畜禽粪便生物炭的二维红外光谱分析
桂向阳1,2, 刘晨3, 许吉宏1, 段方蕾1, 方书伟1, 李飞跃1,*
1.安徽科技学院资源与环境学院, 安徽 凤阳 233100
2.上海交通大学中英国际低碳学院, 上海 200240
3.上海第二工业大学环境与材料工程学院, 上海 201209
*通讯联系人 e-mail: lifeiyue0523@163.com

作者简介: 桂向阳, 1995年生, 安徽科技学院资源与环境学院硕士研究生 e-mail: 304134321@qq.com

摘要

生物炭作为一种新型的环境功能材料, 由于其特殊的结构和理化性质, 具有修复污染、 改良土壤、 固碳减排等多重环境效益, 越来越受到人们的关注, 已经成为一个新的研究热点。 生物炭表面官能团的种类及分布, 决定其理化性质。 而热解温度和原料类型是影响生物炭官能团的两个重要因素。 以鸡粪、 牛粪和猪粪为原料, 在不同温度(200~700 ℃)下热解制备畜禽粪便生物炭。 在一维红外光谱分析的基础上, 通过傅里叶变换红外光谱技术结合二维相关分析方法, 揭示畜禽粪便生物炭的官能团随着热解温度升高的变化规律, 为更好的研究生物炭的构效关系提供理论基础。 结果表明: 随着热解温度的升高, 在3 600~2 800 cm-1波段, 3种畜禽粪便生物炭的官能团变化主要体现在醇、 酚—OH和脂肪族—CH2峰减弱, 其中, 鸡粪生物炭的醇、 酚—OH和脂肪族—CH2峰变化最明显。 此外, —OH峰的变化强度高于—CH2并且先于—CH2变化。 在1 800~800 cm-1波段, 鸡粪、 牛粪和猪粪生物炭受响应的特征峰数目分别为6, 5和6个, 各特征峰的变化强度顺序分别为鸡粪生物炭: 芳环C=C>C=O>C—O/Si—O>P—O/C—H; 牛粪生物炭: C=O>C=C和P—O/C—H>C—OH和C—O/Si—O; 猪粪生物炭: C=O>C=C>COO->羧酸C—OH>P—O/C—H=C—O/Si—O。 3种畜禽粪便生物炭的部分官能团变化相似: 即C=O和C=C官能团均发生了断裂及重组, 且C=O断裂先于C=C。 然而, 不同粪便生物炭的官能团变化也存在差异, 主要表现在鸡粪生物炭中C—O/Si—O变化先于P—O/C—H, 牛粪生物炭中P—O/C—H变化先于C—O/Si—O, 猪粪生物炭中C—O/Si—O与P—O/C—H变化同时发生。

关键词: 畜禽粪便; 生物炭; 红外光谱; 二维相关分析; 官能团
中图分类号:O657.3 文献标志码:A
Two-Dimensional Perturbation Correlation Infrared Spectroscopy Analysis of Animal Manure Biochar
GUI Xiang-yang1,2, LIU Chen3, XU Ji-hong1, DUAN Fang-lei1, FANG Shu-wei1, LI Fei-yue1,*
1. College of Resource and Environment, Anhui Science and Technology University, Fengyang 233100, China
2. China-UK Low Carbon College, Shanghai Jiaotong University, Shanghai 200240, China
3. School of Environmental and Materials Engineering, Shanghai Second Polytechnic University, Shanghai 201209, China
*Corresponding author
Abstract

More attention has been paid to biochar which has been hot research area as a novel environmental functional material due to its special structure and physicochemical properties, it has multiple environmental benefits such as pollution remediation, soil improvement, carbon sequestration and emission reduction. The physicochemical properties of biochar were determined by the types and distribution of surface functional groups. Pyrolysis temperature and material type are two important factors affecting functional groups of biochar. In this article, animal manure biochars were prepared under different temperature from 200 ℃ to 700 ℃ with the materials of chicken manure, dairy manure and pig manure. On the basis of traditional one-dimensional infrared spectrum analysis, the change tendency of functional groups of animal manure biochars with the increase of pyrolysis temperature was revealed using the method of Fourier transform infrared spectroscopy (FTIR) combined with two-dimensional correlation analysis, providing a theoretical basis for better research on structure-activity relationship of biochar. The results showed that with the increase of pyrolysis temperature, Which variations mainly existed in the unceasing weakening of alcohols and phenols —OH peak and aliphatic —CH2 peak of manure biochar, among which the changes of chicken manure were the most obvious. Furthermore, the changing intensity of —OH was higher than —CH2 and was removed before —CH2 in the range of 3 600~2 800 cm-1. The number of auto-peaks was 6, 5, 6 respectively in the two-dimensional perturbation correlation infrared spectroscopy of chicken, dairy and pig manure biochar. The changing intensity of each auto-peak was (chicken manure biochar) aromatic C=C> C=O >Si—O/C—O>P—O/C—H; (dairy manure biochar) C=O > C=C and P—O/C—H>C—OH, C—O/Si—O; (pig manure biochar) C=O > C=C >COO->Carboxylic C—OH> P—O/C—H=C—O /Si—O. The groups of C=O and C=C broke and reconstituted, and C=O broke ahead of C=C. However, there were differences in the functional groups of different animal manure biochars, mainly the change of C—O/Si—O was prior to P—O/C—H in chicken manure biochar, while P—O/C—H in dairy manure biochar changed earlier than C—O/Si—O, and C—O/Si—O and P—O/C—H in pig manure biochar occurred simultaneously in the range of 1 800~800 cm-1.

Keyword: Animal manure; Biochar; Fourier transform infrared spectroscopy; Two-dimensional correlation analysis; Functional groups
引言

我国是畜禽养殖大国, 每年畜禽粪便排放量高达3.8× 109 t, 然而有效处理率却不足50%[1]。 若粪污随意排放, 其携带的病原菌、 滤液和异味等会给人类生存和环境带来巨大危害[2]。 传统的处理畜禽粪便的方法主要有: 堆肥发酵[3]、 生产沼气[4]等, 然而这些处理方法面临二次污染问题。 因此, 亟需研发畜禽粪便处理新技术。

近年来, 通过低温热解方法将畜禽粪便制备成生物炭, 不仅可以减少粪便直接排放带来的环境危害, 还可以提高废物利用效率[5], 受到广泛关注。 此外, 生物炭施入土壤后, 还能够提高土壤肥力, 改善土壤环境, 促进作物增产[6]。 大量研究表明, 热解温度和原材料是影响生物炭性质最关键的因素[7], 特别是生物炭官能团的变化, 直接决定生物炭在环境中的功能和效应。 因而, 明确生物炭官能团随热解温度的变化规律, 为人工设计“ 功能性” 生物炭提供理论基础, 具有十分重要的研究意义。

傅里叶变换红外光谱技术被广泛应用于生物炭表面官能团的表征与分析。 郑庆福等研究了随热解温度升高玉米秸秆的官能团变化, 发现高温使得— OH, — CH3, — CH2和— C=O之间发生了缔合或脱落, 从而使得— OH, — CH2以及— C=O等的吸收峰均有减弱的趋势, 反之, 芳香族基团吸收峰逐渐增强。 李飞跃等[8, 9]在研究不同热解温度下核桃壳和花生壳的结构变化时也发现随着温度的升高生物炭的烷基, 羧基和羟基等基团均逐渐减少, 芳香化程度逐渐升高。 这些研究集中体现在各官能团的自身变化特性上, 不能将官能团之间的变化联系起来, 导致其应用时具有局限性[10]

Noda等[11, 12]弥补了上述的缺陷, 提出二维红外相关光谱技术(two-dimensional perturbation correlation infrared spectroscopy, 2D-PCIS), 不仅可以提高红外光谱图识别, 还可以将分子内或分子间的相互作用及变化趋势很好的显现出来[12], 受到越来越多研究者的青睐。 刘标等[10]分析了低灰竹材250~600 ℃热解过程中焦炭化学成分及结构演变过程, 发现氢键与烷基C— H变化趋势一致, 并早于自由羟基与芳烃C— H发生变化。 醚类C— O— C、 酮类— — C=O、 羟基、 烷基等不断脱除, 期间伴随着官能团间的相互转化, 芳环稠环化大量发生。 Chen等[13]利用2D-PCIS研究了烟茎炭在250~950 ℃的结构变化, 结果表明, 脱甲基早于脱羟基作用, C=C和C=O的拉伸早于C— H的变化。 然而, 将2D-PCIS应用于不同温度下制备的粪便生物炭官能团变化的研究还很少。 为此, 本文采用2D-PCIS, 探究不同热解温度(200~700 ℃)下畜禽粪便(鸡粪, 牛粪, 猪粪)生物炭的官能团变化规律, 为全面认识生物炭的构-效关系积累数据支持。

1 实验部分
1.1 畜禽粪便原材料

鸡粪、 牛粪和猪粪原料取自安徽省凤阳县某养殖场。 采集来的畜禽粪便原料置于烘箱中105 ℃烘干, 然后粉碎过60目筛, 混匀后装袋备用。

1.2 生物炭的制备

称取一定质量的鸡粪、 牛粪和猪粪原料分别装于50 mL坩埚中, 填满压实后, 用锡纸封口, 并用注射器针头在锡箔纸上扎上均匀的小孔(利于热解过程中气体的排出), 置于马弗炉中, 当热解温度到达设定温度时(200~700 ℃)保留2 h, 关闭马弗炉, 待其冷却后, 获得的黑色固体粉末即为生物炭, 鸡粪、 牛粪和猪粪生物炭分别记为JF200— JF700(字母代表生物质种类, 数字代表热解温度), NF200— NF700和ZF200— ZF700, 生物炭过100目筛, 混匀后装袋备用。

1.3 方法

将生物炭粉末与溴化钾粉末(光谱纯)按质量比1:100混合, 研磨均匀, 压制成片, 再置于傅里叶红外变换光谱仪(Nicolet is10)中进行分析测试。 光谱扫描范围400~4 000 cm-1, 分辨率4 cm-1, 扫描次数64次。 将红外光谱数据导入2D-shige软件进行分析, 最后将分析结果导入Origin 2017软件绘制二维红外相关同异步谱图和自动峰强度图。 各官能团的相对峰强度变化为Iv/Imax, 其中Iv表示该官能团在相应温度下的特征峰强度, Imax为该特征官能团在所有热解温度下的最大吸收峰强度[14]

1.4 二维红外相关光谱

Noda[11, 12]提出了二维相关光谱的数值计算。 首先, 明确一种扰动(温度, 时间, 浓度等)所引起的光谱强度y(v, t)在时间TminTmax之间变化, 其中vt分别代表波数和扰动。 该变化可通过式(1)表示

y̅(v, t)=y(v, t)-y̅(v)TmintTmax0其他(1)

其中, y̅(v)定义为

y̅(v)=1Tmax-TminTminTmaxy(v, t)dt(2)

因此, 二维相关光谱X(v1, v2)可以被表示为

X(v1, v2)=< y˙(v1, t)·y˙(v2, t)> (3)

其中X(v1, v2)表示在固定的外部扰动区间内, 不同的波长v1v2上所测量到的光谱强度变化 y˙(v1, t)。 符号< > 表示用于比较光谱强度 y˙(v1, t)变化的相关性函数。 为了更好地运用于数值计算, Noda将式(3)改为复数形式

X(v1, v2)=Φ(v1, v2)+(v1, v2)(4)

根据Noda规则, 同步谱中位于对角线上的峰称为自相关峰, 一般为正值, 其峰值大小可用来反映受外界扰动的作用下该基团受影响程度, 表明不同波数Φ (v1, v2)之间的基团变化; 非对角线上的峰称为交叉峰, 其值有正有负, 交叉峰为正表明两不同波数光谱强度的相似性变化, 即同增或同减, 交叉峰为负表明光谱强度变化的不一致性。 异步谱Ψ (v1, v2)代表了光谱强度的相异性变化, 用以区分谱峰来源和区别被掩盖的重叠峰; 此外, 如表1所示, 交叉峰还可以用来反映不同波数基团的变化先后顺序, 当同步谱和异步谱中交叉峰信号一致时, 表明v1对应的基团变化先于v2, 不一致时v2先于v1, 当同步谱交叉峰存在但异步谱为0时, v1v2同时变化, 当仅存在异步交叉峰时, v1v2的变化无法确定。

表1 谱峰分析规则 Table 1 Noda’ s rules for peaks analysis
2 结果与讨论
2.1 生物炭的一维红外谱图分析

图1中为3种生物炭的一维红外谱图。 其主要的特征峰如下: 3 600~3 300 cm-1处的吸收峰为醇酚羟基存在的缔合— OH伸缩振动, 随着热解温度升高, 3种生物炭的羟基峰强度均显著下降, 这可能是由于3种生物炭在热解过程中伴随着结合水的脱离, 导致羟基键逐渐断裂[15]。 2 970~2 860 cm-1处为脂肪族C— H的伸缩振动或— CH2烷烃的对称、 反对称伸缩[16], 当热解温度超过400 ℃时, 3种生物炭中该峰强度均明显减弱直至消失, 可见, 随着热解温度升高, 3种粪便生物炭的脂肪族烷基烃逐渐消失, 粪便生物炭的芳香性增强; 1 750~1 600 cm-1处为羧基、 醛、 酮和醌等的C=O伸缩振动[8], 其峰强度也随着温度升高而逐渐减弱, 其原因可能是随着热解温度的升高, 会伴随更多的生物油和热解气体的产出, 而C=O的不饱和状态较易生成气体或液体副产物, 因此C=O键显著减少; 1 100~1 030 cm-1处为多糖物质C— O键或硅酸盐杂质的Si— O键伸缩振动[17, 18]; 1 500~1 400 cm-1处的C=C伸缩振动以及1 000~700 cm-1波段的C— H伸缩振动[16], 表明3种生物炭中均含有一定的芳香结构, 一维体系中其变化不明显。

图1 生物炭的红外光谱图
(a): 鸡粪; (b): 牛粪; (c): 猪粪
Fig.1 Infrared spectra of biochars
(a): Chicken manure; (b): Dairy manure; (c): Pig manure

2.2 生物炭的2D-PCIS分析

为了更好的分析生物炭官能团受温度的影响, 利用一维红外数据进行相关分析并扩展到二维, 得到二维红外相关同步、 异步谱图及自动峰强度图。 由图1所知, 3种生物炭的红外光谱图变化主要发生在3 600~2 800及1 800~800 cm-1波段, 因此, 将FTIR光谱分别在3 600~2 800及1 800~800 cm-1波段做二维相关分析。

2.2.1 3种生物炭在3 600~2 800 cm-1波段的2D-PCIS分析

图2为JF, NF和ZF在3 600~2 800 cm-1波段的二维红外相关谱图。 JF出现2个自相关峰, 分别位于3 375和2 930 cm-1处; NF出现3个自相关峰, 分别位于3 415, 3 045和2 930 cm-1(肩峰)处, 而在一维图谱中[图1(b)]3 045 cm-1处的基团未被识别出来, 可见, 2D-PCIS更加有利于将NF与其他粪便区分出来; ZF出现2个自相关峰, 分别位于3 350和2 927 cm-1处。 不同波数对应的基团归属如表2 所示。 3种生物炭在3 600~2 800 cm-1波段的自动峰变化强度高低, 均为— OH伸缩振动> — CH2的C— H反对称伸缩。 此外, NF出现的烯烃— C=H伸缩振动介于— OH伸缩振动和— CH2的C— H反对称伸缩之间。

图2 生物炭在3 600~2 800 cm-1波段的相关谱图
(a): 鸡粪生物炭; (b): 牛粪生物炭; (c): 猪粪生物炭左: 同步谱图; 中: 异步谱图; 右: 自相关峰强度图
Fig.2 Correlation spectra of biochars in 3 600~2 800 cm-1 band
(a): Chicken manure biochar; (b): Cow dung biochar; (c): Pig manure biocharleft: synchronous spectrum; middle: asynchronous spectrum; right: auto-correlation peak intensity graph

表2 生物炭谱峰的基团归属 Table 2 Group identification of biochar peak

3种生物炭的— OH和— CH2在同步谱图中形成的交叉峰Φ (3 375, 2 930), Φ (3 415, 2 930)和Φ (3 350, 2 927)均为正值, 表明羟基和亚甲基基团呈一致性变化, 即随着热解温度增加, 羟基和亚甲基峰强度减弱, 这和图1分析结果相符。 异步谱图中3种生物炭— OH和— CH2形成的交叉峰Ψ (3 375, 2 930)、 Ψ (3 415, 2 930)和Ψ (3 350, 2 927)均为正值, 与同步谱图中的符号相同, 表明随着热解温度的升高, 畜禽粪便生物炭脱羟基先于脱亚甲基。 此外, NF在Φ (3 045, 2 930)处的交叉峰为正, 而Ψ (3 045, 2 930)处的交叉峰符号为负, 表明亚甲基变化先于烯烃

图4(a, b, c)分别为JF, NF和ZF的羟基和亚甲基的相对峰强度变化趋势图, 3种生物炭的羟基均有显著降低。 同时, JF的羟基降低程度明显高于NF和ZF, 这可能由于NF和ZF中含有的部分伯醇基(C— OH)在生物质水分脱除过程中具有一定贡献[19], 因此, 自由羟基的下降程度相对JF较低。 此外, JF和NF分别在400和500 ℃有上升的趋势, 这可能是因为该温度下生物炭芳环发生断裂导致更多的羟基暴露出来。 JF的亚甲基变化趋势和羟基相同, 随着热解温度升高呈指数式降低。 不同的是, 随着热解温度升高, NF和ZF的亚甲基先增加, 均在300 ℃达到最大, 随着热解温度进一步升高, 吸收峰强度逐渐下降。

图4 相对峰强度变化趋势图Fig.4 Variation in relative intensities of groups

2.2.2 3种生物炭在1 800~800 cm-1波段的2D-PCIS分析

图3为JF, NF和ZF在1 800~801 cm-1波段的二维红外相关谱图。 从自相关峰强度图中可以看出JF升温过程中出现6个主要自相关峰, 分别位于1 605, 1 434, 1 085, 1 045, 957和880 cm-1; NF升温过程中出现5个主要自相关峰, 分别位于1 614, 1 440, 1 235, 1 108和965 cm-1处; ZF出现了6个自相关峰, 分别位于1 655, 1 448, 1 380, 1 265, 1 038和976 cm-1处。 不同波数对应的基团归属如表2所示。 3种生物炭在升温过程中引起的官能团变化相似, 但各官能团的变化强度有所不同, 其中, JF官能团变化强度依次为芳环C=C> C=O > Si— O或C— O> P— O> C— H; NF为C=O > C=C和P— O/C— H> C— OH和C— O/Si— O; ZF为C=O > C=C> COO-> 羧酸C— OH> P— O/C— H, C— O/Si— O。

图3 生物炭在1 800~800 cm-1范围内的相关谱图
(a): 鸡粪生物炭; (b): 牛粪生物炭; (c): 猪粪生物炭; 左: 同步谱图; 中: 异步谱图; 右: 自相关峰强度图
Fig.3 Correlation spectra of biochars in 1 800~800 cm-1 band
(a): Chicken manure biochar; (b): Cow dung biochar; (c): Pig manure biochar left: synchronous spectrum; middle: asynchronous spectrum; right: auto-correlation peak intensity graph

同步谱图中交叉峰符号可以得出官能团之间是否存在相同或相异的变化趋势, 由图3(左)可知, JF的C=C, C=C, C— O/Si— O峰与其P— O/C— H的峰强变化相反。 NF的C=C, C=O, C— O/Si— O和羧酸C— OH峰的变化与其P— O/C— H峰变化相反。 ZF的C=C, C=O, COO-和羧酸C— OH峰与P— O/C— H, C— O/Si— O的峰强变化相反。 根据图3同步谱和异步谱图中的交叉峰符号, 得出生物炭官能团变化的先后顺序。 JF中C— O/Si— O最先发生变化, 其次是C=O断裂, 再是C=C断裂, 最后发生C— H/P— O的变化; NF中C=O断裂先于P— O/C— H先于C— OH先于C— O/Si— O, 最后是C=C断裂; ZF先发生C=O断裂, 其次是C=C , COO-和羧酸C— OH的变化, 最后是C— O/Si— O和P— O/C— H的增加。

图4(d, e, f)分别为JF, NF和ZF的C=O, C=C等特征峰的相对峰强度变化趋势图, 与亚甲基的C— H变化相反, JF的芳环C— H逐渐升高, NF和ZF先减低, 在300 ℃后逐渐增加, 芳环稠化开始发生。 在整个热解过程中, JF和NF的羧基C=O在300 ℃前均有上升, 这可能是在低温热解下部分含氧官能团向羧基C=O转化, 之后被取代、 聚合而逐渐降低[20]。 相应的, JF和NF的C— O官能团在300 ℃之后均有增加, 在500 ℃之后因C— O开始生成醚类结构而降低[9]。 ZF在不断升温过程中, C=O持续降低, 其C— O含量在300 ℃之后也在持续增加。 此外, JF的芳环结构在500 ℃遭到显著破坏, 而NF和ZF在300 ℃就已经被破坏。

3 结论

(1) 3种生物炭在3 600~2 800 cm-1波段的二维红外相关光谱中羟基吸收峰最强, 亚甲基次之。 此外, 羟基和亚甲基变化趋势一致(均降低), 但随着热解温度升高, 脱羟基先于脱亚甲基。

(2) 3种生物炭在1 800~800 cm-1波段的二维红外相关光谱中的相关峰主要由C=O, C=C, C— O, P— O/C— H等的变化引起, 其中C=O和C=C基团变化较为强烈, 且3种生物炭C=O断裂均先于C=C。 JF的C— O/Si— O变化先于P— O/C— H, NF的P— O/C— H变化先于C— O/Si— O, 而ZF的C— O/Si— O与P— O/C— H变化同时发生。

(3) 随着热解温度的升高, 3种生物炭的羟基逐渐降低。 JF的亚甲基含量于200 ℃开始逐渐降低, 而C=O在300 ℃时才发生明显的缔合破裂, C— O/Si— O于300 ℃后出现升高并于500 ℃显著降低, 同时, 其芳香化程度逐渐增强, 芳环结构于500 ℃后破坏显著。 而NF和ZF的亚甲基、 C=C, C— OH先升高, 300 ℃达到最大, 之后逐渐减弱, 它们的P— O/C— H则先降低, 300 ℃之后逐渐增大。 ZF和NF的C=O分别从200和300 ℃发生断裂。

参考文献
[1] WANG Tie-ling(王铁灵). China Animal Industry(中国畜牧业), 2016, (21): 38. [本文引用:1]
[2] Cantrell K B, Hunt P G, Uchimiya M, et al. Bioresource Technology, 2012, 107(2): 419. [本文引用:1]
[3] Jiang Z Y, Yin W, Sui-Cheng H E. Phosphate & Compound Fertilizer, 2008, 23(4): 37. [本文引用:1]
[4] ZHANG Tian, BU Mei-dong, GENG Wei(张田, 卜美东, 耿维). Chinese Journal of Ecology(生态学杂志), 2012, 31(5): 1241. [本文引用:1]
[5] Ro K S, Cantrell K B, Hunt P G. Industrial & Engineering Chemistry Research, 2010, 49(20): 10125. [本文引用:1]
[6] Arif M, Ali K, Jan M T, et al. Field Crops Research, 2016, 195: 28. [本文引用:1]
[7] Tang J, Li X, Luo Y, et al. Chemosphere, 2016, 152: 399. [本文引用:1]
[8] LI Fei-yue, TAO Jin-guo, WANG Jian-fei, et al(李飞跃, 陶进国, 汪建飞, ). Chinese Journal of Environmental Engineering(环境工程学报), 2017, (6): 390. [本文引用:2]
[9] LI Fei-yue, WANG Jian-fei, XIE Yue, et al(李飞跃, 汪建飞, 谢越, ). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2015, 31(4): 266. [本文引用:2]
[10] LIU Biao, DONG Ya-bin, HE Tao, et al(刘标, 董亚斌, 何涛, ). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2016, 32(10): 232. [本文引用:2]
[11] Noda I. Vibrational Spectroscopy, 2012, 60(5): 146. [本文引用:2]
[12] Noda I. Journal of Molecular Structure, 2006, 799(1): 2. [本文引用:3]
[13] Chen H, Lin G, Wang X, et al. Journal of Renewable and Sustainable Energy, 2016, 8(1): 013112. [本文引用:1]
[14] Chen Y, Liu B, Yang H, et al. Fuel, 2014, 137: 41. [本文引用:1]
[15] Keiluweit M, Nico P S, Johnson M G, et al. Environmental Science & Technology, 2010, 44(4): 1247. [本文引用:1]
[16] Yang H, Yan R, Chen H, et al. Fuel, 2007, 86(12): 1781. [本文引用:2]
[17] YU Zi-xuan, YANG Jing-jing, WANG Yu-yan, et al(于子旋, 杨静静, 王语嫣, ). Chinese Journal of Applied Ecology(应用生态学报), 2016, 27(6): 2015. [本文引用:1]
[18] Meng J, Feng X, Dai Z, et al. Environmental Science & Pollution Research International, 2014, 21(11): 7035. [本文引用:1]
[19] Xin S, Yang H, Chen Y, et al. Journal of Analytical and Applied Pyrolysis, 2015, (16): 263. [本文引用:1]
[20] Harvey O R, Herbert B E, Kuo L J, et al. Environmental Science & Technology, 2012, 46(19): 10641. [本文引用:]
[21] XIAO He-lan, SUN Su-qin, ZHOU Qun, et al(肖和兰, 孙素琴, 周群, ). Journal of Atomic and Molecular Physics(原子与分子物理学报), 2003, 20(2): 211. [本文引用:1]
[22] Jiang W, Saxena A, Song B, et al. Langmuir, 2004, 20(26): 11433. [本文引用:1]