地震相干属性计算方法及系统技术方案

公开了一种地震相干属性计算方法及系统。该方法可以包括：步骤1：根据地震资料，获得工区内的采样点分布图；步骤2：针对每个采样点，重复执行子步骤201‑203：子步骤201：针对采样点设定采样点坐标系和菱形边界图；子步骤202：以采样点为旋转中心，对菱形边界图进行旋转，获得菱形边界图内的有效点；子步骤203：设定时窗，根据菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定采样点的最大相干值；步骤3：根据每个采样点的最大相干值，获得工区的相干属性图。本发明专利技术能有效进行地震属性的计算，具有更好的断层方向性识别能力，分辨率较高，抗噪性较好，为精细的油藏描述、断裂特征刻画提供有效保障。

Method and System for Calculating Coherent Attributes of Earthquakes

A method and system for calculating seismic coherence attributes are disclosed. The method can include: step 1: obtain the distribution map of sampling points in the work area according to seismic data; step 2: repeat the sub-step 201 203: sub-step 201: set the coordinate system of sampling points and diamond boundary map for sampling points; step 202: rotate the diamond boundary map with sampling points as the rotation center to obtain the validity of the diamond boundary map. Substep 203: Set the time window, calculate the coherence value corresponding to each rotation angle according to the data in the time window of the valid points in the diamond boundary map, and then determine the maximum coherence value of the sampling point; Step 3: According to the maximum coherence value of each sampling point, obtain the coherence attribute map of the work area. The invention can effectively calculate seismic attributes, has better fault orientation recognition ability, higher resolution and better noise resistance, and provides effective guarantee for fine reservoir description and fracture characterization.

【技术实现步骤摘要】
地震相干属性计算方法及系统
本专利技术涉及地震勘探和开发领域，更具体地，涉及一种地震相干属性计算方法及系统。
技术介绍
地震勘探是寻找油气最有效的方法，野外采集的地震数据，经过室内精细的常规处理，得到了叠后数据体。为了油藏描述的需要，就要进行地震属性的计算，而大部分地震属性提取方式都是一种规则的计算，对地下地质条件变化小的情况下是可以接受的。当地质条件变化剧烈(各项异性)的情况下，就需要一种能够适合各种情况的地震相干属性计算方法。因此，有必要开发一种地震相干属性计算方法及系统。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种地震相干属性计算方法及系统，其能有效进行地震属性的计算，具有更好的断层方向性识别能力，分辨率较高，抗噪性较好，为精细的油藏描述、断裂特征刻画提供有效保障。根据本专利技术的一方面，提出了一种地震相干属性计算方法。所述方法可以包括：步骤1：根据地震资料，获得工区内的采样点分布图；步骤2：针对每个采样点，重复执行子步骤201-203：子步骤201：针对所述采样点设定采样点坐标系和菱形边界图，其中所述采样点作为所述采样点坐标系的原点和所述菱形边界图的长短轴交点，所述菱形边界图的长轴与所述采样点坐标系的X轴重合，所述菱形边界图的短轴与所述采样点坐标系的Y轴重合；子步骤202：以所述采样点为旋转中心，对所述菱形边界图进行旋转，获得所述菱形边界图内的有效点；子步骤203：设定时窗，根据所述菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定所述采样点的最大相干值；步骤3：根据每个采样点的最大相干值，获得所述工区的相干属性图。优选地，获得所述菱形边界图内的有效点包括：将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点。优选地，还包括：设置所述菱形边界图的外接圆的外接正方形为扫描范围。优选地，获得所述菱形边界图内的有效点包括：在所述扫描范围内，将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点，并记录所述有效点的编号。优选地，子步骤203包括：记录每个旋转角度的有效点的编号，根据每个旋转角度的有效点的编号，获取所述有效点的时窗内的数据，进而计算每个旋转角度对应的相干值。优选地，根据以下公式(1)计算相干值：其中，C3为相干值，C为有J道每道有N个样点组成的地震数据体，DN*J为协方差矩阵，其中DN*J为参与计算的由J道N个样点构成的地震子体，是DN*J的转置，Tr(C)为矩阵的迹，代表协方差矩阵的能量，λ1是最大特征值，代表占优的能量，cjj是矩阵对角线上的元素，λj是矩阵非负特征值，J为参与计算的地震道数，j为参与计算的地震道编号。优选地，所述菱形边界图的长短轴之比为5:2。根据本专利技术的另一方面，提出了一种地震相干属性计算系统，可以包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：步骤1：根据地震资料，获得工区内的采样点分布图；步骤2：针对每个采样点，重复执行子步骤201-203：子步骤201：针对所述采样点设定采样点坐标系和菱形边界图，其中所述采样点作为所述采样点坐标系的原点和所述菱形边界图的长短轴交点，所述菱形边界图的长轴与所述采样点坐标系的X轴重合，所述菱形边界图的短轴与所述采样点坐标系的Y轴重合；子步骤202：以所述采样点为旋转中心，对所述菱形边界图进行旋转，获得所述菱形边界图内的有效点；子步骤203：设定时窗，根据所述菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定所述采样点的最大相干值；步骤3：根据每个采样点的最大相干值，获得所述工区的相干属性图。优选地，获得所述菱形边界图内的有效点包括：将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点。优选地，根据以下公式(1)计算相干值：其中，C3为相干值，C为有J道每道有N个样点组成的地震数据体，DN*J为协方差矩阵，其中DN*J为参与计算的由J道N个样点构成的地震子体，是DN*J的转置，Tr(C)为矩阵的迹，代表协方差矩阵的能量，λ1是最大特征值，代表占优的能量，cjj是矩阵对角线上的元素，λj是矩阵非负特征值，J为参与计算的地震道数，j为参与计算的地震道编号。本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。图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.一种地震相干属性计算方法，包括：步骤1：根据地震资料，获得工区内的采样点分布图；步骤2：针对每个采样点，重复执行子步骤201‑203：子步骤201：针对所述采样点设定采样点坐标系和菱形边界图，其中所述采样点作为所述采样点坐标系的原点和所述菱形边界图的长短轴交点，所述菱形边界图的长轴与所述采样点坐标系的X轴重合，所述菱形边界图的短轴与所述采样点坐标系的Y轴重合；子步骤202：以所述采样点为旋转中心，对所述菱形边界图进行旋转，获得所述菱形边界图内的有效点；子步骤203：设定时窗，根据所述菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定所述采样点的最大相干值；步骤3：根据每个采样点的最大相干值，获得所述工区的相干属性图。

【技术特征摘要】
1.一种地震相干属性计算方法，包括：步骤1：根据地震资料，获得工区内的采样点分布图；步骤2：针对每个采样点，重复执行子步骤201-203：子步骤201：针对所述采样点设定采样点坐标系和菱形边界图，其中所述采样点作为所述采样点坐标系的原点和所述菱形边界图的长短轴交点，所述菱形边界图的长轴与所述采样点坐标系的X轴重合，所述菱形边界图的短轴与所述采样点坐标系的Y轴重合；子步骤202：以所述采样点为旋转中心，对所述菱形边界图进行旋转，获得所述菱形边界图内的有效点；子步骤203：设定时窗，根据所述菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定所述采样点的最大相干值；步骤3：根据每个采样点的最大相干值，获得所述工区的相干属性图。2.根据权利要求1所述的地震相干属性计算方法，其中，获得所述菱形边界图内的有效点包括：将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点。3.根据权利要求2所述的地震相干属性计算方法，还包括：设置所述菱形边界图的外接圆的外接正方形为扫描范围。4.根据权利要求3所述的地震相干属性计算方法，其中，获得所述菱形边界图内的有效点包括：在所述扫描范围内，将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点，并记录所述有效点的编号。5.根据权利要求4所述的地震相干属性计算方法，其中，子步骤203包括：记录每个旋转角度的有效点的编号，根据每个旋转角度的有效点的编号，获取所述有效点的时窗内的数据，进而计算每个旋转角度对应的相干值。6.根据权利要求1所述的地震相干属性计算方法，其中，根据以下公式(1)计算相干值：其中，C3为相干值，C为有J道每道有N个样点组成的地震数据体，DN*J为协方差矩阵，其中DN*J为参与计算的由J道N个样点构成的地震子体，是DN*J的转置，Tr(C)为矩阵的迹，代表...

【专利技术属性】
技术研发人员：吴永栓，唐金良，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11