Geological interpretation of resistivity models from marine controlled-source electromagnetic (CSEM) and magnetotelluric (MT) data for hydrocarbon exploration and reservoir monitoring can be problematic due to structural complexity and low resistivity contrasts in sedimentary units typically found in new frontier areas. It is desirable to reconstruct three-dimensional (3D) resistivity structures that are consistent with seismic images and geological expectations of the subsurface in order to reduce uncertainty in the evaluation of petroleum ventures. Structural similarity is achieved by promoting a cross-gradient constraint between external seismically-derived gradient fields and the inversion resistivity model. The gradient fields come from coherency weighted structure tensors computed directly from the seismic volume. Consequently, structural similarity is obtained without the requirement for any horizon interpretation or picking, thus reducing significantly the complexity and effort. We demonstrate the effectiveness of this approach using CSEM, MT, and seismic data from a structurally complex fold-thrust belt in offshore northwest Borneo.

Geologic interpretation of 3D anisotropic resistivity models from conventional marine controlled-source electromagnetic (CSEM) data inversion faces difficulties in low-resistivity contrast sediments and structurally complex environments that typify the new frontiers for hydrocarbon exploration. Currently, the typically reconstructed horizontal resistivity [Formula: see text] and vertical resistivity [Formula: see text] models often have conflicting depth structures that are difficult to explain in terms of subsurface geology, and the resulting resistivities may not be close to the true formation resistivities required for estimating reservoir parameters. We have investigated the concept that an objective geologically oriented or structurally tailored inversion can be achieved by requiring that the cross-product of the gradient of horizontal resistivity and the gradient of the vertical resistivity is equal to zero at significant geologic boundaries. We incorporate this boundary-shape criterion in our 3D inverse problem formulations, implemented within nonlinear model-space and conjugate-gradient contexts, for cases in which a priori calibration data from wells and/or seismically derived subsurface boundaries are available and for cases in which these are lacking. The resulting fit-for-purpose solutions serve to better analyze the peculiarity of a given data set. We applied these algorithms to synthetic and field CSEM data sets representing a fold-thrust environment with low-resistivity and low-contrast sediments. The resulting [Formula: see text] and [Formula: see text] models from cross-gradient joint inversion of synthetic data of appropriate frequency bandwidth without a priori information are structurally similar and consistent with the test models, whereas those from the inversions of band-limited field data are consistent with the available seismic and resistivity well-log data. This particular approach will thus be useful for lithologic correlation in frontier regions with limited a priori information using broadband CSEM data. For these band-limited field data, we found that the anisotropic bulk resistivities of the low-contrast sediments are better determined by incorporating a priori calibration data from triaxial resistivity logs and seismic horizons.

The focus of hydrocarbon exploration has now moved into frontier regions where structural complexity, heterogeneous overburden, and hydrocarbon system fundamentals are significant challenges requiring an integrated exploration approach. Three-dimensional controlled-source electromagnetic (CSEM) anisotropic resistivity imaging is emerging as a technique to combine with seismic imaging in such regions. However, the typically reconstructed horizontal resistivity [Formula: see text] and vertical resistivity [Formula: see text] models often have conflicting depth structures that are difficult to explain in terms of subsurface geology. It is highly desirable to reduce ambiguity or subjectivity in depth interpretation of [Formula: see text] and [Formula: see text] models and also achieve comparability with other coincidentally located subsurface models. We have developed a workflow for integrating information from seismic well-based inversion, interpreted seismic horizons, and resistivity well logs in a cross-gradient-guided simultaneous 3D CSEM inversion for geologically realistic [Formula: see text] and [Formula: see text] models whose parameter estimates for a selected reservoir interval can then be better optimized to aid reservoir characterization. We developed our workflow using exploration data from a complex fold-thrust belt. We found that the integrated cross-gradient approach led to [Formula: see text] and [Formula: see text] models that have a common depth structure, are consistent with seismic and resistivity logs, and are hence less ambiguous for geologic interpretation and reservoir parameter estimation.

Industry practices for near-surface analysis indicate difficulties in coping with the increased number of channels in seismic acquisition systems, and new approaches are needed to fully exploit the resolution embedded in modern seismic data sets. To achieve this goal, we have developed a novel surface-consistent refraction analysis method for low-relief geology to automatically derive near-surface corrections for seismic data processing. The method uses concepts from surface-consistent analysis applied to refracted arrivals. The key aspects of the method consist of the use of common midpoint (CMP)-offset-azimuth binning, evaluation of mean traveltime and standard deviation for each bin, rejection of anomalous first-break (FB) picks, derivation of CMP-based traveltime-offset functions, conversion to velocity-depth functions, evaluation of long-wavelength statics, and calculation of surface-consistent residual statics through waveform crosscorrelation. Residual time lags are evaluated in multiple CMP-offset-azimuth bins by crosscorrelating a pilot trace with all the other traces in the gather in which the correlation window is centered at the refracted arrival. The residuals are then used to build a system of linear equations that is simultaneously inverted for surface-consistent shot and receiver time shift corrections plus a possible subsurface residual term. All the steps are completely automated and require a fraction of the time needed for conventional near-surface analysis. The developed methodology was successfully performed on a complex 3D land data set from Central Saudi Arabia where it was benchmarked against a conventional tomographic work flow. The results indicate that the new surface-consistent refraction statics method enhances seismic imaging especially in portions of the survey dominated by noise.

No abstract