CN115201925A - Two-dimensional spectrum identification method, device, equipment and medium based on nuclear magnetic resonance - Google Patents
Two-dimensional spectrum identification method, device, equipment and medium based on nuclear magnetic resonance Download PDFInfo
- Publication number
- CN115201925A CN115201925A CN202210581764.9A CN202210581764A CN115201925A CN 115201925 A CN115201925 A CN 115201925A CN 202210581764 A CN202210581764 A CN 202210581764A CN 115201925 A CN115201925 A CN 115201925A
- Authority
- CN
- China
- Prior art keywords
- relaxation
- fluid component
- nuclear magnetic
- magnetic resonance
- determining
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 102
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 73
- 239000012530 fluid Substances 0.000 claims abstract description 135
- 238000009826 distribution Methods 0.000 claims abstract description 39
- 230000010354 integration Effects 0.000 claims abstract description 19
- 238000005084 2D-nuclear magnetic resonance Methods 0.000 claims description 4
- 238000004043 dyeing Methods 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 abstract description 5
- 230000000877 morphologic effect Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000655 nuclear magnetic resonance spectrum Methods 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/32—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electron or nuclear magnetic resonance
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses a two-dimensional spectrum identification method and device based on nuclear magnetic resonance, electronic equipment and a storage medium, wherein the method comprises the following steps: acquiring peak positions of all fluid components in a relaxation total signal of a nuclear magnetic resonance two-dimensional spectrum; determining the boundary position of the relaxation distribution range of each fluid component by taking each peak position as a central point, and determining the relaxation signal of each fluid component according to the peak position and the boundary position; determining the relaxation signal intensity and the relaxation total signal intensity of each fluid component based on a preset integration method; determining the saturation of each fluid component based on the ratio of the relaxation signal intensity of each fluid component to the intensity of the relaxation total signal. The invention solves the technical problems of low component decomposition efficiency and poor precision based on a two-dimensional spectrum in the prior art.
Description
Technical Field
The invention relates to the technical field of nuclear magnetic resonance logging two-dimensional spectrums, in particular to a method and a device for identifying a two-dimensional spectrum based on nuclear magnetic resonance, electronic equipment and a storage medium.
Background
Nmr logs are produced by measuring nuclear hydrogen signals in formations to obtain relaxation signals that contain various fluid components in the formation. But because of the transverse relaxation time T of oil, gas and water in the reservoir 2 There is some overlap, resulting in the conventional one-dimensional NMR logging method passing only T 2 The relaxation distribution has certain difficulty in identifying the type of fluid in a reservoir, so that a learner introduces a two-dimensional nuclear magnetic resonance technology in the wave spectrum into the field of well logging. Two-dimensional NMR logging methods are compared to one-dimensional NMR logging methods by adding another dimension, such as the longitudinal relaxation time T 1 And diffusion coefficient D or magnetic field intensity G, and the like, thereby obtaining a nuclear magnetic resonance two-dimensional spectrum containing richer formation information. With the continuous development of complex reservoirs such as shale oil, compact oil and the like, the application of two-dimensional nuclear magnetic resonance logging in the field of oil exploration is more and more extensive, and the relaxation characteristics of different fluids are different, so that the qualitative identification and saturation calculation of the fluids can be carried out through the relaxation distribution characteristics in a two-dimensional spectrum.
At present, fluid identification in a two-dimensional spectrum mainly includes manually dividing relaxation distribution ranges of different fluid components by referring to a fluid identification plate or according to personal experience, and performing fluid identification and saturation calculation of a nuclear magnetic resonance logging two-dimensional spectrum of the same batch.
However, in reservoirs with complex fluid properties, the nuclear magnetic resonance two-dimensional spectrum fluid relaxation characteristics are also complex, and the fluid component identification error condition exists only by manually dividing each fluid component, so that the accuracy of the nuclear magnetic resonance logging fluid identification is influenced. Meanwhile, a large amount of nuclear magnetic resonance two-dimensional spectrum data are acquired by nuclear magnetic resonance logging, the efficiency of identifying fluid components in the two-dimensional spectrum manually is low, and the precision is low due to insufficient related interpretation experience of logging interpreters.
Disclosure of Invention
The invention aims to overcome the technical defects, provides a method, a device, electronic equipment and a storage medium for identifying components based on a nuclear magnetic resonance two-dimensional spectrum, and solves the technical problems of low component identification precision and low efficiency caused by identifying the components by relying on a fluid identification plate or experience of workers in the prior art.
In order to achieve the technical purpose, the invention adopts the following technical scheme:
in a first aspect, the invention provides a two-dimensional spectrum identification method based on nuclear magnetic resonance, which comprises the following steps:
acquiring peak positions of all fluid components in a relaxation total signal of a nuclear magnetic resonance two-dimensional spectrum;
determining the boundary position of the relaxation distribution range of each fluid component by taking each peak position as a central point, and determining the relaxation signal of each fluid component according to the peak position and the boundary position;
determining the relaxation signal intensity and the relaxation total signal intensity of each fluid component based on a preset integration method;
determining the saturation of each fluid component based on the ratio of the relaxation signal intensity of each fluid component to the intensity of the relaxation total signal.
In some implementations, the obtaining the peak positions of the respective fluid components in the relaxation total signal of the nuclear magnetic resonance two-dimensional spectrum includes:
marking amplitude points in the nuclear magnetic resonance two-dimensional spectrum;
scanning the amplitude points to obtain a first amplitude point, wherein the first amplitude point is larger than or equal to the value of any adjacent amplitude point;
determining a position of the first amplitude point in a two-dimensional spectrum;
and determining the peak position according to the amplitude difference relationship between the amplitude point of the two-dimensional spectrum and the marked amplitude point.
In some implementations, the determining the boundary position of the relaxation distribution range of each fluid component by taking the position of each peak as a central point respectively comprises:
inputting a nuclear magnetic resonance two-dimensional spectrum relaxation signal into a preset Roberts operator to obtain two-dimensional spectrum data only containing boundary points and marking the boundary points;
based on the boundary points, searching corresponding first boundary points along left and right adjacent point values of the dyeing point area by adopting a preset iteration method;
and smoothly connecting the first boundary points in a serial connection mode to obtain the boundary position of the relaxation distribution range of each fluid component.
In some implementations, the determining relaxation signals for individual fluid components based on the peak locations and boundary locations includes:
and resolving the relaxation signals of the fluid components from the total relaxation signal according to the boundary position and the peak position.
In some implementations, the determining the relaxation signal intensities and the relaxation total signal intensities for the respective fluid components based on a predetermined integration method includes:
acquiring the sum of the amplitudes of relaxation signals by adopting a preset area integration method, and determining the strength of the relaxation signals;
and acquiring the amplitude sum of the relaxation total signals by adopting a preset area integration method, and determining the strength of the relaxation total signals.
In some implementations, the determining the saturation of each fluid component according to a preset proportional relationship based on the relaxation signal intensity and the relaxation total signal intensity of each fluid component includes:
determining a preset proportional relation between the strength of the relaxation signals and the strength of the relaxation total signals as follows:
wherein S is n Is the saturation of the nth fluid component, F is the intensity of the total relaxation signal of the two-dimensional spectrum, F n Is the relaxation signal intensity of the nth fluid component.
In some implementations, after determining the peak position according to the position difference relationship between the amplitude point of the two-dimensional spectrum and the marked amplitude point, the method further includes:
and eliminating false peak positions in the peak positions by adopting a preset peak searching method according to the peak positions.
In a second aspect, the present invention further provides a two-dimensional spectrum identification device based on nuclear magnetic resonance, including:
a peak position acquisition unit for acquiring peak positions of respective fluid components in the relaxation total signal of the nuclear magnetic resonance two-dimensional spectrum;
a relaxation signal acquisition unit, configured to determine a boundary position of a relaxation distribution range of each fluid component by using each peak position as a central point, and determine a relaxation signal of each fluid component according to the peak position and the boundary position;
a relaxation signal intensity acquisition unit for determining the intensity of the relaxation signal and the intensity of the total relaxation signal of each fluid component based on a preset integration method;
a saturation determining unit, configured to determine a saturation of each of the fluid components based on a ratio of the relaxation signal intensity of each of the fluid components to the intensity of the total relaxation signal.
In a third aspect, the present invention further provides an electronic device, including: a processor and a memory;
the memory has stored thereon a computer readable program executable by the processor;
the processor, when executing the computer readable program, implements the steps in the two-dimensional spectrum identification method based on nuclear magnetic resonance as described above.
In a fourth aspect, the present invention also provides a computer-readable storage medium storing one or more programs, which are executable by one or more processors to implement the steps of the method for two-dimensional spectrum identification based on nuclear magnetic resonance as described above.
Compared with the prior art, the nuclear magnetic resonance two-dimensional spectrum identification method, the nuclear magnetic resonance two-dimensional spectrum identification device, the electronic equipment and the storage medium provided by the invention have the advantages that the peak positions of all fluid components in the relaxation total signal of the nuclear magnetic resonance two-dimensional spectrum are firstly obtained; then, respectively taking the peak positions as central points, determining the boundary positions of relaxation distribution ranges of all the fluid components, and determining relaxation signals of all the fluid components according to the peak positions and the boundary positions; then acquiring the relaxation signal intensity and the relaxation total signal intensity of each fluid component based on a preset integration method according to the relaxation signal and the relaxation total signal of each fluid component; finally, determining the saturation degree of each fluid component based on the ratio relation of the relaxation signal intensity of each fluid component in the relaxation total signal intensity; the efficiency and accuracy of fluid component identification are improved.
Drawings
FIG. 1 is a flow chart of an embodiment of a two-dimensional spectrum identification method based on nuclear magnetic resonance provided by the present invention;
fig. 2 is a schematic flowchart of an embodiment of step S101 in the two-dimensional spectrum identification method based on nuclear magnetic resonance provided in the present invention;
fig. 3 is a schematic flowchart of an embodiment of step S102 in the method for two-dimensional spectrum identification based on nuclear magnetic resonance provided in the present invention;
fig. 4 is a schematic flowchart of an embodiment of step S103 in the method for two-dimensional spectrum identification based on nuclear magnetic resonance provided in the present invention;
FIG. 5 is a schematic diagram of an embodiment of a two-dimensional spectrum identification apparatus based on nuclear magnetic resonance provided by the present invention;
fig. 6 is a schematic operating environment diagram of an embodiment of an electronic device provided in the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
According to the two-dimensional nuclear magnetic resonance-based two-dimensional spectrum identification method, device and equipment or computer-readable storage medium, two-dimensional nuclear magnetic resonance logging greatly improves the accuracy of fluid identification and saturation calculation in a reservoir through two-dimensional information. Due to the different information of the characterized fluids, the represented two-dimensional nuclear magnetic logging method is also different. Such as characterization bankA diffusion-relaxation (T2, D) two-dimensional NMR logging method of the transverse relaxation time T2 and the diffusion coefficient of the reservoir fluid, a (T2, T1) two-dimensional NMR logging method for representing the transverse relaxation time T2 and the longitudinal relaxation time T1 of the reservoir fluid, and the like. Although reservoir information represented by different two-dimensional nuclear magnetic resonance logging methods is different, the nuclear magnetic resonance two-dimensional spectra representing different information have consistent data structures. The invention is described by taking a commonly used (T2, T1) two-dimensional nuclear magnetic resonance as an example, and the echo amplitude b obtained by the (T2, T1) two-dimensional nuclear magnetic resonance logging measurement is Is shown as
Wherein b is is As the waiting time is T W,s Amplitude of an ith echo train of the time echo train; f. of jr For a longitudinal relaxation time of T 1r And transverse relaxation time is T 2j Nuclear magnetic signal intensity of time; epsilon is Is noise; m is transverse relaxation time T 2j The number of (2); n is the longitudinal relaxation time T 1r The number of the cells. By comparing the echo amplitudes b is Performing multi-echo-string joint inversion to obtain nuclear magnetic resonance two-dimensional spectrum f jr 。
The method, apparatus, device or computer readable storage medium of the present invention may be integrated with the above system or may be relatively independent.
Fig. 1 is a flowchart of a two-dimensional spectrum identification method based on nuclear magnetic resonance according to an embodiment of the present invention, and please refer to fig. 1, the two-dimensional spectrum identification method based on nuclear magnetic resonance includes:
s101, acquiring peak positions of all fluid components in a relaxation total signal of a nuclear magnetic resonance two-dimensional spectrum;
s102, determining the boundary position of the relaxation distribution range of each fluid component by taking each peak value position as a central point, and determining the relaxation signal of each fluid component according to the peak value position and the boundary position;
s103, determining the relaxation signal intensity and the relaxation total signal intensity of each fluid component based on a preset integration method;
and S104, determining the saturation of each fluid component based on the ratio of the relaxation signal intensity of each fluid component in the relaxation total signal intensity.
In step S101, it is critical to scan the entire nmr two-dimensional spectral data to identify the center of gravity of each fluid component. The echo information obtained by the nuclear magnetic resonance logging measurement comprises echo information of various fluids in the reservoir pores, and the relaxation characteristics in the nuclear magnetic resonance two-dimensional spectrum obtained by inverting the echo information are the linear superposition results of the relaxation characteristics of the various fluids in the reservoir pores. Meanwhile, the nuclear magnetic relaxation characteristics of different substances are different due to different physical properties, so that the fluid composition in the reservoir can be judged through the relaxation characteristics in the two-dimensional spectrum. Identifying relaxation distribution peak points in a two-dimensional spectrum by using a morphological method, and taking the peak points as the gravity centers of two-dimensional nuclear magnetic resonance relaxation distributions of fluid components, wherein one peak point corresponds to one fluid component; wherein, the morphological method is to use a special structural element to extract the characteristic signal in the relaxation total signal of the nuclear magnetic resonance two-dimensional spectrum so as to further obtain the peak signal of each component.
In step S102, the nuclear magnetic relaxation distribution characteristics (relaxation signals) of the single fluid component are distributed in a two-dimensional gaussian normal distribution in the nuclear magnetic resonance two-dimensional spectrum. The relaxation distribution in the nuclear magnetic resonance two-dimensional spectrum is formed by linearly superposing a plurality of two-dimensional Gaussian distributions, so that the relaxation distribution boundary of the fluid components has certain step property, and after the boundary point of the fluid components is determined, the boundary point around the peak position of the fluid components is connected, so that the relaxation distribution range of the fluid components is obtained.
In step S103, the intensity of the relaxation signal and the intensity of the total relaxation signal are obtained by a preset integration method, specifically, the area formed by the relaxation signal and the coordinate axis of each fluid component and the area surrounded by the total relaxation signal and the coordinate axis are obtained.
In the embodiment, the relaxation signals of the respective fluid components are decomposed from the relaxation total signal by acquiring the peak and boundary positions of the relaxation signals of the respective fluid components, the decomposition is independent of the fluid identification layout, the decomposition efficiency is high, and the decomposition is independent of the experience of workers and the decomposition accuracy is high.
In some embodiments, referring to fig. 2, the obtaining the peak positions of the fluid components in the relaxation total signal of the nuclear magnetic resonance two-dimensional spectrum includes:
s201, marking amplitude points in a nuclear magnetic resonance two-dimensional spectrum;
s202, scanning the amplitude points to obtain a first amplitude point, wherein the first amplitude point is larger than or equal to the value of any adjacent amplitude point;
s203, determining the position of the first amplitude point in the two-dimensional spectrum;
s204, determining the peak position according to the amplitude difference relation between the amplitude point of the two-dimensional spectrum and the marked amplitude point.
In this embodiment, the NMR spectrum is essentially a function f describing the characteristic distribution of the relaxation of the fluid jr Of a two-dimensional matrix, distribution function f jr The nuclear magnetic signal intensity is discretely distributed in a coordinate system consisting of m multiplied by n grids in a logarithmic uniform scale mode, the numerical value of each data point represents the nuclear magnetic signal intensity, wherein a peak point (j, r) in a two-dimensional spectrum is defined as an element which is greater than or equal to the signal intensity in the surrounding grids in the two-dimensional spectrum grid, and the peak position is determined by adopting a morphological method, namely:
first, a labeled two-dimensional spectrum J is constructed jr ，J jr ＝f jr -1. To mark two-dimensional spectrum J jr Scanning is carried out, and a first amplitude point meeting the following conditions is reserved:
J(j,r)≥J(j,r+1)&&J(j,r)≥J(j,r-1)&&J(j,r)≥J(j+1,r)&&J(j,r)≥J(j-1,r)
then, a queue is initialized (only front-end delete, back-end insert is allowed), and if J (J, r) is not equal to 0 and the neighbor is equal to 0, then J (J, r) is added to the queue, storing the position of J (J, r) in the queue.
Finally, the labeled two-dimensional spectrum J is updated jr . And circulating the queue, wherein the circulated point is the head of the queue, and the point is removed from the queue when being circulated. If any adjacent point of J (J, r) corresponding to the head of the queue is 0 and the point isIf a point is not 0 in the two-dimensional spectrum f (J, r), then the amplitude value of f (J, r) at that point is assigned to J (J, r) while the neighboring point is added to the queue. Looping to the stop with the queue empty. And (3) obtaining a two-dimensional matrix only containing the peak value position by subtracting the two-dimensional spectrum from the marked two-dimensional spectrum, screening the peak value position after marking the peak value position, dividing effective signal intervals, namely the peak value in the area is an effective peak value point, wherein part of the peak values are possibly caused by instrument measurement errors and inversion algorithm errors and do not accord with actual fluid relaxation characteristic distribution, such as noise. If the two peak points are close to each other in the peak point clock generated by peak searching, comparing the signal amplitude values corresponding to the two points, and taking the peak value with larger amplitude as an effective peak point. And finally, taking the peak point position capable of representing the relaxation characteristics of the reservoir fluid components as the gravity center of the fluid relaxation distribution.
In some embodiments, referring to fig. 3, said determining the boundary position of the relaxation distribution range of each fluid component by taking each peak position as a central point respectively comprises:
s301, inputting a nuclear magnetic resonance two-dimensional spectrum relaxation signal into a preset Roberts operator to obtain two-dimensional spectrum data only containing boundary points, and marking the boundary points;
s302, based on the boundary points, searching corresponding first boundary points along left and right adjacent point values of a dyeing point area by a preset iteration method;
and S303, smoothly connecting the first boundary points in a serial connection mode to obtain the boundary position of each fluid component relaxation distribution range.
In step S301, it is necessary to determine the transverse relaxation time T of each fluid component 2 Central spread width of peak and longitudinal relaxation time T 1 The relaxation distribution range of each fluid component can be divided by the spread width of the peak, but because the fluid properties in an actual reservoir are complex, the two-dimensional relaxation distribution of the fluid components does not completely conform to normal distribution and is difficult to express through a physical formula, and the actual relaxation distribution characteristics of the fluid have certain step property at the boundary, the Roberts operator for image edge detection is adopted to detect the relaxation of the fluidA Yu distribution boundary; the Robert operator is an operator detected by using a local difference operator, and a boundary with step property is identified by local difference, which is defined as follows:
the Roberts operator is divided into horizontal and vertical directions as shown in the following equation:
in this embodiment, two-dimensional spectrum data f (j, r) is input, and two-dimensional spectrum data g (j, r) only including boundary data points is output after convolution and operation are performed by using a Roberts operator, that is, the boundary points mark special values, and other points are 0; then, the peak position marked in the nuclear magnetic resonance two-dimensional spectrum is used as an anchor point for determining the distribution range of the single-component fluid, namely, a certain peak point is used as the center, a dyeing area is continuously and iteratively expanded towards an adjacent point (different values are marked at different peak positions), and if the values of a point to be dyed, a left adjacent point and a right adjacent point of the point are different, the boundary corresponding to the peak point is marked; thus, the boundary points corresponding to the peak points are marked one by one. Smoothly connecting the various boundary points, i.e. the range of the relaxation distribution of the fluid components.
In some embodiments, determining the relaxation signals for each fluid component from the peak locations and the boundary locations comprises:
and resolving the relaxation signals of the fluid components from the total relaxation signal according to the boundary position and the peak position.
In this embodiment, since the relaxation signals of the respective components are normally distributed, after the boundary position and the peak position are determined, a unique relaxation signal, that is, the relaxation signals of the respective components can be determined.
In some embodiments, referring to fig. 4, the determining the relaxation signal intensities and the relaxation total signal intensities for the respective fluid components based on a predetermined integration method includes:
s401, acquiring the sum of the amplitudes of relaxation signals by adopting a preset area integration method, and determining the strength of the relaxation signals;
s402, acquiring the sum of the amplitudes of the total relaxation signals by adopting a preset area integration method, and determining the strength of the total relaxation signals.
In the present embodiment, the intensity of the relaxation signal and the intensity of the total relaxation signal respectively represent the sum of the amplitudes of the respective component relaxation signals and the sum of the amplitudes of the total relaxation signal; when the signal intensity in the fluid component relaxation distribution range is calculated by adopting an area integration method, firstly, the fluid component relaxation distribution range is regarded as a polygon, a ray is sent from a point to be scanned to any direction, whether the ray is in the polygon is judged according to the number of intersection points of the ray and the polygon, and if the number of the intersection points is an odd number, a grid point to be scanned is positioned in the fluid relaxation distribution range; if the number of the intersection points is even number or no intersection point, the grid point to be scanned is positioned outside the fluid relaxation distribution range, if the grid point is in the range, the data values corresponding to the grid point are accumulated, if the grid point is not in the range, skipping is carried out, and after scanning is finished, the fluid component relaxation signal intensity F can be obtained n (ii) a Then, the full spectrum of the two-dimensional spectrum is scanned again, and the total signal intensity F of the two-dimensional spectrum, namely the total accumulated value of the full spectrum signal intensity, is calculated.
In some embodiments, the determining the saturation of each fluid component based on a ratio of the relaxation signal intensity of each fluid component in the intensity of the relaxation total signal includes:
determining a preset proportional relation between the strength of the relaxation signal and the strength of the relaxation total signal as follows:
wherein S is n Is the saturation of the nth fluid component, F is the intensity of the total relaxation signal of the two-dimensional spectrum, F n Is the relaxation signal intensity of the nth fluid component.
In this embodiment, the saturation of each fluid component is obtained by calculating a ratio of the relaxation signal intensity of the fluid component to the total spectrum signal intensity, and determining the content of each fluid component according to the ratio of the fluid component to the total component.
Based on the above two-dimensional spectrum identification method based on nuclear magnetic resonance, an embodiment of the present invention further provides a two-dimensional spectrum identification apparatus 500 based on nuclear magnetic resonance, referring to fig. 5, where the two-dimensional spectrum identification apparatus 500 based on nuclear magnetic resonance includes a peak position obtaining unit 510, a relaxation signal obtaining unit 520, a relaxation signal strength obtaining unit 530, and a saturation determining unit 540.
A peak position obtaining unit 510, configured to obtain peak positions of respective fluid components in a relaxation total signal of a nuclear magnetic resonance two-dimensional spectrum based on a preset morphological method;
a relaxation signal acquiring unit 520, configured to determine a boundary position of a relaxation distribution range of each fluid component by taking each peak position as a central point, and determine a relaxation signal of each fluid component according to the peak position and the boundary position;
a relaxation signal intensity acquiring unit 530 configured to acquire a relaxation signal intensity of each fluid component and an intensity of a relaxation total signal based on a preset integration method according to the relaxation signal and the relaxation total signal of each fluid component;
and a saturation determining unit 540 configured to determine the saturation of each fluid component according to a preset proportional relationship based on the strength of the relaxation signal and the strength of the total relaxation signal of each fluid component.
As shown in fig. 6, based on the above two-dimensional spectrum identification method based on nuclear magnetic resonance, the present invention also provides an electronic device, which may be a mobile terminal, a desktop computer, a notebook, a palm computer, a server, or other computing devices. The electronic device includes a processor 610, a memory 620, and a display 630. Fig. 6 shows only some of the components of the electronic device, but it is to be understood that not all of the shown components are required to be implemented, and that more or fewer components may be implemented instead.
The storage 620 may in some embodiments be an internal storage unit of the electronic device, such as a hard disk or a memory of the electronic device. The memory 620 may also be an external storage device of the electronic device in other embodiments, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), etc. provided on the electronic device. Further, the memory 620 may also include both internal storage units of the electronic device and external storage devices. The memory 620 is used for storing application software installed in the electronic device and various data, such as program codes for installing the electronic device. The memory 620 may also be used to temporarily store data that has been output or is to be output. In an embodiment, the memory 620 stores a two-dimensional spectrum identification program 640 based on nuclear magnetic resonance, and the two-dimensional spectrum identification program 640 based on nuclear magnetic resonance can be executed by the processor 610, so as to implement the two-dimensional spectrum identification method based on nuclear magnetic resonance according to the embodiments of the present application.
The processor 610 may be a Central Processing Unit (CPU), microprocessor or other data Processing chip in some embodiments, and is configured to execute program codes stored in the memory 620 or process data, such as performing a two-dimensional spectrum identification method based on nuclear magnetic resonance, or the like.
The display 630 may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, an OLED (Organic Light-Emitting Diode) touch panel, or the like in some embodiments. The display 630 is used for displaying the information of the nuclear magnetic resonance-based two-dimensional spectrum identification device and displaying a visual user interface. The components 610-630 of the electronic device communicate with each other via a system bus.
Of course, it will be understood by those skilled in the art that all or part of the processes of the methods of the above embodiments may be implemented by a computer program instructing relevant hardware (such as a processor, a controller, etc.), and the program may be stored in a computer readable storage medium, and when executed, the program may include the processes of the above method embodiments. The storage medium may be a memory, a magnetic disk, an optical disk, etc.
The above-described embodiments of the present invention should not be construed as limiting the scope of the present invention. Any other corresponding changes and modifications made according to the technical idea of the present invention should be included in the protection scope of the claims of the present invention.
Claims (10)
1. A two-dimensional spectrum identification method based on nuclear magnetic resonance is characterized by comprising the following steps:
acquiring peak positions of all fluid components in a relaxation total signal of a nuclear magnetic resonance two-dimensional spectrum;
determining the boundary position of the relaxation distribution range of each fluid component by taking each peak position as a central point, and determining the relaxation signal of each fluid component according to the peak position and the boundary position;
determining the relaxation signal intensity and the relaxation total signal intensity of each fluid component based on a preset integration method;
determining the saturation of each fluid component based on the ratio of the relaxation signal intensity of each fluid component to the intensity of the relaxation total signal.
2. The method for identifying two-dimensional spectra based on nuclear magnetic resonance according to claim 1, wherein the step of obtaining the peak position of each fluid component in the relaxation total signal of the nuclear magnetic resonance two-dimensional spectra comprises the following steps:
marking amplitude points in the nuclear magnetic resonance two-dimensional spectrum;
scanning the amplitude points to obtain a first amplitude point, wherein the first amplitude point is larger than or equal to the value of any adjacent amplitude point;
determining a position of the first amplitude point in a two-dimensional spectrum;
and determining the peak position according to the amplitude difference relationship between the amplitude point of the two-dimensional spectrum and the marked amplitude point.
3. The method for identifying two-dimensional spectra based on nuclear magnetic resonance according to claim 2, wherein the determining the boundary position of the relaxation distribution range of each fluid component with the position of each peak as the center point comprises:
inputting a nuclear magnetic resonance two-dimensional spectrum relaxation signal into a preset Roberts operator to obtain two-dimensional spectrum data only containing boundary points and marking the boundary points;
based on the boundary points, searching corresponding first boundary points along left and right adjacent point values of a dyeing point area by adopting a preset iteration method;
and smoothly connecting the first boundary points in a serial mode to obtain the boundary position of the relaxation distribution range of each fluid component.
4. The two-dimensional spectrum identification method based on nuclear magnetic resonance according to claim 3, wherein the determining relaxation signals of each fluid component according to the peak position and the boundary position comprises:
and resolving the relaxation signals of the fluid components from the total relaxation signal according to the boundary position and the peak position.
5. The two-dimensional spectrum identification method based on nuclear magnetic resonance according to claim 1, wherein the determining of the relaxation signal intensity and the relaxation total signal intensity of each fluid component based on a preset integration method comprises:
acquiring the sum of the amplitudes of relaxation signals by adopting a preset area integration method, and determining the strength of the relaxation signals;
and acquiring the amplitude sum of the relaxation total signals by adopting a preset area integration method, and determining the strength of the relaxation total signals.
6. The nmr two-dimensional spectrum identification method according to claim 1, wherein determining the saturation of each fluid component based on the ratio of the relaxation signal intensity of each fluid component to the total relaxation signal intensity comprises:
determining a preset proportional relation between the strength of the relaxation signal and the strength of the relaxation total signal as follows:
wherein S is n Is the saturation of the nth fluid component, F is the intensity of the total relaxation signal of the two-dimensional spectrum, F n Is the relaxation signal intensity of the nth fluid component.
7. The method for recognizing a two-dimensional spectrum based on nuclear magnetic resonance according to claim 2, wherein after determining the peak position according to the position difference relationship between the amplitude point of the two-dimensional spectrum and the marked amplitude point, the method further comprises:
and eliminating false peak positions in the peak positions by adopting a preset peak searching method according to the peak positions.
8. A two-dimensional spectrum recognition device based on nuclear magnetic resonance is characterized by comprising:
a peak position acquisition unit for acquiring the peak position of each fluid component in the relaxation total signal of the nuclear magnetic resonance two-dimensional spectrum;
a relaxation signal acquisition unit, configured to determine a boundary position of a relaxation distribution range of each fluid component by using each peak position as a central point, and determine a relaxation signal of each fluid component according to the peak position and the boundary position;
a relaxation signal intensity acquisition unit for determining the intensity of the relaxation signal and the intensity of the total relaxation signal of each fluid component based on a preset integration method;
a saturation determination unit for determining the saturation of each fluid component based on the ratio of the relaxation signal intensity of each fluid component to the intensity of the relaxation total signal.
9. An electronic device, comprising: a processor and a memory;
the memory has stored thereon a computer readable program executable by the processor;
the processor, when executing the computer readable program, implements the steps of the method for two-dimensional spectrum identification based on nuclear magnetic resonance according to any one of claims 1-7.
10. A computer-readable storage medium storing one or more programs, the one or more programs being executable by one or more processors to perform the steps of the method for two-dimensional nuclear magnetic resonance spectroscopy identification according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210581764.9A CN115201925B (en) | 2022-05-26 | 2022-05-26 | Nuclear magnetic resonance-based two-dimensional spectrum identification method, device, equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210581764.9A CN115201925B (en) | 2022-05-26 | 2022-05-26 | Nuclear magnetic resonance-based two-dimensional spectrum identification method, device, equipment and medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115201925A true CN115201925A (en) | 2022-10-18 |
| CN115201925B CN115201925B (en) | 2024-08-20 |
Family
ID=83576862
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210581764.9A Active CN115201925B (en) | 2022-05-26 | 2022-05-26 | Nuclear magnetic resonance-based two-dimensional spectrum identification method, device, equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115201925B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118191955A (en) * | 2024-02-02 | 2024-06-14 | 中海石油(中国)有限公司湛江分公司 | Nuclear magnetic resonance T2 spectrum component decomposition method, device, equipment and medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20210132250A1 (en) * | 2019-11-04 | 2021-05-06 | Exxonmobil Upstream Research Company | Estimating Pore and Fluid Characteristic Properties in Rock Samples Using Nuclear Magnetic Resonance Analyses |
| CN112799141A (en) * | 2021-01-05 | 2021-05-14 | 中国石油天然气集团有限公司 | Rapid two-dimensional nuclear magnetic resonance logging signal processing and T1T2 spectrum inversion method |
-
2022
- 2022-05-26 CN CN202210581764.9A patent/CN115201925B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20210132250A1 (en) * | 2019-11-04 | 2021-05-06 | Exxonmobil Upstream Research Company | Estimating Pore and Fluid Characteristic Properties in Rock Samples Using Nuclear Magnetic Resonance Analyses |
| CN112799141A (en) * | 2021-01-05 | 2021-05-14 | 中国石油天然气集团有限公司 | Rapid two-dimensional nuclear magnetic resonance logging signal processing and T1T2 spectrum inversion method |
Non-Patent Citations (1)
| Title |
|---|
| 李鹏举;张智鹏;姜大鹏;: "核磁共振测井流体识别方法综述", 测井技术, no. 05, 20 October 2011 (2011-10-20) * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118191955A (en) * | 2024-02-02 | 2024-06-14 | 中海石油(中国)有限公司湛江分公司 | Nuclear magnetic resonance T2 spectrum component decomposition method, device, equipment and medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115201925B (en) | 2024-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Giagopoulos et al. | Structural health monitoring and fatigue damage estimation using vibration measurements and finite element model updating | |
| RU2391686C2 (en) | Computer-based formation and inspection of training modes designed for multipoint geostatistical analysis | |
| Huseby et al. | Alternative environmental contours for structural reliability analysis | |
| Men’Shchikov et al. | A multi-scale, multi-wavelength source extraction method: getsources | |
| CN110794458B (en) | Gas content detection method and device based on time-frequency analysis and storage medium | |
| Adler et al. | Estimating thresholding levels for random fields via Euler characteristics | |
| BR112019011281A2 (en) | dive wave lighting using migrated collections | |
| CN116360000B (en) | Aeromagnetic detection method for crust material structure | |
| CN115201925B (en) | Nuclear magnetic resonance-based two-dimensional spectrum identification method, device, equipment and medium | |
| US9581707B2 (en) | System for automated identification of surfaces for building of geologic hydrodynamic model of oil and gas deposit by seismic data | |
| Kiar et al. | Comparing perturbation models for evaluating stability of neuroimaging pipelines | |
| CN113591606A (en) | Method and device for identifying hidden diseases of asphalt pavement, electronic equipment and storage medium | |
| CN109100803B (en) | Method and device for determining micro-fracture | |
| Lu et al. | Uncertainty of 1‐D fracture intensity measurements | |
| Lavrentiadis et al. | Generation of surface‐slip profiles in the wavenumber domain | |
| CN113818867B (en) | Method, system, medium, equipment and application for constructing pseudo capillary pressure curve | |
| RU2591239C1 (en) | Systems and methods for assessment of capability of collector system | |
| Spöck et al. | Incorporating covariance estimation uncertainty in spatial sampling design for prediction with trans-Gaussian random fields | |
| Schlegel et al. | Towards visual debugging for multi-target time series classification | |
| RU2530324C2 (en) | Method for determining position of marker depth coordinates when constructing geological model of deposit | |
| Akerjord et al. | Assessing mixing models within a common framework | |
| CN117910883B (en) | Logging data quality detection method, device, computing equipment and storage medium | |
| CN117951476A (en) | Construction method of lithology recognition model of shale oil reservoir and lithology recognition method | |
| CN116072232B (en) | Method, device, equipment and storage medium for determining relative permeability curve | |
| CN117310452B (en) | Method, device, computer equipment and storage medium for determining electromagnetic signal leakage |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |