Marine vertical cable seismic (VCS) is a promising survey technique for submarine complex structure imaging and reservoir monitoring, which uses vertical arrays of hydrophones deployed near the seafloor to record seismic wavefields in a quiet environment. Recently, we developed a new type of distributed VCS system for exploration and development of natural gas hydrates preserved in shallow sediments under the seafloor. Using this system and air-gun sources, we accomplished a 3D VCS yield data acquisition for gas hydrates exploration in the Shenhu area, South China Sea. In view of the characteristics of VCS geometry, we implement reverse time migration (RTM) on a common receiver gather to obtain high-resolution images of marine sediments. Due to the unique acquisition method, it is asymmetrical for the reflection path between the sources and the receivers in the VCS survey. Therefore, we apply accurate velocity analysis to common scatter point (CSP) gathers generated from common receiver gathers instead of the conventional velocity analysis based on common depth point gathers. RTM with this reliable velocity model results in high-resolution images of submarine hydrate-bearing sediments in deep water conditions. The RTM imaging section clearly shows the bottom simulating reflector (BSR) and also the reflection characteristics of the hydrate-bearing sediments filled with consolidated hydrates. Moreover, its resolution is relative to that of acoustic logging curves from the nearby borehole, and this imaging section is well consistent with the synthetic seismogram trace generated by the logging data. All these results reveal that VCS is a great potential technology for exploration and production of marine natural gas hydrates.
Reverse time migration (RTM) is based on the two-way wave equation, so its imaging results obtained by conventional zero-lag cross-correlation imaging conditions contain a lot of low-wavenumber noises. So far, the wavefield decomposition method based on the Poynting vector has been developed to suppress these noises; however, this method also has some problems, such as unstable calculation of the Poynting vector, low accuracy of wavefield decomposition, and poor effect of large-angle migration artifacts suppression. This article introduces the optical flow vector method to RTM to realize high-precision wavefield decomposition for both the source and receiver wavefields and obtains four directions of wavefields: up-, down-, left-, and right-going. Then, the cross-correlation imaging sections of one-way propagation components of forward- and back-propagated wavefields are optimized and stacked. On this basis, the reflection angle of each imaging point is calculated based on the optical flow vector, and an attenuation factor related to the reflection angle is introduced as the weight to generate the optimal stack images. The tests of theoretical model and field marine seismic data illustrate that compared with the conventional RTM with wavefield decomposition based on the Poynting vector, the angle-weighted RTM with wavefield decomposition based on the optical flow vector proposed in this article can achieve wavefield decomposition for both the source and receiver wavefields and calculate the reflection angle of each imaging point more accurately and stably. Moreover, the proposed method adopts angle weighting processing, which can further eliminate large-angle migration artifacts and effectively improve the imaging accuracy of RTM.
Vertical seismic profiling (VSP) is an effective technique to provide high-resolution seismic images of the reservoir area. However, the quality of the images is limited by the poor illumination of primary reflection wave. In conventional VSP imaging, only the upgoing primaries are used. Adding free-surface–related multiples into the imaging process can significantly improve the coverage of the illuminated area. Conventional migration methods using multiples need the complex process of multiple prediction. Data-to-data migration (DDM) is an effective imaging technique for multiples in which the recorded data is migrated directly. To improve the imaging quality of DDM in VSP imaging, we propose separating the wavefield into downgoing and upgoing components using Hilbert transform when reverse-time migration (RTM) is implemented in DDM, and the inverse-scattering imaging condition is further applied to the decomposed wavefields. The proposed method eliminates low-frequency noises and false images generated from the conventional cross-correlation imaging condition, and further enhance the illumination in the VSP imaging. Synthetic examples and application to a walkaway field data demonstrate that it can attenuate the noise and improve the imaging resolution effectively. By using DDM with inverse scattering imaging condition and wavefield decomposition based on Hilbert transform, VSP imaging using free-surface–related multiples becomes a practical complement for conventional VSP imaging.