留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

沉淀法回收生物油高温馏分中的酚类物质

赵欣 李凯 李明 马善为 朱锡锋

赵欣, 李凯, 李明, 马善为, 朱锡锋. 沉淀法回收生物油高温馏分中的酚类物质[J]. 机械工程学报, 2016, 44(2): 201-208.
引用本文: 赵欣, 李凯, 李明, 马善为, 朱锡锋. 沉淀法回收生物油高温馏分中的酚类物质[J]. 机械工程学报, 2016, 44(2): 201-208.
ZHAO Xin, LI Kai, LI Ming, MA Shan-wei, ZHU Xi-feng. Separating phenols from high-temperature fraction of bio-oil by precipitation[J]. JOURNAL OF MECHANICAL ENGINEERING, 2016, 44(2): 201-208.
Citation: ZHAO Xin, LI Kai, LI Ming, MA Shan-wei, ZHU Xi-feng. Separating phenols from high-temperature fraction of bio-oil by precipitation[J]. JOURNAL OF MECHANICAL ENGINEERING, 2016, 44(2): 201-208.

沉淀法回收生物油高温馏分中的酚类物质

基金项目: 

国家重点基础研究发展规划 973计划, 2013CB228103

中国科学院重点部署项目 KGZD-EW-304-3

详细信息
  • 中图分类号: TK6

Separating phenols from high-temperature fraction of bio-oil by precipitation

More Information
  • 摘要: 为了研究金属离子沉淀法对生物油高温馏分中酚类物质的回收, 提高生物油中化学物质的利用率, 利用气相色谱-质谱联用仪 (GC-MS) 从NaOH试剂浓度、反应温度和反应时间三个方面研究了镁离子对生物油中酚类物质的回收效果。结果表明, 镁离子与酚类物质形成了不溶物, 而且不同浓度的氢氧化钠溶液 (1.0-4.0 mol/L)、不同的反应温度 (25-85℃) 以及不同的反应时间 (5-35 min) 对实验结果有着不同程度的影响。结果表明, 在反应温度为25℃、氢氧化钠浓度为2.5 mol/L, 反应时间在20 min时为最佳反应条件。在此条件下, 对生物油高温馏分中对乙基苯酚的回收率可达34.97%。

     

  • 生物油制备技术又称生物质热解液化技术,是指在缺氧状态和中温 (500℃左右) 条件下,生物质快速受热分解,热解气经快速冷凝后主要获得液体产物 (生物油) 以及一部分固体产物 (炭粉) 和气体产物 (燃气) 的热化学转化过程[1]。然而,由于生物油具有黏度高、水分含量大、酸性强以及热稳定性差等特性,致使其品质较低,应用受到一定限制[2]

    生物油的化学组成复杂,除了含有20%-30%的水和少量炭粉等固体颗粒之外,还含有数以百计的酸、醇、酮、醛、酚、醚、酯、糖、呋喃等有机化合物。目前,已检测出的就达400多种[3]。生物油分离制取化学品是生物油高值化的一个重要研究方向。对生物油进行分离提取精细化学品,不仅能提高生物油的利用价值,还能提高生物质热解制油技术的竞争力,因而引起了世界各国学者的密切关注。目前,从生物油中分离获得的产品,主要有酚类[4]、酸类[5]、糖类[6]和醛类[7, 8]等。

    其中,酚类物质主要来源于生物质中木质素的热解。木质素具有天然的酚类结构,在热解过程中可以产生多种酚类化合物,如愈创木酚、二甲氧基苯酚及其酚类的烷基化的衍生物等[9, 10]。目前,从生物油中分离出的酚类主要有苯酚、愈创木酚、4-甲基愈创木酚、2,6-二甲氧基苯酚等,大多是采用碱化-酸化与有机溶剂萃取相结合的手段实现的。比如Amen-Chen等[4]用碱化-酸化和有机溶剂萃取出桉树木焦油真空蒸馏馏分中的酚类。Li等[11]利用了蒸馏、碱化-酸化的工艺,即先获得富含酚类的馏分,再通过碱化-酸化、二氯甲烷萃取、减压蒸馏得到纯净的酚类。目前,从生物油中分离酚类仍存在有机溶剂用量过大、易造成环境污染、回收困难等问题。因此,亟待寻找一种高效、低成本且环保的分离方法,以大规模回收生物油中酚类物质。

    Ge等[12]利用金属离子与酚类在碱性条件下会生成沉淀这一特点,研究探索了钡离子沉淀法对低温热解焦油170-220℃高温馏分中酚类的回收过程,获得了理想的结果。Wang等[13]则通过向生物油中加入氢氧化钙,与生物油中的乙酸反应生成钙离子后,再和酚类生成络合物沉淀,从而达到分离酚类物质的目的。这种方法工艺简单、成本低且对环境的污染小,是一种有前景的分离生物油中酚类物质的方法。基于此,实验利用酚类会与大多数金属在碱性条件下产生不可溶络合物,而不溶物又会在酸性条件下释放出酚类物质这一特点[14],研究了镁离子对生物油高温馏分中酚类物质的回收,主要考察了镁离子应用在此方法的可行性。此外,实验还探讨了不同的操作参数 (NaOH的浓度、反应时间、反应温度) 对实验的影响,试图确定一个最佳工艺参数。

    实验所用生物油是以稻壳为原料,由每小时处理量为120kg物料的自热式生物质热解液化中试装置制备得到。

    对乙基苯酚 (97%) 购买于阿拉丁试剂,硫酸镁 (98%)、片状氢氧化钠 (96.0%)、硫酸 (98%)、乙酸乙酯 (99.5%) 均购买于国药集团化学试剂有限公司。

    实验中酚类的提取过程见图 1。在室温下,将300g的生物油进行常压蒸馏,得到60g的高温馏分 (170-250℃)。量取1mL的高温馏分放入小烧杯中,边搅拌边缓慢加入9mL的饱和硫酸镁溶液,持续搅拌30min。然后边搅拌边缓慢加入8mL不同浓度的NaOH溶液,作为pH值调节剂使得溶液pH值大于7,保证金属离子与酚类沉淀反应的顺利进行。淡黄色的不溶物在滴加NaOH的过程中开始析出,然后混合物用中速滤纸在真空泵下进行过滤,之后将不溶物放入真空干燥箱中在105℃下放置4h,然后称重。将干燥过后的不溶物溶于10mL的2mol/L H2SO4中,30min后加入20mL的乙酸乙酯,用于萃取酚类物质。酚类物质可以通过对乙酸乙酯层进行旋蒸得到,而乙酸乙酯则可循环利用。

    图  1  沉淀法回收生物油中酚类物质的工艺路线示意图
    Figure  1.  Extraction procedure of phenolic compounds from bio-oil

    其中,NaOH的浓度分别设为1.0、1.5、2.0、2.5、3.0、3.5、4.0mol/L,不溶物析出时的反应时间分别设为5、10、15、20、25、30、35min,不溶物析出时的反应温度分别设为25、35、45、55、65、75、85℃。

    溶液pH值由上海虹益仪器仪表有限公司pHB-8型笔式pH计进行测量。色谱柱为TR-35MS毛细管柱 (30m×0.25mm×0.25μm),分流比为1:80。GC程序升温条件为:起始温度40℃保持5min,以4.0℃/min的升温速率升温到180℃保持2min,再以10.0℃的升温速率升温至280℃保持2min,终温为280℃,离子源温度为280℃。根据NIST谱库并参考已知的生物油组分,确定样品的化学组成,每次实验重复三次并取平均值。采用美国热电仪器公司Nicolet 8700型傅里叶变换红外光谱仪对不溶物官能团进行分析。

    采用外标法对对乙基苯酚进行定量分析,用乙酸乙酯作溶剂,分别配置不同浓度的对乙基苯酚溶液 (0.20%、0.26%、0.33%、0.39%、0.46%)。研究中对乙基苯酚的萃取率 (%) 为其在乙酸乙酯相中含量 (g) 与生物油中含量 (g) 的百分比。

    由于生物油的成分复杂,直接提取生物油中高附加值化学品的难度较大。而酚类物质的沸点大多处于170-250℃,考虑到实验的目标产物为酚类,故可以先对生物油进行常压蒸馏[15],这样得到的高温馏分主要成分即为酚类物质,从而有利于对其进行进一步分离和提取,实验将高温馏分进行GC/MS分析,然后将得到色谱图中主要峰所对应的化合物及其相对峰面积见表 1

    表  1  生物油高温馏分中酚类物质的化学成分
    Table  1.  Chemical compositions of phenolic compounds from high temperature fractions of bio-oil
    No.Chemical compoundMolecular formulaMolecular weight /(g·mol-1)StructurePeak area /%
    1phenolC6H6O944.47
    2phenol, 2-methyl-C7H8O1085.67
    3phenol, 3-methyl-C7H8O1086.56
    4phenol, 2-methoxy-C7H8O21243.35
    52, 5-xylenolC8H10O1221.88
    6phenol, 4-ethyl-C8H10O12221.5
    7phenol, 2-methoxy-4-methyl-C8H10O21382.86
    81, 2-benzenediolC6H6O21103.41
    9phenol, 3-(1-methylethyl)-C9H10O1343.05
    10phenol, 3-methyl-4-ethyl-C9H10O1341.58
    11phenol, 4-isopropyl-C9H12O1363.89
    12guaiacol, 4-ethyl-C9H10O21503.75
    131, 2-benzenediol, 4-methyl-C7H8O21242.1
    14phenol, 4-sec-butyl-C10H14O1502.37
    15phenol, 2, 6-dimethoxy-C8H10O31543.15
    16guaiacol, 5-allyl-C10H12O21642.51
    17phenol, 2-methoxy-4-propyl-C10H14O21662.94
    18catechol, 4-ethyl-C8H10O21382.2
    19guaiacol, 4-propenyl-C10H12O21642.04
    下载: 导出CSV 
    | 显示表格

    表 1可知,高温馏分中主要包含了苯酚、间甲基苯酚、愈创木酚、对乙基苯酚等多种酚类物质。其中,对乙基苯酚、苯酚、愈创木酚含量较多,具有较高的回收价值,而实验主要研究高温馏分中对乙基苯酚的分离回收效果。

    此不溶物的样品是在反应温度为25℃、氢氧化钠浓度为2.5mol/L,反应时间20min的最佳反应条件下制备得到。对此不溶物进行FT-IR分析,得到不溶物的红外光谱谱图见图 2

    图  2  不溶物的红外光谱谱图
    Figure  2.  FT-IR spectrum of the complex

    图 2可知,在3400-3500cm-1存在羟基峰。如果酚类物质与Mg2+反应生成的不溶物是以盐的形式存在而不是络合物,则在不溶物的红外谱图中是不会出现羟基峰的。

    然后把不溶物溶于10mL的H2SO4(2.0mol/L),再用20mL的乙酸乙酯进行萃取,最后用GC-MS进行分析。图 3为乙酸乙酯层的离子色谱图。

    图  3  乙酸乙酯层的总离子色谱图
    Figure  3.  Total ion chromatogram of ethyl acetate layer

    表 2图 3中主要峰所对应的酚类化合物及其相对峰面积。由图 3表 2可知,用乙酸乙酯萃取得到的主要是酚类物质,其中, 对乙基苯酚的相对峰面积最大 (21.75%),值得注意的是萃取得到的酚类中还是有一些其他蒸馏成分,这可能是因为这些物质附着在不溶物上,而且很难用乙醇清洗干净,从而混入了酚类物质中。同时也证明了不溶物可以通过加入硫酸的方法得到对乙基苯酚,这进一步说明了实验的可行性。

    表  2  乙酸乙酯层中的酚类化合物
    Table  2.  Phenolic compounds identified in the ethyl acetate layer
    No.Chemical compoundMolecular formulaMolecular weight /(g·mol-1)StructurePeak area /%
    1phenolC6H6O945.48
    2phenol, 3-methyl-C7H8O1086.47
    3guaiacolC7H8O21244.03
    4phenol, 4-ethyl-C8H10O12221.75
    5guaiacol, 4-methylC8H10O21383.74
    6catecholC6H6O21103.19
    7phenol, 3-(1-methylethyl)-C9H10O1343.95
    8phenol, 2-ethyl-5-methyl-C9H12O1362.09
    9phenol, 4-propyl-C9H12O1364.01
    10guaiacol, 4-ethyl-C9H10O21505.6
    11catechol, 4-methyl-C8H10O21381.77
    12phenol, 4-sec-butyl-C10H14O1501.48
    13phenol, 2, 6-dimethoxy-C8H10O31542.78
    14guaiacol, 5-allyl-C10H12O21561.85
    15phenol, 2-methoxy-4-propyl-C10H14O21663.96
    16catechol, 4-ethyl-C8H10O21381.56
    17guaiacol, 4-propenyl-C9H12O1362.1
    下载: 导出CSV 
    | 显示表格

    pH值对不溶物的形成具有重要的影响。为了研究pH值对对乙基苯酚萃取率的影响,将NaOH浓度分别设为1.0、1.5、2.0、2.5、3.0、3.5和4.0mol/L,同时将反应温度和反应时间分别设定为25℃和20min。图 4为NaOH浓度对于对乙基苯酚萃取率的影响。由图 4可知,随着NaOH浓度从1.0mol/L增加到2.5mol/L,对乙基苯酚的萃取率则从12.49%迅速增加到37.28%。而当NaOH浓度大于2.5mol/L时,对乙基苯酚的萃取率几乎不再增加。根据溶液中不溶物达到溶解平衡时的溶度积方程式Ksp{Mg (OH)2}=c(Mg2+c(OH-)2,上述现象可能是因为,当NaOH浓度小于2.5mol/L时,c(Mg2+)c(OH-)2并没有达到氢氧化镁的溶度积,此时沉淀主要是酚类与镁离子形成的络合物;而当NaOH浓度大于2.5mol/L时,c(Mg2+)c(OH-)2超过了氢氧化镁的溶度积,此时氢氧化镁开始沉淀,并成为主要的沉淀物,萃取率几乎不再发生变化。同时发现,此时溶液的黏度也不断增加,这可能是产生了过多的氢氧化镁沉淀导致的,这不仅增加了下一步分离的难度,也使得一部分高温馏分附着在了沉淀上,从而导致反应难以进行。综上所述,2.5mol/L的NaOH浓度被选为进一步实验的最佳条件。

    图  4  在不同NaOH浓度下对乙基苯酚的萃取率和pH值
    Figure  4.  4-ethylphenol extraction rate and pH value at various concentrations of sodium hydroxide solution

    为了研究反应温度对于对乙基苯酚萃取率的影响,把反应温度分别设为25、35、45、55、65、75和85℃。图 5为在不同温度下对乙基苯酚的萃取率和pH值。由图 5可知,温度对于对乙基苯酚的萃取率具有较大影响,较低的反应温度有利于酚类与金属离子沉淀的形成。当反应温度从25℃升高到85℃时,萃取率从34.97%降低为30.46%。这可能是因为不溶物的溶解度随温度的升高而升高造成的[13],而且随着温度的升高,酚类还易受热氧化[16],所以在高温时沉淀加速溶解,萃取率降低。综上所述,从实验成本和操作难度上综合考虑,将反应温度25℃作为反应的最佳温度。

    图  5  在不同温度下对乙基苯酚的萃取率和pH值
    Figure  5.  4-ethylphenol extraction rate and pH value at various temperatures

    为了研究反应时间对萃取率的影响,先将NaOH浓度和反应温度分别设为2.5mol/L和25℃。图 6为在不同时间下对乙基苯酚的萃取率和pH值。由图 6可知,当反应时间小于20min时,萃取率随反应时间的延长而显著增加,而当反应时间延长到20min时,萃取率变化极小。此类金属离子与酚类的络合反应都需要一定的反应时间才能达到平衡[14, 16, 17],而且酚类是一种强还原性物质,易被空气氧化,在碱性介质中不稳定,应尽量避免酚类在碱性介质中长时间存留。由图 6还可知,大约20min时酚类物质与金属离子的络合反应已经反应完全,再延长时间已经没有意义。综上所述,将反应时间20min作为最佳反应时间。

    图  6  在不同时间下对乙基苯酚的萃取率和pH值
    Figure  6.  4-ethylphenol extraction rate and pH value at various reaction times

    以生物油高温馏分为研究对象,研究了Mg2+对于其中酚类物质的回收,探讨了NaOH浓度、反应温度和反应时间对对乙基苯酚萃取率的影响。结果表明,利用Mg2+分离生物油高温馏分中的酚类物质是一种简单有效的方法;最佳的反应温度、反应时间和NaOH浓度分别为25℃、20min和2.5mol/L。在上述条件下,对乙基苯酚的萃取率可达34.97%。

  • 图  沉淀法回收生物油中酚类物质的工艺路线示意图

    Figure  1.  Extraction procedure of phenolic compounds from bio-oil

    图  不溶物的红外光谱谱图

    Figure  2.  FT-IR spectrum of the complex

    图  乙酸乙酯层的总离子色谱图

    Figure  3.  Total ion chromatogram of ethyl acetate layer

    图  在不同NaOH浓度下对乙基苯酚的萃取率和pH值

    Figure  4.  4-ethylphenol extraction rate and pH value at various concentrations of sodium hydroxide solution

    图  在不同温度下对乙基苯酚的萃取率和pH值

    Figure  5.  4-ethylphenol extraction rate and pH value at various temperatures

    图  在不同时间下对乙基苯酚的萃取率和pH值

    Figure  6.  4-ethylphenol extraction rate and pH value at various reaction times

    表  1  生物油高温馏分中酚类物质的化学成分

    Table  1.   Chemical compositions of phenolic compounds from high temperature fractions of bio-oil

    No.Chemical compoundMolecular formulaMolecular weight /(g·mol-1)StructurePeak area /%
    1phenolC6H6O944.47
    2phenol, 2-methyl-C7H8O1085.67
    3phenol, 3-methyl-C7H8O1086.56
    4phenol, 2-methoxy-C7H8O21243.35
    52, 5-xylenolC8H10O1221.88
    6phenol, 4-ethyl-C8H10O12221.5
    7phenol, 2-methoxy-4-methyl-C8H10O21382.86
    81, 2-benzenediolC6H6O21103.41
    9phenol, 3-(1-methylethyl)-C9H10O1343.05
    10phenol, 3-methyl-4-ethyl-C9H10O1341.58
    11phenol, 4-isopropyl-C9H12O1363.89
    12guaiacol, 4-ethyl-C9H10O21503.75
    131, 2-benzenediol, 4-methyl-C7H8O21242.1
    14phenol, 4-sec-butyl-C10H14O1502.37
    15phenol, 2, 6-dimethoxy-C8H10O31543.15
    16guaiacol, 5-allyl-C10H12O21642.51
    17phenol, 2-methoxy-4-propyl-C10H14O21662.94
    18catechol, 4-ethyl-C8H10O21382.2
    19guaiacol, 4-propenyl-C10H12O21642.04
    下载: 导出CSV

    表  2  乙酸乙酯层中的酚类化合物

    Table  2.   Phenolic compounds identified in the ethyl acetate layer

    No.Chemical compoundMolecular formulaMolecular weight /(g·mol-1)StructurePeak area /%
    1phenolC6H6O945.48
    2phenol, 3-methyl-C7H8O1086.47
    3guaiacolC7H8O21244.03
    4phenol, 4-ethyl-C8H10O12221.75
    5guaiacol, 4-methylC8H10O21383.74
    6catecholC6H6O21103.19
    7phenol, 3-(1-methylethyl)-C9H10O1343.95
    8phenol, 2-ethyl-5-methyl-C9H12O1362.09
    9phenol, 4-propyl-C9H12O1364.01
    10guaiacol, 4-ethyl-C9H10O21505.6
    11catechol, 4-methyl-C8H10O21381.77
    12phenol, 4-sec-butyl-C10H14O1501.48
    13phenol, 2, 6-dimethoxy-C8H10O31542.78
    14guaiacol, 5-allyl-C10H12O21561.85
    15phenol, 2-methoxy-4-propyl-C10H14O21663.96
    16catechol, 4-ethyl-C8H10O21381.56
    17guaiacol, 4-propenyl-C9H12O1362.1
    下载: 导出CSV
  • [1] 朱锡锋.生物油制备技术[M].北京:化学工业出版社, 2013.

    ZHU Xi-feng.Preparation and application of bio-oil[M].Beijing:Chemical Industry Presss, 2013.
    [2] OASMAA A, CZERNIK S.Fuel oil quality of biomass pyrolysis oils-state of the art for the end user[J].Energy Fuels, 1999, 13:914-921. doi: 10.1021/ef980272b
    [3] BRIDGWATER A V, PEACOCKE G V C.Fast pyrolysis processes for biomass[J].Renew Sust Energy Rev, 2000, 4:1-73. doi: 10.1016/S1364-0321(99)00007-6
    [4] AMEN-CHEN C, PAKDEL H, ROY C.Separation of phenols from Eucalyptus wood tar[J].Biomass Bioenergy, 1997, 13(1/2):25-37.
    [5] MAHFUD F H, VAN GEEL F P, VENDERBOSCH R H, HEERES H J.Acetic acid recovery from fast pyrolysis oil:An exploratory study on liquid-liquid reactive extraction using aliphatic tertiary amines[J].Sep Sci Technol, 2008, 43(11/12):3056-3074.
    [6] GANDER M, RAPP K M, SCHIWECK H.Process for preparing 1, 6-β-anhydroglucopyranose (levoglucosan) in high purity:US, 5023330[P].1991-06-11.
    [7] STRADAL J A, UNDERWOOD G L.Process for producing hydroxyacetaldehyde:US, 5252188[P].1993-10-12.
    [8] STRADAL J A, UNDERWOOD G L.Process for producing hydroxyacetaldehyde:US, 5393542[P].1995-02-28.
    [9] DINESH M, CHARLES U P, PHILIP H S.Pyrolysis of wood/biomass for bio-oil:A critical review[J].Energy Fuels, 2006, 20:848-889. doi: 10.1021/ef0502397
    [10] EFFENDI A, GERHAUSER H, BRIDGWATER A V.Production of renewable phenolic resins by thermochemical conversion of biomass:A review[J].Renew Sust Energy Rev, 2008, 12:2092-2116. doi: 10.1016/j.rser.2007.04.008
    [11] LI J, WANG C, YANG Z.Production and separation of phenols from biomass-derived bio-petroleum[J].J Anal Appl Pyrolysis, 2010, 89(2):218-224. doi: 10.1016/j.jaap.2010.08.004
    [12] GE Y Z, JIN H.Recovery process for phenolic compounds from coal-derived oils by ions of solublemetal salts[J].Fuel, 1996, 75:1681-1683. doi: 10.1016/0016-2361(95)00276-6
    [13] WANG D, LI D B, LIU Y Q, LV D C, YE Y Y, ZHU S J, ZHANG B B.Study of a new complex method for extraction of phenolic compounds from bio-oils[J].Sep Purif Technol, 2014, 134:132-138. doi: 10.1016/j.seppur.2014.07.033
    [14] 葛宜掌, 金红, 李斌.Ba2+离子沉淀法回收酚类化学过程的研究[J].燃料化学学报, 1996, 24(3):266-270. http://www.cnki.com.cn/Article/CJFDTOTAL-RLHX603.014.htm

    GE Yi-zhang, JIN Hong, LI Bin.Study on recovery of phenolic compounds using Ba2+ as precipitant[J].J Fuel Chem Technol, 1996, 24(3):266-270. http://www.cnki.com.cn/Article/CJFDTOTAL-RLHX603.014.htm
    [15] JEWEL A.CAPUNITAN, SERGIO C.CAPAREDA.Characterization and separation of corn stover bio-oil by fractional distillation[J].Fuel, 2013, 112:60-73. doi: 10.1016/j.fuel.2013.04.079
    [16] 明珠, 吴云海, 阿依妮尕尔·艾尔肯.茶汤与金属离子络合反应的研究[J].环境科技, 2015, 28(4):5-11. http://youxian.cnki.com.cn/yxdetail.aspx?filename=JSHJ2015072100F&dbname=CAPJ2015

    MING Zhu, WU Yun-hai, ARKEN·AYINIGAR.Study on complexation reaction of tea infusion with metal ions[J].Environ Sci Technol, 2015, 28(4):5-11. http://youxian.cnki.com.cn/yxdetail.aspx?filename=JSHJ2015072100F&dbname=CAPJ2015
    [17] GE Y, JIN H.Recovery of phenols from coal tar and waste water by precipitation[J].J China Coal Soc, 1995, 20:545-548.
  • 加载中
图(6) / 表(2)
计量
  • 文章访问数:  88
  • HTML全文浏览量:  116
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-09-18
  • 修回日期:  2015-11-19

目录

/

返回文章
返回