韓忠 劉樂(lè)軍 袁星芳 王文瑾 曲洪飛

摘? 要：焦家斷裂帶金礦密集，相鄰礦山開(kāi)采相互干擾明顯，單一礦坑的涌水量研究工作已經(jīng)不能滿(mǎn)足礦山生產(chǎn)需求。通過(guò)建立多礦山聯(lián)合開(kāi)采的焦家控礦構(gòu)造帶地下水流數(shù)值模型，模擬預(yù)測(cè)干擾開(kāi)采條件下不同開(kāi)采中段的礦坑涌水量及隨時(shí)間的變化規(guī)律，為礦山開(kāi)采制定防治水方案提供依據(jù)。模型中采用barrier障礙邊界模擬控礦構(gòu)造阻水帶對(duì)地下水水平流的阻礙作用，在傳統(tǒng)水位及觀測(cè)孔擬合的基礎(chǔ)上，進(jìn)行礦坑涌水量擬合校正，進(jìn)一步提高模型模擬精度。

關(guān)鍵詞：焦家斷裂;數(shù)值模擬;礦坑涌水;礦山聯(lián)合開(kāi)采;GMS

中圖分類(lèi)號(hào)：P641.4+1 文獻(xiàn)標(biāo)志碼：A? ? ? ? ?文章編號(hào)：2095-2945（2019）11-0117-03

Abstract： The gold deposits in Jiaojia fault zone are dense， and the mining of adjacent mines interferes with each other obviously. The research work of water inrush from a single mine can no longer meet the needs of mine production. By establishing a numerical model of groundwater flow in Jiaojia ore-controlling structural zone of multi-mine joint mining， the water inflow of pits in different mining sections under disturbed mining conditions and its variation with time are simulated and predicted. It provides the basis for the formulation of water prevention and control scheme in mining. In the model， the barrier boundary is used to simulate the hindrance of the ore-controlling structural water blocking zone to the horizontal flow of groundwater. On the basis of the traditional water level and observation hole fitting， the mine water inflow is fitted and corrected， and the simulation accuracy of the model is further improved.

Keywords： Jiaojia fault; numerical simulation; mine water gushing; combined mining; GMS

1 研究背景

礦坑涌水是困擾礦山開(kāi)采的重要隱患[1]，礦坑涌水量預(yù)測(cè)是礦山開(kāi)采的先行步驟和必要環(huán)節(jié)，可為礦山開(kāi)采設(shè)計(jì)中制定防、治水方案提供重要依據(jù)[2]。

對(duì)于開(kāi)采條件較簡(jiǎn)單的單個(gè)礦區(qū)，涌水量預(yù)測(cè)多采用水文地質(zhì)比擬法[3]、大井法[4]、水均衡法[5]、涌水量曲線方程法[6]、基于裘布依穩(wěn)定井流理論的解析法以及數(shù)值法[7]。本文采用數(shù)值法進(jìn)行焦家控礦構(gòu)造帶礦坑水模擬研究。采用國(guó)際通用的可視化軟件GMS（Groundwater Modeling System）[8]，建立焦家控礦構(gòu)造帶地下水流數(shù)值模型，預(yù)測(cè)礦山聯(lián)合開(kāi)采條件下不同時(shí)段、不同開(kāi)采深度的礦坑涌水量，為后期礦山開(kāi)采提供技術(shù)指導(dǎo)。

2 研究區(qū)概況

焦家斷裂帶位于煙臺(tái)萊州市境內(nèi)，南鄰青島，東接威海。大地構(gòu)造位置位于沂沭斷裂帶北段東側(cè)，膠北隆起西緣，南接膠萊凹陷。焦家斷裂帶屬壓扭性斷裂帶，走向10°～30°，傾向北西，傾角29°～56°，規(guī)模較大，影響范圍較廣，不僅控制著金礦體的分布與埋藏，同時(shí)也控制著地下水的富集、運(yùn)移及水力特性。

斷裂帶中間部位分布連續(xù)的斷層泥和兩側(cè)的碎裂巖具有良好的隔水性。斷裂帶兩側(cè)受構(gòu)造影響，構(gòu)造裂隙發(fā)育，是地下水賦存和運(yùn)動(dòng)的有利場(chǎng)所，是較好的富水部位。下盤(pán)的基巖構(gòu)造裂隙含水層（帶）是礦坑充水的直接含水層，也是金礦體的主要賦存部位。頂部為分布比較連續(xù)的隔水帶，下部是以二長(zhǎng)花崗巖為主體的隔水底板，與其它含水層之間的水力聯(lián)系較弱。

3 水文地質(zhì)概念模型

3.1 模型范圍及邊界條件

根據(jù)斷裂帶及附近區(qū)域的水文地質(zhì)條件，確定模型范圍及邊界條件：北部邊界為海岸線，設(shè)為定水頭邊界。南部邊界為米山-金華山形成的天然地下水分水嶺，設(shè)定為流量邊界。東西部邊界分別為諸流河和朱橋河，設(shè)為定水頭邊界。

3.2 斷裂帶水力性質(zhì)概化

在模擬區(qū)中部地區(qū)北東-西南向斷裂發(fā)育，主斷裂中心發(fā)育有連續(xù)穩(wěn)定主裂面，兩側(cè)為碎裂巖，該部位導(dǎo)水性和富水性很差，阻水性強(qiáng)，對(duì)地下水流具有控制和阻礙作用。采用MODFLOW模塊下的barrier障礙邊界對(duì)斷層進(jìn)行處理。

3.3 地下水流動(dòng)系統(tǒng)概化

礦坑地下水流從空間上看是以垂向運(yùn)動(dòng)為主、水平運(yùn)動(dòng)為輔，地下水為三維地下水流。地下水系統(tǒng)的輸入、輸出隨時(shí)間、空間變化，地下水系統(tǒng)為非穩(wěn)定的分布參數(shù)系統(tǒng)。地下水系統(tǒng)參數(shù)、補(bǔ)排項(xiàng)隨空間變化，體現(xiàn)了系統(tǒng)的非均質(zhì)性，所以參數(shù)概化為非均質(zhì)各向同性。

綜上所述，將研究區(qū)地下水系統(tǒng)概化成非均質(zhì)各向同性、空間雙層結(jié)構(gòu)、非穩(wěn)定流三維地下水流系統(tǒng)。

4 地下水流數(shù)值模型

4.1 數(shù)學(xué)模型

對(duì)研究區(qū)的地下水系統(tǒng)進(jìn)行抽象概化，在建立水文地質(zhì)概念模型的基礎(chǔ)上，進(jìn)一步轉(zhuǎn)化為數(shù)學(xué)模型。根據(jù)Darcy定律、質(zhì)量守恒定律以及能量守恒定律，并且不考慮水密度變化的情況下，得到地下水流動(dòng)偏微分方程：

式中：Kxx，Kyy，Kzz-滲透系數(shù)在x，y，z方向上的分量，我們假定滲透系數(shù)的主軸方向與坐標(biāo)軸的方向一致，量綱為（LT-1）;h-含水層的水位標(biāo)高（L）;W-單位體積流量（T-1），用以代表流進(jìn)匯或來(lái)自源的水量;Ss-對(duì)非承壓含水層為給水度，對(duì)承壓含水層為貯水率（L-1）;t-時(shí)間（T）。

4.2 地下水流模擬

4.2.1 空間剖分

根據(jù)研究區(qū)域大小及計(jì)算精度要求，在平面上將研究區(qū)剖分為60行×60列，垂向上分為第四系含水層和構(gòu)造裂隙含水層兩層。其中第一層第四系含水層有效單元642個(gè)，第二層構(gòu)造基巖裂隙含水層有效單元1641個(gè)。

4.2.2 時(shí)間離散和初始條件的確定

模擬期確定為2015年1月1日～2016年12月30日，將每個(gè)自然月作為一個(gè)應(yīng)力期，共24個(gè)應(yīng)力期。在建立擬合期初始流場(chǎng)時(shí)采用各含水層水井水位標(biāo)高及各長(zhǎng)觀孔的初始水位標(biāo)高，然后值賦給模型參與迭代計(jì)算。

4.2.3 源匯項(xiàng)的處理

對(duì)源匯項(xiàng)的處理，主要包括兩大類(lèi)。一類(lèi)是以含水層面狀補(bǔ)給率的形式給出，如降水入滲補(bǔ)給;另一類(lèi)是以點(diǎn)井的形式給出，如礦坑排水、山前側(cè)向徑流、河流滲流補(bǔ)給等。這兩類(lèi)源匯項(xiàng)均分配在活動(dòng)單元格上參與計(jì)算。

4.3 模型的識(shí)別與驗(yàn)證

4.3.1 觀測(cè)孔擬合

通過(guò)反復(fù)調(diào)整參數(shù)和均衡量，識(shí)別水文地質(zhì)條件，確定模型結(jié)構(gòu)、參數(shù)和均衡要素。選取有代表性的觀測(cè)點(diǎn)進(jìn)行擬合，計(jì)算各擬合點(diǎn)計(jì)算水位與觀測(cè)水位的絕對(duì)誤差在0.2～0.4之間。

長(zhǎng)觀孔擬合曲線圖中圓點(diǎn)為預(yù)測(cè)值，倒三角為長(zhǎng)觀孔的實(shí)際值。由擬合結(jié)果可知，大部分?jǐn)M合點(diǎn)的平均絕對(duì)誤差小于0.5m，所建立的模型滿(mǎn)足模擬精度要求。

4.3.2 礦坑涌水量擬合

為進(jìn)一步提高模型可靠性，提高礦坑涌水量預(yù)測(cè)精度，本次通過(guò)以往大量的涌水量觀測(cè)數(shù)據(jù)，對(duì)礦區(qū)開(kāi)拓巷道水平的實(shí)際涌水量與預(yù)測(cè)涌水量進(jìn)行擬合對(duì)比（表1）。由擬合結(jié)果可知，礦坑涌水量計(jì)算值相對(duì)誤差小于1.7%，擬合效果較為理想。

表1 各礦區(qū)涌水量擬合對(duì)比表（單位m3/d）

綜上所述，所建立的地下水流數(shù)值模型符合研究區(qū)水文地質(zhì)條件，基本反映了地下水系統(tǒng)的動(dòng)態(tài)及涌水量變化特征，故可利用模型進(jìn)行地下水水位、水量的預(yù)測(cè)預(yù)報(bào)。

5 礦坑涌水量預(yù)測(cè)

5.1 涌水量預(yù)測(cè)模型

本次礦坑涌水量預(yù)測(cè)分為兩部分：一是不同開(kāi)采中段涌水量預(yù)測(cè);二是各個(gè)開(kāi)采中段涌水量隨時(shí)間的動(dòng)態(tài)變化趨勢(shì)預(yù)測(cè)。針對(duì)以上兩種需求，將源匯項(xiàng)中的礦坑涌水量用排水溝的形式表達(dá)，排水底板標(biāo)高設(shè)置為要預(yù)測(cè)涌水量的開(kāi)采水平，模擬礦井隨開(kāi)采深度的增加涌水量動(dòng)態(tài)變化趨勢(shì)。

5.2 不同開(kāi)采中段涌水量預(yù)測(cè)

在當(dāng)前礦井開(kāi)采深度的基礎(chǔ)上，考慮未來(lái)開(kāi)采深度不斷加大條件下，預(yù)測(cè)不同開(kāi)采中段礦坑涌水量（圖2）。預(yù)測(cè)成果表明，隨著開(kāi)采深度的增加，各礦區(qū)涌水量隨之增加，增幅有差異。開(kāi)采至-2000中段，新城礦區(qū)、焦家礦區(qū)、寺莊礦區(qū)及望兒山礦區(qū)涌水量分別達(dá)到4565m3/d、29034m3/d、21140m3/d、45803m3/d，其中望兒山礦區(qū)的涌水量最大，新城礦區(qū)的涌水量最小，增加幅度也最小。

圖2 各礦區(qū)不同開(kāi)采中段預(yù)測(cè)涌水量曲線

可見(jiàn)，隨著金礦的開(kāi)采，礦坑涌水量將成為制約礦山開(kāi)采的重要因素。礦山生產(chǎn)過(guò)程中宜根據(jù)不同開(kāi)采中段的涌水量制定相應(yīng)的防、治水措施，確保安全生產(chǎn)。

5.3 涌水量動(dòng)態(tài)變化預(yù)測(cè)

假定5年開(kāi)采一個(gè)水平中段，20年可以開(kāi)采4個(gè)水平中段，利用上述地下水流數(shù)值模型對(duì)焦家成礦帶開(kāi)展未來(lái)20年開(kāi)采這4個(gè)水平的變化趨勢(shì)進(jìn)行預(yù)測(cè)，圖3為焦家礦區(qū)最后一個(gè)中段涌水量動(dòng)態(tài)預(yù)測(cè)結(jié)果。

由涌水量預(yù)測(cè)動(dòng)態(tài)曲線圖可見(jiàn)，各礦區(qū)的涌水量季節(jié)性變化明顯，總體夏季呈上升趨勢(shì)，冬季呈下降趨勢(shì)。涌水量總體在某一固定值附近波動(dòng)，最后趨于水平，說(shuō)明隨著開(kāi)采的進(jìn)行涌水量最終趨于穩(wěn)定值。

6 結(jié)論

（1）聯(lián)合開(kāi)采條件下，焦家控礦構(gòu)造帶礦坑涌水量開(kāi)采至-2000中段，新城金礦、焦家金礦、寺莊金礦、望兒山金礦礦坑涌水量將分別達(dá)到4565m3/d、29034m3/d、21140m3/d、45803m3/d。

（2）受河流通過(guò)導(dǎo)水?dāng)嗔褜?duì)礦坑水的補(bǔ)給作用影響，寺莊礦坑涌水量隨季節(jié)性變化不明顯，涌水量較為穩(wěn)定。其它礦區(qū)礦坑涌水量具有明顯的季節(jié)變化特征，隨著礦山的開(kāi)采，需根據(jù)涌水量變化特征及時(shí)調(diào)整疏干排水方案。

（3）barrier障礙邊界在處理控礦構(gòu)造阻水帶對(duì)地下水水平流的阻礙作用方面，模擬效果良好。通過(guò)礦坑涌水量擬合校正，提高模型預(yù)測(cè)的準(zhǔn)確性。

參考文獻(xiàn)：

[1]韓玉杰，陶月贊，周蜜，等.基于數(shù)值法和解析法的礦坑涌水量預(yù)測(cè)對(duì)比分析[J].水利科技及與經(jīng)濟(jì)，2014，10（20）：7-8.

[2]戴巖柯，崔世新，張坤.水均衡法和數(shù)值模擬法在礦坑深部涌水量預(yù)測(cè)中的比較-以西石門(mén)鐵礦為例[J].地下水，2010，1（32）：24-26+48.

[3]李小勇，林堅(jiān)，邱鳳，等.水文地質(zhì)比擬法在礦坑涌水量預(yù)測(cè)中的應(yīng)用[J].資源環(huán)境與工程，2014，28（01）：66-68.

[4]馬洪超，林立新.大井法預(yù)測(cè)礦坑涌水量[J].采礦技術(shù)，2009，9（02）：53-54+79.

[5]郭映忠.復(fù)雜巖溶礦床延深水平涌水量的預(yù)測(cè)[J].地球與環(huán)境，2005，33（3）：19-28.

[6]王世玲，楊作舟.遼寧磚廟硼礦區(qū)二人溝礦段的礦坑涌水量預(yù)測(cè)方法及預(yù)測(cè)結(jié)果[J].化工地質(zhì)，1994，16（4）：269-275.

[7]薛禹群.水文地質(zhì)學(xué)的數(shù)值法[M].北京：煤炭工業(yè)出版社，1986.

[8]賀國(guó)平，張彤，趙月芬，等.GMS數(shù)值建模方法研究綜述[J].地下水，2007（03）：32-35+38.