基于支持向量机代理模型的油藏注采参数多目标优化方法与流程

本发明属于油藏模型模拟生产技术领域，具体涉及一种基于替代模型的油藏注采参数多目标优化方法。

背景技术：

油藏注水生产作为一种重要的提高采收率的方法已在世界大多数油田广泛使用，选取合适的注采参数对于高效开发油藏至关重要。

现有的常用注采参数优化方法是采用数值模拟技术或采用正交设计实验方法通过数值模拟技术或采用数值模拟加优化算法的方式对参数进行研究满足具体的开发要求,这种方法在注采参数优化过程中步骤繁琐、所需资料较多且需要花费大量的时间，同时，这种优化结果通常是静态的，得不到动态的优化结果。

技术实现要素：

针对上述针对水驱油藏生产过程中合适的注采参数选取难的问题，本发明提供一种基于支持向量机代理模型的油藏注采参数多目标优化方法，期许在最短的时间内按照目标函数的要求找出最优的方案。

具体的技术方案为：

基于支持向量机代理模型的油藏注采参数多目标优化方法，包括以下步骤：

(1)优化变量x的确立；将需要确定注采参数的n口注采井的生产时间分为l个时间段，从而建立起该优化设计问题需要优化的参数x，x为l*n维向量，确定各维度上优化变量的取值范围；

(2)在x各维度的取值范围内采用随机抽样的方法，得到r个优化变量x的样本，组成样本空间x，x＝[x1，x2，...，xr]；

(3)将样本空间x中的各优化变量xi；(i∈1，2，...r)带入油藏数值模拟器，求取初始优化变量对应的累积采油量np和净现值npv作为目标函数值；

(4)通过优化变量样本及样本对应的目标函数值，建立ls-svm替代模型，模型建立过程中，以模型和真实值之间误差为目标函数，采用粒子群算法寻优算法，获得ls-svm替代模型的参数；

ls-svm模型描述如下：

式中：f(x)为最优分类函数，αi为系数，xi为支持向量，k(xi，x)为核函数，为阈值，n为支持向量的个数；

(5)基于ls-svm替代模型生成大小为m的初始注采参数种群后，选取合适的pareto等级和拥挤距，然后采用带精英策略的非支配排序多目标优化遗传算法nsga-ii对注采参数进行优化，得出pareto解集；

(6)根据目标函数的需求在pareto解集中找出最优化的油藏注采参数。

本发明提供的基于支持向量机代理模型的油藏注采参数多目标优化方法，提出了以净现值和累产油量为目标函数的多目标优化注采参数设计方法，在保证代理模型预测精度的同时提高了效率。

附图说明

图1为本发明的流程图；

图2为实施例的区块渗透率分布图；

图3为实施例的区块孔隙度分布图；

图4为本发明实施例和对比方法两种方法pareto解集结果对比图；

图5a为实施例基础方案注采参数；

图5b为实施例方案1注采参数；

图5c为实施例方案2注采参数；

图5d为实施例方案3注采参数；

图6为实施例四种方案净现值；

图7为实施例四种方案现金流；

图8为实施例四种方案累产油量对比；

图9a原始方案生产时间150天含油饱和度分布；

图9b原始方案生产时间300天含油饱和度分布；

图9c原始方案生产时间450天含油饱和度分布；

图9d原始方案生产时间600天含油饱和度分布；

图10a优化方案生产时间150天含油饱和度分布；

图10b优化方案生产时间300天含油饱和度分布；

图10c优化方案生产时间450天含油饱和度分布；

图10d优化方案生产时间600天含油饱和度分布。

具体实施方式

结合附图说明本发明的具体技术方案。

本发明的方法如图1所示，将该方法应用于某区块两注两采模型，与eclipse油藏数值模拟软件结合带精英策略的非支配排序多目标优化遗传算法nsga-ⅱ相比较展现出了明显的优势。

选取研究区块油藏的两注两采模型作为优化的研究对象，生产井和注水井分布在方形油藏的四个角，以150天为一个阶段生产4个阶段共模拟600天的生产。区块含油饱和度为80％，渗透率和孔隙度分布如图2和图3所示。先用数值模拟软件计算出一定的样本，采用支持向量机形成替代数值模拟软件的代理模型，分别用代理模型和数值模拟软件结合nsga-ⅱ算法。

基于上述两种方法以净现值和累产油量为目标函数的油藏生产井井底压力和注水井注入量的pareto最优解集如图4所示，两种方案所用时间如下表所示。由下图可知两种方案所算出来的pareto最优解集在误差范围内重合，由此证明了基于粒子群算法的最小二乘支持向量机替代模型的多目标优化算法的有效性；表1反应了样本数量为50，做100代遗传优化两种方法所需要的数模次数和总时间，相比较于数值模拟方法，替代模型节省了近500倍的时间。由此证明了基于粒子群算法的最小二乘支持向量机替代模型的多目标优化算法的优越性。

表1两种方法所用时间对比表

在pareto最优解集中选取净现值最大、净现值趋近于8.5*10⁶元累产油量趋近于8.5*10⁴m³和累产油量最大的方案作为优化方案1、2、3，以井底压力15mpa和注入量300m³/天作为基础方案来对优化结果进行分析。

四种方案的注采参数分别如图5a到图5d所示。

四种方案的净现值和现金流随生产时间变换图如6和图7所示。以净现值最大这一单目标来说，优化方案1为最大净现值方案，按照优化方案生产会在600天内获得最大效益。

四种方案的累产油量随生产时间的变化如图8所示。以累产油量最大这一单目标来说，优化方案3为最大累产油量方案，按照优化方案生产会在600天内获得最大的产油量。

根据pareto最优解集的特性，对于多目标优化，净现值和累产油量这两个目标函数同时都很大是不可能的。从上述分析可以看出优化方案2的净现值和累产油量都相对较大，对于多目标优化是比较合适的解。相比较于没有做优化的方案来说，其可以很好的凸显出多目标优化方法的优势。基础方案和优化方案2的每个阶段含油饱和度分布如图9a到图9d、图10a到图10d所示。

由图9a到图9d、图10a到图10d可以看出当采用优化方案2生产相比于基准方案剩余油饱和度会更低，能达到更高的净现值和产出更多的油。