Regularization of the AVO inverse problem by means of a ...
Regularization of the AVO inverse problem by means of a ...
Regularization of the AVO inverse problem by means of a ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
CHAPTER 2. <strong>AVO</strong> MODELING 14<br />
Employing equations (2.5) in equation (2.20) and applying differentiation, we have<br />
2ρ1β 2 ω<br />
1p cos θi<br />
A0 − 2ρ1β<br />
α1<br />
2 1p cos θi<br />
α1<br />
ω<br />
A1 − ρ1β1(1 − 2β 2 1p 2 ) ω<br />
B1<br />
β1<br />
= ρ2β 2 ω<br />
2p cos ϕt A2 + ρ2β2(1 − 2β 2 2p 2 ) ω<br />
B2. (2.21)<br />
β2<br />
The four equations (2.8) ,(2.10), (2.17), and (2.21) are <strong>the</strong> final result after imposing <strong>the</strong><br />
boundary conditions. But <strong>the</strong>y are expressed in terms <strong>of</strong> <strong>the</strong> potential amplitudes. As <strong>the</strong><br />
reflection and transmission coefficients are defined as ratios <strong>of</strong> <strong>the</strong> displacement amplitude,<br />
we can rewrite <strong>the</strong> four equations as<br />
where<br />
−αpRpp − cos ϕrRps + α2pTpp + cos ϕtTps = αp,<br />
cos θiRpp − β1pRps + cos θtTpp − β2pTps = cos θi,<br />
−ρ1α1(1 − 2β 2 1p 2 )Rpp − 2ρ1β 2 1p cos ϕrRps + ρ2α2(1 − 2β 2 2p 2 )Tpp<br />
−2ρ2β 2 2p cos ϕtTpp = ρ1α1(1 − 2β 2 1p 2 ),<br />
2ρ1β 2 1p cos θiRpp + ρ1β1(1 − 2β 2 1p 2 )Rps + ρ2β 2 2p cos ϕtTpp<br />
Rpp = A1<br />
A0<br />
+ρ2β2(1 − 2β 2 2p 2 )Tps = 2ρ1β 2 1p cos θi, (2.22)<br />
, Rps = α1<br />
B1<br />
β1 A0<br />
, Tpp = α1<br />
A2<br />
α2 A0<br />
, Tps = α1<br />
B2<br />
β2 A0<br />
which are <strong>the</strong> reflection and transmission coefficients. The elegant way <strong>of</strong> expressing <strong>the</strong><br />
Zoeppritz equation is in matrix form. In matrix form, equation (2.22) takes <strong>the</strong> form<br />
MX = Y, (2.23)<br />
where<br />
⎛<br />
⎜<br />
M = ⎜<br />
⎝<br />
−αp<br />
cos θi<br />
−ρ1α1(1 − 2β<br />
− cos ϕr<br />
−β1p<br />
α2p<br />
cos θt<br />
cos ϕt<br />
−β2p<br />
2 1p2 ) −2ρ1β2 1p cos ϕr ρ2α2(1 − 2β2 2p2 ) −2ρ2β2 2p cos ϕt<br />
2ρ1β2 1p cos θi ρ1β1(1 − 2β2 1p2 ) ρ2β2 2p cos ϕt 2ρ2β2 ⎞<br />
⎟<br />
⎠<br />
2p cos ϕt<br />
,<br />
⎛<br />
⎜<br />
X = ⎜<br />
⎝<br />
Rpp<br />
Rps<br />
Tpp<br />
Tps<br />
⎞<br />
⎟<br />
⎠<br />
⎛<br />
⎜<br />
, Y = ⎜<br />
⎝<br />
αp<br />
cos θi<br />
ρ1α1(1 − 2β 2 1p 2 )<br />
2ρ1β 2 1p cos θi<br />
⎞<br />
,<br />
α2<br />
⎟<br />
⎠ .