Full-Waveform Inversion for Imaging Faulted Structures: A Case Study from the Japan Trench Forearc Slope

Full-waveform inversion (FWI) of limited-offset marine seismic data is a challenging task due to the lack of refracted energy and diving waves from the shallow sediments, which are fundamentally required to update the long-wavelength background velocity model in a tomographic fashion. When these events are absent, a reliable initial velocity model is necessary to ensure that the observed and simulated waveforms kinematically fit within an error of less than half a wavelength to protect the FWI iterative local optimization scheme from cycle skipping. We use a migration-based velocity analysis (MVA) method, including a combination of the layer-stripping approach and iterations of Kirchhoff prestack depth migration (KPSDM), to build an accurate initial velocity model for the FWI application on 2D seismic data with a maximum offset of 5.8 km. The data are acquired in the Japan Trench subduction zone, and we focus on the area where the shallow sediments overlying a highly reflective basement on top of the Cretaceous erosional unconformity are severely faulted and deformed. Despite the limited offsets available in the seismic data, our carefully designed workflow for data preconditioning, initial model building, and waveform inversion provides a velocity model that could improve the depth images down to almost 3.5 km. We present several quality control measures to assess the reliability of the resulting FWI model, including ray path illuminations, sensitivity kernels, reverse time migration (RTM) images, and KPSDM common image gathers. A direct comparison between the FWI and MVA velocity profiles reveals a sharp boundary at the Cretaceous basement interface, a feature that could not be observed in the MVA velocity model. The normal faults caused by the basal erosion of the upper plate in the study area reach the seafloor with evident subsidence of the shallow strata, implying that the faults are active.

Numerical Investigation of Wellbore Stability in Deepwater Shallow Sediments

An elaborate poro-elastoplastic numerical model has been developed in this paper to explore the stability characteristics of wellbore in shallow sediments of deepwater oil/gas wells. The combined Drucker-Prager/cap plasticity model is employed to characterize the mechanical behavior of the weakly consolidated or unconsolidated shallow sediments, by which both plastic compaction deformation and plastic shear deformation can be considered. Possible penetration of drilling fluid into the formation and its coupling to deformation have also been accounted for in the model. Using this model, deformation, stress evolution, and failure characteristics of the formation around the wellbore are analyzed in detail. Results presented in this paper demonstrate the necessity of considering the plastic compaction capability of the formation during the wellbore stability analysis of shallow sediments in deepwater. For mud pressures lower than the in situ horizontal stress, excessive wellbore shrinkage may occur if the mud pressure is too low, which, however, can be effectively mitigated through properly increasing the mud pressure even fluid penetration into the near-wellbore region may occur. It is also evidenced that, if penetration of drilling fluid into the formation is prevented, fracturing of the wellbore will not occur even the mud pressure is very high. Instead, the wellbore will expand substantially due to plastic compaction, and the deformed wellbore radius could be several times larger than the original value. However, if drilling fluid can penetrate into the formation, high pore pressure will develop within the near-wellbore region, resulting in tensile hoop stress at the wellbore and thus fracturing of the wellbore along the radial direction. The numerical results and implications in this paper are anticipated to be beneficial for the drilling operation in the shallow portion of deepwater oil/gas wells.

Full Waveform Inversion for Imaging Faulted Structures; a Case Study from Japan Trench Forearc Slope

Full waveform inversion (FWI) of limited-offset marine seismic data is a challenging task due to the lack of refracted energy and diving waves from the shallow sediments, which are fundamentally required to update the long-wavelength background velocity model through a tomographic fashion. When these events are absent, a reliable initial velocity model is necessary to assure that the observed and simulated waveforms kinematically fit within an error less than half a wavelength to protect the FWI iterative local optimization scheme from cycle skipping. We use a migration-based velocity analysis (MVA) method, including a combination of the layer stripping approach and iterations of Kirchhoff prestack depth migration (KPSDM) to build an accurate initial velocity model for the FWI application on 2D seismic data with a maximum offset of 5.8 km. The data is acquired in the Japan Trench subduction zone, and we focus on the area where the shallow sediments overlying a highly reflective basement on top of the Cretaceous erosional unconformity are severely faulted and deformed. Despite the limited offsets available in the seismic data, our carefully designed workflow for data preconditioning, initial model building, and waveform inversion provides a velocity model which could improve the depth images down to a depth of almost 3.5 km. We present several quality control measures to assess the reliability of the resulting FWI model, including ray path illuminations, sensitivity kernels, reverse time migration (RTM) images, and KPSDM common image gathers. A direct comparison between the FWI and MVA velocity profiles reveals a sharp boundary at the Cretaceous basement interface, a feature which could not be observed in the MVA velocity model. The normal faults caused by the basal erosion of the upper plate in the study area reach the seafloor with evident subsidence of the shallow strata, implying that the faults are active.

Temporal Temperature Distribution in Shallow Sediments of a Large Shallow Lake and Estimated Hyporheic Flux Using VFLUX 2 Model

Identifying and quantifying exchange flux across sediment-water interface is crucial when considering water and nutrient contributions to a eutrophic lake. In this study, observed temporal temperature distributions in shallow sediment of Lake Taihu (Eastern China) based on three-depth sensors at 14 sites throughout 2016 were used to assess temporal water exchange patterns. Results show that temporal temperature in shallow sediments differed with sampling sites and depths and the temperature amplitudes also clearly shrunk as the offshore distance increasing. Exchange fluxes estimated using the VFLUX 2 model based on temperature amplitude show that alternating-direction temporal flow exists in the eastern zone of Lake Taihu with averages of −13.0, −0.6, and 3.4 mm day−1 (negative represents discharging into the lake) at three nearshore sites (0.5, 2.0, and 6.0 km away from the shoreline, respectively). Whereas downwelling flow occurred throughout almost the entire year with averages of 37.7, 23.5, and 6.6 mm day−1 at the three southern nearshore sites, respectively. However, upwelling flow occurred throughout almost the entire year and varied widely in the western zone with averages of −74.8, 45.9, and −27.0 mm day–1 and in the northern zone with averages of −76.2, −55.3, and −51.1 mm day−1. The estimated fluxes in the central zone were relatively low and varied slightly during the entire year (−15.1 to 22.5 mm day−1 with an average of −0.7 mm day−1). Compared with the sub sensor pair (at 5 and 10 cm), the estimated hyporheic fluxes based on the top sensor pair (at 0 and 5 cm) varied within wider ranges and exhibited relatively larger values. Effects of upwelling flow at the western and northern zones need to be paid attention to on nearshore water quality particularly during winter and spring seasons. Estimated flow patterns at the four zones summarily reflect the seasonal water interaction near the sediment surface of Lake Taihu and are beneficial to improve its comprehensive management. Thermal dispersivity usually used for estimating the thermal diffusivity is more sensitive for upward hyporheic flux estimating even if with a low flux. Temperature amplitude ratio method can be used to estimate the exchange flux and suitable for low flux conditions (either upwelling or downwelling). A better evaluation of the exchange flux near inclined nearshore zones might need an optimized installation of temperature sensors along with the potential flow path and/or a vertical two-dimensional model in the future.