[发明专利]一种采用地震属性识别断层的方法

摘要

本发明是采用地震属性识别断层的方法，根据地震数据的信噪比确定标准化高阶统计量计算中的相邻地震道数目和地震道横向组合模式，在纵向时窗内确定统计时窗，将地层倾角最大处相邻地震道的时间延迟量作为扫描时窗长度，根据地震道数目和组合模式、统计时窗和扫描时窗长度计算标准化高阶统计量，根据地震道横向组合道数和模式将统计量极大值作为提取属性，进行时间切片及剖面成图识别断层空间展布。本发明对地震资料的随机噪音有较强的压制作用，特别对地层倾角引起的不连续性干扰明显得到压制，属性抗噪性能强，对断裂系统引起的地层不连续性特征识别效果好。

主权项

一种采用地震属性识别断层的方法，特点是通过以下步骤实现：1）采集地震数据，经常规处理后得到叠后地震数据体；2）根据地震数据的信噪比确定标准化高阶计算量中的相邻地震道数目，若地震数据信噪比低则相邻地震道数目采用5道或9道，若地震数据信噪比高则相邻地震道数目采用3道；3）根据相邻地震道数目确定地震道横向组合模式；所述的确定地震道横向组合模式是：当采用3道计算时，地震道横向组合模式为：目标点与目标点所在的纵/横测线上2个相邻点；当采用5道计算时，地震道横向组合模式为：目标点与目标点所在的纵/横测线上4个相邻点；当采用9道计算时，地震道横向组合模式为：目标点与目标点周围的8个相邻点；4）在纵向时窗内确定统计时窗，所述的统计时窗为一个完整地震波形长度的1.5到3倍；5）将地层倾角最大处相邻地震道的时间延迟量作为扫描时窗长度；6）根据地震道数目和组合模式，根据统计时窗长度和扫描时窗长度，采用下式计算标准化高级统计量： C HOSN ( t , x , y , τ x , τ y ) = Σ k = - w k = w u 0 2 u x u y Σ k = - w k = + w [ u 0 u x ] 2 × Σ k = - w k = + w [ u 0 u y ] 2 3 × Σ k = - w k = + w [ u 0 ] 4 × Σ k = - w k = + w [ u x u y ] 2 6 - - - ( 6 ) u0=u(t+k,x,y)ux＝u(t+k+τx,x+dx,y)uy＝u(t+k+τy,x,y+dy) C ( t , x , y ) = max τ x , τ y ( C HOSN ( t , x , y , τ x , τ y ) ) 上式中，t表示目标点时间，x表示目标点横坐标，y表示目标点纵坐标；k为统计时窗变量，τx为x方向相邻道对目标道的时间延迟量，τy为y方向相邻道对目标道的时间延迟量，[‑w,w]为统计时窗长度；u0表示目标地震道，ux表示x方向相邻地震道，uy表示y方向相邻地震道；CHOSN（t,x,y,τx,τy）为对应于每个时间延迟量的统计量值，C（t,x,y）为目标点的标准化高阶统计量属性值；7）根据标准化高阶统计量计算中地震道的横向组合道数及组合模式得到的统计量极大值来提取属性值：3道组合得到一个统计量极大值C(t,x,y)，该值作为最终输出的属性值；5道组合得到两个统计量极大值C1（t,x,y）和C2（t,x,y），将这个两个极大值的算数平均值(C1（t,x,y）+C2（t,x,y）)/2作为最终输出的属性值；9道组合得到4个统计量极大值C1（t,x,y），C2（t,x,y）…C4（t,x,y）,将这四个极大值的算数平均值(C1（t,x,y）+C2（t,x,y）+…C4（t,x,y）)/4作为最终输出的属性值；8）对属性体进行横向时间切片及纵向剖面成图，识别断层在横向及纵向的空间展布。