Complex Analysis of Kinematics for Non-Hyperbolic Moveout Corrections

Abstract

A computer-implemented method for processing data includes receiving a collection of traces corresponding to signals received over time at multiple locations due to reflection of seismic waves from subsurface structures. A measure of correlation among the traces as is computed a function of a set of wavefront parameters, which determine respective moveout corrections to be applied in aligning the traces. A matrix having at least three dimensions is generated, wherein the elements of the matrix include the computed measure of the correlation. Using the matrix, values of the wavefront parameters are identified automatically or interactively along the time axis or along selected horizons to maximize the measure of the correlation, and a seismic image of the subsurface structures is generated by aligning and integrating the traces using the moveout corrections that are determined by the identified values of the wavefront parameters.

1 . A computer-implemented method for processing data, comprising:
 receiving a collection of traces corresponding to signals received over time at multiple locations due to reflection of seismic waves from subsurface structures;   computing a measure of correlation among the traces as a function of a set of wavefront parameters, which determine respective moveout corrections to be applied in aligning the traces;   generating a matrix having at least three dimensions, wherein one of the dimensions corresponds to propagation times of the seismic waves and at least two of the dimensions correspond respectively to at least two of the wavefront parameters, and the elements of the matrix comprise the computed measure of the correlation;   identifying, using the matrix, values of the wavefront parameters that maximize the measure of the correlation for each of a plurality of propagation times, respectively; and   generating a seismic image of the subsurface structures by aligning and integrating the traces using the moveout corrections that are determined by the identified values of the wavefront parameters.   
   
   
       2 . The method according to  claim 1 , wherein receiving the collection of traces comprises receiving a supergather of the traces, comprising multiple gathers of the traces having different, respective sources and receivers in a two-dimensional or three-dimensional acquisition geometry, and wherein identifying the values comprises processing the measure of the correlation computed over the multiple gathers in order to identify the values of the wavefront parameters that maximize the measure of the correlation. 
   
   
       3 . The method according to  claim 2 , wherein receiving the supergather comprises selecting multiple gathers, each comprising a respective set of the traces for which the sources and receivers define a different acquisition geometry, and wherein generating the seismic image of the subsurface structures comprises computing respective seismic images in two or three dimensions corresponding to the gathers. 
   
   
       4 . The method according to  claim 1 , wherein the moveout corrections are determined by an approximation that is selected from a set of approximations consisting of hyperbolic, non-hyperbolic and parabolic approximations. 
   
   
       5 . The method according to  claim 1 , wherein the at least two of the wavefront parameters comprise an angular parameter and at least one radius of curvature of a wavefront outgoing from a subsurface reflector. 
   
   
       6 - 9 . (canceled) 
   
   
       10 . The method according to  claim 1 , wherein the at least two of the wavefront parameters comprise an angular parameter and at least one of a wavefront curvature parameter and a velocity parameter. 
   
   
       11 . The method according to  claim 1 , wherein identifying the values comprises finding one or more respective maxima of the computed measure of the correlation for each of the propagation times, and wherein aligning the traces comprises aligning the traces according to the values of the wavefront parameters corresponding to the respective maxima. 
   
   
       12 . The method according to  claim 11 , wherein finding the one or more respective maxima comprises sorting the maxima in order of the computed measure of the correlation at each of the maxima. 
   
   
       13 . (canceled) 
   
   
       14 . The method according to  claim 1 , wherein identifying the values comprises displaying an output that shows the computed measure of the correlation as a function of the propagation times and at least one of the wavefront parameters, and receiving an input from a user that identifies one or more maxima of the computed measure of the correlation. 
   
   
       15 . The method according to  claim 14 , wherein receiving the input comprises receiving a selection by the user of a corridor containing a plurality of the maxima at different propagation times and corresponding values of the at least one of the wavefront parameters. 
   
   
       16 . The method according to  claim 15 , wherein the maxima in the corridor correspond to waves reflected respectively from successive subsurface layers. 
   
   
       17 . The method according to  claim 1 , and comprising defining a spatial horizon within a region of the seismic image, and identifying, using the matrix, values of the wavefront parameters that maximize the measure of the correlation for each of a plurality of locations along the spatial horizon, respectively. 
   
   
       18 . Apparatus for processing data, comprising:
 an interface, which is coupled to receive a collection of traces corresponding to signals received over time at multiple locations due to reflection of seismic waves from subsurface structures; and   a signal processor, which is configured to compute a measure of correlation among the traces as a function of a set of wavefront parameters, which determine respective moveout corrections to be applied in aligning the traces, to generate a matrix having at least three dimensions, wherein one of the dimensions corresponds to propagation times of the seismic waves and at least two of the dimensions correspond respectively to at least two of the wavefront parameters, and the elements of the matrix comprise the computed measure of the correlation, to identify, using the matrix, values of the wavefront parameters that maximize the measure of the correlation for each of a plurality of propagation times, respectively, and to generate a seismic image of the subsurface structures by aligning and integrating the traces using the moveout corrections that are determined by the identified values of the wavefront parameters.   
   
   
       19 . The apparatus according to  claim 18 , wherein the collection of the traces comprises a supergather of the traces, comprising multiple gathers of the traces having different, respective sources and receivers in a two-dimensional or three-dimensional acquisition geometry, and wherein the signal processor is configured to process the measure of the correlation computed over the multiple gathers in order to identify the values of the wavefront parameters that maximize the measure of the correlation. 
   
   
       20 . The apparatus according to  claim 18 , wherein the moveout corrections are determined by an approximation that is selected from a set of approximations consisting of hyperbolic, non-hyperbolic and parabolic approximations. 
   
   
       21 . The apparatus according to  claim 18 , wherein the at least two of the wavefront parameters comprise an angular parameter and at least one radius of curvature of a wavefront outgoing from a subsurface reflector. 
   
   
       22 . The apparatus according to  claim 18 , wherein the at least two of the wavefront parameters comprise an angular parameter and at least one of a wavefront curvature parameter and a velocity parameter. 
   
   
       23 . The apparatus according to  claim 18 , wherein the signal processor is configured to find one or more respective maxima of the computed measure of the correlation for each of the propagation times, and to align the traces according to the values of the wavefront parameters corresponding to the respective maxima. 
   
   
       24 . The apparatus according to  claim 18 , wherein the signal processor is configured to display an output that shows the computed measure of the correlation as a function of the propagation times and at least one of the wavefront parameters, and to receive an input from a user that identifies one or more maxima of the computed measure of the correlation. 
   
   
       25 . A computer software product, comprising a computer-readable medium in which program instructions are stored, which instructions, when read by a computer, cause the computer to receive a collection of traces corresponding to signals received over time at multiple locations due to reflection of seismic waves from subsurface structures, to compute a measure of correlation among the traces as a function of a set of wavefront parameters, which determine respective moveout corrections to be applied in aligning the traces, to generate a matrix having at least three dimensions, wherein one of the dimensions corresponds to propagation times of the seismic waves and at least two of the dimensions correspond respectively to at least two of the wavefront parameters, and the elements of the matrix comprise the computed measure of the correlation, to identify, using the matrix, values of the wavefront parameters that maximize the measure of the correlation for each of a plurality of propagation times, respectively, and to generate a seismic image of the subsurface structures by aligning and integrating the traces using the moveout corrections that are determined by the identified values of the wavefront parameters.