基于球面波近似算法的地震射线追踪改进

射线追踪是一种重要的地震数据处理工具,而地震波旅行时的计算是射线理论的核心问题。平面波近似旅行时计算方法精度低,且对于介质速度的适用性弱,因此本文在精度较高的球面波近似旅行时计算方法的基础上,对射线追踪算法提出了改进。数值模拟结果表明了该算法的正确性,能够有效地处理速度畸变界面,在不影响计算精度的情况下提高计算效率。     

基于改进的双线性旅行时插值的三维射线追踪

旅行时射线追踪的精度和效率一直是影响地震层析成像质量和效果的关键因素。这里在三维近地表速度结构层析成像应用中,结合旅行时插值算法与最短路径算法,提出了一种三维初至波射线追踪算法。通过对双线性旅行时插值算法进行改进,并运用判定条件与简化插值计算公式进行快速计算,减少了插值次数,降低了运算量,同时也保证了较高的射线精度,有效地解决了三维射线追踪算法的计算低效率,射线精度不高的问题。     

不同插值方法绘制断面图效果分析

介绍了绘制断面图中的各种插值方法,具体包括线性插值和分段非线性插值;分析了目前使用的线性插值方法绘制断面图的不足:精度一般、不光滑;提出了用分段非线性插值绘制断面图的方法;分段非线性插值精度高、光滑性好,所以分段非线性插值绘制断面图能很好地弥补线性插值的不足,而且它在计算机处理水文数据领域有不错的应用前景。     

旅行时线性插值射线追踪提高计算精度和效率的改进方法

旅行时线性插值（LTI）射线追踪算法是基于线性假设的,在向前处理过程中仅用按行或按列扫描的方法来计算节点旅行时没有考虑逆向传播射线,导致其计算精度与网格剖分大小有关,在处理复杂介质时会使得追踪出来的射线路径不一定满足最短旅行时。因此,笔者提出了两点改进措施：在向前处理时需采用全方位循环的方法来计算节点最小旅行时;在网格边界加入次生节点。模型试算结果表明：采用全方位循环的LTI法考虑了来自各个方向的射线,可提高其对复杂模型的适应能力;在节点间距相同的情况下,网格边界插入次生节点的LTI法较传统的LTI法计算精度至少可以提高一个数量级,同时,计算速度也更快;随着节点间距剖分的越精细,计算耗时下降也越明显,计算速度较传统的方法可提高n~10n倍。     

复杂海底各种地震波的射线追踪与运动学特征

为了实现适应崎岖海底、大陡坡、海底火山等复杂海底地质条件且灵活、稳定、精度高的射线追踪方法,并基于该方法详细分析复杂海底条件下各种地震波型的运动学特征,综合多种算法的优势,实现了一种快速推进迎风插值射线追踪方法。首先,采用混合网格法剖分复杂海底地质模型;其次,通过融入迎风差分思想的线性插值策略来构建精度高、无条件稳定且灵活的局部走时和射线路径计算公式;然后,综合应用这些公式和多级次快速推进法,灵活计算整个模型各种类型地震波的走时,并基于逆向追踪方法计算射线路径;最后,对该方法进行了精度分析,发现其能够获得相对高的走时和射线路径计算精度,且反射波的计算精度远高于入射波。此外,计算实例分析显示,初至波中富含折射波和陡倾构造的反射波在很大偏移距还能被接收,崎岖海底各种波型传播路径复杂;基于此提出加大采集排列长度和采用直达波走时可提高复杂构造成像质量等复杂海底地震数据采集与处理方面的思考与建议。     

桥梁检测中非结构化网格层析正演算法研究

在桥梁检测中弹性波层析正演是关键步骤。常规的层析正演中通常采用射线追踪算法,由于受规则矩形网格的影响,存在射线不能覆盖的异常区,射线旅行时误差曲线波动较大,进而导致累积误差较大,不利于进一步数值分析与校正。这里利用自适应三角网剖分技术,结合线性插值射线旅行时算法,分别对规则匀速模型、非规则匀速和变速模型进行试算,其射线路径分布结果和旅行时误差曲线结果,表明了该算法的正确性与可靠的计算精度,印证其算法对桥梁内部缺陷的检测是可靠的。     

利用线性旅行时插值射线追踪 计算近地表模型初至波走时

利用线性旅行时插值射线追踪对近地表模型进行正演计算,可以快捷、准确地获得初至波走时和射线路径。由于该算法计算的初至波不局限于折射波,因此很好地解决了浅层折射勘探中的低速"隐蔽层"问题;而且,由于该算法是基于网格划分和线性插值,因此它不仅可以追踪任意复杂介质的初至波,而且可以使得追踪的初至波射线路径逼于真实,避免了同类算法直接连接网络节点形成射线路径的缺陷(路径过于弯折,计算走时偏大)。将LTI算法同其他几种算法的追踪结果进行的对比和分析表明,LTI算法在计算初至波走时和射线路径方面较其他算法更为精准、稳定,是一种有效的射线追踪方法。     

一种变网格差分的快速行进法

地震波旅行时精度直接影响着地震反演、叠前偏移成像、层析成像等各领域研究成果的可靠性,因此,研究如何提高地震波旅行时精度是很有意义的。在双重网格技术的基础上,引入一种基于变网格差分格式的快速行进法(FMM)计算地震波旅行时,通过正演模拟均匀模型、存在高速异常体模型、Marmousi模型来分析变网格FMM的优势及适用性。研究结果表明:均匀模型背景下,变网格FMM与双重网格FMM优势相当,但是在存在高速异常体模型背景下,双重网格FMM可能违背波前扩展的规律,从而导致较大的误差,而变网格FMM则不存在这样的问题; Marmousi模型试算验证了变网格FMM能适应各种复杂模型。因此,该方法是一种有效提高走时计算精度和效率的方法,不仅增强了FMM法的适用性,而且扩展了变网格技术的应用范围。     

基于S变换和瞬时旅行时属性的自动初至拾取方法研究

地震波走时拾取是地震资料前期处理过程中的一个重要步骤,同时也是后期层析成像方法中计算速度的重要依据。本文提出了一种新的自动初至拾取的方法,利用地震信号的瞬时旅行时属性进行地震波自动初至拾取,瞬时旅行时属性的计算包括两个独立但相关联的部分,首先通过S变换将地震信号转换到时频域,然后利用映射算子将时频域表征映像回时间域得到瞬时旅行时,从而实现自动拾取初至。理论模型和实际地震资料处理的结果表明,该方法可以有效地拾取地震波初至时间。     

地震波走时计算的逆风差分算法

快速精确地旅行时计算在地震资料的叠前偏移与层析反演中起重要作用。利用逆足差分格式求解程函方程，首先在波阵面上寻找全局极小点，然后计算该点周围各点的地震波初至走时，同时，采用局部算法以确保计算的稳定性。该算法对于任意复杂地质模型都具有较高的精度。     

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

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

一种快速有效的地震道空间内插方法

地震道空间内插所面临的主要问题是空间假频。为此,Ｓｐｉｔｚ（１９９１年）提出了一种ｔ－ｘ域道内插方法可以很有效地克服空间假频,其数学基础是单位步长预测滤波技术。最近,在Ｓｐｉｔｚ工作的基础上,Ｐｏｒｓａｎｉ（１９９９年）提出了一种半步长预测滤波技术,它可以将原来用于求解预测滤波器系数的两组复杂的复线性系统方程简化为一组更为简单的复线性系统方法,从而大大降低了运算工作量。合成数据及实际地震记录的试处     

投影解析速度建模

当今地球动力学的研究已逐步引起重视 ,尤其是高新探测技术的发展 ,将会大大加快其研究步伐。地球物理勘探作为重要的探测技术之一 ,已在该领域做出了重要贡献 ,比如能源探测、环境保护、防灾减灾工作中都需要物探技术。文中就地震勘探领域中地震波成像前沿技术———深度偏移的核心问题———速度模型的建立展开深入的研究。在分析现有方法缺陷的基础上 ,提出投影解析速度建模方法。该方法建立在两个重要事实依据之上 :(1)地震波旅行时误差是由描述地下反射界面的深度误差和描述地层的速度误差共同引起的 ;(2 )在反演中要保证零炮检距旅行时不变。在局部平面波前假设和拉东投影原理下 ,导出由旅行时误差计算界面深度误差和速度误差的解析公式。该方法可保证界面深度和层速度同步快速收敛。与现有方法相比 ,该解析法可以明显提高计算效率 ,且应用效果相当明显。     

基于三次插值的大地电磁自适应有限元二维正演模拟

大地电磁(MT)数值模拟中通常使用有限单元法,通过伽辽金(Galerkin)法将微分方程转化为与其等价的泛函形式,对泛函求取极值并在单元上定义插值基函数,得到节点上电磁场值的线性方程组,最终形成大型复对称稀疏矩阵。要达到较高的有限元计算精度,一般采用密集的网格或高次插值的方法,这样做大大的减慢了正演的速度。结合两者的优点利用三次插值和h-型自适应相结合的有限元法来实现MT的正演算法。首先从一个粗网格出发并利用三次插值,通过后验误差估计方法局部加密网格,在计算量较小的情况获得较高的计算精度。这种方法可以针对目标区域和介质分界面发生突变处进行网格加密,不需要全局加密网格。最后通过对国际标准模型COMMEMI-2D1的模拟,分别比较二次插值与三次插值的自适应网格数量和数值模拟结果,证明了三次插值自适应有限元算法的可行性。     

一种加权抛物Radon变换地震波场重建方法

基于带限正反抛物Radon变换的基本原理和计算方法,笔者给出了一种加权抛物Radon变换（PRT）算法,该方法可应用于叠前地震波场恢复和重建,如缺失偏移距地震数据的插值外推和反假频重采样,其权系数完全由数据域来确定,也就是说在迭代计算过程中,该权系数也在变化。对于均匀采样的地震数据,笔者给出了近偏移距外推和波场的反假频重采样的计算结果,试算表明笔者给出的加权PRT方法可有效地应用于地震波场重建中,并具有算法稳健、计算精度高和效率高的特点。     

地震波混合阶褶积算法模拟

基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。