趙春虎，虎維岳，靳德武

(1.中國煤炭科工集團(tuán)西安研究院有限公司，陜西 西安 710054；2.陜西省煤礦水害防治技術(shù)重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710177)

西部干旱礦區(qū)采煤引起潛水損失量的定量評(píng)價(jià)方法

趙春虎1,2，虎維岳1,2，靳德武1,2

(1.中國煤炭科工集團(tuán)西安研究院有限公司，陜西 西安 710054；2.陜西省煤礦水害防治技術(shù)重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710177)

針對(duì)西部干旱礦區(qū)煤炭資源開發(fā)對(duì)地下水?dāng)_動(dòng)的相關(guān)研究定量化程度不高問題，通過地下水動(dòng)力學(xué)分析方法，得出采掘擾動(dòng)形成的覆巖垮落帶、裂隙帶、彎曲帶及地面沉陷區(qū)的空間尺度是控制地下水損失的主要因素，其擾動(dòng)面積越大、冒裂帶高度越高、彎曲帶保護(hù)層滲透能力越強(qiáng)則地下水滲漏量越大；提出了覆巖導(dǎo)水裂縫邊界化，彎曲帶滲透能力分區(qū)，地面沉陷區(qū)重新剖分的采掘擾動(dòng)影響下地下水系統(tǒng)數(shù)值評(píng)價(jià)模型的構(gòu)建方法；并通過構(gòu)建案例礦井的評(píng)價(jià)模型，計(jì)算得出神東礦區(qū)補(bǔ)連塔煤礦由于煤礦開采引起松散層潛水損失量為1.90×104m3/d。

地下水；松散含水層；定量評(píng)價(jià)；數(shù)值模擬；導(dǎo)水裂縫帶

地下水是西部干旱地區(qū)水資源主要組成部分，具有極其重要的資源功能和生態(tài)功能。而隨著煤炭資源開發(fā)西進(jìn)戰(zhàn)略的實(shí)施，我國煤炭基地分布與水資源豐富程度已呈逆向分布。尤其西部干旱礦區(qū)在規(guī)?；?、現(xiàn)代化的煤炭資源開發(fā)過程中不可避免地造成了對(duì)地下水影響和破壞，加劇了地區(qū)資源性缺水與生態(tài)環(huán)境破壞，是西部煤炭工業(yè)可持續(xù)發(fā)展面臨的主要問題之一。由于能源基礎(chǔ)結(jié)構(gòu)的不同，國外對(duì)地下水水動(dòng)力擾動(dòng)問題的研究極少，主要以礦區(qū)水污染風(fēng)險(xiǎn)評(píng)價(jià)、關(guān)閉礦井修復(fù)、復(fù)墾以及相關(guān)法規(guī)制定等為主要研究方向[1-4]。國內(nèi)眾多學(xué)者采用不同技術(shù)方法針對(duì)西部地區(qū)煤炭開采造成的地下水位下降、地下水水量損失等水動(dòng)力響應(yīng)問題進(jìn)行了研究[5-6]，但是尚未形成公認(rèn)的評(píng)判方法，直到2011年2月，環(huán)保部第一次正式頒布了《環(huán)境影響評(píng)價(jià)技術(shù)導(dǎo)則-地下水環(huán)境》，對(duì)于指導(dǎo)我國開展地下水資源與環(huán)境保護(hù)工作具有重要的現(xiàn)實(shí)意義。

目前水文地質(zhì)界常用大井法、降落漏斗法[7-8]，以及煤炭生產(chǎn)管理部門的噸煤排水量法等均以宏觀分析為主，評(píng)價(jià)的定量程度不足。近年來數(shù)值法在水資源評(píng)價(jià)、預(yù)測和管理中得到迅速的發(fā)展，使地下水資源研究從傳統(tǒng)的研究方法轉(zhuǎn)到模型研究，大大提高了地下水資源評(píng)價(jià)定量化程度。如劉懷忠[9]、張偉[10]、李瑩[11]、白樂[12]、吳喜軍[13]等基于數(shù)值模擬手段將礦井涌水量或圍巖滲透性乘以一定的比例系數(shù)(如1.2,1.5倍等)，通過變化該類比例系數(shù)來預(yù)測地下水動(dòng)態(tài)變化，達(dá)到定量評(píng)價(jià)采煤對(duì)地下水影響程度的目的，其中涌水量系數(shù)和圍巖滲透系數(shù)的給定具有較強(qiáng)的主觀性，忽視了采煤導(dǎo)致地下水響應(yīng)機(jī)制問題。本文將采礦覆巖破壞與地下水系統(tǒng)模擬研究結(jié)合，通過分析采煤對(duì)地下水?dāng)_動(dòng)的控制因素，探索數(shù)值化的處理方法，建立定量化程度高的地下水?dāng)?shù)值評(píng)價(jià)模型，為我國西部干旱礦區(qū)煤-水資源協(xié)調(diào)開發(fā)提供科學(xué)依據(jù)。

1 地下水環(huán)境擾動(dòng)機(jī)理

1.1 采動(dòng)覆巖的破壞規(guī)律

采煤活動(dòng)導(dǎo)致圍巖體含、導(dǎo)水能力的變異是地下水水動(dòng)力變化的根本原因。目前，采礦與地下水研究的相關(guān)學(xué)者、技術(shù)人員對(duì)采動(dòng)覆巖破壞有較為一致的認(rèn)識(shí)，以經(jīng)典的采礦“覆巖分帶”[14]理論為研究基礎(chǔ)，認(rèn)為煤層開采后采動(dòng)覆巖由下至上依次分為垮落帶、裂隙帶和彎曲帶，同時(shí)在地表形成地面沉陷區(qū)。因而，典型采動(dòng)覆巖擾動(dòng)(圖1)可以總結(jié)為“三帶一區(qū)”的變形損傷規(guī)律。

圖1 采動(dòng)覆巖影響示意

1.2 對(duì)地下水影響機(jī)制分析

陜北與神東煤炭基地位于我國西北干旱地區(qū)，松散層地下潛水是區(qū)域惟一具有供水意義和重要生態(tài)價(jià)值的水資源[15]。而采礦擾動(dòng)使賦存地下水的巖層發(fā)生 “三帶一區(qū)”損傷變形，由此產(chǎn)生的松散層地下水向采掘空間的滲漏流失。因此，須以采礦覆巖破壞為研究基礎(chǔ)，分析采礦引起松散含水層地下水的滲漏損失問題。

煤層開采導(dǎo)致的地下水損失(滲漏)量主要是由松散層地下水從垂向上向采掘空間的滲漏量Q1、冒裂帶周邊揭露含水層側(cè)向排泄量Q2和地面沉陷積水引起無效蒸發(fā)量Q3三部分組成。

(1)地下水滲漏量(Q1)

深埋煤層開采冒裂帶一般未發(fā)育至松散含水層底部，當(dāng)余留一定的基巖保護(hù)層厚度時(shí)，為了便于分析擾動(dòng)因素對(duì)地下水滲漏的影響，認(rèn)為采掘擾動(dòng)形成的導(dǎo)水裂縫為一倒梯形，如圖2所示，定義覆巖體頂部截面積為A(m2)；冒裂帶(導(dǎo)水裂縫帶)發(fā)育高度為L1(m)；松散孔隙含水層與煤層之間基巖厚度為M(m)；其滲透系數(shù)為K(m/d)；松散含水層與保護(hù)層界面地下水水壓力為P(MPa)；u為保護(hù)層中地下水的實(shí)際流速；r為地下水容重(N/m3)。

圖2 地下水環(huán)境擾動(dòng)示意

定義煤層頂板為“0”基準(zhǔn)面，根據(jù)經(jīng)典地下水動(dòng)力學(xué)分析方法，得導(dǎo)水裂縫帶頂部(L1高度處)的水頭值為

H1=L1+P/r+u2/2g

由于井下排水，冒裂帶內(nèi)直接與大氣聯(lián)通，孔隙壓力為大氣壓，即P≈0，則

H1=L1+u2/2g

保護(hù)層頂部(L2高度處)的水頭值為

H2=L2+P/r+u2/2g

則松散含水層向采掘空間滲漏量Q1為

由于L2-L1=M-L1，則

根據(jù)式(2)可以看出，滲漏量Q1與采動(dòng)導(dǎo)水裂縫帶高度L1、彎曲帶滲透能力K、導(dǎo)水裂縫帶揭露范圍A以及含水層水壓P呈正相關(guān)關(guān)系，與余留的保護(hù)層厚度(M-L1)呈負(fù)相關(guān)。即擾動(dòng)揭露的面積越大、冒裂帶高度越高、彎曲帶保護(hù)層滲透能力越強(qiáng)則地下水滲漏量越大。

(2)側(cè)向排泄量(Q2)

淺埋煤層開采冒裂帶一般直接發(fā)育至松散含水層內(nèi)或地表，即保護(hù)層厚度為0，采掘空間頂部含水層地下水直接疏干。地下水資源流失以側(cè)向排泄為主，根據(jù)基于穩(wěn)定流分析的“大井法”，含水層地下水向采掘空間的側(cè)向排量泄Q2為

式中，K為滲透系數(shù)，m/d；H為水頭高度，m；S為由于礦井排水而產(chǎn)生的水位降深值，m；R0為引用影響半徑，m；r0為假想“大井”的半徑(即為導(dǎo)水裂縫揭露區(qū)域的引用半徑)，m。

(3)地下水無效蒸發(fā)量(Q3)

深埋煤層開采一般在地面形成沉陷區(qū)，當(dāng)松散層地下水埋深較淺時(shí)可形成地面積水，地下水由潛水蒸發(fā)轉(zhuǎn)化成蒸發(fā)能力極強(qiáng)的水面蒸發(fā)，則產(chǎn)生了地下水無效蒸發(fā)量(Q3)，定義松散層潛水的蒸發(fā)強(qiáng)度為E1，水面蒸發(fā)為E2，假設(shè)地面已形成地面沉陷積水，積水面積為A2(m2)，則

根據(jù)式(4)可以看出，無效蒸發(fā)量Q3與地區(qū)水面蒸發(fā)能力E2和積水面積A2正相關(guān)。從采掘擾動(dòng)角度分析，積水面積與采掘煤層具體地質(zhì)條件、含水層水文地質(zhì)條件相關(guān)，即采煤厚度越大、面積越大、水位埋深越小則更易形成地面積水。

綜上以上3種形式松散層地下水損失情況，導(dǎo)水裂縫帶影響的高度、平面范圍、彎曲帶保護(hù)層厚度及滲透能力、地面沉陷范圍等是控制地下水資源漏失的主要因素，即采掘擾動(dòng)形成的“三帶一區(qū)”的空間尺度是地下水?dāng)_動(dòng)程度的控制因素。

2 采煤對(duì)地下水?dāng)_動(dòng)的評(píng)價(jià)模型構(gòu)建

由前面分析可知，采掘擾動(dòng)形成的導(dǎo)水裂縫帶、彎曲帶、地面沉陷區(qū)的空間尺度是地下水損失的控制要素，下面就“三帶一區(qū)”在地下水系統(tǒng)數(shù)值模型數(shù)值化處理技術(shù)展開研究。

(1)導(dǎo)水裂縫帶數(shù)值化處理技術(shù)

從地下水系統(tǒng)角度分析，當(dāng)導(dǎo)水裂縫發(fā)育到某一高度(層位)，在冒裂帶發(fā)育范圍內(nèi)的含水層會(huì)被直接溝通，采空區(qū)頂部含水層地下水體沿垂向采動(dòng)裂縫迅速涌入采空區(qū)，在無垂向無越流或大氣降水直接補(bǔ)給時(shí)頂部含水層被迅速疏干。

如圖3所示，在導(dǎo)水裂縫帶外圍邊界處，形成了切割含水層的地下水排泄“帶”，地下水沿排泄“帶”涌入采空區(qū)形成了較為穩(wěn)定的礦井正常涌水，造成地下水水位在排泄“帶”處下降至該含水層底板，其水位高度等于該處的含水層底板標(biāo)高。

圖3 地下水與導(dǎo)水裂縫帶關(guān)系示意

在導(dǎo)水裂縫帶上方外圍邊界處，當(dāng)導(dǎo)水裂縫未揭露上覆的含水層，上部的含水層在垂向上可能形成穩(wěn)定的越流補(bǔ)給，由于導(dǎo)水裂縫直接與大氣聯(lián)通，其孔隙水壓力為大氣壓力，即水頭高度Hc近似等于該處的位置水頭標(biāo)高(冒裂帶高程Hm)為

Hc=Hm

從水文地質(zhì)條件概化的角度，導(dǎo)水裂縫切割含水層的接觸帶構(gòu)成了地下水系統(tǒng)的“內(nèi)邊界”，可以概化為地下水運(yùn)動(dòng)的一類水頭邊界條件，且水頭邊界水頭值即為該處含水層底板標(biāo)高。以此為約束條件，在地下水系統(tǒng)數(shù)值模型中將導(dǎo)水裂縫切割含水層的接觸帶數(shù)值化處理成地下水系統(tǒng)的一類內(nèi)部邊界，通過識(shí)別含水層底板標(biāo)高，來定義水頭邊界的水位值，以計(jì)算模擬含水層地下水的變化特征。

(2)“彎曲帶”覆巖滲透能力分析與數(shù)值化處理技術(shù)

淺埋煤層開采導(dǎo)水裂縫一般直接發(fā)育至地表，覆巖為垮落帶與裂隙帶“兩帶型”采動(dòng)裂縫，深埋煤層開采易形成典型的“三帶一區(qū)型”結(jié)構(gòu)，彎曲帶雖未形成具有較強(qiáng)導(dǎo)水能力的擾動(dòng)裂縫，地應(yīng)力狀態(tài)變化使覆巖發(fā)生不同程度與形式變形。如圖4所示，位于采掘空間上方彎曲帶的壓縮區(qū)其巖體滲透能力相對(duì)減小，兩側(cè)拉伸區(qū)巖體滲透能力相對(duì)減小。根據(jù)西部礦區(qū)相關(guān)研究成果，文獻(xiàn)[16]在神東礦區(qū)利用地球物理手段探測出采前采后彎曲帶覆巖含水性整體較為穩(wěn)定，文獻(xiàn)[17-18]利用原位鉆孔壓水試驗(yàn)得出位于采動(dòng)拉伸區(qū)的離石組黃土滲透能力較采前增大了1～2個(gè)數(shù)量級(jí)，而同位于壓縮區(qū)與拉伸區(qū)的保德組紅土，其采前采后滲透能力相當(dāng)?shù)慕Y(jié)論?？梢姀澢鷰B透能力變異對(duì)地下水漏失影響小。

圖4 彎曲帶滲透能力分區(qū)

在地下水系統(tǒng)建模時(shí)，由于采掘擾動(dòng)后彎曲帶覆巖仍為連續(xù)介質(zhì)，可以根據(jù)滲透能力在采掘前后的變化趨勢，進(jìn)行滲透能力參數(shù)分區(qū)，即劃定滲透能力變化趨勢一致的區(qū)段，進(jìn)行參數(shù)分區(qū)，分區(qū)內(nèi)的滲透能力參數(shù)進(jìn)行重新賦值，達(dá)到采掘擾動(dòng)下滲透能力變化的數(shù)值化處理目的。

(3)地面沉陷區(qū)數(shù)值化處理技術(shù)

地面沉陷降低了潛水水位埋深，易形成地面積水，從而加大了地下水資源的無效蒸發(fā)，因而在建模過程中應(yīng)予以考慮。如圖5所示，根據(jù)地面沉陷預(yù)測或?qū)崪y結(jié)果，通過改變地面高程值，對(duì)模型的上邊界進(jìn)行重新剖分。

圖5 地面沉陷數(shù)值化處理示意

需要注意的是在重新改變了上邊界條件后，須對(duì)模型進(jìn)行試運(yùn)行，分析沉陷區(qū)最低點(diǎn)的地面標(biāo)高Hmin與潛水位標(biāo)高Hw的關(guān)系，當(dāng)HminHw，地下水出露于地表形成積水，此時(shí)需對(duì)地面進(jìn)行蒸發(fā)能力的重新分區(qū)，積水區(qū)范圍按照水面蒸發(fā)處理，其他區(qū)域仍按照潛水蒸發(fā)處理；當(dāng)Hmin≥Hw，地下水未出露，全區(qū)仍為潛水蒸發(fā)，但是由于埋深的降低，潛水蒸發(fā)量可能增大。

綜上，以地下水系統(tǒng)數(shù)值模擬為基礎(chǔ)，通過將采煤擾動(dòng)形成的覆巖導(dǎo)水裂縫邊界化、彎曲帶滲透能力分區(qū)，地面沉陷區(qū)重新剖分的數(shù)值化處理手段是地下水系統(tǒng)數(shù)值評(píng)價(jià)模型構(gòu)建的基本思路。

3 評(píng)價(jià)案例

3.1 案例礦井概況

補(bǔ)連塔井田位于西部神東礦區(qū)中部，主采煤層埋深小于150 m，地表松散層厚度約為30 m，采用綜合機(jī)械化采煤工藝。如圖1所示，井田內(nèi)松散層類孔隙水主要接受大氣降水的入滲補(bǔ)給，在溝谷地下水出露于地表，形成常年或季節(jié)性的地表水體(補(bǔ)連溝、活雞兔等)，烏蘭木倫河是區(qū)內(nèi)所有地表水的排泄區(qū)。

如圖6所示，模型在垂向上按照巖性變化剖分為4層，研究區(qū)水文地質(zhì)單元較為完整，分別以常年流水的烏蘭木倫河、呼和烏蘇溝、活雞兔溝為模型的一類水頭邊界。

圖6 研究區(qū)(補(bǔ)連塔井田)范圍示意

3.2 采掘擾動(dòng)影響下地下水系統(tǒng)模擬

據(jù)礦井實(shí)測，補(bǔ)連塔井田采動(dòng)覆巖為典型的“兩帶”破壞形式，不存在彎曲帶滲透能力變異與地面沉陷積水問題，因此在評(píng)價(jià)模型中定義導(dǎo)水裂縫接觸帶為模型的“內(nèi)邊界”，概化為含水層地下水運(yùn)動(dòng)的一類水頭邊界條件，松散含水層水頭邊界水頭值為該處松散含水層底板標(biāo)高，降雨補(bǔ)給以研究區(qū)平水年數(shù)據(jù)為準(zhǔn)。

(1)水位

由圖7(a)分析可知(藍(lán)色箭頭為地下水流向，紅色箭頭為地下水轉(zhuǎn)地表水流向)，煤層未開采前，表現(xiàn)出大氣降水和地下水補(bǔ)給地表水的三水轉(zhuǎn)化特征，烏蘭木倫河為井田的地下水排泄點(diǎn)。局部由于地形控制，以地表分水嶺為界限地下水分別向兩側(cè)的補(bǔ)連溝、呼和烏素溝和活雞兔溝等地表水體排泄。

圖7 地下水流場變化特征

補(bǔ)連塔井田規(guī)?；_采開始于2000年左右，如圖7(b)所示，2012年井田基本以1-2煤回采完畢的三、四盤區(qū)為中心，形成較為明顯的降落漏斗，水位下降了20 m以上，采空區(qū)頂部含水層直接被疏干(黃色區(qū)域)。

(2)地下水均衡影響

表1中，未采掘前大氣降水是地下水補(bǔ)給的主要來源，側(cè)向補(bǔ)給次之，河流僅在局部地段補(bǔ)給地下水存在補(bǔ)給。地下水向河流等地表水體的排泄是主要排泄形式，蒸發(fā)排泄次之，由于模擬區(qū)水文地質(zhì)單元較為完整，地下水向區(qū)外側(cè)向排泄量為0。

表1 天然條件下地下水均衡

Table 1 Groundwater balance under natural conditions

源匯項(xiàng)數(shù)值占比/%合計(jì)降雨補(bǔ)給量/(104m3·d-1)6.7171.69河流補(bǔ)給量/(104m3·d-1)0.707.469.36側(cè)向補(bǔ)給量/(104m3·d-1)1.9520.85蒸發(fā)排泄量/(104m3·d-1)3.0732.829.36向河流排泄量/(104m3·d-1)6.2967.18

如圖8所示，隨著井田采掘范圍逐年增加，受采掘擾動(dòng)范圍擴(kuò)大，散含水層向采掘區(qū)段排泄的地下水水量逐年增大，以2012年為例，地下水向采掘區(qū)段的側(cè)向排泄量(即地下水的漏失量)為1.9×104m3/d，較之未采掘前，其排泄量主要是襲奪了河流排泄減量和潛水蒸發(fā)的減量。

圖8 采掘影響下水均衡變化

(3)對(duì)比分析

2012年補(bǔ)連塔井田四盤區(qū)的實(shí)測礦井正常排水量約為0.852×104m3/d(不含采空積水量)，由于礦井排水量不包括補(bǔ)給至采空區(qū)積水和井下復(fù)用等水量，其結(jié)果相對(duì)偏小。根據(jù)井田水文地質(zhì)特征，通過大井法計(jì)算地下水向采掘區(qū)的側(cè)向排泄量，其中松散含水層K取2.385 m/d，水頭高度、水位降深均取含水層厚度值30 m，F(xiàn)取截止2012年已采掘的三四盤區(qū)面積約為2.49×107m2，計(jì)算可得松散含水層向采掘區(qū)的地下水排泄量為4.4×104m3/d，由于大井法為四面匯水條件，無限含水層，規(guī)則影響半徑導(dǎo)致其計(jì)算結(jié)果明顯偏大，因此數(shù)值評(píng)價(jià)模型計(jì)算結(jié)果(1.9×104m3/d)能正確、定量地反映出采掘擾動(dòng)引起地下水的損失量。

4 結(jié) 語

(1)“三帶一區(qū)”是地下水賦存介質(zhì)擾動(dòng)的主要表現(xiàn)形式：由采動(dòng)“覆巖分帶”理論可知，采掘擾動(dòng)形成了覆巖的垮落帶、裂隙帶及彎曲帶。

(2)“三帶一區(qū)”空間尺度是控制地下水損失的主要因素：通過地下水動(dòng)力學(xué)分析可知，擾動(dòng)影響的面積越大、冒裂帶高度越高、彎曲帶保護(hù)層得到滲透能力越強(qiáng)則地下水滲漏量越大。

(3)提出了“三帶一區(qū)”模型化處理技術(shù)：在地下水系統(tǒng)數(shù)值模型中將導(dǎo)水裂縫切割含水層的接觸帶數(shù)值化處理成地下水系統(tǒng)的一類內(nèi)部邊界；根據(jù)彎曲帶滲透能力在采掘前后的變化趨勢，進(jìn)行滲透能力參數(shù)分區(qū)；根據(jù)地面沉陷區(qū)預(yù)測或?qū)崪y結(jié)果，通過改變地面高程值，對(duì)模型的上邊界進(jìn)行重新剖分。

(4)通過構(gòu)建神東礦區(qū)補(bǔ)連塔井田地下水?dāng)_動(dòng)的計(jì)算機(jī)評(píng)價(jià)模型，模擬得出由于采掘擾動(dòng)使井田地下水水位下降20 m以上(至2012年)，地下水資源損失量可達(dá)1.90×104m3/d。

[1] Chugh Y P,Behum P T.Coal waste management practices in the USA:An overview[J].International Journal of Coal Science & Technology,2014,1(2):163-176.

[2] Skousen J,Zipper C E.Post-mining policies and practices in the Eastern USA coal region[J].International Journal of Coal Science & Technology,2014,1(2):135-151.

[3] 胡振琪,龍精華,王新靜.論煤礦區(qū)生態(tài)環(huán)境的自修復(fù)、自然修復(fù)和人工修復(fù)[J].煤炭學(xué)報(bào),2014,39(8):1751-1757. Hu Zhenqi,Long Jinghua,Wang Xinjing.Self-healing,natural restoration and artificial restoration of ecological environment for coal mining[J].Journal of China Coal Society,2014,39(8):1751-1757.

[4] Manca P P,Desogus P,Orrù G.The reuse of abandoned Acquaresi mine voids for storage of the Masua flotation tailings[J].International Journal of Coal Science & Technology,2014,1(2):213-220.

[5] 范立民,馬雄德,蔣輝,等.西部生態(tài)脆弱礦區(qū)礦井突水潰沙危險(xiǎn)性分區(qū)[J].煤炭學(xué)報(bào),2016,41(3):531-536. Fan Limin,Ma Xiongde,Jiang Hui,et al.Risk evaluation on water and sand inrush in ecologically fragile coal mine[J].Journal of China Coal Society,2016,41(3):531-536.

[6] 王力,衛(wèi)三平,王全九.榆神府煤田開采對(duì)地下水和植被的影響[J].煤炭學(xué)報(bào),2008,33(12):1408-1414. Wang Li,Wei Sanping,Wang Quanjiu.Effect of coal exploitation on groundwater and vegetation in the Yushenfu Coal Mine[J].Journal of China Coal Society,2008,33(12):1408-1414.

[7] 潘天杭.用降落漏斗體積法求潛水含水層中抽水時(shí)之影響半徑[J].水文地質(zhì)工程地質(zhì),1957(9):41-42. Pan Tianhang.The influence radius of the water pumping in the aquifer by the method of the volume of the descent funnel method[J].Hydrogeology and Engineering Geology,1957(9):41-42.

[8] 張子文.穩(wěn)定降落漏斗中心最大下降及其體積計(jì)算方法探討[J].工程勘察,1989(3):34-36. Zhang Ziwen.Discussion on the maximum descent and its volume calculation method of the stable descend funnel center in[J].Geotechnical Investigation & Surveying,1989(3):34-36.

[9] 劉懷忠.煤礦開采對(duì)礦區(qū)地下水系統(tǒng)擾動(dòng)的定量評(píng)價(jià)研究[D].徐州：中國礦業(yè)大學(xué),2011. Liu Huaizhong.Quantitative evaluation of groundwater system disturbance in mining area[D].Xuzhou:China University of Mining and Technology,2011.

[10] 李楊.淺埋煤層開采覆巖移動(dòng)規(guī)律及對(duì)地下水影響研究[D].北京：中國礦業(yè)大學(xué)(北京),2012. Li Yang.Study on the movement law of overlying strata in shallow coal seam mining and its influence on groundwater[D].Beijing：China University of Mining and Technology(Beijing),2012.

[11] 李瑩.陜北煤炭分布區(qū)地下水資源與煤炭開采引起的水文生態(tài)效應(yīng)[D].西安：長安大學(xué),2008. Li Ying.The hydrological and ecological effects caused by groundwater resources and coal mining in the coal distribution area of Northern Shaanxi Province[D].Xi’an:Chang’an University,2008.

[12] 白樂,李懷恩,何宏謀,等.煤礦開采區(qū)地表水-地下水耦合模擬[J].煤炭學(xué)報(bào),2015,40(4):931-937. Bai Le,Li Huaien,He Hongmou,et al.Integrated simulation of surface water and groundwater in a high intensive coal mining area[J].Journal of China Coal Society,2015,40(4):931-937.

[13] 吳喜軍.煤炭開采地區(qū)河道徑流變化與生態(tài)基流研究-以陜北窟野河流域?yàn)槔齕D].西安：西安理工大學(xué),2013. Wu Xijun.Study on river runoff change and ecological base flow in coal mining area[D].Xi’an:Xi’an University of Technology,2013.

[14] 錢鳴高,石平五,許家林.礦山壓力與巖層控制[M].徐州:中國礦業(yè)大學(xué)出版社,2010:98-101. Qian Minggao,Shi Pingwu,Xu Jialin.Mining pressure and strata control[M].Xuzhou:China University of Mining and Technology Press,2010:98-101.

[15] 王雙明,黃慶享,范立民,等.生態(tài)脆弱礦區(qū)含(隔)水層特征及保水開采分區(qū)研究[J].煤炭學(xué)報(bào),2010,35(1):7-14. Wang Shuangming,Huang Qingxiang,Fan Limin,et al.Study on overburden aquclude and water protection mining regionazation in the ecological fragile mining area[J].Journal of China Coal Society,2010,35(1):7-14.

[16] Du W,Peng S,Zhu G,et al.Time-lapse geophysical technology-based study on overburden strata changes induced by modern coal mining[J].International Journal of Coal Science & Technology,2014,1(2):184-191.

[17] 隋旺華.開采沉陷土體變形工程地質(zhì)研究[M].徐州:中國礦業(yè)大學(xué)出版社,1999. Sui Wanghua.The engineering geology of the soil deformation of mining subsidence[M].Xuzhou:China University of Mining and Technology Press,1999.

[18] 李濤.陜北煤炭大規(guī)模開采含隔水層結(jié)構(gòu)變異及水資源動(dòng)態(tài)研究[D].徐州:中國礦業(yè)大學(xué),2012. Li Tao.Study on the structure variation and water resources of the large-scale mining of coal in the coal of Northern Shaanxi[D].Xuzhou:China University of Mining and Technology,2012.

Method of quantitative evaluation on amount of groundwater loss from uncon-fined aquifer caused by mining disturbance in the arid area of Western China

ZHAO Chun-hu1,2,HU Wei-yue1,2,JIN De-wu1,2

(1.Xi’anResearchInstituteofChinaCoalTechnology&EngineeringGroupCorp,Xi’an710054,China;2.ShaanxiKeyLaboratoryofPreventionandControlTechnologyforCoalMineWaterHazard,Xi’an710077,China)

The quantitative degree of coal resources exploitation on water environment disturbance is insufficient in the coal mining area of western China.In this paper,it concludes that the spatial scale is the main factor to control the loss of groundwater by the method of groundwater dynamics analysis,which includes the caving zone,the fracture zone,the bending zone and the ground subsidence area.The method of numerical evaluation model of groundwater system is proposed,which is based on the boundary of the hydraulic conductivity of overlying rock,the partition of bending zone,and the re-subdivision of surface subsidence area.According to the quantitative evaluation of groundwater dynamic system disturbance in the Bulianta coal mine,the simulation results show that the loss of groundwater resources is about 1.90×104m3/d due to the mining disturbance.

groundwater;loose aquifer;quantitative evaluation;numerical simulation；fractured zone

10.13225/j.cnki.jccs.2016.5001

2016-07-11

2016-11-09責(zé)任編輯:許書閣

國家自然科學(xué)基金資助項(xiàng)目(41402220)；國家“十三五”重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(2016YFC0501102)；中煤科工西安研究院自籌基金資助項(xiàng)目(2014MS009)

趙春虎(1980—),男,陜西扶風(fēng)人,副研究員,博士。E-mail:zhaochunhu@cctegxian.com

TD823

A

0253-9993(2017)01-0169-06

趙春虎，虎維岳，靳德武.西部干旱礦區(qū)采煤引起潛水損失量的定量評(píng)價(jià)方法[J].煤炭學(xué)報(bào),2017,42(1):169-174.

Zhao Chunhu,Hu Weiyue,Jin Dewu.Method of quantitative evaluation on amount of groundwater loss from unconfined aquifer caused by mining disturbance in the arid area of Western China[J].Journal of China Coal Society,2017,42(1):169-174.doi:10.13225/j.cnki.jccs.2016.5001