US20230266498A1 - Full wellbore pressure calculation method, device and computer readable storage medium for fracturing horizontal well - Google Patents
Full wellbore pressure calculation method, device and computer readable storage medium for fracturing horizontal well Download PDFInfo
- Publication number
- US20230266498A1 US20230266498A1 US17/890,248 US202217890248A US2023266498A1 US 20230266498 A1 US20230266498 A1 US 20230266498A1 US 202217890248 A US202217890248 A US 202217890248A US 2023266498 A1 US2023266498 A1 US 2023266498A1
- Authority
- US
- United States
- Prior art keywords
- fracture
- pressure
- matrix
- flow
- unit
- 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.)
- Granted
Links
- 238000003860 storage Methods 0.000 title claims abstract description 13
- 238000004364 calculation method Methods 0.000 title claims abstract description 10
- 238000004519 manufacturing process Methods 0.000 claims abstract description 53
- 230000000977 initiatory effect Effects 0.000 claims abstract description 25
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 18
- 230000008859 change Effects 0.000 claims abstract description 12
- 206010017076 Fracture Diseases 0.000 claims description 124
- 208000010392 Bone Fractures Diseases 0.000 claims description 107
- 239000011159 matrix material Substances 0.000 claims description 73
- 239000007788 liquid Substances 0.000 claims description 35
- 208000013201 Stress fracture Diseases 0.000 claims description 33
- 239000012530 fluid Substances 0.000 claims description 13
- 230000005484 gravity Effects 0.000 claims description 13
- 239000012071 phase Substances 0.000 claims description 11
- 230000035699 permeability Effects 0.000 claims description 8
- 230000009977 dual effect Effects 0.000 claims description 5
- 239000007791 liquid phase Substances 0.000 claims description 5
- 238000013178 mathematical model Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000000704 physical effect Effects 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 238000009826 distribution Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G01V99/005—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Definitions
- the invention relates to the technical field of numerical simulation of oil and gas reservoirs, in particular to a method, device and computer-readable storage medium for calculating a full wellbore pressure for fracturing horizontal wells.
- Unconventional shale gas reservoirs, tight gas reservoirs and other original reservoirs generally have low porosity and permeability, and low natural depletion production.
- the production of each artificial fracture is equal, and the pressure loss in the wellbore is ignored, and the fracturing horizontal well is regarded as having infinite conductivity.
- the actual situation is that due to the heterogeneity of the reservoir and the interference between fractures in the production process, the production of each artificial fracture is not equal.
- the flow in the horizontal section of the wellbore belongs to variablemass multi-phase flow, which will generate frictional pressure drop, gravity pressure drop, and accelerated pressure drop caused by fluid confluence in fractures.
- the pressure change in the wellbore will in turn affect the fluid confluence of each fracture, resulting in the change of the productivity of each fracture. Therefore, it is necessary to study the coupled formation wellbore flow model of fracturing horizontal wells to achieve accurate prediction of fracture production, reservoir pressure variation and the pressure distribution characteristics of the full wellbore in fracturing horizontal wells.
- the purpose of the present invention is to provide a method for calculating the full wellbore pressure of a fracturing horizontal well.
- the technical solution of the present invention is as follows.
- a method for calculating a full wellbore pressure of a fracturing horizontal well includes the following steps:
- S1 collect wellbore structural parameters, fluid physical properties, reservoir parameters and production history data at different time steps, and establish a fully implicit numerical model of formation flow based on embedded discrete fracture model (EDFM).
- EDFM embedded discrete fracture model
- a device for calculating the full wellbore pressure of a fracturing horizontal well includes a processor, an acquisition module for acquiring initial data for calculating the full wellbore pressure of a fracturing horizontal well, an output module for outputting calculation results, and a storage module; a program for calculating the full wellbore pressure of a fracturing horizontal well that can be run on the processor is stored thereon, and the program for calculating the full wellbore pressure of a fracturing horizontal well is executed by the processor to realize the above-mentioned method steps.
- a computer readable storage medium having stored therein program code executable by a processor, the computer readable storage medium including a plurality of instructions configured to cause the processor to perform the described wellbore pressure calculation method of fracturing horizontal well.
- the present invention has the following beneficial effects:
- FIG. 1 is a schematic diagram of the discrete spatialgrid of the gas well based on EDFM in the embodiment of the present invention
- FIG. 2 is reservoir pressure distribution at different production times of the gas well in the embodiment of the present invention.
- FIG. 3 is the prediction curve of the full wellbore pressure distribution at different production times of the gas well according to the embodiment of the present invention
- FIG. 4 is a flow chart of the method for calculating the full wellbore pressure of a fracturing horizontal well according to an embodiment of the present invention
- FIG. 4 A is a flow chart of establishment of a formation flow fully implicit numerical model according to an embodiment of the present invention
- FIG. 4 B is a flow chart of obtaining a wellbore pressure variation of each fracturing section of the horizontal well according to an embodiment of the present invention.
- FIG. 5 is a device for calculating the full wellbore pressure of a fracturing horizontal well for the embodiment of the present invention.
- FIG. 4 To calculate the full wellbore pressure of fracturing horizontal wells, the flow chart of this calculation method can be referred to FIG. 4 , including the following steps:
- S1 collect wellbore structural parameters, fluid physical properties, reservoir parameters and production history data at different time steps, and establish a fully implicit numerical model of formation flow based on embedded discrete fracture model (EDFM);
- EDFM embedded discrete fracture model
- EDFM embedded discrete fracture model
- q l,mm represents the flow exchange between connected matrix grids, the unit is 1/day.
- the finite difference method is used to derivative expansion of the flow exchange term ⁇ q l,mm ⁇ q l,fm on the left side of the above seepage equation, if there is flow exchange between i and j grids, the fully implicit iterative form of the flow exchange between the two grids is.
- n+1 represents the next time step;
- v represents the v-th iteration step;
- ⁇ p lmi , ⁇ p lmj represents the iterative pressure variables in the i-th and j-th grids of the matrix system;
- a lmi a lmj representing coefficients of iterative pressure variables in the i-th and j-th grids of the matrix system;
- ⁇ S lmi , ⁇ S lmj represents the iterative saturation variable in the i-th and j-th grids of the matrix system;
- b lmi , b lmj representing coefficients the iterative saturation variables in the i-th and j-th grids of the matrix system;
- d lmij is the known quantity of flow exchange between the i-th and j-th grids of the matrix system.
- mm represents the coefficient matrix in the matrix flow equation
- ⁇ m represents the unknown quantity matrix in the matrix flow equation
- d m is the matrix of known quantities in the matrix flow equation. Therefore, the fully implicit solution matrix based on the finite difference method for micro-fractures, artificial fractures, the flow exchange terms of each system and the internal flow of the wellbore can be deduced as follows:
- ff, FF, and ww represent the coefficient matrices in the microfractures, artificial fractures, and wellbore flow equations, respectively;
- mf, fm represent the flow exchange coefficient matrices between the matrix and the microfracture system;
- fF, Ff represent the microfractures and artificial fractures;
- Fw represents the flow exchange coefficient matrix between the artificial fracture and the wellbore;
- d f , d F , d w represent the known quantity matrices in the microfracture, artificial fracture, and wellbore flow equations, respectively;
- ⁇ f, ⁇ F, and ⁇ w represent the unknown quantity matrices in the microfracture, artificial fracture, and wellbore flow equations, respectively.
- Equation (8) the gas production and liquid production of each fracture, the grid block pressure where each fracture is located, and the initial bottom-hole flow pressure can be obtained.
- v sl represents the apparent flow rate of the liquid phase, the unit is m/s;
- v sg represents the apparent flow rate of the gas phase, the unit is m/s;
- q sc represents the fracture gas production, the unit m 3 /s;
- q w represents the fracture fluid production, the unit is m 3 /s;
- B g represents the volume coefficient of the gas phase, dimensionless;
- B w represents the liquid phase volume coefficient, dimensionless;
- A represents the cross-sectional area of the wellbore, the unit is m 2 .
- H l exp [ ( c 1 + c 2 ⁇ sin ⁇ ⁇ + c 3 ⁇ sin 2 ⁇ ⁇ + c 4 ⁇ N l 2 ) ⁇ N g ⁇ v c 5 N lv c 6 ] ( 10 )
- N lv , N gv , N l in the above formula are as follows:
- ⁇ l reprents the liquid density the unit is kg/m 3 ;
- ⁇ represents the gas-liquid interfacial tension, the unit is N/m;
- ⁇ l reprents the liquid viscosity the unit is Pa ⁇ s.
- (dp/dL) represents the pressure gradient, the unit is MPa/m; dL represents the length of the horizontal section between adjacent fractures, the unit is m; ⁇ m represents the mixed density, the unit is kg/m 3 ; G represents the acceleration of gravity, the unit is m/s 2 ; ⁇ represents the angle between the wellbore and the horizontal direction; f m represents the two-phase friction coefficient, dimensionless; v m represents the mixing speed, the unit is m/s; D represents the diameter of the pipe, the unit is m.
- N F represents the total number of the fractures
- ⁇ L i represents the distance from the i-th fracture initiation point to the bottom-hole of the well, the unit is m
- ⁇ p wi represents the horizontal wellbore pressure change from the bottom-hole to the i-th fracture, the unit is MPa.
- step S5 bring the bottom-hole flow pressure and the bottom-hole pressure of each fracture initiation point into the formation flow fully implicit numerical model based on the embedded discrete fracture model (EDFM) in step S2, and iteratively cycle steps S2-S4 until the iterative variables converge, obtain the bottom-hole flow pressure and bottom-hole pressure of each fracture initiation point at this time step;
- EDFM embedded discrete fracture model
- a device for calculating the full wellbore pressure of a fracturing horizontal well includes a processor, an acquisition module for acquiring initial data for calculating the full wellbore pressure of a fracturing horizontal well, an output module and a storage module for outputting calculation results; a program for calculating the full wellbore pressure of a fracturing horizontal well is stored on the storage module, which can be run on the processor, and the program for calculating the full wellbore pressure of a fracturing horizontal well is executed by the processor to implement the above method step.
- a computer readable storage medium having stored therein program code executable by a processor, the computer readable storage medium including a plurality of instructions configured to cause the processor to perform the described fracturing horizontal well process. wellbore pressure calculation method.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a full wellbore pressure calculation method, device and computer-readable storage medium for fracturing horizontal wells, includes steps: collect basic parameters, establish fully implicit numerical model of formation flow; solve the full implicit numerical model under an inner boundary condition of a constant gas production rate; calculate wellbore pressure change of each fracturing section of horizontal well; calculate the bottom-hole pressure of each fracture initiation point; repeat steps until the variables converge, the bottom-hole flow pressure and the bottom-hole pressure at each fracture initiation point are obtained at this time step; the wellbore pressure and wellhead casing pressure are calculated at this time step. The method is for predicting the variation law of full wellbore pressure and wellhead casing pressure by using production data in the absence of bottom-hole pressure test data, and has practical value for accurate prediction of production performance of gas reservoir fracturing horizontal wells.
Description
- The invention relates to the technical field of numerical simulation of oil and gas reservoirs, in particular to a method, device and computer-readable storage medium for calculating a full wellbore pressure for fracturing horizontal wells.
- Unconventional shale gas reservoirs, tight gas reservoirs and other original reservoirs generally have low porosity and permeability, and low natural depletion production. In order to obtain industrial production capacity, it is necessary to generate multiple artificial fractures through long horizontal section plus hydraulic fracturing technology to improve the seepage capacity of unconventional reservoirs. In the past research, it is generally believed that the production of each artificial fracture is equal, and the pressure loss in the wellbore is ignored, and the fracturing horizontal well is regarded as having infinite conductivity. However, the actual situation is that due to the heterogeneity of the reservoir and the interference between fractures in the production process, the production of each artificial fracture is not equal. In addition, due to the confluence of reservoir fluid, the flow in the horizontal section of the wellbore belongs to variablemass multi-phase flow, which will generate frictional pressure drop, gravity pressure drop, and accelerated pressure drop caused by fluid confluence in fractures. The pressure change in the wellbore will in turn affect the fluid confluence of each fracture, resulting in the change of the productivity of each fracture. Therefore, it is necessary to study the coupled formation wellbore flow model of fracturing horizontal wells to achieve accurate prediction of fracture production, reservoir pressure variation and the pressure distribution characteristics of the full wellbore in fracturing horizontal wells.
- In view of the above reasons, the purpose of the present invention is to provide a method for calculating the full wellbore pressure of a fracturing horizontal well. The technical solution of the present invention is as follows.
- A method for calculating a full wellbore pressure of a fracturing horizontal well, includes the following steps:
- S1, collect wellbore structural parameters, fluid physical properties, reservoir parameters and production history data at different time steps, and establish a fully implicit numerical model of formation flow based on embedded discrete fracture model (EDFM).
- S2, solve the formation flow fully implicit numerical model based on the embedded discrete fracture model (EDFM) under an inner boundary condition of a constant gas production rate at a certain time step, and obtain the gas production, liquid production, grid block pressure and initial well bottom flow pressure of each fracture.
- S3, bring the gas production volume, liquid production volume and the pressure of the grid block where the fracture initiation point is located into the pressure drop model considering the gravity loss, frictional resistance loss and fracture convergence loss, and obtain the wellbore pressure variation of each fracturing section of the horizontal well.
- S4, calculate the well bottom-hole pressure of each fracture initiation point from the initial well bottom-hole flow pressure and the wellbore pressure change of each fracturing section.
- S5, bring the bottom-hole flow pressure and the bottom-hole pressure of each fracture initiation point into the formation flow fully implicit numerical model based on the embedded discrete fracture model (EDFM) in the step S2, and iteratively cycle steps S2-S4 until the iterative variables converge, obtain the bottom-hole flow pressure and bottom-hole pressure of each fracture initiation point at this time step.
- S6, bring the bottom-hole flow pressure and the total gas production and liquid production at the time step into the pressure drop model considering gravity loss and frictional resistance loss in the inclined section and vertical well section, and obtain each section wellbore pressure and wellhead casing pressure through the iterative cycle step by step.
- S7, for each time step, repeat steps S2-S6 to obtain the full wellbore pressure at different production times.
- A device for calculating the full wellbore pressure of a fracturing horizontal well, includes a processor, an acquisition module for acquiring initial data for calculating the full wellbore pressure of a fracturing horizontal well, an output module for outputting calculation results, and a storage module; a program for calculating the full wellbore pressure of a fracturing horizontal well that can be run on the processor is stored thereon, and the program for calculating the full wellbore pressure of a fracturing horizontal well is executed by the processor to realize the above-mentioned method steps.
- A computer readable storage medium having stored therein program code executable by a processor, the computer readable storage medium including a plurality of instructions configured to cause the processor to perform the described wellbore pressure calculation method of fracturing horizontal well.
- Compared with the prior art, the present invention has the following beneficial effects:
- 1. based on the dual medium-embedded discrete fracture coupling model, a fully implicit numerical model for the seepage flow of fracturing horizontal wells in unconventional tight gas reservoirs and shale gas reservoirs is established, which realizes the accurate characterization of multi-medium flow laws in unconventional gas reservoirs;
- 2. considering the pressure drop, frictional resistance and gravity pressure drop caused by the mass-varying flow of formation fluid into the wellbore, a wellbore pressure drop model for fracturing horizontal wells in unconventional gas reservoirs with gas-liquid two-phase flow was established, and the wellbore pressure drop model for fracturing horizontal wells in unconventional gas reservoirs was established. The full wellbore pressure distribution characteristics are simulated and predicted, and in the absence of bottom-hole pressure test data, the production data can be used to predict the change law of the whole wellbore pressure and wellhead casing pressure, which saves the test cost and has important practical significance.
-
FIG. 1 is a schematic diagram of the discrete spatialgrid of the gas well based on EDFM in the embodiment of the present invention; -
FIG. 2 is reservoir pressure distribution at different production times of the gas well in the embodiment of the present invention; -
FIG. 3 is the prediction curve of the full wellbore pressure distribution at different production times of the gas well according to the embodiment of the present invention; -
FIG. 4 is a flow chart of the method for calculating the full wellbore pressure of a fracturing horizontal well according to an embodiment of the present invention;FIG. 4A is a flow chart of establishment of a formation flow fully implicit numerical model according to an embodiment of the present invention; andFIG. 4B is a flow chart of obtaining a wellbore pressure variation of each fracturing section of the horizontal well according to an embodiment of the present invention. -
FIG. 5 is a device for calculating the full wellbore pressure of a fracturing horizontal well for the embodiment of the present invention. - In order to have a clearer understanding of the technical features, purposes and beneficial effects of the present invention, the technical solutions of the present invention are now described in detail below, but should not be construed as limiting the scope of implementation of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
- To calculate the full wellbore pressure of fracturing horizontal wells, the flow chart of this calculation method can be referred to
FIG. 4 , including the following steps: - S1, collect wellbore structural parameters, fluid physical properties, reservoir parameters and production history data at different time steps, and establish a fully implicit numerical model of formation flow based on embedded discrete fracture model (EDFM);
- This example takes an actual shale gas well as an example, and the specific basic data are as follows:
-
TABLE 1 Basic data of an actual shale gas well Numerical Numerical Parameter value Parameter value Outer zone volume 1400*400*50 Formation 40 (m3) pressure(MPa) Length of main 200 Daily gas 15 fracture(m) production(10,000 cubic meters/day) Horizontal section 1000 Main fracture 500 length(m) permeability(mD) The number of 11 Micro-fracture 0.005 main fracture permeability(mD) Matrix 0.0001 Matrix porosity 5% permeability(mD) Stress Sensitivity 0.1 Micro-fracture 1% coefficient(MP−1) porosity Sleeve size(m) 0.1214 Shape factor(m−2) 0.01 - The establishment of the formation flow fully implicit numerical model based on the embedded discrete fracture model (EDFM) includes the following steps:
- S11, establish a formation flow mathematical model based on the embedded discrete fracture model, use the dual medium model to describe the flow exchange between the matrix and micro-fractures, and use the embedded discrete fracture model to describe the flow exchange between the matrix and artificial fractures, based on the mass conservation equation, the seepage differential equation is established with Darcy's law, and the formation flow mathematical model of the embedded discrete fracture model is as follows:
- Matrix system seepage equation:
-
- The seepage equation of the micro-fracture system:
-
- Seepage equation of artificial fracture system:
-
- The supplementary equation is as follows:
-
- Wherein km and kf represent the permeability of matrix, micro-fracture, and artificial fracture system, respectively, the unit is D; krlm and krlf represent the relative permeability of matrix, micro-fracture and artificial fracture system, the unit is D, wherein l=g,w; μlm, μlf, μlF represent the fluid viscosity of matrix, micro-fracture and artificial fracture system, the unit is Pa·s, wherein l=g,w; Blm, Blf, BlF represent the fluid volume coefficients of the matrix, micro-fracture, and artificial fracture systems, respectively, dimensionless, wherein l=g,w; μm, μf, μF represent the porosity of the matrix, micro-fractures, and artificial fracture systems, respectively, dimensionless; Slm, Slf, SlF represent saturation, dimensionless, wherein l=g,w; t represents production time, the unit is day; plm, plf, plF represent the two-phase pressure of matrix, micro-fracture and artificial fracture system, the unit is MPa, wherein l=g,w; rl represents the phase density, the unit is kg/m3, wherein l=g,w; G represents the acceleration of gravity, the unit is m/s2; z represents vertical depth, the unit is m; ql,fm represents the flow exchange between matrix and micro-fractures based on the dual medium model, the unit is 1/day, wherein l=g,w; α represents the shape factor between matrix and micro-fractures, the unit is m−2; ql,F represents the flow exchange between the grids after an artificial fracture is divided into multiple grids by the matrix grid, the unit is 1/day; ql,FF represents the flow exchange between different artificial fracture grids, the unit is 1/day; ql,Ff represents the flow exchange between the artificial fracture grid and the matrix grid, the unit is 1/day; ql,well represents the flow exchange between the artificial fracture grid and the connected well grid, the unit is 1/day; GF represents geometric parameters for the connection pairs of grids in the same fracture; GFF represents the geometric parameter of the connection pairs of intersected fracture grids between different fractures; GFf is the geometric parameter of the connection pair between the matrix grid and intersected fracture grid; Gwell represents geometric parameters for connected pairs between fracture grid and intersected well grid.
- S12, use a three-dimensional cubic grid to discretize the reservoir in space, based on EDFM, the fractures are embedded in the three-dimensional cubic grid in a two-dimensional plane.
- S13, use the finite difference method to establish the formation flow fully implicit numerical model based on the embedded discrete fracture model (EDFM), and solve the formation flow fully implicit numerical model; include the following steps:
- For the seepage equation (1) of the matrix system, the flow exchange relationship between the matrix grids on the left can be written in the form of flow:
-
- In the formula: ql,mm represents the flow exchange between connected matrix grids, the unit is 1/day.
- The finite difference method is used to derivative expansion of the flow exchange term Σql,mm−ql,fm on the left side of the above seepage equation, if there is flow exchange between i and j grids, the fully implicit iterative form of the flow exchange between the two grids is.
-
q lij,mm n+1 =a lmi δp lmi +a lmj δp lmj +b lmi δS lmwi +b lmj δS lmwj +d lmij (6) - In the formula: n+1 represents the next time step; v represents the v-th iteration step; δplmi, δplmj represents the iterative pressure variables in the i-th and j-th grids of the matrix system; almi, almj representing coefficients of iterative pressure variables in the i-th and j-th grids of the matrix system; δSlmi, δSlmj represents the iterative saturation variable in the i-th and j-th grids of the matrix system; blmi, blmj representing coefficients the iterative saturation variables in the i-th and j-th grids of the matrix system; dlmij is the known quantity of flow exchange between the i-th and j-th grids of the matrix system.
- The mass accumulation term
-
- on the right side of the above seepage equation also derives and expands the unknowns, the coefficients of the unknowns are added to the coefficient matrix, and the known quantities are brought into the known quantity matrix to obtain:
-
mm×δm=d m (7) - In the formula: mm represents the coefficient matrix in the matrix flow equation; δm represents the unknown quantity matrix in the matrix flow equation; dm is the matrix of known quantities in the matrix flow equation. Therefore, the fully implicit solution matrix based on the finite difference method for micro-fractures, artificial fractures, the flow exchange terms of each system and the internal flow of the wellbore can be deduced as follows:
-
- Wherein ff, FF, and ww represent the coefficient matrices in the microfractures, artificial fractures, and wellbore flow equations, respectively; mf, fm represent the flow exchange coefficient matrices between the matrix and the microfracture system; fF, Ff represent the microfractures and artificial fractures; Fw represents the flow exchange coefficient matrix between the artificial fracture and the wellbore; df, dF, dw represent the known quantity matrices in the microfracture, artificial fracture, and wellbore flow equations, respectively; δf, δF, and δw represent the unknown quantity matrices in the microfracture, artificial fracture, and wellbore flow equations, respectively.
- By solving Equation (8), the gas production and liquid production of each fracture, the grid block pressure where each fracture is located, and the initial bottom-hole flow pressure can be obtained.
- S2, solve the formation flow fully implicit numerical model based on the embedded discrete fracture model (EDFM) under an inner boundary condition of a constant gas production rate at a certain time step, and obtain the gas production, liquid production, grid block pressure and initial bottom-hole flow of each fracture; when this step is performed, in the first iteration step, the pressure drop change in the horizontal section is not considered; in the subsequent iteration steps, the pressure drop change in the horizontal section is considered.
- S3, bring the gas production, liquid production, and the pressure of the grid block where the fracture initiation point is located into the pressure drop model considering the gravity loss, frictional resistance loss and fracture convergence loss, and obtain the wellbore pressure change of each fracturing section of the horizontal well. Specifically include the following steps:
- S31, calculate the apparent flow rate of gas and liquid at the initiation point of each fracture from the gas production and liquid production of each fracture, and bring it into the Mukherjee-Brill (M-B) liquid holdup formula to calculate the liquid holdup at each fracture initiation point:
-
v sg =q sc B g /A -
v sl =q w B w /A (9) - Wherein vsl represents the apparent flow rate of the liquid phase, the unit is m/s; vsg represents the apparent flow rate of the gas phase, the unit is m/s; qsc represents the fracture gas production, the unit m3/s; qw represents the fracture fluid production, the unit is m3/s; Bg represents the volume coefficient of the gas phase, dimensionless; Bw represents the liquid phase volume coefficient, dimensionless; A represents the cross-sectional area of the wellbore, the unit is m2.
-
- In the formula: Hl represents the liquid holdup, dimensionless; c1=−0.380113; c2=−0.129875; c3=−0.119788; c4=2.343227; c5=0.475686; c6=0.288657; θ represents the angle between the wellbore and the horizontal direction, the unit is degree (°).
- The Nlv, Ngv, Nl in the above formula are as follows:
-
- In the formula: μl reprents the liquid density, the unit is kg/m3; δ represents the gas-liquid interfacial tension, the unit is N/m; μl reprents the liquid viscosity, the unit is Pa·s.
- S32, calculate the density ρm of the mixture, friction coefficient fm from the liquid holdup at each fracture initiation point and the grid block pressure, and the pressure gradient equation considering gravity loss, frictional resistance loss and fracture convergence loss is established as follows:
-
- In the formula: (dp/dL) represents the pressure gradient, the unit is MPa/m; dL represents the length of the horizontal section between adjacent fractures, the unit is m; ρm represents the mixed density, the unit is kg/m3; G represents the acceleration of gravity, the unit is m/s2; θ represents the angle between the wellbore and the horizontal direction; fm represents the two-phase friction coefficient, dimensionless; vm represents the mixing speed, the unit is m/s; D represents the diameter of the pipe, the unit is m.
- S33, combine the pressure gradient equation (12) and the distance ΔLi from the i-th fracture initiation point to the bottom-hole (i.e, the first fracture initiation point), the wellbore pressure changes in each fracturing section of the horizontal well are obtained as follows:
-
- In the formula: NF represents the total number of the fractures; ΔLi represents the distance from the i-th fracture initiation point to the bottom-hole of the well, the unit is m; Δpwi represents the horizontal wellbore pressure change from the bottom-hole to the i-th fracture, the unit is MPa.
- S4, calculate the bottom-hole pressure of each fracture initiation point from the initial bottom-hole flow pressure and the wellbore pressure change of each fracturing section;
- S5, bring the bottom-hole flow pressure and the bottom-hole pressure of each fracture initiation point into the formation flow fully implicit numerical model based on the embedded discrete fracture model (EDFM) in step S2, and iteratively cycle steps S2-S4 until the iterative variables converge, obtain the bottom-hole flow pressure and bottom-hole pressure of each fracture initiation point at this time step;
- S6, bring the bottom-hole flow pressure and the total gas production and liquid production at the time step into the pressure drop model considering gravity loss and frictional resistance loss in the inclined section and vertical well section, and obtain the wellbore pressure and wellhead casing pressure of each section through the iterative cycle section by section.
- S7, for each time step, repeat steps S2-S6 to obtain the full wellbore pressure under different production times and draw a pressure distribution curve. For specific results, see
FIG. 2 andFIG. 3 . - Referring to
FIG. 5 , a device for calculating the full wellbore pressure of a fracturing horizontal well, includes a processor, an acquisition module for acquiring initial data for calculating the full wellbore pressure of a fracturing horizontal well, an output module and a storage module for outputting calculation results; a program for calculating the full wellbore pressure of a fracturing horizontal well is stored on the storage module, which can be run on the processor, and the program for calculating the full wellbore pressure of a fracturing horizontal well is executed by the processor to implement the above method step. - A computer readable storage medium having stored therein program code executable by a processor, the computer readable storage medium including a plurality of instructions configured to cause the processor to perform the described fracturing horizontal well process. wellbore pressure calculation method.
- The above is not intended to limit the present invention in any form. Although the present invention has been disclosed through the above-mentioned embodiments, it is not intended to limit the present invention. Any person skilled in the art, within the scope of the technical solution of the present invention, When the technical contents disclosed above can be used to make some changes or modifications to equivalent embodiments with equivalent changes, any simple modifications or equivalents to the above embodiments according to the technical essence of the present invention do not depart from the content of the technical solution of the present invention. Changes and modifications still fall within the scope of the technical solutions of the present invention.
Claims (5)
1. A method for calculating a full wellbore pressure of a fracturing horizontal well, wherein the method comprises the following steps:
S1, collect the wellbore structural parameters, fluid physical properties, reservoir parameters and historical production data at different time steps of the gas well, and establish a formation flow fully implicit numerical model based on an embedded discrete fracture model;
S2, calculate the formation flow fully implicit numerical model based on the embedded discrete fracture model under an inner boundary condition of a constant gas production rate at a certain time step, to obtain a gas production, liquid production, grid block pressure and initial bottom-hole flow pressure of each fracture;
S3, establish a pressure drop model considering the gravity loss, frictional resistance loss and fracture convergence loss by combining the gas production, liquid production, and the pressure of the grid block where the fracture initiation point is located, to obtain a wellbore pressure variation of each fracturing section of the horizontal well;
S4, calculate the bottom-hole pressure of each fracture initiation point from the initial bottom-hole flow pressure and the wellbore pressure change of each fracturing section;
S5, perform the formation flow fully implicit numerical model based on the embedded discrete fracture model in step S2 by combining the bottom-hole flow pressure and the bottom-hole pressure of each fracture initiation point, and perform iterative cycles of the steps S2-S4 until the variables converge, to obtain the bottom-hole flow pressure and bottom-hole pressure at each fracture initiation point under this time step;
S6, bring the bottom-hole flow pressure and the total gas production and liquid production at the time step into the pressure drop model considering gravity loss and frictional resistance loss in the inclined section and vertical well section, to obtain wellbore pressure and wellhead casing pressure of each section through the iterative cycle at each section;
S7, for each time step, repeat steps S2-S6 to obtain the full wellbore pressure at different production times.
2. The full wellbore pressure calculation method of fracturing horizontal well according to claim 1 , wherein establishment of the formation flow fully implicit numerical model based on the embedded discrete fracture model described in the step S1 comprises the following steps:
S11, establish a formation flow mathematical model based on the embedded discrete fracture model, in which the dual medium model is used to describe the flow exchange between the matrix and micro-fractures, and the embedded discrete fracture model is used to describe the flow exchange between the matrix and artificial fractures, based on mass conservation equation and Darcy's law to establish seepage differential equation, the formation flow mathematical model of the embedded discrete fracture model is as follows:
matrix system seepage equation:
micro-fracture system seepage equation:
artificial fracture system seepage equation:
supplementary equation is as follows:
wherein km and kf represent the permeability of matrix, and micro-fracture, respectively, the unit is D; l=g or w, where g represents gas phase and w represents liquid phase; krlm and krlf represent the relative permeability of matrix, and micro-fracture, the unit is D; μlm, μlf, μlF represent the fluid viscosity of matrix, micro-fracture and artificial fracture system, the unit is Pa·s; Blm, Blf, BlF represent the fluid volume coefficients of the matrix, micro-fracture, and artificial fracture systems, respectively, dimensionless; φm, φf, φF represent the porosity of the matrix, micro-fractures, and artificial fracture systems, respectively, dimensionless; Slm, Slf, SlF represent saturation, dimensionless; t represents production time, the unit is day; plm, plf, plF represent the two-phase pressure of matrix, micro-fracture and artificial fracture system, the unit is MPa; ql,fm represents the flow exchange between matrix and micro-fractures based on the dual medium model, the unit is 1/day, wherein l=g,w; α represents the shape factor between matrix and micro-fractures, the unit is m−2; ql,F represents the flow exchange between the grids after an artificial fracture is divided into multiple grids by the matrix grid, the unit is 1/day; ql,FF represents the flow exchange between different artificial fracture grids, the unit is 1/day; ql,Ff represents the flow exchange between the artificial fracture grid and the matrix grid, the unit is 1/day; ql,well represents the flow exchange between the artificial fracture grid and the connected well grid, the unit is 1/day; GF represents geometric parameters for the connection pairs of grids in the same fracture, dimensionless; GFF represents the geometric parameter of the connection pairs of intersected fracture grids between different fractures, dimensionless; GFf is the geometric parameter of the connection pair between the matrix grid and intersected fracture grid, dimensionless; Gwell represents geometric parameters for connected pairs between fracture grid and intersected well grid, dimensionless;
S12, use a three-dimensional cubic grid to perform a space grid discretization for the reservoir, and based on the embedded discrete fracture model, the fracturing fractures are embedded in the three-dimensional cubic grid in a two-dimensional plane;
S13, using a finite difference method to establish the formation flow fully implicit numerical model based on the embedded discrete fracture model, comprising steps:
for the seepage equation (1) of the matrix system, the flow exchange relationship between the matrix grids on the left is written in the flow form:
in the formula: ql,mm is the flow exchange between connected matrix grids, the unit is 1/day;
the finite difference method is used to derive and expand the flow exchange term on the left of the above seepage equation, if there is a flow exchange between i and j grids, the fully implicit iterative form of the flow exchange between the two grids is:
q lij,mm n+1 =a lmi δp lmi +a lmj δp lmj +b lmi δS lmwi +b lmj δS lmwj +d lmij (6)
q lij,mm n+1 =a lmi δp lmi +a lmj δp lmj +b lmi δS lmwi +b lmj δS lmwj +d lmij (6)
in the formula: n+1 representing the next time step; δplmi, δplmj representing the iterative pressure variables in the i-th and j-th grids of the matrix system; almi, almj representing coefficients of iterative pressure variables in the i-th and j-th grids of the matrix system; δSlmi, δSlmj representing the iterative saturation variable in the i-th and j-th grids of the matrix system; blmi, blmj representing coefficients the iterative saturation variables in the i-th and j-th grids of the matrix system; dlmij representing the known quantity of flow exchange between the i-th and j-th grids of the matrix system;
the mass accumulation term on the right of the above seepage equation also derives and expands the unknowns, the coefficients of the unknowns are added to the coefficient matrix, and the known quantities are brought into the known quantity matrix to obtain:
mm×δm=d m (7)
mm×δm=d m (7)
in the formula: mm representing the coefficient matrix in the matrix flow equation; δm representing the unknown quantity matrix in the matrix flow equation; dm representing the matrix of known quantities in the matrix flow equation;
in the same way obtaining the coefficient matrix ff, FF, ww in the micro-fracture, artificial fracture and wellbore flow equation; the flow exchange coefficient matrix mf between the matrix and the micro-fracture system; the flow exchange coefficient matrix fF between the micro-fracture and the artificial fracture; flow exchange coefficient matrix Fw between artificial fracture and wellbore; known quantity matrix df, dF, dw in micro-fracture, artificial fracture, and wellbore flow equation; unknown matrix δf, δF, δw in micro-fractures, artificial fractures, and wellbore flow equations;
establish a fully implicit solver matrix based on the finite difference method:
equation (8) is solved to obtain the gas production and liquid production of each fracture, the grid block pressure where each fracture is located, and the initial bottom-hole flow pressure.
3. The full wellbore pressure calculation method of fracturing horizontal well according to claim 2 , it is characterized in that, the step S3 comprises the following steps:
S31, calculate apparent flow rates of gas and liquid at the initiation point of each fracture from the gas production and liquid production of each fracture, and combine the apparent flow rates of gas and liquid into the Mukherjee-Brill liquid holdup formula to calculate a liquid holdup at each fracture initiation point:
v sg =q sc B g /A
v sl =q w B w /A (9)
v sg =q sc B g /A
v sl =q w B w /A (9)
wherein vsl represents the apparent flow rate of the liquid phase, the unit is m/s; vsg representing the apparent flow rate of the gas phase, the unit is m/s; qsc represents the fracture gas production, the unit m3/s; qw represents the fracture fluid production, the unit is m3/s; Bg represents the volume coefficient of the gas phase, dimensionless; Bw represents the liquid phase volume coefficient, dimensionless; A represents the cross-sectional area of the wellbore, the unit is m2;
in the formula: Hl represents the liquid holdup, dimensionless; c1=−0.380113; c2=−0.129875; c3=−0.119788; c4=2.343227; c5=0.475686; c6=0.288657; θ represents the angle between the wellbore and the horizontal direction, the unit is °;
the Nlv, Ngv, Nl in the above formula are as follows:
in the formula: ρl represents the liquid density, the unit is kg/m3; δ represents the gas-liquid interfacial tension, the unit is N/m; μl represents the liquid viscosity, the unit is Pa·s;
S32, calculate the density ρm of the mixture, friction coefficient fm from the liquid holdup at each fracture initiation point and the grid block pressure, and the pressure gradient equation considering gravity loss, frictional resistance loss and fracture convergence loss is established as follows:
wherein dp/dL represents the pressure gradient, the unit is MPa/m; dL represents the length of the horizontal section between adjacent fractures, the unit is m; ρm represents the mixed density, the unit is kg/m3; G represents the acceleration of gravity, the unit is m/s2; θ represents the angle between the wellbore and the horizontal direction; fm represents the two-phase friction coefficient, dimensionless; vm represents the mixing speed, the unit is m/s; D represents the diameter of the pipe, the unit is m;
S33, combine the pressure gradient equation (12) and the distance ΔLi from the i-th fracture initiation point to the bottom-hole, to obtain the wellbore pressure changes in each fracturing section of the horizontal well as follows:
in the formula: NF represents the total number of fractures; ΔLi represents the distance from the i-th fracture initiation point to the bottom-hole of the well, the unit is m; Δpwi represents the horizontal wellbore pressure change from the bottom-hole to the i-th fracture, the unit is MPa.
4. A device for calculating the full wellbore pressure of a fracturing horizontal well, it is characterized in that, comprising:
processor, and
acquisition module configured to acquire the initial data for calculating the full wellbore pressure of the fracturing horizontal well;
storage module having stored thereon a program operable on the processor for calculating the full wellbore pressure of a fracturing horizontal well, the program for calculating the full wellbore pressure of a fracturing horizontal well being processed by the processor implements the steps of the method according to any one of claims 1 -3 ;
output module configured to output calculation results.
5. A non-transitory tangible computer readable storage medium, in which program code executable by a processor is stored, it is characterized in that, the computer readable storage medium comprises a plurality of instructions configured to cause the processor to perform the method for calculating the full wellbore pressure of a fractured horizontal well of any one of claims 1 -3 .
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210160050.0 | 2022-02-22 | ||
CN202210160050.0A CN114580100B (en) | 2022-02-22 | 2022-02-22 | Method and device for calculating full wellbore pressure of fractured horizontal well and computer readable storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
US11733425B1 US11733425B1 (en) | 2023-08-22 |
US20230266498A1 true US20230266498A1 (en) | 2023-08-24 |
Family
ID=81774718
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US17/890,248 Active US11733425B1 (en) | 2022-02-22 | 2022-08-17 | Full wellbore pressure calculation method, device and computer readable storage medium for fracturing horizontal well |
Country Status (2)
Country | Link |
---|---|
US (1) | US11733425B1 (en) |
CN (1) | CN114580100B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117556744B (en) * | 2024-01-11 | 2024-03-26 | 成都理工大学 | Method for predicting productivity of water-gas well produced by abnormal high-pressure low-permeability gas reservoir |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105840187B (en) * | 2016-06-03 | 2018-12-25 | 陕西延长石油(集团)有限责任公司研究院 | Compactness reservoir horizontal well staged fracturing Productivity |
CN107622328B (en) * | 2017-09-22 | 2019-11-05 | 西南石油大学 | A kind of more cluster pressure break horizontal well production prediction methods of shale gas reservoir segmentation |
EP3768939A4 (en) * | 2018-03-21 | 2021-12-29 | Resfrac Corporation | Systems and methods for hydraulic fracture and reservoir simulation |
CN109025942B (en) * | 2018-08-09 | 2021-08-17 | 西南石油大学 | Yield calculation method for irregular and multi-crack fracturing of tight gas reservoir inclined shaft |
CN110056336B (en) * | 2019-05-31 | 2020-01-07 | 西南石油大学 | Automatic diagnosis method for shale gas fracture network fracturing construction pressure curve |
CN110206522B (en) * | 2019-06-10 | 2019-12-20 | 西南石油大学 | Shale gas reservoir fracturing horizontal well fracturing fluid flowback simulation method |
CN112434426B (en) * | 2020-11-26 | 2022-09-27 | 北京科技大学 | Shale gas multistage fracturing horizontal well step gradient pressure drop development method and device |
CN112253103B (en) * | 2020-11-27 | 2021-08-31 | 西南石油大学 | Shale gas reservoir fracturing horizontal well yield prediction method based on random fracture model |
CN113076676B (en) * | 2021-01-19 | 2022-08-02 | 西南石油大学 | Unconventional oil and gas reservoir horizontal well fracture network expansion and production dynamic coupling method |
-
2022
- 2022-02-22 CN CN202210160050.0A patent/CN114580100B/en active Active
- 2022-08-17 US US17/890,248 patent/US11733425B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN114580100B (en) | 2022-09-16 |
US11733425B1 (en) | 2023-08-22 |
CN114580100A (en) | 2022-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110334431B (en) | Single-well control reserve calculation and residual gas analysis method for low-permeability tight gas reservoir | |
CN113826099B (en) | Flow simulation and transient well analysis method based on generalized pipe flow seepage coupling | |
US10908317B2 (en) | Method for experimentally determining conductivity distribution of acid-fractured fracture of carbonate rock oil and gas reservoir | |
Sun et al. | Investigating the effect of improved fracture conductivity on production performance of hydraulically fractured wells: Field-case studies and numerical simulations | |
Warpinski | Measurement of width and pressure in a propagating hydraulic fracture | |
Wang et al. | Estimation of the REV size and equivalent permeability coefficient of fractured rock masses with an emphasis on comparing the radial and unidirectional flow configurations | |
US20230229830A1 (en) | Method for coupling hydraulic fracture network extension and production performance of horizontal well in unconventional oil and gas reservoir | |
CN109522634B (en) | Numerical analysis method for compact gas multistage volume fracturing horizontal well | |
CN107545113A (en) | The complicated seam net of untraditional reservoir hydraulic fracturing forms process analogy method | |
WO2016192077A1 (en) | Method for establishing and solving numerical well-testing model of horizontal well for tight gas hydraulic fracturing | |
CN108319756A (en) | A kind of compact reservoir volume fracturing seam net extended simulation and characterizing method | |
Perera et al. | Numerical simulation of gas flow through porous sandstone and its experimental validation | |
Raghavan | Well-test analysis for multiphase flow | |
Yu et al. | The anisotropic seepage analysis of water-sealed underground oil storage caverns | |
CN104895550A (en) | Tight gas fracturing horizontal well numerical value well testing model building and solving method | |
US9418184B2 (en) | Determining flow through a fracture junction in a complex fracture network | |
Ahn et al. | Modeling of hydraulic fracture propagation in shale gas reservoirs: a three-dimensional, two-phase model | |
Zhao et al. | Simulation of a multistage fractured horizontal well in a tight oil reservoir using an embedded discrete fracture model | |
US20230266498A1 (en) | Full wellbore pressure calculation method, device and computer readable storage medium for fracturing horizontal well | |
Behmanesh | Rate-transient analysis of tight gas condensate and black oil wells exhibiting two-phase flow | |
Xu et al. | Production performance of horizontal wells with dendritic-like hydraulic fractures in tight gas reservoirs | |
Akbarnejad-Nesheli et al. | Relating fracture network characteristics to shale gas reserve estimation | |
CN105931125A (en) | Method for predicting yield of compact oil staged multi-cluster volume fracturing horizontal well | |
Liu et al. | Manufacturing method of large-scale fractured porous media for experimental reservoir simulation | |
Lu et al. | Influence factors of unpropped fracture conductivity of shale |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
AS | Assignment |
Owner name: SOUTHWEST PETROLEUM UNIVERSITY, CHINA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ZHANG, RUIHAN;ZHAO, YULONG;ZHENG, JIAN;AND OTHERS;REEL/FRAME:060838/0383 Effective date: 20220803 |
|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO SMALL (ORIGINAL EVENT CODE: SMAL); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |