王 海,王永剛,張 雁,孫 浩,苗賀朝
(1.中煤科工集團(tuán)西安研究院有限公司,陜西 西安 710077;2.中國礦業(yè)大學(xué) 資源與地球科學(xué)學(xué)院,江蘇 徐州 221116;3.內(nèi)蒙古平莊煤業(yè)(集團(tuán))有限責(zé)任公司,內(nèi)蒙古 赤峰 024050)
我國噸煤開采產(chǎn)生礦井水為1.87 m3,每年產(chǎn)生礦井水約6.88×109m3[1]。因此,礦井排水對(duì)礦區(qū)水資源和生態(tài)環(huán)境具有極為不利的影響,尤其是對(duì)于生態(tài)環(huán)境脆弱的露天礦區(qū)。蒙東地區(qū)氣候干旱,降水量少,水資源匱乏,生態(tài)環(huán)境脆弱,煤炭高強(qiáng)度開采引起地下水系統(tǒng)的補(bǔ)、徑、排條件發(fā)生劇烈變化,產(chǎn)生松散層水位大幅下降、地表河流量衰減、生態(tài)環(huán)境破壞等一系列難題[2-4]。由于本區(qū)域的露天礦生態(tài)閾值較低,抗擾動(dòng)能力差,傳統(tǒng)的疏排水方式將大范圍和大幅度地疏降地下水,導(dǎo)致大面積的喬、灌、草等植被衰敗減少,草場(chǎng)退化,加速了荒漠化進(jìn)程[2-3]。為有效減少礦坑疏排水量,抑制生態(tài)脆弱露天礦區(qū)松散層水位急劇下降的趨勢(shì),保護(hù)露天礦區(qū)并不豐富的地下水資源,防止因缺水而導(dǎo)致生態(tài)環(huán)境進(jìn)一步惡化[5-6],近年來,帷幕截水技術(shù)被引入到露天礦防治水工作中,通過截水帷幕優(yōu)化、轉(zhuǎn)變露天礦地下水控制方式,阻斷礦坑和周邊區(qū)域松散層含水層水力聯(lián)系,保護(hù)第四系松散層水資源,減少露天礦開采活動(dòng)對(duì)區(qū)域水環(huán)境影響,保障礦山的安全生產(chǎn)[7-8]。
截水帷幕在內(nèi)蒙古扎尼河和元寶山露天煤礦取得良好的應(yīng)用效果,降低了礦坑疏干水量,減少了第四系松散層水對(duì)礦坑的補(bǔ)給,帷幕內(nèi)外側(cè)水位也發(fā)生相應(yīng)的變化[9-10]。但生態(tài)脆弱露天礦區(qū)截水帷幕下松散層水位的演化特征和變化規(guī)律尚不清晰,且相關(guān)研究文獻(xiàn)極少?;诖?,筆者以內(nèi)蒙古扎尼河和元寶山露天煤礦實(shí)施的2 個(gè)截水帷幕為研究對(duì)象,分析松散層水位演化規(guī)律,以期為生態(tài)脆弱區(qū)礦山水資源保護(hù)提供參考。
我國西部富煤缺水區(qū)約占國土總面積的 46.1%,該地區(qū)大部分處于缺水帶、少水帶和過渡帶,淡水資源僅占我國的10.2%[11]。生態(tài)脆弱露天礦區(qū)地層結(jié)構(gòu)如圖1 所示,煤層上方為粉砂質(zhì)泥巖和炭質(zhì)泥巖隔水層、黏土隔水層,含水層為砂卵石松散層。露天礦區(qū)煤炭開采導(dǎo)致砂卵石松散層水源源不斷補(bǔ)給礦坑。
圖1 生態(tài)脆弱區(qū)露天煤礦地層Fig.1 Strata of open-pit mining area in ecologically fragile areas
蒙東礦區(qū)位于內(nèi)蒙古東部,屬干旱半干旱氣候,年降水量小于400 mm,生態(tài)環(huán)境極其脆弱,元寶山、伊敏、大雁、霍林河等內(nèi)蒙古東部煤田的開發(fā)形成了多個(gè)大型露天礦區(qū)[12-14]。蒙東礦區(qū)第四系松散含水層為砂卵石層,滲透性好,粒徑1~10 cm,次圓狀,分選好,含少量細(xì)砂,滲透系數(shù)為80~700 m/d,水位埋藏較淺,一般在0~16 m,大多接受地表河流補(bǔ)給[12],是區(qū)域內(nèi)各露天礦生產(chǎn)過程中疏排水的主要補(bǔ)給水源。
如圖2 所示,露天煤礦開采煤層埋藏淺,充水水源固定,主要采用疏排水方法進(jìn)行礦坑水疏降或疏干,具體方式有第四系疏干孔、煤層疏干孔、泄水巷、強(qiáng)排孔、地面截水溝等[7]。
圖2 露天煤礦地下水控制方式[7]Fig.2 Groundwater control methods in open-pit coal mines[7]
采用預(yù)先疏干方式降低地下水位時(shí),應(yīng)根據(jù)采剝推進(jìn)和開拓階段強(qiáng)度要求,結(jié)合水文地質(zhì)條件,確定預(yù)先疏干時(shí)間和水位的降低深度。
采用地面垂直降水孔法的第四系和煤層疏干孔適用于滲透系數(shù)大于2 m/d 的含水層;采用水平放水孔法的第四系和煤層疏干孔適用于地下水補(bǔ)給條件較差、降低露天煤礦邊坡地下水壓或排放含水層殘余水的條件;明溝和暗溝等地面截水溝法適用于水文地質(zhì)條件簡(jiǎn)單、埋深較淺、厚度較小且產(chǎn)狀較穩(wěn)定的松散含水層;泄水巷自流排水適用于水文地質(zhì)條件復(fù)雜、水力聯(lián)系小的多含水層,或含水層厚度、水壓及透水性變化較大、埋藏較深且不適用降水孔疏干的情況[15]。
例如,元寶山露天煤礦采用以“群井疏干為主、坑下明排為輔、基巖疏干巷道為補(bǔ)充”的聯(lián)合疏干系統(tǒng)以及地面防排水系統(tǒng)。群井疏干系統(tǒng)采用環(huán)形加線形的布設(shè)方式,目的在于將地下水疏降工作超前采掘進(jìn)度,保證生產(chǎn)正常進(jìn)行和可持續(xù)發(fā)展。集水坑相當(dāng)于一個(gè)大口徑的集中疏干井,并以排水溝作為擴(kuò)大影響范圍的一種方法解決大氣降水、采場(chǎng)殘余水,同時(shí)起到預(yù)先疏干降水的作用。泄水巷道位于煤層頂界面下方2~3 m 的煤層中,坡度5°,下部近剝蝕區(qū)處落平,形成主水倉,在主水倉上部施工排水井。通過多種地下水控制方法,將露天煤礦地下水位控制在采掘工作面以下,保證安全、高效采掘。
露天開采總體上引起礦區(qū)周邊地下水位下降,本文以蒙東地區(qū)伊敏、東明和元寶山露天煤礦為例,分析露天開采對(duì)松散層水位的影響。
伊敏露天煤礦疏干導(dǎo)致礦坑及周圍地下水位持續(xù)下降,水位埋深由0~16 m 大幅下降至 0~112 m,湖泊總面積由采動(dòng)前的6.94 km2縮小為1.12 km2,面積減少率達(dá)84%,井田范圍內(nèi)所有湖泊、泉水全部干涸,濕地消失[12-13]。東明露天煤礦疏排水量1.493×105m3/d,疏干了礦坑周邊大范圍松散層水,中心水位下降幅度在100 m 以上[5]。元寶山露天煤礦早期疏干水量約4.676×105m3/d,目前疏干水量1.8×105m3/d 左右,地下水位最大降深49 m。
蒙東礦區(qū)具有特殊的自然、生態(tài)和地質(zhì)環(huán)境,煤層采動(dòng)過程改變了含水層的賦存狀態(tài),地下水賦存、補(bǔ)給、循環(huán)模式發(fā)生改變,松散層水位下降,濕生植被向旱生植被演變,加劇土壤鹽漬化和荒漠化問題,嚴(yán)重影響生態(tài)環(huán)境質(zhì)量[11]。
截水帷幕有別于現(xiàn)有的露天礦區(qū)疏降或疏干方法,可保護(hù)礦區(qū)松散層水資源,恢復(fù)礦區(qū)周邊松散層地下水位。
通過探明煤礦區(qū)水文地質(zhì)工程地質(zhì)條件,厘清礦坑水的補(bǔ)給水源和補(bǔ)給通道,確定適宜的截水帷幕路線和墻體空間結(jié)構(gòu),確保帷幕底部進(jìn)入穩(wěn)定隔水層、頂部高于歷史水位,因地制宜地選擇溝槽式截水帷幕或鉆孔注漿帷幕方式及其配套防滲材料,通過構(gòu)筑連續(xù)的帷幕墻體攔截松散層水側(cè)向補(bǔ)給礦坑。
如圖3 所示,將第四系疏干孔、煤層疏干孔、泄水巷、強(qiáng)排孔、地面截水溝等礦坑水疏降或疏干控制礦區(qū)地下水方式改變?yōu)榻厮∧环绞剑瑢⑺缮铀c礦坑隔離,僅有少量松散層水經(jīng)過截水帷幕滲透進(jìn)入礦區(qū)。截水帷幕外側(cè)松散層水位逐漸抬升、恢復(fù),而截水帷幕內(nèi)側(cè)水位迅速下降。
圖3 截水帷幕作用下松散層水位Fig.3 Water level of the unconsolidated strata under the condition of water cutoff curtain
將圖3 簡(jiǎn)化為圖4 所示的截水帷幕內(nèi)外側(cè)松散層水位計(jì)算模型。
圖4 截水帷幕內(nèi)外松散層水位計(jì)算簡(jiǎn)化模型Fig.4 Simplified model for calculating water level of loose layer inside and outside the water cutoff curtain
對(duì)于帷幕外側(cè)的松散層,單寬流量q為:
式中:K1為截水帷幕外側(cè)含水層的滲透系數(shù);h1為距離帷幕外側(cè)L1處松散層水位;h2為帷幕外側(cè)松散層水位;L1為水位h1、h2之間的距離。
將式(1)轉(zhuǎn)化為:
通過式(2)可推算出距離帷幕外側(cè)L1處的松散層水位。
同理,對(duì)于截水帷幕,單寬流量q為:
式中:K2為截水帷幕的滲透系數(shù);h3為帷幕內(nèi)側(cè)松散層水位;D為截水帷幕的厚度。
將式(3)轉(zhuǎn)化為:
同理,對(duì)于帷幕內(nèi)側(cè)的松散層,單寬流量q為:
式中:K3為截水帷幕內(nèi)側(cè)含水層的滲透系數(shù);h4為距離帷幕內(nèi)側(cè)L2處松散層水位;L2為水位h3、h4之間的距離。
將式(5)轉(zhuǎn)化為:
當(dāng)截水帷幕內(nèi)外側(cè)地層滲透系數(shù)一致時(shí),即K1=K3,則式(6)變?yōu)椋?/p>
當(dāng)q、K1、K2、K3、D、L1、L2、h2(或h3)值已知時(shí),可得出截水帷幕內(nèi)外側(cè)不同位置的松散層水位。
假定截水帷幕內(nèi)外側(cè)松散層滲透系數(shù)相同,即K1=K3,帷幕滲透系數(shù)K2為松散層滲透系數(shù)K1的0.45%,帷幕截水效率為80%,采用定水頭補(bǔ)給方式,分別應(yīng)用式(1)-式(7)計(jì)算截水帷幕作用下的不同位置松散層水位,得到如圖5 所示的計(jì)算結(jié)果。無截水帷幕條件下,松散層水位隨著與定水頭補(bǔ)給源的距離增大而不斷降低;有截水帷幕時(shí),帷幕外側(cè)水位抬升明顯,而帷幕內(nèi)側(cè)水位急劇降低,帷幕內(nèi)外水位差極大,充分體現(xiàn)了截水帷幕的正向作用。
圖5 有無截水帷幕下松散層水位變化的計(jì)算結(jié)果Fig.5 Calculation results of water level changes of the unconsolidated strata with or without the water cutoff curtain
2.2.1 物理模型
根據(jù)露天煤礦水文地質(zhì)條件建立如圖6 所示的松散層截水帷幕模型,在松散層模型中增加不同長度的垂向隔水板模擬不同圈閉范圍的截水帷幕,記錄帷幕兩側(cè)水位變化數(shù)據(jù),分析帷幕構(gòu)筑前后第四系松散層水位變化規(guī)律。
圖6 松散層截水帷幕模型縱向剖面Fig.6 Longitudinal profile of the water cutoff curtain model in the unconsolidated strata
松散層截水帷幕模型由3 部分組成,分別為水源倉、滲流倉和排泄倉。模型尺寸為3 000 mm(長)×1 500 mm(寬)×600 mm(高),其中水源倉與排泄倉長度均為300 mm,滲流倉長度為2 400 mm。水源倉用于模擬流場(chǎng)補(bǔ)給源并保持穩(wěn)定水頭。靠近滲流倉的一側(cè)布置滲流孔,通過滲流孔向滲流倉滲水,外側(cè)垂向布置溢出孔保持穩(wěn)定水位。供水水源采用飲用水,使用流量表監(jiān)測(cè)補(bǔ)給流量[16]。
如圖7 所示,滲流倉模擬截水帷幕下松散層水滲流過程,滲流倉中充填滲透材料,底部安設(shè)水位觀測(cè)孔,通過玻璃管連接至滲流倉外壁進(jìn)行水位觀測(cè),滲流倉中部采用不透水拼接擋板模擬截水帷幕,通過改變擋板的長度,模擬不同過水?dāng)嗝娴臐B流過程,擋板底部與滲流倉底部接觸部位進(jìn)行隔水處理[16]。排泄倉收集滲流水,在排泄倉底部布置滲流水收集孔,將滲流水引至儲(chǔ)水容器并進(jìn)行滲流量觀測(cè)。
圖7 松散層截水帷幕試驗(yàn)?zāi)P虵ig.7 Test model of the water cutoff curtain in the unconsolidated strata
2.2.2 水位模擬結(jié)果
如圖8 所示,假設(shè)截水帷幕總長度為L,通過圖7所示的松散層截水帷幕試驗(yàn)?zāi)P?,模擬帷幕長度為0、0.38L、0.51L、0.65L、0.79L、0.86L、0.93L和1.00L等不同條件下的帷幕內(nèi)外側(cè)松散層水位變化過程。截水帷幕位于1 200 mm 位置處,從位于0 mm 位置的水源倉進(jìn)水開始直至流場(chǎng)穩(wěn)定,水源倉初始水位為500 mm,觀測(cè)整個(gè)流場(chǎng)變化過程,2 400 mm 位置留設(shè)排泄倉。
圖8 不同類型截水帷幕作用下的松散層水位變化Fig.8 Changes of water level of the unconsolidated strata under different water cutoff curtain conditions
由圖8 可知,無帷幕時(shí),松散層水位連續(xù)變化,從水源倉至排泄倉范圍的松散層水位逐漸降低,水位由500 mm 降至87 mm,截水帷幕位置無明顯水位差。帷幕長度0.38L時(shí),截水帷幕外側(cè)松散層水位由221 mm抬升至223 mm,而截水帷幕內(nèi)側(cè)松散層水位由209 mm下降至206 mm,截水帷幕位置松散層水位差略有顯現(xiàn)。帷幕長度0.65L時(shí),截水帷幕外側(cè)松散層水位進(jìn)一步抬升至230 mm,外側(cè)松散層水位得到一定保護(hù),而帷幕內(nèi)側(cè)松散層水位降至195 mm,下降明顯,帷幕所在位置的松散層水位差為35 mm。帷幕長度0.93L時(shí),截水帷幕外側(cè)松散層水位抬升至258 mm,帷幕外側(cè)松散層水位得到較好保護(hù),而帷幕內(nèi)側(cè)松散層水位下降至143 mm,帷幕所在位置的松散層水位差較大。當(dāng)截水帷幕形成圈閉時(shí),即帷幕長度1.00L時(shí),帷幕外側(cè)松散層水位迅速抬升至462 mm,與水源倉水位基本持平,外側(cè)松散層水位得到有效保護(hù),帷幕內(nèi)側(cè)松散層水位下降至17.6 mm,帷幕所在位置的松散層水位差巨大。
由圖5 和圖8 截水帷幕下松散層水位變化理論計(jì)算和物理模擬結(jié)果可見,理論計(jì)算結(jié)果和物理模擬結(jié)果的規(guī)律一致,距離截水帷幕較遠(yuǎn)的松散層水位變化速率為0.01~0.03 mm/mm,較為穩(wěn)定,靠近截水帷幕的松散層水位變化相對(duì)劇烈。滲流倉末端松散層穩(wěn)定水位隨著截水帷幕長度的減少呈逐漸抬升趨勢(shì),靠近截水帷幕的松散層水位變化幅度比遠(yuǎn)離截水帷幕的水位變化幅度劇烈。受截水帷幕影響,在帷幕靠近水源倉一側(cè)會(huì)形成高水位區(qū),在帷幕靠近排泄倉一側(cè)會(huì)形成低水位區(qū),且隨著帷幕長度減小,高(低)水位區(qū)域逐漸減小,且水位差在逐漸降低[17-19]。
截水帷幕作用下典型礦區(qū)松散層水位演化規(guī)律分析選擇蒙東地區(qū)的扎尼河和元寶山露天煤礦為研究對(duì)象,分別闡述礦區(qū)和截水帷幕概況,分析帷幕截水效果以及松散層水位演化特征。
3.1.1 截水帷幕概況
扎尼河露天煤礦砂卵石層滲透性好,滲透系數(shù)達(dá)80~180 m/d,礦坑疏排降水破壞周邊草原生態(tài)環(huán)境、造成水資源浪費(fèi)。為減少露天煤礦礦坑疏排水量、保護(hù)草原水資源和生態(tài)環(huán)境,采用截水帷幕替代傳統(tǒng)的露天煤礦疏排降水方式[6-7,16]。如圖9 所示,扎尼河露天煤礦截水帷幕分別應(yīng)用了防滲膜、抗?jié)B混凝土地下連續(xù)墻、超高壓角域變速射流注漿、鉆孔咬合樁等工藝,形成了“頂不越、底不漏、兩端不繞、接頭不滲”的5 815 m松散層地下截水帷幕[7]。
3.1.2 松散層水位演化規(guī)律
如圖9 所示,在扎尼河露天煤礦截水帷幕內(nèi)側(cè)布設(shè)5 個(gè)松散層水位觀測(cè)孔,編號(hào)為K1-K5,在帷幕外側(cè)布設(shè)4 個(gè)松散層水位觀測(cè)孔,編號(hào)K6-K9。扎尼河露天煤礦截水帷幕工程自2018 年7 月開工,如圖10所示,截水帷幕外側(cè)K6-K9 孔處松散層初始水位高程分別為617.95、618.68、618.41、617.86 m,帷幕內(nèi)側(cè)松散層水位觀測(cè)孔K1-K5 的水位分別為618.39、618.06、616.02 和617.80 m。
圖9 扎尼河露天煤礦截水帷幕及水位觀測(cè)孔平面分布Fig.9 Schematic plan of the water cutoff curtain at Zanihe open-pit coal mine
由圖10 和表1 可知,隨著扎尼河露天煤礦截水帷幕完成比例的增加,帷幕外側(cè)松散層水位逐步抬升,而帷幕內(nèi)側(cè)松散層水位不斷下降。帷幕工程完成100%時(shí),K6—K9 孔水位分別變化了2.55、1.93、3.16、3.13 m,K1—K5 孔的水位分別變化了-6.55、-3.49、-4.02、-5.51 和-6.53 m。在帷幕完成后,帷幕外側(cè)松散層水位仍在抬升,內(nèi)側(cè)水位將進(jìn)一步降低,帷幕內(nèi)外側(cè)水頭差增大。
圖10 扎尼河露天煤礦松散層水位變化Fig.10 Changes of water level of the unconsolidated strata inside and outside the curtain at Zanihe open-pit coal mine
表1 扎尼河露天煤礦松散層水位變化規(guī)律Table 1 Changes of water level of the unconsolidated strata inside and outside the curtain at Zanihe open-pit coal mine
3.2.1 截水帷幕概況
元寶山露天煤礦礦坑主要充水水源為第四系松散層水,占礦坑疏排水量95%左右,地層主要為圓礫夾砂、卵石層、黃土及風(fēng)化砂巖、泥巖,圓礫石粒徑為50~150 mm,滲透系數(shù)17~700 m/d。歷史上礦坑的疏排水量每天達(dá)數(shù)十萬立方米,疏排水量以第四系松散含水層的靜儲(chǔ)量為主,動(dòng)態(tài)補(bǔ)給量為輔。隨著靜儲(chǔ)量的持續(xù)疏放,動(dòng)態(tài)補(bǔ)給量占比逐年增大,礦坑疏排水量呈逐年下降并漸趨穩(wěn)定。
為減少礦坑疏排水量,擬沿現(xiàn)有礦坑一級(jí)邊坡平臺(tái)建造礦坑邊界全封閉截水帷幕,在所有來水邊界進(jìn)行截水,可最大程度減少礦坑疏排水量。如圖11 所示,先期在礦坑南側(cè)開展長度1 369 m 的帷幕試驗(yàn)工程,深度16~43 m。
3.2.2 松散層水位演化規(guī)律
如圖11 所示,在元寶山露天煤礦截水帷幕試驗(yàn)工程外側(cè)布設(shè)4 個(gè)松散層水位觀測(cè)孔,編號(hào)Z1-Z4,在帷幕內(nèi)側(cè)布設(shè)2 個(gè)松散層水位觀測(cè)孔,分別為Z5、Z6。隨著截水帷幕長度的增加,帷幕外側(cè)松散層水位不斷抬升,而帷幕內(nèi)側(cè)水位逐漸降低,帷幕內(nèi)外側(cè)水位變幅與帷幕完成率呈正相關(guān)關(guān)系(圖12)。
圖11 元寶山露天煤礦截水帷幕及水位觀測(cè)孔平面布置Fig.11 Arrangement plan of the water cutoff curtain at Yuanbaoshan open-pit coal mine
圖12 元寶山露天煤礦帷幕內(nèi)外側(cè)松散層水位變化Fig.12 Changes of water level of the unconsolidated strata inside and outside the curtain at Yuanbaoshan open-pit coal mine
由表2 可知,元寶山露天煤礦截水帷幕試驗(yàn)工程完成7 d 時(shí),帷幕外側(cè)松散層水位觀測(cè)孔Z1-Z4 的水位分別較原水位升高9.12、6.48、8.18 和5.55 m;帷幕內(nèi)側(cè)松散層水位觀測(cè)孔Z5、Z6 的水位分別較原水位降低3.15、5.82 m。
表2 元寶山露天礦松散層水位變化規(guī)律Table 2 Changes of water level of the unconsolidated strata inside and outside the curtain at Yuanbaoshan open-pit coal mine
a.我國蒙東地區(qū)生態(tài)環(huán)境脆弱,水資源匱乏,露天開采引起地下水系統(tǒng)的補(bǔ)、徑、排條件發(fā)生劇烈變化,導(dǎo)致松散層水位大幅度下降、地表生態(tài)環(huán)境破壞。截水帷幕改變了傳統(tǒng)的地下水控制方式,可減少礦坑疏排水量,抑制煤礦區(qū)松散層水位急劇下降的趨勢(shì),保護(hù)煤礦區(qū)松散層地下水資源。
b.理論計(jì)算、物理模擬結(jié)果與扎尼河和元寶山露天煤礦截水帷幕現(xiàn)場(chǎng)監(jiān)測(cè)結(jié)果一致,截水帷幕阻隔內(nèi)外側(cè)松散層水力聯(lián)系,靠近排泄一側(cè)形成低水位區(qū),帷幕內(nèi)側(cè)水位逐漸降低;在帷幕靠近水源一側(cè)形成高水位區(qū),帷幕外側(cè)水位不斷抬升,且帷幕內(nèi)外側(cè)松散層水位差與截水帷幕完成率呈正相關(guān)關(guān)系;當(dāng)帷幕完成后,帷幕外側(cè)水位仍在抬升,內(nèi)側(cè)水位進(jìn)一步降低,帷幕內(nèi)外側(cè)水位差增大。
c.截水帷幕內(nèi)外側(cè)松散層水位解析計(jì)算和物理模擬可分析簡(jiǎn)單地質(zhì)條件的松散層水位變化特征,指導(dǎo)截水帷幕位置的選擇和松散層水位保護(hù)區(qū)域的確定,下一步將進(jìn)一步開展復(fù)雜地質(zhì)條件下松散層水位計(jì)算與模擬研究。