存在流体介质夹层时瑞利波频散曲线的完整求解

应用介质分界面上位移应力连续性的变化关系,采用δ矩阵法和固体层状介质中的瑞利波频散方程,推导了含流体介质夹层情况下瑞利波的频散方程。利用推导的频散方程对理论模型进行模拟计算,并与以往文献计算方法进行比较,模拟结果吻合较好,验证了公式的正确性;同时,模拟计算结果也说明了本文公式的普遍适用性和计算结果的完整性。根据文中对4个理论模型的模拟结果来看,与固体层状介质的情况相比,当存在流体夹层时,瑞利波频散曲线模式中多出了一个文中定义的0阶模式,在今后的多道瞬态瑞利波法应用中,应该可以将其看成是流体介质夹层存在的标志性曲线。     

基于F-K变换的井下多道瑞利波频散曲线提取

井下瑞利波勘探中频散曲线的提取多采用基于两道之间互相关的表面波谱分析法,此方法仅能提取基阶模式的瑞利波频散曲线.为了能够获得多阶模式频散曲线,提高煤矿井下瑞利波勘探的探测距离和精度,运用基于F-K变换的多道瑞利波频散曲线提取方法处理井下地震数据.针对井下地震数据的特点对方法进行了改进,在利用合成数据验证方法的有效性后,将其应用于工作面内部构造探测中,并从中获取了具有多阶模式的瑞利波频散曲线,从而证实了当煤层中存在断层时多阶模式瑞利波的存在.     

基于MPI的面波有限差分正演模拟

近年来,瑞利波波形反演技术因其避开了常规频散曲线计算,直接进行波场计算和反演不再受水平层状介质理论假设的限制,得到广大学者的高度重视。但瑞利波波形反演过程中需要不断进行波场正演和逆推计算。另外,由于浅地表速度较小,模拟计算时需要较小的网格间距才能避免数值频散,这无疑大大增加了正演模拟的计算量。对于这一问题,通常采用并行化设计来提高正演模拟的计算效率。本文基于消息传递接口（MPI）并行有限差分算法,以区域分解思路将模型区间分解成若干子区域,各区域互相通信,共同完成对模型的正演计算。并详细给出了区域分解、坐标转换、区域通信、波场合并等并行方案中的具体实现方法和实现步骤。通过对弹性模型、Kelvin黏弹性模型和标准线弹性固体（SLS）黏弹性模型不同并行方案的计算结果进行分析,验证了本文并行方案的可行性和有效性。并行计算结果表明,与单处理器计算时间相比,增加处理器数目可以明显减少计算时间,但随着处理器数目的增加,不同处理器之间的通信时间也增大;因此,并行时需要选择合适的处理器数目。对于黏弹性介质模型,SLS黏弹性模型的并行计算效率优于Kelvin黏弹性模型。     

Matlab环境下瑞利波有限差分正演与曲线绘制

以瑞利波频散方程为出发点的Abo-Zena传递矩阵等方法,只能研究层状介质中瑞利波的传播特性,对于非层状介质,传递矩阵无能为力。因此,不得不考虑有限元、有限差分等方法。根据弹性动力学方程,采用交错网格有限差分方法对均匀弹性半空间介质进行全波场模拟,并在Matlab环境下实现编程计算,再现了瑞利波在近地表的传播状态,实现了地震剖面曲线的正振幅充填问题,从而获得了类似野外地震勘探的剖面记录。模拟结果表明,采用一阶差分格式所得到的地震记录存在较严重的数值频散,在采用有限差分法研究瑞利波"之"字形成因时,须尽量压制数值频散现象,才能获得更接近实际情况的频散曲线。      

基于贝叶斯理论面波频散曲线随机反演

面波频散曲线反演是获得地下横波速度结构的重要地球物理方法。常规基于迭代最小二乘等线性反演方法依赖于初始模型,且存在多极值、容易陷入局部最小、反演精度低等问题。基于贝叶斯理论的随机反演方法是一种可以融合先验信息的非线性反演方法,该方法无需人为给定初始模型,仅利用先验信息对模型进行随机采样,根据概率分布筛选接受合适的后验概率密度估计结果,可达到对细节信息的准确估计。本文针对瑞利面波频散曲线,提出了基于GPR数据先验资料约束的贝叶斯马尔科夫蒙特卡洛(MCMC)随机反演方法,通过随机改变模型参数并计算其频散曲线与实际频散曲线的似然函数来选择是否接受新的模型参数,不断重复此过程,最终得到与实际频散曲线拟合效果最佳的最优解以及横波速度解的后验概率密度分布。通过理论模型以及实际数据反演测试,验证了该方法与常规无约束的随机反演相比,可以有效地提高反演速度和反演精度。     

基于多模态分离的面波谱分析方法

面波频谱分析法(SASW) 作为一种横向分辨率较高的瑞雷面波勘探方法, 由于计算得到的基阶面波频散曲线存在较大误差以及无法获得高阶模式频散曲线而在应用中受到限制.通过应用Fourier变换方法(FT法) 分离提取面波各模态数据, 进而对分离后的各模态数据利用SASW法计算频散曲线.通过模型实例分析得出: (1) 利用SASW计算基阶面波频散曲线时必须分离得到基阶面波数据, 然后计算才能得到正确的结果; (2) 基于多模态分离的SASW法可以计算得到面波高阶模式的频散曲线.      

基于Python语言的瑞利面波频散反演实行方案

本文基于Python语言构建了瑞利面波频散反演工作流程：(1)通过pysurf96软件包实现水平层状模型的频散曲线正演;(2)建立描述频散曲线拟合度的目标函数;(3)利用scikit-opt软件包中的启发式算法实现频散曲线反演。提出并解决了调用函数过程中所遇到的问题,实验测试表明反演结果可靠,同时具备一定的计算效率。由此实现了搭建基于Python语言的面波频散反演地下层状结构研究平台,为其他研究者依靠开源软件进行反演运算提供了方法支持。最后,通过利用彭一波研究海拉尔盆地背景噪声时提取的频散曲线对地壳以及上地幔结构进行了反演,得到不错的效果。     

多道瞬态面波在复杂地形条件下岩层划分中的应用研究

在复杂地表条件下,利用瑞雷面波进行勘探时,需要考虑地形对瑞雷波场传播地影响。针对倾斜地表两层地质模型,通过设计不同展布方向的检波器排列进行瑞雷波信号采集,利用高精度交错网格有限差分法,讨论经过快速傅里叶变换后的地震记录F-V频率谱,通过对比信号频率谱与理论多阶频散曲线,分析在不同地形条件下瑞雷波传播能量的变化特征,并提取叠加频散曲线;运用遗传算法对采集到的地震记录进行基阶与高阶瑞雷波联合反演得到目的层厚与初始模型层厚进行误差对比;并对理论正演初始模型的频散曲线与各地质模型利用F-K法提取的频散曲线进行全局拟合误差计算。分析表明:倾斜地形进行多道瞬态瑞雷面波勘探,相对于下部激发、上部激发会产生高阶频散趋势,因此对下部激发的地震记录进行反演产生的厚度误差相对较小。     

基于小波域最小平方滤波的多尺度自适应全波形反演

常规的全波形反演对初始模型和低频数据的依赖性较强,反演精确度经常会受到"跳周"现象的严重影响。将小波域最小平方滤波引入到全波形反演中,并利用小波变换的多尺度特性,有效地提高了反演过程的稳定性,减小反演受到"跳周"现象影响的可能性。小波域最小平方滤波具有更高的精度和更好的性能,利用该算法调整模拟波场的相位,从而改善模拟波场和观测波场之间的相位差异,并由此构造一个新的目标函数。同时利用小波变换的多尺度特性将地震数据分解为不同频带数据,实现多尺度反演。数值模拟实验结果表明,基于小波域最小平方滤波的多尺度自适应全波形反演,对初始模型和低频数据的依赖程度降低,较常规全波形反演具有更好的稳定性,受到"跳周"现象影响的可能性大大降低。     

基于遗传算法的地震面波和大地电磁数据一维联合反演算法研究

利用快速标量传递算法进行了一维层状模型的地震面波频散曲线正演计算,得到了基阶频散曲线;利用传统方法对一维层状模型进行了大地电磁正演计算,得到大地电磁的视电阻率曲线数据和相位数据。利用遗传算法对合成数据分别进行了单独反演和联合反演,得到拟合差分布图来分析两种数据的兼容性并研究联合反演的可行性。通过对比单独反演和联合反演的结果,可以看出地震面波频散曲线和大地电磁数据的联合反演结果较之单方法的反演结果要更加精确,收敛速度更快,两种数据的耦合情况较好,同时联合反演能够有效地限制多解性。     

从地震层析速度确定原位岩土力学参数

地震纵横波层析成像速度与岩土弹性力学参数有关,多道面波分析（MASW）采用类似二维反射地震处理,得到共中点互相关道集（CMPCC）提取面波频散曲线,由面波频散曲线反演二维速度结构,在岩土结构分层等方面有良好的效果。通过同一剖面的折射地震波层析成像得到纵波速度结构,对岩土层的动弹模量和泊松比等力学参数进行分析,为下一步地质工程设计及岩土力学试验提供参考。     

地震资料面波的一种特殊处理技术

作者在本文中对地震资料的面波作了如下研究：（1）引用地震波复速度及地震波能量吸收的形式来描述粘弹性介质；（2）应用目前国外粘弹性介质面波频散参数计算的快速算法，计算出了多阶振型地震面波的各阶频散参数；（3）运用哈特莱变换计算出地震资料中的面波波形。     

瑞利波频散成像方法的实现及成像效果对比研究

瑞利波勘探方法广泛应用在各种工程勘察领域中,其中快速准确的瑞利波频散成像方法是进行瑞利波勘探方法的关键因素之一。这里通过理论数据详细介绍了三种频散成像方法的实现过程,并深入分析了不同方法的频散成像效果,通过对比发现,F-K变换法和τ-p变换法虽可以清晰地看出瑞利波能量的分布特点,但并不利于频散曲线的提取,而相移法则具有更好的频散成像效果,可以得到更高频率和更高阶次的瑞利波频散曲线。     

储层中的导波

储层常常表现为低速,在井间地震中往往可以观察到导波的存在,其主要特点是在储层附近能量集中,而在围岩中则能量迅速衰减,同时波形具有频散特征。导波特性对储层内部结构与岩性十分敏感,它有可能作为储层描述和采油过程监测的有效手段,本文针对三层简单波导模型进行了频散方程及曲线的计算,并进行了人工合成槽波记录的研究。     

煤层中瑞利型槽波的频散特性

频散是槽波的最大的特征,基于此推导三层对称模型的瑞利型槽波的频散特征方程并绘制出频散曲线,对频散曲线进行研究,讨论地层参数与埃里相频率的关系(影响频散的因素)。通过研究表明,利用已知频散曲线可以进行地层厚度、速度等参数的反演和估算。     

利用正演模拟实现面波衰减

瑞雷波是分布在自由表面附近的P-SV面波,具有强振幅、低频、频散的特征,常用于工程物探,而在油气勘探及深部探测中常作为一种"干扰"波,会影响地震资料的信噪比,需要将面波净化处理。为实现面波衰减,本文将工程物探方法和石油物探技术相结合,通过面波频散曲线反演获得高精度的近地表模型,在三维近地表横波速度模型基础上进行面波正演,利用面波模型自适应匹配相减的方法将其减去,在不伤及有效波的情况下将面波的影响降至最低,一定程度上避免了常规面波衰减方法伤及有效波、损失低频弱信息、导致资料有效频带变窄、对地震构造解释及反演造成负面影响等缺点。最后应用实际数据验证了本研究方法的有效性。     

基于最佳小波基的地震面波插值方法

在利用实际地震数据中的面波反演近地表横波速度的过程中,若道间距较大、空间采样率不足,则会产生空间假频现象,从而降低频率速度谱的信噪比,影响频散曲线提取的精度以及反演效果,因此需要针对面波进行插值处理。文中提出了一种基于最佳小波基的地震面波插值方法,通过理论分析和实验误差对比在地震数据处理常用的众多小波基中选出适用于插值处理的最佳小波基bior6.8,提高了插值精度。针对面波同向轴为线性且斜率较大的特点,文中首先采用线性动校正的方法对面波进行拉平处理,再进行小波变换插值,最后进行反线性动校正恢复面波。通过对理论模型与实际资料进行插值处理验证了本文方法的有效性,插值后的面波记录波形恢复较好,显著提高了频率速度谱的信噪比,有效解决了面波数据空间采样率不足引起的假频问题。     

地震层析成像方法综述

地震层析成像作为一种有效还原地下介质速度模型的方法,为全波形反演提供了可靠的初始速度模型,从原始的射线层析到相移旅行时层析和瞬时旅行时层析,实现了地震波传播的有限频特性;从声波方程到弹性波方程,从各向同性介质到VTI,TTI介质,实现了对真实地下介质情况的模拟。减缓层析反演的病态性也一直是研究热点,常用的方法有正则化,用高斯束层析的敏感核代替传统的射线层析敏感核等。此外,为了避免使成像结果的精度依赖于共成像道集上反射位的真实深度,角度域双差分反射层析可以稳定有效地收敛到精确的偏移速度模型。如今,层析成像逐步向各向异性介质过渡,使用的数据从VSP到WVSP,从单一波形到多种波形联合反演发展,然而,分辨率和计算效率的相关问题仍然需要得到关注。     

工程地质中瑞利波法勘探的理论初探

本文以水平层状液体为模型,推导出瑞利面波的频散方程,进而作了多层介质的模型试算并绘出了瑞利波波速与波长的关系曲线.该曲线对应于弹性分界面上的异常特征与实测曲线的异常特征吻合,可用于说明层状介质中瑞利面波勘探的机制.文中还以理论模型试算结果为依据,对当前面波勘探中的半波长解释方法作了评述,并解释了固体层状介质中面波曲线的异常特征.     

WD智能天然源面波数据采集处理系统及其应用试验

我国自主研发的WD智能天然源面波数据采集处理系统,改变了迄今为止现场仪器只采集天然微振动波形记录的现状。该系统利用大数据计算、筛选、叠加、归一化等技术,创先性研制出现场仪器在采集微震信号的同时,可直接显示勘探成果的面波频散曲线;信号采集过程中,无需震源,无需人工控制;装置系统安置在地面,10~30min可获得上百米深度地层的面波频散曲线,得知土层软硬。在不同岩性地区的试验结果证实,该系统能快速、准确地进行地层划分,操作便捷,适应性强,为天然源面波勘探广泛应用于工程勘查提供了技术保障。