Search

We present an integrated algorithm for high‐resolution crustal‐scale imaging utilizing long‐offset wide‐angle ocean‐bottom seismometer (OBS) data and short‐offset multichannel streamer (MCS) data. The algorithm adopts a two‐step imaging strategy, initially using the OBS data to enhance deep structure imaging and capture long‐wavelength features of the migration velocity. Subsequently, the MCS data are migrated to recover detailed short‐wavelength components and shallow structures using the migration velocity model modified by the OBS result. Both steps employ the least‐squares reverse time migration (LSRTM) method with appropriate regularization techniques. The algorithm is implemented using the Bregmanized operator splitting (BOS) approach, known for its efficiency and adaptability in handling non‐smooth regularization, such as total‐variation (TV) constraints. To improve computational efficiency, compressed sensing is employed to store incident wavefields at half their actual size and reconstruct them accurately when required. Convergence is expedited through the use of preconditioners and the Anderson acceleration method. The proposed algorithm is validated through large‐scale numerical examples, demonstrating its robustness even with an initially imprecise velocity model. Results showcase improved resolution and accuracy achieved through the integration of OBS and MCS data. This study offers a comprehensive framework for crustal‐scale imaging, addresses computational complexities, and provides enhanced imaging capabilities. Plain Language Summary: We present a method for high‐resolution crustal‐scale imaging using a combination of long‐offset wide‐angle data from stationary receivers deployed on the seafloor with short‐offset data from floating receivers. Our method is structured as a two‐step imaging algorithm that exploits the unique strengths of each data set. In the first step, the stationary recorded data is used to recover the long‐wavelength features. Subsequently, we integrate the floating receiver data with the preceding image to recover the short‐wavelength components, adding details. The efficacy of our method relies on the utilization of least‐squares reverse time migration, which is enhanced by incorporating appropriate regularization techniques. This enables us to address critical challenges such as memory overhead, convergence speed, and image stability. These key challenges are addressed by incorporating three advanced strategies: compressed sensing theory, Anderson acceleration, and total variation regularization. Accordingly, we enhance computational efficiency, expedite convergence, and elevate imaging quality, which helps provide accurate, high‐resolution images of the subsurface that facilitate better geological interpretation of deep structures and therefore improve our understanding of the target's regional geodynamic context. Moreover, this algorithm allows imaging with less initial knowledge about the image and more robustness to the frequency range. Key Points: Enhanced crustal‐scale seismic imaging from academic dataOcean bottom seismometer and multi‐channel streamer data integrationHigh‐resolution least‐squares reverse time migration [ABSTRACT FROM AUTHOR]

Imaging via pre-stack depth migration (PSDM) of reflection towed-streamer multichannel seismic (MCS) data at the scale of the whole crust is inherently difficult. This is because the depth penetration of the seismic wavefield is controlled, firstly, by the acquisition design, such as streamer length and air-gun source configuration, and secondly by the complexity of the crustal structure. Indeed, the limited length of the streamer makes the estimation of velocities from deep targets challenging due to the velocity–depth ambiguity. This problem is even more pronounced when processing 2-D seismic data due to the lack of multi-azimuthal coverage. Therefore, in order to broaden our knowledge about the deep crust using seismic methods, we present the development of specific imaging workflows that integrate different seismic data. Here we propose the combination of velocity model building using (i) first-arrival tomography (FAT) and full-waveform inversion (FWI) of wide-angle, long-offset data collected by stationary ocean-bottom seismometers (OBSs) and (ii) PSDM of short-spread towed-streamer MCS data for reflectivity imaging, with the former velocity model as a background model. We present an application of such a workflow to seismic data collected by the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) and the Institut Français de Recherche pour l'Exploitation de la Mer (IFREMER) in the eastern Nankai Trough (Tokai area) during the 2000–2001 Seize France Japan (SFJ) experiment. We show that the FWI model, although derived from OBS data, provides an acceptable background velocity field for the PSDM of the MCS data. From the initial PSDM, we refine the FWI background velocity model by minimizing the residual move-outs (RMOs) picked in the pre-stack-migrated volume through slope tomography (ST), from which we generate a better-focused migrated image. Such integration of different seismic datasets and leading-edge imaging techniques led to greatly improved imaging at different scales. That is, large to intermediate crustal units identified in the high-resolution FWI velocity model extensively complement the short-wavelength reflectivity inferred from the MCS data to better constrain the structural factors controlling the geodynamics of the Nankai Trough. [ABSTRACT FROM AUTHOR]