陸 嘉，張春暉，何緒文，梁 鵬，李時蓓，劉俐媛，胡喚雨

（1.中國礦業(yè)大學(xué) （北京），北京 100083；2.環(huán)境保護部環(huán)境工程評估中心，北京 100012；3.中國科學(xué)院高能物理研究所，北京 100049）

隨著我國一系列鼓勵煤層氣開發(fā)政策、規(guī)劃的實施，區(qū)域煤層氣規(guī)?；蜕虡I(yè)化開發(fā)將逐步展開［1］。煤層氣是一種潔凈的高熱值能源資源，對其地面開發(fā)和綜合利用不僅可以減少傳統(tǒng)資源的用量，進而減少溫室氣體排放，還有利于改善煤礦安全生產(chǎn)條件和周圍生態(tài)環(huán)境［2］，但同時煤層氣開發(fā)中的環(huán)境問題也逐步顯現(xiàn)。煤層氣在采氣過程中，會有部分水氣被帶入地表，稱作采出水。由于大量的煤層水與煤層氣共存，煤層氣的開采必須通過排水、降壓、解吸實現(xiàn)。煤層排水水質(zhì)和水量因各地區(qū)水文地質(zhì)條件不同而不同，通常，煤層氣采出水具有較高的堿性［3］，含有較多的鈉、鋇、重碳酸巖、鐵，并有較強的導(dǎo)電性［4］，且其中富含雜質(zhì)氣體和有毒、有害物質(zhì)，這些物質(zhì)在煤層和煤層氣開采過程中，將會污染大氣和地下水［5］。采出水通常會被送至采氣井附近的帶有底部防滲措施的臨時集水池，待積滿后由專用車輛運輸至污水處理廠集中處理。但在實際調(diào)查中發(fā)現(xiàn)，部分井場未建防滲收集池，僅用泥沙圍建臨時集水池收集采出水，這樣會造成采出水回滲入地下，如果采出水中含有污染物則可能對地下水環(huán)境造成污染。

可溶性污染物是地下水污染的主要原因。根據(jù)滲流理論，考慮地下水中水頭變化的影響，利用有限元方法，對可溶性物質(zhì)在地下水中的運移進行數(shù)值模擬，得到了地下水壓力水頭和可溶性污染物在地下水中的遷移變化規(guī)律，為可溶性污染物對地下水的污染研究和預(yù)測提供重要的分析數(shù)據(jù)［6］。目前常用的地下水?dāng)?shù)值模擬軟件主要有Visual MODFLOW［7－8］、GMS［9］以 及 Feflow［10］，其 中Visual MODFLOW和GMS使用的數(shù)值求解方式為有限差分法，F(xiàn)eflow使用的數(shù)值求解方式為有限元法［11］。因有限差分法易于和物理模型直接聯(lián)系，運算比較方便，因此在滲流問題中，有限差分法用的更為廣泛，因此選用有限差分法來對采出水對地下水影響進行研究［12］。

1 研究區(qū)概況

1.1 地層特征

本研究區(qū)域地處淺山丘陵溝壑區(qū)，地貌形態(tài)分為河谷區(qū)和梁狀黃土丘陵區(qū)兩大類。區(qū)內(nèi)出露地層由老到新為：太古界、元古界、寒武系、奧陶系、石炭系、二疊系、三疊系和第四系。區(qū)域地形總的趨勢是西北高東南低，最高點海拔1432m；最低點海拔400m。一般標(biāo)高在600～800m之間，相對高差1032m。

1.2 含水層特征

河谷區(qū)地勢平坦，分布有第四系松散層的孔隙含水層，溝谷兩側(cè)出露二疊系石盒子地層，基巖裸露，地勢高陡，在上石盒子組下部有一層砂巖含水層，為下石盆子組下段砂巖含水層，其富水性較其他含水層相對較強，水質(zhì)良好，河谷區(qū)居民多以自然村為單位就近建設(shè)供水設(shè)施開采該含水層中的地下水作為生活用水。

地下水按其巖性及儲水空間可分為：第四系松散巖類孔隙水、石炭二疊系砂（灰）巖裂隙水及奧灰?guī)r巖溶裂隙水三種類型?？傮w而言，第四系松散巖類含水不豐富，石炭二疊系富水性、透水性不強，奧灰?guī)r巖溶水富水性、透水性強，但極不均一。根據(jù)收集資料，區(qū)內(nèi)地下水含水層劃分見表1。

表1 研究區(qū)域地下水含水層劃分表

2 模型概化

本研究利用數(shù)值模擬的方法對污染物進入地下水后的污染范圍及遷移路徑進行預(yù)測。模擬區(qū)域為陜西某煤層氣開采井田，模擬使用Visual MODFLOW軟件。通過使用MODFLOW模塊構(gòu)建地下水流動模型，使用MT3D模塊構(gòu)建污染物運移模型，對區(qū)域地下水流場及污染羽進行模擬。

該水文地質(zhì)單元地處某褶斷帶，模擬區(qū)域北邊界以F22斷層為邊界，為隔水邊界。東邊界，南邊界分別依托F1、F14斷層設(shè)置，為隔水邊界。據(jù)調(diào)查資料顯示，區(qū)域北部的河流水位常年與第四系潛水含水層的水位基本一致，因此可判定為一類邊界，南側(cè)的支流為季節(jié)性河流，概化為河流邊界，模型西邊界為二類邊界，模型上邊界接受降雨補給。區(qū)域地質(zhì)概況如圖1所示。

從空間上看，評價區(qū)地下水流整體上以水平運動為主，地下水系統(tǒng)符合質(zhì)量守恒定律和能量守恒定律；在常溫常壓下地下水運動符合達(dá)西定律；地下水系統(tǒng)的輸入輸出隨時間、空間變化，但地下水位動態(tài)變化受降水影響明顯，因此將地下水視為非穩(wěn)定流；參數(shù)隨空間變化，體現(xiàn)了系統(tǒng)的非均質(zhì)性；在水平方向上，參數(shù)沒明顯的方向性，為各向同性。

綜上所述，研究區(qū)可概化成非均質(zhì)各向同性、空間三維結(jié)構(gòu)、非穩(wěn)定地下水流系統(tǒng)，即地下水系統(tǒng)的概念模型。模型概化為四層，其中第一層概化為第四系松散層的孔隙潛水含水層；第二層概化為第四系松散層含水層以下，二疊系砂巖裂隙承壓含水層以上的三疊系泥巖相對隔水層；第三和第四層分別對應(yīng)上石盒子組上段砂巖裂隙含水層和上石盒子組下段砂巖裂隙含水層。由于附近村落的水源井的取水點位于上石盒子組下段砂巖含水層，因此我們主要關(guān)注污染物到達(dá)該含水層后的運動動態(tài)。

圖1 區(qū)域地質(zhì)概況圖

3 模擬結(jié)果與討論

通過輸入?yún)?shù)搭建地下水流模型，為了確保模型準(zhǔn)確無誤，需要通過模型識別來校驗參數(shù)的準(zhǔn)確性。在模型識別的過程中，以模型構(gòu)建當(dāng)日的600天前作為初始時間，輸入這600天的歷史觀測值，運行模型，可以得到模型計算值，繼而比較觀測值與計算值的誤差。以第四含水層為例，以模型構(gòu)建當(dāng)天作為比對時間點，模型識別前與識別后的水位擬合情況分別如圖2和圖3所示。

圖2 第四含水層水位模擬誤差分析圖（模型識別前）

通過圖2和圖3可以看出，在模型識別后各觀測點的觀測水位值與計算水位值擬合情況較好，通過比對五個觀測井的最大水位誤差僅為0.8m，可認(rèn)為該水流模型基本正確，繼而在水流模型的基礎(chǔ)上建立污染物運移模型，對污染情況進行預(yù)報。通過對該井田采出水水質(zhì)進行采樣分析，得到的結(jié)果如表2所示。

圖3 第四含水層水位模擬誤差分析圖（模型識別后）

表2 陜西省某氣田采出水水質(zhì)報告

從表2可以看出，氯化物的濃度嚴(yán)重超標(biāo)，因此選擇氯化物作為預(yù)測因子。在構(gòu)建污染物運移模型時，模擬起始日期為模型構(gòu)建當(dāng)日，模擬時間為3000天，運行模型，得出未來3000天內(nèi)4個時段的污染羽圖（圖4，圖4中五個柱狀體為水源井）。

圖4 不同預(yù)測時間點的污染羽圖

研究結(jié)果表明，3000天內(nèi)采出水持續(xù)滲入到地下水中，基本不會對周圍村民取水用水造成影響。分析原因有如下兩個：第一，污染源與水源井的距離相對較遠(yuǎn)，污染物到達(dá)水源井需要很長的時間；第二，第二含水層至第四含水層的滲透系數(shù)均很小，延緩了污染物的滲透速度，有效地減緩了污染物在含水層中的污染進程。由此可見，按最不利情況考慮，在連續(xù)采氣3000天內(nèi)，仍然沒有修建防滲集水池收集采出水的情況下，含高濃度氯化物的廢水持續(xù)滲入到地下水中，基本不會對周圍村民取水用水造成影響。

4 結(jié)論

本文以氯化物為預(yù)測因子，利用地下水?dāng)?shù)值模擬的方法，使用地下數(shù)值模擬軟件Visual MODFLOW搭建研究區(qū)域的水文地質(zhì)模型，對煤層氣采出水滲入地下水后其中所含污染物的遷移路徑和污染范圍進行模擬。結(jié)果表明，因該井場距離周邊水源井距離較遠(yuǎn)，且含水層的滲透性較差，減緩了污染物的下滲速度，污染物基本不會對周圍居民飲水健康造成影響。

［1］胡連伍，陳海霞，楊衛(wèi)國.地面抽采煤層氣環(huán)境影響和管理對策探討［J］.能源環(huán)境保護，2009，23（3）：38-42.

［2］楊海亮.煤層氣地面開發(fā)對環(huán)境的影響及防護對策［J］.能源環(huán)境保護，2009，23（3）：26-30.

［3］劉妍萍，王斌，劉新春等.UF-RO處理煤層氣田產(chǎn)出水的工程應(yīng)用［J］.環(huán)境工程，2011，29（6）：1-4.

［4］張衛(wèi)東，郭敏，閆愛華.美國煤層采出水處理技術(shù)現(xiàn)狀［J］.中國煤層氣，2010，7（3）：11-16.

［5］嚴(yán)群，王慶偉，李金海.煤層氣開采對環(huán)境的影響［J］.中國煤炭地質(zhì)，2008，20（12）：57-59.

［6］劉勇.可溶性污染物在地下水中運移的數(shù)值模擬［J］.資源環(huán)境與工程，2008，22（3）：359-361.

［7］R.K.Tiwary，R.Dhakate，V.AnandaRao，et al.Assessment and prediction of contaminant migration in ground water from chromite waste dump［J］.Environ Geol，2005（48）：420-429.

［8］Tian Gan，Zhao Chunhu.Applications of MODFLOW in quantitative analysis coal seam floor water inrush condition［C］. 2011IEEE 3rd International Conference on Communication Software and Networks （ICCSN 2011），2011：711-714.

［9］戴永旺，常云龍.基于GMS軟件的地下水?dāng)?shù)值模擬應(yīng)用研究［J］.商品與質(zhì)量：學(xué)術(shù)觀察，2012（11）：75，84.

［10］李雪艷，梁立章，聶大鵬.基于FEFLOW傍河水源地地下水?dāng)?shù)值模擬［J］.東北水利水電，2012（8）：22-25.

［11］童彥釗，吳平，韓強強.地下水?dāng)?shù)值模擬軟件簡述［J］.科技信息，2011（29）：48，27.

［12］張宏仁，李俊亭.有限差分法與有限單元法在滲流問題中的對比［J］.水文地質(zhì)工程地質(zhì)，1979（2）：50-54.