CN116610907A - Gear vibration signal characteristic extraction method based on variational modal decomposition - Google Patents
Gear vibration signal characteristic extraction method based on variational modal decomposition Download PDFInfo
- Publication number
- CN116610907A CN116610907A CN202310612456.2A CN202310612456A CN116610907A CN 116610907 A CN116610907 A CN 116610907A CN 202310612456 A CN202310612456 A CN 202310612456A CN 116610907 A CN116610907 A CN 116610907A
- Authority
- CN
- China
- Prior art keywords
- value
- signal
- matrix
- vibration signal
- singular value
- 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
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 24
- 238000000605 extraction Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 38
- 238000001228 spectrum Methods 0.000 claims abstract description 19
- 208000036075 Autosomal dominant tubulointerstitial kidney disease Diseases 0.000 claims abstract description 17
- 101001133056 Homo sapiens Mucin-1 Proteins 0.000 claims abstract description 17
- 102100034256 Mucin-1 Human genes 0.000 claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 13
- 230000035772 mutation Effects 0.000 claims abstract description 13
- 241000254173 Coleoptera Species 0.000 claims abstract description 8
- 210000003608 fece Anatomy 0.000 claims abstract description 8
- 230000004927 fusion Effects 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000002708 enhancing effect Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 claims 1
- 238000012216 screening Methods 0.000 abstract description 3
- 239000013598 vector Substances 0.000 description 14
- 238000003745 diagnosis Methods 0.000 description 3
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
-
- 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
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
A gear vibration signal characteristic extraction method based on variation modal decomposition relates to a gear vibration signal characteristic extraction method. The method comprises the following steps: step one, signal acquisition; step two, constructing a Hankel matrix of the signal, calculating a singular value kurtosis differential spectrum of the Hankel matrix, and positioning a maximum mutation point r of the differential spectrum; step three, signal denoising is carried out according to the effective order of the r positioning singular value matrix; calculating the decomposition layer number K of the VMD by using the proposed optimal K value algorithm; fifthly, screening IMF components by using a correlation and peak value principle; step six, optimizing MCKD parameters by using a dung beetle algorithm, and carrying out pulse enhancement on signals by using the optimized MCKD; and seventhly, feature extraction is achieved by using fuzzy entropy, approximate entropy and sample entropy. The invention provides a gear vibration signal characteristic extraction method which can extract the characteristics of gear vibration signals more completely. The method is suitable for extracting the characteristics of the gear vibration signals.
Description
Technical field:
the invention relates to a vibration signal characteristic extraction method.
Background field:
when the state of the gear is evaluated, because the environmental factors often cause that other interference signals are contained in the obtained gear vibration signals, some weak fault signals can be covered in noise signals, so that fault information contained in the gear vibration signals is not obvious, further, the feature extraction of the gear fault information is incomplete, and the fault diagnosis after the fault diagnosis is affected.
The invention comprises the following steps:
in order to reduce the influence of noise interference on gear fault information feature extraction and improve the accuracy of gear state evaluation, the invention provides a gear vibration signal feature extraction method based on variation modal decomposition.
The gear vibration signal characteristic extraction based on variation modal decomposition is carried out according to the following steps:
step 1, collecting gear vibration signals;
step 2, constructing a Hankel matrix of the signal, calculating a singular value kurtosis difference spectrum of the Hankel matrix, and positioning a maximum mutation position r of the difference spectrum;
step 3, reserving the first r singular values in the Hankel matrix, setting zero for other singular values, and carrying out signal denoising;
step 4, calculating the VMD decomposition layer number by using an optimal K value algorithm;
step 5, using the correlation and peak value index as the basis to screen IMF components and reconstructing signals;
step 6, selecting the optimal parameters in the MCKD by using a dung beetle optimization algorithm, and carrying out pulse enhancement on the reconstructed signal by using the optimized MCKD;
and 7, extracting features by using fuzzy entropy, approximate entropy and sample entropy to obtain a fusion feature matrix.
Further, in the second step, discretizing the collected continuous signals and constructing a Hankel matrix, wherein the matrix expression is as follows:
in formula (1), x (i) is an element in the vibration signal sequence;
further, in the second step, a singular value decomposition Hankel matrix is used and a singular value differential spectrum is calculated, so that a maximum mutation r of the singular value differential spectrum is obtained, and a singular value decomposition expression is as follows:
A=UΣV T (2)
in formula (2), U is the left singular vector of matrix a, Σ is the singular value of matrix a, and V is the right singular vector of matrix a.
Further, the calculation method of the optimal K value algorithm in the fourth step is as follows: singular value decomposition is carried out on a Hankel matrix constructed by the original signal, a singular value kurtosis difference spectrum of the Hankel matrix is calculated, a kurtosis difference spectrum mutation position is located, and the mutation position is used as an initial K value. And decomposing the original vibration signal by using an initial K value, calculating the center frequency distance between adjacent IMF components, adding 1 to the initial K value if all the initial K values are larger than 0.1, and stopping until the distance between any adjacent IMFs is smaller than 0.1, wherein the K value is the optimal K value of the VMD.
Further, the IMF screening method in the fifth step comprises the following steps: and calculating the correlation between the IMF and the original signal and the kurtosis value of the IMF, selecting the IMF with the correlation coefficient larger than the average value and the kurtosis value larger than 3, and carrying out signal reconstruction by using the screened IMF.
Further, the fusion mode of the feature vectors in the step six is as follows: dividing the reconstruction signal into n segments, calculating fuzzy entropy, approximate entropy and sample entropy of each segment, taking entropy value as characteristic of the reconstruction signal, fusing entropy values of all signal segments, for example, the characteristic vector of A is A= [ a ] 1 … a n ]The feature vector of B is b= [ B ] 1 … b n ]The C eigenvector is c= [ C 1 … c n ]So the fusion feature vector is f= [ a ] 1 … a n b 1 … b n c 1 … c n ]。
The principle of the invention is as follows:
the method comprises the steps of constructing a Hankel matrix of a vibration signal during feature extraction, decomposing the Hankel matrix of the vibration signal by singular values, and calculating a singular value kurtosis difference spectrum; positioning the maximum mutation position r of kurtosis spectrum, setting all elements after the r order of the diagonal line of the singular value matrix to zero, and using the singular value matrix to reconstruct vibration signals so as to finish signal denoising; calculating a preset mode number K of the VMD by using a proposed optimal K value algorithm; screening IMF by using cross correlation and kurtosis as indexes to reconstruct signals; introducing a dung beetle algorithm to optimize parameters L and M in the MCKD, and performing pulse enhancement on a reconstruction signal by using the optimized MCKD; features are extracted from different angles by using fuzzy entropy, approximate entropy and sample entropy, and fusion feature vectors are formed.
The invention has the beneficial effects that:
according to the invention, the mutation position of the singular value kurtosis difference spectrum is obtained according to the Hankel matrix of the original signal, the SVD is used for removing noise in the original signal, the proposed optimal K value algorithm is used for calculating the preset mode number K in the variation mode decomposition algorithm, and the problem that the setting of the preset mode number K depends on experience judgment in the past is solved. The invention carries out signal reconstruction according to the correlation and kurtosis value, and screens useful IMF component signals. According to the invention, an optimal combination of parameters L and M in the MCKD is found by using a dung beetle optimization algorithm, and pulse components in the reconstructed signal are enhanced by using the optimized MCKD, so that weak features can be identified during feature extraction. The invention uses three modes to extract the features from different angles to manufacture the fusion feature vector, thereby solving the problem of incomplete feature extraction caused by the single mode of feature extraction and ensuring the quality of the extracted fault features.
Description of the drawings (last):
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a reconstructed signal time domain diagram;
fig. 3 is a signal feature fusion graph.
The specific embodiment is as follows:
the invention will be further described in detail with reference to the drawings and examples below for better explaining the objects, technical solutions and advantages of the invention. The specific embodiments described herein are merely illustrative of the invention and are not intended to be limiting.
As described in connection with fig. 2 and 3, the processed reconstructed signal is more distinct than the original signal characteristics.
The gear vibration signal characteristic extraction method based on variation modal decomposition in the embodiment is carried out according to the following steps:
step one: a gearbox system built in a laboratory collects signals;
step two: constructing a vibration signal Hankel matrix, calculating singular values of the Hankel matrix by singular value decomposition, and locating a maximum mutation r by a singular value peak differential spectrum, wherein the Hankel matrix and the singular value decomposition expression are as follows:
A=UΣV T (2)
in formula (2), U is the left singular vector of matrix a, Σ is the singular value of matrix a, and V is the right singular vector of matrix a.
Step three: performing signal denoising processing by using SVD according to r;
the signal denoising method in the third step comprises the following steps: setting elements after r orders on diagonal lines in the singular value matrix decomposed by SVD to zero, and carrying out signal reconstruction by using the singular value matrix after processing.
Step four: calculating the VMD decomposition layer number by using an optimal K value algorithm;
the optimal K value algorithm in the fourth step is as follows: setting the maximum mutation r of the original signal singular value peak differential spectrum as an initial K value, carrying out VMD decomposition on the denoising signal by using the initial K value, calculating the center frequency distance between adjacent IMFs, adding 1 to the K value when all the center frequency distances are larger than 0.1, and repeating the steps until the step is finished when the center frequency distance between any adjacent IMFs is smaller than 0.1, and obtaining the optimal K value of the VMD.
Step five: using a correlation and kurtosis principle to screen IMF components with correlation larger than the mean value and kurtosis larger than 3 for signal reconstruction;
step six: optimizing MCKD parameters by using a dung beetle algorithm and carrying out pulse enhancement on the reconstructed signals;
the MCKD parameter optimization algorithm in the step six is as follows: randomly initializing 100L, M combinations, using the sample entropy of the MCKD enhanced signal as an objective function, finding the L, M combination that minimizes the sample entropy of the MCKD enhanced result over 15 iterations, and pulse enhancing the reconstructed signal using set L, M.
Step seven: extracting a reconstructed signal characteristic by using a sample entropy, a fuzzy entropy and an approximate entropy to obtain a fusion characteristic;
the method comprises the steps of carrying out equidistant segmentation on a reconstructed signal, dividing the reconstructed signal into n sections, carrying out feature extraction on each section by using sample entropy, fuzzy entropy and approximate entropy, and splicing the extracted features to obtain a fusion feature vector.
Experimental results show that noise in an original signal can be removed by SVD denoising, the optimal K value algorithm can obtain reasonable preset mode numbers, the dung beetle optimization algorithm can effectively jump out the limitation of local optimization and find the optimal MCKD parameter combination, the optimized MCKD enhances weak pulse signals in a reconstructed signal, feature line vectors are extracted in different modes by using fuzzy entropy, approximate entropy and sample entropy and combined into a mixed feature vector matrix, features are extracted from different scales, the obtained feature vectors are more complete, and support is provided for gear fault diagnosis.
Claims (7)
1. The gear vibration signal characteristic extraction method based on variation modal decomposition is characterized by comprising the following specific processes:
step 1, a laboratory collects vibration signals of a gear box;
step 2, constructing a Hankel matrix of the signal, calculating a singular value kurtosis difference spectrum of the Hankel matrix, and positioning a maximum mutation position r of the difference spectrum;
step 3, reserving the first r singular values in the Hankel matrix, setting zero for other singular values, and carrying out signal denoising;
step 4, calculating the VMD decomposition layer number by using an optimal K value algorithm;
step 5, using the correlation and peak value index as the basis to screen IMF components and reconstructing signals;
step 6, selecting the optimal parameters in the MCKD by using a dung beetle optimization algorithm, and carrying out pulse enhancement on the reconstructed signal by using the optimized MCKD;
and 7, extracting features by using fuzzy entropy, approximate entropy and sample entropy to obtain a fusion feature matrix.
2. The method for extracting the gear vibration signal characteristics of the variation modal decomposition according to claim 1, wherein in the step 2, the vibration signal is constructed as a Hankel matrix, a Hankel singular value matrix is obtained by using a singular value decomposition algorithm, a singular value kurtosis differential spectrum is calculated, and a maximum mutation r of the differential spectrum is positioned.
3. The method for extracting the gear vibration signal characteristics by the variation modal decomposition according to claim 1, wherein in the step 3, the effective order of a singular value matrix is determined through the mutation position of the singular value kurtosis spectrum, the singular value after the r order of a diagonal line of the singular value matrix is set to zero, and the processed singular value matrix is used for obtaining the vibration signal after denoising.
4. The method for extracting the gear vibration signal characteristics of the variation modal decomposition according to claim 1, wherein step 4 provides an optimal K value calculation method of VMD: taking the maximum mutation position of the singular value kurtosis difference spectrum of the signal as an initial K value, carrying out VMD decomposition on the signal, calculating the center frequency distance between the decomposed adjacent IMF components, adding 1 to the K value when the center frequency distances of all the IMFs are larger than 0.1, and repeating the steps until the center distance of any adjacent IMF is smaller than 0.1, and obtaining the optimal K value at the moment.
5. The method for extracting the vibration signal characteristics of the gear decomposed by the variation mode according to claim 1, wherein the IMF components decomposed in the step 5 are subjected to signal reconstruction by using correlation and peak values as indexes, and the IMF components with correlation coefficients larger than the average value and kurtosis larger than 3 are screened.
6. The method for extracting features of vibration signals of a gear decomposed by a variation mode according to claim 1, wherein in the step 6, a dung beetle optimization algorithm is used to find the best combination of L and M in the range that the filter length L belongs to [100,500] and the shift number M belongs to [1,7], so that the effect of enhancing weak pulses in the vibration signals by MCKD is the best.
7. The method for extracting vibration signal characteristics of a gear using a variation modal decomposition according to claim 1, wherein in step 7, characteristics of the vibration signal are extracted from different angles using three different signal entropies of fuzzy entropy, approximate entropy, and sample entropy.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310612456.2A CN116610907B (en) | 2023-05-27 | 2023-05-27 | Gear vibration signal characteristic extraction method based on variational modal decomposition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310612456.2A CN116610907B (en) | 2023-05-27 | 2023-05-27 | Gear vibration signal characteristic extraction method based on variational modal decomposition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116610907A true CN116610907A (en) | 2023-08-18 |
CN116610907B CN116610907B (en) | 2024-03-19 |
Family
ID=87683254
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310612456.2A Active CN116610907B (en) | 2023-05-27 | 2023-05-27 | Gear vibration signal characteristic extraction method based on variational modal decomposition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116610907B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117235643A (en) * | 2023-09-06 | 2023-12-15 | 南京林业大学 | Early weak fault diagnosis method for rolling bearing |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3511862A1 (en) * | 2018-01-12 | 2019-07-17 | Qognify Ltd. | System and method for dynamically ordering video channels according to rank of abnormal detection |
CN111413089A (en) * | 2020-04-08 | 2020-07-14 | 北华大学 | Gear fault diagnosis method based on combination of VMD entropy method and VPMCD |
WO2021056727A1 (en) * | 2019-09-27 | 2021-04-01 | 山东科技大学 | Joint noise reduction method based on variational mode decomposition and permutation entropy |
CN113408458A (en) * | 2021-06-29 | 2021-09-17 | 北京交通大学 | Circuit breaker fault diagnosis method based on vibration signal characteristic entropy fusion |
CN114077847A (en) * | 2021-10-27 | 2022-02-22 | 河北建投能源投资股份有限公司 | Variable working condition rolling bearing fault diagnosis method, system, equipment and medium |
CN115452381A (en) * | 2022-09-29 | 2022-12-09 | 湘潭大学 | Rolling bearing fault diagnosis method under strong noise interference based on multi-modal decomposition and long-short term memory network |
CN115730199A (en) * | 2022-11-10 | 2023-03-03 | 天地(常州)自动化股份有限公司北京分公司 | Method and system for noise reduction and fault feature extraction of vibration signal of rolling bearing |
EP4162865A1 (en) * | 2021-10-06 | 2023-04-12 | Origin Wireless, Inc. | System and method for radio-based sleep monitoring |
-
2023
- 2023-05-27 CN CN202310612456.2A patent/CN116610907B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3511862A1 (en) * | 2018-01-12 | 2019-07-17 | Qognify Ltd. | System and method for dynamically ordering video channels according to rank of abnormal detection |
WO2021056727A1 (en) * | 2019-09-27 | 2021-04-01 | 山东科技大学 | Joint noise reduction method based on variational mode decomposition and permutation entropy |
CN111413089A (en) * | 2020-04-08 | 2020-07-14 | 北华大学 | Gear fault diagnosis method based on combination of VMD entropy method and VPMCD |
CN113408458A (en) * | 2021-06-29 | 2021-09-17 | 北京交通大学 | Circuit breaker fault diagnosis method based on vibration signal characteristic entropy fusion |
EP4162865A1 (en) * | 2021-10-06 | 2023-04-12 | Origin Wireless, Inc. | System and method for radio-based sleep monitoring |
CN114077847A (en) * | 2021-10-27 | 2022-02-22 | 河北建投能源投资股份有限公司 | Variable working condition rolling bearing fault diagnosis method, system, equipment and medium |
CN115452381A (en) * | 2022-09-29 | 2022-12-09 | 湘潭大学 | Rolling bearing fault diagnosis method under strong noise interference based on multi-modal decomposition and long-short term memory network |
CN115730199A (en) * | 2022-11-10 | 2023-03-03 | 天地(常州)自动化股份有限公司北京分公司 | Method and system for noise reduction and fault feature extraction of vibration signal of rolling bearing |
Non-Patent Citations (3)
Title |
---|
QINGHUA WANG等: "Utilizing SVD and VMD for Denoising Non-Stationary Signals of Roller Bearings", SENSORS (BASEL), vol. 22, no. 1, pages 1 - 17 * |
任良等: "基于SSA-VMD-MCKD的强背景噪声环境下滚动轴承故障诊断", 振动与冲击, vol. 42, no. 3, pages 220 - 222 * |
陈剑等: "基于SVD-VMD和SVM滚动轴承故障诊断方法", 电子测量与仪器学报, vol. 36, no. 1, pages 221 - 223 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117235643A (en) * | 2023-09-06 | 2023-12-15 | 南京林业大学 | Early weak fault diagnosis method for rolling bearing |
CN117235643B (en) * | 2023-09-06 | 2024-03-22 | 南京林业大学 | Early weak fault diagnosis method for rolling bearing |
Also Published As
Publication number | Publication date |
---|---|
CN116610907B (en) | 2024-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101303764B (en) | Method for self-adaption amalgamation of multi-sensor image based on non-lower sampling profile wave | |
CN116610907B (en) | Gear vibration signal characteristic extraction method based on variational modal decomposition | |
CN110361778B (en) | Seismic data reconstruction method based on generation countermeasure network | |
CN105913393A (en) | Self-adaptive wavelet threshold image de-noising algorithm and device | |
CN110865357A (en) | Laser radar echo signal noise reduction method based on parameter optimization VMD | |
KR20170056687A (en) | Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation | |
CN112364291B (en) | Empirical mode decomposition method and device for pre-filtering extreme point optimization set | |
CN110490219B (en) | Method for reconstructing seismic data by U-net network based on texture constraint | |
CN110568483A (en) | Automatic evaluation method for seismic linear noise suppression effect based on convolutional neural network | |
CN112686096A (en) | Rolling bearing fault diagnosis method based on multi-scale diffusion entropy and VPMCD | |
CN100417191C (en) | Method of reducing noise for combined images | |
CN117390380B (en) | Data analysis method in oil-residue separation system | |
CN112419203B (en) | Diffusion weighted image compressed sensing recovery method and device based on countermeasure network | |
CN113917490A (en) | Laser wind finding radar signal denoising method and device | |
CN116699526A (en) | Vehicle millimeter wave radar interference suppression method based on sparse and low-rank model | |
CN114779238A (en) | Life detection method and system based on radar signals | |
CN114386452B (en) | Nuclear power circulating water pump sun gear fault detection method | |
CN111242873A (en) | Image denoising method based on sparse representation | |
CN111598822A (en) | Image fusion method based on GFRW and ISCM | |
CN110673138A (en) | Ground penetrating radar image processing method based on singular value decomposition and fuzzy C mean value method | |
CN112380932B (en) | Vibration signal characteristic value selection method and elevator health state evaluation or fault diagnosis method | |
CN117473414A (en) | Bearing fault position identification method based on low-noise time-frequency image | |
CN111709116A (en) | Blind signal decomposition method based on similarity measurement | |
CN115563480A (en) | Gear fault identification method for screening octave geometric modal decomposition based on kurtosis ratio coefficient | |
CN112285793B (en) | Magnetotelluric denoising method and system |
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 |