CN105913494A - Multi-scale fracture fine geological modeling and value simulation method and device - Google Patents
Multi-scale fracture fine geological modeling and value simulation method and device Download PDFInfo
- Publication number
- CN105913494A CN105913494A CN201610192042.9A CN201610192042A CN105913494A CN 105913494 A CN105913494 A CN 105913494A CN 201610192042 A CN201610192042 A CN 201610192042A CN 105913494 A CN105913494 A CN 105913494A
- Authority
- CN
- China
- Prior art keywords
- level
- grid
- crack
- model
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a multi-scale fracture fine geological modeling and value simulation method and a device and belongs to the technical field of oil reservoir value simulation. The method comprises steps that a corner-point grid and the fracture network data are acquired; fractures are screened level by level; an unstructured basic grid and an unstructured value simulation grid are generated, matrix attribute mapping of the corner-point grid is carried out; III-level fractures are equivalent to the basic grid and are further equivalent to the value simulation grid; all conductivities of the final value simulation grid are calculated by utilizing the global coarsening technology. According to the method, fractures with different properties can be processed with pertinence level by level, on the condition that value simulation precision is guaranteed, the value simulation grid quantity is reduced, and the value simulation calculation speed is improved.
Description
Technical field
The present invention relates to reservoir description and Research Numerical Simulation Techique field, particularly relate to a kind of many chis
Spend crack fine geology modeling and method for numerical simulation and device.
Background technology
Fractured Reservoir is a kind of Oil Reservoir Types being widely present, and according to incompletely statistics, the most verifies ground
In matter reserves, fractured reservoir account for more than 28%.The modeling of Fractured Reservoir and digital-to-analogue are the most all
It is the big Major Difficulties of in the industry, is mainly reflected in the substrate crack with different stage at yardstick, seepage flow
On the greatest differences of conducting power and other flow behaviors.
For the characteristic in crack, primarily now there are double medium model and two kinds of methods of discrete fractures model
Carry out fracture to be described, and provide model and parameter for follow-up numerical simulation.Both models are respectively
There is different features and adaptability: double medium model advantage is stability and high efficiency, but precision is relatively low,
Often it is only applicable to the region that small fractures is grown, is not suitable for that the flowing of full oil reservoir is had significant impact
The situation in large scale crack；And the feature of discrete fractures model is fracture portrays finely, numerical simulation
Precision is high, and causing yet with high resolution cannot the oil reservoir grown of efficient simulation small fractures.
Fig. 1 i.e. shows a Fractured Reservoir (left), and correspondence double medium model (in) with from
Dissipate the schematic diagram of fractured model (right).
But, owing to fracture reservoir is usually present the crack of different origins, yardstick and the regularity of distribution,
Cause above two model all cannot reach the requirement in precision and efficiency simultaneously, it is therefore necessary to for
A kind of new modeling of multi-scale facture oil reservoir development and numerical simulation technology.
Due to the existence in little yardstick crack, for for the modeling of multi-scale facture oil reservoir and numerical simulation,
Still use double medium model outside Present Domestic more.That is, yardstick or the flow behavior in crack are not differentiated between,
The Fracture System unification of full oil reservoir is calculated numerical simulator parameter by specific method: crack
Interporosity flow coefficient between porosity and permeability, and substrate and crack grid, also referred to as form factor.
Wherein, the calculation of double medium model parameter has a lot, has traditional based on crack geometry
The computational methods of statistics, have the computational methods calculated based on local flow, and also one the most more
Popular method of roughening based on flowing.No matter use which kind of method, the most inevitably exist dual
The limitation of dielectric model.Such as, flowing method of roughening is more grown at fracture or presents netted point
When cloth, it is possible to obtain the double medium model that precision is higher, but be sporadicly distributed in crack or
The active region of some large fracture of person, it is likely that cause bigger error.Fig. 2 i.e. illustrates
Applicable application double medium model (left) and be suitable for the situation of application discrete fractures model (right).
Therefore, in prior art, double medium model technology has a disadvantage in that
1) different scale crack separately cannot be processed；
2) distribution fragmentary for crack or the active region of some large fracture, precision is relatively low.
In sum, the existence of multi-scale facture, result in the simulation of discrete fractures model value and calculate effect
Cannot adapt in rate, be typically only capable to use double medium model to be modeled and numerical simulation.But by
In double medium model for the most coarse description in crack and equivalent process, precision can be had a greatly reduced quality,
For the situation of some major fracture control area crude oil flow it is possible that bigger error.
Summary of the invention
The technical problem to be solved in the present invention is to provide one can take into account computational accuracy and computational efficiency
The modeling of multi-scale facture fine geology and method for numerical simulation and device.
For solving above-mentioned technical problem, the present invention provides technical scheme as follows:
On the one hand, it is provided that a kind of multi-scale facture fine geology modeling and method for numerical simulation, including:
Step 1: loading data, described data include Corner-point Grids model, host properties model, split
Seam network model and crack attribute model；
Step 2: according to fracture parameters, is divided into tri-ranks of I, II and III by crack, is respectively
Discrete fractures level, double porosity media level and EFFECTIVE MEDIUM level, wherein, the fracture scale maximum that I level is corresponding,
Fracture scale corresponding to II level is medium, and fracture scale corresponding to III level is minimum；
Step 3: according to Corner-point Grids model and I level and II level Fracture System geological information, generate I
Level and the unstructured grid G of II level discrete fracturesf, and the host properties of Corner-point Grids is mapped wherein,
Based on grid；
Step 4: according to Corner-point Grids model and I level Fracture System geological information, generate I level discrete
The unstructured grid G in crackc, as digital-to-analogue grid；
Step 5: by III level crack equivalence to GfSubstrate grid in, calculate its equivalent porosity and oozing
Rate thoroughly；
Step 6: II level crack is corresponded to digital-to-analogue grid GcIn, become splitting in double medium model
Seam medium grid；According to overall situation method of roughening, calculate GcThe parameters for numerical simulation of grid, the knot obtained
Fruit is discrete fractures+dual media mixed model.
On the other hand, it is provided that a kind of multi-scale facture fine geology modeling and numerical simulator obtain dress
Put, including:
Load-on module: for loading data, described data include Corner-point Grids model, host properties mould
Type, fracture network model of analogous outcrop and crack attribute model；
Crack diversity module: for according to fracture parameters, crack is divided into tri-levels of I, II and III
Not, respectively discrete fractures level, double porosity media level and EFFECTIVE MEDIUM level, wherein, corresponding the splitting of I level
Seam yardstick is maximum, and fracture scale corresponding to II level is medium, and fracture scale corresponding to III level is minimum；
First mess generation module: for several with I level and II level Fracture System according to Corner-point Grids model
What information, generates I level and the unstructured grid G of II level discrete fracturesf, and by the base of Corner-point Grids
Matter attribute map wherein, based on grid；
Second mess generation module: according to Corner-point Grids model and I level Fracture System geological information, raw
Become the unstructured grid G of I level discrete fracturesc, as digital-to-analogue grid；
First Equivalent Calculation module: for by III level crack equivalence to GfSubstrate grid in, calculate
Its equivalent porosity and permeability；
Second Equivalent Calculation module: for II level crack is corresponded to digital-to-analogue grid GcIn, become dual
Fracture medium grid in dielectric model；According to overall situation method of roughening, calculate GcThe numerical simulation of grid
Parameter, the result obtained is discrete fractures+dual media mixed model.
The method have the advantages that
In such scheme, the crack differing size has carried out sizing screening process.
(1) little yardstick (III level) crack is made EFFECTIVE MEDIUM process, the numerical simulation after equivalence
In model, III level crack is without additional values simulation lattice, it is to avoid, enormous amount less to yardstick
The crack that data uncertainty is bigger simultaneously carries out high-resolution modeling and model, greatly reduces meter
Calculate cost, numerical computations efficiency can be increased substantially on the premise of losing few precision.
(2) use dual media to process in mesoscale (II level) crack, use overall situation roughening flow process,
Can ensure that the equivalent precision in II crack, in final numerical simulator, II level crack is by dual Jie
Matter grid characterizes, it is only necessary to increase little numerical simulation grid.
(3) large scale (I level) crack being carried out discrete fractures process, it is right accurately to describe and simulate
The crack that oil reservoir plays an important role, maintains the high accuracy of conventional discrete fractured model.
To sum up, the inventive method uses separately process for the crack of different scale and reliability, has pin
Employing difference numerical simulation strategy to property, it is ensured that numerical stability, improves numerical computations
Efficiency.
Accompanying drawing explanation
Fig. 1 is Fractured Reservoir (left) in prior art, and correspondence double medium model (in)
Schematic diagram with discrete fractures model (right)；
Fig. 2 is to be suitable for application double medium model (left) in prior art and be suitable for application discrete fractures
The schematic diagram of the Fractured Reservoir of model (right)；
Fig. 3 is to obtain unstructured grid G in the present inventionfBased on the schematic diagram of grid；
Fig. 4 is to obtain unstructured grid G in the present inventioncSchematic diagram as digital-to-analogue grid；
The porosity attribute schematic diagram that Fig. 5 obtains when being and read Corner-point Grids and attribute thereof in the present invention；
Fig. 6 is according to the fracture network geometric data read in the present invention, carries out the signal of crack classification
Figure, wherein upper is III level crack, and lower-left is II level crack, and bottom right is: I level crack；
Fig. 7 is to generate model according to the Corner-point Grids information read in the step 3 of the present invention and step 4
The schematic diagram on border；
Fig. 8 is schematic diagram based on I+II level discrete fractures subdivision unstructured grid Gf in the present invention；
Fig. 9 be in the present invention from Corner-point Grids to the schematic diagram of Gf mesh mapping host properties, wherein
It it is porosity attribute shown in figure；
Figure 10 is schematic diagram based on I level discrete fractures subdivision unstructured grid Gc in the present invention；
Figure 11 is the schematic diagram of III level fracture parameter calculation, and the rightest figure is that the equivalence of III level crack is oozed
Rate thoroughly；
Figure 12 is the schematic diagram of mixed model roughening (obtaining the parameters such as conductivity)；
Figure 13 is the test example using the inventive method, and upper figure is contrast model numerical simulation
Result saturation distribution figure (the numerical simulation calculation time is 1619 seconds), figure below is the mixing of the present invention
Grid values analog result saturation distribution figure (the numerical simulation calculation time is 88 seconds), both errors
It is only 5%, and calculates the most original 1/18；
Figure 14 is the flow process principle schematic of the inventive method.
Detailed description of the invention
For making the technical problem to be solved in the present invention, technical scheme and advantage clearer, below will knot
Conjunction drawings and the specific embodiments are described in detail.
Generally reservoir also exists the crack that size is different with character, such as:
(1) pressure break major fracture general size is big, flow conductivity is strong, convection cell seepage flow and reservoir exploitation
Impact big, from distribution from the point of view of, this kind of crack negligible amounts, present the state of sparse distribution.By
In wall scroll crack, to seepage effect, greatly and distribution is the most sparse, for this kind of crack, uses traditional double
Dense media model has bigger error use, uses discrete fractures model more to have the most effectively.
(2) some crack general size such as pressure-break net and intrinsic fracture seam net are medium, flow conductivity
The strongest, distribution is interspersed in netted, applicable double medium model is simulated, according to discrete
Fractured model, can cause number of grid too much thus to calculate the problem that speed is the slowest.
(3) there is also some microcracks in reservoir, these fractue spacings are the most intensive, yet with
Flaw size is the least, might not be formed and effectively stitch net, therefore this kind of crack had both been not suitable for discrete splitting
Slit die type is also also not suitable for using double medium model to be simulated, and its attribute can be passed through ad hoc approach,
It is imparted in substrate.
As can be seen here, single double medium model or discrete fractures model is used effectively to simulate
Different types of crack in reservoir, present invention employs the thinking of three grades of crack classifications, it is simply that by these
The crack that character is different carries out sizing screening, and then uses distinct methods to set up model, and finally unification is arrived
One mixed model carries out numerical simulation.
On the one hand, the present invention provides a kind of multi-scale facture fine geology modeling and method for numerical simulation,
As shown in figs. 3-14, comprise the following steps:
Step 1: loading (reading) data, these data include Corner-point Grids model, host properties mould
Type, fracture network model of analogous outcrop and crack attribute model；
For the source of each data, it is respectively described below:
(1) Corner-point Grids model
Corner-point Grids model typically by three-dimensional geological modeling Software Create, and can export to Eclipse angle
Dot grid file (* .GRDECL) formatted file, conventional three-dimensional geological modeling software have Petrel,
SKUA/GOCAD, RMS etc..
(2) host properties model
Host properties model is set up after Corner-point Grids model generation, generally also soft at three-dimensional geological modeling
Part completes, it is possible to complete by specialty geology statistical softwares such as GsLib.Set up every attribute
Model is required for hard data (hard data) conduct input data, and hard data can be divided into according to its source
Core data, log data, geological data etc., under the constraint of hard data, by specific attribute
Modeling method sets up the attribute model that grid is corresponding, it may be assumed that give attribute for each grid units.Often
Attribute modeling method have Method of Stochastic (such as Gauss sequential simulation method), interpolation method (as
Common kriging method) and evaluation of expression method etc..This host properties model is based on angle in (1)
Dot grid model, including porosity, permeability and three attributes of net-gross ratio of substrate.
(3) fracture network model of analogous outcrop
Fracture network model of analogous outcrop is typically set up by specialty crack modeling software, and common crack modeling software has
Fraca, FracMan etc., going deep into and promoting along with fractured model research, many three-dimensional geological modelings
Software (such as Petrel, GOCAD etc.) also supports that fracture network models.Derive after modeling and preserve
Surface file (* .ts) form for GOCAD.
(4) crack attribute model
Crack attribute model is generally set up together with fracture network model of analogous outcrop, can be by setting up fracture network model of analogous outcrop
Software automatically calculate, it is possible to be manually entered.Generally comprise the aperture in crack, permeability.
Step 2: crack classification
According to fracture parameters, crack is divided into I, II and III tri-ranks, respectively discrete fractures
Level, double porosity media level and EFFECTIVE MEDIUM level, wherein, the fracture scale maximum that I level is corresponding, II level is right
The fracture scale answered is medium, the fracture scale minimum that III level is corresponding；
Specifically, can be according to the yardstick (such as length, height and aperture) in crack, the leading of crack
The data such as the distribution situation in stream ability and crack, are divided into three ranks by crack:
(1) discrete fractures level (I level)
This grade of fracture scale maximum, fracture condudtiviy is strong, and wall scroll crack convection cell seepage effect is maximum,
Final digital-to-analogue model exists with discrete fractures；
(2) double porosity media level (II level)
This grade of fracture scale is medium, and connection situation in crack at the same level is preferable, wall scroll crack convection cell seepage flow shadow
Ringing medium, in final numerical simulator, the fracture medium with double medium model exists；
(3) EFFECTIVE MEDIUM level (III level)
This grade of fracture scale minimum, wall scroll crack convection cell seepage effect is minimum, at final digital-to-analogue model
In, its to the effect of seepage flow by equivalence to substrate grid.
Concrete classification step can be:
(1) I level crack
According to I level-II level crack set in advance value parameter, including critical length L1, critical angle
d1, Critical permeability k1, then in fracture network F, subset
FI=f ∈ F | Lf>L1And df>d1And kf>k1It it is I level crack；
(2) III level crack
According to II level set in advance-III level crack value parameter, including critical length L2, critical open
Degree d2, Critical permeability k2, then in fracture network F, subset
FIII=f ∈ F | Lf<L2Or df<d2Or kf<k2It it is III level crack；
(3) II level crack
Subset F of fracture network FII=F-FI-FIIIFor II level crack.
Step 3: generate basic grid Gf
According to Corner-point Grids model and I level and II level Fracture System geological information, generate I level and II level
The unstructured grid G of discrete fracturesf, and the host properties of Corner-point Grids is mapped wherein, as base
Plinth grid；
In this step, according to Corner-point Grids model and I+II level Fracture System geological information, subdivision I+II
The unstructured grid G of level discrete fracturesfGrid (as shown in Figure 3) based on,.The most permissible
Including following sub-step:
(1) model boundary is generated
According to input angle dot grid, generate model external envelope face, as unstructured grid GfBorder；
(2) subdivision unstructured grid Gf
Mesh generation work can be completed by Grid Generation Softwares such as Triangle, TetGen or CGAL；
(3) host properties is mapped
For GfIn arbitrary grid g, find Corner-point Grids corresponding to grid according to its center point coordinate,
Grid g is given one by one by the host properties of this Corner-point Grids.
This grid has possessed all parameters of numerical simulation, it practice, this basic grid model passes exactly
System discrete fractures model, compared with the hybrid grid model that this patent proposes, this basic grid model net
Lattice number is many, and numerical simulation calculation cost is higher.Therefore in the art of this patent, this basic grid model
It is only used as a basis of follow-up work, is not directly used for numerical simulation calculation.In follow-up work, III
Level fracture parameter calculation and mixed model roughening step are all carried out based on this grid.
Step 4: generate digital-to-analogue grid Gc
According to Corner-point Grids model and I level Fracture System geological information, generate the non-of I level discrete fractures
Structured grid Gc, as digital-to-analogue grid；
This step is similar with step 3, and difference is the most only I level crack to be carried out unstructured grid
Subdivision, and without carrying out the mapping (carrying out in mixed model roughening step) of host properties.
This step obtains I level discrete fractures unstructured grid Gc(as shown in Figure 4), as finally
Digital-to-analogue grid.Owing to final numerical simulator system of parameters is obtained by hybrid grid roughening, therefore this step
Without carrying out attribute mapping in Zhou.
Concrete, this step can include following sub-step:
(1) model boundary is generated
According to input angle dot grid, generate model external envelope face, as unstructured grid GcBorder；
(2) subdivision unstructured grid Gc
Mesh generation work can be completed by Grid Generation Softwares such as Triangle, TetGen or CGAL.
This grid is as the carrier of final hybrid grid, it may be assumed that the numerical simulator finally given is base
On this grid.
Step 5:III level fracture parameter calculation
By III level crack equivalence to GfSubstrate grid in, calculate its equivalent porosity and permeability；
In the step for of, based on basic grid G in step 3f, by III level crack equivalence to Gf's
In substrate grid.Why by III level fracture parameters equivalence to basic grid GfRather than digital-to-analogue grid Gc,
It is because the former and has possessed all Numerical-Mode information in addition to III level crack, and the latter also lacks
II level crack and other host properties information, therefore by III level crack Equivalent Calculation to basic grid Gf,
Final hybrid grid roughening step just can obtain the hybrid grid model possessing complete information.
In III level fracture parameter calculation, porosity can use Traditional calculating methods equivalence (seeing below formula),
Permeability preferably employs Local Coarsening method based on flowing equivalence.I.e. for GfIn arbitrary grid g:
(1) traversal FIII, determine all III level cracks fallen within grid
(2) revising matrix porosity, computing formula is
Being wherein porosity for φ, V is volume；
(3) revising matrix permeability, computing formula is
Wherein Q is the total flow being perpendicular to permeability direction net cross-section, and d is on permeability direction
Gridding length, A is for being perpendicular to permeability direction net cross-section average area, and Δ P is barometric gradient
Component on permeability direction.
Owing to this rank fracture scale is little and uncertain high, therefore it is made without that the overall situation is discrete to be split
Seam characterizes, only need to be at unstructured grid GfEquivalent Calculation is carried out in single grid block.
Equivalent parameters calculates complete, then III level crack can be rejected from fracture network and be cast out.
Step 6: mixed model is roughened
II level crack is corresponded to digital-to-analogue grid GcIn, become the fracture medium net in double medium model
Lattice；According to overall situation method of roughening, calculate GcThe parameters for numerical simulation of grid, the result obtained is discrete
Crack+dual media mixed model.
Abovementioned steps has generated a set of basic grid G possessing complete digital-to-analogue informationf, its shortcoming is net
Lattice number is numerous, and the calculation cost of numerical simulation is the highest.Additionally, have also obtained final number in step 4
The unstructured grid G that value simulation usesc, digital-to-analogue grid GcG relativelyfFor want much less, but mesh
Till before, it contains only the information in I level crack.Therefore, the target of this step is exactly by basic grid
GfInformation (II level crack and other host properties), by method of roughening, be imparted to digital-to-analogue grid Gc
In.
In order to realize above-mentioned target, the step for use based on flowing overall method of roughening, obtain
All parameters of whole numerical simulator, work is divided into two parts, and one is to be corresponded in II level crack
GcIn grid, becoming the fracture medium grid in double medium model, two is to calculate GcThe numerical value of grid
Analog parameter, the result obtained is discrete fractures+dual media mixed model.Concrete steps are such as:
(1) traversal FII, determine and fall at GcAll II level cracks within each grid；
(2) calculating fracture network checkerwork cell porosity, computing formula is
Being wherein porosity for φ, V is volume,
(3) according to overall situation method of roughening, conductivity T it is calculated respectivelyMM、TFFAnd TMF, this
In T be conductivity, M and F refers to substrate and crack grid respectively.Noting, F has I and II here
The crack of two ranks, therefore, conductivity can be further subdivided into TMM、
So far, obtain based on GcThe final parameters for numerical simulation of grid:
(1) porosity of numerical simulation grid；
(2) the control point degree of depth of numerical simulation grid；
(3) volume of numerical simulation grid；
(4) the conductivity UNICOM table of numerical simulation grid.
Figure 13 is the test example using the inventive method, and upper figure is contrast model numerical simulation
Result saturation distribution figure (the numerical simulation calculation time is 1619 seconds), figure below is the mixing of the present invention
Grid values analog result saturation distribution figure (the numerical simulation calculation time is 88 seconds), both errors
It is only 5%, and calculates the most original 1/18.
On the other hand, corresponding with above-mentioned method, it is fine that the present invention also provides for a kind of multi-scale facture
Geologic modeling and numerical simulation device, including:
Load-on module: for loading data, described data include Corner-point Grids model, host properties mould
Type, fracture network model of analogous outcrop and crack attribute model；
Crack diversity module: for according to fracture parameters, crack is divided into tri-levels of I, II and III
Not, respectively discrete fractures level, double porosity media level and EFFECTIVE MEDIUM level, wherein, corresponding the splitting of I level
Seam yardstick is maximum, and fracture scale corresponding to II level is medium, and fracture scale corresponding to III level is minimum；
First mess generation module: for several with I level and II level Fracture System according to Corner-point Grids model
What information, generates I level and the unstructured grid G of II level discrete fracturesf, and by the base of Corner-point Grids
Matter attribute map wherein, based on grid；
Second mess generation module: according to Corner-point Grids model and I level Fracture System geological information, raw
Become the unstructured grid G of I level discrete fracturesc, as digital-to-analogue grid；
First Equivalent Calculation module: for by III level crack equivalence to GfSubstrate grid in, calculate
Its equivalent porosity and permeability；
Second Equivalent Calculation module: for II level crack is corresponded to digital-to-analogue grid GcIn, become dual
Fracture medium grid in dielectric model；According to overall situation method of roughening, calculate GcThe numerical simulation of grid
Parameter, the result obtained is discrete fractures+dual media mixed model.
Further, described crack diversity module preferably includes:
First classification submodule: according to I level-II level crack set in advance value parameter, including critical
Length L1, critical angle d1, Critical permeability k1, then in fracture network F, subset FI=
{f∈F|Lf>L1And df>d1And kf>k1It it is I level crack；
Second classification submodule: according to II level set in advance-III level crack value parameter, including facing
Boundary's length L2, critical angle d2, Critical permeability k2, then in fracture network F, subset
FIII=f ∈ F | Lf<L2Or df<d2Or kf<k2It it is III level crack；
3rd classification submodule: subset F of fracture network FII=F-FI-FIIIFor II level crack.
Further, described first mess generation module preferably includes:
First model boundary generates submodule: for according to input angle dot grid, generate model external envelope
Face, as unstructured grid GfBorder；
First subdivision submodule: for subdivision unstructured grid Gf；
Mapping submodule: for for GfIn arbitrary grid g, according to its center point coordinate find net
The Corner-point Grids that lattice are corresponding, gives grid g one by one by the host properties of this Corner-point Grids.
Further, described second mess generation module preferably includes:
Second model boundary generates submodule: for according to input angle dot grid, generate model external envelope
Face, as unstructured grid GcBorder；
Second subdivision submodule: for subdivision unstructured grid Gc。
Further, in described first Equivalent Calculation module, during III level fracture parameter calculation, porosity
Can use Traditional calculating methods equivalence, permeability preferably employs Local Coarsening method etc. based on flowing
Effect.
The above is the preferred embodiment of the present invention, it is noted that general for the art
For logical technical staff, on the premise of without departing from principle of the present invention, it is also possible to make some changing
Entering and retouch, these improvements and modifications also should be regarded as protection scope of the present invention.
Claims (10)
1. a multi-scale facture fine geology modeling and method for numerical simulation, it is characterised in that bag
Include:
Step 1: loading data, described data include Corner-point Grids model, host properties model, split
Seam network model and crack attribute model；
Step 2: according to fracture parameters, is divided into tri-ranks of I, II and III by crack, is respectively
Discrete fractures level, double porosity media level and EFFECTIVE MEDIUM level, wherein, the fracture scale maximum that I level is corresponding,
Fracture scale corresponding to II level is medium, and fracture scale corresponding to III level is minimum；
Step 3: according to Corner-point Grids model and I level and II level Fracture System geological information, generate I
Level and the unstructured grid G of II level discrete fracturesf, and the host properties of Corner-point Grids is mapped wherein,
Based on grid；
Step 4: according to Corner-point Grids model and I level Fracture System geological information, generate I level discrete
The unstructured grid G in crackc, as digital-to-analogue grid；
Step 5: by III level crack equivalence to GfSubstrate grid in, calculate its equivalent porosity and oozing
Rate thoroughly；
Step 6: II level crack is corresponded to digital-to-analogue grid GcIn, become splitting in double medium model
Seam medium grid；According to overall situation method of roughening, calculate GcThe parameters for numerical simulation of grid, the knot obtained
Fruit is discrete fractures+dual media mixed model.
Multi-scale facture fine geology the most according to claim 1 modeling and method for numerical simulation,
It is characterized in that, described step 2 includes:
Step 21: according to I level-II level crack set in advance value parameter, including critical length L1、
Critical angle d1, Critical permeability k1, then in fracture network F, subset
FI=f ∈ F | Lf>L1And df>d1And kf>k1It it is I level crack；
Step 22: according to II level set in advance-III level crack value parameter, including critical length L2、
Critical angle d2, Critical permeability k2, then in fracture network F, subset
FIII=f ∈ F | Lf<L2Or df<d2Or kf<k2It it is III level crack；
Step 23: subset F of fracture network FII=F-FI-FIIIFor II level crack.
Multi-scale facture fine geology the most according to claim 2 modeling and method for numerical simulation,
It is characterized in that, described step 3 includes:
Step 31: according to input angle dot grid, generates model external envelope face, as unstructured grid
GfBorder；
Step 32: subdivision unstructured grid Gf；
Step 33: for GfIn arbitrary grid g, according to its center point coordinate find grid corresponding
Corner-point Grids, gives grid g one by one by the host properties of this Corner-point Grids.
Multi-scale facture fine geology the most according to claim 3 modeling and method for numerical simulation,
It is characterized in that, described step 4 includes:
Step 41: according to input angle dot grid, generates model external envelope face, as unstructured grid
GcBorder；
Step 42: subdivision unstructured grid Gc。
Multi-scale facture fine geology the most according to claim 4 modeling and method for numerical simulation,
It is characterized in that, in described step 5, during III level fracture parameter calculation, porosity uses traditional calculations
Method equivalence, permeability uses Local Coarsening method based on flowing equivalence.
6. a multi-scale facture fine geology modeling and numerical simulation device, it is characterised in that bag
Include:
Load-on module: for loading data, described data include Corner-point Grids model, host properties mould
Type, fracture network model of analogous outcrop and crack attribute model；
Crack diversity module: for according to fracture parameters, crack is divided into tri-levels of I, II and III
Not, respectively discrete fractures level, double porosity media level and EFFECTIVE MEDIUM level, wherein, corresponding the splitting of I level
Seam yardstick is maximum, and fracture scale corresponding to II level is medium, and fracture scale corresponding to III level is minimum；
First mess generation module: for several with I level and II level Fracture System according to Corner-point Grids model
What information, generates I level and the unstructured grid G of II level discrete fracturesf, and by the base of Corner-point Grids
Matter attribute map wherein, based on grid；
Second mess generation module: according to Corner-point Grids model and I level Fracture System geological information, raw
Become the unstructured grid G of I level discrete fracturesc, as digital-to-analogue grid；
First Equivalent Calculation module: for by III level crack equivalence to GfSubstrate grid in, calculate
Its equivalent porosity and permeability；
Second Equivalent Calculation module: for II level crack is corresponded to digital-to-analogue grid GcIn, become dual
Fracture medium grid in dielectric model；According to overall situation method of roughening, calculate GcThe numerical simulation of grid
Parameter, the result obtained is discrete fractures+dual media mixed model.
Multi-scale facture fine geology the most according to claim 6 modeling and numerical simulation device,
It is characterized in that, described crack diversity module includes:
First classification submodule: according to I level-II level crack set in advance value parameter, including critical
Length L1, critical angle d1, Critical permeability k1, then in fracture network F, subset FI=
{f∈F|Lf>L1And df>d1And kf>k1It it is I level crack；
Second classification submodule: according to II level set in advance-III level crack value parameter, including facing
Boundary's length L2, critical angle d2, Critical permeability k2, then in fracture network F, subset
FIII=f ∈ F | Lf<L2Or df<d2Or kf<k2It it is III level crack；
3rd classification submodule: subset F of fracture network FII=F-FI-FIIIFor II level crack.
Multi-scale facture fine geology the most according to claim 7 modeling and numerical simulation device,
It is characterized in that, described first mess generation module includes:
First model boundary generates submodule: for according to input angle dot grid, generate model external envelope
Face, as unstructured grid GfBorder；
First subdivision submodule: for subdivision unstructured grid Gf；
Mapping submodule: for for GfIn arbitrary grid g, according to its center point coordinate find net
The Corner-point Grids that lattice are corresponding, gives grid g one by one by the host properties of this Corner-point Grids.
Multi-scale facture fine geology the most according to claim 8 modeling and numerical simulation device,
It is characterized in that, described second mess generation module includes:
Second model boundary generates submodule: for according to input angle dot grid, generate model external envelope
Face, as unstructured grid GcBorder；
Second subdivision submodule: for subdivision unstructured grid Gc。
Multi-scale facture fine geology the most according to claim 9 modeling and numerical simulation device,
It is characterized in that, in described first Equivalent Calculation module, during III level fracture parameter calculation, porosity is adopted
With Traditional calculating methods equivalence, permeability uses Local Coarsening method based on flowing equivalence.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610192042.9A CN105913494B (en) | 2016-03-30 | 2016-03-30 | The modeling of multi-scale facture fine geology and method for numerical simulation and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610192042.9A CN105913494B (en) | 2016-03-30 | 2016-03-30 | The modeling of multi-scale facture fine geology and method for numerical simulation and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105913494A true CN105913494A (en) | 2016-08-31 |
| CN105913494B CN105913494B (en) | 2018-09-14 |
Family
ID=56745324
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610192042.9A Expired - Fee Related CN105913494B (en) | 2016-03-30 | 2016-03-30 | The modeling of multi-scale facture fine geology and method for numerical simulation and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105913494B (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106934185A (en) * | 2017-04-27 | 2017-07-07 | 中国石油大学(华东) | A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid |
| CN107064193A (en) * | 2017-02-28 | 2017-08-18 | 中国石油大学(华东) | A kind of slit formation compact reservoir microcrack graded quantitative characterizing method |
| CN107704685A (en) * | 2017-10-10 | 2018-02-16 | 中国石油大学(北京) | A kind of Meshing Method and device |
| CN107832482A (en) * | 2017-08-14 | 2018-03-23 | 中国石油化工股份有限公司 | Compact reservoir multi-scale facture network modelling and analogy method |
| CN108627882A (en) * | 2017-03-24 | 2018-10-09 | 中国石油化工股份有限公司 | A method of establishing parallel fracture equivalent medium mode |
| CN108710734A (en) * | 2018-05-04 | 2018-10-26 | 南京特雷西能源科技有限公司 | Method for numerical simulation and device based on mesh adaption encryption and coarsening technique |
| CN109558631A (en) * | 2018-10-24 | 2019-04-02 | 北京大学 | A kind of densification oil-gas reservoir multi-scale facture medium automatic history matching method |
| CN110069791A (en) * | 2018-01-22 | 2019-07-30 | 中国石油化工股份有限公司 | A kind of modified laboratory experiment analogy method of scale |
| CN110263434A (en) * | 2019-06-20 | 2019-09-20 | 中国石油大学(华东) | A kind of flow unit method for numerical simulation based on multiple dimensioned mixed finite element |
| CN111324947A (en) * | 2020-01-21 | 2020-06-23 | 中国石油化工股份有限公司 | Oil reservoir numerical simulation non-uniform angular point grid optimization method based on attribute distribution |
| CN111400932A (en) * | 2020-04-09 | 2020-07-10 | 上海勘测设计研究院有限公司 | Method and device for generating grid in water environment simulation, storage medium and terminal |
| CN111898251A (en) * | 2020-07-06 | 2020-11-06 | 中国海洋石油集团有限公司 | Envelope surface construction method based on geological prototype model |
| CN115345079A (en) * | 2022-10-20 | 2022-11-15 | 中科数智能源科技(深圳)有限公司 | Oil reservoir numerical simulation data processing method and system of 2.5D unstructured network |
| CN117390998A (en) * | 2023-12-08 | 2024-01-12 | 西安石油大学 | Local grid encryption method and device for embedded discrete fracture model |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130096889A1 (en) * | 2011-10-12 | 2013-04-18 | IFP Energies Nouvelles | Method for generating a fractured reservoir mesh with a limited number of nodes in the matrix medium |
| CN104730596A (en) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | Discrete fracture modeling method based on multiscale factor restraint |
| CN104750896A (en) * | 2013-12-31 | 2015-07-01 | 中国石油化工股份有限公司 | Numerical simulation method of fractured-vug carbonate reservoir |
| CN104933208A (en) * | 2014-02-21 | 2015-09-23 | 中国石油化工股份有限公司 | Different-dimension crack numerical reservoir simulation comprehensive processing method |
| CN105160134A (en) * | 2015-09-30 | 2015-12-16 | 中国石油天然气股份有限公司 | Mixed medium simulation method and device for oil-gas flow in multiple media of tight reservoir |
| CN105205273A (en) * | 2015-09-30 | 2015-12-30 | 中国石油天然气股份有限公司 | Method and device for simulating oil-gas relay flow in multiple media of tight reservoir |
| CN105260543A (en) * | 2015-10-19 | 2016-01-20 | 中国石油天然气股份有限公司 | Multi-medium oil-gas flow simulation method and device based on double-hole model |
-
2016
- 2016-03-30 CN CN201610192042.9A patent/CN105913494B/en not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130096889A1 (en) * | 2011-10-12 | 2013-04-18 | IFP Energies Nouvelles | Method for generating a fractured reservoir mesh with a limited number of nodes in the matrix medium |
| CN104750896A (en) * | 2013-12-31 | 2015-07-01 | 中国石油化工股份有限公司 | Numerical simulation method of fractured-vug carbonate reservoir |
| CN104933208A (en) * | 2014-02-21 | 2015-09-23 | 中国石油化工股份有限公司 | Different-dimension crack numerical reservoir simulation comprehensive processing method |
| CN104730596A (en) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | Discrete fracture modeling method based on multiscale factor restraint |
| CN105160134A (en) * | 2015-09-30 | 2015-12-16 | 中国石油天然气股份有限公司 | Mixed medium simulation method and device for oil-gas flow in multiple media of tight reservoir |
| CN105205273A (en) * | 2015-09-30 | 2015-12-30 | 中国石油天然气股份有限公司 | Method and device for simulating oil-gas relay flow in multiple media of tight reservoir |
| CN105260543A (en) * | 2015-10-19 | 2016-01-20 | 中国石油天然气股份有限公司 | Multi-medium oil-gas flow simulation method and device based on double-hole model |
Non-Patent Citations (3)
| Title |
|---|
| 张允 等: "复杂介质油藏多尺度数值模拟研究", 《中国力学大学-2015》 * |
| 肖阳 等: "缝洞型碳酸盐岩油气藏数值模拟技术研究", 《新疆石油天然气》 * |
| 郎晓玲 等: "基于DFN离散裂缝网络模型的裂缝性储层建模方法", 《北京大学学报（自然科学版）》 * |
Cited By (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107064193A (en) * | 2017-02-28 | 2017-08-18 | 中国石油大学(华东) | A kind of slit formation compact reservoir microcrack graded quantitative characterizing method |
| CN107064193B (en) * | 2017-02-28 | 2018-05-22 | 中国石油大学(华东) | A kind of slit formation compact reservoir microcrack graded quantitative characterizing method |
| CN108627882A (en) * | 2017-03-24 | 2018-10-09 | 中国石油化工股份有限公司 | A method of establishing parallel fracture equivalent medium mode |
| CN106934185A (en) * | 2017-04-27 | 2017-07-07 | 中国石油大学(华东) | A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid |
| CN106934185B (en) * | 2017-04-27 | 2019-08-06 | 中国石油大学(华东) | A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid |
| CN107832482B (en) * | 2017-08-14 | 2021-03-23 | 中国石油化工股份有限公司 | Compact reservoir multi-scale fracture network modeling and simulation method |
| CN107832482A (en) * | 2017-08-14 | 2018-03-23 | 中国石油化工股份有限公司 | Compact reservoir multi-scale facture network modelling and analogy method |
| CN107704685B (en) * | 2017-10-10 | 2020-04-17 | 中国石油大学(北京) | Mesh division method and device |
| CN107704685A (en) * | 2017-10-10 | 2018-02-16 | 中国石油大学(北京) | A kind of Meshing Method and device |
| CN110069791A (en) * | 2018-01-22 | 2019-07-30 | 中国石油化工股份有限公司 | A kind of modified laboratory experiment analogy method of scale |
| CN110069791B (en) * | 2018-01-22 | 2021-11-02 | 中国石油化工股份有限公司 | Indoor experiment simulation method for scale correction |
| CN108710734B (en) * | 2018-05-04 | 2022-12-20 | 特雷西能源科技(杭州)有限公司 | Numerical simulation method and device based on grid adaptive encryption and coarsening technology |
| CN108710734A (en) * | 2018-05-04 | 2018-10-26 | 南京特雷西能源科技有限公司 | Method for numerical simulation and device based on mesh adaption encryption and coarsening technique |
| CN109558631A (en) * | 2018-10-24 | 2019-04-02 | 北京大学 | A kind of densification oil-gas reservoir multi-scale facture medium automatic history matching method |
| CN110263434A (en) * | 2019-06-20 | 2019-09-20 | 中国石油大学(华东) | A kind of flow unit method for numerical simulation based on multiple dimensioned mixed finite element |
| CN111324947A (en) * | 2020-01-21 | 2020-06-23 | 中国石油化工股份有限公司 | Oil reservoir numerical simulation non-uniform angular point grid optimization method based on attribute distribution |
| CN111400932A (en) * | 2020-04-09 | 2020-07-10 | 上海勘测设计研究院有限公司 | Method and device for generating grid in water environment simulation, storage medium and terminal |
| CN111400932B (en) * | 2020-04-09 | 2022-09-16 | 上海勘测设计研究院有限公司 | Method and device for generating grid in water environment simulation, storage medium and terminal |
| CN111898251A (en) * | 2020-07-06 | 2020-11-06 | 中国海洋石油集团有限公司 | Envelope surface construction method based on geological prototype model |
| CN115345079A (en) * | 2022-10-20 | 2022-11-15 | 中科数智能源科技(深圳)有限公司 | Oil reservoir numerical simulation data processing method and system of 2.5D unstructured network |
| CN115345079B (en) * | 2022-10-20 | 2023-01-10 | 中科数智能源科技(深圳)有限公司 | Oil reservoir numerical simulation data processing method and system of 2.5D unstructured network |
| CN117390998A (en) * | 2023-12-08 | 2024-01-12 | 西安石油大学 | Local grid encryption method and device for embedded discrete fracture model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105913494B (en) | 2018-09-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105913494A (en) | Multi-scale fracture fine geological modeling and value simulation method and device | |
| CN110286421B (en) | Method for modeling natural fractures of tight sandstone reservoir | |
| US6078869A (en) | Method and apparatus for generating more accurate earth formation grid cell property information for use by a simulator to display more accurate simulation results of the formation near a wellbore | |
| US12061306B2 (en) | Constructing structural models of the subsurface | |
| CN106150477B (en) | A kind of method of the single well controlled reserves of determining fracture-pore reservoir | |
| CN106227957A (en) | The method of equivalent fissure modeling | |
| GB2563241A (en) | Hydraulic fracturing simulation | |
| US10309216B2 (en) | Method of upscaling a discrete fracture network model | |
| CN107194104B (en) | The method for calculating inanimate matter porosity in shale reservoir | |
| CN109856674A (en) | Engineering dessert evaluating method and device | |
| CN104933208B (en) | A kind of different scale Fractured Reservoir numerical simulation integrated conduct method | |
| CN103376463B (en) | A kind of Inverse modeling method based on faults control | |
| CN106886046A (en) | Method and device for determining available reserves of unproductive blocks of fracture-cavity gas reservoir | |
| CN114429525A (en) | Fractured reservoir multi-scale fracture simulation method and computer readable storage medium | |
| EP3020913A1 (en) | Method of managing petro-chemical reservoir production and program product therefor | |
| CN107507179B (en) | Rock-soil mass quantitative analysis method based on GOCAD | |
| CN106597542B (en) | Reservoir characteristic parameter prediction method and device | |
| AU2011360602B2 (en) | Computerized method for the estimation of a value for at least a parameter of a hydrocarbon-producing region, for planning the operation and operating the region | |
| CN114970235A (en) | Geomechanical heterogeneous-anisotropic modeling method for fractured reservoir | |
| CN108595782B (en) | Calculation method for mass transfer between matrix and cracks in discrete cracks | |
| Kumar et al. | Comparative analysis of dual continuum and discrete fracture simulation approaches to model fluid flow in naturally fractured, low-permeability reservoirs | |
| CN107832482A (en) | Compact reservoir multi-scale facture network modelling and analogy method | |
| CN115984503B (en) | Geological profile generation method, system, electronic equipment and medium | |
| CN105068141A (en) | Method and device for determining karst reservoir distribution of fracture system | |
| CN103425881A (en) | Method for certainty numerical simulation of crack medium seismic wave response |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180914 Termination date: 20190330 |