Abstract: A method for performing pre-stack time migration of seismic data in a region where the subsurface seismic wave propagation velocity varies in both the vertical direction and a horizontal direction. The migration is performed in the wave number—frequency domain and involves using travel time maps to find a linear fitting of the migration phase shift for a given ray parameter and at a given vertical depth versus horizontal position. The slope and intercept parameters from the linear fit are used to adjust a known pre-stack time migration equation, with the accuracy of the linear approximation requiring mild horizontal velocity variation.

Abstract: A method of analyzing dip in a seismic data volume in which a horizontal gradient is calculated in a first direction in the seismic data volume. A vertical gradient is calculated at data locations in the seismic data volume corresponding to the locations at which the horizontal gradient was calculated. Dip is calculated in the first direction from the horizontal gradient in the first direction and the vertical gradient. Repetition of the process for the entire seismic data volume results in a dip volume.

Abstract: The invention is a method for analyzing spatially-varying noise in seismic data. Transitions between data values at adjacent data locations in a seismic data set are represented by Markov chains. Transition probability matrices are constructed from the Markov chains. Data values are predicted from the calculated transition probabilities. Noise values are determined from the predicted data values.

Abstract: A method for predicting the connectivity of seismic objects identified by seed detection or similar methods.

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: A method of calculating a throw volume corresponding to a seismic data volume. A range of time shifts and a search direction for the seismic data volume are selected. A data location separation and a vertical time window are also selected. A cross-correlation is calculated between data values corresponding to first and second data locations separated by the data location separation and symmetrically located in the search direction on each side of a target data location. The cross-correlation is calculated throughout the vertical time window for each time shift in the range of time shifts. The time shift corresponding to the maximum calculated cross-correlation is stored in the throw volume.

Abstract: A method for assessing the suitability of seismic data for quantitative amplitude analysis, where the concern is excessive residual normal moveout (RNMO). The invention uses a near offset stack and a far offset stack, the time difference between the two, a mute pattern, a reflection shape assumption for the RNMO, and a waveform and frequency for the far stack traces to generate a formula that estimates far stack amplitude error caused by RNMO. The formula can be used to compensate the far stack amplitude where the error is not so great as to require reprocessing of the data. The method can also be applied to interpreted amplitude maps.

Abstract: Seismic horizons are transferred between two three-dimensional seismic data volumes which cover the same area. The seismic data volumes are time-aligned to generate a time-shift volume. A seismic horizon is selected and the time-shift volume is applied to the seismic horizon, generating a time-shifted seismic horizon.

Abstract: A method for seismic exploration using nonlinear conversions between electromagnetic and seismic energy, with particular attention to the electromagnetic source waveform used. According to the invention, seismic returns from a source waveform are correlated with a reference waveform, with both waveforms custom designed to minimize both correlation side lobes and interference from linear electroseismic effects. A waveform element is selected which may be sequenced by a binary or similar digital code embodying the desired custom design to generate an input sweep with the needed depth penetration and noise suppression. Correlation of the seismic response with the reference waveform in a data processing step mathematically aggregates the seismic response from the input sweep into a single wavelet. Preferred binary digital codes include prescribed variations of maximal length shift-register sequences. Also, an apparatus for generating the desired waveforms.

Abstract: A method for predicting hydrocarbon-bearing zones and estimating rock properties by analyzing fluids trapped in the pore spaces or adsorbed on the surfaces of rock samples. The trapped gases are removed under vacuum and analyzed by a mass spectrometer. Data peaks corresponding to petroleum constituent molecules provide an indication of presence and abundance of hydrocarbons. A decrease of the count rate over time is used to estimate permeability and other rock properties. Concentration ratios for selected constituents indicate oil quality and depth of the oil-water contact.

Abstract: AVO anomalies are classified in near-offset and far-offset seismic data volumes, by first calculating a plurality of initial AVO seismic attributes representative of the offset seismic data volumes. A probabilistic neural network is constructed from the calculated initial AVO seismic attributes. AVO anomaly classifications are calculated in a portion of the offset seismic data volumes. The preceding steps are repeated until the calculated AVO anomaly classifications in the portion of the offset seismic data volumes are satisfactory. AVO anomaly classifications are calculated throughout the offset seismic data volumes using the constructed probabilistic neural network.

Abstract: An apparatus and method for suppressing multiples in the collection of marine seismic data comprises at least one source positioned in the body of water; at least one receiver positioned in the body of water below the air-water interface and near the sources; a bubble diffuser positioned in the body of water so that the bubbles emitted by the bubble diffuser are positioned between the receivers and the air-water interface, wherein the emitted bubbles provide high acoustic reflection and substantially suppress specular reflection of seismic waves; and a control for activating the bubble diffuser during the collection of seismic data by the sources and receivers.

Abstract: A method for surface estimation of reservoir properties, wherein location of and average earth resistivities above, below, and horizontally adjacent to the subsurface geologic formation are first determined using geological and geophysical data in the vicinity of the subsurface geologic formation. Then dimensions and probing frequency for an electromagnetic source are determined to substantially maximize transmitted vertical and horizontal electric currents at the subsurface geologic formation, using the location and the average earth resistivities. Next, the electromagnetic source is activated at or near surface, approximately centered above the subsurface geologic formation and a plurality of components of electromagnetic response is measured with a receiver array. Geometrical and electrical parameter constraints are determined, using the geological and geophysical data.

Abstract: A method for training a probabilistic neural network to map seismic attributes or similar quantities.

Abstract: A method for moveout-correcting CDP gathers of seismic data that takes nonhyperbolic moveout at large offsets into account. The entire range of offsets is divided into smaller ranges. Within each smaller range, a curve is fitted to the reflection time vs. offset data, and these separate curves are made to be continuous at range boundaries. The curves may be the hyperbolic or quartic options available in commercial seismic data processing packages. Additional computation time required due to lack of a single velocity function can be compensated for by interpolation techniques. The method has application to, inter alia, AVO analysis.

Abstract: A method for constructing velocity models for stacking seismic data, or, more generally, a method for determining reflection geometries with seismic gathers through the use of trace-to-trace coherency coupled with global editing capability to separate velocity events corresponding to primary reflections from multiple reflections and other noise. A velocity model can be fit to the edited velocity events and used for stacking of seismic data. Alternatively, the measured reflection geometry has other potential uses including providing data to be input to a tomographic velocity updating procedure that produces a migration velocity model.

Abstract: A method for seismic exploration using conversions between electromagnetic and seismic energy, with particular attention to the electromagnetic source waveform used. According to the invention, source waveforms are correlated with reference waveforms selected to minimize correlation side lobes. Line power at 60 Hz may be used to provide a waveform element which may be sequenced by a binary code to generate an extended source waveform segment with minimal correlation side lobes. Preferred binary codes include Golay complementary pairs and maximal length shift-register sequences.

Abstract: A method for correcting seismic amplitudes for geometrical artifacts resulting in non-uniform illumination of reflector surfaces. The illumination is calculated for each offset as a function of spatial position. The amplitudes are scaled by a scale factor which may be the reciprocal of the illumination. Data considered to be poor quality due to low illumination may be discarded, and the amplitudes scaled further to compensate for the discarded data.

Abstract: A method for compensating for amplitude loss in a set of seismic traces, where the amplitude loss is caused by shallow attenuation anomalies. Attenuation factors are calculated from first-arrival amplitudes for long-offset seismic data traces. Those factors are calculated by numerical methods, using matrix inversion and iteration. Those factors are used to adjust upward the amplitudes of all later-arriving seismic traces. Frequency dependent attenuation is treated by dividing the seismic traces into frequency ranges using band-pass filters, then using the above-described method within each frequency range after which the adjusted amplitudes are recombined. Amplitude variations caused by inconsistencies in the data acquisition system, such as a source strength variation, are also treated by a similar approach.

Abstract: A method for calculating seismic velocity for migration purposes as a function of subsurface spatial position that gives the seismic processing analyst direct control of the resulting migration velocity model, and includes a means to ensure that the model is consistent with RMS velocity sweeps computed from available surface seismic or other data. The method involves generating average velocities from the model to compare to the measured RMS velocity data, and comparing the model predictions to actual data on a 3-D interactive visual display allowing quick and easy model adjustment. Other available measured data besides surface seismic data maybe used to develop the model.