Method for processing borehole seismic data

Abstract

Seismic data recorded by subsurface seismic sensors placed in a borehole, such as an oil or gas well, are transformed via a process of upward wavefield propagation to pseudo-receivers at the surface of the earth. The seismic data thus transformed can be processed as though the data had been recorded by the pseudo-receivers at the surface rather than in the borehole where the data were actually recorded. This method accurately accounts for seismic source statics, anisotropy, and all velocity effects between the real receivers in the borehole and the pseudo-receivers at the surface of the earth.

Claims

1. A method, comprising: placing a plurality of receivers in at least one borehole; wherein receivers of the plurality of receivers are vertically spaced from each other; defining, by a processor, a 3D Vertical Seismic Profile (3D VSP) source point map for a 3D VSP survey of an area of an earth surface; wherein the 3D VSP source point map comprises a plurality of seismic source points; initiating a first seismic source corresponding to a first seismic source point; initiating at least one second seismic source corresponding to a plurality of second seismic source points; wherein the first seismic source is located at a first location that is at or near the earth surface of the area and the at least one second seismic source is located at a second location that is at or near the earth surface of the area; receiving, by the processor, first seismic data caused by the first seismic source and recorded by the plurality of receivers; wherein the first seismic data comprises actual seismograms representative of a first wavefield associated with the first seismic source; creating, by the processor, each respective pseudo-receiver of a plurality of pseudo-receivers at each respective location in the 3D VSP source point map; wherein each respective pseudo-receiver corresponds to each respective receiver of the plurality of receivers; wherein each respective location of each respective pseudo-receiver is related to a direction of the first wavefield; wherein each respective location of each respective pseudo-receiver is adjacent to at least one second seismic source point of the plurality of second seismic source points; determining, by the processor, each respective seismic travel time between each respective pseudo-receiver and each corresponding receiver of the plurality of receivers; and continuing, by the processor, the first wavefield of the first seismic source to the plurality of pseudo-receivers, by generating, for each respective pseudo-receiver, based on each respective seismic travel time, each respective synthetic seismogram from each respective actual seismogram recorded by each respective receiver of the plurality of receivers to form a plurality of synthetic seismograms of the first wavefield of the first seismic source that is utilized as a plurality of as-if surface-recorded seismograms.

2. The method of claim 1, wherein the determining each respective seismic travel time further comprises: interpolating measured first arrival times at each respective receiver from at least two second seismic source points adjacent to each respective pseudo-receiver.

3. The method of claim 1, wherein the determining each respective seismic travel time is based on a finite difference method.

4. The method of claim 1, wherein the determining each respective seismic travel time is based on ray tracing.

5. The method of claim 1, wherein each of the plurality of receivers is a multi-component geophone.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The present invention can be more readily understood in conjunction with the accompanying drawings, in which:

(2) FIG. 1 is a cross-sectional diagram of a portion of the earth's surface showing a number of source-receiver pairs in a common source gather.

(3) FIG. 2 is a cross-sectional diagram similar to FIG. 1 showing a number of source-receiver pairs in a common mid-point gather.

(4) FIG. 3 is a graph showing a synthetic common seismic depth point (CDP) gather with reflections that have not been corrected for normal moveout (NMO).

(5) FIG. 4 is a graph showing the synthetic common depth point (CDP) gather in FIG. 3 after correction for NMO.

(6) FIG. 5 is a cross-sectional diagram showing a number of source-receiver pairs in a common source gather recorded with seismometers in a borehole.

(7) FIG. 6 is a graph showing a synthetic VSP common source point gather.

(8) FIG. 7 is a cross-sectional diagram similar to FIG. 5 showing upward continuation of the wavefield to the surface of the earth.

(9) FIG. 8 is a cross-sectional diagram similar to FIG. 7 showing upward continuation of the wavefield to the surface of the earth by way of finite difference methods.

(10) FIG. 9 is a flow chart of the process for upward continuing borehole seismic data by measured, interpolated or computed travel time delays.

(11) FIG. 10 is a surface map showing seismic source point locations for 3D VSP as circles. The well location is depicted as a diamond.

(12) FIG. 11 is a cross-sectional diagram similar to FIG. 5 showing wavefield continuation via Huygens' principle using known or computed travel times.

(13) FIG. 12 is a surface map showing pseudo-receiver positions, p, for a source point, S.

(14) FIG. 13 is a graph of a synthetic common source gather generated by finite difference methods.

(15) FIG. 14 is a graph of an upward continued common source gather done using Huygens' principle and Equation 1.

(16) FIG. 15 is a graph of the CDP stack created by processing upward continued data as surface seismic data.

(17) FIG. 16 is a flow diagram of the process of upward continuing borehole seismic data by finite difference methods.

(18) FIG. 17 is a graph of a synthetic VSP common source point gather generated by finite difference methods.

(19) FIG. 18 is a graph of the upward continued common source gather shown in FIG. 17 after upward continuation by finite difference methods.

(20) FIG. 19 is a graph of the CDP stack created by processing upward continued data as surface seismic data. The data were upward continued by finite difference methods.

DETAILED DESCRIPTION OF THE INVENTION

(21) The present method begins with seismic data recorded as in FIG. 5 in which seismic receivers (seismometers, hydrophones, geophones, etc) are placed in a borehole such as an oil or gas well and seismic sources are initiated at or near the surface of the earth. Seismic waves propagate from the seismic source into the earth. Seismic energy propagates to the borehole receivers whereupon the seismic receivers generate a signal in response to the seismic energy. The seismometer signal is then typically transmitted to a recording instrument at the surface of the earth and recorded on an electronic medium such as magnetic tape, optical device, or hard disk.

(22) Some of the propagating energy travels directly from the seismic source to the receivers and provides a measure of the travel time from the seismic source position to the receiver position. Other energy reflects from subsurface interfaces and is also recorded at the receivers in the borehole. Other energy propagates and is recorded as converted waves and refractions and other energy propagates away from the borehole and is not recorded.

(23) After the data traces are recorded on electronic media, the data are processed to obtain information about the earth. FIG. 6 shows a common source point gather from a synthetic offset VSP dataset. Important components of the wavefield in offset VSP data are shown in this synthetic offset VSP gather. The first arrivals are the first energy that arrives at the receiver directly from the source. The time of the first arrival is a key element of the method described herein. The other important component is the reflections that are recorded at geophones in the borehole. The reflections are propagated to the surface of the earth via upward continuation or application of Huygens' principle and known travel times from the surface of the earth to receivers in the subsurface.

(24) A key step in the present process is to continue the VSP wavefield upward to the surface of the earth in order to represent the data as it would have been recorded at the surface of the earth rather than in the borehole. FIG. 7 depicts upward continuation of the wavefield to the surface of the earth. The dataset recorded in the borehole is thus transformed into a dataset that mimics the case in which the data had been recorded at the surface of the earth. The dashed arrows depict the wavefield being continued to the surface of the earth. After transformation of the data into surface seismic data, the data can be processed using surface seismic processing methods, including common mid-point sorting and the NMO equation. FIG. 7 shows pseudo-receivers P.sub.1 through P.sub.n at the surface of the earth. The term pseudo-receiver is used herein to indicate the location of a seismometer to which data was propagated from a real receiver position in a borehole.

(25) Upward continuation of the wavefield by Huygens' principle can also be modeled via finite difference or finite element methods as shown in FIG. 8. Finite difference propagation essentially propagates each amplitude recorded at each receiver through a velocity field to pseudo-receivers at the surface of earth. The data can then be processed as surface seismic data. Finite difference propagation of the wavefield is currently prohibitively expensive for common practice, though it does work and should be considered to be one possible embodiment of the present invention.

(26) To use the finite difference method the data processor must provide a velocity model of the earth through which to propagate the seismic amplitudes. The extent to which the upward continued data can be used for successfully obtaining an image is related to the accuracy of the velocity model. While finite difference methods can be used for upward continuation of the wavefield, the following alternative method is preferred for the following reasons: (1) the alternative method does not require the provision of a velocity model, but incorporates all anisotropic and statics factors in the upward continuation; and (2) the alternative method requires many orders of magnitude less computer time to compute the upward continued seismic data than finite difference or finite element methods.

(27) FIG. 9 is a flow diagram of the alternative, preferred process for upward continuing borehole seismic data by measured, interpolated, or computed travel time delays. This method of upward continuing the data to the surface of the earth begins with picking the first arrival times for source points in the survey. FIG. 10 shows a 3D VSP source point map for a hypothetical 3D VSP survey. Seismic source locations at the surface of the earth are depicted as filled circles. The well location at the surface of the earth is shown as a diamond near the center of the source points.

(28) Borehole seismic data is initially acquired as previously described (step 61 in FIG. 9) and shown for example in FIG. 5. The travel time from each source point location to every receiver in the borehole can be determined from the first arrival times on the shot records (step 62 in FIG. 9). The travel time from a point that is not on a real source point can be estimated by interpolating travel times based on nearby source points. Such a dataset would constitute a 3D VSP. The first arrival times are picked for as many source points as is required to accurately interpolate the first arrival time for any point at the surface of the earth within the survey area where source points are present. In normal surveys, all first arrivals are picked for each source point. The first arrival time information is stored. The first arrival information for each source-receiver pair must be stored in such a way that the source coordinates (X, Y, and Elevation) can be associated with the first arrival time and the specific receiver of the source-receiver pair. For example the data can be stored in columns as follows.

(29) TABLE-US-00001 SRC_X SRC_Y SRC_EL RNUM FB_PICK 100. 200. 344. 1 221.4 100. 200. 344. 2 225.6 100. 200. 344. 3 229.9 200. 350. 317. 1 188.0 200. 350. 317. 2 194.1 200. 350. 317. 3 199.4
Where SRC_X is the source point X-coordinate, SRC_Y is the source point Y-coordinate, SRC_EL is the source point elevation. RNUM is the receiver number or some number that uniquely identifies each receiver in the borehole, and FB_PICK is the first break pick time. In this example FB_PICK would be stored in milliseconds.

(30) The next step in the preferred method of upward continuation is to actually upward continue the data to the surface of the earth using Huygens' principle. As was discussed above, Huygens' principle states that a wavefront at one time and place is the summation of a set of virtual source points when the wavefront was at another location at another time. Berryhill, Wave-Equation Datuming, Geophysics, vol. 44, no. 8, pp. 1329-1344 (1979) noted that continuation of a wavefield from one location to another through time could be achieved if travel times between the wavefront and the point to which the wave was being continued could be obtained. Berryhill's function for continuation is:
S(t)=Wi(tTi)*f(t)(Equation 1)

(31) Where S(t) is the seismogram as a function of time, t, at the location to which the wavefield was propagated, Wi is the seismogram at the i.sup.th location from which the wavefield was propagated, Ti is the seismic travel time for a wave between the position S and the position Wi, and f(t) is a filter that, while recommended by Berryhill, is optional. FIG. 11 is a cross-sectional diagram illustrating wavefield continuation via Huygens' principle using known or computed travel times. The wavefield can be continued to pseudo-receivers at the surface of the earth in a much faster way than finite difference upward continuation. Berryhill discusses a method of fast wavefield continuation by application of Huygens' principle if the seismic travel time between a real receiver and a pseudo-receiver are known. In the case of 3D VSP data or offset VSP data, the travel times from the surface of the earth to the receivers in the borehole are either known from direct first break pick times (as depicted by the dashed lines) or can be computed if the observed travel times are not known. Thus the wavefield as it would have been recorded at the pseudo-receivers can be quickly and accurately approximated by summing wavefield amplitudes recorded in the borehole after application of the time delay observed on first breaks, as previously described with regard to Equation 1. We refer to the upward continued data as a Huygens Stack.

(32) In particular, the present method creates pseudo-receiver positions at the surface of the earth for each of the source points (step 63 in FIG. 9). For example, FIG. 12 shows the map view of a line of pseudo-receiver positions, p, that extend from the receiver well in a direction that is co-linear with the source, S, and the top of the well location. The travel times from each position, p, to the receivers in the receiver well can be interpolated from first arrival travel times measured for source points (filled circles) that are nearby the position of each point, p. The synthetic seismogram at each location, p, is computed by Equation 1 where the time delay, Ti, which is the only unknown on the right hand side of Equation 1, is either known from measured first arrivals at each receiver location in the borehole, or can be accurately interpolated from the first arrival picks from source points that surround the pseudo-receiver position (step 64 in FIG. 9). The resulting synthetic seismograms can then be processed using conventional techniques as if they were actual surface seismic data (step 65 in FIG. 9).

(33) The seismograms created for the pseudo-receivers at the surface of the earth have the following attributes: a) They can be used as seismograms that were recorded at the surface of the earth and thus standard surface seismic data processing methods can be applied to the data. b) The seismograms have all static time delays that affect the travel time from the seismic source to the seismic receivers. Thus, surface-consistent data processing will yield consistent source and receiver statics. c) The reflection point between the source and receiver can be estimated based on the source-receiver position and the dip of beds in the same way that it is done for normal surface seismic data. This is one of the most important aspects of this technology given that estimating the reflection point for geophones in the borehole is a function of the unknown velocity field and the source and receiver position. When the receiver is near the same elevation as the source, the reflection point can initially be assumed to be half way between the source and receiver position and the unknown velocity field can be determined by standard NMO analysis.
Travel times delays that are applied as the Ti term in Equation 1 are not restricted to measured first arrival travel times. The value used in Ti can be computed from ray tracing or other travel time estimation techniques, such as picking first breaks from data produced from finite difference modeling.

(34) Pseudo-receivers can be placed in arbitrary positions relative to the real source and receiver position if the travel time term in Equation 1 can be either measured or estimated. The geophysical relevance of pseudo-receiver placement is dependent on the dataset and the geophysical method that is being applied. The geophysical relevance of the receiver position is related to the direction of propagating waves that are recorded in the borehole. Multi-component geophones that are typically used in modern VSP recording can be used to determine the direction of wavefield propagation and then be applied to determining the optimal placement of pseudo-receivers.

(35) Example of Upward Continuation by Travel Times. FIG. 13 shows a synthetic common source gather generated by the finite difference software to demonstrate the present method. The gather in FIG. 13 is one of a set of eleven finite difference shot records with varying horizontal source-receiver offsets that were computed to demonstrate the methods described. The upward continued version of this shot is shown in FIG. 14. Upward continuation was done by using Huygens' principle and Equation 1 described above. Source and receiver coordinates were assigned to the data traces for all upward continued common source gathers. Common depth point numbers were assigned to the traces of the twelve upward continued shots (not shown) via the standard method in surface seismic data processing. The data were sorted by CDP gather and normal moveout velocities were picked and applied to the dataset. The two-way time image (see FIG. 15) was then created by stacking the moveout-corrected CDPs. Travel times to the receivers were symmetric on either side of the well because the velocity model was one-dimensional and the receiver well was vertical.

(36) Method Using Finite Difference Upward Continuation. For the reasons mentioned above, upward continuation via finite difference is inferior to upward continuation by the summing method in Equation 1. Finite difference upward continuation methods for VSP processing have been reduced to practice as a result of our work. In spite of its inferiority, this invention should be understood to include the finite difference technique as a way of upward continuation of VSP data to the surface of the earth or to an elevation near the source elevation for the purpose of seismic reflection imaging of borehole seismic data.

(37) FIG. 16 is a diagram of the processes using finite difference upward continuation rather than the preferred method of upward continuation via first arrival times. Borehole seismic data is initially acquired as previously described (step 161 in FIG. 16). A set of pseudo-receiver positions are determined at the surface of the earth for each of the source points (step 162 in FIG. 16). As in the Huygens Stack method described above, the VSP common source gathers are filtered so as to extract the upward traveling seismic waves for upward continuation. The wavefield was then propagated through a velocity model to pseudo-receivers at an elevation that was equal to the seismic data processing datum elevation that was selected for processing (step 163 in FIG. 16). The elevation selected was the mean source elevation for the survey though any elevation can be selected and this invention should not be restricted to using the mean source elevation. As in the use of the stack method of Equation 1, the receiver coordinates were assigned to the pseudo-receivers at the surface of the earth and the data traces were input to a surface seismic data processing flow.

(38) FIG. 17 shows a finite difference common source gather. FIG. 18 is the shot record after upward continuation to the elevation of the source by finite difference methods. The trace with the smallest horizontal offset from the source point is shown on the right. Increasing horizontal source-receiver offset is to the left. FIG. 19 is the image produced after processing the synthetic VSP data after upward continuation with finite difference methods. As before, resulting synthetic seismograms can be processed using conventional techniques as if they were actual surface seismic data (step 164 in FIG. 16).

(39) General Application of Wave Equation Approximations. Upward continuation of VSP data by the Huygens Stack method or finite difference methods are all approximations of the general wave equation. All of these methods can be used to predict a wavefield at the surface of the earth based on seismic data that has been recorded in a borehole, which can then be treated as surface seismic data. It is possible that other approximations of the general wave equation could be employed for this purpose. The present invention should be construed to include other approximations of the general wave equations used to predict the wavefield at the surface of the earth, so that the wavefield can be treated as surface seismic data for purposes of subsequent seismic data processing.

(40) The above disclosure sets forth a number of embodiments of the present invention described in detail with respect to the accompanying drawings. Those skilled in this art will appreciate that various changes, modifications, other structural arrangements, and other embodiments could be practiced under the teachings of the present invention without departing from the scope of this invention as set forth in the following claims.