京津冀城市群国土空间适宜性格局特征及分区研究.pdf

第38 卷第4期2023年10 月文章编号：2 0 96-47 49(2 0 2 3)0 4-0 46 9-12引文格式：徐磊，张世珍，胡浩，杨静园，张长春，马立军.京津冀城市群国土空间适宜性格局特征及分区研究.林业与生态科学,2 0 2 3,38 4)：46 9-48 0XU Lei,ZHANG Shizhen,HU Hao,YANG Jingyuan,ZHANG Changchun,MA Lijun.Regional characteristics andregionalization of land space suitability in Beijing-Tianjin-Hebei urban agglomerationJJ.Forestry and EcologicalSciences,2023,38(4):469-480.林业与生态科学FORESTRY ANDECOLOGICAL SCIENCESDOI:10.13320/ki.hjfor.2023.0058Vol.38 No.4Oct.2 0 2 3京津冀城市群国土空间适宜性格局特征及分区研究徐磊1.2，张世珍1，胡浩1，杨静园，张长春1.2,马立军1（1河北农业大学国土资源学院，河北保定0 7 10 0 1;2 河北农业大学乡土文化与乡村治理研究中心，河北保定0 7 10 0 1；3河北农业大学资源与环境科学学院，河北保定0 7 10 0 1)摘要：科学实现空间适宜性分区是国土空间优化的基础。为了优化国土空间适宜性格局，基于网格尺度构建了生态保护功能一农业生产功能一城镇建设功能的城市群国土空间适宜性评价指标体系，采用多因素综合评价法和标准显示性比较优势指数模型（NRCA),得到京津冀城市群国土空间适宜性格局特征，给出了6 种功能类型区的分区优化建议。研究表明：（1)京津冀城市群国土空间适宜性具有明显的空间异质性特征。生态保护功能低适宜性和一般适宜地区集中于京津冀城市群西北部和东南部的平原地区；农业生产功能低适宜性地区主要分布在坝上高原、中部平原地区；燕山、太行山山脉以北城镇建设功能适宜性较低。（2)划分京津冀城市群国土空间为生态保护功能优势区、农业生产功能优势区、城镇建设功能优势区、生态保护一农业生产功能优势区、生态保护一城镇建设功能优势区、农业生产一城镇建设功能优势区，并提出各类型功能区国土空间开发利用与优化的策略。关键词：国土空间；开发适宜性；比较优势；分区优化；京津冀城市群中图分类号：F205文献标志码：ARegional characteristics and regionalization of land space suitabilityin Beijing-Tianjin-Hebei urban agglomerationXU Leil2,ZHANG Shizhen,HU Hao,YANG Jingyuan,ZHANG Changchun2,MA Lijun(1 College of Land and Resources,Hebei Agricultural University,Baoding O7100l,China;2 Local Culture and Rural Governance Research Center,Hebei Agricultural University,Baoding 07100l,China;3 College of Resources and Environmental Sciences,Hebei Agricultural University,Baoding 07100l,China)Abstract:Scientific space suitability zoning was the basis of territorial space optimization.Inorder to optimize the spatial suitability of urban agglomeration,an evaluation index system收稿日期：2 0 2 3-0 4-10基金项目：2 0 2 2 年度河北省社会科学发展研究课题（2 0 2 2 0 30 30 8 7）；河北省高等学校人文社会科学研究项目（SQ2023207）；河北省省属高等学校基本科研业务费研究项目（KY2022100）；河北省高等学校人文社会科学重点研究基地经费资助。第一作者：徐磊（198 5-)，男，河北唐山人，讲师，研究方向为土地资源利用与区域经济发展。通信作者：马立军（198 0-），男，河北承德人，副教授，研究方向为区域经济与土地管理。470of ecological protection function,agricultural production function and urban constructionfunction was constructed based on grid scale.Multi-factor comprehensive evaluation methodand standard display comparative advantage index model(NRCA)were used to obtain thespatial suitability characteristics of the Beijing-Tianjin-Hebei urban agglomeration.The opti-mization suggestions of six functional areas were given.The results showed that:(1)Thespatial suitability of the Beijing-Tianjin-Hebei urban agglomeration had obvious spatial heter-ogeneity.The areas with low suitability and general suitability for ecological protection wereconcentrated in the plain areas in northwest and southeast of the Beijing-Tianjin-Hebei urbanagglomeration.The areas with low suitability of agricultural production function were main-ly distributed in Bashang plateau and central plain.The suitability of urban constructionfunction in the north of Yanshan and Taihang Mountains was low.(2)The territorial spaceof the Beijing-Tianjin-Hebei urban agglomeration was divided into ecological protection func-tional advantages,agricultural production functional advantages,urban construction func-tional advantages,ecological protection-agricultural production functional advantages,eco-logical protection-urban construction functional advantages,agricultural production-urbanconstruction functional advantages,and put forward the development and utilization and op-timization strategies of each type of functional areas of the territorial space.Keywords:territorial space;development suitability;comparative advantage;partition optimi-zation;Beijing-Tianjin-Hebei urban agglomeration我国人口数目庞大，土地资源有限，保障农业生产空间是实现粮食安全的重要物质基础 1。同时，近年来从中央到地方都把生态文明建设纳人重要议程，着力加快生态空间绿色发展。由于国土空间规划不合理和城市空间扩张管理不当，城镇规模迅速扩张，建设用地不断扩张，耕地锐减，生态问题频发，国土空间规划呕待完善。揭示区域国土空间适宜性的格局特征、科学优化国土空间是解决人地矛盾、实现生态宜居的重要路径，也是保障我国社会经济可持续发展的重要手段 2。京津冀城市群是华北地区经济规模最大的城市群，识别国土空间适宜性分布规律，优化国土空间开发格局是实现京津冀地区协同发展的重要举措，是带动北方经济社会深化发展的必经之路。国内外学者关于国土空间适宜性相关研究较早，研究视角也不断创新，从单一指标评价到多维度发展，单一指标维度有从生境适宜性评价、人居环境适宜性评价和土地开发利用适宜性评价等方面开展研究，多指标维度有从“城镇一农业一生态”空间和从“三生”空间视角进行国土空间适宜性评价 3-7 。基于网格单元尺度，以“农业一生态一建设”功能的空间载体作为评价对象，进行量化识别国土空间适林业与生态科学宜性并优化空间布局的研究还处于薄弱环节。在研究内容上，囊括了理论内涵、指标体系构建、表征和测算、三生空间演变预测、国土空间分区优化等方面 8-12 。在研究区域上，多从市域、县域等范围选取,随着研究深入，关于省域、海岛等特殊地形区域尺度的研究分析也逐步增多 13-1。在研究方法上，运用较多的有层次分析法、专家打分法和GIS空间分析等，近年来，贝叶斯网络、物元模型、综合评价等模型构建也日益常见 15-18 。此外，国内外学者基于生态位、系统科学等理论，从县域、市域、省域、城市群和特殊地貌山区等研究尺度入手，运用指标体系法、三维魔方法和比较优势指数法等方法进行了分区优化，并提出相应对策与建议 9-10,12,19-2 4。国土空间开发适宜性评价作为国土空间规划的重要性基础工作，已有研究多从落实国土空间规划的市、县尺度人手，评价“三生”空间适宜性，对城市群进行适宜性评价并分区优化的研究还较为少见。特别是基于网格单元，从生态、农业、建设功能空间分类的角度，量化识别国土空间适宜性并优化城市群国土空间布局的研究更为缺乏。综上所述，以京津冀城市群为研究区域，以生态保护功能、农业生产功能和城镇建设功能为视角，采取多因素综合评价第38 卷第4期法综合评价国土空间适宜性水平，利用ArcGIS软件空间聚类分析揭示京津冀城市群国土空间适宜性的空间集聚特征，运用NRCA模型提出国土空间开发适宜性分区方案，为区域协调可持续发展提供参考。1研究区概况及数据来源1.1研究区概况京津冀城市群被称为中国的“首都经济圈”，包括北京、天津两大直辖市和河北省的石家庄市、保定市、唐山市、张家口市、承德市、秦皇岛市、廊坊市、沧州市、衡水市、邢台市和邯郸市共11个地级市。地理位置为E1132711950,N360542 40 ，城市群总面积为2 1.8 万km。京津冀城市群地处华北平原，境内河流水系众多，生物矿产资源丰富，地貌类型多样，涵盖高原、山地、丘陵、盆地、平原和洼地，地势整体上呈由西北向东南倾斜，山地约占区域总面积的一半。2 0 2 1年京津冀城市群国内生产总值约96 355.95亿元，占全国的8.4%，是中国北方经济规模最大的地区。1.2数据来源数据来源主要包括：2 0 18 年京津冀土地利用现状数据，空间分辨率为30 m，以Land TM和Land-sat-8遥感数据为主要信息源，利用ENVI与Arc-GIS软件解译所得;数字高程数据(DEM)源于地理空间数据云平台（http：/w w w.g s c l o u d.c n/）；植被覆盖度基础数据为MODIS13Q1，源于陆地产品处理分发数据中心（https:/lpdaac.usgs.gov)；土壤属性数据采用寒区旱区科学数据中心（http：/）；京津冀道路数据由高德地图提供；水系数据由国家青藏高原科学数据中心()提供；人口数据的原始数据来自2019年中国县域统计年鉴和各省市的统计年鉴。本研究考虑到研究区域范围、数据类型特征、可计算性及空间分辨率等因素，所有数据统一转换为高斯一克里格2 0 0 0 国家大地坐标系，通过栅格重采样转换为30 0 0 mX3000m栅格进行计算，研究区共分为2 39 9 1个评价单元。徐磊，等：京津冀城市群国土空间适宜性格局特征及分区研究W,=0.5XM+0.5XN;式中：W,为指标i的权重；M,为层次分析法测算指标权重；N,为独立性权重系数法测算指标权重。2.1.2国土空间适宜性评价指标选取主要包括生态保护功能适宜性、农业生产功能适宜性、城镇建设功能适宜性3方面。(1)生态保护功能适宜性。选取生境质量指数、生态系统服务价值、植被覆盖度、水系密度、生态结构指数和土壤可蚀性表征国土空间生态保护功能适宜性。生态系统服务价值参照谢高地单位面积生态系统服务价值当量表测算 2 6 。水系密度为评价单元内的水系总长度与网格面积的比值，能够反映区域水资源的短缺程度，以水系矢量数据为基础，基于ArcGIS的标识和统计工具测算。土壤可蚀性参考生态保护红线生态功能评价技术指南的土壤可蚀性因子K测算。生态结构指数能够反映区域的防风固沙和水土保持情况，参照生态环境状况评价计算规范测算，见公式2。ESI=C,X0.25+CgX0.25+R,X0.5式中：ESI为生态结构指数；Cf、Cg 分别表示林草地覆盖率、草地覆盖率；R，表示水域湿地面积占比。参照现有研究，综合考虑京津冀城市群研究区概况和专家意见，本研究以旱地、城镇用地、工业用地、公路、铁路、农村居民点用地和未利用地作为威胁因子 2 7 。关于威胁因子参数设置见表1。4712研究方法2.1国土空间适宜性评价体系构建2.1.1国土空间适宜性评价指标体系石的评价单元较多，确定权重应综合考虑主观层面指标的重要程度以及数据的客观规律，采用层次分析法和独立性权重系数法实现主客观相结合计算权重。层次分析法将目标划分为目标层、准则层和方案层，基于指标的重要性确定专家打分表，从生态保护功能、农业生产功能和城镇建设功能3方面构建判断矩阵以确定权重。独立性权重系数法是依据指标间相关程度确定权重的客观赋值法，基于回归分析确定指标的共线性强弱 2 5。综合权重计算公式为：(1)(2)研究涉及472威胁因子Threat factor旱地城镇用地农村居民点用地未利用地公路铁路工业用地生境质量是生态环境对人类的适宜程度，根据已有研究，土地利用程度越强对生境的威胁越大 2 8 。以土地利用现状数据为基础，设定各土地利用类型的土地利用类型生境适宜度旱地Land-useHabitattypesuitability水田0.38旱地0.37有林地0.80灌木林0.86疏林地0.80其它林地0.00高覆盖度草地0.73中覆盖度草地0.68低覆盖度草地0.55河渠0.84湖泊1.00水库坑塘0.88滩涂0.83滩地0.33城镇用地0.00农村居民点0.00其它建设用地0.00沙地0.01盐碱地0.01沼泽地0.60裸土地0.01裸岩石质地0.24林业与生态科学表1威胁因子表Table 1Table of threat factors权重Weight0.70.70.70.30.60.60.9表2 土地利用类型对威胁因子敏感程度Table 2Sensitivity of land use types to threat factors城镇用地IndustrialDry landUrban land0.400.530.100.620.530.640.430.680.500.770.000.000.460.560.520.500.460.530.450.740.650.850.530.820.320.950.030.950.000.000.000.000.000.250.100.100.100.100.470.400.050.100.030.03第38 卷最大威胁距离空间衰退类型Maximum threat distanceSpatial decay type5指数5指数5指数5指数5线性5线性5指数生境适宜度和土地利用类型对威胁因子敏感程度，基于InVEST模型测算生境质量指数。生境适宜度参数设置和土地利用类型对威胁因子敏感程度见表2。工业用地未利用地公路铁路HighroadRailwayland0.530.560.700.660.750.000.540.550.590.800.850.700.870.630.400.750.330.000.000.470.000.00农村居民点用地Ruralresidential land0.330.300.300.250.430.430.400.470.680.780.000.000.300.350.380.480.300.370.570.600.550.600.580.620.700.700.270.270.000.000.000.000.000.000.200.200.200.200.470.470.100.100.100.10Unusedland0.550.000.540.450.600.600.550.600.700.000.000.000.460.700.430.000.450.000.750.500.720.000.750.700.870.750.300.400.000.000.000.000.050.000.100.000.100.000.500.000.050.000.030.00第4期(2)农业生产功能适宜性。采用坡度和土壤的理化性质测算农业生产功能适宜性。土壤理化性质包含土壤pH值、土壤剖面深度、土壤容重、土壤孔隙度和土壤有机碳，评价单元土壤理化性质分值采用分区统计测算获得。坡度以DEM高程为原始数据，在ArcGIS基于坡度处理获得。(3)城镇建设功能适宜性。路网密度为评价单元内全部道路长度与评价单元面积的比值。交通干线可达性为评价单元中心点与1到4级非高速公路的距离。采用京津冀地区乡镇人口及乡镇行政区面积测算人口密度，以各乡镇点人口密度为中心进行Table 3 Weight and classification of evaluation indexes of suitability of ecological protection function in territorial space国土空间功能Territorial spacefunction指标指标名称IndexIndex name指标层次分析法权重独立性权重系数法综合权重指标100分级80赋值604020依据层次分析法和独立性权重系数法主客观相结合计算所得权重，参照已有研究，结合自然间断点徐磊，等：京津冀城市群国土空间适宜性格局特征及分区研究表3国土空间生态保护功能适宜性评价指标权重及指标分级Ecological protection function生境质量生态系统指数服务价值/元VegetationHabitatEcosystemquality indexservicevalue0.060.060.040.040.050.050.672.360.530.671.002.360.460.530.391.000.310.460.120.390.310.801.470.710.800.671.470.570.710.390.670.350.570.290.390.350.230.130.150.230.130.160.030.150.160.190.010.030.190.210.21表4国土空间农业生产功能适宜性评价指标权重及指标分级Table 4 Weight and classification of evaluation indexes of suitability of agricultural production function in territorial space国土空间功能Territorial space function指标指标名称IndexIndex name指标层次分析法权重独立性权重系数法综合权重农业生产功能Agricultural production function土壤剖面平均坡度/土壤pHAverageSoil pHslope0.050.050.050.040.050.05土壤容重/深度/cm(kg m-3)Soil profileSoil bulkdepthdensity0.050.050.060.060.050.05土壤有机碳/土壤孔(g kg-2)度/%Soil organicSoil porositycarbon0.050.090.050.060.050.08474国土空间功能Territorial spacefunction指标指标名称IndexIndex name指标100分级80赋值604020国土空间城镇建设功能适宜性评价指标权重及表5国土空间城镇建设功能适宜性评价指标权重及指标分级Table 5Weight and classification of evaluation indexes of functional suitability of urban construction in territorial space国土空间功能Territorial spacefunction指标IndexIndex name指标层次分析法权重独立性权重系数法综合权重指标100分级80赋值6040202.2国土空间适宜性评价计算模型多因素综合评价法能够综合反映多个指标的适宜性分值，采用多因素综合评价法分别测算生态保护功能适宜性、农业生产功能适宜性和城镇建设功能适宜性分值，依据自然间断点分级法和分位数法将国土空间适宜性划分为低适宜、一般适宜、中等适宜和高适宜4类 18 18F,=ZW:X G1式中：F，为第个评价单元的生态保护功能适宜性、农业生产功能适宜或城镇建设功能适宜性分值；G，为第个评价单元下第i项指标的分级赋值得分。林业与生态科学表4(续）农业生产功能Agricultural production function土壤剖面平均坡度/土壤pHAverageSoil pHslope6.57.535.56.5387.58.58158.525距离乡镇地形起距离/m指标名称伏度/Topographicrelief0.050.040.0525第38 卷土壤容重/土壤有机碳/土壤孔深度/cm(kg:m-3)Soil profileSoil bulkdepthdensity116.281.26106.90116.281.261.28100.05106.901.281.3184.05100.051.311.341.34指标分级，见表5。城镇建设功能Urban construction function交通干线路网密度/可达性Distance(km:km=2)from townsRoad densityand villages0.050.070.066500(3)(g:kg-2)隙度/%Soil organicSoil porositycarbon52.121.3150.7352.120.891.3149.9550.730.720.8949.4549.950.610.7249.453.568001.443.5680020000.671.44200035000.100.6735006.50065002.3国土空间适宜性冷热点分析冷热点分析是通过高值和低值在空间上的聚类，采用ArcGIS软件冷热点分析工具，揭示国土空间适宜性的空间特征 2 3。通过检验统计量得分、概率P值和置信区间识别国土空间适宜性的热点和冷点区域。得分表示标准差的倍数，P值即为概率,对国土空间适宜性显著性进行评估，P值越小，表明结果越显著；显著冷热点区域、冷热点区域和次冷热点区域的置信度分别为9 9%、9 5%和90%。2.4比较优势功能指数比较优势指数不受时空限制，已应用于资源利用、贸易往来、功能评价等多个领域 2 0.2 4。已有研究测度方法中标准显示性比较优势指数(NRCA)不Surfaceroughness0.060.060.060.050.060.056009.493806009.499.771303809.77 9.91481309.919.949.94第4期受时间和空间的约束能跨区域比较，能够准确反映地理单元在某一功能上的比较优势，本研究采用标准显示性比较优势指数识别国土空间优势功能 2 0-2 1。具体计算公式为：NRCA;,=X;/X,X,X:/X,X,式中：X,为i网格单元的i功能值；X，为所有网格单元的全部功能值之和；X，表示所有网格单元j功能值总和；X，表示i网格单元所有值总和。当NRCA0,该网格单元此功能存在比较优势，反之，则说明其不存在比较优势，该功能竞争力较弱。徐磊，等：京津冀城市群国土空间适宜性格局特征及分区研究N4753结果与分析3.1京津冀城市群国土空间适宜性格局特征以京津冀城市群为研究区域，从生态保护功能、(4)农业生产功能和城镇建设功能出发，选取交通干线可达性、生态系统服务价值、生境质量指数、土壤孔隙度、土壤pH及距乡镇距离指标进行可视化分析，进一步探究京津冀城市群国土空间适宜性格局。京津冀城市群国土空间适宜性指标的格局特征见图1。NN0-2.9052906-70037004-13320km24327-4622113321-2432670140（a）交通干线可达性0-4110412093209330-15100Jkm21200-5420015200-2110070140（b）生态系统服务价值NN高：1Jkm低：070140（c）生境质量指数N0.00%0.01%-49.02%49.03%51.31%151.32%54.63%kn54.64%62.72%70140(d)土壤孔隙度Figure 1 Spatial suitability index of Beijing-Tianjin-Hebei urban agglomeration根据京津冀城市群国土空间适宜性评价指标体系，采用多因素综合评价法分别测算适宜性分值，在ArcGIS软件中采用自然间断点分级法和分位数法0.0-5.05.1-6.86.9-7.5ukm7.6-8.0070140(e)土壤pH图1京津冀城市群国土空间适宜性指标图将国土空间适宜性评价结果进行可视化，将评价结果划分为4个等级，见图2。36m-3382m3383m-6238m6239m-10185m-km10186m-19119m8.1-9.270140(f）距离乡镇距离19120m-37485m476林业与生态科学NN第38 卷N低适宜性一般适宜中等适宜高适宜性-km070140（a）生态保护功能适宜性Figure 2 Spatial suitability evaluation results of Bejing-Tianjin-Hebei urban agglomeration京津冀城市群国土空间生态保护适宜性由西北向东南大致呈“低一一般一中一高一中一一般一低”分布。国土空间生态保护高适宜性地区集中分布在燕山山脉、太行山脉，该区域海拔相对较高，包含5478个网格单元，地貌类型以山地为主，位于河流上游的太行山是天然屏障；此外，该类型区域是滦河等多条河流的发源地，水资源丰富，滋养了山、林、田多种生态要素，影响和调节着该区域的气候变化。生态保护功能一般适宜和低适宜性地区主要分布在京津冀城市群西北部和中南部地区，尤其是中部和东南部的河北平原地区，呈镶嵌分布，包含1347 8个网格单元，主要由古黄河、海河、滦河等水系冲积而成，限制了生态保护。京津冀城市群国土空间农业生产适宜性具有明显的空间异质性。农业生产适宜性高值区主要分布在太行山东麓、燕山南麓和山麓平原地区，包含5888个网格单元，地貌类型以平原、丘陵和盆地为主，水系分布有洋河流域和桑干河流域，农业生产水、肥条件较好，为发挥生产功能提供了便利条件。农业生产一般适宜地区和中等适宜地区混合分布在京津冀城市群中部和南部地区，呈面状分布，该区域干旱缺水，农业生产条件较为贫。农业生产适宜性低值区分布在京津冀城市群西北部，主要分布在低适宜性一般适宜中等适宜-km高适宜性070140（b）农业生产功能适宜性图2 京津冀城市群国土空间适宜性评价结果坝上高原，零星分布在城市群中部和西南部地区，包含56 93个网格单元，地貌类型以丘陵为主，海拔较高，环境恶劣，可利用性较低。京津冀城市群国土空间城镇建设适宜性低值区和高值区以燕山山脉和太行山山脉为界分布在南北两端。城镇建设低适宜性地区主要分布在燕山和太行山2 大山脉以北地带，包含548 7 个网格单元，海拔相对较高，气候较为寒冷，人口密度低，社会经济发展相对缓慢，较难发挥资源优势进行城镇建设。城镇建设一般适宜地区呈网络状分布在京津冀城市群内，在低适宜性地区和高适宜地区交界处集中分布，共包含58 40 个网格单元，在区位上起到过渡作用。城镇建设中高适宜地区集中分布在京津冀城市群中部和东南部，由城市中心呈放射状分布，该区域地势优越，气候环境较好，经济水平较高，基础设施相对完备，适宜城镇建设。3.2京津冀城市群国土空间适宜性空间关联特征基于国土空间“生态保护功能、农业生产功能和城镇建设功能”适宜性水平的正态统计量依据自然断点法将其从高到低分层，生成京津冀城市群国土空间适宜性水平冷热点分布图，见图3。表明京津冀城市群国土空间适宜性水平指数存在明显的空间关联特征，高值区和低值区在空间上都趋于集聚。一低适宜性一般适宜中等适宜-km高适宜性070140（c）城镇建设功能适宜性第4期徐磊，等：京津冀城市群国土空间适宜性格局特征及分区研究477NNN显著热点区热点区域次热点区域不显著区域次冷点区域070140km（a）生态保护功能适宜性Figure 3 Cold and hot spot analysis of spatial suitability of Beijing-Tianjin-Hebei urban agglomeration生态保护功能适宜性显著热点区成片分布在京津冀城市群北部和西部边缘地区，呈块状分布，为燕山山脉和太行山山脉交叉处，水源涵养较好，有利于国家实施各项农业、生态项目，主要涉及城市有承德市、秦皇岛市、张家口市、保定市、石家庄市、邢台市、邯郸市的西部部分地区。显著冷点区和冷点区集中分布在京津冀城市群中南部和西北部。生态保护功能适宜性不显著区域处于显著热点和显著冷点区间隙，视觉上形成明显分界。农业生产功能适宜性热点区域和冷点区域图斑分布都相对扩散，呈块状零星分布。就城市分布而言，显著热点区和热点区域主要分布在张家口市中南部、承德市北部、北京市、天津市和秦皇岛市部分地区，区域地貌类型主要为丘陵和盆地，土壤肥沃，适宜植物生长。显著冷点区呈条带状分布，主要涉及城市有张家口市、保定市、廊坊市、石家庄市和邯郸市。不显著区域广泛分布在京津冀城市群除显著冷热点区域以外，图斑连片分布。京津冀城市群城镇建设功能适宜性显著热点区集中分布在河北省中南部平原地区和北京、天津2 个特大城市，并向北部显著热点区热点区城次热点区域不显著区域次冷点区域冷点区域显著冷点区图3京津冀城市群国土空间适宜性冷热点分析图显著热点区热点区域次热点区域不显著区域次冷点区域070（b）农业生产功能适宜性盆地延伸零星分布。显著冷点区和冷点区域主要分布在燕山和太行山、京津冀北部和西部山地、丘陵地区，图斑均相对集中。城镇建设功能不显著区域分布相对较小，处于冷点区和热点区过渡地带。城镇建设功能适宜性总体表现为东南高西北低，东南部资源丰富，交通便利，为工业发展和城镇建设提供了有利条件。3.3京津冀城市群国土空间适宜性分区运用标准显性比较优势指数，从理论上来说，一个网格区域会有某1种或者2 种显性优势功能,其他功能竞争力较弱。基于网格单元尺度的各项功能适宜性评价结果，应用NRCA模型计算各功能比较优势指数，依据比较优势指数的正负将京津冀城市群国土空间分区优化。由于研究尺度选取的是网格尺度，单元数目众多，为了便于直观分类和可视化表达，根据各功能比较优势指数正负情况对应的优势功能分区,应用数理统计法中的判别分析对其NR-CA值进行赋值，并在ArcGIS中采用自然间断点法将比较优势功能相同的网格单元归属为同一类型区。关于NRCA值赋值表及优势功能分区见表6。140km冷点区域显著冷点区0（c）城镇建设功能适宜性70140-km冷点区域显著冷点区478生态保护功能比较优势Comparativeadvantage of ecological advantage of agriculturalprotection function+一注：由于计算结果数值较小且繁多，为便于分析，将结果的正负情况罗列，根据结果不同组合应用判别分析进行16 赋值，得到对应的优势功能分区。得到京津冀城市群6 种功能类型区，将其按照单一或组合“比较优势功能十优势区”进行命名，分别为生态保护功能优势区、生态保护一农业生产功能优势区、生态保护一城镇建设功能优势区、农业生产功能优势区、农业生产一城镇建设功能优势区和城镇建设功能优势区，分析不同功能特征的聚类区域，针对6 类优势功能分区提出以下优化方案，见图4。生态保护一农业生产功能优势区生态保护一一城镇建设功能优势区农业生产功能优势区农业生产一一城镇建设功能优势区城镇建设功能优势区km70140图4京津冀城市群国土空间适宜性分区图Figure 4 Spatial suitability zoning map ofBeijing-Tianjin-Hebei urban agglomeration生态保护功能优势区：该区域涵盖452 4个网格单元，面积占比约为18.8 6%，主要集中在京津冀城市群北部和西部山脉水系处，多山地，不适宜大规模的农业耕作和城镇建设。如承德市依托燕山山脉形成天然生态屏障，并拥有滦河和伊逊河等众多水林业与生态科学表6 NRCA值赋值表及优势功能分区Table 6 NRCA value assignment table and advantage function partition农业生产功能比较优势城镇建设功能比较优势ComparativeComparative advantageof urban constructionproduction functionfunction一一一+资源，是得天独厚的避暑胜地；该区域在不触碰生态红线的前提下，可利用其优越的自然条件和地域人文特色，打造休闲旅游城市，缓解京津大都市发展压力。京津冀城市群西部涉及到的城市有保定市、石家庄市、邢台市和邯郸市，凭借太行山脉和百花山、狼牙山、白石山等名山以及拒马河、沙河、黄壁庄水库等水系，可大力发展生态旅游产业，作为实验点开展自然观光、生物研学等活动，借助新媒体直播等方N式，助力当地绿色发展。城镇建设功能优势区：该区域涵盖536 4个网格单元，面积占比约为2 2.36%，多分布在京津大城市和京广线沿线。北京市由于经济、人才、教育等要素高度集聚，资源承载压力大，该区可优化城市规划，响应政策积极进行功能转移。廊坊市承接中心大都市，拥有得天独厚的城镇建设优势，应着力高新产业和设施建设，向京津大城市引流人才技术，打造邻京津快速发展新区。雄安新区应发挥首都和