Abstract: A method for identifying and suppressing water column reverberations (“multiple reflections”) in two-component ocean bottom seismic data is disclosed. The method involves processing the hydrophone (P) data and the geophone (Z) data separately to produce two stacked images of the subsurface (21). Analyzing the stacked P-image and the stacked Z-image together can be used to identify multiple reflections (22). Analyzing the stacked /?-image and the stacked Z-image together with an image of the subsurface created from hydrophone and geophone data combined in the usual way (PZ-image) (20) can be used to identify residual multiples in the PZ image (23). The stacked P and Z images can be combined using an existing PZ combination technique to suppress multiples (24). The efficiency of the PZ combination technique at suppressing multiples is increased because of the higher signal-to-noise ratio of stacked data.

Abstract: Method and apparatus for producing a bubble curtain with a diversity of bubble diameters for purposes such as suppressing surface-related multiple reflections in marine seismic surveys. Bubble generating elements are used that combine porous wall material with discrete holes.

Abstract: Method for reducing air wave and/or magnetotelluric noise in controlled source electromagnetic surveying by either shielding the source (61) from the air interface, shielding the receivers from downward traveling electromagnetic energy, or by employing a second source (62) to preferentially cancel the air wave (and MT) part of the signal, or a combination of the preceding.

Abstract: Method for determining best and worst cases for values of model parameters such as porosity and shale volume fraction generated by non-unique matrix inversion of physical data such as seismic reflection amplitudes. The matrix is diagonalized, and then orthonormal basis vectors associated with insignificant diagonal elements are used to generate upper and lower bounds on the solution. Best and worst case solutions are determined as linear combinations of the null basis vectors, where the expansion coefficients are determined by making a best fit to the upper and lower bounds.

Abstract: Method for efficient inversion of controlled-source electromagnetic survey data to obtain a resistivity model of the subsurface of the survey area. The method extracts the dimensions and location of sub-surface structures as they may be revealed by existing seismic or other available high resolution survey data from the subsurface area (33). This structure geometry information is used to construct a discretization (grid, or mesh) for the inversion computation (34) that is different from the mesh used for the forward modeling calculations (32) in that (a) it has fewer and hence larger cells; and (b) the cells honor the assumed structural information. The inversion need only extract resistivity information (35), the geometry of the resistive structures being specified by the inversion mesh.

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: One or more techniques for determining time-varying stress and strain fields within a subsurface region include integrating a seismic model (110) of a reservoir within the subsurface region with a geomechanical model (140) of the subsurface region. An estimate of the time-varying stress and strain fields within the subsurface region during production of the reservoir are determined, wherein the estimate is based on the integration of the seismic model with the geomechanical model. The integration of the seismic model with the geomechanical model can be used to predict the feasibility of passive seismic monitoring for a reservoir within a subsurface region (170).

Abstract: Method for identifying, determining and correcting source-related phase errors in data from a controlled source electromagnetic survey by using data from ordinary survey receivers, i.e. without benefit of source monitoring data. Abrupt anomalies indicating source malfunctions are identified (71) in the time domain by plotting time intervals between neighboring zero crossings or by zero-lag cross correlation between consecutive bins of receiver data, and the amount of the time error (73) can be determined by performing cross correlation between two bins on either side of an anomaly. In the frequency domain, transmitter anomalies can be identified by looking for discontinuities in plots of phase vs. offset, and the corrective phase shift can be determined by matching the phase on one side of the anomaly to that on the other side. A global time/phase shift (76) can be determined by using phase frequency-scaling behavior at near offsets.

Abstract: The method for suppressing noise in the controlled seismic source electromagnetic survey data based on the frequency content of the noise (51) contained in the data. The invention recognizes that some data variations across bins cannot be attributed to resistivity variations within the earth. This variation across the bins constitutes a model of noise in such surveys, and the invention mitigates noise that obeys this model.

Abstract: This invention relates generally to a method of simulating the signal of an electromagnetic source using one or more dipole sources. In the method a dipole source is located at an excitation location corresponding to a segment of the electromagnetic source to be simulated. The dipole source is activated, and an electromagnetic signal recorded at one or more receiver locations. This process is repeated for additional excitation locations corresponding to additional segments of the electromagnetic source. The data from the sequence of dipole source excitation locations is processed to determine the simulated signal of the electromagnetic source.

Abstract: A method is disclosed for combining seismic data sets. This method has application in merging data sets of different vintages, merging data sets collected using different acquisition technologies, and merging data sets acquired using different types of sensors, for example merging hydrophone and geophone measurements in ocean bottom seismic data.

Abstract: Method for conducting a seismic survey (23) wherein the acquisition parameters are modified as the survey progresses using information (24) from sensors located to monitor structures and fixtures located in the survey area. The sensor output is compared to acceptable levels, and acquisition parameters such as the seismic source strength are adjusted (26) to prevent exceeding the acceptable levels. An automated arbitration scheme (22) is used to resolve conflicting system goals within survey resource constraints.

Abstract: A structured computer-implemented method (1000) based on graphical user interfaces for interpretation and mapping of data from a controlled-source electromagnetic survey, featuring capability to store electromagnetic data in layers (1004), each layer having the same internal structure for ease of comparison, editing, display and other manipulation by an assortment of software tools (1005) useful to an interpreter. Thus, different layers might contain actual data (1001) from different surveys and simulated results (1002) based on different resistivity models or inversion results (1003), all pertaining to the same subterranean survey region.

Abstract: Method for generating a three-dimensional resistivity data volume for a subsurface region from an initial resistivity model and measured electromagnetic field data from an electromagnetic survey of the region, where the initial resistivity model is preferably obtained by performing multiple ID inversions of the measured data [100]. The resulting resistivity depth profiles are then registered at proper 3D positions [102]. The 3D electromagnetic response is simulated [106] assuming the resistivity structure is given by the initial resistivity model. The measured electromagnetic field data volume is scaled by the simulated results [108] and the ratios are registered at proper 3D positions [110] producing a ratio data volume [112]. A 3D resistivity volume is then generated by multiplying the initial resistivity volume by the ratio data volume (or some function of it), location-by location [114]. A related method emphasizes deeper resistive anomalies over masking effects of shallow anomalies.

Abstract: Method and computer program for accepting controlled-source electromagnetic (“CSEM”) source and receiver data (40) as time series, transforming these data into the time-frequency domain, and reducing these data and survey metadata to a form suitable for interpretation or inversion. The invention includes: a number of processing tools or programs (30), each designed to take a specific action on CSEM data or metadata, combine data types in some way, and/or provide a visual representation of data; a Graphical User Interface (32) to specify the action of specific tools (34) on specific data, supply parameters to tools, and monitor progress of the processing project; and a specified common internal data format, so that processing tools may be applied in various orders (36) during different processing flows and processed CSEM data can be passed on to interpretation or inversion systems.

Abstract: Embodiments of the invention relate to a process for providing a pumpable hydrate slurry in a hydrocarbon pipeline fluid mixture having a water-cut greater than about 50 volume percent. In one or more embodiments the process comprises treating the fluid mixture with an anti-agglomerant and adding water to the fluid mixture in an amount sufficient to lower the gas-water ratio sufficiently to achieve a pumpable hydrate slurry. Also disclosed are methods for producing hydrocarbons utilizing a process for providing a pumpable hydrate slurry in a hydrocarbon pipeline fluid mixture having a water -cut greater than about 50 volume percent.

Abstract: A method, apparatus and computer program for improving the signal-to-noise ratio of a signal S(t), S(t) containing Signal and noise, are disclosed. A measurement of S(t) at a frequency-of-interest is obtained. Noise measurements of S(t) at one or more noise frequencies where the Signal portion of S(t) is expected to be small are obtained. The noise at the frequency-of-interest is estimated using the noise measurements at the one or more noise frequencies. The estimated noise is subtracted from the measurement of S(t) at the frequency-of-interest.

Abstract: The invention is a method for simulating sandstone deposition. The sandstone is simulated by estimating the grain size distribution and mineral composition of grains in the sandstone, simulating sedimentation of grains from the grain size distribution and mineral composition of the grains, simulating compaction of the grains, and simulating cementation of the grains. Properties of the sandstone such as porosity and permeability may be estimated from the simulated sandstone. The method permits multiple mineralogies to be simulated during the burial history of sedimentation, compaction and cementation.

Abstract: Method for conducting an efficient and interpretable controlled-source electromagnetic reconnaissance survey for buried hydrocarbons. While a part of the survey area is being set up for measurement and data are being acquired, data from a nearby part of the survey area, surveyed just previously, are being rapidly processed and analyzed. If the analysis shows resistive anomalies of interest in a portion of a survey area, a fine-grid survey is quickly designed for that portion, and that survey is conducted next before moving source and receivers to a more distant part of the survey area.