A kind of diaclase three-dimensional modeling method restored based on construction face geometry
Technical field
The present invention relates to geological technique fields more particularly to a kind of diaclase three-dimensional restored based on construction face geometry to build
Mould method.
Background technology
Natural formation crack is a kind of spy increasing substantially petro-gas permeation relatively conventional in all kinds of oil and gas reservoirs
Different reservoir medium, it is controlled by tectonic stress, strain variation during tectonic movement, in oil-gas reservoir part or greatly
The intrinsic fracture that subregion occurs.Natural formation crack has the phase time mixes, system is more, distributed more widely, feature is complicated etc.
Feature, characterization and understanding to natural formation crack always are one of the problem during production of hydrocarbons.
The method of natural formation crack characterization is more, there is field section and core analysis method, well log interpretation method, tectonic stress
And tectonic curvaceous, seismic description method, stochastic simulation method etc..Field section and core analysis method are the methods by core description
System, density, inclination angle, the orientation etc. in qualitative characterization's natural formation crack are that acquisition natural formation crack attributive character is most direct
Method.Well log interpretation method is the method that the method explained by imaging logging or conventional logging characterizes natural formation crack, energy
Completely realize the description to full near wellbore natural formation crack attributive character and characterization.Tectonic stress and tectonic curvaceous are
The probability difference in natural formation crack is generated based on Various Tectonic Stresses state or curvature induction and is realized to natural formation crack
The characterization of distribution.Seismic description method is realized using the seismic properties sensitive to natural formation crack to natural formation crack
Prediction.Stochastic simulation method realizes the prediction to natural formation crack using the method for stochastic simulation.
It is field section and core analysis method, well log interpretation method, tectonic stress and tectonic curvaceous, seismic description method, random
The natural formations such as simulation crack characterizing method has the advantage of oneself, also there is the deficiency of itself.Field section and core analysis
Method is acquisition natural formation crack attributive character most straightforward approach, but can only qualitative characterization.Well log interpretation method can be completely real
The now description to full near wellbore natural formation crack attributive character and characterization, but be difficult to realize inter-wells predicting.Tectonic stress and
The prediction of natural formation fractue spacing between well may be implemented in tectonic curvaceous, stochastic simulation method, but one is that can not obtain natural structure
Crack figure information is made, second is that the accuracy of prediction result is difficult to ensure.Earthquake (such as ant body) describes rule due to resolution ratio
The reason of can predict large-scale tomography, it is difficult to predict small natural formation crack, while can only also provide distributed intelligence, can not
Obtain occurrence information.
Invention content
It is an object of the invention to solve the problems of the above-mentioned prior art, provide a kind of based on the recovery of construction face geometry
Diaclase three-dimensional modeling method.
A kind of diaclase three-dimensional modeling method restored based on construction face geometry, is included the following steps:
Step 1:Establish three-dimensional construction FEM layer model;
Step 2:Establish the recovery of construction level geometry and diaclase occurrence index collection of illustrative plates;
Step 3:Establish diaclase occurrence index model;
Step 4:Rock core/explanation imaging logging and extraction diaclase characteristic parameter are described;
Step 5:Establish diaclase system pattern and diaclase system strength model;
Step 6:Establish diaclase discrete model;
Step 7:Establish diaclase equivalent model.
Further, method as described above, in step 2, the construction level geometry recovery is to three-dimensional in step 1
It constructs the deformation layer plane of FEM layer model and fault plane carries out de-shirredization and tomographyization is gone to handle respectively, eliminate later structural and live
The structural deformation of movable property life, restores the deformation behaviour of paleostructure activity early period of origination sedimentary formation, discloses paleostructure motion process
In stress situation occurred and stratum distribution of force situation.
Further, method as described above, in step 2, the diaclase occurrence index is several by construction level
What what recovery process was extracted, diaclase occurrence index be according to level restore in front and back each grid deformed region area with
The ratio between area of plane after recovery, i.e.,
DIi=(Ai-Afi)/Ai
Wherein, DIiFor diaclase occurrence index, i is grid serial number, and Af is distortion of the mesh amount, and A is grid area.
Further, method as described above, establishing diaclase occurrence index model described in step 3 includes:In step
A cutoff value is set on the two diaclase occurrence index collection of illustrative plates, and the region less than or equal to the cutoff value is set as 1, is more than
The region of the cutoff value is set as 0, you can establishes and forms the diaclase occurrence index mould that spatial domain range occurs for indication structure crack
Type.
Further, method as described above, it is according to build joint that diaclase system pattern is established described in step 5
Orientation and inclination angle divide crack system to establish diaclase system pattern;
The diaclase system strength model is to seek each system structure by the diaclase system pattern of establishing
The product of fracture spacing and width is made to build the intensive variable of build joint.
Further, method as described above, it is to be sent out with diaclase that diaclase discrete model is established described in step 6
Raw exponential model is that the spatial domain of diaclase interpolation constrains, and the diaclase system strength model by foundation is primary input, profit
It can be established to form diaclase discrete model with reservoir modeling software.
Further, method as described above, establishing diaclase equivalent model described in step 7 includes:
Physical property is oozed by test development build joint and the hole of the core sample of agensis build joint, is developed by seeking respectively
The average pore and mean permeability of build joint and the core sample of agensis build joint, and the flat of build joint rock sample will be developed
Ratio between porosity and the average pore and mean permeability of the core sample of mean permeability and agensis build joint
Relationship is equivalent to the proportionate relationship of Reservoir Tectonic Fracture porosity, permeability respectively between matrix porosity and permeability, is based on
The proportionate relationship can be established to form structure on the basis of reservoir matrix porosity, penetration rate model using reservoir modeling software
Make crack equivalent model.
The method of the present invention absorb natural formation crack characterization field section and core analysis method, well log interpretation method and with
On the basis of machine simulation advantage, restores by construction level geometry, diaclase occurrence index model is built, as spatial domain
Natural formation crack 3 d-dem model is established in constraint, realizes every group of natural formation fractue spacing position, development scale, extension side
The three-dimensional visualization of position and cutting relationship;Establish natural formation crack three dimensional equivalent model, every group of natural structure of quantitatively characterizing simultaneously
Make the porosity and permeability of reservoir at fractue spacing position.The modeling result of acquisition has not only reproduced natural structure in three dimensions
The distribution scale and feature in crack are made, also providing natural formation crack 3-D quantitative hole for production simulation oozes attribute model, from
And it is that solid foundation has been established in the scientific and reasonable formulation of production of hydrocarbons scheme.The main contents of the present invention include following content.
(1) three-dimensional construction FEM layer model
Combined method, which is shaken, by well identifies that stratum and fracture system, input commercialization modeling software establish three-dimensional construction level
Model is the mature technology in current oil and gas reservoir representational field, establishes the three-dimensional construction FEM layer model of completion and is needed as the present invention
The original input data wanted.
(2) construction level geometry restores and diaclase occurrence index collection of illustrative plates is established
Construction level geometry restores to be to carry out pleat respectively with fault plane to the deformation layer plane of three-dimensional construction FEM layer model
Wrinkleization and go tomographyization handle, eliminate later structural activity generate structural deformation, restore paleostructure activity early period of origination deposition
The deformation behaviour on stratum discloses stress situation occurred and stratum distribution of force situation in paleostructure motion process.According to simple shear
Angle principle, engagement formation analysis and balanced section technique, can realize that construction level geometry restores in high quality.
Diaclase occurrence index extracts to obtain by construction level geometry recovery process, diaclase occurrence index
It is to restore the ratio between deformed region area and the area of plane after recovery in front and back each grid according to level, i.e.,
DIi=(Ai-Afi)/Ai
Wherein DIiFor diaclase occurrence index, i is grid serial number, and Af is distortion of the mesh amount, and A is grid area.
It can be seen that diaclase occurrence index indicates the severe degree of structural deformation, value is bigger, and structural deformation is got over
Weak, value is smaller, and structural deformation is more violent;And the severe degree of structural deformation is directly proportional to diaclase odds, construction
Deformation is more violent, and diaclase odds is bigger.It tints by the size of diaclase occurrence index value on each grid,
Diaclase occurrence index value is smaller, and color is whiter, thus establishes and obtains diaclase occurrence index collection of illustrative plates;And color blackening
Degree then reflect diaclase generation strength information, black region represent build joint occur intensity it is larger, white area
It is smaller to represent build joint generation generation intensity.
(3) foundation of diaclase occurrence index model
Diaclase occurrence index collection of illustrative plates specifies the possibility of build joint generation, but can not limit the sky of build joint generation
Domain position, for this reason, it may be necessary to a cutoff value is set on diaclase occurrence index collection of illustrative plates, it will be less than or equal to the area of the cutoff value
Domain is set as 1, and the region more than the cutoff value is set as 0, you can the construction that foundation forms indication structure crack generation spatial domain range is split
Stitch occurrence index model.And cutoff value selection will be limited with that can reflect the crack that the overwhelming majority may occur, it can be with three-dimensional structure
It is selected after making the fracture system comparison on FEM layer model.
(4) core description/image logging interpretation and diaclase characteristic parameter extraction
Field section and core analysis method are to obtain the attributive character such as natural formation fracture spacing, inclination angle, orientation, width most
If direct method can completely identify the fractue spacing feature of individual well section in conjunction with image logging interpretation.
By core description/image logging interpretation, extract and in Revision area along the system of the diaclase of well profile
The attributive character such as feature, and density, width, inclination angle, the orientation of each system, to build the input data of diaclase modeling
Library.
(5) foundation of diaclase system pattern and strength model
Crack system is divided according to the orientation of build joint and inclination angle, seeks multiplying for each system diaclase density and width
It is soft to be constituted the input modeling of input data body by the intensive variable of product structure build joint for the individual well build joint intensive variable of each system
It is that constraints establishes each system structure using diaclase occurrence index model using suitable stochastic simulation algorithm in part
The strength model of seam characterizes the spatial distribution characteristic of each system build joint growth strength.
(6) foundation of diaclase discrete model
It is constrained by the spatial domain of diaclase interpolation of diaclase exponential model, passes through the diaclase strength model of foundation
It for primary input, can establish to form diaclase discrete model using ripe reservoir modeling business software, to accurately exist
Three dimensions reproduces the distributing position of each crack piece in the system feature in natural formation crack, and each system, orientation, inclines
Angle and form.
(7) foundation of diaclase equivalent model
Physical property is oozed in the hole that the different core sample of development degree is stitched by test configuration, and extraction respectively is obtained developing constructing and be split
Average pore, the mean permeability value of seam and agensis diaclase sample are constructed using development diaclase and agensis
Average pore, the mean permeability value of crack sample are quotient respectively, and using quotient as diaclase with matrix porosity and ooze
Proportionate relationship between saturating rate, maturation is utilized based on the proportionate relationship on the basis of reservoir matrix porosity, penetration rate model
Reservoir modeling business software can establish to form diaclase equivalent model, including equivalent porosity threedimensional model, I, J, K tri-
The equivalent permeability threedimensional model in a direction.
Description of the drawings
Fig. 1 is that the present invention is based on the diaclase three-dimensional modeling method flow charts that construction face geometry restores;
Fig. 2 a are five peaks of East Sichuan Basin of the embodiment of the present invention-No. 1 bedding plane structure crack of Longma small stream group shale gas reservoir
Occurrence index collection of illustrative plates;
Fig. 2 b are five peaks of East Sichuan Basin of the embodiment of the present invention-No. 5 bedding plane structure cracks of Longma small stream group shale gas reservoir
Occurrence index collection of illustrative plates;
Five peaks of Fig. 3 a East Sichuan Basins of the embodiment of the present invention-floor diaclase of Longma small stream group shale gas reservoir 1 occurs
Exponential model;
Five peaks of Fig. 3 b East Sichuan Basins of the embodiment of the present invention-floor diaclase of Longma small stream group shale gas reservoir 5 occurs
Exponential model;
Fig. 4 a are five peaks of East Sichuan Basin of the embodiment of the present invention-Longma small stream group shale gas reservoir core observation diaclase
System characteristic pattern;
Fig. 4 b are five peaks of East Sichuan Basin of the embodiment of the present invention-Longma small stream group shale gas reservoir diaclase system feature
Change individual well sectional view;
Fig. 5 a are five peaks of East Sichuan Basin of the embodiment of the present invention-floor A high angles of Longma small stream group shale gas reservoir 1 construction
Intensity illustraton of model;
Fig. 5 b are five peaks of East Sichuan Basin of the embodiment of the present invention-floor B low angles of Longma small stream group shale gas reservoir 1 construction
Intensity illustraton of model;
Fig. 6 a are five peaks of East Sichuan Basin of the embodiment of the present invention-Longma small stream group shale gas reservoir A high angle natural formations
Crack 3 d-dem illustraton of model;
Fig. 6 b are five peaks of East Sichuan Basin of the embodiment of the present invention-Longma small stream group shale gas reservoir B low angle natural formations
Crack 3 d-dem illustraton of model;
Fig. 7 a are five peaks of East Sichuan Basin of the embodiment of the present invention-equivalent hole of Longma small stream group shale gas reservoir natural formation
Spend three-dimensional model diagram;
Fig. 7 b are five peaks of East Sichuan Basin of the embodiment of the present invention-directions Longma small stream group shale gas reservoir natural formation I etc.
Imitate permeability three-dimensional model diagram;
Fig. 7 c are five peaks of East Sichuan Basin of the embodiment of the present invention-directions Longma small stream group shale gas reservoir natural formation J etc.
Imitate permeability three-dimensional model diagram;
Fig. 7 d are five peaks of East Sichuan Basin of the embodiment of the present invention-directions Longma small stream group shale gas reservoir natural formation K etc.
Imitate permeability three-dimensional model diagram.
Specific implementation mode
To make the object, technical solutions and advantages of the present invention clearer, the technical solution in the present invention is carried out below
It clearly and completely describes, it is clear that described embodiments are some of the embodiments of the present invention, instead of all the embodiments.Base
Embodiment in the present invention, those of ordinary skill in the art obtained without creative efforts it is all its
His embodiment, shall fall within the protection scope of the present invention.
The present invention is primarily based on the origin mechanism that natural formation fracture development is controlled by tectonic stress/strain field, in three-dimensional
On the basis of constructing FEM layer model, carry out construction level geometry and restore, extract diaclase occurrence index collection of illustrative plates, establishes construction and split
Occurrence index model is stitched, the spatial domain range of diaclase modeling is formed;Then pass through core description/imaging logging, extraction construction
FRACTURE CHARACTERISTICS parameter builds diaclase system pattern, establishes diaclase system strength model, discloses diaclase development
Spatial intensity distribution;It constrains for spatial domain according to diaclase occurrence index model, is built using diaclase system strength model again
It is vertical to form diaclase discrete model;Fracture porosity, the permeability information for finally utilizing core analysis to obtain, establish and complete structure
Make crack equivalent model.
Detailed technical solution flow diagram is as shown in Figure 1.
(1) three-dimensional construction FEM layer model
Combined method, which is shaken, by well identifies that stratum and fracture system, input commercialization modeling software establish three-dimensional construction level
Model is the mature technology in current oil and gas reservoir representational field, establishes the three-dimensional construction FEM layer model of completion and is needed as the present invention
The original input data wanted.
(2) construction level geometry restores and diaclase occurrence index collection of illustrative plates is established
Construction level geometry restores to be to carry out pleat respectively with fault plane to the deformation layer plane of three-dimensional construction FEM layer model
Wrinkleization and go tomographyization handle, eliminate later structural activity generate structural deformation, restore paleostructure activity early period of origination deposition
The deformation behaviour on stratum discloses stress situation occurred and stratum distribution of force situation in paleostructure motion process.According to simple shear
Angle principle, engagement formation analysis and balanced section technique, can realize that construction level geometry restores in high quality.
Diaclase occurrence index extracts to obtain by construction level geometry recovery process, diaclase occurrence index
It is to restore the ratio between deformed region area and the area of plane after recovery in front and back each grid according to level, i.e.,
DIi=(Ai-Afi)/Ai
Wherein DIiFor diaclase occurrence index, i is grid serial number, and Af is distortion of the mesh amount, and A is grid area.
It can be seen that diaclase occurrence index indicates the severe degree of structural deformation, value is bigger, and structural deformation is got over
Weak, value is smaller, and structural deformation is more violent;And the severe degree of structural deformation is directly proportional to diaclase odds, construction
Deformation is more violent, and diaclase odds is bigger.It tints by the size of diaclase occurrence index value on each grid,
Diaclase occurrence index value is smaller, and color is whiter, thus establishes and obtains diaclase occurrence index collection of illustrative plates;And color blackening
Degree then reflect diaclase generation strength information, black region represent build joint occur intensity it is larger, white area
It is smaller (Fig. 2 a, Fig. 2 b) to represent build joint generation generation intensity.
Fig. 2 a, 2b are set forth five peaks of East Sichuan Basin-Longma small stream group gas reservoir and construct No. 1 and No. 5 face warp layer by layer
Construction level geometry restores the diaclase occurrence index collection of illustrative plates that extraction obtains.The diaclase occurrence index of two-layer structure level
Collection of illustrative plates is in addition to locally there is nuance, and essential characteristic is almost the same, and the areas Ji Gai diaclase is prone in peripheral edge margin, in
Between region only part be sporadicly dispersed with diaclase.
(3) foundation of diaclase occurrence index model
Diaclase occurrence index collection of illustrative plates specifies the possibility of build joint generation, but can not limit the sky of build joint generation
Domain position, for this reason, it may be necessary to a cutoff value is set on diaclase occurrence index collection of illustrative plates, it will be less than or equal to the area of the cutoff value
Domain is set as 1, and the region more than the cutoff value is set as 0, you can the construction that foundation forms indication structure crack generation spatial domain range is split
Stitch occurrence index model (Fig. 3 a, Fig. 3 b).And cutoff value selection will be limited with that can reflect the crack that the overwhelming majority may occur, it can
It is selected after being compared with the fracture system constructed with three-dimensional on FEM layer model.Fig. 3 a, Fig. 3 b are that selection cut-off is 0.99968 acquisition
Five peaks of East Sichuan Basin-Longma small stream group shale gas reservoir diaclase occurrence index model, the model substantially cover in area
Possible diaclase development area.
(4) core description/image logging interpretation and diaclase characteristic parameter extraction
Field section and core analysis method are to obtain the attributive character such as natural formation fracture spacing, inclination angle, orientation, width most
If direct method can completely identify the fractue spacing feature of individual well section in conjunction with image logging interpretation.
By core description/image logging interpretation, extract and in Revision area along the system of the diaclase of well profile
The attributive character such as feature, and density, width, inclination angle, the orientation of each system, to build the input data of diaclase modeling
Library.
Fig. 4 a, Fig. 4 b are respectively to be identified by five peaks of East Sichuan Basin-Longma small stream group shale gas reservoir typical core
High angle, two groups of diaclases of low angle, respectively describe the attributive character such as density, width, inclination angle, the orientation of diaclase
And be labeled on individual well sectional view, form the knowledge base of research area's internal fissure modeling.
(5) foundation of diaclase system pattern and strength model
Crack system is divided according to the orientation of build joint and inclination angle, seeks multiplying for each system diaclase density and width
It is soft to be constituted the input modeling of input data body by the intensive variable of product structure build joint for the individual well build joint intensive variable of each system
It is that constraints establishes each system structure using diaclase occurrence index model using suitable stochastic simulation algorithm in part
The strength model of seam characterizes the spatial distribution characteristic of each system build joint growth strength.
Fig. 5 a, Fig. 5 b are to establish five peaks of the East Sichuan Basin-Longma small stream group shale completed using sequence Gauss method
No. 1 floor diaclase strength model of gas reservoir, which characterizes No. 1 floor height angle respectively and low angle amounts to two groups of diaclases
Intensity distribution feature.
(6) foundation of diaclase discrete model
It is constrained by the spatial domain of diaclase interpolation of diaclase exponential model, passes through the diaclase strength model of foundation
It for primary input, can establish to form diaclase discrete model using ripe reservoir modeling business software, to accurately exist
Three dimensions reproduces the distributing position of each crack piece in the system feature in natural formation crack, and each system, orientation, inclines
Angle and form.Fig. 6 a, Fig. 6 b are to establish five peaks of the East Sichuan Basin-No. 1 floor height angle of Longma small stream group gas reservoir and low obtained
Angle amounts to two groups of diaclase discrete models, characterizes the distributing position of high angle and low angle diaclase, orientation, inclination angle
And form.
(7) foundation of diaclase equivalent model
Physical property is oozed in the hole that the different core sample of development degree is stitched by test configuration, and extraction respectively is obtained developing constructing and be split
Average pore, the mean permeability value of seam and agensis diaclase sample are constructed using development diaclase and agensis
Average pore, the mean permeability value of crack sample are quotient respectively, and using quotient as diaclase with matrix porosity and ooze
Proportionate relationship between saturating rate is based on the proportionate relationship, be based on the proportionate relationship, in reservoir matrix porosity, penetration rate model
On the basis of, it can establish to form diaclase equivalent model, including equivalent porosity using ripe reservoir modeling business software
Threedimensional model, the equivalent permeability threedimensional model in tri- directions I, J, K.
Fig. 7 a- Fig. 7 d are to establish five peaks of the East Sichuan Basin-equivalent hole of Longma small stream group gas reservoir diaclase obtained
Threedimensional model is spent, the equivalent permeability threedimensional model in tri- directions I, J, K is realized to the gas reservoir three dimensions build joint hole
The quantitatively characterizing of tri- directional structure vectorical structure seam permeabilities of degree and I, J, K.
What the present invention solved natural formation crack 3 d-dem model and equivalent model establishes problem, in Fractured oil gas
To be obtained in reservoir characterization it is extensive promote and apply, at present at China's East Sichuan Basin First Five-Year Plan peak-Longma small stream group gas reservoir
It is applied during the foundation of natural formation crack 3-D quantitative geological model, brings good economic results in society.
The present invention carries out construction level geometry and restores on the basis of three-dimensional constructs FEM layer model, extraction diaclase hair
Raw index collection of illustrative plates, establishes diaclase occurrence index model, forms the spatial domain range of diaclase modeling;Then it is retouched by rock core
/ imaging logging is stated, diaclase characteristic parameter is extracted, diaclase system pattern is built, establishes diaclase system intensity mould
Type discloses the spatial intensity distribution of diaclase development;It is again space constraint according to diaclase occurrence index model, utilizes structure
Crack system strength model is made to establish to form diaclase discrete model;Ooze letter in the crack hole finally obtained using core analysis
Breath is established and completes diaclase equivalent model.
Finally it should be noted that:The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although
Present invention has been described in detail with reference to the aforementioned embodiments, it will be understood by those of ordinary skill in the art that:It still may be used
With technical scheme described in the above embodiments is modified or equivalent replacement of some of the technical features;
And these modifications or replacements, various embodiments of the present invention technical solution that it does not separate the essence of the corresponding technical solution spirit and
Range.