Method and device for determining full-size distribution curve of pore radius of rock

Abstract

The invention discloses a method and device for determining a rock pore radius full-size distribution curve, and the method comprises the steps: obtaining a rock sample, and pretreating the rock sample to obtain a standard rock core, the pretreatment at least comprising oil washing treatment and / or drying treatment; determining gas logging pore information of the standard rock core, and acquiring rock core scanning data of the standard rock core through computed tomography equipment; processing the standard rock core into a saturated water rock core, and obtaining a transverse relaxation time map of the saturated water rock core through magnetic resonance equipment; a rock pore radius distribution curve corresponding to the rock sample is determined according to the rock core scanning data, the gas logging pore information and the transverse relaxation time atlas, the precision and reliability of pore structure analysis are improved, and a scientific basis is provided for optimizing a water drive scheme and improving the use degree of crude oil in tiny pores.

Description

Technical Field

[0001] The present invention relates to the technical field of geological data processing, and particularly relates to a method and device for determining the full-size distribution curve of rock pore radius. Background Art

[0002] In the fields of geology and petroleum engineering, accurately characterizing the pore structure of carbonate reservoirs, especially pores at the micron scale and smaller, is crucial for improving water flooding recovery efficiency and crude oil production.

[0003] In related technologies, pore structure analysis methods, such as mercury intrusion porosimetry and gas adsorption methods, although they can provide certain pore information, have problems such as complex operation, sample destruction, and limited data accuracy. Especially under the conditions of complex pore structures, it is difficult to comprehensively and accurately reflect the pore distribution characteristics of rocks. Summary of the Invention

[0004] The present invention provides a method and device for determining the full-size distribution curve of rock pore radius to solve the problem that it is difficult to reflect the pore distribution characteristics of rocks in related technologies under the condition of complex pore structures.

[0005] According to one aspect of the present invention, a method for determining the full-size distribution curve of rock pore radius is provided, including:

[0006] Obtain a rock sample, and perform pretreatment on the rock sample to obtain a standard core, wherein the pretreatment includes at least oil washing treatment and / or drying treatment;

[0007] Determine the gas-measured pore information of the standard core, and obtain the core scan data of the standard core through a computed tomography device;

[0008] Process the standard core into a water-saturated core, and obtain the transverse relaxation time spectrum of the water-saturated core through a magnetic resonance device;

[0009] Determine the rock pore radius distribution curve corresponding to the rock sample according to the core scan data, the gas-measured pore information, and the transverse relaxation time spectrum.

[0010] According to another aspect of the present invention, a device for determining the full-size distribution curve of rock pore radius is provided, including:

[0011] A standard core acquisition module, configured to obtain a rock sample, and perform pretreatment on the rock sample to obtain a standard core, wherein the pretreatment includes at least oil washing treatment and / or drying treatment;

[0012] A standard core analysis module, configured to determine the gas-measured pore information of the standard core, and obtain the core scan data of the standard core through a computed tomography device;

[0013] A saturated water core analysis module for processing the standard core into a saturated water core and obtaining a transverse relaxation time spectrum of the saturated water core through a magnetic resonance device;

[0014] A pore radius distribution determination module for determining a rock pore radius distribution curve corresponding to the rock sample according to the core scan data, the gas logging pore information, and the transverse relaxation time spectrum.

[0015] According to another aspect of the present invention, there is provided an electronic device, which includes:

[0016] At least one processor; and

[0017] A memory communicatively connected to the at least one processor; wherein,

[0018] The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for determining the full-size distribution curve of the rock pore radius according to any embodiment of the present invention.

[0019] According to another aspect of the present invention, there is provided a computer-readable storage medium storing computer instructions for causing a processor to implement the method for determining the full-size distribution curve of the rock pore radius according to any embodiment of the present invention when executed.

[0020] According to another aspect of the present invention, there is provided a computer program product including a computer program which implements the method for determining the full-size distribution curve of the rock pore radius according to any embodiment of the present invention when executed by a processor.

[0021] In the technical solution of the embodiment of the present invention, by obtaining a rock sample and performing pretreatment on the rock sample to obtain a standard core, since the pretreatment at least includes oil washing treatment and / or drying treatment, impurities and moisture in the rock sample can be effectively removed, the credibility of experimental data can be improved, and reliable data support can be provided for subsequent processing of the rock sample; then, determining the gas logging pore information of the standard core, and obtaining the core scan data of the standard core through a computed tomography device, the internal structure and gas logging pore information of the rock sample can be obtained, and the accuracy and reliability of pore structure analysis can be improved; then, by processing the standard core into a water-saturated core and obtaining the transverse relaxation time spectrum of the water-saturated core through a magnetic resonance device, the distribution information of the fluid in the rock pores can be obtained, and the accuracy of pore structure and fluid property analysis can be improved; finally, by determining the rock pore radius distribution curve corresponding to the rock sample according to the core scan data, the gas logging pore information, and the transverse relaxation time spectrum, multiple data sources can be integrated, the accuracy of rock pore structure analysis can be improved, and the problem that it is difficult to reflect the pore distribution characteristics of the rock in the case of a complex pore structure in the related art is solved. It not only improves the accuracy and reliability of pore structure analysis, but also provides a scientific basis for optimizing the water flooding scheme and improving the recovery degree of crude oil in micro pores.

[0022] It should be understood that the content described in this part is not intended to identify the key or important features of the embodiments of the present invention, nor is it used to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention, and for those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0024] Figure 1 is a flowchart of a method for determining the full-size distribution curve of rock pore radius according to Embodiment 1 of the present invention;

[0025] Figure 2 is a flowchart of a method for determining the full-size distribution curve of rock pore radius according to Embodiment 2 of the present invention;

[0026] Figure 3 is a flowchart of a method for determining the full-size distribution curve of rock pore radius according to Embodiment 3 of the present invention;

[0027] Figure 4It is a schematic structural diagram of a device for determining the full-size distribution curve of rock pore radius according to Embodiment 3 of the present invention;

[0028] Figure 5 It is a schematic diagram of computer tomography equipment data processing and digital core provided by an embodiment of the present invention;

[0029] Figure 6 It is a schematic diagram of the CT pore radius distribution curve of the computer tomography equipment provided by an embodiment of the present invention;

[0030] Figure 7 It is a schematic diagram of the relaxation time and the cumulative distribution curve of pore radius provided by an embodiment of the present invention;

[0031] Figure 8 It is a schematic diagram of the conversion coefficient fitting curve provided by an embodiment of the present invention;

[0032] Figure 9 It is a schematic diagram of the converted nuclear magnetic resonance T2 map provided by an embodiment of the present invention;

[0033] Figure 10 It is a schematic diagram of the full-size pore radius distribution curve provided by an embodiment of the present invention;

[0034] Figure 11 It is a schematic diagram of the comparison between the mercury intrusion pore radius distribution and the full-size pore radius distribution provided by an embodiment of the present invention;

[0035] Figure 12 It is a schematic structural diagram of an electronic device for implementing the method for determining the full-size distribution curve of rock pore radius according to an embodiment of the present invention. Detailed implementation manners

[0036] In order to enable those skilled in the art to better understand the solution of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0037] It should be noted that the terms "first", "second", etc. in the description, claims and the above-mentioned drawings of the present invention are used to distinguish similar objects, and do not necessarily need to describe a specific order or sequence. It should be understood that such used data can be interchanged under appropriate circumstances, so that the embodiments of the present invention described here can be implemented in an order other than those illustrated or described here. In addition, the terms "comprising" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device comprising a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these processes, methods, products or devices.

[0038] It should be noted that the modifications of "one" and "multiple" mentioned in this disclosure are illustrative rather than restrictive. Those skilled in the art should understand that unless clearly specified otherwise in the context, it should be understood as "one or more".

[0039] The names of the messages or information exchanged between multiple devices in the embodiments of this disclosure are only for illustrative purposes and do not limit the scope of these messages or information.

[0040] It can be understood that before using the technical solutions disclosed in the embodiments of this disclosure, the types, usage scopes, usage scenarios, etc. of the personal information involved in this disclosure should be informed to users and user authorization should be obtained through appropriate means in accordance with relevant laws and regulations.

[0041] For example, when responding to receiving an active request from a user, a prompt message is sent to the user to clearly prompt the user that the operation requested by the user will require obtaining and using the user's personal information. Thus, the user can autonomously choose whether to provide personal information to software or hardware such as an electronic device, application program, server or storage medium that performs the operations of the technical solutions of this disclosure according to the prompt message.

[0042] As an optional but non-limiting implementation manner, the manner of sending a prompt message to the user in response to receiving an active request from the user can be, for example, in the form of a pop-up window, and the prompt message can be presented in text in the pop-up window. In addition, the pop-up window can also carry a selection control for the user to choose "agree" or "disagree" to provide personal information to the electronic device.

[0043] It can be understood that the above process of notifying and obtaining user authorization is only illustrative and does not limit the implementation manners of this disclosure. Other manners that meet relevant laws and regulations can also be applied to the implementation manners of this disclosure.

[0044] It is understandable that the data involved in this technical solution (including but not limited to the data itself, the acquisition or use of data) should comply with the requirements of corresponding laws, regulations and related provisions.

[0045] Embodiment 1

[0046] Figure 1 FIG. 1 is a flowchart of a method for determining a full-size distribution curve of rock pore radius provided by Embodiment 1 of the present invention. This embodiment is applicable to determining the pore structure of a rock sample. This method can be executed by a device for determining the full-size distribution curve of rock pore radius, and the device for determining the full-size distribution curve of rock pore radius can be implemented in the form of hardware and / or software. Optionally, it is implemented by an electronic device, and the electronic device can be a mobile terminal, a PC or a server, etc. As Figure 1 shown, the method may specifically include:

[0047] S110. Obtain a rock sample, and perform preprocessing on the rock sample to obtain a standard core, where the preprocessing at least includes oil washing treatment and / or drying treatment.

[0048] Among them, the rock sample can be understood as a small piece of rock collected from a rock reservoir for analysis. The rock sample can provide information on aspects such as formation composition, structure and physicochemical properties, and has a guiding role in fields such as oil exploration. Exemplarily, the core sample can be an artificial core or a natural core. The geometric parameters of the core sample can be a length of 5 cm, a diameter of 2.5 cm or 3.8 cm, and no specific limitation is made here. The preprocessing can be understood as a processing method adopted to improve the accuracy and reliability of the analysis results before specifically analyzing or detecting the rock sample, so as to eliminate factors that may interfere with the experimental results, such as pollutants and moisture. Exemplarily, the preprocessing at least can include oil washing treatment and / or drying treatment, etc. The oil washing treatment can be understood as an operation of removing oil stains or other organic substances in the rock sample by using a solvent for cleaning. It can be understood that some rock samples (especially samples from oil and gas fields) may contain oil stains or other organic substances, and the oil washing treatment can remove the surface oil stains, etc., to avoid affecting the subsequent analysis results. The drying treatment can be understood as an operation of heating the rock sample to remove the moisture in the rock sample. It can be understood that the presence of moisture in the rock sample may cause deviations in detection results (such as density and porosity). The standard core can be understood as a rock sample that meets certain specification requirements obtained after performing oil washing treatment and / or drying treatment on the rock sample. The standard core usually has unified characteristics such as size and shape, which is convenient for data comparison and communication between different laboratories. Exemplarily, the core sample can be subjected to oil washing and / or drying treatment according to a preset dry core treatment standard.

[0049] S120. Determine the gas logging pore information of the standard core, and obtain the core scan data of the standard core through a computed tomography (CT) device.

[0050] Among them, the computed tomography (CT) device can be understood as a technical device that uses a radiation source such as X-rays to scan an object from multiple angles and forms an image of the internal structure of the object through computer reconstruction technology. The CT scan can be used to non-destructively obtain detailed structural information inside a rock sample, such as pore distribution, etc. The core scan data can be understood as the digital information of the internal structure of the standard core obtained through the CT device. The core scan data can be a two-dimensional slice image or a three-dimensional reconstruction model, which can intuitively display the complex structural characteristics inside the standard core and provide an important basis for subsequent analysis. The gas logging pore information can be understood as the relevant parameters of the pore space in the measured rock sample, which is used to map the openness of the pores inside the standard rock and its influence on fluid transport performance. Exemplarily, the gas logging pore information may include, but is not limited to, parameters such as rock porosity and permeability, etc. The pore volume at different pore radii of the core can be calculated through the gas logging pore information of the standard core.

[0051] In an alternative implementation, the CT scan resolution of the core can reach the minimum resolution that the CT device can achieve under the standard core size, so as to ensure that the smallest pores of the standard core can be scanned as much as possible under the constraints of this CT device.

[0052] S130. Process the standard core into a water-saturated core, and obtain the transverse relaxation time spectrum of the water-saturated core through a magnetic resonance device.

[0053] In geology and petroleum engineering, processing the standard core into a water-saturated state can better simulate the conditions in the water flooding process in underground reservoirs. The water-saturated core can be understood as the standard core in a water-saturated state. The water-saturated state can be understood as a state in which all pores in the standard core are completely filled with water. Since the nuclear magnetic resonance signal mainly comes from hydrogen atoms in the fluid, the water-saturated core can be scanned by nuclear magnetic resonance. At this time, all the nuclear magnetic resonance signals come from the water in the pores with different radii of the core, so that the pore radius distribution of the core sample can be qualitatively described.

[0054] Exemplarily, the pores in the core can be completely filled with water by specific methods (such as vacuum suction or pressure saturation, etc.) to make the standard core in a water-saturated state. The Transverse Relaxation Time Spectrum (T2 spectrum) can be understood as a distribution map of the transverse relaxation time (T2) of the fluid in the water-saturated core obtained by a Nuclear Magnetic Resonance (NMR) device. Exemplarily, in rock sample analysis, the transverse relaxation time spectrum can provide information about the pore size, connectivity, and fluid type (such as the ratio of free water to bound water, etc.) of the rock. By analyzing the transverse relaxation time spectrum of the water-saturated core, the pore structure of the rock and its influence on fluid flow behavior can be deeply understood.

[0055] In an alternative embodiment, the process of treating the water-saturated core can perform the configuration of formation water and the operation of vacuum pumping and pressurizing the core with saturated water according to a preset standard for core saturation water treatment. Then, the experimental parameters of the nuclear magnetic resonance instrument are set for the water-saturated core, and the transverse relaxation time spectrum is measured.

[0056] S140. Determine the rock pore radius distribution curve corresponding to the rock sample according to the core scan data, the gas logging pore information, and the transverse relaxation time spectrum.

[0057] Among them, the rock pore radius distribution curve can be understood as a chart showing the distribution of pores of different sizes in the rock, used to reflect the quantity ratio of pores of different sizes in the rock sample. In the embodiments of the present disclosure, the pore radius cumulative distribution curve can be plotted with the pore radius and the cumulative frequency as the horizontal and vertical coordinates respectively. Among them, the cumulative frequency corresponding to the pore radius in the pore radius cumulative distribution curve can be understood as a parameter used to characterize the cumulative volume of pores less than and equal to this pore radius. For example, the cumulative frequency corresponding to the pore radius can specifically be the ratio of the cumulative volume of pores less than and equal to this pore radius to the total pore volume of the core sample. The pore volume can be determined based on the core porosity and the volume of the core.

[0058] Specifically, the first pore radius distribution curve can be determined according to the core scan data and the gas logging pore information, the second pore radius distribution curve can be determined according to the first pore radius distribution curve and the transverse relaxation time spectrum, and the rock pore radius distribution curve corresponding to the rock sample can be determined according to the first pore radius distribution curve and the second pore radius distribution curve.

[0059] Further, determining the first pore radius distribution curve according to the core scan data and the gas logging pore information may include: determining a scanned pore radius distribution curve according to the core scan data, and correcting the core scan data according to the gas logging pore information to obtain the first pore radius distribution curve of the rock sample.

[0060] The technical solution of the embodiment of the present invention obtains a rock sample and preprocesses the rock sample to obtain a standard core. Since the preprocessing at least includes oil washing treatment and / or drying treatment, impurities and moisture in the rock sample can be effectively removed, the credibility of experimental data can be improved, and reliable data support can be provided for subsequent processing of the rock sample. Then, the gas logging pore information of the standard core is determined, and the core scan data of the standard core is obtained by a computed tomography device, so that the internal structure and gas logging pore information of the rock sample can be obtained, and the accuracy and reliability of pore structure analysis can be improved. Then, by processing the standard core into a water-saturated core and obtaining the transverse relaxation time spectrum of the water-saturated core by a magnetic resonance device, the distribution information of the fluid in the rock pores can be obtained, and the accuracy of pore structure and fluid property analysis can be improved. Finally, by determining the rock pore radius distribution curve corresponding to the rock sample according to the core scan data, the gas logging pore information, and the transverse relaxation time spectrum, multiple data sources can be integrated, the accuracy of rock pore structure analysis can be improved, and the problem that it is difficult to reflect the pore distribution characteristics of the rock in the case of a complex pore structure in the related art is solved. It not only improves the accuracy and reliability of pore structure analysis, but also provides a scientific basis for optimizing the water flooding scheme and improving the production degree of crude oil in micro pores.

[0061] Embodiment 2

[0062] Figure 2The flowchart of a method for determining the full-size distribution curve of rock pore radius provided in the second embodiment of the present invention. The technical solution of this embodiment is a further refinement on the basis of the above embodiment, aiming at how to determine the rock pore radius distribution curve corresponding to the sample according to the core scanning data, the gas logging pore information, and the transverse relaxation time spectrum. Optionally, determining the rock pore radius distribution curve corresponding to the sample according to the core scanning data, the gas logging pore information, and the transverse relaxation time spectrum includes: determining the tortuosity, pore fractal dimension, and scanning pore radius distribution curve of the core sample according to the core scanning data, and determining the first pore radius distribution curve according to the scanning pore distribution curve and the gas logging pore information; converting the transverse relaxation time spectrum into a second pore radius distribution curve according to the tortuosity, the pore fractal dimension, and the first pore radius distribution curve; and determining the rock pore radius distribution curve corresponding to the rock sample according to the first pore radius curve and the second pore radius distribution curve. For the specific implementation manner, reference can be made to the description of this embodiment. Among them, the technical features that are the same as or similar to those of the foregoing embodiment will not be described again. As Figure 2 shown, the method may specifically include:

[0063] S210. Obtain a rock sample, and perform preprocessing on the rock sample to obtain a standard core, where the preprocessing includes at least oil washing treatment and / or drying treatment.

[0064] S220. Determine the gas logging pore information of the standard core, and obtain the core scanning data of the standard core through a computed tomography device.

[0065] S230. Process the standard core into a water-saturated core, and obtain the transverse relaxation time spectrum of the water-saturated core through a magnetic resonance device.

[0066] S240. Determine the tortuosity, pore fractal dimension, and scanning pore radius distribution curve of the core sample according to the core scanning data, and determine the first pore radius distribution curve according to the scanning pore distribution curve and the gas logging pore information.

[0067] Among them, the tortuosity can be understood as representing the degree of tortuosity of the flow path of the fluid in the pores of the rock sample. It is a dimensionless parameter used to quantify the complexity of the pore structure. Exemplarily, the greater the tortuosity, the more complex the pores of the rock sample, and the greater the resistance of the fluid passing through the pores of the rock sample. The pore fractal dimension can be understood as a parameter used to describe the complexity of the rock pore structure (such as roughness, etc.). The pore fractal dimension reflects the irregularity of the pore shape and distribution, usually between 2 and 3. It can be understood that a higher pore fractal dimension indicates a more complex and irregular pore structure. The scanned pore radius distribution curve can be understood as a curve of the pore radius distribution of the rock sample obtained based on CT scan data, used to reflect the proportion of the number of pores of different sizes inside the rock sample, and can be extracted from the CT scan image through image processing techniques. The first pore radius distribution curve can be understood as a distribution curve graph of the pore radius determined based on core scan data and gas logging pore information.

[0068] Based on the above solution, optionally, the determining the first pore radius distribution curve according to the scanned pore distribution curve and the gas logging hollow information includes: constructing a core three-dimensional model according to the core scan data, and determining the tortuosity, pore fractal dimension, and scanned pore radius distribution curve of the core sample according to the core three-dimensional model; correcting the scanned pore radius distribution curve according to the gas logging pore information to obtain the first pore radius distribution curve.

[0069] Among them, the core three-dimensional model can be understood as a three-dimensional geometric model of the rock sample constructed by computer modeling technology based on core scan data, used to visually display the pore structure and distribution inside the rock sample, and provide a visualization tool for analyzing the physical properties of the rock sample. The scanned pore radius distribution curve is used to characterize the mapping relationship between the pore radius and the cumulative frequency, and is used to reflect the proportion of the number of pores of different sizes inside the rock sample.

[0070] In order to improve the accuracy of the core scan data in subsequent data analysis, an optional implementation manner further includes: preprocessing the core scan data before constructing the core three-dimensional model according to the core scan data, where the preprocessing may include parameter setting, noise reduction, threshold segmentation, etc.

[0071] The above solution constructs a core three-dimensional model by combining core scan data, determines the tortuosity, pore fractal dimension, and scanned pore radius distribution curve, then corrects the scanned pore radius distribution curve, and finally obtains the first pore radius distribution curve, which not only improves the accuracy of pore structure analysis, but also can more truly reflect the internal characteristics of the rock sample, providing a reliable scientific basis for oil and gas resource development.

[0072] Due to the limited resolution of CT, it is impossible to calculate the pore volume corresponding to less than the minimum resolution. Without correction, the cumulative frequency in the cumulative distribution curve of pore radius starts from 0, without considering the influence of pores smaller than the minimum resolution. Nuclear magnetic resonance is a characterization of the total pore radius of the core. If the two curves are fitted, the error will be very large. Therefore, it is necessary to calculate the cumulative volume of pores smaller than the minimum resolution of the digital core through gas logging core porosity, and correct the pore radius distribution curve of CT.

[0073] On the basis of the above scheme, optionally, the correcting the scanned pore radius distribution curve according to the gas logging pore information to obtain the first pore radius distribution curve includes: determining the gas logging cumulative frequency corresponding to the core sample according to the gas logging pore information, and correcting the cumulative frequency corresponding to the pore radius in the scanned pore radius distribution curve according to the gas logging cumulative frequency to obtain the first pore radius distribution curve.

[0074] Among them, the gas logging cumulative frequency can be understood as the cumulative frequency distribution of pores with different pore sizes in the total pores determined according to the gas logging pore information. Exemplarily, the gas logging cumulative frequency can be less than or equal to the proportion of pores with a certain pore size in the total pores, and can be expressed as a percentage. The gas logging cumulative frequency can provide statistical data on the pore size distribution.

[0075] Correcting the cumulative frequency corresponding to the pore radius in the scanned pore radius distribution curve according to the gas logging cumulative frequency to obtain the first pore radius distribution curve. Specifically, for the same pore radius, the difference between the gas logging cumulative frequency corresponding to the pore radius and the cumulative frequency corresponding to the pore radius in the scanned pore radius distribution curve is calculated to obtain the frequency difference, that is, the cumulative frequency corresponding to the cumulative volume of pores that cannot be recognized by less than the CT minimum resolution. On this basis, the cumulative frequency in the scanned pore radius distribution curve is recalculated. For example, the cumulative frequency in the scanned pore radius distribution curve can be increased by the frequency difference to obtain the first pore radius distribution curve.

[0076] In practical applications, the first pore radius distribution curve and the second pore radius distribution curve are generally represented by a series of scatter points and are not continuous functions. Therefore, the cumulative frequency at any pore radius r p (i) of the first pore radius distribution curve is S(i), and the transverse relaxation time T 2 (i) of the second pore radius distribution curve corresponding to it may not have a value. At this time, interpolation can be used for calculation. Specifically, in the region greater than the CT minimum resolution, any pore radius r p (i) with a cumulative frequency of S(i) is selected, and S(i) is used to interpolate the T 2 spectrum to obtain the value corresponding to any pore radius r p(i) corresponding T 2 (i).

[0077] S250. Convert the transverse relaxation time spectrum into a second pore radius distribution curve according to the tortuosity, the pore fractal dimension, and the first pore radius distribution curve.

[0078] In an embodiment of the invention, specifically, according to the tortuosity of the three-dimensional core model, the pore fractal dimension, the pore radius in the first pore radius distribution curve, and a preset conversion function, determine a target conversion coefficient for converting the transverse relaxation time into a pore radius and a power exponent corresponding to the target conversion coefficient; substitute the target conversion coefficient and the power exponent into the preset conversion function to obtain a target conversion function, and determine the pore radius corresponding to the transverse relaxation time according to the target conversion function, the tortuosity, the fractal dimension, and the transverse relaxation time, so as to obtain a second pore radius distribution curve.

[0079] Among them, the preset conversion function can be understood as a function for converting the transverse relaxation time into a pore radius. The target conversion coefficient can be understood as the coefficient used in the preset conversion function to adjust the proportional relationship in the conversion process. The power exponent can be understood as the exponent used in the preset conversion function to reflect the non-linear relationship between the pore radius and the transverse relaxation time. The target conversion function can be understood as the conversion function obtained by substituting the target conversion coefficient and the power exponent into the preset conversion function. The target conversion function can be used to convert the transverse relaxation time into a pore radius. The second pore radius distribution curve can be understood as a pore radius distribution curve obtained by converting the transverse relaxation time spectrum, and the conversion process is determined according to the tortuosity, the pore fractal dimension, and the transverse relaxation time.

[0080] On the basis of the above solution, specifically, the preset conversion function is:

[0081]

[0082] Among them, r p is the pore radius; C is the conversion coefficient; τ is the tortuosity; D is the pore fractal dimension; n is the power exponent; T 2 is the transverse relaxation time.

[0083] On the basis of the above solution, specifically, the transverse relaxation time can be determined based on the following formula:

[0084]

[0085] Among them, D is the pore fractal dimension; τ is the tortuosity; r p is the pore radius; n is the power exponent; ρ 2is the transverse surface relaxation rate; F s is the pore shape factor.

[0086] Among them, the transverse surface relaxation rate can be understood as representing the transverse relaxation rate per unit pore surface area, which is used to reflect the interaction strength between the pore surface and fluid molecules and affects the relaxation time of T2. To simplify the usage method of the preset conversion function, the logarithm can be taken on both sides of the preset conversion function so as to calculate the values of the conversion coefficient C and the power exponent n by the least squares method. Specifically, the preset conversion function after taking the logarithm is:

[0087]

[0088] In the above solution, by combining the tortuosity, pore fractal dimension and the first pore radius distribution curve of the core three-dimensional model with the transverse relaxation time spectrum, and using the preset conversion function to determine the target conversion coefficient and power exponent corresponding to the core sample, the relaxation time is accurately converted into the pore radius, and the second pore radius distribution curve is generated, which improves the accuracy and reliability of the pore structure analysis and provides a scientific basis for optimizing the water flooding scheme and improving the production degree of crude oil in micro pores.

[0089] S260. Determine the rock pore radius distribution curve corresponding to the rock sample according to the first pore radius curve and the second pore radius distribution curve.

[0090] Based on the above solution, optionally, the determining the rock pore radius distribution curve corresponding to the rock sample according to the first pore radius curve and the second pore radius distribution curve includes: determining the critical pore radius according to the scanning resolution corresponding to the core scanning data, and splicing the curve segment of the first pore radius curve that is greater than or equal to the critical pore radius with the curve segment of the second pore radius curve that is less than the critical pore radius to obtain the rock pore radius distribution curve corresponding to the rock sample.

[0091] Among them, the critical pore radius can be understood as the pore radius threshold determined according to the scanning resolution of the core scanning data. The critical pore radius is used to distinguish pores of different sizes to ensure the accuracy of the pore radius distribution curve. Exemplarily, the minimum resolution that can be measured by a CT device can be used as the critical pore radius, such as 10 μm.

[0092] In the technical solution of the embodiment of the present invention, first, a three-dimensional core model is constructed based on core scanning data to determine the tortuosity, pore fractal dimension, and the distribution curve of scanned pore radius. Then, the gas logging pore information is used to correct the distribution curve of scanned pore radius to obtain the first pore radius distribution curve. Next, using a preset conversion function and nuclear magnetic resonance technology, the transverse relaxation time spectrum is converted into a second pore radius distribution curve. Finally, according to the resolution of the core scanning data, the critical pore radius is determined, and the first pore radius distribution curve and the second pore radius distribution curve are spliced to generate the pore radius distribution curve of the rock sample, which not only improves the accuracy and reliability of pore structure analysis but also provides a scientific basis for optimizing the water flooding scheme and improving the production degree of crude oil in micro pores.

[0093] Embodiment III

[0094] Figure 3 FIG. is a flowchart of a method for determining the full-size distribution curve of rock pore radius provided by Embodiment III of the present invention, which is an alternative embodiment of the above embodiment. The method includes:

[0095] First, the rock sample is subjected to oil washing treatment and / or drying treatment to obtain a standard core.

[0096] Then, the treated core is placed in a CT scanning device for scanning to obtain core scanning data, which needs to reach the minimum resolution that the CT device can achieve under the size of this kind of core to ensure that the smallest pores of the core can be scanned as much as possible under the constraint of this kind of CT device.

[0097] Next, the core scanning data is imported into relevant software to construct a digital core. The reconstruction of the digital core needs to go through steps such as parameter setting, noise reduction, threshold segmentation (dividing pores and skeletons), and three-dimensional reconstruction, and finally form a three-dimensional core model, as Figure 5 shown, and determine the tortuosity, fractal dimension, and the cumulative distribution curve of core pore radius of the three-dimensional core model. The CT pore radius distribution curve is as Figure 6 shown.

[0098] Then, simulate the formation water according to the formation water ion composition table of the target oilfield, and then perform a vacuum saturation water operation on the dry core. Place the core saturated with formation water in a nuclear magnetic resonance device to measure the nuclear magnetic resonance T2 spectrum of the core in the saturated water state.

[0099] Then, solve the conversion coefficient using the CT pore radius distribution curve and the nuclear magnetic resonance T2 spectrum. The specific operation is as follows:

[0100] 1) Plot the cumulative distribution curves of CT pore radius and nuclear magnetic resonance T2 relaxation time. The results are as Figure 7As shown, due to the limitation of the resolution of CT equipment, pores smaller than the minimum resolution cannot be identified, and it is necessary to correct the CT pore radius cumulative distribution curve (the first pore radius distribution curve) by combining core gas logging porosity data;

[0101] 2) In the region larger than the CT minimum resolution, the cumulative frequency at pore radius r p (i) is S(i), and the T 2 spectrum is interpolated using S(i) to obtain T p corresponding to any pore radius r 2 (i);

[0102] 3) On the basis of the original method of converting mercury injection and nuclear magnetic resonance by non-linear power function, quantitative indexes for describing the irregularity (tortuosity τ) and roughness (fractal dimension D) of cores (porous media) in the petroleum and natural gas industries are introduced to correct the conversion relationship:

[0103]

[0104] Let C = (ρ 2 F s ) 1n , then formula (3) can be transformed into:

[0105]

[0106] Taking the natural logarithm of both sides of formula (4), we can get:

[0107]

[0108] Among them, C = (ρ 2 F s ) 1n ; r p is the pore radius; C is the conversion coefficient; τ is the tortuosity; D is the pore fractal dimension; n is the power exponent; T 2 is the transverse relaxation time; F s is the pore shape factor; ρ 2 is the transverse surface relaxation rate.

[0109] According to the least square method principle, parameter fitting calculation is carried out on formula (5) to obtain the C and n values with the minimum error, where the fractal dimension and tortuosity have been obtained in the above steps. As Figure 8 shown, the relaxation time in the T2 spectrum can be converted into the pore radius distribution. From the fitting results, the test data and the fitting line are in good agreement in morphology and the fitting accuracy is high, which proves the rationality of the fitting method. Substituting the solved C and n values into formula (5), the relaxation time in the T 2 spectrum can be converted into the pore radius distribution, that is, the converted nuclear magnetic T2 map, asFigure 9 as shown

[0110] Finally, taking the 10-μm resolution of the CT device as the demarcation point, the CT pore radius distribution curve is then spliced together with the pore radius distribution curve less than 10 μm after the NMR conversion, and at the same time, the repeated pore radius distributions near the junction point are removed. Finally, the full-size distribution curve of the rock pore structure is obtained, realizing the quantitative characterization of the full-size distribution of the reservoir rock pore radius, as Figure 10 shown

[0111] In order to verify the accuracy of the NMR T2 spectrum conversion and the full-size characterization method of the rock pore structure proposed in this experiment, the present invention uses the same core for mercury injection experiment, and compares the pore radius distribution curve of the mercury injection core with the full-size distribution curve of the rock pore radius in this embodiment in the same coordinate system, as Figure 11 shown. It can be clearly obtained that the trends of the core pore radius distribution curves measured by the two different methods are the same, which proves the accuracy of the full-size characterization method of the rock pore structure proposed in the present invention. It can characterize the micro-pores and even nano-pores that cannot be measured by the mercury injection method, and has strong guiding significance for the development of complex core carbonate reservoirs.

[0112] The technical solution of the embodiment of the present invention is as follows: First, the rock sample is washed with oil and dried according to the national standard to obtain a standard core. Then, the internal structure data of the core is obtained by using a high-resolution CT scanning device, and a three-dimensional model of the core is constructed to determine the tortuosity, fractal dimension and pore radius cumulative distribution curve. Next, the T 2 relaxation time spectrum of the core in the saturated water state is measured by nuclear magnetic resonance technology, and combined with the CT pore radius distribution curve and the gas logging porosity data, the conversion coefficient is solved by least square fitting to convert the T 2 relaxation time into the pore radius distribution. Finally, taking the minimum resolution (10 μm) of the CT device as the demarcation point, the CT pore radius distribution curve is spliced with the pore radius distribution curve after the NMR conversion, and the repeated part is removed to generate the full-size pore radius distribution curve of the rock sample, which not only improves the accuracy and reliability of the pore structure analysis, but also provides a scientific basis for optimizing the water flooding scheme and improving the production degree of crude oil in micro-pores, significantly improving the development efficiency and recovery rate of carbonate reservoirs.

[0113] Embodiment 4

[0114] Figure 4 is a schematic structural diagram of a device for determining the full-size distribution curve of the rock pore radius provided in Embodiment 3 of the present invention. As Figure 4 shown, the device includes: a standard core acquisition module 410, a standard core analysis module 420, a saturated water core analysis module 430, and a pore radius distribution determination module 440.

[0115] Among them, the standard core acquisition module 410 is used to obtain a rock sample and preprocess the rock sample to obtain a standard core. Among them, the preprocessing at least includes oil washing treatment and / or drying treatment; the standard core analysis module 420 is used to determine the gas logging pore information of the standard core, and obtain the core scan data of the standard core through a computed tomography device; the saturated water core analysis module 430 is used to process the standard core into a saturated water core and obtain the transverse relaxation time spectrum of the saturated water core through a magnetic resonance device; the pore radius distribution determination module 440 is used to determine the rock pore radius distribution curve corresponding to the rock sample according to the core scan data, the gas logging pore information, and the transverse relaxation time spectrum.

[0116] In the technical solution of the embodiment of the present invention, the standard core acquisition module 410 is used to obtain a rock sample and preprocess the rock sample to obtain a standard core. Since the preprocessing at least includes oil washing treatment and / or drying treatment, impurities and moisture in the rock sample can be effectively removed, the credibility of experimental data can be improved, and reliable data support can be provided for subsequent processing of the rock sample; then, the standard core analysis module 420 is used to determine the gas logging pore information of the standard core, and the core scan data of the standard core is obtained through a computed tomography device, so that the internal structure and gas logging pore information of the rock sample can be obtained, and the accuracy and reliability of pore structure analysis can be improved; then, the saturated water core analysis module 430 is used to process the standard core into a saturated water core and obtain the transverse relaxation time spectrum of the saturated water core through a magnetic resonance device, so that the distribution information of the fluid in the rock pores can be obtained, and the accuracy of pore structure and fluid property analysis can be improved; finally, the pore radius distribution determination module 440 is used to determine the rock pore radius distribution curve corresponding to the rock sample according to the core scan data, the gas logging pore information, and the transverse relaxation time spectrum, which can integrate multiple data sources, improve the accuracy of rock pore structure analysis, solve the problem in the related art that it is difficult to reflect the pore distribution characteristics of the rock in the case of complex pore structures, not only improve the accuracy and reliability of pore structure analysis, but also provide a scientific basis for optimizing the water flooding scheme and improving the recovery degree of crude oil in micro pores.

[0117] On the basis of the above solution, optionally, the pore radius distribution determination module may include: a first pore radius distribution curve determination sub-module, a second pore radius distribution curve determination sub-module, and a rock pore radius distribution curve determination sub-module.

[0118] Among them, the first pore radius distribution curve determination sub-module is used to determine the tortuosity, pore fractal dimension, and scanning pore radius distribution curve of the core sample according to the core scanning data, and determine the first pore radius distribution curve according to the scanning pore distribution curve and the gas logging hollow information; the second pore radius distribution curve determination sub-module is used to convert the transverse relaxation time spectrum into the second pore radius distribution curve according to the tortuosity, the pore fractal dimension, and the first pore radius distribution curve; the rock pore radius distribution curve determination sub-module is used to determine the rock pore radius distribution curve corresponding to the rock sample according to the first pore radius curve and the second pore radius distribution curve.

[0119] On the basis of the above solution, optionally, the first pore radius distribution curve determination sub-module may include: a scanning pore radius distribution curve determination unit, and a first pore radius distribution curve acquisition unit.

[0120] Among them, the scanning pore radius distribution curve determination unit is used to construct a core three-dimensional model according to the core scanning data, and determine the tortuosity, pore fractal dimension, and scanning pore radius distribution curve of the core sample according to the core three-dimensional model, where the scanning pore radius distribution curve is used to characterize the mapping relationship between the pore radius and the cumulative frequency; the first pore radius distribution curve acquisition unit is used to correct the scanning pore radius distribution curve according to the gas logging pore information to obtain the first pore radius distribution curve.

[0121] On the basis of the above solution, optionally, the first pore radius distribution curve acquisition unit is specifically used to: determine the gas logging cumulative frequency corresponding to the core sample according to the gas logging pore information, and correct the cumulative frequency corresponding to the pore radius in the scanning pore radius distribution curve according to the gas logging cumulative frequency to obtain the first pore radius distribution curve.

[0122] On the basis of the above solution, optionally, the second pore radius distribution curve determination sub-module may include: an exponential power determination unit, and a second pore radius distribution curve acquisition unit.

[0123] Among them, the power exponent determination unit is configured to determine a target conversion coefficient for converting the transverse relaxation time into a pore radius and a power exponent corresponding to the target conversion coefficient according to the tortuosity of the three-dimensional core model, the pore fractal dimension, the pore radius in the first pore radius distribution curve, and a preset conversion function; the second pore radius distribution curve acquisition unit is configured to substitute the target conversion coefficient and the power exponent into the preset conversion function to obtain a target conversion function, and determine a pore radius corresponding to the transverse relaxation time according to the target conversion function, the tortuosity, the fractal dimension, and the transverse relaxation time, so as to obtain a second pore radius distribution curve.

[0124] Based on the above solution, optionally, the preset conversion function is:

[0125]

[0126] where r p is the pore radius; C is the conversion coefficient; τ is the tortuosity; D is the pore fractal dimension; n is the power exponent; T 2 is the transverse relaxation time.

[0127] Based on the above solution, optionally, the rock pore radius distribution curve determination sub-module is specifically configured to: determine a critical pore radius according to the scanning resolution corresponding to the core scanning data, and splice a curve segment in the first pore radius curve that is greater than or equal to the critical pore radius with a curve segment in the second pore radius curve that is less than the critical pore radius to obtain the rock pore radius distribution curve corresponding to the rock sample.

[0128] The device for determining the full-size distribution curve of the rock pore radius provided by the embodiments of the present invention can execute the method for determining the full-size distribution curve of the rock pore radius provided by any embodiment of the present invention, and has corresponding functional modules and beneficial effects for executing the method.

[0129] Embodiment 5

[0130] Figure 12 FIG. shows a schematic structural diagram of an electronic device 10 that can be used to implement the embodiments of the present invention. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workbenches, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processors, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are only examples and are not intended to limit the implementation of the present invention described and / or claimed herein.

[0131] As shown Figure 12 in FIG. 1, the electronic device 10 includes at least one processor 11 and a memory communicatively connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc. The memory stores a computer program executable by the at least one processor. The processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 into the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0132] A plurality of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a magnetic disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0133] The processor 11 can be various general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of the processor 11 include but are not limited to a central processing unit (CPU), a graphics processing unit (GPU), various dedicated artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as a method for determining the full-size distribution curve of rock pore radii.

[0134] In some embodiments, a method for determining the full-size distribution curve of rock pore radii can be implemented as a computer program tangibly embodied in a computer-readable storage medium, such as the storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed onto the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the method for determining the full-size distribution curve of rock pore radii described above can be executed. Alternatively, in other embodiments, the processor 11 can be configured to execute a method for determining the full-size distribution curve of rock pore radii by any other suitable means (e.g., by means of firmware).

[0135] The various embodiments of the systems and techniques described above in this specification can be implemented in digital electronic circuitry, integrated circuit systems, field programmable gate arrays (FPGA), application specific integrated circuits (ASIC), application specific standard products (ASSP), systems on a chip (SOC), complex programmable logic devices (CPLD), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include: being implemented in one or more computer programs that are executable and / or interpretable on a programmable system including at least one programmable processor, which can be a special-purpose or general-purpose programmable processor that receives data and instructions from, and transmits data and instructions to, a storage system, at least one input device, and at least one output device.

[0136] The computer programs for implementing the methods of the present invention can be written in any combination of one or more programming languages. These computer programs can be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus, such that the computer programs, when executed by the processor, cause the functions / operations specified in the flowchart and / or block diagram to be implemented. The computer programs can be executed entirely on the machine, partly on the machine, as a stand-alone software package partly on the machine and partly on a remote machine or entirely on the remote machine or server.

[0137] In the context of the present invention, a computer-readable storage medium can be a tangible medium that can contain, or store a computer program for use by or in connection with an instruction execution system, apparatus, or device. The computer-readable storage medium can include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. Alternatively, the computer-readable storage medium can be a machine-readable signal medium. More specific examples of the machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0138] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or a trackball) through which the user can provide input to the electronic device. Other kinds of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0139] The systems and techniques described herein can be implemented in a computing system including backend components (e.g., as a data server), or a computing system including middleware components (e.g., an application server), or a computing system including frontend components (e.g., a user computer having a graphical user interface or a web browser through which the user can interact with an implementation of the systems and techniques described herein), or a computing system including any combination of such backend components, middleware components, or frontend components. The components of the system can be interconnected to each other by digital data communication in any form or medium (e.g., a communication network). Examples of communication networks include: local area network (LAN), wide area network (WAN), blockchain network, and the Internet.

[0140] A computing system can include a client and a server. The client and the server are generally far from each other and usually interact through a communication network. The client-server relationship is created by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or a cloud host, which is a host product in the cloud computing service system and solves the defects of difficult management and weak business scalability existing in traditional physical hosts and VPS services.

[0141] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps recited in the present invention can be executed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of the present invention can be achieved, and no limitation is made herein.

[0142] The above specific embodiments do not constitute a limitation on the protection scope of the present invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims

1. A method for determining a full-size distribution curve of rock pore radius, characterized in that: include: Obtaining a rock sample, and pre-treating the rock sample to obtain a standard core, wherein the pre-treatment at least includes oil washing treatment and / or drying treatment; Determine the gas logging pore information of the standard core, and obtain core scanning data of the standard core by means of a computer tomography device; Processing the standard core into a water-saturated core, and obtaining a transverse relaxation time spectrum of the water-saturated core by a magnetic resonance device; A rock pore radius distribution curve corresponding to the rock sample is determined according to the core scanning data, the gas logging pore information and the transverse relaxation time spectrum.

2. The method according to claim 1, characterized in that The step of determining the rock pore radius distribution curve corresponding to the rock sample according to the core scanning data, the gas logging pore information and the transverse relaxation time spectrum includes: Determine the tortuosity, pore fractal dimension and scanning pore radius distribution curve of the core sample according to the core scanning data, and determine a first pore radius distribution curve according to the scanning pore distribution curve and the gas logging pore information; converting the transverse relaxation time spectrum into a second pore radius distribution curve according to the tortuosity, the pore fractal dimension and the first pore radius distribution curve; A rock pore radius distribution curve corresponding to the rock sample is determined according to the first pore radius curve and the second pore radius distribution curve.

3. The method according to claim 2, characterized in that The determining of a first pore radius distribution curve according to the scanning pore distribution curve and the gas measurement pore information comprises: Constructing a three-dimensional core model according to the core scanning data, and determining the tortuosity, pore fractal dimension and scanning pore radius distribution curve of the core sample according to the three-dimensional core model, wherein the scanning pore radius distribution curve is used to characterize the mapping relationship between pore radius and cumulative frequency; The scanning pore radius distribution curve is corrected according to the gas measurement pore information to obtain a first pore radius distribution curve.

4. The method according to claim 3, characterized in that The step of correcting the scanned pore radius distribution curve according to the gas-measured pore information to obtain a first pore radius distribution curve includes: The gas logging cumulative frequency corresponding to the core sample is determined according to the gas logging pore information, and the cumulative frequency corresponding to the pore radius in the scanning pore radius distribution curve is corrected according to the gas logging cumulative frequency to obtain a first pore radius distribution curve.

5. The method according to claim 2, characterized in that: The converting the transverse relaxation time spectrum into a second pore radius distribution curve according to the tortuosity, the pore fractal dimension and the first pore radius distribution curve comprises: Determine a target conversion coefficient for converting the transverse relaxation time into the pore radius and a power exponent corresponding to the target conversion coefficient according to the tortuosity of the three-dimensional core model, the pore fractal dimension, the pore radius in the first pore radius distribution curve, and a preset conversion function; The target conversion coefficient and the power exponent are substituted into the preset conversion function to obtain a target conversion function, and the pore radius corresponding to the transverse relaxation time is determined according to the target conversion function, the tortuosity, the fractal dimension and the transverse relaxation time to obtain a second pore radius distribution curve.

6. The method according to claim 5, characterized in that The preset conversion function is: Among them, r p is the pore radius; C is the conversion coefficient; τ is the tortuosity; D is the pore fractal dimension; n is the power exponent; T2 is the transverse relaxation time.

7. The method according to claim 2, characterized in that The step of determining a rock pore radius distribution curve corresponding to the rock sample according to the first pore radius curve and the second pore radius distribution curve includes: The critical pore radius is determined according to the scanning resolution corresponding to the core scanning data, and the curve segment in the first pore radius curve that is greater than or equal to the critical pore radius is spliced ​​with the curve segment in the second pore radius curve that is less than the critical pore radius to obtain a rock pore radius distribution curve corresponding to the rock sample.

8. A device for determining a full-size distribution curve of rock pore radius, characterized in that: include: A standard core acquisition module is used to acquire rock samples and pre-treat the rock samples to obtain standard cores, wherein the pre-treatment at least includes oil washing treatment and / or drying treatment; A standard core analysis module, used to determine the gas logging pore information of the standard core, and to obtain core scanning data of the standard core through a computer tomography device; A water-saturated core analysis module, used for processing the standard core into a water-saturated core, and obtaining a transverse relaxation time spectrum of the water-saturated core through a magnetic resonance device; The pore radius distribution determination module is used to determine the rock pore radius distribution curve corresponding to the rock sample according to the core scanning data, the gas logging pore information and the transverse relaxation time spectrum.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for determining the full-size distribution curve of rock pore radius as described in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for determining the rock pore radius full-size distribution curve according to any one of claims 1 to 7 when executed.

Images