Four-dimensional ground stress coupling simulation method based on unified software framework

Abstract

The invention relates to the technical field of petroleum and natural gas engineering numerical simulation, and discloses a unified software framework-based four-dimensional ground stress coupling simulation method, which comprises the following steps of: S1, loading data to establish a four-dimensional ground stress coupling simulation project; s2, establishing a crustal stress model based on a finite element grid system; s3, establishing an oil reservoir numerical simulation model based on the angular point grid; and S4, controlling the ground stress coupling simulation process based on the fluid-structure interaction console. According to the invention, an integrated process from a static geologic model to numerical reservoir simulation and then to dynamic ground stress coupling analysis is realized, and an established unified coupling framework solves the problem that an existing numerical reservoir simulator and ground stress simulation software are mutually independent; the problems that fluid flow and crustal stress field dynamic coupling and cross-platform software data transmission are difficult to achieve are solved, and therefore simulation efficiency and result accuracy are improved. Technical support is provided for unconventional oil and gas reservoir well pattern optimization, refracturing design and development scheme adjustment.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of numerical simulation of oil and gas engineering, and particularly relates to a four-dimensional geostress coupling simulation method based on a unified software framework. BACKGROUND

[0002] The stress state of the formation has a significant control effect on the hydraulic fracture propagation in the process of unconventional oil and gas reservoir reconstruction, directly affects the rationality of well pattern deployment and fracturing parameter design, and increases the difficulty of reservoir plan adjustment. In addition, the change of the reservoir geostress will also cause the change of the pore structure of the rock, and then affect the flow characteristics of the fluid in the porous medium. At present, the established three-dimensional geological model of the reservoir, the rock mechanics model and the finite element model reflecting the geostress field based on the finite element method mainly reflect the static characteristics of the reservoir, and it is difficult to effectively depict the time and space dynamic changes of the geostress field caused by the engineering operation process such as oil and gas reservoir exploitation, infill drilling and repeated fracturing. Obviously, the traditional single physical field simulation method has been difficult to meet the description needs of the dynamic response characteristics of the reservoir geostress in the process of unconventional oil and gas development.

[0003] At present, the reservoir production dynamic simulation mainly uses reservoir numerical simulators such as Eclipse, Petrel RE, tNavigator and CMG, and the geostress field analysis usually relies on finite element mechanics simulation software such as Abaqus and Comsol. The separation type simulation method of the reservoir numerical simulation and the stress simulation has the following disadvantages: firstly, it is impossible to realize the real-time dynamic coupling between the fluid flow and the change of the geostress field; secondly, it is difficult to accurately represent the comprehensive influence of the stress field change on the physical properties of the reservoir rock and the fluid seepage behavior; thirdly, the data interfaces between different simulation software are not unified, and manual intervention is frequent, resulting in low simulation efficiency.

[0004] Therefore, it is urgent to develop a four-dimensional geostress coupling simulation method with a unified software framework, which can organically combine the reservoir numerical simulator and the finite element mechanics simulator. The method needs to establish a dynamic coupling mechanism between the fluid flow and the geostress field, and realize the automatic iterative solution of the two. Through such a coupling simulation method, the dynamic evolution law of the reservoir stress field in the process of oil and gas exploitation can be more accurately predicted, and the influence of the reservoir stress field on the production performance of the oil and gas well can be evaluated, thereby providing reliable technical support for the efficient development of unconventional oil and gas reservoirs. SUMMARY

[0005] The application aims to provide a four-dimensional geostress coupling simulation method based on a unified software framework, aiming at the independence of existing reservoir numerical simulators and finite element mechanics simulation software, the real-time dynamic coupling between fluid flow and geostress field change cannot be realized, the comprehensive influence of stress field change on reservoir rock physical properties and fluid seepage behavior cannot be accurately characterized, and the data interface between different simulation software is not unified, manual intervention is frequent, and the simulation efficiency is low, aiming to provide a unified software framework and data interface, complete the docking of different reservoir numerical simulators and geostress simulation software, and realize the establishment of a four-dimensional geostress fluid-solid coupling unified framework.

[0006] To achieve the above-mentioned purpose, the application adopts a four-dimensional geostress coupling simulation method based on a unified software framework, comprising the following steps: S1: According to the engineering data file, a four-dimensional geostress fluid-solid coupling simulation engineering is established; S2: According to the four-dimensional geostress fluid-solid coupling simulation engineering, a geostress model is established based on a finite element grid system; S3: According to the four-dimensional geostress fluid-solid coupling simulation engineering, an oil reservoir numerical simulation model is established based on an angle point grid; S4: According to the geostress model and the oil reservoir numerical simulation model, the geostress coupling simulation process is controlled based on a fluid-solid coupling control console.

[0007] Further, in the step S1, the following steps are specifically included: S11: Load the engineering data file: load the model basic data, including the geological model angle point grid system, the reservoir rock mechanics parameters, the reservoir pore depth data, the initial stress field data, the model boundary conditions, the model production system, and load the coupling configuration file; S12: Angle point grid system conversion: convert the geological model angle point grid system into a finite element grid system, and obtain the node and element correspondence relationship of the angle point grid system and the finite element grid system; S13: Calculate the interpolation operator between grids: based on the correspondence relationship of the angle point grid system and the finite element grid system, calculate the physical field interpolation operator between the two kinds of grid systems, including the node physical field interpolation operator and the element physical field interpolation operator; S14: Establish a unified coupling framework: including docking external geostress simulation software and numerical simulation software, establishing an input and output file directory, preparing initial data for the coupling process, allocating memory for the coupling process and performing memory initialization.

[0008] Further, in the step S2, the following steps are specifically included: S21: Geostress model grid data processing: according to the finite element grid system converted from the angle point grid, the geostress model grid data is established, and the grid data is stored in the format required by the geostress simulator; S22: Geostress model grid parameter interpolation calculation: interpolate the model rock mechanics parameters and initial physical field data into the geostress model finite element grid system according to the physical field interpolation operator; S23: Geostress model initial and boundary condition processing: convert the boundary data and initial data loaded in step S11 into the boundary conditions and initial conditions of the geostress model, and store them into the data format required by the geostress simulator; S24: Create geostress solver model data file: load the configuration file data in step S11, obtain the related setting parameters of the geostress simulator, and store the geostress model data into the task file format required by the geostress simulator.

[0009] Further, in the step S3, the following steps are specifically included: S31: Load the model data required by the reservoir numerical simulator: based on the corner point grid system and model data loaded in step S11, establish a reservoir numerical simulation model, including model grid data, initial conditions, boundary conditions, and reservoir numerical simulator configuration parameters; S32: Create a reservoir numerical simulator model data file: store all data of the reservoir numerical simulation model into a task file required by the reservoir numerical simulator.

[0010] Further, in the step S4, the following steps are specifically included: S41: Load the initial data of the time step t, including the stress field, the porosity and permeability field, the pressure field, and the saturation field; S42: Submit a task to the reservoir numerical simulator to simulate the reservoir dynamics at the current time step; S43: Obtain the simulation results of the reservoir simulator, and obtain the pore pressure field data after simulation; S44: Transfer the pore pressure field data to the geostress model through the fluid-solid coupling console; S45: Submit a task to the geostress simulator to simulate the geostress field under the current reservoir dynamics; S46: Obtain the simulation results of the geostress simulator, and update the reservoir porosity and permeability model according to the stress field data; the expression of the porosity and permeability model under the influence of geostress effect is:

[0011] wherein, , porosity of two time steps, , stress field of two time steps, , pressure field of two time steps, , , respectively represent the compressibility coefficients of rock (including fluid), fluid mixture phase and rock matrix, is the compressibility coefficient of rock (not including fluid), , respectively represent the current and initial pore permeability, , respectively represent the current and initial porosity, is the power index of the Kozney-Carman equation; S47: next time step t = t +1, when t < t n , update the initial data of the time step t , repeat steps S41-S46; when t >= t n , the fluid-solid iterative coupling is ended, the coupling simulation results are output, the memory is cleared, and the program exits.

[0012] The application further provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to realize the steps of the four-dimensional ground stress coupling simulation method based on the unified software framework.

[0013] The application further provides a computer device, which comprises a memory, a processor and a computer program stored in the memory and executable on the processor, and the processor realizes the steps of the four-dimensional ground stress coupling simulation method based on the unified software framework when executing the computer program.

[0014] The application further provides a computer program product, which comprises a computer program, and the computer program is executed by a processor to realize the steps of the four-dimensional ground stress coupling simulation method based on the unified software framework.

[0015] The four-dimensional ground stress coupling simulation method based on the unified software framework has the following beneficial effects: The application provides a four-dimensional ground stress coupling simulation method based on a unified software framework. BRIEF DESCRIPTION OF DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without any creative effort on the basis of these drawings.

[0017] Figure 1 Fig. 1 is a flow chart of the four-dimensional ground stress coupling simulation method based on the unified software framework of the present application; Figure 2 Fig. 2 is a schematic diagram of the step S1 processing process of the present application; Figure 3 Fig. 3 is a four-dimensional ground stress unified coupling framework provided by the present application; Figure 4 Fig. 4 is a schematic diagram of the step S2 processing process of the present application; Figure 5 Fig. 5 is a schematic diagram of the step S3 processing process of the present application; Figure 6 Fig. 6 is a schematic diagram of the step S4 processing process of the present application; Figure 7 Fig. 7 is a schematic diagram of the Abaqus software and tNavigator software fluid-solid coupling mode provided by the present application; Figure 8 Fig. 8 is a conversion result of the corner point grid and the finite element grid provided by the present application; Figure 9 Fig. 9 is a pore volume and conductivity distribution result at different time provided by the present application; Figure 10The application provides a coupling result and a comparison error distribution of a tNavigator commercial software result; Figure 11 The application provides a daily oil production and average pressure difference curve considering the effect of ground stress; Figure 12 The application provides a structural block diagram of the computer device. DETAILED DESCRIPTION

[0018] In order to facilitate understanding of the four-dimensional ground stress coupling simulation method based on the unified software framework provided by the application, the above technical solutions of the application are described in detail through specific embodiments.

[0019] Referring to Figure 1 , in view of the independence of the existing reservoir numerical simulator and the finite element mechanics simulation software, real-time dynamic coupling between fluid flow and ground stress field change cannot be realized, the comprehensive influence of stress field change on reservoir rock physical properties and fluid seepage behavior cannot be accurately represented, and the data interface between different simulation software is not unified, manual intervention is frequent, and the simulation efficiency is low, the application aims to provide a unified software framework and data interface, complete the connection of different reservoir numerical simulators and ground stress finite element software, and realize the establishment of a four-dimensional ground stress coupling simulation unified framework.

[0020] The method comprises the following steps: S1: establishing a four-dimensional ground stress fluid-solid coupling simulation project according to an engineering data file; S2: establishing a ground stress model based on a finite element grid system according to the four-dimensional ground stress fluid-solid coupling simulation project; S3: establishing a reservoir numerical simulation model based on a corner point grid according to the four-dimensional ground stress fluid-solid coupling simulation project; S4: controlling a ground stress coupling simulation process based on a fluid-solid coupling control console according to the ground stress model and the reservoir numerical simulation model.

[0021] Referring to Figure 2 , the step S1 specifically comprises the following steps: S11: loading an engineering data file: loading model basic data, including a geological model corner point grid system, reservoir rock mechanics parameters, reservoir pore depth data, initial stress field data, model boundary conditions, model production system, and loading a coupling configuration file; S12: corner point grid system conversion: converting the geological model corner point grid system into a finite element grid system, and obtaining the corresponding relationship between the corner point grid system and the finite element grid system nodes and units; S13: Calculate the interpolation operator between the grids: based on the correspondence between the corner point grid system and the finite element grid system, calculate the interpolation operator between the two grid systems, including the node physical field interpolation operator and the element physical field interpolation operator; S14: Establish a unified coupling framework: including interfacing with external stress simulation software and numerical simulation software, establishing input and output file directories, preparing initial data for the coupling process, allocating memory for the coupling process and initializing the memory.

[0022] Specifically, the unified coupling framework provides a unified interface through the coupling console to interface different stress simulators and numerical simulators, where the numerical simulator specifically refers to a simulation software for implementing specific business functions. In addition to implementing the functions in S11, S12, and S13, the main functions of the coupling console include memory space allocation, submission of computing tasks, task state monitoring, data transfer during the coupling process, input and output files, log recording, etc., as well as implementing a stress-sensitive model to transfer the stress effect to the numerical solver (see Figure 3 ).

[0023] Referring to Figure 4 , in step S2, specifically comprising the following steps: S21: Grid data processing of the stress model: according to the finite element grid system converted from the corner point grid, establish the grid data of the stress model, and store the grid data into the format required by the stress simulator; S22: Interpolation calculation of the stress model grid parameters: according to the physical field interpolation operator, interpolate the model rock mechanics parameters and the initial physical field data into the finite element grid system of the stress model; S23: Initial and boundary condition processing of the stress model: convert the boundary data and initial data loaded in step S11 into the boundary conditions and initial conditions of the stress model, and store them into the data format required by the stress simulator; S24: Create a stress solver model data file: load the configuration file data in step S11, obtain the related setting parameters of the stress simulator, and store the stress model data into the task file format required by the stress simulator.

[0024] Referring to Figure 5 , in step S3, specifically comprising the following steps: S31: Load the model data required by the reservoir numerical simulator: based on the corner point grid system and the model data loaded in step S11, establish a reservoir numerical simulation model, including model grid data, initial conditions, boundary conditions, reservoir numerical simulator configuration parameters, etc.; S32: Create a model data file for the reservoir numerical simulator: store all data of the reservoir numerical simulation model into the task file required by the numerical simulator.

[0025] Referring to Figure 6 , the step S4 specifically comprises the following steps: S41: loading initial data of time step t, including stress field, porosity and permeability field, pressure field, saturation field, etc. S42: submitting a task to the reservoir numerical simulator to simulate the reservoir performance of the current time step; S43: obtaining the simulation result of the reservoir simulator, and obtaining the data of the simulated pore pressure field; S44: passing the data of the pore pressure field to the geostress model through the fluid-solid coupling console; S45: submitting a task to the geostress simulator to simulate the geostress field under the current reservoir performance; S46: obtaining the simulation result of the geostress simulator, and updating the reservoir porosity and permeability model according to the stress field data; the expression of the porosity and permeability model under the influence of the geostress effect is:

[0026] wherein, , porosity of two time steps, , stress field of two time steps, , pressure field of two time steps, , , compressibility coefficient of rock (including fluid), fluid mixture phase and rock matrix, respectively, compressibility coefficient of rock (excluding fluid), , current and initial porosity and permeability, respectively, , current and initial porosity, respectively, power index of the Kozney-Carmen equation; S47: next time step t = t +1, when t < t n , the initial data of the time step t is updated, and the steps S41-S46 are repeated; when t >= t n , the fluid-solid iterative coupling is ended, the coupling simulation result is output, the memory is cleared, and the program is exited.

[0027] Referring to Figure 7 to Figure 11The application further provides a four-dimensional geostress coupling simulation method based on a unified software framework in a case application of tNavigator reservoir numerical simulator and Abaqus finite element simulation software. The tNavigator reservoir numerical simulator calculates the seepage response (seepage pressure and fluid mass) under a black oil model, and re-calculates and updates the reservoir mechanical response (reservoir stress and strain) according to obtained reservoir pore volume and permeability change data. The Abaqus finite element software solves the elastic deformation of a porous medium through an internally established pore pressure model to simulate the change of the reservoir geostress field. In the fluid-solid coupling simulation of tNavigator and Abaqus, tNavigator calculates the reservoir pore pressure field of each time step based on a flow model and transmits the result to Abaqus through a coupling console; Abaqus solves the corresponding reservoir stress field based on a finite element model. Both sides perform mesh conversion, model identification and data transmission through the coupling console to realize automatic interactive coupling iteration based on the pore permeability model.

[0028] Specifically, the fluid-solid coupling mode (see Figure 7 ) calculates from the initial geostress balance, and then in each time step, fluid flow simulation and rock stress simulation are alternately performed: first, the flow equation is solved under a fixed stress field to update the pore pressure, then the new pressure is transmitted to the mechanics model as a load to update the stress field, and finally the dynamically changing porosity is calculated according to the difference between the new stress and the pressure, and the iteration of the next time step is entered. The coupling console realizes the influence of the dynamically changing geostress field on the pore flow field by superimposing the geostress change effect on the keyword MULTPV in tNavigator.

[0029] Based on the tNavigator reservoir numerical simulation software and the Abaqus finite element software, the corner point grid model and the finite element tetrahedral grid model (see Figure 8 ) are obtained through the coupling console mesh conversion module, the model has one injection well and three production wells, the simulation time step is 144 minutes, and the total simulation time is 5 days. The pore volume and the conductivity distribution at different times in the coupling simulation result are shown in Figure 9 . By changing the time step (T=2, 20, 39), the relative errors of the reservoir pore volume, pressure and saturation are all less than 2%, which indicates that the numerical simulation result of the model has strong robustness and good convergence, and the framework is suitable for practical application of numerical simulation (see Figure 10 ). Figure 11For whether to consider the daily oil production and reservoir average pressure change curve brought by stress field change, it can be seen that the stress effect has obvious influence on oil production and average pressure field. At the same time, through comparison with the simulation result of internal coupling of tNavigator, it can be seen that the result obtained by the unified coupling framework established by the application is in good agreement with it, the difference between the two models in describing the reservoir pressure dynamic change and oil production is small, which shows that the established fluid-solid coupling simulation workflow is reliable and accurate.

[0030] The embodiment provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to realize the steps of the four-dimensional ground stress coupling simulation method based on the unified software framework. The storage medium can be a disc, an optical disc, a read-only memory (ROM), a random access memory (RAM), a flash memory, a hard disk drive (HDD) or a solid-state drive (SSD), etc. The storage medium can also include a combination of the above-mentioned types of memories.

[0031] The embodiment provides a computer device, which includes a memory, a processor and a computer program stored in the memory and capable of running on the processor, and the processor realizes the steps of the four-dimensional ground stress coupling simulation method based on the unified software framework when executing the program.

[0032] As Figure 12As shown, the computer device 120 can include at least one processor 121, such as a central processing unit (CPU), at least one communication interface 123, a memory 124, and at least one communication bus 122. The communication bus 122 is configured to enable communication between these components. The communication interface 123 can include a display, a keyboard, and optionally the communication interface 123 can further include a standard wired interface, a wireless interface. The memory 124 can be a high-speed random access memory (RAM), and can also be a non-volatile memory, such as at least one disk memory. The memory 124 can also be at least one storage device located away from the aforementioned processor 121. The memory 124 stores an application program, and the processor 121 invokes the program code stored in the memory 124 to execute any of the above method steps. The communication bus 122 can be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus, etc. The communication bus 122 can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 12The bus 123 is used to connect the above-mentioned elements in the system 120, and only one bus is represented, but it does not mean that there is only one bus or only one type of bus. Among them, the memory 124 can include volatile memory such as random-access memory (RAM); the memory can also include non-volatile memory such as flash memory, a hard disk drive (HDD) or a solid-state drive (SSD); the memory 124 can also include a combination of the above-mentioned types of memory. Among them, the processor 121 can be a central processing unit (CPU), a network processor (NP), or a combination of CPU and NP. The processor 121 can further include a hardware chip. The hardware chip can be an application-specific integrated circuit (ASIC), a programmable logic device (PLD), or a combination thereof. The PLD can be a complex programmable logic device (CPLD), a field-programmable gate array (FPGA), a generic array logic (GAL), or any combination thereof. Alternatively, the memory 124 is also used to store program instructions. The processor 121 can invoke the program instructions to implement the four-dimensional ground stress coupling simulation method based on the unified software framework as in the embodiment.

[0033] The embodiment provides a computer program product, including a computer program, which, when executed by a processor, implements the steps of the four-dimensional ground stress coupling simulation method based on the unified software framework.

[0034] The above only discloses a preferred embodiment of the present application, of course, cannot limit the scope of the present application, those skilled in the art can understand that the implementation of all or part of the above-mentioned processes, and according to the equivalent changes of the claims of the present application, still belong to the scope of the present application.

Claims

1. A four-dimensional geostress coupling simulation method based on a unified software framework, characterized in that, Comprise the following steps: S1: according to the engineering data file, establish four-dimensional geostress fluid-solid coupling simulation engineering; S2: according to four-dimensional geostress fluid-solid coupling simulation engineering, based on finite element grid system to establish geostress model; S3: according to four-dimensional geostress fluid-solid coupling simulation engineering, based on the corner point grid to establish reservoir numerical simulation model; S4: according to the geostress model and reservoir numerical simulation model, based on fluid-solid coupling control platform to control geostress coupling simulation process.

2. The four-dimensional geostress coupling simulation method based on a unified software framework according to claim 1, characterized in that, In step S1, as follows: S11: load engineering data file: load model basic data, including geologic model corner point grid system, reservoir rock mechanics parameters, reservoir pore depth data, initial stress field data, model boundary conditions, model production system, load coupling configuration file; S12: corner point grid system conversion: convert the geologic model corner point grid system into the finite element grid system, obtain the corresponding relationship between the corner point grid system and the finite element grid system nodes and elements; S13: calculate the interpolation operator between grids: based on the corresponding relationship between the corner point grid system and the finite element grid system, calculate the physical field interpolation operator between the two kinds of grid systems, including node physical field interpolation operator and element physical field interpolation operator; S14: establish a unified coupling framework: including docking external geostress simulation software and numerical simulation software, establishing input and output file directory, coupling process preparing initial data, allocating memory for coupling process and performing memory initialization.

3. The four-dimensional geostress coupling simulation method based on a unified software framework of claim 2, wherein, In step S2, as follows: S21: geostress model grid data processing: according to the finite element grid system converted from the corner point grid to establish the geostress model grid data, store the grid data into the format required by the geostress simulator; S22: geostress model grid parameter interpolation calculation: according to the physical field interpolation operator, interpolate the model rock mechanics parameters and initial physical field data into the finite element grid system of the geostress model; S23: geostress model initial and boundary condition processing: convert the boundary data and initial data loaded in step S11 into the boundary conditions and initial conditions of the geostress model, and store them into the data format required by the geostress simulator; S24: create geostress solver model data file: load the configuration file data in step S11, obtain the relevant setting parameters of the geostress simulator, store the geostress model data into the task file format required by the geostress simulator.

4. The four-dimensional geostress coupling simulation method based on a unified software framework of claim 3, wherein, In step S3, as follows: S31: load the model data required by the reservoir numerical simulator: based on the corner point grid system and model data loaded in step S11, establish the reservoir numerical simulation model, including model grid data, initial conditions, boundary conditions, and reservoir numerical simulator configuration parameters; S32: create reservoir numerical simulator model data file: store all data of the reservoir numerical simulation model into the task file required by the reservoir numerical simulator.

5. The four-dimensional geostress coupling simulation method based on a unified software framework of claim 1, wherein, In step S4, as follows: S41: load the initial data of time step t, including stress field, porosity and permeability field, pressure field, saturation field; S42: submit the task to the reservoir numerical simulator, simulate the reservoir dynamics of the current time step; S43: obtain the simulation results of the reservoir simulator, obtain the pore pressure field data after simulation; S44: passing the pore pressure field data to the geostress model through the fluid-solid coupling console; S45: submitting a task to the geostress simulator to simulate the geostress field under the current reservoir dynamic; S46: obtaining the simulation result of the geostress simulator, and updating the reservoir porosity-permeability model according to the stress field data; the porosity-permeability model under the influence of the geostress effect is expressed as: where, , denote the porosity at two time steps, , denote the stress field at two time steps, , denote the pressure field at two time steps, , , denote the compressibility coefficients of the rock (including the fluid), the fluid mixture phase and the rock matrix, respectively, denote the compressibility coefficient of the rock (not including the fluid), , denote the current and initial pore permeability, respectively, , denote the current and initial porosity, respectively, denote the power exponent of the Kozeny-Carman equation; S47: next time step t t t t n update the initial data of time step t , repeat steps S41-S46; when t t n the fluid-structure iterative coupling ends, output the coupling simulation results, clear the memory, and the program exits.​​​​ 6. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to implement the steps of the four-dimensional geostress coupling simulation method based on the unified software framework as claimed in any one of claims 1-5.

7. A computer device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor executes the computer program to implement the steps of the four-dimensional geostress coupling simulation method based on the unified software framework as claimed in any one of claims 1-5.

8. A computer program product comprising a computer program, characterized in that, The computer program is executed by the processor to implement the steps of the four-dimensional geostress coupling simulation method based on the unified software framework as claimed in any one of claims 1-5.

Images