B uilding high-resolution models of several physical properties of the subsurface by multiparameter full waveform inversion (FWI) of multicomponent data will be a challenge for seismic imaging for the next decade. The physical properties, which govern propagation of seismic waves in visco-elastic media, are the velocities of the P-and S-waves, density, attenuation, and anisotropic parameters. Updating each property is challenging because several parameters of a different nature can have a coupled effect on the seismic response for a particular propagation regime (from transmission to reflection). This is generally referred to as trade-off or crosstalk between parameters. Moreover, different parameter classes can have different orders of magnitude or physical units and footprints of different strength in the wavefield, which can make the inversion poorly conditioned if it is not properly scaled. These difficulties raise the issue of a suitable parameterization for multiparameter FWI, where the term parameterization must be understood as a set of independent parameter classes that fully describe the subsurface properties. Many combinations of parameters can be viewed and this choice is not neutral as the parameterization controls the trade-off between parameters and the local resolution with which they can be reconstructed. Once this parameterization is selected, the subset of parameter classes in the parameterization, that can be reliably updated during the inversion, must be identified to avoid overparameterization of the optimization problem. The purpose of this tutorial is to provide a comprehensive overview of the promise, pitfalls, and open questions underlying multiparameter FWI. We first review the main FWI ingredient that controls the tradeoff between parameters, namely the radiation pattern of the so-called virtual sources, and some tools for analyzing these trade-offs. Then, we present some illustrative examples of multiparameter FWI, which should provide some guidelines to choose a suitable parameterization for FWI in visco-acoustic anisotropic media. We conclude by proposing a data-driven and model-driven workflow for visco-elastic anisotropic FWI of multicomponent marine data, which has been inspired by a real data case study from Valhall.

[1] Splay faults within accretionary complexes are commonly associated with the updip limit of the seismogenic zone. Prestack depth migration of a multichannel seismic line across the north Ecuador-south Colombia oceanic margin images a crustal splay fault that correlates with the seaward limit of the rupture zone of the 1958 (Mw 7.7) tsunamogenic subduction earthquake. The splay fault separates 5-6.6 km/s velocity, inner wedge basement rocks, which belong to the accreted Gorgona oceanic terrane, from 4 to 5 km/s velocity outer wedge rocks. The outer wedge is dominated by basal tectonic erosion. Despite a 3-km-thick trench fill, subduction of 2-km-high seamount prevented tectonic accretion and promotes basal tectonic erosion. The low-velocity and poorly reflective subduction channel that underlies the outer wedge is associated with the aseismic, décollement thrust. Subduction channel fluids are expected to migrate upward along splay faults and alter outer wedge rocks. Conversely, duplexes are interpreted to form from and above subducting sediment, at $14-to 15-km depths between the overlapping seismogenic part of the splay fault and the underlying aseismic décollement. Coeval basal erosion of the outer wedge and underplating beneath the apex of inner wedge control the margin mass budget, which comes out negative. Intraoceanic basement fossil listric normal faults and a rift zone inverted in a flower structure reflect the evolution of the Gorgona terrane from Cretaceous extension to likely Eocene oblique compression. The splay faults could have resulted from tectonic inversion of listric normal faults, thus showing how inherited structures may promote fluid flow across margin basement and control seismogenesis.

International audienceRifting between Arabia and Somalia started around 35 Ma followed by spreading at 17.6 Ma in the eastern part of the Gulf of Aden. The first-order segment between Alula-Fartak and Socotra-Hadbeen fracture zones is divided into three second-order segments with different structure and morphology. Seismic reflection data were collected during the Encens Cruise in 2006 on the northeastern margin. In this study, we present the results of Pre-Stack Depth Migration of the multichannel seismic data from the western segment, which allows us to propose a tectono-stratigraphic model of the evolution of this segment of the margin from rifting to the present day. The chronological interpretation of the sedimentary sequences is mapped out within relation to the onshore observations and existing dating. After a major development of syn-rift grabens and horsts, the deformation localized where the crust is the thinnest. This deformation occurred in the distal margin graben (DIM) at the northern boundary of the ocean-continent transition (OCT) represented by the OCT ridge. At the onset of the OCT formation differential uplift induced a submarine landslide on top of the deepest tilted block and the crustal deformation was restricted to the southern part of the DIM graben, where the continental break-up finally occurred. Initial seafloor spreading was followed by post-rift magmatic events (flows, sills and volcano-sedimentary wedge), whose timing is constrained by the analysis of the sedimentary cover of the OCT ridge, correlated with onshore stratigraphy. The OCT ridge may represent exhumed serpentinized mantle intruded by post-rift magmatic material, which modified the OCT after its emplacement