Abstract: Techniques for generating a geological model from 3D seismic data and rock property data are disclosed. Rock property data and 3D seismic data are received. Based on the rock property data and the 3D seismic data, an adaptive geological model is generated. The adaptive geological model includes a characteristic geological property. Synthetic seismic data is generated from a first region of interest of the adaptive geological model. The synthetic seismic data is adapted to facilitate a comparison between the first region of interest and a corresponding region of interest of the received 3D seismic data. The characteristic geological property is adjusted until the comparison indicates a result that is within a predetermined threshold region of the corresponding value from the rock properties. A validated geologic model is then generated.

Abstract: The present invention relates to a power supply device for supplying power, from a battery provided at an electric vehicle, to various types of components, the power supply device enabling current sensing to be performed so as to confirm that an electric vehicle charging facility and an electric vehicle are connected to each other through a connector, and reducing electro-magnetic interference (EMI) noise. The objective of the power supply device of the present invention is to: provide a stable voltage to electric vehicle application components; and reduce the EMI noise according to battery charging, self-switching of the power supply device and the like.

Abstract: Methods of analyzing and optimizing a seismic survey design are described. Specifically, the sampling quality is analyzed as opposed to the overall quality of the whole survey. This allows for analysis of the impact of the offsets, obstacles, and other aspects of the survey on the sampling quality, which will improve the ability to compress the resulting data and minimize acquisition footprints.

Abstract: A system and methods for analyzing seismic data are provided herein. The method includes identifying, via a computing device, a representation of a seismic data set (1802) and determining a number of feature descriptors corresponding to each of a number of aggregates within the representation (1804). The method also includes identifying a query relating to the representation and one or more vocabulary definitions relating to the query (1806), analyzing the representation to compute a likelihood that each of the aggregates satisfies the query (1808), and returning a result of the query (1810).

Abstract: Disclosed are an automatic focus identification method and a system for a Karst cave reservoir. According to the method, velocity replacement is performed based on a pre-stack time migration velocity model to obtain a plurality of replacement migration velocity volumes; pre-stack time migration is performed based on the plurality of replacement migration velocity volumes to obtain a plurality of post-stack seismic data volumes; a three dimensional cube window is selected based on a seismic response characteristic of the Karst cave reservoir; sparse statistical data of a seismic trace or sparse statistical data of a seismic trace attribute of each of the post-stack seismic data volumes in the three dimensional cube window are obtained respectively; and an optimized migration velocity and an optimized seismic trace are obtained through optimization determination.

Abstract: One embodiment of the present disclosure includes a method for processing seismic data comprising the steps of receiving data representing seismic energy gathered from a formation by a plurality of seismic receivers, wherein the data include primary and multiple data. A copy of the received data is created and compensated to reduce amplitude attenuation effects due to transmission and absorption losses. A multiple prediction algorithm is applied to the received and compensated data to obtain a multiple data prediction. The multiple data prediction is subtracted from the received data to obtain primary data. The primary data is processed to reduce attenuation effects in the received data.

Abstract: The method and system for measuring low-noise acoustical impulse responses at high sampling rates of the present invention utilizes two exponential sine sweeps (ESSs) to measure the impulse responses. The first ESS is a quick sweep up to the Nyquist frequency to provide an estimate of the system response and sample the ambient noise. This measurement is used to algorithmically determine an appropriate pass-band of the system. A second, slower sweep through the pass-band alone is then executed and a corresponding band-pass filter is applied to the resulting output signal to suppress noise. The result is a measured impulse response with an improved signal-to-noise ratio and a much-reduced pre-response.

Abstract: A sweep generator is employed to generate a sweep to be used by a seismic vibrator device for generating a desired target output spectrum, wherein the frequency sweep is designed so as to comply with one or more constraints imposed by the seismic vibrator device and/or imposed by the environment in which the device is to be used. In one embodiment, a sweep generator determines a sweep for achieving a desired target output spectrum by a given seismic vibrator device in compliance with at least a pump flow constraint imposed by the seismic vibrator device. In another embodiment, a sweep generator determines a sweep for achieving a desired target output spectrum by a given seismic vibrator device in compliance with all of multiple operational constraints of the seismic vibrator device, such as both mass displacement and pump flow constraints. Environmental constraints may also be accounted for in certain embodiments.

Abstract: A series of scans is generated for a subsurface and a derivative image is created using the series of subsurface images. One or more tests are performed on the derivative image, and a subsurface object is detected based on the one or more tests. A sensor is configured to generate a series of scans for a subsurface and a processor is coupled to the sensor. The processor is configured to execute stored program instructions that cause the processor to generate a series of images of the subsurface using the series of scans, create a derivative image using the series of subsurface images, perform one or more tests on the derivative image, and detect a subsurface object based on the one or more tests.

Abstract: A voice signal processor detects background noise sections to reflect characteristics of the background noise on the Wiener filter coefficient to be used for suppressing noise components of input voice signals. In the voice signal processor, directivity signal generators form directivity signals having a directivity pattern. The directivity signals are used by a coherence calculator to obtain coherence, which is in turn used by a targeted voice section detector to detect a targeted voice section. A background noise section detector detects background noise sections containing no voice signal. When a background noise section is detected, a WF adapter uses characteristics of background noise in the detected temporal section to calculate a new WF coefficient.

Abstract: A technique includes receiving data corresponding to at least one measurement of energy produced by a shot of an energy source; and transforming the received data to a domain that is sensitive to event direction. The technique includes comparing the measurement(s) with a reference based at least in part on the transformation; and based at least in part on the comparison, attenuating noise from the measurement(s).

Abstract: There is provided a method of interpolating wave signals, particularly seismic signals acquired through a seismic survey, using the steps of obtaining time series of measured wave signals; and selecting iteratively basis functions to represent said measured signals with a basis function being fully defined by n parameters, wherein in each iteration one or more basis functions are combined with the selected basis functions such that the residual between the measured signals and a representation of the measured signal by the combined basis functions is minimized at the locations of the measurement at each iteration.

Abstract: A technique includes processing seismic data indicative of samples of at least one measured seismic signal in a processor-based machine to, in an iterative process, determine basis functions, which represent a constructed seismic signal. The technique includes in each iteration of the iterative process, selecting another basis function of the plurality of basis functions. The selecting includes based at least in part on the samples and a current version of the constructed seismic signal, determining a cost function; and interpreting the cost function based at least in part on a predicted energy distribution of the constructed seismic signal to select the basis function.

Abstract: Adaptive filtering method to remove ground roll from seismic data. In an M channel adaptive filter, weights Wi are set using an adaptive algorithm based on seeking the minimum in the partial differential of cost function J. The cost function includes an expansion of the primary trace d into d=dg+?d (where: dg is ground roll contribution and ?d=dsig+dran, where dsig is the reflected signal component and dran is a random noise component) and a corresponding expansion of the reference x into x=xg+?x (where xg is ground roll contribution and ?x=xsig+xran; where xsig is a reflected signal component and xran is a random noise component). The delta components are included in the denominator of cost function J to provide an optimal solution of the filter coefficients biased by the reflection signal and random noise is removed.

Abstract: A technique includes spatially filtering a signal that is derived from a seismic acquisition. The filtering is associated with a filter length, and the filtering includes varying the filter length with frequency. The filtering may be used in connection with adaptive noise attenuation, which is applied to decomposed subbands. Furthermore, the filtering may be applied during the reconstruction of the signal from the subbands.

Abstract: The present invention relates to a method of filtering seismic data for noise attenuation. An embodiment of the present invention provides a method of processing seismic data in which the seismic data is transformed into an f-x domain using a discrete Fourier transform and is then filtered at each discrete frequency using an infinite impulse response (HR) filter.

Abstract: Methods and systems utilizing seismic sensors configured or designed for use in seismic signal detection are provided so as to reduce the occurrence of spurious responses of the sensors. A method of seismic surveying using a seismic sensor may include the steps of deploying the seismic sensor at a location for seismic signal detection and acquiring seismic signals. The seismic signals may include high frequency response signals containing spurious response signals at an identifiable bandwith. The method may further include applying spurious response cancellation based on the bandwidth location of the spurious response signals and generating modified seismic waveforms having extended frequency bandwidth.

Abstract: A frequency-domain based echo and NEXT canceller is provided. The canceller uses log2 encoding to precondition the error signal representing the echo. An improved gradient constraint is applied on at least a portion of a full weight vector in a least-mean-square algorithm. The least-mean-square algorithm is used to compute filter coefficients. The filter coefficients are multiplied by a frequency-domain data vector using a frequency-domain multiplier to generate frequency-domain output vector.

Abstract: A sweep generator is employed to generate a sweep to be used by a seismic vibrator device for generating a desired target output spectrum, wherein the frequency sweep is designed so as to comply with one or more constraints imposed by the seismic vibrator device and/or imposed by the environment in which the device is to be used. In one embodiment, a sweep generator determines a sweep for achieving a desired target output spectrum by a given seismic vibrator device in compliance with at least a pump flow constraint imposed by the seismic vibrator device. In another embodiment, a sweep generator determines a sweep for achieving a desired target output spectrum by a given seismic vibrator device in compliance with all of multiple operational constraints of the seismic vibrator device, such as both mass displacement and pump flow constraints. Environmental constraints may also be accounted for in certain embodiments.

Abstract: The invention is a method to predict surface-wave waveforms (306) and subtract them (307) from seismic data. Prediction is done by estimating a set of surface-consistent components (transfer functions in the frequency domain or impulse responses in time domain) that best represent changes in the waveforms for propagation along the surface from source to receiver (303). The prediction uses a mathematical expression, or model, of the earth's filtering effects, both amplitude and phase, as a function of frequency. The desired surface-consistent components are model parameters, and model optimization is used to solve for the surface-consistent components. The surface-consistent components may include filter transfer functions for each source location, each receiver location, and for propagation (302) through each region (301) of the surface that exhibits lateral variation.

Abstract: Measured seismic data is received from a seismic sensor. Rotation data is also received, where the rotation data represents rotation with respect to at least one particular axis. The rotation data is combined, using adaptive filtering, with the measured seismic data to attenuate at least a portion of a noise component from the measured seismic data.

Abstract: A method for seismic event mapping includes adaptively velocity filtering seismic signals recorded at selected positions. The velocity filtered signals are transformed into a domain of possible spatial positions of a source of seismic events. An origin in spatial position and time of at least one seismic event is determined from space and time distribution of at least one attribute of the transformed seismic data.

Abstract: A method and apparatus for use in seismic prospecting are disclosed. The method includes: attenuating coherent noise in a seismic data set generated from a sweep signal; and correlating the coherent noise-attenuated seismic data set to the sweep signal. The apparatus includes a program storage medium encoded with instructions that, when executed by a computing device perform the method and a computer programmed to perform the method.

Abstract: Method for processing seismic data sets corresponding to the same zone, some of which are acquired by means of receivers at the ocean bottom, while others are acquired by means of surface receivers, characterized in that at least one seismic data set that is an estimation of the component common to at least two seismic data sets corresponding respectively to one and to the other of these two acquisition modes is determined.

Abstract: A method waveform processing technique utilizing signal coherence of the array data for processing signals having poor signal-to-noise ratio. Raw waveform data is first transformed into f-k (frequency-wavenumber) domain. A coherence function is then calculated and convolved with the data in the f-k domain, which effectively suppresses non-coherent signals in the data. For the remaining coherent data, the unwanted part is muted and the wanted part is retained and inverse-transformed to yield the coherence-filtered array waveform data. After this processing, small signals that are hidden in the original data are extracted with much enhanced coherence. Subsequent processing of the data yields reliable information about formation acoustic property.

Abstract: A method and apparatus for use in seismic prospecting are disclosed. The method includes: attenuating coherent noise in a seismic data set generated from a sweep signal; and correlating the coherent noise-attenuated seismic data set to the sweep signal. The apparatus includes a program storage medium encoded with instructions that, when executed by a computing device perform the method and a computer programmed to perform the method.

Abstract: A statistical value update unit calculates the fluctuation of measurement data. A filtering processing unit extracts a normal white noise component from the fluctuation of the measurement data using an adaptive lattice filter. A statistical test unit determines whether or not the variance of the normal white noise component is out of a predetermined scope in reference distribution. A change decision unit detects a stationary change of the status of a target system based on a detection ratio of an outlier.

Abstract: A method and apparatus for use in seismic prospecting are disclosed. The method includes: attenuating coherent noise in a seismic data set generated from a sweep signal; and correlating the coherent noise-attenuated seismic data set to the sweep signal. The apparatus includes a program storage medium encoded with instructions that, when executed by a computing device perform the method and a computer programmed to perform the method.

Abstract: To generate an absorption parameter model, estimated values of an effective absorption parameter are received, where the estimated effective absorption parameter values represent absorption encountered by a seismic wave in a subterranean structure. Based on the estimated effective absorption parameter values, an absorption parameter model is generated that varies absorption parameter values along at least one dimension of the subterranean structure.

Abstract: Methods, algorithms, software, architectures, systems and circuits for targeting certain dominant error types in an adaptive FIR filter and/or signal equalizer. The method and algorithm generally include processing a data sequence in accordance with the adaptive algorithm to produce a processed data sequence, filtering the data sequence to generate a filtered data term for the adaptive algorithm, generating a filtered error term for the adaptive algorithm from at least the processed data sequence, and updating the adaptive algorithm in response to the filtered data and error terms. The architectures generally include an equalizer configured to equalize and/or filter a data sequence and provide an equalized data output, a first filter configured to receive the data sequence and generate a filtered data term for the adaptive algorithm, and an error term circuit configured to receive the equalized data output and provide a filtered error term for the adaptive algorithm.

Abstract: A method and apparatus for seismic trace matching to well logs.

Abstract: Coherent wave noise energy is removed from seismic data by modeling both the P-wave primary energy and the coherent wave noise energy. The P-wave primary energy is modeled first and then subtracted from the input data. The data with the P-wave primary energy removed is used as the input for coherent wave energy removal. The coherent wave energy is modeled and subtracted from the original input data, i.e. the data input into P-wave primary removal. This leaves a dataset with P-wave primary energy and noise energy not related to coherent waves. This method can be utilized to remove all types of coherent noise with a velocity difference to the desired P-wave primary energy or with a different type of moveout (change of time of arrival with source-receiver distance) such as, for example, linear moveout.

Abstract: This invention relates in general to vibroseis and, more specifically, but not by way of limitation, to the enhancement and/or signal strength optimization of low frequency content of seismic signals for use in surveying boreholes and/or subsurface earth formations. In embodiments of the present invention, physical properties of a seismic vibrator may be analyzed and used to provide for determination of a driving force necessary to drive a reaction mass to produce a sweep signal with enhanced low frequency content for injection into the ground for vibroseis. In certain aspects, the physical properties may be considered independent of any geophysical properties related to operation of the seismic vibrator.

Abstract: A system and method is provided for characterizing earth formations. In one embodiment, the method includes passing a logging tool through a borehole and repeatedly: (a) triggering an acoustic wave generator; (b) recording acoustic waveforms received by receivers in the logging tool; (c) determining a time semblance of the recorded acoustic waveforms; and (d) smoothing the time semblance. In a different embodiment, a phase semblance of the recorded acoustic waveforms is determined and smoothed. The smoothing may be performed using an adaptive wavelet denoising technique or an adaptive moving average filter technique. In each case the average time or frequency spacing between semblance peaks is preferably determined and used to adapt the smoothing operation in a manner that varies with the slowness value s.

Abstract: A seismic signal is represented as a combination of a geologic signal and a noise signal. The representation of the seismic signal is decomposed into a linear combination of orthonormal components. One or more components are identified in which the geologic signal and the noise signal are uncorrelated. The noise signal is expanded as a product of an unknown noise amplitude modulation signal and a known noise spatial periodicity signal at the identified components. The expansion is solved for an estimate of the noise amplitude modulation signal at the identified components. The noise signal is estimated by multiplying the estimated noise amplitude modulation signal by the noise spatial periodicity signal. The estimated noise signal can be removed from the seismic signal to generate an estimate of the geologic signal.

Abstract: Near-offset and far-offset seismic data volumes are time-aligned by first selecting a plurality of time shifts. The near-offset and far-offset seismic data volumes are cross-correlated at the plurality of time shifts. An initial time-shift volume and a maximum correlation volume are created from the maximal cross-correlations at the plurality of time shifts. Areas of high time shift from the initial time-shift volume and areas of low cross-correlation from the maximum correlation volume are determined. The determined areas of high time shift and low cross-correlation are filtered from the initial time-shift volume, generating a filtered time-shift volume. The filtered time-shift volume is applied to the far-offset seismic volume to generate a time-aligned far-offset volume.

Abstract: This invention unifies a set of statistical signal processing, neuromorphic systems, and microelectronic implementation techniques for blind separation and recovery of mixed signals. A set of architectures, frameworks, algorithms, and devices for separating, discriminating, and recovering original signal sources by processing a set of received mixtures and functions of said signals are described. The adaptation inherent in the referenced architectures, frameworks, algorithms, and devices is based on processing of the received, measured, recorded or otherwise stored signals or functions thereof. There are multiple criteria that can be used alone or in conjunction with other criteria for achieving the separation and recovery of the original signal content from the signal mixtures. The composition adopts both discrete-time and continuous-time formulations with a view towards implementations in the digital as well as the analog domains of microelectronic circuits.

Abstract: A method relating to filtering coherent noise and interference from seismic data by constrained adaptive beamforming is described using a constraint design methodology which allows the imposition of an arbitrary predesigned quiescent response on the beamformer. The method also makes sure that the beamformer response in selected regions of the frequency-wavenumber space is entirely controlled by this quiescent response, hence ensuring signal preservation and robustness to perturbations. Built-in regularization brings an additional degree of robustness. Seismic signals with arbitrary spectral content in the frequency-wavenumber domain are preserved, while coherent noise and interference that is temporally and spatially nonstationary is adaptively filtered. The approach is applicable to attenuation of all types of coherent noise in seismic data including swell-noise, bulge-wave noise, ground-roll, air wave, seismic vessel and rig interference, etc. It is applicable to both linear or areal arrays.

Abstract: An error-suppression signal measurement system and method therefor is provided. The system transmits a test signal from a first probe, through a device under test, and into a second probe. The probes extract normalization signals from reference signals therein, exchange specific ones of the normalization signals, and combine the normalization signals with data signals derived from the test signal to form receiver signals. The probes propagate the receiver signals to a receiver, where the signals are gain-ranged, digitized, normalized, and compensated for phase-noise.

Abstract: A method is disclosed for attenuating ground roll in seismic data signals. Seismic data signals are obtained from nearby seismic sensors. A difference filter is created that incorporates an estimate of ground roll differential move-out between the seismic sensors. The difference filter is then applied to the seismic data signals to produce a ground roll attenuated seismic data signal. The method can be used with one-pass and multiple-pass filter operations and can be used in connection with pairs of seismic data signals or with groups of seismic data signals from an areal pattern of nearby seismic sensors. The seismic data signals may be appropriately pre-processed and factors in addition to the ground roll differential move-out can be incorporated when creating the difference filter. The inventive method can be implemented in the time domain or in the frequency domain.

Abstract: A method of acquiring and processing seismic data involving deploying a number of seismic sensors, actuating a seismic source, receiving seismic signals produced by the seismic source using the seismic sensors, calculating offset distances (44) between the seismic source and the seismic sensors, and (42) producing spatially filtered seismic data from the received seismic signals. The inventive method may be used effectively with receiver arrays having smaller areal footprints and fewer sensors per receiver station location than conventional seismic data acquisition systems. It also allows ground roll and random noise to be effectively attenuated as the seismic data is spatially resampled or prepared for spatial resampling.

Abstract: A method relating to filtering coherent noise and interference from seismic data by constrained adaptive beamforming is described using a constraint design methodology which allows the imposition of an arbitrary predesigned quiescent response on the beamformer. The method also makes sure that the beamformer response in selected regions of the frequency-wavenumber space is entirely controlled by this quiescent response, hence ensuring signal preservation and robustness to perturbations. Built-in regularization brings an additional degree of robustness. Seismic signals with arbitrary spectral content in the frequency-wavenumber domain are preserved, while coherent noise and interference that is temporarily and spatially nonstationary is adaptively filtered. The approach is applicable to attenuation of all types of coherent noise in seismic data including swell-noise, bulge-wave noise, ground-roll, air wave, seismic vessel and rig interference, etc. It is applicable to both linear and a real arrays.

Abstract: A method for predicting seismic events wherein measurable and calculable parameters relating to alterations in the shape of the geoid, such as centripetal and tidal gravitational effects, and gravitational anomalies related to a buildup of energy at any given point in the earth are incorporated with observations from at least two previous, subsequent seismic events, to calculate the energy buildup required to result in a future seismic event.

Abstract: A marine seismic signal is transformed from time domain into frequency domain and represented by matrix D. The marine data signal is truncated in time, transformed into the frequency domain and represented by matrix D.sub.T Eigenvalue decomposition D.sub.T=S.LAMBDA.S.sup.-1 of matrix D.sub.T is computed. Matrix product D S is computed and saved in memory. Conjugate transpose [D S]* is computed and saved in memory. Matrix product [D S]*[D S] is computed and saved in memory. Matrix inverse S.sup.-1 is computed and saved in memory. Conjugate transpose (S.sup.-1)* is computed and saved in memory. Matrix product S.sup.-1 (S.sup.-1)* is computed and saved in memory. An initial estimate of the source wavelet w is made. Source wavelet w is optimized by iterating the steps of computing the diagonal matrix [I-w.sup.-1 .LAMBDA.], computing matrix inverse [I-w.sup.-1 .LAMBDA.].sup.-1, computing conjugate transpose [(I-w.sup.-1 .LAMBDA.).sup.-1 ]*, retrieving matrix products [D S]*[D S] and S.sup.-1 (S.sup.

Abstract: A method and computer system for deriving and applying normal moveout corrections to seismic survey traces is disclosed. In each common midpoint (CMP) gather, a stacking velocity function of travel time is first determined for near-offset traces. For each reflective event, a time window about a depth point is selected. Time anisotropy and velocity anisotropy values are then determined for the depth point, for example by way of a semblance analysis applied to envelopes of the traces over the full offset range. Normal moveout corrections are then determined by a weighted sum expression of travel time (zero-offset time plus moveout correction) The weighted sum adds a first travel time expression corresponding to a Dix equation expression with a second travel time expression corresponding to a skewed hyperbolic expression that includes the time and velocity anisotropy values, with the weighting dependent upon the offset of the particular trace.

Abstract: A method for filtering noise from discrete noisy seismic signals in a time interval ([1, . . . , T]) is provided. The method comprising the steps of using the noisy seismic signals for determining at least one reference channel (x(t)) for the time interval as an estimate of the noise; determining coefficients for M temporally local filters (w(i, t)), the filters forming a filter bank, and M being a number equal to or larger than two, by minimizing a cost function (J(t)) representing a measure of the error of the output of the filter bank, and applying the filter bank to at least one estimate (x(t)) to determine M filtered estimates of the noise. The filtered estimates are multipled with temporal windows (h(i, t)).

Abstract: The present invention relates to a filtering method for discriminating elliptical waves among other waves propagating in a material medium (such as a geologic formation for example) by a combined processing of the components measured along several axes of the waves received by a multi-axis receiver. It comprises detection by at least one receiver (R) (consisting for example of one or several directional detectors such as geophones or accelerometers coupled with said formation) of the wave components in at least two orthogonal directions. A determination (by a processing assembly (1)) of the wavelet transforms of the signals is produced by this receiver (R) in response to the waves received and selection of the amplitudes of the signals resulting from this wavelet transformation is performed according to a criterion based on the ratio of the respective amplitudes thereof. The method can be used for the filtering of surface or tube waves within the scope of a seismic exploration.

Abstract: A method for attenuating source-generated surface wavetrains in a set of seismic data traces. According to the method, the seismic data traces are Fourier transformed with respect to time to determine the frequency components of the surface wavetrain. A phase matching operation is then performed using estimates of the wavenumbers for each of the frequency components to approximately align the surface wavetrain in the seismic data traces. Next, a spatial low-pass filter is used to remove the approximately aligned surface wavetrain. Finally, the phase matching is removed and an inverse Fourier transformation is performed to return the seismic data traces to the time domain.

Abstract: The invention relates to a seismic prospection method in which a seismic disturbance is generated in the sub-soil and sensors are used to pick up seismic data, the seismic data being processed to deduce useful information concerning the geology of the sub-soil, the seismic data containing a signal y.sub.0 that is to be isolated, said signal y.sub.0 being present in the seismic data embedded in an additive signal. The method comprises the steps consisting in: 1) for .alpha. given integer p, computing a prediction error filter a having p coefficients and providing best cancellation of the signal y.sub.0 ; 2) applying the self-deconvoluted prediction error filter to the seismic data to obtain filtered data in which the signal y.sub.0 is absent; and 3) subtracting the filtered data from the initial data to obtain processed data containing the signal y.sub.0 without the additive signal.

Abstract: A computer-implemented method and system of correcting seismic survey data for the effects of NMO stretch is disclosed. The method operates upon a CDP gather of seismic survey data after normal move-out correction (NMO) has been applied. A semblance analysis is used to derive a stretch coefficient .kappa. profile over the gather, as a function of time and offset. The values of the stretch coefficient profile .kappa. are used in producing a time-varying filter that is applied, preferably by way of a time-domain multiplication, to incremental windows in time for each trace, with the products added to one another to generate a corrected trace. After all of the traces in the gather are corrected, an event-weighted filter is applied to the traces by way of a correlation operation, to remove high-frequency artifacts in the corrected traces. The corrected and filtered traces are then ready for stacking and conventional processing.