CN103823244A - Magnetic resonance three-component noise removing device and noise removing method - Google Patents

Magnetic resonance three-component noise removing device and noise removing method Download PDF

Info

Publication number
CN103823244A
CN103823244A CN201410088813.0A CN201410088813A CN103823244A CN 103823244 A CN103823244 A CN 103823244A CN 201410088813 A CN201410088813 A CN 201410088813A CN 103823244 A CN103823244 A CN 103823244A
Authority
CN
China
Prior art keywords
signal
magnetic resonance
coil
component
noise
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
Application number
CN201410088813.0A
Other languages
Chinese (zh)
Other versions
CN103823244B (en
Inventor
林婷婷
张扬
林君
万玲
张思远
杜官峰
蒋川东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN201410088813.0A priority Critical patent/CN103823244B/en
Publication of CN103823244A publication Critical patent/CN103823244A/en
Application granted granted Critical
Publication of CN103823244B publication Critical patent/CN103823244B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

The invention relates to a magnetic resonance three-component noise removing device. The magnetic resonance three-component noise removing device is formed through a structure. The structure is characterized in that an upper computer is connected with a transmitting coil through a control system, a transmitting circuit and a resonance matching capacitor; the transmitting circuit is connected with a power supply; the control system is connected with a receiving coil and a reference coil respectively through an A/D (Analog to Digital) collection and amplification circuit and a conditioning circuit. The magnetic resonance three-component noise removing device is applicable to a detection area with complex and changeable environment noise and uneven space distribution of the noise; due to the adoption of relevance between x-component signals, y-component signals and z-component signals which are sensed by coils paved in the same detection place, the z-component noise is removed, then filtered and processed to obtain reliable magnetic resonance signals, so that the wiring space is greatly saved and the processing result is flexible and accurate; in particular, when sharp-peak noise with great amplitude exists, the noise removing effect is excellent. Through the magnetic resonance three-component noise removing device, the influence of strong electromagnetic interference near the cities, the villages and the water-deficient areas is eliminated, and the problem of limited application of a nuclear magnetic resonance underground water detection method is solved.

Description

Magnetic resonance three-component noise cancellation apparatus and noise-eliminating method
Technical field:
The present invention relates to a kind of geophysical prospecting method and data interpretation thereof, in particular for gathering nuclear magnetic resonance (Nuclear Magnetic Resonance, NMR) three-component underground water detection method and utilizing magnetic resonance three-component to eliminate the data processing method of noise.
Background technology:
Nuclear magnetic resonance underground water detection method has as a kind of method of direct detection underground water that resolution is high, efficiency is high, quantity of information is abundant and the advantage such as uniqueness of solution.But due to the signal that instrument receives very faint (receive and lie prostrate level), highly sensitive receiving trap is very easily subject to the interference of ambient noise composition.Chinese population is intensive, electrical network prosperity, line of electric force and power equipment are covered with cities and towns and village, electromagnetic interference (EMI) situation is serious, and water-deficient area major part near cities and village etc. these strong jamming areas, this is restricted the application of nuclear magnetic resonance underground water detection method, sometimes even can not obtain reliable NMR signal.Noise to received signal carries out effective filtering and becomes accurate inversion interpretation hydrogeological parameter result important prerequisite to realize the reliable extraction of NMR signal.
That CN101251606A discloses is a kind of " in Testing of Feeble Signals instrument in useful signal frequency band industrial frequency harmonic codan ", is the compensating circuit that utilizes frequency, phase place, amplitude comprehensively to follow the tracks of industrial frequency harmonic in useful signal frequency band to cause interference.The connecting and disconnecting of master microprocessor control signal sensor and impedance matching network, impedance matching network is connected with programmable amplifier through prime amplifier, broadband filter, output signal envelope detecting device by control line and data bus be connected from microprocessor, programmable amplifier connects DDS signal generator through phaselocked loop.This circuit effectively suppresses the industrial frequency harmonic in useful signal frequency band in detecting instrument to be disturbed, follow the tracks of in real time frequency, phase place and the amplitude of this interference, by its filtering from receive signal, and do not affect other useful signal, do not affect the useful signal frequency content overlapping with this interfering frequency.
CN102053280A discloses a kind of " with nuclear magnetic resonance underground water detection system and the detection method of reference coil ", by the Full wave shape data of noise signal in the NMR signal in multi-channel A/D collecting unit synchronous acquisition transmitting/receiving coil and reference coil, the noise signal gathering by computing reference coil with and the maximum correlation of magnetic resonance signal, realize the laying of reference coil optimum position and quantity, the in the situation that of signal and noise statistics the unknown, adopt variable step size adaptive algorithm, offset to greatest extent transmitting/receiving coil and obtain the noise in NMR signal, realize the extraction that the complicated very noisy of many field sources disturbs lower and magnetic resonance signal, improve equipment anti-interference ability, for searching ground water provides a kind of reliable sniffer and method.
The industrial frequency harmonic codan of foregoing invention can only solve the filtering of the some industrial frequency harmonics in useful signal frequency band, helpless for electromagnetic interference (EMI) such as other industrial frequency harmonics and unusual noises; Require noise signal and magnetic resonance signal in reference coil to have good correlativity with the underground water detection method of reference coil, thereby distance between reference coil and main coil is at least greater than two times of wire diameters, increase wiring space, and due to the unevenness of complicacy, polytrope and the space distribution of neighbourhood noise, sometimes can not obtain desirable denoising effect.
For the noise problem in magnetic resonance signal, by laying special shape coil, adopt methods such as increasing method, block opposition method, sinusoidal opposition method and the notch filter of signal stacking fold to be also used to the filtering of electromagnetic noise, these methods have obtained good effect under certain precondition, but when noise intensity is large, interference component is complicated, signal center frequency and power frequency disturb frequency difference hour, are difficult to obtain good filter effect.
Summary of the invention:
Object of the present invention is just for above-mentioned the deficiencies in the prior art, provide one to be applicable to neighbourhood noise complexity, the magnetic resonance three-component noise-eliminating method of test operation limited space, thereby the noise to received signal carries out effective filtering to realize the reliable extraction of NMR signal, and accurate inversion interpretation hydrogeological parameter is had to vital role and beyond thought effect.
The object of the invention is to be achieved through the following technical solutions:
Magnetic resonance three-component noise cancellation apparatus, through control system 2, radiating circuit 3 with join humorous electric capacity 5 and be connected with transmitting coil 6 by host computer 1, radiating circuit 3 is connected with power supply 4, and control system 2 connects and composes respectively at receiving coil 7, reference coil 8 and reference coil 9 through A/D collection 12, amplifying circuit 11 and modulate circuit 10.
Centered by point of fixity O along the orthogonal laying reference coil 8 of X, Y, Z axis, reference coil 9 and receiving coil 7.
Host computer 1 arranges emission parameter to control system 2, comprises and excites duration and stimulating frequency.The sine pulse equating with local Larmor frequency to transmission frequency in transmitting coil 6 by radiating circuit 3.NMR signal is spatially made up of three components of x, y, z, receives nuclear magnetic resonance z component signal with receiving coil 7, receives nuclear magnetic resonance x component with reference coil 8, receives nuclear magnetic resonance y component signal with reference coil 9.Wherein, receiving coil 7, reference coil 8 and reference coil 9 are three coils of orthogonal laying centered by a certain point of fixity O.The signal receiving passes through after signal conditioning circuit 10 and amplifying circuit 11, gathered by A/D acquisition module 12, finally pass to host computer 1, save data also starts to process this three components, utilize magnetic resonance three-component to eliminate noise with reference to noise-eliminating method, obtain reliable NMR signal.
Magnetic resonance three-component, with reference to noise-eliminating method, comprises the following steps:
In a, Ce district, select to survey place, lay transmitting coil 6;
B, take vertical direction as normal vector lay receiving coil 7, take North and South direction as normal vector lay reference coil 8, take east-west direction as normal vector lay reference coil 9.Receiving coil 7, reference coil 8 and reference coil 9 be orthogonal laying on locus centered by O point, and transmitting coil 6 must be arranged and be laid on same place with receiving coil 7, reference coil 8 and reference coil 9;
C, first measure noise one time, open magnetic resonance three component seismic data acquisition system, host computer 1 sends control word to control system 2, comprises the information that excites duration and stimulating frequency.Control system 2 produces control signal radiating circuit 3 is controlled, and radiating circuit 3 does not produce transmitter current to transmitting coil 6, and this is emitted as pseudo-transmitting.
D, receiving coil 7, reference coil 8 and reference coil 9 carry out after signal condition and three-channel amplifying circuit 11 amplify being gathered by A/D acquisition module 12 through three-channel modulate circuit 10 by receiving signal after puppet transmitting;
E, acquiescence North and South direction are x direction of principal axis, and the signal that reference coil 8 receives is x component noise; East-west direction is y direction of principal axis, and the signal that reference coil 9 receives is y component noise; Vertical direction is z direction of principal axis, and the signal that receiving coil 7 receives is z component noise;
F, measurement magnetic resonance signal, now host computer 1 sends control word to control system 2 again, comprises the information that excites duration and stimulating frequency.Control system 2 produces control signal radiating circuit 3 is controlled, and radiating circuit 3 produces by joining humorous electric capacity 5 exchange current that frequency is local Larmor frequency in transmitting coil 6;
G, receiving coil 7, reference coil 8 and reference coil 9 are gathered by A/D acquisition module 12 after sensing signal after three-channel modulate circuit 10 carries out signal condition and 11 amplifications of three-channel amplifying circuit, and the signal now receiving is respectively magnetic resonance x component signal, magnetic resonance y component signal and magnetic resonance z component signal;
H, with the handling object of the clock signal of three components recording simultaneously, time-domain filtering algorithm is proposed:
h z(t)=[h x(t)*X(t)]+[h y(t)*Y[t)]
Wherein h x, h yand h zfor x, y and z component noise, X and Y are x and y component signal, * is convolution algorithm;
I, obtain z component noise according to above-mentioned time-domain filtering algorithm, from z component signal, remove the reliable magnetic resonance signal that is that z component noise obtains;
J, change transmitter current size, obtain different pulses apart from different depth of stratums is surveyed, and repeats c-i process, obtains second group of magnetic resonance signal after processing;
K, can repeat above-mentioned j process, repeatedly survey, contrast repeatedly processing the magnetic resonance signal obtaining, obtain accurate information.
Beneficial effect: magnetic resonance three-component noise-eliminating method disclosed by the invention is applicable to the complicated and changeable and neighbourhood noise of neighbourhood noise inhomogeneous search coverage in space distribution, utilize the coil-induced x of same detection place laying and the correlativity of y component signal and z component signal, eliminate after z component noise, filtering processing obtains reliable magnetic resonance signal, wiring space but also result are not only significantly saved accurate flexibly, especially in the time that the larger spike noise of amplitude exists, de-noising effect is fine.The present invention has got rid of the impact of the suffered strong electromagnetic of near cities and village and water-deficient area, has solved nuclear magnetic resonance underground water detection method and has applied a limited difficult problem.
Accompanying drawing explanation:
Fig. 1 is magnetic resonance three component seismic data acquisition system structured flowchart
Fig. 2 is that coil is laid spatial relation figure
Fig. 3 is magnetic resonance three-component direction schematic diagram
1 host computer, 2 control system, 3 radiating circuits, 4 power supplys, 5 join humorous electric capacity, 6 transmitting coils, 7 receiving coils, 8 reference coils, 9 reference coils, 10 modulate circuits, 11 amplifying circuits, 12A/D acquisition module.
Embodiment:
Be described in further detail below in conjunction with drawings and Examples:
As shown in Figure 1, the exchange current that this system is local Larmor frequency by transmitting coil 6 to underground transmission frequency, produces alternating magnetic field.Under the exciting of the alternating magnetic field forming in ground, make hydrogen proton in underground water produce energy level transition, a large amount of hydrogen protons transit on high level.When removing supply current, these high level hydrogen protons are just got back to low-lying level state gradually, discharge a large amount of quantum of energys with Larmor frequency, in receiving coil 7, reference coil 8 and reference coil 9, induce MRS signal.
As shown in Figure 2, the particular orientation of laying due to receiving coil 7, reference coil 8 and reference coil 9, can be by the x, y, z component signal that receives MRS signal of selectivity, utilize the correlativity of x and y component signal and z component signal, time-domain filtering algorithm in use magnetic resonance three-component noise-eliminating method is with reference to eliminating after z component noise, and filtering processing can obtain reliable magnetic resonance signal.
Time-domain filtering algorithm in magnetic resonance three-component noise-eliminating method is calculated:
Set up h x(t)=[h x(t) * X (t)]+discrete data model of [hy (t) * Y (t)]:
h x 0 h y 0 h x - 1 h y - 1 h x - 2 h y - 2 . . . h x - n + 1 h y - n + 1 h x 1 h y 1 h x 0 h y 0 h x - 1 h y - 1 . . . h x - n + 2 h y - n + 2 h x 2 h y 2 h x 1 h y 1 h x 0 h y 0 . . . h x - n + 2 h y - n + 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . h x - 1 h y - 1 . . . . . . . . . h x 0 h y 0 . . . . . . . . . . . h x m - n h y m - n . . . . . . . . . . . h x m - n h y m - n . . . . . . . . . . . h x m - n h y m - n . . . . . . . . . . . . . . . . h x m - 1 h y m - 1 h x m - 2 h y m - 2 . . . . . h x m - n h y m - n X 0 Y 0 X 1 Y 1 . X n - 1 Y n - 1 = h z 0 h z 1 h z 2 h z 3 . . . . . . . . h z m - 1
Wherein, h x, h y, h z, the subscript of X and Y refers to time samples, and n is the number of samples of impulse response, and m is output h znumber of samples.
Variable in this model has comprised in the past and sampling number in the future, realizes by Computer Simulation.
Calculate result and need to estimate the solution X of shock response equation by model 0, X 2x n-2and Y 0, Y 2y n-1, filter coefficient depends on stratal configuration rather than x and y component.Remainder error r icomputing formula:
r i = h z i - h z i p
Wherein, p is estimated value; The computing formula of average error σ between measured value and estimated value:
σ = ( 1 m - n Σ i = 1 m r i 2 ) 1 / 2
Wherein, m is data amount check, and n is filter width; The poor computing formula of z component noise criteria:
S = ( 1 m - 1 Σ i = 1 m { h z i - h ‾ z } 2 ) 1 / 2
Wherein, z component noise mean value.
Snr computation formula:
SNR = σ S = ( 1 m - n Σ i = 1 m ( h z i - h z i p ) 1 1 m - 1 Σ i = 1 m ( h z i - h ‾ z ) 2 ) 1 / 2
Calculate the signal quality of Analysis signal-to-noise ratio (SNR) after with reference to de-noising.
Magnetic resonance three-component noise-eliminating method, comprises following order and step:
A, as shown in Figure 1, connect instrument, host computer 1 is connected with control system 2, control system 2 is connected with radiating circuit 3, radiating circuit 3 is connected with transmitting coil 6 by joining humorous electric capacity 5, and radiating circuit 3 is connected with power supply 4, and control system 2 gathers 12 with A/D and is connected, A/D gather three passages of 12 by amplifying circuit 11 and modulate circuit 10 respectively with receiving coil 7, reference coil 8 is connected with reference coil 9;
In b, Ce district, select to survey place, lay transmitting coil 6;
C, as shown in Figure 2, lays receiving coil 7 take vertical direction as normal vector, lay reference coil 8 take North and South direction as normal vector, lays reference coil 9 take east-west direction as normal vector.Receiving coil 7, reference coil 8 and reference coil 9 be orthogonal laying on locus centered by O point, and transmitting coil 6 must be arranged and be laid on same place with receiving coil 7, reference coil 8 and reference coil 9;
D, first measure noise one time, open magnetic resonance three component seismic data acquisition system, host computer 1 sends control word to control system 2, comprises the information that excites duration and stimulating frequency.Control system 2 produces control signal radiating circuit 3 is controlled, and radiating circuit 3 does not produce transmitter current to transmitting coil 6, and this is emitted as pseudo-transmitting.
E, receiving coil 7, reference coil 8 and reference coil 9 are undertaken after signal condition and three-channel amplifying circuit 11 amplify being gathered by A/D acquisition module 12 by receiving after puppet transmitting after signal through three-channel modulate circuit 10;
F, as shown in Figure 3, acquiescence North and South direction is x direction of principal axis, the signal that reference coil 8 receives is x component noise; East-west direction is y direction of principal axis, and the signal that reference coil 9 receives is y component noise; Vertical direction is z direction of principal axis, and the signal that receiving coil 7 receives is z component noise;
G, measurement magnetic resonance signal, now host computer 1 sends control word to control system 2 again, comprises the information that excites duration and stimulating frequency.Control system 2 produces control signal radiating circuit 3 is controlled, and radiating circuit 3 produces by joining humorous electric capacity 5 exchange current that frequency is local Larmor frequency in transmitting coil 6;
H, receiving coil 7, reference coil 8 and reference coil 9 are gathered by A/D acquisition module 12 after sensing signal after three-channel modulate circuit 10 carries out signal condition and 11 amplifications of three-channel amplifying circuit, and the signal now receiving is respectively magnetic resonance x component signal, magnetic resonance y component signal and magnetic resonance z component signal;
I, with the handling object of the clock signal of three components recording simultaneously, time-domain filtering algorithm is proposed:
h z(t)=[h x(t)*X(t)]+[h y(t)*Y(t)]
Wherein h x, h yand h zfor x, y and z component noise, X and Y are x and y component signal, * is convolution algorithm;
J, obtain z component noise according to above-mentioned time-domain filtering algorithm, from z component signal, remove the reliable magnetic resonance signal that is that z component noise obtains;
K, change transmitter current size, obtain different pulses apart from different depth of stratums is surveyed, and repeats d-j process, obtains second group of magnetic resonance signal after processing;
L, can repeat above-mentioned k process, repeatedly survey, contrast repeatedly processing the magnetic resonance signal obtaining, obtain accurate information.
Embodiment 1
Magnetic resonance three-component noise-eliminating method is a kind of noise-eliminating method of measuring magnetic resonance three component signal and utilizing time-domain filtering algorithm elimination noise, and the exchange current that this system is local Larmor frequency by transmitting coil 6 to underground transmission frequency, produces alternating magnetic field.Under the exciting of the alternating magnetic field forming in ground, make hydrogen proton in underground water produce energy level transition, a large amount of hydrogen protons transit on high level.When removing supply current, these high level hydrogen protons are just got back to low-lying level state gradually, discharge a large amount of quantum of energys with Larmor frequency, in receiving coil 7, reference coil 8 and reference coil 9, induce MRS signal.Acquiescence North and South direction is x direction of principal axis, east-west direction is y direction of principal axis, vertical direction is z direction of principal axis, the particular orientation of laying due to receiving coil 7, reference coil 8 and reference coil 9, can be by the x, y, z component signal that receives MRS signal of selectivity, utilize the correlativity of x and y component signal and z component signal, the time-domain filtering algorithm in use magnetic resonance three-component noise-eliminating method is with reference to eliminating after z component noise, and filtering processing can obtain reliable magnetic resonance signal.
Host computer 1 is connected with control system 2, control system 2 is connected with radiating circuit 3, radiating circuit 3 is connected with transmitting coil 6 by joining humorous electric capacity 5, radiating circuit 3 is connected with power supply 4, control system 2 gathers 12 with A/D and is connected, A/D gathers 12 and is connected with amplifying circuit 11, and amplifying circuit 11 is connected with modulate circuit 10, and modulate circuit 10 is connected with receiving coil 7, reference coil 8 and reference coil 9.
Magnetic resonance three-component noise-eliminating method, comprises following order and step:
In a, Ce district, select to survey place, lay transmitting coil 6;
B, take vertical direction as normal vector lay receiving coil 7, take North and South direction as normal vector lay reference coil 8, take east-west direction as normal vector lay reference coil 9.Receiving coil 7, reference coil 8 and reference coil 9 be orthogonal laying on locus centered by O point, and transmitting coil 6 must be arranged and be laid on same place with receiving coil 7, reference coil 8 and reference coil 9;
C, first measure noise one time, open magnetic resonance three component seismic data acquisition system, host computer 1 sends control word to control system 2, comprises the information that excites duration and stimulating frequency.Control system 2 produces control signal radiating circuit 3 is controlled, and radiating circuit 3 does not produce transmitter current to transmitting coil 6, and this is emitted as pseudo-transmitting.
D, receiving coil 7, reference coil 8 and reference coil 9 are undertaken after signal condition and three-channel amplifying circuit 11 amplify being gathered by A/D acquisition module 12 by receiving after puppet transmitting after signal through three-channel modulate circuit 10;
E, acquiescence North and South direction are x direction of principal axis, and the signal that reference coil 8 receives is x component noise; East-west direction is y direction of principal axis, and the signal that reference coil 9 receives is y component noise; Vertical direction is z direction of principal axis, and the signal that receiving coil 7 receives is z component noise;
F, measurement magnetic resonance signal, now host computer 1 sends control word to control system 2 again, comprises the information that excites duration and stimulating frequency.Control system 2 produces control signal radiating circuit 3 is controlled, and radiating circuit 3 produces by joining humorous electric capacity 5 exchange current that frequency is local Larmor frequency in transmitting coil 6;
G, receiving coil 7, reference coil 8 and reference coil 9 are gathered by A/D acquisition module 12 after sensing signal after three-channel modulate circuit 10 carries out signal condition and 11 amplifications of three-channel amplifying circuit, and the signal now receiving is respectively magnetic resonance x component signal, magnetic resonance y component signal and magnetic resonance z component signal;
H, with the handling object of the clock signal of three components recording simultaneously, time-domain filtering algorithm is proposed:
h z(t)=[h x(t)*X(t)]+[h y(t)*Y(t)]
Wherein h x, h yand h zfor x, y and z component noise, X and Y are x and y component signal, * is convolution algorithm;
I, obtain z component noise according to above-mentioned time-domain filtering algorithm, from z component signal, remove the reliable magnetic resonance signal that is that z component noise obtains;
J, change transmitter current size, obtain different pulses apart from different depth of stratums is surveyed, and repeats c-i process, obtains second group of magnetic resonance signal after processing;
K, repeat above-mentioned j process, carry out 16 times and survey, 16 groups of detection datas are processed.
Embodiment 2
Magnetic resonance three-component noise-eliminating method is a kind of reference noise-eliminating method of measuring magnetic resonance three component signal and utilizing time-domain filtering algorithm elimination noise, the exchange current that this system is local Larmor frequency by transmitting coil 6 to underground transmission frequency, produces alternating magnetic field.Under the exciting of the alternating magnetic field forming in ground, make hydrogen proton in underground water produce energy level transition, a large amount of hydrogen protons transit on high level.When removing supply current, these high level hydrogen protons are just got back to low-lying level state gradually, discharge a large amount of quantum of energys with Larmor frequency, in receiving coil 7, reference coil 8 and reference coil 9, induce MRS signal.Acquiescence North and South direction is x direction of principal axis, east-west direction is y direction of principal axis, vertical direction is z direction of principal axis, the particular orientation of laying due to receiving coil 7, reference coil 8 and reference coil 9, can be by the x, y, z component signal that receives MRS signal of selectivity, utilize the correlativity of x and y component signal and z component signal, the time-domain filtering algorithm in use magnetic resonance three-component noise-eliminating method is with reference to eliminating after z component noise, and filtering processing can obtain reliable magnetic resonance signal.
Host computer 1 is connected with control system 2, control system 2 is connected with radiating circuit 3, radiating circuit 3 is connected with transmitting coil 6 by joining humorous electric capacity 5, radiating circuit 3 is connected with power supply 4, control system 2 gathers 12 with A/D and is connected, A/D gathers 12 and is connected with amplifying circuit 11, and amplifying circuit 11 is connected with modulate circuit 10, and modulate circuit 10 is connected with receiving coil 7, reference coil 8 and reference coil 9.
Magnetic resonance three-component noise-eliminating method, comprises following order and step:
In a, Ce district, select to survey place, lay transmitting coil 6;
B, take vertical direction as normal vector lay receiving coil 7, take North and South direction as normal vector lay reference coil 8, take east-west direction as normal vector lay reference coil 9.Receiving coil 7, reference coil 8 and reference coil 9 be orthogonal laying on locus centered by O point, and transmitting coil 6 must be arranged and be laid on same place with receiving coil 7, reference coil 8 and reference coil 9;
C, first measure noise one time, open magnetic resonance three component seismic data acquisition system, host computer 1 sends control word to control system 2, comprises the information that excites duration and stimulating frequency.Control system 2 produces control signal radiating circuit 3 is controlled, and radiating circuit 3 does not produce transmitter current to transmitting coil 6, and this is emitted as pseudo-transmitting.
D, receiving coil 7, reference coil 8 and reference coil 9 are undertaken after signal condition and three-channel amplifying circuit 11 amplify being gathered by A/D acquisition module 12 by receiving after puppet transmitting after signal through three-channel modulate circuit 10;
E, acquiescence North and South direction are x direction of principal axis, and the signal that reference coil 8 receives is x component noise; East-west direction is y direction of principal axis, and the signal that reference coil 9 receives is y component noise; Vertical direction is z direction of principal axis, and the signal that receiving coil 7 receives is z component noise;
F, measurement magnetic resonance signal, now host computer 1 sends control word to control system 2 again, comprises the information that excites duration and stimulating frequency.Control system 2 produces control signal radiating circuit 3 is controlled, and radiating circuit 3 produces by joining humorous electric capacity 5 exchange current that frequency is local Larmor frequency in transmitting coil 6;
G, receiving coil 7, reference coil 8 and reference coil 9 are gathered by A/D acquisition module 12 after sensing signal after three-channel modulate circuit 10 carries out signal condition and 11 amplifications of three-channel amplifying circuit, and the signal now receiving is respectively magnetic resonance x component signal, magnetic resonance y component signal and magnetic resonance z component signal;
H, with the handling object of the clock signal of three components recording simultaneously, time-domain filtering algorithm is proposed:
h z(t)=[h x(t)*X(t)]+[h y(t)*Y(t)]
Wherein h x, h yand h zfor x, y and z component noise, X and Y are x and y component signal, * is convolution algorithm;
I, obtain z component noise according to above-mentioned time-domain filtering algorithm, from z component signal, remove the reliable magnetic resonance signal that is that z component noise obtains;
J, change transmitter current size, obtain different pulses apart from different depth of stratums is surveyed, and repeats c-i process, obtains second group of magnetic resonance signal after processing;
K, repeat above-mentioned j process, carry out 24 times and survey, 24 groups of detection datas are processed.

Claims (3)

1. a magnetic resonance three-component noise cancellation apparatus, it is characterized in that, through control system (2), radiating circuit (3) with join humorous electric capacity (5) and be connected with transmitting coil (6) by host computer (1), radiating circuit (3) is connected with power supply (4), and control system (2) gathers (12), amplifying circuit (11) and modulate circuit (10) through A/D and connects and composes with receiving coil (7), reference coil (8) and reference coil (9) respectively.
2. magnetic resonance three-component noise cancellation apparatus according to claim 1, is characterized in that, centered by point of fixity O along the orthogonal laying reference coil of X, Y, Z axis (8), reference coil (9) and receiving coil (7).
3. a magnetic resonance three-component noise-eliminating method, is characterized in that, comprises the following steps:
In a, Ce district, select to survey place, lay transmitting coil 6;
B, take vertical direction as normal vector lay receiving coil (7), North and South direction be normal vector lay reference coil (8), take east-west direction as normal vector lay reference coil (9); The orthogonal laying on locus centered by O point of receiving coil (7), reference coil (8) and reference coil (9), and transmitting coil (6) is laid on same place with receiving coil (7), reference coil (8) and reference coil (9);
C, first measure noise one time, open magnetic resonance three component seismic data acquisition system, host computer (1) sends control word to control system (2), comprise the information that excites duration and stimulating frequency, control system (2) produces control signal radiating circuit (3) is controlled, radiating circuit (3) does not produce transmitter current to transmitting coil (6), and this is emitted as pseudo-transmitting;
D, receiving coil (7), reference coil (8) and reference coil (9) are gathered the signal receiving after puppet transmitting after three-channel modulate circuit (10) carries out signal condition and three-channel amplifying circuit (11) amplification by A/D acquisition module (12);
E, acquiescence North and South direction are x direction of principal axis, and the signal that reference coil (8) receives is x component noise; East-west direction is y direction of principal axis, and the signal that reference coil (9) receives is y component noise; Vertical direction is z direction of principal axis, and the signal that receiving coil (7) receives is z component noise;
F, measurement magnetic resonance signal, now host computer (1) sends control word to control system (2) again, comprise the information that excites duration and stimulating frequency, control system (2) produces control signal radiating circuit (3) is controlled, and radiating circuit (3) produces by joining humorous electric capacity (5) exchange current that frequency is local Larmor frequency in transmitting coil (6);
G, receiving coil (7), reference coil (8) and reference coil (9) are gathered by A/D acquisition module (12) after sensing signal after three-channel modulate circuit (10) carries out signal condition and three-channel amplifying circuit (11) amplification, and the signal now receiving is respectively magnetic resonance x component signal, magnetic resonance y component signal and magnetic resonance z component signal;
H, with the handling object of the clock signal of three components recording simultaneously, time-domain filtering algorithm is proposed:
h z(t)=[h x(t)*X(t)]+[h y(t)*Y(t)]
Wherein h x, h yand h zfor x, y and z component noise, X and Y are x and y component signal, * is convolution algorithm;
I, obtain z component noise according to above-mentioned time-domain filtering algorithm, from z component signal, remove the reliable magnetic resonance signal that is that z component noise obtains;
J, change transmitter current size, obtain different pulses apart from different depth of stratums is surveyed, and repeats c-i process, obtains second group of magnetic resonance signal after processing;
K, repeat above-mentioned j process, carry out multi collect, the data of multi collect are processed to the magnetic resonance signal obtaining and contrast, eliminate the false and retain the true, obtain accurate information.
CN201410088813.0A 2014-03-11 2014-03-11 Magnetic resonance three-component noise removing device and noise removing method Active CN103823244B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410088813.0A CN103823244B (en) 2014-03-11 2014-03-11 Magnetic resonance three-component noise removing device and noise removing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410088813.0A CN103823244B (en) 2014-03-11 2014-03-11 Magnetic resonance three-component noise removing device and noise removing method

Publications (2)

Publication Number Publication Date
CN103823244A true CN103823244A (en) 2014-05-28
CN103823244B CN103823244B (en) 2015-04-22

Family

ID=50758391

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410088813.0A Active CN103823244B (en) 2014-03-11 2014-03-11 Magnetic resonance three-component noise removing device and noise removing method

Country Status (1)

Country Link
CN (1) CN103823244B (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104297807A (en) * 2014-09-06 2015-01-21 吉林大学 Magnetic resonance imaging device for underground disaster water source detection and detection and imaging method
CN104777442A (en) * 2015-04-07 2015-07-15 吉林大学 MRS (magnetic resonance sounding) FID (frequency identity) signal noise inhibition method
CN105223622A (en) * 2015-09-24 2016-01-06 吉林大学 Antisaturation high-gain adaptive cancellation type ground nuclear magnetic resonance signal amplifying apparatus
CN105652328A (en) * 2016-03-09 2016-06-08 吉林大学 Automatic resonance matching device for nuclear magnetic resonance water detector
CN106019393A (en) * 2016-07-04 2016-10-12 吉林大学 Method for detecting ground nuclear magnetic off-resonance in Larmor frequency unknown condition
CN107533117A (en) * 2015-04-30 2018-01-02 皇家飞利浦有限公司 Method and apparatus for the magnetic resonance imaging with RF noises
CN108828679A (en) * 2018-06-20 2018-11-16 吉林大学 A kind of multi-point remote manipulation nuclear magnetic resonance detection device and method based on 4G network
CN109100813A (en) * 2018-08-14 2018-12-28 吉林大学 A method of it is filtered based on collaboration and eliminates spike noise in ground nuclear magnetic resonance data
CN109143389A (en) * 2018-08-01 2019-01-04 吉林大学 A kind of three-dimensional industrial frequency interference source for nuclear-magnetism quantifies orienting device and measurement method
CN111999774A (en) * 2020-08-12 2020-11-27 吉林大学 Ground-well combined three-dimensional magnetic resonance depth measurement device and method
CN116449438A (en) * 2023-06-16 2023-07-18 中国地质大学(武汉) Transient electromagnetic method reference denoising device and method
WO2023226116A1 (en) * 2022-05-27 2023-11-30 浙江大学 Method for assessing intrinsic time domain stability of radio-frequency receiving coil in fmri

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5578920A (en) * 1994-07-14 1996-11-26 Spectrospin Ag Compensation of interfering fields in NMR measurements in the earth's magnetic field
CN102053280A (en) * 2010-11-10 2011-05-11 吉林大学 Nuclear magnetic resonance ground water detection system with reference coils and detection method
CN202075424U (en) * 2011-03-03 2011-12-14 西安理工大学 Transient electromagnetic method central loop three-component and overlap loop receiving device
CN103033849A (en) * 2012-12-15 2013-04-10 吉林大学 Multichannel nuclear magnetic resonance underground water detecting instrument and field work method thereof

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5578920A (en) * 1994-07-14 1996-11-26 Spectrospin Ag Compensation of interfering fields in NMR measurements in the earth's magnetic field
CN102053280A (en) * 2010-11-10 2011-05-11 吉林大学 Nuclear magnetic resonance ground water detection system with reference coils and detection method
CN202075424U (en) * 2011-03-03 2011-12-14 西安理工大学 Transient electromagnetic method central loop three-component and overlap loop receiving device
CN103033849A (en) * 2012-12-15 2013-04-10 吉林大学 Multichannel nuclear magnetic resonance underground water detecting instrument and field work method thereof

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JUAN PLATA 等: "MRS experiments in a noisy area of a detrital aquifer in the south of Spain", 《JOURNAL OF APPLIED GEOPHYSICS》, 31 December 2002 (2002-12-31), pages 83 - 94 *
田宝凤 等: "基于参考线圈和变步长自适应的磁共振信号噪声压制方法", 《地球物理学报》, vol. 55, no. 7, 31 July 2012 (2012-07-31), pages 2462 - 2472 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104297807B (en) * 2014-09-06 2017-04-26 吉林大学 Magnetic resonance imaging device for underground disaster water source detection and detection and imaging method
CN104297807A (en) * 2014-09-06 2015-01-21 吉林大学 Magnetic resonance imaging device for underground disaster water source detection and detection and imaging method
CN104777442A (en) * 2015-04-07 2015-07-15 吉林大学 MRS (magnetic resonance sounding) FID (frequency identity) signal noise inhibition method
CN104777442B (en) * 2015-04-07 2017-06-16 吉林大学 A kind of nuclear magnetic resonance depth measurement FID signal noise suppressing method
CN107533117B (en) * 2015-04-30 2020-05-01 皇家飞利浦有限公司 Method and apparatus for magnetic resonance imaging with RF noise
CN107533117A (en) * 2015-04-30 2018-01-02 皇家飞利浦有限公司 Method and apparatus for the magnetic resonance imaging with RF noises
CN105223622A (en) * 2015-09-24 2016-01-06 吉林大学 Antisaturation high-gain adaptive cancellation type ground nuclear magnetic resonance signal amplifying apparatus
CN105652328A (en) * 2016-03-09 2016-06-08 吉林大学 Automatic resonance matching device for nuclear magnetic resonance water detector
CN106019393A (en) * 2016-07-04 2016-10-12 吉林大学 Method for detecting ground nuclear magnetic off-resonance in Larmor frequency unknown condition
CN108828679A (en) * 2018-06-20 2018-11-16 吉林大学 A kind of multi-point remote manipulation nuclear magnetic resonance detection device and method based on 4G network
CN109143389A (en) * 2018-08-01 2019-01-04 吉林大学 A kind of three-dimensional industrial frequency interference source for nuclear-magnetism quantifies orienting device and measurement method
CN109143389B (en) * 2018-08-01 2019-07-12 吉林大学 A kind of three-dimensional industrial frequency interference source for nuclear-magnetism quantifies orienting device and measurement method
CN109100813B (en) * 2018-08-14 2019-07-12 吉林大学 A method of it is filtered based on collaboration and eliminates spike noise in ground nuclear magnetic resonance data
CN109100813A (en) * 2018-08-14 2018-12-28 吉林大学 A method of it is filtered based on collaboration and eliminates spike noise in ground nuclear magnetic resonance data
CN111999774A (en) * 2020-08-12 2020-11-27 吉林大学 Ground-well combined three-dimensional magnetic resonance depth measurement device and method
WO2023226116A1 (en) * 2022-05-27 2023-11-30 浙江大学 Method for assessing intrinsic time domain stability of radio-frequency receiving coil in fmri
CN116449438A (en) * 2023-06-16 2023-07-18 中国地质大学(武汉) Transient electromagnetic method reference denoising device and method
CN116449438B (en) * 2023-06-16 2023-09-12 中国地质大学(武汉) Transient electromagnetic method reference denoising device and method

Also Published As

Publication number Publication date
CN103823244B (en) 2015-04-22

Similar Documents

Publication Publication Date Title
CN103823244B (en) Magnetic resonance three-component noise removing device and noise removing method
CN102053280B (en) Nuclear magnetic resonance ground water detection system with reference coils and detection method
CN102426393B (en) Electric prospecting method
CN102565862B (en) Gradient measurement method of transient electromagnetic response signal and observation device thereof
CN103955003B (en) Noise suppression method for superconducting transient electromagnetic application
CN102590869B (en) Artificial field source frequency domain electrical prospecting method and prospecting system
CN103207413B (en) Electrical prospecting device and system
CN203502602U (en) Nuclear magnetic resonance detection device eliminating power frequency harmonic interference
CN110703316B (en) Optical fiber ground seismic wave detection method and system
CN109085653A (en) A kind of detection method of geology of deep part, sulfide ore body resource
CN110187394B (en) Method and device for acquiring formation resistivity anisotropy by double-field source electromagnetic sounding method
CN109490966A (en) A kind of crustal magnetotelluric measurement system
CN104297807A (en) Magnetic resonance imaging device for underground disaster water source detection and detection and imaging method
CN107942397A (en) With the magnetic resonance multi-channel detection method and device of prepolarizing field enhancing signal amplitude
CN102707323A (en) Controllable source audio-frequency magnetic field sounding method for geological exploration
CN104035137A (en) Underground full-space transient electromagnetic detecting instrument and detection method
CN103941298A (en) Transient electromagnetic instrument and hydro-geologic prospecting method for mine
CN103852794A (en) Hydrocarbon polluted shallow groundwater magnetic resonance detection device and hydrocarbon polluted shallow groundwater magnetic resonance detection method
CN102520451B (en) High-precision double-way ground voice frequency electromagnetic physical geography prospecting instrument
CN106291719A (en) A kind of array artificial source's field frequency depth detecting method
Liu et al. Apsu: a wireless multichannel receiver system for surface nuclear magnetic resonance groundwater investigations
CN103486909B (en) Low-frequency electromagnetic induction secondary field imaginary component mine detection device and detection method
CN114217354A (en) Electromagnetic data acquisition system and method based on optical fiber electromagnetic sensor
CN111856601B (en) Distributed magnetic resonance underground water detection device and detection method
Lin et al. Acquisition and correction for early surface nuclear magnetic resonance signal based on multi-turn small air-core coil sensors

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant