一种叠后数据的振幅能量异常现象的检测方法
【专利摘要】本发明提供了一种叠后数据的振幅能量异常现象的检测方法，属于石油勘探资料处理【技术领域】。所述方法包括：将叠后数据按照预定的时窗进行划分；在每个时窗下获得每个地震道的振幅能量值，并在每个时窗下获得全工区地震道的振幅能量值，获取每个地震道的振幅能量值相对于全工区地震道的振幅能量值的能量偏差值，将能量偏差值作为每个地震道的能量属性；将每个地震道的能量属性与预定门限值的上下限进行比较，则确定是否发生及发生哪种振幅能量异常。本发明能够定性、定量地对振幅能量异常现象进行客观评价，并且能够较好地识别出叠后数据体中的整体及局部能量异常现象，具有较高的检测准确率以及较快的检测效率。
【专利说明】-种叠后数据的振幅能量异常现象的检测方法
【技术领域】
[0001] 本发明设及一种叠后数据的振幅能量异常现象的检测方法，属于石油勘探资料处 理【技术领域】。
【背景技术】
[0002] 在地震资料处理过程中，在叠后纯波数据和成果数据剖面上经常会出现一些振 幅能量变化不均匀的现象，结合图6所示，主要表现为W下几个方面；弱能量分界线、强能 量分界线、局部的强能量、单点的野值、空间方向能量突变等等，该些异常现象的主要来自 于：
[0003] 1、采集原始记录中存在的振幅异常
[0004] 2、老资料由于存放的磁带老化造成的异常振幅
[0005] 3、叠加剖面不能完全覆盖特定部位存在的振幅异常
[0006] 4、处理参数运用不当所造成的振幅异常
[0007] 由于目前常用的处理系统中（如；Omega、CGG、GeoEast等）均无法针对叠后数据 中上述能量异常现象进行快速检测，所W长期W来主要由人工逐炮逐道翻阅剖面，依靠肉 眼来对叠加数据体进行分析，W识别上述能量异常现象，该种方法主要有W下弊端：
[000引 1、无法进行定量评价
[0009] 2、很大程度依靠个人经验，主观性较强
[0010] 3、容易造成漏检
[0011] 4、随着高密度采集技术的发展，海量数据的呈现，仅依靠人工进行识别，
[0012] 使处理质量的质控工作强度睹增，压力巨大
[0013] 5、不能满足目前资料处理工程化质量控制的需求
【发明内容】
[0014] 本发明针对【背景技术】中存在的问题，提出了一种叠后数据的振幅能量异常现象的 检测方法，从而能够W较高的检测准确率W及较快的检测效率对振幅能量异常现象进行检 测。
[0015] 本发明提供的技术方案包括：
[0016] 一种叠后数据的振幅能量异常现象的检测方法，包括：
[0017] 将所述叠后数据中的每个地震道按照预定的时窗大小、重叠百分比或时窗间增量 划分时窗；
[0018] 在每个所述时窗下获得每个地震道的振幅能量值，并在每个所述时窗下获得全工 区地震道的振幅能量值，获取所述每个地震道的振幅能量值相对于所述全工区地震道的振 幅能量值的能量偏差值，将所述能量偏差值作为每个地震道的能量属性；
[0019] 将所述每个地震道的能量属性与预定口限值的上下限进行比较，当所述每个地震 道的能量属性大于所述预定口槛值的上限时，则确定发生强振幅能量异常，当所述每个地 震道的能量属性小于所述预定口槛值的下限时，则确定发生弱振幅能量异常。
[0020] 在本发明所述的叠后数据的振幅能量异常现象的检测方法中，所述每个地震道的 振幅能量值通过最大绝对振幅计算方法、均方根振幅计算方法和/或平均绝对振幅计算方 法获得。
[0021] 在本发明所述的叠后数据的振幅能量异常现象的检测方法中，所述最大绝对振幅 计算方法通过W下公式计算获得所述每个地震道的振幅能量值：
[0022]
【权利要求】
1. 一种叠后数据的振幅能量异常现象的检测方法，其特征在于，包括： 将所述叠后数据中的每个地震道按照预定的时窗大小、重叠百分比或时窗间增量划分 时窗； 在每个所述时窗下获得每个地震道的振幅能量值，并在每个所述时窗下获得全工区地 震道的振幅能量值，获取所述每个地震道的振幅能量值相对于所述全工区地震道的振幅能 量值的能量偏差值，将所述能量偏差值作为每个地震道的能量属性； 将所述每个地震道的能量属性与预定门限值的上下限进行比较，当所述每个地震道的 能量属性大于所述预定门槛值的上限时，则确定发生强振幅能量异常，当所述每个地震道 的能量属性小于所述预定门槛值的下限时，则确定发生弱振幅能量异常。
2. 如权利要求1所述的方法，其特征在于，所述每个地震道的振幅能量值通过最大绝 对振幅计算方法、均方根振幅计算方法和/或平均绝对振幅计算方法获得。
3. 如权利要求2所述的方法，其特征在于，包括：所述最大绝对振幅计算方法通过以下 公式计算获得所述每个地震道的振幅能量值：
其中，i表示在每个预划分的时窗下每个地震道的振幅能量值，ie(1，N)， jG(1，M)，N表不地震道数，M表不时窗个数，i表不地震道序号，j表不时窗序号，t表不时 窗的起始时间，n表示时窗内采样点个数，a(ik)表示第i个地震道的第k个样点的振幅。
4. 如权利要求2所述的方法，其特征在于，所述均方根振幅计算方法通过以下公式计 算获得所述每个地震道的振幅能量值：
其中，i表示在每个预划分的时窗下每个地震道的振幅能量值，ie(1，N)， jG(1，M)，N表不地震道数，M表不时窗个数，i表不地震道序号，j表不时窗序号，t表不时 窗的起始时间，n表示时窗内采样点个数，a(ik)表示第i个地震道的第k个样点的振幅。
5. 如权利要求2所述的方法，其特征在于，所述平均绝对振幅计算方法通过以下公式 计算获得所述每个地震道的振幅能量值：
其中，i表示在每个预划分的时窗下每个地震道的振幅能量值，ie(1，N)， jG(1，M)，N表不地震道数，M表不时窗个数，i表不地震道序号，j表不时窗序号，t表不时 窗的起始时间，n表示时窗内采样点个数，a(ik)表示第i个地震道的第k个样点的振幅。
6. 如权利要求1至5任意一项所述的方法，其特征在于，所述方法还包括： 在每个所述时窗下获得每条主测线和每条联络线的振幅能量值，将所述振幅能量值与 所述全工区地震道的振幅能量值的偏差值作为每条主测线和每条联络线的能量属性； 将所述每条主测线和每条联络线的能量属性通过二维连线的方式绘制主测线和联络 线的振幅属性图，根据所述振幅属性图检测线间能量差异。
7. 如权利要求6所述的方法，其特征在于，确定所述叠后数据是否发生振幅能量异常 现象包括：根据交互可视化图件对所述叠后数据的振幅能量属性进行识别，以确定是否发 生振幅能量异常现象。
8. 如权利要求7所述的方法，其特征在于，所述叠后数据的振幅能量属性通过以下公 式获得：
其中，Su表示每个所述地震道在每个所述时窗的振幅能量属性值，$表示在每个所 述时窗下每个地震道的振幅能量值，^表示在每个所述时窗下所有地震道的振幅能量值， iG(1，N)，jG(1，M)，N表示地震道数，M表示时窗个数，i表示地震道序号，j表示时窗序 号。
9. 如权利要求8所述的方法，其特征在于，所述根据交互可视化图件对所述叠后数据 的振幅能量属性进行识别包括： 通过正态分布图将所述振幅能量属性值划分为若干区间，根据所述每个地震道在每个 时窗的振幅能量属性值的分布对所述叠后数据的振幅能量属性进行识别； 通过二维连线图显示的每条主测线的振幅属性及每条联络线的振幅属性对所述叠后 数据的振幅能量属性进行识别； 通过振幅平面属性图对所述叠后数据的振幅能量属性进行识别。
【文档编号】G01V1/28GK104502963SQ201410785466
【公开日】2015年4月8日 申请日期:2014年12月17日 优先权日:2014年12月17日
【发明者】恽春华, 王克斌, 梁桂美, 周兴元, 李道善, 刘志山 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司