【技术实现步骤摘要】
一种层状介质中的变网格八角汇集射线追踪方法
[0001]本专利技术涉及射线追踪
，更具体地说，涉及一种层状介质中的变网格八角汇集射线追踪方法。
技术介绍
[0002]地震波动的数值模拟方法主要分为两类，一类是波动方程数值模拟法，另一类是基于高频近似的波场计算方法，波动方程数值模拟法包含丰富的波场信息，但方法易引起干扰波，计算速度慢对计算机要求高，射线追踪方法是高频近似波场计算方法，能够直观反映地震波的几何传播路径，计算速度快，在正演模拟、叠前深度偏移成像、层析成像、微地震监测等领域占有重要地位。传统的射线追踪方法为试射法和弯曲法，其存在计算量大、隐影区射线少、难以处理较强速度变化介质、不易找到全局最小走时等问题。近年来基于费马原理和惠更斯
‑
夫列涅尔原理，发展了程函方程有限差分法、最短路径法、慢度匹配法、走时插值法和波前重建法等方法。其中慢度匹配法计算量大，最短路径法路径不完全满足斯奈尔定律，走时插值法的走时线性假定会产生误差，导至路径不够精确，且射线逆向传播时，无法求得最小走的路径，波前重建法存在非...
【技术保护点】
【技术特征摘要】
1.一种层状介质中的变网格八角汇集射线追踪方法，其特征在于：包括以下步骤：S1：进行走时计算时，首先需要对地质模型进行数字参数化，采用的模型参数化采用分块
‑
节点模型，依据速度特征先对模型进行分块，参照图片1中分布块为分块结果，分块边界即为速度界面；S2：然后在每一个分块中用节点对分块进行参数化，每一分块的节点间距可以不同，即网格为变网格，图1中X位置为节点，节点的再次剖分，使一分块分为了很多节点组成的更小分块，每一小分块的速度用其左上角节点值标记，此标记方式利于后续的走时计算；S3：图2中A为待求走时的计算节点，B
‑
J为其周围节点，周围节点不需要考虑走时是否已知；S4：对图2中数个三角形域在考虑局域计算的情况下进行计算；S5；在计算机中计算实现。2.根据权利要求1所述的一种层状介质中的变网格八角汇集射线追踪方法，其特征在于：步骤S3还包括以下步骤：S301：不采用vadale的扩展方阵方式求解的，而是把B
‑
J所围A的区域分划为八个三角形；S302：计算A点走时时，需要同时考虑八个三角形的走时情况，图2中标记为1的箭头表示从边缘透射，运用计算公式进行计算。3.根据权利要求1所述的一种层状介质中的变网格八角汇集射线追踪方法，其特征在于：步骤S302中还包括以下步骤：S3021：计算为公式(1)为：S3022：标记为2的箭头表示从角点直达或衍射，其计算公式为(2)：S3023：标记为3的箭头表示滑行或首波，其计算公式为(3)：4.根据权利要求1所述的一种层状介质中的变网格八角汇集射线追踪方法，其特征在于：步骤S302中，三个计算公式只是表述了三角形ABC的情况，其中上式中t
A1
,t
A2
,t
A3
为三个不同射线路径所求取的走时值，A点最终走时应取三者中最小值，t
B
,t
...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。