A method, system, and program product for pressure-holding coring nuclear porosimetry

Abstract

The present application belongs to the technical field of nuclear magnetic resonance logging, and specifically discloses a pressure-maintaining coring nuclear magnetic porosity measurement method, system and program product. Through a corresponding convex optimization process, the amount of pollution signals of the pressure-maintaining liquid is estimated and the T2 relaxation distribution of the core is corrected, so that the pressure-maintaining liquid relaxation component in the original signal is removed without affecting the original relaxation signal of the core, and the measured T2 relaxation distribution is close to the true distribution of the core. By canceling the pressure-maintaining liquid signal component in the measurement signal, the present application can effectively reduce the interference of the pressure-maintaining liquid signal on the core porosity measurement during pressure-maintaining coring, improve the core porosity measurement accuracy, and further realize in-situ high-precision porosity measurement of the pressure-maintaining coring sample by nuclear magnetic logging.

Description

Technical Field

[0001] This invention belongs to the field of nuclear magnetic resonance logging technology, specifically relating to a method, system, and program product for measuring nuclear magnetic porosity using pressure-maintaining coring. Background Technology

[0002] Nuclear magnetic resonance logging requires precise measurement of core porosity to determine the content of oil and gas reservoirs. Conventional coring methods have unavoidable problems such as gas escape and changes in the physicochemical environment. Pressure-holding coring can better preserve the original state of the core, but the pressure-holding fluid in the pressure-holding cylinder can interfere with the porosity measurement.

[0003] Existing well-side measurement techniques typically employ freezing and then extracting cored samples for measurement. This method cannot completely prevent gas escape and changes in the physicochemical environment. On the other hand, if an in-situ measurement scheme is adopted, it is impossible to avoid the interference of the NMR signal of the pressure-holding liquid on porosity measurement. Summary of the Invention

[0004] The purpose of this invention is to provide a method, system, and program for measuring the porosity of nuclear magnetic resonance imaging (NMR) samples obtained by pressure-holding core sampling, in order to solve the aforementioned problems existing in the prior art.

[0005] To achieve the above objectives, the present invention adopts the following technical solution: Firstly, a method for measuring nuclear magnetic resonance porosity using pressure-holding coring is provided, comprising: The magnetic resonance T2 relaxation signal of the core sample was acquired by the nuclear magnetic resonance measurement sequence and the cross-sectional image was constructed using the magnetic resonance T2 relaxation signal; Using the center of the cross-sectional image as the center, select the central region with a set radius and the outer ring region outside the central region; The average gray level of the central region is used as the first reference signal value, and the average gray level of the outer ring region is used as the second reference signal value. Obtain the set prior T2 relaxation time distribution vector, and use the prior T2 relaxation time distribution vector, the first reference signal value, and the second reference signal value to determine the amount of pollution signal that satisfies the set first convex optimization condition; The corrected T2 relaxation time distribution vector that satisfies the set second convex optimization condition is determined by using the pollution signal quantity, the first reference signal value, and the second reference signal value. The total porosity is calculated based on the corrected T2 relaxation time distribution vector, and the core pore distribution parameters are calculated using the corrected T2 relaxation time distribution vector, the total porosity, and the set positive correlation ratio. Construct a core pore distribution curve based on the core pore distribution parameters, and output the overall porosity and core pore distribution curves.

[0006] In one possible design, the method further includes: Obtain the prior relaxation signal value of the pressure-holding coring fluid used in the pressure-holding coring sample, and use the prior relaxation signal value as the second reference signal value.

[0007] In one possible design, determining the amount of contamination signal satisfying the set first convex optimization condition using the prior T2 relaxation time distribution vector, the first reference signal value, and the second reference signal value includes: Substituting the prior T2 relaxation time distribution vector, the first reference signal value, and the second reference signal value into the first convex optimization condition, the pollution signal quantity satisfying the first convex optimization condition is obtained. The first convex optimization condition is:

[0008] in, To contaminate the signal quantity, The value obtained from the characterization calculation is the smallest. Let T2 be the prior relaxation time distribution vector. The first reference signal value, The second reference signal value, For the set T2 attenuation matrix, To satisfy convexity and about Decreasing loss function.

[0009] In one possible design, the loss function Set as

[0010] in, This is the first regularization hyperparameter set.

[0011] In one possible design, determining the corrected T2 relaxation time distribution vector that satisfies the set second convex optimization condition using the contamination signal quantity, the first reference signal value, and the second reference signal value includes: Substituting the pollution signal quantity, the first reference signal value, and the second reference signal value into the second convex optimization condition, the corrected T2 relaxation time distribution vector satisfying the second convex optimization condition is obtained. The second convex optimization condition is:

[0012] in, To correct the T2 relaxation time distribution vector, The value obtained from the characterization calculation is the smallest. The first reference signal value, The second reference signal value, For the set T2 attenuation matrix, The second regularization hyperparameter is set. This is the third regularization hyperparameter that is set.

[0013] In one possible design, the calculation of overall porosity based on the corrected T2 relaxation time distribution vector includes: The components of the corrected T2 relaxation time distribution vector are added together to obtain the total T2 relaxation distribution S1; Substitute the sum of T2 relaxation distributions S1 into the preset overall porosity calculation formula to obtain the overall porosity P. The overall porosity calculation formula is P = (S1 / S0) × P0, where S0 is the set standard value of the sum of relaxation distributions and P0 is the set standard porosity.

[0014] In one possible design, the calculation of core pore distribution parameters using the corrected T2 relaxation time distribution vector, total porosity, and a set positive correlation ratio includes: Substitute each component of the corrected T2 relaxation time distribution vector into the preset core porosity scatter plot formula K. i = (x i The calculation is performed in (×P×β) / S1 to obtain the scatter plot parameters of core porosity corresponding to each component, where x i To correct the corresponding components of the T2 relaxation time distribution vector, where i is the component index and K is the component number. i Here are the corresponding scatter plot parameters for core porosity, and β is the set positive correlation ratio value. Core porosity distribution parameters are obtained by combining the scattered parameters of porosity from each core sample.

[0015] Secondly, a pressure-holding coring nuclear magnetic resonance porosity measurement system is provided, comprising a signal acquisition unit, a region selection unit, a signal determination unit, a first optimization unit, a second optimization unit, a porosity calculation unit, and a measurement output unit, wherein: The signal acquisition unit is used to acquire the magnetic resonance T2 relaxation signal of the core sample based on the nuclear magnetic resonance measurement sequence, and to construct a cross-sectional image using the magnetic resonance T2 relaxation signal; The region selection unit is used to select a central region with a set radius and an outer ring region outside the central region, with the center of the cross-sectional image as the center. The signal determination unit is used to take the average gray value of the image in the central region as the first reference signal value and the average gray value of the image in the outer ring region as the second reference signal value. The first optimization unit is used to obtain the set prior T2 relaxation time distribution vector, and use the prior T2 relaxation time distribution vector, the first reference signal value and the second reference signal value to determine the amount of pollution signal that satisfies the set first convex optimization condition. The second optimization unit is used to determine the corrected T2 relaxation time distribution vector that satisfies the set second convex optimization condition by using the pollution signal quantity, the first reference signal value and the second reference signal value. The porosity calculation unit is used to calculate the overall porosity based on the corrected T2 relaxation time distribution vector, and to calculate the core pore distribution parameters using the corrected T2 relaxation time distribution vector, the overall porosity, and the set positive correlation ratio. The measurement output unit is used to construct a core pore distribution curve based on the core pore distribution parameters, and output the overall porosity and core pore distribution curve.

[0016] Thirdly, a pressure-holding coring nuclear magnetic resonance porosity measurement system is provided, comprising: Memory, used to store instructions; A processor is configured to read instructions stored in the memory and execute the method described in any one of the first aspects above, according to the instructions.

[0017] Fourthly, a computer-readable storage medium is provided, on which instructions are stored, which, when executed on a computer, cause the computer to perform any of the methods described in the first aspect. A computer program product is also provided, which, when executed on a computer, performs any of the methods described in the first aspect.

[0018] Beneficial effects: By canceling the pressure-holding fluid signal component in the measurement signal, this invention can effectively reduce the interference of the pressure-holding fluid signal on the core porosity measurement during pressure-holding coring, improve the accuracy of core porosity measurement, and thus realize in-situ high-precision porosity measurement of pressure-holding cored samples near nuclear magnetic resonance wells. Attached Figure Description

[0019] 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.

[0020] Figure 1 This is a schematic diagram of the steps in the method of Embodiment 1 of the present invention; Figure 2 This is a schematic diagram of the system configuration in Embodiment 2 of the present invention; Figure 3 This is a schematic diagram of the system configuration in Embodiment 3 of the present invention. Detailed Implementation

[0021] It should be noted that the descriptions of these embodiments are intended to aid in understanding the invention and do not constitute a limitation thereof. The specific structural and functional details disclosed herein are merely for describing exemplary embodiments of the invention. However, the invention may be embodied in many alternative forms and should not be construed as being limited to the embodiments described herein.

[0022] It should be understood that, unless otherwise explicitly specified and limited, the term "connection" should be interpreted broadly. For example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium; it can be a connection within two components. Those skilled in the art can understand the specific meaning of the above terms in the embodiments according to the specific circumstances.

[0023] Specific details are provided in the following description to provide a complete understanding of the exemplary embodiments. However, those skilled in the art will understand that the exemplary embodiments can be implemented without these specific details. For example, apparatus may be shown in block diagrams to avoid obscuring the examples with unnecessary details. In other embodiments, well-known processes, structures, and techniques may be omitted with non-essential details to avoid obscuring the embodiments.

[0024] Example 1: This embodiment provides a method for measuring porosity in pressure-holding coring nuclear magnetic resonance imaging, such as... Figure 1 As shown, the method includes the following steps: S1. The magnetic resonance T2 relaxation signal of the core sample was acquired based on the nuclear magnetic resonance measurement sequence, and the cross-sectional image was constructed using the magnetic resonance T2 relaxation signal.

[0025] In practice, NMR sequences can be used to acquire the T2 relaxation signals of cored samples obtained during pressure holding. NMR logging typically uses CPMG pulse sequences, and the measured data is the T2 relaxation signal of the excited object over time. Based on this signal, the content of components with different T2 decay rates within the object can be determined. Spatially resolution NMR sequences can be used for acquisition. The results obtained from spatially resolution NMR sequences contain a set of T2 relaxation image signals of the object changing over time. Each frame represents the relaxation state of the object's cross-section at the current moment, containing richer information than a standard CPMG sequence. After acquiring the T2 relaxation signals of the cored samples, cross-sectional images can be constructed using these signals. Typically, a non-uniform Fourier transform or CT projection method corresponding to the k-space acquisition trajectory of the pulse sequence is used to reconstruct the cross-sectional image from the k-space NMR signal. For spatially resolution NMR sequences, each echo time corresponds to one cross-sectional image.

[0026] S2. Using the center of the cross-sectional image as the center, select the central region with a set radius and the outer ring region outside the central region.

[0027] In practice, after obtaining the cross-sectional image, a central region with a set radius can be selected with the center of the cross-sectional image as the center. Step S1 uses an NMR measurement sequence with spatial resolution to acquire T2 relaxation signals, and it is also necessary to select the region outside the central region that contains the pressure-holding fluid signal as the outer ring region.

[0028] S3. Take the average gray value of the central region as the first reference signal value, and take the average gray value of the outer ring region as the second reference signal value.

[0029] In practice, the grayscale value of each pixel in the central region is determined. Then, the average grayscale value of the central region is used as the first reference signal value sig, and the average grayscale value of the outer ring region is used as the second reference signal value b0, which is the NMR T2 relaxation signal of the pressure-holding fluid. Alternatively, the prior relaxation signal value of the pressure-holding core sample used can be directly obtained, and this prior relaxation signal value can be used as the second reference signal value b0.

[0030] S4. Obtain the set prior T2 relaxation time distribution vector, and use the prior T2 relaxation time distribution vector, the first reference signal value, and the second reference signal value to determine the amount of contamination signal that satisfies the set first convex optimization condition.

[0031] In specific implementation, firstly, the prior T2 relaxation time distribution vector is obtained. This prior T2 relaxation time distribution vector can be an empirical vector pre-defined based on prior knowledge. Then, the prior T2 relaxation time distribution vector, the first reference signal value, and the second reference signal value are substituted into the first convex optimization condition for calculation, resulting in the pollution signal quantity that satisfies the first convex optimization condition. The first convex optimization condition is...

[0032] in, This is a contamination signal; The value obtained from the characterization calculation is the smallest; The prior T2 relaxation time distribution vector; This is the first reference signal value; This is the second reference signal value; Let A = exp(-t_i / T2_j) be the set T2 attenuation matrix, which is the forward discrete Laplace transform, converting the pore component distribution into a time-domain T2 attenuation signal. t_i represents the echo time series, and T2_j represents the discretized T2 relaxation time sampling point sequence. This linear transformation can be regarded as an attenuation matrix. To satisfy convexity and about A decreasing loss function, exemplarily, the loss function Set as

[0033] in, This is the first regularization hyperparameter set. Since the pressure-holding fluid signal is at a higher level than the core signal, this loss function is used in the optimization fit to... The goal is to favor larger values ​​to improve the accuracy of the optimization.

[0034] S5. Determine the corrected T2 relaxation time distribution vector that satisfies the set second convex optimization condition using the pollution signal quantity, the first reference signal value, and the second reference signal value.

[0035] In practice, the pollution signal quantity, the first reference signal value, and the second reference signal value can be substituted into the second convex optimization condition for calculation to obtain the corrected T2 relaxation time distribution vector that satisfies the second convex optimization condition. The second convex optimization condition is:

[0036] in, To correct the T2 relaxation time distribution vector, The value obtained from the characterization calculation is the smallest. The first reference signal value, The second reference signal value, For the set T2 attenuation matrix, The second regularization hyperparameter is set. This is the third regularization hyperparameter that is set.

[0037] S6. Calculate the overall porosity based on the corrected T2 relaxation time distribution vector, and calculate the core pore distribution parameters using the corrected T2 relaxation time distribution vector, the overall porosity, and the set positive correlation ratio.

[0038] In practice, the components of the corrected T2 relaxation time distribution vector can be added together to obtain the total T2 relaxation distribution S1. Then, the total T2 relaxation distribution S1 is substituted into a preset overall porosity calculation formula to obtain the overall porosity P. The overall porosity calculation formula is P = (S1 / S0) × P0, where S0 is the set standard value of the total relaxation distribution, and P0 is the set standard porosity. For example, the total T2 relaxation distribution S0 of a copper sulfate solution standard sample with a determined overall porosity P0 can be used as the standard value. For pure solution samples, P0 = 100%.

[0039] After calculating the overall porosity P, the components of the corrected T2 relaxation time distribution vector can be substituted into the preset core porosity scatter plot formula K. i = (x iThe calculation is performed in (×P×β) / S1 to obtain the scatter plot parameters of core porosity corresponding to each component, where x i To correct the corresponding components of the T2 relaxation time distribution vector, where i is the component index and K is the component number. i Here, β represents the corresponding scatter plot parameter of core porosity, and β is a set positive correlation ratio value. Then, the core porosity distribution parameter is obtained by combining the scatter plot parameters of each core.

[0040] S7. Construct a core pore distribution curve based on the core pore distribution parameters, and output the overall porosity and core pore distribution curves.

[0041] In practice, a scatter plot of core pore distribution can be constructed using core pore distribution parameters, and curve fitting can be performed on the scatter plot to obtain a core pore distribution curve. Finally, the overall porosity and core pore distribution curve are output to achieve pressure-controlled core sampling and nuclear magnetic resonance porosity measurement.

[0042] Typically, during pressure-holding coring, the core sample is stored in a pressure-holding cylinder filled with pressure-holding fluid. Therefore, the signal measured by NMR is the sum of the interference signal from the pressure-holding fluid and the pore fluid in the core sample, resulting in an inflated NMR relaxation signal. Furthermore, abnormal interference peaks appear in the pressure-holding fluid signal region during T2 relaxation (which may overlap with the core signal and cannot be directly removed), leading to an inflated porosity measurement. The method described in this embodiment can partially or completely remove the pressure-holding fluid relaxation component from the original signal without affecting the original relaxation signal of the core sample. This ensures that the measured T2 relaxation distribution closely approximates the true distribution of the core sample, resulting in a more accurate porosity measurement and improved core porosity measurement precision. Ultimately, this enables in-situ, high-precision porosity measurement of pressure-holding cored samples near the well site using NMR.

[0043] Example 2: This embodiment provides a pressure-holding coring nuclear magnetic resonance porosity measurement system, such as... Figure 2 As shown, it includes a signal acquisition unit, a region selection unit, a signal determination unit, a first optimization unit, a second optimization unit, a porosity calculation unit, and a measurement output unit, wherein: The signal acquisition unit is used to acquire the magnetic resonance T2 relaxation signal of the core sample based on the nuclear magnetic resonance measurement sequence, and to construct a cross-sectional image using the magnetic resonance T2 relaxation signal; The region selection unit is used to select a central region with a set radius and an outer ring region outside the central region, with the center of the cross-sectional image as the center. The signal determination unit is used to take the average gray value of the image in the central region as the first reference signal value and the average gray value of the image in the outer ring region as the second reference signal value. The first optimization unit is used to obtain the set prior T2 relaxation time distribution vector, and use the prior T2 relaxation time distribution vector, the first reference signal value and the second reference signal value to determine the amount of pollution signal that satisfies the set first convex optimization condition. The second optimization unit is used to determine the corrected T2 relaxation time distribution vector that satisfies the set second convex optimization condition by using the pollution signal quantity, the first reference signal value and the second reference signal value. The porosity calculation unit is used to calculate the overall porosity based on the corrected T2 relaxation time distribution vector, and to calculate the core pore distribution parameters using the corrected T2 relaxation time distribution vector, the overall porosity, and the set positive correlation ratio. The measurement output unit is used to construct a core pore distribution curve based on the core pore distribution parameters, and output the overall porosity and core pore distribution curve.

[0044] Example 3: This embodiment provides a pressure-holding coring nuclear magnetic resonance porosity measurement system, such as... Figure 3 As shown, at the hardware level, it includes: The data interface is used to establish data communication between the processor and external data terminals; Memory, used to store instructions; The processor is used to read the instructions stored in the memory and execute the pressure-holding core sampling nuclear magnetic porosity measurement method in Embodiment 1 according to the instructions.

[0045] Optionally, the system also includes an internal bus, through which the processor, memory, and data interface can be interconnected. This internal bus can be an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus, or an EISA (Extended Industry Standard Architecture) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc.

[0046] The memory may include, but is not limited to, random access memory (RAM), read-only memory (ROM), flash memory, first-in-first-out (FIFO) memory, and / or first-in-last-out (FILO) memory. The processor may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it may also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0047] Example 4: This embodiment provides a computer-readable storage medium storing instructions. When these instructions are executed on a computer, the computer performs the pressure-holding core sampling NMR porosity measurement method described in Embodiment 1. The computer-readable storage medium refers to a data storage medium, which may include, but is not limited to, floppy disks, optical disks, hard disks, flash memory, USB flash drives, and / or Memory Sticks. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices.

[0048] This embodiment also provides a computer program product that, when run on a computer, executes the pressure-holding coring NMR porosity measurement method described in Embodiment 1. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device.

[0049] Finally, it should be noted that the above description is merely a preferred embodiment of the present invention and is 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 method for measuring nuclear magnetic resonance porosity using pressure-holding coring, characterized in that, include: The magnetic resonance T2 relaxation signal of the core sample was acquired based on the NMR measurement sequence with spatial resolution, and the cross-sectional image was constructed using the magnetic resonance T2 relaxation signal. Using the center of the cross-sectional image as the center, select the central region with a set radius and the outer ring region outside the central region; The average gray level of the central region is used as the first reference signal value, and the average gray level of the outer ring region is used as the second reference signal value. Obtain the set prior T2 relaxation time distribution vector, and use the prior T2 relaxation time distribution vector, the first reference signal value, and the second reference signal value to determine the amount of pollution signal that satisfies the set first convex optimization condition, including: Substituting the prior T2 relaxation time distribution vector, the first reference signal value, and the second reference signal value into the first convex optimization condition, the pollution signal quantity satisfying the first convex optimization condition is obtained. The first convex optimization condition is: in, To contaminate the signal quantity, The value obtained from the characterization calculation is the smallest. Let T2 be the prior relaxation time distribution vector. The first reference signal value, The second reference signal value, For the set T2 attenuation matrix, To satisfy convexity and about Decreasing loss function, set as in, This is the first regularization hyperparameter set; The corrected T2 relaxation time distribution vector that satisfies the set second convex optimization condition is determined by using the pollution signal quantity, the first reference signal value, and the second reference signal value. The total porosity is calculated based on the corrected T2 relaxation time distribution vector, and the core pore distribution parameters are calculated using the corrected T2 relaxation time distribution vector, the total porosity, and the set positive correlation ratio. Construct a core pore distribution curve based on the core pore distribution parameters, and output the overall porosity and core pore distribution curves.

2. The method for measuring nuclear magnetic resonance porosity using pressure-holding coring according to claim 1, characterized in that, The step of determining the corrected T2 relaxation time distribution vector that satisfies the set second convex optimization condition using the pollution signal quantity, the first reference signal value, and the second reference signal value includes: Substituting the pollution signal quantity, the first reference signal value, and the second reference signal value into the second convex optimization condition, the corrected T2 relaxation time distribution vector satisfying the second convex optimization condition is obtained. The second convex optimization condition is: in, To correct the T2 relaxation time distribution vector, The value obtained from the characterization calculation is the smallest. The first reference signal value, The second reference signal value, For the set T2 attenuation matrix, The second regularization hyperparameter is set. This is the third regularization hyperparameter that is set.

3. The method for measuring nuclear magnetic resonance porosity using pressure-holding coring according to claim 1, characterized in that, The calculation of total porosity based on the corrected T2 relaxation time distribution vector includes: The components of the corrected T2 relaxation time distribution vector are added together to obtain the total T2 relaxation distribution S1; Substitute the sum of T2 relaxation distributions S1 into the preset overall porosity calculation formula to obtain the overall porosity P. The overall porosity calculation formula is P = (S1 / S0) × P0, where S0 is the set standard value of the sum of relaxation distributions and P0 is the set standard porosity.

4. The method for measuring nuclear magnetic resonance porosity in pressure-holding coring according to claim 3, characterized in that, The calculation of core pore distribution parameters using the corrected T2 relaxation time distribution vector, total porosity, and a set positive correlation ratio includes: Substitute each component of the corrected T2 relaxation time distribution vector into the preset core porosity scatter plot formula K. i = (x i The calculation is performed in (×P×β) / S1 to obtain the scatter plot parameters of core porosity corresponding to each component, where x i To correct the corresponding components of the T2 relaxation time distribution vector, where i is the component index and K is the component number. i Here are the corresponding scatter plot parameters for core porosity, and β is the set positive correlation ratio value. Core porosity distribution parameters are obtained by combining the scattered parameters of porosity from each core sample.

5. A pressure-holding coring NMR porosity measurement system, used to implement the pressure-holding coring NMR porosity measurement method according to any one of claims 1-4, characterized in that, It includes a signal acquisition unit, a region selection unit, a signal determination unit, a first optimization unit, a second optimization unit, a porosity calculation unit, and a measurement output unit, wherein: The signal acquisition unit is used to acquire the magnetic resonance T2 relaxation signal of the core sample based on the nuclear magnetic resonance measurement sequence, and to construct a cross-sectional image using the magnetic resonance T2 relaxation signal; The region selection unit is used to select a central region with a set radius and an outer ring region outside the central region, with the center of the cross-sectional image as the center. The signal determination unit is used to take the average gray value of the image in the central region as the first reference signal value and the average gray value of the image in the outer ring region as the second reference signal value. The first optimization unit is used to obtain a predetermined prior T2 relaxation time distribution vector, and to determine the amount of contamination signal that satisfies the predetermined first convex optimization condition using the prior T2 relaxation time distribution vector, the first reference signal value, and the second reference signal value, including: Substituting the prior T2 relaxation time distribution vector, the first reference signal value, and the second reference signal value into the first convex optimization condition, the pollution signal quantity satisfying the first convex optimization condition is obtained. The first convex optimization condition is: in, To contaminate the signal quantity, The value obtained from the characterization calculation is the smallest. Let T2 be the prior relaxation time distribution vector. The first reference signal value, The second reference signal value, For the set T2 attenuation matrix, To satisfy convexity and about Decreasing loss function, set as in, This is the first regularization hyperparameter set; The second optimization unit is used to determine the corrected T2 relaxation time distribution vector that satisfies the set second convex optimization condition by using the pollution signal quantity, the first reference signal value and the second reference signal value. The porosity calculation unit is used to calculate the overall porosity based on the corrected T2 relaxation time distribution vector, and to calculate the core pore distribution parameters using the corrected T2 relaxation time distribution vector, the overall porosity, and the set positive correlation ratio. The measurement output unit is used to construct a core pore distribution curve based on the core pore distribution parameters, and output the overall porosity and core pore distribution curve.

6. A pressure-holding coring nuclear magnetic resonance porosity measurement system, characterized in that, include: Memory, used to store instructions; A processor is configured to read instructions stored in the memory and execute the pressure-holding coring nuclear magnetic resonance porosity measurement method according to any one of claims 1-4.

7. A computer program product, characterized in that, When the computer program product is run on a computer, it executes the pressure-holding coring nuclear magnetic resonance porosity measurement method according to any one of claims 1-4.

Images