USH482H - Seismic migration method - Google Patents
Seismic migration method Download PDFInfo
- Publication number
- USH482H USH482H US06/920,552 US92055286A USH482H US H482 H USH482 H US H482H US 92055286 A US92055286 A US 92055286A US H482 H USH482 H US H482H
- Authority
- US
- United States
- Prior art keywords
- trace
- ref
- seismic
- traces
- migration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000013508 migration Methods 0.000 title claims abstract description 34
- 230000005012 migration Effects 0.000 title claims abstract description 29
- 238000012545 processing Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 abstract description 8
- 230000009466 transformation Effects 0.000 abstract description 7
- 238000000844 transformation Methods 0.000 abstract description 3
- 238000003672 processing method Methods 0.000 abstract 1
- 238000007796 conventional method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
Definitions
- This invention is a method for processing seismic data. More particularly, the invention is an improved seismic migration technique of the type in which seismic data are Fourier-Transformed (or subjected to cosine or sine transformations), the coefficients of the transformed data are thereafter processed ("migrated"), and the processed coefficients are thereafter inverse-Fourier-transformed (or subjected to inverse-cosine or inverse-sine transformations).
- the objective of all migration techniques is to generate a display of processed seismic data in which the position of each seismic reflection indicated by the display accurately represents the position within the earth of whatever caused the reflection.
- KF-migration In one conventional class of seismic migration techniques, known as "KF-migration" (or “FK-migration”) techniques, the seismic data to be processed are first Fourier-transformed, then the Fourier-coefficients are migrated, and thereafter, the migrated coefficients are inverse-Fourier-transformed.
- KF-migration methods are described in the following papers: R. H. Stolt, "Migration by Fourier Transform," Geophysics, V. 43, No. 1 (February, 1978), pp. 23-48; J. Gazdag, “Wave Equation Migration with the Phase-Shift Method," Geophysics, V. 43, No. 7 (December, 1978), pp. 1342-1351; and J. Gazdag, "Wave Equation Migration with the Accurate Space Derivative Method," Geophysical Prospecting, V. 28, (1980), pp. 60-70.
- a drawback of conventional KF-migration techniques of the type described in the Stolt paper is that they require that the relationship between depth of a reflector in the earth and the seismic wave traveltime to and from the reflector must be represented as a single fixed ratio. In other words, those methods require that the velocity of seismic waves in the earth be represented as a constant.
- the Gazdag (1978) technique has the disadvantage that it is time-consuming and expensive to implement, for a number of reasons, including the reason that the number of operations required to migrate each seismic trace in accord with the Gazdag (1978) technique is proportional to N 2 , where N is the number of data points (or samples) comprising the trace.
- the Gazdag (1980) technique too, has the disadvantage of being time-consuming and expensive to implement.
- the present invention represents an improvement over conventional KF-migration techniques because it facilitates efficient migration of seismic data without requiring that the velocity of the associated seismic waves be represented as a constant. Rather, the invention requires only that the earth be represented as a sequence of uniform horizontal layers, each layer having a characteristic velocity, or ratio of depth increment to traveltime.
- the number of operations required to migrate a single trace in accord with the present invention is proportional to N log N, where N is the number of data points (or samples) comprising the trace.
- the earth is modeled as a stack of M horizontal layers, each having, in general, a different seismic wave velocity.
- the sixth step of the preferred embodiment (“recursive KF migration”) is a sequence of operations performed M-1 times in succession (as an index j assumes values from 2 through M) for each transformed seismic trace.
- Recursive KF migration reassigns new values to the data points comprising the traces to approximate changes in the wave field as the associated source-receiver pair is "moved" downward into the earth. No interpolation is required in performing the recursive KF migration step of the invention.
- each data point of the trace is reassigned a single new value (that may be stored in a single computer storage location), and some of these new values are stored while the rest of these new values are discarded.
- the stored values remaining after all M-1 iterations are the result of applying the recursive KF migration operation to the relevant trace.
- steps 2 and 9 are omitted.
- the output traces produced in accord with any embodiment of the invention may be displayed in form meaningful to an interpreter using any of several suitable conventional methods and means.
- FIG. 1 is a diagram showing an unprocessed set of seismic traces, intermediate states of representative ones of these traces as they are processed in accord with the inventive method, and a set of traces in processed form after the invention has been applied thereto.
- FIG. 2 is a diagram showing the iterations performed in performing the step of the invention known as "recursive KF migration”.
- the seismic data to be processed in accordance with the preferred embodiment will be described with reference to Cartesian coordinate system (x,y,z), where the x-axis and y-axis are horizontal, and the z-axis is vertical.
- the seismic data are a set of seismic traces g(x,y,t i ). There is one trace for each horizontal position (x,y).
- An example of a set of traces to be processed in accord with the invention is set A of traces shown in FIG. 1. In FIG. 1, set A includes twenty traces. It should be appreciated that a more typical set of traces to be processed in accord with the invention will be a 256 trace by 512 trace array of 131,072 traces, or an even larger array of traces.
- the parameters t i for each trace comprise a set of N time instants separated by sampling interal ⁇ t.
- each stacked data point is associated with a single position (x,y) and a particular sample time t i
- trace g(x,y,t i ) is a stacked seismic trace associated with single position (x,y).
- the first step of the preferred embodiment of the invention is performed as follows. For each time sample (t i ), add together the depth values specified for each (x,y) mesh point and divide the sum by the number of addends. In other words, for each sample t i , generate ##EQU1## where the set of all mesh points (x,y) is ⁇ x mn , y mn :1 ⁇ m ⁇ N x , 1 ⁇ n ⁇ N y ⁇ , and z imn is the member of set z i associated with mesh point x mn , y mn . This produces the average time-depth relation, a set of depth values (z i av ) corresponding one to each time sample of the seismic traces.
- Each seismic trace of the input data is associated with its individual time-depth relationship ⁇ t i , z i :1 ⁇ i ⁇ N ⁇ defining the time and depth for each sample of the trace.
- the second step of the preferred embodiment is conforming the input seismic traces to the reference time-depth relationship. This second step includes the sub-step of assigning each sample a modified time t' i according to the scheme
- the set of conformed input traces ⁇ (i-1) ⁇ t, f i ⁇ may be represented as f[(m-1) ⁇ x, (n-1) ⁇ y, (i-1) ⁇ t]; where 1 ⁇ m ⁇ N x , 1 ⁇ n ⁇ N y , 1 ⁇ i ⁇ N.
- the Discrete Fourier Transform in any of steps 3, 4, and 5 may be performed using any of several suitable conventional implementations thereof. See, for example, G. D. Bergland, "A Guided Tour of the Fast Fourier Transform," IEEE Spectrum, 6, pp 41-52 (July, 1969) for a discussion of the Discrete Fast Fourier Transform.
- steps two, three and four are represented in FIG. 1 by generation of all the traces d[(m-1) ⁇ k x , (n-1) ⁇ k y , (i-1) ⁇ t] represented by single trace B from set A.
- Step five is represented in FIG. 1 by transformation of all traces represented by trace B into set of transformed traces c[(m-1) ⁇ k x , (n-1) ⁇ k y , (i-1) ⁇ w] represented by single signal C.
- the set of signals represented by signal C is next transformed into the set of numerical signals a[(m-1) ⁇ k x , (n-1) ⁇ k y , (i-1) ⁇ t] according to the "Recursive KF Migration" method described below.
- This recursive KF migration method (step six of the preferred embodiment) is a key element of all embodiments of the invention, and is represented in FIG. 1 by transformation of single signal C into signal D. In the FIG. 1 example, should it be understood that during step six, the recursive KF migration operation will be repeated to process each of the signals (all represented by signal C) resulting from performance of step five.
- the set of signals a is subjected to a Discrete Inverse Fourier Transform on index n to become set of signals h[(m-1) ⁇ k x , (n-1) ⁇ y, (i-1) ⁇ t].
- the Discrete Inverse Fourier Transform of this step seven (and of step eight below) may be performed using any of several suitable conventional implementations thereof.
- the set of numerical signals h is subjected to a Discrete Inverse Fourier Transform on index m to become set of signals q[(m-1) ⁇ x, (n-1) ⁇ y, (i-1) ⁇ t].
- the resulting set q of numerical values is a set of migrated conformed output traces.
- To select a pair of particular values (m,n) is to specify a particular migrated conformed output trace ⁇ (i-1) ⁇ t, q i :1 ⁇ i ⁇ N ⁇ associated with the time-depth relationship at mesh point ((m-1) ⁇ x, (n-1) ⁇ y).
- Performance of all nine steps of the preferred embodiment of the invention results in generation of a migrated output trace ⁇ (i-1) ⁇ t,p i :1 ⁇ i ⁇ N ⁇ for each of the mesh points (x,y) associated with an input trace.
- An example of a set of migrated output traces generated in accord with the invention is set E of traces shown in FIG. 1.
- the output traces produced in accord with the invention may be displayed in form meaningful to an interpreter using any of several suitable and well-known display methods and means.
- Step 5 precedes set 4 or step 3;
- Discrete Cosine Transforms or Discrete Sine Transforms are employed where Discrete Fourier Transforms have been indicated (any of several suitable conventional implementations of the Discrete Cosine Transforms or Discrete Sine Transform may be employed).
- Discrete Fourier Transform (or “Discrete Inverse Fourier Transform”) should be understood to denote any discrete transform selected from the group including Discrete Fourier Transforms, Discrete Cosine Transforms, and Discrete Sine Transforms (and the corresponding discrete inverse transformation).
- the sixth step of the preferred embodiment (“Recursive KF Migration") will next be described.
- the reference time-depth relation is defined by the set of numbers ⁇ t j ref , z j ref :1 ⁇ j ⁇ M, M ⁇ 2 ⁇ .
- the recursive KF migration step is a sequence of operations carried out M-1 times in succession, as index j assumes values from 2 through M. This process is carried out separately for every possible combination of m and n.
- Each performance of recursive KF migration on a trace reassigns a new value to each sample comprising the trace, to approximate changes in the seismic wave field as the source-receiver pair (or pairs) associated with index pair n,n are "moved" downward into the earth.
- the earth is modeled as a stack of M horizontal layers, each having, in general, a different seismic wave velocity, v j (to be defined precisely below).
- FIG. 2 shows the iteration cycles performed during application of recursive KF migration to a set of traces.
- the first operation is known as downward continuation.
- this sub-step is bypassed.
- the set of numbers c' is subjected to a Discrete Inverse Fourier Transform on index i' to become b[(m-1) ⁇ k x ,(n-1) ⁇ k y ,(i'-1) ⁇ t], 1 ⁇ i' ⁇ N.
- a[(m-1) ⁇ k x ,(n-1) ⁇ k y ,(i-1) ⁇ t] b[(m-1) ⁇ k.sub.x, (n-1) ⁇ k y ,(i'-1) ⁇ t], 1+Int[t j-1 ref / ⁇ t] ⁇ i ⁇ 1+Int[t j ref / ⁇ t].
- b[(m-1) ⁇ k x ,(n-1) ⁇ k y , (i-1) ⁇ t]] are stored.
- V max be defined as the maximum of the ratio (z j ref -z j-1 ref )/(t j ref -t j-1 ref ), considering all permissible values of j.
- V max be defined as the maximum of the ratio (z j ref -z j-1 ref )/(t j ref -t j-1 ref ), considering all permissible values of j.
- the seismic data to be processed in accord with the invention may be stacked or unstacked.
- the preferred embodiment described above is suitable for processing zero-offset (i.e., stacked) data.
- each input trace g(x s , y s , x R , y R , t i ) will be associated with both a source position (x s , y s ) and a receiver position (x R , y R ); not with a single position (x,y) representing both source and receiver positions.
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A seismic data processing method in which seismic traces (each trace containing a signal resulting from reflection of a seismic signal at a location within the earth, and each trace being associated with at least one point in a two-dimensional spatial grid (x,y)) are subjected to Fourier transformations, the coefficients of the Fourier-Transformed traces are subjected to a recursive KF migration operation, and the migrated traces are thereafter inverse-Fourier-transformed. When displayed, the processed seismic data accurately represents the position within the earth of whatever caused the reflection. The method may be employed to process stacked seismic traces, each associated with a single point (x,y) in the grid, or may be employed to process unstacked seismic traces, each associated with both a seismic source location (xs,y5) and a different seismic receiver location (xr,yr) in the grid. In performing the method, the earth is modeled as a stack of M horizontal layers each characterized by a seismic wave velocity. The recursive KF migration step is iterated M-1 times for each trace, where part of the output of each iteration is stored and part discarded.
Description
This invention is a method for processing seismic data. More particularly, the invention is an improved seismic migration technique of the type in which seismic data are Fourier-Transformed (or subjected to cosine or sine transformations), the coefficients of the transformed data are thereafter processed ("migrated"), and the processed coefficients are thereafter inverse-Fourier-transformed (or subjected to inverse-cosine or inverse-sine transformations).
There are several types of conventional techniques for migrating seismic data. The objective of all migration techniques is to generate a display of processed seismic data in which the position of each seismic reflection indicated by the display accurately represents the position within the earth of whatever caused the reflection.
In one conventional class of seismic migration techniques, known as "KF-migration" (or "FK-migration") techniques, the seismic data to be processed are first Fourier-transformed, then the Fourier-coefficients are migrated, and thereafter, the migrated coefficients are inverse-Fourier-transformed. Examples of such KF-migration methods are described in the following papers: R. H. Stolt, "Migration by Fourier Transform," Geophysics, V. 43, No. 1 (February, 1978), pp. 23-48; J. Gazdag, "Wave Equation Migration with the Phase-Shift Method," Geophysics, V. 43, No. 7 (December, 1978), pp. 1342-1351; and J. Gazdag, "Wave Equation Migration with the Accurate Space Derivative Method," Geophysical Prospecting, V. 28, (1980), pp. 60-70.
A drawback of conventional KF-migration techniques of the type described in the Stolt paper is that they require that the relationship between depth of a reflector in the earth and the seismic wave traveltime to and from the reflector must be represented as a single fixed ratio. In other words, those methods require that the velocity of seismic waves in the earth be represented as a constant.
The Gazdag (1978) technique has the disadvantage that it is time-consuming and expensive to implement, for a number of reasons, including the reason that the number of operations required to migrate each seismic trace in accord with the Gazdag (1978) technique is proportional to N2, where N is the number of data points (or samples) comprising the trace. The Gazdag (1980) technique, too, has the disadvantage of being time-consuming and expensive to implement.
The present invention represents an improvement over conventional KF-migration techniques because it facilitates efficient migration of seismic data without requiring that the velocity of the associated seismic waves be represented as a constant. Rather, the invention requires only that the earth be represented as a sequence of uniform horizontal layers, each layer having a characteristic velocity, or ratio of depth increment to traveltime. The number of operations required to migrate a single trace in accord with the present invention is proportional to N log N, where N is the number of data points (or samples) comprising the trace.
The preferred embodiment of the invention includes the following nine steps:
1. Determining the reference time-depth relation;
2. Conforming each seismic trace of the input data to the reference time-depth relation;
3. Performing a Discrete Fourier Transform over the x-coordinate of each conformed trace;
4. Performing a Discrete Fourier Transform over the y-coordinate of each conformed trace;
5. Performing a Discrete Fourier Transform on the time coordinate of each conformed trace;
6. Performing the process known as "Recursive KF Migration" on the data after they have been processed in accordance with steps 1-5;
7. Performing a Discrete Inverse Fourier Transform over the y-coordinate of each migrated data trace;
8. Performing a Discrete Inverse Fourier Transform over the x-coordinate of each migrated data trace; and
9. After processing the data in steps 1-8, restoring the processed data's individual time-depth relations.
The earth is modeled as a stack of M horizontal layers, each having, in general, a different seismic wave velocity. The sixth step of the preferred embodiment ("recursive KF migration") is a sequence of operations performed M-1 times in succession (as an index j assumes values from 2 through M) for each transformed seismic trace. Recursive KF migration reassigns new values to the data points comprising the traces to approximate changes in the wave field as the associated source-receiver pair is "moved" downward into the earth. No interpolation is required in performing the recursive KF migration step of the invention. Rather, during each of the M-1 iterations for each trace, each data point of the trace is reassigned a single new value (that may be stored in a single computer storage location), and some of these new values are stored while the rest of these new values are discarded. The stored values remaining after all M-1 iterations are the result of applying the recursive KF migration operation to the relevant trace.
In alternative embodiments, steps 2 and 9 are omitted. The output traces produced in accord with any embodiment of the invention may be displayed in form meaningful to an interpreter using any of several suitable conventional methods and means.
FIG. 1 is a diagram showing an unprocessed set of seismic traces, intermediate states of representative ones of these traces as they are processed in accord with the inventive method, and a set of traces in processed form after the invention has been applied thereto.
FIG. 2 is a diagram showing the iterations performed in performing the step of the invention known as "recursive KF migration".
The steps of the preferred embodiment will first be described in detail with reference to FIG. 1. Variations on this preferred embodiment will than be described. Thereafter, the sixth step of the preferred embodiment (known as "recursive KF migration") will be described in greater detail with reference to FIG. 2.
The seismic data to be processed in accordance with the preferred embodiment will be described with reference to Cartesian coordinate system (x,y,z), where the x-axis and y-axis are horizontal, and the z-axis is vertical. The seismic data are a set of seismic traces g(x,y,ti). There is one trace for each horizontal position (x,y). An example of a set of traces to be processed in accord with the invention is set A of traces shown in FIG. 1. In FIG. 1, set A includes twenty traces. It should be appreciated that a more typical set of traces to be processed in accord with the invention will be a 256 trace by 512 trace array of 131,072 traces, or an even larger array of traces. The parameters ti for each trace comprise a set of N time instants separated by sampling interal Δt. Each set ti may be described using the notation: {ti =iΔt:1≦i≦N}. It is within the scope of the invention to process data associated with only a single value of x (or y), or with many values of x and y (i.e., the method may be performed to process 2-D or 3-D data). For simplicity, it will be assumed for purposes of describing the preferred embodiment, that each trace represents the stacked electrical output of a number of acoustic receivers. For each instant ti, the outputs of the receivers are stacked together so that each stacked data point is associated with a single position (x,y) and a particular sample time ti, and trace g(x,y,ti) is a stacked seismic trace associated with single position (x,y). For a given trace, the sum of the receiver outputs at instant ti is represented by the numerical value gi =g(x,y,ti).
Despite the actual positions of the seismic source and receiver (or receivers) associated with each trace g(x,y,ti), the trace may be regarded as representing the output of an acoustic receiver at position (x,y) over an interval of time (0≦t≦NΔt) resulting from actuation of a seismic source at the same position (x,y) at t=0. It is presumed that each trace is indicative of at least one signal resulting from reflection by a subterranean reflector of the seismic wave proceeding outward from the source. It is further presumed that there exists a known or assumed relationship between the arrival time (ti) of each reflection at mesh point (x,y) and the depth (zi) of any real or hypothetical reflecting object directly below mesh point (x,y). The relationship between time and depth specified for each such regularly spaced time instant may be described using the notation {ti, zi :1≦i≦N}.
The first step of the preferred embodiment of the invention is performed as follows. For each time sample (ti), add together the depth values specified for each (x,y) mesh point and divide the sum by the number of addends. In other words, for each sample ti, generate ##EQU1## where the set of all mesh points (x,y) is {xmn, ymn :1≦m≦Nx, 1≦n≦Ny }, and zimn is the member of set zi associated with mesh point xmn, ymn. This produces the average time-depth relation, a set of depth values (zi av) corresponding one to each time sample of the seismic traces. From the resulting set of time-depth pairs {ti =ti av, zi av :1≦i≦N} select the reduced set {tj ref, zj ref :1≦j≦M, M<N}, constituting the reference time-depth relation, to satisfy the rules t1 =t1 ref =t1 ref =z1 av, |ti av -[tj-i ref +(zi av -zj-1 ref)(tj ref -tj-1 ref)/(zj ref -zj-1 ref)]|<E, where tj-1 ref <ti av ≦tj ref, and E is a selected acceptable tolerance. The feasibility of the invention proceeds from the observation that M may be of order 10 while N is of order 1000, if E is of order 10Δt.
Each seismic trace of the input data is associated with its individual time-depth relationship {ti, zi :1≦i≦N} defining the time and depth for each sample of the trace. The second step of the preferred embodiment is conforming the input seismic traces to the reference time-depth relationship. This second step includes the sub-step of assigning each sample a modified time t'i according to the scheme
t'.sub.i =t.sub.j-1.sup.ref +(t.sub.j.sup.ref -t.sub.j-1.sup.ref)(z.sub.i -z.sub.j-1.sup.ref)/(z.sub.j.sup.ref -z.sub.j-1.sup.ref),
where zj-1 ref ≦zi ≦zj ref. The resulting seismic trace {t'i, gi } is not sampled at regular time intervals. The second sub-step of the second inventive step is to construct from each trace {t'i, gi } the regularly sampled conformed trace {(i-1)Δt, fi } according to the scheme fi =gj-1 +(gj -gj-1) ((i-1)Δt-t'j-1)/(t'j -t'j-1), where t'j-1 ≦(i-1)Δt≦t'j.
The set of conformed input traces {(i-1)Δt, fi } may be represented as f[(m-1)Δx, (n-1)Δy, (i-1)Δt]; where 1≦m≦Nx, 1≦n≦Ny, 1≦i≦N. The third step of the preferred embodiment is to subject this set of Nx Ny N numbers to a Discrete Fourier Transform on index m to produce set of values e[(m-1)Δkx, (n-1)Δy, (i-1)Δt], where Δkx =1/(Nx Δx) and 1≦m≦Nx.
In the fourth step of the preferred embodiment, the set of signals e is subjected to a Discrete Fourier Transform on index n to become set of signals d[m-1)Δkx, (n-1)Δky, (i-1)Δt], where Δky =1/(Ny Δy) and 1≦n≦Ny.
In the fifth step of the preferred embodiment, the set of signals d is subjected to a Discrete Fourier Transform on index i to become set of signals c[(m-1)Δkx, (n-1)Δky, (i-1)Δw] where Δw=1/(N Δt) and 1≦i≦N. The Discrete Fourier Transform in any of steps 3, 4, and 5 may be performed using any of several suitable conventional implementations thereof. See, for example, G. D. Bergland, "A Guided Tour of the Fast Fourier Transform," IEEE Spectrum, 6, pp 41-52 (July, 1969) for a discussion of the Discrete Fast Fourier Transform.
The steps two, three and four are represented in FIG. 1 by generation of all the traces d[(m-1)Δkx, (n-1)Δky, (i-1)Δt] represented by single trace B from set A. Step five is represented in FIG. 1 by transformation of all traces represented by trace B into set of transformed traces c[(m-1)Δkx, (n-1)Δky, (i-1)Δw] represented by single signal C.
The set of signals represented by signal C is next transformed into the set of numerical signals a[(m-1)Δkx, (n-1)Δky, (i-1)Δt] according to the "Recursive KF Migration" method described below. This recursive KF migration method (step six of the preferred embodiment) is a key element of all embodiments of the invention, and is represented in FIG. 1 by transformation of single signal C into signal D. In the FIG. 1 example, should it be understood that during step six, the recursive KF migration operation will be repeated to process each of the signals (all represented by signal C) resulting from performance of step five.
In the seventh step of the preferred embodiment, the set of signals a is subjected to a Discrete Inverse Fourier Transform on index n to become set of signals h[(m-1)Δkx, (n-1)Δy, (i-1)Δt]. The Discrete Inverse Fourier Transform of this step seven (and of step eight below) may be performed using any of several suitable conventional implementations thereof.
In the eighth step of the preferred embodiment, the set of numerical signals h is subjected to a Discrete Inverse Fourier Transform on index m to become set of signals q[(m-1)Δx, (n-1)Δy, (i-1)Δt]. The resulting set q of numerical values is a set of migrated conformed output traces. To select a pair of particular values (m,n) is to specify a particular migrated conformed output trace {(i-1)Δt, qi :1≦i≦N} associated with the time-depth relationship at mesh point ((m-1)Δx, (n-1)Δy).
In the ninth step of the preferred embodiment, to compensate for the effect of step two described above, each sample is assigned a modified time ti according to the scheme ti =tj-1 +(tj -tj-1)(z'i -zj-1)/(zj -zj-1) where zj-1 <z'i ≦zj ; and {tj, zj :1≦j≦N} is the time-depth relationship at the mesh point; and z'i =zj-1 ref +(zj ref -zj-1 ref)((i-1)Δt-tj-1 ref)/(tj ref -tj-1 ref), with tj-1 ref <(i-1)Δt≦tj ref. The resulting migrated trace (ti, qi) is not sampled at regular time intervals. From it is constructed the regularly sampled migrated output trace {(i-1)Δt, pi } according to the scheme pi =qj-1 +(qj -qj-1)((i-1)Δt-tj-1)/(tj -tj-1). Steps seven, eight, and nine are represented in FIG. 1 by transformation of the set of signals represented by single signal D into set E, and will be carried out to process all the signals produced as a result of step six.
Performance of all nine steps of the preferred embodiment of the invention results in generation of a migrated output trace {(i-1)Δt,pi :1≦i≦N} for each of the mesh points (x,y) associated with an input trace. An example of a set of migrated output traces generated in accord with the invention is set E of traces shown in FIG. 1. The output traces produced in accord with the invention may be displayed in form meaningful to an interpreter using any of several suitable and well-known display methods and means.
The data processing steps described above should be translated into a series of specific computer instructions in a manner that will be apparent to those ordinarily skilled in the art of computer programming for seismic data processing. The method steps preferably should be performed using any of many commercially available, appropriately programmed digital computers.
It will be apparent to those of ordinary skill in the field of seismic data processing that numerous variations on the above-described preferred embodiment may be implemented. The inventors consider that it is still within the scope of the invention to vary the above-described method as follows:
(a) Step 5 precedes set 4 or step 3;
(b) Steps 3 and 4 are interchanged;
(c) Steps 7 and 8 are interchanged;
(d) Steps 2 and 9 are omitted;
(e) Nonlinear interpolation is performed where, for simplicity, linear interpolation has been indicated; or
(f) Discrete Cosine Transforms or Discrete Sine Transforms are employed where Discrete Fourier Transforms have been indicated (any of several suitable conventional implementations of the Discrete Cosine Transforms or Discrete Sine Transform may be employed).
Other variations may be within the scope of the invention. Throughout this application (including in the claims), the phrase "Discrete Fourier Transform" (or "Discrete Inverse Fourier Transform") should be understood to denote any discrete transform selected from the group including Discrete Fourier Transforms, Discrete Cosine Transforms, and Discrete Sine Transforms (and the corresponding discrete inverse transformation).
The sixth step of the preferred embodiment ("Recursive KF Migration") will next be described. In this description, the symbol WN ij will represent the complex number cos Q+√-1 sin Q, where Q=2πij/N. The reference time-depth relation is defined by the set of numbers {tj ref, zj ref :1≦j≦M, M≧2}. The recursive KF migration step is a sequence of operations carried out M-1 times in succession, as index j assumes values from 2 through M. This process is carried out separately for every possible combination of m and n. Each performance of recursive KF migration on a trace (identified by index pair m,n) reassigns a new value to each sample comprising the trace, to approximate changes in the seismic wave field as the source-receiver pair (or pairs) associated with index pair n,n are "moved" downward into the earth. The earth is modeled as a stack of M horizontal layers, each having, in general, a different seismic wave velocity, vj (to be defined precisely below). FIG. 2 shows the iteration cycles performed during application of recursive KF migration to a set of traces.
The first operation is known as downward continuation. In the first execution, when j=2, this sub-step is bypassed. In subsequent executions, each number c[m-1)Δkx, (n-1)Δky, (i-1)Δw] as it currently exists is multiplied by the corresponding factor WN pq where q=Int[(tj-1 ref -tj-2 ref)/Δt], understanding that the symbol Int[ ] indicates that the expression within the brackets is to be rounded to the nearest integer, and p=[(-1)2 -k2 Vj-1 ]1/2, Vj-1 =(zj-1 ref -zj-2 ref)/(tj-1 ref -tj-2 ref), k2 =(m-1)2 (Δkx /Δw)2 +(n-1)2 (Δky /Δw)2, provided that i≧1+kVj-1 ; otherwise, c[(m- 1)Δkx, (n-1)Δky, (i-1)Δw]=0.
The next operation is known as relinearization. For each sample index i, a quantity i' is computed according to i'Int[1+{(i-1)2 -k2 Vj 2 1/2 ], 1+kVj≦ i≦N/2+1, or i'=Int[N+1-{(N+1-i)2 -k2 Vj 2 }1/2 ], N/2+1<i≦N-kVj. The objective is to define c'[(m-1)Δkx, (n-1)Δky, (i'-1)Δw]=c[(m-1)Δkx,(n-1)Δky, (i-1)Δw]. If some particular integer value i' is not among those defined in this manner, the corresponding value of c' is taken to be zero. Various interpolation methods are permissible to distribute any sample of c among several samples of c' . The restriction that must be observed is that ##EQU2##
In the next operation, the set of numbers c' is subjected to a Discrete Inverse Fourier Transform on index i' to become b[(m-1)Δkx,(n-1)Δky,(i'-1)Δt], 1≦i'≦N.
In the next operation, one stores the partial result: a[(m-1)Δkx,(n-1)Δky,(i-1)Δt]=b[(m-1)Δk.sub.x, (n-1)Δky,(i'-1)Δt], 1+Int[tj-1 ref /Δt]≦i<1+Int[tj ref /Δt]. Thus, for each iteration associated with a particular value of j, some but not all of the values b[(m-1)Δkx,(n-1)Δky, (i-1)Δt]] are stored.
With reference to FIG. 2, the order in which different values of m and n are employed in performing recursive KF migration is immaterial. It is necessary that each possible combination of m and n be employed once and only once.
In the above description of the recursive KF migration step, let Vmax be defined as the maximum of the ratio (zj ref -zj-1 ref)/(tj ref -tj-1 ref), considering all permissible values of j. Before carrying out the recursive KF migration operation, it may be useful to set c[(m-1)Δkx, (n-1)Δky,(i-1)Δw]=0 for 1≦i≦kVmax /sin dmig, where k2 =(m-1)2 (Δkx /Δw)2, dmig is arbitrary but restricted by the requirement 0< sin dmig ≦1. It may furthermore be useful, before thus modifying the set of numbers c, to generate and set aside the set of numerical signals
c.sub.save (i)=c[(m-1)Δk.sub.x,(n-1)Δk.sub.y,(i-1)Δw]
for (kVmax /tan dnoise)≦i-1≦(kVmax /tan dmig) where dnoise is arbitrary but restricted by the requirement tan dmig ≦tan dnoise. Then, during recursive KF migration, each time c' is computed, csave is added to c', sample by sample, over the range of samples for which csave is defined.
Numerous other variations on the inventive embodiment described above are with the scope of the invention. For example, the seismic data to be processed in accord with the invention may be stacked or unstacked. The preferred embodiment described above is suitable for processing zero-offset (i.e., stacked) data. With unstacked data, each input trace g(xs, ys, xR, yR, ti) will be associated with both a source position (xs, ys) and a receiver position (xR, yR); not with a single position (x,y) representing both source and receiver positions. Accordingly, where unstacked data are to be processed in accordance with the invention, modifications to the above-described must be made in a manner that will be apparent to one of ordinary skill in the seismic data processing art. In particular, the relations employed for downward continuation and relinearization will change to reflect what is known to those of ordinary skill in the art as the "Double Square Root" equation.
The above description is merely illustrative of the invention. It is contemplated that various changes in the details of the structures and methods described may be within the scope of the invention as defined by the appended claims.
Claims (4)
1. A method for processing a set of seismic traces, each trace corresponding to at least one location in a two-dimensional spatial grid (x,y), and each trace including a set of N time samples corresponding to time instants {ti =iΔt:1≦i≦N}, given a time-depth relationship {ti,zi :1≦i≦N} for each trace, said method including the steps of:
(a) determining a reduced set of time-depth pairs {tj ref,zj ref :1≦j≦M, M<N}, where t1 =t1 av =t1 ref, z1 ref =z1 av, |ti av -[tj-1 ref +(zi av -zj-1 ref)(tj ref -tj-1 ref)/(zj ref -zj-1 ref)]|<E, where tj-1 ref <ti av ≦tj ref, zi av is determined by adding together the depth values zi for each trace and dividing by the number of traces, and E is a selected acceptable tolerance on the order of about 10Δt;
(b) subjecting each trace to Discrete Fourier Transforms on the spatial coordinates corresponding to the trace, and on time index i;
(c) performing recursive KF migration on the traces after they have been transformed in step (b); and
(d) subjecting each migrated trace generated in step (c) to Discrete Inverse Fourier Transforms on the spatial coordinates corresponding thereto.
2. The method of claim 1, where each trace is a stacked trace g(x,y,ti), corresponding to a single point (xmn, ymn) in grid (x,y), and (xmn, ymn) is determined by the position coordinates of the seismic sources and receivers associated with the individual traces comprising the stacked trace.
3. The method of claim 1, where each trace is an unstacked trace g(xs, ys, xR, yR, ti) corresponding to a seismic source position (xs, ys) and to an acoustic receiver position (xR, yR).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US06/920,552 USH482H (en) | 1986-10-20 | 1986-10-20 | Seismic migration method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US06/920,552 USH482H (en) | 1986-10-20 | 1986-10-20 | Seismic migration method |
Publications (1)
Publication Number | Publication Date |
---|---|
USH482H true USH482H (en) | 1988-06-07 |
Family
ID=25443938
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US06/920,552 Abandoned USH482H (en) | 1986-10-20 | 1986-10-20 | Seismic migration method |
Country Status (1)
Country | Link |
---|---|
US (1) | USH482H (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4980866A (en) * | 1989-11-30 | 1990-12-25 | Conoco Inc. | Common offset depth migration with residual moveout correction |
US4992996A (en) * | 1990-01-31 | 1991-02-12 | Conoco Inc. | Interval velocity analysis and depth migration using common reflection point gathers |
US5128899A (en) * | 1991-07-16 | 1992-07-07 | Conoco Inc. | Method for sharpening prestack depth migration images by removing residual moveout from common image point gathers before stacking |
US5138584A (en) * | 1989-09-06 | 1992-08-11 | Chevron Research & Technology Company | Migration of seismic turning waves |
WO1993013434A1 (en) * | 1991-12-20 | 1993-07-08 | Schlumberger Technology Corporation | Method of seismic migration using a massively parallel computer |
US5586026A (en) * | 1991-12-20 | 1996-12-17 | Schlumberger Technology Corporation | Method of seismic migration using a massively parallel computer |
US5629904A (en) * | 1994-11-30 | 1997-05-13 | Paradigm Geophysical, Ltd. | Migration process using a model based aperture technique |
US5995904A (en) * | 1996-06-13 | 1999-11-30 | Exxon Production Research Company | Method for frequency domain seismic data processing on a massively parallel computer |
US6049759A (en) * | 1998-01-16 | 2000-04-11 | Bp Amoco Corporation | Method of prestack 3-D migration |
US20110093203A1 (en) * | 2009-10-21 | 2011-04-21 | Chevron U.S.A. Inc. | System and method for clustering arrivals of seismic energy to enhance subsurface imaging |
-
1986
- 1986-10-20 US US06/920,552 patent/USH482H/en not_active Abandoned
Non-Patent Citations (6)
Title |
---|
Berkhout et al., "Recursive Migration in Three Dimensions," Geophysical Prospecting, vol. 29, #5, pp. 758-781, Oct. '81. |
Berryhill, J. R., "Wave-Equation Datuming," Geophysics, vol. 44, #8, Aug. 1979. |
Dubrulle, A. A., "On Numerical Methods for Migration in Layered Media," Geophysical Prospecting, vol. 31, pp. 237-264, 1983. |
J. Gazdag, "Wave Equation Migration with the Accurate Space Derivative Method", Geophysical Prospecting, vol. 28, pp. 60-70 (1980). |
J. Gazdag, "Wave Equation Migration with the Phase-Shift Method," Geophysics, vol. 43, No. 7, pp. 1342-1351 (Dec., 1978). |
R. H. Stolt, "Migration by Fourier Transform," Geophysics, vol. 43, No. 1, pp. 23-48 (Feb., 1978). |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5138584A (en) * | 1989-09-06 | 1992-08-11 | Chevron Research & Technology Company | Migration of seismic turning waves |
US4980866A (en) * | 1989-11-30 | 1990-12-25 | Conoco Inc. | Common offset depth migration with residual moveout correction |
US4992996A (en) * | 1990-01-31 | 1991-02-12 | Conoco Inc. | Interval velocity analysis and depth migration using common reflection point gathers |
US5128899A (en) * | 1991-07-16 | 1992-07-07 | Conoco Inc. | Method for sharpening prestack depth migration images by removing residual moveout from common image point gathers before stacking |
WO1993013434A1 (en) * | 1991-12-20 | 1993-07-08 | Schlumberger Technology Corporation | Method of seismic migration using a massively parallel computer |
US5349527A (en) * | 1991-12-20 | 1994-09-20 | Schlumberger Technology Corporation | Method of seismic time migration using a massively parallel computer |
AU670256B2 (en) * | 1991-12-20 | 1996-07-11 | Geco A.S. | Method of seismic migration using a massively parallel computer |
US5586026A (en) * | 1991-12-20 | 1996-12-17 | Schlumberger Technology Corporation | Method of seismic migration using a massively parallel computer |
US5629904A (en) * | 1994-11-30 | 1997-05-13 | Paradigm Geophysical, Ltd. | Migration process using a model based aperture technique |
US5995904A (en) * | 1996-06-13 | 1999-11-30 | Exxon Production Research Company | Method for frequency domain seismic data processing on a massively parallel computer |
US6049759A (en) * | 1998-01-16 | 2000-04-11 | Bp Amoco Corporation | Method of prestack 3-D migration |
US20110093203A1 (en) * | 2009-10-21 | 2011-04-21 | Chevron U.S.A. Inc. | System and method for clustering arrivals of seismic energy to enhance subsurface imaging |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2318224C (en) | Seismic signal processing method and apparatus | |
CA2126325C (en) | Method of seismic migration using a massively parallel computer | |
US5995904A (en) | Method for frequency domain seismic data processing on a massively parallel computer | |
Symes et al. | Inversion of reflection seismograms by differential semblance analysis: Algorithm structure and synthetic examples 1 | |
US5995905A (en) | Source signature determination and multiple reflection reduction | |
US4380059A (en) | F-K Filtering of multiple reflections from a seismic section | |
US4628492A (en) | Method of avoiding aliasing in slant stacking of seismic data | |
US5500832A (en) | Method of processing seismic data for migration | |
US20090292476A1 (en) | Method of seismic data interpolation by projection on convex sets | |
USH482H (en) | Seismic migration method | |
Hindriks et al. | Reconstruction of 3-D seismic signals irregularly sampled along two spatial coordinates | |
US5079703A (en) | 3-dimensional migration of irregular grids of 2-dimensional seismic data | |
EP1198776A1 (en) | Seismic signal processing method and apparatus for generating a cube of variance values | |
US8213260B2 (en) | Fast residual migration of seismic data through parsimonious image decomposition | |
CN110058303A (en) | Acoustic anisotropy reverse-time migration mixed method | |
Van Borselen et al. | Removal of surface-related wave phenomena—The marine case | |
AU707757B2 (en) | Method for frequency domain seismic data processing on a massively parallel computer | |
US5586026A (en) | Method of seismic migration using a massively parallel computer | |
US8010293B1 (en) | Localized seismic imaging using diplets | |
US7069200B2 (en) | Method intended to obtain reflection travel times from an interpretation of migrated cylindrical wave seismic data | |
US6032104A (en) | Method of processing seismic signals so as to reduce the effects of dip in the analysis of moveouts | |
Bube et al. | Simultaneous determination of source wavelet and velocity profile using impulsive point-source reflections from a layered fluid | |
GB2082771A (en) | F-K Geophysical operations including filtering of seismic records | |
Bregman et al. | A noniterative procedure for inverting plane‐wave reflection data at several angles of incidence using the Riccati equation | |
AU712392B2 (en) | Geophysical attribute migration method for a medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: EXXON PRODUCTION RESEARCH COMPANY, A CORP OF DE Free format text: ASSIGNMENT OF ASSIGNORS INTEREST.;ASSIGNORS:BERRYHILL, JOHN R.;GONZALEZ, RUTH;KIM, YOUNG C.;REEL/FRAME:004647/0736 Effective date: 19861015 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |