一种确定地震子波的方法
【专利摘要】本发明是石油物探地震资料处理中确定地震子波的方法，在原有的目标函数中增加了确定地震子波能量最小的约束项，考虑了地震道中的噪音对确定地震子波精度的影响，通过寻找地震道残差能量与地震子波能量的最佳“平衡点”，确定出反射系数褶积矩阵奇异值的最佳截断阶数，通过对反射系数褶积矩阵奇异值截断来减弱地震道的噪音影响，确定地震子波。本发明不影响到较大奇异值对子波确定的贡献，与常用方法相比更有利于减弱地震道中噪音对子波确定精度的影响，有利于地震子波确定精度的提高。
【专利说明】一种确定地震子波的方法
【技术领域】
[0001]本发明涉及石油勘探开发技术，具体是地震资料处理、解释中的一种确定地震子波的方法。
【背景技术】
[0002]地震子波的确定是地震勘探中一个非常关键的问题。在正演问题中，需要通过波动方程或褶积模型结合地震子波来形成正演模拟地震数据。在反演和反褶积问题中，也需要通过地震道提取一个地震子波，不同的地震子波得到不同的反演结果。
[0003]地震子波提取的基础是褶积模型，它包括地震子波、反射系数、含有噪声的地震道数据。地震子波提取方法是利用测井资料计算出反射系数，然后结合地震道由褶积模型求出地震子波。目前常用的井震联合子波确定方法是基于合成道与地震道残差能量(以下简称地震道残差能量)最小原理:
【权利要求】
1.一种确定地震子波的方法，特点是通过以下步骤实现: 1)采集测井数据和地震道数据以及先验地震子波； 2)利用先验地震子波的时间域反射系数制作合成记录，在时窗范围内对时间域反射系数进行重新标定； 3)给定目标地震子波的长度，确定地震道数据时窗和反射系数时窗范围； 4)构建反射系数褶积矩阵； 5)构建地震道列向量X; 6)按照以下公式确定目标地震子波目标函数；
2.根据权利要求1的方法，特点是步骤I)所述的测井数据主要包括:井名、该井所在位置的大地坐标、经过合成记录准确标定后的时间域反射系数。
3.根据权利要求1的方法，特点是步骤I)所述的地震道数据，是地震采集的一道或井点附近的多道地震数据。
4.根据权利要求1的方法，特点是步骤I)所述的地震道数据，是指要确定地震子波目标层段时窗内的地震数据；地震道数据的时窗顶和底要偏离地震道的波峰和波谷，选择弱反射位置。
5.根据权利要求1的方法，特点是步骤I)所述的先验地震子波是先前用于合成记录标定所已知的地震子波或理论子波。
6.根据权利要求1的方法，特点是步骤2)中的重新标定是通过对合成道和各个井旁道数据相关分析得到反射系数和各个井旁道最佳匹配位置。
7.根据权利要求1的方法，特点是步骤2)中的时窗范围是指在给定地震时窗范围内正负30ms范围内移动时间域反射系数。
8.根据权利要求1的方法，特点是 步骤3)所述的目标地震子波长度根据实际地震资料确定，浅层地震子波长度短，深层地震子波长度长；当时间域反射系数不能满足步骤3)时窗范围的需要，减小子波长度，子波长度的最小长度，必须包括子波的主瓣和旁瓣的完整性。
9.根据权利要求1的方法，特点是步骤3)所述的目标地震子波长度为60--120ms之间。
10.根据权利要求1的方法，特点是步骤3)所述的地震道时窗要大于地震子波长度的I~1.5倍，反射系数时窗比地震道时窗范围上下各扩展半个地震子波长度。
11.根据权利要求1的方法，特点是步骤7)中不同k阶截奇异值取值范围是O< k < n，η为地震子波样点个数。
12.根据权利要求1的方法，特点是步骤7)中地震子波W(k)能量公式是根据目标地震子波目标函数公式得到。
【文档编号】G01V1/28GK103576192SQ201210273225
【公开日】2014年2月12日 申请日期:2012年8月2日 优先权日:2012年8月2日
【发明者】徐广民 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司