王仕昌,邱愷毅

(山東省物化探勘查院,山東 濟(jì)南 250013)

濟(jì)寧鐵礦位于山東省濟(jì)寧市兗州區(qū)、任城區(qū)、汶上縣交界地帶，面積約46km2。濟(jì)寧鐵礦賦存于新太古代濟(jì)寧群變質(zhì)巖中，埋深多在1000m以下，經(jīng)勘查查明鐵礦石資源量約18.34億t[1]。因埋藏深，開(kāi)采難度大，為開(kāi)采可行性研究提供水文地質(zhì)數(shù)據(jù)，施工了3個(gè)水文地質(zhì)孔，進(jìn)行了抽水試驗(yàn)，求取了水文地質(zhì)參數(shù)，采用數(shù)值模擬法[2-3]，對(duì)礦坑涌水量進(jìn)行了預(yù)測(cè)[4-5]。

1 礦區(qū)水文地質(zhì)條件

礦床頂板為寒武紀(jì)饅頭組地層，礦床周邊為濟(jì)寧群變質(zhì)巖地層，巖石堅(jiān)硬、完整，富水性、導(dǎo)水性較差，礦床底部仍為濟(jì)寧群變質(zhì)巖地層[6-8]，隨著深度增加富水性及導(dǎo)水性趨于微弱，礦床底部可視為相對(duì)隔水層。

2 數(shù)值模擬預(yù)測(cè)礦坑涌水量

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

2.1.1 模型模擬范圍

模型平面范圍選擇依據(jù)①濟(jì)寧礦區(qū)屬于湖東汶泗河沖洪積平原中部水文地質(zhì)單元；②鐵礦資源儲(chǔ)量范圍遠(yuǎn)小于礦區(qū)范圍，可以保證模擬范圍在抽水試驗(yàn)時(shí)影響不到模擬邊界；③能將區(qū)域內(nèi)切割深入礦區(qū)的F1斷裂較完整地包含于模型模擬范圍內(nèi)，方便模擬斷層對(duì)抽水試驗(yàn)的影響[9]。

模型垂直范圍模擬含礦層位濟(jì)寧群變質(zhì)巖裂隙含水層，模型上頂界為YK1第3次抽水試驗(yàn)起點(diǎn)-1113m，下底界定為YK1第6次抽水試驗(yàn)底界下移100m，即-1613m處。

2.1.2 地下水動(dòng)力場(chǎng)及流動(dòng)特征

濟(jì)寧群變質(zhì)巖裂隙水補(bǔ)給、徑流、排泄條件差，根據(jù)鉆孔鉆探情況，地下水賦存于千枚巖破碎帶中。千枚巖結(jié)構(gòu)密實(shí)，孔隙、裂隙不發(fā)育[10-11]。根據(jù)抽水試驗(yàn)，斷層上盤(pán)灰?guī)r、白云巖裂隙巖溶含水層的地下水通過(guò)F1導(dǎo)水?dāng)鄬恿飨驖?jì)寧群變質(zhì)巖裂隙含水層。地下水運(yùn)動(dòng)為三維流；地下水水位是隨時(shí)間變化的，由于介質(zhì)的非均勻性造成水文地質(zhì)參數(shù)隨空間變化，體現(xiàn)了系統(tǒng)的非均質(zhì)性。初始流場(chǎng)[12]根據(jù)抽水試驗(yàn)的靜止水位給出，概化模擬含水層天然狀態(tài)下水量保持穩(wěn)定。

2.1.3 水文地質(zhì)參數(shù)及其靈敏度分析

該次研究對(duì)象為深層的濟(jì)寧群變質(zhì)巖裂隙承壓含水層，含水層的水文地質(zhì)參數(shù)主要為滲透系數(shù)[13-14]。根據(jù)YK1鉆孔四次抽水試驗(yàn)擬合求得的濟(jì)寧群變質(zhì)巖含水層的滲透系數(shù)分別為0.0055m/d，0.0083m/d，0.000021m/d，0.0046m/d。

水文地質(zhì)參數(shù)的不確定性[15]往往導(dǎo)致計(jì)算水頭或流量的不確定性，影響礦區(qū)涌水量的可靠性。因此應(yīng)用靈敏度分析，對(duì)識(shí)別驗(yàn)證后的水文地質(zhì)參數(shù)進(jìn)行不確定性分析[8]，為研究礦區(qū)涌水量提供依據(jù)。

靈敏度可度量一種因子的變化對(duì)另一種因子的影響程度，某一模型因變量對(duì)模型輸入?yún)?shù)的靈敏度可用因變量對(duì)該輸入?yún)?shù)的偏導(dǎo)數(shù)來(lái)表示,即:

Xi,k=?yi/?ak

式中:Xi,k—模型的因變量對(duì)第i觀測(cè)點(diǎn),第k參數(shù)的靈敏度系數(shù)。由于不同的參數(shù)單位不同，這樣不同參數(shù)的靈敏度系數(shù)就無(wú)可比性。此時(shí),需對(duì)上式進(jìn)行標(biāo)準(zhǔn)化處理,即:

根據(jù)上式計(jì)算得出靈敏度系數(shù)(表1)，從表1可以看出，滲透系數(shù)的變化對(duì)模型的運(yùn)行結(jié)果有一定影響。

表1 靈敏度系數(shù)

2.2 數(shù)值模型

2.2.1 數(shù)學(xué)模型

對(duì)于建立的礦區(qū)水文地質(zhì)概念模型，可用如下微分方程的定解問(wèn)題來(lái)描述：

式中：Ω—滲流區(qū)域；h—地下水系統(tǒng)的水位標(biāo)高(m)；K—含水介質(zhì)的水平滲透系數(shù)(m/d)；Kz—含水介質(zhì)垂向滲透系數(shù)(m/d)；ε—含水層的源匯項(xiàng)(1/d)；Ss—含水層儲(chǔ)水率(1/m)。

2.2.2 子程序包的選擇

運(yùn)用基于有限差分法的GMS程序求解以上微分方程的定解問(wèn)題。根據(jù)礦區(qū)的水文地質(zhì)結(jié)構(gòu)等特征，選擇GMS中相應(yīng)的子程序包來(lái)實(shí)現(xiàn)地下水流的模擬。

(1)LPF-層狀特征流量程序包。該功能相當(dāng)于計(jì)算單元間滲流子程序包，它對(duì)各個(gè)模型參數(shù)作出了相應(yīng)的假定，選定該程序包可以通過(guò)參數(shù)分區(qū)方式賦值參數(shù)。

(2)WEL1-井流子程序包(WEL)。模型約定，在每個(gè)應(yīng)力期，以指定流量從含水層抽水或向含水層注水；井流量不受井所在計(jì)算單元的大小及水頭影響；負(fù)的流量值表示抽水井，而正的流量值則表示注水井。

(3)GHB-通用水頭子程序包。將流向通用水頭邊界的流量項(xiàng)加進(jìn)有限差分方程組。

2.2.3 數(shù)值模型

(1)空間離散[16]

將模擬區(qū)剖分為104行，81列，共計(jì)4層，各層均采用100m×100m的剖分格式。4層有效單元格共33696個(gè)。

(2)模擬期確定

該次地下水?dāng)?shù)值模擬模型的模擬期為YK1第一次抽水試驗(yàn)開(kāi)始至最后一次抽水試驗(yàn)結(jié)束水位恢復(fù)。在抽水期和水位恢復(fù)期，每小時(shí)作為一個(gè)應(yīng)力期，分層抽水的間隔時(shí)期，時(shí)間節(jié)點(diǎn)以小時(shí)的倍數(shù)逐步過(guò)渡增長(zhǎng)：2h,4h,8h……共計(jì)233個(gè)應(yīng)力期。

(3)定解條件的處理

初始條件：由于缺乏地下水流場(chǎng)，加之模擬含水層為深層承壓水，循環(huán)條件差，遂根據(jù)分層抽水試驗(yàn)的靜止水位分別給出統(tǒng)一水位作為每層初始流場(chǎng)。

邊界條件：全部人為邊界處理為通用水頭邊界。

2.2.4 模型擬合情況

根據(jù)YK1鉆孔的第3～6次分層抽水試驗(yàn)的情況，擬合結(jié)果如圖1～圖4所示。

圖1 YK1第3次抽水試驗(yàn)擬合圖

圖2 YK1第4次抽水試驗(yàn)擬合圖

圖3 YK1第5次抽水試驗(yàn)擬合圖

圖4 YK1第6次抽水試驗(yàn)擬合圖

根據(jù)擬合圖看出，YK1鉆孔的第3～5次抽水試驗(yàn)的擬合效果均較好，第6次抽水試驗(yàn)的擬合有所偏差，分析原因是將模擬的含水層底板放低100m所致。

2.3 礦井涌水量預(yù)測(cè)

2.3.1 根據(jù)等效滲透系數(shù)原理合并模型為單層

礦層含水層為濟(jì)寧群變質(zhì)巖裂隙含水層1層，根據(jù)等效滲透系數(shù)原理將原本垂向劃分為4層的模型合并為1層。

滲流方向與層面水平時(shí)，設(shè)其平行于層面的滲透系數(shù)為KP，當(dāng)單一分層的滲透系數(shù)Ki和厚度Mi已知時(shí)，可求出KP。當(dāng)水流平行于層面時(shí)，通過(guò)層狀含水層總的單寬流量q等于各分層的單寬流量之和，總厚度M等于各分層厚度之和。對(duì)于每一分層而言，水力坡度J均為ΔH/l，因此，每一分層的流量為：

如果用等效的均質(zhì)含水層代替層狀巖層，顯然等效層的厚度等于層狀巖層的總厚度，并且在同一水力坡度ΔH/l作用下應(yīng)當(dāng)有相同的流量q：

由此得：

因而求得平行于層面方向的等效滲透系數(shù)為：

濟(jì)寧群變質(zhì)巖裂隙含水層的最終水文地質(zhì)參數(shù)分區(qū)如圖5所示。

圖5 滲透系數(shù)分區(qū)圖

2.3.2 涌水量預(yù)測(cè)

(1)疏干排水量

一般情況下，對(duì)于給定的疏干降深，可以分別求出不同的疏干流量所對(duì)應(yīng)的疏干時(shí)間，然后作出不同疏干流量的降深與時(shí)間關(guān)系曲線(圖6)和指定疏干水平的流量與時(shí)間關(guān)系曲線(圖7～圖10)，據(jù)此選擇最佳的疏水量與疏干時(shí)間[17-19]。

圖6 不同疏干流量的降深與時(shí)間關(guān)系曲線圖

圖7 開(kāi)采水平-1200m的Q-t曲線

圖8 開(kāi)采水平-1300m的Q-t曲線

圖9 開(kāi)采水平-1400m的Q-t曲線

圖10 開(kāi)采水平-1500m的Q-t曲線

由圖7～圖10可以看出，開(kāi)采水平為-1200m，-1300m，-1400m，-1500m所對(duì)應(yīng)的最佳疏干時(shí)間為20.00h,25.23h,31.70h,39.73h；疏干排水量為3250m3/d。

(2)正常涌水量

將疏干工作面設(shè)置為定水頭，用穩(wěn)定流場(chǎng)求出礦井涌水量(表2)。

表2 正常涌水量預(yù)測(cè)

將模型各參數(shù)分別增減10%后進(jìn)行穩(wěn)定流模型正常涌水量的預(yù)測(cè)，可以得到正常涌水量的1個(gè)區(qū)間值(表3)。

表3 正常涌水量區(qū)間值

(3)最大涌水量

最大涌水量的預(yù)測(cè)考慮到礦體埋藏深度大，結(jié)合該次勘查階段的精度低于C級(jí)、高于D級(jí)，而且未做群孔抽水試驗(yàn)，根據(jù)《地下水資源分類分級(jí)標(biāo)準(zhǔn)》(GB15218-94)附錄C[20]，C級(jí)計(jì)算的允許誤差為±35%， D級(jí)計(jì)算的允許誤差為±50%，因此為提高安全程度，按D級(jí)允許誤差±50%計(jì)算，最大涌水量為正常涌水量的1.5倍計(jì)算，結(jié)果見(jiàn)表4。

表4 礦井最大涌水量

3 結(jié)語(yǔ)

濟(jì)寧鐵礦埋藏于濟(jì)寧群變質(zhì)巖中，埋藏深度在1000m以下，該次通過(guò)數(shù)值模擬預(yù)測(cè)了-1200m,-1300m,-1400m,-1500m開(kāi)采水平的正常礦坑涌水量分別為998.56m3/d,1079.27m3/d，1159.98m3/d,1240.68m3/d，但由于濟(jì)寧鐵礦賦存條件復(fù)雜，礦坑涌水量的預(yù)測(cè)較難，其預(yù)測(cè)數(shù)值與實(shí)際數(shù)值的一致性有待于進(jìn)一步探討研究。