陳紅梅，杭艷紅，楊 林，韓東薇，石 磊，寧 靜

(東北農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，哈爾濱 150030)

0 引 言

水土資源是一個(gè)復(fù)雜的系統(tǒng)，同時(shí)是人類生存發(fā)展的基礎(chǔ)性、戰(zhàn)略性資源[1]。隨著工業(yè)化和城鎮(zhèn)化發(fā)展，其帶來的水土資源污染和短缺問題已成為人類面臨的巨大挑戰(zhàn)[2]。因此，對(duì)水土資源承載力進(jìn)行合理有效地評(píng)價(jià)，了解水土資源的利用結(jié)構(gòu)和分布規(guī)律，深入分析水土資源承載力時(shí)空變化的影響因素，有助于合理調(diào)整水資源供給和土地利用轉(zhuǎn)變，從而促進(jìn)區(qū)域社會(huì)、經(jīng)濟(jì)和環(huán)境協(xié)調(diào)發(fā)展。

國(guó)內(nèi)外學(xué)者對(duì)水土資源承載力進(jìn)行了大量的研究[3-7]，但大多數(shù)的研究都將水資源和土地資源分別作為研究對(duì)象，將水土資源作為整體系統(tǒng)的研究相對(duì)較少[8]。將水土資源作為系統(tǒng)的研究方式主要是通過構(gòu)建水土資源承載力評(píng)價(jià)指標(biāo)體系，采用模糊綜合評(píng)價(jià)法、因子分析、決策分析等方法進(jìn)行評(píng)價(jià)。例如，施開放等結(jié)合可拓學(xué)和熵權(quán)理論，構(gòu)建熵權(quán)可拓決策模型，對(duì)重慶三峽庫(kù)區(qū)水土資源承載力進(jìn)行定量評(píng)價(jià)[9]。周智等將水土資源系統(tǒng)劃分為水資源、土地資源和水土資源耦合3個(gè)子系統(tǒng)，利用因子分析法，對(duì)全國(guó)31個(gè)省水土資源承載力進(jìn)行綜合評(píng)價(jià)[10]。任守德等提出從微觀和宏觀尺度進(jìn)行農(nóng)業(yè)水土資源承載力評(píng)價(jià)的新思路，并對(duì)三江平原農(nóng)墾建三江分局農(nóng)業(yè)水土資源承載力進(jìn)行實(shí)證研究。南彩艷等利用改進(jìn)的SPA水土資源承載力綜合評(píng)價(jià)模型，對(duì)關(guān)中地區(qū)影響水土資源承載力進(jìn)行評(píng)價(jià)和分級(jí)[11]。拓學(xué)森從仿真模擬角度利用系統(tǒng)動(dòng)力學(xué)模型模擬了民勤縣水土資源綜合承載力變化趨勢(shì)[12]。呂添貴等采用信息熵權(quán)法確定耦合評(píng)價(jià)指標(biāo)體系，對(duì)鄱陽(yáng)湖生態(tài)經(jīng)濟(jì)區(qū)水土資源進(jìn)行評(píng)價(jià)[13]；同時(shí)，研究角度也存在較大差異，王薇等從水土資源的空間匹配角度，探討其匹配程度的區(qū)域差異，提出水土資源科持續(xù)利用對(duì)策[14]。張曉青、張衍廣等從預(yù)測(cè)角度，求取了區(qū)域未來水土資源的綜合承載力[15]。呂添貴等將水土資源與經(jīng)濟(jì)發(fā)展進(jìn)行耦合分析[13]。

考慮到水土資源系統(tǒng)具有的動(dòng)態(tài)性和復(fù)雜性，以及子系統(tǒng)間相互作用關(guān)系，同時(shí)綜合黑龍江省的社會(huì)經(jīng)濟(jì)和水土資源狀況，本文采用改進(jìn)的DPSIR水土資源承載力指標(biāo)體系，利用時(shí)序全局因子分析方法，對(duì)黑龍江省水土資源承載力進(jìn)行綜合評(píng)價(jià)，深入分析各個(gè)子系統(tǒng)和各地區(qū)水土資源承載力演化規(guī)律，探討水土資源綜合承載力空間分異特征以及影響因素，摸清黑龍江省的水土資源承載力狀況，從而為提高水土資源利用效率，制定有效水土資源保護(hù)和治理措施提供科學(xué)依據(jù)。

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

1.1 研究區(qū)概況

黑龍江省位于中國(guó)東北部，地理坐標(biāo)121°11′～135°5′E、43°25′～53°33′N，西部屬松嫩平原，東北部為三江平原，北部、東南部為山地，平原海拔50～200 m，年平均氣溫在-5～5 ℃之間。其境內(nèi)有黑龍江、松花江和烏蘇里江三大水系，同時(shí)有興凱湖、鏡泊湖、連環(huán)湖和五大連池4處較大湖泊及星羅棋布的泡沼。全省多年平均地表水資源量為686.0億m3，多年平均地下水資源量為297.44億m3，扣除兩者之間重復(fù)計(jì)算量173.14億m3，全省多年平均水資源量為810.3億m3，人均水量2 160 m3，均低于全國(guó)平均水平[16]。土地總面積為4 730萬(wàn)hm2，占全國(guó)土地總面積的4.9%，其中農(nóng)用地面積3 950.2萬(wàn)hm2，占全省土地總面積的83.5%，是我國(guó)重要的農(nóng)業(yè)大省。2003年實(shí)施東北振興戰(zhàn)略以來，黑龍江省城鎮(zhèn)化步伐不斷加快，水土資源的空間配置也發(fā)生了顯著變化，因此，有必要對(duì)東北振興戰(zhàn)略實(shí)施以來黑龍江省水土資源承載進(jìn)行綜合評(píng)價(jià)，同時(shí)分析其空間分異特征，為保護(hù)和合理有效地利用水土資源提供理論依據(jù)。

1.2 數(shù)據(jù)來源

研究數(shù)據(jù)主要來源于《黑龍江省統(tǒng)計(jì)年鑒》(2005-2014)、《中國(guó)城市統(tǒng)計(jì)年鑒》(2005-2014)、《中國(guó)水利統(tǒng)計(jì)年鑒》(2005-2014)、《黑龍江省水資源公報(bào)》(2005-2014)等統(tǒng)計(jì)資料以及相關(guān)政府報(bào)告資料。

2 研究方法

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

水土資源是一個(gè)復(fù)雜的綜合系統(tǒng)，其結(jié)構(gòu)和組成要素具有復(fù)雜性和動(dòng)態(tài)性，而DPSIR模型是一種在環(huán)境系統(tǒng)中廣泛使用的評(píng)價(jià)指標(biāo)體系概念模型，它將表征一個(gè)自然系統(tǒng)的評(píng)價(jià)指標(biāo)分成驅(qū)動(dòng)力(Driving Force)、壓力( Pressure)、狀態(tài)(State)、影響(Impact)、響應(yīng)(Response)5個(gè)子系統(tǒng)。進(jìn)一步，從經(jīng)濟(jì)社會(huì)發(fā)展需求、環(huán)境質(zhì)量提高、資源合理匹配等發(fā)展問題的角度，選取21項(xiàng)代表性指標(biāo)，構(gòu)建黑龍江省水土資源承載力綜合評(píng)價(jià)指標(biāo)體系(表1)。

表1 水土資源承載力綜合評(píng)價(jià)指標(biāo)體系Tab.1 Indexes system、expression and significance of soil and water resources carrying capacity

續(xù)表1 水土資源承載力綜合評(píng)價(jià)指標(biāo)體系

2.2 時(shí)序全局因子分析模型的建立

水土資源承載力系統(tǒng)具有顯著的動(dòng)態(tài)性，是多維數(shù)據(jù)、多種因素相互影響制約的復(fù)雜系統(tǒng)，同時(shí)，數(shù)據(jù)的積累表征了我們擁有大量的按時(shí)間順序排列的平面數(shù)據(jù)表數(shù)列，構(gòu)成了水土資源綜合承載力的數(shù)據(jù)匣，稱為時(shí)序立體數(shù)據(jù)表[17]。時(shí)序全局因子分析方法將多維的動(dòng)態(tài)時(shí)序立體數(shù)據(jù)表通過全局的主成分公因子變換到相同的主超平面上，使得各年份公因子具有相同的構(gòu)成，根據(jù)不同的時(shí)序變化，數(shù)據(jù)進(jìn)行交換、組合。對(duì)于具有動(dòng)態(tài)性特征的水土資源綜合系統(tǒng)，通過時(shí)序全局因子分析方法，構(gòu)成了整個(gè)系統(tǒng)的統(tǒng)一簡(jiǎn)化子空間，保證了水土資源承載力綜合系統(tǒng)分析的統(tǒng)一性、整體性。

具體步驟如下：

(1)時(shí)序全局立體數(shù)據(jù)表的構(gòu)建。若以n個(gè)城市作為水土資源承載力的樣本點(diǎn)，選取p個(gè)指標(biāo)，以x1，x2，…，xp為變量指標(biāo)，就構(gòu)成了時(shí)序立體數(shù)據(jù)表K，

K={Xt∈Rn×p，t=1,2，…，T} (t為時(shí)間跨度)

在t時(shí)刻數(shù)據(jù)表中的Xt可表示為：

(2)指標(biāo)數(shù)據(jù)的無(wú)量綱標(biāo)準(zhǔn)化處理。水土資源綜合承載力評(píng)價(jià)指標(biāo)數(shù)據(jù)的量綱不同，各個(gè)指標(biāo)數(shù)據(jù)值相差很大，考慮數(shù)據(jù)的特點(diǎn)和分析目的，為了消除指標(biāo)受量綱影響而獲得的不可公度性，對(duì)評(píng)價(jià)指標(biāo)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，采用兩種方法。

對(duì)于指標(biāo)值越大越好的數(shù)據(jù)：

X*ij=[xij-min(xj)]/[max(xj)-min(xj)]

(i=1,2,3,…，n；j=1,2,3，…,p)

對(duì)于指標(biāo)值越大越好的數(shù)據(jù)：

X*ij=[max(xj)-xij]/[max(xj)-min(xj)]

(i=1,2,3,…，n；j=1,2,3，…,p)

其中，max(xj)、min(xj)分別為指標(biāo)的最大值與最小值。由此無(wú)量綱標(biāo)準(zhǔn)化后的決策矩陣記為X*，X*=(x*ij)n×p。

(3)因子模型適用性分析。利用KMO對(duì)評(píng)價(jià)的指標(biāo)變量進(jìn)行分析，確定其是否有較強(qiáng)的相關(guān)性，以判斷其是否適合使用因子分析方法進(jìn)行分析。

(4)計(jì)算X*矩陣的協(xié)方差陣R，求R的m個(gè)特征值λ1，λ2，…，λm(m

(5)利用最大方差法進(jìn)行因子旋轉(zhuǎn)，求得全局主公因子表達(dá)式：

Y1=u11x1+u21x2+…+um1xm

Y2=u12x1+u22x2+…+um2xm

…

Ym=u1mx1+u2mx2+…+ummxm

Y1，Y2，…，Ym表示第 1、2至m個(gè)全局主公因子。

(6)構(gòu)造綜合承載力得分函數(shù)，計(jì)算各子系統(tǒng)承載力分值。根據(jù)主公因子貢獻(xiàn)率，構(gòu)造綜合承載力得分函數(shù)：F=(λ1Y1+λ2Y2+…+λmYm)/(λ1+λ2+…+λm)。

3 結(jié)果分析

3.1 水土資源承載力評(píng)價(jià)結(jié)果分析

以黑龍江省地級(jí)市為基本單元，對(duì)該區(qū)域12個(gè)地區(qū)的數(shù)據(jù)根據(jù)因子分析的原理和步驟，運(yùn)用SPSS軟件進(jìn)行統(tǒng)計(jì)處理，其KMO檢驗(yàn)值為：0.694，表明適宜進(jìn)行因子分析。據(jù)表2和旋轉(zhuǎn)后的因子載荷矩陣，特征值大于1的6個(gè)公因子累積貢獻(xiàn)率為77.158，且其載荷系數(shù)較高，選取為公因子。

表2 時(shí)序全局因子解釋的總方差情況Tab.2 The total variance explained in the case of the timing factor globally

3.2 各子系統(tǒng)變化趨勢(shì)分析

將各地區(qū)子系統(tǒng)值進(jìn)行求和得出黑龍江省各子系統(tǒng)值，進(jìn)一步對(duì)各子系統(tǒng)的變化趨勢(shì)進(jìn)行分析(圖1)，可以看出各子系統(tǒng)在水土資源系統(tǒng)中具有明顯的層次，從高往低依次為影響系統(tǒng)、驅(qū)動(dòng)力系統(tǒng)、壓力系統(tǒng)和響應(yīng)系統(tǒng)，而狀態(tài)系統(tǒng)為負(fù)值，表明由黑龍江省的社會(huì)經(jīng)濟(jì)“驅(qū)動(dòng)力”帶來的發(fā)展，對(duì)水土資源系統(tǒng)帶來了巨大的“影響”，同時(shí)由于水土資源利用方式的調(diào)整或經(jīng)濟(jì)政策的改變，即“響應(yīng)”系統(tǒng)的變化，從而使水土資源承受了相對(duì)較少的“壓力”，因此，其負(fù)向的狀態(tài)值說明水土資源的利用向著更為合理的方向發(fā)展。

圖1 各子系統(tǒng)發(fā)展趨勢(shì)Fig.1 The trends of each subsystem

3.3 各地區(qū)水土資源承載力演化規(guī)律

3.3.1水土資源承載力時(shí)間變化規(guī)律

黑龍江省水土資源承載力綜合指數(shù)如圖2所示，各綜合承載力指數(shù)為各分項(xiàng)承載力指數(shù)加權(quán)后的總和，水土資源承載力綜合指數(shù)最高為1，指數(shù)值越高，說明承載力越大。黑龍江省各地區(qū)水土資源承載力綜合指數(shù)波動(dòng)較大，總體呈先上升后下降再上升的N型曲線特征，其整體呈增長(zhǎng)趨勢(shì)，將其變化階段劃分為：增長(zhǎng)期(2004-2009年)、衰落期(2009-2011年)、反彈期(2011-2013年)。①在增長(zhǎng)期，大部分地區(qū)水土資源承載力綜合指數(shù)呈上升趨勢(shì)，其水土資源承載力在逐步增強(qiáng)，然而，哈爾濱、黑河、佳木斯在2004年承載力綜合指數(shù)較高，之后緩慢發(fā)展，主要由于2004年?yáng)|北西部發(fā)生百年一遇春夏連旱，導(dǎo)致其他地區(qū)水土資源承載力降低。②衰落期間佳木斯、大慶、雙鴨山、雞西、七臺(tái)河、鶴崗和綏化，水土資源綜合承載力指數(shù)整體呈下降趨勢(shì)，說明水土資源匹配狀態(tài)開始惡化，哈爾濱、黑河、伊春發(fā)展較為平穩(wěn)。③在反彈期，大部分地區(qū)水土資源承載力恢復(fù)增長(zhǎng)趨勢(shì)，只有綏化持續(xù)下降，進(jìn)一步發(fā)現(xiàn)綏化水土資源壓力系統(tǒng)的作用遠(yuǎn)大于驅(qū)動(dòng)力系統(tǒng)，其影響系統(tǒng)呈現(xiàn)負(fù)值，可以說明綏化的水土資源利用方式存在嚴(yán)重的問題。

圖2 各地區(qū)水土資源承載力變化態(tài)勢(shì)Fig.2 Carrying capacity of soil and water resources situation changes of Regions

3.3.2水土資源承載力時(shí)間變異特征

圖3 黑龍江省水土資源承載力區(qū)域變異系數(shù)Fig.3 The regional variation coefficient of soil and water resources capacity in Heilongjiang Province

3.4 黑龍江省水土資源承載力空間分異特征分析

為了探討各地區(qū)水土資源的分異特征，采用系統(tǒng)聚類方法對(duì)2004年、2009年、2013年各地區(qū)綜合評(píng)價(jià)結(jié)果進(jìn)行聚類分析，并將12地級(jí)市分為3個(gè)梯度，在空間上進(jìn)行表達(dá)(如圖3)。

2004年各地區(qū)水土資源承載力具有較大差異，4個(gè)農(nóng)業(yè)分區(qū)[18]內(nèi)各有一個(gè)高值聚集點(diǎn)分別是黑河、佳木斯、大慶和牡丹江，中等梯度為齊齊哈爾、伊春和哈爾濱地區(qū)，其他地區(qū)為低等梯度主要分布在三江平原和綏化地區(qū)。2009年大部分地區(qū)處于高等梯度，僅有東部煤電基地的雙鴨山、七臺(tái)河和雞西地區(qū)，以及綏化地區(qū)處于中等梯度。2013年與2009年相比，綏化從中等梯度退到低等梯度，鶴崗從高等梯度退到中等梯度?？傮w而言，2004-2009年，水土資源承載力得到較大提高，高值區(qū)域以綏化地區(qū)為中心，向各個(gè)方向擴(kuò)張，低值區(qū)分布在東部且呈現(xiàn)明顯縮小態(tài)勢(shì)，2009-2013年，部分地區(qū)呈現(xiàn)下降態(tài)勢(shì)，但水土資源承載力整體發(fā)展趨于良性，呈上升的發(fā)展態(tài)勢(shì)。

圖4 黑龍江省各市水土資源綜合承載力 空間分布變化態(tài)勢(shì)(2004、2009、2013年)Fig.4 Comprehensive changes in the spatial distribution of soil and water resources carrying capacity situation of cities in Heilongjiang Province(2004,2009,2013)

4 結(jié) 語(yǔ)

(1)水土資源承載力系統(tǒng)是由各個(gè)子系統(tǒng)共同作用形成的合力系統(tǒng)。黑龍江省的水土資源系統(tǒng)中各子系統(tǒng)的作用形成明顯的層次，表明黑龍江省的社會(huì)經(jīng)濟(jì)帶來的發(fā)展，給水土資源系統(tǒng)帶來影響，同時(shí)由于水土資源利用方式的改變，減少了這種壓力，從而使水土資源利用方式向著更合理的方向發(fā)展。

(2)東北振興戰(zhàn)略實(shí)施以來，黑龍江省各地區(qū)水土資源承載力的變化呈現(xiàn)“N”字型態(tài)勢(shì)，其演化階段可以分為：增長(zhǎng)期、衰落期和反彈期。水土資源狀況在發(fā)展的初期表現(xiàn)出良好的態(tài)勢(shì)，但隨著社會(huì)經(jīng)濟(jì)的發(fā)展對(duì)其造成破壞，然而由于水土資源利用方式和經(jīng)濟(jì)政策調(diào)整水土資源狀態(tài)仍然可以恢復(fù)發(fā)展的活力。同時(shí)，水土資源承載力與各地區(qū)水土資源利用方式的差異性存在較強(qiáng)的關(guān)聯(lián)性，如綏化水土資源承載力持續(xù)下降，其水土資源壓力系統(tǒng)的作用遠(yuǎn)大于驅(qū)動(dòng)力系統(tǒng)，其影響系統(tǒng)呈現(xiàn)負(fù)值，可以說明綏化的水土資源利用方式存在嚴(yán)重的問題。

(3)水土資源承載力具有明顯的空間分異特征，黑龍江省西部、北部和南部水土資源承載力在不斷地提高，而東部煤炭基地建設(shè)，致使水土資源壓力增加，進(jìn)而水土資源承載力相對(duì)較弱，總體而言，近年來黑龍江省水土資源承載力發(fā)展整體趨于良性，呈上升態(tài)勢(shì)。

[1] 姜秋香,付 強(qiáng),王子龍,等. 三江平原水土資源空間匹配格局[J]. 自然資源學(xué)報(bào),2011,26(2):270-277.

[2] 楊 宇,劉 毅,金鳳君,等. 塔里木河流域綠洲城鎮(zhèn)發(fā)展與水土資源效益分析[J]. 地理學(xué)報(bào),2012,67(2):157-168.

[3] Joardar S D. Carrying capacities and standards as bases towards urban Infra structure planning in India: a case of urban water supply and sanitation[J]. Habitat International, 1998,22(3):327-337.

[4] Oh K, Jeong Y, Lee D, et al. Determining development density using the urban carrying capacity assessment system[J]. Landscape and Urban Planning, 2005,73(1):1-15.

[5] 張永勇,夏 軍,王中根. 區(qū)域水資源承載力理論與方法探討[J]. 地理科學(xué)進(jìn)展,2007,26(2):126-132.

[6] 王長(zhǎng)建,張小雷,杜宏茹,等. 開都河-孔雀河流域水資源承載力水平的綜合評(píng)價(jià)與分析[J]. 冰川凍土，2012,34(4):990-998.

[7] 岳曉燕,宋伶英. 土地資源承載力研究方法的回顧與展望[J]. 水土保持研究,2008,15(1):254-257.

[8] 任守德,付 強(qiáng),王 凱. 區(qū)域水土資源承載力研究進(jìn)展[J]. 農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究,2010,26(2):222-226.

[9] 施開放,刁承泰,孫秀鋒. 基于熵權(quán)可拓決策模型的重慶三峽庫(kù)區(qū)水土資源承載力評(píng)價(jià)[J]. 環(huán)境科學(xué)學(xué)報(bào),2013,33(2):609-616.

[10] 周 智,黃 英,黃 娟. 水土資源承載力區(qū)域差異與影響因素[J]. 水土保持通報(bào),2015,35(2):344-349.

[11] 南彩艷,粟曉玲.基于改進(jìn)SPA的關(guān)中地區(qū)水土資源承載力綜合評(píng)價(jià)[J]. 自然資源學(xué)報(bào),2012,27(1):104-114.

[12] 拓學(xué)森,陳興鵬,薛冰.民勤縣水土資源承載力系統(tǒng)動(dòng)力學(xué)仿真模型研究[J].干旱區(qū)資源與環(huán)境,2006,20(6):78-83.

[13] 呂添貴,吳次芳,游和遠(yuǎn). 鄱陽(yáng)湖生態(tài)經(jīng)濟(jì)區(qū)水土資源與經(jīng)濟(jì)發(fā)展耦合分析及優(yōu)化路徑[J]. 中國(guó)土地科學(xué),2013,27(9):3-10.

[14] 王 薇,呂寧江,王 昕,等. 黃河三角洲水土資源空間匹配格局探析[J]. 水資源與水工程學(xué)報(bào),2014,25(2):66-70.

[15] 張曉青,李玉江.山東省水土資源承載力空間結(jié)構(gòu)研究[J].資源科學(xué),2006,28(2):13-21.

[16] 李德一,張樹文.黑龍江省水資源與社會(huì)經(jīng)濟(jì)發(fā)展協(xié)調(diào)度評(píng)價(jià)[J].干旱區(qū)資源與環(huán)境,2010,24(4):8-11.

[17] 任若恩,王惠文.多元統(tǒng)計(jì)數(shù)據(jù)分析----理論、方法、實(shí)例[M].北京:國(guó)防工業(yè)出版社，1997:168-186．

[18] 彭立學(xué),楊 蓓,劉志強(qiáng). 黑龍江省農(nóng)業(yè)資源環(huán)境分區(qū)預(yù)警與農(nóng)業(yè)發(fā)展對(duì)策的探討[J]. 科學(xué)新聞,2003,(9):27-28.