Multiparameter full-waveform inversion of 3D ocean-bottom cable data from the Valhall field

被引:0
作者
Kamath, Nishant [1 ,2 ]
Brossier, Romain [1 ]
Metivier, Ludovic [1 ,3 ]
Pladys, Arnaud [1 ]
Yang, Pengliang [1 ]
机构
[1] Univ Grenoble Alpes, ISTerre, F-38400 Grenoble, France
[2] CGG, Crawley, England
[3] Univ Grenoble Alpes, CNRS, LJK, F-38041 Grenoble, France
关键词
FREQUENCY-DOMAIN; MISFIT FUNCTIONS; DIRECT SOLVER; PART; EFFICIENT; REVERSE; DENSITY; ATTENUATION; PROPAGATION; ALGORITHM;
D O I
10.1190/GEO2019-0705.1
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Full-waveform inversion (FWI) applications on 3D oceanbottom cable (OBC) data from the Valhall oil field in the North Sea have demonstrated the importance of appropriately accounting for attenuation. The Valhall field contains unconsolidated shallow sediments and a low-velocity anomaly in its center - indicative of gas clouds - which have a significant attenuation imprint on the data. Our challenge is to perform time-domain viscoacoustic 3D FWI, which requires more sophisticated tools than in the frequency domain wherein attenuation can be incorporated in a straightforward manner. The benefit of using a viscoacoustic, instead of a purely acoustic, modeling engine is illustrated. We have determined that, in the frequency band used (2.5-7.0 Hz), it is better to reconstruct the velocity only keeping the attenuation lixed because simultaneous inversion of the velocity and quality factor Q does not provide reliable Q updates. We develop an efficient time-domain workflow combining a random source decimation algorithm, modeling using standard linear solid mechanisms, and wavefield preconditioning. Our results are similar to those obtained from state-of-the-art frequency-domain algorithms, at a lower computational cost compared to conventional check-pointing techniques. We clearly illustrate the improvement in terms of imaging and data fit achieved when accounting for attenuation.
引用
收藏
页码:B15 / B35
页数:21
相关论文
共 50 条
[1]   Acoustic approximations for processing in transversely isotropic media [J].
Alkhalifah, T .
GEOPHYSICS, 1998, 63 (02) :623-631
[2]  
Amestoy P., 2015, 85 ANN INT M SEG, P3606, DOI [10.1190/segam2015-5811693.1, DOI 10.1190/SEGAM2015-5811693.1]
[3]   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
[4]  
Anderson JE, 2012, GEOPHYSICS, V77, pS93, DOI [10.1190/geo2011-0114.1, 10.1190/GEO2011-0114.1]
[5]   Time-lapse seismic imaging using regularized full-waveform inversion with a prior model: which strategy? [J].
Asnaashari, A. ;
Brossier, R. ;
Garambois, S. ;
Audebert, F. ;
Thore, P. ;
Virieux, J. .
GEOPHYSICAL PROSPECTING, 2015, 63 (01) :78-98
[6]  
Barkved O., 2003, 65 ANN INT C EXH EAG
[7]   MODELING OF A CONSTANT-Q - METHODOLOGY AND ALGORITHM FOR AN EFFICIENT AND OPTIMALLY INEXPENSIVE VISCOELASTIC TECHNIQUE [J].
BLANCH, JO ;
ROBERTSSON, JOA ;
SYMES, WW .
GEOPHYSICS, 1995, 60 (01) :176-184
[8]   Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion [J].
Brossier, Romain ;
Operto, Stephane ;
Virieux, Jean .
GEOPHYSICS, 2009, 74 (06) :WCC105-WCC118
[9]  
Casula G., 1992, Bollettino di Geofisica Teorica ed Applicada, V34, P235
[10]   A NONREFLECTING BOUNDARY-CONDITION FOR DISCRETE ACOUSTIC AND ELASTIC WAVE-EQUATIONS [J].
CERJAN, C ;
KOSLOFF, D ;
KOSLOFF, R ;
RESHEF, M .
GEOPHYSICS, 1985, 50 (04) :705-708