CN102508291A - Transverse variable speed small scale body reflection coefficient formula application method - Google Patents
Transverse variable speed small scale body reflection coefficient formula application method Download PDFInfo
- Publication number
- CN102508291A CN102508291A CN2011103144896A CN201110314489A CN102508291A CN 102508291 A CN102508291 A CN 102508291A CN 2011103144896 A CN2011103144896 A CN 2011103144896A CN 201110314489 A CN201110314489 A CN 201110314489A CN 102508291 A CN102508291 A CN 102508291A
- Authority
- CN
- China
- Prior art keywords
- small scale
- scale body
- reflection coefficient
- formula
- lateral 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004088 simulation Methods 0.000 claims abstract description 35
- 238000004458 analytical method Methods 0.000 claims abstract description 19
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 18
- 238000013461 design Methods 0.000 claims description 8
- 241000209094 Oryza Species 0.000 claims description 7
- 235000007164 Oryza sativa Nutrition 0.000 claims description 7
- 235000009566 rice Nutrition 0.000 claims description 7
- 239000011435 rock Substances 0.000 claims description 6
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 239000004576 sand Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 238000012797 qualification Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 101100136092 Drosophila melanogaster peng gene Proteins 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 125000001183 hydrocarbyl group Chemical group 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a transverse variable speed small-scale volume reflection coefficient formula application method, which is characterized in that a double-layer medium model is established, and numerical simulation is carried out by utilizing a finite difference full wave field numerical simulation algorithm; AVO analysis is carried out on the obtained numerical simulation result; establishing a theoretical expression of a transverse variable speed small-scale body reflection coefficient formula; and on the basis of a theoretical expression, determining the undetermined coefficient according to the AVO analysis result of the obtained numerical simulation result, and finally obtaining the transverse variable speed small-scale body reflection coefficient formula. The method for solving the width and the elastic parameter of the transverse variable-speed small scale body through the formula is simple and practical, and can effectively describe the relation between the change of the reflection coefficient and the width and the elastic parameter of the scale body.
Description
Technical field
The present invention relates to the exploration of oil and gas field field, particularly a kind of lateral velocity variation small scale body reflection coefficient formula application process.
Background technology
Early 1980s Ostrander (1982) at first proposes to utilize reflection coefficient to change the identification gas sand with incident angle; He notices: the gas sand reflection amplitude increases with offset distance; The water bearing sand reflection amplitude reduces with the increase of offset distance, and the technology that hydro carbons detects has been enriched in the discovery of this phenomenon.The working mark of Ostrander is used the appearance of AVO (Amplitude Versus Offset) technology really; Over more than 30 year; The AVO technology constantly becomes the focus of people's research, and SEG (The Leading Edge) magazine was once repeatedly published special edition AVO theory, method, technology are discussed.Gassaway etc. are used for the inverting Poisson ratio to AVO (Amplitude Versus Offset) information, directly explain lithology and oil gas.
The AVO technology is to come a direct predicting oil and a technology estimating earth's crust lithologic parameter according to amplitude with the underground lithology that Changing Pattern reflected of geophone offset and the character of pore fluid thereof.Sun Peng is far away to be waited sorting out based on the approximate various reflection coefficient approximation methods of Zoeppritz and contrasting, and present most AVO research is being similar to the Zeoppriz equation under uniform dielectric or anisotropic medium prerequisite.But the reservoir heterogeneity of most of hydrocarbon-bearing pool of China is very strong, how to predict effectively that the nonuniformity reservoir is a key link of oil-gas exploration and development.Existing method for predicting reservoir is based on HORIZONTAL LAYERED MEDIUM WITH HIGH ACCURACY and uniform dielectric model mostly, can not satisfy the demand of heterogeneous reservoir prediction.Thereby, must carry out deep research from situation of change with the medium elastic parameter to the reflection coefficient of heterogeneous reservoir.
Under the situation of horizontal homogeneous or large scale lateral velocity variation, Zeoppritz equation and simplified style thereof can be described the reflection coefficient on stratum well.Under the situation of small scale lateral velocity variation; Whether existing reflection coefficient expression formula is suitable for; How the truth of underground medium could be described effectively? Current in the world seldom to the research of reflection coefficient approximate expression under the small scale lateral velocity variation situation; And mainly from the anisotropy angle, Zeoppritz equation approximate expression is revised, can not be applied among actual the AVO analysis and prestack inversion according to anisotropic each parameter.How horizontal small scale defines, how slewing range defines, and laterally small scale speed change body is much to the reflective index impacts degree, and laterally how approximate formula starts with under the small scale speed change situation, and these a series of problems all are badly in need of solving.
To the problems referred to above; Set up the small scale body of different scale; Change the speed in the small scale body,, the reflection coefficient of horizontal small scale speed change body is analyzed with all-wave field The Forward Modeling; Derive the relation of reflection coefficient and yardstick body size and lithologic parameter, promoted the development of heterogeneous reservoir forecasting techniques.
Summary of the invention
Fundamental purpose of the present invention is to solve the problem that exists in the prior art, and a kind of lateral velocity variation small scale body reflection coefficient formula application process is provided.
The objective of the invention is to be achieved through following technical proposals:
Lateral velocity variation small scale body reflection coefficient formula application process is characterized in that, comprises the steps:
Set up the two-layered medium model, in its layer dielectric, design the small scale body, and the width and the elastic parameter of this small scale body are set;
Seismic wavelet type and dominant frequency parameter are set, the wave detector exhaust position is set, carry out numerical simulation through finite difference all-wave Flow Field Numerical modeling algorithm;
The numerical simulation result that is obtained is carried out AVO to be analyzed;
Set up the theoretical expression of lateral velocity variation small scale body reflection coefficient formula: Δ R=aw Δ im+b Δ im+cw+d; Wherein, Δ R is the reflection coefficient changing value, and w is a yardstick body width, and Δ im is the impedance contrast of impedance of yardstick body and country rock, and a, b, c, d are undetermined coefficient;
On the basis of said theoretical expression, according to the AVO analysis result of the numerical simulation result that is obtained, use least square method and find the solution undetermined coefficient a, b, c, d in the said theoretical expression, obtain lateral velocity variation small scale body reflection coefficient formula.
Also include following reflection coefficient formula correction step:
Said two-layered medium model is provided with the width and the elastic parameter of different small scale bodies; And repeat the said step of asking for to undetermined coefficient a, b, c, d; Obtain new one group of expression formula undetermined coefficient a ', b ', c ', d ', said lateral velocity variation small scale body reflection coefficient formula is revised with one group of new expression formula undetermined coefficient a ', b ', c ', the d ' that is obtained.
Said seismic wavelet type is a Ricker wavelet, and its dominant frequency is 30Hz; Said wave detector is placed on from destination layer 50-100 rice.
Also comprise the steps:
Utilize the lateral velocity variation small scale body reflection coefficient formula that is obtained that small scale body width and elastic parameter are predicted.
Lateral velocity variation small scale body reflection coefficient formula application system is characterized in that, comprising: two-layered medium model building module, finite difference all-wave Flow Field Numerical simulation calculation module, AVO analysis module, theoretical expression are set up module, coefficient The Fitting Calculation module;
Said two-layered medium model building module is in order to set up the two-layered medium model; Design small scale body in the layer dielectric of said two-layered medium model; This small scale body comprises width and elastic parameter at least;
Said finite difference all-wave Flow Field Numerical simulation calculation module, in order to according to seismic wavelet type and dominant frequency parameter are set, and the wave detector exhaust position, carry out numerical simulation and obtain numerical simulation result;
Said AVO analysis module is analyzed in order to the numerical simulation result that is obtained is carried out AVO;
Said theoretical expression is set up module, in order to set up the theoretical expression of lateral velocity variation small scale body reflection coefficient formula: Δ R=aw Δ im+b Δ im+cw+d; Wherein, Δ R is the reflection coefficient changing value, and w is a yardstick body width, and Δ im is the impedance contrast of impedance of yardstick body and country rock, and a, b, c, d are undetermined coefficient;
Said coefficient The Fitting Calculation module; AVO analysis result in order to the numerical simulation result that obtained according to the AVO analysis module; Use least square method and find the solution undetermined coefficient a, b, c, d in the said theoretical expression; And, obtain lateral velocity variation small scale body reflection coefficient formula with this undetermined coefficient a, b, c, d substitution theoretical expression.
Also be provided with the formula correcting module;
Said formula correcting module; In order to the width and the elastic parameter of said two-layered medium model building module Small and Medium Sized body are made amendment, and lateral velocity variation small scale body reflection coefficient formula is revised according to one group of new expression formula undetermined coefficient a ', b ', c ', the d ' that said coefficient The Fitting Calculation module is obtained.
In the said finite difference all-wave Flow Field Numerical simulation calculation module, said seismic wavelet type is a Ricker wavelet, and its dominant frequency is 30Hz; Said wave detector is placed on from destination layer 50-100 rice.
Small scale body parameter prediction module also is set;
Said small scale body parameter prediction module is predicted small scale body width and elastic parameter in order to utilize the lateral velocity variation small scale body reflection coefficient formula that is obtained.
Pass through the embodiment of the invention; Can generate and make up lateral velocity variation small scale body reflection coefficient formula; The method of width and elastic parameter of asking for lateral velocity variation small scale body through this formula is simple and practical, can effectively describe reflection coefficient and change the relation with yardstick body width and elastic parameter.
Description of drawings
Accompanying drawing described herein is used to provide further understanding of the present invention, constitutes the application's a part, does not constitute qualification of the present invention.In the accompanying drawings:
Fig. 1 is a two-layered medium small scale body illustraton of model;
Fig. 2 be in the yardstick body during velocity variations reflection coefficient with the graph of relation of incident angle;
Fig. 3 be in the yardstick body during velocity variations reflection coefficient change the synoptic diagram that concerns with offset distance;
Fig. 4 be in the different in width yardstick body during velocity variations reflection coefficient change the synoptic diagram that concerns with offset distance;
Fig. 5 is that small scale body width and reflection coefficient difference concern synoptic diagram;
Fig. 6 is the application process process flow diagram of lateral velocity variation small scale body reflection coefficient formula;
Fig. 7 is the application system structural representation of lateral velocity variation small scale body reflection coefficient formula.
Embodiment
For making the object of the invention, technical scheme and advantage clearer,, the present invention is explained further details below in conjunction with embodiment and accompanying drawing.At this, exemplary embodiment of the present invention and explanation thereof are used to explain the present invention, but not as to qualification of the present invention.
In general, major part all is to carry out AVO with thick-layer homogeneous model to analyze now, and reflection coefficient can be asked for through approximate equation.In fact underground medium in the horizontal speed possibly have variation, if ignore it, just possibly cause and can not carry out accurately predicting reservoir.Forefathers are less for the reflection coefficient research of lateral velocity variation body, and therefore reflection coefficient approximate formula under the horizontal small scale speed change of the Shang Weiyou situation does not see ripe model yet.The present invention comes the variation of simulate formation medium velocity through design small scale body; The size that changes the small scale body is simulated the size of anomalous body; Realized simplified models like this, through add in the uniform dielectric yardstick body (as, wavelength, 1/2 wavelength, 1/4 wavelength or 1/8 wavelength etc.); Study the influence rule of lateral velocity variation to reflection coefficient through the full wave equation numerical simulation technology, the final generation makes up lateral velocity variation small scale body reflection coefficient formula.
Fig. 1 model medium parameter at the middle and upper levels is velocity of longitudinal wave 2192m/s, shear wave velocity 818m/s, density 2.160g/cm
3, the layer dielectric parameter is velocity of longitudinal wave 1543m/s, shear wave velocity 951.5m/s, density 1.973g/cm
3, small scale body width is a wavelength, the velocity of longitudinal wave value becomes 1550m/s, 1570m/s, 1580m/s, 1590m/s from 1543m/s successively in the yardstick body.
Fig. 2 is that reflection coefficient is with the relation curve of incident angle during velocity variations in the yardstick body, and as can be seen from the figure along with near the increase reflection coefficient bounce or jump small scale body of speed in the yardstick body, the absolute value of reflection coefficient is reducing.Be the influence of velocity variations in the further contrast small scale body to reflection coefficient; Do before changing and the residual error after changing; As shown in Figure 3, black, red, blue among the figure, purple is followed successively by the residual error of yardstick body reflection coefficient when being 1550m/s, 1570m/s, 1580m/s, 1590m/s and during 1543m/s rice, and as can be seen from the figure; Along with the increase of speed, the rate of change of reflection coefficient curve is in obvious increase.
The width that reduces the yardstick body successively is to 1/2 λ, 1/4 λ, and 1/8 λ, model medium parameter is identical with Fig. 1.In the yardstick body under the friction speed situation reflection coefficient basic consistent with the incident angle change curve with Fig. 2 trend, still owing to the reducing of yardstick body, reflection coefficient dies down at small scale body place.But still can find out the increase along with speed in the small scale body, near the absolute value of the reflection coefficient the small scale body position is reducing.Be better explanation small scale body to the influence of reflection coefficient, make equally in the yardstick body after the velocity variations that the result is as shown in Figure 4 with unchanged poor.From Fig. 4 a, b, c, contrast can find out that the yardstick body is more little, under the condition of identical filling velocity, more little to the influence of reflection coefficient.With width is that the yardstick body of 1/2 λ is an example, when yardstick is 1/2 λ the variable quantity maximal value within 0.003, and the yardstick body to be λ change maximal value about 0.005, reduce 60%, the size of visible yardstick body is more tangible to the variable effect of reflection coefficient.
With different scale body width is transverse axis, and reflection coefficient is changed to the longitudinal axis, projects to above-mentioned model value corresponding in the rectangular coordinate system, and the corresponding value of different impedances is being carried out linear fit (like Fig. 5).Can find out that from fitting formula changing with yardstick body width for the reflection coefficient that antibody is corresponding on the same group mutually has good linear relationship, definite coefficient minimum value of reflection fitting degree also reaches 99.6.%, and obviously its linear relation is extraordinary.Equally, for the yardstick body of same widths, the impedance contrast that its yardstick body is corresponding is also linear with reflection coefficient.Therefore, we recognize that the reflection coefficient variation that the small scale body causes is the corresponding width of small scale body and the function of respective impedance.Based on above-mentioned analysis and experiment, we have designed the application process of following lateral velocity variation small scale body reflection coefficient formula of the present invention.
Fig. 6 is the application process process flow diagram of this lateral velocity variation small scale body reflection coefficient formula.As shown in the figure, the application process of this lateral velocity variation small scale body reflection coefficient formula comprises the steps:
Step 1 is set up the two-layered medium model, design small scale body in its layer dielectric, and width and the elastic parameter (said elastic parameter is the broad sense title of multiple parameter types such as the velocity of longitudinal wave that comprises medium, shear wave velocity, density) of this small scale body are set;
Step 2 is provided with seismic wavelet type and dominant frequency parameter, and the wave detector exhaust position is set, and carries out numerical simulation through finite difference all-wave Flow Field Numerical modeling algorithm;
Step 4, set up the theoretical expression of lateral velocity variation small scale body reflection coefficient formula:
ΔR=awΔim+bΔim+cw+d
Wherein, Δ R is the reflection coefficient changing value, and w is a yardstick body width, and Δ im is the impedance contrast of impedance of yardstick body and country rock, and a, b, c, d are undetermined coefficient;
Step 5; On the basis of above-mentioned theory expression formula; According to the AVO analysis result of the numerical simulation result that step 3 obtained, use least square method and find the solution undetermined coefficient a, b, c, d in the said theoretical expression, obtain lateral velocity variation small scale body reflection coefficient formula.
Through above-mentioned reflection coefficient formula applying step, we have obtained lateral velocity variation small scale body reflection coefficient formula, and then utilize this formula can predict the width and the elastic parameter of lateral velocity variation small scale body, for AVO analyzes and the AVO inverting provides solid foundation.The method of width and elastic parameter of asking for lateral velocity variation small scale body through this formula is simple and practical, can effectively describe reflection coefficient and change the relation with yardstick body width and elastic parameter.This method has been filled up the blank that the AVO technology is studied under horizontal small scale speed change situation.
The inventor has carried out goodness of fit detection to the lateral velocity variation small scale body reflection coefficient formula that said method obtained.Calculate the goodness of fit through the goodness of fit and confirm that coefficient is 0.99.Thereby the fitting effect that can confirm the reflection coefficient formula that said method obtains is extraordinary.So just can prove that utilizing this formula is reasonable, believable to the reflection coefficient amplitude of variation calculated value that different scale body width and impedance cause.
In order to make the reflection coefficient formula that obtains through said method more accurately more effective, the present invention has also designed the step of as follows the reflection coefficient formula of above-mentioned acquisition being revised:
Step 6; The two-layered medium model of step 1 is provided with the width and the elastic parameter of different small scale bodies; And repeat new one group of expression formula undetermined coefficient a ', b ', c ', the d ' of above-mentioned steps 1 to step 5 acquisition, with one group of new expression formula undetermined coefficient a ', b ', c ', the d ' that is obtained lateral velocity variation small scale body reflection coefficient formula is revised.
We can revise this lateral velocity variation small scale body reflection coefficient formula based on the small scale body of limited group of different parameters combination repeatedly through this step, thereby obtain computing formula more accurately and effectively.
In addition, in order to improve the resolution of analog computation, to seismic wavelet type and the wavelet dominant frequency of being selected for use in the said step 2, and the wave detector exhaust position that is designed also has requirement.In the present embodiment, what specifically adopt is thunder gram (Ricker) wavelet, and its dominant frequency is 30Hz.And in order to eliminate whole the extending influence of seismic event, wave detector is placed on from destination layer 50-100 rice.Can guarantee that so promptly being recorded into ejected wave writes down reflection wave again, estimate that the slowness direction vector of each ripple calculates the AVO curve then.
In addition, this lateral velocity variation small scale body reflection coefficient formula application process can also comprise the steps:
Step 7 utilizes the lateral velocity variation small scale body reflection coefficient formula that is obtained that small scale body width and elastic parameter are predicted.
Through this step, we can utilize the width and the elastic parameter of the prediction lateral velocity variation small scale body that obtains that calculates, for AVO analyzes and the AVO inverting provides solid foundation.
Fig. 7 is the structural representation of the application system of this lateral velocity variation small scale body reflection coefficient formula.As shown in the figure; The application system of this lateral velocity variation small scale body reflection coefficient formula comprises: two-layered medium model building module, finite difference all-wave Flow Field Numerical simulation calculation module, AVO analysis module, theoretical expression are set up module, coefficient The Fitting Calculation module.
Said two-layered medium model building module is in order to set up the two-layered medium model; Design small scale body in the layer dielectric of said two-layered medium model; This small scale body comprises width and elastic parameter at least;
Said finite difference all-wave Flow Field Numerical simulation calculation module, in order to according to seismic wavelet type and dominant frequency parameter are set, and the wave detector exhaust position, carry out numerical simulation and obtain numerical simulation result;
Said AVO analysis module is analyzed in order to the numerical simulation result that is obtained is carried out AVO;
Said theoretical expression is set up module, in order to set up the theoretical expression of lateral velocity variation small scale body reflection coefficient formula:
ΔR=awΔim+bΔim+cw+d
Wherein, Δ R is the reflection coefficient changing value, and w is a yardstick body width, and Δ im is the impedance contrast of impedance of yardstick body and country rock, and a, b, c, d are undetermined coefficient;
Said coefficient The Fitting Calculation module; AVO analysis result in order to the numerical simulation result that obtained according to the AVO analysis module; Use least square method and find the solution undetermined coefficient a, b, c, d in the said theoretical expression; And, obtain lateral velocity variation small scale body reflection coefficient formula with this undetermined coefficient a, b, c, d substitution theoretical expression.
In order to make the reflection coefficient formula that obtains through said method more accurately more effective, in the reflection coefficient formula application system that the present invention designed, also be provided with the formula correcting module.
Said formula correcting module; In order to the width and the elastic parameter of said two-layered medium model building module Small and Medium Sized body are made amendment, and lateral velocity variation small scale body reflection coefficient formula is revised according to one group of new expression formula undetermined coefficient a ', b ', c ', the d ' that said coefficient The Fitting Calculation module is obtained.
We can revise this lateral velocity variation small scale body reflection coefficient formula based on the small scale body of limited group of different parameters combination repeatedly through this formula correcting module, thereby obtain computing formula more accurately and effectively.
In addition, in order to improve the resolution of analog computation, in the said finite difference all-wave Flow Field Numerical simulation calculation module, what specifically adopt is thunder gram (Ricker) wavelet, and its dominant frequency is 30Hz.And in order to eliminate whole the extending influence of seismic event, said wave detector is placed on from destination layer 50-100 rice.
In addition, in the reflection coefficient formula application system that the present invention designed, small scale body parameter prediction module is set also.
Said small scale body parameter prediction module is predicted small scale body width and elastic parameter in order to utilize the lateral velocity variation small scale body reflection coefficient formula that is obtained.
In sum, the present invention studies the influence rule of lateral velocity variation to reflection coefficient through the full wave equation numerical simulation technology, and final the generation makes up lateral velocity variation small scale body reflection coefficient formula.The method of width and elastic parameter of asking for lateral velocity variation small scale body through this formula is simple and practical, can effectively describe reflection coefficient and change the relation with yardstick body width and elastic parameter.Persons skilled in the art any not creative transformation of doing under this design philosophy all should be regarded as within protection scope of the present invention.
Claims (8)
1. lateral velocity variation small scale body reflection coefficient formula application process is characterized in that, comprises the steps:
Set up the two-layered medium model, in its layer dielectric, design the small scale body, and the width and the elastic parameter of this small scale body are set;
Seismic wavelet type and dominant frequency parameter are set, the wave detector exhaust position is set, carry out numerical simulation through finite difference all-wave Flow Field Numerical modeling algorithm;
The numerical simulation result that is obtained is carried out AVO to be analyzed;
Set up the theoretical expression of lateral velocity variation small scale body reflection coefficient formula: Δ R=aw Δ im+b Δ im+cw+d; Wherein, Δ R is the reflection coefficient changing value, and w is a yardstick body width, and Δ im is the impedance contrast of impedance of yardstick body and country rock, and a, b, c, d are undetermined coefficient;
On the basis of said theoretical expression, according to the AVO analysis result of the numerical simulation result that is obtained, use least square method and find the solution undetermined coefficient a, b, c, d in the said theoretical expression, obtain lateral velocity variation small scale body reflection coefficient formula.
2. lateral velocity variation small scale body reflection coefficient formula application process as claimed in claim 1 is characterized in that, also includes following reflection coefficient formula correction step:
Said two-layered medium model is provided with the width and the elastic parameter of different small scale bodies; And repeat the said step of asking for to undetermined coefficient a, b, c, d; Obtain new one group of expression formula undetermined coefficient a ', b ', c ', d ', said lateral velocity variation small scale body reflection coefficient formula is revised with one group of new expression formula undetermined coefficient a ', b ', c ', the d ' that is obtained.
3. lateral velocity variation small scale body reflection coefficient formula application process as claimed in claim 1 is characterized in that said seismic wavelet type is a Ricker wavelet, and its dominant frequency is 30Hz; Said wave detector is placed on from destination layer 50-100 rice.
4. lateral velocity variation small scale body reflection coefficient formula application process as claimed in claim 1 is characterized in that, also comprises the steps:
Utilize the lateral velocity variation small scale body reflection coefficient formula that is obtained that small scale body width and elastic parameter are predicted.
5. lateral velocity variation small scale body reflection coefficient formula application system; It is characterized in that, comprising: two-layered medium model building module, finite difference all-wave Flow Field Numerical simulation calculation module, AVO analysis module, theoretical expression are set up module, coefficient The Fitting Calculation module;
Said two-layered medium model building module is in order to set up the two-layered medium model; Design small scale body in the layer dielectric of said two-layered medium model; This small scale body comprises width and elastic parameter at least;
Said finite difference all-wave Flow Field Numerical simulation calculation module, in order to according to seismic wavelet type and dominant frequency parameter are set, and the wave detector exhaust position, carry out numerical simulation and obtain numerical simulation result;
Said AVO analysis module is analyzed in order to the numerical simulation result that is obtained is carried out AVO;
Said theoretical expression is set up module, in order to set up the theoretical expression of lateral velocity variation small scale body reflection coefficient formula: Δ R=aw Δ im+b Δ im+cw+d; Wherein, Δ R is the reflection coefficient changing value, and w is a yardstick body width, and Δ im is the impedance contrast of impedance of yardstick body and country rock, and a, b, c, d are undetermined coefficient;
Said coefficient The Fitting Calculation module; AVO analysis result in order to the numerical simulation result that obtained according to the AVO analysis module; Use least square method and find the solution undetermined coefficient a, b, c, d in the said theoretical expression; And, obtain lateral velocity variation small scale body reflection coefficient formula with this undetermined coefficient a, b, c, d substitution theoretical expression.
6. lateral velocity variation small scale body reflection coefficient formula application system as claimed in claim 5 is characterized in that, also is provided with the formula correcting module;
Said formula correcting module; In order to the width and the elastic parameter of said two-layered medium model building module Small and Medium Sized body are made amendment, and lateral velocity variation small scale body reflection coefficient formula is revised according to one group of new expression formula undetermined coefficient a ', b ', c ', the d ' that said coefficient The Fitting Calculation module is obtained.
7. lateral velocity variation small scale body reflection coefficient formula application system as claimed in claim 5 is characterized in that in the said finite difference all-wave Flow Field Numerical simulation calculation module, said seismic wavelet type is a Ricker wavelet, and its dominant frequency is 30Hz; Said wave detector is placed on from destination layer 50-100 rice.
8. lateral velocity variation small scale body reflection coefficient formula application system as claimed in claim 5 is characterized in that, small scale body parameter prediction module also is set;
Said small scale body parameter prediction module is predicted small scale body width and elastic parameter in order to utilize the lateral velocity variation small scale body reflection coefficient formula that is obtained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110314489.6A CN102508291B (en) | 2011-10-17 | 2011-10-17 | Transverse variable speed small scale body reflection coefficient formula application method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110314489.6A CN102508291B (en) | 2011-10-17 | 2011-10-17 | Transverse variable speed small scale body reflection coefficient formula application method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102508291A true CN102508291A (en) | 2012-06-20 |
CN102508291B CN102508291B (en) | 2014-01-15 |
Family
ID=46220396
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110314489.6A Active CN102508291B (en) | 2011-10-17 | 2011-10-17 | Transverse variable speed small scale body reflection coefficient formula application method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102508291B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104820235A (en) * | 2015-05-07 | 2015-08-05 | 信真维超能源科技(北京)有限公司 | Method for solving velocity reflectivity and density reflectivity |
CN107356972A (en) * | 2017-06-28 | 2017-11-17 | 中国石油大学(华东) | A kind of imaging method of anisotropic medium |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0175321A2 (en) * | 1984-09-17 | 1986-03-26 | Phillips Petroleum Company | Seismic data processing |
SU1408395A1 (en) * | 1986-03-04 | 1988-07-07 | Киевское геофизическое отделение Украинского научно-исследовательского геологоразведочного института | Seismic prospecting method |
WO1997041456A1 (en) * | 1996-04-29 | 1997-11-06 | The Trustees Of Columbia University In The City Of New York | Method for inverting reflection trace data from 3-d and 4-d seismic surveys |
CN102116869A (en) * | 2011-02-12 | 2011-07-06 | 中国石油大学(华东) | High-precision prestack domain least square migration seismic imaging technology |
-
2011
- 2011-10-17 CN CN201110314489.6A patent/CN102508291B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0175321A2 (en) * | 1984-09-17 | 1986-03-26 | Phillips Petroleum Company | Seismic data processing |
SU1408395A1 (en) * | 1986-03-04 | 1988-07-07 | Киевское геофизическое отделение Украинского научно-исследовательского геологоразведочного института | Seismic prospecting method |
WO1997041456A1 (en) * | 1996-04-29 | 1997-11-06 | The Trustees Of Columbia University In The City Of New York | Method for inverting reflection trace data from 3-d and 4-d seismic surveys |
CN102116869A (en) * | 2011-02-12 | 2011-07-06 | 中国石油大学(华东) | High-precision prestack domain least square migration seismic imaging technology |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104820235A (en) * | 2015-05-07 | 2015-08-05 | 信真维超能源科技(北京)有限公司 | Method for solving velocity reflectivity and density reflectivity |
CN104820235B (en) * | 2015-05-07 | 2017-08-29 | 信真维超能源科技(北京)有限公司 | The decomposition method of speed reflectivity and density reflectivity |
CN107356972A (en) * | 2017-06-28 | 2017-11-17 | 中国石油大学(华东) | A kind of imaging method of anisotropic medium |
CN107356972B (en) * | 2017-06-28 | 2019-09-06 | 中国石油大学(华东) | A kind of imaging method of anisotropic medium |
Also Published As
Publication number | Publication date |
---|---|
CN102508291B (en) | 2014-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103454685B (en) | Method and device for predicting sand body thickness by utilizing logging constrained wave impedance inversion | |
CA2920499C (en) | Stratigraphic function | |
CN105353412A (en) | Calculating method and system of well-to-seismic integration average speed field | |
CN103245970B (en) | Pre-stack seismic wide angle retrieval method | |
CN104597490A (en) | Multi-wave AVO reservoir elastic parameter inversion method based on precise Zoeppritz equation | |
CN103733089B (en) | For including the system and method that the underground of uncertainty estimation characterizes | |
CN104142516B (en) | Method for predicting thickness of thin single sand bed | |
CN102841378B (en) | Method for predicting reservoir stratum by seismic inversion data | |
US9952341B2 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
CN102116869A (en) | High-precision prestack domain least square migration seismic imaging technology | |
CN102109611A (en) | Fast and convenient method for predicting high-quality petroleum reservoir in virtue of seism attributes | |
CN102455438A (en) | Method for predicting volume of carbonate rock fractured cave type reservoir | |
CN105607120A (en) | Time-shifting-logging-based method for building initial model with seismic facies constraint | |
CN104155687A (en) | Phase control post-stack acoustic wave impedance inversion method | |
CN105372705A (en) | Stratum slicing method based on multi-wave data | |
CN104570066A (en) | Construction Method of Seismic Inversion Low-Frequency Model | |
CN106920176A (en) | Mining area scale mineral resource estimation method and system | |
CN102156297A (en) | Fluid substitution method based on sandstone reservoir post-stack seismic data | |
CN109188520A (en) | Thin reservoir thickness prediction method and device | |
CN103245972A (en) | Method for determining complex geologic structure in two-dimensional space | |
CN105242307A (en) | Complex carbonate stratum earthquake porosity obtaining method and apparatus | |
CN112147684A (en) | Method and device for characterizing activity intensity of co-deposition fault | |
CN102866422B (en) | A kind of depth domain geological entity model generation method | |
CN104216013B (en) | The method of C3 coherent body based on wide-azimuth data | |
CN102508291B (en) | Transverse variable speed small scale body reflection coefficient formula application method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |