張皓輝，史俊祥，姜永海，賈永鋒，李秀金

1. 北京化工大學(xué)化學(xué)工程學(xué)院，北京 100029

2. 中國環(huán)境科學(xué)研究院，國家環(huán)境保護(hù)地下水污染模擬與控制重點(diǎn)實(shí)驗(yàn)室，北京 100012

近年來，被譽(yù)為“工業(yè)血液”的石油工業(yè)迅猛發(fā)展，為經(jīng)濟(jì)社會發(fā)展做出了重大貢獻(xiàn). 但石油產(chǎn)品由于開采、運(yùn)輸、貯存和使用過程中的不規(guī)范操作和意外情況的發(fā)生不斷泄漏進(jìn)入自然環(huán)境，嚴(yán)重污染土壤、地表水和地下水，導(dǎo)致健康風(fēng)險升高，對人體健康產(chǎn)生不可忽視的影響，石油污染已經(jīng)成為全球范圍內(nèi)的環(huán)境問題. 石油烴污染物具有毒性、疏水性，且遷移能力較強(qiáng)，被許多國家列為優(yōu)先控制污染物. 它們不僅會吸附在土壤介質(zhì)中，破壞土壤有機(jī)質(zhì)組成和結(jié)構(gòu)的平衡，同時還會進(jìn)入地表水和地下水，造成水環(huán)境質(zhì)量惡化，毒害水生生物. 大量具有“三致”危害的石油烴類可能通過食物鏈的層級傳導(dǎo)、呼吸作用和皮膚接觸等途徑進(jìn)入人體，危害人體健康. 因此，需要進(jìn)行場地地下水石油烴健康風(fēng)險評估.

土著微生物會通過自身代謝作用降解石油烴，因此土著微生物群落的生物降解通常被認(rèn)為是石油烴自然衰減的主要機(jī)制，可以有效控制其健康風(fēng)險. 微生物群落長時間暴露在石油烴中會改變它們的活動、結(jié)構(gòu)和多樣性，由于不同污染場地污染物種類和環(huán)境因子的差異，土著微生物種類也會有所不同. 石油烴的生物可利用度和化學(xué)性質(zhì)的復(fù)雜程度會影響微生物的降解能力. 歐美國家在20世紀(jì)80年代逐漸建立了針對各自國家實(shí)際情況的場地風(fēng)險評估體系. 中國科學(xué)院南京土壤研究所陳夢舫研究團(tuán)隊(duì)針對我國污染場地實(shí)際情況，自主開發(fā)出我國首款污染場地健康與環(huán)境風(fēng)險評估軟件HERA (Health and Environmental Risk Assessment). HERA包含多種模型和美國環(huán)境保護(hù)局發(fā)布的土壤污染物毒理性參數(shù)數(shù)據(jù)庫，評估結(jié)論對場地是否會給人體健康帶來潛在風(fēng)險具有重要的理論和現(xiàn)實(shí)意義.

該文以天津市某污染場地為研究對象，選取石油烴為關(guān)注污染物，利用HERA軟件，開展地下水石油烴健康風(fēng)險評估，并進(jìn)行地下水微生物群落特征分析，以期為污染場地土壤和地下水中石油烴的存在水平、分布特征和健康風(fēng)險評估，以及評價微生物對石油烴的降解能力提供參考.

1 材料與方法

1.1 場地概況

場地位于天津市海河沿岸，20世紀(jì)90年代作為農(nóng)業(yè)用地受化肥及農(nóng)藥施用污染，21世紀(jì)初引進(jìn)電子廠、配件加工廠、精密電子公司等，原料儲存運(yùn)輸過程及生產(chǎn)加工過程中對場地造成污染. 通過污染識別判斷，地塊內(nèi)潛在污染源主要包括氯代烴類和苯系物等揮發(fā)性有機(jī)物，總石油烴和多環(huán)芳烴等半揮發(fā)性有機(jī)物，汞、砷、銅等重金屬，以及有機(jī)氯農(nóng)藥等. 目前，各工廠已拆除，該場地基本為堆土，僅場地南側(cè)有一食品加工廠. 場地周邊有人居住，定義為敏感用地，因此可能通過經(jīng)口攝入、皮膚接觸和吸入地下水蒸氣等暴露途徑對人體健康產(chǎn)生一定影響.

場地包氣帶厚度為0.45～4.76 m，巖性以素填土、雜填土為主. 包氣帶在場地內(nèi)連續(xù)穩(wěn)定存在，防污性能較弱. 場地內(nèi)潛水含水層地下水水位埋深在0.45～3.76 m之間，平均水位埋深為2.51 m，水位標(biāo)高在0.25～1.37 m 之間，平均水位標(biāo)高為0.856 m. 場地地下水由西南向東北流動，工作區(qū)水力坡度約為0.56‰.

1.2 樣品采集

圖 1 研究區(qū)監(jiān)測井點(diǎn)位分布Fig.1 The location of monitoring wells in study area

地塊采用專業(yè)判斷布點(diǎn)法，共布設(shè)地下水監(jiān)測井5個(A1、A2、A3、A4、A5，見圖1)，以污染物的主要泄漏點(diǎn)A5監(jiān)測井為中心，按地下水流向分布布設(shè)上游A3監(jiān)測井、下游A4監(jiān)測井和側(cè)向A1、A2監(jiān)測井. 按照《地下水環(huán)境監(jiān)測技術(shù)規(guī)范》(HJ 164-2020)采集地下水，放入帶有干冰的保溫箱(≤4 ℃)運(yùn)回實(shí)驗(yàn)室，進(jìn)行石油烴濃度分析和微生物群落分析.

1.3 石油烴濃度分析

根據(jù)《水質(zhì)可萃取性石油烴(C～C)的測定氣相色譜法》(HJ 894-2017)檢測石油烴濃度. 地下水經(jīng)二氯甲烷進(jìn)行萃取、濃縮、凈化、定容后，用氣相色譜-火焰離子化檢測器(GC-FID)(Agilent 7890B，安捷倫，美國)進(jìn)行純化和分析. 進(jìn)樣溫度設(shè)置320 ℃，色譜柱流速為2.0 mL/min. 柱箱保持1 min初始溫度60 ℃，以8 ℃/min升至290 ℃，再以30 ℃/min升至320 ℃后保持7 min. FID檢測器設(shè)置為330 ℃，氫氣流量為40.0 mL/min，空氣流量為350.0 mL/min，尾吹氣流量為30.0 mL/min. 進(jìn)樣不分流，進(jìn)樣體積為1.0 μL，0.75 min后以30:1比例分流. 根據(jù)色譜圖組分保留時間，計算目標(biāo)化合物總峰面積，由外標(biāo)法得出目標(biāo)化合物總濃度. 石油烴餾分的碳數(shù)與沸點(diǎn)有關(guān)，沸點(diǎn)與色譜圖的保留時間有關(guān)，根據(jù)保留時間將石油烴餾分分段積分，石油烴分段為C～C、C～C、C～C、C～C.

1.4 微生物高通量測序

將2 L地下水樣中微生物用0.22 μm濾膜濾出，至濾膜表面有可見覆蓋物，使用FastDNASpin Kit for Soil (MP Biomedicals, 美國)提取總DNA，使用引物338F (5′-ACTCCTACGGGAGGCAGCAG-3′ )和806R(5′-GGACTACHVGGGTWTCTAAT-3′ )對16S rRNA基因V3～V4可變區(qū)進(jìn)行PCR擴(kuò)增. 利用Illumina公司的Miseq PE300平臺進(jìn)行測序. 使用UPARSE軟件(http://drive5.com/uparse，version 7.1)，根據(jù)97%的相似度對序列進(jìn)行OTU聚類并剔除嵌合體. 利用RDP classifier(http://rdp.cme.msu.edu，version 2.2)對每條序列進(jìn)行物種分類注釋，比對Silva 16S rRNA數(shù)據(jù)庫(v138)，設(shè)置比對閾值為70%.

1.5 健康風(fēng)險評估

利用HERA軟件進(jìn)行2個層次的健康風(fēng)險評估：首先選擇合適的場地模型，設(shè)置參數(shù)，計算原場地下水中石油烴污染物的篩選值；然后，結(jié)合暴露途徑，調(diào)整場地概念模型，計算石油烴污染物在原場條件下的危害商，對地下水石油烴進(jìn)行原場健康風(fēng)險評估.

2 結(jié)果與討論

2.1 地下水石油烴分布特征

各監(jiān)測井地下水石油烴濃度如表1所示. 由表1可見，A3監(jiān)測井石油烴濃度(0.06 mg/L)最低，A5監(jiān)測井石油烴濃度(6.61 mg/L)最高，A1～A4監(jiān)測井石油烴濃度差別不顯著，A5監(jiān)測井石油烴濃度約是其他4個監(jiān)測井的60～100倍. A3監(jiān)測井石油烴濃度最低可能是距離污染源較遠(yuǎn)所致，污染物濃度隨遷移距離增加而降低，且地下水由西南向東北流動，A3監(jiān)測井位于地下水流向上游區(qū)，所以石油烴濃度最低.A5監(jiān)測井位于多污染源交匯處，受多個污染物泄漏點(diǎn)影響，因此石油烴濃度最高. 5個監(jiān)測井石油烴分布呈由低碳數(shù)到高碳數(shù)逐漸升高的趨勢(見表1). 已有研究表明，石油烴碳數(shù)越多，結(jié)構(gòu)越復(fù)雜，越不容易降解，對人體危害越大.

表 1 地下水石油烴濃度分布Table 1 Groundwater petroleum hydrocarbon concentration

2.2 地下水石油烴健康風(fēng)險評估

篩選值評估

確定污染物和暴露途徑. 根據(jù)場地污染物分析結(jié)果，確定石油烴(C～C)為目標(biāo)污染物，影響對象為場地周邊范圍內(nèi)的成人和兒童. 具體暴露途徑包括土壤中經(jīng)口攝入、皮膚接觸、吸入室內(nèi)土壤顆粒物、吸入室外土壤顆粒物、吸入表層土壤室外蒸汽、吸入下層土壤室外蒸汽、吸入下層土壤室內(nèi)蒸汽和土壤淋濾到地下水、地下水中吸入地下水室外蒸汽和吸入地下水室內(nèi)蒸汽.

模型構(gòu)建. 選用的模型有C-RAG、ASTM、Johnson-Ettinger等. 模型參數(shù)從場地信息調(diào)查報告中獲取，主要包括地下水中污染物濃度、地下水埋深、地下水混合區(qū)厚度、含水層有效孔隙度等水文地質(zhì)參數(shù)，以及暴露期、暴露頻率等受體暴露參數(shù)，其他參數(shù)采用《污染場地風(fēng)險評估技術(shù)導(dǎo)則》(HJ 25.3-2019)中的推薦值.

篩選值分析. 目前頒布的《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T 14848-2017)中，沒有石油烴類限值說明，《地表水質(zhì)量標(biāo)準(zhǔn)》(GB 3838-2002)Ⅳ類標(biāo)準(zhǔn)中石油類限值為0.5 mg/L. 使用HERA軟件綜合分析石油烴分段毒性和理化性質(zhì)，對各分段石油烴計算篩選值(見表2). 經(jīng)過比較，選用HERA結(jié)果作為篩選值. 將各點(diǎn)位石油烴濃度數(shù)據(jù)與篩選值進(jìn)行對比發(fā)現(xiàn)，A1～A4監(jiān)測井的分段石油烴濃度均低于篩選值，所以健康風(fēng)險可以忽略；A5監(jiān)測井中C～C濃度低于篩選值，但C～C、C～C、C～C濃度均高于篩選值，因此A5監(jiān)測井的石油烴存在一定的環(huán)境健康風(fēng)險，需要利用HERA軟件進(jìn)行第二層次分析，并計算出石油烴污染物在所選暴露途徑下的危害商，以確定石油烴的健康風(fēng)險.

表 2 健康風(fēng)險評估HERA分析結(jié)果Table 2 HERA analysis results of health risk assessment

危害商評估

調(diào)整模型. 進(jìn)行危害商計算時，需調(diào)整模型參數(shù)、修改暴露途徑、更新場地概念模型、輸入污染物濃度，對場地地下水石油烴污染物的致癌風(fēng)險與非致癌風(fēng)險進(jìn)行分析.

危害商分析. 《建設(shè)用地土壤污染風(fēng)險評估技術(shù)導(dǎo)則》(HJ 25.3-2019)規(guī)定，危害商是污染物每日攝入劑量與參考劑量的比值，用于表征人體經(jīng)單一途徑暴露于非致癌污染物而受到危害的水平，單一污染物的總危害商由各途徑危害商加和得到，不同暴露途徑的危害商按照HJ 25.3-2019附錄中公式〔A.13～A.20，見式(1)～(6)〕由HERA軟件自動計算得出. 對受體人群的健康影響，包括致癌風(fēng)險和非致癌風(fēng)險，根據(jù)計算結(jié)果發(fā)現(xiàn)，此次風(fēng)險評估可不考慮選定暴露途徑下的致癌風(fēng)險，非致癌危害商如表2所示，單一污染物的非致癌可接受危害商為1. 由表2可見，非致癌風(fēng)險隨碳鏈長度的增加而增高. A1～A4監(jiān)測井的總石油烴和分段石油烴的危害商均小于1，健康風(fēng)險為可接受水平；A5監(jiān)測井C～C的危害商小于1，但是C～C、C～C、C～C的危害商均高于1，總石油烴危害商為44.990，是可接受值的45倍，表示A5監(jiān)測井的風(fēng)險較高，對人體健康可能造成損害.

式中：I OVER為吸入室外空氣中來自地下水的氣態(tài)污染物對應(yīng)的地下水暴露量(致癌效應(yīng))，L/(kg·d)；VF為地下水中污染物擴(kuò)散進(jìn)入室外空氣的揮發(fā)因子，L/m； D AIR為兒童空氣呼吸量，m/d； DAIR為成人空氣呼吸量，m/d； EFO為兒童的室外暴露頻率，d/a； EFO為 成人的室外暴露頻率，d/a； ED為兒童暴露期，a； B W為 兒童體質(zhì)量，kg； E D為成人暴露期，a；BW為成人體質(zhì)量，kg；A T為致癌效應(yīng)平均時間，d.

式中：I OVER為吸入室外空氣中來自地下水的氣態(tài)污染物對應(yīng)的地下水暴露量(非致癌效應(yīng))，L/(kg·d)；AT為非致癌效應(yīng)平均時間，d.

式中： I IVER為吸入室內(nèi)空氣中來自地下水的氣態(tài)污染物對應(yīng)的地下水暴露量(致癌效應(yīng))，L/(kg·d)；VF為地下水中污染物擴(kuò)散進(jìn)入室內(nèi)空氣的揮發(fā)因子，L/m； EFI為 兒童的室內(nèi)暴露頻率，d/a； EFI為成人的室內(nèi)暴露頻率，d/a.

式中， I IVER為吸入室內(nèi)空氣中來自地下水的氣態(tài)污染物對應(yīng)的地下水暴露量(非致癌效應(yīng))，L/(kg·d).

式中：C GWER為飲用受影響地下水對應(yīng)的地下水的暴露量(致癌效應(yīng))，L/(kg·d)； G WCR為兒童每日飲水量，L/d；G WCR為 成人每日飲水量，L/d； E F為兒童暴露頻率，d/a； E F為成人暴露頻率，d/a.

式中，C GWER為飲用受影響地下水對應(yīng)的地下水的暴露量(非致癌效應(yīng))，L/(kg·d).

2.3 地下水微生物分析

微生物門水平和屬水平分析

圖 2 各監(jiān)測井微生物群落門水平的絕對豐度Fig.2 Absolute abundance of microbial community phylum horizontal of all monitoring wells heat map

微生物群落結(jié)構(gòu)的變化可以反映微生物的分布情況. 基于16SrRNA的微生物群落高通量測序結(jié)果，探討了地下水微生物種類豐富程度與石油烴健康風(fēng)險之間的關(guān)系. A3監(jiān)測井地下水樣品未用0.22 μm濾膜濾出足量的微生物，因此缺失其微生物數(shù)據(jù)，其他4個監(jiān)測井地下水中細(xì)菌門分類水平上優(yōu)勢菌群的絕對豐度見圖2. 由圖2可見，A1監(jiān)測井中細(xì)菌門分類水平上絕對豐度由高到低依次為Patescibacteria、變形菌門、厚壁菌門、放線菌門和擬桿菌門，A2監(jiān)測井中依次為厚壁菌門、變形菌門和擬桿菌門，A4監(jiān)測井中依次為變形菌門、Dependentiae、Patescibacteria和厚壁菌門，A5監(jiān)測井中依次為厚壁菌門、擬桿菌門、Patescibacteria和變形菌門. A1～A5(除A3外)監(jiān)測井中均含有豐富的變形菌門和厚壁菌門，與已有研究發(fā)現(xiàn)變形菌門和厚壁菌門是石油烴降解優(yōu)勢類群的結(jié)果一致. 擬桿菌門是一種發(fā)酵和硫酸鹽還原微生物，能降解各種碳?xì)浠衔锊a(chǎn)生短鏈有機(jī)酸、乙酸和H. 4個監(jiān)測井中其他種類細(xì)菌含量較少，可能與這些菌類的耐受性較低有關(guān).

4個監(jiān)測井(A1、A2、A4、A5)地下水中細(xì)菌屬分類水平上優(yōu)勢菌群的絕對豐度位如圖3所示. 由圖3可見：A1和A4監(jiān)測井中優(yōu)勢菌屬是未知分類屬；A2監(jiān)測井中優(yōu)勢菌屬為；A5監(jiān)測井含有的分枝桿菌屬()和湖水水桿狀菌屬()，二者都有降解石油烴的功能. A1和A4監(jiān)測井細(xì)菌屬水平種類較多，A2和A5監(jiān)測井細(xì)菌屬水平種類較少，且4個監(jiān)測井的優(yōu)勢菌屬有所差異，相關(guān)性不明顯，可能與各監(jiān)測井中其他種類污染物有關(guān). 細(xì)菌的差異性表明，石油烴對具有石油烴降解能力的細(xì)菌具有選擇性. 微生物群落的演替是環(huán)境選擇的結(jié)果，但微生物群落的具體演替規(guī)律還有待進(jìn)一步研究.

細(xì)菌α多樣性分析

不同監(jiān)測井中微生物的α多樣性指數(shù)如表3所示. Chao指數(shù)和Shannon-Wiener指數(shù)分別指示物種豐富度和多樣性，數(shù)值越高表明樣品的物種豐富度和多樣性越高. 用Simpson指數(shù)評價物種均勻性，當(dāng)Simpson指數(shù)越接近1時，說明微生物群落分布越均勻. A5監(jiān)測井Chao指數(shù)和Shannon-Wiener指數(shù)最高(見表3)，說明該點(diǎn)物種豐富度最高； Simpson指數(shù)最低，說明該點(diǎn)微生物群落分布種類多樣. 而A1、A2、A4監(jiān)測井的Simpson指數(shù)更接近1(見表3)，說明微生物群落均一性高，但Chao指數(shù)和Shannon-Wiener指數(shù)均較低，說明A2、A4監(jiān)測井物種豐富度均較低. 因此，一定范圍內(nèi)石油烴濃度的增加，利于石油烴降解功能微生物的生長，使降解菌豐富度升高.微生物群落結(jié)構(gòu)的多樣性和豐富度與石油烴的降解有關(guān)，A5監(jiān)測井的微生物豐富度最高，也對應(yīng)A5監(jiān)測井的石油烴降解能力較強(qiáng)(見圖4). 4個監(jiān)測井含有豐富的變形菌門和厚壁菌門(見圖4)，二者相對豐度占細(xì)菌群落總數(shù)的50%以上，有利于石油烴的生物降解. 研究表明，油水界面上微生物相互作用具有混合表面潤濕性，油水相互作用會受到微生物生物量的顯著影響，并且總生物量與油水乳化液質(zhì)量成正比. 污染環(huán)境中物種豐富度和多樣性水平的提高有助于微生物種群特殊功能的實(shí)現(xiàn)，利于污染物衰減.

圖 3 各監(jiān)測井微生物群落屬水平的絕對豐度Fig.3 Absolute abundance of microbial community gene horizontal of all monitoring wells heat map

表 3 監(jiān)測井中微生物的α多樣性指數(shù)Table 3 The α diversity index of microorganism of all monitoring wells

2.4 微生物生化作用對石油烴健康風(fēng)險的影響

風(fēng)險管控的原理指采取各種措施和方法削減、清除、固定污染源，保護(hù)風(fēng)險受體，控制其健康風(fēng)險. 本地微生物菌群在石油烴的生物降解和生物修復(fù)過程中發(fā)揮了重要作用. 石油烴污染使環(huán)境發(fā)生了微妙的變化，表現(xiàn)出對石油烴降解菌的選擇性，可能減少一些對石油烴不耐受微生物的數(shù)量，增加具有石油烴降解功能微生物數(shù)量，群落組成向石油烴優(yōu)勢降解菌種群非均勻轉(zhuǎn)變. 這種微生物群落動態(tài)變化可以提高對石油烴的生物降解性能，從而提高自然衰減效率，所以適量濃度的石油烴會富集石油烴降解功能菌，從而正向刺激石油烴的自然衰減效率.各種細(xì)菌、真菌參與石油烴降解，可通過自身代謝作用，將石油烴污染物質(zhì)降解成HO、CO等無毒無害物質(zhì). 一些細(xì)菌可以產(chǎn)生生物表面活性劑來正向輔助降解石油烴，生物表面活性劑可以乳化石油，使大分子油團(tuán)分散成小分子油滴，包裹小油滴進(jìn)入水相，增加石油烴與微生物細(xì)胞膜的接觸面積，從而提高降解效率. Varjani等報道了假單胞菌產(chǎn)生的生物表面活性劑(鼠李糖脂)能有效去除石油烴污染物. 這些細(xì)菌是有效的石油烴修復(fù)劑. 自然環(huán)境中石油烴降解功能菌僅占群落總數(shù)的不到1%，石油烴污染環(huán)境中石油烴降解功能菌比例可能在10%左右. 多數(shù)石油烴降解菌都有負(fù)責(zé)降解目標(biāo)底物的酶，由于缺少某些酶，很少有微生物能夠降解多種類的復(fù)雜石油烴，所以微生物群落結(jié)構(gòu)的穩(wěn)定性和多樣性是影響石油烴生物降解性能的重要因素，協(xié)同微生物群落呈現(xiàn)出更多的優(yōu)勢，如更多分解代謝基因和共代謝作用，利于石油烴污染物的自然衰減作用，可以降低石油烴的健康風(fēng)險.

圖 4 微生物物種關(guān)系Circos圖Fig.4 Microorganism species relationship Circos diagram

3 結(jié)論

a) A3監(jiān)測井石油烴濃度(0.06 mg/L)最低，A5監(jiān)測井石油烴濃度(6.61 mg/L)最高，A1～A4監(jiān)測井石油烴濃度差別不顯著，A5監(jiān)測井石油烴濃度約為其他監(jiān)測井的60～100倍.

b)利用HERA軟件以石油烴為目標(biāo)污染物進(jìn)行健康風(fēng)險評估，結(jié)果表明，該場地C～C、C～C的篩選值均為0.190 mg/L，C～C、C～C的篩選值均為0.142 mg/L. A1～A4監(jiān)測井的健康風(fēng)險在可接受水平. A5監(jiān)測井總石油烴危害商為44.990，是可接受值的45倍，表示A5監(jiān)測井的風(fēng)險較高，對場地周邊人群的健康損害風(fēng)險較大. 在同一場地條件下，石油烴濃度越高，其健康風(fēng)險越高. 建議根據(jù)不同監(jiān)測井的健康風(fēng)險選用合適的修復(fù)技術(shù)對場地進(jìn)行經(jīng)濟(jì)有效的修復(fù). 例如風(fēng)險較高的監(jiān)測井采取抽出處理技術(shù)、化學(xué)氧化技術(shù)、可滲透反應(yīng)墻技術(shù)等，風(fēng)險在可接受范圍內(nèi)的監(jiān)測井采取監(jiān)測自然衰減技術(shù).

c) 4個監(jiān)測井的微生物群落組成略有差異，A1、A2、A4、A5監(jiān)測井均含有變形菌門和厚壁菌門，且在A5監(jiān)測井中絕對豐度較高，說明一定范圍內(nèi)較高的石油烴污染會刺激有石油烴降解功能的微生物生長，如變形菌門、厚壁菌門和擬桿菌門，對污染物有更好的抵抗能力和降解能力，從而正向提高對石油烴的生物降解能力，降低石油烴健康風(fēng)險. 不同種類石油烴降解菌對石油烴的降解機(jī)制是今后的研究重點(diǎn).