地震波初至拾取的分形研究.pdfVIP

现代地质 第 16 卷第2 期 GEOSCIENCE Vo 1.16 No.2 2002 年 6 月 Jun. 2002 地震波初至拾取的分形研究 曾富英，李敏锋，中 维 (中国地质大学国土资源与高新技术研究中心，北京 100083) 摘要:在地震数据处理过程中，传统的方法都是以振幅、频率、相位等的变化为依据来拾取 地震初至波O 作者突破了这些传统依据，根据分形的原理，用时窗的变化来反映地震道曲线 的局部形态特征，通过一元线性回归模型统计出曲线的形态特征参数一一分数维，然后由分 数维的变化来判断初至地震波，得出了一种高效的拾取初至地震波的新方法。针对该方法的 计算过程开发了处理软件模块。理论的模拟证明该方法非常适用于目前越来越庞大而复杂的 地震勘探数据资料的处理。通过实际的试算表明软件模块的适用性强。该方法精度高，计算 速度极快，抗噪能力强，在低信噪比的情况下也能快速精确地拾取初至地震波。 关键询:分数维;地震数据处理;振幅;频率;相位 中图分类号 P631. 4 文献标识码 A 文章编号 1000 - 8527(2002)02 - 0209 - 05 Box Counting 等[2.3] 。本文介绍分形维的基本概 。 寻| 念及计算处理的过程，通过编程在计算机上实现了 在地震勘探中，确定初至波的旅行时是非常重 分维数的计算。将该方法实际应用于青海某工区 耍的，特别是在折射波和面波勘探中，通常用初至 中，取得了很好的效果。 波的旅行时来计算表层的速度和层厚等介质参数。 1 地震波初至拾取的分形维方法原理 为了提高计算过程的精度和操作效率，人们一般用 地震初至波自动拾取技术来做这项工作。到目前为 1. 1 分形维 止，已提出了许多具体的方法。但是，随着地震勘 在传统几何学(欧几里得空间)中的维数只能 探资料的越来越庞大，计算越来越复杂，这些地震 是整数，而在分形几何学中维数可以是任意正实 初至波自动拾取的方法就显得有些不足了。一般说 数[4] 。分形维针对的研究对象则是空间中的分形 来，地震初至波的拾取就是要确定地震道上纯噪音 体，即非线性系统中的不规则几何形体。其研究方 信号和噪音与地震叠加信号之间的分界时刻。通 法是描述研究对象内部的不均匀性，刻画研究对象 常，人们在人工处理过程中，都是根据地震信号的 层次结构的整体数量特征。对于二维 Euclidean 空 振幅、频率或相位等的变化来确定这一时刻;还有 间，自仿射定义为 f ( 缸， kHy) 与 f (x , y) 是 一些方法，除了比较同一道上的振幅、频率、相位 统计自相似的[5.6] 其中 H 叫作 Hausdorff 测度， 等的变化以外，还将相邻道之间的这些信号作一比 h 是比例系数。根据上面的叙述很自仿射分形体就 较。但是，所有这些方法对信噪比比较高的地震信 可以用分数维来描述其不规则程度[2] 。对于地震 号较为有效，而对低信噪比的信号，其精度就很低 道的抽样序列来说，将其作元量纲图形化处理可得 了[1] 。 到一条自仿射曲线。因此，地震信号可以用分数维 近年来，为了解决这一问题，应用了分形维方 来描述其不规则程度。而且，分数维在地震有效信 法。关于用分形维方法实现这一拾取过程，一些专 号到来前后是有变化的，这成为利用该方法判断初 家提出了一些有效的途径:如 Divider Method