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

直达波旅行时线性插值算法黄翼坚;朱光明;敦敏【摘要】文中首先简要介绍了初至波旅行时LTI算法,并分析了该法的计算结果包含非直达波旅行时的原因,由此提出通过对射线方向进行控制实现直达波旅行时计算的方法.直达波旅行时计算步骤与初至波旅行时LTI算法的计算步骤基本一样,只有两个不同之处:其一,旅行时赋初值不同,对震源点以外的节点,初至波旅行时LTI算法赋予大干该节点初至波旅行时的值,而直达波旅行时LTI算法则赋予大干该节点直达波旅行时的值;其二,在每次更新节点旅行时的过程中,一旦发现有波从某个低速地层传向某个高速地层,则在以后的计算中,判断射线是否是从该高速层指向该低速层,如果是则舍去该射线的旅行时.平层速度模型和具有起伏速度界面的多层模型算例说明,本文提出的直达波旅行时LTI迭代算法简单而有效,反射波射线追踪算例说明基于该算法可以进行大角度反射波射线追踪.【期刊名称】《石油地球物理勘探》【年(卷),期】2010(045)002【总页数】5页(P225-229)【关键词】初至波旅行时;直达波旅行时;线性插值;算法;费马原理;数值模拟【作者】黄翼坚;朱光明;敦敏【作者单位】长安大学地质工程与测绘学院,陕西西安,710054;长安大学地质工程与测绘学院,陕西西安,710054;东方地球物理公司物探技术研究中心,河北涿州,072751【正文语种】中文复杂速度模型中地震波旅行时的计算是基尔霍夫积分法叠前深度偏移及层析成像的一个关键步骤。

自20世纪80年代末以来,国内外许多学者研究了初至波旅行时的计算方法。

例如,Vidale[1]提出了初至波旅行时程函方程有限差分解法;Moser[2]提出了最短路径射线追踪及初至波旅行时计算方法;A sawaka和 Kawanaka[3]提出了初至波旅行时线性插值(L TI)算法;刘洪等 [4]提出了界面网全局旅行时计算技术;王华忠等[5]提出了动态规划法旅行时计算方法;赵连锋等[6]提出了有序波前重建法射线追踪技术;秦义龙等[7]利用单频双程波动方程计算初至波旅行时。

旅行时线性插值射线追踪提高计算精度和效率的改进方法的报
告，600字
旅行时线性插值射线追踪技术是一种用于优化三维图形处理的技术，它能够在计算精度和效率方面取得令人满意的结果。

本文将介绍如何改进这一技术，来提高计算精度和效率。

首先，我们将介绍如何改进射线追踪技术，以增强其效率和精度。

一种方法是利用隐式函数曲面的求解，即通过划分曲面计算其法向量，从而替换曲面离散点之间的线性插值，减少不必要的计算量。

此外，我们还可以使用并行算法来加速射线追踪技术。

例如，可以使用多线程并行技术来加速射线追踪算法以及更复杂的加权函数曲面，从而提高算法的性能。

其次，为了提高射线追踪技术的精度，可以使用深度学习技术。

深度学习技术可以以许多方式改进射线追踪技术，例如使用深度神经网络来快速替换步骤中的点，从而大大减少误差。

此外，深度学习技术还可以为射线追踪技术带来附加的视觉信息，从而提高模型的精度和可读性。

最后，为了提高射线追踪技术的效率，可以使用图形处理器（GPU）。

GPU可以有效利用其并行处理能力来提高射线追
踪技术的效率，降低维数并引入缓存能够加快运算速度，从而提高渲染效率。

总而言之，旅行时线性插值射线追踪技术的改进方法包括利用隐式函数曲面的求解，使用并行算法进行加速，使用深度学习技术提高精度，以及使用GPU加速渲染的方法。

通过这些改
进方法，可以显著提高射线追踪技术的精度和效率，为三维图形处理注入新的活力。

线性走时插值射线追踪算法的改进卢江波;方志【摘要】在LTI(Linear Travel-time Interpolation)射线追踪算法基础上提出的扩张收缩扫描算法能正确追踪直达波、绕射波和回波的射线路径，但其存在计算效率低、收敛速度慢的问题。

采用交叉扫描方式对扩张收缩扫描算法进行改进，并由此提出了基于交叉扫描方式的扩张收缩扫描改进算法。

理论分析及数值模拟结果表明：改进算法在保留了原扩张收缩扫描算法所有优点的同时，具有更高的计算效率；当模型网格尺寸划分较细时，改进算法在计算效率上的优势更为显著。

%The extension & compaction scanning algorithm based on LTI (Linear Travel-time Interpo-lation)ray tracing algorithm can properly trace the ray path of the direct wave,diffraction wave and echo 's,but it has the following problems:low computational efficiency and slow convergencespeed.Aiming at these problems and using the cross-scan mode to improve the algorithm,an improved extension &compac-tion scanning algorithm based on the cross-scan mode was presented.Theoretical analysis and the results of numerical calculation show that the improved algorithm not only maintains all the advantages of the ex-tension & compaction scanning algorithm,but also has a higher computational efficiency.When the model is mesh more refined,the advantage of the improved algorithm is more significant in computational efficien-cy.【期刊名称】《湖南大学学报（自然科学版）》【年(卷),期】2014(000)001【总页数】6页(P39-44)【关键词】射线追踪;LTI算法;交叉扫描;扩张;收缩扫描法;改进算法【作者】卢江波;方志【作者单位】湖南大学土木工程学院，湖南长沙 410082;湖南大学土木工程学院，湖南长沙 410082【正文语种】中文【中图分类】TU317射线追踪技术在地震层析成像以及混凝土超声波射线层析成像等领域具有重要作用.目前常用射线追踪方法主要有有限差分解程函方程法［1－2］、最短路径法［3－4］以及LTI（Linear Travel－time Inter－polation）射线追踪算法［5－6］等.实验表明［5］，LTI算法在走时计算以及射线路径追踪上比其它方法（如有限差分解程函方程法）更为快速、精确.LTI射线追踪算法由Asakawa等人提出［5］，该算法以走时线性变化为前提，分两步进行.第一步，向前计算最小走时：先将模型划分为若干规则的单元，并将各单元边界划分为若干节段，然后，根据走时线性变化的假定以及走时最小原理（Fermat原理），求出从发射点到接收点的最小走时；第二步，根据得到的最小走时，反向追踪射线路径.但是该算法在向前计算走时时，没有考虑射线的逆向传播，不能追踪回波［9］，进而影响射线追踪的精度.针对这一问题，不少学者提出了改进算法，张东等人提出了循环计算LTI改进算法［7］，王浩全提出了交叉扫描LTI改进算法［8］，黄靓等人提出了扩张－收缩扫描改进算法［9－10］.在考虑射线的逆向传播时，文献［7－8］采用从发射点所在列向外逐列逆向扫描的方式，文献［9－10］则采用从边界列（行）向发射点所在列（行）进行逆向扫描的方式.在考虑射线的逆向传播上，文献［9－10］提出的改进算法更为合理，但是该算法存在计算效率低、收敛速度慢的问题.因为在LTI算法中，接收点要得到其最小走时，需其理论射线路径与模型单元边界所有相交节点都已得到其最小走时，这些相交节点的最小走时可能来自列扫描，也可能来自行扫描.根据本文作者对文献［5，9－10］算法步骤的理解以及对文献［9－10］中数值算例的研究，发现文献［9－10］不仅增加了从边界进行逆向传播的射线，即文献［9－10］中的收缩扫描过程，而且将文献［5］中既考虑邻列也考虑来自邻行入射射线的逐列扫描方式改为逐列（行）扫描只考虑邻列（行）入射射线的方式，这点可由文献［9］算例模型1的计算结果得出：若文献［9］的逐列（行）扫描过程既考虑邻列又考虑邻行的入射射线，则模型1经一次扩张扫描即得接收点的最小走时，此时相对误差应为0.247‰，而不应为按行列分开扫描方式经一次扩张扫描得到的计算结果3.3‰①.文献［9］采用行列分开扫描的方式使得文献［9－10］在能正确追踪回波的同时，接收点为获取其最小走时也进行了较多的无效扫描，降低了算法的计算效率.由以上分析可知，若采用行列交叉扫描的方式对扩张收缩扫描算法进行改进，将有效提高扩张收缩扫描算法的计算效率.同时，由于绕射波以及回波的存在，文献［5］算法步骤3对每列都进行水平边界节点最小走时搜索，其意义并不明确.此外，由于扩张收缩扫描算法相比文献［5］增加了收缩扫描过程，在逐列扫描过程中能够考虑上行或下行首波的最小走时，所以文献［5］中的逐行扫描可以省略.最后，由于存在收缩扫描过程，在计算竖向边界各节点最小走时时，若按文献［5］算法步骤5的计算方法，则存在着重复无效扫描.因此，也有必要对文献［5］的扩张扫描过程进行简化和改进.基于此，本文在前人研究的基础上提出了基于交叉扫描方式的扩张－收缩扫描新算法，该算法以扩张－收缩扫描算法为基础，结合文献［5］在逐列扫描过程中进行交叉扫描的思想，改进扩张－收缩扫描算法中逐列扫描的具体计算方法，提高了算法的计算效率，减少了迭代次数；同时，保留了原扩张收缩扫描算法的所有优点.1 LTI射线追踪算法的基本方程如图1所示，射线通过单元下边界的AB节段到达C节点，射线与AB的交点为D.A点及B点走时分别为TA和TB，节段AB长为L，单元慢度为s（速度的倒数），C点距A点的水平及竖向距离均为已知，分别为x，y.建立以A点为原点的局部坐标系，确定A，B，C，D点的局部坐标.现推导C点走时以及交点D 距A 点长度r的计算公式［5－6］.由线性追踪算法的基本假设可得D点的走时：根据D点的走时，结合单元慢度以及各点的局部坐标等条件，可得C点走时：将式（1）代入式（2）得根据费马原理，TC对r的一阶偏导数应当满足等于零的条件，即（设ΔT＝TB－TA）当L2s2－ΔT2 ＞0时，解方程（4）可得① 注：事实上，模型1接收点在按行列分开扫描的方式下，经扩张－收缩－再扩张后能够得到其最小走时，文献［9］给出的是一次扩张扫描结果，这应该是文献［9］作者的疏忽.若r≥0且r≤L，则若r＜0或r＞L，则计算r＝0和r＝L时的TC值，并取两者较小值作为最终TC值. 当L2s2－ΔT2≤0时，TC的计算方法与r＜0或r＞L时的情况一样.图1 经过节段AB到达C点的射线路径图Fig.1 The graph of ray path from segment ABto C在求得r值后，可以通过A点在整体坐标系中的值推出D点在整体坐标系中的值.在这里，定义D点为C点的次级源.2 扩张收缩扫描算法存在的问题及改进为说明扩张收缩扫描算法［9－10］存在的问题，以图2模型为例，模型尺寸为3m×5m，单元边长及速度分别为1m和3 500m／s；单元边界划分为2个节段；发射点S位于2号单元下边界的中点，接收点为a；线段为a节点的理论射线路径，其与模型各单元边界的交点分别为d，c和b.由理论射线路径可以看出，节点d，b只能通过扩张扫描的行扫描得到其最小走时，而节点a，c则需通过扩张扫描的列扫描得到其最小走时，所以接收点a要得到其最小走时，则模型的第3列必须先完成一次行列扫描，再对第4列进行一次列扫描.而扩张收缩扫描算法［9－10］在其扩张过程采用行列分开扫描的方式，所以a 节点在第一次扩张扫描中不能得到其最小走时.扩张－收缩扫描算法的第一次扩张扫描过程如图2所示（假定发射点所在列各节点的走时计算已经完成），其中图2（a）为扩张扫描的列扫描结果，图2（b）为完成所有列扫描后再进行行扫描的结果.从图2（b）可以看出a节点没有得到其最小走时.对于a节点，此次扩张扫描为无效扫描.同时，从图2（b）还可以看出，除a节点外还有其它节点也进行了无效扫描.易知，当模型网格划分得更多时，将有更多节点经历无效扫描.为了减少无效扫描，提高算法的计算效率.本文采用行列交叉扫描的方式对算法［9－10］进行改进，即在扩张（收缩）扫描过程，行列扫描不再分开进行，而是在其逐列（行）扫描时就同时进行行和列的扫描.仍以图2所示模型为例，采用改进后的算法对其进行一次扩张扫描，同样假定发射点所在列各节点的走时计算已经完成.图2 扩张收缩扫描算法的扩张扫描过程及结果Fig.2 Process and results of expansion scan of expansion－contraction scanning algorithm图3（a）为改进算法在扩张扫描时对第3列进行交叉扫描后的结果，可看出，b 和c节点在这次行列交叉扫描过程中都得到了最小走时，进而在第4列进行交叉扫描后，a节点能得到其最小走时（图3（b））.图3 改进算法的扩张扫描过程及结果Fig.3 Scanning process and results of improved algorithm对比图2（b）和图3（b），可以看出，改进算法进行的无效扫描更少，具有更高的计算效率，同时，由于改进算法仅改变了扫描顺序，两种算法的单次扩张扫描计算量不变.同时，对比文献［5］的一次扩张扫描过程，可以看出，本文提出的交叉扫描算法仅保留了逐列扫描过程中的行扫描，去掉了逐列扫描后附加的逐行扫描过程，并对文献［5］的逐列扫描过程进行了简化.需要说明的是，本文的交叉扫描与扩张收缩扫描在概念上是不同的，前者针对的是一次扩张收缩扫描中行列的计算次序，而后者针对的是整个单元网格的计算次序，一个为微观操作，另一个为宏观操作.下文中，除“扩张扫描”和“收缩扫描”中的“扫描”为宏观操作外，其余“扫描”均指微观操作.3 改进算法的基本步骤改进算法仍分两步进行，第一步向前计算走时，第二步，向后追踪射线路径，其中第二步与文献［9－10］相同，不赘述.改进算法的基本步骤如下：1）如图4所示，模型为3×3的网格，发射点S位于第2行第2列.以发射点所在的行和列为轴，将网格划分为四个象限，其中发射点所在行和列为各象限共有部分.图4（a）和图4（b）中的阴影部分分别为模型的第一和第二象限.图4 模型网格及象限划分示意图Fig.4 Model grid and quadrant division schematic diagram2）计算发射点S所在单元边界上各节点的走时，并记录次级源.假定此模型在各单元的边界上均只划分两个节段（图5），则各单元均有8个节点.根据单元节点以及发射点在整体坐标系中的坐标，可求得发射点S所在单元边界上各节点的走时，并将S点记为各节点的次级源.3）计算发射点S所在列所有节点的走时并记录次级源.首先计算图6中与发射点单元上边界相连的EJLG单元各节点的走时.以I点的走时计算为例（仅考虑通过下边界GE到达I节点的射线），易知，满足I 点最小走时要求的射线可能来自GE中的任一节段.此时，根据第1节给出的计算公式和计算方法分别计算出射线通过GF节段和FE节段时I点的最小走时，取两个最小走时的较小值作为I点的最小走时，并记录相应的次级源.以同样步骤求出EJLG单元其它节点的走时，完成该单元的计算.这种通过单元下边界节点走时计算其它节点走时的过程，称为向上扫描.其它通过上边界、左边界、右边界的情况分别称为向下扫描、向右扫描以及向左扫描.图5 改进算法的基本步骤2Fig.5 The basic steps 2of improved algorithm图6 改进算法的基本步骤3Fig.6 The basic steps 3of improved algorithm然后，以同样的方式向上逐个计算其它单元的节点走时，并记录相应的次级源.最后，从发射点单元开始，向下扫描，最终完成发射点S所在列所有节点走时计算以及次级源记录.4）扩张扫描.对步骤1划分的四个象限采用不同的扫描策略，具体如下.象限Ⅰ：从靠近发射点所在列的右侧列开始，向右逐列进行向右向上交叉扫描，直至模型的最后一列.以靠近发射点所在列的右侧列为例说明交叉扫描算法：如图7所示，先对此列单元进行向右扫描，然后再向上扫描，在扫描过程中，若扫描得到的节点走时比原节点走时小，则更新此节点的走时，并记录对应的次级源，否则，原节点走时以及相应次级源保持不变.象限Ⅱ，Ⅲ，Ⅳ的扫描也均从发射点所在列开始，向左（象限Ⅱ，Ⅲ）或向右（象限Ⅳ）逐列进行向左向上（象限Ⅱ）、向左向下（象限Ⅲ）以及向右向下（象限Ⅳ）交叉扫描，直至模型的边界.图7 改进算法的基本步骤4Fig.7 The basic steps 4of improved algorithm5）收缩扫描.四个象限的收缩扫描均从模型的竖向边界开始，向左（象限Ⅰ，Ⅳ）或向右（象限Ⅱ，Ⅲ）逐列进行向左向下（象限Ⅰ）、向左向上（象限Ⅳ）、向右向下（象限Ⅱ）或向右向上（象限Ⅲ）交叉扫描，直至发射点所在列.6）重复步骤4）和5），直到模型中所有节点的走时均不再改变为止.至此得到所有节点的最小走时以及相应的次级源.对比文献［5］的逐列扫描过程，本文省去了对每列都进行水平边界节点最小走时的搜索，原因如前言所述，在绕射波以及回波存在的情况下，搜索水平边界节点最小走时的意义并不明确.同时，为了确定逐列扫描过程中行扫描的顺序，本文参考了均匀介质模型各节点的理论射线路径，如发射点右上方单元各节点的最小走时均来自各自单元的左边界或下边界，据此将模型划分为四个象限，这也与仅改变文献［9］行列扫描顺序的结果一致.最后，本文按象限划分扫描区域进行逐列扫描的方式较文献［5］的逐列扫描过程更简洁，程序编制也更为容易.4 数值算例为了验证本文提出的改进算法在计算效率上的优越性，给出三个数值算例.考虑到改进算法仅改变了扫描顺序，不改变单次扩张收缩扫描计算量，故选取计算收敛所需迭代次数为对比参数.算例1 如图8，模型尺寸为3m×5m，各单元边长及速度分别为1m和3 500m／s，单元边界划分为2个节段，发射点S位于模型的下边界，且距模型左边界1.5m.采用扩张－收缩扫描算法［9－10］及本文改进算法进行计算.结果显示：扩张－收缩扫描算法需进行4次迭代才收敛，而改进算法只需要2次迭代即收敛.图8 算例1计算模型Fig.8 Calculation model of example 1算例2 将算例1模型细划为30×50的网格，单元边长为0.1m.如图9所示，发射点位置不变.单元边界划分为2个节段.计算结果显示：扩张收缩扫描算法需进行31次迭代，而改进算法仍只需2次迭代.图9 算例2计算模型Fig.9 Calculation model of example 2现将1 500号单元在两种算法下的扫描过程及结果列于表1，其中，由于算法收敛的条件，原算法及改进算法的最后一次迭代结果均与前一次的迭代结果相同.表1 1 500号单元在两种算法下的扫描过程及结果Tab.1 Scanning process and results of No.1 500under the two algorithms?算例1及算例2的计算结果表明：对于均质模型，改进算法相比于原算法具有更快的收敛速度，且当模型网格尺寸划分较细时，改进算法在计算效率上的优势更显著.算例3 如图10，模型尺寸为20m×20m，单元边长为1m；单元速度有两种，其中白色单元的速度为3 500m／s，黑色单元的速度为100m／s，单元边界划分为2个节段，发射点S和接收点R分别位于模型的左下角和下边界距发射点6m位置.计算结果显示：改进算法只需7次迭代即收敛，而原算法则需17次迭代；同时，接收点R在两种算法下的计算最小走时均为8.439 6ms，与理论最小走时相差0.19%.图10 算例3计算模型Fig.10 Calculation model of example 3算例3的计算结果表明：改进算法保留了扩张收缩扫描算法能正确处理射线逆向传播的优点，并且具有更快的收敛速度.5 结语理论分析以及数值算例表明，本文提出的基于交叉扫描方式的扩张－收缩扫描算法，不仅具有原扩张－收缩扫描算法的所有优点，而且在不增加单次扩张－扫描计算量的前提下，通过改变算法的扫描顺序，提高了算法的计算效率，加快了算法的收敛速度.特别当模型网格划分数较多时，改进算法在计算效率方面的优势更为显著.参考文献［1］VIDALE J.Finite－difference calculationoftravel times［J］.Bulletin of the Seismological Society of America，1988，78（6）：2062－2076.［2］QIN F，LUO Y，OLSEN K B，et al.Finite－differencesolution of the eikonal equation alongexpandingwavefronts［J］.Geophysics，1992，57（3）：478－487.［3］MOSERT J.Shortest path calculation of seismic rays［J］.Geophysics，1991，56（1）：59－67.［4］刘洪，孟凡林，李幼铭.计算最小走时和射线路径的界面网全局方法［J］.地球物理学报，1995，38（6）：823－832.LIU Hong，MENG Fan－lin，LI You －ming.The interface grid method for seeking global minimum travel－time and the correspondent ray path ［J］，Acta Geophysica Sinica，1995，38（6）：823－832.（In Chinese）［5］ASAKAWA E，KAWANAKA T.Seismic ray tracing using linear traveltimeinterpolation［J］.Geophysical Prospecting，1993，41（1）：99－111.［6］赵改善，郝守玲，杨尔皓，等.基于旅行时线性插值的地震射线追踪算法［J］.石油物探，1998，37（2）：14－24.ZHAO Gai－shan，HAO Shou－lin，YANG Er－hao，et al.Seismic ray tracing algorithm based on the linear traveltimeinterpolation［J］.Geophysical Prospecting for Petroleum，1998，37（2）：14－24.（In Chinese）［7］张东，谢宝莲，杨艳，等.一种改进的线性走时插值射线追踪算法［J］.地球物理学报，2009，52（1）：200－205.ZHANG Dong，XIE Bao－lian，YANG Yan，et al.A ray tracing method based on improved linear traveltimeinterpolation［J］.Chinese Journal of Geophysics，2009，52（1）：200－205.（In Chinese）［8］WANG Hao－quan.An improved method of linear travel－time interpolation ray tracing algorithm［J］.Acta Physica Polonica－Series A，2010，118（4）：521.［9］黄靓，黄政宇.线性插值射线追踪的改进方法［J］.湘潭大学学报：自然科学版，2002，24（4）：105－108.HUANG Liang，HUANG Zheng－yu.An improved method of linear interpolation ray tracing［J］.Journal of Xiangtan University：Natural Science Edition，2002，24（4）：105－108.（In Chinese）［10］黄靓.混凝土超声波层析成像的理论方法和试验研究［D］.长沙：湖南大学土木工程学院，2008：33－35.HUANG Liang.Methodology and experiment research on concrete ultrasonic computerized tomography［D］.Changsha：College of Civil Engineering，Hunan University，2008：33－35.（InChinese）。

起伏地表下基于抛物插值的最短路径射线追踪桑运云;孙军晓;焦淑萍;金艳萍;陈院生【期刊名称】《石油物探》【年(卷),期】2014(053)002【摘要】基于射线理论的初至波走时层析静校正是解决复杂地表静校正问题的实用方法,其核心技术是射线追踪.基于抛物插值的最短路径射线追踪方法利用最短路径、动态网格和抛物旅行时插值实现了初至波射线路径的追踪和初至走时的准确求取,但起伏地表条件下起伏界面的存在给射线追踪带来了许多问题.为此,对基于抛物插值的最短路径射线追踪算法进行了改进,将起伏地表及地下界面精确地投影到部分网格节点上,实现了起伏地表下的初至波路径的射线追踪.通过模型和实际资料试算,证明了改进后的算法适用于起伏地表条件下的近地表层析速度建模和静校正量计算.【总页数】7页(P142-148)【作者】桑运云;孙军晓;焦淑萍;金艳萍;陈院生【作者单位】中国石油天然气集团公司东方地球物理勘探有限责任公司研究院华北分院,河北任丘062552;中国石油天然气集团公司东方地球物理勘探有限责任公司研究院华北分院,河北任丘062552;中国石油天然气集团公司东方地球物理勘探有限责任公司研究院华北分院,河北任丘062552;中国石油天然气集团公司东方地球物理勘探有限责任公司研究院处理中心,河北涿州072751;中国石油天然气集团公司东方地球物理勘探有限责任公司研究院华北分院,河北任丘062552【正文语种】中文【中图分类】P631.4【相关文献】1.起伏地表地震波旅行时混合网格线性插值射线追踪计算方法 [J], 王琦;朱盼;叶佩;李勤;李庆春2.抛物旅行时插值最短路径射线追踪 [J], 桑运云;李振春;张凯3.起伏地表下初值射线追踪的实现 [J], 岳玉波;孙建国;杨昊;韩复兴4.起伏地表条件下各向异性地震波最短路径射线追踪 [J], 赵后越;张美根5.对旅行时进行抛物型插值的地震射线追踪方法 [J], 张赛民;周竹生;陈灵君;周惠群因版权原因，仅展示原文概要，查看原文内容请购买。

复杂地表初至波层析反演静校止应用研究[摘要]常规的方法在复杂近地表条件下很难求得正确的低速带速度和静校正量。

从理论模型上得出层析反演静校正通过非线性算法反演出准确的近地表速度模型，求出准确的静校正量。

通过本研究结果可看出层析反演静校正明显改善了复杂地区（沙漠、山地、黄土塬等）地剖面的质量，证明其是解决复杂地表地区静校正问题的一条比较有效的途径。

[关键词]复杂地表层析反演静校正中图分类号：tf046.6 文献标识码：a 文章编号：1009-914x （2013）11-0162-011 引言随着勘探的不断深入，地震勘探由平原地区逐步转移到山地、沙漠及黄土塬等地区，而这些探区地表条件极为复杂，静校正问题尤为严重。

这些地区地形起伏大，表层岩性变化非常剧烈，低降速带厚度变化大，激发和接收条件复杂，近地表条件纵、横向千差万别，导致地震反射资料不能准确成像，也造成地下构造发生扭曲。

静校正问题是制约地表复杂地区油气勘探的关键性问题。

本次研究通过辽河外围沙漠、山地地区和黄土塬等地震资料特点，选定初至波层析反演静校正的方法，解决地表复杂地区地震资料的静校正问题。

2 层析反演静校正方法层析静校正技术是一种利用单炮初至进行近地表速度反演的方法。

在层析反演中，将地质模型假设由速度单元组成，每个单元是常速，单元之间的速度不同。

首先给定一个初始的速度模型，通过射线追踪计算初至时间，它与实际旅行时的差被用来计算速度模型的修正量，模型修改后，再计算基于新的速度模型的初至旅行时，最终构成了一个迭代过程。

当正演旅行时和实际初至时间之差小于某个阈值时，就得到了最终的速度分布。

层析反演采用正反演迭代的方法，可根据初至时间重构速度场，并充分利用更多的地球物理信息（如直达波、回折波或折射波等初至信息）得到可靠的结果。

当低降速带厚度、速度横向变化非常剧烈，或者纵向速度分布异常时，层析静校正是一种十分有效的方法。

层析反演静校正技术不需要输入风化层速度，它通过运算可以较准确的得到近地表风化层的速度，另外，层析静校正可以让整个偏移距范围内的所有初至参与计算，无需人工将初至分组，所需参数量较少。

吸收介质中射线路径,走时与振幅计算——基于一阶扰动方法姜盼盼;孙建国;石秀林;徐杨杨【摘要】实际地层由于波前扩散、地层吸收和散射等原因发生衰减,为得到更加详细的地下信息和分辨率更高的图像,基于波前构建法的射线追踪,实现了吸收介质中对射线路径,走时与振幅的计算.计算表明吸收介质中的波动属性是复值;当Q为5～20时,衰减是可见的,明显的;而当Q＞30时,衰减微弱,为弱衰减介质.对吸收介质中的射线追踪进行计算,可提高地震波的分辨率,得到更清晰的地震记录.【期刊名称】《世界地质》【年(卷),期】2016(035)001【总页数】6页(P250-255)【关键词】衰减;复值;波前构建;吸收介质;射线追踪【作者】姜盼盼;孙建国;石秀林;徐杨杨【作者单位】吉林大学地球探测科学与技术学院,长春130026;国土资源部应用地球物理综合解释理论实验室,长春130026;吉林大学地球探测科学与技术学院,长春130026;国土资源部应用地球物理综合解释理论实验室,长春130026;吉林大学地球探测科学与技术学院,长春130026;国土资源部应用地球物理综合解释理论实验室,长春130026;吉林大学地球探测科学与技术学院,长春130026;国土资源部应用地球物理综合解释理论实验室,长春130026【正文语种】中文【中图分类】P631.44随着中国油气工业的发展，油气勘探难度越来越大。

地震勘探作为现阶段油气勘探的主要方法，其油气勘探效果主要取决于人们在方法理论研究中所用的介质模型与实际地层的符合程度。

在常规的地震资料处理中，一般把地层介质看成是均匀和完全弹性，忽略了地震波在地层中传播时的衰减效应。

在这样的假设下，理论上的研究结果与地震波在地层中传播时观测到的实际记录在频率及振幅上有很大的差别[1]。

研究表明实际地层更接近于黏弹性介质，由于波前扩散、地层吸收和散射等原因，在其中传播的地震波会被吸收而导致能量衰减，频率越高，衰减越严重。