CN116662894A - Fluid pipeline leakage diagnosis method based on intuitionistic fuzzy weighted least square twin multi-class support vector machine - Google Patents
Fluid pipeline leakage diagnosis method based on intuitionistic fuzzy weighted least square twin multi-class support vector machine Download PDFInfo
- Publication number
- CN116662894A CN116662894A CN202310630259.3A CN202310630259A CN116662894A CN 116662894 A CN116662894 A CN 116662894A CN 202310630259 A CN202310630259 A CN 202310630259A CN 116662894 A CN116662894 A CN 116662894A
- Authority
- CN
- China
- Prior art keywords
- original vibration
- vibration characteristic
- sample
- classification
- samples
- 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
- 238000003745 diagnosis Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000012530 fluid Substances 0.000 title claims abstract description 28
- 238000012706 support-vector machine Methods 0.000 title claims abstract description 19
- 238000013145 classification model Methods 0.000 claims abstract description 38
- 238000005457 optimization Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000017105 transposition Effects 0.000 claims description 6
- 230000006870 function Effects 0.000 description 20
- 230000002159 abnormal effect Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000013136 deep learning model Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 231100000206 health hazard Toxicity 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
- F17D—PIPE-LINE SYSTEMS; PIPE-LINES
- F17D5/00—Protection or supervision of installations
- F17D5/02—Preventing, monitoring, or locating loss
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
- G06N3/0442—Recurrent networks, e.g. Hopfield networks characterised by memory or gating, e.g. long short-term memory [LSTM] or gated recurrent units [GRU]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
The invention relates to a fluid pipeline leakage diagnosis method based on an intuitionistic fuzzy weighted least square twin multi-category support vector machine, which comprises the following steps: obtaining vibration data of the fluid pipeline in a normal state and a leakage state; extracting characteristic parameters of vibration data from two aspects of time domain and waveform to generate an original vibration characteristic sample set; based on an original vibration characteristic sample set, a classification contribution value and an intra-class weight of the original vibration characteristic sample are calculated by using a weighting technology based on k nearest neighbor and an intuitional fuzzy theory, a decision function of an LST-KSVC classification model is reconstructed to generate a pipeline leakage diagnosis model, the original vibration characteristic sample is used as input of the pipeline leakage diagnosis model, the pipeline leakage diagnosis model is trained by using a grid search method with the highest diagnosis precision of the pipeline leakage diagnosis model as an optimization target, and the vibration characteristic sample to be tested is input into the trained pipeline leakage diagnosis model to judge whether leakage occurs in a fluid pipeline.
Description
Technical Field
The invention belongs to the field of pipeline leakage detection, and particularly relates to a fluid pipeline leakage diagnosis method based on an intuitionistic fuzzy weighted least square twin multi-class support vector machine.
Background
The fluid pipeline is an important component of the urban energy transportation system, and is used for transporting important resources such as water, natural gas, petroleum and the like, and the running state of the fluid pipeline is closely related to urban development and resident life. However, due to factors such as out-of-service, chemical corrosion, aging, weld defects and the like, the fluid pipeline in the pipe network system inevitably receives damage to a certain extent, and leakage accidents occur, so that not only is economic direct loss caused, but also potential environmental pollution and public health hazard are caused. Therefore, the suspected leakage hidden danger of the fluid pipeline is diagnosed in time, and the method has important significance for avoiding serious leakage accidents and ensuring the safe operation of the pipeline.
Currently, machine learning methods such as artificial neural networks, deep learning models, support vector machines and the like have been used in the field of pipeline leakage diagnosis. Because the pipeline leakage diagnosis method based on the neural network model needs a large number of samples for model training, and the pipeline is in a normal running state most of the time, and the number of leakage samples is relatively small, the classification methods such as a support vector machine, a twin support vector machine, a least square twin support vector machine and the like which are suitable for small sample learning are more suitable for pipeline leakage diagnosis. The least square twin multi-class support vector machine (Least Squared Twin K-class Support Vector Classification; LST-KSVC) is a novel classification method, and utilizes a one-to-one complementary structure to evaluate all training samples, so that the high training efficiency of the least square twin support vector machine is maintained, and the classification accuracy is higher. Due to the presence of noise disturbances in the pipeline environment, some outliers are unavoidable in the leaky samples. When the LST-KSVC classification model is trained, abnormal values in a training sample can mislead the generation of a classification hyperplane, so that the accuracy rate of pipeline leakage diagnosis is lower.
Disclosure of Invention
In order to solve the problems, the invention provides a fluid pipeline leakage diagnosis method based on an intuitionistic fuzzy weighted least square twin multi-category support vector machine, which considers local field information among data points and improves generalization performance while reducing abnormal value influence by utilizing intuitionistic fuzzy numbers, and specifically comprises the following steps:
s1: obtaining vibration data of a fluid pipeline in a normal state and a leakage state, marking class labels to generate an original vibration sample data set, and extracting time domain characteristics and waveform characteristics of the original vibration sample data from two aspects of time domain and waveform to generate characteristic parameters of the original vibration sample data; normalizing the characteristic parameters of the original vibration sample data to obtain an original vibration characteristic sample set;
s2: calculating classification contribution values of the original vibration characteristic samples by utilizing an intuitional fuzzy theory according to the distance from the original vibration characteristic samples to the class center of the original vibration characteristic samples, the total number of the original vibration characteristic samples and the number of the original vibration characteristic samples under each class; calculating the intra-class weight of the original vibration characteristic sample according to the original vibration characteristic sample and k original vibration characteristic samples adjacent to the original vibration characteristic sample under the category to which the original vibration characteristic sample belongs;
s3: reconstructing a decision function of the LST-KSVC classification model according to the intra-class weight of the original vibration characteristic sample and the classification contribution value of the original vibration characteristic sample to generate a pipeline leakage diagnosis model;
s4: taking an original vibration characteristic sample as input of a pipeline leakage diagnosis model, and training the pipeline leakage diagnosis model by using a grid search method with the highest diagnosis precision of the pipeline leakage diagnosis model as an optimization target;
s5: and acquiring a vibration characteristic sample to be detected of the fluid pipeline, inputting the vibration characteristic sample to be detected into a trained pipeline leakage diagnosis model, and judging whether the fluid pipeline is leaked or not.
Further, the normalizing the characteristic parameters of the original vibration sample data includes:
wherein x is ij A j-th element representing a characteristic parameter of the i-th raw vibration sample data,is x ij Corresponding normalized value, < >>And->The minimum value and the maximum value in the j-th element of the characteristic parameters of all the original vibration sample data are respectively represented.
Further, the calculating the classification contribution value of the original vibration characteristic sample by using the intuitionistic fuzzy theory comprises:
the intuitive ambiguity is defined as θ= (μ) θ ,v θ ) Membership value mu θ ∈[0,1]Non-membership value v θ ∈[0,1]，0≤μ θ +v θ ≤1；
Normalizing the distance from the original vibration characteristic sample to the class center of the original vibration characteristic sample to serve as the membership degree of the original vibration characteristic sample;
taking the ratio of the number of samples under the category to which the original vibration characteristic sample belongs to the total number of samples in the original vibration characteristic sample set as the non-membership degree of the original vibration characteristic sample, wherein θ represents the intuitive fuzzy number of the original vibration characteristic sample; mu (mu) θ Representing membership value of original vibration characteristic sample; v θ A non-membership value representing an original vibration characteristic sample;
calculating a classification contribution value S (theta) of the original vibration characteristic sample according to the membership value and the non-membership value of the original vibration characteristic sample;
further, the intra-class weights of the raw vibration feature samples include:
wherein sigma represents a kernel parameter,intra-class weights, x, representing the ith raw vibration signature sample under the c-th class i Is x j Is denoted as x i Is the same as the jth original vibration characteristic sample x j One of the k nearest samples; x is x i Representing the ith original vibration signature sample, x i And x j The category to which the method belongs is category c.
Further, the decision function of the intuitive fuzzy weighted LST-KSVC classification model comprises:
defining an original vibration feature sample set d= { (x) 1 ,y 1 ),(x 2 ,y 2 ),...,(x i ,y i ),...,(x m ,y m )}，x i ∈R n Representing the original vibration characteristic sample, y i Representing sample x i Category labels of (c); given a sample matrixAnd->The method comprises the steps that the method consists of original vibration characteristic samples corresponding to category labels of +1, -1 and 0, wherein n represents the number of elements in the original vibration characteristic samples;
for a linear intuitive fuzzy weighted LST-KSVC classification model, two classification hyperplanes are constructed:
x T w 1 +b 1 =0 and x T w 2 +b 2 ＝0
Wherein x represents an original vibration characteristic sample, T represents transposition, and w 1 、w 2 And b 1 、b 2 The weights and biases of the two hyperplanes are respectively;
determining two classification hyperplanes according to the intra-class weights of the original vibration characteristic samples and the classification contribution values of the original vibration characteristic samples:
s.t.-(X 2 w 1 +e 2 b 1 )+ξ＝e 2
-(X 3 w 1 +e 3 b 1 )+δ＝e 3 (1-ε)
s.t.(X 1 w 2 +e 1 b 2 )+η＝e 1
-(X 3 w 2 +e 3 b 2 )+γ＝e 3 (1-ε)
wherein ρ is 1 And ρ 2 Representing a matrix consisting of intra-class weights of the original vibration feature samples under the two class labels; s is S 1 And S is 2 Representing a matrix composed of classification contribution values of the original vibration feature samples under the two category labels; ζ, δ, η, γ represent relaxation variables; c (C) 1 、C 2 、C 3 And C 4 Representing penalty parameters; e, e 1 、e 1 And e 3 For dimension m 1 、m 2 And m 3 Is the full 1 vector of (2); epsilon is a positive parameter selected a priori;
the decision function for the linear intuitive fuzzy weighted LST-KSVC classification model is as follows:
where f (x) represents the decision function of the intuitive fuzzy weighted LST-KSVC classification model.
Further, the decision function of the intuitive fuzzy weighted LST-KSVC classification model comprises:
defining an original vibration feature sample set d= { (x) 1 ,y 1 ),(x 2 ,y 2 ),...,(x i ,y i ),...,(x m ,y m )}，x i ∈R n Representing the original vibration characteristic sample, y i Representing sample x i Category labels of (c); given a sample matrixAnd->The method comprises the steps that the method consists of original vibration characteristic samples corresponding to category labels of +1, -1 and 0, wherein n represents the number of elements in the original vibration characteristic samples;
for a nonlinear intuitive fuzzy weighted LST-KSVC classification model, two classification hyperplanes are constructed:
K(x T ,A T )w 1 +b 1 =0 and K (x T ,A T )w 2 +b 2 ＝0
Wherein K (,) represents a kernel function, a= [ X ] 1 ；X 2 ]X represents the original vibration characteristic sample, T represents the transposition, and w 1 、w 2 And b 1 、b 2 Respectively representing the weights and the biases of the two hyperplanes;
determining two classification hyperplanes according to the intra-class weights of the original vibration characteristic samples and the classification contribution values of the original vibration characteristic samples:
s.t.-(X 2 w 1 +e 2 b 1 )+ξ＝e 2
-(X 3 w 1 +e 3 b 1 )+δ＝e 3 (1-ε)
s.t.(X 1 w 2 +e 1 b 2 )+η＝e 1
-(X 3 w 2 +e 3 b 2 )+γ＝e 3 (1-ε)
wherein ρ is 1 And ρ 2 Representing a matrix consisting of intra-class weights of the original vibration feature samples under the two class labels; s is S 1 And S is 2 Representing a matrix composed of classification contribution values of the original vibration feature samples under the two category labels; ζ, δ, η, γ represent relaxation variables; c (C) 1 、C 2 、C 3 And C 4 Representing penalty parameters; e, e 1 、e 1 And e 3 For dimension m 1 、m 2 And m 3 Is the full 1 vector of (2); epsilon is a positive parameter selected a priori; k represents a kernel function; sigma is a kernel parameter.
The decision function for the nonlinear intuitive fuzzy weighted LST-KSVC classification model is as follows:
where f (x) represents the decision function of the intuitive fuzzy weighted LST-KSVC classification model.
Further, the training of the pipeline leakage diagnostic model includes:
for the pipeline leakage diagnosis model to be a linear intuitionistic fuzzy weighting LST-KSVC classification model, inputting an original vibration characteristic sample into the pipeline leakage diagnosis model to determine parameters { C } by a grid search method i Optimal combination of i=1, 2,3,4} solves the optimization problem of the linear intuitionistic fuzzy weighted LST-KSVC model under the optimal parameter combination, thereby determining w 1 ,b 1 And w 2 ,b 2 Is worth training a pipeline leakage diagnosis model;
for the pipeline leakage diagnosis model to be a nonlinear intuitionistic fuzzy weighted LST-KSVC classification model, inputting an original vibration characteristic sample into the pipeline leakage diagnosis model to determine parameters { C } by a grid search method i Optimal combination of i=1, 2,3,4} and σ; determining two optimal non-parallel classification hypersurface K (x T ,C T )w 1 +b 1 =0 and K (x T ,C T )w 2 +b 2 =0 resulted in a trained pipeline leakage diagnostic model.
The invention has at least the following beneficial effects
According to the invention, a membership value and a non-membership value are allocated to each training sample by utilizing the intuitional fuzzy number, so that the classification contribution value of all training samples to the classification model is determined, and the objective function of the LST-KSVC model is improved by utilizing the value, so that the corresponding weight is given to all training samples, and the influence of the abnormal value on the accuracy of the classification model is reduced. In addition, the intuitive fuzzy weighting LST-KSVC model provided by the invention integrates local domain information among data samples based on the weighting technology of k nearest neighbors, so that the generalization performance of the model is improved. The provided intuitive fuzzy weighting LST-KSVC model reduces the influence of abnormal values while retaining higher training efficiency and classification accuracy of the LST-KSVC, and further improves the classification accuracy and generalization performance of the model.
Drawings
Fig. 1 is an overall flow chart of the present invention.
Detailed Description
Other advantages and effects of the present invention will become apparent to those skilled in the art from the following disclosure, which describes the embodiments of the present invention with reference to specific examples. The invention may be practiced or carried out in other embodiments that depart from the specific details, and the details of the present description may be modified or varied from the spirit and scope of the present invention. It should be noted that the illustrations provided in the following embodiments merely illustrate the basic idea of the present invention by way of illustration, and the following embodiments and features in the embodiments may be combined with each other without conflict.
Referring to fig. 1, the present invention provides a fluid pipeline leakage diagnosis method based on an intuitionistic fuzzy weighted least square twin multi-category support vector machine, comprising:
s1: obtaining vibration data of a fluid pipeline in a normal state and a leakage state, marking class labels to generate an original vibration sample data set, and extracting time domain characteristics and waveform characteristics of the original vibration sample data from two aspects of time domain and waveform to generate characteristic parameters of the original vibration sample data; normalizing the characteristic parameters of the original vibration sample data to obtain an original vibration characteristic sample set;
the marking of the category label comprises: vibration data is labeled according to prior knowledge of experts in the field, wherein a leakage state is labeled 1, a normal state is labeled-1, and an abnormal state is labeled 0.
The time domain features include: root mean square, effective value, energy, average value, etc.; the waveform features include: pulse factor, valley degree factor, kurtosis, peak value coefficient, etc.
S2: calculating classification contribution values of the original vibration characteristic samples by utilizing an intuitional fuzzy theory according to the distance from the original vibration characteristic samples to the class center of the original vibration characteristic samples, the total number of the original vibration characteristic samples and the number of the original vibration characteristic samples under each class; calculating the intra-class weight of the original vibration characteristic sample according to k adjacent original vibration characteristic samples under the category to which the original vibration characteristic sample belongs;
the center of the category of the original feature sample is expressed as: if a class of samples is completely evenly distributed within the circle, then ideally the class center of the class is almost at the center of the circle. The center of the actual sample distribution is already obtained by some mature formulas, so that only the center of the category is simply mentioned in the document, and detailed description is omitted.
The k adjacent original vibration feature samples under the category to which the original vibration feature sample belongs are expressed as:
and taking the original vibration characteristic samples as target samples, and taking K original vibration characteristic samples closest to the target samples in the category of the target samples as K adjacent original vibration characteristic samples.
S3: reconstructing a decision function of the LST-KSVC classification model according to the intra-class weight of the original vibration characteristic sample and the classification contribution value of the original vibration characteristic sample to generate a pipeline leakage diagnosis model;
s4: taking an original vibration characteristic sample as input of a pipeline leakage diagnosis model, and training the pipeline leakage diagnosis model by using a grid search method with the highest diagnosis precision of the pipeline leakage diagnosis model as an optimization target;
s5: and acquiring a vibration characteristic sample to be detected of the fluid pipeline, inputting the vibration characteristic sample to be detected into a trained pipeline leakage diagnosis model, and judging whether the fluid pipeline is leaked or not.
Further, the normalizing the characteristic parameters of the original vibration sample data includes:
wherein x is ij Representing the ith raw vibration sample dataIs a j-th element of the characteristic parameter of (c),is x ij Corresponding normalized value, < >>And->The minimum value and the maximum value in the j-th element of the characteristic parameters of all the original vibration sample data are respectively represented.
Further, the calculating the classification contribution value of the original vibration characteristic sample by using the intuitionistic fuzzy theory comprises:
the intuitive ambiguity is defined as θ= (μ) θ ,v θ ) Membership value mu θ ∈[0,1]Non-membership value v θ ∈[0,1]，0≤μ θ +v θ ≤1；
Normalizing the distance from the original vibration characteristic sample to the class center of the original vibration characteristic sample to serve as the membership degree of the original vibration characteristic sample;
taking the ratio of the number of samples under the category to which the original vibration characteristic sample belongs to the total number of samples in the original vibration characteristic sample set as the non-membership degree of the original vibration characteristic sample, wherein θ represents the intuitive fuzzy number of the original vibration characteristic sample; mu (mu) θ Representing membership value of original vibration characteristic sample; v θ A non-membership value representing an original vibration characteristic sample;
calculating a classification contribution value S (theta) of the original vibration characteristic sample according to the membership value and the non-membership value of the original vibration characteristic sample;
further, the intra-class weights of the raw vibration feature samples include:
wherein sigma represents a kernel parameter,intra-class weights, x, representing the ith raw vibration signature sample under the c-th class i Is x j Is denoted as x i Is the same as the jth original vibration characteristic sample x j One of the k nearest points; x is x i Representing the ith original vibration signature sample, x i And x j The category to which the method belongs is category c.
Further, the decision function of the intuitive fuzzy weighted LST-KSVC classification model comprises:
defining an original vibration feature sample set d= { (x) 1 ,y 1 ),(x 2 ,y 2 ),...,(x i ,y i ),...,(x m ,y m )}，x i ∈R n Representing the original vibration characteristic sample, y i Representing sample x i Category labels of (c); given a sample matrixAnd->The vibration detection method comprises the steps of respectively forming original vibration characteristic samples corresponding to category labels of "+1", "-1" and "0", wherein n represents the number of elements in the original vibration characteristic samples;
for a linear intuitive fuzzy weighted LST-KSVC classification model, two classification hyperplanes are constructed:
x T w 1 +b 1 =0 and x T w 2 +b 2 ＝0
Wherein x represents an original vibration characteristic sample, T represents transposition, and w 1 、w 2 And b 1 、b 2 The weights and biases of the two hyperplanes are respectively;
determining two classification hyperplanes according to the intra-class weights of the original vibration characteristic samples and the classification contribution values of the original vibration characteristic samples:
s.t.-(X 2 w 1 +e 2 b 1 )+ξ＝e 2
-(X 3 w 1 +e 3 b 1 )+δ＝e 3 (1-ε)
s.t.(X 1 w 2 +e 1 b 2 )+η＝e 1
-(X 3 w 2 +e 3 b 2 )+γ＝e 3 (1-ε)
wherein ρ is 1 And ρ 2 Representing a matrix consisting of intra-class weights of the original vibration feature samples under the two class labels; s is S 1 And S is 2 Representing a matrix composed of classification contribution values of the original vibration feature samples under the two category labels; ζ, δ, η, γ represent relaxation variables; c (C) 1 、C 2 、C 3 And C 4 Representing penalty parameters; e, e 1 、e 1 And e 3 For dimension m 1 、m 2 And m 3 Is the full 1 vector of (2); epsilon is a positive parameter selected a priori;
the decision function for the linear intuitive fuzzy weighted LST-KSVC classification model is as follows:
where f (x) represents the decision function of the intuitive fuzzy weighted LST-KSVC classification model.
Further, the decision function of the intuitive fuzzy weighted LST-KSVC classification model comprises:
defining an original vibration feature sample set d= { (x) 1 ,y 1 ),(x 2 ,y 2 ),...,(x i ,y i ),...,(x m ,y m )}，x i ∈R n Representing the original vibration characteristic sample, y i Representing sample x i Category labels of (c); given a sample matrixAnd->The method comprises the steps that the method consists of original vibration characteristic samples corresponding to category labels of +1, -1 and 0, wherein n represents the number of elements in the original vibration characteristic samples;
for a nonlinear intuitive fuzzy weighted LST-KSVC classification model, two classification hyperplanes are constructed:
K(x T ,A T )w 1 +b 1 =0 and K (x T ,A T )w 2 +b 2 ＝0
Wherein K (,) represents a kernel function, a= [ X ] 1 ；X 2 ]X represents the original vibration characteristic sample, T represents the transposition, and w 1 、w 2 And b 1 、b 2 Respectively representing the weights and the biases of the two hyperplanes;
determining two classification hyperplanes according to the intra-class weights of the original vibration characteristic samples and the classification contribution values of the original vibration characteristic samples:
s.t.-(X 2 w 1 +e 2 b 1 )+ξ＝e 2
-(X 3 w 1 +e 3 b 1 )+δ＝e 3 (1-ε)
s.t.(X 1 w 2 +e 1 b 2 )+η＝e 1
-(X 3 w 2 +e 3 b 2 )+γ＝e 3 (1-ε)
wherein ρ is 1 And ρ 2 Representing a matrix consisting of intra-class weights of the original vibration feature samples under the two class labels; s is S 1 And S is 2 Representing a matrix composed of classification contribution values of the original vibration feature samples under the two category labels; ζ, δ, η, γ represent relaxation variables; c (C) 1 、C 2 、C 3 And C 4 Representing penalty parameters; e, e 1 、e 1 And e 3 For dimension m 1 、m 2 And m 3 Is the full 1 vector of (2); epsilon is a positive parameter selected a priori; k represents a kernel function; sigma is a kernel parameter.
The decision function for the nonlinear intuitive fuzzy weighted LST-KSVC classification model is as follows:
where f (x) represents the decision function of the intuitive fuzzy weighted LST-KSVC classification model.
Further, the training of the pipeline leakage diagnostic model includes:
for the pipeline leakage diagnosis model to be a linear intuitionistic fuzzy weighting LST-KSVC classification model, inputting an original vibration characteristic sample into the pipeline leakage diagnosis model to determine parameters { C } by a grid search method i Optimal combination of i=1, 2,3,4} solves the optimization problem of the linear intuitionistic fuzzy weighted LST-KSVC model under the optimal parameter combination, thereby determining w 1 ,b 1 And w 2 ,b 2 Is worth training a pipeline leakage diagnosis model;
LST-KSVC classification for nonlinear intuitionistic fuzzy weighting of pipeline leakage diagnostic modelsThe model is characterized in that an original vibration characteristic sample is input into a pipeline leakage diagnosis model, and parameters { C } are determined through a grid search method i Optimal combination of i=1, 2,3,4} and σ; determining two optimal non-parallel classification hypersurface K (x T ,C T )w 1 +b 1 =0 and K (x T ,C T )w 2 +b 2 =0 resulted in a trained pipeline leakage diagnostic model.
According to the invention, a membership value and a non-membership value are allocated to each training sample by utilizing the intuitional fuzzy number, so that the classification contribution value of all training samples to the classification model is determined, and the objective function of the LST-KSVC model is improved by utilizing the value, so that the corresponding weight is given to all training samples, and the influence of the abnormal value on the accuracy of the classification model is reduced. In addition, the intuitive fuzzy weighting LST-KSVC model provided by the invention integrates local domain information among data samples based on the weighting technology of k nearest neighbors, so that the generalization performance of the model is improved. The provided intuitive fuzzy weighting LST-KSVC model reduces the influence of abnormal values while retaining higher training efficiency and classification accuracy of the LST-KSVC, and further improves the classification accuracy and generalization performance of the model.
Finally, it is noted that the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications and equivalents may be made thereto without departing from the spirit and scope of the present invention, which is intended to be covered by the claims of the present invention.
Claims (7)
1. A fluid conduit leak diagnosis method based on an intuitive fuzzy weighted least squares twinning multi-class support vector machine, comprising:
s1: obtaining vibration data of a fluid pipeline in a normal state and a leakage state, marking class labels to generate an original vibration sample data set, and extracting time domain characteristics and waveform characteristics of the original vibration sample data from two aspects of time domain and waveform to generate characteristic parameters of the original vibration sample data; normalizing the characteristic parameters of the original vibration sample data to obtain an original vibration characteristic sample set;
s2: calculating classification contribution values of the original vibration characteristic samples by utilizing an intuitional fuzzy theory according to the distance from the original vibration characteristic samples to the class center of the original vibration characteristic samples, the total number of the original vibration characteristic samples and the number of the original vibration characteristic samples under each class; calculating the intra-class weight of the original vibration characteristic sample according to the original vibration characteristic sample and k original vibration characteristic samples adjacent to the original vibration characteristic sample under the category to which the original vibration characteristic sample belongs;
s3: reconstructing a decision function of the LST-KSVC classification model according to the intra-class weight of the original vibration characteristic sample and the classification contribution value of the original vibration characteristic sample to generate a pipeline leakage diagnosis model;
s4: taking an original vibration characteristic sample as input of a pipeline leakage diagnosis model, and training the pipeline leakage diagnosis model by using a grid search method with the highest diagnosis precision of the pipeline leakage diagnosis model as an optimization target;
s5: and acquiring a vibration characteristic sample to be detected of the fluid pipeline, inputting the vibration characteristic sample to be detected into a trained pipeline leakage diagnosis model, and judging whether the fluid pipeline is leaked or not.
2. The fluid pipeline leakage diagnosis method based on the intuitionistic fuzzy weighted least square twin multi-class support vector machine according to claim 1, wherein the normalizing the characteristic parameters of the original vibration sample data comprises:
wherein x is ij A j-th element representing a characteristic parameter of the i-th raw vibration sample data,is x ij Corresponding normalized value, < >>And->The minimum value and the maximum value in the j-th element of the characteristic parameters of all the original vibration sample data are respectively represented.
3. The fluid conduit leakage diagnosis method based on the intuitionistic fuzzy weighted least square twin multi-class support vector machine according to claim 1, wherein the calculating the classification contribution value of the original vibration characteristic sample by using the intuitionistic fuzzy theory comprises:
the intuitive ambiguity is defined as θ= (μ) θ ,v θ ) Membership value mu θ ∈[0,1]Non-membership value v θ ∈[0,1]，0≤μ θ +v θ ≤1；
Normalizing the distance from the original vibration characteristic sample to the class center of the original vibration characteristic sample to serve as the membership degree of the original vibration characteristic sample;
taking the ratio of the number of samples under the category to which the original vibration characteristic sample belongs to the total number of samples in the original vibration characteristic sample set as the non-membership degree of the original vibration characteristic sample, wherein θ represents the intuitive fuzzy number of the original vibration characteristic sample; mu (mu) θ Representing membership value of original vibration characteristic sample; v θ A non-membership value representing an original vibration characteristic sample;
calculating a classification contribution value S (theta) of the original vibration characteristic sample according to the membership value and the non-membership value of the original vibration characteristic sample;
4. the method for fluid conduit leakage diagnosis based on an intuitionistic fuzzy weighted least squares twinning multi-class support vector machine of claim 1, wherein the intra-class weights of the raw vibration signature samples comprise:
wherein sigma represents a kernel parameter,intra-class weights, x, representing the ith raw vibration signature sample under the c-th class i Is x j Is denoted as x i Is the same as the jth original vibration characteristic sample x j One of the k nearest samples; x is x i Representing the ith original vibration signature sample, x i And x j The category to which the method belongs is category c.
5. The method for fluid conduit leakage diagnosis based on an intuitive fuzzy weighted least squares twinning multi-class support vector machine of claim 1, wherein the decision function of the intuitive fuzzy weighted LST-KSVC classification model comprises:
defining an original vibration feature sample set d= { (x) 1 ,y 1 ),(x 2 ,y 2 ),...,(x i ,y i ),...,(x m ,y m )}，x i ∈R n Representing the original vibration characteristic sample, y i Representing sample x i Category labels of (c); given a sample matrixAndrespectively consists of original vibration characteristic samples corresponding to category labels of "+1", "-1" and "0", wherein n represents the number of elements in the original vibration characteristic samplesAn amount of;
for a linear intuitive fuzzy weighted LST-KSVC classification model, two classification hyperplanes are constructed:
x T w 1 +b 1 =0 and x T w 2 +b 2 ＝0
Wherein x represents an original vibration characteristic sample, T represents transposition, and w 1 、w 2 And b 1 、b 2 The weights and biases of the two hyperplanes are respectively;
determining two classification hyperplanes according to the intra-class weights of the original vibration characteristic samples and the classification contribution values of the original vibration characteristic samples:
s.t.-(X 2 w 1 +e 2 b 1 )+ξ＝e 2
-(X 3 w 1 +e 3 b 1 )+δ＝e 3 (1-ε)
s.t.(X 1 w 2 +e 1 b 2 )+η＝e 1
-(X 3 w 2 +e 3 b 2 )+γ＝e 3 (1-ε)
wherein ρ is 1 And ρ 2 Representing a matrix consisting of intra-class weights of the original vibration feature samples under the two class labels; s is S 1 And S is 2 Representing a matrix composed of classification contribution values of the original vibration feature samples under the two category labels; ζ, δ, η, γ represent relaxation variables; c (C) 1 、C 2 、C 3 And C 4 Representing penalty parameters; e, e 1 、e 1 And e 3 For dimension m 1 、m 2 And m 3 Is the full 1 vector of (2); epsilon is a positive parameter selected a priori;
the decision function for the linear intuitive fuzzy weighted LST-KSVC classification model is as follows:
where f (x) represents the decision function of the intuitive fuzzy weighted LST-KSVC classification model.
6. The method for fluid conduit leakage diagnosis based on an intuitive fuzzy weighted least squares twinning multi-class support vector machine of claim 1, wherein the decision function of the intuitive fuzzy weighted LST-KSVC classification model comprises:
defining an original vibration feature sample set d= { (x) 1 ,y 1 ),(x 2 ,y 2 ),...,(x i ,y i ),...,(x m ,y m )}，x i ∈R n Representing the original vibration characteristic sample, y i Representing sample x i Category labels of (c); given a sample matrixAndthe method comprises the steps that the method consists of original vibration characteristic samples corresponding to category labels of +1, -1 and 0, wherein n represents the number of elements in the original vibration characteristic samples;
for a nonlinear intuitive fuzzy weighted LST-KSVC classification model, two classification hyperplanes are constructed:
K(x T ,A T )w 1 +b 1 =0 and K (x T ,A T )w 2 +b 2 ＝0
Wherein K (,) represents a kernel function, a= [ X ] 1 ；X 2 ]X represents the original vibration characteristic sample, T represents the transposition, and w 1 、w 2 And b 1 、b 2 Respectively representing the weights and the biases of the two hyperplanes;
determining two classification hyperplanes according to the intra-class weights of the original vibration characteristic samples and the classification contribution values of the original vibration characteristic samples:
s.t.-(X 2 w 1 +e 2 b 1 )+ξ＝e 2
-(X 3 w 1 +e 3 b 1 )+δ＝e 3 (1-ε)
s.t.(X 1 w 2 +e 1 b 2 )+η＝e 1
-(X 3 w 2 +e 3 b 2 )+γ＝e 3 (1-ε)
wherein ρ is 1 And ρ 2 Representing a matrix consisting of intra-class weights of the original vibration feature samples under the two class labels; s is S 1 And S is 2 Representing a matrix composed of classification contribution values of the original vibration feature samples under the two category labels; ζ, δ, η, γ represent relaxation variables; c (C) 1 、C 2 、C 3 And C 4 Representing penalty parameters; e, e 1 、e 1 And e 3 For dimension m 1 、m 2 And m 3 Is the full 1 vector of (2); epsilon is a positive parameter selected a priori; k represents a kernel function; sigma is a kernel parameter;
the decision function for the nonlinear intuitive fuzzy weighted LST-KSVC classification model is as follows:
where f (x) represents the decision function of the intuitive fuzzy weighted LST-KSVC classification model.
7. The fluid pipeline leakage diagnosis method based on the intuitive fuzzy weighted least square twin multi-class support vector machine according to claim 5 or 6, wherein the training of the pipeline leakage diagnosis model comprises:
for the pipeline leakage diagnosis model to be a linear intuitionistic fuzzy weighting LST-KSVC classification model, inputting an original vibration characteristic sample into the pipeline leakage diagnosis model to determine parameters (C) by a grid search method i Optimal combination of i=1, 2,3,4} solves the optimization problem of the linear intuitionistic fuzzy weighted LST-KSVC model under the optimal parameter combination, thereby determining w 1 ,b 1 And w 2 ,b 2 Is worth training a pipeline leakage diagnosis model;
for the pipeline leakage diagnosis model to be a nonlinear intuitionistic fuzzy weighted LST-KSVC classification model, inputting an original vibration characteristic sample into the pipeline leakage diagnosis model to determine parameters { C } by a grid search method i Optimal combination of i=1, 2,3,4} and σ; determining two optimal non-parallel classification hypersurface K (x T ,C T )w 1 +b 1 =0 and K (x T ,C T )w 2 +b 2 =0 resulted in a trained pipeline leakage diagnostic model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310630259.3A CN116662894B (en) | 2023-05-31 | 2023-05-31 | Fluid pipeline leakage diagnosis method based on intuitionistic fuzzy weighted least square twin multi-class support vector machine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310630259.3A CN116662894B (en) | 2023-05-31 | 2023-05-31 | Fluid pipeline leakage diagnosis method based on intuitionistic fuzzy weighted least square twin multi-class support vector machine |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116662894A true CN116662894A (en) | 2023-08-29 |
| CN116662894B CN116662894B (en) | 2024-03-26 |
Family
ID=87723687
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310630259.3A Active CN116662894B (en) | 2023-05-31 | 2023-05-31 | Fluid pipeline leakage diagnosis method based on intuitionistic fuzzy weighted least square twin multi-class support vector machine |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116662894B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103472867A (en) * | 2013-09-22 | 2013-12-25 | 浙江大学 | Pesticide production waste liquid incinerator temperature optimization system and method of support vector machine |
| CN104502103A (en) * | 2014-12-07 | 2015-04-08 | 北京工业大学 | Bearing fault diagnosis method based on fuzzy support vector machine |
| WO2016061586A1 (en) * | 2014-10-17 | 2016-04-21 | Cireca Theranostics, Llc | Methods and systems for classifying biological samples, including optimization of analyses and use of correlation |
| CN105572572A (en) * | 2015-12-19 | 2016-05-11 | 渤海大学 | WKNN-LSSVM-based analog circuit fault diagnosis method |
| US11354599B1 (en) * | 2021-06-25 | 2022-06-07 | BRYTER GmbH | Methods and systems for generating a data structure using graphical models |
| CN115690677A (en) * | 2022-10-26 | 2023-02-03 | 浙江工业大学 | Abnormal behavior identification method based on local intuition fuzzy support vector machine |
| CN115978466A (en) * | 2022-12-27 | 2023-04-18 | 重庆市荣冠科技有限公司 | Fluid pipeline leakage detection method based on class imbalance improved twin support vector machine |
-
2023
- 2023-05-31 CN CN202310630259.3A patent/CN116662894B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103472867A (en) * | 2013-09-22 | 2013-12-25 | 浙江大学 | Pesticide production waste liquid incinerator temperature optimization system and method of support vector machine |
| WO2016061586A1 (en) * | 2014-10-17 | 2016-04-21 | Cireca Theranostics, Llc | Methods and systems for classifying biological samples, including optimization of analyses and use of correlation |
| CN104502103A (en) * | 2014-12-07 | 2015-04-08 | 北京工业大学 | Bearing fault diagnosis method based on fuzzy support vector machine |
| CN105572572A (en) * | 2015-12-19 | 2016-05-11 | 渤海大学 | WKNN-LSSVM-based analog circuit fault diagnosis method |
| US11354599B1 (en) * | 2021-06-25 | 2022-06-07 | BRYTER GmbH | Methods and systems for generating a data structure using graphical models |
| CN115690677A (en) * | 2022-10-26 | 2023-02-03 | 浙江工业大学 | Abnormal behavior identification method based on local intuition fuzzy support vector machine |
| CN115978466A (en) * | 2022-12-27 | 2023-04-18 | 重庆市荣冠科技有限公司 | Fluid pipeline leakage detection method based on class imbalance improved twin support vector machine |
Non-Patent Citations (2)
| Title |
|---|
| 安悦瑄等: ""孪生支持向量机综述"", 《计算机科学》, vol. 45, no. 11 * |
| 张丹: ""一种新的直觉模糊最小二乘支持向量机"", 《西安电子科技大学学报》, vol. 49, no. 5 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116662894B (en) | 2024-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111709448B (en) | Mechanical fault diagnosis method based on migration relation network | |
| CN112016251B (en) | Nuclear power device fault diagnosis method and system | |
| Kang et al. | Method of assessing the state of a rolling bearing based on the relative compensation distance of multiple-domain features and locally linear embedding | |
| CN106874957A (en) | A kind of Fault Diagnosis of Roller Bearings | |
| CN102944418A (en) | Wind turbine generator group blade fault diagnosis method | |
| CN110110768A (en) | Fault Diagnosis of Roller Bearings based on Concurrent Feature study and multi-categorizer | |
| CN108596260A (en) | A kind of grid leakage loss localization method and device | |
| CN110595778B (en) | Wind turbine generator bearing fault diagnosis method based on MMF and IGRA | |
| CN112507479B (en) | Oil drilling machine health state assessment method based on manifold learning and softmax | |
| CN110299192A (en) | A kind of environmental suitability evaluation method of firearm components and its composite material, high molecular material | |
| CN108171432B (en) | Ecological risk evaluation method based on multi-dimensional cloud model-fuzzy support vector machine | |
| Zhang et al. | A hybrid forecasting system with complexity identification and improved optimization for short-term wind speed prediction | |
| CN111160457A (en) | Turboshaft engine fault detection method based on soft class extreme learning machine | |
| CN116542170A (en) | Drainage pipeline siltation disease dynamic diagnosis method based on SSAE and MLSTM | |
| CN115978466A (en) | Fluid pipeline leakage detection method based on class imbalance improved twin support vector machine | |
| CN111695288A (en) | Transformer fault diagnosis method based on Apriori-BP algorithm | |
| CN116662894B (en) | Fluid pipeline leakage diagnosis method based on intuitionistic fuzzy weighted least square twin multi-class support vector machine | |
| Jun et al. | Convolutional neural network for burst detection in smart water distribution systems | |
| CN114444663A (en) | Water supply pipe network leakage detection and positioning method based on time convolution network | |
| CN114021610A (en) | Fan fault recognition model training method and system based on transfer learning | |
| Shan et al. | Semisupervised fault diagnosis of gearbox using weighted graph-based label propagation and virtual adversarial training | |
| CN103400021A (en) | KDSNE hydroelectric generating set noise source detection method based on heavy tailed distribution | |
| Dang et al. | seq2graph: discovering dynamic dependencies from multivariate time series with multi-level attention | |
| CN115165366A (en) | Variable working condition fault diagnosis method and system for rotary machine | |
| Rituraj et al. | A Comprehensive Investigation into the Application of Convolutional Neural Networks (ConvNet/CNN) in Smart Grids |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |