深部暫停開采礦井采空區(qū)積水過程分析

易偉功

(長沙礦山研究院有限責(zé)任公司， 湖南 長沙 410012)

0 引言

隨著深部資源枯竭而礦區(qū)暫停開采或者礦井停采閉坑，可能會出現(xiàn)如水體污染、地面塌陷、土地破壞、采空區(qū)積水等一系列地質(zhì)災(zāi)害和環(huán)境問題，分析深部暫停開采礦井的采空區(qū)積水過程，對礦井淺部殘留礦體的回收安全具有重要意義[1?2]。國內(nèi)外學(xué)者對煤礦廢棄礦井和采空區(qū)積水的系統(tǒng)模型和積水過程進行了大量的研究。周建軍等[3]通過對廢棄礦井的含水介質(zhì)及其地下水運動規(guī)律進行分析，建立了水流運動概念模型，并提出數(shù)值模擬方法和技術(shù)，為廢棄礦井的恢復(fù)及災(zāi)害治理提供了依據(jù)；張壯路等[4-5]通過對唐村煤礦水文地質(zhì)條件和礦井廢棄后地下水動力模式分析, 分析了水位回彈過程中礦井的容水空間及其分布特點, 分析了廢棄礦井水位回彈的時空過程和最終動態(tài)平衡模式。YOUNGER P.L.等[6]依據(jù)廢棄礦井中井巷工程和采空區(qū)中含水介質(zhì)和水動力場特征的區(qū)別，將礦井含水區(qū)域劃分為紊流管道區(qū)、低滲透性層流區(qū)和高滲透性層流區(qū)，并建立了3種廢棄礦井采空區(qū)積水系統(tǒng)模型；BANKS D.等[7?8]引入有效面積的概念，建立了積水體積可變、水頭可變的地下水回灌模型，并對觀測井中的水位變化進行了成功預(yù)測。獅鳳山銅礦深部資源暫無開采價值，故決定暫停九中段以下所有開采活動，集中回收九中段以上邊角礦體，由于深部排水系統(tǒng)拆除后，地下水會隨著回升，本文以獅鳳山銅礦為研究背景，在綜合考慮國內(nèi)外研究現(xiàn)狀的基礎(chǔ)上，根據(jù)獅鳳山銅礦礦區(qū)水文地質(zhì)條件，分析獅鳳山銅礦深部暫停開采后采空區(qū)的積水過程。

1 深部暫停開采后采空區(qū)積水體積

根據(jù)礦山現(xiàn)有采礦和地質(zhì)資料整理分析，采用回采空間法針對井下主要巷道和采空區(qū)的體積計算出積水體積。積水體積（容積）是指積水標(biāo)高h以下的礦井充水體積的總量。礦井容水空間（介質(zhì)）由井巷工程、采空區(qū)及裂隙帶、未采動巖石裂隙系統(tǒng)組成。其中井巷和采空區(qū)及裂隙帶空間由開采活動形成，是礦井容納積水的主要空間，本次主要是以采空區(qū)容積和開拓巷道的體積之和作為計算的積水體積。

根據(jù)資料，獅山礦段1965年水文地質(zhì)勘探時期測定的地下水位標(biāo)高為1230m～1228m，但由于礦區(qū)幾十年對地下水的排泄和居民用水對地下水的消耗等原因，其積水標(biāo)高在短時間內(nèi)不會高于綠汁江最低侵蝕面（標(biāo)高在1205m左右），根據(jù)對礦區(qū)采空區(qū)調(diào)查及統(tǒng)計結(jié)果，十中段至十八中段的空區(qū)體積為1566.9萬m3，其中采空區(qū)體積為1508.0萬m3，巷道體積為58.9萬m3，見表1。

根據(jù)現(xiàn)場對采空區(qū)的調(diào)查發(fā)現(xiàn)，礦區(qū)中大部分采空區(qū)已被碎石充填，礦區(qū)中采空區(qū)的充填率約為90%，碎石充填體的空隙率可按45%計算，經(jīng)綜合計算，采空區(qū)容積為采空區(qū)體積的50.5%，因此，十中段至十八中段的總積水體積為820.5萬m3，其中采空區(qū)容積為761.6萬m3，巷道積水體積為58.9萬m3，見表2。

表1 各中段體積計算統(tǒng)計

表2 各中段積水體積計算統(tǒng)計

2 積水過程分析

2.1 積水恢復(fù)過程分析

礦井停止排水后，首先是井下最深部積水，積水會以較快的速度升高，達到一定高度后，在水壓差作用下通過不同形式向低標(biāo)高相鄰區(qū)段滲透和徑流，隨著積水高度增加，水力梯度逐漸減小，導(dǎo)水通道增多和滲透面積增大，水位上升的速度越來越慢，最終達到礦區(qū)內(nèi)地下水動態(tài)平衡狀態(tài)。

礦區(qū)采空區(qū)積水的主要來源是昆陽群綠汁江組含水層的補給量約33000m3/d，以及礦區(qū)上部垮落區(qū)雨季補給量約928m3/d。將九中段以下排水系統(tǒng)撤離后，地下水會隨著上升，在一定時間內(nèi)會形成一定標(biāo)高之下的積水靜貯量。由于礦床東部和西南部均為隔水邊界，且昆陽群綠汁江組鳳山段是礦區(qū)的唯一含水層，因此，礦井積水在一定標(biāo)高之下滲透途徑主要有：昆陽群綠汁江組鳳山段白云巖含水層的自然滲透；采礦作用后礦巖柱松動裂隙帶和其他裂隙的滲透；極端條件下軟巖潰破引起的滲透。

礦井在積水初期其補給量遠大于滲透量，積水速率較快，隨著時間的推移，導(dǎo)水通道和滲透的面積也相應(yīng)的增多，積水速率會越來越小，最終當(dāng)?shù)V井補給量與滲漏量相當(dāng)時，礦井積水將會趨于一種動態(tài)的平衡。

2.2 積水恢復(fù)標(biāo)高計算

獅鳳山礦區(qū)主要含水層為潛水含水層，在九中段以下礦井排水設(shè)備拆除后，礦區(qū)地下水只能通過含水層的自然滲透。隨著時間的推移，礦區(qū)地下水的最終潛水位與礦區(qū)原始潛水位基本持平或者略低于原始潛水位，但由于礦井開采的過程中人為對地下水進行了大量排泄，消耗了大量的地下水，整個地下水流場都發(fā)生了很大的變化，在較短時間內(nèi)地下水很難恢復(fù)到原始水位。故本次計算的積水標(biāo)高是在一定時間周期內(nèi)的一個相對值，主要是以礦井內(nèi)滲出和涌入的水量相等時經(jīng)歷的時間值，即基本達到地下水相對動態(tài)平衡狀態(tài)時的相對積水標(biāo)高值作為積水恢復(fù)標(biāo)高。

本文根據(jù)水均衡原理和達西定律計算積水恢復(fù)標(biāo)高，積水恢復(fù)過程中拐點標(biāo)高是指礦井地下水的補給量和排泄量基本相等時，積水高度相對穩(wěn)定的一個動態(tài)平衡，其計算公式如下：

根據(jù)資料分析，目前白云巖水源滲漏穩(wěn)定補給量在33000m3/d，而其他上部垮塌區(qū)的補給量約928m3/d，根據(jù)地下水運動規(guī)律，在較低積水標(biāo)高時白云巖水源滲透補給量可近似保持不變，但二者之和不大于33928m3/d，為了安全起見取33928m3/d，因此有：

式中，dl為含水層厚度，m；k為含水層滲透系數(shù)，m/d；ω(h)為過水?dāng)嗝婷娣e，m2；h為水頭差，m。

九中段以下最大的積水面積為99646m2，假設(shè)拐點積水高度為H，則根據(jù)積分計算可得：

lnH=5.87，則H=354.25m≈354m。

按十八中段最低底板標(biāo)高為745m計算，則積水恢復(fù)過程中拐點標(biāo)高約為1099m，即基本接近于鳳山十一中段（1108m）標(biāo)高。

2.3 水位恢復(fù)時間計算

利用礦山現(xiàn)有資料，統(tǒng)計分析了礦山在開采過程中歷史水位的變化情況（見圖1）。通過圖1可知：礦區(qū)地下水位在礦山開采的初期，由于礦坑涌水量較小，其水位保持較平緩的波動并緩慢下降，隨著礦山開采深度向下延伸，揭露的導(dǎo)水裂隙增多，加大地下水壓，增加相應(yīng)的井下涌水量，地下水位的下降趨勢也明顯地增大，隨著較長時間疏排礦坑涌水，水位曲線逐漸平緩。

圖1 礦區(qū)水位歷時變化曲線

根據(jù)資料統(tǒng)計九中段以下各中段的積水面積（見表3），依據(jù)地下水的運動規(guī)律和礦區(qū)開采過程中的實際情況，在較低標(biāo)高位置其涌水量基本穩(wěn)定，隨著積水高度增高，滲水面積逐漸增加，其礦坑水位回升上漲的速率也呈現(xiàn)出逐漸減少的規(guī)律。利用曲線擬合法確定出積水高度h與時間t的關(guān)系，得出：

其相關(guān)系數(shù)為0.96，積水回升歷時變化曲線如圖2所示。

圖2 積水回升歷時變化曲線

深部礦床暫停開采，排水設(shè)備全部拆除后，根據(jù)圖2可計算得出，礦坑地下水回升到拐點標(biāo)高1099m時需要的時間約為558d；隨著時間的推移，礦區(qū)地下水位將繼續(xù)緩慢上漲，上漲至當(dāng)?shù)厍治g基準面（標(biāo)高為1205m），水位恢復(fù)需要的時間為1650d。另外，以鳳山十八中段最低底板標(biāo)高（745m）為基礎(chǔ)，計算出地下水位恢復(fù)至各中段水平所需時間，見表4。

表3 各中段積水面積

表4 恢復(fù)至各中段水平時間計算

3 結(jié)論

（1）礦井停止排水后礦區(qū)中的水動力學(xué)特征是一個集層流和紊流為一體的復(fù)雜流體。礦床內(nèi)含水介質(zhì)空隙具有多級次的特征，地下水在介質(zhì)中的流動具有時間性較強及階段性較分明的運動規(guī)律。

（2）本次統(tǒng)計計算了十中段至十八中段的積水體積，總積水體積為820.5萬m3，其中采空區(qū)容積為761.6萬m3，巷道積水體積為58.9萬m3。

（3）運用水均衡原理和達西定律相結(jié)合的穩(wěn)定流計算方法，計算得出礦坑地下水回升過程中拐點標(biāo)高約為1099m，基本接近于鳳山十一中段（1108m）標(biāo)高；并計算得出恢復(fù)到該標(biāo)高所需要的時間約為558d；礦坑地下水位恢復(fù)到當(dāng)?shù)厍治g基準面（1205m）時，所需要的時間約為1650d。根據(jù)礦床開采期間地下水位下降的曲線與暫停開采地下水位上升的曲線基本相符的特征，進一步論證計算結(jié)果基本符合實際。