更全的杂志信息网

超精煤热解行为及动力学研究

更新时间:2009-03-28

0 引 言

目前,煤炭资源仍是世界能源消费的主体[1],对不同品级煤炭能源的高效利用始终是研究的热点。近年来通过可循环有机溶剂对劣质煤进行热熔萃取处理[2]成为煤炭资源综合利用的研究热点。Takanohashi等[3]使用1-甲基萘对不同煤阶煤在不同温度条件下进行处理,发现高温条件有利于提高产率,而且低阶煤处理后具有更高的热塑性。Shui等[4]在室温条件下使用CS2/NMP混合有机溶液处理非炼焦煤,有效提高了其黏结性指数,进行混煤炼焦时其焦炭强度显著提高。通过可循环有机溶剂对煤进行高温萃取,萃取混合液经分离有机溶剂后所得的萃取物称为超精煤(hyper coal,HPC)[5]。HPC具有低灰、低硫的特点,作为清洁能源使用可以有效减少污染物排放,减少对反应容器的损害。樊丽华等[5]对无灰煤的热解行为进行研究,发现无灰煤中中小型分子增多,热解初始温度降低,热失重量增加。目前对于HPC的研究主要集中在热熔萃取的机理研究和较好炼焦能力的应用上,而对于热解过程尤其是热解动力学研究相对较少。煤的热解是煤炭气化与燃烧等煤资源转换方式的初始阶段[6],并且对产物反应活性和转化有重要影响[7]。热解过程的内在机理通过宏观反应动力学的方式表现出来[8]。在热重分析仪器上采用程序升温法研究煤热解过程中的失重现象是一种成熟的手段,其具有试验量小,操作周期短,获得信息多等特点[9]。在对热解过程失重数据进行处理时,将其看为一个整体并不能真正反映煤复杂的热解过程,因此对热解过程分阶段研究可以更好地描述热解过程[10]。平传娟等[7]利用热天平试验测试数据对混煤热解动力学进行研究,发现对混煤热解分段拟合可以很好地描述不同温度段热解机制。

本文采用升温热重分析法在氩气下,对原煤和超精煤热解特性和动力学进行研究。根据热解特征参数对热解过程进行分段处理。采用分段模型拟合法研究确定热解不同阶段的动力学模型及其对应的热解动力学参数。

1 试 验

1.1 样品制备

以岢岚(KL)煤为试验用煤。使用可循环有机溶剂(N-甲基吡咯烷酮,NMP)在高压反应釜内对KL煤粉进行热萃取处理。试验所选高压反应釜有效容积为1 L,按煤粉与溶剂比为8 g∶400 mL,准确称取相应量的KL煤粉和NMP溶剂加入高压反应釜。首先向反应釜内以400 mL/min的速率通入高纯氮气20 min,排除空气后关闭进气口和出气口,此时釜内为大气压力。设置程序控制反应釜升温至350℃后恒温1 h。整个试验过程反应釜内桨叶以100 r/min转速搅拌,使煤粉与有机溶剂两相充分接触。在350℃下保温1 h后关闭控温程序,取下加热套空冷至室温,取出反应物。使用布氏漏斗和滤纸分离萃取混合液,获得液相和残渣。将滤液通过真空旋转蒸发仪回收有机溶剂和析出固体产物。固体产物经酒精和去离子水反复冲洗后放入真空干燥箱(80℃)真空干燥12 h,干燥后产物即为超精煤。将原煤和超精煤分别标记为KL-RAW和KL-HPC。相应煤样的工业及元素分析见表1。样品均研磨至粒度<0.074 mm,以备热重分析使用。

3.1.2 膀胱穿孔 膀胱穿孔是TURBT需要特别警惕的并发症,术前应常规进行CT/MRI扫描,了解肿瘤浸润深度和基底大小。电切过深、膀胱过度充盈及闭孔反射是穿孔的主要原因,电切时膀胱灌入液体不能太多,以免膀胱过度膨胀,使膀胱壁变得太薄而容易穿孔;电切侧壁肿瘤时需警惕闭孔反射发生,切除侧壁肿瘤前,可适当加深麻醉予以肌松剂甚至闭孔神经阻滞,适当充盈膀胱,甚至请手术助手适当固定同侧下肢,电切时尽量采取间歇式触发电切模式,均可能降低由于闭孔反射导致的膀胱穿孔发生。此外,术中应仔细辨认结构,一旦切除组织底部见到脂肪组织时,提示已经穿孔,应立即停止这一区域的电切。

 

表1 煤样的工业及元素分析Table 1 Proximate and ultimate analysis of coal samples %

  

注:∗表示根据差减法所得数据。

 

煤样 工业分析元素分析Mad Vad Aad FCad Cad Had O∗ad Nad Sad KL-RAW 2.14 32.20 10.37 55.29 73.38 4.73 7.48 1.59 0.96 KL-HPC 2.11 37.52 1.20 59.17 79.13 5.27 6.85 4.14 0.53

1.2 试验仪器及过程

试验设备为美国TA仪器的SDT-600系列热重分析仪。试验过程:温度范围为室温至1 500℃,升温速率0.1~50 K/min,样品容量200 mg,天平灵敏度 0.1 μg,DTA 灵 敏 度 0.001 K。 每 次 取 样 品10.0 mg(粒径<0.074 mm)。 反应温度为室温至1 000℃,采用高纯氩(99.999%)为载气,流量为100 mL/min,升温速率分别为 5、10、20 K/min。

1.3 动力学模型

煤热解动力学的研究目的在于定量表征反应过程及确定反应过程中遵循的最佳机理函数,求出动力学参数和,算出速率常数k,提出模拟TA(thermal analysis)曲线的反应速率表达式。通常煤热解过程表示为

面对小学教育多媒体教学的新形势,我们应该首先加强教师思想教育,加强党员和教师的先锋模范作用,因为人们的思想决定了动机、行为和创造力,决定了教育的方向和质量。其次是管理上要跟得上时代,对于小学英语教师的培养要遵守人才发展规律。再次是加强教师的教研改革活动,可以给英语课堂注入新的活力和魅力,甚至可以形成一定的校园氛围,从而全面提高教学质量。还有要适度对教师进行中长期进修培训,让教师有针对性、实效性、全身心的投入到培训中去。

 

3)分段模型拟合法可以很好地解释煤热解过程机理。原煤的低温和中温段的热解机制均为一维扩散模型,高温段热解机制为二级化学反应模型;超精煤的3个温度段分别为一维扩散模型、三维扩散模型和二级化学反应模型。升温速率的提高对原煤热解的中温段影响较大,而对超精煤热解的低温过程影响较大;同时加热速率的提高对超精煤3个温度段的活化能总和并无影响。

 

对于非等温、非均相反应的动力学方程为

 

[5]樊丽华,杜敬文,梁英华,等.无灰煤的热解行为及其在配煤中的添加效果[J].煤炭科学技术,2017,45(3):185-190.FAN Lihua,DU Jingwen,LIANG Yinghua,et al.Pyrolysis behavior of hypercoal and its adding effect in blending coal[J].Coal Science and Technology,2017,45(3):185-190.

根据热重分析数据,热解转化率α

当衰减值分别为20dB和30dB时,图6和图7均有准确测量数据,图7和图8没有测量数据,说明文献[1]和文献[2]的测量方法可准确测量fc=|2ft±fj|类型互调发射,但受到频谱分析仪动态范围的限制,无法对fc=|3ft±fj|类型的互调发射进行测量.文献[3]的测量结果与本文的测量结果一致.

 

式中,m为样品质量,mg;m0为样品初始质量,mg;mf为样品反应终止质量,mg。

因为公司的拍摄人员不懂专业,不知道需要的是哪一种效果,所以请公司拍摄时一定要有同专业的教师来协助,一个人操作、一个人看着拍摄画面,确保拍摄角度和画面达到想要的效果。拍摄顺序不一定按照实验操作的先后顺序,要考虑实验的时间长短、拍摄过程是否需要等待、等待时间长短等,用最经济的时间完成拍摄。

f(α)为热解反应动力学机理函数,表示物质反应速率与反应转换率α之间所遵循的某种函数关系。相应的积分形式定义为

 

其中,P(u)为温度积分,其没有解析解,只有数值解。动力学方程的求解最终归结在对P(u)的处理上,本文采用Coasts-Redfern法[11]P(u)进行近似,即

 

将式(6)代入式(5)得

 

对式(7)两边取自然对数,得到Coats-Redfern方程为

 

对于一般的反应温度区与大部分的E值来说所以式(8)变为

对于占全国学钢琴人数99.9%的人来说,都是属于非专业钢琴学习,那么,非专业钢琴学习应该用什么标准来学习?是否“业余”就代表可以随便放任呢?我认为,“好”的标准是统一的。如果用学习英语来比喻,可能更容易明白。许多学生准备出国,会被要求考托福、雅思……这些考核不是代表英语专业,但是这些非专业性质的考核都有非常规范和严格的系统与标准。钢琴也是一样,钢琴学习应该是很规范而且系统的,应该是循序渐进,踏踏实实,一步一个脚印的。那么钢琴课应该怎么上,怎样才是更负责的,更有效果的呢?

2.5.2.2 农业防治首先进行合理轮作,不能和瓜类、向日葵、茄科等作物轮作,可以和禾本科、大豆、甜菜等轮作。其次是进行深耕秋翻,将种子翻入20cm以下土壤深处,可以减轻食葵列当种子的萌发。最后就是及时铲除列当,在列当开花结籽前,结合除草将列当铲除。

 

对于正确的反应机制 f(α)或 G(α),,因此对1/T是一条直线,可判定反应模型是否正确,常见反应机理函数的微分和积分形式见表2[12-13],包括扩散模型,随机成核模型、收缩核模型和反应级数模型。

 

表2 常用气固反应动力学函数[12-13]Table 2 Basic gas-solid reaction kinetics function[12-13]

  

函数 符号 反应机理 微分形式 积分形式D1 一维扩散 (1/2)α-1α2扩散(Dm) D2 二维扩散 [-ln(1)]-1α+(1)ln(1)D3 三维扩散 (3/2)[(1)-1/3-1]-1 1-2α/3-(1)2/3随机成核(Am) A2m=2 2(1)[-ln(1)]1/2 [-ln(1)]1/2 A3m=3 3(1)[-ln(1)]2/3 [-ln(1)]1/3 S2 收缩圆柱体(面积) 2(1)1/2 1-(1)1/2收缩核(Sm) S3 收缩球状(体积) 3(1)2/3 1-(1)1/3 F1n=1 1-ln(1)化学反应(Fm) F2n=2 (1)2 (1)-1-1 F3n=3 (1)3 [(1-α)-2-1]/2

2 试验结果与讨论

2.1 不同样品热解特性

KL-RAW和KL-HPC煤在20℃/min升温下热解的TG和DTG曲线如图1所示。可以看出,2种煤的热解失重过程具有类似的特征,失重速率先增大后减小。所以煤的热解过程主要分为3个阶段:结晶水和吸附气体脱除阶段,热解气体和半焦生成阶段以及半焦缩聚焦化阶段。从DTG曲线可知,KL-RAW煤热解过程仅在500℃附近存在1个失重速率较大的失重峰,为热解过程热解气体大量生成阶段;KL-HPC煤热解过程除在与KL-RAW煤相近温度存在1个失重速率较大的失重峰外,在较低温度存在2个失重峰。超精煤制备过程中,高温下有机溶剂对煤结构产生溶胀作用,使煤结构疏松,制取的KL-HPC空隙中吸附有少量NMP有机溶剂;另一方面,KL煤粉经有机溶剂高温萃取后,大分子聚合结构打开、断裂,产生较多的芳香自由基碎片,主要以中、小分子量的状态存在。导致超精煤中易挥发组分含量增大,因此在KL-HPC煤粉热解过程DTG图的较低温度存在2个失重峰。

  

图1 KL-RAW和KL-HPC在20℃/min升温下TG和DTG曲线Fig.1 TG and DTG curves of pyrolysis of KL-RAW and KL-HPC at 20℃/min

1)原煤和超精煤热解过程主要分为3个阶段:结晶水和脱气阶段,热解气体和半焦生成阶段以及半焦缩聚焦化阶段。尽管原煤的最大热解失重速率比超精煤大,超精煤具有更高的热解温度区间,造成更大的最终失重率。

 

表3 不同样品热解特征参数Table 3 Characteristic parameters of different samples pyrolysis reaction

  

煤种 Ts/℃Tmax/℃Tf/ ℃(dw/dT)max/(%·℃-1)Tf-Ts KL-RAW 328 473 612 0.249 0 284 KL-HPC 183 467 752 0.186 8 568

2.2 加热速率对热解影响

不同升温速率下样品热解失重过程如图2所示,其特征参数见表4。可以看出,升温速率的提高对KL-RAW和KL-HPC热解过程的影响类似,随着升温速率的提高,热解失重曲线向高温区移动,对应的TsTmTf值均增大,表现出热滞后性[9]。但是相应的增幅KL-RAW较大,因此加热速率对其影响较大。升温速率的提高对样品最大热解速率的影响不大,主要是因为其表达为单位温度的质量变化分数(%/℃),若将其单位转化为单位时间的质量变化分数(%/min),可以看出随着升温速率的提高单位时间内失重量显著增多,这样的单位表达更具有实际意义。

  

图2 不同升温速率TG和DTG曲线Fig.2 TG and DTG curves of pyrolysis at different heating rates

 

表4 不同升温速率条件下热解特征参数Table 4 Characteristic parameters of pyrolysis reaction at different heating rates

  

煤样 β/(℃·min-1) Ts/℃ Tmax/℃ Tf/℃(dw/dT)max/(%·℃-1) /(%·min-1)304 442 580 0.251 1 1.255 5 KL-RAW 10 306 456 605 0.250 2 2.502 0 20 328 473 612 0.249 0 4.980 0增幅 24 31 32 5 5 166 445 725 0.184 0 0.920 0 KL-HPC 10 173 458 742 0.184 4 1.844 0 20 183 467 752 0.186 8 3.376 0增幅 17 22 27

由于煤热解过程的复杂性,不同温度阶段内其热解过程存在差异。一般地,根据热解过程和质量变化速率的不同,将整个热解温度范围划分为不同阶段。本文主要采用TsTmax为节点,将整个热解温度区间划分为低温、中温和高温段,见表5,此区间的划分同样也是动力学参数求解的依据。

 

表5 不同升温速率温度范围Table 5 Temperature ranges at different heating rates

  

煤样 β/(℃·min-1)温度/℃低温段 中温段 高温段5 100~304 304~442 442~1 000 KL-RAW 10 100~306 306~456 456~1 000 20 100~328 328~473 473~1 000 5 100~166 166~445 445~1 000 KL-HPC 10 100~173 173~458 458~1 000 20 100~183 183~467 467~1 000

2.3 动力学参数的求解

采用模型拟合法对煤热解动力学研究时,将整个热解过程视为一个整体有时不能很好地分析热解过程的机理,因此应对热解反应过程依据表5进行分段处理,对不同速率条件下热解过程每个阶段分别进行模型拟合。根据表2中所列动力学机理函数,结合式(9)进行线性拟合,热解过程的机制为拟合过程中相关系数R2最大的机理函数所描述的动力学模型。根据最佳机理函数的拟合直线所对应的斜率和截距分别计算得到活化能E和指前因子A,即可获得“热解动力学三因子”。结合表6,分析KL-RAW和KL-HPC煤在升温速率10℃/min下热解不同温度阶段的动力学参数,可以看出,KLRAW煤的第2个阶段的活化能最高,其次是第3阶段,而第1个阶段的活化能最低。这主要是因为煤的热解过程主要发生在第2阶段,而第1阶段主要发生结晶水的脱除和脱气,随着温度的升高,第3阶段主要发生半焦缩聚焦化过程,这与前文所阐述的热解过程相符。对于KL-HPC,第1阶段活化能最高,随着反应的进行活化能逐渐降低,这可能主要是在制备KL-HPC时存在未完全分离出去的有机溶剂导致的活化能升高。

不同升温速率下热解过程的最佳动力学机理函数见表6,其中的动力学参数可以很好地反映出随升温速率的提高,所出现的“热滞后”现象。在不同升温速率条件下,相应阶段的最佳动力学机理函数并不发生变化,对于不同升温速率条件下KL-RAW煤热解过程,低温段和中温段热解机制均为一维扩散模型,高温段热解机制为二级化学反应模型;对于KL-HPC煤热解过程,3个温度阶段热解机制分别为一维扩散模型、三维扩散模型、二级化学反应模型。对于KL-HPC煤,随着升温速率的变化,每个温度阶段的活化能发生变化,但活化能总和基本保持不变,这与文献[15]中结论一致。不同温度区间活化能的变化情况可以看出,升温速率增大对KLRAW煤中温段和KL-HPC煤低温段影响较大。

3 结 论

2种样品的热解特征参数见表3。其中,热解开始温度Ts为转化率随温度变化中,转化率5%与50%的连线与横坐标的交点所对应的温度;热解最大质量变化速率温度Tmax为热解速率随温度变化中,质量变化速率最大的点对应的温度,对应热解最大质量变化速率为(dw/dT)max;热解终温Tf根据文献定义为 Tf=2Tmax-Ts[14]。可以看出,虽然 KLHPC煤的最大质量变化速率低于KL-RAW煤,但是KL-HPC煤的起始热解温度更低,热解终温更高,导致其有更宽的热解温度区间,因此KL-HPC在整个热解过程中表现出更大的质量变化速率。

因此,正确发挥社会公众监督的作用,特别是新闻监督的作用,既符合我国宪法的立宪价值,又对于弱势群体实现健康权等其他基本权利具有重要意义。但是即便如此,我国目前并未针对舆论监督或者公众监督进行完善专门立法,这着实令人感到遗憾。

 

表6 不同升温速率条件下KL-RAW和KL-HPC热解动力学参数Table 6 Kinetics parameters of the pyrolysis of KL-RAW and KL-HPC at different heating rates

  

煤样 β/(℃·min-1) T/℃ 机理 Ea/(kJ·mol-1) A/min-1R2 100 ~304 D1 19.32 2.55×10-3 0.947 1 5 304 ~442 D1 92.18 4.08×104 0.943 3 442 ~1 000 F2 28.22 3.41 0.981 7 100 ~306 D1 15.46 9.57×10-4 0.950 6 KL-RAW 10 306 ~456 D1 87.11 1.26×104 0.934 8 456 ~1 000 F2 35.71 13.0 0.988 3 100 ~328 D1 16.58 9.52×10-4 0.937 2 20 328 ~473 D1 101.71 1.32×105 0.947 9 473 ~1 000 F2 43.43 47.0 0.988 1 100 ~166 D1 96.62 3.03×107 0.974 7 166 ~445 D3 45.18 3.54 0.991 2 445 ~1 000 F2 23.60 1.92 0.982 8 100 ~173 D1 69.09 5.33×103 0.903 6 KL-HPC 10 173 ~458 D3 46.68 4.28 0.989 0 458 ~1 000 F2 31.55 8.53 0.992 7 100 ~183 D1 50.93 16.3 0.919 7 20 183 ~467 D3 48.50 5.38 0.986 9 467 ~1 000 F2 39.98 37.1 0.989 7 5

基于质量作用定律和Arrhenius方程,均相反应的动力学方程为

参考文献(References):

[1]樊东黎.世界能源现状和未来[J].金属热处理,2011,36(10):119-131.FAN Dongli.World energy:Present and future[J].Heat Treatment of Metals,2011,36(10):119-131.

[2]赵雪飞,吴凯,赖仕全,等.热解温度条件下溶剂萃取煤技术及应用[J].辽宁科技大学学报,2012,35(2):113-117.ZHAO Xuefei,WU Kai,LAI Shiquan,et al.Application of solvent extraction coal technology under pyrolytic temperature condition[J].Journal of University of Science and Technology Liaoning,2012,35(2):113-117.

[6]SOLOMON P R,FLETCHER T H,PUGMIRE R J.Progress in coal pyrolysis[J].Fuel,1993,72(5):587-597.

2)随着升温速率的提高,热解失重曲线向高温区移动,对应的TsTmaxTf值均增大,表现出热滞后性,且对原煤影响较大。以TsTmax为节点可以将热解过程划分为3个阶段:低温、中温和高温段。划分同样适用热解动力学求解过程。

[4]SHUI H,YE W,WANG Z,et al.Hydrothermal treatment of a subbituminous coal and its use in coking blends[J].Energy&Fuels,2013,27(1):138-144.

式中,α为转化率;A为指前因子,min-1;E为活化能,kJ/mol;R为气体常数;T为反应温度,K。

[3]TAKANOHASHI T,SHISHIDO T,KAWASHIMA H,et al.Characterisation of hyper coals from coals of various ranks[J].Fuel,2008,87(4-5):592-598.

原始图像为RGB图像,后续识别与计算需要在灰度图像基础上完成,并且图像中必然存在噪声。在裂缝识别中,灰度值差异为重要基础,即在裂缝部位,其图像灰度值可能产生跳变,出现图像边缘。对于裂缝图像,它是对裂缝实际状态进行可视化表达,需要对其予以灰度转化和除噪,同时进行像素率求解,只有这样,才可以对裂缝实际状况进行定量描述。

菌落总数检测,参照GB4789.2-2010《食品微生物学检验-菌落总数确定》。大肠杆菌总数检测,参照GB4789.3-2010《食品微生物学检验-大肠菌群计数》。

[7]平传娟,周俊虎,程军,等.混煤热解反应动力学特性研究[J].中国电机工程学报,2007,27(17):6-10.PING Chuanjuan,ZHOU Junhu,CHENG Jun,et al.Research on the pyrolysis kinetics of blended coals[J].Proceeding of the CSEE,2007,27(17):6-10.

[8]赵毅,汪黎东,王小明,等.烟气脱硫产物-亚硫酸钙非催化氧化的宏观反应动力学研究[J].中国电机工程学报,2005,25(8):116-120.ZHAO Yi,WANG Lidong,WANG Xiaoming,et al.Study on the macroscopical uncatalyzed oxidation kinetics of desulfurization residual product-calcium sulftte[J].Proceeding of the CSEE,2005,25(8):116-120.

[9]DONG K S,SANG S P,YONG T K,et al.Study of coal pyrolysis by thermo-gravimetric analysis(TGA)and concentration measurements of the evolved species[J].Journal of Analytical&Applied Pyrolysis,2011,92(1):209-216.

[10]DU R,WU K,ZHANG L,et al.Thermal behavior and kinetic study on the pyrolysis of Shenfu coal by sectioning method[J].Journal of Thermal Analysis&Calorimetry,2016,125(2):959-966.

[11]COATS A W,REDFERN J P.Kinetic parameters from thermogravimetric data[J].Nature,1964,201:68-69.

依托成都市医疗废物处理中心项目停车场边坡工程项目,综合计算结果、理论和经验对边坡稳定性进行评价,最后在合理评价的基础上,设计了合理的支护方式,大大提高了场地的利用率。依托边坡工程治理设计综合采用了防滑桩、挡土板以及截水沟等进行联合支护,并采用ANSYS-FLAC 3D联合仿真对支护结构防滑坡效果进行了模拟,结果表明本文采用的联合支护方式效果较好;后期检测数据尚未发现地表隆起、桩的变形,周边地基土变形等现象,目前使用正常,表明了设计方案合理有效。

[12]NIU Z,LIU G,HAO Y,et al.Investigation of mechanism and kinetics of non-isothermal low temperature pyrolysis of perhydrous bituminous coal by in-situ FTIR[J].Fuel,2016,172:1-10.

[13]VLAEV L T,MARKOVSKA I G,LYUBCHEV L A.Non-isothermal kinetics of pyrolysis of rice husk[J].Thermochimica Acta,2003,406(1/2):1-7.

[14]朱丰.COREX熔融气化炉内块煤裂解机理及喷煤燃烧行为研究[D].重庆:重庆大学,2013:11.

第三个坚持,就是岗位分析,对学生就业的主要岗位临床、社区、家庭护理进行分析,提炼出所需要的知识、技能和态度,确定课程标准。

[15]何佳佳,邱朋华,吴少华.升温速率对煤热解特性影响的TG/DTG 分析[J].节能技术,2007,25(4):321-325.HE Jiajia,QIU Penghua,WU Shaohua.Study on the effects of heating-up speed to coal pyrolysis with TG/DTG analysis[J].Energy Conservation Technology,2007,25(4):321-325.

 
郭龙飞,左海滨,赵骏,王亚杰,王静松
《洁净煤技术》 2018年第02期
《洁净煤技术》2018年第02期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息 粤ICP备2023046998号