CN104680025B - 基于遗传算法极限学习机的抽油机参数优化方法 - Google Patents

基于遗传算法极限学习机的抽油机参数优化方法 Download PDF

Info

Publication number
CN104680025B
CN104680025B CN201510109219.XA CN201510109219A CN104680025B CN 104680025 B CN104680025 B CN 104680025B CN 201510109219 A CN201510109219 A CN 201510109219A CN 104680025 B CN104680025 B CN 104680025B
Authority
CN
China
Prior art keywords
mtd
mrow
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
Application number
CN201510109219.XA
Other languages
English (en)
Other versions
CN104680025A (zh
Inventor
辜小花
李太福
梅青平
裴仰军
王坎
周伟
杨永龙
曹旭鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Daqing Suofelectronic Technology Development Co ltd
Original Assignee
Chongqing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN201510109219.XA priority Critical patent/CN104680025B/zh
Publication of CN104680025A publication Critical patent/CN104680025A/zh
Application granted granted Critical
Publication of CN104680025B publication Critical patent/CN104680025B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明提供一种基于遗传算法极限学习机的抽油机参数优化方法，包括：1)确定抽油机的生产效率影响因素和性能变量、2)获得其样本数据、3)对载荷数据降维、4)由非载荷变量与载荷新主元构建网络输入变量集、5)构建输入变量的样本值、6)归一化处理、7)选取训练和测试样本、8)构建前馈神经网络并训练、9)构造产液量偏好函数、10)构建评价决策变量个体支配关系的适应度函数、11)计算系统环境变量的平均值、12)利用决策变量构建初始种群并生成精英解种群、13)进行遗传迭代计算，得到第二代精英种群和父代种群、14)迭代循环，最终得到优化后的生产效率影响因素值。优化后，可保证在产液量基本固定的情况下耗电量最小。

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)在归一化后样本集中，选取样本集前组样本作为网络训练样本集则样本集剩余组作为测试样本集其中其中：训练样本集和测试样本集分别为：
8)构建3层前馈神经网络，采用所述训练输入样本对该网络进行训练，利用遗传算法训练网络权值W、权值β、隐含层阈值b，从而确定网络的参数，最终得到可计算出任一组输入观测值x1,x2,...,xM所对应的输出预测值的前馈神经网络：
9)针对产液量y1构造其偏好函数h＝h(y1)，该偏好函数为U形曲线，在其整个定义域上二阶可导且二阶导数恒大于零，将产液量y1的值划分为好、较好、一般、较差和极差5个区域，并通过所述偏好函数将该五个区域量化为数值h＝h(y1)；其中，产液量y1处于好区域对应偏好函数最小值；
10)构建评价决策变量个体支配关系的适应度函数，系统的性能变量选取产液量(y1)、耗电量(y2)，结合步骤S9构建的产液量偏好函数，得到适应度函数如下：
11)计算抽油机工艺系统环境变量的平均值，以作为优化决策参数时的环境状态；
12)利用决策变量x1,x2构建多目标优化的初始种群PD的个体，即设置决策变量的上下限xmin、xmax,即x1min≤x1≤x1max，x2min≤x2≤x2max，利用决策变量(x1,x2)的K对数据构建初始种群PD，即初始化种群PD,令其为第一代父代种群；生成一个空的精英解种群设置精英个体个数为设置最大遗传代数GEN＝100；
13)进行第一次遗传迭代计算，并得到第二代精英种群A2、第二代父代种群P2。具体步骤如下：
①个体强度求取，将第一代父代种群与精英解种群组合成种群Rt，即求取种群Rt中每个个体的原始适应度函数值，并比较个体之间的相互支配关系；定义变量R(c)为种群Rt中第c个个体强度，即第c个个体可以被种群Rt其他个体支配的数量；其中个体Rt(c)原始适应度函数值求取过程如下：通过种群个体Rt(c)与环境状态变量平均值组建输入样本计算样本Xc原始适应度函数值并作为个体Rt(c)的原始适应度函数值；
②个体密度求取，利用个体Rt(i)与种群Rt中第b个邻近个体的距离值则个体Rt(c)密度函数
③求个体的适应值，将上述所求个体Rt(c)的强度R(c)和所求个体Rt(c)的密度值D(c)的相加作为个体Rt(c)的适应值；
④在种群Rt将所有的非支配个体全部放入精英种群A2，但要保持精英种群个体数为此时存在三种情况：A2中个体数为则不需要在操作，如果A2中个体数小于则需要在种群Rt剩余个体中选取适应值较小的个体放入A2中，保持A2个体数为如果A2中个体数大于则需要在A2个体中比较个体的密度值D(c)，将密度值较大的个体剔除，以保持A2个体数为
⑤将A2中个体放入交配池中进行遗传操作得到第二代父代种群
⑥将第二代父代种群P2与第二代精英种群A2组合，并重复①～⑤过程，直至gen＝GEN，输出精英种群AGEN,将AGEN的个体作为优化结果；
14)将优化后的决策变量，以及环境变量的平均值带入建立工艺过程模型，计算优化后的决策变量的系统性能，该优化后的决策变量取值可保证在固定产液量的情形下，耗电量降低。
本发明的有益效果是，利用遗传算法优化的极限学习机算法(GAELM)建立油田机采过程的高精度模型，并利用建立油田机采过程的高精度模型，并利用智能特性的进化算法改进强度Pareto进化算法(Strength Pareto evolutionary algorithm 2,SPEA2)对建立的模型进行搜索，探寻抽油机生产过程中最佳工艺决策参数，给出面向节能降耗的抽油机生产最佳参数，指导生产。通过选择优化后的运行参数，可以使抽油机在运行过程中保证在产液量基本固定且具有最小偏好值的情况下，耗电量最小，从而可以降低油田生产成本并提高油田生产效率。
附图说明
图1示出了抽油机的工作模型；
图2示出了本发明一个实施例所述的基于遗传算法极限学习机的抽油机参数优化方法的流程图；
图3示出了本发明一个实施例中的前馈神经网络的结构；
图4示出了本发明一个实施例中的产液量偏好函数的图形；
图5示出了利用本发明的方法训练神经网络所产生的产液量训练效果图；
图6示出了利用本发明的方法训练神经网络所产生的耗电量训练效果图；
图7示出了测试样本产液量预测效果图；
图8示出了测试样本耗电量预测效果图；
图9示出了偏好函数的图形；
图10示出了产液量偏好值与耗电量的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。
步骤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：由影响因素观测变量集合中非载荷变量与载荷新主元观测变量集合{αz1,αz2,...,αzd}构建网络输入变量集合：并令输入变量集合为：{x1,x2,x3,...,xM}，即，
步骤S5：构建输入变量集合{x1,x2,x3,...,xM}观测样本值：
其中，x1～x2为决策变量，x3～xM为新的环境变量。
步骤S6：对得到的训练输入样本X、输出样本Y进行归一化处理，得到新的训练输入矩阵输出矩阵
在一个实施例中，所述归一化处理的算法如下：
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}
于是得到：
步骤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个体L＝50，设置最大遗传代数GEN＝100。
①构建遗传算法优化网络权值阈值的适应值函数：
③计算第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：βbest＝(Hbest)-1·Ytrain从而得到最佳模型。
上述选择操作可以包括：在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：构建评价决策变量个体支配关系的适应度函数，系统的性能变量选取产液量(y1)、耗电量(y2)，结合步骤S9构建的产液量偏好函数，得到适应度函数如下：
步骤S11：计算抽油机工艺系统环境变量的平均值，以作为优化决策参数时的环境状态。
计算环境参数平均值具体算法如下：
其中N为该环境变量训练样本数量。
步骤S12：利用决策变量x1,x2构建多目标优化的初始种群PD的个体，即设置决策变量的上下限xmin、xmax,即x1min≤x1≤x1max，x2min≤x2≤x2max，利用决策变量(x1,x2)的K对数据构建初始种群PD，即初始化种群PD,令其为第一代父代种群；生成一个空的精英解种群设置精英个体个数为设置最大遗传代数GEN＝100。
步骤S13：进行第一次遗传迭代计算，并得到第二代精英种群A2、第二代父代种群具体步骤如下：
①个体强度求取，将第一代父代种群与精英解种群组合成种群Rt，即求取种群Rt中每个个体的原始适应度函数值，并比较个体之间的相互支配关系；定义变量R(c)为种群Rt中第c个个体强度，即第c个个体可以被种群Rt其他个体支配的数量；其中个体Rt(c)原始适应度函数值求取过程如下：通过种群个体Rt(c)与环境状态变量平均值组建输入样本计算样本Xc原始适应度函数值并作为个体Rt(c)的原始适应度函数值；
②个体密度求取，利用个体Rt(i)与种群Rt中第b个邻近个体的距离值则个体Rt(c)密度函数
③求个体的适应值，将上述所求个体Rt(c)的强度R(c)和所求个体Rt(c)的密度值D(c)的相加作为个体Rt(c)的适应值；
④在种群Rt将所有的非支配个体全部放入精英种群A2，但要保持精英种群个体数为此时存在三种情况：A2中个体数为则不需要在操作，如果A2中个体数小于则需要在种群Rt剩余个体中选取适应值较小的个体放入A2中，保持A2个体数为如果A2中个体数大于则需要在A2个体中比较个体的密度值D(c)，将密度值较大的个体剔除，以保持A2个体数为
⑤将A2中个体放入交配池中进行遗传操作得到第二代父代种群
⑥将第二代父代种群与第二代精英种群A2组合，并重复①～⑤过程，直至gen＝GEN，输出精英种群AGEN,将AGEN的个体作为优化结果。
步骤S14：将优化后的决策变量，以及环境变量的平均值带入建立工艺过程模型，计算优化后的决策变量的系统性能，该优化后的决策变量取值可保证在固定产液量的情形下，耗电量降低。
下面以大港油田港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＝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对数据构建初始种群P，设置种群大小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.75，此时产液量可能为49.51t/d或者45.08t/d。此时计算得到在保持定量式采油工作状态下，耗电量下降16.57％。系统效率比(产液量/耗电量)可提高8％。达到了节能设计目的。说明该方法有效。
尽管已经结合详细示出并描述的优选实施例公开了本发明，但是本领域技术人员应当理解，对于上述本发明所提出的基于遗传算法极限学习机的抽油机参数优化方法，还可以在不脱离本发明内容的基础上做出各种改进。因此，本发明的保护范围应当由所附的权利要求书的内容确定。

Claims (9)

1.一种基于遗传算法极限学习机的抽油机参数优化方法，包括如下步骤：
1)确定抽油机采油过程生产效率影响因素构成效率观测变量集合其中α1，α2为决策变量，α3～α146载荷数据环境变量，为其他环境变量，选取抽油机系统的性能变量构成性能观测变量集合：{y1,y2,y3,…yl}；
2)获得所述生产效率影响因素和系统性能变量的样本数据，得到效率影响因素样本矩阵α和性能样本矩阵Y：
<mrow> <mi>&amp;alpha;</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>12</mn> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mn>1</mn> <mi>N</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>22</mn> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;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>&amp;alpha;</mi> <mrow> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mi>N</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;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> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> </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>
其中为效率影响因素个数，N为样本个数，αik表示第i个效率影响因素变量的第k个观测值，i＝1,2,...,M；k＝1,2,...,N；
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>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>11</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>12</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>1</mn> <mi>N</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>21</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>22</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;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>&amp;alpha;</mi> <mrow> <mi>z</mi> <mi>d</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mi>d</mi> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mi>d</mi> <mi>N</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
4)由影响因素观测变量集合中非载荷变量与载荷新主元观测变量集合{αz1,αz2,...,α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> <mo>...</mo> </mtd> <mtd> <msub> <mi>X</mi> <mi>N</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>12</mn> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mn>1</mn> <mi>N</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mn>22</mn> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mn>2</mn> <mi>N</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mn>147</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mn>147</mn> <mo>,</mo> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;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>&amp;alpha;</mi> <mrow> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mi>N</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>11</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>12</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;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>&amp;alpha;</mi> <mrow> <mi>z</mi> <mi>d</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mi>d</mi> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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> <mn>...</mn> </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> </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> </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> <mo>...</mo> </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> </mrow>
9)针对产液量y1构造其偏好函数h＝h(y1)，该偏好函数为U形曲线，在其整个定义域上二阶可导且二阶导数恒大于零，将产液量y1的值划分为好、较好、一般、较差和极差5个区域，并通过所述偏好函数将该五个区域量化为数值h＝h(y1)；其中，产液量y1处于好区域对应偏好函数最小值；
10)构建评价决策变量个体支配关系的适应度函数，系统的性能变量选取产液量(y1)、耗电量(y2)，结合步骤S9构建的产液量偏好函数，得到适应度函数如下：
<mrow> <mi>o</mi> <mi>b</mi> <mi>j</mi> <mi>F</mi> <mi>u</mi> <mi>n</mi> <mn>1</mn> <mrow> <mo>(</mo> <mi>X</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mn>2</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>h</mi> <mrow> <mo>(</mo> <msup> <mi>g</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <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> <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> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>;</mo> </mrow>
11)计算抽油机工艺系统环境变量的平均值，以作为优化决策参数时的环境状态；
12)利用决策变量x1,x2构建多目标优化的初始种群PD的个体，即设置决策变量的上下限xmin、xmax,即x1min≤x1≤x1max，x2min≤x2≤x2max，利用决策变量(x1,x2)的K对数据构建初始种群PD，即初始化种群PD,令其为第一代父代种群；生成一个空的精英解种群设置精英个体个数为设置最大遗传代数GEN＝100；
13)进行第一次遗传迭代计算，并得到第二代精英种群A2、第二代父代种群具体步骤如下：
①个体强度求取，将第一代父代种群与精英解种群组合成种群Rt，即求取种群Rt中每个个体的原始适应度函数值，并比较个体之间的相互支配关系；定义变量R(c)为种群Rt中第c个个体强度，即第c个个体可以被种群Rt其他个体支配的数量；其中个体Rt(c)原始适应度函数值求取过程如下：通过种群个体Rt(c)与环境状态变量平均值组建输入样本计算样本Xc原始适应度函数值并作为个体Rt(c)的原始适应度函数值；
②个体密度求取，利用个体Rt(i)与种群Rt中第b个邻近个体的距离值则个体Rt(c)密度函数
③求个体的适应值，将上述所求个体Rt(c)的强度R(c)和所求个体Rt(c)的密度值D(c)的相加作为个体Rt(c)的适应值；
④在种群Rt将所有的非支配个体全部放入精英种群A2，但要保持精英种群个体数为此时存在三种情况：A2中个体数为则不需要在操作，如果A2中个体数小于则需要在种群Rt剩余个体中选取适应值较小的个体放入A2中，保持A2个体数为如果A2中个体数大于则需要在A2个体中比较个体的密度值D(c)，将密度值较大的个体剔除，以保持A2个体数为
⑤将A2中个体放入交配池中进行遗传操作得到第二代父代种群
⑥将第二代父代种群与第二代精英种群A2组合，并重复①～⑤过程；直至gen＝GEN，输出精英种群AGEN,将AGEN的个体作为优化结果；
14)将优化后的决策变量，以及环境变量的平均值带入建立工艺过程模型，计算优化后的决策变量的系统性能，该优化后的决策变量取值可保证在固定产液量的情形下，耗电量降低。
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；
⑥由大到小依次排列特征值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>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>11</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>12</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>1</mn> <mi>N</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>21</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mn>22</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;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>&amp;alpha;</mi> <mrow> <mi>z</mi> <mi>d</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mi>d</mi> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>z</mi> <mi>d</mi> <mi>N</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
5.如权利要求1所述的基于遗传算法极限学习机的抽油机参数优化方法，其中，
所述归一化处理的算法如下：
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <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> <mo>(</mo> <mrow> <msub> <mover> <mi>x</mi> <mo>~</mo> </mover> <mi>max</mi> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>~</mo> </mover> <mi>min</mi> </msub> </mrow> <mo>)</mo> <mo>&amp;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>max</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>min</mi> </mrow> </msub> </mrow> </mfrac> <mo>+</mo> <msub> <mover> <mi>x</mi> <mo>~</mo> </mover> <mi>min</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>M</mi> <mo>;</mo> <mi>k</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>N</mi> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <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> <mo>(</mo> <mrow> <msub> <mover> <mi>y</mi> <mo>~</mo> </mover> <mi>max</mi> </msub> <mo>-</mo> <msub> <mover> <mi>y</mi> <mo>~</mo> </mover> <mi>min</mi> </msub> </mrow> <mo>)</mo> <mo>&amp;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>min</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> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>....</mn> <mo>,</mo> <mi>l</mi> <mo>;</mo> <mi>k</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>N</mi> </mrow> </mtd> </mtr> </mtable> </mfenced>
其中:为设定输入变量归一化后数据范围的最大值、最小值；
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>&amp;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>&amp;equiv;</mo> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> </mrow> 4
<mrow> <mover> <mi>X</mi> <mo>~</mo> </mover> <mo>=</mo> <mo>&amp;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>&amp;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> <mn>...</mn> </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>&amp;beta;</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>&amp;beta;</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;beta;</mi> <mn>12</mn> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;beta;</mi> <mrow> <mn>1</mn> <mi>l</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;beta;</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>&amp;beta;</mi> <mn>22</mn> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;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>&amp;beta;</mi> <mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;beta;</mi> <mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>&amp;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>&amp;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>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> </munderover> <msub> <mi>&amp;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>&amp;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>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> </munderover> <msub> <mi>&amp;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>&amp;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>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> </munderover> <msub> <mi>&amp;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>&amp;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> <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> <mo>...</mo> </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>
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>&amp;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>&amp;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>&amp;lsqb;</mo> <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> <mo>...</mo> <mo>,</mo> <msubsup> <mi>H</mi> <msub> <mi>s</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> </msubsup> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> <mo>=</mo> <mo>&amp;lsqb;</mo> <mi>g</mi> <mrow> <mo>(</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mover> <mi>N</mi> <mo>~</mo> </mover> </munderover> <msubsup> <mi>W</mi> <mrow> <mn>1</mn> <mi>m</mi> </mrow> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> </msubsup> <mo>&amp;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> <mo>)</mo> </mrow> <mi>g</mi> <mrow> <mo>(</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mover> <mi>N</mi> <mo>~</mo> </mover> </munderover> <msubsup> <mi>W</mi> <mrow> <mn>1</mn> <mi>m</mi> </mrow> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> </msubsup> <mo>&amp;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> <mo>)</mo> </mrow> <mo>...</mo> <mi>g</mi> <mrow> <mo>(</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mover> <mi>N</mi> <mo>~</mo> </mover> </munderover> <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>&amp;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> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </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：βbest＝(Hbest)-1·Ytrain从而得到最佳模型。
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> <msub> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>3</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>M</mi> </mrow>
其中N为该环境变量训练样本数量。
CN201510109219.XA 2015-03-12 2015-03-12 基于遗传算法极限学习机的抽油机参数优化方法 Active CN104680025B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510109219.XA CN104680025B (zh) 2015-03-12 2015-03-12 基于遗传算法极限学习机的抽油机参数优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510109219.XA CN104680025B (zh) 2015-03-12 2015-03-12 基于遗传算法极限学习机的抽油机参数优化方法

Publications (2)

Publication Number Publication Date
CN104680025A CN104680025A (zh) 2015-06-03
CN104680025B true CN104680025B (zh) 2017-08-25

Family

ID=53315058

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510109219.XA Active CN104680025B (zh) 2015-03-12 2015-03-12 基于遗传算法极限学习机的抽油机参数优化方法

Country Status (1)

Country Link
CN (1) CN104680025B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106529042B (zh) * 2016-11-14 2019-12-31 重庆科技学院 基于计算智能的油田机采参数动态演化建模与优化方法
CN106502096B (zh) * 2016-11-14 2019-06-14 重庆科技学院 基于偏好多目标优化的油田机采过程决策参数优化方法
CN108804720B (zh) * 2017-04-26 2021-09-14 重庆科技学院 基于改进无迹Kalman滤波与RBF神经网络的抽油机故障诊断方法
CN107346459B (zh) * 2017-05-22 2020-09-18 天津科技大学 一种基于遗传算法改进的多模式污染物集成预报方法
CN109376934A (zh) * 2018-10-31 2019-02-22 广东兴发铝业有限公司 一种基于大数据的挤压机的能耗预测方法
CN111832799A (zh) * 2020-05-25 2020-10-27 山东电力工程咨询院有限公司 一种燃煤发电机组节能改造性能优化方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102982373A (zh) * 2012-12-31 2013-03-20 山东大学 一种混合svm回归算法的oin神经网络训练方法
CN103914622A (zh) * 2014-04-04 2014-07-09 清华大学 一种化学品泄漏快速预测预警应急响应决策方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4475625B2 (ja) * 2003-05-27 2010-06-09 キヤノン株式会社 画像処理装置及び画像処理方法
CN103543889B (zh) * 2012-07-17 2017-12-29 深圳达沃斯光电有限公司 单层电容式触摸屏制作方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102982373A (zh) * 2012-12-31 2013-03-20 山东大学 一种混合svm回归算法的oin神经网络训练方法
CN103914622A (zh) * 2014-04-04 2014-07-09 清华大学 一种化学品泄漏快速预测预警应急响应决策方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《Application of a single-speed three-power energy-saving motor in pumping unit》;Wang H F，et al。;《Electric Machine & Control Application》;20101030;第50-54页 *
《多层前馈神经网络改进算法及其应用》;宋宜斌;《计算机工程》;20030831;第109-111页 *
《面向抽油机节能的ＧＲＮＮ过程建模及工艺参数优化》;辜小花，等;《重庆大学学报》;20130630;第36卷(第6期);第130-136页 *

Also Published As

Publication number Publication date
CN104680025A (zh) 2015-06-03

Similar Documents

Publication Publication Date Title
CN104680025B (zh) 基于遗传算法极限学习机的抽油机参数优化方法
CN105045941B (zh) 基于无迹卡尔曼滤波的抽油机参数优化方法
CN105930901B (zh) 一种基于rbpnn的变压器故障诊断方法
CN105046326B (zh) 基于功图主元分析的抽油机参数优化方法
CN106951983A (zh) 基于采用多父代遗传算法的人工神经网络的喷射器性能预测方法
CN103198354B (zh) 一种利用bp神经网络与非支配排序遗传算法nsga2实现油田抽油机采油节能增产优化方法
CN107909206A (zh) 一种基于深层结构循环神经网络的pm2.5预测方法
CN106980897A (zh) 一种基于变学习率的bp人工神经网络的喷射器性能参数预测方法
CN108009639A (zh) 一种基于ga-bp神经网络算法的城市生态建设评价方法
CN104680023B (zh) 基于多目标决策的抽油机参数优化方法
CN104680257B (zh) 面向抽油机工作过程的智能优化方法
CN107016469A (zh) 电力负荷预测方法
CN106067066A (zh) 基于遗传算法优化装袋算法的电力变压器故障诊断方法
CN105931116A (zh) 基于深度学习机制的自动化信用评分系统及方法
CN107657243A (zh) 基于遗传算法优化的神经网络雷达一维距离像目标识别方法
CN109063247A (zh) 基于深度信念网络的滑坡灾害预报方法
CN110222844A (zh) 一种基于人工神经网络的压缩机性能预测方法
Sahu et al. Economic load dispatch in power system using genetic algorithm
CN106447133A (zh) 一种基于深度自编码网络的短期电力负荷预测方法
CN109492748A (zh) 一种基于卷积神经网络的电力系统的中长期负荷预测模型建立方法
CN106126607A (zh) 一种面向社交网络的用户关系分析方法
CN105574589A (zh) 基于小生境遗传算法的变压器油色谱故障诊断方法
CN106096723A (zh) 一种基于混合神经网络算法的用于复杂工业产品性能评估方法
CN101719193B (zh) 一种制动闸寿命的预测方法
CN104680259B (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: 20230417

Address after: No. 112 Yuanwang Street, Industrial Development Community, Ranghulu District, Daqing City, Heilongjiang Province, 163712

Patentee after: Daqing suofelectronic Technology Development 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