基于微波法的生物质热解油改性

第44卷第2期东北林业大学学报Vol.44 No.22016年2月JOURNAL OF NORTHEAST FORESTRY UNIVERSITYFeb. 2016基于微波法的生物质热解油改性’)崔勇常建民李本王文亮虞羽翔任学勇(北京林业大学.北京，10083)摘要使用微波辐射的方法对生物质热解油进行改性研究,考查了微波功率、微波作用时间和pH值等因素对改性效果的影响。结果表明,改性后生物质热解油的固体颗粒质量分数和黏度降低、甲醛消耗量提高,并得出最优工艺为pH值7.5左右左右，微波功率300 w,微波作用时间4 min;通过GC-MS检测发现，生物质热解油改性前后主要化合物种类基本相同,微波改性可以使酚类、酮类、醛类的相对峰面积增加。关键词生物质;热解油;微波;改性Biomass Pyrolysis Oil Modified by Microwave Method//Cui Yong, Chang Jianmin, Li Ben, Wang Wenliang, YuYuxiang, Ren Xueyong( Beijing Forestry University, Bejjing 100083，P. R. China)//Jourmal of Northeast Forestry Univer-sity ,2016,44(2):70-74.We used microwave to modily biomass pyrolysis oil, and studied the effects of microwave power, microwave process-ing time, and pH value. Biomass pyrolysis oil modifed by microwave has lower solid content and viscosity, and higherformaldehyde consumption. The optimal technological parameters were obtained as adjusting pH value to about 7.5, using300 W of microwave processing for 4 min. By GC - MS analysis, main compound species has the same line after microwaveprocessing in the biomass pyrolysis oil, but the contents of phenols, ketones and aldehydes were increased.Keywords Biomass; Pyrolysis oil; Microwave ; Modification生物质热解油(简称生物油)是指生物质原料分发生反应,是导致生物油品质不佳的重要原因。(如木粉、农作物秸秆树皮、果壳等)在隔绝空气的生物油的自身缺陷,使其直接利用变得很困难。近条件下被迅速加热到一定温度热裂解成低分子有机年来,人们尝试使用物理和化学方法改善生物油的物蒸气,再将其冷凝液化后得到的棕黑色液体,具有性能。价格低廉、原料可再生、反应活性较高等特点,被认微波是波长范围为1~1 000 mm的电磁波,在为是一- 种潜在的新能源!。在木材加工行业中,生电磁波谱中位于红外线和无线电波之间,具有波、粒物油苯酚、甲醛三者共缩聚合成的生物油酚醛树二相性,同时具有一-定的能量。目前,微波技术在通脂,已经被证明可以替代传统酚醛树脂应用到人造信、军事、医疗以及科研等诸多领域都得到了广泛应板的生产中。用。微波是通过热效应和非热效应两种方式对化学生物油的化学组成十分复杂,包括苯酚和愈创反应产生影响[。1。研究表明，微波对木质纤维的物木酚等酚类物质、乙酸等酸类物质、乙醛等醛类物理化学特性具有改性作用，被认为是一种颇具潜力质,以及脂类衍生物和木质素低聚物等(2。研究表的活化方法[7-8]。明,生物质基本上都是由纤维素、木质素、半纤维索本研究针对生物油在制备生物油酚醛树脂过程构成。生物油是以上几类组分热解和解聚过程的产中与甲醛化学反应活性不足的问题,采用微波辐射物,因此不同原料和热解工艺产出的生物油在化学的方法对其进行处理,使生物油中未分解的低聚物组成上基本相似”。由于在热解过程中并未达到进一步解聚为小分子类物质,从而达到提高其化学热力学平衡,所以生物油的物理化学性质并不稳定，反应活性的目的。.主要表现为其黏度随贮存期的延长而逐渐升高,同1材料与方法时还有pH值低、热值低、固体杂质含量高等缺点[4。生物油中含有20%~ 30%的木质素低聚1.1材料物[5],其稳定性差,易与生物油中活泼性不饱和组氢氧化钠(分析纯,广东汕头市西陇化工厂有限公司)、甲醇(分析纯,国药集团化学试剂有限公1)国家高技术研究发展计划( 2012AA101808-06)司).乙醇(分析纯,95%,国药集团化学试剂有限公第- -作者简介:崔勇.男, 1984年7月生,北京林业大学材料科司)、盐酸(分析纯ρ=1.19g/mL,国药集团化学试.与技术学院，博士研究生, E-mail :232550184@ qq.com。通信作者:常建民,北京林业大学材料科学与技术学院教授。E剂有限公司) ;浓盐酸(分析纯,37% ,北京化学试剂收稿日期:2015年8月17日。公司);盐酸羟胺(分析纯,北京化学试剂公司);蒸责任编辑:戴芳天。馏水(分析纯,中国科学院半导体所);卡尔费休试第2期崔勇,等;基于微波法的生物质热解油改性71剂( KFR-02型，天津市科密欧化学试剂有限公司);通过滤纸孔隙为3 μm的过滤装置预滤一遍后,再生物油(采用山东泰然新能源有限公司生产的松木通过滤纸孔隙为0.1 μm的滤纸;然后用乙醇溶液洗热解油,颜色为棕褐色)。涤该过滤装置直到滤液变得纯净,并把带有滤渣的1.2 设备滤纸放在烘箱中干燥30 min后,在干燥器中冷却至表1主要仪器和设备室温后称质量并计算滤渣质量。设生物油的质量设备名称规格型号产地m,滤渣的质量为m2 ,则固体颗粒质量分数为m2与美的微波炉EM7KCCW3-NR广东m;的比值。PB-10赛多利斯科学仪器有限公司卡尔费休水分测定仪ZDJ-IS北京先锋威驱技术开发公司1.3.7黏度测定数显恒温水浴锅NH-4金坛荣华仪器制造有限公司将盛有试样的容器(直径不小于6 cm)置于25电子分析天平BS124SC的恒温水浴中,温度达到平衡后,再将选择好的转循环水式真空泵SHB-3 .北京神泰伟业仪器设备有限公司旋转蒸发仪RE-52A上海亚荣生化仪器厂子垂直浸人到试样中心,并使液面达到容器标线,打旋转黏度计NDJ-SS上海昌吉地质仪器有限公司电热鼓式干燥箱101-0AB北京利康达圣科技发展有限公司开旋转黏度计,待测定的黏度值稳定后读取数值。气相色谱-质谱联用仪Agilen16980/5973N 美国Agilen公司1.3.8 甲醛消耗量测定1.3方法测定生物油与甲醛反应过程中消耗甲醛的量能1.3.1 微波改性够直观地体现出生物油的反应活性"。取10 g热解在室温下,利用磁力搅拌器将生物油样品搅拌油和30 g过量的甲醛溶液,在圆底烧瓶中混合,放入均匀后静置一段时间待用。称取一定质量的生物搅拌转子,用质量分数40%的NaOH溶液调节pH值油,加入一定质量的质量分数为40%的NaOH溶液至9 ,将圆底烧瓶置于80 C水浴中反应90 min。调节其pH后,放入微波炉中，在- -定的微波功率下取1g反应后的待测试样放人250 mL的烧杯处理一段时间。中,加人50mL甲醇,用0.9mol/L的HCl溶液调节1.3.2微波功率的确定pH值至3.5;随后加入25mL质量分数10%的盐酸微波的作用功率会影响木质素低聚物的解聚效羟铵溶液，搅拌10 min;最后用0.54 mol/L的NaOH果,本实验将微波的作用功率依次设置为100、150、标准液滴定该试样,使其pH值恢复至3.5以上,并200、250.300、350W,分别对6组pH值为7.32的生用相同步骤做空白实验(不加试样)的测定。物油处理4 min。10g生物油可以消耗甲醛的质量:按照公式计1.3.3 微波作用时间的确定算生物油与甲醛实际反应过程中消耗的甲醛量,反微波的作用时间会影响木质素低聚物的解聚效应中消耗的甲醛质量=反应加人的甲醛质量-溶液果,本实验将微波的作用时间依次设置为1.2、3、4、中游离 甲醛的质量。5.6min,分别对6组pH值为7.32的生物油使用_3e(V,-Vo)100%;300 W功率处理。n2M2=M. -wXm。1.3.4 pH 值的确定利用pH计测得本实验所用的松木生物油pH式中:10为游离甲醛质量分数(%);m,为投料质量值为3.45,呈强酸性。由于木质素低聚物的解聚会;(g) ;m2为取样质量(g) ;M,为反应加入的甲醛质受到pH值的影响,因此本实验通过向生物油中加量(g) ;M2为反应消耗的甲醛质量(g);c为NaOH人一定质量的质量分数为40%的Na0H溶液将pH标准溶液的浓度(mol/L);V,为待测试样所消耗分别调节为3.45.5.55 、6.43 .7.32 8.42 9.51、10.31,NaOH溶液的体积(mL);V为空白试样所消耗在微波功率为300w,微波作用时间为4min的条件NaOH溶液的体积(mL)。下处理。1.3.9主要化学组分 分析1.3.5含水率测定采用美国PerkinElmer公司生产的气相色谱一使用ZDJ-1S型卡尔费休含水率测定仪,以卡质谱联用仪(GC-MS)对生物油的组分进行分析。尔试剂为标准液,甲醇为滴定溶剂。通过标准液消试验前,生物油经离心机以转速5 000 r/min预处理耗体积计算样品含水量[91。20min,然后进行真空抽滤处理,除去混在生物油中1.3.6固体颗粒质量分 数测定的杂质。最后用分液漏斗按V(生物油) : V(二氯采用乙醇溶解法来测量生物油中固体颗粒物的甲烷)=1 : 10将两者混合，待溶液分层,取下层清质量分数[10],将约2g生物油样品溶解在200mL无液作为待测试样。水乙醇溶剂中,m(样品) : m(溶剂)约为1 : 100,先气相色谱条件:色谱柱DB-5 MS柱,尺寸30.0072东北林业大学学报第44卷mx0.25 mmx0.25 μm。程序升温:初始温度60 C ,当微波功率在250 W以下时,生物油的黏度呈保温2 min;以6 C/min的速率升温至220 C ,保温现下降趋势;当超过250W时，生物油的黏度开始增3min;以20C/min的速率升温至280C,保温5加。黏度降低可以增强生物油的流动性,有利于分min。载气为氦气;进样口温度250C ;分流比100: 1;散和混合。微波对大分子有分解作用,可以使生物油取0.5 μL生物油样品注人GC-MS进样器。中未分解的木质素等低聚物发生解聚,变成流动性较质谱条件:电离源EI,电子轰击能量70eV;接好的小分子物质,从而降低生物油的黏度;随着功率口温度260C,离子源温度220C;质谱扫描范围的增加,微波的加热作用逐渐增强,体系内温度升高,m/z=29~500;溶剂延迟时间3min。会促进活性基团之间发生聚合反应,使黏度增加。按照上述试验条件进行测试，得到对应峰的总当微波功率从100W增加到300W时，甲醛的离子流TIC谱图。根据谱库检索以及文献资料,采消耗量呈增加的趋势,此后小幅下降,这表明在一定用峰面积归一法,分析确定松木生物油微波改性前范围内提高微波作用功率有利于提高生物油的反应后的主要化学成分及其相对峰面积变化。活性。综合考虑微波功率对生物油以上各项基本性质的2结果与分析影响,本实验认为将微波功率设置在300 W较合理。2.1微波功率对生物油性质的影响2.2微波作用时间对生物油性质的影响微波功率对生物油性质的影响如表2所示。经微波作用时间对生物油性质的影响如表3所.过微波处理的生物油含水率普遍略高于生物油原示。在规定功率下,微波处理6 min的时间范围内,油,一方面在调节pH时会引人水分导致含水率升生物油含水率呈现先升高再降低的趋势,但均低于高,另一方面微波处理过程中,生物油中的某些组分30% ,波动范围较小。随着微波作用时间的延长,微发生脱水反应也会产生少量水分。生物油在水中的波对生物油的加热作用持续增强,体系内温度逐渐溶解能力有限,如果含水率过高,会破坏体系中的微升高,部分活性基团之间会发生缩合脱水,进而使生乳液结构,造成水相和有机相的分离,这将对生物油物油中水分含量增加;如果体系内温度过高,又会加的后续利用很不利。研究认为，生物油的含水率控快水分及小分子物质的挥发,使含水率略有降低。制在30%以内时,其具有较好的均匀性和稳定表3微波作用时间对生 物油性质的影响性[12]。可以发现，微波功率的增加并没有引起含水微波作用含水固体颗粒质黏度/甲醛消耗量/率的显著变化,在本实验选取的微波功率范围内,含时间/min率/%量分数/%mPa.smol. g~'水率均低于30%,并且生物油具有较好的均一性，026.510.3550.006326.820.280.0075未出现分层现象。27.130.224.0081表2微波功率对生 物油性质的影响27.320.180.008828.580.150.0090微波功29.010.140.009率/W率/% .mPa.s mol●g'28.410.130.009 2生物油固体颗粒质量分数随着微波作用时间的1027.010.3180.007增加呈下降趋势,且初期下降幅度明显，并逐渐减15026.810.230.007 520027.800.210.0078缓。微波加热作用随时间延长而增强，固体颗粒物25028.020.008330028.520.0091被加热到一定程度时会膨胀破裂,变成更小颗粒或3500.120.0079。者分解。任何生物油中都含有一定的固体物质,这些物在6min的微波作用时间内,生物油的黏度呈质主要是炭颗粒以及热解过程中的流化床料等。这现逐渐降低的趋势,表明微波作用时间延长有利于些固体颗粒物可能会成为加速生物油老化反应的催生物油黏度的降低,微波作用带来的温度升高以及化剂,对生物油的储存及使用不利B。随着微波功对低聚物的降解作用是引起上述现象的原因。但是率的增加,固体颗粒质量分数呈现下降趋势,且前期随着微波作用时间的持续增加，降黏作用不断减弱，降低幅度大于后期。微波对生物油中的颗粒具有加这是由于温度升高促进了体系内的聚合反应。热效果,可以使其膨胀发生破裂,进而降低固体颗粒生物油甲醛消耗量随着微波作用时间的延长呈的直径。但是,随着微波强度的逐渐增加,这种效应现出先缓慢增加后基本保持不变的趋势,这表明适并非无限增强,过高的功率强度反而会使生物油内当地增加微波作用时间有利于提高生物油的活性,部发生其他副反应。但是过长的微波作用时间对生物油反应活性的提高第2期.崔勇,等:基于微波法的生物质热解油改性意义不大。moL左右的甲醛,相当于反应能力提高了54%,表综合考虑微波作用时间对生物油以上各项基本明微波改性后生物油的反应活性得到了较大改善。性质的影响,本实验认为将作用时间定为4min较表5重复试验结果合理。含水固体颗粒质黏度/ 甲 醛消耗量/编号2.3 pH 值对生物油性质的影响率/%量分数/% mPa.smol.g!27.20.16310.0110pH值对生物油性质的影响如表4所示。本实228.20.15330.0095验使用的松木生物油原油含水率为26%,在利用27.80.170.009 8NaOH溶液调节pH值的过程中会引入外来水分,随26.20.50500.0065着生物油pH值增高其含水率不断增加。在微波处2.5 化学组分对比分析理的情况下,pH值的增加有利于降低生物油中固体根据总离子流TIC谱图,通过检索谱库及查阅颗粒质量分数;但是当pH增加到- -定值时,固体颗文献资料,采用峰面积归一法,得到落叶松生物油微波改性前后的主要化学组成，如表6、表7所示。可粒质量分数降低的趋势明显减缓。以看出,通过GC-MS对微波改性前后的落叶松生表4 pH 值对生物油性质的影响pH值固体颗粒质度/甲醛消耗量物油化学成分进行分析,主要检测出酚类、酮类、醛率/9量分数/%mPa.sol.g'类、酸类物质,由于离子总图存在部分无法确定的谱3.4526.210.25420.0075峰,表中所列成分相对峰面积的总和小于100%。5.5526.950.21100.008 I改性前落叶松生物油中所含成分相对峰面积由大到6.4327.510.18s50.00857328.530.0090小依次为:酮类(21.05%,以羟基丙酮为主)、醛类8.4229.710.0089(17.73%,以羟基乙醛为主)、酸类( 16.71%,以乙酸9.5130.820.13300.008810.32_31.61.13290.008 7为主)、酚类(16.52%, 以4-甲基愈创木酚为主)。生物油的黏度随着pH值的增加,呈持续降低微波改性后的落叶松生物油主要成分的相对峰面积的趋势,使整个体系的流动性增加。微波的降解作则发生了一定变化，相对峰面积由大到小依次为:酮用和调节pH值时水分的引人，都是引起黏度降低类(24.97% ,以羟基丙酮为主)、酚类(22.79%,以4-甲的原因。基愈创木酚和4-乙基愈创木酚为主)、醛类(20.82%,生物油甲醛消耗量随着pH值的增加,呈现先以羟基乙醛为主)、酸类( 16.66% ,以乙酸为主)。增加后降低的趋势。研究表明,提高NaOH的用量,表6未改性生物油主要化学组成可以有效地促进生物油中与苯酚结构类似的酚类化号类别比例化合物名称相对峰面积/%合物与甲醛的反应活性,并有利于高分子聚合物的酚类(16.52%) 4-甲基愈创木酚2.75形成14)。但是当pH值继续增加时,甲醛消耗量有4-丙烯基-2-甲氧基苯酚4-乙基愈创木酚.03所降低,这可能是由于含水率过高，影响了生物油的2-甲氧基苯酚.88反应活性。2.6-二甲基苯酚.584-甲基邻奉二酚.48综合考虑pH值对生物油以上各项基本性质的4-烯丙基-2,6-二甲氧基苯酚.18影响,本实验认为将pH值调节至7.5左右(弱碱邻苯二酚0.97性)进行处理较合理。对甲苯酚0.782.4工艺条件优化及结果1丁香酚0.53综合以上试验结果,得出微波改性生物油的最123.4-二甲基苯酚0.52佳工艺条件为:生物油的pH值调节至7.5左右(弱2-甲酚碱性),微波功率为300w,微波作用时间为4min。14酮类(21.05%) 羟基丙酮15.984-羟基-3-甲氧基苯丙酮3.96在上述最佳工艺条件范围内进行3次重复性试验，163-甲基-1.2-环戍二酮0.69试验结果见表5。可知,与未经改性的生物油相比,173,4-二羟基-3-环丁烯-1,2-二酮0.42在最佳工艺条件下获得的微波改性生物油的含水率18醛类(17.73%)羟基乙醛14.6119丙醛3.12稳定在28%左右,满足含水率低于30%的使用要:20酸类(16.71%)_ 乙13.13求;改性生物油的固体颗粒质量分数和黏度降低明乙酰氧基乙酸58显,这有利于提高生物油的均- -性以及流动性;甲醛通过微波改性后,生物油中的酚类、酮类、醛类消耗量的实验显示,1 g未改性生物油的甲醛消耗量组分的相对峰面积均有所增加,而酸类组分相对峰为0.0065moL,而1g改性生物油可以消耗0.001面积基本不变。其中，酚类、醛类和酮类组分相对峰74东北林业大学学报第44卷面积的增加,说明在微波作用下生物油中木质素、纤GC-MS分析表明,微波改性前后落叶松生物油维素和半纤维素的未完全热裂解产物发生了断键和的主要成分都是由酚类、酮类、醛类、酸类等几种物质解聚反应。组成;但是,通过微波改性可以使酚类酮类、醛类组表7微波改性生物油主要化学组成分相对峰面积增加,而酸类组分相对峰面积变化不序号类别比例化合物名称相对峰面积/%大,表明在微波作用下生物油中木质素纤维素和半酚类(279%) 4-甲 基愈创木酚3.51纤维素的未完全热裂解产物发生了断键和解聚反应。4-乙基愈创木酚3.324-甲基邻苯二酚2.42参考文献4-丙烯基-2-甲氧基苯酚1.92[1]王树荣 ,骆仲泱,谭洪,等.生物质热裂解生物油特性的分析研2,6-二甲基苯酚1.85究[J].工程热物理学报,2004 ,25( 6) : 1049- 1052.2-甲氧基苯酚1.72[2] MARCELOE D, ANDREC V, YVESS, et al. Coprocessing of邻苯二酚1.45oxygenated biomass compounds and hydrocarbons for the produc-tion of sustainable fuel[ J]. ChemSuChem, 2008,1(3): 179-181.4-嫣丙基-2.6-二甲氧基苯酚1.18 .[3] 武景丽,汪从伟,阴秀丽，等生物油分离方法的研究进展[J].对甲苯酚1.05石油化工,2008,37(1):95-99.4-乙基邻苯二酚0.92 :[4] PEACOCKE G v C, RUSSELL P A, JENKINSJ D, et al. Physi-丁香酚0.91cal proproperties of flash pyrolysis liquids[ J]. Biomass and Bioener-苯酚gy, 194,7:169-177CHAAALA A., PAKDEL H, etal. Character-[5] CARCIA-PEREZ M，3,4-二甲基苯酚0.74 :ization of bio-oil in chemical families[ J]. Biomass and Bioenergy,4-丙基间苯二酚0.2007 ,31(4)222-242.2-甲酚0.27[6] BANIK s, BANDYOPADHYAY s, GANGULY s. Bioffets of5 酮类(24.97%) 羟基丙酮15.82microwave-a brief review [J]. Bioresource Technology .2003 , 874-羟基-3-甲氧基苯丙酮3.43(2):155-159.[7] TSENGKH, YONGFS, YUT Y. A study of biomass maerial182-丁酮2.28pretreatment by microwave-based heating[ J]. Advanced Materials3-甲基-1,2-环戍二酮2.13Research ,2011 ,347-353 :2991-2995.203,4-二羟基-3-环丁烯-1,2-二酮1.31[8] SAPCI Z. The efct of microwave pretreatment on bio-gas produc-21醛类( 20.82%)羟基乙醛14.56tion from agricultural straws[ J ] . Bioresource Technology，2013,128(1) :487-494.丙醛4.31[9]国家标准化管理委员会.原油水含量测定卡尔 .费休库伦滴23香草醛1.95定法:CB/T 1146- 200[ S].北京:中国标准出版社,2009.24 酸类(16.6%6) 乙酸13.04[10] OASMA A A, MEIER D. Norms and standards for fast pyrolysis25亚油酸3.62liquids: I. Round robin test[ J]. Joumal of Analytical and Ap-plied Pyrolysis，2005 ,73(2) :323-334.3结论[11] 国家标准化管理委 员会.木材胶黏剂及其树脂检验方法:GB/T 14074- -2006[ S].北京:中国标准出版社,2006.[12] RADLEIN D. Study of levoglucosan production-a review. Fast py-综合单因素实验结果,得出了微波改性生物油rolysis of biomass: a hand book[ M]. Newbury: CPL Press,的优选工艺为:生物油的pH值调节至7.5左右(弱[13] CHIARAMONTI D, OASMAA A, SOLANTAUSTA Y. Power碱性),微波功率为300W,微波作用时间为4min。Beneraton.using asipytolysisugucs2007, 11(6): 1056- 1086.经此工艺处理的生物油,其固体颗粒质量分数和黏度[14]范东斌尿素 、生物油一苯酚一甲醛共缩聚树脂的合成、结构明显降低,甲醛消耗量提高了54%,反应活性得到了与性能研究[ Dj.北京:北京林业大学,2009.较大提高,同时含水率低于30% ,外观均一-稳定。(上接69页)2009:238-271.[9] 全国塑料标准化技术委员会方法和产品分会.塑料拉伸性能的[15] Bledzki A K, Faruk 0. Creep and impact properties of wood f-测定:第3部分薄膜和薄片 的试验条件:GB/T 1040.3- -2006ber-polypropylene composites: inluence of temperature and mois-[S].北京:中国标准出版社,2007.ture content[ J]. Composites Science & Technology , 2004, 64(5) :693- 700.[10] 全国人造 板标准化技术委员会,木塑地板:GB/T 24508- 2009赵永生，朱复华，薛平,等.木粉对PVC木塑复合材料力学性[S].北京:中国标准出版社,2010.[11] ANATOLE A Klyosov. Wood-plastic composites[M]. New Jersey:能影响1[].现代塑料加工应用2005,17(6):12-15[17] 杜虎虎,李涛 ,王伟宏,等.木粉含量对木粉- -HDPE 复合材料Wiley-Interscience ,2007 :403- 405.[12] 任宁，刘一 星,巩翠芝.木材微观构造与拉伸断裂的关系[J].物理力学性能的影响[J].西南林业大学学报，2013,33(1):东北林业大学学报,2008 ,36<2):33-35.13] 李清玲,张晓文,裴鑫杰.等高密度聚乙烯的拉伸性能影响因[18]刘庄，黄旭江，张翔,等HDPE微发泡塑木复合材料的制备及其性能研究[J].化工新型材料,2014.42(3):71-73.素分析[J].炼油与化T,2012,23(1):16-18.焦书科.橡胶化学与物理导论[M].北京:化学工业出版社,