CN105334239B - Multi-dimensional nmr fluid components content measuring method and device - Google Patents
Multi-dimensional nmr fluid components content measuring method and device Download PDFInfo
- Publication number
- CN105334239B CN105334239B CN201510766342.9A CN201510766342A CN105334239B CN 105334239 B CN105334239 B CN 105334239B CN 201510766342 A CN201510766342 A CN 201510766342A CN 105334239 B CN105334239 B CN 105334239B
- Authority
- CN
- China
- Prior art keywords
- sample
- echo train
- amplitude
- representing
- train signal
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
Abstract
The embodiment of the present invention provides a kind of multi-dimensional nmr fluid components content measuring method and device.This method includes：First pulse train is applied to sample, so that sample produces the first echo string signal, obtains the amplitude of the first echo string signal；Second pulse train is applied to sample, so that sample produces the second echo string signal, obtains the amplitude of the second echo string signal；The amplitude of amplitude and the second echo string signal to the first echo string signal carries out the content that joint inversion obtains each component in sample.The embodiment of the present invention is by obtaining amplitude of the sample to echo string signal caused by two pulse trains difference, and two echo string signals are carried out with the content that joint inversion obtains each component in sample, the measurement number for applying pulse train in the sample compared to prior art is 5 10, improves the detection speed of the content of each component in sample.
Description
Technical Field
The embodiment of the invention relates to the field of petroleum exploration, in particular to a method and a device for measuring the component content of a multi-dimensional nuclear magnetic resonance fluid.
Background
In the field of oil exploration, a nuclear magnetic resonance apparatus is an important analysis testing tool, a pulse sequence generated by the nuclear magnetic resonance apparatus is applied to a sample, so that the sample generates an echo train signal for the pulse sequence, and the content of each component in the sample can be analyzed and obtained according to the amplitude of the echo train signal.
In the prior art, a spectrometer of a nuclear magnetic resonance spectrometer generates a pulse sequence, a sample is placed in a probe of the nuclear magnetic resonance spectrometer, an upper computer of the nuclear magnetic resonance spectrometer controls the pulse sequence generated by the spectrometer to be applied to the sample so as to enable the sample to generate echo train signals, 5-10 times of measurement is needed, joint inversion is carried out on the echo train signals to obtain a two-dimensional nuclear magnetic resonance spectrum, and the content of each component in the sample is obtained according to analysis of the two-dimensional nuclear magnetic resonance spectrum.
Because the number of the pulse sequences applied to the sample is 5-10 or more, the sample generates an echo train signal for each pulse sequence, and a two-dimensional nuclear magnetic resonance spectrum can be obtained only by joint inversion of a plurality of echo train signals, so that the measurement speed of the content of each component in the sample is low.
Disclosure of Invention
The embodiment of the invention provides a method and a device for measuring the component content of a multi-dimensional nuclear magnetic resonance fluid, which are used for improving the speed of measuring the content of each component in a sample.
One aspect of an embodiment of the present invention is to provide a method for measuring a component content of a multi-dimensional nuclear magnetic resonance fluid, including:
applying a first pulse sequence to a sample to enable the sample to generate a first echo train signal, and acquiring the amplitude of the first echo train signal;
applying a second pulse sequence to the sample to enable the sample to generate a second echo train signal, and acquiring the amplitude of the second echo train signal;
performing joint inversion on the amplitude of the first echo train signal and the amplitude of the second echo train signal to obtain the content of each component in the sample;
wherein the first pulse sequence or the second pulse sequence affects different physical properties of the sample, including longitudinal relaxation time, transverse relaxation time, and diffusion coefficient.
Another aspect of the embodiments of the present invention is to provide a multi-dimensional nmr fluid component content measuring apparatus, including:
a pulse sequence applying module, configured to apply a first pulse sequence to a sample, so that the sample generates a first echo train signal; applying a second pulse sequence to the sample to cause the sample to generate a second echo train signal;
the amplitude acquisition module is used for acquiring the amplitude of the first echo train signal; obtaining the amplitude of the second echo string signal;
the joint inversion module is used for performing joint inversion on the amplitude of the first echo string signal and the amplitude of the second echo string signal to obtain the content of each component in the sample;
wherein the first pulse sequence or the second pulse sequence affects different physical properties of the sample, including longitudinal relaxation time, transverse relaxation time, and diffusion coefficient.
According to the method and the device for measuring the content of the components in the multi-dimensional nuclear magnetic resonance fluid, provided by the embodiment of the invention, the two pulse sequences are sequentially applied to the sample, one of the two pulse sequences influences different physical properties of the sample, the amplitude of echo train signals respectively generated by the sample on the two pulse sequences is obtained, and the combined inversion is carried out on the amplitudes of the two echo train signals to obtain the content of each component in the sample, so that the measurement frequency is 5-10, and the speed of measuring the content of each component in the sample is improved.
Drawings
FIG. 1 is a flow chart of a method for measuring the component content of a multi-dimensional NMR fluid according to an embodiment of the invention;
FIG. 2 is a schematic diagram of a fast measurement T2-D pulse sequence provided by an embodiment of the present invention;
FIG. 3 is a schematic diagram of a two-shot T1-T2 pulse sequence provided by an embodiment of the present invention;
FIG. 4 is a schematic diagram of a single fast measurement T1-T2 pulse sequence according to an embodiment of the present invention;
FIG. 5 is a two-dimensional NMR T provided by an embodiment of the invention 2 -a D-atlas model;
FIG. 6 shows a graph of T according to another embodiment of the present invention 2 -a schematic representation of a D-atlas forward-evolving echo train;
FIG. 7 shows a two-dimensional NMR T according to another embodiment of the invention 2 -D inversion map;
fig. 8 is a structural diagram of a component content measuring apparatus for a multi-dimensional nmr fluid according to an embodiment of the invention.
Detailed Description
FIG. 1 is a flow chart of a method for measuring the component content of a multi-dimensional NMR fluid according to an embodiment of the invention; FIG. 2 is a schematic diagram of a fast measurement T2-D pulse sequence provided by an embodiment of the present invention; FIG. 3 is a schematic diagram of a two-shot T1-T2 pulse sequence provided by an embodiment of the present invention; FIG. 4 is a schematic diagram of a single fast measurement T1-T2 pulse sequence according to an embodiment of the present invention; FIG. 5 shows a two-dimensional NMR T 2 -a D-atlas model. The embodiment of the invention provides a method for measuring the component content of a multidimensional nuclear magnetic resonance fluid, aiming at the problems that the number of pulse sequences applied to a sample is 5-10 or more, the sample generates an echo train signal for each pulse sequence, and a two-dimensional nuclear magnetic resonance spectrum can be obtained only by joint inversion of a plurality of echo train signals, so that the measuring speed of the content of each component in the sample is slow, and the method comprises the following specific steps:
step S101, a first pulse sequence is applied to a sample to enable the sample to generate a first echo train signal, and the amplitude of the first echo train signal is obtained;
the sample is placed in advance in the probe of the NMR spectrometer, and the magnet in the probe generates B 0 A magnetic field of B 0 The magnetic field polarizes the hydrogen atoms in the sample, so that the spin orientation of the hydrogen atoms is transited from a disordered state before being placed in the magnetic field to an ordered state, and energy level transition is generated, and the hydrogen atoms in the sample are completely polarized to form a macroscopic magnetization vector along with the extension of the time for placing the sample in the probe, namely the extension of the polarization time.
The spectrometer of the nuclear magnetic resonance apparatus generates a pulse sequence having the same frequency as the spin precession, and the upper computer of the nuclear magnetic resonance apparatus controls the application of the pulse sequence generated by the spectrometer to the sample to cause the sample to generate an echo train signal. In the embodiment of the present invention, only two pulse sequences, i.e., a first pulse sequence and a second pulse sequence, are applied to the sample, specifically, the first pulse sequence is the sequence a in fig. 2, and the second pulse sequence is the sequence b in fig. 2; or the first pulse sequence is a sequence a in fig. 3, and the second pulse sequence is b sequence in fig. 3; or the first pulse sequence is the DE sequence in fig. 4 and the second pulse sequence is the CPMG sequence in fig. 4.
The first pulse sequence applied to the sample by the NMR spectrometer is the sequence a in FIG. 2, and the magnetic field of the sequence a in FIG. 2 has no gradient, then the sequence a in FIG. 2 influences the transverse relaxation time T of the sample 2 Generating a first echo train signal for the first pulse sequence by the sample, detecting and receiving the first echo train signal by the nuclear magnetic resonance spectrometer, and simultaneously obtaining the amplitude of the first echo train signal; the second pulse sequence applied to the sample by the NMR spectrometer is the sequence b in FIG. 2, and the gradient G exists in the magnetic field of the sequence b in FIG. 2, then the sequence b in FIG. 2 influences the transverse relaxation time T of the sample 2 And a diffusion coefficient D.
If the first pulse sequence applied to the sample by the NMR spectrometer is the sequence a in FIG. 3, and the sequence a in FIG. 3 influences the longitudinal relaxation time T of the sample 1 And transverse relaxation time T 2 The second pulse sequence applied by the NMR spectrometer to the sample is the sequence b in FIG. 3, and the sequence b in FIG. 3 influences the transverse relaxation time T of the sample 2 Then, the two-dimensional nuclear magnetic resonance T is obtained by the method of the embodiment of the invention 1 -T 2 And (4) mapping.
Similarly, if the first pulse sequence applied to the sample by the NMR spectrometer is the DE sequence in FIG. 4, and the DE sequence in FIG. 4 affects the longitudinal relaxation time T of the sample 1 And transverse relaxation time T 2 The second pulse sequence applied to the sample by the NMR spectrometer is the CPMG sequence in FIG. 4, and the CPMG sequence in FIG. 4 influences the transverse relaxation time T of the sample 2 Then a two-dimensional NMR T will be obtained by the method of the embodiment of the invention 1 -T 2 And (4) mapping.
A first pulse sequence is firstly applied to a sample so as to enable the sample to generate a first echo train signal, and the amplitude of the first echo train signal is obtained.
Step S102, a second pulse sequence is applied to the sample, so that the sample generates a second echo train signal, and the amplitude of the second echo train signal is obtained;
after step S101, a second pulse sequence is applied to the sample, so that a second echo train signal is generated by the sample, and an amplitude of the second echo train signal is obtained.
Before the first pulse sequence is applied to the sample, the hydrogen atoms in the sample are completely polarized to form a macroscopic magnetization vector, when the first pulse sequence is applied to the sample, the spins of the hydrogen atoms in the sample are converted from an equilibrium state to a non-equilibrium state after the complete polarization, the magnetization vector is switched, a specific 90-degree pulse pulls the macroscopic magnetization vector to a transverse plane, and then the half-echo time tau is CPMG The internal phase dispersion occurs, 180-degree pulses reunite the magnetization vectors of the dispersed phase to finally form echo, and the 180-degree pulses are continuously applied to the magnetization vectors after the phase dispersion to generate a CPMG echo train. The spins of the hydrogen atoms in the sample transition from a non-equilibrium state to an equilibrium state for a period of time after the first pulse sequence is ended and reach an equilibrium state before the second pulse sequence is applied.
Similarly, when the second pulse sequence is applied to the sample, the spins of the hydrogen atoms in the sample are changed from the equilibrium state to the non-equilibrium state after complete polarization, the magnetization vector is switched, the macroscopic magnetization vector is switched to the transverse plane by the specific 90-degree pulse, and the following half-echo time tau is obtained CPMG Generating a dispersed phase inside, reuniting the magnetization vector of the dispersed phase by 180-degree pulse, finally forming echo, continuously applying 180-degree pulse to the dispersed magnetization vector to generate a CPMG echo train, and converting the spin of hydrogen atoms in the sample from a non-equilibrium state to an equilibrium state in a period of time after the second pulse sequence is finished. The difference is that in this case the second pulse sequence is influenced by a magnetic field gradient.
Step S103, performing joint inversion on the amplitude of the first echo train signal and the amplitude of the second echo train signal to obtain the content of each component in the sample;
wherein the first pulse sequence or the second pulse sequence affects different physical properties of the sample, including longitudinal relaxation time, transverse relaxation time, and diffusion coefficient.
And performing joint inversion on the amplitude of the acquired first echo string signal and the amplitude of the acquired second echo string signal to obtain a two-dimensional nuclear magnetic resonance spectrum, and obtaining the content of each component in the sample according to the two-dimensional nuclear magnetic resonance spectrum. The sample has three important physical properties, namely the longitudinal relaxation time T 1 Transverse relaxation time T 2 And a diffusion coefficient D.
For example, the sample generates a second echo train signal for the second pulse sequence, and the nmr detects and receives the second echo train signal, and obtains an amplitude of the second echo train signal; performing joint inversion on the amplitude of the first echo train signal and the amplitude of the second echo train signal to obtain a two-dimensional nuclear magnetic resonance T shown in FIG. 7 2 -D map according to the two-dimensional NMR T 2 The D-profile allows the content of the components in the sample to be obtained by means of the prior art.
According to the embodiment of the invention, two pulse sequences are applied to the sample in sequence, one of the two pulse sequences influences different physical properties of the sample, the amplitudes of echo train signals generated by the sample on the two pulse sequences respectively are obtained, and the amplitudes of the two echo train signals are subjected to joint inversion to obtain the content of each component in the sample, compared with the prior art that the number of the pulse sequences applied to the sample is 5-10, the speed of measuring the content of each component in the sample is improved.
On the basis of the above embodiment, the first pulse sequence affects the transverse relaxation time of the sample, and the second pulse sequence affects the transverse relaxation time and the diffusion coefficient of the sample.
Echo train signal in the embodiment of the present invention, by analysing two-dimensional nuclear magnetic resonance T 1 -T 2 Atlas or two-dimensional NMR T 2 The content of each component in the sample can be obtained by a D map。
According to the embodiment of the invention, different physical characteristics of the sample can be influenced by applying different two pulse sequences to the sample in sequence, so that different two-dimensional nuclear magnetic resonance maps can be obtained, the content of each component in the sample can be obtained according to different two-dimensional nuclear magnetic resonance maps, and methods for measuring the content of each component in the sample are enriched.
Based on the above embodiments, the amplitude M of the first echo train signal C (nT E ) Expressed as formula (1):
wherein, 0<n<N c ，N c Represents the total number of the first echo train signals, N represents N c An nth one of the first echo train signals, T E Representing the echo interval, N representing the number of points of transverse relaxation time, N 1 Number of points, T, representing diffusion coefficient 2i For transverse relaxation time distribution, f ij Denotes the content of a component of the component type i, j, A 1 Representing a kernel function matrix corresponding to the first echo train signal;
amplitude M of the second echo train signal D (mT E ) Expressed as formula (2):
wherein, 0<m<N d ，N d Representing the total number of said second echo train signals, m representing N d M-th of the second echo train signals, gamma represents gyromagnetic ratio, G represents static magnetic field gradient, and D j Denotes the jth diffusion coefficient, A 2 Representing a kernel function matrix corresponding to the second echo train signal;
performing joint inversion on the amplitude of the first echo train signal and the amplitude of the second echo train signal to obtain the content of each component in the sample, including:
performing joint inversion according to the formula (1) and the formula (2) to obtain a formula (3):
wherein M is C An amplitude M representing the first echo train signal C (nT E )，M D An amplitude M representing the second echo train signal D (mT E )；
Obtaining the content f of the component with the component type i, j according to the formula (3) ij 。
As the first pulse sequence applied to the sample by the nmr according to the above embodiment is the sequence a in fig. 2, since the first pulse sequence is a discrete sequence, the first echo train signal generated by the sample to the first pulse sequence is also a discrete signal, and the amplitude of the first echo train signal obtained by the nmr is also a discrete value, which is specifically expressed by formula (1):
wherein, 0<n<N c ，N c Represents the total number of the first echo train signals, N represents N c An nth one of the first echo train signals, T E Representing the echo interval, N representing the number of points of transverse relaxation time, N 1 Number of points, T, representing diffusion coefficient 2i For transverse relaxation time distribution, f ij Denotes the content of a component of component type i, j, A 1 Representing a kernel function matrix corresponding to the first echo train signal.
The second pulse sequence applied to the sample by the nmr is the sequence b in fig. 2, and since the second pulse sequence is a discrete sequence, the second echo train signal generated by the sample to the second pulse sequence is also a discrete signal, and the amplitude of the second echo train signal obtained by the nmr is also a discrete value, which is specifically expressed as formula (2):
wherein, 0<m<N d ，N d Representing the total number of said second echo train signals, m representing N d M-th second echo train signal among the second echo train signals, gamma represents a gyromagnetic ratio, G represents a static magnetic field gradient, and D j Denotes the jth diffusion coefficient, A 2 Representing a kernel function matrix corresponding to the second echo train signal.
Performing joint inversion on the formula (1) and the formula (2) to obtain a formula (3):
wherein M is C An amplitude M representing the first echo train signal C (nT E )，M D An amplitude M representing the second echo train signal D (mT E )；
In formula (3), the amplitude M of the first echo train signal C (nT E ) Can be detected to obtain C For a known quantity, the amplitude M of the second echo train signal is likewise D (mT E ) Can also be detected to obtain D A kernel function matrix A corresponding to the first echo train signal 1 A kernel function matrix A corresponding to the second echo train signal 2 Are all known, the content f of the component type i, j is obtained according to the formula (3) ij 。
Preferably, the first and second liquid crystal materials are,B＝S*V T (ii) a Performing singular value decomposition on A to obtain A = U S V T (ii) a According to A = U S V T ，B＝S*V T And formula (3) to obtain Bf ij = M; according to B x f ij = M obtaining content f of component type i, j ij 。
In the embodiment of the present invention, it is,a is an ultra-large kernel function, which makes the inversion speed slow. In order to improve the inversion rate, singular value decomposition is carried out on the kernel function A by adopting an SVD (singular value decomposition) method to obtain A = U S V T The embodiment of the inventionB＝S*V T A = U × S × V T ， B＝S*V T Substituting equation (3) to obtain B f ij = M, wherein B and M are known quantities in terms of B f ij = M obtaining content f of component type i, j ij 。
The embodiment of the invention provides a formula for carrying out joint inversion on the amplitude values of two echo train signals to obtain the content of each component in a sample, and the calculation precision of the content of each component in the sample is improved.
FIG. 6 is a schematic diagram of a forward echo train according to another embodiment of the present invention; FIG. 7 shows a two-dimensional NMR T according to another embodiment of the invention 2 -D inversion map; shown in FIG. 6 is a graph according to T shown in FIG. 5 2 -obtaining a first echo train signal and a second echo train signal by forward modeling of a D-atlas model, and inverting the first echo train signal and the second echo train signal shown in fig. 6 to obtain T shown in fig. 7 2 -D inversion map. T can be seen by comparing FIGS. 5 and 7 2 -D atlas model and T 2 Similarity of inversion maps of-D, illustrating the multidimensional NMR provided by the embodiments of the inventionThe fluid component content measurement method is correct.
Fig. 8 is a structural diagram of a component content measuring apparatus for a multi-dimensional nmr fluid according to an embodiment of the invention. The multi-dimensional nuclear magnetic resonance fluid component content measuring device provided by the embodiment of the present invention can execute the processing flow provided by the multi-dimensional nuclear magnetic resonance fluid component content measuring method embodiment, as shown in fig. 8, the multi-dimensional nuclear magnetic resonance fluid component content measuring device 90 includes a pulse sequence applying module 91, an amplitude obtaining module 92, and a joint inversion module 93, where the pulse sequence applying module 91 is configured to apply a first pulse sequence to a sample, so that the sample generates a first echo train signal; applying a second pulse sequence to the sample to cause the sample to generate a second echo train signal; the amplitude obtaining module 92 is configured to obtain an amplitude of the first echo train signal; obtaining the amplitude of the second echo string signal; the joint inversion module 93 is configured to perform joint inversion on the amplitude of the first echo train signal and the amplitude of the second echo train signal to obtain the content of each component in the sample; wherein the first pulse sequence or the second pulse sequence affects different physical properties of the sample, including longitudinal relaxation time, transverse relaxation time, and diffusion coefficient.
According to the embodiment of the invention, two pulse sequences are applied to the sample in sequence, one of the two pulse sequences influences different physical properties of the sample, the amplitudes of echo train signals generated by the sample on the two pulse sequences respectively are obtained, and the amplitudes of the two echo train signals are subjected to joint inversion to obtain the content of each component in the sample, compared with the prior art that the number of the pulse sequences applied to the sample is 5-10, the speed of measuring the content of each component in the sample is improved.
On the basis of the above embodiment, the first pulse sequence affects the transverse relaxation time of the sample, and the second pulse sequence affects the transverse relaxation time and the diffusion coefficient of the sample.
The first pulse sequence affects a longitudinal relaxation time and a transverse relaxation time of the sample, and the second pulse sequence affects a transverse relaxation time of the sample.
Amplitude M of the first echo train signal C (nT E ) Expressed as formula (1):
wherein 0<n<N c ，N c Represents the total number of the first echo train signals, N represents N c An nth one of the first echo train signals, T E Representing the echo interval, N representing the number of points of transverse relaxation time, N 1 Number of points, T, representing diffusion coefficient 2i For transverse relaxation time distribution, f ij Denotes the content of a component of component type i, j, A 1 Representing a kernel function matrix corresponding to the first echo train signal;
amplitude M of the second echo train signal D (mT E ) Expressed as formula (2):
wherein 0<m<N d ，N d Representing the total number of said second echo train signals, m representing N d M-th of the second echo train signals, gamma represents gyromagnetic ratio, G represents static magnetic field gradient, and D j Denotes the jth diffusion coefficient, A 2 Representing a kernel function matrix corresponding to the second echo train signal;
the joint inversion module 93 is specifically configured to perform joint inversion according to the formula (1) and the formula (2) to obtain a formula (3):
wherein M is C An amplitude M representing the first echo train signal C (nT E )，M D An amplitude M representing the second echo train signal D (mT E ) (ii) a Obtaining the content f of the component with the component type i, j according to the formula (3) ij 。
B＝S*V T (ii) a The joint inversion module 93 is further configured to perform singular value decomposition on a to obtain a = U × S × V T (ii) a According to A = U S V T ，B＝S*V T And formula (3) to obtain Bf ij = M; according to B x f ij = M obtaining content f of component type i, j ij 。
The device for measuring the component content of the multi-dimensional nuclear magnetic resonance fluid provided by the embodiment of the present invention may be specifically configured to perform the method embodiment provided in fig. 1, and specific functions are not described herein again.
According to the embodiment of the invention, different physical characteristics of the sample can be influenced by applying different two pulse sequences to the sample in sequence, so that different two-dimensional nuclear magnetic resonance spectrums can be obtained, the content of each component in the sample can be obtained according to different two-dimensional nuclear magnetic resonance spectrums, and the method for measuring the content of each component in the sample is enriched; the formula for carrying out joint inversion on the amplitudes of the two echo train signals to obtain the content of each component in the sample is provided, and the calculation accuracy of the content of each component in the sample is improved.
In summary, in the embodiment of the present invention, two pulse sequences are sequentially applied to a sample, and one of the two pulse sequences affects different physical properties of the sample, the amplitudes of echo train signals respectively generated by the sample on the two pulse sequences are obtained, and the amplitudes of the two echo train signals are jointly inverted to obtain the content of each component in the sample, so that compared with the prior art in which the number of pulse sequences applied to the sample is 5-10, the speed of measuring the content of each component in the sample is increased; different physical characteristics of the sample can be influenced by applying two different pulse sequences to the sample in sequence, so that different two-dimensional nuclear magnetic resonance maps can be obtained, the content of each component in the sample can be obtained according to the different two-dimensional nuclear magnetic resonance maps, and methods for measuring the content of each component in the sample are enriched; the formula for carrying out joint inversion on the amplitudes of the two echo train signals to obtain the content of each component in the sample is provided, and the calculation accuracy of the content of each component in the sample is improved.
In the embodiments provided in the present invention, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, or in a form of hardware plus a software functional unit.
The integrated unit implemented in the form of a software functional unit may be stored in a computer readable storage medium. The software functional unit is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) or a processor (processor) to execute some steps of the methods according to the embodiments of the present invention. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
It is obvious to those skilled in the art that, for convenience and simplicity of description, the foregoing division of the functional modules is merely used as an example, and in practical applications, the above function distribution may be performed by different functional modules according to needs, that is, the internal structure of the device is divided into different functional modules to perform all or part of the above described functions. For the specific working process of the device described above, reference may be made to the corresponding process in the foregoing method embodiment, which is not described herein again.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and these modifications or substitutions do not depart from the spirit of the corresponding technical solutions of the embodiments of the present invention.
Claims (4)
1. A method for measuring the component content of a multi-dimensional nuclear magnetic resonance fluid is characterized by comprising the following steps:
applying a first pulse sequence to a sample to enable the sample to generate a first echo train signal, and acquiring the amplitude of the first echo train signal;
applying a second pulse sequence to the sample to enable the sample to generate a second echo train signal, and acquiring the amplitude of the second echo train signal;
performing joint inversion on the amplitude of the first echo train signal and the amplitude of the second echo train signal to obtain the content of each component in the sample;
wherein the first pulse sequence or the second pulse sequence affects different physical properties of the sample, including longitudinal relaxation time, transverse relaxation time, and diffusion coefficient;
the first pulse sequence affects the transverse relaxation time of the sample, and the second pulse sequence affects the transverse relaxation time and diffusion coefficient of the sample;
amplitude M of the first echo train signal C (nT E ) Expressed as formula (1):
wherein N is more than 0 and less than N c ，N c Representing the total number of the first echo train signals, N representing N c N-th of the first echo train signals, T E Representing the echo interval, N representing the number of points of transverse relaxation time, N 1 Number of points, T, representing diffusion coefficient 2i For transverse relaxation time distribution, f ij Denotes the content of a component of component type i, j, A 1 Representing a kernel function matrix corresponding to the first echo train signal;
amplitude M of the second echo train signal D (mT E ) Expressed as formula (2):
wherein m is more than 0 and less than N d ，N d Representing the total number of said second echo train signals, m representing N d M-th of the second echo train signals, gamma represents gyromagnetic ratio, G represents static magnetic field gradient, and D j Denotes the jth diffusion coefficient, A 2 Representing a kernel function matrix corresponding to the second echo train signal;
performing joint inversion on the amplitude of the first echo train signal and the amplitude of the second echo train signal to obtain the content of each component in the sample, including:
performing joint inversion according to the formula (1) and the formula (2) to obtain a formula (3):
wherein, M C Representing the amplitude M of the first echo train signal C (nT E )，M D An amplitude M representing the second echo train signal D (mT E )；
Obtaining the content f of the component with the component type i, j according to the formula (3) ij 。
2. The method of claim 1,B＝S*V T ；
the method further comprises the following steps:
performing singular value decomposition on A to obtain A = U S V T ；
According to A = U S V T ，B＝S*V T And formula (3) to obtain B f ij ＝M；
According to B x f ij = M obtaining content f of component type i, j ij 。
3. A multi-dimensional nuclear magnetic resonance fluid component content measuring device is characterized by comprising:
a pulse sequence applying module, configured to apply a first pulse sequence to a sample so that the sample generates a first echo train signal; applying a second pulse sequence to the sample to cause the sample to generate a second echo train signal;
an amplitude acquisition module, configured to acquire an amplitude of the first echo train signal; acquiring the amplitude of the second echo string signal;
the joint inversion module is used for performing joint inversion on the amplitude of the first echo string signal and the amplitude of the second echo string signal to obtain the content of each component in the sample;
wherein the first pulse sequence or the second pulse sequence affects different physical properties of the sample, including longitudinal relaxation time, transverse relaxation time, and diffusion coefficient;
the first pulse sequence affects the transverse relaxation time of the sample, and the second pulse sequence affects the transverse relaxation time and diffusion coefficient of the sample;
amplitude M of the first echo train signal C (nT E ) Expressed as formula (1):
wherein N is more than 0 and less than N c ，N c Representing the total number of the first echo train signals, N representing N c N-th of the first echo train signals, T E Representing the echo interval, N representing the number of points of transverse relaxation time, N 1 Number of points, T, representing diffusion coefficient 2i For transverse relaxation time distribution, f ij Denotes the content of a component of the component type i, j, A 1 Representing a kernel function matrix corresponding to the first echo train signal;
amplitude M of the second echo train signal D (mT E ) Expressed as formula (2):
wherein m is more than 0 and less than N d ，N d Representing the total number of said second echo train signals, m representing N d The mth second echo train signal in the second echo train signals, gamma tableShowing gyromagnetic ratio, G showing static magnetic field gradient, D j Denotes the jth diffusion coefficient, A 2 Representing a kernel function matrix corresponding to the second echo train signal;
the joint inversion module is specifically configured to perform joint inversion according to a formula (1) and a formula (2) to obtain a formula (3):
wherein M is C Representing the amplitude M of the first echo train signal C (nT E )，M D An amplitude M representing the second echo train signal D (mT E ) (ii) a Obtaining the content f of the component with the component type i, j according to the formula (3) ij 。
4. The multi-dimensional NMR fluid component content measuring apparatus according to claim 3,B＝S*V T ；
the joint inversion module is further used for performing singular value decomposition on A to obtain A = U S V T (ii) a According to A = U × S V T ，B＝S*V T And formula (3) to obtain Bf ij = M; according to B x f ij = M obtaining content f of component type i, j ij 。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510766342.9A CN105334239B (en) | 2015-11-11 | 2015-11-11 | Multi-dimensional nmr fluid components content measuring method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510766342.9A CN105334239B (en) | 2015-11-11 | 2015-11-11 | Multi-dimensional nmr fluid components content measuring method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105334239A CN105334239A (en) | 2016-02-17 |
| CN105334239B true CN105334239B (en) | 2018-02-02 |
Family
ID=55284892
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510766342.9A Active CN105334239B (en) | 2015-11-11 | 2015-11-11 | Multi-dimensional nmr fluid components content measuring method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105334239B (en) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106908370B (en) * | 2017-03-22 | 2019-08-06 | 中国石油天然气股份有限公司 | Seep determination method for parameter and device in sandstone hole |
| CN109270107B (en) * | 2017-03-27 | 2021-10-08 | 北京青檬艾柯科技有限公司 | Multi-dimensional nuclear magnetic resonance measurement method |
| US10724975B2 (en) * | 2017-08-11 | 2020-07-28 | Schlumberger Technology Corporation | Apparatus and methods for determining properties of liquid-bearing solids using nuclear magnetic resonance |
| US10890685B2 (en) * | 2017-08-11 | 2021-01-12 | Schlumberger Technology Corporation | Apparatus and methods for determining properties of hydrogen-containing samples using nuclear magnetic resonance |
| CN107525553B (en) * | 2017-09-19 | 2019-09-06 | 中国石油天然气股份有限公司 | A kind of method and device of determining heterogeneous fluid composition flow rate |
| CN107966465B (en) * | 2017-11-23 | 2019-12-03 | 中国石油大学(华东) | A kind of acquisition of rock core nuclear magnetic signal and inversion method based on three-dimensional pulse train |
| CN108005633A (en) * | 2017-11-27 | 2018-05-08 | 中石化石油工程技术服务有限公司 | The two-dimensional NMR Logging observation mode parameter determination method of tight gas reservoir |
| CN108426909B (en) * | 2018-02-11 | 2019-12-27 | 中国石油大学(北京) | Method and device for synchronously and rapidly measuring nuclear magnetic resonance relaxation time and flow velocity |
| CN108872291B (en) * | 2018-04-25 | 2019-12-27 | 中国石油大学(北京) | Rapid localized multidimensional Laplace nuclear magnetic resonance method and device |
| CN108680212B (en) * | 2018-05-18 | 2020-01-07 | 中国石油天然气股份有限公司 | Multiphase flow magnetic resonance flowmeter calibration device and water content and flow rate calibration method thereof |
| CN109507221B (en) * | 2018-12-03 | 2019-12-13 | 中国科学院生物物理研究所 | Multidimensional nuclear magnetic resonance molecular diffusion coupling imaging method |
| CN111721795B (en) * | 2020-06-29 | 2021-05-11 | 无锡鸣石峻致医疗科技有限公司 | Material measurement method and system based on nuclear magnetic resonance system |
| CN112683940B (en) * | 2020-12-24 | 2022-07-19 | 潍坊欣龙生物材料有限公司 | Method for rapidly determining content of flame retardant in inorganic flame-retardant cellulose fiber |
| CN112903567B (en) * | 2021-04-15 | 2022-06-17 | 中国石油大学(北京) | Aperture measurement method, device, equipment and storage medium |
-
2015
- 2015-11-11 CN CN201510766342.9A patent/CN105334239B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105334239A (en) | 2016-02-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105334239B (en) | Multi-dimensional nmr fluid components content measuring method and device | |
| Buonincontri et al. | MR fingerprinting with simultaneous B1 estimation | |
| CN108027409B (en) | Time domain MRI | |
| Kiselev | Fundamentals of diffusion MRI physics | |
| CN109270107B (en) | Multi-dimensional nuclear magnetic resonance measurement method | |
| Aganj et al. | Reconstruction of the orientation distribution function in single‐and multiple‐shell q‐ball imaging within constant solid angle | |
| Drobnjak et al. | The matrix formalism for generalised gradients with time-varying orientation in diffusion NMR | |
| Lawrenz et al. | A tensor model and measures of microscopic anisotropy for double-wave-vector diffusion-weighting experiments with long mixing times | |
| US10168444B2 (en) | Data processing with magnetic resonance tool | |
| US20130325408A1 (en) | Two dimensional nmr of diffusion and relaxation for material characterization | |
| AU773421B2 (en) | Methods and apparatus for mapping internal and bulk motion of an object with phase labeling in magnetic resonance imaging | |
| JP7218288B2 (en) | How to perform diffusion-weighted magnetic resonance measurements on a sample | |
| Awojoyogbe et al. | Mathematical concept of the Bloch flow equations for general magnetic resonance imaging: A review | |
| CN109254255A (en) | Nuclear magnetic resonance lateral relaxation time spectrometry method, apparatus and readable storage medium storing program for executing | |
| Liu et al. | A novel pulse sequence and inversion algorithm of three-dimensional low field NMR technique in unconventional resources | |
| Bi et al. | Stabilization of parameter estimates from multiexponential decay through extension into higher dimensions | |
| US10871539B2 (en) | Determination of a joint probability distribution of radius and length of anisotropic pores from double pulsed field gradient MRI data | |
| CN109507221B (en) | Multidimensional nuclear magnetic resonance molecular diffusion coupling imaging method | |
| Sigmund et al. | Multiple echo diffusion tensor acquisition technique | |
| Özarslan et al. | Characterizing magnetic resonance signal decay due to Gaussian diffusion: the path integral approach and a convenient computational method | |
| CN107727678A (en) | A kind of intrinsic mode coupling process of NMR relaxation height | |
| WO1994003823A1 (en) | A single shot magnetic resonance method to measure diffusion, flow and/or motion | |
| US10302733B2 (en) | NMR spin-echo amplitude estimation | |
| Lamy et al. | Sycomore-an MRI simulation toolkit | |
| Sørgård et al. | Investigating pore geometry in heterogeneous porous samples using spatially resolved G0-Δχapp and G0-Δν correlations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |