1、基于马尔科夫状态转移模型的天气衍生品定价杨刚 杨徐摘 要 引入马尔科夫状态转移（MRS）模型拟合长沙市每日平均气温变化，利用最大期望算法估计马尔科夫状态转移模型参数，通过误差分析得到了最佳 MRS 模型.基于最佳的 MRS 模型，采用无套利定价原理定价气温衍生品，并利用蒙特卡罗方法得到了取暖指数（HDD）欧式看涨期权的数值解.实证结果表明，五状态的 MRS 模型对长沙市每日平均气温变化的拟合效果明显优于其他的 MRS 模型，它使得气温衍生品定价结果相比以前的方法更为精确.关键词 数量经济学;最优气温模型;EM 算法;无套利定价原理;天气衍生品;MRS 模型中图分类号 F840.67文献标识码
2、AAbstract The Markov regime-switching model is introduced to fit thevariation of Changsha daily average temperature.The parameters of Markovregime-switching（MRS）models are estimated by Expectation-Maximizationalgorithm.The optimal MRS model is derived by error analysis.On the basis ofthe optimal MRS
3、 model，no-arbitrage pricing principle is applied to pricingthe temperature derivatives.Numerical solution to the HDD European calloption is determined by Monte Carlo method.As the empirical result is shown，the fitting effect of five-state-MRS model for the variation of Changsha dailyaverage temperat
4、ure is superior to other MRS models.Hence，the pricingresult of temperature derivatives is more accurate than ever.Key words quantitative economics;optimal temperature model;EM algorithm;no-arbitrage pricing principle;weather derivatives;MRS model1 引 言自 1997 年 12 月 84 個国家在日本京都签订京都协议书以来，人们日益认识到气候变化对生产
5、生活的严重影响，并且把该年定为天气衍生品交易元年.2014 年 8 月，国务院发布了关于加快发展现代保险服务业的若干意见，指出要“探索天气指数等新兴产品和服务”.天气风险通常分为巨灾天气风险和非巨灾天气风险.对于巨灾风险，许多学者已进行了深入研究.杨刚等（2008）1通过 Esscher 变换对巨灾超额损失再保险进行定价，沈明轩和何朝林（2012）2假定巨灾指数服从分数跳扩散情形，利用保险精算方法定价巨灾期权.天气衍生品是一种基于天气指数的未定权益，它通常用于管理因不利天气变化造成的经济损失，常见标的包括温度、降雨量、降雪量和霜冻天数等天气指标.由于天气衍生品的标的不可交易，天气衍生品市场属于
6、不完全市场，因此不能应用经典的 Black-Scholes 期权定价方法对天气衍生品进行定价.现有的天气衍生品定价方法种类多样，主要包含精算定价方法、无套利定价方法、效用无差异定价方法和均衡定价方法（Lee 和 Oren，（2010）3）.传统的精算定价方法通常基于精算等价原则，往往采用真实概率测度得到贴现价格，由于天气衍生品市场存在显著的风险市场价格，故容易导致较大的定价误差.效用无差异定价方法和均衡定价方法通常以效用函数为基础，由于对个人效用函数的选取较为主观，所以容易产生较大的定价误差.选取风险中性定价方法确定气温衍生品的无套利价格比较合理，它是一种灵活而又便于实施的金融衍生品定价方法.
7、在无套利定价方法中，Alaton 等（2002）4首次引入布朗运动驱动的 O-U 过程对气温衍生品进行无套利定价.王明亮等（2015）5利用布朗运动驱动的 O-U 过程拟合北京市每日平均气温变化，通过蒙特卡洛方法定价气温衍生品.王晶（2016）6利用布朗运动驱动的 O-U 过程拟合北京市每日平均气温变化，通过有限差分方法得到气温期权的数值结果.胡亚茹（2019）7采用布朗运动驱动的 O-U 过程拟合哈尔滨、北京等五个城市的每日平均气温变化，通过蒙特卡洛模拟得到气温期权价格.Alaton 等（2002）4采用常数方差刻画气温的波动，Benth 等（2007）8采用截断的傅里叶级数描述方差变化的时
8、变结果，Li（2018）9对布朗运动驱动的 O-U 模型进行改进，采用偏微分方程（PDE）定价天气衍生品.Benth 等（2005）10发现气温残差并不完全符合正态分布的假定，这使得气温衍生品在定价过程中存在较大的误差，所以他们采用 lvy 过程驱动的 O-U 过程来拟合气温过程的变化轨迹.Elias 等（2014）11首次将两状态马尔科夫状态转移模型（Markov regime-switching，MRS）模型应用于气温衍生品定价.随后，Evarest 等（2017）12对两状态MRS 模型进行拓展，提出时变波动率 MRS 模型.Gyamerah 等（2018）13在 Evarest 等（2
9、017）12的工作上对 MRS 模型进一步修正，得到 lvy 过程驱动的两状态 MRS 模型.Xiong 和 Mamon（2018）14采用受隐马尔可夫链调制的 O-U 过程对天气衍生品进行定价.MRS 模型在一定程度上提升了天气衍生品定价的精确性，同时也增加了计算的复杂性.通过对不同气温模型对比分析，最终采用 MRS 模型对气温衍生品进行无套利定价.受 Evarest 等（2017）12和 Gyamerah 等（2018）13提出的两状态 MRS 模型启发，进行了两方面的拓展，一是将时变波动率均值回复 MRS 模型中基本状态的时变波动率拓展到带有常数指数时变波动率的 MRS 模型;二是将两状
10、态的 MRS 模型拓展到五状态的 MRS模型.同时，基于最优的 MRS 模型对取暖指数（heating degree days，HDD）期权进行无套利定价，并通过蒙特卡罗方法得到了 HDD 期权价格的数值解.2 每日平均气温分解和马尔科夫状态转移模型通过分解得到每日平均气温的残差，并用 MRS 模型拟合每日平均气温残差的变化.3 气温模型数值分析首先描述數据来源和数据特征，然后对方程（2）的参数使用非线性最小二乘法估计.最大期望算法（Expectation-Maximization algorithm，EM）以极大似然估计方法为基础，较为简便也便于操作，算法收敛的稳定性强，采用 EM 算法对
11、MRS 模型的参数进行估计应该是可行的.基于所有的 MRS 模型，拟合气温数据并使用误差分析，挑选最优的 MRS模型以便于下一节对气温衍生品进行定价.3.1 数据来源与数据特征描述气温数据来源于美国海洋和大气管理局（NOAA）美国海洋和大气管理局官网（https：/www.ncdc.noaa.gov/data-access/quick-links#dsi-3505）.将气温衍生品应用于农业天气风险管理之前，需要对气温衍生品进行定价.由于气温衍生品用于管理农业天气风险可能存在地理基差风险，而减少基差风险的简便方式是选取覆盖地点较近的气温衍生品合约.因此选取长沙站点作为湖南省具有代表性的气象站点，
12、时间区间为 2016 年 1 月 1 日到 2019 年 12 月 31 日，样本数量为 1461 个.在此研究的气温衍生品是以长沙市每日平均气温数据为标的的气温衍生品.需要注意的是，基于长沙市每日平均气温指数的气温衍生品的价格数值并不一定与采用其他地区的气温数据得到的价格数值一致，原因是不同地区的气温数据存在差异，但仍然可以采用与此相同的模型、方法和步骤得到其他地区的气温衍生品的价格.由于样本量较大，可能存在缺失数据的情形.对缺失数据的处理，详见 Alexandridi 和 Zapranis（2012）16.表 1 中描述性统计呈现长沙每日平均气温特征，表明样本数据不服从正态分布，并且波动性
13、较大.3.2 确定性函数的参数估计采用非线性最小二乘法得到方程（2）的最佳参数拟合值见表 2.R2 为 84.79%，表明除e3 以外确定性成分平均解释了 84.79%的气温变化 e3 的 P 值很大而参数值很小，所以将其剔除.在图 1 中，拟合的确定性成分总体上与真实数据相吻合.在图 2 中，气温残差有较为显著的均值回复特征，残差序列的波动性呈现强烈的“跳跃”性，表明适合引入 MRS 模型拟合气温残差.5 结 论通过将 Evarest 等（2017）12和 Gyamerah 等（2018）13MRS 模型的时变波动率拓展到带常数指数的时变波动率 MRS 模型，将两状态的 MRS 模型拓展到五
14、状态的 MRS 模型，发现 MRS7 模型对长沙市每日平均气温残差的拟合效果相对较好.五状态 MRS 模型优于其他的 MRS 模型.在气温期权定价的数值结果中，存在气温风险市场价格和不存在气温风险市场价格两种情形下得到的期权价格相差较大，在气温期权定价时，需要考虑风险市场价格，这与 Xiong 和 Mamon（2018）14的研究结果一致.在高状态 MRS 模型的拓展中，并未使用 lvy 过程作为“跳跃”状态构建 MRS 模型，这可能使得某些极端异常气温难以被模型拟合，未来可以构建常数指数时变波动率均值回复过程为基本状态、lvy 过程为“跳跃”状态的 MRS 模型.参考文献1 杨刚，刘再明，欧
15、阳资生，等.隐含风险中性分布在巨灾超额损失再保险定价中的应用J.经济数学，2008，25（4）：331-337.2 沈明轩，何朝林.分数跳扩散环境下的巨灾期权定价J.经济数学，2012，29（3）：78-81.3 LEE Y，OREN S S.A multi-period equilibrium pricing model of weatherderivativesJ，Energy Systems，2010，1（1）：3-30.4 ALATON P，DJEHICHE B，STILLBERGER D.On modelling and pricingweather derivativesJ.App
16、lied Mathematical Finance，2002，9（1）：1-20.5 王明亮，何建敏，陈百硕，等.时变 O-U 模型在气温预测及气温期货定价中的适应性研究基于北京市 1951-2012 年的日平均气温数据J.中国管理科学，2015，23（2）：44-49.6 王晶.天气衍生品定价模型的构建J.统计与决策，2016，32（15）：152-155.7 胡亚茹.基于 O-U 模型的气温衍生品定价研究J.武汉金融，2019，36（4）：31-35，71.8 BENTH F E，ALTYT-BENTH J.The volatility of temperature and pricing
17、of weather derivativesJ.Quantitative Finance，2007，7（5）：553-561.9 LI P.Pricing weather derivatives with partial differential equationsof the Ornstein-Uhlenbeck processJ.Computers&Mathematics withApplications，2018，75（3）：1044-1059.10BENTH F E，ALTYT-BENTH J.Stochastic modelling of temperaturevariations
18、with a view towards weather derivativesJ.Applied MathematicalFinance，2005，12（1）：53-85.11ELIAS R S，WAHAB M I M，FANG L.A comparison of regime-switchingtemperature modeling approaches for applications in weather derivativesJ.European Journal of Operational Research，2014，232（3）：549-560.12EVAREST E，BERNT
19、SSON F，SINGULL M，et al.Regime switching modelson temperature dynamicsJ.International Journal of Applied Mathematics andStatistics，2017，56（2）：19-36.13GYAMERAH S A，NGARE P，IKPE D.Regime-switching temperature dynamicsmodel for weather derivativesJ.International Journal of StochasticAnalysis，2018，32（1）：
20、1-15.14XIONG H，MAMON R.Putting a price tag on temperatureJ.Computational Management Science，2018，15（2）：259-296.15XIONG H，MAMON R.A self-updating model driven by a higher-orderhidden Markov chain for temperature dynamicsJ.Journal of ComputationalScience，2016，17（1）：47-61.16ALEXANDRIDIS A K，ZAPRANIS A
21、D.Weather derivatives：modeling and pricing weather-related riskM.New York：Springer Science&Business Media，2012.17HAMILTON J D.Analysis of time series subject to changes in regimeJ.Journal of Econometrics，1990，45（1/2）：39-70.18XIONG H，MAMON R.A higher-order Markov chain-modulated model forelectricity spot-price dynamicsJ.Applied Energy，2019，233（1）：495-515.