工业水处理, 2020, 40(3): 43-48 doi: 10.11894/iwt.2019-0056

试验研究

好氧颗粒污泥处理印染废水的效能及其微生物特征

李黔花,1,2, 李志华,1, 岳秀3, 于广平2, 王慧娟4

Efficiency and microbial characteristics of aerobic granular sludge for textile dyeing wastewater treatment

Li Qianhua,1,2, Li Zhihua,1, Yue Xiu3, Yu Guangping2, Wang Huijuan4

通讯作者: 李志华,博士,教授, E-mail:lizhihua@xauat.edu.cn

收稿日期: 2019-12-27  

基金资助: 国家自然科学基金.  51878539
国家自然科学基金.  51808157
广东省科技计划项目.  2016B0909180369

Received: 2019-12-27  

Fund supported: 国家自然科学基金.  51878539
国家自然科学基金.  51808157
广东省科技计划项目.  2016B0909180369

作者简介 About authors

李黔花(1992-),硕士E-mail:709729228@qq.com , E-mail:709729228@qq.com

摘要

研究了好氧颗粒污泥处理印染废水的效果,并对系统中微生物菌群进行了分析。结果表明:颗粒污泥长时间运行会发生解体,但具有自修复潜力。该体系的COD、氨氮、总磷、色度的平均去除率分别为92.53%、56.65%、95.12%、69.03%。微生物测序结果显示:颗粒污泥解体前后的微生物群落发生了显著变化,体系中的优势菌门为变形菌门,其中磁螺菌属Magnetospirillum对染料降解和COD去除起到了重要的作用。

关键词: 好氧颗粒污泥 ; 印染废水 ; 生物降解 ; 微生物菌群

Abstract

The performance of aerobic granular sludge on textile dyeing wastewater treatment was studied, and the population structure of microorganisms in the system was analyzed. The results showed that granular sludge disintegrated during long-term operation, but had self-repairing potential. The average removal rates of COD, ammonia nitrogen, total phosphorus and chroma were 92.53%, 56.65%, 95.12% and 69.03%, respectively. Microbial sequencing results showed that the microbial community changed significantly before and after the disintegration of granular sludge. The dominant bacteria in the system were Proteobacteria, and the Magnetospirillum played an important role in dye degradation and COD removal.

Keywords: aerobic granular sludge ; textile dyeing wastewater ; biodegradation ; microbial population

PDF (0KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

李黔花, 李志华, 岳秀, 于广平, 王慧娟. 好氧颗粒污泥处理印染废水的效能及其微生物特征. 工业水处理[J], 2020, 40(3): 43-48 doi:10.11894/iwt.2019-0056

Li Qianhua. Efficiency and microbial characteristics of aerobic granular sludge for textile dyeing wastewater treatment. Industrial Water Treatment[J], 2020, 40(3): 43-48 doi:10.11894/iwt.2019-0056

纺织印染废水一直以排放量大、水质复杂、色度高、处理难度大等特点而成为废水治理工艺研究的重点和难点1-3。目前印染废水的处理方法通常是生物法和物化法结合,单独使用任何传统处理方法都不能有效地去除纺织废水中存在的各种污染物。好氧颗粒污泥是一种微生物自凝聚形成的特殊生物膜4,由外向内延伸依次聚集生长好氧、兼性厌氧、厌氧微生物,特殊的层状分布能使好氧颗粒污泥在一个生物池中同时对废水中各类特征污染物进行去除,生物量高,剩余污泥量少,且对有毒有害物质和高负荷废水的冲击耐受力强5-6。因此,好氧颗粒污泥技术是一种潜在的高效废水生物处理技术。由于好氧颗粒污泥技术应用于印染废水处理方面的研究较少,因此笔者在单独的好氧条件下,采用低高径比(3.9)的SBR反应器,接种剩余污泥,在有机负荷(OLR)为3 kg/(m3·d)的模拟印染废水中驯化好氧颗粒污泥,通过测定COD、氨氮(NH4+-N)、总磷(TP)、色度及颗粒特性等指标,评估好氧颗粒污泥对印染废水的处理效能,同时对系统中不同时期的微生物菌群结构进行分析,探究生物转化机理,通过试验进一步论证好氧颗粒污泥技术应用到印染废水处理中的可行性。

1 材料与方法

1.1 实验装置

试验反应器采用SBR的运行方式,主体由圆柱形有机玻璃制成,内径100 mm,有效高度为390 mm,总有效容积为3.0 L,排水率为50%,即每周期进水1.5 L。每天3个运行周期,每周期8 h,HRT为16 h,水力负荷为0.57 m3/(m2·d)。其中进水5 min,曝气450 min,沉淀时间8 min,排水2 min,静置时间15 min。进水泵为蠕动泵(BT100-2J),转速为93.5 r/min;曝气采用鼓风曝气,曝气量为60 L/h,由气体流量计控制;有1#、2#、3#共3个排水口,出水按照顺序依次排出,反应温度为室温。试验装置采用PLC自动控制。

1.2 接种污泥

实验接种污泥取自广州市某污水处理厂二沉池的普通活性污泥。接种前去除污泥表面杂质,将污泥静置倒掉上清液后,剩下的浓缩污泥闷曝48 h,消耗污泥中的有机物,接种污泥的MLSS为8 973 mg/L,SVI为83%,取1.5 L接种污泥倒入SBR反应器中,加入模拟废水至3 L,反应器中的MLSS为4 487 mg/L。

1.3 模拟废水组分

实验用水采用人工配水(各成分单位以mg/L计),以乙酸钠和葡萄糖作为碳源,氯化铵作为氮源,磷酸二氢钾为磷源,根据COD:N:P=100:5:1进行配比。进水基质分为两个阶段,第一阶段(第1天—第23天),以葡萄糖为单一碳源,葡萄糖2 000;第二阶段(第24天—实验结束),以乙酸钠和葡萄糖作为碳源,其中乙酸钠1 280,葡萄糖1 000。其他成分都相同,包括NH4Cl 382,K2HPO4 48,MgSO4·7H2O 120,CaCl2·H2O 160,FeSO4·7H2O 80,NaHCO3 400,活性黑KN-B染料20。每升水中加入0.8 mL的微量元素溶液,微量元素成分为:ZnSO4·7H2O 120,CuSO4·5H2O 30,KI 180,MnCl2·4H2O 120,CoCl2·6H2O 150,H3BO4 150。进水COD 1 704~2 310 mg/L、NH4+-N 85.2~114.8 mg/L、TP 12.94~16.3 mg/L、pH 7.0~7.5。

1.4 分析项目与检测方法

实验中水质测定均采用国家标准方法7,其中NH4+-N测定采用水杨酸试剂法,COD测定采用重铬酸钾法,TP用哈希分光光度法;pH采用便携式pH计(PJBJ-260,上海雷磁公司)测定,污泥微观结构采用扫描电子显微镜观测(S-3000N,Hitachi公司)。

1.5 基于Illumina平台的16S rDNA宏基因组测序

1.5.1 污泥微生物选取与制备

分别选取系统中好氧颗粒污泥部分解体期(第115天)和颗粒修复后期(第128天以后)的污泥进行研究,污泥样品的制备方法为:两个样品均选自系统中的混合污泥,在3 000 r/min条件下离心3 min,去除上清液,取浓缩污泥50 mL,备用。

1.5.2 高通量测序方法

将采集的微生物样品,送至检测机构(广州锐博生物科技有限公司)进行16S rDNA宏基因组测序。

2 结果和讨论

2.1 好氧颗粒污泥颗粒化

本实验在模拟废水中共运行144 d,整个实验过程分为第Ⅰ阶段驯化期(第1天—第48天),第Ⅱ阶段稳定期(第49天—第113天),第Ⅲ阶段解体期(第114天—第125天)与第Ⅳ阶段修复期(第126天—第144天)4个阶段。图1是用数码相机拍摄的颗粒污泥不同时期的形态特征。

图1

图1   不同阶段的颗粒污泥形态


在驯化期的第1天—第23天,采用葡萄糖为单一碳源,系统中有大量的丝状菌繁殖,并且产生部分大而疏松的颗粒污泥,但该污泥结构不稳定,在高剪切力作用下被破碎成絮状,而后丝状菌以颗粒碎片方式黏附形成沉降性差的菌胶团。此时,系统排泥量大,导致污泥浓度下降,为保证系统的运行稳定,将部分排出的污泥进行回流。第24天—实验结束,进水基质改为混合碳源,据J. H. Tay等8研究,由葡萄糖培养的颗粒污泥以丝状菌为主,由乙酸钠培养的颗粒污泥以杆菌为主。为使颗粒污泥变得密实,第24天—实验结束采用乙酸钠和葡萄糖联合培养,引起体系中的丝状菌减少,颗粒污泥逐渐变得密实。从第49天开始,系统中开始出现黑色、粒径2~5 mm的规则密实颗粒,水质处理效果稳定,即认为好氧颗粒污泥进入稳定培养期。经稳定运行113 d后,部分颗粒有解体现象,通过逐步排出解体污泥,微生物又以破碎颗粒为核心聚集生成新颗粒(第Ⅳ阶段修复期),表明好氧颗粒污泥有自修复潜力。

2.2 好氧颗粒污泥的微观结构

还对好氧颗粒污泥的微观结构进行了观察。可以看到,稳定期的颗粒污泥轮廓清晰,结构完整但较为疏松,另外颗粒污泥表面空隙较大,有机物和染料可通过这些孔洞进入AGS内部9,通过孔隙进行物质交换,从而维持好氧颗粒污泥的稳定结构。同时在高倍的扫描电镜下,观察到微生物的聚集状态,主要以杆菌为主,丝状菌缠绕在颗粒中构成骨架10。实验证明,用乙酸钠和葡萄糖联合培养的好氧颗粒污泥,避免了丝状菌的大量繁殖,从而使颗粒结构更稳定。

2.3 印染废水中水质处理效果

印染废水中水质运行情况如图2所示,整个过程分为4个阶段。

图2

图2   反应体系中水质运行情况


图2可看出,第Ⅰ阶段水质波动较大。污泥处于接种初期,微生物量较高,对COD的去除效果明显,第18天时,由于丝状颗粒增大后破碎,导致系统出水SS增高,COD去除率下降;系统中NH4+-N的去除能力有限,在第1天—第18天,由于自养菌与异养菌存在基质竞争关系,异养菌在基质充足的条件下占据优势,自养菌对氨氮的降解逐渐减小,NH4+-N的去除率逐渐下降;出水中TP质量浓度由3.548 mg/L逐渐下降到0.541 mg/L;驯化期的颗粒结构不稳定,因此印染废水中色度的去除波动较大。第42天—第44天时,系统中pH升高到8.74,且有许多小泡沫浮在表面,导致出水COD由170.67 mg/L上升到621.5 mg/L,可能由于曝气过大破坏了系统中CO2-CO32-缓冲平衡体系,CO2逸出使得pH升高,对COD的降解造成冲击,而对NH4+-N和TP的影响较小。

第Ⅱ阶段好氧颗粒污泥形态结构完整,水质处理效果比较稳定。颗粒表面聚集生长好氧微生物,颗粒内部的缺氧/厌氧环境为脱氮除磷提供有利条件,COD的平均去除率为92.53%,NH4+-N的平均去除率达到56.65%,TP的最高去除率达到98.62%,平均去除率为95.12%,色度最大的去除率达到82.92%,平均去除率为69.03%。本研究中COD的去除率高于R. D. G. Franca等11的研究结果,该研究在6 h的厌氧-好氧循环下,COD的去除率为80%~90%,且证明了染料的存在及其潜在的有毒的中间产物,不会对碳底物的摄取产生负面影响。M. Sarvajith等12也得出类似的结论,好氧颗粒污泥上的生物过程(C,N和P的去除)不受偶氮染料负荷的影响。好氧颗粒污泥存在好氧、缺氧、厌氧区,可以在一个系统中对印染废水进行有效降解。本试验中反应器的曝气系统属于上向流,颗粒在反应器底部会逐渐沉积,形成局部厌氧区使污泥老化,在第75天—第80天这段时间,通过排出反应器底部沉积污泥,系统中微生物量减少,COD、NH4+-N和TP的去除受到一定的影响,其中TP的去除率波动最大,最低为82.26%,可能由于排出的污泥中带走部分聚磷菌,使出水TP的浓度变高,去除率下降。

第Ⅲ阶段,由于颗粒粒径增大,颗粒内部会由于传质受阻产生一个厌氧的核心,从而导致颗粒污泥解体13。从整体运行趋势上看,解体期系统对COD和NH4+-N的降解影响较小,出水COD维持在106.33~212.67 mg/L;出水中NH4+-N从40 mg/L升高到46.25 mg/L;而在好氧的反应系统中,颗粒内部的厌氧环境被破坏,TP和染料的色度脱除情况明显受到冲击,TP去除率从97.08%下降到90.73%;色度的去除率由82.92%降为60.80%。

第Ⅳ阶段,微生物以解体污泥为支架,重新聚合生成了新颗粒。系统对COD、NH4+-N、TP、色度的平均去除率分别为:92.59%、61.58%、96.77%、74.49%。这一结果体现好氧颗粒污泥法对染料废水稳定去除的技术优势。

2.4 微生物菌群的宏基因组16S rDNA测序分析

2.4.1 微生物Alpha多样性分析

Alpha多样性是用于分析样品内的微生物群落多样性。Alpha多样性有多种衡量指标,如Chao1、Shannon、Simpson等,这些指数数值越大,表明物种丰度越大14。通常用稀释曲线来评价测序量是否足以覆盖所有类群,并间接反映样品中物种的丰富程度。当曲线趋于平缓或者达到平台期时则认为测序深度已基本覆盖样品中的所有物种。

本研究选取2个污泥样品,分别为系统运行第115天(1号)和第128天(2号)的污泥,代表着解体期与颗粒修复期。微生物Alpha多样性相关的各项指标如表1所示。

表1   微生物样品的Alpha多样性

样本Chaol指数Shannon指数Simpson指数OUT数目Tags数目样本覆盖度/%
1号806.144.280.7967573 72899.59
2号561.982.470.5341964 18399.67

新窗口打开| 下载CSV


结果显示本研究中样品曲线基本趋于平缓,测序结果占整个基因组的99%,样品覆盖度高,证明样品测序量基本能够反映出样品的物种丰度。样品2即修复期颗粒的各项指数都低于解体期的样品1,说明好氧颗粒污泥经解体修复后,物种的多样性减少。

2.4.2 微生物群落结构特征

对宏基因组16S rDNA序列在微生物门分类层面上的比对分析,发现在模拟印染废水培养的好氧颗粒污泥系统中,微生物菌群结构比较单一,但各物种的丰度较高,富集的优势菌群对系统运行起到了重要的作用。在两个污泥样品中,除2.5%的未识别菌种外,微生物共检测出19种菌门,解体期时物种丰度大于1%的菌门有7种,包括变形菌门Proteobacteria(81.19%)、拟杆菌门Bacteroidetes(9.63%)、厚壁菌门Firmicutes(2.08%)、疣微菌门Verrucomicrobia(1.98%)、TM7(1.31%)、放线菌门Actinobacteria(1.18%)、绿弯菌门Chloroflexi(1%)。其中变形杆菌Proteobacteria、拟杆菌Bacteroidetes和厚壁菌门Firmicutes,是偶氮染料中常见的微生物15。变形杆菌门Proteobacteria为系统中的绝对优势菌门,Y. M. Kolekar等16研究显示α-,β-和γ-proteobacterial菌群参与偶氮染料脱色。颗粒污泥修复后,微生物菌群发生较大的变化,物种丰度大于1%的菌门有2种,变形杆菌Proteo-bacteria(90.83%)丰度比解体时期提高了9.64%,而拟杆菌Bacteroidetes(5.99%)相对解体期下降了3.64%,表明变形菌门Proteobacteria对颗粒稳定性有重要意义。其次,TM7在污泥解体时丰度为1.31%,而在颗粒修复后这个菌门消失,P. Hugenholtz等17研究发现,TM7会引起污泥膨胀,因此本研究中颗粒污泥解体可能与TM7菌门有关。

从微生物的属分类水平进行分析,结果如表2所示。在属分类水平上解体期污泥与颗粒修复后期的污泥存在显著差异,变形菌门Proteobacteria的α-proteobacteria纲下属红螺菌科Rhodospirillaceae的Magnetospirillum菌属为系统的绝对优势菌属。Bin Hou等18研究显示,红螺菌科Rhodospirillaceae参与偶氮染料的脱色与生物转化,在本研究中,红螺菌科的相对丰度远高于研究中的数量,相对丰度从46.36%提升到72.92%,说明解体颗粒修复后,对印染废水中有机物和色度的去除有良好的效果,COD和色度的平均去除率分别为92.59%和74.49%。在解体期系统中的脱氮除磷功能主要由DechloromonasAcinetobacter19-20完成,这两种菌的相对丰度分别为9.59%和2.37%,王启镔等〔21〕研究表明,在较高的有机负荷〔>0.91 kg/(m3·d)〕条件下有利于丝状细菌Acinetobacter的生长。随着颗粒逐渐修复,DechloromonasAcinetobacter菌群减少,而β-proteobacteria纲中的红环菌科Rhodocyclaceae的KD1-23型菌属在颗粒修复后相对丰度从0.06%提升到2.69%,菌群结构发生了变化。

表2   反应体系内微生物的菌群演化

微生物菌属解体期丰度/%修复期丰度/%
Bacteroidetes;c_〔Saprospirae〕;o_〔Saprospirales〕;f_Chitinophagaceae;g_Unnamed1.010.31
Bacteroidetes;c_〔Saprospirae〕;o_〔Saprospirales〕;f_Saprospiraceae;g_Unnamed3.511.16
Proteobacteria;c_Alphaproteobacteria;o_BD7-3;f_Unnamed;g_Unnamed4.431.61
Proteobacteria;c_Alphaproteobacteria;o_Rhizobiales;f_Hyphomicrobiaceae;g_Unnamed1.410.47
Proteobacteria;c_Alphaproteobacteria;o_Rhizobiales;f_Hyphomicrobiaceae;g_Hyphomicrobium1.280.4
Proteobacteria;c_Alphaproteobacteria;o_Rhizobiales;f_Methylocystaceae;g_Pleomorphomonas7.17.74
Proteobacteria;c_Alphaproteobacteria;o_Rhodospirillales;f_Rhodospirillaceae;g_Magnetospirillum46.3672.92
Proteobacteria;c_Betaproteobacteria;o_Unnamed;f_Unnamed;g_Unnamed1.17
Proteobacteria;c_Betaproteobacteria;o_Rhodocyclales;f_Rhodocyclaceae;g_Dechloromonas9.590.29
Proteobacteria;c_Betaproteobacteria;o_Rhodocyclales;f_Rhodocyclaceae;g_KD1-232.69
Proteobacteria;c_Gammaproteobacteria;o_Pseudomonadales;f_Moraxellaceae;g_Acinetobacter2.370.54
TM7;c_TM7-3;o_;f_Unnamed;g_Unnamed1.240.14

新窗口打开| 下载CSV


3 结论

(1)好氧颗粒污泥于第49天在模拟印染废水中驯化成功,3 kg/(m3·d)条件下培养的颗粒污泥粒径较大(2~5 mm),长时间运行会使颗粒发生解体,但好氧颗粒污泥表现出自修复潜力。

(2)该体系的COD、氨氮、总磷、色度的平均去除率分别为92.53%、56.65%、95.12%、69.03%,好氧颗粒污泥能在一个污泥系统中,对印染废水进行有效降解。

(3)本研究中的微生物优势菌门是变形菌门Proteobacteria,其中α-Proteobacteria纲中的磁螺菌属Magnetospirillum对染料降解和COD去除起到了重要的作用。

参考文献

申洋洋, 刘锐, 徐灿灿, .

印染及染料行业废水生物处理系统中的AOX污染研究

[J]. 环境科学, 2015, 36 (9): 3304- 3310.

URL     [本文引用: 1]

陈红, 李响, 薛罡, .

当前印染废水治理中的关键问题

[J]. 工业水处理, 2015, 35 (10): 16- 19.

DOI:10.11894/1005-829x.2015.35(10).016     

王俊峰, 赵英武, 毛燕芳.

我国印染废水处理概况及研究进展

[J]. 中国环保产业, 2012, (4): 30- 33.

DOI:10.3969/j.issn.1006-5377.2012.04.009      [本文引用: 1]

Jiang H L , Tay J H , Maszenan A M , et al.

Bacterial diversity and function of aerobic granules engineered in a sequencing batch reactor for phenol degradation

[J]. Applied and Environmental Microbiology, 2004, 70 (11): 6767- 6775.

DOI:10.1128/AEM.70.11.6767-6775.2004      [本文引用: 1]

苏海佳, 王陆玺, 邓爽, .

好氧颗粒污泥技术及研究进展

[J]. 化工进展, 2016, 35 (6): 1914- 1922.

URL     [本文引用: 1]

Zheng Yuming , Yu Hanqing , Liu Shuangjiang , et al.

Formation and instability of aerobic granules under high organic loading conditions

[J]. Chemosphere, 2006, 63 (10): 1791- 1800.

DOI:10.1016/j.chemosphere.2005.08.055      [本文引用: 1]

国家环保总局.水和废水监测分析方法[M]. 4版.北京:中国环境科学出版社, 2002: 200-284.

[本文引用: 1]

Tay J H , Liu Q S , Liu Y .

Characteristics of aerobic granules grown on glucose and acetate in sequential aerobic sludge blanket reactors

[J]. Environmental Technology, 2002, 23 (8): 931- 936.

DOI:10.1080/09593332308618363      [本文引用: 1]

王建龙, 张子健, 吴伟伟.

好氧颗粒污泥的研究进展

[J]. 环境科学学报, 2009, (3): 449- 473.

URL     [本文引用: 1]

Beun J J , Hendriks A , van Loosdrecht M C M , et al.

Aerobic granulation in a sequencing batch reactor

[J]. Water Research, 1999, 33 (10): 2283- 2290.

DOI:10.1016/S0043-1354(98)00463-1      [本文引用: 1]

Franca R D G , Vieira A , Mata A M T , et al.

Effect of an azo dye on the performance of an aerobic granular sludge sequencing batch reactor treating a simulated textile wastewater

[J]. Water Research, 2015, 85, 327- 336.

DOI:10.1016/j.watres.2015.08.043      [本文引用: 1]

Sarvajith M , Reddy G K K , Nancharaiah Y V .

Textile dye biodecolourization and ammonium removal over nitrite in aerobic granular sludge sequencing batch reactors

[J]. Journal of Hazardous Materials, 2018, 342, 536- 543.

DOI:10.1016/j.jhazmat.2017.08.064      [本文引用: 1]

Adav S S , Lee D J , Lai J Y .

Proteolytic activity in stored aerobic granular sludge and structural integrity

[J]. Bioresource Technology, 2009, 100 (1): 68- 73.

DOI:10.1016/j.biortech.2008.05.045      [本文引用: 1]

岳秀, 唐嘉丽, 于广平, .

双氧水协同生化法强化处理印染废水

[J]. 环境科学, 2017, 38 (9): 3769- 3780.

URL     [本文引用: 1]

Chaudhari A U , Paul D , Dhotre D , et al.

Effective biotransformation and detoxification of anthraquinone dye reactive blue 4 by using aerobic bacterial granules

[J]. Water Research, 2017, 122, 603- 613.

DOI:10.1016/j.watres.2017.06.005      [本文引用: 1]

Kolekar Y M , Nemade H N , Markad V L , et al.

Decolorization and biodegradation of azo dye, reactive blue 59 by aerobic granules

[J]. Bioresour. Technol., 2012, 104, 818- 822.

DOI:10.1016/j.biortech.2011.11.046      [本文引用: 1]

Hugenholtz P , Tyson G W , Webb R I , et al.

Investigation of candidate division TM7, a recently recognized major lineage of the domain bacteria with no known pure-culture representatives

[J]. Applied and Environmental Microbiology, 2001, 67 (1): 411- 419.

DOI:10.1128/AEM.67.1.411-419.2001      [本文引用: 1]

Hou Bin , Sun Jian , Hu Yongyou .

Effect of enrichment procedures on performance and microbial diversity of microbial fuel cell for Congo red decolorization and electricity generation

[J]. Applied Microbiology & Biotechnology, 2011, 90 (4): 1563- 1572.

URL     [本文引用: 1]

Ferrera I , Sánchez O .

Insights into microbial diversity in wastewater treatment systems:How far have we come?

[J]. Biotechnology Advan ces, 2016, 34 (5): 790- 802.

DOI:10.1016/j.biotechadv.2016.04.003      [本文引用: 1]

王启镔, 苑泉, 宫徽, .

SBR系统在低浓度污水条件下培养的好氧颗粒污泥特性及微生物分析

[J]. 环境工程学报, 2018, 12 (11): 3043- 3052.

DOI:10.12030/j.cjee.201805080      [本文引用: 1]

/