# Difference between revisions of "Water reverberations"

m |
|||

(4 intermediate revisions by the same user not shown) | |||

Line 1: | Line 1: | ||

+ | <languages/> | ||

+ | <translate> | ||

+ | </translate> | ||

{{publication | {{publication | ||

| name = Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing | | name = Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing | ||

Line 9: | Line 12: | ||

| frompg = 193 | | frompg = 193 | ||

| topg = 216 | | topg = 216 | ||

− | | author = [[Enders A. Robinson]] and [[Sven Treitel]] | + | | author = [[Special:MyLanguage/Enders A. Robinson|Enders A. Robinson]] and [[Special:MyLanguage/Sven Treitel|Sven Treitel]] |

| doi = 10.1190/1.9781560801610 | | doi = 10.1190/1.9781560801610 | ||

| pdf = | | pdf = | ||

Line 15: | Line 18: | ||

| isbn = 9781560801481 | | isbn = 9781560801481 | ||

}} | }} | ||

− | The water-reverberation problem in marine seismic operations can be described as follows: Consider just two interfaces — the water surface and the water bottom — with a separation given by the two-way traveltime ''T''. The reflectivity is <math>\{ \varepsilon _0 ,\;0,\;0,\;...,\;0,\;\varepsilon _T \} </math>. For definiteness, let us assume that ''T'' = 3. Then the reflectivity is <math>\{ \varepsilon _0 ,\;0,\;0,\;\varepsilon _3 \} </math>. The two-sided autocorrelation of the reflectivity is | + | <translate> |

+ | <!--T:1--> | ||

+ | The water-reverberation problem in marine seismic operations can be described as follows: Consider just two interfaces — the water surface and the water bottom — with a separation given by the two-way traveltime ''T''. The reflectivity is <math>\{ \varepsilon _0 ,\;0,\;0,\;...,\;0,\;\varepsilon _T \} </math>. For definiteness, let us assume that ''T'' = 3. Then the reflectivity is <math>\{ \varepsilon _0 ,\;0,\;0,\;\varepsilon _3 \} </math>. The two-sided autocorrelation ''g''<sub>''t''</sub> of the reflectivity is | ||

+ | </translate> | ||

+ | <translate><!--T:2--> | ||

{{NumBlk|:|<math>\begin{align} | {{NumBlk|:|<math>\begin{align} | ||

\left(g_{3}{ ,0,0,\ }g_0{ ,0,0,\ }g_{3}\right)=\left({\varepsilon }_0{\varepsilon }_{2}{ \ ,\ 0,0,\ }{\varepsilon }^{2}_0+{\varepsilon }^{2}_{2}{ ,0,0,\ }{\varepsilon }_0{\varepsilon }_{3}\right) . | \left(g_{3}{ ,0,0,\ }g_0{ ,0,0,\ }g_{3}\right)=\left({\varepsilon }_0{\varepsilon }_{2}{ \ ,\ 0,0,\ }{\varepsilon }^{2}_0+{\varepsilon }^{2}_{2}{ ,0,0,\ }{\varepsilon }_0{\varepsilon }_{3}\right) . | ||

\end{align}</math>|{{EquationRef|16}}}} | \end{align}</math>|{{EquationRef|16}}}} | ||

+ | <!--T:3--> | ||

First, take the zero-lag term and the right-hand side of this autocorrelation, and then replace the zero-lag term by unity. The result is | First, take the zero-lag term and the right-hand side of this autocorrelation, and then replace the zero-lag term by unity. The result is | ||

+ | </translate> | ||

+ | <translate><!--T:4--> | ||

{{NumBlk|:|<math>\begin{align} | {{NumBlk|:|<math>\begin{align} | ||

\left({ 1,0,0,\ }g_{3}\right)=\left({ 1,\ 0,0,\ }{\varepsilon }_0{\varepsilon }_{3}\right) . | \left({ 1,0,0,\ }g_{3}\right)=\left({ 1,\ 0,0,\ }{\varepsilon }_0{\varepsilon }_{3}\right) . | ||

\end{align}</math>|{{EquationRef|17}}}} | \end{align}</math>|{{EquationRef|17}}}} | ||

+ | <!--T:5--> | ||

As a ''Z''-transform, this result is | As a ''Z''-transform, this result is | ||

+ | </translate> | ||

+ | <translate><!--T:6--> | ||

{{NumBlk|:|<math>\begin{align} | {{NumBlk|:|<math>\begin{align} | ||

{ 1+}g_{3}Z^{3}={1+}{\varepsilon }_0{\varepsilon }_{3}Z^{3}. | { 1+}g_{3}Z^{3}={1+}{\varepsilon }_0{\varepsilon }_{3}Z^{3}. | ||

Line 34: | Line 47: | ||

|{{EquationRef|18}}}} | |{{EquationRef|18}}}} | ||

+ | <!--T:7--> | ||

We will show now that this ''Z''-transform represents the feedback loop that produces a reverberating wavetrain. | We will show now that this ''Z''-transform represents the feedback loop that produces a reverberating wavetrain. | ||

+ | <!--T:8--> | ||

The surface (the water-air interface) has the reflection coefficient <math>{\varepsilon }_0</math>. Because the surface is a strong reflector, this reflection coefficient is close to unity in magnitude. The reflection coefficient for a wave striking the water surface from below is <math>-{\varepsilon }_0</math>. | The surface (the water-air interface) has the reflection coefficient <math>{\varepsilon }_0</math>. Because the surface is a strong reflector, this reflection coefficient is close to unity in magnitude. The reflection coefficient for a wave striking the water surface from below is <math>-{\varepsilon }_0</math>. | ||

+ | <!--T:9--> | ||

The water-bottom interface also is a strong reflector, with the reflection coefficient <math>\varepsilon _3 </math>. Thus, the water layer acts as an imperfect energy trap in which a seismic pulse is reflected successively between its two interfaces. Energy that bounces back and forth between two interfaces is called a reverberation. | The water-bottom interface also is a strong reflector, with the reflection coefficient <math>\varepsilon _3 </math>. Thus, the water layer acts as an imperfect energy trap in which a seismic pulse is reflected successively between its two interfaces. Energy that bounces back and forth between two interfaces is called a reverberation. | ||

+ | <!--T:10--> | ||

Let us examine the dynamics of this water reverberation. Our approach is applicable to any reverberation occurring between two interfaces, not just to the water reverberation occurring between the water-surface layer and the water-bottom layer encountered in marine seismic work. | Let us examine the dynamics of this water reverberation. Our approach is applicable to any reverberation occurring between two interfaces, not just to the water reverberation occurring between the water-surface layer and the water-bottom layer encountered in marine seismic work. | ||

+ | <!--T:11--> | ||

In summary, the source and receiver locations are just below the water surface. Source energy travels down to the water bottom, from which it is reflected. The reflected upgoing energy proceeds toward the water surface, where it is recorded as a primary reflection. However, this upgoing energy also is reflected from the surface and then continues to bounce back and forth in the water layer. These multiple reflections within the water layer make up the water-layer reverberation. Such reverberations are generally undesirable because they obscure reflections from deeper horizons. Figure 8 shows the raypaths for the water reverberation. (For clarity, the raypaths have been drawn as slanting lines, although in our normal-incidence model, they are perpendicular to the two interfaces.) | In summary, the source and receiver locations are just below the water surface. Source energy travels down to the water bottom, from which it is reflected. The reflected upgoing energy proceeds toward the water surface, where it is recorded as a primary reflection. However, this upgoing energy also is reflected from the surface and then continues to bounce back and forth in the water layer. These multiple reflections within the water layer make up the water-layer reverberation. Such reverberations are generally undesirable because they obscure reflections from deeper horizons. Figure 8 shows the raypaths for the water reverberation. (For clarity, the raypaths have been drawn as slanting lines, although in our normal-incidence model, they are perpendicular to the two interfaces.) | ||

+ | <!--T:12--> | ||

[[file:Ch08_fig8.png|thumb|{{figure number|8.}} The downward pass through the water layer.]] | [[file:Ch08_fig8.png|thumb|{{figure number|8.}} The downward pass through the water layer.]] | ||

+ | <!--T:13--> | ||

In our example, the integer ''T'' = 3 represents the two-way traveltime parameter in the water layer. The source is a unit impulse at ''A'', which is a point just below the surface. The pulse travels downward in the water and arrives at the water bottom. There, it is reflected. The resulting upgoing pulse at ''B'' has value <math>{\varepsilon }_{3}</math> and arrives at the water surface ''C'' at time 3. The receiver records this upgoing pulse as the primary reflection. However, this upgoing pulse then is reflected downward. In the case of an upcoming incident wave, the water surface has a reflection coefficient of –<math>{\varepsilon }_{3}</math>, so the downgoing pulse at ''D'' has the value <math>-{\varepsilon }_0{\varepsilon }_{3}</math>. This downgoing pulse travels downward and is reflected from the water bottom at ''E''. It returns to the surface at ''F''. The receiver records this upgoing pulse as the first multiple reflection, with the value <math>{\varepsilon }_{3}\left(-{\varepsilon }_0{\varepsilon }_{3}\right)</math> occurring at time 2''T'' = 6. This repeats to produce the second multiple reflection with the value <math>{\varepsilon }_{3}{\left(-{\varepsilon }_0{\varepsilon }_{3}\right)}^{2}</math> occurring at time 3''T'' = 9. The process keeps repeating itself, each time producing another multiple reflection. Because we omit the shot, we write the trace as | In our example, the integer ''T'' = 3 represents the two-way traveltime parameter in the water layer. The source is a unit impulse at ''A'', which is a point just below the surface. The pulse travels downward in the water and arrives at the water bottom. There, it is reflected. The resulting upgoing pulse at ''B'' has value <math>{\varepsilon }_{3}</math> and arrives at the water surface ''C'' at time 3. The receiver records this upgoing pulse as the primary reflection. However, this upgoing pulse then is reflected downward. In the case of an upcoming incident wave, the water surface has a reflection coefficient of –<math>{\varepsilon }_{3}</math>, so the downgoing pulse at ''D'' has the value <math>-{\varepsilon }_0{\varepsilon }_{3}</math>. This downgoing pulse travels downward and is reflected from the water bottom at ''E''. It returns to the surface at ''F''. The receiver records this upgoing pulse as the first multiple reflection, with the value <math>{\varepsilon }_{3}\left(-{\varepsilon }_0{\varepsilon }_{3}\right)</math> occurring at time 2''T'' = 6. This repeats to produce the second multiple reflection with the value <math>{\varepsilon }_{3}{\left(-{\varepsilon }_0{\varepsilon }_{3}\right)}^{2}</math> occurring at time 3''T'' = 9. The process keeps repeating itself, each time producing another multiple reflection. Because we omit the shot, we write the trace as | ||

+ | </translate> | ||

+ | <translate><!--T:14--> | ||

{{NumBlk|:|<math>\begin{align} | {{NumBlk|:|<math>\begin{align} | ||

− | h={\ }\left({ 0,0,0,\ }{\varepsilon }_{3}{ ,0,0,\ }{\varepsilon }_{3}\left(-{\varepsilon }_0{\varepsilon }_{3}\right){ ,\ 0,0,\ }{\varepsilon }_{1}{\left(-{\varepsilon }_0{\varepsilon }_{3}\right)}^{2}{ \ ,\ 0,0,\ }\dots \right) , | + | h={\ }\left({ 0,0,0,\ }{\varepsilon }_{3}{ ,0,0,\ }{\varepsilon }_{3}\left(-{\varepsilon }_0{\varepsilon }_{3}\right){ ,\ 0,0,\ }{\varepsilon }_{1}{\left(-{\varepsilon }_0{\varepsilon }_{3}\right)}^{2}{ \ ,\ 0,0,,\ }\dots \right) , |

\end{align}</math>|{{EquationRef|19}}}} | \end{align}</math>|{{EquationRef|19}}}} | ||

+ | <!--T:15--> | ||

where the first 0 is at time 0, the second 0 is at time 1, the third 0 is at time 2, the <math>\varepsilon _3 </math> is at time 3, and so on. In other words, the shot occurs at time 0, the first nonzero value <math>{\varepsilon }_{3}</math> (the primary) occurs at time ''T'' = 3, the second nonzero value (the first multiple) occurs at time 6, the third nonzero value (the second multiple) occurs at time 9, and so on. The successive nonzero values (after the primary) are separated by the time parameter ''T'' = 3. For that reason, ''T'' is called the cycle time of the reverberation. Thus, the ''Z''-transform of the impulsive trace is | where the first 0 is at time 0, the second 0 is at time 1, the third 0 is at time 2, the <math>\varepsilon _3 </math> is at time 3, and so on. In other words, the shot occurs at time 0, the first nonzero value <math>{\varepsilon }_{3}</math> (the primary) occurs at time ''T'' = 3, the second nonzero value (the first multiple) occurs at time 6, the third nonzero value (the second multiple) occurs at time 9, and so on. The successive nonzero values (after the primary) are separated by the time parameter ''T'' = 3. For that reason, ''T'' is called the cycle time of the reverberation. Thus, the ''Z''-transform of the impulsive trace is | ||

+ | <!--T:16--> | ||

<math>H\left(Z\right)={\varepsilon }_{3}Z^{3}+{\varepsilon }_{3}\left(-{\varepsilon }_0{\varepsilon }_{3}\right)Z^{6}+{\varepsilon }_{3}{\left(-{\varepsilon }_0{\varepsilon }_{3}\right)}^{2}Z^{9}+\dots</math> | <math>H\left(Z\right)={\varepsilon }_{3}Z^{3}+{\varepsilon }_{3}\left(-{\varepsilon }_0{\varepsilon }_{3}\right)Z^{6}+{\varepsilon }_{3}{\left(-{\varepsilon }_0{\varepsilon }_{3}\right)}^{2}Z^{9}+\dots</math> | ||

+ | </translate> | ||

+ | <translate><!--T:17--> | ||

{{NumBlk|:|<math>\begin{align} | {{NumBlk|:|<math>\begin{align} | ||

− | ={\varepsilon }_{3}Z^{3}\left[{1}+\left(-{\varepsilon }_0{\varepsilon }_{3}\right)Z^{3}+{\left(-{\varepsilon }_0{\varepsilon }_{3}\right)}^{2}Z^{6}+\dots \right]=\frac{{\varepsilon }_{3}Z^{3}}{{1}+{\varepsilon }_0{\varepsilon }_{3}Z^{3}}. | + | ={\varepsilon }_{3}Z^{3}\left[{1}+\left(-{\varepsilon }_0{\varepsilon }_{3}\right)Z^{3}+{\left(-{\varepsilon }_0{\varepsilon }_{3}\right)}^{2}Z^{6}+\dots \right]=\frac{{\varepsilon }_{3}Z^{3}} |

+ | {{1}+{\varepsilon }_0{\varepsilon }_{3}Z^{3}}. | ||

\end{align}</math> | \end{align}</math> | ||

|{{EquationRef|20}}}} | |{{EquationRef|20}}}} | ||

+ | <!--T:18--> | ||

The geometric series within square brackets is convergent and summable because the autocorrelation coefficient <math>g_3 = \varepsilon _0 \;\varepsilon _3</math> is less than one in magnitude. Now comes an important point. We recognize that the multiple-generating factor <math>{\varepsilon }_0{\varepsilon }_{3}</math> that occurs in the above expression for the trace ''X'' is the autocorrelation coefficient <math>g_{3}</math>. Thus, the ''Z''-transform becomes | The geometric series within square brackets is convergent and summable because the autocorrelation coefficient <math>g_3 = \varepsilon _0 \;\varepsilon _3</math> is less than one in magnitude. Now comes an important point. We recognize that the multiple-generating factor <math>{\varepsilon }_0{\varepsilon }_{3}</math> that occurs in the above expression for the trace ''X'' is the autocorrelation coefficient <math>g_{3}</math>. Thus, the ''Z''-transform becomes | ||

+ | </translate> | ||

+ | <translate><!--T:19--> | ||

{{NumBlk|:|<math>\begin{align} | {{NumBlk|:|<math>\begin{align} | ||

− | H\left(Z\right)=\frac{{\varepsilon }_{3}Z^{3}}{{1}+g_{3}Z^{3}}. | + | H\left(Z\right)=\frac{{\varepsilon }_{3}Z^{3}} |

+ | {{1}+g_{3}Z^{3}}. | ||

\end{align}</math>|{{EquationRef|21}}}} | \end{align}</math>|{{EquationRef|21}}}} | ||

+ | <!--T:20--> | ||

We therefore have shown that the primary reflection is represented by the feedforward term <math>{\varepsilon }_{3}Z^{3}</math>, whereas the reverberation is represented by the feedback term <math>1+g_{3}Z^{3}</math> (Figure 9). | We therefore have shown that the primary reflection is represented by the feedforward term <math>{\varepsilon }_{3}Z^{3}</math>, whereas the reverberation is represented by the feedback term <math>1+g_{3}Z^{3}</math> (Figure 9). | ||

[[file:Ch08_fig9.png|thumb|{{figure number|9.}} Block diagram of the feedback system that generates reverberation.]] | [[file:Ch08_fig9.png|thumb|{{figure number|9.}} Block diagram of the feedback system that generates reverberation.]] | ||

__NOTOC__ | __NOTOC__ | ||

− | == Continue reading == | + | |

+ | == Continue reading == <!--T:21--> | ||

+ | |||

+ | <!--T:22--> | ||

{| border="1" cellpadding="5" cellspacing="0" align="center" | {| border="1" cellpadding="5" cellspacing="0" align="center" | ||

|- | |- | ||

Line 76: | Line 111: | ||

! style="background: #E5D75F;" | Next section | ! style="background: #E5D75F;" | Next section | ||

|- | |- | ||

− | | align="center" | [[Synthetic seismogram without multiples]] | + | | align="center" | [[Special:MyLanguage/Synthetic seismogram without multiples|Synthetic seismogram without multiples]] |

− | | align="center" | [[Synthetic seismogram with multiples]] | + | | align="center" | [[Special:MyLanguage/Synthetic seismogram with multiples|Synthetic seismogram with multiples]] |

|- | |- | ||

! style="background: #426580; color: white;" | Previous chapter | ! style="background: #426580; color: white;" | Previous chapter | ||

Line 89: | Line 124: | ||

|} | |} | ||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | == External links == | + | == Also in this chapter == <!--T:23--> |

− | {{search}} | + | |

+ | |||

+ | <!--T:24--> | ||

+ | * [[Special:MyLanguage/Introduction - Chapter 8|Introduction]] | ||

+ | * [[Special:MyLanguage/Polarity|Polarity]] | ||

+ | * [[Special:MyLanguage/Reflection coefficients and transmission coefficients|Reflection coefficients and transmission coefficients]] | ||

+ | * [[Special:MyLanguage/Ghost reflection|Ghost reflection]] | ||

+ | * [[Special:MyLanguage/Layer-cake model|Layer-cake model]] | ||

+ | * [[Special:MyLanguage/Synthetic seismogram without multiples|Synthetic seismogram without multiples]] | ||

+ | * [[Special:MyLanguage/Synthetic seismogram with multiples|Synthetic seismogram with multiples]] | ||

+ | * [[Special:MyLanguage/Examples|Examples]] | ||

+ | * [[Special:MyLanguage/Small and white reflection coefficients|Small and white reflection coefficients]] | ||

+ | * [[Appendix H: Exercises|Appendix H: Exercises]] | ||

+ | |||

+ | |||

+ | == External links == <!--T:25--> | ||

+ | |||

+ | </translate> | ||

+ | {{search}} | ||

+ | <translate> | ||

+ | <!--T:26--> | ||

[[Category:Digital Imaging and Deconvolution]] | [[Category:Digital Imaging and Deconvolution]] | ||

[[Category:Publications]] | [[Category:Publications]] | ||

[[Category:Synthetics]] | [[Category:Synthetics]] | ||

+ | </translate> |

## Latest revision as of 09:58, 4 May 2021

Series | Geophysical References Series |
---|---|

Title | Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing |

Author | Enders A. Robinson and Sven Treitel |

Chapter | 8 |

DOI | http://dx.doi.org/10.1190/1.9781560801610 |

ISBN | 9781560801481 |

Store | SEG Online Store |

The water-reverberation problem in marine seismic operations can be described as follows: Consider just two interfaces — the water surface and the water bottom — with a separation given by the two-way traveltime *T*. The reflectivity is . For definiteness, let us assume that *T* = 3. Then the reflectivity is . The two-sided autocorrelation *g*_{t} of the reflectivity is

**(**)

First, take the zero-lag term and the right-hand side of this autocorrelation, and then replace the zero-lag term by unity. The result is

**(**)

As a *Z*-transform, this result is

**(**)

We will show now that this *Z*-transform represents the feedback loop that produces a reverberating wavetrain.

The surface (the water-air interface) has the reflection coefficient . Because the surface is a strong reflector, this reflection coefficient is close to unity in magnitude. The reflection coefficient for a wave striking the water surface from below is .

The water-bottom interface also is a strong reflector, with the reflection coefficient . Thus, the water layer acts as an imperfect energy trap in which a seismic pulse is reflected successively between its two interfaces. Energy that bounces back and forth between two interfaces is called a reverberation.

Let us examine the dynamics of this water reverberation. Our approach is applicable to any reverberation occurring between two interfaces, not just to the water reverberation occurring between the water-surface layer and the water-bottom layer encountered in marine seismic work.

In summary, the source and receiver locations are just below the water surface. Source energy travels down to the water bottom, from which it is reflected. The reflected upgoing energy proceeds toward the water surface, where it is recorded as a primary reflection. However, this upgoing energy also is reflected from the surface and then continues to bounce back and forth in the water layer. These multiple reflections within the water layer make up the water-layer reverberation. Such reverberations are generally undesirable because they obscure reflections from deeper horizons. Figure 8 shows the raypaths for the water reverberation. (For clarity, the raypaths have been drawn as slanting lines, although in our normal-incidence model, they are perpendicular to the two interfaces.)

In our example, the integer *T* = 3 represents the two-way traveltime parameter in the water layer. The source is a unit impulse at *A*, which is a point just below the surface. The pulse travels downward in the water and arrives at the water bottom. There, it is reflected. The resulting upgoing pulse at *B* has value and arrives at the water surface *C* at time 3. The receiver records this upgoing pulse as the primary reflection. However, this upgoing pulse then is reflected downward. In the case of an upcoming incident wave, the water surface has a reflection coefficient of –, so the downgoing pulse at *D* has the value . This downgoing pulse travels downward and is reflected from the water bottom at *E*. It returns to the surface at *F*. The receiver records this upgoing pulse as the first multiple reflection, with the value occurring at time 2*T* = 6. This repeats to produce the second multiple reflection with the value occurring at time 3*T* = 9. The process keeps repeating itself, each time producing another multiple reflection. Because we omit the shot, we write the trace as

**(**)

where the first 0 is at time 0, the second 0 is at time 1, the third 0 is at time 2, the is at time 3, and so on. In other words, the shot occurs at time 0, the first nonzero value (the primary) occurs at time *T* = 3, the second nonzero value (the first multiple) occurs at time 6, the third nonzero value (the second multiple) occurs at time 9, and so on. The successive nonzero values (after the primary) are separated by the time parameter *T* = 3. For that reason, *T* is called the cycle time of the reverberation. Thus, the *Z*-transform of the impulsive trace is

**(**)

The geometric series within square brackets is convergent and summable because the autocorrelation coefficient is less than one in magnitude. Now comes an important point. We recognize that the multiple-generating factor that occurs in the above expression for the trace *X* is the autocorrelation coefficient . Thus, the *Z*-transform becomes

**(**)

We therefore have shown that the primary reflection is represented by the feedforward term , whereas the reverberation is represented by the feedback term (Figure 9).

## Continue reading

Previous section | Next section |
---|---|

Synthetic seismogram without multiples | Synthetic seismogram with multiples |

Previous chapter | Next chapter |

Wavelets | Wavelet Processing |

## Also in this chapter

- Introduction
- Polarity
- Reflection coefficients and transmission coefficients
- Ghost reflection
- Layer-cake model
- Synthetic seismogram without multiples
- Synthetic seismogram with multiples
- Examples
- Small and white reflection coefficients
- Appendix H: Exercises