董 冬，羅 毅，王麗宸，邰 宇，顧康康①

(1.安徽建筑大學(xué)建筑與規(guī)劃學(xué)院，安徽 合肥 230022；2.安徽省國土空間規(guī)劃與生態(tài)研究院，安徽 合肥 230022)

新安江西起安徽省黃山市休寧縣，東至浙江省杭州市，是長江流域重要的水源地和生態(tài)屏障[1]。新安江水質(zhì)關(guān)乎流域全局，具有重要意義。20世紀(jì)末，新安江上游城鎮(zhèn)化和工業(yè)化進(jìn)程加快，加重了流域資源環(huán)境壓力，也嚴(yán)重影響區(qū)域人居環(huán)境質(zhì)量。2012年，新安江流域生態(tài)補(bǔ)償機(jī)制正式實(shí)施，這是習(xí)近平總書記親自倡導(dǎo)和推動的全國首個跨省流域生態(tài)補(bǔ)償試點(diǎn)。經(jīng)過3輪試點(diǎn)[2]，新安江-千島湖生態(tài)補(bǔ)償試驗(yàn)區(qū)(簡稱試驗(yàn)區(qū))于2020年正式開始建設(shè)。作為命運(yùn)共同體，試驗(yàn)區(qū)通過禁止和限制上游重點(diǎn)生態(tài)功能區(qū)的開發(fā)，達(dá)到保障并改善下游生態(tài)安全的目的；下游地區(qū)對上游地區(qū)開展生態(tài)環(huán)境保護(hù)與修復(fù)、放棄發(fā)展機(jī)會和治理農(nóng)村面源污染等方面進(jìn)行價值補(bǔ)償[3]。生態(tài)補(bǔ)償作為一項(xiàng)重要的經(jīng)濟(jì)、社會和環(huán)境政策，已經(jīng)拓展到生態(tài)環(huán)境保護(hù)與治理、經(jīng)濟(jì)發(fā)展、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整、村民就業(yè)與生活水平提高等領(lǐng)域。這些措施改變了鄉(xiāng)村地區(qū)生態(tài)環(huán)境和生產(chǎn)、生活方式，進(jìn)而對鄉(xiāng)村人居環(huán)境系統(tǒng)產(chǎn)生影響。試驗(yàn)區(qū)是皖浙兩省生態(tài)環(huán)境同治的“生態(tài)文明創(chuàng)建共同體”，也是協(xié)同推進(jìn)鄉(xiāng)村人居環(huán)境整治與發(fā)展的重要空間單元。筆者研究對該區(qū)域鄉(xiāng)村人居環(huán)境質(zhì)量進(jìn)行綜合評價，探索關(guān)鍵制約因素并構(gòu)建反饋機(jī)制，這對優(yōu)化區(qū)域鄉(xiāng)村人居環(huán)境質(zhì)量和生態(tài)補(bǔ)償政策具有重要現(xiàn)實(shí)意義。

鄉(xiāng)村人居環(huán)境指鄉(xiāng)村地域內(nèi)村民進(jìn)行生產(chǎn)生活、休憩游樂和社會交往所需的物質(zhì)和非物質(zhì)環(huán)境[4]。人居環(huán)境質(zhì)量是衡量村民生活品質(zhì)和鄉(xiāng)村區(qū)域進(jìn)步的重要載體[5]，鄉(xiāng)村人居環(huán)境評價一直是國內(nèi)外學(xué)者們關(guān)注的焦點(diǎn)。從研究內(nèi)容來看，國外研究主要集中在鄉(xiāng)村聚落形態(tài)演替[6-7]、鄉(xiāng)村人居環(huán)境質(zhì)量[6,8]和鄉(xiāng)村可持續(xù)發(fā)展等方面[9-11]。國內(nèi)人居環(huán)境研究起源于吳良鏞院士提出“人居環(huán)境科學(xué)”理論[12]，主要聚焦于相關(guān)理論體系研究[13-15]、鄉(xiāng)村空間聚落分異[16-18]、鄉(xiāng)村人居環(huán)境演化[19-21]和鄉(xiāng)村環(huán)境整治與優(yōu)化[22-23]等方面。其中，鄉(xiāng)村人居環(huán)境評價研究多為對復(fù)雜鄉(xiāng)村人居環(huán)境從單一角度進(jìn)行綜合評價[24-25]，鮮有從時間和空間多角度開展研究。從研究尺度來看，現(xiàn)有鄉(xiāng)村人居環(huán)境質(zhì)量評價研究尺度多聚焦于生態(tài)敏感區(qū)、省域、市域和縣域尺度[26-28]，科學(xué)評估跨區(qū)域生態(tài)補(bǔ)償試驗(yàn)區(qū)人居環(huán)境的研究鮮見報道。從研究方法來看，多采用層次分析法、模糊綜合層次法、熵值法和結(jié)構(gòu)方程模型等方法進(jìn)行綜合評價[29-31]，且評價指標(biāo)體系特色不夠突出[32]?；诖?，以試驗(yàn)區(qū)為研究區(qū)域，選取生態(tài)補(bǔ)償機(jī)制確立之前的2011年以及3輪試點(diǎn)收官之年(2014、2017和2019年)的相關(guān)研究數(shù)據(jù)，構(gòu)建具有地域特征的評價指標(biāo)體系，采用熵權(quán)TOPSIS法和時間地理加權(quán)回歸模型(GTWR)相結(jié)合的方法評價鄉(xiāng)村人居環(huán)境質(zhì)量，并理清其時間演變規(guī)律、時空分異特征及影響機(jī)制，以期為區(qū)域鄉(xiāng)村人居環(huán)境質(zhì)量的優(yōu)化提升和生態(tài)補(bǔ)償政策的制定提供建議。

1 材料與方法

1.1 研究區(qū)概況與數(shù)據(jù)來源

新安江干流全長約為359 km。2012年以來，新安江流域生態(tài)保護(hù)補(bǔ)償機(jī)制啟動實(shí)施并經(jīng)過3輪試點(diǎn)(2012—2014、2015—2017和2018—2020年)。試驗(yàn)區(qū)包括安徽省黃山市全境和宣城市績溪縣以及浙江省杭州市淳安縣和建德市(圖1)，共計10個縣(區(qū)、市)。該區(qū)域?qū)儆趤啛釒Ъ撅L(fēng)性氣候區(qū)，四季分明，雨水充沛。全區(qū)位于長江中下游平原，地形復(fù)雜多樣，以丘陵山地為主，森林覆蓋率較高，地表水系發(fā)達(dá)，河網(wǎng)湖泊密布。各縣(區(qū)、市)在生態(tài)環(huán)境現(xiàn)狀、人口分布情況、經(jīng)濟(jì)水平和基礎(chǔ)設(shè)施完善程度等方面差距較大。

圖1 新安江流域地理范圍與區(qū)位

數(shù)據(jù)資料主要包括2011、2014、2017和2019年《中國縣域統(tǒng)計年鑒》《黃山統(tǒng)計年鑒》《杭州統(tǒng)計年鑒》《宣城統(tǒng)計年鑒》以及《新安江流域上下游橫向生態(tài)補(bǔ)償試點(diǎn)績效評估報告(2012—2017年)》，各市統(tǒng)計局、環(huán)保局、住建局和水利局等部門提供的相關(guān)數(shù)據(jù)以及相關(guān)規(guī)劃、衛(wèi)星影像數(shù)據(jù)和實(shí)地調(diào)研資料。

1.2 研究方法

1.2.1指標(biāo)體系的構(gòu)建

鄉(xiāng)村人居環(huán)境是一項(xiàng)涉及經(jīng)濟(jì)、社會、資源和生態(tài)環(huán)境等方面的復(fù)雜系統(tǒng)，且內(nèi)部存在著物質(zhì)、能量和信息等要素的流動與交互，與鄉(xiāng)村空間尺度上的生態(tài)空間、生產(chǎn)空間和生活空間的合理利用相契合[33]。生態(tài)補(bǔ)償機(jī)制的實(shí)施使鄉(xiāng)村人居環(huán)境賴以依存的社會、經(jīng)濟(jì)和自然環(huán)境發(fā)生較大變化。生態(tài)補(bǔ)償區(qū)政府對現(xiàn)實(shí)環(huán)境作出響應(yīng)，通過政策調(diào)控以及實(shí)施“農(nóng)藥替代花費(fèi)工程”“綠色特色農(nóng)業(yè)發(fā)展工程”“農(nóng)村環(huán)境整治工程”“畜禽規(guī)模養(yǎng)殖工程”等“十大工程”，直接或間接改善鄉(xiāng)村人居環(huán)境[34]。構(gòu)建生態(tài)補(bǔ)償試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境體系時充分體現(xiàn)生態(tài)補(bǔ)償機(jī)制的影響以及生態(tài)功能區(qū)人居環(huán)境的復(fù)雜性和主體功能。在指標(biāo)選擇上，參照以往研究成果[32,35]并遵循數(shù)據(jù)科學(xué)性、對比性和可獲得性等原則，著重體現(xiàn)生態(tài)補(bǔ)償機(jī)制對生態(tài)空間、生活空間和生產(chǎn)空間質(zhì)量要素的影響，在初擬評價指標(biāo)基礎(chǔ)上，結(jié)合10位長期從事村莊規(guī)劃、鄉(xiāng)村振興研究專家意見，構(gòu)建了包含鄉(xiāng)村生態(tài)環(huán)境(B1)、鄉(xiāng)村生產(chǎn)環(huán)境(B2)和鄉(xiāng)村生活環(huán)境(B3)3個1級指標(biāo)和27個2級指標(biāo)的試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量綜合評價指標(biāo)體系(表1)。

表1 生態(tài)補(bǔ)償試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量評價指標(biāo)體系

1.2.2熵權(quán)TOPSIS模型

(1)

(2)

式(1)～(2)中，rj為指標(biāo)j原值；n為指標(biāo)數(shù)。

各縣(區(qū)、市)鄉(xiāng)村人居環(huán)境與正負(fù)理想解的貼近度(Ci)計算公式為

Ci取值范圍為0～1，Ci越接近1，說明鄉(xiāng)村人居環(huán)境質(zhì)量越好。

1.2.3時空地理加權(quán)回歸(GTWR)模型

采用GTWR模型與熵權(quán)TOPSIS法計算的試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量及其子系統(tǒng)得分進(jìn)行回歸分析，測度鄉(xiāng)村人居環(huán)境質(zhì)量及其影響因素的時空分異規(guī)律。GTWR模型是以地理加權(quán)回歸模型為基礎(chǔ)，將時間維度數(shù)據(jù)納入其中構(gòu)成的時空三維分析模型[37]，計算公式為

(3)

式(3)中，yi為鄉(xiāng)村人居環(huán)境質(zhì)量綜合評價得分，i=1，2，…，n；ui、vi和ti分別為點(diǎn)i的經(jīng)度、緯度和時間；β0(ui，vi，ti)為截距項(xiàng)；βj(ui，vi，ti)為影響因素j在點(diǎn)i上的回歸系數(shù)；xij為影響因素j在點(diǎn)i上的值；εi為隨機(jī)誤差。

2 結(jié)果與分析

2.1 鄉(xiāng)村人居環(huán)境質(zhì)量時間演變特征分析

根據(jù)人居環(huán)境質(zhì)量評價指標(biāo)體系和權(quán)重，采用熵權(quán)TOPSIS法和GTWR模型進(jìn)行綜合評價計算，繪制2011、2014、2017和2019年試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量以及鄉(xiāng)村生態(tài)環(huán)境、鄉(xiāng)村生產(chǎn)環(huán)境和鄉(xiāng)村生活環(huán)境3個子系統(tǒng)平均值和標(biāo)準(zhǔn)偏差的變化趨勢(圖2)。圖2顯示，隨著新安江流域生態(tài)補(bǔ)償機(jī)制制度的演進(jìn)，鄉(xiāng)村人居環(huán)境綜合質(zhì)量呈不斷上升態(tài)勢，但各子系統(tǒng)得分變化趨勢差異較大。鄉(xiāng)村生態(tài)和鄉(xiāng)村生活環(huán)境得分由生態(tài)補(bǔ)償機(jī)制實(shí)施前的0.485和0.478分別上升到0.573和0.499，而鄉(xiāng)村生長環(huán)境呈現(xiàn)下降趨勢。各縣(區(qū)、市)鄉(xiāng)村人居環(huán)境綜合質(zhì)量和生態(tài)環(huán)境質(zhì)量的標(biāo)準(zhǔn)偏差不斷增加。通過3輪新安江流域生態(tài)補(bǔ)償試點(diǎn)，全面推進(jìn)了新安江流域各縣(區(qū)、市)經(jīng)濟(jì)、社會、產(chǎn)業(yè)的發(fā)展。雖然試驗(yàn)區(qū)各縣(區(qū)、市)之間鄉(xiāng)村人居環(huán)境質(zhì)量差距在不斷擴(kuò)大，但試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量總體上不斷提升。

圖2 2011—2019年新安江-千島湖生態(tài)補(bǔ)償試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量變化

新安江生態(tài)補(bǔ)償機(jī)制通過3輪9年試點(diǎn)，各縣(區(qū)、市)鄉(xiāng)村人居環(huán)境質(zhì)量總體上趨于向好，但在不同試點(diǎn)階段間存在一定波動(圖3)，2011—2019年各縣(區(qū)、市)鄉(xiāng)村人居環(huán)境質(zhì)量增幅在前5位的分別為祁門縣(12.69%)、建德市(12.67%)、徽州區(qū)(7.91%)、休寧縣(7.48%)和淳安縣(7.14%)，這些地區(qū)均位于新安江-千島湖源頭和腹地。

圖3 各縣(區(qū)、市)鄉(xiāng)村人居環(huán)境綜合質(zhì)量比較

第1輪試點(diǎn)(2012—2014年)期間，除屯溪區(qū)鄉(xiāng)村人居環(huán)境綜合質(zhì)量有所下降外，其他縣(區(qū)、市)綜合質(zhì)量得分均有不同程度增加。其中，淳安縣和建德市增幅最明顯。這是由于隨著上下游補(bǔ)償協(xié)議的簽訂，生態(tài)環(huán)境修復(fù)與保護(hù)投入增加，鄉(xiāng)村人居環(huán)境綜合質(zhì)量逐漸提升。而屯溪區(qū)為黃山市主城區(qū)，受生態(tài)補(bǔ)償機(jī)制的影響較小。

第2輪試點(diǎn)(2015—2017年)期間，新安江流域生態(tài)補(bǔ)償機(jī)制逐步完善，生態(tài)保護(hù)成效明顯，生態(tài)效益迅速轉(zhuǎn)化為經(jīng)濟(jì)效益與社會效益，休寧縣、徽州區(qū)和黟縣鄉(xiāng)村人居環(huán)境綜合質(zhì)量大幅提升?？兿h和歙縣經(jīng)濟(jì)基礎(chǔ)相對薄弱，基礎(chǔ)設(shè)施處于相對落后水平，在第2輪生態(tài)補(bǔ)償試點(diǎn)期間，績溪縣和歙縣因加強(qiáng)生態(tài)環(huán)境保護(hù)而導(dǎo)致自身發(fā)展機(jī)會成本增加，所以績溪縣和歙縣鄉(xiāng)村人居環(huán)境總體質(zhì)量反而有所下降，得分分別下降0.01和0.02。受第2輪中央財政資金退坡的影響，淳安縣鄉(xiāng)村人居環(huán)境總體質(zhì)量下降明顯。第3輪試點(diǎn)(2018—2020年)期間，皖浙兩省實(shí)施橫向生態(tài)補(bǔ)償2+2億元，研究區(qū)鄉(xiāng)村人居環(huán)境綜合質(zhì)量出現(xiàn)明顯波動，其中，祁門縣、淳安縣和建德市上升明顯，上升幅度分別為9.42%、5.10%和6.96%，而黟縣出現(xiàn)一定幅度下降，降幅達(dá)3.85%。

2.2 鄉(xiāng)村人居環(huán)境質(zhì)量空間分布特征分析

為進(jìn)一步揭示試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量空間分布特征，采用ArcGIS 10.2軟件自然斷點(diǎn)法將2011—2019年研究區(qū)各縣(區(qū)、市)鄉(xiāng)村人居環(huán)境及其子系統(tǒng)評價結(jié)果分為高、較高、中等、較低和低水平5個等級，并將分類結(jié)果進(jìn)行可視化表達(dá)(圖4)。

圖4顯示，從鄉(xiāng)村人居環(huán)境綜合質(zhì)量空間分布來看，試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境及其子系統(tǒng)得分的空間分布較為分散，具有一定梯度性，在整體上呈現(xiàn)東西高、中間低的空間分布特征，形成以淳安縣和建德市為核心的鄉(xiāng)村人居環(huán)境質(zhì)量高水平集聚區(qū)，休寧縣-黃山區(qū)南北連片分布的較高水平區(qū)，屯溪區(qū)和祁門縣屬于中等水平區(qū)，黟縣屬于較低水平區(qū)，歙縣和績溪縣屬于低水平區(qū)。從所屬行政區(qū)來看，各地級市鄉(xiāng)村人居環(huán)境質(zhì)量由高到低依次為杭州市(建德市和淳安縣)、黃山市(黃山區(qū)、休寧縣、屯溪區(qū)、祁門縣、黟縣、徽州區(qū)和歙縣)和宣城市(績溪縣)。

從各子系統(tǒng)空間分布來看，鄉(xiāng)村人居環(huán)境高、較高和中等水平集聚區(qū)與鄉(xiāng)村生態(tài)環(huán)境高、較高和中水平集聚區(qū)基本吻合〔圖4(a)～(b)〕，說明鄉(xiāng)村生態(tài)環(huán)境是影響鄉(xiāng)村人居環(huán)境質(zhì)量的主要因素。從鄉(xiāng)村生態(tài)環(huán)境子系統(tǒng)空間分布來看，黃山區(qū)依托黃山風(fēng)景名勝區(qū)等優(yōu)勢生態(tài)資源，鄉(xiāng)村生態(tài)環(huán)境質(zhì)量屬于高水平區(qū)。休寧縣、淳安縣和建德市鄉(xiāng)村生態(tài)環(huán)境質(zhì)量屬于較高水平區(qū)。生態(tài)補(bǔ)償政策的實(shí)施極大改善了上游地區(qū)生態(tài)環(huán)境，緩解了下游地區(qū)環(huán)境保護(hù)和生態(tài)修復(fù)壓力，鄉(xiāng)村生態(tài)環(huán)境質(zhì)量水平較高。從鄉(xiāng)村生產(chǎn)環(huán)境子系統(tǒng)空間分布來看，淳安縣鄉(xiāng)村生產(chǎn)環(huán)境質(zhì)量屬于高水平區(qū)。淳安縣近年來將礦泉水加工、農(nóng)副產(chǎn)品生產(chǎn)與加工同電子商務(wù)相結(jié)合；同時，圍繞千島湖風(fēng)景區(qū)大力發(fā)展旅游產(chǎn)業(yè)，探索生態(tài)經(jīng)濟(jì)化、經(jīng)濟(jì)生態(tài)化的有效路徑?？兿h和祁門縣生態(tài)本底優(yōu)良，但經(jīng)濟(jì)基礎(chǔ)薄弱，生態(tài)補(bǔ)償政策限制了污染型工礦產(chǎn)業(yè)的發(fā)展，發(fā)展條件受到約束，鄉(xiāng)村生產(chǎn)環(huán)境水平較低。從鄉(xiāng)村生活環(huán)境子系統(tǒng)空間分布來看，建德市、黃山區(qū)、淳安縣、歙縣和休寧縣鄉(xiāng)村生活環(huán)境質(zhì)量在中等水平之上。其中，建德市和淳安縣經(jīng)濟(jì)較為發(fā)達(dá)且生態(tài)文明建設(shè)起步早，生活保障及基礎(chǔ)設(shè)施較為完善；黃山區(qū)和休寧縣由于政府補(bǔ)貼較多及鄉(xiāng)村人居環(huán)境整治力度大，鄉(xiāng)村生活環(huán)境質(zhì)量較高。

圖4 2011—2019年鄉(xiāng)村人居環(huán)境質(zhì)量及其子系統(tǒng)空間分布

2.3 鄉(xiāng)村人居環(huán)境時空演變的影響機(jī)制分析

2.3.1影響因素的選擇與檢驗(yàn)

鄉(xiāng)村人居環(huán)境時空演變受到眾多因素綜合影響，尤其是新安江生態(tài)補(bǔ)償機(jī)制對鄉(xiāng)村人居環(huán)境的影響機(jī)制尚不明晰。為探索生態(tài)補(bǔ)償機(jī)制背景下各因素對鄉(xiāng)村人居環(huán)境質(zhì)量差異的影響，基于鄉(xiāng)村生態(tài)、生產(chǎn)和生活因素3個方面選擇生態(tài)補(bǔ)償、生態(tài)環(huán)境、鄉(xiāng)村人口、經(jīng)濟(jì)、基礎(chǔ)設(shè)施和公共服務(wù)設(shè)施等方面因素，通過相關(guān)性檢驗(yàn)和復(fù)共線性檢驗(yàn)(表2)，最終篩選出14個指標(biāo)。鄉(xiāng)村生態(tài)因素包括生態(tài)補(bǔ)償指數(shù)(X1)、新安江流域水質(zhì)(X2)、化肥施用強(qiáng)度(X3)和農(nóng)藥施用強(qiáng)度(X4)；鄉(xiāng)村生產(chǎn)因素包括鄉(xiāng)村從業(yè)人員比例(X5)、鄉(xiāng)村人均農(nóng)林牧漁業(yè)總產(chǎn)值(X6)、單位GDP能耗(X7)、鄉(xiāng)村旅游業(yè)收入比例(X8)和農(nóng)業(yè)用電量比例(X9)；鄉(xiāng)村生活因素包括鄉(xiāng)村人口占比(X10)、鄉(xiāng)村人均可支配收入(X11)、鄉(xiāng)村醫(yī)療設(shè)施水平(X12)、鄉(xiāng)村通訊設(shè)施水平(X13)和鄉(xiāng)村低保人數(shù)比例(X14)。

表2顯示，14個因素P值均<0.05，且方差膨脹因子值(VIF)均小于4，遠(yuǎn)小于7.5，表明各解釋變量之間多重共線性現(xiàn)象較弱，可以進(jìn)行回歸分析[36]。時間加權(quán)回歸(TWR)模型、地理加權(quán)回歸(GWR)模型和時空地理加權(quán)回歸(GTWR)模型分別側(cè)重于時間作用、空間作用和時空共同作用，因此各模型擬合結(jié)果存在差異。由表3可知，各模型解釋力由強(qiáng)到弱依次為GTWR、GWR、TWR和全局普通線性回歸(Globe-OLS)模型。其中，GTWR模型R2及調(diào)整R2均為0.98，回歸模型擬合度高；GTWR模型AICc和RSS值最小，觀測數(shù)據(jù)解釋力最強(qiáng)；GTWR模型Bandwidth值也最小，模型平滑程度最高。因此，綜合考慮了時空共同作用的GTWR模型為該研究最優(yōu)選擇。

表2 影響因素相關(guān)性檢驗(yàn)和復(fù)共線性檢驗(yàn)結(jié)果

表3 Globe-OLS、TWR、GWR和GTWR模型回歸屬性對比

2.3.2影響因子空間分異分析

采用GTWR模型得到2011、2014、2017和2019年研究區(qū)鄉(xiāng)村人居環(huán)境影響參數(shù)估計結(jié)果，僅對2019年各影響因素參數(shù)估計值空間分異(表4)進(jìn)行分析。其中，生態(tài)補(bǔ)償指數(shù)(X1)、新安江流域水質(zhì)(X2)、鄉(xiāng)村人均農(nóng)林牧漁業(yè)總產(chǎn)值(X6)、鄉(xiāng)村人均可支配收入(X11)、鄉(xiāng)村醫(yī)療設(shè)施水平(X12)和鄉(xiāng)村低保人數(shù)比例(X14)是主要影響因素。

表4 鄉(xiāng)村人居環(huán)境質(zhì)量各影響因素參數(shù)估計值

從鄉(xiāng)村生態(tài)因素來看，生態(tài)補(bǔ)償指數(shù)和新安江流域水質(zhì)均為正值，說明新安江生態(tài)補(bǔ)償政策的實(shí)施明顯改善了新安江流域水質(zhì)和生態(tài)環(huán)境質(zhì)量，進(jìn)而促進(jìn)鄉(xiāng)村人居環(huán)境質(zhì)量提升?；屎娃r(nóng)藥施用強(qiáng)度均為負(fù)值，但數(shù)值較小，說明化肥和農(nóng)藥會在一定程度上破壞生態(tài)環(huán)境，從而影響鄉(xiāng)村人居環(huán)境質(zhì)量提升。從鄉(xiāng)村生產(chǎn)因素來看，除單位GDP能耗外，鄉(xiāng)村從業(yè)人員比例、鄉(xiāng)村人均農(nóng)林牧漁業(yè)總產(chǎn)值、鄉(xiāng)村旅游業(yè)收入比例和農(nóng)業(yè)用電量比例均為正值，可以促進(jìn)鄉(xiāng)村人居環(huán)境質(zhì)量提升。其中，鄉(xiāng)村人均農(nóng)林牧漁業(yè)總產(chǎn)值和農(nóng)業(yè)用電量比例影響較大，說明該區(qū)域大力發(fā)展生態(tài)農(nóng)業(yè)、綠色農(nóng)業(yè)可以促進(jìn)經(jīng)濟(jì)發(fā)展，進(jìn)而促進(jìn)鄉(xiāng)村人居環(huán)境質(zhì)量提升。從鄉(xiāng)村生活因素來看，鄉(xiāng)村低保人數(shù)比例是明顯短板，鄉(xiāng)村人口占比、鄉(xiāng)村人均可支配收入、鄉(xiāng)村醫(yī)療設(shè)施水平和鄉(xiāng)村通訊設(shè)施水平均為正值，其中，醫(yī)療衛(wèi)生水平和通訊水平改善可以顯著提升鄉(xiāng)村人居環(huán)境質(zhì)量。

2.3.3影響因子時序波動分析

為進(jìn)一步分析影響鄉(xiāng)村人居環(huán)境質(zhì)量各影響因素參數(shù)估計值的時序波動，對2011、2014、2017和2019年各變量對鄉(xiāng)村人居環(huán)境質(zhì)量影響參數(shù)的平均值進(jìn)行分析(圖5)。從鄉(xiāng)村生態(tài)因素(X1～X4)來看，隨著新安江流域生態(tài)補(bǔ)償機(jī)制的推進(jìn)，生態(tài)文明建設(shè)水平也逐漸提高，生態(tài)補(bǔ)償指數(shù)影響參數(shù)值正向提升，由2011年的9.25提升至2019年的10.43。

X1為生態(tài)補(bǔ)償指數(shù)，X2為新安江流域水質(zhì)，X3為化肥施用強(qiáng)度，X4為農(nóng)藥施用強(qiáng)度，X5為鄉(xiāng)村從業(yè)人員比例，X6為鄉(xiāng)村人均農(nóng)林牧漁業(yè)總產(chǎn)值，X7為單位GDP能耗，X8為鄉(xiāng)村旅游業(yè)收入比例，X9為農(nóng)業(yè)用電量比例，X10為鄉(xiāng)村人口占比，X11為鄉(xiāng)村人均可支配收入，X12為鄉(xiāng)村醫(yī)療設(shè)施水平，X13為鄉(xiāng)村通訊設(shè)施水平，X14為鄉(xiāng)村低保人數(shù)比例。

化肥和農(nóng)藥施用強(qiáng)度對鄉(xiāng)村人居環(huán)境質(zhì)量的抑制效果逐年提升，在生態(tài)環(huán)境持續(xù)向好的情況下，破壞生態(tài)環(huán)境的行為會明顯影響鄉(xiāng)村人居環(huán)境的提升。從鄉(xiāng)村生產(chǎn)因素(X5～X9)來看，隨著生態(tài)文明建設(shè)與供給側(cè)改革的推進(jìn)，由依靠增加生產(chǎn)投入來擴(kuò)大生產(chǎn)規(guī)模的粗放型經(jīng)濟(jì)增長方式轉(zhuǎn)變?yōu)橐揽靠萍紕?chuàng)新和高素質(zhì)勞動者來增數(shù)提質(zhì)的集約型經(jīng)濟(jì)增長方式。因此，單位GDP能耗的影響由正效應(yīng)轉(zhuǎn)為負(fù)效應(yīng)，2011年波峰值為0.02，而2019年波谷值為-0.77，總跨度達(dá)0.79。2011—2019年農(nóng)業(yè)用電量的影響由負(fù)效應(yīng)轉(zhuǎn)為正效應(yīng)，這說明隨著皖浙兩省之間合作的深入，新理念、新技術(shù)和新設(shè)備在農(nóng)業(yè)領(lǐng)域投入使用。農(nóng)業(yè)用電量增加反映生態(tài)農(nóng)業(yè)、綠色農(nóng)業(yè)和現(xiàn)代化農(nóng)業(yè)投入增加，促進(jìn)了農(nóng)村經(jīng)濟(jì)發(fā)展，進(jìn)而促進(jìn)鄉(xiāng)村人居環(huán)境提升。其他影響因素均為正效應(yīng)，且參數(shù)估計值呈現(xiàn)下降趨勢。從鄉(xiāng)村生活因素(X10～X14)來看，2011—2019年鄉(xiāng)村人口占比的影響由負(fù)效應(yīng)轉(zhuǎn)為正效應(yīng)，2011年波谷值為-0.41，2019年波峰值為0.21，總跨度為0.62。研究區(qū)城鄉(xiāng)二元結(jié)構(gòu)明顯，城鄉(xiāng)之間差距大。近年來，隨著生態(tài)補(bǔ)償機(jī)制、精準(zhǔn)扶貧和鄉(xiāng)村振興等政策的提出，村民逐漸成為鄉(xiāng)村建設(shè)與治理改善的主體，這可以促進(jìn)鄉(xiāng)村經(jīng)濟(jì)水平和人居環(huán)境的提升。鄉(xiāng)村人均可支配收入和鄉(xiāng)村低保人數(shù)比例的影響呈下降趨勢，鄉(xiāng)村醫(yī)療設(shè)施水平和鄉(xiāng)村通訊設(shè)施水平的影響呈上升趨勢。

2.3.4鄉(xiāng)村人居環(huán)境演變影響機(jī)制分析

試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量時空分異的形成是一個不斷發(fā)展且復(fù)雜多變的過程。經(jīng)過3輪試點(diǎn)，研究區(qū)鄉(xiāng)村人居環(huán)境賴以依存的生態(tài)、社會、經(jīng)濟(jì)和生活環(huán)境均發(fā)生較大變化。地方政府通過政策調(diào)控、生態(tài)治理等措施，直接或間接改善了鄉(xiāng)村人居環(huán)境。皖浙兩省通過構(gòu)建縱橫結(jié)合的補(bǔ)償機(jī)制，形成了共謀共治共建共享的聯(lián)動格局。在生態(tài)補(bǔ)償機(jī)制的影響下，鄉(xiāng)村人居環(huán)境發(fā)展對補(bǔ)償制度、政策及項(xiàng)目的依賴性不斷增強(qiáng)。生態(tài)、生產(chǎn)和生活環(huán)境受生態(tài)補(bǔ)償政策影響發(fā)生改變，進(jìn)而導(dǎo)致人居環(huán)境的空間差異進(jìn)一步擴(kuò)大。鄉(xiāng)村生態(tài)因素是試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量時空分異的主導(dǎo)因素，也是其形成與演化的“基石”，鄉(xiāng)村生產(chǎn)因素是試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量時空分異的內(nèi)在驅(qū)動因素，村民對于美好鄉(xiāng)村生活向往與現(xiàn)實(shí)環(huán)境的差距是提升鄉(xiāng)村人居環(huán)境質(zhì)量的直接驅(qū)動因素(圖6)。各方面因素的相互影響、相互作用，形成了試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量時空分布基本格局。

圖6 新安江-千島湖生態(tài)補(bǔ)償試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量時空分異影響機(jī)制

農(nóng)戶空間行為是影響鄉(xiāng)村人居環(huán)境的外在表現(xiàn)，在不同程度上影響了鄉(xiāng)村人居環(huán)境演化。生態(tài)補(bǔ)償機(jī)制實(shí)施之前，農(nóng)戶對生態(tài)環(huán)境保護(hù)的認(rèn)知不足，日常行為對人居環(huán)境質(zhì)量存在一定程度的破壞，人居環(huán)境質(zhì)量發(fā)展較慢。隨著生態(tài)補(bǔ)償機(jī)制的不斷推進(jìn)，農(nóng)戶空間行為在政策的約束和引導(dǎo)下不斷發(fā)生轉(zhuǎn)型調(diào)整。農(nóng)村生產(chǎn)、生活和生態(tài)空間格局趨于優(yōu)化，鄉(xiāng)村人居環(huán)境質(zhì)量不斷提高。隨著新安江-千島湖生態(tài)補(bǔ)償制度進(jìn)入長效階段，準(zhǔn)確把握試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境的影響機(jī)制，落實(shí)生態(tài)補(bǔ)償政策對農(nóng)戶空間行為的引導(dǎo)與規(guī)范，有利于進(jìn)一步共建生態(tài)文明共同體和一體化發(fā)展格局，共同促進(jìn)生態(tài)補(bǔ)償試驗(yàn)區(qū)人居環(huán)境快速發(fā)展。

3 討論

筆者從生態(tài)、生產(chǎn)和生活3個維度構(gòu)建指標(biāo)體系對試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境進(jìn)行評價，探討試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境時空演變過程，但在以下2個方面研究尚不充分：(1)研究區(qū)域范圍尚未完全明確。筆者研究區(qū)域?yàn)樾掳步饔?0個縣(區(qū)、市)，而政府部門尚未完全確定新安江-千島湖生態(tài)補(bǔ)償試驗(yàn)區(qū)邊界，其中，黃山區(qū)、休寧縣、黟縣、祁門縣、績溪縣和建德市的部分地區(qū)屬于試驗(yàn)區(qū)。但該研究使用上述6個縣(區(qū))全域的行政范圍和統(tǒng)計數(shù)據(jù)進(jìn)行制圖表達(dá)和分析。(2)鄉(xiāng)村人居環(huán)境各要素之間的交互影響機(jī)制分析尚不充分。雖然筆者采用時空地理加權(quán)回歸模型對鄉(xiāng)村人居環(huán)境影響要素及影響機(jī)制進(jìn)行分析，但受限于指標(biāo)數(shù)量及統(tǒng)計數(shù)據(jù)缺失等原因，未能充分研究各要素之間的交互影響。如何分析各影響要素對鄉(xiāng)村人居環(huán)境的交互作用，進(jìn)而從更深層次揭示鄉(xiāng)村人居環(huán)境發(fā)展對影響因素的響應(yīng)機(jī)制，是今后研究的重點(diǎn)方向。

4 結(jié)論與建議

4.1 結(jié)論

采用熵權(quán)TOPSIS法和時空地理加權(quán)回歸模型對2011—2019年新安江-千島湖生態(tài)補(bǔ)償試驗(yàn)區(qū)所在10個縣(區(qū)、市)鄉(xiāng)村人居環(huán)境質(zhì)量進(jìn)行時空分異及影響機(jī)制研究。得出以下結(jié)論：

(1)試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量的時間演變特征。2011—2019年試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境綜合得分以及鄉(xiāng)村生態(tài)環(huán)境和鄉(xiāng)村生活環(huán)境得分分別由生態(tài)補(bǔ)償機(jī)制實(shí)施前2011年的1.454、0.485和0.478上升到2019年的1.541、0.573和0.499，試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境及其子系統(tǒng)質(zhì)量總體呈上升趨勢。

(2)鄉(xiāng)村人居環(huán)境的空間分布特征。從鄉(xiāng)村人居環(huán)境質(zhì)量空間分布來看，試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境及其子系統(tǒng)質(zhì)量空間分布較為分散，具有一定梯度性，在整體上呈現(xiàn)東西高、中部低的空間分布特征。形成了以淳安縣和建德市為核心的鄉(xiāng)村人居環(huán)境質(zhì)量高水平區(qū)集聚區(qū)，休寧縣-黃山區(qū)南北連片分布的較高水平區(qū)，屯溪區(qū)和祁門縣屬于中等水平區(qū)，黟縣屬于較低水平區(qū)，歙縣和績溪縣屬于低水平區(qū)。從所屬行政區(qū)來看，各地級市鄉(xiāng)村人居環(huán)境質(zhì)量由高到低依次為杭州市(建德市和淳安縣)、黃山市(黃山區(qū)、休寧縣、屯溪區(qū)、祁門縣、黟縣、徽州區(qū)和歙縣)和宣城市(績溪縣)。

(3)鄉(xiāng)村人居環(huán)境質(zhì)量影響因子分析。鄉(xiāng)村人居環(huán)境時空演變受到多重因素的綜合影響，生態(tài)環(huán)境因素是主導(dǎo)因素，生產(chǎn)和生活因素是內(nèi)在驅(qū)動因素，政策因素是外部直接因素。3大因素相互影響形成了試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量時空分異格局。其中，生態(tài)補(bǔ)償指數(shù)、新安江流域水質(zhì)、鄉(xiāng)村人均農(nóng)林牧漁業(yè)總產(chǎn)值、鄉(xiāng)村人均可支配收入、鄉(xiāng)村醫(yī)療設(shè)施水平和鄉(xiāng)村低保人數(shù)比例6個因素是研究期內(nèi)驅(qū)動試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量上升的主要影響因子，且影響力隨著時間推移呈上升趨勢。

4.2 建議

(1)貫徹落實(shí)新安江流域生態(tài)補(bǔ)償制度，優(yōu)化鄉(xiāng)村生態(tài)環(huán)境。新安江流域生態(tài)補(bǔ)償機(jī)制不僅給流域各縣(區(qū)、市)帶來了生態(tài)、生產(chǎn)和生活等方面效益，同時也是生態(tài)文明建設(shè)的積極探索和有效實(shí)踐。隨著生態(tài)補(bǔ)償政策3輪試點(diǎn)的順利完成，研究區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量不斷提升。僅依靠新安江上游限制經(jīng)濟(jì)發(fā)展和開展生態(tài)修復(fù)，不利于全流域的統(tǒng)籌和諧及可持續(xù)發(fā)展。應(yīng)繼續(xù)貫徹落實(shí)新安江流域生態(tài)補(bǔ)償制度，充分考慮新安江上游地區(qū)生態(tài)環(huán)境保護(hù)成本、治理成本和經(jīng)濟(jì)發(fā)展成本，提高生態(tài)補(bǔ)償資金標(biāo)準(zhǔn)。應(yīng)落實(shí)“清單制”管理模式，提高相關(guān)驗(yàn)收標(biāo)準(zhǔn)，建立相關(guān)監(jiān)督檢查機(jī)制，以提高鄉(xiāng)村生態(tài)環(huán)境質(zhì)量，從而促進(jìn)試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量提高。

(2)大力發(fā)展生態(tài)經(jīng)濟(jì)，改善鄉(xiāng)村生產(chǎn)環(huán)境。應(yīng)充分利用黃山風(fēng)景區(qū)、徽州傳統(tǒng)村落和千島湖等優(yōu)勢自然資源及底蘊(yùn)深厚的徽州文化，發(fā)展生態(tài)旅游業(yè)和生態(tài)產(chǎn)業(yè)等多種經(jīng)營模式。充分挖掘水源地綠色發(fā)展新動能，更好地把生態(tài)效益轉(zhuǎn)化為經(jīng)濟(jì)效益和社會效益。積極引入民間組織、資本、技術(shù)和服務(wù)，在農(nóng)業(yè)生產(chǎn)中減少農(nóng)藥化肥等污染物使用量，構(gòu)建綠色農(nóng)業(yè)生產(chǎn)體系，優(yōu)化鄉(xiāng)村生產(chǎn)環(huán)境質(zhì)量，從而促進(jìn)試驗(yàn)區(qū)鄉(xiāng)村人居環(huán)境質(zhì)量提高。

(3)完善鄉(xiāng)村基礎(chǔ)建設(shè)，提高鄉(xiāng)村生活環(huán)境。各地經(jīng)濟(jì)水平、發(fā)展模式、人口分布和歷史文化都不同，應(yīng)結(jié)合當(dāng)?shù)噩F(xiàn)狀，將鄉(xiāng)村人居環(huán)境整治行動與鄉(xiāng)村經(jīng)濟(jì)發(fā)展水平相適應(yīng)。在新安江流域上游經(jīng)濟(jì)相對薄弱的地區(qū)，政府應(yīng)加大生態(tài)補(bǔ)償資金投入，加強(qiáng)基礎(chǔ)設(shè)施和鄉(xiāng)村居民生活保障設(shè)施建設(shè)；而在經(jīng)濟(jì)條件較好的地區(qū)，應(yīng)推動供給側(cè)改革，加快經(jīng)濟(jì)結(jié)構(gòu)轉(zhuǎn)型升級，大力提高鄉(xiāng)村生活環(huán)境質(zhì)量，促進(jìn)生態(tài)補(bǔ)償試驗(yàn)區(qū)人居環(huán)境快速發(fā)展。