馬素英 孫梅英 白振江 張爽

摘要：針對南水北調(diào)中線受水區(qū)外調(diào)水及當?shù)氐乇硭畞硭^程的不均勻性及來水頻率的不確定性、供用水系統(tǒng)的復(fù)雜性（多水源、多工程、多用戶）等特點，以河北省南水北調(diào)受水區(qū)為研究對象，以多目標群決策理論和大系統(tǒng)分解協(xié)調(diào)理論為指導(dǎo)，以優(yōu)化與模擬為主要技術(shù)，構(gòu)建了基于MIKE BASIN的1+7受水區(qū)水資源優(yōu)化調(diào)配模型，制定了基于不同年份、不同保證率、不同節(jié)水水平、不同用水戶需求條件的受水區(qū)水資源優(yōu)化調(diào)配方案。在2020年供受水區(qū)來水頻率組合（引江水P=50%、當?shù)氐乇硭甈=50%）條件下，對基本節(jié)水及高效節(jié)水的優(yōu)化調(diào)配方案優(yōu)化效果進行了分析，以期在充分利用水資源、滿足需水要求的情況下，尋求多水源與多用戶之間最優(yōu)的配水方案，實現(xiàn)區(qū)域綜合效益最大化，促進區(qū)域經(jīng)濟可持續(xù)發(fā)展。

關(guān)鍵詞：南水北調(diào)受水區(qū)；水資源優(yōu)化調(diào)配；多水源；多用戶；多目標

中圖分類號：TV213 文獻標志碼：A 文章編號：

16721683（2018）05006611

Study on the unified allocation of water resources in the water receiving area of the Middle Route of the SouthtoNorth Water Diversion Project in Hebei

MA Suying1，SUN Meiying1，BAI Zhenjiang2，ZHANG Shuang1

（

1.Hebei Institute of Water Resources，Shijiazhuang 050051，China；2.Hebei Branch of Construction and Administration Bureau of SouthtoNorth Water Diversion Middle Route Project，Shijiazhuang 050035，China）

Abstract：

Targeting the inhomogeneity of inflow process and uncertainty of inflow frequency of transferred water and local surface water in the water receiving area of the Middle Route of the SouthtoNorth Water Transfer Project as well as the complexity of the water supply system （multisource，multiproject，multiuser），we selected a water receiving area of the SouthtoNorth Water Transfer Project in Hebei province as the research object.With the guidance of multitarget group decision theory and large system decomposition coordination theory，we constructed a 1+7 optimal allocation model of water resources based on MIKE BASIN，and formulated an optimal allocation scheme for water resources in the water receiving area based on different years，different guarantee rates，different water saving levels，and different water users′ demand conditions.Under the condition of inflow frequency combination in 2020 （transferred water P=50% and local surface water P=50%），we analyzed the effects of the optimal allocation schemes for basic water saving and efficient water saving，with a view to finding the optimal water distribution scheme between multiple water sources and multiple users under the premise of utilizing water resources fully and satisfying the water demand.Exploring the optimal water distribution scheme between multiple water sources and multiple users will maximize the comprehensive benefits of the region and promote regional sustainable development.

Key words：

SouthtoNorth Water Transfer Area；optimal allocation of water resources；multiple water sources；multiple users；multiple objectives

南水北調(diào)工程是迄今為止世界上最大的水利工程，是我國水資源優(yōu)化配置的重大戰(zhàn)略性基礎(chǔ)設(shè)施[1]，其目的是用長江的水資源供給京、津、冀、豫及沿線城市工業(yè)及城鎮(zhèn)生活用水， 使其將本地水還給農(nóng)業(yè)，同時有效改善當?shù)厣鷳B(tài)環(huán)境。2014年12月12日中線一期工程干線全線通水，分配給河北304億m3（多年平均）。受配套工程建設(shè)進度等因素影響，截止2018年5月1日河北省累計利用引江水188億m3，一定程度上緩解了河北省南水北調(diào)中線受水區(qū)水資源短缺狀況。到2020年，引江水利用量將達到304億m3。不過，供水區(qū)不同年份來水過程在年內(nèi)、年際分配上存在不均勻性，而受水區(qū)當?shù)氐乇硭煌攴菀泊嬖诓痪鶆蛐裕盟畱粜杷侩S著經(jīng)濟的發(fā)展和生活水平的提高也存在不確定性，同時受水區(qū)水源系統(tǒng)十分復(fù)雜，既有引江水、引黃水、當?shù)氐乇硭?、地下水及非常?guī)水等多水源系統(tǒng)，也有引江供水工程、當?shù)氐乇硭┧こ?、地下水供水工程及污水處理回用工程等供水工程系統(tǒng)，還有生產(chǎn)、生活、環(huán)境等多用戶的復(fù)雜用水系統(tǒng)，由此受水區(qū)形成了一個多水源、多工程、多用戶的復(fù)雜水資源系統(tǒng)[2]。在此情況下，受水區(qū)必須進行水資源統(tǒng)一的優(yōu)化配置與合理調(diào)度，利用多種供水工程把多種水源調(diào)配給多端用水戶，保證這個復(fù)雜多端的水資源系統(tǒng)正常、有序、高效運行，才能實現(xiàn)“引江水高效利用、地表水調(diào)節(jié)利用、地下水控制利用”的高效用水目標。

目前我國針對南水北調(diào)受水區(qū)水資源調(diào)配的研究逐漸由 單水源單用戶調(diào)配向多水源、多用戶、多目標聯(lián)合優(yōu)化調(diào)配問題發(fā)展。閆寶偉等[3]運用Copula方法理論對水源區(qū)與各受水區(qū)的豐枯遭遇進行了研究，康玲等[4]運用貝葉斯網(wǎng)絡(luò)理論對南水北調(diào)中線工程水源區(qū)和4個受水區(qū)的降水豐枯遭遇進行了深入的分析研究，梁鐘元等[5]分析了氣候變化對南水北調(diào)中線受水區(qū)徑流量影響，上述研究為南水北調(diào)工程水資源調(diào)度研究奠定了基礎(chǔ)；廖衛(wèi)紅等[6]研制了南水北調(diào)中線受水區(qū)簡易地下水模型，程雙虎等[7]利用FEFLOW模型軟件建立了河北省南水北調(diào)受水區(qū)地下水模型，為科學、準確評價南水北調(diào)中線受水區(qū)地下水資源量奠定了基礎(chǔ)；劉昌明等[8]對南水北調(diào)中線主要城市節(jié)水潛力進行了分析，并提出了節(jié)水對策。朱永楠等[9]構(gòu)建了能夠全面反映區(qū)域用水結(jié)構(gòu)和全社會全過程用水效率的評價指標與評價標準，對南水北調(diào)受水區(qū)節(jié)水潛力進行了分析。上述研究為科學預(yù)測南水北調(diào)受水區(qū)各行業(yè)需水量奠定了基礎(chǔ)；耿六成等[10]開展了河北省南水北調(diào)受水區(qū)水資源合理配置研究，提出了科學的運行管理規(guī)則。吳澤宇等[11]建立了多水源大系統(tǒng)水資源配置模擬模型，開展了中線工程總干渠至各省市的分水量配置方法研究，實現(xiàn)了單一目標的多水源配置。張平等[12]建立了南水北調(diào)中線受水區(qū)的水資源優(yōu)化配置模型，王海云等[13]采用AHPLP法建立了河北省受水區(qū)水資源優(yōu)化配置模型，進行了省級層面行業(yè)間的多目標、多水源優(yōu)化配置。上述研究已充分考慮到水資源優(yōu)化配置的多優(yōu)化目標問題，同時也嘗試了多水源在多用水戶間的聯(lián)合調(diào)配，但研究層面均停留在?。ㄊ校┘?，無法結(jié)合省（市）級所轄的供水目標實際供水條件，做到真正意義上的水資源優(yōu)化配置。本文充分考慮了供水及節(jié)水的工程、政策、經(jīng)濟條件可行性，圍繞河北省南水北調(diào)受水區(qū)多水源聯(lián)合條件下各供水目標、各行業(yè)間的水資源統(tǒng)一調(diào)配問題進行了深入探索。

1.1 建模方法

本次研究采用優(yōu)化技術(shù)與模擬技術(shù)有機耦合的方法進行水資源優(yōu)化調(diào)配。選擇MIKE BASIN軟件模擬河北省南水北調(diào)中線受水區(qū)水資源系統(tǒng)。MIKE BASIN是由丹麥水利研究院（DHI）開發(fā)的[14]、完全與ArcGIS整合的一個用于綜合性水資源管理和規(guī)劃的決策支持工具，是一個可在流域（區(qū)域）尺度基礎(chǔ)上，解決水量的優(yōu)化配置、用水戶連接、水庫調(diào)度規(guī)則及水質(zhì)模擬等問題的綜合性水資源模擬軟件。數(shù)據(jù)信息和模擬結(jié)果通過MIKE BASIN在ArcGIS的平臺直觀展示，并可進行空間和時間尺度上的分析。

本研究引入多目標優(yōu)化算法，在MIKE BASIN水資源數(shù)學模型基礎(chǔ)上，通過模型調(diào)度/運行規(guī)則、內(nèi)置非線性優(yōu)化算法及基于軟件計算引擎進行的二次開發(fā)建立優(yōu)化模型，把水資源系統(tǒng)模擬規(guī)則概化為該算法能識別的決策變量，然后把該算法與模擬模型進行嵌套調(diào)用， 實現(xiàn)模擬與優(yōu)化的耦合。模擬與優(yōu)化技術(shù)耦合的方法詳見圖1。

1.2 建模思路

南水北調(diào)受水區(qū)水資源優(yōu)化調(diào)配系統(tǒng)是一個規(guī)模龐大、結(jié)構(gòu)復(fù)雜、功能綜合、影響因素眾多的大系統(tǒng)[15]，采用傳統(tǒng)單一整體模型求解往往難以克服系統(tǒng)復(fù)雜、維數(shù)（階數(shù)）眾多、求解困難等缺陷[16]。故采用分解協(xié)調(diào)等大系統(tǒng)優(yōu)化技術(shù)，建立河北省南

水北調(diào)中線受水區(qū)水資源優(yōu)化調(diào)配模型，將復(fù)雜的大[CM（22]系統(tǒng)分解為若干子系統(tǒng)，以解決優(yōu)化配置中更宏觀

的問題，使得建模及求解過程較容易。模型包含兩個級別，第一級（省級總模型），以市級計算單元為基礎(chǔ)，進行整個受水區(qū)的優(yōu)化調(diào)配；第二級（市級子模型），以縣級計算單元為基礎(chǔ)，進行各個子模型的優(yōu)化調(diào)配。模型運行時，在給定的供、需水條件下，先進行第二級（市級子模型）的優(yōu)化模擬，并將局部優(yōu)化結(jié)果g（X）反饋給第一級（省級總模型）?？偰Ｐ透鶕?jù)各子模型的反饋結(jié)果，計算總系統(tǒng)的優(yōu)化結(jié)果f（M）；同時將優(yōu)化結(jié)果在各子模型重新進行分配，并再次進行各子模型的優(yōu)化模擬，得到一組新的局部優(yōu)化結(jié)果和總系統(tǒng)優(yōu)化結(jié)果，如此反復(fù)進行，直到求得全局優(yōu)化結(jié)果。建模思路見圖2。

1.3 系統(tǒng)概化

根據(jù)建模思路，將受水區(qū)水資源系統(tǒng)設(shè)計為三層譜系結(jié)構(gòu)：第一層為河北省南水北調(diào)中線受水區(qū)整體協(xié)調(diào)層，第二層為市級子系統(tǒng)自控層（亦為市級計算單元），第三層為縣級基礎(chǔ)層（亦為縣級計算單元，即基本計算單元）。

受水區(qū)水資源系統(tǒng)構(gòu)成復(fù)雜，其組成要素主要包括點、線和面等信息。為將水資源系統(tǒng)中的各類元素有機聯(lián)系起來，實現(xiàn)對水資源系統(tǒng)網(wǎng)絡(luò)的描述，故在模型中將用水戶、水源、河流節(jié)點、分水口等概化為點，并將其分為兩類用水節(jié)點和供水節(jié)點；將河流、渠道等概化為線；將河網(wǎng)、受水區(qū)概化為面，即計算單元。經(jīng)概化后：系統(tǒng)節(jié)點704個，其中用水節(jié)點329個（生活用水節(jié)點110個，工業(yè)用水節(jié)點111個，環(huán)境用水節(jié)點108個）；供水節(jié)點375個（引江水供水節(jié)點112個，地表水供水節(jié)點38個，地下水供水節(jié)點112個，再生水供水節(jié)點112個，海水淡化水供水節(jié)點1個）。計算單元共119個，其中市級計算單元7個，縣級計算單元112個。詳見表1。

1.4 模型建立

1.4.1 目標函數(shù)

本次研究確定水資源優(yōu)化調(diào)配的總目標為，在滿足社會環(huán)境、生態(tài)環(huán)境和技術(shù)經(jīng)濟等方面限制條件下，追求經(jīng)濟效益、社會效益和環(huán)境效益的權(quán)衡最優(yōu)[17]。依據(jù)《河北省水資源綜合規(guī)劃》、《河北省實行最嚴格水資源管理制度實施方案》、《河北省南水北調(diào)中線配[HJ2.25mm]套工程規(guī)劃》、《河北省水中長期供求規(guī)劃》等成果中涉及到的優(yōu)水優(yōu)用原則、水資源優(yōu)化配置原則、南水北調(diào)中線工程建設(shè)目標、地下水保護政策等，在明確受水區(qū)供水水源及用水戶的供水、用水優(yōu)先級順序基礎(chǔ)上（詳見143運行規(guī)則），結(jié)合河北省實際工程條件及現(xiàn)狀用水習慣，最終確定了優(yōu)化目標的優(yōu)先級。模型目標函數(shù)設(shè)置詳見表2。

社會效益目標是一個不易度量的目標，僅以水資源對社會影響的角度考慮，可以認為缺水量大小或缺水程度直接影響到社會的發(fā)展和安定，是社會效益的一個側(cè)面反映，在需水一定的條件下，供水利用量最大即缺水量最小。河北省屬水資源嚴重短缺、地下水長期超采地區(qū)，生態(tài)環(huán)境需水長期被擠占，而南水北調(diào)工程的主要目的就是解決我國北方的上述問題。故本次研究確定社會效益目標為供水利用量最大及地下水開采量最小，其中供水利用量最大目標包括中水回用水利用量最大、海水淡化水水利用量最大、引江水利用量最大及當?shù)氐乇硭昧孔畲?。目標函?shù)如下：

式中：MaxF1，c-i 至MaxF4，c-i為求市級子模型中第1至4個目標函數(shù)最大值，其中F1至F4表示第1至4個目標函數(shù)；c-i 表示第i 個市級計算單元下屬縣級計算單元集合；YjRA，c-i為第i個市級子模型中第j個縣級計算單元的中水回用水利用量，YRA表示中水回用水利用量，j表示第j個縣級計算單元；MaxF1，m為求省級子模型中第1個目標函數(shù)最大值，m表示市級計算單元集合；MinF5，c-i為求市級子模型中第5個目標函數(shù)最小值；YDW、YDA、YSA、YGA分別為海水淡化水供水利用量、引江水供水利用量、當?shù)氐乇硭┧昧?、地下水供水利用量?/p>

經(jīng)濟效益目標用一定運行規(guī)則條件下供水帶來的直接凈效益最大表示，為最大區(qū)域經(jīng)濟效益與最低供水成本之差。

對于環(huán)境效益目標，考慮河北省實際，依據(jù)《河北省水中長期供求規(guī)劃》，全省近期“生態(tài)環(huán)境需水量”只能是通過“以供定需”方式基本達到“最小生態(tài)環(huán)境需水量”（詳見142約束條件），而對地下水的有效保護則通過取多年平均淺層地下水可開采量（在經(jīng)濟合理，技術(shù)可行且利用后不會造成地下水位持續(xù)下降、水質(zhì)惡化、海水入侵、地面沉降等環(huán)境地質(zhì)問題和不對生態(tài)環(huán)境造成不良影響的情況下，允許從地下含水層中取出的最大水量）作為地下水開采量上限值加以體現(xiàn)（詳見142約束條件）。故本次充分考慮優(yōu)水優(yōu)用配置原則，優(yōu)先配置處理達標的再生水作為生態(tài)環(huán)境用水水源，將廢污水回用量最大和環(huán)境缺水量最小作為環(huán)境目標。

廢污水回用量最大：

1.4.2 約束條件

為使水資源開發(fā)利用與生態(tài)環(huán)境保護相協(xié)調(diào)，使水資源利用促進經(jīng)濟社會發(fā)展，謀求水資源凈效益最大[18]，在進行優(yōu)化模擬計算過程中要考慮多方面的約束條件，其中主要包括行政約束、市場約束、生存需求約束、工程約束等。具體表達式詳見表3。

行政約束主要包括引江水分水量約束、供水目標供水保證程度約束。因中線工程分配給河北省的分水指標是一定的，故引江水利用量必須在此分水指標范圍內(nèi)，而對于受水區(qū)的112個供水目標，雖然亦有分水指標控制，但是考慮到河北省水資源的統(tǒng)一優(yōu)化管理調(diào)度，允許調(diào)配給受水區(qū)內(nèi)任意供水目標。至于供水目標保證程度約束，則根據(jù)用水戶（供水目標）的性質(zhì)不同來確定：對于生活用水戶來說供水保證程度需達到95%以上；工業(yè)用水戶供水保證程度需達到90%以上，最低需水量取需水量的90%；對于一些特殊地區(qū)（特區(qū)、工業(yè)園區(qū)）如雄安新區(qū)等，則依據(jù)正式出臺的相關(guān)政策文件、規(guī)劃及方案等，設(shè)置其不同供水目標的需水保證程度約束條件。

市場約束可解釋為由于水價的變化引起的用水戶需水量變化[19]，而形成的用水戶用水量約束。對于生活用水戶來說，其提出的需水量為其在現(xiàn)狀水價條件下需水量，當水價發(fā)生變化時，用水戶對水的需求量則隨之發(fā)生變化，變?yōu)橐?guī)劃水價條件下的需水量，而經(jīng)優(yōu)化調(diào)配后的生活用水戶利用水量需小于此值。同樣，對于工業(yè)用水戶來說，經(jīng)優(yōu)化調(diào)配后的工業(yè)用水戶利用水量需小于規(guī)劃水價條件下的需水量。

生存條件約束是從滿足人類基本需求角度，要求水資 源開發(fā)利用首先必須要解決人類生存問題，發(fā)展帶來的社會凈效益必須滿足人類最低生活標準的需要[20]，主要包括生活基本需水約束、社會凈效益約束、生態(tài)補水約束等。生活基本需水約束按兩種情景考慮：一是各用水戶根據(jù)自身需水要求上報的需水量，即基本節(jié)水條件下（方案一）需水量；二是本次研究參照《河北省實行最嚴格水資源管理制度實施方案》中“三條紅線”用水指標，確定的各用水戶需水量，即高效節(jié)水條件下（方案二）需水量。社會凈效益約束必須滿足人類最低生活標準的需要[20]，本次以河北省最低工資標準代替，同時采用萬元工業(yè)增加值取水量折算為工業(yè)最小需水量。另外考慮到電力行業(yè)是關(guān)系國計民生和社會發(fā)展的重要基礎(chǔ)產(chǎn)業(yè)，其用水需優(yōu)先滿足，故亦將充分考慮節(jié)能減排措施條件的電力需水量作為工業(yè)需水量的下限值。工業(yè)需水量上限值，同樣按兩種情景考慮，同生活基本需水約束，不再贅述。生態(tài)補水約束即環(huán)境用水量需滿足環(huán)境最低需水量，本次以現(xiàn)狀環(huán)境需水量代替。

工程約束主要包括水庫庫容約束、河/渠（管）道流量約束、地表水供水約束、地下水供水約束、中水回用水及海水淡化水約束，其中后四個約束包括供水量及供水能力兩方面約束。

1.4.3 運行規(guī)則

水資源配置模擬模型的核心是其在優(yōu)化模型基礎(chǔ)上確定的配置規(guī)則和運行規(guī)則，配置與運行規(guī)則直接影響模擬模型對實際情況的反映程度。在河北省南水北調(diào)受水區(qū)水資源系統(tǒng)模擬運算中，各供水節(jié)點與用水節(jié)點均需要設(shè)置具體的運行規(guī)則及優(yōu)先順序，故本次運行規(guī)則分為兩類：供水節(jié)點（供水水源）運行規(guī)則和用水節(jié)點（用水戶）運行規(guī)則。

供水節(jié)點（供水水源）運行規(guī)則。參考《河北省水資源綜合規(guī)劃》、《河北省實行最嚴格水資源管理制度實施方案》及《河北省南水北調(diào)中線配套工程規(guī)劃》，河北省南水北調(diào)受水區(qū)各水源使用的優(yōu)先次序是：非常規(guī)水源（淡化海水、中水回用水）優(yōu)先配置和使用；外調(diào)水優(yōu)先于水庫供水首先直供，有余水時存蓄于地下水庫或地表水庫；當?shù)氐乇韽搅骼脙?yōu)先于水庫供水，水庫供水除結(jié)合汛前泄水供水外，主要在汛后配合外調(diào)水用于城市、工業(yè)等；淺層地下水配合外調(diào)水和地表水參與調(diào)蓄、利用。

用水節(jié)點（用水戶）運行規(guī)則。規(guī)則可分為兩種：按用水類別可將用水節(jié)點分為生活用水戶、工業(yè)用水戶、環(huán)境用水戶，其優(yōu)先級順序為生活、工業(yè)、環(huán)境；按用水區(qū)域分的用水節(jié)點（用水戶）運行規(guī)則，是在按用水類別分的用水節(jié)點（用水戶）運行規(guī)則的基礎(chǔ)上設(shè)定的主要用水計算單元間供水量調(diào)配過程中的運行規(guī)則。按用水區(qū)域可將用水節(jié)點分為地下水資源嚴重匱乏地區(qū)用水單元（缺水率>50%）、嚴重缺水用水單元（30<缺水率≤50%）、一般缺水用水單元（0<缺水率≤30%）。涉及到計算單元間供水量調(diào)配的水源主要是引江水、地表水，優(yōu)先接受引江水調(diào)配水量的用水單元，依次為：地下水資源嚴重匱乏地區(qū)用水單元、嚴重缺水用水單元、一般缺水用水單元；按用水節(jié)點重要程度可將用水節(jié)點分為一類目標（國家級特區(qū)、設(shè)區(qū)市、重點工業(yè)區(qū)）、二類目標（縣級市）、三類目標（縣城）。優(yōu)先接受引江水和地表水調(diào)配水量的用水單元依次為一類目標（設(shè)區(qū)市、重點工業(yè)區(qū)）、二類目標（縣級市）、三類目標（縣城）。

2 模型求解

模型求解時，首先在供水預(yù)測與需水預(yù)測的基礎(chǔ)上，充分考慮基于決策者、用水戶及水市場而制定的約束條件，確定水資源優(yōu)化調(diào)配模型的供水方案集合和需水方案集合，形成水資源供需方案組合。其次，根據(jù)既定的目標函數(shù)及運行規(guī)則，采用順序決策法對確定的若干級優(yōu)化目標進行求解，高優(yōu)先權(quán)的目標比低優(yōu)先權(quán)的目標優(yōu)先得到滿足，每一個最優(yōu)化的目標均遵從模型約束限制，以確保其他具有更高優(yōu)先級的目標不受到負面影響。同一優(yōu)先級的多個優(yōu)化目標，則采用遺傳算法進行求解。通過優(yōu)化和模擬模型的相互反饋和不斷調(diào)整，得出不同水資源供需方案組合下的推薦調(diào)配方案集合。最后，將推薦方案提交至意見征詢系統(tǒng)，根據(jù)決策者、用水戶及水市場對推薦方案的反應(yīng)，調(diào)整調(diào)配模型的約束條件、運行規(guī)則及目標函數(shù)，重新運行優(yōu)化調(diào)配模型系統(tǒng)，得出推薦方案。這一過程一直進行到意見征詢系統(tǒng)無信息反饋，決策者滿意為止。最終的推薦方案即決策方案?？梢哉f本模型的求解過程是一個“滾動修正，不斷尋優(yōu)”的過程。受水區(qū)水資源優(yōu)化調(diào)配模型求解思路詳見圖4。

3 結(jié)果分析

本優(yōu)化調(diào)配模型可根據(jù)設(shè)定的水平年、節(jié)水條件（基本節(jié)水方案、高效節(jié)水方案）、配置方式等要素，或任意用水戶、任意決策者、任意水源、任意供水工程等水資源影響因素的不同，生成無數(shù)組優(yōu)化調(diào)配方案，本次以2020年供受水區(qū)來水頻率組合50%+50%調(diào)配方案（供受水區(qū)來水頻率均為P=50%、不同節(jié)水水平條件下2020年受水區(qū)優(yōu)化調(diào)配方案，下同）為例，對計算結(jié)果進行分析。

3.1 2020年供受水區(qū)來水頻率組合50%+50%基本節(jié)水方案調(diào)配結(jié)果

引江水利用情況分析。由圖5和表4可知，在2020年供受水區(qū)來水頻率均為P=50%、城市各業(yè)用水均為基本節(jié)水條件下，經(jīng)過優(yōu)化調(diào)配，受水區(qū)利用引江水量總計298億m3，利用率100%。其中邯鄲36億m3、[HJ2.05mm]廊坊23億m3、石家莊81億m3、保定54億m3、邢臺34億m3、滄州42億m3、衡水28億m3。

減少地下水開采情況分析。經(jīng)過優(yōu)化調(diào)配，與無引江水相比減少地下水開采量266億m3，其中邯鄲26億m3、廊坊23億m3、[HJ2.08mm]石家莊71億m3、保定50億m3、邢臺30億m3、滄州42億m3、衡水24億m3。

水資源供需平衡分析。全省南水北調(diào)受水區(qū)總供水量為443億m3、需水量499億m3，經(jīng)過優(yōu)化調(diào)配，受水區(qū)不同供水水源得到充分利用，置換出（剩余）水量32億m3，（其中地下水21億m3、非 常規(guī)水11億m3，可用于環(huán)境與農(nóng)業(yè)），缺水88億m3，缺水率176%。從行業(yè)來看，生活用水優(yōu)先全部得到滿足，工業(yè)需水量及環(huán)境需水量分別缺水82億m3和缺水06億m3，缺水率分別為292%和136%。從地區(qū)來看，邯鄲缺水04億m3，缺水率59%； 廊坊缺水02億m3，缺水率74%；石家莊缺水01億m3，缺水率09%；保定缺水01億m3，缺水率14%；邢臺市缺水05億m3，缺水率86%；滄州市缺水64億m3，缺水率538%；衡水市缺水11億m3，缺水率256%?？梢?020年保持現(xiàn)狀節(jié)水力度，不足以支撐受水區(qū)的水資源供需平衡。

3.2 2020年供受水區(qū)來水頻率組合50%+50%高效節(jié)水方案調(diào)配結(jié)果[BT）][HJ]

引江水利用情況分析。由圖6和表5可知，在2020年供受水區(qū)來水頻率均為P=50%、城市各業(yè)用水均為高效節(jié)水條件下，經(jīng)過優(yōu)化調(diào)配，受水區(qū)利用引江水量總計298億m3，利用率100%。其中邯鄲27億m3、廊坊22億m3、石家莊71億m3、保定53億m3、邢臺33億m3、滄州69億m3、衡水23億m3。

減少地下水開采情況分析。經(jīng)過優(yōu)化調(diào)配，減少地下水開采[HJ2.1mm]量249億m3，其中邯鄲13億m3、廊坊23億m3、石家莊53億m3、保定47億m3、邢臺25億m3、滄州69億m3、衡水19億m3。

水資源供需平衡分析。全省南水北調(diào)受水區(qū)總供水量為443億m3、需水量406億m3，經(jīng)過優(yōu)化調(diào)配，受水區(qū)不同供水水源得到充分利用，置換出（剩余）水量49億m3（其中地下水38億m3、非常規(guī)水11億m3，可用于環(huán)境與農(nóng)業(yè)），缺水12億m3，缺水率30%。從行業(yè)來看，生活用水優(yōu)先全部得到滿足，工業(yè)需水量及環(huán)境需水量分別缺水09億m3 和缺水03億m3，缺水率分別為44%和77%。從地區(qū)來看，其中邯鄲、廊坊、石家莊、邢臺、衡水

不缺水；保定缺水01億m3，缺水率14%；滄州市缺水11億m3，缺水率120%?？梢娡ㄟ^加大節(jié)水力度，大部分設(shè)區(qū)市可以實現(xiàn)受水區(qū)的水資源供需平衡。

4 結(jié)語

本文以河北省南水北調(diào)中線受水區(qū)為研究區(qū)，以基礎(chǔ)數(shù)據(jù)需求分析為起點，以多目標群決策理論為基礎(chǔ)，以優(yōu)化與模擬為主要技術(shù)，構(gòu)建了基于MIKE BASIN的水資源優(yōu)化調(diào)配模型（WOMMB），解決了多水源、多區(qū)域、多層次、多用戶水資源統(tǒng)一調(diào)配問題，模型可根據(jù)規(guī)劃年不同保證率、不同節(jié)水水平條件、不同用水戶需求，提出綜合效益最優(yōu)的受水區(qū)水資源優(yōu)化調(diào)配方案，為驗證調(diào)配方案的合理性及調(diào)配效果的顯著性，以2020年供受水區(qū)來水頻率組合50%+50%調(diào)配方案為例，采用受水區(qū)總?cè)彼颗c缺水率、生活及工業(yè)供水保證率、置換（剩余）地表水與地下水量、引江水利用量與利用率、減少地下水開采量等指標，對優(yōu)化調(diào)配前后結(jié)果進行對比。結(jié)果顯示，5項指標均較優(yōu)化調(diào)配前效果顯著，各項指標達到最優(yōu)，調(diào)配結(jié)果合理。

（1）缺水量及缺水率經(jīng)優(yōu)化調(diào)配后，基本節(jié)水條件下受水區(qū)缺水量減少了04億m3，缺水率降低08%，為176%；高效節(jié)水條件下受水區(qū)缺水量減少了30億m3，缺水率降低74%，為30%。

（2）生活、工業(yè)供水保證率經(jīng)優(yōu)化調(diào)配后，生活供水保證率達到100%；工業(yè)供水保證率，基本節(jié)水條件下提高14%，達到708%，高效節(jié)水條件下提高150%，達到956%，實現(xiàn)了模型的既定目標。

（3）地表水與地下水置換（剩余）量經(jīng)優(yōu)化調(diào)配后，地表水剩余量進一步得到充分利用，基本節(jié)水和高效節(jié)水條件下，地表水剩余量置換（剩余）量分別從04億m3、21億m3降低至0。地下水置換（剩余）量優(yōu)化調(diào)配前后未發(fā)生變化。

（4）引江水利用量與利用率經(jīng)優(yōu)化調(diào)配后，引江水進一步得到充分利用，引江水利用量達到上限298億m3，利用率達到100%。

（5）減少地下水開采量經(jīng)優(yōu)化調(diào)配后，與無引江水對比，基本節(jié)水和高效節(jié)水條件下，分別較優(yōu)化調(diào)配前減少地下水開采量04億m3和30億m3。

本文通過河北省南水北調(diào)中線受水區(qū)水資源統(tǒng)一調(diào)配方案研究，提出了綜合效益最優(yōu)的受水區(qū)水資源優(yōu)化調(diào)配方案，可實現(xiàn)“引江水高效利用、地表水調(diào)節(jié)利用、地下水控制利用”的高效用水目標，從水資源的開發(fā)利用及保護等方面為水資源管理部門提供決策支持，為實現(xiàn)水資源的可持續(xù)利用及經(jīng)濟社會的可持續(xù)發(fā)展提供保障。

參考文獻（References）：

[1] 高文君.南水北調(diào)中線工程建設(shè)對河南省經(jīng)濟的促進作用[J].法制與社會，2013（9）：214215，224.（GAO W J.Promoting effect of middle route project of southtonorth water diversion on Henan economy[J].Legal System and Society，2013（9）：214215，224.（in Chinese）） DOI：10.19387/j.cnki.10090592.2013.09.114.

[2] 高亮，張玲玲.區(qū)域多水源多用戶水資源優(yōu)化配置研究[J].節(jié)水灌溉，2015（3）：3841.（GAO L，ZHANG L L.Study on Optimal Allocation of Regional Water Resources with Multiple Resources and Multiple Users[J].Water Saving Irrigation，2015（3）：3841.（in Chinese））

[3] 閆寶偉，郭生練，肖義.南水北調(diào)中線水源區(qū)與受水區(qū)降水豐枯遭遇研究[J].水利學報，2007（10）：11781185.（YAN B W，GUO S L，XIAO Y.Synchronousasynchronous encounter probability of richpoor precipitation between water source area and water receiving areas in the Middle Route of SouthtoNorth Water Transfer Project[J].Journal of Hydraulic Engineering，2007（10）：11781185.（in Chinese）） DOI：10.3321/j.issn：05599350.2007.10.005.

[4] 康玲，何小聰，熊其玲.基于貝葉斯網(wǎng)絡(luò)理論的南水北調(diào)中線工程水源區(qū)與受水區(qū)降水豐枯遭遇風險分析[J].水利學報，2010，41（8）：908913.（KANG L，HE X C，XIONG Q L.Risk analysis for precipitation richpoor encounter between source area and receiving area of the Middle Route of SouthtoNorth Water Transfer Project Based on Bayesnet theory[J].Journal of Hydraulic Engineering，2010，41（8）：908913.（in Chinese）） DOI：10.13243/j.cnki.slxb.2010.08.005.

[5] 梁鐘元，王海潮，雷曉輝.氣候變化對南水北調(diào)中線受水區(qū)徑流量影響[J].南水北調(diào)與水利科技，2014，12（2）：137141，145.（LIANG Z Y，WANG H CH，LEI X H.Impact of climate change on the intake area runoff in the MiddleRoute of SouthtoNorth Water Diversion Project[J].SouthtoNorth Water Transfers and Water Science & Technology，2014，12（2）：137141，145.（in Chinese）） DOI：10.13476/j.cnki.nsbdqk.2014.02.032.

[6] 廖衛(wèi)紅，雷曉輝，蔣云鐘.簡易地下水模型EasyGWM在南水北調(diào)中線受水區(qū)的應(yīng)用[J].南水北調(diào)與水利科技，2011，9（4）：812.（LIAO W H，LEI X H，JIANG Y ZH.Application of Easy Groundwater Model in the Reception Area of the Middle Route of SouthtoNorth Water Diversion Project in China[J].SouthtoNorth Water Transfers and Water Science & Technology，2011，9（4）：812.（in Chinese）） DOI：10.3724/SP.J.1201.2011.04008.

[7] 程雙虎，李明良，張俊芝，等.河北省南水北調(diào)受水區(qū)地下水模型研究[J].北京水務(wù)，2016（3）：38.（CHENG SH H，LI M L，ZHANG J Z，et al.Application of FEFLOW in groundwater modeling in intake area of the SouthtoNorth water transfer project in Hebei province[J].Beijing Water，2016（3）：38.（in Chinese）） DOI：10.3969/j.issn.16734637.2016.03.002.

[8] 劉昌明，左建兵.南水北調(diào)中線主要城市節(jié)水潛力分析與對策[J].南水北調(diào)與水利科技，2009，7（1）：17.（LIU CH M，ZUO J B.Water saving potential analysis and countermeasures for major cities for the Midline of SouthtoNorth Water Transfer Project（SNWTP）[J].SouthtoNorth Water Transfers and Water Science & Technology，2009，7（1）：17.（in Chinese）） DOI：10.3969/j.issn.16721683.2009.01.001.

[9] 朱永楠，王慶明，任靜，等.南水北調(diào)受水區(qū)節(jié)水指標體系構(gòu)建及應(yīng)用[J].南水北調(diào)與水利科技，2017，15（6）：187195.（ZHU Y N，WANG Q M，REN J，et al.The construction and application of water saving index system in the SouthtoNorth Water Transfer Project[J].SouthtoNorth Water Transfers and Water Science & Technology，2017，15（6）：187195.（in Chinese）） DOI：10.13476/j.cnki.nsbdqk.2017.06.027.

[10] 耿六成，馮戰(zhàn)洪，趙玲.南水北調(diào)受水區(qū)水資源合理配置研究[J].河北水利水電技術(shù)，2003（1）：4041.（GENG L CH，F(xiàn)ENG ZH H，ZHAO L.Research on rational allocation of water resources in Water Area of SouthtoNorth Water Diversion[J].HeBei Water Resources And Hydropower Engineering，2003（1）：4041.（in Chinese））

[11] 吳澤宇，蔣為群，肖萬格，等.南水北調(diào)中線受水區(qū)水資源配置[J].人民長江，2007（5）：101103.（WU Z Y，JIANG W Q，XIAO W G，et al.Water resources allocation of water service area of Middle Route Project of SN Water Transfer[J].

Yangtze River，2007（5）：101103.（in Chinese）） DOI：10.3969/j.issn.10014179.2007.05.045.

[12] 張平，鄭垂勇.南水北調(diào)中線受水區(qū)水資源優(yōu)化配置研究[J].統(tǒng)計與決策，2006（06）：3841.（ZHANG P，ZHENG CH Y.Optimal allocation of water resources in the receiving area of the Middle Route of the SouthtoNorth Water Transfer Project[J].Statistics and Decision，2006（6）：3841.（in Chinese））

[13] 王海云，高太忠，高京，等.基于AHPLP法的南水北調(diào)中線水資源優(yōu)化配置——以河北省受水區(qū)為例[J].長江流域資源與環(huán)境，2007（5）：588592.（WANG H Y，GAO T ZH，GAO J，et al.Water resources optimal allocation in the Middle Line of SouthToNorth Water Transfer Project using AHPLP——A case of Hebei Province[J].Resources and Environment in the Yangtze Basin，2007（5）：588592.（in Chinese）） DOI：10.3969/j.issn.10048227.2007.05.008.

[14] 盧書超.基于MIKE BASIN的石羊河流域水資源管理模型研究[D].北京：清華大學，2016.（LU S C.Research on simulation model of Water Resources Management Based on MIKE BASIN for Shiyanghe River Basin[D].Beijing：Tsinghua University，2016.（in Chinese））

[15] 趙璧奎.城市原水系統(tǒng)水質(zhì)水量聯(lián)合調(diào)度優(yōu)化方法及應(yīng)用研究[D].北京：華北電力大學，2013.（ZHAO B K.Study on optimization method and application of water quality and quantity joint scheduling of urban raw water system[D].Beijing：North China Electric Power University，2013.（in Chinese））

[16] 畢方明.空間信息處理的P2P分層網(wǎng)絡(luò)模型的研究[D].北京：中國礦業(yè)大學，2010.（BI F M.Research on P2P hierarchical network model of spatial information processing[D].Beijing：China University of Mining，2010.（in Chinese））

[17] 李曉娜.不確定環(huán)境下城市需水量預(yù)測及多水源聯(lián)合供水調(diào)度研究[D].邯鄲：河北工程大學，2018.（LI X N.Prediction of urban water demand under uncertain environment and joint water supply scheduling with multiple water sources[D].Handan：Hebei University of Engineering，2018.（in Chinese））

[18] 龍藝.沱江流域水資源合理配置及開發(fā)利用研究[J].考試周刊，2011（46）：238240.（LONG Y.Study on Rational Allocation and Exploitation of Water Resources in Tuojiang River Basin[J].KaoShiZhouKan，2011（46）：238240.（in Chinese））

[19] 汪黨獻.水資源需求分析理論與方法研究[D].北京：中國水利水電科學研究院，2002.（WANG D X.Study on The theory and method of water resource demand analysis[D].Beijing：China Institute of Water Resources and Hydropower Research，2002.（in Chinese））

[20] 邵東國，賀新春，黃顯峰，等.基于凈效益最大的水資源優(yōu)化配置模型與方法[J].水利學報，2005（9）：10501056.（SHAO D G，HE Q CH，HUANG X F，et al.Optimal water resources deployment model based on maximal net benefit[J].Journal of Hydraulic Engineering，2005（9）：10501056.（in Chinese）） DOI：10.3321/j.issn：05599350.2005.09.005.