基于分层介质的地震走时参数化网格模型构建方法及应用

技术领域

本发明涉及禁核试核查监测技术领域，具体涉及一种基于分层介 质的地震走时参数化网格模型构建方法及应用。

背景技术

走时定量描述了波在介质中的传播时间。事件源产生的地震波， 经地壳、地幔等介质传播到台站，反映在地震台站探测设备生成的连 续波形数据中。走时估计是地震数据处理和事件判断的重要基础，主 要可分为基于高频近似射线理论的最短路径方法、基于程函方程的数 值解方法、基于惠更斯原理的波前构建法、基于频率域波动方程方法 等。这些研究方法在计算精度、稳定性、对复杂介质的适应性上各有 权衡，但普遍存在计算效率低的问题。鉴于到走时运算是波形数据处 理中频繁运算，这些研究的运算效率难以满足监测实时处理需求，通 常面向事件后精细分析、重定位等处理。

基于射线理论，可能事件源k与台站α

和全球大多数监测机构一样，考虑数据处理的实时性需求，我国 禁核试核查国家数据中心(NDC)采用主流的IASP91简化模型，导 致走时估计普遍存在偏差，甚至部分台站对敏感地区存在严重偏差。 一些近期研究引入了对介质异构性的考虑，实际应用表明能有效降低 走时估计偏差。但此类方法的运算效率通常面向事件后精细分析、重 定位等处理。高性能地给出低偏差的走时估计，是改进实时监测应用 场景中事件定位的迫切需求。

发明内容

(一)要解决的技术问题

本发明提供一种地震走时参数化建模方法及应用，针对以下问题 的改进：

(1)现有考虑介质特性的精确地震走时运算实现复杂度高，难 以适用于实时监测处理；

(2)三维介质模型难以应用于实时监测处理场景；

(3)如何使用位置已知的标准事件修正参数化模型；

(4)集成走时参数化模型到实时监测应用场景。

(二)技术方案

为了解决上述技术问题，本发明提供了一种基于分层介质的地震 走时参数化网格模型构建方法，包括以下步骤：

(1)基于特定的水平分层介质模型，离线建立全球网格剖分； 其中采用基于矩形、三角镶嵌的均匀网格剖分，以及使用Voronoi图 这些不规则形状的非均匀网格剖分的方法之一进行全球网格剖分；

(2)离线参数化各个网格节点，建立特定于位置固定监测台站 的节点走时估计值；

其中使用走时估计方法，通过预先估算节点到特定台站的走时， 并参数化节点位置存储，为每个节点关联了一个数组数据结构，每个 数组保存一维的数据值，规模对应于监测台网中的固定台站数量，并 支持多种数值类型存储；

(3)离线建立网格节点索引；

其中通过索引将数据结构关联到三维空间节点，集成空间索引， 支持从空间位置、深度、震相、台站多维度对百万数量级的节点关联 数据进行查询和修改；

(4)提供自然邻点插值算法实现方法，针对非线性最小二乘反 演定位这些对偏导连续性有要求的应用场景，在水平和垂向上提供核 函数连续可微的自然邻点插值函数；

(5)针对给定固定台网离线构建走时参数化网格模型，模型提 供在线的网格模型检索接口，走时估计过程中使用插值方法给出对空 间任意一点到给定台站的震相走时的估计值。

优选地，步骤1中，对水平分层介质模型确定的三维空间做网格 剖分，选定水平剖分粒度，垂向深度参考层次化介质模型的深度设定， 由此形成分层的全球网格剖分。

优选地，步骤2中，对于剖分给出的每个网格节点，按确定台网 中台站的数量，分配数组的数据结构规模，数组中的每个元素，顺次 存储节点到确定台站的走时估计，在面向定位运算时，额外顺次存储 作为射线参数的慢度信息。

优选地，步骤2中，对于落在距离台站15-20度范围内的节点， 调用SALSA3D方法计算该节点处的走时值，写入该节点处数组的相 应位置；

对于落在距离台站0-15度范围内的节点，调用RSTT方法计算 节点处的走时值，写入节点处数组的相应位置；

对于不满足RSTT计算条件的节点，调用IASP91速度表做插值 结算节点处的走时值，写入节点处数组的相应位置。

优选地，步骤3离线建立网格节点索引时，对所形成的数组集合， 建立空间索引，以便形成二进制的走时参数化网格模型，通过节点数 据访问接口，提供根据空间位置、深度、台站和震相的多维度数据检 索和修改。

优选地，步骤4中，在水平和垂向上提供核函数连续可微的自然 邻点插值函数时，插值权重使用基于空间距离的Sibson权重。

优选地，步骤5中，分别对浅源近震震相Pn、Pg、Sn和Lg震 相构建模型文件，设定模型覆盖范围是每个台站周边20度以内；

加载震相特定的二进制走时参数化网格模型文件到内存；通过数 据访问接口，获取模型内空间任一点的紧邻参数化节点数值；

做自然邻点插值运算，将给出的拟合结果作为走时估计值返回。

优选地，所述走时参数化网格模型M形式化表示为六元组：

其中，Ω是层次介质模型，描述全球三维介质特性；Θ(n)是网格 划分方法，将三维空间划分为n个网格节点；集合Α表示由m个台站 构成的台网，其中i∈{1，m}表示集合A中的特定台站；T(k，i)是走时 估计方法，计算节点k到台站i的走时；Ψ

从模型的具体实现上看，走时参数化网格模型M定义了离散节 点的空间位置k、节点上的参数向量v(k)和插值函数

M＝{v(k)＝

←Ω,i∈{1,m},k∈{1,n}}( 2)

建立特定于m个台站的参数化网格模型M，首先根据物理介质Ω 给出层次化网格剖分方法Θ(n)，形成n个节点的集合，每个节点k∈ {1,n}维护一个m维向量v(k)＝

优选地，台站i在空间位置r处的走时为：

其中，f

本发明还提供了一种所述的方法在禁核试核查监测技术领域中 的应用。

(三)有益效果

本发明提出的地震-声学走时参数化网格模型构建方法，基于分 层介质物理模型构建全球三维空间网格，整合先进算法和标准事件离 线估算各网格节点到台站的走时，进而建立台站特定的走时参数化模 型，在实现提高走时估计精度的同时，解决了实时监测中难以考虑复 杂介质等因素改进走时估计的问题。本发明的主要优点是：

(1)实现了一种适用于多种介质模型、走时算法、空间内插和 节点检索的走时参数化网格建模方法，针对定位运算引入连续可微的 自然近邻基函数支持水平和垂向内插，集成空间索引技术支持高效的 多维度网格节点检索。

(2)面对禁核试核查监测的典型应用场景，构建了针对特定台 网的浅源近震震相走时参数化模型，在满足监测数据处理实时性要求 的同时，提高了地震数据自动处理的走时运算精度，降低了定位运算 误差。

附图说明

图1为本发明的参数化网格模型的建模和应用流程图；

图2为本发明的参数化网格模型中节点处走时估计的计算流程图。

具体实施方式

为使本发明的目的、内容、和优点更加清楚，下面结合附图和实 施例，对本发明的具体实施方式作进一步详细描述。

本发明提供了一种地震走时参数化模型的建模方法，并集成应用 于实时监测处理。包括以下步骤：

(1)基于特定的水平分层介质模型，离线建立全球网格剖分。 本发明支持多种全球网格剖分方法，包括基于矩形、三角镶嵌的均匀 网格剖分，以及使用Voronoi图等不规则形状的非均匀网格剖分等。

(2)离线参数化各个网格节点，建立特定于位置固定监测台站 的节点走时估计值。本发明支持使用现有走时估计方法，通过预先估 算节点到特定台站的走时，并参数化节点位置存储。为每个节点关联 了一个数组数据结构。每个数组保存一维的数据值，规模对应于监测 台网中的固定台站数量，并支持多种数值类型存储。

(3)离线建立网格节点索引，提高模型检索效率，并优化模型 的规模。通过索引将数据结构关联到三维空间节点，集成高效的空间 索引，支持从空间位置、深度、震相、台站等多维度对百万数量级的 节点关联数据进行高效查询和修改；

(4)提供自然邻点插值算法实现方法，针对非线性最小二乘反 演定位等对偏导连续性有要求的应用场景，在水平和垂向上提供核函 数连续可微的自然邻点插值函数。

(5)实际应用中，针对给定固定台网离线构建走时参数化网格 模型，模型提供在线的网格模型检索接口，走时估计过程中使用插值 方法给出对空间任意一点到给定台站的震相走时的估计值。

其中，本发明提出一种地震走时参数化网格模型的通用形式化表 示。从模型的组成看，所述走时参数化网格模型M可形式化表示为 六元组：

其中，Ω是层次介质模型，描述全球三维介质特性；Θ(n)是网格 划分方法，将三维空间划分为n个网格节点；集合A表示由m个台站 构成的台网，其中i∈{1,m}表示集合A中的特定台站；T(k,i)是走时 估计方法，计算节点k到台站i的走时；Ψ

从模型的具体实现上看，走时参数化网格模型M定义了离散节 点的空间位置k、节点上的参数向量v(k)和插值函数

M＝{v(k)＝

←Ω,i∈{1,m},k∈{1,n}}

建立特定于m个台站的参数化网格模型M，首先根据物理介质Ω 给出层次化网格剖分方法Θ(n)，形成n个节点的集合。每个节点k∈ {1,n}维护一个m维向量v(k)＝

本发明还支持定位应用和标准事件修正的空间插值。离散网格点 模型参数化通过定义离散节点的位置、节点上的参数值和插值函数来 描述介质模型。模型中任意一点的参数值可以通过节点上的值和插值 函数获得，即台站i在空间位置r处的走时为：

其中，f

具体地，本发明实施例提供一种地震走时参数化网格模型建模方 法和应用，结构如图1，具体如下：

(1)构建走时参数化网格模型，离线计算过程为：

(1.1)建立全球网格剖分：引入特定的层次化介质模型，例如 NOAA的etopo1数据，构建全球三维空间；对介质模型确定的三维 空间做网格剖分，选定水平剖分粒度(例如在经度和维度方向均为2 度)，垂向深度参考层次化模型深度设定，由此形成分层的全球网格剖分；

(1.2)离线参数化各网格节点：对于剖分给出的每个网格节点， 按确定台网中台站的数量，分配数组的数据结构规模。数组中的每个 元素，顺次存储该节点到确定台站的走时估计。在面向定位运算时， 额外顺次存储作为射线参数的慢度信息。

(1.2.1)对于落在距离台站15-20度范围内的节点，调用 SALSA3D方法计算该节点处的走时值，写入该节点处数组的相应 位置。

(1.2.2)对于落在距离台站0-15度范围内的节点，调用RSTT 方法计算节点处的走时值，写入节点处数组的相应位置。

(1.2.3)对于不满足RSTT计算条件的节点，调用IASP91 速度表做插值结算节点处的走时值，写入节点处数组的相应位置。

(1.3)离线建立网格节点索引：对所形成的数组集合，建立空 间索引，形成二进制的走时参数化网格模型。通过节点数据访问接口， 提供根据空间位置、深度、台站和震相的多维度数据检索和修改。

(1.4)提供自然邻点插值：在水平和垂向上提供核函数连续可 微的自然邻点插值函数，插值权重使用基于空间距离的Sibson权重。

(2)参考图2，应用走时参数化网格模型，在线估计空间任一 点到特定台站的走时；

(2.1)分别对浅源近震震相Pn、Pg、Sn和Lg震相构建模型文 件，设定模型覆盖范围是每个台站周边20度以内。

(2.2)，加载震相特定的二进制走时参数化网格模型文件到内存； 通过数据访问接口，获取模型内空间任一点的紧邻参数化节点数值；

(2.3)做自然邻点插值运算，将给出的拟合结果作为走时估计 值返回。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领 域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以 做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。