Full-size aperture characterization method and device
By employing Gaussian function decomposition and pore size transformation, the problem of low accuracy in NMR T2 spectroscopy was solved, enabling full-size pore size characterization and improving the accuracy of pore size spectroscopy and evaluation of pore structure in strongly heterogeneous reservoirs.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-09-03
- Publication Date
- 2026-03-10
AI Technical Summary
In existing technologies, the pore size spectrum obtained after pore size conversion of the nuclear magnetic resonance T2 spectrum has low accuracy, which cannot meet the needs of practical applications and cannot accurately evaluate the pore structure of strongly heterogeneous reservoirs.
The core NMR T2 spectrum was decomposed using Gaussian function decomposition. The number and position parameters of the Gaussian function were determined by the multi-derivative peak finding method. The spectrum was decomposed into multiple sets of component NMR T2 spectra with different pore size ranges. After pore size conversion, the full-size pore size spectrum was finally formed.
It improves the accuracy of pore size spectrum, realizes full-size pore size characterization, and can more accurately evaluate the pore structure of strongly heterogeneous reservoirs.
Smart Images
Figure CN121633166A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of petroleum geology, and particularly relates to a full-size pore size characterization method and device. BACKGROUND
[0002] This section is intended to provide background or context to the embodiments of the application recited in the claims. The description herein does not constitute admission that the prior art is prior art nor does it constitute an admission of any description in this section as prior art to an application described herein and / or in another application also owned by the applicant of the present application.
[0003] Strongly heterogeneous reservoirs have the characteristics of multiple types of reservoir spaces and complex pore structures, which contain clay minerals and intergranular pores, intragranular pores and organic matter pores, and the pore diameters of different reservoir spaces differ greatly. Therefore, accurate evaluation of the pore structure of strongly heterogeneous reservoirs is of great significance to fill point optimization and productivity evaluation. There are many experimental methods for core pore structure evaluation, including scanning electron microscopy, casting thin section, nitrogen adsorption, high-pressure mercury injection, nuclear magnetic resonance, etc. The prior art also proposes to combine multiple pore size characterization methods, but the pore size spectrum obtained by pore size conversion of the nuclear magnetic T2 spectrum has low precision and the pore size is not comprehensive, which cannot meet the needs of practical applications. Therefore, there is an urgent need for a method to solve the above problems. SUMMARY
[0004] The embodiments of the present application provide a full-size pore size characterization method to improve the precision of the size pore size spectrum and realize pore size characterization under full size. The method comprises:
[0005] Obtaining a core nuclear magnetic T2 spectrum obtained by core nuclear magnetic resonance experiment;
[0006] Decomposing the core nuclear magnetic T2 spectrum by using multiple Gaussian functions to obtain multiple groups of component nuclear magnetic T2 spectra of different pore size intervals; wherein the process of decomposing the core nuclear magnetic T2 spectrum into multiple groups of component nuclear magnetic T2 spectra of different pore size intervals by using Gaussian functions is as follows: determining multiple Gaussian functions, the total number and position parameters of the multiple Gaussian functions by using a multiple-order derivative peak finding method; characterizing the core nuclear magnetic T2 spectrum according to the multiple Gaussian functions, the total number and the position parameters of the multiple Gaussian functions; decomposing the characterized core nuclear magnetic T2 spectrum to obtain multiple component nuclear magnetic T2 spectra; dividing the multiple component nuclear magnetic T2 spectra into different pore size intervals according to the position parameters to obtain multiple groups of component nuclear magnetic T2 spectra of different pore size intervals;
[0007] Converting the multiple groups of component nuclear magnetic T2 spectra of different pore size intervals to obtain multiple groups of component pore size spectra;
[0008] Combining the multiple groups of component pore size spectra to form a full-size pore size spectrum.
[0009] The embodiment of the present application also provides a full-size aperture characterization device for improving the accuracy of the size aperture spectrum and realizing aperture characterization under full size.
[0010] The nuclear magnetic T2 spectrum acquisition module is configured to acquire a core nuclear magnetic T2 spectrum obtained through a core nuclear magnetic resonance experiment.
[0011] The component nuclear magnetic T2 spectrum acquisition module is configured to decompose the core nuclear magnetic T2 spectrum into a plurality of component nuclear magnetic T2 spectra of different aperture intervals by using a plurality of Gaussian functions.
[0012] The aperture conversion module is configured to convert the plurality of component nuclear magnetic T2 spectra of different aperture intervals to obtain a plurality of component aperture spectra.
[0013] The full-size aperture spectrum determination module is configured to combine the plurality of component aperture spectra to form a full-size aperture spectrum.
[0014] The embodiment of the present application also provides a computer device, which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, and the processor realizes the full-size aperture characterization method when executing the computer program.
[0015] The embodiment of the present application also provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to realize the full-size aperture characterization method.
[0016] The embodiment of the present application also provides a computer program product, which comprises a computer program, and the computer program is executed by a processor to realize the full-size aperture characterization method.
[0017] In this embodiment of the invention, the core NMR T2 spectrum obtained from the core NMR experiment is acquired; the core NMR T2 spectrum is decomposed using multiple Gaussian functions to obtain multiple sets of component NMR T2 spectra for different pore size intervals; the process of decomposing the core NMR T2 spectrum into multiple sets of component NMR T2 spectra for different pore size intervals using Gaussian functions is as follows: Multiple Gaussian functions, the total number of Gaussian functions, and their position parameters are determined using the multi-derivative peak finding method; the core NMR T2 spectrum is characterized according to the multiple Gaussian functions, the total number of Gaussian functions, and their position parameters; the characterized core NMR T2 spectrum is decomposed to obtain multiple component NMR T2 spectra; the multiple component NMR T2 spectra are divided into intervals for different pore sizes according to the position parameters to obtain multiple sets of component NMR T2 spectra for different pore size intervals; the multiple sets of component NMR T2 spectra for different pore size intervals are subjected to pore size conversion to obtain multiple component pore size spectra; the multiple component pore size spectra are combined to form a full-size pore size spectrum. In the above process, this embodiment of the invention utilizes a Gaussian function to decompose the core NMR T2 spectrum into multiple component NMR T2 spectra for different pore size ranges, characterizing the core NMR T2 spectrum. The characterized core NMR T2 spectrum is then decomposed to obtain multiple component NMR T2 spectra. These multiple component NMR T2 spectra for different pore size ranges are then converted to form a full-size pore size spectrum. This improves the accuracy of the pore size spectrum and achieves full-size pore size characterization. Attached Figure Description
[0018] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. In the drawings:
[0019] Figure 1 This is a flowchart of the full-size aperture characterization method in an embodiment of the present invention;
[0020] Figure 2 These are the T2 NMR spectra of 36 shale rocks in this embodiment of the invention;
[0021] Figure 3 This is a diagram showing the decomposition results of the nuclear magnetic resonance spectrum of a certain shale in an embodiment of the present invention;
[0022] Figure 4 This is a graph showing the relationship between the geometric mean of the first component spectrum and the average pore size of nitrogen adsorption in an embodiment of the present invention.
[0023] Figure 5 This is a graph showing the relationship between the geometric mean of T2 and the average pore size of nitrogen adsorption in an embodiment of the present invention.
[0024] Figure 6This is a schematic diagram of the full-size aperture characterization device in an embodiment of the present invention. Detailed Implementation
[0025] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the embodiments of the present invention will be further described in detail below with reference to the accompanying drawings. Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.
[0026] Figure 1 This is a flowchart of a full-size aperture characterization method in an embodiment of the present invention. The method includes:
[0027] Step 101: Obtain the core nuclear magnetic resonance T2 spectrum obtained from the core nuclear magnetic resonance experiment;
[0028] Step 102 involves decomposing the core NMR T2 spectrum using multiple Gaussian functions to obtain multiple sets of component NMR T2 spectra for different pore size intervals. The process of decomposing the core NMR T2 spectrum into multiple sets of component NMR T2 spectra for different pore size intervals using Gaussian functions is as follows: Multiple Gaussian functions, their total number, and position parameters are determined using the multi-derivative peak-finding method; the core NMR T2 spectrum is characterized based on the multiple Gaussian functions, their total number, and position parameters; the characterized core NMR T2 spectrum is then decomposed to obtain multiple component NMR T2 spectra; and the multiple component NMR T2 spectra are divided into intervals for different pore sizes based on the position parameters to obtain multiple sets of component NMR T2 spectra for different pore size intervals.
[0029] Step 103: Perform aperture conversion on the component NMR T2 spectra of multiple groups with different aperture ranges to obtain multi-component component aperture spectra;
[0030] Step 104: Combine the multi-component pore size spectra to form a full-size pore size spectrum.
[0031] Each step will be explained in detail below.
[0032] In step 101, the core NMR T2 spectrum obtained from the core NMR experiment is acquired. In step 102, the core NMR T2 spectrum is decomposed using multiple Gaussian functions to obtain multiple sets of component NMR T2 spectra for different pore size intervals. The process of decomposing the core NMR T2 spectrum into multiple sets of component NMR T2 spectra for different pore size intervals using Gaussian functions is as follows: using the multi-derivative peak finding method, multiple Gaussian functions, the total number of Gaussian functions, and their position parameters are determined; the core NMR T2 spectrum is characterized according to the multiple Gaussian functions, the total number of Gaussian functions, and their position parameters; the characterized core NMR T2 spectrum is then decomposed to obtain multiple component NMR T2 spectra; the multiple component NMR T2 spectra are divided into intervals for different pore sizes according to the position parameters to obtain multiple sets of component NMR T2 spectra for different pore size intervals.
[0033] In one embodiment, the T2 NMR spectrum of the core is characterized based on multiple Gaussian functions, the total number of Gaussian functions, and position parameters, including:
[0034] By taking the second derivative of the Gaussian function, and based on the relationship between the second derivative of the Gaussian function and the amplitude and scale parameters, the amplitude and scale parameters are calculated.
[0035] The T2 NMR spectrum of the core was characterized based on amplitude parameters, scale parameters, multiple Gaussian functions, the total number of multiple Gaussian functions, and location parameters.
[0036] In a specific embodiment, based on the detection principle of nuclear magnetic resonance (NMR) experiments and the statistical principle of NMR data inversion, the aforementioned core NMR T2 spectrum can be considered as a linear superposition of multiple independent normally distributed curves. Therefore, the core NMR T2 spectrum can be decomposed into multiple component NMR T2 spectra with different pore size ranges using the Gaussian function decomposition method. In signal processing, signals conforming to a normal distribution can be characterized using a Gaussian function; therefore, the NMR T2 spectrum can be considered as a combined spectrum resulting from the superposition of multiple Gaussian function spectra. Based on the above analysis, the functional relationship for characterizing the NMR T2 spectrum using a Gaussian function can be determined as follows:
[0037]
[0038] In the formula, H is the original NMR T2 spectrum; i is the NMR inversion point number; j is the Gaussian function number; N is the number of NMR inversion points; M is the total number of Gaussian functions; x is the NMR inversion point value; u is the location parameter, which is a normal distribution symmetrical about the axis of symmetry x = u; σ is the scale parameter, which is the standard deviation of the normal distribution; a is the amplitude parameter, representing the signal amplitude when x = u; and G is the Gaussian function equation.
[0039] Therefore, the essence of Gaussian function decomposition of NMR T2 spectra is to determine the number M of Gaussian functions, the amplitude parameter a, the position parameter u, and the scale parameter σ of each Gaussian function, so as to minimize the error between the NMR T2 spectrum and the combined Gaussian function spectrum. Thus, in one embodiment, the core NMR T2 spectrum is characterized according to the following functional relationship:
[0040]
[0041] Where i is the NMR inversion point number, j is the Gaussian function number, N is the number of NMR inversion points, M is the total number of Gaussian functions, x is the NMR inversion point value, u is the location parameter, which is a normal distribution symmetrical about the axis of symmetry x = u, σ is the scale parameter, which is the standard deviation of the normal distribution; a is the amplitude parameter, representing the signal amplitude when x = u; and f(x) represents the core NMR T2 spectrum.
[0042] In one embodiment, the second derivative of the Gaussian function is taken, and the amplitude and scale parameters are calculated based on the relationship between the second derivative of the Gaussian function and the amplitude and scale parameters, including:
[0043]
[0044] Where G is the Gaussian function equation, x is the NMR inversion point value, u is the location parameter, the normal distribution is symmetrical about x = u, σ is the scale parameter, which is the standard deviation of the normal distribution; and a is the amplitude parameter, representing the signal amplitude when x = u.
[0045] In one embodiment, the method of finding peaks using multiple derivatives is used to determine multiple Gaussian functions, the total number of Gaussian functions, and their location parameters, including:
[0046] The constraints were determined using the peak-finding method based on multiple derivatives.
[0047] Based on the constraints, the total number of Gaussian functions and their positional parameters are solved.
[0048] In one embodiment, the total number and position parameters of the Gaussian function are solved based on constraints, including:
[0049] The constraints are as follows:
[0050]
[0051] Where f(x) represents the T2 NMR spectrum of the core, and x is the NMR inversion point value;
[0052] Collect each NMR inversion point that simultaneously satisfies the above constraints, count the total number of NMR inversion points, and determine the total number of NMR inversion points as the total number of Gaussian functions;
[0053] For each NMR inversion point value that simultaneously satisfies the above constraints, determine it as the location parameter of the corresponding Gaussian function.
[0054] In a specific embodiment, for core or downhole NMR, the core NMR T2 spectrum is formed by the linear summation of multiple cross-over superimposed Gaussian functions. The signal amplitude of the combined NMR T2 spectrum at x = u is not equal to the signal function amplitude of the component spectrum at that point, leading to inaccurate scale parameters. Characterizing the core NMR T2 spectrum requires determining 2M nonlinear variables (M position parameters and M scale parameters) and M linear variables (M amplitude parameters). The Nelder-Mead simplex method is used iteratively to solve for the nonlinear variables. During the iteration process, the least squares method is simultaneously used to determine the linear variables and the error term. The solution for the linear variables determines the decomposition of the core NMR T2 spectrum according to the following optimized functional relationship:
[0055]
[0056] In the formula, G is the Gaussian function equation, N is the number of NMR inversion points, M is the total number of Gaussian functions, A is the amplitude parameter of the Gaussian function to be solved, and f is the T2 NMR spectrum of the core.
[0057] At this point, the least squares method can be used to solve for the linear variables and calculate the error term ε:
[0058] A = (G T G) -1 G T f
[0059] ε=f-GA
[0060] At this time, Transformed into an optimization problem involving 2M nonlinear variables:
[0061]
[0062] Based on the above optimization solution, the component NMR T2 spectrum is obtained.
[0063] In step 103, the component NMR T2 spectra of multiple groups with different pore sizes are converted by pore size conversion to obtain multi-component component pore size spectra.
[0064] In one embodiment, multiple sets of component NMR T2 spectra from different pore size ranges are subjected to pore size conversion to obtain multi-component component pore size spectra, including:
[0065] The geometric mean of the component NMR T2 spectra of multiple groups with different aperture ranges is compared with a preset threshold.
[0066] When the geometric mean of the component NMR T2 spectrum is less than a preset threshold, the component NMR T2 spectrum is converted based on the pore size conversion relationship of the nitrogen adsorption pore size spectrum to obtain the converted component pore size spectrum; the pore size conversion equation of the nitrogen adsorption pore size spectrum is determined based on the nitrogen adsorption pore size spectrum obtained by nitrogen adsorption experiment.
[0067] When the geometric mean of the component NMR T2 spectrum is greater than a preset threshold, the component NMR T2 spectrum is converted based on the NMR T2 spectrum aperture conversion relationship of the high-pressure mercury intrusion pore spectrum to obtain the converted component aperture spectrum; the NMR T2 spectrum aperture conversion equation of the high-pressure mercury intrusion pore spectrum is determined based on the high-pressure mercury intrusion pore spectrum obtained by high-pressure mercury intrusion experiment.
[0068] In a specific embodiment, the process of pore size conversion of component NMR T2 spectra based on high-pressure mercury intrusion pore size spectrum and nitrogen adsorption pore size spectrum is as follows:
[0069] Using the average pore size d1 of the nitrogen adsorption pore size spectrum of all core samples as the dependent variable and the geometric mean T2LM of the component NMR T2 spectra of all core samples (T2LM < 10 ms) as the independent variable, the pore size conversion relationship based on the nitrogen adsorption pore size spectrum of the NMR T2 spectrum was determined:
[0070] d1 = f1(T2lm) where (T2lm < 10ms);
[0071] Using the average pore size d2 of the high-pressure mercury intrusion pore size spectrum of all core samples as the dependent variable, and the geometric mean T2LM (T2LM > 10 ms) of the component NMR T2 spectra of all core samples as the independent variable, the pore size conversion relationship based on the high-pressure mercury intrusion pore size spectrum of NMR T2 spectra was determined:
[0072] d2 = f2(T2lm) where (T2lm > 10ms);
[0073] By combining the pore size conversion relationships of NMR T2 spectra based on nitrogen adsorption pore size spectra and those based on high-pressure mercury intrusion pore size spectra, the pore size conversion relationships of component NMR T2 spectra are determined.
[0074] The full-scale pore size conversion relationship at the core scale was determined by combining nitrogen adsorption pore size spectra, high-pressure mercury intrusion pore size spectra, and nuclear magnetic resonance (NMR) T2 spectra. For macropores with pore sizes greater than 50 nm (T2LM greater than 10 ms), a NMR pore size conversion method was established by combining high-pressure mercury intrusion and NMR data. Simultaneously, for mesopores with pore sizes less than 50 nm (T2LM less than 10 ms), a NMR pore size conversion method was established by combining nitrogen adsorption and NMR data.
[0075] d1 = f1(T2lm) where (T2lm < 10ms)
[0076] d2 = f2(T2lm) where (T2lm > 10ms);
[0077] In step 104, the multi-component pore size spectra are combined to form a full-size pore size spectrum.
[0078] In a specific embodiment, the above-mentioned NMR T2 spectrum Gaussian function decomposition method and component NMR T2 spectrum pore size conversion equation are applied to NMR logging to achieve full-size pore size characterization of strongly heterogeneous reservoirs.
[0079] Taking a specific embodiment as an example, this invention is applied to the characterization of spring size and pore size in a shale reservoir of a basin. For example... Figure 2The T2 NMR spectra of 36 shale samples were analyzed. The shale T2 spectra were predominantly bimodal and trimodal. The left peak had a transverse relaxation time range of 0.05–3 ms, with a large proportion of porosity components, primarily representing NMR signals from bound water in clay or organic matter pores. The right peak (including the second or third peak) had a transverse relaxation time range of 5–1000 ms, with a smaller proportion of porosity components, primarily representing NMR signals from capillary-bound water or intergranular pores of clastic particles. Gaussian decomposition of the shale T2 spectra revealed two or three peaks in all samples. Figure 3 The diagram shows the decomposition results of the nuclear magnetic resonance spectrum of shale, as follows: Figure 3 The results of NMR spectral decomposition for a typical shale rock are presented. Through spectral decomposition, the T2 NMR spectrum of this shale rock was decomposed into three component spectra. The geometric mean values of the three component spectra are 0.43 ms, 23.6 ms, and 117 ms, respectively, with porosity components of 4.93%, 0.53%, and 0.53%. The original NMR porosity is 5.98%, and the combined NMR porosity is 5.88%, indicating that this method has good application results in actual shale core analysis.
[0080] Based on the relationship between the average pore size of nitrogen adsorption in 12 shale samples and the geometric mean of the first component spectrum and the geometric mean of the T2 spectrum, the average pore size of nitrogen adsorption is directly proportional to the geometric mean of the first component spectrum, and the average pore size increases with the increase of the geometric mean of the first component spectrum. However, the correlation between the average pore size of nitrogen adsorption and the geometric mean of the T2 spectrum is not obvious, indicating that nitrogen adsorption only characterizes the pore structure of mesopores. Meanwhile, the T2 NMR spectrum characterizes the pore structure of all pores in the shale; therefore, only at the mesopore scale do the geometric mean of the NMR spectrum and the average pore size of nitrogen adsorption show a good correspondence. Figure 4 This is a graph showing the relationship between the geometric mean of the first component spectrum and the average pore size of nitrogen adsorption. Figure 5 This is a graph showing the relationship between the geometric mean of T2 and the average pore size of nitrogen adsorption.
[0081] This invention also provides a full-size aperture characterization device, as described in the following embodiments. Since the principle by which this device solves the problem is similar to that of the full-size aperture characterization method, the implementation of this device can be found in the implementation of the method, and repeated details will not be elaborated further. Figure 6 This is a schematic diagram of a full-size aperture characterization device according to an embodiment of the present invention. The device includes:
[0082] The nuclear magnetic resonance T2 spectrum acquisition module 601 is used to acquire the core nuclear magnetic resonance T2 spectrum obtained from the core nuclear magnetic resonance experiment;
[0083] The component NMR T2 spectrum acquisition module 602 for different pore size ranges is used to decompose the core NMR T2 spectrum using multiple Gaussian functions to obtain multiple sets of component NMR T2 spectra for different pore size ranges. The process of decomposing the core NMR T2 spectrum into multiple sets of component NMR T2 spectra for different pore size ranges using Gaussian functions is as follows: Multiple Gaussian functions, the total number of Gaussian functions, and their position parameters are determined using the multi-derivative peak-finding method; the core NMR T2 spectrum is characterized based on the multiple Gaussian functions, the total number of Gaussian functions, and their position parameters; the characterized core NMR T2 spectrum is then decomposed to obtain multiple component NMR T2 spectra; and the multiple component NMR T2 spectra are divided into different pore size ranges based on the position parameters to obtain multiple sets of component NMR T2 spectra for different pore size ranges.
[0084] Aperture conversion module 603 is used to convert the component NMR T2 spectra of multiple sets of different aperture ranges to obtain multi-component component aperture spectra;
[0085] The full-size pore size spectrum determination module 604 is used to combine multi-component pore size spectra to form a full-size pore size spectrum.
[0086] In one embodiment, the component NMR T2 spectrum acquisition module 602 for different aperture ranges is specifically used for:
[0087] By taking the second derivative of the Gaussian function, and based on the relationship between the second derivative of the Gaussian function and the amplitude and scale parameters, the amplitude and scale parameters are calculated.
[0088] The T2 NMR spectrum of the core was characterized based on amplitude parameters, scale parameters, multiple Gaussian functions, the total number of multiple Gaussian functions, and location parameters.
[0089] In one embodiment, the T2 NMR spectrum of the core was characterized according to the following functional relationship:
[0090]
[0091] Where i is the NMR inversion point number, j is the Gaussian function number, N is the number of NMR inversion points, M is the total number of Gaussian functions, x is the NMR inversion point value, u is the location parameter, which is a normal distribution symmetrical about the axis of symmetry x = u, σ is the scale parameter, which is the standard deviation of the normal distribution; a is the amplitude parameter, representing the signal amplitude when x = u; and f(x) represents the core NMR T2 spectrum.
[0092] In one embodiment, the component NMR T2 spectrum acquisition module 602 for different aperture ranges is specifically used for:
[0093]
[0094] Where G is the Gaussian function equation, x is the NMR inversion point value, u is the location parameter, the normal distribution is symmetrical about x = u, σ is the scale parameter, which is the standard deviation of the normal distribution; and a is the amplitude parameter, representing the signal amplitude when x = u.
[0095] In one embodiment, the component NMR T2 spectrum acquisition module 602 for different aperture ranges is specifically used for:
[0096] The constraints were determined using the peak-finding method based on multiple derivatives.
[0097] Based on the constraints, the total number of Gaussian functions and their positional parameters are solved.
[0098] In one embodiment, the component NMR T2 spectrum acquisition module 602 for different aperture ranges is specifically used for:
[0099] The constraints are as follows:
[0100]
[0101] Where f(x) represents the T2 NMR spectrum of the core, and x is the NMR inversion point value;
[0102] Collect each NMR inversion point that simultaneously satisfies the above constraints, count the total number of NMR inversion points, and determine the total number of NMR inversion points as the total number of Gaussian functions;
[0103] For each NMR inversion point value that simultaneously satisfies the above constraints, determine it as the location parameter of the corresponding Gaussian function.
[0104] In one embodiment, the aperture conversion module 603 is specifically used for:
[0105] The geometric mean of the component NMR T2 spectra of multiple groups with different aperture ranges is compared with a preset threshold.
[0106] When the geometric mean of the component NMR T2 spectrum is less than a preset threshold, the component NMR T2 spectrum is converted based on the pore size conversion relationship of the nitrogen adsorption pore size spectrum to obtain the converted component pore size spectrum; the pore size conversion equation of the nitrogen adsorption pore size spectrum is determined based on the nitrogen adsorption pore size spectrum obtained by nitrogen adsorption experiment.
[0107] When the geometric mean of the component NMR T2 spectrum is greater than a preset threshold, the component NMR T2 spectrum is converted based on the NMR T2 spectrum aperture conversion relationship of the high-pressure mercury intrusion pore spectrum to obtain the converted component aperture spectrum; the NMR T2 spectrum aperture conversion equation of the high-pressure mercury intrusion pore spectrum is determined based on the high-pressure mercury intrusion pore spectrum obtained by high-pressure mercury intrusion experiment.
[0108] This invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the above-described full-size aperture characterization method.
[0109] This invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described full-size aperture characterization method.
[0110] This invention also provides a computer program product, which includes a computer program that, when executed by a processor, implements the above-described full-size aperture characterization method.
[0111] In this embodiment of the invention, the core NMR T2 spectrum obtained from the core NMR experiment is acquired; the core NMR T2 spectrum is decomposed using multiple Gaussian functions to obtain multiple sets of component NMR T2 spectra for different pore size intervals; the process of decomposing the core NMR T2 spectrum into multiple sets of component NMR T2 spectra for different pore size intervals using Gaussian functions is as follows: Multiple Gaussian functions, the total number of Gaussian functions, and their position parameters are determined using the multi-derivative peak finding method; the core NMR T2 spectrum is characterized according to the multiple Gaussian functions, the total number of Gaussian functions, and their position parameters; the characterized core NMR T2 spectrum is decomposed to obtain multiple component NMR T2 spectra; the multiple component NMR T2 spectra are divided into intervals for different pore sizes according to the position parameters to obtain multiple sets of component NMR T2 spectra for different pore size intervals; the multiple sets of component NMR T2 spectra for different pore size intervals are subjected to pore size conversion to obtain multiple component pore size spectra; the multiple component pore size spectra are combined to form a full-size pore size spectrum. In the above process, this embodiment of the invention utilizes a Gaussian function to decompose the core NMR T2 spectrum into multiple component NMR T2 spectra for different pore size ranges, characterizing the core NMR T2 spectrum. The characterized core NMR T2 spectrum is then decomposed to obtain multiple component NMR T2 spectra. These multiple component NMR T2 spectra for different pore size ranges are then converted to form a full-size pore size spectrum. This improves the accuracy of the pore size spectrum and achieves full-size pore size characterization.
[0112] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0113] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.
[0114] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.
[0115] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.
[0116] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above descriptions are merely specific embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.
Claims
1. A full-size aperture characterization method, characterized in that, The method comprises the following steps: obtaining a core nuclear magnetic T2 spectrum obtained by a core nuclear magnetic experiment; decomposing the core nuclear magnetic T2 spectrum into a plurality of component nuclear magnetic T2 spectra of different pore size intervals by using a plurality of Gaussian functions, wherein the decomposition of the core nuclear magnetic T2 spectrum into the plurality of component nuclear magnetic T2 spectra of different pore size intervals by using the plurality of Gaussian functions comprises the following steps: determining the plurality of Gaussian functions, the total number and the position parameters of the plurality of Gaussian functions by using a multi-order derivative peak searching method; characterizing the core nuclear magnetic T2 spectrum according to the plurality of Gaussian functions, the total number and the position parameters of the plurality of Gaussian functions, decomposing the characterized core nuclear magnetic T2 spectrum to obtain the plurality of component nuclear magnetic T2 spectra; and dividing the plurality of component nuclear magnetic T2 spectra into different pore size intervals according to the position parameters to obtain the plurality of component nuclear magnetic T2 spectra of different pore size intervals; performing pore size conversion on the plurality of component nuclear magnetic T2 spectra of different pore size intervals to obtain a plurality of component pore size spectra; combining the plurality of component pore size spectra to form a full-size pore size spectrum.
2. The method of claim 1, wherein, The characterization of the core nuclear magnetic T2 spectrum according to the plurality of Gaussian functions, the total number and the position parameters of the plurality of Gaussian functions comprises the following steps: performing second-order derivation on the Gaussian function, and calculating the amplitude parameter and the scale parameter according to the correlation between the second-order derivative of the Gaussian function and the amplitude parameter and the scale parameter; characterizing the core nuclear magnetic T2 spectrum according to the amplitude parameter, the scale parameter, the plurality of Gaussian functions, the total number and the position parameters of the plurality of Gaussian functions.
3. The method of claim 2, wherein, The core nuclear magnetic T2 spectrum is characterized according to the following equation: wherein i is a nuclear magnetic inversion point distribution number, j is a Gaussian function number, N is a nuclear magnetic inversion point number, M is a total number of Gaussian functions, x is a nuclear magnetic inversion point value, u is a position parameter, a normal distribution is symmetrically distributed about the symmetric axis x=u, σ is a scale parameter, which is a standard deviation of the normal distribution, a is an amplitude parameter, which represents a signal amplitude when x=u, and f(x) represents a core nuclear magnetic T2 spectrum.
4. The method of claim 2, wherein, The second-order derivation is performed on the Gaussian function, and the amplitude parameter and the scale parameter are calculated according to the correlation between the second-order derivative of the Gaussian function and the amplitude parameter and the scale parameter, comprising: wherein G is a Gaussian function equation, x is a nuclear magnetic inversion point value, u is a position parameter, a normal distribution is symmetrically distributed about the symmetric axis x=u, σ is a scale parameter, which is a standard deviation of the normal distribution, and a is an amplitude parameter, which represents a signal amplitude when x=u.
5. The method of claim 1, wherein, The plurality of Gaussian functions, the total number and the position parameters of the plurality of Gaussian functions are determined by using a multi-order derivative peak searching method, comprising: determining a constraint condition by using the multi-order derivative peak searching method; solving the total number and the position parameters of the Gaussian functions based on the constraint condition.
6. The method of claim 5, wherein, The solving of the total number and the position parameters of the Gaussian functions based on the constraint condition comprises: the constraint condition is: wherein f(x) represents a core nuclear magnetic T2 spectrum, and x is a nuclear magnetic inversion point value; collecting each nuclear magnetic inversion point satisfying the above constraint condition, counting the total number of the nuclear magnetic inversion points, and determining the total number of the nuclear magnetic inversion points as the total number of the Gaussian functions; respectively determining each nuclear magnetic inversion point value satisfying the above constraint condition as the position parameter of the corresponding Gaussian function.
7. The method of claim 1, wherein, The component nuclear magnetic T2 spectra of multiple groups of different pore size intervals are subjected to pore size conversion to obtain multiple groups of component pore size spectra, including: The geometric mean of the component nuclear magnetic T2 spectra of multiple groups of different pore size intervals is compared with a preset threshold value; When the geometric mean of the component nuclear magnetic T2 spectrum is less than the preset threshold value, the component nuclear magnetic T2 spectrum is subjected to pore size conversion based on a nuclear magnetic T2 spectrum pore size conversion relationship of a nitrogen adsorption pore size spectrum to obtain a converted component pore size spectrum; the nuclear magnetic T2 spectrum pore size conversion relationship of the nitrogen adsorption pore size spectrum is determined according to a nitrogen adsorption pore size spectrum obtained by using a nitrogen adsorption experiment; When the geometric mean of the component nuclear magnetic T2 spectrum is greater than the preset threshold value, the component nuclear magnetic T2 spectrum is subjected to pore size conversion based on a nuclear magnetic T2 spectrum pore size conversion equation of a high-pressure mercury intrusion pore size spectrum to obtain a converted component pore size spectrum; the nuclear magnetic T2 spectrum pore size conversion equation of the high-pressure mercury intrusion pore size spectrum is determined according to a high-pressure mercury intrusion pore size spectrum obtained by using a high-pressure mercury intrusion experiment.
8. A full-size aperture characterization device, characterized in that, including: The core nuclear magnetic T2 spectrum acquisition module is configured to acquire a core nuclear magnetic T2 spectrum obtained by a core nuclear magnetic resonance experiment; The component nuclear magnetic T2 spectrum acquisition module is configured to decompose the core nuclear magnetic T2 spectrum by using multiple Gaussian functions to obtain multiple groups of component nuclear magnetic T2 spectra of different pore size intervals; the process of decomposing the core nuclear magnetic T2 spectrum into multiple groups of component nuclear magnetic T2 spectra of different pore size intervals by using the Gaussian functions includes the following steps: determining multiple Gaussian functions, a total number and position parameters of the multiple Gaussian functions by using a multiple-order derivative peak searching method; characterizing the core nuclear magnetic T2 spectrum according to the multiple Gaussian functions, the total number and the position parameters of the multiple Gaussian functions; and decomposing the characterized core nuclear magnetic T2 spectrum to obtain multiple component nuclear magnetic T2 spectra; and dividing the multiple component nuclear magnetic T2 spectra into different pore size intervals according to the position parameters to obtain multiple groups of component nuclear magnetic T2 spectra of different pore size intervals; The pore size conversion module is configured to perform pore size conversion on the multiple groups of component nuclear magnetic T2 spectra of different pore size intervals to obtain multiple groups of component pore size spectra; The full-size pore size spectrum determination module is configured to combine the multiple groups of component pore size spectra to form a full-size pore size spectrum.
9. The apparatus of claim 8, wherein, The component nuclear magnetic T2 spectrum acquisition module is specifically configured to: perform second-order derivation on the Gaussian functions, calculate the amplitude parameters and the scale parameters according to a relationship between the second-order derivatives of the Gaussian functions and the amplitude parameters and the scale parameters; and characterize the core nuclear magnetic T2 spectrum according to the amplitude parameters, the scale parameters, the multiple Gaussian functions, a total number and position parameters of the multiple Gaussian functions.
10. The apparatus of claim 9, wherein, The core nuclear magnetic T2 spectrum is characterized according to the following function relationship: wherein i is a nuclear magnetic inversion point distribution number, j is a Gaussian function number, N is a nuclear magnetic inversion point number, M is a total number of Gaussian functions, x is a nuclear magnetic inversion point value, u is a position parameter, a normal distribution is symmetrically distributed about the symmetric axis x=u, σ is a scale parameter, which is a standard deviation of the normal distribution; a is an amplitude parameter, which represents a signal amplitude when x=u; and f(x) represents the core nuclear magnetic T2 spectrum.
11. The apparatus of claim 9, wherein, The component nuclear magnetic T2 spectrum acquisition module is specifically configured to: Wherein, G is a Gaussian function equation, x is a nuclear magnetic inversion point value, u is a position parameter, a normal distribution is symmetrically distributed with x=u as a symmetric axis, and sigma is a scale parameter, which is a standard deviation of a normal distribution; a is an amplitude parameter, which represents a signal amplitude when x=u.
12. The apparatus of claim 8, wherein, The component nuclear magnetic T2 spectrum acquisition module of different pore size intervals is specifically used for: The constraint condition is determined by a multi-order derivative peak finding method; Based on the constraint condition, the total number of Gaussian functions and the position parameters are solved.
13. The apparatus of claim 12, wherein, The component nuclear magnetic T2 spectrum acquisition module of different pore size intervals is specifically used for: The constraint condition is: Wherein, f(x) represents a core nuclear magnetic T2 spectrum, and x is a nuclear magnetic inversion point value; Each nuclear magnetic inversion point value that meets the above constraint condition is determined as a position parameter of a corresponding Gaussian function. The pore size conversion module is specifically used for:
14. The apparatus of claim 8, wherein, Comparing the geometric mean value of a plurality of groups of component nuclear magnetic T2 spectra of different pore size intervals with a preset threshold value; When the geometric mean value of the component nuclear magnetic T2 spectrum is less than the preset threshold value, the component nuclear magnetic T2 spectrum is subjected to pore size conversion based on a nuclear magnetic T2 spectrum pore size conversion relationship of a nitrogen adsorption pore size spectrum, to obtain a converted component pore size spectrum; the nuclear magnetic T2 spectrum pore size conversion equation of the nitrogen adsorption pore size spectrum is determined according to a nitrogen adsorption pore size spectrum obtained by using a nitrogen adsorption experiment; When the geometric mean value of the component nuclear magnetic T2 spectrum is greater than the preset threshold value, the component nuclear magnetic T2 spectrum is subjected to pore size conversion based on a nuclear magnetic T2 spectrum pore size conversion relationship of a high-pressure mercury intrusion pore size spectrum, to obtain a converted component pore size spectrum; the nuclear magnetic T2 spectrum pore size conversion equation of the high-pressure mercury intrusion pore size spectrum is determined according to a high-pressure mercury intrusion pore size spectrum obtained by using a high-pressure mercury intrusion experiment. The processor executes the computer program to realize the method of any one of claims 1 to 7.
15. A computer device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The computer readable storage medium stores a computer program, and the computer program is executed by the processor to realize the method of any one of claims 1 to 7.
16. A computer-readable storage medium, characterized in that, The computer program product includes a computer program, and the computer program is executed by the processor to realize the method of any one of claims 1 to 7.
17. A computer program product, characterised in that,