臥龍湖礦太灰含水層水文地質(zhì)求參

王 濤

中煤第三建設(shè)集團安徽地質(zhì)測量工程有限責任公司，安徽宿州，234000

臥龍湖礦太灰含水層水文地質(zhì)求參

王 濤

中煤第三建設(shè)集團安徽地質(zhì)測量工程有限責任公司，安徽宿州，234000

通過對臥龍湖煤礦101和102采區(qū)進行放水試驗，得到了放水第一、二階段主要太灰觀測孔降深與流量歷時曲線以及水位恢復階段觀測孔降深與流量的歷時曲線。根據(jù)放水試驗結(jié)果與采區(qū)水文地質(zhì)條件，建立水文地質(zhì)概念模型并進行水文地質(zhì)單元劃分，采用FEFLOW模擬軟件進行模型的建立、識別、對比和驗證，求得了較為合理的各水文地質(zhì)單元地下水動力學參數(shù)，進而揭示出了礦區(qū)內(nèi)101與102采區(qū)兩個太灰含水層的高水位分布區(qū)。

放水試驗；太灰含水層；數(shù)值模擬；臥龍湖礦

臥龍湖煤礦深部煤層開采受到下伏太原群巖溶含水層(簡稱“太灰”)和奧陶系巖溶含水層(簡稱“奧灰”)的嚴重威脅。為了查明臥龍湖礦101與102采區(qū)太灰含水層的富水性、補給來源、補給方式及可能的補給地段，進而得出采區(qū)灰?guī)r水流場分布特征，開展井下放水試驗是一種有效手段。FEFLOW是由德國WASY水資源規(guī)劃系統(tǒng)研究所研制開發(fā)的基于有限單元法的地下水模型軟件包，是迄今為止功能最為齊全的地下水模擬軟件包之一。該軟件提供圖形人機對話功能、具備地理信息系統(tǒng)數(shù)據(jù)接口、能夠自動產(chǎn)生空間各種有限元網(wǎng)、具有空間參數(shù)區(qū)域化、快速精確的數(shù)值算法和先進的圖形視覺化技術(shù)等特點[1]。將放水試驗與FEFLOW模擬同時應(yīng)用于采區(qū)水文地質(zhì)條件分析,是煤礦防治水研究工作中的有益嘗試。

1 研究區(qū)水文地質(zhì)概況

101與102采區(qū)位于臥龍湖井田西南部，采區(qū)西部及南部以10煤層露頭為界，東部以F9斷層為界，北部以F8斷層為界。101與102采區(qū)上限為煤層露頭線，煤層底板標高-310 m左右，下限煤層底板標高為-550 m，采區(qū)整體形態(tài)呈梯形，北寬南窄，南北長約5000 m，東西寬600～1400 m，面積約5 km2。兩采區(qū)內(nèi)地表無基巖出露，沖積層厚330～355 m，由西南至東北逐漸變厚，區(qū)內(nèi)地勢平坦。

四含水因與基巖直接接觸，構(gòu)成四含水與煤系含水層水、灰?guī)r含水層水成為互補關(guān)系，煤系含水層與下伏灰?guī)r含水層之間有相對隔水層存在，本身聯(lián)系不密切，但采區(qū)由于F8、F9及BF34等落差大于50米斷層的存在，造成灰?guī)r與煤系地層直接接觸，煤系含水層與灰?guī)r含水層為互補關(guān)系[2]。

2 研究區(qū)放水試驗

此次研究的主要對象為太原組上段灰?guī)r水。由于其為高承壓含水層，井下也不具備足夠的排水能力和完備的測流、測壓裝備，故本次試驗按照非穩(wěn)定流方法，采用井下放水，井上、井下同步觀測的方法。為充分暴露礦井水文地質(zhì)條件，首先102-觀3開始放水，其他鉆孔作為觀測孔。持續(xù)放水一段時間后，增加102-放5為放水孔。放水試驗歷時15天，第1d12:00之前開始原始水位觀測，12:00時后打開101-觀3開始放水，由于效果不明顯，于第4d11:00打開102-放5;第10d18:00關(guān)閉兩個放水孔，開始水位恢復。放水第一、二階段主要太灰觀測孔降深與流量歷時曲線如圖1所示，水位恢復階段太灰觀測孔降深歷時曲線如圖2所示。

圖1 放水第一、二階段太灰觀測孔降深與流量歷時曲線

圖2 水位恢復階段太灰觀測孔降深歷時曲線

3 水文地質(zhì)概念模型與數(shù)學模型

3.1 水文地質(zhì)概念模型

根據(jù)臥龍湖煤礦水文地質(zhì)條件，結(jié)合太灰101-觀2的觀測結(jié)果，可以初步推斷F9斷層以西的塊段為一相對獨立的水文地質(zhì)單元。東側(cè)以落差0～70 m的F9斷層為區(qū)外相隔，北側(cè)以落差0～90 m的F8斷層為區(qū)外相隔，南側(cè)與西側(cè)以井田淺部10煤層露頭區(qū)為界。采區(qū)主要水害威脅水源太灰含水層包括4個薄層灰?guī)r，各含水層間無穩(wěn)定隔水層，水體聯(lián)系密切，將其視為一非均質(zhì)、各向異性統(tǒng)一含水體。該區(qū)域的北部邊界定義為流量邊界，邊界的具體流量根據(jù)北部太灰2001-觀1在放水過程中降深與流場分布規(guī)律進行分段估計。F9斷層以南，101采區(qū)東翼軌道巷的太灰101-觀1在整個放水試驗過程中壓力不變，表明F9斷層阻隔太灰水流動，F(xiàn)9斷層為隔水斷層。研究區(qū)域的F8斷層落差較大，在計算區(qū)域內(nèi)延長不大，定義為隔水邊界。在F8與F9斷層之間的邊界無構(gòu)造分布，直接與外界相同，定義為

圖3 101與102采區(qū)太灰水文地質(zhì)概念模型示意圖

流量邊界。區(qū)內(nèi)太灰水主要通過兩個放水孔及102采區(qū)局部泄水點等通道進行排泄。

根據(jù)整個101與102采區(qū)放水試驗結(jié)果，結(jié)合兩采區(qū)水文地質(zhì)條件與構(gòu)造條件，將計算區(qū)域進行了水文地質(zhì)參數(shù)分區(qū)：一共劃分出6個水文地質(zhì)參數(shù)區(qū)，并在局部復雜區(qū)域定義了緩沖區(qū)域(圖3)。

3.2 數(shù)學模型

數(shù)學模型是水文地質(zhì)概念模型的數(shù)學語言表達形式，它是根據(jù)地下水動力學原理，由一個泛定方程加定解條件(邊界條件和初始條件)組成。根據(jù)上述水文地質(zhì)概念模型，即得二維非均質(zhì)、各向異性、非穩(wěn)定承壓地下水流的數(shù)學模型，可用下列方程式來描述，即：

H(x,y,t)=H1(x,y)∈Γ1

=-q(x,y,t) (x,y)∈Γ2,t>0

4 數(shù)值模型的建立、識別與對比驗證

4.1 數(shù)值模型的建立

圖4 研究區(qū)有限元計算剖分網(wǎng)格

通過對研究區(qū)內(nèi)所有的水文資料進行整理，將研究區(qū)進行參數(shù)分區(qū)[5-6]。采用三角形單元將兩個研究區(qū)進行三角剖分，將101與102研究區(qū)分為2238個單元，2360個節(jié)點(圖4)。

4.2 數(shù)值模型識別與驗證

根據(jù)已掌握的礦井水文地質(zhì)條件，將模型所需全部數(shù)據(jù)輸入模擬計算程序，通過不斷調(diào)整模型參數(shù)，修正模型的各部分，完成數(shù)值模型對客觀水文地質(zhì)原型較為逼真的擬合。擬合過程采用人工試算與機器調(diào)試相結(jié)合的方法，即先人為給定未知量(參數(shù)或水量)進行正向計算，求得目標函數(shù)，并不斷地修改未知量，重復進行正演計算，直至求得的目標函數(shù)滿足誤差為止，此時的未知量即是所要求的參數(shù)或水量。再給定未知量的約束條件，讓機器自動尋優(yōu)，不斷地解正問題。求得目標函數(shù)達極小值的未知量，即得到所求的參數(shù)。然后結(jié)合礦井的水文地質(zhì)條件，進行綜合分析，確定較合理的一組。

根據(jù)放水試驗，本次模擬時間段為第4d15：00～第10d18：00～第15d9:00兩個時間段。從圖5中第10d18：00的時刻水位分布模擬與實測對比圖可以看出，模擬結(jié)果基本與觀測結(jié)果相符合，從圖中可以明顯地看出兩個高水位分布區(qū)與兩個降落漏斗的分布區(qū)。從圖6中第15d9:00的時刻水位分布模擬與實測對比圖可以看出，水位恢復速度較慢，高水位分布區(qū)水位變化很小，少數(shù)觀測孔水位在此時間段內(nèi)恢復到初始水位。盡管部分觀測孔上水位恢復模擬不太理想，但大體上水位恢復階段的模擬還是符合放水試驗反映的客觀實際。

圖5 第10d 18：00時刻水位分布模擬與實測等值線對比圖

太灰水文觀測孔水位變化曲線如圖7所示，從擬合曲線上可看出，放水試驗中兩水位觀測孔模擬水位與觀測水位基本一致，變化趨勢大體相符。

4.3 水文地質(zhì)參數(shù)分區(qū)結(jié)果

101與102采區(qū)水文地質(zhì)參數(shù)識別結(jié)果如表1所示，單元分區(qū)內(nèi)的水文地質(zhì)參數(shù)基本與采用s-lg(t/r2)直線圖解法一致。

表1 101與102采區(qū)水文地質(zhì)參數(shù)識別結(jié)果

5 結(jié) 論

根據(jù)放水試驗，結(jié)合FEFLOW數(shù)值模擬，建立了臥龍湖煤礦101與102采區(qū)太灰含水層水文地質(zhì)概念模型，合理確定了研究區(qū)水文地質(zhì)單元，求得各水文地質(zhì)單元地下水動力學參數(shù)，揭示了101與102采區(qū)兩個太灰含水層高水位分布區(qū)。數(shù)值模擬結(jié)果與實測結(jié)果基本一致。

[1]張國偉,孫亞軍,徐智敏,等.FEFLOW在區(qū)域地下水流場演化數(shù)值模擬中的應(yīng)用[J].地下水,2010,32(1):3-6

[2]桂和榮,陳陸望.礦區(qū)地下水水文地球化學演化與識別[M].北京:地質(zhì)出版社,2007:57-61

[3]任印國,柳華武,李明良.石家莊市東部平原FeFLow地下水數(shù)值模擬與研究[J].水文,2009,29(5):59-62

[4]冶雪艷,杜新強.FeFlow在紅興水庫滲控方案確定中的應(yīng)用[J].人民黃河,2009,31(9):116-118

[5]林壢,楊峰,崔亞莉等.FeFLow在模擬大區(qū)域地下水流中的特點[J].北京水務(wù),2007,(1):43-46

[6]李守波,趙傳燕,馮兆東.黑河下游地下水波動帶地下水時空分布模擬研究：FeFLow模型應(yīng)用[J].干旱區(qū)地理,2009,32(3):391-395

(責任編輯：汪材印)

10.3969/j.issn.1673-2006.2016.12.030

2016-06-12

王濤(1987-)，安徽宿州人，主要研究方向：工程地質(zhì)與水文地質(zhì)。

P641

A

1673-2006(2016)12-0108-04