基于GIS的三江平原耕地利用效益研究

張文琦+姜博+羅沖+初楠臣+李曉慶+王玉迪

摘要：基于經(jīng)濟(jì)、社會(huì)與生態(tài)視角利用多因素綜合評(píng)價(jià)模型并結(jié)合ArcGIS空間分析手段探究2005—2014年三江平原耕地利用效益時(shí)間動(dòng)態(tài)變化與各縣（市、區(qū)）耕地利用效益空間變異的特征。結(jié)果表明，2005—2014年三江平原耕地利用效益整體呈上升趨勢(shì)，生態(tài)效益在2005—2009年呈下降趨勢(shì)，隨著農(nóng)業(yè)資源利用效率的提高，2009年以后耕地利用生態(tài)效益基本穩(wěn)定；現(xiàn)代農(nóng)業(yè)綜合改革試驗(yàn)區(qū)的建立，帶來三江平原耕地利用效益空間格局的重塑，東北部低平原區(qū)耕地利用綜合效益高于西南部丘陵地區(qū)；未來三江平原耕地利用效益還有較大的提升空間，實(shí)現(xiàn)農(nóng)業(yè)現(xiàn)代化的同時(shí)應(yīng)更加注重耕地利用生態(tài)效益與經(jīng)濟(jì)、社會(huì)效益的協(xié)調(diào)發(fā)展。

關(guān)鍵詞：耕地利用效益；時(shí)間變化；空間分異；GIS技術(shù)；三江平原；經(jīng)濟(jì)效益；社會(huì)效益；生態(tài)效益；綜合效益

中圖分類號(hào)：F301.21 文獻(xiàn)標(biāo)志碼： A文章編號(hào)：1002-1302（2017）23-0243-05

國(guó)外對(duì)耕地利用效益研究開始的比較早，波蘭農(nóng)業(yè)地理學(xué)家Kostowicki根據(jù)農(nóng)業(yè)類型學(xué)委員會(huì)判定的原則和方法把世界農(nóng)業(yè)劃分為6個(gè)一級(jí)類型、25個(gè)二級(jí)類型和93個(gè)三級(jí)類型[7]；Herdt等提出從生態(tài)效益、社會(huì)效益和經(jīng)濟(jì)效益3個(gè)方面來衡量農(nóng)業(yè)系統(tǒng)的可持續(xù)發(fā)展能力[8-9]。國(guó)內(nèi)學(xué)者傅伯杰等認(rèn)為應(yīng)依據(jù)土地可持續(xù)利用的特點(diǎn)，從生態(tài)、社會(huì)、經(jīng)濟(jì)角度建立土地可持續(xù)利用的評(píng)價(jià)指標(biāo)體系[10]；黃秋昊等指出土地可持續(xù)利用的評(píng)價(jià)指標(biāo)應(yīng)該具有時(shí)間序列性與空間敏感性[11]；林庶民等分別采用主成分分析、聚類分析和相關(guān)性分析等方法進(jìn)行實(shí)證研究，揭示了區(qū)域耕地利用效益的時(shí)空差異和特征[12-14]。

但目前我國(guó)耕地利用效益的研究主要針對(duì)耕地質(zhì)量下降明顯的一些人口密集地區(qū)，對(duì)縣域耕地利用效益時(shí)空分異的研究比較少，針對(duì)三江平原耕地利用效益的時(shí)間動(dòng)態(tài)變化規(guī)律和空間分異的研究目前還鮮有報(bào)道。本研究在參考已有研究成果的基礎(chǔ)上，建立三江平原耕地利用效益評(píng)價(jià)指標(biāo)體系，應(yīng)用地理信息系統(tǒng)（geographic information system，GIS）技術(shù)對(duì)三江平原耕地的經(jīng)濟(jì)效益、社會(huì)效益、生態(tài)效益與綜合效益的時(shí)空變化特征進(jìn)行分析，以期為三江平原耕地資源合理利用及綜合效益的提高提供參考。

1研究區(qū)域概況

三江平原位于黑龍江省東部，小興安嶺以東，長(zhǎng)白山以北，介于43°49′55″N～48°27′20″N、129°11′20″E～135°5′10″E之間。主要包括黑龍江、松花江、烏蘇里江沖擊形成的低平地，以及完達(dá)山和南烏蘇里江支流域興凱湖交匯形成的沖積平原。包含佳木斯市、雙鴨山市、鶴崗市、雞西市、七臺(tái)河市所屬的21個(gè)縣（市、區(qū)）和牡丹江所屬的穆棱市、哈爾濱所屬的依蘭縣，總面積108 900 km2，占黑龍江總面積的23.04%。該區(qū)域是溫帶濕潤(rùn)、半濕潤(rùn)的大陸性季風(fēng)氣候，夏季炎熱濕潤(rùn)，冬季寒冷干燥，四季鮮明，1月平均氣溫低于-18 ℃，7月平均氣溫21～22 ℃，年降水量500～650 mm。區(qū)域內(nèi)土壤土層相對(duì)較厚，質(zhì)地主要是以壤土和壤質(zhì)粘土為主，土壤有機(jī)質(zhì)含量相對(duì)較高，已墾耕地土壤以草甸土、白疆土、黑土和暗棕壤為主。三江平原具有典型的河谷平原特點(diǎn)，平均海拔為50～60 m，平原面積大約占總面積的61.2%，是我國(guó)重要的農(nóng)業(yè)區(qū)和商品糧生產(chǎn)基地。

2研究方法

2.1耕地利用效益評(píng)價(jià)指標(biāo)體系的構(gòu)建

在科學(xué)性、客觀性、系統(tǒng)性和可操作性等原則的指導(dǎo)下從經(jīng)濟(jì)、社會(huì)與生態(tài)視角綜合測(cè)度三江平原耕地利用效益。經(jīng)濟(jì)效益從耕地資金、技術(shù)等投入與生產(chǎn)成本降低、糧食產(chǎn)量提高角度判斷耕地的利用效益；社會(huì)效益從宏觀社會(huì)效應(yīng)角度分析耕地利用對(duì)社會(huì)環(huán)境系統(tǒng)作用的效益；生態(tài)效益從生態(tài)環(huán)境角度闡釋在耕地利用過程中對(duì)生物、土壤、水資源等生態(tài)環(huán)境造成的影響[15]。結(jié)合三江平原的實(shí)際情況，經(jīng)因子分析篩選出13個(gè)指標(biāo)，構(gòu)建三江平原耕地利用效益評(píng)價(jià)指標(biāo)體系（表1），該體系由目標(biāo)層A、準(zhǔn)則層B和指標(biāo)層C構(gòu)成，將耕地利用綜合效益定為目標(biāo)層，經(jīng)濟(jì)效益、社會(huì)效益和生態(tài)效益設(shè)為準(zhǔn)則層，各單項(xiàng)指標(biāo)設(shè)為指標(biāo)層。

2.2耕地利用效益評(píng)價(jià)指標(biāo)權(quán)重的確定

運(yùn)用層次分析法（AHP）從主觀視角分析各項(xiàng)指標(biāo)對(duì)耕地利用綜合效益的作用強(qiáng)度[16]。咨詢有代表性的學(xué)者專家，比較所列指標(biāo)的重要程度并逐層進(jìn)行判斷打分，進(jìn)而構(gòu)建判斷矩陣，設(shè)為A=（aij）n×n，其中aij=1/aji（i≠j，i、j=1，2，…，n）；aij=1（i=j）。判斷矩陣中標(biāo)度值通常依據(jù)比較標(biāo)度法得到，2個(gè)指標(biāo)i、j同等重要，標(biāo)度值為1，微弱重要為3，明顯重要為5，突出重要為7，絕對(duì)重要為9，相鄰判斷的中間值為2、4、6、8。依據(jù)方根法計(jì)算最大特征根所對(duì)應(yīng)的特征向量，進(jìn)而得到各項(xiàng)指標(biāo)的重要性權(quán)值，經(jīng)過層次排序和一致性檢驗(yàn)，得出三江平原耕地利用效益各評(píng)價(jià)指標(biāo)的權(quán)重Wi，如表2所示。

2.3數(shù)據(jù)處理

以各項(xiàng)指標(biāo)的原始數(shù)值為基礎(chǔ)，通過標(biāo)準(zhǔn)化處理消除量綱，得出標(biāo)準(zhǔn)化數(shù)值，以此作為指標(biāo)體系評(píng)價(jià)標(biāo)準(zhǔn)。采用極差變換法對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。

正向指標(biāo)：yij=xij-bijaij-bij；逆向指標(biāo)：yij=aij-xijaij-bij。

式中：yij為無量綱處理后的數(shù)值；xij指j年第i項(xiàng)指標(biāo)的實(shí)際值；aij為第i項(xiàng)指標(biāo)的最大值；bij為第i項(xiàng)指標(biāo)的最小值。考慮到后期繪圖便于觀察，將標(biāo)準(zhǔn)化處理數(shù)值轉(zhuǎn)化為指標(biāo)得分：Qi=yij×10+100。指標(biāo)得分按照正向指標(biāo)和逆向指標(biāo)分別計(jì)算，得分與耕地利用效益均呈正相關(guān)，即Qi越大，耕地利用效益越高，反之亦然。

研究數(shù)據(jù)來源于2006—2015年《黑龍江省統(tǒng)計(jì)年鑒》[17]和各地方統(tǒng)計(jì)年鑒，部分?jǐn)?shù)據(jù)來源于我國(guó)經(jīng)濟(jì)與社會(huì)發(fā)展統(tǒng)計(jì)數(shù)據(jù)庫(kù)和黑龍江省統(tǒng)計(jì)公報(bào)，評(píng)價(jià)方法使用的部分?jǐn)?shù)據(jù)是根據(jù)各項(xiàng)指標(biāo)由原始數(shù)據(jù)計(jì)算整理得到。endprint

2.4耕地利用效益評(píng)價(jià)模型

根據(jù)各指標(biāo)數(shù)據(jù)標(biāo)準(zhǔn)化處理后計(jì)算得分和指標(biāo)綜合權(quán)重，采用多因素綜合評(píng)價(jià)法建立耕地利用效益評(píng)價(jià)模型。P=∑13i=1QiWi；E=∑4i=1QiBi；S=∑9i=5QiBi；F=∑13i=10QiBi。其中，P為綜合效益指數(shù)；E為經(jīng)濟(jì)效益指數(shù)；S為社會(huì)效益指數(shù)；F為生態(tài)效益指數(shù)；Qi為第i個(gè)指標(biāo)的得分值；Wi為各指標(biāo)的組合權(quán)重值，∑Wi=1；Bi為單因子權(quán)重，∑Bi=1。

3結(jié)果與分析

3.1三江平原耕地利用效益時(shí)序特征分析

由圖1可知，2005—2014年三江平原耕地利用社會(huì)效益、經(jīng)濟(jì)效益和綜合效益的時(shí)序特征大致相同，曲線重合度較高，呈波動(dòng)上升趨勢(shì)，2010、2013年是2個(gè)轉(zhuǎn)折點(diǎn)；而生態(tài)效益曲線呈下降趨勢(shì)，且下降幅度較大。

耕地利用經(jīng)濟(jì)效益是耕地綜合利用效益的核心，經(jīng)濟(jì)效益的提高有助于社會(huì)需求滿足度的提高和耕地生態(tài)環(huán)境的改善。三江平原耕地利用經(jīng)濟(jì)效益在2005—2014年這10年內(nèi)增長(zhǎng)幅度較大，2010、2013年出現(xiàn)轉(zhuǎn)折點(diǎn)，均比上一年小幅度下降。分析影響三江平原耕地利用經(jīng)濟(jì)效益的影響因素可以看出，糧食單位面積總產(chǎn)值、農(nóng)業(yè)機(jī)械化程度始終處于增長(zhǎng)態(tài)勢(shì)，技術(shù)效率的進(jìn)步釋放多余勞動(dòng)力，使得單位面積勞動(dòng)力逐漸減少。由于2010年三江平原腹地遭遇百年一遇的早春低溫寡照雪災(zāi)致使農(nóng)時(shí)約延后15 d，嚴(yán)重影響糧食產(chǎn)量，2013年8月份，受降水偏多和主要江河水位上漲的影響，局部農(nóng)田出現(xiàn)澇情，糧食生產(chǎn)受到直接威脅，致使2010、2013年糧食播面單產(chǎn)下降，這2個(gè)年份的耕地利用經(jīng)濟(jì)效益出現(xiàn)轉(zhuǎn)折。

根據(jù)耕地利用社會(huì)效益評(píng)價(jià)指數(shù)變化，2005—2014年三江平原耕地利用社會(huì)效益基本呈穩(wěn)步上升趨勢(shì)，間或有小的回落。隨著社會(huì)經(jīng)濟(jì)的發(fā)展，人均糧食產(chǎn)量、社會(huì)需求滿足度、人均農(nóng)業(yè)國(guó)內(nèi)生產(chǎn)總值、農(nóng)民人均純收入和勞動(dòng)力轉(zhuǎn)移指數(shù)均有所增加；三江平原地廣人稀且人口增長(zhǎng)較為緩慢，人均耕地面積基本趨于穩(wěn)定。另外，由于自然災(zāi)害，2010年、2013年人均糧食產(chǎn)量與社會(huì)需求滿足度有所回落，導(dǎo)致三江平原耕地利用社會(huì)效益出現(xiàn)小幅度下降。

耕地利用的生態(tài)效益關(guān)系到耕地利用的可持續(xù)發(fā)展，2005—2014年三江平原的耕地利用生態(tài)效益的評(píng)價(jià)值與經(jīng)濟(jì)、社會(huì)和綜合效益的相反，出現(xiàn)了明顯下降又趨于平穩(wěn)的態(tài)勢(shì)。2005—2009年，三江平原耕地資源生態(tài)環(huán)境狀況相對(duì)較好，前期生態(tài)效益指數(shù)明顯高于經(jīng)濟(jì)效益、社會(huì)效益和綜合效益，但糧食增產(chǎn)的同時(shí)伴隨著化肥、農(nóng)藥投入量的加大，致使這一階段的耕地利用生態(tài)效益逐年下降。2009—2014年，三江平原耕地生態(tài)效益評(píng)價(jià)值趨于平穩(wěn)，萬元產(chǎn)值能耗的逐步降低緩解了耕地利用生態(tài)效益下降的態(tài)勢(shì)，這種態(tài)勢(shì)從一定意義上表明耕地利用不可持續(xù)的危險(xiǎn)有所控制。

三江平原耕地利用的綜合效益總體呈上升趨勢(shì)，2005—2007年生態(tài)效益的下降一定程度上影響了耕地利用綜合效益，加之經(jīng)濟(jì)效益和社會(huì)效益提升緩慢，2007年耕地利用綜合效益達(dá)到10年間最低水平。2008年以來，三江平原加快推進(jìn)黑龍江、松花江、烏蘇里江三江連通工程，實(shí)現(xiàn)江河、航運(yùn)和水系貫通，發(fā)揮三江平原水利工程綜合效益，隨著水利設(shè)施的修建、農(nóng)業(yè)科技水平的改善與農(nóng)民農(nóng)業(yè)素養(yǎng)的提高，其耕地利用綜合效益在整體上也處于上升態(tài)勢(shì)。

3.2三江平原耕地利用效益空間變異

基于三江平原2014年縣域耕地利用效益指數(shù)，利用ArcGIS 10.0軟件的專題圖顯示功能對(duì)三江平原耕地利用效益屬性數(shù)據(jù)生成空間分布圖，研究經(jīng)濟(jì)效益、社會(huì)效益、生態(tài)效益和綜合效益在地域空間上的分布差異（圖2）。

三江平原耕地利用經(jīng)濟(jì)效益指數(shù)具有明顯的空間集聚特征，松花江沿線平原區(qū)的佳木斯和雙鴨山地區(qū)受到周邊城市群經(jīng)濟(jì)發(fā)展的輻射，整體上經(jīng)濟(jì)利用效益較好，而三江平原的南部地區(qū)，由于人口密集，耕地面積較少，耕地利用經(jīng)濟(jì)效益指數(shù)相對(duì)較低。經(jīng)濟(jì)效益評(píng)價(jià)值得分小于104.00的縣（市、區(qū)）有8個(gè)，經(jīng)濟(jì)效益評(píng)價(jià)值得分在104.01～105.00之間的縣（市、區(qū)）有8個(gè)，經(jīng)濟(jì)效益評(píng)價(jià)值得分在105.01～106.99之間的縣（市、區(qū)）有7個(gè)，表明各區(qū)域之間的耕地利用經(jīng)濟(jì)效益差異較大，各縣（市、區(qū)）經(jīng)濟(jì)耕地利用經(jīng)濟(jì)效益程度參差不齊。三江地區(qū)各縣（市、區(qū)）糧食播面單產(chǎn)平均值為 6.21 t/hm2，大于平均水平的縣（市、區(qū)）有11個(gè)，2014年三江平原單位面積總產(chǎn)值平均值為 21 821.91元/hm2，大于平均水平的縣（市、區(qū)）有9個(gè)，均占45%左右。農(nóng)業(yè)機(jī)械化程度大于平均值的縣（市、區(qū)）有7個(gè)，且空間分布差異較大，這主要是受地形和農(nóng)業(yè)現(xiàn)代化水平的影響。

三江平原耕地利用社會(huì)效益評(píng)價(jià)值得分高于103.00的地區(qū)主要集中在三江平原北部的富錦市、同江市和撫遠(yuǎn)縣，以及三江平原邊緣區(qū)的虎林市和依蘭縣，這些縣（市、區(qū)）地廣人稀，人均耕地面積交廣，人均糧食產(chǎn)量較高，耕地利用社會(huì)效益評(píng)價(jià)值得分也相對(duì)較高。三江平原人均耕地面積約為全國(guó)平均水平的5倍，各縣（市、區(qū)）平均值為0.53 hm2，人均糧食產(chǎn)量平均值為3 164.49 kg，超過平均值的縣（市、區(qū)）有8個(gè)，人均農(nóng)業(yè)國(guó)內(nèi)生產(chǎn)總值和農(nóng)民人均純收入平均值分別為11 514.33、10 244.72 元， 高于平均水平的縣（市、 區(qū)）分別有9、13個(gè)，分別占40%、56%。表明三江平原水土資源豐富，農(nóng)業(yè)產(chǎn)值相對(duì)較高，但農(nóng)民收入差距較大，收入高的農(nóng)民集中在佳木斯、鶴崗、雙鴨山、雞西和七臺(tái)河等5個(gè)市轄區(qū)及其周邊縣（市、區(qū)），這些地區(qū)的勞動(dòng)力轉(zhuǎn)移指數(shù)也比較高。

三江平原耕地利用生態(tài)效益評(píng)價(jià)值得分為小于105.00、105.01～107.00、107.01～109.13的縣（市、區(qū)）分別有6、10、7個(gè)。高值區(qū)主要分布在東部平原區(qū)，耕地利用的生態(tài)效益評(píng)價(jià)值得分大于107.01的有撫遠(yuǎn)縣、虎林市、穆棱市、友誼縣、雞東縣和樺南縣，這些地區(qū)土壤肥沃、農(nóng)業(yè)現(xiàn)代化程度高，化肥、農(nóng)藥的利用率高，對(duì)水體、土壤污染較小。低值區(qū)則集中在七臺(tái)河及其周邊地區(qū)，影響這些縣（市、區(qū)）生態(tài)效益的主要因素是耕地灌溉面積小、化肥施用量大以及能耗相對(duì)較高，其中勃利縣、樺川縣、集賢縣、七臺(tái)河市市轄區(qū)和同江市位于后5位。但從整體上看，三江平原地區(qū)積極響應(yīng)國(guó)家發(fā)展農(nóng)業(yè)建設(shè)高水平農(nóng)業(yè)示范區(qū)的號(hào)召，政府加大農(nóng)業(yè)水利等基礎(chǔ)設(shè)施投入，耕地有效灌溉面積比率逐年上升，這也使三江平原耕地利用生態(tài)效益下降趨勢(shì)有所緩解，耕地利用可持續(xù)度提高。endprint

三江平原耕地利用綜合效益評(píng)價(jià)值呈東西高、中間低的空間分布形態(tài)，這一方面受到地形地理差異的影響，另一方面是由于各縣（市、區(qū)）農(nóng)業(yè)經(jīng)濟(jì)和技術(shù)發(fā)展水平的差異所致。撫遠(yuǎn)縣的耕地利用綜合效益評(píng)價(jià)值最高，為106.93，富錦市、虎林市和雞東縣位于第二梯隊(duì)，這些耕地利用綜合效益評(píng)價(jià)值較高的地區(qū)主要位于三江平原的東北部，該地區(qū)耕地開發(fā)時(shí)間較短，耕地質(zhì)量較高，再加上黑龍江省農(nóng)墾系統(tǒng)歸屬其中，機(jī)械化程度較高，整體上耕地利用綜合效益較好。三江平原耕地利用綜合效益評(píng)價(jià)值得分小于103.00的有七臺(tái)河市市轄區(qū)和勃利縣，其余縣（市、區(qū)）耕地利用綜合效益評(píng)價(jià)值得分均在103.01～105.00之間，表明三江平原大部分縣（市、區(qū)）的耕地利用綜合效益均處于中等水平，還有較大的提升空間。

4結(jié)論與討論

三江平原耕地利用經(jīng)濟(jì)效益、社會(huì)效益和綜合效益在時(shí)間序列上曲線重合度較高，表現(xiàn)為波動(dòng)上升，期間2010、2013年遭受自然災(zāi)害，耕地利用綜合效益評(píng)價(jià)值出現(xiàn)小幅度回落，生態(tài)效益評(píng)價(jià)值在2005—2009年持續(xù)下降，2009年以后變化態(tài)勢(shì)逐漸趨于平穩(wěn)。10年間，耕地利用效益的提高主要依賴于糧食產(chǎn)量、農(nóng)業(yè)機(jī)械化程度、農(nóng)民人均純收入的增長(zhǎng)與萬元產(chǎn)值能耗的降低。

三江平原耕地利用效益空間變異較大，呈明顯的東西高、中間低的空間分布形態(tài)，各縣（市、區(qū)）的農(nóng)業(yè)生產(chǎn)水平不均衡發(fā)展，耕地利用效益高值區(qū)集中在三江平原東北部低平原區(qū)，耕地利用綜合效益較高的地區(qū)有撫遠(yuǎn)市、富錦市、虎林市和雞東縣。對(duì)耕地利用效益空間差異影響較大的指標(biāo)有農(nóng)業(yè)機(jī)械化程度、萬元產(chǎn)值能耗、農(nóng)民人均純收入，表明三江平原耕地利用技術(shù)、勞動(dòng)力資源的區(qū)域分配不平衡和農(nóng)民的貧富差距較大。

未來三江平原在大力發(fā)展農(nóng)業(yè)的同時(shí)更應(yīng)重視耕地利用綜合效益的提高，根據(jù)三江平原各縣（市、區(qū)）耕地利用效益的空間差異，制定相應(yīng)的區(qū)域規(guī)劃，引導(dǎo)各區(qū)域經(jīng)濟(jì)、社會(huì)和生態(tài)效益協(xié)調(diào)發(fā)展。依靠科技進(jìn)步增加糧食播面單產(chǎn)，提高農(nóng)業(yè)生產(chǎn)效率，減小農(nóng)業(yè)機(jī)械化水平和能耗的空間差異，在耕地利用過程中控制農(nóng)藥、化肥施用量的增長(zhǎng)，增加耕地有效灌溉面積，提高耕地資源的有效利用率。

運(yùn)用ArcGIS技術(shù)可以直觀地反映三江平原各縣（市、區(qū)）耕地利用效益水平的相對(duì)空間分布，可為建立耕地利用效益動(dòng)態(tài)監(jiān)測(cè)平臺(tái)提供理論依據(jù)。此外，由于客觀條件的局限性，在今后的研究中，還須進(jìn)一步優(yōu)化耕地利用效益評(píng)價(jià)指標(biāo)體系，減小評(píng)價(jià)結(jié)果與實(shí)際值的偏差。

參考文獻(xiàn)：

[1]劉瓊峰，李明德，段建南，等. 基于GIS的湖南省耕地利用效益時(shí)空變異研究[J]. 經(jīng)濟(jì)地理，2013，33（9）：142-147.

[2]李俏，谷國(guó)鋒，姚麗，等. 黑龍江墾區(qū)耕地利用效率變化分析[J]. 干旱區(qū)資源與環(huán)境，2016，30（4）：30-35.

[3]李瑞華，李明秋. 河南省耕地利用效益時(shí)空分異及驅(qū)動(dòng)因素研究[J]. 地域研究與開發(fā)，2010，29（1）：109-112.

[4]劉吉平，趙丹丹，田學(xué)智，等. 1954—2010年三江平原土地利用景觀格局動(dòng)態(tài)變化及驅(qū)動(dòng)力[J]. 生態(tài)學(xué)報(bào)，2014，34（12）：3234-3244.

[5]王娟，于鳳榮，孫曉兵，等. 三江平原耕地后備資源評(píng)價(jià)及空間特征分析[J]. 現(xiàn)代化農(nóng)業(yè)，2016（5）：59-62.

[6]宋開山，劉殿偉，王宗明，等. 1954年以來三江平原土地利用變化及驅(qū)動(dòng)力[J]. 地理學(xué)報(bào)，2008，63（1）：93-104.

[7]Todkari G U，Suryawanshi S P，Suryawanshi M V，et al. Agriculture landuse pattern in solapur district of Maharashtra[J]. International Journal of Agriculture Sciences，2010，2（2）：1-8.

[8]Herdt A R，Byeongsu Kim A，Taton T A. Encapsulated magnetic nanoparticles as supports for proteins and recyclable biocatalysts[J]. Bioconjugate Chemistry，2007，18（1）：183-189.

[9]Vitousek P M，Naylor R，Crews T，et al. Nutrient imbalances in agricultural development[J]. Science，2009，324（5934）：1519-1520.

[10]傅伯杰，陳利頂，馬克明. 黃土丘陵區(qū)小流域土地利用變化對(duì)生態(tài)環(huán)境的影響——以延安市羊圈溝流域?yàn)槔齕J]. 地理學(xué)報(bào)，1999，54（3）：241-246.

[11]黃秋昊，蔡運(yùn)龍. 國(guó)內(nèi)幾種土地利用變化模型述評(píng)[J]. 中國(guó)土地科學(xué)，2005，19（5）：25-30.

[12]林庶民. 耕地利用效益時(shí)空差異：基于山東省的實(shí)證研究[J]. 內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2007，9（2）：97-100.

[13]牛海鵬，張安錄. 耕地利用生態(tài)社會(huì)效益測(cè)算方法及其應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26（5）：316-323.

[14]陳美球，黃靚，蔡海生，等. 鄱陽湖區(qū)土地健康評(píng)價(jià)[J]. 自然資源學(xué)報(bào)，2004，19（2）：170-175.

[15]宋戈，梁海鷗，林佳，等. 黑龍江省墾區(qū)耕地利用綜合效益評(píng)價(jià)及驅(qū)動(dòng)力分析[J]. 經(jīng)濟(jì)地理，2010，30（5）：835-840.

[16]李豐玉，董子銘. 基于層次分析法（AHP）的休閑農(nóng)業(yè)產(chǎn)業(yè)集群競(jìng)爭(zhēng)力評(píng)價(jià)指標(biāo)體系[J]. 江蘇農(nóng)業(yè)科學(xué)，2014，42（12）：484-486.

[17]黑龍江省統(tǒng)計(jì)局，國(guó)家統(tǒng)計(jì)局黑龍江調(diào)查總隊(duì). 黑龍江統(tǒng)計(jì)年鑒2006—2015[M]. 北京：中國(guó)統(tǒng)計(jì)出版社，2006—2015.江蘇農(nóng)業(yè)科學(xué)2017年第45卷第23期周鴻文，孫艷偉，呂文星，等. 基于SWAT模型的青海引黃灌區(qū)耗水系數(shù)模擬[J]. 江蘇農(nóng)業(yè)科學(xué)，2017，45（23）：248-251.endprint