晏蔚楠，夏建國，2，陳星燎，熊冰瑤

（1.四川農(nóng)業(yè)大學(xué)管理學(xué)院，四川成都611130；2.四川農(nóng)業(yè)大學(xué)資源學(xué)院，四川成都611130；3.四川農(nóng)業(yè)大學(xué)農(nóng)學(xué)院，四川成都611130）

耕地是人類賴以生存的基本資源和條件，是構(gòu)成糧食綜合生產(chǎn)能力最基本的要素之一，其存量與質(zhì)量是保證我國糧食安全的重要方面［1－2］。耕地的保護也由數(shù)量管理向數(shù)量、質(zhì)量并重管理發(fā)展［3］，使得如何加強耕地質(zhì)量保護、穩(wěn)定提高產(chǎn)能成為近年來的研究熱點。耕地產(chǎn)能通過理論產(chǎn)能、可實現(xiàn)產(chǎn)能和實際產(chǎn)能3種層次的表達，衡量了耕地的質(zhì)量以及耕地的“隱形”價值［4］，對其研究不僅能摸清耕地的數(shù)量、質(zhì)量狀況，同時也可以為國家實施耕地保護戰(zhàn)略、確保糧食安全提供依據(jù)。

國外對耕地產(chǎn)能研究主要從作物生產(chǎn)潛力入手，采用生產(chǎn)力數(shù)值模型預(yù)測作物產(chǎn)量潛力進而對土地生產(chǎn)力進行評價［5－6］。而我國自1996年開展農(nóng)用地分等工作以來，普遍采用基于農(nóng)用地分等成果的產(chǎn)能核算方法［7－9］，產(chǎn)能成果也被廣泛應(yīng)用于基本農(nóng)田劃定［10－12］、耕地保護級別評定［13］、耕地占補平衡、土地空間管制［14－15］、土地整理工作以及潛力測算等方面［16－17］。但耕地產(chǎn)能多以市縣為單位，集中于東中部區(qū)域，西南部省域范圍內(nèi)的空間分異特征研究相對較少，而省域尺度是國家糧食安全與耕地保護關(guān)注的重點［18］，因此本研究結(jié)合四川省實際情況，分析四川省空間分異特征和集聚情況，并探討區(qū)域耕地產(chǎn)能的主導(dǎo)因子，以期為提高耕地產(chǎn)能、制訂差異化的耕地保護策略，實現(xiàn)省域耕地可持續(xù)發(fā)展和空間宏觀調(diào)控提供依據(jù)。

1 研究區(qū)域概況

四川省地處長江上游，介于 97°21′～108°33′E、26°03′～34°19′N之間，南北跨度916 km，東西跨度為 1 062 km，處于第1級青藏高原和第2級長江中下游平原的過渡帶，地形呈西高東低的特點。氣候條件差異顯著，東部屬于亞熱帶季風(fēng)氣候區(qū)，西部屬于高原高寒氣候區(qū)。區(qū)域內(nèi)農(nóng)業(yè)利用歷史悠久，耕地資源質(zhì)量較好，耕地集中分布于7～11等。

結(jié)合四川省氣候、地形地貌特征，將研究區(qū)分為5個二級指標(biāo)區(qū)，分別是成都平原區(qū)、盆地丘陵低山區(qū)、盆周山地區(qū)、川西南山地區(qū)和川西北高原山地區(qū)，耕地總面積3 960 448 hm2（表 1）。

2 研究方法與數(shù)據(jù)

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

研究數(shù)據(jù)主要來源于2013年四川省農(nóng)用地分等更新調(diào)查數(shù)據(jù)庫、2011—2013年四川省農(nóng)業(yè)統(tǒng)計年鑒、2013年四川省統(tǒng)計年鑒和各縣理論產(chǎn)量樣點資料。

2.2 研究方法

2.2.1 耕地產(chǎn)能核算 3種層次的產(chǎn)能采用基于耕地分等與農(nóng)業(yè)統(tǒng)計的產(chǎn)能核算方法［7，19］，其中耕地理論產(chǎn)能以農(nóng)用地分等成果、各縣樣地資料為基礎(chǔ)，按《農(nóng)用地質(zhì)量分等規(guī)程》中標(biāo)準(zhǔn)糧換算系數(shù)和標(biāo)準(zhǔn)耕作制度，將指定作物單產(chǎn)轉(zhuǎn)化為標(biāo)準(zhǔn)糧單產(chǎn)。以二級指標(biāo)區(qū)為單位，建立樣本標(biāo)準(zhǔn)糧單產(chǎn)與相應(yīng)地塊的自然等指數(shù)函數(shù)：

式中：yi表示第i個單元標(biāo)準(zhǔn)糧理論單產(chǎn)；Ri表示第i個分等單元自然等指數(shù)；a、b表示回歸系數(shù)值。

表1 四川省二級指標(biāo)區(qū)耕地狀況

按照理論單產(chǎn)模型，將二級指標(biāo)內(nèi)的耕地自然等指數(shù)帶入模型，獲取理論單產(chǎn)，進而核算出耕地理論產(chǎn)能：

式中：WFi表示耕地理論產(chǎn)能；YFi表示耕地理論單產(chǎn)；Si表示第i個分等單元耕地面積。

耕地可實現(xiàn)單產(chǎn)將3年內(nèi)指定作物的最高產(chǎn)量核算成標(biāo)準(zhǔn)糧產(chǎn)量，構(gòu)建標(biāo)準(zhǔn)糧產(chǎn)量與利用等指數(shù)間的函數(shù)：

式中：y′i表示第i個單元標(biāo)準(zhǔn)糧可實現(xiàn)單產(chǎn)；R′i表示第i個分等單元利用等指數(shù)；c、d表示回歸系數(shù)值。

依據(jù)可實現(xiàn)產(chǎn)能核算模型，將二級指標(biāo)區(qū)的利用等指數(shù)帶入模型，獲取耕地可實現(xiàn)單產(chǎn)，然后，核算出耕地可實現(xiàn)產(chǎn)能：

式中：W′Fi表示耕地可實現(xiàn)產(chǎn)能；Y′Fi表示耕地可實現(xiàn)單產(chǎn)；Si表示第i個分等單元耕地面積。

耕地實際產(chǎn)能以2013年農(nóng)業(yè)統(tǒng)計年鑒為準(zhǔn)，將二級指標(biāo)區(qū)內(nèi)的指定作物換算為標(biāo)準(zhǔn)糧，作為耕地實際單產(chǎn)，用耕地實際單產(chǎn)乘以耕地面積即為耕地實際產(chǎn)能。

2.2.2 GIS－ESDA分析法 GIS－ESDA分析法是指在整體把握耕地產(chǎn)能空間分布狀況的情況下，進一步發(fā)現(xiàn)空間集聚和空間異質(zhì)現(xiàn)象，為功能分區(qū)劃分和宏觀調(diào)控提供依據(jù)。其中ESDA技術(shù)的核心內(nèi)容是空間自相關(guān)測度。本研究以全局Morans I指數(shù)來進行測算全局空間自相關(guān)，Morans I計算公式如下：

式中：N表示空間單元數(shù)目為觀測值，為Xi的平均值。W（i，j）為研究范圍內(nèi)空間單元i與空間單元j（i，j＝1，2，3，…，n）的空間連接矩陣。本研究所使用的空間連接矩陣為Rook方式，一階鄰接關(guān)系，即以1和0表示i與j的相鄰關(guān)系，1表示i與j相鄰，0表示不相鄰，定義W（i，j）＝0，依此得到一個N維的矩陣W（n，n）。

同時采用LISA統(tǒng)計指標(biāo)來描述局部區(qū)域具有一定相似程度的數(shù)值在空間上的聚集程度，其計算公式為：

式中：N、Xi、Xj、W（i，j）含義同公式（5）。

Local Morans I采用條件隨機方法獲得偽顯著水平P值進行判斷，并形成LISA聚類圖，可視化的表達區(qū)域上高值聚集（H－H型）、低值聚集（L－L型）、高低聚集（H－L型）和低高聚集（L－H型）4種聚集模式。

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

基礎(chǔ)數(shù)據(jù)處理采用Excel 2010，耕地產(chǎn)能空間分布采用ArcGIS軟件完成，耕地產(chǎn)能空間集聚采用OpenGeoDa軟件完成，影響因子分析采用SPSS軟件完成。

3 結(jié)果與分析

3.1 耕地產(chǎn)能數(shù)量分異特征

按照上述耕地產(chǎn)能核算方法，核算出四川省耕地產(chǎn)能（表 2），結(jié)果表明：

（1）耕地單產(chǎn)最高區(qū)域位于成都平原區(qū)、盆地丘陵低山區(qū)，2個二級指標(biāo)區(qū)耕地單產(chǎn)均值相差較小，其中盆地丘陵低山區(qū)耕地單產(chǎn)最大。成都平原區(qū)的理論單產(chǎn)、可實現(xiàn)單產(chǎn)和實際單產(chǎn)分別比盆地丘陵低山區(qū)低 2.1、0.36、082 t／hm2。可見，成都平原區(qū)、盆地丘陵低山區(qū)作為四川省傳統(tǒng)農(nóng)業(yè)資源最優(yōu)區(qū)域，耕地單產(chǎn)高，區(qū)域農(nóng)業(yè)發(fā)展的自然條件和農(nóng)業(yè)配套設(shè)施投入差距小，耕地單產(chǎn)均值也相差較小；但在市場條件作用下，成都平原區(qū)耕地等別略高于盆地丘陵低山區(qū)，但區(qū)域耕地多用于種植經(jīng)濟作物，而盆地丘陵低山區(qū)作為四川省糧、經(jīng)、豬、工貿(mào)綜合大區(qū)，耕地仍多用于種植水稻、玉米等糧食作物，因此盆地丘陵低山區(qū)的耕地單產(chǎn)大于成都平原區(qū)。

（2）盆周山地區(qū)、川西南山地區(qū)單產(chǎn)均值相差較小，但與成都平原區(qū)和盆地丘陵低山區(qū)相比，各層次單產(chǎn)相差較大，差值范圍介于5～10 t／hm2。主要是盆周山地區(qū)和川西南山地區(qū)作為平壩、淺丘向深丘、山地的過渡帶，氣候生產(chǎn)潛力、地形條件、經(jīng)濟發(fā)展水平和農(nóng)業(yè)配套設(shè)施差距小，但相較于成都平原區(qū)和盆地丘陵低山區(qū)差距還較大。

（3）耕地單產(chǎn)最低的區(qū)域位于川西北高原山地區(qū)，區(qū)域理論單產(chǎn)、可實現(xiàn)單產(chǎn)和實際單產(chǎn)與盆地丘陵低山區(qū)相比，相差 25.20、23.09、20.88 t／hm2。一方面是因為川西北的地形以高原、山地為主，海拔高，光溫生產(chǎn)潛力低，農(nóng)業(yè)利用水平差；另一方面是區(qū)域內(nèi)種植的糧食作物多以玉米、馬鈴薯、青稞為主，且青稞僅種植于海拔3 000 m的地方，按照農(nóng)用地分等規(guī)程中四川省標(biāo)準(zhǔn)糧核算方式來看，馬鈴薯并未納入核算范圍。

表2 四川省二級指標(biāo)區(qū)各層次耕地產(chǎn)能單產(chǎn)

（4）從耕地總產(chǎn)看，不同區(qū)域耕地總產(chǎn)差距較大，成都平原區(qū)與盆地丘陵低山區(qū)相比，理論總產(chǎn)、可實現(xiàn)總產(chǎn)和實際總產(chǎn)分別相差 5 728.89×104、5 069.73×104、4 627.21×104t，耕地產(chǎn)能的變化差異與耕地面積變化一致（圖1），可見保有相對數(shù)量的耕地面積是保證糧食安全的基礎(chǔ)。而對比耕地理論總產(chǎn)、可實現(xiàn)總產(chǎn)和實際總產(chǎn)（圖1）發(fā)現(xiàn)，除盆地丘陵低山區(qū)外，其余指標(biāo)區(qū)3種層次耕地總產(chǎn)相差較小，說明四川省耕地整體利用強度大，而盆地丘陵低山區(qū)作為糧食綜合大區(qū)耕地產(chǎn)能還有較大的提升空間。

3.2 耕地產(chǎn)能空間分異特征

為進一步分析省域尺度不同層次耕地產(chǎn)能的空間分異特征，結(jié)合耕地產(chǎn)能數(shù)量，采用ArcGIS，對省域尺度范圍內(nèi)不同層次的耕地單產(chǎn)進行普通克里格空間插值，按自然斷點法劃分5個［20－21］等級，得到四川省不同層次耕地單產(chǎn)分布圖（圖2）和總產(chǎn)分布圖（圖3）。

耕地單產(chǎn)整體呈東高西低分布，其理論單產(chǎn)和可實現(xiàn)單產(chǎn)的分異特征和等級區(qū)間范圍基本一致，實際單產(chǎn)分布更為復(fù)雜；西部與理論單產(chǎn)、可實現(xiàn)單產(chǎn)的分布特征基本一致，但在東部成都平原區(qū)和盆地丘陵低山區(qū)則呈環(huán)狀分布，從成都市主城區(qū)向外實際單產(chǎn)不斷增加，介于 20 920.91～69 329.89 kg／hm2之間（圖2）。耕地單產(chǎn)東高西低的分布特征與氣候、地形及農(nóng)業(yè)利用水平相關(guān)；而實際單產(chǎn)東西部分布差異，可能是因為西部地區(qū)由于自然條件和交通的限制，受社會經(jīng)濟和農(nóng)業(yè)利用水平影響較小，與理論單產(chǎn)和可實現(xiàn)單產(chǎn)分異基本一致，而東部地區(qū)受城市擴張影響，實際單產(chǎn)呈現(xiàn)環(huán)狀分布的特點。

耕地總產(chǎn)的分異特征與耕地單產(chǎn)相似，表現(xiàn)為東部高、西部低的特點，理論總產(chǎn)和可實現(xiàn)總產(chǎn)分異基本一致，實際總產(chǎn)分布更為復(fù)雜（圖3）。實際總產(chǎn)在川西南山地區(qū)的分布南北差異較大，米易縣、仁和區(qū)、東區(qū)、西區(qū)、西昌市的耕地實際總產(chǎn)高于川西南山地區(qū)的其余區(qū)縣，介于 71.02×104～100.02×104t之間，可見在實際農(nóng)業(yè)生產(chǎn)過程中，自然條件相對優(yōu)越的河谷地農(nóng)業(yè)利用水平更高，耕地產(chǎn)能潛力更大。

3.3 耕地產(chǎn)能空間自相關(guān)性

GIS分析是一種松散的空間分布［21］，為進一步了解不同層次耕地產(chǎn)能的集中、分散程度，更好地分析省域耕地產(chǎn)能的空間分異特征，本研究采用ESDA方法揭示四川省耕地產(chǎn)能的空間相互作用，為統(tǒng)籌區(qū)域發(fā)展策略提供依據(jù)。

3.3.1 全局空間自相關(guān)分析 利用OpenGeoDa空間分析軟件，采用Rook一階鄰接定義權(quán)重矩陣，通過999次隨機檢驗，獲得偽顯著性P＜0.05的不同層次耕地產(chǎn)能全局Morans I指數(shù)值（表3）。

結(jié)果表明：四川省耕地單產(chǎn)、耕地總產(chǎn)的分布具有較強的空間正相關(guān)性，空間集聚程度由大到小分別是理論產(chǎn)能＞可實現(xiàn)產(chǎn)能＞實際產(chǎn)能，表明受社會經(jīng)濟活動的影響，耕地單產(chǎn)的空間集聚程度在減弱。對比耕地單產(chǎn)和耕地總產(chǎn)的Global Morans I指數(shù)可見，3種層次耕地總產(chǎn)的 Global Morans I指數(shù)值小于耕地單產(chǎn)，說明省域耕地產(chǎn)能的空間集聚性受耕地面積影響較大。

3.3.2 局部空間自相關(guān)分析 基于各縣（市、區(qū)）不同層次耕地產(chǎn)能核算結(jié)果，結(jié)合ArcGIS 10.0生成縣（市、區(qū)）中心點圖，利用OpenGeoDa的空間相關(guān)分析得到四川省耕地產(chǎn)能的LISA聚類圖，揭示了省域耕地產(chǎn)能局部空間的分異特征（表4、圖4）。

表3 四川省不同層次耕地產(chǎn)能全局Morans I指數(shù)值

表4 局部空間自相關(guān)類型縣級個數(shù)統(tǒng)計個

結(jié)果表明：（1）省域耕地產(chǎn)能的集聚效性較低，且耕地總產(chǎn)的集聚性小于耕地單產(chǎn)。從表4中可以看出，集聚區(qū)縣占全省的比例均小于50%，理論單產(chǎn)所占比例最大，為45.60%，實際單產(chǎn)所占比例最小，為25.27%，說明四川省農(nóng)業(yè)發(fā)展的協(xié)調(diào)性較差。（2）從省域集聚分布（圖4）看，L－L聚集區(qū)呈現(xiàn)較強的一致性，集中分布于西部川西北高原山地區(qū)。H－H聚集區(qū)呈現(xiàn)較大差異，耕地總產(chǎn)集中分布于盆地丘陵低山區(qū)；耕地單產(chǎn)分布復(fù)雜，理論單產(chǎn)分布在盆地丘陵低山區(qū)、成都平原區(qū)的涪城區(qū)、游仙區(qū)、江陽區(qū)、廣漢市和盆周山地區(qū)的興文縣、萬源市；可實現(xiàn)單產(chǎn)分布在盆地丘陵低山區(qū)、成都平原區(qū)的成都市三圈層、涪城區(qū)、游仙區(qū)、羅江縣、江陽區(qū)、廣漢市；實際單產(chǎn)分布在成都平原區(qū)的溫江區(qū)、雙流區(qū)、龍泉驛區(qū)以及盆地丘陵低山區(qū)的金堂縣、簡陽市、射洪縣、嘉陵區(qū)等地。L－H聚集區(qū)差異更大，實際單產(chǎn)分布在威遠(yuǎn)縣和榮縣；理論總產(chǎn)和可實現(xiàn)總產(chǎn)分布在通川區(qū)、華鎣市、羅江縣、彭山縣和江陽區(qū)；實際總產(chǎn)分布在通川區(qū)、順慶區(qū)、涪城區(qū)、雙流區(qū)和青神縣。通過對比LISA聚類圖（圖4）與耕地產(chǎn)能分布圖（圖2、圖3）得出，四川省耕地產(chǎn)能的空間集中分布帶來了集聚效應(yīng)，但由于下墊面和人為活動的影響，空間集中分布對H－H型聚集效應(yīng)影響較小。

3.4 耕地產(chǎn)能主導(dǎo)因子分析

參照農(nóng)用地質(zhì)量分等規(guī)程中的評價因子，結(jié)合四川省實際情況，本研究選取了14個參評因子分析影響耕地單產(chǎn)的主導(dǎo)因子。選取的自然因素有：X1代表年均氣溫，X2代表年降水量，X3代表日照時數(shù)，X4代表地形坡度，X5代表土壤質(zhì)地，X6代表地表裸露程度，X7代表土層厚度，X8代表土壤有機質(zhì)含量；社會經(jīng)濟因素有：X9代表化肥施用量，X10代表農(nóng)業(yè)機械總動力，X11代表有效灌溉面積，X12代表排水條件，X13代表人均收入，X14代表農(nóng)村勞動力人口，探討了二級指標(biāo)區(qū)各參評因子與不同層次耕地單產(chǎn)的相關(guān)性（表5），結(jié)果表明：

（1）盆西平原區(qū)理論單產(chǎn)與年降水量、地表裸露程度和排水條件顯著相關(guān)；可實現(xiàn)單產(chǎn)與年平均氣溫、年降水量和人均收入顯著相關(guān)；實際單產(chǎn)則與年降水量、地形坡度、日照時數(shù)和農(nóng)村勞動人口、化肥施用量、有效灌溉面積顯著相關(guān)。3種層次耕地單產(chǎn)均與年降水量顯著相關(guān)，說明保證水源是盆西平原區(qū)提高耕地產(chǎn)能的必要條件。

（2）盆地丘陵低山區(qū)理論單產(chǎn)與土層厚度、土壤有機質(zhì)含量和年均降水量顯著相關(guān)；可實現(xiàn)單產(chǎn)與地表裸露程度、土壤有機質(zhì)含量和年降水量顯著相關(guān)；實際單產(chǎn)與土壤質(zhì)地、有機質(zhì)含量、日照時數(shù)、年均氣溫、農(nóng)村勞動人口、農(nóng)業(yè)機械化、化肥施用量、有效灌溉面積顯著相關(guān)。3種層次耕地單產(chǎn)均與有機質(zhì)含量顯著相關(guān)，說明提高土壤養(yǎng)分是盆地丘陵低山區(qū)提升耕地產(chǎn)能的重要因素，而農(nóng)村勞動人口與實際單產(chǎn)顯著相關(guān)，則說明盆地丘陵低山區(qū)農(nóng)業(yè)現(xiàn)代化發(fā)展程度較低。

（3）盆周山地區(qū)理論單產(chǎn)與地形坡度極顯著相關(guān)；可實現(xiàn)單產(chǎn)與排水條件、地形坡度、農(nóng)業(yè)機械化、化肥使用量、有效灌溉面積以及年均降水量顯著相關(guān)；實際單產(chǎn)與排水條件、地形坡度、農(nóng)村機械、化肥施用量、年均降水相關(guān)性強。3種層次耕地單產(chǎn)均與地形坡度顯著相關(guān)，說明地形坡度是盆周山地區(qū)農(nóng)業(yè)發(fā)展的主要限制因素。

（4）川西南山地區(qū)耕地理論單產(chǎn)與排水條件、日照時數(shù)和年均氣溫顯著相關(guān)；可實現(xiàn)單產(chǎn)與排水條件、地形坡度、日照時數(shù)、年均氣溫、人均收入顯著相關(guān)；實際單產(chǎn)則與土層厚度、人均收入、農(nóng)業(yè)機械、化肥施用量和有效灌溉面積顯著相關(guān)。

（5）川西北高原山地區(qū)理論單產(chǎn)與年均氣溫顯著相關(guān)；可實現(xiàn)單產(chǎn)與土壤質(zhì)地、化肥施用量和年均氣溫相關(guān)性強；實際單產(chǎn)則與年均氣溫、土層厚度、農(nóng)業(yè)機械總動力、化肥使用量和有效灌溉面積顯著相關(guān)。3種層次耕地單產(chǎn)均與氣溫顯著相關(guān)，主要是因為川西北高原山地區(qū)垂直地帶性明顯，氣溫差異顯著，因此發(fā)展川西北高原山地區(qū)農(nóng)業(yè)重點考慮海拔和氣候條件。

4 結(jié)論

四川省不同區(qū)域耕地產(chǎn)能數(shù)量差異大，二級指標(biāo)區(qū)耕地產(chǎn)能由大到小為：盆地丘陵低山區(qū)＞成都平原區(qū)＞盆周山地區(qū)＞川西南山地區(qū)＞川西北高原山地區(qū)。但耕地單產(chǎn)成都平原區(qū)和盆地丘陵低山區(qū)、盆周山地區(qū)和川西南山地區(qū)兩兩對比單產(chǎn)均值相差較??；耕地總中的理論總產(chǎn)、可實現(xiàn)總產(chǎn)和實際總產(chǎn)，除盆地丘陵低山區(qū)外3種層次數(shù)量差異較小，說明四川省耕地利用強度大，盆地丘陵低山區(qū)是發(fā)展農(nóng)業(yè)的潛力區(qū)域。

四川省耕地產(chǎn)能總體分布特征是東部高于西部，耕地理論產(chǎn)能與可實現(xiàn)產(chǎn)能分布規(guī)律基本一致，耕地實際產(chǎn)能分布復(fù)雜。實際產(chǎn)能西部川西北高原山地區(qū)與理論產(chǎn)能和可實現(xiàn)產(chǎn)能分布基本一致，大體處于一個等級區(qū)間；東部盆地丘陵地山區(qū)與成都平原區(qū)呈環(huán)狀分布，川西南山地區(qū)中的米易縣、仁和區(qū)、東區(qū)、西區(qū)、西昌市的實際總產(chǎn)高于東南部其他區(qū)縣，盆周山地區(qū)分異特征不明顯。

四川省各層次耕地產(chǎn)能均存在比較明顯的空間正相關(guān)性，空間集聚程度由大到小分別是理論產(chǎn)能＞可實現(xiàn)產(chǎn)能＞實際產(chǎn)能，但省域耕地產(chǎn)能集聚效應(yīng)較低，且耕地總產(chǎn)的集聚性小于耕地單產(chǎn)。集聚分布表現(xiàn)為：L－L聚集區(qū)呈現(xiàn)較強的一致性，H－H聚集區(qū)和L－H聚集區(qū)呈現(xiàn)較大差異。四川省耕地產(chǎn)能的空間集中分布帶來了集聚效應(yīng)，但由于下墊面和人為活動的影響，空間集中分布對H－H型聚集效應(yīng)影響較小。

四川省不同區(qū)域影響耕地產(chǎn)能的主導(dǎo)因子不同，3種層次耕地單產(chǎn)在成都平原區(qū)、盆地丘陵低山區(qū)、盆周山地區(qū)、川西北高原分別主要受降水、土壤有機質(zhì)、地形坡度、氣溫的影響。5個二級指標(biāo)區(qū)耕地實際單產(chǎn)均與農(nóng)業(yè)機械化、化肥施用量、有效灌溉面積顯著相關(guān)，因此在實際農(nóng)業(yè)生產(chǎn)過程，提高耕地產(chǎn)能要重視農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè)和農(nóng)業(yè)機械化發(fā)展。

本研究從宏觀上對四川省耕地理論產(chǎn)能、可實現(xiàn)產(chǎn)能和實際產(chǎn)能的分異特征進行了分析，同時對不同指標(biāo)區(qū)提升耕地產(chǎn)能的主導(dǎo)因子進行了探討，但就耕地產(chǎn)能分異特征的原因未進行系統(tǒng)分析，還需深入研究。

參考文獻：

［1］封志明，李香蓮.耕地與糧食安全戰(zhàn)略：藏糧于土，提高中國土地資源的綜合生產(chǎn)能力［J］.地理學(xué)與國土研究，2000，16（3）：1－5.

［2］孔祥斌.區(qū)域耕地數(shù)量控制與質(zhì)量提升實踐［M］.北京：科學(xué)出版社，2014.

［3］楊邦杰，鄖文聚，程 鋒.論耕地質(zhì)量與產(chǎn)能建設(shè)［J］.中國發(fā)展，2012，12（1）：1－6.

［4］伍育鵬，鄖文聚，鄒 如.基于耕地質(zhì)量評價的耕地產(chǎn)能核算方法與分布［J］.農(nóng)業(yè)工程學(xué)報，2008，24（增刊1）：85－89.

［5］Acemoglu D，Johnson S，Robinson J.Institutions as the fundamental cause of long－run growth［M］.［S.l.］：NBER Working Paper，2004.

［6］Miller S，Upadhyay M.Total factor productivity and the convergence hypothesis［J］.Journal of Management，2002，24（2）：267－286.

［7］鄖文聚，王洪波，王國強，等.基于農(nóng)用地分等與農(nóng)業(yè)統(tǒng)計的產(chǎn)能核算研究［J］.中國土地科學(xué)，2007，21（4）：32－37.

［8］錢鳳魁，王秋兵，邊振興.遼寧省農(nóng)用地產(chǎn)能核算研究［J］.土壤通報，2014，45（3）：543－549.

［9］喻 婧.基于RS和農(nóng)用地分等的武穴市耕地生產(chǎn)能力核算研究［D］.武漢：華中師范大學(xué)，2014.

［10］潘洪義，蔣貴國，何 偉.基于農(nóng)用地產(chǎn)能核算成果基本農(nóng)田劃定研究——以安縣為例［J］.中國農(nóng)學(xué)通報，2012，28（8）：160－165.

［11］湯思陽，雷國平，張 慧，等.基于耕地質(zhì)量評價的基本農(nóng)田劃定方法［J］.水土保持研究，2013，20（1）：173－176，181.

［12］孔祥斌，靳 京，劉 怡.基于農(nóng)用地利用等別的基本農(nóng)田保護區(qū)劃定［J］.農(nóng)業(yè)工程學(xué)報，2008，24（10）：46－51.

［13］劉華賓，張俊梅，門明新，等.基于產(chǎn)能的耕地資源保護級別評定方法研究——以河北省盧龍縣為例［J］.自然資源學(xué)報，2009，24（8）：1325－1333.

［14］曹小曙，李 濤，黃曉燕.耕地綜合產(chǎn)能指數(shù)的構(gòu)建及其實踐意義［J］.廣東土地科學(xué)，2013，12（1）：25－29.

［15］陳亞恒，張 利，張俊梅，等.基于糧食生產(chǎn)能力的耕地非農(nóng)可轉(zhuǎn)量測算模型［J］.自然資源學(xué)報，2013，28（3）：450－458.

［16］劉巧芹，李子君，吳克寧，等.中國耕地整理潛力測算方法研究綜述［J］.資源開發(fā)與市場，2013，29（2）：127－130，148.

［17］林和明，肖北生，陳旭飛，等.粵西南丘陵山區(qū)耕地占補平衡中的質(zhì)量與產(chǎn)能變化研究［J］.廣東農(nóng)業(yè)科學(xué)，2013，40（12）：202－206，215.

［18］張紅富，周生路，吳紹華，等.省域尺度耕地產(chǎn)能空間分異規(guī)律及其影響因子［J］.農(nóng)業(yè)工程學(xué)報，2010，26（8）：308－314.

［19］中華人民共和國國土資源部.農(nóng)用地質(zhì)量分等規(guī)程：TD／T 1004—2012［S］，2012：18－19.

［20］許 妍，吳克寧，程先軍，等.東北地區(qū)耕地產(chǎn)能空間分異規(guī)律及產(chǎn)能提升主導(dǎo)因子分析［J］.資源科學(xué)，2011，33（11）：2030－2040.

［21］李子良，王樹濤，張 利，等.經(jīng)濟快速發(fā)展地區(qū)耕地生產(chǎn)能力空間格局［J］.農(nóng)業(yè)工程學(xué)報，2010，26（11）：323－331.

［22］馬榮華，黃杏元，朱傳耿.用ESDA技術(shù)從GIS數(shù)據(jù)庫中發(fā)現(xiàn)知識［J］.遙感學(xué)報，2002，6（2）：102－107.