CN109143397A - Carbonate reservoir fracture hole charges recognition methods and system - Google Patents
Carbonate reservoir fracture hole charges recognition methods and system Download PDFInfo
- Publication number
- CN109143397A CN109143397A CN201710506999.0A CN201710506999A CN109143397A CN 109143397 A CN109143397 A CN 109143397A CN 201710506999 A CN201710506999 A CN 201710506999A CN 109143397 A CN109143397 A CN 109143397A
- Authority
- CN
- China
- Prior art keywords
- effective medium
- modulus
- wave impedance
- rock
- fracture hole
- 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
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Abstract
Disclose a kind of carbonate reservoir fracture hole charges recognition methods and system.This method may include: to obtain the elasticity modulus of Rock Matrix based on well-log information；Elasticity modulus based on Rock Matrix obtains the elasticity modulus of dry rock；Based on the elasticity modulus of the dry rock, EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance are obtained；Based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, carbonate reservoir fracture hole filler component is identified；The empirical relation of the fracture hole pack completeness and density established based on well-log information obtains carbonate reservoir fracture hole filling rate.The method can identify carbonate reservoir fracture hole filler component, obtain carbonate reservoir fracture hole filling rate, provide foundation to measure oil gas reservoir ability and the percolation ability of fracture hole.
Description
Technical field
The present invention relates to oil gas field of geophysical exploration, fill more particularly, to a kind of carbonate reservoir fracture hole
Object recognition methods and system.
Background technique
System in Tahe Oilfield is the Paleozoic marine facies source elephant based on ORDOVICIAN CARBONATE ancient karst oil reservoir, cave and its week
The crack of edge is its major reservoir spaces, however and not all fracture hole be effective open fracture and corrosion hole, wherein past
Toward various charges are mingled with, carbonate reservoir fracture hole charges property is studied, the oil gas for measuring fracture hole stores up
Collection ability and percolation ability have great significance.
The research achievement of forefathers focuses mostly in the research of effective open fracture and corrosion hole, to filling fracture hole, especially
The Study of recognition of charges property is less in fracture hole.Charges Quality Research key is to identify the type and filling of charges
Degree.Common fracture hole charges have shale, flour sand, organic matter, heavy mineral and the dust that collapses, corn spar etc., most of
Charges ingredient is more complex, is the mixture of many kinds of substance.And filling operation can be in identification Dong Ding, hole bottom and the base for filling thickness
It is quantitatively sought on plinth, full-filling, half filling and unfilled three classes can also be roughly divided into.
Identify that charges property can be divided into Logging Identification Method and seismic identification using geophysical method.The former
Basic means be to analyze the corresponding well logging sound of different charges using rock core and Image Logging Data calibration conventional logging information
Feature is answered, sensitive log parameter is extracted, by crossplot analysis, determines the door of sensitivity log identification different type charges
Limit value.But cross-plot, while with simple and fast advantage, different charges types is often difficult on cross plot
Completely separable, recognition effect is simultaneously inaccurate.Some scholars propose charges Logging Identification Method neural network based, seek to stitch
The mapping relations of hole charges feature and well logging information, the quality of recognition result are largely determined by the survey of selected sample
Well feature vector.The eigenvector recognition effect for meeting the more typical response characteristic of each filling type can protrude very much, and and allusion quotation
Type response differs biggish log value, and then recognition effect is bad.Therefore the precision of model identification is largely also dependent on
The accuracy of well-log information and to it is each filling type recognition capability；Opposite logging method, before seismic recognition belongs to really brill
Prediction, can provide the spatial distribution characteristic of Packing character solution cavity, find unfilled advantageous earthquake characteristic portion, instruct well location portion
Administration.The charges of seismic method identification at present mainly on the basis of carrying out the summary of solution-cavity filling object geology, are taken off by typical well
Show and establish the geological model of solution cavity and different charges actual parameters, simulation actual seismic acquisition parameter carries out forward model system
Make, analyze influence of the different Packing characters to Seismic reflection character, charges characterization description is carried out using seismic facies.But based on just
Drill to establish true geological model not a duck soup, model can only underground medium carried out a degree of simplification, fluid type,
The variation of the parameters such as saturation degree, reservoir structure may influence than charges property to analog result it is bigger, be difficult in model
It totally disappeared the difference except these factors and real medium, affect the precision using earthquake phase method identification charges.
Therefore, it is necessary to develop a kind of high-precision forecast carbonate reservoir fracture hole charges component and filling rates of capableing of
Carbonate reservoir fracture hole charges recognition methods and system.
Summary of the invention
The invention proposes a kind of carbonate reservoir fracture hole charges recognition methods and system, the method can be identified
Carbonate reservoir fracture hole filler component obtains carbonate reservoir fracture hole filling rate, for the oil gas reservoir energy for measuring fracture hole
Power and percolation ability provide foundation.
To achieve the goals above, according to an aspect of the invention, it is proposed that a kind of carbonate reservoir fracture hole charges
Recognition methods, comprising:
Based on well-log information, the elasticity modulus of Rock Matrix is obtained；
Elasticity modulus based on Rock Matrix obtains the elasticity modulus of dry rock；
Based on the elasticity modulus of the dry rock, EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance are obtained；
Based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identify that carbonate reservoir fracture hole is filled out
Fill object component；
The empirical relation of the fracture hole pack completeness and density established based on well-log information obtains the filling of carbonate reservoir fracture hole
Rate.
Preferably, described to be based on well-log information, it obtains the elasticity modulus of Rock Matrix: based on well-log information, passing through
Voight-Reuss-Hill model calculates, and obtains the elasticity modulus of the Rock Matrix, specific formula is
Wherein, MHFor the elasticity modulus of Rock Matrix, fiIndicate the volume content of mineral composition ingredient, MiIndicate mineral group
At the elasticity modulus of ingredient；Voigt and Reuss model each provides the upper limit M of equivalent Modulus of RocksVWith lower limit MR, to upper and lower
Limit carries out arithmetic mean and obtains the elasticity modulus of Rock Matrix.
Preferably, the elasticity modulus based on Rock Matrix, obtains the elasticity modulus of dry rock: blowhole is equivalent
For the ideal ellipsoid hole with single aspect ratio, Rock Matrix is added in dry equivalent hole using differential equivalent medium mode
In, obtain the elasticity modulus of the dry rock, specific formula are as follows:
Kd(0)=Km,μd(0)=μm (4)
Wherein, KdTo do rock volume modulus, μdFor dry Shear Modulus of Rock in Situ, KmFor Rock Matrix bulk modulus, μmFor rock
Ground mass matter modulus of shearing, KfFor pore-fluid bulk modulus, μfFor pore-fluid modulus of shearing, y is porosity, and P and Q are shape
The factor is equivalent aspect ratio αeFunction, for dry hole, KfAnd μfTake 0.
Preferably, the elasticity modulus based on the dry rock, obtains EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM is vertical
Wave impedance includes: to be calculated by Gassmann equation, obtains EFFECTIVE MEDIUM elasticity modulus, passes through the EFFECTIVE MEDIUM elasticity modulus
Obtain the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, specific formula are as follows:
μs=μd (6)
Wherein, KsTo be EFFECTIVE MEDIUM bulk modulus, μsFor EFFECTIVE MEDIUM modulus of shearing, KdTo do rock volume modulus, μd
For dry Shear Modulus of Rock in Situ, KfFor pore-fluid bulk modulus, KmFor matrix minerals bulk modulus, φ is porosity.
Wherein, Ip is EFFECTIVE MEDIUM p-wave impedance, and Is is EFFECTIVE MEDIUM S-wave impedance, and ρ is density, by formula (6) and
Formula (7) obtains the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance.
Specifically, it is established between charges property and geophysical parameters by formula (5), formula (6) and formula (7)
Mathematical relationship, to utilize seismic data direct inversion charges component, identification charges type provides bridge.
Preferably, described to be based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identify carbonate
Rock reservoir fracture hole filler component includes: by the EFFECTIVE MEDIUM S-wave impedance, the EFFECTIVE MEDIUM p-wave impedance and well logging
Density, the inverting under Bayes's inverting framework, inverting solution cavity section mineral constituent parameter identify that the carbonate reservoir fracture hole is filled out
Fill object component.
Preferably, the empirical relation of the fracture hole pack completeness established based on well-log information and density, obtains carbonate rock
Reservoir fracture hole filling rate includes: public according to the recurrence for having solution cavity interpretation results building solution-cavity filling rate between density of logging well
Formula obtains the spatial distribution of solution-cavity filling rate by the carbonate reservoir fracture hole filler calculation.
According to another aspect of the present invention, a kind of carbonate reservoir fracture hole charges identifying system is provided, comprising:
Based on well-log information, the unit of the elasticity modulus of Rock Matrix is obtained；
Elasticity modulus based on Rock Matrix obtains the unit of the elasticity modulus of dry rock；
Based on the elasticity modulus of the dry rock, the list of EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance is obtained
Member；
Based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identify that carbonate reservoir fracture hole is filled out
Fill the unit of object component；
The empirical relation of the fracture hole pack completeness and density established based on well-log information obtains the filling of carbonate reservoir fracture hole
The unit of rate.
Preferably, the elasticity modulus based on Rock Matrix, the unit for obtaining the elasticity modulus of dry rock includes: by rock
Stone hole is equivalent to the ideal ellipsoid hole with single aspect ratio, is added dry equivalent hole using differential equivalent medium mode
Enter in Rock Matrix, obtain the elasticity modulus of the dry rock, specific formula are as follows:
Kd(0)=Km,μd(0)=μm (4)
Wherein, KdTo do rock volume modulus, μdFor dry Shear Modulus of Rock in Situ, KmFor Rock Matrix bulk modulus, μmFor rock
Ground mass matter modulus of shearing, KfFor pore-fluid bulk modulus, μfFor pore-fluid modulus of shearing, y is porosity, and P and Q are shape
The factor is equivalent aspect ratio αeFunction, for dry hole, KfAnd μfTake 0.
Preferably, the elasticity modulus based on the dry rock, obtains EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM is vertical
The unit of wave impedance includes: to be calculated by Gassmann equation, obtains EFFECTIVE MEDIUM elasticity modulus, passes through the EFFECTIVE MEDIUM bullet
Property modulus obtain the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, specific formula are as follows:
μs=μd (6)
Wherein, KsTo be EFFECTIVE MEDIUM bulk modulus, μsFor EFFECTIVE MEDIUM modulus of shearing, KdTo do rock volume modulus, μd
For dry Shear Modulus of Rock in Situ, KfFor pore-fluid bulk modulus, KmFor matrix minerals bulk modulus, φ is porosity.
Wherein, Ip is EFFECTIVE MEDIUM p-wave impedance, and Is is EFFECTIVE MEDIUM S-wave impedance, and ρ is density, is obtained by formula (7)
Take the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance.
It is preferably based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identification carbonate rock storage
The unit of layer fracture hole filler component includes: by the EFFECTIVE MEDIUM S-wave impedance, the EFFECTIVE MEDIUM p-wave impedance and survey
Well density, the inverting under Bayes's inverting framework, inverting solution cavity section mineral constituent parameter identify the carbonate reservoir fracture hole
Filler component.
The beneficial effects of the present invention are: conventional carbonate rock rock section petrophysical model is improved, is allowed to suitable
For fracture-cavity type carbonate reservoir, the mathematical relationship between charges property and geophysical parameters is established, identifies carbonate
Rock reservoir fracture hole filler component, while the empirical relation for the fracture hole pack completeness and density established based on well-log information, obtain carbon
Carbonate Reservoir fracture hole filling rate.
Methods and apparatus of the present invention has other characteristics and advantages, these characteristics and advantages are attached from what is be incorporated herein
It will be apparent in figure and subsequent specific embodiment, or will be in the attached drawing and subsequent specific implementation being incorporated herein
It is stated in detail in example, these the drawings and specific embodiments are used together to explain specific principle of the invention.
Detailed description of the invention
Exemplary embodiment of the present is described in more detail in conjunction with the accompanying drawings, of the invention is above-mentioned and other
Purpose, feature and advantage will be apparent, wherein in exemplary embodiments of the present invention, identical reference label is usual
Represent same parts.
Fig. 1 shows carbonate reservoir fracture hole charges recognition methods process according to an embodiment of the invention
Figure.
Fig. 2 shows identification carbonate reservoir fracture hole filler component according to an embodiment of the invention modelings to imitate
Fruit schematic diagram.
Fig. 3 shows the carbonate reservoir fracture hole filling rate of one embodiment of the present of invention and the cross plot of its density.
Specific embodiment
The present invention will be described in more detail below with reference to accompanying drawings.Although showing the preferred embodiment of the present invention in attached drawing,
However, it is to be appreciated that may be realized in various forms the present invention and should not be limited by the embodiments set forth herein.On the contrary, providing
These embodiments are of the invention more thorough and complete in order to make, and can will fully convey the scope of the invention to ability
The technical staff in domain.
Embodiment 1
In this embodiment, carbonate reservoir fracture hole charges recognition methods according to the present invention, comprising: based on well logging
Data obtains the elasticity modulus of Rock Matrix；Elasticity modulus based on Rock Matrix obtains the elasticity modulus of dry rock；It is based on
The elasticity modulus of the dry rock obtains EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance；Based on the EFFECTIVE MEDIUM
S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance identify carbonate reservoir fracture hole filler component；It is built based on well-log information
The empirical relation of vertical fracture hole pack completeness and density obtains carbonate reservoir fracture hole filling rate.
The embodiment obtains the elasticity modulus of Rock Matrix elasticity modulus and dry rock by logging quality, seeks equivalent Jie
Matter S-wave impedance and EFFECTIVE MEDIUM p-wave impedance are based on EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identification
The empirical relation of carbonate reservoir fracture hole filler component, the fracture hole pack completeness and density established based on well-log information obtains carbon
Carbonate Reservoir fracture hole filling rate provides foundation to measure oil gas reservoir ability and the percolation ability of fracture hole.
The following detailed description of the specific steps of carbonate reservoir fracture hole charges recognition methods according to the present invention.
Based on well-log information, the elasticity modulus of Rock Matrix is obtained.
In one example, described to be based on well-log information, it obtains the elasticity modulus of Rock Matrix: based on well-log information, leading to
The calculating of Voight-Reuss-Hill model is crossed, obtains the elasticity modulus of the Rock Matrix, specific formula is
Wherein, MHFor the elasticity modulus of Rock Matrix, fiIndicate the volume content of mineral composition ingredient, MiIndicate mineral group
At the elasticity modulus of ingredient；Voigt and Reuss model each provides the upper limit M of equivalent Modulus of RocksVWith lower limit MR, to upper and lower
Limit carries out arithmetic mean and obtains the elasticity modulus of Rock Matrix.
Specifically, Rock Matrix mineral can be calcite, dolomite, quartz, clay etc., and elasticity modulus includes bulk modulus
K, other elasticity modulus parameters such as modulus of shearing μ, Voigt and Reuss model each provide the upper limit M of equivalent Modulus of RocksV
With lower limit MR, arithmetic mean is carried out to bound and obtains MH。
Elasticity modulus based on Rock Matrix obtains the elasticity modulus of dry rock.
In one example, based on the elasticity modulus of Rock Matrix, the elasticity modulus of dry rock is obtained: by blowhole etc.
Effect is the ideal ellipsoid hole with single aspect ratio, and rock base is added in dry equivalent hole using differential equivalent medium mode
In matter, the elasticity modulus of the dry rock, specific formula are obtained are as follows:
Kd(0)=Km,μd(0)=μm (4)
Wherein, KdTo do rock volume modulus, μdFor dry Shear Modulus of Rock in Situ, KmFor Rock Matrix bulk modulus, μmFor rock
Ground mass matter modulus of shearing, KfFor pore-fluid bulk modulus, μfFor pore-fluid modulus of shearing, y is porosity, and P and Q are shape
The factor is equivalent aspect ratio αeFunction, for dry hole, KfAnd μfTake 0.
Specifically, equivalent aspect ratio αeIt is known as pore components for the width and diameter ratio of hole, it is to define hole
The parameter of shape.When carrying out numerical simulation, the porosity type of rock complexity can be reduced to the equivalent hole of single shape,
Shape is indicated with equivalent aspect ratio.P and Q is the function of pore components, and different aspect ratios corresponds to different P, Q values, is implied
Influence of the pore shape for rock elastic property, therefore can be described as form factor.
Based on the elasticity modulus of the dry rock, EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance are obtained.
In one example, the elasticity modulus based on the dry rock obtains EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM
P-wave impedance includes: to be calculated by Gassmann equation, obtains EFFECTIVE MEDIUM elasticity modulus, passes through the EFFECTIVE MEDIUM springform
Amount obtains the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, specific formula are as follows:
μs=μd (6)
Wherein, KsTo be EFFECTIVE MEDIUM bulk modulus, μsFor EFFECTIVE MEDIUM modulus of shearing, KdTo do rock volume modulus, μd
For dry Shear Modulus of Rock in Situ, KfFor pore-fluid bulk modulus, KmFor matrix minerals bulk modulus, φ is porosity.
Wherein, Ip is EFFECTIVE MEDIUM p-wave impedance, and Is is EFFECTIVE MEDIUM S-wave impedance, and ρ is density, is obtained by formula (7)
Take the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance.
Specifically, input parameter of the porosity as petrophysical model, by the total pore space in well log interpretation performance data
Line of writing music provides.
Specifically, the petrophysical model of fracture-cavity type carbonate reservoir is established, the model to be made to be applicable not only to carbon
The rock section (being free of solution cavity) of Carbonate Reservoir, will also be suitable for solution-cavity filling object.Existing model is built both for rock section
Vertical, do not consider the influence of solution cavity and its charges, model be used for solution cavity section reluctantly, often generates biggish error,
The subsequent prediction result that solution cavity section can only be corrected by debugging model parameter, and debugging process has no theoretical foundation.One kind compared with
It for popular rock section modeling method is modified to classical Xu-White model, in view of carbonate rock and sand shale in rock
Property, the difference in porosity type, the mineral constituent in master mould is substituted for the main mine object space of carbonate rock by quartz, clay
Xie Shi and dolomite, in addition, porosity type is divided into rigid hole, intergranular in order to characterize the pore system of carbonate rock complexity
Hole and crack.However, this revised Xu-White model, the especially division of its porosity type are not appropriate for simulating
Based on mostly being filled with sand shale, the relatively single Carbonate Karst Cave section of porosity type.In order to enhance the applicability of model, introduce
Complicated porosity type is equivalent to aspect ratio be α by this concept of porosityeSingle pore morphology, construct fracture hole type
The petrophysical model of carbonate reservoir.
Based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identify that carbonate reservoir fracture hole is filled out
Fill object component.
In one example, it is based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identifies carbonic acid
Rock salt reservoir fracture hole filler component includes: by the EFFECTIVE MEDIUM S-wave impedance, the EFFECTIVE MEDIUM p-wave impedance and survey
Well density, the inverting under Bayes's inverting framework, inverting solution cavity section mineral constituent parameter identify the carbonate reservoir fracture hole
Filler component.
Specifically, pass through the EFFECTIVE MEDIUM S-wave impedance, the EFFECTIVE MEDIUM p-wave impedance and well logging density and fracture hole
Type In Carbonate Rock physical relations { Ip, Is, ρ }=f (Vcl, Vqu ...), under Bayes's inverting framework, inverting solution cavity section
The mineral constituents parameters such as shale content Vcl, sandy content Vqu realize the identification of charges lithology.
The empirical relation of the fracture hole pack completeness and density established based on well-log information obtains the filling of carbonate reservoir fracture hole
Rate.
In one example, the empirical relation of the fracture hole pack completeness and density established based on well-log information, obtains carbonate
Rock reservoir fracture hole filling rate includes: public according to the recurrence for having solution cavity interpretation results building solution-cavity filling rate between density of logging well
Formula obtains the spatial distribution of solution-cavity filling rate by the carbonate reservoir fracture hole filler calculation.
Specifically, the regression formula for constructing solution-cavity filling rate fd between density p of logging well according to existing solution cavity interpretation results
Fd=f (ρ) calculates the spatial distribution of solution-cavity filling rate using the density item in inverting solution cavity section mineral constituent parameter.
Specifically, the charges recognition methods based on petrophysical model and seismic inversion can quantitatively calculate charges component
And filling operation, circleization different lithology, the spatial distribution filled in various degree provide well location portion to the validity for determining fracture hole
Management side case has certain directive significance.
Specifically, the empirical relation of the fracture hole pack completeness and density established based on well-log information, by the density of seismic inversion
Calculate fracture hole pack completeness.Pack completeness is on the basis of well-log information counts the empirical relation of pack completeness and density, by earthquake
The density of inverting is sought.
Embodiment
Fig. 1 shows carbonate reservoir fracture hole charges recognition methods process according to an embodiment of the invention
Figure.Fig. 2 shows identification carbonate reservoir fracture hole filler component according to an embodiment of the invention modeling effects to show
It is intended to.Fig. 3 shows the carbonate reservoir fracture hole filling rate of one embodiment of the present of invention and the cross plot of its density.
By taking the ORDOVICIAN CARBONATE oil reservoir of system in Tahe Oilfield work area as an example, Development Practice shows that Caves Reservoir Body is it
Main Reservoir Body type is influenced by charges type in solution cavity scale, solution cavity and filling operation, different Filling Characteristics it is molten
Hole type Reservoir Body volume variance is larger, and identification solution-cavity filling feature has great importance before boring.
As shown in Figure 1, the rock physics mineralogical composition and content that are obtained according to interpretation of logging data, use Voight-
The elasticity modulus of Reuss-Hill model calculating Rock Matrix；Blowhole is equivalent to the ideal ellipsoid with single aspect ratio
Hole calculates the elasticity modulus of dry rock with differential equivalent medium mode；EFFECTIVE MEDIUM elasticity modulus is calculated with Gassmann equation；
To work area seismic data cube carry out prestack elastic parameter inversion, according to EFFECTIVE MEDIUM elasticity modulus establish rock physics relations,
Under Bayes's inverting framework, inverting solution cavity section mineral constituent parameter realizes the identification of charges lithology；According to existing solution cavity solution
The regression formula that fruit building solution-cavity filling rate is interpreted between density of logging well, using close in inverting solution cavity section mineral constituent parameter
It spends item and calculates solution-cavity filling rate.
The committed step of the process is accurately to construct the petrophysical model for being suitable for fracture-cavity type carbonate reservoir, root
Mathematical relationship between the charges lithology provided according to model and geophysical parameters, seeks filling by Bayes's inverting framework
Object component is filled out, identifies charges type.On the other hand, solution cavity section pack completeness data and its corresponding well logging in existing well are counted
Density can establish the regression relation between pack completeness and density, and then seek solution cavity section using density item seismic inversion and fill
Fill out rate data volume.
The precision of seismic data inverting mineral constituent is largely dependent on built petrophysical model in the suitable of the work area
Ying Xing.We select four mouthfuls in work area typical wells to test model, and test results are shown in figure 2.Black curve in figure
The calculated velocity of longitudinal wave of model and Well-Log Acoustic Velocity are respectively indicated with Grey curves, the black stripe by ordinate indicates molten
The corresponding depth bounds of hole section.Rock core information shows as full-filling or half based on showing in well that solution cavity section is mostly filled with sand shale
Filling.As can be seen that the prediction velocity of longitudinal wave of four mouthfuls of wells surveys velocity of longitudinal wave in rock section and solution cavity Duan Junneng quite well,
Predict that error is smaller, precision meets production requirement, and it is good to illustrate that the petrophysical model of this patent building has in the work area
Adaptability lays a solid foundation for the subsequent physical parameter inverting based on model.
In addition, carrying out data mining to solution cavity section interpretation results existing in work area, solution-cavity filling rate and corresponding depth are found
Well logging density degree of correlation highest at degree, as shown in Figure 3.48 pack completeness data points are had collected altogether, with corresponding density data
Intersection analysis is carried out, related coefficient illustrates that density data can be used for calculating the pack completeness of work area solution cavity, pass through prestack up to 0.95
After seismic inversion obtains density data body, the spatial distribution of pack completeness can be estimated based on the regression formula in figure.Fig. 3 is based on real
The relationship of the well log interpretation achievement of border carbonate reservoir and the pack completeness that counts and its density.Purpose is built according to real data
The empirical relation of vertical pack completeness and density, to go out density for seismic inversion and then predict that pack completeness is established using density
Basis.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
Embodiment 2
Embodiment according to the present invention provides a kind of carbonate reservoir fracture hole charges identifying system, comprising:
Based on well-log information, the unit of the elasticity modulus of Rock Matrix is obtained；
Elasticity modulus based on Rock Matrix obtains the unit of the elasticity modulus of dry rock；
Based on the elasticity modulus of the dry rock, the list of EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance is obtained
Member；
Based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identify that carbonate reservoir fracture hole is filled out
Fill the unit of object component；
The empirical relation of the fracture hole pack completeness and density established based on well-log information obtains the filling of carbonate reservoir fracture hole
The unit of rate.
In one example, the elasticity modulus based on Rock Matrix obtains the unit packet of the elasticity modulus of dry rock
It includes: blowhole is equivalent to the ideal ellipsoid hole with single aspect ratio, it will be dry etc. using differential equivalent medium mode
It imitates hole to be added in Rock Matrix, obtains the elasticity modulus of the dry rock, specific formula are as follows:
Kd(0)=Km,μd(0)=μm (4)
Wherein, KdTo do rock volume modulus, μdFor dry Shear Modulus of Rock in Situ, KmFor Rock Matrix bulk modulus, μmFor rock
Ground mass matter modulus of shearing, KfFor pore-fluid bulk modulus, μfFor pore-fluid modulus of shearing, y is porosity, and P and Q are shape
The factor is equivalent aspect ratio αeFunction, for dry hole, KfAnd μfTake 0.
In one example, the elasticity modulus based on the dry rock obtains EFFECTIVE MEDIUM S-wave impedance and equivalent
The unit of medium p-wave impedance includes: to be calculated by Gassmann equation, EFFECTIVE MEDIUM elasticity modulus is obtained, by described equivalent
Dielectric resilient modulus obtains the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, specific formula are as follows:
μs=μd(6)
Wherein, KsTo be EFFECTIVE MEDIUM bulk modulus, μsFor EFFECTIVE MEDIUM modulus of shearing, KdTo do rock volume modulus, μd
For dry Shear Modulus of Rock in Situ, KfFor pore-fluid bulk modulus, KmFor matrix minerals bulk modulus, φ is porosity.
Wherein, Ip is EFFECTIVE MEDIUM p-wave impedance, and Is is EFFECTIVE MEDIUM S-wave impedance, and ρ is density, is obtained by formula (7)
Take the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance.
In one example, it is based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identifies carbonic acid
The unit of rock salt reservoir fracture hole filler component includes: to be hindered by the EFFECTIVE MEDIUM S-wave impedance, the EFFECTIVE MEDIUM longitudinal wave
Anti- and well logging density, the inverting under Bayes's inverting framework, inverting solution cavity section mineral constituent parameter identify the carbonate rock storage
Layer fracture hole filler component.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.The selection of term used herein, purport
In the principle, practical application or improvement to the technology in market for best explaining each embodiment, or make the art
Other those of ordinary skill can understand each embodiment disclosed herein.
Claims (10)
1. a kind of carbonate reservoir fracture hole charges recognition methods, comprising:
Based on well-log information, the elasticity modulus of Rock Matrix is obtained；
Elasticity modulus based on Rock Matrix obtains the elasticity modulus of dry rock；
Based on the elasticity modulus of the dry rock, EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance are obtained；
Based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, carbonate reservoir fracture hole filler is identified
Component；
The empirical relation of the fracture hole pack completeness and density established based on well-log information obtains carbonate reservoir fracture hole filling rate.
2. carbonate reservoir fracture hole charges recognition methods according to claim 1, wherein described based on well logging money
Material, obtains the elasticity modulus of Rock Matrix: based on well-log information, being calculated by Voight-Reuss-Hill model, described in acquisition
The elasticity modulus of Rock Matrix, specific formula are
Wherein, MHFor the elasticity modulus of Rock Matrix, fiIndicate the volume content of mineral composition ingredient, MiIndicate mineral composition at
The elasticity modulus divided；Voigt and Reuss model each provides the upper limit M of equivalent Modulus of RocksVWith lower limit MR, to bound into
Row arithmetic mean obtains the elasticity modulus of Rock Matrix.
3. carbonate reservoir fracture hole charges recognition methods according to claim 1, wherein described to be based on Rock Matrix
Elasticity modulus, obtain the elasticity modulus of dry rock: blowhole being equivalent to the ideal ellipsoid hole with single aspect ratio, is made
Dry equivalent hole is added in Rock Matrix with differential equivalent medium mode, obtains the elasticity modulus of the dry rock, tool
Body formula are as follows:
Kd(0)=Km,μd(0)=μm (4)
Wherein, KdTo do rock volume modulus, μdFor dry Shear Modulus of Rock in Situ, KmFor Rock Matrix bulk modulus, μmFor rock base
Matter modulus of shearing, KfFor pore-fluid bulk modulus, μfFor pore-fluid modulus of shearing, y is porosity, and P and Q are form factor
It is equivalent aspect ratio αeFunction, for dry hole, KfAnd μfTake 0.
4. carbonate reservoir fracture hole charges recognition methods according to claim 1, wherein described to be based on the dry rock
The elasticity modulus of stone, obtains EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance includes: to be calculated by Gassmann equation,
EFFECTIVE MEDIUM elasticity modulus is obtained, the EFFECTIVE MEDIUM S-wave impedance and described etc. is obtained by the EFFECTIVE MEDIUM elasticity modulus
Imitate medium p-wave impedance, specific formula are as follows:
μs=μd (6)
Wherein, KsTo be EFFECTIVE MEDIUM bulk modulus, μsFor EFFECTIVE MEDIUM modulus of shearing, KdTo do rock volume modulus, μdIt is dry
Shear Modulus of Rock in Situ, KfFor pore-fluid bulk modulus, KmFor matrix minerals bulk modulus, φ is porosity；
Wherein, Ip is EFFECTIVE MEDIUM p-wave impedance, and Is is EFFECTIVE MEDIUM S-wave impedance, and ρ is density, obtains institute by formula (7)
State EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance.
5. carbonate reservoir fracture hole charges recognition methods according to claim 1, wherein described based on described equivalent
Medium S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identification carbonate reservoir fracture hole filler component includes: to pass through institute
State EFFECTIVE MEDIUM S-wave impedance, the EFFECTIVE MEDIUM p-wave impedance and well logging density, the inverting under Bayes's inverting framework, inverting
Solution cavity section mineral constituent parameter identifies the carbonate reservoir fracture hole filler component.
6. carbonate reservoir fracture hole charges recognition methods according to claim 1, wherein described to be based on well-log information
The fracture hole pack completeness of foundation and the empirical relation of density, obtaining carbonate reservoir fracture hole filling rate includes: according to existing solution cavity
The regression formula that interpretation results construct solution-cavity filling rate between density of logging well, passes through the carbonate reservoir fracture hole filler
Calculation obtains the spatial distribution of solution-cavity filling rate.
7. a kind of carbonate reservoir fracture hole charges identifying system, comprising:
Based on well-log information, the unit of the elasticity modulus of Rock Matrix is obtained；
Elasticity modulus based on Rock Matrix obtains the unit of the elasticity modulus of dry rock；
Based on the elasticity modulus of the dry rock, the unit of EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance is obtained；
Based on the EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, carbonate reservoir fracture hole filler is identified
The unit of component；
The empirical relation of the fracture hole pack completeness and density established based on well-log information obtains carbonate reservoir fracture hole filling rate
Unit.
8. carbonate reservoir fracture hole charges identifying system according to claim 7, wherein described to be based on Rock Matrix
Elasticity modulus, the unit for obtaining the elasticity modulus of dry rock includes: the reason for being equivalent to blowhole to have single aspect ratio
Think ellipsoid hole, dry equivalent hole is added in Rock Matrix using differential equivalent medium mode, obtains the dry rock
Elasticity modulus, specific formula are as follows:
Kd(0)=Km,μd(0)=μm (4)
Wherein, KdTo do rock volume modulus, μdFor dry Shear Modulus of Rock in Situ, KmFor Rock Matrix bulk modulus, μmFor rock base
Matter modulus of shearing, KfFor pore-fluid bulk modulus, μfFor pore-fluid modulus of shearing, y is porosity, and P and Q are form factor
It is equivalent aspect ratio αeFunction, for dry hole, KfAnd μfTake 0.
9. carbonate reservoir fracture hole charges identifying system according to claim 7, wherein described to be based on the dry rock
The elasticity modulus of stone, the unit for obtaining EFFECTIVE MEDIUM S-wave impedance and EFFECTIVE MEDIUM p-wave impedance includes: by the side Gassmann
Journey calculates, and obtains EFFECTIVE MEDIUM elasticity modulus, obtains the EFFECTIVE MEDIUM S-wave impedance by the EFFECTIVE MEDIUM elasticity modulus
And the EFFECTIVE MEDIUM p-wave impedance, specific formula are as follows:
μs=μd (6)
Wherein, KsTo be EFFECTIVE MEDIUM bulk modulus, μsFor EFFECTIVE MEDIUM modulus of shearing, KdTo do rock volume modulus, μdIt is dry
Shear Modulus of Rock in Situ, KfFor pore-fluid bulk modulus, KmFor matrix minerals bulk modulus, φ is porosity；
Wherein, Ip is EFFECTIVE MEDIUM p-wave impedance, and Is is EFFECTIVE MEDIUM S-wave impedance, and ρ is density, obtains institute by formula (7)
State EFFECTIVE MEDIUM S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance.
10. carbonate reservoir fracture hole charges identifying system according to claim 7, wherein be based on equivalent Jie
The unit of matter S-wave impedance and the EFFECTIVE MEDIUM p-wave impedance, identification carbonate reservoir fracture hole filler component includes: logical
Cross the EFFECTIVE MEDIUM S-wave impedance, the EFFECTIVE MEDIUM p-wave impedance and density of logging well, the inverting under Bayes's inverting framework,
Inverting solution cavity section mineral constituent parameter, identifies the carbonate reservoir fracture hole filler component.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710506999.0A CN109143397B (en) | 2017-06-28 | 2017-06-28 | Carbonate reservoir fracture-cave filling identification method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710506999.0A CN109143397B (en) | 2017-06-28 | 2017-06-28 | Carbonate reservoir fracture-cave filling identification method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109143397A true CN109143397A (en) | 2019-01-04 |
| CN109143397B CN109143397B (en) | 2020-05-19 |
Family
ID=64805636
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710506999.0A Active CN109143397B (en) | 2017-06-28 | 2017-06-28 | Carbonate reservoir fracture-cave filling identification method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109143397B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112649870A (en) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | Method and system for determining mineral elastic parameters in rock physical modeling |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103064118A (en) * | 2013-01-08 | 2013-04-24 | 西南石油大学 | Method of acoustic logging and quantifying cavern filling degree |
| CN104007485A (en) * | 2014-04-16 | 2014-08-27 | 孙赞东 | Method for establishing complex pore shale rock physical model |
| CN104834007A (en) * | 2015-05-04 | 2015-08-12 | 中国石油天然气股份有限公司 | Method for calculating filling degrees of carbonate fracture-cave reservoir during seismic inversion process |
| CN106772608A (en) * | 2017-02-16 | 2017-05-31 | 甘肃省地震局 | The elastic impedance and Generalized Fluid factor-analysis approach of equivalent aperture gap fracture medium |
-
2017
- 2017-06-28 CN CN201710506999.0A patent/CN109143397B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103064118A (en) * | 2013-01-08 | 2013-04-24 | 西南石油大学 | Method of acoustic logging and quantifying cavern filling degree |
| CN104007485A (en) * | 2014-04-16 | 2014-08-27 | 孙赞东 | Method for establishing complex pore shale rock physical model |
| CN104834007A (en) * | 2015-05-04 | 2015-08-12 | 中国石油天然气股份有限公司 | Method for calculating filling degrees of carbonate fracture-cave reservoir during seismic inversion process |
| CN106772608A (en) * | 2017-02-16 | 2017-05-31 | 甘肃省地震局 | The elastic impedance and Generalized Fluid factor-analysis approach of equivalent aperture gap fracture medium |
Non-Patent Citations (2)
| Title |
|---|
| 滕龙: "贝叶斯框架下联合AVA反演与统计岩石物理的储层参数估计", 《石油地球物理勘探》 * |
| 陈怀震: "基于等效各向异性和流体替换的地下裂缝地震预测方法", 《中国科学》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112649870A (en) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | Method and system for determining mineral elastic parameters in rock physical modeling |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109143397B (en) | 2020-05-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104516017B (en) | A kind of carbonate rock petrophysical parameter seismic inversion method | |
| CN103256046B (en) | Unconventionaloil pool hides method and the device that horizontal well stitches the simulation of long fracturing parameter entirely | |
| CN104950331B (en) | A kind of porosity of sand mud reservoir and the earthquake prediction method of shale content | |
| CN106368691B (en) | Three-dimensional abnormal pore pressure prediction method based on rock physics seismic information | |
| CN107121699A (en) | A kind of sedimentary facies identification method under earthquake phase control | |
| CN105445791A (en) | Stratum aperture pressure prediction method based on variety earthquake attributes | |
| CN106842301B (en) | A kind of quantitative judge and prediction technique of tuffaceous sandstone Favorable Reservoir | |
| CN103135135A (en) | Method and device used for hydrocarbon quantitative prediction and based on unconsolidated sandstone model | |
| CN107203005A (en) | A kind of method that quantification calculates crack characterising parameter | |
| CN103867198B (en) | Method for distinguishing formation density of carbonatite natural gas reservoir | |
| CN107678072B (en) | Based on magnetic force, earthquake, the united igneous reservoirs prediction technique of drilling well | |
| CN105629308A (en) | Phase control heterogeneous mechanical parameter crustal stress prediction method | |
| Sohail et al. | An evaluation of empirical and rock physics models to estimate shear wave velocity in a potential shale gas reservoir using wireline logs | |
| CN105840189A (en) | Shale gas rock physics and logging method | |
| CN105629325A (en) | Foreland basin alluvial fan fine description and prediction method | |
| Massaro et al. | Analyzing a suitable elastic geomechanical model for Vaca Muerta Formation | |
| Mandal et al. | Geomechanical appraisal and prospectivity analysis of the Goldwyer shale accounting for stress variation and formation anisotropy | |
| CN109709608A (en) | The mixed product hydrocarbon source rock TOC content of one kind and lithology compound mensuration method and application | |
| CN105804732A (en) | Oil-gas distribution prediction method based on relative resistivity stochastic simulation inversion | |
| CN109143397A (en) | Carbonate reservoir fracture hole charges recognition methods and system | |
| Qian et al. | A rock physics driven Bayesian inversion for TOC in the Fuling shale gas reservoir | |
| Paiva et al. | Comparison of methodologies to estimate the clay content–A case study in the Roncador Field, Campos Basin | |
| Jiang et al. | An efficient rock physics scheme for estimating crack density and fluid saturation of shale gas reservoir | |
| Himawan et al. | Integration of Petrophysical Analysis and Elastic Log Properties as an Input to Optimize Development Wells Target in Unique Globigerina Limestone Gas Reservoir in Madura Strait | |
| Kemper et al. | From inversion results to reservoir properties |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |