[发明专利]一种混合二维地震波走时计算方法

摘要

本发明公开了一种混合二维地震波走时计算方法，包括以下步骤：读入相关参数、速度模型；沿着震源向不同方向发出射线并计算中心射线信息；使用波前构建法计算射线范围内网格点的地震波走时；对全空间网格点走时属性进行分类并以此为依据建立初始窄带；使用快速推进法计算剩余网格点地震波走时。本发明通过使用窄带技术连接了快速推进法与波前构建法，通过波前构建法提高了震源附近较小区域的地震波走时计算精度，从而提高了快速推进法在剩余区域网格节点的计算精度，实现了一种兼顾计算效率与计算精度的二维地震波走时计算。

主权项

1.一种混合二维地震波走时计算方法，其特征在于，包括以下步骤：步骤1：读入相关参数文件、速度模型，其中，所述参数文件包含速度模型的网格点数、网格间距和震源位置；步骤2：从炮点沿着不同方向发射射线并计算中心射线信息；步骤3：使用波前构建法计算射线覆盖范围内网格节点走时；步骤4：对所有网格点进行分类，将它们分为接受点、窄带点或远离点，即如果一个网格点的地震波走时已经计算过，并且它周围点的走时全部计算过，则这个点为接受点；如果一个网格点的地震波走时计算过，而它周围点的走时并非全部计算过，则这个点为窄带点；如果一个点的走时没有被计算过，则这个点为远离点；步骤5：将所有窄带点移入窄带内，构建初始窄带；步骤6：扩展窄带，直至窄带为空；在这一过程中，网格点走时是通过迎风差分法求解二维程函方程获得的，所述二维程函方程为：|▽τ|＝s其中，τ为地震波走时，s为慢度，▽为梯度符号，上述公式中梯度项的迎风差分表达式为： | ▿ τ | = max ( D i , j - x τ , 0 ) 2 + min ( D i , j + x τ , 0 ) 2 + max ( D i , j - z τ , 0 ) 2 + min ( D i , j + z τ , 0 ) 2 1 / 2 ]]>其中，

下载完整专利技术内容需要扣除积分，VIP会员可以免费下载。 免登录下载普通用户下载升级VIP会员，免费下载

该专利技术资料仅供研究查看技术是否侵权等信息，商用须获得专利权人授权。该专利全部权利属于西南交通大学，未经西南交通大学许可，擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作，请联系【客服】

本文链接：http://www.vipzhuanli.com/patent/201810077621.8/，转载请声明来源钻瓜专利网。