地震记录的处理方法技术

提供一种地震记录的处理方法，包括：选择多道地震记录并将每道地震记录变换至复赛谱域以确定每道地震记录的复赛谱序列；从所述地震记录的输入子波的复赛谱序列确定所述输入子波的振幅谱和输入子波趋势拟合曲线；对输入子波的振幅谱进行傅里叶反变换，以确定所述输入子波；计算所述多道地震记录的信噪比谱，并基于输入子波趋势拟合曲线，参考所述信噪比谱的曲线来获得期望输出子波的振幅谱；对期望输出子波的振幅谱进行傅里叶反变换以确定所述期望输出子波；通过将输入子波与期望输出子波进行匹配滤波，以获得反褶积滤波器；基于褶积模型对地震记录进行反褶积以获得反射系数序列，将反褶积滤波器与反射系数序列进行褶积以获得优化的地震记录。

【技术实现步骤摘要】

本专利技术涉及一种，更具体地讲，涉及一种在信噪比谱约束下进行复赛谱域反褶积处理地震记录的方法。
技术介绍
近年来国内油气藏勘探重点由构造油气藏转向了岩性油气藏。精细勘探目标由大的构造圈闭转向了小构造、小断层控制下的构造岩性圈闭，这对地震资料分辨率提出了更高的要求。谱模拟、谱白化等常规反褶积是实际地震资料处理中经常采用的提高地震资料分辨率的方法，这些方法是建立在地震子波的振幅谱可以表示为下面光滑的解析函数表示的基础上的。一般而言，实际地震记录的子波很难由统一的解析表达式进行刻画和描述，导致地震子波振幅谱的估算存在误差；另外，多项式L是个很难确定的用户参数；而且，对于含噪地震记录，常规反褶积方法的稳定性不够理想，特别是当地震记录的信噪比较低时，常规反褶积方法具有噪声放大效应，严重降低了谱模拟之后地震记录的信噪比。Oppenheim首先提出同态反褶积方法，该方法采用非线性滤波方法来分离地震子波和反射系数序列，对地震子波不作最小相位假设，利用平均复赛谱求取地震子波，从而得到反射系数序列。在实际应用中，由于受到噪声、反射系数序列非随机性等因素的影响，利用同态反褶积方法求取的子波稳定性较差。同态滤波的应用效果在很大程度上还取决于地震子波和反射系数序列在复赛谱的可分离性，尽管两者 在复赛谱中占据不同的区域，但要设计一个线性滤波器使二者完全分离，经过长期的实验几乎是不可能的。因此后来又发展了对数谱平均法、随机窗法等，目的之一是压制噪声的影响，另外也是为了在复赛谱中压制反射系数序列，增强地震子波的可分离性。但是，这种做法很难消除噪声和反射系数序列对地震子波相位谱的影响。由于地震资料自身的高频段能量较弱，并且受噪声影响严重，单纯依靠地震的高频信息进行反褶积处理容易造成在提高地震资料分辨率的同时降低地震资料的信噪比，不利于地震资料的下一步解释。因此，地震资料处理人员除了关心分辨率本身之外，更加关心的是地震资料在提高分辨率处理之后能否较好地保持地震资料的信噪比，但常规的反褶积方法往往都无法有效质控提高分辨后地震资料的信噪比。
技术实现思路
本专利技术的一方面在于通过将地震记录变换至复赛谱域来使地震记录中的地震输入子波与反射系数序列分离，从而消除噪声和反射系数序列对地震输入子波的影响。本专利技术的另一目的在于通过计算地震记录的信噪比谱，在信噪比谱约束下对地震记录中的输入子波振幅谱进行调整，从而在较好地保持地震数据的信噪比的情况下提高地震记录的分辨率。本专利技术提供一种，所述处理方法包括以下步骤:(1)选择多道地震记录并将每道地震记录变换至复赛谱域以确定每道地震记录的复赛谱序列，其中，Xi (t)表示所述多道地震记录中的第i道地震记录，A⑴表示地震记录Xi (t)的复赛谱序列；(2)从地震记录Xi (t)的复赛谱序列足(0中提取地震记录Xi (t)的输入子波Wi (t)的复赛谱序列4(0，并从所述复赛谱序列示(0确定所述输入子波Wi (t)的振幅谱Wi (f)和输入子波趋势拟合曲线；(3)对输入子波Wi(t)的振幅谱^&)进行傅里叶反变换，以确定所述输入子波Wi (t) ; (4)计算所述多道地震记录的信噪比谱，并基于输入子波趋势拟合曲线，参考所述信噪比谱的曲线来获得期望输出子波Wi' (t)的振幅谱w/ (f) ;(5)对期望输出子波Wi' (t)的振幅谱Wi' (f)进行傅里叶反变换以确定所述期望输出子波w/ (t);(6)通过将步骤(3)中确定的输入子波Wi (t)与步骤(5)中确定的期望输出子波w/ (t)进行匹配滤波，以获得反裙积滤波器djt) ; (7)基于裙积模型Xi(^t)=Wi (thrjt)对地震记录Xi (t)进行反褶积以获得反射系数序列A (t)，将反褶积滤波器djt)与反射系数序列& (t)进行裙积以获得优化的地震记录⑴，其中，Yi (t) =(I^t) ^ri (t), *为裙积算子。可选地，确定输入子波Wi (t)的振幅谱&江)的步骤可包括:使用预先设置的低通滤波器对地震记录Xi (t)的复赛谱序列A⑴进行低通滤波，以提取输入子波Wi (t)的复赛谱序列A (0;将输入子波Wi (t)的复赛谱序列记(0变换至频域来获得输入子波Wi (t)的振幅谱 Wi(f)。可选地，确定输入子波趋势拟合曲线的步骤可包括:计算输入子波Wi(t)的振幅谱Wi (f)在不同频率段上的均值；利用样条拟合的方法将计算出的多个均值拟合成输入子波趋势拟合曲线。可选地，计算信噪比谱的步骤可包括:通过公式本文档来自技高网...

【技术保护点】
一种地震记录的处理方法，所述处理方法包括以下步骤：（1）选择多道地震记录并将每道地震记录变换至复赛谱域以确定每道地震记录的复赛谱序列，其中，xi(t)表示所述多道地震记录中的第i道地震记录，表示地震记录xi(t)的复赛谱序列；（2）从地震记录xi(t)的复赛谱序列中提取地震记录xi(t)的输入子波wi(t)的复赛谱序列并从所述复赛谱序列确定所述输入子波wi(t)的振幅谱wi(f)和输入子波趋势拟合曲线；（3）对输入子波wi(t)的振幅谱wi(f)进行傅里叶反变换，以确定所述输入子波wi(t)；（4）计算所述多道地震记录的信噪比谱，并基于输入子波趋势拟合曲线，参考所述信噪比谱的曲线来获得期望输出子波wi′(t)的振幅谱wi′(f)；（5）对期望输出子波wi′(t)的振幅谱wi′(f)进行傅里叶反变换以确定所述期望输出子波wi′(t)；（6）通过将步骤（3）中确定的输入子波wi(t)与步骤（5）中确定的期望输出子波wi′(t)进行匹配滤波，以获得反褶积滤波器di(t)；（7）基于褶积模型xi(t)=wi(t)*ri(t)对地震记录xi(t)进行反褶积以获得反射系数序列ri(t)，将反褶积滤波器di(t)与反射系数序列ri(t)进行褶积以获得优化的地震记录yi(t)，其中，yi(t)=di(t)*ri(t)，*为褶积算子。FDA0000418395820000011.jpg,FDA0000418395820000012.jpg,FDA0000418395820000013.jpg,FDA0000418395820000014.jpg...

【技术特征摘要】
1.一种地震记录的处理方法，所述处理方法包括以下步骤: (1)选择多道地震记录并将每道地震记录变换至复赛谱域以确定每道地震记录的复赛谱序列，其中，Xi (t)表示所述多道地震记录中的第i道地震记录，无⑴表示地震记录Xi (t)的复赛谱序列； (2)从地震记录Xi(t)的复赛谱序列天的中提取地震记录Xi (t)的输入子波Wi (t)的复赛谱序列，并从所述复赛谱序列》⑴确定所述输入子波Wi⑴的振幅谱Wi (f)和输入子波趋势拟合曲线； (3)对输入子波Wi(t)的振幅谱&&)进行傅里叶反变换，以确定所述输入子波Wi (t)； (4)计算所述多道地震记录的信噪比谱，并基于输入子波趋势拟合曲线，参考所述信噪比谱的曲线来获得期望输出子波Wi' (t)的振幅谱Wi' (f)； (5)对期望输出子波w/(t)的振幅谱w/ (f)进行傅里叶反变换以确定所述期望输出子波Wi' (t); (6)通过将步骤(3)中确定的输入子波Wi(t)与步骤(5)中确定的期望输出子波Wi' (t)进行匹配滤波，以获得反褶积滤波器Cli (t)； (7)基于褶积模型Xi(t)=Wi(t)*ri(t)对地震记录Xi(t)进行反褶积以获得反射系数序列rjt)，将反褶积滤波器ddt)与反射系数序列ri(t)进行褶积以获得优化的地震记录Yi (t)，其中，Yi (t) =Cli (t)*ri(t)，* 为裙积算子。2.根据权利要求1所述的处理方法，其...

【专利技术属性】
技术研发人员：张华，何光明，陈爱萍，金德刚，罗红明，曹中林，尹陈，康亮，
申请(专利权)人：中国石油集团川庆钻探工程有限公司地球物理勘探公司，
类型：发明
国别省市：