一种从叠后地震数据采集脚印衰减的方法
【专利摘要】本发明提供一种从叠后地震数据采集脚印衰减的方法，所述方法包括：将叠后地震数据划分为多个时窗，对每个时窗执行以下步骤，对时窗内的叠后地震数据进行截断奇异值分解，得到分解后的地震数据；计算分解后的地震数据在不同秩值下的平均功率谱；对叠后地震数据进行多项式拟合，利用多项式拟合的结果合成地震纯波模型数据；计算地震纯波模型数据的平均功率谱；利用地震纯波模型数据的平均功率谱和分解后的地震数据在不同秩值下的平均功率谱，计算分解后的地震数据在不同秩值下的信噪比谱；将计算的信噪比谱之中的最大信噪比谱所对应的秩值作为地震数据重构的最佳秩值；利用最佳秩值对叠后地震数据进行重构，得到重构后的地震数据。
【专利说明】一种从叠后地震数据采集脚印衰减的方法
【技术领域】
[0001]本发明属于石油勘探地震资料处理领域，特别涉及一种从叠后地震数据采集脚印衰减的方法。
【背景技术】
[0002]随着岩性油气藏和隐蔽油气藏勘探对地震成像分辨率要求的提高，对地震资料的细节刻画要求也更高了，不仅不规则的采集和空间采样不足会产生采集脚印，规则的三维观测系统、规则的空间采样同样也会产生采集脚印，它们都影响地质目标的高精度、高质量地震成像，影响AVO属性的分析和速度分析，这种“采集脚印”又称“采集痕迹”，一般表现在三维叠加数据上。
[0003]目前针对叠后地震数据采集脚印衰减的常用方法有:道混合法、F-K滤波法、三维真振幅倾角滤波法、针对时间切片的F-KX-KY法、自适应滤波法等，虽然这些方法都能不同程度地衰减叠后地震数据中的采集脚印，但处理后的时间切片上往往会以模糊横向的振幅变化为代价。例如，F-KX-KY陷波滤波法利用采集脚印在时间切片上的X、Y方向上存在等间隔的现象，与时间方向的50Hz工业电干扰在频率域表现为一样的脉冲，采集脚印会在频率切片的KX、KY域也表现为脉冲，则可在F-KX-KY域用陷波滤波来消除采集脚印，但由于采集脚印在X、Y方向的频率切片上存在数个等间隔脉冲，需要一定的频率计算宽度和个数，由于无法准确计算出其间隔宽度，为了衰减掉采集脚印，就会给出较大的计算宽度和个数，从而造成部分有效信号被当作噪声脉冲衰减掉，使得采集脚印衰减后的时间切片横向振幅较模糊。
[0004]对叠后地震数据进行采集脚印衰减，获得高信噪比地震数据是高精度地震解释的关键，但是目前常用的采集脚印衰减方法在处理地震数据过程中无法对其信噪比进行衡量和质控，使得对地震资料衰减程度不同，其噪声压制效果不同，往往得不到最佳信噪比的地震数据。例如，常用的SVD截断奇异值分解压制采集脚印法利用分解后对角线上的元素，将其作为奇异值，并根据大的奇异值代表着有效信号，小的奇异值代表着噪声，对大的奇异值进行重构从而获得采集脚印衰减后的地震数据，但是该方法没有一个有效判断奇异值大小的根据，以及大奇异值的个数，使得不同个数的奇异值重构得到不同信噪比的地震数据，无法对处理过程中的地震数据进行有效衡量和质控。
【发明内容】
[0005]本发明的目的为利用截断奇异值分解重构衰减采集脚印技术获得具有最佳信噪比的地震数据，使重构后的地震数据采集脚印与有效信号得到最佳分离，且时间切片上振幅较清晰，较大程度上改善地震数据去噪效果。
[0006]本发明提供一种从叠后地震数据采集脚印衰减的方法，所述方法包括:将叠后地震数据划分为多个时窗，对每个时窗执行以下步骤，对时窗内的叠后地震数据进行截断奇异值分解，得到分解后的地震数据；计算分解后的地震数据在不同秩值下的平均功率谱；对叠后地震数据进行多项式拟合，利用多项式拟合的结果合成地震纯波模型数据；计算地震纯波模型数据的平均功率谱；利用地震纯波模型数据的平均功率谱和分解后的地震数据在不同秩值下的平均功率谱，计算分解后的地震数据在不同秩值下的信噪比谱；将计算的信噪比谱之中的最大信噪比谱所对应的秩值作为地震数据重构的最佳秩值；利用最佳秩值对叠后地震数据进行重构，得到重构后的地震数据。
[0007]此外，所述相邻两个时窗的移动间隔为半个时窗。
[0008]此外，计算分解后的地震数据在不同秩值下的信噪比谱的公式为，
【权利要求】
1.一种从叠后地震数据采集脚印衰减的方法，所述方法包括:将叠后地震数据划分为多个时窗，对每个时窗执行以下步骤， (1)对时窗内的叠后地震数据进行截断奇异值分解，得到分解后的地震数据； (2)计算分解后的地震数据在不同秩值下的平均功率谱； (3)对叠后地震数据进行多项式拟合，利用多项式拟合的结果合成地震纯波模型数据； (4)计算地震纯波模型数据的平均功率谱； (5)利用地震纯波模型数据的平均功率谱和分解后的地震数据在不同秩值下的平均功率谱，计算分解后的地震数据在不同秩值下的信噪比谱； (6)将计算的信噪比谱之中的最大信噪比谱所对应的秩值作为地震数据重构的最佳秩值； (7)利用最佳秩值对叠后地震数据进行重构，得到重构后的地震数据。
2.如权利要求1所述的方法，其中，所述相邻两个时窗的移动间隔为半个时窗。
3.如权利要求1所述的方法，其中，计算分解后的地震数据在不同秩值下的信噪比谱的公式为，
【文档编号】G01V1/36GK103558636SQ201310554758
【公开日】2014年2月5日 申请日期:2013年11月8日 优先权日:2013年11月8日
【发明者】张华 , 何光明, 陈爱萍, 金德刚, 罗红明, 曹中林, 尹陈, 康亮 申请人:中国石油集团川庆钻探工程有限公司地球物理勘探公司