地震相干属性计算方法及系统与流程
本发明涉及地震勘探和开发领域，更具体地，涉及一种地震相干属性计算方法及系统。
背景技术：
地震勘探是寻找油气最有效的方法，野外采集的地震数据，经过室内精细的常规处理，得到了叠后数据体。为了油藏描述的需要，就要进行地震属性的计算，而大部分地震属性提取方式都是一种规则的计算，对地下地质条件变化小的情况下是可以接受的。当地质条件变化剧烈(各项异性)的情况下，就需要一种能够适合各种情况的地震相干属性计算方法。因此，有必要开发一种地震相干属性计算方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现要素：
本发明提出了一种地震相干属性计算方法及系统，其能有效进行地震属性的计算，具有更好的断层方向性识别能力，分辨率较高，抗噪性较好，为精细的油藏描述、断裂特征刻画提供有效保障。
根据本发明的一方面，提出了一种地震相干属性计算方法。所述方法可以包括：步骤1：根据地震资料，获得工区内的采样点分布图；步骤2：针对每个采样点，重复执行子步骤201-203：子步骤201：针对所述采样点设定采样点坐标系和菱形边界图，其中所述采样点作为所述采样点坐标系的原点和所述菱形边界图的长短轴交点，所述菱形边界图的长轴与所述采样点坐标系的x轴重合，所述菱形边界图的短轴与所述采样点坐标系的y轴重合；子步骤202：以所述采样点为旋转中心，对所述菱形边界图进行旋转，获得所述菱形边界图内的有效点；子步骤203：设定时窗，根据所述菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定所述采样点的最大相干值；步骤3：根据每个采样点的最大相干值，获得所述工区的相干属性图。
优选地，获得所述菱形边界图内的有效点包括：将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点。
优选地，还包括：设置所述菱形边界图的外接圆的外接正方形为扫描范围。
优选地，获得所述菱形边界图内的有效点包括：在所述扫描范围内，将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点，并记录所述有效点的编号。
优选地，子步骤203包括：记录每个旋转角度的有效点的编号，根据每个旋转角度的有效点的编号，获取所述有效点的时窗内的数据，进而计算每个旋转角度对应的相干值。
优选地，根据以下公式(1)计算相干值：
其中，c3为相干值，c为有j道每道有n个样点组成的地震数据体，
优选地，所述菱形边界图的长短轴之比为5:2。
根据本发明的另一方面，提出了一种地震相干属性计算系统，可以包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：步骤1：根据地震资料，获得工区内的采样点分布图；步骤2：针对每个采样点，重复执行子步骤201-203：子步骤201：针对所述采样点设定采样点坐标系和菱形边界图，其中所述采样点作为所述采样点坐标系的原点和所述菱形边界图的长短轴交点，所述菱形边界图的长轴与所述采样点坐标系的x轴重合，所述菱形边界图的短轴与所述采样点坐标系的y轴重合；子步骤202：以所述采样点为旋转中心，对所述菱形边界图进行旋转，获得所述菱形边界图内的有效点；子步骤203：设定时窗，根据所述菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定所述采样点的最大相干值；步骤3：根据每个采样点的最大相干值，获得所述工区的相干属性图。
优选地，获得所述菱形边界图内的有效点包括：将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点。
优选地，根据以下公式(1)计算相干值：
其中，c3为相干值，c为有j道每道有n个样点组成的地震数据体，
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的地震相干属性计算方法的步骤的流程图。
图2示出了根据本发明的一个实施例的菱形边界图在采样点坐标系中旋转θ角度的示意图。
图3示出了根据本发明的一个实施例的二维地震测线1的叠后剖面的示意图。
图4示出了根据本发明的一个实施例的二维地震测线1常规相干剖面的示意图。
图5示出了根据本发明的一个实施例的二维地震测线1菱形相干剖面的示意图。
图6示出了根据本发明的一个实施例的三维地震1200ms切片的示意图。
图7示出了根据本发明的一个实施例的三维地震1200ms常规相干切片的示意图。
图8示出了根据本发明的一个实施例的三维地震1200ms菱形相干切片的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的地震相干属性计算方法的步骤的流程图。
在该实施例中，根据本发明的地震相干属性计算方法可以包括：步骤1：根据地震资料，获得工区内的采样点分布图；步骤2：针对每个采样点，重复执行子步骤201-203：子步骤201：针对采样点设定采样点坐标系和菱形边界图，其中采样点作为采样点坐标系的原点和菱形边界图的长短轴交点，菱形边界图的长轴与采样点坐标系的x轴重合，菱形边界图的短轴与采样点坐标系的y轴重合；子步骤202：以采样点为旋转中心，对菱形边界图进行旋转，获得菱形边界图内的有效点；子步骤203：设定时窗，根据菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定采样点的最大相干值；步骤3：根据每个采样点的最大相干值，获得工区的相干属性图。
在一个示例中，获得菱形边界图内的有效点包括：将当菱形边界图旋转角度θ时处于菱形边界图内的所有其他采样点标记为有效点。
在一个示例中，还包括：设置菱形边界图的外接圆的外接正方形为扫描范围。
在一个示例中，获得菱形边界图内的有效点包括：在扫描范围内，将当菱形边界图旋转角度θ时处于菱形边界图内的所有其他采样点标记为有效点，并记录有效点的编号。
在一个示例中，子步骤203包括：记录每个旋转角度的有效点的编号，根据每个旋转角度的有效点的编号，获取有效点的时窗内的数据，进而计算每个旋转角度对应的相干值。
在一个示例中，根据以下公式(1)计算相干值：
其中，c3为相干值，c为有j道每道有n个样点组成的地震数据体，
在一个示例中，菱形边界图的长短轴之比为5:2。
具体地，地震波场是一种矢量波场，地震相干属性计算又是油藏描述的基础，通常情况下包括整道提取，沿目的层开时窗提取，按直角、十字排列等简单的空间组合方式等。无论采用单道、多道空间组合，还是体元的方式，实际都是一种规则的计算方式，都能满足均匀地质体解释要求；但对于较为复杂的地质体，理论上都会存在局限性。为了使局限性降到最低，采用可以任意旋转、长短轴可以变化调整的菱形地震属性图进行计算，利用菱形对角线具有的方向性来判断特殊地质体的方向、范围等特征。
图2示出了根据本发明的一个实施例的菱形边界图在采样点坐标系中旋转θ角度的示意图。
根据地震资料，获得工区内的采样点分布图；针对每个采样点，重复执行子步骤201-203：子步骤201：针对采样点设定采样点坐标系和菱形边界图，如图2所示，其中采样点作为采样点坐标系的原点和菱形边界图的长短轴交点，长短轴之比为5:2菱形边界图的长轴与采样点坐标系的x轴重合，菱形边界图的短轴与采样点坐标系的y轴重合，若地震资料为二维地震剖面，则x轴表示共深度点道集方向，y轴表示采样个数，若地震资料为三维地震剖面，则x轴表示主测线方向，y轴表示横向测线方向，设置菱形边界图的外接圆的外接正方形为扫描范围，本领域技术人员可以根据具体情况设定菱形边界图的长轴与短轴长度，包括采样点分布情况、结果精度要求等；子步骤202：以采样点为旋转中心，对菱形边界图进行旋转，在扫描范围内，将当菱形边界图旋转角度θ时处于菱形边界图内的所有其他采样点标记为有效点，并记录有效点的编号；子步骤203：记录每个旋转角度的有效点的编号，设定时窗，根据每个旋转角度的有效点的编号，获取有效点的时窗内的数据，根据公式(1)计算每个旋转角度对应的相干值，进而确定采样点的最大相干值，其中，本领域技术人员可以根据具体情况设定时窗长度；步骤3：根据每个采样点的最大相干值，获得工区的相干属性图。
本方法能有效进行地震属性的计算，具有更好的断层方向性识别能力，分辨率较高，抗噪性较好，为精细的油藏描述、断裂特征刻画提供有效保障。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
根据地震资料，获得工区内的采样点分布图；针对每个采样点，重复执行子步骤201-203：子步骤201：针对采样点设定采样点坐标系和菱形边界图，其中采样点作为采样点坐标系的原点和菱形边界图的长短轴交点，长短轴之比为5:2菱形边界图的长轴与采样点坐标系的x轴重合，菱形边界图的短轴与采样点坐标系的y轴重合，若地震资料为二维地震剖面，则x轴表示共深度点道集方向，y轴表示时间方向，若地震资料为三维地震剖面，则x轴表示inline方向，y轴表示crossline方向，设置菱形边界图的外接圆的外接正方形为扫描范围；子步骤202：以采样点为旋转中心，对菱形边界图进行旋转，在扫描范围内，将当菱形边界图旋转角度θ时处于菱形边界图内的所有其他采样点标记为有效点，并记录有效点的编号；子步骤203：记录每个旋转角度的有效点的编号，设定时窗，根据每个旋转角度的有效点的编号，获取有效点的时窗内的数据，根据公式(1)计算每个旋转角度对应的相干值，进而确定采样点的最大相干值度；步骤3：根据每个采样点的最大相干值，获得工区的相干属性图。
图3示出了根据本发明的一个实施例的二维地震测线1的叠后剖面的示意图；图4示出了根据本发明的一个实施例的二维地震测线1常规相干剖面的示意图；图5示出了根据本发明的一个实施例的二维地震测线1菱形相干剖面的示意图，其分辨率与边界刻画比图4精细。
图6示出了根据本发明的一个实施例的三维地震1200ms切片的示意图；图7示出了根据本发明的一个实施例的三维地震1200ms常规相干切片的示意图；图8示出了根据本发明的一个实施例的三维地震1200ms菱形相干切片的示意图，在分辨率及边界刻画能力上有较大的提升。
综上所述，本发明能有效进行地震属性的计算，具有更好的断层方向性识别能力，分辨率较高，抗噪性较好，为精细的油藏描述、断裂特征刻画提供有效保障。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例，提供了一种地震相干属性计算系统，可以包括：存储器，存储有计算机可执行指令；处理器，处理器运行存储器中的计算机可执行指令，执行以下步骤：步骤1：根据地震资料，获得工区内的采样点分布图；步骤2：针对每个采样点，重复执行子步骤201-203：子步骤201：针对采样点设定采样点坐标系和菱形边界图，其中采样点作为采样点坐标系的原点和菱形边界图的长短轴交点，菱形边界图的长轴与采样点坐标系的x轴重合，菱形边界图的短轴与采样点坐标系的y轴重合；子步骤202：以采样点为旋转中心，对菱形边界图进行旋转，获得菱形边界图内的有效点；子步骤203：设定时窗，根据菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定采样点的最大相干值；步骤3：根据每个采样点的最大相干值，获得工区的相干属性图。
在一个示例中，获得菱形边界图内的有效点包括：将当菱形边界图旋转角度θ时处于菱形边界图内的所有其他采样点标记为有效点。
在一个示例中，还包括：设置菱形边界图的外接圆的外接正方形为扫描范围。
在一个示例中，获得菱形边界图内的有效点包括：在扫描范围内，将当菱形边界图旋转角度θ时处于菱形边界图内的所有其他采样点标记为有效点，并记录有效点的编号。
在一个示例中，子步骤203包括：记录每个旋转角度的有效点的编号，根据每个旋转角度的有效点的编号，获取有效点的时窗内的数据，进而计算每个旋转角度对应的相干值。
在一个示例中，根据以下公式(1)计算相干值：
其中，c3为相干值，c为有j道每道有n个样点组成的地震数据体，
在一个示例中，菱形边界图的长短轴之比为5:2。
本发明能有效进行地震属性的计算，具有更好的断层方向性识别能力，分辨率较高，抗噪性较好，为精细的油藏描述、断裂特征刻画提供有效保障。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
- 还没有人留言评论。精彩留言会获得点赞！