# User:Ageary/Chapter 13

Problems-in-Exploration-Seismology-and-their-Solutions.jpg | |

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

Title | Problems in Exploration Seismology and their Solutions |

Author | Lloyd P. Geldart and Robert E. Sheriff |

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

ISBN | ISBN 9781560801153 |

Store | SEG Online Store |

## Contents

- 1 Chapter 13 Specialized techniques
- 2 13.1 S-wave conversion in marine surveys
- 3 13.2 Equally inclined orthogonal geophones
- 4 13.3 Guided (Channel) waves; Normal-mode propagation
- 5 13.4 Vertical seismic profiling
- 6 13.5 Effect of velocity change on VSP traveltime
- 7 13.6 Mapping the vertical flank of a salt dome
- 8 13.7 Poission’s ratio from P- and S-wave traveltimes

## Chapter 13 Specialized techniques

## 13.1 S-wave conversion in marine surveys

**13.1a In a marine survey, the water depth is 100 m and a reflector is 3 km below the seafloor. Use Figure 13.1a to determine the optimum range of offsets for S-wave generation. Take Poisson’s ratio just below the seafloor as 0.35 and the P-wave velocity as 2.8 km/s. The velocity in sea water is 1.5 km/s.**

**Solution**

Referring to Figure 13.1a(i), we see that the offset for the 3-km reflector is

**(**)

where is the angle of incidence of the P-wave at the bottom of the water layer, and is the angle of refraction of the converted S-wave.

To find for a given , we need the S-wave velocity . Using equation (1,8) in Table 2.2a, we find that for , so . To get the angles of incidence and of refraction , we have from equation (3.1a),

The optimum angles of incidence for S-wave conversion for are obtained by interpolating between the curves for and 3.0 in Figure 13.1a(ii). This gives a range of about to for . The corresponding values of are and . Equation (13.1a) now gives offsets of 4 and 11 km for and .

**13.1b Most marine S-wave surveys that wish to deal with S-waves utilize conversion at the reflector and recording with three-component geophones laid on the seafloor (OBC, ocean-bottom cables), so that only one mode conversion is involved. Assume a ray leaving an airgun source at to the vertical in water 100 m deep, an increase in P-wave velocities from 1.5 km/s at the seafloor to 3.0 km/s at a reflector 3 km below the sea floor (average velocity in the sediments of 2.25 km/s). Take in the sediments as 0.3. What will be the source-geophone offset?**

**Solution**

The P-wave gives an incident angle at the reflector of . Using equation 1,8 in Table 2.2a, the S-wave velocity here is about 0.53 or 1.59 km/s. Hence the angle of reflection is . The average direction of the P-wave ray in the sediments is . With constant in the sediments, is constant and, hence, . Hence, the average direction of the S-wave is . The offset is thus

This is a very reasonable offset.

## 13.2 Equally inclined orthogonal geophones

**Determine the inclination angles for three orthogonal geophones that are equally inclined to the vertical.**

**Solution**

The three orthogonal geophones define a rectangular coordinate system with the -, -, and -axes along the geophones axes. We take a straight line through the origin of this system such that it is equally inclined to each axis, and rotate the coordinate system so that the straight line is vertical. The direction cosines of the vertical line are and . Because [Sheriff and Geldart, problem (15.9a)],

Thus the geophones must be inclined to the vertical.

## 13.3 Guided (Channel) waves; Normal-mode propagation

**13.3a In Figure 13.3a the first arrival () has traveled at the velocity 2.7 km/s; find the water depth.**

**Background**

A wave guide is a layer in which a wave can propagate with little loss of energy. In a water layer nearly total reflection occurs at both boundaries, at the surface because of the very large impedance contrast and at the bottom reflection beyond the critical angle. The phase is inverted at the surface because the reflectivity is nearly , but not at the seafloor (beyond the critical angle) until the angle becomes very large.

Figure 13.3b(i) shows waves bouncing back and forth in a wave guide. For certain incident angles and frequencies , they interfere constructively. In Figure 13.3b(ii), is a wavefront traveling upward at the angle . The previous cycle of a parallel wavefront that passed through the position earlier followed paths such as *EFGH* and *BDA* and thus coincides with . Clearly . , , so . Taking into account the phase reversal at the surface the condition for constructive interference is

**(**)

Writing for the phase velocity, the frequencies that are reinforced are

**(**)

In addition to the upgoing waves parallel to and in Figure 13.3b(iii), a second downgoing set and will combine with the set to build up the energy traveling in the direction along the wave guide. Energy travels from to in the time that wavefront moves to , so the phase velocity of the energy traveling along , , is

**(**)

Since both and are functions of , is dispersive with a group velocity given by equation (2.7a):

**(**)

The derivative is always negative for a water channel wave (Figure 13.3c), so .

If the wave-guide effect had been due to reflection beyond the critical angle at both boundaries, as with a low-velocity coal seam, the phase would not have been reversed at either boundary and equations (13.3a,b) would reduce to

**(**)

**Solution**

Equation (13.3c) gives

and we use equations (13.3b,e) with to get :

**13.3b What frequency is reinforced when ?**

**Solution**

We use equation (13.3b) with , , :

## 13.4 Vertical seismic profiling

**13.4a (a) A source is offset 1000 m from a vertical well in which a geophone is suspended, and a horizontal reflector is present at a depth of 2000 m. Calculate the traveltimes for the reflection when the geophone is 800, 1200, and 1600 m deep if the velocity =3000 m/s**

**Solution**

Using a coordinate system with origin at the source and the -axis positive downward, the image point (see problem 4.1) in Figure 13.4a is at (0, 4000) and the geophone coordinates are (1000, ); hence,

For 800 m depth ; for 1200 m, ; for 1600 m, .

**13.4b Repeat for reflectors dipping , where the reflector intersects the well at the same point as in part (a) (see Figure 13.4b).**

**Solution**

i) __Dip down toward the well__

ii) __Dip down toward source __

Coordinates of image are (, 4170),

Table 13.4a shows calculated traveltimes for , 1200, and 1600 m.

Depth(m) | 800 | 1200 | 1600 | ||

time (s) | 1.09 | 0.97 | 0.86 | ||

time (s) | 1.14 | 1.00 | 0.87 |

800 | 958 | 799 | 1.114 | |

1200 | 937 | 1198 | 0.985 | |

1600 | 916 | 1598 | 0.857 |

**13.4c By how much do the values in part (a) change if the well deviates by towards the source (see Figure 13.4c)?**

**Solution**

Well inclination changes the geophone coordiantes from (1000, ) where the in-hole depth , 1200, 1600, to where

where subscript denotes values for the deviated well. The image point remains at (0, 4000), so the traveltime is

**(**)

Substituting the values of , we get the results in Table 13.4b. is the difference in traveltime from that calculated in part (a).

## 13.5 Effect of velocity change on VSP traveltime

**13.5a A source is offset 3.00 km west of a vertical well in which a geophone is suspended. There is a vertical north-south fault 0.80 km west of the well with velocities and 3.00 km/s west and east of the fault, respectively. A horizontal reflector is present west of the fault at a depth of 2.00 km. Find the reflection traveltime for geophone depths of 0.60 and 1.20 km.**

**Solution**

This problem must be solved by trial-and-error methods, either graphical or numerical. We have chosen the latter, because it is easier, quicker, and more accurate than graphical methods.

Figure 13.5a shows the geometry. The horizontal reflector at a depth of 2 km has an image point 4.00 km below the source. A typical ray is shown going from the image point to a geophone located in the well at a depth .

The angle of incidence and of refraction where the ray passes through the fault are related by Snell’s law: that is, . We have a further constraint in that the vertical distance between the image and the geophone is . These relations fix :

- ;
- .

We denote the left-hand side of (b) by . For the geophone at depth 0.60 km, the bending of the ray will have a small effect because of the shallowness and proximity to the fault, so we take as our first approximation the straight path,

To allow for some bending, we try and get . Substituting theses values in relation (b) we get

This should be , so we must increase by a small amount; taking gives and . We now find the traveltime:

For , . We know that will be less than for the shallower geophone, so we take and get , . Taking , , . The change of increased by 0.18 m; we need an increase of 0.10, so we take , , and . The traveltime is

**13.5b What is the deepest geophone location for which the reflection can be recorded?**

**Solution**

The deepest geophone position is that for which the ray passes through the intersection of the fault and reflector. For this point,

## 13.6 Mapping the vertical flank of a salt dome

**13.6a Two media of velocities and are separated by a vertical plane interface. A source is located at on the surface of the high-velocity medium and a geophone at point in a well in the other medium. If is the traveltime of a wave from to along the path SAG where is on the salt/sediment interface, discuss the locus of points located at the intersection of arcs centered at and having radii and , respectively, where .**

**Solution**

Any path where the interface is tangent to the locus curve satisfies the observed traveltime. To draw the locus we try a series of points such that . If a number of loci can be drawn for various geophone (or source) loctions, then their common tangent must locate the interface.

While this problem illustrates the concept, it can also be applied in three dimensions to allow for situations where and do not lie in the same plane.

500 | 0.44 | 1250 | 0.52 | 2000 | 0.6 |

750 | 0.46 | 1500 | 0.56 | 2250 | 0.67 |

1000 | 0.49 | 1750 | 0.60 | 2500 | 0.70 |

**13.6b An outcropping salt dome has roughly vertical flanks. A source is located on the salt and a geophone is suspended in a vertical well in the sediments 1600 m from the source point. Determine the outline of the salt dome from the data in Table 13.6a. Take the velocities in the salt and adjacent sediments as 5.00 and 3.00 km/s.**

**Solution**

The usual method of resolving this problem is to prepare a vertical section through the well W and source point S, then draw a series of circles centered at S with radii equal to the distances traveled in salt for convenient time intervals such as 0.2 s, 0.4 s, etc., the circles being labeled with the time value. A second set of concentric circles using the sediment velocity is drawn on a transparency; the center of these circles is then placed over a geophone position, and intersections of the two sets of circles are marked wherever the time values of the two radii add up to give the traveltime to the geophone.

An alternative method is to dispense with the circles and use only the radii. The centers of the circles and the extremities of the radii are marked along the edges of two narrow strips of light cardboard (one for salt, one for sediments), the time values being marked as before (the marks and the centers should be on opposite sides of the two strips to permit accurate determinations of the intersections). Map pins can be used to attach the zero points of the two strips to the vertical section, one at the source point and the other at the geophone location in the well. Intersections are found as before and marked directly on the section. After curves have been plotted for each geophone position, the flank is outlined by the curve through the apices of the curves. The results are shown in Figure 13.6a.

## 13.7 Poission’s ratio from P- and S-wave traveltimes

**Find Poisson’s ratio for the five events in Figure 13.7a.**

**Solution**

Poisson’s ratio can be obtained from the ratio using equation (10.2) Table 2.2a. Since , , we can get from the traveltimes and then get using equation (10,2) (see Table 2.2a), that is,

Measurements and calculations are listed in Table 13.7a.

Event | |||||
---|---|---|---|---|---|

1 | 0.39 | 0.96 | 2.46 | 6.05 | 0.40 |

2 | 0.62 | 1.50 | 2.42 | 5.86 | 0.40 |

3 | 0.74 | 1.65 | 2.23 | 4.97 | 0.37 |

4 | 0.80 | 1.83 | 2.29 | 5.24 | 0.38 |

5 | 0.99 | 2.16 | 2.18 | 4.75 | 0.37 |