基于PCA的龙口市土地可持续利用评价*

苏晨晨，周 奥，潘玉翠，王瑷玲

(山东农业大学资源与环境学院，泰安 271018)

0 引言

土地资源是人类赖以生存发展的基础，土地资源可持续利用是区域可持续发展的重要保障，土地资源的可持续利用程度直接关系到区域的社会经济发展与生态文明建设[1]。21世纪以来，人口快速增长、经济迅猛发展，人地矛盾日益突出，生态环境退化[2]，如何实现土地可持续利用和社会经济可持续发展成为社会热点。土地可持续性利用评价是土地可持续利用的基础，只有了解土地资源可持续利用现状及变化，发现存在的问题，才能科学地制定土地可持续利用政策和区域发展规划与可持续发展战略，实现区域社会、经济与生态的协调发展。

随着持续利用土地思想的提出，学者们广泛开展了土地可持续利用有关研究，主要从土地可持续利用评价指标体系和评价方法两个方面展开。评价指标体系构建较多借鉴了PSR体系[3]、FELSM框架[4-5]以及社会—经济—生态三分体系[6-8]包括城镇化水平、人口密度、人均收入、地区生产总值、林地覆盖率、化肥、农药使用量等指标。我国人多地少、耕地相对不足，土地可持续利用应考虑节约集约利用土地，在追求社会经济发展的同时，尽量减少建设用地的消耗，实现人口、经济发展与土地资源消耗的协调，但目前的评价指标体系没有考虑。评价方法有综合评价法[2-4]、三角模型评价法[5,9]、灰色关联分析法[10]、生态足迹评价法[11]等，方法较先进，但较多评价指标难以避免相互之间共线性，指标权重确定多采用层次分析法[1，]，带有一定的主观性。主成分分析法是采用数理统计方法，将多个指标转化为几个综合指标进行评价的方法，避免了评价指标间的共线性和人为确定指标权重的主观性，在土地可持续利用评价中应用较少。龙口市是山东省经济第一县市，其土地利用是否可持续直接影响其社会经济发展的可持续，对其他县级市的可持续发展也有启示作用。因此文章选取龙口市为研究区域，引入土地节约集约利用评价指标体系中能够反应社会经济发展与建设用地消耗协调程度的人口、经济用地弹性系数，运用主成分分析法，对2004—2015年龙口市土地可持续利用状况进行动态评价分析。

1 研究区概况与数据来源

1.1 研究区概况

龙口市位于东经120°13′～120°44′，北纬37°27′～37°47′，地处山东东北部、胶东半岛的西北部，辖8个镇、5个街道办事处、1个省级经济开发区、1个省级高新技术产业园区、1个滨海旅游度假区，总面积901km2，总人口69.79万人。2015年龙口市地区生产总值达1 041.306 7亿元，其中二、三产业产值1 005.151 2亿元，占总产值的96.53%，人均地区生产总值14.919 6万元。2015年全国县域经济与县域基本竞争力百强县中龙口市排名第十位、山东省第一位。龙口市社会经济和城镇化水平发展迅速， 2004—2015年全市地区生产总值平均年增长率为24.47%，城镇化水平由2004年的26.21%上升到2015年的59.61%。龙口市在发展经济的同时，比较重视生态环境保护，获得国家可持续发展实验区、全国生态保护与建设示范区等称号，但还是对土地资源利用与生态环境产生了一些不利影响，人类活动与土地资源利用的矛盾比较突出，对实现区域可持续性发展造成了一定威胁。因此，进行龙口市土地可持续利用评价，揭示其可持续利用水平及变化，对其社会经济可持续发展具有重要意义。

1.2 数据来源

数据主要包括社会经济数据和土地利用数据。基于数据可获取性，研究时段确定为2004—2015年。社会经济数据主要来源于《龙口市统计年鉴(2004—2015年)》《烟台市统计年鉴(2004—2015年)》和《山东省统计年鉴(2004—2015年)》。土地利用数据主要来源于龙口市国土资源局的土地利用变更调查数据。因2009年“二调”前后的土地利用分类体系发生变化，为保持土地利用数据的一致性，将2009年以前的土地利用数据以“二调”土地利用分类进行调整、统一数据口径。

2 主成分分析原理与步骤

2.1 主成分分析原理

主成分分析法(Principal component analysis，简称PCA)是利用降维的思想，在损失很少信息的前提下把多个指标转化为几个综合指标的统计方法[13]。在进行土地可持续利用评价时，为更全面地体现土地可持续性利用，需要考虑社会、经济和生态3类因素，每类因素下又包括多个因子，因子间可能有信息重叠，引起指标之间的共线性。主成分分析通过将原始指标的线性组合来解释原始指标的绝大部分信息，且对指标赋权，计算综合评价值，可降低评价指标的数据维数，提高评价效率[14-15]和评价结果的客观性。

2.2 主成分分析步骤

(1)建立m个年份n个指标的原始数据矩阵X。

X=(xij)m*n，i=1， 2…，m，j=1， 2…，n； 其中，i代表研究年份，j代表评价指标。

(2)标准化处理。为消除指标量纲差异的影响，采用标准差标准化方法，对原始值进行标准化处理，计算公式如下：

(1)

(3)计算指标的相关系数矩阵Rn*n，求它的特征值λ、特征向量μ。

根据特征方程|R-λ|=0计算特征值λ，同时求得特征向量μ。

(4)计算方差贡献率与累积方差贡献率。取特征值大于1、累积方差贡献率在80%以上的作为主成分。

(5)计算主成分载荷。因子载荷量是主成分与原始指标的相关系数。

(6)计算各主成分得分:

(2)

(7)以特征值方差贡献率为权数，采用线性回归方法，计算综合得分，进行综合评价:

Zi=∑FiPn

(3)

(4)

式(3)、(4)中，Zi为第i年的综合得分，Pn为第n个主成分因子的方差贡献率；Z′为百分制化后的综合分值，Zmax为百分制化前各因素得分的最大值，Zmin百分制化前各因素得分的最小值。

3 土地可持续利用评价指标体系构建

3.1 土地可持续利用内涵分析

土地可持续利用思想最早于1990年国际土地持续利用系统研讨会上提出，但对于具体内涵，还没有统一的定义。1987年世界环境与发展委员会在《我们共同的未来》报告中提出可持续发展包括生态、经济、社会3个关键要素[3]。1993年FAO颁布的《持续土地管理评价大纲》认为土地可持续利用包括生产性、安全性、保护性、可行性和接受性等5个方面[16]。也有学者提出土地持续利用就是协调人类当代与后代之间的经济、社会与环境方面的需求，同时维持与加强土地资源质量[17]。中国人多地少、耕地相对不足，土地应节约集约利用，尤其建设用地扩张应与社会经济发展协调。综合可持续发展和土地可持续利用观点，考虑土地节约集约利用，该文认为土地可持续利用是指通过土地资源的节约集约、合理利用，土地利用各子系统相互协调，促进社会、经济、生态可持续发展。

3.2 土地可持续利用评价指标体系构建

构建评价指标体系是科学评价的关键。土地节约集约利用应成为土地可持续利用评价的重要内容，土地可持续利用评价应考虑人口增加、经济发展对建设用地占用的弹性，包括人口—建设用地弹性和经济—建设用地弹性。人口—建设用地弹性系数是在分析一定时期区域人口、建设用地的变化量、变化幅度相关关系的基础上，判断区域人口变化与建设用地增加的协调程度，而经济—建设用地弹性系数则是在分析区域地区生产总值、建设用地的变化量、变化幅度的基础上判断经济变化与建设用地增加的匹配程度，二者在一定程度上反映了建设用地消耗与自身社会、经济的发展协调程度。因此在土地可持续利用内涵界定基础上，参考有关土地可持续利用评价和土地节约集约利用评价等评价指标体系，并考虑相关指标数据的可获取性，结合龙口市特点，引入建设用地弹性系数，从社会发展、经济发展和生态环境3个方面，选择城镇化水平等6个社会因子、地区生产总值等7个经济因子和林地覆盖率等6个生态因子，构建土地可持续利用评价指标体系(表1)。

指标体系中，城镇化水平、人口密度、人口用地弹性系数、教育财政支出比例等反映了社会发展状况。地区生产总值、固定资产投资、经济用地弹性系数、经济增长率、收入等反映了经济发展状况。林地覆盖率、城镇生活污水处理率、单位耕地化肥负荷等反映了生态环境状况。社会发展、经济发展和生态环境状况共同体现了区域土地资源的可持续利用水平。与FAO土地可持续利用的生产性、安全性、保护性、可行性和接受性相比，社会发展指标主要体现了可行性、接受性，经济发展指标体现了生产性，生态环境体现了安全性和保护性。该研究创新提出的用地弹性指标，也体现了土地可持续利用的保护性、接受性。

4 基于PCA的土地可持续利用评价

土地利用可持续包括社会发展可持续、经济发展可持续和生态环境可持续3方面。为了了解各因素可持续水平，首先分社会发展、经济发展和生态环境3个因素，分别评价2004—2015年龙口市土地利用可持续水平。然后综合评价2004—2015年龙口市土地可持续利用水平。

4.1 指标标准化

通过数据收集与初步计算，得到2004—2015年各指标原始值。采用式(1)进行指标标准化，得到各指标标准化值(表2)。

表1 龙口市土地可持续利用评价指标体系

目标层准则层指标层指标计算方法指标属性土地可持续利用社会发展城镇化水平(C1)城镇人口/总人口+人口密度(C2)总人口/土地总面积+人口—建设用地弹性系数(C3)[(t年总人口-t-1年总人口)/t-1年总人口]/[(t年城乡建设用地面积-t-1年城乡建设用地面积)/t-1年城乡建设用地面积]+人均社会消费品总额(C4)社会消费品零售总额/总人口+教育财政支出比例(C5)教育财政支出/财政总支出+就业人口比重(C6)就业人口/总人口+经济发展固定资产投资(C7)—+地区生产总值(C8)—+经济—建设用地弹性系数(C9)[(t年地区生产总值-t-1年地区生产总值)/t-1年地区生产总值]/[(t年建设用地面积-t-1年建设用地面积)/t-1年建设用地面积]+经济增长率(C10)(t年地区生产总值-t-1年地区生产总值)/t-1年地区生产总值+财政总收入(C11)—+农村居民人均纯收入(C12)—+城镇居民人均纯收入(C13)—+生态环境城镇生活污水处理率(C14)—+林地覆盖率(C15)林地面积/土地总面积+单位耕地化肥负荷(C16)化肥折纯量/耕地面积-单位耕地农药负荷(C17)农药使用量/耕地面积-单位耕地薄膜负荷(C18)农用塑料薄膜使用量/耕地面积-工业废气排放达标率(C19)—+

表2 2004—2015年龙口市指标标准化值

指标200420052006200720082009201020112012201320142015C1-2.55-0.39-0.36-0.36-0.48-0.480.660.660.670.810.960.86C2-1.62-0.260.641.651.320.94-1.12-0.74-0.53-0.560.34-0.05C30.112.611.41-0.22-0.87-0.46-0.44-0.44-0.41-0.46-0.23-0.60C4-1.27-1.09-0.97-0.80-0.56-0.310.010.310.630.991.331.73C5-1.27-1.72-1.45-0.24-0.130.110.330.410.990.971.001.00C6-1.33-0.91-0.87-1.15-0.70-0.380.181.071.051.490.840.70C7-1.39-1.16-0.97-0.88-0.52-0.120.360.410.410.871.281.70C8-1.48-1.23-0.96-0.66-0.47-0.240.090.410.721.061.311.46C91.061.871.30-0.86-0.79-0.80-0.94-1.040.440.47-0.50-0.22C101.041.751.400.10-1.20-0.21-0.240.64-0.52-0.64-0.99-1.12C11-1.19-1.05-0.91-0.73-0.58-0.44-0.160.230.621.091.371.75C12-1.29-1.13-0.97-0.77-0.54-0.36-0.060.360.761.151.551.28C13-1.34-1.15-1.01-0.79-0.53-0.29-0.030.440.851.231.301.32C14-2.57-0.53-1.220.030.080.580.600.600.680.690.380.69C15-1.15-0.92-0.80-0.87-0.84-0.890.011.051.071.091.111.13C160.630.29-0.10-0.98-2.670.12-0.180.410.360.391.320.41C171.100.670.861.050.44-0.07-1.37-0.08-2.310.02-0.28-0.02C180.580.130.39-0.95-0.08-2.35-0.770.310.290.260.481.73C190.370.370.370.370.370.37-0.65-3.040.370.370.370.37

4.2 分因素评价

利用SPSS20.0进行主成分分析。对标准化数据分析计算，得到各因素的特征值、方差贡献率和累积方差贡献率，按照累积方差贡献率≥80%、特征值>1的原则，提取主成分(表3)。

表3 分因素特征值与方差贡献率

由表3可知，社会发展因素前两项成分的累积方差贡献率为83.389%，超过80%，将其作为主成分因子； 经济发展因素第一个成分的方差贡献率达到81.856%，将其作为主成分因子； 生态环境因素前3项成分的累积方差贡献率为81.221%，但第3项成分特征值小于1，考虑到主成分需涵盖原始数据的绝大部分信息，所以将前3项作为主成分因子。

选取主成分后，计算各因素因子载荷矩阵(表4)。

表4 分因素因子载荷矩阵

社会发展12经济发展1生态环境123C10.8480.04C70.981C140.734-0.5620.095C2-0.1260.949C80.992C150.9480.1930.051C3-0.658-0.348C9-0.501C160.4720.5790.080C40.950-0.044C10-0.819C17-0.7750.313-0.290C50.9720.141C110.970C180.2870.8100.021C60.937-0.255C120.979C19-0.4340.0790.891C130.980

表4表明，社会发展因子中，城镇化水平C1、人口—建设用地弹性系数C3、人均社会消费品总额C4、教育财政支出比例C5和就业人口比重C6在主成分1中占较高的负荷，而人口密度的影响相对较弱，这表明城镇化水平、居民生活水平的提高对社会可持续发展的影响是突出的。经济发展因子中，固定资产投资C7、地区生产总值C8、经济增长率C10、财政总收入C11和农村居民人均纯收入C12、城镇居民人均纯收入C13有较高的负荷，而经济用地弹性系数的影响相对较小，这表明经济可持续发展对土地投入、产出水平的依赖较高。生态环境因子中，城镇生活污水处理率C14、林地覆盖率C15以及单位耕地农药负荷C17在主成分1中的作用较明显，单位耕地化肥负荷C16、单位耕地薄膜负荷C18在主成分2中有较高负荷，这表明环境污染的处理能力和环境资源的压力水平对生态环境资源持续发展有重要影响。

4.3 综合评价

与分因素评价原理相同，利用SPSS20.0 软件进行处理，根据累积方差贡献率≥80%、特征值>1的原则提取出前4项主成分，其特征值、方差贡献率、累积方差贡献率以及因子载荷矩阵见表5。

表5 综合因素因子载荷矩阵

评价指标主成分1主成分2主成分3主成分4评价指标主成分1主成分2主成分3主成分4C10.857-0.011-0.0660.420C120.9780.1140.1340.000C2-0.121-0.7150.4320.336C130.9860.1010.0980.005C3-0.6140.5390.1070.525C140.801-0.378-0.1430.393C40.9730.1180.16-0.028C150.9360.295-0.073-0.009C50.960-0.207-0.016-0.091C160.3160.749-0.1080.100C60.9300.212-0.1650.071C17-0.650-0.0180.298-0.207C70.9750.0550.078-0.014C180.2730.6490.333-0.365C80.9910.0370.1050.014C19-0.210-0.0700.8600.118C9-0.4810.7310.3380.197特征值11.5962.9461.5341.025C10-0.7550.490-0.3090.225方差贡献率/(%)61.03115.5048.0765.394C110.9570.1530.219-0.040累积方差贡献率/(%)61.03176.53684.61190.005

图1 2004—2015年龙口市土地可持续利用分值

4.3 可持续分值计算

根据各指标的特征向量和标准化数据，采用式(2)计算因素各主成分得分。然后根据因素主成分得分，以因素内各主成分的方差贡献率作为权重，采用式(3)计算因素综合得分。为更直观地反映龙口市评价时期内土地可持续利用状况，将因素综合得分采用式(4)进行百分制化，结果见图1。

5 结果分析

5.1 社会发展可持续水平

图1表明， 2004—2015年龙口市土地社会可持续性水平逐年提高。土地社会可持续水平大致经历3个阶段： 2004—2007年，社会发展可持续水平增长较快， 2004年社会发展可持续水平相对最低， 2007年达到46.57分，年增长率为11.64%； 2007—2009年，社会发展可持续水平发展迟缓，年均增长率为3.83%， 2009年达到58.05分； 2009以后进行快速增长阶段，年增长率为10.32%， 2015年达到相对最高。

分析原始数据可知， 2004—2015年龙口市人均社会消费品消费总额、教育财政支出比例、就业人口比重指标值整体上逐年增加， 2015年龙口市人均社会消费水平达5.571 6万元，教育财政支出比例29.41%，就业人口比重达56.96%，表明龙口市居民消费能力水平逐渐提高，人口就业比率有所提高，且政府大力支持教育事业的发展，使得龙口市土地社会可持续利用水平不断提高。而人口用地弹性系数由2007年的0.30，下降到2009年-0.01，人口呈现下降的趋势， 2007—2009年人口减少了1 340人，人口变化速度明显低于城乡建设用地增长速度，人口变化与城乡建设用地增长不协调，从而造成2007—2009年土地社会发展可持续水平增长缓慢。另外，随着城镇化进程的加快，龙口市城市人口不断增长、农村人口不断减少，但城镇、村庄用地面积却均逐年增加，农村出现“一户多宅”、宅基地闲置等现象，土地集约利用水平低，导致了土地社会可持续发展受到了限制。

5.2 经济发展可持续水平

龙口市土地经济发展可持续水平整体呈上升趋势。2004、2005年土地经济发展可持续水平相对较低，分值分别为0.24分、0分。2005—2008年经济发展可持续水平快速增长，由0分增长至41.35分，年增长率达到10.34%。2008—2009年经济可持续发展平衡，增长不明显。2009年后又进入快速增长时期，年增长率为8.06%， 2015年达到相对最高。

分析发现， 2004年以来龙口市经济持续发展，固定资产投资不断增加，地区生产总值逐年提高， 2015年龙口市固定资产投资总额与地区生产总值分别达到594.785 1亿元、1 041.306 7亿元，整体的经济发展带动了居民收入的提高和土地经济可持续利用水平提高。2008、2009年受国际金融危机和国内市场影响，龙口市固定资产投资、地区生产总值、财政总收入等均增长较小，导致这2年土地经济可持续水平增长不明显。随着经济发展，建设用地占用增加，并且建设用地增长快于经济增长，尤其2004、2011年经济用地弹性系数分别为14.27、2.51，建设用地的消耗与经济的产出不协调。2011年后情况有所改善。

5.3 生态环境可持续水平

龙口市土地生态环境可持续水平总体呈上升趋势，但期间有较大波动。2004年生态发展可持续分值为26.05分， 2005年有所增加，达到33.77分， 2005—2008年逐渐下降， 2008年达到相对最低值。2009—2012年呈快速上升趋势，年增长率达20%， 2012年达到相对最高值。2013年有所回落， 2014年开始又逐渐上升， 2015年土地生态环境可持续分值为97.23分。

2004—2015年，龙口市城镇生活污水处理率、林地覆盖率总体呈提高趋势， 2015年分别达到100%、55.60%，工业废气排放达标率除2010、2011年为99.70%、99.00%，其余各年均达到100%，这为土地生态环境可持续利用提供了基本保障。但2006年，城镇生活污水处理率较2005年有较大下降，由87.89%下降为80.99%， 2005—2008年影响生态环境的负向指标单位耕地化肥增加较快， 2008年达到最大值2771.27kg/hm2，对生态环境带来较大影响，造成土地生态可持续利用水平降低。2013年土地生态可持续发展水平下降同样与单位耕地农药负荷有关系，可见其对生态环境的影响比较突出。

5.4 综合可持续利用水平

龙口市土地综合可持续利用水平总体呈上升趋势，土地资源可持续利用状况不断改善。龙口市土地可持续利用水平大致分为两个阶段：第一阶段(2004—2009年)，稳定发展阶段，由图1可知， 2004—2009年土地可持续利用水平稳定增长，年均增长率为6.07%， 2009年土地可持续利用水平达到36.40分； 第二阶段(2009—2015年)快速上升阶段，其中2009—2012年综合得分由36.40上升到83.78，上升速度幅度较快， 2012年后上升趋势有所减缓。

2004年以来，龙口市人们社会生活水平不断提高，经济发展成果显著，生态保护意识逐步提高，不断加强生态环境保护建设，从而提高了土地可持续利用整体水平。但龙口市还存在人口增长与城乡建设用地增长不协调、建设用地消耗与经济产出不协调、单位耕地农药化肥负担较高等问题。今后需在城乡建设用地增长、耕地农药化肥施用等加强管制，从而不断提高土地可持续利用水平，促进社会、经济与生态协调发展。

6 结论与建议

6.1 结论

(1)该文以山东省龙口市为研究区，创新引入建设用地弹性系数，构建包括社会发展、经济发展和生态环境3个因素、19个因子的评价指标体系，基于主成分分析法(PCA)，客观评价2004—2015年龙口市土地可持续利用状况与变化。

(2)2004年以来龙口市土地利用可持续发展综合水平逐年提高，其中社会发展、经济发展可持续水平持续上升，生态环境可持续水平整体呈上升趋势但波动较大。到2015年，龙口市土地社会、经济可持续利用水平达到相对最高值，发展态势良好，生态环境可持续利用水平有待提高。

(3)人口变化与城乡建设用地的增加不协调、经济增长与建设用地的消耗不协调以及耕地环境污染可能成为目前龙口市土地可持续利用的制约因素。

6.2 建议

(1)强化“一户一宅”制度，加强农村宅基地登记管理，开展村庄用地整治及城乡建设用地增减挂钩，对闲置宅基地、闲置房屋进行整治，加大村庄建设用地的“减量”力度，提高村庄用地集约度。

(2)继续加大单位面积土地的投入程度，提高土地的产出效率，适应经济发展新常态，转变经济发展方式，全面优化产业结构，加快经济的发展速度。

(3)加强生态环境保护建设，倡导发展生态农业，合理进行化肥、农药、塑料薄膜的使用，提高化肥农药利用率，鼓励农民增施有机肥，减少对耕地资源的污染，进一步提高生态环境可持续性。