CN106502096B - 基于偏好多目标优化的油田机采过程决策参数优化方法 - Google Patents
基于偏好多目标优化的油田机采过程决策参数优化方法 Download PDFInfo
- Publication number
- CN106502096B CN106502096B CN201611001501.7A CN201611001501A CN106502096B CN 106502096 B CN106502096 B CN 106502096B CN 201611001501 A CN201611001501 A CN 201611001501A CN 106502096 B CN106502096 B CN 106502096B
- Authority
- CN
- China
- Prior art keywords
- variable
- state
- variables
- covariance
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000003921 oil Substances 0.000 claims abstract description 46
- 238000004519 manufacturing process Methods 0.000 claims abstract description 38
- 238000003062 neural network model Methods 0.000 claims abstract description 37
- 239000007788 liquid Substances 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 239000010721 machine oil Substances 0.000 claims abstract description 7
- 238000010606 normalization Methods 0.000 claims abstract description 4
- 239000002245 particle Substances 0.000 claims description 47
- 238000005070 sampling Methods 0.000 claims description 25
- 210000002569 neuron Anatomy 0.000 claims description 21
- 238000000605 extraction Methods 0.000 claims description 16
- 230000007613 environmental effect Effects 0.000 claims description 15
- 238000009826 distribution Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000002068 genetic effect Effects 0.000 claims description 7
- 238000013528 artificial neural network Methods 0.000 claims description 6
- 238000000513 principal component analysis Methods 0.000 claims description 5
- 239000013598 vector Substances 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 238000012952 Resampling Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 230000001351 cycling effect Effects 0.000 claims 1
- 238000012163 sequencing technique Methods 0.000 claims 1
- 238000011084 recovery Methods 0.000 abstract description 6
- 238000005265 energy consumption Methods 0.000 abstract description 2
- 239000012530 fluid Substances 0.000 description 7
- 230000000694 effects Effects 0.000 description 3
- 241001331845 Equus asinus x caballus Species 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 239000000725 suspension Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供的基于偏好多目标优化的油田机采过程决策参数优化方法，包括：确定油田机采油过程中的效率影响因素和性能变量；对样本中的载荷变量进行降维处理构建新样本，并归一化新样本；基于归一化后的新样本构建神经网络模型；利用ST‑UPFNN算法估计神经网络模型中由权值阈值所组成的状态变量的最优状态；并利用最优状态变量重构更新后的神经网络模型获得油田机采油过程模型；构建实际产液量的偏好函数；利用多目标进化算法对决策参量各自的上下限进行优化；将优化后的决策变量，带入油田机采油过程模型，计算优化后的决策变量的系统性能的平均值，与实际样本的系统性能的平均值进行比较。利用本发明可以提高油田机采油的生产效率，降低能耗。
Description
技术领域
本发明涉及油田机采技术领域，更为具体地，涉及一种基于偏好多目标优化的油田机采过程决策参数优化方法。
背景技术
油田机采油是一种机械采油方式，主要由电动机、地面传动设备和井下抽油设备三部分组成。油田机采油过程主要分为上、下两个冲程，上冲程，即驴头悬点向上运动，需提起抽油杆柱和液柱，电动机需消耗大量的能量；下冲程，即驴头悬点向下运动，油田机杆柱转拉动对电动机做功。在杆柱上下运动过程中，液柱负载发生周期性变化，使得油田机系统在电机做功、传动装置等方面能耗较大，以致系统工作效率低下。
发明内容
鉴于上述问题，本发明的目的是提供一种基于偏好多目标优化的油田机采过程决策参数优化方法，以解决上述背景技术所提出的问题。
本发明提供的基于偏好多目标优化的油田机采过程决策参数优化方法，包括：
步骤S1：确定油田机采油过程中的效率影响因素,构成效率观测变量集合{x1,x2,x3,L xn}；以及，选取油田机工艺系统的性能变量，构成性能观测变量集合{y1,y2}；
其中，x1为冲次决策变量，x2为有效冲程决策变量，x3～x5分别为计算泵效环境变量、含水率环境变量，平均功率因数环境变量，x6～xn均为载荷环境变量；性能观测变量的个数l＝2，y1为日产液量，y2为日耗电量；
步骤S2：根据效率观测变量集合{x1,x2,x3,L xn}和性能观测变量集合{y1,y2}，采集通过ST-UPFNN算法构建神经网络模型的观测变量的样本值矩阵[x1，x2L xn，y1，y2]；其中，
设定采样周期为T，在采集观测变量的过程中，如果采样周期小于T，对T周期内的样本求取平均值以作为该T周期的样本[I,Y]；如果采样周期大于T，剔除采集到的观测变量，将样本中的I作为输入样本，将样本中的Y作为输出样本；
步骤S3：利用主元分析算法对载荷环境变量进行降维，构建新的载荷主元变量{Lz1,Lz2,...,Lzd}；
其中，构建新的载荷主元变量{Lz1,Lz2,...,Lzd}为d个载荷主元分量，每个载荷主元分量的维度与所述样本[I,Y]的数量相同；
步骤S4：重新组合非载荷变量与d个载荷主元分量，构建新的输入样本I1，并对新的输入样本I1和输出样本Y进行归一化，获得归一化后的样本其属于[-1,1]；其中，非载荷变量包括冲次决策变量x1、有效冲程决策变量x2、计算泵效环境变量x3、含水率环境变量x4、平均功率因数环境变量x5；
步骤S5：基于归一化后的样本构建神经网络模型和神经网络模型的初始状态变量X，以及，将归一化后的样本中的作为神经网络模型的输入，将归一化后的样本中的作为神经网络模型的输出；
其中，神经网络模型为：
其中，Ik为训练样本的矢量样本值，并作为神经网络模型的输入，为网络输入层到隐含层的神经元的连接权值，为网络输入层到隐含层的神经元的阈值，为隐含层到网络输出层的神经元的连接权值，为隐含层到网络输出层的神经元的阈值，其中，i＝1，2…S0；j＝1，2…S1；k＝1，2…S2；S0为网络输入层的神经元的数量，S1为网络隐含层的神经元的数量，S2为网络输出层的神经元的数量；
初始状态变量X为：
步骤S6：利用ST-UPFNN算法估计所述神经网络模型的最优状态变量；
步骤S7：将最优状态变量作为所述神经网络模型的和重构神经网络表达式，获得油田机采油过程模型；
步骤S8：构建日产液量y1的偏好函数perfc(y)；
步骤S9：利用NSGA-Ⅱ算法对日耗电量y2的决策参量各自的上下限进行优化；
步骤S10：将优化后的决策变量结合环境变量，带入由ST-UPFNN算法所建立的油田机采油过程模型，在油田机采油过程模型的基础之上进行优化，得到优化后的决策变量的系统性能的平均值，与实际样本的系统性能的平均值进行比较，如果优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值，利用优化后的决策变量对实际生产进行指导；否则重复上述步骤S1-S9，直至优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值为止。
本发明提供的基于偏好多目标优化的油田机采过程决策参数优化方法，通过ST-UPFNN算法挖掘油田机的生产规律，并利用多目标进化算法优化油田机生产过程决策参量，提高油田机的生产效率。
附图说明
通过参考以下结合附图的说明及权利要求书的内容，并且随着对本发明的更全面理解，本发明的其它目的及结果将更加明白及易于理解。在附图中：
图1为主原分量的贡献率柱状图；
图2为产液量的偏好函数图；
图3为产液量偏好值与耗电量的Pareto解集图；
图4为产液量实际值与耗电量的Pareto解集图。
具体实施方式
名称解释
ST-UKFNN：Strong TrackUnscented Kalman FilterNeural Network，强追踪无迹卡尔曼滤波神经网络；
ST-UPFNN：Strong TrackUnscented Particle FilterNeural Network，强追踪无迹粒子滤波神经网络，其将ST-UKFNN、粒子滤波(Particle Filter)、BP神经网络相结合。
本发明提供的基于偏好多目标优化的油田机采过程决策参数优化方法，包括：
步骤S1：确定油田机采油过程中的效率影响因素，构成效率观测变量集合{x1,x2,x3,L xn}；以及，选取油田机工艺系统的性能变量，构成性能观测变量集合{y1,y2}。
其中，x1为冲次决策变量，x2为有效冲程决策变量，x3～x5分别为计算泵效环境变量、含水率环境变量，平均功率因数环境变量，x6～xn均为载荷环境变量；性能观测变量的个数l＝2，y1为日产液量，y2为日耗电量。
本发明中，选取性能影响因素与性能指标如表1所示：
表1
| 变量类型 | 变量名称 |
| 决策变量 | 冲次 |
| 决策变量 | 有效冲程 |
| 环境变量 | 计算泵效 |
| 环境变量 | 含水率 |
| 环境变量 | 平均功率因数 |
| 环境变量 | 载荷 |
| 输出变量 | 日产液量 |
| 输出变量 | 日耗电量 |
步骤S2：根据效率观测变量集合{x1,x2,x3,L xn}和性能观测变量集合{y1,y2}，采集通过ST-UPFNN算法构建神经网络模型的观测变量的样本值矩阵[x1，x2L xn，y1，y2]。
设定采样周期为T，在采集观测变量的过程中，如果采样周期小于T，对T周期内的样本求取平均值以作为该T周期的样本[I,Y]，即[I,Y]通过神经网络模型的观测变量的样本值矩阵[x1，x2L xn，y1，y2]经平均值后得到；如果采样周期大于T，说明存在样本不足的现象，直接剔除采集到的观测变量，将样本中的I作为输入样本，将样本中的Y作为输出样本。
样本[I,Y]如表2所示：
表2
| 参量 | 1 | 2 | 3 | …… | 256 |
| 冲次/(min<sup>-1</sup>) | 3.12 | 3.12 | 3.1 | …… | 4.19 |
| 有效冲程/(m) | 3.49 | 3.51 | 3.5 | …… | 3.43 |
| 平均功率因数 | 0.26 | 0.25 | 0.25 | …… | 0.67 |
| 计算泵效/％ | 92.08 | 95.82 | 94.51 | …… | 94.05 |
| 含水率/％ | 93.7 | 93.7 | 93.5 | …… | 93.4 |
| 载荷1/(KN) | 29.47 | 29.72 | 29.59 | …… | 33.99 |
| 载荷2/(KN) | 30.75 | 30.74 | 30.66 | …… | 34.23 |
| …… | …… | …… | …… | …… | …… |
| 载荷144/(KN) | 29.88 | 30.13 | 30.03 | …… | 33.59 |
| 产液量/(t/d) | 40.95 | 41.02 | 41.25 | …… | 54.91 |
| 耗电量/(kw·h) | 129 | 130.1 | 129 | …… | 154.1 |
步骤S3：利用主元分析算法对载荷环境变量进行降维，构建新的载荷主元变量{Lz1,Lz2,...,Lzd}。
其中，构建新的载荷主元变量{Lz1,Lz2,...,Lzd}为d个载荷主元分量，每个载荷主元分量的维度与所述样本[I,Y]的数量相同；
本发明采用示功图描绘数据的144个载荷点作为部分环境变量进行建立神经网络模型，利用144维数据建立神经网络模型为参数维度灾难。故而利用主元分析算法(Principal ComponentAnalysis，PCA)对载荷环境变量进行降维处理，构建新的载荷主元变量，新的载荷主元变量构成的集合:{Lz1,Lz2,...,Lzd}，其为d个载荷主元分量，每个主元分量维度与样本[X,Y]的数量相同。令功图数据为：设置样本累计贡献率precent＝0.90；如图1所示，得到前5个主元分量的贡献率以及累计贡献率。
故此，取前2个主元分量B1、B2作为载荷环境变量的特征变量，其部分值如下表所示：
表3部分主元分量数据
| B1 | B2 |
| 399.31 | 23.85 |
| 399.28 | 25.86 |
| 399.11 | 24.69 |
| 401.31 | 25.52 |
| …… | …… |
| 454.67 | 42.03 |
| 448.12 | 37.89 |
步骤S4：重新组合非载荷变量与d个载荷主元分量，构建新的输入样本I1，并对新的输入样本I1和输出样本Y进行归一化，获得归一化后的样本其属于[-1,1]；其中，非载荷变量包括冲次决策变量x1、有效冲程决策变量x2、计算泵效环境变量x3、含水率环境变量x4、平均功率因数环境变量x5。
步骤S5：基于归一化后的样本构建神经网络模型和神经网络模型的初始状态变量X，以及，将归一化后的样本中的作为神经网络模型的输入，将归一化后的样本中的作为神经网络模型的输出；
其中，构建的神经网络模型为：
其中，Ik为训练样本的矢量样本值，并作为神经网络模型的输入，为网络输入层到隐含层的神经元的连接权值，为网络输入层到所述隐含层的神经元的阈值，为隐含层到网络输出层的神经元的连接权值，为隐含层到网络输出层的神经元的阈值，其中，i＝1，2…S0；j＝1，2…S1；k＝1，2…S2；S0为网络输入层的神经元的数量，S1为网络隐含层的神经元的数量，S2为网络输出层的神经元的数量；
构建的初始状态变量为：
步骤S6：利用ST-UPFNN算法估计神经网络模型的状态变量X，以得到最优状态变量，完成所建模型的权值阈值更新，使得所得到的模型更符合实际生产过程。
利用ST-UPFNN算法估计神经网络模型的最优状态变量的过程，包括：
步骤S61：针对粒子滤波器设置粒子的数目N，并以x0为均值,P0为方差进行正态分布采样，得到初始粒子集并将初始粒子集中的每个粒子的权值均设为1/N；记粒子x0为k＝0时刻状态；
步骤S62：在获取(k+1)时刻的观测变量值后，为归一化样本中第(1)组样本性能观测变量，利用ST-UKFNN算法对每个粒子进行状态估计，得到最优状态估计值和协方差
利用ST-UKFNN算法对每个粒子进行状态估计的过程如下：
步骤S621：对初始状态变量X进行Sigma采样，获得2n+1个采样点，初始化控制2n+1个采样点的分布状态参数α、待选参数κ，以及非负权系数β，对所述初始状态变量X的Sigma采样如下：
步骤S622：计算每个采样点的权重，每个采样点的权重如下：
其中，Wc为计算状态变量的协方差的权重，Wm为计算状态估计和观测预测时的权重，是的第一列，是的第一列；
步骤S623：通过离散时间非线性系统的状态方程将每个采样点的k时刻的最优状态变量的状态估计变换为(k+1)时刻的状态变量的状态估计以及，通过合并(k+1)时刻的状态估计的向量，获得(k+1)时刻的状态变量的状态先验估计和协方差Pk+1|k；其中，
(k+1)时刻的状态变量的状态估计为：
其中，为k时刻的最优状态估计，wk为过程噪声，其协方差矩阵Qk为cov(wk,wj)＝Qkδkj，
(k+1)时刻的状态变量的状态先验估计为：
(k+1)时刻的状态变量的协方差Pk+1|k为：
步骤S624：通过离散时间非线性系统的观测方程建立(k+1)时刻的状态变量的状态估计和(k+1)时刻的观测预测的联系：
其中，νk为观测噪声，其协方差矩阵Rk为cov(vk,vj)＝Rkδkj，
步骤S625：通过估计(k+1)时刻的观测预测的向量，获得(k+1)时刻的先验观测预测并根据先验观测预测估计(k+1)时刻的观测预测的协方差
(k+1)时刻的先验观测预测为：
(k+1)时刻的观测预测的协方差为：
其中，在此处引入强追踪算法，即渐消因子λk+1增强模型的追踪能力以提高模型精度；
其中，β为弱化因子，β≥1；
步骤S626：计算(k+1)时刻的状态变量的状态先验估计与(k+1)时刻的先验观测预测之间的协方差
协方差为：
步骤S627：通过建立协方差和预测协方差的关系，以更新(k+1)时刻的状态变量的状态估计和协方差，分别获得(k+1)时刻的最优状态估计值和协方差
建立的协方差和协方差的关系为：
其中，Kk+1为增益矩阵，以及更新后的(k+1)时刻的状态变量的状态估计协方差Pk+1为：
将更新后的(k+1)时刻的状态变量的状态估计和协方差Pk+1分别作为(k+1)时刻的最优状态估计值和协方差
步骤S63：将最优状态估计值和协方差作为粒子的重要性密度函数进行抽样，得到新粒子由所有新粒子组成的粒子集中的每个新粒子的正态分布概率密度值如下：
其中，p为每个新粒子的条件概率，randnorm为正态分布随机误差，正态分布密度函数：x、μ、σ分别为正态分布函数的三个变量；
对于x,μ,σ分别与一一对应；
对于x,μ,σ分别与一一对应；
对于x,μ,σ分别与一一对应。
步骤S64：对新粒子的权值进行更新，并进行归一化处理。
权值更新公式为：
权值归一化公式为：
步骤S65：根据粒子权值和重采样策略对粒子集进行重采样，从而获取新粒子集并求取新粒子集中每个新粒子的状态估计值
设变量u，令取u1∈(0,1)
步骤S66：以粒子的数目N作为循环次数循环步骤S61-步骤S65的计算过程，将最后一次估计得到系统状态变量作为利用ST-UPFNN算法估计得到的所述神经网络模型的最优状态变量；其中，将新粒子的状态估计值作为本时刻的最优估计赋给进行下一时刻的状态估计。
神经网络模型的最优状态变量的结构参数如下：
步骤S7：将最优状态变量作为神经网络模型的和重构神经网络表达式，获得油田机采油过程模型。
步骤S8：构建日产液量y1的偏好函数perfc(y)。
在系统工艺参数优化计算中，考虑对不同参数具有不同的喜好程度，利用物理规划构建系统偏好函数。设定产液量最优值为y1best，设定值为ybest，在设定值ybest周围某一邻域范围[ybest-△y,ybest+△y]内波动为非常满意(HD)，且在[ybest-△y-△y1,ybest-△y],[ybest+△y,ybest+△y+△y1]内为满意(D)，依次得到可接受(T)，不满意(U)和非常不满意(HU)，对应的偏好值区间用[0,2]，[2,4]，[4,6]，[6,8]，[8,10]表示。
假定将所有样本的平均日产液量作为给定产液量及极好程度的偏好值(47.38)。同时设定所有日产液量数据的最小值(40.22)和最大值(56.92)作为不可接受域的临界值。故而设计偏好程度区间为：[0,2],[2,4],[4,6],[6,8],[8,10]等，且设计的偏好区间边界值所对应的实际产液量区间的边界值如表4所示，偏好函数如图2所示。
表4偏好函数的边界值对应表
拟合得到日产液量的偏好函数为：
步骤S9：利用NSGA-Ⅱ算法(Non-dominated Sorting Genetic Algorithm-Ⅱ，多目标优化算法)对日耗电量y2的决策参量各自的上下限进行优化。
优化的过程，包括：
步骤S91：通过决策变量个体P＝[x1x2L xn]的非支配性比较寻找最佳个体；其中，个体非支配性由适应度函数值以及个体拥挤度共同决定。
多目标优化适应度函数为：
其中，是由ST-UPFNN算法构建的油田机采油过程模型：
结合向最小值方向搜索最值的特点，将部分求最大值的性能变量函数取反，从而得到多目标优化适应度函数。
本发明中，由于优化计算的过程中，得到的产液量越接近最佳值越好；耗电量越低越好，故而将净化气产量的预测函数取反函数作为该性能变量的多目标优化适应度函数。
步骤S92：计算油田机工艺系统的环境变量的平均值：
其中，环境变量包括计算泵效环境变量x3、含水率环境变量x4、平均功率因数环境变量x5，N为环境变量的输入样本的数量。
表5环境变量平均值表
步骤S93：利用决策参量(x1,x2)构建父代种群P和子代种群Q，其中，
其中，K1为父代种群P中的个体(1≤m≤L)的数量；K2为子代种群Q中的个体(1≤m≤K)的数量；L为初始化的种群样本数量，L＝50；GEN为最大遗传代数，GEN＝100。
步骤S94：根据决策参量的上下限xi,min≤xi≤xi,max(i＝1,2,L,n)初始化父种群P。
其中，初始化父种群P的过程为：从x1的取值范围内随机取值赋予 从x2的取值范围内，随机取值赋予
步骤S95：进行第一次遗传迭代GEN＝1，从父代种群P中选出任意对个体，对每对个体(1≤m≤K,1≤n≤K)进行遗传交叉、变异计算，并将计算结果赋予子代种群Q中相应的一对个体将父代种群P与子代种群Q进行合并得到种群R。
步骤S96：从种群R中选择精英个体作为下一代种群；其中，将种群R中所有个体按照支配性层级进行排序，并计算所有个体的拥挤度，按照种群R中的个体按照层级从低到高的顺序从非支配集中选择个体作为精英个体，而对于同一层级的非支配集，按照拥挤度从大到小的顺序选择个体作为精英个体，保持种群个体为L。
步骤S97：循环GEN次步骤S94～步骤S96，得到第GEN代种群作为优化结果输出；其中，
得到Pareto解集，产液量偏好函数与耗电量Pareto解集如图3所示，产液量实际值与耗电量的Pareto解集如图4所示。
由优化所得Pareto解集分析可知优化前后效果对比如表6所示：
表6优化参数对应目标函数与生产实际输出对比表
优化平均产耗比提升约为3％，达到了节能增效的优化效果，说明本次结果有效。
步骤S10：将优化后的决策变量结合环境变量，带入ST-UPFNN算法建立的油田机采油过程模型，计算优化后的决策变量的系统性能，与实际样本的系统性能的平均值进行比较，如果优化后的决策变量的系统性能大于实际样本的系统性能的平均值，利用优化后的决策变量对实际生产进行指导；否则重复上述步骤S1-S9，直至优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值为止。
产液量越接近最优值、耗电量越低则效果越优越好。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (5)
1.一种基于偏好多目标优化的油田机采过程决策参数优化方法，包括：
步骤S1：确定油田机采油过程中的效率影响因素,构成效率观测变量集合{x1,x2,x3,Lxn}；以及，选取油田机工艺系统的性能变量，构成性能观测变量集合{y1,y2}；
其中，x1为冲次决策变量，x2为有效冲程决策变量，x3～x5分别为计算泵效环境变量、含水率环境变量，平均功率因数环境变量，x6～xn均为载荷环境变量；性能观测变量的个数l＝2，y1为日产液量，y2为日耗电量；
步骤S2：根据效率观测变量集合{x1,x2,x3,L xn}和性能观测变量集合{y1,y2}，采集通过ST-UPFNN算法构建神经网络模型的观测变量的样本值矩阵[x1，x2L xn，y1，y2]；其中，
设定采样周期为T，在采集观测变量的过程中，如果采样周期小于T，对T周期内的样本求取平均值以作为该T周期的样本[I,Y]；如果采样周期大于T，剔除采集到的观测变量；其中，将样本中的I作为输入样本，将样本中的Y作为输出样本；
步骤S3：利用主元分析算法对载荷环境变量进行降维，构建新的载荷主元变量{Lz1,Lz2,...,Lzd}；
其中，构建新的载荷主元变量{Lz1,Lz2,...,Lzd}为d个载荷主元分量，每个载荷主元分量的维度与所述样本[I,Y]的数量相同；
步骤S4：重新组合非载荷变量与d个载荷主元分量，构建新的输入样本I1，并对新的输入样本I1和输出样本Y进行归一化，获得归一化后的样本其属于[-1,1]；其中，非载荷变量包括冲次决策变量x1、有效冲程决策变量x2、计算泵效环境变量x3、含水率环境变量x4、平均功率因数环境变量x5；
步骤S5：基于所述归一化后的样本构建神经网络模型和所述神经网络模型的初始状态变量X，以及，将所述归一化后的样本中的作为所述神经网络模型的输入，将所述归一化后的样本中的作为所述神经网络模型的输出；
其中，所述神经网络模型为：
其中，Ik为训练样本的矢量样本值，并作为所述神经网络模型的输入，为网络输入层到隐含层的神经元的连接权值，为网络输入层到所述隐含层的神经元的阈值，为所述隐含层到网络输出层的神经元的连接权值，为所述隐含层到所述网络输出层的神经元的阈值，其中，i＝1，2…S0；j＝1，2…S1；k＝1，2…S2；S0为所述网络输入层的神经元的数量，S1为所述网络隐含层的神经元的数量，S2为所述网络输出层的神经元的数量；
所述初始状态变量X为：
步骤S6：利用ST-UPFNN算法估计所述神经网络模型的最优状态变量；
步骤S7：将所述最优状态变量作为所述神经网络模型的和重构神经网络表达式，获得油田机采油过程模型；
步骤S8：构建日产液量y1的偏好函数perfc(y)；
步骤S9：利用NSGA-Ⅱ算法对日耗电量y2的决策参量各自的上下限进行优化；
步骤S10：将优化后的决策变量结合环境变量，带入由ST-UPFNN算法所建立的油田机采油过程模型，在所述油田机采油过程模型的基础之上进行优化，计算优化后的决策变量的系统性能的平均值，与实际样本的系统性能的平均值进行比较，如果优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值，利用优化后的决策变量对实际生产进行指导；否则重复上述步骤S1-S9，直至优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值为止。
2.如权利要求1所述的基于偏好多目标优化的油田机采过程决策参数优化方法，所述步骤S6包括：
步骤S61：针对粒子滤波器设置粒子的数目N，并以x0为均值,P0为方差进行正态分布采样，得到初始粒子集并将所述初始粒子集中的每个粒子的权值均设为1/N；记粒子x0为k＝0时刻状态；
步骤S62：在获取k+1时刻的观测变量值后，为归一化样本中第(1)组样本性能观测变量，利用ST-UKFNN算法对每个粒子进行状态估计，得到最优状态估计值和协方差其中，k＝0；
步骤S63：将所述最优状态估计值和协方差作为粒子的重要性密度函数进行抽样，得到新粒子由所有新粒子组成的粒子集中的每个新粒子的正态分布概率密度值如下：
其中，p为每个新粒子的条件概率，randnorm为正态分布随机误差，正态分布密度函数：x、μ、σ分别为正态分布的三个变量；
步骤S64：对新粒子的权值进行更新，并进行归一化处理；其中，
权值更新公式为：
权值归一化公式为：
步骤S65：根据粒子权值和重采样策略对粒子集进行重采样，从而获取新粒子集并求取新粒子集中每个新粒子的状态估计值
设变量u，令取u1∈(0,1)
步骤S66：以粒子的数目N作为循环次数循环步骤S61-步骤S65的计算过程，将最后一次估计得到系统状态变量作为利用ST-UPFNN算法估计得到的所述神经网络模型的最优状态变量；其中，将新粒子的状态估计值作为本时刻的最优估计赋给进行下一时刻的状态估计。
3.如权利要求2所述的基于偏好多目标优化的油田机采过程决策参数优化方法，利用ST-UKFNN算法对每个粒子进行状态估计的过程，包括：
步骤S621：对所述初始状态变量X进行Sigma采样，获得2n+1个采样点，初始化控制2n+1个采样点的分布状态参数α、待选参数κ，以及非负权系数β，对所述初始状态变量X的Sigma采样如下：
步骤S622：计算每个采样点的权重，每个采样点的权重如下：
其中，Wc为计算状态变量的协方差的权重，Wm为计算状态估计和观测预测时的权重，是的第一列，是的第一列；
步骤S623：通过离散时间非线性系统的状态方程将每个采样点的k时刻的最优状态变量的状态估计变换为(k+1)时刻的状态变量的状态估计以及，通过合并(k+1)时刻的状态估计的向量，获得(k+1)时刻的状态变量的状态先验估计和协方差Pk+1|k；其中，
所述状态估计为：
其中，为k时刻的最优状态估计，wk为过程噪声，其协方差矩阵Qk为cov(wk,wj)＝Qkδkj，
所述状态先验估计为：
所述状态变量的协方差Pk+1|k为：
步骤S624：通过离散时间非线性系统的观测方程建立(k+1)时刻的状态变量的状态估计和(k+1)时刻的观测预测的联系：
其中，νk为观测噪声，其协方差矩阵Rk为cov(vk,vj)＝Rkδkj，
步骤S625：通过估计(k+1)时刻的观测预测的向量，获得(k+1)时刻的先验观测预测并根据先验观测预测估计(k+1)时刻的观测预测的协方差
(k+1)时刻的先验观测预测以及观测预测的协方差分别为：
步骤S626：计算(k+1)时刻的状态变量的状态先验估计与(k+1)时刻的先验观测预测之间的协方差所述协方差为：
步骤S627：通过建立协方差和协方差的关系，更新(k+1)时刻的状态变量的状态估计和协方差，分别获得(k+1)时刻的最优状态估计值和协方差
4.如权利要求3所述的基于偏好多目标优化的油田机采过程决策参数优化方法，其中，建立的协方差和协方差的关系为：
其中，Kk+1为增益矩阵，以及更新后的(k+1)时刻的状态变量的状态估计协方差Pk+1为：
将更新后的(k+1)时刻的状态变量的状态估计和协方差Pk+1分别作为(k+1)时刻的最优状态估计值和协方差
5.如权利要求1所述的基于偏好多目标优化的油田机采过程决策参数优化方法，其中，步骤S9包括：
步骤S91：通过决策变量个体P＝[x1x2L xn]的非支配性比较寻找最佳个体；其中个体非支配性由多目标适应度函数值以及个体拥挤度共同决定，多目标优化适应度函数为：
其中，是由ST-UPFNN算法构建的油田机采油过程模型：
步骤S92：计算油田机工艺系统的环境变量的平均值：
其中，所述环境变量包括计算泵效环境变量x3、含水率环境变量x4、平均功率因数环境变量x5，N为所述环境变量的输入样本的数量；
步骤S93：利用决策参量(x1,x2)构建父代种群P和子代种群Q，其中，
其中，K1为父代种群P中的个体的数量，1≤m≤L；K2为子代种群Q中的个体的数量，1≤m≤K；L为初始化的种群样本数量，L＝50；GEN为最大遗传代数，GEN＝100；
步骤S94：根据决策参量的上下限xi,min≤xi≤xi,max，初始化父种群P；其中，i＝1,2,L,n；以及，初始化父种群P的过程为：从x1的取值范围内随机取值赋予1≤m≤K，从x2的取值范围内，随机取值赋予1≤m≤K；
步骤S95：进行第一次遗传迭代GEN＝1，从父代种群P中选出任意对个体，对每对个体进行遗传交叉、变异计算，1≤m≤K，1≤n≤K，并将计算结果赋予子代种群Q中相应的一对个体将父代种群P与子代种群Q进行合并得到种群R：
步骤S96：从种群R中选择精英个体作为下一代种群；其中，将种群R中所有个体按照支配性层级进行排序，并计算所有个体的拥挤度，按照种群R中的个体按照层级从低到高的顺序从非支配集中选择个体作为精英个体，而对于同一层级的非支配集，按照拥挤度从大到小的顺序选择个体作为精英个体，保持种群个体为L；
步骤S97：GEN＝GEN+1，循环100次步骤S94～步骤S96，得到第GEN代种群作为优化结果输出；其中，
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611001501.7A CN106502096B (zh) | 2016-11-14 | 2016-11-14 | 基于偏好多目标优化的油田机采过程决策参数优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611001501.7A CN106502096B (zh) | 2016-11-14 | 2016-11-14 | 基于偏好多目标优化的油田机采过程决策参数优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106502096A CN106502096A (zh) | 2017-03-15 |
| CN106502096B true CN106502096B (zh) | 2019-06-14 |
Family
ID=58324461
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611001501.7A Active CN106502096B (zh) | 2016-11-14 | 2016-11-14 | 基于偏好多目标优化的油田机采过程决策参数优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106502096B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107862129B (zh) * | 2017-11-03 | 2021-02-02 | 哈尔滨工业大学 | 一种基于moead的偏差区间偏好引导多目标决策优化方法 |
| CN107748937B (zh) * | 2017-11-03 | 2021-08-13 | 哈尔滨工业大学 | 一种基于moead的比例区间偏好引导多目标决策优化方法 |
| CN109085752B (zh) * | 2018-03-09 | 2020-09-29 | 重庆科技学院 | 基于角度支配关系的铝电解偏好多目标优化算法 |
| CN110598163A (zh) * | 2019-09-09 | 2019-12-20 | 北京邮电大学 | 基于双精英差分进化的表面活性剂驱动态优化方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104680259A (zh) * | 2015-03-12 | 2015-06-03 | 天津市万众科技发展有限公司 | 抽油机运行参数的节能优化方法 |
| CN104680025A (zh) * | 2015-03-12 | 2015-06-03 | 重庆科技学院 | 基于遗传算法极限学习机的抽油机参数优化方法 |
| CN104680023A (zh) * | 2015-03-12 | 2015-06-03 | 重庆科技学院 | 基于多目标决策的抽油机参数优化方法 |
| CN104680257A (zh) * | 2015-03-12 | 2015-06-03 | 天津市万众科技发展有限公司 | 面向抽油机工作过程的智能优化方法 |
| CN105045941A (zh) * | 2015-03-13 | 2015-11-11 | 重庆科技学院 | 基于无迹卡尔曼滤波的抽油机参数优化方法 |
| CN105046326A (zh) * | 2015-03-13 | 2015-11-11 | 重庆科技学院 | 基于功图主元分析的抽油机参数优化方法 |
-
2016
- 2016-11-14 CN CN201611001501.7A patent/CN106502096B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104680259A (zh) * | 2015-03-12 | 2015-06-03 | 天津市万众科技发展有限公司 | 抽油机运行参数的节能优化方法 |
| CN104680025A (zh) * | 2015-03-12 | 2015-06-03 | 重庆科技学院 | 基于遗传算法极限学习机的抽油机参数优化方法 |
| CN104680023A (zh) * | 2015-03-12 | 2015-06-03 | 重庆科技学院 | 基于多目标决策的抽油机参数优化方法 |
| CN104680257A (zh) * | 2015-03-12 | 2015-06-03 | 天津市万众科技发展有限公司 | 面向抽油机工作过程的智能优化方法 |
| CN105045941A (zh) * | 2015-03-13 | 2015-11-11 | 重庆科技学院 | 基于无迹卡尔曼滤波的抽油机参数优化方法 |
| CN105046326A (zh) * | 2015-03-13 | 2015-11-11 | 重庆科技学院 | 基于功图主元分析的抽油机参数优化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106502096A (zh) | 2017-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106502096B (zh) | 基于偏好多目标优化的油田机采过程决策参数优化方法 | |
| CN105045941B (zh) | 基于无迹卡尔曼滤波的抽油机参数优化方法 | |
| CN105046326B (zh) | 基于功图主元分析的抽油机参数优化方法 | |
| Zha et al. | Forecasting monthly gas field production based on the CNN-LSTM model | |
| CN106529042B (zh) | 基于计算智能的油田机采参数动态演化建模与优化方法 | |
| CN106600059B (zh) | 基于改进rbf神经网络的智能电网短期负荷预测方法 | |
| CN103345559B (zh) | 铝电解过程电解槽工艺能耗的动态演化建模方法 | |
| CN106773667B (zh) | 基于无迹粒子滤波神经网络的油田机采参数建模方法 | |
| CN111931418B (zh) | 一种基于进化神经网络的马蹄焰玻璃窑炉的热效率预测方法 | |
| Wang et al. | Parameter estimation of a bioreaction model by hybrid differential evolution | |
| CN105787259A (zh) | 一种多元气象因素对负荷变化影响相关性的分析方法 | |
| CN106951995A (zh) | 一种超高压输电电场极限学习机预测多目标优化屏蔽方法 | |
| CN104680023B (zh) | 基于多目标决策的抽油机参数优化方法 | |
| CN104680025A (zh) | 基于遗传算法极限学习机的抽油机参数优化方法 | |
| CN114862015A (zh) | 一种基于FEN-ConvLSTM模型的台风风速智能预测方法 | |
| CN104680257A (zh) | 面向抽油机工作过程的智能优化方法 | |
| CN105426965A (zh) | 一种应用于多目标区间遗传算法的排序方法 | |
| Bi et al. | Self-adaptive Teaching-learning-based Optimizer with Improved RBF and Sparse Autoencoder for Complex Optimization Problems | |
| Irwandi et al. | Performance analysis of subtractive clustering algorithm in determining the number and position of cluster centers | |
| CN106530130A (zh) | 油田机采过程动态演化建模与节能优化方法 | |
| CN104680259A (zh) | 抽油机运行参数的节能优化方法 | |
| CN107766887A (zh) | 一种局部加权的不完整数据混杂聚类方法 | |
| CN110175698B (zh) | 一种基于改进粒子群优化bp神经网络的水泵设备状态预测方法 | |
| D’Ambrosio et al. | Meta-model assisted evolutionary optimization of cellular automata: An application to the SCIARA model | |
| CN110909922B (zh) | 一种水资源效率检测及预测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |