易偉功

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

0 引言

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

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

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

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

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

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

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

2 積水過(guò)程分析

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

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

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

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

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

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

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

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

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

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

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

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

2.3 水位恢復(fù)時(shí)間計(jì)算

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

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

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

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

圖2 積水回升歷時(shí)變化曲線

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

表3 各中段積水面積

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

3 結(jié)論

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

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

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