摘要：基于地理统计和空间计量实证方法，研究中原经济区2000年和2013年的全局与局域空间就业格局演变。结果表明，从全局看，2000年和2013年均为“背空”就业格局，除初始的3个中型就业聚集中心外，未形成大型就业聚集中心，就业核单元年均复合增长率为1.37%。从局域看，中原各市2000年存在显著的东进效应和不明显的南下效应，2013年则存在显著、势均的北上、西进效应。因此，中原经济区应加强地区协同合作来改善就业格局。
关键词：中原经济区;卡斯特经-纬度空间扩张模型;就业格局
中原经济区涵盖5省30个市，人口1.79亿(侯永志， 刘云中，2013)，区域就业格局决定地区繁荣或衰落。实际上，地区就业格局早已引起了学界关注。如Blanchard等(1992)、Carlino等(2002)发现美国地区就业呈现经济繁荣时劳动力流入，衰退时流出;Giuliano等(2007)发现Los Angeles region from 1980， 1990 and 2000， we identify employment，洛杉矶主要就业区呈现内部稳定、向外扩散趋势，需要用动态模型进行刻画。此外Holm和Qstergaard (2013)、Mameli 等(2014)、Dutsch和Struck(2014)等都通过描述、计量等方法揭示就业格局演变。国内学者大多仍沿传统就业转移主线研究就业转移，但也有学者关注就业格局。如殷李松、卓翔芝(2010)采用生物系统论研究方法研究发现东部正遭遇“就业难与招工荒”并存;蔡窻(2013)从农民工和城镇职工两大群体角度，揭示了劳动力市场上的结构性失业风险;牟宇峰(2014)运用产业格局分析法研究发现长三角核心区的就业格局正向多中心演变;罗奎等(2014)采用地理加权回归发现中国地级市就业以大城市为主导、快速增长并向周边扩散。虽然文献研究成果丰富，但至今没有文献对中原经济区的就业格局展开研究，有鉴于此，本研究综合地理统计、核密度分析其格局演变，为科学决策提供依据。[1]
一、研究方法与数据来源
(一) 研究方法
1、全局空间就业格局研究方法：核密度分析制图法
核密度分析属于非参数检验方法，目的是对给定数据x1，x2……xn，估计出概率密度函数：
[^][p(x)]=K()(1)
式(1)中：h为带宽(步长)，主要是为了控制核密度估计的平滑程度， 一般可由软件自动选择。n为观测值个数，K(u)为核密度函数，其种类很多，常用的有均匀核、二次核、三角核、三权核等，ARCGIS软件中以Silverman(1986)描述的二次核函数为基础。其公式如下：
k(u)=(1-u2)I(|u|?1)(2)
2、 识别局部空间就业格局特征的空间计量法：卡斯特经-纬度空间扩张模型
对空间就业格局，来自解释变量的低阶近邻效应将高于高阶近邻，局部空间就业格局特征难以从全局就业格局中推导出来。对此要采用局部空间计量扩张模型进行分析，本研究采用卡斯特经-纬度空间扩张模型，从解释变量沿经度方向(南至北)与纬度方向(东至西)在近邻域上分解空间依赖或外部效应，模型形式如下：
Li=α+Pi+β0+PiZXi βxi+PiZYi βyi+εi(3)
ε～N(0，σ2In)(4)
式(3)、(4)中i代表地市，L代表就业，P代表人口，ZX代表经度，ZY代表纬度，参数β0、常数项α代表非扩张变量影响，βxi、βyi则代表解释变量在经-纬度空间中对就业的扩张效应，βxi为沿经度(由南至北)的扩张效应，即北上效应，若为负，则是南下效应，βyi代表沿纬度(由西至东)扩张效应，即东进效应，若为负，则是西进效应，ε为正态分布。 [2]
(二) 数据来源
本文最小统计单元为地级市，2000年就业数据来自全国第五次人口普查数据(国务院人口普查办公室，2003)，2013年就业数据来自省级统计2014年鉴，地区面积数据来自全国第六次人口普数据(国务院人口普查办公室，2012);经纬度数据来自国家基础地理信息系统数据库。
二、实证分析
(一) 中原经济区全局空间就业格局演变分析
对全局空间就业格局演变，采用核密度分析法(搜索半径100公里)估值就业密度，形成平滑的空间曲面。栅格�色越深说明就业密度越高。若出现2个以上核单元形成深色的圆形或带状区域，则表明出现了区域性就业聚集。对核单元，色彩加深核在成长，反之衰减。[3]
面向正北方位，中原经济区全局空间就业格局为：①2000年为“背空”就业格局，南阳、信阳、运城与三门、阜阳、邢台与邯郸、亳州至淮南连成就业带，开封与周口为空心，历经13年，至2013年仍为“背空”格局，不过，2013年相比2000年各地级的就业核有所成长。②各单元就业核分散成长，除初始的3个中型就业聚集中心，即邢台―邯郸、郑州、漯河―周口―驻马店，并未形成区域性的大型就业集聚中心。③无论是核单元，还是3个中型就业中心成长十分缓慢， 2000年为314.42人/平方公里，2013年为384.13人/平方公里，年复合增长率为1.37%。由此可见，13年来中原各市单元就业核成长缓慢，成长差异微小，没有出现大的就业中心。
(二)中原经济区局域空间就业格局演变分析
1、实证结果
卡斯特经-纬度空间扩张模型中，均以就业为因变量、人口为解释变量，借助MATLAB估计模型得到基准模型为：
式(5)、(6)中参数含义为：①常数项代表是否累积了充足资本保持就业的增长。②人口前的系数项为非扩张的就业效应，由区域内的资源禀赋状况引起。正号为增溢，表明人口增长有相应资本深化，引起就业增长，反之亦然。③为人口沿经纬度方向对就业的扩张效应。经度方向从南至北，正号为北上效应，负号为南下效应;纬度方向从西至东，正号为东进效应，负号则为西进效应。若大量地级市扩张效应同向，表明受到共同干扰，并且对增长的核单元不利，因为资本资源充裕，劳动力被吸出不利于经济成长，反之亦然。[4]
2、实证分析
(1)中原经济区局域空间就业格局基准效应演变
式(5)表明2000年中原地区就人口而言：①资本广化能为就业提供支撑(214.3293);②受资本深化制约，就业存在衰减(-1.7467)，不过影响不显著;③人口沿经度方向对就业的基准扩张效应为负，即存在南下效应(0.0072)，不过不显著;④人口沿纬度方向对就业的基准扩张效应为正，即存在东进效应(0.0053)且显著。由此中原在2000年受到东部经济的强大干扰，以及不显著的南部经济的干扰，使中原各市就业核的增长缓慢且无法进行核聚集，没有出现大的就业中心，只有邢台-邯郸、郑州、漯河―周口―驻马店3个中型就业中心。
