中国东部沿海四大城市群协同创新效率综合测度及影响因素研究

孙振清，李欢欢，刘保留

(天津科技大学 经济与管理学院，天津 300222)

0 引言

城市群作为国家实现新型城镇化和工业化发展中较高形式的产物，逐渐取代中国传统省域经济发展格局，成为未来地区经济发展的重要趋势[1]。在政策指引和现实需求的双重推动下，城市群经济发展呈现经济增速快、覆盖面广和区域联动性强等特征[2]，发展不均衡问题仍然存在。其中，东部沿海城市群因资源禀赋、创新能力和区位优势，成为中国对外开放程度和经济发展水平较高的典型城市群，经济发展也呈现区域化特征。单个主体或者城市的技术创新能力可能无法有效应对瞬息万变的竞争环境，城市群间的协同创新则是实现区域一体化的关键，有助于资源在城市间各生产环节协同整合[3]，促进各城市优势互补。因此，提升东部沿海城市群协同创新能力，促进区域间技术和知识流动，实现创新资源整合，强化协同创新多维溢出效应，对于提升技术创新水平和经济发展潜力有着至关重要的作用。另外，在中国创新资源紧缺的情况下，如何有效解决效率问题，优化创新要素配置，对于推动国民经济发展具有重要现实意义。为此，本文立足于协同发展角度，对东部沿海城市群(京津冀、山东半岛、长三角和珠三角)区域协同创新效率进行测度，并分析其内在机理和影响因素，既符合国家发展的战略需求，也为其它城市群实现区域经济和资源协同发展提供借鉴。

1 文献综述

查阅相关文献发现，针对协同创新的研究主要集中在协同创新效率测度以及影响因素识别、作用等方面，具体来说，主要对协同创新内涵、动因、模式以及作用机制进行探讨。协同创新内涵方面，吴钊阳等[4]认为，协同创新是企业、政府和知识研究机构等为实现重大科技创新而开展的大跨度创新组织模式。动因方面，学者们从外部资源获取、企业绩效提高和风险共担等角度进行探讨，发现企业采用协同创新模式可以进行信息流和知识传递、分享等多维互动，加速企业技术创新能力和绩效提升[5]。模式方面，学者们主要从研发协同、专利转让以及技术创新等层面展开分析，如何郁冰[6]提出“战略—知识—组织”三重互动的产学研协同模式，为实现产学研协同创新奠定了理论基础；Bjerregaard[7]认为，研发协同有利于外部知识和技术整合利用，降低成本和风险，促进科研合作和信息共享，从而实现企业协同创新发展。作用机制方面，杨晓昕等[8]通过探讨一些地区推进的“创新驿站”发现，新型区域协同创新机制需要构建集研究型大学、中介机构和企业以及政府“四位一体”的区域协同创新体系。协同创新效率测度以及影响因素方面，主要从相同创新主体和不同创新主体两个角度加以分析：相同创新主体角度侧重于高校、企业和科研院所等机构，如尹洁等[9]采用问卷调查方式，分析高校科研创新团队知识共享机制下创新团队内部知识共享绩效影响因素；Kim等[10]以韩国中小型企业为案例研究对象，采用多元回归方法证实，中小型企业开放式创新活动对其创新产出的影响存在差异性；温珂等[11]以中国101家公立研究院所为样本，发现协同创新能力与绩效之间存在正向促进关系。不同创新主体角度侧重于模型建立和度量，如段云龙等[12]运用超效率DEA模型和网络DEA模型测算各省市产学研之间的协同创新效率，发现区域共生系统下各主体创新效率呈现不均衡状态；蔡翔等[13]通过构建三螺旋强度模型，实证分析大学—企业—政府协同创新效率及影响因素。此外，Fuentes等[14]、Li等[15]、叶斌等[16]从不同角度对区域协同创新进行了研究。

随着空间经济学的兴起和发展，学者们开始关注区域协同创新外溢问题。如赵增耀[17]、白俊红[18]、刘军[19]和孙大明等[20]从空间维度厘清区域协同创新空间外溢效应和机制，重点考察其对区域创新绩效、产业集聚和升级的影响，但基于城市群视角的研究较少。以往学者往往以省际(区域)尺度为研究对象，忽略了不同地区发展内质的不均匀性。本文以城市群为研究对象，更能突出不同城市间的内在关联性。与以往研究相比，本文可能的创新贡献在于：

(1)在研究内容方面，将京津冀、山东半岛、长三角和珠三角地区统一纳入分析框架，重点考察不同地区协同创新效率变化趋势，弥补了传统研究未考虑城市群发展对区域协同创新效率影响的不足。

(2)在变量测度方面，采用三阶段DEA模型测算2005—2018年中国四大城市群协同创新效率，对处于生产前沿的决策单元进行分析。

(3)在影响因素方面，本文采用面板Tobit回归模型实证分析不同城市群协同创新效率的主要影响因素，以期为地区因地制宜采取相关政策，实现协同发展提供依据和思路。

2 研究方法与数据说明

2.1 研究方法

2.1.1 三阶段DEA模型

传统DEA模型未考虑环境因素和随机噪声对决策单元效率评价的影响，导致测算结果有所偏差，而三阶段DEA在测算过程中能够消除自身不可控因素对效率的影响。

第一阶段：传统DEA模型。Charnes等[21]最早提出被称为数据包络分析的方法，用于评价相同部门间的相对有效性，因而被称为DEA有效，该模型即CCR模型。从生产函数角度看，这一模型是用来研究具有多个投入，特别是具有多个产出“生产部门”的理想且卓有成效的方法。DEA模型分为投入导向和产出导向；Banker等[22]提出BCC模型，测度规模报酬及可变条件下决策单元有效性问题，将综合技术效率(crste)分解为规模效率(scal)和纯技术效率(vrste)，具体关系为crste=scal×vrste。本文选择BCC模型对研究区域内X个城市的初始投入与产出进行测算。

(1)

其中，i=1，2，...，m；r=1,2，...，s。k=1,2，...，n为决策单元数，m、s分别代表投入指标和产出指标个数，xik(i=1，2，...，m)为第k个决策单元第i个投入要素，yik(r=1,2，...，s)为第k个测算单元第r个产出要素，θ为决策单元有效值。若θ=1且s+=s-=0，则决策单元DEA有效；若θ=1且s+≠0或s-≠0，则决策单元弱DEA有效；若θ<1，则决策单元非DEA有效。

第二阶段：利用相似SFA回归剔除环境因素和统计噪声。第一阶段的DEA模型无法消除不可控因素对效率值的影响，因此，第二阶段通过引进SFA模型，将第一阶段的松弛变量分解成包含环境因子、随机因素和管理无效率的3个自变量函数。以投入导向为基础构建类似SFA回归函数。

Sni=f(Zi;βn)+vni+μni;i=1,2,...,N;n=1,2,...,N

(2)

为剔除环境因素和随机因素对效率测度的影响，将所有决策单元调置于相同的外部环境下，基于最有效预测单元，以投入量为基础对上述公式进行如下调整：

(3)

第三阶段：调整后的DEA模型。将调整后的投入产出变量再次带入第一阶段的传统DEA模型，测算各决策单元第三阶段的协同创新效率。此时，效率已经剔除环境因素和随机因素的影响，得到的效率值是相对真实准确的。

2.1.2 Tobit回归模型

Tobit回归模型是基于极大似然估计法建立的，可以较好地规避参数不一致、有偏性等问题。因此，采用因变量受限的截断Tobit回归模型对不同城市群协同创新效率影响因素进行分析，具体模型形式设定如下：

(4)

2.2 指标选取与说明

从经济学角度而言，从事某项生产活动的基本要素是劳动和资本，协同创新活动亦是如此。基于协同创新效率的基本概念以及现有研究成果，考虑协同创新主体之间的关系，本文从投入和产出两个层面出发，选择投入产出指标如下：

投入指标方面，考虑到区域协同创新主体构成主要是高校、科研机构、企业和政府，从资本和人力劳动两个方面进行考量，在人力劳动方面参考廖名岩等[23]的做法，采用地区总R&D人员减去企业独立R&D人员、高校应用研究R&D人员和科研机构应用研究R&D人员的值衡量。对于资本投入，参考谢子远等[24]的研究，选取高校科技活动中来自企业的经费与科研机构科技活动中来自企业的经费之和作为协同创新效率投入要素指标。对于项目要素投入，由于企业是协同创新主体，在协同创新活动中占据主导地位[25]。因此，协同创新项目要素投入指标为企业R&D项目总数减去企业 R&D独立项目数。

产出指标方面，协同创新活动的目的是为获得知识产出或者具有市场前景的新产品或者创新技术。因此，本文主要从经济收益和创新成果两个方面对产出指标进行考量。经济收益方面，选取科技创新研发的新产品销售收入；创新成果方面选取科技创新活动的成果展示专利数，由于所申请的专利并不一定全部通过审查，选择专利授权数量更能反映出一个地区的创新水平。因此，选取新产品销售收入和专利授权数作为区域协同创新产出指标。

第二阶段主要计算对沿海四大城市群协同创新效率产生显著影响、但样本无法主观控制和改变的要素，即外部环境变量，主要包括研发支持力度、产业结构、区域经济发展水平、对外开放程度及政府财政支持力度。具体来说，研发支持力度对区域创新水平有极大的支持作用，以研究和试验发展为主的资金支持是重要方式，因而选择研究和试验发展(R&D)经费投资强度衡量；区域第三产业占比往往更能体现一个区域的产业升级水平，本文以第三产业与第二产业占比即产业结构高级化测度；区域经济环境采取人均GDP衡量，反映区域经济发展水平；对外开放可以引进国外先进知识与创新技术，提升科技创新水平、促进要素流动，进而有利于协同创新发展水平提升，本文采取城市进出口额占城市GDP的比重衡量；政府作为区域协同创新主体之一，其制定的相关政策以及提供的创新资金会对区域协同创新效率产生影响，因而选取财政科技支出占财政总支出的比重衡量政府财政支持力度。具体指标见表1。

表1 区域协同创新效率测算相关指标

2.3 研究区域与数据来源

本文借鉴中国城市群现有研究成果[2]，选取京津冀城市群、山东半岛城市群、长江三角洲城市群和珠江三角洲城市群4个沿海城市群为研究对象，包含8个省份30个城市，如表2所示。

表2 东部沿海城市群范围及选取依据

鉴于数据可得性与准确性，样本时间跨度为2005—2018年，数据来源于《中国统计年鉴》《中国环境统计年鉴》《中国科技统计年鉴》《中国城市统计年鉴》以及相关地市国民经济和社会发展公告。针对缺失数据，根据实际情况采取外延法和均值法计算得出。

3 结果分析

3.1 第一阶段传统DEA测算结果

第一阶段运用BCC模型，采用DEAP2.1软件包对2005—2018年我国东部沿海城市群30个城市协同创新效率进行测度，分别得到综合效率、纯技术效率和规模效率，表3为各项效率的平均值。

表3 第一阶段2005—2018年四大城市群协同创新效率平均值

在不考虑外部环境与随机误差影响的条件下，东部沿海四大城市群30个城市中有7个城市在14年中处于随机前沿面，分别为北京、沧州、烟台、苏州、上海、广州和中山，其它城市的综合效率与最优值之间仍存在差距。总体来看，珠三角地区协同创新的综合效率值最高，京津冀和长三角地区次之，山东半岛的综合效率值最低，说明珠三角和京津冀地区将资源投入转化为协同效益的能力较强。纯技术效率呈现出与综合效率相同的变化趋势，表明珠三角和京津冀地区的投入要素配置合理且利用效率高。如表4所示，从时间变化上看，四大城市群的综合效率和纯技术效率呈现波动上升趋势，且在研究区间内有显著变化差异。其中，长三角和珠三角地区的综合效率呈现先下降后上升的“V”型演变格局，在2008年分别达到最低值0.787和0.820后稳步上升。截至2018年，两个地区的上升幅度分别为7.8%和8.4%，说明二者资源配置逐渐趋于合理。四大城市群的规模效率一直比较平稳，说明东部沿海四大城市群协同创新发展集聚水平较高。

表4 第一阶段四大城市群协同创新效率时序变化

3.2 第二阶段：相似SFA分析

第一阶段利用Deap2.1软件测算效率值可以得到各城市投入要素的松弛变量，投入松弛变量会受外部环境因素、随机干扰项和内部管理无效率等方面影响。本阶段以投入松弛变量为被解释变量，以5个环境变量为解释变量，采用相似SFA回归模型，运用Frontier4.1分别估算2005—2018年环境变量对理想投入与实际投入差额的影响。受篇幅受限，选取2018年的回归结果进行分析，如表5所示。

表5 第二阶段SFA回归结果

表5回归分析中，环境变量的SFA回归系数大多通过了1%的显著性水平检验且γ值趋近于1，说明在混合误差项中，管理无效率项对投入松弛变量的影响比随机因素的影响更大，进一步表明模型设定具有一定的可靠性，外部环境因素对城市群协同创新投入效率存在影响。因此，对原始投入值作出调整是合理且必要的，具体结果分析如下：

研发支持力度对劳动、资本、项目投入松弛变量的影响为负且通过1%的显著性水平检验，说明随着研发支持力度加强，各城市在科技人员、研究经费及项目上的投入可以有效转化为专利和财政收入。目前，国家不断加大新产品研发和科技创新支持力度，致力于推动龙头企业技术进步，扩大科研专项资金规模，使得城市创新水平不断提升。随着区域协同发展政策深入推进，创新要素在城市间充分流动所产生的协同效应使相同的要素投入获取更多的创新产出，因此研发力度增大可以减少劳动力与资金投入冗余。

产业结构对资本、项目投入松弛变量的影响回归系数为负，说明产业结构高级化水平越高，劳动力、资本及项目投入冗余越少，越有利于提高城市协同创新效率。我国正处于产业结构转型的关键期，以科学技术创新为主的新型第三产业正在崛起。第三产业比重增加会减少劳动力、资本及项目投入冗余量，避免浪费现象发生。同一相关关系下，产业结构与资本、项目投入松弛变量的回归系数大于前者与劳动力投入松弛变量的回归系数，充分说明资本和项目利用效率受第三产业比重影响更大，与我国发展阶段相符，第三产业增加会带来协同创新效率提升。

区域经济水平与劳动、资本和项目投入松弛变量正相关，表明城市人均GDP对劳动力投入松弛变量的影响不显著，而人均GDP提高会导致资本和项目投入冗余增加，从而降低城市协同创新效率。究其原因在于人均GDP增加对资本利用效率的作用存在门槛效应[26]，只有达到门槛值时，人均GDP增加才会促进资本利用效率提升，进而提升城市协同创新效率。

对外开放程度与劳动、资本和项目投入松弛变量的回归系数为负且均通过显著性检验，表明对外开放程度提升会吸引外国创新资源流入，加速创新领域发展，避免国内传统工艺造成的资源浪费，使得劳动力、资本和项目投入冗余减少，有利于节约劳动力和资本。

政府财政支持力度与劳动、资本和项目投入松弛变量的回归系数为正且通过1%的显著性水平检验。政府作为区域协同创新主体之一，在提升创新效率方面有强烈的意愿，但投入要素在各主体间的利用率并不高，导致部分投入要素未被充分利用，进而增加了投入冗余。

综上所述，外部环境因素会使不同城市在不同环境下的协同创新效率产生较大变差，故需要对原来的投入变量进行调整，剔除环境变量中的随机因素，保证各城市处于相同的外部环境，进而得出更加准确的城市协同创新效率。

3.3 第三阶段：调整后的DEA计算结果

将第二阶段调整后的投入数据与原始产出数据再次导入传统DEA模型进行分析，得到第三阶段东部沿海30个城市协同创新综合效率值、纯技术效率值和规模效率值，以及与第一阶段相比的改进情况，测算结果如表6所示。

表6 第三阶段2005-2018年四大城市群协同创新效率平均值

在剔除外部环境因素和随机误差因素后，东部沿海30个城市协同创新效率的综合效率值由0.833升至0.836，纯技术效率值由0.907降至0.899，规模效率值由0.915升至0.930。可以看出，综合效率与规模效率变化趋势一致，说明综合效率的变化主要取决于技术创新水平。调整后的规模效率值高于调整前水平，表明城市间协同合作会形成一定的规模效应，进而提升规模效率值。4个城市群中，不同城市协同创新效率出现了较大变动，除山东半岛综合效率值较第一阶段有所上升外，其它3个城市群的综合效率值呈现下降趋势，说明城市群协同创新效率值受其所处外部环境的影响。如图1所示，就综合效率而言，北京、烟台、苏州、广东、中山5个城市处于前沿面，较第一阶段减少2个，沧州和上海由于受外部环境的影响综合效率值有所下降，退出效率前沿面。整体上看，综合效率值上升的有14个城市，下降的有11个城市。

图1 四大沿海城市群各城市协同创新综合效率、纯技术效率与规模效率演变

对比第三阶段与第一阶段改进比例可以看出，天津、唐山、威海等14个城市的综合效率值增加，说明外部环境和随机干扰项对上述城市协同创新效率存在较为不利的影响。廊坊、南通、杭州等11个城市综合效率值呈现下降趋势，说明上述城市协同创新效率受外部环境的干扰较大，进而暴露出城市自身创新水平相关问题。此外，北京、广州等5个城市的改进比例不变，表明外部环境因素以及随机干扰项对协同创新投入转化利用率的影响较小，城市协同创新水平较高。

3.4 东部沿海城市协同创新效率及时空特征分析

为深入分析四大城市群之间协同创新效率变化趋势，绘制综合效率、纯技术效率与规模效率变化趋势图(见图2)。总体来看，四大城市群的协同创新效率呈现波动上升趋势。从城市群层面看，四大城市群的协同创新效率存在一定差异，呈现珠三角>京津冀>山东半岛>长三角的发展格局，即珠江三角城市群协同创新效率最高，这与珠三角率先通过改革开放获得的有利环境与机遇密切相关。珠三角凭借国家政策支持以及自身高新技术产业优势，在区域协同创新发展中取得重要突破，京津冀城市群在“京津冀协同一体化”发展政策的推动下，城市发展处于较高水平。不同于以往研究[23]，长三角城市群协同创新效率在四大城市群处于较落后，可能是由于较高水平的技术创新尚未得到有效转化。协同创新属于科技创新活动范畴[28]，易受区域经济发展水平、国家政策等外在因素影响[29]。最早提出协同发展政策的京津冀城市群前期协同创新效率呈现上升趋势，其它3个区域同时期协同创新效率则呈现平缓甚至下降趋势。之后，随着区域协同发展、国家创新驱动发展等一系列战略实施，推动传统发展方式向依靠技术进步、知识积累以及高素质劳动力转变，四大城市群协同创新效率开始稳步上升。从图2可以看出，综合效率值与纯技术效率的波动整体上趋于一致，规模效率值趋于平稳。由此可以认为，四大城市群协同创新发展规模为城市群协同创新发展提供了有利条件。

图2 城市群协同创新综合效率、纯技术效率与规模效率演变趋势

为了更直观地分析沿海四大城市群内城市协同创新效率时空演变特征，采用Arc GIS10.2中的自然断点法将协同创新效率分为3类，得到2005年、2011年、2018年城市综合效率值分级分类情况，如表7所示。

由表7可见，京津冀城市群协同创新效率空间分布呈现以北京、天津为中心，并向周边辐射的递进式演变格局，秦皇岛市由2005年的Ⅲ级上升为2011年的Ⅱ级，且其与大部分城市均位于Ⅰ-Ⅱ级；山东半岛城市群主要表现为沿海城市协同创新效率高于内陆城市的空间分布格局，高值区位于东部沿海一带，青岛市和潍坊市呈现等级上升趋势，威海市则呈现先上升后下降的演变格局；长山角城市群内部协同创新效率下降的城市较多，呈现高值区由沿海向内陆减少的总局势，展现出沿海城市可以有效吸引国外先进技术与管理模式的优势；珠三角城市群协同创新效率整体水平较高，截至2018年，Ⅰ级城市数量上升至5个，但惠州市由Ⅱ级下降至Ⅲ级，协同创新发展水平较不稳定。

表7 四大城市群协同创新综合效率分类等级

3.5 Tobit回归的结果分析

运用Stata14.0软件，采用面板Tobit模型对我国沿海四大城市群协同创新效率影响因素进行分析，回归结果见表8。

表8 Tobit面板回归结果

表8是采用Tobit模型的面板回归结果，可以看出，就研发投入支持力度而言，珠三角地区研发投入水平更高，说珠三角地区创新实力明显强于其它3个城市群，山东半岛研发实力和创新水平仍有待提升；产业结构对四大城市群均具有负向影响，表明地区产业结构调整和优化仍面临较大压力，不同城市群的产业结构仍处于不平衡状态；对于人均GDP的影响，长三角和珠三角地区由于地区经济发展水平较高，对区域协同创新效率的积极作用更加显著；对外开放程度和地方财政支出对区域协同创新效率的影响存在明显差异，京津冀和山东半岛城市群由于对外开放程度、财政支出水平较低，并不能有效提高协同创新效率水平，长三角和珠三角地区凭借地理优势和经济实力对协同创新效率具有显著促进作用。

4 结语

4.1 研究结论

本文研究时间跨度为2005-2018年，选取中国东部沿海四大城市群30个城市相关数据，运用三阶段DEA模型测度其协同创新效率并进行比较分析，得到以下结论：

(1)在剔除外部环境与随机因素的影响后，四大城市群协同创新效率均发生了较为明显的变动，证明采用三阶段DEA测度协同创新效率可得到较为客观的结果。调整后，除山东半岛城市群外，其它三大城市群的平均综合效率值均较第一阶段下降，且纯技术效率的变动趋势与综合效率相同，区域协同创新效率提升的关键在于创新技术利用率提升。

(2)外部环境因素和随机因素对协同创新效率具有显著影响，区域经济水平和政府财政支持力度提升会为城市技术创新提供一定资金支持，但是发展动力不足或技术转化率低会对城市协同创新发展产生不利影响；研发支持力度、产业结构和对外开放水平可以有效提升劳动力与资本利用效率，加上国家推动技术创新协同发展的相关政策支持，会对城市群协同创新发展具有正向影响。

(3)我国四大城市群协同创新效率均呈上升趋势，但内部城市等级存在空间不平衡性。受不同时期国家政策以及自身发展水平影响，四大城市群协同创新效率呈现出珠三角>京津冀>山东半岛>长三角的发展格局。城市群内协同创新效率呈现出核心城市以及沿海周边城市较优，内陆城市较低的时空分布特征，除处于前沿面的5个城市外，其它城市协同创新效率均有较大提升空间。

(4)研发支持力度和区域经济发展水平提高有助于四大城市群协同创新效率提升。其中，珠三角系数最高说明其创新实力最强；对外开放度对京津冀和山东半岛协同创新效率具有抑制作用，进一步说明引进先进国外技术对本区域创新发展具有重要作用。

4.2 政策建议

(1)依据四大城市群协同创新效率测度结果，建议采取如下差别化措施：对于珠三角、京津冀城市群,应更加注重技术创新与制度创新，不断强化城市群内部与城市群之间创新主体协同合作；对于山东半岛城市群,在加大投入的同时，应建立科学化协同模式与合理有效的监督管理制度，通过加大协同创新投入增加创新产出；长三角则应该在保证投入不变的情况下提高要素利用率，切实提高城市协同创新效率。

(2)适度增加科研政策及资金支持力度，提高自身研发能力。通过竞争择优吸引高素质人才，注重技术积累以及对共享知识的消化与吸收。

(3)提升产业结构高级化水平，积极推动以技术创新为主的新兴产业发展，逐步实现产业结构调整，借助创新驱动发展战略，提升城市群整体协同创新效率。

(4)提高城市开放程度，增强外资吸引能力，一方面，政府应合理减少外资准入限制，增加外资投资领域与机会；另一方面，要制定相关财政支持政策，优化投资运营环境，为外资投入提供支持与保障。