工业机器人应用、劳动保护与异质性技能劳动力就业
2022-08-04胡嘉琪
明 娟，胡嘉琪
(广东工业大学 经济与贸易学院，广东 广州 510520)
一、引言
近年来，我国工业机器人应用规模呈现“井喷式增长”态势。国际机器人联合会(International Federation of Robotics，IFR)全球工业机器人报告显示，截至2019年，我国工业机器人保有量以78.3万台的数量稳居世界第一，约占世界工业机器人总保有量的29%。工业机器人广泛应用，其可能带来的“机器替代劳动”的技术性失业风险也引起广泛关注。部分研究证实：工业机器人应用对中国制造业就业总量产生冲击。工业机器人具有鲜明的技术进步特征，对就业可能产生偏向性影响，导致就业极化，即工业机器人使用会促进高技能岗位和低技能岗位增加，而常规型的中等技能岗位减少。不过，对于工业机器人应用引发的就业结构极化特征，实证证据存在较大争议。部分来自制造业上市公司微观数据的研究证实，机器人应用对不同技能劳动力就业的影响具有显著差异，将促进高、低技能劳动力的就业，导致中国劳动力就业结构整体上呈现两极化特征，不过部分研究也指出我国部分地区出现单极化趋势，甚至就业反极化效应。虽然对于是否出现就业两极化态势，实证研究结果并不一致，但一个相对一致的结论是，工业机器人的使用对就业技能结构带来一定冲击。
当前中国劳动力市场供需不匹配问题依然突出，“招工难”与“就业难”并存的矛盾凸显，而工业机器人的大规模应用势必进一步加剧这一就业结构性矛盾。如何有效化解人工智能给技术变革脆弱群体带来的不利影响，提升劳动力市场对技术冲击的适应性也成为当前加快发展新一代人工智能技术需要思考的问题。而劳动保护可能是解决这一问题的重要渠道和手段，在劳动保护程度较高的地区，企业往往受制于高昂的诉讼成本和严格的监管约束，不能随意调整或解除与员工签订的劳动合同，这可以在一定程度上缓解机器人应用对劳动力市场的负面冲击。但现有研究更多关注工业机器人应用的就业影响，而忽视了劳动保护对就业负面冲击的缓冲作用。
基于此，本文将工业机器人、劳动保护和就业技能结构纳入统一的分析框架，探讨人工智能技术下我国不同技能劳动力就业的影响趋势，试图厘清人工智能和就业技能结构在劳动保护政策体系下的互动关系，并针对性地提出相关劳动保护举措的优化路径和政策建议。研究致力于回答三个问题：工业机器人应用是否引发就业技能结构极化？劳动保护政策是否能够有效缓解工业机器人应用带来的就业负面冲击？不同劳动保护政策发挥的调节作用是否存在异质性影响？以上问题的探讨有助于捕捉劳动保护政策体系可能存在的发展困境，促进我国就业保护机制的长期改革，实现技术变革和劳动力市场的动态适配和协调发展。
本文可能的创新之处在于：一是现有研究大多利用国际机器人联合会(IFR)提供的分行业工业机器人数据，通过Bartik工具变量法计算各行业机器人渗透率与各省份分行业劳动就业份额，以此衡量省级层面的工业机器人应用水平。相较之下，本文直接采用工业机器人进口金额作为工业机器人应用水平的测度指标，能更精确衡量我国省级层面的机器人应用水平，有效地减少加权谬误，从而保证估计结果的可靠性。二是已有文献大多关注工业机器人应用的就业影响，随着人工智能技术不断推进，系统考察机器人技术对就业技能结构的影响以及劳动保护在其中的缓冲效应具有重要意义。有鉴于此，本文从工会力量和职业培训的双重视角出发，检验劳动保护是否缓解工业机器人对中低技能劳动力就业的负面冲击，扩展了机器人应用对劳动力市场就业领域的文献研究，丰富了后续的研究视角。三是研究发现，工业机器人应用导致我国就业技能结构呈现单极化趋势。工会力量和职业培训在“机器换人”与就业技能结构的调整中发挥了积极作用，证实了劳动保护在“机器换人”中的缓冲作用，为政府制定人工智能产业政策提供经验证据，对就业结构性矛盾下实现稳就业和更充分就业目标也有一定借鉴意义。
二、理论分析与研究假设
1. 机器人应用与就业技能结构
工业机器人应用对劳动力市场影响的分析框架大多基于任务制模型展开。任务制模型将工作任务划分为认知和手工任务、常规和非常规任务，并考虑机器人与劳动力在不同任务中的比较优势，以探究机器人对就业总量的影响。在任务制分析框架下，工业机器人应用对就业产生替代和补偿两种效应：一方面，工业机器人能够替代一系列重复性、程序化、标准化的常规任务，减少相应工作任务的劳动需求，产生就业替代效应。另一方面，工业机器人应用可扩大生产规模并创造新的工作任务，进而增加劳动需求，实现就业补偿。就业替代效应和补偿效应往往同时存在，而工业机器人应用对就业总量影响的净效应，要取决于两者之间相互冲减程度，在不同行业、地区和不同技能劳动者群体之间，这一影响存在较大的不确定性，研究结论也不尽一致，甚至相反。
然而工业机器人鲜明的偏向型技术进步特征，决定了工业机器人应用的益处不能均等地渗透到所有劳动人口，即机器人对劳动力市场影响会出现较大的群体异质性，因此人工智能的发展不可避免地会对就业结构产生冲击。特别是随着产业自动化和智能化技术改造的加速，工作岗位的生命周期不断缩短，落后岗位消失和新岗位创造的速度持续加快，使得技能成为技术变革中“任务—劳动力”匹配的关键。这种匹配互动关系不仅加速了不同技能劳动力在就业市场中的比重变动，而且加剧了就业极化趋势的形成，这一特征可能有两种态势：一是技能偏向型技术进步使高技能劳动力与技术进步的匹配度更高，因此信息化与自动化技术进步会降低中低技能劳动力的就业水平，并提升高技能劳动力的就业水平，形成就业单极化趋势；二是程序偏向型技术进步下常规型任务通常由中等技能劳动者执行，而对非常规型任务主要由低技能和高技能劳动者执行。工业机器人应用对执行常规型任务的中等技能劳动力产生挤出效应，而对执行非常规型任务的低技能和高技能劳动力产生互补效应，使得就业结构呈现两极化特征。
已有研究证实，工业机器人在中国的应用会产生就业结构极化现象，但由于对低技能劳动力就业的冲击存在较大不确定性，对于就业技能结构是两极化还是单极化，结论并不一致。如王永钦和董雯基于中国行业机器人应用数据和制造业上市公司微观数据，利用Bartik工具变量的因果关系识别策略分析机器人对劳动力市场的影响，研究结果表明，工业机器人应用在促进高技能和低技能劳动力就业的同时，挤压了部分中等技能劳动力的就业空间，进而引致就业技能结构的两极化。韩民春等学者利用2013—2017年中国286个地级市的面板数据，以科技活动人员和非科技活动人员作为高低技能劳动力的划分标准，估计工业机器人对制造业就业结构影响时发现，制造业的机器人应用对高技能劳动者的就业有促进作用，但对低技能劳动力产生了就业挤出效应，此时制造业就业结构朝着单极化方向调整。而孙早和侯玉琳利用2001—2015年全国省级面板数据，构造工业智能化指标，探讨工业智能化对就业结构的影响，研究指出，工业智能化替代初、高中学历劳动力的同时，增加了专科及以上与小学及以下学历劳动力的就业，使全国范围内劳动力就业结构整体呈现两极化态势。而在东南发达地区，工业智能化和过高生活成本的交互作用挤出低技能劳动力，致使发达地区就业结构呈现单极化趋势。
基于此，假设如下：
H1：工业机器人应用引发就业技能结构呈现极化特征。
2. 劳动保护在机器人应用对就业技能结构冲击中的缓冲作用
工业机器人应用规模增大给劳动力就业结构带来巨大改变，随之热议的“技术性失业”浪潮触发公众对失业风险的担忧，而加大劳动保护力度可能有助于劳动者应对技术性失业风险。其作用机制主要表现在三个方面。
第一，劳动保护直接提高了企业的雇佣成本和解雇成本，限制了企业调节劳动力数量的手段和力度，使得企业在进行雇佣决策时采取相对谨慎的态度，进一步规范企业用工行为，使得企业无法轻易调整或解除与员工签订的劳动合同，一定程度上缓解了机器人应用所带来的就业冲击；第二，为了维护劳动力市场的稳定性，政府很可能会通过强化劳动保护的方式来影响企业劳动力调整。如政府可以通过约谈或行业协会沟通等方式来劝诫企业有限度地裁减低技能劳动力，也可采取税收优惠、补贴等政策工具影响企业雇佣策略；第三，考虑到保障就业的需要，政府鼓励企业在发展智能机器人的同时，引导企业在利润分配时更注重工人利益，增强对职工的岗位保护和职业培训，在企业内部转移吸收“富余”劳动力，从企业内部缓解智能机器人应用对就业的不利影响。
进一步来看，劳动保护对异质性技能劳动群体所发挥的作用也存在差异。对于低技能劳动力而言，劳动保护导致雇佣弹性不足，可能会增加企业雇佣成本，进而刺激企业转而使用临时工、劳务派遣、共享员工等灵活用工形式来替代低技能用工，增加其就业难度，最终造成低技能劳动力利益减损。而对于高技能劳动力而言，劳动保护强化了企业对高技能劳动力的偏好。以职业培训为例，高技能劳动力具有较强的学习能力和适应能力，往往能更快理解新的思想理念、更易掌握新的技术手段以及适应新的工作环境。企业雇佣高技能劳动力可以节约更多的培训投资和调整成本以适应机器人等新一轮自动化技术。因此，高技能劳动力在职业培训的“加持”作用下，会赢得企业更强的偏好、拥有更好的就业机会。
基于此，做出如下假设：
H2：劳动保护有助于缓解工业机器人对就业技能结构的冲击，但对不同技能水平劳动力的影响有显著差异。
三、计量模型设定、变量选取与数据
1. 模型设定
本文采用2012—2020年中国省级层面数据，主要讨论工业机器人应用是否导致我国就业技能结构出现极化趋势，以及劳动保护能否有效缓解机器人应用对就业技能结构的负面冲击。基准计量模型设定如下：
=+++++
(1)
其中，下标和分别表示省份和年份。为被解释变量，衡量不同技能劳动力的就业水平。为关键解释变量，衡量工业机器人应用水平。代表一系列省级层面的控制变量，包含人力资本投资水平、产业结构、贸易开放度和城镇化水平。另外，模型还控制了省份固定效应和年份固定效应，表示随机扰动项。估计参数用以刻画工业机器人应用对不同技能劳动力就业水平的影响，若显著为正，则表明工业机器人应用会提高该类技能劳动力的就业水平；若显著为负，则表明工业机器人应用会降低该类技能劳动力的就业水平。
为探讨劳动保护是否有助于缓解工业机器人对就业技能结构的冲击，进一步明确劳动保护在其中的调节效应，构建劳动保护的调节效应模型如下：
=+++*++++
(2)
其中，表示省份在第年的劳动保护强度，采用工会力量和职业培训两种工具手段。
2. 变量选取和数据说明
(1)不同技能劳动力的就业水平()。人工智能技术的迅速发展，机器人等自动化设备的采用，在创造新岗位的同时摒弃落后岗位，使得不同技能劳动力的就业水平产生差异。本文采用受教育程度刻画其技能水平，将小学及以下的劳动力定义为低技能劳动力，初中和高中的劳动力定义为中等技能劳动力，大学专科及以上的劳动力定义为高技能劳动力。为反映不同受教育程度就业人员受工业机器人影响的程度差异，参照孙早和侯玉琳学者的做法，将不同受教育程度就业人员占总就业的比重作为被解释变量。
(2)工业机器人应用()。尽管近年来中国工业机器人本土企业市场占有率一路攀升，但与国际企业相比，我国工业机器人国产化率依然偏低。《中国工业机器人产业发展白皮书(2020)》显示，以世界机器人“四大家族”(ABB、库卡、发那科、安川)为代表的国外企业占中国工业机器人市场份额合计超过60%。考虑到省级层面数据的可得性与中国工业机器人市场的实际情况，本文采用工业机器人进口金额来衡量工业机器人的应用水平。根据《中国商品贸易数据库》基于HS2012八位编码体系提供的数据，我国工业机器人具体包括以下8类：喷涂机器人(84248920)、搬运机器人(84289040)、未列名的多功能工业机器人(84795010)、其他未列名的工业机器人(84795090)、集成电路工厂专用的自动搬运机器人(84864031)、电阻焊接机器人(85152120)、电弧焊接机器人(85153120)和激光焊接机器人(85158010)。本文使用的各省份年度工业机器人进口金额是上述八类机器人月度进口金额的加总。
(3)劳动保护()。本文将从工会力量()和职业培训()的角度来探讨劳动保护在机器人应用对就业技能结构冲击中的调节效应。工会力量方面，本文采用各省份工会集体劳动争议人数作为工会力量对应的代理变量。职业培训方面，本文采用各省份工会职业培训机构数来衡量职业培训的覆盖程度，并对其作对数化处理。
(4)控制变量。①人力资本投资水平()，采用各省份一般公共预算教育经费占一般公共预算支出的比重来测度；②产业结构()，采用各省份第三产业增加值占GDP的比重来衡量；③贸易开放度()，采用各省份进出口总额占GDP的比重来表示；④城镇化水平()，采用各省份年末城镇人口占总人口的比重来反映。
由于西藏、青海、港澳台地区工业机器人进口金额数据缺失严重，本文选取2012—2020年我国除此之外的29个省级行政单位的面板数据作为初始样本，并且根据数据处理原则，删除在样本期间数据缺失的部分。同时，为了保证研究数据真实可靠，本文未对缺失值作插补处理，因此本文实证研究有效样本量为225。除工业机器人应用数据外，本文其他变量的数据来源主要包括《中国统计年鉴》、《中国区域经济统计年鉴》、《中国劳动统计年鉴》和《中国人口和就业统计年鉴》。变量描述性统计见表1。
表1 变量的描述性统计
从表1可以看出，就业劳动力的技能以中等技能为主，占61.32%，而高技能劳动力和低技能劳动力份额近似，分别占20.33%和18.35%。全国历年机器人应用水平差异较大，最高达到7.05亿美元，最低为0，主要类型机器人年度平均进口金额为0.52亿美元。
四、实证结果及分析
1. 基准回归结果
基准回归根据公式(1)进行，检验工业机器人应用对就业技能结构的影响。面板估计模型的选择根据Hausman检验结果采用固定效应模型，具体回归结果见表2。结果显示，在控制省份固定效应和年份固定效应的条件下，工业机器人应用导致我国就业技能结构呈现单极化的趋势。具体来说，工业机器人应用促进了高技能劳动力的就业，同时挤出中等技能劳动力，但对低技能劳动力的影响并不显著，与已有的研究结论相似。
表2 基准回归结果
由于区域间要素禀赋和经济发展水平的不同，无论是工业机器人应用水平还是就业技能结构，都存在明显的区域异质性特点。因此，工业机器人应用对就业技能结构的影响可能存在区域差异。研究在基准回归的基础上按照三大经济区的分类标准将总样本划分为东部、中部和西部三组子样本进行回归，结果如表3所示。
从表3中可以看出，工业机器人应用对三大经济区的就业技能结构的影响存在明显差异。其中，东部地区的结果与全样本的估计结果基本一致，工业机器人使用显著提高了对大学专科及以上受教育程度劳动者的就业水平，同时降低了高中及以下受教育程度劳动者的就业水平，工业机器人对就业技能结构的影响呈现偏向高技能型单极化的趋势。西部地区的工业机器人应用在增加低技能劳动力就业的同时，挤压中高技能劳动力的就业空间，对中等技能劳动力产生更为明显的替代效应，就业技能结构呈现偏向低技能型单极化的趋势。而在中部地区，工业机器人应用促进高、低技能劳动力就业的同时，对中等技能劳动力形成了就业挤出效应，呈现两极化趋势，虽然作用效果并不显著(如图1所示)。
表3 工业机器人应用对劳动力技能结构的区域异质性影响
图1 工业机器人应用对劳动力技能结构影响的区域差异说明：横轴代表不同技能水平的劳动力，纵轴表示工业机器人应用变量系数，每条直线加粗部分为模型检验显著的部分。
上述异质效应可能归因于以下两点。
一是东部地区的人工智能基础技术和机器人发展规模位于全国前列。特别是长三角、珠三角和京津冀地区作为我国机器人产业园的集聚地，其建成与建设中的机器人产业园已达到46个，占全国机器人产业园总数的70.77%。而东部地区通过机器人产业体系的不断完备、国际机器人龙头企业的引进合作和核心技术的不断突破等途径，实现其在产业结构水平、集聚水平、规模效应、发展环境和创新能力上的优势，进而改变社会就业结构和企业内部劳动者的供求结构，东部地区机器人应用促进了高技能劳动力就业，但降低了中等技能劳动力的就业水平。因此，东部地区工业机器人应用对就业技能结构的影响呈现偏向高技能型单极化的趋势。
二是西部地区的产业发展具有劳动密集型和资本密集型特征，机器人应用提升了生产效率并降低了生产成本，使得企业扩大再生产，刺激了服务和手工劳动任务生产的产品需求，进而促进低技能劳动力就业。但西部地区机器人产业发展起步较晚，机器人产业发展环境受限，高技能劳动力流失现象仍然明显。中部地区作为老牌装备制造业基地，则侧重于雇佣生产关键基础零部件等重复性、标准化工作的低技能劳动力，机器人应用对中等技能劳动力的负面冲击和对高技能劳动力的正面创造效应尚未显现，所以中部地区机器人应用对就业技能结构的影响不显著。
2. 内生性处理
(1)IV-2SLS估计。工业机器人应用对就业技能结构的基准回归结果可能会受到内生性问题的干扰。一方面，就业技能结构影响企业投入工业机器人的生产决策。一般来说，新技术变革对人才需求旺盛的地区往往会产生剧烈影响，而后者也常常偏向于借助工业机器人等高端技术设备来实现工业再升级和智能化转型。另一方面，基准模型已纳入人力资本投资、产业结构、贸易开放等控制变量，但受限于数据可获得性，仍可能存在变量遗漏问题。
为克服模型双向因果和遗漏变量导致的内生性问题，本文选择滞后一期的劳动年龄人口比重作为工具变量，并使用二阶段最小二乘法(2SLS)进行估计。本文选择滞后一期的15—64岁的劳动力人口占总人口比重作为工具变量主要基于两点：一是我国劳动供给总量减少，人口结构老龄化趋势带来的劳动力供给紧缺可能倒逼企业进行技术创新，促使企业更多地使用工业机器人替代人力劳动工作，以弥补用工缺口所带来的负面影响。因此，滞后一期的劳动年龄人口比重与工业机器人应用之间必然存在强相关关系。二是劳动年龄人口比重会直接引致劳动供给结构的变化，但与就业技能结构并没有直接关联，即滞后一期的劳动年龄人口比重与就业技能结构正交，满足外生性的要求。
IV-2SLS估计结果如表4所示，第一阶段回归统计值为35.97，工具变量不存在弱工具变量问题。此外，第一阶段回归系数显著为负，说明劳动年龄人口比重越低则越倾向于采用自动化和智能化机器人。第二阶段估计结果表明，虽然基准回归可能低估工业机器人应用对就业技能结构的影响，但工业机器人应用导致就业技能结构呈现单极化特征的基本结论依然稳健。
表4 工业机器人应用对劳动力技能结构影响的IV-2SLS估计结果
(2)双重差分倾向得分匹配(PSM-DID)估计。工业机器人应用导致我国就业技能结构呈现单极化的特征，这一估计也可能存在样本自选择问题，即并非工业机器人应用导致高技能劳动力的就业水平提高和中等技能劳动力的就业水平降低，而是广泛应用工业机器人生产的省份本身就是具有高技能偏向型特征的省份。
为了克服不同省份对不同技能劳动力就业水平可能存在的系统性差异，本文采用PSM-DID的方法进行稳健性检验，并将《中国制造2025》这一国家行动纲领作为外生政策冲击。之所以如此，其原因在于《中国制造2025》明确提出了我国实施制造强国战略第一个十年的行动纲领，并将“高档数控机床和机器人”作为大力推动的十大重点领域之一，由此以工业机器人为代表的智能制造飞速发展。基于此，本文将2015年视为冲击发生的年份，同时借鉴杨光和侯钰的做法，设定2012—2020年人均工业机器人进口金额的平均数大于等于157(单位：万美元/百万人)的省份为实验组，否则视为对照组，计量模型如下：
=+*++++++
(3)
在使用双重差分法时，依据上述标准设置准自然实验分组虚拟变量，对实验组赋值为1，对控制组赋值为0。同时考虑到《中国制造2025》这一行动纲领的颁布时间，将准自然实验分期虚拟变量在2015年及以后的年份赋值为1，2015年之前赋值为0。因此，所有的样本可以分为四组：行动纲领颁布之前的控制组(=0，=0)、行动纲领颁布之后的控制组(=0，=1)、行动纲领颁布之前的实验组(=1，=0)和行动纲领颁布之后的实验组(=1，=1)。
要过滤样本自选择的影响，需要保证实验组和对照组之间不存在显著的系统性差异。本文采用logit模型来获得倾向得分，然后运用核匹配(Kernel Matching)对得分值进行匹配。而协变量检验结果显示，匹配后实验组和对照组各个协变量之间没有明显差异，说明采用PSM-DID方法进行稳健性检验是合适的。
PSM-DID的回归结果如表5所示，在剔除了样本的自选择效应后，工业机器人应用仍然形成就业技能结构呈现偏向高技能型单极化的特征，进一步支撑了前面回归分析的结论。
表5 工业机器人应用对劳动力技能结构影响的PSM-DID估计结果
3. 劳动保护的调节机制
工业机器人具有明显的技术进步偏向特征，可能加重中低技能劳动力的技术性失业风险，如前结论显示，工业机器人应用同时降低了低技能和中等技能劳动力的就业水平，而加强劳动保护力度可能有助于缓解这一冲击。因此，研究尝试将工业机器人、劳动保护与就业技能结构纳入统一的分析框架，检验劳动保护是否缓解工业机器人对中低技能劳动者就业的负面影响。具体估计采用公式(2)进行，对于劳动保护，我们选择工会组织和职业培训的覆盖程度作为其代理变量，依据如下：一是工会组织在维护劳动者权益和调解劳资纠纷方面发挥着重要作用。一方面，工会组织通过监督劳动合同的履行、订立包括防止企业强制性裁员条款的工会协议，提高了企业随意招募和恶意遣散员工的成本，使得企业在进行雇佣决策时采取相对谨慎的态度，进而规范了企业用工行为。另一方面，工会利用“集体发声”的方式，推进集体合同和集体协商，调解企业在机器换人过程中可能出现的劳动争议，一定程度上缓解了机器人应用所带来的就业冲击。二是机器人应用正在重塑工作所需要的技能，伴随机器人的大规模推广应用，职业所需的专业知识与技能将变化的越来越快。职业技能培训有助于劳动者实现技能提升，进而获得更为稳定的就业机会和可持续发展的能力，以实现广泛而有效的人机协作。
劳动保护调节效应的估计结果见表6。从中可以看出，工会力量和机器人应用的交互项系数同时对低技能和中等技能劳动力显著。更直观地，如图2(a)、(b)的交互项边际效应图所示，可以发现，随着工会力量的加强，工业机器人应用对低技能劳动者的替代效应逐渐增强，而对中等技能劳动者的替代效应逐渐减弱。
表6 工业机器人应用、劳动保护与劳动力技能结构的调节效应检验
图2 交互项边际效应的变化趋势
工会力量增强了机器人应用对低技能劳动者的替代效应。工会力量增强，可能在用工、解雇、工资福利等方面强化了劳动保护的作用，不仅直接增加了企业的用工成本，提高了企业使用劳动要素的价格，提高企业雇佣成本，还加大了企业调整雇佣规模和雇佣工资的难度。用工成本粘性降低了企业的劳动需求弹性，一定程度上增加了企业的雇佣调整成本。以工会力量为代表的劳动保护加强提高了企业的雇佣成本和调整成本，更加刺激企业采用机器资本来替代人工，而这种替代效应对低技能劳动力影响更大，因为低技能劳动力的解雇风险成本相对较低，对企业生产经营影响不大。
随着工会组织力量增强，工业机器人对中等技能劳动力的替代效应有所减缓。这可能是因为工会组织强化了劳动力在工作场所中的谈判力量及议价权，中等技能劳动力更愿意通过调解、仲裁等方式来处理“机器换人”所带来的劳资争议。此外，当企业采用机器人替代中等技能劳动力以执行常规型生产任务时，工会组织也会积极采取措施以应对中等技能劳动力的失业风险。工会通过集体协商的途径与企业管理层达成一致，利用企业内部劳动力流转机制，为劳动力提供在岗培训、转岗和再就业的机会，使中等技能劳动力在被替代和失业之间得到缓冲。此外，相对于低技能劳动力，中等技能劳动力仍是企业的核心竞争力，受到解雇的威胁相对较小。
职业培训的调节效应结果显示，职业培训的调节效应同时对低技能和高技能劳动力有显著影响，降低了低技能劳动力的就业水平，促进了高技能劳动力就业增加。而图2(c)、(d)的交互项边际效应进一步表明，随着职业培训覆盖程度的提高，工业机器人应用对低技能劳动力的替代效应和高技能劳动力的补偿效应逐渐增强。
职业培训增强了工业机器人应用对低技能劳动力的替代效应，不利于低技能劳动力就业，可能的原因在于：低技能劳动力人力资本积累不足，学习能力有限，需要花费较多的时间和精力接受系统的技能培训才能适应工业机器人的应用，引发“腐蚀效应”，降低了企业对低技能劳动力的有效需求。这进一步印证了新技术变革速度往往超过技能调整速度，被替代的低技能劳动力通过培训重新进入就业市场或转移到高技能岗位时面临较大的阻力，低技能劳动力工作能力与现有工作技能要求的结构性矛盾短期内难以消除。而随着职业培训覆盖程度的提高，工业机器人应用对高技能劳动力的边际效应逐渐增强，这一结果表明：接受职业培训促使高技能劳动力在劳动力市场中拥有更多的就业机会，而职业培训在“工业机器人应用—技能供求变迁—就业结构优化”的动态调整机制中发挥着积极作用，推动就业技能结构向高级演化。
五、结论与政策建议
基于2012—2020年中国省级面板数据，本文探讨了工业机器人应用对异质性技能劳动力就业的影响，并考察了劳动保护在“机器换人”中的缓冲作用。研究结果证实，工业机器人应用使得我国就业技能结构呈现出单极化特征，即工业机器人应用在促进高技能劳动力就业的同时，挤出中等技能劳动力，但对低技能劳动力就业的影响并不显著。在运用工具变量法和PSM-DID方法检验后，以上研究结果依然稳健。进一步的区域异质性分析表明，东部地区工业机器人应用对就业技能结构的影响呈现偏向高技能型单极化的趋势，西部地区的工业机器人应用在增加低技能劳动力就业的同时，挤压中高技能劳动力的就业空间，对中等技能劳动力产生更为明显的替代效应，使得就业技能结构呈现偏向低技能型单极化特征。而在中部地区，机器人应用对异质性技能劳动力就业的影响尚未显现，对就业技能结构的影响并不显著。
劳动保护的调节效应分析进一步显示，工会力量增强了机器人应用对低技能劳动力的替代效应，这可能与工会力量增强提高了企业的雇佣成本和调整成本，导致企业采取反雇佣策略，增加机器资本替代人工有关。与此同时，工会组织力量减缓了工业机器人对中等技能劳动力的替代效应，可能的解释是当企业采用机器人替代中等技能劳动力时，工会通过集体协商的途径与企业管理层达成一致，利用企业内部劳动力流转机制，为被替代的劳动者提供转岗培训和再就业机会，从而减缓了机器人应用对中等技能劳动力的负面冲击。而职业培训在强化工业机器人应用对高技能劳动力促进作用的同时，也加剧了工业机器人对低技能劳动力的替代作用。其原因可能是低技能劳动力人力资本积累不足，学习能力有限，需要花费较多的时间和精力接受系统的技能培训才能适应工业机器人的应用，而高技能劳动力具备更强的学习能力和环境适应能力。因此职业培训降低了企业对低技能劳动力的偏好，但赋能了高技能劳动力。
基于以上结论，本文得到如下政策启示。
第一，稳定重点群体就业，强化技能人才支撑。研究发现，工业机器人应用对中等技能劳动力有显著的替代效应，而对高技能劳动力有显著的正向影响。因此，为打好稳就业保就业攻坚战，实现就业质量的持续提升，政府要着力于完善现代化劳动保护体系，加强职业技能培训力度。完善企业工会体系，健全工会组织架构，畅通劳动者利益诉求渠道，充分发挥工会组织对劳动力市场的稳定作用。同时，要优化职业技能培训和转岗培训机制，建立终身学习体系，提高不同技能劳动者对新经济的适应能力，引导劳动者与工业机器人的互补式发展。尤其应关注中低技能劳动力在工作场所中的弱势处境，为中低技能劳动群体提供有效而持续的培训机制，着力推动我国就业技能结构整体升级。
第二，充分考虑地区异质性特征，制定适配产业扶持政策。工业机器人应用对就业技能结构的影响存在显著的区域异质性。其中东部地区的工业机器人应用对就业技能结构的影响呈现偏向高技能型的单极化趋势，而西部地区工业机器人的发展对当地的劳动者技能结构反而呈现偏向低技能型单极化特征，这与各区域的要素禀赋密不可分。因此，各区域应根据当地的发展阶段和就业现状，推出因地制宜、全面协调的机器人发展战略，积极探索与之适配的产业转型升级路径，不应盲目进行资本扩张和推广机器人使用，而应注重提高资本、技术与劳动力的融合水平。
(致谢：感谢广东制造业大数据创新研究中心、工业大数据战略决策实验室对本研究的支持。)
