Subsalt imaging in northern Germany using multiphysics (magnetotellurics, gravity, and seismic)

Imaging subsalt is still a challenging task in oil and gas exploration. We have used magnetotellurics (MT) to improve the integration of seismic and gravity data to image the Wedehof salt dome, located in the Northern German Basin. High-density natural field source broadband MT data were acquired and enhanced the definition of the top and overhanging salt structures in addition to imaging the salt dome root. Salt boundaries show strong resistivity contrasts with the surrounding sediments and thus represent a good target for electromagnetic measurements, especially for top salt and salt flanks imaging. With integrated 3D gravity modeling focusing on the salt dome’s flanks at intermediate depths, an improved model was achieved. The new model provided sound input to a follow-up seismic depth migration that led to an improved imaging of the subsalt target proven by subsequent exploration drilling. The integrated interpretation of MT, gravity, and seismic combines the strengths of the different physics, thus increasing imaging reliability and reducing exploration drilling risks. Using a conservative workflow that included a feasibility study with field noise evaluation and careful acquisition parameter testing prior to survey start, a broadband array data acquisition, and advanced processing, the survey area's severe cultural noise issues could be overcome.

Near-Surface Multi-Point Vibration Location Method Based on Seismic Depth Migration

This paper presents a method for multi-point vibration location near-surface. Unlike traditional source location technologies, it is not using travel times of seismic waves for positioning calculation, but according to the entire seismic data record, using the wave equation migration method to calculate the source location. Similar to exploration seismic, the records from a survey line within a certain period of time are data volumes with dimensions of time and ground coordinates. Based on the data, combined with surface seismic wave propagation characteristics, by using the improved seismic depth migration algorithm, the vibration energy will return to the starting position where exactly the source location is. The method can solve the problem of location calculation error by using traditional method when several vibration at the same time or continuous vibration occurs at some point.

Reducing two-way splitting error of FFD method in dual domains

The Fourier finite-difference (FFD) method is very popular in seismic depth migration. But its straightforward 3D extension creates two-way splitting error due to ignoring the cross terms of spatial partial derivatives. Traditional correction schemes, either in the spatial domain by the implicit finite-difference method or in the wavenumber domain by phase compensation, lead to substantially increased computational costs or numerical difficulties for strong velocity contrasts. We propose compensating the two-way splitting error in dual domains, alternately in the spatial and wavenumber domains via Fourier transform. First, we organize the expanded square-root operator in terms of two-way splitting FFD plus the usually ignored cross terms. Second, we select a group of optimized coefficients to maximize the accuracy of propagation in both inline and crossline directions without yet considering the diagonal directions. Finally, we further optimize the constant coefficient of the compensation part to further improve the overall accuracy of the operator. In implementation, the compensation terms are similar to the high-order corrections of the generalized-screen method, but their functions are to compensate the two-way splitting error rather than the expansion error. Numerical experiments show that optimized one-term compensation can achieve nearly perfect circular impulse responses and the propagation angle with less than 1% error for all azimuths is improved up to 60° from 35°. Compared with traditional single-domain methods, our scheme can handle lateral velocity variations (even for strong velocity contrasts) much more easily with only one additional Fourier transform based on the two-way splitting FFD method, which helps retain the computational efficiency.