基于Visual MODFLOW的煤礦開采對地下水影響的數值模擬研究

趙曉玲，王　博

（安徽省地質環境監測總站，安徽合肥230000）

基于Visual MODFLOW的煤礦開采對地下水影響的數值模擬研究

趙曉玲*，王博

（安徽省地質環境監測總站，安徽合肥230000）

為了研究該煤礦開采對研究區內地下水的影響，通過對研究區的水文地質條件的分析，構建了水文地質概念模型，并采用Visual MODFLOW軟件就開采3煤層對地下水的影響進行了數值模擬研究。結果表明：3煤層開采后，潛水含水層的地下水位變化幅度較小，最大水位降深為6m；基巖裂隙含水層的地下水位降幅較大，最大水位降深為25m。

煤礦開采；含水層；地下水位；數值模擬

隨著我國經濟和社會的高速發展，對能源的需求也日益增長，而在我國的能源結構中，煤炭資源占比超過7成，煤炭資源的開發為經濟和社會的發展提供了能源和動力。但是在煤炭資源開發的過程中，也產生了采空區塌陷，地下水漏斗、水質污染和生態環境惡化等不良影響和危害。某煤礦所處地區降雨量小、蒸發量大，地下水資源匱乏，區內工農業生產用水、人民生活用水和生態用水大部分來源于地下水。煤礦開采與工農業、生活以及生態用水之間存在一定矛盾，亟待解決。因此，為了研究該煤礦開采對地下水的影響［1-2］，首先對該區的水文地質條件進行了分析，并構建了該區的水文地質概念模型，然后采用Visual MODFLOW軟件就該煤礦開采3煤層對其上覆含水層地下水系統的影響進行了數值模擬研究，對煤礦的合理開采和所處地區地下水資源保護具有一定的指導意義。

1　水文地質概況

圖1　研究區位置示意圖

本次研究的煤礦位于安徽省淮北市城東北的杜集區石臺鎮境內（圖1），地處淮北平原北部，地形平坦，地表高程在35～37.6m之間；地貌類型沖洪積平原地貌單元；研究區屬北方型大陸氣候與濕潤性氣候之間的季風氣候區，氣候溫和，日照充足，春秋季明顯短于冬夏季，春夏季節盛行東南風，炎熱潮濕多雨，秋冬季盛行東北風，寒冷干燥少雨。

研究區的地下水為類型主要為第四系松散巖類孔隙水和基巖裂隙水。松散巖類孔隙水由細砂、粉砂層組成、存在砂質粘土夾層，總厚33～38m，平均為35m，較穩定。地表以下4～5m為潛水。砂層含水組鉆孔單位涌水量為q=0.24mL/（s·m），滲透系數k=1～2m/d，水質類型為HCO3--CI--Ca2+-Mg2+型。基巖裂隙水主要為下石盒子組中、細砂巖孔隙裂隙含水層，厚度15～45m，一般在25m左右，該含水組砂巖裂隙不發育，富水性較弱。鉆孔單位涌水量q=0.000085～0.0624mL/（s·m），滲透系數k=0.00047～0.202m/d，水質類型屬于CI-—Na+或CO23--C1--Na+型。

2　地下水模型的建立

2.1水文地質概念模型

根據研究區水文氣象、水文地質特征、主要含水層的流場特征以地層結構的分析結果，并依照主要含水層的滲透性能、地下水的遷移特征及水力聯系等，對實際的水文地質條件進行概化［3-5］，建立水文地質概念模型（圖2）。模擬區地下水流方向自北東向南西，因此將東部與北部設置為補給邊界，將南部和西部設置為排泄邊界，主要含水層的邊界劃為二類邊界。潛水的自由面作為上部邊界，底部邊界則為二疊系下二疊統山西組泥巖。

圖2　水文地質概念模型示意圖

2.2數學模型

根據水文地質模型，將實際地下水流概化為非均質各向同性三維非穩定地下水系統，可用如下偏微分方程來描述［6-8］：

式中：Kx、Ky、Kz——x、y、z方向的滲透系數；

h——含水層水頭；

ε——源匯項；

μ——潛水含水層潛水面上的重力給水度；

h0——含水層的初始水頭；

Γ0——潛水面邊界；

Γ1——滲流區第二類邊界；

Kn——邊界法向方向的滲透系數；

Ω——計算區范圍；

p——潛水面上的降水入滲和蒸發；

2.3數值模型結構及其識別與檢驗

2.3.1數值模型結構

根據研究區水文地質概念模型和數學模型，采用Visual MODFLOW 4.6構建數值模型并進行求解。研究區長約12.5km，寬約12.2km，總面積約152.2km2。計算前先對研究區進行網格剖分，網格間距為100m，在X、Y方向上分別剖分136×133，總網格共18088個。并根據模擬范圍將其范圍外地網格設定為不活動單元格。根據地下水的長期觀測資料，選取2014年3月～2015年2月的水位資料，應力期定為1個月，時間步長根據模型自動控制。初始地下水位數據采用2014年3月的統一觀測資料，利用Kriking插值法獲得初始水位等值線。

2.3.2模型的識別與檢驗

采用試估校正法對模型進行識別與檢驗，將滲透系數、貯水率、入滲補給系數和邊界條件等［9-10］輸入軟件中，根據模擬計算的水位值與實際觀測值進行對比，根據擬合情況不斷修正各參數，直到結果符合誤差的要求。根據研究區內的水位觀測的位置，本次選取J1113、J1307、J913、觀36等4個觀測點，在2014年3月～2015 年2月的地下水位觀測值與模擬值進行對比，分別計算各觀測點的模擬水位與實測水位之間的誤差，均滿足要求，擬合結果比較符合實際。各觀測孔的實測地下水位與模擬地下水位的擬合曲線如圖3所示。

2.3.3識別檢驗后的水文地質參數

通過對研究區實測水位與模擬水位的擬合、對地下水均衡的分析、對邊界流量和水文地質參數的調節，得出較為符合實際情況的水文地質參數，如表1所示。

3　煤層開采對地下水的影響

圖3　觀測孔實測與模擬地下水位的擬合曲線

表1　研究區識別檢驗后的水文地質參數

煤層開采過程中，對地下水采取排水疏干措施，會導致地下水流場發生明顯的改變，部分區域地下水甚至會被疏干，會導致各含水層出現水頭差并進一步加大。若各含水層之間無隔水層，則對淺層地下水的水位影響顯著，對煤礦的開采也有著巨大的影響；若各含水層之間存在弱透水地層，則兩者之間的滲透系數決定了越流量的大小，進而決定了淺層含水層水位；若各含水層之間存在連續完整且滲透系數很小的隔水層，煤層開采對淺層含水層影響很小。

研究區的潛水含水層與基巖裂隙含水層之間存在相對隔水層，該隔水層厚度不均勻，有多處該層缺失，因此，潛水含水層與基巖裂隙含水層中的地下水可以互相補給。基巖裂隙含水層是研究區內3煤層的直接充水水源，該煤層開采時對地下水的疏干作用，會導致基巖裂隙含水層中的地下水位大幅降低，進而與上覆潛水含水層之間形成較大水頭差，潛水含水層會對基巖裂隙含水層進行補給，影響潛水含水層中的地下水位。

本次根據地下水數值模型對3煤層開采后的潛水含水層和基巖裂隙含水層的地下水水位的變化情況進行的模擬，并進行了分析預測。3開采后各含水層的水位等值線和水位降深如圖4所示。

由圖4a和圖4b可以看出：3煤層開采后潛水含水層的水位等值線有所降低，但是變化幅度較小，最大水位降深為6m，對潛水含水層的影響不大；在開采疏干作用下，局部地段潛水含水層水位已近徹底疏干，快達到潛水含水層底板，而多部分地段的地下水位下降幅度很小，造成差別的原因主要在于含水層的滲透性能和給水能力有較大差異，在透水性和給水能力較強的地段，疏干水量大于補給量，地下水位下降，在透水性和給水能力較弱的地段，地下水位下降幅度不明顯。

圖4　3開采后各含水層的水位等值線和水位降深圖

由圖4c和圖4d可以看出：3煤層開采后基巖裂隙含水層的水位等值線降幅較大，最大水位降深為25m；造成上述差別的原因是潛水含水層與基巖裂隙含水層之間的相對隔水層有缺失，基巖裂隙含水層是3煤層的直接充水水源，煤層開采導致基巖裂隙含水層中的地下水位大幅降低，潛水含水層對其進行補給，進而導致局部地段潛水含水層的地下水位降幅較大。

4　結論

（1）采用Visual MODFLOW軟件對研究區煤層開采對含水層地下水位進行數值模擬分析，其模擬結果與實際觀測結果較為接近，基本反映了研究區的分布情況。

（2）采用試估校正法對模型進行識別與檢驗，通過對研究區實測水位與模擬水位的擬合、對地下水均衡的分析、對邊界流量和水文地質參數的調節，根據擬合情況不斷修正各參數，得出較為符合實際情況的水文地質參數。對研究區內的水位觀測點1個水文年內的地下水位觀測值與模擬值進行對比，擬合結果比較符合實際。

（3）3煤層開采后，在開采疏干作用下，局部地段潛水含水層水位已近徹底疏干，快達到潛水含水層底板，潛水含水層的水位等值線有所降低，但是變化幅度較小，最大水位降深為6m，對潛水含水層的影響不大；而3煤層開采后基巖裂隙含水層的水位等值線降幅較大，最大水位降深為25m。

［1］付佳妮，孫建明，林青.青島市黃島區潮河水源地地下水數值模擬［J］.山東國土資源，2015，31（7）：47-49.

［2］張文喆，王錦國，徐爍.徐州奎河淺層地下水數值模擬［J］.科學技術與工程，2013，13（13）：3793-3797.

［3］邱康利，李旭，賀江輝，等.Visual MODFLOW在地下水數值模擬中的應用［J］.煤炭與工化工，2014，37（6）：54-56.

［4］李志邦，張永波.Visual MODFLOW在煤礦開采地下水數值模擬中的應用［J］.礦業安全與環保，2014，41（4）：63-65.

［5］寇文杰.格爾木河流域地下水數值模擬［J］.水文地質工程地質，2013，40（1）：34-40.

［6］陳社明，王威，劉宏偉，等.河北省某擴建煤礦地下水流場變化數值模擬［J］.人民黃河，2015，37（10）：70-73.

［7］高翠萍，杜新強，郭中小，等.樺南縣地下水數值模擬研究［J］.人民黃河，2014，36（1）：78-81.

［8］孫琦，高為超，陳劍杰，等.西北某干旱區多級儲水洼地地下水數值模擬［J］.地下水，2012，34（1）：38-41.

［9］范堯，趙正國，黃旭，等.基于Visual MODFLOW的騰龍臺地地下水數值模擬研究［J］.水利與建筑工程學報，2015，13（4）：195-200.

［10］涂亮，宋漢周.基于Visual MODFLOW的黑河中游地下水流數值模擬［J］.勘察科學技術，2009（2）：19-23，33.

TD82

A

1004-5716（2016）07-0124-05

2016-04-02

2016-04-06

趙曉玲（1983-），女（漢族），河南禹州人，工程師，現從事水工環地質工作。