CN100507509C - 基于主成分分析和支持向量机的油气水多相流流型识别方法 - Google Patents
基于主成分分析和支持向量机的油气水多相流流型识别方法 Download PDFInfo
- Publication number
- CN100507509C CN100507509C CNB2006100170900A CN200610017090A CN100507509C CN 100507509 C CN100507509 C CN 100507509C CN B2006100170900 A CNB2006100170900 A CN B2006100170900A CN 200610017090 A CN200610017090 A CN 200610017090A CN 100507509 C CN100507509 C CN 100507509C
- Authority
- CN
- China
- Prior art keywords
- vector machine
- pressure
- signal
- flow type
- support vector
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Complex Calculations (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种基于主成分分析和支持向量机的油气水多相流流型识别方法,其特点是:利用水平管路、取压环、取压管、压力、差压变送器、数据采集卡和计算机对不同流型的流动参数信号进行采集;将采集到的压力和不同取压间距的压差波动信号进行小波包去除噪声预处理,再对上述信号进行经验模式分解,得到的信号固有模态函数,将这些固有模态函数构成一个特征矩阵X,利用主成分分析得到流型的特征向量;再将这些特征作为支持向量机的输入样本,应用支持向量机完成从特征空间到流型空间的映射,最终实现流型识别。所采用的特征向量融合了多个传感器的信息,能够全面反映流型的信息,所应用的支持向量机同神经网络分类器相比,具有学习速度快、分类能力强等特点。
Description
技术领域
本发明涉及多相流测量技术领域,尤其涉及一种基于主成分分析和支持向量机的油气水多相流流型识别方法。
背景技术
传统的流型识别方法有两大类:一类是采用实验方法做出流型图;另一类是根据对流型转变机理的分析得到的转变准则关系式,利用现场的流动参数来确定流型。由于在生产实际中流动参数往往是需要解决的问题,因此传统的方法无法得到广泛应用。
流型的现代测量方法从工作原理上可分为直接测量和间接测量方法。常用的直接测量方法有目测法、高速摄影法等,该方法具有一定的主观性,对于多个观察者可能得出不同的结论。间接测量方法是利用测量的仪器设备测量的压力、压差和含气率等波动信号,通过对信号进行分析、提取特征,结合神经网络等模式识别技术来识别流型。该方法有两个关键问题:一是流型特征向量的提取;二是网络模型的选取,其中特征提取尤为关键。现有的方法主要是利用单一的压力或压差波动信号来提取信号的单一特征,导致提取的特征无法全面反映流型的信息,致使流型识别率不高。
发明内容
本发明的目的在于克服上述现有技术的缺陷,提出一种流型识别准确性、可靠性高,通用性强的基于主成分分析和支持向量机的油气水多相流流型识别方法。
本发明的目的是由以下技术方案来实现的:一种基于主成分分析和支持向量机的油气水多相流流型识别方法,其特征在于,它包括以下步骤:
(1)对不同流型的流动参数信号进行采集:利用取压间距分别为5倍、10倍和15倍管径D;采样频率为256Hz;采样点数为1024点;1个压力变送器和3个不同取压间距的差压变送器对流动参数信号进行实时采集;
(2)去除噪声预处理:利用小波包对压力和压差波动信号进行去除噪声处理。
(a)对1个压力波动和3个不同取压间距的压差波动信号S1、S2、S3和S4按照下面的递归式(1)进行4层小波包分解:
g(k)=(-1)kh(1-k) (2)
(b)对分解的各个频带分别设定阈值,采用如下的浮动阈值形式
其中σ为噪声强度,用中位数MAD(υk)/0.6745来估计,υk为小波包分解系数,n为信号的长度。
(c)用非线性函数ξ=sgn(υk)(|υk|-tk)作用于信号的小波包分解系数υk,得到系数υk,即将浮动阈值以上的系数减去阈值后并保留下来。
(3)在对去除噪声后的1个压力和3个压差波动信号 进行经验模式分解,提取流型的特征向量。EMD方法从本质上讲是对一个信号进行平稳化处理,其结果是将信号中不同尺度的波动或趋势逐级分解出来,产生一系列具有不同特征尺度的数据序列,每一个序列称为一个固有模态函数IMF。IMF分量必须满足下面两个条件:一是其极点数和零点数相同(或最多相差一个),二是其上下包络线关于时间轴局部对称。
EMD方法步骤如下:假设任何信号都由不同的IMF组成,每个IMF可以是线性的,也可以是非线性的,这样任何一个信号就可以分解为有限个IMF之和,则IMF可以按以下方法“筛选”(Sift):
(a)确定信号x(t)的所有局部极值点,将所有极大值点用三次样条线连接起来形成上包络线,将所有极小值点用三次样条线连接起来形成下包络线,这两条包络线包络了所有的信号数据。
(b)将两条包络线的均值记为μ1,求出
y1(t)=x(t)-μ1 (4)
(c)判断y1(t)是否为IMF,若y1(t)不满足IMF条件,则将y1(t)作为原始数据,此时,记y1(t)=c1(t),则c1(t)为信号x1(t)的第一个IMF分量,它代表信号x1(t)中最高频率的分量;
(d)将c1(t)从x1(t)中分离出来,即得到一个去掉高频分量的差值信号r1(t),即有
r1(t)=x1(t)-c1(t) (5)
将r1(t)作为原始数据,重复步骤(1)、(2)和(3),得到第二个IMF分量c2(t),重复n次,得到n个IMF分量。这样就有
当cn(t)或rn(t)满足给定的终止条件(通常使rn(t)成为一个单调函数)时,循环结束,由(2)式和(3)式可得到:
式中,rn(t)为残余函数,代表信号的平均趋势。而各IMF分量c1(t),c2(t),Λ,cn(t)分别包含了信号从高到低不同频率段的成分,每一频率段所包含的频率成分都是不同的,且随信号本身的变化而变化。
(e)特征向量的提取。1个压力波动信号和3个压差波动信号的n个IMF分量分别记为c11(t),c12(t),Λ,、c21(t),c22(t),Λ,、c31(t),c32(t),Λ,、c41(t),c42(t),Λ,。将这些IMF分量组成一个特征矩阵X,即
式中:X为p×n的特征矩阵,p为采样点数,n=n1+n2+n3+n4为4个原始信号的所有IMF分量的个数,xi(i=1,2,Λ,n)为p维向量。
(f)对矩阵X用下式变换为相关矩阵,即主成分矩阵R,并计算出相关矩阵R的m个特征值。
①将原始变量进行标准化
各指标量纲往往不同,分析时,不同的量纲和数量级会引出新的问题,故要标准化,使标准化的变量平均值为0,方差为1,即
②计算X的相关矩阵R
标准化后,相关阵和协方差阵相同,且R的子分量可表示为
③求相关矩阵R的特征值及特征向量
计算出各特征值(λ1,λ2,Λ,λn)和各特征向量(u1,u2,Λ,un)。
④选择主成分
计算各主成分的贡献率及累计贡献率,主成分的个数视具体问题而定,本发明取累计贡献率为90%所对应的m个变量作为主成分。
⑤计算得分
T=[μ1,μ2,Λ,μm] (11)
则T为不同流型的特征向量。
(4)支持向量机分类器。采用支持向量机来实现流型的识别。水平管内油气水多相流流型有4种,即分层流、间歇流、泡状流和环状流。根据上述分类情况,采用“一对多”的策略,构造4个二值分类器,对每一个流型构造一个支持向量机,4个支持向量机的输出是一个4维向量,每一个分量表示样本是否对应该流型,训练好的支持向量机可实现对流型的识别。
本发明基于主成分分析和支持向量机的油气水多相流流型识别方法的优点在于:
1.对1个压力和3个差压变送器的信号进行融合,提取的特征能全面反映流型的信息,因此,识别准确性和可靠性高。
2.采用支持向量机来实现流型的识别,计算速度快,有利于流型的在线识别。
附图说明
图1是基于主成分分析和支持向量机的油气水多相流流型识别方法流程图。
图2是基于主成分分析和支持向量机的油气水多相流流型信号采集装置的结构示意图。
具体实施方式
下面利用附图和实施例对本发明基于主成分分析和支持向量机的油气水多相流流型识别方法作进一步说明。
如图2所示,基于主成分分析和支持向量机的油气水多相流流型信号采集装置,具有多相流水平管路16,在多相流水平管路16上套置的取压环1a、2b、3c、4d、5e、6f、7g,取压环1a与取压管8a连接,取压环2b、3c与取压管8b连接,取压环4d、5e与取压管8c连接,取压环6f、7g与取压管8d连接,取压管8a与压力变送器9a连接,取压管8b、8c、8d分别与差压变送器10b、11c、12d连接,压力变送器9a、差压变送器10b、11c、12d均与数据采集卡13电连接,数据采集卡13与计算机14电连接。所述取压环2b、3c的间距为5倍管径D,取压环4d、5e的间距为10倍管径D,取压环6f、7g的间距为15倍管径D。其中取压环1a距离入口300mm,取压环1a与取压环2b相隔200mm,取压环3c与取压环4d相隔200mm,取压环5e与取压环6f相隔200mm。压力变送器9a、差压变送器10b、11c、12d型号为PD-23的输出标准信号与型号为IMP3595的数据采集卡13相连,数据采集卡13与计算机14相连,由计算机14完成压差信号的数据处理,进而实现对油气水多相流流型信号的采集。计算机14还可连接打印机15。基于主成分分析和支持向量机的油气水多相流流型信号采集装置的软件程序是依据自动检测技术及计算机数据处理技术编制的,为本领域技术人员所熟悉的技术。
如图1所示,本发明基于主成分分析和支持向量机的油气水多相流流型识别方法包括下述步骤:
(1)对不同流型的流动参数信号进行采集:在基于主成分分析和支持向量机的油气水多相流流型信号采集装置中,利用取压间距分别为5倍、10倍和15倍管径D;采样频率为256Hz;采样点数为1024点,采样时间为16s;用1个压力变送器9a和3个不同取压间距的差压变送器10b、11c、12d对流动参数信号进行实时采集;
(2)应用小波理论对采集的压差波动信号进行去除噪声预处理:(a)对1个压力波动和3个不同取压间距的压差波动信号S1、S2、S3和S4按照下面的递归式(1)进行4层小波包分解:
g(k)=(-1)kh(1-k) (2)
(b)对分解的各个频带分别设定阈值,采用如下的浮动阈值形式
其中σ为噪声强度,用中位数MAD(υk)/0.6745来估计,υk为小波包分解系数,n为信号的长度。
(c)用非线性函数ξ=sgn(υk)(|υk|-tk)作用于信号的小波包分解系数υk,得到系数υk,即将浮动阈值以上的系数减去阈值后并保留下来。
(d)由系数υk重建原信号,从而达到去除噪声的目的。
(3)提取流型的特征:运用EMD对采集到的压差波动信号x(t)按如下步骤进行分解:
(a)确定信号x(t)的所有局部极值点,将所有极大值点用三次样条线连接起来形成上包络线,将所有极小值点用三次样条线连接起来形成下包络线,这两条包络线包络了所有的信号数据。
(b)将两条包络线的均值记为μ1,求出
y1(t)=x(t)-μ1 (4)
(c)判断y1(t)是否为IMF,若y1(t)不满足IMF条件,则将y1(t)作为原始数据,此时,记y1(t)=c1(t),则c1(t)为信号x1(t)的第一个IMF分量,它代表信号x1(t)中最高频率的分量;
(d)将c1(t)从x1(t)中分离出来,即得到一个去掉高频分量的差值信号r1(t),即有
r1(t)=x1(t)-c1(t) (5)
将r1(t)作为原始数据,重复步骤(1)、(2)和(3),得到第二个IMF分量c2(t),重复n次,得到n个IMF分量。这样就有
当cn(t)或rn(t)满足给定的终止条件(通常使rn(t)成为一个单调函数)时,循环结束,由(2)式和(3)式可得到:
式中,rn(t)为残余函数,代表信号的平均趋势。而各IMF分量c1(t),c2(t),Λ,cn(t)分别包含了信号从高到低不同频率段的成分,每一频率段所包含的频率成分都是不同的,且随信号本身的变化而变化。
(e)特征向量的提取。1个压力波动信号和3个压差波动信号的n个IMF分量分别记为c11(t),c12(t),Λ,、c21(t),c22(t),Λ,、c31(t),c32(t),Λ,、c41(t),c42(t),Λ,。将这些IMF分量组成一个特征矩阵X,即
式中:X为p×n的特征矩阵,p为采样点数,n=n1+n2+n3+n4为4个原始信号的所有IMF分量的个数。
(f)对矩阵X用下式变换为相关矩阵,即主成分矩阵R,并计算出相关矩阵R的m个特征值。
①将原始变量进行标准化
各指标量纲往往不同,分析时,不同的量纲和数量级会引出新的问题,故要标准化,使标准化的变量平均值为0,方差为1,即
②计算X的相关矩阵R
标准化后,相关阵和协方差阵相同,且R的子分量可表示为
③求相关矩阵R的特征值及特征向量
计算出各特征值(λ1,λ2,Λ,λn)和各特征向量(u1,u2,Λ,un)。
④选择主成分
计算各主成分的贡献率及累计贡献率,主成分的个数视具体问题而定,本发明取累计贡献率为90%所对应的m个变量作为主成分。
⑤计算得分
T=[μ1,μ2,Λ,μm] (11)
则T为不同流型的特征向量。
(4)采用支持向量机来实现流型的识别:水平管内油气水多相流流型有4种,即分层流、间歇流、泡状流和环状流。根据上述分类情况,采用“一对多”的策略,构造4个二值分类器,分别表示为SVM1、SVM2、SVM3、SVM4。如图2所示。对每一个流型构造一个支持向量机,每个支持向量机的输出yi∈{-1,+1},(i=1,2,3,4),则4个支持向量机的输出是一个4维向量,每一个分量表示样本是否对应该流型,(1,-1,-1,-1)代表分层流,(-1,1,-1,-1)代表间歇流,(-1,-1,1,-1)代表泡状流,(-1,-1,-1,1)代表环状流。用120个(4种流型各30个样本)对支持向量机进行训练,通过比较,确定本发明的核函数为径向基核函数,即 式中σ取0.05。训练好的支持向量机可实现对流型的识别。再用200个样本进行测试,结果如表1所示(部分测试样本)。
表1
Claims (5)
1.一种基于主成分分析和支持向量机的油气水多相流流型识别方法,其特征在于,将主成分分析方法与支持向量机识别技术相结合,自动对油气水多相流管流流型进行识别,其方法包括下述步骤:
1)对不同流型的流动参数信号进行采集;
2)对压力和压差波动信号进行小波包去除噪声预处理;
3)对去除噪声后的压力和压差波动信号进行经验模式分解,将分解得到的各个信号的固有模式函数组成一个特征矩阵,再利用主成分分析方法提取流型的特征向量;
4)将经验模式分解的能量特征向量作为支持向量机的输入样本,输入与输出的关系由支持向量机来完成,并构造4个二值分类器,由径向基核函数支持向量机完成对流型的识别。
2.根据权利要求1所述的基于主成分分析和支持向量机的油气水多相流流型识别方法,其特征是:所述的步骤1)是利用取压间距分别为5倍、10倍和15倍管径(D);采样频率为256Hz;采样点数为1024点;1个压力变送器和3个不同取压间距的差压变送器对流动参数信号进行实时采集。
3.根据权利要求1所述的基于主成分分析和支持向量机的油气水多相流流型识别方法,其特征是:所述的步骤2)的小波包去除噪声预处理,
(a)对1个压力波动和3个不同取压间距的压差波动信号S1、S2、S3和S4按照下面的递归式(1)进行4层小波包分解:
g(k)=(-1)kh(1-k) (2)
(b)对分解的各个频带分别设定阈值,采用如下的浮动阈值形式:
其中σ为噪声强度,用中位数MAD(υk)/0.6745来估计,υk为小波包分解系数,n为信号的长度;
(c)用非线性函数ξ=sgn(υk)(|υk|-tk)作用于信号的小波包分解系数υk,得到系数υk,即将浮动阈值以上的系数减去阈值后并保留下来;
(d)由系数υk重建原信号,从而达到去除噪声。
4.根据权利要求1所述的基于主成分分析和支持向量机的油气水多相流流型识别方法,其特征是:所述的步骤3)流型特征向量的提取,
(1)利用经验模式分解技术对1个压力波动信号和3个压差波动信号进行分解,得到n个IMF分量分别记为c11(t),c12(t),Λ,c21(t),c22(t),Λ,c31(t),c32(t),Λ,c41(t),c42(t),Λ,将这些IMF分量组成一个特征矩阵X,即
式中:X为p×n的特征矩阵,p为采样点数,n=n1+n2+n3+n4为4个原始信号的所有IMF分量的个数,xi(i=1,2,Λ,n)为p维向量;
(2)将原始变量进行标准化:各变量的不同量纲不能直接计算,故要标准化,使标准化的变量平均值为0,方差为1,即:
(3)计算X的相关矩阵R:标准化后,相关阵和协方差阵相同,且R的子分量可表示为:
(4)求相关矩阵R的特征值及特征向量:计算出各特征值(λ1,λ2,Λ,λn)和各特征向量(u1,u2,Λ,un);
(5)选择主成分:计算各主成分的贡献率及累计贡献率,主成分的个数视具体问题而定,本发明取累计贡献率为90%所对应的m个变量作为主成分;
T=[μ1,μ2,Λ,μm] (9)
则T为不同流型的特征向量。
5.根据权利要求1所述的基于主成分分析和支持向量机的油气水多相流流型识别方法,其特征是:所述的步骤4)是采用支持向量机来实现流型的识别,水平管内油气水多相流流型有4种,即分层流、间歇流、泡状流和环状流,采用“一对多”的策略,构造4个二值分类器,对每一个流型构造一个支持向量机,4个支持向量机的输出是一个4维向量,每一个分量表示样本是否对应该流型,训练好的支持向量机可实现对流型的识别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100170900A CN100507509C (zh) | 2006-08-01 | 2006-08-01 | 基于主成分分析和支持向量机的油气水多相流流型识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100170900A CN100507509C (zh) | 2006-08-01 | 2006-08-01 | 基于主成分分析和支持向量机的油气水多相流流型识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1904581A CN1904581A (zh) | 2007-01-31 |
CN100507509C true CN100507509C (zh) | 2009-07-01 |
Family
ID=37673886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100170900A Expired - Fee Related CN100507509C (zh) | 2006-08-01 | 2006-08-01 | 基于主成分分析和支持向量机的油气水多相流流型识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100507509C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269684A (zh) * | 2011-06-22 | 2011-12-07 | 浙江大学 | 小管道液液两相流流型辨识系统及方法 |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102841979A (zh) * | 2012-06-26 | 2012-12-26 | 广州市新之地环保产业有限公司 | 一种基于pca-ga-svr的污泥回流量预测方法 |
CN103257181B (zh) * | 2013-05-24 | 2015-08-26 | 中国矿业大学 | 一种高速运动状态下金属板材缺陷检测装置与方法 |
CN103674478B (zh) * | 2013-12-11 | 2016-01-06 | 中国石油大学(华东) | 低含气率的气液两相流流型识别方法 |
CN103967478B (zh) * | 2014-05-21 | 2017-10-27 | 北京航空航天大学 | 一种基于电导探针的垂直井流型辨识方法 |
CN104698837B (zh) * | 2014-12-11 | 2017-10-20 | 华侨大学 | 一种时变线性结构工作模态参数识别方法、装置及应用 |
CN106972966A (zh) * | 2017-03-28 | 2017-07-21 | 国网辽宁省电力有限公司阜新供电公司 | 一种通信网络中网络流量异常的实时检测方法 |
CN107991706B (zh) * | 2017-11-24 | 2020-07-31 | 华北科技学院 | 基于小波包多重阈值和改进经验模态分解的煤层水力压裂微震信号联合降噪方法 |
CN108388848B (zh) * | 2018-02-07 | 2022-02-22 | 西安石油大学 | 一种多尺度油气水多相流动力学特性分析方法 |
CN109538944B (zh) * | 2018-12-03 | 2020-07-07 | 北京无线电计量测试研究所 | 一种管道泄漏检测方法 |
CN110044533B (zh) * | 2019-04-11 | 2020-10-20 | 大连理工大学 | 一种基于emd能量熵与支持向量机相结合的螺栓预紧力监测方法 |
CN111503010B (zh) * | 2019-09-18 | 2022-03-18 | 山东建筑大学 | 一种离心泵宽带辐射噪声信号特征提取方法 |
-
2006
- 2006-08-01 CN CNB2006100170900A patent/CN100507509C/zh not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
基于支持向量机和小波包能量特征的气液两相流流型识别方法. 孙斌,周云龙.中国电机工程学报,第25卷第17期. 2005 |
基于支持向量机和小波包能量特征的气液两相流流型识别方法. 孙斌,周云龙.中国电机工程学报,第25卷第17期. 2005 * |
油气水多相流流型智能识别系统的设计与实现. 吴浩江,周芳德.西安交通大学学报,第34卷第3期. 2000 |
油气水多相流流型智能识别系统的设计与实现. 吴浩江,周芳德.西安交通大学学报,第34卷第3期. 2000 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269684A (zh) * | 2011-06-22 | 2011-12-07 | 浙江大学 | 小管道液液两相流流型辨识系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1904581A (zh) | 2007-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100507509C (zh) | 基于主成分分析和支持向量机的油气水多相流流型识别方法 | |
CN111882446B (zh) | 一种基于图卷积网络的异常账户检测方法 | |
CN101984340B (zh) | 一种基于人工智能的软性磨粒两相湍流流型识别方法 | |
CN108595414B (zh) | 基于源汇空间变量推理的土壤重金属企业污染源识别方法 | |
CN109102005A (zh) | 基于浅层模型知识迁移的小样本深度学习方法 | |
CN104807589B (zh) | 一种集输-立管系统内气液两相流流型的在线识别方法 | |
CN100485359C (zh) | 基于信息融合的气液两相流流型识别方法及流型信号采集装置 | |
CN105389480A (zh) | 多类不平衡基因组学数据迭代集成特征选择方法及系统 | |
CN103412557A (zh) | 一种适于非线性过程在线监控的工业故障检测与诊断方法 | |
CN111965121B (zh) | 一种基于光谱分析技术的煤炭焦渣的快速分类检测方法 | |
CN104688252A (zh) | 采用方向盘转角信息的驾驶人疲劳状态检测方法 | |
CN103914705A (zh) | 一种基于多目标免疫克隆的高光谱图像分类和波段选择的方法 | |
CN116049637A (zh) | 基于流动噪声信号的气液两相流流型识别方法 | |
CN115795380A (zh) | 一种基于孤立森林和加权随机森林的烟气制酸数据清洗及优化方法 | |
CN106650801A (zh) | 一种基于gps数据的多类型车辆分类方法 | |
CN105823634A (zh) | 基于时频关联向量卷积玻尔兹曼机的轴承损伤识别方法 | |
CN106326335A (zh) | 一种基于显著属性选择的大数据归类方法 | |
CN106055645A (zh) | 一种用于高维数据分析的维度重要性评估方法 | |
CN113283467A (zh) | 一种基于平均损失和逐类选择的弱监督图片分类方法 | |
CN117475191A (zh) | 一种特征对齐域对抗神经网络的轴承故障诊断方法 | |
CN115374858B (zh) | 基于混合集成模型的流程工业生产品质的智能诊断方法 | |
CN107122394A (zh) | 异常数据检测方法和装置 | |
CN116933151A (zh) | 一种基于闪锌矿微量元素判别矿床类型方法 | |
CN116610990A (zh) | 一种基于特征空间差异化的断路器隐患识别方法及装置 | |
CN117235489A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090701 Termination date: 20120801 |