基于局部优化粒子群算法的背景噪声反演

从背景噪声中提取瑞雷波频散曲线并通过反演获得地下横波速度结构已被广泛应用于大尺度的地下结构探测和小尺度的工程勘探中.基于频散函数的反演目标函数可以有效解决多阶模频散曲线联合反演的模式误判问题,然而其广泛分布的局部极值导致更为严重的多解性,在大范围的参数搜索空间下很难获得最优解,需要搭配全局搜索性能强的优化算法.本文提出局部优化粒子群算法(PSOG),通过粒子迭代过程中引入局部优化方法提高种群多样性,避免陷入局部极值并加快收敛速度.为验证新算法的有效性,结合基于久期函数的目标函数对理论合成数据进行反演,结果表明,局部优化粒子群算法比传统算法的稳定性与准确性都有显著提高.处理了上海苏州河地区的背景噪声数据,成功地对古河道切割造成的软弱层进行成像.PSOG算法与新型反演目标函数的结合在背景噪声勘探的工程应用上具有巨大潜力.

Optimization method of fi rst-arrival waveform inversion based on the L-BFGS algorithm

The fi rst arrival waveform inversion(FAWI)has a strong nonlinearity due to the objective function using L2 parametrization.When the initial velocity is not accurate,the inversion can easily fall into local minima.In the full waveform inversion method,adding a cross-correlation function to the objective function can eff ectively reduce the nonlinearity of the inversion process.In this paper,the nonlinearity of this process is reduced by introducing the correlation objective function into the FAWI and by deriving the corresponding gradient formula.We then combine the first-arrival wave travel-time tomography with the FAWI to form a set of inversion processes.This paper uses the limited memory Broyden-Fletcher-Goldfarb-Shanno(L-BFGS)algorithm to improve the computational effi ciency of inversion and solve the problem of the low effi ciency of the FAWI method.The overthrust model and fi eld data test show that the method used in this paper can eff ectively reduce the nonlinearity of inversion and improve the inversion calculation effi ciency at the same time.

Reflection full-waveform inversion using a modified phase misfit function

Reflection full-waveform inversion （RFWI） updates the low- and high- wavenumber components, and yields more accurate initial models compared with conventional full-waveform inversion （FWI）. However, there is strong nonlinearity in conventional RFWI because of the lack of low-frequency data and the complexity of the amplitude. The separation of phase and amplitude information makes RFWI more linear. Traditional phase-calculation methods face severe phase wrapping. To solve this problem, we propose a modified phase-calculation method that uses the phase-envelope data to obtain the pseudo phase information. Then, we establish a pseudophase-information-based objective function for RFWI, with the corresponding source and gradient terms. Numerical tests verify that the proposed calculation method using the phase-envelope data guarantees the stability and accuracy of the phase information and the convergence of the objective function. The application on a portion of the Sigsbee2A model and comparison with inversion results of the improved RFWI and conventional FWI methods verify that the pseudophase-based RFWI produces a highly accurate and efficient velocity model. Moreover, the proposed method is robust to noise and high frequency.

Back analysis for soil slope based on measuring inclination data

Based on the analysis of several objective functions,a new method was proposed.Firstly,the feature of the inclination curve was analyzed.On this basis,the soil could be divided into several blocks with different displacements and deformations.Then,the method of the soil division was presented,and the characteristic of single soil block was studied.The displacement of the block had two components:sliding and deformation.Moreover,a new objective function was constructed according to the deformation of the soil block.Finally,the sensitivities of the objective functions by traditional method and the new method were calculated,respectively.The result shows that the new objective function is more sensitive to mechanical parameters and the inversion result is close to that obtained by the large direct shear apparatus.So,this method can be used in slope back analysis and its effectiveness is proved.

Study of frequency domain full waveform inversion based on Huber norm and L-BFGS algorithm

Full waveform inversion( FWI) is a high resolution inversion method,which can reveal detailed information of the structure and lithology under complex geological background. It is limited by many kinds of noises when the method applied to the real seismic data. Based on Huber function criterion,the objective function combinates the anti-noise of L1 norm and the stability of L2 norm in theory,the authors derive the gradient formula of the Huber function by using L-BFGS algorithm for FWI. The new method is proved by synthetic seismic data with the Gaussian noise and the impulse noise. Numerical test results show that L-BFGS algorithm is applied to the frequency domain FWI with the convergence speed and high calculation accuracy,and can effectively reduce computer memory usage; and the Huber function is more robust and stable than L2 norm even with the noises.

Robust elastic impedance inversion using L1-norm misfit function and constraint regularization

The classical elastic impedance （EI） inversion method, however, is based on the L2-norm misfit function and considerably sensitive to outliers, assuming the noise of the seismic data to be the Guassian-distribution. So we have developed a more robust elastic impedance inversion based on the Ll-norm misfit function, and the noise is assumed to be non-Gaussian. Meanwhile, some regularization methods including the sparse constraint regularization and elastic impedance point constraint regularization are incorporated to improve the ill-posed characteristics of the seismic inversion problem. Firstly, we create the Ll-norm misfit objective function of pre-stack inversion problem based on the Bayesian scheme within the sparse constraint regularization and elastic impedance point constraint regularization. And then, we obtain more robust elastic impedances of different angles which are less sensitive to outliers in seismic data by using the IRLS strategy. Finally, we extract the P-wave and S-wave velocity and density by using the more stable parameter extraction method. Tests on synthetic data show that the P-wave and S-wave velocity and density parameters are still estimated reasonable with moderate noise. A test on the real data set shows that compared to the results of the classical elastic impedance inversion method, the estimated results using the proposed method can get better lateral continuity and more distinct show of the gas, verifying the feasibility and stability of the method.

地震勘探全波形反演的应用与发展分析

本文首先对20世纪80年代发展起来的全波形反演应用及其在勘探地球物理领域的发展进行了分析;其次,面对定量化、精细化的地震勘探要求,提出了将地震勘探全波形反演与其他数据处理环节或处理技术相结合的研究设想,并展望了全波形反演的发展趋势;最后,论述了全波形反演研究中地震波场数值模拟、反演初始速度模型获取、目标函数形式选择、寻优算法启用及各向异性介质中的应用等关键问题,并总结了通过Laplace域的全波形反演获取反演初始速度模型、结合射线追踪并充分发挥并行计算之于波动方程方法来模拟地震波场的巨大优势,及灵活选用反演目标函数形式和寻优算法更新速度模型参数来加快全波形反演方法的实用化进程。