进口贸易、技术溢出的环境效应与工业大气污染——基于省级面板数据的中介路径分析

丁日佳 ， 张诗悦

（中国矿业大学（北京） 管理学院，北京 100083）

改革开放以来，中国工业化发展取得了举世瞩目的成绩。然而，中国以重工业为主的“先污染、后治理”和“高污染、高耗能、低质量、低产出”的粗放型和外延式发展模式对生态、资源和环境产生了不良的影响，尤其是生产和加工过程中排放的大量烟尘、粉尘、硫氧化物、氮氧化物等有害气体给人们身体健康和日常生活带来了极大的危害[1]。如果不采取合理的治理措施来遏制环境污染，则大气质量可能无法跨越改善的拐点[2]。人们的健康福祉也可能随之无法得到有效保障。因此，工业大气污染的治理势在必行，而技术创新和技术进步是缓解大气污染的根本推动力。“十二五”和“十三五”期间，各级政府积极提倡工业领域加强技术改造升级，研发节能环保的新技术、新工艺和新装备，发挥科技创新在工业绿色发展中的引领作用[3]。大气污染的治理技术不仅来源于本土的自主研发，而且依赖于国外先进的清洁生产技术。在全球价值链分工和生产网络体系中，拥有绝对优势研发资本的发达国家占据国际技术创新的制高点，掌握保障“蓝天白云”的绿色科技。国务院发布的《大气污染防治行动计划》中提出要“加强大气污染治理先进技术、管理经验等方面的国际交流与合作。”国务院发布的《“十三五”生态环境保护规划》中提出要“促进生态环保理念、管理制度政策、环保产业技术等方面的国际交流合作。”《工业绿色发展规划（2016—2020年）》指出要“强化绿色科技国际合作”。可见，国外先进的绿色技术溢出是治理工业大气污染不可忽略的途径。

目前，就多数发展中国家而言，FDI技术溢出和国际贸易技术溢出成为全球间技术转移和技术扩散的主要途径[4]。相比之下，囿于知识产权保护、外部环境和资金条件等因素限制，利用进口贸易来引进和吸收国外先进技术无疑是更为直接和有效的方式[5]。据海关统计，1998年和2019年中国进口贸易总值分别约为1.16万亿元和14.31万亿元，二十年增长了约12.34倍。在进口商品结构中，中国的环境产品技术创新能力较弱。一些环保、节能的生产工艺使用装备和环境监测设备等，例如：大气、水和土壤等监测设备和核心技术，基本依赖于进口；一些技术含量较高的中间品，例如太阳能级的多晶硅、汽车发动机和计算机芯片等也需要大量进口[6]。不可否认，进口贸易技术溢出对中国大气环境质量改善产生了深远影响。然而，受国际贸易市场和进口国内部发展水平的影响，进口环境产品所带来的技术溢出是否一直促减工业大气污染排放？如何有效利用进口贸易渠道带来的技术材料和研发资本治理工业大气污染？基于此，本文通过探索技术溢出的环境效应内容，从理论解析和实证分析双重角度，验证进口贸易技术溢出对工业大气污染排放强度的内在机理和三类中介效应，以期对中国工业大气污染治理和绿色发展有所启发。

一、文献综述

进口贸易技术溢出不仅关系到全要素生产率和经济增长，还会影响到产业结构、能源效率等，同时其影响方向存在“抑制论”“促进论”以及“双刃剑论”三种截然不同的观点。一般情况下，物化在进口贸易产品的有形知识溢出和非物化型无形知识转移对技术接收国的经济增长和技术进步产生一定的外部性作用[7]。Coe与Helpman最早运用回归法验证了国外R&D溢出效应(CH模型)，随后证明了工业化国家的R&D资本存量每增加1%，发展中国家的TFP就会增加0.06%[8-9]。此后，投入密度[10]、私人品和公共品[11]、人口数量[12]、资本品和中间品贸易总量[13]等概念被引入，反复修正CH模型的不足之处，且均证实了发达国家的进口贸易技术溢出对发展中国家的TFP和经济增长产生了正效应。然而，受本国贸易结构、吸收能力、区域差异、国外知识产权保护和技术垄断等因素影响，进口贸易技术溢出对TFP和经济增长产生不确定性作用，例如全国范围内，进口贸易技术溢出刺激经济增长和全要素生产效率的提高，然而东、中、西部地区因其异质性区域条件限制和吸收能力程度，未必都能获取积极的技术溢出效应[14-15]；从进口产品的传导途径来看，进口资本品主要包括成套的机器设备和运输设备，由于其技术复杂性和国外知产权保护，可能难以汲取先进的技术[16]。此外，进口贸易渠道的技术溢出红利也蔓延到产业结构[17]、能源效率[18-19]、创新研发[20]等要素。制造业结构优化升级、全要素能源效率提升、科技创新产出都与进口贸易技术溢出的积极效应息息相关。

贸易自由化和经济全球化使得环境污染成为国际性问题。国际贸易对生态环境的影响错综复杂，众多学者辨析两者之间是存在“污染光环”假说，还是存在“污染避难所”和“向底线赛跑”假说？以及国际贸易规模效应是占据主导地位加重了环境污染，还是技术效应和结构效应抵消了规模效应进而提高了环境质量？进一步，延伸至进口贸易技术溢出与大气污染的关系的相关探索却很少，目前主要聚焦于国际贸易与大气污染、国际贸易技术溢出与温室气体两个层面的机理内容。一方面，最为经典的国际贸易与大气污染研究是Grossman和Krueger[21]提出的国际贸易环境效应，论证了经济增长对二氧化硫等污染排放物的影响。随后，研究领域集中在贸易自由化如何通过收入效应、技术效应、规模效应、结构效应和规制效应等途径对空气质量产生作用，并测算了贸易产品所隐含的大气污染[22-23]。另一方面，对外贸易技术溢出不一定具有节能减排效应。虽然国际前沿的环保技术和绿色技术会通过国际贸易渠道转移至技术落后国家，在一定程度上降低了温室气体排放[24]，但是进出口市场价格与产品质量迫使企业面临利润压缩和改进能效成本增大的双重竞争压力，从而减少进口贸易技术溢出的环境福利[25]。此外，进口国的内资企业实际需求状况[26]、环境政策和规制水平[27]、模仿和吸收能力特征[28]，也导致技术溢出对温室气体的减排作用具有不确定性。

现有文献表明，进口贸易产品可以为一国带来显著的技术溢出效应，从而进口贸易技术溢出对其经济规模发展、科技创新、要素结构调整和能源效率改善等产生了重要影响，同时国际贸易一般通过对经济活动产生规模效应、结构效应和技术效应来影响大气污染。然而，鲜有文献研究进口贸易技术溢出在大气污染治理过程中的作用与地位。基于此，本文边际贡献在于：（1）研究视角上，首创性探索资本品和中间品的进口贸易技术溢出与工业大气污染排放强度的关系；（2）研究方法上，以工业经济规模扩张效应、要素禀赋结构优化效应和工业技术进步效应作为中介切入点，运用广义空间两阶段最小二乘法(GS2SLS)和逐步回归法具体探析进口贸易技术溢出对工业大气污染排放的传导路径，以及不同程度技术溢出的环境效应如何作用于污染减排。

二、理论分析与研究假设

根植于众多学者的研究成果，进口贸易结构中的资本品和中间品蕴含着较高的技术含量和技术优势，技术溢出效应显著。再结合国际贸易环境效应理念、技术溢出的外部性和相关理论分析，针对工业领域提出了进口贸易技术溢出的环境效应内涵，认为国外先进的R&D活动通过扩大工业经济规模、优化工业结构和提升工业技术进步分别带来经济规模扩张效应、要素禀赋结构优化效应、工业技术进步效应，对工业大气污染排放产生异质性影响趋势。具体分析思路如图1所示。

图1 进口贸易技术溢出对工业大气污染的中介路径图

（一）经济规模扩张效应

内生增长理论认为，技术进步是内生的，是经济增长的内在动力，使得经济能够不依赖于外力推动实现持续增长。国外R&D知识通过国际贸易技术溢出渠道间接影响进口国的技术进步，进而促进经济发展，也可作为国外技术溢出直接影响其经济水平[29]。这就意味着，进口贸易技术溢出有助于推动工业经济增长，强化工业企业经济规模效应。然而，经济规模的扩大是否起到节能减排作用尚未定论。一方面，就供给侧而言，生产规模的扩张可能带来规模经济、提高生产效率、降低生产成本，以较少的投入换取更多的产出，加速生产分工趋于专业化和合理化，生产集约效应优势明显，进而节约能源和资源的消耗，缓解大气污染；或者带来了规模不经济，在产业结构和污染排放系数不变的情况下，生产要素投入过多，导致产出成本有所增加，加大了对能源和自然资源的消耗力度，加重了空气污染物的排放。另一方面，就需求侧而言，经济发展初期既是工业规模扩充跃进阶段，也是污染型产品大量需求的时期。较低的经济发展水平促使消费者更多倾向于温饱型物质需求。为迎合市场发展趋势和人民生活基本需要，伴随工业经济集聚规模扩张、大幅度提升非清洁型生产要素的投入，如果此时技术进步和政府监管无法有效遏制工业规模扩张带来的环境污染，则技术溢出导致的规模效应会促进大气污染增加。后期随着经济质量的提升，当满足物质需求时，越来越多的居民偏好于绿色消费，更愿意选取技术含量高和环境质量检测严格的产品。市场的绿色刚需也使得规模扩张后的工业企业投入新技术和新工艺，加大对资源节约型和环境友好型工业品的生产和加工，从而在一定程度上降低大气污染的排放。由此，本文提出如下假设:

假设1a.进口贸易技术溢出通过扩大工业企业经济规模加剧了工业大气污染排放；

假设1b.进口贸易技术溢出通过扩大工业企业经济规模降低了工业大气污染排放。

（二）要素禀赋结构优化效应

钱纳里工业化阶段理论认为，后发国家的产业发展遵循劳动密集型产业转向资本密集型产业再转至技术密集型产业的变动规律。其中，技术进步是推动整个产业结构变动，迈向高科技、高附加值和绿色环保领域的重要助推器[30]。现阶段，中国是以劳动密集型行业和非技术密集型行业为主的工业结构，由于严重依赖于非清洁能源，存在落后的行业产能和结构分布不均衡现象，导致工业生产运营过程中产生了大量污染废弃物[31]。进口贸易技术溢出作为国外先进技术的扩散形式，在一定程度上能够促进工业生产结构从劳动密集型、资源密集型和资本密集型转变为技术密集型，促使工业生产可能性前沿外移，降低了高耗能、高污染行业份额，使其向绿色化、集约化和清洁化状态转换。（1）优化高碳密集型和高能耗密集型的行业结构。目前诸多发达国家积极开发和使用可再生能源和新兴能源，致力于研发绿色低碳产品。将进口环保材料、治污设备、环境监测设备等带来的新能源技术和低碳技术投入工业生产使用，有利于提高能源利用效率，促使高碳行业转向低碳行业，逐步优化石油加工炼焦、有色金属冶炼、电力热力等高耗能行业结构，降低高污染行业的密集度。（2）提升技术密集型行业规模。发展中国家往往热衷于进口发达国家的技术密集型产品，例如：机械设备、计算机芯片、航空发动机、医疗器械等。这些进口商品一般具有高附加值、高技术化和高集约化特征，带来了高新技术的溢出和普及，促使产业重心向技术密集型行业转移。而技术密集型行业具有低污染和低排放性质，其行业比例的增加无疑有助于减少大气污染。（3）升级原有工业结构，促进新兴环保产业的诞生。原有的产业部门依托进口贸易技术溢出的要素禀赋，通过诸如改造原有的环保生产技术和产品工艺、摒弃高污染型结构、聚焦国际创新技术，以较少的资源投入，获取清洁型高产出，促使传统工业部门与时俱进、绿色发展。进口贸易带来的节能环保技术和清洁生产技术等也会掀起新兴环保产业的浪潮。高端的节能减排技术溢出会促使工业领域不断探索降污排污的生产方式和管理运作理念。在充分掌握国外高精尖设备和零部件等核心技术基础上，追踪和加强自身绿色环保理念和原创设计，同时伴随市场环保产品需求大幅度提升，促使工业领域新兴环保产业的诞生，进而起到节能减排作用，降低大气污染。由此，本文提出如下假设:

假设2.进口贸易技术溢出通过调整和优化工业结构来降低工业大气污染排放。

（三）工业技术进步效应

国际贸易技术溢出的途径可分为静态途径和动态途径[32]。静态途径是指通过进口资本品和增强中间品质量和种类直接影响了技术转移。动态途径是指国际贸易通过示范与模仿效应、人员培训效应和关联效应等间接影响技术溢出。联合静态效应和动态效应，进口贸易技术溢出发挥的技术进步效应能够在一定程度上降低工业大气污染。从静态角度来看，中间品种类的增加能够带来多样化的技术溢出[33]，使得最终产品部门从中选择更加符合绿色要求的产品，增强最终产品的生产效率。中间品投入质量的提升与其所带来的生产力成正比[34]，高质量的中间品往往具备高效清洁型技术性质，因此中间品投入质量越高就会带来更多的绿色生产力和环保技术溢出。另外，资本品内含有技术知识，直接投入使用或投入生产过程的进口治污机器设备和清洁型运输设备可有效提高生产效率和能源效率，改造和升级企业现有的生产技术，增加更多的绿色技术溢出。从动态角度来看，在“干中学”效应推动下，针对进口贸易带来的清洁型技术溢出，工业企业通过逆向工程来获取环保技术的构建原理，通过培训优质的人力资源来获取绿色技术服务和使用方法等隐性知识[35]，促进工业企业的技术进步，推动节能减排。此外，进口商在实施贸易交易活动时与产业链的上游和中下游企业发生前向和后向关联效应带动了产业间的绿色技术溢出。通过产业链溢出渠道，国内工业供应商等也会获得更为直接的先进技术，最终降低大气污染。

然而，长期的进口贸易技术溢出可能会产生环境回弹效应和抑制中国自主创新现象，未必一直促减空气污染。一方面，技术进步作为提高能源效率的重要途径并不是总以降低污染排放为目标。杰文斯悖论（Jevons Paradox）详细剖析了环境回弹效应的发生机理，并指出技术进步会带来经济规模的扩充，刺激经济增长从而导致更多能源消费来替代劳动力和生产活动中的资本投资。同时，技术进步带来的能源效率的提升也降低了单位生产成本和能源价格，为追求利润最大化，企业必然增加更大范围的能源消费以促进产出[36]。因此，能源消费量的增加也伴随着环境污染再次反弹，恶化了空气质量[37]。另一方面，技术溢出同时对本土自主创新具有技术互补效应和技术替代效应[38]。受技术差距、专利垄断、吸收能力等因素制约，大量进口贸易技术溢出的涌进冲击了进口国市场份额，市场竞争带来的高额研发成本和技术挤出效应只能使少数内资企业真正掌握国际前沿的节能减排核心技术，很容易形成“技术锁定”困局。长期循环的贸易模式可能会引发国内企业的研发惰性，削弱了工业企业开发清洁型能源技术的绿色竞争优势，阻碍其自主创新步伐，无法从根本上解决中国工业大气污染问题。由此，本文提出如下假设:假设3a.短期而言，进口贸易技术溢出通过增强进口国工业企业技术进步来降低工业大气污染排放；假设3b.长期而言，进口贸易技术溢出可能会形成回弹效应和阻滞自主创新局面，从而增加了工业大气污染排放。

三、模型设定、指标测算和数据来源

（一）模型设定

1.空间自相关检验

大气污染具有物理流动性和空间扩散性特征，且区域间的经济网络发展联系紧密，有必要考虑空间溢出效应。空间相关性检验是构建空间计量模型的重要前提。如果通过了该检验，则证实工业大气污染排放强度存在空间依赖性，可以构建空间计量模型。Moran’I检验是目前应用最为广泛的空间相关性检验方法。全局Moran’I统计量取决于反映观测区域间的地理关系强度的空间权重矩阵，表明空间邻接或空间邻近的区域单元属性值的相似程度。全局Moran’I值用如下公式计算

其中，xi和xj分别为第i个地区和第j地区工业大气污染排放强度的水平；x为xi或xj的平均值；n为地区总数；wij为空间权重矩阵中W第i行第j列的元素，采用省域i和j间的地理距离空间权重矩阵表示，具体如式（2）所示。Moran’I统计量取值一般在—1～1之间，大于0表示工业大气污染排放强度空间正相关，小于0表示工业大气污染排放强度空间负相关，等于0表示区域间工业大气污染排放强度不相关。越接近1表示区域间工业大气污染排放强度差异较小或空间集聚性较强，容易出现大气污染排放强度的区域高值集聚和低值集聚；越接近—1，则反之；接近0，则表示工业大气污染排放强度空间不相关。

2.空间中介效应模型设定

如果上述的莫兰指数通过显著性检验，则表明可进一步构建空间计量模型。多数文献采用最大似然估计(MLE)来构建空间模型，但是前提条件是必须满足独立同分布与正态分布的假定才可达到最有效率的估计。相比之下，广义空间两阶段最小二乘法(GS2SLS)放宽了假设限定，即使存在异方差与非正态分布，依然保持一致估计[39]。同时，依据上述理论分析可知，受不同的中介路径影响，进口贸易技术溢出可能在不同的阶段对工业大气污染排放存在异质性作用，由此引入进口贸易技术溢出的二次项和三次项来分析与工业大气污染排放的非线性关系更为合理。因此，本文运用GS2SLS方法先检验两者是否存在非线性关系，并选择一阶空间滞后项作为工具变量。构建如下基准模型

为消除异方差和保证数据平稳性，将各变量加以对数化表示。其中，i和t为各省市和年份；SO2为工业大气污染排放强度；W为空间权重矩阵；importit为进口贸易技术溢出，分别由资本品技术溢出(cap)和中间品技术溢出(inte)构成；γ、α、β为估计参数；εit表示误差项。

进一步，基于理论分析和研究假设、中介效应思想和逐步回归法[40]，可以有效阐述进口贸易技术溢出通过何种路径影响工业大气污染的排放强度。为有效验证各变量间的非线性关系，借鉴Gan等[41]研究思路：第一步以上述基准回归式（3）检验技术溢出对工业大气污染排放强度是否显著，即α1、α2和α3的显著性；如果显著，则进行第二步，即分别以经济规模、要素禀赋结构和技术进步三类环境效应为因变量，进口贸易技术溢出为解释变量，检验技术溢出对中介变量的回归系数α11、α21和α31是否显著，如式（4）所示；第三步，再次以工业大气污染排放强度作为被解释变量，中介变量加入基准回归式（3），检验三类中间机制和技术溢出对排放强度的影响，如式（5）所示。如果式（3）～式（5）的技术溢出和中介变量的回归系数全部显著，则表明存在部分中介效应；如果式（5）中的技术溢出回归系数不显著，则表明存在完全中介效应

其中，mediateij表示中介变量，包括工业经济规模（scale）、要素禀赋结构（struc）和工业技术进步（tech）。

（二）指标测算

1.核心解释变量

资本品(cap)和中间品(inte)进口贸易技术溢出构成核心解释变量。联合国BEC分类法(广义分类法)根据商品的最终用途，把国际贸易结构分成资本品、中间品和消费品三部分。相比之下，资本品和中间品的技术含量较高，技术溢出效应更为明显。众多研究倾向于采用资本品和中间品技术溢出来探析对TFP的影响。为此，本文利用LP的方法测算两种贸易技术溢出[10]，具体如下所示

其中，t为时间；i代表30个省市；j表示9个进口国家。中国的进口贸易伙伴国主要集中于OCED国家，尤其与G7国家贸易往来密切。发达国家是国际先进技术的领先者，工业化国家的主要研发活动也集中在G7国家。结合各省市进口国别的贸易结构现状，最终选取进口贸易伙伴为日本、新加坡、韩国、英国、德国、法国、意大利、加拿大和美国[14][16][28]。

Capitalijt和Intermediateijt分别表示在t时期、i省市从j国进口的资本品和中间品贸易额，采用BEC分类法中的代码41和521作为资本品，代码42和代码53表示中间品[16]。由于缺乏各省市的中间品和资本品进口贸易额数据，因此本文首先计算每年各省市在j国的进口贸易总额占全国在j国的进口贸易总额的比例；然后，利用这个比例分别乘以全国在j国进口的资本品和中间品贸易总额。 此外，GDPjt表示t时期j国的产出水平，利用各国的GDP平减指数将现价GDP折算为2000年不变价格；RDjt表示t时期j国的R&D资本存量。各国的研发支出存量参考Griliches的永续盘存法[42]。其中，研发资本的折旧率，一般设定为5%[8]。

为便于下文实证分析，承接上述测算，表1刻画了30个省市统计期内的资本品和中间品进口贸易技术溢出的均值。从表1可知，内陆经济较为发达或位于沿海地带的区域，如北京、天津、上海、浙江、广东等，均超过100亿元的技术溢出，而经济发展较为落后，位于中、西部地区的进口贸易技术溢出较低，表明区域间因贸易结构和经济水平的差异导致不同程度的贸易技术溢出。

表1 各省市资本品和中间品进口贸易技术溢出的均值 单位：亿元

2.被解释变量

工业大气污染物排放主要是指工业废气排放。工业废气又包括二氧化硫、烟尘、粉尘和氮氧化物等。在衡量工业大气污染时，相关研究直接选取工业废气排放量或具有代表性的工业二氧化硫作为其代理变量，或者间接构造污染指数和污染强度加以综合度量。鉴于不同的大气污染指标选取会影响到研究问题本质和实证分析结果的可靠性，同时较为完整和连续的工业大气污染数据主要以二氧化硫和烟尘粉尘为主，而工业二氧化硫对空气和降水等环境监测指标影响巨大，也是困扰居民身体健康的重要大气污染物[43]。因此，在遵循以往学者的测度方式基础上，本文采取工业二氧化硫排放强度(SO2)作为工业大气污染的代理变量。具体算法为工业大气污染等于工业二氧化硫排放量除以实际GDP总量。

3.中介变量

工业经济规模(scale)：众多学者衡量企业经济规模主要从投入参数和产出参数两个方面入手。一般情况下，人力资本和物质资本等作为企业规模投入要素。产品产量、产值、营业收入和利润总额等可表示企业规模产出要素[44]。鉴于数据的连贯性和一致性，本文利用人均规模以上工业企业的工业总产值来表示经济规模效应，并利用各地区工业产品出厂价格指数将其折算为2000年不变价格。

要素禀赋结构(struc)：依据行业性质和结构，工业内各类行业可划分为技术密集型和非技术密集型。通常，资源型密集行业、资本密集型行业和劳动密集型行业构成非技术密集型工业结构。由于本文年份跨越维度较长，为保证各年间的行业统计口径一致性，遵循《国民经济行业分类标准GB/T 4754—2019》文件内容，对工业的27个行业按要素密集度进行分类。各行业的工业销售产值能够体现报告期内工业企业销售的工业产品总量，有效反映行业的构成要素和产出规模；同时也能够保证数据的连贯性和可获取性，因此选取该指标作为衡量标准。此外，通过理论分析可知，进口贸易技术溢出既可以优化和调整非技术密集型行业，也可以提升技术密集型行业份额。为对比技术溢出优化和改善结构效应的程度，借鉴产业结构高级化指标测算方式，最终运用各省市的技术密集型行业的工业销售产值占非技术密集型行业的工业销售产值的比例来测量工业结构的高级化水平。

工业技术进步(tech)：国内技术进步一般来源于内部的自主创新和技术引进，本文重点实证外部的技术溢出作用于工业内部技术进步。测算技术进步的研究方法主要包括创新要素的投入、创新要素的产出和集成两者的数据包络分析法。综合来看，以人力资本、R&D经费等创新投入可能会出现资源利用效率较低的现实情况而无法准确衡量技术进步[45]。以DEA为首的投入产出算法虽可精准地诠释全要素生产率的内涵，反映技术溢出效应，也是众多学者热衷的测量方法[14-16]，但是目前难以获取口径一致的工业劳动力全时当量和工业固定资产存量数据。从产出角度看，发明专利数是被授权和认可的技术创新，更能真实地反映工业企业技术进步和“干中学”效应。因此，本文采用发明专利数作为工业技术进步效应的代理变量。

4.控制变量

防止因遗漏变量造成的样本误差，参照相关理论与实证分析，选定如下变量来控制模型：能源结构(es)：工业生产运作和加工中的大量化石能源消耗是造成工业大气污染排放的主要来源，本文采用煤炭消费占能源消费总量的比例加以表示；环境规制强度(er)：政府适度和持续的环境监管和约束，有助于实现“波特假说”，提升污染防控治理水平，本文利用工业污染治理完成投资额占工业增加值的比例进行考察；城市化水平(ubran)：城市化发展带来了工业化水平提升，城市化进程以粗放式或集约型模式推进对污染排放至关重要，本文运用人均铺装道路面积为城市化进程的替代变量；劳动力成本(labor)：劳动力价格的高低突显着企业生产成本和劳动力素质程度，是诱发产业流动与转换，以及环境污染的重要渠道[46-47]，本文采用城镇单位就业人员平均工资衡量。

（三）数据来源

本文选取中国2000—2017年30个省市（除西藏及港、澳、台地区）作为研究区域。2000—2015年的工业二氧化硫数据来源于《中国统计年鉴》，2016年和2017年的数据来源于《中国城市统计年鉴》和各省市统计年鉴，部分缺失数据利用ARIAM模型进行预测和估算；关于进口贸易技术溢出数据、公式中各类数据主要来源于世界银行数据库、OECD数据库、UN Comtrade数据库、中国各省市统计年鉴(2001—2018年)；中介变量数据主要来源于《中国工业统计年鉴》和国家统计局网站，部分缺失数据利用插值法加以补充；其他控制变量数据来源于《中国能源统计年鉴》《中国统计年鉴》《中国环境年鉴》。本文利用stata15.1进行模型验证。

四、实证分析

各模型回归结果显示，R2值均在0.9之上表明模型拟合优度较好。工业大气污染排放强度的莫兰指数和空间滞后项均在1%水平下显著，验证了各区域的大气污染排放强度具有明显的空间集聚效应和溢出效应。由于测算的拐点值较多，为确保拐点值的精准度和后续结果分析的准确性，将回归系数保留至五位小数。

（一）空间自相关检验

诸多研究表明，工业污染物具有显著的空间依赖性。因此，本文根据式（1）和式（2）采取空间自相关指标，即莫兰指数来检验工业大气污染排放强度是否具有空间相关性。由表2可知，工业二氧化硫排放强度的莫兰指数在0.048～0.102区间变化，P值在0～0.01间变动，全部通过1%显著性水平检验，说明工业二氧化硫排放强度的空间自相关性和集聚性较强，省市间的污染排放强度存在正相关性。

表2 工业二氧化硫排放强度的莫兰指数

（二）基准回归结果

表3对应了式（1）的基准回归结果。表3中模型1和模型2显示资本品和中间品技术溢出的三次项系数全部在1%水平下显著，验证了进口贸易技术溢出与工业大气污染排放具有“正N形”关系，即进口贸易技术溢出在不同阶段促进或抑制了工业二氧化硫。“正N形”曲线是否存在拐点，还需通过进一步计算方可得知。以模型1的资本品技术溢出的拐点测算为例，具体的拐点求值过程如下所示

表3 基准回归结果

原回归方程见式（7），可以通过求导数的形式对该方程求出极值点，对式（7）中的ln cap进行求导，计算结果如下

在控制变量方面，以模型1为例，煤炭燃烧为主的能源结构、城市化发展、环境规制均加剧了工业大气污染排放，影响系数分别为0.065 20、0.334 28和0.045 38，并在10%、1%、10%水平下显著，而劳动力成本的提升减缓了工业二氧化硫污染，影响系数为0.057 11，但未通过显著性水平检验。原因可能为煤炭是工业生产和运作的主要基础原材料，煤炭的燃烧和加工释放了大量二氧化硫；交通密度和城市化水平的提升推动了工业化发展，城市化和工业化建设过程中对能源和资源的消费需求加重了工业污染的排放强度；即使工业污染治理用于废气投资的力度和比例逐年上升，但是环境规制并未有效改善工业大气污染排放强度。一方面可能是治理投资成效往往具有较长滞后性，在短时间内难以显现规制成果；另一方面废气治理投资额占工业总产值比重不足，无法满足经济社会发展对环境污染治理投资的客观要求；劳动力成本和质量越高代表该区域的经济水平和基础设施越完善，存在较高素质的人力资本水平，相应地环保意识越强，减污排污设备越齐全。受篇幅所限，表4和表5的控制变量回归结果均做省略。

表4 进口贸易技术溢出对中介环境效应影响的回归结果

表5 直接效应和中介效应的回归结果

（三）中介效应回归结果

基于上述的显著性水平验证，可进一步检验进口贸易技术溢出对工业大气污染排放的直接效应和中介效应。表4和表5对应式（3）和式（4）。由表4可知，控制不同的中介变量影响效应后，进口贸易技术溢出与工业大气污染仍呈“正N形”特征，直接效应依旧显著，说明存在部分中介效应，分析内容如下所示。

首先，资本品和中间品贸易技术溢出为全国多数省市带来了工业规模扩张效应。虽然模型3和模型6反映资本品和中间品技术溢出与工业经济规模具有“倒N形”特征，但是纵观样本期间内技术溢出的均值（如表1所示），仅有北京、天津、辽宁、上海、江苏、浙江、山东和广东八个地区的资本品技术溢出，外加福建九个省市的中间品技术溢出超越第二个拐点，其余省市均在两个拐点之间。因此，大约73%省市的进口贸易技术溢出扩充了该区域工业企业的经济规模。而较为发达的经济地带受技术溢出的规模扩张效应影响并不明显，可能由于其经济发展迅速，高度的经济集聚性已成饱和状态，同时这些省市多数位于沿海地区，外商直接投资业务往来频繁，更多地受到FDI技术溢出带来的经济作用。进一步，模型9和模型12的工业经济规模回归系数为正，在1%水平下显著，这说明在工业结构和污染排放系数不变的情况下，工业企业非清洁型生产要素投入过多，产出成本有所增加，加大了对能源和自然资源的消耗力度。同时市场绿色刚需和居民绿色观念尚未成熟，更多注重污染型产品消费，未能倒逼工业企业生产环境友好型工业品，加重了空气污染物的排放。综合来看，目前全国约73%省市的进口贸易技术溢出通过扩张工业经济规模提升了工业大气污染排放力度，在一定程度上验证了假设1a的合理性。

其次，资本品和中间品贸易技术溢出起到了优化和改善要素禀赋结构的作用。经实证检测，技术溢出的三次项和二次项均不显著，只有一次项系数在1%水平下显著，模型4和模型7的变量影响系数分别为0.120 39和0.191 98，这表明贸易技术溢出积极推动了全国工业结构高级化发展，提升技术密集型行业份额，增强节能环保生产性服务业发展，降低高碳高耗能型行业比例，淘汰过剩产能，促使非技术密集型行业转型至资本知识密集型行业。与此同时，模型10和模型13展示要素禀赋结构的回归系数为−0.097 61和−0.112 82，表明高级化的工业结构提升了行业间的协同发展程度，推进高污染行业治理力度和升级改造，推动绿色环保产业市场化、集约化和专业化发展，有效促进二氧化硫减排效应。整体来看，全国30个省市的进口贸易技术溢出通过调整和优化工业结构降低了工业二氧化硫排放强度，验证了假设2的合理性。

再次，资本品和中间品贸易技术溢出能够带来工业技术进步效应。经实证检验，虽然贸易技术溢出与工业技术进步呈“倒N形”非线性关系，但是同样对比表1中的统计期内的技术溢出均值，模型5和模型8展现30个省市全部超过第一个拐点，但未到达第二个拐点，验证了样本期间内全国的贸易技术溢出对工业技术进步产生了积极效应，发挥“进口干中学”作用，推动了工业企业的科技创新。此外，模型11和模型14的工业技术进步回归系数为−0.079 18和−0.068 51，节能减排作用显著。正如理论分析可知，资本品和中间品技术含量较高，可直接发挥技术溢出效应来降低工业空气污染。同时，借助示范与模仿效应、人员培训效应和产业关联效应等动态路径，技术溢出带来的技术进步效应融入工业企业人力资源、创新成果、产业联动等要素，有助于降低空气污染排放。因此，进口贸易技术溢出通过提升工业技术进步降低了大气污染排放，验证了假设3a的合理性。

（四）实证结果分析

表5显示，添加中介变量后，对比模型1和模型4，各拐点值发生了移动与变化。表6比较了基准回归表3和表5中各列回归方程中拐点值变动情况。从模型9和模型12可知，添加工业经济规模变量后，资本品和中间品的技术溢出的两个拐点值分别有所上升和下降，延展了作用于污染排放强度的“N形”关系中的第一段和第三段正效应。在模型10和模型13中，添加工业结构高级化变量后，两个拐点值分别下降和上升，延长了“N形”形状中的第二段距离，增强了技术溢出的负效应。从模型11和模型14可知，控制了工业技术进步变量，资本品技术溢出的两个拐点均下降，而中间品技术溢出的两个拐点值先下降和后上升，这表明资本品技术溢出先后增强了第二段负效应和第三段的正效应，而中间品技术溢出延伸了第二段负效应的距离。

表6 拐点变化趋势

承接上述“N形”拐点值变动情况，再结合前文分析的进口贸易技术溢出环境效应的影响方向可知，在较低的进口贸易技术溢出水平下，未跨越第一个拐点时，工业规模扩张效应占据主导地位，增加了大气污染排放强度。随着技术溢出水平不断渗入工业企业，当跨越了第一个拐点后，位于两个拐点之间时，要素禀赋结构优化效应和工业技术进步效应显现，并超越规模扩张效应，进而降低污染排放。这表明经济发展与技术溢出程度相辅相成，遵循经济增长与环境污染的“规模→技术+结构”路径。最后处于较高水平的进口贸易技术溢出时，即跨越第二个拐点后，资本品技术溢出再次加剧大气污染排放，依据拐点变化情况，主要由技术进步效应作用导致。一方面，技术进步引致的资本化扩充增强了生产规模，再次引发经济规模扩张效应对环境污染的正外部性，继而呈现环境回弹趋势，抵消了技术进步产生的节能环保效应。另一方面，受技术差距、专利垄断、吸收能力等因素制约，不易直接获取资本品内含的技术知识。即使有效汲取了资本品技术溢出，将二次创新成果应用于企业生产实践也具有较长的滞后性。为节约研发成本和快速获取最终工业品，往往技术依赖效应和替代效应在工业企业中应运而生，制约自主研发的内在动力，无法从根源上缓解工业大气污染排放现状。与之不同的是跨越第二个拐点后，致使中间品技术溢出促增排放强度的原因可能只受环境回弹效应的影响，因为中间品作为组装生产的关键材料和零配件，在加工过程中易于获取嵌套在内的核心技术，耗费的时间和成本较少，从而减弱技术锁定效应带来的阻滞自主创新的影响。以上分析验证了假设3b的合理性。

依据上述分析内容，不同程度的进口贸易技术溢出水平通过不同的中介效应来影响工业大气污染排放。因此，依据基准回归模型的拐点值，以样本期间均值和2017年数值为标准，将30个省市划分为不同的类别，探讨各省市处于何种阶段的进口贸易技术溢出，通过何种效应来影响工业大气污染排放，如图2～图5①图2～图5是在国家测绘地理信息局标准地图服务网站（http://bzdt.ch.mnr.gov.cn/index.html）地理底图基础上，添加相关信息制作，审图号：GS（2016）2923号。所示。

图2 按均值分类各区域的资本品技术溢出水平

图3 按均值分类各区域的中间品技术溢出水平

图4 2017年各区域的资本品技术溢出水平

图5 2017年各区域的中间品技术溢出水平

图2和图3结果显示，以均值为研究视角，约11个省市拥有较低的进口贸易技术溢出水平，因工业经济规模扩张效应加剧当地企业空气污染物排放。其余19个省市与之相反，受结构效应和技术效应影响，资本品和中间品技术溢出有效降低了空气污染。图4和图5结果显示，截至2017年，拥有低水平中间品技术溢出的省市包括贵州、甘肃、青海、宁夏和新疆，资本品在此基础上外加山西和云南。仅广东省的资本品技术溢出超越第二拐点，其余省市均处于两个拐点值之间。综合来看，首先，作为实行改革开放最早的东部地区，集群了环渤海经济圈、长三角地区和珠三角地区等高新技术示范区，具有优越的人力资源和研发资金等要素。这些优质的资本为快速吸收和积累国外技术溢出提供了有利条件，使其迅速跨越规模扩张带来的污染正效应。然而，随着贸易规模和进口总量的扩大，广东省的资本品边际减排效应却呈现递减趋势。原因可能在于，自主创新体系和机制不健全、研发投入结构不合理和廉价的劳动力无法短时间吸收资本品核心技术，促使清洁技术研发本土化难以实现，造成技术依赖效应明显。其次，多数中部地区的进口贸易技术溢出也发挥了节能减排作用。在中部崛起战略的推动下，该区域的经济发展和贸易活动取得突破性进展，且煤炭石油等能源资源禀赋条件较好，使得技术溢出在优化工业能源生产结构和提供新能源技术更容易发挥支撑作用。再次，西部地区因闭塞的地理形势和落后的经济发展致使技术溢出的工业规模扩张效应发挥主导作用。虽然西部大开发战略积极推进西部地带融入国际分工和开拓贸易市场，但是粗放型经济增长方式和劣势的技术研发水平无法短时间激发技术溢出的排污效应，而是刺激工业企业投入更多的生产要素和能源资源，加剧当地大气污染。

（五）稳健性检验

为验证中介效应模型的稳健性，本文采用替换空间权重矩阵和改变工具变量方法。将文中的地理空间权重矩阵替换为经济地理空间权重，同时选择最高二阶的空间滞后项作为工具变量。参考Gan[41]具体做法，空间权重的设置如式（8）所示

限于篇幅，具体回归结果在此不予显示。实证结果与上述各变量的影响方向基本保持一致，回归系数有所差异。首先，在基准回归结果中，资本品技术溢出的一次项、二次项和三次项系数分别为5.959 18、−0.444 26和0.010 77，拐点值分别为10.91亿元和804.03亿元。中间品技术溢出的回归系数分别为5.616 22、−0.411 38和0.009 68，拐点值分别为9.56亿元和2 108.19亿元，仍验证了进口贸易技术溢出与工业大气污染排放强度的“正N形”关系。其次，中介效应回归模型显示，进口贸易技术溢出同样在一定程度上扩张了工业经济规模，优化了要素禀赋结构和提升了工业技术进步。添加三个中介变量后，进口贸易技术溢出的系数仍旧显著，证明存在中介效应。最后，对比添加中介变量后的拐点值变化情况（如表7所示），变动趋势与表6显示的结果保持一致，再次验证了实证模型的稳健性。

表7 更换空间权重矩阵和工具变量后的拐点变化趋势

五、结论与建议

本文利用2000—2017年中国30个省市和九个OECD发达国家的面板数据，采用GS2SLS和逐步回归法，以贸易环境效应为空间传导路径，探析资本品和中间品的贸易技术溢出与工业大气污染排放强度的中介效应，得到以下结论：

1.进口贸易技术溢出能够带来工业经济规模扩张效应、要素禀赋结构优化效应和工业技术进步效应三类环境效应，在一定程度上通过扩大工业规模提升工业二氧化硫排放强度，通过优化工业结构和推动技术进步降低工业大气污染。

2.资本品和中间品技术溢出与工业大气污染排放强度呈“N形”特征。在未跨越第一个拐点时，低水平的技术溢出带来的“工业经济规模扩张效应>要素禀赋结构优化效应+工业技术进步效应”，致使工业大气污染显著加剧。在两个拐点之间，技术溢出带来的“工业经济规模扩张效应<要素禀赋结构优化效应+工业技术进步效应”发挥了节能减排作用，减缓了工业大气污染排放。跨越第二个拐点后，高水平的资本品技术溢出带来的环境回弹效应和阻滞自主创新作用，以及中间品技术溢出带来的环境回弹效应再次加重了工业大气污染。

3.从区域角度来看，均值比较下，位于中西部地区的11个省市拥有低水平进口贸易技术溢出，受经济规模扩张效应影响，加剧了大气污染排放，其余省市则处于较高水平进口贸易技术溢出带来的结构效应和技术进步效应阶段。直至2017年，五个省域中间品技术溢出和七个省市的资本品技术溢出仍旧以工业规模扩张效应为主，其余省市则全部进入节能减排层面，只有广东省因环境回弹效应、技术替代和锁定效应再次加剧大气污染。区域的经济水平、吸收能力、技术研发、市场要素流动、地理位置等差异性导致各省市的进口贸易技术溢出水平不同，进而受到异质性环境效应的影响。

根据本文研究结果，提出以下几点建议：

1.推动工业经济规模扩张与大气环境质量的协调发展。实证分析显示，进口贸易技术溢出受规模扩张效应影响，会加剧空气污染排放强度。因此，从理念上，工业企业应当摒弃传统的“先污染、后治理”的粗放型经济发展思想，而是在工业规模扩张初期坚持绿色发展和可持续发展观念，追求工业经济和生态环境的高质量发展，形成大气污染治理与经济规模扩张“协同—耦合—并进”的思想。从执行力上，工业企业应当积极汲取和模仿国外先进绿色科技，将进口贸易技术溢出融入到国内技术应用，重新协调技术、资本和劳动力等生产要素比例，提升自身的管理水平和科技水平，加强精细分工，最终以较少的生产要素换取较高的利润，形成规模经济效应，降低空气污染排放。政府和相关机构应当加强对资源和能源利用的监管，采取税收和法律等政策和制度约束工业企业浪费资源和违法排污等行为。此外，政府、工业企业和第三方数据平台应当联合打造“绿色工业互联网”模式。在规模扩张时期，严格监测工业企业的大气污染排放和处理等实时动态情况，并优化生产工艺，调整高污染行业结构，以环保大数据的手段助力两者协调发展。

2.加强进口贸易技术溢出与自主创新有机结合。自主创新是解决中国工业大气污染排放的根源措施，对进口贸易技术溢出既要积极汲取和引进以获取较高的环境效益，也要实现其与自主创新的和谐发展。对自主研发创新的工业企业，政府应给予最大力度上的政策扶持，通过深化金融市场改革，完善环境规制和知识产权保护政策等加快“贸易技术溢出—消化吸收—自主研发—再创新”的良性生态链。同时，工业企业应当抓住全球数字经济发展机遇，通过进口贸易渠道加强国际技术合作与交流，并采取产、学、研等合作方式积极进行二次研发，加大力度建设智慧环境创新工程、清洁替代创新工程等，加强外资技术和内资技术的创新互补，最终形成适合本土环境需要的绿色技术。

3.因地制宜地扩大进口贸易规模。东部地区如江苏、浙江，中部地区如湖南、湖北等地区现阶段拥有中等水平的国外技术溢出，并未跨越第二个拐点，结构优化和技术进步效应显著。因此，今后针对这些区域，应适度扩大进口贸易规模，防止高度的技术溢出引发工业实体经济的环境回弹效应和阻滞自主创新现象；而西部地区如新疆和甘肃仍停留在较低的技术溢出程度，受规模扩张效应致使工业大气污染排放加剧。因此，针对这些区域，应当高度扩大进口规模。一方面，需要政府的支持和引导。政府应当结构性下调关税，降低绿色产品的进口税率，放宽对环保治污装备、新能源材料和节能与环保服务等市场准入限制。同时，应适当简化环保产品技术的进口流程，并完善相应的公共服务和基础设施建设。另一方面，当地的工业企业要积极抓住扩大进口贸易机会和响应政府的扶持政策。坚持对外开放政策、“一带一路”政策、中部崛起和西部大开发战略，并积极参加国际进口贸易博览会和相关展会，加强国际合作和交流。基于此，东部和中部地区应当注重进口贸易的方向、质量和效益，持续优化进口产品的结构；而西部地区在此基础上应当加大进口贸易的规模，尤其是含有高技术的资本品和中间品的进口规模；同时，要加强建设区域间协同化和一体化发展的联动机制，西部地区应积极学习和引入东、中部地区先进的自由贸易试验区和进口贸易创新示范区的经验，加强区域间的知识溢出和转移，最终形成内部辐射和外部联动的进口贸易扩展格局。