Stratum deformation analysis method, device and system based on fluid-solid coupling calculation

Abstract

The invention relates to a stratum deformation analysis method, device and system based on fluid-solid coupling calculation. The method comprises the following steps: acquiring geological parameters of a target area, and constructing a geometric model of the target area based on the geological parameters; the geometric model comprises a stratum solid model and a fluid flow model; based on the stratum solid model and the fluid flow model, setting fluid-solid coupling boundary conditions at the interface of fluid and solid; based on a preset stratum stress field initial condition and a fluid-solid coupling boundary condition of the target area, solving the fluid flow control equation to obtain a first calculation result; taking the pressure distribution data as a boundary condition of a stratum solid model, solving a solid mechanical equilibrium equation, and calculating to obtain a second calculation result; target parameters in the fluid flow model are updated based on the deformation data and the stress distribution data, and the stratum solid model is updated based on the pressure distribution data. According to the scheme, the accuracy of stratum deformation analysis is improved.

Description

Technical Field

[0001] This disclosure relates to the field of formation deformation monitoring technology, and in particular to a formation deformation analysis method, device and system based on fluid-structure interaction calculation. Background Technology

[0002] In related technologies, formation fluid injection is a common engineering operation in fields such as oil and gas extraction and geological engineering. For example, it involves injecting water into the formation to maintain formation pressure and improve recovery rates. In the scenario of co-injection into two formations, changes in bottomhole injection pressure significantly affect formation deformation, which in turn affects the fluid flow characteristics. This fluid-solid interaction is complex. Traditional formation deformation analysis methods often consider fluids and solids separately, ignoring the coupling effect between them, leading to significant deviations between the analysis results and actual conditions. Although some fluid-solid coupling analysis methods have been applied to formation deformation research, existing methods are not accurate enough in simulating the dynamic interaction between bottomhole injection pressure and formation deformation when dealing with co-injection into two formations. They cannot comprehensively and accurately reflect the complexities of actual engineering situations and fail to meet the needs of engineering design and safety assessment. Summary of the Invention

[0003] To overcome the problems existing in related technologies, this disclosure provides a method, apparatus and system for analyzing formation deformation based on fluid-structure interaction calculations.

[0004] According to a first aspect of the present disclosure, a method for analyzing formation deformation based on fluid-structure interaction calculations is provided, comprising:

[0005] Geological parameters of the target area are obtained, and a geometric model of the target area is constructed based on the geological parameters; the geometric model includes a solid formation model and a fluid flow model.

[0006] Based on the formation solid model and fluid flow model, fluid-structure interaction boundary conditions are set at the fluid-solid interface; the fluid-structure interaction boundary conditions include stress equilibrium conditions and flow rate matching conditions.

[0007] Based on the preset initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target area, the fluid flow control equation is solved to obtain the first calculation result; the first calculation result includes the pressure distribution data and velocity distribution data of the fluid in the formation pores of the target area at the current moment;

[0008] Using the pressure distribution data as the boundary condition of the formation solid model, the solid mechanics equilibrium equations are solved to obtain a second calculation result; the second calculation result includes the deformation data and stress distribution data of the formation under fluid pressure.

[0009] The target parameters in the fluid flow model are updated based on the deformation data and stress distribution data. The formation solid model is updated based on the pressure distribution data. The process of reverting to the steps of solving the fluid flow control equation based on the preset initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target region, and obtaining the pressure distribution data and velocity distribution data of the fluid in the formation pores of the target region at the current moment, continues until both the first calculation result and the second calculation result satisfy the preset conditions. The target parameters include the pore structure information and permeability of the formation.

[0010] According to a second aspect of the present disclosure, a formation deformation analysis apparatus based on fluid-structure interaction calculation is provided, comprising:

[0011] A construction unit is used to acquire geological parameters of a target area and construct a geometric model of the target area based on the geological parameters; the geometric model includes a solid formation model and a fluid flow model;

[0012] The setting unit is used to set the fluid-solid coupling boundary conditions at the fluid-solid interface based on the formation solid model and the fluid flow model; the fluid-solid coupling boundary conditions include stress balance conditions and flow rate matching conditions.

[0013] The solving unit is used to solve the fluid flow control equations based on the preset initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target area, and obtain a first calculation result; the first calculation result includes the pressure distribution data and velocity distribution data of the fluid in the formation pores of the target area at the current moment;

[0014] The solving unit is also used to use the pressure distribution data as the boundary condition of the formation solid model, solve the solid mechanics equilibrium equation, and calculate the second calculation result; the second calculation result includes the deformation data and stress distribution data of the formation under fluid pressure;

[0015] The update unit is used to update the target parameters in the fluid flow model based on the deformation data and stress distribution data, update the formation solid model based on the pressure distribution data, and return to the steps of solving the fluid flow control equation based on the preset initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target region, obtaining the pressure distribution data and velocity distribution data of the fluid in the formation pores of the target region at the current time, until the first calculation result and the second calculation result both meet the preset conditions; the target parameters include the pore structure information and permeability of the formation.

[0016] According to a third aspect of the present disclosure, an electronic device includes: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the method as described in any one of the first aspects.

[0017] According to a fourth aspect of the present disclosure, a computer-readable storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the method as described in any one of the first aspects.

[0018] According to a fifth aspect of the present disclosure, a computer program product is provided, including a computer program that, when executed by a processor, implements the method as described in any one of the first aspects.

[0019] The technical solutions provided by the embodiments of this disclosure can include the following beneficial effects: Geological parameters of the target area are obtained, and a geometric model of the target area is constructed based on these parameters; the geometric model includes a formation solid model and a fluid flow model; based on the formation solid model and the fluid flow model, fluid-solid coupling boundary conditions are set at the fluid-solid interface; based on the preset initial conditions of the formation stress field and the fluid-solid coupling boundary conditions of the target area, the fluid flow control equation is solved to obtain a first calculation result; pressure distribution data is used as the boundary conditions of the formation solid model to solve the solid mechanics equilibrium equation, and a second calculation result is obtained; the target parameters in the fluid flow model are updated based on deformation data and stress distribution data, and the formation solid model is updated based on pressure distribution data. This takes into account the coupling effect between the fluid and solid, achieves accurate simulation of the target area, and improves the accuracy of formation deformation analysis.

[0020] It should be understood that the above general description and the following detailed description are exemplary and explanatory only, and are not intended to limit this disclosure. Attached Figure Description

[0021] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention.

[0022] Figure 1 This is a flowchart illustrating a formation deformation analysis method based on fluid-structure interaction calculations, according to an exemplary embodiment.

[0023] Figure 2 This is a block diagram illustrating a formation deformation analysis device based on fluid-structure interaction calculations, according to an exemplary embodiment.

[0024] Figure 3 This is a block diagram illustrating an apparatus for a formation deformation analysis method based on fluid-structure interaction calculations, according to an exemplary embodiment. Detailed Implementation

[0025] Exemplary embodiments will now be described in detail, examples of which are illustrated in the accompanying drawings. When the following description relates to the drawings, unless otherwise indicated, the same numerals in different drawings denote the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present invention. Rather, they are merely examples of apparatuses and methods consistent with some aspects of the invention as detailed in the appended claims.

[0026] The terminology used in this disclosure is for the purpose of describing particular embodiments only and is not intended to be limiting of the present disclosure. The singular forms “a” and “the” as used in this disclosure and the appended claims are also intended to include the plural forms, unless the context clearly indicates otherwise.

[0027] It should be understood that although the terms first, second, third, etc., may be used to describe various information in embodiments of this disclosure, such information should not be limited to these terms. These terms are only used to distinguish information of the same type from one another. For example, first information may also be referred to as second information without departing from the scope of embodiments of this disclosure, and similarly, second information may also be referred to as first information. Depending on the context, the words “if” and “suppose” as used herein may be interpreted as “when”, “when”, or “in response to a determination”.

[0028] Furthermore, various forms of processes shown in the embodiments of this disclosure can be used to reorder, add, or delete steps. For example, the steps described in this application can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution disclosed in this disclosure can be achieved, and no limitation is imposed herein.

[0029] Figure 1 This is a flowchart illustrating a formation deformation analysis method based on fluid-structure interaction calculations, according to an exemplary embodiment. Figure 1 As shown, it should be noted that the formation deformation analysis method based on fluid-structure interaction calculation in this embodiment of the present disclosure is applied to a formation deformation analysis device based on fluid-structure interaction calculation. Figure 1 As shown, the method may include the following steps:

[0030] Step 101: Obtain the geological parameters of the target area and construct a geometric model of the target area based on the geological parameters.

[0031] The geometric model includes a formation solid model and a fluid flow model.

[0032] In one embodiment, core samples can be obtained through geological drilling, and lithological composition can be analyzed using X-ray diffraction (XRD). Combined with well logging data (such as sonic logging and resistivity logging), formation thickness and porosity distribution can be determined. The collected discrete data is then imported into specialized data processing software (such as Petrel) for interpolation processing to generate a continuous formation solid model.

[0033] In addition, 3D modeling software can be used to construct a geometric model of two strata based on the processed stratigraphic data and the actual spatial morphology of the strata. During the modeling process, the boundary shapes and interlayer distribution of the strata are accurately depicted. Mechanical parameters such as elastic modulus and Poisson's ratio are assigned to the stratigraphic model based on rock mechanics experimental results; seepage parameters such as permeability are determined based on Darcy's flow experiments or empirical formulas and then input into the model. The finite element method is used, employing tetrahedral or hexahedral meshes to discretize the stratigraphic model. Mesh refinement is applied to key areas such as fluid-solid interfaces and regions of lithological variation to ensure the accuracy of model calculations.

[0034] In this embodiment, a fluid flow model can be established, treating the injected fluid as a compressible Newtonian fluid and considering parameters such as fluid density and viscosity. Based on Darcy's law, and combined with the formation's permeability and porosity, governing equations describing the fluid flow in the formation pores are established to calculate parameters such as fluid velocity and pressure distribution.

[0035] Furthermore, based on Darcy's law, and combined with the mass conservation equation, momentum conservation equation (in Darcy's law form), and energy conservation equation, we can establish the governing equations describing fluid flow in formation pores:

[0036]

[0037] Where ρ is the fluid density, φ is the porosity, and t is the time. This is the Darcy velocity vector.

[0038] The equation of motion is Where k is the permeability tensor, μ is the fluid viscosity, p is the pressure, g is the gravitational acceleration, and D is the elevation (depth, which can be positively downward).

[0039] Step 102: Based on the formation solid model and the fluid flow model, set the fluid-solid coupling boundary conditions at the fluid-solid interface.

[0040] Among them, fluid-structure interaction boundary conditions include stress equilibrium conditions and flow matching conditions.

[0041] In one embodiment, the stress equilibrium condition can be: at the fluid-solid interface, according to the principles of continuum mechanics, the fluid pressure p and the solid surface stress σ are balanced. nSatisfying p = σ n , where σ n This represents the normal stress on the solid surface. Through the coupling interface of the finite element software, fluid pressure is applied as a surface load to the corresponding boundaries of the solid model.

[0042] In some embodiments of this disclosure, the stress balance condition includes the fluid pressure being equal to the solid surface stress.

[0043] In some embodiments of this disclosure, traffic matching conditions include:

[0044] When a solid undergoes deformation, the rate of change of pore volume and the fluid flow rate satisfy the following formula:

[0045]

[0046] Where φ is porosity. Let Q be the pore volume change rate, and Q be the fluid flow rate.

[0047] In one embodiment, the flow matching condition can be: establishing a relationship between the fluid flow rate and the pore volume change rate caused by solid deformation. When the solid deforms, the porosity φ changes. This porosity change can be programmed to update the permeability and feed it back into the fluid flow model, ensuring that the fluid flow rate matches the solid deformation.

[0048] Step 103: Based on the preset initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target area, solve the fluid flow control equations to obtain the first calculation results.

[0049] The first calculation result includes the pressure distribution data and velocity distribution data of the fluid in the formation pores of the target area at the current moment.

[0050] In some embodiments of this disclosure, step 103 may specifically include:

[0051] Based on the initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target area, the fluid flow control equations are solved using the finite volume method or the finite element method to obtain the first calculation results.

[0052] In one embodiment, the fluid flow control equations can be solved using the finite volume method or the finite element method. Taking the finite element method as an example, the fluid control equations are discretized in space and transformed into a system of linear algebraic equations K. p =f, where K is the stiffness matrix, p is the pressure vector, and f is the load vector. The system of equations is solved using an iterative algorithm (such as the conjugate gradient method) to obtain the pressure and velocity distribution data of the fluid in the formation pores at the current moment.

[0053] Step 104: Using the pressure distribution data as the boundary condition of the formation solid model, solve the solid mechanics equilibrium equations to obtain the second calculation result.

[0054] The second calculation result includes deformation data and stress distribution data of the formation under fluid pressure.

[0055] In some embodiments of this disclosure, step 104 may specifically include:

[0056] Pressure distribution data is used as the boundary condition of the formation solid model. The solid mechanical equilibrium equation is discretized using the finite element method. The discretized solid mechanical equilibrium equation is solved using the Newton-Raphson iterative method, and the second calculation result is obtained.

[0057] In one embodiment, the pressure distribution obtained from fluid calculations can be applied as a surface load to the formation solid model, and the solid mechanics equilibrium equation ▽·σ+f can be solved in conjunction with the initial stress field. b =0, where σ is the stress tensor, f b The force is a volume force. The equations are discretized using the finite element method, and the nonlinear equations are solved using the Newton-Raphson iterative method to calculate the deformation and stress distribution data of the formation under fluid pressure.

[0058] Step 105: Update the target parameters in the fluid flow model based on deformation data and stress distribution data, update the formation solid model based on pressure distribution data, return to execute the initial conditions of the formation stress field and fluid-structure interaction boundary conditions based on the preset target area, solve the fluid flow control equation, and obtain the pressure distribution data and velocity distribution data of the fluid in the formation pores of the target area at the current moment, until the first calculation result and the second calculation result both meet the preset conditions.

[0059] The target parameters include the pore structure information and permeability of the formation.

[0060] In one embodiment, the aforementioned preset conditions, namely the calculation results of the fluid flow control equation and the solid mechanics equilibrium equation, are such that the changes in two consecutive iterations are both less than their respective preset convergence thresholds.

[0061] In one embodiment, the change in porosity φ can be calculated based on the formation displacement obtained from solid calculations, using an empirical formula. Where C p Δε is the porosity compressibility coefficient. v This represents volumetric strain. The permeability k is updated based on porosity changes using the formula... Where C is a constant. The updated permeability is fed back into the fluid flow model to recalculate the fluid flow state. Simultaneously, the new fluid pressure distribution is applied to the solid model again for the next round of solid calculations. A convergence threshold is set (e.g., pressure change less than 10). -3 MPa, displacement change less than 10 -5 m) When the changes in the calculation results of the fluid and solid in two consecutive iterations are less than this threshold, the fluid-structure interaction system is considered to have reached a stable state.

[0062] In some embodiments of this disclosure, multi-condition simulation and parameter analysis can be performed. That is, by changing key factors such as bottom-hole injection pressure, injected fluid flow rate, and formation physical parameters, simulation calculations can be conducted under various operating conditions. The patterns and characteristics of formation deformation under different operating conditions are analyzed, and the influence mechanism of bottom-hole injection pressure on formation deformation is studied. For example, the evolution of formation displacement and stress concentration areas with changes in injection pressure is analyzed. Through parameter sensitivity analysis, the key factors affecting formation deformation and their degree of influence are determined, providing data support and theoretical basis for engineering design.

[0063] Furthermore, multiple operating conditions can be designed, such as changing the injection pressure (e.g., increasing from 10MPa to 20MPa in 2MPa increments) and the injection fluid flow rate (e.g., from 50m³ / h to 20MPa). 3 / d at 10m 3 The gradient of / d increases to 100m 3 The simulation and data analysis are performed using a combination of parameters such as (d) and formation elastic modulus (fluctuating within 20% of the actual value range). Fluid-structure interaction calculations are conducted for each working condition to obtain data on formation deformation (such as displacement and strain), stress distribution, and fluid pressure distribution. Statistical methods and numerical analysis techniques are used to analyze the patterns and characteristics of formation deformation under different working conditions. The sensitivity coefficients of the parameters are calculated. Where y is the formation deformation index (such as maximum displacement) x i To determine the influencing parameters, we need to identify the key factors affecting formation deformation and their degree of influence.

[0064] In some embodiments of this disclosure, the calculated results of formation deformation, stress distribution, and fluid pressure distribution can be visualized, presenting the analysis results in intuitive graphical and chart formats to facilitate understanding and analysis by engineering technicians. Simultaneously, the calculation results are compared and verified with field monitoring data or physical experimental results to evaluate the accuracy and reliability of the analysis method. If deviations exist, the causes are analyzed, and the model and calculation parameters are optimized and adjusted to further improve the accuracy of the analysis method.

[0065] In this embodiment, by establishing a refined fluid-structure interaction model and employing an iterative calculation method, the dynamic interaction between bottom-hole injection pressure and formation deformation under a two-layer formation injection scenario can be accurately simulated. Compared with traditional methods, this significantly improves the accuracy and reliability of the analysis results, providing a more realistic reference for engineering design. Furthermore, the influence of multiple factors on formation deformation can be comprehensively considered. Through multi-condition simulation and parameter analysis, the laws and characteristics of formation deformation under different conditions can be studied in depth, providing engineers with comprehensive formation deformation analysis results. This helps to better understand the formation deformation mechanism and formulate reasonable engineering measures. Additionally, based on accurate analysis results, engineers can optimize formation fluid injection schemes, rationally control bottom-hole injection pressure and flow rate, avoid engineering safety problems caused by excessive formation deformation, such as formation fracturing and ground subsidence, while improving resource extraction efficiency and reducing engineering costs.

[0066] According to the formation deformation analysis method based on fluid-structure interaction calculation proposed in this disclosure, the geological parameters of the target area are obtained, and a geometric model of the target area is constructed based on the geological parameters. The geometric model includes a formation solid model and a fluid flow model. Based on the formation solid model and the fluid flow model, fluid-structure interaction boundary conditions are set at the fluid-solid interface. Based on the preset initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target area, the fluid flow control equation is solved to obtain a first calculation result. The pressure distribution data is used as the boundary condition of the formation solid model to solve the solid mechanics equilibrium equation to obtain a second calculation result. The target parameters in the fluid flow model are updated based on the deformation data and stress distribution data, and the formation solid model is updated based on the pressure distribution data. This method considers the coupling effect between fluid and solid, achieves accurate simulation of the target area, and improves the accuracy of formation deformation analysis.

[0067] Figure 2 This is a block diagram illustrating a formation deformation analysis device based on fluid-structure interaction calculations, according to an exemplary embodiment. (Refer to...) Figure 2 The device includes a construction unit 201, a setting unit 202, a solving unit 203, and an updating unit 204.

[0068] The construction unit 201 is used to acquire geological parameters of the target area and construct a geometric model of the target area based on the geological parameters; the geometric model includes a solid formation model and a fluid flow model.

[0069] The setting unit 202 is used to set the fluid-solid coupling boundary conditions at the fluid-solid interface based on the formation solid model and the fluid flow model; the fluid-solid coupling boundary conditions include stress balance conditions and flow matching conditions.

[0070] Solving unit 203 is used to solve the fluid flow control equation based on the preset initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target area, and obtain a first calculation result; the first calculation result includes the pressure distribution data and velocity distribution data of the fluid in the formation pores of the target area at the current moment;

[0071] The solver unit 203 is also used to use the pressure distribution data as the boundary condition of the formation solid model, solve the solid mechanics equilibrium equation, and calculate the second calculation result; the second calculation result includes the deformation data and stress distribution data of the formation under fluid pressure;

[0072] The update unit 204 is used to update the target parameters in the fluid flow model based on the deformation data and stress distribution data, update the formation solid model based on the pressure distribution data, and return to the steps of solving the fluid flow control equation based on the preset initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target region, and obtaining the pressure distribution data and velocity distribution data of the fluid in the formation pores of the target region at the current time, until the first calculation result and the second calculation result both meet the preset conditions; the target parameters include the pore structure information and permeability of the formation.

[0073] In some embodiments of this disclosure, the stress balance condition includes the fluid pressure being equal to the solid surface stress.

[0074] In some embodiments of this disclosure, the traffic matching conditions include:

[0075] When a solid undergoes deformation, the rate of change of pore volume and the fluid flow rate satisfy the following formula:

[0076]

[0077] Where φ is porosity. Let Q be the pore volume change rate, and Q be the fluid flow rate.

[0078] In some embodiments of this disclosure, the solving unit 203 may be specifically used to: solve the fluid flow control equations using the finite volume method or the finite element method based on the preset initial conditions of the formation stress field and the fluid-structure interaction boundary conditions of the target area, and obtain the first calculation result.

[0079] In some embodiments of this disclosure, the solving unit 203 may be specifically used to: use the pressure distribution data as the boundary condition of the formation solid model, discretize the solid mechanical equilibrium equation using the finite element method, solve the discretized solid mechanical equilibrium equation using the Newton-Raphson iteration method, and calculate the second calculation result.

[0080] Regarding the apparatus in the above embodiments, the specific manner in which each module performs its operation has been described in detail in the embodiments related to the method, and will not be elaborated upon here.

[0081] The formation deformation analysis device based on fluid-structure interaction calculation proposed in the embodiments of this disclosure is a method.

[0082] The technical solutions provided by the embodiments of this disclosure can include the following beneficial effects: Geological parameters of the target area are obtained, and a geometric model of the target area is constructed based on these parameters; the geometric model includes a formation solid model and a fluid flow model; based on the formation solid model and the fluid flow model, fluid-solid coupling boundary conditions are set at the fluid-solid interface; based on the preset initial conditions of the formation stress field and the fluid-solid coupling boundary conditions of the target area, the fluid flow control equation is solved to obtain a first calculation result; pressure distribution data is used as the boundary conditions of the formation solid model to solve the solid mechanics equilibrium equation, and a second calculation result is obtained; the target parameters in the fluid flow model are updated based on deformation data and stress distribution data, and the formation solid model is updated based on pressure distribution data. This takes into account the coupling effect between the fluid and solid, achieves accurate simulation of the target area, and improves the accuracy of formation deformation analysis.

[0083] Figure 3 This is a block diagram illustrating an apparatus for a formation deformation analysis method based on fluid-structure interaction calculations, according to an exemplary embodiment. For example, apparatus 300 may be an electronic device, such as a mobile phone, computer, digital broadcasting terminal, messaging device, game console, tablet device, medical device, fitness equipment, personal digital assistant, etc.

[0084] Reference Figure 3 The device 300 may include one or more of the following components: a processing component 302, a memory 304, a power component 306, a multimedia component 308, an audio component 310, an input / output (I / O) interface 312, a sensor component 314, and a communication component 316.

[0085] Processing component 302 typically controls the overall operation of device 300, such as operations associated with display, telephone calls, data communication, camera operation, and recording. Processing component 302 may include one or more processors 320 to execute instructions to perform all or part of the steps of the methods described above. Furthermore, processing component 302 may include one or more modules to facilitate interaction between processing component 302 and other components. For example, processing component 302 may include a multimedia module to facilitate interaction between multimedia component 308 and processing component 302.

[0086] Memory 304 is configured to store various types of data to support the operation of device 300. Examples of this data include instructions for any application or method operating on device 300, contact data, phonebook data, messages, pictures, videos, etc. Memory 304 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk.

[0087] The power supply component 306 provides power to the various components of the device 300. The power supply component 306 may include a power management system, one or more power sources, and other components associated with generating, managing, and distributing power to the device 300.

[0088] Multimedia component 308 includes a screen that provides an output interface between the device 300 and the user. In some embodiments, the screen may include a liquid crystal display (LCD) and a touch panel (TP). If the screen includes a touch panel, the screen may be implemented as a touchscreen to receive input signals from the user. The touch panel includes one or more touch sensors to sense touches, swipes, and gestures on the touch panel. The touch sensors may sense not only the boundaries of the touch or swipe action but also the duration and pressure associated with the touch or swipe operation. In some embodiments, multimedia component 308 includes a front-facing camera and / or a rear-facing camera. When the device 300 is in an operating mode, such as a shooting mode or a video mode, the front-facing camera and / or the rear-facing camera may receive external multimedia data. Each front-facing camera and rear-facing camera may be a fixed optical lens system or have focal length and optical zoom capabilities.

[0089] Audio component 310 is configured to output and / or input audio signals. For example, audio component 310 includes a microphone (MIC) configured to receive external audio signals when device 300 is in an operating mode, such as call mode, recording mode, and voice recognition mode. The received audio signals may be further stored in memory 304 or transmitted via communication component 316. In some embodiments, audio component 310 also includes a speaker for outputting audio signals.

[0090] I / O interface 312 provides an interface between processing component 302 and peripheral interface modules, such as keyboards, click wheels, buttons, etc. These buttons may include, but are not limited to, home buttons, volume buttons, power buttons, and lock buttons.

[0091] Sensor assembly 314 includes one or more sensors for providing status assessments of various aspects of device 300. For example, sensor assembly 314 may detect the on / off state of device 300, the relative positioning of components such as the display and keypad of device 300, changes in the position of device 300 or a component of device 300, the presence or absence of user contact with device 300, the orientation or acceleration / deceleration of device 300, and temperature changes of device 300. Sensor assembly 314 may include a proximity sensor configured to detect the presence of nearby objects without any physical contact. Sensor assembly 314 may also include a light sensor, such as a CMOS or CCD image sensor, for use in imaging applications. In some embodiments, sensor assembly 314 may also include an accelerometer, a gyroscope, a magnetometer, a pressure sensor, or a temperature sensor.

[0092] Communication component 316 is configured to facilitate wired or wireless communication between device 300 and other devices. Device 300 can access wireless networks based on communication standards, such as WiFi, 2G, or 3G, or combinations thereof. In one exemplary embodiment, communication component 316 receives broadcast signals or broadcast-related information from an external broadcast management system via a broadcast channel. In one exemplary embodiment, communication component 316 also includes a near-field communication (NFC) module to facilitate short-range communication. For example, the NFC module may be implemented based on radio frequency identification (RFID) technology, Infrared Data Association (IrDA) technology, ultra-wideband (UWB) technology, Bluetooth (BT) technology, and other technologies.

[0093] In an exemplary embodiment, the apparatus 300 may be implemented by one or more application-specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field-programmable gate arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic components to perform the methods described above.

[0094] In an exemplary embodiment, a non-transitory computer-readable storage medium including instructions is also provided, such as a memory 304 including instructions, which can be executed by a processor 320 of the device 300 to perform the above-described method. For example, the non-transitory computer-readable storage medium may be a ROM, random access memory (RAM), CD-ROM, magnetic tape, floppy disk, and optical data storage device, etc.

[0095] In an exemplary embodiment, a computer program product is also provided, including a computer program that implements the above-described method when executed by the processor 320 of the device 300.

[0096] Other embodiments of the invention will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein. This disclosure is intended to cover any variations, uses, or adaptations of the invention that follow the general principles of the invention and include common knowledge or customary techniques in the art not disclosed herein. The specification and examples are to be considered exemplary only, and the true scope and spirit of the invention are indicated by the following claims.

[0097] It should be understood that the present invention is not limited to the precise structure described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from its scope. The scope of the invention is limited only by the appended claims.

Claims

1. A method of formation deformation analysis based on fluid-structure coupling calculation, characterized by, The method comprises the following steps: obtaining geological parameters of a target area and constructing a geometric model of the target area based on the geological parameters; the geometric model comprises a formation solid model and a fluid flow model; setting a fluid-structure coupling boundary condition at a fluid-solid interface based on the formation solid model and the fluid flow model; the fluid-structure coupling boundary condition comprises a stress balance condition and a flow matching condition; solving a fluid flow control equation based on a preset initial condition of a formation stress field of the target area and the fluid-structure coupling boundary condition to obtain a first calculation result; the first calculation result comprises pressure distribution data and flow velocity distribution data of a fluid in a formation pore of the target area at a current time; using the pressure distribution data as a boundary condition of the formation solid model to solve a solid mechanics balance equation and obtain a second calculation result; the second calculation result comprises deformation data and stress distribution data of the formation under the action of fluid pressure; updating a target parameter in the fluid flow model based on the deformation data and the stress distribution data, updating the formation solid model based on the pressure distribution data, and returning to execute the step of solving the fluid flow control equation based on the preset initial condition of the formation stress field of the target area and the fluid-structure coupling boundary condition to obtain the pressure distribution data and the flow velocity distribution data of the fluid in the formation pore of the target area at the current time until the first calculation result and the second calculation result both satisfy a preset condition; the target parameter comprises pore structure information and permeability of the formation.

2. The fluid-solid coupling computation-based formation deformation analysis method of claim 1, wherein, the stress balance condition comprises that the fluid pressure is equal to the surface stress of the solid.

3. The fluid-structure interaction computation-based formation deformation analysis method of claim 1, wherein, the flow matching condition comprises: in the case that the solid deforms, the pore volume change rate and the fluid flow satisfy the following formula: where φ is the porosity, is the rate of change of pore volume, Q is the fluid flow rate.

4. The fluid-solid coupling computation-based formation deformation analysis method of claim 1, wherein, the step of solving the fluid flow control equation based on the preset initial condition of the formation stress field of the target area and the fluid-structure coupling boundary condition to obtain the first calculation result comprises: solving the fluid flow control equation based on the preset initial condition of the formation stress field of the target area and the fluid-structure coupling boundary condition by using a finite volume method or a finite element method to obtain the first calculation result.

5. The fluid-structure interaction computation-based formation deformation analysis method of claim 1, wherein, the step of using the pressure distribution data as the boundary condition of the formation solid model to solve the solid mechanics balance equation and obtain the second calculation result comprises: using the pressure distribution data as the boundary condition of the formation solid model, discretely processing the solid mechanics balance equation by using a finite element method, solving the discretely processed solid mechanics balance equation by using a Newton-Raphson iteration method, and obtaining the second calculation result.

6. A formation deformation analysis device based on fluid-structure coupling calculation, characterized by, The method comprises the following steps: a construction unit is configured to obtain geological parameters of a target area and construct a geometric model of the target area based on the geological parameters; the geometric model comprises a formation solid model and a fluid flow model; a setting unit is configured to set a fluid-structure coupling boundary condition at a fluid-solid interface based on the formation solid model and the fluid flow model; the fluid-structure coupling boundary condition comprises a stress balance condition and a flow matching condition; The solving unit is configured to solve a fluid flow control equation based on preset initial conditions of a stress field of the target area and the fluid-solid coupling boundary condition to obtain a first calculation result, wherein the first calculation result comprises pressure distribution data and flow velocity distribution data of the fluid in the pore of the target area at a current time point; The solving unit is further configured to solve a solid mechanics equilibrium equation based on the pressure distribution data as a boundary condition of the formation solid model to obtain a second calculation result, wherein the second calculation result comprises deformation data and stress distribution data of the formation under the action of the fluid pressure; The updating unit is configured to update a target parameter in the fluid flow model based on the deformation data and the stress distribution data, update the formation solid model based on the pressure distribution data, and return to execute the step of solving the fluid flow control equation based on the preset initial conditions of the stress field of the target area and the fluid-solid coupling boundary condition to obtain the pressure distribution data and the flow velocity distribution data of the fluid in the pore of the target area at the current time point until the first calculation result and the second calculation result both satisfy preset conditions. The target parameter comprises pore structure information and permeability of the formation.

7. An electronic device, comprising: The method comprises: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method according to any one of claims 1 to 5.

8. 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 method according to any one of claims 1 to 5.

9. A computer program product comprising a computer program, characterized in that, The computer program is executed by the processor to implement the method according to any one of claims 1 to 5.

Images