CN105046326A - 基于功图主元分析的抽油机参数优化方法 - Google Patents
基于功图主元分析的抽油机参数优化方法 Download PDFInfo
- Publication number
- CN105046326A CN105046326A CN201510111157.6A CN201510111157A CN105046326A CN 105046326 A CN105046326 A CN 105046326A CN 201510111157 A CN201510111157 A CN 201510111157A CN 105046326 A CN105046326 A CN 105046326A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- mover
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000005086 pumping Methods 0.000 title claims abstract description 40
- 238000005457 optimization Methods 0.000 title claims abstract description 16
- 238000010586 diagram Methods 0.000 title claims abstract description 15
- 238000000513 principal component analysis Methods 0.000 title claims abstract description 14
- 238000004519 manufacturing process Methods 0.000 claims abstract description 32
- 238000013528 artificial neural network Methods 0.000 claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 238000012549 training Methods 0.000 claims abstract description 18
- 238000010606 normalization Methods 0.000 claims abstract description 17
- 239000007788 liquid Substances 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 54
- 230000008569 process Effects 0.000 claims description 22
- 210000002569 neuron Anatomy 0.000 claims description 17
- 230000007613 environmental effect Effects 0.000 claims description 15
- 238000001914 filtration Methods 0.000 claims description 10
- 238000000605 extraction Methods 0.000 claims description 9
- 230000009467 reduction Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 7
- 230000002068 genetic effect Effects 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000004087 circulation Effects 0.000 abstract 1
- 238000000638 solvent extraction Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000005265 energy consumption Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000004134 energy conservation Methods 0.000 description 2
- 239000000725 suspension Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000012847 principal component analysis method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于功图主元分析的抽油机参数优化方法,包括:1)确定抽油机的效率影响因素和性能变量、2)获得其样本数据、3)对载荷数据进行降维、4)由非载荷变量与载荷新主元构建网络输入变量、5)构建输入变量的样本值、6)归一化处理、7)构建前馈神经网络、8)利用无迹卡尔曼滤波对网络进行训练、9)构造父代和子代种群;10)对父代个体作遗传变异计算,以产生子代个体;11)对父代和子代个体求适应度函数;12)将父代和子代个体划分到层级不同的非支配集中;13)从这些非支配集中选择个体构成新的父代种群,循环10)-13)多次,得到优化后的效率影响因素值。优化后,可得到产液量最大时,耗电量最小。
Description
技术领域
本发明属于采油领域,具体涉及一种基于功图主元分析的抽油机参数优化方法。
背景技术
抽油机采油作为一种机械采油方式,主要由电动机、地面传动设备和井下抽油设备三部分组成,如图1所示。抽油机的整个采油过程主要分为上下两个冲程:上冲程,即驴头悬点向上运动,提起抽油杆柱和井下抽油设备,此过程中电动机需消耗大量的能量;下冲程,即驴头悬点向下运动,抽油机的抽油杆柱电动机做功。在抽油杆柱上下运动过程中,电动机的负载发生周期变化。抽油机的运行参数的选择对整个抽油机系统的能量消耗影响很大。为了使抽油机采油生产过程既能完成预定的产液量,又能使抽油机生产过程的耗电量最低,需要对抽油机运行参数进行节能优化。
发明内容
本发明是为了解决现有技术中存在的上述技术问题而做出,其目的在于提供一种基于功图主元分析的抽油机参数优化方法,以保证抽油机的生产状态最佳,从而达到减少能耗,提高系统效率的目的。
为了实现上述目的,本发明提供一种基于功图主元分析的抽油机参数优化方法,该方法包括的步骤如下:
1)确定抽油机采油过程生产效率影响因素构成效率观测变量集合其中α1,α2为决策变量,α3~α146载荷数据环境变量,为其他环境变量,选取抽油机系统的性能变量构成性能观测变量集合:{y1,y2,y3,…yl};
2)获得所述生产效率影响因素和系统性能变量的样本数据,得到效率影响因素样本矩阵α和性能样本矩阵Y:
其中为效率影响因素个数,N为样本个数,αik表示第i个效率影响因素变量的第k个观测值,i=1,2,...,M;k=1,2,...,N;
3)利用主元分析算法对载荷数据进行降维处理,从而构建新的载荷主元变量矩阵:
4)由影响因素观测变量集合中非载荷变量与载荷新主元观测变量集合{αz1,αz2,...,αzd}构建网络输入变量集合:并令输入变量集合为:{x1,x2,x3,...,xM},即,
5)构建输入变量集合{x1,x2,x3,...,xM}观测样本值:
其中,x1~x2为决策变量,x3~xM为新的环境变量;
6)对得到的训练输入样本X、输出样本Y进行归一化处理,得到新的训练输入矩阵输出矩阵
7)构建三层前馈神经网络,其输入变量集为输出变量集为隐含层神经元个数为s1,输入层、隐含层、输出层通过权值、阈值进行连接,并且该神经网络的输入输出函数表达式为:
式中函数F(X)为S型函数;
8)利用无迹卡尔曼滤波对所述前馈神经网络进行训练,得到该神经网络的结构参数值,该训练过程包括:
①将所述神经网络中的所有权值和阈值组成状态变量I:
其中,M为输入层神经元数,s1为隐层神经元数,l为输出层神经元数,输入层至隐层神经元的连接权值为阈值为隐层至输出层的连接权值为阈值为I中的元素个数为n;设定非线性方程:
其中,函数表达式参考步骤7),为K时刻的神经网络归一化输入样本,令ωk=0,vk=0,为神经网络归一化输出样本;
②设定无迹卡尔曼计算过程中控制采样点的分布状态参数a、待选参数κ,以及非负权系数β;
③计算2n+1个采样点σ点以及σ点的相应权重,其中n为状态矩阵的I维度,λ=a2(n+κ)-n,2n+1个采样点计算如下:
每个采样点的权值计算如下:
④计算σ点的一步状态预测及状态变量协方差Pk+1|k;
⑤计算输出的一步提前预测以及协方差
⑥进行滤波更新获取新的状态矩阵、协方差矩阵、增益矩阵:
⑦对获取的新样本数据重新进行②~⑥步骤,直至所有样本对状态矩阵、协方差矩阵、增益矩阵进行了更新,从而得到适应于所有样本状态矩阵;
⑧对最后一组样本得到状态矩阵I,作为网络训练得到的权值和阈值;
⑨在得到网络参数各层权值、阈值后,确定所述前馈神经网络的函数模型为:
9)利用决策变量(x1,x2)构建多目标优化过程父代种群PD,
其中,父代种群PD中的个体的数量为K,并从x1的取值范围x1,min≤x1≤x1,max内随机取值赋予从x2的取值范围x2,min≤x2≤x2,max内随机取值赋予从而对父代种群PD进行初始化;
10)从父代种群PD中选出任意对个体,对于每对个体 进行遗传交叉计算或变异计算,并将计算结果赋予子代种群QD中相应的一对个体
11)将父代种群PD与子代种群QD进行合并得到种群R=PD∪QD,即有 将种群R的每个个体与环境变量平均值i=3,…,M合成输入样本 并计算相应的适应度函数 其中,函数h(x)=-x;
12)将种群R的所有个体所对应的适应度函数相互进行比较,将种群R的所有个体划分到具有不同层级的非支配集中,其中,对于层级较低的非支配集中的任一个体所对应的适应度函数objFun(Xs)和层级较高的非支配集中的任一个体所对应的适应度函数objFun(Xt)来说,均不存在且而对于同一层级的非支配集中的任两个个体来说,该两个不等式中至少有一个不成立;
13)按照层级从低到高的顺序从所述非支配集中选择K个个体,将选择出的K个个体的值赋予父代种群PD中的个体,并执行步骤10)-步骤13)的过程GEN次,GEN为预先确定的循环次数,最终得到优化后的L组决策变量将优化后的决策变量以及所述环境变量的平均值构成优化后的输入样本 这K个样本保证了在产液量最大,耗电量最小。
本发明的有益效果是,利用功图主元分析方法进行降维,以简化计算过程,利用无迹卡尔曼滤波神经网络(UKFNN)建立油田机采过程实施动态演化高精度模型,并利用带精英策略的快速非支配排序遗传算法对建立的模型进行搜索,探寻抽油机生产过程中最佳工艺决策参数,给出面向节能降耗的抽油机生产最佳参数指导生产,从而达到节能降耗目的。通过选择优化后的运行参数,可以使抽油机在运行过程中保证在产液量基本固定的情况下,耗电量最小,从而可以降低油田生产成本并提高油田生产效率。
附图说明
图1示出了抽油机的工作模型;
图2示出了本发明一个实施例所述的基于功图主元分析的抽油机参数优化方法的流程图;
图3示出了本发明一个实施例中的前馈神经网络的结构;
图4示出了测试样本产液量预测效果图;
图5示出了测试样本耗电量预测效果图;
图6示出了产液量偏好值与耗电量的pareto解集关系。
具体实施方式
在下面的描述中,出于说明的目的,为了提供对一个或多个实施例的全面理解,阐述了许多具体细节。然而,很明显,也可以在没有这些具体细节的情况下实现这些实施例。在其它例子中,为了便于描述一个或多个实施例,公知的结构和设备以方框图的形式示出。
图2是流程图,示出了本发明的一个实施例所述的基于功图主元分析的抽油机参数优化方法。如图2所示,本发明所述的基于功图主元分析的抽油机参数优化方法包括如下步骤:
步骤S1:确定抽油机采油过程生产效率影响因素构成效率观测变量集合其中α1,α2为决策变量,α3~α146载荷数据环境变量,为其他环境变量,选取抽油机系统的性能变量构成性能观测变量集合:{y1,y2,y3,…yl}。
在本发明的一个实施例中,选取决策变量α1为冲次、决策变量α2为有效冲程、α3~α146为载荷1至载荷144,其余环境变量包括:理论排量、功率因数、有功功率、无功功率、含水率中的一个或多个变量;选取抽油机生产过程性能变量y1为产液量、y2为耗电量。
步骤S2:获得所述生产效率影响因素和系统性能变量的样本数据,得到效率影响因素样本矩阵α和性能样本矩阵Y:
其中为效率影响因素个数,N为样本个数,αik表示第i个效率影响因素变量的第k个观测值,i=1,2,...,M;k=1,2,...,N。
设所述决策变量、环境变量和性能变量的观测值采集周期的最大值为tmax,则这些变量中的任一变量的样本取为tmax时间内该变量的观测值的平均值。
步骤S3:利用主元分析算法对载荷数据进行降维处理,从而构建新的载荷主元变量。本发明中采用示功图描绘数据的144个载荷点做为部分环境变量进行建模。然而利用144维数据建模易产生维度灾难。故而利用主元分析算法对载荷数据进行降维处理。
在一个实施例中,利用主元分析算法对载荷数据进行降维处理的步骤可以包括:
①设置样本累计贡献率precent=0.95;
②获取载荷数据每个Lk具有第k观测变量的N个观测数据,3≤k≤146;
③求出数据平均值并利用原始数据减去均值得到
④计算协方差矩阵
⑤计算协方差矩阵的特征值E1,E2,...,E144与特征向量EV1,EV2,...,EV144;
⑥由大到小依次排列特征值E′1,E'2,...,E'M,对应特征向量为EV′1,EV′2,...,EV′144,按特征值大小顺序取前d个特征值的特征向量构成矩阵[EV1',EV2',...,EVd'],此时其中d<144;特征向量代表原数据的分布方向,其对应的特征值越大,则该向量越重要(即为主元);其对应的特征值越小,则该向量越次要。
⑦由[EV′1,EV′2,...,EV′d]与原始样本求取载荷新的主元,其新载荷主元观测变量构成集合:{αz1,αz2,...,αzd},其为d个新变量,且每个变量为N个观测值构成的新主元矩阵:
步骤S4:由影响因素观测变量集合中非载荷变量与载荷新主元观测变量集合{αz1,αz2,...,αzd}构建网络输入变量集合:并令输入变量集合为:{x1,x2,x3,...,xM},即,
步骤S5:构建输入变量集合{x1,x2,x3,...,xM}观测样本值:
其中,x1~x2为决策变量,x3~xM为新的环境变量。
步骤S6:对得到的训练输入样本X、输出样本Y进行归一化处理,得到新的训练输入矩阵输出矩阵
在一个实施例中,所述归一化处理的算法如下:
其中:为设定输入变量归一化后数据范围的最大值、最小值;
xik为归一化前的第i个输入变量第k个样本值;
为归一化后第i个输入变量第k个样本值;
xi,min=min{xik|1≤k≤N}
xi,max=max{xik|1≤k≤N}
为设定输出变量归一化后数据范围的最大值、最小值;
yjk为归一化前第j个输出变量的第k个采集样本值;
为归一化后第j个输出变量的第k个值;
yj,max=max{yjk|1≤k≤N}
yj,min=min{yjk|1≤k≤N}
于是得到:
步骤S7:构建三层前馈神经网络,其输入变量集为输出变量集为隐含层神经元个数为s1,输入层、隐含层、输出层通过权值、阈值进行连接,并且该神经网络的输入输出函数表达式为:
式中函数F(X)为S型函数。图3示出了本发明的一个实施例所使用的前馈神经网络的结构。
步骤S8:利用无迹卡尔曼滤波对所述前馈神经网络进行训练,得到该神经网络的结构参数值,具体说,该训练过程包括:
①将所述神经网络中的所有权值和阈值组成状态变量I:
其中,M为输入层神经元数,s1为隐层神经元数,l为输出层神经元数,输入层至隐层神经元的连接权值为阈值为隐层至输出层的连接权值为阈值为I中的元素个数为n;设定非线性方程:
其中,函数表达式参考步骤S7,为K时刻的神经网络输入样本,令ωk=0,vk=0,为神经网络输出样本。
②设定无迹卡尔曼计算过程中控制采样点的分布状态参数a、待选参数κ,以及非负权系数β。
③计算2n+1个σ点(即采样点,一个采样点即为一组I值)以及σ点(采样点)的相应权重,其中n为状态矩阵的的I维度,λ=a2(n+κ)-n,2n+1个采样点计算如下:
每个采样点的权值计算如下:
④计算σ点的一步状态预测及状态变量协方差Pk+1|k;
⑤计算输出的一步提前预测以及协方差
⑥进行滤波更新获取新的状态矩阵、协方差矩阵、增益矩阵:
⑦对获取的新样本数据重新进行②~⑥步骤,直至所有样本对状态矩阵、协方差矩阵、增益矩阵进行了更新,从而得到适应于所有样本状态矩阵。
⑧对最后一组样本得到状态矩阵I,作为网络训练得到的权值和阈值。
⑨在得到网络参数各层权值、阈值后,确定所述前馈神经网络的函数模型为:
步骤S9:利用决策变量(x1,x2)构建父代种群PD,
其中,父代种群PD中的个体的数量为K,并从x1的取值范围x1,min≤x1≤x1,max内随机取值赋予从x2的取值范围x2,min≤x2≤x2,max内随机取值赋予从而对父代种群PD进行初始化。
步骤S10:从父代种群PD中选出任意对个体,对于每对个体 进行遗传交叉计算或变异计算,并将计算结果赋予子代种群QD中相应的一对个体
在一个实施例中,每对个体之间的遗传交叉计算公式为:
随机数α∈[0,1];
每个个体的变异计算公式为:
随机数β∈[0,1]。
步骤S11:将父代种群PD与子代种群QD进行合并得到种群R=PD∪QD,即有 将种群R的每个个体与环境变量平均值i=3,…,M合成输入样本 并计算相应的适应度函数 其中函数h(x)=-x,这是因为由于优化计算过程中,得到的产液量值越大,耗电量越小,其性能变量越好。
所述环境参数的平均值的计算公式为:
步骤S12:将种群R的所有个体所对应的适应度函数相互进行比较,将种群R的所有个体划分到具有不同层级的非支配集中,其中,对于层级较低的非支配集中的任一个体所对应的适应度函数objFun(Xs)和层级较高的非支配集中的任一个体所对应的适应度函数objFun(Xt)来说,均不存在且而对于同一层级的非支配集中的任两个个体来说,该两个不等式中至少有一个不成立。
步骤S13:按照层级从低到高的顺序从所述非支配集中选择K个个体,将选择出的K个个体的值赋予父代种群PD中的个体,返回步骤S10。
在按照层级从低到高的顺序从所述非支配集中选择K个个体时,对于同一层级非支配集中的个体,选择个体拥挤度ds较大的个体,所述个体拥挤度ds的计算方法为:
对当前种群R中所有个体所对应的适应度函数值objFun(Xs)中的按从小到大的顺序排序,另外,对所有objFun(Xs)中的按从小到大的顺序排序,令每次排序的第一个和最后一个个体的拥挤距离为无穷大,种群个体的拥挤度ds为
分别为在所述所排的序列中值的后一个值和前一个值;
分别为在所述所排的序列中值的后一个值和前一个值;
步骤S14:循环执行步骤10)-步骤13)的过程GEN次,GEN为预先确定的循环次数,最终得到优化后的K组决策变量将优化后的决策变量以及所述环境变量的平均值构成优化后的输入样本 这K个样本保证了在产液量最大时,耗电量最小。
下面以大港油田港510-3抽油机为实验对象采用本发明的方法进行优化。
确定抽油机采油过程生产效率影响因素构成效率观测变量集合其中α1,α2为决策变量冲次、有效冲程,α3~α146功图载荷数据,α147为环境变量平均功率因数,α148为环境变量平均有功功率,α149为环境变量平均无功功率,α150为环境变量理论排量、α151为环境变量计算泵效、α152为环境变量含水率。选择抽油机生产过程的产液量y1、耗电量y2作为建立抽油机系统模型的性能变量{y1,y2}。
采集油井2013年12月26日至2014年12月23日数据。对所有变量按照24小时为采集时间间隔求取24小时数据的平均值,并作为该变量样本。采集数据部分地示于表1和表2。
表1
表2
利用PCA输入变量降维处理,得到新主元数据示于表3。
表3
构建新的采集建模样本数据[X,Y]。共获得可建模数据256组,示于表4中。
表4
将匹配后的数据与载荷新主元数据进行匹配,并进行归一化处理。可建模样本归一化后部分数据示于表5中。
表5
通过无迹卡尔曼滤波对前馈神经网络的权值、阈值进行估计,将神经网络权值、阈值作为无迹卡尔曼滤波的状态变量,神经网络的输出作为无迹卡尔曼滤波的测量变量,从而得到抽油机工艺系统的精确模型。设置网络隐含层节点数为s1=5。设置隐层神经元输出函数为F(x)=1/(exp(-x)+1),得到:
输入层到隐层的权值w1:
隐层神经元阈值b1:
隐层到输出层权值w2:
输出层阈值b2:
利用上述权值和阈值,构建数学模型,选择可建模样本中最后60组数据作为测试样本进行验证。图4示出了测试样本产液量预测效果图;图5示出了测试样本耗电量预测效果图。
由模型的相对误差可知,建模效果较好,随着样本的不断训练,模型精度越来越高,符合动态建模的特性。
利用决策变量(x1,x2)的K对数据构建初始种群PD,设置种群大小50,即初始化的种群样本数量为K=50;决策变量冲次(x1)的优化范围2.0≤x1≤4.5;决策变量有效冲程(x2)的范围3.0≤x2≤3.6。设置最大遗传代数GEN=100。
计算环境参数平均值具体算法如下:
具体的环境参数平均值示于表6中。
表6
计算得到的产液量偏好值与耗电量的pareto解集关系示于图6中。
Pareto前沿所对应的数值解集示于表7中。
表7
在优化过程采用向最小值优化方法策略。故而,产液量值取反优化,向最小值方向计算,从而得到pareto前沿。
将优化后的决策变量,以及环境变量的平均值带入建立工艺过程模型,计算优化后的决策变量的系统性能,与实际样本的系统性能平均值进行比较,在保证固定的产液量的前提下,耗电量降低,则上述方法有效。去产液量偏好值解集的平均值求取实际对应的产液量大小。其对比结果示于表8。
表8
由表8可知,产液量得到提高,耗电量得到下降,系统效率得到了提高。达到了节能增产的目的。说明该方法有效。
尽管已经结合详细示出并描述的优选实施例公开了本发明,但是本领域技术人员应当理解,对于上述本发明所提出的基于功图主元分析的抽油机参数优化方法,还可以在不脱离本发明内容的基础上做出各种改进。因此,本发明的保护范围应当由所附的权利要求书的内容确定。
Claims (8)
1.一种基于功图主元分析的抽油机参数优化方法,包括如下步骤:
1)确定抽油机采油过程生产效率影响因素构成效率观测变量集合其中α1,α2为决策变量,α3~α146载荷数据环境变量,α147~αM为其他环境变量,选取抽油机系统的性能变量构成性能观测变量集合:{y1,y2,y3,…yl};
2)获得所述生产效率影响因素和系统性能变量的样本数据,得到效率影响因素样本矩阵α和性能样本矩阵Y:
其中为效率影响因素个数,N为样本个数,αik表示第i个效率影响因素变量的第k个观测值,i=1,2,...,M;k=1,2,...,N;
3)利用主元分析算法对载荷数据进行降维处理,从而构建新的载荷主元变量矩阵:
4)由影响因素观测变量集合中非载荷变量与载荷新主元观测变量集合{αz1,αz2,...,αzd}构建网络输入变量集合:并令输入变量集合为:{x1,x2,x3,...,xM},即,
5)构建输入变量集合{x1,x2,x3,...,xM}观测样本值:
其中,x1~x2为决策变量,x3~xM为新的环境变量;
6)对得到的训练输入样本X、输出样本Y进行归一化处理,得到新的训练输入矩阵输出矩阵
7)构建三层前馈神经网络,其输入变量集为输出变量集为隐含层神经元个数为s1,输入层、隐含层、输出层通过权值、阈值进行连接,并且该神经网络的输入输出函数表达式为:
式中函数F(X)为S型函数;
8)利用无迹卡尔曼滤波对所述前馈神经网络进行训练,得到该神经网络的结构参数值,该训练过程包括:
①将所述神经网络中的所有权值和阈值组成状态变量I:
其中,M为输入层神经元数,s1为隐层神经元数,l为输出层神经元数,输入层至隐层神经元的连接权值为 阈值为 隐层至输出层的连接权值为 阈值为 I中的元素个数为n;设定非线性方程:
其中,函数表达式参考步骤7),为K时刻的神经网络输入样本,令ωk=0,vk=0,为神经网络输出样本;
②设定无迹卡尔曼计算过程中控制采样点的分布状态参数a、待选参数κ,以及非负权系数β;
③计算2n+1个采样点(σ点)以及σ点的相应权重,其中n为状态矩阵的I维度,λ=a2(n+κ)-n,2n+1个采样点计算如下:
每个采样点的权值计算如下:
④计算σ点的一步状态预测及状态变量协方差Pk+1|k;
⑤计算输出的一步提前预测以及协方差
⑥进行滤波更新获取新的状态矩阵、协方差矩阵、增益矩阵:
⑦对获取的新样本数据重新进行②~⑥步骤,直至所有样本对状态矩阵、协方差矩阵、增益矩阵进行了更新,从而得到适应于所有样本状态矩阵;
⑧对最后一组样本得到状态矩阵I,作为网络训练得到的权值和阈值;
⑨在得到网络参数各层权值、阈值后,确定所述前馈神经网络的函数模型为:
9)利用决策变量(x1,x2)构建多目标优化父代种群PD
其中,父代种群PD中的个体的数量为K,并从x1的取值范围x1,min≤x1≤x1,max内随机取值赋予从x2的取值范围x2,min≤x2≤x2,max内随机取值赋予从而对父代种群PD进行初始化;
10)从父代种群PD中选出任意对个体,对于每对个体 进行遗传交叉计算或变异计算,并将计算结果赋予子代种群QD中相应的一对个体
11)将父代种群PD与子代种群QD进行合并得到种群R=PD∪QD,即有 将种群R的每个个体与环境变量平均值i=3,…,M合成输入样本 并计算相应的适应度函数 其中,函数h(x)=-x;
12)将种群R的所有个体所对应的适应度函数相互进行比较,将种群R的所有个体划分到具有不同层级的非支配集中,其中,对于层级较低的非支配集中的任一个体所对应的适应度函数objFun(Xs)和层级较高的非支配集中的任一个体所对应的适应度函数objFun(Xt)来说,均不存在且而对于同一层级的非支配集中的任两个个体来说,该两个不等式中至少有一个不成立;
13)按照层级从低到高的顺序从所述非支配集中选择K个个体,将选择出的K个个体的值赋予父代种群PD中的个体,并执行步骤10)-步骤13)的过程GEN次,GEN为预先确定的循环次数,最终得到优化后的K组决策变量将优化后的决策变量以及所述环境变量的平均值构成优化后的输入样本 这K个样本保证了在产液量最大,耗电量最小。
2.如权利要求1所述的基于功图主元分析的抽油机参数优化方法,其中,
所述决策变量α1为冲次、决策变量α2为有效冲程、α3~α146为载荷1至载荷144,其余环境变量包括:理论排量、功率因数、有功功率、无功功率、含水率中的一个或多个变量;所述抽油机生产过程性能变量y1为产液量、y2为耗电量。
3.如权利要求1所述的基于功图主元分析的抽油机参数优化方法,其中,
设所述决策变量、环境变量和性能变量的观测值采集周期的最大值为tmax,则这些变量中的任一变量的样本取为tmax时间内该变量的观测值的平均值。
4.如权利要求1所述的基于功图主元分析的抽油机参数优化方法,其中,
利用主元分析算法对载荷数据进行降维处理的步骤包括:
①设置样本累计贡献率precent=0.95;
②获取载荷数据每个Lk具有第k观测变量的N个观测数据,3≤k≤146;
③求出数据平均值并利用原始数据减去均值得到
④计算协方差矩阵
⑤计算协方差矩阵的特征值E1,E2,...,E144与特征向量EV1,EV2,...,EV144;
⑥由大到小依次排列特征值E′1,E'2,...,E'M,对应特征向量为EV′1,EV′2,...,EV′144,按特征值大小顺序取前d个特征值的特征向量构成矩阵[EV1',EV2',...,EVd'],此时 其中d<144;
⑦由[EV′1,EV′2,...,EV′d]与原始样本求取载荷新的主元,其新载荷主元观测变量构成集合:{αz1,αz2,...,αzd},其为d个新变量,且每个变量为N个观测值构成的新主元矩阵:
5.如权利要求1所述的基于功图主元分析的抽油机参数优化方法,其中,
所述归一化处理的算法如下:
i=1,2,...,M;k=1,2,...,N
j=1,2,....,l;k=1,2,...,N
其中:为设定输入变量归一化后数据范围的最大值、最小值;
xik为归一化前的第i个输入变量第k个样本值;
为归一化后第i个输入变量第k个样本值;
xi,min=min{xik|1≤k≤N}
xi,max=max{xik|1≤k≤N}
为设定输出变量归一化后数据范围的最大值、最小值;
yjk为归一化前第j个输出变量的第k个采集样本值;
为归一化后第j个输出变量的第k个值;
yj,max=max{yjk|1≤k≤N}
yj,min=min{yjk|1≤k≤N}
于是得到:
6.如权利要求1所述的基于功图主元分析的抽油机参数优化方法,其中,
多目标优化的每对个体之间的遗传交叉计算公式为:
随机数α∈[0,1];
每个个体的变异计算公式为:
随机数β∈[0,1]。
7.如权利要求1所述的基于功图主元分析的抽油机参数优化方法,其中,
所述环境参数的平均值的计算公式为:
8.如权利要求1所述的基于功图主元分析的抽油机参数优化方法,其中,
在按照层级从低到高的顺序从所述非支配集中选择K个个体时,对于同一层级非支配集中的个体,选择个体拥挤度ds较大的个体,所述个体拥挤度ds的计算方法为:
对当前种群R中所有个体所对应的适应度函数值objFun(Xs)中的按从小到大的顺序排序,另外,对所有objFun(Xs)中的按从小到大的顺序排序,令每次排序的第一个和最后一个个体的拥挤距离为无穷大,种群个体 的拥挤度ds为
分别为在所述所排的序列中值的后一个值和前一个值;
分别为在所述所排的序列中值的后一个值和前一个值;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510111157.6A CN105046326B (zh) | 2015-03-13 | 2015-03-13 | 基于功图主元分析的抽油机参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510111157.6A CN105046326B (zh) | 2015-03-13 | 2015-03-13 | 基于功图主元分析的抽油机参数优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105046326A true CN105046326A (zh) | 2015-11-11 |
CN105046326B CN105046326B (zh) | 2019-04-19 |
Family
ID=54452855
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510111157.6A Active CN105046326B (zh) | 2015-03-13 | 2015-03-13 | 基于功图主元分析的抽油机参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105046326B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106502096A (zh) * | 2016-11-14 | 2017-03-15 | 重庆科技学院 | 基于偏好多目标优化的油田机采过程决策参数优化方法 |
CN106529042A (zh) * | 2016-11-14 | 2017-03-22 | 重庆科技学院 | 基于计算智能的油田机采参数动态演化建模与优化方法 |
CN106773667A (zh) * | 2016-11-14 | 2017-05-31 | 重庆科技学院 | 基于无迹粒子滤波神经网络的油田机采参数建模方法 |
CN106777466A (zh) * | 2016-11-14 | 2017-05-31 | 重庆科技学院 | 基于st‑upfnn算法的高含硫天然气净化工艺的动态演化建模方法 |
CN106777465A (zh) * | 2016-11-14 | 2017-05-31 | 重庆科技学院 | 高含硫天然气净化工艺动态演化建模与节能优化方法 |
CN107229970B (zh) * | 2017-06-27 | 2018-08-17 | 重庆科技学院 | 共享直饮水水质自适应动态自学习在线监测系统 |
CN108678728A (zh) * | 2018-05-15 | 2018-10-19 | 北京理工大学 | 一种基于k-means的油井参数分析组合算法 |
CN109935304A (zh) * | 2019-03-15 | 2019-06-25 | 重庆科技学院 | 一种基于用户体验测试系统的自闭症拥抱机智能设计建模与决策参数优化方法 |
CN113153267A (zh) * | 2021-03-11 | 2021-07-23 | 中国石油大学(北京) | 一种抽油机井示功图多尺度归一化方法及系统 |
CN116756513A (zh) * | 2023-07-05 | 2023-09-15 | 大庆石油管理局有限公司 | 一种塔架式抽油机数据滤波算法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7295690B2 (en) * | 2003-05-27 | 2007-11-13 | Canon Kabushiki Kaisha | Image processing method and apparatus |
CN102982373A (zh) * | 2012-12-31 | 2013-03-20 | 山东大学 | 一种混合svm回归算法的oin神经网络训练方法 |
CN103061749A (zh) * | 2013-01-25 | 2013-04-24 | 东北大学 | 一种有杆泵抽油井井下动态液位软测量方法 |
CN103544389A (zh) * | 2013-10-18 | 2014-01-29 | 丽水学院 | 基于故障树和模糊神经网络的汽车起重机故障诊断方法 |
-
2015
- 2015-03-13 CN CN201510111157.6A patent/CN105046326B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7295690B2 (en) * | 2003-05-27 | 2007-11-13 | Canon Kabushiki Kaisha | Image processing method and apparatus |
CN102982373A (zh) * | 2012-12-31 | 2013-03-20 | 山东大学 | 一种混合svm回归算法的oin神经网络训练方法 |
CN103061749A (zh) * | 2013-01-25 | 2013-04-24 | 东北大学 | 一种有杆泵抽油井井下动态液位软测量方法 |
CN103544389A (zh) * | 2013-10-18 | 2014-01-29 | 丽水学院 | 基于故障树和模糊神经网络的汽车起重机故障诊断方法 |
Non-Patent Citations (1)
Title |
---|
李宾元: "抽油机四连杆机构参数的多目标函数优化的探讨", 《西南石油学院学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529042B (zh) * | 2016-11-14 | 2019-12-31 | 重庆科技学院 | 基于计算智能的油田机采参数动态演化建模与优化方法 |
CN106777466B (zh) * | 2016-11-14 | 2020-06-05 | 重庆科技学院 | 基于st-upfnn算法的高含硫天然气净化工艺的动态演化建模方法 |
CN106773667A (zh) * | 2016-11-14 | 2017-05-31 | 重庆科技学院 | 基于无迹粒子滤波神经网络的油田机采参数建模方法 |
CN106777466A (zh) * | 2016-11-14 | 2017-05-31 | 重庆科技学院 | 基于st‑upfnn算法的高含硫天然气净化工艺的动态演化建模方法 |
CN106777465A (zh) * | 2016-11-14 | 2017-05-31 | 重庆科技学院 | 高含硫天然气净化工艺动态演化建模与节能优化方法 |
CN106777465B (zh) * | 2016-11-14 | 2020-06-09 | 重庆科技学院 | 高含硫天然气净化工艺动态演化建模与节能优化方法 |
CN106529042A (zh) * | 2016-11-14 | 2017-03-22 | 重庆科技学院 | 基于计算智能的油田机采参数动态演化建模与优化方法 |
CN106502096B (zh) * | 2016-11-14 | 2019-06-14 | 重庆科技学院 | 基于偏好多目标优化的油田机采过程决策参数优化方法 |
CN106502096A (zh) * | 2016-11-14 | 2017-03-15 | 重庆科技学院 | 基于偏好多目标优化的油田机采过程决策参数优化方法 |
CN107229970B (zh) * | 2017-06-27 | 2018-08-17 | 重庆科技学院 | 共享直饮水水质自适应动态自学习在线监测系统 |
CN108678728A (zh) * | 2018-05-15 | 2018-10-19 | 北京理工大学 | 一种基于k-means的油井参数分析组合算法 |
CN109935304A (zh) * | 2019-03-15 | 2019-06-25 | 重庆科技学院 | 一种基于用户体验测试系统的自闭症拥抱机智能设计建模与决策参数优化方法 |
CN113153267A (zh) * | 2021-03-11 | 2021-07-23 | 中国石油大学(北京) | 一种抽油机井示功图多尺度归一化方法及系统 |
CN113153267B (zh) * | 2021-03-11 | 2023-04-18 | 中国石油大学(北京) | 一种抽油机井示功图多尺度归一化方法及系统 |
CN116756513A (zh) * | 2023-07-05 | 2023-09-15 | 大庆石油管理局有限公司 | 一种塔架式抽油机数据滤波算法 |
Also Published As
Publication number | Publication date |
---|---|
CN105046326B (zh) | 2019-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105045941B (zh) | 基于无迹卡尔曼滤波的抽油机参数优化方法 | |
CN105046326A (zh) | 基于功图主元分析的抽油机参数优化方法 | |
CN109492822B (zh) | 空气污染物浓度时空域关联预测方法 | |
CN102622418B (zh) | 一种基于bp神经网络的预测装置及设备 | |
CN103235877B (zh) | 机器人控制软件模块划分方法 | |
CN110348713A (zh) | 一种基于关联性分析和数据挖掘的台区线损计算方法 | |
CN107657243A (zh) | 基于遗传算法优化的神经网络雷达一维距离像目标识别方法 | |
CN106055791B (zh) | 基于预估校正算法的飞行器全局气动优化方法 | |
CN104523264B (zh) | 一种心电信号处理方法 | |
CN106529166A (zh) | 一种基于maepso算法的区域水资源优化配置方法 | |
CN111985610A (zh) | 一种基于时序数据的抽油机井泵效预测系统和方法 | |
CN105184416A (zh) | 基于粒子群优化反向传播神经网络的脉动风速预测方法 | |
CN104680025A (zh) | 基于遗传算法极限学习机的抽油机参数优化方法 | |
CN106355003B (zh) | 基于t分布的马尔科夫链蒙特卡洛自动历史拟合方法及系统 | |
CN103226741A (zh) | 城市供水管网爆管预测方法 | |
CN106778846A (zh) | 一种基于支持向量机的风电场风速预测方法 | |
CN108445752A (zh) | 一种自适应选择深度特征的随机权神经网络集成建模方法 | |
CN105760673A (zh) | 一种河流相储层地震敏感参数模板分析方法 | |
CN102622515A (zh) | 一种天气预测方法 | |
CN104680023B (zh) | 基于多目标决策的抽油机参数优化方法 | |
CN109492748A (zh) | 一种基于卷积神经网络的电力系统的中长期负荷预测模型建立方法 | |
CN104778338A (zh) | 一种低能耗分子蒸馏过程给定值的优化方法 | |
CN104050547A (zh) | 一种油田开发规划方案非线性优选决策方法 | |
CN103793887A (zh) | 基于自适应增强算法的短期电力负荷在线预测方法 | |
CN105631554A (zh) | 一种基于时间序列的油井油液含水率多模型预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201204 Address after: 151400 west of Anton 1st Avenue and south of Anton 10th Road, Hadaqi industrial corridor, Anda City, Suihua City, Heilongjiang Province Patentee after: HEILONGJIANG RONGZE PETROLEUM EQUIPMENT Co.,Ltd. Address before: 401331 Chongqing city Shapingba District hogye University City, Chongqing University of Science and Technology Patentee before: Chongqing University of Science & Technology |