On the role of density and attenuation in three-dimensional multiparameter viscoacoustic VTI frequency-domain FWI: an OBC case study from the North Sea

被引:59
作者
Operto, S. [1 ]
Miniussi, A. [2 ]
机构
[1] Univ Cote Azur, CNRS, IRD, OCA, F-06905 Geoazur, Sophia Antipoli, France
[2] Univ Cote Azur, OCA, F-06394 Nice 4, France
关键词
Body waves; Computational seismology; Controlled source seismology; Inverse theory; Seismic anisotropy; Wave propagation; WAVE-FORM INVERSION; BOTTOM-CABLE DATA; PART; FINITE-DIFFERENCE; DIRECT SOLVER; DIFFRACTION TOMOGRAPHY; VELOCITY; VALHALL; PROPAGATION; SENSITIVITY;
D O I
10.1093/gji/ggy103
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
3-D frequency-domain full waveform inversion (FWI) is applied on North Sea wide-azimuth ocean-bottom cable data at low frequencies (<= 10 Hz) to jointly update vertical wave speed, density and quality factor Q in the viscoacoustic VTI approximation. We assess whether density and Q should be viewed as proxy to absorb artefacts resulting from approximate wave physics or are valuable for interpretation in the presence of soft sediments and gas cloud. FWI is performed in the frequency domain to account for attenuation easily. Multiparameter frequency-domain FWI is efficiently performed with a few discrete frequencies following a multiscale frequency continuation. However, grouping a few frequencies during each multiscale step is necessary to mitigate acquisition footprint and match dispersive shallow guided waves. Q and density absorb a significant part of the acquisition footprint hence cleaning the velocity model from this pollution. Low Q perturbations correlate with low-velocity zones associated with soft sediments and gas cloud. However, the amplitudes of the Q perturbations show significant variations when the inversion tuning is modified. This dispersion in the Q reconstructions is however not passed on the velocity parameter suggesting that cross-talks between first-order kinematic and second-order dynamic parameters are limited. The density model shows a good match with a well log at shallow depths. Moreover, the impedance built a posteriori from the FWI velocity and density models shows a well-focused image with however local differences with the velocity model near the sea bed where density might have absorbed elastic effects. The FWI models are finally assessed against time-domain synthetic seismogram modelling performed with the same frequency-domain modelling engine used for FWI.
引用
收藏
页码:2037 / 2059
页数:23
相关论文
共 55 条
[1]  
Aki K., 2002, QUANTITATIVE SEISMOL
[2]   Fast 3D frequency-domain full-waveform inversion with a parallel block low-rank multifrontal direct solver: Application to OBC data from the North Sea [J].
Amestoy, Patrick ;
Brossier, Romain ;
Buttari, Alfredo ;
L'Excellent, Jean-Yves ;
Mary, Theo ;
Metivier, Ludovic ;
Miniussi, Alain ;
Operto, Stephane .
GEOPHYSICS, 2016, 81 (06) :R363-R383
[3]   Full waveform inversion for seismic velocity and anelastic losses in heterogeneous structures [J].
Askan, Aysegul ;
Akcelik, Volkan ;
Bielak, Jacobo ;
Ghattas, Omar .
BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA, 2007, 97 (06) :1990-2008
[4]  
Bai J., 2014, SOC EXPLORATION GEOP, P922
[5]  
Barkved O., 2010, 80th SEG Annual International Meeting, Expanded Abstracts, P925
[6]  
Boiero D., 2013, The Leading Edge, V32, P638, DOI DOI 10.1190/TLE32060638.1
[7]  
Boiero D., 2014, 76 EAGE C EXH AMST W
[8]  
Borisov D., 2014, 76 ANN EAGE M AMST
[9]   Fast full waveform inversion with source encoding and second-order optimization methods [J].
Castellanos, Clara ;
Metivier, Ludovic ;
Operto, Stephane ;
Brossier, Romain ;
Virieux, Jean .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2015, 200 (02) :718-742
[10]  
Dupuy B., 2016, LEADING EDGE, V35, P1146