黃耀文，張曉平

（東華理工大學(xué)測(cè)繪工程學(xué)院，江西·南昌 330013）

國家重點(diǎn)生態(tài)功能區(qū)指對(duì)國家，區(qū)域等存在生態(tài)安全意義的限制開發(fā)區(qū)，即在水源和水土保護(hù)，荒漠化防治，生物多樣性維護(hù)等四個(gè)方面具有關(guān)鍵作用的區(qū)域限制開發(fā)。自《國家重點(diǎn)生態(tài)功能區(qū)轉(zhuǎn)移支付辦法》等相關(guān)文件實(shí)施之后，轉(zhuǎn)移支付成為我國重點(diǎn)生態(tài)功能區(qū)補(bǔ)償?shù)某Ｓ梅绞胶捅匾侄?。生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估是重點(diǎn)生態(tài)功能區(qū)轉(zhuǎn)移支付績(jī)效評(píng)價(jià)和激勵(lì)約束的重要理論依據(jù)，是確立生態(tài)補(bǔ)償標(biāo)準(zhǔn)的價(jià)值基礎(chǔ)根本[1]。

迄今，對(duì)于生態(tài)系統(tǒng)服務(wù)價(jià)值的研究眾多。從研究?jī)?nèi)容上看，有浙江省土地利用景觀格局演變對(duì)生態(tài)系統(tǒng)服務(wù)的影響過程及機(jī)制研究[2]，有多山地丘陵的生態(tài)脆弱敏感區(qū)域土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響研究[3]；有結(jié)合當(dāng)量因子法、空間自相關(guān)分析、熱點(diǎn)分析等方法研究洞庭湖生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化，直接使用地理探測(cè)器工具探究生態(tài)系統(tǒng)服務(wù)價(jià)值空間分異的驅(qū)動(dòng)因素[4]；有分析長(zhǎng)江經(jīng)濟(jì)帶周邊省份近年來的生態(tài)系統(tǒng)服務(wù)價(jià)值特征并采用灰色預(yù)測(cè)模型對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行預(yù)測(cè)[5]；有利用CA-Markov模型預(yù)測(cè)未來不同情景下白洋淀流域土地利用變化情況，估算各情景下生態(tài)系統(tǒng)服務(wù)價(jià)值，并對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行時(shí)空分析[6]。從研究方法上看，主要有兩種方法，一是功能價(jià)值法，如運(yùn)用市場(chǎng)價(jià)值法、造林成本法、影子工程法、費(fèi)用替代法等方法估算甘肅祁連山國家自然保護(hù)區(qū)山丹馬場(chǎng)草地生態(tài)系統(tǒng)的服務(wù)價(jià)值[7]；二是當(dāng)量因子法，國內(nèi)學(xué)者謝高地等[8]結(jié)合我國實(shí)際情況在Costanza等[9]研究的基礎(chǔ)上進(jìn)行修改，提出適合我國生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值表，并在后續(xù)進(jìn)行了改進(jìn)[10]。該改進(jìn)成果在土地利用變化與生態(tài)系統(tǒng)服務(wù)價(jià)值的研究中廣泛應(yīng)用。從研究尺度上來看，涵蓋國家級(jí)[11]、省級(jí)[12]、市級(jí)[13]、城市群[14]、流域[15-16]、區(qū)域[17-18]等空間尺度。近年，有關(guān)研究拓展到微觀尺度，以格網(wǎng)為研究單元展開分析[19-20]。格網(wǎng)尺度的研究打破了行政區(qū)的研究范式[21]，可以更有效地獲取土地利用信息，更準(zhǔn)確地判斷生態(tài)系統(tǒng)服務(wù)價(jià)值空間變化關(guān)系。然而在格網(wǎng)尺度上，關(guān)于重點(diǎn)生態(tài)功能區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值量化分析較少。

江西省大余縣屬于南嶺山地森林及生物多樣性生態(tài)功能區(qū)，是重要的水源涵養(yǎng)區(qū)。隨著工業(yè)化城鎮(zhèn)化的推進(jìn)，生態(tài)系統(tǒng)功能退化。因此，研究大余縣生態(tài)系統(tǒng)服務(wù)價(jià)值對(duì)生態(tài)保護(hù)具有重要意義。研究以大余縣為研究區(qū)，依據(jù)多時(shí)段土地利用覆被數(shù)據(jù)，修正生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表，在格網(wǎng)尺度分析區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值空間變化特征，以期為國家重點(diǎn)生態(tài)功能區(qū)的生態(tài)政策制定提供參考。

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

大余縣地處南嶺山地森林及生物多樣性生態(tài)功能區(qū)，主體功能定位是贛、粵、湘、桂4省重要的水源涵養(yǎng)區(qū)，生態(tài)功能區(qū)位十分重要。

全縣位于江西省的西南角，與湖南、廣東交界，是江西前往廣東的重要通道之一，享有江西南邊門戶的稱謂。全境三面環(huán)山，東部平坦的西高東低地勢(shì)，呈現(xiàn)為東西長(zhǎng)，南北寬的長(zhǎng)條狀（圖1）。全縣共轄11個(gè)鄉(xiāng)鎮(zhèn)，其中山地、丘陵約占81.62%，平原約占18.38%。大余縣四季分明，溫和濕潤，屬于亞熱帶季風(fēng)氣候。

圖1 研究區(qū)地貌形態(tài)Fig.1 Geomorphic morphology of the study area

研究中涉及的各期土地覆蓋數(shù)據(jù)（30 m年度土地覆蓋數(shù)據(jù)集，總體準(zhǔn)確率達(dá)80%）來源于武漢大學(xué)遙感信息處理研究所（http://irsip.whu.edu.cn/）。為便于生態(tài)系統(tǒng)服務(wù)價(jià)值分析，將大余縣土地利用類型劃分為草地、耕地、建設(shè)用地、林地、水域和未利用地（圖2）。水田旱地比例來源于大余縣自然資源局土地變更調(diào)查數(shù)據(jù)。糧食作物播種面積、產(chǎn)量糧食作物和總產(chǎn)值來源于2020年《大余年鑒》，2020年《贛州統(tǒng)計(jì)年鑒》和國家統(tǒng)計(jì)局（https://data.stats.gov.cn/）。

圖2 大余縣各期土地利用現(xiàn)狀圖Fig.2 Current land use map for different periods in Dayu County

2 研究方法

2.1 生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估方法

以謝高地等[5]改進(jìn)后當(dāng)量表為基礎(chǔ)，并根據(jù)徐麗芬等[22]的地區(qū)修正公式（1）進(jìn)行研究區(qū)地域修正。

式中：θ為研究區(qū)修正因子，Q和Q0分別為大余縣和全國單位面積平均產(chǎn)量，Ei表示i種地類修正后當(dāng)量因子，Ei0表示謝高地等確定的同種地類當(dāng)量因子。

2019年大余縣耕地單位面積糧食產(chǎn)量為5.606噸/公頃；與此同時(shí)，全國耕地單位面積糧食產(chǎn)量為5.720噸/公頃；根據(jù)公式(1)得到大余縣生態(tài)服務(wù)當(dāng)量修正因子為0.98。根據(jù)贛州統(tǒng)計(jì)年鑒，2019年的糧食作物總產(chǎn)值為23786萬元，糧食總產(chǎn)量89613噸，估算大余縣2019年糧食價(jià)格為2654.30元/噸，現(xiàn)有單位面積農(nóng)田提供的食物生產(chǎn)服務(wù)經(jīng)濟(jì)價(jià)值是在沒有人為干預(yù)情況下的自然生態(tài)系統(tǒng)提供的經(jīng)濟(jì)價(jià)值的7倍[8]，由此估算大余縣1個(gè)生態(tài)服務(wù)價(jià)值量因子的經(jīng)濟(jì)價(jià)值量為2125.72元。

通過地區(qū)修正系數(shù)修訂得到大余縣各地類單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)（表1）。其中，耕地中水田旱地的比例取自2013至2018年土地變更調(diào)查數(shù)據(jù)中水田旱地的平均比例，為0.905：0.095；水域采用朱增云等[23]的研究：取水系和濕地的均值；林地、草地、未利用地和建設(shè)用地借鑒冉玉菊等[14]的處理方法。

表1 大余縣各地類單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)Table 1 Coefficient of ecosystem services value per unit area over Dayu County

生態(tài)系統(tǒng)服務(wù)價(jià)值估算公式如下：

式中：ESV為生態(tài)系統(tǒng)服務(wù)價(jià)值，Ak是地類k的面積，VCk是生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)，ESVf為第f項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值，VCfk是第k地類的第f項(xiàng)服務(wù)價(jià)值系數(shù)。

2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值空間聚集特征分析

空間自相關(guān)是評(píng)估空間分異規(guī)律的常用方法。其中莫蘭指數(shù)是研究的一項(xiàng)常用指標(biāo)，分為全局Moran′s I和局部Moran′s I。全局Moran′s I用于檢驗(yàn)大余縣ESV是否存在空間相關(guān)性以及相關(guān)程度有多大，取值范圍在[-1,1]，當(dāng)Moran′s I＞ 0時(shí)，呈正相關(guān)；當(dāng)Moran′s I＜ 0時(shí)，為負(fù)相關(guān)；當(dāng)Moran′s I= 0時(shí)，空間分布呈現(xiàn)隨機(jī)分布狀態(tài)。局部Moran′s I主要通過GeoDa空間數(shù)據(jù)分析軟件實(shí)現(xiàn)，局部Moran′s I可以進(jìn)一步確定大余縣的ESV空間集聚區(qū)或孤立區(qū)所在的位置。其中，LISA聚類圖能清晰反映地域?qū)傩缘年P(guān)聯(lián)性[24]。

3 結(jié)果分析

3.1 大余縣土地利用變化

大余縣主要土地利用類型為林地，耕地和建設(shè)用地。研究期間，林地、水域和未利用地的面積均減少，其中林地面積減少最多，為2289.48公頃；建設(shè)用地、耕地和草地面積均增加，建設(shè)用地面積增加最多，達(dá)1536.02公頃（表2）。

表2 2000-2020年大余縣土地利用變化Table 2 Land use change in Dayu County from 2000 to 2020

通過計(jì)算大余縣土地利用轉(zhuǎn)移矩陣[25]，分析各個(gè)地類的轉(zhuǎn)移方向和數(shù)量。2000-2005年，耕地主要轉(zhuǎn)向林地和建設(shè)用地，分別為809.26公頃（占總轉(zhuǎn)出量的72.62%）和190.17公頃（占總轉(zhuǎn)出量的17.07%）；林地轉(zhuǎn)向耕地1180.82公頃，占總轉(zhuǎn)出量的99.00%；草地、建設(shè)用地，水域和未利用地轉(zhuǎn)向其他地類轉(zhuǎn)換的面積較小；2005-2010年、2010-2015年及2015-2020年三期的轉(zhuǎn)換關(guān)系與第一期類似，均以耕地、林地、建設(shè)用地之間轉(zhuǎn)換為主（表3）。總體上，大余縣各地類的總轉(zhuǎn)出量為逐期遞增，符合社會(huì)不斷發(fā)展，人類活動(dòng)更加頻繁的規(guī)律；建設(shè)用地、林地、以及耕地三者之間的轉(zhuǎn)換是由于建設(shè)用地?cái)U(kuò)張占用耕地，林地開發(fā)用于耕地占補(bǔ)平衡；2005-2010年的耕地和水域的轉(zhuǎn)換關(guān)系與2015-2020年的關(guān)系剛好相反，前期由于耕地經(jīng)濟(jì)效益小于漁業(yè)效益，利益驅(qū)使導(dǎo)致耕地轉(zhuǎn)向坑塘，或因水利工程等基礎(chǔ)設(shè)施建設(shè)占用耕地，而后期水域向耕地轉(zhuǎn)移的原因可能是河湖范圍收縮導(dǎo)致部分河湖界內(nèi)的灘涂等地塊被開墾成耕地。

表3 2000-2020年大余縣土地利用轉(zhuǎn)移矩陣 Table 3 Land use transfer matrix in Dayu County from 2000 to 2020（hm2）

3.2 大余縣生態(tài)系統(tǒng)服務(wù)價(jià)值變化

根據(jù)敏感性指數(shù)方法[15,26]，對(duì)各個(gè)地類的生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)上下調(diào)整50%，計(jì)算大余縣2000-2020年各個(gè)地類ESV的敏感性指數(shù)。研究期內(nèi)各個(gè)地類的敏感性系數(shù)均小于1，意味著大余縣的總ESV對(duì)于生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)缺乏彈性，表明研究采用的生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)適合研究區(qū)的情況，生態(tài)服務(wù)價(jià)值結(jié)果可信。

（1）不同土地類型的ESV變化

根據(jù)大余縣當(dāng)量表和各個(gè)時(shí)段各地類的面積，計(jì)算大余縣的生態(tài)系統(tǒng)服務(wù)價(jià)值。大余縣2000年、2005年、2010年、2015年和2020年的生態(tài)系統(tǒng)服務(wù)價(jià)值分別是47.92億元、47.72億元、47.62億元、47.35億元和46.43億元，總體上呈逐年遞減趨勢(shì)（表4）。林地的ESV最高，占總ESV超過94%，且林地ESV減少最快，是總ESV減少的主要原因之一；耕地和水域的ESV較高，耕地ESV近期增長(zhǎng)較快，水域ESV波動(dòng)較大，研究期間呈減少趨勢(shì)；草地ESV總體上保持相對(duì)穩(wěn)定的狀態(tài)；因未利用地面積較少，各期未利用地的ESV對(duì)于大余縣總ESV而言幾乎可以忽略不計(jì)；建設(shè)用地持續(xù)加速擴(kuò)張，導(dǎo)致減少用地ESV大量減少，是大余縣總ESV不斷減少的另外一個(gè)主因。

表4 2000-2020年大余縣土地利用類型ESV及變化Table 4 ESV and changes of land use types in Dayu County from 2000 to 2020

（2）生態(tài)系統(tǒng)單項(xiàng)服務(wù)價(jià)值變化

大余縣生態(tài)系統(tǒng)服務(wù)以調(diào)節(jié)服務(wù)為主，其次為支付服務(wù)，兩者價(jià)值分別占68.5%、23%左右。2000-2020年，調(diào)節(jié)服務(wù)減少8120.44萬元，供給服務(wù)減少4431.31萬元，支付服務(wù)減少1911.74萬元，文化服務(wù)價(jià)值減少423.58萬元。除食物生產(chǎn)價(jià)值增長(zhǎng)外，其他二級(jí)生態(tài)系統(tǒng)服務(wù)類型價(jià)值均減少，其中水資源供給、氣體調(diào)節(jié)、氣候調(diào)節(jié)、凈化環(huán)境、水文調(diào)節(jié)減少規(guī)模超過1500萬元。研究期間，因建設(shè)用地及耕地增加，水資源供給服務(wù)價(jià)值變化率達(dá)到62.038%，其他單項(xiàng)生態(tài)系統(tǒng)服務(wù)類型價(jià)值變化率不大（表5）。

表5 2000-2020年大余縣單項(xiàng)ESV變化Table 5 Changes of individual ESV in Dayu County from 2000 to 2020

（3）基于格網(wǎng)尺度的生態(tài)系統(tǒng)服務(wù)價(jià)值空間聚集特征變化

為了確定大余縣總ESV的空間分異和聚類特征，參考相關(guān)格網(wǎng)構(gòu)建的文獻(xiàn)[12,13,21,27]以及多次實(shí)驗(yàn)的基礎(chǔ)上，通過ArcGIS10.5軟件在大余縣行政區(qū)上構(gòu)建500m×500m的漁網(wǎng)并計(jì)算各個(gè)格網(wǎng)的生態(tài)系統(tǒng)服務(wù)價(jià)值，利用GeoDa軟件在格網(wǎng)上建立空間鄰接權(quán)重，進(jìn)行空間自相關(guān)分析。

大余縣總ESV的全局Moran′s I均為正值，在GeoDa中隨機(jī)化999次置換后，P值均為0.001，Z得分均大于2.58，說明大余縣總ESV在空間上呈聚集分布，表現(xiàn)出顯著的空間正相關(guān)性。全局Moran′s I呈現(xiàn)先升后降再升的趨勢(shì)，且在2020年達(dá)到最大值，為0.637，表明大余縣總ESV在格網(wǎng)尺度下的空間集聚性有逐漸增強(qiáng)趨勢(shì)（表6）。

表6 2000-2020年大余縣ESV全局空間自相關(guān)對(duì)比Table 6 Comparison of global spatial autocorrelation of ESV in Dayu County from 2000 to 2020

為進(jìn)一步確定大余縣總ESV聚集區(qū)域的位置和空間差異特征，利用LISA分布圖（圖3）可以清楚表達(dá)局部地區(qū)的空間自相關(guān)。2000-2020年，大余縣總ESV的高高（HH）集聚區(qū)主要分布在大余縣西側(cè)山區(qū)，以及東部北側(cè)和南側(cè)的山區(qū)，符合大余縣三面環(huán)山、林地包夾的地理環(huán)境特點(diǎn)；大余縣總ESV的低低（LL）集聚區(qū)主要分布在大余縣東部中間，該區(qū)域地勢(shì)平坦，有利于人類活動(dòng)，是建設(shè)用地的集中分布區(qū)；低高（LH）和高低（HL）集聚區(qū)比較少，且呈零星分布。總體上說明大余縣總ESV的空間差異小。

圖3 2000-2020年大余縣總ESV的LISA分布圖Fig.3 LISA distribution map of total ESV in Dayu County from 2000 to 2020

將LISA分布數(shù)據(jù)導(dǎo)入轉(zhuǎn)移矩陣中，得到各期格網(wǎng)變化（表7）。從ESV聚集區(qū)格網(wǎng)數(shù)統(tǒng)計(jì)來看，研究期間，大余縣生態(tài)系統(tǒng)服務(wù)價(jià)值HH集聚區(qū)減少最多，LL集聚區(qū)增加最多，與土地利用類型變化特征吻合。2000-2005年和2005-2010年期間，均為不顯著（NS）區(qū)域轉(zhuǎn)出，且多數(shù)轉(zhuǎn)為HH集聚區(qū)；2010-2015年和2015-2020年期間，均為HH集聚區(qū)轉(zhuǎn)出，多數(shù)轉(zhuǎn)為不顯著區(qū)域。

表7 大余縣ESV各聚集區(qū)格網(wǎng)數(shù)變化Table 7 Variation of grid number in ESV aggregation areas in Dayu County

4 結(jié)論

根據(jù)國家重點(diǎn)生態(tài)功能區(qū)江西省大余縣2000-2020年土地利用覆蓋數(shù)據(jù)，運(yùn)用轉(zhuǎn)移矩陣分析土地利用變化特征，以改進(jìn)后的當(dāng)量表為基準(zhǔn)進(jìn)行地區(qū)修訂，估算了大余縣多期生態(tài)系統(tǒng)服務(wù)價(jià)值，并在500m×500m格網(wǎng)尺度下，研究大余縣生態(tài)系統(tǒng)服務(wù)價(jià)值空間聚集特征，得到以下結(jié)論：

（1）研究期間大余縣林地、水域和未利用地面積呈減少態(tài)勢(shì)，建設(shè)用地、草地和耕地呈增加態(tài)勢(shì)，且林地和建設(shè)用地面積變化明顯；林地、建設(shè)用地、耕地和水域之間轉(zhuǎn)換頻繁，以林地轉(zhuǎn)變成耕地和建設(shè)用地，耕地轉(zhuǎn)變成建設(shè)用地和水域?yàn)橹?。作為國家重點(diǎn)生態(tài)功能區(qū)，經(jīng)濟(jì)發(fā)展受到限制，一方面，應(yīng)轉(zhuǎn)變傳統(tǒng)的發(fā)展模式，發(fā)展生態(tài)旅游業(yè)，加強(qiáng)林地和耕地保護(hù)，嚴(yán)格執(zhí)行征占用林地制度，守住林地資源紅線和耕地紅線；另一方面要嚴(yán)格控制建設(shè)用地的規(guī)模，充分挖掘城鎮(zhèn)內(nèi)部建設(shè)用地潛力，推進(jìn)城鄉(xiāng)增減掛鉤，引導(dǎo)國家重點(diǎn)生態(tài)功能區(qū)綠色發(fā)展。

（2）大余縣生態(tài)系統(tǒng)服務(wù)價(jià)值總體上呈現(xiàn)下降的趨勢(shì)，林地減少和建設(shè)用地?cái)U(kuò)張是生態(tài)系統(tǒng)服務(wù)價(jià)值減少的重要因素，尤其林地對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值有著顯著放大作用，是主要貢獻(xiàn)因子和敏感因子。各單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值量的比例結(jié)構(gòu)比較穩(wěn)定，一級(jí)類中調(diào)節(jié)服務(wù)價(jià)值最高，供給服務(wù)價(jià)值最小。從二級(jí)類來看，氣體和水文調(diào)節(jié)價(jià)值在大余縣生態(tài)系統(tǒng)服務(wù)價(jià)值中占據(jù)了主導(dǎo)地位。研究期間水資源供給功能服務(wù)價(jià)值大量減少，今后應(yīng)嚴(yán)格保護(hù)生態(tài)公益林及水源涵養(yǎng)林等生態(tài)環(huán)境脆弱地區(qū)的林地，加快實(shí)施荒山造林、封山育林等重點(diǎn)工程，進(jìn)一步鞏固退耕還林成果，擴(kuò)大有林地面積，從而避免快速城鎮(zhèn)化和頻繁人為活動(dòng)對(duì)區(qū)域生態(tài)系統(tǒng)造成負(fù)面影響。

（3）大余縣總ESV在空間上呈聚集分布，表現(xiàn)出顯著的空間正相關(guān)性。且莫蘭指數(shù)總體上呈上升趨勢(shì)，表明大余縣生態(tài)系統(tǒng)服務(wù)價(jià)值空間集聚性在進(jìn)一步增強(qiáng)。2000-2020年，大余縣總ESV的高高集聚區(qū)主要分布在大余縣山區(qū)丘陵，低低集聚區(qū)主要分布在大余縣東部盆地，低高和高低集聚區(qū)面積較少，且呈零星分布。對(duì)于高高集聚區(qū)內(nèi)生態(tài)環(huán)境要繼續(xù)維持，強(qiáng)化聚集區(qū)之間的連通性；低低聚集區(qū)周邊的區(qū)域是重點(diǎn)生態(tài)環(huán)境監(jiān)測(cè)區(qū)，需要引導(dǎo)其向高高集聚區(qū)轉(zhuǎn)變。

以改進(jìn)后的當(dāng)量因子為基礎(chǔ)，結(jié)合大余縣的實(shí)際進(jìn)行修訂，實(shí)現(xiàn)了區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值更準(zhǔn)確的測(cè)算。在格網(wǎng)尺度下分析大余縣的生態(tài)系統(tǒng)服務(wù)價(jià)值空間變化，突破了把行政區(qū)當(dāng)作評(píng)價(jià)單元的研究模式，顯化了大余縣的生態(tài)系統(tǒng)服務(wù)價(jià)值空間關(guān)系。本研究也存在一些不足：一是研究以遙感影像解譯結(jié)果（一級(jí)地類）為基礎(chǔ)，今后可采用高分辨率遙感影像解譯二級(jí)地類，研究結(jié)果將更精確；二是研究時(shí)間跨度20年，只取2019年糧食價(jià)格作為計(jì)算依據(jù)，若能考慮時(shí)間價(jià)值的影響，更能凸顯生態(tài)系統(tǒng)服務(wù)價(jià)值的規(guī)律；三是本研究?jī)H在格網(wǎng)尺度下分析重點(diǎn)生態(tài)功能區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值變化還不夠，今后考慮景觀格局等其他變化對(duì)價(jià)值的影響。