原杰輝

（煤炭工業(yè)太原設(shè)計研究院， 山西 太原 030001）

經(jīng)驗交流

MODFLOW水流數(shù)值模型在礦山環(huán)境影響評價中的應用

原杰輝

（煤炭工業(yè)太原設(shè)計研究院， 山西 太原 030001）

以某鐵礦為例，通過對MODFLOW在礦山開發(fā)活動中的應用進行研究，對研究區(qū)礦山開采疏排礦坑水過程進行了水流數(shù)值模擬，預測礦山開采活動對周邊地下水環(huán)境的影響，并結(jié)合礦區(qū)水文地質(zhì)條件提出礦區(qū)開發(fā)過程中針對水環(huán)境保護可采取的保護措施，旨在對類似礦山開發(fā)對水環(huán)境的影響分析工作起到一定參考作用。

MODFLOW 數(shù)值模擬 環(huán)境保護 水資源

MODFLOW是現(xiàn)今世界范圍內(nèi)應用最廣泛的地下水流模擬程序，國內(nèi)引進應用也很多，在礦山防治水方面具有廣泛的應用[1-2]。MODFLOW模擬地下水的方法，是用多層的長方形網(wǎng)格刻畫三維含水層系統(tǒng)，輸入含水層參數(shù)，然后對每個單元格建立非穩(wěn)定流的有限差分方程進行數(shù)值求解[3]。Visual MODFLOW作為其前后數(shù)據(jù)處理的軟件，可以在Windows系統(tǒng)中運行，是目前較主流的MODFLOW運行載體。

研究區(qū)地處呂梁山脈東麓的狐偃山一帶，屬中低山區(qū)。區(qū)內(nèi)溝谷切割，地貌類型復雜，地表大部分被第四系黃土覆蓋，植被發(fā)育、荊棘叢生。礦區(qū)地勢總體南高北低，海拔高度在（+1694～+1918）m之間，研究區(qū)內(nèi)相對高差達到224m。區(qū)內(nèi)溫差大、蒸發(fā)強烈、日照時間長，礦區(qū)及礦床水文地質(zhì)條件較為復雜。

1 研究區(qū)地質(zhì)與水文地質(zhì)條件簡介

研究區(qū)處于狐偃山勘探區(qū)東南分水嶺北側(cè)，區(qū)內(nèi)植被覆蓋茂盛，溝谷中基巖裸露，地形高差相對較大。區(qū)域出露地層有：奧陶系、石炭系、二疊系、三疊系和新生界。地質(zhì)構(gòu)造方面位于上白泉橢圓形不規(guī)則環(huán)形構(gòu)造的南部，南岔背斜（軸向NNW－SSE）南段的東翼，礦區(qū)內(nèi)巖層傾斜130°～250°，傾角0°～10°，沿層面有波狀起伏現(xiàn)象。本礦區(qū)以北地表零星出露有二長斑巖，自上而下可見二長斑巖脈順層侵入于太原組內(nèi)、斑狀二長巖侵入于奧陶系中統(tǒng)峰峰組地層中，此外奧陶系中統(tǒng)上馬家溝組二段下尚見有霓輝二長斑巖侵入。

礦區(qū)內(nèi)共圈出隱伏磁鐵礦體5個：1、2號礦體賦存于奧陶系中統(tǒng)峰峰組與（巨）斑狀二長巖的接觸帶上，屬熱液－接觸交代矽卡巖型磁鐵礦；3、4號礦體賦存于本溪組底部的山西式鐵礦或鐵鋁巖巖體層位；28號礦體賦存于上馬家溝三段的大理巖構(gòu)造裂隙中，屬充填交代型磁鐵礦。

根據(jù)地下水的賦存條件和水力特征，區(qū)內(nèi)地下水含水層主要為：松散巖類孔隙水、碎屑巖夾碳酸鹽巖類巖溶裂隙水和碳酸鹽巖巖溶水。地下水的補給主要為大氣降水沿各類巖層的風化帶下滲。礦區(qū)范圍內(nèi)全部被第四系黃土覆蓋，基巖含水層隱伏發(fā)育于第四系黃土層之下。礦床產(chǎn)于分水嶺地段，礦床周邊沖溝發(fā)育，地形有利于自然排泄。據(jù)上白泉礦區(qū)普查地質(zhì)報告資料，當?shù)貛r溶水位標高為1 410m，部分礦體位于巖溶水位以下，具有較高的靜水壓力。

2 地下水流數(shù)值模型的建立

2.1 建立水文地質(zhì)概念模型

鐵礦開采對礦體涉及的石炭二疊系碎屑巖夾碳酸鹽巖類巖溶裂隙水含水層、奧陶系碳酸鹽巖巖溶水含水層、侵入巖風化裂隙水含水層造成疏排影響。根據(jù)地質(zhì)資料，礦體以上并沒有發(fā)育連續(xù)穩(wěn)定的含水層，多為淺部基巖風化裂隙發(fā)育后形成的分散含水層。評價將其假設(shè)為具有統(tǒng)一水力聯(lián)系的淺部基巖風化裂隙水含水層，作為評價的目標含水層[4]。

根據(jù)研究區(qū)地下水流向整體上自南向北的現(xiàn)狀，取目標含水層地下水流場南側(cè)上游地表分水嶺作為南邊界；結(jié)合地形取流場下游1.7～2.1 km作為北邊界。這兩個邊界均處理為二類流量邊界。模擬區(qū)東、西邊界均取地形高點連線圈定，處理為流量為零通量的隔水邊界。

評價區(qū)目標含水層地下水流從空間上假設(shè)以水平運動為主，地下水系統(tǒng)符合質(zhì)量守恒定律和能量守恒定律：未開采階段，目標含水層在水平方向上接受南側(cè)地下水徑流補給，在垂直方向上接受降水入滲補給，模型內(nèi)部整體上向北方向徑流，在模型北邊界處向下游方向排泄。開采階段，模型于采區(qū)抽水形成降落漏斗，目標含水層水平方向接受上游及周邊地下水徑流補給，在垂直方向上接受降水入滲補給，模型內(nèi)部局部向降落漏斗中心方向徑流，以礦坑水的形式排泄。由于目標含水層地下水位埋藏整體相對較深，不考慮蒸發(fā)排泄。

2.2 建立MODFLOW水流數(shù)值模型

將1：50 000的鐵礦地形圖導入模型作為計算模擬區(qū)的剖分底圖，模型識別區(qū)域南北長約6.8 km，東西長約7.3 km。將模型區(qū)在水平方向上剖分為100×100的微小單元格，在垂直方向上按水文地質(zhì)概念模型設(shè)計層次剖分。

選擇初始水位校核為模型識別與參數(shù)校正階段，將水文地質(zhì)調(diào)查和抽水試驗得到的水文地質(zhì)參數(shù)、邊界條件、水頭初始條件作為模型調(diào)參的初始值，運行計算模型，將實測水位和計算水位進行擬合分析，如果計算水位與實測水位相差很大，則根據(jù)參數(shù)變化范圍和實際水位差值，重新給定一組參數(shù)，再迭代計算，直至二者擬合較好為止。

通過調(diào)參計算，實際水位和計算水位等值線的水位擬合小于0.5m的絕對誤差占已知水位的85%以上，擬合結(jié)果較好，表明所建的模型、對水文地質(zhì)條件的概化、邊界條件的確定都與研究區(qū)域?qū)嶋H情況吻合較好，因此，校正后的模型可以滿足預測評價的要求。

3 地下水流數(shù)值模擬預測與分析

根據(jù)礦區(qū)開拓布置，開發(fā)利用方案根據(jù)礦體分布及賦存形態(tài)、地表地形地質(zhì)，確定采用主豎井開拓系統(tǒng)，主豎井開拓、斜風井回風，對礦區(qū)內(nèi)1、2、4號礦體進行開采，在模型中將礦山開采期間產(chǎn)生的礦坑涌水視為穩(wěn)定連續(xù)的抽水過程，模擬階段依開采順序在開采水平內(nèi)設(shè)置抽水孔，將正常涌水量以抽水井的形式排泄出去。

模擬預測結(jié)果顯示，礦山開采后，形成以擬采礦體中底板較低點（1、2礦體底板）為中心的地下水位降水漏斗，漏斗中心位于1、2號礦體底板低點附近，最大水位降深約為42.72m，往周邊方向水位降深逐漸變小。地下水位受影響面積約為1.79 km2。

根據(jù)預測結(jié)果，分析本研究區(qū)內(nèi)礦山開采對地下水的影響具有如下特點。

1）礦山開采過程中，地下水降落漏斗隨時間由小變大，直至一定時間段后補排平衡達到穩(wěn)定；同一時間節(jié)點內(nèi)地下水降深隨著與礦體距離增大而減小，結(jié)合礦體等高線可以看出礦體周邊地下水位降至最低開采高度之下。

2）隨著開采礦體的轉(zhuǎn)移，抽水位置發(fā)生變化，地下水降落幅度與影響范圍隨之變化。

3）鐵礦開采對地下水位的影響主要是抽排礦坑水造成的地下水位下降，其下降中心為開采礦體位置，在水平方向上影響礦體周邊地下水位，本礦開采時段內(nèi)礦體及周邊地下水位整體出現(xiàn)下降。

4）在礦體開采結(jié)束后，受影響區(qū)域內(nèi)地下水位在降落漏斗中心轉(zhuǎn)移后緩慢上升。由于區(qū)域降水量較小、水資源不豐富，且基巖風化裂隙水含水層滲透性較差、富水性很弱，水資源更新補給緩慢，礦區(qū)范圍地下水位恢復速度緩慢。

4 提出水環(huán)境影響保護措施

4.1 開采期間涵養(yǎng)水土，及時進行生態(tài)恢復

開采礦體使礦區(qū)及周邊地區(qū)地下水位下降，加劇水土流失，建議開采時一方面要嚴格實施分區(qū)開采，另一方面及時進行水土保持工作，涵養(yǎng)水土，降低鐵礦開采對淺部地下水資源和生態(tài)環(huán)境的影響。

4.2 建立地下水觀測網(wǎng)系統(tǒng)

結(jié)合觀測區(qū)地質(zhì)、水文地質(zhì)條件，建立地下水動態(tài)觀測網(wǎng)，以掌握地下水位動態(tài)變化規(guī)律，有效預測疏干涌水量，指導疏干工作。若實際開發(fā)中造成區(qū)域地下水水位嚴重下降，建設(shè)單位應及時組織水文地質(zhì)專家查找原因，針對性地制定工程防止措施和配套補救措施，對可能造成的不良影響給以經(jīng)濟補償，并根據(jù)項目可能誘發(fā)的環(huán)境水文地質(zhì)問題制定相應的監(jiān)測方案。

4.3 對水資源綜合利用

本礦開采對礦體圍巖含水層的破壞不可避免，該部分水資源主要以礦坑水的方式產(chǎn)生，排水均進入礦坑水處理站經(jīng)混凝、沉淀、過濾處理工藝處理后盡量尋求回用途徑，減少礦坑水排放。

4.4 制訂供水預案

以預測結(jié)果為基礎(chǔ)，針對受影響范圍內(nèi)的飲用水井，制訂由礦方提供水源和供水的預案，結(jié)合地下水動態(tài)觀測系統(tǒng)，一旦發(fā)現(xiàn)居民用水受到影響，立即啟動供水預案，保障居民用水不受影響。

5 結(jié)語

文章通過對MODFLOW在礦山開發(fā)活動中的應用進行研究，對研究區(qū)礦山開采疏排礦坑水過程進行了水流數(shù)值模擬，預測礦山開采活動對周邊地下水環(huán)境的影響，據(jù)此提出環(huán)境影響防范措施，對于類似礦山開發(fā)對水環(huán)境的影響分析工作具有一定參考作用。

[1] 馮更辰，郝俊杰，譚俊，等.Visual Modflow模型在白澗鐵礦區(qū)礦井涌水量預測中的應用[J].中國巖溶，2011（3）：271-277.

[2] 王潤堂.石硐溝銀多金屬礦床水文地質(zhì)條件分析[J].甘肅冶金，2008（2）：24-26.

[3] 王慶永，賈忠華，劉曉峰，等.VisualMODFLOW及其在地下水模擬中的應用[J].水資源與水工程學報，2007（5）：90-92.

[4] 薛禹群，謝春紅.地下水數(shù)值模擬[M].北京：科學出版社，2007.

App lication of MODFLOW Flow Numerical M odel in M ine Environmental Im pact Assessment

Yuan Jiehui
（Taiyuan Design and Research Institute of Coal Industry,Taiyuan Shanxi 030001）

Taking an iron mine as example,based on the research of the application of MODFLOW in mine development activities,the flow numerical simulation of drainage process ofmine drainage is conducted to predict the impact of mining activities on surrounding groundwater environment.Combined with the hydrogeological conditions,this paper puts forward protectivemeasures for water environment in developing process of themine so as to provide reference for the analysis of the impactofmine developmenton water environment.

MODFLOW;numerical simulation;environmental protection;water resources

P641.7

A

1672-1152（2017）05-0127-03

10.16525/j.cnki.cn14-1167/tf.2017.05.48

2017-10-23

原杰輝（1985—），男，碩士研究生，畢業(yè)于吉林大學，工程師，主要從事環(huán)境影響評價和水資源論證方面的工作。

（編輯：王瑾）