靳德武,周振方,趙春虎,馮龍飛,許 峰

(1.中煤科工集團西安研究院有限公司 礦山災(zāi)害防治與環(huán)境治理技術(shù)研發(fā)中心，陜西 西安 710054; 2.陜西省煤礦水害防治技術(shù)重點實驗室，陜西 西安 710077; 3.煤炭科學(xué)研究總院，北京 100013; 4.中煤科工集團西安研究院有限公司 水文地質(zhì)研究所，陜西 西安 710054)

我國侏羅系可采煤炭資源占煤炭總儲量的67%以上，主要分布在干旱半干旱、生態(tài)環(huán)境脆弱的西部地區(qū)[1]，先期煤炭資源以淺部開發(fā)為主[2]。煤層開采破壞覆巖結(jié)構(gòu)，形成采動導(dǎo)水裂隙，造成含水層、特別是淺表層水資源大量漏失，加劇了植被退化、土壤沙化等生態(tài)環(huán)境惡化進程[3-6]。

為化解西部礦區(qū)煤炭資源開發(fā)與水資源、生態(tài)環(huán)境保護之間的矛盾，20世紀(jì)90年代初期，王雙明、范立民團隊首先提出了“保水采煤”觀點[7-9]。圍繞這一觀點，眾多學(xué)者和專家開展了大量的科學(xué)研究和工程實踐，研究成果主要集中在以薩拉烏蘇組含水層為保水對象、以古近系紅黏土層為隔水關(guān)鍵層的保水采煤分區(qū)、導(dǎo)水裂隙帶控制方面[3,4,10-22]。保水采煤分區(qū)方面，王雙明等根據(jù)煤層頂板覆巖組合特征，將榆神府礦區(qū)劃分成自然保水開采區(qū)、可控保水開采區(qū)、保水限采區(qū)和無水開采區(qū)4種保水開采分區(qū)[3];孫亞軍等根據(jù)覆巖隔水層發(fā)育及埋藏條件，提出了>70 m厚基巖分布區(qū)上覆雙隔水層結(jié)構(gòu)、單隔水層結(jié)構(gòu)以及隔水層缺失結(jié)構(gòu)保水開采模式[16];根據(jù)植被與地下水關(guān)系，馬雄德等將榆神礦區(qū)劃分成植被約束區(qū)、地下水約束區(qū)和無約束區(qū)3種保水開采分區(qū)[23]。不管是哪種分區(qū)方法，受保護含水層水資源漏失通道(導(dǎo)水裂隙帶)的發(fā)育規(guī)律與控制技術(shù)是實現(xiàn)保水采煤的研究重點[10-11,24-25]。由此，黃慶享通過物理模擬和地裂縫實測分析，提出隔水巖組“上行裂隙”、“下行裂隙”之間的導(dǎo)通性決定了受保護含水層水量漏失程度的觀點，建立了隔水巖組隔水性判據(jù)[11];李文平等研究揭示了隔水關(guān)鍵層三趾馬紅土具有采動破裂后主動彌合、自修復(fù)的功能，從地質(zhì)結(jié)構(gòu)組成、水-土相互作用及其流變特性方面，闡述了其隔水性自恢復(fù)機理[14-15];之后，多位學(xué)者結(jié)合工程實踐，探討了控制煤層采高、充填采煤技術(shù)、隔水關(guān)鍵層注漿改造技術(shù)、含水層水轉(zhuǎn)移存儲和轉(zhuǎn)移存貯技術(shù)等[1,5,14-19,22,26-28]，形成了圍繞導(dǎo)水裂隙帶控制為主的保水采煤技術(shù)體系，為西部干旱礦區(qū)煤—水協(xié)調(diào)開發(fā)、生態(tài)環(huán)境保護提供了理論和技術(shù)參考。

近30 a來，圍繞“保水采煤”取得了豐碩的理論和技術(shù)成果，但保水采煤的根本是受保護含水層水資源量的保護，相關(guān)研究還較為少見，部分學(xué)者僅從礦區(qū)或井田尺度計算了采煤引起的水資源損失量[29-30]，以采煤工作面為水資源損失量計算的基本單元，定量計算工作面采動水量損失過程，是推動保水采煤理論定量化研究的關(guān)鍵。

筆者以榆神礦區(qū)典型覆巖結(jié)構(gòu)組合下煤層開采水量漏失過程為研究對象，通過覆巖破壞數(shù)值模擬，歸納不同的水量損失模式，構(gòu)建采動過程水量損失動力學(xué)模型，計算頂板水疏放和無疏放兩種條件下水量損失過程強度以及水資源損失總量，構(gòu)建頂板水疏放條件下工作面采動水量漏失模型，定量計算薩拉烏蘇組含水層和風(fēng)化基巖含水層水量損失比例。

1 研究區(qū)概況及覆巖破壞特征

陜北煤炭基地侏羅系上組煤開采覆巖結(jié)構(gòu)組合以沙層-土層-風(fēng)化層-基巖層類型為主，約占全區(qū)的65%,主要分布于榆神礦區(qū)[12]。研究區(qū)某煤礦主采延安組3-1煤層，平均煤厚3.0 m，平均埋深115 m，屬典型淺埋煤層。煤層覆巖結(jié)構(gòu)從上到下依次為新生代薩拉烏蘇組含水層、紅黏土隔水層、侏羅紀(jì)直羅組風(fēng)化基巖含水層以及正常基巖層。根據(jù)覆巖厚度分布特征，將井田開采煤層覆巖結(jié)構(gòu)概化為6種工況(表1)。

根據(jù)煤層覆巖組合類型，通過改變松散沙層、黏土層以及風(fēng)化基巖層厚度，設(shè)置30 m開挖步距，利用UDEC軟件分別模擬分析了工作面推采過程煤層頂板導(dǎo)水裂隙帶發(fā)育規(guī)律(圖1,表2)。模擬結(jié)果顯示，工作面推采至180 m以后，導(dǎo)水裂隙帶發(fā)育范圍達到穩(wěn)定;基巖段厚度與導(dǎo)水裂隙帶最大高度的大小關(guān)系，直接決定了煤層采動導(dǎo)水裂隙帶(上行裂隙)能否與下行裂隙溝通;厚基巖分布區(qū)，導(dǎo)水裂隙帶發(fā)育最大高度在45～50 m,模擬結(jié)果與現(xiàn)場實測45.72 m基本吻合。

表1 研究區(qū)3-1煤覆巖組合厚度特征Table 1 Thickness characteristics of overburden combination of No.3-1 coal seam in study area

注:表中巖層厚度是根據(jù)井田水文地質(zhì)補勘報告進行的覆巖特征概化，并非嚴(yán)格的實際厚度。

圖1 不同覆巖組合煤層開采導(dǎo)水裂隙帶發(fā)育高度范圍Fig.1 Development height range of water conducting fracture zone in coal seam mining with different overburden rock combinations simulated by the UDEC software

總的來看，根據(jù)煤層開采上行裂隙和下行裂隙發(fā)育特征，頂板含水層水量損失可分為2種模式，當(dāng)上行和下行裂隙在隔水黏土層未能溝通，頂板含水層水量損失主要來源于單一的風(fēng)化基巖含水層(模式1);兩種裂隙突破隔水黏土層后相互溝通，形成良好的垂向過水通道，水量損失來源于風(fēng)化基巖和松散層砂巖復(fù)合含水層(模式2)。

2 水量損失動力學(xué)模型構(gòu)建

2.1 頂板水疏放工作面推采水量損失過程動力學(xué)模型構(gòu)建

一般，工作面采前均經(jīng)過頂板水預(yù)疏放[31]，根據(jù)前期研究成果，煤層開采引起的頂板水量損失主要來源于靜儲量釋放和采空區(qū)大井動態(tài)水量損失，兩種形式的水量損失過程可分別用一階指數(shù)衰減模型(式(1))和一階指數(shù)恢復(fù)模型表征(式(2))，相應(yīng)工作面頂板水量損失過程可用二階動力系統(tǒng)模型描述(式(3))[32]。

表2 工作面推采過程覆巖導(dǎo)水裂隙帶發(fā)育最大高度Table 2 Maximun altitude of water-conducting fracture zone in overburden rock with mining m

工作面頂板水量損失過程動力系統(tǒng)模型表述為

Qj=A1exp(-k1at/b)+C1

(1)

Qd=C2-A2exp(-k2at/b)

(2)

Q=Qj+Qd=A1exp(-k1at/b)-

A2exp(-k2at/b)+C

(3)

其中，Qj和Qd分別為at/b位置頂板靜儲量釋放水量和采空區(qū)大井動態(tài)損失量，m3/h;at/b為推采t時刻工作面推采走向長度與傾向?qū)挾鹊谋戎?簡稱“走傾比”);k1和k2分別為靜儲量釋水衰減系數(shù)和大井動態(tài)損失增長系數(shù)，m3/(h·1);C=C1+C2，C1和C2分別為頂板靜儲量釋放水量動態(tài)穩(wěn)定值和大井水量損失動態(tài)穩(wěn)定值，m3/h;根據(jù)模型意義，容易推導(dǎo)A1=Qj0-Qje，Qj0為推采直接頂初次垮落時刻位置釋放靜儲量，可表述為Qj0=bvhc，v為直接頂初次垮落位置推采速度，m/h;hc為直接頂初次垮落位置殘余水頭高度，m。根據(jù)模型意義，初始時刻至直接頂初次垮落時刻位置，頂板等效于無破壞，大井動態(tài)水量損失為0，at/b趨于無窮大時，Qd=C2，據(jù)此，容易推導(dǎo)A2>C2;A1和A2分別為工作面推采靜儲量釋放和動態(tài)補給初始參量。

工作面推采頂板動、靜儲量疊加水量損失過程，一般呈現(xiàn)兩種變化趨勢:當(dāng)采動導(dǎo)水裂隙揭露覆巖單一含水層時，地下水損失量一般呈現(xiàn)先增加至峰值，后逐步衰減至動態(tài)穩(wěn)定水量的趨勢;當(dāng)采動導(dǎo)水裂隙揭露覆巖復(fù)合含水層時，地下水水量損失過程呈現(xiàn)逐步增加至動態(tài)穩(wěn)定或基本趨于動態(tài)穩(wěn)定的變化趨勢。據(jù)此，工作面推采頂板含水層水量損失過程一般變化規(guī)律如圖2所示。

圖2 工作面推采頂板水量損失組成及變化過程Fig.2 Composition and change of water-inflow from roof aquifer with the working face mining process

2.2 無疏放狀態(tài)工作面推采頂板含水層水量損失模型構(gòu)建

工作面推采頂板水量損失規(guī)律是推采速度、頂板覆巖組合條件以及含水層水文地質(zhì)特征共同作用的結(jié)果[32-33]。根據(jù)采空區(qū)實測水量損失數(shù)據(jù)擬合得到的曲線方程僅僅反映的是頂板水疏放之后的水量損失過程變化規(guī)律，對于不同頂板水疏放程度工作面推采水量損失量的求解意義不大。無疏放狀態(tài)水量損失過程方程對于計算頂板水頭控制疏降工作面水量損失強度具有重要意義，等效于現(xiàn)狀頂板水疏放條件下水量損失過程與前期鉆孔疏放水量疊加的結(jié)果。對于同一工作面，無疏放和有疏放兩個水量損失過程具有相同的初始和動態(tài)平衡位置，以及對應(yīng)相同的特征水量損失強度(圖3)。據(jù)此，2種疏放狀態(tài)之間的水量損失強度定量關(guān)系可以用定積分方程表述(式(4))。

(4)

式中，QT為天然無疏放狀態(tài)工作面推采水量損失曲線方程;QS為現(xiàn)狀疏放狀態(tài)工作面推采水量損失曲線方程;ΔQ為鉆孔疏放水總量;ate為工作面推采水量損失動態(tài)平衡位置;te為工作面推采水量損失動態(tài)平衡時間。

圖3 工作面推采水量損失定積分概念模型Fig.3 Definite integral computing model for water resources loss with the working face mining process

3 水量損失模型求解及應(yīng)用

3.1 頂板水疏放工作面推采覆巖含水層水量損失規(guī)律

按照工作面覆巖結(jié)構(gòu)特征進行劃分，M1，M2，M3工作面屬于模式1，N1，N2，N3，N4，N5，N6，N7工作面屬于模式2，工作面推采頂板水量損失過程表現(xiàn)出了顯著差異。單一風(fēng)化基巖充水條件，水量損失強度呈現(xiàn)初期增加至峰值水量，而后隨推采逐步下降至動態(tài)平衡的趨勢(圖4);與之相比，風(fēng)化基巖與松散層復(fù)合充水條件，煤層上覆直羅組風(fēng)化基巖相對較薄，煤層開采導(dǎo)水裂隙發(fā)育至松散層，工作面推采水量損失強度呈現(xiàn)逐步增加的趨勢，但增速逐步下降，基本趨近動態(tài)穩(wěn)定(圖5)。

圖4 模式1工作面推采水量損失過程擬合曲線Fig.4 Variations of water resources loss with the working face mining process described using the double-exponential decay model (Mode 1)

基于麥克勞林公式指數(shù)多項式轉(zhuǎn)化為非線性多項式近似計算的基本思想[32,34]，求解了2種模式工作面覆巖含水層水量損失過程曲線方程相應(yīng)參數(shù)(表3)。兩種模式相比，頂板靜儲量衰減速率和采空區(qū)大井動態(tài)損失強度增長速率存在數(shù)量級的差異，指示模式2采空區(qū)水量損失強度的持續(xù)增加主要受到復(fù)合含水層上段松散層潛水持續(xù)疊加補給影響;而模式1單一含水層充水，前期逐段靜儲量釋水量較大，控制著涌水過程逐步增加的趨勢，后期隨著靜儲量釋水量逐步下降，動態(tài)損失水量開始占據(jù)主導(dǎo)地位，控制著水量損失過程緩慢增加的趨勢，這種增速與前期靜儲量的逐段集中釋放相比，基本可以忽略，這與吳吉春等的研究結(jié)果一致[35]。相同模式，特征水量損失強度及發(fā)生位置存在差異，且二者存在正比例關(guān)系，指示導(dǎo)水裂隙波及層段巖層綜合富水性是控制工作面推采頂板水量損失過程的主要因素，這與其他學(xué)者研究結(jié)果是一致的[36-37]。

3.2 無疏放狀態(tài)工作面推采頂板水量損失規(guī)律

選取有鉆孔疏放水觀測記錄的2個工作面(M2和N5工作面)，進行2種模式無疏放狀態(tài)工作面推采頂板水量損失過程線求解。M2工作面疏放水總量128.90萬m3，根據(jù)式(4)列出實際疏放與無疏放工作面推采頂板水量損失過程線之間的定積分關(guān)系式:

圖5 模式2工作面水量損失過程擬合曲線Fig.5 Variations of water resources loss with the working face mining process described using the double-exponential decay model (Mode 2)

工作面A1/(m3·h-1)A2/(m3·h-1)k1/(m3·(h·1)-1)k2/(m3·(h·1)-1)Qe/(m3·h-1)ate/bQm/(m3·h-1)atm/bR2M1518.53 992.30 3.19 14.71 46.67 2.323 268.87 0.189 0.79模式1M2473.37 608.99 6.48 12.73 91.83 1.131 180.70 0.148 0.64M3(a)523.03 475.71 0.12 3.28 42.17 28.500 488.18 1.017 0.84N1385.20 560.03 0.20 0.20 180.00 29.300? ——0.91N2450.20 533.03 1.00 0.92 115.00 6.850? ——0.24N3405.20 443.13 0.20 0.20 160.00 21.600? ——0.22模式2N4583.60 672.93 0.31 0.37 170.00 21.550? ——0.52N5565.60 729.59 0.20 0.22 188.00 12.220? ——0.74N6463.60 666.23 0.30 0.26 290.00 26.600? ——0.89N7175.20 490.37 0.15 0.10 390.00 68.800? ——0.83

注:(a)為M3工作面推采水量損失強度觀測樣本較少，代表性弱;*為曲線斜率等于0.05時近似為水量損失平衡位置。

(5)

式中，QTM為天然無疏放狀態(tài)M2工作面推采水量損失曲線方程;QSM為現(xiàn)狀疏放狀態(tài)M2工作面推采水量損失曲線方程。

列出現(xiàn)狀疏放狀態(tài)M2工作面推采頂板水量損失過程線方程，以及無疏放狀態(tài)M2工作面推采頂板水量損失過程線一般模型:

QSM=473.37 exp(-6.48at/b)-

608.99 exp(-12.73at/b)+91.83

(6)

QTM=A1TMexp(-k1TMat/b)-A2TMexp(-k2TMat/b)+

CTM

(7)

式中，A1TM和A2TM分別為無疏放狀態(tài)M2工作面推采靜儲量和動態(tài)補給初始參量;k1TM和k2TM分別為靜儲量釋水衰減系數(shù)和采空區(qū)大井動態(tài)損失增長系數(shù);CTM為無疏放狀態(tài)M2工作面推采頂板靜儲量釋放水量動態(tài)穩(wěn)定值和大井水量損失動態(tài)穩(wěn)定值的疊加。

對式(5)進行化簡，變形得到關(guān)系式:

(8)

根據(jù)模型概化結(jié)果，對于M2工作面，實際疏放和無疏放工作面推采水量損失過程線存在相同的初值、動態(tài)平衡值，以及水量損失峰值發(fā)生位置，據(jù)此，列出上述特征值數(shù)學(xué)表達式:

QTM|at/b=0=QSM|at/b=0

(9)

QTM|at/b=1.131=QSM|at/b=1.131

(10)

(QTM|at/b=0.148)′=0

觀念是支配人的行為的主觀意識。在科技成果向教學(xué)資源轉(zhuǎn)化之前，科技成果的擁有者要樹立正確的觀念。只有充分認識到科技成果教學(xué)轉(zhuǎn)化的重要意義，才能為科技成果教學(xué)轉(zhuǎn)化提供動力，這種動力就是科技成果擁有者的內(nèi)在動力。但是，無論是社會還是高校內(nèi)部，都沒有重視科技成果教學(xué)轉(zhuǎn)化，對其重要意義的認識也不足，其原因主要有以下幾個方面。

(11)

經(jīng)化簡、變形得到上述條件的對應(yīng)方程:

A1TM-A2TM=-135.62

(12)

A1TMexp(-1.131k1TM)-A2exp(-1.131k2TM)=0

(13)

-A1TMk1TMexp(-0.148k1TM)+

A2TMk2TMexp(-0.148k2TM)=0

(14)

聯(lián)立式(8)，(12)～(14)，求解得到M2工作面無疏放狀態(tài)推采過程水量損失過程線模型(式(15)，圖6)。

QTM=3 786.96exp(-4.527at/b)-3 922.58×

exp(-10.356at/b)+91.83

(15)

頂板水預(yù)疏放顯著降低了工作面推采過程水量損失峰值強度，從M2工作面來看，峰值水量從1 182.56 m3/h減至180.70 m3/h，減少了84.72%(圖6)，無疏放狀態(tài)頂板水靜儲量的集中釋放導(dǎo)致了水量損失陡增陡降過程。從安全角度來看，頂板水預(yù)疏放起到了顯著的“削峰”作用，使工作面推采水量變化過程變得相對平緩。從水資源損失方面來看，在煤層開采不可避免破壞覆巖含水層時，頂板預(yù)疏放部分水資源能夠避免直接進入采空區(qū)受到機械殘油污染的風(fēng)險，從而得到轉(zhuǎn)移利用、轉(zhuǎn)移存儲或轉(zhuǎn)移存貯[4,16]。

圖6 M2工作面推采水量損失過程線Fig.6 Process line of water loss with the M2 working face mining

與M2工作面不同的是，頂板水疏放后，模式2 N5工作面水量損失過程呈現(xiàn)逐步增加的趨勢，無疏放狀態(tài)過程線求解無法準(zhǔn)確給出水量損失峰值發(fā)生位置，以及對應(yīng)數(shù)學(xué)表達式。以往研究成果顯示，覆巖水文地質(zhì)條件相似工作面推采頂板水量損失過程峰值水量及相應(yīng)位置主要受控于推采速度和導(dǎo)水裂隙帶發(fā)育高度[32]。覆巖破壞數(shù)值模擬結(jié)果指示模式2導(dǎo)水裂隙帶發(fā)育在走傾比在0.5～0.6位置達到穩(wěn)定，這與模式一導(dǎo)水裂隙帶發(fā)育穩(wěn)定位置類似;同時，N5和M2工作面初期推采速度分別是177.6 m/月和167.8 m/月，基本一致?？梢酝茰y，無疏放狀態(tài)2個工作面推采峰值水量損失位置走傾比基本接近(<0.2)，筆者通過假定N5工作面水量損失峰值強度走傾比分別在0.2,1.0,3.0三個位置，進行相應(yīng)近似過程線求解，找出水量損失峰值強度與發(fā)生位置之間的數(shù)學(xué)近似關(guān)系。

實測N5工作面頂板水預(yù)疏放水量256.29萬m3，參照M2工作面水量損失過程線求解方法，得到3種假設(shè)峰值位置水量損失過程線方程(式(16)～(18)，圖7)。

QTN1=4 524.80exp(-0.838at/b)-

4 688.79exp(-15.886at/b)+188.00

(16)

QTN2=3 673.80exp(-0.510at/b)-

3 738.79exp(-1.834at/b)+188.00

(17)

QTN3=1 414.00exp(-0.179at/b)-

1 577.99exp(-0.821at/b)+188.00

(18)

其中，QTN1,QTN2,QTN3分別為水量損失峰值強度在0.2，1.0和3.0三個位置條件下天然無疏放狀態(tài)N5工作面推采水量損失曲線方程。

圖7 N5工作面推采水量損失過程線Fig.7 Process line of water loss with the N5 working face mining

3種假設(shè)位置對應(yīng)水量損失峰值強度分別為3 819.05，1 796.74，880.07 m3/h，可以看出，瞬時水量損失峰值強度隨著其發(fā)生位置的推移不斷減小，且呈現(xiàn)指數(shù)下降關(guān)系(圖7)。工作面推采前期，峰值水量發(fā)生位置越小，水量集中損失強度越顯著，由于水量損失總量相同，推采后期涌水量呈現(xiàn)相反的特點，即水量損失峰值強度發(fā)生位置越大，衰減后期頂板水量損失瞬時值越大。按照研究區(qū)峰值水量發(fā)生位置<0.2的規(guī)律，無疏放條件下N5工作面推采頂板水量損失峰值強度>3 819.05 m3/h。

3.3 工作面推采覆巖含水層水量損失分析

利用定積分方法求解了無疏放狀態(tài)水量損失過程線方程，得到了工作面推采水量損失峰值。單一充水水源條件M2工作面和復(fù)合充水水源條件N5工作面，頂板水頭分別平均疏降68.42%,66.67%(圖8)，水量損失峰值分別下降84.72%,95.08%，疏放水后，工作面推采過程頂板水量分別損失74.97,94.82萬m3，進入采空區(qū)水量分別減少41.84%,72.99%，結(jié)果指示鉆孔預(yù)疏放對于復(fù)合充水水源條件工作面推采頂板涌水“削峰”效果更為顯著。在風(fēng)化基巖承壓含水層儲水量相似的條件下，采煤分別造成M2和N5工作面頂板含水層總的水量損失203.87,351.11萬m3。與N5工作面相比，M2工作面開采，古近系黏土層至少阻隔了140萬m3松散層水進入采空區(qū)。

在工作面推采水量損失過程強度和損失總量計算結(jié)果基礎(chǔ)上，根據(jù)工作面覆巖與采動裂隙空間結(jié)構(gòu)特征，構(gòu)建了M2和N5工作面采動頂板水量損失強度數(shù)學(xué)和數(shù)值模型(圖9)，進一步計算頂板水疏放后單位走向長度工作面覆巖含水層水量損失強度。

根據(jù)地下水動力學(xué)基本原理，導(dǎo)水裂隙帶外圍非采動擾動區(qū)為達西滲流區(qū)域(式(19))。導(dǎo)水裂隙帶范圍以內(nèi)含水層由于被迅速疏干，形成導(dǎo)水裂隙帶與含水層空間接觸曲面，曲面水壓力p為大氣壓，定義為0(式(20))，基于Comsol Multiphysics軟件進行了工作面回采過程單位走向長度覆巖地下水損失過程穩(wěn)定流模擬。

非采動擾動區(qū)Darcy線性滲流方程:

(19)

導(dǎo)水裂隙帶與含水層空間接觸曲面:

p=0

(20)

式中，u為流速;p為水壓;k′為含水層滲透率。

單位走向長度水量損失強度薩拉烏蘇組潛水和直羅組風(fēng)化基巖承壓水占比計算結(jié)果指示，2種覆巖組合模式相比，采動導(dǎo)水裂隙揭露不同的覆巖含水層，水量損失規(guī)律差異明顯，采動導(dǎo)水裂隙揭露的隔水保護層厚度越大，松散含水層失水強度與占比越高(表4)。工作面推采水量損失總量、水量損失強度計算結(jié)果，揭示了研究區(qū)古近系黏土隔水層對于保護薩拉烏蘇組松散層潛水的重要意義。

表4 單位走向長度水量損失強度計算結(jié)果Table 4 Characteristic intensity of water loss with mining

注:水量損失強度指水文地質(zhì)特征概化后每百米覆巖含水層水量損失強度。

4 結(jié) 論

(1)基于工作面推采頂板水量損失二階動力系統(tǒng)模型及麥克勞林公式，求解了頂板水實際疏放條件下單一和復(fù)合含水層水量損失10個工作面的曲線方程，計算了水量損失特征值及發(fā)生位置，結(jié)果指示導(dǎo)水裂隙帶波及層段巖層綜合富水性是水量損失變化過程的主控因素。

(2)構(gòu)建了無疏放狀態(tài)工作面推采頂板水量損失數(shù)學(xué)模型，利用定積分方法確定了M2和N5工作面無疏放條件水量損失過程方程，結(jié)果顯示鉆孔預(yù)疏放使工作面水量損失峰值顯著下降，為減少水資源污染，實施水資源轉(zhuǎn)移存儲方案提供了可能。

(3)利用定積分模型、含水層水量損失強度數(shù)值模型分別計算了工作面推采頂板水資源損失總量、薩拉烏蘇組潛水和直羅組風(fēng)化基巖承壓水損失強度占比，定量分析了榆神礦區(qū)古近系黏土層保護薩拉烏蘇組潛水的重要意義。