大汶河流域生態(tài)補(bǔ)償經(jīng)濟(jì)環(huán)境的效應(yīng)分析

摘 要：本文基于2004—2023年大汶河流域內(nèi)各地市的面板數(shù)據(jù)，以該流域生態(tài)補(bǔ)償實(shí)際情況為研究對(duì)象，運(yùn)用雙重差分模型探究流域生態(tài)補(bǔ)償機(jī)制的環(huán)境、經(jīng)濟(jì)效應(yīng)。研究發(fā)現(xiàn)，大汶河生態(tài)補(bǔ)償對(duì)該流域的環(huán)境狀況改善和整體經(jīng)濟(jì)發(fā)展都存在正面影響，但其影響強(qiáng)度隨著政策實(shí)施周期推進(jìn)逐漸降低，應(yīng)當(dāng)延續(xù)并優(yōu)化此類流域生態(tài)補(bǔ)償政策，推動(dòng)大汶河流域可持續(xù)發(fā)展。

關(guān)鍵詞：流域生態(tài)補(bǔ)償；經(jīng)濟(jì)效應(yīng)；環(huán)境效益；雙重差分模型

中圖分類號(hào)：TP 79 文獻(xiàn)標(biāo)志碼：A

隨著經(jīng)濟(jì)發(fā)展與生態(tài)保護(hù)之間的矛盾凸顯，以促進(jìn)流域高質(zhì)量發(fā)展、上下游流域協(xié)同發(fā)展的生態(tài)補(bǔ)償制度逐漸清晰。黨中央和國(guó)務(wù)院對(duì)流域生態(tài)保護(hù)的關(guān)注度日益增加，并將健全我國(guó)生態(tài)保護(hù)補(bǔ)償機(jī)制作為重要任務(wù)。狹義上，流域生態(tài)補(bǔ)償就是由下游生態(tài)受益區(qū)居民對(duì)流域內(nèi)上游經(jīng)濟(jì)主體的生態(tài)保護(hù)行為進(jìn)行激勵(lì)和補(bǔ)貼，以達(dá)到生態(tài)保護(hù)的目的，我國(guó)的相關(guān)學(xué)術(shù)探索工作起源于20世紀(jì)80年代。從生態(tài)補(bǔ)償?shù)膬?nèi)涵來(lái)說(shuō)，門寶輝等[1]認(rèn)為流域生態(tài)補(bǔ)償是指對(duì)因人類活動(dòng)而產(chǎn)生的流域生態(tài)環(huán)境負(fù)外部效應(yīng)進(jìn)行補(bǔ)償?shù)男袨?。牛坤玉等[2]提出生態(tài)補(bǔ)償也就是生態(tài)系統(tǒng)付費(fèi)服務(wù)。是人類從自然界獲得的福利。從生態(tài)補(bǔ)償標(biāo)準(zhǔn)來(lái)看，羅萬(wàn)云等[3]提出在生態(tài)公平視域下理清權(quán)責(zé)關(guān)系，確定流域補(bǔ)償標(biāo)準(zhǔn)。鄧夢(mèng)華等[4]基于學(xué)習(xí)模型XGBoost-SHAP，將流域居民補(bǔ)償意愿納入其中。綜上，學(xué)者們對(duì)流域生態(tài)補(bǔ)償進(jìn)行規(guī)范化研究，但該研究還有待進(jìn)一步深入。本文以大汶河流域的生態(tài)補(bǔ)償機(jī)制為研究對(duì)象，為流域協(xié)調(diào)發(fā)展提供依據(jù)。

1 模型構(gòu)建與變量選取

本文采用雙重差分法（DID）對(duì)大汶河流域生態(tài)補(bǔ)償政策進(jìn)行評(píng)估。雙重差分法（DID）雙重差分法最早由ASHENFELTER引入經(jīng)濟(jì)學(xué)研究中，是基于各地市面板數(shù)據(jù)進(jìn)行機(jī)制評(píng)估的方法，本文建立的模型如公式（1）所示。

Yit=β0+β1didit+β2x+δi+θt+εit （1）

式中：Yit為表征環(huán)境效應(yīng)、經(jīng)濟(jì)效應(yīng)的被解釋變量；didit為核心解釋變量；X為表征環(huán)境效應(yīng)、經(jīng)濟(jì)效應(yīng)的控制變量組；δi為個(gè)體效應(yīng)；θt為時(shí)間效應(yīng)；εit為隨機(jī)干擾項(xiàng)；β1為所求政策效應(yīng)值，如果系數(shù)符號(hào)為正，就說(shuō)明生態(tài)補(bǔ)償機(jī)制將改善環(huán)境，促進(jìn)經(jīng)濟(jì)發(fā)展，如果系數(shù)符號(hào)為負(fù)，就說(shuō)明生態(tài)補(bǔ)償機(jī)制將不利于環(huán)境改善和經(jīng)濟(jì)發(fā)展；β2為控制變量組系數(shù)。

1.1 被解釋變量

廢水排放量：減少流域內(nèi)污染物排放，促進(jìn)生態(tài)系統(tǒng)良性發(fā)展是流域生態(tài)補(bǔ)償?shù)闹匾康?。本文通過(guò)大汶河流域各地市的年工業(yè)廢水排放總量來(lái)衡量[4]。

地區(qū)生產(chǎn)總值：GDP是國(guó)民經(jīng)濟(jì)核算的核心指標(biāo)，本文通過(guò)引用大汶河流域地市生產(chǎn)總值來(lái)反映流域內(nèi)生態(tài)補(bǔ)償政策實(shí)施后帶來(lái)的經(jīng)濟(jì)效果影響。

1.2 核心解釋變量

本文在模型中引入地區(qū)與時(shí)間虛擬變量的交互項(xiàng)，來(lái)檢驗(yàn)試驗(yàn)組與對(duì)照組由大汶河生態(tài)補(bǔ)償機(jī)制而引起的流域內(nèi)經(jīng)濟(jì)和生態(tài)環(huán)境的變化情況。地區(qū)虛擬變量（treat）試驗(yàn)組取1，對(duì)照組取0；時(shí)間虛擬變量將大汶河流域機(jī)制實(shí)施的時(shí)間（t≥2008）取1，機(jī)制實(shí)施之前的年份（t＜2008）取0。

1.3 控制變量

產(chǎn)業(yè)結(jié)構(gòu)：區(qū)域產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整是區(qū)域內(nèi)經(jīng)濟(jì)發(fā)展水平和生態(tài)環(huán)境變化的重要體現(xiàn)指標(biāo)之一。城市現(xiàn)代化的發(fā)展過(guò)程當(dāng)中，第二、三產(chǎn)業(yè)的發(fā)展情況對(duì)經(jīng)濟(jì)進(jìn)步和環(huán)境變化起到重要影響，因此本文取流域內(nèi)各地市的第三產(chǎn)業(yè)和第二產(chǎn)業(yè)增加值比重的對(duì)數(shù)來(lái)表示。

科技支出：創(chuàng)新會(huì)帶來(lái)科技水平的進(jìn)步，是發(fā)展的源泉和第一動(dòng)力，而科技進(jìn)步對(duì)轉(zhuǎn)換城市經(jīng)濟(jì)增長(zhǎng)動(dòng)力、提高各產(chǎn)業(yè)水資源利用率和提升環(huán)境污染治理能力等方面均能起到積極作用，促進(jìn)生態(tài)環(huán)境與經(jīng)濟(jì)建設(shè)可持續(xù)發(fā)展。因此，本文以各地區(qū)年政府科技支出的對(duì)數(shù)來(lái)表示。

城鎮(zhèn)化率：改革開放以來(lái)，我國(guó)的城鎮(zhèn)化率大大提高，包括大汶河流域在內(nèi)，為滿足日益增長(zhǎng)的人口對(duì)生態(tài)資源的需要，人類社會(huì)的發(fā)展在一定程度上對(duì)該流域生態(tài)環(huán)境施加了更大的壓力。因此本文以各地市的城鎮(zhèn)常住人口占地區(qū)常住人口的比例來(lái)表示。

教育支出：教育支出是影響該地區(qū)市民的受教育水平的重要指標(biāo)，受教育程度越高，越能夠?qū)斫馍鷳B(tài)保護(hù)的重要性并積極配合政府相關(guān)的環(huán)境政策和環(huán)保號(hào)召。因此本文取教育和文化財(cái)政支出的對(duì)數(shù)來(lái)表示。

生態(tài)支出：生態(tài)支出是政府購(gòu)買的重要組成部分，主要用于政府生態(tài)修復(fù)工程的投資、相關(guān)生態(tài)產(chǎn)品與服務(wù)的購(gòu)買支出、有關(guān)生態(tài)金融衍生品和政策的補(bǔ)貼與補(bǔ)助等方面。因此本文取各地區(qū)歷年的農(nóng)業(yè)、林業(yè)和水資源這三類環(huán)境保護(hù)支出總和的對(duì)數(shù)來(lái)表示。

2 大汶河流域生態(tài)補(bǔ)償政策效應(yīng)評(píng)估

2.1 描述性統(tǒng)計(jì)

表1為大汶河流域生態(tài)補(bǔ)償機(jī)制的環(huán)境效應(yīng)和經(jīng)濟(jì)效應(yīng)的相關(guān)變量的描述性統(tǒng)計(jì)。樣本總數(shù)為270，其中2004—2023年該地區(qū)工業(yè)廢水排放量平均數(shù)為7543萬(wàn)t，標(biāo)準(zhǔn)差為5755萬(wàn)t，最小值1210萬(wàn)t，中位數(shù)為5772萬(wàn)t，最大值為21212萬(wàn)t，說(shuō)明該流域工業(yè)廢水排放量分布不均衡，也表明少部分地區(qū)工業(yè)較發(fā)達(dá)，大部分地區(qū)工業(yè)化水平較落后。2004—2023年該地區(qū)生產(chǎn)總值平均數(shù)為2828億元，標(biāo)準(zhǔn)差為3259億元，最小為值223.9億元，中位數(shù)為2233億元，最大值為22901億元，說(shuō)明經(jīng)濟(jì)水平發(fā)展不均衡，生產(chǎn)總值最小地區(qū)只占生產(chǎn)總值最大地區(qū)的10%左右。2004—2023年該地區(qū)產(chǎn)業(yè)結(jié)構(gòu)中第三產(chǎn)業(yè)和第二產(chǎn)業(yè)增加值的比重平均數(shù)為0.810，標(biāo)準(zhǔn)差為0.330，最小值為0.380，最大值為1.680，說(shuō)明該流域內(nèi)不同地市在產(chǎn)業(yè)發(fā)展方向有所側(cè)重，并未全部發(fā)展工業(yè)或服務(wù)業(yè)，產(chǎn)業(yè)發(fā)展結(jié)構(gòu)相對(duì)合理。

2.2 平行趨勢(shì)檢驗(yàn)

雙重差分法的應(yīng)用前提是試驗(yàn)組與對(duì)照組在機(jī)制的實(shí)施前具有相同的發(fā)展趨勢(shì)[5]。因此，在機(jī)制效應(yīng)評(píng)估前，本文構(gòu)建模型，如公式（2）所示。

（2）

式中：Yit為表征環(huán)境效應(yīng)、經(jīng)濟(jì)效應(yīng)的被解釋變量，分別代入廢水排放量和地區(qū)生產(chǎn)總值進(jìn)行檢驗(yàn)；β0為生態(tài)補(bǔ)償系數(shù)；δi為個(gè)體效應(yīng)；θt為時(shí)間效應(yīng)；εit為隨機(jī)干擾項(xiàng)；縱坐標(biāo)為流域污水排放量和地區(qū)生產(chǎn)總值回歸系數(shù)；prej和lastl分別為政策實(shí)施前、后的時(shí)間虛擬變量。

如果大汶河流域生態(tài)補(bǔ)償機(jī)制的實(shí)施年份為2008年，那么pre2代表政策實(shí)施的兩年前，即2006年；同理，post2為政策實(shí)施后第二年，即2010年。結(jié)果如圖1、圖2所示，在政策執(zhí)行前，對(duì)生態(tài)補(bǔ)償機(jī)制系數(shù)β的置信區(qū)間均為零，這說(shuō)明在政策執(zhí)行前，資料處理試驗(yàn)組與對(duì)照組的廢水排放量和地區(qū)生產(chǎn)總值并未出現(xiàn)明顯區(qū)別，因此可以將其視為平行。在2008年流域相關(guān)保護(hù)政策正式執(zhí)行后，政策系數(shù)β的置信區(qū)間在往后連續(xù)的四期內(nèi)都不為零，說(shuō)明政府制定并執(zhí)行的生態(tài)政策已經(jīng)對(duì)現(xiàn)實(shí)情況中的處理廢水排放量和地區(qū)生產(chǎn)總值等相關(guān)研究要素形成了一定限制與影響，該結(jié)果滿足平行趨勢(shì)假設(shè)成立所需要求，因此可以運(yùn)用雙重差分法展開進(jìn)一步的分析研究。

2.3 雙重差分檢驗(yàn)結(jié)果

本文借助Stata17軟件進(jìn)行基本回歸分析。雙重查分檢驗(yàn)結(jié)果見表2。其中（1）、（3）列為廢水排放量環(huán)境效應(yīng)評(píng)估模型，（2）、（4）列為地區(qū)生產(chǎn)總值為經(jīng)濟(jì)效應(yīng)評(píng)估模型。并且（1）列為對(duì)照組，（3）列為試驗(yàn)組，試驗(yàn)組相較對(duì)照組增加了上文所描述與統(tǒng)計(jì)的若干解釋變量。從該模型檢驗(yàn)的結(jié)果可知，生態(tài)補(bǔ)償機(jī)制的環(huán)境效應(yīng)系數(shù)為-4.273，在1%的水平上顯著為負(fù)，表明其在減少該流域工業(yè)廢水排放量上起到積極作用，即減少了該流域工業(yè)廢水排放總量，顯著改善了流域水資源的生態(tài)環(huán)境。同時(shí)由表2可得，生態(tài)補(bǔ)償機(jī)制經(jīng)濟(jì)效應(yīng)系數(shù)為4.323，在1％的水平上顯著為正，表明生態(tài)補(bǔ)償政策對(duì)流域的經(jīng)濟(jì)發(fā)展起到了明顯作用，即制定實(shí)施相關(guān)政策顯推動(dòng)了GDP等經(jīng)濟(jì)指標(biāo)增長(zhǎng)，提高了該區(qū)域的整體經(jīng)濟(jì)發(fā)展水平。

2.4 穩(wěn)健性檢驗(yàn)

為進(jìn)一步驗(yàn)證結(jié)果的準(zhǔn)確性，本文采用滯后一期實(shí)施有關(guān)流域生態(tài)政策的檢驗(yàn)方法對(duì)上述結(jié)果進(jìn)行穩(wěn)健性檢驗(yàn)?；貧w結(jié)果見表3，其中（1）列代表廢水排放量，（2）列代表生產(chǎn)總值。從表3中可得，就業(yè)率在1%的水平上對(duì)廢水排放量指標(biāo)顯著為負(fù)，對(duì)地區(qū)生產(chǎn)總值指標(biāo)顯著為正，說(shuō)明該政策的實(shí)施通過(guò)在生態(tài)保護(hù)行業(yè)吸納就業(yè)人口，達(dá)到了明顯降低了污水排放量和提高經(jīng)濟(jì)社會(huì)發(fā)展的目的，同時(shí)該流域生態(tài)政策的環(huán)境機(jī)制效應(yīng)系數(shù)為-2.403，仍在1%的水平上顯著為負(fù)，即延后實(shí)施一期比前一期工業(yè)污水排放量明顯下降，經(jīng)濟(jì)機(jī)制效應(yīng)系數(shù)為2.126，在5%的水平上顯著為正，即延后實(shí)施一期比前一期經(jīng)濟(jì)發(fā)展水平明顯上升，與本文雙重差分法的檢驗(yàn)結(jié)果保持一致。

3 結(jié)語(yǔ)

本文經(jīng)過(guò)相關(guān)研究結(jié)論如下：大汶河流域生態(tài)補(bǔ)償政策實(shí)施后起到了改善當(dāng)?shù)厣鷳B(tài)環(huán)境的作用，并且?guī)?dòng)了當(dāng)?shù)亟?jīng)濟(jì)的發(fā)展；但是大汶河流域生態(tài)補(bǔ)償機(jī)制的環(huán)境效應(yīng)隨時(shí)間發(fā)展其相關(guān)性逐漸減弱，最明顯的就是政策實(shí)施后的第三年效果最顯著。

大汶河流域的生態(tài)補(bǔ)償政策本身和其他方面仍存在不足。1）大汶河流域生態(tài)補(bǔ)償?shù)膮⑴c主體基本是政府，這在一定程度上帶來(lái)了轉(zhuǎn)移支付的財(cái)政壓力和政府部門為實(shí)現(xiàn)政策目標(biāo)而進(jìn)行人員配置等環(huán)節(jié)的執(zhí)行成本，在后疫情時(shí)代對(duì)財(cái)政相對(duì)緊張的各級(jí)政府來(lái)說(shuō)是很難繼續(xù)保證政策實(shí)施成效的。2）大汶河流域的經(jīng)濟(jì)效益提升更多來(lái)自社會(huì)生產(chǎn)力的發(fā)展和產(chǎn)業(yè)升級(jí)，對(duì)生態(tài)環(huán)境本身經(jīng)濟(jì)價(jià)值的探索還停留在初級(jí)階段，沒(méi)有很好地發(fā)揮大汶河流域生態(tài)的經(jīng)濟(jì)效應(yīng)。3）大汶河生態(tài)補(bǔ)償機(jī)制的手段不夠豐富，基本依靠較直接的政府轉(zhuǎn)移支付和政府投資來(lái)推動(dòng)生態(tài)修復(fù)和拉動(dòng)經(jīng)濟(jì)發(fā)展，在更復(fù)雜的產(chǎn)權(quán)交易方面缺少探索和實(shí)踐，很難形成由市場(chǎng)為主導(dǎo)驅(qū)動(dòng)的生態(tài)產(chǎn)業(yè)鏈。

參考文獻(xiàn)

[1]門寶輝，劉菁蘋.基于最小數(shù)據(jù)法的黃河中游延安市生態(tài)補(bǔ)償標(biāo)準(zhǔn)研究[J].華北水利水電大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，43（5）：101-108.

[2]牛坤玉，劉靜，郭靜利，等.農(nóng)業(yè)生態(tài)補(bǔ)償：內(nèi)涵、要素特征與政策創(chuàng)設(shè)[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃，2024（2）：1-10.

[3]羅萬(wàn)云，周楊，王小娟.生態(tài)公平視域下額爾齊斯河流域生態(tài)補(bǔ)償標(biāo)準(zhǔn)及空間選擇研究[J].生態(tài)學(xué)報(bào)，2024，（21）：1-16.

[4]鄧夢(mèng)華，張?zhí)焓妫愜婏w.基于XGBoost-SHAP模型的太湖流域居民生態(tài)補(bǔ)償支付意愿影響因素研究[J].水利經(jīng)濟(jì)，2024，42（2）：44-50.

[5]于冰，史穎.新安江生態(tài)補(bǔ)償環(huán)境經(jīng)濟(jì)效應(yīng)及影響機(jī)制研究[J].生態(tài)經(jīng)濟(jì)，2023，39（7）：157-164.