生态系统服务(ecological services)是人类能够从生态环境中得到的各种恩惠，包括有形物质产品和无形服务两方面，可分为供给服务、调节服务、支持服务和文化服务，而这些是人类生存和发展的必需品和生活质量的保证^[1-3]。生态系统服务价值(ecological service values，ESV)评估是用货币的形式评估生态系统的能力，能直观反映生态系统服务效益的变化^[4-5]。现在生态服务价值的研究已成为生态经济学和环境经济学的研究热点，并且国内外对生态服务价值评估的方法已逐步成型^[6-7]。国内外学者已对城市^[8]、森林^[9]、草地^[10]、流域^[11]、湿地等^[12]不同生态系统的生态服务价值进行了大量深入研究。土地利用/土地覆盖变化(LUCC)是人类社会经济活动行为与自然生态过程交互衔接的纽带，能够引起生态服务价值的改变^[13-14]。土地利用与生态系统服务相互影响，相互制约。土地利用活动影响着各类生态系统类型、面积以及空间分布的变化，是生态系统服务功能变化的重要驱动力^{[13, 15]}。现有的土地利用变化与生态服务价值的相关研究多是对于土地利用变化下生态服务价值时空演变趋势的研究，主要分析土地利用由某一地类转换为另一地类时其本身变化带来的生态服务价值变化。而忽略了这一土地类型的变化对周围土地生态服务价值也存在着影响，不同的土地类型对周围土地生态服务价值影响是不同的，所以当土地类型发生变化时，会连带影响其周围土地的生态服务价值。本研究将土地类型转换对其周围土地生态服务价值产生影响的现象称为土地利用变化与生态服务价值的遥相关。分析这种遥相关，对研究地利用变化与生态服务价值的关系、优化土地结构、制定生态补偿机制以及促进社会经济与生态环境的协调可持续发展等具有重要意义。为此，本文以湖南省韶山市为研究对象，通过遥感解译提取1985年及2020年韶山市土地利用数据。基于该数据研究土地利用变化对周围土地生态服务价值的影响，提出土地利用变化与生态服务价值的遥相关概念，构建土地利用变化与生态服务价值的遥相关计算模型，分析遥相关影响因子，探讨土地利用变化对周围土地生态服务价值的影响，为土地利用变化与生态服务价值的相关研究提供依据，也为生态服务价值评估体系的确定提供参考。

1 研究区概况

韶山市位于湖南省湘潭市，距长沙市约120 km。地理坐标为东经112°23′52″—112°38′13″，北纬27°51′40″—28°1′53″。该市现辖2镇2乡(清溪镇、银田镇、韶山乡及杨林乡)，总人口约12万，行政区域面积约247 km²。该区属于湘中低山丘陵区，地势西高东低，地貌轮廓以韶峰山脉和韶河、石狮江两水为骨架，构成西部隆起，往东及东南倾斜的地势(图 1)。该区地处亚热带湿润气候区，四季分明，冬冷夏热、夏热期长，严寒期短，年均温16.7 ℃，年均降水1 358 mm。在社会经济发展的过程中，韶山市立足于旅游发展，是国家AAAAA级旅游景区、全国爱国主义教育基地和中国优秀旅游城市。

2 数据处理与研究方法 2.1 数据获取与数据处理

选择韶山市1985和2020年遥感卫星数据(Google Earth下载)为基础数据。运用ENVI 5.2和ArcGIS 10.2软件对研究区不同时期的遥感影像进行图像镶嵌、裁剪及几何校正等预处理，参照国家土地利用分类/土地覆盖研究系统及研究区的特点，将研究区土地利用类型划分为林地、耕地、草地、水体、建设用地和未利用地等6类^[16]。在此基础上运用监督分类法对遥感影像进行判读解译，由此得到1985年和2020年的土地利用数据(附图7)。粮食产量、粮食价格、旅人综合收入及旅游人次等数据主要来自源于1985—2020年《韶山市国民经济和社会发展统计公报》《湖南省统计年鉴》及《湘潭市统计年鉴》。

2.2 土地利用变化分析方法

(1) 土地利用转移矩阵：是马尔科夫模型在土地利用变化方面的应用，是分析研究区各土地利用类型间相互转换情况的一种比较好的方法。

(2) 土地利用动态度：主要反映某研究时间段内土地利用类型数量的变化情况^[17]。单一土地动态度可表示研究区在一定时间段内某种土地类型的数量变化情况，不仅反映出其变化的剧烈程度，还反映了土地利用变化速度，其计算公式^[18]为：

$ K=\frac{U_{b}-U_{a}}{U_{a}} \times \frac{1}{T} \times 100 \% $

(1)

式中：K为研究时间段内某一土地利用类型动态度; U_a, U_b分别为研究时间段初期和末期该土地类型面积; T为研究时段, 设定为年。

2.3 生态服务价值测算方法

1997年Contanza等^[19]提出的生态服务价值评价体系是目前最为广泛的方法。而谢高地等^[20]依据Contanza的评估方法，从我国实际情况出发，制定了中国陆地生态服务价值当量因子表。并确定1个生态服务价值当量等于当年全国平均粮食单位面积市场价值的1/7^[21]，计算公式为：

$ \mathrm{VC}=\frac{1}{7} \times T_{a} \times T_{b} $

(2)

式中：VC表示单位当量因子的价值; T_a，T_b分别表示为研究区平均粮食产量和平均粮食价格。

本文采用2015年最新改进的单位面积生态服务价值当量^[22]，结合韶山市风景名胜区的实地情况，主要参考风景名胜区的旅游特性对中国陆地生态系统单位面积生态服务价值表中的文化服务指标进行调整。对于建设用地的生态服务价值，许多学者都根据谢高地等的研究成果将建设用地的生态服务价值赋值为0^[23]，本文将采用这一赋值作为建设用地的值。由韶山市国民经济和社会发展统计公报及相关论文等资料可以得到研究时间段1985—2020年韶山市平均粮食产量6 210.9 kg/hm²和平均粮食价格3.69元/kg。由公式(2)计算出韶山市一个当量生态服务价值为3 274元/hm²。由此得到韶山市各土地利用类型的单位面积生态服务价值表(表 1)。

利用遥感解译提取韶山市的土地利用数据，参考Costanza^[19]的生态服务价值计算公式，结合韶山市实地情况修改后的谢高地等提出中国陆地生态系统单位面积生态服务价值表，来计算韶山市不同地类的生态服务价值，其计算公式为：

$ \mathrm{ESV}_{i}=\sum A_{i} \cdot \mathrm{VC}_{i} $

(3)

式中：ESV_i表示为i地类的生态服务价值; A_i表示i地类的面积; VC_i表示i类地类的单位面积生态服务价值。

2.4 土地利用变化与生态服务价值的遥相关测算方法 2.4.1 土地利用变化与生态服务价值的遥相关概念

所有土地利用类型都不是独立存在的，它们相互作用相互影响。土地类型发生变化，不仅影响了本身土地的生态服务价值，还对其周围土地的生态服务价值产生了影响。例如，某地类由水体转换为林地，对周围土地生态服务价值的影响由水体的影响转换为林地的影响，水体消失其周围耕地、草地等可能出现因干旱而导致减产、干枯等现象，进而影响土地的生态服务价值；林地出现固化土地，减少水土流失，与周围耕地、草地等构成新的气体及气候调节，进而影响其生态服务价值。土地利用通过地类的转换进而影响其周围土地生态服务价值，称这种关系为遥相关。基于此提出遥相关概念：即所有土地利用类型都不是独立存在的，他们相互影响相互作用，每一土地类型都对周围土地生态服务价值产生影响，而当土地利用类型发生转换时，其对周围土地生态服务价值的影响也发生变化，称这种现象为土地利用变化与生态服务价值的遥相关。

2.4.2 遥相关计算模型构建

基于生态服务的基本定义，综合考虑土地利用变化与土地生态服务价值变化的关系、土地利用类型转换情况及相隔距离等，构建土地利用变化与生态服务价值的遥相关计算模型：土地利用变化与生态服务价值的遥相关影响价值=影响因子×周围土地生态服务价值。

(1) 影响因子。遥相关的影响因子是反映土地利用变化对周围土地生态服务价值的影响程度，与土地利用的转换、面积及距离等有关。影响因子的计算主要依据转换前后的土地类型对其周围土地提供的生态服务价值当量及周围土地本身的生态服务价值当量。在谢高地等提出的生态服务中，主要由调节服务这一服务功能对外构成循环并产生影响，由此选取调节服务功能的价值当量作为各土地类型对周围土地提供的生态服务价值当量，以此为基础数据得到各个距离下的土地类型对周围土地提供的生态服务价值当量。由于土地类型面积的大小对周围土地的影响不同，所以选取间隔面积斑块个数为距离单位，计量单位为个，距离为间隔0，1，2，…，n个面积斑块。其中面积斑块是指i类土地利用类型转换为j类土地利用类型的转换面积，各地类转换的面积斑块各不相同。如某一水体转换为林地的面积大小为X m²，则该地类转换的面积斑块大小为X m²，其距离选取为间隔0，1，2，…，n个面积斑块(当距离单位由面积斑块个数转换成米时，一个面积斑块代表的距离大约为$\sqrt X $ m)。由此构建土地利用变化与生态服务价值的遥相关影响因子计算公式：

$ \mu=\left(\mathrm{ESVE}_{j_{n}}-\mathrm{ESVE}_{i_{n}}\right) \cdot \mathrm{ESVE}^{-1} $

(4)

式中：μ表示i类土地利用类型转换为j类土地利用类型在距离为n斑块(土地类型变化的面积斑块)下的影响因子; ESVE_jn，ESVE_in分别为距离为n斑块下的转换后和转换前土地类型对其他土地提供的生态服务价值当量; ESVE表示为各土地利用类型的平均生态服务价值当量。

(2) 遥相关模型。遥相关模型主要阐述土地利用变化对周围土地生态服务价值的影响，由遥相关的影响因子和周围土地的生态服务价值构建其计算公式：

$ \mathrm{RESV}=\sum\left(\mu \cdot \mathrm{ESV}_{n}\right) $

(5)

式中：RESV表示土地利用变化对周围土地生态服务价值的遥相关影响值; μ表示i地类转换为j地类在距离为n个的土地利用变化对生态服务价值的影响因子; ESV_n表示在距离为n个面积斑块下的周围土地的平均生态服务价值, 由i地类转换为j地类在距离为n个面积斑块下周围土地的面积和平均的生态服务价值当量(由表 1得到), 据生态服务价值计算公式(3)计算得到。

土地利用转换对本身土地生态服务价值的影响因子为1，即μ=1，ESV_n由该土地面积和生态服务价值当量计算得到，此时得到的生态服务价值当量为该土地类型转换而带来的生态服务价值。而当将距离足够远时，土地利用转换对本身土地生态服务价值的影响因子为0，即μ=0，土地利用转换对生态服务价值的影响为0。由此可知该公式的设定基本符合要求。

3 结果与分析 3.1 土地利用变化分析

1985—2020年韶山市土地利用发生了明显变化，通过表 2对土地利用变化和转换趋势进行分析。从各地类的转换数量来看，林地主要转出为草地和耕地，其中大约有4%的林地转换为草地，3%的林地转换为耕地，林地的转入来源以耕地为主。耕地和林地的转换幅度较大，在耕地的转出中有48%的都转出为林地，这主要由于研究区的退耕还林政策。草地占总面积的29%，是韶山市风景区第二大土地利用类型，主要转换为林地和耕地，其转入来源主要以耕地为主，其中30%的耕地转入草地，还有未利用地，这说明韶山市的景区规划加大了对土地的利用。受韶山市风景名胜区地形地貌的影响，建设用地的面积所占比例较小，约为3.6%，其转换的主要来源以耕地和草地为主。其他土地类型，如林地、水体等由于其开发难度和生态保护问题，转入建设用地的面积比较少。从单一土地利用动态来看，由于时间跨度比较长，各地类的动态度差异不大。其中为正的有林地、草地和建设用地，研究期内其面积净增。增加幅度比较大的为林地，表明韶山市风景名胜区退耕还林和保护生态林业方面功效显著；其次为建设用地，表明韶山市风景名胜区在城市扩张和景区建设等方面加大了对土地的利用。单一土地利用动态度为负的有耕地、水体和未利用地，各土地类型面积有不同程度的减少，其中减少幅度最大的为未利用地，主要由于城市建设的侵占。

3.2 生态服务价值分析

根据生态服务价值计算公式(3)，计算得到韶山市1985年、2020年的生态服务价值(表 3)。由表 3可知1985—2020年韶山市的生态服务价值是增加的，由12.82亿元上升到17.11亿元，35 a里生态服务价值增长了4.69亿元。其中林地生态服务价值最高，随着研究期内韶山市退耕还林政策的实施，林地面积迅速增加，面积比例从25%上升到45%，生态服务价值比例从51%上升到68%；其次为草地，面积比例约为25%，生态服务价值的比例在23%左右。林地和草地这两者对生态服务价值的贡献率超过了80%，这主要是由于林地和草地在韶山市风景名胜区的生态系统中占有重要地位。而耕地的面积则随着“退耕还林”措施的实施迅速减少，面积比例由42%下降到20%，其生态服务价值也在急剧降低，比例从14%下降到5%；同上，水体和未利用地比例就比较小。地类相互转换的过程，必然会改变土地原有的生态结构和功能，进而影响其生态服务价值。

通过附图8分析生态服务价值变化的空间分布，可以看出生态服务价值剧烈增加主要发生在城市东北方向和东南方向，这是由于退耕还林政策，林地和耕地进行了大量转换；生态服务价值急剧降低主要发生在城市中心和西南方向，这是由于城市建设和景区规划，建筑用地急剧增加。总体来说为了韶山市的可持续发展，应重点保护对生态服务价值贡献率较大的林地、草地和耕地等土地利用类型。

3.3 土地利用变化与生态服务价值的遥相关分析 3.3.1 遥相关影响因子分析

根据公式(4)计算各土地类型转换的遥相关影响因子(表 4)，由于数据较多选取关键数据(各地类转换的遥相关临界距离下的影响因子数据)列举。韶山市各地类转换对其周围土地生态服务价值的影响不同，但各地类转换都在其他地类转换为水体这一地类时影响因子达到最大正值，其次为林地，主要由于水体和林地在生态系统中占有重要地位，其中水体对调节服务功能的贡献最大，林地次之。当土地类型转换为建设用地时影响因子皆为最小负值，主要由于建设用地对生态服务价值的影响最小。

各地类转换的影响因子都随着距离的增加而逐渐减小，直到减小为0。当土地类型转换为林地、水体时其影响因子在距离大约为60，100个面积斑块左右减小为0，土地类型转换为草地、耕地、建设用地及未利用地时其影响因子在距离大约为50个面积斑块减小为0。例如，面积为X m²的土地由其他地类转换为林地或水体时，其影响因子在距离大约为60 $\sqrt X $ m或100 $\sqrt X $ m时减小为0，也就是说这种土地类型转换对周围土地生态服务价值的影响距离大约为60 $\sqrt X $ m或100 $\sqrt X $ m；同样面积为X m²的土地由其他地类转换为草地、耕地、建设用地或未利用地时，其影响因子在距离大约为50 $\sqrt X $ m时减小为0，也就是说这种地类转换对周围土地生态服务价值的影响距离为大约为50 $\sqrt X $ m。水体的影响距离最大，这是由于水体调节服务功能最为强劲，其次为林地，最后草地、耕地、建设用地及未利用地的影响距离基本持平，这是由于草地和耕地的调节服务功能差别不大，故他们的影响距离基本相同，而建设用地和未利用地则是因为未利用地的调节服务功能非常小与建设用地的调节服务功能(0)差别不大，当其他地类转换为建设用地或未利用地时，其对周围土地生态服务价值的影响为其他地类调节服务功能价值当量的负值，其平均影响值与草地和耕地的影响值的绝对值基本相同，故建设用地及未利用地的影响距离与草地和耕地基本持平。因此为了保证研究区生态循环的健康，应注重水资源和森林资源的管理。

3.3.2 遥相关影响分析

通过遥相关计算模型(公式5)，得到韶山市土地利用变化对周围土地生态服务价值的遥相关影响价值表(表 5)。由于数据较多，本文选取遥相关影响价值较大的前10个数据列举。土地利用变化对周围土地生态服务价值的遥相关影响值为4.67亿元，所以1985—2020年韶山市生态服务价值由12.66亿元最终上升到22.01亿元，总体增加了9.35亿元。其中土地类型转换为林地、草地和水体时，遥相关影响值为正，主要由于林地、草地和水体的调节服务功能对周围土地生态系统的贡献率较大。

土地类型转换为耕地、建设用地及未利用地时，遥相关影响值为负，主要由于耕地、建设用地及未利用地的调节服务功能对周围土地生态系统的贡献率较小。其中林地的遥相关影响值最大，为4.78亿元，主要是由于退耕还林政策，大量耕地转入林地，林地面积迅速增加，且对调节服务功能的贡献率较大；建设用地的遥相关影响值最小，为-0.24亿元，主要由于建设用地面积变化比较小，对调节服务功能的贡献率也最小。各地类转换的交互影响，不仅影响其本身土地的生态服务价值，还影响其周围土地的生态服务价值，所以在规划土地利用结构时，应综合考虑各土地类型转换的交互影响，制定合理的土地利用政策，实现韶山市经济建设与生态环境的协调可持续发展。

4 结论

(1) 通过对土地利用分析，研究期内韶山市土地利用面积发生不同程度的变化，土地利用转换主要集中在林地、耕地和草地之间。林地和耕地的转换最为剧烈，林地面积急剧增加，耕地面积急剧降低。林地、草地和建设用地的土地动态度为正，耕地、水体和未利用地的土地动态度为负。说明韶山市风景名胜区在城市扩张和景区建设等方面加大了对土地的利用，应根据韶山市生态环境和社会需求合理制定土地利用政策，加强对土地利用的宏观管理，实现土地利用结构的优化。

(2) 通过对生态服务价值分析，1985—2020年韶山市的生态服务价值由12.66亿元上升到17.34亿元，35 a里生态服务价值增长了4.68亿元。林地的生态服务价值最大，约占总值的68%，且随林地面积的扩张而增加。由于韶山市城市生态文明建设力度的加大和“退耕还林”政策的进一步实施，韶山市林地及草地用地面积不断增加，促使韶山市生态服务价值增加。表明林地是维持韶山市风景名胜区生态系统服务价值的关键所在，需加强森林资源管理，促进森林系统的生态效应的有效发挥，实现当地经济建设和生态环境的和谐发展。

(3) 土地利用变化与生态服务价值的遥相关分析，研究期内韶山市遥相关影响值为4.67亿元，所以韶山市生态服务价值由12.66亿元最终上升到22.01亿元，总体增加了9.35亿元。土地类型转换为林地、草地和水体时，遥相关影响值为正，转换为耕地、建设用地及未利用地时，遥相关影响值为负。地类转换的遥相关影响距离各不相同，水体的影响距离最大，其次为林地，而其余地类的影响距离基本持平。表明水体的对调节服务功能的贡献率最大，为了保证研究区生态循环的健康，应注重水资源管理。各土地类型相互作用相互影响，土地利用变化不仅影响其本身土地的生态服务价值，还影响其周围土地的生态服务价值。所以在规划土地结构时，应结合各地类转换的交互影响，合理调整土地资源，制定适宜的土地利用政策，实现当地经济与环境的协调可持续发展。

对于土地利用转换对周围土地生态服务价值的影响本文只进行了初步的探讨，还有很多问题需要进一步的研究分析。例如，构建土地利用变化与生态服务价值的遥相关计算模型时，未考虑周围土地的具体利用类型、土地利用类型转换复杂的叠加状态及土地利用变化与生态服务价值的遥相关计算模型的具体验证等。