CN116542037A - Numerical simulation method and device for fluid channeling law of first interface of well cementing cement sheath - Google Patents
Numerical simulation method and device for fluid channeling law of first interface of well cementing cement sheath Download PDFInfo
- Publication number
- CN116542037A CN116542037A CN202310465253.5A CN202310465253A CN116542037A CN 116542037 A CN116542037 A CN 116542037A CN 202310465253 A CN202310465253 A CN 202310465253A CN 116542037 A CN116542037 A CN 116542037A
- Authority
- CN
- China
- Prior art keywords
- model
- cement sheath
- interface
- well cementation
- crossflow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000004568 cement Substances 0.000 title claims abstract description 120
- 239000012530 fluid Substances 0.000 title claims abstract description 87
- 230000005465 channeling Effects 0.000 title claims abstract description 84
- 238000004088 simulation Methods 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000008859 change Effects 0.000 claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims description 24
- 239000000463 material Substances 0.000 claims description 16
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000002347 injection Methods 0.000 claims description 7
- 239000007924 injection Substances 0.000 claims description 7
- 238000011065 in-situ storage Methods 0.000 claims description 5
- 239000011396 hydraulic cement Substances 0.000 claims description 4
- 238000010206 sensitivity analysis Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 2
- 239000007789 gas Substances 0.000 description 20
- 238000010586 diagram Methods 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 239000011148 porous material Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000007789 sealing Methods 0.000 description 2
- 241000218631 Coniferophyta Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000002343 natural gas well Substances 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002341 toxic gas Substances 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Geometry (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Computing Systems (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a numerical simulation method and a device for a cement sheath first interface channeling law, wherein the method comprises the following steps: acquiring well cementation data of a well cementation cement sheath to be tested; establishing a cross flow model of a first interface of the well cementing cement sheath according to the well cementing data, wherein the cross flow model comprises a geometric model and a grid model; determining a crossflow rule theoretical formula of a crossflow model according to the geometric model; dividing the grid number of the grid model and setting the cross flow condition of the cross flow model; and obtaining a numerical simulation result of the height change rule of the first interface crossflow according to the crossflow model. According to the invention, the fluid channeling model of the first interface of the well cementation cement sheath is constructed, the structures of the well cementation cement sheath are simulated by the geometric model, the stress state of the structures of the well cementation cement sheath is modeled by the grid model, and the fluid channeling rule of the first interface of the well cementation cement sheath under the conditions of high temperature and high pressure is simulated and analyzed with high precision.
Description
Technical Field
The invention relates to the technical field of well cementation, in particular to a numerical simulation method and device for a first interface channeling law of a well cementation cement sheath and electronic equipment.
Background
Natural gas is a clean and environment-friendly excellent energy source, and along with the collection and development of natural gas, weather heat exploitation extends towards unconventional high-temperature high-pressure complex geological conditions. However, in the well cementation process of the natural gas well under complex geological conditions, the cementing quality of the sleeve, the cement sheath and the stratum is difficult to ensure, the well cementation tightness is easily damaged due to the influence of the later operation and high-temperature high-pressure conditions, the sealing failure of the cement sheath is a current ubiquitous problem, and the phenomenon of ring blank belt pressure is caused after the sealing failure of the interface of the cement sheath.
If the first interface of the cement ring (the cementing interface of the pipe sleeve and the cement ring) is damaged, the gas of the stratum is easy to flow into the first interface of the cement ring due to the action of high temperature and high pressure, and can be quickly diffused upwards at the interface in a short time to generate a cross flow phenomenon, so that the safety problems of unstable reservoir pressure, stranding among reservoirs with different pressure, toxic gas outflow and the like can be caused. In order to ensure the safety of a gas well in exploitation operation, analysis and research on a gas channeling rule when a first interface is damaged are needed.
The conventional analysis of the well cementation cement sheath fluid of the high-temperature high-pressure gas well is divided into a theoretical analysis method and an indoor simulation experiment method, but both methods have certain limitations, the theoretical analysis method is very difficult to analyze the first interface fluid of the well cementation cement sheath under the high-temperature high-pressure condition, and the first interface fluid of the well cementation cement sheath is difficult to accurately analyze; the indoor simulation experiment method has higher cost and safety risk, and the accuracy of the obtained result is lower because the indoor simulation experiment method is difficult to duplicate complex geological conditions such as high temperature and high pressure.
Disclosure of Invention
In view of the foregoing, it is necessary to provide a numerical simulation method, a device and an electronic device for a cement sheath first interface channeling law of well cementation, which are used for solving the technical problem in the prior art that complex conditions are difficult to be re-carved, so that analysis accuracy of the cement sheath first interface channeling law of well cementation under high temperature and high pressure conditions is low.
In order to solve the problems, the invention provides a numerical simulation method for a cement sheath first interface channeling law, which comprises the following steps:
acquiring well cementation data of a well cementation cement sheath to be tested;
establishing a cross flow model of a first interface of the well cementing cement sheath according to the well cementing data, wherein the cross flow model comprises a geometric model and a grid model;
determining a crossflow rule theoretical formula of a crossflow model according to the geometric model;
dividing the grid number of the grid model and setting the cross flow condition of the cross flow model;
and obtaining a numerical simulation result of the height change rule of the first interface crossflow according to the crossflow model.
Further, obtaining well cementation data of a well cementation cement sheath to be tested includes:
acquiring material parameters and well cementation interface parameters of a well cementation cement sheath to be tested, wherein the material parameters comprise the geometric dimension, young modulus and Poisson's ratio of a well cementation material; the well cementation interface parameters comprise shear strength, hydraulic cement strength, rigidity coefficient and critical fracture energy.
Further, establishing a fluid channeling model of a first interface of the well cementing cement sheath according to the well cementing data, including:
building a geometric model of a well cementation cement sheath according to the well cementation data;
dividing the geometric model into a plurality of basic stress analysis units by a plurality of hexahedral meshes, and constructing a mesh model of the well cementation cement sheath.
Further, determining a crossflow law theoretical formula of a crossflow model according to the geometric model, including:
analyzing the state of the fluid channeling of the cement sheath to be tested based on the Navigator equation and the laminar flow assumption, and determining the flow law of the fluid channeling in the cement sheath;
determining a flow velocity distribution type in the annular gap based on a slip-free boundary condition of the annular gap and according to a flow law of a cross flow in the well cementation cement ring;
the flow velocity in the annular gap is distributed and integrated to obtain a flow expression of the annular gap, and a flow approximation expression in the annular gap is obtained according to the law of three times;
obtaining a flow rule expression of the microcosmic-layer gas in the annular space based on the Fuchshimer equation;
and obtaining a crossflow rule theoretical formula of the crossflow model according to the flow velocity distribution type, the flow similarity expression and the flow rule expression in the annular gap.
Further, the mesh number division of the mesh model includes:
setting an error standard threshold, carrying out grid sensitivity analysis on the grid model, determining the minimum grid number of which the error does not exceed the error standard threshold in the process of calculating the cross flow height rule numerical simulation, and carrying out grid number division according to the minimum grid number and a preset division principle.
Further, setting a crossflow condition of the crossflow model includes:
setting initial conditions of a model, wherein the initial conditions comprise applying an in-situ stress field to a model body of the channeling model, and defining an initial damage unit and an injection point boundary;
setting boundary conditions of the model, wherein the boundary conditions comprise displacement constraints of the geometric model in the three-dimensional direction;
based on the ground stress balance analysis step, eliminating the initial displacement of the model body in the channeling model, and applying the initial stress of the cement sheath interface;
and determining the ground stress according to the initial stress and the crossflow model, and applying the ground stress to the crossflow model.
Further, obtaining a numerical simulation result of a height change rule of the crossflow according to the crossflow model, including:
determining an influence parameter influencing the height of the fluid channeling of a first interface of the well cementing ring according to a fluid channeling rule theoretical formula of the fluid channeling model;
and obtaining a numerical simulation result of the influence rule of the influence parameter on the height of the crossflow according to the crossflow model.
Further, obtaining a numerical simulation result of the influence rule of the influence parameter on the height of the crossflow according to the crossflow model, including:
and taking any parameter in the influence parameters as a single variable to carry out simulation input to obtain a corresponding crossflow height, and obtaining a numerical simulation result of the influence rule of the influence parameters on the crossflow height according to the crossflow height.
The invention also provides a numerical simulation device for the cement sheath first interface channeling law, which comprises:
the data acquisition unit is used for acquiring well cementation data of the well cementation cement sheath to be tested;
the model construction unit is used for building a cross flow model of a first interface of the well cementing cement sheath according to the well cementing data, and the cross flow model comprises a geometric model and a grid model;
the model analysis unit is used for determining a crossflow rule theoretical formula of the crossflow model according to the geometric model;
the condition setting unit is used for dividing the grid number of the grid model and setting the cross flow condition of the cross flow model;
and the simulation result output unit is used for obtaining a numerical simulation result of the height change rule of the first interface crossflow according to the crossflow model.
The invention also provides an electronic device comprising a memory and a processor, wherein,
the memory is used for storing programs;
the processor is coupled with the memory and is used for executing the program stored in the memory so as to realize the steps in the numerical simulation method for the cement sheath first interface channeling law of any one of the above.
Compared with the prior art, the beneficial effects of adopting the embodiment are as follows: in the numerical simulation method of the fluid channeling law of the first interface of the well cementation cement sheath, firstly, well cementation data of the well cementation cement sheath to be tested are obtained; then establishing a cross flow model of a first interface of the well cementing cement sheath according to the well cementing data, wherein the cross flow model comprises a geometric model and a grid model; determining a crossflow rule theoretical formula of a crossflow model according to the geometric model; dividing the grid number of the grid model and setting the cross flow condition of the cross flow model; and finally, obtaining a numerical simulation result of the height change rule of the first interface crossflow according to the crossflow model. In summary, the invention realizes the technical effect of high-precision simulation analysis on the fluid channeling law of the first interface of the well cementation cement sheath under the conditions of high temperature and high pressure by constructing the fluid channeling model of the first interface of the well cementation cement sheath and simulating each structure of the well cementation cement sheath by using the geometric model and modeling the stress state of each structure of the well cementation cement sheath by using the grid model.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the description of the embodiments will be briefly described below, it being evident that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of an embodiment of a method for simulating a fluid channeling law value of a first interface of a well cementation cement sheath;
FIG. 2 is a schematic diagram of a geometric model and a corresponding mesh model of a sleeve portion according to an embodiment of the present invention;
FIG. 3 is a schematic structural diagram of an embodiment of a numerical simulation device for a fluid channeling law of a first interface of a well cementation cement sheath;
fig. 4 is a schematic structural diagram of an electronic device provided by the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. It will be apparent that the described embodiments are only some, but not all, embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
It should be understood that the drawings of the schematic drawings are not drawn to scale. A flowchart, as used in this disclosure, illustrates operations implemented according to some embodiments of the present invention. It should be appreciated that the operations of the flow diagrams may be implemented out of order and that steps without logical context may be performed in reverse order or concurrently. Moreover, one or more other operations may be added to or removed from the flow diagrams by those skilled in the art under the direction of the present disclosure.
Some of the block diagrams shown in the figures are functional entities and do not necessarily correspond to physically or logically separate entities. These functional entities may be implemented in software or in one or more hardware modules or integrated circuits or in different networks and/or processor systems and/or microcontroller systems.
Reference herein to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment may be included in at least one embodiment of the invention. The appearances of such phrases in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. Those of skill in the art will explicitly and implicitly appreciate that the embodiments described herein may be combined with other embodiments.
Fig. 1 is a flow chart of an embodiment of a method for simulating a fluid channeling law of a first interface of a well cementation cement sheath, as shown in fig. 1, where the method for simulating the fluid channeling law of the first interface of the well cementation cement sheath comprises:
s101, obtaining well cementation data of a well cementation cement sheath to be tested;
s102, establishing a cross flow model of a first interface of the well cementing cement sheath according to the well cementing data, wherein the cross flow model comprises a geometric model and a grid model;
s103, determining a crossflow rule theoretical formula of a crossflow model according to the geometric model;
s104, dividing the grid number of the grid model and setting the cross flow condition of the cross flow model;
s105, obtaining a numerical simulation result of the height change rule of the first interface crossflow according to the crossflow model.
Specifically, in the numerical simulation method of the first interface channeling law of the well cementation cement sheath provided by the invention, firstly, well cementation data of the well cementation cement sheath to be tested are obtained; then establishing a cross flow model of a first interface of the well cementing cement sheath according to the well cementing data, wherein the cross flow model comprises a geometric model and a grid model; determining a crossflow rule theoretical formula of a crossflow model according to the geometric model; dividing the grid number of the grid model and setting the cross flow condition of the cross flow model; and finally, obtaining a numerical simulation result of the height change rule of the first interface crossflow according to the crossflow model. In summary, the invention realizes the technical effect of high-precision simulation analysis on the fluid channeling law of the first interface of the well cementation cement sheath under the conditions of high temperature and high pressure by constructing the fluid channeling model of the first interface of the well cementation cement sheath and simulating each structure of the well cementation cement sheath by using the geometric model and modeling the stress state of each structure of the well cementation cement sheath by using the grid model.
In a specific embodiment of the present invention, obtaining well cementation data of a well cementation cement sheath to be tested includes:
acquiring material parameters and well cementation interface parameters of a well cementation cement sheath to be tested, wherein the material parameters comprise the geometric dimension, young modulus and Poisson's ratio of a well cementation material; the well cementation interface parameters comprise shear strength, hydraulic cement strength, rigidity coefficient and critical fracture energy.
Specifically, the material parameters and the well cementation interface parameters of the well cementation cement sheath well cementation material to be tested can be obtained through measuring, testing, inquiring the material specification and the like, wherein the well cementation material comprises three parts of a casing, a cement sheath and a stratum, and the material parameters comprise the geometric dimension, the Young modulus and the conifer ratio; the parameters of the well cementation interface, such as shear strength, hydraulic cement strength, rigidity coefficient and critical fracture energy.
It will be appreciated that the material parameters and cementing interface parameters of the cementing cement sheath, including but not limited to those described above, may be derived as desired.
In a specific embodiment of the present invention, establishing a fluid channeling model of a first interface of a cementing cement sheath according to the cementing data includes:
building a geometric model of a well cementation cement sheath according to the well cementation data;
dividing the geometric model into a plurality of basic stress analysis units by a plurality of hexahedral meshes, and constructing a mesh model of the well cementation cement sheath.
It can be understood that the constructed mesh model of the well cementation cement sheath is a channeling model of the first interface of the well cementation cement sheath.
Specifically, in the process of building a cement sheath, as shown in fig. 2, fig. 2 is a schematic structural diagram of a geometric model and a corresponding grid model of a casing portion according to an embodiment of the present invention. Firstly, building a geometric model of a well cementing cement ring according to geometric dimensions and other material parameters in well cementing data, and simultaneously, in order to simulate complex conditions of the well cementing cement ring in a fluid channeling process, using hexahedral meshes as basic stress analysis units to carry out stress analysis on the fluid channeling model, and dividing the geometric model into a mesh model consisting of a plurality of hexahedral meshes.
In a specific embodiment of the present invention, determining a theoretical formula of a cross-flow rule of a cross-flow model according to the geometric model includes:
analyzing the state of the fluid channeling of the cement sheath to be tested based on the Navigator equation and the laminar flow assumption, and determining the flow law of the fluid channeling in the cement sheath;
determining a flow velocity distribution type in the annular gap based on a slip-free boundary condition of the annular gap and according to a flow law of a cross flow in the well cementation cement ring;
the flow velocity in the annular gap is distributed and integrated to obtain a flow expression of the annular gap, and a flow approximation expression in the annular gap is obtained according to the law of three times;
obtaining a flow rule expression of the microcosmic-layer gas in the annular space based on the Fuchshimer equation;
and obtaining a crossflow rule theoretical formula of the crossflow model according to the flow velocity distribution type, the flow similarity expression and the flow rule expression in the annular gap.
Specifically, when a crossflow law theoretical formula of a crossflow model is determined, firstly, according to an established geometric model and based on a Navier-Stokes model equation (Navistos equation, a motion equation describing conservation of viscous incompressible fluid momentum) and laminar flow assumption, the flow law of gas at a cement sheath interface is analyzed. From the analysis it can be assumed that the pressure gradient in the annular gap is uniformly distributed parallel to the x-axis, the flow rate only exists in the x-direction, and the Navier-Stokes equation is simplified as:
wherein u is x (r) represents the flow rate of the gas in the annulus as the radius r varies,represents the pressure gradient and μ represents the gas viscosity.
The laminar flow is a flow state of the fluid, and the laminar flow is assumed to be a flow state of the fluid when the fluid flows in the pipe at a low speed, and the fluid moves smoothly and linearly in a direction parallel to the pipe axis. The flow rate of the fluid is greatest at the center of the tube and smallest at the wall of the tube.
Based on laminar flow assumptions, the fluid has a radius r within the annulus 1 And annular gap outer radius r 2 The slip-free boundary condition exists, namely the speed or the relative speed of the fluid on the wall surface is 0, and the speed distribution in the annular space can be calculated:
integrating the velocity distribution in the annular space to obtain the gas flow Q in the annular space x ：
And the approximate value of the gas flow in the annular space is obtained through three times of law calculation:
wherein the aperture h of the annular gap is (r) 2 -r 1 ) The inner diameter perimeter w of the annular gap is 2 pi r 1 The third law refers to that for a circular ring-shaped pore with a small pore diameter, the pore can be approximately regarded as rectangular calculation, and for a sleeve size and a gap size under typical field conditions, the error of the third law is less than 1 percent.
For the flow law that the gas follows in the annulus, i.e. there is a viscous (darcy) flow and also a non-linear (inertial) flow, it can be determined by the forshheimer equation (the foshimer equation, which is considered to be one of the most efficient formulas for calculating the non-darcy percolation problem), which describes the flow law that the gas follows in the annulus from a microscopic level:
where Q is the volumetric flow rate, K is the permeability, A is the cross-sectional area involved in the flow, μ is the fluid viscosity, β is the inertia coefficient, ρ is the fluid density.
According to the flow velocity distribution type, the flow approximate expression and the flow law in the annular gap, a channeling law theoretical formula for the steady-state isothermal state airflow can be obtained:
where M is the molecular weight of the gas, L is the length of the cross flow, z is the gas compression factor, R is the universal gas constant, T is the temperature, and subscripts u and d represent upstream and downstream, respectively.
In a specific embodiment of the present invention, the mesh number division of the mesh model includes:
setting an error standard threshold, carrying out grid sensitivity analysis on the grid model, determining the minimum grid number of which the error does not exceed the error standard threshold in the process of calculating the cross flow height rule numerical simulation, and carrying out grid number division according to the minimum grid number and a preset division principle.
Specifically, when a geometric model is divided by using a plurality of hexahedral meshes to build a mesh model, the division of the number of meshes has a larger influence on the calculation of the model, the more the number of meshes is, the more accurate the calculation result is, but the longer the calculation time is. Therefore, when dividing the grid number, an error standard threshold is set, and the minimum grid number of which the error of the calculated result does not exceed the error standard threshold is determined through grid sensitivity analysis. And according to the obtained minimum grid number and the importance degree of each part of the cross-flow model, carrying out grid division on each part. For example, the sleeve is the part of the analysis that is the focus of the analysis, and the mesh density of the division of this part is greater.
In a specific embodiment of the present invention, setting a crossflow condition of the crossflow model includes:
setting initial conditions of a model, wherein the initial conditions comprise applying an in-situ stress field to a model body of the channeling model, and defining an initial damage unit and an injection point boundary;
setting boundary conditions of the model, wherein the boundary conditions comprise displacement constraints of the geometric model in the three-dimensional direction;
based on the ground stress balance analysis step, eliminating the initial displacement of the model body in the channeling model, and applying the initial stress of the cement sheath interface;
and determining the ground stress according to the initial stress and the crossflow model, and applying the ground stress to the crossflow model.
Specifically, a cross-flow condition process of a cross-flow model is set, and initial conditions of the model are set first, wherein the initial conditions comprise applying an in-situ stress field to the whole model, and defining an initial damage unit and an injection point boundary. The in-situ stress field can be obtained by inquiring geological data, the initial damage unit refers to a damaged position of the first interface, the cohesive energy of the initial damage unit is set to be zero, and the boundary of the injection point is the initial boundary position of the injection point of the cross flow gas of the damaged position.
And setting boundary conditions, wherein the boundary conditions comprise displacement constraint conditions in all directions in the three-dimensional coordinate axis.
According to the Geostatic analysis step, initial displacement of the casing-cement sheath-formation is eliminated and initial stress of the cement sheath interface is applied. And determining the ground stress applied to the model according to the initial stress, the well cementation structure in the channeling model and the like.
The ground stress means an internal stress effect caused by a crustal substance due to a geologic structure or the like, and the ground stress motion balance analysis step Geostatic is a method for balancing ground stress in geologic analysis.
In a specific embodiment of the present invention, obtaining a numerical simulation result of a height change rule of a cross flow according to the cross flow model includes:
determining an influence parameter influencing the height of the fluid channeling of a first interface of the well cementing ring according to a fluid channeling rule theoretical formula of the fluid channeling model;
and obtaining a numerical simulation result of the influence rule of the influence parameter on the height of the crossflow according to the crossflow model.
Specifically, according to the theoretical formula of the obtained fluid channeling law, the influence parameters influencing the fluid channeling height of the first interface of the well cementing cement sheath are analyzed and determined. The crossflow model simulates the interface characteristics and the interface crossflow characteristics of a first interface through a Cohesive unit (Cohesive force unit) to obtain a crossflow rule that crossflow gas flows into a Cohesive viscous layer after the crossflow gas breaks through an injection point, and obtains a numerical simulation result of the influence rule of each influence parameter on the crossflow height through numerical simulation.
In a specific embodiment of the present invention, obtaining a numerical simulation result of the influence rule of the influence parameter on the height of the cross-flow according to the cross-flow model includes:
and taking any parameter in the influence parameters as a single variable to carry out simulation input to obtain a corresponding crossflow height, and obtaining a numerical simulation result of the influence rule of the influence parameters on the crossflow height according to the crossflow height.
Specifically, in one embodiment of the invention, three influencing factors of fluid viscosity, interface cementing strength and fluid channeling flow are taken as examples, and the influence relation of the fluid viscosity, the interface cementing strength and the fluid channeling flow on the first interface fluid channeling height of the well cementation cement sheath is simulated through a single variable value.
Wherein, the viscosity change of the fluid-channeling gas is larger under the high-temperature and high-pressure condition, when the influence of the fluid viscosity on the fluid-channeling height is simulated, the interface cementing strength is set to be 0.001MPa, the fluid-channeling time is set to be 0.7s, and the results are shown in Table 1.
TABLE 1 influence of fluid viscosity on fluid height of channeling
As can be seen from table 1, the cement sheath first interface fluid height decreases with increasing fluid viscosity over the same period of time.
After the well cementation cement paste is solidified to form a cement ring, the temperature and pressure in the sleeve are changed due to the influence of various working conditions of subsequent well completionComplicated, the cementing quality is deteriorated, and the interfacial cementing strength is reduced. To study the effect of the bond strength on the height of the fluid channeling, six bond strengths of 0.001MPa, 0.1MPa, 0.2MPa, 0.3MPa, 0.4MPa and 0.5MPa were set, and the fluid viscosity was set to 1X 10 -5 Pa.s, the channeling flow rate is set to be 1 multiplied by 10 -5 m 3 The cross-flow time was 0.7s, and the results were shown in Table 2.
TABLE 2 influence of interfacial bond strength on fluid-channeling height
As can be seen from table 2, the cross-flow height increases exponentially with decreasing interfacial bond strength at the same time and cross-flow rate, indicating that interfacial bond strength plays an important role in preventing cross-flow.
In addition, in order to study the influence of the cross-flow on the cross-flow height, a fluid viscosity of 1×10 was set -5 And under the condition of Pa.s and 0.001MPa of interface cementing strength, the cross flow height results of different cross flow flows are shown in table 3.
TABLE 3 influence of crossflow flow flux on crossflow height
As can be seen from table 3, under the condition that the fluid viscosity and the interfacial bond strength are unchanged, the influence of the fluid channeling flow on the fluid channeling height is larger, and the fluid channeling height increases in an exponential function form along with the increase of the fluid channeling.
The numerical simulation results of the three comprehensive influence parameters can show that the fluid property, the interface cementing strength and the fluid channeling flow rate have important influence on the fluid channeling height of the first interface of the well cementation cement sheath, and the influence relationship between the interface cementing strength and the fluid channeling flow rate is in an exponential function form.
In order to better implement the numerical simulation method of the cement sheath first interface channeling law of the cement sheath first interface in the embodiment of the invention, on the basis of the numerical simulation method of the cement sheath first interface channeling law of the cement sheath first interface, the invention also provides a numerical simulation device 300 of the cement sheath first interface channeling law of the cement sheath first interface, as shown in fig. 3, comprising:
the data acquisition unit 301 is configured to acquire well cementation data of a well cementation cement sheath to be tested;
the model building unit 302 is configured to build a fluid-channeling model of the first interface of the well cementing cement sheath according to the well cementing data, where the fluid-channeling model includes a geometric model and a grid model;
a model analysis unit 303, configured to determine a crossflow rule theoretical formula of a crossflow model according to the geometric model;
a condition setting unit 304, configured to divide the grid number of the grid model and set a cross flow condition of the cross flow model;
and the simulation result output unit 305 is configured to obtain a numerical simulation result of the height change rule of the first interface cross flow according to the cross flow model.
The numerical simulation device 300 for the fluid-filled cement sheath first interface fluid-channeling law provided in the foregoing embodiment may implement the technical solution described in the foregoing embodiment of the method for simulating the fluid-channeling law of the fluid-filled cement sheath first interface, and the specific implementation principle of each module or unit may refer to the corresponding content in the foregoing embodiment of the method for simulating the fluid-channeling law of the fluid-filled cement sheath first interface, which is not described herein again.
The invention further provides an electronic device 400 based on the numerical simulation method of the fluid channeling law of the first interface of the well cementation cement sheath, as shown in fig. 4, fig. 4 is a schematic structural diagram of an embodiment of the electronic device provided by the invention, the electronic device 400 comprises a processor 401, a memory 402 and a computer program stored in the memory 402 and capable of running on the processor 401, and the numerical simulation method of the fluid channeling law of the first interface of the well cementation cement sheath is realized when the processor 401 executes the program.
As a preferred embodiment, the electronic device further includes a display 403, configured to display that the processor 401 performs the above-mentioned numerical simulation method process of the fluid channeling law of the first interface of the well cementing ring.
The processor 401 may be an integrated circuit chip with signal processing capability. The processor 401 may be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), etc.; also digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC). The disclosed methods, steps, and logic blocks in the embodiments of the present invention may be implemented or performed. A general purpose processor may also be a microprocessor or the processor may be any conventional processor or the like.
The Memory 402 may be, but is not limited to, a random access Memory (Random Access Memory, RAM), a Read Only Memory (ROM), a Secure Digital (SD Card), a Flash Card (Flash Card), etc. The memory 402 is configured to store a program, and the processor 401 executes the program after receiving an execution instruction, and the method for defining a flow disclosed in any one of the foregoing embodiments of the present invention may be applied to the processor 401 or implemented by the processor 401.
The display 403 may be an LED display, a liquid crystal display, a touch display, or the like. The display 403 is used to display various information on the electronic device 400.
It is to be appreciated that the configuration shown in fig. 4 is merely a schematic diagram of one configuration of the electronic device 400, and that the electronic device 400 may include more or fewer components than shown in fig. 4. The components shown in fig. 4 may be implemented in hardware, software, or a combination thereof.
The present invention is not limited to the above-mentioned embodiments, and any changes or substitutions that can be easily understood by those skilled in the art within the technical scope of the present invention are intended to be included in the scope of the present invention.
Claims (10)
1. The numerical simulation method for the fluid channeling law of the first interface of the well cementation cement sheath is characterized by comprising the following steps:
acquiring well cementation data of a well cementation cement sheath to be tested;
establishing a cross flow model of a first interface of the well cementing cement sheath according to the well cementing data, wherein the cross flow model comprises a geometric model and a grid model;
determining a crossflow rule theoretical formula of a crossflow model according to the geometric model;
dividing the grid number of the grid model and setting the cross flow condition of the cross flow model;
and obtaining a numerical simulation result of the height change rule of the first interface crossflow according to the crossflow model.
2. The method for numerical simulation of the fluid channeling law of the first interface of the well cementation cement sheath according to claim 1, wherein the obtaining the well cementation data of the well cementation cement sheath to be tested comprises:
acquiring material parameters and well cementation interface parameters of a well cementation cement sheath to be tested, wherein the material parameters comprise the geometric dimension, young modulus and Poisson's ratio of a well cementation material; the well cementation interface parameters comprise shear strength, hydraulic cement strength, rigidity coefficient and critical fracture energy.
3. The method for numerical simulation of the fluid channeling law of the first interface of the well cementation cement sheath according to claim 1, wherein the establishing a fluid channeling model of the first interface of the well cementation cement sheath according to the well cementation data comprises:
building a geometric model of a well cementation cement sheath according to the well cementation data;
dividing the geometric model into a plurality of basic stress analysis units by a plurality of hexahedral meshes, and constructing a mesh model of the well cementation cement sheath.
4. The method for modeling the fluid-channeling law value of a first interface of a well cementation cement sheath according to claim 1, wherein determining a fluid-channeling law theoretical formula of a fluid-channeling model according to the geometric model comprises:
analyzing the state of the fluid channeling of the cement sheath to be tested based on the Navigator equation and the laminar flow assumption, and determining the flow law of the fluid channeling in the cement sheath;
determining a flow velocity distribution type in the annular gap based on a slip-free boundary condition of the annular gap and according to a flow law of a cross flow in the well cementation cement ring;
the flow velocity in the annular gap is distributed and integrated to obtain a flow expression of the annular gap, and a flow approximation expression in the annular gap is obtained according to the law of three times;
obtaining a flow rule expression of the microcosmic-layer gas in the annular space based on the Fuchshimer equation;
and obtaining a crossflow rule theoretical formula of the crossflow model according to the flow velocity distribution type, the flow similarity expression and the flow rule expression in the annular gap.
5. The method for numerical simulation of the fluid channeling law of the first interface of the well cementation cement sheath according to claim 1, wherein the step of dividing the grid model into the number of grids comprises the steps of:
setting an error standard threshold, carrying out grid sensitivity analysis on the grid model, determining the minimum grid number of which the error does not exceed the error standard threshold in the process of calculating the cross flow height rule numerical simulation, and carrying out grid number division according to the minimum grid number and a preset division principle.
6. The method for numerical simulation of a fluid channeling law of a first interface of a well cementation cement sheath according to claim 1, wherein the setting of fluid channeling conditions of the fluid channeling model comprises:
setting initial conditions of a model, wherein the initial conditions comprise applying an in-situ stress field to a model body of the channeling model, and defining an initial damage unit and an injection point boundary;
setting boundary conditions of the model, wherein the boundary conditions comprise displacement constraints of the geometric model in the three-dimensional direction;
based on the ground stress balance analysis step, eliminating the initial displacement of the model body in the channeling model, and applying the initial stress of the cement sheath interface;
and determining the ground stress according to the initial stress and the crossflow model, and applying the ground stress to the crossflow model.
7. The method for numerically simulating the cross flow law of the first interface of the well cementation cement sheath according to claim 1, wherein the step of obtaining the numerical simulation result of the cross flow height change law according to the cross flow model comprises the following steps:
determining an influence parameter influencing the height of the fluid channeling of a first interface of the well cementing ring according to a fluid channeling rule theoretical formula of the fluid channeling model;
and obtaining a numerical simulation result of the influence rule of the influence parameter on the height of the crossflow according to the crossflow model.
8. The method for numerically simulating the law of fluid channeling on a first interface of a well cementation cement sheath according to claim 7, wherein the obtaining, according to the fluid channeling model, a numerical simulation result of the law of influence of the influencing parameter on the fluid channeling height comprises:
and taking any parameter in the influence parameters as a single variable to carry out simulation input to obtain a corresponding crossflow height, and obtaining a numerical simulation result of the influence rule of the influence parameters on the crossflow height according to the crossflow height.
9. The utility model provides a well cementation cement sheath first interface channeling law numerical simulation device which characterized in that includes:
the data acquisition unit is used for acquiring well cementation data of the well cementation cement sheath to be tested;
the model construction unit is used for building a cross flow model of a first interface of the well cementing cement sheath according to the well cementing data, and the cross flow model comprises a geometric model and a grid model;
the model analysis unit is used for determining a crossflow rule theoretical formula of the crossflow model according to the geometric model;
the condition setting unit is used for dividing the grid number of the grid model and setting the cross flow condition of the cross flow model;
and the simulation result output unit is used for obtaining a numerical simulation result of the height change rule of the first interface crossflow according to the crossflow model.
10. An electronic device comprising a memory and a processor, wherein,
the memory is used for storing programs;
the processor is coupled to the memory, and is configured to execute the program stored in the memory, so as to implement the steps in the method for simulating the fluid channeling law value of the first interface of the cement sheath of the well cementation cement sheath according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310465253.5A CN116542037A (en) | 2023-04-26 | 2023-04-26 | Numerical simulation method and device for fluid channeling law of first interface of well cementing cement sheath |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310465253.5A CN116542037A (en) | 2023-04-26 | 2023-04-26 | Numerical simulation method and device for fluid channeling law of first interface of well cementing cement sheath |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116542037A true CN116542037A (en) | 2023-08-04 |
Family
ID=87442876
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310465253.5A Pending CN116542037A (en) | 2023-04-26 | 2023-04-26 | Numerical simulation method and device for fluid channeling law of first interface of well cementing cement sheath |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116542037A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117174218A (en) * | 2023-11-02 | 2023-12-05 | 济南大学 | Cement strength conditional probability distribution estimation method and system based on flow model |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106596400A (en) * | 2016-12-29 | 2017-04-26 | 中国地质大学（北京） | Evaluation device and experimental method for cementation quality of second interface for well cementation |
| CN110424947A (en) * | 2019-07-31 | 2019-11-08 | 长江大学 | A method of measurement the first interface microannulus size of well cementation |
| CN113756744A (en) * | 2021-09-29 | 2021-12-07 | 中海石油(中国)有限公司 | Numerical simulation method for sealing failure of deep-water high-temperature and high-pressure well casing-cement sheath-stratum system |
| CN113898312A (en) * | 2021-09-29 | 2022-01-07 | 中海石油(中国)有限公司 | Numerical simulation method and system for annular oil-gas channeling of well cementation of deepwater high-temperature and high-pressure well and storage medium |
| CN115455770A (en) * | 2022-09-14 | 2022-12-09 | 东北石油大学 | Oil-gas field wellbore structure stress fine simulation method based on substructure |
-
2023
- 2023-04-26 CN CN202310465253.5A patent/CN116542037A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106596400A (en) * | 2016-12-29 | 2017-04-26 | 中国地质大学（北京） | Evaluation device and experimental method for cementation quality of second interface for well cementation |
| CN110424947A (en) * | 2019-07-31 | 2019-11-08 | 长江大学 | A method of measurement the first interface microannulus size of well cementation |
| CN113756744A (en) * | 2021-09-29 | 2021-12-07 | 中海石油(中国)有限公司 | Numerical simulation method for sealing failure of deep-water high-temperature and high-pressure well casing-cement sheath-stratum system |
| CN113898312A (en) * | 2021-09-29 | 2022-01-07 | 中海石油(中国)有限公司 | Numerical simulation method and system for annular oil-gas channeling of well cementation of deepwater high-temperature and high-pressure well and storage medium |
| CN115455770A (en) * | 2022-09-14 | 2022-12-09 | 东北石油大学 | Oil-gas field wellbore structure stress fine simulation method based on substructure |
Non-Patent Citations (2)
| Title |
|---|
| JOHN C. STORMONT ET AL: "Gas flow through cement-casing microannuli under varying stress conditions", GEOMECHANICS FOR ENERGY AND THE ENVIRONMENT, 26 December 2017 (2017-12-26), pages 1 - 12 * |
| 曲路路: "高温高压井套管-水泥环-地层组合体应力变化规律研究", 中国优秀硕士学位论文全文数据库工程科技Ⅰ辑, no. 02, 15 February 2021 (2021-02-15), pages 019 - 160 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117174218A (en) * | 2023-11-02 | 2023-12-05 | 济南大学 | Cement strength conditional probability distribution estimation method and system based on flow model |
| CN117174218B (en) * | 2023-11-02 | 2024-01-26 | 济南大学 | Cement strength conditional probability distribution estimation method and system based on flow model |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111794727B (en) | Pump injection frequency selection method and device for pulse circulation hydraulic fracturing | |
| CN116542037A (en) | Numerical simulation method and device for fluid channeling law of first interface of well cementing cement sheath | |
| CN105913494A (en) | Multi-scale fracture fine geological modeling and value simulation method and device | |
| CN110671104A (en) | Interpretation method of interwell parameters of fracture-cavity type oil reservoir based on interference well testing interpretation | |
| CN106709219A (en) | Area initial stress field inversion method and device under complicated geological conditions | |
| Qadeer et al. | Correcting oil/water relative permeability data for capillary end effect in displacement experiments | |
| CN108227032B (en) | Determination method and device for ground stress | |
| CN116401966A (en) | Simulation method for time-varying slurry diffusion process of porous fractured rock mass | |
| CN104990853B (en) | The Forecasting Methodology of the full rank permeability tensor of porous media | |
| CN107169227B (en) | A kind of the coarse grid analogy method and system of staged fracturing horizontal well | |
| CN116882221B (en) | Geothermal exploitation numerical simulation method and system based on three-dimensional fracture type thermal storage model | |
| Ferreira et al. | Improved SWMM modeling for rapid pipe filling incorporating air behavior in intermittent water supply systems | |
| Kolditz et al. | Quasi-steady-state strategy for numerical simulation of geothermal circulation in hot dry rock fractures | |
| Shi et al. | Large eddy simulations of model-scale turbulent atmospheric boundary layer flows | |
| CN110377990B (en) | Method for solving surrounding rock stress and displacement of shallow tunnel in stratum containing cavity under action of gravity | |
| CN110857625B (en) | Method, device and equipment for acquiring reserve of closed water body fracture-cave unit of carbonate rock | |
| CN112377184B (en) | Physical property parameter analysis method and device for beaded fracture-cavity type carbonate reservoir | |
| CN106930759A (en) | A kind of Forecasting Methodology of low permeability sandstone reservoir production capacity | |
| CN115563894A (en) | Permeability prediction method for complex layered and spliced digital core | |
| Guo et al. | Finite-volume multi-stage scheme for advection-diffusion modeling in shallow water flows | |
| CN114462278A (en) | Method and system for analyzing production induced stress before shale gas reservoir repeated fracturing | |
| CN107965308B (en) | Method and device for determining water yield of single well | |
| Zheng et al. | Hydro-mechanical modeling of impermeable discontinuity in rock by extended finite element method | |
| CN117421939B (en) | Shale oil fracture system simulation agent method based on track piecewise linearization | |
| RU2745142C1 (en) | Method and system for modeling hydraulic fractures of a formation of finite/infinite conductivity and lateral-longtitudinal location relative to a horizontal wellbore |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |