用于识别断层的相干增强处理方法
【专利摘要】一种用于识别断层的相干增强处理方法，通过确定目的层段在地震数据体上的双程反射时间范围；对目的层段对应的原始地震资料进行构造导向滤波，再进行相干运算，得到构造导向滤波本征值相干数据体；对数据体进行反射强度交流分量的运算，得到交流分量数据体，将数据体正负对调，得到正负值对调的相干反射强度交流分量数据体；将数据体进行数据体充零运算，将表现为断层的低相干的正值全部保留，得到正值数据体；将相干反射强度交流分量的正值数据体进行蚂蚁追踪运算，得到蚂蚁追踪数据体；将蚂蚁追踪数据体与原始地震进行融合，生成融合数据体；利用在原始地震数据体上解释的层位数据在融合数据体上提取沿层切片，得到断层与原始地震平面图。
【专利说明】
用于识别断层的相干增强处理方法
【技术领域】
[0001]本发明涉及地球物理勘探方法，是一种用于识别断层的相干增强的处理方法。
【背景技术】
[0002]相干是识别断层有效的一种地球物理方法，目前的相干有道相关相干、方差、相似性相干、本征值相干等相干方法。但无论哪种相干方法，在地震资料品质差的情况下，都不能获得理想的效果。为了突出相干效果，目前流行一些增强相干的方法，包括在相干运算之前对地震资料进行叠后处理，包括优势频带滤波、构造导向滤波等，虽然相干得到了增强，但有时效果仍不理想。目前还有的方法是在相干之后进行数据处理，包括对相干数据体进行蚂蚁追踪等，但蚂蚁追踪受制于相干效果的影响，地层因素影响较重，所以蚂蚁追踪的结果看起来过于破碎，除了断层因素外地层因素也包含在内，很难区分哪些是断层，哪些是地层。
【发明内容】
[0003]为解决上述问题，本发明目的在于提供一种用于识别断层的相干增强的处理方法，以此，有效解决上述地震资料品质较差的问题。
[0004]本发明所提供的用于识别断层的相干增强处理方法，具体包含如下步骤:步骤一:通过钻井分层、声波时差测井曲线、密度测井曲线制作合成地震记录，并进行井震层位标定，通过研究区内已知井资料建立井震关系，确定目的层段在地震数据体上的双程反射时间范围；步骤二:将目的层段对应的原始地震资料置于所述双程反射时间范围内，以所述双程反射时间范围为时窗对目的层段对应的原始地震资料进行构造导向滤波，得到导向滤波地震资料数据体一；步骤三:对导向滤波地震资料数据体一进行相干运算，得到构造导向滤波相干数据体二 ；步骤四:对构造导向滤波相干数据体二进行反射强度交流分量的运算，得到相干体反射强度交流分量数据体三，将相干体反射强度交流分量数据体三正负对调，得到正负值对调的相干反射强度交流分量数据体四；步骤五:将正负值对调的相干反射强度交流分量数据体四进行数据体赋值运算，将表现为断层的低相干的正值全部保留，得到相干反射强度交流分量的正值数据体五；步骤六:将相干反射强度交流分量的正值数据体五进行蚂蚁追踪运算，得到相干反射强度交流分量正值体的蚂蚁追踪数据体六；步骤七:将相干反射强度交流分量正值体的蚂蚁追踪数据体六与原始地震进行融合，生成融合数据体七；步骤八:利用在原始地震数据体上解释的层位数据在融合数据体七上提取沿层切片，得到断层与原始地震平面图。
[0005]在上述实施例中，所述导向滤波的算法优选的还包含:中值滤波或均值滤波或主分量滤波。
[0006]在上述实施例中，所述步骤三中的相干算法优选的还包含:本征值相干或相关或相似性相干。
[0007]在上述实施例中，优选的所述相干算法还包含:通过5点相干或7点相干或9点相干，得到构造导向滤波本征值相干数据体二。
[0008]在上述实施例中，所述步骤四中优选的还包含:数据体三在数值上是正值和负值相间的数据体，低相干的断层表现为负值，高相干的地层信息表现为正值。
[0009]在上述实施例中，所述数据体四中优选的还包含:断层表现为正值，地层信息表现为负值。
[0010]在上述实施例中，所述步骤六中的蚁群算法优选的还包含:非侵略性蚁群算法或侵略性算法，得到相干反射强度交流分量正值体的蚂蚁追踪数据体六。
[0011]在上述实施例中，所述步骤七中优选的还包含:融合时蚂蚁追踪数据体范围选择为 _0.7-1 ο
[0012]在上述实施例中，优选的还包含:所述原始地震数据体选择为8位数据体，生成融合数据体七。
[0013]本发明的有益技术效果在于:在地震资料品质较差的情况下，通过本发明能够得到更好的识别断层的数据体，与实际钻井情况吻合较好。
【专利附图】
【附图说明】
[0014]此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中:
[0015]图1为利用钻井分层、声波测井曲线、密度测井曲线、制作的合成地震记录和层位标定；
[0016]图2为构造导向滤波前后地震剖面的对比图，左图为原始地震剖面，右图为构造导向滤波剖面；
[0017]图3为构造导向滤波相干数据体的时间切片；
[0018]图4为构造导向滤波相干数据体的反射强度交流分量数据体的时间切片；
[0019]图5为将正负对调的构造导向滤波相干数据体的反射强度交流分量数据体负数赋值为O的数据体的时间切片；
[0020]图6为将图4负数赋值为O的数据体进行蚂蚁追踪后的数据体的时间切片；
[0021]图7蚂蚁追踪数据体与原始数据体融合的数据体的时间切片。
【具体实施方式】
[0022]为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
[0023]请参考图1至图7所示，本发明所提供的用于识别断层的相干增强的处理方法具体如下:
[0024]步骤一:通过钻井分层、声波时差测井曲线、密度测井曲线制作合成地震记录，并进行井震层位标定，通过研究区内已知井资料建立井震关系，确定目的层段在地震数据体上的双程反射时间范围；例如:通过钻井分层、声波时差测井曲线、密度测井曲线制作合成地震记录，并进行井震层位标定，通过研究区内已知井资料建立井震关系，确定目的层段在地震数据体上的双程反射时间范围，根据图1所示的标定结果，目的层段的双程反射时间确定为2400-2900毫秒。
[0025]步骤二:将目的层段对应的原始地震资料置于所述双程反射时间范围内，以所述双程反射时间范围为时窗对目的层段对应的原始地震资料进行构造导向滤波，得到导向滤波地震资料数据体一；例如:将目的层段对应的原始地震资料置于2400-2900ms内，以该时间范围为时窗对原始地震资料进行构造导向滤波，为了取得最佳效果，导向滤波的算法可采用主分量滤波，得到导向滤波后地震资料数据体一，图2所示为构造导向滤波后的地震数据体一的地震剖面。
[0026]在上述实施例中，所述导向滤波的算法采用中值滤波或均值滤波或主分量滤波均可得到数据体一；其中主分量滤波效果较好，但本发明并不限制该导向滤波仅仅为主分量滤波。
[0027]步骤三:对导向滤波地震资料数据体一进行相干运算，得到构造导向滤波相干数据体二 ；例如:对导向滤波地震资料数据体一进行相干运算，为了得到最佳效果，相干算法可采用本征值相干，为了获得较快的运算速度，本征值相干选择5点相干，得到构造导向滤波本征值相干数据体二，图3所示为得到的数据体二的2800ms时间切片，由于地震资料品质差的原因，此时效果比较杂乱，断层不清晰。
[0028]在上述实施例中，所述步骤三中的相干算法采用本征值相干、相关、相似性相干都可得到需要的数据体二 ；其中本征值相干算法效果较好，但本发明并不限制该相干算法仅仅为本征值相干。
[0029]该本征值相关算法采用5点相干或7点相干或9点相干，得到构造导向滤波本征值相干数据体二 ；其中5点想干算法获得的数据体二效果较好。
[0030]步骤四:对构造导向滤波相干数据体二进行反射强度交流分量的运算，得到相干体反射强度交流分量数据体三，将相干体反射强度交流分量数据体三正负对调，得到正负值对调的相干反射强度交流分量数据体四；例如:对构造导向滤波本征值相干数据体二进行反射强度交流分量的运算，得到相干体反射强度交流分量数据体三，数据体三在数值上是正值和负值相间的数据体，低相干的断层表现为负值，高相干的地层信息表现为正值；将数据体三乘以-1得到正负符号对调的相干体反射强度交流分量数据体四，数据体中低相干的断层表现为正值，高相干的地层信息表现为负值。图4为数据体四的2800ms时间切片。
[0031]步骤五:将正负值对调的相干反射强度交流分量数据体四进行数据体赋值运算，将表现为断层的低相干的正值全部保留，得到相干反射强度交流分量的正值数据体五；例如:将相干体反射强度交流分量数据体四进行数据体赋值运算，将表现为地层因素的高相干的负数值全部赋值为0，将表现为断层的低相干的正值全部保留，得到相干反射强度交流分量的正值数据体五，图5为数据体五的2800ms时间切片，高相干的因素全部被充零。
[0032]步骤六:将相干反射强度交流分量的正值数据体五进行蚂蚁追踪运算，得到相干反射强度交流分量正值体的蚂蚁追踪数据体六；例如:将相干反射强度交流分量的正值数据体五进行蚂蚁追踪运算，为了获得最佳效果，运算方法采用非侵略性蚁群算法，得到相干反射强度交流分量正值体的蚂蚁追踪数据体六，图6为数据体六的2800ms时间切片，此时断层非常清晰。
[0033]在上述实施例中，所述步骤六中的蚁群算法采用非侵略性蚁群算法和侵略性算法可得到需要的蚂蚁追踪数据体六；其中非侵略性蚁群算法效果较好，但本发明并不限制该蚁群算法仅仅为非侵略性蚁群算法。
[0034]步骤七:将相干反射强度交流分量正值体的蚂蚁追踪数据体六与原始地震进行融合，生成融合数据体七；例如:将相干反射强度交流分量正值体的蚂蚁追踪数据体六与原始地震进行融合，为了获得较好的效果，融合时蚂蚁追踪数据体范围选择为-0.7-1，为了获得较快的运算速度，原始地震数据体选择为8位数据体，生成融合数据体七。
[0035]最后，利用在原始地震数据体上解释的层位数据在融合数据体七上提取沿层切片，得到断层与原始地震平面图，图7为平面图结果，可以用于断层解释。
[0036]通过上述实施例，在地震资料品质较差的情况下，通过本发明能够得到更好的识别断层的数据体，与实际钻井情况吻合较好。
[0037]以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【权利要求】
1.一种用于识别断层的相干增强处理方法，其特征在于，所述处理方法包含: 步骤一:通过钻井分层、声波时差测井曲线、密度测井曲线制作合成地震记录，并进行井震层位标定，通过研究区内已知井资料建立井震关系，确定目的层段在地震数据体上的双程反射时间范围； 步骤二:将目的层段对应的原始地震资料置于所述双程反射时间范围内，以所述双程反射时间范围为时窗对目的层段对应的原始地震资料进行构造导向滤波，得到导向滤波地震资料数据体一； 步骤三:对导向滤波地震资料数据体一进行相干运算，得到构造导向滤波相干数据体-* ; 步骤四:对构造导向滤波相干数据体二进行反射强度交流分量的运算，得到相干体反射强度交流分量数据体三，将相干体反射强度交流分量数据体三正负对调，得到正负值对调的相干反射强度交流分量数据体四； 步骤五:将正负值对调的相干反射强度交流分量数据体四进行数据体赋值运算，将表现为地层因素的高相干的负值全部赋值为零，将表现为断层的低相干的正值全部保留，得到相干反射强度交流分量的正值数据体五； 步骤六:将相干反射强度交流分量的正值数据体五进行蚂蚁追踪运算，得到相干反射强度交流分量正值体的蚂蚁追踪数据体六； 步骤七:将相干反射强度交流分量正值体的蚂蚁追踪数据体六与原始地震进行融合，生成融合数据体七； 步骤八:利用在原始地震数据体上解释的层位数据在融合数据体七上提取沿层切片，得到断层与原始地震平面图。
2.根据权利要求1所述的用于识别断层的相干增强处理方法，其特征在于，所述导向滤波的算法包含:中值滤波或均值滤波或主分量滤波。
3.根据权利要求1所述的用于识别断层的相干增强处理方法，其特征在于，所述步骤三中的相干算法包含:本征值相干或相关或相似性相干。
4.根据权利要求3所述的用于识别断层的相干增强处理方法，其特征在于，所述相干算法包含:通过5点相干或7点相干或9点相干,得到构造导向滤波本征值相干数据体二。
5.根据权利要求1所述的用于识别断层的相干增强处理方法，其特征在于，所述步骤四中还包含:数据体三在数值上是正值和负值相间的数据体，低相干的断层表现为负值，高相干的地层信息表现为正值。
6.根据权利要求5所述的用于识别断层的相干增强处理方法，其特征在于，所述数据体四中，断层表现为正值，地层信息表现为负值。
7.根据权利要求1所述的用于识别断层的相干增强处理方法，其特征在于，所述步骤六中的蚁群算法包含:非侵略性蚁群算法或侵略性算法，得到相干反射强度交流分量正值体的蚂蚁追踪数据体六。
8.根据权利要求1所述的用于识别断层的相干增强处理方法，其特征在于，所述步骤七中还包含:融合时蚂蚁追踪数据体范围选择为-0.7-1 ο
9.根据权利要求8所述的用于识别断层的相干增强处理方法，其特征在于，所述原始地震数据体选择为8位数据体，生成融合数据体七。
【文档编号】G06F19/00GK104459801SQ201410756927
【公开日】2015年3月25日 申请日期:2014年12月10日 优先权日:2014年12月10日
【发明者】陈志刚, 马辉, 唐建超, 王亚玲, 潘良云 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司