一种初至波迭代拾取的方法技术

本发明专利技术是石油物探资料处理的初至波迭代拾取的方法，计算每一道每个采样点处的能量比值，得到能量比值曲线，确定初至预测值，以此为起点找到其中最大振幅值，计算相邻道的初至时间差的绝对值和的平均值，取一道和前后两道初至时间，排列计算可信度系数，对可信度系数为0的初至利用前后相邻4个可靠初至进行拉格朗日多项式拟合插值，得到初至预测值。本发明专利技术对于初至波形变化较大、各种波相互干扰，同时外界干扰等的影响的复杂地区，可以获得准确的初至时间，解决了初至拾取方法的难题。

【技术实现步骤摘要】

【技术保护点】
１．一种初至波迭代拾取的方法，特点是通过以下技术方案实现：１）激发并采集记录单炮地震数据；２）利用以下公式计算每一道每个采样点处的能量比值：（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）Ｒ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｒ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＝（／ｍｏ）?（ｍｆｒａｃ）?（ｍｒｏｗ）?（ｍｓｕｐ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｆｒａｃ）?（ｍｉ）Ｎ（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｏ）＋（／ｍｏ）?（ｍｓｕｐ）?（ｍｉ）Ａ（／ｍｉ）?（ｍｎ）２（／ｍｎ）?（／ｍｓｕｐ）?（ｍｉ）Ｗ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｆｒａｃ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｎ）４（／ｍｎ）?（／ｍｓｕｐ）?（ｍｏ）＆ｔｉｍｅｓ；（／ｍｏ）?（ｍｓｕｐ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｉ）Ｎ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｎ）２（／ｍｎ）?（／ｍｓｕｐ）?（／ｍｒｏｗ）?（ｍｓｕｐ）?（ｍｉ）ｒ（／ｍｉ）?（ｍｎ）２（／ｍｎ）?（／ｍｓｕｐ）?（／ｍｆｒａｃ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｎ）２（／ｍｎ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｒｏｗ）?（／ｍａｔｈ）（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｕｎｄｅｒｏｖｅｒ）?（ｍｉ）＆Ｓｉｇｍａ；（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｐ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（ｍｉ）ｒ（／ｍｉ）?（／ｍｕｎｄｅｒｏｖｅｒ）?（ｍｓｕｂｓｕｐ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｉ）ｐ（／ｍｉ）?（ｍｎ）２（／ｍｎ）?（／ｍｓｕｂｓｕｐ）?（／ｍｒｏｗ）?（／ｍａｔｈ）（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）Ｎ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｕｎｄｅｒｏｖｅｒ）?（ｍｉ）＆Ｓｉｇｍａ；（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｐ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｉ）ｒ（／ｍｉ）?（／ｍｒｏｗ）?（ｍｒｏｗ）?（ｍｉ）ｒ（／ｍｉ）?（ｍｏ）＋（／ｍｏ）?（ｍｉ）Ｗ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｕｎｄｅｒｏｖｅｒ）?（ｍｓｕｂｓｕｐ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｉ）ｐ（／ｍｉ）?（ｍｎ）２（／ｍｎ）?（／ｍｓｕｂｓｕｐ）?（／ｍｒｏｗ）?（／ｍａｔｈ）（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）Ａ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｆｒａｃ）?（ｍｎ）１（／ｍｎ）?（ｍｉ）ｎ（／ｍｉ）?（／ｍｆｒａｃ）?（ｍｕｎｄｅｒｏｖｅｒ）?（ｍｉ）＆Ｓｉｇｍａ；（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｒ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（ｍｉ）ｎ（／ｍｉ）?（／ｍｕｎｄｅｒｏｖｅｒ）?（ｍｏ）｜（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｉ）ｒ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ）｜（／ｍｏ）?（／ｍｒｏｗ）?（／ｍａｔｈ）式中：每一道的第一个采样点到当前采样点的长度为长时窗，ｒ为长时窗的结束样点，ｎ为一道的采样点数，Ｎ为当前点之后Ｗ个采样点的能量之和，Ｗ为短时窗的计算点数，Ａ为所用道的所有点振幅绝对值和的平均值，Ｍ为所用道的长时窗内所有采样点的能量之和；３）取任一道，利用步骤２）中公式（２）计算该道每个采样点处的能量比值，得到能量比值曲线，在能量比值曲线上确定门槛值对应的时间点，在对应的原始地震道上该时间点为初至预测值；４）以初至预测值为起点，长度缺省３０个采样点为第一时窗，找到其中最大振幅值，修正波峰点的能量参数与最大振幅的乘积，作为振幅门槛值，以初至预测值为起点，长度是第一时窗的１／２为第二时窗，第二时窗内第一个振幅大于振幅门槛值的位置，就是所要找的初至波峰点，如果在第二个时窗内找不到满足条件的振幅大于振幅门槛值的位置，则该道初至不存在；５）计算相邻道的初至时间差的绝对值和的平均值式中：ｔ１，ｔ２，．．．，ｔｍ是第一道到最后一道的初至时间，ｍ为一炮的总道数；以五道为一组，任一道初至时间ｔｉ，ｉ为道号，ｉ大于等于１小于等于ｍ，满足下面条件中一个，该道初至的可信度系数为１，如果均不满足，该道初至的可信度系数为０：第一个条件：ｊ＝ｉ－４，ｉ－３，ｉ－２，ｉ－１；第二个条件：ｊ＝ｉ－３，ｉ－２，ｉ－１，ｉ＋１；第三个条件：ｊ＝ｉ－２，ｉ－１，ｉ＋１，ｉ＋２；第四个条件：ｊ＝ｉ－１，ｉ＋１，ｉ＋２，ｉ＋３；第五个条件：ｊ＝ｉ＋１，ｉ＋２，ｉ＋３，ｉ＋４；ｔｊ为ｔｉ相邻道的初至时间，ｊ为道号，ｊ大于等于１小于等于ｍ，门槛值ｋ为大于０小于等于２，缺省情况下为１。６）每炮除...

【技术特征摘要】

【专利技术属性】
技术研发人员：许银坡，杨海申，蒋先艺，郭敏，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：11