水資源污染負荷強度的灰色有效性預測模型及其應用

王 晗，張 峰，薛惠鋒

(1.中國航天系統(tǒng)科學與工程研究院，北京 100048；2.山東理工大學 管理學院，山東 淄博 255012)

隨著全面建成小康社會步伐的推進，社會居民的生活質(zhì)量持續(xù)提升，對高質(zhì)量的水資源需求程度也愈發(fā)明顯，而基于中國過去長期水資源利用效率偏低、水體污染治理投入不足、水資源總量及其時空分布不均等系列限制性因素，現(xiàn)有水質(zhì)狀況已難以滿足水資源可持續(xù)開發(fā)利用的實際需求。對此，國家不斷完善水污染防治制度，試圖構(gòu)建包括水資源污染負荷在內(nèi)的更加精準的水環(huán)境質(zhì)量及污染物排放監(jiān)測體系。

鑒于水資源污染負荷預測在水污染防治工作中的重要性，學者們對其研究的思路也逐步拓寬，縱觀目前關(guān)于水資源污染負荷的預測方法則主要集中于線性回歸、BP神經(jīng)網(wǎng)絡和系統(tǒng)動力仿真等。Sehgal等[1]利用線性擬合的方式對水污染指標趨勢進行預判，經(jīng)過其實例驗證發(fā)現(xiàn)基于線性回歸的方程擬合雖然操作簡便，但無法適應水資源污染負荷多指標精確預測的需求。相似的觀點還出現(xiàn)在Mouatadid等[2]對城市水量需求及水質(zhì)變化的預測中。Djerbouai等[3]采取神經(jīng)網(wǎng)絡的方法構(gòu)建水資源污染排放量的非線性擬合模型，結(jié)果發(fā)現(xiàn)隨著神經(jīng)網(wǎng)絡訓練次數(shù)的增加，其誤差參數(shù)控制效果呈非確定性，尤其是容易陷入局部極值。張巖祥等[4]引入層次分析與模糊函數(shù)建立水質(zhì)綜合評價模型，認為對于水質(zhì)預測的關(guān)鍵在于評價指標權(quán)系數(shù)的合理確定。張峰等[5]基于復合系統(tǒng)的視角，利用系統(tǒng)動力學仿真方法檢驗并預測水資源與水環(huán)境復雜影響因素的變化趨勢，操作過程較為繁雜。由上可見，現(xiàn)階段應用在水量與水質(zhì)方面的傳統(tǒng)預測分析方法多難以滿足水資源污染負荷預測的高精度要求。近年來，組合預測方法在多領(lǐng)域中開始應用[6]，并體現(xiàn)出較高預測精度，但根據(jù)相關(guān)領(lǐng)域的現(xiàn)有研究成果，不難發(fā)現(xiàn)構(gòu)成其組合預測模型的各類單個模型的類別、方法及數(shù)量還需要進行多維度的考察。而鄧聚龍教授提出的灰色系統(tǒng)理論，為分析信息不完備系統(tǒng)提供了一種良好的建模方法，其包含的灰色矩陣在處理有限樣本的集成預測問題上已被學者們驗證了其具有低誤差的建模優(yōu)勢[7]，但在水資源污染負荷預測方面的應用成效還需要深入研討。

綜上，本文基于水資源污染負荷強度預測的現(xiàn)實需求，選取工業(yè)單位產(chǎn)值化學需氧量排放量為表征指標，利用灰色系統(tǒng)理論中GM(1,1)模型、Verhulst模型和SCGM(1,1)c模型的優(yōu)點及有效度測算方法，分別建立水資源污染負荷強度的單預測模型及GM-Verhulst-SCGM組合預測模型，并在驗證其誤差可接受的基礎(chǔ)上進一步對水資源污染負荷強度的演變態(tài)勢進行預測分析。

1 模型構(gòu)建

利用灰色系統(tǒng)理論的白化函數(shù)、累積運算等解決水資源污染負荷強度分析中易出現(xiàn)的隨機變化要素和數(shù)據(jù)樣本規(guī)模的非一致性問題[8]，具體過程如下。

1.1 GM(1,1)模型

GM(1,1)模型是現(xiàn)階段灰色系統(tǒng)理論建模應用中較為廣泛與經(jīng)典的模型之一。其基本的建模理論與實現(xiàn)過程如下：

(1)

式中：a、ε是待識參數(shù)。

定義2：參數(shù)向量為∏=[au]，同時：

Yn=[x(0)(2),x(0)(3),…,x(0)(N)]，及

利用Least square方法，計算∏=(ΓTΓ)-1(ΓTYn)，得到x(1)(t)白化微分方程的解：

(2)

要取得原數(shù)據(jù)序列的預測序列，還需要對以上式子進行逆向累減計算，即：

(3)

1.2 Verhulst模型

Verhulst模型常用于描述分析具有“S”飽和生物特性的增長型模型，該方法目前應用所涉領(lǐng)域相對較多，包括自然資源開發(fā)消耗、企業(yè)成本控制、生態(tài)治理投資等諸多方面[9]。本文研究發(fā)現(xiàn)水資源污染負荷強度滿足“S”型增長的特性，可采用該模型對水資源污染負荷強度的趨勢預測。

利用Least square取得Φ=(ZTZ)-1(ZTY)，描述Verhulst模型為：

(4)

1.3 灰色系統(tǒng)云SCGM(1,1)c模型

SCGM(1,1)c模型是依據(jù)系統(tǒng)云理論對GM(1,1)模型的發(fā)展而得，其關(guān)鍵的計算環(huán)節(jié)是積分變換及數(shù)據(jù)關(guān)聯(lián)計算。在該模型的基礎(chǔ)上，進一步提煉出單因子灰色系統(tǒng)云SCGM(1,h)c，當h=1即為SCGM(1,1)c模型。此模型的優(yōu)勢在于能夠積分累計計算最大限度地挖掘數(shù)據(jù)內(nèi)在信息，同時相比之下，其建模所需的最小樣本數(shù)據(jù)、計算冗余度等相對較低等[10]。模型構(gòu)建過程如下：

定義4：水資源污染負荷強度的原序列描述為：X(0)={X(0)(1),X(0)(2),…,X(0)(n)}。其中，n指序列長度。

(5)

其中

(6)

(7)

模型(7)單次響應函數(shù)為：

(8)

其中：

(9)

(10)

(11)

(12)

1.4 基于預測有效度的GM-Verhulst-SCGM預測模型

水資源污染負荷強度的GM-Verhulst-SCGM組合預測模型的總體形式可以描述為：

(13)

(14)

Si=1-δi

(15)

(16)

1.5 水資源污染負荷強度預測模型精度評價指標

2 水資源污染負荷強度實證預測

相比傳統(tǒng)的線性回歸擬合預測方法，灰色預測模型具有對于訓練樣本的需求量相對較少及精度高等特點，本文為驗證 組合預測模型在水資源污染負荷強度實際預測的可行性與有效性，以工業(yè)單位產(chǎn)值化學需氧量排放量2004-2016年歷史數(shù)據(jù)為算例樣本，并分別運用 模型、 模型和 模型對工業(yè)單位產(chǎn)值化學需氧量排放量指標進行計算，同時觀測各預測模型下其參數(shù)的變化走向，數(shù)據(jù)主要源于《中國統(tǒng)計年鑒》(2005-2017年)、《中國水資源統(tǒng)計公報》(2004-2017年)。在上述基礎(chǔ)上，綜合運用水資源污染負荷強度預測模型精度評價指標將各水資源污染負荷強度預測模型計算結(jié)果與 組合預測模型進行對比，其中設置模型擬合數(shù)據(jù)區(qū)間為2004-2013年，模型驗證數(shù)據(jù)區(qū)間為2014-2016年。

2.1 單模型預測

(1)GM(1,1)模型。

(2)Verhulst模型。

(3)SCGM(1,1)c模型。

利用各模型對2004-2013年工業(yè)單位產(chǎn)值化學需氧量排放量樣本數(shù)據(jù)進行線性擬合，按照統(tǒng)計指標對其分別進行評價，見表1。

2.2 GM-Verhulst-SCGM組合模型預測

將單工業(yè)單位產(chǎn)值化學需氧量排放量預測值代入上述公式，取得水資源污染負荷強度組合預測值，結(jié)果見表1，其誤差對比情況見圖1所示。

表1 各模型預測擬合結(jié)果

表2 各模型預測結(jié)果

圖1 模型預測誤差值對比

2.3 水資源污染負荷強度預測效果評價

按照各水資源污染負荷強度單個及組合預測模型，分別對工業(yè)單位產(chǎn)值化學需氧量排放2014-2016年時序數(shù)據(jù)進行動態(tài)估計，各預測模型的預估效果利用MRE、MAE雙項評價指標進行對比檢驗，以確保其準確性，結(jié)果見表2。按照MRE與MAE的統(tǒng)計值，可發(fā)現(xiàn)在各單工業(yè)單位產(chǎn)值化學需氧量排放量預測模型預估中，GM(1,1)模型呈現(xiàn)出了較小的統(tǒng)計指標值，為單模型預測中預測效果較佳模型。其次預測效果由優(yōu)至劣分別為Verhulst模型、SCGM(1,1)c模型。而在整體效測度果上，GM-Verhulst-SCGM組合模型盡管在MRE值(0.001 430)上要略高于GM(1,1)模型(0.012 933)，但是均低于Verhulst模型和SCGM(1,1)c模型(其測算值依次為0.026 938、0.036 396)，同時，GM-Verhulst-SCGM組合模型的MAE值為各組模型中的最低值(0.002 233)。說明GM-Verhulst-SCGM組合預測的誤差落于預測可控范圍之內(nèi)，更能夠逼近工業(yè)單位產(chǎn)值化學需氧量排放量的實際值，同比各個單工業(yè)單位產(chǎn)值化學需氧量排放量預測模型的預測效果更具針對性，且精度實現(xiàn)了有效提升。

2.4 水資源污染負荷強度預測及分析

通過對工業(yè)單位產(chǎn)值化學需氧量排放量歷史數(shù)據(jù)的擬合預測驗證，可知其灰色組合預測模型在應用于水資源污染負荷強度預測中具有較高的實效性。據(jù)此，本文選取GM-Verhulst-SCGM組合模型預估工業(yè)單位產(chǎn)值化學需氧量排放量2017-2021年期間的變化趨勢，見圖2。

圖2 2017-2020年排污強度趨勢預測結(jié)果

根據(jù)預測結(jié)果可知，工業(yè)單位產(chǎn)值化學需氧量排放量在未來一段時間內(nèi)的變化整體上呈穩(wěn)步下降的基本態(tài)勢，且從排放強度的下降速率來看，2017-2020年期間其落幅相對顯著，從1.281 4 萬t/千億元下降到1.191 1 萬t/千億元，這說明近年來關(guān)于推進水污染防治的相關(guān)政策措施相對有效，在控制水污染強度方面可取得良好的成效。但從2020-2021年期間的工業(yè)單位產(chǎn)值化學需氧量排放量變化情況來看，其下降速率出現(xiàn)了明顯的放緩態(tài)勢，由此形成了其排放強度的“拐點”，該現(xiàn)象在較大程度上可能是受前期水污染治理強度的影響，即在“高投入、高污染、高排放”的傳統(tǒng)發(fā)展模式下，水污染程度愈發(fā)突出，而近年來通過加快對重點污染源等進行水污染調(diào)控倒逼，促使水污染治理總體成效上出現(xiàn)了顯著性的扭轉(zhuǎn)，但由此在一定程度上也削減了后期所需要水污染治理規(guī)模，并逐步形成水資源污染負荷強度控制的“換擋期”，由高減速向穩(wěn)步趨緩轉(zhuǎn)變，使其表現(xiàn)形態(tài)上出現(xiàn)了“拐點”，但這也表征著水污染防治由“淺水區(qū)”向“深水區(qū)”的發(fā)展，其深層次治理的質(zhì)量要求更為嚴格。因此，在加強重點污染源的整治倒逼的同時，也要逐步提高對中輕度污染源的治理引導與水質(zhì)“紅線”監(jiān)管力度，并根據(jù)水污染防治各階段出現(xiàn)的實際問題，及時通過采取制度約束、技術(shù)推進、社會激勵等相關(guān)綜合性措施進行系統(tǒng)調(diào)控。

3 結(jié)束語

基于水資源污染負荷強度的預測需求，本文以2004-2016年期間工業(yè)單位產(chǎn)值化學需氧量排放量歷史數(shù)據(jù)為樣本，利用灰色系統(tǒng)理論構(gòu)建其 模型、 模型和 模型，并在此基礎(chǔ)上選取預測有效度方法建立其灰色 組合預測模型。通過實證研究發(fā)現(xiàn)，其灰色組合預測模型可集成各單預測模型的優(yōu)勢，具有更加理想的預測效果，能夠滿足水資源污染負荷強度的高精度預測需求。而通過對水資源污染負荷強度未來趨勢的預測，發(fā)現(xiàn)在未來一段時間內(nèi)水資源污染負荷強度整體上呈現(xiàn)為持續(xù)下降的態(tài)勢，印證了現(xiàn)階段水污染防治措施有效性與可持續(xù)性，而在2020年存在出現(xiàn)水資源污染負荷強度退落的“拐點”，即由高減速向穩(wěn)步趨緩轉(zhuǎn)變的“換擋期”，但這也預示著水污染防治將由“淺水區(qū)”向“深水區(qū)”的轉(zhuǎn)變?？傮w而言，灰色組合預測模型雖在計算復雜度上較單預測模型有一定提升，但可利用MATLAB等現(xiàn)有計算手段較易實現(xiàn)，具有較強的可操作性與靈活性。而通過對水資源污染負荷強度的預測驗證，也可為相關(guān)資源領(lǐng)域的預測分析提供方法借鑒。