© 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved.
第1期
肖 坤等 : 矿产资源需求的计量分析
27
2 矿产资源需求的模型分析
211 相关研究评价及说明
该领域学者大多运用规范的分析法 , 在纯理论上 研究了矿产资源的开发利用水平 , 动态地说明了为实 现帕累托最优或是潜在帕累托最优 , 矿产资源在若干 时期内的最优开发利用水平 , 为矿产资源的供需确定 理论上的标准和指导 。然而他们不能在实际中说明和 预测矿业活动的规模等状况 , 在制定相关矿产资源规 划时 , 往往也只是根据本地矿产资源的实际条件和往 年供需情况 , 结合相关经济政策和总体规划 , 并考虑 国内国际一般经验 , 初略地估算今后一段时间内矿产 资源的开发和需求用量 , 给出或正或负的供需缺口 。 他们的研究重点一般是当前和过去的情况 , 对未来的 预测并没有较好的令人信服的说明 。
注 : 数据来源于《湖北统计年鉴》、《武汉四十年》、《武汉统计年鉴》、《武汉五十年》。Iron 表示生铁消耗量 , Oil 表示原油消耗量 , Coal 表示 煤炭消耗量 , 下同 。
收稿日期 : 2003 - 11 - 19 ; 作者简介 : 肖坤 ,男 ,中国地质大学 (武汉) 硕士研究生 ,研究方向为产业组织 、环境管理 。 ① 本文受到 “武汉市矿产资源总体规划”项目 (99046031) 资助 。
表 3 矿产资源消耗量差分形式的单位根检验
变量 检验形式 ( С, T , K) ADF 检验统计量
临界值
△2Coal
( c , n ,1)
- 31210
- 21756 3 3 3
△Iron
( c , n ,2)
- 41378
- 41326 3
△Oil
( c , n ,1)
- 81140
- 41221 3
注 : △2Coal 为煤炭 ( Coal) 的二阶差分过程 ; 通常 , “ 3 ”表示 在 1 %的显著水平下 , “ 3 3 ”表示在 5 %的显著水平下 , “ 3 3 3 ”表 示在 10 %的显著水平下 。 21313 VAR 方程的构建
上述三者的 ADF 检验值均小于临界值 , 拒绝单 位根假设 , 因此都为平稳过程 , 此时可以分别求出各 自的向量自回归方程 ( VAR) 。同样观察它们的自相 关系数和偏自相关系数 , 初步确定 △2Coal 的自回归 部分是 5 阶的 , △Iron 的自回归部分是 5 阶的 , △Oil 的自回归部分是 2 阶的 , 实际操作中适当增加或减少 AR 项也是必要的 。构建 VAR 方程 , 从统计检验的 角度考虑 , 略去一些不显著的项 , 有如下形式 (见参 考文献 1) :
1 矿产资源储量及消耗的基本情况
武汉市目前已发现矿种 33 种 , 矿点 154 处 , 潜 在经济价值 8 400 多亿元 。其中冶金辅料和建筑材料 储量最大 , 拥有全国最大的溶剂石灰岩 、白云岩 、石 英砂岩基地 , 主要分布于江夏乌龙泉和黄枝山一带 , 储量 4 000 多万 t ; 黄陂区黄子山 、丰山等地有大量 白云岩分布 ; 石英石矿见于江夏区八分山 , 蔡甸区侏 儒山 和 横 山 等 地 , 含 二 氧 化 硅 97 % 以 上 , 储 量 近 3 000 万 t ;膨润土矿储量 1118 多亿 t , 居全国第一位 ; 优质矿泉水 、热泉已发现 8 例 ; 此外发现金矿点 11 处 , 铜矿点 9 处 , 锰矿 、磷矿各 2 处 , 并有 5 处发现 石油 、天然气的 “油气显示”。由于早期实际统计工 作的不足 , 一些数据是不可考的 , 特别是用量少的小 宗矿产 。因此不得不将其忽略处理 , 而主要集中于用 量大的战略性资源 : 原油 、煤炭 、生铁的消耗统计量 (见表 1) 。
解释变量 (滞后 2 阶)
解释变量 (滞后 3 阶)
解释变量 (滞后 4 阶)
- 01432 ( - 11633) d1 w = 21199
F = 617
解释变量 (滞后 5 阶)
- 01767 ( - 21387)
C
- 181222 ( - 11611)
△Iron
R 2 = 01963
01337 (21002) d1 w = 21534
从武汉市矿产资源储量 、分布 、开发和利用的具 体情况来看 , 关系到后续发展的矿产资源供需形式并 不乐观 , 供需缺口将继续增大 , 政府及经济管理部门 必须对此作出积极的应对 , 因此矿产资源供需形势预 测是一项重要的经济预测工作 。本文的意图即在于运 用科学的经济预测方法对武汉市矿产资源需求形势作 出合理的判断和预测 , 并提出相关的政策建议 。
表 1 武汉市主要矿产消耗量
(单位 : 万 t)
年份 主要
矿产 消耗量
Iron
1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 465 46612 474 46319 50513 53417 52911 54419 55316 60118 618144 623104 647181 682105
注 : (1) 本表的检验及下述各表的统计或检验由 Eviews310 完成 。 (2) 检验形式 ( C , T , K) 分别表示单位根检验方程 , 包括常数项 , 时间趋势和滞后阶数 , N 是指不包括 C 或 T , 加入滞后项是为了使残 差项为白噪声 , 滞后项阶数由 AIC 准则确定 , Δ表示差分算子 。 21312 时间序列的平稳化
表 2 矿产资源消耗量的单位根检验
变量
检验形式 ( C , T , K)
ADF 检验统计量 临界值 (5 %)
Coal
( c , n , 1)
- 11176
- 31148
Oil
( c , n , 1)
- 11056
Hale Waihona Puke - 31148212 本文的研究方法
Iron
( c , t , 1)
- 21789
第 6 卷第 1 期 资源 ·产业
2004 年 2 月 RESOU RCES & INDUSTRIES
资源评价
Vol16 No11 Feb12004
矿产资源需求的计量分析
———武汉市主要矿产资源需求的实例分析 ①
肖 坤 , 杨树旺
© 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved.
28
资 源 ·产 业
2004 年
变量 △2Coal
解释变量 (滞后 1 阶)
- 01497 ( - 21674) R 2 = 01870
表 4 矿产资源消耗量差分形式的向量自回归方程
01449 (31523) F = 191513
01691 (41969)
11065 (71598)
- 191912 ( - 21383)
△Oil
- 11496 ( - 61055) R 2 = 01930
- 11591 ( - 51664) d1 w = 11426
F = 171804
21289 (21650)
的时序模型是对时间序列的严谨的分析方法 。 21311 时间序列的平稳性检验
对于本文的时间序列数据 , 首先将考虑上述数据 的平稳性 , 在这里使用的是单位根检验 (DF) 。由于 在 DF 检验中不能保证方程中的残差为白噪声过程 , 通常对 DF 检验进行扩展 , 形成 Augmented Dickey Fuller 检验 (ADF) 。以下是这三种矿产品的 ADF 检 验结果 :
161477 21828
注 : R 2 为方程拟合优度的判定系数 , d1 w 为检验序列相关性的统计量 , F 为检验总显著性的统计量 , C 为方程常数项 。
213 现代经济计量模型的构建
本文所选用的时序模型区别于传统的基于某种经 济理论的结构模型 。随着现代计量经济理论和建模方 法的发展 , 传统的方法越来越显现出不足 。首先 , 传 统的方法没有考虑到所选数据 , 特别是时间序列的稳 定性问题 , 对于非平稳过程生成的时间序列数据 , 传 统的计量经济方法学是无能为力的 , 因为作为推断和 检验理论基础的中心极限定理不再适用 , 并且有时还 会导致伪回归 。如太阳黑子和失业人口 , 即使他们并 不相关却往往能得到显著的相关系数 , 因为这两者的 时间序列数据有明确的趋势 , 而这种趋势在回归估计 中又起着主要的作用 。其次 , 传统的建模方法实际上 是一种依赖于先验的经济理论的理论驱动方法 , 它只 能起到检验理论的作用而难以发现新的理论 , 现代的 时序建模理论则充分描述数据本身的特征 , 是一种数 据驱动的方法 。再次 , 传统的建模方法也即结构建模 法 , 在检验被拒绝后要根据所反映的问题进行重新设 定和修改 , 直到通过检验 , 这也被称为隐式建模设 计 , 而这种策略显然是孤立和盲目的 。Pyrrho 的理论 表明 , 当建模者想使某些解释变量的系数达到某预期 值或拟合优度达到某一水平时 , 他总能通过向模型中 引入适当的变量来实现这种愿望 。这也说明了隐式建 模的随意性 , 容易将研究者引入歧途 。因此本文采用
(中国地质大学人文与经济学院 , 湖北 武汉 430074)
摘 要 : 本文介绍了武汉市矿产资源存储消耗的现状 ,结合相关研究成果 ,主要应用计量经济学的方法 ,指出 了传统建模方法的不足并运用现代建模方法 ,对矿产资源 (主要指煤炭 、原油和生铁三种战略性矿产资源) 消耗 的时间序列建立了 VAR 模型 。最后运用该模型对武汉市矿产资源需求作了预测并提出有关政策建议 。 关键词 : 矿产资源 ; 时序模型 ; 预测 中图分类号 : F41611 ; O21211 文献标识码 : A 文章编号 : 1006 - 9399 (2004) 01 - 0026 - 04