王榮森，吳 蓉

(中國(guó)科學(xué)技術(shù)大學(xué) 管理學(xué)院，合肥 230026)

1 水環(huán)境賬戶計(jì)算模型

本研究采取2013年張義等人改進(jìn)模型的水污染生態(tài)足跡計(jì)算方法，舍棄水資源產(chǎn)量因子，不同的是，本研究單獨(dú)考慮水環(huán)境狀況，并在其總氮、總磷、化學(xué)需氧量三大污染物的基礎(chǔ)上，添加氨氮污染物，具體方法如下：

EFT=max(EFN，EFP，EFCOD，EFNH3N)

EFN=CN/PN

EFP=CP/PP

(1)

EFCOD=CCOD/PCOD

EFNH3N=CNH3N/PNH3N

上式中，EFT為總的水污染生態(tài)足跡(hm2)，EFN為總氮污染水生態(tài)足跡，EFP為總磷污染水生態(tài)足跡，EFCOD為有機(jī)物水生態(tài)足跡，EFNH3N為氨氮污染水生態(tài)足跡，CN、CP、CCOD、CNH3N分別為廢水中總氮、總磷、有機(jī)物和氨氮的排放量(t)；PN、PP、PCOD、PNH3N為水域消納總氮、總磷、有機(jī)物和氨氮的全球平均能力(t/hm2)，由于幾種污染物在環(huán)境影響中的重疊作用，故取其中的最大值作為最終的水污染生態(tài)足跡。以下是承載力模型：

CCT=min(CCN,CCP,CCCOD,CCNH3N)

CCN=0.88×(QWP×UN/PN)

CCP=0.88×(QWP×UP/PP)

(2)

CCCOD=0.88×(QWP×UCOD/PCOD)

CCNH3N=0.88×(QWP×UNH3N/PNH3N)

QWP=QW-(QR×K)

其中，CCT為總的水環(huán)境承載力(hm2)，CCN、CCP、CCCOD、CCNH3N依次是氮、磷、有機(jī)物、氨氮水環(huán)境承載力；QWP為研究區(qū)用于消納污染物的水量(m3)，UN、UP、UCOD、UNH3N為Ⅲ類水質(zhì)標(biāo)椎中基本能維持水體生態(tài)服務(wù)功能的氮、磷、有機(jī)物和氨氮含量上限(mg/L)；承載力中須扣除12%用于生物多樣性保護(hù)，故有0.88的系數(shù)[1]；QW是水資源總量(m3)，QR為地區(qū)取水量(m3)，K為地區(qū)綜合耗水率(可由各省水資源公報(bào)得到)；由于任一污染物超標(biāo)都會(huì)造成水體生態(tài)系統(tǒng)功能的破壞，所以最終的水環(huán)境生態(tài)承載力取氮、磷、有機(jī)物和氨氮水環(huán)境承載力中的最小值；水環(huán)境賬戶盈余ES為水環(huán)境承載力與水污染生態(tài)足跡之差；水環(huán)境生態(tài)壓力指數(shù)IP為水污染生態(tài)足跡與水環(huán)境承載力之比，如下式：

ES=CCT-EFT

IP=EFT/CCT

(3)

2 數(shù)據(jù)來(lái)源

Ⅲ類水質(zhì)標(biāo)椎中基本能維持水體生態(tài)服務(wù)功能的氮、磷、有機(jī)物和氨氮含量上限分別為1.0 mg/L、0.2 mg/L、20 mg/L和1.0 mg/L[2]，水資源全球平均生產(chǎn)能力為3 140 m3/hm2，通過(guò)單位換算得到水域消納總氮、總磷、有機(jī)物和氨氮的全球平均能力分別為0.003 145 t/hm2、0.000 629 t/hm2、0.062 893 t/hm2和 0.003 145 t/hm2，其他數(shù)據(jù)主要來(lái)源于2013-2017年的中國(guó)統(tǒng)計(jì)年鑒、各省份統(tǒng)計(jì)年鑒、水資源公報(bào)，部分?jǐn)?shù)據(jù)根據(jù)原始數(shù)據(jù)計(jì)算得來(lái)。

3 全國(guó)水環(huán)境賬戶情況

根據(jù)全國(guó)除西藏和港澳臺(tái)以外的30個(gè)省份的水污染生態(tài)足跡的計(jì)算結(jié)果，得到2013-2017年的各省份水環(huán)境盈余/赤字情況，由于大部分地區(qū)處于水環(huán)境赤字狀況，所以僅列出盈余省份，具體情況如下：

2013年：青海、貴州、廣西、云南、江西、海南。

2014年：青海、貴州、廣西、云南、江西、福建。

2015年：青海、貴州、廣西、云南、江西、浙江、福建。

2016年：新疆、青海、貴州、四川、廣西、云南、吉林、安徽、江西、湖北、湖南、福建、海南。

2017年：新疆、青海、四川、貴州、廣西、云南、江西、湖南、福建、海南。

除上述所列省份之外，其余省份均處于赤字狀態(tài)。2013-2017年，只有青海、貴州、廣西、云南、江西5省在各年均處于水環(huán)境盈余狀態(tài)，可見(jiàn)我國(guó)水環(huán)境污染形勢(shì)十分嚴(yán)峻。2015年之后，處于水環(huán)境盈余的省份明顯增加，說(shuō)明近兩年地方水污染狀況得到一定程度的改善。根據(jù)水污染生態(tài)足跡、水環(huán)境承載力數(shù)據(jù)和公式(3)，分別得到西、中、東部地區(qū)的人均水環(huán)境壓力指數(shù)平均值，如圖1所示：

圖1 分地區(qū)人均水環(huán)境生態(tài)壓力指數(shù)Fig.1 Ecological pressure index of water environment per capita by region

從圖1可見(jiàn)，東部人均水環(huán)境生態(tài)壓力最大，中部次之，西部最小。東部地區(qū)2013-2014年壓力指數(shù)由20.91陡然增加至84.47，2014年之后又迅速回落，2016年下降到最低點(diǎn)；中部和東部地區(qū)變化較為平穩(wěn)，2016年之前中部的人均壓力指數(shù)均位于西部地區(qū)之上，2017年兩者趨于一致，西、中部地區(qū)在這5年間的人均壓力指數(shù)呈平穩(wěn)波動(dòng)下降趨勢(shì)，東部地區(qū)波動(dòng)較大，且在2017年有翹尾，但總體上處于下降態(tài)勢(shì)。綜合來(lái)看，全國(guó)人均水環(huán)境壓力指數(shù)在2015年之后處于降低狀態(tài)。

4 分地區(qū)水環(huán)境賬戶狀況

4.1 西部地區(qū)

西部地區(qū)的水污染生態(tài)足跡與水環(huán)境生態(tài)承載力的具體狀況如圖2所示：

圖2 西部地區(qū)水污染生態(tài)足跡與水環(huán)境生態(tài)承載力Fig.2 Ecological footprint of water pollution and ecological carrying capacity of water environment in Western China

2013-2015年西部地區(qū)的人均總污染足跡遠(yuǎn)大于人均承載力，處于水環(huán)境赤字狀態(tài)，但自從2015年以后，人均水污染生態(tài)足跡迅速俯沖至水環(huán)境承載力以下，使之轉(zhuǎn)為水環(huán)境盈余并維持這一狀態(tài)。在此期間，人均水環(huán)境承載力先下降后上升，波動(dòng)較為平穩(wěn)，而水污染生態(tài)足跡變動(dòng)較大，可見(jiàn)西部的水環(huán)境盈虧主要取決于總的水污染生態(tài)足跡情況。

4.2 中部地區(qū)

中部地區(qū)的水污染生態(tài)足跡與水環(huán)境生態(tài)承載力的具體狀況如圖3所示：

圖3 中部地區(qū)水污染生態(tài)足跡與水環(huán)境生態(tài)承載力Fig.3 Ecological footprint of water pollution and ecological carrying capacity of water environment in Central China

中部地區(qū)的水環(huán)境除了2016年以外，其余年份均處于赤字狀態(tài)，且水環(huán)境承載力有下降趨勢(shì)，2017年恢復(fù)水環(huán)境赤字狀態(tài)。2013-2015年水污染足跡最高達(dá)到人均承載力的3倍，2015年之后人均總污染足跡顯著下降，2016年達(dá)到水環(huán)境盈余，但2017年由于承載力有所下降，又恢復(fù)赤字狀態(tài)，但比之前的赤字狀態(tài)緩解許多。

4.3 東部地區(qū)

東部地區(qū)的水污染生態(tài)足跡與水環(huán)境生態(tài)承載力的具體狀況如圖4所示。東部地區(qū)的水污染生態(tài)足跡與水環(huán)境生態(tài)承載力往相反方向波動(dòng)，形態(tài)上近似對(duì)稱圖形，2013-2015年水污染足跡緩慢上升，與此同時(shí)人均總承載力緩慢下降，2015年赤字達(dá)到峰值；2016年人均總污染足跡與人均總承載力之間的落差迅速收窄，赤字達(dá)到最小，之后赤字又逐漸擴(kuò)大。2015年國(guó)家通過(guò)《水污染防治行動(dòng)計(jì)劃》，簡(jiǎn)稱“水十條”，強(qiáng)制規(guī)定各省切實(shí)加強(qiáng)水污染防治力度，逾期未達(dá)標(biāo)的，決定取消其加劇水污染物排放的建設(shè)項(xiàng)目審批資格，并依照有關(guān)規(guī)定采取撤銷園區(qū)資格等懲罰措施[3]。由全國(guó)各地區(qū)的人均水環(huán)境走勢(shì)圖可看出，2015年以前人均水污染足跡總體處于高位，但2015年以后都急劇下降，大大緩解水環(huán)境赤字狀況，可見(jiàn)“水十條”對(duì)全國(guó)的水污染防治效果顯著，但中、東部地區(qū)在2016年以后又有恢復(fù)赤字的狀況出現(xiàn)，這種現(xiàn)象表明建立水污染防治的長(zhǎng)效機(jī)制任重道遠(yuǎn)。

圖4 東部地區(qū)水污染生態(tài)足跡與水環(huán)境生態(tài)承載力Fig.4 Ecological footprint of water pollution and ecological carrying capacity of water environment in eastern China

5 綠色因子分析

為進(jìn)一步分析人均水環(huán)境壓力指數(shù)與各地區(qū)綠色發(fā)展因素之間的關(guān)系，本研究參考《中國(guó)省級(jí)綠色經(jīng)濟(jì)指標(biāo)體系》和《中國(guó)區(qū)域產(chǎn)業(yè)綠色發(fā)展指數(shù)評(píng)價(jià)》，并結(jié)合與水環(huán)境賬戶有關(guān)的原則，綜合羅列了影響地方水環(huán)境壓力指數(shù)的22項(xiàng)綠色指標(biāo)，并分別與水環(huán)境壓力指數(shù)進(jìn)行相關(guān)性分析，以2017年數(shù)據(jù)為例，具體結(jié)果如表1。

單從相關(guān)性來(lái)看，pearson系數(shù)在0.5以上的只有人均壽命、平均受教育年限[4]、人均地區(qū)生產(chǎn)總值、第三產(chǎn)業(yè)勞動(dòng)生產(chǎn)率和人均可支配收入5項(xiàng)因子，即在綠色因子中，與人均水環(huán)境生態(tài)壓力指數(shù)有相對(duì)較強(qiáng)的線性相關(guān)性。代表第一、二產(chǎn)業(yè)發(fā)展水平的第一、二產(chǎn)業(yè)勞動(dòng)生產(chǎn)率沒(méi)有通過(guò)顯著性檢驗(yàn)，創(chuàng)新因子(13～15)、環(huán)境污染因子(16～18)、環(huán)境治理因子(20～22)以及人均耕地面積均未通過(guò)檢驗(yàn)，說(shuō)明這些因子與人均水環(huán)境壓力指數(shù)線性無(wú)關(guān)，但不排除有非線性關(guān)系。

對(duì)以上通過(guò)相關(guān)性檢驗(yàn)的10個(gè)指標(biāo)(1～10)進(jìn)行標(biāo)準(zhǔn)化處理后(消除量綱效應(yīng))[5]進(jìn)行因子分析，KMO值為0.801，巴特利特檢驗(yàn)結(jié)果顯著，故變量適合做因子分析。提取4個(gè)主成分，累計(jì)方差貢獻(xiàn)率為87.094%，即所提取的4個(gè)主成分可以解釋總體信息的87.094%，旋轉(zhuǎn)后的因子載荷矩陣如表2所示。

表1 人均水環(huán)境壓力指數(shù)與綠色發(fā)展因子相關(guān)性分析結(jié)果Tab.1 Analysis results of correlation between per capita water environmental pressure index and green development factors

表2 旋轉(zhuǎn)后的因子載荷矩陣Tab.2 Factor load matrix after rotation

根據(jù)各主成分在不同變量上的載荷大小進(jìn)行分組，并根據(jù)各因子結(jié)構(gòu)特征重新命名，如表3所示：

表3 因子分組結(jié)構(gòu)Tab.3 Factor grouping structure

6神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)分析

鑒于人工神經(jīng)網(wǎng)絡(luò)處理非線性問(wèn)題的優(yōu)越性，現(xiàn)運(yùn)用SPSS中神經(jīng)網(wǎng)絡(luò)的徑向基函數(shù)功能，選取2017年30個(gè)省份的數(shù)據(jù)，將地區(qū)發(fā)展水平、資源承載度、資源利用效率、政府支持度的標(biāo)準(zhǔn)化得分作為協(xié)變量，人均水環(huán)境壓力指數(shù)為因變量，直觀呈現(xiàn)四大綠色影響因素與水環(huán)境壓力指數(shù)之間存在的復(fù)雜非線性關(guān)系，并進(jìn)行因子相對(duì)重要性分析，結(jié)果如下：

圖5 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)值散點(diǎn)圖Fig.5 Scatter diagram of neural network prediction value

圖6 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)值殘差圖Fig.6 Residual map of neural network prediction value

神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)值的散點(diǎn)圖如圖5所示。從散點(diǎn)圖可以看出人均壓力指數(shù)的預(yù)測(cè)值基本呈現(xiàn)線性趨勢(shì)。神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)值的殘差圖如圖6所示。模型的殘差基本在0附近波動(dòng)，說(shuō)明模型效果好。自變量重要性如表4所示。未正態(tài)化前的變量重要性依次是地區(qū)發(fā)展水平0.494、資源承載度0.178、政府支持度0.167、資源利用效率0.161，正態(tài)化后的變量重要性依次是地區(qū)發(fā)展水平100%、資源承載度36%、政府支持度33.9%、資源利用效率32.6%。圖7是神經(jīng)網(wǎng)絡(luò)自變量正態(tài)化后的變量重要性條形圖?？梢?jiàn)地區(qū)發(fā)展水平在評(píng)價(jià)人均水環(huán)境壓力指數(shù)的過(guò)程中重要性占比最高，資源承載度、政府支持度、資源利用效率重要性均為30%左右，約占地區(qū)發(fā)展水平的1/3。這種分布特征為后續(xù)制定具體、有傾向化的政策方針提供了方向和依據(jù)。

表4 自變量的重要性Tab.4 Importance of independent variable

圖7 變量相對(duì)重要性條形圖Fig.7 Bar chart of relative importance of variables

7 結(jié)論

本研究通過(guò)計(jì)算2013-2017年全國(guó)30個(gè)省份的水污染生態(tài)足跡，表明大部分地區(qū)的水環(huán)境處于不安全的狀態(tài)，這一狀況在2015年以后有明顯好轉(zhuǎn)，但中東部地區(qū)有恢復(fù)赤字趨勢(shì)。因子分析綜合相關(guān)變量信息，并將輸出作為預(yù)測(cè)水環(huán)境壓力指數(shù)的參考變量，神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)的結(jié)果表明，要從根本上改善全國(guó)的水環(huán)境安全狀況，需從以下方面著手：

地區(qū)發(fā)展水平。目前地區(qū)的高度發(fā)展仍然是以水環(huán)境的惡化為代價(jià)的，第三產(chǎn)業(yè)勞動(dòng)生產(chǎn)率與人均水環(huán)境生態(tài)壓力指數(shù)成正相關(guān)，說(shuō)明僅調(diào)整產(chǎn)業(yè)結(jié)構(gòu)向第三產(chǎn)業(yè)轉(zhuǎn)移并不能有效緩解地區(qū)水環(huán)境壓力，而應(yīng)將綠色產(chǎn)業(yè)作為第三產(chǎn)業(yè)中最重要的經(jīng)濟(jì)增長(zhǎng)點(diǎn)，從根本上解決經(jīng)濟(jì)效益與水環(huán)境保護(hù)的矛盾，提升全國(guó)水環(huán)境安全水平。人均壽命、平均受教育年限與水環(huán)境壓力指數(shù)正相關(guān)，說(shuō)明人口素質(zhì)中的水環(huán)境保護(hù)意識(shí)有待加強(qiáng)，為培養(yǎng)居民的綠色用水習(xí)慣，應(yīng)將與水環(huán)境保護(hù)有關(guān)的知識(shí)培訓(xùn)納入國(guó)民教育體系，樹(shù)立防治水污染與節(jié)約水資源的綠色消費(fèi)理念。加速淘汰落后產(chǎn)能，提高項(xiàng)目環(huán)評(píng)的審核標(biāo)準(zhǔn)，嚴(yán)格規(guī)定用水及排放必須達(dá)到行業(yè)先進(jìn)水準(zhǔn)，方可授予開(kāi)發(fā)經(jīng)營(yíng)權(quán)益。整合政府部門(mén)與銀行業(yè)的信息資源，全面建立企業(yè)環(huán)境信用評(píng)價(jià)制度，并根據(jù)企業(yè)環(huán)評(píng)狀況實(shí)現(xiàn)貸款發(fā)放階梯化，放寬有利于水生態(tài)保護(hù)的企業(yè)貸款限制，嚴(yán)格限制環(huán)境違法企業(yè)貸款，營(yíng)造綠色信貸環(huán)境。另一方面，應(yīng)大力發(fā)展環(huán)保產(chǎn)業(yè)，盡快實(shí)現(xiàn)水環(huán)境保護(hù)體系化、服務(wù)化、產(chǎn)業(yè)化，打造一批現(xiàn)代化的環(huán)保服務(wù)型標(biāo)桿企業(yè)。

資源承載度。在人均水資源量一定的情況下，節(jié)約水資源對(duì)水環(huán)境安全起到至關(guān)重要的作用。將用水計(jì)劃管理落到實(shí)處，規(guī)定用水設(shè)施必須通過(guò)節(jié)水標(biāo)準(zhǔn)，否則不予投入使用。引進(jìn)新型適用的節(jié)水裝備和器具，全面更換市場(chǎng)中不符合節(jié)水標(biāo)準(zhǔn)的用水設(shè)備，并禁止其再次生產(chǎn)、出售。以政策補(bǔ)貼的方式激勵(lì)用戶更換落后的耗水型用具，選購(gòu)使用節(jié)水型水嘴、便器水箱等環(huán)保產(chǎn)品。充分發(fā)揮森林資源在涵養(yǎng)水源方面所起的作用，加大退耕還林、還草、還濕力度，不定期更新危害生態(tài)修復(fù)工作的商業(yè)開(kāi)發(fā)負(fù)面清單，扶持花卉苗木、森林旅游休閑等新型業(yè)態(tài)的適度規(guī)模經(jīng)營(yíng)，依托綠色資源，因地制宜。

政府支持度。加大政府資金投入，設(shè)立水環(huán)境保護(hù)專項(xiàng)基金，為實(shí)現(xiàn)項(xiàng)目開(kāi)發(fā)經(jīng)營(yíng)各環(huán)節(jié)的水環(huán)境監(jiān)控提供資金支持，將事故應(yīng)對(duì)轉(zhuǎn)為事前預(yù)防。參考以往環(huán)境違法典型案例，加緊推進(jìn)水環(huán)境保護(hù)立法工作，確定環(huán)境違法訴訟制度并建立咨詢委員會(huì)，追溯監(jiān)管部門(mén)、排污企業(yè)在水環(huán)境違法中的責(zé)任歸屬。企業(yè)為支持水環(huán)境保護(hù)政策所依賴進(jìn)口的核心原材料和零部件，應(yīng)予以減免關(guān)稅，將高耗水、高水污染產(chǎn)品納入消費(fèi)稅征收范圍。影響周邊居民水環(huán)境的大型建設(shè)項(xiàng)目需舉行公開(kāi)聽(tīng)證會(huì)，限期研究并解決群眾反映強(qiáng)烈的問(wèn)題。

資源利用效率。為有效降低單位地區(qū)生產(chǎn)總值水耗，減少工業(yè)固體廢棄物中有害物質(zhì)所造成的水體污染，應(yīng)加大包括水循環(huán)在內(nèi)的資源綜合利用領(lǐng)域的經(jīng)濟(jì)投入，加緊相關(guān)重大工藝及材料的技術(shù)突破。對(duì)自主研發(fā)循環(huán)利用產(chǎn)品并探索市場(chǎng)化路徑的企業(yè)實(shí)行降低準(zhǔn)入，稅收優(yōu)惠和政策補(bǔ)貼，充分調(diào)動(dòng)資本市場(chǎng)的投資熱情。