Full-Text Search:
Home|Journal Papers|About CNKI|User Service|FAQ|Contact Us|中文
《Chinese Journal of Geophysics》 2015-01
Add to Favorite Get Latest Update

Frequency multiscale full-waveform velocity inversion

ZHANG Wen-Sheng;LUO Jia;TENG Ji-Wen;Academy of Mathematics and Systems Science,Institute of Computational Mathematics and Scientific/Engineering Computing, State Key Laboratory of Scientific and Engineering Computing,Chinese Academy of Sciences;Institute of Geology and Geophysics,Chinese Academy of Sciences;  
Inverse problems arise naturally in oil geophysical exploration.In order to detect the structures of oil/gas objects,one of the most important step in seismic data processing is that we need to know the velocity information of subsurface.This belongs to scope of seismic inversion.There are various seismic inversion methods to achieve the goal,for example the seismic tomography based on ray tracing and the conventional velocity analysis based on CDP stack.Currently,the popular method is the full-waveform inversion(FWI)method based on wave equations.The FWI method uses the prestack data rather than the poststack data as an input and has high accuracy to image complicated velocity structures.The inversion can be implemented in the time domain or in the frequency domain.While there are many works of frequency-domain FWI methods,here we focus on the time-domain FWI method and try to inverse complicated velocity structures.Mathematically,the prestack data can be obtained by the forward modeling of wave equations.The seismic waves generated by a Ricker wavelet source propagate through the underground media.When the waves encounter the interfaces or inhomogeneous objects with different physical parameters,they are reflected back to the surface.These reflected waves are received by the geophones on the surface and used for inversion or other seismic data processing.This process is usually modeled by wave equations such as the acoustic wave equation if we assume the underground media is acoustic.In this paper,we investigate the full-waveform velocity inversion based on the 2-D acoustic wave equation and the corresponding computations are completed in the time domain.The forward problem is to solve the acoustic wave equation numerically and obtain the prestack data.Many numerical methods can be used to solve the forward problem,for example,the finite element method and the finite volume method and so on.In this paper we adopt the finite difference method because of easy programming and its high efficiency.Since the computational domain is a bounded rectangular domain and boundary reflections may devastate inversion results,it is necessary to use the absorbing boundary conditions and eliminate boundary reflections.The boundary conditions we adopted are the second-order Clayton′s absorbing boundary conditions.The computational discrete schemes for the equation and boundary conditions are of second-order accuracy both in time and space.The inversion is required to solve a nonlinear least-square problem.It is an iterative minimization process between the synthetic data and the observed data.Aiming at the difficulty which the waveform inversion is easy to fall into local extreme points,we propose the"stepwise inversion"strategy.First we use the wavefield on the lower frequency scale to obtain a reasonable initial model.Then we use the wavefield on the other high frequency scales to do inversion step by step.In order to use the wavefield information fully,the information at large scale contain that at low scale.Moreover,the inversion result at the previous low-frequency scale is chosen as the initial model for the next inversion at the large-frequency scale.The optimization method in inversion is the L-BFGS method.The full-waveform inversion is a typical large-scale scientific computational problem and the computations are implemented based on the MPI parallel algorithm on the PC cluster.The theoretical formulae and algorithms,including the finite-difference forward modeling,velocity model correction,gradient computation and corresponding algorithm,are expounded and derived in detail.In numerical computations,the full-waveform inversion based on the MPI parallel computations for the Marmousi model are completed.The Marmousi model is a complicated benchmark 2-D model which is usually used for testing the ability of migration and inversion methods.In the process of"stepwise inversion",the data at five different scales,i.e.0~5Hz,0~10 Hz,0~25 Hz,0~35 Hz,0~60 Hz,are selected.The initial model at first scale is a linear velocity model which is reasonable in practical case.The iteration inversion results at the first scale(i.e.0~5Hz)show that the inversion gets better as the iteration number increases till some iteration number such as the 50 th iteration.This phenomenon can be observed by the variations of object function.Then we stop the inversion and the corresponding inversion result is used for the initial velocity model for the inversion at the second cale(i.e.0~10Hz).The following procedure is similar till the inversion at the fifth scale(i.e.0~60Hz)is completed.At each scale,the inversion iteration number may be set 50 as more inversion iterations can not decline the value of object function.For comparisons,the inversion result without using the stepwise strategy is also given and shows that the iteration inversion converges to a wrong local extreme result.The detail comparisons between numerical and exact results at a fixed CDP are presented.The computations are completed on the PC cluster.The parallel efficiency is high and the scalability is about 0.9.The full waveform inversion is an iterative process of residual minimization between synthetic data and the known records.The inversion is easy to fall into local minimum points.We develop the sequential inversion method based on the inversion with the data at different frequency scales.The main idea is that the inversion result at low frequency scale is chosen as the initial guessing model for the inversion at the next high frequency scale.This strategy effectively solves the problem of inversion divergence when the initial value is far from the true solution.The detailed descriptions including theoretical formula and the corresponding algorithms are given or derived in this paper.Numerical computations for the complex structure model named Marmousi model are carried out.Relative good inversion results are yielded.Many computations show that the method is effective and has high robustness to the initial model.The full-waveform inversion based on wave equations is a typical large scale scientific computational problem. The implementation based on the MPI algorithm improves the computational efficiency greatly,which provides the basis for further application to real data.
【Fund】: 国家自然科学基金项目(11471328);; 973项目(2010CB731503)资助
【CateGory Index】: P315
Download(CAJ format) Download(PDF format)
CAJViewer7.0 supports all the CNKI file formats; AdobeReader only supports the PDF format.
【Citations】
Chinese Journal Full-text Database 2 Hits
1 YANG WEN-CAI(Institute of Geophysical & Geochemlcal Exploration Ministry of Geology and Mineral Resources, Langjang 102849);NONLINEAR CHAOTIC INVERSION OF SEISMIC TRACES:(I)THEORY AND NUMERICAL EXPERIMENTS[J];Chinese Journal of Geophysics;1993-02
2 YANG WEN-CAI(Institute of Geophysical and Ceochemical Exploration, Ministry of Geology and Mineral Resources, Langiang 102849);NONLINEAR CHAOTIC INVERSION OF SEISMIC TRACES: II. LYAPUNOV EXPONENTS AND ATTRACTORS[J];Chinese Journal of Geophysics;1993-03
【Co-citations】
Chinese Journal Full-text Database 10 Hits
1 Chen Zheng;Judgement on the existence of a common quadratic Lyapunov function for a kind of systems[J];Journal of Anhui Vocational College of Electronics & Information Technology;2006-04
2 SUN Xiao-jun,JIAO Jian-min,HE Jun-hong(Department of Mathematics,Baoji College of Arts and Science,Baoji 721013,China);A genetic accelerated trust region algorithm for solving optimization problems[J];Journal of Anhui University(Natural Sciences);2008-03
3 ZHANG Jia-xin1,2,DUAN Fu-jian1(1.School of Computational Science and Mathematics,Guilin University of Electronic Technology,Guilin 541004,China;2.School of Science,Anhui Science and Technology University,Fengyang 233100,China);A new method of SQP-filter with NCP function[J];Journal of Anhui University(Natural Sciences Edition);2010-05
4 GAO Hong1,MENG Shu2,ZHA Wei-min2,LI Ling-chun1,ZHANG Hai-tao1(1.Anhui Key Laboratory of Advanced Numerical Control & Servo Technology Anhui Ploytechnic University,Wuhu 241000,China; 2.Wuhu Bureau of Quality and Technical Supervision,Wuhu 241000,China);Configuration optimization of physical prototype for 6-3-3 parallel mechanism[J];Journal of Anhui Polytechnic University;2011-03
5 LI Bingfang,HU Shuiling(Department of Basic Course,Shaanxi Railway Institute,Weinan 714000,Shaanxi,China);Analysis on the Properties of Multi-Variables Convex Function[J];Journal of Ankang University;2011-03
6 YAN Feng,LI Shi-guo(School of Anshan University of Science and Technology,Anshan 114044,China);Design and optimization of automobile seats[J];Journal of Anshan University of Science and Technology;2004-05
7 LI Jiang-hai, SUN Qin (School of Aeronautics, Northwestern Polytechnical University, Xi'An 710072,China);Study of optimization algorithm for stealth performance of SRAM[J];Ordnance Material Science and Engineering;2004-02
8 CUI Peng ZHANG Li'ang JI Lijiu (School of Electronic Engineering & Computer Science,Peking University,Beijing,100871)GAO Li (School of Mathematical Sciences,Peking University,Beijing,100871);Comparisons of Two Smoothing Algorithms for Optimization in One-Dimensional Space[J];Acta Scicentiarum Naturalum Universitis Pekinesis;2003-05
9 XUE Yi College of Applied Sciences, Beijing Polytechnic University, Beijing 100022, China);Newton Method for Solving Non-negative Constraint Optimization Problems[J];Journal of Beijing Polytechnic University;2002-03
10 Ge Xinsheng (Division of Basic Courses, Beijing Institute of Machineny, Beijing 100085);OPTIMAL ATTITUDE CONTROL OF SPACECRAFT WITH STRETCHING SOLAR ARRAYS[J];JOURNAL OF BEIJING INSTITUTE OF MACHINERY;1998-04
China Proceedings of conference Full-text Database 10 Hits
1 LIU Yu,HAN Feng,DONG Nan,LU Xi-cheng,LEI Ming (Northwest Institute of Nuclear Technology,Xi'an 710024,China);The Confidence Interval Estimation of the Security Probability of Cylindrical Explosion Vessel[A];[C];2011
2 YAO Zhiying,LIU Dong 1.The Second Artillery Engineering College,Xi'an 710025,P.R.China E-mail:littlepigyao@163.com 2.Xi'an Institute of Surveying and Mapping,Xi'an 710054,P.R.China E-mail:ld_yzy@yahoo.com.cn;Iterative Filtering Based on Kazakov Linearization[A];[C];2010
3 YAO Hai-xiang(School of Informatics,Guangdong University of Foreign Studies,Guangzhou 510006,China);Portfolio Selection Based on Nonparametric Estimation and Quadric Utility Maximization Framework[A];[C];2011
4 YANG Xiao-chun,LI Xiao-fan,ZHANG Mei-gen (Institute of Geology & Geophysics,Chinese Academy of Science,Beijing 100101,China);Some Mathematical Prolems in Nonlinear Inversion of Seismic Wave[A];[C];2003
5 ZHANG Xue-mei~1,SUN Ruo-mei~1,YANG Hui~1,TENG Ji-wen~(1,2) (1.InsLiLuLe of Geology and Geophysics,Chinese Academy of Sciences,Beijing 100029,China; 2.Jilin University,Changchun 130026,China);Application of improved simulated annealing algorithm in the study of depth of interfaces in the crust and upper mantle[A];[C];2009
6 Zeng Rui.Liu Hong.Qin Yueshuang.Niu YanJiang. Institute of Geology and Geophysies.Chinese Academy of Sciences.Beijing 100029.China;Seismic inversion based on local wave decomposition and constrained interpolation.[A];[C];2009
7 LiYanPing, He Yu (School of Computer Science and Technology, ) (Beijing Technology and Business University, 100037);Finding a global extremum of a continuous function with VPM[A];[C];2006
8 Yi-Yong Pan~1 Ping-Qi Pan~1 1 Department of Mathematics,Southeast University,Nanjing 210096;Diagonal Second Order Quasi-Newton Method for Unconstrained Optimization[A];[C];2008
9 Zhao-Hui Jia Qin Ni Nanjing University of Aeronautics and Astronautics,Nanjing 210016;An Algorithm of Moving Asymptotes for Solving Unconstrained Optimization Problems[A];[C];2008
10 Tang Mingyun Academy of Mathematics and Systems Science,CAS,Beijing 100190;A Trust-Region-Newton Method for Unconstrained Optimization[A];[C];2008
【Secondary Citations】
Chinese Journal Full-text Database 6 Hits
1 Yang Wencai (Beijng Computer Center of Geology, Beijng, 100083);METHODS AND PROBLEMS IN SEISMIC INVERSION OF NONLINEAR WAVE EQUATIONS[J];Progress In Geophysics;1992-01
2 YANG WEN-CAI(Beijing Computer Center, Ministry of Geology and Mineral Resources);A GENERALIZED INVERSION TECHNIQUE FOR POTENTIAL FIELD DATA PROCESSING[J];Chinese Journal of Geophysics;1986-03
3 YANG WEN-CAI JIAO FU-GUANG(Beijing Computer Center, The Ministry of Geology and Mineral Resources);VELOCITY IMAGING FROM REFLECTION SEISMIC DATA BY JOINT INVERSION TECHNIQUES[J];Chinese Journal of Geophysics;1987-06
4 LIU SHI-DA(Department of Geophysics, Peking University, Beijing);EARTH SYSTEM MODELING AND CHAOTIC TIME SERIES[J];Chinese Journal of Geophysics;1990-02
5 YANG WEN-CAI(Institute of Geophysical & Geochemlcal Exploration Ministry of Geology and Mineral Resources, Langjang 102849);NONLINEAR CHAOTIC INVERSION OF SEISMIC TRACES:(I)THEORY AND NUMERICAL EXPERIMENTS[J];Chinese Journal of Geophysics;1993-02
6 YANG WEN-CAI(Institute of Geophysical and Ceochemical Exploration, Ministry of Geology and Mineral Resources, Langiang 102849);NONLINEAR CHAOTIC INVERSION OF SEISMIC TRACES: II. LYAPUNOV EXPONENTS AND ATTRACTORS[J];Chinese Journal of Geophysics;1993-03
©2006 Tsinghua Tongfang Knowledge Network Technology Co., Ltd.(Beijing)(TTKN) All rights reserved