一种复杂分散断块油田群投产顺序的决策方法

摘要

本发明属于断块油田开发的技术领域，具体的涉及一种复杂分散断块油田群投产顺序的决策方法。该决策方法通过在复杂分散断块油田群中选择典型断块油田建立油藏数值模拟模型，预测累计产油量。基于相关性分析方法优选投产顺序决策参数，结合插值方法，形成投产顺序决策参数对应累计产油量的数据样本。建立各时间节点下基于帝国竞争算法优化神经网络的产量预测模型，运用插值方法形成能够预测开发年限内每年累计产油量的预测模型。建立经济计算模型，计算各实际断块油田的NPV，实现复杂分散断块油田群投产顺序决策。

权利要求

1.一种复杂分散断块油田群投产顺序的决策方法，其特征在于，包括以下步骤：(1)首先基于复杂分散断块油田群中断块油田包括的构造特征、原油性质及天然能量大小在内的相关油田属性特征，选择典型断块油田建立典型油藏数值模拟模型；然后在复杂分散断块油田群的地质参数界限范围内，通过改变所建典型油藏数值模拟模型中的相关地质参数值，预测断块油田开发年限内的累计产油量；

(2)首先基于相关性分析方法，分析复杂分散断块油田群地质参数与开发年限内累计产油量的相关性，将相关性最大的地质参数作为复杂分散断块油田群投产顺序决策参数；

然后从开发年限中选择时间节点，建立复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下的累计产油量数据样本，通过插值方法计算，形成复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下的累计产油量关系曲线，扩充复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下累计产油量的数据样本；

(3)基于步骤(2)所扩充的复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下累计产油量数据样本，建立开发年限内各时间节点对应的基于帝国竞争算法优化BP神经网络的累计产油量预测模型，预测开发年限内每个时间节点下的累计产油量；所述步骤(3)中开发年限内各时间节点对应的基于帝国竞争算法优化BP神经网络的累计产油量预测模型的建立步骤如下：

①首先将复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下的累计产油量数据样本按比例划分为训练集与测试集，其中训练集所占比例至少70％；然后从中选择一个时间节点下的复杂分散断块油田群投产顺序决策参数与累计产油量数据样本；

②构建选择的时间节点下的BP神经网络模型，具体包括以下步骤：a.初步构建选择的时间节点下的BP神经网络模型结构：设计选择的时间节点下BP神经网络模型为由1个输入层，1个隐含层和1个输出层组成的3层回归网络结构；将选择的时间节点下训练集中的投产顺序决策参数作为BP神经网络输入层，选择的时间节点下训练集中对应的累计产油量作为输出层；隐含层节点数确定方法由公式(1)所示：式(1)中，hidden为隐含层节点数；n为输入层投产顺序决策参数个数；m为1～10之间的常数；

b.选择传递函数：隐含层和输出层的传递函数分别选择双极性S型函数和线性传递函数，由公式(2)～(3)所示：

foutput(x)＝x        (3)式(2)～(3)中，fhidden(x)为隐含层传递函数；foutput(x)为输出层传递函数，x为BP神经网络模型各层输出计算结果；

c.初始化选择的时间节点下BP神经网络模型：初始化权值和阈值，设定权值和阈值在的正态分布范围内随机取值，形成选择的时间节点下初始BP神经网络模型；

d.数据样本归一化：数据样本的归一化公式如公式(4)所示：式(4)中，l’为归一化计算后的数据；l为选择的时间节点下投产顺序决策参数与累计产油量形成的训练集；lmin和lmax为训练集的最小值和最大值；

e.确定选择的时间节点下BP神经网络模型训练过程评价指标：定义MS函数作为每个时间节点下BP神经网络模型训练过程的评价指标，计算公式如公式(5)所示：式(5)中，Ri为选择时间节点下训练集中累计产油量数据，Pi为模型预测的累计产油量数据，N为选择时间节点下训练集的样本容量；

③利用帝国竞争算法优化选择时间节点下BP神经网络模型，包括以下步骤：a.编码选择时间节点下BP神经网络模型的权值和阈值成为初始国家：利用帝国竞争算法优化选择时间节点下BP神经网络初始化的权值和阈值，结合上述②中b，将公式(5)改写成与权值和阈值对应的函数关系，如公式(6)所示；将MS函数作为帝国竞争算法的目标函数，将选择时间节点下BP神经网络中权值和阈值编码组合成一个实数码串，形成初始国家位置；

式(6)中，Rj为选择时间节点下训练集中累计产油量数据，n为输入层向量个数，μ为隐含层输出向量个数，m为输出层输出向量个数，ωei表示输入层投产顺序决策参数到隐含层连接权值，ωie表示隐含层到输出层累计产油量的连接权值，θe表示隐含层节点阈值，θj表示输出层节点阈值，fhidden(x)是隐含层传递函数，foutput(x)是输出层传递函数；

b.随机生成Nall个权值和阈值组成的实数码串，产生Nall个初始国家，利用③中a所述的目标函数MS计算国家的权力值，将Nall个初始国家实数码串按权力值大小排序，权力值大的前Np个实数码串作为初始强势国家，剩余的Nl个实数码串作为初始殖民地；

c.定义国家标准化权力值计算公式，如公式(7)所示，对初始强势国家标准化权力值的大小与初始殖民地数量Nl的乘积取整，完成初始殖民地对初始国家的分配，形成初始帝国；

式(7)中，Pn为第n个强势国家的标准化权力值，MSn为第n个强势国家的权力值；

d.通过帝国内作为殖民地的实数码串向作为强势国家的实数码串移动以更新选择时间节点下BP神经网络的权值和阈值，再返回选择时间节点下BP神经网络中计算目标函数MS，更新殖民地的权力值，若出现某个殖民地权力值大于该殖民地所属帝国中强势国家权力值，则将该殖民地与强势国家互换；

e.设定改革概率β控制改革的殖民地数量Nmn，如公式(8)～(9)所示；将满足②中c所述的随机权值和阈值构成的实数码串与改革殖民地的实数码串进行交换，再次更新改革殖民地的权值和阈值，改革殖民地的权力值即再次更新，若出现某个殖民地权力值大于该殖民地所属帝国中强势国家权力值，则将该殖民地与强势国家互换；

Ncn＝round{Pn·Nl}                (8)Nmn＝Ncn×β              (9)式(8)～(9)中，Ncn为n个强势国家所分配到的殖民地的数量；Nmn为某帝国中进行改革的殖民地数量；β为改革概率；

f.计算帝国总权力值，其计算公式为帝国中强势国家权力值和平均殖民地权力值之和，如公式(10)所示；根据各帝国总权力值的大小，将总权力值最小帝国的殖民地向其他帝国分配，分配概率计算公式如公式(11)所示；

式(10)中，Tn为第n个帝国总权力值，η为影响因子，0<η<1；

式(11)中，Ptn为第n个帝国分配最弱帝国殖民地的概率；

所有帝国的占领概率形成一个向量P，构造一个和向量P同维的随机向量R，R满足0～1上均匀分布，向量P与向量R相减得到向量D，计算公式如公式(12)～(14)所示；向量D中元素数值最大的帝国将分配最弱殖民地，使得选择时间节点下BP神经网络模型的权值和阈值接近最优解；

P＝[P1,P2,L,Pn]   (12)R＝[r1,r2,L,rn]ri～U(0,1)   (13)D＝P‑R   (14)

g.消灭失去殖民地的帝国，即抛弃使选择时间节点下BP神经网络模型预测误差较大的权值和阈值；

h.重复③中d～g，只剩一个帝国即只剩最优权值和阈值，或达到最大迭代次数，则终止算法，将最优权值和阈值代入至选择时间节点下BP神经网络模型中，构建选择时间节点下帝国竞争算法优化BP神经网络累计产油量预测模型；

2

④设定平均绝对百分比误差MAPE，相关系数R作为评价指标，用于评价选择时间节点下帝国竞争算法优化BP神经网络累计产油量预测模型的可靠性，计算公式如公式(15)～(16)所示；

y y

式(15)～(16)中，Ri为选择时间节点下测试集中累计产油量数据，Pi 为训练后模型预测测试 集累计产油量数据，N为样本数量；

⑤将选择时间节点下测试集中投产顺序决策参数作为输入层，通过选择时间节点下的帝国竞争算法优化BP神经网络累计产油量预测模型预测该时间节点下的累计产油量，计算2

MAPE和R ；优选出的帝国竞争算法优化BP神经网络累计产油量预测模型可靠的标准为MAPE2

小于10％，R大于0.9；

⑥选择其他时间节点下的复杂分散断块油田群投产顺序决策参数与累计产油量数据样本，重复①～⑤，建立各类典型油藏开发年限内每个时间节点下的帝国竞争算法优化BP神经网络累计产油量预测模型；

(4)运用插值方法，根据步骤(3)所得预测结果，计算不同典型油藏各时间节点与累计产油量预测结果形成的对应关系曲线，形成开发年限内每年累计产油量的预测模型并验证精度；

(5)考虑复杂分散断块油田群投资因素，计算投资成本；结合经济参数，建立NPV经济计算模型；基于开发年限内每年累计产油量计算NPV；

(6)利用所建开发年限内每年累计产油量的预测模型预测实际复杂分散断块油田群中各断块油田开发年限内每年的累计产油量，并依据所建NPV经济计算模型计算NPV；根据NPV大小，完成复杂分散断块油田群投产顺序的决策。

2.根据权利要求1所述复杂分散断块油田群投产顺序的决策方法，其特征在于，所述步骤(1)中复杂分散断块油田群中的断块油田数量大于30；所选择的典型断块油田数量为3～

6个；复杂分散断块油田群的地质参数包括含油面积、油层厚度、平均渗透率、平均孔隙度、平均含油饱和度、油藏原始压力、水体倍数、原油粘度、原油密度、饱和压力、地层温度和原始气油比；典型油藏数值模拟模型中相关地质参数值的改变次数至少200次。

3.根据权利要求1所述复杂分散断块油田群投产顺序的决策方法，其特征在于，所述步骤(2)中相关性分析方法为相关系数法、信息熵和灰色关联法；其中相关系数法中包括皮尔森系数、斯皮尔曼系数以及肯德尔相关性系数；插值方法为分段三次赫米特插值法；从开发年限中所选择的各个时间节点之间的间隔小于等于6年；最终复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下的累计产油量数据样本数量至少10000组。

4.根据权利要求1所述复杂分散断块油田群投产顺序的决策方法，其特征在于，所述步骤(4)中插值方法为分段三次赫米特插值法。

5.根据权利要求1所述复杂分散断块油田群投产顺序的决策方法，其特征在于，所述步骤(4)中验证精度的具体操作如下：采用测试集中不同典型断块油藏数值模拟计算结果作为验证样本，选择所述帝国竞争算法优化BP神经网络累计产油量预测模型建立的步骤④中2

平均绝对百分比误差MAPEH，相关系数RH作为评价指标，计算公式如公式(17)～(18)所示；

2

定义MAPEH小于10％，RH大于0.9即验证模型精度在有效范围内；

H H

式(17)～(18)中，Ri为插值计算形成的累计产油量数据，Pi 为验证样本中累计产油量数据，NH为验证样本容量。

6.根据权利要求1所述复杂分散断块油田群投产顺序的决策方法，其特征在于，所述步骤(5)中复杂分散断块油田群投资因素包括总井数、累计产油量、边际成本以及与原油中心处理站的地面距离权重；通过多元线性回归方法计算投资成本，计算公式如公式(19)所示：CAPEXCPF＝e+aωCPF+bn+cOilproduct+dωborder          (19)式(19)中，CAPEXCPF为地面总投资，ωCPF为断块油田与CPF的地面距离权重，n为总井数，Oilproduct为累计产油量，ωborder为边际成本即计算的断块与大断块的距离，e、a、b、c、d为多元回归系数。

7.根据权利要求1所述复杂分散断块油田群投产顺序的决策方法，其特征在于，所述步骤(5)中NPV计算公式如公式(20)～(24)所示：式(20)中，t为时间步长，单位为年；T表示总生产时间，单位为年；i为年折现率，单位为％；CPFt为时间t的运营现金流；CEX为投资金额；

运营现金流CPFt计算公式如下：

CPFt＝INCOME‑OPEX                    (21)式(21)中，INCOME为收入，OPEX为操作总投资；

其中收入INCOME计算公式如下：式(22)中，Rit为国内企业收入比例；λ为原油的商品率； 表示t年的累计产油量；Poil表示石油销售价格；

操作总投资OPEX计算公式如下：

式(23)中，Re为国内企业投资比例，n表示总井数，μs表示固定操作费，μv表示可变操作费，a表示操作费上涨率；

投资金额CEX计算公式为：

CEX＝Re[n(CEXO+CEXD)(1+θt)+CEXCPF]        (24)式(24)中，n表示总井数；CEXO表示单井采油投资；CEXD表示单井钻井投资；θt为弃置率，表示当年产油量达到了总体的75％以上时， θt为5％，反之θt为零；CEXCPF为地面总投资。

说明书

一种复杂分散断块油田群投产顺序的决策方法

技术领域

[0001] 本发明属于断块油田开发的技术领域，具体的涉及一种复杂分散断块油田群投产顺序的决策方法。

背景技术

[0002] 复杂断块油田分布广泛，在油气资源的勘探与开发中占据重要地位，通常在一定构造背景下若干个具有独立水动力系统、面积大小各异的含油断块形成一个复杂分散的断
块油田群。复杂分散断块油田群具有多含油层系、多油藏类型、多压力系统以及流体性质多
变等特性，同时由于断层发育，断块油田群各层油气水界面不统一，油、气、水三者关系复
杂。仅断块油田群自身复杂多样性的特点便为如何科学、客观地决策开发投产顺序带来了
诸多难点。此外，针对该类复杂分散断块油田群的开发，还需要兼顾地面上较小断块油田到
大断块油田以及原油中心处理站远近等因素所带来的投资差异，更加大了复杂分散断块油
田群的投产顺序决策难度。

[0003] 目前，对于复杂分散断块油田群投产顺序的决策方法较少，主要以层次分析法、综合模糊评判法等主观性较强的决策方法为主。该类方法通过选择复杂分散断块油田群中断
块油田的指标参数，进而计算各断块油田的综合评价系数，以此作为投产顺序决策依据，此
类方法难以体现复杂分散断块油田群的投资特点。若采用传统油藏数值模拟方法虽然能够
较为客观地预测每个断块油田的开发效果，但复杂分散断块油田群中断块油田数量众多，
按照传统油藏数值模拟方法需要对每个断块油田开展数值模拟研究工作，工作量巨大，实
际操作难度高，难以实现。

[0004] 因此亟需提供一种基于复杂分散断块油田群特点的投产顺序决策方法，既能够尽可能全面地考虑诸多因素，满足科学、客观性；又具有高效的实际可操作性。

发明内容

[0005] 本发明的目的在于针对目前复杂分散断块油田群投产顺序无法客观高效决策的现状而提供一种复杂分散断块油田群投产顺序的决策方法，该方法能够高效准确地完成对
复杂分散断块油田群投产顺序的决策，对复杂分散断块油田群动用开发具有重要意义。

[0006] 本发明的技术方案为：一种复杂分散断块油田群投产顺序的决策方法，包括以下步骤：

[0007] (1)首先基于复杂分散断块油田群中断块油田的包括构造特征、原油性质及天然能量大小在内的相关油田属性特征，选择典型断块油田建立典型油藏数值模拟模型；然后
在复杂分散断块油田群的地质参数界限范围内，通过改变所建典型油藏数值模拟模型中的
相关地质参数值，预测断块油田开发年限内的累计产油量；

[0008] (2)首先基于相关性分析方法，分析复杂分散断块油田群地质参数与开发年限内累计产油量的相关性，将相关性最大的地质参数作为复杂分散断块油田群投产顺序决策参
数；然后从开发年限中选择时间节点，建立复杂分散断块油田群投产顺序决策参数与开发
年限内每个时间节点下的累计产油量数据样本，通过插值方法计算，形成复杂分散断块油
田群投产顺序决策参数与开发年限内每个时间节点下的累计产油量关系曲线，扩充复杂分
散断块油田群投产顺序决策参数与开发年限内每个时间节点下累计产油量的数据样本；

[0009] (3)基于步骤(2)所扩充的复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下累计产油量数据样本，建立开发年限内各时间节点对应的基于帝国竞争算
法优化BP神经网络的累计产油量预测模型，预测开发年限内每个时间节点下的累计产油
量；

[0010] (4)运用插值方法，根据步骤(3)所得预测结果，计算不同典型油藏各时间节点与累计产油量预测结果形成的对应关系曲线，形成开发年限内每年累计产油量的预测模型并
验证精度；

[0011] (5)考虑复杂分散断块油田群投资因素，计算投资成本；结合经济参数，建立NPV经济计算模型；基于开发年限内每年累计产油量计算NPV；

[0012] (6)利用所建开发年限内每年累计产油量的预测模型预测实际复杂分散断块油田群中各断块油田开发年限内每年的累计产油量，并依据所建NPV经济计算模型计算NPV；根
据NPV大小，完成复杂分散断块油田群投产顺序的决策。

[0013] 所述步骤(1)中复杂分散断块油田群中的断块油田数量大于30；所选择的典型断块油田数量为3～6个；复杂分散断块油田群的地质参数包括含油面积、油层厚度、平均渗透
率、平均孔隙度、平均含油饱和度、油藏原始压力、水体倍数、原油粘度、原油密度、饱和压
力、地层温度和原始气油比；典型油藏数值模拟模型中相关地质参数值的改变次数至少200
次。

[0014] 所述步骤(2)中相关性分析方法为相关系数法、信息熵和灰色关联法；其中相关系数法中包括皮尔森系数、斯皮尔曼系数以及肯德尔相关性系数；插值方法为分段三次赫米
特插值法；从开发年限中所选择的各个时间节点之间的间隔小于等于6年；最终复杂分散断
块油田群投产顺序决策参数与开发年限内每个时间节点下的累计产油量数据样本数量至
少10000组。

[0015] 所述步骤(3)中开发年限内各时间节点对应的基于帝国竞争算法优化BP神经网络的累计产油量预测模型的建立步骤如下：

[0016] ①首先将复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下的累计产油量数据样本按比例划分为训练集与测试集，其中训练集所占比例至少70％；然
后从中选择一个时间节点下的复杂分散断块油田群投产顺序决策参数与累计产油量数据
样本；

[0017] ②构建选择的时间节点下的BP神经网络模型，具体包括以下步骤：

[0018] a.初步构建选择的时间节点下的BP神经网络模型结构：设计选择的时间节点下BP神经网络模型为由1个输入层，1个隐含层和1个输出层组成的3层回归网络结构；将选择的
时间节点下训练集中的投产顺序决策参数作为BP神经网络输入层，选择的时间节点下训练
集中对应的累计产油量作为输出层；隐含层节点数确定方法由公式(1)所示：

[0019]

[0020] 式(1)中，hidden为隐含层节点数；n为输入层投产顺序决策参数个数；m为1～10之间的常数；

[0021] b.选择传递函数：隐含层和输出层的传递函数分别选择双极性S型函数和线性传递函数，由公式(2)～(3)所示：

[0022]

[0023] foutput(x)＝x        (3)

[0024] 式(2)～(3)中，fhidden(x)为隐含层传递函数；foutput(x)为输出层传递函数，x为BP神经网络模型各层输出计算结果；

[0025] c.初始化选择的时间节点下BP神经网络模型：初始化权值和阈值，设定权值和阈值在 的正态分布范围内随机取值，形成选择的时间节点下初始BP神经网络模型；

[0026] d.数据样本归一化：数据样本的归一化公式如公式(4)所示：

[0027]

[0028] 式(4)中，l’为归一化计算后的数据；l为选择的时间节点下投产顺序决策参数与累计产油量形成的训练集；lmin和lmax为训练集的最小值和最大值；

[0029] e.确定选择的时间节点下BP神经网络模型训练过程评价指标：定义MS函数作为每个时间节点下BP神经网络模型训练过程的评价指标，计算公式如公式(5)所示。

[0030]

[0031] 式(5)中，Ri为选择时间节点下训练集中累计产油量数据，Pi为模型预测的累计产油量数据，N为选择时间节点下训练集的样本容量；

[0032] ③利用帝国竞争算法优化选择时间节点下BP神经网络模型，包括以下步骤：

[0033] a.编码选择时间节点下BP神经网络模型权值和阈值成为初始国家：利用帝国竞争算法优化选择时间节点下BP神经网络初始化的权值和阈值，结合上述②中b，将公式(5)改
写成与权值和阈值对应的函数关系，如公式(6)所示；将MS函数作为帝国竞争算法的目标函
数，将选择时间节点下BP神经网络中权值和阈值编码组合成一个实数码串，形成初始国家
位置；

[0034]

[0035] 式(6)中，Rj为选择时间节点下训练集中累计产油量数据，n为输入层向量个数，μ为隐含层输出向量个数，m为输出层输出向量个数，ωei表示输入层投产顺序决策参数到隐
含层连接权值，ωie表示隐含层到输出层累计产油量的连接权值，θe表示隐含层节点阈值，
θj表示输出层节点阈值，fhidden(x)是隐含层传递函数，foutput(x)是输出层传递函数；

[0036] b.随机生成Nall个权值和阈值组成的实数码串，产生Nall个初始国家，利用③中a所述的目标函数MS计算国家的权力值，将Nall个初始国家实数码串按权力值大小排序，权力值
大的前Np个实数码串作为初始强势国家，剩余的Nl个实数码串作为初始殖民地；

[0037] c.定义国家标准化权力值计算公式，如公式(7)所示，对初始强势国家标准化权力值的大小与初始殖民地数量Nl的乘积取整，完成初始殖民地对初始国家的分配，形成初始
帝国；

[0038]

[0039] 式(7)中，pn为第n个强势国家的标准化权力值，MSn为第n个强势国家的权力值；d.通过帝国内作为殖民地的实数码串向作为强势国家的实数码串移动以更新选择时间节点
下BP神经网络的权值和阈值，再返回选择时间节点下BP神经网络中计算目标函数MS，更新
殖民地的权力值，若出现某个殖民地权力值大于该殖民地所属帝国中强势国家权力值，则
将该殖民地与强势国家互换；

[0040] e.设定改革概率β控制改革的殖民地数量Nmn，如公式(8)～(9)所示；将满足②中c所述的随机权值和阈值构成的实数码串与改革殖民地的实数码串进行交换，再次更新改革
殖民地的权值和阈值，改革殖民地的权力值即再次更新，若出现某个殖民地权力值大于该
殖民地所属帝国中强势国家权力值，则将该殖民地与强势国家互换；

[0041] Ncn＝round{Pn·Nl}           (8)

[0042] Nmn＝Ncn×β             (9)

[0043] 式(8)～(9)中，Ncn为n个强势国家所分配到的殖民地的数量；Nmn为某帝国中进行改革的殖民地数量；β为改革概率；

[0044] f.计算帝国总权力值，其计算公式为帝国中强势国家权力值和平均殖民地权力值之和，如公式(10)所示；根据各帝国总权力值的大小，将总权力值最小帝国的殖民地向其他
帝国分配，分配概率计算公式如公式(11)所示；

[0045]

[0046] 式(10)中，Tn为第n个帝国总权力值，η为影响因子，0<η<1；

[0047]

[0048] 式(11)中，Ptn为第n个帝国分配最弱帝国殖民地的概率；

[0049] 所有帝国的占领概率形成一个向量P，构造一个和向量P同维的随机向量R，R满足0～1上均匀分布，向量P与向量R相减得到向量D，计算公式如公式(12)～(14)所示；向量D中
元素数值最大的帝国将分配最弱殖民地，使得选择时间节点下BP神经网络模型的权值和阈
值接近最优解；

[0050] P＝[P1,P2,L,Pn]  (12)

[0051] R＝[r1,r2,L,rn] ri～U(0,1)  (13)

[0052] D＝P‑R  (14)

[0053] g.消灭失去殖民地的帝国，即抛弃使选择时间节点下BP神经网络模型预测误差较大的权值和阈值；

[0054] h.重复③中d～g，只剩一个帝国即只剩最优权值和阈值，或达到最大迭代次数，则终止算法，将最优权值和阈值代入至选择时间节点下BP神经网络模型中，构建选择时间节
点下帝国竞争算法优化BP神经网络累计产油量预测模型；

[0055] ④设定平均绝对百分比误差MAPE，相关系数R2作为评价指标，用于评价选择时间节点下帝国竞争算法优化BP神经网络累计产油量预测模型的可靠性，计算公式如公式(15)
～(16)所示；

[0056]

[0057]

[0058] 式(15)～(16)中， 为选择时间节点下测试集中累计产油量数据， 为训练后模型预测验证集累计产油量数据，N为样本数量；

[0059] ⑤将选择时间节点下测试集中投产顺序决策参数作为输入层，通过选择时间节点下的帝国竞争算法优化BP神经网络累计产油量预测模型预测该时间节点下的累计产油量，
2
计算MAPE和R ；优选出的帝国竞争算法优化BP神经网络累计产油量预测模型可靠的标准为
2
MAPE小于10％，R大于0.9；

[0060] ⑥选择其他时间节点下的复杂分散断块油田群投产顺序决策参数与累计产油量数据样本，重复①～⑤，建立各类典型油藏开发年限内每个时间节点下的帝国竞争算法优
化BP神经网络累计产油量预测模型。

[0061] 所述步骤(4)中插值方法为分段三次赫米特插值法。

[0062] 所述步骤(4)中验证精度的具体操作如下：采用测试集中不同典型断块油藏数值模拟计算结果作为验证样本，选择所述帝国竞争算法优化BP神经网络累计产油量预测模型
2
建立的步骤④中平均绝对百分比误差MAPEH，相关系数RH 作为评价指标，计算公式如公式
2
(17)～(18)所示；定义MAPEH小于10％，RH大于0.9即验证模型精度在有效范围内；

[0063]

[0064]

[0065] 式(17)～(18)中， 为插值计算形成的累计产油量数据， 为验证样本中累计产油量数据，NH为验证样本容量。

[0066] 所述步骤(5)中复杂分散断块油田群投资因素包括总井数、累产油量、边际成本以及与原油中心处理站的地面距离权重；通过多元线性回归方法计算投资成本，计算公式如
公式(19)所示：

[0067] CAPEXCPF＝e+aωCPF+bn+cOilproduct+dωborder  (19)

[0068] 式(19)中，CAPEXCPF为地面总投资，ωCPF为断块油田与CPF的地面距离权重，n为井数，Oilproduct为累计产油量，ωborder为边际成本即计算的断块与大断块的距离，e、a、b、c、d
为多元回归系数。

[0069] 所述步骤(5)中NPV计算公式如公式(20)～(24)所示：

[0070]

[0071] 式(20)中，t为时间步长，单位为年；T表示总生产时间，单位为年；i为年折现率，单位为％；CPFt为时间t的运营现金流；CEX为投资金额；

[0072] 运营现金流CPFt计算公式如下：

[0073] CPFt＝INCOME‑OPEX               (21)

[0074] 式(21)中，INCOME为收入，OPEX为操作总投资；

[0075] 其中收入INCOME计算公式如下：

[0076]

[0077] 式(22)中，Rit为国内企业收入比例；λ为原油的商品率； 表示t年的累计产油量；Poil表示石油销售价格；

[0078] 操作总投资OPEX计算公式如下：

[0079]

[0080] 式(23)中，Re为国内企业投资比例，n表示总井数，μs表示固定操作费，μv表示可变操作费表示可变操作费，a表示操作费上涨率；

[0081] 投资金额CEX计算公式为：

[0082] CEX＝Re[n(CEXO+CEXD)(1+θt)+CEXCPF]   (24)

[0083] 式(24)中，n表示总井数；CEXO表示单井采油投资；CEXD表示单井钻井投资；θt为弃置率，表示当该年产油量达到了总体的75％以上θt为5％，反之θt为零；CEXCPF为地面总投资。

[0084] 本发明的有益效果为：本发明通过在复杂分散断块油田群中选择典型断块油田建立油藏数值模拟模型，预测累计产油量。基于相关性分析方法优选投产顺序决策参数，结合
插值方法，形成投产顺序决策参数对应累计产油量的数据样本。建立各时间节点下基于帝
国竞争算法优化神经网络的产量预测模型，运用插值方法形成能够预测开发年限内每年累
计产油量的预测模型。建立经济计算模型，计算各实际断块油田的NPV，实现复杂分散断块
油田群投产顺序决策。相较于现有技术，具有如下优势：

[0085] 1、本发明突破了现有断块油田群投产顺序决策方法较强主观性的局限，通过油藏数值模拟方法，建立复杂分散断块油田群投产顺序决策数据样本，为本方法提供了客观基
础。通过帝国竞争算法优化BP神经网络模型的方法代替了对每个断块油田的数值模拟研
究，大大节省了工作量与计算成本，能够高效准确地完成对复杂分散断块油田群投产顺序
的决策。

[0086] 2、本发明通过在开发年限内设置时间节点，结合插值方法，有效扩充了复杂分散断块油田群投产顺序决策参数与累计产油量数据样本，能够较好地提高累计产油量预测模
型训练精度并防止过拟合等情况，为后续理论研究提供了充分的数据支持。

[0087] 3、本发明基于帝国竞争算法优化BP神经网络模型的权值和阈值，改善了BP神经网络模型自身训练精度低、易局部最小化等缺陷，提高了BP神经网络模型预测累计产油量的
精度和效率，增加了复杂分散断块油田群投产顺序决策方法的准确性和可靠性。

[0088] 4、本发明针对复杂分散断块油田群地面投资差异性，提出了复杂分散断块油田群地面投资计算方法，有效结合了复杂分散断块油田群开发特色，形成了NPV计算方法，为复
杂分散断块油田群投产顺序决策提供了可靠的参考依据。

附图说明

[0089] 图1为本发明的方法流程图。

[0090] 图2为实施例1中的典型边水单斜稀油油藏数值模型示例图。

[0091] 图3为实施例1中的平均渗透率插值扩充效果示例图。

[0092] 图4为实施例1中地面投资模型回归R2图。

[0093] 图5为实施例1中预测实际断块油田累计产油量效果图。

[0094] 图6为实施例1中对某复杂分散断块油田群102个断块油田投产顺序排序结果。

具体实施方式

[0095] 下面结合附图，对本发明的技术方案进行详细的说明。

[0096] 实施例1

[0097] 为使本发明的目的、技术方案和优点更加清楚，下面对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基
于本发明中的实施例，本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实
施例，都属于本发明保护的范围。

[0098] 以某断块油田群为例对所提出的技术进行验证，该断块油田群中的断块油田数量为102个。

[0099] 所述复杂分散断块油田群投产顺序的决策方法，包括以下步骤：

[0100] (1)首先基于复杂分散断块油田群中断块油田的包括构造特征、原油性质及天然能量大小在内的相关油田属性特征，选择4个典型断块油田建立典型油藏数值模拟模型，分
别为边水单斜典型稀油油藏、含气顶断垒典型稀油油藏、岩性典型稀油油藏、边水断垒典型
稠油油藏；然后在复杂分散断块油田群的包括含油面积、油层厚度、平均渗透率、平均孔隙
度、平均含油饱和度、油藏原始压力、水体倍数、原油粘度、原油密度、饱和压力、地层温度和
原始气油比的地质参数界限范围内，通过改变所建典型油藏数值模拟模型中的相关地质参
数值300次，以预测断块油田开发年限35年内的累计产油量；

[0101] (2)首先基于斯皮尔曼系数法的相关性分析方法，分析复杂分散断块油田群地质参数与开发年限内累计产油量的相关性，筛选出复杂分散断块油田群投产顺序决策参数为
平均孔隙度、平均渗透率、平均含油饱和度、油层厚度、含油面积、原油粘度、水体倍数；然后
以5年的间隔选择时间节点分别为开发5年、开发10年、开发15年、开发20年、开发25年、开发
30年、开发35年共7个时间节点，建立复杂分散断块油田群投产顺序决策参数与开发年限内
每个时间节点下的累计产油量数据样本，通过分段三次赫米特插值法计算，形成复杂分散
断块油田群投产顺序决策参数与开发年限内每个时间节点下的累计产油量关系曲线，扩充
复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下累计产油量的数据
样本，数量为14000组；

[0102] (3)基于步骤(2)所扩充的复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下累计产油量数据样本，建立开发年限内各时间节点对应的基于帝国竞争算
法优化BP神经网络的累计产油量预测模型，预测开发年限内每个时间节点下的累计产油
量；其中累计产油量预测模型的建立步骤如下：

[0103] ①首先将复杂分散断块油田群投产顺序决策参数与开发年限内每个时间节点下的累计产油量数据样本按比例划分为训练集与测试集，其中训练集所占比例为70％，测试
集所占比例为30％；然后从中选择一个时间节点下的复杂分散断块油田群投产顺序决策参
数与累计产油量数据样本；

[0104] ②构建选择的时间节点下的BP神经网络模型，具体包括以下环节：

[0105] a.初步构建选择的时间节点下的BP神经网络模型结构：设计选择的时间节点下BP神经网络模型为由1个输入层，1个隐含层和1个输出层组成的3层回归网络结构；将选择的
时间节点下训练集中的投产顺序决策参数作为BP神经网络输入层，选择的时间节点下训练
集中对应的累计产油量作为输出层；隐含层节点数确定方法由公式(1)所示：

[0106]

[0107] 式(1)中，hidden为隐层节点数；n为输入层投产顺序决策参数个数；m为1～10之间的常数；通过计算确定隐含层节点数为3个。

[0108] b.选择传递函数：隐含层和输出层的传递函数分别选择双极性S型函数和线性传递函数，由公式(2)～(3)所示：

[0109]

[0110] foutput(x)＝x        (3)

[0111] 式(2)～(3)中，fhidden(x)为隐含层传递函数；foutput(x)为输出层传递函数，x为BP神经网络模型各层输出计算结果；

[0112] c.初始化选择的时间节点下BP神经网络模型：初始化权值和阈值，设定权值和阈值在 的正态分布范围内随机取值，形成选择的时间节点下初始BP神经网络模型；

[0113] d.数据样本归一化：数据样本的归一化公式如公式(4)所示：

[0114]

[0115] 式(4)中，l’为归一化计算后的数据；l为选择的时间节点下投产顺序决策参数与累计产油量形成的训练集；lmin和lmax为训练集的最小值和最大值；

[0116] e.确定选择的时间节点下BP神经网络模型训练过程评价指标：定义MS函数作为每个时间节点下BP神经网络模型训练过程的评价指标，计算公式如公式(5)所示。

[0117]

[0118] 式(5)中，Ri为选择时间节点下训练集中累计产油量数据，Pi为模型预测的累计产油量数据，N为选择时间节点下训练集的样本容量；

[0119] ③利用帝国竞争算法优化选择时间节点下BP神经网络模型，包括以下环节：

[0120] a.编码选择时间节点下BP神经网络模型权值和阈值成为初始国家：帝国竞争算法优化选择时间节点下BP神经网络初始化的权值和阈值，结合上述②中b，将公式(5)改写成
与权值和阈值对应的函数关系，如公式(6)所示；将MS函数作为帝国竞争算法的目标函数，
将选择时间节点下BP神经网络中权值和阈值编码组合成一个实数码串，形成初始国家位
置；

[0121]

[0122] 式(6)中，Rj为选择时间节点下训练集中累计产油量数据，n为输入层向量个数，μ为隐含层输出向量个数，m为输出层输出向量个数，ωei表示输入层投产顺序决策参数到隐
含层连接权值，ωie表示隐含层到输出层累计产油量的连接权值，θe表示隐含层节点阈值，
θj表示输出层节点阈值，fhidden(x)是隐含层传递函数，foutput(x)是输出层传递函数；

[0123] b.随机生成40个权值和阈值组成的实数码串，产生40个初始国家，③中a所述的目标函数MS计算国家的权力值，将40个初始国家实数码串按权力值大小排序，选择初始优势
国家数量为10个，初始殖民地数量为30个；

[0124] c.定义国家标准化权力值计算公式，如公式(7)所示，对初始强势国家标准化权力值的大小与初始殖民地数量Nl的乘积取整，完成初始殖民地对初始国家的分配，形成初始
帝国；

[0125]

[0126] 式(7)中，pn为第n个强势国家的标准化权力值，MSn为第n个强势国家的权力值；d.通过帝国内作为殖民地的实数码串向作为强势国家的实数码串移动以更新选择时间节点
下BP神经网络的权值和阈值，再返回选择时间节点下BP神经网络中计算目标函数MS，更新
殖民地的权力值，若出现某个殖民地权力值大于该殖民地所属帝国中强势国家权力值，则
将该殖民地与强势国家互换；

[0127] e.设定改革概率β控制改革的殖民地数量Nmn，如公式(8)～(9)所示，其中β为0.2；将满足②中c所述的随机权值和阈值构成的实数码串与改革殖民地的实数码串进行交换，
再次更新改革殖民地的权值和阈值，改革殖民地的权力值即再次更新，若出现某个殖民地
权力值大于该殖民地所属帝国中强势国家权力值，则将该殖民地与强势国家互换；

[0128] Ncn＝round{Pn·Nl}            (8)

[0129] Nmn＝Ncn×β              (9)

[0130] 式(8)～(9)中，Ncn为n个优势国家所分配到的殖民地的数量；Nmn为某帝国中进行改革的殖民地数量；β为改革概率；

[0131] f.计算帝国总权力值，其计算公式为帝国中强势国家权力值和平均殖民地权力值之和，如公式(10)所示；根据各帝国总权力值的大小，将总权力值最小帝国的殖民地向其他
帝国分配，分配概率计算公式如公式(11)所示；

[0132]

[0133] 式(10)中，Tn为第n个帝国总权力值，η为影响因子，0<η<1；

[0134]

[0135] 式(11)中，Ptn为第n个帝国分配最弱帝国殖民地的概率；

[0136] 所有帝国的占领概率形成一个向量P，构造一个和向量P同维的随机向量R，R满足0～1上均匀分布，向量P与向量R相减得到向量D，计算公式如公式(12)～(14)所示；向量D中
元素数值最大的帝国将分配最弱殖民地，使得选择时间节点下BP神经网络模型的权值和阈
值接近最优解；

[0137] P＝[P1,P2,L,Pn]  (12)

[0138] R＝[r1,r2,L,rn] ri～U(0,1)  (13)

[0139] D＝P‑R  (14)

[0140] g.消灭失去殖民地的帝国，即抛弃使选择时间节点下BP神经网络模型预测误差较大的权值和阈值；

[0141] h.重复③中d～g，至达到最大迭代次数200次，则终止算法，将最优权值和阈值代入至选择时间节点下BP神经网络模型中，构建选择时间节点下帝国竞争算法优化BP神经网
络累计产油量预测模型；

[0142] ④设定平均绝对百分比误差MAPE，相关系数R2作为评价指标，用于评价选择时间节点下帝国竞争算法优化BP神经网络累计产油量预测模型的可靠性，计算公式如公式(15)
～(16)所示；

[0143]

[0144]

[0145] 式(15)～(16)中， 为选择时间节点下测试集中累计产油量数据，Pi为训练后模型预测验证集累计产油量数据，N为样本数量；

[0146] ⑤将选择时间节点下测试集中投产顺序决策参数作为输入层，通过选择时间节点下的帝国竞争算法优化BP神经网络累计产油量预测模型预测该时间节点下的累计产油量，
2
计算MAPE和R ；优选出的帝国竞争算法优化BP神经网络累计产油量预测模型可靠的标准为
2
MAPE小于10％，R大于0.9；

[0147] ⑥选择其他时间节点下的复杂分散断块油田群投产顺序决策参数与累计产油量数据样本，重复①～⑤，建立各类典型油藏开发年限内每个时间节点下的帝国竞争算法优
化BP神经网络累计产油量预测模型。四类典型油藏对应各时间节点帝国竞争算法优化BP神
经网络累计产油量预测模型共28组，所有模型的MAPE范围为2.19％～8.18％，R2范围为
0.9872～0.9963，符合优选要求。

[0148] (4)运用分段三次赫米特插值法，根据步骤(3)所得预测结果，计算不同典型油藏各时间节点与累计产油量预测结果形成的对应关系曲线，形成开发年限内每年累计产油量
的预测模型并验证精度；其中验证精度的具体操作如下：随机从4类典型断块油藏数值模拟
计算结果中每类选取50组数据作为验证样本，共200组数据，选择所述帝国竞争算法优化BP
2
神经网络累计产油量预测模型建立的步骤④中平均绝对百分比误差MAPEH，相关系数RH 作
2
为评价指标，计算公式如公式(17)～(18)所示；定义MAPEH小于10％，RH大于0.9即验证模型
精度在有效范围内；

[0149]

[0150]

[0151] 式(17)～(18)中， 为插值计算形成的累计产油量数据， 为验证样本中累计产2
油量数据，NH为验证样本容量；累计产油量预测模型的MAPE范围为5.16％～9.82％，R范围
为0.9172～0.9863，符合精度要求，验证了4类典型断块油田每年累计产油量的预测模型精
度。

[0152] (5)考虑复杂分散断块油田群投资因素，包括总井数、累产油量、边际成本以及与原油中心处理站的地面距离权重，通过多元线性回归方法计算投资成本，计算公式如下：

[0153] CEXCPF＝0.0076+0.0027ωCPF+0.00063n+0.0897Oilproduct+0.00023ωborder

[0154] 式中，CAPEXCPF为地面总投资，ωCPF为断块油田与CPF的地面距离权重，n为井数，Oilproduct为累计产油量，ωborder为边际成本即计算的断块与大断块的距离，e、a、b、c、d为多
元回归系数；

[0155] 结合经济参数，建立NPV经济计算模型；基于开发年限内每年累计产油量计算NPV；NPV计算公式如公式(20)～(24)所示：

[0156]

[0157] 式(20)中，t为时间步长，1年；T表示总生产时间，35年；i为年折现率，2％；CPFt为时间t的运营现金流；CEX为投资金额；

[0158] 运营现金流CPFt计算公式如下：

[0159] CPFt＝INCOME‑OPEX               (21)

[0160] 式(21)中，INCOME为收入，OPEX为操作总投资；

[0161] 其中收入INCOME计算公式如下：

[0162]

[0163] 式(22)中，Rit为0.48；λ为98％； 表示t年的累计产油量；Poil为53.48美元/桶；操作总投资OPEX计算公式如下：

[0164]

[0165] 式(23)中，Re为0.76，n表示井数，μs为15.6万美元/桶，μv为5.45美元/桶，a为2％；

[0166] 投资金额CEX计算公式为：

[0167] CEX＝Re[n(CEXO+CEXD)(1+θt)+CEXCPF]    (24)

[0168] 式(24)中，n表示井数；CEXO为0.5百万美元；CEXD为3百万美元；θt为弃置率，表示当该年产油量达到了总体的75％以上θt为5％，反之θt为零；CEXCPF为地面总投资；

[0169] (6)利用所建开发年限内每年累计产油量的预测模型预测实际复杂分散断块油田群中各断块油田开发年限内每年的累计产油量，并依据所建NPV经济计算模型计算NPV；根
据NPV大小，完成复杂分散断块油田群投产顺序的决策。