2024年1月3日发(作者:出数学试卷答案用什么APP)
第36卷第6期Vol.
36
No.
6统计与信息论坛JOURNAL
OF
STATISTICS
AND
INFORMATION2021年6月Jun.2021【资源与环境统计】结构突变下国际油价波动与中国石油行业
发展的多重联动关系吕政1,胡晨沛2(1中央财经大学
统计与数学学院,北京102206
\'.国家统计局
国际统计信息中心,北京100826)摘要:国际油价突变和中国原油对外依存度攀升增加了国际石油市场与中国石油行业联系的复杂度,廓
清两者关系有助于提升中国在国际石油市场的话语权。为此,厘清了国际油价波动与中国石油行业股指变
动之间的传导机制。在实证检验部分,应用Bai-Perron内生结构突变检验方法诊断突变点,在均值方程
XVAR模型中将结构突变虚拟变量设置为外生变量,以此为基础构建BEEK、CCC、DCC等多元GARCH模
型,对国际油价波动与中国石油开采、石油加工、石油贸易行业股指收益之间的波动溢出性、静态相关性和动
态相关性等多重关系提供新的经验证据)研究发现:国际油价波动对中国石油行业股票价格指数变化起先
导作用;中国石油开采、石油加工行业股指收益率与国际油价变动率存在双向的波动溢出效应;中国石油开
采、石油加工、石油贸易行业股指收益率与国际油价变动率之间的动态相关系数均呈现明显的时变特征。中
国作为发展中的世界第二大经济体,应抓住国际油价低位徘徊的有利窗口期,提升买方市场的定价话语权)关键词:国际油价;石油行业;波动溢出性;时变性;多元GARCH中图分类号:F224. 0
:
F832.
5 文献标志码:A
文章编号#007
—
3116(2021)06
—
0073
—
11引用格式:吕政,胡晨沛.结构突变下国际油价波动与中国石油行业发展的多重联动关系统计与信息论
坛,2021,36(6):on Form:LYU
Zheng,
HU
Chen-pei.
The
Mutual Linkage
between
International
Oil
Price
Fluctuation
and
the
Development
of
Chinese
Oil
Industry
under
the
Structural
Change[J$.
Journal
of
Statistics
and
Information,
2021
36(6):73-83&一、引言2020年中国原油对外依存度攀升至72.7%,伴随国内人民生活水平持续提高以及工业化继续推进,对
化石能源的需求还将不断增加。《世界能源展望2017之中国特别报告》预测,中国将于2030年超越美国成
为世界第一石油消费国,并且在2040年占据全球近30%的国际石油贸易,成为仅次于美国的第二大市场。
显然,随着中国对国际石油市场的依赖不断加深,中国对国际油价波动的反应也将愈发敏感。在短期,国际
油价下跌有利于中国经济边际企稳,同时有利于减少外汇支出,平衡国际收支,但不排除通货紧缩的风险。
2014年,国际油价从每桶110美元暴跌至每桶30美元,此后两年在低位徘徊。2016年以来,国际油价筑底
企稳、震荡上行。2020年,受俄罗斯与OPEC组织减产合作破裂以及新冠肺炎疫情全球爆发综合影响,加剧
石油供给需求不平衡,导致原油价格再度大幅下跌。在国际石油贸易中,油价暴跌常常会将巨额财富从石油
出口国转移到石油进口国,中国作为全球最大的石油净进口国,国际油价持续走低对宏观经济运行的影响更收稿日期:2020
—07
—28;修复日期!020
—12
—15基金项目:全国统计科学研究一般项目“大数据背景下组合风险的统计估计及应用研究\"(2018LY41);贵外|省教育厅科技人
才成长项目“考虑非同步交易影响的金融高频协方差阵的估计及应用”(黔教合KY字〔2018860);中央财经大学
研究生科研创新基金资助项目“制造业集聚与生产率的非线性关系研究一#以运输条件为视角\"(20192Y011)
作者简介:吕政(通讯作者),男,浙江嘉兴人,博士生,研究方向:数量经济学理论与应用;胡晨沛,男,浙江永嘉人,中级统计师,研究方向:世界经济)73
统计与信息论坛深。值得关注的是,与发达国家的石油企业不同,中石化、中石油和中海油等国有企业主导了中国石油行业,
国内石油企业垄断程度高、竞争能力较弱,对国际油价信号的反应灵敏度较差,国际石油市场和中国石油行
业之间存在何种联动关系,已成为广受关注的重要课题)石油资源作为中国石油行业生产经营的重要投入要素,其价格波动将直接影响石油行业的经营绩效)
在实体经济与资本市场联系愈发密切的当下,石油板块股票价格指数是反映石油企业实际运营情况的重要
综合性指标,那么国际油价波动对中国石油行业的冲击必然会体现在石油板块的股票指数上。Broadstock
等的研究表明,国际油价波动对中国股票市场的冲击不容小觑1。基于现有研究证据,本文着重探讨国际石
油价格波动与中国石油行业股指变动之间的关系,结合国际油价曾出现暴涨暴跌的特点,在实证分析中引入
结构突变虚拟变量,以提高实证检验的准确性,并进一步应用多个多元波动率建模工具,分别考察国际油价
波动与中国石油开采、石油加工、石油贸易行业股指收益之间的波动溢出性、静态相关性和动态相关性,借此
分析国际油价波动与中国石油行业发展之间的多重联动关系,以期为中国石油行业发展建言献策。二、文献综述尽管中国股票市场经过近30年的快速发展,已成为全球金融市场的重要组成部分,但与发达国家的金
融市场相比尚处在成长期。石油作为重要战略资源,其价格变动必然会波及中国股票市场,那么油价波动的
冲击幅度有多大?相反,中国股票市场对国际油价是否也存在价格信息外溢?以上问题越来越受到学界和
业界的关注。国外学者较早关注到石油价格波动对股票市场的冲击,但早期研究成果结论差异较大。Hamilton首次
发现国际油价反向影响美国股市,引起了大量学者对这一问题的持续深入研究Ferson等将石油价格视
为风险因子,发现油价对全球多个股市产生了不同程度的影响3。发达国家和发展中国家的股市在面对国
际石油价格波动时,资产价格的波动幅度往往不同。Park等研究发现,石油价格波动对发展中国家股票收
益率的冲击较发达国家更显著4。之所以如此,是因为非石油输出的发展中国家,经济起飞阶段的发展模式
较为单一,对石油的依赖程度远远高于发达国家。Nandha等细化了股票市场的行业,发现除石油、天然气
以及采掘行业外,石油冲击对大部分行业都存在负面影响。国际油价波动对股票市场的影响通常具有杠
杆效应,即油价上涨对石油净进口国股市冲击的影响远比油价下跌大得多。2010年中国成为世界第二大经
济体后,经济继续保持较快增长,与此同时,石油对外依存度以及金融市场对外开放程度也在提高,国际油价
波动对中国经济的影响加深。对此,国内学者就国际油价波动对中国股市的冲击同样做了广泛且深入的研
究,许金华和范英采取GARCH模型和多因素动态模型,观察了金融市场驱动因素和供需基本面在石油定
价机制演化路径中的作用,发现油价波动能够被金融市场驱动因素解释归。杨胜刚等应用VAR-
TGARCH-AGDCC多元模型,从动态相关性的角度检验了大宗商品市场与股票市场之间的联系,发现油市
和股市的联动性正不断增强7。李素芳等构建贝叶斯机制转换协整模型考察了国际石油价格与亚太股票市
场之间的非线性动态关系,发现两者存在显著的非对称效应王朝阳等应用BEKK模型观察了国际油价
与中国新能源股票价格之间的波动溢出效应9。综上,现有的研究成果得到了国际油市和国内股市之间存在双向互动关系的结论,但还存在以下缺陷:
其一,在研究对象上,大多考察国际油价对中国股市整体的影响,但事实上,股市中与石油关系微弱的行业受
到油价波动的影响极其有限,以至于用股票价格综合指数作为油价冲击对象并不精确;其二,在研究方法上,
主要运用VAR模型和多元GARCH模型讨论价格波动之间的互动关系,但大多未考虑当金融时间序列存
在结构突变时,可能造成的实证结果失真。鉴于此,本文以结构突变为视角,将结构突变虚拟变量作为外生变量引入传统的VAR模型,在XVAR
模型的基础上,应用静态的BEKK模型、CCC模型和动态的DCC模型,分别检验国际石油价格波动和中国
石油行业股指变动之间的波动溢出、静态相关和动态相关关系。与以往研究相比,本文的边际贡献包括以下
三个方面:第一,指标选取方面,由于石油板块股指受国际油价波动的冲击最为直接,因此选取中国石油开
采、石油加工和石油贸易行业股指衡量国内石油行业经营状况,观察国际油价波动对中国石油行业在资本市
场上的冲击;第二,研究方法方面,借助Bai-Perron内生结构突变检验方法诊断国际油价的结构突变点,将
74
吕政,胡晨沛:结构突变下国际油价波动与中国石油行业发展的多重联动关系突变点设为外生的虚拟变量,进而构建包含结构突变外生变量的XVAR模型,并在此基础上建立多个多元
波动率模型,增强了实证结果的可信度;第三,理论分析方面,厘清国际油价与国内石油板块间的传导机制,
并从需求信息的角度考虑了中国股市对国际油市的反馈,进一步丰富了能源金融领域的研究成果。三、作用机制分析中国已成为全球最大的石油净进口国,石油能源作为战略资源的重要性不言而喻。韩立岩等的研究表
明,国际油价信号与中国实体经济及资本市场运行具有长期均衡关系[10]。众所周知,国际油价波动将直接
影响国内石油行业的经营业绩,而石油行业的股票指数又是资本市场上反映石油企业实际经营绩效的综合
性指标,那么国际油价波动对石油行业的冲击自然会在石油板块股指上有所反映。石油行业之间具有成本
接近和产品同质的特征,业绩易受共同因素影响,股价变动存在较明显的板块联动效应,进而造成股票市场
整体的次生冲击。根据石油产业链的划分方式,石油产业链上游行业从事石油的勘探开采,中游行业负责石油的提炼加
工、下游行业进行石油销售等业务。当国际油价下跌时,中国石油开采行业在采掘成本上不再具有竞争力,
因此利润减少业绩恶化,然而对中下游的石油加工、石油贸易行业带来了利好,国际油价下跌降低了加工环
节的成本,贸易行业又能够在一定程度上转嫁油价波动的风险。总体而言,国际石油价格下落会降低石油上
游行业业绩,提高中下游行业业绩,并且这些市场价格信号会在股票指数上加以体现。对于与石油产品联系
紧密的消费投资行业而言,如果石油价格下跌,将降低石油产品支出成本,增加企业的现金流,有利于股价上
涨,这里的石油产品使用行业包括:食品饮料业、餐饮旅游业、纺织服装业、建筑建材业、电子信息业、交运设
备业、机械设备业、家用电器业、轻工制造业、医药生物业、公用事业、交通运输业、金融服务业等。这些行业
在股票市场中所占比重大,国际油价下跌容易产生联动效应,促使股票市场总体价格水平上涨。反之,国际
油价上涨会对上述行业造成相反的影响。特别在石油产品生产成本普遍提高时,极易引起成本推动型通货
膨胀,倘若中央政府为缓解通膨压力,执行紧缩性货币政策,又会带来货币供给下降,整体经济增长减速,最
终导致股市大盘整体下跌。值得关注的是,柳明等的研究表明,考虑到中国巨额的石油进口需求,中国股票
市场的运行状况同样会反作用于国际石油市场,从而实现中国股市对国际油市的反馈[11因此,不排除中
国石油行业对国际石油市场价格信息反馈的可能。虽然绝大多数研究认为,国际油价波动会对中国股票市场价格变动产生影响。但是,由郭国峰和郑召锋
提供的经验证据却显示,国际油价波动对中国股市整体的影响在统计意义上并不显著[12],该研究结论意味
着,股票市场中与石油关系微弱的行业对油价冲击的响应较为有限,用股市综合指数作为国际油价的冲击对
象并不精确。考虑到石油板块股票指数受国际油价波动的冲击最为直接,因此本文根据石油产业链的上下
游关系,选取石油开采、石油加工、石油贸易行业作为国内石油产业链上中下游行业的代表,以检验国际石油
市场与中国石油板块价格信息传递的多重联动关系。四、样本数据和研究方法(一)样本数据本文采用投资型行业分类标准反映中国石油产业链各行业所具有的投资价值,选择申万三级行业指数
中的石油开采(OD)、石油加工(OP)、石油贸易(OT)股票指数,刻画中国石油上下游行业股指的整体变动趋
势,国际石油价格指标采用英国北海布伦特原油期货价格(BRENT)和美国西德克萨斯轻质原油期货价格
(WTI),其中WTI原油期货价格用于稳健性检验。鉴于2009年的石油对外依存度首破国际警戒线50%,
故搜集的样本为从2010年1月1日至2020年6月30日的日度数据。考虑到国际原油价格以美元计价,而
中国申万行业股票指数采用人民币计价,为保证可比性,原油期货价格利用美元兑人民币汇率折算。中国石
油行业股票指数和美元兑人民币汇率数据均来自WIND数据库,国际原油期货价格来自美国能源情报署。
匹配数据时剔除上海证券交易所与纽约商品交易所之间交易日期不一致的交易日,使样本数据的交易日期
保持一致,最后匹配到2 446个数据。为减弱异方差性,对股票指数和石油价格均做了自然对数变换,用
LNOD、LNOP和LNOT表示中国石油开采、石油加工和石油贸易行业股指,用LNBRENT和LNWTI表示
75
统计与信息论坛国际石油价格①。将股票指数对数值的一阶差分记为股指收益率ROD、ROP和ROT,石油价格对数值的一
阶差分记为国际石油价格变动率RBRENT、RWTI。(二)研究方法对于时间序列模型而言,保证估计系数稳定至关重要,如果某个时间序列本身存在结构突变,但未加以
考虑,会大大降低模型的准确度和拟合效果。早期,邹至庄检验能够在一定程度上识别时间序列是否存在结
构变动,但该检验方法要求事先确定转折点位置,并且转折点有且仅有一个。此后,如何准确鉴定时间序列
的突变点,一直是学术界的热点和难点,在此期间,结构突变检验方法经历了一系列改进。Bai等构造的检
验方法,不仅能够有效诊断时间序列是否存在结构突变,还克服了以往研究对突变点个数的限制,目前被认
为是检验结构突变较为客观、准确的方法[13]。
Bai-Perron内生结构突变检验方法的具体思路如下:假设某个5期的时间序列数据存在加个潜在突变点(产生加+
1个分割区间),并且该数据的生成过程
(DGP)如下:$,
=
\')
+
\'/,
+u,(t=5]—1
+1,…,5,
\'
=
1,…,.+
1)
(1)其中,\"为被解释变量,解释变量由\",和两部分组成,\",表示系数未发生改变的变量,-表示系数发生改
变的变量,和&为相应的系数向量,%,为残差项。51,兀,…,丁”是.个未知的结构突变点。为了便于定
义结构突变点的日期,设5=0,5”+1
=
5。Bai-Perron结构突变检验分为三步。第一步,针对式(1)中每个可能的分割点(51,52,-,5„
),利用普
。通最小二乘法(OLS),计算出0和&的估计值,并得到对应的残差平方和:St
(51,5Z,…,5.
)
=
min\'
\'
(y,
—
x)
—
)2第二步,比较不同分割方式得到的残差平方和,取最小残差平方和的分割:
(51,52
,…,5.)
=
argminS5(51,5Z
,…,5.)第三步,对时间序列的生成过程是否发生结构突变进行显著性检验。倘若时间序列数据确实存在结构突变,但建模时却被忽略,得到的实证结果显然是有偏误的。本文首先
基于Bai-Perron内生结构突变检验方法捕捉金融时间序列中可能存在的突变点,将结构突变虚拟变量设为
外生变量,构建包含外生变量的XVAR模型。结构突变虚拟变量用DT,表示,当,<
结构突变点时,DT,=
0,反之,DT,
=
1。与VAR模型相比,包含突变的XVAR模型可以更真实地估计出均值方程,在此基础上提
取残差并做标准化处理,对国际石油价格与中国石油行业股指之间互动关系的估计更为可靠。由前述传导机制分析可知,一方面,国际石油价格波动会对中国股票市场产生影响,另一方面,中国股票
市场同样能够对国际石油价格产生一定反作用,其中,股市中的石油板块受到的冲击和反馈更为明显。分析
跨市场价格波动信息的传导特征,对于分散、防范金融风险具有重要意义。从统计学的角度看,各个市场价
格收益率的方差反映了价格波动的信息,价格波动信息的传导体现在收益率条件的二阶矩中,而多元
GARCH模型正是研究跨市场风险交叉传递的理想计量经济工具。本文应用多元GARCH模型中的
BEKK、CCC、DCC模型对多变量之间的方差协方差矩阵展开分析,这三个模型分别捕捉了价格波动信息在
波动溢出性、静态相关性和动态相关性方面的特征,后续的实证检验包含三个部分,接下来分别介绍BEKK、
CCC、DCC模型的特点及其在本文中的作用:第一部分,应用Engle等提出的MGARCH-BEKK模型,检验国际油价与中国石油行业股指之间的波
动溢出性[14模型的基本表达式为:H,
=
CC
+
+其中,P和g分别表示影响本期残差项条件方差的ARCH项最优滞后期数和GARCH项最优滞后期数,H,
为残差#的条件方差协方差矩阵,C、A,、B,均为\"阶方阵。C为下三角常数矩阵,代表常数项,和B,分①2020年4月20
E],美国西德克萨斯轻质原油期货价格降为一36.
98美元每桶,无法进行对数化处理,本文将其视为异常
值在样本中剔除。76
吕政,胡晨沛:结构突变下国际油价波动与中国石油行业发展的多重联动关系别反映了滞后第i期的ARCH项和滞后第i期的GARCH项对本期条件方差的影响。BEKK模型的优点
在于,保证矩阵H,满足正定条件的同时,较好地刻画了多变量条件方差和协方差交互作用的关系,为解释
各变量之间的联合波动提供了依据,但其经济意义不如单变量GARCH模型+15。借助BEKK模型,可以检
验国际石油价格与中国石油行业股指之间价格波动传导的方向和路径,即BEKK模型用于判断是否存在以
及存在何种波动溢出性,较为全面地反映了跨市场的风险传递过程。第二部分,应用Bollerslev提出的MGARCH-CCC模型,估计国际石油价格与中国石油行业股指之间
不随时间变化的静态相关性[16])关于CCC模型的设定为:H9=D*D9
=
%,
D
=
diag(K11,,i,K12,)
)
(5)(6)h,,,t
—
c,
2
\'l#2,—,
+〉\'Mih,,,—,
8
&(7)其中,H,是\"个时间序列的常条件协方差矩阵*是k_k阶的常系数相关矩阵,D,为是\"阶标准差对角矩
阵是条件方差。与BEKK模型相比,CCC模型的形式更简单,不需要估计太多参数。CCC模型的特点
在于,其估计得到的跨市场价格波动的条件相关系数为常数,因此可以比较整个样本区间内,国际石油价格
与中国石油行业股指之间价格波动传导的相关关系强弱程度,具有较高的经济意义。但是,CCC模型无法
刻画相关系数的时变特征[17第三部分,应用Engle提出的MGARCH-DCC模型,考察国际油价与中国石油上下游行业股指之间随
时间变化的动态相关性[18类似地,H,=DtRtDt是Engle假设的残差条件协方差矩阵表达形式,矩阵H,
中的R,
=
%,,),(”,=
1,为\"阶条件相关系数矩阵,矩阵H中的D,同式(6),是\"阶标准差对角矩阵,条件
方差同式(7)。记标准化的残差向量为%,
=(M1.t,MZ.t,—
%„,),满足ut=D71/#t
,并进一步将动态相关结
构改写成:R,
=
Q;
tQQ-t
其中,Q-是对角元素为槡8二的对角矩阵,Q,为\"阶对称正定矩阵,当滞后期数为1时,Q,形式如下:(8)Qt
=
(1:a:^)Q
+
awt-1M/t-1
+)Qt:1
(9)其中,/为残差的无条件方差矩阵,*和0为动态条件相关系数,*表示随机扰动对相关关系影响的大小,0反
映动态相关关系持续性的大小。根据*和0的显著性水平,判断中国石油开采、加工、贸易行业与国际石油
市场之间的相关性是否具有时变性。与同为检验价格波动相关性的CCC模型相比,DCC模型侧重于刻画
多变量之间随时间变化的相关关系,直观地反映了价格波动相关性的时变特征,但不如CCC模型得到的常
相关系数简单明了。在现实经济中,对于多变量价格波动的常相关约束可能并不成立[19]。例如,个别时间
点上的突发风险事件会对国际油价带来巨大冲击,在该冲击下,国际油市和中国股市石油板块之间的相关性
极可能发生变动。因为金融时间序列的相关性信息能够较好的被1阶的波动相关性所解释,所以在实证检
验部分本文选取BEKK(1,1)、CCC(1,1)和DCC(1,1)模型进行多元波动率建模。五、实证检验(-)Bai-Perron内生结构突变检验判断时间序列是否存在结构突变,主要观察Bai-Perron检验方法提供的SupFT(k)统计量和Dmax统
计量,进一步借助SupFT
(L
+
1
|F)统计量估计突变点个数。对BRENT、WTI、OD、OP和OT序列的
Bai-Perron检验结果见表1。由表1可知,仅有BRENT和WTI序列的SupFT
(1)统计量以及更高阶的
SupFr(k)统计量均显著,同时相应的UDmax统计量和WDmax统计量也在1
%显著性水平下显著,意味着
BRENT和WTI序列均存在结构突变。关于OD、OP和OT序列,虽然更高阶的SupFT(k)统计量.UDmax
统计量和WDmax统计量均在1
%显著性水平下显著,拒绝不存在结构突变的原假设,但是对应的SupFr(1)
统计量并不显著,暗示OD、OP和OT序列是否存在结构突变有待进一步考察,经检验这3个时间序列并不
存在结构突变点。检验结果表明,BRENT序列和WTI序列均发生了
1次结构突变,结构突变点较为接近,
分别为2014年11月4日和2014年11月13日。77
统计与信息论坛表1Bai-Perron检验统计量Bai-Perron内生结构突变检验BRENT25.
710---26.
260---31.
143---相关统计量及突变点WTI序列55.612---OD0.0365.91727.422---OPOT2.7596.539-10.759---10.292---SupFT(k)统计量SupFr
(1)SupFT
⑵SupF^
(3)SupFT
(4)SupFT
(5)2.2475.34311.673
---9.329---69.910---47.418---32.542---39.362---60.
616---44.
943---13.003---21.917---27.422---8.147---4.634--10.759--Dmax统计量UDmax60.
616---120.
351***69.910---98.525---55.612---14.983
--20.392---WDmaxSupFT(L+1
|L)统计量SupFT
(1
0)SupFT(2|1)突变点54.858
---20.434---25.710---0.0362.2472.7591.7782014/11/041.75312014/11/13注!---分别表示在1%、5%、10%的水平下显著)下同)检验得到的结构突变点具有深厚的政治经济背景。国际油价的两个代理变量在样本期内的结构突变点
分别位于2014年11月13日和2014年11月4日,事实上,2014年下半年,国际油价大幅度下跌。从需求层
面看,当年全球经济恢复远低于预期,加之中国经济增速放缓,使得全球对石油的总需求下降。供给层面的
冲击主要来自于两个方面:一方面,美国以页岩油气为代表的新能源繁荣兴起,打压了石油价格;另一方面,
OPEC成员国的地缘政治危机得到一定程度缓解,逐渐恢复原油正常供应。在全球石油需求量不足的情况
下,供给量却在不断增加,供大于求是石油价格暴跌的根本原因。该突变时点与经济现实吻合,说明结构突
变检验结果是可信的。(二)构建均值方程XVAR模型1.单位根检验为避免建模过程中出现伪回归,需要对各变量做单位根检验,经Bai-Perron内生结构突变检验可知,
BRENT和WTI序列存在结构突变,因此本文对这两个时间序列及其对数收益率序列进行突变点单位根检
验,采用AO检验方法,该方法假设突变的发生是即时的,瞬间的)经检验,BRENT的对数序列LNBRENT
接受零假设,是一个具有结构突变点的非平稳序列,一阶差分后的对数收益率序列RBRENT拒绝零假设,
表明LNBRENT为1阶单整序列,记为I(1),
WTI的对数序列LNWTI及其对数收益率序列RWTI均拒绝
存在单位根的零假设,属于具有结构突变点的平稳序列。未发生结构突变的国内石油行业股指序列采用一
般的单位根检验方法,本文ADF检验的滞后项阶数由SC准则确定。考虑到样本容量和滞后阶数的变化都
会在一定程度上影响ADF检验结果,提高犯第二类错误的可能,所以补充了
PP检验和KPSS检验作参考。
经检验,三种单位根检验的结果完全一致,LNOD、LNOP和LNOT序列均为1阶单整序列I(1)。2
Johansen
检验考虑到国际油价和国内石油行业股指4个时间序列具有同阶单整关系,在该情形下,单整变量的线性组
合很可能是平稳的,不排除多个时间序列之间存在长期均衡关系的可能。因此有必要进行协整检验确认长
期均衡关系是否存在,这将为后文选择向量自回0模型(VAR)或向量误差修正模型(VECM)作为均值方程
提供重要依据。需要说明,协整检验的滞后阶数为VAR模型的最优滞后阶数减1。本文综合比较了多项评
价指标,最终确定VAR模型的最优滞后阶数为5阶,协整检验的滞后阶数为4阶。Johansen协整检验结果
显示,样本区间内国际油价和国内石油行业股指之间不存在长期的协整关系,因此不适合建立VECM模型。
下面就LNBRENT、LNOD、LNOP和LNOT变量的差分序列RBRENT、ROD、ROP和ROT进行格兰杰因
果检验,进一步考察四者之间的因果关系。3.格兰杰因果关系检验国际油价变动率、中国石油开采、石油加工、石油贸易行业股指收益率的格兰杰因果关系检验结果见
表2。表2显示,从选取的样本期看,在5
%的显著性水平下,国际油价波动是国内石油开采、石油加工、石油贸
易行业股指收益变动的格兰杰原因。此外,国内石油开采、石油贸易行业股指收益变动和国际油价波动互78
吕政,胡晨沛:结构突变下国际油价波动与中国石油行业发展的多重联动关系为因果关系,石油开采行业股指收益变动是石油加工、石油
贸易行业的单向格兰杰原因,石油加工行业股指收益变动
表2格兰杰因果关系检验零假设F统计量2.9028.2561.5995.301P值0.0130.0000.1570.0000.0090.0110.645是石油贸易行业的单向格兰杰原因。从检验结果看,国际
油价波动作为中国石油板块股价变动的格兰杰原因,基本
ROD不是RBRENT的格兰杰原因RBRENT不是ROD的格兰杰原因ROP不是RBRENT的格兰杰原因符合经济事实)中国石油行业股指收益率对国际油价变动
率的因果关系,仅石油加工行业在统计意义上不显著,反映
RBRENT不是ROP的格兰杰原因ROT不是RBRENT的格兰杰原因RBRENT不是ROT的格兰杰原因ROP不是ROD的格兰杰原因出现阶段中国石油加工行业尚不能对国际油价造成影响,
而石油加工、石油贸易行业股指收益变动能够在一定程度
上反作用于国际油价)中国上游石油开采行业股指收益变
3.1002.9670.6722.4291.711ROD不是ROP的格兰杰原因0.0330.1290.001动是中下游石油加工、石油贸易行业的单向格兰杰原因,中
ROT不是ROD的格兰杰原因ROD不是ROT的格兰杰原因游石油加工行业股指收益变动是下游石油贸易行业的单向
格兰杰原因,也符合石油产业链从上至下的价格信息传递
规律。4.提取均值方程残差4.3101.117ROT不是ROP的格兰杰原因ROP不是ROT的格兰杰原因0.3490.0212.648注:选取的滞后阶数为5阶)VAR模型抛开了经验意义上的因果判断,从样本数据自身出发,寻找变量间的互动关系。因此,本文参
考白万平等的做法,将VAR模型设为均值方程[20]。为更好地反映结构突变带来的经济影响,本文在VAR
模型的基础上引入外生的结构突变虚拟变量,进一步构建XVAR模型。多元波动率建模的前提是提取均值
方程的残差,因此本文对残差进行了一系列检验。经检验,4个残差序列的偏度均小于0,峰度均大于3JB
检验的P值几乎为0,说明4个序列均不服从正态分布。ROP和ROT残差序列Q(8)、Q(16)和Q(24)的P
值都较大,接受残差#不存在自相关的原假设,而RBRENT、ROD残差序列的Q(8)、Q(16)和Q(24)统计量
分别在不同程度的水平下显著,说明残差存在一定程度的自相关。4个残差平方序列Qz
(8),QZ
(16)和
Qz(24)的P值均小于1%,表明均值方程的残差平方序列#2存在显著的自相关。综合残差项的LM检验和
ARCH检验结果,RBRENT、ROD的残差序列存在一定程度的自相关,而ROP和ROT的残差序列不存在
自相关,以上4个残差序列均存在明显的ARCH效应,该检验结果说明,接下来构建的多元GARCH模型
是合理的。(三)波动溢出性估计以往文献在运用BEKK模型时,未考虑序列本身存在的结构突变,提取VAR模型标准化后的残差直接
建模,可能导致实证结果失真。对此,本文尝试改进,将结构突变的影响纳入模型,以增加实证结果的可靠
性。根据前述分析,在波动溢出性检验部分,将国际油价变动率、中国石油开采、石油加工和石油贸易股指收
益率的多元波动率模型设为MGARCH-BEKK(1,1),国际油价与中国石油行业的波动溢出性分析基于A、
B矩阵元素的估计结果,C为常数矩阵,模型估计结果见表3。表3中A矩阵的主对角线元素s代表#序列
自身滞后1期残差平方对自身本期条件方差的影响,非主对角线元素5代表#序列滞后1期残差平方对j
序列本期收益率条件方差的影响。同理,B矩阵中的元素反映了条件异方差的影响。从表3可以看出,A.B
矩阵的对角线元素绝大多数在1%的显著性水平下显著,说明4个标准化残差序列具有明显的条件异方差
特征,即油市和股市的波动均存在明显的聚集性,当期的波动会随滞后一期的扰动和波动的增大而增大。观
察表3中非主对角元素可知,除了部分系数不显著外,绝大多数非对角元素上的系数均在1%的显著性水平
下显著。BEKK模型的估计结果初步表明,价格波动的信号在国际石油价格和中国石油行业股指之间具有
传递性,但考虑到BEKK模型估计系数的经济意义不强,需要进一步给出BEKK模型波动溢出效应的Wald
检验结果。BEKK模型的Wald检验结果显示了
4个时间序列波动溢出的方向,检验结果见表4。由表4可知,国
际油价与石油开采行业的联合检验表明,在1%的显著性水平下,拒绝原假设L15=M”=0、L51=M51=0,意味
着国际油价与中国石油开采行业存在双向波动溢出效应。国际油价与石油加工行业的联合检验表明,在
5%的显著性水平下,拒绝原假设L13=Mi3=0、31=仇1=0,意味着国际油价与中国石油加工行业存在双向波
动溢出效应。国际油价与石油贸易行业的联合检验表明,在统计意义上不拒绝原假设L14=Mm=0、S1=仇1
79
统计与信息论坛=0,意味着国际油价与中国石油贸易行业之间不存在显著的波动溢出效应。类似地,容易得到中国石油行
业上下游各环节之间存在单向波动溢出效应的结论。上述结果表明,国际油价对中国石油行业的上游、中游
行业存在较为明显的波动溢出效应,且中国石油开采、石油加工行业股指收益变化会引起国际油价波动,但
石油贸易行业股指收益变化对国际油价的影响并不明显)表3波动溢出效应模型估计结果待估矩阵0.050**-(0.
011)参数矩阵:0.143***
(0.
013):0.069***
(0.011):0.056***
(0.019):0.040***
(0.014)0.036** (0.017)0.002(0.012)0.005(0.010):0.053*
(0.032):0.032***
(0.011)0.058***
(0.013)0.166***
(0.017):0.085*** (0.027):0.038*** (0.016)0.339*** (0.019)0.044** (0.022):0.013(0.012)0.004(0.014)A:0.070***
(0.013)0.016(0.014)0.015(0.012)0.993***
(0.001)0.032***
(0.004)0.013(0.021):0.014(0.022)0.216*** (0.020):0.056*** (0.015)0.006*** (0.002)0.921*** (0.007)0.018(0.013)0.001(0.001):0.022***
(0.003)0.992***
(0.004):0.001(0.001):0.006* (0.003)0.017*** (0.006)0.961*** (0.006):0.008***
(0.003)0.007(0.005)0.017*** (0.004)0.00(0.002)LogLikelihood注:小括号内为标准课)下同):0.005*
(0.003):10448.263(四)静态相关性估计基于标准化残差建立MGARCH-
CCC(1,1)模型,以考察国际油价波动和
波动溢出关系表4波动溢出效应检验Wald检验中国石油行业股指变动之间的常相关关
RBRENT对ROD的检验国际油价与石油开采H
:
L12
=
&12
=
0
02
ROD
RBRENT
的
验系,估计结果见表5。由表5可知,4个
时间序列两两之间的常相关系数均在
1%的统计意义上显著。具体地,国际油
H0
:21
=
&21=
0
02
=
20.
423***=
89.
940
***国
油价与
油RBRENT
ROP
的
验H0
!13=M13=002
ROP
RBRENT
的
验H0
!31=仇1=002
=
8.
948**价和中国
油开
、
油
、
油
行业的相关系数分别为0.
128、0.
131、
0.108(
油开
行业和
油
、
油国
油价与
油=
15.
485
***RBRENT
ROT
的
验H0
:
L14
=
&14
=
002
ROT
RBRENT
的
验
H0
L41=M41=0
02
=
1
694=
2.
431行
业
的
关
数
分
0.724、
0.524(
油
行业和
油
行业的相关系数为0.683。该结果表明,国内
ROD
ROP
的
验ROP
ROD
的
验油开
与
油H0
L23=M23=0
02
H0
L35=M35=002
=
46.
029***=
3
697石油产业链之间的相关性较国际油价与
国内石油行业的相关性更高,基本符合
经济现实。从CCC模型关于静态相关
的
计
(国
油价石油加工与石油贸易油开
与
油ROD
ROT
的
验日0
:
L
=
&24
=
002
ROT
ROD
的
验H0
L45=M45=0
=
3
67302 =
24.
396***ROT
ROD
的
验ROD
ROT
的
验H0
:
&34
=加4
=
002
H0
L43=M43=005与中国
油
业链
游、
中
游的
油开、
油
行业
指
益变
的
关
性较强,达到0.13左右,而国际石油价
RWTIRWTIROD=
11
660=3.595表5静态相关关系模型估计结果RODROPROT与中国
油
业链下游的
油
行业
指
益变
的
关
10.128***
(0.020)10.724***
(0.010)0.524***
(0.014)弱,大致为0.10。因此,位于中国石油
行业上游、中游的石油上市公司需要密切ROP0.131
***
(0.020)0.108***
(0.020)10.683*** (0.011)ROT1关注国际油价走势,尽可能降低国际油价冲击风险。80
吕
政,胡晨沛:结构突变下国际油价波动与中国石油行业发展的多重联动关系(五)动态相关性估计建
准化残差项的四元MGARCH-DCC(1,1)模型,研究国际石油价格与中国石油行业的动态关联
性。国际油价、石油开采、石油加工和石油贸易行业之间的相关性是否随时间变化而变化取决于DCC模型a和仔系数的显著性(度量随机
表6动态相关关系模型估计结果MGARCH部分参数估计扰
关关
的大小,仔i=
16.664--(0.002)6.651--(.008)0.
947
--(.008)反映
关关
续性的长
i
=26.621---(0.
004)6.696---(0.
010)i
=30.
011---(0
.
0
0
3)0.
0
53
---(0
.
0
0
6)i=40.
019---(0
.
0
0 4)短,估计
表6。从表6可,与)的估计值均
异于
0,且两者之和小于1,符合DCC
0
.
052---(0
.
0
0
7)0.929
---(0
.
0
10
)0.
882---(0.
012)0
.
935
---(0.
007)模型的约束条件。)=0.
973接于1,表明
参数
DCC部分参数估计当期
?Log Likelihood0.
021
---
(0.
003)0.
973---
(0.
004)异方
要来源于上一期的波—16436.612准
间的相关系数具有极强的持续性。油行业股指收益率间的动态条件相关性,将动态的相关系数绘制成图1。从图1可以看
出,自2010年以来,国
市场的供求关
需求
油价格和中国石油行业股指之间的价
形势,以及
信号传
明显的联动关系。,导
界石油2014年年底,受中美欧等主要经济
亚和埃及等中东
油供给增加,新
不平衡,国际油价与石油开采行业的相关系数骤减,这也和检验
的突变时间点
J。近期,由于俄
与OPEC组织的减
,最终导
经济运行
油价格暴跌,国
作破裂引
来全球原油增加,股市中能源
油市场振荡。由于市场不确
的投资者避险
增强,转向配置安全 高的
如
等,投资者的 于
行为,在
冲击时,石油行业股指收益率的相关
下降。值得关注的是,在2020年初国际油关性均大幅下滑,但国际油价波动与国
油于石油开采行业更大,说明石油产业链不同行业2010-01-12/2020-06-30价暴跌的时点上,国
、石油
油价格和国
油行业股指的
行业股指收益变动的相关系数下
对油价冲击的响应具有异质性。0.60.40.20-0.2国际石油价格与国内石油开采行业股指的动态相关系数
1
月
12(2010)
1
月
04(2011)
1
月
04(2012)
1
月
04(2013)
1
月
02(2014)
1月
05(2015)
1
月
04(2016)
1
月
03(2017)
1
月
02(2018)
1
月
02(2019)
1
月
02(2020)0.6
-国际石油价格与国内石油加工行业股指的动态相关系数2010-01-12/2020-06-30-0.2
-1
月
12(2010)
1
月
04(2011)
1
月
04(2012)
1
月
04(2013)
1
月
02(2014)
1
月
05(2015)
1
月
04(2016)
1月03(2017)
1
月
02(2018)
1
月
02(2019)
1
月
02(2020)国际石油价格与国内石油贸易行业股指的动态相关系数0.40.22010-01-12/2020-06-3001
月
12(2010)
1
月
04(2011)
1
月
04(2012)
1
月
04(2013)
1
月
02(2014)
1
月
05(2015)
1
月
04(2016)
1月03(2017)
1
月
02(2018)
1
月
02(2019)
1
月
02(2020)图1国际油价与中国石油行业动态条件相关性走势(六)稳健性检验国
(WTI)
重要的原油期货合约以英国北海布伦特原油期货(BRENT)和美国西德克萨斯轻质原油期货,两种期
间的关联
强。英国北海布伦特原油期
的
量约占
81
统计与信息论坛市场总量的60%,是世界上流动性最好的期货品种。美国西德克萨斯轻质原油价格是全球原油定价的基
准,其价格走势是其他原油价格的风向标。为说明本文实证检验结果的可靠程度,本文将WTI油价作为国
际油价的代理变量,对国际石油市场和中国石油行业的互动关系进行稳健性检验。得到的结论如下:一是WTI油价的结构突变点在2014年11月13日,与BRENT油价的突变点较为接
近,说明两者具有相似的波动轨迹;二是WTI油价变动率同样与中国石油开采、石油贸易行业股指收益率存
在双向因果关系;三是WTI油价与中国石油行业之间同样存在明显的波动溢出效应;四是得到了
WTI油
价、中国石油开采、石油加工和石油贸易行业之间的静态和动态相关系数,国际油价和国内石油行业的动态
相关系数虽然不高,但大体呈正相关变动关系。可以认为,国际石油市场和中国石油行业之间的波动溢出
性、静态相关性和动态相关性在不同代理变量下是稳健的。六、结论本文以结构突变为视角,对2010年1月1日至2020年6月30日国际石油价格与中国石油行业之间的
多重联动关系提供了新的经验证据。选取英国北海布伦特原油期货价格为国际原油市场代表,将中国石油
产业链上中下游行业在资本市场上的代理变量设为石油开采、石油加工和石油贸易行业股票指数,并以西德
克萨斯轻质原油期货价格作为国际油价的其他度量进行了稳健性检验。Bai-Perron内生结构突变检验结果表明,英国北海布伦特原油期货价格和西德克萨斯轻质原油期货价
格的结构突变点均在2014年11月,国际原油供大于求是当时油价暴跌的主要原因。从统计因果关系看,国
际油价波动是国内石油开采、石油加工、石油贸易行业收益率变动的单向格兰杰原因,意味着国际油价波动
对中国石油行业经营绩效起先导作用。中国上游石油开采行业是中下游石油加工、石油贸易行业的单向格
兰杰原因,说明石油行业具有成本接近和产品同质的特征,业绩易受共同因素冲击,股票价格走势趋同,存在
较明显的板块联动效应。将结构突变因素引入均值方程,在此基础上建立多元波动率模型,能够更加准确地
捕捉到国际油价与中国石油行业的多重关系。BEKK模型的估计结果表明,中国石油开采、石油加工行业
与国际油价具有双向的波动溢出效应。综合CCC模型和DCC模型的估计结果可以看出,中国石油开采、石
油加工、石油贸易行业股指收益率与国际油价变动率之间大致保持正相关关系,且具有明显的时变特征,该
变动与国际油价波动具有紧密联系。国际石油价格高位跳水对石油行业的运营带来了潜在隐患,中国石油行业由最具代表性的中石油、中石
化和中海油三大国有企业主导,如何规避油价波动风险是一个不得不面对的重要问题。国际石油市场供大
于求的局面不可逆转,世界能源市场逐渐转变为买方市场,中国作为发展中的世界第二大经济体,应抓住有
利的窗口期,在国际石油贸易谈判中应抓牢主动权,积极提升买方市场的话语权。同时,中国作为世界最大
的石油净进口国,巨额石油需求对全球石油定价有举足轻重的影响,中国应积极推进国际石油交易使用人民
币计价与结算,不仅有利于提高人民币补充全球流动性的潜力,更有助于扩大中国的国际影响力。参考文献:[1,
Broadstock
D
C,Cao
H,Zhang
D.
Oil
Shocks
and
Their
Impact
on
Energy
Related
Stocks
in
China
[J,.
EnergyEconomics,2012,34(6)
:
18881895.[2,
Hamilton
J
D.
Oil
and
the
Macroeconomy
Since
World
War
II
[J,
Journal
of
Political
Economy,1983,91(2)
:228-248.
[3,
Ferson
W
E,
Harvey
C
R.
Sources
of
Risk
and
Expected
Returns
in
Global
Equity
Markets
[J,.
Journal
of
Banking
&Finance,1994,18(4):775-803[4,
Park
J
,Ratti
R
A.
Oil
Price
Shocks
and
Stock
Markets in
the
US
and
13
European
Countries
[J,.
Energy
Economics,2008,30(5):2587-2608[5,
Nandha
M,
Faff
R.
Does
Oil
Move
Equity
Prices?
A
Global
View
[J,
Energy
Economics
,008,30(3):986-997.许金华,范英•金融危机对石油市场定价机制的影响分析[J,数理统计与管理,2014,33(5):790801.杨胜刚,成博.股票市场与大宗商品市场互动特征比较研究当代财经,2014(6)5766.李素芳,朱慧明,李荣.基于贝叶斯机制转换协整模型的石油一#股市非对称效应研究中国管理科学2015,23(9):
46-5482
吕政,胡晨沛:结构突变下国际油价波动与中国石油行业发展的多重联动关系+
,王朝阳,陈宇峰,金曦.国际油价对中国新能源市场的传导效应研究数量经济技术经济研究,2018,35(4):131-146.
[10,韩立岩,甄贞,蔡立新.国际油价的长短期影响因素+,中国管理科学,2017,25(8):687
8+
1,柳明,杨运泽,孙文鑫.成品油调价对能源类股票收益率的影响一#基于事件分析的方法+,经济评论,2013(3):
109-117[12,郭国峰,郑召锋.国际能源价格波动对中国股市的影响#—基于计量模型的实证检验+,中国工业经济,2011(6):26-35[13]
Bai
J,Perron
P.
Computation
and
Analysis
of
Multiple
Structural
Change
Models
[J,.
Journal
of
Applied
Econometrics,
2003
18(1):1-22[14,
Engle
R
F,Kroner
K
F.
Multivariate
Simultaneous
Generalized
ARCH
[J,
Econometric
Theory,1995,11(1)
:
122-150.
+
5,雷钦礼,陈晓蒙.中、印、美股市联动差异性研究统计与信息论坛,2017,32(4):4150.[16,
Bollerslev
T.
Modelling
the
Coherence
in
Short-run
Nominal
Exchange
Rates:
A
Multivariate
Generalized
ARCH
Model
[J,
TheReviewofEconomicsandStatistics,1990,72(3)
498-505[7,刘丽萍,马丹,唐晓彬.基于高维数据的改进CCC-GARCH模型的估计及应用统计与信息论坛,2016,31(9):2228
[18,
Engle
R
Dynamic
Conditional
Correlation:
A
Simple
Class
of
Multivariate
Generalized
Autoregressive
Conditional
Heteroskedasticity
Models
[J,. Journal
of
Business
&
Economic
Statistics,
2002,20(3)
:
339-350.[19,唐晓彬,董曼茹,乔天立,等.中国资本市场资产价格波动的动态关联性检验[J,统计与信息论坛,2020,35(1):5363.
[0,白万平,吕政,刘丽萍.石油产业链股指收益率动态关系研究一一基于石油开采、加工和贸易行业数据分析价格理论
与实践,2018(7):
Mutual
Linkage
between
International
Oil
Price
Fluctuation
and
the
Development
of
Chinese
Oil
Industry
under
the
Structural
ChangeLYU
Zheng1
,
HU
Chen-pei2(1.
School
of
Statistics
and
Mathematics,Central
University
of
Finance
and
Economics,Beijing
102206
,China;
2InternationalStatisticalInformationCenter,NationalBureauofStatistics,Beijing100826,China)Abstract:The
sudden
change
in
international
oil
prices
and
the
rising
dependence
of
China\'s
crude
oil
on
foreign
sources
have
increased
the
complexity
of
the
international
oil
market
?s connection
with
the
Chinese
oil
industry
and
clarifying
the
relationship
between
the
two
will
help
enhance
China\'s
voice
in
the
end,thetransmission
mechanism
betweenfluctuationsininternationaloil
prices
and
changes
in
China\'s
oil
industry
stock
indexes
has
been
clarified.
In
the
empirical
test
part,
the
Bai-Perron
endogenous
structural
mutation
test
method
is
used
to
diagnose
the
mutation point,
and
the
structural
mutation dummy
variable
is
set
as
an
exogenous
variable
in
the
mean
equation
XVAR
model
,nd
amultivariate
GARCH
model
such
as
BEEK,
CCC,
DCC,tructed
on
this
basis,
provide
new
empiricalevidenceforthemultiplerelationshipsbetweeninternationaloilpricefluctuationsandChina\'soil
exploration,oilprocessing,andoiltradingindustrystockindexreturns,suchasvolatilityspilovers,static
correlations,earcharrivesatthefolowingconclusions:Internationaloil
pricefluctuationsplayaleadingroleinthechangesin
China\'sa
two-wayvolatilityspiloverefectbetweenChina\'soilexplorationandoilprocessingindustrystockindex
amiccorrelationcoeficientsbetweentherateof
returnofthestockindexintheChina\'soilextraction,oilprocessing,andoiltradingindustriesandtherate
ofchangeininternationaloilpricesalshow econd-largest
economyinthedeveloping
world,Chinashouldseizethefavorable
window
period
wheninternationaloil
pricesarehoveringatalowleveltoenhancethebuyer\'ds:internationaloilpricespetroleumindustry
volatilityspilovertime-varying\'multivariate
GARCH(责任编辑:于茜)83
更多推荐
石油,行业,中国,国际,油价,波动
发布评论