南水北調(diào)東線調(diào)水仿真模型研究

南水北調(diào)東線調(diào)水仿真模型研究

1北銀東線工程的主要特點1.1洪水河流水系銀川項目是中國銀川項目的重要組成部分。東線工程主要供水區(qū)域為黃淮海平原和山東半島,主要目標(biāo)為解決天津、河北、山東半島和沿線城市生活和工業(yè)用水,并適當(dāng)兼顧環(huán)境、農(nóng)業(yè)和其它用水。該調(diào)水工程自長江干流引水,利用京杭大運河以及與其平行的河道輸水。黃河以南地勢北高南低,設(shè)有13級泵站,抽水總揚程65m,連通洪澤湖、駱馬湖、下級湖、上級湖、東平湖等湖泊作為調(diào)節(jié)水庫,在山東省位山附近黃河河底建穿黃隧洞,調(diào)水過黃河。黃河以北地勢南高北低,可順大運河或新擴建河道自流到天津,利用千頃洼、大屯、大浪淀、北大港4座平原水庫進(jìn)行水量調(diào)蓄。從長江到天津北大港輸水主干線長1156km,其中黃河以南646km,穿黃17km,黃河以北493km;分干線總長795km,其中黃河以南629km,黃河以北166km。此外,為解決山東半島嚴(yán)重缺水的狀況,需從東平湖開辟膠東輸水分干線,向濟南、青島、煙臺、威海等城市供水,輸水線路總長701km。1.2兩側(cè)地區(qū)特別大清河地區(qū)南水北調(diào)東線工程供水范圍為:(1)江蘇省里下河地區(qū)以外的蘇北地區(qū)和里運河?xùn)|西兩側(cè)地區(qū),安徽省蚌埠市、淮北市以東沿淮、沿新汴河、沿高郵湖地區(qū)和山東省南四湖、東平湖地區(qū);(2)山東半島;(3)黃河以北大清河淀東清南天津平原,河北省黑龍港東部平原,運東平原和山東省徒駭馬頰河平原。南水北調(diào)東線工程調(diào)水規(guī)模見表1。2水生百合工程的模擬模型結(jié)構(gòu)2.1用仿真技術(shù)模擬水系統(tǒng)建設(shè)南水北調(diào)東線工程是一個復(fù)雜的,具有多水源、多供水對象、多湖泊的水資源系統(tǒng),其特點主要表現(xiàn)在:(1)工程供水范圍較大;(2)工程供水對象多;(3)現(xiàn)有工程與規(guī)劃工程結(jié)合;(4)湖泊供水、抽水及棄水相結(jié)合;(5)多種水源聯(lián)合調(diào)度。這些特點的存在,使得模型建立的難度大大增加,對于這樣一個系統(tǒng),采用一般的數(shù)學(xué)優(yōu)化技術(shù)建模,要做很大程度的簡化和數(shù)學(xué)處理,會造成許多重要信息的遺失和模型的失真,計算結(jié)果與實際情況會產(chǎn)生較大出入。仿真技術(shù)能直觀、形象地模擬水資源系統(tǒng)中來水、蓄水、供水、輸水等一系列運行過程以及水庫、湖泊和洼淀間的水力關(guān)系,具有很強的實用性。因此,本次研究應(yīng)用了仿真技術(shù)對南水北調(diào)東線工程的水量調(diào)度過程進(jìn)行模擬,在擬定的調(diào)水目標(biāo)下,根據(jù)水區(qū)水利工程特點,兼顧調(diào)水區(qū)和受水區(qū)的利益,確定合理的水量調(diào)度原則,建立南水北調(diào)東線工程仿真模型,通過仿真模型運行計算來模擬東線工程調(diào)水量在時間和空間的調(diào)度分配過程。2.2修正穿黃水量南水北調(diào)東線工程調(diào)水系統(tǒng)分為黃河以南和黃河以北兩個子系統(tǒng),以穿黃水量作為連接這兩個系統(tǒng)的紐帶。整個仿真模型的計算思路則是依據(jù)這個基本結(jié)構(gòu),逐時段進(jìn)行水量調(diào)配計算,先由黃河以南子系統(tǒng)開始,自南向北逐個湖泊進(jìn)行水量調(diào)節(jié)計算,并且在滿足系統(tǒng)各項要求的條件下,按照最大的抽水能力向北抽水,得到出東平湖的最大水量,扣去輸送山東半島的水量后,即得到穿黃水量。將穿黃水量輸入黃河以北子系統(tǒng),在黃河以北子系統(tǒng)中,聚合計算出在位山處的黃河以北總需水量,將其與穿黃水量進(jìn)行比較,對穿黃水量進(jìn)行修正,然后將修正后的穿黃水量作為黃河以南子系統(tǒng)的輸出,黃河以北子系統(tǒng)對修正后的穿黃水量進(jìn)行分區(qū)的水量配置。黃河以南子系統(tǒng)依據(jù)穿黃水量的修正量,自北向南進(jìn)一步修正各湖泊的抽水量、棄水量、湖泊蓄水量和抽江水量。由此,完成該時段的水量調(diào)配計算。如此循環(huán),可模擬出各年和整個水文系列的抽江、穿黃、供需水量狀況的東線工程調(diào)水過程的一系列過程和數(shù)據(jù)。2.3需水面積的確定南水北調(diào)東線工程系統(tǒng)復(fù)雜,是一個多水源、多用戶、多湖泊洼淀的水資源工程系統(tǒng)。因此,在建模中根據(jù)天然水系的分布和實際供水體系的設(shè)置,將這一水資源系統(tǒng)以黃河為界,分為黃河以南和黃河以北兩個子系統(tǒng),穿黃水量作為兩個子系統(tǒng)的控制變量。在黃河以南子系統(tǒng)中,洪澤湖、駱馬湖、下級湖、上級湖和東平湖5個湖泊作為東線工程的調(diào)蓄水庫,在系統(tǒng)概化時,以這5個湖泊為中心,突出湖泊的調(diào)節(jié)能力,反映各湖泊間的水量補償和調(diào)配作用,以及各個湖泊間存在的上下級間的水力關(guān)系,將各級抽水泵站群復(fù)合為湖泊的出、入湖站,同時,忽略輸水河道的調(diào)蓄作用,將其簡化為單一的輸水河道,將5個湖泊串聯(lián)在一起,形成東線工程黃河以南的水資源子系統(tǒng)。南水北調(diào)東線工程的黃河以北為受水地區(qū),根據(jù)該地區(qū)的供水特點,劃分為6個供水分區(qū):山東省的聊城區(qū)、德州區(qū)、河北省的衡水區(qū)、景縣區(qū)、滄州區(qū)和天津區(qū),這6個分區(qū)分別通過6個分水口與輸水干線相連接。在本區(qū)內(nèi)有4座規(guī)劃的調(diào)蓄工程:北大港水庫、大浪淀水庫、千頃洼水庫和大屯水庫。將每一個供水分區(qū)與其相應(yīng)的調(diào)蓄工程一起概化成一個計算單元,來進(jìn)行水量供需分析。本次模擬中考慮到要對山東半島的城市用水進(jìn)行補充,故建模時,將其概化為一個計算單元。南水北調(diào)東線工程系統(tǒng)的概化見圖1。2.4智慧湖泊群水資源配置南水北調(diào)東線工程仿真模型主要由黃河以南計算模塊、黃河以北計算模塊和計算成果瀏覽模塊等組成。各個模塊根據(jù)其功能的不同和計算要求不同,又分別由湖泊群水量計算模塊、水量供需分析模塊、穿黃水量修正模塊、黃河以北配水模塊、黃河以北水庫補水模塊、統(tǒng)計計算模塊組成,各模塊相互間的關(guān)系見南水北調(diào)東線工程調(diào)水仿真模型結(jié)構(gòu)圖2。黃河以南計算模塊包括了湖泊群水量計算模塊、水量供需分析模塊、穿黃水修正模塊和統(tǒng)計計算模塊。本模塊主要對黃河以南5個湖泊進(jìn)行水量平衡計算,分別對5個湖泊所在分區(qū)的工農(nóng)業(yè)用水進(jìn)行供需分析。對計算成果進(jìn)行統(tǒng)計計算,通過成果瀏覽模塊顯示計算結(jié)果。黃河以北計算模塊包括了水量分配子模塊、水庫補水子模塊和計算成果統(tǒng)計子模塊。本模塊主要是先聚合計算黃河以北的需水量,然后根據(jù)黃河以南模塊提供的穿黃水量,在黃河北水量分配子模塊中,計算本區(qū)可以使用的水量也就是實際穿黃水量,同時向黃河以南反饋對穿黃水量的修正量。并按照利益共享的原則,將實際穿黃水量分配到各個區(qū),進(jìn)行各區(qū)的水量供需平衡計算,再根據(jù)黃河以北的補庫要求,運用黃河北水庫補水子模塊,將穿黃水量扣除黃河北的需水量后多余的水,補充到黃河北的水庫中,進(jìn)行水庫補水調(diào)蓄作用的計算,同時也計算當(dāng)穿黃水量不能滿足黃河北需調(diào)水量時,由水庫的蓄水來補充不能滿足的部分。計算成果瀏覽模塊的功能是根據(jù)用戶的要求,從仿真模型計算成果庫中調(diào)出所需瀏覽的仿真模型計算結(jié)果,以報表的形式在計算機屏幕上顯示,并提供了打印機輸出的功能。2.5量調(diào)配原則東線工程仿真模型不僅要模擬東線工程各個湖泊間的水力聯(lián)系,并且要通過選擇合理的水量調(diào)配原則,來對地區(qū)之間、湖泊之間、用水部門之間的水量進(jìn)行協(xié)調(diào),以求達(dá)到水資源合理配置的目的。因此,東線工程仿真模型,要滿足東線工程調(diào)水系統(tǒng)中的各種約束條件,主要包括水量平衡約束、湖泊蓄水約束、河道輸水約束、抽水站規(guī)模約束和供用水約束等。(1)水質(zhì)vi費即湖泊時段末庫容等于時段初庫容加來水、抽進(jìn)水量減去湖泊的蒸發(fā)、滲漏損失、供水量、抽出水量和棄水量。其數(shù)學(xué)表達(dá)式為:式中:Vi,j—湖泊蓄水量;Qi,j—湖泊來水量;P1i,j、P2i,j—湖泊抽進(jìn)水量、抽出水量;Ei,j—湖泊蒸發(fā)損失;Ri,j—供水量;Y1i,j—上一湖泊棄水量;Yi,j—本湖泊棄水量。(2)湖泊的蓄水是有限的湖泊蓄水要低于等于蓄水上限,高于等于蓄水下限。(3)流量站和站群的限制抽水量要小于等于抽水泵站的抽水能力。(4)河流的流量受到限制通過河道輸水的流量要低于等于河道輸水能力。(5)供需水量限制供給生活和工業(yè)、農(nóng)業(yè)的水量要小于相應(yīng)的需水量。(6)不接受額外的限制上述約束中的各變量均為非負(fù)。2.6確保湖泊供水量(1)尊重已用水源現(xiàn)實,調(diào)水工程的興建以不影響原有灌區(qū)合理用水為原則。各地對當(dāng)?shù)厮蓛?yōu)先使用,多余水量由工程統(tǒng)一調(diào)度。(2)黃河以南地區(qū)城市生活、工業(yè)和航運用水,農(nóng)業(yè)規(guī)劃灌區(qū)用水,都要達(dá)到各自供水保證率要求。(3)湖泊調(diào)度圖包含兩條控制線:一是北調(diào)水量控制線(v2),其作用是控制本級湖泊向上級湖泊的北調(diào)抽出水量,以保證本級湖泊供水對象的用水要求;在一般情況下,當(dāng)?shù)陀诒闭{(diào)控制線時,停止向北調(diào)水,二是用水控制線(v1),由于湖泊供水對象具有不同的供水保證率要求,有必要設(shè)置用水控制線以保證其不同要求;由于工業(yè)生活用水的保證率要高于農(nóng)業(yè)灌溉保證率,因此,當(dāng)湖泊蓄水低于用水控制線(v1)時,則停止向農(nóng)業(yè)供水,以保證本片下一時段工業(yè)和生活用水。(4)為了保證黃河以北城市工業(yè)生活用水,設(shè)置最小向北調(diào)水量,在湖泊停止向北調(diào)水時,利用調(diào)水工程新增的裝機向北調(diào)水。最小向北抽水量小于等于東線各級新增抽水能力。(5)當(dāng)穿黃水量大于等于黃河以北的總需水量時,黃河以北各區(qū)的需水全都滿足,仍有多余水時,則需要修正穿黃水量。當(dāng)穿黃水量小于黃河以北總需水時,各區(qū)先用水庫的水供各區(qū)的需水,并重新計算黃河以北的總需調(diào)水后,再與穿黃水量比較,修正各區(qū)的供水量。(6)當(dāng)考慮用穿黃水量補水庫蓄水時,穿黃水首先要在滿足黃河以北各區(qū)的需水后,如仍有多余的水,才能進(jìn)行補庫。當(dāng)補庫后仍有多余水量時,則需要修正穿黃水量。3模擬模型的應(yīng)用在東部工程3.1工程供水區(qū)分區(qū)南水北調(diào)東線工程供水范圍包括黃河以南、山東半島和黃河以北3大片,涉及江蘇、安徽、山東、河北和天津等5省市。根據(jù)東線工程供水的實際情況和仿真模型計算的要求,進(jìn)一步將黃河以南地區(qū)分為5個供水分區(qū),黃河以北地區(qū)分為6個供水分區(qū),山東半島為1個分區(qū),在仿真模型中共有12個計算分區(qū)。分區(qū)見表2。3.2固定組合方案按照南水北調(diào)東線工程調(diào)水仿真模型的功能,計算方案的選擇有多種組合,但主要取決于基本數(shù)據(jù)的詳細(xì)程度,表3列出了所有可能的組合方案。根據(jù)目前所掌握的東線基本數(shù)據(jù)情況,在上述南水北調(diào)東線工程調(diào)水計算方案參數(shù)的組合中,按照各湖的規(guī)劃來水和抽水工程規(guī)模情況,考慮南水北調(diào)東線工程第一、二、三期不同的抽水規(guī)模,不同水平年工農(nóng)業(yè)的需水要求,以及模型計算采用的不同供水方式,選擇、組合成計算方案。3.3計算方案的參數(shù)描述3.3.1北魯西南、東北部南水北調(diào)東線第一期工程供水范圍包括了江蘇蘇北、安徽東北和山東的魯西南、膠東和魯北地區(qū)。南水北調(diào)東線第二期工程供水范圍在第一期的基礎(chǔ)上向北伸延,增加了河北和天津。南水北調(diào)東線工程第三期的供水范圍同第二期。各期抽水規(guī)模見表1。3.3.2限制了產(chǎn)水及需調(diào)水量2010水平年和2030水平年的城市、工農(nóng)業(yè)需水量是以淮委《南水北調(diào)東線工程規(guī)劃——水量調(diào)配》報告中預(yù)測的總量為控制,以及各分區(qū)工農(nóng)業(yè)布局和國民經(jīng)濟發(fā)展的情況分配到各分區(qū)。在東線第一、二期工程中,黃河以南地區(qū)采用的是2010水平年的城市、工農(nóng)業(yè)需水量,為225.88億m3,在東線第三期工程中采用的是2030水平年的城市、工農(nóng)業(yè)需水量,為249.2億m3;黃河以北地區(qū)需調(diào)水量在東線第一、二、三期工程中都不相同;山東半島需調(diào)水量在東線第一、二期工程中采用了2010水平年的需調(diào)水量,為8.76億m3,在第三期工程中采用了2030年的需調(diào)水量,為21.29億m3。3.3.3北方單段水質(zhì)調(diào)度系統(tǒng)在模型計算時有兩種供水方式可以選擇,第一種稱之為按需供水方式,是指模型在東線工程調(diào)水計算過程中,以滿足北方本時段的需調(diào)水量為目標(biāo),進(jìn)行水量調(diào)度。第二種供水方式稱之為補庫供水方式,是指模型在東線調(diào)水計算過程中,不僅考慮本時段黃河以北的需調(diào)水量,同時也考慮了北方水庫的蓄水狀況,當(dāng)黃河以南地區(qū)可以有多余水北調(diào),并且同時北方水庫有庫容可以蓄水時,則提前調(diào)水至水庫。3.3.4農(nóng)業(yè)用水控制線的設(shè)置在仿真模型的湖泊調(diào)節(jié)計算中,對黃河以南地區(qū)各湖泊的水量調(diào)度采用了兩條調(diào)度控制線:第一條是為了保證當(dāng)?shù)毓┧O(shè)置的控制湖泊向北抽水的北調(diào)控制線。計算時當(dāng)湖泊水位低于北調(diào)控制線時,為了保證當(dāng)?shù)氐挠盟?一般將停止向其上游湖泊抽水。第二條是為了滿足黃河以南各計算分區(qū)的城市、工農(nóng)業(yè)用水不同供水保證率的要求,而設(shè)置的農(nóng)業(yè)供水控制線。計算時當(dāng)湖泊的水位低于供水控制線時,本時段將停止向當(dāng)?shù)剞r(nóng)業(yè)供灌溉水,以保證當(dāng)?shù)叵乱粫r段城市生活、工業(yè)供水量。3.4資源供需平衡應(yīng)用南水北調(diào)東線工程仿真模型分別按照按需供水方式進(jìn)行調(diào)水計算和按照補庫供水方式進(jìn)行調(diào)水計算,根據(jù)各湖泊42a的長系列規(guī)劃來水進(jìn)行計算得到了多個方案的調(diào)水工程多年平均、不同頻率年抽江量、湖泊入湖和穿黃水量,沿線各地區(qū)水資源供需平衡計算結(jié)果。表4給出了多年平均抽江量、湖泊入湖和穿黃水量。從模型計算結(jié)果來看,東線第一期工程(方案1)總供水量為192.5億m3,黃河以南供水量為178.98億m3,其中城市生活、工業(yè)和航運供水量的供水保證率達(dá)到100%,農(nóng)業(yè)供水量供水保證率達(dá)到86.7%以上;黃河以北供水保證率達(dá)到96.8%以上。滿足了兼顧調(diào)水區(qū)和受水區(qū)的利益的要求,在對8個典型枯水年計算結(jié)果分析發(fā)現(xiàn),在枯水年汛前的4—6月份,東線工程調(diào)過黃河水量是無法滿足黃河以北需調(diào)水量。因此,在枯水年東線工程調(diào)水時要考慮這個因素。從模型計算結(jié)果分析來看