CN110080745A - Separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS and device - Google Patents
Separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS and device Download PDFInfo
- Publication number
- CN110080745A CN110080745A CN201910408401.3A CN201910408401A CN110080745A CN 110080745 A CN110080745 A CN 110080745A CN 201910408401 A CN201910408401 A CN 201910408401A CN 110080745 A CN110080745 A CN 110080745A
- Authority
- CN
- China
- Prior art keywords
- model
- straight well
- separate stratum
- stratum fracfturing
- fracture
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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 DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Abstract
The present invention provides a kind of separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS and devices, this method comprises: establishing the reservoir model of separate stratum fracfturing straight well according to the oil deposit parameter of separate stratum fracfturing straight well；According to the fracture parameters of separate stratum fracfturing straight well, the fisstured flow model of separate stratum fracfturing straight well is established；Grid dividing is carried out to the crack of separate stratum fracfturing straight well, obtains multiple crack grids；To each crack grid, the considerations of obtaining the Pressure Drop Model at the fracture network center of a lattice, establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow numerical model, obtain the oil reservoir and flow in fracture coupling model of separate stratum fracfturing straight well；According to the oil reservoir of separate stratum fracfturing straight well and flow in fracture coupling model, the straight well capacity of separate stratum fracfturing is obtained.The present invention can predict separate stratum fracfturing straight well production capacity, and calculating speed is fast, and consider the timeliness of flowing and fracture condudtiviy of the fluid in crack.
Description
Technical field
The present invention relates to the prediction technique of oil reservoir production capacity more particularly to a kind of separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS and dresses
It sets.
Background technique
For fine and close oil as one of unconventional petroleum resources, widely distributed and potentiality to be exploited is huge.But since its reservoir is fine and close,
Permeability is low, poor connectivity, and development difficulty is very big.Especially for multiple reservoirs are contained, reservoir properties and heterogeneity are poor
Compact oil reservoir mostly uses the mining method of staged fracturing vertical well to be exploited to improve production efficiency.But the production practice
Primary construction cannot press off all oil reservoirs, and multiple pressing crack construction can increase considerably individual well pressure break cost again.It is best to obtain
Effect of increasing production uses separate stratum fracfturing --- and longitudinal direction multilayered reservoir serious to heterogeneity carries out separate stratum fracfturing, and is directed to each storage
Layer combines its physical property to distinguish pressure break, to improve fracturing effect.
Currently, for separate zone whole production oil reservoir production capacity research mainly for commingling production straight well (non-pressure break), and for point
Lamination splits well capacity analytic modell analytical model and although calculates fairly simple with form but seldom consider flowing of the fluid in crack, only
It is that each laminar flow is entered to the flow in crack to be simply superimposed, obtains total flow；Although existing numerical model can align
Still modeling process is complicated for the characteristics of really describing separate stratum fracfturing straight well, needs to carry out grid dividing, crack to reservoir and crack
With reservoir differing greatly due to size, tend not to using a nested grid (grid of same specification), so in grid division
When generally require to be encrypted near crack, in order to guarantee model and simulate precision, generally also need to carry out logarithm
Encryption, therefore calculating speed is slow.On the other hand the existing productivity model for separate stratum fracfturing straight well, seldom consideration fracture guide
Can ability change with time this actual conditions, cause calculated result and actual production result with biggish discrepancy.
Summary of the invention
The embodiment of the present invention proposes a kind of separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS, to predict that separate stratum fracfturing straight well produces
Can, calculating speed is fast, and considers flowing and fracture condudtiviy timeliness of the fluid in crack, this method comprises:
According to the oil deposit parameter of separate stratum fracfturing straight well, the reservoir model of separate stratum fracfturing straight well is established；
According to the fracture parameters of separate stratum fracfturing straight well, the fisstured flow model of separate stratum fracfturing straight well is established；
Grid dividing is carried out to the crack of separate stratum fracfturing straight well, obtains multiple crack grids；
To each crack grid, according to the reservoir model of separate stratum fracfturing straight well, obtain at the fracture network center of a lattice
Pressure Drop Model；
To each crack grid, calculated according to the fisstured flow model of separate stratum fracfturing straight well and fracture condudtiviy timeliness
Model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow numerical model；
According to the considerations of Pressure Drop Model and all slits grid at all slits net center of a lattice when fracture condudtiviy
The fisstured flow numerical model of the discretization of effect property obtains the oil reservoir and flow in fracture coupling model of separate stratum fracfturing straight well；
According to the oil reservoir of separate stratum fracfturing straight well and flow in fracture coupling model, the straight well capacity of separate stratum fracfturing is obtained.
The embodiment of the present invention proposes a kind of separate stratum fracfturing straight well capability forecasting device, to predict that separate stratum fracfturing straight well produces
Can, calculating speed is fast, and considers flowing and fracture condudtiviy timeliness of the fluid in crack, which includes:
Reservoir model establishes module, for the oil deposit parameter according to separate stratum fracfturing straight well, establishes separate stratum fracfturing straight well
Reservoir model；
Fisstured flow model building module establishes separate stratum fracfturing straight well for the fracture parameters according to separate stratum fracfturing straight well
Fisstured flow model；
Grid dividing module carries out grid dividing for the crack to separate stratum fracfturing straight well, obtains multiple crack grids；
First crack grid computing module is used for each crack grid, according to the flow through oil reservoir mould of separate stratum fracfturing straight well
Type obtains the Pressure Drop Model at the fracture network center of a lattice；
Second crack grid computing module is used for each crack grid, according to the fisstured flow mould of separate stratum fracfturing straight well
Type and fracture condudtiviy timeliness computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization
Fisstured flow numerical model；
Oil reservoir and flow in fracture coupling model obtain module, for according to the pressure drop mould at all slits net center of a lattice
The considerations of type and all slits grid fracture condudtiviy timeliness discretization fisstured flow numerical model, be layered
The oil reservoir and flow in fracture coupling model of pressure break straight well；
AOF calculation module obtains layering pressure for the oil reservoir and flow in fracture coupling model according to separate stratum fracfturing straight well
Split straight well capacity.
The embodiment of the present invention also proposed a kind of computer equipment, including memory, processor and storage are on a memory
And the computer program that can be run on a processor, the processor realize above-mentioned separate stratum fracfturing when executing the computer program
Straight well PRODUCTION FORECASTING METHODS.
The embodiment of the present invention also proposed a kind of computer readable storage medium, the computer-readable recording medium storage
There is the computer program for executing above-mentioned separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS.
In embodiments of the present invention, according to the oil deposit parameter of separate stratum fracfturing straight well, the oil reservoir for establishing separate stratum fracfturing straight well seeps
Flow model；According to the fracture parameters of separate stratum fracfturing straight well, the fisstured flow model of separate stratum fracfturing straight well is established；It is straight to separate stratum fracfturing
The crack of well carries out grid dividing, obtains multiple crack grids；To each crack grid, seeped according to the oil reservoir of separate stratum fracfturing straight well
Flow model obtains the Pressure Drop Model at the fracture network center of a lattice；To each crack grid, according to the crack of separate stratum fracfturing straight well
Flow model in porous media and fracture condudtiviy timeliness computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness
The fisstured flow numerical model of discretization；According to the Pressure Drop Model and all slits grid at all slits net center of a lattice
Consider the fisstured flow numerical model of the discretization of fracture condudtiviy timeliness, obtains oil reservoir and the crack of separate stratum fracfturing straight well
Coupled with Flow model；According to the oil reservoir of separate stratum fracfturing straight well and flow in fracture coupling model, the straight well capacity of separate stratum fracfturing is obtained.
In embodiments of the present invention, crack carries out grid dividing, without carrying out grid dividing, modeling process letter to oil reservoir
Single, calculating speed is fast, and to each crack grid, when according to the fisstured flow model and fracture condudtiviy of separate stratum fracfturing straight well
Effect property computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow Numerical-Mode
Type, it is contemplated that flowing and fracture condudtiviy timeliness of the fluid in crack, computational accuracy are high.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.In the accompanying drawings:
Fig. 1 is the flow chart of separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS in the embodiment of the present invention；
Fig. 2 is the physical model diagram of separate stratum fracfturing straight well in the embodiment of the present invention；
Fig. 3 is the discrete schematic diagram in single crack in the embodiment of the present invention；
Fig. 4 is the conflux schematic diagram of the separate stratum fracfturing straight well near wellbore in the embodiment of the present invention；
Fig. 5 is the schematic diagram of oil reservoir and crack coupling model in the embodiment of the present invention；
Fig. 6 is the capability forecasting comparative result figure using method proposed by the present invention and Eclipse software；
Fig. 7 is to obtain capability forecasting result and practical comparison of production figure using method proposed by the present invention；
Fig. 8 is the AOF calculation result in the embodiment of the present invention under different coefficient of permeability variation；
Fig. 9 is the structural schematic diagram of separate stratum fracfturing straight well capability forecasting device in the embodiment of the present invention.
Specific embodiment
Understand in order to make the object, technical scheme and advantages of the embodiment of the invention clearer, with reference to the accompanying drawing to this hair
Bright embodiment is described in further details.Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
Fig. 1 is the flow chart of separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS in the embodiment of the present invention, as shown in Figure 1, this method
Include:
Step 101, according to the oil deposit parameter of separate stratum fracfturing straight well, the reservoir model of separate stratum fracfturing straight well is established；
Step 102, according to the fracture parameters of separate stratum fracfturing straight well, the fisstured flow model of separate stratum fracfturing straight well is established；
Step 103, grid dividing is carried out to the crack of separate stratum fracfturing straight well, obtains multiple crack grids；
Step 104, to each crack grid, according to the reservoir model of separate stratum fracfturing straight well, the crack grid is obtained
Center at Pressure Drop Model；
Step 105, to each crack grid, when according to the fisstured flow model and fracture condudtiviy of separate stratum fracfturing straight well
Effect property computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow Numerical-Mode
Type；
Step 106, it is led according to crack the considerations of the Pressure Drop Model and all slits grid at all slits net center of a lattice
The fisstured flow numerical model of the discretization of stream ability timeliness obtains the oil reservoir and flow in fracture coupled mode of separate stratum fracfturing straight well
Type；
Step 107, according to the oil reservoir of separate stratum fracfturing straight well and flow in fracture coupling model, the production of separate stratum fracfturing straight well is obtained
Energy.
In embodiments of the present invention, crack carries out grid dividing, without carrying out grid dividing, modeling to oil reservoir
Process is simple, and calculating speed is fast, and to each crack grid, according to the fisstured flow model and fracture guide of separate stratum fracfturing straight well
Ability timeliness computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow number
It is worth model, it is contemplated that flowing and fracture condudtiviy timeliness of the fluid in crack, computational accuracy are high.
In one embodiment, according to the oil deposit parameter of separate stratum fracfturing straight well, the flow through oil reservoir mould of separate stratum fracfturing straight well is established
Type may include:
According to the oil deposit parameter of separate stratum fracfturing straight well, the elastic instability seepage flow number of each layer of separate stratum fracfturing straight well is established
Learn model；
According to elastic instability seepage experiment, the face source of the unit strength of oil reservoir exists where obtaining separate stratum fracfturing straight well
Pressure-drop model at arbitrary point.
When it is implemented, the oil deposit parameter of the separate stratum fracfturing straight well may include: each layer of Reservoir Permeability, each
The oil reservoir conductivity of layer, reservoir fluid viscosity, each layer of reservoir pore degree and oil reservoir system compressibility.
When it is implemented, be primarily based on the semi-analytic method of flow through oil reservoir problem solving, provide it is assumed hereinafter that:
First, the separate stratum fracfturing straight well is in a homogeneous anisotropic oil reservoir, and the reservoir properties of each layer are different；
Second, the fluid individual event in oil reservoir is micro- compressible；
Third is carved at the beginning, and the pressure of same layer oil reservoir is equal everywhere, is each layer of initial reservoir pressure, different
The initial reservoir pressure of layer is different；
4th, have in the oil reservoir and the only described separate stratum fracfturing straight well, the vertical section of the separate stratum fracfturing straight well is pressed
It splits, it is each to be laminated out a crack perpendicular to the separate stratum fracfturing straight well；
5th, there are interlayer between each layer oil reservoir, reservoir and crack in each layer are not communicated with each other, and flow event is not mutual
It influences.
6th, each layer fluid all first flows into crack from stratum, then flows into separate stratum fracfturing straight well along crack.
7th, the pit shaft of separate stratum fracfturing straight well is water conservancy diversion, and each layer bottom pressure is identical, is separate stratum fracfturing straight well inside crack
One-dimensional Darcy Flow, consider the epidermis that generates near separate stratum fracfturing straight well due to conflux.
Fig. 2 is the physical model diagram of separate stratum fracfturing straight well in the embodiment of the present invention, as shown in Fig. 2, the separate stratum fracfturing is straight
Well includes 4 layers, and there are interlayers between every layer.
On the basis of assumed above, according to the oil deposit parameter of separate stratum fracfturing straight well, separate stratum fracfturing straight well below is established
Each layer of elastic instability seepage experiment:
Wherein, x, y, z are respectively three directions of the separate stratum fracfturing straight well on 3-D walls and floor；
kx,ky, kzThe permeability in three directions respectively on 3-D walls and floor；
ηx,ηy, ηzRespectively correspond to x, the conductivity in tri- directions y, z, unit m2/ s can characterize the unit time
The speed that upper pressure wave is propagated in this direction；
μ is the viscosity of fluid, unit Pas；
For reservoir pore degree；
ctFor the system compressibility of oil reservoir, unit Pa-1；
P is the pressure of somewhere at a time in stratum, unit Pa.
The primary condition of formula (1.1) is as follows:
P (x, y, z, 0)=pIi (1.2)
Wherein, pIiFor the initial reservoir pressure of tier I in stratum, unit Pa.
The boundary condition of formula (1.1) is as follows:
The internal boundary condition of formula (1.1) is as follows:
Wherein,
Q (t) is the production capacity of separate stratum fracfturing straight well at a time, unit m3/s；
H is the depth on the direction z, unit m；
rwFor the radius of separate stratum fracfturing straight well, unit m；
K is the permeability k in three directions on 3-D walls and floorx,ky, kzGeometrical mean.
In one embodiment, according to elastic instability seepage experiment, layering pressure can be obtained using method of Source
The pressure-drop model of the face source of the unit strength of oil reservoir at any point where splitting straight well.
Formula (1.1) be one have time variable boundary condition elastic instability Seepage problems, method of Source by
Gringarten Ramey is proposed, can be used for solving problem above, when it is implemented, can be calculated using following formula, is obtained
Intensity is the face source (x of q (τ)i,yi,zi) at (x, y, z) generate pressure-drop model:
Wherein, △ p is the face source (x that intensity is q (τ)i,yi,zi) pressure drop at (x, y, z)；
Sj(x, y, z) is the face source of unit intensity；
xw, yw, zwRespectively SjThree values of (x, y, z) in the centre coordinate of three-dimensional coordinate；
xf, zfRespectively SjThe length and height of (x, y, z)；
xe, ye, zeOil reservoir is respectively indicated in the corresponding size in three directions.
According to formula (1.5) it is found that multiple face sources can be calculated simultaneously to generation at (x, y, z) according to principle of stacking
Pressure drop, calculation formula are as follows:
When it is implemented, the fracture parameters of the separate stratum fracfturing straight well may include: fracture length, fracture condudtiviy and
The fracture condudtiviy coefficient of variation.
When it is implemented, the seepage flow in crack can since the height and width in crack are smaller with respect to for the length in crack
To regard the one-dimensional Darcy Flow in fracture length direction as.
In one embodiment, according to the fracture parameters of separate stratum fracfturing straight well, the fisstured flow mould of separate stratum fracfturing straight well is established
Type may include:
Consider otherness of the fracture condudtiviy in spatial distribution, and the insertion due to proppant, broken and hole
The timeliness of the low caused fracture condudtiviy of gap pressure drop at any time, obtains in the crack that oil reservoir considered below is flowed into crack
One-dimensional Diffusion Equation it is as follows:
Wherein, qF(x, t) is the fluid volume that some crack grid is flowed into the unit time, unit m3/s；
kF(x, t) is the permeability in t moment crack at x, unit m2；
ydFor slit width, unit m；
CFtFor the system compressibility in crack, unit 1/pa；
For the porosity in crack.
Cinco etc. points out that the compressibility of man-made fracture can be ignored, because of the volume very little in crack.Therefore the stream in crack
It is dynamic may be considered it is incompressible.Under this assumption, above-mentioned formula can simplify, and obtain following separate stratum fracfturing straight well
Fisstured flow model:
In one embodiment, to each crack grid, according to the reservoir model of separate stratum fracfturing straight well, the crack is obtained
Pressure Drop Model at net center of a lattice may include:
To each crack grid, according to the pressure of the face source of the unit strength of oil reservoir where separate stratum fracfturing straight well at any point
Computation model is dropped, obtains pressure drop of each crack grid at the fracture network center of a lattice respectively；
Pressure drop summation of all slits grid at the fracture network center of a lattice is obtained at the fracture network center of a lattice
The computation model of pressure drop.
When it is implemented, not remembering since the width in crack can be ignored compared with height, length, each crack grid
Face source can be regarded as, according to formula (1.5), the intensity in same layer when the source of crack grid j is unit intensity,
The pressure drop of generation at the center of crack grid i can be calculated with following formula:
Since each interlayer is independent of each other, when many cracks produce simultaneously, the pressure of each crack grid is solely dependent upon
The pressure drop that crack grid in this layer generates it distinguishes use space principle of stacking to every crack, according to formula (1.10)
By pressure drop summation of each crack grid at the fracture network center of a lattice, the center (x of crack grid i is obtainedi,yi,zi) at pressure
It is reduced to:
Wherein, InIndicate that the crack grid number that crack divides in tier I is n；
qFjFor the flow of j-th of crack grid, unit m3/s；
p((xj,yj,zj)；(xj,yj,zj)；T) center (x for being crack grid jj,yj,zj) at pressure, it is by same
Influence when other crack grids produce in crack；
MijFor the pressure drop of the generation at the center of crack grid i.
The pressure drop at any crack grid element center can be obtained using Duhamel principle for Variable Mass Flow are as follows:
Reasonable time step-length △ t is chosen, is write as the number about time step in conjunction with formula (1.11) and formula (1.12)
It is worth integrated form, should has to each crack grid i:
Mij(0)=0 (1.14)
Mij(tk-tk-1)=Mij(Δt) (1.15)
Δpi=pini-pi(tk) (1.16)
tk=k × Δ t (1.17)
Wherein, m represents time step number；
K is total time step number of required calculating.
The pressure drop at any time fracture network center of a lattice can be calculated according to formula (1.13).
In one embodiment, to each crack grid, according to the fisstured flow model and fracture guide of separate stratum fracfturing straight well
Ability timeliness computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow number
It is worth model, may include:
The experiment of crack long-term flow conductivity is carried out, determines the fracture condudtiviy coefficient of variation；
According to the fracture condudtiviy coefficient of variation, fracture condudtiviy timeliness computation model is obtained；
The discretization of the crack grid is established according to the fisstured flow model of separate stratum fracfturing straight well to each crack grid
Fisstured flow numerical model；
To each crack grid, according to the fisstured flow numerical model and fracture guide energy of the discretization of the crack grid
Power timeliness computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow numerical value
Model.
When it is implemented, multiple crack long-term flow conductivity experiment can be carried out, the fracture condudtiviy coefficient of variation is determined,
Then according to the fracture condudtiviy coefficient of variation, following fracture condudtiviy timeliness computation model is obtained:
[kFyd]i,n=(kF0yd0)i,n(1-βlg(n/86400+1)) (1.18)
Wherein, β is fracture condudtiviy or coefficient of permeability variation, and zero dimension, value is generally 0.2~0.3；
kFoydoFor incipient crack flow conductivity, unit m2·m；
N is moment, unit s.
Each crack grid seeps the crack for the separate stratum fracfturing straight well that formula (1.9) indicate using centered Finite Difference Methods
Flow model carries out sliding-model control, obtains the fully implicit scheme format of formula (1.9) are as follows:
To the isometric net in the rectangular coordinate system after grid dividing, there are following relationships:
Simultaneous formula (1.19) and formula (1.20), the fisstured flow numerical model of available following discretization:
Wherein, [kF yd]i,nFor the fracture condudtiviy of moment n crack grid i, unit m2·m；
The flow velocity of the fluid of grid i, unit m are flowed into for moment n3/s。
Simultaneous formula (1.21) and formula (1.18), the considerations of obtaining the following crack grid fracture condudtiviy timeliness
The fisstured flow numerical model of the discretization of property:
Fig. 3 is the discrete schematic diagram in single crack in the embodiment of the present invention, as shown in figure 3, crack can be divided into multiple split
Grid is stitched, the fisstured flow model of the discretization of separate stratum fracfturing straight well is calculated based on the grid of crack.
When fracture network center of a lattice is connected with separate stratum fracfturing straight well, as shown in figure 4, Fig. 4 is in the embodiment of the present invention
The conflux schematic diagram of separate stratum fracfturing straight well near wellbore, when the flowing that fracture network center of a lattice flows into separate stratum fracfturing straight well is treated as
One-dimensional Darcy Flow, and consider that issuable epidermis building-up effect near separate stratum fracfturing straight well is connected with separate stratum fracfturing straight well
The considerations of crack grid fracture condudtiviy timeliness discretization fisstured flow numerical model it is as follows:
Wherein, scFor the epidermis as caused by conflux effect, zero dimension；
pwfFor flowing bottomhole pressure (FBHP)；
For the flow at n moment at i-th of grid, m3/s。
In one embodiment, according to the Pressure Drop Model and all slits grid at all slits net center of a lattice the considerations of
The fisstured flow numerical model of the discretization of fracture condudtiviy timeliness obtains the oil reservoir and flow in fracture of separate stratum fracfturing straight well
Coupling model may include:
To each crack grid, according to the Pressure Drop Model at the fracture network center of a lattice, the fracture network center of a lattice is obtained
The pressure at place；With the pressure at the fracture network center of a lattice, the considerations of replacing the crack grid fracture condudtiviy timeliness
Pressure in the fisstured flow numerical model of discretization, the considerations of obtaining the replaced crack grid fracture condudtiviy timeliness
The fisstured flow numerical model of the discretization of property；
By all slits net center of a lattice Pressure Drop Model and replaced all slits grid the considerations of fracture guide
The fisstured flow numerical model of the discretization of ability timeliness forms the oil reservoir and flow in fracture coupled mode of separate stratum fracfturing straight well
Type.
When it is implemented, at using the oil reservoir net center of a lattice at the available oil reservoir of formula (1.16) and crack interface
Pressure p1,p2,...,pi,...,pNThe relationship between flow.
Assuming that the discretization of the considerations of formula (1.23) denotable crack grid fracture condudtiviy timeliness is split
Stitching the pressure in seepage flow numerical model is pF1,pF2,...,pFi,...,pFN, from oil reservoir net center of a lattice into crack grid
Flow q at the heartF1,qF2,......,qFN.Due to the continuity of the pressure and flow that go out in interface, oil reservoir net center of a lattice
Locate the equal p of pressure at pressure and fracture network center of a latticei=pFi。
Therefore with the pressure p at the fracture network center of a lattice1,p2,...,pi,...,pN, replace pF1,pF2,...,
pFi,...,pFN, the considerations of obtaining the replaced crack grid fracture condudtiviy timeliness discretization fisstured flow number
It is worth model to get to about qF1,qF2,......,qFN, pF1,pF2,......,pFNTotal 2N formula, qF1,qF2,......,
qFN, pF1,pF2,......,pFNThe equation group of 2N of composition closing formula composition, the as oil reservoir of separate stratum fracfturing straight well with split
Coupled with Flow model is stitched, Fig. 5 is the schematic diagram of oil reservoir and crack coupling model in the embodiment of the present invention, as shown in figure 5, the oil reservoir
It include 4 layers of stratum, 4 cracks with crack coupling model.
In one embodiment, according to the oil reservoir of separate stratum fracfturing straight well and flow in fracture coupling model, it is straight to obtain separate stratum fracfturing
Well capacity may include:
Obtain the production capacity of each crack grid in oil reservoir and flow in fracture coupling model；
Production capacity summation to all slits grid, obtains the straight well capacity of separate stratum fracfturing.
A specific embodiment is given below, illustrates the concrete application of separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS of the invention.
Firstly, according to each layer of Reservoir Permeability, each layer of oil reservoir conductivity, reservoir fluid viscosity, each layer
The oil deposit parameter of reservoir pore degree and oil reservoir system compressibility hierarchical pressure break straight well is established and is layered shown in formula (1.1)
The reservoir model of pressure break straight well.
Then the reservoir model of the separate stratum fracfturing straight well according to shown in formula (1.1) obtains shown in formula (1.5)
Intensity is the face source (x of q (τ)i,yi,zi) at (x, y, z) generate pressure-drop model.
Grid dividing is carried out to 4 cracks of separate stratum fracfturing straight well, it is assumed that every layer of each crack is divided into 8 fracture networks
Lattice, i.e., include 32 crack grids altogether, and oil reservoir and the flow in fracture coupling model of the separate stratum fracfturing straight well eventually formed then wrap
32 × 2=64 formula is included, i.e. unknown quantity is that pressure is p at fracture network center of a latticeF1,pF2,...,pFi,...,pFNWith oil reservoir
To the flow q at fracture network center of a lattice at net center of a latticeF1,qF2,......,qFN64 closing formula compositions of composition
The expression formula of equation group, equation group can write a Chinese character in simplified form are as follows:
A(64×64)X(641)=B(64×1) (1.24)
Wherein, A is the coefficient matrix of unknown quantity, specific as follows:
Matrix A is divided into 8 block matrix, I is unit matrix, in representation formula (1.13) at the net center of a lattice of induced cleavage
Pressure pF1,pF2,......,pFNThe coefficient of front, available:
MIIt (I=1,2,3,4) is at tier I oil reservoir net center of a lattice to the flow q at fracture network center of a latticeF1,
qF2,......,qFNThe coefficient of front, available in conjunction with formula (1.13):
PI(I=1,2,3,4) pressure p at the net center of a lattice of tier I induced cleavage is indicatedF1,pF2,......,pFNFront
Coefficient, available in conjunction with formula (1.2):
Q in formula (1.22)F1,qF2,…qFNCorresponding matrix is
X is the column vector being made of 64 unknown quantitys: X(64×1)=[X1X2X3X4], concrete form is as follows:
XI(16×1)=[p1 p2 ...... p8 q1 q2 ...... q8] I=1,2,3,4
B is the constant term of oil reservoir and flow in fracture coupling model right end: B(64×1)=[B1 B2 B3 B4], concrete form is such as
Under:
Wherein, T indicates transposition；
dpiRepresent the right-hand vector of formula (1.13).
It should be noted that dp when first time stepi=0, remaining time step dpi(i > 1) only need to be by preceding i- when calculating
The flow of the corresponding crack grid of 1 time step is brought formula (1.13) into and can be calculated.
The 4 matrix M finally obtained according to 64 formula in oil reservoir and flow in fracture coupling modelI、PI, A, B, simultaneous
4 matrixes, solve system of equation (1.24), so that it may obtain the flow q of each crack gridF1,qF2,…qFNAnd pressure pF1,
pF2,......,pFN。
By all slits grid qF1,qF2,…qFNFlow be added, can get the straight well capacity of separate stratum fracfturing.
It, will be by the way that the present invention be counted in order to verify the validity of separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS proposed by the present invention
The result that the result and business numerical simulation software of calculation calculate compares.When comparing, used parameter such as 1 institute of table
Show.Business numerical simulation software selects Eclipse, the side which is combined using local mesh reflnement and equivalent conductivity
Method, to construct the tiny grid that can simulate true flow in fracture, the encryption rule of grid is logarithm encryption.
1 separate stratum fracfturing straight well capability forecasting calculating parameter of table
Fig. 6 is the capability forecasting comparative result figure using method proposed by the present invention and Eclipse software, and Fig. 7 is to use
Method proposed by the present invention obtains capability forecasting result and practical comparison of production figure, and the present invention mentions it can be seen from Fig. 6 and Fig. 7
Method out and Eclipse software and the production capacity goodness of fit all with higher, to verify effectiveness of the invention and accurate
Property.
Method proposed by the present invention is compared with business numerical simulation softwares such as Eclipse, and advantage is, proposed by the present invention
Method when solving flow through oil reservoir using by source function based on analytic solutions, only flow in fracture use numerical value
Solution, therefore the model of this paper only needs fracture to carry out grid dividing, without carrying out grid dividing, a side to entire oil reservoir
The scale of grid has greatly been lacked in face, on the other hand also substantially increases the convergence during model solution, therefore can obtain
Obtain calculating speed more faster than numerical model.
Fig. 8 is the AOF calculation in the embodiment of the present invention under different coefficient of permeability variation as a result, as shown in figure 8, to dividing
Lamination is split for straight well, and the variation of fracture condudtiviy is affected to production capacity, is considered as crack when carrying out capability forecasting
Flow conductivity changes with time, to more accurately predict production capacity, can be tested by long-term flow conductivity in actual production
Experiment determines that fracture condudtiviy changes with time rule, method Accurate Prediction proposed by the present invention is then recycled to be layered
The production history situation of fractured well.
In embodiments of the present invention, according to the oil deposit parameter of separate stratum fracfturing straight well, the oil reservoir for establishing separate stratum fracfturing straight well seeps
Flow model；According to the fracture parameters of separate stratum fracfturing straight well, the fisstured flow model of separate stratum fracfturing straight well is established；It is straight to separate stratum fracfturing
The crack of well carries out grid dividing, obtains multiple crack grids；To each crack grid, seeped according to the oil reservoir of separate stratum fracfturing straight well
Flow model obtains the Pressure Drop Model at the fracture network center of a lattice；To each crack grid, according to the crack of separate stratum fracfturing straight well
Flow model in porous media and fracture condudtiviy timeliness computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness
The fisstured flow numerical model of discretization；According to the Pressure Drop Model and all slits grid at all slits net center of a lattice
Consider the fisstured flow numerical model of the discretization of fracture condudtiviy timeliness, obtains oil reservoir and the crack of separate stratum fracfturing straight well
Coupled with Flow model；According to the oil reservoir of separate stratum fracfturing straight well and flow in fracture coupling model, the straight well capacity of separate stratum fracfturing is obtained.
In embodiments of the present invention, crack carries out grid dividing, without carrying out grid dividing, modeling process letter to oil reservoir
Single, calculating speed is fast, and to each crack grid, when according to the fisstured flow model and fracture condudtiviy of separate stratum fracfturing straight well
Effect property computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow Numerical-Mode
Type, it is contemplated that flowing of the fluid in crack, computational accuracy are high.
In addition, the embodiment of the present invention considers otherness of the fracture condudtiviy in spatial distribution and time-varying at any time
Property, it can more accurately predict production capacity, can determine fracture condudtiviy by long-term flow conductivity test experiments in actual production
Change with time rule, then recycles method Accurate Prediction separate stratum fracfturing well capacity situation of change proposed by the present invention.
Based on the same inventive concept, the embodiment of the invention also provides a kind of separate stratum fracfturing straight well capability forecasting devices, such as
Described in following implementation.Since the principle that these are solved the problems, such as is similar to separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS, device
Implementation may refer to the implementation of method, repeat place and do not repeating.
Fig. 9 is the structural schematic diagram of separate stratum fracfturing straight well capability forecasting device in the embodiment of the present invention, as shown in figure 9, should
Device includes:
Reservoir model establishes module 901, and for the oil deposit parameter according to separate stratum fracfturing straight well, it is straight to establish separate stratum fracfturing
The reservoir model of well；
It is straight to establish separate stratum fracfturing for the fracture parameters according to separate stratum fracfturing straight well for fisstured flow model building module 902
The fisstured flow model of well；
Grid dividing module 903 carries out grid dividing for the crack to separate stratum fracfturing straight well, obtains multiple fracture networks
Lattice；
First crack grid computing module 904 is used for each crack grid, according to the flow through oil reservoir of separate stratum fracfturing straight well
Model obtains the Pressure Drop Model at the fracture network center of a lattice；
Second crack grid computing module 905 is used for each crack grid, according to the fisstured flow of separate stratum fracfturing straight well
Model and fracture condudtiviy timeliness computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness it is discrete
The fisstured flow numerical model of change；
Oil reservoir and flow in fracture coupling model obtain module 906, for according to the pressure drop at all slits net center of a lattice
The considerations of model and all slits grid fracture condudtiviy timeliness discretization fisstured flow numerical model, divided
The oil reservoir and flow in fracture coupling model of straight well are split in lamination；
AOF calculation module 907 is layered for the oil reservoir and flow in fracture coupling model according to separate stratum fracfturing straight well
The straight well capacity of pressure break.
In conclusion in embodiments of the present invention, according to the oil deposit parameter of separate stratum fracfturing straight well, establishing separate stratum fracfturing straight well
Reservoir model；According to the fracture parameters of separate stratum fracfturing straight well, the fisstured flow model of separate stratum fracfturing straight well is established；To point
The crack that straight well is split in lamination carries out grid dividing, obtains multiple crack grids；To each crack grid, according to separate stratum fracfturing straight well
Reservoir model, obtain the Pressure Drop Model at the fracture network center of a lattice；It is straight according to separate stratum fracfturing to each crack grid
The fisstured flow model and fracture condudtiviy timeliness computation model of well, the considerations of establishing the crack grid fracture condudtiviy
The fisstured flow numerical model of the discretization of timeliness；According to the Pressure Drop Model at all slits net center of a lattice and all split
The fisstured flow numerical model of the discretization of the considerations of stitching grid fracture condudtiviy timeliness, obtains the oil of separate stratum fracfturing straight well
Hiding and flow in fracture coupling model；According to the oil reservoir of separate stratum fracfturing straight well and flow in fracture coupling model, it is straight to obtain separate stratum fracfturing
Well capacity.In embodiments of the present invention, crack carries out grid dividing, without carrying out grid dividing to oil reservoir, builds
Mold process is simple, and calculating speed is fast, and to each crack grid, is led according to the fisstured flow model of separate stratum fracfturing straight well and crack
Stream ability timeliness computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow
Numerical model, it is contemplated that flowing of the fluid in crack, computational accuracy are high.
In addition, the embodiment of the present invention considers otherness of the fracture condudtiviy in spatial distribution and time-varying at any time
Property, it can more accurately predict production capacity, can determine fracture condudtiviy by long-term flow conductivity test experiments in actual production
Change with time rule, then recycles method Accurate Prediction separate stratum fracfturing well capacity situation of change proposed by the present invention.
It should be understood by those skilled in the art that, the embodiment of the present invention can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the present invention, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects
Describe in detail it is bright, it should be understood that the above is only a specific embodiment of the present invention, the guarantor being not intended to limit the present invention
Range is protected, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in this
Within the protection scope of invention.
Claims (10)
1. a kind of separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS characterized by comprising
According to the oil deposit parameter of separate stratum fracfturing straight well, the reservoir model of separate stratum fracfturing straight well is established；
According to the fracture parameters of separate stratum fracfturing straight well, the fisstured flow model of separate stratum fracfturing straight well is established；
Grid dividing is carried out to the crack of separate stratum fracfturing straight well, obtains multiple crack grids；
To each crack grid, according to the reservoir model of separate stratum fracfturing straight well, the pressure at the fracture network center of a lattice is obtained
Model drops；
To each crack grid, mould is calculated according to the fisstured flow model of separate stratum fracfturing straight well and fracture condudtiviy timeliness
Type, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow numerical model；
According to fracture condudtiviy timeliness the considerations of Pressure Drop Model and all slits grid at all slits net center of a lattice
Discretization fisstured flow numerical model, obtain separate stratum fracfturing straight well oil reservoir and flow in fracture coupling model；
According to the oil reservoir of separate stratum fracfturing straight well and flow in fracture coupling model, the straight well capacity of separate stratum fracfturing is obtained.
2. separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS as described in claim 1, which is characterized in that according to separate stratum fracfturing straight well
Oil deposit parameter establishes the reservoir model of separate stratum fracfturing straight well, comprising:
According to the oil deposit parameter of separate stratum fracfturing straight well, the elastic instability seepage flow mathematical modulo of each layer of separate stratum fracfturing straight well is established
Type；
According to elastic instability seepage experiment, the face source of the unit strength of oil reservoir is any where obtaining separate stratum fracfturing straight well
Pressure-drop model at point.
3. separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS as claimed in claim 2, which is characterized in that according to elastic instability seepage flow
Mathematical model, the pressure of the face source of the unit strength of oil reservoir at any point where obtaining separate stratum fracfturing straight well using method of Source
Computation model drops.
4. separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS as claimed in claim 2, which is characterized in that each crack grid, root
According to the reservoir model of separate stratum fracfturing straight well, the Pressure Drop Model at the fracture network center of a lattice is obtained, comprising:
To each crack grid, according to the pressure drop meter of the face source of the unit strength of oil reservoir where separate stratum fracfturing straight well at any point
Model is calculated, obtains pressure drop of each crack grid at the fracture network center of a lattice respectively；
To pressure drop summation of all slits grid at the fracture network center of a lattice, the pressure drop at the fracture network center of a lattice is obtained
Computation model.
5. separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS as described in claim 1, which is characterized in that each crack grid, root
According to the fisstured flow model and fracture condudtiviy timeliness computation model of separate stratum fracfturing straight well, the considerations of establishing the crack grid
The fisstured flow numerical model of the discretization of fracture condudtiviy timeliness, comprising:
The experiment of crack long-term flow conductivity is carried out, determines the fracture condudtiviy coefficient of variation；
According to the fracture condudtiviy coefficient of variation, fracture condudtiviy timeliness computation model is obtained；
Splitting for the discretization of the crack grid is established according to the fisstured flow model of separate stratum fracfturing straight well to each crack grid
Stitch seepage flow numerical model；
To each crack grid, when according to the fisstured flow numerical model and fracture condudtiviy of the discretization of the crack grid
Effect property computation model, the considerations of establishing the crack grid fracture condudtiviy timeliness discretization fisstured flow Numerical-Mode
Type.
6. separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS as described in claim 1, which is characterized in that according to all slits grid
The considerations of Pressure Drop Model and all slits grid at center fracture condudtiviy timeliness discretization fisstured flow numerical value
Model obtains the oil reservoir and flow in fracture coupling model of separate stratum fracfturing straight well, comprising:
To each crack grid, according to the Pressure Drop Model at the fracture network center of a lattice, obtain at the fracture network center of a lattice
Pressure；With the pressure at the fracture network center of a lattice, the considerations of replacing the crack grid fracture condudtiviy timeliness it is discrete
Pressure in the fisstured flow numerical model of change, the considerations of obtaining the replaced crack grid fracture condudtiviy timeliness
The fisstured flow numerical model of discretization；
By all slits net center of a lattice Pressure Drop Model and replaced all slits grid the considerations of fracture condudtiviy
The fisstured flow numerical model of the discretization of timeliness forms the oil reservoir and flow in fracture coupling model of separate stratum fracfturing straight well.
7. separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS as described in claim 1, which is characterized in that according to separate stratum fracfturing straight well
Oil reservoir and flow in fracture coupling model obtain the straight well capacity of separate stratum fracfturing, comprising:
Obtain the production capacity of each crack grid in oil reservoir and flow in fracture coupling model；
Production capacity summation to all slits grid, obtains the straight well capacity of separate stratum fracfturing.
8. a kind of separate stratum fracfturing straight well capability forecasting device characterized by comprising
Reservoir model establishes module, for the oil deposit parameter according to separate stratum fracfturing straight well, establishes the oil of separate stratum fracfturing straight well
Hide flow model in porous media；
Fisstured flow model building module establishes splitting for separate stratum fracfturing straight well for the fracture parameters according to separate stratum fracfturing straight well
Stitch flow model in porous media；
Grid dividing module carries out grid dividing for the crack to separate stratum fracfturing straight well, obtains multiple crack grids；
First crack grid computing module, for being obtained to each crack grid according to the reservoir model of separate stratum fracfturing straight well
Obtain the Pressure Drop Model at the fracture network center of a lattice；
Second crack grid computing module, for each crack grid, according to the fisstured flow model of separate stratum fracfturing straight well and
Fracture condudtiviy timeliness computation model, the considerations of establishing the crack grid discretization of fracture condudtiviy timeliness split
Stitch seepage flow numerical model；
Oil reservoir and flow in fracture coupling model obtain module, for according to the Pressure Drop Model at all slits net center of a lattice, and
The considerations of all slits grid fracture condudtiviy timeliness discretization fisstured flow numerical model, obtain separate stratum fracfturing it is straight
The oil reservoir and flow in fracture coupling model of well；
It is straight to obtain separate stratum fracfturing for the oil reservoir and flow in fracture coupling model according to separate stratum fracfturing straight well for AOF calculation module
Well capacity.
9. a kind of computer equipment including memory, processor and stores the meter that can be run on a memory and on a processor
Calculation machine program, which is characterized in that the processor is realized described in any one of claim 1 to 7 when executing the computer program
Method.
10. a kind of computer readable storage medium, which is characterized in that the computer-readable recording medium storage has perform claim
It is required that the computer program of any one of 1 to 7 the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910408401.3A CN110080745B (en) | 2019-05-16 | 2019-05-16 | Method and device for predicting productivity of vertical well through separate-layer fracturing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910408401.3A CN110080745B (en) | 2019-05-16 | 2019-05-16 | Method and device for predicting productivity of vertical well through separate-layer fracturing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110080745A true CN110080745A (en) | 2019-08-02 |
| CN110080745B CN110080745B (en) | 2023-04-07 |
Family
ID=67420514
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910408401.3A Active CN110080745B (en) | 2019-05-16 | 2019-05-16 | Method and device for predicting productivity of vertical well through separate-layer fracturing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110080745B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111749687A (en) * | 2020-07-23 | 2020-10-09 | 中海石油国际能源服务(北京)有限公司 | Method, device and equipment for determining main force horizon of multilayer oil reservoir and storage medium |
| CN112818503A (en) * | 2019-11-15 | 2021-05-18 | 中国石油天然气股份有限公司 | Method for determining moving speed of temporary plugging ball |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110307227A1 (en) * | 2010-06-15 | 2011-12-15 | Schlumberger Technology Corporation | Characterizing flow production |
| CN104533370A (en) * | 2014-11-06 | 2015-04-22 | 中国石油大学(北京) | Oil deposit, crack and shaft fully-coupled simulating method of fractured horizontal well |
| CN105089595A (en) * | 2015-05-27 | 2015-11-25 | 中国石油天然气股份有限公司 | Oil reservoir numerical simulation method and device under the action of horizontal fracturing fracture diversion |
| US20160010444A1 (en) * | 2014-07-14 | 2016-01-14 | Saudi Arabian Oil Company | Methods, systems, and computer medium having computer programs stored thereon to optimize reservoir management decisions |
| US20160305238A1 (en) * | 2015-02-06 | 2016-10-20 | Schlumberger Technology Corporation | Average/initial reservoir pressure and wellbore efficiency analysis from rates and downhole pressures |
| CN106545336A (en) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | Consider the Productivity of tight gas reservoir seepage flow mechanism |
| CN107762498A (en) * | 2017-09-27 | 2018-03-06 | 中国地质调查局油气资源调查中心 | A kind of pressure analysis method in the area of tight gas reservoir straight well volume fracturing two |
| CN107956463A (en) * | 2016-10-18 | 2018-04-24 | 中国石油化工股份有限公司 | A kind of ultra-low permeability reservoir straight well omnibearing stereo fracturing process |
| CN108344853A (en) * | 2018-01-24 | 2018-07-31 | 西南石油大学 | Plane heterogeneity constant volume dry gas reservoir straight well absolute open flow (AOF) test method |
| CN108959767A (en) * | 2018-07-02 | 2018-12-07 | 中国地质大学（北京） | A kind of narrow river channel type gas reservoir difference well type condensate injury method for numerical simulation |
| CN109025942A (en) * | 2018-08-09 | 2018-12-18 | 西南石油大学 | A kind of irregular cracky Production rate method of tight gas reservoir Deviated Well Fracturing |
| CN109614736A (en) * | 2018-12-24 | 2019-04-12 | 中国海洋石油集团有限公司 | A kind of determination method of the exploitation of offshore viscous oil field steam soak directional well and straight well production capacity multiple |
-
2019
- 2019-05-16 CN CN201910408401.3A patent/CN110080745B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110307227A1 (en) * | 2010-06-15 | 2011-12-15 | Schlumberger Technology Corporation | Characterizing flow production |
| US20160010444A1 (en) * | 2014-07-14 | 2016-01-14 | Saudi Arabian Oil Company | Methods, systems, and computer medium having computer programs stored thereon to optimize reservoir management decisions |
| CN104533370A (en) * | 2014-11-06 | 2015-04-22 | 中国石油大学(北京) | Oil deposit, crack and shaft fully-coupled simulating method of fractured horizontal well |
| US20160305238A1 (en) * | 2015-02-06 | 2016-10-20 | Schlumberger Technology Corporation | Average/initial reservoir pressure and wellbore efficiency analysis from rates and downhole pressures |
| CN105089595A (en) * | 2015-05-27 | 2015-11-25 | 中国石油天然气股份有限公司 | Oil reservoir numerical simulation method and device under the action of horizontal fracturing fracture diversion |
| CN106545336A (en) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | Consider the Productivity of tight gas reservoir seepage flow mechanism |
| CN107956463A (en) * | 2016-10-18 | 2018-04-24 | 中国石油化工股份有限公司 | A kind of ultra-low permeability reservoir straight well omnibearing stereo fracturing process |
| CN107762498A (en) * | 2017-09-27 | 2018-03-06 | 中国地质调查局油气资源调查中心 | A kind of pressure analysis method in the area of tight gas reservoir straight well volume fracturing two |
| CN108344853A (en) * | 2018-01-24 | 2018-07-31 | 西南石油大学 | Plane heterogeneity constant volume dry gas reservoir straight well absolute open flow (AOF) test method |
| CN108959767A (en) * | 2018-07-02 | 2018-12-07 | 中国地质大学（北京） | A kind of narrow river channel type gas reservoir difference well type condensate injury method for numerical simulation |
| CN109025942A (en) * | 2018-08-09 | 2018-12-18 | 西南石油大学 | A kind of irregular cracky Production rate method of tight gas reservoir Deviated Well Fracturing |
| CN109614736A (en) * | 2018-12-24 | 2019-04-12 | 中国海洋石油集团有限公司 | A kind of determination method of the exploitation of offshore viscous oil field steam soak directional well and straight well production capacity multiple |
Non-Patent Citations (5)
| Title |
|---|
| SHENGYONG HU: "Methane Extraction from Abandoned Mines by Surface Vertical Wells: A Case Study in China", 《HINDAWI》 * |
| WU ZHENGBIN: "A novel model and sensitive analysis for productivity estimate of nitrogen assisted cyclic steam stimulation in a vertical well", 《INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER》 * |
| 史晓东: "致密油直井多层缝网压裂产能预测方法", 《特种油气藏》 * |
| 李准等: "考虑储层应力敏感效应的体积压裂水平井瞬态压力分析", 《东北石油大学学报》 * |
| 李纪: "SDN气藏气井生产动态分析及合理工作制度研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112818503A (en) * | 2019-11-15 | 2021-05-18 | 中国石油天然气股份有限公司 | Method for determining moving speed of temporary plugging ball |
| CN112818503B (en) * | 2019-11-15 | 2022-11-01 | 中国石油天然气股份有限公司 | Method for determining moving speed of temporary plugging ball |
| CN111749687A (en) * | 2020-07-23 | 2020-10-09 | 中海石油国际能源服务(北京)有限公司 | Method, device and equipment for determining main force horizon of multilayer oil reservoir and storage medium |
| CN111749687B (en) * | 2020-07-23 | 2023-11-21 | 中海石油国际能源服务(北京)有限公司 | Multi-layer oil reservoir principal force horizon determination method, device, equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110080745B (en) | 2023-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xu et al. | Development of an embedded discrete fracture model for field-scale reservoir simulation with complex corner-point grids | |
| CN104533370B (en) | Pressure break horizontal well oil reservoir, crack, pit shaft coupled model method | |
| Yushi et al. | Numerical investigation into the influence of bedding plane on hydraulic fracture network propagation in shale formations | |
| CN105260543B (en) | Multi-dielectric oil gas flow simulating method and device based on double porosity model | |
| Gong et al. | Upscaling discrete fracture characterizations to dual-porosity, dual-permeability models for efficient simulation of flow with strong gravitational effects | |
| Wan et al. | Semi-analytical well model of horizontal wells with multiple hydraulic fractures | |
| CN104112057B (en) | A kind of large scale fractured reservoir method for numerical simulation | |
| CN105089595B (en) | Numerical reservoir simulation method and device under the effect of horizontal fractures fracture guide | |
| CN105160134B (en) | The blending agent analogy method and device that oil gas flows in compact reservoir multi-dielectric | |
| CN109241588A (en) | A kind of analogy method of the monolete extension based on quasi-continuous geomechanics model | |
| Zhao et al. | Simulation of a multistage fractured horizontal well in a tight oil reservoir using an embedded discrete fracture model | |
| CN109829217A (en) | Pressure break Fractured Reservoir productivity simulation method and device | |
| Xia et al. | An enriched approach for modeling multiscale discrete-fracture/matrix interaction for unconventional-reservoir simulations | |
| Liu et al. | Manufacturing method of large-scale fractured porous media for experimental reservoir simulation | |
| CN110080745A (en) | Separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS and device | |
| Jinzhou et al. | Prediction of pore pressure–induced stress changes during hydraulic fracturing of heterogeneous reservoirs through coupled fluid flow/geomechanics | |
| Liu et al. | A simplified and efficient method for water flooding production index calculations in low permeable fractured reservoir | |
| Yu et al. | A sand-arch stability constrained dynamic fractal wormhole growth model for simulating cold heavy-oil production with sand | |
| KR101506321B1 (en) | Method for analyzing multi phase and heat flow of fluids in reservoir and Recording media therefor | |
| Mei et al. | Fractal analysis of shape factor for matrix-fracture transfer function in fractured reservoirs | |
| CN107169227B (en) | A kind of the coarse grid analogy method and system of staged fracturing horizontal well | |
| Guo et al. | Theoretical and experimental determination of proppant crushing rate and fracture conductivity | |
| CN107832482B (en) | Compact reservoir multi-scale fracture network modeling and simulation method | |
| Wang et al. | Development and application of a streamline micellar/polymer simulator | |
| Mi et al. | Hydraulic fracture optimization based on discrete fracture network model in shale reservoirs |
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 |