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

Abstract

The application discloses a method and device for determining a full-size distribution curve of rock pore radius, and the method comprises the following steps: obtaining a rock sample, pretreating the rock sample to obtain a standard core, wherein the pretreatment at least comprises oil washing treatment and / or drying treatment; determining gas-measured pore information of the standard core, and obtaining core scanning data of the standard core through a computer tomography device; treating 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; and determining a rock pore radius distribution curve corresponding to the rock sample according to the core scanning data, the gas-measured pore information and the transverse relaxation time spectrum. The method not only improves the accuracy and reliability of pore structure analysis, but also provides a scientific basis for optimizing water flooding schemes and improving the producing degree of crude oil in micro pores.

Description

TECHNICAL FIELD

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

[0002] In the field of geology and petroleum engineering, accurately characterizing the pore structure of carbonate reservoirs, especially micron-sized and smaller pores, is crucial for improving waterflood recovery and oil production.

[0003] In related technologies, pore structure analysis methods such as mercury injection and gas adsorption can provide some pore information, but have problems such as complex operation, sample destruction, and limited data accuracy. Especially under complex pore structure conditions, it is difficult to comprehensively and accurately reflect the pore distribution characteristics of rocks. SUMMARY

[0004] The present application provides a method and device for determining a full-size distribution curve of rock pore radius to solve the problem that related technologies cannot reflect the pore distribution characteristics of rocks under complex pore structure conditions.

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

[0006] Obtaining a rock sample, pre-treating the rock sample to obtain a standard core, wherein the pre-treatment at least includes oil washing treatment and / or drying treatment;

[0007] Determining gas logging pore information of the standard core, and obtaining core scanning data of the standard core through a computer tomography device;

[0008] 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;

[0009] Determining a 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.

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

[0011] A standard core obtaining module is configured to obtain a rock sample, pre-treat the rock sample to obtain a standard core, wherein the pre-treatment at least includes oil washing treatment and / or drying treatment;

[0012] A standard core analysis module is configured to determine gas logging pore information of the standard core, and obtain core scanning data of the standard core through a computer tomography device;

[0013] a saturated water core analysis module configured to process the standard core into a saturated water core and acquire a transverse relaxation time profile of the saturated water core by a magnetic resonance device;

[0014] a pore radius distribution determination module configured to determine a rock pore radius distribution curve corresponding to the rock sample according to the core scanning data, the gas logging porosity information and the transverse relaxation time profile.

[0015] According to another aspect of the present application, there is provided an electronic device, comprising:

[0016] at least one processor; and

[0017] a memory connected to the at least one processor in communication; 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 to enable the at least one processor to perform the method for determining a rock pore radius full-size distribution curve according to any one of the embodiments of the present application.

[0019] According to another aspect of the present application, there is provided a computer readable storage medium storing computer instructions for enabling a processor to perform the method for determining a rock pore radius full-size distribution curve according to any one of the embodiments of the present application when executed by the processor.

[0020] According to another aspect of the present application, there is provided a computer program product comprising a computer program for enabling a processor to perform the method for determining a rock pore radius full-size distribution curve according to any one of the embodiments of the present application when executed by the processor.

[0021] The technical scheme of the embodiment of the present application is characterized in that a rock sample is obtained, and the rock sample is pretreated to obtain a standard core, and the pretreatment at least includes oil washing treatment and / or drying treatment, so that impurities and moisture in the rock sample can be effectively removed, the reliability of experimental data is improved, and reliable data support is provided for subsequent processing of the rock sample; then, the gas logging pore information of the standard core is determined, and the core scanning data of the standard core is obtained by a computer tomography device, so that the internal structure and the gas logging pore information of the rock sample can be obtained, and the accuracy and reliability of the pore structure analysis are improved; then, the standard core is processed into a saturated water core, and the transverse relaxation time spectrum of the saturated water core is obtained by a magnetic resonance device, so that the distribution information of fluid in the rock pore can be obtained, and the accuracy of the pore structure and fluid property analysis is improved; finally, the 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, so that multiple data sources can be comprehensively used, the accuracy of the rock pore structure analysis is improved, and the problem that the pore distribution characteristics of the rock are difficult to reflect in the related art under the condition of a complex pore structure is solved, and not only the accuracy and reliability of the pore structure analysis are improved, but also a scientific basis is provided for optimizing a water drive scheme and improving the producing degree of crude oil in a small pore.

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

[0023] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed in the embodiment description will be briefly introduced as follows. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0024] Figure 1 is a flow chart of a method for determining a rock pore radius full-size distribution curve according to the first embodiment of the present application;

[0025] Figure 2 is a flow chart of a method for determining a rock pore radius full-size distribution curve according to the second embodiment of the present application;

[0026] Figure 3 is a flow chart of a method for determining a rock pore radius full-size distribution curve according to the third embodiment of the present application;

[0027] Figure 4It is a structural schematic view of a rock pore radius full-size distribution curve determination device provided according to the third embodiment of the present application.

[0028] Figure 5 A computer tomography device data processing and digital core schematic view is provided for the embodiment of the present application.

[0029] Figure 6 A computer tomography device CT pore radius distribution curve schematic view is provided for the embodiment of the present application.

[0030] Figure 7 A relaxation time and pore radius cumulative distribution curve schematic view is provided for the embodiment of the present application.

[0031] Figure 8 A conversion coefficient fitting curve schematic view is provided for the embodiment of the present application.

[0032] Figure 9 A converted nuclear magnetic resonance T2 spectrum schematic view is provided for the embodiment of the present application.

[0033] Figure 10 A full-size pore radius distribution curve schematic view is provided for the embodiment of the present application.

[0034] Figure 11 A mercury injection pore radius distribution and full-size pore radius distribution comparison schematic view is provided for the embodiment of the present application.

[0035] Figure 12 It is a structural schematic view of an electronic device for implementing a rock pore radius full-size distribution curve determination method. DETAILED DESCRIPTION

[0036] In order to make the personnel in the art better understand the present application scheme, the technical solutions in the embodiments of the present application will be described clearly and completely in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by the person skilled in the art without creative labor should belong to the scope of protection of the present application.

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

[0038] It should be noted that the "one" or "multiple" modification mentioned in the present disclosure is illustrative but not limiting, and those skilled in the art should understand that unless the context clearly indicates otherwise, it should be understood as "one or more".

[0039] The names of the messages or information exchanged between the plurality of devices in the embodiments of the present disclosure are only for illustrative purposes, and are not intended to limit the scope of the messages or information.

[0040] It can be understood that before using the technical solutions disclosed in the embodiments of the present disclosure, the type, use range, use scenario, etc. of the personal information involved in the present disclosure should be informed to the user and the authorization of the user should be obtained through appropriate means according to relevant laws and regulations.

[0041] For example, in response to receiving the active request of the user, the user is sent prompt information to explicitly prompt the user that the operation requested to be performed will require obtaining and using the personal information of the user. Thus, the user can voluntarily choose whether to provide personal information to the software or hardware such as electronic devices, application programs, servers or storage media that perform the operation of the technical solutions of the present disclosure according to the prompt information.

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

[0043] It can be understood that the above notification and user authorization process is only illustrative and does not limit the implementation manner of the present disclosure, and other manners that meet the relevant laws and regulations can also be applied to the implementation manner of the present disclosure.

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

[0045] Embodiment one

[0046] Figure 1 A flowchart of a method for determining a rock pore radius full-size distribution curve is provided for the first embodiment of the present application. The present embodiment can be applicable to determining the pore structure of a rock sample. The method can be executed by a rock pore radius full-size distribution curve determination device, which can be realized in the form of hardware and / or software. Optionally, the device can be realized by an electronic device, which can be a mobile terminal, a PC terminal, a server, or the like. As shown in FIG. 1, the method can specifically include the following steps. Figure 1

[0047] S110, 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.

[0048] 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 about the composition, structure, and physical and chemical properties of the formation, which can guide the field of oil exploration. The core sample can be a man-made core or a natural core. The geometric parameters of the core sample can be a length of 5 cm and a diameter of 2.5 cm or 3.8 cm, which is not limited herein. The pre-treatment can be understood as a treatment method taken before a specific analysis or detection of the rock sample to improve the accuracy and reliability of the analysis results, so as to remove factors that may interfere with the experimental results, such as pollutants and moisture. The pre-treatment can at least include oil washing treatment and / or drying treatment. The oil washing treatment can be understood as an operation of removing oil stains or other organic matter 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) can contain oil stains or other organic matter, and the oil washing treatment can remove the oil stains on the surface to avoid affecting the subsequent analysis results. The drying treatment can be understood as an operation of heating the rock sample to remove moisture in the rock sample. It can be understood that the moisture in the rock sample may cause deviations in the detection results (such as density and porosity). The standard core can be understood as a rock sample that meets certain specification requirements after being treated by oil washing and / or drying. The standard core usually has uniform size, shape, and other characteristics, which facilitates the comparison and exchange of data between different laboratories. The core sample can be treated according to the preset dry core treatment standard.​

[0049] S120. Determine the gas porosity information of the standard core and acquire the core scanning data of the standard core using a computed tomography (CT) scanner.

[0050] The computed tomography (CT) scanner can be understood as a technical device that uses X-rays or other radiation sources to scan an object from multiple angles and then uses computer reconstruction technology to form an image of the object's internal structure. The CT scan can be used to non-destructively obtain detailed structural information about the interior of a rock sample, such as pore distribution. The core scan data can be understood as digital information about the internal structure of a standard rock core obtained through CT scanning. The core scan data can be two-dimensional slice images or three-dimensional reconstructed models, which can intuitively display the complex structural features inside the standard rock core, providing important basis for subsequent analysis. The gas-measured porosity information can be understood as the relevant parameters of the pore space in the rock sample, used to map the degree of openness of the pores inside the standard rock and its impact on fluid transport performance. For example, the gas-measured porosity information may include, but is not limited to, parameters such as rock porosity and permeability. The pore volume at different pore radii can be calculated using the gas-measured porosity information of the standard rock core.

[0051] In one alternative implementation, the core CT scan resolution can reach the minimum resolution achievable by the CT equipment for a standard core size, so as to ensure that the minimum porosity of the standard core under the constraints of the CT equipment can be scanned as much as possible.

[0052] S130. The standard core is processed into a water-saturated core, and the transverse relaxation time spectrum of the water-saturated core is obtained by magnetic resonance equipment.

[0053] In geology and petroleum engineering, treating standard core samples to a water-saturated state allows for better simulation of water-driven oil recovery processes in underground reservoirs. A water-saturated core can be understood as a standard core in a state of complete water saturation, where all pores are completely filled with water. Since nuclear magnetic resonance (NMR) signals primarily originate from hydrogen atoms in the fluid, NMR scans can be performed on water-saturated cores. In this case, the NMR signals originate entirely from water within different pore radii of the core sample, thus enabling a qualitative description of the pore radius distribution.

[0054] Exemplarily, the standard core can be saturated with water by a specific method (e.g., vacuum suction or pressure saturation, etc.) to make the standard core in a saturated water state. The transverse relaxation time spectrum (T2 spectrum) can be understood as a transverse relaxation time (T2) distribution diagram of fluid in the saturated water core obtained by a nuclear magnetic resonance (NMR) device. Exemplarily, in rock sample analysis, the transverse relaxation time spectrum can provide information about the rock pore size, connectivity, and fluid type (e.g., the proportion of free water and bound water, etc.). By analyzing the transverse relaxation time spectrum of the saturated water core, the pore structure of the rock and its influence on fluid flow behavior can be understood in depth.

[0055] In an alternative embodiment, the process of the saturated water core can be carried out according to the preset core saturated water treatment standard to configure the formation water and the core vacuum pumping and pressurized saturated water operation. Then, the saturated water core is subjected to the NMR instrument experimental parameter setting and the transverse relaxation time spectrum measurement.

[0056] S140, 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.

[0057] The rock pore radius distribution curve can be understood as a graph showing the distribution of pores of different sizes in the rock, which is used to reflect the number proportion 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. The cumulative frequency corresponding to the pore radius in the pore radius cumulative distribution curve can be understood as a parameter for characterizing the cumulative volume of pores smaller than and equal to the pore radius, for example, the cumulative frequency corresponding to the pore radius can be the ratio of the cumulative volume of pores smaller than and equal to the 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 scanning 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 scanning data and the gas logging pore information can include: determining a scanning pore radius distribution curve according to the core scanning data, and correcting the core scanning data according to the gas logging pore information to obtain the first pore radius distribution curve of the rock sample.

[0060] The technical scheme of the embodiment of the present application, by obtaining a rock sample, pre-treating the rock sample to obtain a standard core, since the pre-treatment at least includes oil washing treatment and / or drying treatment, can effectively remove impurities and moisture in the rock sample, improve the reliability of experimental data, and provide reliable data support for subsequent processing of the rock sample; then, determining the gas logging pore information of the standard core, and obtaining the core scanning data of the standard core by a computer tomography device, can obtain the internal structure and gas logging pore information of the rock sample, improve the accuracy and reliability of the pore structure analysis; then, by treating the standard core into a saturated water core and obtaining the transverse relaxation time spectrum of the saturated water core by a magnetic resonance device, the distribution information of fluid in the rock pore can be obtained, and the accuracy of the pore structure and fluid property analysis is improved; finally, by 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, multiple data sources can be integrated, the accuracy of the rock pore structure analysis is improved, and the problem that in the related art, under the condition of complex pore structure, it is difficult to reflect the pore distribution characteristics of the rock is solved, not only the accuracy and reliability of the pore structure analysis are improved, but also a scientific basis for optimizing the water drive scheme and improving the producing degree of crude oil in the micro-pore is provided.

[0061] Embodiment two

[0062] Figure 2A flowchart of a method for determining a full-size distribution curve of a rock pore radius is provided for the second embodiment of the present application. The technical solution of the present embodiment is a further refinement of the method for determining the rock pore radius distribution curve corresponding to the sample pattern based on the core scanning data, the gas-measured pore information, and the transverse relaxation time map, based on the above-mentioned embodiment. Alternatively, the method for determining the rock pore radius distribution curve corresponding to the sample pattern based on the core scanning data, the gas-measured pore information, and the transverse relaxation time map comprises: determining the tortuosity, pore fractal dimension, and scanning pore radius distribution curve of the core sample based on the core scanning data, determining a first pore radius distribution curve based on the scanning pore distribution curve and the gas-measured pore information, converting the transverse relaxation time map into a second pore radius distribution curve based on 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 based on the first pore radius curve and the second pore radius distribution curve. The specific implementation can be seen from the description of the present embodiment. Wherein, the same or similar technical features as the foregoing embodiments are not described herein. As shown in Figure 2 The method can specifically include:

[0063] S210, 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.

[0064] S220, determining the gas-measured pore information of the standard core, and obtaining the core scanning data of the standard core by a computer tomography device.

[0065] S230, treating the standard core into a saturated water core, and obtaining the transverse relaxation time map of the saturated water core by a magnetic resonance device.

[0066] S240, determining the tortuosity, pore fractal dimension, and scanning pore radius distribution curve of the core sample based on the core scanning data, and determining a first pore radius distribution curve based on the scanning pore distribution curve and the gas-measured pore information.

[0067] The tortuosity can be understood as a dimensionless parameter representing the degree of tortuosity of the fluid flow path in the rock sample pores, and is used to quantify the complexity of the pore structure. For example, the greater the tortuosity, the more complex the rock sample pores, and the greater the resistance of the fluid passing through the rock sample pores. The pore fractal dimension can be understood as a parameter for describing the complexity of the rock pore structure, such as roughness. The pore fractal dimension reflects the irregularity of the pore shape and distribution, and is usually between 2 and 3. It can be understood that a higher pore fractal dimension indicates a more complex and irregular pore structure. The scanning pore radius distribution curve can be understood as a curve of the pore radius distribution of the rock sample obtained based on the CT scanning data, which is used to reflect the proportion of different size pores inside the rock sample, and can be extracted from the CT scanning image by image processing technology. The first pore radius distribution curve can be understood as a distribution curve of the pore radius determined based on the core scanning data and the gas logging pore information.

[0068] Based on the above scheme, optionally, the first pore radius distribution curve is determined based on the scanning pore distribution curve and the gas logging pore information, comprising: constructing a core three-dimensional model based on the core scanning data, determining the tortuosity, pore fractal dimension and scanning pore radius distribution curve of the core sample based on the core three-dimensional model; correcting the scanning pore radius distribution curve based on the gas logging pore information to obtain the first pore radius distribution curve.

[0069] The core three-dimensional model can be understood as a three-dimensional geometric model of the rock sample constructed based on the core scanning data by computer modeling technology, which is 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 scanning pore radius distribution curve is used to represent the mapping relationship between the pore radius and the cumulative frequency, and is used to reflect the proportion of different size pores inside the rock sample.

[0070] In order to improve the accuracy of the core scanning data in subsequent data analysis, an optional embodiment, before constructing the core three-dimensional model based on the core scanning data, further comprising: preprocessing the core scanning data, wherein the preprocessing can include parameter setting, noise reduction and threshold segmentation, etc.

[0071] The above scheme combines the core scanning data to construct the core three-dimensional model, and determines the tortuosity, pore fractal dimension and scanning pore radius distribution curve, and then corrects the scanning pore radius distribution curve to finally obtain the first pore radius distribution curve, which not only improves the accuracy of the pore structure analysis, but also more truly reflects the internal characteristics of the rock sample, and provides a reliable scientific basis for oil and gas resource development.

[0072] Due to the limited resolution of CT, the pore volume smaller than the minimum resolution cannot be calculated. If not corrected, the cumulative frequency in the cumulative distribution curve of pore radius starts from 0, without considering the influence of the pores smaller than the minimum resolution. NMR is a representation of the whole pore radius of the core. If the two curves are fitted, the error is large. Therefore, the cumulative volume of the pores smaller than the minimum resolution of the digital core is calculated by the gas logging core porosity, and the pore radius distribution curve of the CT is corrected.

[0073] On the basis of the above scheme, optionally, the correcting the scanning pore radius distribution curve according to the gas logging pore information to obtain a first pore radius distribution curve, comprising: determining a 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 scanning pore radius distribution curve according to the gas logging cumulative frequency to obtain the first pore radius distribution curve.

[0074] The gas logging cumulative frequency can be understood as the cumulative frequency distribution of pores of different pore sizes in the total pores according to the gas logging pore information. For example, the gas logging cumulative frequency can be the proportion of pores with a certain pore size in the total pores, which can be expressed in percentage. The gas logging cumulative frequency can provide statistical data about the pore size distribution.

[0075] 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 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 scanning pore radius distribution curve is obtained, that is, the cumulative frequency corresponding to the cumulative volume of the pores smaller than the minimum resolution of the CT. On this basis, the cumulative frequency in the scanning pore radius distribution curve is recalculated, for example, the cumulative frequency in the scanning pore radius distribution curve can be increased by the frequency difference value 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 scattered points, not continuous functions. Therefore, the cumulative frequency of the first pore radius distribution curve at any pore radius r p The transverse relaxation time T2(i) corresponding to S(i) in the second pore radius distribution curve does not necessarily have a value at this time, which can be calculated by interpolation. Specifically, in the region larger than the minimum resolution of the CT, the cumulative frequency S(i) at the pore radius r p (i) is selected, and S(i) is used to interpolate the T2 spectrum to obtain T2(i) corresponding to the pore radius r p (i).

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

[0078] In the embodiments of the application, specifically, a target conversion coefficient for converting the transverse relaxation time into the pore radius and a power index corresponding to the target conversion coefficient can be determined according to the tortuosity, the pore fractal dimension, the pore radius in the first pore radius distribution curve and a preset conversion function of the core three-dimensional model; the target conversion coefficient and the power index are brought 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.

[0079] The preset conversion function can be understood as a function for converting the transverse relaxation time into the pore radius. The target conversion coefficient can be understood as a coefficient used in the preset conversion function for adjusting the proportional relationship in the conversion process. The power index can be understood as an index used in the preset conversion function for reflecting the nonlinear relationship between the pore radius and the transverse relaxation time. The target conversion function can be understood as a conversion function obtained by bringing the target conversion coefficient and the power index into the preset conversion function. The target conversion function can be used to convert the transverse relaxation time into the pore radius. The second pore radius distribution curve can be understood as a pore radius distribution curve obtained by converting the transverse relaxation time map, 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 scheme, specifically, the preset conversion function is:

[0081]

[0082] wherein r p is the pore radius; C is the conversion coefficient; τ is the tortuosity; D is the pore fractal dimension; n is the power index; and T2 is the transverse relaxation time.

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

[0084]

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

[0086] The transverse surface relaxation rate can be understood as a transverse relaxation rate per unit pore surface area, and is used to reflect the interaction strength between the pore surface and fluid molecules, and has an influence on the T2 relaxation time. In order to simplify the use of the preset conversion function, the preset conversion function can be taken as a logarithm on both sides, so as to calculate the conversion coefficient C and the power index n by the least square method. Specifically, the preset conversion function after taking the logarithm is:

[0087]

[0088] The above scheme combines the tortuosity, pore fractal dimension and first pore radius distribution curve of the three-dimensional core model with the transverse relaxation time map, determines the target conversion coefficient and power index corresponding to the core sample by using the preset conversion function, accurately converts the relaxation time into the pore radius, generates the second pore radius distribution curve, improves the accuracy and reliability of the pore structure analysis, and provides a scientific basis for optimizing the water drive scheme and improving the producing degree of crude oil in the micro-pore.

[0089] S260, 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.

[0090] On the basis of the above scheme, optionally, the 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, comprising: determining a critical pore radius according to the scanning resolution corresponding to the core scanning data, splicing the curve segment greater than or equal to the critical pore radius in the first pore radius curve with the curve segment less than the critical pore radius in the second pore radius curve, to obtain the rock pore radius distribution curve corresponding to the rock sample.

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

[0092] The technical scheme of the embodiment of the present application firstly constructs a core three-dimensional model through core scanning data, determines tortuosity, pore fractal dimension and scanning pore radius distribution curve, and corrects the scanning pore radius distribution curve by using gas measurement pore information to obtain a first pore radius distribution curve; then, the preset conversion function and nuclear magnetic resonance technology are used to convert the transverse relaxation time map into a second pore radius distribution curve; finally, the critical pore radius is determined according to the resolution of the core scanning data, the first pore radius distribution curve and the second pore radius distribution curve are spliced, and the pore radius distribution curve of the rock sample is generated, which not only improves the accuracy and reliability of the pore structure analysis, but also provides a scientific basis for optimizing the water drive scheme and improving the producing degree of crude oil in the micro-pore.

[0093] Embodiment three

[0094] Figure 3 The flowchart of the method for determining the full-size pore radius distribution curve of the rock provided in the third embodiment of the present application is an optional embodiment of the above-mentioned embodiment. The method comprises:

[0095] Firstly, 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 into a CT scanning device for scanning to obtain core scanning data, which needs to reach the minimum resolution that can be achieved by the CT device under the size of the core to ensure that the smallest pore of the core under the constraint of the CT device can be scanned as much as possible.

[0097] Then, the core scanning data is imported into related software to construct a digital core, and the digital core needs to go through parameter setting, noise reduction, threshold segmentation (pore and skeleton division) and three-dimensional reconstruction steps, and finally forms a core three-dimensional model, as shown in Figure 5 , and determines the tortuosity, fractal dimension and core pore radius cumulative distribution curve of the core three-dimensional model, and the CT pore radius distribution curve is as shown in Figure 6 .

[0098] Then, according to the target oilfield formation water ion composition table, the simulated formation water is configured, and then the dry core is subjected to vacuum saturation water operation, and the core saturated with formation water is placed into a nuclear magnetic resonance device to measure the nuclear magnetic resonance T2 map of the core in the saturated water state.

[0099] Then, the conversion coefficient is solved by using the CT pore radius distribution curve and the nuclear magnetic resonance T2 map, and the specific operation is as follows:

[0100] 1) Draw the cumulative distribution curves of the CT pore radius and the nuclear magnetic resonance T2 relaxation time, and the results are as shown in Figure 7As shown, due to the limitation of CT device resolution, the pores smaller than the minimum resolution part cannot be identified, and it is necessary to combine the core gas measured porosity data to correct the CT pore radius cumulative distribution curve (the first pore radius distribution curve);

[0101] 2) In the area greater than the CT minimum resolution, the cumulative frequency of the pore radius r p (i) is S(i), and the T2 spectrum is interpolated using S(i) to obtain the corresponding T2(i) corresponding to any pore radius r p (i);

[0102] 3) On the basis of the original mercury injection and nuclear magnetic resonance nonlinear power function conversion method, the quantitative indicators for quantitatively describing the irregularity (tortuosity τ) and roughness (fractal dimension D) of the core (porous medium) in the oil and gas industry are introduced to correct the conversion relationship:

[0103]

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

[0105]

[0106] Take the natural logarithm of both sides of formula (4), and the following formula (5) can be obtained:

[0107]

[0108] Where C=(ρ2F 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 index; T2 is the transverse relaxation time; F s is the pore shape factor; and ρ2 is the transverse surface relaxation rate.

[0109] According to the least square method, the parameters of formula (5) are fitted and calculated to obtain the minimum error C and n values, wherein the fractal dimension and tortuosity have been obtained in the above steps. As shown in FIG. 4, the relaxation time in the T2 spectrum can be converted into the pore radius distribution, and from the fitting results, the test data and the fitting line have good consistency in shape, the fitting accuracy is high, and the rationality of the fitting method is proved. The C and n values obtained are brought into formula (5), and the relaxation time in the T2 spectrum can be converted into the pore radius distribution, that is, the converted nuclear magnetic T2 spectrum, as shown in FIG. 5. Figure 8 Figure 9

[0110] ​​Finally, taking the CT device 10 μm resolution as a demarcation point, then the CT pore radius distribution and the pore radius distribution curve after the nuclear magnetic resonance conversion are spliced together, and the repeated pore radius distribution near the junction point is removed, and finally the full-size distribution curve of the rock pore structure is obtained, realizing the full-size distribution quantitative characterization of the reservoir rock pore radius distribution, as shown in Figure 10 .

[0111] In order to verify the accuracy of the nuclear magnetic resonance T2 spectrum conversion and the full-size characterization method of the rock pore structure proposed in the experiment, the present application uses the same core to carry out the mercury injection experiment, and compares the mercury injection core pore radius distribution curve with the full-size distribution curve of the rock pore radius in the embodiment in the same coordinate system, as shown in Figure 11 It can be obviously obtained that the core pore radius distribution curves measured by the two different methods have the same trend, which proves the accuracy of the full-size characterization method of the rock pore structure proposed in the present application. The micron pores and even nanometer pores that cannot be measured by the mercury injection method can be characterized, which has a strong guiding significance for the development of complex core carbonate reservoirs.

[0112] The technical scheme of the embodiment of the present application first carries out oil washing and drying treatment on the rock sample according to the national standard to obtain a standard core, then uses a high-resolution CT scanning device to obtain the internal structure data of the core and constructs a three-dimensional model of the core to determine the tortuosity, fractal dimension and pore radius cumulative distribution curve, then measures the T2 relaxation time spectrum of the core in the saturated water state by the nuclear magnetic resonance technology, and combines the CT pore radius distribution curve and the gas measured porosity data to solve the conversion coefficient by the least square method, converts the T2 relaxation time into the pore radius distribution, finally, taking the minimum resolution (10 μm) of the CT device as the demarcation point, splicing the CT pore radius distribution and the pore radius distribution curve after the nuclear magnetic resonance conversion, removing the repeated part, and generating 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 drive scheme and improving the producing degree of the crude oil in the micro-pore, and significantly improves the development efficiency and recovery of the carbonate reservoir.

[0113] Embodiment four

[0114] Figure 4 It is a structure schematic view of a rock pore radius full-size distribution curve determination device provided by the embodiment three of the present application. As shown in Figure 4 The device comprises 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] The standard core obtaining module 410 is configured to obtain a rock sample, and pre-process the rock sample to obtain a standard core, wherein the pre-processing at least includes oil washing processing and / or drying processing.

[0116] The technical scheme of the embodiment of the present application can effectively remove impurities and moisture in the rock sample through the standard core obtaining module 410, which pre-processes the rock sample to obtain a standard core, and can improve the reliability of experimental data and provide reliable data support for subsequent processing of the rock sample, because the pre-processing at least includes oil washing processing and / or drying processing.

[0117] In the above scheme, the pore radius distribution determining module can comprise a first pore radius distribution curve determining sub-module, a second pore radius distribution curve determining sub-module, and a rock pore radius distribution curve determining sub-module.

[0118] The first pore radius distribution curve determination submodule is configured to determine tortuosity, pore fractal dimension, and a 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 information of the empty core.

[0119] On the basis of the above scheme, the first pore radius distribution curve determination submodule can optionally include a scanning pore radius distribution curve determination unit and a first pore radius distribution curve acquisition unit.

[0120] The scanning pore radius distribution curve determination unit is configured to construct a core three-dimensional model according to the core scanning data, and determine tortuosity, pore fractal dimension, and a scanning pore radius distribution curve of the core sample according to the core three-dimensional model, wherein the scanning pore radius distribution curve is used to represent a mapping relationship between a pore radius and an accumulated frequency; and the first pore radius distribution curve acquisition unit is configured to correct the scanning pore radius distribution curve according to the gas logging information of the empty core to obtain a first pore radius distribution curve.

[0121] On the basis of the above scheme, the first pore radius distribution curve acquisition unit can be specifically configured to determine a gas logging accumulated frequency corresponding to the core sample according to the gas logging information of the empty core, and correct the accumulated frequency corresponding to the pore radius in the scanning pore radius distribution curve according to the gas logging accumulated frequency to obtain the first pore radius distribution curve.

[0122] On the basis of the above scheme, the second pore radius distribution curve determination submodule can optionally include a power exponent determination unit and a second pore radius distribution curve acquisition unit.

[0123] The power index determination unit is configured to determine a target conversion coefficient for converting the transverse relaxation time into the pore radius and a power index corresponding to the target conversion coefficient according to the tortuosity of the core three-dimensional model, the pore fractal dimension, the pore radius in the first pore radius distribution curve, and a preset conversion function.

[0124] Optionally, the preset conversion function is as follows:

[0125]

[0126] wherein r p is the pore radius; C is the conversion coefficient; τ is the tortuosity; D is the pore fractal dimension; n is the power index; and T2 is the transverse relaxation time.

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

[0128] The rock pore radius full-size distribution curve determination device provided in the embodiments of the present application can execute the rock pore radius full-size distribution curve determination method provided in any of the embodiments of the present application, and has the corresponding function modules and beneficial effects of the execution method.

[0129] Embodiment five

[0130] Figure 12 A structural schematic diagram of an electronic device 10 that can be used to implement embodiments of the present application is shown. The electronic device is intended to represent various forms of digital computers, such as laptops, desktops, tablets, personal digital assistants, servers, blade servers, mainframes, and other appropriate computers. The electronic device can also represent various forms of mobile devices, such as personal digital assistants, cellular telephones, smartphones, wearable devices (e.g., headsets, glasses, watches, etc.), and other similar computing devices. The components shown here, their connections and relationships, and their functions, are meant to be examples only, and are not intended to limit the implementations of the present application described and / or claimed in this document.

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

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

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

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

[0135] The various embodiments of the systems and techniques described above can be implemented in digital electronic circuitry, integrated circuitry, a field programmable gate array (FPGA), an application specific integrated circuit (ASIC), a system on a chip (SOC), a load programmable logic device (CPLD), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include implementation 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 special or general purpose, coupled to receive data and instructions from, and to transmit data and instructions to, a storage system, at least one input device, and at least one output device.

[0136] Computer programs used to implement the processes of the application 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 program

[0137] In the context of the present application, a computer-readable storage medium can be a tangible medium that can contain or store computer programs for use by or in connection with an instruction execution system, apparatus, or device. Computer-readable storage media can include, but are not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium can be a machine-readable signal medium. More specific examples of the machine-readable storage medium will include one or more lines of electrical connections, portable computer disks, hard disk drives, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or Flash memory), optical fibers, portable compact disc read-only memories (CD-ROMs), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0138] To provide for interaction with a user, the systems and techniques described here 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) by which the user can provide input to the electronic device. Other kinds of devices can be used to provide for interaction with a user as well; for example, 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, speech, or tactile input.

[0139] The systems and techniques described here can be implemented in a computing system that includes a back end component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front end component (e.g., a user computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such back end, middleware, or front end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0140] The computing system can include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other. A server can be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system, to solve the defects of large management difficulty and weak business scalability in traditional physical host and VPS service.

[0141] It should be understood that the various forms of flow shown above can be re-ordered, added to, or deleted from without departing from the scope of the present disclosure. For example, the steps recited in the present disclosure can be executed in parallel, executed in sequence, or executed in a different order, as long as the desired results of the present disclosure are achieved, and the present disclosure is not limited herein.

[0142] The specific embodiments described above are not intended to be limiting, and persons skilled in the art will appreciate that various modifications, combinations, sub-combinations and alternatives can be made to the specific embodiments without departing from the spirit and principles of the disclosure. Accordingly, the disclosure is not limited to the specific embodiments described above, but only by the scope of the appended claims.

Claims

1. A method for determining the full-size distribution curve of rock pore radius, characterized in that, include: Obtain rock samples and pre-treat the rock samples to obtain standard rock cores, wherein the pre-treatment includes at least oil washing and drying. The gas porosity information of the standard core was determined, and the core scanning data of the standard core was obtained by a computed tomography (CT) scanner. The standard core was processed into a water-saturated core, and the transverse relaxation time spectrum of the water-saturated core was obtained by magnetic resonance imaging. The rock pore radius distribution curve corresponding to the rock sample is determined based on the core scanning data, the gas porosity information, and the transverse relaxation time spectrum. The step of determining the rock pore radius distribution curve corresponding to the rock sample based on the core scanning data, the gas porosity information, and the transverse relaxation time spectrum includes: The tortuosity, pore fractal dimension, and scanning pore radius distribution curve of the core sample are determined based on the core scanning data. The first pore radius distribution curve is determined based on the scanning pore distribution curve and the gas pore information. The transverse relaxation time profile is converted into a second pore radius distribution curve based on the tortuosity, the pore fractal dimension, and the first pore radius distribution curve. The rock pore radius distribution curve corresponding to the rock sample is determined based on the first pore radius curve and the second pore radius distribution curve; The step of converting the transverse relaxation time profile into a second pore radius distribution curve based on the tortuosity, the pore fractal dimension, and the first pore radius distribution curve includes: Based on the tortuosity of the three-dimensional core model, the fractal dimension of the pores, the pore radius in the first pore radius distribution curve, and the preset conversion function, the target conversion coefficient for converting the transverse relaxation time into the pore radius and the power exponent corresponding to the target conversion coefficient are determined. The target transformation coefficient and the power exponent are substituted into the preset transformation function to obtain the target transformation function. Based on the target transformation function, the tortuosity, the fractal dimension and the lateral relaxation time, the pore radius corresponding to the lateral relaxation time is determined to obtain the second pore radius distribution curve. The preset conversion function is: ； in, Where is the pore radius; These are the conversion factors; For tortuosity; The fractal dimension of the pores; It is a power exponent; This refers to the lateral relaxation time; The step of determining the rock pore radius distribution curve corresponding to the rock sample based on the first pore radius curve and the second pore radius distribution curve includes: The critical pore radius is determined based on the scanning resolution corresponding to the core scanning data. The curve segments of the first pore radius curve that are greater than or equal to the critical pore radius are spliced ​​with the curve segments of the second pore radius curve that are less than the critical pore radius to obtain the rock pore radius distribution curve corresponding to the rock sample.

2. The method according to claim 1, characterized in that, The step of determining the first pore radius distribution curve based on the scanned pore distribution curve and the gas-measuring pore information includes: A three-dimensional model of the core is constructed based on the core scanning data. The tortuosity, pore fractal dimension, and scanning pore radius distribution curve of the core sample are determined based on the three-dimensional model of the core. 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 based on the gas pore information to obtain a first pore radius distribution curve.

3. The method according to claim 2, characterized in that, The step of correcting the scanned pore radius distribution curve based on the gas pore information to obtain a first pore radius distribution curve includes: The gas measurement cumulative frequency corresponding to the core sample is determined based on the gas measurement porosity information, and the cumulative frequency corresponding to the pore radius in the scanned pore radius distribution curve is corrected based on the gas measurement cumulative frequency to obtain the first pore radius distribution curve.

4. A device for determining the 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. The pre-treatment includes at least oil washing and drying. The standard core analysis module is used to determine the gas porosity information of the standard core and to acquire the core scanning data of the standard core using a computed tomography (CT) scanner. The saturated water core analysis module 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. A pore radius distribution determination module is used to determine the rock pore radius distribution curve corresponding to the rock sample based on the core scanning data, the gas pore information and the transverse relaxation time spectrum. The pore radius distribution determination module includes: The first pore radius distribution curve determination submodule is used to determine the tortuosity, pore fractal dimension and scanning pore radius distribution curve of the core sample based on the core scanning data, and to determine the first pore radius distribution curve based on the scanning pore distribution curve and the gas pore information. The second pore radius distribution curve determination submodule is used to convert the transverse relaxation time spectrum into a second pore radius distribution curve based on the tortuosity, the pore fractal dimension, and the first pore radius distribution curve. The rock pore radius distribution curve determination submodule is used to determine the rock pore radius distribution curve corresponding to the rock sample based on the first pore radius curve and the second pore radius distribution curve; The second pore radius distribution curve determination submodule includes: The power exponent determination unit is used to determine the target conversion coefficient for converting the transverse relaxation time into the pore radius and the power exponent corresponding to the target conversion coefficient based on the tortuosity of the core three-dimensional model, the pore fractal dimension, the pore radius in the first pore radius distribution curve, and the preset conversion function. The second pore radius distribution curve acquisition unit is used to input the target conversion coefficient and the power exponent into the preset conversion function to obtain the target conversion function, and to determine the pore radius corresponding to the transverse relaxation time based on the target conversion function, the tortuosity, the fractal dimension and the transverse relaxation time to obtain the second pore radius distribution curve. The preset conversion function is: ， in, Where is the pore radius; These are the conversion factors; For tortuosity; The fractal dimension of the pores; It is a power exponent; This refers to the lateral relaxation time; The rock pore radius distribution curve determination submodule is specifically used for: The critical pore radius is determined based on the scanning resolution corresponding to the core scanning data. The curve segments of the first pore radius curve that are greater than or equal to the critical pore radius are spliced ​​with the curve segments of the second pore radius curve that are less than the critical pore radius to obtain the rock pore radius distribution curve corresponding to the rock sample.

5. An electronic device, characterized in that, The electronic device includes: At least one processor; and A memory communicatively connected to the at least one processor; wherein, The memory stores a computer program that can be executed by the at least one processor, the computer program being executed by the at least one processor to enable the at least one processor to perform the method for determining the full-size distribution curve of rock pore radius as described in any one of claims 1-3.

6. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions that, when executed by a processor, implement the method for determining the full-size distribution curve of rock pore radius as described in any one of claims 1-3.

Images