One of the most important applications of gravity surveys in regional geophysical studies is determining the depth to basement. Conventional methods of solving this problem are based on the spectrum and/or Euler deconvolution analysis of the gravity field and on parameterization of the earth's subsurface into prismatic cells. We have developed a new method of solving this problem based on 3D Cauchy-type integral representation of the potential fields. Traditionally, potential fields have been calculated using volume integrals over the domains occupied by anomalous masses subdivided into prismatic cells. This discretization can be computationally expensive, especially in a 3D case. The technique of Cauchy-type integrals made it possible to represent the gravity field and its gradients as surface integrals. In this approach, only the density contrast surface between sediment and basement needed to be discretized for the calculation of gravity field. This was especially significant in the modeling and inversion of gravity data for determining the depth to the basement. Another important result was developing a novel method of inversion of gravity data to recover the depth to basement, based on the 3D Cauchy-type integral representation. Our numerical studies determined that the new method is much faster than conventional volume discretization method to compute the gravity response. Our synthetic model studies also showed that the developed inversion algorithm based on Cauchy-type integral is capable of recovering the geometry and depth of the sedimentary basin effectively with a complex density profile in the vertical direction.

This paper presents a linear edge-based finite element method for numerical modeling of 3D controlledsource electromagnetic data in an anisotropic conductive medium. We use a nonuniform rectangular mesh in order to capture the rapid change of diffusive electromagnetic field within the regions of anomalous conductivity and close to the location of the source. In order to avoid the source singularity, we solve Maxwell's equation with respect to anomalous electric field. The nonuniform rectangular mesh can be transformed to hexahedral mesh in order to simulate the bathymetry effect. The sparse system of finite element equations is solved using a quasi-minimum residual method with a Jacobian preconditioner. We have applied the developed algorithm to compute a typical MCSEM response over a 3D model of a hydrocarbon reservoir located in both isotropic and anisotropic mediums. The modeling results are in a good agreement with the solutions obtained by the integral equation method.

No abstract

In induced-polarization (IP) surveys, the raw data are usually distorted significantly by the presence of electromagnetic (EM) interferences, including cultural noise. Several methods have been proposed to improve the signal-to-noise ratio of these data. However, signal processing in an electromagnetically noisy environment is still a challenging problem. We have determined a new and simple technique based on the analysis of the correlation between the measured potential and the injected primary current signals. This processing is applied to the data acquired using a new frequency-domain IP method called the spread-spectrum induced-polarization (SSIP) approach. In this approach, we use a pseudorandom m-sequence (also called the maximum length sequence) for the injected primary current. One of the advantages of this sequence is to be essentially spectrally flat in a given frequency range. Therefore, complex resistivity can be determined simultaneously at various frequencies. A new SSIP data set is acquired in the vicinity of Baiyin mine, Gansu Province, China. The correlation between potential difference and transmitting current signals for each period can be used to assess data quality. Only when the correlation coefficient between the two signals is greater than 0.5 can the SSIP data be used for subsequent processing and tomography. We determine what threshold value should be used for the correlation coefficient to extract high-quality apparent complex resistivity data and eliminate EM-contaminated data. We then compare the pseudosections with and without using the correlation analysis. When the correlation analysis is used, the noisy data are filtered out, and the target anomaly obtained through tomography is clearly enhanced. The inversion results of the apparent complex resistivity (amplitude and phase) for the survey area are consistent with some independent geologic and drilling information regarding the position of the ore body demonstrating the effectiveness of the approach.