A shale nuclear magnetic resonance aperture calibration method, device, equipment and medium

Abstract

The present application discloses a shale nuclear magnetic resonance pore size calibration method, device, equipment, and medium, relating to the technical field of shale reservoir pore characterization. The method comprises the following steps: preparing a shale core sample into a nuclear magnetic resonance column sample, drying the nuclear magnetic column sample under preset high temperature conditions, measuring the original nuclear magnetic resonance spectrum of the processed column sample, saturating the processed column sample with water under preset high pressure conditions, and measuring the target nuclear magnetic resonance spectrum of the saturated column sample; obtaining a shale pore nuclear magnetic resonance curve, drying the saturated column sample under preset high temperature conditions, performing a nitrogen adsorption experiment on the obtained powder sample, and calculating a nitrogen adsorption pore size distribution curve of the powder sample; obtaining a nuclear magnetic resonance porosity component curve and a nitrogen adsorption porosity component curve; performing nonlinear fitting on the curves to determine the corresponding relationship between nuclear magnetic resonance relaxation time and nitrogen pore size, and calibrating the shale nuclear magnetic resonance pore size based on the corresponding relationship. This method achieves pore size calibration.

Description

Technical Field

[0001] The present invention relates to the technical field of shale reservoir pore characterization, and in particular to a shale nuclear magnetic resonance pore size calibration method, device, equipment and medium. Background Art

[0002] Shale reservoir evaluation is very important for shale gas exploration and development. Since the pore size distribution range of shale reservoirs spans a wide range, with pores ranging from nanometer to micrometer scale, it is very difficult to achieve full pore size characterization of shale reservoirs. Therefore, the detailed characterization of shale pore structure has become the key to shale gas exploration and development.

[0003] The commonly used technologies for shale reservoir pore characterization include field emission scanning electron microscopy, gas adsorption, high-pressure mercury injection, nuclear magnetic resonance, etc. Currently, these experimental technologies all have problems. The pores characterized are mostly concentrated in a certain scale, and it is difficult to achieve full-pore shale pore characterization. Nuclear magnetic resonance is an emerging non-destructive experimental technology in recent years. It mainly uses the behavior of water molecules in a magnetic field to infer the pore structure of shale. When liquid molecules are in the pores, their magnetic resonance signals will be affected by the pore morphology, size and surface properties. Compared with other experimental techniques, its advantages are mainly reflected in the wide range of pores characterized and the experiment is non-destructive to the sample. However, nuclear magnetic resonance experiments are conducted through The relationship between relaxation time and pore size can be used to characterize the distribution characteristics of shale pores. The relaxation time is calibrated to the pore size, thus achieving the quantitative characterization of the multi-scale pore structure of shale reservoirs.

[0004] At present, NMR is usually combined with other technologies (such as shale column saturation-centrifugation experiment technology) to achieve NMR analysis of shale reservoirs. The calibration of relaxation time and pore size is mainly calculated by physical principle formula or nuclear magnetic resonance The curve shape is matched with other curves to achieve this, among which the method of calculating the physical principle formula is mainly through nuclear magnetic resonance The empirical formula for spectrum and aperture is calculated in conjunction with the formula for pressure and aperture. Due to the errors in the empirical formula, this method often results in large errors and inaccurate calibration results. However, due to the differences in measurement ranges between the various methods, peak values ​​can vary significantly. Direct curve matching results in highly subjective calibration results and is not applicable to a wide range of conditions. Summary of the Invention

[0005] In view of this, the purpose of the present invention is to provide a shale nuclear magnetic resonance pore size calibration method, device, equipment and medium, which can solve the problem of directly calibrating nuclear magnetic resonance based on nitrogen adsorption curve. The problem of large aperture error of spectrum curve is solved by improving calibration accuracy. The specific solution is as follows:

[0006] In a first aspect, the present application discloses a shale nuclear magnetic resonance pore size calibration method, comprising:

[0007] The shale core sample is prepared into a nuclear magnetic column sample, the nuclear magnetic column sample is dried under a preset high temperature condition, and the original nuclear magnetic resonance of the corresponding column sample after treatment is measured. The treated column sample is saturated with water under a preset high pressure condition, and the target nuclear magnetic resonance of the corresponding saturated column sample is measured. Atlas;

[0008] Based on the original NMR Spectra and the target NMR Shale pore nuclear magnetic resonance curve, drying the saturated column sample under a preset high temperature condition to obtain a powder sample, and performing a nitrogen adsorption experiment on the powder sample to calculate the nitrogen adsorption pore size distribution curve of the powder sample;

[0009] Based on the shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve; the nuclear magnetic resonance porosity component curve is a curve determined based on the nuclear magnetic resonance relaxation time and the nuclear magnetic resonance porosity component; the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size and the nitrogen adsorption porosity component;

[0010] The nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve are nonlinearly fitted to determine the corresponding relationship between the nuclear magnetic resonance relaxation time and the nitrogen pore size, and the nuclear magnetic resonance pore size of the shale is calibrated according to the corresponding relationship.

[0011] Optionally, the original nuclear magnetic resonance Spectra and the target NMR Shale pore nuclear magnetic resonance Curves, including:

[0012] The original NMR Spectrum and the target NMR The difference between the spectra is determined to be the shale pore NMR curve.

[0013] Optionally, performing a nitrogen adsorption experiment on the powder sample to calculate a nitrogen adsorption pore size distribution curve of the powder sample includes:

[0014] A nitrogen adsorption experiment is performed on the powder sample, and a nitrogen adsorption pore size distribution curve of the powder sample is obtained by BET specific surface area detection method and BJH method.

[0015] Optionally, the shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve, including:

[0016] NMR of the shale pores The values ​​corresponding to the curve are accumulated to obtain a first overall accumulated value;

[0017] Determine the shale pore NMR a first accumulated value of each data point corresponding to the curve, and dividing the first accumulated value by the first overall accumulated value to obtain an NMR porosity component of each data point;

[0018] determining the NMR porosity component curve according to each NMR porosity component and the NMR relaxation time;

[0019] Accumulating the values ​​corresponding to the nitrogen adsorption pore size distribution curve to obtain a second overall accumulated value;

[0020] determining a second cumulative value for each data point corresponding to the nitrogen adsorption pore size distribution curve, and dividing the second cumulative value by the second overall cumulative value to obtain a nitrogen adsorption porosity component for each data point;

[0021] The nitrogen adsorption porosity component curve is determined according to each of the nitrogen adsorption porosity components and the nitrogen pore diameter.

[0022] Optionally, before performing nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve, the method further includes:

[0023] The nitrogen adsorption data points are internally interpolated according to the number of NMR data points so that the shale pore NMR The data points of the curve correspond one-to-one with the data points of the nitrogen adsorption pore size distribution curve.

[0024] Optionally, before performing nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve, the method further includes:

[0025] The conversion coefficient of the nonlinear fitting is determined based on the plurality of the nuclear magnetic resonance porosity components and the plurality of the nitrogen adsorption porosity components by using a preset conversion coefficient determination formula; the preset conversion coefficient determination formula is:

[0026] ;

[0027] in, is the NMR porosity component; is the nitrogen adsorption porosity component; a and b are the conversion coefficients.

[0028] Optionally, performing nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve includes:

[0029] The nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve are nonlinearly fitted based on the conversion coefficient.

[0030] In a second aspect, the present application discloses a shale nuclear magnetic resonance aperture calibration device, comprising:

[0031] The spectrum measurement module is used to prepare the shale core sample into a nuclear magnetic column sample, dry the nuclear magnetic column sample under a preset high temperature condition, and measure the original nuclear magnetic resonance of the corresponding column sample after treatment. The treated column sample is saturated with water under a preset high pressure condition, and the target nuclear magnetic resonance of the corresponding saturated column sample is measured. Atlas;

[0032] Curve calculation module for calculating the original NMR Spectra and the target NMR Shale pore nuclear magnetic resonance curve, drying the saturated column sample under a preset high temperature condition to obtain a powder sample, and performing a nitrogen adsorption experiment on the powder sample to calculate the nitrogen adsorption pore size distribution curve of the powder sample;

[0033] Curve acquisition module for shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve; the nuclear magnetic resonance porosity component curve is a curve determined based on the nuclear magnetic resonance relaxation time and the nuclear magnetic resonance porosity component; the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size and the nitrogen adsorption porosity component;

[0034] The pore size calibration module is used to perform nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve to determine the corresponding relationship between the nuclear magnetic resonance relaxation time and the nitrogen pore size, and calibrate the nuclear magnetic resonance pore size of the shale according to the corresponding relationship.

[0035] In a third aspect, the present application discloses an electronic device, comprising:

[0036] Memory, used to store computer programs;

[0037] A processor is used to execute the computer program to implement the aforementioned shale nuclear magnetic resonance pore size calibration method.

[0038] In a fourth aspect, the present application discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, the aforementioned shale nuclear magnetic resonance aperture calibration method is implemented.

[0039] The present invention first prepares a shale core sample into a nuclear magnetic column sample, and then performs a drying treatment on the nuclear magnetic column sample under a preset high temperature condition, and measures the original nuclear magnetic resonance of the corresponding column sample after treatment. The treated column sample is saturated with water under a preset high pressure condition, and the target nuclear magnetic resonance of the corresponding saturated column sample is measured. Spectrum; Based on the original NMR Spectra and the target NMR Shale pore nuclear magnetic resonance curve, drying the saturated column sample under a preset high temperature condition to obtain a powder sample, performing a nitrogen adsorption experiment on the powder sample to calculate the nitrogen adsorption pore size distribution curve of the powder sample; then based on the shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve; the nuclear magnetic resonance porosity component curve is a curve determined based on the nuclear magnetic resonance relaxation time and the nuclear magnetic resonance porosity component; the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size and the nitrogen adsorption porosity component; finally, the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve are nonlinearly fitted to determine the correspondence between the nuclear magnetic resonance relaxation time and the nitrogen pore size, and the nuclear magnetic resonance pore size of the shale is calibrated according to the correspondence. It can be seen that the present application obtains the nuclear magnetic resonance porosity component curve and the nitrogen adsorption porosity component curve through multiple experiments, and determines the correspondence between the nuclear magnetic resonance relaxation time and the nitrogen pore size by fitting the two curves, thereby realizing the precise calibration of the nuclear magnetic resonance pore size by nitrogen adsorption and improving the accuracy of pore size calibration. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are merely embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0041] Figure 1 This is a flow chart of a shale nuclear magnetic resonance pore size calibration method disclosed in this application;

[0042] Figure 2 A sample nuclear magnetic resonance disclosed in this application Schematic diagram of the curve;

[0043] Figure 3 This is a schematic diagram of the nitrogen adsorption pore size distribution of a sample disclosed in this application;

[0044] Figure 4 This is a schematic diagram of nitrogen adsorption porosity component distribution disclosed in this application;

[0045] Figure 5 This is a schematic diagram of the distribution of nuclear magnetic porosity components disclosed in this application;

[0046] Figure 6 This is a schematic diagram of fitting of nuclear magnetic resonance and nitrogen adsorption porosity components disclosed in this application;

[0047] Figure 7 This is a schematic diagram of a nuclear magnetic resonance pore size distribution disclosed in this application;

[0048] Figure 8 This is a schematic structural diagram of a shale nuclear magnetic resonance aperture calibration device disclosed in this application;

[0049] Figure 9 This is a structural diagram of an electronic device disclosed in this application. DETAILED DESCRIPTION

[0050] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0051] At present, NMR is usually combined with other technologies (such as shale column saturation-centrifugation experiment technology) to achieve NMR analysis of shale reservoirs. The calibration of relaxation time and pore size is mainly calculated by physical principle formula or nuclear magnetic resonance The curve shape is matched with other curves to achieve this, among which the method of calculating the physical principle formula is mainly through nuclear magnetic resonance The empirical formula of spectrum and pore size is calculated together with the formula of pressure and pore size. Due to the error problem of the empirical formula, this method usually has a large error and the calibration result is inaccurate. However, the method of peak matching of curve morphology has a large difference in peak value due to the difference in measurement range of each method. Direct curve morphology matching will result in a highly subjective calibration result and is not applicable to a large range. In order to solve the above technical problems, the present application discloses a shale nuclear magnetic resonance pore size calibration method, device, equipment and medium, which can solve the problem of directly calibrating nuclear magnetic resonance based on nitrogen adsorption curve. The problem of large aperture error of spectrum curve is solved to improve calibration accuracy.

[0052] See also Figure 1 As shown, an embodiment of the present invention discloses a shale nuclear magnetic resonance aperture calibration method, comprising:

[0053] Step S11: Prepare the shale core sample into a nuclear magnetic column sample, dry the nuclear magnetic column sample under a preset high temperature condition, and measure the original nuclear magnetic resonance of the corresponding column sample after treatment. The treated column sample is saturated with water under a preset high pressure condition, and the target nuclear magnetic resonance of the corresponding saturated column sample is measured. Atlas.

[0054] In this embodiment, the present application first prepared fresh shale core samples, prepared them into standard nuclear magnetic column samples (cylinders with a length of 5 cm and a diameter of 2.5 cm), dried the column samples at high temperature (60°C) for 24 hours, and measured their original nuclear magnetic resonance (NMR) After that, the column sample was taken out for high-pressure saturation, that is, the column sample was saturated with water under the preset high-pressure conditions, saturated with formation water at 20 MPa pressure for 48 hours, and its nuclear magnetic resonance after saturation was measured. Spectrum, and then obtain the original NMR Spectra and targeted NMR of the column sample after saturation treatment Atlas.

[0055] Step S12: Based on the original nuclear magnetic resonance Spectra and the target NMR Shale pore nuclear magnetic resonance curve, drying the saturated water column sample under a preset high temperature condition to obtain a powder sample, and performing a nitrogen adsorption experiment on the powder sample to calculate the nitrogen adsorption pore size distribution curve of the powder sample.

[0056] In this embodiment, the original NMR Spectra and targeted NMR of the column sample after saturation treatment After the spectrum, the water-saturated NMR Spectra minus the original NMR shallow shale pore nuclear magnetic resonance Curve. The original NMR Atlas Targeted NMR Atlas The difference between the two values ​​is determined as the shale pore NMR curve .

[0057] ;

[0058] Then, the column sample after saturated NMR was dried at high temperature (60°) for 48 hours and processed into a powder sample with a mesh size of 50 to 10. The powder sample was subjected to a nitrogen adsorption experiment to calculate the nitrogen adsorption pore size distribution curve of the powder sample. Specifically, the nitrogen adsorption pore size distribution curve was calculated using the BET method (Formula 1) and the BJH method (Formula 2):

[0059] ; (Formula 1)

[0060] Where V is the nitrogen adsorption capacity (cm 3 / g); is the relative pressure, dimensionless; C is the BET constant, dimensionless; is the monolayer molecular adsorption capacity (cm 3 / g).

[0061] ; (Formula 2)

[0062] in, is the small change in adsorption amount related to pore size (cm 3 / g); The maximum adsorption capacity (cm 3 / g); is the density of nitrogen (cm 3 / g); is the experimental pressure (MPa); is the experimental pressure (MPa).

[0063] Step S13: Based on the shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve; the nuclear magnetic resonance porosity component curve is a curve determined based on the nuclear magnetic resonance relaxation time and the nuclear magnetic resonance porosity component; the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size and the nitrogen adsorption porosity component.

[0064] In this embodiment, shale pore nuclear magnetic resonance curve and nitrogen adsorption pore size distribution curve, based on shale pore nuclear magnetic resonance The corresponding NMR porosity component curve and nitrogen adsorption porosity component curve are obtained from the obtained shale pore NMR. The curve and the nitrogen adsorption pore size distribution curve are processed, and the two sets of curve values ​​are accumulated respectively to obtain their final overall cumulative value, and the accumulated value of each data point is divided by the final overall cumulative value to obtain its individual data porosity components. Among them, for the accumulated value of each data point, for example, if there are 10 data points, the accumulated value of the 3rd data point is the value corresponding to the third data point plus the values ​​corresponding to the first two data points; the accumulated value of the 4th data point is the value corresponding to the 4th data point plus the values ​​corresponding to the first 3 data points. The overall cumulative value is the accumulated value of 10 data points. It should be noted that the nuclear magnetic resonance porosity component curve here is a curve determined based on the nuclear magnetic resonance relaxation time (horizontal axis) and the nuclear magnetic resonance porosity component (vertical axis); the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size (horizontal axis) and the nitrogen adsorption porosity component (vertical axis). Specifically:

[0065] ;

[0066] in, is the NMR porosity component; For nuclear magnetic resonance Spectral signal intensity.

[0067] ;

[0068] in, is the nitrogen adsorption porosity component; is the nitrogen adsorption pore volume.

[0069] Step S14: performing nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve to determine the corresponding relationship between the nuclear magnetic resonance relaxation time and the nitrogen pore size, and calibrating the nuclear magnetic resonance pore size of the shale according to the corresponding relationship.

[0070] In this embodiment, due to the nuclear magnetic resonance The frequency of data acquisition of the curve is high, and the data points are more, while the data points of nitrogen adsorption are fewer. Therefore, the exponential fitting in the origin software is used to internally interpolate the data points of nitrogen adsorption between 2 and 50 nm (the data of the mesopore part). The number of curve data points is adjusted to match the nitrogen data points. The NMR full pore size data is then fitted with the mesopore data (N>30) to calibrate the NMR pore size distribution data. The conversion coefficient for the nonlinear fit is determined based on several NMR porosity components and several nitrogen adsorption porosity components using a preset conversion coefficient determination formula; the preset conversion coefficient determination formula is:

[0071] ;

[0072] in, is the NMR porosity component; is the nitrogen adsorption porosity component; a and b are conversion coefficients. Based on the conversion coefficients, the nitrogen adsorption porosity component curve and the NMR porosity component curve are nonlinearly fitted. Where:

[0073] ;

[0074] ;

[0075] is the NMR transverse relaxation time, d is the aperture; is the pore shape factor; is the surface area of ​​fluid-containing pores; wei surface relaxation rate; The fluid pore volume is used to clarify the correspondence between NMR relaxation time and nitrogen pore size. The NMR pore size of shale is then calibrated based on this correspondence, achieving the correspondence and conversion between NMR relaxation time and pore size.

[0076] In summary, the present invention first prepares the shale core sample into a nuclear magnetic column sample, and then dry the nuclear magnetic column sample under a preset high temperature condition, and measures the original nuclear magnetic resonance of the corresponding column sample after treatment. The treated column sample is saturated with water under a preset high pressure condition, and the target nuclear magnetic resonance of the corresponding saturated column sample is measured. Spectrum; Based on the original NMR Spectra and the target NMR Shale pore nuclear magnetic resonance curve, drying the saturated column sample under a preset high temperature condition to obtain a powder sample, performing a nitrogen adsorption experiment on the powder sample to calculate the nitrogen adsorption pore size distribution curve of the powder sample; then based on the shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve; the nuclear magnetic resonance porosity component curve is a curve determined based on the nuclear magnetic resonance relaxation time and the nuclear magnetic resonance porosity component; the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size and the nitrogen adsorption porosity component; finally, the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve are nonlinearly fitted to determine the correspondence between the nuclear magnetic resonance relaxation time and the nitrogen pore size, and the nuclear magnetic resonance pore size of the shale is calibrated according to the correspondence. It can be seen that the present application obtains the nuclear magnetic resonance porosity component curve and the nitrogen adsorption porosity component curve through multiple experiments, and determines the correspondence between the nuclear magnetic resonance relaxation time and the nitrogen pore size by fitting the two curves, thereby realizing the precise calibration of the nuclear magnetic resonance pore size by nitrogen adsorption and improving the accuracy of pore size calibration.

[0077] Next, the specific process of the shale NMR aperture calibration method will be described in detail with reference to specific embodiments.

[0078] In a specific embodiment, for the calibration of pore diameter in shallow shale (current burial depth less than 2000m), five fresh shale core samples were first selected and the same sample was prepared into a standard NMR column sample with a length of 5cm and a diameter of 2.5cm. The column sample was then dried at high temperature (60°C) for 24 hours and its original NMR was measured. The column sample was then taken out for high-pressure saturation. The formation water was saturated at 20 MPa for 48 hours, and the nuclear magnetic resonance (NMR) after saturation was measured. Spectrum, and then obtain Figure 2 NMR images of the five samples shown Spectrum. Using NMR after saturation of water Spectra minus the original NMR shallow shale pore nuclear magnetic resonance curve.

[0079] Afterwards, the column sample after saturated NMR was dried at high temperature (60°) for 48 hours and processed into 50-10 mesh powder samples for nitrogen adsorption experiment. The BET method and BJH method were used to calculate the following: Figure 3 Its nitrogen adsorption pore size distribution curve is shown.

[0080] The obtained NMR and nitrogen adsorption data were processed according to the following formula. The curve and nitrogen adsorption pore size distribution data were processed as follows Figure 4 and Figure 5 Porosity component data shown.

[0081] ;

[0082] in, is the NMR porosity component; For nuclear magnetic resonance Spectral signal intensity.

[0083] ;

[0084] in, is the nitrogen adsorption porosity component; is the nitrogen adsorption pore volume.

[0085] However, since the nitrogen adsorption data is mainly for the pores in the mesopore part, the pore size distribution range of the control is mainly in the range of 2~50nm, and since the pore size distribution range of the nuclear magnetic resonance curve is relatively large, linear fitting is not applicable, so nonlinear fitting is used. The data in the nitrogen pore size range of 2~50nm are intercepted respectively, and the porosity component corresponding to the nuclear magnetic resonance porosity component is selected according to the porosity component data of nitrogen adsorption. According to the number of NMR data points, the data points of nitrogen adsorption are internally interpolated to make them correspond. Then, the porosity components of the two methods are used as the basis to clarify the correspondence between NMR relaxation time and nitrogen pore size, such as Figure 6 As shown. The conversion coefficient a is determined to be 50.537 and b is 0.9106. Finally, the following is obtained: Figure 7 The calibrated NMR pore size distribution diagram is shown.

[0086] It can be seen that this application obtains the NMR porosity component curve and the nitrogen adsorption porosity component curve through multiple experiments, and determines the correspondence between the NMR relaxation time and the nitrogen pore size by fitting the two curves, thereby realizing the precise calibration of the NMR pore size by nitrogen adsorption and improving the accuracy of the pore size calibration.

[0087] See also Figure 8 As shown, an embodiment of the present invention discloses a shale nuclear magnetic resonance aperture calibration device, comprising:

[0088] The spectrum measurement module 11 is used to prepare the shale core sample into a nuclear magnetic column sample, dry the nuclear magnetic column sample under a preset high temperature condition, and measure the original nuclear magnetic resonance of the corresponding column sample after treatment. The treated column sample is saturated with water under a preset high pressure condition, and the target nuclear magnetic resonance of the corresponding saturated column sample is measured. Atlas;

[0089] Curve calculation module 12, for calculating the original nuclear magnetic resonance Spectra and the target NMR Shale pore nuclear magnetic resonance curve, drying the saturated column sample under a preset high temperature condition to obtain a powder sample, and performing a nitrogen adsorption experiment on the powder sample to calculate the nitrogen adsorption pore size distribution curve of the powder sample;

[0090] Curve acquisition module 13, for obtaining the shale pore NMR curve based on the shale pore NMR curve. The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve; the nuclear magnetic resonance porosity component curve is a curve determined based on the nuclear magnetic resonance relaxation time and the nuclear magnetic resonance porosity component; the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size and the nitrogen adsorption porosity component;

[0091] The pore size calibration module 14 is used to perform nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve to determine the correspondence between the nuclear magnetic resonance relaxation time and the nitrogen pore size, and calibrate the nuclear magnetic resonance pore size of the shale according to the correspondence.

[0092] In summary, the present invention first prepares the shale core sample into a nuclear magnetic column sample, and then dry the nuclear magnetic column sample under a preset high temperature condition, and measures the original nuclear magnetic resonance of the corresponding column sample after treatment. The treated column sample is saturated with water under a preset high pressure condition, and the target nuclear magnetic resonance of the corresponding saturated column sample is measured. Spectrum; Based on the original NMR Spectra and the target NMR Shale pore nuclear magnetic resonance curve, drying the saturated column sample under a preset high temperature condition to obtain a powder sample, performing a nitrogen adsorption experiment on the powder sample to calculate the nitrogen adsorption pore size distribution curve of the powder sample; then based on the shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve; the nuclear magnetic resonance porosity component curve is a curve determined based on the nuclear magnetic resonance relaxation time and the nuclear magnetic resonance porosity component; the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size and the nitrogen adsorption porosity component; finally, the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve are nonlinearly fitted to determine the correspondence between the nuclear magnetic resonance relaxation time and the nitrogen pore size, and the nuclear magnetic resonance pore size of the shale is calibrated according to the correspondence. It can be seen that the present application obtains the nuclear magnetic resonance porosity component curve and the nitrogen adsorption porosity component curve through multiple experiments, and determines the correspondence between the nuclear magnetic resonance relaxation time and the nitrogen pore size by fitting the two curves, thereby realizing the precise calibration of the nuclear magnetic resonance pore size by nitrogen adsorption and improving the accuracy of pore size calibration.

[0093] In some specific embodiments, the curve calculation module 12 may specifically include:

[0094] A curve determining unit is used to convert the original nuclear magnetic resonance Spectrum and the target NMR The difference between the spectra is determined to be the shale pore NMR curve.

[0095] In some specific embodiments, the curve calculation module 12 may specifically include:

[0096] The curve calculation unit is used to perform a nitrogen adsorption experiment on the powder sample and calculate the nitrogen adsorption pore size distribution curve of the powder sample by using the BET specific surface area detection method and the BJH method.

[0097] In some specific embodiments, the curve acquisition module 13 may specifically include:

[0098] The first accumulating unit is used to measure the shale pore nuclear magnetic resonance The values ​​corresponding to the curve are accumulated to obtain a first overall accumulated value;

[0099] The NMR porosity component acquisition unit is used to determine the shale pore NMR a first accumulated value of each data point corresponding to the curve, and dividing the first accumulated value by the first overall accumulated value to obtain an NMR porosity component of each data point;

[0100] a first curve determining unit, configured to determine the NMR porosity component curve according to each NMR porosity component and the NMR relaxation time;

[0101] a second accumulating unit, configured to accumulate the values ​​corresponding to the nitrogen adsorption pore size distribution curve to obtain a second overall accumulated value;

[0102] a nitrogen adsorption porosity component acquisition unit, configured to determine a second cumulative value of each data point corresponding to the nitrogen adsorption pore size distribution curve, and divide the second cumulative value by the second overall cumulative value to acquire a nitrogen adsorption porosity component of each data point;

[0103] The second curve determining unit is configured to determine the nitrogen adsorption porosity component curve according to each of the nitrogen adsorption porosity components and the nitrogen pore diameter.

[0104] In some specific embodiments, the device may further include:

[0105] The difference module is used to perform internal interpolation on the nitrogen adsorption data points according to the number of NMR data points so that the shale pore NMR The data points of the curve correspond one-to-one with the data points of the nitrogen adsorption pore size distribution curve.

[0106] In some specific embodiments, the device may further include:

[0107] A conversion coefficient determination module is used to determine the conversion coefficient of the nonlinear fitting based on the plurality of nuclear magnetic resonance porosity components and the plurality of nitrogen adsorption porosity components using a preset conversion coefficient determination formula; the preset conversion coefficient determination formula is:

[0108] ;

[0109] in, is the NMR porosity component; is the nitrogen adsorption porosity component; a and b are the conversion coefficients.

[0110] In some specific embodiments, the aperture calibration module 14 may specifically include:

[0111] A nonlinear fitting unit is used to perform nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve based on the conversion coefficient.

[0112] Furthermore, the embodiment of the present application also discloses an electronic device, Figure 9 This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content in the diagram should not be considered as any limitation to the scope of application of the present application.

[0113] Figure 9This is a schematic diagram of the structure of an electronic device 20 provided in an embodiment of the present application. The electronic device 20 may include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is used to store a computer program, which is loaded and executed by the processor 21 to implement the relevant steps of the shale NMR aperture calibration method disclosed in any of the aforementioned embodiments. Furthermore, the electronic device 20 in this embodiment may be a computer.

[0114] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and the external device. The communication protocol it follows is any communication protocol that can be applied to the technical solution of this application and is not specifically limited here; the input and output interface 25 is used to obtain external input data or output data to the outside world. Its specific interface type can be selected according to specific application needs and is not specifically limited here.

[0115] In addition, the memory 22 as a carrier for resource storage can be a read-only memory, random access memory, disk or optical disk, etc. The resources stored thereon can include an operating system 221, a computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0116] The operating system 221 is used to manage and control the hardware devices and computer program 222 on the electronic device 20. The operating system 221 can be Windows Server, NetWare, Unix, Linux, etc. In addition to including a computer program capable of performing the shale nuclear magnetic resonance aperture calibration method performed by the electronic device 20 as disclosed in any of the aforementioned embodiments, the computer program 222 can further include computer programs capable of performing other specific tasks.

[0117] Furthermore, this application discloses a computer-readable storage medium for storing a computer program. When executed by a processor, the computer program implements the aforementioned shale NMR aperture calibration method. The specific steps of this method can be found in the aforementioned embodiments and will not be further elaborated here.

[0118] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from the other embodiments. Reference can be made to the descriptions of the identical or similar parts between the various embodiments. For the devices disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple, and the relevant parts can be referred to the descriptions of the methods.

[0119] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the components and steps of each example according to their functions. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians may use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0120] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in random access memory (RAM), internal memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0121] Finally, it should be noted that, in this document, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or device. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article, or device comprising the element.

[0122] The above is a detailed introduction to the technical solution provided by the present application. Specific examples are used herein to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method of the present application and its core idea. At the same time, for those skilled in the art, according to the ideas of the present application, there may be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present application.

Claims

1. A shale nuclear magnetic resonance aperture calibration method, characterized in that: include: The shale core sample is prepared into a nuclear magnetic column sample, the nuclear magnetic column sample is dried under a preset high temperature condition, and the original nuclear magnetic resonance of the corresponding column sample after treatment is measured. The treated column sample is saturated with water under a preset high pressure condition, and the target nuclear magnetic resonance of the corresponding saturated column sample is measured. Atlas; Based on the original NMR Spectra and the target NMR Shale pore nuclear magnetic resonance curve, drying the saturated column sample under a preset high temperature condition to obtain a powder sample, and performing a nitrogen adsorption experiment on the powder sample to calculate the nitrogen adsorption pore size distribution curve of the powder sample; Based on the shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve; the nuclear magnetic resonance porosity component curve is a curve determined based on the nuclear magnetic resonance relaxation time and the nuclear magnetic resonance porosity component; the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size and the nitrogen adsorption porosity component; The nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve are nonlinearly fitted to determine the corresponding relationship between the nuclear magnetic resonance relaxation time and the nitrogen pore size, and the nuclear magnetic resonance pore size of the shale is calibrated according to the corresponding relationship.

2. The shale nuclear magnetic resonance aperture calibration method according to claim 1, characterized in that: Based on the original nuclear magnetic resonance Spectra and the target NMR Shale pore nuclear magnetic resonance Curves, including: The original NMR Spectrum and the target NMR The difference between the spectra is determined to be the shale pore NMR curve.

3. The shale nuclear magnetic resonance aperture calibration method according to claim 1, characterized in that: The performing of a nitrogen adsorption experiment on the powder sample to calculate a nitrogen adsorption pore size distribution curve of the powder sample comprises: A nitrogen adsorption experiment is performed on the powder sample, and a nitrogen adsorption pore size distribution curve of the powder sample is obtained by BET specific surface area detection method and BJH method.

4. The shale nuclear magnetic resonance aperture calibration method according to claim 1, characterized in that: The shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve, including: NMR of the shale pores The values ​​corresponding to the curve are accumulated to obtain a first overall accumulated value; Determine the shale pore NMR a first accumulated value of each data point corresponding to the curve, and dividing the first accumulated value by the first overall accumulated value to obtain an NMR porosity component of each data point; determining the NMR porosity component curve according to each NMR porosity component and the NMR relaxation time; Accumulating the values ​​corresponding to the nitrogen adsorption pore size distribution curve to obtain a second overall accumulated value; determining a second cumulative value for each data point corresponding to the nitrogen adsorption pore size distribution curve, and dividing the second cumulative value by the second overall cumulative value to obtain a nitrogen adsorption porosity component for each data point; The nitrogen adsorption porosity component curve is determined according to each of the nitrogen adsorption porosity components and the nitrogen pore diameter.

5. The shale nuclear magnetic resonance aperture calibration method according to claim 1, characterized in that: Before performing nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve, the method further includes: The nitrogen adsorption data points are internally interpolated according to the number of NMR data points so that the shale pore NMR The data points of the curve correspond one-to-one with the data points of the nitrogen adsorption pore size distribution curve.

6. The shale nuclear magnetic resonance aperture calibration method according to any one of claims 1 to 5, characterized in that: Before performing nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve, the method further includes: The conversion coefficient of the nonlinear fitting is determined based on the plurality of the nuclear magnetic resonance porosity components and the plurality of the nitrogen adsorption porosity components by using a preset conversion coefficient determination formula; the preset conversion coefficient determination formula is: ； in, is the NMR porosity component; is the nitrogen adsorption porosity component; a and b are the conversion coefficients.

7. The shale nuclear magnetic resonance aperture calibration method according to claim 6, characterized in that: The performing nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve comprises: The nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve are nonlinearly fitted based on the conversion coefficient.

8. A shale nuclear magnetic resonance aperture calibration device, characterized in that: include: The spectrum measurement module is used to prepare the shale core sample into a nuclear magnetic column sample, dry the nuclear magnetic column sample under a preset high temperature condition, and measure the original nuclear magnetic resonance of the corresponding column sample after treatment. The treated column sample is saturated with water under a preset high pressure condition, and the target nuclear magnetic resonance of the corresponding saturated column sample is measured. Atlas; Curve calculation module for calculating the original NMR Spectra and the target NMR Shale pore nuclear magnetic resonance curve, drying the saturated column sample under a preset high temperature condition to obtain a powder sample, and performing a nitrogen adsorption experiment on the powder sample to calculate the nitrogen adsorption pore size distribution curve of the powder sample; Curve acquisition module for shale pore nuclear magnetic resonance The curve and the nitrogen adsorption pore size distribution curve respectively obtain the corresponding nuclear magnetic resonance porosity component curve and nitrogen adsorption porosity component curve; the nuclear magnetic resonance porosity component curve is a curve determined based on the nuclear magnetic resonance relaxation time and the nuclear magnetic resonance porosity component; the nitrogen adsorption porosity component curve is a curve determined based on the nitrogen pore size and the nitrogen adsorption porosity component; The pore size calibration module is used to perform nonlinear fitting on the nitrogen adsorption porosity component curve and the nuclear magnetic resonance porosity component curve to determine the corresponding relationship between the nuclear magnetic resonance relaxation time and the nitrogen pore size, and calibrate the nuclear magnetic resonance pore size of the shale according to the corresponding relationship.

9. An electronic device, characterized in that: include: Memory, used to store computer programs; A processor, configured to execute the computer program to implement the shale nuclear magnetic resonance aperture calibration method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that Used to store a computer program; wherein, when the computer program is executed by a processor, the shale nuclear magnetic resonance aperture calibration method according to any one of claims 1 to 7 is implemented.

Images