（河北地質(zhì)大學(xué)，河北 石家莊 050031）

邢臺(tái)市高新技術(shù)產(chǎn)業(yè)開(kāi)發(fā)區(qū)（以下簡(jiǎn)稱高開(kāi)區(qū)）位于河北省南部，太行山東麓，邢臺(tái)市區(qū)東南部。隨著經(jīng)濟(jì)的發(fā)展，礦山能源開(kāi)采與環(huán)境發(fā)展之間的矛盾日趨尖銳。正確分析井下水補(bǔ)給量和排泄量，為能源的可持續(xù)利用及科學(xué)管理提供依據(jù)。

國(guó)內(nèi)外學(xué)者對(duì)邢臺(tái)地區(qū)礦井水文均衡情況進(jìn)行了大量的研究。邵愛(ài)軍等(1996)[2]用水量均衡法確定了騰發(fā)量，根據(jù)土壤水動(dòng)力學(xué)原理，推導(dǎo)了某些礦山難以測(cè)定的均衡要素的計(jì)算公式；許勝寶等(2005)[3]對(duì)湘西礦床富集區(qū)地下水進(jìn)行了均衡分析,并構(gòu)建了均衡模型,對(duì)井下水的補(bǔ)給進(jìn)行了評(píng)價(jià)；陳引峰等(2004)[4]運(yùn)用地下水動(dòng)態(tài)資料,對(duì)欒城縣井下水均衡進(jìn)行了初步分析和研究；Richard(1956)[5]首次提出零通量面法，其根本是包氣帶中的均衡量的測(cè)定，隨后便應(yīng)用到大量的研究中。Cooper(1990)[6]應(yīng)用零通量面法對(duì)英格蘭白堊系礦層井下水進(jìn)行了補(bǔ)給量的研究，證明年均補(bǔ)給量占降水量的27%。此外，水位動(dòng)態(tài)法適用于潛水含水層，Healy and Cook（2002）[7]認(rèn)為只要確定井下水位上升量和含水層給水度，就可以利用該方法計(jì)算出補(bǔ)給量。Healy(2007)[8]提出井下水均衡的本質(zhì)是質(zhì)量守恒原理，進(jìn)行均衡情況分析，不受時(shí)空尺度的限制,對(duì)流域的水均衡提出了水量平衡方程式。

1 研究區(qū)概況

研究區(qū)屬于暖溫帶半濕潤(rùn)季風(fēng)性大陸氣候，多年平均降水量500.4mm，多年平均徑流深17.9mm多年平均水面蒸發(fā)量1230.9mm。研究區(qū)水文分類(lèi)主要分淺層和深層井水，淺層水為孔隙潛水。深層水以內(nèi)丘至邢臺(tái)弧形斷裂為界，以西為巖溶水，以東為孔隙水。淺層含水層與深層含水層之間由10m～40m左右的粘性土隔離，基本沒(méi)有水力聯(lián)系。本次進(jìn)行淺層水水均衡分析。

研究區(qū)地下水主要接受大氣降水、灌溉回歸和地下水側(cè)向徑流補(bǔ)給，然后由西向東徑流，在東部邊界排出區(qū)外。淺層水文地質(zhì)圖見(jiàn)圖1。

2 水均衡分析

本次研究礦區(qū)水文均衡區(qū)為邢臺(tái)市開(kāi)發(fā)區(qū)行政界內(nèi)所包括的范圍，面積較大。

2.1 均衡方程

根據(jù)均衡區(qū)各地層水文補(bǔ)給和排泄條件，地下水補(bǔ)給包括：大氣降水入滲補(bǔ)給、井灌回歸補(bǔ)給量、側(cè)向徑流補(bǔ)給。井下水排泄包括：人工開(kāi)采（包括礦山開(kāi)采及礦區(qū)周邊生活用水開(kāi)采）、側(cè)向徑流排泄。

根據(jù)研究區(qū)均衡項(xiàng)，建立均衡方程式如下：

式中：Q降滲—大氣降水入滲補(bǔ)給量；Q灌—井灌回歸補(bǔ)給量；Q側(cè)補(bǔ)—井下水側(cè)向徑流補(bǔ)給量；Q開(kāi)采—孔隙水開(kāi)采量；Q側(cè)排—井下水側(cè)向徑流排泄量；Q儲(chǔ)變—井下水儲(chǔ)變量。

2.2 均衡項(xiàng)計(jì)算

根據(jù)《河北省邢臺(tái)市區(qū)水資源中長(zhǎng)期供求計(jì)劃報(bào)告》，確定均衡計(jì)算參數(shù)，降水入滲系數(shù)粘性土α=0.238、砂類(lèi)土α=0.446；潛水給水度μ=0.06；灌溉回歸系數(shù)β=0.12；滲透系數(shù)K=55m/d。

2.2.1 井下水補(bǔ)給量的計(jì)算

（1）降水入滲補(bǔ)給量

計(jì)算公式：

式中：Q降-降水入滲補(bǔ)給量（104m3/a）；P-年降水量(mm/a)；α-降水入滲系數(shù)；F-計(jì)算區(qū)面積（km2），粘性土面積37km2、砂類(lèi)土面積21km2。

（2）井灌回歸量

圖1 淺層孔隙水文地質(zhì)圖

計(jì)算公式：

式中：Q井-井灌回歸量（104m3/a）；Q灌開(kāi)-井灌開(kāi)采量（104m3/a）；β-井灌回歸系數(shù)，取0.12。

（3）側(cè)向流入量

計(jì)算公式：

式中：Q側(cè)入-側(cè)向流入量（104m3/a）；K-計(jì)算斷面含水層滲透系數(shù)（m/d）；M-計(jì)算斷面含水層厚度（m）；B-計(jì)算斷面寬度（m）；I-垂直計(jì)算斷面的水力坡度。

根據(jù)研究區(qū)水文地質(zhì)條件，水力坡度取0.0033；計(jì)算斷面按實(shí)際寬度取7.0km；補(bǔ)給斷面的含水層平均厚度取19m。

經(jīng)計(jì)算斷面?zhèn)认驈搅餮a(bǔ)給量為881.1×104m3/a。

綜上，均衡區(qū)內(nèi)淺層井下水年均總補(bǔ)給量為1892.0×104m3。

2.2.2 井下水量計(jì)算

（1）人工開(kāi)采量

地下水開(kāi)采包括：工業(yè)開(kāi)采和生活用水開(kāi)采。工業(yè)開(kāi)采量和生活用水開(kāi)采量依據(jù)實(shí)測(cè)統(tǒng)計(jì)值，農(nóng)業(yè)開(kāi)采量依據(jù)定額法計(jì)算值。該區(qū)人工開(kāi)采地下水總量為164.7×104m3/a。

（2）井下水側(cè)向徑流排泄量

依據(jù)計(jì)算公示（4），均衡區(qū)東部邊界為排泄斷面，其水力坡度取平均值0.0017，斷面寬度6.0km，計(jì)算斷面徑流排泄量為389.1×104m3/a。

綜上，均衡區(qū)內(nèi)井下水年均總排泄量為193.8×104m3。

2.2.3 儲(chǔ)變量的計(jì)算

均衡區(qū)內(nèi)孔隙水儲(chǔ)變量計(jì)算公式為：

式中：Q儲(chǔ)變-地下水儲(chǔ)存量的變化量（104m3）；μ-水位變動(dòng)帶的給水度；F-均衡區(qū)面積（km2）；ΔH—均衡時(shí)間內(nèi)年均水位變差（m）。

根據(jù)研究區(qū)水文地質(zhì)條件，潛水層給水度0.06，均衡區(qū)面積58km2，均衡期內(nèi)研究區(qū)地下水水位年平均下降0.19m。

均衡期內(nèi)井下水儲(chǔ)存量的變化量年均減少66.1×104m3。

3 結(jié)論

在分析研究礦區(qū)水文地質(zhì)條件的基礎(chǔ)上，總結(jié)了淺礦層水文補(bǔ)給和排泄項(xiàng)。計(jì)算得出均衡區(qū)內(nèi)靜下水年均總補(bǔ)給量為182.0×104m3，年均總排泄量為193.8×104m3，排泄量大于補(bǔ)給量，為負(fù)均衡，年均均衡差為-50.8×104m3。均衡期內(nèi)地下水年均儲(chǔ)存量變化量減少11.1×104m3，與均衡差比較接近。