Full-Waveform Inversion Using a Learned Regularization

被引:7
作者
Sun, Pengpeng [1 ]
Yang, Fangshu [2 ]
Liang, Hongxian [1 ]
Ma, Jianwei [3 ]
机构
[1] SINOPEC Shengli Oilfield Co, Geophys Res Inst, Dongying 257022, Shandong, Peoples R China
[2] Harbin Inst Technol, Dept Math, Weihai 264209, Peoples R China
[3] Peking Univ, Sch Earth & Space Sci, Beijing 100871, Peoples R China
来源
IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING | 2023年 / 61卷
基金
中国国家自然科学基金;
关键词
Mathematical models; Data models; Optimization; Convolutional neural networks; Inverse problems; Propagation; Imaging; Deep convolutional neural network (CNN) denoiser; full-waveform inversion (FWI); learned regularization; transfer learning; NEURAL-NETWORKS; PLAY PRIORS; ALGORITHM; GRADIENT;
D O I
10.1109/TGRS.2023.3322964
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Full-waveform inversion (FWI) is an efficient technique for capturing the subsurface physical features by iteratively minimizing the misfit between simulated and observed seismograms. As such a problem is ill-posed, a significant ingredient for a satisfactory solution is to incorporate desirable priors. Most traditional regularized FWI approaches suffer from the nonadaptiveness and sensitivity of regularizers. Deep-learning-assisted inversion methods can use pretrained priors or parametric network priors as regularizers. We develop a novel FWI method based on a physics-constrained iterative algorithm with a learned regularization (FWIPLR). The introduced framework easily enables high-quality inversion by integrating two complementary terms: the physical constraints of the imaging system characterized by its forward model and a priori knowledge of the expected results characterized by a deep convolutional neural network (CNN)-based regularizer. In particular, the advanced CNN denoiser, which corresponds to an implicit regularization term, is first trained with large-scale natural images and then fine-tuned with small-scale geological images. Such a transfer learning strategy seems appealing as it makes FWIPLR more generic to distinct geological models. To stabilize the optimization, we use the spectral normalization instead of batch normalization to impose Lipschitz constraint on the networks. We validate the method effectiveness on Marmousi, overthrust, and 2004 BP models. The comparable results illustrate the ability of FWIPLR for producing high-resolution subsurface structures compared with total variation (TV)-regularized FWI, conventional denoiser, and other CNN denoiser-regularized FWI; in particular, when the initial model is inaccurate, data lose low-frequency information or model contains high-contrast media.
引用
收藏
页数:15
相关论文
共 71 条
[1]   Deep Learning for Seismic Inverse Problems: Toward the Acceleration of Geophysical Analysis Workflows [J].
Adler, Amir ;
Araya-Polo, Mauricio ;
Poggio, Tomaso .
IEEE SIGNAL PROCESSING MAGAZINE, 2021, 38 (02) :89-119
[2]   Full waveform inversion by proximal Newton method using adaptive regularization [J].
Aghamiry, H. S. ;
Gholami, A. ;
Operto, S. .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2021, 224 (01) :169-180
[3]   NTIRE 2017 Challenge on Single Image Super-Resolution: Dataset and Study [J].
Agustsson, Eirikur ;
Timofte, Radu .
2017 IEEE CONFERENCE ON COMPUTER VISION AND PATTERN RECOGNITION WORKSHOPS (CVPRW), 2017, :1122-1131
[4]  
Aminzadeh F., 1994, The Leading Edge, V13, P949, DOI DOI 10.1190/1.1437054
[5]  
[Anonymous], 2017, 87 ANN INT M SEG, DOI DOI 10.1190/SEGAM2017-17627643.1
[6]  
[Anonymous], 2005, P 67 EAGE C EXH
[7]  
[Anonymous], 2004, Finite Difference Schemes and Partial Differential Equations
[8]  
[Anonymous], 1994, The Leading Edge, DOI [DOI 10.1190/1.1437051, 10.1190/1.1437051]
[9]   Deep learning-driven velocity model building workflow [J].
Araya-Polo M. ;
Farris S. ;
Florez M. .
Leading Edge, 2019, 38 (11) :872A1-872A9
[10]   Regularized seismic full waveform inversion with prior model information [J].
Asnaashari, Amir ;
Brossier, Romain ;
Garambois, Stephane ;
Audebert, Francois ;
Thore, Pierre ;
Virieux, Jean .
GEOPHYSICS, 2013, 78 (02) :R25-R36