Simulation method, device, equipment and medium for predicting flow and deformation of porous medium

Abstract

The application discloses a simulation method, device and equipment for predicting the flow and deformation of a porous medium and a medium, and the method comprises the following steps: selecting a calculation domain, generating a corresponding calculation grid, and initializing the porosity of the calculation domain; determining a time step and a discrete format, setting an uneven factor to solve a control equation set, wherein the control equation set comprises a continuity equation, a fluid dynamics equation, a skeleton mechanics equation and a porosity updating equation; and analyzing and displaying the solving result. The application provides a simulation method for predicting the flow and deformation of a porous medium, which is independent of assumptions and empirical parameters. The method is based on the basic principle of micro fluid-solid coupling action, adopts a porous medium fluid-solid coupling body dynamics equation, avoids the dependence of a traditional modeling method on volume averaging assumptions and the dependence on a multiphase action model, has cross-scale applicability, and provides a new theoretical basis and a prediction means for the technical development of a porous medium multiphase dynamics and a multi-field coupling problem.

Description

Technical Field

[0001] This application belongs to the field of numerical simulation technology of porous media, specifically relating to a simulation method, apparatus, equipment and medium for predicting the flow and deformation of porous media, which can be used to solve multiphase dynamics and multifield coupling problems of porous media. Background Technology

[0002] In both natural and engineering fields, numerous problems relate to fluid flow within deformable porous media. Accurate interphase coupling is crucial for predicting fluid flow and solid deformation, including applications such as subsurface reservoirs, hydraulic fracturing, slope instability, submarine landslides, and biomedical flows. Predicting the dynamics of fluids and solids accurately presents significant challenges, especially when dealing with complex and deformable geometries.

[0003] In specific fields such as geotechnical mechanics, multi-component media are often treated as continuous phases represented by statistical characteristics. In macroscopic constitutive models, the behavior of discrete particles, gases, and liquids is comprehensively considered. This treatment may limit the generality of the method and the accuracy of the corresponding predictions. This prompts us to develop an alternative method that should consider the behavior of each phase as much as possible based on fundamental principles.

[0004] In multiphase modeling, the Eulerian continuum model is the most widely used continuous phase model. However, because the fluid quantities defined in the computational domain lack applicability in the actual solid region, the model needs to be modified based on specific assumptions such as volume averaging and equilibrium assumptions. This often leads to different versions of the governing equations and the introduction of unclosed terms, relying on empirical models and further introducing uncertainty. Overall, the modeling is complex, and accuracy and universality are difficult to guarantee.

[0005] Therefore, there is an urgent need for a numerical simulation method for the flow and deformation of porous media that does not rely on assumptions and empirical parameters, so as to provide an accurate and reliable theoretical basis and prediction means for the development of technologies related to multiphase dynamics and multi-field coupling problems of porous media. Summary of the Invention

[0006] The purpose of this application is to provide a simulation method, apparatus, device and medium for predicting the flow and deformation of porous media, so as to facilitate the accurate solution of multiphase dynamics and multifield coupling problems of porous media, without relying on assumptions and empirical parameters.

[0007] The first aspect of this application provides a simulation method for predicting the flow and deformation of porous media, the method comprising:

[0008] Select the computational domain, generate the corresponding computational grid, and initialize the porosity of the computational domain;

[0009] A discretization method is determined and the governing equations are discretized. The time step and discretization scheme are determined, and a non-uniformity factor is set to solve the governing equations. The governing equations include a continuity equation, a fluid dynamics equation, a skeletal mechanics equation, and a porosity update equation; wherein, the fluid dynamics equation is:

[0010]

[0011]

[0012]

[0013] In the formula, α These are the coefficients of the transient terms in the fluid dynamics equations; β These are the convection coefficients in the fluid dynamics equations; γ These are the source term coefficients in the fluid dynamics equations; These are source terms in the fluid dynamics equations; The pore velocity at the next time step to be solved; t For time; For dynamic pressure; The viscosity of the fluid; g It is the acceleration due to gravity; For fluid density; Density of solid; φ Porosity of porous media; The non-uniformity factor is defined as the ratio of the resistance to fluid flow through an actual porous medium to the resistance to fluid flow through a uniform porous medium with the same porosity; σ s For solid stress; For time step; The velocity of the solid at the current time; The velocity of the fluid at the current time;

[0014] The mechanical equations of the skeleton are:

[0015]

[0016]

[0017]

[0018]

[0019] In the formula, λ represents the transient term coefficients of the skeleton mechanics equations; λ represents the source term coefficients of the skeleton mechanics equations. These are source terms in the solid dynamics equations; The velocity of the solid; The velocity of the fluid at the next time step;

[0020] The solution results are analyzed and presented.

[0021] In some embodiments, based on the principle of momentum conservation in the interaction between incompressible fluid and porous framework, fluid dynamics equations and framework mechanics equations describing the fluid-structure interaction in a single field are derived.

[0022] In some embodiments, the continuity equation is:

[0023]

[0024] The above is the continuity equation.

[0025] In some embodiments, the porosity update equation is:

[0026]

[0027] The above is the porosity update equation.

[0028] In some embodiments, the discretization method includes the finite volume method, the finite difference method, and the finite element method.

[0029] In some embodiments, the setting non-uniformity factor includes:

[0030] For simulation of homogeneous porous media, the inhomogeneity factor Take 1;

[0031] For simulation of non-uniform porous media, the inhomogeneity factor The change in velocity of the solid skeleton per unit time Δ u pore The velocity change Δ of the solid skeleton under the condition of uniform porous media with the same porosity u macro The ratio is obtained, that is =Δ u pore / Δ u macro Alternatively, it can be derived from existing penetration rate information.

[0032] In some embodiments, the flow parameters are first obtained by solving the continuity equation and the fluid dynamics equation, then the skeleton mechanics equation is solved using the forward Euler time-stepping scheme, and finally the porosity update equation is solved.

[0033] According to a second aspect of this application, a simulation apparatus for predicting the flow and deformation of porous media is provided. The apparatus applies the simulation method for predicting the flow and deformation of porous media as described in any one of the first aspects, and includes a preprocessing unit, a solution unit, and a result processing unit.

[0034] The preprocessing unit is used to select the computational domain, generate the corresponding computational grid, and initialize the porosity of the computational domain.

[0035] The solution unit is used to determine the discretization method and discretize the governing equations, determine the time step and discretization scheme, and set the inhomogeneity factor to solve the governing equations. The governing equations include continuity equations, fluid dynamics equations, skeletal mechanics equations, and porosity update equations; wherein the fluid dynamics equations are:

[0036]

[0037]

[0038]

[0039] In the formula, α These are the coefficients of the transient terms in the fluid dynamics equations; β These are the convection coefficients in the fluid dynamics equations; γ These are the source term coefficients in the fluid dynamics equations; These are source terms in the fluid dynamics equations; The pore velocity at the next time step to be solved; t For time; For dynamic pressure; The viscosity of the fluid; g It is the acceleration due to gravity; For fluid density; Density of solid; φ Porosity of porous media; The non-uniformity factor is defined as the ratio of the resistance to fluid flow through an actual porous medium to the resistance to fluid flow through a uniform porous medium with the same porosity; σ s For solid stress; For time step; The velocity of the solid at the current time; The velocity of the fluid at the current time;

[0040] The mechanical equations of the skeleton are:

[0041]

[0042]

[0043]

[0044]

[0045] In the formula, λ represents the transient term coefficients of the skeleton mechanics equations; λ represents the source term coefficients of the skeleton mechanics equations. These are source terms in the solid dynamics equations; The velocity of the solid; The velocity of the fluid at the next time step;

[0046] The results processing unit is used to analyze and display the solution results.

[0047] According to a third aspect of this application, a computer device is provided, comprising: a processor and a memory, the memory storing a program or instructions executable on the processor, wherein the program or instructions, when executed by the processor, implement the steps of the simulation method for predicting the flow and deformation of porous media as described in any one of the first aspects.

[0048] According to a fourth aspect of this application, a computer-readable storage medium is provided having a program or instructions stored thereon, which, when executed by a processor, implement the steps of the simulation method for predicting the flow and deformation of porous media as described in any one of the first aspects.

[0049] In summary, compared with the prior art, the above-described technical solutions conceived in this application can achieve the following beneficial effects:

[0050] This application provides a simulation method for the flow and deformation of porous media that does not rely on assumptions and empirical parameters. Based on the fundamental principle of microscopic fluid-structure interaction, this method employs the fluid-structure interaction dynamics equations for porous media, avoiding the reliance on volume averaging assumptions and multiphase interaction models found in traditional modeling methods. It also possesses cross-scale applicability. This method provides a new theoretical foundation and prediction tool for the technological development of multiphase dynamics and multi-field coupling problems in porous media. Attached Figure Description

[0051] Figure 1 A schematic flowchart illustrating a simulation method for predicting the flow and deformation of porous media provided in an embodiment of this application;

[0052] Figure 2 A solid distribution variation map predicted by the DBS method is provided in an embodiment of this application;

[0053] Figure 3 A solid distribution variation diagram predicted by the method of this application is provided for an embodiment of this application;

[0054] Figure 4 A block diagram of a simulation device for predicting the flow and deformation of porous media provided in an embodiment of this application;

[0055] Figure 5 This is a schematic diagram of the hardware structure of a computer device provided in an embodiment of this application. Detailed Implementation

[0056] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application. All other embodiments obtained by those skilled in the art based on the embodiments provided in this application without inventive effort are within the scope of protection of this application.

[0057] Obviously, the accompanying drawings described below are merely some examples or embodiments of this application. Those skilled in the art can apply this application to other similar scenarios based on these drawings without any inventive effort. Furthermore, it is understood that although the efforts made in this development process may be complex and lengthy, for those skilled in the art related to the content disclosed in this application, any changes to design, manufacturing, or production based on the technical content disclosed in this application are merely conventional technical means and should not be construed as insufficient disclosure of the content of this application.

[0058] In this application, the reference to "embodiment" means that a specific feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of this application. The appearance of this phrase in various places in the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment that is mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described in this application may be combined with other embodiments without conflict.

[0059] Unless otherwise defined, the technical or scientific terms used in this application shall have the ordinary meaning understood by one of ordinary skill in the art to which this application pertains. The terms “a,” “an,” “an,” “the,” and similar words used in this application do not indicate quantity limitation and may indicate singular or plural. The terms “comprising,” “including,” “having,” and any variations thereof used in this application are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or device that includes a series of steps or modules (units) is not limited to the listed steps or units, but may also include steps or units not listed, or may include other steps or units inherent to these processes, methods, products, or devices. The terms “connected,” “linked,” “coupled,” and similar words used in this application are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. “Multiple” used in this application refers to two or more. “And / or” describes the relationship between related objects, indicating that three relationships may exist; for example, “A and / or B” can represent: A alone, A and B simultaneously, and B alone. The character " / " generally indicates that the preceding and following objects are in an "or" relationship. The terms "first," "second," and "third" used in this application are merely to distinguish similar objects and do not represent a specific ordering of the objects.

[0060] This application provides a simulation method, apparatus, device, and medium for predicting the flow and deformation of porous media. The method is based on the fundamental principle of microscopic fluid-structure interaction and adopts the fluid-structure interaction dynamics equation of porous media. It avoids the reliance on volume averaging assumptions and multiphase interaction models of traditional modeling methods and has cross-scale applicability.

[0061] Figure 1 This is a schematic flowchart illustrating the simulation method for predicting the flow and deformation of porous media according to an embodiment of this application. Figure 1 As shown, the method includes the following steps:

[0062] S101: Preprocessing. Select the computational domain, generate the corresponding computational mesh, and initialize the porosity of the computational domain.

[0063] S102: Computation and Solution. Select a suitable discretization method and discretize the corresponding governing equations. Determine the solution time step and discretization scheme, and set the inhomogeneity factor to solve the governing equation set. This set of governing equations includes continuity equations, fluid dynamics equations, skeletal mechanics equations, and porosity update equations. The discretization scheme refers to using a discretization method to transform continuous mathematical equations into discretized algebraic equations for numerical solution by computer. Discretization methods include the finite volume method, finite difference method, and finite element method.

[0064] In some embodiments, the following continuity equation is employed:

[0065]

[0066] The above is the expression for the continuity equation.

[0067] Based on the principle of momentum conservation in the interaction between incompressible fluid and porous framework, the fluid dynamics equations and framework mechanics equations describing the fluid-structure interaction in a single field are derived. The fluid dynamics equations are as follows:

[0068]

[0069]

[0070]

[0071] In the formula, α These are the coefficients of the transient terms in the fluid dynamics equations; β These are the convection coefficients in the fluid dynamics equations; γ These are the source term coefficients in the fluid dynamics equations; These are source terms in the fluid dynamics equations; The pore velocity at the next time step to be solved; t For time; For dynamic pressure; The viscosity of the fluid; g It is the acceleration due to gravity; For fluid density; Density of solid; φ Porosity of porous media; The non-uniformity factor is defined as the ratio of the resistance to fluid flow through an actual porous medium to the resistance to fluid flow through a uniform porous medium with the same porosity; σ s For solid stress; For time step; The velocity of the solid at the current time; The velocity of the fluid at the current time.

[0072] The mechanical equations of the skeleton are:

[0073]

[0074]

[0075]

[0076]

[0077] In the formula, λ represents the transient term coefficients of the skeleton mechanics equations; λ represents the source term coefficients of the skeleton mechanics equations. These are source terms in the solid dynamics equations; The velocity of the solid; Let be the velocity of the fluid in the next time step.

[0078] The following porosity update equation is adopted:

[0079]

[0080] The above is the porosity update equation.

[0081] The solution sequence for the governing equations is as follows: first, solve for the flow parameters based on the continuity equation and the fluid dynamics equations; second, solve for the skeleton mechanics equations using the forward Euler time-step scheme; and finally, solve for the porosity update equations. The forward Euler time-step scheme is used when solving the skeleton mechanics equations.

[0082] The following discussion focuses on the unevenness factor. The settings are explained below:

[0083] For simulation of homogeneous porous media, the inhomogeneity factor Take 1;

[0084] For simulation of non-uniform porous media, the inhomogeneity factor The change in velocity of the solid skeleton per unit time Δ u pore The change in velocity Δ of the solid skeleton under the condition of uniform porous media with the same porosity u macro The ratio is obtained, that is =Δ u pore / Δ u macro Alternatively, it can be derived from existing penetration rate information.

[0085] S103: Post-processing. This involves visualizing and analyzing the calculation results, displaying images, and performing secondary analysis.

[0086] The derivation processes of the fluid dynamics equations and the skeleton mechanics equations are as follows:

[0087] (1) The fluid flow process is described by the following equation:

[0088]

[0089] In the formula, The force exerted by the porous media solid framework on the fluid is expressed as:

[0090]

[0091] right Perform the following transformation:

[0092]

[0093] Among them, the first term in the square brackets can be combined with the first term on the left side of the fluid flow process description equation, the second term is a known term, and the third term, the rate of change of the solid skeleton velocity, needs to be obtained by combining the skeleton mechanical equation.

[0094] (2) Similarly, the mechanical behavior of the skeleton is described by the following equations:

[0095]

[0096] In the formula, For the force exerted by the fluid on the porous medium solid framework, according to the basic principle of momentum conservation, the expression is:

[0097]

[0098] right Perform the following transformation:

[0099]

[0100] (3) The last term in the square brackets of the transformation formula is replaced with an equivalent term from the equation describing the mechanical behavior of the skeleton, and then... Expressions to represent the It is possible to obtain data that does not contain any unknown parameters related to solids. Expression. This Substituting these equations into the fluid flow process description equations, we obtain the final fluid dynamics equations:

[0101]

[0102]

[0103]

[0104] (4) After solving the fluid dynamics equations, all fluid-related parameters are known quantities. In the transformation formula It does not contain any unknown fluid-related parameters and can be directly combined with the equations describing the skeleton's mechanical behavior to obtain the corresponding skeleton mechanical equations:

[0105]

[0106]

[0107]

[0108]

[0109] This completes the derivation of the fluid dynamics equations and the skeletal mechanics equations.

[0110] The embodiments of this application provide a numerical simulation method for the flow and deformation of porous media that does not rely on assumptions and empirical parameters. This method provides a new theoretical basis and prediction means for the development of technologies related to multiphase dynamics and multi-field coupling problems of porous media.

[0111] The simulation method for predicting the flow and deformation of porous media presented in this application can be used to solve complex problems in many fields, such as the analysis of the coupled process of shale gas fracturing fluid-rock fracture propagation, the stability assessment of dam seepage-stress coupling, the early warning of landslides caused by rainfall-induced slope seepage-deformation, the prediction of coupled disasters of groundwater-surrounding rock in tunnel water and mud inrush, the coupled analysis of reservoir pressure drop-ground settlement in oil and gas extraction, the prediction of high-speed railway subgrade seepage-filler compression settlement, the safety assessment of groundwater-container corrosion coupling in nuclear waste disposal, and the stability analysis of foundation pit dewatering seepage-pit wall deformation.

[0112] To verify the ability of this application to predict flow and deformation, a simulation case of fluid flow causing deformation of a porous structure is considered. The system consists of a rectangular container with dimensions of 50 × 30 cm (discretized into 500 × 300 elements), with the top open to the atmosphere. The container is filled with a porous medium with a porosity of 0.4 (local porosity fluctuates randomly within ±0.05), a solid particle density of 2650 kg / m³, and a permeability of [missing information]. Water was injected at a rate of 6.5 mL / s through an injection well at the bottom boundary, and the porous medium was considered as a plastic solid with a yield stress of 0.22 m² / s². Under the combined action of viscous fluid flow, gravity, and buoyancy, the fluid surged, causing the plastic porous medium to move, ultimately leading to phenomena such as medium rupture, crack propagation, and surface bulging.

[0113] The above problem was solved using both the current method and the traditional Darcy-Brinkman–Stokes (DBS) method, and the simulation results of the two methods were compared. Figure 2 and Figure 3 This is a map showing the changes in solid distribution caused by water injection (crack initiation, cracking, and surface uplift), predicted by the DBS method and the current method, respectively. For example... Figure 2 and Figure 3As shown, although the governing equations of this application are completely different from the DBS equations, the simulation results of the two for the complex flow and deformation processes of porous media are highly consistent. Both methods successfully captured key processes such as crack propagation, hydraulic fracturing and surface uplift, indicating that the current method is accurate and feasible in dealing with complex fluid-structure interaction problems.

[0114] In terms of computational efficiency, the method of this invention demonstrates outstanding performance: thanks to the optimized design of interphase exchange terms, its total computation time is reduced by 17% compared to the DBS method, and the speed of each iteration is significantly improved. Especially after reaching steady-state flow, the computation time per time step is only 50% of that of the DBS method, which means that the current method will achieve a higher speedup for simulating long-term physical processes. These advantages further validate the practical value of the current method in handling flow and deformation problems in porous media.

[0115] It should be noted that the steps shown in the above process or in the flowchart of the accompanying figures can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.

[0116] This application also provides a simulation device for predicting the flow and deformation of porous media. These devices are used to implement the above embodiments and preferred embodiments, and will not be repeated as described previously. Although the devices described in the following embodiments are preferably implemented in software, hardware implementation or a combination of software and hardware is also possible and contemplated.

[0117] like Figure 4 As shown, the simulation device for predicting the flow and deformation of porous media in this application includes a preprocessing unit 201, a solution unit 202, and a result processing unit 203.

[0118] The preprocessing unit 201 is used to select the computational domain, generate the corresponding computational grid, and initialize the porosity of the computational domain.

[0119] Solver unit 202 is used to determine the discretization method and discretize the governing equations, determine the time step and discretization scheme, and set the inhomogeneity factor to solve the governing equations. The governing equations include a continuity equation, a fluid dynamics equation, a skeletal mechanics equation, and a porosity update equation; wherein, the fluid dynamics equation is:

[0120]

[0121]

[0122]

[0123] In the formula, α These are the coefficients of the transient terms in the fluid dynamics equations; β These are the convection coefficients in the fluid dynamics equations; γ These are the source term coefficients in the fluid dynamics equations; These are source terms in the fluid dynamics equations; The pore velocity at the next time step to be solved; t For time; For dynamic pressure; The viscosity of the fluid; g It is the acceleration due to gravity; For fluid density; Density of solid; φ Porosity of porous media; The non-uniformity factor is defined as the ratio of the resistance to fluid flow through an actual porous medium to the resistance to fluid flow through a uniform porous medium with the same porosity; σ s For solid stress; For time step; The velocity of the solid at the current time; The velocity of the fluid at the current time;

[0124] The mechanical equations of the skeleton are:

[0125]

[0126]

[0127]

[0128]

[0129] In the formula, λ represents the transient term coefficients of the skeleton mechanics equations; λ represents the source term coefficients of the skeleton mechanics equations. These are source terms in the solid dynamics equations; The velocity of the solid; The velocity of the fluid at the next time step;

[0130] The result processing unit 203 is used to analyze and display the solution results.

[0131] In addition, combined Figure 1 The simulation method for predicting the flow and deformation of porous media described in this application can be implemented by a computer device. Figure 5 This is a schematic diagram of the hardware structure of a computer device according to an embodiment of this application. Figure 5 As shown, the device may include a processor 301 and a memory 302 storing computer program instructions.

[0132] Specifically, the processor 301 may include a central processing unit (CPU), an application specific integrated circuit (ASIC), or one or more integrated circuits that can be configured to implement the embodiments of this application.

[0133] Memory 302 may include a large-capacity memory for data or instructions. For example, and not limitingly, memory 302 may include a hard disk drive (HDD), a floppy disk drive, a solid-state drive (SSD), flash memory, an optical disk drive, a magneto-optical disk drive, magnetic tape, or a Universal Serial Bus (USB) drive, or a combination of two or more of these. Where appropriate, memory 302 may include removable or non-removable (or fixed) media. Where appropriate, memory 302 may be internal or external to a data processing device. In a particular embodiment, memory 302 is non-volatile memory. In a particular embodiment, memory 302 includes read-only memory (ROM) and random access memory (RAM). Where appropriate, the ROM may be a mask-programmed ROM, a programmable read-only memory (PROM), an erasable read-only memory (EPROM), an electrically erasable read-only memory (EEPROM), an electrically alterable read-only memory (EAROM), or flash memory, or a combination of two or more of these. Where appropriate, the RAM can be Static Random-Access Memory (SRAM) or Dynamic Random-Access Memory (DRAM). DRAM can be Fast Page Mode Dynamic Random-Access Memory (FPMDRAM), Extended Data Out Dynamic Random-Access Memory (EDODRAM), Synchronous Dynamic Random-Access Memory (SDRAM), etc.

[0134] The memory 302 can be used to store or cache various data files that need to be processed and / or communicated, as well as possible computer program instructions executed by the processor 301.

[0135] The processor 301 reads and executes computer program instructions stored in the memory 302 to implement any of the simulation methods for predicting the flow and deformation of porous media in the above embodiments.

[0136] In some embodiments, the computer device may further include a communication interface 303 and a bus 300. For example, Figure 5 As shown, the processor 301, memory 302, and communication interface 303 are connected through bus 300 and complete communication with each other.

[0137] The communication interface 303 is used to enable communication between the various modules, devices, units, and / or equipment in the embodiments of this application. The communication interface 303 can also enable data communication with other components such as external devices, image / data acquisition devices, databases, external storage, and image / data processing workstations.

[0138] Bus 300 includes hardware, software, or both, that couples components of a computer device together. Bus 300 includes, but is not limited to, at least one of the following: data bus, address bus, control bus, expansion bus, and local bus. For example, and not as a limitation, bus 300 may include an Accelerated Graphics Port (AGP) or other graphics bus, an Extended Industry Standard Architecture (EISA) bus, a Front Side Bus (FSB), a Hyper Transport (HT) interconnect, an Industry Standard Architecture (ISA) bus, an InfiniBand interconnect, a Low Pin Count (LPC) bus, a memory bus, a MicroChannel Architecture (MCA) bus, a Peripheral Component Interconnect (PCI) bus, a PCI-Express (PCI-X) bus, a Serial Advanced Technology Attachment (SATA) bus, a Video Electronics Standards Association Local Bus (VLB) bus, or other suitable buses, or a combination of two or more of these. Where appropriate, bus 300 may include one or more buses. Although specific buses are described and illustrated in embodiments of this application, this application contemplates any suitable bus or interconnection.

[0139] The computer device can execute the simulation method for predicting the flow and deformation of porous media as described in the embodiments of this application, thereby achieving a combination of Figure 1 The simulation method described is used to predict the flow and deformation of porous media.

[0140] Furthermore, in conjunction with the simulation methods for predicting the flow and deformation of porous media in the above embodiments, this application embodiment can provide a computer-readable storage medium for implementation. This computer-readable storage medium stores computer program instructions; when executed by a processor, these computer program instructions implement any of the simulation methods for predicting the flow and deformation of porous media in the above embodiments.

[0141] It should be noted that the technical features of the above embodiments can be combined arbitrarily. For the sake of brevity, not all possible combinations of the technical features in the above embodiments have been described. However, as long as the combination of these technical features does not contradict each other, it should be considered within the scope of this specification. In addition, depending on the implementation needs, the various steps / components described in this application can be broken down into more steps / components, or two or more steps / components or parts of steps / components can be combined into new steps / components to achieve the purpose of this application.

[0142] It will be readily understood by those skilled in the art that the above-described embodiments merely illustrate several implementation methods of this application, and while the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this patent application should be determined by the appended claims.

Claims

1. A simulation method for predicting the flow and deformation of porous media, characterized in that, The method includes: Select the computational domain, generate the corresponding computational grid, and initialize the porosity of the computational domain; A discretization method is determined and the governing equations are discretized. The time step and discretization scheme are determined, and a non-uniformity factor is set to solve the governing equations. The governing equations include a continuity equation, a fluid dynamics equation, a skeletal mechanics equation, and a porosity update equation; wherein, the fluid dynamics equation is: In the formula, α These are the coefficients of the transient terms in the fluid dynamics equations; β These are the convection coefficients in the fluid dynamics equations; γ These are the source term coefficients in the fluid dynamics equations; These are source terms in the fluid dynamics equations; The pore velocity at the next time step to be solved; t For time; For dynamic pressure; The viscosity of the fluid; g It is the acceleration due to gravity; For fluid density; Φ represents the solid density; Φ represents the porosity of the porous medium. The non-uniformity factor is defined as the ratio of the resistance to fluid flow through an actual porous medium to the resistance to fluid flow through a uniform porous medium with the same porosity; σ s For solid stress; For time step; The velocity of the solid at the current time; The velocity of the fluid at the current time; The mechanical equations of the skeleton are: In the formula, λ represents the transient term coefficients of the skeleton mechanics equations; λ represents the source term coefficients of the skeleton mechanics equations. These are source terms in the solid dynamics equations; The velocity of the solid; The velocity of the fluid at the next time step; The solution results are analyzed and presented.

2. The simulation method for predicting the flow and deformation of porous media according to claim 1, characterized in that, Based on the principle of momentum conservation in the interaction between incompressible fluid and porous framework, we derive the fluid dynamics equations and framework mechanics equations that describe the fluid-structure interaction in a single field.

3. The simulation method for predicting the flow and deformation of porous media according to claim 1, characterized in that, The continuity equation is: The above is the continuity equation.

4. The simulation method for predicting the flow and deformation of porous media according to claim 1, characterized in that, The porosity update equation is as follows: The above is the porosity update equation.

5. The simulation method for predicting the flow and deformation of porous media according to claim 1, characterized in that, The discretization methods include the finite volume method, the finite difference method, and the finite element method.

6. The simulation method for predicting the flow and deformation of porous media according to claim 1, characterized in that, The non-uniformity factor includes: For simulation of homogeneous porous media, the inhomogeneity factor Take 1; For simulation of non-uniform porous media, the non-uniformity factor The change in velocity of the solid skeleton per unit time Δ u pore The velocity change Δ of the solid skeleton under the condition of uniform porous media with the same porosity u macro The ratio is obtained, that is =Δ u pore / Δ u macro Alternatively, it can be derived from existing penetration rate information.

7. The simulation method for predicting the flow and deformation of porous media according to claim 1, characterized in that, First, the continuity equation and fluid dynamics equation are solved to obtain the flow parameters. Then, the skeletal mechanics equation is solved using the forward Euler time-step scheme. Finally, the porosity update equation is solved.

8. A simulation device for predicting the flow and deformation of porous media, characterized in that, The device employs the simulation method for predicting the flow and deformation of porous media as described in any one of claims 1 to 7, and includes a preprocessing unit, a solution unit, and a result processing unit. The preprocessing unit is used to select the computational domain, generate the corresponding computational grid, and initialize the porosity of the computational domain. The solution unit is used to determine the discretization method and discretize the governing equations, determine the time step and discretization scheme, and set the inhomogeneity factor to solve the governing equations. The governing equations include continuity equations, fluid dynamics equations, skeletal mechanics equations, and porosity update equations; wherein the fluid dynamics equations are: In the formula, α These are the coefficients of the transient terms in the fluid dynamics equations; β These are the convection coefficients in the fluid dynamics equations; γ These are the source term coefficients in the fluid dynamics equations; These are source terms in the fluid dynamics equations; The pore velocity at the next time step to be solved; t For time; For dynamic pressure; The viscosity of the fluid; g It is the acceleration due to gravity; For fluid density; Φ represents the solid density; Φ represents the porosity of the porous medium. The non-uniformity factor is defined as the ratio of the resistance to fluid flow through an actual porous medium to the resistance to fluid flow through a uniform porous medium with the same porosity; σ s For solid stress; For time step; The velocity of the solid at the current time; The velocity of the fluid at the current time; The mechanical equations of the skeleton are: In the formula, λ represents the transient term coefficients of the skeleton mechanics equations; λ represents the source term coefficients of the skeleton mechanics equations. These are source terms in the solid dynamics equations; The velocity of the solid; The velocity of the fluid at the next time step; The results processing unit is used to analyze and display the solution results.

9. A computer device, characterized in that, include: A processor and a memory, wherein the memory stores a program or instructions that can run on the processor, and when the program or instructions are executed by the processor, implement the steps of the simulation method for predicting the flow and deformation of porous media as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, It stores a program or instructions that, when executed by a processor, implement the steps of the simulation method for predicting the flow and deformation of porous media as described in any one of claims 1 to 7.

Images