應(yīng)對畜禽養(yǎng)殖跨界污染的流域水資源生態(tài)補償量研究
——以國家生態(tài)補償示范區(qū)九洲江流域為例

趙素芹,孫 翔,*, 侯東林,朱 燚

1 廣西大學(xué)資源環(huán)境與材料學(xué)院, 南寧 530004 2 生態(tài)環(huán)境部環(huán)境與經(jīng)濟(jì)政策研究中心, 北京 100029

隨著九洲江流域生豬等養(yǎng)殖業(yè)的迅猛發(fā)展,畜禽養(yǎng)殖污染日趨嚴(yán)重,對當(dāng)?shù)亓饔蛩h(huán)境造成巨大壓力,畜禽養(yǎng)殖業(yè)成為流域污染的主要來源之一。流域水資源生態(tài)補償是平衡上下游發(fā)展與保護(hù)關(guān)系、解決跨區(qū)域污染的有效途徑,亦是環(huán)境經(jīng)濟(jì)學(xué)、環(huán)境管理學(xué)領(lǐng)域研究的熱點問題。自20世紀(jì)90年代末生態(tài)補償理論被引入到流域管理領(lǐng)域,國內(nèi)外專家、學(xué)者對流域生態(tài)補償進(jìn)行了積極的研究探索[1- 3]。目前,國內(nèi)外學(xué)者針對流域生態(tài)補償?shù)难芯恐饕性谏鷳B(tài)補償理論內(nèi)涵[4- 6]、補償模式[7- 9]、補償標(biāo)準(zhǔn)[10- 13]、支付意愿[14- 19]等方面。流域生態(tài)補償機(jī)制理論的研究主要從跨省流域生態(tài)補償、省內(nèi)跨市流域生態(tài)補償、市域內(nèi)跨縣流域生態(tài)補償?shù)?個層次展開,剖析涵義、原則、標(biāo)準(zhǔn)、驅(qū)動力與障礙、政策服務(wù)等[20]。在補償模式與途徑設(shè)計方面,主要提出政府主導(dǎo)型和市場交易型模式,且提出了各自應(yīng)用的條件與準(zhǔn)則,剖析了現(xiàn)有模式的局限問題,建立了生態(tài)補償制度選擇和操作范式[21]?，F(xiàn)階段,流域生態(tài)補償標(biāo)準(zhǔn)的測算主要從投入成本、環(huán)境效益、補償意愿等角度切入[22],主要方法包括條件價值法、機(jī)會成本法、價值和費用分析法等[23],這些核算方法為流域生態(tài)補償標(biāo)準(zhǔn)的確定提供依據(jù)。而有關(guān)養(yǎng)殖污染主控流域管理方面研究主要集中在養(yǎng)殖污染治理成本效益、養(yǎng)殖戶偏好、管理模式、養(yǎng)殖模式、受償意愿等[24- 26]。有些學(xué)者探討了畜禽養(yǎng)殖治理和管理模式優(yōu)化相結(jié)合的定量研究[27]；孫若梅[28]探討了畜禽養(yǎng)殖業(yè)生態(tài)補償內(nèi)涵的界定,并提出具體的生態(tài)補償措施；張郁等[29]分析了生態(tài)補償情景下家庭資源稟賦對養(yǎng)豬戶的環(huán)境行為影響；賓幕容等[24]以湘江流域367位生豬養(yǎng)殖戶的實地調(diào)查數(shù)據(jù)為基礎(chǔ),對農(nóng)戶生豬養(yǎng)殖污染治理意愿的影響因素及層次結(jié)構(gòu)進(jìn)行了分析；昌苗苗和陳洋[30]對九洲江流域生態(tài)補償實踐進(jìn)展進(jìn)行了介紹；王西琴等[31]以九洲江為研究區(qū)域,計算了限制生豬養(yǎng)殖的機(jī)會成本與環(huán)境效益,建立了基于畜禽養(yǎng)殖模式轉(zhuǎn)變的生態(tài)補償空間優(yōu)化模型；唐潥等[32]以九洲江流域生豬養(yǎng)殖戶為研究對象,采用意愿調(diào)查法分別給出禁養(yǎng)與升級改造受償意愿。這些理論與實踐成果為本研究的開展奠定了很好的方法學(xué)基礎(chǔ),但在養(yǎng)殖污染主控流域水資源生態(tài)補償量的分析框架與技術(shù)方法上仍有較大的發(fā)展空間。本文以國家生態(tài)補償示范區(qū)九洲江流域作為典型案例,基于Ⅲ類水質(zhì)目標(biāo)要求計算污染物削減目標(biāo)量,設(shè)置禁養(yǎng)區(qū)關(guān)閉拆遷養(yǎng)豬場、高架床養(yǎng)殖模式改造、養(yǎng)殖糞尿廢水集中處理三種典型控制方式,以污染物削減目標(biāo)為約束條件,基于費效比構(gòu)建線性規(guī)劃模型,求算用于水污染控制的最小政府補償總額,并以取水量、用水量、GDP為主要約束條件,建立上下游地方政府及中央政府生態(tài)補償金額分配模型。以期為流域生態(tài)補償政策的制定提供參考和借鑒,對其他流域生態(tài)補償量的確定提供方法學(xué)參考。

1 研究區(qū)域概況

九洲江是一條跨粵桂兩省(區(qū))的獨流入海河流,發(fā)源于廣西陸川縣沙坡鎮(zhèn)秦鏡村的文龍徑,流經(jīng)廣西玉林市陸川、博白兩縣共10個鎮(zhèn)進(jìn)入廣東湛江市的鶴地水庫,最后注入北部灣,全長168 km,總集雨面積3396 km2。鶴地水庫位于湛江廉江市河唇鎮(zhèn),是九洲江流域下游湛江市重要的飲用水源地,水面面積共122 km2。近年來,以生豬養(yǎng)殖主控的畜禽養(yǎng)殖業(yè)是流域上游廣西玉林市的主要發(fā)展產(chǎn)業(yè),養(yǎng)殖源成為九洲江流域水體污染的主要來源,養(yǎng)殖污染排放中化學(xué)需氧量、氨氮、總氮、總磷占總量比例分別為60.69%、57.83%、66.00%、78.51%[25]?；浌饍墒?區(qū))黨委、政府和社會各界高度重視九洲江流域水環(huán)境安全,2015年3月,粵桂兩省聯(lián)合出臺《九洲江流域水環(huán)境補償實施方案》；同年9月,中共中央、國務(wù)院明確將九洲江流域環(huán)境綜合治理列為國家《生態(tài)文明體制改革總體方案》跨地區(qū)生態(tài)補償三個試點之一。2016年3月,兩省簽訂《九洲江流域上下游橫向生態(tài)補償協(xié)議》,在兩省及中央的共同努力下,流域水環(huán)境已經(jīng)得到初步改善?；诹饔蚋鬣l(xiāng)鎮(zhèn)生豬養(yǎng)殖現(xiàn)狀,欲通過3個方面進(jìn)行流域畜禽污染治理工作(圖1)：一是根據(jù)《粵桂兩省九洲江水污染防治規(guī)劃》,劃定禁養(yǎng)區(qū)(干流沿岸500 m及支流沿岸200 m范圍內(nèi))、限養(yǎng)區(qū)(干流500—2000 m及支流沿岸200—2000 m范圍內(nèi))及適養(yǎng)區(qū)(其他區(qū)域),實施生豬養(yǎng)殖精細(xì)化分區(qū)管理；二是推廣生態(tài)養(yǎng)殖新模式,實現(xiàn)高架床養(yǎng)殖；三是進(jìn)行養(yǎng)殖廢棄物處理設(shè)施建設(shè),實現(xiàn)廢棄物集中處理。因此,如何用最小的生態(tài)補償穩(wěn)定達(dá)到水質(zhì)要求,取得環(huán)境效益目標(biāo)的同時兼顧九洲江上游廣西玉林片區(qū)為保護(hù)水質(zhì)所作的經(jīng)濟(jì)犧牲,是九洲江流域畜禽養(yǎng)殖污染治理和生態(tài)補償政策制定面臨的重要問題。

圖1 九洲江流域養(yǎng)殖區(qū)域劃分及主要污染處理管控方式Fig.1 Spatial distribution of pig farms, and major pollution treatment methods for livestock and poultry in Jiuzhou River Basin

2 流域生態(tài)補償量計算方法

流域生態(tài)補償額度的研究,從最初定性確定補償理論機(jī)制,到定量研究具體補償量,方法逐漸明確科學(xué),也更加具有合理性。當(dāng)前,投入、效益、補償意愿等3方面是絕大多數(shù)流域生態(tài)補償額度測算的主題方向。在投入方面,通過計算為保護(hù)或改善水資源環(huán)境狀況,所進(jìn)行的各項投入,包括保護(hù)環(huán)境而喪失的機(jī)會成本。主要方法包括機(jī)會成本法[33]、生態(tài)重建成本分?jǐn)偡╗34]、損失價值核算法[35]等；在效益方面,通過估算保護(hù)環(huán)境投入在經(jīng)濟(jì)、社會、生態(tài)等方面產(chǎn)生的外部效應(yīng)來確定補償值,主要有生態(tài)系統(tǒng)服務(wù)價值法、水質(zhì)水量法、污染權(quán)等[36]方法；在補償意愿方面,通過對消費者進(jìn)行直接調(diào)查和詢問,了解消費者的對于改善或保護(hù)環(huán)境的支付意愿,主要有條件價值法、博弈法等[37- 38]。

補償資金額度的測算是實現(xiàn)生態(tài)補償?shù)那疤?也是生態(tài)補償?shù)年P(guān)鍵環(huán)節(jié)?；诂F(xiàn)有研究成果,本研究基于水質(zhì)水量的跨界區(qū)域生態(tài)補償量計算方法,結(jié)合九洲江流域的主要特點,建立上下游地區(qū)及中央生態(tài)補償標(biāo)準(zhǔn),推動流域生態(tài)補償具體實施與貫徹執(zhí)行。

2.1 水環(huán)境目標(biāo)約束下的流域生態(tài)補償量

2.1.1流域水環(huán)境污染物削減目標(biāo)總量計算

水環(huán)境容量是在特定水環(huán)境目標(biāo)約束下,綜合考慮水體自凈能力下能夠容納的最大排污量[34- 35]。模型選取COD、總磷、總氮為指標(biāo),在Ⅲ類水質(zhì)要求下,基于水環(huán)境容量為約束條件計算畜禽污染物削減目標(biāo)總量。以下為削減目標(biāo)總量模型函數(shù)：

Q削=Q入河-W

(1)

Q入河=365Ciq1(T1r1+T2r2)

(2)

W=Q0(Cs-C0)+KVCs

(3)

式中,Q削為污染物削減目標(biāo)量；Q入河為畜禽養(yǎng)殖污染物現(xiàn)狀排放入河量；W為水環(huán)境容量；Ci為畜禽污水污染物濃度(見表1),mg/L；T1、T2為禁養(yǎng)區(qū)、限養(yǎng)區(qū)生豬存欄量,頭；q1為傳統(tǒng)養(yǎng)殖模式日排水量,0.03 t 頭-1d-1(實地養(yǎng)殖場調(diào)研求均值所得,方差δ2為0.001)；r1、r2為禁養(yǎng)區(qū)、限養(yǎng)區(qū)污染物入河系數(shù)(根據(jù)文獻(xiàn)查閱[39-41],流域內(nèi)禁養(yǎng)區(qū)范圍污染物入河系數(shù)在0.60—0.70,限養(yǎng)區(qū)范圍內(nèi)污染物入河系數(shù)在0.30—0.60?；诰胖藿饔?qū)嵉卣{(diào)查(圖2),禁養(yǎng)區(qū)范圍內(nèi)養(yǎng)殖密度大,取最大入河系數(shù)0.70,限養(yǎng)區(qū)范圍內(nèi)取均值0.50)；Q0為流域多年平均徑流量,9.28億m3/a；Cs為污染物的環(huán)境標(biāo)準(zhǔn)值(Ⅲ類)；C0為污染物的環(huán)境背景值,這里忽略本底值的影響,假設(shè)C0= 0；V為流域水資源體積,11.05億m3；K為降解系數(shù),根據(jù)文獻(xiàn)現(xiàn)有研究[42-43],確定COD降解系數(shù)為0.20 d-1；TP、TN降解系數(shù)為0.10 d-1。

表1 生豬養(yǎng)殖廢水中污染物質(zhì)量濃度

備注：數(shù)據(jù)來源于《畜禽養(yǎng)殖業(yè)污染治理工程技術(shù)規(guī)范(HJ497—2009)》; COD:化學(xué)需氧量, Chemical Oxygen Demand; TP: 總磷, Total Phosphorus; TN: 總氮, Total Nitrogen

圖2 九洲江流域調(diào)研路線圖Fig.2 Survey roadmap in Jiuzhou River Basin

2.1.2污染物削減分配優(yōu)化模型構(gòu)建

在污染物削減目標(biāo)確定的基礎(chǔ)上,總削減量的分配采取整型規(guī)劃及優(yōu)先次序遞推法,禁養(yǎng)區(qū)內(nèi)禁養(yǎng)清拆為第一步；限養(yǎng)區(qū)處理方式為第二步,并設(shè)高架床養(yǎng)殖模式改造和污染物集中處理兩種情景。按照先后順序?qū)B(yǎng)區(qū)、限養(yǎng)區(qū)進(jìn)行可削減量估算,比較總目標(biāo)削減量和各部分可削減量,若禁養(yǎng)區(qū)可滿足總削減量要求,則完成削減目標(biāo)；若禁養(yǎng)區(qū)不能滿足總削減量要求,則遞推至限養(yǎng)區(qū),直到累積可削減量達(dá)到總目標(biāo)削減量值為止[44](圖3)。

圖3 污染物削減量分配模型Fig.3 Optimization model of pollutant reduction allocation

2.1.3生態(tài)補償量計算

限養(yǎng)區(qū)內(nèi)有兩種畜禽養(yǎng)殖污染治理方案,不同治理方案的治理成本存在差異,因此本文基于限養(yǎng)區(qū)不同污染治理方案設(shè)置不同成本計算情景,從而選取最小生態(tài)補償量。根據(jù)玉林市環(huán)保局出臺的《粵桂兩省九洲江水污染防治規(guī)劃》,其限養(yǎng)區(qū)污染治理分高架床養(yǎng)殖模式改造、養(yǎng)殖糞污集中式處理點建設(shè)兩種方案,因此生態(tài)補償成本計算可設(shè)置兩種情景,即“禁養(yǎng)+高架床生態(tài)養(yǎng)殖改造”為情景一,“禁養(yǎng)+集中式處理設(shè)施建設(shè)”為情景二,以最小生態(tài)補償金額為目標(biāo)函數(shù),以污染物削減目標(biāo)量為主要約束條件,得到目標(biāo)函數(shù)如下：

MinY={Y1,Y2}

(4)

(5)

情景二：Y2=F1T1+(Q削-Q禁)F3/Ci

(6)

Q禁=365T1q1Cir1

(7)

Q高架床=365(q1-q2)Cir2

(8)

式中,Y為最小目標(biāo)生態(tài)補償量,元；Y1、Y2為不同情景下的生態(tài)補償成本,元；Q禁為禁養(yǎng)區(qū)內(nèi)可削減污染物量,t；Q高架床為高架床養(yǎng)殖模式改造可削減污染物量,t 頭-1a-1;M為高架床養(yǎng)殖模式豬舍容積比,1.10 m2/頭(兩次實地調(diào)研取均值所得,方差δ2為0.11)；q2為高架床養(yǎng)殖模式下生豬養(yǎng)殖排水量,0.0018 t 頭-1d-1(兩次調(diào)研取均值所得,方差δ2為0.0002)；F1、F2分別為禁養(yǎng)區(qū)內(nèi)清拆禁養(yǎng)、限養(yǎng)區(qū)內(nèi)高架床養(yǎng)殖改造治理策略下的補償金額,分別為500.00、440.00元/ m2,數(shù)據(jù)參考于當(dāng)?shù)卣卫硪?guī)劃文本報告[45],且實地調(diào)查進(jìn)行驗證；F3為集中式養(yǎng)殖糞污處理模式下補償金額,100.00元/t,數(shù)據(jù)來自博世科陸川縣固廢新能源項目實地考察,此項目是流域上游地區(qū)唯一的養(yǎng)殖糞污資源化處理的項目點,數(shù)據(jù)參考于當(dāng)?shù)卣畬Υ隧椖康难a償標(biāo)準(zhǔn)。

2.2 九洲江流域生態(tài)補償量分擔(dān)模型構(gòu)建

圍繞九洲江水環(huán)境污染的現(xiàn)狀和特點,粵桂雙方多次協(xié)商。2014年8月,廣東省人民政府和廣西壯族自治區(qū)人民政府聯(lián)合簽署《九洲江流域跨界水環(huán)境保護(hù)合作協(xié)議》；2016年3月21 日,在環(huán)保部、財政部的指導(dǎo)下,廣東、廣西兩省區(qū)在福建省簽訂了《九洲江流域上下游橫向生態(tài)補償協(xié)議》。兩省雖有較好的合作基礎(chǔ),然而九洲江流域上下游及中央補償比例還是存在紛爭。本研究結(jié)合流域上下游用水量以及經(jīng)濟(jì)承受能力,引入中央政府財政縱向轉(zhuǎn)移支付,建立地區(qū)及中央生態(tài)補償成本分?jǐn)偰Ｐ?以期為九洲江養(yǎng)殖水污染生態(tài)補償?shù)膶嵤┨峁┛茖W(xué)建議。具體模型如下：

(9)

(10)

Y廣東=Y-Y中央-Y廣西

(11)

式中,Y為流域生態(tài)補償量,元；Y中央為中央政府縱向財政轉(zhuǎn)移支付,元；Y廣西為廣西省生態(tài)補償成本分?jǐn)偭?元；Y廣東為廣東省生態(tài)補償成本分?jǐn)偭?元；Q廣西為上游廣西取水量,m3；Q供水量為九洲江流域供水量,m3；G上游為廣西玉林市GDP總值(2016年),億元；G上下游為上下游(廣西玉林、廣東湛江)GDP總值,億元。

2.3 參數(shù)靈敏度計算

考慮到模型參數(shù)的取值容易受到外界環(huán)境的影響,且具有不確定性。本文利用Crystal Ball軟件進(jìn)行模型參數(shù)的靈敏度分析。首先采用蒙特卡羅模擬方法對模型參數(shù)的概率分布函數(shù)進(jìn)行大量的重復(fù)抽樣,代入式(2)、式(3)、式(5)、式(6)建立生態(tài)補償量的分布函數(shù),從而獲得模型中重要的數(shù)學(xué)特征,如數(shù)學(xué)期望、方差、區(qū)間估計等。在此步驟基礎(chǔ)上,利用方差分析法對參數(shù)進(jìn)行方差分析,根據(jù)參數(shù)輸入對生態(tài)補償量輸出的秩相關(guān)系數(shù)的平方歸一化后所占的比例來確定參數(shù)的靈敏度[46]。具體公式如下所示：

(12)

(13)

式中,γk為k參數(shù)變量的方差貢獻(xiàn)率；αk0表示k參數(shù)變量對生態(tài)補償量的秩相關(guān)系數(shù)；n表示樣本數(shù)；Rik表示參數(shù)變量k第i次抽樣在樣本中的排序；Ri0表示第i次抽樣對應(yīng)的生態(tài)補償量的結(jié)果排序。

3 結(jié)果分析

3.1 污染物削減目標(biāo)總量

根據(jù)公式(3)計算可知,COD、TP、TN水環(huán)境容量分別為20682.00、186.93、934.65 t/a。2016年九洲江流域生豬養(yǎng)殖頭數(shù)為311.00萬頭,其中禁養(yǎng)區(qū)范圍內(nèi)有66.90萬頭,限養(yǎng)區(qū)生豬存欄量為244.10萬頭。在Ⅲ類水質(zhì)要求下,流域COD、TN、TP污染物削減目標(biāo)總量分別為378752.98、13951.70 t和2161.60 t。

3.2 生態(tài)補償量及分?jǐn)偙壤?/h3>

根據(jù)公式(5)、(6)計算可知,其情景一(“禁養(yǎng)+高架床生態(tài)養(yǎng)殖改造”模式)測算出生態(tài)補償金額為15.03億元；情景二(“禁養(yǎng)+集中式污水處理設(shè)施建設(shè)”模式)測算出生態(tài)補償金額為15.55億元。按照成本最小化補償原則確定流域生態(tài)補償量為15.03億元。九洲江流域多年平均總供水量為13.97億m3,其向廣西、廣東境內(nèi)供水量分別為3.00億m3和10.97億m3。據(jù)2016年統(tǒng)計公報顯示,廣西玉林市GDP為1553.91億元,廣東湛江市GDP為2584.78億元。經(jīng)計算,中央、廣西、廣東分擔(dān)率分別為0.25、0.22和0.53,中央、廣西、廣東生態(tài)補償分擔(dān)金額分別為3.76、3.31、7.96億元。

3.3 參數(shù)靈敏度分析

本文采用Crystal Ball軟件進(jìn)行蒙特卡羅模擬和靈敏度分析,模型的輸入設(shè)置了高架床養(yǎng)殖容積比(k1)、限養(yǎng)區(qū)污水入河系數(shù)(k2)、高架床養(yǎng)殖模式日排水量(k3)及傳統(tǒng)養(yǎng)殖模式日排水量(k4)等4個參數(shù)變量。預(yù)測變量為生態(tài)補償量。采用蒙特卡羅法進(jìn)行5000次重復(fù)抽樣,得到生態(tài)補償量預(yù)測值得直方圖。95%置信區(qū)間水平下,情景一生態(tài)補償量置信區(qū)間為12.51—17.42億元,平均值為14.95億元(圖4)；情景二生態(tài)補償量置信區(qū)間為9.34—21.93億元,平均值為15.53億元(圖4)。生態(tài)補償量預(yù)測值的置信區(qū)間包含了上述計算結(jié)果,代表了4個模型參數(shù)不確定條件下的所有可能結(jié)果。情景一下的生態(tài)補償量大于15.00億元的確定性內(nèi)概率為45.95%；情景二下的生態(tài)補償量大于15.00億元的確定性內(nèi)概率為52.66%(圖5)。說明情景一下的流域水污染治理在大部分情況下較情景二的方案更為經(jīng)濟(jì),與本文上述生態(tài)補償量具體化計算結(jié)果一致。

圖4 蒙特卡羅模擬生態(tài)補償量直方圖Fig.4 The histograms of ecological compensation by Monte Carlo simulation

圖5 兩種情景下生態(tài)補償成本對比圖Fig.5 Contrast figure of ecological compensation with two scenarios

情景一下引起生態(tài)補償量變化的靈敏參數(shù)為高架床養(yǎng)殖容積比(k1)、限養(yǎng)區(qū)污水入河系數(shù)(k2)、高架床養(yǎng)殖模式日排水量(k3),方差貢獻(xiàn)率分別為79.50%、18.90%和1.60%；情景二下引起生態(tài)補償量變化的靈敏參數(shù)為限養(yǎng)區(qū)污水入河系數(shù)(k2)、傳統(tǒng)養(yǎng)殖模式日排水量(k4),方差貢獻(xiàn)率依次為88.80%、11.20%(圖6)?？傮w來看,情景一下的高架床養(yǎng)殖容積比(k1)參數(shù)對流域生態(tài)補償量最敏感,其余參數(shù)影響較少。情景二模型中涉及兩個參數(shù),其限養(yǎng)區(qū)污水入河系數(shù)(k2)為主要影響生態(tài)補償量的參數(shù)因子。

圖6 參數(shù)敏感性分析Fig.6 The Monte Carlo sensitivity analysis for parameters

綜上所述,不同流域水污染治理方案的選擇下,其需要的生態(tài)補償成本存在差異,影響其生態(tài)補償成本的參數(shù)貢獻(xiàn)率也各有不同。本文通過模擬計算不確定條件下的生態(tài)補償量的所有可能值以及進(jìn)行模型參數(shù)靈敏度分析,為決策者提供了更全面的信息。

4 結(jié)論與展望

4.1 結(jié)論

流域生態(tài)補償是一個跨行政區(qū)域的綜合生態(tài)環(huán)境問題。生態(tài)補償量的計算和測定是流域生態(tài)補償?shù)那疤?也是決定能否順利實施補償?shù)年P(guān)鍵環(huán)節(jié)。針對九洲江流域水生態(tài)補償政策的問題與不足,本文嘗試構(gòu)建以水質(zhì)達(dá)標(biāo)為前提的、以最小投入為目標(biāo)的、兼顧責(zé)任與公平的跨區(qū)域流域上下游政府水環(huán)境資源生態(tài)補償量計算模型,主要結(jié)論如下：

在Ⅲ類水質(zhì)目標(biāo)要求下,計算得出水環(huán)境治理生態(tài)補償資金最小投入為15.03億元。其禁養(yǎng)區(qū)內(nèi)生豬養(yǎng)殖拆遷補償為3.35億元,限養(yǎng)區(qū)內(nèi)高架床生態(tài)養(yǎng)殖模式改造補償為11.68億元。中央及廣西廣東地方政府生態(tài)補償分?jǐn)偨痤~分別為3.76、3.30億元及7.96億元。此外,模型參數(shù)的蒙特卡羅模擬和靈敏度分析結(jié)果顯示,95%置信水平下,情景一生態(tài)補償量置信區(qū)間為12.51—17.42億元,平均值為14.95億元；情景二生態(tài)補償量置信區(qū)間為9.34—21.93億元,平均值為15.53億元。情景一模型中高架床養(yǎng)殖容積比(k1)、限養(yǎng)區(qū)污水入河系數(shù)(k2)、高架床養(yǎng)殖模式日排水量(k3)等參數(shù)對生態(tài)補償量具有一定的影響,方差貢獻(xiàn)比分別為79.50%、18.90%和1.60%；情景二模型中限養(yǎng)區(qū)污水入河系數(shù)(k2)參數(shù)對生態(tài)補償量最敏感,方差貢獻(xiàn)率達(dá)到88.80%,其傳統(tǒng)養(yǎng)殖模式日排水量(k4)參數(shù)對生態(tài)補償量影響較小,方差貢獻(xiàn)率為11.20%?？梢婌`敏度分析結(jié)果代表了模型中參數(shù)不確定條件下生態(tài)補償量的所有可能值,為決策者提供了更全面的信息。

4.2 不足與展望

通過上文研究分析可以看出,該文以保持水質(zhì)投入的成本為補償金額,構(gòu)建了基于水質(zhì)水量的生態(tài)補償量測算模型,使得補償量得到了合理量化,但如何形成一整套完備的評價理論、指標(biāo)體系、實施原則與計算方法等許多科學(xué)問題依然存在模糊和不確定性,生態(tài)補償成本構(gòu)成方面仍然有待更深入的探討和分析,以避免統(tǒng)計的不完全。此外,養(yǎng)殖規(guī)模,養(yǎng)殖空間差異等也是生態(tài)補償研究中重要內(nèi)容,生態(tài)補償標(biāo)準(zhǔn)在不同地區(qū),不同養(yǎng)殖規(guī)模與方式存在差異性,與地區(qū)養(yǎng)殖特點有一定的相關(guān)性,在今后補償標(biāo)準(zhǔn)制定中還應(yīng)考慮到此異質(zhì)性。

對于未來還可開展的研究工作,一方面,本文未考慮生態(tài)補償方式的可行性,本文生態(tài)補償模式是資金的縱向轉(zhuǎn)移和橫向支付,由于存在諸多的制度限制與障礙,現(xiàn)實中單一的資金補償模式無法滿足實際需要,也不容易操作,因此,還可考慮到補償方式的不同,采取多種方式優(yōu)勢互補,將產(chǎn)業(yè)扶持、生態(tài)旅游、綠色農(nóng)業(yè)等新式產(chǎn)業(yè)納入到生態(tài)補償機(jī)制中。另一方面,畜禽養(yǎng)殖過程中會使用到一定比例的抗生素藥劑,畜禽糞便污水中會攜帶一部分抗生素殘留,未來的生態(tài)補償研究還應(yīng)適當(dāng)?shù)目紤]到抗生素等難降解污染物的控制。因此該研究的生態(tài)補償模型根據(jù)實際情況還需要進(jìn)一步的修正。建立適合我國流域的生態(tài)補償還需各利益相關(guān)方的共同參與,各種制度的完善,還有待進(jìn)一步的研究。