張 弘

(遼寧省丹東水文局，遼寧 丹東118001)

1 水資源配置目的

海水入侵3個嚴(yán)重區(qū)域分別位于三大灌溉體系的海岸附近。本文配置水資源主要是針對海水入侵嚴(yán)重地區(qū)的用水量與時空分布，結(jié)合水動力學(xué)，根據(jù)線性規(guī)劃模型進行計算。

三大海水入侵嚴(yán)重區(qū)中從左至右，第1個區(qū)域為孤山鎮(zhèn)劉大房附近，該區(qū)域為鹽堿化沼澤地，農(nóng)業(yè)生產(chǎn)少，且地勢較低，與大洋河河口有直接的水力聯(lián)系。第2個區(qū)域是椅圈鎮(zhèn)于家村，為農(nóng)業(yè)區(qū)域。第3個區(qū)域是東港市新興區(qū)葦場和勝利村，一個在東港市區(qū)，為工業(yè)用水區(qū)，一個為農(nóng)業(yè)用水區(qū)［1］。對水資源進行配置以防治海水入侵，主要分為3個目的層，最低目的層:一般在枯水期，要保證海水入侵不加劇，水資源利用得到基本滿足;中間目的層:平水期，要保證農(nóng)業(yè)用水區(qū)、工業(yè)(城區(qū))用水區(qū)海水入侵在合理范圍，水資源利用得到滿足;高級目的層:一般在豐水期，在滿足用水的前提下，最大能力地回灌地下水，增加入海水量。同時，由于海水入侵的滯后性，保護效果也會有一定的滯后［2］。

2 水資源開發(fā)利用現(xiàn)狀

2.1 水資源開發(fā)利用率

研究區(qū)2012年水資源總量為12.34 億m3，總供水量7.59 億m3，耗水5.39 億m3，水資源開發(fā)利用率達到61.4%，消耗率達到71%，針對水資源總量，耗水比率為43.6%。另外研究區(qū)每年都從土門子水庫、鴨綠江提水利用，實際開發(fā)利用率要低于此值。

2.2 用水結(jié)構(gòu)分析

2012年用水結(jié)構(gòu)主要為農(nóng)業(yè)用水、工業(yè)用水和生活用水，其中生活用水有城市生活用水和農(nóng)村生活用水。農(nóng)業(yè)灌溉用水6.59 億m3，占總用水量87%，工業(yè)用水0.64 億m3，占總用水量8.4%，生活用水0.35 億m3，占總用水量4.6%。詳見表1。

表1 2012年用水結(jié)構(gòu)表

3 水資源開發(fā)利用預(yù)測

3.1 中長期水資源需水預(yù)測

3.1.1 預(yù)測方法

預(yù)測城鎮(zhèn)生活需水、農(nóng)村人畜需水、農(nóng)業(yè)灌溉需水采用定額法。即通過分析確定不同水平年的各種需水定額，乘以相應(yīng)的發(fā)展指標(biāo)，即得到不同水平年的需水量。

計算公式為:

式中:W 是需水量;Q 是用水定額;C 發(fā)展指標(biāo)，i 是不同的需水類別(城鎮(zhèn)生活、農(nóng)村人畜、農(nóng)業(yè)灌溉)。

預(yù)測工業(yè)各部門的需水采用萬元產(chǎn)值耗水量的方法，通過對萬元產(chǎn)值耗水量分析，確定合理的重復(fù)利用率，然后利用萬元產(chǎn)值用水定額推算不同水平年的工業(yè)需水量，公式如下:

式中:q 為萬元產(chǎn)值用水定額;G 為萬元產(chǎn)值;k 為重復(fù)利用系數(shù)。

3.1.2 預(yù)測的依據(jù)

本次規(guī)劃需水預(yù)測主要依據(jù)是:東港市《國民經(jīng)濟和社會發(fā)展第十一五年計劃和2010年遠景目標(biāo)綱要》、市計生委制定的《人口規(guī)劃》。水資源的合理利用是關(guān)系到國民經(jīng)濟健康發(fā)展的主要因素之一，無論現(xiàn)在還是將來，本市各部門的各個行業(yè)都要本著計劃用水、節(jié)約用水的原則，千方百計降低用水定額，提高單位水資源的效益。因此在需水預(yù)測中要充分考慮到科技進步、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整、行業(yè)技術(shù)改造、先進的節(jié)水技術(shù)和節(jié)水措施推廣使用等方面節(jié)水因素［3－4］。

根據(jù)本市的實際情況，確定城鎮(zhèn)生活用水和農(nóng)村生活用水保證率95%，農(nóng)業(yè)灌溉用水保證率為75%，工業(yè)用水保證率為95%。

3.1.3 預(yù)測結(jié)果

需水總量預(yù)測(如表2，圖1、圖2)

圖1 2020年需水結(jié)構(gòu)圖

圖2 2030年需水結(jié)構(gòu)圖

表2 需水總量預(yù)測表 萬m3

從圖中看出，農(nóng)業(yè)依然是主要用水戶，分別高達87%和85%，生活用水量逐步增加，工業(yè)用水量有減少趨勢［5］。

3.2 中長期水資源供水結(jié)構(gòu)與供需平衡分析

3.2.1 供水結(jié)構(gòu)

根據(jù)東港市水資源規(guī)劃，地下水控制開采量，2020年和2030年保持4200 萬方水不變。根據(jù)保證率，地表水、地下水供水量如下表3，其中地表水所占比例高達93%。

表3 供水結(jié)構(gòu)表

3.2.2 供需平衡分析

從4 表中可以看出:2020年缺水9 819.3 萬m3，其中鴨干和鴨～大缺水12 366 萬m3，大洋河余水2 546 萬m3。

2030年缺水11 792.7 萬m3，其中鴨干和鴨～大缺水14 270 萬m3，大洋河余水2 477.3 萬m3。

表4 供需平衡分析表 萬m3

3.3 節(jié)水預(yù)測

3.3.1 農(nóng)業(yè)節(jié)水

改造現(xiàn)有水利工程，完善配套工程與調(diào)水工程，節(jié)水改造工程完成后，灌溉定額平均降低256 m3/hm2，灌區(qū)每年節(jié)省水量1.43 億m3。

3.3.2 工業(yè)節(jié)水

工業(yè)節(jié)水預(yù)測根據(jù)定額計算，2020年定額比2010年 少10 方/萬元，392 億 元 產(chǎn) 值 節(jié) 約1 960 萬m3水。

4 水資源配置方案

基于保證海水入侵程度最小思想的水資源優(yōu)化配置模型，以入海水量最大化為目標(biāo)函數(shù)，以生存條件約束和承載能力約束組成約束空間，構(gòu)造出水資源配置方案生成模型［6］。

模型數(shù)學(xué)表述如下:

4.1 目的函數(shù)

入海水資源量最大化，即地下水入海量和河口水入海量最大:

式中:W 為入海水量;Q地為地下水初水量;W地引為地下水引水量;W農(nóng)用為農(nóng)業(yè)用水量;Q河為河流初水量;W河引為河流水引水量;W工用為工業(yè)用水量;W生用為生活用水量，1.2 是地下水防治海水入侵效應(yīng)系數(shù)，0.5、0.1、0.6、0.75 為用水回歸系數(shù)。

4.2 約束條件

4.2.1 生存條件約束

生存條件約束主要有4 條:

1)Q河－ W河引+0.1W農(nóng)用+ 0.6W工用+ 0.6W生用≥Min，W生態(tài)為保證生態(tài)需水量要求。

2)Q庫≥W庫引為保證水庫可引水量。

3)Q河≥W河引為保證河水可引水量。

4)Q地≥W地引為保證地下水可引水量。4.2.2 承載能力約束

W河引+W地引+W庫引≥W農(nóng)用+W工用+W生用為保證引水量大于用水量。

(Q河－W河引)×(C標(biāo)－C本底)≥W污為保證水質(zhì)安全(納污能力)。

5 結(jié) 論

結(jié)合2020年、2030年預(yù)測的供需分析，預(yù)測7.3 億m3和7.5 億m3的需水量和6.3 億的供水量，并計劃對灌區(qū)續(xù)建配套與節(jié)水改造，新建鴨綠江引水工程，大洋河新建攔河閘及改建沿河提水站工程等，可 供 水 量 達 到7.54 億m3，2010年 余 水0.24 億m3，2030年余水0.04 億m3，水資源總量基本滿足要求［7－8］。

根據(jù)地下水的1.2 作用系數(shù)，一般只關(guān)注河水約束條件，鴨綠江的生態(tài)需水和納污能力滿足約束條件，故在充分引用水庫水后不足部分引用河水。

本模型以防治海水入侵為唯一目的，在輸水線路的經(jīng)濟效益沒有考慮，又由于本地區(qū)水資源豐富，一般以需定供，三水配置約束條件一般處于閑置狀態(tài)。

［1］王濤，周旭東，李晶. 大連地區(qū)海水入侵成因分析及防治對策研究［J］. 東北水利水電，2008(10).

［2］陳鴻漢，王新民，張永祥，任仲宇. 濰河下游地區(qū)海咸水入侵動態(tài)三維數(shù)值模擬分析［J］. 地學(xué)前緣，2000年S2期.

［3］豐華麗，王超，李劍超. 干旱區(qū)流域生態(tài)需水量估算原則分析［J］. 環(huán)境科學(xué)與技術(shù)，2002(01).

［4］高學(xué)平，楊建華，涂向陽，張明星. 帷幕防治海水入侵的數(shù)值模擬研究［J］. 海洋環(huán)境科學(xué)，2006，25(04):55－58.

［5］劉青勇，耿樹德，董廣清，等. 王河流域防治海水入侵的攔蓄源及調(diào)度技術(shù)研究［J］. 水利水電技術(shù)，1997，28(08):18 －21.

［6］王國利，周惠成. 大連市地下水資源優(yōu)化管理模型［J］.大連理工大學(xué)學(xué)報，2001，41(01):112 －115.

［7］于百榮，馮厚群. 龍口市北馬河流域緩解海水入侵試驗研究［J］. 水利水電技術(shù)，1996(09):5 －10.

［8］張展羽，高玉芳，李龍昌，徐征和. 沿海缺水罐區(qū)水資源優(yōu)化調(diào)配藕合模型［J］. 水力學(xué)報，2006，37(10):1246－1252.