郭?，? 魏立新, 耿孝恒, 賈新磊, 霍洪俊, 王鵬華, 黃春峰, 唐善法*

(1.長(zhǎng)江大學(xué)石油工程學(xué)院，武漢 430100；2.濱州學(xué)院化工與安全學(xué)院，濱州 256600；3.常州大學(xué)石油工程學(xué)院，常州 213164)

含油污泥是原油開(kāi)采儲(chǔ)運(yùn)過(guò)程中產(chǎn)生的油、水等的混合廢物，會(huì)影響人體健康和環(huán)境質(zhì)量[1]。因此，必須對(duì)其資源化無(wú)害化處理，使其處理殘余物達(dá)到國(guó)家環(huán)保要求。而現(xiàn)有其他資源化處理方法則均需對(duì)含油污泥進(jìn)行必要的油泥分離或高溫?zé)峤?，造成分離過(guò)程復(fù)雜、成本較高，且存在二次污染風(fēng)險(xiǎn)。

生物電化學(xué)系統(tǒng)是利用微生物催化氧化有機(jī)物并將有機(jī)質(zhì)中化學(xué)能轉(zhuǎn)變成電能的新型系統(tǒng)[2-3]。但受多種因素的影響系統(tǒng)輸出電壓及功率密度始終不高[4-5]。沉積型生物電化學(xué)系統(tǒng)，發(fā)現(xiàn)其啟動(dòng)時(shí)間短，產(chǎn)電效能高，它不需要質(zhì)子交換膜，結(jié)構(gòu)相對(duì)簡(jiǎn)單，將陽(yáng)極插入相對(duì)缺氧的被修復(fù)底物中，陰極放在上層富含氧氣或氧化物的陰極液中，通過(guò)外接電路就可進(jìn)行電子的傳遞[6-7]。同時(shí)，在底物中的陽(yáng)極會(huì)為其提供連續(xù)且穩(wěn)定的受體，加快有機(jī)物的代謝過(guò)程，從而實(shí)現(xiàn)其底泥修復(fù)的功能[8]。為此，將該沉積型生物電化學(xué)系統(tǒng)運(yùn)用于油泥處理，既能利用微生物降解油泥中的石油等有機(jī)物，同時(shí)還能輸出電能，實(shí)現(xiàn)保護(hù)環(huán)境和產(chǎn)電節(jié)能的雙贏。

陽(yáng)極生物膜的電活性會(huì)直接影響生物電化學(xué)系統(tǒng)的產(chǎn)電和降解性能[9]。有研究表明，以丙酸為陽(yáng)極底物陽(yáng)極生物膜優(yōu)勢(shì)菌群為壁厚菌門(mén)，而以丙酸、丁酸和葡萄糖為底物的生物電化學(xué)系統(tǒng)優(yōu)勢(shì)菌群為β-變形菌門(mén)(Proteobacteria)，尤以葡萄糖底物的微生物群落豐富度最高。陽(yáng)極膜微生物組成和豐度會(huì)因陽(yáng)極底物不同而表現(xiàn)出明顯差異[10-11]。

鑒于此，構(gòu)筑了以含油污泥為底物的單室無(wú)膜沉積型生物電化學(xué)系統(tǒng)，并探討了該系統(tǒng)的產(chǎn)電效能及陽(yáng)極膜微生物特征。通過(guò)系統(tǒng)輸出電壓、功率密度曲線、極化曲線及循環(huán)伏安測(cè)試，考察系統(tǒng)的產(chǎn)電性能；同時(shí)，采用生物宏基因組分類測(cè)序分析產(chǎn)電穩(wěn)定陽(yáng)極碳?xì)稚虾统跏己臀勰嗟奈⑸锶郝浣M成及豐度，以期了解以含油污泥生物電化學(xué)系統(tǒng)的陽(yáng)極膜菌群變化特征。

1 實(shí)驗(yàn)

1.1 生物電化學(xué)系統(tǒng)的構(gòu)筑及運(yùn)行

實(shí)驗(yàn)采用的是單室無(wú)膜的沉積型生物電化學(xué)系統(tǒng)，有效體積為2 L(圖1)。含油污泥來(lái)自勝利油田，含油率為40.71%，含水率為10.70%。含油污泥與培養(yǎng)基，以9∶1的比例均勻混合后，取1 L沉積于系統(tǒng)底部，上部為陰極液[12]，體積為1 L。陰陽(yáng)兩電極均采用經(jīng)預(yù)處理的圓形石墨氈(厚度為10 mm，直徑分別為80 mm)，陽(yáng)極埋于含油污泥中，陰極浮于陰極液表面且與空氣相接觸，兩極間距為12.5 cm。陽(yáng)極和陰極通過(guò)銅導(dǎo)線與外電阻(500 Ω)相連，外電阻兩端連接數(shù)據(jù)采集系統(tǒng)，與電腦相連導(dǎo)出數(shù)據(jù)，數(shù)據(jù)采集器每隔10 min記錄一次。裝置于恒溫培養(yǎng)箱[(31±1) ℃]中進(jìn)行，電池運(yùn)行期間定期添加陰極緩沖液。

圖1 含油污泥沉積型生物電化學(xué)系統(tǒng)示意圖Fig.1 Schematic diagram of oily sludge bioelectrochemical system

1.2 產(chǎn)電性能測(cè)試方法

1.2.1 輸出電壓、功率密度及極化曲線

系統(tǒng)產(chǎn)電過(guò)程中，輸出電壓、功率密度等可用于評(píng)價(jià)系統(tǒng)產(chǎn)電效能。系統(tǒng)輸出電壓可由數(shù)據(jù)采集器直接獲得。由P=UI/A，I=U/RA得出功率密度和電流密度，其中，I表示電流密度(mA·m-2)；U表示電壓(mV)；R表示外接電阻(Ω)；P表示陽(yáng)極功率密度(mW·m-2)；A表示陽(yáng)極有效面積(m2)。

在系統(tǒng)發(fā)電穩(wěn)定期間，改變外接電阻大小，變化范圍0～9 000 Ω，當(dāng)0～10 Ω時(shí)，調(diào)節(jié)間隔為1 Ω，10～100 Ω時(shí)，調(diào)節(jié)間隔10 Ω，100～900 Ω時(shí)調(diào)節(jié)間隔為100 Ω，1 000～9 000 Ω時(shí)，調(diào)節(jié)間隔為1 000 Ω；測(cè)量時(shí)每改變一次電阻待MFC穩(wěn)定10 min后記錄數(shù)據(jù)。測(cè)得數(shù)據(jù)計(jì)算電流、電流密度及功率密度并做出系統(tǒng)的功率密度曲線與極化曲線。

1.2.2 循環(huán)伏安曲線測(cè)試

系統(tǒng)的循環(huán)伏安(CV)曲線通過(guò)電化學(xué)工作站(CHI660E，上海辰華儀器有限公司)進(jìn)行測(cè)試。采用三電極體系：對(duì)電極采用鉑絲電極，Ag/AgCl 電極(+0.195 V，標(biāo)準(zhǔn)氫電極電位)為參比電極，對(duì)系統(tǒng)中的陽(yáng)極膜進(jìn)行CV循環(huán)伏安掃描時(shí)，陽(yáng)極為工作電極，電位掃描的范圍為-1.5～1.5 V，掃描速度為10 mV/s。

1.3 陽(yáng)極膜菌群變化特征

1.3.1 DNA提取

采用生物宏基因組分類測(cè)序的方法分析陽(yáng)極碳?xì)旨俺跏加湍嘀械奈⑸锓N群結(jié)構(gòu)。分別提取初始油泥和產(chǎn)電穩(wěn)定期陽(yáng)極碳?xì)治⑸顳NA。

1.3.2 聚合酶鏈?zhǔn)椒磻?yīng)

對(duì)提取的DNA進(jìn)行聚合酶鏈?zhǔn)椒磻?yīng)(polymerase chain reaction, PCR)擴(kuò)增，選取 V4-V5 區(qū)作為擴(kuò)增引物515F(5＇-GTGCCAGCMGCCGCGG-3＇)和907R(5＇-CCGTCAATTCMTTTRAGTTT-3＇)。

1.3.3 操作分類單元的獲得

通過(guò)Illumina MiseqTM平臺(tái)測(cè)得序列，通過(guò)拼接及篩選后，在相似度97%下，將數(shù)據(jù)聚成操作分類單元(operational taxonomic units，OTU)。

1.3.4 生物群落組成及豐度分析

獲得的OTU通過(guò)核糖體數(shù)據(jù)庫(kù)(ribosomal database project, RDP)數(shù)據(jù)庫(kù) (http://rdp.cme.msu.edu/)中的Classifer程序進(jìn)行檢索分類，按照文獻(xiàn)[13]的分析方法分析群落的組成、相對(duì)豐度。

2 結(jié)果與討論

2.1 系統(tǒng)產(chǎn)電性能分析

2.1.1 輸出電壓、功率密度及極化曲線

系統(tǒng)輸出電壓隨時(shí)間變化如圖2所示，系統(tǒng)電壓在0～12 d逐漸上升； 12～20 d，輸出電壓基本穩(wěn)定在300 mV左右，期間最大產(chǎn)電壓為Vmax=320.7 mV；20～30 d，電壓逐漸下降。圖3為12～20 d即系統(tǒng)產(chǎn)電穩(wěn)定期的極化曲線和功率密度曲線。可知，隨著外電阻的減小，輸出的電流密度先增大后減小，且輸出電壓逐漸減小。由功率密度-電流密度曲線(功率密度曲線)可得，當(dāng)電流密度為183.13 mA·m-2時(shí)，輸出的功率密度最大為Pmax=3 353.7 mW·m-2，此時(shí)外接電阻為500 Ω左右。研究表明，電子從陽(yáng)極向陰極的流動(dòng)會(huì)受外阻的控制，因此外阻值會(huì)影響系統(tǒng)產(chǎn)電性能[14]。由極化曲線可知，系統(tǒng)符合電流密度與輸出電壓之間呈直線關(guān)系的規(guī)律(V=-1.036i+324.96，相關(guān)系數(shù)0.98)。根據(jù)電流與電壓間關(guān)系式V=ε-ir(ε為電動(dòng)勢(shì)，V為輸出電壓，i為電流，r為內(nèi)阻)，可得電池電動(dòng)勢(shì)為324.96 mV，電池內(nèi)阻值為509.73 Ω。由此可以得出，系統(tǒng)的內(nèi)阻與功率輸出最大值時(shí)的外阻值(500 Ω)相接近。

圖2 系統(tǒng)輸出電壓隨時(shí)間變化Fig.2 The output voltage change of the system with time

圖3 極化曲線和功率密度曲線Fig.3 Polarization and power density curves

2.1.2 循環(huán)伏安曲線分析

循環(huán)伏安曲線(CV)能確定活性微生物氧化還原反應(yīng)有機(jī)物的中心電位，同時(shí)明確陽(yáng)極膜菌群的催化能力大小、電子傳遞能力及極限電流的大小等[15]。圖4為系統(tǒng)陽(yáng)極生物膜不同產(chǎn)電時(shí)期的CV曲線。可以看出，在產(chǎn)電初期(0 d)，陽(yáng)極生物膜循環(huán)伏安曲線光滑，表明氧化還原反應(yīng)不明顯。在產(chǎn)電穩(wěn)定期(15 d)，CV曲線呈“S”形狀，說(shuō)明產(chǎn)電微生物的催化作用與有機(jī)物的降解具有耦合效果，與之前報(bào)道結(jié)果一致[16]。氧化還原峰分別在0 V和-0.5 V，極限電流值為0.12 A/cm2，表明該陽(yáng)極生物膜具有較強(qiáng)電活性。系統(tǒng)產(chǎn)電30 d時(shí)，CV曲線光滑即無(wú)氧化還原峰。

圖4 陽(yáng)極循環(huán)伏安曲線Fig.4 Cyclic voltammetry curve for anode electrode

2.2 陽(yáng)極膜菌群變化特征

2.2.1 多樣性指數(shù)

初始油泥(YN1)和系統(tǒng)處于產(chǎn)電穩(wěn)定期的碳?xì)?TZ)的多樣性指數(shù)如表1所示，在樣本覆蓋度(coverage)為0.96條件下，YN1的Seq 數(shù)、AEC指數(shù)、Chao1分別為55 020、50 398.64、22 327.10均高于TZ(50 149、31 064.92、15 617.76)，這表明系統(tǒng)處于產(chǎn)電穩(wěn)定期的陽(yáng)極膜比對(duì)應(yīng)初始油泥的微生物菌群多樣性降低，這與前人報(bào)道具有一致性[17]。

表1 油泥及陽(yáng)極膜中微生物菌群多樣性指數(shù)

2.2.2 群落結(jié)構(gòu)組成、相對(duì)豐度及差異性分析

圖5是系統(tǒng)處于產(chǎn)電穩(wěn)定期的碳?xì)?TZ)和原始油泥(YN1)菌群組成及相對(duì)豐度圖(圖5)，并進(jìn)一步對(duì)比兩者菌群存在的差異性(圖6)。由圖5(a)可知， TZ、YN1的菌群組成主要有變形菌門(mén)(Proteobacteria)(46.39%、30.94%)、厚壁菌門(mén)(Firmicutes)(23.84%、45.5%)、綠彎菌門(mén)(Chloroflexi)(12.07%、9.29%)、擬桿菌門(mén)(Bacteroidetes)(6.94%、4.94%)。從門(mén)水平來(lái)看， Proteobacteria為T(mén)Z的優(yōu)勢(shì)菌群，這與以其他有機(jī)質(zhì)為陽(yáng)極底物的生物電化學(xué)系統(tǒng)具有一致性，F(xiàn)irmicutes為YN1優(yōu)勢(shì)菌群，其能在胞外進(jìn)行電子傳遞，是產(chǎn)電微生物群落的重要組成[17]。結(jié)合圖6(a)可知，YN1上Firmicutes相對(duì)豐度增大為T(mén)Z的1.5倍。研究表明，Bacteroidetes 和 Chloroflexi 都包含有多種產(chǎn)電菌[18]。

圖5 初始油泥(YN1)及碳?xì)?TZ)微生物菌群組成及相對(duì)豐度對(duì)比Fig.5 Composition and abundance of microbial communities at different levels on the anode surface and Initial oil sludge

圖6 油泥(YN1)及碳?xì)?TZ)微生物菌群組成差異性圖Fig.6 Difference diagram of microbial flora composition of sludge (YN1) and carbon felt (TZ)

在綱水平上，相對(duì)豐度大于5%的菌群如圖5(b)所示，TZ和YN1中均為梭菌綱(Clostridia)(21.99%，40.04%)、α-變形菌綱(Alphaproteobacteria)(20.12%，14.94%)、厭氧繩菌綱(Anaerolineae)(18.81%，13.33%)和γ-變形菌綱(Gammaproteobacteria)(11.56%，9.03%)。YN1中Clostridia的相對(duì)豐度是TZ的1.8倍，有研究表明，梭菌綱(Clostridias)可通過(guò)胞外傳遞電子[19]。YN1中Alphaproteobacteria、Anaerolineae 和Gammaproteobacteria的相對(duì)豐度均低于TZ。Anaerolineae 綱屬于Chloroflexi門(mén)，是由Yamada 等通過(guò) 16SrRNA 基因分析，構(gòu)建系統(tǒng)進(jìn)化樹(shù)而發(fā)現(xiàn)并命名，最適宜在 37 ℃左右、嚴(yán)格厭氧條件下生存，生長(zhǎng)最適 pH 為 7.0[20]，因此，廣泛存在于系統(tǒng)中。

在目水平上，相對(duì)豐度大于5%的菌群如圖5(c)所示。TZ中有梭菌目(Clostridiales)(22.01%)、厭氧繩菌目(Anaerolineales)(11.72%)、擬桿菌目(Bacteroidales)(5.63%)、紅育菌目(Rhodobacterales)(6.2%)、Rhizobiales(9.99%)、假單胞菌目(Pseudomonadales)(13.26%)。其中，相對(duì)豐度與YN1相差最大的是Pseudomonadales。有研究表明，Pseudomonadales分泌的介體既能提高其自身的產(chǎn)電降解能力，也能為其他菌屬提供電子傳遞的中介體[21]。

在屬水平上[圖5(d)、圖6(b)]，TZ與YN1差異性較大主要集中在Proteiniclasticum、Pseudomonas。Proteiniclasticum在YN1(39.43%)中的相對(duì)豐度是TZ(15.83%)的2.5倍，Proteiniclasticum屬于厚壁菌門(mén)。

3 結(jié)論

(1)構(gòu)筑了以含油污泥為陽(yáng)極底物的單室無(wú)膜沉積型生物電化學(xué)系統(tǒng)，研究發(fā)現(xiàn)系統(tǒng)最大輸出功率、產(chǎn)出電壓為3 353.7 mW·m-2、320.7 mV。陽(yáng)極生物膜具有較強(qiáng)電活性，循環(huán)伏安曲線呈“S”形狀，氧化還原峰分別在-0.5 V和0 V，極限電流值為0.12 A·cm-2。

(2)系統(tǒng)產(chǎn)電穩(wěn)定期陽(yáng)極膜(TZ)相比于相應(yīng)的初始沉積物(YN1)的菌群多樣性降低，且兩者微生物組成及相對(duì)豐度有較大差別，在TZ表面的微生物主要為Proteobacteria，YN1占優(yōu)勢(shì)的微生物為 Firmicutes。該技術(shù)既能高效、清潔資源化利用含油污泥，又能有效避免油泥分離且輸出電能，具有良好的產(chǎn)業(yè)化前景。