关于能源需求量预测常用模型及方法

1 关于能源需求量预测常用模型及方法

预测不仅要对事物未来的走势和可能出现的情景作定性分析，还要对某些关键的构成作定量判断，由于人们认识事物的能力有限，由于许多事物在已经经历过的过程中尚未能展现出许多重要因素的影响，即使对认识到的影响因素也难以有很深入的了解，因而人们的预测很难与未来的发展完全吻合，能够从大趋势上做到分析基本正确，就是成绩不小。一般说来，定性的预测，特别是定量的预测所涉及的对象越简单的、涉及未来的时间越短，预测的准确程度可能高些；涉及的事物越复杂、影响因素越多，涉及未来的时间越长，预测的准确性越差，这使许多定量的预测也仅带有半定量性质。在这种情况下，我们应把它看成是一种很可能出现的方案。针对这种情况，从预测上说，就要求：（1）把复杂事物的预测简化成相对单一的模式，或从多个相对简单的角度做出预测并加以综合；（2）随着发展，不断做出新的预测，这种新的预测不但可以收纳许多原来未料及或未重视的新的影响因素，而且可以用这一段时间已经发生的事实去修正原来的预测公式或参数。预测的基本思路是由过去和现在推断未来，由已知推测未知。因此，可将预测方法大致分为两类：由事物本身的发展过程推测未来和由已知的或设定的条件推导未来的油气发展状况：具体方法也可归纳为直接预测和间接预测。这里仅列举直接预测法。

该类方法从事物本身的历史和现状来认识影响它们的各种因素，利用这些因素的可能变化，推测未来一段时间内的变化趋势乃至变化量。换言之，对油气未来的预测主要依靠油气本身的变化。从逻辑上，把影响油气的各种因素起作用的综合结果——储量和产量（在一定程度上也表现为品质和空间分布）作为最直接的预测资料，将其作为最基本的依据来得出预测结论。这种预测，可以是简单的外预期，也可以考虑到多变量变化趋势的适应预期和合理预期。（1）客观存在受过去和现状影响太深，如果未来的基本框架与它们一致，这个预测的可靠性很大。而如果未来的影响因素发生了相当大的变化，这个预测结果可能产生相当大的偏差。而如果未来与历史有某种质的差异，预测结果有可能谬误。油气预测的对象是有生命周期的，处于不同生命阶段间的发展趋势可以是不同的。不注意这一点，简单的依据历史和现状，预测可能有趋势性错误。（2）在强调油气本身的因素时，往往对外部因素的影响考虑不周。如在对排放物有法律强制性要求的情况下，企业宁肯冒增加成本的风险，也要用清洁能源，这时对天然气的需求量就会加大。对此，单纯由过去和现在的增长趋势进行外推就难以预料。

对油气储、产量预测常用的方法归类的简单介绍，以利于后面选取几种预测方法，并在讨论后给出预测的推荐值。

（一）趋势分析法

根据历史上，特别是近期的发展趋势推测未来的一种方法，主要采用的是时间序列分析方法。时间序列分析方法是将历史资料按时间排列并对它进行分析，找出它随时间变化的趋势与规律，并用数学关系式把它表示出来。然后根据这个关系式对今后进行预测。该方法是进行短期预测（例如短期油气的生产量和消费量预测）的一种常用方法。其主要优点是灵活 2 简便，短期预测较准确。主要缺点是没有揭示影响事物变化的原因，仅从时序数据中寻找变化规律，因而不适用于长期的、特别是其发展阶段不同变化过程的预测。此外，对处于大起大伏的发育不成熟的事物，则预测效果也差。

趋势分析方法分确定性方法和随机性方法。

1、确定性趋势分析模型

（1）滑动平均预测模型

假设某预测量时间序列数据为1x，2x，…，tx，那么第t+1时期的预测值为1tx=（tx＋1tx＋K＋1Ntx）／N，即以N个实际值的平均值作为预测值，N的大小依时间序列数据显示的规律而选定。若以加权平均数为预测值，x＝（0atx＋1a1tx＋K＋1Na1Ntx）／N，称为加权滑动平均预测模型，。为加权因子，也依时间序列数据显示的规律而选定，满足（0a＋1a＋K＋1Na）／N＝1。若对实际时间序列数据进行一次滑动平均，得到滑动平均数序列，再对该平均数序列运用滑动平均预测模型，求得第t+1时期的预测值，称为二次滑动平均预测模型。

（2）指数平滑预测模型

以t时期实际值tx和预测值tx的加权平均值作为第t+1时期的预测值，1tx＝ttxaax)1(，0

（3）分解预测模型

将时间序列数据进行分解，求出趋势因素（ty）、季节因素（ts）、循环因素（tc）和随机因素（t），分别预测各种因素的预测值1ty、1ts和1t，利用模型1tx＝1ty＊1ts＊1t，求得第t+1时的预测值。在分解模型中，滑动平均是常用的分解长期趋势因素的方法。分解预测模型适用于季度能源消费量的预测。

2、随机性趋势分析模型

（1）自回归模型

tptptttaKaa2211，称为p阶自回归模型，简记为AR（p）。1a、2a、K、pa为自回归参数，t为随机误差项。利用回归分析技术估计参数，得到参数估计 3 量1a、2a、K、pa，于是第t+1时期的预测值为112111ptptttaKaa。

（2）自回归滑动平均模型

qtqtttptptttKaKaa22112211称为（p，q）阶自回归滑动平均模型，简记为ARMA（p，q）。1a、2a、K、pa、1、2、K、q，分别为自回归和滑动平均参数，由回归分析技术估计得到其估计值。

（二）生命模型拟合法

建立一个模型去逼近实际数据序列的过程称为拟合。判断拟合模型与实际体系的符合程度有多种标准，例如最小二乘法、最大离差为最小的最大最小原则、最大似然性、最小绝对偏差等，其中根据最小二乘法原理用拟合模型方程去逼近实际观察值是一种最常见的拟合方法。

根据拟合模型函数形式的不同，拟合形式有多种多样，诸如线性拟合、多项式拟合、对数拟合、指数拟合及生命模型拟合等。其中生命模型拟合是油气储量、产量预测的重要方法。生命模型函数是一类收敛型数学模型，常用的有翁氏生命旋回模型、逻辑斯谛Logistic）模型、龚帕兹（Compertz）模型等，它适用于石油、天然气、煤炭等勘探开发过程中类似生命总量有限体系的发展、消亡过程的描述和预测。它通过对过去油气储量、产量的变化趋势进行拟合，获取有关模型参数，再用相应的模型公式预测未来储量、产量的发展趋势。

1、翁氏生命旋回模型

翁文波教授创建的生命旋回理论认为，一切事物的发展变化可以划分为兴起、成长、达到鼎盛，然后逐渐衰亡4个阶段，从而形成一个完整的生命周期。油、气、煤等化石能源是一种不可再生资源，资源总量是有限的，就勘探开发的全过程来说，其发展规律总体上是符合生命旋回理论。

翁文波认为这种旋回的发展合乎泊松分布函数，就把它称为“泊松旋回”。“泊松旋回”问世后，立即得到石油界的广泛关注。1986年，石油工业部石油勘探开发科学研究院成立了项目组进行专题预测。从事数学地质研究的项目负责人高级工程师赵旭东利用半年时间，对国内外150多个油气田年产量和最终可采储量进行预测。他惊奇地发现，这些油气田的实际产量与预测值的相关系数，绝大多数都大于0.9，而正规开采的油气田的相关系数，一般都在0.95以上。这一系数的重要意义在于，能够根据油气田的以往产量非常准确地预测出最终可采储量。特别是可以发现不少油气田在原来计算时，没有包括进去的潜在可采储量，而且这部分可采储量往往是相当可观的。

赵旭东在实际预测中，发现“泊松旋回”的命名不当。不久，赵旭东写了一篇题目为“用Weng旋回模型对生命总量有限体系的预测”的文章，在中国科学院《科学通报》上发表。从此，国内外即正式称这种生命旋回为翁氏旋回或翁氏生命旋回。

该模型有3个假设： 4 （1）有一个起始点0x，在该点前系统x（t）不存在；

（2）过程开始后，x（t）的增长速率与系统x（t）本身成正比；

（3）增长率有一个比例因子使系统在某一时刻“b”达到有限极值时增长率为0。当tb时，增长率为负。这个因子可假设为（t/b—1）。

根据第（2）、（3）假设，可以建立下列微分方程：

)1)(()(tbtxdttdx

解上述方程并考虑第一个假设，则生命旋回模型可表达为：

 t0 x teat x t 0x t)(0t-b00不存在tx

式中：a，b为拟合系数；0x为预测初始时间；x为预测时间。根据最小二乘法原理，利用历史数据可以求出拟合系数a和b。

x（t）的总和是收敛的，从t=0至无限的积分为：

)1()(00badteatdttxstb

因为x（t）的总和是收敛的，这一模型只能用于总量有限的体系，如矿产资源的产量和储量等。

2、逻辑斯谛（Logistic）模型

逻辑斯谛（Logistic）预测模型是一个属于增长类型的模型。逻辑斯谛并非是人名，而只是Logistic英文一词的中文译名，它表示增长曲线的数学逻辑推理之义。美国哈伯特（Hubbert）于1962年首次提出逻辑斯谛曲线（10gistic curve）预测方法，因此，该法又被称为哈伯特模型。该模型已被用在勘探开发阶段的新老油气田可采储量的预测。

逻辑斯谛模型的微分方程为：

)1(byaydtdy （1）

式中：y为模型函数；t为时间变量；a、b为模型参数。

对（1）式积分后可得：

y=atceb1 （2）

（2）式就是逻辑斯谛模型的表达式。由（2）式可以看出，当a<0，t→∞时，y趋向于它的极值b。

由于油气田的累积产量随时间的变化属于增长曲线的类型。因此，为了预测油气田的可 5 采储量，可把（2）式改写为如下表达式：

atRpceNN1 （3）

式中：丹F为油气田的累积产量；NR为油气田的最终可采储量；t为开发年限。

将（3）式对时间，求导数可到预测油气产量（q）的关系式：

2)1(atatRceeacNq （4）

当（4）式的dtdq=0时，可求得油气田最高年产量发生时间（mt）的表达式：

atm1lnc （5）

把（5）式分别代人（4）、（3）式，得到油气田最高年产量（maxq）和最高年产量发生时的累积产量（pmN）的表达式：

RaNq41max

RpmNN21

逻辑斯谛模型预测油气田的产量与时间的变化关系时，它是一个带峰值的函数，而且最高年产量（峰值）刚好发生在累积采出可采储量50%的时间。因此，对于那些采出可采储量50%左右进入递减阶段的油气田，利用该模型可以得到相当满意的预测结果。否则，只能对油气田的可采储量进行有效预测。