CN1904581A - Oil gas water multiphase flow type identification method based on main component analysis and supporting vector machine - Google Patents
Oil gas water multiphase flow type identification method based on main component analysis and supporting vector machine Download PDFInfo
- Publication number
- CN1904581A CN1904581A CN 200610017090 CN200610017090A CN1904581A CN 1904581 A CN1904581 A CN 1904581A CN 200610017090 CN200610017090 CN 200610017090 CN 200610017090 A CN200610017090 A CN 200610017090A CN 1904581 A CN1904581 A CN 1904581A
- Authority
- CN
- China
- Prior art keywords
- vector machine
- pressure
- signal
- flow
- flow type
- 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.)
- Granted
Links
Images
Landscapes
- Complex Calculations (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention relates to an oil/gas/water multi-phase flow type identification method based on main constituent analysis and supporting vector machine. The feature is that: taking collection to flow parameter signal of different flow type by using horizontal tubes, pressure sampling ring, pressure sampling tube, pressure difference transmitter, data collection card and computer; taking denoising processing to pressure difference fluctuation signal between the collected different pressure sampling distance; taking empiric model decomposition to the signal to gain the fixed modal function, and forming a matrix X from the modal functions, using the main constituent analysis to gain the eigenvector of the flow type; using the supporting vector machine to accomplish the mapping from feature space to flow type space to realize the final flow type identification. The eigenvector fusing the information from plural sensors and would fully reflect the information of flow type. It has the feature of rapid learning speed and strong classifying ability.
Description
Technical field
The present invention relates to the multiphase flow measurement technical field, relate in particular to a kind of oil gas water multiphase flow type identification method based on principal component analysis (PCA) and support vector machine.
Background technology
Traditional flow type identification method has two big classes: a class is to adopt experimental technique to make flow regime map; Another kind of is the transformation criterion relational expression that the analysis according to the convection transformation mechanism obtains, and utilizes on-the-spot flow parameter to determine flow pattern.Because flow parameter needs the problem that solves often in production reality, therefore traditional method can't be used widely.
The modern surveying method of flow pattern can be divided into direct measurement and indirect measurement method on principle of work.Direct measuring method commonly used has ocular estimate, high-speed photography method etc., and this method has certain subjectivity, may draw different conclusions for a plurality of observers.Indirect measurement method is the fluctuation signals such as pressure, pressure reduction and void fraction that utilize the instrument and equipment measurement of measuring, and by signal being analyzed, extracted feature, discerns flow pattern in conjunction with neural network isotype recognition technology.This method has two key issues: the one, and the extraction of flow pattern proper vector; The 2nd, the choosing of network model, wherein feature extraction is particularly crucial.Existing method mainly is to utilize single pressure differential pressure fluctuation signal to extract the single feature of signal, and the feature that causes extracting can't reflect the information of flow pattern comprehensively, causes the flow pattern discrimination not high.
Summary of the invention
The objective of the invention is to overcome the defective of above-mentioned prior art, propose a kind of flow pattern identification accuracy, reliability height, the oil gas water multiphase flow type identification method based on principal component analysis (PCA) and support vector machine of highly versatile.
The objective of the invention is to be realized by following technical scheme: a kind of oil gas water multiphase flow type identification method based on principal component analysis (PCA) and support vector machine is characterized in that it may further comprise the steps:
(1) the flow parameter signal of different flow patterns is gathered: utilize the pressure spacing to be respectively 5 times, 10 times and 15 times of caliber D; Sample frequency is 256Hz; Sampling number is 1024 points; The differential pressure transmitter of 1 pressure unit and 3 different pressure spacings is gathered in real time to the parameter signal that flows;
(2) remove the noise pre-service: utilize wavelet packet that pressure and pressure-difference fluctuation signal are removed noise processed.
(a) to the pressure-difference fluctuation signal S of 1 pressure surge and 3 different pressure spacings
1, S
2, S
3And S
4Carry out 4 layers of WAVELET PACKET DECOMPOSITION according to following recursion (1):
H in the formula (k) is the Hi-pass filter group; G (k) is the low-pass filter group; u
0(t)=and (t), be scaling function; u
1(t)=and ψ (t), be wavelet function; The yardstick of k for decomposing.H (k) and g (k) satisfy orthogonality relation:
g(k)=(-1)
kh(1-k) (2)
(b) each frequency band difference setting threshold to decomposing adopts following floating threshold form
Wherein σ is a noise intensity, with median MAD (υ
k)/0.6745 (υ
kBe the WAVELET PACKET DECOMPOSITION coefficient) estimate that n is the length of signal.
(c) with nonlinear function ζ=sgn (υ
k) (| υ
k|-t
k) act on the WAVELET PACKET DECOMPOSITION coefficient υ of signal
k, obtain coefficient υ
k, soon the above coefficient of floating threshold deducts after the threshold value and remains.
(d) by coefficient υ
kRebuild original signal, thereby reach the removal noise, the signal that obtains removing behind the noise is S
1', S
2', S
3' and S
4'.
(3) to 1 pressure and 3 pressure-difference fluctuation signal S after removing noise
1', S
2', S
3', S
4' carry out empirical mode decomposition, extract the proper vector of flow pattern.The EMD method is a signal to be carried out tranquilization handle in essence, consequently the fluctuation or the trend of different scale in the signal are decomposed out step by step, produce a series of data sequences with different characteristic yardstick, each sequence is called an intrinsic mode function IMF.The IMF component must satisfy following two conditions: the one, and its limit number and number at zero point identical (or differing at most), the 2nd, envelope is local symmetrical about time shaft up and down for it.
The EMD method step is as follows: suppose that any signal all is made up of different IMF, each IMF can be linear, also can be non-linear, any one signal just can be decomposed into limited IMF sum like this, and then " screening " (Sift) by the following method for IMF:
(a) determine all Local Extremum of signal x (t), all maximum points are coupled together with the cubic spline line form the coenvelope line, all minimum points are coupled together with the cubic spline line form the lower envelope line, these two envelope envelopes all signal datas.
(b) average with two envelopes is designated as μ
1, obtain
y
1(t)=x(t)-μ
1 (4)
(c) judge y
1(t) whether be IMF, if y
1(t) do not satisfy the IMF condition, then with y
1(t) as raw data, at this moment, note y
1(t)=c
1(t), c then
1(t) be signal x
1(t) first IMF component, its representation signal x
1(t) component of highest frequency in;
(d) with c
1(t) from x
1(t) separate in, promptly obtain a difference signal r who removes high fdrequency component
1(t), promptly have
r
1(t)=x
1(t)-c
1(t) (5)
With r
1(t) as raw data, repeating step (1), (2) and (3) obtain second IMF component c
2(t), repeat n time, obtain n IMF component.So just have
Work as c
n(t) or r
n(t) satisfy given end condition and (make r usually
n(t) become a monotonic quantity) time, loop ends can be obtained by (2) formula and (3) formula:
In the formula, r
n(t) be remaining function, the average tendency of representation signal.And each IMF component c
1(t), c
2(t), Λ, c
n(t) comprised the signal composition of different frequency section from high to low respectively, the frequency content that each frequency band comprised all is different, and changes with the variation of signal itself.
(e) extraction of proper vector.N IMF component of 1 Pressure Fluctuation Signal and 3 pressure-difference fluctuation signals is designated as c respectively
11(t), c
12(t), Λ, c
1n1(t), c
21(t), c
22(t), Λ, c
2n2(t), c
31(t), c
32(t), Λ, c
3n3(t), c
41(t), c
42(t), Λ, c
4n4(t).These IMF components are formed this matrix of spy X, promptly
In the formula: X is the eigenmatrix of p * n, and p is a sampling number, n=n
1+ n
2+ n
3+ n
4Be the number of all IMF components of 4 original signals, x
i(i=1,2, Λ n) is the p dimensional vector.
(f) matrix X is transformed to correlation matrix with following formula, i.e. major component matrix R, and calculate m the eigenwert of correlation matrix R.
1. original variable is carried out standardization
Each index dimension is often different, and during analysis, the different dimensions and the order of magnitude can be drawn new problem, so want standardization, making standardized average of variable is 0, and variance is 1, promptly
i=1,2,Λ,n;j=1,2,Λ,p (9)
In the formula: x
Ij *Be the observed quantity after the standardization, x
iBe mean value,
Be standard deviation.
2. calculate the correlation matrix R of X
After the standardization, relevant battle array is identical with covariance matrix, and the subcomponent of R can be expressed as
In the formula: x
Ti *, x
Ij *Be respectively capable i row of t and j column element in the eigenmatrix.
3. ask eigenwert and the proper vector of correlation matrix R
Calculate each eigenwert (λ
1>λ
2>Λ>λ
n) and each proper vector (u
1, u
2, Λ, u
n).
4. select major component
Calculate the contribution rate and the contribution rate of accumulative total of each major component, the number of major component is decided on particular problem, and it is that 90% a pairing m variable is as major component that the present invention gets contribution rate of accumulative total.
5. count the score
Calculate each principal component scores according to the major component expression formula.In principal component analysis (PCA),
Be the contribution rate of major component, ratio
The percentage of the primary data information (pdi) of k major component representative before the reflection.Order
T=[μ
1,μ
2,Λ,μ
m] (11)
Then T is the proper vector of different flow patterns.
(4) support vector machine classifier.Adopt support vector machine to realize the identification of flow pattern.Oil gas water multiphase flow type has 4 kinds, i.e. stratified flow, intermittent flow, bubble flow and annular flow in the horizontal tube.According to above-mentioned classification situation, adopt the strategy of " one-to-many ", construct 4 two-value sorters, to support vector machine of each flow pattern structure, the output of 4 support vector machine is one 4 dimensional vectors, whether each representation in components sample is to should flow pattern, and the support vector machine that trains can realize the identification of convection.
The advantage that the present invention is based on the oil gas water multiphase flow type identification method of principal component analysis (PCA) and support vector machine is:
1. the signal to 1 pressure and 3 differential pressure transmitters merges, and the feature of extraction can reflect the information of flow pattern comprehensively, therefore, and identification accuracy and reliability height.
2. adopt support vector machine to realize the identification of flow pattern, computing velocity is fast, helps the ONLINE RECOGNITION of flow pattern.
Description of drawings
Fig. 1 is based on the oil gas water multiphase flow type identification method process flow diagram of principal component analysis (PCA) and support vector machine.
Fig. 2 is based on the structural representation of the oil gas water multiphase flow type signal pickup assembly of principal component analysis (PCA) and support vector machine.
Embodiment
Utilize drawings and Examples that the oil gas water multiphase flow type identification method that the present invention is based on principal component analysis (PCA) and support vector machine is described further below.
As shown in Figure 2, oil gas water multiphase flow type signal pickup assembly based on principal component analysis (PCA) and support vector machine, has the horizontal pipeline 16 of polyphasic flow, the pressure ring 1a that on the horizontal pipeline 16 of polyphasic flow, is nested with, 2b, 3c, 4d, 5e, 6f, 7g, pressure ring 1a is connected with pressure pipe 8a, pressure ring 2b, 3c is connected with pressure pipe 8b, pressure ring 4d, 5e is connected with pressure pipe 8c, pressure ring 6f, 7g is connected with pressure pipe 8d, pressure pipe 8a is connected with pressure unit 9a, pressure pipe 8b, 8c, 8d respectively with differential pressure transmitter 10b, 11c, 12d connects, pressure unit 9a, differential pressure transmitter 10b, 11c, 12d all is electrically connected with data collecting card 13, and data collecting card 13 is electrically connected with computing machine 14.The spacing of described pressure ring 2b, 3c is 5 times of caliber D, and the spacing of pressure ring 4d, 5e is 10 times of caliber D, and the spacing of pressure ring 6f, 7g is 15 times of caliber D.Wherein pressure ring 1a is apart from inlet 300mm, pressure ring 1a and the pressure ring 2b 200mm of being separated by, pressure ring 3c and the pressure ring 4d 200mm of being separated by, pressure ring 5e and the pressure ring 6f 200mm of being separated by.Pressure unit 9a, differential pressure transmitter 10b, 11c, 12d model are that outputting standard signal and the model of PD-23 is that the data collecting card 13 of IMP3595 links to each other, data collecting card 13 links to each other with computing machine 14, finish the data processing of pressure difference signal by computing machine 14, and then realize collection the oil gas water multiphase flow type signal.Computing machine 14 also can connect printer 15.According to Automatic Measurement Technique and microcomputer data processing establishment, be the technology that those skilled in the art were familiar with based on the software program of the oil gas water multiphase flow type signal pickup assembly of principal component analysis (PCA) and support vector machine.
As shown in Figure 1, the oil gas water multiphase flow type identification method that the present invention is based on principal component analysis (PCA) and support vector machine comprises the steps:
(1) the flow parameter signal of different flow patterns is gathered: in oil gas water multiphase flow type signal pickup assembly, utilize the pressure spacing to be respectively 5 times, 10 times and 15 times of caliber D based on principal component analysis (PCA) and support vector machine; Sample frequency is 256Hz; Sampling number is 1024 points, and the sampling time is 16s; Differential pressure transmitter 10b, 11c, 12d with 1 pressure unit 9a and 3 different pressure spacings gather in real time to the parameter signal that flows;
(2) use wavelet theory the pressure-difference fluctuation signal of gathering is removed the noise pre-service: (a) to the pressure-difference fluctuation signal S of 1 pressure surge and 3 different pressure spacings
1, S
2, S
3And S
4Carry out 4 layers of WAVELET PACKET DECOMPOSITION according to following recursion (1):
H in the formula (k) is the Hi-pass filter group; G (k) is the low-pass filter group; u
0(t)=and (t), be scaling function; u
1(t)=and ψ (t), be wavelet function; The yardstick of k for decomposing.H (k) and g (k) satisfy orthogonality relation:
g(k)=(-1)
kh(1-k) (2)
(b) each frequency band difference setting threshold to decomposing adopts following floating threshold form
Wherein σ is a noise intensity, with median MAD (υ
k)/0.6745 (υ
kBe the WAVELET PACKET DECOMPOSITION coefficient) estimate that n is the length of signal.
(c) with nonlinear function ζ=sgn (υ
k) (| υ
k|-t
k) act on the WAVELET PACKET DECOMPOSITION coefficient υ of signal
k, obtain coefficient υ
k, soon the above coefficient of floating threshold deducts after the threshold value and remains.
(d) by coefficient υ
kRebuild original signal, thereby reach the purpose of removing noise.
(3) feature of extraction flow pattern: utilization EMD decomposes as follows to the pressure-difference fluctuation signal x (t) that collects:
(a) determine all Local Extremum of signal x (t), all maximum points are coupled together with the cubic spline line form the coenvelope line, all minimum points are coupled together with the cubic spline line form the lower envelope line, these two envelope envelopes all signal datas.
(b) average with two envelopes is designated as μ
1, obtain
y
1(t)=x(t)-μ
1 (4)
(c) judge y
1(t) whether be IMF, if y
1(t) do not satisfy the IMF condition, then with y
1(t) as raw data, at this moment, note y
1(t)=c
1(t), c then
1(t) be signal x
1(t) first IMF component, its representation signal x
1(t) component of highest frequency in;
(d) with c
1(t) from x
1(t) separate in, promptly obtain a difference signal r who removes high fdrequency component
1(t), promptly have
r
1(t)=x
1(t)-c
1(t) (5)
With r
1(t) as raw data, repeating step (1), (2) and (3) obtain second IMF component c
2(t), repeat n time, obtain n IMF component.So just have
Work as c
n(t) or r
n(t) satisfy given end condition and (make r usually
n(t) become a monotonic quantity) time, loop ends can be obtained by (2) formula and (3) formula:
In the formula, r
n(t) be remaining function, the average tendency of representation signal.And each IMF component c
1(t), c
2(t), Λ, c
n(t) comprised the signal composition of different frequency section from high to low respectively, the frequency content that each frequency band comprised all is different, and changes with the variation of signal itself.
(e) extraction of proper vector.N IMF component of 1 Pressure Fluctuation Signal and 3 pressure-difference fluctuation signals is designated as c respectively
11(t), c
12(t), Λ, c
1n1(t), c
21(t), c
22(t), Λ, c
2n2(t), c
31(t), c
32(t), Λ, c
3n3(t), c
41(t), c
42(t), Λ, c
4n4(t).These IMF components are formed this matrix of spy X, promptly
In the formula: X is the eigenmatrix of p * n, and p is a sampling number, n=n
1+ n
2+ n
3+ n
4It is the number of all IMF components of 4 original signals.
(f) matrix X is transformed to correlation matrix with following formula, i.e. major component matrix R, and calculate m the eigenwert of correlation matrix R.
1. original variable is carried out standardization
Each index dimension is often different, and during analysis, the different dimensions and the order of magnitude can be drawn new problem, so want standardization, making standardized average of variable is 0, and variance is 1, promptly
In the formula: x
Ij *Be the observed quantity after the standardization, x
iBe mean value,
Be standard deviation.
2. calculate the correlation matrix R of X
After the standardization, relevant battle array is identical with covariance matrix, and the subcomponent of R can be expressed as
In the formula: x
Ti *, x
Ij *Be respectively capable i row of t and j column element in the eigenmatrix.
3. ask eigenwert and the proper vector of correlation matrix R
Calculate each eigenwert (λ
1>λ
2>Λ>λ
n) and each proper vector (u
1, u
2, Λ, u
n).
4. select major component
Calculate the contribution rate and the contribution rate of accumulative total of each major component, the number of major component is decided on particular problem, and it is that 90% a pairing m variable is as major component that the present invention gets contribution rate of accumulative total.
5. count the score
In principal component analysis (PCA),
Be the contribution rate of major component, ratio
The percentage of the primary data information (pdi) of k major component representative before the reflection.Order
T=[μ
1,μ
2,Λ,μ
m] (11)
Then T is the proper vector of different flow patterns.
(4) adopt support vector machine to realize the identification of flow pattern: oil gas water multiphase flow type has 4 kinds, i.e. stratified flow, intermittent flow, bubble flow and annular flow in the horizontal tube.According to above-mentioned classification situation, adopt the strategy of " one-to-many ", construct 4 two-value sorters, be expressed as SVM1, SVM2, SVM3, SVM4 respectively.As shown in Figure 2.To support vector machine of each flow pattern structure, the output y of each support vector machine
i∈ 1 ,+1}, (i=1,2,3,4), then the output of 4 support vector machine is one 4 dimensional vectors, and whether each representation in components sample is to should flow pattern, (1 ,-1 ,-1 ,-1) represents stratified flow, (1,1 ,-1 ,-1) represents intermittent flow, (1 ,-1,1 ,-1) represents bubble flow, (1 ,-1 ,-1,1) represents annular flow.With 120 (4 kinds of each 30 samples of flow pattern) support vector machine is trained,, determine that kernel function of the present invention is radially basic kernel function, promptly by relatively
σ gets 0.05 in the formula.The support vector machine that trains can realize the identification of convection.Test result's (partial test sample) as shown in table 1 again with 200 samples.
Table 1
Test sample book | Output | Flow pattern | ||||||
0.47561 | 0.29038 | 0.11616 | 0.08258 | 0.02529 | 0.00434 | 0.00473 | (-1,-1,-1,1) | Annular flow |
0.48657 | 0.28142 | 0.10764 | 0.07964 | 0.02637 | 0.00538 | 0.00456 | (-1,-1,-1,1) | Annular flow |
0.48841 | 0.27038 | 0.11524 | 0.08357 | 0.02735 | 0.00469 | 0.00504 | (-1,-1,-1,1) | Annular flow |
0.57743 | 0.31873 | 0.06160 | 0.01991 | 0.00848 | 0.00563 | 0.00469 | (-1,-1,1,-1) | Bubble flow |
0.57628 | 0.31772 | 0.06036 | 0.01814 | 0.01584 | 0.00433 | 0.00297 | (-1,-1,1,-1) | Bubble flow |
0.57628 | 0.31772 | 0.06036 | 0.01814 | 0.01584 | 0.00433 | 0.00297 | (-1,-1,1,-1) | Bubble flow |
0.69828 | 0.11918 | 0.09776 | 0.06454 | 0.01191 | 0.004421 | 0.00229 | (-1,1,-1,-1) | Intermittent flow |
0.69828 | 0.11918 | 0.09776 | 0.06454 | 0.01191 | 0.004421 | 0.00229 | (-1,1,-1,-1) | Intermittent flow |
0.68550 | 0.23403 | 0.05095 | 0.02419 | 0.00263 | 0.001564 | 0.00115 | (-1,1,-1,-1) | Intermittent flow |
0.47509 | 0.28644 | 0.13284 | 0.07738 | 0.01721 | 0.005168 | 0.00364 | (1,-1,-1,-1) | Stratified flow |
0.42308 | 0.21099 | 0.16566 | 0.10569 | 0.074812 | 0.01078 | 0.00551 | (1,-1,-1,-1) | Stratified flow |
0.41356 | 0.23654 | 0.15324 | 0.12341 | 0.06244 | 0.01025 | 0.00631 | (1,-1,-1,-1) | Stratified flow |
Claims (5)
1. oil gas water multiphase flow type identification method based on principal component analysis (PCA) and support vector machine, it is characterized in that, principal component analytical method is combined with the support vector machine recognition technology, automatically oil gas water multiphase pipe stream flow pattern is discerned, its method comprises the steps:
1) the flow parameter signal of different flow patterns is gathered;
2) pressure and pressure-difference fluctuation signal are carried out the pre-service of wavelet packet removal noise;
3) pressure and the pressure-difference fluctuation signal of removing behind the noise carried out empirical mode decomposition, the natural mode function of each signal that decomposition is obtained is formed an eigenmatrix, utilizes principal component analytical method to extract the proper vector of flow pattern again;
4) with the energy feature vector of the empirical mode decomposition input sample as support vector machine, input is finished by support vector machine with the relation of output, and constructs 4 two-value sorters, is finished the identification of convection by radially basic kernel function support vector machine.
2. the oil gas water multiphase flow type identification method based on principal component analysis (PCA) and support vector machine according to claim 1 is characterized in that: described step 1 is to utilize the pressure spacing to be respectively 5 times, 10 times and 15 times of calibers (D); Sample frequency is 256Hz; Sampling number is 1024 points; The differential pressure transmitter of 1 pressure unit and 3 different pressure spacings is gathered in real time to the parameter signal that flows.
3. the oil gas water multiphase flow type identification method based on principal component analysis (PCA) and support vector machine according to claim 1 is characterized in that: the wavelet packet of described step 2 is removed the noise pre-service,
(a) to the pressure-difference fluctuation signal S of 1 pressure surge and 3 different pressure spacings
1, S
2, S
3And S
4Carry out 4 layers of WAVELET PACKET DECOMPOSITION according to following recursion (1):
H in the formula (k) is the Hi-pass filter group, and g (k) is the low-pass filter group, u
0(t)= (t) is a scaling function, u
1(t)=and ψ (t) is a wavelet function, the yardstick of k for decomposing, and h (k) and g (k) satisfy orthogonality relation:
g(k)=(-1)
kh(1-k) (2)
(b) each frequency band difference setting threshold to decomposing, adopt following floating threshold form:
Wherein σ is a noise intensity, with median MAD (υ
k)/0.6745 (υ
kBe the WAVELET PACKET DECOMPOSITION coefficient) estimate that n is the length of signal;
(c) with nonlinear function ξ=sgn (υ
k) (| υ
k|-t
k) act on the WAVELET PACKET DECOMPOSITION coefficient υ of signal
k, obtain coefficient υ
k, soon the above coefficient of floating threshold deducts after the threshold value and remains;
(d) by coefficient υ
kRebuild original signal, thereby reach the removal noise.
4. the oil gas water multiphase flow type identification method based on principal component analysis (PCA) and support vector machine according to claim 1 is characterized in that: the extraction of described step 3 flow pattern proper vector,
(1) utilizes the empirical mode decomposition technology that 1 Pressure Fluctuation Signal and 3 pressure-difference fluctuation signals are decomposed, obtain n IMF component and be designated as c respectively
11(t), c
12(t), Λ, c
1n1(t), c
21(t), c
22(t), A, c
2n2(t), c
31(t), c
32(t), Λ, c
3n3(t), c
41(t), c
42(t), Λ, c
4n4(t); These IMF components are formed this matrix of spy X, promptly
In the formula: X is the eigenmatrix of p * n, and p is a sampling number, n=n
1+ n
2+ n
3+ n
4Be the number of all IMF components of 4 original signals, x
i(i=1,2, Λ n) is the p dimensional vector;
(2) original variable is carried out standardization: the different dimensions of each variable can not directly be calculated, so want standardization, making standardized average of variable is 0, and variance is 1, that is:
In the formula: x
Ij *Be the observed quantity after the standardization, x
iBe mean value,
Be standard deviation;
(3) the correlation matrix R of calculating X: after the standardization, relevant battle array is identical with covariance matrix, and the subcomponent of R can be expressed as:
In the formula: x
Ti *, x
Ij *Be respectively capable i row of t and j column element in the eigenmatrix;
(4) ask eigenwert and the proper vector of correlation matrix R: calculate each eigenwert (λ
1>λ
2>Λ>λ
n) and each proper vector (u
1, u
2, Λ, u
n);
(5) select major component: calculate the contribution rate and the contribution rate of accumulative total of each major component, the number of major component is decided on particular problem, and it is that 90% a pairing m variable is as major component that the present invention gets contribution rate of accumulative total;
(6) count the score: calculate each principal component scores according to the major component expression formula, in principal component analysis (PCA),
Be the contribution rate of major component, ratio
The percentage of the primary data information (pdi) of k major component representative before the reflection, order:
T=[μ
1,μ
2,Λ,μ
m] (9)
Then T is the proper vector of different flow patterns.
5. the oil gas water multiphase flow type identification method based on principal component analysis (PCA) and support vector machine according to claim 1, it is characterized in that: described step 4 is to adopt support vector machine to realize the identification of flow pattern, oil gas water multiphase flow type has 4 kinds in the horizontal tube, it is stratified flow, intermittent flow, bubble flow and annular flow, adopt the strategy of " one-to-many ", construct 4 two-value sorters, to support vector machine of each flow pattern structure, the output of 4 support vector machine is one 4 dimensional vectors, whether each representation in components sample is to should flow pattern, and the support vector machine that trains can realize the identification of convection.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100170900A CN100507509C (en) | 2006-08-01 | 2006-08-01 | Oil gas water multiphase flow type identification method based on main component analysis and supporting vector machine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100170900A CN100507509C (en) | 2006-08-01 | 2006-08-01 | Oil gas water multiphase flow type identification method based on main component analysis and supporting vector machine |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1904581A true CN1904581A (en) | 2007-01-31 |
CN100507509C CN100507509C (en) | 2009-07-01 |
Family
ID=37673886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100170900A Expired - Fee Related CN100507509C (en) | 2006-08-01 | 2006-08-01 | Oil gas water multiphase flow type identification method based on main component analysis and supporting vector machine |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100507509C (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102841979A (en) * | 2012-06-26 | 2012-12-26 | 广州市新之地环保产业有限公司 | Sludge reflux prediction method based on PCA-GA-SVR (principal component analysis-genetic algorithm-support vector regression) |
CN103257181A (en) * | 2013-05-24 | 2013-08-21 | 中国矿业大学 | Metal sheet defect detection device and method in high-speed motion state |
CN103674478A (en) * | 2013-12-11 | 2014-03-26 | 中国石油大学(华东) | Low-gas-content gas-liquid two phase flow pattern identifying method |
CN103967478A (en) * | 2014-05-21 | 2014-08-06 | 北京航空航天大学 | Method for identifying vertical well flow patterns based on conducting probe |
CN104698837A (en) * | 2014-12-11 | 2015-06-10 | 华侨大学 | Method and device for identifying operating modal parameters of linear time-varying structure and application of the device |
CN106972966A (en) * | 2017-03-28 | 2017-07-21 | 国网辽宁省电力有限公司阜新供电公司 | The real-time detection method of exception of network traffic in a kind of communication network |
CN107991706A (en) * | 2017-11-24 | 2018-05-04 | 华北科技学院 | Combine noise-reduction method based on wavelet packet multiple thresholds and the coal seam hydraulic fracture microseismic signals for improving empirical mode decomposition |
CN108388848A (en) * | 2018-02-07 | 2018-08-10 | 西安石油大学 | A kind of multiple dimensioned oil-gas-water three-phase flow mechanical characteristic analysis method |
CN109282965A (en) * | 2018-11-06 | 2019-01-29 | 中国海洋石油集团有限公司 | It is a kind of to collect the quick identification device of nocuousness flow pattern and method in defeated standpipe |
CN109538944A (en) * | 2018-12-03 | 2019-03-29 | 北京无线电计量测试研究所 | A kind of pipeline leakage detection method |
CN110044533A (en) * | 2019-04-11 | 2019-07-23 | 大连理工大学 | A kind of tools for bolts ' pretension force monitoring methods combined based on EMD Energy-Entropy with support vector machines |
CN111503010A (en) * | 2019-09-18 | 2020-08-07 | 山东建筑大学 | Method for extracting broadband radiation noise signal characteristics of centrifugal pump |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269684A (en) * | 2011-06-22 | 2011-12-07 | 浙江大学 | Small-diameter pipeline liquid-liquid two-phase flow flow pattern identification system and method |
-
2006
- 2006-08-01 CN CNB2006100170900A patent/CN100507509C/en not_active Expired - Fee Related
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102841979A (en) * | 2012-06-26 | 2012-12-26 | 广州市新之地环保产业有限公司 | Sludge reflux prediction method based on PCA-GA-SVR (principal component analysis-genetic algorithm-support vector regression) |
CN103257181A (en) * | 2013-05-24 | 2013-08-21 | 中国矿业大学 | Metal sheet defect detection device and method in high-speed motion state |
CN103257181B (en) * | 2013-05-24 | 2015-08-26 | 中国矿业大学 | Sheet metal defect detecting device and method under a kind of high-speed motion state |
CN103674478A (en) * | 2013-12-11 | 2014-03-26 | 中国石油大学(华东) | Low-gas-content gas-liquid two phase flow pattern identifying method |
CN103674478B (en) * | 2013-12-11 | 2016-01-06 | 中国石油大学(华东) | The Method for Discriminating Gas-liquid Two Phase Flow of low air void |
CN103967478B (en) * | 2014-05-21 | 2017-10-27 | 北京航空航天大学 | A kind of peupendicular hole meteor trail echoes method based on conducting probe |
CN103967478A (en) * | 2014-05-21 | 2014-08-06 | 北京航空航天大学 | Method for identifying vertical well flow patterns based on conducting probe |
CN104698837A (en) * | 2014-12-11 | 2015-06-10 | 华侨大学 | Method and device for identifying operating modal parameters of linear time-varying structure and application of the device |
CN106972966A (en) * | 2017-03-28 | 2017-07-21 | 国网辽宁省电力有限公司阜新供电公司 | The real-time detection method of exception of network traffic in a kind of communication network |
CN107991706A (en) * | 2017-11-24 | 2018-05-04 | 华北科技学院 | Combine noise-reduction method based on wavelet packet multiple thresholds and the coal seam hydraulic fracture microseismic signals for improving empirical mode decomposition |
CN107991706B (en) * | 2017-11-24 | 2020-07-31 | 华北科技学院 | Coal seam hydraulic fracturing microseismic signal combined noise reduction method based on wavelet packet multiple threshold and improved empirical mode decomposition |
CN108388848A (en) * | 2018-02-07 | 2018-08-10 | 西安石油大学 | A kind of multiple dimensioned oil-gas-water three-phase flow mechanical characteristic analysis method |
CN108388848B (en) * | 2018-02-07 | 2022-02-22 | 西安石油大学 | Multi-scale oil-gas-water multiphase flow mechanics characteristic analysis method |
CN109282965A (en) * | 2018-11-06 | 2019-01-29 | 中国海洋石油集团有限公司 | It is a kind of to collect the quick identification device of nocuousness flow pattern and method in defeated standpipe |
CN109538944A (en) * | 2018-12-03 | 2019-03-29 | 北京无线电计量测试研究所 | A kind of pipeline leakage detection method |
CN109538944B (en) * | 2018-12-03 | 2020-07-07 | 北京无线电计量测试研究所 | Pipeline leakage detection method |
CN110044533A (en) * | 2019-04-11 | 2019-07-23 | 大连理工大学 | A kind of tools for bolts ' pretension force monitoring methods combined based on EMD Energy-Entropy with support vector machines |
CN110044533B (en) * | 2019-04-11 | 2020-10-20 | 大连理工大学 | Bolt pretightening force monitoring method based on combination of EMD (empirical mode decomposition) energy entropy and support vector machine |
CN111503010A (en) * | 2019-09-18 | 2020-08-07 | 山东建筑大学 | Method for extracting broadband radiation noise signal characteristics of centrifugal pump |
Also Published As
Publication number | Publication date |
---|---|
CN100507509C (en) | 2009-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1904581A (en) | Oil gas water multiphase flow type identification method based on main component analysis and supporting vector machine | |
CN1811367A (en) | Integrated supporting vector machine mixed intelligent diagnosing method for mechanical fault | |
Kafeel et al. | An expert system for rotating machine fault detection using vibration signal analysis | |
CN1904580A (en) | Gas liquid two phase flow type identification method based on information coalescence and flow type signal collection device | |
CN101059796A (en) | Two-stage combined file classification method based on probability subject | |
CN1177407A (en) | Method and system for velocity-based head writing recognition | |
CN101078776A (en) | Nuclear logging data multi-measure processing method | |
CN111965121B (en) | Rapid classification detection method for coal coke slag based on spectral analysis technology | |
CN110702648B (en) | Fluorescent spectrum pollutant classification method based on non-subsampled contourlet transformation | |
CN107169497A (en) | A kind of tumor imaging label extracting method based on gene iconography | |
CN113192076B (en) | MRI brain tumor image segmentation method combining classification prediction and multi-scale feature extraction | |
CN1916940A (en) | Template optimized character recognition method and system | |
CN101034438A (en) | Intelligent locating method face for micro-fluidic chip | |
CN1763740A (en) | Info web piece extracting method and device | |
CN110443303B (en) | Image segmentation and classification-based intelligent identification method for coal-rock microcomponents | |
CN116049637A (en) | Gas-liquid two-phase flow pattern recognition method based on flow noise signals | |
CN1084502C (en) | Method and device for recognition of similar writing | |
CN114694746A (en) | Plant pri-miRNA coding peptide prediction method based on improved MRMD algorithm and DF model | |
CN110609936A (en) | Intelligent classification method for fuzzy address data | |
CN102902984A (en) | Remote-sensing image semi-supervised projection dimension reducing method based on local consistency | |
CN108573285A (en) | A kind of experimental parameter clustering method towards dredging energy consumption | |
CN1214362C (en) | Device and method for determining coretative coefficient between signals and signal sectional distance | |
CN108763289B (en) | Massive heterogeneous sensor format data analysis method | |
CN1563940A (en) | Method of recogniting flowing type of level pipe gas/liquid 2-phase flow based on Hilbert-Huang conversion | |
CN104036527A (en) | Human motion segmentation method based on local linear embedding |
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 |