【技术实现步骤摘要】
本专利技术属于地球物理勘探，涉及地震资料处理技术,具体地说是一种适用于黄土塬区的叠前振幅补偿方法。
技术介绍
1、鄂尔多斯盆地黄土塬区表层大面积覆盖第四系巨厚黄土、黄土疏松干燥对地震波的吸收衰减极为严重，造成黄土塬区地震勘探资料频带窄、低信噪比的现状；黄土塬区地震资料的提频主要依靠叠前反褶积，李庆忠在《地震勘探分辨率与信噪比谱的关系》一文中指出：“如果信噪比谱很理想，即有效频宽(即信噪比大于1的那些频率成份)很宽，那么可以把子波的视周期通过反褶积将其压缩成接近一个窄脉冲，分辨率可以提高”，而目前鄂尔多斯黄土山地区低信噪比、频宽窄的原始资料难以获得较好的反褶积效果(如附图1)；因此需对资料进行叠前振幅补偿，提升高频有效反射的能量，拓宽频带。
2、常规的振幅补偿方法在应用时利用全频带的单炮进行补偿效果的质控，而全频带的单炮主要反映了低频信息时间域的能量变化，忽略了黄土对高频成份的吸收衰强于低频，高频衰减远快于低频的特点，往往低频补偿较好、而高频补偿严重不足(如附图1b，大于25hz以上剖面深层出现能量补偿不足的现象)，难以拓宽有效频带，反褶积的效果不理想,所以需要分频的振幅处理。
3、目前在公开的文献和报告中针对时间─频率─振幅相关的振幅补偿方法主要有反q滤波、时频分析法、反q偏移法，以上方法在黄土塬地区的应用效果不理想，这是因为：
4、反q滤波及反q偏移法需要一个高精度的q值，目前鄂尔多斯盆地尝试采用表层q(点：双井微测井，面：层析反演速度场)加深层q(井资料约束的叠加速度场) 的方法建立v-q
5、时频分析法，目前较常用的分频振幅补偿通常采用余弦函数在频率域构造分频滤波函数对数据分频，而后对各个频段吸收衰减函数进行多项式拟合，重构数据，得到最终补偿输出结果。由于该方法采用等间隔分频的方法，在频带划分时有可能把吸收衰减较少段的低频与吸收衰减强的高频分到一段，造成统计不准确；另外针对分频的资料自动提取补偿因子，黄土塬低信噪比资料在补偿因子的求取时由于随机噪音的混入，导致补偿因子提取的不准确；谱白化也可以进行频率-振幅的处理，但其主要方法是对期望的频带进行时变增益的处理，对其他频带进行压制，保真保幅性差。
技术实现思路
1、为解决现有技术中存在的以上不足，本专利技术旨在提供一种适用于黄土塬区的叠前振幅补偿方法，以达到有效拓宽频带，通过反褶积达到提升分辨率的目的。
2、为实现上述目的，本专利技术所采用的技术方案如下：一种适用于黄土塬区的叠前振幅补偿方法，包括以下步骤：
3、步骤1：对地震数据进行去噪；
4、步骤2：对去噪的数据进行叠加；
5、步骤3：依据叠加剖面的品质对数据进行分区或不分区；
6、步骤4：各分区数据进行叠加剖面的频谱分析；
7、步骤5：依据各分区数据的叠加剖面频谱分析结果对频带进行划分；
8、步骤6：确定各分区的低频补偿因子；
9、步骤7：计算各分区高频各分频段的补偿因子增量；
10、步骤8：各分区域、分频数据进行分频的振幅补偿；
11、步骤9：数据重构。
12、作为本专利技术的限定，步骤3中：根据叠加剖面的分辨率以及三维的数据，结合剖面提取的属性图、近地表结构图，进行数据分区。
13、作为本专利技术的进一步限定，步骤4中：
14、固定剖面分析时窗(t2-t1)，且(t2-t1)>300ms；
15、将(t2-t1)作为一个传播周期，利用公式1,计算出振幅经过(t2-t1)的衰减系数表达式coef(t2-t1)的值，
16、
17、式中，n2为t2时的补偿因子，n1为t1时刻的补偿因子。
18、作为本专利技术的再进一步限定，步骤5中：固定时窗(t2-t1)进行频谱分析后，频谱上能量处于保持状态的低频范围作为无吸收衰减的标准，从能量开始下降段作为要增加补偿的段；
19、在进行振幅补偿的同时，对高频段的补偿进行补偿控制，对叠加剖面的高通滤波扫描，预估可以进行能量恢复的最大有效频率，针对大于该频率的频段不做补偿。
20、作为本专利技术的更进一步限定，步骤5中：
21、在频段划分时，将频谱上每下降固定的分贝β，划为一个频段；
22、根据公式2，计算β：
23、
24、a1、a2为相邻波峰的振幅值，其计算公式为：
25、a2＝a1·e-αr·ω(t)—————————————————公式3
26、式中r为传播距离，α为吸收衰减系数，ω(t)为波动函数；如果是球面波，在a1前乘以1/r。
27、作为本专利技术的再进一步限定，黄土山地区取β＝10分贝。
28、作为本专利技术的进一步限定，步骤6中：对分离出作为补偿标准的低频部分，利用吸收衰减补偿公式1，结合剖面或单炮补偿的能量曲线,测试出低频基础补偿因子n0。
29、作为本专利技术的更进一步限定，步骤7中：(n1-n0)＝0.6，(n2-n1)＝0.6， (n3-n2)＝0.6。
30、作为本专利技术的其它限定，步骤5至步骤9的表达公式为：
31、xi(t)＝x(t)×hi(t)———————————————公式4
32、
33、
34、hi(t)为第i个滤波器的滤波算子，x(t)为原始记录，xi(t)为分频数据，wi为加权因子，zi(t)为补偿后的分频数据，coef(xi(t))为分频数据振幅随时间变化的吸收衰减系数表达式，t为频谱分析的时窗大小，nxi为各分频数据对应的吸收衰减补偿参数，y(t)为补偿后的重构记录。
35、由于采用了上述的技术方案，本专利技术与现有技术相比，所取得的有益效果是：
36、(1)本专利技术是针对黄土塬区地震资料横向变化快、对地震资料高频成份吸收衰减剧烈的特点(如附图3)，提出的一种基于高低频吸收衰减差异性的叠前振幅补偿技术，根据低频吸收衰减程度较小的特点(如附图3)，在时间域对低频简单进行补偿即可获得较好的能量补偿效果(如附图2中小于10hz的低频补偿效果)，可将低频作为没有近地表吸收衰减的模型标准，高频成份向该标准方向进行能量的补偿，即可实现解决近地表对高频吸收衰减的问题，以此实现在保真保幅的基础上有效拓宽频带，通过反褶积达到提升分辨率的效果。
37、(2)本专利技术在叠前可以对黄土塬资料不同频带进行有效的振幅补偿，保真度高、保幅性好，可有效的拓宽频带、提升高频反射的能量，通过反褶积提频最终得到高分辨率的地震剖面(如附图8,9,10)；可以时本文档来自技高网...
【技术保护点】
1.一种适用于黄土塬区的叠前振幅补偿方法，其特征在于包括以下步骤：
2.根据权利要求1所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤3中：根据叠加剖面的分辨率以及三维的数据，结合剖面提取的属性图、近地表结构图，进行数据分区。
3.根据权利要求1或2所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤4中：
4.根据权利要求3所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤5中：固定时窗(t2-t1)进行频谱分析后，频谱上能量处于保持状态的低频范围作为无吸收衰减的标准，从能量开始下降段作为要增加补偿的段；
5.根据权利要求4所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤5中：
6.根据权利要求5所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤5中：黄土山地区取β＝10分贝。
7.根据权利要求4-6中任意一项所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤6中：对分离出作为补偿标准的低频部分，利用吸收衰减补偿公式1，结合剖面或单炮补偿的能量曲线,测试出低
8.根据权利要求7所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤7中：(n1-n0)＝0.6，(n2-n1)＝0.6，(n3-n2)＝0.6。
9.根据权利要求1、2、4、5、6、8中任意一项所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤5至步骤9的表达公式为：
...【技术特征摘要】
1.一种适用于黄土塬区的叠前振幅补偿方法，其特征在于包括以下步骤：
2.根据权利要求1所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤3中：根据叠加剖面的分辨率以及三维的数据，结合剖面提取的属性图、近地表结构图，进行数据分区。
3.根据权利要求1或2所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤4中：
4.根据权利要求3所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤5中：固定时窗(t2-t1)进行频谱分析后，频谱上能量处于保持状态的低频范围作为无吸收衰减的标准，从能量开始下降段作为要增加补偿的段；
5.根据权利要求4所述的一种适用于黄土塬区的叠前振幅补偿方法，其特征在于步骤5...
【专利技术属性】
技术研发人员：王春，林彧涵，宋方述，
申请(专利权)人：中国石油天然气集团有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。