盧小蘭

摘要依據(jù)PRED框架，建立資源環(huán)境壓力和資源環(huán)境承載力評(píng)價(jià)指標(biāo)體系，利用“縱橫向拉開檔次法”評(píng)價(jià)2006-2011年間省域資源環(huán)境承載情況，基于空間統(tǒng)計(jì)方法分析資源環(huán)境壓力、承載力和承載率的空間相關(guān)性。結(jié)果表明：研究期間，省域資源環(huán)境壓力由西向東、由北到南依次遞減，資源環(huán)境承載力和承載率呈現(xiàn)出由東部沿海向西北部梯狀下降模式。此外，資源環(huán)境壓力、承載力和承載率均存在較強(qiáng)的全局正相關(guān)，呈現(xiàn)出“中心—外圍”格局。

關(guān)鍵詞資源環(huán)境壓力資源環(huán)境承載力縱橫向拉開檔次法Moran'sI指數(shù)空間統(tǒng)計(jì)

中圖分類號(hào)：F062.1 文獻(xiàn)標(biāo)識(shí)碼：Ａ

資源環(huán)境承載力概念是衡量人類經(jīng)濟(jì)社會(huì)活動(dòng)與資源環(huán)境相互關(guān)系的科學(xué)指標(biāo)，是人類可持續(xù)發(fā)展的度量和管理決策的重要依據(jù)。國(guó)內(nèi)外關(guān)于資源環(huán)境壓力的定量評(píng)價(jià)提出了一些指標(biāo)和方法，如生態(tài)足跡法、STIRPAT模型法等。關(guān)于資源環(huán)境承載力的研究主要涉及承載體、承載對(duì)象和承載力等三要素。毛漢英，余丹林等（2001a，2001b） 建立由壓力類、承壓類和區(qū)際交流共27項(xiàng)指標(biāo)構(gòu)成的指標(biāo)體系，綜合利用狀態(tài)空間法、系統(tǒng)動(dòng)力學(xué)、回歸預(yù)測(cè)、灰色預(yù)測(cè)、差分遞歸預(yù)測(cè)等多種方法，構(gòu)建了環(huán)渤海地區(qū)區(qū)域承載力的SD模型，利用DYNAMO軟件進(jìn)行預(yù)測(cè)和模擬。葉明霞，羅國(guó)云（2009）利用矢量法對(duì)長(zhǎng)江中上游地區(qū)資源環(huán)境承載力進(jìn)行實(shí)證分析。唐凱等（2012）采用集對(duì)分析法對(duì)長(zhǎng)株潭城市群的資源環(huán)境綜合承載力進(jìn)行的測(cè)算。

通過(guò)對(duì)相關(guān)文獻(xiàn)的閱讀，目前資源環(huán)境承載情況評(píng)價(jià)有以下特點(diǎn)：權(quán)重設(shè)定缺乏客觀性；評(píng)價(jià)缺乏連續(xù)性；資源環(huán)境的空間相關(guān)性未考慮。

由于區(qū)域是由人口、資源、環(huán)境、發(fā)展（population，resources，environment，development，PRED）等4個(gè)子系統(tǒng)組成的巨系統(tǒng)，本文依據(jù)PRED框架，建立資源環(huán)境壓力和資源環(huán)境承載力評(píng)價(jià)指標(biāo)體系，利用“縱橫向拉開檔次法”對(duì)2006-2011年間中國(guó)30個(gè)省域資源環(huán)境壓力、資源環(huán)境承載力評(píng)價(jià)，進(jìn)一步計(jì)算資源環(huán)境承載力等的Moran'sI指數(shù)，對(duì)資源環(huán)境承載情況進(jìn)行空間統(tǒng)計(jì)研究。

一、區(qū)域資源環(huán)境壓力、資源環(huán)境承載力和資源環(huán)境承載率的科學(xué)內(nèi)涵

區(qū)域資源環(huán)境壓力是指由于資源消耗巨大和生態(tài)環(huán)境破壞嚴(yán)重等原因，同時(shí)人類為謀求自身發(fā)展等原因而對(duì)承載體施加的壓力。

區(qū)域資源環(huán)境承載力是指，“在一定的時(shí)期和一定尺度區(qū)域范圍內(nèi)，在確保資源合理開發(fā)利用和生態(tài)環(huán)境良性循環(huán)的條件下，資源環(huán)境能夠承載的人口數(shù)量及相應(yīng)的經(jīng)濟(jì)社會(huì)總量的能力”。

區(qū)域資源環(huán)境承載率是指，以一個(gè)比研究區(qū)域更大的參考區(qū)域作為對(duì)比區(qū)域，計(jì)算出相對(duì)于對(duì)比區(qū)域的各研究區(qū)域的相對(duì)資源環(huán)境承載力，即資源環(huán)境承載率。

二、區(qū)域資源環(huán)境壓力、資源環(huán)境承載力評(píng)價(jià)指標(biāo)體系

構(gòu)建資源環(huán)境壓力和資源環(huán)境承載力指標(biāo)體系時(shí)，應(yīng)體現(xiàn)科學(xué)性、代表性、可行性和綜合性等原則。指標(biāo)體系如下表1。

以上指標(biāo)數(shù)據(jù)來(lái)源于中國(guó)統(tǒng)計(jì)信息網(wǎng)和中經(jīng)網(wǎng)統(tǒng)計(jì)數(shù)據(jù)庫(kù)等。由于西藏自治區(qū)的部分統(tǒng)計(jì)數(shù)據(jù)缺失，評(píng)價(jià)對(duì)象未包括西藏。

三、研究方法

（一）縱橫向拉開檔次評(píng)價(jià)法。

為了充分利用時(shí)序立體數(shù)據(jù)信息，指標(biāo)權(quán)重可由縱橫向拉開檔次法確定。具體步驟如下：

第1步，數(shù)據(jù)無(wú)量綱化處理。對(duì)原始數(shù)據(jù)進(jìn)行無(wú)量綱化處理，得到時(shí)序立體數(shù)據(jù)及無(wú)量綱數(shù)據(jù)矩陣Ak（k=1，2，…N）。

第2步，確定權(quán)重。由無(wú)量綱數(shù)據(jù)矩陣Ak（k=1，2，…N）求對(duì)應(yīng)的對(duì)稱矩陣Hk=AAk（k=1，2，…N），求和得到H=，計(jì)算H的最大特征值對(duì)應(yīng)的特征向量或求解最優(yōu)化問(wèn)題，對(duì)特征向量或規(guī)劃問(wèn)題最優(yōu)解進(jìn)行歸一化處理得到的向量即為權(quán)重w。

第3步，計(jì)算綜合評(píng)價(jià)值。由yi（tk）=，計(jì)算各樣本在各時(shí)序點(diǎn)上的綜合評(píng)價(jià)值。

（二）空間統(tǒng)計(jì)的理論與方法。

空間統(tǒng)計(jì)研究的空間效應(yīng)包括空間相關(guān)性和空間差異性。空間差異性是指地理空間上存在發(fā)達(dá)地區(qū)和落后地區(qū)、核心地區(qū)和邊緣地區(qū)等經(jīng)濟(jì)地理結(jié)構(gòu)，從而導(dǎo)致資源環(huán)境壓力和資源環(huán)境承載力存在較大的空間差異性。此外，由于相鄰地區(qū)資源環(huán)境稟賦的相似性和空間上的觀測(cè)值缺乏獨(dú)立性，省級(jí)區(qū)域資源環(huán)境壓力和資源環(huán)境承載力同時(shí)存在空間相關(guān)性，可用Moran'sI指數(shù)來(lái)測(cè)算和檢驗(yàn)空間相關(guān)性。

全域的Moran'sI定義如下：

其中Yi，表示i第個(gè)區(qū)域的觀測(cè)值，wij為采用鄰近標(biāo)準(zhǔn)的空間權(quán)值矩陣。

由于全局Moran'sI估計(jì)的是抵消后的剩余值，無(wú)法揭示每一區(qū)域的局域空間關(guān)聯(lián)效應(yīng)，可用局域Moran'sI散點(diǎn)圖直觀地描述省級(jí)區(qū)域資源環(huán)境壓力和資源環(huán)境綜承載力的局域相關(guān)性。

四、省域資源環(huán)境承載情況的實(shí)證分析

（一）省域資源環(huán)境壓力、承載力評(píng)價(jià)及承載率計(jì)算。

由資源環(huán)境壓力類和資源環(huán)境承載力類總和對(duì)稱矩陣計(jì)算得到相對(duì)權(quán)重wy和wc：

wy=（0.090，0.189，0.153，0.185，0.153，0.094，0.135）

wc=（-0.309，-0.198，-0.159，0.061，0.297，0.061，0.330，0.239，-0.196，0.365，0.288，0.406，0.398）

將標(biāo)準(zhǔn)化處理后的和上述wy，wc分別代入=，得到省級(jí)區(qū)域資源環(huán)境壓力、承載力和承載率，如下表2（2007-2010年數(shù)據(jù)略）。

（二）省域資源環(huán)境承載情況分析。

1、省級(jí)區(qū)域資源環(huán)境壓力分析。

利用聚類分析將區(qū)域資源環(huán)境壓力分為3類。2011年省域資源環(huán)境壓力空間分布如圖1所示。資源環(huán)境壓力嚴(yán)重區(qū)域包括山西、內(nèi)蒙古、廣西、貴州、青海、寧夏、新疆等西部及西南部地區(qū)。資源環(huán)境壓力居中型區(qū)域有河北、遼寧、吉林、黑龍江、安徽、江西、河南、湖北、湖南、海南、重慶、四川、云南、陜西和甘肅等地區(qū)。這些區(qū)域主要分布于中部及東北部地區(qū)。資源環(huán)境壓力較輕區(qū)域有北京、天津、上海、江蘇、浙江、福建、山東和廣東，主要分布于東部沿海地區(qū)。

2、資源環(huán)境承載力分析。

2011年省域資源環(huán)境承載力空間分布如圖2所示。資源環(huán)境承載力強(qiáng)區(qū)域主要有北京、天津、上海、江蘇、浙江等東南部沿海地區(qū)。其中，北京和上海2006-2011年間一直穩(wěn)居資源環(huán)境承載力的超強(qiáng)位置?？傮w來(lái)說(shuō)，省域資源環(huán)境承載力的第一集團(tuán)區(qū)域已形成。資源環(huán)境承載力居中型區(qū)域主要有福建、山東、廣東、廣西、河北、遼寧、安徽、江西、河南、湖北、湖南、海南、重慶等地區(qū)。資源環(huán)境承載力弱的區(qū)域有山西、內(nèi)蒙古、貴州、青海、寧夏、新疆、吉林、黑龍江、四川、云南、陜西、甘肅。

3、資源環(huán)境承載率分析。

2011年省域資源環(huán)境承載率空間分布如圖3所示。資源環(huán)境承載能力富余區(qū)域主要有北京、天津、上海、江蘇、浙江和廣東等東南部沿海地區(qū)。資源環(huán)境承載力適中的區(qū)域主要有河北、遼寧、吉林、黑龍江、安徽、福建、江西、山東、河南、湖北、湖南、海南、重慶、四川、云南和陜西等地區(qū)。資源環(huán)境承載力超載的區(qū)域主要有山西、內(nèi)蒙古、廣西、貴州、甘肅、青海、寧夏、新疆等8個(gè)省域。

五、省域資源環(huán)境壓力、承載力和承載率的空間相關(guān)性分析

（一）空間權(quán)值矩陣選擇及全局空間相關(guān)分析。

分別采用rook鄰近、k-nearest鄰近和d距離鄰近的空間權(quán)值矩陣計(jì)算Moran'sI值，由于選取空間權(quán)值矩陣Wr1效果顯著，下文采用Wr1空間權(quán)值矩陣分析省域資源環(huán)境壓力、綜合承載力和承載率的空間相關(guān)關(guān)系。

全局空間自相關(guān)計(jì)算結(jié)果表明，資源環(huán)境壓力（承載力和承載率）的空間分布都是集聚的，而且這種集聚程度均呈進(jìn)一步加深趨勢(shì)。

表3資源環(huán)境壓力、承載力及承載率的全局Moran'sI值

（二）局域空間相關(guān)分析。

1、資源環(huán)境壓力的局域空間相關(guān)分析。

2006年及2011年省域資源環(huán)境承載力的Moran'sI散點(diǎn)圖（圖略）顯示，2006年云南、新疆、陜西、山西、青海、寧夏、內(nèi)蒙古、貴州、廣西和甘肅等10省域位于第一象限，其資源環(huán)境壓力呈現(xiàn)出HH型集聚態(tài)勢(shì)，而浙江、山東、湖北、河南、北京、天津、海南、福建、安徽、上海、江蘇和河北等12個(gè)省域的資源環(huán)境壓力表現(xiàn)為L(zhǎng)L型集聚，LH型離群地區(qū)包括四川、遼寧、吉林、湖南、黑龍江、重慶和廣東等7個(gè)省域，江西省屬于HL型離群地區(qū)。與2006年相比，2011年有5個(gè)省區(qū)的資源環(huán)境壓力的Moran散點(diǎn)圖所在象限發(fā)生改變，其他25個(gè)省區(qū)所在象限均沒(méi)有發(fā)生改變。

2、資源環(huán)境承載力的局域空間相關(guān)分析。

2006年及2011年省域資源環(huán)境承載力的Moran'sI散點(diǎn)圖表明，2006年度浙江、山東、湖南、北京、天津、海南、福建、安徽、上海、江蘇、廣東和河北等12省市位于第一象限，其資源環(huán)境承載力呈現(xiàn)出HH型集聚態(tài)勢(shì)，而云南、新疆、四川、陜西、山西、青海、寧夏、內(nèi)蒙古、遼寧、吉林、湖北、黑龍江、河南、貴州、廣西和甘肅等16個(gè)省域的資源環(huán)境承載力表現(xiàn)為L(zhǎng)L型集聚，江西為L(zhǎng)H型離群地區(qū)，重慶為HL型離群地區(qū)。2011年資源環(huán)境承載力除極少數(shù)地區(qū)（湖南和海南從HH型轉(zhuǎn)變?yōu)長(zhǎng)L型，河北由HH型轉(zhuǎn)變?yōu)長(zhǎng)H型）局域關(guān)系改變外，其它27個(gè)地區(qū)的空間格局均未改變。

3、資源環(huán)境承載率的局域空間相關(guān)分析。

2006年及2011年省域資源環(huán)境承載率的Moran'sI散點(diǎn)圖顯示，2006年浙江、山東、北京、天津、海南、福建、安徽、上海、江蘇等9省市位于第一象限，其資源環(huán)境承載率呈現(xiàn)出HH型集聚態(tài)勢(shì)，而云南、新疆、四川、陜西、山西、青海、寧夏、內(nèi)蒙古、吉林、湖南、黑龍江、河南、貴州、廣西、甘肅等15個(gè)省域的資源環(huán)境承載率表現(xiàn)為L(zhǎng)L型集聚，LH型空間離群地區(qū)包括江西和河北等2個(gè)省域，HL型空間離群地區(qū)包括遼寧、湖北、重慶和廣東等4個(gè)省域。相對(duì)于2006年，2011年承載率除極少數(shù)地區(qū)空間局域關(guān)系改變外，其它26個(gè)地區(qū)的局域空間格局均未改變。

六、結(jié)論

文章基于“縱橫向拉開檔次法”對(duì)省域資源環(huán)境壓力和承載力進(jìn)行評(píng)價(jià)，對(duì)資源環(huán)境承載情況進(jìn)行空間相關(guān)性分析，借助ArcView地理信息系統(tǒng)將分析結(jié)果“視圖化”，結(jié)果表明：

2006-2011年間，中國(guó)30個(gè)省域的資源環(huán)境壓力由西向東、由北到南依次遞減，造成西部、西南部、東北部及中部的資源環(huán)境壓力較大的主要原因是資源利用粗放，技術(shù)水平相對(duì)落后，東部沿海地區(qū)資源環(huán)境壓力較小的主要因素是區(qū)域經(jīng)濟(jì)的發(fā)展水平較高，資源環(huán)境技術(shù)較為先進(jìn)。省域資源環(huán)境承載力由東向西、由南到北依次遞減，制約西部地區(qū)資源環(huán)境承載力的主要因素是區(qū)域經(jīng)濟(jì)發(fā)展水平較低，限制中部資源環(huán)境承載力的主要原因是自然資源的相對(duì)缺乏。省域資源環(huán)境承載率的空間分布亦呈現(xiàn)出比較明顯的由沿海向內(nèi)陸至西部梯狀下降的趨勢(shì)，區(qū)域差異較為明顯。

以Wrl為權(quán)值矩陣的資源環(huán)境壓力、承載力和承載率均存在著較強(qiáng)的全局空間正相關(guān)性，且空間自相關(guān)性較為顯著。中國(guó)省域資源環(huán)境壓力、承載力和承載率空間分布不均衡，存在著較大的地帶性差異，資源環(huán)境壓力、承載力和承載率均顯現(xiàn)高值與高值地區(qū)集聚，低值與低值地區(qū)集聚的特點(diǎn)，即呈現(xiàn)出“中心—外圍”的格局?！?/p>

（作者：江漢大學(xué)商學(xué)院金融與國(guó)際經(jīng)濟(jì)貿(mào)易系，講師，統(tǒng)計(jì)學(xué)博士，主要從事經(jīng)濟(jì)應(yīng)用統(tǒng)計(jì)研究）

注釋：

毛漢英，余丹林，區(qū)域承載力定量研究方法探討，地球科學(xué)進(jìn)展，2001第4期，第549-555頁(yè).

毛漢英，余丹林，環(huán)渤海地區(qū)區(qū)域承載力研究，地理學(xué)報(bào)，2001第3期，第363-371頁(yè).

葉明霞，羅國(guó)云，長(zhǎng)江中上游地區(qū)資源環(huán)境承載力的實(shí)證分析，華東經(jīng)濟(jì)管理，2009年第3期，第1-4頁(yè).

唐凱等，基于集對(duì)分析法的長(zhǎng)株潭城市群資源環(huán)境承載力評(píng)價(jià)，國(guó)土資源科技管理，2012年第1期，第46-53頁(yè).

郭亞軍，綜合評(píng)價(jià)理論、方法及應(yīng)用，科學(xué)出版社，2008年6月，第67-70；107-125頁(yè).

Elhorst．J．P．SpecificationandEstimationofSpatialPanelDataModel，InternationalRegionalScienceReview.2003（26）：2444-2468.