Abstract: Seismic image filtering machines, systems, program products, and computer implemented methods are provided to generate a filtered seismic image responsive to filtered seismic image data generated by attenuating coherent seismic noise from surface waves of an unfiltered wavefield constructed from unfiltered seismic image data through a single downward extrapolation of the unfiltered wavefield using a plurality of nonstationary convolution operators to perform localized filtering at each of a plurality of spatial locations of the unfiltered wavefield. Various embodiments, for example, can beneficially handle strong lateral velocity variations thus making various embodiments effective tools to remove complicated coherent seismic noise which is typically in the form of exponentially decaying evanescent waves.

Abstract: Seismic data acquired by independent simultaneous sweeping (ISS®) techniques are processed is to attenuate random uncompressed cross-talk signals and improve the resolution of the pre-stack migrated time image. A frequency-varying mean filter is applied on cross-spread offset-azimuth gathers of the data. The frequency-space domain filter may vary its window size according to the characteristics of the cross-talk.

Abstract: Acquisition of data by managing crosstalk interference with sector designs and unique sweeps is conducted and the resultant data are processed in 3D common receiver domain to attenuate crosstalk noise while preserving the signals for high source and receiver density acquisition designs. High-amplitude spectral amplitudes are attenuated and inter-ensemble statics or structural time delays are applied to achieve optimum filter performance. If the spectral amplitudes have been attenuated to a level consistent with non-simultaneous acquisition, conventional surface consistent processing can be performed to correct for statics and amplitude variations. A 3-point filter in different frequency bands may then be applied to remove any remaining residual crosstalk noise.

Abstract: Disclosed are techniques for acquiring, processing and presenting data produced by one or more devices and meters in the field. One or more servers, such as Supervisory Control And Data Acquisition (SCADA) servers, gather and store data from one or more field devices. Remote polling interfaces (RPI) provide a secure transport from the SCADA servers to processing and storage elements of the present embodiments. Data gathered from a plurality of heterogeneous SCADA servers are normalized and consolidated into one or more databases for subsequent processing and presentation, thereby collecting and unifying data from a plurality of heterogeneous sources into a device-agnostic and unified presentation for user viewing and interaction.

Abstract: Beamed data can be obtained from shot gather data, and debeamed data can be obtained from beamed data. Shot gather data for a geophysical volume of interest in a first domain and/or beamed data in a third domain may be received. The first domain has a component relating to seismic wave propagation time and a spatial component relating to lateral spacing. The third domain has a component relating to local plane wave arrival time at a beam center surface location and a component relating to plane wave arrival directions. Data may be transformed between the first domain and a second domain. The second domain has a frequency component and a spatial component corresponding to the spatial component of the first domain. A forward or reverse transform may be applied to transform shot gather data in the first domain to/from beamed data in the third domain.

Abstract: Seismic survey data is converted into a subsurface data volume that users can view to understand and analyze the subsurface structures in the survey region. At least some of the disclosed systems and methods employ three-dimensional reverse time migration using wave equations adapted for use in formations having tilted transverse isotropy. Relative to existing methods, the disclosed systems and methods rely on fewer approximations and suffer fewer limitations on the circumstances in which they can be employed. Moreover, because the disclosed wave equations are derived from Hooke's law (and consequently they operate on fields that are directly tied to physical quantities), they demonstrate an increased stability relative to existing methods. Survey data analysts employing the disclosed systems and methods should obtain better images of the subsurface and be better able to identify reservoirs and deposits for commercial exploitation.

Abstract: The instant invention is a method for numerically propagating waves or solving wave equations on a digital computer. This invention can be used to compute highly accurate solutions to the wave equation, in fact in some cases it computes the analytical solution, something previously considered impossible. The instant method can also propagate waves that are not described by differential equations, such as anisotropic scalar waves. The invention has the advantages that it is computationally efficient, accurate, and flexible. Of importance is the ability to propagate waves that simulate the P-wave arrivals in both isotropic and anisotropic media with a scalar as opposed to a vector equation.

Abstract: A method for obtaining a volumetric image of a geological structure in the subsoil includes generating a plurality of seismic waves by plural seismic sources, acquiring signals or traces obtained by processing the seismic waves reflected by the geological structure for creating a seismic record or shot gather for the activation position of each seismic source, migrating, through a numerical processing, each shot gather to create a respective partial volumetric image of the geological structure, summing the single partial volumetric images to obtain the overall volumetric image of the geological structure, the summing of the single partial volumetric images being obtained by grouping the seismic sources into groups, grouping the groups of seismic sources into super-groups, and effecting N migrations for each super-group.

Abstract: A method is disclosed for automatically extending interpreter horizon picks over a wider area of traces in such manner that the automatically generated picks are very similar to picks that an interpreter would pick manually. The method applies optical filters to seismic sections to determine the intrinsic orientation of seismic events. Seismic orientation is captured in the Orientation Vector Field, which is then used to guide the picking process.

Abstract: To perform a marine survey, a streamer is deployed into a body of water, where the streamer has plural first sensors to perform a subterranean survey. Indications are received from second sensors in corresponding sections of the streamer, where the indications are regarding which sections are in the body of water. Information is updated regarding which sections of the streamer are in the body of water in response to the received indications.

Abstract: Methods and computing systems are disclosed to enhance imaging of acquired data. In one embodiment, a method is performed that includes receiving acquired data that corresponds to the medium; computing a first wavefield by injecting a noise; and computing the cumulative illumination by auto-correlating the first wavefield.

Abstract: The invention is a method of constructing an image representative of a heterogeneous medium by a procedure of joint inversion of seismic data represented on different time scales having an application for underground reservoir exploration. Sequential inversion of the seismic data is performed so as to determine a first model from the seismic data expressed in a first time scale, and a second model from the seismic data expressed in the second time scale. A scale factor t1(t0) allowing synthetic data described in a first time scale to be expressed in a second time scale is defined by a differential equation relating the traveltime variations of a seismic wave to the models resulting from the sequential inversion. Finally, joint inversion wherein a cost function using the scale factor is minimized is carried out so as to evaluate a difference between the synthetic data and the seismic data expressed in another time scale.

Abstract: Computing device, computer instructions and method for denoising seismic data recorded with seismic receivers. The method includes receiving the seismic data recorded with the seismic receivers, wherein the seismic data is recorded in a time-space domain; applying with a computing device a high resolution transform to the seismic data in the time-space domain to obtain transformed seismic data in a different domain than the time-space domain, such that the method is amplitude preserving; determining regions of noise and regions of true signal in the transformed seismic data; scaling down the regions of noise; and reverse transforming the transformed seismic data to the time-space domain.

Abstract: A method of processing stratigraphic data comprising a plurality of stratigraphic features, such as horizon surfaces, within a geological volume is provided. The method includes the steps of: extending a plurality of spaced sampling traces through the volume to traverse the stratigraphic features; and assigning the stratigraphic features respective relative geological ages. On each sampling trace, the relative geological age of each stratigraphic feature traversed by the sampling trace in the direction from geologically younger to geologically older stratigraphic features is increased in relation to the relative geological age of its preceding stratigraphic feature, under the condition that each stratigraphic feature takes the same relative geological age across all the sampling traces by which it is traversed.

Abstract: A method for fault and fracture identification based on seismic data representing a geological section using dispersion properties of reflected and diffracted waves. The method includes scanning N parameters associated with the seismic data. The array includes the coordinate axes of the angle of emergence (?), the radius of curvature of the wave front (R) and either time or depth samples. The method also includes processing the N parameters, generating a new image having a cross-sectional shape associated with one of the reflected and diffracted waves, calculating parameters DS and LS, evaluating DS for the case of fracture characterization, and comparing, for the case of fault identification, parameter LS with a threshold value defining the type of wave as one of the reflected and diffracted wave, the cross-sectional shape being substantially circular for the reflected wave, and being elliptical for the diffracted wave.

Abstract: A computer-implemented method for processing data includes accepting a first collection of traces corresponding to signals received over time due to reflection of seismic waves from subsurface structures. A Radon transform is defined with respect to a set of wavefront parameters of the seismic waves. The transform defines a summation of amplitudes of the seismic waves over trajectories defined by the wavefront parameters. The Radon transform is applied to the first collection of traces, so as to convert the first collection into a multidimensional data array that is defined as a function of at least two of the wavefront parameters. The multidimensional data array is processed to produce a second collection of traces having an improved imaging quality with respect to the first collection. The second collection of traces is processed to generate a seismic image of the subsurface structures at the improved imaging quality.

Abstract: A method for stratigraphic analysis of seismic data, known as GeoLeaf, comprises: selecting a seismic data volume comprising a plurality of samples; selecting a number of horizons (TH, BH1 TPH, BPH, SMH) within the seismic data volume; selecting a scale of interest within the seismic data volume, which ranges from 500 m up to the size of the entire seismic data volume; applying an attribute function to at least one of the selected number of horizons (TH, BH, TPH, BPH, SMH) at the selected scale of interest; and displaying a stratigraphic attribute, which represents a geologic feature, such as a seismic fades attribute (FM), an unconformity attribute (Zmax/Z) and/or a thinning attribute (Z1/Z2).

Abstract: A visualization method for the computation and display of seismic attributes from time-frequency spectral analysis enables high resolution investigation of 3D seismic data for the exploration of oil and gas. First, high frequency resolution time derivative or space derivative of the amplitude and phase spectra of time-windowed signals are computed. Second, the derivatives are transformed into seismic attributes that are suitable for seismic data processing or interpretation. Applications for the method are illustrated with graphical screen 3D volume displays of dip and azimuth, curvature and faults. Aside from providing high resolution these displays may improve the productivity of a seismic interpreter or data processor.

Abstract: A method and visualization apparatus for spectral analysis of time-and-space varying signals enables high resolution investigation of 3D seismic data for the exploration of oil and gas. The method extrapolates multi-resolution short windows into an average long window then computes its FFT. Extrapolation uses the continuity relationship between data inside and outside of short windows. Applications of the method are illustrated with graphical screen 3D volume displays of amplitude spectra, dip and azimuth, curvature and faults (figure below). Aside from high resolution these displays improve the productivity of a seismic interpreter.

Abstract: A method and system for processing synchronous array seismic data includes acquiring synchronous seismic data from a plurality of sensors to obtain synchronized array measurements. A reverse-time data propagation process is applied to the synchronized array measurements to obtain dynamic particle parameters associated with subsurface locations. A maximum energy density imaging condition is applied to the dynamic particle parameters to obtain imaging values associated with subsurface locations. Subsurface positions of energy sources are located from the relative maximum of a plurality of the imaging values associated with subsurface locations.

Abstract: An apparatus and method for mapping a water pathway. A group of sensors can be employed for detecting one or more features associated with the water pathway in a direction of flow through the water pathway. A buoyant vessel maintains the sensors, and the sensors assist in compiling data indicative of the detected features. The velocity in the direction of flow through the water pathway can be then determined based on data indicative of the detected feature(s). A mapping of the water pathway can be thereafter generated utilizing the velocity with respect to the features detected by the sensors.

Abstract: Embodiments of the invention provide a system and method for partitioning data for modelling a geological structure including laterally partitioning the data into multiple columns, where each column may be substantially laterally centered about a dual fiber and vertically partitioning each column at each of multiple intersection points of multiple surfaces and the dual fiber about which the column may be substantially laterally centered.

Abstract: A system and method for generating images of a subsurface region of interest is provided. In one embodiment, a computer-implemented method of generating images related to a subsurface region of interest includes: accessing, via a central processing unit (CPU), seismic data and an earth model related to the subsurface region of interest; forward propagating a source wavefield using the earth model at a first time interval via at least one external co-processor coupled to the CPU; transferring, at a second time interval, the forward propagated source wavefield to the CPU for compression and external storage; backward propagating the seismic data at the first time interval via the external co-processor to derive backward propagated receiver wavefield; and transferring, at the second time interval, the backward propagated receiver wavefield to the CPU.

Abstract: A step frequency inverse synthetic aperture radar (ISAR) includes a transmitter configured to transmit a transmission pulse at a transmission frequency to a near earth object (NEO), the transmission frequency having a frequency range comprising a starting frequency, an ending frequency, and a step size; a receiver configured to receive a pulse response from the NEO, the pulse response corresponding to the transmission pulse; and a computer configured to determine a 3-dimensional image of the interior of the NEO from the pulse response.

Abstract: Method, system, and computer program product for tracing seismic sections to convert to digital format are provided. The method includes: receiving a scan of a seismic section having multiple overlapping traces measured at spot points; splitting the traces into segments; associating the segments with appropriate baselines; estimating the traces between segments using logical rules; and deriving amplitudes of the traces.

Abstract: Method for obtaining rock parameters such as porosity and vshale directly from inversion of seismic data corresponding to a single trace location. This method is distinguished from existing methods that obtain elastic properties from inversion of seismic data, then relate the elastic parameters to rock lithology parameters such as porosity or vshale because it is accomplished in one step, can incorporate anisotropy and does not require multiple trace locations for stability. The data are separated into partial stacks, and a wavelet is specified for each stack. A set of linearized equations are constructed relating seismic reflectivity to changes in elastic parameters, and another set of linearized equations is constructed relating the changes in elastic parameters to the lithologic parameters. The linearized reflectivity equations are combined with the linearized rock physics equations, convolved with the specified wavelets, and equated to the seismic data.

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: Embodiments of the invention provide a system and method for converting coordinate systems for representing image data such as for example seismic data, including accepting a first set of seismic data, mapping the first set of seismic data to a second set of seismic data, where the dimensionality of the second set of seismic data is less than the dimensionality of the first set of seismic data, and generating image data by processing the second set of seismic data.

Abstract: The invention relates to acquiring seismic data in either land or marine environments, but typically marine environments where a pulse-type source is fired in a distinctive composite pulse like a distinctive rumble. In a preferred embodiment, a number of pulse-type seismic sources, sometimes called an array, are fired in a distinctive composite pulse to be able to identify within the returning wavefield the energy resulting from the composite pulse. Firing the pulse-type sources creates an identifiable signature so that two or more marine seismic acquisition systems with source arrays can be acquiring seismic data concurrently and the peak energy delivered into the water will be less, which will reduce the irritation of seismic data acquisition to marine life.

Abstract: A method for performing chrono-stratigraphic interpretation of a subterranean formation. The method includes obtaining a seismic volume containing stratigraphic features of the subterranean formation deformed by structural events, performing structural restoration of the seismic volume to generate a restored seismic volume by removing deformation due to the structural events, performing a chrono-stratigraphic interpretation based on the restored seismic volume to generate chrono-stratigraphic objects each associated with a respective relative geologic age, and displaying the chrono-stratigraphic objects in a chrono-stratigraphic space according to the respective relative geologic age of each of the stratigraphic objects.

Abstract: A method for performing seismic interpretation of a subterranean formation by enabling dual-domain interpretation of seismic features in the present day depth domain and simultaneously in a structurally restored “mapped” seismic domain. Specifically, seismic interpretation is performed on structurally restored seismic volumes while concurrently viewing the interpretation results in the structural domain. This increases the interpretation confidence by improved correlation of structural deformed events to their pre-structurally deformed geometry. The method includes the ability to progressively remove structural deformation from the seismic volume, corresponding to moving back in geologic time. The interpretation can be performed in either domain (present-day structure or structurally restored to a defined geologic event), and the interpretation will be mapped to any other computed geologic age or present structural age.

Abstract: A method and system of detecting and mapping a subsurface hydrocarbon reservoir includes determining ratio data for a plurality of orthogonal spectral components of naturally occurring low frequency background seismic data. The ratio data may be compared, plotted, contoured and displayed as a subsurface hydrocarbon reservoir map or a hydrocarbon potential map. The ratio data may represent a vertical spectral component of the seismic data over a horizontal spectral component of the seismic data. The subsurface hydrocarbon reservoir map may include contouring the ratio data over a geographical area associated with the seismic data.

Abstract: According to a preferred aspect of the instant invention, there is provided herein a system and method for Interpolation of seismic data with a POCS (projection onto convex sets) algorithm that can produce high quality interpolation results, at a reduced computational cost. In particular, optimizing the threshold calculations, allowing for aliased data, allowing the algorithm to skip unneeded Fourier transforms, and parallelization makes the method more practical and robust. Robust interpolation of aliased data may be done by limiting the calculations of the beginning iterations to the lowest temporal and spatial frequencies included in the data, then gradually allowing higher temporal and spatial frequencies into the following iterations.

Abstract: A method and system of detecting and mapping a subsurface hydrocarbon reservoir includes determining ratio data for a plurality of orthogonal spectral components of naturally occurring low frequency background seismic data. The ratio data may be compared, plotted, contoured and displayed as a subsurface hydrocarbon reservoir map or a hydrocarbon potential map. The ratio data may represent a vertical spectral component of the seismic data over a horizontal spectral component of the seismic data. The subsurface hydrocarbon reservoir map may include contouring the ratio data over a geographical area associated with the seismic data.

Abstract: A method and system for processing synchronous array seismic data includes acquiring synchronous passive seismic data from a plurality of sensors to obtain synchronized array measurements. A reverse-time data propagation process is applied to the synchronized array measurements to obtain a plurality of dynamic particle parameters associated with subsurface locations. Imaging conditions are applied to obtain imaging values that may be summed or stacked to obtain a time reverse image attribute. A volume of imaging values may be scaled by a non-signal noise function to obtain a modified image that is compensated for noise effects.

Abstract: A method and system for processing synchronous array seismic data includes acquiring synchronous passive seismic data from a plurality of sensors to obtain synchronized array measurements. A reverse-time data propagation process is applied to the synchronized array measurements to obtain a plurality of dynamic particle parameters associated with subsurface locations. Imaging conditions are applied to obtain imaging values that may be summed or stacked to obtain a time reverse image attribute. A volume of imaging values may be scaled by a non-signal noise function to obtain a modified image that is compensated for noise effects.

Abstract: Systems, program product, and methods of suppressing residual water bottom energy in seismic data, are provided. An example of a system, program product, and method can be applied to post-stack datasets and can combine multi-channel deconvolution with novel sorting keys to efficiently identify and suppress residual water bottom energy in common depth point (CDP) stacked seismic data, thereby increasing the resolving power of seismic data leading to an improved interpretation of seismic signals reflected from oil reservoirs.

Abstract: A system, which may be and/or may include a computer system, and method for identifying a fault patch position from a seismic data volume. The system and method may include defining an initial active surface in the seismic data volume; identifying one or more active surfaces by minimizing the value of a function for energy; tracking back the active surfaces to find a set of control points; and producing the fault patch position formed by the set of control points.

Abstract: Seismic data are assembled and stored for a set of cross-correlation lag times to form an array of common image gathers over depth levels of interest. The two dimensional gathers assembled over different lag times form a three-dimensional cube of common image data. The data are analyzed to determine the travel time shift required to equalize upgoing and downgoing wavefields. Events in the common image gathers are then modeled using Green's functions to generate a data set representing the data resulting from processing had a precise velocity model been obtainable from the seismic data. The generated data are then processed with inversion techniques to form a velocity model for seismic data analysis.

Abstract: A method for processing seismic data. The method includes plotting a travel-time curve and one or more actual times that each seismic receiver receives one or more seismograms on a graph. The travel-time curve represents an approximate time that each seismic receiver in a seismic survey area receives the seismograms. The method may further include determining a residual static shift for each seismogram and using the residual static shift to estimate a time shift correction. The residual static shift represents a difference between an actual time in which a seismic receiver received a seismogram and a time on the travel-time curve that corresponds to the actual time. The time shift correction corrects a timing error in the time indicated by the travel-time curve.

Abstract: A long-range lightning detection and characterization system and method. Electromagnetic radiation produced by a lightning strike is sensed at a sensing location and a measured signal representative of the strike is produced that defines an amplitude versus time. A set of reference data containing waveforms comparable to the measured signal is established. The set of reference data defines a set of reference amplitudes versus time, representative of one or more predetermined reference lightning strikes at various predetermined distances, where for each of the distances, the reference lightning strikes have associated therewith various predetermined day/night percentages. The measured signal is compared with the reference data, the comparison including finding a correlation between the measured signal and a correlating waveform in the reference data.

Abstract: Application of nonlinear resonance interferometry is introduced as a new geophysical approach to improve predictability in characterization of subsurface microseismic event analysis and propagation of fracture. In contrast to reflection methods that remove random information noise, nonlinear resonance interferometry exploits the full microseismic acquisition spectrum. In some examples, systems and techniques implement novel computational interactions between acquired microseismic wavefield attributes and a nonlinear system in software to amplify distortions in microseismic noise and exploits injection of synthetic noise, in software format, to fracture events.

Abstract: To perform a marine survey, a streamer is deployed into a body of water, where the streamer has plural first sensors to perform a subterranean survey. Indications are received from second sensors in corresponding sections of the streamer, where the indications are regarding which sections are in the body of water. Information is updated regarding which sections of the streamer are in the body of water in response to the received indications.

Abstract: A method and system for processing synchronous array seismic data includes acquiring synchronous passive seismic data from a plurality of sensors to obtain synchronized array measurements. A reverse-time data propagation process is applied to the synchronized array measurements to obtain a plurality of dynamic particle parameters associated with subsurface locations. Imaging conditions are applied to obtain imaging values that may be summed or stacked to obtain a time reverse image attribute. A volume of imaging values may be scaled by a non-signal noise function to obtain a modified image that is compensated for noise effects.

Abstract: A method and system for processing synchronous array seismic data includes acquiring synchronous passive seismic data from a plurality of sensors to obtain synchronized array measurements. A reverse-time data propagation process is applied to the synchronized array measurements to obtain a plurality of dynamic particle parameters associated with subsurface locations. Imaging conditions are applied to obtain imaging values that may be summed or stacked to obtain a time reverse image attribute. A volume of imaging values may be scaled by a non-signal noise function to obtain a modified image that is compensated for noise effects.

Abstract: To determine fault transmissivity in a subterranean reservoir, fault structures in the subterranean reservoir are identified, and elastic impedance data is computed from well log data and seismic data. Normalization parameters that have a predetermined orientation with respect to the identified fault structures are computed, and elastic impedance flux values are normalized with respect to the normalization parameters, where the normalized elastic impedance flux values are indicative of transmissivities of fault structures in the subterranean structure.

Abstract: To perform a marine survey, a streamer is deployed into a body of water, where the streamer has plural first sensors to perform a subterranean survey. Indications are received from second sensors in corresponding sections of the streamer, where the indications are regarding which sections are in the body of water. Information is updated regarding which sections of the streamer are in the body of water in response to the received indications.

Abstract: Method, apparatus, computer program and data carrier for modeling a multidimensional, heterogeneous structure using a digital processing unit, by means of a grid built up of stacks of multidimensional cells. A cell is bounded by boundary surfaces and represents at least one property of the structure. The composition of a stack, the position of the stacks in the grid and the properties represented by cells are stored in memory means by the processing unit. For modeling a discontinuity, a respective stack is divided into two or more substacks separated by separating surfaces. The composition and the position of the substacks in a respective stack are stored in the memory means. The modeling method enables an accurate and fast computation of properties of the modeled structure or space.

Abstract: A system and method monitor a hydrocarbon reservoir for drainage in volumetric three dimensions. Monitoring between wells is imperative for optimum reservoir management and is achieved by mapping the hydrocarbon fluid pathways in a producing reservoir. Unlike conventional 4D or time-lapse reflection seismic imaging systems that use a controlled active seismic source and records reflected seismic energy at receivers, the system and method exploit the minute vibrations, or micro-earthquakes generated in the reservoir layers that are induced by fluid movement. These microseisms are detected as the fluids move in the reservoir.

Abstract: Described herein is a method for optimizing a plurality of calibration maps for an algorithm of estimation of a control quantity of an internal combustion engine, each of the maps comprising a plurality of calibration values of said control quantity estimated by said algorithm. The optimization method comprises measuring the control quantity, estimating the control quantity, and individually optimizing each calibration map based on the measured control quantity and the estimated control quantity.