A geomechanical simulation method, device, equipment and storage medium

Abstract

The application discloses a geomechanical simulation method and device, equipment and a storage medium, and relates to the technical field of oil and gas field development. The method comprises the following steps: acquiring basic application parameters of a porous medium in a target region in a process of underground fluid injection and production; determining equivalent volume force and equivalent area force when the porous medium produces elastic medium deformation in the process of underground fluid injection and production according to the basic application parameters; performing convolution solution on the equivalent volume force and the equivalent area force based on a Mindlin basic solution to obtain a deformation characteristic parameter; and guiding an oil and gas reservoir development strategy in the process of underground fluid injection and production according to the deformation characteristic parameter. Through the technical scheme, the deformation characteristic parameter determined in the process of underground fluid injection and production is not affected by the grid range and the degree of detail outside the reservoir, the calculation advantage of the far field such as the ground surface is obvious, and the deformation characteristic parameter is not limited by the shape of the reservoir and the pressure distribution state. The reservoir of any shape and the internal pressure of any distribution state can be simulated.

Description

Technical Field

[0001] This invention relates to the field of oil and gas field development technology, and in particular to a geomechanical simulation method, apparatus, equipment and storage medium. Background Technology

[0002] Subsurface fluid injection and production can easily induce corresponding mechanical problems, such as surface uplift and deformation, damage to caprock integrity, and fault activation. Previous studies have extensively investigated reservoir deformation during fluid injection and production, but currently, only near-field and far-field numerical simulation methods can accurately predict surface deformation in reservoirs with complex shapes. The target reservoirs for subsurface fluid injection and production projects such as carbon dioxide geological storage are generally buried at depths below 500m. The thickness of the strata from the reservoir to the surface is hundreds or even thousands of times greater than the reservoir thickness. To accurately simulate surface deformation, the model must include a considerable area of ​​the platform below the reservoir (at a depth where the rebound effect is negligible) and all rock layers from the reservoir to the surface. For such a large model, finely meshing the entire model as in traditional reservoir flow simulations would be impossible.

[0003] Existing research generally achieves this by refining the near-field reservoir grid and coarsening the far-field grid. The near-field and far-field numerical simulation method, represented by the coupling of ECLIPSE (a reservoir numerical simulation software) and VISAGE (a geomechanical simulation software), achieves the solution by refining the near-field reservoir seepage calculation grid and coarsening the far-field mechanical calculation grid. However, this method suffers from large far-field grids and insufficient accuracy in far-field calculations, such as those at the surface. Furthermore, due to the limitation of finite grids, it cannot fully consider the problems of semi-infinite strata, thus limiting its practical guidance.

[0004] Therefore, how to provide a solution to the above-mentioned technical problems is a problem that needs to be solved by those skilled in the art. Summary of the Invention

[0005] In view of this, the purpose of this invention is to provide a geomechanical simulation method, apparatus, equipment, and storage medium that can perform geomechanical simulations during underground fluid injection and extraction, ensuring that the calculation results are not affected by the range and fineness of the grid outside the reservoir, and accurately predicting the ground subsidence during the underground fluid injection and extraction process. The specific solution is as follows:

[0006] In a first aspect, this application discloses a geomechanical simulation method, including:

[0007] Obtain basic application parameters of porous media during the underground fluid injection and production process in the target area;

[0008] The equivalent volume force and equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process are determined using the aforementioned basic application parameters.

[0009] Based on the Mindlin fundamental solution, the equivalent volume force and the equivalent area force are convolved to obtain the deformation characteristic parameters;

[0010] The deformation characteristic parameters are used to guide the oil and gas reservoir development strategy during the underground fluid injection and production process.

[0011] Optionally, the acquisition of basic application parameters of porous media during the underground fluid injection and production process in the target area includes:

[0012] The bulk modulus of the skeleton material, the apparent bulk modulus of the rock, the pore fluid pressure of the porous medium, the porosity of the rock, the density of the rock skeleton material, and the average density and depth data of the fluid mixture in the pores are obtained during the underground fluid injection and extraction process in the target area.

[0013] Optionally, the equivalent volumetric force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined using the basic application parameters, including:

[0014] Using a preset equivalent volume force formula, the equivalent volume force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined by the basic application parameters.

[0015] The preset equivalent volume force formula is:

[0016]

[0017] Among them, K s K is the bulk modulus of the skeleton material. p ρ is the apparent bulk modulus of the rock, p is the pore fluid pressure, and ρ s ρ is the density of the rock skeleton material. l The average density of the fluid mixture in the pores is given, D is the depth data, g is the gravitational acceleration, and i represents the different axis directions of the three-dimensional Cartesian coordinate system.

[0018] Optionally, the equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined using the basic application parameters, including:

[0019] Using a preset equivalent area force formula, the equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined by the basic application parameters.

[0020] The formula for the equivalent area force is:

[0021] Among them, K s K is the bulk modulus of the skeleton material.p Let Ω represent the apparent bulk modulus of the rock, and let Ω represent the interior of the closed boundary of the target region. n represents the outside of the closed boundary of the target region. Ω Let be the unit in-line normal vector inside the closed boundary at the closed boundary Γ.

[0022] Optionally, the step of convolving the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain deformation characteristic parameters includes:

[0023] Based on the Mindlin fundamental solution, the equivalent volume force and the equivalent area force are convolved to obtain the stress, strain and displacement of the target region.

[0024] Optionally, the step of convolving the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain deformation characteristic parameters includes:

[0025] The stress is determined by integrating the pre-set stress determination formula using the Mindlin fundamental solution, the equivalent volume force, and the equivalent area force, based on the superposition principle. The pre-set stress determination formula is as follows: in, The stress Green's function, determined based on the Mindlin fundamental solution, is used to determine the stress at point x = (x1, x2, x3) under stress x. k The unit force applied in the direction causes a Δσ at point ω = (ω1, ω2, ω3). ij ;

[0026] The strain is determined by integrating the Mündlin fundamental solution, the equivalent volume force, and the equivalent area force using a preset strain determination formula, based on the superposition principle. The preset strain determination formula is as follows: in, The strain Green's function, determined based on the Mindlin fundamental solution, is used to determine the strain at point x = (x1, x2, x3) under stress x. k The unit force in the direction causes a Δε at point ω = (ω1, ω2, ω3). ij ;

[0027] The displacement is determined by integrating the preset displacement determination formula using the Mindlin fundamental solution, the equivalent volume force, and the equivalent area force, based on the superposition principle. The preset displacement determination formula is as follows: in, The displacement Green's function, determined based on the Mindlin fundamental solution, is used to determine the stress x at the point x = (x1, x2, x3). i A unit force acting in a direction causes x at point ω = (ω1, ω2, ω3) to be... j Displacement in the direction;

[0028] in, Δp=p t＝T -p t ＝0 ;

[0029] Optionally, after determining the equivalent volume force and equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process using the basic application parameters, the method further includes:

[0030] The equivalent volume force and the equivalent area force are discretized.

[0031] Accordingly, the convolution solution based on the Mindlin fundamental solution to obtain the deformation characteristic parameters includes:

[0032] Based on the fundamental solution of Mindlin, the equivalent volume force and the equivalent area force after discretization are convolved to obtain the deformation characteristic parameters after discretization.

[0033] Secondly, this application discloses a geomechanical simulation device, comprising:

[0034] The parameter acquisition module is used to acquire the basic application parameters of porous media during the underground fluid injection and extraction process in the target area.

[0035] The equivalent force determination module is used to determine the equivalent volume force and equivalent area force when the porous medium undergoes elastic media deformation during the underground fluid injection and extraction process, based on the basic application parameters.

[0036] The deformation characteristic parameter determination module is used to perform convolution between the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain the deformation characteristic parameters.

[0037] The strategy guidance module is used to guide the oil and gas reservoir development strategy during the underground fluid injection and production process based on the deformation characteristic parameters.

[0038] Thirdly, this application discloses an electronic device comprising a processor and a memory; wherein the memory is used to store a computer program, which is loaded and executed by the processor to implement the geomechanical simulation method as described above.

[0039] Fourthly, this application discloses a computer-readable storage medium for storing a computer program; wherein the computer program, when executed by a processor, implements the geomechanical simulation method as described above.

[0040] This application provides a geomechanical simulation method, comprising: acquiring basic application parameters of porous media during underground fluid injection and production in a target area; determining, through the basic application parameters, the equivalent volume force and equivalent area force when the porous media undergoes elastic deformation during the underground fluid injection and production process; performing convolution on the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain deformation characteristic parameters; and guiding oil and gas reservoir development strategies during the underground fluid injection and production process based on the deformation characteristic parameters. It can be seen that this application equates the mechanical changes of porous media during underground fluid injection and production to a continuous medium problem subjected to certain volume forces and area forces. In other words, by using equivalent volume forces and equivalent area forces, the rock deformation during underground fluid injection and production is transformed into the deformation of a general elastic medium subjected to equivalent volume forces and equivalent area forces. A semi-analytical solution is then performed based on the Mindlin fundamental solution and convolution method, avoiding insufficient computational accuracy caused by excessively large far-field grids. This achieves accurate calculations across the entire spatial domain, making the calculation results unaffected by the range and fineness of the grid outside the reservoir. It exhibits significant advantages in far-field calculations, such as those at the surface, and is not limited by reservoir shape or pressure distribution. It can calculate reservoirs of any shape and with any internal pressure distribution. Furthermore, the calculated deformation characteristic parameters can be used to predict surface subsidence, guiding oil and gas reservoir development strategies during underground fluid injection and production.

[0041] In addition, the geomechanical simulation device, equipment and storage medium provided in this application correspond to the above-mentioned geomechanical simulation method and have the same effect. Attached Figure Description

[0042] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0043] Figure 1 This is a flowchart of a geomechanical simulation method disclosed in this application;

[0044] Figure 2 This is a flowchart of a specific geomechanical simulation method disclosed in this application;

[0045] Figure 3 This is a flowchart of a specific geomechanical simulation method disclosed in this application;

[0046] Figure 4 This is a flowchart of a specific geomechanical simulation method disclosed in this application;

[0047] Figure 5This is a schematic diagram of the stress-strain calculation results for a target analysis region disclosed in this application.

[0048] Figure 6 This is a schematic diagram of the structure of a geomechanical simulation device disclosed in this application;

[0049] Figure 7 This is a structural diagram of an electronic device disclosed in this application. Detailed Implementation

[0050] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0051] Currently, analytical studies on the mechanical problems caused by underground fluid injection and production only focus on reservoirs with regular morphology and uniform pressure drop. Numerical studies suffer from problems such as large far-field grids, low computational accuracy, and inability to fully consider semi-infinite strata, thus having limited guiding significance for practical applications.

[0052] Therefore, this application provides a geomechanical simulation scheme that can perform geomechanical simulation during the underground fluid injection and extraction process, so that the calculation results are not affected by the range and fineness of the grid outside the reservoir, and accurately predict the ground subsidence during the underground fluid injection and extraction process.

[0053] This invention discloses a geomechanical simulation method, see [link to relevant documentation]. Figure 1 As shown, the method includes:

[0054] Step S11: Obtain the basic application parameters of porous media during the underground fluid injection and extraction process in the target area.

[0055] In this embodiment of the application, in order to perform high-precision geomechanical simulation of the underground fluid injection and extraction process, the rock deformation during the underground fluid injection and extraction process is transformed into the deformation of a general elastic medium subjected to equivalent force based on the principle of equivalent force. Therefore, some basic application parameters of the porous medium in the target area during the underground fluid injection and extraction process are first obtained to calculate the equivalent force.

[0056] For example, the acquisition of basic application parameters of porous media during the underground fluid injection and extraction process in the target area may include: acquiring the bulk modulus of the skeleton material, the apparent bulk modulus of the rock, the pore fluid pressure of the porous media, the rock porosity, the density of the rock skeleton material, and the average density and depth data of the fluid mixture in the pores during the underground fluid injection and extraction process in the target area.

[0057] Step S12: Determine the equivalent volume force and equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process using the basic application parameters.

[0058] In this embodiment, after obtaining the basic application parameters of the porous medium, the equivalent forces in the underground fluid injection and production process can be calculated, namely, the equivalent volume force and the equivalent area force. Thus, during underground fluid injection and production, such as in the development of oil and gas reservoirs, the mechanical changes of saturated porous media can be equated to a continuous medium problem subjected to certain volume forces and area forces. The deformation law of the porous medium skeleton is equivalent to general elastic deformation under the action of equivalent volume forces, and the influence of the reservoir sealing boundary on formation deformation is equivalent to general elastic deformation under the application of equivalent area forces on the virtual surface of the sealing boundary.

[0059] Step S13: Based on the Mindlin fundamental solution, perform convolution between the equivalent volume force and the equivalent area force to obtain the deformation characteristic parameters.

[0060] This application provides a novel high-precision geomechanical simulation method. Based on Mindlin's fundamental solution and convolution method, it performs semi-analytical solving, avoiding insufficient computational accuracy caused by excessively large far-field grids and achieving accurate calculation across the entire spatial domain. For example... Figure 2 The diagram shows the overall implementation process, in which the equivalent volume force and the equivalent area force are convolved to obtain the deformation characteristic parameters as the stress, strain and displacement of the target region.

[0061] Step S14: Guide the oil and gas reservoir development strategy during the underground fluid injection and production process based on the deformation characteristic parameters.

[0062] In this embodiment, the deformation characteristic parameters obtained are mainly stress, strain and displacement. These parameters can be used to perform fluid-structure interaction calculations during geomechanical simulation, thereby guiding the development strategy of oil and gas reservoirs during underground fluid injection and production.

[0063] It should be noted that the difficulty in actual reservoir research during underground fluid injection and production lies in the irregular reservoir shape and uneven pressure changes. Actual pore pressure variations can only be obtained through discrete numerical simulation; therefore, numerical calculations are essential for actual reservoir research. Consequently, it is necessary to discretize the parameters in the overall process. Specifically, the equivalent volume force and the equivalent area force are discretized; correspondingly, the convolution of the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain deformation characteristic parameters includes: convolving the discretized equivalent volume force and the discretized equivalent area force based on the Mindlin fundamental solution to obtain the discretized deformation characteristic parameters.

[0064] This application provides a geomechanical simulation method, comprising: acquiring basic application parameters of porous media during underground fluid injection and production in a target area; determining, through the basic application parameters, the equivalent volume force and equivalent area force when the porous media undergoes elastic deformation during the underground fluid injection and production process; performing convolution on the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain deformation characteristic parameters; and guiding oil and gas reservoir development strategies during the underground fluid injection and production process based on the deformation characteristic parameters. It can be seen that this application equates the mechanical changes of porous media during underground fluid injection and production to a continuous medium problem subjected to certain volume forces and area forces. In other words, by using equivalent volume forces and equivalent area forces, the rock deformation during underground fluid injection and production is transformed into the deformation of a general elastic medium subjected to equivalent volume forces and equivalent area forces. A semi-analytical solution is then performed based on the Mindlin fundamental solution and convolution method, avoiding insufficient computational accuracy caused by excessively large far-field grids. This achieves accurate calculations across the entire spatial domain, making the calculation results unaffected by the range and fineness of the grid outside the reservoir. It exhibits significant advantages in far-field calculations, such as those at the surface, and is not limited by reservoir shape or pressure distribution. It can calculate reservoirs of any shape and with any internal pressure distribution. Furthermore, the calculated deformation characteristic parameters can be used to predict surface subsidence, guiding oil and gas reservoir development strategies during underground fluid injection and production.

[0065] Accordingly, in one specific implementation, see [link to relevant documentation]. Figure 3 As shown, step S12 may include:

[0066] Step S121: Using a preset equivalent volume force formula, determine the equivalent volume force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process through the basic application parameters.

[0067] In this embodiment, the equivalent volume force when the porous medium undergoes elastic deformation is determined using a preset equivalent volume force formula.

[0068] The preset equivalent volume force formula is:

[0069] The parameters are pre-acquired basic application parameters, K. s K is the bulk modulus of the skeleton material. p ρ is the apparent bulk modulus of the rock, p is the pore fluid pressure, and ρ s ρ is the density of the rock skeleton material. l The average density of the fluid mixture in the pores is given, D is the depth data, g is the gravitational acceleration, and i represents the different axis directions of the three-dimensional Cartesian coordinate system.

[0070] Step S122: Using the preset equivalent area force formula, determine the equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process through the basic application parameters.

[0071] In this embodiment, the equivalent area force when the porous medium undergoes elastic deformation is determined using a preset equivalent area force formula.

[0072] The formula for the equivalent area force is:

[0073] Among them, K s K is the bulk modulus of the skeleton material. p Let Ω represent the apparent bulk modulus of the rock, and let Ω represent the interior of the closed boundary of the target region. n represents the outside of the closed boundary of the target region. Ω Let be the unit in-line normal vector inside the closed boundary at the closed boundary Γ.

[0074] Furthermore, in one feasible implementation, see [link to relevant documentation]. Figure 4 As shown, step S13 includes:

[0075] Step S131: By using the preset stress determination formula, the stress is determined by integrating the Mindlin fundamental solution, the equivalent volume force, and the equivalent area force according to the superposition principle.

[0076] In this embodiment, using the Mindlin fundamental solution for a semi-infinite body subjected to concentrated forces, a complete geomechanical solution, i.e., the stress, can be obtained by integrating according to the superposition principle using a preset stress determination formula. In engineering, we are more concerned with the changes in the problem, specifically the stress relative to the initial state. Therefore, the preset stress determination formula is:

[0077] in, The stress Green's function, determined based on the Mindlin fundamental solution, is used to determine the stress at point x = (x1, x2, x3) under stress x. k The unit force applied in the direction causes a Δσ at point ω = (ω1, ω2, ω3). ij Its analytical expression is:

[0078]

[0079]

[0080]

[0081]

[0082]

[0083]

[0084]

[0085]

[0086]

[0087]

[0088]

[0089]

[0090]

[0091]

[0092]

[0093]

[0094]

[0095]

[0096] In addition, among them Δp=p t＝T -p t ＝0 G is the shear modulus of the rock; v is the Poisson's ratio of the rock; x = (x1, x2, x3) is the point of application of the force; ω = (ω1, ω2, ω3) is the target calculation point;

[0097] Step S132: By using the preset strain determination formula, the strain is determined by integrating the Mindlin fundamental solution, the equivalent volume force, and the equivalent area force according to the superposition principle.

[0098] In this embodiment, using the Mindlin fundamental solution for a semi-infinite body subjected to concentrated forces, a complete geomechanical solution, i.e., strain, can be obtained by integrating according to the superposition principle using a preset strain determination formula. In engineering, we are more concerned with the changes in the problem, specifically the strain relative to the initial state. Therefore, the preset stress determination formula is:

[0099] in, The strain Green's function, determined based on the Mindlin fundamental solution, is used to determine the strain at point x = (x1, x2, x3) under stress x. k The unit force in the direction causes a Δε at point ω = (ω1, ω2, ω3). ij Its analytical expression is:

[0100] Where i = 1, 2, 3.

[0101] Furthermore, the volumetric strain produced by the rock can be calculated as follows:

[0102] in, The volumetric strain Green's function, determined based on the Mindlin fundamental solution, is used to determine the strain at point x = (x1, x2, x3) subjected to x. i The volumetric strain caused by a unit force in a given direction at point ω = (ω1, ω2, ω3) is expressed analytically as follows:

[0103]

[0104]

[0105]

[0106] In addition, among them Δp=p t＝T -p t ＝0 G is the shear modulus of the rock; v is the Poisson's ratio of the rock; x = (x1, x2, x3) is the point of application of the force; ω = (ω1, ω2, ω3) is the target calculation point;

[0107] Step S133: By using the preset displacement determination formula, the displacement is determined by integrating the Mindlin fundamental solution, the equivalent volume force, and the equivalent area force according to the superposition principle.

[0108] In this embodiment, using the Mindlin fundamental solution for a semi-infinite body subjected to concentrated forces, a complete geomechanical solution, i.e., the displacement, can be obtained by integrating according to the superposition principle using a preset displacement determination formula. In engineering, we are more concerned with the changes in the problem relative to the initial stress; that is, the preset displacement determination formula is:

[0109] in, The displacement Green's function, determined based on the Mindlin fundamental solution, is used to determine the stress x at the point x = (x1, x2, x3). iA unit force acting in a direction causes x at point ω = (ω1, ω2, ω3) to be... j The displacement in the direction; its analytical expression is:

[0110]

[0111]

[0112]

[0113]

[0114]

[0115]

[0116]

[0117] In addition, among them Δp=p t＝T -p t ＝0 G is the shear modulus of the rock; v is the Poisson's ratio of the rock; x = (x1, x2, x3) is the point of application of the force; ω = (ω1, ω2, ω3) is the target calculation point;

[0118] It should be pointed out that,

[0119] Where, Δp=p t＝T -p t ＝0 , During underground fluid injection and production, the fluid density varies significantly with pore pressure and fluid composition, while the rock skeleton density changes very little. Therefore, Δρ is neglected. s After Δf i e That is in, No fluid extraction occurs outside the reservoir's sealed boundary, and the change in pore pressure outside the sealed boundary is negligible. therefore It can be simplified to:

[0120] In this embodiment, the difficulty in studying actual reservoirs lies in their irregular shapes and uneven pressure variations. Since actual pore pressure variations can only be obtained through discrete numerical simulation, the equivalent volumetric force can be discretized as follows:

[0121]

[0122] Where a, b, and c are the grid cell numbers in dimensions x1, x2, and x3, respectively. The grid cell numbers increase in the same direction as the positive coordinate axes. Δx i It is the size of the i-dimensional grid cell, and the gradient value of the reservoir boundary is calculated using the one-sided difference method.

[0123] The equivalent area force can be discretized as follows:

[0124]

[0125] It should be noted that, in cases other than those described above, the equivalent area force is zero.

[0126] Correspondingly, the discrete stress is:

[0127]

[0128] Discrete strain is:

[0129]

[0130] Discrete volume strain is:

[0131]

[0132] The discrete displacement is:

[0133]

[0134] For example, consider a complex-shaped actual brine layer, such as an initial formation pressure of 17.2 MPa, a water layer thickness of 54 m, a water layer depth of 1472–1687 m, a shear modulus of G = 4.878 GPa, a Poisson's ratio of v = 0.23, Kp = 7.407 GPa, K0 = 400 GPa, an average porosity of 0.23, an average initial permeability of 15 mD, a simulation period of 15 years, and an average pore pressure rising to 29.8 MPa. Using the calculation method in this application embodiment, by inputting the porosity, saturation, phase fluid density, pore pressure increment, and three-dimensional coordinates of each grid, the equivalent volume force and equivalent area force of each grid can be calculated. Then, by selecting a target region, the deformation characteristic parameters of the target region can be calculated. Figure 5 The figure shows the calculated stress and strain in a target area using a specific geological profile as the model.

[0135] Accordingly, this application also discloses a geomechanical simulation device, see [link to relevant documentation]. Figure 6 As shown, the device includes:

[0136] Parameter acquisition module 11 is used to acquire basic application parameters of porous media during the underground fluid injection and extraction process in the target area;

[0137] The equivalent force determination module 12 is used to determine the equivalent volume force and equivalent area force when the porous medium undergoes elastic media deformation during the underground fluid injection and extraction process, based on the basic application parameters.

[0138] The deformation characteristic parameter determination module 13 is used to perform convolution between the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain the deformation characteristic parameters.

[0139] The strategy guidance module 14 is used to guide the oil and gas reservoir development strategy during the underground fluid injection and production process based on the deformation characteristic parameters.

[0140] For more detailed information on the working process of each of the above modules, please refer to the relevant content disclosed in the foregoing embodiments, which will not be repeated here.

[0141] Therefore, the above-described scheme in this embodiment includes: obtaining basic application parameters of porous media during the underground fluid injection and production process in the target area; determining the equivalent volume force and equivalent area force when the porous media undergoes elastic deformation during the underground fluid injection and production process using the basic application parameters; performing convolution on the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain deformation characteristic parameters; and guiding the oil and gas reservoir development strategy during the underground fluid injection and production process based on the deformation characteristic parameters. It is evident that this application equates the mechanical changes of porous media during the underground fluid injection and production process to a continuous medium problem subjected to certain volume forces and area forces. In other words, by using equivalent volume forces and equivalent area forces, the rock deformation during underground fluid injection and production is transformed into the deformation of a general elastic medium subjected to equivalent volume forces and equivalent area forces. A semi-analytical solution is then performed based on the Mindlin fundamental solution and convolution method, avoiding insufficient computational accuracy caused by excessively large far-field grids. This achieves accurate calculations across the entire spatial domain, making the calculation results unaffected by the range and fineness of the grid outside the reservoir. It exhibits significant advantages in far-field calculations, such as those at the surface, and is not limited by reservoir shape or pressure distribution. It can calculate reservoirs of any shape and with any internal pressure distribution. Furthermore, the calculated deformation characteristic parameters can be used to predict surface subsidence, guiding oil and gas reservoir development strategies during underground fluid injection and production.

[0142] Furthermore, embodiments of this application also disclose an electronic device, Figure 7 This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content of the diagram should not be construed as limiting the scope of this application.

[0143] Figure 7This is a schematic diagram of the structure of an electronic device 20 provided in an embodiment of this application. Specifically, 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 stores a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the geomechanical simulation method disclosed in any of the foregoing embodiments. Alternatively, the electronic device 20 in this embodiment may specifically be a computer.

[0144] 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 external devices, and the communication protocol it follows can be any communication protocol applicable to the technical solution of this application, and is not specifically limited here; the input / output interface 25 is used to acquire external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs, and is not specifically limited here.

[0145] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk, or optical disk, etc. The resources stored on it can include an operating system 221, computer programs 222, and data 223, etc. The data 223 can include various types of data. The storage method can be temporary storage or permanent storage.

[0146] The operating system 221 is used to manage and control the various hardware devices on the electronic device 20 and the computer program 222, which may be Windows Server, Netware, Unix, Linux, etc. In addition to including computer programs capable of performing the geomechanical simulation method executed by the electronic device 20 as disclosed in any of the foregoing embodiments, the computer program 222 may further include computer programs capable of performing other specific tasks.

[0147] Furthermore, this application also discloses a computer-readable storage medium, which includes random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disks, magnetic disks, optical disks, or any other form of storage medium known in the art. The computer program, when executed by a processor, implements the aforementioned geomechanical simulation method. Specific steps of this method can be found in the corresponding content disclosed in the foregoing embodiments, and will not be repeated here.

[0148] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since it corresponds to the method disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to in the method section.

[0149] The steps of the geomechanical simulation or algorithm described in conjunction with the embodiments disclosed herein can be implemented directly using hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.

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

[0151] The above provides a detailed description of the geomechanical simulation method, apparatus, equipment, and storage medium provided by the present invention. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A geomechanical simulation method, characterized in that, include: Obtain basic application parameters of porous media during the underground fluid injection and production process in the target area; The equivalent volume force and equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process are determined using the aforementioned basic application parameters. Based on the Mindlin fundamental solution, the equivalent volume force and the equivalent area force are convolved to obtain the deformation characteristic parameters; The deformation characteristic parameters are used to guide the oil and gas reservoir development strategy during the underground fluid injection and production process; The basic application parameters of porous media obtained during the underground fluid injection and production process in the target area include: The bulk modulus of the skeleton material, the apparent bulk modulus of the rock, the pore fluid pressure of the porous medium, the rock porosity, the density of the rock skeleton material, and the average density and depth data of the fluid mixture in the pores are obtained during the underground fluid injection and extraction process in the target area. The equivalent volume force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined using the aforementioned basic application parameters, including: Using a preset equivalent volume force formula, the equivalent volume force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined by the basic application parameters. The preset equivalent volume force formula is: ; in, Indicates the equivalent volume force in The directional components, where i represents different axes of the three-dimensional Cartesian coordinate system; The apparent bulk modulus of the rock is given by [the value of the value]. Where p is the bulk modulus of the skeleton material, ρ is the pore fluid pressure, and D is the depth data. where g is the acceleration due to gravity; The average density of the fluid mixture in the pores. The porosity of the rock is... The density of the rock skeleton material; The equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined using the aforementioned basic application parameters, including: Using a preset equivalent area force formula, the equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined by the basic application parameters. The formula for the equivalent area force is: ; in, For equivalent area force in directional components, This refers to the interior of the closed boundary of the target region. This indicates the outside of the closed boundary of the target area. Indicates closed boundary Pore ​​pressure on the outside, Indicates closed boundary The pore pressure on the inside, The interior of the closed boundary is within the closed boundary. The component of the unit in-place normal vector at a given location in the i-th direction.

2. The geomechanical simulation method according to claim 1, characterized in that, The method of convolving the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain deformation characteristic parameters includes: Based on the Mindlin fundamental solution, the equivalent volume force and the equivalent area force are convolved to obtain the stress, strain and displacement of the target region.

3. The geomechanical simulation method according to claim 2, characterized in that, The method of convolving the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain deformation characteristic parameters includes: The stress is determined by integrating the pre-set stress determination formula using the Mindlin fundamental solution, the equivalent volume force, and the equivalent area force, based on the superposition principle. The pre-set stress determination formula is as follows: ;in, For the point Stress at the point, , The stress Green's function, determined based on the Mindlin fundamental solution, is used to determine the stress Green's function in... Point received Point of action of unit force in the direction place ; The strain is determined by integrating the Mündlin fundamental solution, the equivalent volume force, and the equivalent area force using a preset strain determination formula, based on the superposition principle. The preset strain determination formula is as follows: ;in, For the point Strain at the point, The strain Green's function, determined based on the Mindlin fundamental solution, is used to determine the strain in... Point received Point of action of unit force in the direction place ; The displacement is determined by integrating the preset displacement determination formula using the Mindlin fundamental solution, the equivalent volume force, and the equivalent area force, based on the superposition principle. The preset displacement determination formula is as follows: ;in, For the point Displacement at that point The displacement Green's function, determined based on the Mindlin fundamental solution, is used to determine the displacement in... Point received Point of action of unit force in the direction Place Displacement in the direction; in, , , ; .

4. The geomechanical simulation method according to any one of claims 1 to 3, characterized in that, After determining the equivalent volume force and equivalent area force of the porous medium during elastic deformation in the underground fluid injection and extraction process using the basic application parameters, the method further includes: The equivalent volume force and the equivalent area force are discretized. Accordingly, the convolution solution based on the Mindlin fundamental solution to obtain the deformation characteristic parameters includes: Based on the fundamental solution of Mindlin, the equivalent volume force and the equivalent area force after discretization are convolved to obtain the deformation characteristic parameters after discretization.

5. A geomechanical simulation device, characterized in that, include: The parameter acquisition module is used to acquire the basic application parameters of porous media during the underground fluid injection and extraction process in the target area. The equivalent force determination module is used to determine the equivalent volume force and equivalent area force when the porous medium undergoes elastic media deformation during the underground fluid injection and extraction process, based on the basic application parameters. The deformation characteristic parameter determination module is used to perform convolution between the equivalent volume force and the equivalent area force based on the Mindlin fundamental solution to obtain the deformation characteristic parameters. The strategy guidance module is used to guide the oil and gas reservoir development strategy during the underground fluid injection and production process based on the deformation characteristic parameters. The parameter acquisition module is specifically used for: The bulk modulus of the skeleton material, the apparent bulk modulus of the rock, the pore fluid pressure of the porous medium, the rock porosity, the density of the rock skeleton material, and the average density and depth data of the fluid mixture in the pores are obtained during the underground fluid injection and extraction process in the target area. The equivalence determination module is specifically used for: Using a preset equivalent volume force formula, the equivalent volume force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined by the basic application parameters. The preset equivalent volume force formula is: ; in, Indicates the equivalent volume force in The directional components, where i represents different axes of the three-dimensional Cartesian coordinate system; The apparent bulk modulus of the rock is given by [the value of the value]. Where p is the bulk modulus of the skeleton material, ρ is the pore fluid pressure, and D is the depth data. where g is the acceleration due to gravity; The average density of the fluid mixture in the pores. The porosity of the rock is... The density of the rock skeleton material; Using a preset equivalent area force formula, the equivalent area force when the porous medium undergoes elastic deformation during the underground fluid injection and extraction process is determined by the basic application parameters. The formula for the equivalent area force is: ; in, For equivalent area force in directional components, This refers to the interior of the closed boundary of the target region. This indicates the outside of the closed boundary of the target area. Indicates closed boundary Pore ​​pressure on the outside, Indicates closed boundary The pore pressure on the inside, The interior of the closed boundary is within the closed boundary. The component of the unit in-place normal vector at a given location in the i-th direction.

6. An electronic device, characterized in that, The electronic device includes a processor and a memory; wherein the memory is used to store a computer program, which is loaded and executed by the processor to implement the geomechanical simulation method as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that, Used to store computer programs; wherein the computer programs, when executed by a processor, implement the geomechanical simulation method as described in any one of claims 1 to 4.

Images