1、AVO地震资料的解释IntroductionIn our previous section on rock physics we discussed fluid effects on P and S-wave velocity,and density.We then looked at post-stack inversion,its advantages and limitations.Next,we considered the recording of S-wave data.In this section,we will first review the basic principle
2、s of AVO and show its relationship with P and S-wave velocity.We will then look at the AVO response of two simple models,one wet and one gas-saturated.We will then look at various AVO modeling schemes,including full wave equation modeling and anisotropic modeling.“Bright spots”Recall that in the sec
3、tion on inversion,we showed the“bright spot”shown above,and pointed out that in the 1970s this would have been interpreted as a gas sand.3The AVO methodBut“bright spots”can be caused by lithologic variations as well as gas sands.This lead geophysicists in the 1980s to start looking at pre-stack seis
4、mic data.The amplitude increase with offset shown here is an example of a Class 3 sand,as we will discuss later.4Reflected P-wave=RP(q1)Reflected SV-wave=RS(q1)Transmitted P-wave=TP(q1)Incident P-waveTransmitted SV-wave=TS(q1)VP1,VS1,r1VP2,VS2,r2q11q1q22Consider the interface between two geologic ho
5、rizons of differing P and S-wave velocity and density and an incident P-wave at angle q q1.This will produce P and S reflected and transmitted waves,as shown above.Mode Conversion5But how do we utilize mode conversion?There are actually two ways:(1)Record the converted S-waves using three-component
6、receivers(in the X,Y and Z directions).This was discussed in the last chapter.(2)Interpret the amplitudes of the P-waves as a function of offset,or angle,which contain implied information about the S-waves.This is called the AVO(Amplitude versus Offset)method.In the AVO method,we can make use of the
7、 Zoeppritz equations,or some approximation to these equations,to extract S-wave type information from P-wave reflections at different offsets.Utilizing Mode Conversion6The Zoeppritz EquationsZoeppritz derived the amplitudes of the reflected and transmitted waves using the conservation of stress and
8、displacement across the layer boundary,which gives four equations with four unknowns.Inverting the matrix form of the Zoeppritz equations gives us the exact amplitudes as a function of angle:111112112221122111122111221221112221111221122112cos2sincossin2sin2cos2sin2cos2cos2cos2cos2sinsincossincoscoss
9、incossinqqqrrrrrrrrqqqqqPSPPPSSPSPSPSSPSPSPVVVVVVVVVVVVVVVTTRRThe Aki-Richards EquationThe is a linearized approximation to the.The initial form(Richards and Frasier,1976)separated the velocity and density terms:SSPPPVVcbVVaRrrq)(where:,sin4,sin25.0,cos2122222qqqPSPSVVcVVba.2,2,2,2121212121212tiSSSS
10、SSPPPPPPandVVVVVVVVVVVVqqqrrrrrrWiggins version of Aki-RichardsA totally equivalent form was derived by Wiggins.He separated the equation into three reflection terms,as follows:qqqq222sintansin)(CBARPwhere:curvature.21gradient,2421intercept,21220pPPSSSPSpPpPPVVCVVVVVVVVBVVRArrrrTo see why the A term
11、 in the linearized approximation is approximately equal to the zero-offset reflection coefficient,recall that in Part 2 we showed that:Reflectivity approximation.,212120PPPPPPPPPVZZZ ZZ ZZ=RrThis leads to:.212ln2ln2ln0rrrPPPPPVVVZR10.and,21,sin2tan21,sin8,tan1:where ,)(02222222200rrrrqqqqqDSSSPSPSDS
12、PRVVRVVfVVedfReRdRRFattis version of Aki-RichardsAnother equivalent form was derived by Fatti et al.They separated the equation into three reflection terms,as follows:A Summary of the Aki-Richards Eq.All three forms of the Aki-Richards equation consist of the sum of three terms,each term consisting
13、of a weight multiplied by an elastic parameter(i.e.a function of VP,VS or r).Here is a summary:Note that the weighting terms a,b,c and d,e,f contain the squared VP/VS ratio as well as q.However,in the Wiggins et al.formulation,this term is in the elastic parameter B.cba,EquationWeightsElastic Parame
14、tersAki-RichardsWiggins et al.Fatti et al.SSPPVVVV,rrqqq222sintan,sin,1CBA,00DSPRRRfed,12A physical interpretation of the three equations is as follows:(1)Since the seismic trace consists of changes in impedance rather than velocity or density independently,the original form of the Aki-Richards equa
15、tion is rarely used.(2)The A,B,C formulation of the Aki-Richards equation is very useful for extracting empirical information about the AVO effect(i.e.A,which is called the intercept,B,called the gradient,and C,called the curvature)which can then be displayed or cross-plotted.Explicit information ab
16、out the Vp/Vs ratio is not needed in the weights.(3)The Fatti et al.formulation gives us a way to extract quantitative information about the P and S reflectivity which can then be used for pre-stack inversion.The terms RP0 and RS0 are the linearized zero-angle P and S-wave reflection coefficients.Ph
17、ysical interpretation13Let us now see how to get from the geology to the seismic using the second two forms of the Aki-Richards equation.We will do this by using the two models shown below.Model A consists of a wet,or brine,sand,and Model B consists of a gas-saturated sand.Wet and Gas Models(a)Wet m
18、odel(b)Gas modelVP1,VS1,r r1VP2,VS2,r r2VP1,VS1,r r1VP2,VS2,r r214In the section on rock physics,we computed values for wet and gas sands using the Biot-Gassmann equations.The computed values were:Wet:VP2=2500 m/s,VS2=1250 m/s,r2=2.11 g/cc,s2 =0.33 Gas:VP2=2000 m/s,VS2=1310 m/s,r2=1.95 g/cc,s2=0.12
19、Values for a typical shale are:Shale:VP1=2250 m/s,VS1=1125 m/s,r1=2.0 g/cc,s1=0.33This gives us the following values at the top of the sand/shale zone:Model Values25.0 and ,054.0,105.0,105.02PSSSPPVVVVVVrrWet Sand:328.0 and ,025.0,152.0,118.02PSSSPPVVVVVVrrGas Sand:15Exercise 4-1Compute the intercep
20、t,A,the gradient term,B,and the curvature,C,for the top and base of both the wet model and the gas model,given the parameters on the previous page.Note that the base values are the negative of the top due to symmetry:Wet Model Top:A=B=C=Gas Model Top:A=B=C=Wet Model Base:A=B=C=Gas Model Base:A=B=C=E
21、xercise 4-2For both the next exercise and a later exercise on anisotropy effects in AVO,you will need to compute the following trigonometric functions of four angles.Compute them here for later use.Angle0153045sin2qtan2qsin2q*tan2qExercise 4-3Using the A,B,and C terms for the wet and gas models that
22、 were computed in Exercise 1,work out the values for R(q)at angles of 0o,15o,30o,and 45oin the table below.Then,plot the results on graph paper,with and without the third term,as a function of q.Top:Base:Wet ModelGas ModelAngle(degrees)R(q)R(q)-term3R(q)R(q)-term30153045Wet ModelGas ModelAngle(degre
23、es)R(q)R(q)-term3R(q)R(q)-term30153045Wet Model AVO Curves0.0RP(q)q0.10.2-0.1-0.210o20o30o40o50o19Gas Model AVO Curvesq0.0RP(q)0.10.2-0.1-0.210o20o30o40o50o20Wet Model AVO Curvesq0.0RP(q)10o20o30o40o50o0.10.2-0.1-0.2Wet Sand BaseWet Sand Top21Gas Model AVO Curvesq0.0RP(q)10o20o30o40o50o0.10.2-0.1-0.
24、2Gas Sand TopGas Sand Base22This figure on the right shows the computed AVO curves for the top and base interfaces of the gas sand using all three terms(A,B,and C)in the Aki-Richards equation,and then only the first two terms(A and B).Note the deviation of the two above 25 degrees.Gas Model AVO Curv
25、esThis figure on the right shows the computed AVO curves for the top and base interfaces of the wet sand using all three terms(A,B,and C)in the Aki-Richards equation,and then only the first two terms(A and B).Note the deviation of the two above 25 degrees.Wet Model AVO CurvesParameters for ABC and F
26、atti eqs.Her is a comparison of the results from the ABC and Fatti equations for the top of the sands(because of symmetry in this example,the base of sand values are simply these values multiplied by-1):025.0,063.0,071.0059.0,242.0,071.000DSPRRRCBAGas Sand:Wet Sand:053.0,079.0,079.0054.0,079.0,079.0
27、00DSPRRRCBANote that A and B have the same polarity for the gas sand and opposite polarity for the wet sand,whereas RP0 and RS0 have opposite polarity for the gas sand and the same polarity for the wet sand.The reason for this will be clear later.25Aki-Richards valuesHere are computed values for the
28、 ABC and Fatti versions of the Aki-Richards equation at angles of 0,30 and 60 degrees:0o Gas30o Gas60o Gas0o Wet30o Wet60o Wet1st Term 2nd Term 3rd Term2nd Term 3rd Term1st TermRP(q)ABC MethodFatti MethodAngle/Sand-0.0710.079-0.0710.07900000000-0.0710.079-0.1370.064-0.071-0.0710.0790.079-0.060-0.006
29、-0.095-0.0426x10-5-0.0200.0050.106-0.040-0.002-0.181-0.385-0.133-0.285-0.1250.025-0.0600.1190.318-0.119-0.0610.13826There was a lot of information in the last slide,but the key points are:(1)The individual terms in each approach are different,but the sum is always identical.(2)For an angle of zero d
30、egrees,the second two terms in both methods are equal to zero,and the scalar on the first term in the Fatti method is equal to one.(3)In the ABC method,the first term is always the zero offset reflection coefficient,but this is true only at zero angle in the Fatti method.(4)The third term makes less
31、 of a contribution to the sum in the Fatti method than in the ABC method.(5)The next slides will show the results at all angles.Summary of the ABC and Fatti methods27This figure on the right shows the AVO curves computed using the Zoeppritz equations and the two and three term ABC equation,for the g
32、as sand model.Notice the strong deviation for the two term versus three term sum.Note:On the next four plots,the curves have been calculated as a function of incident angle and scaled to average angle.Zoeppritz vs ABC Gas SandZoeppritzABC method:two termABC method:three term28This figure on the righ
33、t shows the AVO curves computed using the Zoeppritz equations and the two and three term ABC equation,for the wet sand model.Again,notice the strong deviation for the two term versus three term sum.Zoeppritz vs ABC Wet SandZoeppritzABC method:two termABC method:three term29This figure on the right s
34、hows the AVO curves computed using the Zoeppritz equations and the two and three term Fatti equation,for the gas sand model.Notice there is less deviation between the two term and three term sum than with the ABC approach.Zoeppritz vs Fatti Gas SandZoeppritzFatti method:two termFatti method:three te
35、rm30This figure on the right shows the AVO curves computed using the Zoeppritz equations and the two and three term Fatti equation,for the wet sand model.As in the gas sand case,there is less deviation between the two term and three term sum than with the ABC approach.Zoeppritz vs Fatti Wet SandZoep
36、pritzFatti method:two termFatti method:three term31This final computed synthetic seismogram is shown above on the right,where the log curves are on the left.Notice that the sand is thin enough that the wavelets from the top and bottom of the layer“tune”together.The final synthetic seismogram32Ostran
37、ders Paper(1984)was one of the first to write about in gas sands and proposed a simple two-layer model which encased a low impedance,low Poissons ratio sand,between two higher impedance,higher Poissons ratio shales.This model is shown in the next slide.Ostranders model worked well in the Sacramento
38、valley gas fields.However,it represents only one type of (Class 3)and the others will be discussed in the next section.Ostranders ModelOstrander(1984)wrote the classic paper on AVO.His model is shown above.Notice that the model consists of a low acoustic impedance gas sand encased between two shales
39、34Synthetic from Ostranders Model(a)Well log responses for the model.(b)Synthetic seismic.Notice that Ostranders model produces an increase in amplitude on the pre-stack synthetic gather.35AVO Curves from Ostrander(a)Response from top of model to 45o.(Note that the transmitted P-wave amplitude is sh
40、ifted to plot within the data range).(b)Response from base of model to 45o.36Ostranders case study-stackOstranders case study is from the Sacramento basin.The stack above has“bright spots”at locations A,B,and C,but only A and B are due to gas.37Ostranders case study Supergathers from locations A,B,a
41、nd C.Note that locations A and B show amplitude increases with offset but C does not.(A)(B)(C)38Shueys EquationShuey(1985)rewrote the Aki-Richards equation using VP,r r,and s s,writing the basic form the same way:2)1(121)1(2ssssDDAB12122,/:ssssssrrPPPPVVVVDwhereOnly the gradient is different than in
42、 the Aki-Richards expression,and is given by:qqqq222sintansin)(CBARShuey vs Aki-Richards In this course,we have been using a modeled gas sand and wet sand example.Using Shueys equation for this example,we get the following comparison with the answers in exercise 6-1:B(Aki-Richards)B(Shuey)Gas Sand T
43、op:-0.242-0.252 Wet Sand Top:-0.079-0.079Why do we get the same values in the wet case but not in the gas case?Gas Sand ModelAki-Richards vs Shuey-0.250-0.200-0.150-0.100-0.0500.0000.0500.1000.1500.2000.250051015202530354045Angle(degrees)AmplitudeA-R TopShuey TopA-R BaseShuey BaseThis figure showsa
44、comparisonbetween the Aki-Richards and Shuey equationsfor the gas sand we just considered.Shuey vs Aki-Richards Multi-layer AVO modelingMulti-layer modeling in the consists first of creating a stack of N layers,generally using well logs,and defining the thickness,P-wave velocity,S-wave velocity,and
45、density for each layer.You must then decide what effects are to be included in the model:primaries only,converted waves,multiples,or some combination of these.42The Possible Modelled EventsThe following example,taken from Simmons and Backus(AVO Modeling and the locally converted shear wave,Geophysic
46、s 59,p1237,August,1994),illustrates the effect of wave equation modeling.The figure above shows the modelling options.43The Oil Sand ModelSimmons and Backus used the thin bed oil sand model shown above.44Response to various algorithms(A)Primaries-only Zoeppritz,(B)+single leg shear,(C)+double-leg sh
47、ear,(D)+multiples,(E)Wave equation solution,(F)Linearized approximation.Simmons and Backus(1994)45Primary and Converted WavesSimmons and Backus(1994)Zoeppritz,primaries onlyZoeppritz,primaries+single leg conversionsAki-Richards,primariesSingle leg conversions46Logs from a real data exampleThe logs s
48、hown above come from a real data example in the Colony sand that we will look at in the next section.Models from a Real Data Example(a)Full elastic wave.(b)Zoeppritzequation.(c)Aki-Richardsequation.AVO ModelingBased on AVO theory and the rock physics of the reservoir,we can perform AVO modeling,as s
49、hown above.Note that the model result is a fairly good match to the offset stack.P-waveDensityS-wavePoissons ratioSyntheticOffset StackAnisotropy and AVOSo far,we have considered only the isotropic case,in which earth parameters such as velocity do not depend on seismic propagation angle.In the next
50、 few slides,we will discuss anisotropy,in particular the case of Transverse Isotropy with a vertical symmetry axis,or VTI.We will then see how anisotropy affects the AVO response.Finally,we will look at this effect on our original modelIsotropic versus anisotropic velocityAs mentioned,in an isotropi
