趙瑞玨，毛德強，劉再斌，姬中奎，曹祖寶

（1.山東大學(xué)土建與水利學(xué)院，山東 濟(jì)南 250061；2.中煤科工集團(tuán)西安研究院有限公司，陜西 西安 710054；3.長安大學(xué)環(huán)境科學(xué)與工程學(xué)院，陜西 西安 710054）

在煤炭開采的過程中，礦井突水事故是影響安全生產(chǎn)的重要因素之一[1]。為了探明同一含水層內(nèi)部或不同含水層之間的水力聯(lián)系，深入了解含水層富水性以及補給、徑流、排泄條件，通常進(jìn)行同一礦區(qū)大流量、大降深、大范圍的放水試驗[2-3]。查清本井田及相鄰井田的水文地質(zhì)條件，為煤礦安全生產(chǎn)提供技術(shù)保障[4]。

傳統(tǒng)放水試驗分析方法得到的含水層水文地質(zhì)參數(shù)是含水層的等效參數(shù)，不能準(zhǔn)確刻畫參數(shù)的空間異質(zhì)性。邵紅旗等[5]提出了以水量衰減率、單位涌水量的水位降深值及水位恢復(fù)速率為主要指標(biāo)的放水試驗分析方法；曹祖寶等[6]利用“大井法”分析放水試驗，得到的水文地質(zhì)參數(shù)是含水層的等效參數(shù)。

水力層析法是基于抽水或注水試驗的改進(jìn)方法，通過一系列交叉抽水或注水試驗，運用水頭響應(yīng)數(shù)據(jù)，識別含水層水文地質(zhì)參數(shù)的空間分布特征[7]。在同樣抽水井?dāng)?shù)量的情況下，水力層析法能夠獲得更多有關(guān)含水層的水頭響應(yīng)數(shù)據(jù)，分析結(jié)果能夠更加準(zhǔn)確地反映含水層水文地質(zhì)參數(shù)的空間分布特征，求解結(jié)果更加接近實際。蔣立群等[8]基于室內(nèi)非均質(zhì)含水層砂箱實驗，研究了水力層析法和克立金法刻畫非均質(zhì)含水層滲透系數(shù)場的優(yōu)劣；郝永紅等[9-10]利用數(shù)值算例系統(tǒng)分析了水力層析法可準(zhǔn)確獲取含水層的裂隙帶分布信息；Illman 等[11]基于兩次抽水試驗數(shù)據(jù)，運用水力層析法建立了日本Mizunami 地下實驗室水文地質(zhì)參數(shù)的三維分布圖，揭示了斷層導(dǎo)水特征；董艷輝等[12]應(yīng)用水力層析法研究了離散裂隙網(wǎng)絡(luò)的連通性，獲取的滲透特性分布可以更好地預(yù)測滲流規(guī)律；Fischer 等[13]通過結(jié)構(gòu)約束模型推斷了巖溶地區(qū)水力特性的空間分布。但目前為止，還未有利用煤礦放水試驗獲取含水層水文地質(zhì)參數(shù)分布信息的相關(guān)研究。

在本研究中，運用水力層析法分析煤礦井下疊加放水試驗數(shù)據(jù)，反演計算含水層水文地質(zhì)參數(shù)，并根據(jù)互相關(guān)分析具體解釋水力層析法識別含水層參數(shù)空間分布特征的原理。該方法操作簡單，適用性廣，獲取的參數(shù)詳細(xì)具體，可準(zhǔn)確刻畫含水層的非均質(zhì)性，精確定位富水異常區(qū)的位置及范圍，有效降低突水事故風(fēng)險，為煤礦放水試驗數(shù)據(jù)分析及礦井水文地質(zhì)條件評價提供一套新方法。

1 水力層析法原理

同步連續(xù)線性估計算法（Simultaneous Successive Linear Estimator，簡稱SimSLE）發(fā)展成熟，且已得到廣泛應(yīng)用[14-16]。本研究基于SimSLE 算法，結(jié)合煤礦井下疊加放水試驗，具體解釋水力層析法的基本原理。

1.1 地下水流控制方程

地下水流控制方程用于模擬地下水流場變化，是進(jìn)行水文地質(zhì)參數(shù)反演的基礎(chǔ)。對于飽和帶非均質(zhì)各向同性介質(zhì)中的非穩(wěn)定承壓水問題，地下水流控制方程為：

邊界條件和初始條件分別為：

式中：H—水頭/m；

x—二維空間坐標(biāo)/{m,m}；

Q（xp,t）—xp處t時刻的放水孔流量/（m·s-1）；

T（x）—導(dǎo)水系數(shù)/（m2·s-1）；

K（x）—滲透系數(shù)/（m·s-1）；

S（x）—儲水系數(shù)；

H1—Dirichlet 邊界Γ1的水頭/m；

q—Neumann 邊界Γ2單位法向量n方向的流量/（m·s-1）；

H0—含水層初始水頭/m。

1.2 同步連續(xù)線性估計算法

SimSLE是一種基于地質(zhì)統(tǒng)計學(xué)的參數(shù)估計方法[17]。不同于傳統(tǒng)放水試驗分析中的均質(zhì)性假設(shè)，SimSLE將導(dǎo)水系數(shù)和儲水系數(shù)的自然對數(shù)及觀測孔水頭處理成空間隨機(jī)變量：

f'、s'、h'—擾動值。

SimSLE 從帶有權(quán)重的線性估計開始，利用直接測量的含水層水文地質(zhì)參數(shù)值和放水試驗過程中的水頭測量值，計算得到參數(shù)的初始估計值。權(quán)重根據(jù)水頭測量值與水文地質(zhì)參數(shù)的統(tǒng)計屬性，如均值、相關(guān)長度、水頭測量值和參數(shù)之間的互協(xié)方差和自協(xié)方差，利用協(xié)同克里格（cokriging）方法計算得到。在估計過程中，參數(shù)的自協(xié)方差不斷更新，代表估計不確定性的變化[18-19]。以導(dǎo)水系數(shù)為例，初始估計結(jié)果為：

H*—實際水頭測量值；

He—通過地下水流控制方程求得的水頭模擬值；

Nf—參數(shù)測量點的總數(shù)；

Np—放水試驗的總次數(shù)；

Nh（k）—第k次放水試驗中水頭測量點的總數(shù)；

Nt（k,j）—第k次放水試驗中第j個水頭測量點隨時間變化的水頭測量數(shù)據(jù)的總數(shù)；

λ0i—協(xié)同克里格權(quán)重，即第i個參數(shù)測量點的參數(shù)值f*對x0處的參數(shù)估計值的權(quán)重；

μ0kjl—實際水頭測量值H*（k,xj,tl）對參數(shù)估計值的權(quán)重。

將參數(shù)估計值代入地下水流控制方程，計算在水頭測量點不同時刻的模擬值，判斷模擬值與實際測量值之間的差異，若差異沒有達(dá)到指定精度，上述計算過程會重復(fù)迭代進(jìn)行。

1.3 互相關(guān)分析

運用水力層析法時，不同位置的放水孔和觀測孔水頭響應(yīng)數(shù)據(jù)反映了含水層不同區(qū)域的水文地質(zhì)參數(shù)信息。通過分析隨時間變化的水頭測量值與導(dǎo)水系數(shù)和儲水系數(shù)的互相關(guān)性，可以進(jìn)一步解釋水力層析法識別含水層空間異質(zhì)性的原理[20]。處的一階泰勒公式為擾動項的表達(dá)式為：

式中：J—雅可比（Jacobian）系數(shù)。

假設(shè)導(dǎo)水系數(shù)和儲水系數(shù)相互獨立，式（9）兩邊同時右乘f′，并求期望，且利用矩陣形式表達(dá)：

式中：Rhf—h與f的互協(xié)方差矩陣；

Rff—f的自協(xié)方差矩陣。含水層參數(shù)的空間變異性通過自協(xié)方差矩陣中的方差和相關(guān)長度表示；

Rhs—h與s的互協(xié)方差矩陣；

Rss—s的自協(xié)方差矩陣。

h′

式（9）兩邊同時右乘，并求期望，代入式（10）、（11）：

式中：Rhh—h的自協(xié)方差矩陣。

由相關(guān)系數(shù)的定義，代入式（10）、（11）：

任意位置的水頭H（x,t）在均值

式中：ρ—互相關(guān)系數(shù)矩陣；

σ2h、σ2f、σ2s—h、f、s的方差，即自協(xié)方差矩陣Rhh、Rff、Rss的對角線元素。

互相關(guān)系數(shù)矩陣可表征空間任意位置水文地質(zhì)參數(shù)的不確定性對觀測孔水頭測量值的影響?；ハ嚓P(guān)分析是在隨機(jī)背景下進(jìn)行的敏感性分析，采用了地質(zhì)統(tǒng)計學(xué)的概念，通過期望和方差代入水文地質(zhì)參數(shù)的先驗信息（prior information），綜合分析水文地質(zhì)參數(shù)的空間異質(zhì)性對水頭的影響。

互相關(guān)系數(shù)的取值范圍為[-1,1]?；ハ嚓P(guān)系數(shù)大于零，即為正相關(guān)，表現(xiàn)為當(dāng)空間某一位置水文地質(zhì)參數(shù)增大時，水頭測量值相應(yīng)增大；反之，互相關(guān)系數(shù)小于零，即為負(fù)相關(guān)，表現(xiàn)為水文地質(zhì)參數(shù)增大時，水頭測量值相應(yīng)減小。互相關(guān)系數(shù)的絕對值越大，相關(guān)性越強，表明若賦予水文地質(zhì)參數(shù)一個微小的擾動值，水頭測量值相應(yīng)地改變較大，說明對于該處水文地質(zhì)參數(shù)的變化比較敏感，該觀測孔的水頭響應(yīng)數(shù)據(jù)對于反演水文地質(zhì)參數(shù)效果顯著。

2 研究區(qū)概況

研究區(qū)檸條塔井田位于陜西省榆林市境內(nèi)。煤系地層地表露頭近乎水平展布，下伏煤層底板向西傾斜，坡降為4‰～10‰。在以往勘探及礦井建設(shè)和到目前為止的生產(chǎn)過程中，均未發(fā)現(xiàn)較大的斷層（落差大于15 m）或褶皺，僅煤層底板有一些寬緩的波狀起伏[21-22]。

據(jù)地質(zhì)填圖及鉆孔揭露，檸條塔井田地層由老至新（自下而上）分別為：侏羅系中統(tǒng)延安組（J2y）、直羅組（J2z），新近系上新統(tǒng)保德組（N2b），第四系中更新統(tǒng)離 石組（Q2l），第 四系全新統(tǒng)風(fēng) 積層（Q4el）和沖 積層（Q4al）。延安組為主要含煤地層；直羅組為砂巖含水層；保德組為相對隔水層，結(jié)構(gòu)較致密，滲透性較弱。開采延安組2-2煤層的S1210 工作面涌水量接近1 000 m3/h，維持了三個多月，涌水來源判斷為直羅組砂巖含水層。涌水期間，直羅組上覆和下覆地層的觀測孔未發(fā)生相應(yīng)的水頭變化，因此本研究將直羅組含水層概化為二維含水層。

3 放水試驗與數(shù)值建模

以檸條塔煤礦水文地質(zhì)勘探過程中的井下疊加放水試驗為研究內(nèi)容，依據(jù)水力層析技術(shù)，利用基于地質(zhì)統(tǒng)計學(xué)的同步連續(xù)線性估計算法SimSLE，對水頭數(shù)據(jù)進(jìn)行定量解譯，并對侏羅系中統(tǒng)直羅組（J2z）的含水層結(jié)構(gòu)進(jìn)行模擬計算。

3.1 井下疊加放水試驗

在井田南翼較大范圍內(nèi)，工作面推采前方的水文地質(zhì)條件尚未完全探明，需進(jìn)一步查明是否存在富水異常區(qū)。工作面井下施工鉆孔數(shù)量較多，自切眼往外約2 000 m 范圍內(nèi)工作面頂板有幾十個鉆孔（圖1），斜向上全部貫穿直羅組含水層直至保德組黏土層。利用井下鉆孔進(jìn)行采前大流量放水試驗。第四系地面觀測孔（圖1Q5）測量放水試驗過程中離石組含水層的水頭變化，探查工作面涌水與第四系水的水力聯(lián)系，測量數(shù)據(jù)表明放水過程中沒有明顯的水頭變化，故進(jìn)一步確定保德組為相對隔水層，放水試驗與第四系水無水力聯(lián)系。

放水試驗共持續(xù)9 d，放水次序為自南向北，采用疊加放水的方式。根據(jù)放水前對井下所有鉆孔涌水量的排查情況，共選取工作面頂板上30個水量較大的鉆孔作為放水孔。分區(qū)域選取28個井下鉆孔（其中工作面頂板鉆孔17個，巷道外側(cè)鉆孔11個）和6個地面鉆孔作為水頭觀測孔。

圖1 鉆孔分區(qū)布置平面圖Fig.1 Plan of the boreholes zoning scheme

34個觀測孔全程測量水頭。試驗正式開始時，第一階段開啟F3-1、F5、F5-1、F6、F6-1、F7 共6個放水孔閥門（圖1）。第一階段放水持續(xù)2 d。第二階段開啟T6、F8、F9、F10、F11、F12、F12-2 共7個放水孔閥門（圖1）。此時共有13個鉆孔聯(lián)合放水。第二階段放水持續(xù)2 d。第三階段開啟T9、F15、F16、F17-2、F19、F19-1、F20-1、F21-2 共8個放水孔閥門（圖1）。此時共有21個鉆孔聯(lián)合放水。第三階段放水持續(xù)3 d。第四階段開啟YZ10-2、YZ11-1、F22、F22-1、F24、F26-1、F26-2、F28、F29 共9個放水孔閥門（圖1）。此時共有30個鉆孔聯(lián)合放水。第四階段放水持續(xù)2 d。每個階段的穩(wěn)定放水量分別為270，214，167，96 m3/h，最后穩(wěn)定放水總量為747 m3/h。按照1 h時間步長，將采集到的流量數(shù)據(jù)分段累加，并保持水量不變，換算至整小時時刻處。

在放水的四個階段分別選取一個井下觀測孔，及一個地面觀測孔，共5個觀測孔表示放水試驗期間的降深隨時間變化的動態(tài)趨勢（圖2）。位于第一階段放水孔群附近的觀測孔SK14-2 在放水的第一階段便有水頭響應(yīng)，降深隨著放水的進(jìn)行不斷增大；位于第二階段放水孔群附近的觀測孔F9-1，由于第一階段放水影響的含水層范圍有限，故在第一階段該處水頭響應(yīng)較弱，降深較小，第二階段初期降深迅速增大，之后穩(wěn)步上升；第三階段放水孔群附近的觀測孔F18的水頭響應(yīng)降深見圖2；第四階段放水孔群附近的觀測孔YZ10-1的水頭響應(yīng)降深見圖2。在整個放水試驗過程中，直羅組地面觀測孔J12、J17、SK6、SK8、SK11、SYJ的降深較小，均小于1 m，J17的降深見圖2。

圖2 降深隨時間變化的動態(tài)圖Fig.2 Representative temporal drawdown data

3.2 模型建立

概化數(shù)值模型為1 800 m × 3 200 m的二維非均質(zhì)飽和含水層，將含水層模型劃分為90 × 80個單元格，每個單元格為20 m × 40 m，四周為1 000 m的定水頭邊界，初始水頭為1 000 m。定水頭邊界及初始水頭設(shè)置為1 000 m可保證含水層始終維持飽和。水力層析反演計算時需給定水文地質(zhì)參數(shù)的初始場，設(shè)定lnT=1.3，方差為1.0，lnS=-5.0，方差為0.5，x和y方向的相關(guān)長度均為100 m。利用放水試驗獲得的每小時觀測孔水頭數(shù)據(jù)，迭代估計含水層水文地質(zhì)參數(shù)的空間分布特征。選取的18個時刻的水頭數(shù)據(jù)分別位于t=1，6，20，40，50，56，70，80，97，100，110，130，150，169，174，185，195，216 h，保留了水頭的整體趨勢和關(guān)鍵轉(zhuǎn)折點[23-24]。水力層析法的反演算法SimSLE 由亞利桑那大學(xué)Yeh 教授課題組開發(fā)的軟件VSAFT2 實現(xiàn)，下載鏈接為http://tian.hwr.arizona.edu/downloads。

4 結(jié)果與分析

4.1 水力層析反演結(jié)果分析

水文地質(zhì)參數(shù)的反演結(jié)果見圖3。圖3（a）、（b）是僅由第一階段的觀測孔水頭數(shù)據(jù)反演得到的導(dǎo)水系數(shù)和儲水系數(shù)的空間分布，一個靠近涌水位置的高T高S區(qū)初步顯現(xiàn)；圖3（c）、（d）是由第一、二階段的水頭數(shù)據(jù)聯(lián)合反演得到的參數(shù)空間分布；圖3（e）、（f）是由前三個階段的水頭數(shù)據(jù)聯(lián)合反演得到的參數(shù)空間分布；圖3（g）、（h）是由完整的四個階段的水頭數(shù)據(jù)聯(lián)合反演得到的參數(shù)空間分布。

隨著放水孔閥門的逐階段開啟，T和S的分布模式逐漸形成且輪廓清晰。工作面的南部涵蓋強導(dǎo)水區(qū)域，且有可能延伸至模型區(qū)域之外，與大流量、長時間的涌水密切相關(guān)。在巷道內(nèi)側(cè)，強導(dǎo)水區(qū)域向北延伸，繼而出現(xiàn)相對弱導(dǎo)水區(qū)域，并與弱儲水區(qū)域?qū)?yīng)。

實施第四階段放水試驗的目的是刻畫工作面北部甚至更遠(yuǎn)區(qū)域參數(shù)的空間分布。由于放水試驗影響范圍有限，且鉆孔集中分布于巷道附近，估計參數(shù)的空間異質(zhì)性僅體現(xiàn)在巷道內(nèi)側(cè)有限區(qū)域內(nèi)，巷道外側(cè)大部分區(qū)域基本按照參數(shù)均值均勻分布。

圖3中標(biāo)注了T=2 m2/d、T=8 m2/d、S=0.005和S=0.011時的等值線，等值線內(nèi)的區(qū)域分別對應(yīng)相對低T區(qū)、高T區(qū)、低S區(qū)和高S區(qū)。由圖3可推知，工作面南北富水性差異大，涌水區(qū)位于強導(dǎo)水帶上；從涌水區(qū)往北，導(dǎo)水系數(shù)和儲水系數(shù)整體逐漸減?。荒喜繉?dǎo)水系數(shù)和儲水系數(shù)均較大，屬于富水異常區(qū)。據(jù)文獻(xiàn)[25]所述，工作面南北富水性差異大與直羅河沖刷沉積有關(guān)，工作面位于沖刷帶的尖滅部位，因此南北富水性差異大。

4.2 水頭與水文地質(zhì)參數(shù)的互相關(guān)分析

為了直觀地將互相關(guān)系數(shù)的變化規(guī)律展現(xiàn)出來，在放水試驗的各階段分別選取一個觀測孔，計算了觀測孔在各階段初期、中期、末期的水頭與導(dǎo)水系數(shù)和儲水系數(shù)的互相關(guān)系數(shù)。

圖3 水文地質(zhì)參數(shù)的反演結(jié)果Fig.3 Inversion results of the hydrogeological parameters

圖4 水頭與導(dǎo)水系數(shù)的互相關(guān)系數(shù)ρhf時空分布Fig.4 Spatial and temporal distribution of cross-correlation between head and transmissivity

由圖4（a）、（b）、（c）可知，在放水試驗初期（t=10 h），直羅組地面觀測孔SYJ 與放水孔群之間存在明顯的負(fù)相關(guān)性，該負(fù)相關(guān)區(qū)域在整個放水階段持續(xù)存在，但隨著時間的延續(xù)，范圍逐漸縮小，負(fù)相關(guān)性逐漸減弱。在放水試驗中期（t=110 h），該觀測孔的上游（根據(jù)圖中流線所示方向判斷）及放水孔群的上游存在正相關(guān)性，隨著時間的延續(xù)，該區(qū)域范圍逐漸擴(kuò)展，最終連為一個整體，環(huán)繞在負(fù)相關(guān)區(qū)域周圍，正相關(guān)性逐漸加強。

由圖4（d）、（e）、（f）可知，在放水試驗的第一階段，井下觀測孔F9-1水頭未發(fā)生明顯變化，故互相關(guān)分析的初期選取t=80 h?；ハ嚓P(guān)系數(shù)空間分布的發(fā)展趨勢與觀測孔SYJ 類似，由初期的負(fù)相關(guān)發(fā)展為末期的正相關(guān)。

由圖4（g）、（h）、（i）可知，在第三階段放水試驗初期（t=100 h），井下觀測孔F14-2 下游及上游分別存在明顯的負(fù)相關(guān)區(qū)域及正相關(guān)區(qū)域，這兩個區(qū)域在第三、四階段持續(xù)存在，但隨著時間的延續(xù)，負(fù)相關(guān)區(qū)域范圍逐漸縮小，并逐漸減弱；正相關(guān)區(qū)域范圍逐漸擴(kuò)大，且逐漸增強?；ハ嚓P(guān)系數(shù)空間分布的發(fā)展趨勢與從第一階段開始的放水類似，但沒有顯現(xiàn)出放水孔群上游的正相關(guān)區(qū)域。在通常的水頭與導(dǎo)水系數(shù)的互相關(guān)分析中，一般針對一個觀測孔及一個放水孔，由于這里是放水孔群，觀測孔與放水孔之間的負(fù)相關(guān)區(qū)域及放水孔上游的正相關(guān)分布區(qū)域較大，最終導(dǎo)致正相關(guān)區(qū)域僅出現(xiàn)在觀測孔上游。

由圖4（j）、（k）、（l）可知，井下觀測孔W1 與觀測孔F14-2的互相關(guān)系數(shù)空間分布的發(fā)展趨勢類似。

由圖4可知，水頭與導(dǎo)水系數(shù)的互相關(guān)系數(shù)空間分布的發(fā)展趨勢類似，在初期負(fù)相關(guān)性較強，末期正相關(guān)性較強，且相關(guān)區(qū)域范圍擴(kuò)大。故放水試驗?zāi)┢诘乃^數(shù)據(jù)更利于估計導(dǎo)水系數(shù)，且不同位置觀測孔的相關(guān)區(qū)域不同，相近時間的互相關(guān)系數(shù)分布類似，因此可聯(lián)合不同觀測孔在不同時間的水頭數(shù)據(jù)，求取含水層導(dǎo)水系數(shù)的高精度空間分布。

對于相同觀測孔，在同一時間段的水頭與儲水系數(shù)的互相關(guān)系數(shù)矩陣的計算結(jié)果見圖5。水頭與儲水系數(shù)的互相關(guān)性規(guī)律相對簡單，整個放水過程均呈正相關(guān)性，且靠近觀測孔附近的相關(guān)性較強。隨著時間的延續(xù)，正相關(guān)性逐漸減弱。由圖5可知，放水試驗初期的水頭數(shù)據(jù)對于估計儲水系數(shù)比較有效。

圖4、圖5中的虛線為正、負(fù)相關(guān)性較強區(qū)域的輪廓線，直觀表征了該觀測孔的水頭數(shù)據(jù)可精細(xì)刻畫的參數(shù)空間范圍。對于同一觀測孔，圖中標(biāo)注了對應(yīng)時刻的同條流線，以判斷觀測孔和放水孔群各自的上、下游區(qū)域。

綜合水頭與導(dǎo)水系數(shù)和儲水系數(shù)的相關(guān)性分析結(jié)果，各個階段的相關(guān)系數(shù)分布具有相似性，正相關(guān)和負(fù)相關(guān)區(qū)域按照一定的規(guī)律變化。本研究針對每個觀測孔選擇的18個水頭數(shù)據(jù)，集中分布在每個階段的初期和末期，可以更加有效的獲取研究區(qū)的水文地質(zhì)參數(shù)分布信息。

圖5 水頭與儲水系數(shù)的互相關(guān)系數(shù)ρhs時空分布Fig.5 Spatial and temporal distribution of cross-correlation between head and storativity

5 結(jié)論

通過本研究可見，利用水力層析法分析井下疊加放水試驗，是在傳統(tǒng)放水試驗分析基礎(chǔ)上的創(chuàng)新，根據(jù)本次井下疊加放水試驗，可以得出以下結(jié)論：

（1）水力層析法是一種有效的非均質(zhì)含水層參數(shù)識別的新技術(shù)，相較于傳統(tǒng)放水試驗分析方法，它在識別含水層非均質(zhì)特性方面具有優(yōu)勢。

（2）互相關(guān)分析結(jié)果表明每個觀測孔數(shù)據(jù)與水文地質(zhì)參數(shù)都存在相關(guān)性較強的區(qū)域，且初期及末期數(shù)據(jù)分別與儲水系數(shù)和導(dǎo)水系數(shù)相關(guān)性較強，可增加觀測數(shù)據(jù)，中期數(shù)據(jù)適當(dāng)稀疏選取，將多個觀測孔數(shù)據(jù)聯(lián)合起來結(jié)合水力層析技術(shù)，可以較準(zhǔn)確地反演出符合實際情況的參數(shù)分布信息。

（3）本研究將井下放水視為針對含水層的刺激源，聯(lián)合多個觀測孔的水頭響應(yīng)數(shù)據(jù)，反演研究區(qū)域水文地質(zhì)參數(shù)的空間分布，捕捉到了研究區(qū)域的主要地質(zhì)結(jié)構(gòu)特征，刻畫較精細(xì)。針對有突水風(fēng)險的礦井，預(yù)先采用水力層析法分析礦井水文地質(zhì)條件，推采時避開富水異常區(qū)，或者提前采取有針對性的防治水工程措施，可以有效降低突水事故風(fēng)險。