Method, device, equipment and storage medium for characterizing pores in microstructure of volcanic rocks

Abstract

The embodiment of the present application discloses a method, device, equipment and storage medium for characterizing pores in the microstructure of volcanic rocks. The method includes: determining the pore profile of pores in at least one pore image; determining the pore shape of at least one pore based on the pore profile of at least one pore; determining the initial pore radius of at least one pore based on the pore shape of at least one pore; determining the pore radius range based on the initial pore radius of at least one pore; determining the pressure range corresponding to the pore radius range; determining the target relationship curve corresponding to the pressure range from the relationship curve; obtaining the standard relationship curve based on the target relationship curve; determining the target volume value of at least one pore from the standard relationship curve; determining the equivalent pore radius of at least one pore; determining the final pore radius of at least one pore based on the pore shape and equivalent pore radius of at least one pore. The present application can obtain a more accurate pore radius.

Description

Technical Field

[0001] The embodiments of the present application relate to the technical field of rock structure characterization, and in particular to a method, device, equipment and storage medium for characterizing pores in the microstructure of volcanic rocks. Background Art

[0002] In recent years, with the further development of oil and gas exploration, the exploration of conventional oil and gas reservoirs has become increasingly difficult, and the exploration field of oil and gas resources has gradually shifted to unconventional oil and gas.

[0003] So far, volcanic reservoirs, as one of the unconventional oil and gas resources, have made significant progress in exploration and development in my country. Compared with conventional reservoirs, volcanic reservoirs have special geological characteristics, which makes the microscopic pore throat matching relationship of reservoir rocks more complex and difficult to accurately characterize, which brings many difficulties to exploitation and development. In order to accurately characterize the geological characteristics of oil and gas reservoirs and exploit and transform oil and gas reservoirs more efficiently and economically, it is necessary to analyze and study the microscopic heterogeneity of volcanic reservoirs and characterize the rock microstructure of volcanic reservoirs in a detailed and accurate manner. However, how to accurately characterize the microstructure of volcanic reservoirs remains a problem to be solved. Summary of the invention

[0004] The embodiment of the present application provides a method, device, equipment and storage medium for characterizing the pores of the microstructure of volcanic rocks, which can obtain a more accurate pore radius, thereby achieving accurate microstructural characterization of dense volcanic rock samples. The technical solution is as follows:

[0005] On the one hand, an embodiment of the present application provides a method for characterizing pores in a volcanic rock microstructure, the method comprising:

[0006] Acquiring at least one pore image of a first rock sample of a volcanic rock core;

[0007] determining a pore profile of a pore in at least one of the pore images;

[0008] determining a pore shape of at least one of the pores based on a pore profile of at least one of the pores;

[0009] determining an initial pore radius of at least one of the pores based on a pore shape of at least one of the pores;

[0010] determining a pore radius range of the first rock sample based on an initial pore radius of at least one of the pores;

[0011] Based on the pore radius range of the first rock sample, determining a pressure range corresponding to the pore radius range;

[0012] Determining a target relationship curve corresponding to the pressure range from the relationship curve, wherein the relationship curve refers to a relationship curve between pressure and volume obtained by performing a constant-rate mercury injection experiment on a second rock sample of the volcanic rock core;

[0013] Based on the target relationship curve, a standard relationship curve between pressure and volume is obtained;

[0014] determining a target volume value for at least one of the pores from the standard relationship curve;

[0015] determining an equivalent pore radius of at least one of the pores based on a target volume value of at least one of the pores;

[0016] Based on the pore shape of at least one of the pores and the equivalent pore radius, a final pore radius of at least one of the pores is determined.

[0017] On the other hand, an embodiment of the present application provides a volcanic rock microstructure pore characterization device, the device comprising:

[0018] An image acquisition module, configured to acquire at least one pore image of a first rock sample of a volcanic rock core;

[0019] a contour determination module, configured to determine a pore contour of a pore in at least one of the pore images;

[0020] a shape determination module for determining a pore shape of at least one of the pores based on a pore profile of at least one of the pores;

[0021] a radius determination module, configured to determine an initial pore radius of at least one of the pores based on a pore shape of at least one of the pores;

[0022] a range determination module, configured to determine a pore radius range of the first rock sample based on an initial pore radius of at least one of the pores;

[0023] The range determination module is further used to determine a pressure range corresponding to the pore radius range based on the pore radius range of the first rock sample;

[0024] a curve determination module, used to determine a target relationship curve corresponding to the pressure range from a relationship curve, wherein the relationship curve refers to a relationship curve between pressure and volume obtained by performing a constant-rate mercury injection experiment on a second rock sample of the volcanic rock core;

[0025] The curve determination module is further used to obtain a standard relationship curve between pressure and volume based on the target relationship curve;

[0026] a volume determination module, configured to determine a target volume value of at least one of the pores from the standard relationship curve;

[0027] The radius determination module is further used to determine an equivalent pore radius of at least one of the pores based on a target volume value of at least one of the pores;

[0028] The radius determination module is further used to determine a final pore radius of at least one of the pores based on the pore shape of at least one of the pores and the equivalent pore radius.

[0029] On the other hand, an embodiment of the present application provides a computer device, comprising a processor and a memory, wherein the memory stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, the at least one program, the code set or the instruction set is loaded and executed by the processor to implement the volcanic rock microstructure pore characterization method as described in the above aspects.

[0030] On the other hand, an embodiment of the present application provides a computer-readable storage medium, wherein the storage medium stores at least one instruction, at least one program, a code set, or an instruction set, and the at least one instruction, the at least one program, the code set, or the instruction set is loaded and executed by a processor to implement the method for characterizing the pores of the volcanic rock microstructure as described in the above aspects.

[0031] The beneficial effects brought by the technical solution provided by the embodiments of the present application may include:

[0032] By determining the pore shape of each pore, an initial pore radius is determined based on the pore shape; a final standard constant-rate mercury injection curve is determined based on the initial pore radius; and finally, a final pore radius of the pore is determined based on the constant-rate mercury injection curve and the pore shape. In the embodiment of the present application, a scanning electron microscope is used to correct the pore size in the constant-rate mercury injection result to obtain a more accurate pore radius, thereby achieving accurate microstructural characterization of dense volcanic rock samples. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0034] Figure 1 is a flow chart of a method for characterizing pores in a volcanic rock microstructure provided by an embodiment of the present application;

[0035] Figure 2 is a schematic diagram of determining a standard relationship curve provided by an embodiment of the present application;

[0036] Figure 3 is a schematic diagram of determining shape distribution frequency provided by an embodiment of the present application;

[0037] Figure 4 is a schematic diagram of pore size distribution frequency provided by the related art;

[0038] Figure 5 is a schematic diagram of pore size distribution frequency provided in an embodiment of the present application;

[0039] Figure 6 is a schematic diagram of a volcanic rock microstructure pore characterization device provided by an embodiment of the present application;

[0040] Figure 7 It is a schematic diagram of the structure of a computer device provided by an embodiment of the present application. DETAILED DESCRIPTION

[0041] In order to make the objectives, technical solutions and advantages of the present application clearer, the implementation methods of the present application will be further described in detail below with reference to the accompanying drawings.

[0042] In the relevant technologies, there are many technical means for characterizing rock microstructure, which can be divided into qualitative characterization means and quantitative characterization means according to their properties. Qualitative characterization mainly studies the types, morphological distribution and connectivity of pores and throats. The main technical means used are cast thin sections, field emission scanning electron microscope (FESEM) and environment scanning electron microscope (ESEM). Scanning electron microscope and cast thin sections can directly observe the pore throat structure of tight reservoirs, which are mainly used for qualitative evaluation and can also provide quantitative characterization results; however, the resolution of cast thin sections for microscopic pores is low, and the scanning electron microscope will lead to a decrease in the sample analysis field of view as the resolution increases. The extraction of two-dimensional pore throat structure information after subsequent image processing is limited, and a large amount of experimental data is required to ensure that the results are representative. Quantitative characterization is a quantitative analysis and parameter characterization of pore throat size and distribution. The main technical means used at this stage are constant speed pressure pump, high pressure pressure pump and nuclear magnetic resonance. Based on the understanding that mercury is a non-wetting phase for most rock-forming mineral components, the mercury injection method injects mercury into the rock sample at a constant pressure or constant speed. According to the experimental data between the amount of mercury injected and the mercury injection pressure, the pore throat size and distribution curve of the rock sample are calculated. Among them, the constant rate mercury injection method can distinguish the pores and throats of the rock sample, and realize the size and distribution characterization of the pores and throats of the rock sample. However, in data analysis and processing, the pores of the rock sample are generally equivalent to spheres to determine the pore throat radius, which will cause a large difference between the calculated results and the actual situation. Nuclear magnetic resonance technology measures the T2 spectrum curve of the saturated rock sample and then calculates the porosity, pore size and distribution of the rock sample. Compared with other methods, it is faster, more efficient and non-destructive, but this method cannot directly obtain the pore throat radius distribution information, and other quantitative characterization methods are needed.

[0043] Please refer to Figure 1 , which shows a flow chart of a method for characterizing pores in a volcanic rock microstructure provided by an embodiment of the present application. The method can be applied to a computer device with data processing and storage capabilities such as a PC (Personal Computer), for example, the computer device can be a terminal or a computer device. The above method may include the following steps.

[0044] Step 101: Acquire at least one pore image of a first rock sample of a volcanic rock core.

[0045] Volcanic rock is a porous stone formed after a volcanic eruption. In a possible implementation, the volcanic rock in the embodiment of the present application can also be referred to as a volcanic rock reservoir. A core refers to a cylindrical rock sample taken out of a hole using a core ring drill bit and other coring tools according to the needs of geological work or engineering. The core is an important physical material for studying and understanding underground geology and mineral conditions.

[0046] For example, a pore image refers to an image containing pores. Pores refer to larger storage spaces surrounded by rock particles and are basic storage spaces for fluids.

[0047] For example, the technician takes a plug rock sample of a volcanic rock reservoir, washes the salt and dries it according to GB / T 29172-2012 Core Analysis Method, and measures the basic physical properties of the core such as porosity and permeability (the basic physical properties have a reference for the selection of the core and the subsequent processing process). The technician cuts the core into two sections, which are respectively recorded as the first rock sample and the second rock sample, and then dries it again. The technician selects the first rock sample for scanning electron microscopy experiments, and observes the fresh cross-section sample (the first rock sample) using a scanning electron microscope according to the provisions of SY / T5362-2014 Rock Sample Scanning Electron Microscope Analysis Method: The technician first scans and obtains a global image of the first rock sample at a lower magnification (for example, 500 times); then, the global image is divided into different areas, and the magnification of the display microscope is increased (for example, from 500 times to 1000 times), and each pore image in the first rock sample is photographed to obtain at least one pore image of the first rock sample of the volcanic rock core. The technician may upload at least one pore image of the first rock sample of the volcanic rock core to a computer device, so that the computer device processes the at least one pore image.

[0048] Step 102 : determining a pore profile of a pore in at least one pore image.

[0049] In a possible implementation, the pore contours of the pores in the at least one pore image are determined by image processing software, for example, the pore contours of the pores in the at least one pore image are extracted by Photoshop software.

[0050] Step 103: determining a pore shape of at least one pore based on a pore profile of at least one pore.

[0051] The pore shape of a pore may refer to a cross-sectional shape of the pore.

[0052] In a possible implementation, the computer device may determine the pore shape of the pore based on the pore profile characteristics of the pore. For example, the computer device may determine the pore factor of the pore based on the pore profile, and then determine the pore shape based on the pore factor of the pore. How the computer device determines the pore shape of the pore is not described here.

[0053] Step 104 : determining an initial pore radius of at least one pore based on the pore shape of at least one pore.

[0054] The initial pore radius refers to the initial radius determined based on the pore shape.

[0055] In a possible implementation, the pore shape includes any one of the following: a first triangle (trilateral 1), a second triangle (trilateral 2), a third triangle (trilateral 3), a first quadrilateral (quadrilateral 1), a second quadrilateral (quadrilateral 2), a first pentagon (pentagon 1) and a second pentagon (pentagon 2). The first triangle, the second triangle and the third triangle are mutually different triangles, the first triangle and the second triangle are irregular triangles, the third triangle is a regular triangle, the first quadrilateral and the second quadrilateral are two different quadrilaterals, the second quadrilateral is a regular quadrilateral, the first pentagon and the second pentagon are different pentagons, and the second pentagon refers to a regular pentagon (it should be noted that the triangle in the embodiment of the present application refers to a figure with three vertices, the quadrilateral refers to a figure with four vertices, and the pentagon refers to a figure with five vertices).

[0056] When the pore shape is a first triangle, the initial pore radius can be half of the width of the first triangle; when the pore shape is a second triangle, the initial pore radius can be half of the height of the second triangle; when the pore shape is a third triangle, the initial pore radius can be determined by the following formula: initial pore radius = 2*area of ​​the third triangle / perimeter of the third triangle; when the pore shape is a first quadrilateral, the initial pore radius can be half of the height of the first quadrilateral; when the pore shape is a second quadrilateral, the initial pore radius can be determined by the following formula: initial pore radius = 2*area of ​​the second quadrilateral / perimeter of the second quadrilateral; when the pore shape is a first pentagon, the initial pore radius can be half of the width of the first pentagon; when the pore shape is a second pentagon, the initial pore radius can be the area of ​​the second pentagon / perimeter of the second pentagon.

[0057] Step 105: Determine a pore radius range of the first rock sample based on an initial pore radius of at least one pore.

[0058] Exemplarily, the smallest initial pore radius among the initial pore radius of at least one pore is determined as the lower limit of the pore radius range of the first rock sample, and the largest initial pore radius among the initial pore radius of at least one pore is determined as the upper limit of the pore radius range of the first rock sample.

[0059] Step 106: Based on the pore radius range of the first rock sample, determine a pressure range corresponding to the pore radius range.

[0060] In a possible implementation, the pressure value P corresponding to the pore radius r is determined by the following formula: c :

[0061]

[0062] Exemplarily, the computer device can calculate the pressure value corresponding to the minimum initial pore radius and the pressure value corresponding to the maximum initial pore radius, determine the pressure value corresponding to the minimum initial pore radius as the lower limit of the pressure range, and determine the pressure value corresponding to the maximum initial pore radius as the upper limit of the pressure range.

[0063] Step 107: determining a target relationship curve corresponding to the pressure range from the relationship curve.

[0064] In the embodiment of the present application, the relationship curve refers to the relationship curve between pressure and volume obtained by performing a constant rate mercury injection experiment on the second rock sample of the volcanic rock core. Exemplarily, the technician performs a constant rate mercury injection experiment on the second rock sample to obtain a relationship curve between the mercury injection pressure (referred to as pressure) and the mercury injection volume (referred to as volume) of the second rock sample. Exemplarily, the length of the first rock sample is 2 cm, and the length of the second rock sample is 3 cm.

[0065] Mercury is non-wetting to most rock-forming minerals. If pressure is applied to mercury, when the pressure of the injected mercury reaches the capillary pressure of the pore throat, the mercury will overcome the capillary resistance and enter the pore. According to the pore volume percentage and corresponding pressure of the continuously injected mercury, a relationship curve can be drawn. The constant rate mercury injection method injects mercury into the pores of the rock at a relatively low constant flow rate. The mercury injection rate is low and can be considered as a quasi-static process. During the injection of mercury, the pressure changes. By monitoring the pressure rise and fall, the pores and throats inside the rock can be separated. Mercury first enters throat 1. After the pressure rises to a certain value, mercury breaks through throat 1 and enters pore 1. At this time, the pressure drops rapidly. When pore 1 is filled with mercury, the pressure begins to rise. At this time, mercury is entering throat 2. After the pressure rises to a certain value, mercury breaks through throat 2 and enters pore 2. The pressure drops rapidly again. Similarly, mercury fills other throats and other pores in turn until all throats connected to throat 1 and the controlled pores are filled. Only then will mercury flow into the pore units controlled by throats with a smaller radius than throat 1. Repeat the above process. When the pressure reaches the upper limit of the test instrument, the experiment ends. According to the pressure fluctuations during the mercury injection process, information on the pore structure can be obtained, the number of throats and pores can be quantitatively obtained, and the characteristic parameters of the pore structure such as pore radius distribution and throat radius distribution can be obtained.

[0066] Exemplarily, the abscissa of the relationship curve is volume, and the ordinate is pressure. A curve in the relationship curve corresponding to the pressure range is determined as the target relationship curve.

[0067] Step 108: obtaining a standard relationship curve between pressure and volume based on the target relationship curve.

[0068] The target relationship curve is enlarged, combined with the enlarged target relationship curve, and the abnormal points in the target relationship curve are removed to obtain the optimized points. Since the target relationship curve is actually a curve composed of multiple points, the fluctuation of multiple points can be obtained after the target relationship curve is enlarged. The points with abnormal fluctuations in the target relationship curve are removed to obtain the optimized points. The standard relationship curve of pressure and volume is drawn based on the optimized points. The standard relationship curve can also be called the standard "M" curve. Figure 2 As shown, it shows a schematic diagram of determining a standard relationship curve provided by an embodiment of the present application. The computer device determines the target relationship curve based on the relationship curve 210 and the pressure range, assuming that the pore radius range is 0.226 microns-0.533 microns, and determines the pressure range based on the pore radius range, then intercepts the target relationship curve 220 in the relationship curve 210, amplifies the target relationship curve 220, and removes abnormal points in the target relationship curve 220 to obtain optimized points, and draws a standard relationship curve 230 between pressure and volume based on the optimized points. (It should be noted that Figure 2This is merely an example. In possible implementations, the relationship curve may also be other types of curves, which are not limited in the embodiments of the present application).

[0069] Step 109, determining a target volume value of at least one pore from the standard relationship curve.

[0070] For example, Figure 2 As shown, V1 represents a target volume value of a pore, V2 represents a target volume value of a pore, V3 represents a target volume value of a pore, and so on, to obtain a target volume value of at least one pore.

[0071] In a possible implementation, an insertion sort algorithm is used to determine a target volume value of at least one pore from a standard relationship curve.

[0072] Step 110: determining an equivalent pore radius of at least one pore based on a target volume value of at least one pore.

[0073] Exemplarily, the pores are equivalent to spheres, and the equivalent pore radius of at least one pore is determined based on a sphere volume formula and a target volume value.

[0074] The equivalent pore radius refers to the radius determined based on the equivalent pores.

[0075] The formula for the volume of a sphere is Wherein, v represents volume and r represents radius. When V in the above formula is known, r can be inferred inversely, and therefore, based on the target volume value of at least one pore, the equivalent pore radius of at least one pore can be determined.

[0076] Step 111 : determining a final pore radius of at least one pore based on the pore shape and the equivalent pore radius of at least one pore.

[0077] The final pore radius refers to the radius of the pore that is finally determined.

[0078] Exemplarily, the computer device may determine the final pore radius of the at least one pore in the following manner:

[0079] 1. Determine a radius determination method corresponding to a pore shape of at least one pore.

[0080] The radius determination method refers to the method of determining the final pore radius.

[0081] Due to different pore shapes, the methods for determining the final pore radius are different. Therefore, it is necessary to first determine the radius determination method based on the pore shape.

[0082] 2. Based on the radius determination method and the equivalent pore radius, determine the final pore radius of at least one pore.

[0083] In a possible implementation, when the pore shape of the pore is a triangle, the equivalent pore radius of the first value is determined as the final pore radius.

[0084] In the case where the pore shape of the pore is a quadrilateral, the equivalent pore radius of the second value is determined as the final pore radius.

[0085] In the case where the pore shape of the pore is a pentagon, the equivalent pore radius of the third value is determined as the final pore radius.

[0086] Exemplarily, the first value, the second value, and the third value are different values ​​from each other.

[0087] In a possible implementation, when the pore shape is a triangle, the computer device can determine the final pore radius in the following manner:

[0088]

[0089] Among them, R 三边形内切球半径 The final pore radius of a pore with a triangular pore shape, R 等效孔隙半径 represents the equivalent pore radius.

[0090] In a possible implementation, when the pore shape is a quadrilateral, the computer device may determine the final pore radius in the following manner:

[0091]

[0092] Among them, R 四边形内切球半径 The final pore radius of a quadrilateral pore, R 等效孔隙半径 represents the equivalent pore radius.

[0093] In a possible implementation, when the pore shape is a pentagon, the computer device may determine the final pore radius in the following manner:

[0094]

[0095] Among them, R 五边形内切球半径 The final pore radius of the pore with a pentagonal shape, R 等效孔隙半径 represents the equivalent pore radius.

[0096] Volcanic rock reservoirs generally undergo strong compaction and cementation (siliceous cementation, clay cementation), which results in extremely complex pore shape and size variation characteristics; if the microstructure of volcanic rocks is to be accurately characterized, it is necessary to combine the advantages of constant-rate mercury injection experiments and other methods. The embodiment of the present application provides a method for characterizing the microscopic pore structure of volcanic rocks combined with constant-rate mercury injection experiments. For the pore structure of dense volcanic rock samples, a scanning electron microscope is used to correct the pore and throat size in the constant-rate mercury injection structure, thereby accurately characterizing the microstructure of dense volcanic rock samples.

[0097] In summary, in the technical solution provided in the embodiment of the present application, the pore shape of each pore is determined, and the initial pore radius is determined based on the pore shape; the final standard constant-rate mercury injection curve is determined based on the initial pore radius; finally, the final pore radius of the pore is determined based on the constant-rate mercury injection curve and the pore shape. The embodiment of the present application uses a scanning electron microscope to correct the pore size in the constant-rate mercury injection result to obtain a more accurate pore radius, thereby realizing accurate microstructural characterization of dense rock samples of volcanic rocks.

[0098] In an exemplary embodiment, the pore shape of the pores is determined as follows:

[0099] First, based on the pore profile of the at least one pore, a cross-sectional area of ​​the at least one pore and a perimeter of the at least one pore are determined.

[0100] In a possible implementation, the pore contour is processed by ImageJ software to obtain the cross-sectional area and perimeter of at least one pore.

[0101] Second, determining a shape factor of at least one pore based on a cross-sectional area of ​​the at least one pore and a perimeter of the at least one pore.

[0102] In a possible implementation, the computer device determines the shape factor G by the following formula:

[0103]

[0104] Where A represents the cross-sectional area of ​​the pore, and P represents the perimeter of the pore. For example, the unit of the cross-sectional area in the formula may be cm 2 , the unit of circumference can be cm.

[0105] Third, based on the shape factor of the at least one pore, a pore shape of the at least one pore is determined.

[0106] like Figure 3As shown, by dividing the global image 310 of the first rock sample into different regions, each pore image 320 in the first rock sample is obtained by photographing; the pore contour 330 of the pore in the pore image 320 is determined; the pore factor of the pore is determined based on the pore contour 330; the computer device can determine the pore shape of the pore through the shape factor distribution diagram 340 (the shape factor distribution diagram 340 is used to indicate the corresponding relationship between the pore shape and the shape factor). Exemplarily, if the shape factor is located at the intersection of triangle 2 and quadrilateral 1, the shape corresponding to the lower limit closer to the shape factor is selected as the pore shape of the pore. For example, if the shape factor is 0.0218, quadrilateral 1 is determined as the pore shape of the pore, because 0.0281 is closer to 0.0217 (that is, the lower limit of the shape factor corresponding to quadrilateral 1).

[0107] In an exemplary embodiment, after the computer device determines the pore shape of at least one pore, the computer device may further perform the following steps:

[0108] First, based on the pore shape of at least one pore, the number of pores having the same pore shape is determined.

[0109] After the computer device determines the pore shape of each pore, the number of pores with the same pore shape can be determined. For example, the computer device can determine the number of pores with a pore shape of a triangle 1, the number of pores with a pore shape of a quadrilateral 1, the number of pores with a pore shape of a pentagon 1, etc.

[0110] Second, based on the number of pores with the same type of pore shape and the total number of at least one pore, determine the shape distribution frequency corresponding to the pores with different types of pore shapes. The shape distribution frequency is used to indicate the ratio of the number of pores with the same type of pore shape to the total number of at least one pore.

[0111] For example, after the computer device determines the number of pores of a certain type, the ratio of the value to the total number of pores can be determined as the shape distribution frequency of the pores of a certain type. Figure 3 As shown, it shows a schematic diagram 350 of shape distribution frequency provided by an embodiment of the present application.

[0112] In an exemplary embodiment, after the computer device determines the final pore radius of the pore, the following step may be performed: determining the final equivalent radius of the core based on the pore shape and the final pore radius of at least one pore.

[0113] Exemplarily, the computer device may determine the final equivalent radius of the core based on the final pore radius of pores with a triangular pore shape, the final pore radius of pores with a quadrilateral pore shape, and the final pore radius of pores with a pentagonal pore shape.

[0114] In a possible implementation, the computer device may determine the final equivalent radius of the core in the following manner: based on the shape distribution frequencies corresponding to pores of different types of pore shapes and the final pore radii corresponding to pores of different types of pore shapes, the final equivalent radius of the core is determined.

[0115] In a possible implementation, the final equivalent radius of the core can be determined as follows:

[0116] 1. Based on the final pore radius of at least one pore, determine a target pore radius of pores of the same type of pore shape.

[0117] Exemplarily, the final pore radii of the pores of the same type of pore shape are added and averaged to obtain the target pore radius of the pores of the same type of pore shape.

[0118] 2. The sum of the product of the shape distribution frequency corresponding to pores of the same type and the target pore radius is determined as the final equivalent radius of the core.

[0119] by Figure 3 The schematic diagram of shape distribution frequency shown in FIG. 1 is used as an example for description. Assuming that the first rock sample includes pores with pore shapes of triangle 3, quadrilateral 2, pentagon 1 and pentagon 2, the shape distribution frequency of pores with pore shape of triangle 3 is 0.3, the shape distribution frequency of pores with pore shape of quadrilateral 2 is 0.4, the shape distribution frequency of pores with pore shape of pentagon 1 is 0.2, and the shape distribution frequency of pores with pore shape of pentagon 2 is 0.1, assuming that the target pore radius of pores with pore shape of triangle 3 is radius 1, the target pore radius of pores with pore shape of quadrilateral 2 is radius 2, the target pore radius of pores with pore shape of pentagon 1 is radius 3, and the target pore radius of pores with pore shape of pentagon 2 is radius 4. Then the final equivalent radius of the core is (0.3*radius 1+0.4*radius 2+0.2*radius 3+0.1*radius 4).

[0120] For a rock sample of a Permian volcanic reservoir in a certain basin, based on the results of scanning electron microscope experiment and constant rate mercury injection experiment, the corrected volcanic rock pore radius (i.e., final pore radius) and distribution frequency are obtained according to the process proposed in the embodiment of the present application. The pore radius and distribution frequency can be referred to as the pore size distribution frequency. Comparative analysis shows that the pore size distribution frequency intuitively obtained from the scanning electron microscope experiment is consistent with the pore size distribution frequency corrected based on the embodiment of the present application, which shows that the volcanic rock microscopic pore structure characterization method combined with constant rate mercury injection experiment proposed in the embodiment of the present application is feasible and highly operable. Figure 4A schematic diagram 400 showing the uncorrected pore size distribution frequency based on a constant rate mercury injection experiment in the related art is shown. Figure 5 A schematic diagram 500 showing the pore size distribution frequency after correction based on an embodiment of the present application is shown.

[0121] The following is an embodiment of the device of the present application, which can be used to execute the embodiment of the method of the present application. For details not disclosed in the embodiment of the device of the present application, please refer to the embodiment of the method of the present application.

[0122] Please refer to Figure 6 It shows a block diagram of a volcanic rock microstructure pore characterization device provided by an embodiment of the present application. The device has the function of implementing the above method example, and the function can be implemented by hardware or by hardware executing corresponding software. The device 600 may include: an image acquisition module 610, a contour determination module 620, a shape determination module 630, a radius determination module 640, a range determination module 650, a curve determination module 660, and a volume determination module 670.

[0123] An image acquisition module 610 is used to acquire at least one pore image of a first rock sample of a volcanic rock core;

[0124] a contour determination module 620, for determining a pore contour of a pore in at least one of the pore images;

[0125] a shape determination module 630 for determining a pore shape of at least one of the pores based on a pore profile of at least one of the pores;

[0126] A radius determination module 640, configured to determine an initial pore radius of at least one of the pores based on a pore shape of at least one of the pores;

[0127] A range determination module 650, configured to determine a pore radius range of the first rock sample based on an initial pore radius of at least one of the pores;

[0128] The range determination module 650 is further configured to determine a pressure range corresponding to the pore radius range based on the pore radius range of the first rock sample;

[0129] A curve determination module 660, configured to determine a target relationship curve corresponding to the pressure range from a relationship curve, wherein the relationship curve is a relationship curve between pressure and volume obtained by performing a constant-rate mercury injection experiment on a second rock sample of the volcanic rock core;

[0130] The curve determination module 660 is further used to obtain a standard relationship curve between pressure and volume based on the target relationship curve;

[0131] A volume determination module 670, for determining a target volume value of at least one of the pores from the standard relationship curve;

[0132] The radius determination module 640 is further configured to determine an equivalent pore radius of at least one of the pores based on a target volume value of at least one of the pores;

[0133] The radius determination module 640 is further configured to determine a final pore radius of at least one of the pores based on the pore shape of at least one of the pores and the equivalent pore radius.

[0134] In summary, in the technical solution provided in the embodiment of the present application, the pore shape of each pore is determined, and the initial pore radius is determined based on the pore shape; the final standard constant-rate mercury injection curve is determined based on the initial pore radius; finally, the final pore radius of the pore is determined based on the constant-rate mercury injection curve and the pore shape. The embodiment of the present application uses a scanning electron microscope to correct the pore size in the constant-rate mercury injection result to obtain a more accurate pore radius, thereby realizing accurate microstructural characterization of dense rock samples of volcanic rocks.

[0135] In an exemplary embodiment, the shape determination module 630 is used to:

[0136] determining a cross-sectional area of ​​at least one of the pores and a perimeter of at least one of the pores based on a pore profile of at least one of the pores;

[0137] determining a shape factor of at least one of the pores based on a cross-sectional area of ​​at least one of the pores and a perimeter of at least one of the pores;

[0138] Based on a shape factor of at least one of the pores, a pore shape of at least one of the pores is determined.

[0139] In the exemplary embodiment, the radius determination module 640 includes: a mode determination unit and a radius determination unit (not shown in the figure).

[0140] a mode determination unit, configured to determine a radius determination mode corresponding to a pore shape of at least one of the pores;

[0141] The radius determination unit is used to determine a final pore radius of at least one of the pores based on the radius determination method and the equivalent pore radius.

[0142] In an exemplary embodiment, the radius determination unit is used to:

[0143] In the case where the pore shape of the pore is a triangle, determining the equivalent pore radius of the first value as the final pore radius;

[0144] In the case where the pore shape of the pore is a quadrilateral, determining the equivalent pore radius of the second value as the final pore radius;

[0145] In a case where the pore shape of the pore is a pentagon, the equivalent pore radius of the third value is determined as the final pore radius.

[0146] In an exemplary embodiment, the device further includes: a frequency determination module (not shown in the figure).

[0147] Determining the number of pores having the same pore shape based on the pore shape of at least one of the pores;

[0148] Based on the number of pores with the same type of pore shape and the total number of at least one of the pores, the shape distribution frequencies corresponding to the pores with different types of pore shapes are determined, and the shape distribution frequencies are used to indicate the ratio of the number of pores with the same type of pore shape to the total number of at least one of the pores.

[0149] In an exemplary embodiment, the radius determination module 640 is further configured to:

[0150] A final equivalent radius of the core is determined based on the pore shape of at least one of the pores and the final pore radius.

[0151] In an exemplary embodiment, the radius determination module 640 is used to:

[0152] determining a target pore radius for pores of the same type of pore shape based on the final pore radius of at least one of the pores;

[0153] The sum of the products of the shape distribution frequencies corresponding to pores of the same type of pore shape and the target pore radius is determined as the final equivalent radius of the core.

[0154] One point that needs to be explained is that the device provided in the above embodiment only uses the division of the above-mentioned functional modules as an example to implement its functions. In actual applications, the above-mentioned functions can be assigned to different functional modules according to actual needs, that is, the content structure of the device can be divided into different functional modules to complete all or part of the functions described above.

[0155] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0156] Please refer to Figure 7, which shows a schematic diagram of the structure of a computer device 700 provided in one embodiment of the present application. The computer device 700 can be used to implement the volcanic rock microstructure pore characterization method provided in the above embodiment on the computer device side. Specifically:

[0157] The computer device 700 includes a central processing unit (CPU) 701, a system memory 704 including a RAM (Random Access Memory) 702 and a ROM (Read-Only Memory) 703, and a system bus 705 connecting the system memory 704 and the central processing unit 701. The computer device 700 also includes a basic input / output system (I / O system, Input / Output system) 706 that helps to transmit information between various components in the computer, and a large-capacity storage device 707 for storing an operating system 713, application programs 714 and other program modules 715.

[0158] The basic input / output system 706 includes a display 708 for displaying information and an input device 709 such as a mouse and a keyboard for user inputting information. The display 708 and the input device 709 are connected to the central processing unit 701 through an input / output controller 710 connected to the system bus 705. The basic input / output system 706 may also include an input / output controller 710 for receiving and processing inputs from a plurality of other devices such as a keyboard, a mouse, or an electronic stylus. Similarly, the input / output controller 710 also provides output to a display screen, a printer, or other types of output devices.

[0159] The mass storage device 707 is connected to the central processing unit 701 via a mass storage controller (not shown) connected to the system bus 705. The mass storage device 707 and its associated computer readable medium provide non-volatile storage for the computer device 700. That is, the mass storage device 707 may include a computer readable medium (not shown) such as a hard disk or a CD-ROM (Compact Disc Read-Only Memory) drive.

[0160] Without loss of generality, the computer-readable medium may include computer storage media and communication media. Computer storage media include volatile and non-volatile, removable and non-removable media implemented by any method or technology for storing information such as computer-readable instructions, data structures, program modules or other data. Computer storage media include RAM, ROM, EPROM (Erasable Programmable Read-Only Memory), EEPROM (Electrically Erasable Programmable Read-Only Memory), Flash Memory or other solid-state storage devices, CD-ROM, DVD (Digital Versatile Disc) or other optical storage, cassettes, tapes, disk storage or other magnetic storage devices. Of course, those skilled in the art will know that the computer storage medium is not limited to the above. The above-mentioned system memory 704 and mass storage device 707 can be collectively referred to as memory.

[0161] According to various embodiments of the present application, the computer device 700 can also be connected to a remote computer on the network through a network such as the Internet. That is, the computer device 700 can be connected to the network 712 through the network interface unit 711 connected to the system bus 705, or the network interface unit 711 can be used to connect to other types of networks or remote computer systems (not shown).

[0162] The memory further includes one or more programs, which are stored in the memory and configured to be executed by one or more processors. The one or more programs include instructions for implementing the volcanic rock microstructure pore characterization method on the computer device side.

[0163] In an exemplary embodiment, a computer device is also provided, the computer device comprising a processor and a memory, wherein the memory stores at least one instruction, at least one program, code set or instruction set. The at least one instruction, at least one program, code set or instruction set is configured to be executed by one or more processors to implement the volcanic rock microstructure pore characterization method on the computer device side.

[0164] In an illustrative embodiment, a computer-readable storage medium is also provided, in which at least one instruction, at least one program, a code set or an instruction set is stored. When the at least one instruction, the at least one program, the code set or the instruction set is executed by a processor of a computer device, the method for characterizing the pores of the volcanic rock microstructure on the computer device side is implemented.

[0165] Optionally, the computer-readable storage medium may be a ROM (Read-Only Memory), a RAM (Random Access Memory), a CD-ROM (Compact Disc Read-Only Memory), a magnetic tape, a floppy disk, an optical data storage device, or the like.

[0166] In an exemplary embodiment, a computer program product or a computer program is also provided, the computer program product or the computer program comprising computer instructions, the computer instructions being stored in a computer-readable storage medium. A processor of a computer device reads the computer instructions from the computer-readable storage medium, and the processor executes the computer instructions, so that the computer device executes the above-mentioned volcanic rock microstructure pore characterization method on the computer device side.

[0167] The above description is only an exemplary embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method for characterizing pores in volcanic rock microstructure, characterized in that: The method comprises: Acquiring at least one pore image of a first rock sample of a volcanic rock core; determining a pore profile of a pore in at least one of the pore images; determining a pore shape of at least one of the pores based on a pore profile of at least one of the pores; determining an initial pore radius of at least one of the pores based on a pore shape of at least one of the pores; determining a pore radius range of the first rock sample based on an initial pore radius of at least one of the pores; Based on the pore radius range of the first rock sample, a pressure range corresponding to the pore radius range is determined; wherein the pressure value P corresponding to the pore radius r is determined by the following formula: c : Determining a target relationship curve corresponding to the pressure range from the relationship curve, wherein the relationship curve refers to a relationship curve between pressure and volume obtained by performing a constant-rate mercury injection experiment on a second rock sample of the volcanic rock core; Based on the target relationship curve, a standard relationship curve between pressure and volume is obtained; determining a target volume value for at least one of the pores from the standard relationship curve; determining an equivalent pore radius of at least one of the pores based on a target volume value of at least one of the pores; determining a radius determination method corresponding to a pore shape of at least one of the pores; Determining a final pore radius of at least one of the pores based on the radius determination method and the equivalent pore radius; Based on the final pore radius of at least one pore, the target pore radius of pores with the same pore shape is determined; the sum of the products of the shape distribution frequencies corresponding to the pores with the same pore shape and the target pore radius is determined as the final equivalent radius of the core, and the shape distribution frequency is used to indicate the ratio of the number of pores with the same pore shape to the total number of at least one pore, and the pore shapes include triangular pores, quadrilateral pores and pentagonal pores.

2. The method according to claim 1, characterized in that The determining the pore shape of at least one of the pores based on the pore profile of at least one of the pores comprises: determining a cross-sectional area of ​​at least one of the pores and a perimeter of at least one of the pores based on a pore profile of at least one of the pores; determining a shape factor of at least one of the pores based on a cross-sectional area of ​​at least one of the pores and a perimeter of at least one of the pores; Based on a shape factor of at least one of the pores, a pore shape of at least one of the pores is determined.

3. The method according to claim 1, characterized in that The determining of a final pore radius of at least one of the pores based on the radius determination method and the equivalent pore radius comprises: In the case where the pore shape of the pore is a triangle, determining the equivalent pore radius of the first value as the final pore radius; In the case where the pore shape of the pore is a quadrilateral, determining the equivalent pore radius of the second value as the final pore radius; In a case where the pore shape of the pore is a pentagon, the equivalent pore radius of the third value is determined as the final pore radius.

4. The method according to claim 1, characterized in that: After determining the pore shape of at least one of the pores based on the shape factor of at least one of the pores, the method further includes: Determining the number of pores having the same pore shape based on the pore shape of at least one of the pores; Based on the number of pores with the same type of pore shape and the total number of at least one of the pores, shape distribution frequencies corresponding to pores with different types of pore shapes are determined.

5. A volcanic rock microstructure pore characterization device, characterized in that: The device comprises: An image acquisition module, configured to acquire at least one pore image of a first rock sample of a volcanic rock core; a contour determination module, configured to determine a pore contour of a pore in at least one of the pore images; a shape determination module for determining a pore shape of at least one of the pores based on a pore profile of at least one of the pores; a radius determination module, configured to determine an initial pore radius of at least one of the pores based on a pore shape of at least one of the pores; A range determination module is used to determine the pore radius range of the first rock sample based on the initial pore radius of at least one of the pores; wherein the pressure value P corresponding to the pore radius r is determined by the following formula c : The range determination module is further used to determine a pressure range corresponding to the pore radius range based on the pore radius range of the first rock sample; a curve determination module, used to determine a target relationship curve corresponding to the pressure range from a relationship curve, wherein the relationship curve refers to a relationship curve between pressure and volume obtained by performing a constant-rate mercury injection experiment on a second rock sample of the volcanic rock core; The curve determination module is further used to obtain a standard relationship curve between pressure and volume based on the target relationship curve; a volume determination module, configured to determine a target volume value of at least one of the pores from the standard relationship curve; The radius determination module is further used to determine an equivalent pore radius of at least one of the pores based on a target volume value of at least one of the pores; The radius determination module is further used to determine a radius determination method corresponding to the pore shape of at least one of the pores; based on the radius determination method and the equivalent pore radius, determine a final pore radius of at least one of the pores; The device is also used to: determine the target pore radius of pores with the same type of pore shape based on the final pore radius of at least one pore; determine the sum of the products of the shape distribution frequencies corresponding to the pores with the same type of pore shape and the target pore radius as the final equivalent radius of the core, the shape distribution frequency is used to indicate the ratio of the number of pores with the same type of pore shape to the total number of at least one pore, and the pore shapes include triangular pores, quadrilateral pores and pentagonal pores.

6. A computer device, characterized in that: The computer device includes a processor and a memory, wherein the memory stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, the at least one program, the code set or the instruction set is loaded and executed by the processor to implement the volcanic rock microstructure pore characterization method according to any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: The storage medium stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, the at least one program, the code set or the instruction set is loaded and executed by the processor to implement the volcanic rock microstructure pore characterization method according to any one of claims 1 to 4.

Images