谷國政，宋 戈

（東北大學(xué)文法學(xué)院，遼寧 沈陽 110169）

1 引言

耕地資源為人類生存和發(fā)展提供基礎(chǔ)物質(zhì)保障，承載著食物生產(chǎn)、社會(huì)穩(wěn)定、鄉(xiāng)村發(fā)展、農(nóng)耕文明以及生態(tài)保育等重要功能，是區(qū)域高質(zhì)量發(fā)展的基礎(chǔ)性資源和戰(zhàn)略性條件[1-2]。根據(jù)“土地利用行為理論”分析，土地系統(tǒng)發(fā)展方向主要取決于人的直接利用行為，而由其資源特性誘致下的不合理利用及生態(tài)環(huán)境問題又將深刻影響經(jīng)濟(jì)社會(huì)發(fā)展[3]。耕地利用方式、強(qiáng)度、功能及價(jià)值的演變折射了中國鄉(xiāng)村地區(qū)“人地”關(guān)系的動(dòng)態(tài)變化，這也是對(duì)農(nóng)業(yè)、農(nóng)村和農(nóng)民發(fā)展現(xiàn)狀及現(xiàn)實(shí)問題的一種側(cè)面反映。然而，隨著城鄉(xiāng)融合戰(zhàn)略的深入推進(jìn)，城鄉(xiāng)經(jīng)濟(jì)社會(huì)高速發(fā)展背景下區(qū)域空間結(jié)構(gòu)與產(chǎn)業(yè)結(jié)構(gòu)重構(gòu)驅(qū)動(dòng)土地、勞動(dòng)力、資本、技術(shù)、數(shù)據(jù)等要素產(chǎn)生適應(yīng)性調(diào)整和改變，城鎮(zhèn)發(fā)展、農(nóng)業(yè)生產(chǎn)和生態(tài)保護(hù)的空間競爭愈發(fā)強(qiáng)烈。城鄉(xiāng)關(guān)系轉(zhuǎn)型重塑與城市蔓延引發(fā)了耕地非農(nóng)化、非糧化、粗放化、邊際化、逆集約化[4]等問題。相應(yīng)的，耕地資源的稀缺性日益凸顯，公眾對(duì)耕地資源的需求變化推動(dòng)耕地功能從單一生產(chǎn)功能向生產(chǎn)、生活和生態(tài)功能發(fā)展。同時(shí)，不同區(qū)域耕地功能在提供農(nóng)業(yè)生產(chǎn)場所、保障農(nóng)民生活資料、調(diào)節(jié)農(nóng)村生態(tài)環(huán)境等方面存在差異，依據(jù)耕地多功能利用過程中的價(jià)值表現(xiàn)合理配置耕地資源，對(duì)于平衡經(jīng)濟(jì)社會(huì)發(fā)展的空間競爭需求具有積極意義，由此，耕地價(jià)值被拓展為衡量耕地“三生”功能的重要標(biāo)準(zhǔn)[5]。中共二十大報(bào)告提出“全方位夯實(shí)糧食安全根基，全面落實(shí)糧食安全黨政同責(zé)，牢牢守住十八億畝耕地紅線”，為嚴(yán)守耕地紅線、保障國家糧食安全指明了方向。但結(jié)合上述耕地保護(hù)面臨的新形勢(shì)、新問題來看，僅通過數(shù)量保護(hù)與質(zhì)量提升難以應(yīng)對(duì)耕地保護(hù)的重大難題，更難以全面構(gòu)建耕地?cái)?shù)量、質(zhì)量、生態(tài)“三位一體”保護(hù)新格局[6]。對(duì)此，如何基于耕地多功能融合及其價(jià)值增益的協(xié)同管理，以推進(jìn)耕地資源保護(hù)與可持續(xù)利用值得深入研究和思考。

近年來，理論界圍繞耕地功能與價(jià)值研究等取得有益成果。在耕地功能方面，主要集中在內(nèi)涵辨析[1,3]、理論框架[7-8]、指標(biāo)構(gòu)建[9-10]、模式識(shí)別[11]、功能評(píng)價(jià)[2,12]及權(quán)衡/協(xié)同關(guān)系[13]、管理策略[14]等；在耕地價(jià)值方面，主要關(guān)注價(jià)值重建[15]、影響因素[5,16]、核算體系[17]及方法[18-19]、功能價(jià)值感知[20]等。耕地多功能時(shí)空分異特征及功能間互作規(guī)律、耕地資源價(jià)值分布及影響因素等總體態(tài)勢(shì)已然揭示[2,5,12-13,16]，為新時(shí)期耕地保護(hù)研究奠定基礎(chǔ)。眾多研究認(rèn)為耕地價(jià)值源于其自身功能，耕地資源的多功能性決定其多元化價(jià)值趨向[7,14,20]。耕地功能研究實(shí)現(xiàn)了從單一生產(chǎn)功能向生產(chǎn)功能、生態(tài)功能、社會(huì)功能等演進(jìn)[15]，耕地價(jià)值也形成了基于經(jīng)濟(jì)價(jià)值、社會(huì)保障價(jià)值和生態(tài)價(jià)值的內(nèi)涵體系[17]，為耕地功能測度與價(jià)值核算提供了理論基礎(chǔ)[20]。但是，耕地多功能演變與耕地價(jià)值響應(yīng)研究尚處于探索階段，已有成果對(duì)耕地功能與價(jià)值之間響應(yīng)關(guān)系的理論闡釋較弱，亟需開展典型區(qū)域耕地功能與價(jià)值的動(dòng)態(tài)變化研究。同時(shí)，由國際全球環(huán)境變化研究科學(xué)機(jī)構(gòu)IGBP和IHDP共同發(fā)起的全球土地計(jì)劃（GLP）提出將“土地多功能性”作為社會(huì)經(jīng)濟(jì)—自然生態(tài)耦合研究的基礎(chǔ)框架，并倡議推進(jìn)該領(lǐng)域研究以科學(xué)系統(tǒng)地認(rèn)知人類—資源環(huán)境耦合系統(tǒng)的變化機(jī)理與互饋機(jī)制，對(duì)于促進(jìn)生態(tài)環(huán)境治理、實(shí)現(xiàn)資源永續(xù)利用具有積極意義[21]。因此，耕地功能與價(jià)值之間具有哪些關(guān)聯(lián)？其理論依據(jù)是什么？耕地多功能演變究竟如何影響其價(jià)值變化？這些問題亟需回答。

遼寧省是東北地區(qū)沿海經(jīng)濟(jì)大省和農(nóng)業(yè)大省，具備發(fā)展耕地多功能融合與復(fù)合價(jià)值增益的先行條件。但隨著新型城鎮(zhèn)化進(jìn)程加快，全省優(yōu)質(zhì)耕地流失、質(zhì)量退化、生態(tài)安全風(fēng)險(xiǎn)加劇，為糧食可持續(xù)生產(chǎn)帶來威脅。而遼寧省作為糧食主產(chǎn)區(qū)，更應(yīng)充分考慮耕地資源稀缺性，通過耕地多功能利用發(fā)揮其復(fù)合價(jià)值，推進(jìn)區(qū)域土地資源優(yōu)化配置?；谇笆霈F(xiàn)實(shí)背景及相關(guān)研究，本文以遼寧省作為典型研究區(qū)域，嘗試從理論視角辨析耕地功能與價(jià)值的內(nèi)在關(guān)聯(lián)，并對(duì)耕地多功能演變特征及“經(jīng)濟(jì)—社會(huì)—生態(tài)”價(jià)值響應(yīng)過程進(jìn)行實(shí)證分析，以期為深化新時(shí)期耕地保護(hù)轉(zhuǎn)型、完善耕地保護(hù)補(bǔ)償機(jī)制提供支撐。

2 耕地功能與價(jià)值的關(guān)聯(lián)與響應(yīng)

耕地利用系統(tǒng)作為典型的半自然—半人工生態(tài)系統(tǒng)，關(guān)聯(lián)著社會(huì)、經(jīng)濟(jì)、生態(tài)等多個(gè)子系統(tǒng)，承載或孕育的糧食生產(chǎn)、經(jīng)濟(jì)產(chǎn)出、生活保障、生態(tài)維持等功能價(jià)值具有共生性，本質(zhì)上是耕地自身及其外部因素耦合互饋形成的綜合能力[1,13]。依據(jù)耕地資源“二維要素—三重功能—多元價(jià)值”認(rèn)知系統(tǒng)[22]進(jìn)行審視，立地、社會(huì)、經(jīng)濟(jì)、生產(chǎn)、生態(tài)等不同維度的要素是耕地資源的基礎(chǔ)，這些要素以特定比例、分層次耦合形成具有一定秩序的結(jié)構(gòu)序列[23]，并伴隨著不同的產(chǎn)品產(chǎn)出及功能服務(wù)，這是耕地資源保持系統(tǒng)性、整體性、協(xié)同性且具有特定功能和價(jià)值的內(nèi)在根據(jù)（圖1）。

圖1 耕地功能與價(jià)值的內(nèi)在關(guān)聯(lián)Fig.1 Internal relations between cultivated land function and its value

從耕地功能與價(jià)值的內(nèi)在關(guān)聯(lián)來看，生產(chǎn)功能是由耕地的本質(zhì)屬性產(chǎn)生，生態(tài)功能源于耕地的固有屬性，而生活功能依賴于耕地的一般屬性，耕地多重屬性的內(nèi)在關(guān)聯(lián)決定其功能之間的隸屬關(guān)系，最終表現(xiàn)為權(quán)衡/協(xié)同的減損或增益效應(yīng)。同時(shí)，耕地價(jià)值主要源自其稀缺性、有用性和需求性，是居于耕地保護(hù)外部效益核心、難以進(jìn)行市場交易卻又客觀存在的價(jià)值[24]。耕地功能和價(jià)值是耕地利用形態(tài)、人類需求及耕地產(chǎn)出（產(chǎn)品與服務(wù)）的重要紐帶，二者之間的響應(yīng)將深刻影響耕地資源、資產(chǎn)、資本價(jià)值實(shí)現(xiàn)。

耕地資源內(nèi)部要素與結(jié)構(gòu)的層次關(guān)系不斷演化，導(dǎo)致其功能與價(jià)值產(chǎn)生相互依賴、彼此制約的疊加態(tài)勢(shì)。一方面，耕地功能供需受社會(huì)基礎(chǔ)、制度環(huán)境、市場機(jī)制、人類福祉和生物物理過程控制，持續(xù)推動(dòng)不同功能在空間層面呈現(xiàn)此消彼長的動(dòng)態(tài)性變化，使得耕地功能與價(jià)值保持并存且具有空間集聚性分布規(guī)律的狀態(tài)；另一方面，耕地價(jià)值關(guān)聯(lián)耕地資源屬性，是協(xié)調(diào)耕地自然條件、質(zhì)量建設(shè)、產(chǎn)能提升與環(huán)境要素關(guān)系等可持續(xù)利用目標(biāo)層面的內(nèi)在基礎(chǔ)，加之聯(lián)合生產(chǎn)不斷內(nèi)化耕地規(guī)?；?、集約化利用的可持續(xù)性，進(jìn)而使得耕地功能與價(jià)值存在一定的耦合協(xié)調(diào)特征?？梢?，耕地功能與價(jià)值的耦合互饋可以有效反映耕地系統(tǒng)內(nèi)部物質(zhì)、能量、信息交換的結(jié)果及外部效益。例如，在耕地利用實(shí)踐中，無論城鎮(zhèn)化快速發(fā)展地區(qū)還是生態(tài)脆弱地區(qū)，不同程度的土地利用變化均會(huì)影響其地表屬性特征，一旦耕地利用強(qiáng)度逼近資源環(huán)境本底，耕地功能供需失衡容易引發(fā)耕地價(jià)值不協(xié)調(diào)、不匹配、不平衡等問題，進(jìn)而產(chǎn)生耕地細(xì)碎化、“非農(nóng)化”“非糧化”、農(nóng)田生態(tài)系統(tǒng)退化和環(huán)境成本外部化等威脅耕地安全的一系列現(xiàn)實(shí)問題。因此，耕地功能對(duì)耕地價(jià)值的影響是直接性的，價(jià)值對(duì)功能的反饋是間接性的，而且兩者之間存在一對(duì)一、一對(duì)多和多對(duì)多等多維互動(dòng)過程，總體呈現(xiàn)出一種疊加性、連續(xù)性、耦合性且非對(duì)稱的循環(huán)響應(yīng)路徑（圖2）。

圖2 耕地功能與價(jià)值響應(yīng)關(guān)系Fig.2 Response relationship between cultivated land function and its value

綜觀耕地功能與價(jià)值的內(nèi)在關(guān)聯(lián)，結(jié)合二者響應(yīng)關(guān)系的理論意涵，特定的效用或功能是相對(duì)應(yīng)關(guān)聯(lián)價(jià)值的重要源泉，價(jià)值則是功能的外在表現(xiàn)，即耕地功能決定其價(jià)值，而價(jià)值大小取決于各項(xiàng)功能的實(shí)現(xiàn)水平?？傊?，耕地利用系統(tǒng)的構(gòu)成要素、內(nèi)源結(jié)構(gòu)、多重功能等在不同尺度間的時(shí)空轉(zhuǎn)換及耦合互饋共同塑造了耕地價(jià)值特征，表現(xiàn)出一定的順序性、階段性、周期性不平衡等特點(diǎn)。對(duì)此，耕地資源作為農(nóng)業(yè)生產(chǎn)的基本物質(zhì)條件，對(duì)區(qū)域經(jīng)濟(jì)社會(huì)轉(zhuǎn)型以及高質(zhì)量發(fā)展既有直接支撐作用又有間接支撐作用，揭示耕地功能演變對(duì)價(jià)值影響過程，闡明二者之間的作用關(guān)系，有益于完善耕地資源價(jià)值體系、深化新時(shí)期耕地保護(hù)的體制機(jī)制。

3 數(shù)據(jù)來源與研究方法

3.1 研究區(qū)概況

遼寧省位于東北地區(qū)南部，全省共轄14個(gè)地級(jí)市，分別是沈陽、大連、鞍山、撫順、本溪、丹東、錦州、營口、阜新、遼陽、盤錦、鐵嶺、朝陽、葫蘆島。2020年，遼寧省鄉(xiāng)村人口為2 235.5萬人，占全省總?cè)丝诘?2.91%，耕地面積為496.81萬hm2，占全省土地總面積的33.43%，其中80%耕地資源分布于遼中平原區(qū)和遼西北低山丘陵的河谷地帶。2020年全省糧食產(chǎn)量為2 338.8萬t，糧食播種面積為352.72萬hm2，水稻播種面積為52.04萬hm2，小麥播種面積為0.31萬hm2，玉米播種面積為269.93萬hm2，是中國糧食主產(chǎn)區(qū)之一，對(duì)保障國家糧食安全發(fā)揮著“壓艙石”作用。此外，2020年全省農(nóng)業(yè)產(chǎn)值為1 912億元，占第一產(chǎn)業(yè)產(chǎn)值的43.77%，第一產(chǎn)業(yè)從業(yè)人員為711.8萬人，占總就業(yè)人員的31.49%①數(shù)據(jù)來源于《遼寧省統(tǒng)計(jì)年鑒（2020）》，https://tjj.ln.gov.cn/tjj/tjxx/xxcx/tjnj/otherpages/2021/indexch.htm。。

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

本文使用的數(shù)據(jù)來源于2005年、2010年、2015年和2020年《全國農(nóng)產(chǎn)品成本收益資料匯編》、《遼寧省統(tǒng)計(jì)年鑒》與各地市統(tǒng)計(jì)年鑒，以及2004年、2009年、2014年和2019年遼寧省各地市國民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)。其中，農(nóng)業(yè)補(bǔ)貼、糧食市場價(jià)格、城鄉(xiāng)居民社會(huì)養(yǎng)老保險(xiǎn)等數(shù)據(jù)來自遼寧省及各地市相關(guān)職能部門在網(wǎng)上公開公布的政策文件與統(tǒng)計(jì)資料。

3.3 研究方法

3.3.1 評(píng)價(jià)指標(biāo)體系及權(quán)重確定

基于耕地功能與價(jià)值關(guān)聯(lián)的理論分析，結(jié)合相關(guān)文獻(xiàn)[2,7,9,11]和遼寧省耕地資源利用現(xiàn)狀，遵循系統(tǒng)性、代表性、合理性和可行性相統(tǒng)一原則，本文構(gòu)建了耕地多功能評(píng)價(jià)指標(biāo)體系（3項(xiàng)主導(dǎo)類型、12項(xiàng)指標(biāo)層），采用極值標(biāo)準(zhǔn)化與熵權(quán)法確定各指標(biāo)的客觀權(quán)重，具體見表1。

表1 耕地多功能評(píng)價(jià)指標(biāo)體系Tab.1 Evaluation system of cultivated land multifunction

（1）生產(chǎn)功能。主要包括耕地墾殖率、糧食單產(chǎn)、人均糧食保證率和地均種植業(yè)總產(chǎn)值等評(píng)價(jià)指標(biāo)。

（2）生活功能。主要包括人均耕地面積、農(nóng)戶最低生活保障水平、種植業(yè)從業(yè)比重和農(nóng)業(yè)機(jī)械化水平等評(píng)價(jià)指標(biāo)。其中農(nóng)戶最低生活保障水平[12,25]計(jì)算公式如下：

式（1）中：SSg是g市的耕地社會(huì)保障價(jià)值；P1是農(nóng)村人口數(shù)量；G是各市域最低生活保障標(biāo)準(zhǔn)；R1是農(nóng)村人均支配收入；R2是城鎮(zhèn)人均可支配收入。

（3）生態(tài)功能。主要包括地均化肥投施強(qiáng)度、地均農(nóng)藥投施強(qiáng)度、耕地自然固碳能力和農(nóng)田系統(tǒng)生境質(zhì)量等評(píng)價(jià)指標(biāo)。其中，耕地自然固碳能力[26]計(jì)算公式為：

式（2）中：Ac為耕地自然固碳能力；Aci為農(nóng)作物i的碳吸收量；Cpi為農(nóng)作物i的含碳率；Yi為農(nóng)作物i的經(jīng)濟(jì)產(chǎn)量；Wi為農(nóng)作物i的水分系數(shù)；Hi為農(nóng)作物i的經(jīng)濟(jì)系數(shù)；Ri為農(nóng)作物i的根冠比；m為研究區(qū)主要作物類型數(shù)目（m=5）。各種農(nóng)作物固碳量的計(jì)算參數(shù)[27]詳見表2，遼寧省油料作物類型主要為花生，占比達(dá)到89.5%，所以油料的含碳率、水分系數(shù)、經(jīng)濟(jì)系數(shù)以及根冠比均以花生作為替代。

表2 農(nóng)作物固碳量計(jì)算參數(shù)Tab.2 Calculation parameters of crop carbon sequestration

參考相關(guān)研究中的農(nóng)作物種類多樣性計(jì)算方法[1,12]，結(jié)合遼寧省農(nóng)作物種植實(shí)際情況，挑選糧食作物、經(jīng)濟(jì)作物、蔬菜作物及瓜類作物進(jìn)行農(nóng)田系統(tǒng)生境質(zhì)量的測算。

式（3）中：Sim為農(nóng)田系統(tǒng)生境質(zhì)量系數(shù)；Pi為第i類農(nóng)作物播種面積占農(nóng)作物總播種面積的比重；i為農(nóng)作物類型；n為農(nóng)作物類型數(shù)。

3.3.2 評(píng)價(jià)模型

（1）耕地多功能指數(shù)計(jì)算。功能評(píng)價(jià)的計(jì)算公式：

式（4）中：Fi為不同維度耕地多功能指數(shù)；Wij為第i個(gè)評(píng)價(jià)指標(biāo)的權(quán)重；T′ij為i指標(biāo)中第j個(gè)歸一化后的數(shù)值。

（2）耕地價(jià)值核算。耕地綜合價(jià)值（V）是經(jīng)濟(jì)價(jià)值（Ve1）、社會(huì)價(jià)值（Vs）和生態(tài)價(jià)值（Ve2）三者的加和，計(jì)算公式如下：

耕地的經(jīng)濟(jì)價(jià)值指耕地利用過程中以農(nóng)產(chǎn)品收入減掉其生產(chǎn)周期中所有投入成本得到的凈收益[28]。計(jì)算公式如下：

式（6）中：Ve1為耕地經(jīng)濟(jì)價(jià)值；a為耕地的年均純收益（主要包括水稻、小麥、玉米、大豆、番薯）；r為耕地的收益還原率。耕地的收益還原率是指將耕地的純收益還原成耕地經(jīng)濟(jì)價(jià)值的比率，其本質(zhì)是資本投資的未來收益貼現(xiàn)到現(xiàn)在市場情況下的資本總價(jià)值，安全利率為1年期銀行存款利率，即1.5%。

耕地的社會(huì)價(jià)值是指耕地利用的直接性生產(chǎn)功能維持農(nóng)民生存生活、保障國家糧食安全以及社會(huì)穩(wěn)定等方面的外部價(jià)值[5,14,16]，主要由社會(huì)保障價(jià)值與社會(huì)穩(wěn)定價(jià)值構(gòu)成。計(jì)算公式如下：

式（7）—式（9）中：Vs為耕地社會(huì)價(jià)值；Va為單位面積耕地的社會(huì)保障價(jià)值；Vb為單位面積耕地的社會(huì)穩(wěn)定價(jià)值；Va1為單位面積耕地的生活保障價(jià)值；Va2為單位面積耕地的就業(yè)保障價(jià)值；Y1為政府提供的城鎮(zhèn)居民養(yǎng)老保險(xiǎn)金；Y2為政府提供的農(nóng)村居民養(yǎng)老保險(xiǎn)金；S為人均耕地面積；I1為農(nóng)民人均可支配收入；I2為城鎮(zhèn)居民人均可支配收入；U政府提供的城鎮(zhèn)居民失業(yè)保險(xiǎn)金；Vb為單位耕地社會(huì)安全價(jià)值；Qi為i種農(nóng)產(chǎn)品年總產(chǎn)量；Pi為i種農(nóng)產(chǎn)品的平均單價(jià)；A為區(qū)域耕地面積。

耕地生態(tài)價(jià)值是指耕地資源及所依附環(huán)境構(gòu)成的耕地生態(tài)系統(tǒng)承載的氣候調(diào)節(jié)、水土保持、凈化環(huán)境、優(yōu)化生境以及維持生物多樣性等多方面價(jià)值[29-30]，采用生態(tài)價(jià)值當(dāng)量系數(shù)法計(jì)算耕地的生態(tài)價(jià)值。計(jì)算公式如下:

式（10）中：Ve2為耕地生態(tài)價(jià)值；Ea為1單位生物當(dāng)量因子的價(jià)值量；F為單位面積當(dāng)量因子；r為主要農(nóng)作物種類；Pi為第i種農(nóng)作物的平均價(jià)格；qi為第i種農(nóng)作物面積單產(chǎn)；ai為第i種農(nóng)作物播種面積；A為總面積。

3.3.3 敏感度分析

為分析耕地功能對(duì)價(jià)值的影響作用，借鑒敏感度分析模型[31]，計(jì)算公式如下：

式（11）中：δf為遼寧省14地市耕地功能對(duì)其價(jià)值變化的敏感程度，δp、δl、δe為耕地生產(chǎn)功能、生活功能和生態(tài)功能的敏感系數(shù)，分別用Ve1-Fp、Vs-Fl、Ve2-Fe表示；Vt和Vt+1分別為研究基期和末期的耕地價(jià)值；Ft和Ft+1分別為該時(shí)期內(nèi)14地市基期和末期的耕地多功能指數(shù)。若δ＜0，表明耕地價(jià)值與耕地功能呈反向變化，耕地價(jià)值與功能不存在響應(yīng)關(guān)系；若δ＞0，表明耕地價(jià)值與耕地功能呈正向變化，耕地價(jià)值與功能存在響應(yīng)關(guān)系，且δ值越大響應(yīng)程度越高，即耕地功能的較小變化造成其價(jià)值的較大波動(dòng)，具體劃分見表3。

表3 敏感度劃分標(biāo)準(zhǔn)Tab.3 Sensitivity classification standard

4 結(jié)果分析

4.1 耕地多功能演變特征

從時(shí)序上來看，2004—2019年，遼寧省耕地生產(chǎn)功能指數(shù)（Fp）有所上升（圖3）。朝陽、本溪、撫順增幅較大，其功能指數(shù)分別從0.070 6增至0.171 9、0.027 8增至0.059 3、0.040 6增至0.085 8；大連、丹東、營口、阜新、遼陽等增幅均超過50%；沈陽增幅最小，僅增加19.61%?？臻g層面上呈現(xiàn)北高南低趨勢(shì)，高值區(qū)與較高值區(qū)集中在遼北區(qū)域，表現(xiàn)出空間擴(kuò)展態(tài)勢(shì)；中值區(qū)長期聚集在遼中區(qū)域，呈現(xiàn)“俱樂部收斂”特征；低值區(qū)與較低值區(qū)分列遼東、遼西兩側(cè)（圖4）。主要得益于遼寧省各地區(qū)耕地復(fù)種指數(shù)、糧食單產(chǎn)及農(nóng)業(yè)產(chǎn)業(yè)效益得到不同程度提高，而西部低山丘陵區(qū)的耕地面積雖有所增加，但該區(qū)域的生態(tài)脆弱性較高且耕地質(zhì)量較差。此外，沈陽、大連、盤錦、阜新等地存在耕地收縮現(xiàn)象。

2004—2019年，遼寧省耕地生活功能指數(shù)（Fl）有所提升（圖3），但增幅較小。朝陽、大連增幅均超過50%，其功能指數(shù)分別從0.192 5增至0.325 3、0.192 5增至0.295 9；沈陽、本溪、丹東、錦州、盤錦、葫蘆島增幅趨同，處在20%～30%之間；阜新、撫順有式微下降，降幅為3.6%和6.8%?？臻g層面上呈現(xiàn)由“南北高、東西低”轉(zhuǎn)變?yōu)椤澳媳蔽鞲?、東低”格局特征，高值區(qū)自遼北逐步擴(kuò)散外圍區(qū)域，最終散列南、北、西三側(cè)；低值區(qū)與較低值區(qū)始終集中在遼東區(qū)域。其中，2009—2014年的耕地生活功能衰退明顯，除朝陽和葫蘆島之外，其他12地市均出現(xiàn)不同程度降低（圖4）。主要由于遼寧省東部、西部多為丘陵，耕地破碎化程度高且年平均降水量少，中部下遼河平原的農(nóng)業(yè)生產(chǎn)自然條件與耕地利用的匹配程度相對(duì)較高，永久基本農(nóng)田集中分布于此，而且各地農(nóng)村社會(huì)保障水平、農(nóng)業(yè)機(jī)械化程度不斷提高，部分地區(qū)農(nóng)業(yè)勞動(dòng)力梯度析出等。

圖3 2004—2019年遼寧省耕地多功能指數(shù)均值與變異系數(shù)Fig.3 Average value and variation coefficient of cultivated land multifunction index in Liaoning Province from 2004 to 2019

圖4 2004—2019年遼寧省耕地多功能演變格局Fig.4 Evolution pattern of cultivated land multifunction in Liaoning Province from 2004 to 2019

2004—2019年，遼寧省耕地生態(tài)功能指數(shù)（Fe）顯著上升（圖3）。營口、大連、遼陽、朝陽增幅較大，其功能指數(shù)分別從0.040 7增至0.119 1、0.043 8增至0.112 9、0.055 2增至0.131 2、0.080 3增至0.164 1；沈陽、鞍山、撫順、葫蘆島增幅處于60%～70%之間；鐵嶺增幅最小且為52.9%?？臻g層面上呈現(xiàn)“西北高、東南低”趨勢(shì)，高值區(qū)與較高值區(qū)聚集遼北、遼西區(qū)域，具有“擴(kuò)散—收斂—再擴(kuò)散”態(tài)勢(shì)；低值區(qū)與較低值區(qū)轉(zhuǎn)移至遼中南、遼東區(qū)域。其中，2009—2014年，大連、鞍山、本溪、錦州、營口、葫蘆島的耕地生態(tài)功能出現(xiàn)衰退，但降幅未超過20%。此外，從耕地多功能指數(shù)的變異系數(shù)來看，各功能的變異系數(shù)均處于0～0.5區(qū)間，屬于中等變異，并且總體呈現(xiàn)下降趨勢(shì)，表明遼寧省耕地多功能間的空間差異在逐年縮小；而且耕地生產(chǎn)功能與生態(tài)功能指數(shù)的平均值出現(xiàn)交替上升現(xiàn)象，說明耕地多功能在協(xié)同增益（圖4）。主要是近年來遼寧省各地區(qū)的化肥農(nóng)藥減施增效技術(shù)覆蓋面日益擴(kuò)展，以及高標(biāo)準(zhǔn)農(nóng)田建設(shè)、全域土地綜合整治等項(xiàng)目實(shí)施，不斷推進(jìn)“山水林田湖草沙”一體化保護(hù)和修復(fù)，為農(nóng)田生物多樣性保護(hù)創(chuàng)造了有利條件，例如丹東等地實(shí)施了新一輪退耕還林措施。

4.2 耕地價(jià)值演變特征

2004—2019年，遼寧省耕地價(jià)值呈總體上升趨勢(shì)（圖5），耕地社會(huì)價(jià)值占比最多且變化幅度最大，耕地經(jīng)濟(jì)價(jià)值次之，耕地生態(tài)價(jià)值最小。具體來看，地均經(jīng)濟(jì)價(jià)值方面，各節(jié)點(diǎn)年份平均值分別為5.81萬元/hm2、8.58萬元/hm2、10.76萬元/hm2、11.44萬元/hm2，2004年高于平均值的地市為沈陽、大連、鞍山、錦州、營口和盤錦；其中大連最高，阜新最低（圖6）。與2004年相比，2009年整體上升47.9%，年變化率為9.59%；2014年高于平均值的地市為沈陽、大連、鞍山、撫順、錦州、遼陽和朝陽；2019年地均經(jīng)濟(jì)價(jià)值最大是沈陽，最小為阜新，這與水田和旱地的空間分布有關(guān)。2004—2019年，葫蘆島地均經(jīng)濟(jì)價(jià)值的增幅明顯，由5.21萬元/hm2增至11.27萬元/hm2。主要是研究期內(nèi)葫蘆島耕地面積有所增加，糧食生產(chǎn)支持政策實(shí)施與農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè)等促進(jìn)了農(nóng)業(yè)產(chǎn)值與糧食產(chǎn)量的增長。

圖5 2004—2019年遼寧省耕地價(jià)值均值與變異系數(shù)Fig.5 Average value and variation coefficient of cultivated land value in Liaoning Province from 2004 to 2019

地均社會(huì)價(jià)值方面，各節(jié)點(diǎn)年份平均值分別為105.81萬元/hm2、156.58萬元/hm2、139.19萬元/hm2、144.49萬元/hm2，2004年高值集中分布在大連、營口、鞍山和丹東等遼南地區(qū)，低值主要分布在朝陽、阜新和錦州等遼西地區(qū)，與朝陽的建平、朝陽縣以及阜新的彰武、蒙古族自治縣等耕地轉(zhuǎn)為建設(shè)用地有關(guān)；大連最高，阜新最小。2004—2009年上升幅度為47.9%，高低值鑲嵌分布于遼南與遼西地區(qū)；2014年上升幅度略有下降，變化率為17.39%；2019年最大值為營口，最小值為阜新。2004—2019年間地均社會(huì)價(jià)值先增后降再增加，呈“N”型變化趨勢(shì)（圖5—圖6），主要受不同時(shí)期城鎮(zhèn)空間擴(kuò)張占用耕地的影響，即與各地區(qū)的城鎮(zhèn)化發(fā)展進(jìn)程一致。其中，變化幅度最大的是朝陽，由65.62萬元/hm2增至129.26萬元/hm2，與朝陽部分地區(qū)的林地與水域轉(zhuǎn)為耕地有關(guān)，而且該地區(qū)建設(shè)用地規(guī)模及擴(kuò)張速度也較低；鐵嶺增幅最小，由79.59萬元/hm2增至86.88萬元/hm2。

圖6 2004—2019年遼寧省耕地價(jià)值變化態(tài)勢(shì)Fig.6 Evolution pattern of cultivated land value in Liaoning Province from 2004 to 2019

地均生態(tài)價(jià)值方面，各節(jié)點(diǎn)年份平均值分別為4.58萬元/hm2、6.79萬元/hm2、9.60萬元/hm2、11.76萬元/hm2，2004年高于平均值的地市為沈陽、錦州、遼陽、盤錦和鐵嶺，2009年、2014年和2019年均集中在沈陽、錦州、營口、阜新、遼陽、盤錦和鐵嶺。與2004年相比，2009年除朝陽降低0.42萬元/hm2外，其他地市均呈現(xiàn)增長態(tài)勢(shì)（圖6）；2014年高值分布在遼中地區(qū)，低值主要在遼東地區(qū)；2019年最大值為沈陽，最小值為本溪，這與遼中地區(qū)耕地地勢(shì)平緩、水田集中有關(guān)。2004—2019年，地均生態(tài)價(jià)值每年漲幅10%，變化幅度最大的是朝陽，由1.65萬元/hm2增至8.01萬元/hm2，該地區(qū)具有一定的耕地后備資源，大連變化幅度最小且年均變化率為14.4%。此外，從耕地價(jià)值平均值的變異系數(shù)來看（圖5），各類價(jià)值的變異系數(shù)均處于0～1之間的中等變異層級(jí)，耕地生態(tài)價(jià)值的變異系數(shù)始終高于0.5，但總體上呈逐年下降趨勢(shì)，表明遼寧省耕地價(jià)值間的內(nèi)部差異在縮?。欢页鐣?huì)價(jià)值有所波動(dòng)之外，其他兩類價(jià)值的平均值緩慢上升，說明耕地復(fù)合價(jià)值在協(xié)同增益（圖6）。主要是近15年遼寧省各地市的耕地利用動(dòng)態(tài)變化不一，下遼河平原及經(jīng)濟(jì)發(fā)達(dá)地區(qū)的耕地利用變化相對(duì)劇烈，引發(fā)耕地的自然質(zhì)量條件、集中連片程度、基礎(chǔ)設(shè)施等立地環(huán)境條件產(chǎn)生適應(yīng)性變化，同時(shí)也表明氣候變化趨勢(shì)、經(jīng)濟(jì)社會(huì)轉(zhuǎn)型升級(jí)、土地利用結(jié)構(gòu)轉(zhuǎn)變以及農(nóng)業(yè)人口轉(zhuǎn)移等將深刻影響耕地價(jià)值演變?？傊捎诟貎r(jià)值概念內(nèi)涵缺乏共識(shí)、價(jià)值體系及認(rèn)知途徑不夠完整，耕地價(jià)值體系殘缺或?qū)o形中壓縮農(nóng)用地轉(zhuǎn)為建設(shè)用地的綜合成本，抑制耕地保護(hù)中政策激勵(lì)—規(guī)制約束背景下的規(guī)劃引導(dǎo)和市場機(jī)制作用，進(jìn)而導(dǎo)致大量優(yōu)質(zhì)耕地被占用。

4.3 耕地功能與價(jià)值的響應(yīng)分析

從上述分析可見，遼寧省耕地價(jià)值與耕地功能的演變趨勢(shì)接近，一定程度說明耕地功能演變?nèi)允侵萍s其價(jià)值變化的客觀基礎(chǔ)，即二者之間存在響應(yīng)關(guān)系。從圖7可以看出，2004—2019年，耕地經(jīng)濟(jì)價(jià)值對(duì)生產(chǎn)功能的響應(yīng)程度較高，全省14地市均為敏感區(qū)域，丹東、鐵嶺、阜新、葫蘆島為高度敏感區(qū)域，中度及以上的敏感區(qū)域占比64.29%（圖8）。其中，2004—2009年，高度敏感區(qū)域較多且集中在沈陽、錦州、撫順、遼陽、丹東；輕度、中度及較高敏感區(qū)域集中在遼中和遼東地區(qū)；無敏感區(qū)域分布在葫蘆島、朝陽、阜新、鐵嶺。2009—2014年，全省有7個(gè)低度敏感區(qū)域與6個(gè)較低敏感區(qū)域，僅葫蘆島為無敏感區(qū)域，該時(shí)期未出現(xiàn)中度敏感及以上的區(qū)域。2014—2019年，低度敏感區(qū)域較多，占比為57.14%；無敏感區(qū)域分布在大連、盤錦、遼陽、本溪、撫順；較低敏感及以上的區(qū)域僅分布在鐵嶺。近年來遼寧省農(nóng)業(yè)科技推廣促進(jìn)農(nóng)業(yè)現(xiàn)代化工程建設(shè)，土地流轉(zhuǎn)推動(dòng)耕地規(guī)?；图s化利用，糧食播種面積、糧食產(chǎn)量和農(nóng)業(yè)產(chǎn)值穩(wěn)定提升，從而使得耕地經(jīng)濟(jì)價(jià)值對(duì)耕地生產(chǎn)功能保持積極的響應(yīng)關(guān)系。

2004—2019年，耕地社會(huì)價(jià)值對(duì)生活功能的響應(yīng)程度較低，全省有9個(gè)地市為低度敏感區(qū)域，大連、阜新、撫順為無敏感區(qū)域，錦州為較低敏感區(qū)域，僅遼陽是較高敏感區(qū)域（圖7）。其中，2004—2009年，低度敏感區(qū)域占比達(dá)78.57%（圖8）；較低敏感區(qū)域分布在錦州、撫順，僅阜新為高度敏感區(qū)域。2009—2014年，低度敏感區(qū)域占比為71.43%；葫蘆島、錦州、遼陽均為無敏感區(qū)域，高度敏感區(qū)域僅分布在朝陽。2014—2019年，無敏感區(qū)域占比為64.29%，大連、沈陽、鐵嶺、撫順、遼陽均為無敏感區(qū)域，該階段未出現(xiàn)較低敏感及以上的區(qū)域。近年來遼寧省耕地要素投入水平增強(qiáng)、農(nóng)業(yè)機(jī)械化水平提升，對(duì)農(nóng)業(yè)勞動(dòng)力產(chǎn)生一定的替代效應(yīng)，單位面積投入帶來的綜合收益具有顯著的擴(kuò)散效應(yīng)，導(dǎo)致耕地社會(huì)價(jià)值對(duì)生活功能形成消極的響應(yīng)關(guān)系。

2004—2019年，耕地生態(tài)價(jià)值對(duì)生態(tài)功能的響應(yīng)程度一般，全省有10個(gè)地市為低度敏感區(qū)域，朝陽、阜新、撫順為較低敏感區(qū)域，高度敏感區(qū)域僅分布在本溪（圖7）。其中，2004—2009年，低度敏感區(qū)域占比高達(dá)85.71%（圖8）；本溪為高度敏感區(qū)域，朝陽為無敏感區(qū)域。2009—2014年，中度敏感及以上的區(qū)域占比為50%；朝陽、盤錦、丹東為高度敏感區(qū)域，沈陽、阜新、撫順為較高敏感區(qū)域；大連、葫蘆島為低度敏感區(qū)域，僅遼陽為中度敏感區(qū)域。2014—2019年，低度敏感區(qū)域仍高居85.71%；葫蘆島為較低敏感區(qū)域，盤錦為無敏感區(qū)域。近年來遼寧省水田減少、旱地增加，加劇了部分地區(qū)“水土氣熱”條件的不匹配，農(nóng)業(yè)勞動(dòng)力稟賦關(guān)聯(lián)約束和化肥農(nóng)藥不合理利用，對(duì)區(qū)域農(nóng)田生態(tài)系統(tǒng)造成負(fù)面效應(yīng)，進(jìn)而使得耕地生態(tài)價(jià)值對(duì)生態(tài)功能保持一定的響應(yīng)關(guān)系。

圖7 2004—2019年遼寧省耕地功能與價(jià)值的時(shí)空響應(yīng)Fig.7 Spatial-temporal response of cultivated land multifunction and its value in Liaoning Province from 2004 to 2019

圖8 2004—2019年遼寧省耕地功能與價(jià)值響應(yīng)的敏感區(qū)域占比Fig.8 Proportion of sensitive areas of cultivated land function and value response in Liaoning Province from 2004 to 2019

綜合來看，耕地功能演變與價(jià)值響應(yīng)的敏感區(qū)域數(shù)量、空間集聚特征及其敏感程度總體呈現(xiàn)動(dòng)態(tài)變化趨勢(shì)，不同時(shí)期的各類價(jià)值對(duì)功能的響應(yīng)存在明顯分化現(xiàn)象，因此，實(shí)施區(qū)域耕地多功能管理是保障耕地復(fù)合價(jià)值顯化的關(guān)鍵舉措?？傊?，在新型城鎮(zhèn)化、生態(tài)文明建設(shè)、東北全面振興背景下，遼寧省耕地?cái)?shù)量、質(zhì)量、生態(tài)、安全狀況及其利用方式與空間結(jié)構(gòu)均發(fā)生變化，耕地利用強(qiáng)度、規(guī)模化與集約化程度、耕地綜合產(chǎn)能等不斷提高，促進(jìn)了耕地多功能提升及復(fù)合價(jià)值增益，并引致耕地功能與價(jià)值之間產(chǎn)生內(nèi)在響應(yīng)。但不可忽視經(jīng)濟(jì)比較效益與機(jī)會(huì)成本驅(qū)動(dòng)下的林地、草地、水域向耕地轉(zhuǎn)移之后再轉(zhuǎn)為城鎮(zhèn)建設(shè)或工業(yè)用地等，逐漸形成建設(shè)用地蔓延擴(kuò)張的慣性路徑依賴，并成為阻抑耕地功能及其價(jià)值協(xié)同發(fā)展的重要因素。可見，耕地功能與價(jià)值之間存在著特定關(guān)聯(lián)，但并非是同增同減的線性響應(yīng)，而是形成了包含正向與負(fù)向在內(nèi)的耦合互饋過程，即在時(shí)間序列上表現(xiàn)為耦合效應(yīng)的非線性動(dòng)態(tài)變化，在空間層面上呈現(xiàn)出周期波動(dòng)的非均衡特征。

5 結(jié)論與討論

5.1 結(jié)論

本文以遼寧省作為典型研究區(qū)，基于2004—2019年全省14個(gè)地市的相關(guān)統(tǒng)計(jì)數(shù)據(jù)，構(gòu)建耕地多功能評(píng)價(jià)及其價(jià)值核算的指標(biāo)體系，揭示了耕地生產(chǎn)功能、生活功能、生態(tài)功能的演變特征，測算了遼寧省各地市耕地價(jià)值，闡明了耕地多功能演變及“經(jīng)濟(jì)—社會(huì)—生態(tài)”價(jià)值響應(yīng)關(guān)系，主要結(jié)論如下：

（1）2004—2019年，全省耕地多功能演變具有顯著的空間分異特征，呈“北高南低、東西一致”態(tài)勢(shì)，其中耕地生態(tài)功能提升最大，耕地生產(chǎn)功能次之，耕地生活功能增幅最小?？梢姡芯繀^(qū)化肥農(nóng)藥零增長、減量化專項(xiàng)活動(dòng)的實(shí)施成效對(duì)耕地生態(tài)功能增長具有較強(qiáng)推動(dòng)作用。

（2）2004—2019年，各地市耕地價(jià)值呈整體上升趨勢(shì)，綜合價(jià)值的平均值為153.846 8萬元/hm2。其中，耕地社會(huì)價(jià)值最大，平均值為136.514 8萬元/hm2，總體呈“N”型變化態(tài)勢(shì)；耕地經(jīng)濟(jì)價(jià)值次之，平均值為9.149 1萬元/hm2；耕地生態(tài)價(jià)值最小，平均值為8.182 9萬元/hm2，但變化幅度較大。以社會(huì)價(jià)值為主導(dǎo)的耕地價(jià)值體系構(gòu)成凸顯出耕地資源的農(nóng)民生計(jì)保障和社會(huì)穩(wěn)定作用，也側(cè)面反映了惠農(nóng)支農(nóng)政策的穩(wěn)定性和連續(xù)性。

（3）耕地功能與價(jià)值存在一定的響應(yīng)關(guān)系，總體表現(xiàn)出伴隨時(shí)間發(fā)展的階段性變化與空間分異的非均衡性特征。其中，研究區(qū)耕地生產(chǎn)功能與耕地經(jīng)濟(jì)價(jià)值具有積極的響應(yīng)關(guān)系，耕地生活功能與耕地社會(huì)價(jià)值具有消極的響應(yīng)關(guān)系，耕地生態(tài)功能與耕地生態(tài)價(jià)值的響應(yīng)程度處于一般水平。由此說明基于耕地多功能供需差異的統(tǒng)籌調(diào)控將深刻影響耕地復(fù)合價(jià)值實(shí)現(xiàn)及其協(xié)同發(fā)展，各地市應(yīng)進(jìn)一步加強(qiáng)耕地保護(hù)與質(zhì)量建設(shè)以顯化其功能及價(jià)值。

5.2 討論

在當(dāng)前已有研究的基礎(chǔ)上進(jìn)一步探析耕地利用系統(tǒng)內(nèi)部要素聚合、結(jié)構(gòu)均衡、功能協(xié)調(diào)、價(jià)值穩(wěn)定的過程規(guī)律及主控因素，對(duì)區(qū)域耕地資源優(yōu)化調(diào)控具有積極意義。本文基于耕地功能與價(jià)值關(guān)聯(lián)的理論分析，從二者“解耦”視角初步分析了耕地功能與價(jià)值的響應(yīng)關(guān)系，對(duì)豐富和發(fā)展耕地資源功能—價(jià)值互饋研究或具有一定實(shí)踐價(jià)值。同時(shí)，耕地承載著重要的社會(huì)價(jià)值，但是由于耕地價(jià)值內(nèi)涵認(rèn)知不足、價(jià)值體系缺位、非市場價(jià)值與耕保“外部性”孱弱等因素，在比較利益驅(qū)動(dòng)下，容易弱化農(nóng)民種糧成本、惠農(nóng)補(bǔ)貼、綜合收益與相關(guān)主體保護(hù)行為的協(xié)同性，如何有效遏制耕地“非農(nóng)化”、防止耕地“非糧化”問題成為當(dāng)務(wù)之急。因此，在耕地功能與價(jià)值研究的基礎(chǔ)上，探索信息、物質(zhì)、能量耦合互饋下耕地功能—價(jià)值之間的復(fù)雜關(guān)系及響應(yīng)機(jī)理，深入揭示耕地多功能影響價(jià)值的關(guān)鍵過程與驅(qū)動(dòng)機(jī)理，闡釋耕地多功能形態(tài)因素對(duì)區(qū)域耕地價(jià)值影響以及與經(jīng)濟(jì)社會(huì)高質(zhì)量發(fā)展之間的響應(yīng)適應(yīng)研究應(yīng)成為近期關(guān)注的方向。例如，統(tǒng)籌考慮不同外部條件下耕地利用過程的投入/產(chǎn)出、結(jié)構(gòu)形態(tài)、功能強(qiáng)弱變化及價(jià)值演替等因素，強(qiáng)化耕地利用系統(tǒng)要素共享、結(jié)構(gòu)聯(lián)動(dòng)、功能融合與價(jià)值增益，推進(jìn)高標(biāo)準(zhǔn)農(nóng)田建設(shè)、優(yōu)質(zhì)糧食工程、大豆振興計(jì)劃、綠色種養(yǎng)循環(huán)農(nóng)業(yè)試點(diǎn)、秸稈綜合利用、畜禽糞污資源化利用等高效實(shí)施，有助于落實(shí)耕地保護(hù)政策、增強(qiáng)農(nóng)民的可持續(xù)增收能力。此外，考慮到耕地系統(tǒng)內(nèi)部要素、結(jié)構(gòu)、功能與價(jià)值關(guān)聯(lián)具有高度復(fù)雜性和多維性，不同功能與價(jià)值的互作強(qiáng)度、方向和速率存在差異。本文主要基于定量研究視角下的關(guān)聯(lián)規(guī)則分析了耕地功能與價(jià)值之間的關(guān)系，后續(xù)研究將進(jìn)一步使用回歸分析等方法等揭示其驅(qū)動(dòng)因素及互饋機(jī)制。