張義，程戀婷

(河池學(xué)院 經(jīng)濟(jì)與管理學(xué)院，廣西 宜州 546300)

基于改進(jìn)模型的廣西北部灣經(jīng)濟(jì)區(qū)水污染生態(tài)足跡時(shí)空分析

張義，程戀婷

(河池學(xué)院 經(jīng)濟(jì)與管理學(xué)院，廣西 宜州 546300)

當(dāng)前的水污染生態(tài)足跡模型存在較大不足：明顯低估水污染程度；水資源均衡因子值可靠性不高；水資源產(chǎn)量因子多余；無(wú)法區(qū)分水環(huán)境承載力與水資源承載力。為此，我們對(duì)原模型進(jìn)行了如下改進(jìn)：從點(diǎn)源和面源角度核算排放到水體的多種污染物；棄用均衡因子和產(chǎn)量因子；單獨(dú)核算水環(huán)境承載力。接下來(lái)以2008-2013年的廣西北部灣經(jīng)濟(jì)區(qū)為例進(jìn)行了計(jì)算與分析。計(jì)算所需的數(shù)據(jù)來(lái)自相應(yīng)年份的廣西水資源公報(bào)、廣西統(tǒng)計(jì)年鑒、廣西環(huán)保廳內(nèi)部資料。核算結(jié)果表明：(1)經(jīng)濟(jì)區(qū)及其各市均為氮污染水生態(tài)足跡最大，磷污染水生態(tài)足跡次之，有機(jī)物污染水生態(tài)足跡最小。該地區(qū)有機(jī)物污染足跡呈下降趨勢(shì)，2008-2010年最大驅(qū)動(dòng)因子為工業(yè)賬戶(hù)，2011-2013年則為城鎮(zhèn)生活賬戶(hù)；氮/磷污染足跡總體表現(xiàn)為上升，其最大影響因子來(lái)自畜禽養(yǎng)殖賬戶(hù)。各市3種水污染生態(tài)足跡大小排序不一致，其變化趨勢(shì)也各異。(2)經(jīng)濟(jì)區(qū)水環(huán)境承載力先降后升，總體略呈上升趨勢(shì)；內(nèi)部非常不平衡，其中，防城港市該值很大，其他3市則較小。(3)基于水環(huán)境生態(tài)盈虧和水環(huán)境生態(tài)壓力指數(shù)，經(jīng)濟(jì)區(qū)水環(huán)境安全度較差；內(nèi)部差異顯著，其中，防城港市很好，欽州市和北海市較差，南寧市很差。最后，提出了改善該地區(qū)水環(huán)境狀況的建議，同時(shí)指出了改進(jìn)模型的優(yōu)點(diǎn)及有待進(jìn)一步研究的幾個(gè)問(wèn)題。

水污染生態(tài)足跡；水環(huán)境承載力；水環(huán)境生態(tài)盈虧；水環(huán)境生態(tài)壓力指數(shù)；廣西北部灣經(jīng)濟(jì)區(qū)

生態(tài)足跡是由生態(tài)經(jīng)濟(jì)學(xué)家Rees和Wackernagel于20世紀(jì)90年代初提出的衡量自然資源可持續(xù)利用的方法，近年來(lái)在全球范圍內(nèi)得到廣泛應(yīng)用。它的優(yōu)點(diǎn)主要是角度新穎、形象生動(dòng)、容易理解、可操作性強(qiáng)、研究結(jié)果可比性高，但也存在諸多不足[1]73-77。水資源生態(tài)足跡是研究者針對(duì)傳統(tǒng)生態(tài)足跡忽略水域水資源供給和水質(zhì)凈化功能而提出來(lái)的。它通過(guò)建立水資源用地這一新的土地類(lèi)型將人類(lèi)對(duì)上述兩種生態(tài)系統(tǒng)服務(wù)的占用納入生態(tài)足跡核算框架，一般包括淡水生態(tài)足跡和水污染生態(tài)足跡。其中，淡水生態(tài)足跡反映人類(lèi)在生產(chǎn)生活過(guò)程中對(duì)淡水資源的消費(fèi)，該方面的研究成果眾多且其核算方法高度一致[2-4]；水污染生態(tài)足跡反映的則是人類(lèi)活動(dòng)排放的廢棄物對(duì)水環(huán)境凈化功能的占用程度，該方面的研究成果較少，測(cè)算方法各異且存在某些不足之處。洪輝[5]構(gòu)建的水污染生態(tài)足跡模型(使用者相對(duì)較多)的主要問(wèn)題有：僅核算點(diǎn)源污染而未考慮面源污染，并且僅測(cè)算未達(dá)標(biāo)廢污水占用的水資源用地，故明顯低估水污染生態(tài)足跡；水資源均衡因子值計(jì)算方法迥異于其他均衡因子值估值方法；水資源產(chǎn)量因子多余；無(wú)法區(qū)分水環(huán)境承載力與水資源承載力[6]。段錦等[7]、于冰等[8]提出的水污染生態(tài)足跡計(jì)算方法是基于水域?qū)ξ鬯奈{能力，而實(shí)際上應(yīng)核算的是污水中的污染物對(duì)水環(huán)境凈化功能的占用狀況，同時(shí)不同來(lái)源污水中的污染物濃度也是不相同的，因而這種核算精確度較差。白鈺等[9]將水污染物占用生態(tài)足跡的方式分為“污染物處理”和“污染物危害”兩類(lèi)，重復(fù)計(jì)算了水污染對(duì)水生態(tài)系統(tǒng)的占用。另外，以貨幣為媒介將水污染生態(tài)足跡轉(zhuǎn)化為水產(chǎn)品生態(tài)足跡也存在是否適宜的問(wèn)題。閔慶文等[10]明確提出并界定了“基于生態(tài)系統(tǒng)服務(wù)的生態(tài)足跡”概念，并對(duì)其中的污染足跡的基本理論與模型進(jìn)行了闡釋與構(gòu)建，為水污染生態(tài)足跡的研究提供了理論基礎(chǔ)。焦雯珺等[11]利用該污染足跡模型的簡(jiǎn)化版開(kāi)展了案例研究，但從中可以發(fā)現(xiàn)該模型中的某些參數(shù)值(如污染物吸納服務(wù)供給因子)很難得到，因而推廣不易。因此水污染生態(tài)足跡研究有待進(jìn)一步深入。本文在洪輝構(gòu)建的模型的基礎(chǔ)上提出了改進(jìn)的水污染生態(tài)足跡模型，并以之核算了廣西北部灣經(jīng)濟(jì)區(qū)2008-2013年的水污染生態(tài)足跡、水環(huán)境承載力、水環(huán)境生態(tài)赤字及水環(huán)境生態(tài)壓力指數(shù)，同時(shí)分析其時(shí)空特征，以期有益于水生態(tài)足跡研究并為該地區(qū)水資源(環(huán)境)科學(xué)管理與可持續(xù)利用提供參考依據(jù)。

1 研究區(qū)概況

廣西北部灣經(jīng)濟(jì)區(qū)位于廣西南部，主要包括南寧市、北海市、欽州市和防城港市4個(gè)地級(jí)市。區(qū)域范圍東經(jīng)107°20′-109°47′，北緯20°26′-24°03′，東連廣東省及玉林市，南瀕臨北部灣，與海南省隔海相望，西南與越南和崇左市毗鄰，北與貴港市、來(lái)賓市、河池市、百色市等相連。它是我國(guó)西部大開(kāi)發(fā)地區(qū)唯一的沿海區(qū)域，區(qū)位優(yōu)勢(shì)明顯，戰(zhàn)略地位突出。經(jīng)濟(jì)區(qū)總面積4.22萬(wàn)km2，約占廣西總面積的18%；2013年末總?cè)丝? 383.37萬(wàn)，約占廣西同期總?cè)丝诘?6%；2013年末地區(qū)生產(chǎn)總值為4 817.43億元，約占廣西同期地區(qū)生產(chǎn)總值的34%。經(jīng)濟(jì)區(qū)多年平均水資源量為349.6億m3，約占廣西水資源總量的18%。境內(nèi)的水系為珠江流域西江水系和桂南沿海諸河水系，前者的主要河流是郁江、右江、左江、紅水河、清水江，后者主要包括南流江、欽江、大風(fēng)江、茅嶺江、防城河。經(jīng)濟(jì)區(qū)自2008年上升為國(guó)家戰(zhàn)略以來(lái)經(jīng)濟(jì)快速發(fā)展，但與此同時(shí)其水環(huán)境承受的壓力亦顯著增加。

2 研究方法與參數(shù)取值

2.1 改進(jìn)的水污染生態(tài)足跡模型

本模型包括水污染生態(tài)足跡、水環(huán)境承載力、水環(huán)境生態(tài)盈虧和水環(huán)境生態(tài)壓力指數(shù)4個(gè)核算項(xiàng)目。水污染生態(tài)足跡指吸納(稀釋和凈化)一定人口產(chǎn)生的水體污染物所需要的水資源用地面積[6]。它反映的是人類(lèi)活動(dòng)對(duì)水生態(tài)系統(tǒng)水質(zhì)凈化功能的占用情況，其計(jì)算思路如下：先核算將某區(qū)域人類(lèi)排放到水體的污染物稀釋到Ⅲ類(lèi)水質(zhì)標(biāo)準(zhǔn)所需消耗的水資源量(如此處理的理由是Ⅲ類(lèi)水基本能夠維持水體的生態(tài)服務(wù)功能，而低于此標(biāo)準(zhǔn)的水體的生態(tài)服務(wù)功能會(huì)衰退或喪失[12])，然后以之除以水資源全球平均生產(chǎn)能力就可以得到該地區(qū)水污染生態(tài)足跡。本文選取具有代表性的氮、磷和有機(jī)物作為研究?jī)?nèi)容，以其中的最大污染足跡作為最終水污染生態(tài)足跡(因?yàn)檫@三種污染物在環(huán)境影響上具有明顯重疊效應(yīng)[13])?；谏鷳B(tài)足跡的水環(huán)境承載力是指某區(qū)域水體能夠被繼續(xù)使用并仍保持良好生態(tài)系統(tǒng)時(shí)所能提供的用于容納污染物的水資源用地面積，它反映的是水生態(tài)系統(tǒng)水質(zhì)凈化功能的供給能力。兩者之差為水環(huán)境生態(tài)盈虧，用來(lái)反映人類(lèi)對(duì)水環(huán)境凈化功能的占用是否可持續(xù)，盈余為可持續(xù)，赤字為不可持續(xù)。兩者之比為水環(huán)境生態(tài)壓力指數(shù)，可用來(lái)反映研究區(qū)域水環(huán)境承受的壓力狀況或水環(huán)境安全程度，具體劃分標(biāo)準(zhǔn)見(jiàn)表1。

表1 水環(huán)境生態(tài)壓力指數(shù)的等級(jí)劃分標(biāo)準(zhǔn)

注：見(jiàn)參考文獻(xiàn)[14]

改進(jìn)模型可表示為：

(1)

ECWE=0.88×(QWP/PW)

(2)

QWP=WQ-(QR×K)

(3)

EDWE(ESWE)=ECWE-EFWP

(4)

EPIWE=EFWP/ECWE

(5)

式中，EFWP、EFNP、EFPP和EFCODP分別為水污染生態(tài)足跡、氮污染水生態(tài)足跡、磷污染水生態(tài)足跡和有機(jī)物污染水生態(tài)足跡(hm2)；CCOD、CN和CP為分別為研究區(qū)排入水體的有機(jī)物、總氮和總磷量(g)；UN、UP和UCOD分別為Ⅲ類(lèi)水質(zhì)標(biāo)準(zhǔn)中的總氮、總磷和化學(xué)需氧量含量上限(g/m3)；PW為水資源世界平均生產(chǎn)能力(m3/hm2)；ECWE為研究區(qū)水環(huán)境承載力(hm2)；0.88來(lái)自傳統(tǒng)生態(tài)足跡模型，指的是承載力中保留12%用于生物多樣性保護(hù)；QW、QWP和QR分別為研究區(qū)水資源總量、用于接納污染物的水量和取水量(m3)；K為研究區(qū)綜合耗水率；EDWE/ESWE為研究區(qū)水環(huán)境生態(tài)赤字或盈余(hm2)，其中，大于0為盈余狀態(tài)，小于0為赤字狀態(tài)，等于0為臨界狀態(tài)；EPIWE為研究區(qū)水環(huán)境生態(tài)壓力指數(shù)。

本模型舍棄了水資源均衡因子，這是因?yàn)椋壕庖蜃雍茈y獲得其合理取值[15]；而且水資源均衡因子值來(lái)自?xún)r(jià)值量評(píng)價(jià)法，而生態(tài)足跡模型中的其他均衡因子值來(lái)自物質(zhì)量評(píng)價(jià)法，不同方法求出的參數(shù)使其模型的有效性受到明顯影響[6]。

因?yàn)檠芯繀^(qū)水資源總量除以水資源全球平均生產(chǎn)能力得到的即是全球意義上的水資源用地面積，故棄用傳統(tǒng)模型中的水資源產(chǎn)量因子。如果保留它，其后果是：一方面低估了水資源貧乏地區(qū)的水環(huán)境承載力(因水資源產(chǎn)量因子小于1)，另一方面則高估了水資源豐富地區(qū)的水環(huán)境承載力(因水資源產(chǎn)量因子大于1)，因而不能比較準(zhǔn)確地反映研究區(qū)水環(huán)境承載力狀況。

2.2 數(shù)據(jù)來(lái)源和參數(shù)取值

研究區(qū)水資源總量、取水量和綜合耗水率來(lái)自《廣西壯族自治區(qū)水資源公報(bào)》。人口數(shù)據(jù)來(lái)自研究區(qū)統(tǒng)計(jì)年鑒。水資源世界平均生產(chǎn)能力取3 140 m3/hm2[2]。Ⅲ類(lèi)水質(zhì)標(biāo)準(zhǔn)中的化學(xué)需氧量、總氮和總磷含量上限分別取20 mg/L、1.0 mg/L和0.2 mg/L[17]。

受數(shù)據(jù)限制，本研究主要估算研究區(qū)城鎮(zhèn)及農(nóng)村生活污水、工業(yè)廢水、耕地地面徑流和畜禽養(yǎng)殖向水生態(tài)系統(tǒng)排放的化學(xué)需氧量、總氮和總磷量。其中，工業(yè)和城鎮(zhèn)生活化學(xué)需氧量入水量來(lái)自廣西環(huán)保廳內(nèi)部資料，城鎮(zhèn)生活氮/磷入水量及其他污染源產(chǎn)生的污染物的入水量通過(guò)以下方法估算(由于數(shù)據(jù)難以獲得，沒(méi)有估算工業(yè)廢水中的總氮、總磷量)。

研究區(qū)城鎮(zhèn)生活污水中的總氮/總磷排放量估算方法如下：

CP1=P1-E1=(3 650×RT×F1)-TW×(D1-D0)/1 000

(6)

D1=P1/EW×1 000

(7)

式中，CP1為城鎮(zhèn)生活污水總氮/總磷排放量(kg/a)；P1為城鎮(zhèn)生活污水總氮/總磷產(chǎn)生量(kg/a)；E1為污水處理廠(chǎng)去掉的總氮/總磷量(kg/a)；RT為城鎮(zhèn)居民常住人口(104人)；F1為城鎮(zhèn)生活污水總氮/總磷人均產(chǎn)生系數(shù)，具體取值來(lái)自文獻(xiàn)[18]；TW為污水處理廠(chǎng)處理的污水量(m3)，可從研究區(qū)統(tǒng)計(jì)年鑒獲得；DI為污水進(jìn)廠(chǎng)時(shí)總氮/總磷濃度(mg/L)；DO為污水出廠(chǎng)時(shí)總氮/總磷濃度(mg/L)，根據(jù)污水處理廠(chǎng)出水一級(jí)標(biāo)準(zhǔn)的B標(biāo)準(zhǔn)取值[19]，其中總氮濃度取20 mg/L，總磷濃度取1.5 mg/L；EW為城鎮(zhèn)居民生活污水排放量(m3)，可從研究區(qū)統(tǒng)計(jì)年鑒獲得。

研究區(qū)農(nóng)村生活污水中的化學(xué)需氧量、總氮、總磷排放量估算方法如下：

CP2=P2×β1=(3 650×RC×F2)×β1

(8)

式中，CP2為農(nóng)村生活污水污染物(化學(xué)需氧量/總氮/總磷)排放量(kg/a)；P2為農(nóng)村生活污水污染物(化學(xué)需氧量、總氮、總磷)產(chǎn)生量(kg/a)；β1為農(nóng)村生活污染物(化學(xué)需氧量、總氮、總磷)入河系數(shù)，一般取值為0.2-0.5[20]30，本文取0.35；RC為農(nóng)村人口數(shù)(104人)；F2為農(nóng)村生活污水污染物人均產(chǎn)生系數(shù)(g/人·d)，本文分別取16.4 g/人·d(化學(xué)需氧量)、5 g/人·d(總氮)、0.44 g/人·d(總磷)[20]31。

耕地地面徑流中的化學(xué)需氧量、總氮、總磷入河量可根據(jù)式(9)計(jì)算得到：

CP3=(M1×F3)×β2×γ1

(9)

式中，CP3為耕地污染物(化學(xué)需氧量/總氮/總磷)入河量(t/a)；M1為耕地面積(km2)，可從研究區(qū)統(tǒng)計(jì)年鑒獲得；β2為耕地污染物入河系數(shù)，取0.2(文獻(xiàn)[20]中取值范圍的中間值)；γ1為修正系數(shù)[20]31，當(dāng)耕地化肥施用量在25 kg以下時(shí)取0.8-1.0(本文取0.9)，當(dāng)施用量在25-35 kg時(shí)取1.0-1.2(本文取1.1)，當(dāng)施用量在35 kg以上時(shí)取1.2-1.5(本文取1.3)，耕地化肥使用量可從研究區(qū)統(tǒng)計(jì)年鑒獲得；F3為耕地排污系數(shù)，分別取0.067 t/km2(化學(xué)需氧量)[20]31、2.2 t/km2(總氮)[6]、0.16 t/km2(總磷)[6]。

畜禽養(yǎng)殖(本文主要考慮豬、牛、羊和家禽養(yǎng)殖)向水體排放的化學(xué)需氧量、總氮、總磷的估算方法如下[21]：

CP4=(N×F4)×C×β3

(10)

式中，CP4為畜禽養(yǎng)殖向水體排放的化學(xué)需氧量/總氮/總磷(kg/a)；N為畜禽飼養(yǎng)量(萬(wàn)頭、只)，可從研究區(qū)統(tǒng)計(jì)年鑒獲得；F4為畜禽糞便排泄系數(shù)(kg/a);C為單位質(zhì)量畜禽糞便中化學(xué)需氧量/總氮/總磷含量(%);β3為畜禽糞便中化學(xué)需氧量/總氮/總磷入水率(%)，各參數(shù)取值見(jiàn)表2。

表2 畜禽糞便排泄系數(shù)、養(yǎng)分含量和入水率

注：數(shù)據(jù)來(lái)源于文獻(xiàn)[22]第26，77-78頁(yè)。

3 結(jié)果與分析

3.1 廣西北部灣經(jīng)濟(jì)區(qū)水污染生態(tài)足跡時(shí)空分析

3.1.1 廣西北部灣經(jīng)濟(jì)區(qū)水污染生態(tài)足跡動(dòng)態(tài)分析

2008-2013年，經(jīng)濟(jì)區(qū)人均有機(jī)物污染水生態(tài)足跡介于0.307 7-0.459 4 hm2之間，期間呈明顯下降趨勢(shì)(年均降幅為6.60%)(圖1)。這反映出近年來(lái)經(jīng)濟(jì)區(qū)有機(jī)物污染對(duì)水體的壓力大幅下降。進(jìn)一步分析其二級(jí)賬戶(hù)可以發(fā)現(xiàn)：工業(yè)賬戶(hù)值迅速下降(年均降幅高達(dá)13.42%)，城鎮(zhèn)生活賬戶(hù)值也有明顯下降(年均降幅為4.13%)，農(nóng)村生活賬戶(hù)、畜禽賬戶(hù)和耕地賬戶(hù)值變化很小，表明經(jīng)濟(jì)區(qū)有機(jī)物污染水生態(tài)足跡的跌落主要來(lái)自工業(yè)有機(jī)物和城鎮(zhèn)生活有機(jī)物的減排。從均值來(lái)看，城鎮(zhèn)生活賬戶(hù)年均值(0.141 0 hm2)及所占比例(36%左右)最大，緊隨其后的是工業(yè)賬戶(hù)(年均值和比例分別為0.125 2 hm2和32%左右)，畜禽養(yǎng)殖賬戶(hù)次之(年均值和比例分別為0.102 8 hm2和26%左右)，農(nóng)村生活賬戶(hù)(年均值和比例分別為0.018 8 hm2和5%左右)，耕地賬戶(hù)微不足道(年均值和比例分別為0.000 2 hm2和0.06%左右)，表明該足跡的最大驅(qū)動(dòng)因子來(lái)自城鎮(zhèn)生活，其次是工業(yè)和畜禽養(yǎng)殖。但從走向來(lái)看，工業(yè)賬戶(hù)貢獻(xiàn)率呈明顯下降趨勢(shì)，說(shuō)明其驅(qū)動(dòng)力在減弱，尤其是最近3年所占比例由首位降至第3位；城鎮(zhèn)生活賬戶(hù)貢獻(xiàn)率呈上升態(tài)勢(shì)，說(shuō)明其驅(qū)動(dòng)力逐漸增強(qiáng)；畜禽養(yǎng)殖賬戶(hù)比重上升較快，最近3年所占比例已超工業(yè)賬戶(hù)；其他賬戶(hù)變化較小。

圖1 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)人均有機(jī)物污染水生態(tài)足跡

經(jīng)濟(jì)區(qū)人均氮污染水生態(tài)足跡介于1.294 6-1.438 6 hm2之間，其走向?yàn)橄壬蠼?，總體略呈上升態(tài)勢(shì)(年均增幅為1.13%)(見(jiàn)圖2)。這說(shuō)明近年來(lái)經(jīng)濟(jì)區(qū)氮污染對(duì)水體的壓力有所上升。其中，城鎮(zhèn)生活賬戶(hù)和畜禽養(yǎng)殖賬戶(hù)值波動(dòng)中呈上升走向(年均增幅分別為2.36%和1.22%)，農(nóng)村生活賬戶(hù)和耕地賬戶(hù)值變化較小，說(shuō)明經(jīng)濟(jì)區(qū)氮污染水生態(tài)足跡的上升主要來(lái)自城鎮(zhèn)生活和畜禽養(yǎng)殖氮排放的增加。從均值來(lái)看，畜禽養(yǎng)殖賬戶(hù)是主體(年均值和比例分別為0.803 1 hm2和59%左右)，其次是城鎮(zhèn)生活賬戶(hù)(年均值和比例分別為0.261 8 hm2和19%左右)，農(nóng)村生活賬戶(hù)和耕地賬戶(hù)相對(duì)較小(年均值分別為0.151 0 hm2和0.141 8 hm2，所占比例分別為11%和10%左右)，表明該足跡的最大驅(qū)動(dòng)因子來(lái)自畜禽養(yǎng)殖。從時(shí)間序列來(lái)看，各賬戶(hù)貢獻(xiàn)率變化相對(duì)較小，說(shuō)明該足跡的驅(qū)動(dòng)機(jī)制沒(méi)有明顯改變。

圖2 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)人均氮污染水生態(tài)足跡

經(jīng)濟(jì)區(qū)人均磷污染水生態(tài)足跡處于0.935 8-1.045 7 hm2之間，期間表現(xiàn)為先升后降，總體呈增長(zhǎng)趨勢(shì)(年均增幅為1.69%，見(jiàn)圖3)。經(jīng)分析，經(jīng)濟(jì)區(qū)磷污染水生態(tài)足跡的上升也主要來(lái)自城鎮(zhèn)生活賬戶(hù)和畜禽養(yǎng)殖賬戶(hù)。從均值來(lái)看，畜禽養(yǎng)殖賬戶(hù)占絕對(duì)優(yōu)勢(shì)(年均值和比例分別為0.792 2 hm2和79%左右)，城鎮(zhèn)生活賬戶(hù)年均值和比例分別為0.102 5 hm2和10%左右)，農(nóng)村生活賬戶(hù)和耕地賬戶(hù)年均值(分別為0.051 9 hm2和0.051 6 hm2)和所占比例(均為5%左右)較小，表明該足跡的主要驅(qū)動(dòng)因子源自畜禽養(yǎng)殖。各賬戶(hù)貢獻(xiàn)率變化相對(duì)較小，說(shuō)明該足跡的驅(qū)動(dòng)格局也沒(méi)有發(fā)生多大變化。

圖3 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)人均磷污染水生態(tài)足跡

3種水污染生態(tài)足跡中，氮污染水生態(tài)足跡最大，磷污染水生態(tài)足跡次之，有機(jī)物污染水生態(tài)足跡最小，故以氮污染水生態(tài)足跡作為最終水污染生態(tài)足跡。

3.1.2 廣西北部灣經(jīng)濟(jì)區(qū)水污染生態(tài)足跡空間分析

2008-2013年，經(jīng)濟(jì)區(qū)內(nèi)各市人均有機(jī)物污染水生態(tài)足跡介于0.239 7-0.683 9 hm2之間(見(jiàn)圖4)，年度最大最小值比介于1.46-2.19之間，說(shuō)明其空間差異度不大?；谀昃?見(jiàn)圖4)，各市大小排序?yàn)榉莱歉?北海>南寧>欽州。從走向來(lái)看(見(jiàn)圖4)，各市皆呈下降趨勢(shì)(其中，欽州市和北海市年均降幅9%左右，南寧市和防城港市年均降幅約5%)，這說(shuō)明各市有機(jī)污染物排放對(duì)當(dāng)?shù)厮w的壓力日益減輕，主因是其工業(yè)和城鎮(zhèn)生活化學(xué)需氧量排放量降幅非常明顯。進(jìn)一步分析其二級(jí)賬戶(hù)可以發(fā)現(xiàn)(見(jiàn)表3)：防城港市和北海市最大驅(qū)動(dòng)力分別來(lái)自工業(yè)賬戶(hù)和城鎮(zhèn)生活賬戶(hù)(其年均貢獻(xiàn)率分別約為52%和44%)；南寧市和欽州市2008-2010年最大驅(qū)動(dòng)力源自工業(yè)賬戶(hù)(其貢獻(xiàn)率介于36-43%之間)，2011-2013年則為城鎮(zhèn)生活賬戶(hù)(其貢獻(xiàn)率介于38-43%之間)。第二位驅(qū)動(dòng)因子方面，防城港市是城鎮(zhèn)生活賬戶(hù)(其年均貢獻(xiàn)率約為30%左右)；北海市2008-2010年為工業(yè)賬戶(hù)，2011-2013年則為畜禽養(yǎng)殖和工業(yè)賬戶(hù)，其年均貢獻(xiàn)率分別約為27%和29%；南寧市和欽州市2008-2010年為城鎮(zhèn)生活賬戶(hù)(其年均貢獻(xiàn)率分別約為30%和34%)，2011-2013年則為畜禽養(yǎng)殖賬戶(hù)(其年均貢獻(xiàn)率分別約為35%和36%)。從時(shí)間序列來(lái)看(見(jiàn)表3)，南寧市和欽州市的城鎮(zhèn)生活賬戶(hù)和畜禽養(yǎng)殖賬戶(hù)的驅(qū)動(dòng)力正明顯增強(qiáng)而工業(yè)賬戶(hù)的驅(qū)動(dòng)力日漸減弱；北海市城鎮(zhèn)生活賬戶(hù)驅(qū)動(dòng)力弱化明顯而其畜禽養(yǎng)殖賬戶(hù)驅(qū)動(dòng)力日益增強(qiáng)；防城港市的各類(lèi)賬戶(hù)貢獻(xiàn)率變化相對(duì)較小，說(shuō)明其該足跡的驅(qū)動(dòng)機(jī)制沒(méi)有明顯改變。

圖4 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)各市人均有機(jī)物污染水生態(tài)足跡

年份各賬戶(hù)占南寧有機(jī)物污染水生態(tài)足跡比例各賬戶(hù)占北海有機(jī)物污染水生態(tài)足跡比例城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖工業(yè)城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖工業(yè)200829.534.130.0625.1341.1551.203.270.0317.2228.28200930.554.180.0726.0739.1353.613.370.0418.2024.79201031.364.070.0628.5136.0048.833.690.0419.0128.43201142.774.620.0832.6319.9036.075.560.0629.6528.67201239.194.830.0835.6920.2136.405.240.0529.5228.79201342.285.000.0836.3516.2735.375.680.0631.2727.63均值35.954.470.0730.7328.7843.584.470.0424.1427.76年份各賬戶(hù)占防城港有機(jī)物污染水生態(tài)足跡比例各賬戶(hù)占?xì)J州有機(jī)物污染水生態(tài)足跡比例城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖工業(yè)城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖工業(yè)200828.943.340.0411.7855.9033.564.940.0318.7642.71200930.793.300.0511.9953.8733.634.890.0219.3042.17201030.722.810.0411.6254.8034.845.520.0320.7138.90201131.503.430.0515.0949.9238.189.390.0535.9916.37201231.423.410.0616.1448.9639.259.040.0536.7414.92201331.403.830.0616.9447.7740.279.320.0537.6712.67均值30.803.350.0513.9351.8736.627.180.0428.2027.96

2008-2013年，經(jīng)濟(jì)區(qū)內(nèi)各市人均氮污染水生態(tài)足跡介于0.905 9-1.738 8 hm2之間(見(jiàn)圖5)，年度最大最小值比介于1.55-1.77之間，說(shuō)明其空間差異度相對(duì)較小?；谀昃?見(jiàn)圖5)，各市大小排序?yàn)槟蠈?防城港>北海>欽州。從走向來(lái)看(見(jiàn)圖5)，南寧市和欽州市總體呈上升趨勢(shì)，說(shuō)明它們的氮排放對(duì)其水體的壓力在加大；北海市和防城港市大體表現(xiàn)為下降，表明這2個(gè)地區(qū)的氮排放對(duì)其水體的壓力在減輕。進(jìn)一步分析其二級(jí)賬戶(hù)可以發(fā)現(xiàn)(見(jiàn)表4)：畜禽養(yǎng)殖賬戶(hù)是各市氮污染水生態(tài)足跡的最大驅(qū)動(dòng)因子，前者對(duì)后者的年均貢獻(xiàn)率介于50-64%之間。南寧市、北海市和防城港市該足跡的第二位驅(qū)動(dòng)因子是城鎮(zhèn)生活賬戶(hù)(其年均貢獻(xiàn)率介于20-29%之間)；欽州市該足跡的第二位驅(qū)動(dòng)因子則為農(nóng)村生活賬戶(hù)(其年均貢獻(xiàn)率為16%左右)。從時(shí)間序列來(lái)看(見(jiàn)表4)，防城港市的城鎮(zhèn)生活賬戶(hù)的驅(qū)動(dòng)力在減弱同時(shí)其畜禽養(yǎng)殖賬戶(hù)的驅(qū)動(dòng)力在增強(qiáng)；其余地區(qū)的各類(lèi)賬戶(hù)貢獻(xiàn)率變化較小，說(shuō)明它們?cè)撟阚E的驅(qū)動(dòng)機(jī)制沒(méi)有明顯改變。

圖5 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)各市人均氮污染水生態(tài)足跡

年份各賬戶(hù)占南寧氮污染水生態(tài)足跡比例各賬戶(hù)占北海氮污染水生態(tài)足跡比例城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖200817.4711.7211.0759.7325.508.429.1156.98200917.6911.5112.1658.6522.068.729.6759.55201018.4610.1810.9660.3917.689.739.9562.63201126.039.3510.4154.2218.999.289.7262.01201218.979.9311.3159.7818.968.899.5262.63201322.589.6110.5257.2818.459.179.5862.80均值20.2010.3811.0758.3420.289.039.5961.10年份各賬戶(hù)占防城港氮污染水生態(tài)足跡比例各賬戶(hù)占?xì)J州氮污染水生態(tài)足跡比例城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖200836.247.9911.2244.5510.6216.408.2264.76200933.618.0312.6045.7610.0016.117.9365.97201030.237.6413.1149.0310.8416.727.5064.94201123.467.9814.2454.3116.2315.248.4360.10201224.577.4513.6854.2914.8014.848.3162.05201324.637.9713.3254.089.3915.718.6366.26均值28.797.8413.0350.3411.9815.848.1764.01

2008-2013年，經(jīng)濟(jì)區(qū)內(nèi)各市人均磷污染水生態(tài)足跡介于0.791 1-1.168 6 hm2之間(見(jiàn)圖6)，年度最大最小值比介于1.21-1.48之間，說(shuō)明其空間差異度較小?；谀昃?見(jiàn)圖6)，各市大小排序?yàn)槟蠈?北海>欽州>防城港。從走向來(lái)看(見(jiàn)圖6)，南寧市和欽州市總體呈上升趨勢(shì)，防城港市大體表現(xiàn)為下降，北海市則呈波動(dòng)態(tài)勢(shì)。各市磷污染水生態(tài)足跡的驅(qū)動(dòng)機(jī)制與其氮污染水生態(tài)足跡驅(qū)動(dòng)機(jī)制基本相同，不過(guò)前者的畜禽養(yǎng)殖賬戶(hù)貢獻(xiàn)率要明顯大于后者(前者年均貢獻(xiàn)率介于66-85%之間)(見(jiàn)表5)。

圖6 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)各市人均磷污染水生態(tài)足跡

年份各賬戶(hù)占南寧磷污染水生態(tài)足跡比例各賬戶(hù)占北海磷污染水生態(tài)足跡比例城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖20089.775.095.9579.1915.945.494.6973.8920099.644.906.4179.0612.875.524.8376.78201010.484.455.9379.149.815.944.7979.46201115.524.085.6374.7610.495.674.6979.15201210.254.145.8479.7710.375.404.5679.67201312.934.125.5877.3710.065.554.5879.82年均值11.434.465.8978.2211.595.594.6978.13年份各賬戶(hù)占防城港磷污染水生態(tài)足跡比例各賬戶(hù)占?xì)J州磷污染水生態(tài)足跡比例城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖城鎮(zhèn)生活農(nóng)村生活耕地畜禽養(yǎng)殖200826.855.917.6059.644.316.883.1585.66200924.485.858.4161.263.996.683.0086.33201021.235.408.4764.904.376.972.8685.80201115.675.438.8570.046.826.553.3183.33201216.265.048.4770.236.096.273.2184.44201316.245.388.2370.153.756.513.2786.46年均值20.125.508.3466.044.896.643.1385.34

各市3種水污染生態(tài)足跡中均為氮污染足跡最大，磷污染足跡次之，有機(jī)物污染足跡最小，故以其中的氮污染足跡作為最終水污染生態(tài)足跡。

3.2 廣西北部灣經(jīng)濟(jì)區(qū)水環(huán)境承載力、生態(tài)盈虧及生態(tài)壓力時(shí)空分析

3.2.1 廣西北部灣經(jīng)濟(jì)區(qū)水環(huán)境承載力、生態(tài)盈虧及生態(tài)壓力動(dòng)態(tài)分析

2008-2013年，經(jīng)濟(jì)區(qū)人均水環(huán)境承載力介于0.447 6-0.823 5 hm2之間，其走向?yàn)橄冉岛笊?，總體略呈上升趨勢(shì)(期末值相對(duì)期初值增加2.36%，見(jiàn)圖7)。

當(dāng)以化學(xué)需氧量作為評(píng)價(jià)對(duì)象時(shí)，經(jīng)濟(jì)區(qū)除2009年外均表現(xiàn)為水環(huán)境生態(tài)盈虧，其人均盈虧值居于-0.005 2-0.515 8 hm2之間(見(jiàn)圖7)，說(shuō)明該地區(qū)有機(jī)污染物對(duì)水體的壓力基本在其承載力范圍內(nèi)；當(dāng)以氮、磷作為評(píng)價(jià)對(duì)象時(shí)，經(jīng)濟(jì)區(qū)歷年均呈赤字狀態(tài)，其人均赤字值分別介于0.490 1-0.890 9 hm2之間和0.131 4-0.515 1 hm2之間(見(jiàn)圖7)，說(shuō)明該地區(qū)氮、磷排放對(duì)當(dāng)?shù)厮h(huán)境的壓力已超過(guò)其承載力。由于氮污染水生態(tài)足跡是三種污染足跡中的最大值，故最終水環(huán)境生態(tài)盈虧狀況由該評(píng)價(jià)對(duì)象體現(xiàn)。由此可以認(rèn)為，經(jīng)濟(jì)區(qū)歷年水環(huán)境均處于不安全狀態(tài)。其走向與人均水環(huán)境承載力基本一致，這主要是因?yàn)榻?jīng)濟(jì)區(qū)歷年人均水污染生態(tài)足跡變化相對(duì)較小而其歷年人均水環(huán)境承載力變化較大。

圖7 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)人均水環(huán)境承載力、水環(huán)境生態(tài)盈虧和水環(huán)境生態(tài)壓力指數(shù)

當(dāng)以化學(xué)需氧量作為評(píng)價(jià)對(duì)象時(shí)，經(jīng)濟(jì)區(qū)歷年水環(huán)境生態(tài)壓力指數(shù)介于0.373 6-1.011 6之間(見(jiàn)圖7)，其中2009年和2010年經(jīng)濟(jì)區(qū)水環(huán)境處于稍不安全和臨界水平，其余年份則為很安全或較安全水平，說(shuō)明該地區(qū)有機(jī)污染物排放對(duì)當(dāng)?shù)厮w的壓力相對(duì)較??；當(dāng)以氮、磷作為評(píng)價(jià)對(duì)象時(shí)，該指數(shù)均超過(guò)1(見(jiàn)圖7)，說(shuō)明該地區(qū)氮、磷排放使當(dāng)?shù)厮h(huán)境處于不安全狀態(tài)。基于前述，最終水環(huán)境生態(tài)壓力狀況由氮污染足跡估算結(jié)果反映，由此可以認(rèn)為，經(jīng)濟(jì)區(qū)近年來(lái)水環(huán)境安全度較差(較不安全和很不安全年份各占一半)。因?yàn)闅v年人均水污染生態(tài)足跡變幅相對(duì)較小，所以水環(huán)境生態(tài)壓力指數(shù)變化趨勢(shì)主要由人均水環(huán)境承載力決定。

3.2.2 廣西北部灣經(jīng)濟(jì)區(qū)水環(huán)境承載力、生態(tài)盈虧及生態(tài)壓力空間分析

2008-2013年，各市人均水環(huán)境承載力介于0.230 4-3.378 0 hm2之間，其中，防城港市人均水環(huán)境承載力很高(年均值高達(dá)2.7hm2)，而其他3市則相對(duì)低很多(年均值介于0.4-0.7 hm2之間)(見(jiàn)圖8)，說(shuō)明經(jīng)濟(jì)區(qū)內(nèi)部水環(huán)境承載力非常不平衡。從走向來(lái)看，各市波動(dòng)明顯，說(shuō)明其水環(huán)境承載力不穩(wěn)定。從圖8、表6可以看出，各市人均水環(huán)境承載力高低與其經(jīng)濟(jì)社會(huì)發(fā)展水平關(guān)聯(lián)度較小，而與其水資源豐裕度及人口密度相關(guān)性較大。

圖8 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)各市人均水環(huán)境承載力

行政分區(qū)多年平均人均GDP(元/人)多年平均人均水資源量(m3/人)水資源平均生產(chǎn)能力(m3/hm2)南寧1657416766260北海2252018769670防城港23479901512334欽州1033722509834

注：人均GDP和人均水資源量原始數(shù)據(jù)來(lái)自《廣西統(tǒng)計(jì)年鑒》；水資源平均生產(chǎn)能力數(shù)據(jù)來(lái)自《廣西水功能區(qū)修訂報(bào)告》[23]

各市人均水環(huán)境生態(tài)盈虧值處于-1.329 0-2.005 4 hm2之間。其中，防城港市歷年水環(huán)境均表現(xiàn)為生態(tài)盈余，說(shuō)明該地區(qū)水環(huán)境安全狀況很好；其余地區(qū)所有年份均呈赤字水平，說(shuō)明它們的水環(huán)境安全狀況較差(見(jiàn)圖9)。

圖9 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)各市人均水環(huán)境生態(tài)盈虧

各市水環(huán)境生態(tài)壓力指數(shù)居于0.401 8-5.621 7之間(見(jiàn)圖10)?；谠撝笖?shù)，防城港市有3個(gè)年份表現(xiàn)為水環(huán)境很安全水平，2個(gè)年份為較安全水平，其水環(huán)境生態(tài)壓力指數(shù)最大值僅0.809 2，說(shuō)明它的水環(huán)境安全狀況很好，這主要是因?yàn)樗乃h(huán)境承載力很高。欽州市和北海市水環(huán)境安全狀況較差(欽州市僅2年為稍不安全，其余年份為很不安全或較不安全；北海市3年為很不安全，3年為較不安全)，這是因?yàn)椋核鼈兊乃廴旧鷳B(tài)足跡雖相對(duì)較低，但其水環(huán)境承載力更低。南寧市所有年份均為很不安全(其水環(huán)境生態(tài)壓力指數(shù)年均值高達(dá)3.9)，說(shuō)明其水環(huán)境安全狀況很差，這是由于該市水污染生態(tài)足跡較高同時(shí)其水環(huán)境承載力又很低。

圖10 2008-2013年廣西北部灣經(jīng)濟(jì)區(qū)各市水環(huán)境生態(tài)壓力指數(shù)

4 結(jié)論與建議

2008-2013年，廣西北部灣經(jīng)濟(jì)區(qū)及其各市均為氮污染水生態(tài)足跡最大，磷污染水生態(tài)足跡次之，有機(jī)物污染水生態(tài)足跡最小。該地區(qū)有機(jī)物污染足跡呈下降趨勢(shì)，前3年最大驅(qū)動(dòng)因子為工業(yè)賬戶(hù)，后3年則為城鎮(zhèn)生活賬戶(hù)；氮/磷污染足跡總體表現(xiàn)為上升，其最大影響因子來(lái)自畜禽養(yǎng)殖賬戶(hù)。各市3種水污染生態(tài)足跡大小排序不一致，其變化趨勢(shì)也各異，說(shuō)明經(jīng)濟(jì)區(qū)水污染生態(tài)足跡空間差異明顯。經(jīng)濟(jì)區(qū)水環(huán)境承載力先降后升，總體略呈上升趨勢(shì)；內(nèi)部非常不平衡，其中，防城港市該值很大，其他3市則較小?；谒h(huán)境生態(tài)盈虧和水環(huán)境生態(tài)壓力指數(shù)，經(jīng)濟(jì)區(qū)水環(huán)境安全度較差，其內(nèi)部差異顯著：防城港市很好，欽州市和北海市較差，南寧市很差。

為改善經(jīng)濟(jì)區(qū)水環(huán)境狀況，首先應(yīng)大幅降低該地區(qū)畜禽養(yǎng)殖業(yè)氮、磷排放量，這要求：1)控制畜禽養(yǎng)殖業(yè)規(guī)模(尤其是南寧市)；2)合理調(diào)整養(yǎng)殖結(jié)構(gòu)，例如減少牛養(yǎng)殖量，擴(kuò)大羊養(yǎng)殖規(guī)模；3)各地要科學(xué)劃定畜禽養(yǎng)殖禁養(yǎng)區(qū)，并依法關(guān)閉和搬遷禁養(yǎng)區(qū)內(nèi)的畜禽養(yǎng)殖場(chǎng)(小區(qū))和養(yǎng)殖專(zhuān)業(yè)戶(hù)；4)各地要根據(jù)養(yǎng)殖方式、規(guī)模和自身?xiàng)l件采用適宜的防治技術(shù)。其次，應(yīng)提高經(jīng)濟(jì)區(qū)城鎮(zhèn)生活污水污染防治能力。為此，應(yīng)加強(qiáng)城鎮(zhèn)污水處理廠(chǎng)配套管網(wǎng)工程建設(shè)，對(duì)于新建污水處理設(shè)施，必須“廠(chǎng)網(wǎng)并舉，管網(wǎng)先行”；已建和在建的城鎮(zhèn)污水處理廠(chǎng)要加快管網(wǎng)建設(shè)，提高城鎮(zhèn)污水處理率?，F(xiàn)有合流制排水系統(tǒng)應(yīng)加快實(shí)施雨污分流改造，難以改造的，應(yīng)采取截流、調(diào)蓄和治理等措施。再次，加大工業(yè)廢水防治力度：加強(qiáng)造紙、制糖、酒精、淀粉、紡織、制革、有色金屬、電鍍等主要污染企業(yè)的治理，實(shí)行更嚴(yán)格的淘汰制度；新建以接納工業(yè)廢水為主的集中式污水處理廠(chǎng)必須配套建設(shè)脫氮除磷設(shè)施，已建的污水處理廠(chǎng)按新的排放限值進(jìn)行提標(biāo)改造。另外，隨著該地區(qū)農(nóng)民生活水平的提高，農(nóng)村生活污水排放量及其中污染物濃度也隨之上升，故應(yīng)建設(shè)適宜的污水收集系統(tǒng)，采用最佳可行技術(shù)實(shí)現(xiàn)污水的資源化利用。

相對(duì)于傳統(tǒng)水污染生態(tài)模型，改進(jìn)模型的優(yōu)點(diǎn)主要有：比較全面的核算了研究區(qū)水污染生態(tài)足跡；比較準(zhǔn)確的衡量了研究區(qū)水環(huán)境承載力；避免了均衡因子取值困難；水污染生態(tài)足跡只與水環(huán)境承載力比較，解決了賬戶(hù)匯總計(jì)算的問(wèn)題。但該模型仍存在某些不足：核算項(xiàng)目仍不完整。例如，沒(méi)有考慮船舶等流動(dòng)污染源和城鄉(xiāng)垃圾對(duì)水體的污染，也沒(méi)有估算水體對(duì)溫室氣體的吸納功能；由于問(wèn)題的復(fù)雜性以及人力物力財(cái)力不足，本研究?jī)H從稀釋角度估算研究區(qū)水環(huán)境承載力，而實(shí)際上水環(huán)境承載力大小與水域特性(如幾何特征、水文特征、物理自?xún)裟芰?、化學(xué)自?xún)裟芰?、生物降?、水環(huán)境功能要求、污染物質(zhì)特性、排污方式等因素密切相關(guān)[20]5；本研究采用的水環(huán)境生態(tài)壓力等級(jí)標(biāo)準(zhǔn)主要參考他人相關(guān)研究成果，這雖具有一定合理性但有必要研究確定適用于該模型的等級(jí)標(biāo)準(zhǔn)。這些問(wèn)題的解決有待于進(jìn)一步深入研究。

[1]王書(shū)華.區(qū)域生態(tài)經(jīng)濟(jì)——理論、方法與實(shí)踐[M].北京：中國(guó)發(fā)展出版社，2008.

[2]黃林楠，張偉新，姜翠玲，等.水資源生態(tài)足跡計(jì)算方法[J].生態(tài)學(xué)報(bào)，2008，28(3)：1279-1286.

[3]王文國(guó)，何明雄，潘科，等.四川省水資源生態(tài)足跡與生態(tài)承載力的時(shí)空分析[J].自然資源學(xué)報(bào)，2011，26(9)：1555-1565.

[4]冷建飛，張倩.基于生態(tài)足跡模型的水資源可持續(xù)利用研究——以江蘇省為例[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2015，29(7)：54-59.

[5]洪輝.基于生態(tài)足跡法的西安市水資源生態(tài)足跡研究[D].西安：西安建筑科技大學(xué)，2007.

[6]張義，張合平.基于生態(tài)系統(tǒng)服務(wù)的廣西水生態(tài)足跡分析[J].生態(tài)學(xué)報(bào)，2013，33(13)：4111-4124．

[7]段錦，康慕誼，江源.基于淡水資源賬戶(hù)和污染賬戶(hù)的生態(tài)足跡改進(jìn)模型[J].自然資源學(xué)報(bào)，2012，27(6)：953-963.

[8]于冰，徐琳瑜.城市水生態(tài)系統(tǒng)可持續(xù)發(fā)展評(píng)價(jià)——以大連市為例[J].資源科學(xué)，2014，36(12)：2578-2583.

[9]白鈺，曾輝，魏建兵，等.基于環(huán)境污染賬戶(hù)核算的生態(tài)足跡模型優(yōu)化——以珠江三角洲城市群為例[J].應(yīng)用生態(tài)學(xué)報(bào)，2008，19(8)：1789-1796.

[10]閔慶文，焦雯珺，成升魁.污染足跡：一種基于生態(tài)系統(tǒng)服務(wù)的生態(tài)足跡[J].資源科學(xué)，2011，33(2)：195-200.

[11]焦雯珺，閔慶文，李文華，等.基于ESEF的水生態(tài)承載力：理論、模型與應(yīng)用[J].應(yīng)用生態(tài)學(xué)報(bào)，2015，26(4)：1041-1048.

[12]李芬，孫然好，楊麗蓉，等.基于供需平衡的北京地區(qū)水生態(tài)服務(wù)功能評(píng)價(jià)[J].應(yīng)用生態(tài)學(xué)報(bào)，2010，21(5)：1146-1152.

[13]焦雯珺，閔慶文，成升魁，等.基于污染足跡模型的太湖流域水環(huán)境演變的人文驅(qū)動(dòng)力評(píng)估——以江蘇省常州市為例[J].資源科學(xué)，2011，33(2)：223-229.

[14]趙先貴，韋良煥，馬彩虹，等.西安市生態(tài)足跡與生態(tài)安全的動(dòng)態(tài)研究[J].干旱區(qū)資源與環(huán)境，2007，21(1)：1-5.

[15]宋旭光.生態(tài)占用測(cè)度問(wèn)題研究[J].統(tǒng)計(jì)研究，2003，20(2)：44-47.

[16]陳成忠.生態(tài)足跡模型的多尺度分析及其預(yù)測(cè)研究[D].南京：南京師范大學(xué)，2008.

[17]國(guó)家環(huán)境保護(hù)總局.地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB3828-2002)[S].北京：中國(guó)標(biāo)準(zhǔn)出版社，2002.

[18]環(huán)境保護(hù)部華南環(huán)境科學(xué)研究所.生活源產(chǎn)排污系數(shù)及使用說(shuō)明[EB/OL].http：//www.eiafans.com/thread-263005-1-1.html，2011-9-28.

[19]國(guó)家環(huán)境保護(hù)總局.城鎮(zhèn)污水處理廠(chǎng)污水排放標(biāo)準(zhǔn)(GB18918-2002)[S].北京：中國(guó)環(huán)境科學(xué)出版社，2002.

[20]逄勇，陸桂華.水環(huán)境容量計(jì)算理論及應(yīng)用[M].北京：科學(xué)出版社，2010.

[21]張義，張合平，李豐生，等.基于改進(jìn)模型的廣西水資源生態(tài)足跡動(dòng)態(tài)分析[J].資源科學(xué)，2013，35(8)：1601-1610.

[22]國(guó)家環(huán)境保護(hù)總局自然生態(tài)保護(hù)司.全國(guó)規(guī)?；笄蒺B(yǎng)殖業(yè)污染情況調(diào)查及防治對(duì)策[M].北京：中國(guó)環(huán)境科學(xué)出版社，2002.

[23]廣西水利廳.廣西水功能區(qū)劃修訂報(bào)告[R].南寧：廣西水利廳，2013.

[責(zé)任編輯 韋志巧]

Spatio-temporal Analysis of Water Pollution Ecological Footprint of Guangxi North Gulf Economic Zone Based on an Improved Model

ZHANG Yi, CHENG Lianting

(College of Economics and Management, Hechi University, Yizhou， Guangxi 546300, China)

The current water pollution ecological footprint model had a big problem. Thus, this papershows an improved model, and takes the Guangxi North Gulf Economic Zone in 2008-2013 as an example for calculation and analysis. The result showed that: (1) The ecological footprint of the economic zone and its cities were mostly polluted by nitrogen, and the secondly by phosphorus, and organic pollutants were the least. The organic pollution footprint decreased, and the nitrogen/phosphorus pollution footprint showed increase trend. In these cities, the size of the 3 kinds of water pollution ecological footprint did not show consistency, and theirchange trends were also varied. (2) The water environment carrying capacity of the economic zone first dropped and then rose, and overall showed upward slight trends. Its interior was very unbalance, and the water environment carrying capacity of Fangchenggang city was high, but the other 3 cities were very low. (3) Based on water environment ecological surplus/deficit and water environment ecological pressure index, water environment security of the economic zone was poor, but there are significant differences in water environment security among different cities. Water environment security in Fangchenggang city is good, but water environment security in Qinzhou city and Beihai city are worse than that in Fangchenggang city, and the Nanning city is the worst.

water pollutant ecological footprint; water environment; carrying capacity; ecological surplus/deficit; water environment ecological pressure index; Guangxi North Gulf Economic Zone

X824

A

1672-9021(2017)02-0001-13

張義(1978-)，男，湖南湘潭人，河池學(xué)院經(jīng)濟(jì)與管理學(xué)院副教授，博士，主要研究方向：生態(tài)(資源)經(jīng)濟(jì)學(xué)，生態(tài)安全。

廣西哲學(xué)社會(huì)科學(xué)規(guī)劃一般項(xiàng)目(13BJY022)。

2017-01-16