Abstract: A geophysical model of a subsurface region is generated based on seismic data, e.g., seismic reflection data. Migration and seismic inversion are applied to the seismic data to generate estimates of one or more physical or seismic properties of the subsurface region. Seismic inversion, such as spectral shaping inversion, is applied before or after migrating the seismic data through a variety of techniques that each avoid the amplification of dipping energy while optimizing computational efficiency and/or accuracy.

Abstract: A migration method using hybrid one-way and full-way (HOF) wave equation propagation. The HOF method extrapolates seismic wavefields in less complex media with a one-way wave equation propagator and extrapolates seismic wavefields in extremely complex media with a full-way wave equation propagator. For prestack depth migration, the HOF extrapolates source-side and receiver-side wavefields independently. Frequency-space domain and time-space domain imaging conditions are applied to the one-way and the full-way extrapolated wavefields, respectively. A suitable amplitude matching factor is introduced to combine the one-way and full-way images. The HOF method is a cost-effective migration that produces superior image quality with less noises and less computational resources.

Abstract: Computing device and method for processing seismic traces to produce an image of a subsurface area. The method includes receiving a series of seismic traces related to the subsurface area and recorded by one or more seismic receivers, wherein the one or more seismic sources are originally generated by a source; applying a phase encoding function to the series of seismic traces, at least a portion of said seismic traces comprise signals reflected by geological interfaces of the subsurface area; applying a 3 dimensional (3D) harmonic-source reverse time migration of the series of seismic traces encoded with the phase encoding function; computing a forward wavefield by solving a first wave equation; computing a backward wavefield by solving a second wave equation; and cross-correlating the forward wavefield with the backward wavefield to generate an image of the subsurface.

Abstract: Images relating to a subsurface region may be generated based at least in part on a backward propagated source wavefield and a receiver wavefield. A source wavefield may be propagated from an initial wavefield-state forward in time, from an initial time-state to a final time-state, through an earth model associated with the subsurface region. The backward propagated source wavefield may be determined by propagating the source wavefield backward in time, from the final time-state to the initial time-state, through the earth model to reconstruct the initial wavefield-state. The receiver wavefield may be propagated, from the final time-state, through the earth model. The earth model may include at least one boundary region that can be defined as having one or more of absorbing characteristics, boosting characteristics, randomly perturbed characteristics, and/or other characteristics. As such, wavefields may be dampened, amplified, randomly scattered, and/or otherwise altered at the at least one boundary region.

Abstract: A technique includes monitoring acquisition activity of a plurality of seismic vibrators. The technique includes receiving signals from the seismic vibrators during the monitoring. Each of the signals indicates that at least one of the seismic vibrators is available for an associated seismic operation. The technique includes, in response to the signals, regulating the operations based on the monitored acquisition activity.

Abstract: Presented are systems and methods for wireless data acquisition. The wireless data acquisition may involve synchronizing modules within a data acquisition array. The synchronized data acquisition array may be used to facilitate a seismic survey. Synchronization may be facilitated by receipt of a reference time event such that a clock is synchronized based on the reference time event.

Abstract: A synchronization system for a nodal geophysical data recorder includes an electromagnetic transmitter associated with a master recording unit. The master unit includes devices for determining time and geodetic position from an external reference. The transmitter includes a code generator to cause transmission of a time synchronization signal as a coded sequence. The transmitter is configured to induce an electromagnetic field in at least one of subsurface rock formations and a body of water. At least one nodal geophysical data recorder includes at least one geophysical data sensor. The at least one sensor has measurements therefrom stored in a data storage device associated with the recorder wherein the recorder includes a clock for time indexing the stored data measurements. The recorder includes a receiver for detecting and decoding the time synchronization signal in the coded sequence to synchronize the clock with the synchronization signal.

Abstract: The invention is a method of calculating seismic time shifts ?{right arrow over (t)}(t), comprising: conducting a number (n) of three or more seismic surveys with large time differences such as months or years, producing a number (n) of three or more vintages of seismic data d1(t), d2(t), . . .

Abstract: Method for designing a converted mode (PS or SP) seismic survey to accomplish specified vertical and lateral resolution objectives at target depth. An equation (181) is provided for determining the minimum bandwidth required for a desired vertical resolution at a selected scattering angle, as a function of incident and reflected wave velocities, one of which is the P-wave velocity and the other is the S-wave velocity. A second equation (182) is provided for determining migration acceptance angle from the desired vertical and lateral resolutions. Source and receiver apertures may then be determined by ray tracing. Finally, a third equation (183) is provided for the maximum bin size to avoid aliasing, given the migration acceptance angle and a maximum frequency needed to achieve the bandwidth requirement. Source and receiver spacing may then be based on the maximum bin size.

Abstract: A wireless seismic data acquisition unit with a wireless receiver providing access to a common remote time reference shared by a plurality of wireless seismic data acquisition units in a seismic system. The receiver is capable of replicating local version of remote time epoch to which a seismic sensor analog-to-digital converter is synchronized. The receiver is capable of replicating local version of remote common time reference for the purpose of time stamping local node events. The receiver is capable of being placed in a low power, non-operational state over periods of time during which the seismic data acquisition unit continues to record seismic data, thus conserving unit battery power. The system implements a method to correct the local time clock based on intermittent access to the common remote time reference. The method corrects the local time clock via a voltage controlled oscillator to account for environmentally induced timing errors.

Abstract: Seismic systems and methods are provided to synchronize both source and receiver data using inexpensive timers and/or low energy timers to obtain high resolution seismic data.

Abstract: A method (which can be computer implemented) for addressing errors during migration of software applications includes the step of obtaining access to a data repository, which includes a listing of migration decisions and a listing of dependencies associated with the migration decisions. The method further includes the steps of obtaining an indication of a migration error, tracing the migration error to at least a first one of the migration decisions in the data repository, and employing the listing of dependencies to identify at least a second one of the migration decisions, depending on the at least first one of the migration decisions identified in the tracing step, which is impacted by the migration error. Techniques for populating the data repository are also provided.

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

Abstract: A method and system for generating images of a subsurface region of interest. In general, one embodiment of the present invention includes establishing boundary conditions utilizing seismic data and initial conditions which include excitation from source locations in an earth model. Source wavefields are then propagated forward through the earth model to a maximum time, and saved at a plurality of checkpoints sparsely in time and also corresponding boundary values of the source wavefields at each time step are saved. Source wavefields are also propagated backward through the earth model from the maximum time utilizing the plurality of checkpoints when available and the saved boundary values at each time step. Receiver wavefields are propagated backward concurrently through the earth model from the maximum time.

Abstract: 1st order free-surface multiples recorded in VSP data or reverse VSP data are processed using a 3-C 3-D vector migration method to produce an image of the subsurface. This image produces a larger coverage than that obtained in 3-C 3-D processing of reflection data acquired in the VSP.

Abstract: A method for estimating seismic data from sources of noise in the earth. The method includes calculating a first set of Green's functions between a first seismic receiver and each seismic receiver of an array of seismic receivers using interferometry. The array of seismic receivers is disposed around the first seismic receiver, and the first seismic receiver is part of the array of seismic receivers. The method also includes calculating one or more correction factors to correct the first set of Green's functions and calculating a second set of Green's functions between each seismic receiver of the array and a second seismic receiver using interferometry. Here, the first seismic receiver is disposed at a predetermined distance away from the second seismic receiver. Also, the first set and the second set of Green's functions are biased due to non-identical strengths of the sources of noise.

Abstract: A method, apparatus and system for acquiring land seismic data includes acquiring seismic data with a first autonomous seismic data acquisition unit and a second autonomous seismic data acquisition unit wherein each acquisition unit comprises a plurality of digitally controlled temperature-compensated crystal oscillators. Oscillator-based timing signals are acquired that are associated with the plurality of digitally controlled temperature compensated crystal oscillators and a time correction is determined to apply to the seismic data acquired with the first autonomous seismic data acquisition unit. The time correction is determined using the oscillator-based timing signals from the first and second autonomous seismic data acquisition units.

Abstract: A method, apparatus and system for acquiring land seismic data includes acquiring seismic data with a first autonomous seismic data acquisition unit and a second autonomous seismic data acquisition unit wherein each acquisition unit comprises a plurality of digitally controlled temperature-compensated crystal oscillators. Oscillator-based timing signals are acquired that are associated with the plurality of digitally controlled temperature-compensated crystal oscillators and a time correction is determined to apply to the seismic data acquired with the first autonomous seismic data acquisition unit. The time correction is determined using the oscillator-based timing signals from the first and second autonomous seismic data acquisition units.

Abstract: Implementations of various technologies for a method for processing seismic data. In one implementation, the method includes (a) selecting a first trace from a first seismic data set and a second trace from a second seismic data set; (b) extracting one or more features of the same types from the first trace and the second trace; (c) matching the extracted features from the first trace with the extracted features from the second trace; and (d) calculating for a displacement field using one or more of the matching features of the first trace and the second trace.

Abstract: To process seismic data, a combined four-dimensional (4D) binning and regularization procedure is performed on the seismic data, where the combined 4D binning and regularization procedure includes computing measures associated with regularization of the seismic data, computing measures associated with 4D binning, and processing the seismic data according to the regularization measures and 4D binning measures.

Abstract: A method for interpreting seismic data below a salt layer includes depth migrating the seismic data to a bottom of the salt layer. The migrating including generating an initial model of velocities below the salt layer. The initial model is scaled by a plurality of scale factors at at least one image position. At least one of the plurality of scale factors for which a depth migrated image below the salt layer is optimum is selected as the scale factor.

Abstract: A method is disclosed for reflection time shift matching a first and a second set of seismic reflection data (10,30) comprising first and second seismic reflection traces (1,3) with series of generally corresponding seismic reflections (11,31). The second set of seismic data (30) comprises at least one laterally extending series (40) of new seismic events (4) not present in the first set of seismic data (10). Reflection time shifts (22) are calculated as required for matching seismic reflections (31) of the second reflection traces (3) of the second seismic reflection data (30) with corresponding seismic reflections (11) of the first reflection traces (1). The calculation of the time shifts (22) are made by calculating coefficients of basis function estimates such as Legendre polynomials in order to save calculation effort.

Abstract: A network distributed seismic data acquisition system comprises seismic receivers connected to remote acquisition modules, receiver lines, line tap units, base lines, a central recording system and a seismic source event generation unit synchronized to a master clock. One or more high precision clocks is added to the network to correct for timing uncertainty associated with propagation of commands through the network. Seismic receivers and seismic sources are thereby synchronized with greater accuracy than otherwise possible. Timing errors that interfere with the processing of the seismic recordings are greatly reduced, thus enabling an improvement in subsurface geologic imaging.

Abstract: A method for spectrally conditioning surface seismic data. In one implementation, the method may include correcting surface seismic data for distortions due to anomalous spectral amplitudes, thereby generating a first set of corrected data; correcting the first set of corrected data for deterministic distortions, thereby generating a second set of corrected data; correcting the second set of corrected data for spectral distortions due to the seismic waves traveling through the near-surface, thereby generating a third set of corrected data; and correcting the third set of corrected data for spectral distortions due to the seismic waves traveling through deeper strata.

Abstract: A method is used to convert seismic data between bounding upper and lower surfaces in the time domain to the depth domain. An upper depth surface corresponding to the upper time surface is determined, either from well data or by time-depth conversion of the upper time surface. An interval velocity value is determined for each seismic sample between the bounding time surfaces from velocity data generated from well logs. The seismic times between the bounding time surfaces are then converted to the depth domain by adding the thickness of each seismic sample to the upper depth surface, wherein the thickness of each seismic sample is a function of the seismic sample rate and interval velocity value for each seismic sample. The seismic attribute values from the converted seismic data can be interpolated to a linear depth sampling.

Abstract: The technologies described herein include systems and methods for encoding/decoding seismic sources and responses, generating and using of source-side derivatives while also generating and using the conventional source response. Sources in an array may be encoded such that activation of each source in the array constitutes a single spike in a sequence orthogonal to another sequence emitted by another source. The responses to these different sources that are in close spatial proximity can be decoded and separated. Source-side derivatives may be calculated and utilized in various applications in combination with the monopole response from the source location, including source-side deghosting, spatial (horizontal and vertical) interpolation and imaging.

Abstract: 1storder free-surface multiples recorded in VSP data or reverse VSP data are processed using a 3-C 3-D vector migration method to produce an image of the subsurface. This image produces a larger coverage than that obtained in 3-C 3-D processing of reflection data acquired in the VSP.

Abstract: A method and apparatus for predicting a plurality of surface multiples for a plurality of traces in a record of seismic data. In one embodiment, the method includes providing a plurality of target traces at a nominal offset and a nominal azimuth; selecting a plurality of pairs of input traces, wherein the midpoints of the input traces in each pair are separated by half the nominal offset and the azimuth of a line connecting the midpoints of the input traces in each pair is equal to the nominal azimuth; convolving the selected pairs of input traces to generate a plurality of convolutions; and applying a three dimensional operator to the convolutions.

Abstract: Earth analysis methods, subsurface feature detection methods, earth analysis devices, and articles of manufacture are described. According to one embodiment, an earth analysis method includes engaging a device with the earth, analyzing the earth in a single substantially lineal direction using the device during the engaging, and providing information regarding a subsurface feature of the earth using the analysis.

Abstract: Methods of processing seismic data to remove unwanted noise from meaningful reflection signals are provided for. The methods comprise the steps of assembling seismic data into common geometry gathers in an offset-time domain without correcting the data for normal moveout. The amplitude data are then transformed from the offset-time domain to the time-slowness domain using a limited Radon transformation. That is, the Radon transformation is applied within defined slowness limits pmin and pmax, where pmin is a predetermined minimum slowness and pmax is a predetermined maximum slowness. A corrective filter is then applied to the transformed data enhance the primary reflection signal content of the data and to eliminate unwanted noise events. After filtering, the enhanced signal content is inverse transformed from the time-slowness domain back to the offset-time domain using an inverse Radon transformation.

Abstract: A method and system for seismic data processing utilizes azimuthal variations in the velocity of seismic signals. The system and method utilizes a plurality of seismic energy sources that are located at known positions at the surface of the earth. The seismic energy sources generate seismic signals that propagate downward into the earth. Some of the seismic signals are reflected and diffracted by various sub-surface layers and are returned to the surface of the earth. The returned seismic signals are received by a plurality of receivers. The method includes the step of determining the distance from an energy source to an image point. A fast travel time of the seismic signal from the energy source to the image point is determined, and a slow travel time of the seismic signal from the energy source to the image point is determined. The azimuth angle between the energy source and the surface location of the image point is calculated.

Abstract: A method is disclosed for migrating seismic data which includes determining travel time of a compressional wave from a source location to a scatter point, taking into account ray bending. Travel time of a shear wave from the scatter point to a receiver location is determined, taking into account ray bending at the interfaces between subsurface strata. The determined travel times are then used to migrate the seismic data. In one embodiment, the travel times take account of vertically transversely isotropic media with a vertical symmetry axis.

Abstract: A method of seismic processing for obtaining information on the geophysics of the subsoil includes acquiring seismic traces at at least one point on the surface of the subsoil or in the subsoil, the seismic traces corresponding on each occasion to two perpendicular components of a shear wave emitted into the subsoil and reflected by different interfaces therein; applying a succession of transformations (?) at least to a temporal portion of the traces for each of these assumptions, determining the value of a parameter representative of the coherence/similarity between the result traces obtained in this way; and election as a function of the values obtained in this way that one of the hypotheses which is considered as being most representative of the subsoil, the two result traces obtained for said hypothesis being compensated traces of the subsoil birefringency.

Abstract: Time slices of seismic data are transformed from rectangular space-time domain to cylindrical space-time domain. 2-D seismic migration is performed on the transformed data for each radial direction. Slices of the migrated data are inverse transformed back to the rectangular space-time domain, generated migrated 3-D data for generally inhomogeneous media.

Abstract: A method for improving the seismic resolution in the field of petroleum exploration and development, using the steps of: (1) ensuring the depth and the velocity of the low-velocity layer and the high-velocity layer and the characteristic of the frequency spectrum of the seismic wave at the near-surface using the borehole seismic survey; (2) separating the ghost wave from the seismic record in the high-velocity layer using the processing method of the vertical seismic profiling, and improving the SNR of the seismic record of the direct wave at the high-velocity layer using the uphole stacking of the vertical seismic profiling; (3) evaluating the de-convolution operator of the near-surface by the convolution math-physics equation; (4) evaluating the seismic record with high resolution by the de-convolution math-physics equation; (5) providing the band-pass filtering and geologic interpretation to the seismic record with high resolution.

Abstract: A one-return wave equation migration is used to extrapolate both down-going and up-going waves. Followed by a properly designed imaging condition, the partial image contributed form turning waves is correctly reconstructed. Numerical examples show that this method can significantly enhance definition of an overhanging salt boundary and a geological structure with vertical features.

Abstract: An apparatus and method for processing a record of seismic traces. The method includes receiving the record of seismic traces in a time-offset domain, transforming the record of seismic traces to a tau-p domain, applying a zero offset inverse Q algorithm to the record of seismic traces in the tau-p domain, and transforming the record of seismic traces back to the time-offset domain.

Abstract: A technique includes providing a source in a first well and a seismic receiver in a second well to detect a seismic event that is caused by the source. The technique includes referencing clocks in the source and receiver to a common reference time frame and determining a time in the reference time frame at which the seismic source generates the seismic event.

Abstract: A method is disclosed for depth migrating seismic data. The method includes pre-stack time migrating the seismic data to form a stacked, time migrated image. The stacked, time migrated image is demigrated, and post-stack depth migration is then performed on the demigrated image. In some embodiments, the pre-stack time migration and the demigration account for ray bending in vertically transversely isotropic media.

Abstract: Disclosed is a method for locating a seismic event. The method includes processing seismic data from at least one seismic receiver to validate a potential seismic event, computing a signal travel time between at least one node in an area of interest and the at least one seismic receiver, adjusting the seismic data according to the travel time, and identifying a location of the seismic event based on the adjusted seismic data. Systems for locating a seismic event are also disclosed.

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

Abstract: System and method for performing anisotropy corrections of post-imaging seismic data for a subsurface formation. The method may receive seismic data, preferably pre-stack seismic data comprising a plurality of traces, e.g., collected from a plurality of source and receiver locations. The seismic data may be imaged/migrated to produce imaged seismic data, which is organized into an arrangement that preserves aspects of the relative seismic propagation angle in the subsurface. One or more anisotropic parameters and corresponding corrections may be determined by analyzing the organized imaged seismic data. The determined parameters or corrections may be used to correct at least a subset of the imaged seismic data, thereby producing corrected seismic data which is useable in analyzing the formation. The corrected data may then optionally be stacked to produce a collection of corrected stacked traces, and/or analyzed as desired. The pre-stack data and/or the corrected seismic data may optionally be displayed.

Abstract: Method of determining the specular source/receiver pairs associated with seismic images produced from a transformation of multi-offset seismic wave records as a function of time in the depth domain. When using a prestack depth migration algorithm of Kirchhoff type, one also calculates and saves, in addition to the necessary traveltimes of the rays from the sources and receivers of the acquisition device to the image points of the image to be migrated, the line parameters of these rays and the desired quantities to be reached in fine. Then, once the prestack migrated volume produced, the local slopes of seismic events of interest for the image points considered are extracted. Comparison between line parameters saved during migration and the normal to the local slopes at the image points considered shows which ray is specular among those for which information has been saved, and allows to produce then the required quantity, which is then referred to as specular quantity.

Abstract: The present invention provides a method for analyzing traces collected along a plurality of adjacent sail lines in a marine seismic survey area. The method comprises selecting a plurality of trace groups, at least one trace in each group sharing a common midpoint with at least one trace from another group, determining an initial zero-offset travel time for each trace in the trace groups, and generating a plurality of updated zero-offset travel times and a plurality of time corrections for the trace groups using a pre-selected function of the initial zero-offset travel times.

Abstract: A method is used to convert seismic data between bounding upper and lower surfaces in the time domain to the depth domain. An upper depth surface corresponding to the upper time surface is determined, either from well data or by time-depth conversion of the upper time surface. An interval velocity value is determined for each seismic sample between the bounding time surfaces from velocity data generated from well logs. The seismic times between the bounding time surfaces are then converted to the depth domain by adding the thickness of each seismic sample to the upper depth surface, wherein the thickness of each seismic sample is a function of the seismic sample rate and interval velocity value for each seismic sample. The seismic attribute values from the converted seismic data can be interpolated to a linear depth sampling.

Abstract: A system and method for pre-stack imaging system that includes multiple processors, with each processor configured to perform two-way wave equation-based forward wave propagation of a source event to obtain a source wave field using high order spatial finite differences for a shot gather, and to perform two-way wave equation-based backward wave propagation of received waves to obtain a receiver wave field using high order spatial finite differences for the shot gather; and to correlate the source and receiver wave fields to obtain an image region.

Abstract: A method is disclosed for amplitude preserving Kirchhoff time migration. The method includes calculating estimates of the geometrical spreading terms using a constant velocity approximation. Takeoff and emergence angles are accurately calculated based on a model of velocity which varies with respect to travel time. An image can be calculated at at least one travel time using the estimated weights and calculated takeoff and emergence angles.

Abstract: A filter that includes an analyzer, thresholding circuit, and synthesizer. The analyzer generates a low-frequency component signal and a high-frequency component signal from an input signal. The thresholding circuit generates a processed high-frequency signal from the high-frequency component signal, the processed high-frequency signal having an amplitude of zero in those regions in which the high-frequency component signal has an amplitude that is less than a threshold value. The synthesizer generates a filtered signal from input signals that include the low-frequency component signal and the processed high-frequency signal. The filtered signal is identical to the input signal if the threshold value is zero. The analyzer is preferably constructed from a plurality of finite impulse response filters that operate on a small fraction of the input signal at a time.

Abstract: A method is disclosed for seismic imaging of subsurface diffractors. The method includes performing migration velocity analysis on a seismic time record section and depth migrating the time section for offsets exceeding one-half a distance between a seismic energy source and a seismic receiver most distant from the source during acquisition of seismic data used to generate the time record section.

Abstract: A method is disclosed for processing seismic data. The method includes prestack depth migrating the seismic data to generate common image gathers using an initial velocity-depth model. Horizons in the migrated seismic data are selected. Residual migration velocity analysis in the depth-offset domain is performed with respect to each selected horizon, and the velocity-depth model is updated based on the residual migration velocity analysis.