Abstract: Diagonal gather trace interpolation systems and methods are disclosed. In some embodiments, the method includes obtaining seismic traces from a marine seismic survey performed using typically two seismic sources fired in a flip-flop pattern. The seismic traces are processed in common field file trace groups from each cable by performing trace interpolation in opposing diagonal directions. Among other things, diagonal trace interpolation may improve azimuthal regularization. In any event, the disclosed interpolation methods maintain spatial bandwidth increased spatial resolution with increased inline and crossline sampling components. Diagonal gather trace interpolation exploits reduced input trace separation to provide improved stability and detection of a greater range of formation dip angles. After interpolation and regularization, the seismic traces may be imaged and interpreted for improved seismic exploration and monitoring.

Abstract: A method and system of processing seismic data includes acquiring three-component seismic data and combining horizontal components of the three-component seismic data to obtain a merged horizontal component. Frequency spectra are determined for the acquired three-component seismic data. A ratio of a vertical component of the seismic data to the merged horizontal component is determined. A V/H integration-measure is obtained from the integration of the area bounded by a selected minimum-amplitude value and the amplitude values greater than the selected minimum-amplitude value and the V/H integration-measure is stored in a form for display.

Abstract: A method of processing seismic data includes storing in a single file structure a first type of data representing an inline location, a second type of data representing a crossline location, and a third type of data representing at least one of a depth and time, said first, second and third type of data each being associated with a horizon of a poststack seismic trace at a location defined by said inline location and said crossline location; and storing in said single file structure at least a fourth type of data representing an attribute associated with said horizon of said poststack seismic trace at said location.

Abstract: Method for reducing instability and increasing computational efficiency in tomographic inversion for velocity model building. A system of tomographic equations is developed for a uniform grid. A non-uniform parameterization is found for which a linear mapping exists between the space of the uniform grid and the space of the non-uniform grid. The matrix that relates velocity to the tomographic data in the non-uniform representation is then given by the matrix product of the corresponding matrix in the uniform grid representation and the mapping matrix. Inversion can then be performed for the non-uniform parameterization on a smaller, more stable matrix.

Abstract: In a first aspect, a method for use in a time lapse, marine seismic survey includes accessing a set of baseline seismic data; accessing a set of acquired, time lapse seismic data; and interpolating a set of time lapse seismic data from the baseline seismic data and the acquired time lapse seismic data, at least one of the baseline seismic data and the acquired time lapse seismic data being multicomponent data. In other aspects, a program storage medium is encoded with instructions that, when executed by a computing device, perform the above method and a computing apparatus programmed to perform one or more of such methods.

Abstract: Techniques and a system for performing geological interpretation operations in support of energy resources exploration and production perform well log correlation operations for generating a set of graphical data describing the predetermined geological region. The process and system interpret the geological environment of the predetermined geological region from measured surface and fault data associated with the predetermined geological region. Allowing the user to query and filter graphical data representing the predetermined geological region, the method and system present manipulable three-dimensional geological interpretations of two-dimensional geological data relating to the predetermined geological region and provide displays of base map features associated with the predetermined geological region.

Abstract: Solutions to the problem of reversing seismic fault movements are formulated using a model based on elasticity theory, and using finite element and boundary element methods for generating a solution. The solution involves defining slip vectors from known formations in the fault and applying a space constraint restriction to traction values on the fault surface. The method may be applied in either 2D or 3D applications. The method is computationally fast enough to allow interactive fault reversal, and permits experimentation with various unfaulting scenarios so that a geologically acceptable solution is provided.

Abstract: A method for hydraulic fracture imaging consistently integrates ground deformation data and microseismic data. A time sequence of hydraulic fracture images is computed from the combined data by first computing an orientation and location of the fracture. A size of the fracture (e.g., length, width, and opening) and detailed shape is then computed consistent with the computed orientation using a stochastic optimization procedure. This procedure may be performed by selecting a candidate model of the fracture consistent with the computed orientation, computing a predicted surface deformation from the candidate model (e.g., using elastic crack theory), computing a predicted seismic activity from the candidate model (e.g., using Dieterich seismicity rate theory), and comparing the predicted surface deformation and predicted seismic activity with the collected ground deformation data and earthquake locations derived from the collected seismic activity data.

Abstract: A method and system for processing three component seismic data includes determining a covariance data matrix from three component seismic data for each of a plurality of time periods to obtain eigenvectors and eigenvalues. One or more seismic attributes are calculated from the eigenvectors and eigenvalues for each of a plurality of time periods. A stability measure is determined from the calculated seismic attribute for each of the plurality of time periods and stored for display.

Abstract: The analysis of scattering diagrams of the correlation-type representation theorem in inhomogeneous media is improved by the use of virtual events. Virtual events here are events which are not directly recorded in standard seismic data acquisition, but the assumption of their existence permits the construction, of internal multiples with scattering points at the sea surface.

Abstract: A computer implemented method of producing macro-view conclusion data related to a seismic event. The method includes receiving at a processing center a set of data related to a parameter of a seismic event. The parameter is measured using a vibration sensor of a hard disk drive of a computer at a known, fixed location, wherein the vibration sensor produces the set of data. The method further includes processing the set of data at the processing center to produce the macro-view conclusion data. The macro-view conclusion data is then stored in a storage device.

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: The invention relates to a method of performing oilfield operations for an oilfield, the oilfield having a subterranean formation. The method includes collecting a first volume data set of seismic data and a second volume data set of seismic data, co-rendering a visually-melded scene directly from the first volume data set and the second volume data set, displaying the visually-melded scene comprising a visualized geobody, where the visualized geobody corresponds to a portion of the first volume data set and the second volume data set, and selectively adjusting the oilfield operations based on the visualized geobody.

Abstract: The invention is a method for stratigraphic interpretation of a seismic image for determining a sedimentary history of an underground zone having application to petroleum exploration. The seismic reflectors of the image are digitized by means of an applied shape recognition technique. At least one attribute characteristic of the structure of the image in the vicinity of each reflector is calculated for a set of reflectors by means of an image analysis technique. This analysis is based on at least one analysis window whose shape is controlled by the associated reflector. Pertinent reflectors, in the sense of the stratigraphic interpretation, are selected from among this set of reflectors on the basis of this attribute. Finally, the sedimentary history of the underground zone is reconstructed from a stratigraphic interpretation of the reflectors thus selected.

Abstract: A method for analyzing seismic data from a formation. The method may first receive two-dimensional (2-D) seismic data of a formation comprising a number of individual seismic lines acquired over an exploration area and having a large interline spacing. The 2-D seismic data may then be preprocessed to enhance diffracted energy. For each respective hypothetical diffractor location, the method may then search for coherent diffraction arrivals on nearby 2-D seismic lines consistent with the respective hypothetical diffractor location. The method may then store information regarding identified diffractor locations. The method may then create a map based on the identified diffractor locations, wherein the map illustrates areas of high diffraction. The map may then be displayed on a display, wherein the map is useable to assess the formation.

Abstract: Diagonal gather trace interpolation systems and methods are disclosed. In some embodiments, the method includes obtaining seismic traces from a marine seismic survey performed using typically two seismic sources fired in a flip-flop pattern. The seismic traces are processed in common field file trace groups from each cable by performing trace interpolation in opposing diagonal directions. Among other things, diagonal trace interpolation may improve azimuthal regularization. In any event, the disclosed interpolation methods maintain spatial bandwidth increased spatial resolution with increased inline and crossline sampling components. Diagonal gather trace interpolation exploits reduced input trace separation to provide improved stability and detection of a greater range of formation dip angles. After interpolation and regularization, the seismic traces may be imaged and interpreted for improved seismic exploration and monitoring.

Abstract: Information from a network of seismic device is received, and a geographic representation of the network of seismic devices is generated. The geographic representation contains information regarding the seismic devices.

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: Method for identifying one or more fracture clusters in a formation surrounding a reservoir. In one implementation, the method may include generating a P to S image, comparing the P to S image to one or more images from a borehole, and identifying one or more fracture clusters using the P to S image and the borehole images.

Abstract: Systems and methods are disclosed to construct subsurface images from diffuse seismic energy. Various disclosed system embodiments include multiple seismic sensors that each convert received seismic energy into one or more seismic signals. One or more processor combine the seismic signals to determine a subsurface map. As part of determining the map, the processor(s) systematically focus the array of seismic sensors on each bin in the subsurface volume of interest. In this manner each bin becomes a focal point of the array. For each bin, the processor(s) analyze the seismic wave travel time to each seismic sensor and apply a corresponding time shift to align the seismic signals with a uniform travel time. The time-shifted seismic signals are then combined to determine an intensity value for seismic energy radiating from the focal point. A subsurface map can then be derived from the intensity value as a function of position.

Abstract: A method, article of manufacture, and data set for simulating seismic prestack depth migrated images on the basis of a model of a selected GF-node, without the use of real or synthetic recorded data, provides a very efficient and flexible way to calculate a simulated depth migrated image as a function of parameters such as survey, overburden model, pulse, elastic wavefield, and local reflectivity structure. Important information needed for the method is the scattering wavenumber, calculated for example, by ray methods and other equivalent methods. Complex model geometry can be done in 2-D and 3-D.

Abstract: A method and system for characterizing a seismic reflection point. The system receives vibratory seismic data corresponding to a plurality of events, calculates a curvature and a reflection point corresponding to each of the plurality of events, and characterizes reflection points according to their corresponding curvatures.

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: 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 Radon transformation. A corrective filter is then applied to enhance the primary reflection signal content of the data and to eliminate unwanted noise events. The corrective filter has a pass region with a lower pass limit and a higher pass limit. The higher pass limit is set within 15% above the slowness of the primary reflection signals and, preferably, it is more closely set to the slowness of the primary reflection signals. The lower pass limit is also preferably set within 15% below the slowness of the primary reflection signals.

Abstract: Systems and methods for imaging waveform volumes. An image of the waveform volume may be drawn on a display device as a two-dimensional image or a three-dimensional image of a sampling probe and redrawn in real-time at interactive rates using a graphics accelerator or a graphics card. The image of the waveform volume may also include seismic-data traces that are color-filled according to texture coordinates for pixels on the display device that intersect the waveform volume.

Abstract: Method for removing effects of shallow resistivity structures in electromagnetic survey data to produce a low frequency resistivity anomaly map, or alternatively imaging resistivity structures at their correct depth levels. The method involves solving Maxwell's electromagnetic field equations by either forward modeling or inversion, and requires at least two survey data sets, one taken at the source frequency selected to penetrate to a target depth, the other a higher frequency able to penetrate only shallow depths.

Abstract: A method to extract fault surfaces from seismic data (1,3) is disclosed. The method comprises: (a) generating at least two fault sticks (5) from the same fault from at least two slices of the seismic data wherein each slice comprises at least one fault stick from the same fault, (b) constructing an initial three-dimensional fault surface (7) containing the fault sticks, and (c) reconstructing the initial fault surface (9) using a deformable surface model to fit discontinuity or coherency information in the seismic data in an iterative process. Techniques are disclosed for constructing the initial fault surface from interpreter-provided fault nodes, and for performing the deformable surface iteration by defining an energy function for the fault surface and then minimizing the surface energy.

Abstract: A method is disclosed for reconstructing complex wave attributes from limited view measurements of a scattering object. The method involves the analytic continuation of the Fourier transform of the object function into the area in which there is an absence of K-space coverage by requiring objects to be an even function. (It is assumed that physical objects are even functions, and it is this assumption that allows analytic continuation.) When the object function is not centered at the origin, the measurements are shifted to the origin prior to determining the analytic continuation and returned to their original location following analytic continuation.

Abstract: The invention concerns a method for extracting a geological horizon and related properties of a seismic representation, comprising a step (100) which consists in digital modeling with continuous local seismic traces, calculating the optimal offset and defining a conditional neighbourhood of a reference central continuous local seismic trace; a step (101) which consists in defining a two-dimensional matrix whereof the line and column indices correspond to the coordinates of the geophones; a third step (102) which consists in selecting a seed point; a fourth step (103) which consists in determining the point vertically closest to the seed point and a fifth step (104) which consists in assigning to the point P(p,q,t) the value P(p,q,t+hij,pq,k), where hij,pq,k is optimal offset of the neighbouring point P(i,j,k), so as to estimate the related properties of the conditional neighbourhood thereby filling the two-dimensional extraction matrix of step (101).

Abstract: A method and apparatus for locating hydrocarbon deposits includes temperature 16 sensor for measuring earth surface temperature at each of a plurality of locations, and a solar absorption sensor 14 for measuring the absorption coefficient of the earth's surface at each of those plurality of locations. A sky temperature sensor 18 may also be provided. A computer 30 records measurements as a function of location, and outputs a plot of corrected surface temperatures as a function of the geographic location of the measurements. By comparing corrected surface temperature measurements at the plurality of locations, a hydrocarbon deposit may be located and its boundaries determined.

Abstract: A method of locating subsurface hydrocarbon reservoirs includes acquiring seismic data, dividing the seismic data into time windows, applying a data transform to the seismic data having a plurality of components to obtain a smoothed frequency spectrum, selecting one or more local maxima and minima, determining a frequency product as a function of the maxima and minima in a form for display. Lower and higher ratio thresholds may be selected to further discriminate frequency product selection for display. A display or map of frequency products over a survey area may indicate the presence of subsurface fluid reservoirs.

Abstract: The present invention relates to a method of locating oil and gas exploration prospects. A digital model of each of the many geologic strata that compose a basin is created for all or a portion of a basin by digitizing well log data 100, collecting all digitized well log data in a data base, marking common geological time markers throughout the entire basin and visually displaying the basin in a number of useful views 380, 580 to permit the isolation or identification of oil and gas producing geological formations. This allows for the generation of oil and gas drilling prospects 600 through the extrapolation of productive depositional features into un-drilled portions of the basin. The present invention expedites these processes and allows the available well log data to be viewed in ways not previously possible thereby greatly enabling oil and gas exploration and development in a basin.

Abstract: A method for displaying a plurality of measurement data along a borehole trajectory includes displaying a borehole model representing the borehole trajectory; and displaying the plurality of measurement data as at least one layer overlaying the borehole model, wherein the plurality of measurement data are displayed at measurement depths corresponding to measurement depths of the borehole model.

Abstract: A method for locally determining a profile of geological horizons includes a step (100) which consists in digital modeling with continuous local seismic traces, calculating an optimal offset and defining a conditional neighborhood of a reference central continuous local seismic trace; a step (101) which consists in defining residuals relative to the reference central continuous local seismic trace; a third step (102) which consists in minimizing the set of residuals on the conditional neighborhood and a fourth step (103) which consists in determining parametric coefficients corresponding to the minimization of the set of residuals on the conditional neighborhood carried out at step (102).

Abstract: An system and method are disclosed for efficiently approximating analytical circuit device models. A preferred embodiment includes a method for obtaining smooth and accurate approximations of analytical device models, comprising the steps of identifying a first set of measurement units; locating two or more sets of units that neighbor one or more of said measurement units; for each set of the two or more sets of neighbor units, obtaining the union of one or more of said sets of neighbor units and the first set of measurement units; calculating the smoothness of the analytical device model within one or more of said unions; and selecting at least one of said unions within which the analytical device model is the smoothest as the new set of measurement units.

Abstract: A monitoring system for use with a horizontal directional drilling machine and adapted to determine the roll and/or pitch orientation of the downhole tool assembly by detecting signals from a sensor assembly supported by the downhole tool assembly. In a preferred embodiment, the monitoring system uses an antenna assembly to detect a signal transmitted from the sensor assembly. The antenna assembly transmits the detected signal to a processor which determines a calibration factor indicative of the actual orientation of the sensor relative to the known orientation of the downhole tool assembly. The processor then uses the calibration factor and the signal emitted from the sensor, indicative of the actual orientation of the sensor, to determine the actual orientation of the downhole tool assembly.

Abstract: A method for automated extraction of surface primitives from seismic data is presented. A preferred embodiment of the method includes defining, typically with sub-sample precision, positions of seismic horizons through an extrema representation of a 3D seismic input volume; deriving coefficients that represent the shape of the seismic waveform in the vicinity of the extrema positions; sorting the extrema positions into groups that have similar waveform shapes by applying classification techniques with the coefficients as input attributes using unsupervised or supervised classification based on an underlying statistical class model; and extracting surface primitives as surface segments that are both spatially continuous along the extrema of the seismic volume and continuous in class index in the classification volume.

Abstract: Interpretation of a shot gather and a stack of seismic data. In one implementation, the shot gather is sorted into a first component shot gather, a second component shot gather and a third component shot gather. The first component shot gather is represented as a first subimage, the second component shot gather is represented as a second subimage, and the third component is represented as a third subimage. The first subimage is then merged with the second subimage and the third subimage to create a false color image. The shot gather of seismic traces may then be interpreted based on the false color image. The first component may by an x component, the second component may be a y component and the third component may be a z component. The first subimage may be a red or cyan subimage, the second subimage may be a green or yellow subimage and the third subimage may be a blue or magenta subimage.

Abstract: A system for providing irrigation control is provided. The system includes a number of non-local data sources for providing data, a processor and an irrigation system. The processor is configured to receive data from one or more of the non-local data sources and calculate an evapotranspiration (ET) value for an irrigation area that is non-local with respect to the non-local data sources. The processor is further configured to create or alter an irrigation program based on the ET value. The irrigation system is configured to receive the irrigation program from the processor and provide appropriate irrigation control for the irrigation area using the irrigation program.

Abstract: A method for controlled source electromagnetic Earth surveying includes deploying a plurality of electromagnetic sensors in a selected pattern at the top of an area of the Earth's subsurface to be surveyed. At least one of a transient electric field and a transient magnetic field is applied to the Earth in the vicinity of the sensors at a plurality of different positions. At least one of electric field amplitude and magnetic field amplitude at each of the sensors is recorded each time the transient electric field and/or magnetic field is applied. Each recording is adjusted for acquisition geometry. An image is generated corresponding to at least one sensor position using at least two stacked, adjusted recordings.

Abstract: A method designed to combine both pseudo compressional interval velocity and pseudo shear interval velocity from seismic data for predicting fluid type and porosity in oil and gas exploration. The first step is the generation of ?S, pseudo shear reflectivity data, and P, pseudo compressional reflectivity data from the seismic data. The second step of the method is to produce scaled pseudo shear reflectivity ?Vs and scaled pseudo compressional reflectivity ?Vp from both the ?S and P data and then combine and re-project the ?Vs and the ?Vp data. The compression and shear reflectivity data are inverted, scaled to pseudo-velocity and then calibrated to well control. The resulting ?Vp and ?Vs volumes produce fluid volume and porosity volume. The fluid and porosity volumes can be stored in an electronic storage medium for analysis to distinguish between commercial and non-commercial gas accumulations prior to drilling.

Abstract: A method for identifying, imaging and monitoring dry or fluid-saturated underground reservoirs using seismic waves reflected from target porous or fractured layers is set forth. Seismic imaging the porous or fractured layer occurs by low pass filtering of the windowed reflections from the target porous or fractured layers leaving frequencies below low-most corner (or full width at half maximum) of a recorded frequency spectra. Additionally, the ratio of image amplitudes is shown to be approximately proportional to reservoir permeability, viscosity of fluid, and the fluid saturation of the porous or fractured layers.

Abstract: A method for multipath imaging of three-dimensional volumes of seismic data, given the subsurface velocity distribution. Ray tracing is used to generate the volumes of information needed for migration, but using coarse-scaled grids for computational efficiency. Interpolation methods are disclosed that address the specific shot points and receiver locations and that enable a final image on a grid with the desired resolution. Data storage and retrieval techniques are also disclosed.

Abstract: A method for creating a facies probability cube is disclosed. A S-grid containing facies is first created. The grid includes layers of cells and columns of cells. Vertical facies proportion data for the layers of cells is derived from sources such as well data, conceptual vertical geologic sections and graphs describing the proportion of facies found in each of the layers. Areal depocenter maps are created which have defined boundaries to create depocenter regions in which respective facies are likely to occur. Ideally, the boundaries for the depocenter regions for each of the facies are sequentially and independently defined. Horizontal or map facies proportion data are generated for the columns of cells preferably by filtering the depocenter regions so that the proportion of each of the facies ranges from a maximum value to a minimum value.

Abstract: A computer implemented method for interpreting faults from a fault-enhanced 3-D seismic attribute cube. The method includes the steps of extracting faults from a 3-D seismic attribute cube, and of calculating a minimum path value for each voxel of the 3-D seismic attribute cube. A fault network skeleton is extracted from the 3-D seismic attribute cube by utilizing the minimum path values which correspond to voxels within the 3-D seismic attribute cube. The individual fault networks are then labeled, and a vector description of the fault network skeleton is created. The fault network skeleton is subdivided into individual fault patches wherein the individual fault patches are the smallest, non-intersecting, non-bifurcating patches that lie on only one geologic fault. The individual fault patches are then correlated into a representation of geologic faults.

Abstract: A method of measuring a water effect on an underground non-water well includes the steps of assembling a history of water cut (WC) measurements for the well over an extended period of time, determining a water cut rate of change (R) for the well based upon the history of water cut measurements, wherein the water cut rate of change is defined as follows: ? ( WC ) ? ( t ) = R , and identifying whether an anomaly exists based on the determined R.

Abstract: The present invention relates to a method of locating oil and gas exploration prospects. A digital model of each of the many geologic strata that compose a basin is created for all or a portion of a basin by digitizing well log data 100, collecting all digitized well log data in a data base, marking common geological time markers throughout the entire basin and visually displaying the basin in a number of useful views 380, 580 to permit the isolation or identification of oil and gas producing geological formations. This allows for the generation of oil and gas drilling prospects 600 through the extrapolation of productive depositional features into un-drilled portions of the basin. The present invention expedites these processes and allows the available well log data to be viewed in ways not previously possible thereby greatly enabling oil and gas exploration and development in a basin.

Abstract: In one embodiment, a computer-implemented seismic viewing/editing collaboration method includes performing real-time collaborative cursor tracking, copaging, picking, and image manipulation in a distributed-display-processing, peer-to-peer architecture. A parameterized, minimal set of information required to update a display is transferred directly between different clients. A group state containing events generated by different clients is enforced to be synchronized on the different clients.

Abstract: The invention is a method for performing a stratigraphically-based seed detection in a 3-D seismic data volume. The method incorporates criteria that honor the layered nature of the subsurface so that the resulting seismic objects are stratigraphically reasonable. The method may be used to extract from a seismic data volume all seismic objects that satisfy the input criteria. Alternatively, the method may be used to determine the size and shape of a specific seismic object in a seismic data volume.

Abstract: A Method of measuring local similarities between seismic trace cubes (3D survey) obtained from a volume of an underground zone, corresponding to prestack data or to repeated seismic surveys (4D survey). For each point of the volume considered, the method comprises a) extracting, from each seismic trace cube, a volume neighborhood centered on a point, referred to as current point, and consisting of a set of seismic traces in limited number; b) applying an analysis technique referred to as (GPCA) allowing defining synthetic variables; and c) determining a coherence value from the synthetic extracted variables measuring the local similarity between the seismic trace cubes extracted from the volume neighborhood of the current point. The coherence value thus calculated is assigned to the current point. The coherence values of all of the current points form a coherence cube. An application is finer monitoring of the evolution with time of a reservoir under development.