小南海泉域地下水資源管理模型

1、倪付燕1 河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇南京(210098摘要:本文選取北方巖溶地區(qū)的小南海泉域?yàn)檠芯繉?duì)象,將數(shù)值模型與線性規(guī)劃相結(jié)合,建立小南海泉域地下水資源管理模型,在滿足泉域工農(nóng)業(yè)生活用水需求、維持小南海泉流量3.5/s 的前提下,利用單純形法尋求最優(yōu)的地下水開采量。計(jì)算結(jié)果表明,泉域枯水年份地下水已經(jīng)超采,平水年份具有一定的開采潛力,豐水年份可以多開采以豐補(bǔ)歉。 3m 關(guān)鍵詞: 小南海泉域 數(shù)值模型 線性規(guī)劃 管理模型1 引言小南海泉位于河南省安陽市,是安陽市的主要供水水源之一,擔(dān)負(fù)著安陽市區(qū)近一半的供水任務(wù),是安鋼、電廠、化肥廠等企業(yè)以及萬金灌區(qū)的主要水源,同時(shí)

2、也是市區(qū)生活用水的重要后備水源,被稱為安陽的“生命泉”。小南海泉作為泉域地下水的主要排泄方式之一,流量大小與泉域地下水位息息相關(guān)。近年來,由于泉域內(nèi)打井、挖煤、開礦等現(xiàn)象嚴(yán)重,地下水位下降,泉流量逐年減小, 如何確保安陽市的“生命泉”長(zhǎng)盛不衰,并且最大限度地開采地下水滿足泉域內(nèi)的用水需求,是一個(gè)亟待解決的問題1。本文利用數(shù)值模型和線性規(guī)劃相耦合的管理模型建模方法,對(duì)小南海泉域的地下水資源進(jìn)行優(yōu)化管理。2 研究區(qū)概況小南海泉域包括林州市大部,安陽縣的西部以及鶴壁市西北部分地區(qū),面積為934.6。研究區(qū)屬于大陸性半干旱氣候區(qū),多年平均降水量為573.5,年內(nèi)降水集中在79月。研究區(qū)內(nèi)主要河流為洹河

3、,建有舉世聞名的“人工天河”紅旗渠。2km mm 研究區(qū)地勢(shì)總的趨勢(shì)是西高東低,地貌以中低山、丘陵為主。出露的底層有太古界、元古界和古生界寒武系、奧陶系、石炭系、二疊系及新生界第三系、第四系。區(qū)內(nèi)構(gòu)造形跡以斷裂為主,構(gòu)造不甚發(fā)育,除斷層附近外,地層產(chǎn)狀較平緩。泉域內(nèi)巖溶發(fā)育,地表巖溶不太發(fā)育,巖溶形態(tài)以溶溝、溶孔、溶洞為主,地下巖溶以溶蝕裂隙為主,大型溶洞不多,局部可形成強(qiáng)巖溶帶。根據(jù)地下水賦存介質(zhì)的性質(zhì)及其空隙特性,將小南海泉域地下水類型分為松散巖類孔隙含水巖組、碎屑巖類裂隙孔隙含水巖組、基巖裂隙含水巖組及碳酸鹽巖裂隙巖溶含水巖組四種類型(圖1。其中碳酸鹽巖裂隙巖溶含- 1 -本區(qū)巖溶水的補(bǔ)

4、給主要為大氣降水面狀滲入與河流線狀滲漏,其次為渠道滲漏及農(nóng)田灌溉回滲。泉域地下水徑流總體趨勢(shì)為由西向東,最終匯集于南海泉溢流帶。巖溶水主要以泉的形式排泄,其次是工農(nóng)業(yè)開采地下水和煤礦排水。 A 基巖裂隙水B 松散巖類孔隙水C 裸露型碳酸鹽巖裂隙巖溶水D 碎屑巖類裂隙巖溶孔隙水圖1 小南海泉域水文地質(zhì)圖3 區(qū)域數(shù)值模型及水文地質(zhì)參數(shù)計(jì)算泉域內(nèi)的含水層概化為具有同一水力聯(lián)系的非均質(zhì)各向同性的潛水含水層,不考慮地下水垂向運(yùn)動(dòng),地下水流為水平二維流。本研究以2001年12月繪制的地下水流場(chǎng)作為模型的初始流場(chǎng)。泉域內(nèi)邊界條件根據(jù)地質(zhì)條件概化為流量邊界或者隔水邊界。由于工農(nóng)業(yè)開采發(fā)生在整個(gè)泉域范圍內(nèi),工農(nóng)

5、業(yè)開采量概化為面狀單位開采率;根據(jù)泉域內(nèi)煤礦的分布,礦坑排水概化為單個(gè)的開采井。根據(jù)水文地質(zhì)概念模型,泉域內(nèi)的地下水運(yùn)動(dòng)規(guī)律描述如下:(12002(,0(,(,(,0(,(,0y x y t h h h S h b K h b K W x y t T t x x y y h x y t h x y x y t h h b K q x y t x y t Tn =+<=G =- 2 -x、K y分別為x、y方向滲透系數(shù)分量,S y為給水度,W為源匯項(xiàng),包括機(jī)井開采強(qiáng)度、降水補(bǔ)給強(qiáng)度以及河流滲漏補(bǔ)給強(qiáng)度,h為地下水位,q為第二類邊界2上的流量,h0為初始水位,b為隔水層頂板高程。本計(jì)算采用有

6、限差分方法求解方程(1, 利用商業(yè)軟件VisualMODFLOW進(jìn)行模型設(shè)計(jì)和計(jì)算。利用2002年和2003年實(shí)測(cè)水文、氣象、地下水水位和泉流量資料,率定研究區(qū)水文地質(zhì)參數(shù),擬合結(jié)果見圖2。利用2004年和2005年實(shí)測(cè)資料進(jìn)行模型驗(yàn)證,驗(yàn)證結(jié)果見圖3。率定的水文地質(zhì)參數(shù)結(jié)果見圖4-1、4-2、4-3。從圖中可以看出,巖溶區(qū)由于巖性區(qū)域性很強(qiáng),參數(shù)空間分布變化較大,滲透系數(shù)最小0.6m/d,接近阻水巖層,最大在強(qiáng)徑流帶,達(dá)600m/d。分區(qū)降水入滲補(bǔ)給系數(shù)加權(quán)值:松散層覆蓋區(qū)為0.22、裸露區(qū)為0.3552。 圖4-1 降水入滲補(bǔ)給系數(shù)分區(qū)圖圖4-2 滲透系數(shù)分區(qū)圖(單位:m/d- 3 - 圖

7、4-3給水度分區(qū)圖4管理模型4.1 數(shù)學(xué)模型為了進(jìn)行管理模型的計(jì)算,根據(jù)疊加原理,我們把模型(1分解為(2和(3兩個(gè)模型3-5。模型(2為沒有水井抽水,僅僅由于泉域內(nèi)的大氣降水補(bǔ)給、洹河滲漏等不可控因素及非零初始條件和邊界條件引起的地下水運(yùn)動(dòng)問題,其表達(dá)式為:2'00(,2(,0(,0(,(,0(,x y t x y h h h Sy h b K h b K W x y t T t x x y y x y t h x t t h h x y h b K q x y t n =+<=(2 t T 其中是沒有抽水的源匯項(xiàng),模型的解為h 'W H =。模型(3為泉域內(nèi)齊次的初始

8、條件和邊界條件下僅僅由于可控的人工抽水引起的地下水運(yùn)動(dòng)問題,表達(dá)式為:(3 - 4 -202(,(,0(,0(,0(,0(0x y t h h h Sy h b K h b K p x y t x y t T t x x y y h x t t x y t h x y h b K n =+<=(,p x y t h S=模型(1的解為兩個(gè)模型的解的代數(shù)和,h H S =+。4.2響應(yīng)函數(shù)的計(jì)算建立巖溶地下水資源管理模型,有多種方法,其中利用響應(yīng)矩陣法建立線性規(guī)劃管理模型是目前最常用的方法之一。本研究區(qū)含水層參數(shù)不隨時(shí)間變化,并且模型(3是一個(gè)線性時(shí)不變系統(tǒng),所以本區(qū)可以用線性規(guī)劃進(jìn)行地下

9、水資源的優(yōu)化管理。建立此類管理模型,首先要求取單位脈沖響應(yīng)函數(shù)。為了便于管理將泉域地下水開采量分為泉口附近的礦坑排水和泉域內(nèi)的工農(nóng)業(yè)開采兩部分,將管理時(shí)段確定為3個(gè)月,一年4個(gè)時(shí)段。利用模型(3求得單位脈沖響應(yīng)函數(shù)(,i j k ,表示第j 點(diǎn)以單位流量持續(xù)抽水,i 點(diǎn)處在第個(gè)時(shí)段末刻產(chǎn)生的累計(jì)降深與第-1個(gè)時(shí)段末刻產(chǎn)生的累計(jì)降深的差值k k 3。泉口附近的礦坑排水,以10000為單位流量在各個(gè)時(shí)段持續(xù)抽水,計(jì)算得到泉排泄口處單位脈沖響應(yīng)函數(shù)3/m d s (,1,s k (k=1、2、3、4的值。泉域內(nèi)的工農(nóng)業(yè)生活用水開采量,以10000為單位流量在各個(gè)時(shí)段持續(xù)抽水,計(jì)算得到泉排泄口處單位脈沖響應(yīng)函數(shù)3/m d s (,2,s k 的值。4.3管理模型用線性規(guī)劃進(jìn)行本區(qū)地下水資源的優(yōu)化管理。本區(qū)的管理目標(biāo)是在保持小南海泉流量3.5以上3/m s 2,且滿足泉域工農(nóng)業(yè)生活用水量的情況下,如何最大限度地開采地下水滿足日益增長(zhǎng)的用水需求。因此,管理目標(biāo)確定為地下水的總開采量最大,其表達(dá)式為: Max Z=(44211(,k j Q j k =其中為管理時(shí)段。k 根據(jù)年內(nèi)降雨量的大小,選取2002年枯水年、1998年平水年、2003年豐水年為典型年對(duì)泉域地下水資源進(jìn)行優(yōu)化管理。