专利名称：一种高效地表一致性反褶积的方法
技术领域：
本发明涉及地球物理勘探技术，属于地震资料处理过程中提高分辨率的技术范畴，是用于提高地震资料的分辩率和消除地震记录炮点、检波点和CMP之间的差异的一种高效地表一致性反褶积的方法。
背景技术：
在地球物理勘探技术中对记录的数据要进行处理，地表一致性约束下的 反褶积技术是地震资料处理的一种常规处理方法。在地表起伏、表层结构复 杂地区，当震源和接受点位置变化时，激发条件、接收条件等都发生变化， 对地震记录将造成不同影响，地表一致性反褶积可以消除这种影响，反褶积 过程中道间振幅补偿了由震源和接收器耦合引起的变化，同时能有效的提高 地震资料的分辩率。通常的地表一致性反褶积原理是假定地震记朵可以用以 下的公式表示(源自yilmaz， 1987)
xij(t)=sj(t)*hm(t)*en(t)*qi(t)+n(t)(1)
式中xij(t)是地震记录,sj(t)是炮点位置为y的响应，q(t)是检波点位置为i的响应,hm(t)是依赖于炮检距的响应,en(t)与面元有关的响应，n(t))为噪 音响应，这里假定n(t)=0。
方程(1)在频率域则为
xij(ω)= Sj(ω)Hm(ω)En(ω)Qi(ω)( 2 )
如果预测反褶积，只需估算功率谱，对功率谱取对数(为书写方便，下面用xij(ω),Sj(ω), Hm(ω),En(ω)和Qi(ω)表示对应的功率谱)
lnxij(ω)= lnSj(ω)+lnHm(ω)+lnEn(ω)+lnQi(ω) (3 )
应用最小平方方法<formula>complex formula see original document page 5</formula>
式中，lnXij(w)是地震记录的对数功率谱。 对L求偏导
<formula>complex formula see original document page 5</formula>
其中j=l,2，…J， J为炮点总个数；/ = l,2，...A/， M为以米为单位的不同炮 检距的个数；w = l，2,...7V， W为面元个数；口1，2，…/， /为检波点个数， w二lAw，2Aw，…『Aw，『为频率采样点个数。对每个频率采样点都要用高斯一 塞德尔迭代法求解。而后对每一项响应取指数，反傅立叶变换，得到相应的 自相关函数，解托布里兹方程，得到各个分量的反褶积算子。
由此可见，地表一致性反褶积需要很大的设备内存和计算量。随着大面 积三维采集的广泛应用和高密度采集的推广，需要处理的数据量越来越大， 这种比较成熟的方法无法适应大数据量处理。解决办法有采用对输入数据直 接做自相关，然后用三角窗截断自相关函数，经较少点数的傅立叶变换，对 对数功率谱迭代分解，这种方法虽然可节省内存和减少计算量，但反褶积结 果分辩率降低。也有的直接截断功率谱来节省内存和计算量，这种方法容易 造成高频异常。
发明目的
本发明提供一种适合大面积三维采集和高密度采集、既节省计算量又节省 内存空间、效率高效果好的一种高效地表一致性反褶积的方法。
发明内容
本发明提供如下解决方案
先采集叠前炮集或CMP地震记录，用常规方法截取地震记录， 一般选择 地震记录信噪比较高的时窗段，用常规的方法对地震记录段做傅立叶变换，把地震记录从时间域变换到频率域，并计算功率谱，然后对功率谱的每个采 样点值采取查表的方式取对数，得到对数功率谱，把地震记录段变换成复赛 谱，并把复赛谱截断，然后进行地表一致性分解。
本发明还提供如下解决方案
地表一致性反褶积的方法，采用以下步骤
1) 采集叠前炮集或CMP地震记录；
2) 用常规方法截取地震记录， 一般选择地震记录信噪比较高的时窗段；
3) 用常规的方法对地震记录段做傅立叶变换，把地震记录从吋间域变换
到频率域，并计算功率谱；
4) 对功率谱的每个采样点值取对数，得到对数功率谱；
5) 对得到的对数功率谱做反傅立叶变换，得到复赛谱，
6) 取复赛谱零点附近的记录，时间长度与人工填写的自相关长度一致；
7) 在复赛域分解，对复赛谱用高斯一塞德尔方法进行地表一致性分解， 分解后得到各个分量的响应；
8) 对各个分量响应做傅氏变换，得到各个分量的对数功率谱；
9) 对每个分量的对数功率谱分别取指数，得到功率谱；
10) 对功率谱做反傅氏变换，得到各个分量的自相关函数；
11) 分别解各个分量自相关函数组成的托布里兹方程组，求出预测算子；
12) 分别在地震道上应用各个分量的反褶积算子，用常规的方法处理得到 高分辩率的地震资料，并消除各炮点和检波点之间的差别。
所述的信噪比较高的时窗段同相轴的位置在0.8s 3s之间。 所述的各个分量是炮点、检波点、CMP和炮检距。
所述的地震记录段的功率谱取对数，采用造表杳表的方式。 本发明由于在复赛谱域地震子波一般集中在复赛谱的零点附近，所以截断零点附近的复赛谱进行地表一致性分解，这样既充分保留了地震子波的信 息，也节省了内存空间，同时由于复赛谱被截断，所要进行地表一致性分解 的釆样点减少，提高了分解的效率。本发明后处理时内存大约节省到原来的
1/5~1/10，效率大约节省到原来的1/5 1/15。效果要好于常规的地表一致性反褶积。
图1是本发明在单炮上面的应用效果图1 (a)是单炮记录；
图1 (b)是一般的地表一致性反褶积方法；
图l (C)是本发明地表一致性反褶积方法；
图2是地震记录自相关结果图2 (a)是原始地震记录的自相关函数；
图2 (b)是经过地表一致性反褶积后结果图2 (c)是经过本发明处理后结果图。
具体实施方案
本发明先采集叠前炮集或CMP地震记录，用常规方法截取地震id录，一 般选择地震记录信噪比较高的时窗段，用常规的方法对地震记录段做傅立叶 变换，把地震记录从时间域变换到频率域，并计算功率谱，然后对功率谱的 每个采样点值采取査表的方式取对数，得到对数功率谱，把地震记录段变换 成复赛谱，并把复赛谱截断，然后进行地表一致性分解。
采用以下具体步骤
1)采集叠前炮集或CMP地震记录。
可设地震记录为;x(1)， t= 1，..、N， 1，…，N为采样点，
2)用常规方法截取地震记录， 一般选择地震记录信噪比较高的时窗段，如在0.8s 3: 之间。
设截取后的地震记录为x(,)， / = ^,..，仏，W,，^为截取地震记录的起止采 样点；
3)用常规的方法对地震记录段做傅立叶变换，把地震记录从时间域变换 到频率域，并计算功率谱；
如所述的对雄)做傅立叶变换<formula>complex formula see original document page 8</formula>，其中的功率谱为
4) 造对数表，只计算一次对数值，并把自变量和相应的对数值制成对数 表，以后需要计算对数值得时候根据自变量直接査找对数值，可以避免相同 自变量的重复计算，然后对功率谱的每个采样点值采取查表的方式取对数， 得到对数功率谱；设l义(/f的对数功率谱表示为14^M2
5) 对得到的对数功率谱做反傅立叶变换，得到复赛谱，设<formula>complex formula see original document page 8</formula>为复赛谱；
6) 取复赛谱零点附近的记录，时间长度与人工填写的自相关长度一致， 取邻)零点附近的记录邻)，,=1,-.，W';
7) 在复赛域分解，对复赛谱用高斯一塞德尔方法进行地表一致性分解， 分解后得到炮点项、检波点项、CMP项和炮检距项的响应，如；<formula>complex formula see original document page 8</formula>
其中外)为炮点项， )为检波点项，^)为CMP项，/(0为炮检距项，(n)为第n次迭代的结果，(n-1)为第n-1次迭代的结果，，i为炮号，j为炮集i中的所有检波点号；
8) 对各个分量响应做傅氏变换，得到各个分量的对数功率谱如；
<formula>complex formula see original document page 9</formula>为炮点项的对数功率谱 <formula>complex formula see original document page 9</formula>为检波点项的对数功率谱<formula>complex formula see original document page 9</formula>为cmp项的对数功率谱<formula>complex formula see original document page 9</formula>为炮检距项的对数功率谱
9) 对每个分量的对数功率谱分别取指数，得到功率谱，如；
以炮点项为例<formula>complex formula see original document page 9</formula>为炮点项的功率谱
10) 对功率谱做反傅氏变换，得到各个分量的自相关函数，如；
以炮点项为例<formula>complex formula see original document page 9</formula>为炮点项的自相关函数
11) 分别解由炮点、检波点、cmp和炮检距自相关函数组成的托布里兹方 程组，求出预测算子；
如解托布里兹方程组<formula>complex formula see original document page 9</formula>
12)分别在地震道上应用各个分量的反褶积算子。得到高分辩率的地震 资料，并消除各炮点和检波点之间的差别。
本发明实施例
1) 采集叠前炮集地震记录，总炮数为210炮。
2) 选择时窗2000 3000ms对地震记录截断
3) 对地震记录段变换到复赛域，对复赛域地震记录进行地表一致性分
解，得到各个分量的自相关函数。
4) 然后分别解自相关函数组成的托布里兹方程组，得到各个分量的反褶 积算子。
5) 在地震记录上分别应用这些反褶积算子。
图1 (c)是本实施例效果，明显好于图1 (b) —般地表一致性反褶积的 效果。
通常对地表一致性反褶积效果的监控， 一般都是利用自相关函数。本发 明另一实施例处理后对地震道做自相关处理，对自相关后的数据按炮集叠加， 炮集自相关结果，用来判别地表一致性反褶积的效果。
图2是实际地震记录地表一致性反褶积前后自相关的对比。显示的自相 关函数是炮集上各道自相关函数之和。图2 (a)是原始地震记录的自相关函 数，炮间差别很大，分辩率并不高。经过地表一致性反褶积(限频80Hz)后 图2 (b)分辩率有所提高，但是地表一致性并不是很好，而且有高频分量， 复赛谱域地表一致性反褶积结果图2 (c)，各炮之间的差异很小，分辩率提 高。
权利要求
1、一种地表一致性反褶积的方法，采集叠前炮集或CMP地震记录，用常规方法截取地震记录，一般选择地震记录信噪比较高的时窗段，用常规的方法对地震记录段做傅立叶变换，把地震记录从时间域变换到频率域，并计算功率谱，其特征在于对功率谱的每个采样点值采取查表的方式取对数，得到对数功率谱，把地震记录段变换成复赛谱，并把复赛谱截断，然后进行地表一致性分解。
2、 根据权利要求1所述的一种地表一致性反褶积的方法，其特征在于 采用以下步骤1 )采集叠前炮集或CMP地震记录；2) 用常规方法截取地震记录， 一般选择地震记录信噪比较高的时窗段；3) 用常规的方法对地震记录段做傅立叶变换，把地震记录从时间域 变换到频率域，并计算功率谱；4) 对功率谱的每个采样点值采取査表的方式取对数，得到对数功率谱；5) 对得到的对数功率谱做反傅立叶变换，得到复赛谱，6) 取复赛谱零点附近的记录，时间长度与人工填写的自相关长度一致；7) 在复赛域分解，对复赛谱用高斯一塞德尔方法进行地表一致性分 解，分解后得到各个分量的响应；8) 对各个分量响应做傅氏变换，得到各个分量的对数功率谱；9) 对每个分量的对数功率谱分别取指数，得到功率谱；10) 对功率谱做反傅氏变换，得到各个分量的自相关函数；11) 分别解各个分量自相关函数组成的托布里兹方程组，求出预测算子；12) 分别在地震道上应用各个分量的反褶积算子，用常规的方法处理 得到高分辩率的地震资料，并消除各炮点和检波点之间的差别。
3、 根据权利要求2所述的一种地表一致性反褶积的方法，其特征在于步骤2)所述的信噪比较高的时窗段同相轴的位置在0.8s 3s之间。
4、 根据权利要求2所述的一种地表一致性反褶积的方法，其特征在 于所述的各个分量是炮点、检波点、CMP和炮检距。
5、 根据权利要求2所述的一种地表一致性反褶积的方法，其特征在 于所述的地震记录段的功率谱取对数，采用造表査表的方式。
全文摘要
本发明是地震资料处理过程中地表一致性反褶积的方法，采集叠前炮集或CMP地震记录，用常规方法截取地震记录，一般选择地震记录信噪比较高的时窗段，用常规的方法对地震记录段做傅立叶变换，把地震记录从时间域变换到频率域，并计算功率谱，其特征在于对功率谱的每个采样点值采取查表的方式取对数，得到对数功率谱，把地震记录段变换成复赛谱，并把复赛谱截断，然后进行地表一致性分解。本发明充分保留地震子波的信息，节省内存空间，所要进行地表一致性分解的采样点减少，后处理内存节省，可得到高分辩率的地震资料，并消除各炮点和检波点之间的差别。
文档编号G01V1/28GK101201406SQ200610165019
公开日2008年6月18日 申请日期2006年12月12日 优先权日2006年12月12日
发明者任晓乔, 周兴元 申请人:中国石油天然气集团公司;中国石油集团东方地球物理勘探有限责任公司