CN104749617B - A kind of multi-scale facture reservoir forward method for establishing model - Google Patents
A kind of multi-scale facture reservoir forward method for establishing model Download PDFInfo
- Publication number
- CN104749617B CN104749617B CN201310728474.3A CN201310728474A CN104749617B CN 104749617 B CN104749617 B CN 104749617B CN 201310728474 A CN201310728474 A CN 201310728474A CN 104749617 B CN104749617 B CN 104749617B
- Authority
- CN
- China
- Prior art keywords
- data
- reservoir
- seismic
- model
- grid
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to exploration geophysics data processing method field, a kind of multi-scale facture reservoir forward method for establishing model is especially related to.The horizontal non-homogeneous model towards advantage reservoir is set up first with the interval velocity data based on seismic velocity spectrum, p-wave impedance inverting data and lithological inversion data, secondly the geometry in crack is implemented using the earthquake detection method such as earthquake curvature or seismic coherence of different scale, well shake is combined and assigns its geophysical parameters, it is embedded into the horizontal non-homogeneous model being previously obtained, finally sets up the heterogeneous multi-scale facture reservoir model of transverse and longitudinal.The fracture reservoir model that the method is set up, take full advantage of well shake data, the true law of development and distribution characteristics for reflecting crack, fracture reservoir forward simulation is carried out to next step, recognize the wave field response feature and physical attribute feature of fracture reservoir, the precision for improving reservoir prediction provides basis.
Description
Technical field
The present invention relates to exploration geophysics data processing method field, a kind of multi-scale facture reservoir is especially related to
Forward model method for building up.
Background technology
Slit formation hydrocarbon reservoiring is in various rocks such as mud shale, metamorphic rock and carbonate rock.Crack in itself can be with
Oil gas is preserved as reservoir space, the passage of oil-gas migration can also be turned into, be a kind of important reservoir space type.Often at present
The foundation of crack forward model typically has two ways:First, the more commonly used is anisotropy theory and the equivalent reason in crack
By.Think that thin reservoir and fracture reservoir all show as anisotropy under normal conditions.Crack mould is set up by equivalent theory
The relational expression of shape parameter and anisotropic parameters, and then theory is just being drilled using anisotropy carry out numerical simulation.This method master
Consider the characteristic parameters such as Fracture System yardstick, density, portray not fine enough for fracture spaces development models and distribution characteristics.
2nd, crack or fracture swarm Direct Modeling.The method directly sets up fractured model, wants to set up high-precision fractured model and requires
Model meshes are small, and very big requirement is proposed for calculator memory, and computation rate is very slow.Theoretical mould is primarily adapted for use at present
Type research.
The content of the invention
The purpose of the present invention is directed to the problem of prior art presence, there is provided a kind of fracture spaces development models and distribution are special
Levy and portray multi-scale facture reservoir forward method for establishing model that is fine and being easy to operation to calculate.
The purpose of the present invention can be achieved by the following technical measures：
A kind of multi-scale facture reservoir forward method for establishing model, including：
Step 1, is built using the interval velocity data based on seismic velocity spectrum, p-wave impedance inverting data and lithological inversion data
Horizontal non-homogeneous model from facade to advantage reservoir；
Step 2, using the earthquake curvature or seismic coherence data of different scale, the heterogeneous mould of transverse direction obtained to step 1
Type enters row constraint, sets up the heterogeneous multi-scale facture reservoir forward model of transverse and longitudinal.
Refinement scheme based on such scheme includes,
The method that step 1 is set up towards the horizontal non-homogeneous model of advantage reservoir is：First with same mesh parameter
Interval velocity data based on seismic velocity spectrum, p-wave impedance inverting data and lithological inversion data are carried out into gridding, respectively
To the interval velocity grid data based on seismic velocity spectrum, p-wave impedance inverting grid data and lithological inversion grid data；Root again
Go out the calculation relational expression of wave impedance and speed according to Gardner's formula to calculating, and p-wave impedance inverting grid number is utilized by this formula
According to being calculated inverting interval velocity grid data；Then point by point scanning is carried out to lithological inversion grid data, if the mesh point is
Advantage reservoir lithology, then assign inverting interval velocity grid data to the grid；If non-advantage reservoir lithology, then the grid is assigned
Give the interval velocity grid data based on seismic velocity spectrum, after all mesh points are all scanned and assign speed data, obtain one it is new
Grid data body, the data volume is the heterogeneous forward model of transverse direction towards advantage reservoir.
The method that step 2 sets up the heterogeneous multi-scale facture reservoir forward model of transverse and longitudinal is：Using same with step 1
Mesh parameter earthquake curvature or seismic coherence data are carried out into gridding, obtain curvature or relevant grid data；Then over the ground
Shake curvature or seismic coherence grid data carry out point by point scanning, if the mesh point is fracture reservoir, the grid is assigned and is logged well
Curve interpolation speed；If not fracture reservoir, then the heterogeneous forward model speed of transverse direction towards advantage reservoir is assigned；All nets
After lattice point is all scanned and assigns speed data, a new grid data body is obtained, the data volume is multi-scale facture reservoir
Forward model.
Such scheme is further illustrated：
The interval velocity data based on seismic velocity spectrum described in step 1 refer to folded with what is extracted in seismic velocity spectrum data
Based on acceleration, by the formation interval velocity data being calculated；
Described p-wave impedance inverting data refer to utilize geological data, and the stratum obtained with reference to log data inverting is indulged
Acoustic Impedance Data；
Described lithological inversion data are a kind of using WELL LITHOLOGY data and seismic data as constraints, by being deteriorated
Function characterizes lithology data of the lithology obtained by the geostatistical inversion of three-dimensional spatial distribution form.
Earthquake curvature data described in step 2 is based on the fracture obtained by the calculating of seismic event space second dervative
With the seismic attributes data that tomography has preferable response characteristic；
The seismic coherence data be by calculating the correlation between the adjacent seismic channel of geological data, obtain reflecting by
Seismic attributes data of the caused seismic response of construction segmentation, formation lithology change in cross directional variations.
The log interpolation speed data are the interval velocities being calculated using Sonic Logging Data, with geological model
As Lateral constrained condition, the data of reflection formation interval velocity obtained from being calculated by inverse distance weighted interpolation.
Such scheme is further included：
In the calculating process of log interpolation speed, add the crack that operator virtually builds fast using pseudo- well curve
Degree information；Described pseudo- well curve be operator based on the logging character and geological knowledge of apriority, by non-drilling well
Virtual the constructing in position can reflect the position geologic feature（It is preferred that fracture reservoir feature）Pseudo- well curve.It is preferred that crack is stored up
Layer feature is used as the position geologic feature.
The multi-scale facture reservoir forward model set up by the method, takes full advantage of well shake data, and FORWARD AND INVERSE PROBLEMS is melted
Conjunction mode, truly reflects the law of development and distribution characteristics in crack.The foundation of the forward model, crack storage is carried out to next step
Layer forward simulation, recognizes the wave field response feature and physical attribute feature of fracture reservoir, improves the precision of reservoir prediction
Provide the foundation.
Brief description of the drawings
Fig. 1 is a kind of flow chart of specific embodiment of the invention；
Fig. 2 is lithological inversion data generalized section；
Fig. 3 is the interval velocity data generalized section based on seismic velocity spectrum；
Fig. 4 is that p-wave impedance inverting data section is calculated inverting interval velocity data section；
Fig. 5 is the horizontal non-homogeneous model schematic diagram towards advantage reservoir；
Fig. 6 is coherence data schematic diagram；
Fig. 7 is log interpolation speed data generalized section；
Fig. 8 is multi-scale facture reservoir forward model schematic.
Specific embodiment
Embodiment 1, is described further with reference to accompanying drawing 1 to technical scheme first.
Step 1, is built using the interval velocity data based on seismic velocity spectrum, p-wave impedance inverting data and lithological inversion data
Horizontal non-homogeneous model from facade to advantage reservoir.Each data implication and specific method for building up are as described below.
Each data implication：Interval velocity data based on seismic velocity spectrum refer to the superposition to be extracted in seismic velocity spectrum data
Based on speed, by the formation interval velocity data being calculated.
P-wave impedance inverting data refer to utilize geological data, and the stratum p-wave impedance obtained with reference to log data inverting
Data.
Lithological inversion data are a kind of based on geological model, using WELL LITHOLOGY data and seismic data as constraints,
Lithology data of the lithology obtained by the geostatistical inversion of three-dimensional spatial distribution form is characterized by variogram.
Specific method for building up：In step 1, the interval velocity of seismic velocity spectrum will be based on first with same mesh parameter
Data, p-wave impedance inverting data and lithological inversion data carry out gridding, respectively obtain the interval velocity based on seismic velocity spectrum
Grid data, p-wave impedance inverting grid data and lithological inversion grid data.Further according to Gardner's formula（Density=0.31 ×
Speed1/4）Extrapolate calculation relational expression (wave impedance=0.31 × speed of wave impedance and speed5/4), and by this formula using vertical
Wave impedance inversion grid data is calculated inverting interval velocity grid data.Then pointwise is carried out to lithological inversion grid data to sweep
Retouch, if the mesh point is advantage reservoir lithology（According to different geological conditions, by operator according to lithological inversion grid data from
Row formulation standard specify what be advantage reservoir lithology what be non-advantage reservoir lithology）, then inverting interval velocity net is assigned to the grid
Lattice data；If non-advantage reservoir lithology, then the interval velocity grid data based on seismic velocity spectrum is assigned to the grid.All nets
After lattice point is all scanned and assigns speed data, a new grid data body is obtained, the data volume is towards advantage reservoir
Horizontal heterogeneous forward model.
Step 2, the transverse direction towards advantage reservoir obtained to step 1 using earthquake curvature or seismic coherence data is heterogeneous
Forward model carries out crack description, and fracture assigns log interpolation speed data to set up multi-scale facture reservoir forward
Model.Each data implication and specific method for building up are as described below.
Each data implication：Earthquake curvature data be based on seismic event space second dervative calculate obtained by fracture and
Tomography has the seismic attributes data of preferable response characteristic.
Seismic coherence data are by calculating the correlation between the adjacent seismic channel of geological data, obtaining reflecting by constructing
Seismic attributes data of the caused seismic response of segmentation, formation lithology change in cross directional variations.
Log interpolation speed data are the interval velocities being calculated using Sonic Logging Data, using geological model as
Lateral constrained condition, the data of reflection formation interval velocity obtained from being calculated by inverse distance weighted interpolation；Due to log
Longitudinal high-resolution, log interpolation speed can include crack velocity information.In the calculating of log interpolation speed
During, pseudo- well curve can also be used【That is operator based on the logging character and geological knowledge of apriority, by not
Virtual the constructing in the position of drilling well can reflect the position geologic feature（Consider fracture reservoir feature emphatically in the present invention）'s
Well curve】, the crack velocity information that operator fabricates, the multi-scale facture for so finally giving can be added using pseudo- well curve
Reservoir forward model just adds operator's viewpoint.
Specific method for building up：In step 2, using the mesh parameter same with step 1 by earthquake curvature or seismic coherence
Data（The final result that two kinds of data are obtained is equal to, and is voluntarily selected according to actual conditions by operator）Gridding is carried out, is obtained
Curvature or relevant grid data.Then point by point scanning is carried out to earthquake curvature or seismic coherence grid data, if the mesh point is
Fracture reservoir（According to different geological conditions, standard gauge is voluntarily formulated according to earthquake curvature or seismic coherence data by operator
Determine fracture reservoir, multiple standards can be set up and represent multi-scale facture reservoir）, then log interpolation speed is assigned to the grid；
If not fracture reservoir, then the heterogeneous forward model speed of transverse direction towards advantage reservoir is assigned.All mesh points are all scanned simultaneously
After assigning speed data, a new grid data body is obtained, the data volume is multi-scale facture reservoir forward model.
Embodiment 2, with reference to shown in accompanying drawing 1-8, is described in detail below.
Fig. 1 is the flow chart of the specific embodiment of multi-scale facture reservoir forward method for establishing model of the invention.
Step 101, the interval velocity data of seismic velocity spectrum, p-wave impedance inverting number will be based on using same mesh parameter
Gridding is carried out according to, lithological inversion data, the interval velocity grid data based on seismic velocity spectrum, p-wave impedance inverting is respectively obtained
Grid data and lithological inversion grid data.
Step 102, inverting interval velocity data are calculated using p-wave impedance inverting data.
Step 103, point by point scanning is carried out to lithological inversion grid data, judges whether each model meshes is advantage reservoir
Lithology.
If the grid is advantage reservoir lithology, inverting interval velocity grid data is assigned to it；If the grid is not advantage
Reservoir lithology, then assign the interval velocity grid data based on seismic velocity spectrum to it.
All grids all judge to finish, then flow enters 104, that is, obtain that the transverse direction towards advantage reservoir is heterogeneous just to drill mould
Type.
Step 101 can be applicable to 104, set up the horizontal non-homogeneous model towards advantage reservoir.
In step 105, by earthquake curvature or seismic coherence data gridding.
In step 106, point by point scanning is carried out to earthquake curvature or seismic coherence grid data.If fracture reservoir, then
Log interpolation speed is assigned to the grid；If not fracture reservoir, then assign heterogeneous just towards the transverse direction of advantage reservoir
Drill model velocity.
All grids all judge to finish, then flow enters 107, that is, obtain multi-scale facture reservoir forward model.
Fig. 2 is lithological inversion data generalized section；Fig. 3 is that the interval velocity data section based on seismic velocity spectrum is illustrated
Figure；Fig. 4 is the inverting interval velocity data generalized section that step 102 is obtained；Fig. 5 be step 104 obtain towards advantage reservoir
Horizontal non-homogeneous model schematic diagram.
Fig. 6 is coherence data schematic diagram；Fig. 7 is log interpolation speed data generalized section；Fig. 8 is step 106
The multi-scale facture reservoir forward model schematic for obtaining.
The fracture reservoir model that the method is obtained, takes full advantage of well shake data, really reflects the development rule in crack
Rule and distribution characteristics.
Claims (6)
1. a kind of multi-scale facture reservoir forward method for establishing model, it is characterised in that including：
Step 1, by the use of lithological inversion data as constraints, the layer speed that p-wave impedance inverting is obtained is assigned to advantage reservoir
Degree, the interval velocity that seismic velocity spectrum is obtained is assigned to non-advantage reservoir, so sets up heterogeneous just towards the transverse direction of advantage reservoir
Drill model；
Step 2, by the use of the earthquake curvature or seismic coherence data of different scale as constraints, fracture reservoir assigns statistics
The crack speed or acoustic logging interpolation speed for obtaining, the transverse direction that step 1 foundation is assigned to non-fracture reservoir are heterogeneous just
Model velocity is drilled, the heterogeneous multi-scale facture reservoir forward model of transverse and longitudinal is so set up.
2. multi-scale facture reservoir forward method for establishing model according to claim 1, it is characterised in that
The method that step 1 is set up towards the horizontal non-homogeneous model of advantage reservoir is：First with same mesh parameter by base
Interval velocity data, p-wave impedance inverting data and lithological inversion data in seismic velocity spectrum carry out gridding, respectively obtain base
Interval velocity grid data, p-wave impedance inverting grid data and lithological inversion grid data in seismic velocity spectrum；Further according to plus
Moral receives formula to calculating and goes out the calculation relational expression of wave impedance and speed, and utilizes p-wave impedance inverting grid data meter by this formula
Calculation obtains inverting interval velocity grid data；Then point by point scanning is carried out to lithological inversion grid data, if the mesh point is advantage
Reservoir lithology, then assign inverting interval velocity grid data to the grid；If non-advantage reservoir lithology, then base is assigned to the grid
In the interval velocity grid data of seismic velocity spectrum, after all mesh points are all scanned and assign speed data, a new net is obtained
Lattice data volume, the data volume is the heterogeneous forward model of transverse direction towards advantage reservoir；
The method that step 2 sets up the heterogeneous multi-scale facture reservoir forward model of transverse and longitudinal is：Using the net same with step 1
Earthquake curvature or seismic coherence data are carried out gridding by lattice parameter, obtain curvature or relevant grid data；Then it is bent to earthquake
Rate or seismic coherence grid data carry out point by point scanning, if the mesh point is fracture reservoir, log are assigned to the grid
Interpolation speed；If not fracture reservoir, then the heterogeneous forward model speed of transverse direction towards advantage reservoir is assigned；All mesh points
After all scanning and assigning speed data, a new grid data body is obtained, the data volume is multi-scale facture reservoir forward
Model.
3. multi-scale facture reservoir forward method for establishing model according to claim 1 and 2, it is characterised in that
The interval velocity data based on seismic velocity spectrum described in step 1 refer to the superposition speed to be extracted in seismic velocity spectrum data
Based on degree, by the formation interval velocity data being calculated；
Described p-wave impedance inverting data refer to utilize geological data, and the stratum compressional wave resistance obtained with reference to log data inverting
Anti- data；
Described lithological inversion data are a kind of using WELL LITHOLOGY data and seismic data as constraints, by variogram
Characterize lithology data of the lithology obtained by the geostatistical inversion of three-dimensional spatial distribution form.
4. multi-scale facture reservoir forward method for establishing model according to claim 1 and 2, it is characterised in that
Earthquake curvature data described in step 2 is based on the fracture obtained by the calculating of seismic event space second dervative and breaks
Layer has the seismic attributes data of preferable response characteristic；
The seismic coherence data are by calculating the correlation between the adjacent seismic channel of geological data, obtaining reflecting by constructing
Seismic attributes data of the caused seismic response of segmentation, formation lithology change in cross directional variations；
The log interpolation speed data are the interval velocities being calculated using Sonic Logging Data, using geological model as
Lateral constrained condition, the data of reflection formation interval velocity obtained from being calculated by inverse distance weighted interpolation.
5. multi-scale facture reservoir forward method for establishing model according to claim 4, it is characterised in that in log
In the calculating process of interpolation speed, the virtual crack velocity information for building of operator is added using pseudo- well curve；Described pseudo- well
Curve be operator based on the logging character and geological knowledge of apriority, virtually construct energy by the position of non-drilling well
Enough reflect the pseudo- well curve of the position geologic feature.
6. multi-scale facture reservoir forward method for establishing model according to claim 5, it is characterised in that in non-drilling well
Position virtually builds during pseudo- well curve preferred fracture reservoir feature as the position geologic feature.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310728474.3A CN104749617B (en) | 2013-12-26 | 2013-12-26 | A kind of multi-scale facture reservoir forward method for establishing model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310728474.3A CN104749617B (en) | 2013-12-26 | 2013-12-26 | A kind of multi-scale facture reservoir forward method for establishing model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104749617A CN104749617A (en) | 2015-07-01 |
| CN104749617B true CN104749617B (en) | 2017-05-31 |
Family
ID=53589587
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310728474.3A Active CN104749617B (en) | 2013-12-26 | 2013-12-26 | A kind of multi-scale facture reservoir forward method for establishing model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104749617B (en) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105652323B (en) * | 2015-08-11 | 2017-12-29 | 中国石油化工股份有限公司 | A kind of method for predicting reservoir |
| CN106597531B (en) * | 2015-10-16 | 2019-10-29 | 中国石油化工股份有限公司 | The Forward Modeling of the wave field propagation characteristic of shale containing vertical fracture |
| CN106842290B (en) * | 2015-12-04 | 2019-12-06 | 中国石油化工股份有限公司 | plane representation method of torsion structure |
| CN106842326B (en) * | 2015-12-04 | 2020-10-13 | 中国石油化工股份有限公司 | Sand-mud interbed CO2 geological sequestration time-lapse seismic forward modeling method during non-transverse-wave velocity logging |
| CN107843923A (en) * | 2016-09-19 | 2018-03-27 | 中国石油化工股份有限公司 | Curvature analysis method based on earthquake attribute volume |
| CN108020648A (en) * | 2016-10-28 | 2018-05-11 | 核工业北京地质研究院 | A kind of method of rapid preliminary identification rock fracture different spaces distribution characteristics |
| CN106556868B (en) * | 2016-11-01 | 2019-05-07 | 中国石油天然气股份有限公司 | Quantitative identification method and device for groove |
| CN108572392B (en) * | 2017-03-09 | 2020-07-28 | 中国石油化工股份有限公司 | Method for depicting carbonate rock fracture-cave communicating body |
| JP6937391B2 (en) * | 2018-01-31 | 2021-09-22 | 富士フイルム株式会社 | Repair length determination method and repair length determination device |
| CN109696704B (en) * | 2019-01-28 | 2020-05-01 | 中国海洋石油集团有限公司 | Seismic anisotropy delta modeling method based on longitudinal wave impedance constraint |
| CN109799531B (en) * | 2019-03-20 | 2022-04-01 | 成都理工大学 | Fracture reservoir prediction method based on seismic frequency division coherence attribute |
| CN110988988B (en) * | 2019-11-25 | 2020-10-23 | 中国矿业大学(北京) | Seismic wave field simulation method and device based on vertical fracture medium |
| CN113138407B (en) * | 2020-01-20 | 2024-05-03 | 中国石油天然气集团有限公司 | Multi-scale fracture earthquake prediction method and system for deep shale gas |
| CN112068198B (en) * | 2020-08-24 | 2022-03-18 | 西南科技大学 | Crack fracture dimension description method based on seismic wave full waveform characteristics |
| CN112946782B (en) * | 2021-03-15 | 2022-03-25 | 西南石油大学 | Earthquake fine depicting method for dense oil-gas storage seepage body |
| CN113777656A (en) * | 2021-08-25 | 2021-12-10 | 中国石油化工股份有限公司 | Shale reservoir inversion model building method based on structural deformation |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0745870A2 (en) * | 1995-06-01 | 1996-12-04 | Phillips Petroleum Company | Non-unique seismic lithologic inversion for subterranean modeling |
| CN102121995A (en) * | 2010-12-03 | 2011-07-13 | 中国石油天然气集团公司 | Seismic inversion reservoir prediction method for complex structure containing inversion fracture |
| CN102253415A (en) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | Method for establishing earthquake response mode based on fracture equivalent medium model |
-
2013
- 2013-12-26 CN CN201310728474.3A patent/CN104749617B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0745870A2 (en) * | 1995-06-01 | 1996-12-04 | Phillips Petroleum Company | Non-unique seismic lithologic inversion for subterranean modeling |
| CN102121995A (en) * | 2010-12-03 | 2011-07-13 | 中国石油天然气集团公司 | Seismic inversion reservoir prediction method for complex structure containing inversion fracture |
| CN102253415A (en) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | Method for establishing earthquake response mode based on fracture equivalent medium model |
Non-Patent Citations (3)
| Title |
|---|
| 三维速度场建立方法;韩宏伟，等;《油气地质与采收率》;20100131;第17卷(第1期);54-56、61 * |
| 正演模拟技术在先巴扎地区储层预测中的应用;宁松华，等;《石油天然气学报》;20130331;第35卷(第3期);58-60、71 * |
| 玉北地区碳酸盐岩储层地震响应特征研究;张旭光;《石油物探》;20120930;第51卷(第5期);493-501 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104749617A (en) | 2015-07-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104749617B (en) | A kind of multi-scale facture reservoir forward method for establishing model | |
| CN101421640B (en) | Methods and apparatus for predicting the hydrocarbon production of a well location | |
| CN102759745B (en) | Method for forecasting carbonate reservoir based on forward modeling of digital geological outcrop model | |
| CN106368691B (en) | Three-dimensional abnormal pore pressure prediction method based on rock physics seismic information | |
| CN102636812B (en) | Method for obtaining storage space volume of carbonate reservoir | |
| CN104898161B (en) | Effective sandstone predicting method based on logging response simulator | |
| CN107121699A (en) | A kind of sedimentary facies identification method under earthquake phase control | |
| CN105527653B (en) | A kind of virtual log construction method based on geological information | |
| CN106855636A (en) | Based on the prototype geological model Seismic forward method that carbonate reservoir is appeared | |
| CN102937720B (en) | Well control improves the method for seismic data resolution | |
| CN105986815A (en) | Method for recognizing geological sweet spots in shale stratum | |
| CN101634717A (en) | Fine shear-wave (S-wave) impedance access technology based on logging and prestack channel set seismic data | |
| CN103439740B (en) | Method and device for predicting relative impedance based on dipole seismic wavelet multiple integral | |
| CN105089615A (en) | Log data historical retrogression treatment method based on oil reservoir model | |
| CN105629308A (en) | Phase control heterogeneous mechanical parameter crustal stress prediction method | |
| CN104297800B (en) | A kind of from phased prestack inversion method | |
| CN107544093A (en) | The structure interpretation layer depth system compensation method of borehole restraint | |
| CN102288996B (en) | Method for quantitatively predicting gas containing probability by AFI (prestack inversion) | |
| CN106096081B (en) | The estimation method of reserve of fracture hole type bottom water reservoir | |
| Patterson et al. | Integration of microseismic data, completion data, and production data to characterize fracture geometry in the Permian Basin | |
| CN104834007B (en) | Method for calculating filling degree of carbonate fracture-cave reservoir in seismic inversion process | |
| CN102565852B (en) | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir | |
| CN109696704B (en) | Seismic anisotropy delta modeling method based on longitudinal wave impedance constraint | |
| CN110095811A (en) | Method and device for constructing and processing velocity model of paste rock stratum | |
| CN106761647A (en) | A kind of method that plane transforms area after estimation shale reservoir pressure |
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 |