[发明专利]一种确定三维工区地震数据外边界的方法及装置

说明书

本申请提供了一种确定三维工区地震数据外边界的方法及装置，该方法包括：包括根据地震数据建立第一二维数组、第二二维数组及第三二维数组，每一二维数组中的元素与地面面元一一对应；对每一二维数组中的元素赋予初始值；将地震数据中被覆盖的地面面元对应的第一二维数组的元素赋予一个与初始值不同的标记值；根据赋值后的第一二维数组，确定第一维的边界及第二维的边界，并根据第一维的边界对第二二维数组的相应元素赋予标记值，根据第二维的边界对第三二维数组的相应元素赋予标记值；若第二二维数组和第三二维数组中相同坐标的元素值均为标记值，则将标记值填入第一二维数组的对应坐标。本申请可提高地震成像程序运行效率。

技术领域

本发明涉及计算机数据处理技术领域，具体涉及一种确定三维工区地震数据外边界的方法及装置。

背景技术

三维地震勘探技术引入到石油勘探行业以来，已在石油勘探中取得了显著的地质效果，其解决复杂构造地质问题的能力被实际生产所广泛证实。随着勘探技术的不断进步，勘探需求和难度也在不断增加，三维勘探从最初的平原地带逐渐发展到戈壁沙漠、复杂山地、海洋等。最为理想的三维地震工区在大地坐标平面内常常被设计为规则的长方形，但是受到以下因素的影响，野外采集实际得到的地震数据边界范围千奇百怪。

复杂地表因素：在山前带受到山体高度、山前堆积物等因素的影响，在工区边界，地震采集设备无法布设到设计位置，致使地震采集不能完成，地震数据在工区边界形成缺口。

地表建筑物因素：受到农业设施如水库、建筑设施如村庄，海上钻井平台等因素的影响，地震采集工作往往需要劈开这些障碍，这样地震工区会形成中心空白袋。

不同年代三维连片作业因素：有些探区在不同年份采集了多块不同方向、不同大小、不同地质目标的多块三维数据，后续地震数据处理时地质家要求把多块三维进行连片处理，形成如图3所示地震工区。

综合上述几个因素，三维地震数据实际覆盖范围具有以下特点：外边界为不规则多边形，工区内可能会存在一个或多个数据空洞，空洞边界为不规则形状。

在后续的地震数据处理和地质成图过程中，各个步骤对地震数据覆盖范围的要求会有很大的差别。以时间域积分法叠前偏移成像为例，其成像数据体的覆盖范围是由参数(即成像结果最小、最大线号，最小、最大CMP号)确定的，其形状为一标准的长方形。在地震构造解释和成图过程中，只要求提供工区外边界以内成像结果数据，丢弃边界外所有数据，并且内部空洞中的数据必须保留。导致实际测得的地震数据覆盖范围与数据处理所要求的范围不能够很好的匹配，地震成像程序运行效率较低。

发明内容

针对现有技术中的问题，本申请提供一种自动获取三维工区外边界的方法及装置，能够快速、准确的确定出三维地震工区数据外边界的范围。方法成果可用于地震成像方法计算或者数据输出阶段，能够屏蔽工区外边界以外所有面元的计算和数据输出任务，提高地震成像程序运行效率。

为解决上述技术问题，本申请提供以下技术方案：

第一方面，本申请提供一种确定三维工区地震数据外边界的方法，包括：

根据所述地震数据建立第一二维数组、第二二维数组及第三二维数组，每一二维数组中的元素与地面面元一一对应；

对每一二维数组中的元素赋予初始值；

将所述地震数据中被覆盖的地面面元对应的所述第一二维数组的元素赋予一个与初始值不同的标记值；

根据赋值后的所述第一二维数组，确定第一维的边界及第二维的边界，并根据第一维的边界对所述第二二维数组的相应元素赋予标记值，根据第二维的边界对所述第三二维数组的相应元素赋予标记值；

若第二二维数组和第三二维数组中相同坐标的元素值均为所述标记值，则将所述标记值填入所述第一二维数组的对应坐标。