一种多尺度伪弯曲射线追踪方法技术

本发明专利技术公开了一种多尺度伪弯曲射线追踪方法，涉及地球物理勘探技术领域。该方法包括：通过将初始射线路径的所有相邻两节点间的射线路径长度一一与提取的对应尺度的三维速度模型分量的立方体网格的边长进行比较；当该提取的对应尺度的速度模型分量的立方体网格的边长小于相邻两节点间射线路径长度时，将该相邻两节点的中间点作为新节点；采用伪弯曲射线追踪方法对该新节点进行修改，将初始射线路径的两节点和修改后的节点之间的连线作为优化后的射线路径，将该优化后的射线路径作为初始射线路径。通过本发明专利技术的方法使伪弯曲法从技术上进行改进，使其适用于包括速度突变的任意复杂介质中，以充分发挥这种方法本身所具有的高效率高精度的优点。

【技术实现步骤摘要】
一种多尺度伪弯曲射线追踪方法
本专利技术涉及地球物理勘探
，更具体的涉及一种多尺度伪弯曲射线追踪方法。
技术介绍
射线追踪作为一种快速有效的正演模拟技术，在野外数据采集设计、波场分析与识别、静校正、层析成像以及偏移成像等方面都有重要的应用。射线追踪是已知介质和观测系统，求取震源和接收点之间波传播的路径、走时、振幅等问题。传统的射线追踪技术可分为：试射法(或打靶法)、弯曲法和波前重建法等。众多射线追踪方法中弯曲法具有迭代速度快，精度也较高，但不适于连续介质，而伪弯曲法射线追踪考虑了速度梯度，因而能适用于连续速度变化的介质中，却不适应于存在速度突变的介质中。
技术实现思路
本专利技术实施例提供一种多尺度伪弯曲射线追踪方法，用以解决现有技术中伪弯曲法射线追踪不适用于速度突变介质中的问题。本专利技术实施例提供一种多尺度伪弯曲射线追踪方法，包括：步骤1、将预定的第一三维速度模型按立方体网格剖分，进行空间离散化；步骤2、根据预定的炮点坐标和接收点坐标在剖分后的第一三维速度模型中确定初始射线路径；步骤3、设置初始尺度因子k＝1；步骤4、根据公式(1)，从所述剖分后的第一三维速度模型中提取对应尺度的三维速度模型分量Vk；其中，公式(1)中，L为高斯滤波器长度，(x0，y0,z0,)是高斯窗中心点，k为尺度因子，W(k)为高斯滤波器，Vk为高斯滤波器与三维速度模型的卷积结果，V为三维速度模型；步骤5、将初始射线路径的所有相邻两节点间的射线路径长度一一与提取的对应尺度的三维速度模型分量的立方体网格的边长进行比较；步骤6、当所述提取的对应尺度的速度模型分量的立方体网格的边长小于相邻两节点间射线路径长度时，将所述相邻两节点的中间点作为新节点；步骤7、采用伪弯曲射线追踪方法对所述新节点进行修改，将初始射线路径的两节点和修改后的节点之间的连线作为优化后的射线路径，将所述优化后的射线路径作为初始射线路径；步骤8、令尺度因子k加1，并当尺度因子k对应的网格化后的三维速度模型与预定的第二三维速度模型不匹配时，重复步骤4-7；步骤9、输出最终的射线路径，并计算射线走时；其中，所述最终的射线路径为当尺度因子k对应的网格化后的三维速度模型与预定的第二三维速度模型匹配时的初始射线路径。较佳地，所述采用伪弯曲射线追踪方法对所述新节点进行修改具体包括：S1、假设初始射线路径的两节点为Pk-1和Pk+1，Pk-1和Pk+1的连线的中心处插入的新节点为Pmid；S2、通过以下迭代公式组对中间节点Pmid进行修改，修改后的节点为Pk，即优化后的路径节点为Pk-1、Pk和Pk+1；其中，Vmid为中点Pmid处的速度值，为修改后点Pk处的速度梯度向量，为中点Pmid处的速度梯度向量，xk-1为节点Pk-1的横坐标，xk+1为节点Pk+1的横坐标，Γ为方向矢量，n为单位方向矢量，R为中点Pmid到修改后点Pk的距离，c为点Pk-1和点Pk+1处慢度的平均，L＝|xk+1-xmid|，L为节点Pk+1和节点Pmid之间的距离，xmid为节点Pmid的横坐标，xk为节点pk的横坐标。较佳地，所述中点Pmid处的速度值Vmid可根据如下公式(2)计算得到；其中，公式(2)式中，Vmid(x,y,z)为中点Pmid(x,y,z)处的速度，Vmid((i+l,j+m,k+n)表示中点Pmid周围8个网格节点处的速度值；(x,y,z)是中点Pmid的坐标，(xi,yj,zk)分别是网格节点的x坐标，y坐标和z坐标。本专利技术实施例中，通过判断提取的对应尺度的三维速度模型分量的网格大小是否小于初始射线路径；当所述提取的对应尺度的速度模型分量的网格大小小于初始射线路径时，将初始射线路径的中心处插入新节点；采用伪弯曲射线追踪方法对所述新节点进行修改，将初始射线路径的两节点和修改后的节点之间的连线作为优化后的射线路径，将所述优化后的射线路径作为初始射线路径，也即，通过逐级提取，并且比较每一步提取的对应尺度的三维速度模型分量的网格大小与初始射线路径后，对初始射线路径进行加密和优化，直至路径满足设定的终止条件，从而使得伪弯曲法从技术上进行改进，使其适用于包括速度突变的任意复杂介质中，以充分发挥这种方法本身所具有的高效率高精度的优点。附图说明图1为本专利技术实施例提供的一种多尺度伪弯曲射线追踪方法的流程示意图；图2为本专利技术网格剖分示意图；图3为伪弯曲法三点迭代示意图图4为复杂速度模型多尺度伪弯曲法射线追踪过程。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。多尺度伪弯曲射线追踪方法属于地球物理勘探中的正演问题，即通过给定三维速度模型和炮点、接收点排列的位置来得到正确的地震波传播的射线路径，进而得到正确的地震记录。为实现上述目的，本专利技术实施例提供的一种多尺度伪弯曲射线追踪方法的流程示意图，如图1所示，该方法包括：步骤1、将预定的第一三维速度模型按立方体网格剖分，进行空间离散化。其中，按照立方体网格对第一三维速度模型进行剖分，如图1所示。步骤2、根据预定的炮点坐标和接收点坐标在剖分后的第一三维速度模型中确定初始射线路径。步骤3、设置初始尺度因子k＝1。步骤4、根据公式(1)，从该剖分后的第一三维速度模型中提取对应尺度的三维速度模型分量Vk。其中，公式(1)中，L为高斯滤波器长度，(x0，y0,z0,)是高斯窗中心点，k为尺度因子，W(k)为高斯滤波器，Vk为高斯滤波器与三维速度模型的卷积结果，V为三维速度模型。步骤5、将初始射线路径的所有相邻两节点间的射线路径长度一一与提取的对应尺度的三维速度模型分量的立方体网格的边长进行比较。步骤6、当该提取的对应尺度的速度模型分量的立方体网格的边长小于相邻两节点间射线路径长度时，将该相邻两节点的中间点作为新节点。步骤7、采用伪弯曲射线追踪方法对该新节点进行修改，将初始射线路径的两节点和修改后的节点之间的连线作为优化后的射线路径，将该优化后的射线路径作为初始射线路径。其中，该采用伪弯曲射线追踪方法对该新节点进行修改具体包括：S1、假设初始射线路径的两节点为Pk-1和Pk+1，Pk-1和Pk+1的连线的中心处插入的新节点为Pmid。S2、通过以下迭代公式组对中间节点Pmid进行修改，修改后的节点为Pk，即优化后的路径节点为Pk-1、Pk和Pk+1。其中，Vmid为中点Pmid处的速度值，为修改后点Pk处的速度梯度向量，为中点Pmid处的速度梯度向量，xk-1为节点Pk-1的横坐标，xk+1为节点Pk+1的横坐标，Γ为方向矢量，n为单位方向矢量，R为中点Pmid到修改后点Pk的距离，c为点Pk-1和点Pk+1处慢度的平均，L＝|xk+1-xmid|，L为节点Pk+1和节点Pmid之间的距离，xmid为节点Pmid的横坐标，xk为节点Pk的横坐标。另外，插入新节点的方式也称为加密射线节点，如图3所示，假设初始射线路径的两节点为Pk-1和Pk+1，在Pk-1和Pk+1的连线的中心处插入的新节点为Pmid，即通过增加射线中间节点，得到加密后的射本文档来自技高网...

【技术保护点】
一种多尺度伪弯曲射线追踪方法，其特征在于，包括：步骤1、将预定的第一三维速度模型按立方体网格剖分，进行空间离散化；步骤2、根据预定的炮点坐标和接收点坐标在剖分后的第一三维速度模型中确定初始射线路径；步骤3、设置初始尺度因子k＝1；步骤4、根据公式(1)，从所述剖分后的第一三维速度模型中提取对应尺度的三维速度模型分量V

【技术特征摘要】
1.一种多尺度伪弯曲射线追踪方法，其特征在于，包括：步骤1、将预定的第一三维速度模型按立方体网格剖分，进行空间离散化；步骤2、根据预定的炮点坐标和接收点坐标在剖分后的第一三维速度模型中确定初始射线路径；步骤3、设置初始尺度因子k＝1；步骤4、根据公式(1)，从所述剖分后的第一三维速度模型中提取对应尺度的三维速度模型分量Vk；其中，公式(1)中，L为高斯滤波器长度，(x0，y0,z0,)是高斯窗中心点，k为尺度因子，W(k)为高斯滤波器，Vk为高斯滤波器与三维速度模型的卷积结果，V为三维速度模型；步骤5、将初始射线路径的所有相邻两节点间的射线路径长度一一与提取的对应尺度的三维速度模型分量的立方体网格的边长进行比较；步骤6、当所述提取的对应尺度的速度模型分量的立方体网格的边长小于相邻两节点间射线路径长度时，将所述相邻两节点的中间点作为新节点；步骤7、采用伪弯曲射线追踪方法对所述新节点进行修改，将初始射线路径的两节点和修改后的节点之间的连线作为优化后的射线路径，将所述优化后的射线路径作为初始射线路径；步骤8、令尺度因子k加1，并当尺度因子k对应的网格化后的三维速度模型与预定的第二三维速度模型不匹配时，重复步骤4-7；步骤9、输出最终的射线路径，并计算射线走时；其中，所述最终的射线路径为当尺度因子k对应的网格化后的三维速度模型与预定的第二三维速度模型匹配时的初始射线路径。2.如权利要求1所述的多尺度伪弯曲射线追踪方法，其特征在于，所述采用伪弯曲射线追踪方法对所述新节点进行修改具体包括：S1、假设初始射线路径的两节点为Pk-1和Pk+1，Pk-1和Pk+1的连线的中心处插入的新节点为Pmid；S2、通过以下迭代公式...

【专利技术属性】
技术研发人员：张唤兰，
申请(专利权)人：西安科技大学，
类型：发明
国别省市：陕西,61