CN108681653A - 基于动态子空间高阶累积量分析的天然气净化过程异常监测方法 - Google Patents
基于动态子空间高阶累积量分析的天然气净化过程异常监测方法 Download PDFInfo
- Publication number
- CN108681653A CN108681653A CN201810566582.8A CN201810566582A CN108681653A CN 108681653 A CN108681653 A CN 108681653A CN 201810566582 A CN201810566582 A CN 201810566582A CN 108681653 A CN108681653 A CN 108681653A
- Authority
- CN
- China
- Prior art keywords
- order cumulants
- dynamic
- model
- natural gas
- isolated component
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- C—CHEMISTRY; METALLURGY
- C10—PETROLEUM, GAS OR COKE INDUSTRIES; TECHNICAL GASES CONTAINING CARBON MONOXIDE; FUELS; LUBRICANTS; PEAT
- C10L—FUELS NOT OTHERWISE PROVIDED FOR; NATURAL GAS; SYNTHETIC NATURAL GAS OBTAINED BY PROCESSES NOT COVERED BY SUBCLASSES C10G, C10K; LIQUEFIED PETROLEUM GAS; ADDING MATERIALS TO FUELS OR FIRES TO REDUCE SMOKE OR UNDESIRABLE DEPOSITS OR TO FACILITATE SOOT REMOVAL; FIRELIGHTERS
- C10L3/00—Gaseous fuels; Natural gas; Synthetic natural gas obtained by processes not covered by subclass C10G, C10K; Liquefied petroleum gas
- C10L3/06—Natural gas; Synthetic natural gas obtained by processes not covered by C10G, C10K3/02 or C10K3/04
- C10L3/10—Working-up natural gas or synthetic natural gas
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Mathematical Physics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Organic Chemistry (AREA)
- Bioinformatics & Computational Biology (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- General Chemical & Material Sciences (AREA)
- Algebra (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Gas Separation By Absorption (AREA)
Abstract
本发明公开了一种基于动态子空间高阶累积量分析的天然气净化过程异常监测方法,该监测方法按如下步骤进行:(1)选取所监测的高含硫天然气净化脱硫生产过程中的工艺参数;(2)采集天然气净化过程中所述工艺参数的正常工况下数据;(3)对数据进行预处理;(4)构建基于动态子空间高阶累积量分析的模型,计算正常工况下的总体独立分量三阶累积量HS与总体残差三阶累积量HE,并用核密度估计获得其控制限;(5)在线获取高含硫天然气净化脱硫生产过程中所述工艺参数运行工况下数据,重新输入模型计算统计量,判断所得统计量是否超出控制限,若超限则发生异常工况,否则正常。本方法能够及时检测故障发生,有效降低漏报率和误报率。
Description
技术领域
本发明涉及高含硫天然气脱硫生产过程故障检测与诊断技术,具体为一种基于动态子空间高阶累积量分析的天然气净化过程异常监测方法。
背景技术
我国高含硫气藏资源丰富,累计探明储量逾1×1012m3,开发潜力巨大。此类气田酸性组分含量高达普通气田的数十倍,并且所含剧毒成分硫化氢 (Hydrogen sulfide,H2S)对集采设施、管线及净化装置都会造成强烈腐蚀,安全开发风险极高。一旦生产过程出现故障,导致酸性气体泄漏,将会造成严重的人员伤亡和环境污染事故。因此,高含硫天然气净化过程的实时监控对于保证***可靠运行、实现安全生产具有重要意义。
对于此类化工过程,变量统计过程监测(MSPM)方法因无需精准数学模型而在此监控领域备受关注。近几年来,MSPM及其衍生算法也有针对性的解决了关于非线性、动态、非高斯等不同特性的故障监测问题。其中高阶累积量分析(HCA)方法在独立分量分析(ICA)方法的基础上,采用独立分量和残差矩阵的三阶累积量重新构造统计指标进行检测,大大提高了非高斯过程故障的监测率。然而HCA属于静态算法,会因无法有效识别和反馈工况参数上的动态相关性而存在过程监测误报警问题。
发明内容
鉴于上述问题,本发明的目的是提供一种基于动态子空间高阶累积量分析的天然气净化过程异常监测方法,其能够及时检测故障发生,降低误报率和漏报率,消除正常工况下的误报警。
本发明技术方案如下:
一种基于动态子空间高阶累积量分析的天然气净化过程异常监测方法,其特征在于,该监测方法按如下步骤进行:
步骤1:选取所监测的高含硫天然气净化脱硫生产过程中的m个工艺参数,其中,m=26,各工艺参数分别为:x1原料气中H2S含量、x2溶液循环量、 x3贫MDEA溶液入塔压力、x4吸收塔液位、x5吸收塔压差、x6湿净化气压力、 x7湿净化气H2S含量、x8湿净化气CO2含量、x9闪蒸塔压差、x10闪蒸气流量、 x11闪蒸塔压力、x12闪蒸塔液位、x13MDEA活性炭过滤器压差、x14贫富液换热器富液出口温度、x15贫富液换热器贫液出口温度、x16再生塔液位、x17再生塔塔顶温度、x18再生塔上部压力、x19再生塔下部压力、x20重沸器进口半贫液温度、x21重沸器出口半贫液温度、x22重沸器蒸汽压力、x23重沸器蒸汽流量、 x24酸气分离器压力、x25酸气分离器液位、x26酸水回流量;
步骤2:预设采样间隔时间,采集所述工艺参数的正常工况下数据,剔除误差样本后形成训练集x∈Rm×n;
步骤3:对步骤2所得的训练样本数据x∈Rm×n进行均值中心化以及标准化的数据处理;
步骤4:基于步骤3所得的训练样本数据,构建基于动态子空间高阶累积量分析的模型,计算在正常工况下的总体独立分量三阶累积量HS与总体残差三阶累积量HE,并通过核密度估计获得HS与HE统计量的控制限;
步骤5:在线获取高含硫天然气净化脱硫生产过程中的所述工艺参数的运行工况下数据,通过步骤4构建的基于动态子空间高阶累积量分析的模型,计算当前工况下的HS与HE统计量,分析统计量是否超出正常工况下控制限,若超限则发生异常工况,否则正常。
由于采用了上述技术方案,本发明具有如下有益效果:
本发明针对HCA方法未考虑数据间动态时序相关性问题,首先,建立动态子空间模型对数据进行时序上的动态扩展,将天然气净化过程视为动态***,引入工况参数及其时间序列作为输入备选特征,进行时滞阶次分析,得到对应的动态扩展数据,再结合HCA方法验证动态子空间模型在处理净化过程数据动态相关性上的性能,本发明引入动态子空间模型,提高故障样本的处理能力,有效降低漏报率和误报率,实现对实际工业过程的在线监控,有利于保障故障安全。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。
图1为本发明流程图;
图2为吸收塔拦液故障下各变量变化趋势图;
图3为吸收塔拦液故障下统计指标图。
具体实施方式
下面结合附图对本发明的实施例作进一步详细说明。
如图1所示的一种基于动态子空间高阶累积量分析的天然气净化过程异常监测方法按如下步骤进行:
步骤1:选取所监测的高含硫天然气净化脱硫生产过程中的m个工艺参数,其中,m=26,各工艺参数分别为:x1原料气中H2S含量、x2溶液循环量、 x3贫MDEA溶液入塔压力、x4吸收塔液位、x5吸收塔压差、x6湿净化气压力、 x7湿净化气H2S含量、x8湿净化气CO2含量、x9闪蒸塔压差、x10闪蒸气流量、 x11闪蒸塔压力、x12闪蒸塔液位、x13MDEA活性炭过滤器压差、x14贫富液换热器富液出口温度、x15贫富液换热器贫液出口温度、x16再生塔液位、x17再生塔塔顶温度、x18再生塔上部压力、x19再生塔下部压力、x20重沸器进口半贫液温度、x21重沸器出口半贫液温度、x22重沸器蒸汽压力、x23重沸器蒸汽流量、 x24酸气分离器压力、x25酸气分离器液位、x26酸水回流量。26个工艺参数见下表:
表1高含硫天然气净化过程主要操作变量及技术指标
步骤2:预设采样间隔时间,采集所述工艺参数的正常工况下数据,剔除误差样本后形成训练集x∈Rm×n;
步骤3:对步骤2所得的训练样本数据x∈Rm×n行均值中心化以及标准化的数据处理;
本实施例中,以一分钟的采样间隔时间采集正常工况下数据1440组,利用3σ法去除粗大误差样本,获得有效样本1173组。取其中前500组样本作为训练集,剩余673组作为测试集,用于测试算法性能。
步骤4:基于步骤3所得的500组训练样本数据,构建基于动态子空间高阶累积量分析模型,计算在正常工况下的总体独立分量三阶累积量HS与总体残差三阶累积量HE,并通过核密度估计获得HS与HE统计量的控制限。
所述步骤4包括:
步骤4.1:对正常工况下的样本预处理数据集x∈Rm×n,进行自回归模型分析,结合平行分析法利用自相关系数来确定整个模型动态时滞阶次h,动态时滞阶次h的计算方法为:
步骤4.1.1:初始化时滞阶次h,假设h=0;
步骤4.1.2:选择原始矩阵对应u个最大主元;
步骤4.1.3:定义原始矩阵X新的关系数为:
r(h)=(l+1)m-u (10)
步骤4.1.4:将扩展后的动态矩阵x(h)带入计算扩展后新的关系数为:
步骤4.1.5:判断若rnew(h)>0,设定h=h+1,返回步骤4.1.2否则结束h=h-1。
步骤4.2:对预处理数据进行动态扩展获得动态矩阵x(h):
假设过程按多元AR(h)模型表示为:
x(k)=β1x(k-1)+β2x(k-2)+…+βhx(k-h)+εk (12)
其中,x(k)为模型变量,β1,β2,...,βk为模型回归系数,εk为模型随机误差,h 为模型阶次,用x(k),x(k-1),...,x(k-h)组成的增广数据矩阵可以反映变量间动态关系,即:
其中在序列矩阵中,指在k时刻,m维的观测矢量。
步骤4.3:对获得的动态矩阵x(h),运用独立分量分析求得分离矩阵W并重构独立元s矩阵;
x=As,s=Wx (14)
步骤4.4:将整个过程分为主导独立分量和模型残差两部分:
sd=Wdx (15)
其中,sd为根据解混矩阵W行二范数降序选择的前d个主导独立分量矩阵,相应的W的行构成的子矩阵记作Wd,与其对应的混合子矩阵记为Ad,是x 在主导独立分量sd上的投影。
步骤4.5:计算总体独立分量三阶累积量HS并用核密度估计其控制限;
对于过程独立分量空间sd,在采样i处,第p个主导独立分量sp的样本独立分量三阶累积量为:
hsp(i)=sp(i)sp(i-1)sp(i-2)
=wpx(i)wpx(i-1)wpx(i-2) (17)
其中,wp是解混矩阵Wd的第p行,p=1,2,…,d。为了监测全部样本独立分量三阶累积量,进而定义总体独立分量三阶累积量HS为:
其中,mhsp与vhsp分别为hsp的均值与标准差。用核密度估计其控制限记为lHS。
步骤4.6:计算总体残差三阶累积量HE并用核密度估计其控制限记为lHE。
对于模型残差空间e,在采样i处,模型对第q个变量的样本残差三阶累积量为:
heq(i)=eq(i)eq(i-1)eq(i-2)
=lqx(i)lqx(i-1)lqx(i-2) (19)
其中,lq是L的第q行,q=1,2,…,m。为了监测所有样本残差三阶累积量,定义了总体残差三阶累积量HE为:
其中,mheq与vheq分别为heq的均值与标准差。
步骤5:在线获取高含硫天然气净化脱硫生产过程中的所述工艺参数的运行工况下数据,通过步骤4构建的基于动态子空间高阶累积量分析的模型,计算当前工况下的HS与HE统计量,分析统计量是否超出正常工况下控制限,若超限则发生异常工况,否则正常。
在实际工业现场,相比误报警,漏报警更为重要,接下来以天然气净化过程中一个典型故障(拦液)为例验证算法对于故障工况的监测效果。异常工况检测分为离线训练和在线测试两个阶段。首先离线阶段采用1173组有效样本中前500组样本作为训练集,建立高阶累积量分析(HCA)和动态子空间高阶累积量(DHCA)模型,计算置信度为α=99%下的控制限。生成剩余 673组拦液故障下的样本集作为测试集,用于测试算法性能。
本发明通过几组测试得到如下的技术效果:
图2为26个变量在拦液故障下的变化趋势图。故障于测试集的第301个样本点开始引入,其中变量5(吸收塔差压)发生了阶跃;变量10(闪蒸气流量)有明显的波动;变量12(闪蒸塔液位)和变量16(再生塔液位)均逐渐下降;其他各变量维持正常范围。
将测试数据输入模型计算HCA和DHCA算法统计量,并绘制监测统计图,结果如图3所示,图3a为HCA算法,图3bwei DHCA算法。从图3中可以看出,静态HCA模型在故障发生的时间段内,残差空间指标HE不能全面的监测出故障的发生,出现了大面积的漏检情况。而采用了动态矩阵的 DHCA方法,对于故障动态特性有了更好的监测性能,减少了故障数据的漏检率。而对于正常样本,HCA的独立分量空间指标也明显不如DHCA的监测效果好。
表2展示了两种算法检测率的比较。两种算法在独立空间指标HS上都取得了较低漏报率(MAR)。这表明,它们独立空间指标对故障样本具有较好的检测灵敏度。同样的,在残差空间监测指标HE上两种算法也都取得了较好的误检率(FAR),表明残差空间指标对正常样本具有较好的检测灵敏度。尽管如此, DHCA的误简率和漏检率均低于HCA,表明在考虑动态相关性后,算法的监测性能确有提高。另一方面,与HCA相比,DHCA的HS FAR降低了36%。更值得一提的是,DHCA的HE MAR值显现出了比HCA优越很多的效果。综合表明,引入动态子空间模型可以提高故障样本的处理能力,保证故障的安全。
表2 HCA和DHCA故障检测效果对比
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解;其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的范围,其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (3)
1.一种基于动态子空间高阶累积量分析的天然气净化过程异常监测方法,其特征在于,该监测方法按如下步骤进行:
步骤1:选取所监测的高含硫天然气净化脱硫生产过程中的m个工艺参数,其中,m=26,各工艺参数分别为:x1原料气中H2S含量、x2溶液循环量、x3贫MDEA溶液入塔压力、x4吸收塔液位、x5吸收塔压差、x6湿净化气压力、x7湿净化气H2S含量、x8湿净化气CO2含量、x9闪蒸塔压差、x10闪蒸气流量、x11闪蒸塔压力、x12闪蒸塔液位、x13MDEA活性炭过滤器压差、x14贫富液换热器富液出口温度、x15贫富液换热器贫液出口温度、x16再生塔液位、x17再生塔塔顶温度、x18再生塔上部压力、x19再生塔下部压力、x20重沸器进口半贫液温度、x21重沸器出口半贫液温度、x22重沸器蒸汽压力、x23重沸器蒸汽流量、x24酸气分离器压力、x25酸气分离器液位、x26酸水回流量;
步骤2:预设采样间隔时间,采集所述工艺参数的正常工况下数据,剔除误差样本后形成训练集x∈Rm×n;
步骤3:对步骤2所得的训练样本数据x∈Rm×n进行均值中心化以及标准化的数据处理;
步骤4:基于步骤3所得的训练样本数据,构建基于动态子空间高阶累积量分析的模型,计算在正常工况下的总体独立分量三阶累积量HS与总体残差三阶累积量HE,并通过核密度估计获得HS与HE统计量的控制限;
步骤5:在线获取高含硫天然气净化脱硫生产过程中的所述工艺参数的运行工况下数据,通过步骤4构建的基于动态子空间高阶累积量分析的模型,计算当前工况下的HS与HE统计量,分析统计量是否超出正常工况下控制限,若超限则发生异常工况,否则正常。
2.根据权利要求1所述的基于动态子空间高阶累积量分析的天然气净化过程异常监测方法,其特征在于,所述步骤4包括:
步骤4.1:对正常工况下的样本预处理数据集x∈Rm×n,进行自回归模型分析,结合平行分析法利用自相关系数确定整个模型动态时滞阶次h;
步骤4.2:对预处理数据进行动态扩展获得动态矩阵x(h):
过程按多元AR(h)模型表示为:
x(k)=β1x(k-1)+β2x(k-2)+…+βhx(k-h)+εk (1)
其中,x(k)为模型变量,β1,β2,...,βk为模型回归系数,εk为模型随机误差,h为模型阶次,用x(k),x(k-1),...,x(k-h)组成的增广数据矩阵反映变量间动态关系:
其中在序列矩阵中,是指在k时刻,m维的观测矢量;
步骤4.3:对获得的动态矩阵x(h),运用独立分量分析求得分离矩阵W并重构独立元s矩阵;
x=As,s=Wx (3)
步骤4.4:将整个过程分为主导独立分量和模型残差两部分:
sd=Wdx (4)
其中,sd为根据解混矩阵W行二范数降序选择的前d个主导独立分量矩阵,相应的W的行构成的子矩阵记作Wd,与其对应的混合子矩阵记为Ad,是x在主导独立分量sd上的投影;
步骤4.5:计算总体独立分量三阶累积量HS并用核密度估计其控制限;
对于过程独立分量空间sd,在采样i处,第p个主导独立分量sp的样本独立分量三阶累积量为:
hsp(i)=sp(i)sp(i-1)sp(i-2)
=wpx(i)wpx(i-1)wpx(i-2) (6)
其中,wp是解混矩阵Wd的第p行,p=1,2,…,d,为了监测全部样本独立分量三阶累积量,进而定义总体独立分量三阶累积量HS为:
其中,mhsp与vhsp分别为hsp的均值与标准差,用核密度估计其控制限记为lHS。
步骤4.6:计算总体残差三阶累积量HE并用核密度估计其控制限;
对于模型残差空间e,在采样i处,模型对第q个变量的样本残差三阶累积量为:
heq(i)=eq(i)eq(i-1)eq(i-2)
=lqx(i)lqx(i-1)lqx(i-2) (8)
其中,lq是L的第q行,q=1,2,…,m。为了监测所有样本残差三阶累积量,定义了总体残差三阶累积量HE为:
其中,mheq与vheq分别为heq的均值与标准差。
3.根据权利要求2所述的基于动态子空间高阶累积量分析的天然气净化过程异常监测方法,其特征在于,所述步骤4.1计算时滞后阶次h的方法为:
步骤4.1.1:初始化时滞阶次h,假设h=0;
步骤4.1.2:选择原始矩阵对应u个最大主元;
步骤4.1.3:定义原始矩阵X新的关系数为:
r(h)=(l+1)m-u (10)
步骤4.1.4:将扩展后的动态矩阵x(h)带入计算扩展后新的关系数为:
步骤4.1.5:判断若rnew(h)>0,设定h=h+1,返回步骤4.1.2,否则结束h=h-1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810566582.8A CN108681653A (zh) | 2018-06-05 | 2018-06-05 | 基于动态子空间高阶累积量分析的天然气净化过程异常监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810566582.8A CN108681653A (zh) | 2018-06-05 | 2018-06-05 | 基于动态子空间高阶累积量分析的天然气净化过程异常监测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108681653A true CN108681653A (zh) | 2018-10-19 |
Family
ID=63810031
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810566582.8A Pending CN108681653A (zh) | 2018-06-05 | 2018-06-05 | 基于动态子空间高阶累积量分析的天然气净化过程异常监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108681653A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112560924A (zh) * | 2020-12-07 | 2021-03-26 | 华东理工大学 | 基于动态内部慢特征分析的丙烯精馏塔状态监控方法 |
CN114437844A (zh) * | 2020-11-03 | 2022-05-06 | 中国石油化工股份有限公司 | 一种天然气选择性脱氮工艺的参数计算机自动优化方法 |
CN114632409A (zh) * | 2022-02-23 | 2022-06-17 | 国能龙源环保有限公司 | 适用不同浓度的羰基硫和硫化氢脱除***及其使用方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104635724A (zh) * | 2014-12-25 | 2015-05-20 | 重庆科技学院 | 基于动态核独立分量分析的天然气净化过程异常检测方法 |
CN104656635A (zh) * | 2014-12-31 | 2015-05-27 | 重庆科技学院 | 非高斯动态高含硫天然气净化过程异常检测与诊断方法 |
-
2018
- 2018-06-05 CN CN201810566582.8A patent/CN108681653A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104635724A (zh) * | 2014-12-25 | 2015-05-20 | 重庆科技学院 | 基于动态核独立分量分析的天然气净化过程异常检测方法 |
CN104656635A (zh) * | 2014-12-31 | 2015-05-27 | 重庆科技学院 | 非高斯动态高含硫天然气净化过程异常检测与诊断方法 |
Non-Patent Citations (1)
Title |
---|
YOUQING WANG,ET AL.: "Online Monitoring of Multivariate Processes Using Higher-Order", 《INDUSTRIAL & ENGINEERING CHEMISTRT REASEARCH》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114437844A (zh) * | 2020-11-03 | 2022-05-06 | 中国石油化工股份有限公司 | 一种天然气选择性脱氮工艺的参数计算机自动优化方法 |
CN112560924A (zh) * | 2020-12-07 | 2021-03-26 | 华东理工大学 | 基于动态内部慢特征分析的丙烯精馏塔状态监控方法 |
CN112560924B (zh) * | 2020-12-07 | 2024-01-19 | 华东理工大学 | 基于动态内部慢特征分析的丙烯精馏塔状态监控方法 |
CN114632409A (zh) * | 2022-02-23 | 2022-06-17 | 国能龙源环保有限公司 | 适用不同浓度的羰基硫和硫化氢脱除***及其使用方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109522600B (zh) | 基于组合深度神经网络的复杂装备剩余使用寿命预测方法 | |
Du et al. | Fault detection and diagnosis using empirical mode decomposition based principal component analysis | |
CN110207997B (zh) | 基于卷积自编码器的液体火箭发动机故障检测方法 | |
CN105259895B (zh) | 一种工业过程微小故障的检测和分离方法及其监测*** | |
Maurya et al. | A signed directed graph and qualitative trend analysis-based framework for incipient fault diagnosis | |
CN108681653A (zh) | 基于动态子空间高阶累积量分析的天然气净化过程异常监测方法 | |
Tong et al. | Decentralized monitoring of dynamic processes based on dynamic feature selection and informative fault pattern dissimilarity | |
CN103713628B (zh) | 基于符号有向图和数据重构的故障诊断方法 | |
CN104035431B (zh) | 用于非线性过程监控的核函数参数的获取方法和*** | |
CH698898B1 (de) | System und Verfahren zur Zustandsüberwachung eines Anlagensystems. | |
Li et al. | Fault identification in PCA method during sensor condition monitoring in a nuclear power plant | |
CN111780800A (zh) | 一种传感器故障的监测、隔离以及重构方法和*** | |
CN111474919A (zh) | 基于aann网络组的航空发动机控制***传感器故障诊断方法 | |
CN109871000A (zh) | 一种数据驱动的闭环工业过程传感器故障诊断方法 | |
Li et al. | Dynamic time warping based causality analysis for root-cause diagnosis of nonstationary fault processes | |
CN104503436B (zh) | 一种基于随机投影和k近邻的快速故障检测方法 | |
CN113011248A (zh) | 一种基于lstm的核动力装置故障诊断方法 | |
CN103853144B (zh) | 基于采油生产数据的现场传感器故障检测方法 | |
CN107918381A (zh) | 一种基于组合核函数的类均值核主元故障诊断方法 | |
CN107272648B (zh) | 一种工业过程微小故障的分离方法 | |
CN103995985B (zh) | 基于Daubechies小波变换和弹性网的故障检测方法 | |
CN112052496B (zh) | 一种基于var模型的高拱坝谷幅变形影响因素分析***的操作方法 | |
CN206833239U (zh) | 一种基于数据驱动的火电厂控制***故障检测*** | |
CN105279553B (zh) | 一种高加给水***故障程度识别方法 | |
Uren et al. | An integrated approach to sensor FDI and signal reconstruction in HTGRs–Part I: Theoretical framework |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181019 |