2條平行隔水邊界干擾井礦坑涌水量預算方法
——以道真縣洛龍鋁土礦為例

侯江勇，羅陽喬，黃海韻

（貴州省地質礦產勘查開發局117地質大隊，貴州貴陽550018）

2條平行隔水邊界干擾井礦坑涌水量預算方法
——以道真縣洛龍鋁土礦為例

侯江勇*，羅陽喬，黃海韻

（貴州省地質礦產勘查開發局117地質大隊，貴州貴陽550018）

通過對道真縣洛龍鋁土礦區水文地質條件進行概述、分析，劃分了含、隔水巖組，明確了礦床主要充水含水層為二疊系中統棲霞、茅口組（P2q+m）。根據礦區水文地質條件，建立了礦區水文地質概念模型，通過概念模型，建立了2條平行隔水邊界穩定流干擾井計算數學模型，最后將各參數帶入數學模型預算出礦坑涌水量。

涌水量預算方法；2條平行隔水邊界；干擾井

礦區涌水量預算精確度的高低，直接影響到未來礦山的安全度及經濟效益，一直以來，從事礦區水文地質工作的工作者對礦區涌水量預算較為頭疼，預算出的礦坑涌水量與礦山開采后實際涌水量偏差相當大，造成這種偏差實際上是預算時選則的預算方法不妥。采用大井法預算礦坑涌水量時，有部分礦區，在工具書、教課書上往往沒有適合該礦區條件現成的公式，只能自己推導。本人在預算道真縣洛龍鋁土礦礦坑涌水量時，在工具書、教課書、網上都無法找到適合該礦區邊界條件的預算公式，即2條平行隔水邊界的“穩定流潛水完整干擾井”涌水量預算公式。本人采用2條平行隔水邊界單井涌水量預算公式，通過映射原理、疊加原理等推導出了2條平行隔水邊界的“穩定流潛水完整干擾井”涌水量預算公式，解決了礦坑涌水量預算問題。

1　水文地質概況

從區域上看，礦區位于芙蓉江流域大塘向斜水文地質單元中段。礦區總體地勢北高南低，東西高、中間低，主干河流洛龍河從礦區中部穿過，洛龍河為當地侵蝕基準面和集中排泄場所。出露的含水巖組有：三疊系下統夜郎組（T1y）、二疊系上統吳家坪組（P3w）、二疊系上統梁山—茅口組（P2l+q+m），出露的相對隔水巖組有：志留系下統韓家店群（S1Hj）、志留系下統石榴欄、龍馬溪組（S1l+sh）。區內構造中等發育，發育一大塘向斜匯水構造，F1斷裂構造，構造線總體走向近南北向。見圖1。

礦區含礦層為二疊系下統大竹園組（P1d），礦層直接頂板為梁山—茅口組（P2l+q+m），間接頂板為吳家坪組（P3w）、夜郎組（T1y），頂板為碳酸鹽巖，遍布于整個礦區。直接頂板梁山—茅口組巖性主要為灰巖、泥灰巖，底部為0.18～2.71m厚的炭質粘土巖，該層總厚395.05～484.60m，巖溶管道較發育，含裂隙溶洞水及暗河水，富水性強，泉點流量為0.26～385.90L/s，在水文孔中對棲霞、茅口做混合抽水試驗，其單位涌水量平均值1.74L/（s·m），平均滲透系數（K）為0.396m/d，為礦坑直接充水含水層；間接頂板吳家坪組（P3w）、夜郎組（T1y）底部為2～5m厚的碎屑巖相對穩定，起到一定的隔水作用，且離含礦層較遠，對未來礦井充水可能性較小。直接底板主要為韓家店群（S1Hj），巖性為碎屑巖，厚度大于500m，含水性弱，透水性較差，為較穩定的相對隔水層。

2　礦坑涌水量預算

2.1礦坑充水因素

未來礦坑的主要充水水源地表水主要為洛龍河河水，地下水主要為棲霞、茅口巖溶水；充水通道主要為溶洞、暗河、溶蝕裂隙等巖溶管道，頂板采動裂隙以及斷層；充水方式主要為頂板直接充水。

2.2礦坑涌水量預算

根據礦區水文地質條件，采用大井法預算未來礦井的正常涌水量和最大涌水量。經勘查，所圈定的礦體為南北2個不相連的礦體，即Ⅰ號和Ⅱ號礦體，本次預算的礦坑涌水量，Ⅰ號礦體預算到資源量最低標高+280m水平，Ⅱ號礦體預算到資源量最低標高+430m水平。

圖1　水文地質略圖

2.2.1礦區水文地質概念模型

（1）水文地質邊界條件：

底板邊界：以韓家店群（S1Hj）隔水底板為界，礦層上覆直接巖溶充水層棲霞、茅口組（P2q+m）。

側向邊界：東西兩側以韓家店群（S1Hj）與棲霞組（P2q）地層界限為側向隔水邊界。南北為無限邊界。

（2）含水介質特征：主要充水含水層為棲霞、茅口組（P2q+m）巖溶含水巖組，其含水介質為構造裂隙、層間裂隙、溶蝕裂隙組合巖型，含水性不均勻，其含水性可近視為各項同性。

（3）地下水動力類型：天然條件下地下水動力特征為潛水。

（4）地下水徑流場：礦區地下水的運動為有限含水層中潛水自然運動。但在開采條件下，其流速、運動方向、水位都會改變，但其運動仍為潛水運動。

（5）概念模型的確定：由于大井位于東西2條平行隔水邊界內，未來開采條件下，2口大井將形成相互干擾，地下水的運動形成2平行隔水邊界附近地下水向積水構筑物（礦井）的無壓運動。通過測量和計算（見圖2），礦坑離東西2條隔水邊界線的距離小于影響半徑的一半，根據邊界條件的概化，選擇平行帶狀2條隔水邊界的“穩定流潛水完整干擾井”公式預算礦坑涌水量，見圖2。

圖2　概念模型圖

2.2.22條平行隔水邊界穩定流干擾井計算數學模型

由于礦段內充水含水層厚度大，而排水水位降深與含水層厚相差不大，因此，為簡便計算，可將含水系統近似地視為潛水水系統來處理。

對于潛水含水層降深（S）不是線性函數，不能進行簡單的疊加，故有：

Δh2=H02-h2，令u=H02-h2，使u成為線性函數后再進行疊加。采用2條平行隔水邊界單井涌水量預算公式進行疊加

同理：

式中:R1、R2——Ⅰ、Ⅱ號礦體未來開采井影響半徑；

L——2條平行隔水邊界之間的垂直距離；

b——大井距西邊隔水邊界的距離；

r1、r7——Ⅰ、Ⅱ號礦體未來開采大井的半徑；

x——兩大井之間的距離；

Q1、Q2——Ⅰ、Ⅱ號礦體未來開采大井涌水量。

2.2.3計算參數取值

（1）滲透系數（K）：取鉆孔SZK52442抽水試驗成果平均數0.396m/d。

（2）潛水含水層厚度（H）：由23個鉆孔水位平均標高（872.81m）與礦段內含水層最低標高之差確定，礦段內含水層最低標高為280m，由此可得含水層厚度H為592.81m。

圖3　礦坑涌水量預算平面圖

（3）排水降深（S）：由23個鉆孔水位平均標高（872.81m）與礦體最低開采標高確定，Ⅰ、Ⅱ號礦體最低開采標高分別為+280m，+430m，由此得Ⅰ號礦體排水降深為592.81m、Ⅱ號礦體排水降深為442.81m。

（4）大井半徑（r）：將Ⅰ號礦體概化成不規則的圓形，由公式確定，公式中的F為未來開采區的面積；Ⅱ號礦體概化成橢圓形，由公式r7=d1+d2/4確定。d1、d2為橢圓形長軸、短軸的長，在礦坑涌水量圖（圖3）上量取。

Ⅱ號礦體：r7=d1+d2/4=1012m

（5）影響半徑（R）：

（6）引用影響半徑（R0）：

Ⅰ號礦體：R0=R+r0=15041m

Ⅱ號礦體：R0=R+r0=7797m

（7）大井中心點距西邊隔水直線邊界的垂直距離（b）：在礦坑涌水量圖（圖3）上量取，經測量得b= 3622m。

（8）2條隔水直線邊界的垂直距離（L）：在礦坑涌水量圖（圖3）上量取得L=9298m。

（9）2大井之間的距離（x）：在礦坑涌水量圖（圖3）上量取得x=3750m。

2.2.4礦井涌水量預算結果

根據建立的數值模型和水文地質參數，以Ⅰ號礦體最低開采標高+280m水平之上，Ⅱ號礦體最低開采標高+430m水平之上預算礦坑涌水量，將其各參數代入公式（1）、（2）聯立求解關于Q1、Q2的二元一次方程組，其計算結果見表1。

表1　礦井涌水量預算結果統計表

3　結論

以道真縣洛龍鋁土礦為例，根據水文地質特征、邊界條件等建立水文地質概念模型，運用疊加原理建立了數學模型，最綜推導出了2條平行隔水邊界的“穩定流潛水完整干擾井”預算礦坑涌水量的公式，解決了礦坑涌水量預算問題。

［1］侯江勇，向通，等.貴州省道真縣洛龍鋁土礦詳查報告［R］.貴州省地質礦產勘查開發局117地質大隊，2014.

［2］薛群禹，朱學愚，等.地下水動力學［M］.2版.北京：地質出版社，1997.

［3］中國地質調查局.水文地質手冊［M］.北京：地質出版社，2012.

TV131

A

1004-5716（2016）06-0123-05

2016-02-26

2016-02-26

侯江勇（1982-），男（土家族），貴州江口人，工程師，現從事礦區水工環地質災害調查及勘查等工作。