珠三角農(nóng)村產(chǎn)業(yè)發(fā)展的水環(huán)境動(dòng)態(tài)效應(yīng)分析
——以東莞市水鄉(xiāng)片為例

甘建文，劉小玲(華南農(nóng)業(yè)大學(xué)公共管理學(xué)院，廣州 510642)

0 引 言

隨著珠三角地區(qū)農(nóng)村城鎮(zhèn)化進(jìn)程的不斷推進(jìn)，農(nóng)村產(chǎn)業(yè)快速發(fā)展導(dǎo)致污染排放急劇增加，農(nóng)藥化肥使用、鄉(xiāng)鎮(zhèn)工業(yè)廢水排放及生活污水排放等對(duì)水體造成污染，珠三角日益劣化的農(nóng)村水環(huán)境已經(jīng)長期為人詬病。根據(jù)2014年廣州社情民意研究中心民調(diào)數(shù)據(jù)顯示，廣東村民對(duì)水環(huán)境表示“不太滿意”和“不滿意”的受訪者自2012年以來都在三成及以上，其中珠三角村民對(duì)水環(huán)境不滿意程度最高。

從產(chǎn)業(yè)發(fā)展對(duì)水環(huán)境的影響機(jī)制出發(fā)，探討協(xié)調(diào)產(chǎn)業(yè)發(fā)展和水環(huán)境保護(hù)的路徑在學(xué)術(shù)界已經(jīng)取得一定的研究成果。國外研究多從微觀角度切入，分析水環(huán)境某一污染指標(biāo)對(duì)產(chǎn)業(yè)結(jié)構(gòu)或單一產(chǎn)業(yè)變動(dòng)的響應(yīng)，提出流域水環(huán)境治理模型和建議[1,2]。國內(nèi)研究則大多以黃河、淮河、渭河、太湖、長江等主要流域?yàn)檠芯繉?duì)象，通過建立產(chǎn)業(yè)污染排放與流域水質(zhì)變化的定量模型，綜合分析流域產(chǎn)業(yè)發(fā)展與水環(huán)境質(zhì)量變化的內(nèi)在關(guān)系[3-5]。然而，目前大部分的研究傾向于分析城市工業(yè)污水排放與區(qū)域水體水質(zhì)的關(guān)系，而著眼于解決農(nóng)村產(chǎn)業(yè)發(fā)展和農(nóng)村水環(huán)境保護(hù)之間矛盾的研究并不多見。

農(nóng)村產(chǎn)業(yè)發(fā)展作為區(qū)域經(jīng)濟(jì)發(fā)展環(huán)境負(fù)荷的重要影響要素之一，通過調(diào)整產(chǎn)業(yè)規(guī)模、優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)及改進(jìn)生產(chǎn)技術(shù)的方式來實(shí)現(xiàn)農(nóng)村產(chǎn)業(yè)經(jīng)濟(jì)與資源環(huán)境協(xié)調(diào)發(fā)展，已經(jīng)成為全社會(huì)的普遍共識(shí)。產(chǎn)業(yè)生態(tài)化是協(xié)調(diào)產(chǎn)業(yè)發(fā)展和環(huán)境保護(hù)的重要理論支撐。為此，本文在理清產(chǎn)業(yè)生態(tài)化理論內(nèi)涵的基礎(chǔ)上，定性分析珠三角農(nóng)村產(chǎn)業(yè)發(fā)展對(duì)水環(huán)境作用機(jī)理的基礎(chǔ)上，以東莞市水鄉(xiāng)片為實(shí)證區(qū)域，通過基于VAR模型的脈沖響應(yīng)函數(shù)分析方法來重點(diǎn)探討農(nóng)村產(chǎn)業(yè)發(fā)展的水環(huán)境動(dòng)態(tài)效應(yīng)過程，以期為政府制定珠三角農(nóng)村產(chǎn)業(yè)發(fā)展及水環(huán)境保護(hù)政策提供現(xiàn)實(shí)信息支撐。

1 產(chǎn)業(yè)生態(tài)化理論的提出及其內(nèi)涵

現(xiàn)代社會(huì)生態(tài)脅迫型和環(huán)境破壞型產(chǎn)業(yè)經(jīng)濟(jì)增長模式導(dǎo)致了資源枯竭和生態(tài)退化，引起了人類對(duì)產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展方式的反思。1989年，Robert Frosch與 Nicholas E. Gallopoulos 在《制造業(yè)的戰(zhàn)略》一文中[6]提出通過模擬自然生態(tài)系統(tǒng)的物質(zhì)循環(huán)來構(gòu)建工業(yè)生產(chǎn)系統(tǒng)，在系統(tǒng)中生產(chǎn)企業(yè)相互依存、互相聯(lián)系，按照生態(tài)規(guī)律和經(jīng)濟(jì)規(guī)律進(jìn)行生產(chǎn)活動(dòng)，通過閉路循環(huán)過程達(dá)到最大程度利用資源和減少廢棄物排放的目的，實(shí)現(xiàn)經(jīng)濟(jì)與環(huán)境的雙贏。隨后Nicholas等人進(jìn)一步從生態(tài)系統(tǒng)的角度提出產(chǎn)業(yè)生態(tài)系統(tǒng)和產(chǎn)業(yè)生態(tài)學(xué)的概念。1991年美國國家科學(xué)院與貝爾實(shí)驗(yàn)室共同組織產(chǎn)業(yè)生態(tài)學(xué)論壇，對(duì)產(chǎn)業(yè)生態(tài)學(xué)的概念、內(nèi)容和方法以及應(yīng)用前景進(jìn)行全面系統(tǒng)的總結(jié)，認(rèn)為產(chǎn)業(yè)生態(tài)學(xué)是對(duì)各種產(chǎn)業(yè)活動(dòng)及其產(chǎn)品與環(huán)境之間的相互關(guān)系的跨學(xué)科研究[7]。

產(chǎn)業(yè)生態(tài)學(xué)理論是研究以生態(tài)學(xué)原理構(gòu)建產(chǎn)業(yè)發(fā)展的閉路循環(huán)系統(tǒng)，對(duì)整個(gè)生產(chǎn)活動(dòng)的物質(zhì)周期過程加以優(yōu)化協(xié)調(diào)的系統(tǒng)方法集合。產(chǎn)業(yè)生態(tài)學(xué)建立了三次產(chǎn)業(yè)之間的聯(lián)系，通過工藝改進(jìn)、污染減量、產(chǎn)品生命周期分析、代謝分析等技術(shù)和分析工具。產(chǎn)業(yè)生態(tài)化則是指在產(chǎn)業(yè)生態(tài)學(xué)理論指導(dǎo)下產(chǎn)業(yè)發(fā)展的高級(jí)形態(tài)，通過模仿自然生態(tài)系統(tǒng)閉路循環(huán)的模式構(gòu)建產(chǎn)業(yè)系統(tǒng)，按照生態(tài)規(guī)律和經(jīng)濟(jì)規(guī)律安排生產(chǎn)活動(dòng)，以實(shí)現(xiàn)產(chǎn)業(yè)的資源高效利用、環(huán)境污染最小、廢棄物循環(huán)利用和社會(huì)經(jīng)濟(jì)效益最大的目的。產(chǎn)業(yè)生態(tài)化理論為協(xié)調(diào)農(nóng)村產(chǎn)業(yè)發(fā)展和水環(huán)境保護(hù)之間的關(guān)系、制定環(huán)境經(jīng)濟(jì)政策與環(huán)境法律制度等提供重要指導(dǎo)，近年來在清潔生產(chǎn)、生態(tài)產(chǎn)業(yè)集群建設(shè)、整合梳理產(chǎn)業(yè)鏈條、產(chǎn)品廢棄物循環(huán)利用、產(chǎn)業(yè)科技發(fā)展等方面應(yīng)用廣泛并取得了豐富經(jīng)驗(yàn)。

2 東莞市水鄉(xiāng)片農(nóng)村產(chǎn)業(yè)發(fā)展及水環(huán)境概況

東莞市水鄉(xiāng)片在區(qū)位條件、地質(zhì)水文、產(chǎn)業(yè)經(jīng)濟(jì)及社會(huì)發(fā)展等方面具有珠三角農(nóng)村的典型特征，而且其文化、自然、社會(huì)等因素差異較小，為研究區(qū)域污染水平與經(jīng)濟(jì)因素的理想微觀尺度和實(shí)證樣本。東莞市水鄉(xiāng)片位于東莞市西北部，地處珠三角地區(qū)的中心腹地，屬于東江河口三角洲河網(wǎng)平原地區(qū)，在珠三角的水域面積最大，河網(wǎng)密度居于前列[8]。為了研究東莞市水鄉(xiāng)片農(nóng)村產(chǎn)業(yè)發(fā)展對(duì)水環(huán)境的影響情況，本文選取石龍鎮(zhèn)、中堂鎮(zhèn)、望牛墩鎮(zhèn)等九個(gè)鎮(zhèn)為研究對(duì)象，不包括萬江街道及虎門港。

2.1 東莞市水鄉(xiāng)片農(nóng)村產(chǎn)業(yè)發(fā)展概況

從1990-2013年東莞市水鄉(xiāng)片產(chǎn)業(yè)結(jié)構(gòu)變動(dòng)情況上看(圖1)，第一、二、三產(chǎn)業(yè)產(chǎn)值占比結(jié)構(gòu)由1990年的29∶54∶17發(fā)展為2013年的1∶61∶38。在這24年間，非農(nóng)產(chǎn)業(yè)比重呈現(xiàn)逐年上升的趨勢(shì)，其中第二產(chǎn)業(yè)在東莞市水鄉(xiāng)片產(chǎn)業(yè)結(jié)構(gòu)中長期保持主導(dǎo)地位，這種第二產(chǎn)業(yè)產(chǎn)值比重過大的產(chǎn)業(yè)結(jié)構(gòu)特征在當(dāng)前珠三角農(nóng)村產(chǎn)業(yè)發(fā)展過程中普遍存在，由于較城市而言，農(nóng)村地區(qū)主要依靠土地和低水平人力要素投入的擴(kuò)大來實(shí)現(xiàn)農(nóng)村經(jīng)濟(jì)的持續(xù)增長，這一種粗放型的增長方式最終導(dǎo)致環(huán)境的惡化。第三產(chǎn)業(yè)比重與第一產(chǎn)業(yè)呈反向變動(dòng)，而且該地區(qū)的第三產(chǎn)業(yè)成分具有以生活性服務(wù)部門為主的特征。

圖1 1990-2013年東莞市水鄉(xiāng)片產(chǎn)業(yè)結(jié)構(gòu)變動(dòng)情況Fig.1 Industrial structure changes of Dongguan Watery Region during 1990-2013

2.2 東莞市水鄉(xiāng)片水環(huán)境概況

隨著農(nóng)村工業(yè)化和城鎮(zhèn)化進(jìn)程的不斷推進(jìn)，化肥農(nóng)藥過度使用、鄉(xiāng)鎮(zhèn)企業(yè)污水排放及農(nóng)村居民生活廢物的隨意傾倒對(duì)農(nóng)村水環(huán)境造成嚴(yán)重破壞，東莞市水鄉(xiāng)片的產(chǎn)業(yè)發(fā)展及水環(huán)境現(xiàn)狀已經(jīng)成為珠三角過去30年傳統(tǒng)發(fā)展模式的縮影。據(jù)統(tǒng)計(jì)數(shù)據(jù)顯示，2012 年東莞市水鄉(xiāng)地區(qū)COD、SO2排放量占全市的比重高達(dá)91%、52%，其能耗也超過全市的30%，總體水質(zhì)狀況為中度污染。根據(jù)《2013年東莞市環(huán)境狀況公報(bào)》數(shù)據(jù)顯示，在納入監(jiān)測(cè)的水鄉(xiāng)地區(qū)13條主要河流中，Ⅴ類及劣Ⅴ類水質(zhì)河流占比高達(dá)53.8%，水質(zhì)目標(biāo)達(dá)標(biāo)率僅為15.4%(表1)，超標(biāo)項(xiàng)目以氨氮、總磷、化學(xué)需氧量等指標(biāo)為主，污染源具有農(nóng)業(yè)面源污染、工業(yè)廢水污染和生活污水污染的復(fù)合污染特征。

3 研究設(shè)計(jì)

3.1 研究方法與步驟

經(jīng)濟(jì)發(fā)展和生態(tài)環(huán)境之間存在著復(fù)雜的相互依賴、相互影響、相互制約的互動(dòng)關(guān)系，生態(tài)經(jīng)濟(jì)復(fù)合系統(tǒng)內(nèi)部具有其特殊的結(jié)構(gòu)、功能及其運(yùn)動(dòng)規(guī)律，向量自回歸模型(Vector Autoregressive Model, VAR)作為多變量時(shí)間序列的一種方法，同時(shí)能很好地解決內(nèi)生變量偏誤和變量缺損的問題。另外，為了分析一個(gè)變量對(duì)另一個(gè)變量的影響情況，本文采用脈沖響應(yīng)函數(shù)(Impulse response function, IRF)來衡量任意一個(gè)變量的擾動(dòng)對(duì)其他變量當(dāng)前和未來取值的影響。為了消除變量排序?qū)Ψ治鼋Y(jié)果的干擾，本文采用改進(jìn)的廣義脈沖響應(yīng)函數(shù)(Generalized impulse response function, GIRF)進(jìn)行分析。

表1 東莞市地表水功能區(qū)劃水質(zhì)目標(biāo)及其達(dá)標(biāo)情況Tab.1 Surface water quality objectives and compliance

注：①超標(biāo)項(xiàng)目以氨氮、總磷、化學(xué)需氧量等指標(biāo)為主。②數(shù)據(jù)來源：《東莞市環(huán)境狀況公報(bào)2013》、《東莞市地表水功能區(qū)劃成果表》，東莞市環(huán)保局。

本文首先對(duì)東莞市水鄉(xiāng)片1990-2013年的水環(huán)境進(jìn)行評(píng)價(jià)，以此來考察農(nóng)村水環(huán)境污染的變化過程。接著對(duì)觀察變量進(jìn)行平穩(wěn)性和協(xié)整性檢驗(yàn)，確定變量間是否存在長期穩(wěn)定的均衡關(guān)系。最后采用基于VAR模型的脈沖響應(yīng)函數(shù)分析方法來考察變量間的動(dòng)態(tài)影響效應(yīng)。

3.2 變量選取與數(shù)據(jù)來源

珠三角農(nóng)村水環(huán)境受到了農(nóng)業(yè)面源污染、工業(yè)污水和生活廢水的復(fù)合影響，其中以工業(yè)污水和生活廢水為主要污染源[9-11]。在綜合考慮數(shù)據(jù)的可得性和可靠性的基礎(chǔ)上，鑒于本文的主要目的是考察農(nóng)村產(chǎn)業(yè)發(fā)展對(duì)水環(huán)境污染總體水平的影響，此處選取東莞市水鄉(xiāng)片9鄉(xiāng)鎮(zhèn)1990-2013年化肥農(nóng)藥實(shí)際施用量、工業(yè)廢水排放量及生活污水排放量3個(gè)污染物排放總量指標(biāo)來測(cè)度東莞市水鄉(xiāng)片水環(huán)境污染指數(shù)(SWRZS)，同時(shí)該指數(shù)作為實(shí)證分析的因變量 。

產(chǎn)業(yè)發(fā)展作為社會(huì)經(jīng)濟(jì)發(fā)展的主要驅(qū)動(dòng)力，承擔(dān)了自然資源輸入向產(chǎn)品服務(wù)輸出的轉(zhuǎn)換過程，其對(duì)資源環(huán)境的影響是通過規(guī)模效應(yīng)、結(jié)構(gòu)效應(yīng)和技術(shù)效應(yīng)途徑實(shí)現(xiàn)的[12-14]，因此， 衡量農(nóng)村產(chǎn)業(yè)發(fā)展從產(chǎn)業(yè)規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)和技術(shù)進(jìn)步這3個(gè)自變量入手：首先，產(chǎn)業(yè)規(guī)模(CYGM)用地區(qū)生產(chǎn)總值表示，反映農(nóng)村地域范圍內(nèi)三次產(chǎn)業(yè)的整體經(jīng)營活動(dòng)成果；其次，由于產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化是指各產(chǎn)業(yè)部門之間比例的合理化和高級(jí)化，而第三產(chǎn)業(yè)比重的不斷擴(kuò)大是區(qū)域產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化的重要標(biāo)志之一，因此用第三產(chǎn)業(yè)產(chǎn)值占GDP的比重來表示產(chǎn)業(yè)結(jié)構(gòu)(CYJG)的變動(dòng)情況；最后，技術(shù)進(jìn)步(JSJB)用全員勞動(dòng)生產(chǎn)率表示，反映區(qū)域內(nèi)各企業(yè)的生產(chǎn)技術(shù)、經(jīng)營管理、職工技能的綜合水平。

以上所有基礎(chǔ)數(shù)據(jù)均來自于歷年《東莞統(tǒng)計(jì)年鑒》 ，時(shí)序長度為1990-2013年，所有涉及產(chǎn)值的數(shù)據(jù)均以1990年為基期進(jìn)行指數(shù)平減，為了消除和減少數(shù)據(jù)中可能存在的異方差影響，本文實(shí)際分析時(shí)均采用各變量的對(duì)數(shù)值。資料的整理和處理使用Excel2003和Eviews6.0軟件。

4 東莞市水鄉(xiāng)片農(nóng)村產(chǎn)業(yè)發(fā)展與水環(huán)境動(dòng)態(tài)影響機(jī)理

4.1 農(nóng)村水環(huán)境污染綜合指數(shù)測(cè)算及其檢驗(yàn)

4.1.1農(nóng)村水環(huán)境污染綜合指數(shù)測(cè)算

本文選取東莞市水鄉(xiāng)片9鄉(xiāng)鎮(zhèn)1990-2013年化肥農(nóng)藥實(shí)際施用量、工業(yè)廢水排放量及生活污水排放量3個(gè)指標(biāo)來測(cè)度東莞市水鄉(xiāng)片農(nóng)村水環(huán)境污染指數(shù)(SWRZS)。為了消除量綱和量綱單位的不同所帶來的不可公度性，采用式(1)對(duì)原始數(shù)據(jù)進(jìn)行無量綱化處理：

(1)

式中：i為年份；V′i為標(biāo)準(zhǔn)化后的賦值；Vi為各污染指標(biāo)的當(dāng)年值；Vmax為各污染指標(biāo)最大值；Vmin為各污染指標(biāo)的最小值。

根據(jù)珠三角農(nóng)村水環(huán)境污染源特征，結(jié)合東莞水鄉(xiāng)片的產(chǎn)業(yè)結(jié)構(gòu)比例，采用層次分析方法確定各污染指標(biāo)在計(jì)算水環(huán)境污染綜合指數(shù)中所占的權(quán)重(表2)。基于Excel2003對(duì)判斷矩陣進(jìn)行處理并進(jìn)行一致性檢驗(yàn)，得到化肥農(nóng)藥實(shí)際施用量、工業(yè)廢水排放量及生活污水排放量的權(quán)重分別為0.07、0.65與0.28，其中一致性檢驗(yàn)CR=0.056。

表2 水環(huán)境污染綜合指數(shù)測(cè)量指標(biāo)權(quán)重判斷矩陣Tab.2 The power matrix of SWRZS

最后，根據(jù)式(2)計(jì)算歷年水環(huán)境污染綜合指數(shù)：

(2)

式中：Yi為第i年的水環(huán)境污染綜合指數(shù)；j代表污染指標(biāo)種類；Yji為第i年第j種污染指標(biāo)的標(biāo)準(zhǔn)化值；Wj為第j種污染指標(biāo)的權(quán)重值。

水環(huán)境污染綜合指數(shù)如圖2所示。1990-2013年這20余年來，東莞市水鄉(xiāng)片9鄉(xiāng)鎮(zhèn)在農(nóng)村產(chǎn)業(yè)快速發(fā)展的同時(shí)，其水環(huán)境污染水平呈現(xiàn)先增后減的變化趨勢(shì)，符合環(huán)境庫茲涅茲曲線的理論假設(shè)。

圖2 1990-2013年東莞市水鄉(xiāng)片水環(huán)境污染綜合指數(shù)變化趨勢(shì)圖Fig.2 The SWRZS changes of Dongguan Watery Region during 1990-2013

4.1.2平穩(wěn)性檢驗(yàn)與協(xié)整檢驗(yàn)

為了防止偽回歸的產(chǎn)生，對(duì)時(shí)間序列數(shù)據(jù)進(jìn)行平穩(wěn)性檢驗(yàn)。在對(duì)所選取的時(shí)間序列進(jìn)行一階差分后，發(fā)現(xiàn)一階差分序列的ADF值都在1%顯著水平下拒絕原假設(shè)，即都是平穩(wěn)的。因此，CYGM、CYJG、JSJB及SWRZS序列均為一階單整序列，滿足協(xié)整檢驗(yàn)的前提條件，見表3。

表3 lnCYGM、lnCYJG、lnJSJB、lnSWRZS的協(xié)整檢驗(yàn)Tab.3 The cointegration of lnCYGM、lnCYJG、lnJSJB、lnSWRZS

注：* denotes rejection of the hypothesis at the 0.05 level；**MacKinnon-Haug-Michelis (1999) p-values。

在平穩(wěn)性檢驗(yàn)的基礎(chǔ)上還要考慮各序列之間是否存在長期協(xié)整關(guān)系，本文采用 Johansen 極大似然估計(jì)法對(duì)變量進(jìn)行協(xié)整檢驗(yàn)。從表4的檢驗(yàn)結(jié)果可以看出，CYGM、CYJG、JSJB與SWRZS之間有且只有一個(gè)協(xié)整方程，即存在唯一的長期動(dòng)態(tài)均衡關(guān)系，說明農(nóng)村產(chǎn)業(yè)發(fā)展水平的提高與農(nóng)村水環(huán)境污染指數(shù)之間具有長期的比例關(guān)系，即農(nóng)村產(chǎn)業(yè)發(fā)展與水環(huán)境污染息息相關(guān)。

表4 標(biāo)準(zhǔn)化的協(xié)整方程表Tab.4 Standardized cointegration equation

注：模型中括號(hào)內(nèi)為估計(jì)標(biāo)準(zhǔn)差；Log likelihood= 67.896 22。

由表3可以得出，標(biāo)準(zhǔn)化的協(xié)整方程表達(dá)式為(方程下方括號(hào)內(nèi)數(shù)字為標(biāo)準(zhǔn)誤差)：

由協(xié)整方程可以看出，各變量之間具有長期穩(wěn)定的均衡關(guān)系。從各變量的符號(hào)上看，農(nóng)村產(chǎn)業(yè)規(guī)模和產(chǎn)業(yè)結(jié)構(gòu)的發(fā)展在一定程度上加重了水環(huán)境污染，而技術(shù)進(jìn)步因素則可以減輕水環(huán)境的污染程度。由于協(xié)整模型估算而得的系數(shù)是原序列取對(duì)數(shù)之后的結(jié)果，系數(shù)估計(jì)值代表彈性，由此可得水環(huán)境污染的產(chǎn)業(yè)規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)、技術(shù)進(jìn)步彈性分別為1.902 233、0.225 730、-3.909 579，表明產(chǎn)業(yè)規(guī)模擴(kuò)大1%將會(huì)帶來1.9%水環(huán)境污染指數(shù)的上升；第三產(chǎn)業(yè)比重提高1%就會(huì)帶來0.23%水環(huán)境污染指數(shù)的上升；而技術(shù)進(jìn)步因素提高1%就能夠帶來3.91%水環(huán)境污染指數(shù)的下降。

4.2 基于VAR模型的廣義脈沖響應(yīng)分析

為了進(jìn)一步解釋農(nóng)村產(chǎn)業(yè)發(fā)展各自變量對(duì)水環(huán)境污染指數(shù)的動(dòng)態(tài)影響過程，有必要進(jìn)行廣義脈沖響應(yīng)分析，分析結(jié)果如圖3、表5所示。橫軸代表沖擊作用的滯后期數(shù)，縱軸表示因變量對(duì)各自變量的響應(yīng)大小，實(shí)線表示脈沖響應(yīng)曲線，虛線代表兩倍標(biāo)準(zhǔn)差的置信帶。

圖3 lnSWRZS對(duì)lnCYGM、lnCYJG、lnJSJB的脈沖響應(yīng)圖Fig.3 lnSWRZS impulse response to lnCYGM、lnCYJG、lnJSJB

如圖3和表5結(jié)果顯示，在整個(gè)沖擊響應(yīng)期內(nèi)，lnSWRZS對(duì)當(dāng)期lnCYGM一個(gè)單位沖擊的平均響應(yīng)值為0.61%。具體分時(shí)段來看，lnSWRZS對(duì)lnCYGM沖擊的當(dāng)期響應(yīng)值為1.37%，隨后迅速下降，到第2期達(dá)到其負(fù)的最大值-0.56%，并從第3期開始平穩(wěn)上升到第8期后穩(wěn)定在1.00%水平。該沖擊響應(yīng)過程說明：水環(huán)境污染指數(shù)對(duì)產(chǎn)業(yè)規(guī)模一個(gè)標(biāo)準(zhǔn)差沖擊的響應(yīng)總體為正，且存在明顯的時(shí)滯效應(yīng)。即東莞市水鄉(xiāng)片農(nóng)村產(chǎn)業(yè)規(guī)模擴(kuò)大先是在當(dāng)期對(duì)水環(huán)境產(chǎn)生最大的污染作用，隨后產(chǎn)業(yè)規(guī)模擴(kuò)大帶來經(jīng)濟(jì)效益的提高會(huì)對(duì)水環(huán)境產(chǎn)生一定正向效應(yīng)，但這種改善并不是持續(xù)的，緊接著農(nóng)村產(chǎn)業(yè)規(guī)模的不斷擴(kuò)大將導(dǎo)致水環(huán)境污染的進(jìn)一步加劇，且污染作用具有長期性。

lnSWRZS對(duì)當(dāng)期lnCYJG一個(gè)單位沖擊的平均響應(yīng)值為2.97%，在前4期呈現(xiàn)單調(diào)遞增的趨勢(shì)。具體分時(shí)段來看，lnSWRZS對(duì)lnCYJG沖擊的當(dāng)期響應(yīng)值為-8.10%，隨后逐漸增加，到第7期增至最大值5.85%，隨后穩(wěn)定于這一水平。該沖擊響應(yīng)過程說明：水環(huán)境污染指數(shù)對(duì)產(chǎn)業(yè)結(jié)構(gòu)的一個(gè)標(biāo)準(zhǔn)差沖擊的響應(yīng)為正，且表現(xiàn)出可持續(xù)的特征。即東莞市水鄉(xiāng)片農(nóng)村第三產(chǎn)業(yè)比重提高先是在當(dāng)期對(duì)水環(huán)境產(chǎn)生顯著的正向效應(yīng)，隨著第三產(chǎn)業(yè)比重的不斷增加，這種正向作用會(huì)逐漸下降，最終會(huì)導(dǎo)致水環(huán)境污染的加劇。

lnSWRZS對(duì)當(dāng)期lnJSJB一個(gè)單位沖擊的平均響應(yīng)值為2.33%。具體分時(shí)段來看，lnSWRZS對(duì)lnJSJB沖擊的當(dāng)期響應(yīng)值為0.07%，第2期達(dá)到其負(fù)的最大值-4.08%，隨后逐漸平穩(wěn)至-1.39%的水平。該沖擊響應(yīng)過程說明：水環(huán)境污染指數(shù)對(duì)技術(shù)進(jìn)步一個(gè)標(biāo)準(zhǔn)差沖擊的響應(yīng)總體為負(fù)，且表現(xiàn)出可持續(xù)的特征。即東莞市水鄉(xiāng)片農(nóng)村產(chǎn)業(yè)技術(shù)進(jìn)步有利于水環(huán)境質(zhì)量的改善，且該正向效應(yīng)長期有效。

表5 lnSWRZS對(duì)lnCYGM、lnCYJG、lnJSJB的脈沖響應(yīng)函數(shù)值Tab.5 The lnSWRZS impulse response function value tolnCYGM、lnCYJG、lnJSJB

5 結(jié)論與討論

5.1 主要結(jié)論

珠三角農(nóng)村產(chǎn)業(yè)發(fā)展受到城鄉(xiāng)二元結(jié)構(gòu)體制的影響，農(nóng)業(yè)發(fā)展效益低、工業(yè)大而不強(qiáng)、服務(wù)業(yè)單一化已經(jīng)成為珠三角農(nóng)村產(chǎn)業(yè)結(jié)構(gòu)的典型特征，這一種低層次、不均衡的產(chǎn)業(yè)發(fā)展方式最終會(huì)導(dǎo)致環(huán)境的惡化和發(fā)展內(nèi)生動(dòng)力不足。與此同時(shí)，缺乏創(chuàng)新驅(qū)動(dòng)轉(zhuǎn)型的能力也會(huì)致使珠三角農(nóng)村經(jīng)濟(jì)發(fā)展內(nèi)卷化，抵抗環(huán)境變化風(fēng)險(xiǎn)的水平不高。

珠三角農(nóng)村水環(huán)境劣化狀況不容樂觀。以東莞市水鄉(xiāng)片為研究區(qū)域，該片區(qū)水環(huán)境污染綜合指數(shù)在總體上呈現(xiàn)“倒U型”的變化趨勢(shì)，符合環(huán)境庫茲涅茲曲線的理論假設(shè)。農(nóng)村水環(huán)境受到了農(nóng)村產(chǎn)業(yè)規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)和生產(chǎn)技術(shù)的聯(lián)動(dòng)影響，它們之間存在長期穩(wěn)定的協(xié)整關(guān)系，并且具有唯一的長期動(dòng)態(tài)均衡關(guān)系。協(xié)整方程結(jié)果顯示，產(chǎn)業(yè)規(guī)模擴(kuò)大和產(chǎn)業(yè)結(jié)構(gòu)發(fā)展都會(huì)對(duì)水環(huán)境有負(fù)向效應(yīng)，且產(chǎn)業(yè)規(guī)模擴(kuò)大的水環(huán)境劣化效應(yīng)要比產(chǎn)業(yè)結(jié)構(gòu)發(fā)展的水環(huán)境劣化效應(yīng)大，生產(chǎn)技術(shù)的進(jìn)步則具有顯著的正向效應(yīng)。

從動(dòng)態(tài)效應(yīng)角度看，短期內(nèi)產(chǎn)業(yè)規(guī)模和產(chǎn)業(yè)結(jié)構(gòu)對(duì)水環(huán)境起到了雙刃劍的作用，但是從長期趨勢(shì)上看依然會(huì)導(dǎo)致水環(huán)境質(zhì)量的惡化。農(nóng)村產(chǎn)業(yè)規(guī)模增大對(duì)水環(huán)境污染的動(dòng)態(tài)效應(yīng)遵循“深刻加劇-稍微緩解-持續(xù)加劇”的過程。農(nóng)村產(chǎn)業(yè)結(jié)構(gòu)高級(jí)化對(duì)水環(huán)境污染的動(dòng)態(tài)效應(yīng)遵循“短期緩解-逐漸加劇”的過程，主要是因?yàn)橹槿寝r(nóng)村第三產(chǎn)業(yè)以生活性服務(wù)部門為主，尤其是以傳統(tǒng)服務(wù)業(yè)和旅游休閑服務(wù)業(yè)為主，導(dǎo)致了生活污水的增加，加劇水環(huán)境污染。盡管這樣，仍然無法否定產(chǎn)業(yè)結(jié)構(gòu)高級(jí)化對(duì)當(dāng)期水環(huán)境具有較強(qiáng)的正向效應(yīng)。農(nóng)村生產(chǎn)技術(shù)進(jìn)步對(duì)降低水環(huán)境污染仍然有持續(xù)的正向效應(yīng)。

5.2 政策啟示

(1)拉長產(chǎn)業(yè)鏈條，適度發(fā)展產(chǎn)業(yè)集群。對(duì)于珠三角農(nóng)村地區(qū)而言，農(nóng)村產(chǎn)業(yè)分散式經(jīng)營局面導(dǎo)致企業(yè)之間組織化程度低、生產(chǎn)聯(lián)系弱、協(xié)作配套差，無法在空間上產(chǎn)生集聚和整合效應(yīng)，最終造成嚴(yán)重的資源浪費(fèi)和重復(fù)建設(shè)問題。產(chǎn)業(yè)生態(tài)化要求產(chǎn)業(yè)發(fā)展形成閉路循環(huán)系統(tǒng)，對(duì)整個(gè)生產(chǎn)活動(dòng)的物質(zhì)周期過程加以優(yōu)化協(xié)調(diào)。拉長產(chǎn)業(yè)鏈條有利于農(nóng)業(yè)、工業(yè)、服務(wù)業(yè)之間的聯(lián)系和延伸，產(chǎn)生相互補(bǔ)充、資源重復(fù)利用的生產(chǎn)經(jīng)營群體，加強(qiáng)各經(jīng)營群體間相互協(xié)作，形成閉合循環(huán)的生產(chǎn)系統(tǒng)。而且，完整的產(chǎn)業(yè)系統(tǒng)與農(nóng)村水環(huán)境系統(tǒng)緊密結(jié)合起來，形成有機(jī)統(tǒng)一的農(nóng)村水環(huán)境網(wǎng)絡(luò)治理體系，實(shí)現(xiàn)污染減量。發(fā)展適度產(chǎn)業(yè)集群，有利于扭轉(zhuǎn)現(xiàn)有珠三角農(nóng)村產(chǎn)業(yè)獨(dú)立分布現(xiàn)狀，激活農(nóng)村經(jīng)濟(jì)發(fā)展活力，引導(dǎo)技術(shù)、資金、人才等生產(chǎn)要素流向農(nóng)村。

(2)合理開發(fā)布局，調(diào)整優(yōu)化農(nóng)村產(chǎn)業(yè)結(jié)構(gòu)。珠三角農(nóng)村產(chǎn)業(yè)發(fā)展受到城鄉(xiāng)二元結(jié)構(gòu)體制的影響，在產(chǎn)業(yè)轉(zhuǎn)型升級(jí)過程中一直處于被動(dòng)地位，只能通過依靠土地要素、生態(tài)資源及廉價(jià)勞動(dòng)力投入來實(shí)現(xiàn)農(nóng)村經(jīng)濟(jì)的持續(xù)增長。產(chǎn)業(yè)結(jié)構(gòu)高級(jí)化尤其是加快發(fā)展高端、綠色、低碳產(chǎn)業(yè)是產(chǎn)業(yè)生態(tài)化的具體體現(xiàn)，同時(shí)產(chǎn)業(yè)結(jié)構(gòu)高級(jí)化是城鄉(xiāng)統(tǒng)籌縱深發(fā)展的必然趨勢(shì)。因此，珠三角農(nóng)村要充分發(fā)揮嶺南水鄉(xiāng)獨(dú)特的生態(tài)環(huán)境和自然資源優(yōu)勢(shì)，加快發(fā)展高效生態(tài)農(nóng)業(yè)和循環(huán)農(nóng)業(yè)。同時(shí)，要努力發(fā)展環(huán)境友好型工業(yè)，運(yùn)用高新技術(shù)改造和提升傳統(tǒng)行業(yè)，堅(jiān)決淘汰高耗能、高污染行業(yè)，真正實(shí)現(xiàn)清潔高效生產(chǎn)。更為關(guān)鍵的是要優(yōu)化第三產(chǎn)業(yè)成分和布局，積極扶持現(xiàn)代物流部門、生產(chǎn)性服務(wù)部門和公共服務(wù)部門的發(fā)展，優(yōu)先發(fā)展交通運(yùn)輸業(yè)、郵電通訊業(yè)、金融保險(xiǎn)業(yè)、房地產(chǎn)管理業(yè)和社會(huì)福利事業(yè)。

(3)加強(qiáng)技術(shù)創(chuàng)新，推動(dòng)產(chǎn)業(yè)升級(jí)換代。在“退二進(jìn)三”的產(chǎn)業(yè)轉(zhuǎn)移倒逼形勢(shì)下，珠三角農(nóng)村承接了城市能源消耗大、生產(chǎn)工藝落后、生產(chǎn)附加值低的勞動(dòng)密集型或資源密集型產(chǎn)業(yè)。由于農(nóng)村地區(qū)創(chuàng)新意識(shí)不強(qiáng)、創(chuàng)新環(huán)境不優(yōu)導(dǎo)致了農(nóng)村產(chǎn)業(yè)發(fā)展只能簡單復(fù)制城市的生產(chǎn)技術(shù)，重走“先污染，后治理”的發(fā)展老路。產(chǎn)業(yè)升級(jí)換代是珠三角農(nóng)村產(chǎn)業(yè)發(fā)展的出路，產(chǎn)業(yè)技術(shù)創(chuàng)新則是產(chǎn)業(yè)升級(jí)換代的重要驅(qū)動(dòng)力，更是促進(jìn)產(chǎn)業(yè)生態(tài)化、改善農(nóng)村水環(huán)境的有效手段。珠三角農(nóng)村產(chǎn)業(yè)發(fā)展應(yīng)該更多地依賴科學(xué)技術(shù)的創(chuàng)新和發(fā)展，注重高素質(zhì)人才的引進(jìn)和培養(yǎng)。一方面，要?jiǎng)?chuàng)新技術(shù)對(duì)傳統(tǒng)產(chǎn)業(yè)進(jìn)行改造和升級(jí)，比如通過參考丹麥卡倫堡工業(yè)園區(qū)等成功案例，建立共享資源和互換副產(chǎn)品的企業(yè)共生組合，形成工農(nóng)業(yè)相互促進(jìn)的物質(zhì)和能量循環(huán)，從而促進(jìn)傳統(tǒng)產(chǎn)業(yè)向生態(tài)化方向發(fā)展。另一方面，在依托農(nóng)村自然資源和生態(tài)環(huán)境等優(yōu)勢(shì)的基礎(chǔ)上，致力于經(jīng)濟(jì)增長的培育和開發(fā)，通過建立科技創(chuàng)新鼓勵(lì)機(jī)制、完善技術(shù)成果交易平臺(tái)、積極開展產(chǎn)學(xué)研合作等措施促進(jìn)科研成果和創(chuàng)新技術(shù)資本化、產(chǎn)業(yè)化，推進(jìn)生產(chǎn)技術(shù)由末端向中端、尖端前移，增強(qiáng)珠三角農(nóng)村產(chǎn)業(yè)發(fā)展后勁和競(jìng)爭(zhēng)力。

5.3 研究不足與展望

隨著人們對(duì)農(nóng)村發(fā)展及生態(tài)保護(hù)問題的日 益重視，農(nóng)村產(chǎn)業(yè)發(fā)展對(duì)水環(huán)境的影響研究逐漸成為環(huán)境經(jīng)濟(jì)學(xué)領(lǐng)域的關(guān)注熱點(diǎn)，具有一定的前沿性。由于本文分析架構(gòu)是從理論到實(shí)際、再從實(shí)際提升到理論的過程，在宏觀與微觀的轉(zhuǎn)換中可能存在不足。例如，為了提高研究區(qū)域的代表性，本文以東莞水鄉(xiāng)片9個(gè)鄉(xiāng)鎮(zhèn)作為一個(gè)整體來開展研究，對(duì)農(nóng)村產(chǎn)業(yè)發(fā)展及水污染變化的時(shí)間演化過程進(jìn)行詳細(xì)分析，因此缺乏對(duì)各鄉(xiāng)鎮(zhèn)的空間分異特征的深入研究。嘗試通過擴(kuò)大研究區(qū)域來探索農(nóng)村產(chǎn)業(yè)發(fā)展對(duì)水環(huán)境影響的空間分異特征是未來研究的方向。

□

致謝：感謝黃翅清同學(xué)在本文數(shù)據(jù)整理過程中給予的幫助。

[1] Krishna A K, Satyanarayanan M, Govil P K. Assessment of heavy metal pollution in water using multivariate statistical techniques in an industrial area: A case study from Patancheru, Medak District, Andhra Pradesh, India[J]. Journal of Hazardous Materials, 2009,167:366-373.

[2] Schaffner M, Bader H, Scheidegger R. Modeling non-point source pollution from rice farming in the Thachin River Basin[J]. Environment, Development and Sustainability, 2011,(13):403-422.

[3] 王 金, 李進(jìn)華, 陳 來, 等. 巢湖流域產(chǎn)業(yè)結(jié)構(gòu)與水污染程度的關(guān)系研究——基于灰色關(guān)聯(lián)分析法[J]. 資源開發(fā)與市場(chǎng), 2009,25(7):606-609.

[4] 趙小風(fēng), 黃賢金. 產(chǎn)業(yè)結(jié)構(gòu)演變的流域水環(huán)境響應(yīng)研究——以社瀆港流域?yàn)槔齕J]. 環(huán)境污染與防治, 2010,32(1):9-13.

[5] 王 磊, 張 磊, 段學(xué)軍, 等. 江蘇省太湖流域產(chǎn)業(yè)結(jié)構(gòu)的水環(huán)境污染效應(yīng)[J]. 生態(tài)學(xué)報(bào), 2011,31(22):6 832-6 844.

[6] Frosch R A，Gallopoulos N E. Straegies for Manufacturing[J].Scientific American,1989,261(3):94-102.

[7] 傅 沂, 隋廣軍. 生態(tài)管理的產(chǎn)業(yè)生態(tài)基礎(chǔ)研究[J]. 科學(xué)與科學(xué)技術(shù)管理, 2006,(4):86-91.

[8] 歐陽杰. 東莞統(tǒng)籌水鄉(xiāng)地區(qū)潛力旅游資源開發(fā)研究[J]. 旅游縱覽,2014,(2):177-178.

[9] 譚秀娟,鄭欽玉. 我國水資源生態(tài)足跡分析與預(yù)測(cè)[J]. 生態(tài)學(xué)報(bào),2009,(7):3 559-3 568.

[10] 陳百明. 區(qū)域土地可持續(xù)利用指標(biāo)體系框架的構(gòu)建與評(píng)價(jià)[J]. 地理科學(xué)進(jìn)展,2002,(3):204-215

[11] 劉小玲,甘建文. 珠三角地區(qū)水環(huán)境空間分異及其優(yōu)化對(duì)策研究[J]. 中國農(nóng)業(yè)資源與區(qū)劃,2015,(4):1-9.

[12] 趙海霞,曲福田,郭忠興. 環(huán)境污染影響因素的經(jīng)濟(jì)計(jì)量分析——以江蘇省為例[J]. 環(huán)境保護(hù),2006,(4):43-49.

[13] 單瑞峰,孫小銀. 環(huán)境污染區(qū)域差異及其影響因素灰色關(guān)聯(lián)法分析——以山東省為例[J]. 環(huán)境科學(xué)與管理,2008,(10):5-9.

[14] 張彥博,郭亞軍. FDI的環(huán)境效應(yīng)與我國引進(jìn)外資的環(huán)境保護(hù)政策[J]. 中國人口·資源與環(huán)境,2009,(4):7-12.