王 凡
(錦州市水利事務(wù)服務(wù)中心,遼寧 錦州 121000)
水資源是國(guó)家可持續(xù)發(fā)展和人類賴以生存的基礎(chǔ)資源,對(duì)維持地球生命系統(tǒng)發(fā)揮著重要的作用[1-2]。近年來(lái),隨著經(jīng)濟(jì)的快速發(fā)展和社會(huì)的不斷進(jìn)步,水資源供需矛盾不斷突出,水體污染、水環(huán)境惡化、水資源短缺、地下水超采等問(wèn)題逐漸成為制約社會(huì)發(fā)展的瓶頸[3-6]。水資源優(yōu)化配置是解決區(qū)域水資源短缺問(wèn)題的重要途徑,通過(guò)合理分配區(qū)域內(nèi)各種形式的水資源實(shí)現(xiàn)經(jīng)濟(jì)可持續(xù)發(fā)展和水資源的永續(xù)利用。
目前,國(guó)內(nèi)學(xué)者越來(lái)越關(guān)注水資源優(yōu)化配置問(wèn)題,研究方法也日趨成熟,如劉玒玒等以黑河流域?yàn)槔?,利用蟻群算法進(jìn)行了優(yōu)化配置;潘俊等以沈陽(yáng)市為例,利用遺傳算法對(duì)區(qū)域水資源量進(jìn)行了合理的配置[7-10]。曾萌等以廣東省為例,對(duì)水資源配置模型利用魚(yú)群算法進(jìn)行求解;沙金霞等對(duì)不同規(guī)劃水平年的引大濟(jì)湟工程水資源量,運(yùn)用粒子群算法進(jìn)行優(yōu)化配置[11-12]。在優(yōu)化配置水資源時(shí),這些算法大多存在收斂速度滿、易陷入局部最優(yōu)等缺陷。因此,文章以錦州市為例,以生態(tài)、社會(huì)和經(jīng)濟(jì)效益最優(yōu)為目標(biāo),考慮供水能力、需水能力、污染物質(zhì)量濃度、變量非負(fù)等約束條件構(gòu)建優(yōu)化配置模型,并利用具有計(jì)算精度高、全局尋優(yōu)能力強(qiáng)、收斂速度快等優(yōu)點(diǎn)的鯨魚(yú)優(yōu)化算法,求解出2025年和2030年(不同規(guī)劃年)的配置方案,旨在實(shí)現(xiàn)該區(qū)域經(jīng)濟(jì)社會(huì)的協(xié)調(diào)發(fā)展和水資源合理開(kāi)發(fā)利用。
研究發(fā)現(xiàn),鯨魚(yú)大腦特定區(qū)域內(nèi)的梭形細(xì)胞使鯨魚(yú)區(qū)別于其他生物,即能夠發(fā)生類似于人類的判斷、情感行為。2016年,有學(xué)者提出一種WOA鯨魚(yú)優(yōu)化算法,其理論來(lái)源于座頭鯨的捕食行為,即隨機(jī)捕獵、收縮包圍和位置螺旋式更新等,總體上可以劃分成環(huán)繞式、泡泡網(wǎng)式和搜索式3個(gè)捕魚(yú)階段。
環(huán)繞式捕食時(shí)能夠精準(zhǔn)識(shí)別獵物的位置,并快速完成包圍。將目標(biāo)獵物位置作為當(dāng)前最優(yōu)位置,最優(yōu)位置定義之后其他鯨魚(yú)個(gè)體會(huì)不斷更新其位置,即朝著最優(yōu)位置更替,位置更新可表示成:
(1)
(2)
結(jié)合泡泡網(wǎng)捕食行為能夠建立位置螺旋式更新、收縮包圍兩種數(shù)學(xué)模型,具體如下:
1)螺旋式位置更新機(jī)制。首先,計(jì)算最優(yōu)位置與鯨魚(yú)個(gè)體間的距離,構(gòu)建螺旋數(shù)學(xué)模型模擬座頭鯨的螺旋形運(yùn)動(dòng),其表達(dá)式為:
(3)
(4)
一般地,座頭鯨會(huì)按照螺旋形的路徑收縮包圍獵物,對(duì)此可假設(shè)螺旋式位置更新機(jī)制和收縮包圍機(jī)制的更新概率均為50%,由此構(gòu)建的數(shù)學(xué)模型為:
(5)
式中:p為隨機(jī)數(shù),取值區(qū)間為[-1,1]。
(6)
(7)
錦州市地處“遼西走廊”東端,總面積為10301km2。境內(nèi)地勢(shì)東南低、西北高,從西北向東南依次為低山區(qū)、丘陵區(qū)、平原區(qū),其土地結(jié)構(gòu)整體呈“五山一水四分田”。錦州市屬于溫暖帶半濕潤(rùn)大陸性季風(fēng)氣候,多年平均降水量567mm,年均氣溫7.8℃-9.0℃,四季降水不均,夏季降水占全年的60%-70%。境內(nèi)主要有饒陽(yáng)河、大凌河、女兒河、百股河、觀音洞水、小凌河、二郎洞水、小壩溝等河流水系,主要有流量小變化大、含沙量大、凍期長(zhǎng)、含沙量大、利用率不高、地區(qū)差異顯著等特點(diǎn)。近年來(lái),錦州市憑借其優(yōu)越的資源環(huán)境和地理?xiàng)l件快速發(fā)展經(jīng)濟(jì),各行業(yè)用水需求量和水資源開(kāi)發(fā)利用量急劇增大,致使水資源短缺等問(wèn)題。因此,科學(xué)配置水資源對(duì)于促進(jìn)經(jīng)濟(jì)社會(huì)、水資源、水生態(tài)的協(xié)調(diào)發(fā)展具有重要意義。
錦州市主要有再生水、地下水和地表水(錦凌水庫(kù)、凌河水、遼西北供水工程、其他地表水)等供水水源,預(yù)測(cè)的2025年和2030年可供水量如表1所示。
表1 不同規(guī)劃年的可供水量預(yù)測(cè) 萬(wàn)m3
錦州市用水部門主要有生態(tài)用水、第三產(chǎn)業(yè)用水、工業(yè)用水、農(nóng)業(yè)和生活用水。以2020年為現(xiàn)狀基準(zhǔn)年,對(duì)各用水部門2025年、2030年的需水量利用定額法進(jìn)行預(yù)測(cè),如表2所示。
表2 不同規(guī)劃年的需水量預(yù)測(cè) 萬(wàn)m3
2.3.1 目標(biāo)函數(shù)
1)經(jīng)濟(jì)效益目標(biāo)以區(qū)域最大供水效益來(lái)描述,其函數(shù)表達(dá)式為:
(8)
2)社會(huì)效益目標(biāo)以區(qū)域缺水量最小來(lái)衡量,其函數(shù)表達(dá)式為:
(9)
式中:Djk為子區(qū)k內(nèi)用戶j的需水量,萬(wàn)m3。
3)生態(tài)效益目標(biāo)利用最小化學(xué)需氧量(COD)排放量來(lái)表征,其函數(shù)表達(dá)式為:
(10)
式中:pjk、djk為子區(qū)k中用戶j的污水排放系數(shù)和廢水中COD污染物的濃度,mg/L。
2.3.2 約束條件
1)以子區(qū)k中水源i的總供水量不超過(guò)其可供水量作為供水能力約束條件,即:
(11)
2)以水源向用戶分配的水量處于用需水量上下限之間作為需水能力約束條件,即:
(12)
式中:L(k,j)、H(k,j)為子區(qū)k中用戶j的最大和最小需水量。
3)以子區(qū)k中用戶j排放的重要污染物總量不超過(guò)該區(qū)域最大允許排放量D且廢水中的重要污染物COD濃度在國(guó)家允許排放指標(biāo)d0內(nèi)作為廢水污染物質(zhì)量濃度約束條件,即:
(13)
djk≤d0
(14)
4)以各變量不<0作為非負(fù)約束條件,即:
xijk≥0
(15)
2.3.3 參數(shù)確定
1)權(quán)重系數(shù):采用熵權(quán)和層次分析組合賦權(quán)法確定生態(tài)環(huán)境、社會(huì)、經(jīng)濟(jì)3個(gè)效益目標(biāo)的權(quán)重系數(shù)為0.24、0.31、0.45。
2)供水效益系數(shù):結(jié)合錦州市及其周邊地區(qū)實(shí)際情況和居民生活用水優(yōu)先序最高的原則確定生活供水效益系數(shù)取0.05萬(wàn)元/m3,2025年農(nóng)業(yè)、工業(yè)和第三產(chǎn)業(yè)供水效益系數(shù)為0.0028萬(wàn)元/m3、0.0185萬(wàn)元/m3、0.2012萬(wàn)元/m3,2030年農(nóng)業(yè)、工業(yè)和第三產(chǎn)業(yè)供水效益系數(shù)為0.0047萬(wàn)元/m3、0.0260萬(wàn)元/m3、0.1542萬(wàn)元/m3??紤]到居民生活與城市生態(tài)環(huán)境息息相關(guān),確定生態(tài)環(huán)境供水效益系數(shù)取0.05萬(wàn)元/m3。
3)供水費(fèi)用系數(shù):根據(jù)錦州市稅費(fèi)及水價(jià)征收標(biāo)準(zhǔn),合理確定生活、農(nóng)業(yè)、工業(yè)、第三產(chǎn)業(yè)和生態(tài)環(huán)境用水費(fèi)用系數(shù)為2.45元/m3、0.42元/m3、2.80元/m3、3.50元/m3、2.45元/m3。
4)供水次序系數(shù):與其他水資源相比子區(qū)k水源i的優(yōu)先使用程度用供水次序系數(shù)qik來(lái)衡量,即利用下式計(jì)算供水次數(shù)系數(shù):
(16)
式中:nmax、ni為最大水源和水源i的供水次序序號(hào)。
根據(jù)錦州市各用水部門用水的重要程度和各供水水源的實(shí)際情況,確定水源供水次序?yàn)?再生水<地下水<遼西北供水工程<錦凌水庫(kù)<凌河水<其他地表水,以此計(jì)算的供水次序系數(shù)為0.06、0.12、0.15、0.20、0.25、0.30。
5)用水公平系數(shù):結(jié)合錦州市各用水部門的重要程度確定先后用水次序?yàn)檗r(nóng)業(yè)<第三產(chǎn)業(yè)<工業(yè)<生態(tài)<生活用水,應(yīng)用供水次序系數(shù)公式計(jì)算各部門用水公平系數(shù)依次為0.08、0.12、0.21、0.26、0.33。
6)模型參數(shù)設(shè)置:鯨魚(yú)優(yōu)化算法的最大迭代次數(shù)tmax取1000,種群規(guī)模N為150。
2.3.4 模型求解
在合理確定相關(guān)參數(shù)的基礎(chǔ)上,對(duì)水資源優(yōu)化配置模型利用鯨魚(yú)優(yōu)化算法進(jìn)行求解,主要步驟為:
步驟4:經(jīng)更新計(jì)算,確定鯨魚(yú)個(gè)體的適應(yīng)度值并找到全局最優(yōu)個(gè)體。
步驟5:經(jīng)檢驗(yàn),若滿足終止運(yùn)算條件則輸出最優(yōu)適應(yīng)度值及鯨魚(yú)個(gè)體位置;若不滿足終止條件,令t=t+1返回步驟3再次迭代運(yùn)算。
根據(jù)以上參數(shù)和錦州市水資源實(shí)際情況,利用Matlab軟件求解優(yōu)化配置模型確定不同規(guī)劃年的配置結(jié)果,如表3所示。
表3 不同規(guī)劃年的水資源供需平衡結(jié)果 萬(wàn)m3
1)配水量分析。由表3可知,2025年按配水比例確定錦州市各用水部門排序依次為第三產(chǎn)業(yè)、生態(tài)、生活、農(nóng)業(yè)和工業(yè)配水,該排序與需水量預(yù)測(cè)結(jié)果保持一致。根據(jù)配水比例確定的2025年各用水部門排序與2030年保持相同,該排序可以客觀反映區(qū)域用水情況。結(jié)合配水比例可知,2030年的農(nóng)業(yè)配水量低于2025年,而生態(tài)、生活、工業(yè)和第三產(chǎn)業(yè)配水高于2025年,該配水變化趨勢(shì)與需水變化情況基本相符。
2)缺水量分析。對(duì)比錦州市需水量預(yù)測(cè)結(jié)果和水資源優(yōu)化配置結(jié)果發(fā)現(xiàn),該配置方案符合用水公平和可持續(xù)發(fā)展原則,能夠較好地滿足生態(tài)和生活蓄水。優(yōu)化配置后的缺水情況以生產(chǎn)部門為主,從綜合效益目標(biāo)的角度上,在確保工業(yè)和第三產(chǎn)業(yè)用水需求的條件下允許存在少量的缺水,但考慮到農(nóng)業(yè)用水效率較低的實(shí)際情況,故農(nóng)業(yè)部門出現(xiàn)缺水的概率更高??傮w而言,在實(shí)現(xiàn)生態(tài)、經(jīng)濟(jì)和社會(huì)效益目標(biāo)的情況下,經(jīng)優(yōu)化配置錦州市缺水率不超過(guò)10%,可見(jiàn)相關(guān)參數(shù)設(shè)置和優(yōu)化配置模型是科學(xué)合理的。
文章錦州市為例,利用相關(guān)約束條件和生態(tài)環(huán)境、社會(huì)、經(jīng)濟(jì)目標(biāo)構(gòu)建優(yōu)化配置模型,并利用鯨魚(yú)優(yōu)化算法求解模型。結(jié)果發(fā)現(xiàn),錦州市水資源量能夠滿足2025年和2030年的生態(tài)需水、生活需水,生產(chǎn)用水部門缺水率不超過(guò)10%,缺水量較少,該配置方案科學(xué)合理能夠滿足未來(lái)用水需求,為錦州市水資源可持續(xù)利用和科學(xué)管理提供一定指導(dǎo)。在水資源優(yōu)化配置問(wèn)題求解時(shí)鯨魚(yú)優(yōu)化算法具有較好的有效性、適應(yīng)性,能夠有效解決傳統(tǒng)算法存在的問(wèn)題,在優(yōu)化配置水資源等領(lǐng)域具有一定應(yīng)用價(jià)值。