CN104680257B - 面向抽油机工作过程的智能优化方法 - Google Patents
面向抽油机工作过程的智能优化方法 Download PDFInfo
- Publication number
- CN104680257B CN104680257B CN201510108394.7A CN201510108394A CN104680257B CN 104680257 B CN104680257 B CN 104680257B CN 201510108394 A CN201510108394 A CN 201510108394A CN 104680257 B CN104680257 B CN 104680257B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- mover
- mtr
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- External Artificial Organs (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供一种面向抽油机工作过程的智能优化方法,包括:1)确定抽油机的生产效率影响因素和性能变量、2)获得其样本数据、3)对载荷数据降维、4)由非载荷变量与载荷新主元构建网络输入变量集、5)构建输入变量的样本值、6)归一化处理、7)选取训练和测试样本、8)构建前馈神经网络并训练、9)构造产液量偏好函数、10)构造父代和子代种群、11)对父代种群个体进行遗传交叉或变异计算,以产生子代种群个体、12)求个体求适应度函数、13)将个体划分到不同层级的非支配集中、14)从非支配集中选择个体构成新父代种群并循环,最终得到优化后的生产效率影响因素值。优化后,可保证在产液量基本固定的情况下,耗电量最小。
Description
技术领域
本发明属于采油领域,具体涉及一种面向抽油机工作过程的智能优化方法。
背景技术
抽油机采油作为一种机械采油方式,主要由电动机、地面传动设备和井下抽油设备三部分组成,如图1所示。抽油机的整个采油过程主要分为上下两个冲程:上冲程,即驴头悬点向上运动,提起抽油杆柱和井下抽油设备,此过程中电动机需消耗大量的能量;下冲程,即驴头悬点向下运动,抽油机的抽油杆柱电动机做功。在抽油杆柱上下运动过程中,电动机的负载发生周期变化。抽油机的运行参数的选择对整个抽油机系统的能量消耗影响很大。为了使抽油机采油生产过程既能完成预定的产液量,又能使抽油机生产过程的耗电量最低,需要对抽油机运行参数进行节能优化。
发明内容
本发明是为了解决现有技术中存在的上述技术问题而做出,其目的在于提供一种面向抽油机工作过程的智能优化方法,以保证抽油机的生产状态最佳,从而达到减少能耗,提高系统效率的目的。
为了实现上述目的,本发明提供一种面向抽油机工作过程的智能优化方法,该方法包括的步骤如下:
1)确定抽油机采油过程生产效率影响因素构成效率观测变量集合{α1,α2,α3,…αM},其中α1,α2为决策变量,α3~α146载荷数据环境变量,α147~αM为其他环境变量,选取抽油机系统的性能变量构成性能观测变量集合:{y1,y2,y3,…yl};
2)获得所述生产效率影响因素和系统性能变量的样本数据,得到效率影响因素样本矩阵α和性能样本矩阵Y:
其中M为效率影响因素个数,N为样本个数,αik表示第i个效率影响因素变量的第k个观测值,i=1,2,...,M;k=1,2,...,N;
3)利用主元分析算法对载荷数据进行降维处理,从而构建新的载荷主元变量矩阵:
4)由影响因素效率观测变量集合{α1,α2,α3,…αM}中非载荷变量与载荷新主元观测变量集合{αz1,αz2,...,αzd}构建网络输入变量集合:{α1,α2,α147,…,αM,αz1,…,αzd},并令输入变量集合为:{x1,x2,x3,...,xM},即,{α1,α2,α147,…,αM,αz1,…,αzd}={x1,x2,x3,...,xM};
5)构建输入变量集合{x1,x2,x3,...,xM}观测样本值:
其中,x1~x2为决策变量,x3~xM为新的环境变量;
6)对得到的训练输入样本X、输出样本Y进行归一化处理,得到新的训练输入矩阵输出矩阵
7)在归一化后样本集中,选取样本集前组样本作为网络训练样本集则样本集剩余组作为测试样本集其中其中:训练样本集和测试样本集分别为:
8)构建3层前馈神经网络,采用所述训练输入样本对该网络进行训练,利用遗传算法训练网络权值W、权值β、隐含层阈值b,从而确定网络的参数,最终得到可计算出任一组输入观测值x1,x2,...,xM所对应的输出预测值的前馈神经网络:
9)针对产液量y1构造其偏好函数h=h(y1),该偏好函数为U形曲线,在其整个定义域上二阶可导且二阶导数恒大于零,将产液量y1的值划分为好、较好、一般、较差和极差5个区域,并通过所述偏好函数将该五个区域量化为数值h=h(y1);其中,产液量y1处于好区域对应偏好函数最小值;
10)利用决策变量(x1,x2)构建多目标优化父代种群PD,
其中,父代种群PD中的个体的数量为K,并从x1的取值范围x1,min≤x1≤x1,max内随机取值赋予从x2的取值范围x2,min≤x2≤x2,max内随机取值赋予从而对父代种群PD进行初始化;
11)从父代种群PD中选出任意对个体,对于每对个体 进行遗传交叉计算或变异计算,并将计算结果赋予子代种群QD中相应的一对个体
12)将父代种群PD与子代种群QD进行合并得到种群R=PD∪QD,即有将种群R的每个个体与环境变量平均值i=3,…,M合成输入样本并计算相应的适应度函数
13)将种群R的所有个体所对应的适应度函数相互进行比较,将种群R的所有个体划分到具有不同层级的非支配集中,其中,对于层级较低的非支配集中的任一个体所对应的适应度函数objFun1(Xs)和层级较高的非支配集中的任一个体所对应的适应度函数objFun1(Xt)来说,均不存在且而对于同一层级的非支配集中的任两个个体来说,该两个不等式中至少有一个不成立;
14)按照层级从低到高的顺序从所述非支配集中选择K个个体,将选择出的K个个体的值赋予父代种群PD中的个体,并执行步骤11)-步骤14)的过程GEN次,GEN为预先确定的循环次数,最终得到优化后的K组决策变量将优化后的决策变量以及所述环境变量的平均值构成优化后的输入样本这K个样本保证了在产液量基本固定且具有最小偏好值的情况下,耗电量最小。
本发明的有益效果是,利用遗传算法优化的极限学习机算法(GAELM)建立油田机采过程的高精度模型,并利用带精英策略的快速非支配排序遗传算法(NSGA2)对建立的模型进行搜索,探寻抽油机生产过程中最佳工艺决策参数,给出面向节能降耗的抽油机生产最佳参数,指导生产。通过选择优化后的运行参数,可以使抽油机在运行过程中保证在产液量基本固定且具有最小偏好值的情况下,耗电量最小,从而可以降低油田生产成本并提高油田生产效率。
附图说明
图1示出了抽油机的工作模型;
图2示出了本发明一个实施例所述的面向抽油机工作过程的的智能优化方法的流程图;
图3示出了本发明一个实施例中的前馈神经网络的结构;
图4示出了本发明一个实施例中的产液量偏好函数的图形;
图5示出了利用本发明的方法训练神经网络所产生的产液量训练效果图;
图6示出了利用本发明的方法训练神经网络所产生的耗电量训练效果图;
图7示出了测试样本产液量预测效果图;
图8示出了测试样本耗电量预测效果图;
图9示出了偏好函数的图形;
图10示出了产液量偏好值与耗电量的pareto解集关系。
具体实施方式
在下面的描述中,出于说明的目的,为了提供对一个或多个实施例的全面理解,阐述了许多具体细节。然而,很明显,也可以在没有这些具体细节的情况下实现这些实施例。在其它例子中,为了便于描述一个或多个实施例,公知的结构和设备以方框图的形式示出。
图2是流程图,示出了本发明的一个实施例所述的面向抽油机工作过程的智能优化方法。如图2所示,本发明所述的面向抽油机工作过程的智能优化方法包括如下步骤:
步骤S1:确定抽油机采油过程生产效率影响因素构成效率观测变量集合{α1,α2,α3,…αM},其中α1,α2为决策变量,α3~α146载荷数据环境变量,α147~αM为其他环境变量,选取抽油机系统的性能变量构成性能观测变量集合:{y1,y2,y3,…yl}。
在本发明的一个实施例中,选取决策变量α1为冲次、决策变量α2为有效冲程、α3~α146为载荷1至载荷144,其余环境变量包括:理论排量、功率因数、有功功率、无功功率、含水率中的一个或多个变量;选取抽油机生产过程性能变量y1为产液量、y2为耗电量。
步骤S2:获得所述生产效率影响因素和系统性能变量的样本数据,得到效率影响因素样本矩阵α和性能样本矩阵Y:
其中M为效率影响因素个数,N为样本个数,αik表示第i个效率影响因素变量的第k个观测值,i=1,2,...,M;k=1,2,...,N。
步骤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个特征值的特征向量构成矩阵[EV′1,EV′2,...,EV′d],此时其中d<144;特征向量代表原数据的分布方向,其对应的特征值越大,则该向量越重要(即为主元);其对应的特征值越小,则该向量越次要。
⑦由[EV′1,EV′2,...,EV′d]与原始样本求取载荷新的主元,其新载荷主元观测变量构成集合:{αz1,αz2,...,αzd},其为d个新变量,且每个变量为N个观测值构成的新主元矩阵:
步骤S4:由影响因素效率观测变量集合{α1,α2,α3,…αM}中非载荷变量与载荷新主元观测变量集合{αz1,αz2,...,αzd}构建网络输入变量集合:{α1,α2,α147,…,αM,αz1,…,αzd},并令输入变量集合为:{x1,x2,x3,...,xM},即,{α1,α2,α147,…,αM,αz1,…,α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:在归一化后样本集中,选取样本集前组样本作为网络训练样本集则样本集剩余组作为测试样本集其中其中:训练样本集和测试样本集分别为:
步骤S8:构建3层前馈神经网络,采用所述训练输入样本对该网络进行训练,利用遗传算法训练网络权值W、权值β、隐含层阈值b,从而确定网络的参数,最终得到可计算出任一组输入观测值x1,x2,...,xM所对应的输出预测值的前馈神经网络:
具体地,在一个实施例中,如图3所示,所述3层前馈神经网络可以由输入层、隐含层和输出层组成,输入层与隐含层由权值W连接,其Wki表示第i个输入神经元与隐含层第k个神经元的间的连接权值;隐含层与输出层由权值β连接,其βkj表示第k个隐层神经元与输出层第j个神经元的间的连接权值;隐含层的神经元个数为s1,隐含层神经元的阈值为b,bk为第k个隐层神经元阈值;
各权值、阈值具体为:
设隐含层神经元的激活函数为网络的第m组输入样本的输出为:
其中,Wk=[Wk1,Wk2,…,WkM],最终得到可计算出任一组输入观测值x1,x2,...,xM所对应的输出预测值的前馈神经网络:
在步骤S8中,利用遗传算法训练网络权值W、权值β、隐含层阈值b的步骤可以包括:
①利用输入层到隐含层的权值W,隐含层阈值b构建初始种群P,该种群中第i个个体为设置初始种群P个体K=50,设置最大遗传代数GEN=100。
②构建遗传算法优化权值阈值计算的适应值函数objfun2(Pi):
③计算第i个个体的适应度函数值,利用个体和训练样本输入样本计算隐含层神经元的输出矩阵H(i):
由训练样本输入样本对应的输出样本与隐含层的输出矩阵H(i)共同求出隐含层到输出层的权值β(i),(H(i))-1为矩阵H(i)的逆矩阵;在求出个体对应的β(i)后,构建网络预测函数F(i)(Pi,β(i),X);并将测试样本集Xtest中的输入矩阵带入F(i)(Pi,β(i),X)得到预测测试集预测输出故而计算个体Pi适应度函数值:
④对种群P进行选择、交叉、变异操作,得到子群Q。
⑤种群Q进行③、④操作得到种群Q的子群。循环进行以上步骤,得到第GEN代子群,并在第GEN代子群中选择一个最佳个体Pbest,个体Pbest第GEN代子群中计算得到测试样本预测输出值与实际输出值的均方差最小。
⑥由Pbest确定的输入到隐层的权值W,隐层阈值b,计算隐层输出矩阵Hbest。计算计算隐层到输出的权值矩阵输出权值βbest:从而得到最佳模型。
上述选择操作可以包括:在P中K个个体中,按每个个体的概率φ(Pi)大小决定选中的机会,分K次从P中选择K个个体,其中概率
上述交叉操作包括:将选择操作得出的K个个体进行二进制编码,个体中每个变量占据4字节,这样编码后一个个体的字节数为:4*(s1*M+s1);将个体Pi和Pj二进制编码后,将两个个体二进制编码第4·c位以后的编码进行调换,从而得到两个新的个体和
所述变异操作包括:将交叉操作得到的K个个体进行二进制编码,其中某个二进制编码后的个体若第c位为1,则将1变成0;若第c位为0,则将0变成1,,从而得到一个新的个体Qi。
步骤S9:针对产液量y1构造其偏好函数h=h(y1),该偏好函数为U形曲线,在其整个定义域上二阶可导且二阶导数恒大于零,将产液量y1的值划分为好、较好、一般、较差和极差5个区域,并通过所述偏好函数将该五个区域量化为数值h=h(y1);其中,产液量y1处于好区域对应偏好函数最小值。图4示出了本发明的一个实施例所使用的产液量偏好函数的图形。
步骤S10:利用决策变量(x1,x2)构建多目标优化的父代种群PD:
其中,父代种群PD中的个体的数量为K,并从x1的取值范围x1,min≤x1≤x1,max内随机取值赋予从x2的取值范围x2,min≤x2≤x2,max内随机取值赋予从而对父代种群PD进行初始化。
步骤S11:从父代种群PD中选出任意对个体,对于每对个体 进行遗传交叉计算或变异计算,并将计算结果赋予多目标优化的子代种群QD中相应的一对个体
在一个实施例中,每对个体 之间的遗传交叉计算公式为:
随机数α∈[0,1];
每个个体的变异计算公式为:
随机数β∈[0,1]。
步骤S12:将父代种群PD与子代种群QD进行合并得到种群R=PD∪QD,即有将种群R的每个个体与环境变量平均值i=3,…,M合成输入样本并计算相应多目标优化计算适应度函数
所述环境参数的平均值的计算公式为:
步骤S13:将种群R的所有个体所对应的适应度函数相互进行比较,将种群R的所有个体划分到具有不同层级的非支配集中,其中,对于层级较低的非支配集中的任一个体所对应的适应度函数objFun1(Xs)和层级较高的非支配集中的任一个体所对应的适应度函数objFun1(Xt)来说,均不存在且而对于同一层级的非支配集中的任两个个体来说,该两个不等式中至少有一个不成立;
步骤S14:按照层级从低到高的顺序从所述非支配集中选择K个个体,将选择出的K个个体的值赋予父代种群PD中的个体,返回步骤S11。
在按照层级从低到高的顺序从所述非支配集中选择K个个体时,对于同一层级非支配集中的个体,选择个体拥挤度ds较大的个体,所述个体拥挤度ds的计算方法为:
对当前种群R中所有个体所对应的适应度函数值objFun1(Xs)中的按从小到大的顺序排序,另外,对所有objFun1(Xs)中的按从小到大的顺序排序,令每次排序的第一个和最后一个个体的拥挤距离为无穷大,种群个体的拥挤度ds为
分别为在所述所排的序列中值的后一个值和前一个值;
分别为在所述所排的序列中值的后一个值和前一个值;
步骤S15:循环执行步骤S11-步骤S14的过程GEN次,GEN为预先确定的循环次数,最终得到优化后的K组决策变量将优化后的决策变量以及所述环境变量的平均值构成优化后的输入样本这K个样本保证了在产液量基本固定且具有最小偏好值的情况下,耗电量最小。
下面以大港油田港510-3抽油机为实验对象采用本发明的方法进行智能优化。
确定抽油机采油过程生产效率影响因素构成效率观测变量集合{α1,α2,α3,…αM}。其中α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=20,选择可建模样本中的200组作为训练样本。基于遗传算法得到最佳输入层到隐含层权值W和阈值b如下:
图5示出了利用本发明的方法训练神经网络所产生的产液量训练效果图;图6示出了利用本发明的方法训练神经网络所产生的耗电量训练效果图。图7示出了测试样本产液量预测效果图;图8示出了测试样本耗电量预测效果图。
由建模效果可知,训练效果较优,产液量、耗电量训练误差均控制在较小的范围,训练效果较好。由所得的建模与测试样本进行预测,得到测试样本误差百分比也控制较小的范围内,故知模型泛化能力较好。
产液量偏好函数设计,根据偏好函数设计准则,采用3-S类软偏好设计,根据生产者对产液量的多少的喜爱程度用:极好、好、一般、较差、差等五个个等级描述。极好程度利用建模训练集中产液量所有值的平均值作为极好程度的偏好值为:47.3,并利用训练集中产液量最小值37.3、最大值57.3作为不可接受域的临界值。故而设计偏好程度区间为:[0,2],[2,4],[4,6],[6,8],[8,10]等,且设计的偏好程度区间所对应的实际产液量区间为:[37.3,38.4],[38.4,39.6],[39.6,40.9],[40.9,42.8],[42.8,47.3]以及[47.3,51.8],[51.8,53.6],[53.6,55.5],[55.5,56.2],[56.2.57.3]。故而设计产液量偏好函数为h(y1)=0.1×(y1-47.3)2。图9示出了偏好函数的图形;
利用决策变量(x1,x2)的K对数据构建初始种群PD,设置种群大小50,即初始化的种群样本数量为K=50;决策变量冲次(x1)的优化范围2.0≤x1≤4.5;决策变量有效冲程(x2)的范围3.0≤x2≤3.6。设置最大遗传代数GEN=100。
算环境参数平均值具体算法如下:
具体的环境参数平均值示于表6中。
表6
计算得到的产液量偏好值与耗电量的pareto解集关系示于图10中。Pareto前沿所对应的数值解集示于表7中。
表7
在表7中,产液量偏好程度属于[0,1.3]区间,产液量在极好的范围内优化,符合抽油机定量采油的工作制度。
将优化后的决策变量,以及环境变量的平均值带入建立工艺过程模型,计算优化后的决策变量的系统性能,与实际样本的系统性能平均值进行比较,在保证固定的产液量的前提下,耗电量降低,则上述方法有效。去产液量偏好值解集的平均值求取实际对应的产液量大小。其对比结果示于表8。
表8
由表8可知,计算得到偏好程度平均值为0.52,此时产液量可能为49.58t/d或者45.02t/d。此时计算得到在保持定量式采油工作状态下,耗电量下降16.75%。系统效率比(产液量/耗电量)可提高8%。达到了节能设计目的。说明该方法有效。
尽管已经结合详细示出并描述的优选实施例公开了本发明,但是本领域技术人员应当理解,对于上述本发明所提出的面向抽油机工作过程的智能优化方法,还可以在不脱离本发明内容的基础上做出各种改进。因此,本发明的保护范围应当由所附的权利要求书的内容确定。
Claims (10)
1.一种面向抽油机工作过程的智能优化方法,包括如下步骤:
1)确定抽油机采油过程生产效率影响因素构成效率观测变量集合{α1,α2,α3,…αM},其中α1,α2为决策变量,α3~α146载荷数据环境变量,α147~αM为其他环境变量,选取抽油机系统的性能变量构成性能观测变量集合:{y1,y2,y3,…yl};
2)获得所述生产效率影响因素和系统性能变量的样本数据,得到效率影响因素样本矩阵α和性能样本矩阵Y:
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>M</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>M</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&equiv;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mi>M</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>Y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mi>l</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mi>l</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mi>l</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中M为效率影响因素个数,N为样本个数,αik表示第i个效率影响因素变量的第k个观测值,i=1,2,...,M;k=1,2,...,N;yjk为归一化前第j个输出变量的第k个采集样本值,j=1,2,....,l;
3)利用主元分析算法对载荷数据进行降维处理,从而构建新的载荷主元变量矩阵:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>z</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>z</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>z</mi>
<mn>3</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>11</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>12</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>21</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>22</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
载荷新主元观测变量构成集合:{αz1,αz2,...,αzd},其为d个新变量,且每个变量为N个观测值构成的新主元矩阵;
4)由影响因素效率观测变量集合{α1,α2,α3,…αM}中非载荷变量与载荷新主元观测变量集合{αz1,αz2,...,αzd}构建网络输入变量集合:{α1,α2,α147,…,αM,αz1,…,αzd},并令输入变量集合为:{x1,x2,x3,...,xM},即,{α1,α2,α147,…,αM,αz1,…,αzd}={x1,x2,x3,...,xM};
5)构建输入变量集合{x1,x2,x3,...,xM}观测样本值:
<mrow>
<mi>X</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>X</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>X</mi>
<mi>N</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mn>147</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mn>147</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mn>147</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>M</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>M</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>11</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>12</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>31</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>32</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>M</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>M</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>Y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Y</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mi>N</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mi>l</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mi>l</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mrow>
<mi>l</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,x1~x2为决策变量,x3~xM为新的环境变量;
6)对得到的训练输入样本X、输出样本Y进行归一化处理,得到新的训练输入矩阵输出矩阵
<mrow>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>3</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>M</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>&equiv;</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>l</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>&equiv;</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>N</mi>
</msub>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mi>N</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
7)在归一化后样本集中,选取样本集前组样本作为网络训练样本集
则样本集剩余组作为测试样本集其中其中:训练样本集和测试样本集分别为:
<mrow>
<msup>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<msup>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msup>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<msup>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
8)构建3层前馈神经网络,采用所述训练输入样本对该网络进行训练,利用遗传算法训练网络权值W、权值β、隐含层阈值b,从而确定网络的参数,最终得到可计算出任一组输入观测值x1,x2,...,xM所对应的输出预测值的前馈神经网络:
<mrow>
<mover>
<mi>Y</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>g</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>g</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msup>
<mi>g</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mi>l</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>;</mo>
</mrow>
9)针对产液量y1构造其偏好函数h=h(y1),该偏好函数为U形曲线,在其整个定义域上二阶可导且二阶导数恒大于零,将产液量y1的值划分为好、较好、一般、较差和极差5个区域,并通过所述偏好函数将该五个区域量化为数值h=h(y1);其中,产液量y1处于好区域对应偏好函数最小值;
10)利用决策变量(x1,x2)构建多目标优化父代种群PD,
<mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<mo>{</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>m</mi>
<mo>&le;</mo>
<mi>K</mi>
<mo>}</mo>
</mrow>
其中,父代种群PD中的个体的数量为K,并从x1的取值范围x1,min≤x1≤x1,max内随机取值赋予从x2的取值范围x2,min≤x2≤x2,max内随机取值赋予从而对父代种群PD初始化;
11)从父代种群PD中选出任意对个体,对于每对个体 进行遗传交叉计算或变异计算,并将计算结果赋予子代种群QD中相应的一对个体
12)将父代种群PD与子代种群QD进行合并得到种群R=PD∪QD,即有 将种群R的每个个体与环境变量平均值i=3,…,M合成输入样本并计算相应的多目标适应度函数
13)将种群R的所有个体所对应的适应度函数相互进行比较,将种群R的所有个体划分到具有不同层级的非支配集中,其中,对于层级较低的非支配集中的任一个体所对应的适应度函数objFun1(Xs)和层级较高的非支配集中的任一个体所对应的适应度函数objFun1(Xt)来说,均不存在且而对于同一层级的非支配集中的任两个个体来说,该两个不等式中至少有一个不成立;
14)按照层级从低到高的顺序从所述非支配集中选择K个个体,将选择出的K个个体的值赋予父代种群PD中的个体,并执行步骤11)-步骤14)的过程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;
③求出数据平均值并利用原始数据减去均值得到
④计算协方差矩阵
<mrow>
<mi>C</mi>
<mi>o</mi>
<mi>v</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>144</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>144</mn>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>L</mi>
<mi>k</mi>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>L</mi>
<mi>k</mi>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
⑤计算协方差矩阵的特征值E1,E2,...,E144与特征向量EV1,EV2,...,EV144;
⑥由大到小依次排列特征值E1',E'2,...,E'M,对应特征向量为EV1',EV2',...,EV′144,按特征值大小顺序取前d个特征值的特征向量构成矩阵[EV1',EV2',...,EVd'],此时其中d<144;
⑦由[EV1',EV2',...,EVd']与原始样本求取载荷新的主元,其载荷新主元观测变量构成集合:{αz1,αz2,...,αzd},其为d个新变量,且每个变量为N个观测值构成的新主元矩阵:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>z</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>z</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>z</mi>
<mn>3</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>11</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>12</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>21</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>22</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>z</mi>
<mi>d</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
5.如权利要求1所述的面向抽油机工作过程的智能优化方法,其中,
所述归一化处理的算法如下:
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>min</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>min</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
i=1,2,...,M;k=1,2,...,N
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
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}
于是得到:
<mrow>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>3</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>M</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>&equiv;</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>l</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>&equiv;</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>N</mi>
</msub>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mi>N</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>l</mi>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
6.如权利要求1所述的面向抽油机工作过程的智能优化方法,其中,
所述前馈神经网络由输入层、隐含层和输出层组成,输入层与隐含层由权值W连接,其Wki表示第i个输入神经元与隐含层第k个神经元间的连接权值;隐含层与输出层由权值β连接,其βkj表示第k个隐层神经元与输出层第j个神经元间的连接权值;隐含层的神经元个数为s1,隐含层神经元的阈值为b,bk为第k个隐层神经元阈值;
所述各权值、阈值具体为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>W</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>W</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<mn>1</mn>
<mi>M</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>W</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<mn>2</mn>
<mi>M</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mi>M</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&beta;</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&beta;</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&beta;</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&beta;</mi>
<mrow>
<mn>1</mn>
<mi>l</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&beta;</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&beta;</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&beta;</mi>
<mrow>
<mn>2</mn>
<mi>l</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&beta;</mi>
<mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>&beta;</mi>
<mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&beta;</mi>
<mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mi>l</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>b</mi>
<mo>=</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
设隐含层神经元的激活函数为网络的第m组归一化输入样本的归一化预测输出为:
<mrow>
<msub>
<mover>
<mover>
<mi>Y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mrow>
<mi>l</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</munderover>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>k</mi>
<mn>1</mn>
</mrow>
</msub>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>W</mi>
<mi>k</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</munderover>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>k</mi>
<mn>2</mn>
</mrow>
</msub>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>W</mi>
<mi>k</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</munderover>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>k</mi>
<mi>l</mi>
</mrow>
</msub>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>W</mi>
<mi>k</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Wk=[Wk1,Wk2,…,WkM],最终得到可计算出任一组输入观测值x1,x2,...,xM所对应的输出预测值的前馈神经网络:
<mrow>
<mover>
<mi>Y</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>g</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>g</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msup>
<mi>g</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>.</mo>
</mrow>
7.如权利要求6所述的面向抽油机工作过程的智能优化方法,其中,
利用遗传算法训练网络权值W、权值β、隐含层阈值b的步骤包括:
①利用输入层到隐含层的权值W,隐含层阈值b构建初始种群P,该种群中第i个个体为
②构建遗传算法优化的隐层权值阈值的适应值函数:
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
<mi>f</mi>
<mi>u</mi>
<mi>n</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>Y</mi>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msup>
<mi>&beta;</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
<msup>
<mi>X</mi>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>Y</mi>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msup>
<mi>&beta;</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
<msup>
<mi>X</mi>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msub>
</mfrac>
</mrow>
③计算第i个个体的适应度函数值,利用个体和训练样本输入样本计算隐含层神经元的输出矩阵H(i):
<mrow>
<msup>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msubsup>
<mi>H</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>H</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msubsup>
<mi>H</mi>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>&Sigma;</mi>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
</mover>
<msubsup>
<mi>W</mi>
<mrow>
<mn>1</mn>
<mi>m</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>m</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>b</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>&Sigma;</mi>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
</mover>
<msubsup>
<mi>W</mi>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>m</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>b</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>&Sigma;</mi>
<mover>
<mi>N</mi>
<mo>~</mo>
</mover>
</mover>
<msubsup>
<mi>W</mi>
<mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mi>m</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mi>m</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>b</mi>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
由训练样本输入样本对应的输出样本与隐含层的输出矩阵H(i)共同求出隐含层到输出层的权值β(i),(H(i))-1为矩阵H(i)的逆矩阵;在求出个体对应的β(i)后,构建网络预测函数F(i)(Pi,β(i),X);并将测试样本集Xtest中的输入矩阵带入F(i)(Pi,β(i),X)得到预测测试集预测输出故而计算个体Pi适应度函数值:
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
<mi>f</mi>
<mi>u</mi>
<mi>n</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>Y</mi>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mover>
<mi>Y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>)</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>Y</mi>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mover>
<mi>Y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msub>
</mfrac>
<mo>;</mo>
</mrow>
④对种群P进行选择、交叉、变异操作,得到子群Q;
⑤种群Q进行③、④操作得到种群Q的子群,循环进行以上步骤,得到第GEN代子群,并在第GEN代子群中选择一个最佳个体Pbest,个体Pbest第GEN代子群中计算得到测试样本预测输出值与实际输出值的均方差最小;
⑥由Pbest确定的输入到隐层的权值W,隐层阈值b,计算隐层输出矩阵Hbest,计算隐层到输出的权值矩阵输出权值βbest:从而得到最佳模型。
8.如权利要求7所述的面向抽油机工作过程的智能优化方法,其中,
所述选择操作包括:在P中K个个体中,按每个个体的概率φ(Pi)大小决定选中的机会,分K次从P中选择K个个体,其中概率
所述交叉操作包括:将选择操作得出的K个个体进行二进制编码,个体中每个变量占据4字节,这样编码后一个个体的字节数为:4*(s1*M+s1);将个体Pi和Pj二进制编码后,将两个个体二进制编码第4·c位以后的编码进行调换,从而得到两个新的个体和
所述变异操作包括:将交叉操作得到的K个个体进行二进制编码,其中某个二进制编码后的个体若第c位为1,则将1变成0;若第c位为0,则将0变成1,从而得到一个新的个体Qi。
9.如权利要求1所述的面向抽油机工作过程的智能优化方法,多目标优化的每对决策变量个体之间的遗传交叉计算公式为:
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>Q</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>=</mo>
<mn>0.5</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>&alpha;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>Q</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>=</mo>
<mn>0.5</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>&alpha;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
<msub>
<mi>Q</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>=</mo>
<mn>0.5</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>&alpha;</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
<msub>
<mi>Q</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>=</mo>
<mn>0.5</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>&alpha;</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
随机数α∈[0,1];
每个个体的变异计算公式为:
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>Q</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>Q</mi>
<mi>D</mi>
</msub>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
随机数β∈[0,1]。
10.如权利要求1所述的面向抽油机工作过程的智能优化方法,其中,
在按照层级从低到高的顺序从所述非支配集中选择K个个体时,对于同一层级非支配集中的个体,选择个体拥挤度ds较大的个体,所述个体拥挤度ds的计算方法为:
对当前种群R中所有个体所对应的适应度函数值objFun1(Xs)中的按从小到大的顺序排序,另外,对所有objFun1(Xs)中的按从小到大的顺序排序,令每次排序的第一个和最后一个个体的拥挤距离为无穷大,种群个体的拥挤度ds为
<mrow>
<msub>
<mi>d</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>h</mi>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
</msub>
<mo>-</mo>
<mi>h</mi>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
</msub>
</mrow>
<mrow>
<mi>h</mi>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>h</mi>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>min</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
分别为在所述所排的序列中值的后一个值和前一个值;
分别为在所述所排的序列中值的后一个值和前一个值;
<mrow>
<mi>h</mi>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>{</mo>
<mi>h</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>s</mi>
<mo>&le;</mo>
<mn>2</mn>
<mi>K</mi>
<mo>}</mo>
</mrow>
<mrow>
<mi>h</mi>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>{</mo>
<mi>h</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>s</mi>
<mo>&le;</mo>
<mn>2</mn>
<mi>K</mi>
<mo>}</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>{</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>s</mi>
<mo>&le;</mo>
<mn>2</mn>
<mi>K</mi>
<mo>}</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>{</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>s</mi>
<mo>&le;</mo>
<mn>2</mn>
<mi>K</mi>
<mo>}</mo>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510108394.7A CN104680257B (zh) | 2015-03-12 | 2015-03-12 | 面向抽油机工作过程的智能优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510108394.7A CN104680257B (zh) | 2015-03-12 | 2015-03-12 | 面向抽油机工作过程的智能优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104680257A CN104680257A (zh) | 2015-06-03 |
CN104680257B true CN104680257B (zh) | 2018-03-02 |
Family
ID=53315266
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510108394.7A Active CN104680257B (zh) | 2015-03-12 | 2015-03-12 | 面向抽油机工作过程的智能优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104680257B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106502096B (zh) * | 2016-11-14 | 2019-06-14 | 重庆科技学院 | 基于偏好多目标优化的油田机采过程决策参数优化方法 |
CN107886324A (zh) * | 2017-11-29 | 2018-04-06 | 重庆科技学院 | 一种基于5g的城市一卡通小额支付终端控制方法 |
CN109869137B (zh) * | 2017-12-05 | 2021-06-15 | 中国科学院沈阳自动化研究所 | 一种基于流量计和示功图的抽油井定产模式控制方法 |
CN108678728A (zh) * | 2018-05-15 | 2018-10-19 | 北京理工大学 | 一种基于k-means的油井参数分析组合算法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103177155A (zh) * | 2013-02-28 | 2013-06-26 | 重庆科技学院 | 一种基于bp神经网络与spea2算法的油田抽油机采油节能增产优化方法 |
CN103198354A (zh) * | 2013-02-28 | 2013-07-10 | 重庆科技学院 | 一种利用bp神经网络与nsga2实现油田抽油机采油节能增产优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6941289B2 (en) * | 2001-04-06 | 2005-09-06 | Sas Institute Inc. | Hybrid neural network generation system and method |
-
2015
- 2015-03-12 CN CN201510108394.7A patent/CN104680257B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103177155A (zh) * | 2013-02-28 | 2013-06-26 | 重庆科技学院 | 一种基于bp神经网络与spea2算法的油田抽油机采油节能增产优化方法 |
CN103198354A (zh) * | 2013-02-28 | 2013-07-10 | 重庆科技学院 | 一种利用bp神经网络与nsga2实现油田抽油机采油节能增产优化方法 |
Non-Patent Citations (1)
Title |
---|
面向抽油机节能的GRNN过程建模及工艺参数优化;辜小花 等;《重庆大学学报》;20130630;第36卷(第6期);第130-136页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104680257A (zh) | 2015-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104680025B (zh) | 基于遗传算法极限学习机的抽油机参数优化方法 | |
CN106951983A (zh) | 基于采用多父代遗传算法的人工神经网络的喷射器性能预测方法 | |
CN105046326B (zh) | 基于功图主元分析的抽油机参数优化方法 | |
CN104680257B (zh) | 面向抽油机工作过程的智能优化方法 | |
CN105930901B (zh) | 一种基于rbpnn的变压器故障诊断方法 | |
CN105045941A (zh) | 基于无迹卡尔曼滤波的抽油机参数优化方法 | |
CN107909206A (zh) | 一种基于深层结构循环神经网络的pm2.5预测方法 | |
CN109932903A (zh) | 多父代优化网络和遗传算法的风机控制多目标优化方法 | |
CN109634121A (zh) | 基于径向基函数神经网络的多父代遗传算法空气源热泵多目标优化控制方法 | |
CN103177155B (zh) | 一种基于bp神经网络与spea2算法的油田抽油机采油节能增产优化方法 | |
CN106980897A (zh) | 一种基于变学习率的bp人工神经网络的喷射器性能参数预测方法 | |
CN108954680A (zh) | 一种基于运行数据的中央空调能耗预测方法 | |
CN109598092A (zh) | 融合bp神经网络和多父代遗传算法的空气源热泵多目标优化设计方法 | |
CN108009639A (zh) | 一种基于ga-bp神经网络算法的城市生态建设评价方法 | |
CN106067066A (zh) | 基于遗传算法优化装袋算法的电力变压器故障诊断方法 | |
CN107657243A (zh) | 基于遗传算法优化的神经网络雷达一维距离像目标识别方法 | |
CN103198354B (zh) | 一种利用bp神经网络与非支配排序遗传算法nsga2实现油田抽油机采油节能增产优化方法 | |
CN105931116A (zh) | 基于深度学习机制的自动化信用评分系统及方法 | |
CN113988408B (zh) | 基于多目标规划的地铁站内客流疏散人工诱导方法 | |
CN105574589B (zh) | 基于小生境遗传算法的变压器油色谱故障诊断方法 | |
CN104680259B (zh) | 抽油机运行参数的节能优化方法 | |
CN104680023B (zh) | 基于多目标决策的抽油机参数优化方法 | |
CN110222844A (zh) | 一种基于人工神经网络的压缩机性能预测方法 | |
CN106126607A (zh) | 一种面向社交网络的用户关系分析方法 | |
CN109492748A (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 | ||
CP03 | Change of name, title or address |
Address after: 300280 building 7, South Ring Road, Jinnan, Tianjin, 7 Co-patentee after: Chongqing University of Science & Technology Patentee after: Tianjin people's Polytron Technologies Inc Address before: 300280 Tianjin Grand Port Dagang Oilfield port West Transport Shuangfeng District Co-patentee before: Chongqing University of Science & Technology Patentee before: Tianjin Wanzhong Technology Development Co., Ltd. |
|
CP03 | Change of name, title or address |