吉林省經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染水平演進(jìn)關(guān)系的實(shí)證研究

許清濤，高 標(biāo)，房 驕

（1.白城師范學(xué)院 地理科學(xué)學(xué)院，吉林 白城137000；2.東北師范大學(xué) 城市與環(huán)境科學(xué)學(xué)院，長(zhǎng)春130024）

環(huán)境問題的產(chǎn)生、發(fā)展與擴(kuò)大，與人類社會(huì)經(jīng)濟(jì)活動(dòng)密不可分，有效地計(jì)量評(píng)價(jià)經(jīng)濟(jì)發(fā)展和環(huán)境污染是當(dāng)前研究的熱點(diǎn)，環(huán)境庫(kù)茲涅茨曲線假設(shè)（EKC）是描述經(jīng)濟(jì)發(fā)展與環(huán)境污染演替關(guān)系的計(jì)量模型，是研究二者關(guān)系的有效方法［1－2］。1991年，環(huán)境經(jīng)濟(jì)學(xué)家Gross man等將庫(kù)茲涅茨曲線應(yīng)用于環(huán)境問題的研究，發(fā)現(xiàn)環(huán)境質(zhì)量先隨著經(jīng)濟(jì)的增長(zhǎng)而不斷惡化，當(dāng)經(jīng)濟(jì)發(fā)展到一定程度以后，環(huán)境質(zhì)量將隨著經(jīng)濟(jì)的增長(zhǎng)而不斷改善，即曲線呈現(xiàn)“倒U型”，進(jìn)而提出了環(huán)境庫(kù)茲涅茨曲線假設(shè)［3－4］。隨后，EKC在國(guó)內(nèi)外得到了廣泛的實(shí)證研究應(yīng)用［5－7］。隨著實(shí)證研究的不斷深入，一些學(xué)者發(fā)現(xiàn)，在不同的研究區(qū)域，EKC具有直線型、N型或U型等多種形式，并不只是傳統(tǒng)的“倒 U 型”［8－9］。本研究基于吉林省1981—2010年環(huán)境污染指標(biāo)與經(jīng)濟(jì)增長(zhǎng)指標(biāo)數(shù)據(jù)，利用PASW Statistics 18.0統(tǒng)計(jì)軟件進(jìn)行回歸模擬，得到各個(gè)環(huán)境污染指標(biāo)與經(jīng)濟(jì)增長(zhǎng)指標(biāo)的計(jì)量模型，研究分析其EKC曲線特征，結(jié)合灰色關(guān)聯(lián)分析方法探析EKC的影響因素，并提出經(jīng)濟(jì)與環(huán)境協(xié)調(diào)發(fā)展的建議和對(duì)策，以期為政府或相關(guān)部門制定政策與規(guī)劃提供科學(xué)、合理的依據(jù)。

1 研究區(qū)域概況

吉林省位于東北地區(qū)中部，地處北緯40°52′—46°18′，東經(jīng)121°38′—131°19′之間，總面積18.74萬(wàn)k m2，占全國(guó)的1.95%；目前人口為2 749.41萬(wàn)人，占全國(guó)總?cè)丝诘?.04%，近年來(lái)，伴隨著振興東北老工業(yè)基地戰(zhàn)略的大力有效實(shí)施，吉林省在經(jīng)濟(jì)規(guī)模和結(jié)構(gòu)上都有了較大的提高和優(yōu)化，2010年末，全省GDP總量達(dá)8 577.06億元，人均GDP達(dá)到31 306元，三次產(chǎn)業(yè)的結(jié)構(gòu)比例為12.2∶51.5∶36.3，對(duì)經(jīng)濟(jì)增長(zhǎng)的貢獻(xiàn)率分別為3.0%，67.1%和29.9%。吉林省在經(jīng)濟(jì)增長(zhǎng)的同時(shí)，也產(chǎn)生了大量的工業(yè)污染物排放，2010年，吉林省排放38 656萬(wàn)t工業(yè)廢水、8 240億m3工業(yè)廢氣、30.1萬(wàn)t工業(yè)SO2、5.3萬(wàn)t工業(yè)粉塵、13.1萬(wàn)t工業(yè)廢水COD，產(chǎn)生4 642萬(wàn)t工業(yè)固體廢棄物。

2 研究方法

2.1 指標(biāo)選取

選取的指標(biāo)包括環(huán)境污染指標(biāo)和經(jīng)濟(jì)增長(zhǎng)指標(biāo)。針對(duì)環(huán)境污染指標(biāo)，從流量指標(biāo)和存量指標(biāo)兩個(gè)方面考慮［10］。其中，流量指標(biāo)包括工業(yè)廢水排放量、工業(yè)廢氣排放量、工業(yè)固廢物產(chǎn)生量；存量指標(biāo)包括工業(yè)二氧化硫排放量、工業(yè)粉塵排放量和工業(yè)廢水COD排放量。關(guān)于經(jīng)濟(jì)增長(zhǎng)指標(biāo)，由于人均GDP能較好地說(shuō)明經(jīng)濟(jì)增長(zhǎng)水平且有利于不同研究區(qū)域的相互比較，故本研究選取人均GDP為經(jīng)濟(jì)增長(zhǎng)指標(biāo)。本研究收集整理了吉林省1981—2010年的工業(yè)污染物排放及人均GDP相關(guān)數(shù)據(jù)，具體見表1。

2.2 計(jì)量模型構(gòu)建

應(yīng)用PASW Statistics 18.0統(tǒng)計(jì)軟件，以經(jīng)濟(jì)增長(zhǎng)指標(biāo)（人均GDP）為自變量（X），以各環(huán)境污染指標(biāo)為因變量（Y），分別選用線性、二次項(xiàng)、對(duì)數(shù)、立方、指數(shù)分布及Logistic等多種函數(shù)對(duì)表1中的數(shù)據(jù)進(jìn)行回歸模擬，建立工業(yè)污染物排放量與人均GDP的計(jì)量模型，根據(jù)決定系數(shù)和模型檢驗(yàn)參數(shù)進(jìn)行優(yōu)選。

2.3 灰色關(guān)聯(lián)分析

灰色關(guān)聯(lián)分析是根據(jù)因素間動(dòng)態(tài)過程的相似性或相異程度來(lái)衡量因素間發(fā)展態(tài)勢(shì)的一種方法［11－12］。其計(jì)算公式為：

式中：r0i——關(guān)聯(lián)度；N——每個(gè)相關(guān)因素序列的數(shù)據(jù)個(gè)數(shù)；L0i（k）——相關(guān)因素序列k時(shí)刻的關(guān)聯(lián)系數(shù)；Δ0i（k）——相關(guān)因素序列k時(shí)刻的絕對(duì)差；Δmax，Δmin——所有相關(guān)因素序列各個(gè)時(shí)刻絕對(duì)差中的最大值與最小值；ρ——分辨系數(shù)，本研究中取ρ＝0.5。

3 結(jié)果與分析

3.1 計(jì)量模型的選取及分析

3.1.1 計(jì)量模型的選取與設(shè)計(jì) 應(yīng)用PASW Statistics 18.0軟件，針對(duì)多種函數(shù)類型進(jìn)行回歸，根據(jù)決定系數(shù)和檢驗(yàn)參數(shù)優(yōu)選出各個(gè)工業(yè)污染物排放量與人均GDP的計(jì)量模型，除了工業(yè)粉塵排放量與人均GDP符合對(duì)數(shù)函數(shù)關(guān)系以外，其他環(huán)境污染指標(biāo)與人均GDP均符合三次函數(shù)的計(jì)量模型，具體結(jié)果見表2。

3.1.2 EKC曲線特征分析 根據(jù)優(yōu)選的計(jì)量模型，得到各個(gè)環(huán)境指標(biāo)的EKC，如圖1所示。從圖1可以看出，6個(gè)環(huán)境指標(biāo)的曲線形狀各不相同，各具特點(diǎn)。

表1 吉林省1981－2010年環(huán)境污染與經(jīng)濟(jì)增長(zhǎng)指標(biāo)

表2 吉林省人均GDP與污染物排放量的計(jì)量模型

（1）從工業(yè)廢水排放量的EKC可以看到，整個(gè)曲線為三次函數(shù)，整體形狀不符合理想的“倒U型”，而是呈現(xiàn)一個(gè)“正U＋倒U型”，存在兩個(gè)拐點(diǎn)，取計(jì)量模型函數(shù)的導(dǎo)數(shù)Y′＝0，求解得到X1＝11975.949，X2＝25540.427，分別對(duì)應(yīng)的時(shí)間為2004—2005年之間、2008—2009年之間。取計(jì)量模型函數(shù)二階導(dǎo)數(shù)Y″＝0，求解得到X3＝18578.188，對(duì)應(yīng)時(shí)間為2006—2007年之間。以上說(shuō)明：吉林省在1981—2001年，工業(yè)廢水排放量隨人均GDP的增長(zhǎng)而下降，2004后到達(dá)拐點(diǎn)，從2005開始到2010年，整個(gè)曲線符合EKC的先升后降的“倒U型”，2006—2007年之間，工業(yè)廢水排放量由加速增加轉(zhuǎn)變?yōu)闇p速增加，并且在2008—2009年之間到達(dá)拐點(diǎn)。

（2）工業(yè)廢氣排放量的EKC表明，1981—2010年，工業(yè)廢氣排放量一直隨人均GDP的增長(zhǎng)而增長(zhǎng)，針對(duì)計(jì)量模型而言，取其導(dǎo)數(shù)Y′＝0，無(wú)實(shí)數(shù)解，取其二階導(dǎo)數(shù)Y″＝0，得：X3＝4777.382，對(duì)應(yīng)時(shí)間范圍是1995—1996年之間。說(shuō)明雖然在1996年以后，工業(yè)廢氣排放量已經(jīng)由加速增長(zhǎng)轉(zhuǎn)變?yōu)闇p速增長(zhǎng)，但整體持續(xù)增長(zhǎng)的趨勢(shì)并未發(fā)生改變，不存在拐點(diǎn)。如果吉林省未來(lái)不對(duì)現(xiàn)有的經(jīng)濟(jì)發(fā)展方式，特別是能源利用方式進(jìn)行優(yōu)化和改變，工業(yè)廢氣排放量的狀況將進(jìn)一步惡化。

（3）從工業(yè)固體廢棄物產(chǎn)生量的EKC可以看到，吉林省1981—2010年的工業(yè)固體廢棄物產(chǎn)生量基本呈現(xiàn)逐漸上升的態(tài)勢(shì)，大致處于“倒U型”的左半部分，存在一個(gè)拐點(diǎn)，取計(jì)量模型函數(shù)的導(dǎo)數(shù)Y′＝0，求解得到X1＝－1456.279（舍去）、X2＝52723.144，根據(jù)表1中人均GDP變化數(shù)據(jù)，應(yīng)用PASW Statistics 18.0對(duì)人均GDP與年份（以1981年為1，1982年為2，其余年份以此類推）進(jìn)行線性回歸，得到計(jì)量模型為：Y＝251.608exp（0.160X）＋661.259（R2＝0.994、F＝4917.602），在模型中令Y＝52723.144，得X＝33.327，取X＝34，說(shuō)明吉林省如果保持現(xiàn)有的發(fā)展態(tài)勢(shì)，其工業(yè)固廢產(chǎn)生量的EKC將在2014年到達(dá)拐點(diǎn)。取計(jì)量模型函數(shù)二階導(dǎo)數(shù)Y″＝0，求解得到X3＝25633.432，對(duì)應(yīng)時(shí)間為2008—2009年之間，說(shuō)明這一時(shí)間段內(nèi)，工業(yè)固廢排放量由加速增長(zhǎng)轉(zhuǎn)變?yōu)闇p速增長(zhǎng)。

（4）工業(yè)粉塵排放量的EKC呈對(duì)數(shù)函數(shù)，處于逐漸減小的趨勢(shì)，不存在拐點(diǎn)；工業(yè)SO2排放量的EKC呈“正U＋倒U型”，存在兩個(gè)拐點(diǎn)，取導(dǎo)數(shù)Y′＝0，求解得到X1＝6762.671、X2＝23294.039，分別對(duì)應(yīng)的時(shí)間為1999—2000年之間、2007—2008年之間。取計(jì)量模型函數(shù)二階導(dǎo)數(shù)Y″＝0，求解得到X3＝15028.355，對(duì)應(yīng)時(shí)間為2005—2006年之間；工業(yè)廢水COD排放量的EKC呈現(xiàn)不是很明顯的“正U＋倒U型”，存在兩個(gè)拐點(diǎn)，取計(jì)量模型函數(shù)的導(dǎo)數(shù)Y′＝0，求解得到X1＝13907.226、X2＝26386.922，分別對(duì)應(yīng)的時(shí)間為2005—2006年之間、2008—2009年之間。取計(jì)量模型函數(shù)二階導(dǎo)數(shù)Y″＝0，求解得到X3＝20147.074，對(duì)應(yīng)時(shí)間為2007—2008年之間，由此說(shuō)明：吉林省的工業(yè)SO2和工業(yè)廢水COD排放量的EKC在人均GDP達(dá)到一定數(shù)值后，均存在拐點(diǎn)，符合 EKC 理論［13］。

圖1 吉林省1981－2010年人均GDP與工業(yè)污染物排放量擬合曲線

3.2 吉林省EKC影響因素分析

為進(jìn)一步探析吉林省EKC的影響因素，根據(jù)表3中吉林省1981—2010年的數(shù)據(jù)，應(yīng)用模型對(duì)吉林省環(huán)境污染指標(biāo)與各個(gè)影響因素進(jìn)行灰色關(guān)聯(lián)分析，結(jié)果見表4。

表3 吉林省1981－2010年環(huán)境變化的影響因素

表4 吉林省環(huán)境污染指標(biāo)與影響因素關(guān)聯(lián)度

表4表明，吉林省環(huán)境污染的各個(gè)指標(biāo)與環(huán)境變化的影響因素有著密切的聯(lián)系，主要的影響因素有工業(yè)比率、第三產(chǎn)業(yè)比率、單位GDP能耗、環(huán)境科研人數(shù)4個(gè)方面，其余因素影響較小，在此不予詳細(xì)討論。工業(yè)廢水排放量的影響因素依次為工業(yè)比率（0.996）、單 位 GDP 能 耗 （0.989）、第 三 產(chǎn) 業(yè) 比 率（0.976），環(huán)境科研人數(shù)（0.897）；工業(yè)廢氣排放量的影響因素依次為第三產(chǎn)業(yè)比率（0.989）、工業(yè)比率（0.972）、單位 GDP能耗（0.957）、環(huán)境科研人員數(shù)（0.935）；工業(yè)固廢產(chǎn)生量的影響因素依次為第三產(chǎn)業(yè)比率（0.987）、工業(yè)比率（0.984）、單位 GDP能耗（0.971）、環(huán)境科研人員數(shù)（0.907）；工業(yè)粉塵排放量的影響因素依次為單位GDP能耗（0.996）、工業(yè)比率（0.990）、第三產(chǎn)業(yè)比率（0.968）、環(huán)境科研人員數(shù)（0.890）；工業(yè)SO2排放量的影響因素依次為工業(yè)比率（0.997）、單位 GDP能耗（0.988）、第三產(chǎn)業(yè)比率（0.975）、環(huán)境科研人員數(shù)（0.895）；工業(yè)廢水 COD排放量的影響因素依次為工業(yè)比率（0.982）、單位GDP能耗（0.980）、第三產(chǎn)業(yè)比率（0.955）、環(huán)境科研人員數(shù)（0.935）。

綜上可知，吉林省的環(huán)境污染指標(biāo)的變化主要受產(chǎn)業(yè)結(jié)構(gòu)與技術(shù)進(jìn)步的影響。從產(chǎn)業(yè)結(jié)構(gòu)方面來(lái)說(shuō)，工業(yè)比率的增加會(huì)加大環(huán)境污染物的排放量，而第三產(chǎn)業(yè)比率的增加會(huì)在很大程度上減少環(huán)境污染物的排放量［14］。從表1、表3可知，1981—2010年，吉林省的工業(yè)比率與第三產(chǎn)業(yè)比率幾乎是互為反方向發(fā)展的，而環(huán)境污染物排放量隨著工業(yè)比率的增大和第三產(chǎn)業(yè)比率的減小而呈現(xiàn)增長(zhǎng)趨勢(shì)，反之亦成立。自新世紀(jì)，尤其是振興東北老工業(yè)基地戰(zhàn)略實(shí)施以來(lái)，吉林省的經(jīng)濟(jì)得到了快速發(fā)展，工業(yè)與第三產(chǎn)業(yè)的生產(chǎn)總值一直處于迅速增加的狀態(tài)，2010年，工業(yè)產(chǎn)值和第三產(chǎn)業(yè)產(chǎn)值分別達(dá)到3 929.31億元、3 111.12億元，但從表3可見，第三產(chǎn)業(yè)比率卻從2002年的40.83%降到2010年的35.89%，同時(shí)期工業(yè)比率則由34.21%升高到45.33%，2010年的國(guó)家工業(yè)比率和第三產(chǎn)業(yè)比率分別為40.0%和43.2%。因此，在產(chǎn)業(yè)結(jié)構(gòu)方面，吉林省具有非常大的優(yōu)化空間；從技術(shù)進(jìn)步方面來(lái)說(shuō)，1981—2010年，吉林省的單位GDP能耗一直處于下降的態(tài)勢(shì)，2000年以后下降趨勢(shì)變緩，2010年達(dá)到0.96 t標(biāo)準(zhǔn)煤／萬(wàn)元，小于當(dāng)年國(guó)家0.81 t標(biāo)準(zhǔn)煤／萬(wàn)元的單位GDP能耗，存在進(jìn)一步減小的可能；同時(shí)環(huán)境科研人數(shù)的增加也能有效減少環(huán)境污染物的排放量。吉林省在經(jīng)濟(jì)快速發(fā)展的同時(shí)，也帶來(lái)了較大的資源能源消耗及污染物排放量，吉林省要想進(jìn)一步改善環(huán)境質(zhì)量，控制環(huán)境污染水平，必須合理布局產(chǎn)業(yè)，大力發(fā)展第三產(chǎn)業(yè)與高技術(shù)產(chǎn)業(yè)等低消耗、低排放產(chǎn)業(yè)；加強(qiáng)科技創(chuàng)新，打破技術(shù)瓶頸，進(jìn)一步降低單位GDP能耗；提倡發(fā)展循環(huán)經(jīng)濟(jì)技術(shù)，有效控制工業(yè)發(fā)展的過度廢棄物排放。

4 結(jié)論與建議

通過對(duì)吉林省1981—2010年的環(huán)境污染指標(biāo)EKC的變化分析，以及環(huán)境污染指標(biāo)與各影響因素的灰色關(guān)聯(lián)分析，可以得到以下結(jié)論：

（1）吉林省的環(huán)境污染指標(biāo)EKC形狀各不相同，不完全符合典型的“倒U型”特征。其中，工業(yè)廢水排放量的EKC呈現(xiàn)“正U＋倒U型”，在2009年達(dá)到拐點(diǎn)；工業(yè)廢氣排放量的EKC呈連續(xù)增長(zhǎng)態(tài)勢(shì)，無(wú)拐點(diǎn)；工業(yè)固體廢棄物產(chǎn)生量的EKC大致處于“倒U型”的左半部分，拐點(diǎn)在現(xiàn)有數(shù)據(jù)之外，通過人均GDP計(jì)量模型預(yù)測(cè)2014年到達(dá)拐點(diǎn)；工業(yè)粉塵排放量的EKC呈對(duì)數(shù)函數(shù)，呈逐漸減小的趨勢(shì)，不存在拐點(diǎn)；工業(yè)SO2排放量的EKC呈“正U＋倒U型”，2008年達(dá)到拐點(diǎn)；工業(yè)廢水COD排放量的EKC呈現(xiàn)不明顯的“正U＋倒U型”，2009年達(dá)到拐點(diǎn)。伴隨著經(jīng)濟(jì)增長(zhǎng)，除了工業(yè)廢氣排放量呈污染連續(xù)加劇以外，吉林省總的環(huán)境污染處于不斷改善的階段。從這種意義上說(shuō)，環(huán)境庫(kù)茲涅茨曲線假設(shè)成立，即經(jīng)濟(jì)增長(zhǎng)初期會(huì)造成環(huán)境污染水平加劇，但當(dāng)經(jīng)濟(jì)增長(zhǎng)達(dá)到一定的程度以后，環(huán)境污染水平將隨經(jīng)濟(jì)增長(zhǎng)而不斷改善。

（2）吉林省環(huán)境污染指標(biāo)與環(huán)境變化的影響因素有著密切的聯(lián)系，主要的影響因素有工業(yè)比率、第三產(chǎn)業(yè)比率、單位GDP能耗、環(huán)境科研人數(shù)4個(gè)方面。

雖然吉林省的部分環(huán)境污染指標(biāo)的EKC已經(jīng)出現(xiàn)拐點(diǎn)，但像工業(yè)廢水、工業(yè)廢氣和工業(yè)固體廢棄物排放等環(huán)境污染指標(biāo)排放的絕對(duì)量卻一直在增加，為有效遏制吉林省環(huán)境污染物的排放，從根本上解決存在的環(huán)境問題，未來(lái)吉林省應(yīng)該在以下方面采取對(duì)策：優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，合理布局工業(yè)，減少工業(yè)比率，大力發(fā)展第三產(chǎn)業(yè)與高技術(shù)產(chǎn)業(yè)等低消耗、低排放產(chǎn)業(yè)；轉(zhuǎn)換經(jīng)濟(jì)增長(zhǎng)模式，大力發(fā)展循環(huán)經(jīng)濟(jì)，減少能源資源消耗與廢棄物排放；加大科技創(chuàng)新力度，加快科技轉(zhuǎn)化為生產(chǎn)力的速度，加速降低單位GDP能耗；有效加強(qiáng)公民環(huán)境意識(shí)，提高環(huán)境科研人數(shù)及素質(zhì)；加大工業(yè)污染治理投資，豐富資金來(lái)源，監(jiān)管資金有效投入；完善相關(guān)法律法規(guī)，有效實(shí)施環(huán)境保護(hù)政策。

［1］ 宋萍，洪偉，吳承禎，等.福建省污染物排放的環(huán)境Kuznets曲線研究［J］.江西農(nóng)業(yè)大學(xué)學(xué)報(bào)，2009，31（3）：574-578.

［2］ 楊麗霞.浙江省經(jīng)濟(jì)增長(zhǎng)與環(huán)境質(zhì)量演進(jìn)關(guān)系的實(shí)證研究［J］.地理與地理信息科學(xué)，2010，26（1）：78-80.

［3］ Ckossmai G M，Kr ueger A B.Econo mic growth and the environ ment［J］.Quarterly Jour nal of Economics，1995，110（2）：353-377.

［4］ 姚煥玫，唐國(guó)滔，莫?jiǎng)?chuàng)榮，等.基于環(huán)境庫(kù)茲涅茨曲線的經(jīng)濟(jì)增長(zhǎng)與環(huán)境質(zhì)量實(shí)證研究［J］.環(huán)境污染與防治，2010，32（11）：74-77.

［5］ 陳強(qiáng)強(qiáng)，竇學(xué)誠(chéng).甘肅省環(huán)境庫(kù)茲涅茨曲線估計(jì)及其驅(qū)動(dòng)因子分析［J］.干旱區(qū)地理，2011，34（5）：866-873.

［6］ 蘇秋實(shí)，王立本.山東省工業(yè)“三廢”排放的庫(kù)茲涅茨特征研究［J］.水土保持研究，2008，15（3）：258-259.

［7］ 蘇力葉·木沙江，孜比布拉·司馬義，周玄德.吐魯番地區(qū)經(jīng)濟(jì)增長(zhǎng)與環(huán)境質(zhì)量水平的計(jì)量模型研究［J］.水土保持研究，2011，18（3）：55-60.

［8］ Or ubu C O，Omotor D G.Envir on mental quality and economic growth：searching for environmental Kuznets curves for air and water pollutants in Africa［J］.Energy Policy，2011，39（7）：4178-4188.

［9］ 李政，王坎，任津，等.湖北省主要工業(yè)污染物排放的環(huán)境庫(kù)茲涅茨特征研究［J］.環(huán)境科學(xué)與技術(shù)，2011，34（12）：55-59.

［10］ 黃一綏，邱健斌，佘晨興.福建省經(jīng)濟(jì)發(fā)展與工業(yè)污染水平計(jì)量模型研究［J］.環(huán)境污染與防治，2010，32（3）：90-93.

［11］ 于文良，王伯鐸，吳良興.陜西省可持續(xù)發(fā)展能力變化趨勢(shì)和影響因素分析［J］.干旱區(qū)資源與環(huán)境，2009，23（6）：13-18.

［12］ 宮繼萍，潘竟虎，石培基.基于生態(tài)足跡和灰色關(guān)聯(lián)度的甘肅省可持續(xù)發(fā)展研究［J］.水土保持研究，2011，18（2）：198-201.

［13］ 夏自蘭，趙小風(fēng)，王繼軍.江蘇省環(huán)境庫(kù)茲涅茨曲線特征及其成因分析［J］.水土保持研究，2010，17（1）：198-201.

［14］ 柯文嵐，沙景華，閆晶晶.山西省環(huán)境庫(kù)茲涅茨曲線特征及其影響因素分析［J］.中國(guó)人口·資源與環(huán)境，2011，21（12）：389-392.