CN104975851B - For amplitude with the reservoir model optimization method of geophone offset variation road set analysis - Google Patents
For amplitude with the reservoir model optimization method of geophone offset variation road set analysis Download PDFInfo
- Publication number
- CN104975851B CN104975851B CN201410142188.3A CN201410142188A CN104975851B CN 104975851 B CN104975851 B CN 104975851B CN 201410142188 A CN201410142188 A CN 201410142188A CN 104975851 B CN104975851 B CN 104975851B
- Authority
- CN
- China
- Prior art keywords
- model
- parameter
- reservoir model
- avo
- velocity
- 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 is the reservoir model optimization method for amplitude with geophone offset variation road set analysis, establish reservoir model, petrophysical model is demarcated, establish the non-linear relation between seismic response parameter and rock parameter, using calibrated petrophysical model, the velocity of longitudinal wave, shear wave velocity and density of reservoir model are calculated；Velocity of longitudinal wave is multiplied the impedance synthesized with density again, the angular range and angle interval that up-and-down boundary, wavelet and the earthquake of forward modeling of reservoir model are arrived is determined with well-log information, generates AVO prestack trace gather forward modeling result.The present invention provides foundation for the modification of geological model, and the available and closest reservoir model of practical subsurface picture is conducive to oil reservoir development.
Description
Technical field
The present invention relates to the field of seismic exploration of petroleum, more particularly to a kind of amplitude that is used for is with geophone offset variation road set analysis
Reservoir model optimization method.
Background technique
Amplitude with geophone offset variation (AVO) analysis be after bright point technique another item using amplitude information research lithology,
Detect the important method of oil gas.Earthquake reflected amplitude is analyzed with the changing rule of geophone offset (incidence angle) in prestack, into
And determine and covered on reflecting interface, the lithologic character transitivity parameter for the medium that underlies, whereby to the fluid properties in hole in rock
Deduction is made with lithology.
In reflection wave seismic exploration most important hypothesis first is that: the interface between the different medium of two kinds of wave impedance is
Reflecting surface, after the wave impedance when interface both sides determines, the reflection coefficient of reflecting interface has also just uniquely been determined.When according to reflection
Between and echo amplitude, can be with certain main characters of the structural configuration of inverting subsurface interface and medium.However actual seismic is anti-
Penetrating coefficient is the variation with geophone offset and changes, and AVO technology exactly carries out hydrocarbon fluids analysis using this variation.
Zoeppritz equation is the theoretical basis of AVO technical research, but since the equation is more complicated, is never directly answered
With descendant summarizes the reduced form of many Zoeppritz on the basis of Zoeppritz.
Traditional AVO technology, in conjunction with characteristics of reservoirs, is analyzed under different geological conditions using forward modeling AVO phenomenon
The AVO feature that the different fluids such as oil, gas and water generate establishes corresponding AVO mark, in actual earthquake record, Direct Recognition
Lithology and oil-gas possibility.
However, the actual conditions error of ordinary circumstance drag and underground can be bigger, the reason is that reservoir model usually by
Well-log information, seismic data, creation data and core data, are constructed by history matching, only understand ground by well-log information
The actual conditions of lower certain point, although and seismic data is solved to a certain extent to lateral understanding, it is longitudinal
Limited resolution, and in fact, subsurface geology situation be it is complicated, modeling result is often larger with actual conditions error,
Therefore, the AVO result of forward modeling may bring deviation to explanation underground oil gas.
Summary of the invention
Object of the present invention is to propose that a kind of forward modeling result can change for amplitude with geophone offset closer in actual
The reservoir model optimization method of road set analysis.
The present invention is realized by following steps:
1) quasi- by history according to work area seismic interpretation data, log analysis data, creation data and core data
It closes, establishes reservoir model；
Data described in step 1) include: fluid parameter,
Elastic parameter, oil reservoir partitioned parameters, initialization calculating parameter, output control parameter, simulates work at rock properties
Essential information and manufacturing parameter；Wherein:
Fluid parameter includes the ground surface density or severe of oil, gas and water, the reservoir volume ratio of oil, gas, the viscosity of water, body
Product parameter, the compressed coefficient, rock compressibility；
Elastic parameter includes shear wave velocity, Lame constants and Poisson's ratio.
Reservoir model is established described in step 1), is to establish Skeleton Model using seismic interpretation result, in conjunction with log data
And core data, physical property model is established, history matching is carried out on the basis of creation data and obtains reservoir model.
2) petrophysical model is demarcated using work area log analysis data and reservoir model data, establishes seismic response ginseng
Non-linear relation between several and rock parameter；
The petrophysical model is Voigt and Reuss theoretical model or Gassmann theoretical model.
The non-linear relation established between seismic response parameter and rock parameter is to establish acoustic wave parameter, elasticity ginseng
The several and relationship between physical properties of rock, fluid parameter.
3) oil reservoir is calculated using calibrated petrophysical model according to the fluid parameter of reservoir model and elastic parameter
Velocity of longitudinal wave, shear wave velocity and the density of model；Velocity of longitudinal wave is multiplied the impedance synthesized with density again；
The petrophysical model is Voigt and Reuss theoretical model or Gassmann theoretical model.
Described calculate uses gassmann model, and rock volume elasticity modulus K is calculated:
Wherein, KsFor rock particles elasticity modulus；KdFor dry elastic modulus of rock；Φ is porosity；
KfFor elasticity of fluid modulus, calculated by following formula:
Wherein, Kw, Ko, KgIt is the elasticity modulus of water, oil, gas, S respectivelyw,
SoIt is the saturation degree of water, oil, 1-S respectivelyw-SoFor gas saturation；
Velocity of longitudinal wave V of the sound wave in rockpIt is calculated by following formula:
Wherein, μ is Shear Modulus of Rock in Situ；ρ is rock density；Vs is shear wave
Speed.
Wherein, ρw、ρo、ρgBe respectively water, oil,
Air tightness, ρmaFor matrix density.
Resultant impedance are as follows: AI=ρ Vp。
4) angle model that up-and-down boundary, wavelet and the earthquake of forward modeling of reservoir model are arrived is determined using work area well-log information
It encloses and carries out the trace gather forward modeling of AVO prestack using velocity of longitudinal wave, shear wave velocity and density with angle interval, generate AVO prestack trace gather
Forward modeling result；
Before trace gather forward modeling, the acoustic wave parameter and density value surveyed at invalid grid using well-log information are to reservoir model
In velocity of longitudinal wave at invalid grid and density set；
The well-log information is density log curve and acoustic travel time logging curve.
It is described be set as by well-log information density and speed be assigned to it is invalid in the reservoir model of corresponding depth
Influence at grid, to avoid null value of the reservoir model in forward modeling in invalid grid to forward modeling result.
The phase of the wavelet and frequency and the sub-wave phase of actual seismic data are consistent with frequency.
Before trace gather forward modeling, according between the velocity of longitudinal wave, shear wave velocity, density and the angular range of setting and angle
Every calculating reflection coefficient；
The method of the calculating reflection coefficient includes boltfeld or Aki.
5) seismic data of AVO prestack trace gather forward modeling result and actual acquisition is compared, completes oil reservoir if consistent
The optimization of model, such as inconsistent modification are consistent.
The comparison is:
The AVO trace gather for the well point position sometime put, compares with the production logging curve at time point or AVO is folded
The AVO of preceding trace gather forward modeling result and actual seismic acquisition time consistent AVO trace gather and actual acquisition seismic data analyze into
Row comparison.
Described is unanimously:
AVO characteristic meets the rule of the fluid at well point, and the consistent indistinction of AVO characteristic, amplitude increases with the increase of geophone offset
It is big or reduce, in addition occur phase reversion or,
AVO prestack trace gather forward modeling result and actual seismic acquisition time consistent AVO trace gather and actual acquisition earthquake number
According to AVO analysis compare, whether AVO characteristic consistent.
The modification is to the porosity or shale content or net-gross ratio or oil saturation or pressure in reservoir model
Power or temperature are modified.
Optimization method for building up of the present invention for the reservoir model of the road AVO set analysis, provides for the modification of geological model
Foundation, and the modification of conventional geological model is usually the experience according to reservoir engineer, thus there is many uncertainties,
Therefore the available and closest reservoir model of practical subsurface picture of the invention, is conducive to oil reservoir development.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those skilled in the art without any creative labor, can be with root
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is that the reservoir model for the road AVO set analysis of the embodiment of the present invention optimizes the method flow diagram of method for building up；
Fig. 2 is that the reservoir model for the road AVO set analysis of the embodiment of the present invention optimizes the parameter calculating stream of method for building up
Cheng Tu；
Fig. 3 is the model flat distribution map of specific example trial zone of the present invention；
Fig. 4 is that the petrophysical model of specific example trial zone of the present invention demarcates schematic diagram；
Fig. 5,6,7 are respectively the velocity of longitudinal wave, shear wave velocity, close of specific example test block of the present invention reservoir model calculating
Degree；
Fig. 8 is velocity of longitudinal wave and the calculated impedance of density of specific example test block of the present invention；
Fig. 9 and Figure 10 is respectively setting up-and-down boundary and the wavelet signal of specific example trial zone of the present invention reservoir model
Figure；
Figure 11 is the angular range of the setting of the selected well point the L8-182 position of specific example trial zone of the present invention, angle
Interval and the schematic diagram for calculating reflection coefficient method；
Figure 12 is the PP prestack trace gather schematic diagram of the well point the L8-182 position synthesis of specific example trial zone of the present invention；
Figure 13 is the schematic diagram of the AVO prestack road set analysis of the well point the L8-182 position of specific example trial zone of the present invention；
Figure 14 is the actually measured water saturation curve of the well point the L8-182 position of specific example trial zone of the present invention.
Specific embodiment
The embodiment of the present invention provides a kind of optimization method for building up of reservoir model for the road AVO set analysis, comprising: according to
Seismic interpretation data, log analysis data, creation data and core data establish reservoir model by history matching；It utilizes
Work area log analysis data and reservoir model data demarcate petrophysical model, establish seismic response parameter and rock parameter it
Between non-linear relation；It is calculated according to the fluid parameter and elastic parameter of reservoir model using calibrated petrophysical model
Velocity of longitudinal wave, shear wave velocity and the density of reservoir model；Velocity of longitudinal wave is multiplied the impedance synthesized with density again；It utilizes
Work area well-log information determines the angular range and angle interval that up-and-down boundary, wavelet and the earthquake of forward modeling of reservoir model are arrived,
Using velocity of longitudinal wave, shear wave velocity and density, the trace gather forward modeling of AVO prestack is carried out, generates AVO prestack trace gather forward modeling result；By AVO
The seismic data of prestack trace gather forward modeling result and actual acquisition compares, and the optimization of reservoir model is completed if consistent, if not
Consistent modification is consistent.
Following will be combined with the drawings in the embodiments of the present invention, carries out to the technical solution in present example clear, complete
Ground description, it is clear that described embodiment is only a part of example of the present invention, instead of all the embodiments.Based on this hair
Example in bright, every other example obtained by those of ordinary skill in the art without making creative efforts,
It shall fall within the protection scope of the present invention.
Fig. 1 is that the reservoir model for the road AVO set analysis of the embodiment of the present invention optimizes the method flow diagram of method for building up,
Stream is calculated with reference to the parameter that the reservoir model for the road AVO set analysis that Fig. 2 is the embodiment of the present invention optimizes method for building up together
Cheng Tu.As depicted in figs. 1 and 2, the reservoir model for the road AVO set analysis of the invention optimizes method for building up, and core is rock
Pre-stack Forward Modeling based on reservoir model after physical model calibration and calibration, by velocity of longitudinal wave, the density etc. that generate after calibration,
Wave impedance is formed, the up-and-down boundary information and wavelet of reservoir model is set, is finally completed the AVO prestack trace gather of different angle, into
And it is compared and analyzed with well real fluid situation.In Fig. 1 of the present embodiment, the reservoir model for the road AVO set analysis is excellent
Changing method for building up includes:
Step S101 passes through history according to seismic interpretation data, log analysis data, creation data and core data
Fitting, establishes reservoir model.Geologist carries out detailed pool description work to oil field, is based on seismic interpretation data, well logging
Explain data, core data and combination geology man establish three-dimensional geological model to the understanding in oil field.According to digital-to-analogue engineer's
It is required that geology man carries out roughening treatment to three-dimensional geological model, directly digital-to-analogue engineer can be met for digital-to-analogue engineer output
The reservoir simulation model needed.Wherein, reservoir model data include: fluid parameter, elastic parameter, rock properties, oil reservoir subregion
Parameter, initialization calculating parameter, output control parameter, the essential information and manufacturing parameter for simulating work；Wherein, fluid parameter
Ground surface density or severe including oil, gas and water, the reservoir volume ratio of oil, gas, the viscosity of water, volume parameter, the compressed coefficient,
Rock compressibility；Elastic parameter includes shear wave velocity, Lame constants and Poisson's ratio etc..
Step S102 demarcates petrophysical model using work area log analysis data and reservoir model data, establishes ground
Shake the non-linear relation between response parameter and rock parameter.There is no speed and density these values in reservoir model data, therefore
It needs to be fitted using some dynamic parameters of logging data interpretation result well logging position corresponding with reservoir model, such one
A process is referred to as to demarcate.It is then sharp first the speed applications in well logging at the corresponding well point of reservoir model by taking speed as an example
With the corresponding relationship of the dynamic parameter at the well point Chu Yufei of well point, synthesize the speed at non-well point.That is: petrophysical model mark
Surely be will by carrying out comprehensive analysis to various core data and logging informations and seismic data, study of rocks lithology, porosity,
Porosity type, pore-fluid, fluid saturation and frequency etc. are to such as longitudinal wave of resilient nature in rock (P wave) and shear wave (S wave) speed
The influence of degree and decaying, and propose the theory and method using seismic response prediction rock physical property, it is seismic response and storage
The bridge contacted between layer rock parameter is the basic premise for carrying out quantitative reservoir prediction.Characterize the ground of petrophysics feature
Shake parameter mainly has elasticity modulus, density, velocity of longitudinal wave, shear wave velocity of rock etc., they are that we identify lithology and oil gas
Important parameter, and carry out the bridge of quantitative seismic reservoir characterization.
In the present embodiment, the petrophysical model of foundation includes: Voigt and Reuss theoretical model, Gassmann theoretical
Model.
Wherein, the Gassmann theoretical model are as follows:
Rock volume elasticity modulus calculates:
Wherein, Ks is rock particles elasticity modulus；Kd is dry elastic modulus of rock；
Φ is porosity；
Kf is elasticity of fluid modulus, can be calculated by following formula:
Wherein, Kw, Ko, Kg are the elasticity modulus of water, oil, gas respectively, Sw,
So, Sg are the saturation degree of water, oil, gas respectively；
Velocity of longitudinal wave of the sound wave in rock is calculated by following formula:
Wherein, μ is Shear Modulus of Rock in Situ；ρ is rock density, it and the linear pass of saturation degree
System:
Wherein, ρw、ρo、ρgBe respectively water, oil,
Air tightness, ρw、ρo、ρgIt needs to seek dry rock under the conditions of physical parameter to obtain reasonable Vp for matrix density
Skeleton elasticity modulus (Kd).Kd is typically expressed as under Well-Log Acoustic Velocity:In formula, K is
The elasticity modulus of acquisition is calculated under known Vp and Vp/Vs relationship.
Using pressure as the function of Kd, and meet exponential relationship, so that it may which pressure is put into Gassmann relational expression.It can
To simulate more waves and prestack feature and observation attribute by the more waves and prestack characteristic simulation of oil, water layer and grease transition zone
Comparison obtain the quasi- of optimal simulation feature and observational characteristic to optimize petrophysical parameter, such as the parameter of rock matrix
It closes.
Step S103 utilizes calibrated petrophysical model, meter according to the fluid parameter and elastic parameter of reservoir model
Velocity of longitudinal wave, shear wave velocity and the density of reservoir model are calculated, then velocity of longitudinal wave is multiplied the impedance synthesized with density.?
After petrophysical model calibration, then it can pass through the oil of the acquisitions such as the elastic parameter of the fluid parameter of reservoir model and rock synthesis
Hide p-and s-wave velocity and the density etc. of model.It here can be by choosing the generated data of well point position and the speed of well logging actual measurement
Degree, density etc. carry out site-test analysis.Further include the velocity of longitudinal wave and density by synthesizing in the present embodiment, after this step, calculates
The impedance data synthesized out.
Step S104 determines that up-and-down boundary, wavelet and the earthquake of forward modeling of reservoir model are arrived using work area well-log information
Angular range and angle interval, utilize velocity of longitudinal wave, shear wave velocity and density, carry out the trace gather forward modeling of AVO prestack, generate AVO
Prestack trace gather forward modeling result.It is a part since reservoir model is in the longitudinal direction generally directed to some target zone or oily group
Information.Therefore, it is necessary to the velocity of longitudinal waves and density of the invalid grid in coboundary, lower boundary and the model to reservoir model
It is set.The acoustic wave parameter that the synthesis at bottom is pushed up in the acoustic wave parameter and model being normally set up is close.It is normally set up the phase of wavelet
Position and frequency are consistent with the wavelet of actual seismic data.
It further include according to the velocity of longitudinal wave, shear wave velocity, density and setting before the progress AVO prestack trace gather forward modeling
Fixed angular range and angle interval calculation reflection coefficient, the method for calculating reflection coefficient is mainly Zoeppritz equation, however
The accounting equation of Zoeppritz equation is more complicated, therefore the side of the calculating reflection coefficient of the Zoeppritz reduced form provided
Method has: boltfeld or Aki.
Step S105 compares the seismic data of AVO prestack trace gather forward modeling result and actual acquisition, complete if consistent
At the optimization of reservoir model, such as inconsistent modification is consistent.Since the dynamic parameter in reservoir model changes over time, because
This, the dynamic AVO analysis that can be changed over time.The AVO prestack trace gather of synthesis can be with the earthquake of well or actual acquisition
Data compare and analyze:
The first, the AVO trace gather for the well point position sometime put is analyzed, and is carried out with the production logging curve at time point
It is regular to analyze the fluid whether its AVO characteristic meets at well point for comparison；
The second, it synthesizes and is analyzed with the AVO of actual seismic acquisition time consistent AVO trace gather and actual acquisition seismic data
It compares, whether the AVO characteristic for analyzing the two is consistent.
When fluid rule at AVO characteristic and well point is not inconsistent or runs counter to the AVO characteristic of the seismic data of actual acquisition,
Reservoir model is then returned to, the modification of model is carried out, again history matching, to achieve the purpose that the meticulous depiction of reservoir model.
Fig. 3 to Figure 10 is with specific trial zone example, to illustrate the optimization method for building up of reservoir model of the invention
Process.
Specific example shares 156 mouthfuls of wells based on delta deposit, in test block, and area is about 3.6km2, such as Fig. 3 institute
Show, is established by seismic interpretation data, log analysis data, creation data and the core data of trial zone by history matching
The model profile figure of reservoir model, is analyzed by taking reservoir model as an example.Firstly, extracting interval of interest production logging data, shape
At petrophysical model template, as shown in Figure 4.After petrophysical model calibration, calculated using calibration result and reservoir model
Velocity of longitudinal wave, shear wave velocity and density, then velocity of longitudinal wave is multiplied the impedance synthesized with density, Fig. 5,6,7 are respectively
The calculated velocity of longitudinal wave of this Case Experiments On A area reservoir model, shear wave velocity, density.Fig. 8 is this example laboratory area by longitudinal wave speed
Degree and the calculated impedance of shear wave velocity.Fig. 9 and Figure 10 be respectively this Case Experiments On A area reservoir model setting up-and-down boundary and
Wavelet schematic diagram, that is, confirmed the relationship between dynamic parameter and acoustic wave parameter, synthesize entire reservoir model based on this
Acoustic wave parameter field, the acoustic wave parameter of the up-and-down boundary of setting model, and its son is obtained using the actual seismic data of this area
The phase and frequency of wave.Figure 11 is the angular range of the setting of the selected well point the L8-182 position in this Case Experiments On A area, angle
Interval and the schematic diagram for calculating reflection coefficient method set its angular range 0 and arrive as shown, choosing the well point L8-182 position
30 degree, 5 degree are wherein divided between angle, the method for calculating reflection coefficient is Aki approximate equation, and the PP prestack trace gather synthesized is such as
Shown in Figure 12, choose in 925m to 975m(Figure 12 shown in A frame) AVO analysis (shown in Figure 13) is carried out, pass through amplitude and incidence angle
Its slope of crossploting analysis be 405, intercept be -326, it can be deduced that amplitude with incidence angle (i.e. the distance of geophone offset) increase
And reduce, belong to first kind AVO.Analyze result be this section it is aqueous, by the actual well production logging curve (Figure 14) into
Row comparative analysis, the section belong to High water cut.Therefore, the analysis result of AVO be really consistent, it can be deduced that reservoir model
Accurately portray the conclusion of underground fluid situation.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects
Describe in detail it is bright, it should be understood that the above is only a specific embodiment of the present invention, the guarantor being not intended to limit the present invention
Range is protected, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in this
Within the protection scope of invention.
Claims (13)
1. a kind of reservoir model optimization method for changing road set analysis with geophone offset for amplitude, feature is through the following steps that reality
It is existing:
1) it is built according to work area seismic interpretation data, log analysis data, creation data and core data by history matching
Vertical reservoir model；
2) petrophysical model is demarcated using work area log analysis data and reservoir model data, establish seismic response parameter with
Non-linear relation between rock parameter；
3) reservoir model is calculated using calibrated petrophysical model according to the fluid parameter of reservoir model and elastic parameter
Velocity of longitudinal wave, shear wave velocity and density；Velocity of longitudinal wave is multiplied the impedance synthesized with density again；
4) using work area well-log information determine angular range that up-and-down boundary, wavelet and the earthquake of forward modeling of reservoir model are arrived and
Angle interval carries out the trace gather forward modeling of AVO prestack using velocity of longitudinal wave, shear wave velocity and density, generates the trace gather forward modeling of AVO prestack
As a result；
Before trace gather forward modeling, the acoustic wave parameter and density value surveyed at invalid grid using well-log information are to nothing in reservoir model
Velocity of longitudinal wave and density at effect grid are set；
Before trace gather forward modeling, counted according to the velocity of longitudinal wave, shear wave velocity, density and the angular range of setting and angle interval
Calculate reflection coefficient；
5) AVO prestack trace gather forward modeling result and the AVO analysis result of actual acquisition seismic data are compared, it is complete if consistent
It at the optimization of reservoir model, as inconsistent, modifies to the parameter of reservoir model, so that the AVO prestack trace gather forward modeling knot
Fruit is consistent with the AVO of actual acquisition seismic data analysis result.
2. the method according to claim 1, feature is that reservoir model data described in step 2) include: fluid parameter, elasticity ginseng
Number, rock properties, oil reservoir partitioned parameters, initialization calculating parameter, output control parameter, the essential information and life for simulating work
Produce parameter；Wherein:
Fluid parameter includes the ground surface density or severe of oil, gas and water, the reservoir volume ratio of oil, gas, viscosity, the volume ginseng of water
Number, the compressed coefficient；
Elastic parameter includes shear wave velocity, Lame constants and Poisson's ratio.
3. the method according to claim 1, feature, which is to establish reservoir model described in step 1), is built using seismic interpretation result
Vertical Skeleton Model establishes physical property model in conjunction with log data and core data, and history matching is carried out on the basis of creation data
Obtain reservoir model.
4. the method according to claim 1, feature is that petrophysical model described in step 2) is Voigt and Reuss theory mould
Type or Gassmann theoretical model.
5. the method according to claim 1, feature is established between seismic response parameter and rock parameter described in step 2)
Non-linear relation is the relationship established between acoustic wave parameter, elastic parameter and physical properties of rock, fluid parameter.
6. the method according to claim 1, feature is that petrophysical model described in step 3) is Voigt and Reuss theory mould
Type or Gassmann theoretical model；
Described calculate uses Gassmann theoretical model, and rock volume elasticity modulus K is calculated:
Wherein, KsFor rock particles elasticity modulus；KdFor dry elastic modulus of rock；For porosity；
KfFor elasticity of fluid modulus, calculated by following formula:
Wherein, Kw, Ko, KgIt is the elasticity modulus of water, oil, gas, S respectivelyw, SoRespectively
It is the saturation degree of water, oil, 1-Sw-SoFor gas saturation；
Velocity of longitudinal wave V of the sound wave in rockpIt is calculated by following formula:
Wherein, μ is Shear Modulus of Rock in Situ；ρ is rock density；Vs is shear wave velocity；
Wherein, ρw、ρo、ρgIt is that water, oil, gas are close respectively
Degree, ρmaFor matrix density；
Resultant impedance are as follows: AI=ρ Vp。
7. the method according to claim 1, feature is that the step 4) well-log information is density log curve and interval transit time
Log.
8. the method according to claim 1, feature, which is that step 4) is described, to be set as the density and speed assignment in well-log information
At the invalid grid in the reservoir model of corresponding depth, to avoid null value pair of the reservoir model in forward modeling in invalid grid
The influence of forward modeling result.
9. the method according to claim 1, feature is the phase of the step 4) wavelet and the son of frequency and actual seismic data
Wave phase is consistent with frequency.
10. the method according to claim 1, feature be the method that reflection coefficient is calculated described in step 4) include boltfeld or
Aki method.
11. the method according to claim 1, feature is that comparison described in step 5) is:
It analyzes and carries out in the AVO prestack trace gather forward modeling result of actual seismic acquisition time and the AVO of actual acquisition seismic data
Comparison.
12. the method according to claim 1, feature, which is described in step 5), is unanimously:
AVO characteristic meet the fluid at well point rule, the consistent indistinction of AVO characteristic, amplitude increase with the increase of geophone offset or
Person reduces, or the reversion of phase occurs, or,
The AVO of the AVO prestack trace gather forward modeling result of actual seismic acquisition time and actual acquisition seismic data analysis carries out pair
Than AVO characteristic is consistent.
13. the method according to claim 1, feature is that modification described in step 5) is to the porosity or shale in reservoir model
Content or net-gross ratio or oil saturation or pressure or temperature are modified.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410142188.3A CN104975851B (en) | 2014-04-10 | 2014-04-10 | For amplitude with the reservoir model optimization method of geophone offset variation road set analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410142188.3A CN104975851B (en) | 2014-04-10 | 2014-04-10 | For amplitude with the reservoir model optimization method of geophone offset variation road set analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104975851A CN104975851A (en) | 2015-10-14 |
| CN104975851B true CN104975851B (en) | 2019-01-01 |
Family
ID=54272767
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410142188.3A Active CN104975851B (en) | 2014-04-10 | 2014-04-10 | For amplitude with the reservoir model optimization method of geophone offset variation road set analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104975851B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105653815B (en) * | 2016-01-19 | 2019-03-26 | 中国海洋石油集团有限公司 | A kind of quantitative interpretation reservoir fluid location mode based on petrophysical model theory |
| US11567232B2 (en) | 2016-05-27 | 2023-01-31 | Total Sa | Method for providing a calibrated rock-physics model of a subsoil |
| CN110376645A (en) * | 2019-08-02 | 2019-10-25 | 中国石油大学(北京) | A kind of determination method, device and equipment of Thin Sandbody top interface location |
| CN113495293B (en) * | 2020-04-01 | 2023-09-26 | 中国石油天然气股份有限公司 | Reservoir fluid prediction method and device |
| CN113391373B (en) * | 2021-05-10 | 2023-12-12 | 三明学院 | Young modulus calculation method, device, equipment and storage medium |
| CN117572503B (en) * | 2024-01-15 | 2024-03-26 | 湖北神龙工程测试技术有限公司 | Rock wave velocity testing method |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8923093B2 (en) * | 2009-08-25 | 2014-12-30 | Westerngeco L.L.C. | Determining the quality of a seismic inversion |
| CN101872024B (en) * | 2010-06-02 | 2012-07-04 | 中国海洋石油总公司 | Method for carrying out well design by using time-lapse seismic |
| CN102305942B (en) * | 2011-07-15 | 2013-12-04 | 中国石油天然气集团公司 | Nonlinear AVO Fluid Discrimination Method Based on Three Parameters |
| CN103412336B (en) * | 2013-07-22 | 2015-12-09 | 中国石油天然气股份有限公司 | Longitudinal wave velocity prediction method for rock system in heterogeneous oil reservoir |
| CN103527184B (en) * | 2013-10-28 | 2016-04-20 | 北京大学 | A kind of Forecasting Methodology of dolostone reservoirs and system |
-
2014
- 2014-04-10 CN CN201410142188.3A patent/CN104975851B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104975851A (en) | 2015-10-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10983232B2 (en) | Seismic rock physics inversion method based on large area tight reservoir | |
| CA2895549C (en) | Fracturing and reactivated fracture volumes | |
| CN106855636B (en) | The prototype geological model Seismic forward method appeared based on carbonate reservoir | |
| CN104975851B (en) | For amplitude with the reservoir model optimization method of geophone offset variation road set analysis | |
| CN104516017B (en) | A kind of carbonate rock petrophysical parameter seismic inversion method | |
| CN105445791B (en) | A kind of formation pore pressure Forecasting Methodology based on a variety of seismic properties | |
| CN102636812B (en) | Method for obtaining storage space volume of carbonate reservoir | |
| CN104898161B (en) | Effective sandstone predicting method based on logging response simulator | |
| NO20180189A1 (en) | Seismic constrained discrete fracture network | |
| NO347906B1 (en) | Stratigraphic function | |
| CN105182424B (en) | Method and device for quantitatively predicting reservoir porosity based on patch saturation model | |
| WO2017035104A1 (en) | Velocity model seismic static correction | |
| US20140336940A1 (en) | Estimation of q-factor in time domain | |
| CA2940406A1 (en) | Characterizing a physical structure using a multidimensional noise model to attenuate noise data | |
| CN105607120A (en) | Time-shifting-logging-based method for building initial model with seismic facies constraint | |
| CN104675392A (en) | Reservoir lithology identification method based on pre-stack multi-parameter dimensionality reduction | |
| CN107728205B (en) | A kind of Formation pressure prediction method | |
| CN113031068A (en) | Reflection coefficient accurate base tracking prestack seismic inversion method | |
| Bigi et al. | Discrete fracture network of the Latemar carbonate platform | |
| CN104199088A (en) | Incident angle gather extraction method and system | |
| CN107422384B (en) | A method of vertical seismic profile data is imitated using well-log information | |
| US20120099396A1 (en) | System and method for characterization with non-unique solutions of anisotropic velocities | |
| CN109521470B (en) | Method for analyzing influence of geological structure on seismic inversion crack density | |
| CN109239777B (en) | Method for detecting development of tectonic coal by using joint inversion method | |
| CN104316959B (en) | Fluid identification based on equivalent fluid acoustic wave impedance |
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 |