CN106777465A - 高含硫天然气净化工艺动态演化建模与节能优化方法 - Google Patents

Abstract

本发明提供的高含硫天然气净化工艺动态演化建模与节能优化方法，包括选取影响脱硫效率的工艺参数和脱硫单元的性能指标后进行采集形成样本集；对样本集进行归一化形成归一化样本集，并从中选取训练样本和测试样本；基于训练样本构建神经网络模型并确定神经网络模型的初始状态变量；用ST‑UPFNN算法估计神经网络模型的最优状态变量；根据最优状态变量对神经网络模型进行更新；分别构建H₂S浓度和CO₂浓度的偏好函数；利用SPEA‑II算法对H₂S浓度和CO₂浓度的工艺参数各自的上下限进行优化，将优化后的工艺参数带入更新后的神经网络模型，计算优化后的工艺参数的系统性能，与实际样本的系统性能的平均值进行比较。利用本发明可以提高高含硫天然气净化的生产效率。

Description

高含硫天然气净化工艺动态演化建模与节能优化方法

技术领域

本发明涉及高含硫天然气净化技术领域，更为具体地，涉及一种高含硫天然气净化工艺动态演化建模与节能优化方法。

背景技术

高含硫天然气酸性组分含量比常规天然气高出数倍，其脱硫过程胺液循环量大、工艺流程复杂、能耗高。统计表明，脱硫单元能耗占高含硫天然气净化厂总能耗50％以上，其单位综合能耗高达1729.3MJ·t-1，属于高耗能单元。对大型净化厂而言，通过脱硫单元优化可降低能耗5％～10％。此外，高含硫天然气酸性组分浓度高，经过净化后的产品气量相对原料气流量有显著下降。为此，对高含硫天然气脱硫过程进行工艺优化，实现节能降耗，提高产率和气体加工经济效益是十分必要的。

发明内容

鉴于上述问题，本发明的目的是提供一种高含硫天然气净化工艺动态演化建模与节能优化方法，以解决上述背景技术所提出的问题。

本发明提供的高含硫天然气净化工艺动态演化建模与节能优化方法，包括：

步骤S1：选择影响脱硫效率的工艺参数和脱硫单元的性能指标；其中，工艺参数包括进入尾气吸收塔贫的胺液流量x₁、进入二级吸收塔的贫胺液流量x₂、原料气处理量x₃、尾气单元返回脱硫单元的半富胺液流量x₄、一级吸收塔胺液入塔温度x₅、二级吸收塔胺液入塔温度x₆、闪蒸罐压力x₇、一个重沸器的蒸汽消耗量x₈、另一个重沸器的蒸汽消耗量x₉和蒸汽预热器的蒸汽消耗量x₁₀；脱硫单元的性能指标包括净化气中H₂S的浓度和CO₂的浓度

步骤S2：采集预设时间的工艺参数和所述性能指标的数据，剔除误差样本后形成样本集[X,Y]；

步骤S3：对样本集[X,Y]进行归一化，形成归一化样本集取所述归一化样本集中前80％的样本作为训练样本，而剩余的20％样本作为测试样本；

步骤S4：基于所述训练样本训练构建神经网络模型并根据模型确定初始状态变量X，将训练样本中的作为所述神经网络模型的输入，将训练样本中的作为神经网络模型的输出；

其中，神经网络模型为：

其中，I_k为训练样本的矢量样本值，并作为神经网络模型的输入，为网络输入层到隐含层的神经元的连接权值，为网络输入层到隐含层的神经元的阈值，为隐含层到网络输出层的神经元的连接权值，为隐含层到网络输出层的神经元的阈值，其中，i＝1，2…S₀；j＝1，2…S₁；k＝1，2…S₂；S₀为网络输入层的神经元的数量，S₁为网络隐含层的神经元的数量，S₂为网络输出层的神经元的数量；

初始状态变量为：

步骤S5：利用ST-UPFNN算法估计神经网络模型的最优状态变量；

步骤S6：将最优状态变量作为神经网络模型的和对式(1)进行更新，获得权值阈值更新后的神经网络模型；

步骤S7：分别构建H₂S浓度的偏好函数和CO₂浓度的偏好函数；

步骤S8：利用SPEA-II算法对H₂S浓度的偏好函数和CO₂浓度的偏好函数进行多目标极值寻优优化，获得符合生产实际的决策变量；

步骤S9：将优化后的决策变量带入权值阈值更新后的神经网络模型，计算优化后的决策变量的系统性能的平均值，与实际样本的系统性能的平均值进行比较，如果优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值，利用优化后的决策变量对实际生产进行指导；否则重复上述步骤S1-S8，直至优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值为止。

本发明提供的高含硫天然气净化工艺动态演化建模与节能优化方法，能够有效地跟踪装置的整个经济效益最优路线，有效克服过程干扰、设备性能变化、经济效益和生产目标的变化问题。

附图说明

通过参考以下结合附图的说明及权利要求书的内容，并且随着对本发明的更全面理解，本发明的其它目的及结果将更加明白及易于理解。在附图中：

图1为H₂S浓度训练样本的训练效果图；

图2为H₂S浓度训练样本的训练误差图；

图3为H₂S浓度测试样本的测试效果图；

图4为H₂S浓度测试样本的测试误差图；

图5为CO₂浓度训练样本的训练效果图；

图6为CO₂浓度训练样本的训练误差图；

图7为CO₂浓度测试样本的测试效果图；

图8为CO₂浓度测试样本的测试误差图；

图9a-图9b为净化气中H₂S含量和CO₂含量的偏好函数图；

图10为H₂S浓度的偏好值与CO₂浓度的偏好值的Pareto解集图；

图11为H₂S浓度的实际值与CO₂浓度的实际值的Pareto解集图。

具体实施方式

名称解释

ST-UKFNN：Strong TrackUnscented Kalman FilterNeural Network，强追踪无迹卡尔曼滤波神经网络。

ST-UPFNN：Strong TrackUnscented Particle FilterNeural Network，强追踪无迹粒子滤波神经网络，其将ST-UKFNN、粒子滤波(Particle Filter)、BP神经网络相结合。

SPEA-II：Strength pareto evolutionary algorithm-II，改进的强度pareto进化算法。

本发明提供的高含硫天然气净化工艺动态演化建模与节能优化方法，包括：

步骤S1：选择影响脱硫效率的工艺参数和脱硫单元的性能指标；其中，工艺参数包括进入尾气吸收塔贫的胺液流量、进入二级吸收塔的贫胺液流量、原料气处理量、尾气单元返回脱硫单元的半富胺液流量、一级吸收塔胺液入塔温度、二级吸收塔胺液入塔温度、闪蒸罐压力、一个重沸器的蒸汽消耗量、另一个重沸器的蒸汽消耗量和蒸汽预热器的蒸汽消耗量；脱硫单元的性能指标包括净化气中H₂S和CO₂的浓度以及净化气的产量。如表1所示：

表1

步骤S2：采集预设时间的工艺参数和性能指标的数据，剔除误差样本后形成样本集[X，Y]。样本集[X，Y]如下表2所示：

表2

步骤S3：对样本集[X，Y]进行归一化，形成归一化样本集取归一化样本集中前80％的样本作为训练样本，而剩余的20％样本作为测试样本。

步骤S4：基于训练样本构建神经网络模型和神经网络模型的初始状态变量X，以及，将训练样本中的作为神经网络模型的输入即将训练样本中的作为神经网络模型的输出即

其中，初始状态变量为：

步骤S5：利用ST-UPFNN算法估计神经网络模型的最优状态变量。

本发明利用ST-UPFNN算法估计所述神经网络模型的状态变量，以达到连接权值、阈值的不断调整，直到满足要求。将得到的最优状态变量作为上述所建立神经网络模型的连接权值、阈值。需要说明的是，该连接权值、阈值为通过ST-UPFNN算法调整后的连接权值、阈值，也是上述所建立的神经网络模型的全部连接权值与阈值，包括和

利用ST-UPFNN算法估计神经网络模型的最优状态变量的过程包括：

步骤S51：针对粒子滤波器设置粒子的数目N，并对抽样，得到初始粒子集并将所述初始粒子集中的每个粒子的权值均设为1/N。

其中，表示以x₀为均值、P₀为方差的正态分布采。

步骤S52：在获取(k+1)时刻的观测变量值后，利用ST-UKFNN算法对每个粒子进行(k+1)时刻的状态估计，得到最优状态估计值和协方差

利用ST-UKFNN算法对每个粒子进行状态估计的过程，包括：

步骤S521：对所述初始状态变量X进行Sigma采样，获得2n+1个采样点，初始化控制2n+1个采样点的分布状态参数α、待选参数κ，以及非负权系数β，对所述初始状态变量X的Sigma采样如下：

步骤S522：计算每个采样点的权重，每个采样点的权重如下：

其中，W_c为计算状态变量的协方差的权重，W_m为计算状态估计和观测预测时的权重，是的第一列，是的第一列。

步骤S523：通过离散时间非线性系统的状态方程将每个采样点的k时刻的最优状态变量的状态估计变换为(k+1)时刻的状态变量的状态估计以及，通过合并(k+1)时刻的状态估计的向量，获得(k+1)时刻的状态变量的状态先验估计和协方差P_k+1|k；其中，

(k+1)时刻的状态变量的状态估计为：

其中，表示k时刻的最优估计，w_k为过程噪声，其协方差矩阵Q_k为cov(w_k,w_j)＝Q_kδ_kj，

(k+1)时刻的状态变量的状态先验估计为：

(k+1)时刻的状态变量的协方差P_k+1|k为：

步骤S524：通过离散时间非线性系统的观测方程将获得的(k+1)时刻的状态变量的状态估计转化为(k+1)时刻的观测预测

其中，ν_k为观测噪声，其协方差矩阵R_k为cov(v_k,v_j)＝R_kδ_kj，

步骤S525：通过并估计(k+1)时刻的观测预测的向量，获得(k+1)时刻的先验观测预测并根据先验观测预测估计(k+1)时刻观测预测的协方差

(k+1)时刻的先验观测预测为：

(k+1)时刻观测预测的预测协方差为：

其中，在此处引入强追踪算法，即渐消因子λ_k+1增强模型的追踪能力以提高模型精度；

N_k+1＝V_k+1-βR_k+1 (12)

其中，β为弱化因子，β≥1；

步骤S526：计算(k+1)时刻的状态变量的状态先验估计与(k+1)时刻的先验观测预测之间的状态协方差状态协方差为：

步骤S527：通过建立状态协方差和预测协方差的关系，更新(k+1)时刻的状态变量的状态估计和协方差，分别获得(k+1)时刻的最优状态估计值和协方差

其中，建立的状态协方差和预测协方差的关系为：

其中，K_k+1为增益矩阵，以及，

更新后的(k+1)时刻的状态变量的状态估计为：

更新后的(k+1)时刻的状态变量的协方差P_k+1为：

将更新后的(k+1)时刻的状态变量的状态估计和协方差P_k+1分别作为(k+1)时刻的最优状态估计值和协方差

步骤S53：将最优状态估计值和初始化后的协方差作为粒子的重要性密度函数进行抽样，得到新粒子由所有新粒子组成的粒子集中的每个新粒子的正态分布概率密度值如下：

其中，rand为正态分布随机误差

正态分布密度函数为：

其中，N为每个新粒子的正态分布，x、μ、σ分别为正态分布函数的三个变量，p为每个新粒子的条件概率；

式(2)中，x,μ,σ分别一一对应

式(3)中，x,μ,σ分别一一对应

式(4)中，x,μ,σ分别一一对应

步骤S54：对新粒子的权值进行更新，并进行归一化处理；其中，

权值更新公式为：

权值归一化公式为：

步骤S55：根据新粒子归一化处理后的权值和重采样策略对粒子集进行重采样，获取新粒子集并求取新粒子集中每个新粒子的状态估计值

设变量u，令取u₁∈(0,1)

步骤S56：以粒子的数目N作为循环次数循环步骤S51-步骤S55的计算过程，将最后一次估计得到系统状态变量作为利用ST-UPFNN算法估计神经网络模型的最优状态变量；其中，将新粒子的状态估计值作为本时刻的最优估计付给进行下一时刻的状态估计。

神经网络模型的最优状态变量的结构参数如下：

(1)H₂S浓度模型参数值

w₂(1×15)＝[-9.88 12.31 ... -6.53 -13.39] b₂(1×1)＝[0.27]

(2)CO₂浓度模型参数值

w₂(1×15)＝[10.43 -9.29 ... -4.74 -11.14] b₂(1×1)＝[-4.54]

步骤S6：将最优状态变量作为所述神经网络模型的和对式(1)进行更新，获得权值阈值更新后的神经网络模型。

利用测试样本测试权值阈值更新后的神经网络模型的预测精度。H₂S浓度模型的测试结果如图1-图4所示，CO₂浓度模型的测试结果如图5-图8所示。

步骤S7：分别构建H₂S浓度和CO₂浓度的偏好函数。

在系统工艺参数优化计算中，考虑设计者对不同参数具有不同的喜好程度，利用物理规划构建系统偏好函数。设定产液量最优值y_1best，设定值y_best周围某一邻域范围[y_best-Δy,y_best+Δy]内波动为非常满意(HD)，且在[y_best-Δy-Δy₁,y_best-Δy],[y_best+Δy,y_best+Δy+Δy₁]内为满意(D)，依次得到可接受(T)，不满意(U)和非常不满意(HU)，对应的偏好值区间用[0,2]，[2,4]，[4,6]，[6,8]，[8,10]表示。

构建偏好值与实际量的函数关系perf_c(H₂S)和perf_c(CO₂)，且构建偏好值与实际量的函数关系perf_c(H₂S)和perf_c(CO₂)，且与偏好程度数值区间对应的H₂S浓度以及CO₂浓度区间值如下表3所示：

表3 偏好程度数值区间与变量实际值对应区间表

根据表3的区间边界值的对应关系，可以求得H₂S浓度、CO₂浓度的偏好函数表达式如式(28)、式(29)所示。其各自偏好函数图像如图9a-图9b所示。

步骤S8：利用SPEA-II算法对H₂S浓度的偏好函数和CO₂浓度的偏好函数进行多目标极值寻优优化，获得符合生产实际的决策变量。

优化的具体过程为：

步骤S81：通过决策变量个体P＝[x₁ L x₁₀]的非支配性比较寻找最佳个体；

其中个体非支配性由适应度函数值以及个体拥挤度共同决定。

多目标优化适应度函数为：

其中，是由ST-UPFNN算法构建的油田机采油过程模型：

步骤S82：设置决策变量x₁ L x₁₀的上限和下限，其中，x_1min≤x₁≤x_1max，x_2min≤x₂≤x_2max，x_3min≤x₃≤x_3max，x_4min≤x₄≤x_4max，x_5min≤x₅≤x_5max，x_6min≤x₆≤x_6max，x_7min≤x₇≤x_7max，x_8min≤x₈≤x_8max，x_9min≤x₉≤x_9max，x_10min≤x₁₀≤x_10max。

表4 决策变量最值表

步骤S83：利用决策变量的K构建父代种群P，设置父代种群P的样本数量K＝50，设置最大遗传代数GEN＝100。

其中，K为父代种群P中的个体的数量；L为初始化的父代种群P的样本数量L＝50；GEN为最大遗传代数，GEN＝100。

步骤S84：利用决策参量(x₁,x₂)构建父代种群P的个体P_k。

P_k＝[x_1k,x_2k,...,x_10k] (33)

其中，初始化父代种群P，令初始化父代种群P为第一代父代种群P¹，设置第一代父代种群P¹的样本数量为50，设置第一代精英种群设置第一代精英种群中个体的数量为K。

步骤S85：进行第一次遗传迭代(GEN＝1)，得到第二代精英种群A²和第二代父代种群P²。

步骤S85具体包括如下步骤：

步骤S851：将第一代父代种群P¹与第一代精英种群A¹组合成种群R_t：R_t＝P¹U A¹；

定义变量R(i)为种群R_t中个体R_t(i)的强度；其中，R_t(i)为种群R_t中的第i个个体；

将种群R_t中的个体R_t(i)与环境状态变量平均值组建输入样本利用式(30)计算样本X_i的原始适应度函数值，并作为个体R_t(i)的原始适应度函数值，以求取受个体支配的解的个数S(i)，根据受个体支配的解的个数S(i)通过式(35)获得个体R_t(i)的强度值R(i)。

S(i)＝|{j|x^j∈P_k+A¹,xⁱf x^j}| (34)

步骤S852：利用种群R_t中个体R_t(i)与种群R_t中第b个邻近个体R_t(i)的距离值通过式(36)计算个体R_t(i)的密度函数D(c)。

步骤S853：根据种群R_t中个体R_t(i)的强度R(i)与种群R_t中个体R_t(i)的密度值D(c)通过式(38)获得种群R_t中个体R_t(i)的适应值F(i)。

F(i)＝R(i)+D(i) (38)

步骤S854：在保持精英种群A²中个体的数量为的情况下，将种群R_t中所有的非支配个体全部放入精英种群A²中；其中，如果精英种群A²中个体的数量小于比较种群R_t中剩余个体的适应值，选取适应值小的个体放入精英种群A²中，直到精英种群A²中个体的数量为为止；如果精英种群A²中个体的数量大于则在精英种群A²中比较个体的密度值D(c)，将密度值大的个体剔除，直到精英种群A²中个体的数量为为止。

步骤S855：将精英种群A²中的个体放入交配池中进行遗传操作获得第二代父代种群P²。

步骤S86：将第二代父代种群P²与第二代精英解种群A²组合，且GEN＝GEN+1。

步骤S87：循环100次步骤S85～步骤S86，得到精英种群A^GEN，将A^GEN的个体作为优化结果输出。

A^GEN中的个体为优化后的工艺参数x₁-x₁₀。

得到Pareto解集，H₂S浓度的偏好值与CO₂浓度的偏好值的Pareto解集如图10所示H₂S浓度的实际值与CO₂浓度的实际值的Pareto解集如图11所示。

步骤S9：将优化后的决策变量带入权值阈值更新后的神经网络模型，计算优化后的决策变量的系统性能的平均值，与实际样本的系统性能的平均值进行比较，如果优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值，利用优化后的决策变量对实际生产进行指导；否则重复上述步骤S1-S8，直至优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值为止

优化后的决策变量即为优化后的工艺参数x₁-x₁₀。

对比结果如表4所示：

表4 优化前后结果对标表

	H2S浓度(ppm)	CO2浓度(％)
			优化前	0.06	1.16
优化后	2.96	2.3

根据优化前后H₂S浓度和CO₂浓度的实际值可以看出，优化前H₂S浓度平均值和CO₂浓度平均值均远低于国标要求的4ppm和3％；经过优化之后的H₂S浓度平均值和CO₂浓度平均值既满足国标要求又提高产量降低能耗。可有效的提升原料气净化过程产出率。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (6)

1.一种高含硫天然气净化工艺动态演化建模与节能优化方法，包括：

步骤S1：选择影响脱硫效率的工艺参数和脱硫单元的性能指标；其中，所述工艺参数包括进入尾气吸收塔贫的胺液流量x₁、进入二级吸收塔的贫胺液流量x₂、原料气处理量x₃、尾气单元返回脱硫单元的半富胺液流量x₄、一级吸收塔胺液入塔温度x₅、二级吸收塔胺液入塔温度x₆、闪蒸罐压力x₇、一个重沸器的蒸汽消耗量x₈、另一个重沸器的蒸汽消耗量x₉和蒸汽预热器的蒸汽消耗量x₁₀；所述脱硫单元的性能指标包括净化气中H₂S的浓度和CO₂的浓度

步骤S2：采集预设时间的所述工艺参数和所述性能指标的数据，剔除误差样本后形成样本集[X,Y]；

步骤S3：对样本集[X,Y]进行归一化，形成归一化样本集取所述归一化样本集中前80％的样本作为训练样本，而剩余的20％样本作为测试样本；

步骤S4：基于所述训练样本训练构建神经网络模型并根据模型确定初始状态变量X，将所述训练样本中的作为所述神经网络模型的输入，将所述训练样本中的作为所述神经网络模型的输出；

其中，所述神经网络模型为：

$Y=\underset{j=1}{\overset{s_2}{\Σ}}\left(f\right(\underset{i=1}{\overset{S_1}{\Σ}}{w}_{ik}^1{I}_k+b_i^1\left)\right)\·w_{kj}^2+b_j^2---\left(1\right)$

其中，I_k为所述训练样本的矢量样本值，并作为所述神经网络模型的输入，为网络输入层到隐含层的神经元的连接权值，为网络输入层到所述隐含层的神经元的阈值，为所述隐含层到网络输出层的神经元的连接权值，为所述隐含层到所述网络输出层的神经元的阈值，其中，i＝1，2…S₀；j＝1，2…S₁；k＝1，2…S₂；S₀为所述网络输入层的神经元的数量，S₁为所述网络隐含层的神经元的数量，S₂为所述网络输出层的神经元的数量；

所述初始状态变量为：

$X={\left[\begin{array}{cccccccc}{w}_{11}^{1}L& w_{s_0{s}_1}^1& b_1^{1}L& b_{s_1}^1& w_{11}^{2}L& w_{s_1{s}_2}^2& b_1^{2}L& b_{s_2}^2\end{array}\right]}^T---\left(2\right)$

步骤S5：利用ST-UPFNN算法估计所述神经网络模型的最优状态变量；

步骤S6：将所述最优状态变量作为所述神经网络模型的和对式(1)进行更新，获得权值阈值更新后的神经网络模型；

步骤S7：分别构建H₂S浓度的偏好函数和CO₂浓度的偏好函数；

步骤S8：利用SPEA-II算法对H₂S浓度的偏好函数和CO₂浓度的偏好函数进行多目标极值寻优优化，获得符合生产实际的决策变量；

步骤S9：将优化后的决策变量带入权值阈值更新后的神经网络模型，计算优化后的决策变量的系统性能的平均值，与实际样本的系统性能的平均值进行比较，如果优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值，利用优化后的决策变量对实际生产进行指导；否则重复上述步骤S1-S8，直至优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值为止。

2.如权利要求1所述的高含硫天然气净化工艺动态演化建模与节能优化方法，其中，所述步骤S5包括：

步骤S51：针对粒子滤波器设置粒子的数目N，对进行抽样，得到初始粒子集并将所述初始粒子集中的每个粒子的权值均设为1/N；

其中，表示以x₀为均值、P₀为方差的正态分布采样；

步骤S52：在获取(k+1)时刻的观测变量值后，利用ST-UKFNN算法对每个粒子进行(k+1)时刻的状态估计，得到最优状态估计值和协方差

步骤S53：将所述最优状态估计值和所述初始化后的协方差作为粒子的重要性密度函数进行抽样，得到新粒子由所有新粒子组成的粒子集中的每个新粒子的正态分布概率密度值如下：

$X_{k+1}^{\′j}={\hat{X}}_{k+1}^j+rand---\left(3\right)$

其中，rand为正态分布随机误差；

正态分布密度函数：

$p\left(X_{k+1}^{\′j}\right|{\hat{X}}_{k+1}^j,Y_{k+1})=N(X_{k+1}^{\′j},{\hat{X}}_{k+1}^j,P_{k+1}^j)---\left(5\right)$

$p\left(Y_{k+1}\right|X_{k+1}^{\′j})=N(Y_{ek,k+1}^j,h(X_{k+1}^{\′j},I_{k+1}),1)---\left(6\right)$

$p\left(X_{k+1}^{\′j}\right|{\hat{X}}_{k+1}^j)=N(X_{k+1}^{\′j},{\hat{X}}_{k+1}^j,1)---\left(7\right)$

其中，N为每个新粒子的正态分布，x、μ、σ分别为正态分布的三个变量，p为每个新粒子的条件概率；

步骤S54：对新粒子的权值进行更新，并进行归一化处理；其中，

权值更新公式为：

权值归一化公式为：

步骤S55：根据新粒子归一化处理后的权值和重采样策略对粒子集进行重采样，获取新粒子集并求取新粒子集中每个新粒子的状态估计值

设变量u，令取u₁∈(0,1)

$X_{k+1}^j=\left\{\begin{array}{cc}{X}_{k+1}^{\&prime;c}& u\left(j\right)<\underset{i=1}{\overset{c}{\mathrm{\&Sigma;}}}{\mathrm{\&omega;}}_{k+1}^i,c=1,2L,N,\\ X_{k+1}^{\&prime;j},& else\end{array}\right.,j=1,2L,N---\left(10\right)$

${\hat{X}}_{k+1}=\left(\underset{j=1}{\overset{N}{\&Sigma;}}{X}_{k+1}^j\right)/N---\left(11\right)$

步骤S56：以粒子的数目N作为循环次数循环步骤S51-步骤S55的计算过程，将最后一次估计得到的最优状态变量作为利用ST-UPFNN算法估计所述神经网络模型的最优状态变量；其中，将新粒子的状态估计值作为本时刻的最优估计付给进行下一时刻的状态估计。

3.如权利要求2所述的高含硫天然气净化工艺动态演化建模与节能优化方法，其中，利用ST-UKFNN算法对每个粒子进行状态估计的过程，包括：

步骤S521：对所述初始状态变量X进行Sigma采样，获得2n+1个采样点，初始化控制2n+1个采样点的分布状态参数α、待选参数κ，以及非负权系数β，对所述初始状态变量X的Sigma采样如下：

$\left\{\begin{array}{c}\begin{array}{cc}{\hat{X}}_{k|k}^{\left(i\right)}=X_k+\sqrt{(n+\mathrm{\&lambda;})p_k}& i=1:n\\ {\hat{X}}_{k|k}^{\left(i\right)}=X_k-\sqrt{(n+\mathrm{\&lambda;})p_k}& i=n+1:2n\end{array}\\ \mathrm{\&lambda;}=a^2\left(n+\mathrm{\&kappa;}\right)-n\end{array}\right.---\left(12\right)$

步骤S522：计算每个采样点的权重，每个采样点的权重如下：

$\left\{\begin{array}{c}{W}_m^{\left(0\right)}=\mathrm{\&lambda;}/(n+\mathrm{\&lambda;})\\ W_c^{\left(0\right)}=\mathrm{\&lambda;}/(n+\mathrm{\&lambda;})+(1-{\mathrm{\&alpha;}}^2+\mathrm{\&beta;})\\ \begin{array}{cc}{W}_m^{\left(i\right)}=W_c^{\left(i\right)}=\mathrm{\&lambda;}/(2\&times;(n+\mathrm{\&lambda;}\left)\right)& i=1:2n\end{array}\end{array}\right.---\left(13\right)$

其中，W_c为计算状态变量的协方差的权重，W_m为计算状态估计和观测预测时的权重，是的第一列，是的第一列；

步骤S523：通过离散时间非线性系统的状态方程将每个采样点的k时刻的最优状态变量的状态估计变换为(k+1)时刻的状态变量的状态估计并通过合并(k+1)时刻的状态估计的向量，获得(k+1)时刻的状态变量的状态先验估计和协方差P_k+1|k；其中，

所述状态估计为：

$X_{k+1|k}^{\left(i\right)}=F\left({\hat{X}}_{k|k}^{\left(i\right)}\right)+w_k---\left(14\right)$

其中，表示k时刻的最优估计，w_k为过程噪声，其协方差矩阵Q_k为cov(w_k,w_j)＝Q_kδ_kj，

所述状态先验估计为：

${\hat{X}}_{k+1|k}=\underset{i=0}{\overset{2n}{\mathrm{\&Sigma;}}}{W}_m^{\left(i\right)}\&CenterDot;X_{k+1|k}^{\left(i\right)}---\left(15\right)$

所述状态变量的协方差P_k+1|k为：

$P_{k+1|k}=\underset{i=0}{\overset{2n}{\mathrm{\&Sigma;}}}{W}_c^{\left(i\right)}\&CenterDot;\&lsqb;X_{k+1|k}^{\left(i\right)}-{\hat{X}}_{k+1|k}^{\left(i\right)}\&rsqb;{\&lsqb;X_{k+1|k}^{\left(i\right)}-{\hat{X}}_{k+1|k}^{\left(i\right)}\&rsqb;}^T+Q_k---\left(16\right)$

步骤S524：通过离散时间非线性系统的观测方程将获得的(k+1)时刻的状态变量的状态估计与(k+1)时刻的观测预测建立联系：

$Y_{k+1|k}^{\left(i\right)}=g\left(\underset{k=1}{\overset{s_2}{\&Sigma;}}\right(f\left(\underset{j=1}{\overset{S_1}{\&Sigma;}}{w}_{ji}^1{X}_{k+1|k}^{\left(i\right)}+b_j^1\right))\&CenterDot;w_{kj}^2+b_k^2)+v_k---\left(17\right)$

其中，ν_k为观测噪声，其协方差矩阵R_k为cov(v_k,v_j)＝R_kδ_kj，

步骤S525：通过估计(k+1)时刻的观测预测的向量，获得(k+1)时刻的先验观测预测并根据先验观测预测估计(k+1)时刻观测预测的协方差

(k+1)时刻的观测预测为：

${\hat{Y}}_{k+1|k}=\underset{i=0}{\overset{2n}{\&Sigma;}}{W}_m^{\left(i\right)}\&CenterDot;Y_{k+1|k}^{\left(i\right)}---\left(18\right)$

(k+1)时刻观测预测的协方差为：

$P_{y_{k+1}}={\mathrm{\&lambda;}}_{k+1}\underset{i=0}{\overset{2n}{\&Sigma;}}{W}_c^{\left(i\right)}\&CenterDot;(Y_{k+1|k}^{\left(i\right)}-{\hat{Y}}_{k+1|k}){(Y_{k+1|k}^{\left(i\right)}-{\hat{Y}}_{k+1|k})}^T+R_{k+1}---\left(19\right)$

其中，

N_k+1＝V_k+1-βR_k+1 (21)

其中，β为弱化因子，β≥1；

$M_{k+1}=\underset{i=0}{\overset{2n}{\mathrm{\&Sigma;}}}{W}_c^{\left(i\right)}\&CenterDot;(Y_{k+1|k}^{\left(i\right)}-{\hat{Y}}_{k+1|k}){(Y_{k+1|k}^{\left(i\right)}-{\hat{Y}}_{k+1|k})}^T+R_{k+1}---\left(22\right)$

$V_{k+1}=\left\{\begin{array}{cc}{e}_{k+1}{e}_{k+1}^T& k=0\\ \frac{{\mathrm{\&rho;V}}_k+e_{k+1}{e}_{k+1}^T}{1+\mathrm{\&rho;}}& k\&GreaterEqual;1\end{array}\right.,\mathrm{\&rho;}\&Element;(0,1)---\left(23\right)$

$e_{k+1}=Y_{k+1}-\underset{i=0}{\overset{2n}{\&Sigma;}}{W}_m^{\left(i\right)}\&CenterDot;h({\hat{X}}_{k+1|k}^{\left(i\right)},I_k)---\left(24\right)$

步骤S526：计算(k+1)时刻的状态变量的状态先验估计与(k+1)时刻的先验观测预测之间的协方差所述协方差为：

$P_{x_{k+1}{y}_{k+1}}={\mathrm{\&lambda;}}_{k+1}\underset{i=0}{\overset{2n}{\&Sigma;}}{W}_c^{\left(i\right)}\&CenterDot;(X_{k+1|k}^{\left(i\right)}-{\hat{X}}_{k+1|k}){(Y_{k+1|k}^{\left(i\right)}-{\hat{Y}}_{k+1|k})}^T---\left(25\right)$

步骤S527：通过建立协方差和协方差的关系，更新(k+1)时刻的状态变量的状态估计和协方差，分别获得(k+1)时刻的最优状态估计值和协方差

4.如权利要求3所述的高含硫天然气净化工艺动态演化建模与节能优化方法，其中，S527建立的协方差和协方差的关系为：

$K_{k+1}=\frac{P_{x_{k+1}{y}_{k+1}}}{P_{y_{k+1}}}---\left(26\right)$

其中，K_k+1为增益矩阵；以及，

更新后的(k+1)时刻的状态变量的最优状态估计为：

${\hat{X}}_{k+1}={\hat{X}}_{k+1|k}+K_{k+1}(Y_{k+1}-{\hat{Y}}_{k+1|k})---\left(27\right)$

更新后的(k+1)时刻的状态变量的协方差P_k+1为：

$P_{k+1}=P_{k+1|k}-K_{k+1}{\left(P_{x_{k+1}{y}_{k+1}}\right)}^T---\left(28\right)$

将更新后的(k+1)时刻的最优状态变量的状态估计和协方差P_k+1分别作为(k+1)时刻的最优状态估计值和协方差

5.如权利要求1所述的高含硫天然气净化工艺动态演化建模与节能优化方法，其中，步骤S8包括：

步骤S81：通过决策变量个体P＝[x₁ L x₁₀]的非支配性比较寻找最佳个体；其中个体非支配性由适应度函数值以及个体拥挤度共同决定，多目标优化适应度函数为：

$objFun\left(X\right)=\left[\begin{array}{cc}{\mathrm{perf}}_c\left(g^{-1}\right({\hat{y}}_{H_{2}S}\left(f\left(X\right)\right)\left)\right)& g^{-1}\left({\hat{y}}_{{\mathrm{CO}}_2}\right(f\left(X\right)\left)\right)\end{array}\right]---\left(29\right)$

其中，是由ST-UPFNN算法构建的油田机采油过程模型：

$\hat{Y}\left(X\right)=g^{-1}\left(\tilde{Y}\right(\tilde{X}\left)\right)=g^{-1}\left(\tilde{Y}\right(f\left(X\right)\left)\right)={\&lsqb;{\hat{Y}}_{H_{2}S},{\hat{Y}}_{{\mathrm{CO}}_2}\&rsqb;}^T---\left(30\right)$

步骤S82：设置决策变量x₁ L x₁₀的上限和下限，其中，x_1min≤x₁≤x_1max，x_2min≤x₂≤x_2max，x_3min≤x₃≤x_3max，x_4min≤x₄≤x_4max，x_5min≤x₅≤x_5max，x_6min≤x₆≤x_6max，x_7min≤x₇≤x_7max，x_8min≤x₈≤x_8max，x_9min≤x₉≤x_9max，x_10min≤x₁₀≤x_10max；

步骤S83：利用决策变量的K构建父代种群P，设置父代种群P的样本数量K＝50，设置最大遗传代数GEN＝100；

$P=\{(x_1^1,x_2^1,...x_{10}^1),(x_1^2,x_2^2,...x_{10}^2),L(x_1^K,x_2^K,...x_{10}^K)\}---\left(31\right)$

其中，K为父代种群P中的个体的数量；L为初始化的父代种群P的样本数量L＝50；GEN为最大遗传代数，GEN＝100；

步骤S84：利用决策参量(x₁,x₂,…x₁₀)构建父代种群P的个体P_k；

P_k＝[x_1k,x_2k,...,x_10k] (32)

其中，初始化父代种群P，令初始化父代种群P为第一代父代种群P¹，设置第一代父代种群P¹的样本数量为50，设置第一代精英种群设置第一代精英种群中个体的数量为

步骤S85：进行第一次遗传迭代(GEN＝1)，得到第二代精英种群A²和第二代父代种群P²；

步骤S86：将第二代父代种群P²与第二代精英解种群A²组合，且GEN＝GEN+1；

步骤S87：循环100次步骤S85～步骤S86，得到精英种群A_GEN，将A^GEN中的个体作为优化结果输出，A^GEN中的个体为优化后的工艺参数x₁-x₁₀。

6.如权利要求5所述的高含硫天然气净化工艺动态演化建模与节能优化方法，其中，步骤S85包括：

步骤S851：将第一代父代种群P¹与第一代精英种群A¹组合成种群R_t：R_t＝P¹U A¹；

定义变量R(i)为种群R_t中个体R_t(i)的强度；其中，R_t(i)为种群R_t中的第i个个体；

将种群R_t中的个体R_t(i)与环境状态变量平均值组建输入样本利用式(29)计算样本X_i的原始适应度函数值，并作为个体R_t(i)的原始适应度函数值，以求取受个体支配的解的个数S(i)，根据受个体支配的解的个数S(i)通过式(34)获得个体R_t(i)的强度值R(i)；

S(i)＝|{j|x^j∈P_k+A¹,xⁱf x^j}| (33)

$R\left(i\right)=\underset{x^j\&Element;P_k+A^1,x^j{\mathrm{fx}}^i}{\&Sigma;}S\left(j\right)---\left(34\right)$

步骤S852：利用种群R_t中个体R_t(i)与种群R_t中第b个邻近个体R_t(i)的距离值通过式(35)计算个体Rt(i)的密度函数D(c)；

$D\left(c\right)=\frac{1}{{\mathrm{\&sigma;}}_c^b+2}---\left(35\right)$

$b=\sqrt{K+\stackrel{\&OverBar;}{K}}---\left(36\right)$

步骤S853：根据种群R_t中个体R_t(i)的强度R(i)与种群R_t中个体R_t(i)的密度值D(c)通过式(37)获得种群R_t中个体R_t(i)的适应值F(i)；

F(i)＝R(i)+D(i) (37)

步骤S854：在保持精英种群A²中个体的数量为的情况下，将种群R_t中所有的非支配个体全部放入精英种群A²中；其中，如果精英种群A²中个体的数量小于比较种群R_t中剩余个体的适应值，选取适应值小的个体放入精英种群A²中，直到精英种群A²中个体的数量为为止；如果精英种群A²中个体的数量大于则在精英种群A²中比较个体的密度值D(c)，将密度值大的个体剔除，直到精英种群A²中个体的数量为为止；

步骤S855：将精英种群A²中的个体放入交配池中进行遗传操作获得第二代父代种群P²。