首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 171 毫秒
1.
包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.  相似文献   

2.
包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.  相似文献   

3.
瑞雷波频散曲线反演的本质是对目标函数求极值的过程.传统的线性局部反演算法容易陷入局部极小值,增加了反演结果的不确定性.粒子群算法作为一种全局非线性优化手段,能够保证各搜索空间内局部寻优的同时,逼近全局最优,保证迭代反演的收敛性.建立含软夹层型地层模型,正演计算得到理论频散曲线,采用粒子群算法反演得到的横波速度与模型高度吻合.同时,在理论频散曲线中加入10%的高斯白噪声,粒子群算法的反演结果仍然具有很好的可靠程度.通过对实测频散曲线的反演表明:粒子群算法在含软夹层型地层的横波速度探测中相对于最小二乘法优势明显,具有很高的应用价值.  相似文献   

4.
田宵  汪明军  张雄  张伟  周立 《中国地震》2021,37(2):452-462
微地震事件的空间分布可以用来监测水力压裂过程中裂缝的发育情况.因此,震源定位是微震监测中重要的环节.震源定位依赖准确的速度模型,而震源位置和速度模型的耦合易导致线性迭代的同时反演方法陷入局部极小值.邻近算法作为一种非线性全局优化算法,能够最大程度地避免陷入局部最优解.本文将邻近算法应用于单井监测的微震定位和一维速度模型...  相似文献   

5.
基于L2范数的常规全波形反演目标函数是一个强非线性泛函,在反演过程中容易陷入局部极小值.本文提出归一化能量谱目标函数来缓解全波形反演过程中的强非线性问题,同时能够有效地缓解噪声和震源子波不准等因素的影响.能量谱目标函数是通过匹配观测数据与模拟数据随频率分布的能量信息来实现最小二乘反演的,其忽略了地震数据波形与相位变化的细节特征,这在反演的过程中能够有效缓解波形匹配错位等问题.数值测试结果表明,基于归一化能量谱目标函数在构建初始速度模型、抗噪性和缓解震源子波依赖等方面都优于归一化全波形反演目标函数.金属矿模型测试结果表明,即使地震数据缺失低频分量,基于归一化能量谱目标函数的全波形反演方法在像金属矿这样的强散射介质反演问题上同样具有一定的优势.  相似文献   

6.
从MT信号中提取激电信息的研究对深部矿产资源和油气资源的开发具有重要的现实意义.在激电信息提取方法的研究中,过去的理论多以线性反演方法为主,主要存在三个问题:(1)依赖初始模型,(2)容易陷入局部极值,(3)多解性严重.本文通过以下方法解决上述问题:(1)采用差分进化算法对MT信号中的激电信息以反演的方式进行提取,该方法受初始模型影响小,不易陷入局部极值.(2)在差分进化算法的适应度目标函数中引入最小构造约束,并提出对极化率和电阻率分别约束来提高反演结果的稳定性.对极化层处于不同位置的不同地电类型的反演仿真结果表明,本文算法能有效获得电阻率和激电参数的光滑模型,反演结果具有较高的稳定一致性和准确性.加入噪声测试后的实验结果表明,本文算法对高斯白噪声有较强的鲁棒性.与微粒群优化算法(PSO)比较结果表明,本文算法具有更为优越的收敛速度,能够获得更好的反演效果.上述结果表明,用非线性方法提取MT信号激电信息具有积极的意义.  相似文献   

7.
全波形反演可以为叠前深度偏移成像提供更高精度的速度模型,但该方法具有较强的非线性,对初始速度模型的依赖性较强,尤其是在实际应用中,地质条件复杂多变,速度变化不连续,增加了反演非线性程度,常常使反演陷入局部极小值,影响反演的精度.全变差约束在图像去噪领域应用广泛,属于非光滑约束,在去噪过程中能有效的保留图像的不连续界面和边缘信息.本文提出基于Hinge损失函数的垂向全变差约束全波形反演方法,在全变差约束的基础上,利用Hinge损失函数控制模型的更新方向,并使用原-对偶混合梯度算法进行求解,给出这一优化问题的迭代格式,有效提高了对地下不连续界面的重构精度,同时也降低反演对初始速度模型的依赖程度.数值算例证明:与常规全波形反演方法相比,基于全变差约束的全波形反演方法可以有效的重构速度模型中的不连续界面,尤其对高速体边缘的重构效果更明显,但该方法对初始速度模型的依赖性仍然较强;基于Hinge损失函数的垂向全变差约束全波形反演方法降低了对初始速度模型的依赖程度,可以从一个较差的初始模型通过循环迭代的方式最终得到同样精确的速度模型,较好的重构了高速体边缘和不连续界面.  相似文献   

8.
基于多网格的频率域全波形反演(英文)   总被引:2,自引:1,他引:1  
频率域全波形反演虽然克服了时间方向上的局部极小值问题,但是地下介质的复杂性使其在空间域仍然存在局部极小值缺陷。在优化梯度法基础上,本文采用预条件双共轭梯度稳定算法和多重网格方法计算反演中的波场传播和目标函数的梯度,在保证计算速度的同时,减小计算机内存的消耗。频率域波形反演和多重网格的多尺度性质有效改善问题极小值缺陷,加快反演的收敛速度。以局部非均匀的三孔模型和Marmousi模型的数值模拟结果验证了该算法的有效性。  相似文献   

9.
多分量地震数据中低频缺失是弹性波全波形反演中的一大难题,低频的缺失导致全波形反演无法有效恢复介质的长波长成分进而使反演陷入局部极值。为此,本文提出了反演介质纵横波速度长波长分量的弹性波包络反演方法。该方法利用包络算子具有的解调多分量数据中隐含的低频信息的能力,构造多分量地震数据的包络目标函数进行反演,用以恢复地下介质纵横波速度的长波长成分。一系列数值试验表明,即使在多分量地震数据中缺失低频信息、并且初始模型缺少先验信息的情况下,这种弹性波包络反演方法能够有效降低波形反演的非线性,可以为后续的常规弹性波全波形反演或者深度偏移提供足够精确的初始模型,且该方法对横波速度长波长分量的重建尤为有效。Mamousi-2模型的高精度纵横波速度的反演结果表明,利用该方法反演的纵横波速度作为常规弹性波全波形反演的初始模型,可以显著提高反演结果的精度。此外,本文对弹性波包络反演方法的适用性也进行了初步的研究与讨论。  相似文献   

10.
从CSAMT信号中提取激电信息有利于提高频率域电磁法反演与解释的精度.目前的研究多以线性反演方法为主,存在依赖初始模型、易陷入局部极值的问题.针对CSAMT信号IP提取问题的非线性和非凸特征,本文提出了一种基于柯西分布和惯性权重的二阶段最小构造混合蛙跳反演方法来提取IP信息.该方法首先利用柯西算子取代随机算子来提高算法的全局搜索能力,并通过引入混沌震荡惯性权重来均衡进化过程中的个体经验和群体经验,保证算法后期的稳定收敛;然后通过引入第二阶段反演过程来强化极化率对观测数据的影响,同时将正则化参数引入混合蛙跳算法的适应度函数来改善反演的多解性问题;最后利用CPU并行计算加速了算法的模因组搜索过程.反演结果表明,上述方法能够较好地重构地电结构和提取激电信息,在加噪环境下具有较强的鲁棒性.相比其他非线性算法(标准混合蛙跳算法SFLA,差分进化算法DE和粒子群优化算法PSO)的反演结果,本文算法具有更强的全局搜索能力和更高的计算效率,适合对微弱的激电信息进行提取.  相似文献   

11.
反射地震数据的逐层波形反演   总被引:10,自引:2,他引:8       下载免费PDF全文
本文针对层状介质并结合梯度法波形反演,提出逐层波形反演的方法. 首先给出介质扰动响应的概念,并在此基础上分析了梯度法波形反演方法. 波形反演实质上是将实测地震记录和预测地震记录的波形残差信息转化为实际地质模型和预测地质模型的模型残差信息. 波形反演的优点是利用大量振幅相位信息得到高分辨率的反演结果, 其缺点是运行耗时大;当初始模型和实际模型相差较大时,迭代算法容易陷入局部极小点,这是因为目标函数和初始模型同实际模型间的差异是非线性的关系. 逐层波形反演方法是使自上而下每一层的目标函数最小,这样总的目标函数也是最小的. 利用二分法速度扫描确定每一层速度不仅提高了运算速度也避免了迭代算法陷入局部极小点的问题. 结合介质扰动响应和目标函数值变化可以更为准确迅速地确定每一层速度和该层界面位置.  相似文献   

12.
Waveform inversion met severe challenge in retrieving long‐wavelength background structure. We have proposed to use envelope inversion to recover the large‐scale component of the model. Using the large‐scale background recovered by envelope inversion as new starting model, we can get much better result than the conventional full waveform inversion. By comparing the configurations of the misfit functional between the envelope inversion and the conventional waveform inversion, we show that envelope inversion can greatly reduce the local minimum problem. The combination of envelope inversion and waveform inversion can deliver more faithful and accurate final result with almost no extra computation cost compared to the conventional full waveform inversion. We also tested the noise resistance ability of envelope inversion to Gaussian noise and seismic interference noise. The results showed that envelope inversion is insensitive to Gaussian noise and, to a certain extent, insensitive to seismic interference noise. This indicates the robustness of this method and its potential use for noisy data.  相似文献   

13.
Reflection full waveform inversion can update subsurface velocity structure of the deeper part, but tends to get stuck in the local minima associated with the waveform misfit function. These local minima cause cycle skipping if the initial background velocity model is far from the true model. Since conventional reflection full waveform inversion using two‐way wave equation in time domain is computationally expensive and consumes a large amount of memory, we implement a correlation‐based reflection waveform inversion using one‐way wave equations to retrieve the background velocity. In this method, one‐way wave equations are used for the seismic wave forward modelling, migration/de‐migration and the gradient computation of objective function in frequency domain. Compared with the method using two‐way wave equation, the proposed method benefits from the lower computational cost of one‐way wave equations without significant accuracy reduction in the cases without steep dips. It also largely reduces the memory requirement by an order of magnitude than implementation using two‐way wave equation both for two‐ and three‐dimensional situations. Through numerical analysis, we also find that one‐way wave equations can better construct the low wavenumber reflection wavepath without producing high‐amplitude short‐wavelength components near the image points in the reflection full waveform inversion gradient. Synthetic test and real data application show that the proposed method efficiently updates the background velocity model.  相似文献   

14.
改进的模拟退火-单纯形综合反演方法   总被引:15,自引:1,他引:14       下载免费PDF全文
实际中的大量地球物理反演是一个多参数、非线性优化问题,所采用的目标函数,即度量由参数化的理论模型得出的预测值与观测值的吻合程度,往往具有多个局部极值.针对这类问题,本文综合全局反演方法具有的全域搜索能力强、局部方法收敛速度快和“均匀设计”布点效率高的特点,提出了模拟退火-单纯性综合反演方法,并通过一维声波非线性反演验证了这种综合方法的搜索能力和效率.  相似文献   

15.
改进的模拟退火-单纯形综合反演方法   总被引:19,自引:6,他引:13       下载免费PDF全文
实际中的大量地球物理反演是一个多参数、非线性优化问题,所采用的目标函数,即度量由参数化的理论模型得出的预测值与观测值的吻合程度,往往具有多个局部极值.针对这类问题,本文综合全局反演方法具有的全域搜索能力强、局部方法收敛速度快和“均匀设计”布点效率高的特点,提出了模拟退火-单纯性综合反演方法,并通过一维声波非线性反演验证了这种综合方法的搜索能力和效率.  相似文献   

16.
不依赖子波、基于包络的FWI初始模型建立方法研究   总被引:3,自引:3,他引:0       下载免费PDF全文
地震全波形反演(FWI)从理论走向实际面临着诸多难题,其中之一就是需要一个较高精度的初始模型,另一个难题就是需要一个较为精确的震源子波,初始模型和震源子波的准确程度严重影响着全波形反演的最终结果.为此,本文提出了不依赖子波、基于包络的FWI初始模型建立的方法,建立了相应的目标函数,推导出了反演的梯度,给出了伴随震源的表达式,理论上分析了不依赖子波FWI的可行性.在数值试验中,讨论了参考道的选取方式,通过分析归一化目标函数收敛速率,认为近偏移距参考道优于远偏移距参考道,在地震数据含干扰噪音时,平均道作为参考道要优于最小偏移距参考道.通过包络、包络对数、包络平方三种目标函数反演结果的比较,发现包络对数目标函数对深层的反演效果最好.通过不同子波的试验进一步验证了本方法的正确性.  相似文献   

17.
Highly non-linear seismic trace inversion problems can be solved efficiently by an implementation of Tabu Search, a meta-heuristic method related to artificial intelligence. The implementation under consideration is a deterministic, global search that combines the advantages ofa local search, giving a quick descent to local misfit minima, with an ability to cross misfit barriers in the model space. Once Tabu Search has found an area of low misfit, it performs an extensive exploration of its deepest points. This property makes it possible to use Tabu Search for a semiquantitative resolution and uncertainty analysis of the inverse problem.  相似文献   

18.
基于声波方程的井间地震数据快速WTW反演方法   总被引:4,自引:2,他引:4       下载免费PDF全文
WTW(Wave equation traveltime+Waveform inversion)反演是基于波动方程的走时反演(WT反演)和波形反演的联合反演方法.WT反演利用波动方程计算走时和走时关于速度的导数,和传统以射线为基础的走时反演相比,具有不必射线追踪、不必拾取初至、不必高频假设以及初始模型和实际模型差别较大时也能较好收敛等优点,但WT反演与波形反演相比其结果分辨率低.与之互补的是,波形反演的反演结果分辨率高,但是当所给初始模型和实际模型相差太大时,波形反演迭代算法容易陷入局部极小点.可见结合两种方法的WTW反演是一种比较好的联合反演方法.常规WTW迭代算法是首先以WT反演为主反演得到地质模型的整体特征,然后再以波形反演为主反演模型细节,该算法耗时和占用计算机存储空间接近WT反演或波形反演的两倍.为了节省运算耗时和计算机存储空间,往往采取首先单独利用WT反演然后再单独利用波形反演的算法.这样做的缺点是不能紧密结合两种反演方法,使得它们的优缺点在每一次迭代中无法得到互补,从而影响了最终的反演结果.针对以上事实,本文提出一种新的方法实现WTW,使得WTW运算速度和存储空间在任何情况下等同于WT反演或波形反演.模型计算表明新的算法具有更好的收敛性.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号