1、第 卷第期 年月中国矿业 ，收稿日期：责任编辑：刘硕基金项目：国家自然科学基金面上项目“共伴生关键矿产多尺度峰值预测及预警研究”资助（编号：）第一作者简介：周之逸（），男，江苏南通人，硕士研究生，主要从事应用统计与矿产资源经济研究，：。通讯作者简介：徐德义（），男，博士，教授，博士生导师，主要从事应用统计与矿产资源经济研究，：。引用格式：周之逸，徐德义，樊锦绮，等 高质量发展视角下我国矿业经济发展影响因素分析 中国矿业，（）：，（）：文章编号：（）：高质量发展视角下我国矿业经济发展影响因素分析周之逸，徐德义，樊锦绮，龚玉娜（中国地质大学（武汉）经济管理学院，湖北 武汉 ）摘要：高质量发展理念背
2、景下，我国矿业经济的发展是兼顾各方面多要素的均衡发展，各要素既是衡量矿业经济高质量发展的评价指标，也是影响矿业经济发展的重要因素。如何在符合高质量发展理念的框架下抓住矿业经济发展动力，发挥各地区位优势，是促进我国矿业经济发展、保障国家初级产品供给的重要一环。本文从环境、技术、产业、管理个方面选取影响矿业经济效益的因素，依托 年 个省级行政单位的面板数据，采用基于随机森林的多重插补法对数据进行补全，从而建立面板固定效应模型和面板向量自回归模型，从静态和动态的角度研究我国矿业经济的影响因素。研究结果表明：产业因素是推动我国矿业经济发展的主要因素，环境因素对矿业经济存在一定的推动作用，管理因素对矿业
3、经济存在一定抑制作用。在全国层面，应当规范市场环境，引进投资，促进矿业企业健康发展，同时深入贯彻“两山”理论，推进绿色矿山建设，促进环境保护与经济发展协同增效。另外，本文根据各地区矿业经济发展主导因素的不同，将 个省级行政单位分为环保反馈型、管理优先型、产业驱动型、技术支撑型种类型，各地区应当结合自身区位优势发展矿业经济。关键词：高质量发展；矿业经济；多重插补；面板变系数回归；面板向量自回归中图分类号：；文献标识码：，（，（），）：，：，第期周之逸，等：高质量发展视角下我国矿业经济发展影响因素分析 ，：，：；引言我国是矿产资源大国，矿业既是我国的基础工业部门，也是我国经济的重要组成部分，为我国
4、快速发展提供巨大动力，也为经济增长、减贫以及加强与全球经济融合提供机会。但无论是初期对矿产资源的勘探开采，还是后期的冶炼加工，矿业带来经济效益的同时，也对生态环境和经济社会造成破坏。在生态文明建设和“三条红线”硬约束的背景下，我国矿产资源开发与环境保护的矛盾异常突出。一方面，由于环境承载力的限制，对矿产资源开发产生约束；另一方面，由于生态保护逐步落实，良好的环境也会形成反馈，促进矿业经济的发展。但除了生态环境，技术、产业、管理等因素同样也会对矿业的发展施力。一部分学者认为，由于已探明储量充沛，技术已经不是制约矿产可供性的主要因素，社会矛盾和环境限制才是影响矿产可供性的关键因素。另一部分学者认为
5、，科技创新是引领经济高质量发展的核心驱动力，提高科技创新能力是实现我国矿业经济高质量发展的关键；地质矿产科技创新对矿业经济增长的作用在短期存在滞后性，而长期促进作用显著。此外，我国自然资源政策经历了早期发展后，更加注重系统性、协调性、权威性、公正性，更加注重以科学发展观为指导，为国家可持续发展做出重要贡献，因而其对于矿业的规范、引导作用也不容忽视。同时，矿山管理也是影响矿业经济的重要因素，包括美国地质调查局（）、标普全球市场财智（）等在矿山可供性评估中都加入了管理、安全、政治等风险评估考量。“十三五”规划和“十四五”规划中明确提出，要提高矿产资源开发保护水平，发展绿色矿业，建设绿色矿山。绿色矿
6、业是指以可持续发展为理念，包含生态环境、安全、科技创新等内容的矿业发展的高级形态，这要求各地区发展矿业时不仅要尊重自然，也要遵循自然资源价值规律、市场经济规律和人类社会发展规律 。从“九五”到“十四五”期间，我国对矿业发展的要求也从“资源节约、环境友好”转向追求“人与自然和谐”的绿色、高质量发展。目前对于矿业经济的研究主要是基于高质量发展框架进行的评价指标体系构建。当前主流的矿业经济高质量发展指标选取，除了最核心的经济效益和环境因素外，逐步加入了技术创新、产业投资、资源政策、社会就业等方面的指标。高苇等 基于矿业产值、矿业损坏土地面积和矿业从业人数等指标构建指数以衡量矿业绿色水平；任思达等 构
7、建我国矿业经济绿色发展评价指标时，将经济、环境、政策等个方面作为主要考虑因素；黄洁等 构建包含资源、环境、转型发展和安全和谐等个方面的矿业绿色发展指数用于反映我国绿色矿业的发展情况；等 基于驱动压力状态影响响应模型，综合经济、环境、技术、管理等方面，建立了由 个指标组成的评价指标体系，用于评价我国绿色矿山建设现状；等 选取 项指标，构建山西省矿产资源开发和生态环境评价指标体系。基于高质量发展理念，新形势下矿业经济的发展是兼顾各方面多要素的均衡发展，各要素既是衡量矿业经济高质量发展的评价指标，也是影响矿业经济发展的重要因素。而各要素对于矿业的经济效益也产生了不同影响。同时由于各地区的区位条件差异
8、，各要素对于地区矿业经济的影响也不尽相同。如何在新形势下，既遵循高质量发展理念，又能有效促进矿业经济发展，提高矿业经济效益，带动地区经济发展，对促进我国矿业发展、实现我国经济高质量发展目标具有显著意义。因此，本文基于已有的矿业经济高质量发展评价指标体系，从环境、技术、产业、管理个方面选取影响矿业经济效益的类因素，依托 年跨度为 年的中国 个省级行政单位的面板数中 国 矿 业第 卷据，研究矿业经济高质量发展的个要素对于其经济效益的影响。研究设计与数据说明 研究设计与模型设定本文参照高质量发展内涵和矿业经济绿色发展研究，将矿业经济的高质量发展要素分为环境、技术、产业、管理个方面，这个方面既构成了矿
9、业经济高质量发展的综合评价指标，也作为矿业经济发展的投入要素，对矿业经济的发展造成正面或负面的影响。本文搜集到类要素加上矿业经济和地区宏观经济总计个维度、个变量，使用基于随机森林的多重插补填补了缺失值，在确保插补效果的前提下完善了研究数据与资料，为相关领域的后续研究做出一定的贡献，并在此基础上进一步研究矿业经济高质量发展的类要素对于经济效益的影响。静态模型本文选择面板回归模型作为静态研究的基准模型见式（）。（）式中：为地区；为年份；为矿业的经济效益；为矿业的高质量发展要素，包含环境因素、技术因素、产业因素和管理因素类因素；为地区因素；为时间因素；为残差项。地区经济的整体发展往往会带动地区大部分
10、甚至所有产业的共同发展，迁移到本文的研究中，某地区矿业经济的发展，可能是受到整体经济水平的发展带动而非矿业经济高质量发展投入要素的影响。为了避免这一干扰对于本文研究的影响，加入衡量当地经济发展水平的变量 ，以达到在本模型中消除整体经济水平发展对于矿业经济影响的目的。在基准面板回归模型的基础上代入控制变量的思想，依次对类要素进行面板回归，进一步研究剥离了整体经济发展趋势对矿业经济的影响。动态模型生产要素的投入对于产业的产出影响存在复杂的动态关系，并且可能存在滞后效应和内生性问题，使用传统的面板回归模型无法体现与衡量动态性和滞后性，因此本文在静态研究的基础上，建立面板向量自回归（）模型，见式（）。
11、，（）式中：为地区；为年份；为滞后阶数；为 的一阶差分；为残差项。指标选取与数据处理 样本选取与介绍本文参照任思达等 构建的包含经济、环境、政策、技术、产业、社会方面内容的矿业经济绿色发展水平综合评价指标体系，并综合高质量发展内涵，从矿业经济、宏观经济、环境、技术、产业和管理个维度分别收集 个指标，使用基于随机森林的多重插补法填补缺失值，完善相关数据，并在此基础上对每类因素分别进行熵值法计算，对数据进行降维，从而得到每个因素的综合指标，综合指标的数值越大，代表该因素在地区的发展水平较高。考虑到各地区由于经济基础、历史因素、地理条件等区位因素造成的总体状态上的差异，使用绝对数值的变量来代表各地的
12、各指标水平不适于横向比较，因此本文尽可能地采用相对指标进行建模，以便于各地区之间的横向比较（表）。表变量说明 一级指标二级指标单位类型数据来源矿业劳动效率（）正向矿业产值矿业从业人数 中国工业统计年鉴矿业经济矿业产值增长率（）正向 中国工业统计年鉴（）矿业收入利润率（）正向矿业利润总额矿业销售收入 中国工业统计年鉴非油气综合利用产值增长率（）正向 中国矿业年鉴人均国内总产值增长率（）正向 中国矿业年鉴地区经济（）固定资产投资效率（）正向地区生产总值地区固定资产投资额 中国统计年鉴产业升级系数（）正向第一产业占比第二产业占比第三产业占比 中国统计年鉴第期周之逸，等：高质量发展视角下我国矿业经济发
13、展影响因素分析续表一级指标二级指标单位类型数据来源矿业废水排放总量增长率（）负向矿业污染物排放量全国矿业行业污染排放量矿业废气排放量增长率（）负向地区矿业产值全国矿业产值环境矿业二氧化硫排放总量增长率（）负向 中国工业统计年鉴（）矿业粉尘排放量增长率（）负向 中国环境统计年鉴矿业固体废物排放量增长率（）负向本年矿业开采新增占用、损坏土地面积（）负向 数据库本年矿山恢复治理面积（）正向 数据库国家级绿色矿山试点数量占比（）正向国家级绿色矿山试点单位数地区矿山企业数国家级绿色矿山试点单位公布名单统计（年开始统计，之前均记为）技术矿业单位产值能耗（）吨标准煤万元负向 中国能源统计年鉴（）从业人员（）
14、人正向 数据库经费支出（）万元正向矿业就业人数占地区总就业人数比重（）正向 中国统计年鉴矿业固定资产投资金额增长率（）正向国家统计局产业非油气年矿产量增长率（）正向 中国矿业年鉴（）非油气大中型矿山占比（）正向 中国矿业年鉴 国土资源年鉴矿业产值占工业行业比重（）正向 中国工业统计年鉴管理矿产资源勘查开采违法案件查处数量（）件负向（）有效采矿许可证（）个正向 中国矿业年鉴省、市、县级国土资源管理机构数（）个正向 数据来源与处理由于高质量发展最早于 年在党的十九大提出，且受限于矿业产业划分、数据可获取性等因素，本文选用 年我国 个省级行政单位（香港、澳门、台湾除外）的相关年度数据。本文所选数据主
15、要来自 中国统计年鉴 中国矿业年鉴 中国工业统计年鉴 中国环境统计年鉴 以 及 全球统计数据分析平台和各省级行政区的政府工作报告。缺失变量情况如图所示。由图可知，缺失最多的是本年矿业开采新增占用、损坏土地图缺失数据总览图 （注：空白处表示数据缺失）面积，大约缺失，绝大部分变量的缺失数据在 以下。图为各变量缺失情况的相关系数矩阵。鉴于多维度、短序列、缺失规律性不强的缺失情况与特点，本文采用基于随机森林的多重插补法补全了含有缺失值的所有变量。基于随机森林的多重插补法对连续型变量和分类型变量都有很好的插补效果，且无需先验知识。该方法首先基于随机森林预测缺失值，得到多个初始数据集，然后通过多次迭代来减
16、少各初始数据集之间的差异，最终让所有的预测数据集趋向于一个共同的结果。图缺失变量热力图 （注：系数代表两个变量缺失情况的相关系数）中 国 矿 业第 卷插补过程与结果的平均折线图如图所示。由图可知，个初始数据集在经过次迭代后，插补结果都大致在同一趋势下，数据集的均值走向趋于一条直线，变量都趋于均值收敛，因此插补结果有效。实证研究 矿业经济高质量发展要素静态影响分析 平稳性检验若数据不平稳，模型可能产生伪回归等问题，对面板模型的估计结果和脉冲响应分析结果都产生巨大的偏差和影响。传统的面板单位根检验通常使用 检验，但这种检验方法对于短面板数据并不可用。而 检验和 检验则克服了这一缺点，检验更是在检验短面板数据的基础上，放宽了同质性假定，允许单位根在不同的面板单位中变化，更接近客观事实。因此，本文分别用 检验、检验和 检验等种方法对、和 进行检验，检验结果见表。由表可知，可以拒绝非平稳原假设的序列是、和 ，无法拒绝非平稳假设的序列是 、和 。在经过一阶差分后，所有序列都可以拒绝非平稳原假设。图多重插补过程和结果 （注：子图标题表示变量的序号，黑线表示模拟的插补数据）表平稳性检验 序列名称检验方