折射波静校正技术的开发与应用.doc

折射波静校正技术的开发与应用 编写人 周兴海 华北研究中心 二OO四年十一月 折射波静校正技术的开发与应用 项目负责人 周兴海 项目参加人 秋剑霞 屈志华 高现俊 单位负责人 张志国 一、 前言 4 二、概述 4 1、Sditr具有两大功能： 4 2、两种层析成像法 5 3、拾取的问题 5 4、初至叠加 5 5、几何定义问题 5 6、SDITR的优势 6 7、理论依据 6 8、SDITR的环境控制 6 9、SDITR的策略 7 三、关键词 7 四、数据的组织 7 1、输入数据 8 2、输出结果 8 五、处理过程 9 1、重新置位 9 2、计算校正量 11 3、折射波速度分析 13 4、原始道的拾取 15 六、应用实例 15 七、认识与体会 15 一、 前言 近些年来，三维折射波静校正技术发展迅速。几乎成了许多专业地球物理软件包的标志。造成这种趋势的原因之一是与小道距、短炮检距排列的高分辨采集产生的折射波静校正量有关。另外一个原因就是独立地描述同相轴近地表异常的能力。 虽然该模块儿存在多年，历经数个版本，但我们对它的了解仅限于最基本、最简单的功能，其最主要的原因就是华北地区高程起伏不大，低降速带厚度稳定，该区的地震资料静校正问题不是很突出，靠剩余静校正完全可以解决。所以无形中它也就成了制约研究开发这项技术的主要客观因素。 随着华北研究中心处理市场的对外开拓，有更多地表复杂、静校正问题严重的资料需要处理，解决静校正问题自然成了燃眉之急，这也是开发研究该项技术的最主要的原因。 二、概述 1、Sditr具有两大功能： 检查并校正几何定义错误。 计算折射波静校正量。 2、两种层析成像法 两种方法用来计算折射波静校正量。 第一种：全模型层析成像，称作GLI。（归纳线性反演） 第二种：利用旅行时分解和简单的层析成像方法。 3、拾取的问题 两种方法都从时间拾取开始，无论方法多么先进，对于低信噪比的资料来说，其性能也会大打折扣。虽然可以用局部交互拾取的方法代替自动拾取，而这个过程是非常耗时和乏味的。 4、初至叠加 多少年来，大多数折射波静校正量的质量控制是在线性动校的初至叠加的基础上完成的，运用这种方法暗含了一些简单化的假设： a、折射波速度在每个炮点----检波点矢量上被均化，它相当于一个平滑滤波。 b、初至叠加应该变得更好并且拾取时间更接近于拾取量的平均值。 5、几何定义问题 处理周期不仅取决于初至时间拾取，炮点和检波点位置的错误仍然是一个很大的麻烦。最合乎逻辑的办法就是先查出主要的异常炮点或是检波点，一旦这些异常值被剔除，以递减的数量级迭代校正其它量。 6、SDITR的优势 交互应用是能够得到的最好的完成迭代的方法，它还可以减轻管理无数个文件的负担以及等待无休止的批处理工作。 SDITR可以应用源自数据解释的全程交互，并仅限于一个主要的折射层。在完成适当的质量控制后把数据转回到原始数据上应用。此项技术非常有效，校正量或经过滤波的校正量都可以应用或舍弃，数据可以被校正也可以不被校正。同样，用户还可以一直监控每一个记录上的校正结果。 很显然，我们都希望有更大存储量的工作站而不把交互陷入长时间的等待。这里给出了一个优化的数据量。在一次工作任务当中处理50000---100000道。迄今为止，最大的一次处理数据量是350000道。7、理论依据 用代数重建法(Algebraic Reconstruction Techniques)而不是矩阵代数法(Matrix Algebra). 8、SDITR的环境控制 a、位置图 目的就是检查炮点和检波点位置。 b、几何形状选择 要处理的道的系列将由用户决定， 程序提供了三个选件。指环型、线型和变化型。 指环型 这个选件适应于高程变化不大，折射层结构简单的一般情况。 线型 它适合于地形变化、折射层结构复杂的地区。剧烈变化的地形容易引起两条接受线间产生较大的差别。对于这种类型的几何形状，扇形区被线型分量所取代。这个选件仅对重新置位有效。 变化型 这个选件只与来自SDITR的速度分析部分的空变折射波有关。 9、SDITR的策略 工作是应该从重新置位还是从校正量的计算开始？正确的做法是首先检查主要的异常值，去掉这些值以后，用一种递减的数量级迭代其他量。 三、关键词 延迟，风化层速度，替换速度，折射波速度，风化层校正量，高程校正量。（计算方法见图一） 四、数据的组织 1、输入数据 SDITR的启动始于道文件或一个来自SDS(