林甸縣地下水可開采資源量評價

2010-06-26 13:38:32靳和平王晨光

黑龍江水利科技 2010年2期

吳丹,靳和平,王晨光

(林甸縣水務局,黑龍江林甸 162300)

1 水量均衡評價

1.1 均衡法原理

水量均衡法是地下水資源數量評價的基本方法。它的原理基礎是質量守恒定律,實質上是在某一時間段內確定某一地下水系統的補給量和排泄量與地下水儲存量變化間的關系。根據地下水資源可分為補給資源量和儲存資源量兩大類的原則,采用水量均衡法可計算出不同條件下地下水補給資源量,同時對地下水儲存資源的可恢復性和調節功能進行分析,就可以確定出地下水的可開采量,從而達到地下水資源評價的目的。對于一個地下水含水層系統來說,任意一個時間段的補給量與排泄量之差,必然等于該時間段內地下水儲存量的變化。即:

式中:Q1為均衡時段內地下水的天然補給量;Q0為均衡時段內地下水的天然排泄量;u為對潛水為給水度,對承壓水為貯水系數;A為均衡區含水層的分布面積;⊿h為均衡時段內地下水位的變化量。

1.2 均衡方程式

在天然條件下,地下水的補給量和排泄量基本上是均衡的。在開采條件下,由于人為地增加了地下水的排泄量,這必然改變地下水的天然排泄條件,會形成新的均衡方程式:

式中:⊿Q1為開采條件下地下水補給量的增加量;⊿ Q0為開采條件下地下水排泄量的減少量;Q人工為地下水的人工開采量。

由于地下水多年天然補給量和天然排泄量基本相等(Q1≈Q0),因此,開采條件下水均衡式可簡化為:

可見,地下水可開采量由 3個部分組成:開采條件下的補給量增量、天然排泄量的減少量及地下水儲存量的減少量。

2 地下水可開采量計算

2.1 潛水補給量

潛水總補給量為垂向補給量,包括大氣降水入滲補給量、地表水體滲漏補給量、農田灌溉回歸量。

潛水補給量計算式為:

式中:Q降為大氣降水入滲補給量;Q灌為田間灌溉入滲補給量;Q渠為渠系滲漏補給量;Q河為河道滲漏補給量;Q湖為湖泡水滲補給量。

2.2 潛水總排泄量

潛水總排泄量包括蒸發排泄量、人工開采量、向承壓含水層垂向越流排量、側向排泄量。式中:Q潛蒸為潛水蒸發量;Q側排為河流排泄量;Q越為潛水區通過弱透水層向承壓水區的排泄量;Q開采為人工開采潛水量。

2.3 承壓水補給量

承壓水總補給量主要包括越流補給量和側向補給量。承壓水補給量計算式為:

式中:Q側排為承壓含水層側向徑流補給量;Q越為越流補給量。承壓水總排泄量主要為人工開采量和側向排泄量。

式中:Q側排為側向排泄量;Q人工為人工開采消耗量。

3 林甸縣潛水補、排總量評價

3.1 潛水補給總量

3.1.1 大氣降水入滲補給量

式中:Q降為年大氣降水滲入補給量,108m3/a;P為年降水量,m m;a為年降水入滲補給系數;F為計算面積,k m2。

3.1.2 田間灌溉入滲補給量

式中:B渠灌為田間灌溉入滲系數或井灌回歸系數;Q1為渠灌凈水量。

截止 2009年,全縣實有耕地面積12.47萬h m2,其中水田地0.51萬 h m2,旱田11.96萬 h m2。

根據《黑龍江省行業用水定額》中水稻田淺型節水灌溉定額表(P=75%),林甸屬低平原、中等滲漏區,凈灌溉定額為7 695～8 700 m3/h m2,渠系利用系數取 0.6,則毛灌溉定額為12 825～14 490 m3/h m2,取灌溉定額為13 500m3/h m2。由于“五小”工程的灌區水田灌溉用水是短距離取地表水或地下水,所以在計算灌溉入滲補給量時,不計這部分灌溉水量。

3.1.3 渠系滲漏補給量

式中:m為渠系滲漏補給系數 m=r(1—n);r為修正系數,取值范圍一般在 0.3-0.8,取 0.4;n為渠系有效利用系數,取0.60;Q渠首引為渠首年引水量,其中渠系滲漏補給系數為0.16。

3.1.4 湖泡水入滲補給量

采用積水模型法建立入滲補給模型,計算公式為:

式中:K s為湖泡底泥導水率,取 2.1×10-4m/d;H為湖泡平均水深,取0.87m;Z s為潛水平均埋深,取4.60 m;F為湖泡水底面積,k m2;t為計算時段,取 365d。

林甸縣湖泡水體總面積為0.48萬h m2,入滲補給量為2 428萬 m3/a。

3.2 潛水排泄總量評價

3.2.1 潛水蒸發量

潛水蒸發量計算評價采用如下公式:

式中:Q為潛水蒸發量,億m3;E0為潛水接近地面時的蒸發強度,m m/a,一般用 E601型蒸發器的水面蒸發量近似代替,E0可根據蒸發器觀測的水面蒸發量E20估算而得,即E0=a.E20,a為經驗系數,這里取 0.53,即 E0=E20×0.53;⊿為潛水水位平均埋深,m;⊿0為潛水蒸發極限埋深,m;n為無綱量指數,與土壤質有關,一般取值范圍 1≤n≤3;F為計算面積,k m2。

根據本縣實際情況,利用地下水動態觀測資料,按潛水蒸發經驗公式分析計算,具體方法與動態觀測資料分析計算給水度相同。