一种地震记录变时窗自动增益控制方法
【专利摘要】一种地震记录变时窗自动增益控制方法，属于地震资料数据处理领域。该方法首先根据给定的最浅层和最深层的时窗长度，计算出各个时窗的长度，在此基础上计算出各个时窗内各个点的平均振幅值，进而得到每个时窗内每个采样点的均衡值。模型数据的实验结果显示，基于变时窗自动增益控制处理的炮集记录，既可保证浅层的弱反射得到较好的能量增益，又可最大限度的保留深层的动力学特征，是一种较为理想的能量增益处理方法。
【专利说明】
一种地震记录变时窗自动増益控制方法
技术领域
[0001] 本发明属于地震资料数据处理领域，是一种既可增强浅层的弱反射信息又可保证 深层反射动力学特征的有效均衡方法。
【背景技术】
[0002] 由于地震波的几何扩散效应以及在复杂的地质构造条件下地层吸收衰减的影响， 会使得地震记录上浅层的反射波能量较强，而深层的反射波能量较弱，浅层与深层能量差 异较大，并且随着炮检距的增大，道与道之间相同深度也存在能量不均衡的情况。炮集记录 的能量不均衡情况会严重影响深层的成像精度，因此地震数据处理时需要对地震记录进行 能量均衡处理。
[0003] 能量均衡是应用一定的数学处理方法对地震记录进行处理，使得浅、中、深层的反 射同相轴能量达到均衡。当前的能量均衡方法主要有两种:一种是真振幅恢复，其主要处理 步骤包括波前扩散能量补偿、地层吸收能量补偿和地表一致性能量调整等，理论上若能够 得到地下介质的真实速度、Q值等参数，则真振幅恢复处理将可严格地按地震波的动力学特 征对地震记录进行能量均衡处理，但是实际工程中地下介质的速度、Q值等参数都是未知 的，因此真振幅恢复的能量均衡效果受到影响；另一种常用的能量增益方法是自动增益控 制(AGC)，其应用时窗内的平均能量值对时窗内每个点的能量进行均衡，并通过时窗的滑动 处理实现记录内所有点的均衡处理。
[0004] 常规的自动增益控制采用固定时窗，其在时窗选取时遇到困难:若时窗太大，浅层 的弱反射难以得到有效增益，而时窗太小又会破坏深层的动力学特征，不利于后续的地质 解释。因此本文提出了一种基于变时窗的地震记录自动增益控制方法，其在浅层利用小时 窗，深层利用大时窗，时窗随着深度线性增加，以达到浅中深层有效增益的目的。数值试验 结果显示，相比于传统的自动增益控制方法，基于变时窗的自动增益控制方法既可保证浅 层的弱反射得到较好的能量增益，又可最大限度地保留深层的动力学特征，是一种较为理 想的能量增益处理方法。
【发明内容】
[0005] 本发明要解决的技术问题在于提供一种地震记录变时窗自动增益控制方法，以保 证浅层的弱反射得到较好的能量增益，并最大限度的保留深层的动力学特征。对于每炮地 震记录，首先在确定最浅层和最深层的时窗长度的基础上，计算出时窗移动的线性增量;然 后在计算过程中时窗长度由浅至深线性增加，以满足浅层和深层的不同增益需要。
[0006] 本发明采取以下技术方案：
[0007] -种地震记录变时窗自动增益控制方法，其特征在于它具体包括以下步骤：
[0008] (1)设记录长度为T，给定最浅层时窗长度Ti和最深层时窗长度T2,按下列步骤计算 从浅层到深层每个时窗的长度；
[0009 ]①利用式(1)计算每次时窗移动的线性增量1;
[0010]
(1)
[0011]式中，N2和Νι分别是最深层和最浅层时窗中心点的米样点序号；
[0012] ②利用式(2)计算每个时窗的长度；
[0013]
(2)
[0014] 式中，int{}表示取整运算，At为采样间隔，I表示第I个时窗（取值范围为1~N2-Ni)；
[0015] (2)计算地震记录每一道每个采样点的平均振幅值，按以下步骤计算：
[0016] ①在第一个时窗内，基于式(3)计算时窗中心点的平均振幅值；
[0017]
(3)
[0018] 式中，七,表示第j道η时刻的平均振幅值，Xj,k表示第j道k时刻的采样值，m= int {L
[1]/2}为时窗长度的一半；
[0019] ②将时窗向下滑动一个米样点，时窗中心点位置不变，时窗长度调整为L[I]，重复 ①操作，从而得到第二个时窗中心点的平均振幅值。
[0020] ③随着时窗不断向下滑动，最终可得到IV2至T-IV2范围内所有点的平均振幅值；
[0021] ④对于tCh/2的时间采样点，令其等于IV2点的平均振幅值;对于t>T_T2/2的时间 采样点，令其等于T-IV2点的平均振幅值，这样即可得到该道每个采样点的平均振幅值。
[0022] ⑤根据步骤①~④求取地震记录每一道每个采样点的平均振幅值；
[0023] (3)利用式(4)计算地震记录每一道每个采样点的加权系数；
[0024]
(4)
[0025] 式中，参数Q为均衡参数，可控制输出数据的大小，选择1000;
[0026] (4)利用式(5)对每个采样点进行振幅均衡计算，即可得到该地震记录每一道每个 采样点的均衡振幅值；
[0027] Y^ = PhkxJhk 〇 (5)
[0028] 本发明与现有技术相比的有益效果：
[0029] 本发明提出的一种地震记录变时窗自动增益控制方法，在变时窗长度的基础上计 算记录中每道每个采样点的加权系数，然后利用这些加权系数求取记录中各点的均衡振幅 值。相比于传统的自动增益控制方法，基于变时窗的自动增益控制方法，既可保证浅层的弱 反射得到较好的能量增益，又可最大限度地保留深层的动力学特征。
【附图说明】
[0030] 图1水平层状速度模型；
[0031] 图2单炮地震记录；
[0032]图3变时窗自动增益控制处理的单炮地震记录；
[0033]图4 100ms固定时窗自动增益控制处理的单炮地震记录；
[0034]图5 1000ms固定时窗自动增益控制处理的单炮地震记录；
[0035]图6基于变时窗的自动增益控制处理流程图。
【具体实施方式】
[0036]下面通过实施例来对本发明的技术方案作进一步解释，但本发明的保护范围不受 实施例任何形式上的限制。
[0037] 实施例1
[0038] 本发明提出的地震记录变时窗自动增益控制方法，具体的实施过程主要分为以下 四步：（1)在确定地震记录最浅层和最深层的时窗长度的基础上，计算出每个时窗的长度； (2)计算地震记录每一道每个采样点的平均振幅值；（3)计算地震记录每一道每个采样点的 加权系数;（4)计算地震记录每一道每个采样点的均衡振幅值。
[0039] 本论文基于一个7层的水平层状模型（速度模型如图1所示，具体参数见表1，为得 到浅层弱反射信号，第二与第三层以及第三与第四层速度间隔较小），应用有限差分正演模 拟方法模拟得到一个单炮地震记录，然后应用该单炮地震记录进行能量均衡处理实验。该 单炮地震记录共有300道，道间距为5m，采样间隔为0.5ms，记录时间6000ms (单炮记录如图2 所示）。
[0040] 表1速度模型参数表
[0041]
[0042] 下面将详细阐述本发明的具体实施过程：
[0043] (1)给定最浅层时窗长度为Ti = 100ms，最深层时窗长度为T2 = 1000ms，按以下步骤 计算从浅层到深层每个时窗的长度：
[0044] ①利用式(6)计算得到每次移动时窗时时窗长度的线性增量1 = 0.08ms;
[0045]
⑴.
[0046] 式中，他=11000和见=100分别是最深层和最浅层时窗中心点的采样点序号；
[0047] ②利用式(7)计算每个时窗的长度(采样点个数）。
[0048]
(2)
[0049] 式中，int{}表示取整运算，At为采样间隔，I表示第I个时窗（取值范围为1~ 10900)；
[0050] (2)计算地震记录每一道每个采样点的平均振幅值，按以下步骤计算：
[0051] ①在第一个时窗内，基于式(8)计算时窗中心点的平均振幅值。
[0052]
(3)
[0053] 式中，為表示第j道η时刻的平均振幅值，Xj, k表示第j道k时刻的采样值，m = int {L [1]/2}为时窗长度的一半；
[0054] ②将时窗向下滑动一个采样点，时窗中心点位置不变，时窗长度调整为L[I]，重复 ①操作，从而得到第二个时窗中心点的平均振幅值。
[0055] ③随着时窗不断向下滑动，最终可得到50ms至5500ms范围内所有点的平均振幅 值；
[0056] ④对于t〈50ms的时间采样点，令其等于50ms处的平均振幅值;对于t>5500ms的时 间采样点，令其等于5500ms处的平均振幅值，这样即可得到该道每个采样点的平均振幅值。
[0057] ⑤根据步骤①~④求取地震记录每一道每个采样点的平均振幅值。
[0058] (3)利用式(9)计算地震记录每一道每个采样点的加权系数；
[0059]
⑷
[0060] 式中，参数Q为均衡参数，可控制输出数据的大小，选择1000;
[0061] (4)利用式(5)对每个采样点进行振幅均衡计算，即可得到该地震记录每一道每个 采样点的均衡振幅值。
[0062] Y^^=PltxX" (5)
[0063] 以上步骤的流程图如图6所示。
[0064] 为了说明本发明所述方法的有效性，这里与基于传统的自动增益控制方法得到的 单炮地震记录(如图4、5所示)进行了对比。由图4可知，经100ms小时窗自动增益处理后，浅 层的弱反射能量得到较好的增益效果，深层的能量也得到较好的增益。但对比速度模型可 以发现，在深层多出了很多虚假信号(如图4箭头所示），这是因为小时窗包含的地震信号较 少，在深层作用时，把层间能量均衡到相近的级别，破坏了深层所需要的动力学特征，均衡 出了许多虚假同相轴。由图5可知，经过1000ms时窗均衡后，深层能量得到较好恢复并未出 现虚假信号，但浅层的弱反射同相轴却没有得到有效增益(见图5中圆圈标定部分）。由此可 见，虽然大时窗能保证深层的动力学特征，但因为大时窗中包含过多的强能量信号，不利于 浅层的弱反射同相轴的均衡。由图3可知，经变时窗自动增益控制处理后，浅层的弱反射得 到了较好的增益，而深层的能量也得到了恢复，并且由于深层的时窗较大，深层地震波的动 力学特征得到了保证，未出现虚假同相轴，满足了地震记录在浅中深层不同的均衡要求。 [0065]由此可见，基于地震记录变时窗自动增益控制方法处理的炮集记录，在保证浅层 的弱反射得到较好的能量增益的基础上，又最大限度地保留了深层的动力学特征，是一种 理想的炮集记录均衡方法。
【主权项】
1. 一种地震记录变时窗自动增益控制方法，其特征在于它具体包括W下步骤： (1) 设记录长度为T，给定最浅层时窗长度Τι和最深层时窗长度T2，按下列步骤计算从浅 层到深层每个时窗的长度； ① 利用式(1)计算每次时窗移动的线性增量1;(1) 式中，化和化分别是最深层和最浅层时窗中屯、点的采样点序号； ② 利用式(2)计算每个时窗的长度；(2) 式中，int{}表示取整运算，Δ t为采样间隔，I表示第I个时窗，取值范围为1~化-Ni; (2) 计算地震记录每一道每个采样点的平均振幅值，按W下步骤计算： ① 在第一个时窗内，基于式(3)计算时窗中屯、点的平均振幅值；(3) 式中，表示第j道η时刻的平均振幅值，X^k表示第j道k时刻的采样值，m=int{L[I]/ 2}为时窗长度的一半； ② 将时窗向下滑动一个采样点，时窗中屯、点位置不变，时窗长度调整为L[I]，重复①操 作，从而得到第二个时窗中屯、点的平均振幅值； ③ 随着时窗不断向下滑动，最终可得到Ti/2至T-T2/2范围内所有点的平均振幅值； ④ 对于t<Ti/2的时间采样点，令其等于Ti/2点的平均振幅值;对于t〉T-T2/2的时间采样 点，令其等于T-T2/2点的平均振幅值，运样即可得到该道每个采样点的平均振幅值； ⑤ 根据步骤①~④求取地震记录每一道每个采样点的平均振幅值； (3) 利用式(4)计算地震记录每一道每个采样点的加权系数；(4) 式中，参数Q为均衡参数，可控制输出数据的大小，选择1000; (4) 利用式(5)对每个采样点进行振幅均衡计算，即可得到该地震记录每一道每个采样 点的均衡振幅值；(5)
【文档编号】G01V1/36GK105974467SQ201610289295
【公开日】2016年9月28日
【申请日】2016年5月5日
【发明人】夏冬明, 李金山, 宋鹏, 谭军, 姜秀萍, 赵波, 张晓波, 钟梦轩
【申请人】中国海洋大学