劉先彪，臧彥超，李 靜

(1.神華國能寶清煤電化有限公司 朝陽露天煤礦，黑龍江 雙鴨山 155625；2.大地工程開發(fā)(集團(tuán))有限公司，北京 100102)

煤炭資源開發(fā)工程中疏干降水是安全生產(chǎn)的關(guān)鍵環(huán)節(jié)。礦區(qū)疏干排水方案優(yōu)化一直是礦山開發(fā)研究設(shè)計(jì)面臨主要問題之一[1-5]。近幾十年來，地下水?dāng)?shù)值模擬技術(shù)在疏干排水方案分析研究中得到廣泛應(yīng)用。例如，陳時(shí)磊等開展了基于礦井生產(chǎn)進(jìn)度疏干條件下的三維地下水?dāng)?shù)值模擬研究，分析滿足生產(chǎn)進(jìn)度疏干條件礦井涌水量[6]；白樂等研究了煤礦開采區(qū)地表水-地下水轉(zhuǎn)化關(guān)系，開展了地表水-地下水耦合模擬研究[7]；馮翔開展了利用GMS地下水模擬軟件，開發(fā)礦井三維地下水系統(tǒng)模型，進(jìn)行礦井涌水量預(yù)測[8]。對(duì)于復(fù)雜多含水層系統(tǒng)，以及不同疏干排水方式的模擬與方案優(yōu)化尚不多見。

朝陽露天礦區(qū)含水層呈多層結(jié)構(gòu)，滲透性差，不易疏干，地下水位(壓)較高，采掘場出現(xiàn)多次滑坡，為保證露天礦安全生產(chǎn)，迫切需要準(zhǔn)確預(yù)測露天礦疏干水量，采取有效疏干方式降低地下水位。本文對(duì)礦區(qū)水文地質(zhì)資料和疏干排水?dāng)?shù)據(jù)詳細(xì)分析，運(yùn)用Visual Modflow軟件，開發(fā)了礦區(qū)地下水三維數(shù)值模型，提出了地面降水孔降水、采場內(nèi)水平放水孔排水、泄壓坑排水等多種方式聯(lián)合的綜合疏干降水優(yōu)化方案，為礦區(qū)確定治水方案和措施提供科學(xué)依據(jù)[9-16]。

1 礦區(qū)水文地質(zhì)條件

礦區(qū)地處三江平原的南部，沉積有第三系和第四系地層。露天礦開采的對(duì)象為10煤層，依據(jù)含水層和煤層空間分布關(guān)系，本區(qū)含水層分為上部含水層、中部含水層和下部含水層，煤層位于中部和下部弱含水層之間，平均厚度15m。上部含水層組：全區(qū)發(fā)育，巖性以中砂和細(xì)砂為主，厚度變化大，平均厚度11.75m；中部含水層組：全區(qū)發(fā)育，巖性主要為第三系微膠結(jié)砂巖及第四系底部中粗砂，平均厚度21m；下部含水層組：發(fā)育在10煤底板以下，巖性為第三系微膠結(jié)中、細(xì)砂巖和粉砂巖，平均厚度5～6m。礦區(qū)發(fā)育三個(gè)相對(duì)弱含水層，即頂部、中部及底部弱含水層，平均厚度分別6.76m、8m和2～3m。頂部弱含水層裸露地表，中部弱含水層埋深在30m左右，底部弱含水層埋藏在下部含水層與煤層底板之間。弱含水層主要為粘土、砂質(zhì)粘土和弱膠結(jié)泥巖、砂質(zhì)泥巖。上部含水層和中部含水層之間的弱含水層連續(xù)分布，而中、下部含水層之間弱含水層由于發(fā)育有斷層或天窗，含水層存在一定的水力聯(lián)系。

礦區(qū)地下水主要接受大氣降水入滲、地表水和灌溉渠系滲漏、灌溉田間回歸，以及山前側(cè)向徑流補(bǔ)給。天然條件下，地下水由西南向東北方向徑流，消耗于潛水蒸發(fā)和側(cè)向徑流。

2 礦區(qū)地下水?dāng)?shù)值模型

2.1 含水層和邊界條件概化

礦區(qū)研究范圍界定在以礦區(qū)為中心，東和南到山區(qū)山前，西到地表分水嶺、北到平原區(qū)山丘隆起帶，模型區(qū)面積為114.72km2。

依據(jù)巖性與富水性，含水層自上而下概化為7層：第1層為第四系亞粘土、粘土層相對(duì)弱含水層；第2層即第四系中砂細(xì)砂含水層(上含水層)，滲透系數(shù)3.1～4.1m/d；第3層為第四系粘土、砂質(zhì)粘土相對(duì)弱含水層；第4層為第四系底部中粗砂和第三系微膠結(jié)砂巖強(qiáng)含水層(中含水層)，滲透系數(shù)2.31～2.98m/d，含水層與煤層頂板直接接觸；第5層為相對(duì)弱含水層，即煤系地層；第6層為相對(duì)弱含水層；第7層為第三系微膠結(jié)中、細(xì)砂巖強(qiáng)含水層(下含水層)，滲透系數(shù)3.6～5.2m/d，是礦床直接充含水層。各個(gè)分層的頂?shù)装逡罁?jù)鉆孔資料，進(jìn)行空間插值得到。

頂部邊界接收大氣降水和地表水體的入滲補(bǔ)給，采用入滲系數(shù)法計(jì)算。四周邊界的概化為流量邊界。為反映邊界變化的特點(diǎn)概化為混合邊界(MODFLOW中的GHB邊界類型)，邊界流量與邊界內(nèi)外水位及含水層導(dǎo)水性有關(guān)。

礦區(qū)范圍內(nèi)，深埋的基底第三系泥巖為相對(duì)隔水層，底部邊界概化為隔水邊界。

2.2 數(shù)學(xué)模型及求解

含水層概化為非均質(zhì)各項(xiàng)異性，地下水水流概化為三維非穩(wěn)定流，主微分方程：

式中，kxx，kyy，kzz分別為沿主軸方向滲透系數(shù)，m/d；h為測壓水頭，m；Q為單位體積流量(降水或地表水入滲、人工抽水等)，(m3/d)/m3；Ss為單位儲(chǔ)水系數(shù)，m-1；t為時(shí)間，d。

初始及邊界條件：

h(x，y，z，t)∣t=0=h0(x，y，z)

式中，h0為含水層的初始水位分布，m；Γ12為滲流區(qū)域的側(cè)向混合邊界；q(x，y，t)為邊界流量，采用GHB計(jì)算方法得到。

設(shè)計(jì)降水孔，按分布和開采深度(濾水管)布置于各單元，隨水位下降，降水孔抽水量也將衰減。降水孔抽水量與水位以下的濾水管長度和揭露含水層導(dǎo)水性有關(guān)。計(jì)算公式如下：

式中，Qi為在某時(shí)段內(nèi)從第i層進(jìn)入某井的水量；Qt為在該時(shí)段內(nèi)的井的抽水量；Li為第i層中濾水管段長度；Kx為第i層X方向的滲透系數(shù)。

采場排水概化為Drain邊界類型。采坑排泄地下水量計(jì)算方法如下：

Qk=(h-hc)C

式中，Qk為采坑排泄地下水量，m3/d；h為地下水位標(biāo)高，m；hc為坑內(nèi)水位高程，坑內(nèi)無水時(shí)，為坑底高程，m；C為綜合導(dǎo)水系數(shù)，m2/d。

當(dāng)?shù)叵滤坏陀诳拥赘叱虝r(shí)，不再排泄地下水。

計(jì)算區(qū)域采用40m×40m的矩形網(wǎng)格，在開采區(qū)和水文地質(zhì)條件變化較大的地方加密部分網(wǎng)格(20m×20m)，以提高計(jì)算的精度。為了保證收斂的同時(shí)提高計(jì)算速度，對(duì)網(wǎng)格作平滑處理。

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

采用2008年11月26日—12月10日的群孔抽水試驗(yàn)數(shù)據(jù)進(jìn)行模型參數(shù)和邊界條件識(shí)別，抽水試驗(yàn)用162C2和162C3共2個(gè)井抽水，5個(gè)觀測井同步觀測水位。2個(gè)抽水井水位模擬和實(shí)測水位過程曲線如圖1所示。

圖1 2個(gè)抽水井水位模擬和實(shí)測水位過程曲線

運(yùn)用2019年1月1日—4月30日的孔排、坑排水量和觀測井水位數(shù)據(jù)，進(jìn)行模型驗(yàn)證。收集的6個(gè)觀測孔實(shí)測數(shù)據(jù)模擬結(jié)果進(jìn)行置信檢驗(yàn)，結(jié)果如圖2所示，均差-0.76m，相關(guān)系數(shù)為0.976，RMS為8.52%，滿足模型精度要求。模型驗(yàn)證期內(nèi)，地下水總補(bǔ)給量和總排泄量分別為0.91萬m3/d和1.98萬m3/d，靜儲(chǔ)量變化量為-1.07萬m3/d?？着藕涂优欧謩e為1.49萬m3/d和0.36萬m3/d?？偱潘繛?.85萬m3/d，期間消耗靜儲(chǔ)量為-1.07萬m3/d，占總排水量的57.84%。

圖2 計(jì)算和實(shí)測水位擬合曲線

3 降水方案模擬預(yù)測

3.1 降水方案

降水方案主要目的是在充分利用已有降水孔和坑內(nèi)疏干排水設(shè)施的基礎(chǔ)上，通過新增降水設(shè)施，有效降低采掘場內(nèi)地下水位，保證采剝工作順利開展，提高剝離、采煤工作效率。同時(shí)，最大程度降低非工作幫和東、西端幫的地下水位，減少邊幫地下水滲出量，提高邊幫穩(wěn)定性。擬定方案如下：

1)利用89個(gè)已有降水孔，2020年新投入60～90個(gè)設(shè)計(jì)地面降水孔。

2)2020—2021年在工作幫臺(tái)階分階段設(shè)置100～120個(gè)水平放水孔(傾角5°)，疏放采場工作幫的中含水層地下水。

3)隨著開采工作面推進(jìn)，在采場坑底煤層底板標(biāo)高較低地帶，開挖泄壓坑(井)，降低下含水層地下水位(水頭)。

3.2 模擬預(yù)測與方案優(yōu)化

降水方案模擬預(yù)測時(shí)，大氣降水取多年平均降水量，保持不變。用2018年底流場為初始流場，時(shí)間步長為30d，預(yù)測到2023年12月30日。

初始設(shè)計(jì)地面降水孔和水平孔的數(shù)量分別為45個(gè)和80個(gè)，通過逐步增加設(shè)計(jì)的地面降水孔和水平孔的數(shù)量，預(yù)測各含水層水位動(dòng)態(tài)，直到水位降滿足疏干和降壓水位要求，即得到優(yōu)化排水方案。通過模擬確定的排水方案為68個(gè)地面孔和120個(gè)水平孔。

2019年總排水量1.79萬m3/d，2020年啟用北部、東北和西北部67個(gè)降水孔后排水量急劇增加，達(dá)到3.4萬m3/d，之后幾年呈減少趨勢，2023年排水量為2.41萬m3/d，2019—2023年礦區(qū)總排水量過程預(yù)測曲線如圖3所示；相應(yīng)消耗靜儲(chǔ)量由2020年2.53萬m3/d，減少到0.98萬m3/d。排水量組成中，靜儲(chǔ)量占42%，其他來自側(cè)向補(bǔ)給、降水和地表水入滲補(bǔ)給。按分層地下水組成看，排水主要來自中含水層和下含水層，分層含水層排水量變化過程見表1，上含水層、中含水層和下含水層排水量分別占6%、61%和33%。下含水層的排水需要孔排和坑排同時(shí)進(jìn)行，坑排水量主要來自下含水層的泄壓排水。

圖3 2019—2023年礦區(qū)總排水量過程預(yù)測曲線

表1 分層含水層排水量變化過程 萬m3/d

2023年中含水層和下含水層預(yù)測地下水流場如圖4所示。隨著開采向北推進(jìn)，采坑范圍中含水層逐步疏干，外圍水位在40～45m變化。由于外圍的孔排和坑內(nèi)泄壓坑排水，坑內(nèi)下含水層不同年份水位(水頭)能夠控制在煤層底板之下，滿足了煤礦安全水位(水頭)要求，預(yù)測礦區(qū)下含水層水位(頭)標(biāo)高見表2。坑外圍水位(水頭)標(biāo)高控制在15m。

圖4 預(yù)測期末中含水層和下含水層水位標(biāo)高等值線(m)

表2 預(yù)測礦區(qū)下含水層水位(頭)標(biāo)高 m

4 結(jié) 論

1)在礦區(qū)水文地質(zhì)模型基礎(chǔ)上，建立了三維地下水流數(shù)學(xué)模型。運(yùn)用已有疏干降水資料和地下水水位觀測數(shù)據(jù)，進(jìn)行模型參數(shù)與邊界條件識(shí)別與驗(yàn)證以及地下水均衡分析。結(jié)果表明，開發(fā)的模型是可靠的。目前排水條件下，疏干排水主要是來自中含水層和下含水層，分別占61%和33%。地下水均衡為確定疏干降水方案增加了堅(jiān)實(shí)的依據(jù)。

2)運(yùn)用開發(fā)的地下水?dāng)?shù)值模擬模型，預(yù)測分析了考慮“地面孔排—臺(tái)階水平孔—坑底泄壓坑”等綜合排水方案下的礦坑涌水量及疏干時(shí)間。預(yù)測最大排水量為3.9萬m3/d，排水量隨工作面推進(jìn)呈緩慢減小趨勢，排水水源來自周邊地下水側(cè)向補(bǔ)給、消耗靜儲(chǔ)量以及降水地表水入滲補(bǔ)給，其中消耗靜儲(chǔ)量占42%。

3)根據(jù)礦區(qū)水文地質(zhì)條件研究和數(shù)值模擬成果，優(yōu)化了降水方案。模擬模型預(yù)測結(jié)果表明，綜合疏干降水方案是可行的。采用優(yōu)化后的降水方案，地下水位(頭)下降顯著，從而有效保證露天礦正常開采。