CN109933747A - A kind of quasi- phase infiltration calculation method based on fracture hole mode - Google Patents
A kind of quasi- phase infiltration calculation method based on fracture hole mode Download PDFInfo
- Publication number
- CN109933747A CN109933747A CN201910042454.8A CN201910042454A CN109933747A CN 109933747 A CN109933747 A CN 109933747A CN 201910042454 A CN201910042454 A CN 201910042454A CN 109933747 A CN109933747 A CN 109933747A
- Authority
- CN
- China
- Prior art keywords
- quasi
- fracture hole
- hole mode
- mutually
- calculation method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to the invention discloses a kind of quasi- phases for fracture hole mode to seep calculation method comprising following steps collect the basic data of selected typical fracture hole mode first and establish two dimensional model；According to Navier-Stokes equation simulation water drive oil process, obtains different time and walk upper flow parameter；It seeks intending mutually seeping initial value under the mode using meteod of instability according to analog result；This is mutually seeped to the business simulation softward that initial value substitutes into based on Darcy formula to be fitted, obtains the quasi- mutually infiltration of relatively accurate fracture hole mode.
Description
Technical field
The present invention relates to oil-gas field development technical field, specially a kind of quasi- phase based on fracture hole mode seeps calculation method.
Background technique
Fracture hole oil reservoir efficient accurate method for numerical simulation, method one general at present because its strong heterogeneity lacks always
Kind is to simulate fracture hole medium, although this kind of method concept is simply convenient for calculating, simulation process mistake using dual media method
In simply causing degree of verisimilitude lower；Another method is then that entire fracture-cavity units are regarded as to a disposed of in its entirety, using being similar to
The mode of sandstone oil reservoir is fitted, although the analysis of this method is intuitive, fracture-cavity units stronger for heterogeneity are this
The effect that fit approach reaches excessively coincidence, subsequent prediction accuracy are poor.
For this purpose, the present invention proposes a kind of calculation method mutually seeped for more junior unit i.e. fracture hole mode, the seam being calculated
Hole mode mutually seeps characteristic sub-area and is assigned to fracture-cavity units, greatlys save subsequent fitting predicted time, also improves simulation accuracy.
Summary of the invention
The purpose of the present invention is to provide a kind of, and the quasi- phase based on fracture hole mode seeps calculation method, passes through N-S equation calculation
Fracture hole mode water drive oil dynamic parameter, the quasi- phase for recycling dynamic parameter to calculate fracture hole mode seep initial value, which is substituted into number
It is fitted Production development in value simulation softward and obtains relatively accurate quasi- phase infiltration, is suitable for fracture hole reservoir numerical simulation and parameter is excellent
Change.
To achieve the above object, the invention provides the following technical scheme: a kind of quasi- phase based on fracture hole mode seeps calculating side
Method comprising following steps:
S1 collects the basic data of selected typical fracture hole mode and establishes two dimensional model；
S2 obtains different time and walks upper flow parameter according to Navier-Stokes equation simulation water drive oil process；
S3 seeks intending mutually seeping initial value under the mode using meteod of instability according to analog result；
This is mutually seeped the business simulation softward that initial value substitutes into based on Darcy formula and is fitted by S4, is obtained relatively accurate
Fracture hole mode quasi- mutually seep；
S5 converts fracture hole pattern configurations parameter, repeats step S2-S4, obtains the quasi- phase under the fracture hole mode different shape
It seeps.
Preferably, the concrete operations of the step S1 are: according to oil field actual conditions, extracting Reasonable Injection and adopt model, by it
It is mapped to two-dimensional surface, constructs model.
Preferably, the concrete operations of the step S2 are: using Navier-Stokes equation analysis flow regime, use
Level Set Method handles oil water interfacial tension problem, spreads interface problem using Cahn-Hilliard equation control water-oil phase.
Preferably, the concrete operations of the step S3 are: model being equivalent to a rock core, according to equisaturation movement side
Journey and B-L equation calculate water phase relative permeability and oil relative permeability under corresponding saturation degree.
Preferably, the concrete operations of the step S4 are: the mutually infiltration initial value that step S3 is obtained being substituted into business software and is transported
It calculates, analog parameter and fitting parameter are consistent with step S2 operational parameter, and fitting Production development result obtains relatively accurate fracture hole
Mode is quasi- mutually to seep.
Preferably, the concrete operations of the step S5 are: transformation fracture aperture, and the parameters such as solution cavity size proportion repeat to walk
Rapid S2-S4 calculates the quasi- mutually infiltration of the fracture hole mode under different shape.
Detailed description of the invention
Fig. 1 is quasi- phase infiltration calculation method interwell communication mould of certain block T402-TK429C based on fracture hole mode in the present invention
Type；
Fig. 2 is double holes double slit pattern diagram in the present invention；
Fig. 3 is to intend mutually seeping initial value operation result in the present invention；
Fig. 4 is Production development fitting result in the present invention；
Fig. 5 is to intend mutually seeping final result in the present invention；
Fig. 6, which is that different solution cavity radius ratios are corresponding in the present invention, to be intended mutually seeping result；
Fig. 7, which is that different fracture apertures are quasi- in the present invention, mutually seeps result；
Fig. 8, which is that different fracture lengths are quasi- in the present invention, mutually seeps result；
Fig. 9, which is that different fracture dips are quasi- in the present invention, mutually seeps result.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Fig. 1-9 is please referred to, the present invention provides a kind of technical solution: a kind of quasi- phase infiltration calculation method based on fracture hole mode,
It includes the following steps:
Step S1 extracts Reasonable Injection and adopts model, map that two-dimensional surface, construct model according to oil field actual conditions.
Step S2 is to handle oil-water interfaces using Level Set Method using Navier-Stokes equation analysis flow regime
Strain problem spreads interface problem using Cahn-Hilliard equation control water-oil phase.
The quality and momentum of Navier-Stokes equation simulation fluid of the model based on incompressible fluid are transmitted.Consider
The influence of capillary force, it is necessary in a model include surface tension；Modelling is longitudinal profile, summarizes addition gravity in equation
?.Therefore, Navier-Stokes equation is expressed as formula (1), (2):
▽ u=0 (2)
Wherein, ρ represents fluid density, kg/m3；μ indicates fluid kinematic viscosity, Ns/m2；U indicates speed, m/s；P is indicated
Pressure, Pa；G indicates acceleration of gravity, m/s2；FstIndicate the interfacial tension of oil-water interfaces.
It is calculated shown in interfacial tension such as formula (3), (4) using Level Set Method:
Fst=▽ T (3)
T=σ (I- (nnT))δ (4)
Wherein, I indicates known matrix, and n is interface normal, and σ indicates interfacial tension coefficient (N/m), and δ indicates dirac letter
Number is only non-zero number on fluid boundary.When using finite element method Navier-Stokes equation, equation is multiplied
It is then integrated in computational domain with test function.If handled for part, test can will be moved to the derivative of T
In equation, such result is exactly that the integral of entire computational domain adds the integral of graphic limit.
For the diffusion interface of tracking water phase and oily phase, controlled using Cahn-Hilliard equation.Diffusion interface definition is nothing
The region of dimension phase field variable φ from -1 to 1.For convenience of solution, Cahn-Hilliard equation is divided into two equations:
ψ=- ▽ ε2▽φ+(φ2-1)φ (6)
Wherein, u is fluid velocity, m/s；γ is mobility, m3·s/kg；λ is mixed tensor density N；ε is interfacial thickness
Parameter, m；ψ variable becomes phase field master variable.Following equation (7) is by mixed tensor density and interfacial thickness and surface tension coefficient
It is associated:
Usually interfacial thickness parameter can be set as ε=hc/ 2, wherein hcIt is the feature grid ruler in the region of interface process
It is very little.Mobility parameters γ determines the time scale of Cahn-Hilliard diffusion, it is necessary to select with caution, which needs sufficiently large
To guarantee that interfacial thickness is constant, at the same need again it is sufficiently small with guarantee will not extra-inhibitory convective term.Default value γ=ε2Usually
It is applicable initial valuation, needs to obtain the correct pressure change in interface in this model therefore use higher mobility.
Claims (6)
1. a kind of quasi- phase based on fracture hole mode seeps calculation method comprising following steps:
S1 collects the basic data of selected typical fracture hole mode and establishes two dimensional model；
S2 obtains different time and walks upper flow parameter according to Navier-Stokes equation simulation water drive oil process；
S3 seeks intending mutually seeping initial value under the mode using meteod of instability according to analog result；
This is mutually seeped the business simulation softward that initial value substitutes into based on Darcy formula and is fitted by S4, obtains relatively accurate seam
Hole mode is quasi- mutually to seep；
S5 converts fracture hole pattern configurations parameter, repeats step S2-S4, and the quasi- phase obtained under the fracture hole mode different shape is seeped.
2. a kind of quasi- phase based on fracture hole mode according to claim 1 seeps calculation method, it is characterised in that: the step
The concrete operations of S1 are: according to oil field actual conditions, extracting Reasonable Injection and adopt model, map that two-dimensional surface, construct model.
3. a kind of quasi- phase based on fracture hole mode according to claim 1 seeps calculation method, it is characterised in that: the step
The concrete operations of S2 are: using Navier-Stokes equation analysis flow regime, use Level Set Method processing oil-water interfaces
Power problem spreads interface problem using Cahn-Hilliard equation control water-oil phase.
4. a kind of quasi- phase based on fracture hole mode according to claim 1 seeps calculation method, it is characterised in that: the step
The concrete operations of S3 are: model being equivalent to a rock core, according to equisaturation mobile equation and B-L equation, calculates and corresponds to
Water phase relative permeability and oil relative permeability under saturation degree.
5. a kind of quasi- phase based on fracture hole mode according to claim 1 seeps calculation method, it is characterised in that: the step
The concrete operations of S4 are: the mutually infiltration initial value that step S3 is obtained substitute into operation in business software, analog parameter and fitting parameter with
Step S2 operational parameter is consistent, and fitting Production development result obtains the quasi- mutually infiltration of relatively accurate fracture hole mode.
6. a kind of quasi- phase based on fracture hole mode according to claim 1 seeps calculation method, it is characterised in that: the step
The concrete operations of S5 are: transformation fracture aperture, and the parameters such as solution cavity size proportion repeat step S2-S4, and calculating should under different shape
Fracture hole mode is quasi- mutually to seep.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910042454.8A CN109933747A (en) | 2019-01-17 | 2019-01-17 | A kind of quasi- phase infiltration calculation method based on fracture hole mode |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910042454.8A CN109933747A (en) | 2019-01-17 | 2019-01-17 | A kind of quasi- phase infiltration calculation method based on fracture hole mode |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109933747A true CN109933747A (en) | 2019-06-25 |
Family
ID=66985131
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910042454.8A Pending CN109933747A (en) | 2019-01-17 | 2019-01-17 | A kind of quasi- phase infiltration calculation method based on fracture hole mode |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109933747A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110714757A (en) * | 2019-11-26 | 2020-01-21 | 中国石油大学(华东) | Mechanical accuracy measurement method for residual oil of gas-drive starting blind end seam hole body |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102339326A (en) * | 2010-07-16 | 2012-02-01 | 中国石油化工股份有限公司 | Method for analyzing and simulating fluid flow of fracture-cavity oil reservoir |
| CN106951641A (en) * | 2017-03-22 | 2017-07-14 | 中国石油大学(华东) | A kind of method and system of fracture-pore reservoir numerical simulation |
-
2019
- 2019-01-17 CN CN201910042454.8A patent/CN109933747A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102339326A (en) * | 2010-07-16 | 2012-02-01 | 中国石油化工股份有限公司 | Method for analyzing and simulating fluid flow of fracture-cavity oil reservoir |
| CN106951641A (en) * | 2017-03-22 | 2017-07-14 | 中国石油大学(华东) | A kind of method and system of fracture-pore reservoir numerical simulation |
Non-Patent Citations (4)
| Title |
|---|
| AIMIN LV，等: "Characteristics of Oil-Water Relative Permeability and Influence Mechanism in Fractured-Vuggy Medium", 《THE SECOND SREE CONFERENCE ON CHEMICAL ENGINEERING》 * |
| HAOJUN XIE，等: "Coupling of two-phase flow in fractured-vuggy reservoir with fi lling medium", 《DE GRUYTER ACESS》 * |
| 吕爱民，等: "基于Buckley－Leverett理论的典型缝洞油藏一维驱替开发指标计算方法", 《西安石油大学学报(自然科学版)》 * |
| 赵健男，等: "碳酸盐岩缝洞型油藏均衡驱替机理研究", 《IFEDC 2018》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110714757A (en) * | 2019-11-26 | 2020-01-21 | 中国石油大学(华东) | Mechanical accuracy measurement method for residual oil of gas-drive starting blind end seam hole body |
| CN110714757B (en) * | 2019-11-26 | 2023-03-21 | 中国石油大学(华东) | Mechanical accuracy measurement method for residual oil of gas-drive starting blind end seam hole body |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gilman et al. | Improvements in simulation of naturally fractured reservoirs | |
| Qin et al. | A dynamic pore-network model for spontaneous imbibition in porous media | |
| Pariseau et al. | A new model for effects of impersistent joint sets on rock slope stability | |
| Pancharoen et al. | Modeling foam displacement in fractures | |
| CN102915406A (en) | Calculation method for relative permeability curve of oil and water in radial flow condition | |
| Yan et al. | An enhanced discrete fracture network model for multiphase flow in fractured reservoirs | |
| Ewing et al. | Mixed finite element approximation of phase velocities in compositional reservoir simulation | |
| CN113836695B (en) | Oil reservoir numerical simulation method based on gridless connecting element | |
| Cao | Development of a two-phase flow coupled capacitance resistance model | |
| Toll et al. | Soil–atmosphere interactions for analysing slopes in tropical soils in Singapore | |
| CN109958413A (en) | A kind of ultra-high water cut stage Reservoir behavior flow unit division methods | |
| CN107832515B (en) | Method and device for simulating coupling of oil reservoir and shaft | |
| CN109933747A (en) | A kind of quasi- phase infiltration calculation method based on fracture hole mode | |
| Gassara et al. | Equivalence between semi-empirical and population-balance foam models | |
| Zhang et al. | High-order streamline simulation and macro-scale visualization experimental studies on waterflooding under given pressure boundaries | |
| Douglas et al. | A data assimilation enabled model for coupling dual porosity flow with free flow | |
| Koçyigit et al. | Three-dimensional numerical modelling of wind-driven circulation in a homogeneous lake | |
| Zhao et al. | A comprehensive model for calculating relative permeability based on spontaneous imbibition and CT scanning measurement | |
| Zhang et al. | A multiscale mixed finite element method with oversampling for modeling flow in fractured reservoirs using discrete fracture model | |
| Muir Wood et al. | Deposition and self-weight consolidation of a shrinking fill | |
| Fuks et al. | Analysis of travel time distributions for uncertainty propagation in channelized porous systems | |
| Lewis et al. | A finite element analysis of multiphase immiscible flow in deforming porous media for subsurface systems | |
| CN108049861A (en) | Suitable for the bearing calibration of middle and high infiltration oil reservoir sealed coring well fluid saturation | |
| CN110516295A (en) | A kind of calculating CO2Drive the method for numerical simulation of minimum miscibility pressure | |
| CN107169227A (en) | The coarse grid analogy method and system of a kind of staged 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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190625 |
|
| WD01 | Invention patent application deemed withdrawn after publication |