曾一凡,包 函,武 強,孟世豪,華照來,苗彥平,張 曄,卜文揚
(1.中國礦業(yè)大學(北京) 國家煤礦水害防治工程技術研究中心,北京 100083;2.中國礦業(yè)大學(北京) 內(nèi)蒙古研究院,內(nèi)蒙古 鄂爾多斯 017000;3.中國礦業(yè)大學(北京) 地球科學與測繪工程學院,北京 100083;4.礦山水害防治國家礦山安全監(jiān)察局重點實驗室,北京 100083;5.陜西陜煤曹家灘礦業(yè)有限公司,陜西 榆林 719000;6.陜煤集團神木紅柳林礦業(yè)有限公司,陜西 神木 719300)
煤炭是我國的基礎能源和重要工業(yè)原料,是國家繁榮及可持續(xù)發(fā)展的關鍵。據(jù)國家統(tǒng)計局數(shù)據(jù),2022 年,規(guī)模以上工業(yè)原煤產(chǎn)量45.0 億t,創(chuàng)歷史新高。陜西榆神礦區(qū)位于我國西部干旱半干旱地區(qū),生態(tài)環(huán)境脆弱,但其巨大的原煤產(chǎn)量在我國能源安全中占據(jù)了極其重要的位置。榆神礦區(qū)煤與水的賦存空間組合特點是“煤水共生,水在上、煤在下”[1],因此頂板水害威脅成為該區(qū)域煤礦高效生產(chǎn)和生態(tài)環(huán)境保護的主要矛盾。
近年來,榆神礦區(qū)煤炭資源的高強度開采,造成了地下水水位下降,引發(fā)了一系列生態(tài)問題[2-3]。新近系保德組紅土位于煤層上覆基巖與潛水含水層之間,其阻水能力直接關系到煤礦生產(chǎn)安全,更是工農(nóng)業(yè)、生活、生態(tài)用水安全的重要保障。基于榆神礦區(qū)的特殊地位,不少學者針對此地區(qū)的關鍵隔水層開展研究。在推動榆神礦區(qū)保水采煤工作中[4-6],李文平等[7]以榆神礦區(qū)地表的生態(tài)環(huán)境類型為基礎對各礦井展開保水采煤等級劃分。鄧念東等[8]從保水采煤角度將榆神礦區(qū)煤層開采的工程地質(zhì)條件劃分為5 個區(qū)域,并對各區(qū)域開采的難易程度進行分析。馬立強等[9]為解決隔水層厚度小、保水開采難度大的問題,提出一種壁式連采連充的采煤方法。王雙明等[10]點出保水采煤是實現(xiàn)綠色開采的手段之一,也是未來煤炭開采的趨勢。此外,丁仲禮等[11]通過對紅土的土壤學特征進行野外觀察,結(jié)合紅土的沉積學特征及稀土元素分析確認黃土高原地區(qū)紅土為風成成因。賀曉浪等[12]提出了保水安全厚度的概念,并針對陜北小保當煤礦建立了保水安全厚度的計算模型。大量物理相似模擬試驗[13-15]為研究采動作用下紅土隔水機理奠定了基礎。
以上研究對榆神礦區(qū)保水采煤工作起到了推動作用,對榆神礦區(qū)保水采煤及生態(tài)保護作出了巨大貢獻。前人對紅土性質(zhì)已做過大量研究,但是針對紅土正常沉積到紅土尖滅缺失的沉積過渡區(qū)(紅土薄弱區(qū))的研究很少。紅土薄弱區(qū)的紅土阻水性能判定僅以紅土厚度進行線性推斷,其工程性質(zhì)和表征其阻水性能的關鍵因素仍需深入挖掘,對于礦井生產(chǎn)和生態(tài)保護的有效隔水層厚度尚未形成統(tǒng)一的評判方法,臨界厚度目前仍難以準確界定。為緩解高強度采煤與脆弱生態(tài)環(huán)境之間的矛盾,同時也是保障煤炭的安全開采,在收集大量鉆孔數(shù)據(jù)的基礎上,通過在紅土沉積邊緣布設鉆孔和取樣對比,對榆神礦區(qū)保德組紅土的厚度分布、工程特性及阻水性能進行進一步的研究。首先通過巖性描述和測井數(shù)據(jù)重新優(yōu)化紅土厚度數(shù)據(jù),并結(jié)合沉積規(guī)律對紅土層分布進行修正;之后,從微觀尺度對比正常沉積區(qū)和沉積薄弱區(qū)紅土的工程特性。最后,通過計算非完全阻水的紅土薄弱區(qū)臨界厚度,進行整個榆神礦區(qū)紅土阻水性能的劃分。深入研究紅土的工程性質(zhì)和阻水能力將為推動保水工作提供重要的理論依據(jù)。
榆神礦區(qū)為本次研究的主要區(qū)域,同時將研究區(qū)范圍延伸至相鄰的含紅柳林、檸條塔、張家峁煤礦的神南礦區(qū)(圖1)。研究區(qū)位于陜西省榆林市西北部地區(qū),多年年均降水量約為411.2 mm,季節(jié)降水量分布不平衡,降水量由西北向東南遞增,屬于水資源匱乏型的干旱半干旱區(qū)域,是非常典型的生態(tài)脆弱區(qū)。榆神礦區(qū)地表多被新近系、第四系覆蓋,西部是由風力作用形成的沙丘地貌,而東部和北部都是以黃土梁峁為主。地表水系自東向西分屬黃河一級支流窟野河、禿尾河和黃河二級支流榆溪河流域(圖2)。礦區(qū)侏羅系中統(tǒng)延安組中賦存大量煤炭資源,其中主采煤層為2-2煤,且賦存較為穩(wěn)定,除部分區(qū)域發(fā)生自燃造成缺失,大部分區(qū)域可采。
圖1 榆神礦區(qū)地貌及水文特征分布Fig.1 Distribution of landform and hydrological characteristics in Yushen mining area
圖2 榆神礦區(qū)A—A'典型剖面Fig.2 A-A' typical profile in Yushen mining area
根據(jù)地下水的賦存條件和水力特征,可將主采煤層以上的含水層分為3 類[16]。最底層的基巖裂隙含水層主要由侏羅系延安組(J2y)、直羅組(J2z)和安定組(J2a)砂巖組成,結(jié)構(gòu)致密,含水量低;風化基巖孔隙裂隙含水層位于基巖頂部30 m 內(nèi),整個研究區(qū)均有分布,其特點是風化嚴重,滲透性強,孔隙大,為地下水儲存提供了充足的空間?;鶐r含水層的水量往往較少,風化基巖含水層的水量相對豐富,是煤礦開采的主要充水含水層。最后,第四系松散孔隙含水層主要由第四系薩拉烏蘇組(Q3s)和離石組(Q2l)組成,易于接受降雨補給,側(cè)向補給充分。關鍵隔水層是位于第四系含水層和風化基巖含水層之間的新近系保德組紅土(N2b)。
榆神礦區(qū)紅土分布不均勻,厚度及阻水性能有較大差異。本次研究共收集研究區(qū)鉆孔425 個(表1),其中為準確描述陜西曹家灘煤礦12 盤區(qū)西翼紅土層從沉積到尖滅的分布特征,布置了8 條勘探線,共22 個勘探鉆孔,其中紅土層全段取心。本次研究鉆孔數(shù)據(jù)的收集區(qū)域超過5 000 km2,統(tǒng)計時間跨越50 a。但紅土層在不同時期、不同位置的鉆孔描述中并不完全一致,因此通過巖性描述和測井數(shù)據(jù)重新優(yōu)化紅土數(shù)據(jù),并結(jié)合沉積規(guī)律對第四系松散層及紅土層數(shù)據(jù)進行修正。
表1 鉆孔分布統(tǒng)計結(jié)果Table 1 List of borehole distribution
修正后的紅土沉積厚度分布如圖3a 所示,紅土厚度最大可達140 m 以上,但分布極不均勻,存在大面積缺失區(qū),缺失面積達2 597.23 km2,占比達48.80%。在孟家灣、爾林灘、阿包兔、中雞、馬王廟勘察區(qū),金雞灘、杭來灣煤礦等均有較大范圍的缺失。而在曹家灘煤礦附近,紅土厚度達到最大,僅在礦井西南側(cè)存在局部缺失。
圖3 榆神礦區(qū)新近系保德組紅土分布特征Fig.3 Distribution characteristics of Neogene Baode Formation laterite in Yushen mining area
通過對比圖3a、圖3b 可知,紅土沉積厚度與沉積底界高程具有較高吻合性,總體上高程越低,紅土沉積厚度越大。在5.2 Ma B.P.之前,榆神礦區(qū)地勢總體上為東高西低,伴隨著西部六盤山的隆起,發(fā)生了一次反轉(zhuǎn),導致東凹西凸,造成了新近系紅土沉積時西高東低的沉積底界(圖3b)。之后在風力作用下,周邊山區(qū)的巖屑和泥沙輸送到盆地內(nèi),并在干旱環(huán)境下發(fā)生氧化而形成了紅土沉積物[17]。故其沉積厚度呈現(xiàn)由西北向東南逐漸增大的總體趨勢。
東南部低洼的環(huán)境給紅土沉積提供了極佳條件,但隨著河流的演化,在河流侵蝕作用下紅土的分布被進一步打亂,最終形成目前極不規(guī)則的分布特征,更是造成榆神礦區(qū)西側(cè)有大面積紅土缺失區(qū)的直接原因。通過圖3 對比可知,河流侵蝕和沉積底界的高程差異是造成紅土厚度分布不均勻的主要原因。
研究區(qū)內(nèi)紅土顏色呈淺紅棕色。紅土碎屑以石英為主,其次是長石,含少量巖屑,填隙物主要為黏土礦物;碎屑顆粒分選性較好,磨圓度差,呈次棱角狀,基底式膠結(jié),內(nèi)含較多鈣質(zhì)結(jié)核,且隨埋深的增加,結(jié)核密集程度也隨之增加。研究區(qū)內(nèi)紅土鈣質(zhì)結(jié)核多呈星散狀展布,極少成層。紅土層經(jīng)過淋濾,土體中的鈣質(zhì)成分下滲,在原有顆粒周圍或空隙中沉積形成鈣質(zhì)結(jié)核。保德組紅土層沉積相主要為陸相風成或半風成沉積,反映了當時的干旱氣候和強烈的風力作用。通過沉積相分析,紅土層垂向可分為兩段(圖4),上部亞相主要由黃褐色、紅棕色亞砂土及亞黏土組成,這些沉積物是由于風力將周邊山區(qū)的巖屑和泥沙輸送到盆地內(nèi),并在干旱環(huán)境下發(fā)生氧化而形成;下部亞相主要由土紅色黏土組成,富含鈣質(zhì)結(jié)核。這些沉積物可能是由于風力將遠處的火山灰輸送到盆地內(nèi),并在濕潤環(huán)境下發(fā)生水化和碳酸鹽化而形成。為進一步研究其阻水性能,則需要通過對紅土樣品的微觀特征進行測試,開展垂向和水平對比,對其沉積環(huán)境及隔水性能展開分析。
圖4 紅土土樣及內(nèi)部結(jié)核Fig.4 Laterite samples and internal nodules
本節(jié)通過選取同一位置不同埋深的紅土樣品進行粒徑特征和滲透系數(shù)的垂向?qū)Ρ?,從而進一步論述紅土層的沉積過程及隔水性能在垂直方向的差異。
3.1.1粒徑組成
沉積物的粒徑組成可以反映沉積物搬運介質(zhì)類型、搬運方式及沉積環(huán)境等特征,是重要的環(huán)境分析指標,并在紅土環(huán)境分析中受到廣泛應用[18]。榆神礦區(qū)紅土以粉土顆粒和黏土顆粒為主,表現(xiàn)出風成的特點。通過一系列上下段紅土樣品的對比可以發(fā)現(xiàn)(圖5),兩段紅土質(zhì)量分數(shù)最大的均為粉土顆粒,占比分別為50.53%和46.24%。此外相比于上段沉積土樣,下段土樣的小粒徑黏土顆粒相對較多,質(zhì)量分數(shù)達到35.55%??傮w上看,紅土顆粒的集中程度相對較高,分選性較好,上段粒徑整體大于下段粒徑,下段的隔水性能更好。粒徑特征符合風力搬運的沉積描述,只是氣候略有不同,上段搬運營力相對下段略有增強。
圖5 紅土沉積上下段樣品的粒徑特征對比Fig.5 Comparison of particle size characteristics of upper and lower samples of laterite deposition
3.1.2滲透性特征
為直觀掌握紅土隔水能力的強弱,通過自主研發(fā)的新型土體高壓滲流試驗系統(tǒng),開展土樣滲流特性試驗,研究土體在原位恢復狀態(tài)下的滲透系數(shù)變化規(guī)律。首先對土樣施加軸壓和圍壓恢復原位應力,在此基礎上施加上下滲流反壓,測試土樣在不同滲流壓力下滲透水量。由于深層土中充填結(jié)合水,不同于重力水,存在起始梯度問題,所以通過測定滲透水量與時間的變化曲線,計算深層黏性土的起始梯度及符合達西定律的時間,進而得到滲透系數(shù)。
滲透系數(shù)的計算結(jié)果如圖6 所示,通過對比上下段樣品,可見相同位置的樣品隨著沉積深度的增加,滲透系數(shù)大幅減小。說明除了沉積特征以外,下段紅土所受到的壓實作用時間更長,其結(jié)構(gòu)更為緊密,隔水性能更強。此外,在水平方向,不同沉積厚度紅土的滲透系數(shù)也具有一定的規(guī)律性。
圖6 紅土沉積上下段樣品的滲透系數(shù)對比Fig.6 Comparison of permeability coefficients of upper and lower samples of laterite deposits
通過選取不同紅土沉積厚度的下段土樣進行水平對比,通過SEM 掃描電鏡微觀結(jié)構(gòu)和X 射線衍射XRD物質(zhì)組成來進一步分析紅土層的隔水性能。
3.2.1SEM 微觀結(jié)構(gòu)特征
SEM 電子顯微鏡的基本原理是利用高能電子束來照射紅土結(jié)構(gòu)表面,電子束會與紅土樣品物理作用產(chǎn)生不同的檢測信號,再由探測器接收,便可得到紅土試樣的微觀結(jié)構(gòu)信息。
結(jié)果如圖7 所示,在3 萬倍鏡下清晰可見有大量纖維狀交織在一起,其孔隙極小,說明紅土具有極佳的隔水性。通過對比多組樣品及單個樣品的多組SEM影像圖片可以發(fā)現(xiàn),位于正常沉積區(qū)的紅土樣品交織結(jié)構(gòu)緊密,以纖維狀結(jié)構(gòu)(圖7a)為主,而紅土薄弱區(qū)的紅土樣品少見交織結(jié)構(gòu),呈現(xiàn)出零散的片狀結(jié)構(gòu)(圖7b)。
圖7 紅土樣品的微觀結(jié)構(gòu)特征對比Fig.7 Comparison of microstructure characteristics of samples with different laterite deposition thickness
3.2.2地球化學特征
紅土的形成受多種因素控制,分析紅土黏土礦物組成對其形成環(huán)境及隔水能力有重要的作用[19]。紅土中含有較多的黏土礦物,黏土礦物具有遇水膨脹的特性,會形成強大的黏聚力,增強留水能力,大大降低滲透率,且不同的黏土礦物膨脹性不同,因此黏土的成分是評價含、隔水層的重要指標。本文選取不同沉積厚度的紅土樣品,分析不同厚度區(qū)域其成分含量。
X 射線衍射儀是用于分析黏土礦物成分的主要手段,具備方便快捷、精度高的優(yōu)點。依據(jù)設定角度、掃描時間等測試條件下對研究區(qū)紅土樣品進行測試,得到相關的X 衍射圖譜,利用物相檢索,定量計算礦物成分組成。
在樣品礦物組成定性分析的基礎上,采用K值計算法對礦物含量進行相對定量計算[20]。其原理是物相衍射強度和含量呈正相關關系,可通過衍射強度的大小求出混合物中某種物相的衍射的體積分數(shù)或質(zhì)量分數(shù)。即通過不同礦物衍射峰強度值與K值相對比,得到不同礦物在樣品中的相對含量,計算公式如下:
式中:i為需要計算含量的礦物編號;j為參比的礦物編號;n為參比礦物總數(shù);Xi為礦物i的相對含量;I為礦物的衍射強度;K為PDF 卡片中礦物的參比強度。
X 衍射結(jié)果如圖8 所示,研究區(qū)紅土主要由石英、長石、方解石等非黏土礦物及包括伊利石、高嶺石、綠泥石、伊蒙混層在內(nèi)的黏土礦物組成。礦物成分的質(zhì)量分數(shù)見表2,紅土中的黏土含量與紅土沉積厚度呈正相關。紅土厚度為2.92 m 的區(qū)域黏土質(zhì)量分數(shù)僅為12.67%,隨著紅土厚度增大至31.8 m,其黏土質(zhì)量分數(shù)可達23.10%??梢姡S著紅土厚度的增加,起著主要隔水作用的黏土含量會顯著增加。
表2 新近系保德組紅土礦物成分組成Table 2 Mineral composition of N2b laterite
圖8 XRD 定性分析衍射圖譜結(jié)果Fig.8 XRD qualitative analysis of diffraction pattern results
進一步細分黏土中的礦物組成有助于明確紅土的沉積環(huán)境和阻水性能。比如蒙脫石在堿性介質(zhì)條件、溫帶半干旱半濕潤的環(huán)境下更易形成,親水性強,具有較強的膨脹性,從而具備優(yōu)良的隔水性能;伊利石和綠泥石在弱堿性環(huán)境、干旱氣候條件下能較多地保存,是較弱化學風化階段的產(chǎn)物;高嶺石是在溫暖濕潤氣候、酸性介質(zhì)條件下由長石、云母和輝石等礦物經(jīng)較強的化學風化過程形成。黏土礦物隔水能力大小的排序為蒙脫石>伊蒙混層>伊利石>高嶺石[21]。因此,不同的沉積環(huán)境會造成紅土的隔水性能不均一。
從圖9 可知,隨著紅土厚度的增大,蒙脫石、伊蒙混層這類強親水性礦物的質(zhì)量分數(shù)分別從31%和16%上升至37%和27%;伊利石、高嶺土、綠泥石這類相對弱親水性礦物的含量均呈減少趨勢,最大減幅達35.02%。通過分析可知,紅土厚的區(qū)域不僅其黏土礦物含量高,其中親水性礦物的含量也越高,隔水性能隨紅土厚度呈非線性關系。
圖9 不同紅土沉積厚度樣品中各黏土礦物成分對比Fig.9 Comparison of clay mineral composition in samples with different laterite deposition thickness
通過上述研究可以發(fā)現(xiàn),在同一沉積位置不同沉積深度的紅土阻水性能并不一致,同時位于沉積邊緣的紅土阻水能力遠弱于正常沉積紅土。僅用紅土厚度來表征該關鍵隔水層的阻水性能不夠嚴謹,本章旨在通過上述研究成果,進一步量化紅土阻水性能,基于穩(wěn)定達西流速概念進行紅土阻水臨界厚度的計算,并通過多源數(shù)據(jù)融合準確表征紅土薄弱區(qū)的紅土阻水性能,最終得到榆神礦區(qū)的紅土阻水性能分區(qū)。
首先通過上述粒徑組成研究成果,結(jié)合鉆孔測井的連續(xù)曲線確定紅土上下分段的閾值。并結(jié)合對應的數(shù)學模型計算紅土全段的等效阻水性能,為同一沉積位置不同沉積深度阻水能力不一致的問題提供有效解決方案。
(1) 基于測井曲線的上下分段。自然伽馬測井技術通過測量地層中的伽馬射線強度,可以反映地層中放射性元素的含量。自然伽馬測井可以為評價土層性質(zhì)提供重要參考信息,本次研究擬用測井數(shù)據(jù)來區(qū)分同一地層中沉積環(huán)境和古地理信息的差異。因為自然伽馬測量結(jié)果受地殼中放射性核素分布的影響,而這些放射性核素分布的空間變化與土層的密度、含水量等參數(shù)有關。一般來說,對于相同類型的土層,如果自然伽馬值越大,則說明土層密度越小,固結(jié)程度越低;反之,自然伽馬值越小,則表示土層密度越大,固結(jié)程度越高。通過單獨提取多個鉆孔紅土全段的自然伽馬曲線,進行對比分析可以發(fā)現(xiàn)(圖10),紅土層明顯分為上下兩段,上段的自然伽馬平均值大于下段,表明了紅土下段沉積固結(jié)性更好。通過對比大量鉆孔的測井曲線,上段厚度占全段比值在58%~65%,結(jié)合土樣顆粒組成的研究成果,最終設置紅土層沉積上下段比值為6∶4。
圖10 測井曲線分段對比Fig.10 Logging curve section comparison
式中:(i,j,k)為 單元空間位置的三維坐標;Qi,j,k、Qi,j,k+1分別為通過土層上段單元和下段單元的流量,L3/T;Q為通過整個土層的總流量,L3/T;Δh為水頭差,L;A為單元面積,L2;Kv為垂向滲透系數(shù),L/T;M為單元厚度,L。
(2) 基于達西定律的等效滲透性計算模型。由于上下段紅土層所表現(xiàn)出的性質(zhì)具有差異,需要通過建立有效的計算模型得到能夠完整表現(xiàn)紅土層的等效滲透系數(shù)。計算模型如圖11 所示,將整個紅土層分為上下2 個單元,通過達西定律可以獲得分別通過上下單元的流量和通過整個土層的總流量。通過上下2 個單元的流量沒有累積或消耗,同時總水頭差由2 個單元水頭差疊加,代入上下兩段沉積厚度的閾值,并聯(lián)合以下公式進行推導,可以得到等效滲透系數(shù)的計算公式。
圖11 等效滲透系數(shù)數(shù)學概念模型Fig.11 Mathematical conceptual model of equivalent permeability coefficient
等效滲透系數(shù)的計算結(jié)果如圖12 所示,通過Logistic 模型函數(shù)擬合可以得到紅土等效滲透系數(shù)與紅土沉積厚度的關系。
圖12 紅土等效滲透系數(shù)與紅土沉積厚度的關系曲線Fig.12 Equivalent permeability coefficient vs deposition thickness of laterite
通過等效滲透系數(shù)和礦物組成對紅土阻水過程進行詳細描述,最終在榆神礦區(qū)紅土厚度分布基礎上,形成紅土阻水性能的準確劃分。
黏土的滲透分為2 個過程,第一階段可以簡單概括為非飽和黏土的滲流,由于黏土的親水性和膨脹性,低水力梯度下,水流實際過水斷面小、流動途徑彎曲程度大、孔隙通道阻力大等多種原因使得難以形成有效流動;在水力梯度增大,或充分飽水情況下,第二階段可概括為穩(wěn)定的飽和達西流。當在一定水力梯度下,水流僅停留在第一階段時對應的紅土層厚度可成為有效阻水臨界厚度。紅土臨界厚度以內(nèi)會產(chǎn)生有效流動,水力梯度影響下會激發(fā)或襲奪上覆含水層形成越流補給。
通過對榆神礦區(qū)紅土層上下含水層的水位統(tǒng)計,以及采動影響下的水位變化,最大水位差約為45 m。以達到穩(wěn)定達西流速為臨界狀態(tài),可計算紅土有效阻水的臨界厚度為10 m。由此可將榆神礦區(qū)紅土分為完全阻水區(qū)和非完全阻水區(qū)(紅土薄弱區(qū))。在非完全阻水區(qū)中,通過礦物組成含量和等效滲透系數(shù)對紅土阻水性能進行準確劃分,為位于沉積邊緣的紅土阻水能力遠弱于正常沉積紅土且無法準確描述紅土阻水能力與厚度非線性關系的問題提供有效解決方案。最終結(jié)果如圖13 所示,將榆神礦區(qū)紅土阻水性能分為3 個大區(qū),分別為完全阻水區(qū)、紅土薄弱區(qū)和紅土缺失區(qū)。進一步將紅土薄弱區(qū)分為極弱阻水區(qū)、阻水過渡區(qū)和高阻水區(qū)。
圖13 榆神礦區(qū)紅土阻水性能分區(qū)Fig.13 The zoning of water resistance performance of laterite in Yushen mining area
曹家灘煤礦位于榆神礦區(qū)I 期中部(圖13),其工作面西翼位于紅土沉積邊緣,導致東西兩翼截然不同的涌水特征[22]。通過水樣測試和數(shù)據(jù)收集,分析各含水層的水化學特征[23]。進一步將水化學數(shù)據(jù)和紅土分布厚度聯(lián)系可以發(fā)現(xiàn)(圖14),紅土上部的第四系薩拉烏蘇組間接充水含水層中Ca2+毫克當量百分比在全礦區(qū)分布均勻,平均為62.25%;而在紅土下部的風化基巖直接充水含水層中分布不均,在礦區(qū)西部邊界的紅土缺失及薄弱區(qū),Ca2+毫克當量百分比比東區(qū)高近12%。在隔水層厚度為0~10 m 的風化基巖地下水中Ca2+毫克當量百分比為43.33%,而當紅土隔水層厚度大于10 m 時急劇降低,平均為27.62%。由此說明當紅土厚度為0~10 m 時為非完全阻水區(qū),位于紅土薄弱區(qū)的礦區(qū)需要特別關注煤層開采條件下淺表水體是否參與礦井涌水。
圖14 不同紅土厚度下Ca2+毫克當量對比Fig.14 Comparison of Ca2+ content underdifferent laterite thickness
作為典型的地下工程,礦山資源開發(fā)常常維持幾十年。早期以開采效率為主要目的,地下水主要體現(xiàn)的是災害屬性,對礦區(qū)地下水研究的主要目的是保障礦山安全。隨著新形勢下“以水定產(chǎn)”“把水資源作為最大的剛性約束”等政策的轉(zhuǎn)變,以及開采重心轉(zhuǎn)向生態(tài)更加脆弱的西北地區(qū),礦區(qū)地下水的資源屬性逐漸被重視起來。其發(fā)展主要分為以下幾個階段:(1) 礦井水利用構(gòu)想階段:自20 世紀80 年代,國內(nèi)外學者開始對地下資源開發(fā)引起的水資源破壞問題展開研究,Wu Qiang 等[24]宏觀上提出了排水、供水和環(huán)境保護三位一體的礦井涌水處理優(yōu)化方案,針對礦井水的研究從單純的治理轉(zhuǎn)變?yōu)樗廴咎幚砑熬C合利用。(2) 狹義的保水開采階段:該階段是通過研究覆巖及含水層的破壞程度,轉(zhuǎn)變開采方式和方法,調(diào)整煤柱間距和尺寸,進行采空區(qū)充填等技術[25],形成一種具有被動性和局限性的保水開采模式[26]。(3) 廣義的保水開采階段:為了增大開采效率,最大程度保證資源采出率,一方面通過利用或再造關鍵隔水層保護水資源安全,另一方面通過轉(zhuǎn)移地下水至地下采空區(qū),構(gòu)成地下水庫來維持水資源的開發(fā)利用價值。(4) 煤水雙資源協(xié)同開采階段:通過水害防控消除地下水災害屬性的負效應,基于礦井水資源化利用充分挖掘其資源屬性的正效應[27]。以全局視角協(xié)調(diào)資源開發(fā)過程,實現(xiàn)礦井水控制、處理、利用、回灌與生態(tài)保護“五位一體”的優(yōu)化管理[28]。本文對紅土關鍵隔水層的穩(wěn)定性和阻水性能的研究,對于榆神礦區(qū),特別是非完全阻水的紅土薄弱區(qū)礦井實現(xiàn)煤水雙資源協(xié)同共采具有重要意義。
首先,關鍵隔水層工程性質(zhì)和阻水性能的深入分析,是查清紅土薄弱區(qū)礦井充水因素和隱蔽致災地質(zhì)因素的關鍵步驟。通過本文研究可以進一步綜合隔水層的滲透性能、破壞概率、穩(wěn)定范圍等指標,有助于查清充水通道水力性質(zhì),并進行水害風險辨識。
其次,關鍵隔水層工程性質(zhì)和阻水性能的系統(tǒng)評估,是實現(xiàn)水害分源精準防治和水資源綜合利用的必要條件。通過本文研究可以針對紅土薄弱區(qū)推演其應力狀態(tài)、變形特征、裂縫發(fā)育等情況,從而精準識別水害源頭及災變模式。通過選擇最適宜的加固材料、方法和時機,提高隔水層保護措施的可行性和有效性,實現(xiàn)松散層淺表水的防護和治理。
最后,關鍵隔水層研究對于指導保水開采工程實踐有著重要的價值。本研究成果可為煤水雙資源協(xié)同開采提供更準確和科學的評價依據(jù),從而確定最適合的開采方法、參數(shù)和條件,提高開采效率和安全性。根據(jù)紅土隔水層在不同地區(qū)、不同條件下表現(xiàn)出來的工程性質(zhì)和阻水性能差異,分析其影響因素和演化規(guī)律,并給出相應的調(diào)整和改進建議,從而提高煤水雙資源協(xié)同開采在實際應用中的適應性和靈活性,促進該技術在實際應用中的推廣和發(fā)展。
a.通過巖性描述和測井數(shù)據(jù)重新優(yōu)化紅土數(shù)據(jù),榆神礦區(qū)紅土缺失區(qū)面積達48.80%。結(jié)合紅土沉積底界高程,認為河流侵蝕和沉積底界的高程差異是造成紅土厚度分布極不均勻的直接原因。
b.紅土的隔水性能與其厚度之間并非呈線性關系,隨著紅土厚度的減小,其隔水性能急劇減小。通過多種手段的綜合判定,給定紅土層上下2 段的閾值為6∶4;并在厚度分布的基礎上通過等效滲透系數(shù)和礦物組成對紅土阻水性能進行分區(qū)。
c.對紅土關鍵隔水層的穩(wěn)定性和阻水性能的研究,對于榆神礦區(qū),特別是非完全阻水的紅土薄弱區(qū)礦井查清礦井充水因素和隱蔽致災地質(zhì)因素、關鍵隔水層保護及再建和實現(xiàn)煤水雙資源協(xié)同共采提供了重要的理論基礎。