A rock permeability calculation method, device, equipment and medium

Abstract

The present application discloses a rock permeability calculation method, device, equipment and medium, which relates to the field of geomechanics, including: obtaining permeability calculation parameters; wherein the permeability calculation parameters include initial porosity, normal strain parameter and initial permeability; performing apparent strain calculation on the permeability calculation parameters to obtain apparent area strain parameters, and calculating the apparent area strain of each rock and the apparent area strain of each rock skeleton based on the apparent area strain parameters; calculating the porosity of each rock after deformation based on the apparent area strain of each rock and the apparent area strain of each rock skeleton; obtaining capillary parameters after rock deformation, and calculating the permeability of each rock after deformation based on the permeability calculation parameters, the capillary parameters after rock deformation and the porosity of each rock. Through the above technical solution of the present application, the permeability change in different directions can be solved according to the strain field, and the dynamic permeability model can be used to realize the numerical simulation of fluid-solid coupling of anisotropic formations.

Description

Technical Field

[0001] The present invention relates to the field of geomechanics, and in particular to a rock permeability calculation method, device, equipment and medium. Background Art

[0002] Because the arrangement of rock particles during sedimentation and compaction is directional, with the long axes of the particles primarily parallel to the direction of water flow, rocks exhibit anisotropy. Permeability anisotropy has a significant impact on geological activities such as oil and gas field development, groundwater flow, and earthquakes. Therefore, anisotropic permeability has been widely accepted and applied. Changes in underground pore pressure during processes such as underground fluid injection and production can cause changes in the formation stress state and rock deformation. Rock deformation can cause changes in parameters such as rock porosity and permeability, which in turn couples and influences the seepage process. Permeability is the bridge that connects the deformation field and the seepage field. Only by establishing a relationship between permeability and stress or strain fields (i.e., a dynamic permeability model) can the influence of the deformation field on the seepage field be reflected and true seepage coupling analysis be achieved. However, currently mainstream dynamic permeability models do not consider anisotropy.

[0003] From the above, it can be seen that how to solve the permeability changes in different directions based on the strain field and use the dynamic permeability model to realize the numerical simulation of fluid-solid coupling in anisotropic formations is a problem to be solved in this field. Summary of the Invention

[0004] In view of this, the present invention aims to provide a rock permeability calculation method, device, equipment, and medium that can solve the permeability variation in different directions based on the strain field and realize the numerical simulation of fluid-solid coupling in anisotropic formations using a dynamic permeability model. The specific scheme is as follows:

[0005] In a first aspect, the present application discloses a rock permeability calculation method, which is applied to a dynamic permeability calculation model, comprising:

[0006] Obtaining permeability calculation parameters; wherein the permeability calculation parameters include initial porosity, normal strain parameter and initial permeability;

[0007] Performing apparent strain calculation on the permeability calculation parameters to obtain apparent area strain parameters, and calculating the apparent area strain of each rock and the apparent area strain of each rock skeleton based on the apparent area strain parameters;

[0008] Calculating the porosity of each rock after deformation based on the apparent area strain of each rock and the apparent area strain of each rock skeleton;

[0009] The capillary parameters of the rock after deformation are obtained, and the permeabilities of the rocks after deformation are calculated based on the permeability calculation parameters, the capillary parameters of the rock after deformation and the surface porosity of each rock.

[0010] Optionally, before obtaining the permeability calculation parameters, the method further includes:

[0011] The initial dynamic permeability calculation model is configured and analyzed based on the capillary network physical model and according to preset model configuration rules to obtain the dynamic permeability calculation model.

[0012] Optionally, performing apparent strain calculation on the permeability calculation parameter to obtain an apparent area strain parameter includes:

[0013] Determining the rock deformation direction; wherein the rock deformation direction includes a horizontal axis direction, a longitudinal axis direction, and a vertical axis direction;

[0014] Apparent strain calculations are performed on the permeability calculation parameters according to the rock deformation direction to obtain a horizontal axis apparent area strain parameter, a vertical axis apparent area strain parameter, and a vertical axis apparent area strain parameter.

[0015] Optionally, the calculating the apparent area strain of each rock and the apparent area strain of each rock skeleton according to the apparent area strain parameter includes:

[0016] Calculate the apparent area strain of the rock in the transverse-longitudinal section, the apparent area strain of the rock in the transverse-vertical section, and the apparent area strain of the rock in the longitudinal-vertical section according to the apparent area strain parameter of the transverse axis, the apparent area strain parameter of the longitudinal axis, and the apparent area strain parameter of the vertical axis;

[0017] The apparent area strain of the rock skeleton in the transverse-longitudinal section, the apparent area strain of the rock skeleton in the transverse-vertical section and the apparent area strain of the rock skeleton in the longitudinal-vertical section are calculated according to the apparent area strain parameter of the transverse axis, the apparent area strain parameter of the longitudinal axis and the apparent area strain parameter of the vertical axis.

[0018] Optionally, the calculating of each surface porosity of the rock after deformation based on each rock apparent area strain and each rock skeleton apparent area strain respectively includes:

[0019] Calculating the porosity of the rock in the transverse and longitudinal sections after deformation based on the apparent area strain of the rock in the transverse and longitudinal sections and the apparent area strain of the rock skeleton in the transverse and longitudinal sections;

[0020] Calculating the porosity of the rock in the transverse-vertical section after deformation based on the apparent area strain of the rock in the transverse-vertical section and the apparent area strain of the rock skeleton in the transverse-vertical section;

[0021] The porosity of the longitudinal-vertical cross-section of the rock after deformation is calculated based on the apparent area strain of the longitudinal-vertical cross-section rock and the apparent area strain of the longitudinal-vertical cross-section rock skeleton.

[0022] Optionally, the calculating of the permeability of each deformed rock based on the permeability calculation parameter, the capillary parameter of the rock after deformation, and each surface porosity includes:

[0023] Calculate the vertical axis permeability of the rock after deformation based on the initial porosity, the capillary parameter of the rock after deformation, and the surface porosity of the transverse and longitudinal sections in the permeability calculation parameters;

[0024] Calculate the vertical axis permeability of the rock after deformation based on the initial porosity, the capillary parameter of the rock after deformation, and the porosity of the horizontal-vertical cross section in the permeability calculation parameters;

[0025] The transverse axis permeability of the rock after deformation is calculated based on the initial porosity, the capillary parameter of the rock after deformation, and the surface porosity of the longitudinal-vertical cross section in the permeability calculation parameters.

[0026] Optionally, after calculating the permeability of each deformed rock based on the permeability calculation parameter, the capillary parameter of the rock after deformation, and each surface porosity, the method further includes:

[0027] The permeabilities of the rocks after deformation are sent to the client so that the client can analyze the permeabilities of the rocks after deformation to obtain the permeability changes.

[0028] In a second aspect, the present application discloses a rock permeability calculation device, comprising:

[0029] A parameter acquisition module is used to acquire permeability calculation parameters; wherein the permeability calculation parameters include initial porosity, normal strain parameter and initial permeability;

[0030] an apparent strain calculation module, configured to perform apparent strain calculation on the permeability calculation parameters to obtain apparent area strain parameters, and calculate the apparent area strain of each rock and the apparent area strain of each rock skeleton based on the apparent area strain parameters;

[0031] A surface porosity calculation module, configured to calculate the surface porosity of each rock after deformation based on the apparent area strain of each rock and the apparent area strain of each rock skeleton;

[0032] The module for calculating the permeability of deformed rock is used to obtain the capillary parameters of the rock after deformation, and calculate the permeability of each deformed rock based on the permeability calculation parameters, the capillary parameters of the rock after deformation and the surface porosity of each rock.

[0033] In a third aspect, the present application discloses an electronic device, comprising:

[0034] Memory, used to store computer programs;

[0035] A processor is used to execute the computer program to implement the aforementioned rock permeability calculation method.

[0036] In a fourth aspect, the present application discloses a computer storage medium for storing a computer program; wherein, when the computer program is executed by a processor, the steps of the rock permeability calculation method disclosed above are implemented.

[0037] It can be seen that the present application provides a method for calculating rock permeability, including obtaining permeability calculation parameters; wherein the permeability calculation parameters include initial porosity, normal strain parameters, and initial permeability; performing apparent strain calculation on the permeability calculation parameters to obtain apparent area strain parameters, and calculating the apparent area strain of each rock and the apparent area strain of each rock skeleton based on the apparent area strain parameters; respectively calculating the porosity of each rock after deformation based on the apparent area strain of each rock and the apparent area strain of each rock skeleton; obtaining the capillary parameters after rock deformation, and calculating the permeability of each deformed rock based on the permeability calculation parameters, the capillary parameters after rock deformation, and the porosity of each rock. The present application clarifies the physical relationship between the permeability tensor and the strain tensor, establishes a reasonable relationship between the principal value of the permeability tensor and the normal strain, and provides more realistic dynamic parameters for the coupling process of the deformation field and the seepage field. The calculation method of the present application can directly solve the permeability changes in different directions based on the strain field obtained by geomechanical calculation, providing a more reliable dynamic permeability model for the numerical simulation of fluid-solid coupling in anisotropic formations. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are merely embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0039] Figure 1 This is a flow chart of a rock permeability calculation method disclosed in this application;

[0040] Figure 2 This is a schematic diagram of rock area strain disclosed in this application;

[0041] Figure 3 This is a schematic structural diagram of a rock permeability calculation device disclosed in this application;

[0042] Figure 4 This is a structural diagram of an electronic device provided in this application. DETAILED DESCRIPTION

[0043] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0044] Because the arrangement of rock particles during sedimentation and compaction is directional, with the long axis of the particles primarily parallel to the direction of water flow, the rock exhibits anisotropy. Permeability anisotropy has a significant impact on geological activities such as oil and gas field development, groundwater flow, and earthquakes. Therefore, anisotropic permeability has been widely accepted and applied. Changes in underground pore pressure during processes such as underground fluid injection and production can cause changes in the formation stress state, causing rock deformation. Rock deformation can cause changes in parameters such as rock porosity and permeability, which in turn couples and affects the seepage process. Permeability is the bridge through which the deformation field and seepage field interact and couple. Only by establishing a relationship between permeability and stress or strain fields (i.e., a dynamic permeability model) can the influence of the deformation field on the seepage field be reflected and true seepage coupling analysis be achieved. However, current mainstream dynamic permeability models do not consider anisotropy. As can be seen from the above, how to solve the permeability changes in different directions based on the strain field and use dynamic permeability models to achieve numerical simulations of fluid-solid coupling in anisotropic formations is an unresolved problem in this field.

[0045] See also Figure 1 As shown, the embodiment of the present invention discloses a rock permeability calculation method, which is applied to a dynamic permeability calculation model and may specifically include:

[0046] Step S11: Obtain permeability calculation parameters; wherein the permeability calculation parameters include initial porosity, normal strain parameter and initial permeability.

[0047] In this embodiment, before obtaining the permeability calculation parameters, the method further includes: configuring and analyzing the initial dynamic permeability calculation model based on the capillary network physical model and according to the preset model configuration rules to obtain the dynamic permeability calculation model. Specifically, the initial dynamic permeability calculation model is first configured and analyzed, and the specific steps are as follows: (1) The principal value of the permeability tensor (K x ,K y ,K z ) and normal strain (ε x ,ε y ,ε z) in the same direction; (2) the analysis is based on the capillary network physical model; (3) the initial rock porosity is consistent with the surface porosity; (4) the strain follows the convention that tension is positive and compression is negative.

[0048] In addition, the process of determining the initial permeability in the permeability calculation parameters is as follows, based on the definition of permeability and considering that the initial rock porosity and surface area are consistent:

[0049]

[0050] Where φ0 is the initial rock porosity, is the initial permeability in the x direction.

[0051] The initial capillary radius is calculated using the above formula, and then the Poiseuille equation for a circular tube is used to calculate the initial capillary radius to obtain the single capillary flow rate. The cross-sectional surface area ratio is determined, and the initial permeability is calculated based on the single capillary flow rate and the cross-sectional surface area ratio to obtain the initial permeability. Taking the x-direction (i.e., the horizontal axis) as an example, the Poiseuille equation for a single circular tube is:

[0052]

[0053] Where r is the capillary radius, μ is the fluid viscosity, Δp is the pressure difference across the capillary, L is the capillary length, and q is the flow rate of a single capillary.

[0054] The number of capillaries per unit yz cross-sectional area (i.e. longitudinal-vertical) is:

[0055]

[0056] Among them, φ yz is the surface ratio of the yz section.

[0057] Step S12: performing apparent strain calculation on the permeability calculation parameters to obtain apparent area strain parameters, and calculating the apparent area strain of each rock and the apparent area strain of each rock skeleton based on the apparent area strain parameters.

[0058] In this embodiment, the rock deformation direction is determined; wherein the rock deformation direction includes the horizontal axis direction, the vertical axis direction, and the vertical axis direction; the apparent strain calculation is performed on the permeability calculation parameters according to the rock deformation direction to obtain the horizontal axis apparent area strain parameter, the vertical axis apparent area strain parameter, and the vertical axis apparent area strain parameter, and then the apparent area strain of each rock and the apparent area strain of each rock skeleton are calculated based on the apparent area strain parameters.

[0059] The specific process of calculating the apparent area strain of each rock and the apparent area strain of each rock skeleton is as follows: the apparent area strain of the rock in the transverse-longitudinal section, the apparent area strain of the rock in the transverse-vertical section, and the apparent area strain of the rock in the longitudinal-vertical section are calculated according to the apparent area strain parameter of the transverse axis, the apparent area strain parameter of the longitudinal axis, and the apparent area strain parameter of the vertical axis; the apparent area strain of the rock skeleton in the transverse-longitudinal section, the apparent area strain of the rock skeleton in the transverse-vertical section, and the apparent area strain of the rock skeleton in the longitudinal-vertical section are calculated according to the apparent area strain parameter of the transverse axis, the apparent area strain parameter of the longitudinal axis, and the apparent area strain parameter of the vertical axis.

[0060] ε x , ε y , ε z is the apparent area strain parameter of the horizontal axis, the apparent area strain parameter of the vertical axis and the apparent area strain parameter of the vertical axis, ε xy , ε xz , ε yz is the apparent area strain of the rock in the xy section, xz section, and yz section (i.e., the apparent area strain of the rock in the transverse-vertical section, the apparent area strain of the rock in the transverse-vertical section, and the apparent area strain of the rock in the longitudinal-vertical section). The specific calculation process is: xy =ε x +ε y ; ε xz =ε x +ε z ; ε yz =ε y +ε z , is the area strain of the rock skeleton in the xy section, is the area strain of the rock skeleton in the xz section, is the area strain of the rock skeleton in the yz section (i.e., the apparent area strain of the rock skeleton in the transverse-longitudinal section, the apparent area strain of the rock skeleton in the transverse-vertical section, and the apparent area strain of the rock skeleton in the longitudinal-vertical section). The specific calculation process is:

[0061] Step S13: Calculate the porosity of each rock after deformation based on the apparent area strain of each rock and the apparent area strain of each rock skeleton.

[0062] In this embodiment, the porosity of the transverse-longitudinal section of the rock after deformation is calculated based on the apparent area strain of the transverse-longitudinal section rock and the apparent area strain of the transverse-longitudinal section rock skeleton; the porosity of the transverse-vertical section of the rock after deformation is calculated based on the apparent area strain of the transverse-vertical section rock and the apparent area strain of the transverse-vertical section rock skeleton; the porosity of the longitudinal-vertical section of the rock after deformation is calculated based on the apparent area strain of the longitudinal-vertical section rock and the apparent area strain of the longitudinal-vertical section rock skeleton.

[0063] In this embodiment, the relationship between the surface ratio, capillary radius and strain is as follows: Figure 2 As shown, the face rate after deformation is:

[0064]

[0065] right Figure 2 In the case shown, the number of capillaries in the study area remains constant before and after deformation, and the ratio of the pore areas in the cross section is the square of the ratio of the capillary radii:

[0066]

[0067] Step S14: Obtain capillary parameters of the rock after deformation, and calculate the permeability of each deformed rock based on the permeability calculation parameters, the capillary parameters of the rock after deformation, and the surface porosity of each rock.

[0068] In this embodiment, after obtaining the capillary parameters of the rock after deformation, the vertical axis permeability of the deformed rock is calculated based on the initial porosity in the permeability calculation parameters, the capillary parameters of the rock after deformation, and the surface porosity of the transverse-vertical section; the vertical axis permeability of the deformed rock is calculated based on the initial porosity in the permeability calculation parameters, the capillary parameters of the rock after deformation, and the surface porosity of the transverse-vertical section; the transverse axis permeability of the deformed rock is calculated based on the initial porosity in the permeability calculation parameters, the capillary parameters of the rock after deformation, and the surface porosity of the transverse-vertical section.

[0069] In this embodiment, the permeability in the x direction after deformation is:

[0070]

[0071] Among them, R is the radius of the capillary tube after deformation, K x is the permeability in the x direction after deformation.

[0072] The ratio of the permeability in the x direction before and after deformation is:

[0073]

[0074] Substituting the above formula into the square of the capillary radius ratio, we can obtain the relationship between the permeability in the x-direction and the rock strain:

[0075]

[0076] Similarly, the relationship between the y direction (i.e., longitudinal axis direction) and the z direction (i.e., vertical axis direction) can be obtained as follows:

[0077]

[0078]

[0079] Substitute the initial capillary radius, initial porosity, capillary parameters after rock deformation, and vertical-vertical cross-section porosity into the above equation to obtain the horizontal axis permeability. Similarly, the vertical axis permeability and vertical axis permeability are obtained.

[0080] In this embodiment, after calculating the permeability of each deformed rock based on the permeability calculation parameters, the capillary parameters of the rock after deformation, and the surface porosity of each, it also includes: sending the permeability of each deformed rock to the client, so that the client can analyze the permeability of each deformed rock to obtain the permeability change.

[0081] For example, the initial porosity at a point in a saltwater layer is 0.26, and the initial permeability is:

[0082]

[0083] After the fluid is injected into the phase formation, the apparent strain at this point is calculated through geomechanics:

[0084] ε x =-3.9240×10 -4 , ε y =-3.9211×10 -4 , ε z =0.0020

[0085]

[0086] According to the method in this paper:

[0087] ε xy =ε x +ε y =-3.9240×10 -4 -3.9211×10 -4 =-7.8451×10 -4

[0088] ε xz =ε x +ε z =-3.9240×10 -4 +0.0020=1.6076×10 -3

[0089] ε yz =ε y +ε z =-3.9211×10 -4 +0.0020=1.60789×10 -3

[0090]

[0091]

[0092]

[0093]

[0094]

[0095]

[0096]

[0097] Through the examples given above, it can be seen that the permeability in different directions is not only inconsistent in the amplitude of change, but even the direction of change may be inconsistent. In the above case, the permeability in the z direction decreases, while the permeability in the x and y directions increases. If the traditional method is used to calculate only based on volumetric strain, the wrong conclusion that all permeabilities increase will be obtained. Therefore, the calculation method of the present application can directly solve the permeability variation in different directions based on the strain field obtained by geomechanical calculation, providing a more reliable dynamic permeability model for the numerical simulation of fluid-solid coupling in anisotropic strata.

[0098] In this embodiment, permeability calculation parameters are obtained; wherein, the permeability calculation parameters include initial porosity, normal strain parameters, and initial permeability; apparent strain calculation is performed on the permeability calculation parameters to obtain apparent area strain parameters, and the apparent area strain of each rock and the apparent area strain of each rock skeleton are calculated based on the apparent area strain parameters; the porosity of each rock after deformation is calculated based on the apparent area strain of each rock and the apparent area strain of each rock skeleton; the capillary parameters after rock deformation are obtained, and the permeability of each deformed rock is calculated based on the permeability calculation parameters, the capillary parameters after rock deformation, and the porosity of each rock. This application clarifies the physical relationship between the permeability tensor and the strain tensor, establishes a reasonable relationship between the principal value of the permeability tensor and the normal strain, and provides more realistic dynamic parameters for the coupling process of the deformation field and the seepage field. The calculation method of this application can be used to directly solve the permeability changes in different directions based on the strain field obtained by geomechanical calculation, providing a more reliable dynamic permeability model for the numerical simulation of fluid-solid coupling in anisotropic formations.

[0099] See also Figure 3 As shown, the embodiment of the present invention discloses a rock permeability calculation device, which may specifically include:

[0100] The parameter acquisition module 11 is used to acquire permeability calculation parameters; wherein the permeability calculation parameters include initial porosity, normal strain parameter and initial permeability;

[0101] An apparent strain calculation module 12 is used to perform apparent strain calculation on the permeability calculation parameters to obtain apparent area strain parameters, and calculate the apparent area strain of each rock and the apparent area strain of each rock skeleton based on the apparent area strain parameters;

[0102] A surface porosity calculation module 13 is configured to calculate the surface porosity of each rock after deformation based on the apparent area strain of each rock and the apparent area strain of each rock skeleton;

[0103] The deformed rock permeability calculation module 14 is used to obtain the capillary parameters of the rock after deformation, and calculate the permeability of each deformed rock based on the permeability calculation parameters, the capillary parameters of the rock after deformation and the surface porosity of each rock.

[0104] In this embodiment, permeability calculation parameters are obtained; wherein, the permeability calculation parameters include initial porosity, normal strain parameters, and initial permeability; apparent strain calculation is performed on the permeability calculation parameters to obtain apparent area strain parameters, and the apparent area strain of each rock and the apparent area strain of each rock skeleton are calculated based on the apparent area strain parameters; the porosity of each rock after deformation is calculated based on the apparent area strain of each rock and the apparent area strain of each rock skeleton; the capillary parameters after rock deformation are obtained, and the permeability of each deformed rock is calculated based on the permeability calculation parameters, the capillary parameters after rock deformation, and the porosity of each rock. This application clarifies the physical relationship between the permeability tensor and the strain tensor, establishes a reasonable relationship between the principal value of the permeability tensor and the normal strain, and provides more realistic dynamic parameters for the coupling process of the deformation field and the seepage field. The calculation method of this application can be used to directly solve the permeability changes in different directions based on the strain field obtained by geomechanical calculation, providing a more reliable dynamic permeability model for the numerical simulation of fluid-solid coupling in anisotropic formations.

[0105] In some specific embodiments, the parameter acquisition module 11 may specifically include:

[0106] The configuration analysis module is used to configure and analyze the initial dynamic permeability calculation model based on the capillary network physical model and according to preset model configuration rules to obtain the dynamic permeability calculation model.

[0107] In some specific embodiments, the apparent strain calculation module 12 may specifically include:

[0108] A direction determination module is used to determine the rock deformation direction; wherein the rock deformation direction includes a horizontal axis direction, a longitudinal axis direction, and a vertical axis direction;

[0109] The apparent strain calculation module is used to perform apparent strain calculation on the permeability calculation parameters according to the rock deformation direction to obtain the horizontal axis apparent area strain parameter, the vertical axis apparent area strain parameter and the vertical axis apparent area strain parameter.

[0110] In some specific embodiments, the apparent strain calculation module 12 may specifically include:

[0111] a rock apparent area strain calculation module, configured to calculate the rock apparent area strain in the transverse-longitudinal section, the rock apparent area strain in the transverse-vertical section, and the rock apparent area strain in the longitudinal-vertical section based on the transverse axis apparent area strain parameter, the longitudinal axis apparent area strain parameter, and the vertical axis apparent area strain parameter;

[0112] The rock skeleton apparent area strain calculation module is used to calculate the apparent area strain of the rock skeleton in the transverse-longitudinal section, the apparent area strain of the rock skeleton in the transverse-vertical section, and the apparent area strain of the rock skeleton in the longitudinal-vertical section based on the apparent area strain parameter of the transverse axis, the apparent area strain parameter of the longitudinal axis, and the apparent area strain parameter of the vertical axis.

[0113] In some specific embodiments, the face-changing face rate calculation module 13 may specifically include:

[0114] A first deformed porosity calculation module is used to calculate the porosity of the rock after deformation in the transverse and longitudinal sections based on the apparent area strain of the rock in the transverse and longitudinal sections and the apparent area strain of the rock skeleton in the transverse and longitudinal sections;

[0115] A second deformation porosity calculation module is used to calculate the porosity of the rock after deformation in the transverse-vertical section based on the apparent area strain of the rock in the transverse-vertical section and the apparent area strain of the rock skeleton in the transverse-vertical section;

[0116] The third deformed surface porosity calculation module is used to calculate the surface porosity of the longitudinal-vertical section of the rock after deformation based on the longitudinal-vertical section rock apparent area strain and the longitudinal-vertical section rock skeleton apparent area strain.

[0117] In some specific embodiments, the deformed rock permeability calculation module 14 may specifically include:

[0118] A first rock permeability calculation module is used to calculate the vertical axis permeability of the rock after deformation based on the initial porosity, the capillary parameter of the rock after deformation, and the surface porosity of the transverse and longitudinal sections in the permeability calculation parameters;

[0119] A second rock permeability calculation module is used to calculate the longitudinal axis permeability of the rock after deformation based on the initial porosity in the permeability calculation parameters, the capillary parameter of the rock after deformation, and the porosity of the horizontal-vertical cross section;

[0120] The third rock permeability calculation module is used to calculate the transverse axis permeability of the deformed rock based on the initial porosity, the capillary parameters of the rock after deformation, and the vertical-vertical cross-section porosity in the permeability calculation parameters.

[0121] In some specific embodiments, the deformed rock permeability calculation module 14 may specifically include:

[0122] The permeability change determination module is used to send each of the dynamic permeability ratios to the client so that the client can analyze the dynamic permeability ratios to obtain the permeability change.

[0123] Figure 4 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. The electronic device 20 may include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is used to store a computer program, which is loaded and executed by the processor 21 to implement the relevant steps of the rock permeability calculation method performed by the electronic device as disclosed in any of the aforementioned embodiments.

[0124] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and the external device. The communication protocol it follows is any communication protocol that can be applied to the technical solution of this application and is not specifically limited here; the input and output interface 25 is used to obtain external input data or output data to the outside world. Its specific interface type can be selected according to specific application needs and is not specifically limited here.

[0125] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or CD, etc. The resources stored thereon include an operating system 221, a computer program 222 and data 223, etc. The storage method can be temporary storage or permanent storage.

[0126] The operating system 221 is used to manage and control the hardware devices and computer program 222 on the electronic device 20, enabling the processor 21 to calculate and process data 223 in the memory 22. The operating system 221 can be Windows, Unix, Linux, etc. In addition to including computer programs capable of performing the rock permeability calculation method performed by the electronic device 20 as disclosed in any of the aforementioned embodiments, the computer program 222 may further include computer programs capable of performing other specific tasks. Data 223 may include data transmitted by external devices and received by the rock permeability calculation device, as well as data collected by its own input / output interface 25.

[0127] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in a random access memory (RAM), internal memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0128] Furthermore, an embodiment of the present application also discloses a computer-readable storage medium, in which a computer program is stored. When the computer program is loaded and executed by a processor, the steps of the rock permeability calculation method disclosed in any of the aforementioned embodiments are implemented.

[0129] Finally, it should be noted that, in this document, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or device. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article, or device comprising the element.

[0130] The above is a detailed introduction to the rock permeability calculation method, device, equipment and storage medium provided by the present invention. Specific examples are used herein to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. At the same time, for those skilled in the art, according to the ideas of the present invention, there may be changes in the specific implementation methods and application scopes. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A rock permeability calculation method, characterized in that: Applied to dynamic permeability calculation models, including: Obtaining permeability calculation parameters; wherein the permeability calculation parameters include initial porosity, normal strain parameter and initial permeability; Performing apparent strain calculation on the permeability calculation parameters to obtain apparent area strain parameters, and calculating the apparent area strain of each rock and the apparent area strain of each rock skeleton based on the apparent area strain parameters; Calculating the porosity of each rock after deformation based on the apparent area strain of each rock and the apparent area strain of each rock skeleton; The capillary parameters of the rock after deformation are obtained, and the permeabilities of the rocks after deformation are calculated based on the permeability calculation parameters, the capillary parameters of the rock after deformation and the surface porosity of each rock.

2. The rock permeability calculation method according to claim 1, characterized in that: Before obtaining the permeability calculation parameters, the method further includes: The initial dynamic permeability calculation model is configured and analyzed based on the capillary network physical model and according to preset model configuration rules to obtain the dynamic permeability calculation model.

3. The rock permeability calculation method according to claim 1, characterized in that: The performing of apparent strain calculation on the permeability calculation parameter to obtain an apparent area strain parameter includes: Determining the rock deformation direction; wherein the rock deformation direction includes a horizontal axis direction, a longitudinal axis direction, and a vertical axis direction; Apparent strain calculations are performed on the permeability calculation parameters according to the rock deformation direction to obtain a horizontal axis apparent area strain parameter, a vertical axis apparent area strain parameter, and a vertical axis apparent area strain parameter.

4. The rock permeability calculation method according to claim 3, characterized in that: The calculating of the apparent area strain of each rock and the apparent area strain of each rock skeleton according to the apparent area strain parameter includes: Calculate the apparent area strain of the rock in the transverse-longitudinal section, the apparent area strain of the rock in the transverse-vertical section, and the apparent area strain of the rock in the longitudinal-vertical section according to the apparent area strain parameter of the transverse axis, the apparent area strain parameter of the longitudinal axis, and the apparent area strain parameter of the vertical axis; The apparent area strain of the rock skeleton in the transverse-longitudinal section, the apparent area strain of the rock skeleton in the transverse-vertical section and the apparent area strain of the rock skeleton in the longitudinal-vertical section are calculated according to the apparent area strain parameter of the transverse axis, the apparent area strain parameter of the longitudinal axis and the apparent area strain parameter of the vertical axis.

5. The rock permeability calculation method according to claim 4, characterized in that: The method of calculating the porosity of each rock after deformation based on the apparent area strain of each rock and the apparent area strain of each rock skeleton comprises: Calculating the porosity of the rock in the transverse and longitudinal sections after deformation based on the apparent area strain of the rock in the transverse and longitudinal sections and the apparent area strain of the rock skeleton in the transverse and longitudinal sections; Calculating the porosity of the rock in the transverse-vertical section after deformation based on the apparent area strain of the rock in the transverse-vertical section and the apparent area strain of the rock skeleton in the transverse-vertical section; The porosity of the longitudinal-vertical cross-section of the rock after deformation is calculated based on the apparent area strain of the longitudinal-vertical cross-section rock and the apparent area strain of the longitudinal-vertical cross-section rock skeleton.

6. The rock permeability calculation method according to claim 5, characterized in that: The step of calculating the permeability of each deformed rock based on the permeability calculation parameter, the capillary parameter of the rock after deformation, and each surface porosity includes: Calculate the vertical axis permeability of the rock after deformation based on the initial porosity, the capillary parameter of the rock after deformation, and the surface porosity of the transverse and longitudinal sections in the permeability calculation parameters; Calculate the vertical axis permeability of the rock after deformation based on the initial porosity, the capillary parameter of the rock after deformation, and the porosity of the horizontal-vertical cross section in the permeability calculation parameters; The transverse axis permeability of the rock after deformation is calculated based on the initial porosity, the capillary parameter of the rock after deformation, and the surface porosity of the longitudinal-vertical cross section in the permeability calculation parameters.

7. The rock permeability calculation method according to any one of claims 1 to 6, characterized in that: After calculating the permeability of each deformed rock based on the permeability calculation parameter, the capillary parameter of the rock after deformation, and each surface porosity, the method further includes: The permeabilities of the rocks after deformation are sent to the client so that the client can analyze the permeabilities of the rocks after deformation to obtain the permeability changes.

8. A rock permeability calculation device, characterized in that: include: A parameter acquisition module is used to acquire permeability calculation parameters; wherein the permeability calculation parameters include initial porosity, normal strain parameter and initial permeability; an apparent strain calculation module, configured to perform apparent strain calculation on the permeability calculation parameters to obtain apparent area strain parameters, and calculate the apparent area strain of each rock and the apparent area strain of each rock skeleton based on the apparent area strain parameters; A surface porosity calculation module, configured to calculate the surface porosity of each rock after deformation based on the apparent area strain of each rock and the apparent area strain of each rock skeleton; The module for calculating the permeability of deformed rock is used to obtain the capillary parameters of the rock after deformation, and calculate the permeability of each deformed rock based on the permeability calculation parameters, the capillary parameters of the rock after deformation and the surface porosity of each rock.

9. An electronic device, characterized in that: include: Memory, used to store computer programs; A processor, configured to execute the computer program to implement the rock permeability calculation method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that Used to store a computer program; wherein, when the computer program is executed by a processor, the rock permeability calculation method according to any one of claims 1 to 7 is implemented.

Images