以2011—2021年中部6省87个市、州面板数据为样本，集成三阶段数据包络分析、全局及局部空间自相关方法，分析各城市科学技术创新效率的时空演变特征与集聚格局，探究其动态和静态驱动作用机制。结果显示：科学技术创新效率区域差异程度提高，空间分布逐步由“整体北强南弱、局部差异分布”，向省会核心外围转变；空间集聚性增强，形成因地而异的创新驱动机制；为中小微公司可以提供完善的政策与科技支撑服务，是该地区提高科学技术创新效率的重要路径。
区域协同创新模式是凝聚创新驱动理论和创新驱动战略的科学路径，日益成为推动区域经济高度协调化的“中国智慧”。注重创新驱动发展，强化创新体系和创造新兴事物的能力建设，推动科学技术创新和经济社会持续健康发展深层次地融合，塑造更多依靠创新驱动、更多发挥先发优势的引领型发展，是我国“十四五”期间实现创新型国家建设的重要工作方向。党的十九大以来，为实现供给侧改革和区域产业体系升级，各地区科技要素投入愈发强劲，区域创新合作网络愈发密集，形成了以京津冀城市群、长三角城市群、珠三角城市群为核心的科学技术创新高地，雄安新区、粤港澳大湾区、浦东新区、南沙新区和两江新区等一批国家级新区也逐步建设为驱动区域经济高水平发展的重要引擎。然而，科技资金浪费、科技资源管控低效和创新成果转换率低仍是抑制地区创新技术跃升的内在问题，科技资源跨区域配置低和科学技术创新管理经验交流欠缺仍是阻隔现代化区域创新体系建设的掣肘性因素。如何优化科技资源跨区域合理配置和提高科学技术创新效率，成为推进新型城镇化建设和形成以创新为主要引领支撑的经济体系的重点和难点，开展区域科学技术创新效率评价研究尤为重要。
目前研究人员多采用数据包络分析、随机前沿模型、标准差椭圆和趋势面分析等方法，综合研究长三角、京津冀和珠三角等地区科学技术创新效率的空间异质性，对中西部欠发达或发展中地区的相关研究还较少涉及，难以为该区域产业转型升级和创新技术管理上的水准提升提供科学借鉴。
地区经济规模、全球化水平和政府服务能力等城市要素对科学技术创新效率的驱动影响，近年来得到研究人员的广泛关注，普通最小二乘法和Tobit回归模型被运用于相关研究中，但以地理空间要素为基础的驱动机制探讨尚有不足，难以对区域整体科学技术创新效率的持续健康发展提供科学借鉴。基于此，本研究以中部6省87市、州为研究对象，集成三阶段数据包络分析（data envelopment analysis，DEA）、全局自相关和局部自相关方法分析科学技术创新效率的时空演进特征，并基于地理探测器和地理加权回归（geographically weighted regression，GWR）模型探究影响科学技术创新效率分布的动静态驱动机制。
中部6省是实现“中部崛起”战略的主体区域，包括河南省、山西省、湖北省、安徽省、湖南省、江西省共6省87市、州。该地区面积达102.8万km2，占全国陆地面积的10.7%，2021年GDP共计14.59万亿元，占全国的17.6%，具有经济加快速度进行发展的巨大潜力和优势（图1）。近年来，随着科创资源由国际向国内转移，由东部沿海向中西部地区迁移，中部6省的科学技术创新实力和科学技术研发水平持续提升。2021年该地区公共预算支出中科技支出共计894.38亿元，为优化企业创新和产业升级环境奠定了充足的基础。持续增长的科技投入、产学研高新技术资源云集，推动中部6省逐步成为国家科学技术创新试点的重点区域之一。
首先是结合三阶段DEA模型和探索性空间数据分析，研究中部6省科学技术创新效率的动态演化；其次是以地理探测器和地理加权回归模型分析科学技术创新效率分布的驱动机制，重点是运用地理探测器和地理加权回归模型分别展开驱动因素的动态和静态分析。
结合前人研究成果，建立三阶段DEA模型，排除外部环境对效率评价的干扰，提高科学技术创新效率的测度精准度。其中，第一阶段采用纯技术效率（pure technical efficiency，PTE）模型评价“规模报酬不变条件下”各单元科学技术创新效率、纯科学技术创新效率和规模效率；第二阶段采用随机前沿分析 （stochastic frontier approach，SFA）模型，将管理效率、外因和统计随机因素3种对研究单元科学技术创新效率的影响分离；第三阶段将前阶段调整后的产出数据和原始产出数据重新代入PTE模型中，对科学技术创新效率、纯科学技术创新效率和规模效率进行重新测算。三阶段DEA模型的运算步骤不作展示，具体可参考张宁等的研究。我们从投入、产出及环境变量3个维度建立中部6省科学技术创新效率的评价指标体系（表1）。
科学技术创新活动需要大量的资本支撑和智力投入，从资金和人力2个层面选取全社会研究与试验发展经费（research and experimental development，R＆D）的内部支出（X1）、研发人员数全时当量作为科技创新活动的投入指标（X2）；期刊论文、专利数是常见的科研产物，工业公司盈利收入与科学技术创新成效关联，选取专利申请授权数（Y1）、规模以上工业公司总利润（Y2）为技术创新活动的产出指标；科学技术创新活动受区域经济规模、科教研水平等外在外因影响，选取地区生产总值（Z1）、研究机构全时当量（Z2）和教育支出占一般公共预算支出比重（Z3）作为科学技术创新效率的环境变量。
全局自相关分析：选取全局莫兰指数（Global Morans I，Morans I），用于探测中部6省科学技术创新效率的空间关联结构和空间自相关集聚特性：
式中，xi、xj表示第i、j个市、州的科学技术创新效率值；n为参与空间分析的单元数；wij为空间权重矩阵；S0为所有要素标准差之和。
局部自相关分析：采用安瑟伦局部莫兰指数（Anselin Local Morans I，Local Morans I）分析局部区域的空间自相关程度：
式中，Zi、Zj表示观测值的标准化形式；Ii表示局部空间关联程度，当Ii0，表明区域邻近单元相似值集聚，Ii0表明区域邻近单元非相似值集聚；Ii为0表示区域邻近单元值随机分布；Z（Ii）为局部自相关的Z检验。
地理探测器是中国科学院王劲峰团队开发的空间分异与影响因子探测软件，该模型有因子探测、交互作用探测、风险区探测及生态探测4种类型。本研究选取的因子探测揭示驱动因子对中部6省科学技术创新效率的解释程度。模型如下：
式中，d为中部6省科学技术创新效率的驱动因子；Td为d驱动因子对中部6省科学技术创新效率的解释力，取值介于，该值越大，解释能力越强；h为分类数；Nh和N分别为h类和全部区域的单元数，和σ2分别为h类和全部区域的y值的方差。
地理加权回归（GWR）模型考虑要素的非平稳性，能够衡量驱动因素因地区而异的空间差异。本研究利用ArcGIS平台GWR工具进行地理加权回归分析，地理加权回归模型的方程如下：
式中，（ui，vi）是地区i的空间坐标；yi和xi1，xi2，xi3，…，xik分别是因变量y和解释变量x1，x2，x3，…，xk在（ui，vi）处的观测值；εi是随机误差项；系数β（ui，vi）是关于空间位置的k个未知函数。
式中，x是自变量矩阵；xT是自变量转换矩阵；W是n×n的空间权重矩阵，反映了地理空间位置对参数估计的影响。
科技创新效率评价指标数据主要来自于相关年份《中国科技统计年鉴》《中国城市统计年鉴》，并以河南省、山西省、湖北省、安徽省、湖南省、江西省的省域统计年鉴及87市、州统计年鉴、统计公报和政府工作报告等统计资料加以补充。
表2反映各省科学技术创新效率平均值、变异系数和等级序列情况。整体看来，中部6省2011—2021年的科学技术创新效率整体水平相对来说比较稳定（AVG：0.349→0.341），但区域内城市间的差异程度明显提高 （CV：0.465→0.715）。从省域对比看，各省科学技术创新效率的第一梯队比重（TFS）普遍增长，河南省替换山西省成为首位省份。2011—2021年，山西省科学技术创新效率CV值增长0.43，效率均值由0.671骤降至0.108，效率高值城市接连消失。表明该时段山西省科学技术创新资源的投入冗余问题突出，科技资源利用绩效缩水，既缺乏核心带动的高值城市，城市间科学技术创新管理也缺乏良性协作机制；河南省科学技术创新效率均值虽有下降（AVG：0.52→0.458）但在2021年居于首位，城市内部科学技术创新效率差异也比较小（CV2017：0.547），表明河南省保持较优的科技资源配置优化水平，城市间的科创管理模式互有借鉴，实现了整体创新效率排名的稳步上升。
图2反映出中部6省2011—2021年间科学技术创新效率的空间分布逐步由“整体北强南弱、局部破碎分布”向以省会为中心的核心外围模式转变。2011年，湖南省、江西省等南部地区表现孱弱，北部地区形成以晋城市、吕梁市、漯河市、周口市和芜湖市等为代表的创新效率高地，各省域边缘区位尚未形成创新效率规模集聚。到2021年，武汉市、长沙市、郑州市、合肥市、太原市等省会城市，因科创资金利用、科研机构管理优势成为各省的科学技术创新效率中心，并逐步产生创新管理示范效应和溢出效应，带动周边地区形成科学技术创新效率核心外围组团。
图3反映创新技术效率全局Morans I变化。中部6省科学技术创新效率正相关空间集聚性，随着Morans I逐步的提升（0.170→0.225），空间集聚特性持续加强，Z检验值也始终超过临界值，通过95%显著性水平上的统计检验。图4是局部自相关分析下的（local indicators of spatial association，LISA）聚类图，将研究区域分割为高高（HH）、低低（LL）、低高（LH）和高低（HL）4种类型。2011—2021年，HH型和LL型是中部6省科学技术创新效率的主流集聚类型。其中，HH型由破碎布局向组团集聚演化，2011年，该类型零碎布局于临汾市、荆州市、鄂州市和马鞍山市，到2021年，该类型集聚布局于以开封市为代表的河南省东部和以荆州市为代表的湖北省南部等8市。这反映出河南省和湖北省科学技术创新效率的空间溢出效应和城市创新互动机制随时间推移而逐步增强。LL型空间集聚格局与经济发展水平协同，布局于中部6省的边缘地区。2011年，布局于湖南省西部和安徽省东部6市，到2021年，该类型城市集中布局于以湘西州为代表的湖南省西部地区和以大同市为代表的山西省北部地区。这些城市的经济发展水平相对滞后，目前正处于产业转变发展方式与经济转型的接续期，科创资源利用水平较低，形成科学技术创新效率低值集聚区。此外，LH型集聚区维持破碎状态，由布局于恩施州、神农架林区转移至平顶山市、神农架林区和马鞍山市3地，这些地区受自身产业供血能力不强、经济发展水平落后、管理经验贫瘠等因素制约，与周围地区对比形成科学技术创新效率沉降区。
在中部6省发展实际和数据可得性基础上，选取对外开放水平、科技投入程度、产业现代化水平、企业经营效益、市场消费水平等要素作为驱动因子变量，并以进出口贸易总额占GDP比重、科技支出占财政支出比重、第三产业增加值占GDP比重、规模以上工业公司总利润和社会消费品零售总额占GDP比重等指标分别反映上述驱动因子，建立驱动因子指标体系。利用地理探测器因子探测工具，在对各驱动因子变量的自然断点分类的基础上，分析各驱动因子对于中部6省科学技术创新效率分布的解释力（表3）。结果显示，2011年，中部6省科学技术创新效率驱动因子的解释力排序为市场消费水平、产业现代化水平、企业经营效益、对外开放水平、科技投入程度。这说明，市场的消费水平极大程度地刺激了私营企业和公共部门产业技术升级和科技资源管理更新，促进科学技术创新效率的持续发展。到2021年，主导驱动因子由市场消费水平转为企业经营效益，对外开放水平、科技投入程度和产业现代化水平的解释力均呈现出不同程度的增长。这说明中部6省科学技术创新效率空间分布的驱动模式逐步由单要素主导驱动向多元要素联合驱动转变。其中，企业经营利润增长吸引大量产业科技跟随者和模仿者，慢慢地增加企业内部机制管理和科技投入，促进创新效率的持续发展，而地区的国家化分工、对科技开发的重视程度和现代化产业体系建设均以不同程度贡献地区科学技术创新效率变迁。
利用地理加权回归模型分析2021年驱动因子在各地区的非平稳性特征，以各驱动因子为自变量，以科学技术创新效率值作为因变量，建立中部6省科学技术创新效率的GWR模型。回归结果为，Bandwidth（带宽或相邻点数目）为440357.6，AICc值（性能度量值）为−173.2385，Sigma值为0.0811，各项指标均通过显著性检验，且不存在多重共线，模型拟合程度较优。图5反映了各驱动因子回归系数的地理分布。
1）进出口贸易总额占GDP比重的回归系数有正负之分，整体分布由中间地区向南北两侧递增。回归系数高值区集中分布于山西省，低值区布局于湖南省北部、湖北省南部、安徽省南部、江西省北部的部分城市。这主要是由于近年来这些地区的进出口贸易逐渐做大做强，工程机械、有色金属和特色农产品远销东南亚、拉丁美洲和欧洲，但受外贸产品的技术层次较低和生产流程标准化等因素制约，导致进出口贸易本身的技术需求不高，对于科学技术创新效率的增益效果不明显。而山西省作为资源转型大省，亟待建立新的经济增长点，发展外向型经济对于提振山西省产业经济，增强科技投入和管理上的水准有重要促进作用。
2）科技支出占财政支出比重的回归系数总体呈现自北向南递减，高值区分布于山西省南部、河南省西部、湖北省北部5市。表明这些地区近年来一直在优化科研项目经费管理上的水准，坚持推进财政科研项目经费管理改革，科研单位、科研人员和社会企业的创新活力充分释放，为该地区的科学技术创新效率提升奠定了充足基础。
3）第三产业增加值占GDP比重的回归系数为负值，高值区在地理空间上分布较为破碎。这表明第三产业的发展对提升中部6省科学技术创新效率的贡献不足。原因可能在于中部多数城市仍是以生活性服务业为第三产业的主体形式，产品技术层次与东部沿海地区尚存差距，导致科学技术创新的现实需求不高，在某些特定的程度上制约了科学技术创新效率的进步。
4）规模以上工业公司总利润的回归系数自东南向西北递减。高值区布局于江西省东部、安徽省南部，这些城市地处长江经济带，邻近粤港澳大湾区和长三角城市群，具有对外交流的区位优势，近年来逐步打开国际机械设备、轻工纺织市场，获得巨额利润收入，直接影响相关企业对科学技术研发和资金管理的重视程度，驱动区域科学技术创新效率的提升。
5）社会消费品零售总额占GDP比重的回归系数以正向分布为主。空间分布由北向南递减，高值区布局于山西省南部和河南省西北部11市。究其原因主要在于随着该地区经济水平和居民可支配收入的提升，对新型科学技术产品的消费需求慢慢地增强，为扩大科学技术产品市场容量和科技资源的优化配置起到正向推动作用。
科技创新效率涉及科学技术研发、产品生产和资源管理等诸多领域，是影响城市持续包容发展的重要支撑系统。本研究集成多种数理实证模型和空间计量方法，揭示了中部6省科学技术创新效率的时空演变特征及空间分布的内在驱动机制，能为丰富科学技术创新效率研究范畴和空间界限起到一定的贡献作用。
未来中部6省应把科学技术创新的工作重心放在优化科技投入质量和改善企业经营环境等方面。一方面，要加快落实国家创新驱动战略政策，协调好科学技术研发投入和科学技术创新管理服务的均衡关系，不断驱动科学技术创新效率的发展；另一方面，对标东部发达省份，积极融入国家“一带一路”倡议、中部崛起及城市群建设的浪潮，一直在优化市场环境，为中小微公司可以提供完善的政策服务和科技支撑服务，鼓励企业提升自主创造新兴事物的能力和制定产品技术改造计划，以高质量的企业经济发展带动区域创新实力和创新效率的不断的提高。另外，本研究仍存在指标选取不充分、研究尺度偏宏观和机制构建不完善等不足，在研究数据可得条件下，将科学技术创新效率研究深入到中部6省内部城市群、都市圈等多元城市体系，是下一步工作的重点方向。
以中部6省87市、州为研究对象，集成三阶段DEA模型、全局自相关、局部自相关、地理探测器和地理加权回归模型等方法，分析2011—2021年中部6省科学技术创新效率的空间集聚演化及驱动机制，得出以下结论。
1）科学技术创新效率整体水平相对来说比较稳定，区域内部城市的差异程度明显提高，空间分布模式逐步由“整体北强南弱、局部破碎分布”向以省会为中心的核心外围模式转变。
2）科学技术创新效率分布的空间集聚性明显地增强，HH型和LL型是主流集聚类型。其中，HH型由破碎布局向组团集聚演化，向河南省东部和湖北省南部加速布局，LL型空间集聚格局与经济发展水平协同，布局于中部6省的边缘地区，平顶山市、神农架林区和马鞍山市3地形成LH沉降区。
3）驱动因子动态分析显示企业经营效益取代市场消费水平成为主导驱动因子，逐步由单要素主导驱动向多元要素联合驱动转变；地理加权回归分析显示，各驱动因子在不同城市的影响具有非常明显差异，形成因城市而异的科学技术创新效率静态驱动机制。
作者简介：王思成，贵州大学建筑与城市规划学院，讲师，研究方向为城乡韧性ECO规划；高明明（通信作者），贵州大学建筑与城市规划学院，浙江大学建筑工程学院，讲师，研究方向为区域产业经济协调发展与治理。
《科技导报》创刊于1980年，中国科协学术会刊，主要刊登科学前沿和技术热点领域突破性的成果报道、权威性的科学评论、引领性的高端综述，发表促进经济社会持续健康发展、完善科学技术管理、优化科研环境、培育科学文化、促进科学技术创新和科技成果转化的决策咨询建议。常设栏目有院士卷首语、智库观点、科技评论、热点专题、综述、论文、学术聚焦、科学人文等。