适宜性分析是土地资源对某种特殊利用适合程度的确定过程，是生态规划的重要分析手段[1]。随着计算机的普及和GIS技术的发展，适宜性分析已逐渐由人工方式转变为计算机技术支撑，影响因素也由单因子发展到多因子复合因素的综合考虑，适宜性分析先后应用于农场规模、资源分配[2]和生态学的多准则决策[3]。21世纪后，随着计算机专家系统的建设和完善，地理信息专家系统开始应用于土地适宜性分析[4]。本文借鉴了近年来国内的相关研究，如单勇兵[5]等从自然、环境、人文3个方面确定了10个影响因子，并通过取最小值法对徐州市乡村聚落布局进行了适宜性评价；唐常春[6]等遵循“三步骤四原则”，以行政区划为评价单元对长江流域进行了空间开发适宜性分析；祁豫玮[7]等以开发引导和约束为标准建立了二维判别矩阵，并将其与开发强度进行组合判断，从而进一步研究南京市空间开发功能类型。

1 研究区概况与数据来源

1.1 研究区概况

新沂市隶属江苏省徐州市，地势平坦，降水季节变化明显，区域内土壤和水资源总量充足但污染严重，如图1所示。

新沂市地处苏北地区重要的发展轴上，为顺应“一带一路”发展，政府正加紧实现经济发展的战略对接；但因地质活动不稳定导致其对大规模经济和人口的聚集能力不高；受徐连双核结构[8]的影响，其经济发展相对落后。近年来，旅游业开发是新沂市GDP增长的主要部分[9]。另外，城镇等级规模体系的缺陷性使得城镇结构布局不够紧凑[10]，城镇整体功能及其各支撑要素间的联通性被降低；生态环境的破坏和发展的不平衡日益明显，导致其资源环境承载力大大降低且空间分布不均匀。

本工程地下室顶板覆土大部分厚1.5m，标准跨度为8.4m×8.4m，综合考虑经济性且避免以往类似的地下车库工程发生的事故，选择井字梁作为地下室顶板的承重结构，其他特殊区域如坡道、楼梯、设备间和有较大集中荷载的区域按实际情况设置主次梁。经计算，将主次梁的配筋率控制在1.5%～2%，且地下1层梁顶面和底面的纵向钢筋应比计算值增大10%。

图1 新沂市地理位置（审图号：苏S(2016)096）

1.2 数据来源与预处理

本文采用的主要数据包括新沂市地形图、新沂市土地利用现状图、土地资源图、土壤图以及其他资料（交通网络、统计年鉴等）。

在完成评价指标指数计算后，采用互斥矩阵法进行适宜性等级划分。综合考虑资源环境承载力和空间开发潜力，根据适宜性指数建立互斥矩阵（表2）。

2 适宜性分析过程与方法

本文以Geodatabase为基础划分评价单元，并进行空间分析。根据适宜性分区评价的需要，本文首先选取生态重要性、灾害易损性、水环境容量、可建设用地比例、交通可达性以及人口与经济集聚度6个评价指标；再采用Delphi法进行权重计算，并利用互斥矩阵法实现矩阵分类；最后对区域空间开发适宜性进行分析和制图，并结合当地自然经济社会等情况，对比现今的空间开发情况，给出政策性建议，以促进区域生态、经济等的可持续发展。具体流程如图2所示。

图2 适宜性分析流程图

2.1 评价单元的选择

评价单元的选择方法一般包括自然边界法、行政边界法和网格法。自然边界法和行政边界法在结果的完整性和辨识度方面存在一定缺陷；网格法则是根据计算精度确定一定大小的网格单元作为评价单元，评价精度和工作量与网格单元大小成正比。虽然网格法划分评价单元会割裂行政区划和自然边界，但其划分方式灵活简便，单元大小具有一定的可控性。考虑到研究区属于县级市，行政分级较小，且需要打破行政区界进行统筹分析，本文采用网格法划分评价单元。

2.2 指标选择与权重判断

根据一系列指标图层属性值和适宜性指标指数的获取与计算，利用ArcGIS制图分析功能，在ArcMap中进行专题图的制作与输出。本文充分考虑6个指标在研究区域各自的分布情况，尤其是重点影响因子在研究区的分布，以专题图的形式表示[14]；利用ArcGIS空间分析功能，结合由Delphi法计算得到的各单因子权重，判断各单因子影响的重要程度，从而计算空间开发适宜性指数（图3）。其计算公式为：

空间开发潜力包括可建设用地比例、交通可达性和人口与经济集聚度3个指标。可建设用地比例反映了区域开发建设的饱和程度，比例越大，潜力越高。城镇作为交通网络上的节点，交通可达性反映了其交通区位的优劣[12]。新沂市交通区位优越，交通运输主要以汽运、铁路为主，且正处于高铁和水运的发展期，但市内公共交通设施有待进一步完善。人口与经济集聚度分别是指单位面积上居住人口的总量和单位面积上国民生产总值[13]，通过综合人口密度、城镇化体系结构、人均GDP以及建设用地二三产业产出研究区域发展情况。

当然这一切，我都没有告诉黄玲。我们时常一起上班，一起下班走那段昏暗的小巷子。旁边的大楼都已经竣工，这使我想起那次被抢劫的遭遇，对黄玲说，其实如果那天你不出现，我就撒手了，他们只是两个月没领到工资的民工，不是什么坏人。

2.2.2 空间开发潜力

协方差矩阵特征值的大小反映了其对应的特征变量所包含的信息多少。将协方差矩阵的特征值降序排列，取前k个特征值对应的特征向量(kp)，这些特征向量称为主成分，这个过程称为降维[13]。根据特征值大于1的规则，以此来确定k的数值。PCA在分解矩阵时不单独考虑输出矩阵，因此用该方法提取主成分获取的信息比PLS更充分。PCA将原始矩阵X分解为式(9)。

确定权重的方法主要包括Delphi法、两两比较法和熵值确定法。考虑到通用性、计算过程和工作量等因素，本文采用Delphi法确定权重。首先客观地选择数名专家，请他们根据相关资料和指标独立地进行权重判断；然后收集所有判断结果并计算各指标所赋权重的均值和标准差；再将计算结果返还给各位专家，重复上述过程，直到各专家给出的权重判断结果具有较强的一致性为止。利用Delphi法确定评价区域的各等级评价指标权重，构建评价指标如表1所示。

2.2.1 资源环境承载力

2.2.3 权重判断

在设备频率确定后，正确选择比功率值，对满足零件的技术要求是十分重要的。如上所述，频率、比功率确定后，加热速度也就随之而定，经过某一时间后，零件表面一定的厚层就被加热到淬火温度。此时进行淬火冷却，加热淬火层的深度也就完全确定了。实际上，在透入式加热条件下因表面升温速度相对较小，大都是用调整时间的方法来满足要求的加热深度。过多地延长加热时间，将引起表面的严重过热。因此，此时就要重新选择较低的比功率。相反，如要求加热层深度很小时，就需要选用较大的比功率进行加热。

The Current Situation of Agricultural Products Logistics and the Conception of Logistics Systemin Henan Province: based on the Case of Zhengzhou

资源环境承载力包括生态重要性、灾害易损性和水环境容量3个指标。生态越敏感，损后恢复能力就越差，越应重视对其的生态保护。灾害易损性是指某区域易受某灾害影响的客观程度以及该区域人们应对灾害的主观反应能力。新沂市地处郯庐断裂带之上，灾害易损性较高。水环境容量反映了水体对污染的最大承受能力和损后自净能力，其与污染承受力和自净力成正比。新沂市水域面积广阔，可调用水量达14×108 m3，境内骆马湖水域水质可达国家二级标准[11]；水利用方式多为农业灌溉、生活和工业用水，化工业等重金属污染为水体污染日渐严重的主要因素。

2.3 适宜性等级划分

在ArcGIS系统的支撑下，对各图形要素进行数字化处理并配准；运用ArcGIS中Spatial Analyst模块进行距离、密度和坡度分析，获取评价单元的各类指标值；并利用Delphi法进行原始数据标准化处理，计算适宜性指数并分级。

根据资源环境承载力指数大小，将研究区资源环境承载力分为高、较高、中等、较低和低5个等级；根据空间开发潜力指数大小，将研究区空间开发潜力也划分为高、较高、中等、较低和低5个等级；各等级彼此组合，形成25个矩阵单元，并将其归到高、较高、中等、较低和低5类适宜性程度中。资源环境承载力低的区域，无论其空间开发潜力高低，均为低空间开发适宜性；空间开发潜力低的区域，当其资源环境承载力为中等偏低时，其空间开发适宜性较低；资源环境承载力高且空间开发潜力大的区域，其空间开发适宜性也较高；资源环境承载力和空间开发潜力都属于中等偏下水平的区域，其空间开发适宜性较低。

3 研究结果分析

3.1 适宜性分区

从研究区实际出发，在确保分析精度的基础上，根据指标选择的可获性、特殊性和综合性原则，结合分析目的等相关要求，本文从资源环境承载力和空间开发潜力两方面着手，选取生态重要性、灾害易损性、水环境容量、可建设用地比例、交通可达性以及人口与经济集聚度6个指标进行综合评价。

表1 空间开发适宜性分析指标及其权重分配

系统层 因子层 权重 指标层 权重生态重要性 0.200公园 0.020风景名胜区 0.030重要山体丘陵 0.050水源涵养区 0.050饮用水重点保护区 0.050资源环境承载力 灾害易损性 0.150旱涝灾害 0.010地震灾害 0.050地面塌陷 0.040滑坡 0.012泥石流 0.012地裂灾害 0.026水环境容量 0.200水体特征 0.040水质目标 0.025水体通达性 0.050污染物形式 0.050水利用方式 0.035可建设用地比例 0.150城市居住用地 0.030工农建设用地 0.037公共设施用地 0.033军事用地 0.023旅游用地 0.027空间开发潜力 交通可达性 0.150运输成本 0.030道路服务水平 0.027道路等级 0.020距离阻抗 0.043中心区可达性 0.030人口与经济集聚度 0.150人口密度 0.035城镇体系结构 0.028人均GDP 0.052建设用地二三产业产出 0.035

表2 互斥矩阵

资源环境承载力高 较高 中等 较低 低空间开发潜力高 高适宜 较高适宜 较高适宜 中等适宜 低适宜较高 高适宜 较高适宜 中等适宜 中等适宜 低适宜中等 较高适宜 中等适宜 中等适宜 较低适宜 低适宜较低 较低适宜 较低适宜 较低适宜 较低适宜 低适宜低 较低适宜 较低适宜 低适宜 低适宜 低适宜

式中，sj为第j个单元的适宜性指数；wi为第i个指标的权重；xij为第j个单元第i个指标的指数。

从资源环境承载力和空间开发潜力两个方面进行互斥矩阵分类，得到适宜性分区并制作专题图（图4）。由图4可知，高适宜空间开发区主要分布在中部偏北，即新安镇主城区，因为主城区经济发展速度快，空间开发潜力大，人口与经济密集度高，且资源环境承载力较好；受市区辐射影响，北沟镇、唐店镇与市区交界地区空间开发适宜性也较高。就各乡镇而言，由镇中心区开始，呈同心圆状向四周扩散，适宜性程度递减。其中，窑湾、合沟等边缘乡镇受邳州市影响，加之良好的自然条件和丰厚的历史文化底蕴，存在部分中等以上适宜性的区块。然而，由于交通可达性低、乡镇各支撑要素联通困难，农业占比大、可持续意识不高，生态资源的自我再生速度降低、资源环境承载力低，导致其他乡镇空间开发适宜性较低。另外，骆马湖、沭河沿岸以及马陵山地区均为低适宜空间开发区，归因于长期以来的过饱和开发，开发潜力较低，且生态污染情况日益严重。

图3 空间开发适宜性指数分布图

图4 空间开发适宜性分类结果

本文选用的网格大小为200 m×200 m，根据各适宜性等级所占网格数量，确定各等级所占面积的数量指标（表3），可以看出，新沂市中等以上适宜性空间开发区域面积占比较少，约为7 %；其中高适宜空间开发区仅占2 %；较低适宜空间开发区占比最大，约占53 %。

表3 各等级面积统计

适宜性等级 面积/km2 面积占比/%低适宜性 300.84 19较低适宜性 842.96 53中等适宜性 324 21较高适宜性 85.72 5高适宜性 24.92 2

3.2 空间开发区划建议

本文将空间开发适宜性分析结果与区域空间开发强度现状相结合[15]，建立二维判别矩阵（表4），划分为5类空间开发区，分别是立体提升区、经济拓展区、适度开发区、后备发展区和生态养护区；并对研究区进行空间开发区划建议，制作专题图进行开发分区（图5）。

表4 空间开发矩阵法分区

空间开发适宜性高适宜 较高适宜 中等适宜 较低适宜 低适宜空间开发强度现状高 经济拓展区 经济拓展区 后备发展区 后备发展区 生态养护区较高 经济拓展区 经济拓展区 适度开发区 后备发展区 生态养护区中等 立体提升区 经济拓展区 适度开发区 后备发展区 生态养护区较低 立体提升区 经济拓展区 适度开发区 后备发展区 生态养护区低 立体提升区 立体提升区 适度开发区 适度开发区 生态养护区

图5 空间开发区划

对于低适宜空间开发区，无论空间开发强度如何，均划分为生态养护区；对于适宜性高且未充分开发的地区，则划为立体提升区。由图5可知，立体提升区分布在新安镇主城区、北沟和草桥镇，其中主城区可利用交通区位条件、基础设施条件和资金技术等优势，发展精细化工、食品加工、纺织、机械等第二产业，促进市区经济整体多方面提升，完善商旅服务功能；考虑到窑湾、棋盘和高流等良好的生态区位条件和较适宜的空间开发程度，可进行经济拓展开发，在骆马湖等河湖沿岸发展水稻和水产品生产，但同时也应注意生态资源可持续，以旅游业代替高污染工业；而双塘、新店和邵店镇的生态承载力较低，应重视生态养护，可选择性的发展物流业、生态农业等第三产业。

4 结 语

本文在充分认识土地利用适宜性评价重要性的前提下，综合考虑多方面影响因子，选取了6个评价指标；利用Delphi法判断权重，计算适宜性指数；利用互斥矩阵法划分适宜性等级并制图分析，从而将新沂市空间开发划分为高适宜性、较高适宜性、中等适宜性、较低适宜性和低适宜性5类。结合其空间开发现状，本文提出了空间开发区划和发展方式建议，开发强度由高到低分别为立体提升区、经济拓展区、适度开发区、后备发展区和生态养护区。由于本文通过网格法划分评价单元，网格大小和不同网格的邻近性导致评价指标的选择具有一定局限性，且选择过程中忽视了地域基本农田政策的部分影响；同时，Delphi法确定的专家人数有限，通过专家进行权重判断存在一定的主观性，因此如何更为合理地进行评价指标选择和权重判断，进一步提高适宜性分析的准确性有待下一步研究。

为加大产业扶贫力度，《指导意见》提出，多渠道拓宽农产品营销渠道，推动批发市场、电商企业、大型超市等市场主体与贫困村建立长期稳定的产销关系，支持供销、邮政及各类企业把服务网点延伸到贫困村，推广以购代捐的扶贫模式。加快推进“快递下乡”工程，完善贫困地区农村物流配送体系。

参考文献

[1] 许捍卫,马文波,赵相伟,等.地理信息系统教程(2010)[M].北京:国防工业出版社,2010:229

[2] Campbell J C, Radke J, Gless J T, et al. An Application of Linear Programming and Geographic Information Systems:Cropland Allocation in Antigua[J]. Environment and Planning A,1992(24):535-549

[3] Pereira J C, Lucien D. A Multiple Criteria Decision-making Approach to GIS-based Land Suitability Evaluation[J]. International Journal of Geographical Information Systems,1993,7(5):407-424

[4] 史冲.基于GIS的土地适宜性评价研究[J].北京测绘,2013,9(5):38-40

[5] 单勇兵,马晓冬,宣勇.基于GIS 的徐州市乡村聚落空间适宜性分析[J].地域研究与开发,2012,31(6):156-160

[6] 唐常春,孙威.长江流域国土空间开发适宜性综合评价[J].地理学报,2012,67(12):1 587-1 598

[7] 祁豫玮,顾朝林.市域开发空间区划方法与应用:以南京市为例[J].地理研究,2010,29(11):2 035-2 044

[8] 陆玉麒.区域双核结构模式的形成机理[J].地理学报,2002,57(1):85-95

[9] 江苏省统计局.江苏统计年鉴[EB/OL].(2015-02-21)[2016-06-30].http://www.jssb.gov.cn

[10] 马随随.主体功能区导向下新沂市空间结构现状及其优化对策[J].淮海文江,2013,13(3):33-35

[11] 胡晓萌.新沂市乡镇空间组织优化研究[D].徐州:江苏师范大学,2012

[12] 王先鹏,戚红年,曹荣林.基于GIS的县域交通网络可达性探讨[J].地理空间信息,2009,7(5):109-113

[13] 闫庆武,马晓冬,卞正富.基于地统计的徐州市人口密度空间分布研究[J].云南地理环境研究,2007,19(4):13-18

[14] 金志丰,陈雯,孙伟,等.基于土地开发适宜性分区的土地空间配置:以宿迁市区为例[J].中国土地科学,2008,22(9):43-50

[15] 陆树启,陆玉麒,侯立春.影响县域发展的空间结构:以江苏省新沂市为例[J].经济地理,2006,26(6):919-921