第24卷第6期 2015年6月 中 国 矿 业 CHINA MINING MAGAZlNE Vo1.24,No.6 煤炭价格波动特征及市场风险研究 宋建新 (中国矿业大学(北京)管理学院,北京100083) 摘要:煤炭价格的大幅波动给煤炭企业及相关产业带来了巨大的经营风险,本文应用GARCH类 模型对煤炭价格的波动特征进行研究,研究发现:煤炭价格波动存在明显的聚集性和长期记忆性;煤炭市 场存在“高风险、高收益”的特征和非对称效应——利好消息对煤炭价格波动的影响要大于不利消息的影 响。其后,通过构建VaR—GARCH族模型对煤炭市场的风险价值进行测算,并与实际损失相比较,发现广 义误差分布下的VaR-EGARCH(1,1)模型能够更好地描述煤炭市场的风险。 关键词:煤炭价格;GARCH类模型;风险价值VaR 中图分类号:C93 文献标识码:A 文章编号:1004—4051(2015)06—0048—04 Price volatility characteristics and market riskof coal SONG Jian-xin (School of Management,China University of Mining and Technology(Beijing),Beijing 100083,China) Abstract:Thevolatility of coal price has brought great risk to coal enterprises and its downstream firms.Therefore,based on GARCH models,fluctuant features of coal is researched empirically,and we find that the price of COS】has long memory and clustering in volatility,the coal market SHOWS the characteristjcs of high——risk and high—。return and the non-symmetrical effect that favorable news have a larger impact on the market than bad news.Then,the VaR of coal market is calculated by VaR—GARCH models,and by comparing with actua1 1OSS we find that VaR—EGARCH(1,1)model based on GED-distribution can describe the risk of coal market more effectively. Key words:price of coal;GARCH models;value at risk 我国具有“富煤、贫油、少气”的能源资源赋存 特点,这决定了煤炭在我国的基础能源地位,其在 一2011年l1月3日,秦皇岛动力煤(山西优混,5500 大卡)的价格由570元/t上涨到860元/t,增长了 50.9 ;而后,煤炭价格开始步入下跌之途,到2014 次能源消费结构中的比重长期高达70 左右。 近几年,为了应对气候变化、节能减排和大气污染 防治,国家提出了优化能源结构、控制煤炭消费总 年2月8日跌至505元/t,降低了41.3 。煤炭价格 的大幅波动对我国煤炭生产及下游企业都产生了严 重的影响,加大了企业经营风险,不利于我国经济社 会的平稳运行和发展。因此,研究煤炭价格波动特 征,度量波动产生的风险,具有重要的现实意义。 现有对煤炭价格的研究大多集中在价格形成 机制『1]、煤炭价格预测[2]、煤炭价格的影响因素[3]、 量等措施,并大力推广清洁能源的开发和利用,但 由于受到技术、资金等因素的影响,清洁能源产量 短期内仍难以大幅增加。因此,在短期内以煤炭为 主体的能源消费结构不会发生根本性的改变,煤炭 在未来较长的一段时间内仍将是我国重要的能源 资源。另外,受宏观经济调整、供需失衡、低价进口 煤冲击等因素的影响,煤炭价格经历了由大涨到大 跌的巨幅波动。如图1所示,从2009年1月4日到 煤炭价格波动对下游产业和宏观经济的冲击¨4 以 及煤电博弈 ]、煤电联动[7 等方面,很少有文献对煤 炭价格的波动特征和市场风险进行量化研究。因 此,本文应用GARCH类模型对煤炭价格的波动特 征进行实证研究,并运用VaR—GARCH模型对煤炭 收稿日期:2015一Ol一05 作者简介:宋建新(1972一),男,山东济宁人,博士研究生,主要从 事资源经济领域研究。E—mail:s ̄xxj@126.com。 市场的风险价值进行了测算,这可为相关企业和政 府部门进行风险控制提供参考依据。 第6期 宋建新,等:煤炭价格波动特征及市场风险研究 样本周期 图1 2009年1月4日至2014年2月l9日 秦皇岛动力煤价格变化情况 1数据与研究方法 1.1数据 本文所用数据为2009年1月4日至2014年2 月19日的秦皇岛动力煤平仓价(山西优混,5500大 卡),共计1277组,数据来源为Wind资讯。为了方 便计算,我们将煤炭价格收益率定义为对数收益率 (式(1))。 —in(pcoal )一in(pcoal 1) (1) 式中 表示第t期的煤炭价格波动率;pcoal 表示 第£期的煤炭价格。 表1给出了序列 的统计性描述结果,煤炭价 格收益率的均值非常接近于0;偏度大于0,峰度远 大于3,J—B统计量也远大于临界值5.992,因此我 们认为煤炭价格收益率序列rf不服从正态分布,并 且具有尖峰、厚尾的特征。 图2给出了煤炭价格收益率序列r 的走势,可 以发现煤炭价格收益率具有明显的波动聚集特征。 1.2研究方法 为研究煤炭价格的波动特征和市场风险,本文 的模型分为两部分。首先,应用GARCH类模型对 煤炭价格的波动特征进行实证研究,该模型能够检 验煤炭价格波动是否具有长期记忆性、是否存在 “高风险高收益”的特征以及非对称效应等。其次, 将GARcH类模型和VaR方法相结合,构建VaR— GARCH类模型来测算煤炭市场的风险价值。 1.2.1 GARCH类模型 为了刻画资产收益尖峰、厚尾及波动聚集的特 征,Engle提出了自回归条件异方差(ARCH)模型。 但当滞后阶数较大时,无约束的估计常常会违 背参数是非负的限定条件,因此Bollerslec又提出 了广义自回归条件异方差(GARCH)模型。假设有 资产收益率Y ,那么GARCH( ,q)模型可表示为 式(2)、式(3)。式(2)、式(3)分别为均值方程和方 表1 煤炭价格收益率序列rt的统计性描述结果 图2秦皇岛动力煤价格收益率走势图 差方程。 Yt—C+∑Yt一十厶 十 -i十 (2), 2—60+∑ai;l 2 +∑ J;1 (3) 式中:地为残差项; 为 的条件方差;C和 为常数 项。而在资本市场上经常存在“高风险、高收益”的现 象,针对这种现象Engle等又提出了GARCH—M模 型,其方差方程仍为式(3),均值方程变为式(4)。 Y 一∑Y + + (4) 式中,参数.0是指可观测到的预期风险波动对 的 影响程度,它代表了风险和收益之间的一种权衡。 另外,为了描述好消息和坏消息对资产价格冲击的 非对称性,Zakoian等提出了TARCH模型,Nelson等 提出了EGARCH模型。两个模型的均值方程仍都 是式(2),而方差方程分别变为式(5)和(6)。 一∞+∑ 2 +∑ dr2.一J+∑OkfL d i=1 J 1 h=1 (5) in(a ̄)= +∑fliin(aL,) J=1 + J 一E( )I+ ㈤ 式(5)中,d 是一个虚拟变量,当 小于0时,d“ 等于1;否则就等于0,因此只要 不等于0,就存在 5O 中国矿业 第24卷 非对称效应。式(6)中, 是非对称效应参数,当 一0时,好消息和坏消息的冲击是对称的;当 >0 时,坏消息的冲击要小于好消息;当y <0时,坏消 息的冲击就会大于好消息。 1.2.2风险价值VaR及其测算方法 所谓风险价值VaR就是指在一定的置信水平 下,某一资产在未来一段时间内可能遭受的最大损 失,其数学定义为式(7)。 VaR=~inf{y lPr(Ⅱ≤ )>1~a)(7) 式中:Ⅱ表示该资产的未来损益;inf{y l A}表示使 A成立的所有Y组成的集合下界; 是给定的置信 水平。因此,VaR是对应于置信水平a损益分布的 下分位数,它描述了资产的未来可能遭受的最大损 失。以未来资产价值的期望值为参考,VaR的基本 计算公式见式(8)。 VaRr—VO (8) 式中: 为资产期初的价值; 为方差; 为下分位 数;T为持有期。传统的VaR计算方法通常假设资 产收益服从元条件的正态分布,但大量的研究表明 资产收益的分布通常具有明显的尖峰、厚尾特征, 收益的波动也具有聚集效应和非对称效应。因此, 我们用GARCH族模型来描述价格波动的时变性 和聚集效应,构建VaR—GARCH族模型,如式(9) 所示o VaR 一 P 1,t=1,2,…, (9) 式中,P 是t一1时刻资产的价格, ̄/ 是由 GARCH族模型计算得到的资产收益波动的条件标 准差,z。为资产收益服从某种分布时的a分位数。 2 实证分析 2.1数据的平稳性和ARCH效应检验 在构建GARCH类模型之前,要先确定所用数 据是否平稳并且存在ARCH效应。首先,我们采用 ADF单位根检验法来验证煤炭价格收益率序列r 的平稳性。由式(10)可知最大滞后阶数为23,由图 2可以知序列r 不包含截距项和趋势项。ADF检 验结果如表2所示,因此我们拒绝序列r 存在单位 根的原假设,并认为它是平稳的。 一『12*( )i] ㈣ 式中:lOmax为Schwert(1989)建议的最大滞后阶数;N 为样本容量;口为取整符号。 经过多次的试验,参考AIC和sC准则,我们最 终确定煤炭价格收益率序列GARCH模型的均值 方程,见式(11)。 r — 1r 1+ 2 1+ (11) 对均值方程(11)的残差项进行ARCH效应的 拉格朗日乘数LM检验,如表3所示,均值方程(11) 存在高阶的ARCH效应,说明所选数据适合构建 GARCH类模型。 表2煤炭价格波动率序列Rt的ADF单位根检验 注:***表示在1 显著性水平下显著。 表3均值方程残差项的ARCH—LM检验结果 F统计量8.8750…7.4007…6.6100… 户值0.00o0 0.o000 0.0000 注:***表示在1 显著性水平下显著。 2.2 GARCH类模型的构建 为了考察煤炭价格波动的不同特征,我们分别构 建序列 的GARCH(1,1)、GARCH(1,1)一M、 TARCH(1,1)以及EGARCH(1,1)模型,结果见表4。 1)GARCH(1,1)模型中ARCH项、GARCH 项的系数a和 分别为0.1267和0.8269,且均在 1%的显著性水平下通过检验;二者之和为0.9536, 小于1且非常接近于1,表明外部冲击和自身过去 的波动对煤炭价格的影响具有较强的持续性,煤炭 价格波动具有长期记忆性。GARCH项的系数明显 大于ARCH项系数,说明煤炭价格的波动主要是由 自身过去的波动引起的。 2)GARCH(1,1)一M模型中风险溢价系数ID为 0.1540,且非常显著,说明当煤炭市场中的预期风 险增加1 时,收益率就会相应地增加0.1540 ,即 煤炭市场存在“高风险、高收益”的特征。 3)TARCH(1,1)模型和EGARCH(1,1)模型 中代表非对称效应的参数 和7分别为一0.0992 和0.2902,均小于0且都在1 的显著性水平下通 过检验,说明煤炭市场存在非对称效应,即“好消 息”对煤炭价格波动率的影响要大于“坏消息”。 4)对上述四个模型的残差项分别进行ARCH— LM检验,发现F统计量均不显著,即这四个模型均 不再存在ARCH效应,表明拟合效果较好。 2.3 VaR值的测算 广义误差分布GED是一种比正态分布更加灵 活的分布,通过对参数的调整可以拟合不同的分布 第6期 宋建新,等:煤炭价格波动特征及市场风险研究 厚度参数。当 <2时,GED为厚尾分布;当 >2 时,GED为瘦尾分布;当 一2时,GED退化为正态 厂 ' 一 一 町而L L I A I-I ‘ 2 分布。因此,为了更好地拟合煤炭价格收益率波动 的厚尾特征,下面分别应用正态分布假设下的 GARCH(1,1)模型和广义误差GED分布假设下的 … 斗罐 … 系数 GARCH(1。1) Ai EGARCH(1,1)模型来计算煤炭价格收益率序列 的条件方差h ,再将 代入VaR的计算公式(9), 就可得到煤炭市场的风险价值VaR,结果如见5。 表4序列n的GARCH类模型估计结果 模型 GARCH(1,1)一M TARCH(1,1) 0.9448…一 0.9800… 均值方程 2 0.8428…0.9788… 4.05*10 7… 一14.0539…0.4452…0.4090… 0.157l… 方差方程 口 0.8379…0.0992…y 0.29O2…0.6703 0.5569 — ARCH LM 相关统计量 AIC SC —7.6912 7.6670 7.2287 7.2044 —注:***表示在1%显著性水平下显著;**表示在5 水平下显著;*表示在1O%水平下显著。 表5煤炭市场的风险价值VaR测算结果 由表5可知,当置信水平为95%和99 时,正 态分布假设下的VaR—GARCH(1,1)模型的测算值 小于实际损失的比例显著地大于5 9/6和1 ,说明该 模型低估了风险;而广义误差分布下的VaR— EGARCH(1,1)模型的测算值小于实际损失的比例 非常接近5 和1%,说明该模型对煤炭市场风险的 测算比较准确。 3结论及建议 游企业在研究煤炭市场时,不仅要重视供需基本面 分析,还要加强对煤炭历史价格波动特征的研究。 2)煤炭市场存在“高风险、高收益”的特征,当 煤炭市场中的预期风险增加1 时,收益率就会相 应地增加0.1540 。这说明煤炭市场上的交易者 大多根据理性因素进行决策,煤炭市场交易机制建 设比较完善。 3)煤炭市场存在非对称效应,即“好消息”对煤 炭价格波动的影响要大于“坏消息”的影响。因此 应用GARCH类模型对煤炭价格的波动特征 进行实证研究,然后构建VaR—GARCH模型对煤炭 市场风险进行测算,主要得到以下结论。 1)煤炭价格收益率的波动存在明显的长期记 忆性,影响的时间比较长,因此建议煤炭生产及下 建议相关部门加强市场信息监管,防止有不法 分子通过发布不实信息来操纵煤炭价格。 4)煤炭价格收益率具有明显的尖峰、厚尾和波 (下转第56页) 56 1966:6—7. [[ [ [ [ [ 4 ]5 ] 6 ] 7 ] 8 ] 9 ]中国矿业 15一Z3. 第24卷 chard R.Nelson.A Theory of the Low LevelEquilibrium [3] Ri[16]黄子健.“资源诅咒”:测度与破解——一个来自的经验 证据[J].财经大学学报,2013(4):19—27. Elf]邵帅,齐中英.西部地区的能源开发与经济增长一一基于“资 源诅咒”假说的实证分析[J].经济研究,2008(4):147—160. [18]李影.“福音”还是“诅咒”——自然资源与经济增长关系研究 综[J].经济论坛,2009(6):19-22 Trap in Underdeveloped Economies[J].American Economic Review,1956,46:894—908. Debr ̄Ray.DevelopmentEconomics[D_.Princeton University,1988. 吴国宝.对中国扶贫战略的简评[J].中国农村经济.中国农村 经济,1996(8):26—30. 李实,古斯塔夫森.八十年代末中国贫困规模和程度的估计 [J].中国社会科学,1996(6):13—34. [19]崔学锋.“资源诅咒”论不成立[J].经济问题探索,2013(5): 27—31. 郭宏宝,仇伟杰.财政投资对农村脱贫效应的边际递减趋势 及对策口].当代经济科学,2005(5):53—57. [2O]世界银行东亚及太平洋地区扶贫与经济管理局.从贫困地区 到贫困人群:中国扶贫议程的演进中国贫困和不平等问题评 估JR].2009. [21]成升魁,苏筠.西部地区资源开发与消除贫困[c]//西部大开 发科教先行与可持续发展——中国科协2000年学术年会文 集,2000. 刘晓昀,辛贤,毛学峰.贫困地区农村基础设施投资对农户收 入和支出的影响[J].中国农村观察,2003(1):31-36. 于海龙,宋山梅.土地特征与粮食生产、农村贫困关系实证分 析一以贵州喀斯特贫困地区为例[Jl江西农业大学学报:社 会科学版,2009(3):47—53 xed [io] Glyfason,T.,Herhertsson,T.T,Zoega,G.A Mi[22]郑雯雯.资源依赖型区域贫困化增长:形成机理与破解途径 [D].兰州:兰州大学,2008. Blessing:Natural Resources and Economic Growth[J]. Macroeconomic Dynamics,1999(3):204—225. [23]魏燕妮.陕北资源富集型贫困的制度成因与对策研究[D].西 安:西安科技大学,2011. s,G.A.,J.E.Tilton.The Resource Curse[J].Natural [11] DaviResources Forum,2005(29):233—242. s,E,Gerlagh R.Resource Abundance and Economic [12] Papyraki[24] 王雅敬,任晓冬,苏维词,等.贵州喀斯特地区资源性企业参 与扶贫模式研究——以贵州J矿社区发展为例[J].贵州 师范大学学报:自然科学版,2011,29(2):l14一l19. Growth in Unite States[J].European Economic Review,2006 (4):253-282. [25]宋文飞,李国平,韩先锋.矿产资源地贫困理论机理及改 革——基于税制与“双重扭曲”的视角口].制度经济学研 究,2013(1):123—146. [13] 徐康宁,王剑.自然资源丰裕程度与经济发展水平关系的研 究[J].经济研究,2006(1):78—89. [26]王永平,黄海燕.新阶段贵州农村贫困特征与反贫困策略调 整[J].贵州农业科学,2010,38(7):204—208. [14] 徐康宁,邵军.自然禀赋与经济增长:对“资源诅咒”命题的再 检验[J].世界经济,2006(II):38—47. [15] 胡援成,肖德勇.经济发展门槛与自然资源诅咒——基于我 国省际层面的面板数据实证研究[J].管理世界,2007(4): [27] 贵州省国土资源厅.2012贵州省国土资源公报[R].2013. [28]袁国华,资源基础与中国工业化战略[M].武汉:中国地质大 学出版社,1999:35—92. 米米米米米米米米米米米米米米米米米米米米米米米米米米米柴米米柴桨米米米米米米米米米米米米果米米米米 (上接第51页) 动聚集特征,传统的VaR计算方法会低估实际的风 险值。另外,广义误差分布GED假设下的VaR— EGARCH(1,1)模型对煤炭市场风险的测算更加准 的应用口].煤炭学报,2001,26(1):109一u2. [3] 丁志华,赵洁,周梅华.基于VEC模型的煤炭价格影响因素研 究[J].经济问题,2011(3):45—48. [4] 牟敦果,林伯强.中国经济增长、电力消费和煤炭价格相互影 响的事变参数研究[J].金融研究,2012(6):42—53. [5] 王伟,董德彪,赵联振,等.基于投入产出模型的煤价变动对 确。因此,建议煤炭生产及下游企业采用该模型进 行风险价值测算,以提高风险管理能力。 ■■● 国民经济影响分析[J].煤炭科学技术,2012,40(12): 125—128. 参考文献 [1] 张华明,赵国浩.煤炭价格形成机制存在的问题及对策分析 [6] 李丽,杨力,韩静.煤电一体化趋势下煤电联营的博弈分析 口].中国煤炭,201I(I):30—32. r^1 [J].资源科学,2010,32(11):2210—2215. 刘希颖,林伯强.改革能源定价机制以保障可持续发展—— [2] 王立杰,刘志东.经济时间序列分析技术在煤炭价格预测中 以煤电联动为例[J].金融研究,2013(4):112—126.