A reality check on full-wave inversion applied to land seismic data for near-surface modeling

We evaluate the performance of traveltime tomography and full-wave inversion (FWI) for near-surface modeling using the data from a shallow seismic field experiment. Eight boreholes up to 20-m depth have been drilled along the seismic line traverse to verify the accuracy of the P-wave velocity-depth model estimated by seismic inversion. The velocity-depth model of the soil column estimated by traveltime tomography is in good agreement with the borehole data. We used the traveltime tomography model as an initial model and performed FWI. Full-wave acoustic and elastic inversions, however, have failed to converge to a velocity-depth model that desirably should be a high-resolution version of the model estimated by traveltime tomography. Moreover, there are significant discrepancies between the estimated models and the borehole data. It is understandable why full-wave acoustic inversion would fail — land seismic data inherently are elastic wavefields. The question is: Why does full-wave elastic inversion also fail? The strategy to prevent full-wave elastic inversion of vertical-component geophone data trapped in a local minimum that results in a physically implausible near-surface model may be cascaded inversion. Specifically, we perform traveltime tomography to estimate a P-wave velocity-depth model for the near-surface and Rayleigh-wave inversion to estimate an S-wave velocity-depth model for the near-surface, then use the resulting pairs of models as the initial models for the subsequent full-wave elastic inversion. Nonetheless, as demonstrated by the field data example here, the elastic-wave inversion yields a near-surface solution that still is not in agreement with the borehole data. Here, we investigate the limitations of FWI applied to land seismic data for near-surface modeling.

Multidirectional Anisotropic Total Variation and Its Application in the Tomography of the Surrounding Rock of Coal Mining Faces

The computed tomography (CT) reconstruction algorithm is one of the crucial components of the CT system. To date, total variation (TV) has been widely used in CT reconstruction algorithms. Although TV utilizes the a priori information of the longitudinal and lateral gradient sparsity of an image, it introduces some staircase artifacts. To overcome the current limitations of TV and improve imaging quality, we propose a multidirectional anisotropic total variation (MATV) that uses multidirectional gradient information. The surrounding rock of coal mining faces uses principles of tomography similar to those of medical X-rays. The velocity distribution for the surrounding rock can be obtained by the first-arrival traveltime tomography of the transmitted waves in the coal mining face. Combined with the geological data, we can interpret the geological hazards in the coal mining face. To perform traveltime tomography, we first established the objective function of the first-arrival traveltime tomography of the transmitted waves based on the MATV regularization and then used the split Bregman method to solve the objective function. The simulated data and real data show that the MATV regularization method proposed in this paper can better maintain the boundaries of geological anomalies and reduce the artifacts compared with the isotropic total variation regularization method and the anisotropic total variation regularization method. Furthermore, this approach describes the distribution of geological anomalies more accurately and effectively and improves imaging accuracy.

Ray illumination compensation for adjoint-state first-arrival traveltime tomography

First-arrival traveltime tomography is an essential method for obtaining near-surface velocity models. The adjoint-state first-arrival traveltime tomography is appealing due to its straightforward implementation, low computational cost, and low memory consumption. Because solving the point-source isotropic eikonal equation by either ray tracers or eikonal solvers intrinsically corresponds to emanating discrete rays from the source point, the resulting traveltime gradient is singular at the source point, and we denote such a singular pattern the imprint of ray illumination. Because the adjoint-state equation propagates traveltime residuals back to the source point according to the negative traveltime gradient, the resulting adjoint state will inherit such an imprint of ray illumination, leading to singular gradient descent directions when updating the velocity model in the adjoint-state traveltime tomography. To mitigate this imprint, we propose to solve the adjoint-state equation twice but with different boundary conditions: one being taken to be regular data residuals, and the other taken to be ones uniformly, so that we are able to use the latter adjoint state to normalize the regular one and we further use the normalized quantity to serve as the gradient direction to update the velocity model; we call this process the ray-illumination compensation. To overcome the issue of limited aperture, we propose a spatially varying regularization method to stabilize the new gradient direction. A synthetic example demonstrates that the proposed method is able to mitigate the imprint of ray illumination, remove the footprint effect near source points, and provide uniform velocity updates along ray paths. A complex example extracted from the Marmousi2 model and a migration example illustrate that the new method accurately recovers the velocity model, and an offset-dependent inversion strategy can further improve the quality of recovered velocity models.

Adjoint traveltime tomography unravels a scenario of horizontal mantle flow beneath the North China craton

The North China craton (NCC) was dominated by tectonic extension from late Cretaceous to Cenozoic, yet seismic studies on the relationship between crust extension and lithospheric mantle deformation are scarce. Here we present a three dimensional radially anisotropic model of NCC derived from adjoint traveltime tomography to address this issue. We find a prominent low S-wave velocity anomaly at lithospheric mantle depths beneath the Taihang Mountains, which extends eastward with a gradually decreasing amplitude. The horizontally elongated low-velocity anomaly is also featured by a distinctive positive radial anisotropy (VSH > VSV). Combining geodetic and other seismic measurements, we speculate the presence of a horizontal mantle flow beneath central and eastern NCC, which led to the extension of the overlying crust. We suggest that the rollback of Western Pacific slab likely played a pivotal role in generating the horizontal mantle flow at lithospheric depth beneath the central and eastern NCC.

Imaging structure and geometry of slabs in the greater Alpine area – A P-wave traveltime tomography using AlpArray Seismic Network data

We perform a teleseismic P-wave traveltime tomography to examine the geometry and structure of subducted lithosphere in the upper mantle beneath the Alpine orogen. The tomography is based on waveforms recorded at over 600 temporary and permanent broadband stations of the dense AlpArray Seismic Network deployed by 24 different European institutions in the greater Alpine region, reaching from the Massif Central to the Pannonian Basin and from the Po plain to the river Main. Teleseismic traveltimes and traveltime residuals of direct teleseismic P-waves from 331 teleseismic events of magnitude 5.5 and higher recorded between 2015 and 2019 by the AlpArray Seismic Network are extracted from the recorded waveforms using a combination of automatic picking, beamforming and cross-correlation. The resulting database contains over 162.000 highly accurate absolute P-wave traveltimes and traveltime residuals. For tomographic inversion, we define a model domain encompassing the entire Alpine region down to a depth of 600 km. Outside this domain, a laterally homogeneous standard earth model is assumed. Predictions of traveltimes are computed in a hybrid way applying a fast Tau-P method outside the model domain and continuing the wavefronts into the model domain using a fast marching method. For teleseismic inversion, we iteratively invert demeaned traveltime residuals for P-wave velocities in the model domain using a regular discretization with an average lateral spacing of about 25 km and a vertical spacing of 15 km. The inversion is regularized towards an initial model constructed from an a priori model of the crust and uppermost mantle and a standard earth model beneath. The resulting model provides a detailed image of slab configuration beneath the Alpine and Apenninic orogens. Major features are an overturned Adriatic slab beneath the Apennines reaching down to 400 km depth still attached in its northern part to the crust but exhibiting detachment towards the southeast. A fast anomaly beneath the western Alps indicates a short western Alpine slab that ends at about 100 km depth close to the Penninic front. Further to the east and following the arcuate shape of the western Periadriatic Fault System, a deep-reaching coherent fast anomaly with complex interior stucture generally dipping to the SE down to about 400 km suggests a slab of European origin extending eastward to the Giudicarie fault. This slab is detached from overlying lithosphere at its eastern end below a depth of about 100 km. Further to the east, well-separated from the slab beneath the western and central Alps, another deep-reaching, nearly vertically dipping high-velocity anomaly suggests the existence of a slab beneath the Eastern Alps of presumably European origin which is completely detached from the orogenic root. Our image of this slab does not require a polarity switch because of its nearly vertical dip and full detachment from the overlying lithosphere. Fast anomalies beneath the Dinarides are weak and concentrated to the northernmost part and shallow depths. Low-velocity regions surrounding the fast anomalies beneath the Alps to the west and northwest follow the same dipping trend as the overlying fast ones, indicating a kinematically coherent subducting tectosphere in this region. In contrast, low-velocity anomalies to the east suggest asthenospheric upwelling presumably driven by retreat of the Carpathian slab and extrusion of eastern Alpine lithosphere towards the east while low velocities to the south are presumably evidence of asthenospheric upwelling and mantle hydration due to the backarc position behind the European slab.

Low-frequency, long-offset elastic waveform inversion in the context of velocity model building

Classic imaging approaches consist of splitting the earth into background and reflectivity models. When justified, this separation of scale is quite powerful, although this approach relies on some smoothness and weak contrast assumptions. This approach allows for the imaging methods to be based on acoustic wave propagation after having identified the compressional waves through picking or signal processing. Over the past years, wave-equation tomography and waveform inversion approaches have become routine, complementing the classic approaches to derive background models. They do not rely on high-frequency picks, unlike ray-based traveltime tomography, but on low-frequency cross-correlation to define time shifts and on waveform matching. In the presence of large earth parameter contrasts, time shifts and waveforms of compressional waves may depend on elastic parameters when interferences occur within the Fresnel zones. This challenges the recovery of the background model under an acoustic assumption with low-frequency data. Accounting for an elastic propagation in waveform inversion, even in the context of model building, could help to reduce the artifacts seen in acoustic results. A synthetic and a real data example are presented to illustrate the potential benefit of using an elastic waveform inversion approach when inverting long-offset, low-frequency seismic data.