无射线层析静校正方法的应用研究.docVIP

无射线层析静校正方法的应用研究 　　【摘要】无射线追踪层析静校正方法不同于典型的折射波反演。它同时利用了首波和回折波，对其进行空间旅行时分解，再使用Herglotz- Wiecher方程来建立速度模型，进而实现静校正量的计算。并将其应用于实际资料中，得到很好的效果。 　　【关键词】静校正；无射线追踪层析静校正；非显式射线追踪层析静校正 　　1. 概述 　　（1）无射线追踪层析静校正可称之为非显式射线追踪层析静校正法，或折射层析静校正法。Enright等（1992）和Ecoublet等（1995）指出，在层析中避免射线追踪的方法，该方法通过程函方程对空间旅行时进行分解，虽然能很好的适用于VSP和井间地震，但很难用于折射数据中。Zhang和Toksoz（1998）证明使用平均和视慢度旅行时反演比用观测到的旅行时进行直接反演更加稳定。随后Konstantin Osypov在1999年首次提出，它使用空间旅行时分解，但不用程函方程，而使用Herglotz- Wiecher方程来计算速度模型（Aki and Richards，1980），进而实现静校正量的计算。 　　（2）该方法不同于典型的折射波反演，它包含了首波和回折波。第一步，从拾取的初至时间中估算出每个站点的时距函数；第二步，将这些函数转换成速度-深度函数，进而组合成近地表模型；第三步，利用导出的模型计算长波长静校正量；最后，通过对剩余旅行时进行地表一致性分解算出短波长静校正量。 　　（3）其优点是：有效地结合了首波的稳定性和回折波的灵活性，不显式地进行射线追踪和模型网格化，不依赖初始模型，所得反演结果比较稳定，同时还可以利用全偏移距范围的初至信息。 　　图1积分路径 　　2. 无射线追踪层析静校正原理及参数 　　2.1无射线追踪层析静校正原理。 　　无射线追踪层析方法，将旅行时反演分为两个部分： 　　2.1.1第一步，估算视速度和延迟时。将初至旅行时划离散化为一系列偏移距子范围，对每个子范围中的旅行时ti（Rs，Rd） 进行空间分解，i是偏移距子范围索引。 　　2.1.2第二步，利用导出的视速度计算折射波穿透深度Z（VXsd，Rmp） 。这一步用Herglotz-Wiechert积分公式求解： 　　式中，VX·R 表示深度Z处的瞬时层速度，公式假设随偏移距和深度的增加，VX·R 单调递增。 Z（VXsd，Rmp）表示回折波的回折深度或首波的折射深度。所得一系列的Z值形成该位置地下速度-深度函数（即V-Z函数），对于一系列地表控制点，这些V-Z函数便可估算出近地表反演模型。 　　综上即为无射线层析反演的理论基础，它其实是折射波原理的一种推广。在反演的过程中，非显式地利用折射波原理导出速度模型，得到一个无射线追踪的层析重建目标。 　　图3高程静校正叠加剖面 　　图4无射线追踪的层析静校正叠加剖面 　　2.2影响无射线追踪层析静校正的主要参数。 　　在实际资料处理过程中，参数的选取是非常重要的。影响无射线追踪层析静校正技术的主要参数如下： 　　（1）偏移距范围。理论上，无射线追踪层析静校正适用于全偏移距的初至信息，由此能增加模型反演的稳定性。但因近偏移距折射数据常受噪声干扰或者不可用，这种情况下，可充分利用先验信息（如井口时间、测井数据、重力数据和各种地质资料）以便约束表层模型的速度。同时，受计算效率与波的穿透深度限制，利用的最大偏移距也小于实际的最大偏移距。 　　（2）反演层数。反演层数对应速度模型的层数。层数的确定也是影响成像效果的关键因素，选取的范围不能太大，也不能太小。太大时容易造成反演结果的不稳定；太小则难以充分反映地下地层的结构。通常，选择的范围小于5层。 　　（3）风化层速度。通常，可用直达波来估算。不过，直达波在近偏移距处，当近偏移距信息受干扰较大或不可用时，则需提供风化层速度，或者直接利用微测井或小折射数据。 　　（4）速度模型的平滑距离。该参数影响反演结果和短波长分量。因此，应根据不同的地质条件及野外观测系统，选择合适的平滑半径。 　　2.3计算静校正量。 　　（1）在地震资料处理过程中，每个站点的静校正量包括长波长静校正量和短波长静校正量。用反演得到的近地表V-Z模型就可计算激发点和接收点的静校正量，所得量是用于解决长波长问题的。 　　（2）其中，长波长量包括：一是从地表高程校正到浮动基准面（或浮动面），这部分量是由导出的速度计算的垂直旅行时，即高频量；二是从浮动面校正到固定面，此部分量是采用替换速度计算出的旅行时，即低频量。将计算出的校正量进行高低频分离，分离出的高频量应用于静校正中，将其校正到浮动面上。而分离出的低频量先保留着，待叠加之后，再将地震数据校正到固定面，在固定面上进行偏移处理。在实际处理中