As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing i...As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues： first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.展开更多
The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the tim...The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.展开更多
Frequency-domain waveform seismic tomography includes modeling of wave propagation and full waveform inversion of correcting the initial velocity model. In the forward modeling, we use direct solution based on sparse ...Frequency-domain waveform seismic tomography includes modeling of wave propagation and full waveform inversion of correcting the initial velocity model. In the forward modeling, we use direct solution based on sparse matrix factorization, combined with nine-point finite-difference for the linear system of equations. In the waveform inversion, we use preconditioned gradient method where the preconditioner is provided by the diagonal of the approximate Hessian matrix. We successfully applied waveform inversion method from low to high frequency in two sets of Marmousi data. One is the data set generated by frequencydomain finite-difference modeling, and the other is the original Marmousi shots data set. The former result is very close to the true velocity model. In the original shots data set inversion, we replace the prior source with estimated source; the result is also acceptable, and consistent with the true model.展开更多
基于常规模拟退火算法的零偏VSP全波形反演面临着计算量大和耗时长的问题。为此提出了一种不同阶段对应不同扰动模型和退火方式的分段快速模拟退火(segmented fast simulated annealing,SFSA)反演策略,以提高零偏VSP资料全波形反演的效...基于常规模拟退火算法的零偏VSP全波形反演面临着计算量大和耗时长的问题。为此提出了一种不同阶段对应不同扰动模型和退火方式的分段快速模拟退火(segmented fast simulated annealing,SFSA)反演策略,以提高零偏VSP资料全波形反演的效率。在反演前期采用大模型扰动空间和较慢温度衰减速度,充分发挥全局搜索能力,而在后期引入限制因子产生扰动模型,在迭代不断增加的时候逐渐减小模型的扰动空间,同时采用较快的温度衰减速度,有效提高反演的速度,使反演快速收敛到最优解。基于相同的初始温度和马尔可夫链长度,分别利用基于SFSA和非常快速模拟退火(very fast simulated annealing,VFSA)方法进行零偏VSP纵波速度全波形反演测试。结果表明,基于SFSA的反演方法的反演效率提高约50%,在迭代次数更少的条件下能获得更好的反演效果。基于SFSA的零偏VSP全波形反演具有高效和高精度的特点,其反演结果为地震地质层位标定、成果解释及油气预测奠定了基础。展开更多
基金financially supported by the National Natural Science Foundation of China(No.41074075/D0409)the National Science and Technology Major Project(No.2011ZX05025-001-04)
文摘As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues： first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.
基金Supported by the National Natural Science Foundation of China (69983005)
文摘Frequency-domain waveform seismic tomography includes modeling of wave propagation and full waveform inversion of correcting the initial velocity model. In the forward modeling, we use direct solution based on sparse matrix factorization, combined with nine-point finite-difference for the linear system of equations. In the waveform inversion, we use preconditioned gradient method where the preconditioner is provided by the diagonal of the approximate Hessian matrix. We successfully applied waveform inversion method from low to high frequency in two sets of Marmousi data. One is the data set generated by frequencydomain finite-difference modeling, and the other is the original Marmousi shots data set. The former result is very close to the true velocity model. In the original shots data set inversion, we replace the prior source with estimated source; the result is also acceptable, and consistent with the true model.
基金funded by the project‘Study on tomographic velocity inversion method based on P-wave and S-Wave constraints’(33550000-21-FW0399-0009)of State Key Laboratory of Shale Oil and Gas Enrichment Mechanisms and Effective Development and Sinopec Key Laboratory of Seismic Elastic Wave Technology.
文摘基于常规模拟退火算法的零偏VSP全波形反演面临着计算量大和耗时长的问题。为此提出了一种不同阶段对应不同扰动模型和退火方式的分段快速模拟退火(segmented fast simulated annealing,SFSA)反演策略,以提高零偏VSP资料全波形反演的效率。在反演前期采用大模型扰动空间和较慢温度衰减速度,充分发挥全局搜索能力,而在后期引入限制因子产生扰动模型,在迭代不断增加的时候逐渐减小模型的扰动空间,同时采用较快的温度衰减速度,有效提高反演的速度,使反演快速收敛到最优解。基于相同的初始温度和马尔可夫链长度,分别利用基于SFSA和非常快速模拟退火(very fast simulated annealing,VFSA)方法进行零偏VSP纵波速度全波形反演测试。结果表明,基于SFSA的反演方法的反演效率提高约50%,在迭代次数更少的条件下能获得更好的反演效果。基于SFSA的零偏VSP全波形反演具有高效和高精度的特点,其反演结果为地震地质层位标定、成果解释及油气预测奠定了基础。
