CN108363109A - The harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data - Google Patents
The harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data Download PDFInfo
- Publication number
- CN108363109A CN108363109A CN201810029428.7A CN201810029428A CN108363109A CN 108363109 A CN108363109 A CN 108363109A CN 201810029428 A CN201810029428 A CN 201810029428A CN 108363109 A CN108363109 A CN 108363109A
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- atom
- nuclear magnetic
- dictionary
- resonance data
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/14—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electron or nuclear magnetic resonance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Abstract
The invention belongs to ground nuclear magnetic resonance data processing fields, are particularly a kind of harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data, and this method includes acquisition ground nuclear magnetic resonance data and typing；The industrial frequency harmonic dictionary for including be possible to fundamental frequency is established according to the data of typing；Calculate the autopower spectral density of each atom in ground nuclear magnetic resonance data and dictionary；Calculate the cross-power spectral density of each atom in ground nuclear magnetic resonance data and dictionary；The amplitude square related coefficient of each atom in ground nuclear magnetic resonance data and dictionary is calculated according to autopower spectral density and cross-power spectral density；It sums to the corresponding amplitude square related coefficient of each atom, the frequency of the corresponding atom of maximum value is industrial frequency harmonic fundamental frequency, the method of the present invention interferes the efficient of fundamental frequency search to industrial frequency harmonic in the nuclear magnetic resonance data of ground, and accuracy is good, substantially increases the processing speed of ground nuclear magnetic resonance data.
Description
Technical field
The invention belongs to ground nuclear magnetic resonance data processing fields, are particularly a kind of ground nuclear magnetic resonance data
Harmonic wave interference fundamental frequency method for fast searching.
Background technology
Ground nuclear magnetic resonance (surface nuclear magnetic resonance, SNMR) be it is a kind of it is non-intrusive,
The geophysical method of direct detection groundwater occurrence attribute, has many advantages, such as high resolution and quantitative interpretation.But the side SNMR
The signal that method obtains is very faint, is vulnerable to the interference of various ambient noises, low signal-to-noise ratio is still that SNMR often is faced with most
Big challenge.It is the strongest by the interference of the industrial frequency harmonic of the generations such as high-pressure delivery power grid and electrical equipment in SNMR data.For work
Frequency harmonic wave interference, presently the most effective method are industrial frequency harmonic modeling methods.
Industrial frequency harmonic modeling needs to utilize accurate fundamental frequency.In practice, the fundamental frequency f of power frequency0It cannot be guaranteed that being accurate
50Hz, but change within the scope of 50 ± 0.1Hz, therefore f0Search it is particularly significant in SNMR Data processings.On the one hand, f0
Accuracy directly affect the accuracy of data processed result, it is desirable that at least precision to 1mHz.If being all made of f0=50Hz is carried out
Modeling, harmonic error can increase with the increase of order, lead to modeling result mistake.On the other hand, since SNMR tests number
According to have a large amount of single measurements as a result, each measurement result fundamental frequency and differ, f0Search speed will have a direct impact at data
The speed of reason.
Patent CN106772646A discloses a kind of ground nuclear magnetic resonance method for extracting signal.Disappeared by Harmonic Modeling method
Except industrial frequency harmonic noise, the fundamental frequency of industrial frequency harmonic is searched for using the method for variable step scanning, but this method needs repeatedly calculating more
The amplitude and phase of harmonic wave, calculating speed are very slow under a fundamental frequency.
Patent CN105549097A discloses a kind of transient electromagnetic signal power frequency and its harmonic wave interference removing method.Utilize two
The fundamental frequency of point-score iterative search industrial frequency harmonic.This method is gradually reduced and is changed by the mean square deviation minimum point of searching random noise
For section, but the amplitude and phase of harmonic wave under multiple fundamental frequencies equally need to be repeatedly calculated, calculating speed is slow.
Invention content
Technical problem to be solved by the present invention lies in providing, a kind of harmonic wave interference fundamental frequency of ground nuclear magnetic resonance data is fast
Fast searching method.Under the premise of ensureing search precision, accelerate the search speed of fundamental frequency, moreover it is possible to simultaneously scan for two kinds of harmonic wave bases
Frequently, so improve magnetic resonance signal processing speed, be particularly suitable for including the ground magnetic resonance two and three dimensions of mass data
Measure situation.
The invention is realized in this way
A kind of harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data, this method include,
Acquire ground nuclear magnetic resonance data and typing；
The industrial frequency harmonic dictionary for including be possible to fundamental frequency is established according to the data of typing；
Calculate the autopower spectral density of each atom in ground nuclear magnetic resonance data and dictionary；
Calculate the cross-power spectral density of each atom in ground nuclear magnetic resonance data and dictionary；
It is calculated according to autopower spectral density and cross-power spectral density each former in ground nuclear magnetic resonance data and dictionary
The amplitude square related coefficient of son；
It sums to the corresponding amplitude square related coefficient of each atom, the frequency of the corresponding atom of maximum value is industrial frequency harmonic
Fundamental frequency.
Further, each of industrial frequency harmonic dictionary is classified as an atom di, generating function is：
Wherein, n is harmonic wave number, f0iFor industrial frequency harmonic fundamental frequency, the value range of fundamental frequency arrives for 49.9Hz in dictionary
50.1Hz, step-length are 1mHz or 0.1mHz, and t is the sampling time.
Further, amplitude square related coefficient is calculated as：
Wherein, PSS(f) it is each atom d in ground nuclear magnetic resonance dataiAutopower spectral density, PDD(f) it is in dictionary
Each atom diAutopower spectral density, PSD(f) it is each atom d in ground nuclear magnetic resonance data and dictionaryiAlternating power
Spectrum density, CSD(f) it is each atom diCorresponding amplitude square related coefficient.
Fundamental frequency for some European state-owned railroads power transmission lines is about the one third of power frequency, and the present invention may extend to power frequency
Two fundamental frequency method for fast searching when harmonic wave and railway harmonic exist simultaneously.
A kind of harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data,
Acquire ground nuclear magnetic resonance data and typing；
The industrial frequency harmonic dictionary for including be possible to fundamental frequency is established according to the data of typing；
Calculate the autopower spectral density of each atom in ground nuclear magnetic resonance data and dictionary；
Calculate the cross-power spectral density of each atom in ground nuclear magnetic resonance data and dictionary；
It is calculated according to autopower spectral density and cross-power spectral density each former in ground nuclear magnetic resonance data and dictionary
The amplitude square related coefficient of son；
It sums to the corresponding amplitude square related coefficient of each atom, two peak values is obtained, wherein the corresponding frequency of larger peak value
Rate is railway harmonic fundamental frequency, and frequency corresponding compared with small leak is the one third of industrial frequency harmonic fundamental frequency.
Further, each of industrial frequency harmonic dictionary is classified as an atom di, generating function is：
Wherein, n is harmonic wave number, f0iFor industrial frequency harmonic fundamental frequency, the value range of fundamental frequency arrives for 15.6Hz in dictionary
17.6Hz, step-length are 1mHz or 0.1mHz, and t is the sampling time.
Further, amplitude square related coefficient is calculated as：
Wherein, PSS(f) it is each atom d in ground nuclear magnetic resonance dataiAutopower spectral density, PDD(f) it is in dictionary
Each atom diAutopower spectral density, PSD(f) it is each atom d in ground nuclear magnetic resonance data and dictionaryiAlternating power
Spectrum density, CSD(f) it is each atom diCorresponding amplitude square related coefficient.
Compared with prior art, the present invention advantageous effect is：The harmonic wave interference fundamental frequency of ground nuclear magnetic resonance data is quick
Searching method, need not calculate the amplitude and phase of each fundamental frequency harmonics, and then avoid the multiple solution of large-scale equation group.Root
It is different according to ground nuclear magnetic resonance signal and the interatomic amplitude square correlation of different fundamental frequencies, and then realize to industrial frequency harmonic fundamental frequency
Search.This method substantially increases search speed under the premise of ensureing search precision, moreover it is possible to simultaneously scan for two kinds of harmonic wave bases
Frequently so that the two and three dimensions ground nuclear magnetic resonance method data comprising multigroup record are treated as possibility in real time.
Through experiment, the method for the present invention interferes the efficient of fundamental frequency search to industrial frequency harmonic in the nuclear magnetic resonance data of ground, accurate
True property is good, substantially increases the processing speed of ground nuclear magnetic resonance data.
Description of the drawings
Fig. 1 is the flow chart of the harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data；
The amplitude square related coefficient of atom and nuclear magnetic resonance data, (a) fundamental frequency are equal in Fig. 2 dictionaries
(b) fundamental frequency difference 10mHz (c) fundamental frequency differs 20mHz；
The sum of Fig. 3 nuclear magnetic resonance datas amplitude square related coefficient corresponding with atom；
Fig. 4 is there are when two harmonic fundamental frequencies, the sum of nuclear magnetic resonance data amplitude square related coefficient corresponding with atom.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to embodiments, to the present invention
It is further elaborated.It should be appreciated that described herein, specific examples are only used to explain the present invention, is not used to limit
The fixed present invention.
Embodiment 1：
Be 100nV, relaxation time 0.2s, Larmor frequency 2310Hz with one group of initial amplitude, initial phase be π/
The ground nuclear magnetic resonance signal of 3rad；Industrial frequency harmonic fundamental frequency is 50.0125Hz, and the amplitude and phase of 100 harmonic waves are in random point
Cloth；Random noise is that the data instance of 20nV illustrates.
Referring to Fig. 1, step 1, typing ground nuclear magnetic resonance data establish the industrial frequency harmonic dictionary for including be possible to fundamental frequency
The each of D, dictionary D is classified as an atom di, generating function is：
Wherein, n is harmonic wave number, and 100, f are taken in this implementation0iFor industrial frequency harmonic fundamental frequency, the value range of fundamental frequency in dictionary D
For 49.9Hz to 50.1Hz, step-length 0.1mHz, t are the sampling time.
Step 2, each atom d in ground nuclear magnetic resonance data and dictionary is calculated separatelyiAutopower spectral density PSS(f) and
PDD(f), each atom d in ground nuclear magnetic resonance data and dictionaryiCross-power spectral density PSD(f).Autopower spectral density PSS
(f) and PDD(f) and cross-power spectral density PSD(f) calculation formula for being all made of power spectral density below obtains：
The calculation formula of power spectral density is
Rxy(m)=E [xk+my* k]
Wherein, xkAnd ykIt is two groups of time series datas respectively, k and m are moment, RxyFor cross-correlation function, E [] is represented
Mathematic expectaion, y*For the conjugate complex number of y sequences.When x and y respectively represent the atom D in ground nuclear magnetic resonance data S and dictionary
When, PSDFor crosspower spectrum spectrum density.When x and y are the atom D in ground nuclear magnetic resonance data S or dictionary, PSSAnd PDDFor
Auto-power spectrum spectrum density.
Amplitude square related coefficient is calculated according to formula (2)：
Fig. 2 (a) is the width when fundamental frequency of atom is equal with the fundamental frequency that industrial frequency harmonic in nuclear magnetic resonance data interferes in dictionary
Spend squared correlation coefficient；Fig. 2 (b) is that the fundamental frequency of atom is differed with the fundamental frequency that industrial frequency harmonic in nuclear magnetic resonance data interferes in dictionary
Amplitude square related coefficient when 10mHz；Fig. 2 (c) is that industrial frequency harmonic is dry in the fundamental frequency and nuclear magnetic resonance data of atom in dictionary
Amplitude square related coefficient when the fundamental frequency difference 20mHz disturbed.
Step 3, each atom d obtained according to step 2iCorresponding amplitude square related coefficient CSD(f), respectively to CSD
(f) summation that adds up obtains C∑(di), as a result see Fig. 3, wherein C∑(di) maximum value correspondence atom diFrequency f0iAs power frequency is humorous
The fundamental frequency f of wave0=50.0125Hz.
Embodiment 2：
Be 100nV, relaxation time 0.2s, Larmor frequency 2310Hz with one group of initial amplitude, initial phase be π/
The ground nuclear magnetic resonance signal of 3rad；Industrial frequency harmonic fundamental frequency is 50.022Hz, and the amplitude and phase of 100 harmonic waves are in random point
Cloth；Railway harmonic fundamental frequency is 16.602Hz, and the amplitude and phase of 300 harmonic waves are in random distribution；Random noise is the number of 20nV
It is illustrated for.
Step 1, typing ground nuclear magnetic resonance data establish industrial frequency harmonic the dictionary D, dictionary D for including be possible to fundamental frequency
Each be classified as an atom di, generating function is：
Based in ground nuclear magnetic resonance data, the prior information of industrial frequency harmonic and railway line harmonic wave interference, fundamental frequency in dictionary
Value range be 15.6Hz to 17.6Hz, the present embodiment takes step-length 1mHz, and harmonic order n is the sampling time for 300, t.
Step 2, each atom d in ground nuclear magnetic resonance data and dictionary is calculated separatelyiAuto-power spectrum spectrum density PSS(f)
And PDD(f), each atom d in ground nuclear magnetic resonance data and dictionaryiCross-power spectral density PSD(f).It is calculated according to formula (2)
Obtain amplitude square related coefficient：
Step 3, each atom d obtained according to step 2iCorresponding amplitude square related coefficient CSD(f), respectively to CSD
(f) summation that adds up obtains C∑(di), Fig. 4 is as a result seen, wherein larger peak value is the fundamental frequency 16.602Hz of railway harmonic, compared with small peak
It is 16.674 × 3=50.022Hz of industrial frequency harmonic fundamental frequency that value, which is multiplied by 3,.
Claims (6)
1. a kind of harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data, which is characterized in that this method includes,
Acquire ground nuclear magnetic resonance data and typing；
The industrial frequency harmonic dictionary for including be possible to fundamental frequency is established according to the data of typing；
Calculate the autopower spectral density of each atom in ground nuclear magnetic resonance data and dictionary；
Calculate the cross-power spectral density of each atom in ground nuclear magnetic resonance data and dictionary；
Each atom in ground nuclear magnetic resonance data and dictionary is calculated according to autopower spectral density and cross-power spectral density
Amplitude square related coefficient；
It sums to the corresponding amplitude square related coefficient of each atom, the frequency of the corresponding atom of maximum value is industrial frequency harmonic base
Frequently.
2. according to the method for claim 1, which is characterized in that each of industrial frequency harmonic dictionary is classified as an atom di, life
It is at function：
Wherein, n is harmonic wave number, f0iFor industrial frequency harmonic fundamental frequency, the value range of fundamental frequency is 49.9Hz to 50.1Hz in dictionary, is walked
A length of 1mHz or 0.1mHz, t are the sampling time.
3. according to the method for claim 1, which is characterized in that amplitude square related coefficient is calculated as：
Wherein, PSS(f) it is each atom d in ground nuclear magnetic resonance dataiAutopower spectral density, PDD(f) it is each in dictionary
Atom diAutopower spectral density, PSD(f) it is each atom d in ground nuclear magnetic resonance data and dictionaryiCross power spectrum it is close
Degree, CSD(f) it is each atom diCorresponding amplitude square related coefficient.
4. a kind of harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data, which is characterized in that
Acquire ground nuclear magnetic resonance data and typing；
The harmonic wave dictionary for including be possible to fundamental frequency is established according to the data of typing；
Calculate the autopower spectral density of each atom in ground nuclear magnetic resonance data and dictionary；
Calculate the cross-power spectral density of each atom in ground nuclear magnetic resonance data and dictionary；
Each atom in ground nuclear magnetic resonance data and dictionary is calculated according to autopower spectral density and cross-power spectral density
Amplitude square related coefficient；
It sums to the corresponding amplitude square related coefficient of each atom, two peak values is obtained, wherein the corresponding frequency of larger peak value is
Railway harmonic fundamental frequency, frequency corresponding compared with small leak are the one third of industrial frequency harmonic fundamental frequency.
5. according to the method for claim 4, which is characterized in that each of harmonic wave dictionary is classified as an atom di, generate letter
Number is：
Wherein, n is harmonic wave number, f0iFor harmonic fundamental frequency, the value range of fundamental frequency is 15.6Hz to 17.6Hz in dictionary, and step-length is
1mHz or 0.1mHz, t are the sampling time.
6. according to the method for claim 4, which is characterized in that amplitude square related coefficient is calculated as：
Wherein, PSS(f) it is each atom d in ground nuclear magnetic resonance dataiAutopower spectral density, PDD(f) it is each in dictionary
Atom diAutopower spectral density, PSD(f) it is each atom d in ground nuclear magnetic resonance data and dictionaryiCross power spectrum it is close
Degree, CSD(f) it is each atom diCorresponding amplitude square related coefficient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810029428.7A CN108363109B (en) | 2018-01-12 | 2018-01-12 | The harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810029428.7A CN108363109B (en) | 2018-01-12 | 2018-01-12 | The harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108363109A true CN108363109A (en) | 2018-08-03 |
| CN108363109B CN108363109B (en) | 2019-06-21 |
Family
ID=63011337
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810029428.7A Active CN108363109B (en) | 2018-01-12 | 2018-01-12 | The harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108363109B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104898172A (en) * | 2015-05-19 | 2015-09-09 | 吉林大学 | Cross-correlation-based nuclear magnetic resonance full wave signal noise filtering method |
| CN105551501A (en) * | 2016-01-22 | 2016-05-04 | 大连民族大学 | Harmonic signal fundamental frequency estimation algorithm and device |
| CN105785124A (en) * | 2016-03-07 | 2016-07-20 | 国网技术学院 | Method for measuring harmonics and interharmonics of electric power system through spectrum estimation and cross correlation |
| US20160211823A1 (en) * | 2015-01-21 | 2016-07-21 | Qualcomm Incorporated | Integrated filters in output match elements |
| CN106772646A (en) * | 2016-12-19 | 2017-05-31 | 吉林大学 | A kind of ground nuclear magnetic resonance method for extracting signal |
-
2018
- 2018-01-12 CN CN201810029428.7A patent/CN108363109B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160211823A1 (en) * | 2015-01-21 | 2016-07-21 | Qualcomm Incorporated | Integrated filters in output match elements |
| CN104898172A (en) * | 2015-05-19 | 2015-09-09 | 吉林大学 | Cross-correlation-based nuclear magnetic resonance full wave signal noise filtering method |
| CN105551501A (en) * | 2016-01-22 | 2016-05-04 | 大连民族大学 | Harmonic signal fundamental frequency estimation algorithm and device |
| CN105785124A (en) * | 2016-03-07 | 2016-07-20 | 国网技术学院 | Method for measuring harmonics and interharmonics of electric power system through spectrum estimation and cross correlation |
| CN106772646A (en) * | 2016-12-19 | 2017-05-31 | 吉林大学 | A kind of ground nuclear magnetic resonance method for extracting signal |
Non-Patent Citations (1)
| Title |
|---|
| 孙晓东 等: ""复杂噪声背景下谐波信号频率估计新方法"", 《吉林大学学报（工学版）》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108363109B (en) | 2019-06-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8067942B2 (en) | Method for locating phase to ground faults in DC distribution systems | |
| Zhu | Exact harmonics/interharmonics calculation using adaptive window width | |
| Li et al. | Damage localization of stacker’s track based on EEMD-EMD and DBSCAN cluster algorithms | |
| Han et al. | The fault feature extraction of rolling bearing based on EMD and difference spectrum of singular value | |
| Liu et al. | Correlation analysis for spread-spectrum induced-polarization signal processing in electromagnetically noisy environments | |
| CN108345039B (en) | A method of eliminating adjacent frequency harmonic wave interference in ground nuclear magnetic resonance data | |
| Xie et al. | De-noising methods for NMR logging echo signals based on wavelet transform | |
| Chanda et al. | Application of wavelet multiresolution analysis for classification of faults on transmission lines | |
| Shi et al. | Signal extraction using complementary ensemble empirical mode in pipeline magnetic flux leakage nondestructive evaluation | |
| Zhang et al. | A precise and adaptive algorithm for interharmonics measurement based on iterative DFT | |
| Jin et al. | A novel power harmonic analysis method based on Nuttall-Kaiser combination window double spectrum interpolated FFT algorithm | |
| Yang et al. | A noise evaluation method for CSEM in the frequency domain based on wavelet transform and analytic envelope | |
| Zhang et al. | Harmonic detection method based on permutation entropy and variational modal decomposition optimized by genetic algorithm | |
| Chen et al. | A bilateral second-order synchrosqueezing transform and application to vibration monitoring of aerospace engine | |
| Li et al. | Research on power quality disturbance identification and classification technology in high noise background | |
| CN108363109B (en) | The harmonic wave interference fundamental frequency method for fast searching of ground nuclear magnetic resonance data | |
| CN112083509B (en) | Method for detecting induced polarization abnormity in time-frequency electromagnetic method | |
| CN110245321A (en) | Duration power quality disturbances method and device based on match tracing sparse decomposition | |
| CN109541304A (en) | The weak amplitude harmonic detecting method of power grid high order based on six minimum secondary lobe window interpolation | |
| Zhang et al. | Partial discharge signal self-adaptive sparse decomposition noise abatement based on spectral kurtosis and S-transform | |
| Kolobov et al. | The KVVN-7 multifunction digital measuring station for electromagnetic monitoring of seismoactive zones | |
| Zhao et al. | Analysis of transient waveform based on combined short time Fourier transform and wavelet transform | |
| LIN et al. | Segmented time-frequency peak filtering for random noise reduction of MRS oscillating signal | |
| Ortega et al. | Hilbert–Huang transform analysis of storm waves | |
| Liguo et al. | Power systems faults location with traveling wave based on Hilbert-Huang transform |
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 |