CN107133373A - A kind of coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network - Google Patents
A kind of coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network Download PDFInfo
- Publication number
- CN107133373A CN107133373A CN201610111437.1A CN201610111437A CN107133373A CN 107133373 A CN107133373 A CN 107133373A CN 201610111437 A CN201610111437 A CN 201610111437A CN 107133373 A CN107133373 A CN 107133373A
- Authority
- CN
- China
- Prior art keywords
- pit shaft
- msub
- gas reservoir
- pipeline network
- mrow
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000010168 coupling process Methods 0.000 claims abstract description 102
- 230000008878 coupling Effects 0.000 claims abstract description 98
- 238000005859 coupling reaction Methods 0.000 claims abstract description 98
- 239000012071 phase Substances 0.000 claims description 36
- 239000008346 aqueous phase Substances 0.000 claims description 32
- 239000000463 material Substances 0.000 claims description 21
- 238000012797 qualification Methods 0.000 claims description 21
- 239000012530 fluid Substances 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 9
- 230000035699 permeability Effects 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000009795 derivation Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 239000007788 liquid Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification 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/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network, belong to numerical simulation for oil-gas reservoir technical field, solve that existing shale gas well method for numerical simulation solution procedure is unstable and the larger technical problem of error.This method includes：Obtain the parameters for numerical simulation of shale gas reservoir, pit shaft and surface pipeline network；Set up the numerical simulation fully implicit solution coupling model of shale gas reservoir, pit shaft and surface pipeline network；The fully implicit solution coupling model is carried out according to the parameters for numerical simulation to solve the numerical simulation result for obtaining shale gas reservoir, pit shaft and surface pipeline network.
Description
Technical field
The present invention relates to numerical simulation for oil-gas reservoir technical field, specifically, it is related to a kind of shale gas reservoir, pit shaft
And the coupled simulation method of surface pipeline network.
Background technology
In current shale gas Resource Development Process both at home and abroad, numerical simulation means are widely used.Numerical simulation
Application on shale gas well can not only carry out history matching to the conventional creation data of shale gas well, obtain shale
Each regional petrophysics parameter of reservoir and flow conductivity, can also calculate the physical parameter of fluid in shale gas well reservoir,
And the dynamic physical property and long-term production capacity to reservoir are showed and made prediction.Because its applicability is wide, reservoir Dynamic is retouched
The features such as specific, prediction is accurate is stated, numerical simulation means are on the Dynamic profiling and capability forecasting of shale gas well
It is widely used.
However, at present in the numerical simulation of shale gas well, the boundary condition of gas reservoir model is often assumed
To fix flowing bottomhole pressure (FBHP) or sandface flow rate, but actual conditions that this and shale gas are developed differ greatly.True
In situation, because the real-time dynamic change of pit shaft and fluid in surface pipeline network is more violent, cause each parameter in shaft bottom
(flowing bottomhole pressure (FBHP), sandface flow rate) can have a greater change with the time.Therefore to mould in conventional oil reservoir/gas reservoir simulation
The hypothesis of type boundary condition can cause the error of analog result.Although existing half implicit and complete explicit coupling process
The problem of can solving to assume not conforming to the actual conditions in convenient value simulation, but still suffer from convergence difficulties (system equation
It is unstable), the problems such as error is larger.
Therefore, need one kind badly and can stablize, it is accurate gas reservoir model and pit shaft/pipe net leakage rate are carried out and meanwhile calculate and
The analogy method of solution, it is accurate, stability forecast to hydrodynamic and gas reservoir production capacity in shale gas exploitation to meet
Demand.
The content of the invention
It is an object of the invention to provide a kind of coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network, with
Solve that existing shale gas well method for numerical simulation solution procedure is unstable and the larger technical problem of error.
The present invention provides a kind of coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network, and this method includes：
Obtain the parameters for numerical simulation of shale gas reservoir, pit shaft and surface pipeline network；
Set up the numerical simulation fully implicit solution coupling model of shale gas reservoir, pit shaft and surface pipeline network；
The fully implicit solution coupling model solve according to the parameters for numerical simulation and obtains shale gas reservoir, pit shaft
And the numerical simulation result of surface pipeline network.
It is described set up fully implicit solution coupling model the step of include：
Set up shale gas reservoir governing equation, pit shaft governing equation, surface pipeline network governing equation；
Set up the Coupling point governing equation of Coupling point governing equation, pit shaft and the surface pipeline network of shale gas reservoir and pit shaft；
Define the boundary condition of the fully implicit solution coupling model；
It is defined on the qualifications of more new variables when being solved to the fully implicit solution coupling model.
It is described set up shale gas reservoir governing equation the step of include：
Build shale gas reservoir flowing material equilibrium equation：
Build shale gas reservoir absorbing boundary equation：
Wherein, RrwFor the material balance differential equation of aqueous phase in gas reservoir, RrgFor the material balance of gas phase in gas reservoir
The differential equation, ρw、ρgThe respectively gentle phase density of gas reservoir aqueous phase densities, krw、krgRespectively aqueous phase and gas phase
Relative permeability, k is reservoir permeability, μw、μgThe respectively gentle phase viscosity of gas reservoir aqueous viscosity, pw、pg
Respectively Gas Reservoir Water phase pressure and gaseous pressure,For divergence of the vertical direction to all directions, φ is porosity,
Sw、SgThe respectively gentle phase saturation of gas reservoir water phase saturation,Respectively converge or source aqueous phase, gas phase
Flow, RrbFor gas reservoir border residual equation, WIpFor pit shaft index,For sandface flow rate,For pit shaft
The pressure of place grid future time step, pwfFor bottom pressure.
It is described set up pit shaft and surface pipeline network governing equation the step of include：
Build pit shaft/surface pipeline network flowing material equilibrium equation：
Rt/nw=qwi-qw(i-1)；
Rt/ng=qgi-qg(i-1)；
Build pit shaft/surface pipeline network fluid dynamic energy equilibrium equation：
Wherein, Rt/nwFor aqueous phase material balance residual equation, R in pit shaft/surface pipeline networkt/ngFor pit shaft/ground line
Gaseous substance balances residual equation, q in netwi、qw(i-1)Grid i and grid respectively in pit shaft/surface pipeline network
I-1 aqueous phase flow, qgi、qg(i-1)Grid i and grid i-1 aqueous phase flow respectively in pit shaft/surface pipeline network,
Rt/npFor the remainder equation of pit shaft/surface pipeline network kinetic energy balancing equation,Respectively pit shaft/surface pipeline network
Middle grid i and grid i-1 pressure, Δ pnjFor each section of the pressure loss of pit shaft/surface pipeline network, j is i and i-1
The sequence number of the interface location of two grids.
It is described set up Coupling point governing equation the step of include：
Build gas reservoir and pit shaft/pit shaft and the Coupling point equation of surface pipeline network：
Rcp=pn-pwf；
Wherein, RcwIt is that gas reservoir and pit shaft/pit shaft balance coupled wave equation, R with surface pipeline network Coupling point aqueous phase substancecg
It is that gas reservoir and pit shaft/pit shaft balance coupled wave equation, R with surface pipeline network Coupling point gaseous substancecpFor gas reservoir and pit shaft/
Pit shaft is coupled equation with surface pipeline network Coupling point kinetic energy balancing,For the aqueous phase flow of pit shaft/surface pipeline network,For shaft bottom aqueous phase flow, qngFor the gas phase flow rate of pit shaft/surface pipeline network,For shaft bottom gas phase flow rate, pn
For the bottom pressure of wellbore model, pwfFor the bottom pressure of gas reservoir model.
Include the step of the boundary condition of the definition fully implicit solution coupling model：
It is surface pipeline network separator or the pressure of end port to define the border of fully implicit solution coupling model；
The pressure of the surface pipeline network separator or end port is set, definite value is or on the time
Function.
It is described fully implicit solution coupling model is solved the step of include：
Bring the parameters for numerical simulation into the fully implicit solution coupling model；
By each equation in the coupling model respectively to each variable derivation in the coupling model, coupling is set up
The linear equation of model；
It is divided by using the linear equation with the governing equation, is solved according to Newton iteration method and obtain the coupling
The approximation of model root.
It is described fully implicit solution coupling model is solved the step of in also include：
The scope of the approximation of described is limited according to the qualifications of the more new variables, to ensure
State the convergence range of coupling model.
The scope of the approximation of described is limited in the qualifications of more new variables described in the basis
Step includes：
Judge whether the approximation of described exceeds the qualifications scope of the more new variables, if going beyond the scope,
Then the root of the coupling model takes the qualifications value of the more new variables.
The coupled simulation method of shale gas reservoir provided in an embodiment of the present invention, pit shaft and surface pipeline network is considered in page
In rock gas actual production process, the boundary condition of system is the port pressure of surface pipeline network, is setting up model system
During governing equation, the material balance of pit shaft/pipe network and momentum balance equation are added in the governing equation of gas reservoir,
It is achieved thereby that the coupling to gas reservoir model and pipe net leakage rate.Also, in view of including many in coupled system
Individual main body (gas reservoir, pit shaft, surface pipeline network), in the iteration mistake for the governing equation for calculating and solving coupled system
Cheng Zhong, accomplishes synchronous calculating to the governing equation of each main body and solves, ensure that each governing equation
Synchronous convergence, it is to avoid the vibration in convergence process.Meanwhile, gas reservoir governing equation is considered in solution procedure
With the difference of basic physical property in pipe flow control equation, in iteration more new variables, type range of variables is entered
Go limitation, ensure that the convergence range of equation.
This method by setting up the coupling model governing equation of shale gas gas reservoir, pit shaft, the human subject of surface pipeline network three,
And it is multiple ensure the convergent qualifications of this model equation, realize for shale gas reservoir, pit shaft and surface pipeline network
Numerical simulation, overcome in half implicit coupling process and be when being solved caused by the governing equation difference of each main body
The problem of equation of uniting vibrates brought convergence difficulties, also overcomes the problem of error is big in full explicit method, real
Now gas reservoir and pit shaft/surface pipeline network are modeled and solved by the model cootrol equation of complete set.
Other features and advantages of the present invention will be illustrated in the following description, also, partial from specification
In become apparent, or by implement the present invention and understand.The purpose of the present invention and other advantages can pass through
Specifically noted structure is realized and obtained in specification, claims and accompanying drawing.
Brief description of the drawings
, below will be to needed for embodiment description for the technical scheme in the clearer explanation embodiment of the present invention
The accompanying drawing wanted does simple introduction：
Fig. 1 is the schematic flow sheet of coupled simulation method provided in an embodiment of the present invention；
Fig. 2 is the schematic flow sheet provided in an embodiment of the present invention for setting up fully implicit solution coupling model；
Fig. 3 is the schematic flow sheet provided in an embodiment of the present invention solved to fully implicit solution coupling model；
Fig. 4 is the linear equation distribution map of fully implicit solution coupling model provided in an embodiment of the present invention；
Fig. 5 is fully implicit solution model provided in an embodiment of the present invention and half implicit model convergence rate comparison diagram.
Embodiment
Describe embodiments of the present invention in detail below with reference to drawings and Examples, whereby to the present invention how
Application technology means solve technical problem, and reach the implementation process of technique effect and can fully understand and real according to this
Apply.As long as it should be noted that conflict is not constituted, in each embodiment and each embodiment in the present invention
Each feature can be combined with each other, and the technical scheme formed is within protection scope of the present invention.
The coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network that the present invention is provided in implementing, such as Fig. 1 institutes
Show, including：Step 101, step 102 and step 103.Wherein, in a step 101 obtain shale gas reservoir,
The parameters for numerical simulation of pit shaft and surface pipeline network, such as gas reservoir model permeability, porosity, mineshaft diameter, pit shaft
Frictional resistance, geostatic pressure, surface temperature etc..
In a step 102, the numerical simulation fully implicit solution coupling model of shale gas reservoir, pit shaft and surface pipeline network is set up.
As shown in Fig. 2 including the step of coupling model is set up：Step 201 is to step 205.In step
In 201, shale gas reservoir governing equation is set up, in this step, the material balance differential side of gas reservoir aqueous phase is built
Journey, the material balance differential equation of gas reservoir gas phase and gas reservoir border residual equation.
Aqueous phase in shale gas reservoir is built and during the gas phase material balance differential equation, by by the equation of motion (darcy
Law) bring conservation of matter equation into and by obtaining last gas reservoir mould after fluid state equation simplification its coefficient
Type flowing material equilibrium equation：：
Wherein, RrwFor the material balance differential equation of aqueous phase in gas reservoir, RrgFor the material balance of gas phase in gas reservoir
The differential equation, ρw、ρgThe respectively gentle phase density of gas reservoir aqueous phase densities, krw、krgRespectively aqueous phase and gas phase
Relative permeability, k is reservoir permeability, μw、μgThe respectively gentle phase viscosity of gas reservoir aqueous viscosity, pw、pg
Respectively Gas Reservoir Water phase pressure and gaseous pressure,For divergence of the vertical direction to all directions, φ is porosity,
Sw、SgThe respectively gentle phase saturation of gas reservoir water phase saturation,Respectively converge or source aqueous phase, gas phase
Flow.
The absorbing boundary equation of gas reservoir model is the qualifications to flowing bottomhole pressure (FBHP) or sandface flow rate, and shale gas reservoir border is remained
Remaining equation directly carries out absorbing boundary equation foundation, such as constant pressure outer boundary equation (boundary bit by defining boundary types
Put pressure constant) or closed outer boundary equation (boundary position barometric gradient is 0) etc..In the embodiment of the present invention
In, the absorbing boundary equation for setting up gas reservoir model is：
R in formularbFor gas reservoir border residual equation (being based on Peaceman equations), WIpFor pit shaft index,For well
Underflow amount,The pressure of grid future time step, p where pit shaftwfFor bottom pressure.
Further, in a step 102, pit shaft and surface pipeline network governing equation, in this step, structure are set up
Material balance residual equation, the kinetic energy balancing residual equation of the gentle phase of aqueous phase in shaft building cylinder and surface pipeline network.
Wherein, pit shaft/surface pipeline network matter balance equation is：
Rt/nw=qwi-qw(i-1)
Rt/ng=qgi-qg(i-1)
In formula, Rt/nwFor aqueous phase material balance residual equation, R in pit shaft/surface pipeline networkt/ngFor pit shaft/ground line
Gaseous substance balances residual equation, q in netwi、qw(i-1)Grid i and grid respectively in pit shaft/surface pipeline network
I-1 aqueous phase flow, qgi、qg(i-1)Grid i and grid i-1 gas phase flow rate respectively in pit shaft/surface pipeline network.
Stream of the correlation formula to each section in one-dimensional pit shaft/pipe net leakage rate is being utilized during setting up kinetic energy balancing equation
Type is judged, each grid of pit shaft/surface pipeline network is may determine that according to the correlation formula of pit shaft/surface pipeline network
Liquid holdup, the parameter such as the density of fluid, viscosity in each grid is obtained according to liquid holdup, with obtained fluid parameter
Calculate the pressure loss Δ p of each grid in pit shaftnj, further obtain pit shaft/surface pipeline network kinetic energy balancing equation：
In formula：Rt/npFor the remainder equation of pit shaft/surface pipeline network kinetic energy balancing equation,Respectively well
Grid i and grid i-1 pressure in cylinder/surface pipeline network, Δ pnjFor each section of the pressure loss of pit shaft/surface pipeline network, j
For the sequence number of the interface location of two grids of i and i-1.
Further, in step 203, Coupling point governing equation is set up.Coupling point is shaft bottom Coupling point and well
Mouth Coupling point, shaft bottom Coupling point is the Coupling point of shale gas reservoir model and wellbore model, and well head Coupling point is pit shaft
The Coupling point of model and surface pipeline network model.
In this step, shaft bottom Coupling point and well head Coupling point exchange adjacent governing equation parameter (including pressure,
Each phase flow rate), by setting up Coupling point governing equation, realize each main body equation (gas of fully implicit solution coupling model
Hide governing equation, pit shaft governing equation, pipe network governing equation) simultaneous.So that the perimeter strip of gas reservoir model
Part (sandface flow rate of gas reservoir model is pressed with stream) (shaft bottom of wellbore model equal with the boundary condition of wellbore model
Flow is pressed with stream), make the boundary condition of wellbore model equal with the boundary condition of terrestrial network tube model.
In embodiments of the present invention, the Coupling point control of gas reservoir and pit shaft Coupling point, pit shaft and surface pipeline network Coupling point
Equation processed is identical, because model hypothesis wellbore fluids are the quality stream of stationary flow, i.e. pit shaft and surface pipeline network each point
Amount is identical.
Coupling point governing equation is：
Rcp=pn-pwf
R in formulacwIt is that gas reservoir and pit shaft/pit shaft balance coupled wave equation, R with surface pipeline network Coupling point aqueous phase substancecgFor
Gas reservoir balances coupled wave equation, R with pit shaft/pit shaft with surface pipeline network Coupling point gaseous substancecpFor gas reservoir and pit shaft/well
Cylinder is coupled equation with surface pipeline network Coupling point kinetic energy balancing,For the aqueous phase flow of pit shaft/surface pipeline network,
For shaft bottom aqueous phase flow, qngFor the gas phase flow rate of pit shaft/surface pipeline network,For shaft bottom gas phase flow rate, pnFor
The bottom pressure of wellbore model, pwfFor the bottom pressure of gas reservoir model.
Further, in step 204, the boundary condition of fully implicit solution coupling model is defined.In this step,
First, it is surface pipeline network separator or the pressure of end port to define the border of fully implicit solution coupling model.Then, it is right
The pressure of the surface pipeline network separator or end port is provided, becomes definite value or the function on the time.
Can be by ground side by the material balance residual equation of above-mentioned pit shaft and surface pipeline network, kinetic energy balancing residual equation
Boundary's pressure convert is shaft bottom terminal pressure.
Further, in step 205, it is defined on more new variables when being solved to the fully implicit solution coupling model
Qualifications.In this step, define including the pressure in gas reservoir model in grid, each section in conduit flow model
The scope of the parameters such as interior pressure, flow is provided, makes each parameter in fully implicit solution coupling model solution procedure
No more than real physical values scope.
As shown in figure 1, in step 103, fully implicit solution is carried out to the coupling model according to the analog parameter
Difference solves the numerical simulation result for obtaining shale gas reservoir, pit shaft and surface pipeline network.
As shown in figure 3, including the step of being solved to fully implicit solution coupling model：Step 301 is to step
305.In step 301, bring parameters for numerical simulation into fully implicit solution coupling model, linearize coupling model side
Journey, by each equation in coupling model respectively to each variable derivation in coupling model, sets up the line of coupling model
Property equation：
In formula：R is each equation in coupling model, and U is each variable in coupling model, and r is gas reservoir model,
T is wellbore model, and n is surface pipeline network model, RrIncluding Rrw、Rrg、Rrb, RtAnd RnInclude R respectivelyt/nw、
Rt/ng、Rt/np、Rcw、Rcg、Rcp。
In step 302, coupling model linear equation is solved.In this step, linear equation and gas are utilized
Tibetan/pit shaft/surface pipeline network governing equation is divided by, and the approximation for obtaining coupling model root is solved according to Newton iteration method.
In an iterative process, the system variable value of following iteration is the system variable value and linear equation root of current iteration
The sum of approximation.
Solution formula is：
In formula：δU*For the approximation of coupling model linear equation root, U*+1For the system variable of following iteration, R*
For each equation in coupling model, U*For the system variable of current iteration,For the linear equation of coupling model.
After the approximation of coupling model linear equation root is obtained, using the qualifications of more new variables to each ox
The approximation of root obtained by iteration of pausing is judged, it is judged as the reasonability of iteration more new variables, repeatedly
During generation more new variables, the scope of the approximation of root is limited according to the qualifications of more new variables, to ensure
The convergence range of coupling model.I.e. in step 303, judge whether the approximation of gained root exceeds more new variables
Qualifications scope, if the approximation of root exceed qualifications scope, perform step 304, coupled mode
Root (the system variable U of following iteration of type*+1) the qualifications value of more new variables is taken, if the approximation of root is not
Beyond the scope of qualifications, then step 305 is performed, root takes income value.
For example, as the new U of gained*+1Pressure in variable (that is, can not be negative value beyond actual physical scope
/ excessive), then force U*+1Pressure in variable is equal to the maximum of qualifications pressure value, i.e. model specification/most
Small value (can such as set the maximum/minimum value of pressure as：100MPa/0.1MPa).
The coupled simulation method provided using the present invention sets up the fully implicit solution coupled mode of a set of gas reservoir, pit shaft and pipe stream
Type, model is shale gas exploitation model, and including a bite pressure break horizontal well, concrete model parameter is as follows：
Vertical bore model meshes are set to 1 × 1 × 100, and hole diameter is 0.1m；Horizontal segment length 1000m, pressure break 15
Section, crack is double-vane symmetric fracture, and fracture half-length is 150m, and crack is infinite fluid diversion；Geostatic pressure is 2.8MPa.
Gas reservoir model meshes are set to 35 × 35 × 1, and permeability is 0.001mD, and porosity is 0.05, and gas phase is initially satisfied
It is 1 with degree.Gas reservoir top depth is 3000m, and thickness is 50m, and gas reservoir mean temperature is 105 DEG C, reset pressure
For 45MPa.Fig. 4 is the linear equation distribution map for the fully implicit solution coupling model set up, and the matrix lower right corner is pit shaft
/ surface pipeline network element.
Meanwhile, a set of identical model is set up using the implicit coupling process of tradition half, and to the convergence of two methods
Speed is compared, the iterations of each time step of two methods as shown in figure 5, in figure, abscissa
Time step, ordinate is iterations, fully implicit solution coupling model in convergence process, iterations is than half
Implicit model is more stable, and convergence rate is faster.
The coupled simulation method of shale gas reservoir provided in an embodiment of the present invention, pit shaft and surface pipeline network is considered in page
In rock gas actual production process, the boundary condition of system is the port pressure of surface pipeline network, is setting up model system
During governing equation, the material balance of pit shaft/pipe network and momentum balance equation are added in the governing equation of gas reservoir,
It is achieved thereby that the coupling to gas reservoir model and pipe net leakage rate.Also, in view of including many in coupled system
Individual main body (gas reservoir, pit shaft, surface pipeline network), in the iteration mistake for the governing equation for calculating and solving coupled system
Cheng Zhong, accomplishes synchronous calculating to the governing equation of each main body and solves, ensure that each governing equation
Synchronous convergence, it is to avoid the vibration in convergence process.Meanwhile, gas reservoir governing equation is considered in solution procedure
With the difference of basic physical property in pipe flow control equation, in iteration more new variables, type range of variables is entered
Go limitation, ensure that the convergence range of equation.
This method by setting up the coupling model governing equation of shale gas gas reservoir, pit shaft, the human subject of surface pipeline network three,
And it is multiple ensure the convergent qualifications of this model equation, realize for shale gas reservoir, pit shaft and surface pipeline network
Numerical simulation, overcome in half implicit coupling process and be when being solved caused by the governing equation difference of each main body
The problem of equation of uniting vibrates brought convergence difficulties, also overcomes the problem of error is big in full explicit method, real
Now gas reservoir and pit shaft/surface pipeline network are modeled and solved by the model cootrol equation of complete set.
While it is disclosed that embodiment as above, but described content is only to facilitate understand the present invention
And the embodiment used, it is not limited to the present invention.Technology people in any the technical field of the invention
Member, on the premise of spirit and scope disclosed in this invention are not departed from, can implementation formal and details
On make any modification and change, but the scope of patent protection of the present invention still must be with appended claims institute
The scope defined is defined.
Claims (9)
1. a kind of coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network, it is characterised in that including：
Obtain the parameters for numerical simulation of shale gas reservoir, pit shaft and surface pipeline network；
Set up the numerical simulation fully implicit solution coupling model of shale gas reservoir, pit shaft and surface pipeline network；
The fully implicit solution coupling model solve according to the parameters for numerical simulation and obtains shale gas reservoir, pit shaft
And the numerical simulation result of surface pipeline network.
2. the coupled simulation method of shale gas reservoir according to claim 1, pit shaft and surface pipeline network, its
Be characterised by, it is described set up fully implicit solution coupling model the step of include：
Set up shale gas reservoir governing equation, pit shaft governing equation, surface pipeline network governing equation；
Set up the Coupling point governing equation of Coupling point governing equation, pit shaft and the surface pipeline network of shale gas reservoir and pit shaft；
Define the boundary condition of the fully implicit solution coupling model；
It is defined on the qualifications of more new variables when being solved to the fully implicit solution coupling model.
3. the coupled simulation method of shale gas reservoir according to claim 2, pit shaft and surface pipeline network, its
Be characterised by, it is described set up shale gas reservoir governing equation the step of include：
Build shale gas reservoir flowing material equilibrium equation：
Build shale gas reservoir absorbing boundary equation：
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>q</mi>
<mi>w</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>WI</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<msubsup>
<mi>p</mi>
<mi>p</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mi>w</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, RrwFor the material balance differential equation of aqueous phase in gas reservoir, RrgFor the material balance of gas phase in gas reservoir
The differential equation, ρw、ρgThe respectively gentle phase density of gas reservoir aqueous phase densities, krw、krgRespectively aqueous phase and gas phase
Relative permeability, k is reservoir permeability, μw、μgThe respectively gentle phase viscosity of gas reservoir aqueous viscosity, pw、pg
Respectively Gas Reservoir Water phase pressure and gaseous pressure,For divergence of the vertical direction to all directions, φ is porosity,
Sw、SgThe respectively gentle phase saturation of gas reservoir water phase saturation,Respectively converge or source aqueous phase, gas phase
Flow, RrbFor gas reservoir border residual equation, WIpFor pit shaft index,For sandface flow rate,For pit shaft
The pressure of place grid future time step, pwfFor bottom pressure.
4. the coupled simulation method of shale gas reservoir according to claim 3, pit shaft and surface pipeline network, its
Be characterised by, it is described set up pit shaft and surface pipeline network governing equation the step of include：
Build pit shaft/surface pipeline network flowing material equilibrium equation：
Rt/nw=qwi-qw(i-1)；
Rt/ng=qgi-qg(i-1)；
Build pit shaft/surface pipeline network fluid dynamic energy equilibrium equation：
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>t</mi>
<mo>/</mo>
<mi>n</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>p</mi>
<msub>
<mi>n</mi>
<mi>i</mi>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<msub>
<mi>n</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;p</mi>
<mrow>
<mi>n</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
Wherein, Rt/nwFor aqueous phase material balance residual equation, R in pit shaft/surface pipeline networkt/ngFor pit shaft/ground line
Gaseous substance balances residual equation, q in netwi、qw(i-1)Grid i and grid respectively in pit shaft/surface pipeline network
I-1 aqueous phase flow, qgi、qg(i-1)Grid i and grid i-1 aqueous phase flow respectively in pit shaft/surface pipeline network,
Rt/npFor the remainder equation of pit shaft/surface pipeline network kinetic energy balancing equation,Respectively pit shaft/surface pipeline network
Middle grid i and grid i-1 pressure, Δ pnjFor each section of the pressure loss of pit shaft/surface pipeline network, j is i and i-1
The sequence number of the interface location of two grids.
5. the coupled mode of the shale gas reservoir, pit shaft and surface pipeline network according to any one of claim 2 to 4
Plan method, it is characterised in that it is described set up Coupling point governing equation the step of include：
Build gas reservoir and pit shaft/pit shaft and the Coupling point equation of surface pipeline network：
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>c</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>q</mi>
<msub>
<mi>n</mi>
<mi>w</mi>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mrow>
<msub>
<mi>wf</mi>
<mi>w</mi>
</msub>
</mrow>
</msub>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>c</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>n</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mrow>
<msub>
<mi>wf</mi>
<mi>g</mi>
</msub>
</mrow>
</msub>
<mo>;</mo>
</mrow>
Rcp=pn-pwf；
Wherein, RcwIt is that gas reservoir and pit shaft/pit shaft balance coupled wave equation, R with surface pipeline network Coupling point aqueous phase substancecg
It is that gas reservoir and pit shaft/pit shaft balance coupled wave equation, R with surface pipeline network Coupling point gaseous substancecpFor gas reservoir and pit shaft/
Pit shaft is coupled equation with surface pipeline network Coupling point kinetic energy balancing,For the aqueous phase flow of pit shaft/surface pipeline network,For shaft bottom aqueous phase flow, qngFor the gas phase flow rate of pit shaft/surface pipeline network,For shaft bottom gas phase flow rate, pn
For the bottom pressure of wellbore model, pwfFor the bottom pressure of gas reservoir model.
6. the coupled simulation method of shale gas reservoir according to claim 5, pit shaft and surface pipeline network, its
It is characterised by, includes the step of the boundary condition of the definition fully implicit solution coupling model：
It is surface pipeline network separator or the pressure of end port to define the border of fully implicit solution coupling model；
The pressure of the surface pipeline network separator or end port is set, definite value is or on the time
Function.
7. the coupled simulation method of shale gas reservoir according to claim 6, pit shaft and surface pipeline network, its
Be characterised by, it is described fully implicit solution coupling model is solved the step of include：
Bring the parameters for numerical simulation into the fully implicit solution coupling model；
By each equation in the coupling model respectively to each variable derivation in the coupling model, coupling is set up
The linear equation of model；
It is divided by using the linear equation with the governing equation, is solved according to Newton iteration method and obtain the coupling
The approximation of model root.
8. the coupled simulation method of shale gas reservoir according to claim 7, pit shaft and surface pipeline network, its
Be characterised by, it is described fully implicit solution coupling model is solved the step of in also include：
The scope of the approximation of described is limited according to the qualifications of the more new variables, to ensure
State the convergence range of coupling model.
9. the coupled simulation method of shale gas reservoir according to claim 8, pit shaft and surface pipeline network, its
It is characterised by, the scope of the approximation of described is limited in the qualifications of more new variables described in the basis
Processed the step of, includes：
Judge whether the approximation of described exceeds the qualifications scope of the more new variables, if going beyond the scope,
Then the root of the coupling model takes the qualifications value of the more new variables.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610111437.1A CN107133373B (en) | 2016-02-29 | 2016-02-29 | Coupling simulation method for shale gas reservoir, shaft and ground pipe network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610111437.1A CN107133373B (en) | 2016-02-29 | 2016-02-29 | Coupling simulation method for shale gas reservoir, shaft and ground pipe network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107133373A true CN107133373A (en) | 2017-09-05 |
| CN107133373B CN107133373B (en) | 2020-05-15 |
Family
ID=59720944
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610111437.1A Active CN107133373B (en) | 2016-02-29 | 2016-02-29 | Coupling simulation method for shale gas reservoir, shaft and ground pipe network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107133373B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111340298A (en) * | 2020-02-28 | 2020-06-26 | 西安石油大学 | Shale reservoir induced heterogeneity dynamic evolution prediction method |
| CN115860266A (en) * | 2023-02-16 | 2023-03-28 | 东北石油大学三亚海洋油气研究院 | Shale gas/coal bed gas well productivity evaluation method and system and electronic equipment |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7516056B2 (en) * | 2005-04-26 | 2009-04-07 | Schlumberger Technology Corporation | Apparatus, method and system for improved reservoir simulation using a multiplicative overlapping Schwarz preconditioning for adaptive implicit linear systems |
| CN103246820A (en) * | 2013-05-21 | 2013-08-14 | 中国石油大学(北京) | Oil and gas reservoir numerical simulation calculation method |
| CN104239602A (en) * | 2014-07-15 | 2014-12-24 | 浙江科技学院 | Shale gas reservoir characteristic simulation method |
| CN104533370A (en) * | 2014-11-06 | 2015-04-22 | 中国石油大学(北京) | Oil deposit, crack and shaft fully-coupled simulating method of fractured horizontal well |
| CN104881522A (en) * | 2015-05-05 | 2015-09-02 | 中国石油大学(北京) | Fluid potential field based remaining oil representation method and device |
| CN104989351A (en) * | 2015-06-01 | 2015-10-21 | 四川大学 | Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process |
-
2016
- 2016-02-29 CN CN201610111437.1A patent/CN107133373B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7516056B2 (en) * | 2005-04-26 | 2009-04-07 | Schlumberger Technology Corporation | Apparatus, method and system for improved reservoir simulation using a multiplicative overlapping Schwarz preconditioning for adaptive implicit linear systems |
| CN103246820A (en) * | 2013-05-21 | 2013-08-14 | 中国石油大学(北京) | Oil and gas reservoir numerical simulation calculation method |
| CN104239602A (en) * | 2014-07-15 | 2014-12-24 | 浙江科技学院 | Shale gas reservoir characteristic simulation method |
| CN104533370A (en) * | 2014-11-06 | 2015-04-22 | 中国石油大学(北京) | Oil deposit, crack and shaft fully-coupled simulating method of fractured horizontal well |
| CN104881522A (en) * | 2015-05-05 | 2015-09-02 | 中国石油大学(北京) | Fluid potential field based remaining oil representation method and device |
| CN104989351A (en) * | 2015-06-01 | 2015-10-21 | 四川大学 | Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111340298A (en) * | 2020-02-28 | 2020-06-26 | 西安石油大学 | Shale reservoir induced heterogeneity dynamic evolution prediction method |
| CN111340298B (en) * | 2020-02-28 | 2021-11-02 | 西安石油大学 | Shale reservoir induced heterogeneity dynamic evolution prediction method |
| CN115860266A (en) * | 2023-02-16 | 2023-03-28 | 东北石油大学三亚海洋油气研究院 | Shale gas/coal bed gas well productivity evaluation method and system and electronic equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107133373B (en) | 2020-05-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2020224539A1 (en) | Flow simulation and transient well analysis method based on generalized pipe flow seepage coupling | |
| CN111581786B (en) | Well test interpretation method of well test interpretation model for analyzing fracture-cave series-mode double-hole composite reservoir | |
| CN104504230B (en) | Estimation method for recovery ratio and limit drainage radius of low-permeability gas well | |
| CN108729908B (en) | Dense oil flow simulation and permeability prediction method based on pore network model | |
| An et al. | Influence of pore structure parameters on flow characteristics based on a digital rock and the pore network model | |
| CN106599449B (en) | Well test interpretation method for karst cave volume calculation | |
| CN107045671A (en) | Water-producing gas well hydrops Risk Forecast Method | |
| CN111236931B (en) | Method and system for generating unsteady well testing chart of gas reservoir highly-deviated well | |
| WO2016192077A1 (en) | Method for establishing and solving numerical well-testing model of horizontal well for tight gas hydraulic fracturing | |
| CN105574320B (en) | The evaluation method of the effective percolation ability of low permeability sandstone reservoirs | |
| CN106761733A (en) | A kind of horizontal wells in heavy oil reservoir steam soak initial productivity Forecasting Methodology | |
| CN104847314B (en) | HTHP oil gas straight well single-phase flow perforation completion parameter optimization method | |
| CN111428426A (en) | Shale gas multiphase flow simulation method based on lattice Boltzmann | |
| CN113836695B (en) | Oil reservoir numerical simulation method based on gridless connecting element | |
| CN107145671A (en) | A kind of numerical reservoir simulation method and system | |
| CN104989385B (en) | The HTHP oil gas straight well perforating parameter optimization method calculated based on skin factor | |
| CN106951641B (en) | A kind of method and system of fracture-pore reservoir numerical simulation | |
| CN107944620A (en) | A kind of Non-linear of individual well steady state productivity | |
| CN107133373A (en) | A kind of coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network | |
| Śliwa et al. | Theoretical model of borehole heat exchanger | |
| CN112163379B (en) | Unsteady two-phase relative permeability calculation method based on pore network model | |
| CN113919247A (en) | Complex fracture network oil reservoir flow simulation method | |
| Skorpa et al. | Leakages through radial cracks in cement sheaths: Effect of geometry, viscosity, and aperture | |
| CN113255247A (en) | Multi-scale numerical simulation method for oil reservoir in high water-cut stage | |
| CN109933951B (en) | Method for establishing multi-scale and multi-mechanism coupling seepage model of tight oil reservoir volume fracturing horizontal well |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |