馬才學，金 瑩，柯新利，馬艷春

（華中農業(yè)大學公共管理學院，湖北 武漢 430070）

耕地是人類生存發(fā)展至關重要的生產要素，也是糧食生產中不可替代的資源。隨著人們對物質需求的多樣化，耕地功能日益豐富，耕地多功能特征逐步體現(xiàn)。耕地資源作為土地資源的重要組成，在食物供給、生態(tài)環(huán)境維持和社會穩(wěn)定中發(fā)揮著重要功效[1]，同時耕地的休閑和娛樂功能也受到人們的重視[2]。

自20世紀80年代提出農業(yè)多功能以來，農業(yè)多功能性受到了各界認可，從最初局限于自由貿易的探討延伸至農業(yè)政策手段、多功能指標建立并實際運用于國家政策的調整等[3-6]，耕地的多功能性也得到了國內外學者普遍關注。Y Takatori以東京練馬區(qū)為例研究城市農業(yè)多功能發(fā)現(xiàn)農業(yè)多功能可劃分為農業(yè)生產功能、環(huán)境維護功能、災害防御功能、休閑娛樂功能及土地供給功能5類[7]。姜廣輝等通過大量文獻梳理及定量分析法提出耕地多功能可分為基本和衍生兩個層次功能，應以耕地多功能的層次性作為指導制定耕地保護工作方案[8]。宋小青等采用政策反推的耕地功能分析方法構建指標體系，并運用熵權法對1949年以來中國耕地功能進行評價[9]。施園園等通過信息熵模型對北京13個區(qū)縣（不包括中心城區(qū)）耕地多功能與社會經濟的空間相關性進行測算，討論耕地多功能與其經濟社會發(fā)展的互動關系以及發(fā)展階段的協(xié)調模式[10]。

但當前研究也存在不足：（1）層次分析法等指標表征法因其系統(tǒng)體現(xiàn)了相關評價的綜合性、整體性與層次性，以及評價過程與結果的簡單明了、易于公眾認知等優(yōu)點得到廣泛應用，但易受主觀判斷和選擇對結果產生影響[11]；（2）多功能綜合強度的定量評價往往采用多功能指標加權求和來反映，但不同功能指標間高低值的相互補償使得亟待重視的多功能均衡發(fā)展問題被忽略[12]。

因此，本文通過構建耕地多功能評價指標體系，引入全排列多邊形圖示指標法測度耕地多功能實現(xiàn)情況。該評價模型簡單直觀，能有效克服傳統(tǒng)主觀判斷權重的不足，且引入幾何圖形面積和形狀指數(shù)公式對耕地多功能的綜合強度和協(xié)調度進行測算，從耕地多功能強度和協(xié)調度兩方面進行評價，實現(xiàn)耕地多功能評價與可視化表達的集成。通過探究2015年湖北省各縣（區(qū)）耕地多功能發(fā)展現(xiàn)狀，識別耕地多功能強度與協(xié)調度典型模式，厘清湖北省耕地多功能典型模式的空間分布，為耕地多功能利用提供科學依據(jù)。

1 研究區(qū)域及數(shù)據(jù)來源

1.1 研究區(qū)域

本文以湖北省為研究區(qū)域。湖北省位于中國中部，長江中游，其所在的江漢平原是中國重要的糧食生產基地之一。全省復雜多樣的地形地貌形態(tài)，使得耕地利用類型復雜多樣，為調整農業(yè)產業(yè)結構，發(fā)揮耕地多功能性提供了前提條件。截至2015年末，湖北省耕地總資源達3960.48×103hm2，列全國第八，人均耕地面積為1.02畝，低于中國平均水平。2009—2015年，全省耕地累計減少6.3×104hm2，年均減少耕地1.05×104hm2，城鎮(zhèn)化和工業(yè)化快速發(fā)展對建設用地的需求與農業(yè)生產對耕地的保護之間矛盾加劇。隨著湖北省實施的“兩圈一帶”發(fā)展戰(zhàn)略，以“兩圈”為引擎、“一帶”為傳動，所帶來的全省經濟平穩(wěn)較快發(fā)展的同時，耕地資源管理面臨著更大的壓力。

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

本文的數(shù)據(jù)主要包括湖北省各縣（區(qū)）耕地、人口、產業(yè)發(fā)展、就業(yè)和農業(yè)生產等數(shù)據(jù)。人口、產業(yè)發(fā)展、經濟等數(shù)據(jù)來源于《湖北省統(tǒng)計年鑒》（2016年）；耕地、農業(yè)生產、農村人口就業(yè)等數(shù)據(jù)來源于《湖北農村統(tǒng)計年鑒》（2016年）；其他部分數(shù)據(jù)由湖北省各縣域統(tǒng)計年鑒獲得。由于部分區(qū)縣耕地數(shù)量較少和數(shù)據(jù)獲取影響，最終確定湖北省各區(qū)縣耕地多功能評價單元共包括秭歸、嘉魚、竹溪等79個縣（區(qū)）。

2 研究方法

2.1 耕地多功能評價指標體系

已有學者開展了耕地多功能研究并構建相應的評價體系，但耕地功能分類尚未統(tǒng)一。耕地功能即耕地系統(tǒng)提供滿足人類生存與發(fā)展所需產出的能力[13-14]，從系統(tǒng)論角度來看，耕地是由自然生態(tài)子系統(tǒng)和社會經濟子系統(tǒng)耦合而成，耕地功能除包括糧食生產等傳統(tǒng)商品性功能外，還包括調節(jié)氣候、提供休閑觀光的綠色景觀等生態(tài)功能，保持傳統(tǒng)農耕文化、為農民提供社會保障的社會功能以及建設空間儲備等承載功能的非商品性生產功能[15]。參考宋小青采用的政策反推的耕地功能分析方法[9]，結合中國城鄉(xiāng)二元化土地制度的特殊背景，將耕地功能劃分為糧食安全保障、就業(yè)保障、社會安定維護、國民經濟貢獻、家庭經濟貢獻及生態(tài)恢復6類（表1）。

表1 耕地多功能評價指標體系Tab.1 Multi-functional evaluation index system of farmland

2.2 全排列多邊形圖示指標法

目前全排列多邊形圖示指標法應用于生態(tài)城市[17]、土地集約利用[18]、綠地生態(tài)服務效能[19]等多方面評價。其具體定義為[20]：以n個評價指標（標準化）的上限值為半徑構成一個中心正N邊形，指標值連線構成一個不規(guī)則中心N邊形，頂點是n個指標的一個首尾相接的全排列。

指標值標準化采用雙曲線標準化函數(shù)Q（x）：

式（1）中，L、U、T為指標的下限、上限和臨界閾值。Q（x）滿足：Q（x）|x=L= -1，Q（x）|x=T= 0，Q（x）|x=U= 1。標準化函數(shù)Q（x）把位于區(qū)間［L，U］的指標值映射到［-1，+1］區(qū)間。指標由原始的線性變化調整為快—慢—快的非線性變化，臨界值為指標變化速度的轉折點。

因此，第i個單項功能的標準化公式為：式（2）中，fi為指標標準化值，Ni是指標實際值，Ui、Li、Ti分別為指標值的上限值、下限值和閾值。n個指標構成一個中心正n邊形，n個頂點為fi= 1時的值，中心點為fi= -1時的值，中心點到頂點的線段為各指標標準化值所在區(qū)間［-1，+1］，而fi= 0時構成的多邊形為指標的臨界區(qū)。

耕地單項功能值計算公式：

當功能指標值等于指標上限時，耕地功能值達到理論最大值1；當功能指標值等于指標下限時，耕地功能值為理論最小值0，即耕地功能值在［0，1］區(qū)間。借鑒相關研究[21-22]，確定分級評價標準：功能指標值F∈［0，0.33］，功能等級為低度持續(xù)；功能指標值F∈（0.33，0.66］，功能等級為中度持續(xù)；功能指標值F∈（0.66，1］，功能等級為高度持續(xù)。

2.3 耕地多功能綜合強度及協(xié)調度測算

引入多邊形面積（S）表征區(qū)域耕地多功能實現(xiàn)的綜合強度，是該區(qū)域耕地多功能評價構成的全排列多邊形的面積。當所有功能指標值均達上限時，耕地多功能綜合強度達理論最大值S = nsin（360°/n）/2。S越大，表明耕地多功能實現(xiàn)程度越大。

耕地多功能綜合強度函數(shù)計算公式：

引入景觀形狀指數(shù)（L）表征區(qū)域耕地各功能實現(xiàn)的協(xié)調程度，是該區(qū)域耕地多功能評價構成的全排列多邊形與相同面積的正n邊形之間的偏離程度。當所有功能指標值大小相等時，耕地多功能協(xié)調度達到理論最小值1，因此耕地多功能協(xié)調度L≥1。L越大，表明耕地各功能實現(xiàn)程度差異越大，越趨于單一化，反之越趨于多元化。耕地多功能協(xié)調度函數(shù)L（x）：

2.4 耕地多功能的情景模式聚類

應用系統(tǒng)聚類分析法，根據(jù)不同評價單元的耕地6大單項功能發(fā)展水平的組合， 對79個研究對象的耕地多功能評價結果進行聚類分析， 并結合耕地多功能強度和協(xié)調度，分類得到耕地多功能典型模式。

3 結果與分析

3.1 湖北省耕地單項功能空間差異分析

將湖北省79個縣（區(qū)）耕地單項功能值結果作為屬性數(shù)據(jù)，根據(jù)功能等級的劃分規(guī)則，繪制出6大單項功能分級圖（圖1），從空間分析的角度探討耕地各單項功能的空間分布特征。

3.1.1 糧食安全保障功能

由圖1可以看出，位于鄂中部地區(qū)的棗陽市等12個縣（區(qū)）耕地的糧食安全保障功能較高，為高度持續(xù)等級。鄂中部地區(qū)位于長江中下游平原，土壤肥沃、灌溉保證率高、農業(yè)生產條件優(yōu)越，并且湖北省的農產品主產區(qū)也主要分布在該區(qū)域。這12個縣（區(qū)）的耕地普遍具有面積廣、墾殖率高、產量多的特點。而糧食安全保障功能較低的縣（區(qū)）一部分位于鄂西的山區(qū)地帶，由于地形因素，該區(qū)域耕地面積較少且多為坡耕地，生產力較低；另一部分位于湖北省省會武漢市周邊，雖也位于長江中下游平原，但由于接近湖北省經濟發(fā)展中心，該區(qū)域縣（區(qū)）從行政區(qū)劃來看屬于城市發(fā)展新區(qū)，人口密集，建設用地為主要用地類型，耕地面積較小。

3.1.2 就業(yè)保障功能

位于鄂西部地區(qū)的神農架林區(qū)等10個縣（區(qū)）耕地的就業(yè)保障功能較高，為高度持續(xù)等級。這與近年來湖北省實施的“兩圈一帶”發(fā)展戰(zhàn)略有密不可分的關系，湖北省著重在鄂西部地區(qū)打造鄂西生態(tài)文化旅游圈，與武漢城市發(fā)展圈及湖北長江經濟帶相比，該區(qū)域更注重耕地等生態(tài)用地的保育工作。且鄂西部地區(qū)人口總數(shù)相對較少，農業(yè)勞動力占人口總數(shù)比例較高，因此該區(qū)域耕地的就業(yè)保障功能相對較高。而湖北省省會武漢市周邊、省域副中心城市宜昌市、襄陽市周邊等區(qū)域為城市發(fā)展新區(qū)，著重于城市發(fā)展工作，非農勞動力占人口總數(shù)比例相對較高，因此該區(qū)域耕地的就業(yè)保障功能相對較低。

圖1 湖北省耕地單項功能空間分異圖Fig.1 Spatial pattern of farmland function in Hubei Province

3.1.3 社會安定維護功能

該功能空間分異情況呈明顯區(qū)分：鄂中部地區(qū)近40個縣（區(qū)）耕地的社會安定維護功能較高，為高度持續(xù)等級；鄂東南地區(qū)的上10個縣（區(qū)）耕地的社會安定維護功能次之，為中度持續(xù)等級；鄂西部及鄂東北地區(qū)的30多個縣（區(qū)）耕地的社會安定維護功能較低，為低度持續(xù)等級。鄂中部地區(qū)由于較好的區(qū)位條件，例如較優(yōu)的交通運輸條件，使得該區(qū)域農產品市場輻射力強，農產品銷售市場活躍。此外，該區(qū)域人均耕地面積較大且市場環(huán)境良好使農民能獲得穩(wěn)定持續(xù)的收入。鄂西部地區(qū)處于山區(qū)，人均耕地面積較其他區(qū)縣略小，且困難的交通運輸條件使得該區(qū)域城鄉(xiāng)居民收入差距較大。

3.1.4 國民經濟貢獻功能

位于鄂西部地區(qū)的鄖西縣等多個縣（區(qū)）耕地的國民經濟貢獻功能較高，為高度持續(xù)等級。這些區(qū)域的地區(qū)生產總值與其他縣（區(qū)）相比較低，湖北省作為全國糧食主產地之一，各縣（區(qū)）的農業(yè)產值都較高，因此這些區(qū)域的農業(yè)產值占地區(qū)生產總值的比重較大。而屬于湖北省農產品主產區(qū)的縣（區(qū)）其農業(yè)產值在湖北省所有縣（區(qū)）排行中名列前茅，但其地區(qū)生產總值也在湖北省所有縣（區(qū)）排行中排名靠前，因此這些區(qū)域的農業(yè)產值占地區(qū)生產總值的比重次之，為中度持續(xù)等級。

3.1.5 家庭經濟貢獻功能

耕地的家庭經濟貢獻功能較高的區(qū)域大多位于“鄂西生態(tài)文化旅游圈”中，由于農業(yè)觀光園等休閑旅游產業(yè)在這個區(qū)域的發(fā)展，一方面提高了農戶農業(yè)收入水平，另一方面相對較為不利的區(qū)位條件使得該區(qū)域農村勞動力選擇非農活動的機會較低，因此農戶純收入主要來源于農業(yè)收入， 耕地的家庭經濟貢獻為高度持續(xù)等級。而屬于湖北省農產品主產區(qū)的區(qū)域由于農業(yè)生產的逐漸現(xiàn)代化帶來的生產設備的改良，不僅使得農戶收入增加、農業(yè)勞動力流失較小，也使得農戶有更多的精力選擇非農活動，因此該區(qū)域農戶收入結構更為多元化，耕地的家庭經濟貢獻功能為中度持續(xù)等級。

3.1.6 生態(tài)恢復功能

耕地的生態(tài)恢復功能較高的區(qū)域相對比較集中在鄂中的農產品主產區(qū)，這些區(qū)域農業(yè)基礎設施更為先進，與其他區(qū)域相比先一步實現(xiàn)了農業(yè)的現(xiàn)代化、都市化。并且該區(qū)域土壤肥沃、灌溉保證率高、農業(yè)生產條件優(yōu)越的自然條件使得該區(qū)域耕地農藥化肥施用強度相對較低，因此其耕地的生態(tài)恢復功能為高度持續(xù)等級。 耕地的生態(tài)恢復功能較低的區(qū)域相對比較集中于湖北省的山區(qū)或丘陵地區(qū)，由于該區(qū)域坡耕地較多，現(xiàn)代化的農業(yè)技術難以在這些地區(qū)實施，因此該區(qū)域仍舊保留原始的資本密集型生產方式使得該區(qū)域耕地農藥化肥施用強度相對較高，因此其耕地的生態(tài)恢復功能為低度持續(xù)等級。

3.2 湖北省各縣（區(qū)）耕地多功能評價

計算得到湖北省各縣（區(qū)）2015年的耕地多功能綜合強度和協(xié)調度，并將結果用ArcGIS進行空間分析（圖2）。

3.2.1 耕地多功能強度由鄂中向鄂西、鄂東弱化

研究以各縣（區(qū)）耕地單項功能值作為指標值構建出全排列多邊形，并測算其面積作為衡量耕地多功能強度的標準。由圖2可知，鄂中部地區(qū)耕地多功能綜合強度普遍較高，耕地多功能強度較強，鄂西、鄂東部地區(qū)耕地多功能綜合強度普遍較低，耕地多功能強度較弱?？傮w上來看，耕地多功能強度由鄂中向鄂西、鄂東弱化。

3.2.2 耕地多功能協(xié)調度由鄂東向鄂西呈減弱趨勢

研究以各縣（區(qū)）耕地單項功能值作為指標值構建出全排列多邊形，并測算其與相同面積的正多邊形之間的偏離程度作為衡量耕地多功能協(xié)調度的標準。由圖2可知，鄂東、鄂中部地區(qū)耕地多功能協(xié)調度普遍較高，耕地多功能協(xié)調度較好，鄂東部地區(qū)耕地多功能協(xié)調度普遍較低，耕地多功能協(xié)調度較差?？傮w上來看，耕地多功能協(xié)調度由鄂東向鄂西呈減弱趨勢。

圖2 湖北省各縣（區(qū)）耕地多功能綜合強度與協(xié)調度Fig.2 Comprehensive value and coordination degrees of farmland in Hubei Province

3.3 湖北省各縣（區(qū)）耕地多功能典型模式識別

耕地多功能的影響因素較多，加上區(qū)域自身發(fā)展狀況及其發(fā)展目標等因素差異，導致耕地多功能表現(xiàn)模式差異顯著。對各縣（區(qū)）耕地多功能的表現(xiàn)模式進行識別及分析，可以探究適宜湖北省發(fā)展戰(zhàn)略的耕地多功能發(fā)展模式。鑒于此，本文根據(jù)湖北省各縣（區(qū)）耕地多功能的差異特征進行系統(tǒng)聚類分析，將湖北省79個縣（區(qū)）分為5類耕地多功能發(fā)展模式（表2）。

湖北省各縣（區(qū)）的自然條件、經濟社會發(fā)展實力等因素具有明顯差異性，導致不同地區(qū)耕地多功能結構存在顯著差異性，從而形成了多樣化的耕地多功能模式，并呈現(xiàn)出明顯的區(qū)域差異性（圖3）。具體表現(xiàn)為由鄂西至鄂中、鄂東地區(qū)，耕地多功能模式逐漸從生態(tài)優(yōu)先發(fā)展模式轉向多功能協(xié)調發(fā)展模式及多功能強協(xié)調發(fā)展模式，最終轉變?yōu)槎喙δ苋鯀f(xié)調發(fā)展模式。該研究結果與目前湖北省實施的“兩圈一帶”發(fā)展戰(zhàn)略所確定的區(qū)域功能定位相符：鄂東部地區(qū)處于武漢市城市圈范圍，是城市功能拓展區(qū)，耕地多功能綜合強度較弱但各單項功能協(xié)同發(fā)展特征明顯；鄂中部地區(qū)是湖北省農產品主產區(qū)聚集區(qū)域，良好的自然條件、日益現(xiàn)代化的農業(yè)設備，使得該區(qū)域耕地各功能持續(xù)等級均較高并且差異較小，多功能綜合強度較強；鄂西部地區(qū)處于鄂西生態(tài)文化旅游圈范圍，是生態(tài)涵養(yǎng)發(fā)展區(qū)，生態(tài)功能為該區(qū)域的主導功能，其他單項功能發(fā)展滯后，耕地多功能結構相對失衡。

表2 湖北省耕地多功能典型模式Tab2 The typical mode of farmland in Hubei Province

圖3 湖北省耕地多功能模式現(xiàn)狀圖Fig.3 Status of farmland multi-functional modes in Hubei Province

4 結論與討論

研究基于耕地多功能的概念及政策反推的耕地功能分析方法，構建耕地多功能評價體系，引入全排列多邊形圖示指標法，從耕地多功能的綜合強度及協(xié)調度兩方面評價了湖北省各縣（區(qū)）耕地利用多功能情況，并對評價結果聚類得出湖北省各縣（區(qū)）的耕地多功能典型模式。

（1）研究提出的耕地多功能綜合強度及協(xié)調度二維評價模型，能夠具體、客觀地反映研究區(qū)耕地多功能情況及空間分異規(guī)律。結果顯示，耕地單項功能發(fā)展水平及多功能綜合水平在各縣（區(qū)）呈現(xiàn)明顯的差異，總體上表現(xiàn)出東、中、西的區(qū)域分異現(xiàn)狀。耕地多功能綜合強度由鄂中向鄂西、鄂東弱化，耕地多功能協(xié)調度由鄂東向鄂西呈減弱趨勢。湖北省耕地多功能綜合強度與協(xié)調度評價存在顯著差異，這一結果進一步說明了多功能綜合功能的定量評價往往只能顯示出耕地多功能性的一部分，想要更加全面地對耕地多功能性進行評價，引入?yún)f(xié)調度的測定是必要的。

（2）受自然環(huán)境及社會經濟發(fā)展水平在湖北省各縣（區(qū)）空間差異影響，研究區(qū)耕地多功能強度與協(xié)調程度關系的典型模式可劃分為5類：多功能強協(xié)同發(fā)展模式、多功能協(xié)同發(fā)展模式、多功能弱協(xié)同發(fā)展模式、生態(tài)優(yōu)先發(fā)展模式、社會主導發(fā)展模式。在空間分布上具體表現(xiàn)為由鄂西至鄂中、鄂東地區(qū)，耕地多功能模式逐漸從生態(tài)優(yōu)先發(fā)展模式轉向多功能協(xié)調發(fā)展模式及多功能強協(xié)調發(fā)展模式，最終轉變?yōu)槎喙δ苋鯀f(xié)調發(fā)展模式。該研究結果與目前湖北省實施的“兩圈一帶”發(fā)展戰(zhàn)略所確定的區(qū)域功能定位相符，比較客觀地揭示了處于同一發(fā)展時期不同地區(qū)的耕地多功能現(xiàn)狀，對指導不同地區(qū)耕地多功能優(yōu)化調整具有重要作用。

耕地多功能的提出強調了耕地利用的合理方向，也增強了嚴格保護耕地的說服力和行動能力?，F(xiàn)有大量學者運用功能指標體系法對耕地多功能進行評價，例如楊雪[10]、施園園[22]分別運用灰色關聯(lián)度、層次分析法對北京市耕地多功能演變及空間差異進行了研究，發(fā)現(xiàn)近年來北京市耕地各項功能都在不斷增強，但耕地各項功能演變是否協(xié)調并未提及。而本文提出的耕地多功能綜合強度及協(xié)調度二維評價模型，一方面解決了傳統(tǒng)多功能評價易受主觀判斷和選擇對結果產生影響的確定，另一方面解決了多功能綜合功能的定量評價中對多功能均衡發(fā)展忽視的問題，實現(xiàn)耕地多功能評價與可視化表達的集成。為識別耕地多功能發(fā)展現(xiàn)狀提供合理方法，對指導耕地多功能發(fā)展提供了一些依據(jù)，但要想更好指導耕地多功能發(fā)展，僅有本文的研究還是不足的，還應當進一步研究同一地區(qū)在不同時期耕地多功能的變化規(guī)律和模式。城市化發(fā)展導致耕地多功能結構變化與重組，耕地多功能的優(yōu)化調整應當堅持在可持續(xù)發(fā)展的條件下合理滿足城市化發(fā)展的需求，即以實現(xiàn)城市化與耕地多功能的最佳耦合狀態(tài)為發(fā)展目標。因此，如何根據(jù)地區(qū)城市化發(fā)展特征測度耕地多功能最佳結構模式應當成為下一步研究的重點。

［1］ 羅成，蔡銀鶯. 湖北省農產品主產區(qū)耕地資源功能的時空演變［J］ . 經濟地理，2016，36（3）：153 - 161.

［2］ 楊雪，談明洪. 近年來北京市耕地多功能演變及其關聯(lián)性［J］ . 自然資源學報，2014，29（5）：733 - 743.

［3］ Huang J, Tichit M, Poulot M, et al. Comparative Review of Multifuntionality and Ecosystem Services in Sustainable Agriculture［J］ . Journal of Environmental Management，2015（149）：138 - 147.

［4］ 胡偉艷，魏安奇，趙志尚，等. 農地多功能供需錯位與協(xié)同作用研究進展及趨勢［J］ . 中國土地科學，2017，31（3）：89 - 97.

［5］ Burkhard B, Kroll F, Nedkov S, et al. Mapping Ecosystem Service Supply, Demand and Budgets［J］ . Ecological Indicators，2012（21）：17 - 29.

［6］ Van Huylenbroeck G, Vandermeulen V, Mettepenningen E, et al. Multifunctionality of Agriculture: A Review of Definitions,Evidence and Instruments［J］ . Living Rev.Landscape Res，2007（3）：5 - 31.

［7］ Takatori Y. The Multi-Functional System of Urban Agriculture in the Suburbs of Tokyo-A Case Study in the Nishi-Ooizumi Area of Nerima Ward, Tokyo［J］ . Journal of Geography，2010，109（3）：401 - 417.

［8］ 姜廣輝，張鳳榮，孔祥斌，等. 耕地多功能的層次性及其多功能保護［J］ . 中國土地科學，2011，25（8）：42 - 47.

［9］ 宋小青，吳志峰，歐陽竹. 1949年以來中國耕地功能變化［J］ . 地理學報，2014（4）：435 - 447.

［10］ 施園園，趙華甫，鄖文聚，等. 北京市耕地多功能空間分異及其社會經濟協(xié)調模式解釋［J］ . 資源科學，2015，37（2）：247 - 257.

［11］ 吳瓊，王如松，李宏卿，等. 生態(tài)城市指標體系與評價方法［J］ . 生態(tài)學報，2005，25（8）：2090 - 2095.

［12］ 彭建，劉志聰，劉焱序. 農業(yè)多功能性評價研究進展［J］ .中國農業(yè)資源與區(qū)劃，2014，35（6）：1 - 8.

［13］ 宋小青，歐陽竹. 中國耕地多功能管理的實踐路徑探討［J］ . 自然資源學報，2012，27（4）：540 - 551.

［14］ Renetzeder C, Eupen, Mucher C A, et al. Sustainability Impact Assessment of Land Use Changes［J］ . Sustainability Impact Assessment of Land Use Changes，2008，33（4）：159 - 180.

［15］ 宋小青，歐陽竹. 耕地多功能內涵及其對耕地保護的啟示［J］ . 地理科學進展，2012，31（7）：859 - 868.

［16］ 龍方，曾福生. 中國糧食安全的戰(zhàn)略目標與模式選擇［J］ .農業(yè)經濟問題，2008（7）：32 - 38.

［17］ 張偉，張宏業(yè)，王麗娟，等. 生態(tài)城市建設評價指標體系構建的新方法——組合式動態(tài)評價法［J］ . 生態(tài)學報，2014（16）：4766 - 4774.

［18］ 周偉，曹銀貴，喬陸印. 基于全排列多邊形圖示指標法的西寧市土地集約利用評價［J］ . 中國土地科學，2012，26（4）：84 - 90.

［19］ 王洪威，徐建剛，桂昆鵬，等. 城市綠地系統(tǒng)生態(tài)服務效能評價及優(yōu)化研究——以淮安生態(tài)新城為例［J］ . 環(huán)境科學學報，2012，32（4）：1018 - 1024.

［20］ 張路路，鄭新奇，原智遠，等. 基于全排列多邊形綜合圖示法的唐山市土地利用多功能性評價［J］ . 中國土地科學，2016，30（6）：23 - 32.

［21］ 王楓，董玉祥. 廣州市土地利用多功能的空間差異及影響因素分析［J］ . 資源科學，2015，37（11）：2179 - 2192.

［22］ 楊雪，談明洪. 北京市耕地功能空間差異及其演變［J］ .地理研究，2014，33（6）：1106 - 1118.