陶春梅，王建成，溫江淯，崔浩林

(1.江蘇省連云港市水利局,江蘇 連云港 222000；2.南京市水利規(guī)劃設(shè)計(jì)院股份有限公司連云港分公司,江蘇 連云港 222000；3.連云港市市區(qū)水工程管理處,江蘇 連云港 222000;4.連云港市臨洪工程管理處,江蘇 連云港 222000)

當(dāng)前城市發(fā)展迅猛，而排污系統(tǒng)發(fā)展卻相對(duì)滯后，導(dǎo)致湖氮、磷和其他污染物大幅增加，周邊水庫水質(zhì)急劇下降，生物多樣性下降，嚴(yán)重影響生態(tài)環(huán)境和居民生活健康[1- 3]，因此對(duì)此類生態(tài)環(huán)境展開生態(tài)修復(fù)并展開評(píng)價(jià)具有重要意義。當(dāng)前國(guó)內(nèi)外對(duì)次展開了相關(guān)研究[4]。

趙金召等[5]運(yùn)用層次分析法，以復(fù)綠植被生態(tài)的適應(yīng)性、綠化效果、經(jīng)濟(jì)性與應(yīng)用潛力為準(zhǔn)則層，選取11個(gè)指標(biāo)，構(gòu)建高陡巖質(zhì)邊坡礦山生態(tài)修復(fù)綠化效果評(píng)價(jià)體系;李舟等[6]用層次分析法對(duì)銀西高鐵寧夏段干旱區(qū)的植物物種進(jìn)行選擇和評(píng)價(jià)技術(shù)研究;史帥航等[7]以西南地區(qū)某礦產(chǎn)集采區(qū)為例，探索區(qū)內(nèi)土壤重金屬含量及其隨流域遷移規(guī)律，分析區(qū)內(nèi)土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)；王瑾[8]通過分析調(diào)查區(qū)生態(tài)環(huán)境現(xiàn)狀，建立了相關(guān)生態(tài)保護(hù)修復(fù)體制。此外，還有部分學(xué)者基于組合權(quán)重法和定量評(píng)價(jià)方法對(duì)生態(tài)評(píng)價(jià)系統(tǒng)進(jìn)行研究。

從以上文獻(xiàn)分析可以得出，當(dāng)前的生態(tài)評(píng)價(jià)方法大多都是基于復(fù)雜繁瑣的數(shù)據(jù)計(jì)算，在評(píng)價(jià)效率方面十分不利，因此，本文開發(fā)了一種生態(tài)修復(fù)效果自適應(yīng)評(píng)價(jià)算法，以石梁河水庫生態(tài)修復(fù)工程為例，選取了高錳酸鹽，氨氮和總磷三個(gè)分類指標(biāo)，采用決策樹自適應(yīng)評(píng)價(jià)算法，以GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》和修正的Carlson營(yíng)養(yǎng)狀態(tài)指數(shù)為評(píng)價(jià)閾值對(duì)生態(tài)修復(fù)后的石梁河水庫進(jìn)行自動(dòng)評(píng)價(jià)分析(2016—2019)。研究成果可為相關(guān)工程提供參考。

1 工程概況

石梁河水庫90年代曾被水利部列為第二批全國(guó)重點(diǎn)病險(xiǎn)水庫之一，1958—1962年，石梁河水庫建成，樞紐工程包括主壩、副壩各一座，設(shè)計(jì)流量4000m3/s泄洪閘一座，1120千瓦水電站一座，輸水涵洞(閘)6座。石梁河水庫建成后由防洪、灌溉功能演變?yōu)槎喙δ芫C合性水庫。根據(jù)淮委規(guī)劃設(shè)計(jì)研究院1998年提出的《對(duì)石梁河新、老泄洪閘泄量分配的建議》，在正常年和50年一遇洪水、閘上水位分別是24.0和26.08m時(shí)，老閘泄量均為2500m3/s，新閘泄量分別為2500、3500m3/s;在100年一遇洪水、閘上水位26.81m時(shí)，老閘泄量3000m3/s，新閘泄量為3500m3/s;在2000年一遇洪水、閘上水位27.95m時(shí)，老閘泄量5000m3/s，新閘泄量為5131m3/s。

2 生態(tài)修復(fù)措施

2.1 生態(tài)田埂工程

石梁河水庫(又名海陵湖)位于江蘇省東海縣石梁河鎮(zhèn)北側(cè)，本次采用生態(tài)田埂、果園生草技術(shù)，石梁河水庫周邊水系如圖1所示。本次田埂工程主要是改變田埂高程：即加高10～15cm，從而將與降雨相關(guān)的地表徑流控制30～50cm，進(jìn)而從根源上有效防止污染的農(nóng)田通過大部分中小降雨流入水庫。此外，選取部分植物進(jìn)行密集種植，分布在田埂的兩側(cè)，形成隔離帶，這一舉措能夠有效阻截氮、磷養(yǎng)分損失，在產(chǎn)生地表徑流時(shí)可以對(duì)殘留農(nóng)藥進(jìn)行固存，防止向水體遷移。

圖1 石梁河水庫周邊水系圖

2.2 生態(tài)溝渠

田間溝渠同時(shí)具備灌溉輸水和田間排水兩種功能。規(guī)劃基于生態(tài)學(xué)原理，將原有的土質(zhì)或三面光水泥溝渠改造成生態(tài)攔截溝渠。生態(tài)溝渠主要在田塊間進(jìn)行交錯(cuò)分布，由3部分組成：初沉池，泥硬質(zhì)生態(tài)溝框架和生態(tài)植物。其中第一部分初沉池用于沉淀由于徑流帶來的農(nóng)田顆粒物，并進(jìn)行回收處理。初沉池通常安裝在農(nóng)田排水出口處，并連接至生態(tài)溝渠處。第二和第三部分(溝底、溝板)為水泥硬質(zhì)板，并用鉆機(jī)進(jìn)行打孔，用于種植植物。該部分作用主要是形成過濾系統(tǒng)，滯留農(nóng)田徑流污染物，截流降雨形成的徑流以及滲漏水中的氮、磷。

2.3 生態(tài)塘

規(guī)劃區(qū)內(nèi)的部分溝塘被村民用作魚塘或種植水生經(jīng)濟(jì)作物，還有很多坑塘沒有得到充分利用，有些甚至被作為養(yǎng)殖場(chǎng)動(dòng)物糞便、生活垃圾的傾倒場(chǎng)所，嚴(yán)重影響水生態(tài)健康。規(guī)劃采用多塘系統(tǒng)改善水生態(tài)環(huán)境，多塘系統(tǒng)是指由多個(gè)沉降塘或滯留塘組成并共同發(fā)揮生態(tài)凈化功效的系統(tǒng)工程。沉降塘位于系統(tǒng)的前端，滯留塘位于系統(tǒng)的后端，沉降塘深度要大于滯留塘。多塘系統(tǒng)平面布置示意圖如圖2所示。

圖2 多塘系統(tǒng)平面布置示意圖

3 生態(tài)修復(fù)評(píng)價(jià)指標(biāo)

3.1 水質(zhì)質(zhì)量

水質(zhì)質(zhì)量是生態(tài)修復(fù)結(jié)果的中套評(píng)價(jià)指標(biāo)，本規(guī)劃主要以高錳酸鹽指數(shù)指數(shù)、氨氮、總磷作為評(píng)價(jià)項(xiàng)目，以GB 3838—2002為評(píng)價(jià)標(biāo)準(zhǔn)，水質(zhì)評(píng)價(jià)方法采用單因子指數(shù)法評(píng)價(jià)。石梁河水庫及入庫河流水質(zhì)評(píng)價(jià)采用Ⅲ類標(biāo)準(zhǔn)，標(biāo)準(zhǔn)限值見表1。

表1 水質(zhì)評(píng)價(jià)標(biāo)準(zhǔn)限值

3.2 Carlson營(yíng)養(yǎng)狀態(tài)指數(shù)

由于造成庫區(qū)富營(yíng)養(yǎng)化的因素有很多，本文選取修正的Carlson營(yíng)養(yǎng)狀態(tài)指數(shù)對(duì)2016—2019年的綜合營(yíng)養(yǎng)狀態(tài)指數(shù)進(jìn)行比較分析，分級(jí)標(biāo)準(zhǔn)見表2。

表2 湖泊(水庫)營(yíng)養(yǎng)狀態(tài)分級(jí)方法

4 自適應(yīng)評(píng)價(jià)算法

本文采用的自適應(yīng)評(píng)價(jià)算法為機(jī)器學(xué)習(xí)當(dāng)中的決策樹系統(tǒng)，能夠自動(dòng)根據(jù)輸入數(shù)據(jù)和相應(yīng)的閾值門檻進(jìn)行評(píng)價(jià)分析。決策樹模型是監(jiān)督學(xué)習(xí)算法中比較經(jīng)典的模型，會(huì)學(xué)習(xí)并回答將樣本分成不同類別的一系列問題。每個(gè)問題都有一個(gè)二元答案，根據(jù)滿足的條件，樣本將被分配給樹的左分支或右分支，并且分支中還可以包含分支；一旦學(xué)習(xí)了模型，就可以用圖形表示一棵樹的形狀?；跇涞哪Ｐ图瓤梢杂糜诜诸惾蝿?wù)，也可以用于回歸任務(wù)，因此可將它們描述為分類和回歸樹(CART)。其基本算法如圖3所示。

圖3 決策樹基本算法

在樹的構(gòu)建過程中的每個(gè)階段，決策樹算法都會(huì)考慮所有輸入變量，并選擇類別之間可分性較高的變量。該算法將從根節(jié)點(diǎn)開始，在每個(gè)分支上再次查找最適合區(qū)別類別的樣本特征并選作當(dāng)前分支。但是，決定每次劃分的最優(yōu)特征可通過兩種不同的方法來實(shí)現(xiàn)：熵的差異(稱為信息增益)和基尼指數(shù)的差異(稱為基尼增益)。這兩種方法得到的結(jié)果通常非常相似；但基尼指數(shù)(以社會(huì)學(xué)家和統(tǒng)計(jì)學(xué)家的計(jì)算速度要快一些，所以本文采用基尼增益劃分樹。

5 評(píng)價(jià)算法應(yīng)用結(jié)果

如圖4所示2016—2019年石梁河水庫水質(zhì)監(jiān)測(cè)數(shù)據(jù)，即本次分析的輸入數(shù)據(jù)集，輸入的分類指標(biāo)包括高錳酸鹽，氨氮和總磷，分析時(shí)只需將監(jiān)測(cè)數(shù)據(jù)文件轉(zhuǎn)為CSV.格式輸入寫好的程序即可，評(píng)價(jià)分類閾值見表1，水質(zhì)評(píng)價(jià)結(jié)果見表2。根據(jù)輸入數(shù)據(jù)和表2中結(jié)果，石梁河水庫高錳酸鹽指數(shù)濃度和平均濃度分別為2.75～5.20mg/L和4.00mg/L，100.0%達(dá)標(biāo)；氨氮濃度最小0.05mg/L，最大0.71mg/L，平均濃度為0.34mg/L，達(dá)標(biāo)率為100.0%，主要超標(biāo)因子為總磷(0.13)?？偭诐舛葹樽钚?.02mg/L，最大0.18mg/L，最大超標(biāo)倍數(shù)2.50倍，平均濃度為0.06mg/L，達(dá)標(biāo)率為70.7%。

圖4 2016—2019年石梁河水庫水質(zhì)監(jiān)測(cè)數(shù)據(jù)

2016年石梁河水庫水質(zhì)良好，符合地表水Ⅲ類標(biāo)準(zhǔn)。2017年及2018年水質(zhì)符合地表水Ⅳ類標(biāo)準(zhǔn)，為輕度污染，主要超標(biāo)因子為總磷。2019年石梁河水庫水質(zhì)良好，符合地表水Ⅲ類標(biāo)準(zhǔn)。石梁河水庫綜合營(yíng)養(yǎng)狀態(tài)指數(shù)評(píng)價(jià)結(jié)果見表3。由富營(yíng)養(yǎng)化評(píng)價(jià)結(jié)果可知，石梁河水庫2016—2019年綜合營(yíng)養(yǎng)狀態(tài)指數(shù)在富營(yíng)養(yǎng)限值左右波動(dòng)，2016年及2017年為中營(yíng)養(yǎng)，2018年為輕度富營(yíng)養(yǎng)，2019年為中營(yíng)養(yǎng)。由以上評(píng)價(jià)結(jié)果可知，本文采用的決策樹評(píng)價(jià)算法可以直接根據(jù)所選的生態(tài)監(jiān)測(cè)指標(biāo)，用程序直接分析所給閾值后的評(píng)價(jià)類型，避免了人工操作的繁瑣，值得應(yīng)用推廣，部分代碼接口如圖5所示。

表2 石梁河水庫2016—2019年水質(zhì)評(píng)價(jià)結(jié)果

表3 石梁河水庫綜合營(yíng)養(yǎng)狀態(tài)指數(shù)評(píng)價(jià)結(jié)果

圖5 部分代碼接口

6 結(jié)語

當(dāng)前的生態(tài)修復(fù)評(píng)價(jià)過程復(fù)雜繁瑣，在評(píng)價(jià)效率方面十分不利。本文為提高評(píng)價(jià)效率，以石梁河水庫生態(tài)修復(fù)工程為例，選取了3個(gè)分類指標(biāo)，采用決策樹自適應(yīng)評(píng)價(jià)算法對(duì)生態(tài)修復(fù)后的石梁河水庫進(jìn)行自動(dòng)評(píng)價(jià)分析(2016—2019)。結(jié)果表明，2016年石梁河水庫符合地表水Ⅲ類標(biāo)準(zhǔn)，2017年及2018年水質(zhì)符合地表水Ⅳ類標(biāo)準(zhǔn)，2019年符合地表水Ⅲ類標(biāo)準(zhǔn)。此外，石梁河水庫2016年及2017年為中營(yíng)養(yǎng)，2018年為輕度富營(yíng)養(yǎng)，2019年為中營(yíng)養(yǎng)，自動(dòng)評(píng)價(jià)結(jié)果與公式計(jì)算結(jié)果一致。因此，本文采用的評(píng)價(jià)算法可以直接根據(jù)生態(tài)監(jiān)測(cè)指標(biāo)實(shí)現(xiàn)了評(píng)價(jià)自動(dòng)化，準(zhǔn)確度高，值得應(yīng)用推廣。