裴巍，付強，劉東，李天霄，成琨

（東北農(nóng)業(yè)大學(xué)水利與建筑學(xué)院，哈爾濱　150030）

基于改進投影尋蹤模型黑龍江省土地資源生態(tài)安全評價

裴巍，付強*，劉東，李天霄，成琨

（東北農(nóng)業(yè)大學(xué)水利與建筑學(xué)院，哈爾濱150030）

針對Friedman-Tukey投影指標中窗口半徑不易確定，分析目前窗口半徑確定方法不足，提出基于聚類思想窗口半徑確定方法，改進投影尋蹤模型，將其應(yīng)用于黑龍江省各地級市土地資源生態(tài)安全評價，實例分析表明，改進模型優(yōu)于原模型，所得投影點分布與正態(tài)分布差異更大。空間差異分析表明，哈爾濱、齊齊哈爾、七臺河、綏化土地資源生態(tài)安全狀況較差，屬于Ⅳ級，雙鴨山、佳木斯、黑河、雞西、鶴崗和牡丹江屬于Ⅲ級或Ⅱ級，生態(tài)安全狀況適中，土地資源生態(tài)安全狀況較好地區(qū)是大慶、伊春和大興安嶺，屬于Ⅰ級。從自然因素上看，大多數(shù)區(qū)域?qū)儆冖窦壔颌蚣?，表明黑龍江總體自然狀況較好；從經(jīng)濟因素上看，經(jīng)濟狀況差，改善土地資源生態(tài)安全能力有限；從社會因素上看，區(qū)域等級分布均衡，表明社會因素對土地壓力適中，有待改善。

投影尋蹤；窗口半徑；聚類；土地資源；生態(tài)安全；黑龍江

網(wǎng)絡(luò)出版時間2016-7-21 14:12:07［URL］http://www.cnki.net/kcms/detail/23.1391.S.20160721.1412.010.html

裴巍,付強,劉東,等.基于改進投影尋蹤模型黑龍江省土地資源生態(tài)安全評價[J].東北農(nóng)業(yè)大學(xué)學(xué)報,2016,47(7):92-100.

Pei Wei,Fu Qiang,Liu Dong,et al.Assessing land resource ecological security in the Heilongjiang Province based on improved projection pursuit model[J].Journal of Northeast Agricultural University,2016,47(7):92-100.(in Chinese with English abstract)

土地資源是人類賴以生存和發(fā)展的基礎(chǔ)性自然資源，其生態(tài)安全問題是土地規(guī)劃、管理和保護的出發(fā)點和落腳點，關(guān)系到糧食、生態(tài)、經(jīng)濟安全等方面［1］。近年來，水土流失嚴重，區(qū)域內(nèi)生態(tài)環(huán)境惡化。土地資源生態(tài)安全評價可評估區(qū)域土地安全狀態(tài)，為土地生態(tài)安全預(yù)警及調(diào)控，土地規(guī)劃、管理及資源優(yōu)化配置提供理論依據(jù)，保護土地資源生態(tài)安全，實現(xiàn)可持續(xù)發(fā)展。

在評價指標體系選取方面，Hua等基于PSR模型從土地生態(tài)壓力、土地生態(tài)狀態(tài)和人文社會響應(yīng)三個方面構(gòu)建評價指標體系［2-5］；劉勇等從自然、經(jīng)濟和社會系統(tǒng)安全三個方面選取典型指標構(gòu)建土地資源生態(tài)安全評價指標體系［6-7］。在評價方法選擇方面，有綜合指數(shù)法［2，8］，層次分析法［7，9-10］，熵權(quán)法［2，11］，物元分析法［12-13］，主成分分析法［3，14］，主成分聚類分析［15］，模糊綜合評價法［16］，灰色綜合評價法［17］，投影尋蹤評價法［18］等。在模型構(gòu)建方面，宋戈等運用生態(tài)足跡模型研究土地資源生態(tài)安全［19-20］，Zhang等開展景觀生態(tài)學(xué)模型方面研究［21-22］；此外，將新方法和模型應(yīng)用于土地資源生態(tài)安全評價之中，如正態(tài)云模型［23］，Double Land Ecological Security（DLES）方法［24］，突變模型［25］等。同時，隨“3S”技術(shù)發(fā)展，基于“3S”技術(shù)生態(tài)安全評價逐步發(fā)展［26-27］。土地資源生態(tài)安全評價已成為土地生態(tài)安全研究重要內(nèi)容。

投影尋蹤是處理高維非線性、非正態(tài)數(shù)據(jù)有效統(tǒng)計方法，通過極大化或極小化反映數(shù)據(jù)聚類程度投影指標，把高維數(shù)據(jù)投影到低維空間，投影后數(shù)據(jù)可以進行分類或預(yù)測，找到高維數(shù)據(jù)結(jié)構(gòu)特征［28］。投影尋蹤方法中，投影指標選擇是核心問題。Friedman和Tukey提出以樣本投影值標準差和局部密度乘積作為投影指標應(yīng)用，基本思想是投影點總體分散，局部聚集［29-31］。Friedman-Tukey投影指標中，局部密度窗口半徑選擇是核心問題。Friedman等采用R=0.1Sz，Sz為投影值標準差［31］；陳曜等規(guī)定R范圍R∈［rmax，2p］，rmax=max｛r（i，j）｝，p為待分類樣本維數(shù)［32］；張欣莉等提出R取常數(shù)方案［33］；樓文高等給出R范圍R∈［rmax/5，rmax/3］［34］。上述方案對于不同數(shù)據(jù)，確定半徑R會過大或過小，造成分類效果較差。

本文基于Friedman-Tukey投影指標，針對投影數(shù)據(jù)，提出基于聚類思想窗口半徑確定方法，改進投影尋蹤模型，將其應(yīng)用于黑龍江省13個地市土地資源生態(tài)安全評價，確定各地市相對土地資源生態(tài)安全等級，以期為區(qū)域土地資源生態(tài)安全預(yù)警、調(diào)控、土地規(guī)劃管理及資源配置提供理論指導(dǎo)。

1　研究區(qū)域與研究方法

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

黑龍江省位于中國東北部，是我國緯度最高省份。北部和東部以黑龍江及烏蘇里江為界，與俄羅斯相望，西部與內(nèi)蒙古自治區(qū)毗鄰，南部與吉林省接壤；南北相距1 120 km，東西寬930 km。黑龍江省地勢大致西北部、北部和東南部較高，東北部、南部較低。氣候?qū)僦袦貛А⒑疁貛Т箨懶约撅L(fēng)氣候。區(qū)域總面積463 730 km2，其中耕地面積158 660 km2，林地面積218 370 km2，濕地面積86 700 km2，其中天然濕地面積55 600 km2，草地面積20 630 km2。區(qū)域內(nèi)有黑龍江、烏蘇里江、松花江、綏芬河四大水系，有興凱湖、鏡泊湖、連環(huán)湖和五大連池四處較大湖泊。區(qū)域概況見圖1。

黑龍江是農(nóng)業(yè)大省，耕地資源豐富，近年來，區(qū)域內(nèi)過度開墾、采伐、放牧，農(nóng)藥化肥過量施用，耕地長期重用輕養(yǎng)，土壤污染問題嚴重，有機質(zhì)含量下降，地力減退；林地、濕地、草原面積逐漸減少，水土涵養(yǎng)能力下降，土地沙化鹽漬化現(xiàn)象嚴重，水土流失面積逐年增加；土地后備資源減少，土地生態(tài)功能減弱。

本文研究對象是黑龍江省13個地市，具體包括：哈爾濱市、齊齊哈爾市、雞西市、鶴崗市、雙鴨山市、大慶市、伊春市、佳木斯市、七臺河市、牡丹江市、黑河市、綏化市和大興安嶺地區(qū)。所涉及數(shù)據(jù)主要來自于黑龍江省統(tǒng)計年鑒（2014）和黑龍江省農(nóng)墾年鑒（2014）；黑龍江省水利建設(shè)資料（2012）；相關(guān)部門官方網(wǎng)站。

圖1　黑龍江省區(qū)域概況Fig.1Regional map of the Heilongjiang Province

1.2Friedman-Tukey投影尋蹤模型及局部密度窗口半徑確定

1.2.1Friedman-Tukey投影尋蹤模型簡介

Friedman與Tukey提出Friedman-Tukey投影指標［31］，投影指標函數(shù)為：

1.2.2已有局部密度窗口半徑確定方法不足

投影指標中標準差Sz反映投影點數(shù)據(jù)分散程度，Sz越大則數(shù)據(jù)整體分散程度越大。局部窗口密度Dz反映投影點聚集程度，Dz中u（R-r（i，j））主要作用是篩選數(shù)據(jù)，投影尋蹤模型對窗口半徑要求：對于點團內(nèi)投影點，應(yīng)R≥r（i，j），此時u（R-r（i，j））＝1，保留R-（ri，j），即保留點團內(nèi)投影點信息；對于點團與點團之間，應(yīng)有R＜r（i，j），此時u（R-r（i，j））=0，即舍棄R-r（i，j），舍棄點團間投影點信息?？梢奟過大，過多投影點在一個點團內(nèi)，分類效果不明顯；R過小，點團內(nèi)投影點過少，點團過多。

假定在某一方向下投影數(shù)據(jù)為z=｛0.8、0.9、1.0、1.1、1.2；1.8、1.9、2.0、2.1、2.2；2.8、2.9、3.0、3.1、3.2｝，顯然這組數(shù)據(jù)被分成3個點團，z1=｛0.8、0.9、1.0、1.1、1.2｝，z2=｛1.8、1.9、2.0、2.1、2.2｝，z3=｛2.8、2.9、3.0、3.1、3.2｝，滿足點團內(nèi)投影點距離小于點團之間投影點距離，即點團內(nèi)差異均小于點團間差異。把上述數(shù)據(jù)作為測試數(shù)據(jù)。按照投影尋蹤模型對窗口半徑要求，R取值范圍R∈［0.4，0.6），此時對于點團內(nèi)投影點滿足R≥r（i，j），對于點團之間投影點滿足R＜r（i，j）。

①R=0.1Sz，Sz為投影值標準差。對于測試數(shù)據(jù)，樣本間最小距離rmin（i，j）=0.1。計算得到樣本方差Sz=0.8577，R=0.1Sz=0.8577?［0.4，0.6）且R＜rmin（i，j），即對于任意r（i，j），有u（R-r（i，j））=0，此時，每個點團內(nèi)僅一個樣本點，數(shù)據(jù)被分成15個點團，并不合適。以上分析表明，對于某些特定數(shù)據(jù)，R=0.1Sz不能正確分類。

②R∈［rmax，2p］，rmax=max｛r（i，j）｝，p為待分類樣本維數(shù)。對于測試數(shù)據(jù)，rmax=2.4，R＞rmax=2.4，則R?［0.4，0.6），可見窗口半徑過大，導(dǎo)致窗口過寬，所有樣本點集中在一個點團內(nèi)。對于特定數(shù)據(jù)，R∈［rmax，2p］分類效果較差。

③R取常數(shù)方案不合適，由以上分析，可知R作用即篩選數(shù)據(jù)，對于不同數(shù)據(jù)R取值不同。對于一組特定數(shù)據(jù)，理論上應(yīng)存在最優(yōu)值，但該值無法預(yù)知。

④R∈［rmax/5，rmax/3］，R∈［rmax/5，rmax/3］=［0.48，0.8］，與最優(yōu)區(qū)間有交叉部分，未完全包含在最優(yōu)區(qū)間［0.4，0.6）。對于特定數(shù)據(jù)而言，區(qū)間過寬，分類結(jié)果未達到最優(yōu)。這種確定R方法不穩(wěn)定，易受個別極端值影響，當個別極端值出現(xiàn)時，rmax會被無限擴大，窗口半徑會無限變大，窗口過寬，分類效果不明顯。

1.2.3局部密度窗口半徑確定方法

分析發(fā)現(xiàn)，R應(yīng)大于或等于點團內(nèi)部所有點與點間距離，另一方面R應(yīng)小于點團與點團之間所有點與點間距離，實現(xiàn)數(shù)據(jù)篩選。在數(shù)據(jù)整體分散而局部聚集下，Sz和Dz同時達到極大。整體分散可以使Sz達到極大；整體分散、局部聚集可保證每一個點團內(nèi)點盡可能聚集，點團與點團間盡量分散，因為點團與點團間點距離被舍棄，不影響Dz，Dz達到極大。為此，本文提出一種基于聚類思想局部密度窗口半徑確定方法，改進投影尋蹤評價模型。

投影尋蹤模型把p維數(shù)據(jù)｛xij|j=1，2，…，p｝綜合成以a=｛a1，a2，a3，…，ap｝為投影方向一維投影值zi，即：

式中a-單位長度向量。

在綜合投影指標值時，要求投影值zi散布特征應(yīng)為：局部投影點盡可能密集，最好凝聚成若干個點團，整體上點團與點團間盡可能散開。定義投影指標：

Q（a）=SzDz

各項符號意義同前。要求指標值Q（a）越大越好，這是一個以a（jj=1，2，…，p）為優(yōu)化變量復(fù)雜非線性優(yōu)化問題。可以采用基于實數(shù)編碼遺傳算法RGA求解。

將n個樣本做K均值聚類，假定將n分成k類（點團），每個點團中含有樣本數(shù)為x1，x2，…，xk，且x1+x2+…+xk=n，r（i，j）表示投影點之間距離，即r（i，j）=|z（i）-z（j）|。降序排序，記為r（i，j）（k），k=1，2，…，n2，r（i，j）（k）為排序后序號為k個（ri，j）值，則

R=r（i，j）（p）

其中，p=∑xixj（1≤i，j≤k）。投影尋蹤模型對投影值局部密度要求是保留點團內(nèi)距離r（i，j），去除點團之間距離r（i，j）。即對于窗口半徑而言，要求點團與點團之間點保證R-r（i，j）＜0，即這些R-r（i，j）值不計入Dz中，具體個數(shù)為p=∑xixj（1≤i，j≤k），而對于點團內(nèi)部點要保證R-r（i，j）≥0。點團之間距離r（i，j）必定大于點團內(nèi)距離r（i，j），將r（i，j）降序排序后，第p個值r（i，j）（p）即滿足要求。

1.3評價標準及指標體系構(gòu)建

土地資源生態(tài)安全評價指標選擇是區(qū)域土地資源生態(tài)安全重要內(nèi)容之一，評價指標選擇恰當與否關(guān)系評價結(jié)果優(yōu)劣。土地資源生態(tài)安全受到自然、經(jīng)濟、社會等相關(guān)因素影響，在遵循可操作性、系統(tǒng)性、區(qū)域性等原則基礎(chǔ)上，結(jié)合文獻［1-2，4-5，20］，重點結(jié)合黑龍江省區(qū)域特色，從自然子系統(tǒng)、經(jīng)濟子系統(tǒng)、社會子系統(tǒng)三個方面選取16個指標，見表1。結(jié)合國際、國內(nèi)相關(guān)標準和黑龍江省區(qū)域特性，將所有指標分成4個相對等級，等級越高，土地資源生態(tài)風(fēng)險越高，安全性越低。

本文將土地資源生態(tài)安全評價指標分為正向指標和負向指標，正向指標包括：森林覆蓋率（X1）、草地面積比例（X2）、濕地面積比例（X3）、水土協(xié)調(diào)度（X5），這些指標越大表明區(qū)域涵養(yǎng)水源、保持水土能力越強，土地資源生態(tài)安全狀況越好；人均GDP（X6）、農(nóng)民人均純收入（X7）、環(huán)境投資比例（X10）、農(nóng)電集約度（X16），這些指標值越大，表明區(qū)域抵御外來沖擊和壓力能力越強。負向指標包括：水土流失面積比例（X4）、農(nóng)業(yè)經(jīng)濟比重（X8）、單位耕地化肥負荷（X9）、人口密度（X11）、人口自然增長率（X12）、農(nóng)業(yè)人口比重（X13）、城鎮(zhèn)化水平（X14）、單位土地廢水負荷（X15），這些值越大，表明區(qū)域土地受損情況越嚴重，承載壓力或沖擊力能力越差，生態(tài)安全狀況越差。

為解決量綱及數(shù)量級不同，將上述正向、逆向指標作極差標準化［32-33］，標準化后從單個指標上看，指標越大，安全性越高。

表1　黑龍江省土地資源生態(tài)安全評價指標體系Table 1Evaluation indices and relative grades of land resource ecological security in the Heilongjiang Province

2　評價結(jié)果與分析

2.1最佳投影方向及土地資源生態(tài)安全等級確定

2.2評價結(jié)果對比分析

本文在Friedman-Tukey模型基礎(chǔ)上改進窗口半徑確定方法，為從實例上說明本方法有效性，將其與最常用窗口半徑R=0.1Sz對比，R=0.1Sz時Friedman-Tukey模型評價結(jié)果見表3。

由評價結(jié)果可知，兩種窗口半徑確定方式得到最終評價結(jié)果大致相同，說明本文窗口半徑確定方法有效性。為定量分析兩組數(shù)據(jù)差異，采用W檢驗分析兩組數(shù)據(jù)。

W檢驗由夏皮羅和威爾克1965年提出［35］，主要用于小樣本（n≤50）時檢驗數(shù)據(jù)分布與正態(tài)分布差異，W值越小，與正態(tài)分布差異越大。投影尋蹤模型中較為成熟方法是利用數(shù)據(jù)分布與高斯分布差異定義投影指標［31，36］，可用W值判斷最終結(jié)果。

W值具體計算過程見文獻［35］，得到R=0.1Sz時W值為0.9419，采用本文窗口半徑確定方法W值為0.9317，可見本文得到結(jié)果與正態(tài)分布差異較大，結(jié)果相對優(yōu)良。

表2　黑龍江省土地資源生態(tài)安全性評價結(jié)果Table 2Assessment results of land resource ecological security in the Heilongjiang Province

表3　 Friedman-Tukey 模型評價結(jié)果Table 3 Assessment results of Friedman-Tukey model

2.3黑龍江省土地資源生態(tài)安全空間差異分析

為了解不同指標系統(tǒng)對土地資源生態(tài)安全影響，結(jié)合最佳投影方向，從自然子系統(tǒng)、經(jīng)濟子系統(tǒng)、社會子系統(tǒng)三個方面分別計算投影值。以自然子系統(tǒng)為例，將指標X1，X2，X3，X4，X5歸一化后結(jié)合投影方向前五個分量（0.4847，0.2869，0.5807，0.1797，0.1681），計算自然子系統(tǒng)角度投影值，再根據(jù)相應(yīng)分級標準，得到分界界點0.4012，0.6503，0.9042，從自然子系統(tǒng)角度評價結(jié)果，經(jīng)濟子系統(tǒng)和社會子系統(tǒng)計算過程類似，計算結(jié)果如表4所示。

為分析黑龍江省土地資源生態(tài)安全區(qū)域差異，利用ArcGis可視化功能，繪制黑龍江省土地資源生態(tài)安全空間差異分布圖，繪制自然、經(jīng)濟、社會子系統(tǒng)下黑龍江省土地資源生態(tài)安全空間差異分布圖，見圖2?？傮w來看，黑龍江省土地資源生態(tài)安全等級為Ⅲ級。從各個子系統(tǒng)來看，自然子系統(tǒng)和社會子系統(tǒng)等級為Ⅱ級，經(jīng)濟子系統(tǒng)等級為Ⅲ級。從自然因素來看，黑龍江省森林覆蓋率為45.70%，濕地面積比例為12.29%，高于全國平均水平24.83%，5.58%，表明黑龍江省自然狀況相對良好，土地資源生態(tài)安全基礎(chǔ)較好。從社會因素上看，黑龍江省人口密度、城鎮(zhèn)化水平較低，分別為81.08人·km-2，48.39%，低于全國平均水平141.74人·km-2，53.73%，對土地壓力較小。從經(jīng)濟因素上看，黑龍江省總體經(jīng)濟水平較低，人均GDP為38 058元，低于全國均值41 980元，經(jīng)濟對農(nóng)業(yè)依賴較大，農(nóng)業(yè)經(jīng)濟比重為17.50%，高于全國均值10%，對土地壓力較大。

從自然子系統(tǒng)角度來看，Ⅰ級地區(qū)包括大慶、大興安嶺、雙鴨山和伊春4個地區(qū)，Ⅳ級地區(qū)有哈爾濱、齊齊哈爾和七臺河，其余地區(qū)分別為Ⅱ級和Ⅲ級。大興安嶺、伊春、雙鴨山森林覆蓋率均超過60%，伊春和大興安嶺更高達88.34%和81.23%，森林資源豐富，涵養(yǎng)水源和水土保持能較強。大興安嶺、大慶和雙鴨山濕地面積比例均超過20%，濕地具有抵御洪水、蓄洪防旱、控制污染、抵御土壤侵蝕等作用，對土地資源生態(tài)安全起到積極作用。但哈爾濱、齊齊哈爾、七臺河等地土地資源生態(tài)安全風(fēng)險較高，這些地區(qū)在發(fā)展經(jīng)濟同時，開發(fā)耕地資源，對森林、草地等資源有不同程度破壞，影響區(qū)域土地資源生態(tài)安全。

由經(jīng)濟子系統(tǒng)可知，黑龍江經(jīng)濟總體欠發(fā)達，經(jīng)濟對農(nóng)業(yè)依賴程度較強，其中綏化、黑河、大興安嶺農(nóng)業(yè)經(jīng)濟比重超過40%，土地過度開墾和不合理利用導(dǎo)致水土流失問題嚴重，土壤荒漠化現(xiàn)象日愈嚴重，威脅區(qū)域土地資源生態(tài)安全。此外，在追求經(jīng)濟發(fā)展同時，忽略區(qū)域生態(tài)建設(shè)，資源約束趨緊，土地資源生態(tài)系統(tǒng)遭到不同程度破壞。

由社會子系統(tǒng)可知，分布比較均衡，其中Ⅰ級城市有3個，Ⅱ級城市有4個，Ⅲ級城市有2個，Ⅳ城市有5個。除哈爾濱、齊齊哈爾、大慶和綏化外，其余地區(qū)人口密度均低于80人·km-2，對土地壓力較小，一定程度上起到保護耕地、森林等作用。在污染程度上，哈爾濱、大慶、齊齊哈爾單位土地廢水負荷較高，這些城市經(jīng)濟發(fā)展相對較好，工業(yè)“三廢”排放量較大，對土地安全造成壓力。

表4　自然、經(jīng)濟、社會子系統(tǒng)下評價結(jié)果Table 4Assessment results of Natural，Economic，Social subsystem

圖2　黑龍江省土地資源生態(tài)安全等級Fig.2Grades of land resources ecological security Heilongjiang Province

綜合所有因素，哈爾濱、齊齊哈爾、七臺河、綏化屬于Ⅳ級，土地資源安全狀況較差，主要原因是這些地區(qū)經(jīng)濟相對發(fā)達，人口稠密，森林、草地、濕地面積比例較低、城鎮(zhèn)化水平較高。雞西、鶴崗、牡丹江屬于Ⅲ級，雙鴨山、黑河、佳木斯屬于Ⅱ級，這些地區(qū)總體上土地資源生態(tài)安全狀況適中，雖然一些指標不佳，但是其他指標表現(xiàn)較好恢復(fù)性。土地資源生態(tài)安全狀況良好的是大興安嶺、大慶和伊春，屬于第Ⅰ級。大慶草地面積、濕地面積比例明顯高于其他地區(qū)。大興安嶺和伊春森林覆蓋率較高，自然狀況較好，涵養(yǎng)水源、保持水土能力較高，且區(qū)域人口密度較低，地廣人稀，對土地壓力較小。黑龍江作為國家糧食基地，耕地面積和糧食產(chǎn)量均位于全國首位。然而，過度開墾土地嚴重破壞區(qū)域土地生態(tài)，農(nóng)業(yè)經(jīng)濟附加值不高，在經(jīng)濟結(jié)構(gòu)中處于弱勢地位，但作用不可或缺，所以在國家層面上應(yīng)考慮如何平衡和補償農(nóng)業(yè)，降低土地壓力，保護地方土地資源生態(tài)安全。從地方角度看，應(yīng)關(guān)注并解決環(huán)境投資較低、重耕輕養(yǎng)、土地質(zhì)量下降，破壞區(qū)域土地資源生態(tài)安全問題。

3　結(jié)論

a.改進模型中基于聚類思想窗口半徑確定方法，滿足投影尋蹤模型對窗口半徑要求。

b.本文窗口半徑確定方法和R=0.1Sz比較，W檢驗結(jié)果表明改進模型優(yōu)于原模型。

c.從自然、經(jīng)濟、社會子系統(tǒng)角度分別分析黑龍江省各地市土地資源生態(tài)安全狀況，結(jié)果表明黑龍江總體自然狀況良好，自然基礎(chǔ)較好，但黑龍江地區(qū)經(jīng)濟水平較低，經(jīng)濟對農(nóng)業(yè)依賴程度較高，土地壓力較大。個別地區(qū)如哈爾濱、齊齊哈爾、綏化等地土地資源生態(tài)安全狀況較差。

［1］徐美，朱翔，劉春臘.基于RBF湖南省土地生態(tài)安全動態(tài)預(yù)警［J］.地理學(xué)報，2012，67（10）:1411-1422.

［2］Hua Y E，Yan M A，Dong L.Land ecological security assessment for bai autonomous prefecture of dali based using PSR modelwith data in 2009 as case［J］.Energy Procedia，2011，5（22）:2172-2177.

［3］Ying J，Hai-Shan C，Xue-Fei S.Evaluation of land ecological securityinHubeiProvince［C］//ManagementScienceand Engineering（ICMSE），2013 International Conference on.IEEE，2013:2236-2242.

［4］張虹波，劉黎明，張軍連，等.黃土丘陵區(qū)土地資源生態(tài)安全及其動態(tài)評價［J］.資源科學(xué)，2007，29（4）:193-200.

［5］Pei L，Du L，Yue G.Ecological security assessment of Beijing based on PSR model［J］.Procedia Environmental Sciences，2010，2（1）:832-841.

［6］劉勇，劉友兆，徐萍.區(qū)域土地資源生態(tài)安全評價--以浙江嘉興市為例［J］.資源科學(xué)，2004，26（3）:69-75.

［7］Yang S M.Evaluation on ecological security of land resources［J］. Advanced Materials Research，2013，664:129-132.

［8］Wang J G，Zhao Z Z，Zhao G R.Evaluation on land ecological security in Hainan Island，China［J］.Asian Agricultural Research，2009:69-73.

［9］Zhao Y Z，Zou X Y，Cheng H，et al.Assessing the ecological security of Tibetan plateau:Methodology and a case study for Lhaze County［J］.Journal of Environmental Management，2006，80（2）:120-131.

［10］李玉平，蔡運龍.河北省土地生態(tài)安全評價［J］.北京大學(xué)學(xué)報:自然科學(xué)版，2007，43（6）:784-789.

［11］鄭敬剛，許淑娜，喻錚錚，等.中原經(jīng)濟區(qū)核心區(qū)土地資源生態(tài)安全評價［J］.水土保持研究，2013，20（3）:273-277.

［12］黃輝玲，羅文斌，吳次芳，等.基于物元分析土地生態(tài)安全評價［J］.農(nóng)業(yè)工程學(xué)報，2010，26（3）:316-322.

［13］余建，房莉，倉定幫，等.熵權(quán)模糊物元模型在土地生態(tài)安全評價中應(yīng)用［J］.農(nóng)業(yè)工程學(xué)報，2012，28（5）:260-266.

［14］杜忠潮，韓申山.基于主成分分析土地生態(tài)安全評價實證研究——以陜西省10個省轄市為例［J］.水土保持通報，2009，29（6）:198-202.

［15］鮑艷，胡振琪，柏玉，等.主成分聚類分析在土地利用生態(tài)安全評價中應(yīng)用［J］.農(nóng)業(yè)工程學(xué)報，2006，22（8）:87-90.

［16］Su S，Chen X，DeGloria S D，et al.Integrative fuzzy set pair model for land ecological security assessment:a case study of Xiaolangdi Reservoir Region，China［J］.Stochastic Environmental Research and Risk Assessment，2010，24（5）:639-647.

［17］Xi C L，Jing W，Lan B，et al.The pattern of multiple objects and multiple stages grey relationship for land ecological security assessment［J］.Journal of Grey System，2008，20（4）:351-358.

［18］Gao Y，Wu Z，Lou Q，et al.Landscape ecological security assessment based on projection pursuit in Pearl River Delta［J］.Environmental Monitoring and Assessment，2012，184（4）:2307-2319.

［19］宋戈，韓天宇，王越.基于生態(tài)足跡齊齊哈爾市土地承載力研究［J］.東北農(nóng)業(yè)大學(xué)學(xué)報，2014，45（8）:34-40.

［20］Wang H，Sun H，Wang X，et al.Assessment of the ecological security in the three gorges reservoir area by using the ecological footprint method［J］.Journal of Mountain Science，2012，9（6）:891-900.

［21］Zhang Z，Liu S，Dong S.Ecological security assessment of yuan river watershed based on landscape pattern and soil erosion［J］. Procedia Environmental Sciences，2010，2（6）:613-618.

［22］萬利，陳佑啟，譚靖，等.北京郊區(qū)生態(tài)安全動態(tài)評價與分析［J］.地理科學(xué)進展，2009，28（2）:238-244.

［23］張楊，嚴金明，江平，等.基于正態(tài)云模型湖北省土地資源生態(tài)安全評價［J］.農(nóng)業(yè)工程學(xué)報，2013，29（22）:252-258.

［24］Xu L，Yin H，Li Z，et al.Land ecological security evaluation of Guangzhou，China［J］.International Journal of Environmental Research&Public Health，2014，11（10）:10537-10558.

［25］Su S，Li D，Yu X，et al.Assessing land ecological security in Shanghai（China）based on catastrophe theory［J］.Stochastic Environmental Research&Risk Assessment，2011，25（6）:737-746.

［26］喻鋒，李曉兵，王宏，等.皇甫川流域土地利用變化與生態(tài)安全評價［J］.地理學(xué)報，2006，61（6）:645-653.

［27］呂建樹，吳泉源，張祖陸，等.基于RS和GIS濟寧市土地利用變化及生態(tài)安全研究［J］.地理科學(xué)，2012，32（8）:928-935.

［28］Kruscal J B.Toward a practical method which helps uncover the structure of a set of multivariate observations by finding the linear transformation which optimizes a new index of condensation［J］. Statistical Computer.Academic，1969:427-440.

［29］李芳花，王柏，孫艷玲，等.基于微粒群算法的投影尋蹤模型對調(diào)虧灌溉模式評價［J］.東北農(nóng)業(yè)大學(xué)學(xué)報，2013，44（2）:77-81.

［30］李陶，付強，丁紅.基于投影尋蹤的集對分析模型及其在地下水質(zhì)評價中應(yīng)用［J］.東北農(nóng)業(yè)大學(xué)學(xué)報，2010，41（4）:48-52.

［31］Friedman J H，Turkey J W.A projection pursuit algorithm for exploratory data analysis［J］.IEEE Trans On Computer，1974，23（9）:881-890.

［32］陳曜，丁晶，趙永紅.基于投影尋蹤原理的四川省洪災(zāi)評估［J］.水利學(xué)報，2010，41（2）:220-225.

［33］張欣莉，丁晶，李祚泳，等.投影尋蹤新算法在水質(zhì)評價模型中應(yīng)用［J］.中國環(huán)境科學(xué)，2000，20（2）:187-189.

［34］樓文高，喬龍.投影尋蹤分類建模理論新探索與實證研究［J］.數(shù)理統(tǒng)計與管理，2015，34（1）:47-58.

［35］Shapiro S S，Wilk M B.An analysis of variance test for normality（complete samples）［J］.Biometrika，1965，52（4）:591-611.

［36］Hall P.On polynomial-based projection indices for exploratory projection pursuit［J］.The Annals of Statistics，1989，17（2）:589-605.

AssessinglandresourceecologicalsecurityintheHeilongjiang Province based on improved projection pursuit model

PEI Wei,FU Qiang,LIU Dong,LI Tianxiao,CHENG Kun
(School of Water Conservancy and Architecture,Northeast Agricultural University,Harbin150030,China)

For the cutoff radius of Friedman-Tukey projection index was not easy determined,the shortcomings of the existing method for determining the cutoff radius were analyzed.For these shortcomings,the method of determining the cutoff radius based on clustering was proposed.This improved model was applied to assess land resource ecological security of 13 cities located in the Heilongjiang Province.Case analysis showed that the results of the improved model was better than the original model,the difference between the projected point distribution and normal distribution was even greater.Spatial difference analysis reveal that the ecological security of land resources in Harbin, Qigihar,Qitaihe,Suihua were poorer,belong to classⅣ.Jixi,Hegang,Mudanjiang,Heihe,Shuangyashan,Jiamusi belong to class III or II,the ecological security was moderate.The better regions of ecological security of land resources were Yichun,Daqing and Da Hinggan Mountains and belong to classⅠ.In addition,from the perspective of natural factors,most regions belonged to classⅠorⅡ,itshowed that the natural condition of Heilongjiang Province was better;From the perspective of economic factors,the economic situation of Heilongjiang Province was poor,and the ability to improve the ecological security of land resources was limited;From the perspective of social factors,the regional level distribution was more balanced,which showed that social factors were moderate and there was room for improvement.

projection pursuit;cutoff radius;clustering;land resources;ecological security; Heilongjiang

X826

A

1005-9369（2016）07-0092-09

2016-03-19

國家自然科學(xué)基金（51479032，51279031，51579044）；黑龍江省自然科學(xué)基金（E201241）；黑龍江省高校長江學(xué)者后備支持計劃項目；黑龍江省水利廳科技項目（201318，201503）；黑龍江省杰出青年基金（JC201402）

裴?。?985-），男，講師，博士研究生，研究方向為農(nóng)業(yè)水土資源風(fēng)險分析。E-mail:nothing_pw@163.com

付強，教授，博士生導(dǎo)師，研究方向為農(nóng)業(yè)水土資源系統(tǒng)分析、凍融土壤水熱作用機理。E-mail:fuqiang0629@126.com