《一种振幅补偿因子的反演方法.pdf》由会员分享,可在线阅读,更多相关《一种振幅补偿因子的反演方法.pdf(12页完整版)》请在专利查询网上搜索。
1、(10)申请公布号 (43)申请公布日 (21)申请号 201410806081.4 (22)申请日 2014.12.22 G01V 1/28(2006.01) G01V 1/36(2006.01) (71)申请人 中国石油天然气集团公司 地址 100007 北京市东城区东直门北大街 9 号 申请人 中国石油集团东方地球物理勘探有 限责任公司 (72)发明人 张固澜 张建军 李飞 吴俊军 张天仓 (74)专利代理机构 北京三友知识产权代理有限 公司 11127 代理人 王涛 (54) 发明名称 一种振幅补偿因子的反演方法 (57) 摘要 本发明涉及一种振幅补偿因子的反演方法, 包括 : 获取垂。
2、直地震剖面数据、 地面地震数据, 并 获取得到每道垂直地震剖面数据对应的监控子 波记录 ; 获取每道垂直地震剖面数据的初至时间 一, 并根据垂直地震剖面数据、 地面地震数据、 监 控子波记录获取各深度点的第三波场的绝对能 量, 从而确定振幅补偿因子。本发明能够基于反 Q 滤波振幅补偿, 反演振幅补偿因子, 并用于相对保 幅处理。 (51)Int.Cl. (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书3页 说明书7页 附图1页 (10)申请公布号 CN 104502965 A (43)申请公布日 2015.04.08 CN 104502965 A 1/3 页 2 1.一种。
3、振幅补偿因子的反演方法, 其特征在于, 包括 : 获取垂直地震剖面数据、 地面地震数据, 并获取每道垂直地震剖面数据对应的监控子 波记录 ; 获取每道垂直地震剖面数据的初至时间一, 并通过所述初至时间一反演垂直地震剖面 的层速度 ; 根据所述垂直地震剖面数据生成频率波数谱, 并根据所述频率波数谱生成一下行波 ; 在所述下行波中, 根据所述初至时间一设置一时窗, 并将所述时窗外的信号清零, 生成 一第二波场 ; 根据所述第二波场生成一振幅谱, 根据所述振幅谱获取各个深度点和炮点间的平均 Q 值, 根据所述平均 Q 值确定相邻深度点之间的层 Q 值 ; 根据所述垂直地震剖面的层速度和所述层 Q 值。
4、, 生成所述垂直地震剖面的层速度和所 述层 Q 值的关系函数 ; 对所述地面地震数据进行噪音压制, 并利用共中心点道集进行速度谱分析, 获取一层 速度体 ; 根据所述关系函数和所述层速度体, 生成一层 Q 体 ; 在经过噪音压制的地面地震数据中, 选取目标区域不同信噪比的多个单炮数据或共中 心点道集数据 ; 设置不同的振幅补偿增益限, 通过所述层 Q 体和所述层速度体, 对所述单炮数据或共 中心点道集数据进行反 Q 滤波振幅补偿 ; 根据所述经过噪音压制的地面地震数据, 设置一第二时窗, 并确定一优化后的振幅补 偿增益限 ; 根据所述层 Q 体、 所述优化后的振幅补偿增益限、 所述层速度体对所。
5、述多个单炮数据 或共中心点道集数据进行反 Q 滤波振幅补偿 ; 根据反 Q 滤波振幅补偿后的多个单炮数据或共中心点道集数据、 所述第二时窗, 确定 一优化后的层 Q 体 ; 根据所述优化后的振幅补偿增益限、 所述优化后的层 Q 体、 所述层速度体、 对所述第二 波场进行反 Q 滤波振幅补偿, 生成一第三波场 ; 在所述第三波场中, 获取每个深度点的第三波场的相对能量 ; 在所述监控子波记录中, 获取每道监控子波记录的初至时间二, 并根据所述初至时间 二设置一第三时窗, 根据所述第三时窗确定所述每个深度点的激发子波的能量 ; 根据各深度点的所述第三波场的相对能量和所述激发子波的能量, 获取得到各。
6、深度点 的第三波场的绝对能量 ; 根据所述初至时间一、 各深度点的第三波场的绝对能量确定振幅补偿因子。 2.根据权利要求 1 所述的振幅补偿因子的反演方法, 其特征在于, 所述根据所述垂直 地震剖面数据生成频率波数谱, 并根据所述频率波数谱生成一下行波, 包括 : 对所述垂直地震剖面数据进行傅里叶变换, 生成频率波数谱 ; 根据所述频率波数谱进行波场分离, 生成一下行波。 3.根据权利要求 2 所述的振幅补偿因子的反演方法, 其特征在于, 所述根据所述第二 波场生成一振幅谱, 根据所述振幅谱获取各个深度点和炮点间的平均 Q 值, 根据所述平均 Q 权 利 要 求 书 CN 104502965 。
7、A 2 2/3 页 3 值确定相邻深度点之间的层 Q 值, 包括 : 对所述第二波场进行傅里叶变换, 生成一振幅谱 ; 根据所述振幅谱进行谱比法反演, 获取得到各个深度点和炮点间的平均 Q 值 ; 将所述平均 Q 值进行平滑处理, 反演得到相邻深度点之间的层 Q 值。 4.根据权利要求 3 所述的振幅补偿因子的反演方法, 其特征在于, 所述根据所述垂直 地震剖面的层速度和所述层 Q 值, 生成所述垂直地震剖面的层速度和所述层 Q 值的关系函 数, 包括 : 以所述垂直地震剖面层速度为自变量、 层 Q 值为应变量, 利用最小二乘法拟合出垂直 地震剖面层速度和层 Q 值之间的关系函数 : 其中, 。
8、Qi-1,i为所述层 Q 值, V i-1,i为垂直地震剖面层速度, G(Vi-1,i) 为所述关系函数, an 为拟合系数, N 为自然数。 5.根据权利要求 4 所述的振幅补偿因子的反演方法, 其特征在于, 所述根据所述关系 函数和所述层速度体, 生成一层 Q 体, 包括 : 将所述关系函数中的垂直地震剖面层速度替换为所述层速度体, 生成所述层 Q 体。 6.根据权利要求 5 所述的振幅补偿因子的反演方法, 其特征在于, 所述根据所述经过 噪音压制的地面地震数据, 设置一第二时窗, 并确定一优化后的振幅补偿增益限, 包括 : 将所述第二时窗内的地震信号进行傅里叶变换得到第二振幅谱 ; 根据。
9、所述第二振幅谱, 确定所述第二时窗内的地震信号的分辨率 ; 对比反 Q 滤波振幅补偿前后地震信号的分辨率, 若按一增益限进行反 Q 滤波振幅补偿 后的所述第二时窗内的地震信号的分辨率大于等于反 Q 滤波振幅补偿之前的地震信号的 分辨率, 确定该增益限为所述优化后的振幅补偿增益限。 7.根据权利要求 6 所述的振幅补偿因子的反演方法, 其特征在于, 所述根据所述层 Q 体、 所述优化后的振幅补偿增益限、 所述层速度体对所述多个单炮数据或共中心点道集数 据进行反 Q 滤波振幅补偿, 包括 : 根据所述层 Q 体和多个不同的比例系数 k, 确定得到多个不同的层 Q 体 ; 根据所述多个不同的层 Q 。
10、体、 所述优化后的振幅补偿增益限、 所述层速度体对所述多 个单炮数据或共中心点道集数据进行反 Q 滤波振幅补偿。 8.根据权利要求7所述的振幅补偿因子的反演方法, 其特征在于, 所述根据反Q滤波振 幅补偿后的多个单炮数据或共中心点道集数据、 所述第二时窗, 确定一优化后的层 Q 体, 包 括 : 在第二时窗内的地震信号的分辨率达到最高值时, 确定所述第二时窗内的地震信号的 分辨率对应的层 Q 体为所述优化后的层 Q 体。 9.根据权利要求 8 所述的振幅补偿因子的反演方法, 其特征在于, 所述在所述第三波 场中, 获取每个深度点的第三波场的相对能量, 包括 : 在所述第三波场中, 对每个深度点。
11、处的垂直地震剖面数据求取地震信号振幅的平方 和, 再进行开方运算, 生成每个深度点的第三波场的相对能量。 10.根据权利要求 9 所述的振幅补偿因子的反演方法, 其特征在于, 所述在所述监控子 权 利 要 求 书 CN 104502965 A 3 3/3 页 4 波记录中, 获取每道监控子波记录的初至时间二, 并根据所述初至时间二设置一第三时窗, 根据所述第三时窗确定所述每个深度点的激发子波的能量, 包括 : 求取所述第三时窗的地震信号振幅的平方和, 再进行开方操作, 生成所述每个深度点 的激发子波的能量。 11.根据权利要求 10 所述的振幅补偿因子的反演方法, 其特征在于, 所述根据各深度。
12、 点的所述第三波场的相对能量和所述激发子波的能量, 获取得到各深度点的第三波场的绝 对能量, 包括 : 将各深度点的所述第三波场的相对能量除以该深度点处激发子波的能量, 获取得到各 深度点的第三波场的绝对能量。 12.根据权利要求 11 所述的振幅补偿因子的反演方法, 其特征在于, 所述根据所述初 至时间一、 各深度点的第三波场的绝对能量确定振幅补偿因子, 包括 : 将各深度点处的初至时间一进行自然对数运算, 将所获得的自然对数设置为自变量 ; 将各深度点处第三波场的绝对能量进行自然对数运算, 将所获得的自然对数设置为应 变量 ; 将所述自变量和应变量通过最小二乘拟合方法进行直线拟合, 获取得。
13、到拟合后的直线 斜率值的绝对值 ; 确定所述拟合后的直线斜率值的绝对值为所述振幅补偿因子 ; 其中, 进行直线拟合的拟合函数为 : ln A A0-ln t t 为所述初至时间一, A 为所述第三波场的绝对能量, A0为常数, 为所述振幅补偿因 子。 权 利 要 求 书 CN 104502965 A 4 1/7 页 5 一种振幅补偿因子的反演方法 技术领域 0001 本发明涉及地震勘探数据处理技术领域, 特别涉及一种振幅补偿因子的反演方 法。 背景技术 0002 当前, 在地震勘探数据处理技术中, 为了采集垂直地震剖面法的地震资料, 可以将 检波器放置在井中, 直接接收到来自地层的下行波, 从。
14、而可以得到各深度的检波器接收到 的下行初至波的能量, 而在这个过程中存在地震波的能量损失。地震波的能量损失包括有 透射能量损失、 球面扩散能量损失、 以及与地层品质因子 Q 吸收衰减有关的能量损失。其 中, 透射能量损失和球面扩散能量损失与频率无关。 Q吸收衰减造成的能量损失与频率成正 比, 且高频的能量损失比低频更严重。 0003 随着当前对地震波保持振幅处理需求的越来越迫切, 振幅补偿变得尤为重要。目 前对地震信号的振幅补偿, 主要是利用零偏移距 VSP 资料的下行初至波, 统计其全频带的 能量, 并进行与时间有关的指数拟合, 从而得到振幅补偿因子, 并用于地面地震资料振幅补 偿。这种振幅。
15、补偿因子, 其实是透射能量损失, 球面扩散能量损失, 以及地层品质因子 Q 对 地震波的吸收衰减造成的能量损失的一种综合补偿。从理论上来讲, 当不考虑地震资料信 噪比和硬件的截断误差时, 反Q滤波可完全消除由于地层Q吸收对地震波造成的能量损失。 0004 随着 Q 提取技术的越来越成熟, 反 Q 滤波技术已经越来越广泛的应用到地面地震 资料处理当中。但由于反 Q 滤波振幅补偿, 受地震信号的信噪比和硬件的截断误差影响, 会 造成很强的不稳定性, 为控制反Q滤波的稳定性, 当前研究提出了增益控制的反Q滤波补偿 函数, 即反 Q 滤波振幅补偿函数的最大值为增益限。在一定的增益限约束的前提下, 增益。
16、控 制的反 Q 滤波, 只能恢复一部分由于 Q 吸收衰减造成能量损失。此时, 如果用零偏移距 VSP 资料全频信号反演的振幅补偿因子, 对反 Q 滤波后的地面地震资料进行振幅补偿, 势必会 造成能量补偿过度, 达不到地震资料保幅处理的目的。 0005 在当前的技术中, 还难以基于反 Q 滤波振幅补偿, 反演振幅补偿因子, 并用于相对 保幅处理。 发明内容 0006 本发明提供了一种振幅补偿因子的反演方法, 能够为了进行地震资料的振幅恢 复, 优化 Q 值来提高地震资料的分辨率, 而反演出一振幅补偿因子。 0007 为实现上述目的, 本发明提供了一种振幅补偿因子的反演方法, 该方法包括 : 00。
17、08 获取垂直地震剖面数据、 地面地震数据, 并获取每道垂直地震剖面数据对应的监 控子波记录 ; 0009 获取每道垂直地震剖面数据的初至时间一, 并通过所述初至时间一反演垂直地震 剖面的层速度 ; 0010 根据所述垂直地震剖面数据生成频率波数谱, 并根据所述频率波数谱生成一下行 说 明 书 CN 104502965 A 5 2/7 页 6 波 ; 0011 在所述下行波中, 根据所述初至时间一设置一时窗, 并将所述时窗外的信号清零, 生成一第二波场 ; 0012 根据所述第二波场生成一振幅谱, 根据所述振幅谱获取各个深度点和炮点间的平 均 Q 值, 根据所述平均 Q 值确定相邻深度点之间的。
18、层 Q 值 ; 0013 根据所述垂直地震剖面的层速度和所述层 Q 值, 生成所述垂直地震剖面的层速度 和所述层 Q 值的关系函数 ; 0014 对所述地面地震数据进行噪音压制, 并利用共中心点道集进行速度谱分析, 获取 一层速度体 ; 0015 根据所述关系函数和所述层速度体, 生成一层 Q 体 ; 0016 在经过噪音压制的地面地震数据中, 选取目标区域不同信噪比的多个单炮数据或 共中心点道集数据 ; 0017 设置不同的振幅补偿增益限, 通过所述层 Q 体和所述层速度体, 对所述单炮数据 或共中心点道集数据进行反 Q 滤波振幅补偿 ; 0018 根据所述经过噪音压制的地面地震数据, 设置。
19、一第二时窗, 并确定一优化后的振 幅补偿增益限 ; 0019 根据所述层 Q 体、 所述优化后的振幅补偿增益限、 所述层速度体对所述多个单炮 数据或共中心点道集数据进行反 Q 滤波振幅补偿 ; 0020 根据反 Q 滤波振幅补偿后的多个单炮数据或共中心点道集数据、 所述第二时窗, 确定一优化后的层 Q 体 ; 0021 根据所述优化后的振幅补偿增益限、 所述优化后的层 Q 体、 所述层速度体、 对所述 第二波场进行反 Q 滤波振幅补偿, 生成一第三波场 ; 0022 在所述第三波场中, 获取每个深度点的第三波场的相对能量 ; 0023 在所述监控子波记录中, 获取每道监控子波记录的初至时间二,。
20、 并根据所述初至 时间二设置一第三时窗, 根据所述第三时窗确定所述每个深度点的激发子波的能量 ; 0024 根据各深度点的所述第三波场的相对能量和所述激发子波的能量, 获取得到各深 度点的第三波场的绝对能量 ; 0025 根据所述初至时间一、 各深度点的第三波场的绝对能量确定振幅补偿因子。 0026 具体的, 所述根据所述垂直地震剖面数据生成频率波数谱, 并根据所述频率波数 谱生成一下行波, 包括 : 0027 对所述垂直地震剖面数据进行傅里叶变换, 生成频率波数谱 ; 0028 根据所述频率波数谱进行波场分离, 生成一下行波。 0029 具体的, 所述根据所述第二波场生成一振幅谱, 根据所述。
21、振幅谱获取各个深度点 和炮点间的平均 Q 值, 根据所述平均 Q 值确定相邻深度点之间的层 Q 值, 包括 : 0030 对所述第二波场进行傅里叶变换, 生成一振幅谱 ; 0031 根据所述振幅谱进行谱比法反演, 获取得到各个深度点和炮点间的平均 Q 值 ; 0032 将所述平均 Q 值进行平滑处理, 反演得到相邻深度点之间的层 Q 值。 0033 具体的, 所述根据所述垂直地震剖面的层速度和所述层 Q 值, 生成所述垂直地震 剖面的层速度和所述层 Q 值的关系函数, 包括 : 说 明 书 CN 104502965 A 6 3/7 页 7 0034 以所述垂直地震剖面层速度为自变量、 层 Q 。
22、值为应变量, 利用最小二乘法拟合出 垂直地震剖面层速度和层 Q 值之间的关系函数 : 0035 0036 其中, Qi-1,i为所述层Q值, Vi-1,i为垂直地震剖面层速度, G(Vi-1,i)为所述关系函数, an为拟合系数, N 为自然数。 0037 具体的, 所述根据所述关系函数和所述层速度体, 生成一层 Q 体, 包括 : 0038 将所述关系函数中的垂直地震剖面层速度替换为所述层速度体, 生成所述层 Q 体。 0039 具体的, 所述根据所述经过噪音压制的地面地震数据, 设置一第二时窗, 并确定一 优化后的振幅补偿增益限, 包括 : 0040 将所述第二时窗内的地震信号进行傅里叶变。
23、换得到第二振幅谱 ; 0041 根据所述第二振幅谱, 确定所述第二时窗内的地震信号的分辨率 ; 0042 对比反 Q 滤波振幅补偿前后地震信号的分辨率, 若按一增益限进行反 Q 滤波振幅 补偿后的所述第二时窗内的地震信号的分辨率大于等于反 Q 滤波振幅补偿之前的地震信 号的分辨率, 确定该增益限为所述优化后的振幅补偿增益限。 0043 具体的, 所述根据所述层 Q 体、 所述优化后的振幅补偿增益限、 所述层速度体对所 述多个单炮数据或共中心点道集数据进行反 Q 滤波振幅补偿, 包括 : 0044 根据所述层 Q 体和多个不同的比例系数 k, 确定得到多个不同的层 Q 体 ; 0045 根据所述。
24、多个不同的层 Q 体、 所述优化后的振幅补偿增益限、 所述层速度体对所 述多个单炮数据或共中心点道集数据进行反 Q 滤波振幅补偿。 0046 具体的, 所述根据反 Q 滤波振幅补偿后的多个单炮数据或共中心点道集数据、 所 述第二时窗, 确定一优化后的层 Q 体, 包括 : 0047 在第二时窗内的地震信号的分辨率达到最高值时, 确定所述第二时窗内的地震信 号的分辨率对应的层 Q 体为所述优化后的层 Q 体。 0048 具体的, 所述在所述第三波场中, 获取每个深度点的第三波场的相对能量, 包括 : 0049 在所述第三波场中, 对每个深度点处的垂直地震剖面数据求取地震信号振幅的平 方和, 再进。
25、行开方运算, 生成每个深度点的第三波场的相对能量。 0050 具体的, 所述在所述监控子波记录中, 获取每道监控子波记录的初至时间二, 并根 据所述初至时间二设置一第三时窗, 根据所述第三时窗确定所述每个深度点的激发子波的 能量, 包括 : 0051 求取所述第三时窗的地震信号振幅的平方和, 再进行开方操作, 生成所述每个深 度点的激发子波的能量。 0052 具体的, 所述根据各深度点的所述第三波场的相对能量和所述激发子波的能量, 获取得到各深度点的第三波场的绝对能量, 包括 : 0053 将各深度点的所述第三波场的相对能量除以该深度点处激发子波的能量, 获取得 到各深度点的第三波场的绝对能量。
26、。 0054 具体的, 所述根据所述初至时间一、 各深度点的第三波场的绝对能量确定振幅补 偿因子, 包括 : 说 明 书 CN 104502965 A 7 4/7 页 8 0055 将各深度点处的初至时间一进行自然对数运算, 将所获得的自然对数设置为自变 量 ; 0056 将各深度点处第三波场的绝对能量进行自然对数运算, 将所获得的自然对数设置 为应变量 ; 0057 将所述自变量和应变量通过最小二乘拟合方法进行直线拟合, 获取得到拟合后的 直线斜率值的绝对值 ; 0058 确定所述拟合后的直线斜率值的绝对值为所述振幅补偿因子 ; 0059 其中, 进行直线拟合的拟合函数为 : 0060 ln。
27、A A0-lnt 0061 t 为所述初至时间一, A 为所述第三波场的绝对能量, A0为常数, 为所述振幅补 偿因子。 0062 本发明实施例振幅补偿因子的反演方法, 能够为了进行地震资料的振幅恢复, 而 反演出一振幅补偿因子 ; 该反演出的振幅补偿因子与实际的地震资料处理完全结合, 能直 接用于地面地震资料保幅处理中, 反演方法具有很好的稳定性, 方便简洁, 运算效率较高。 附图说明 0063 为了更清楚地说明本发明实施例或现有技术中的技术方案, 下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍, 显而易见地, 下面描述中的附图仅仅是本 发明的一些实施例, 对于本领域普通技术人。
28、员来讲, 在不付出创造性劳动的前提下, 还可以 根据这些附图获得其他的附图。 0064 图 1 为本发明实施例提出的振幅补偿因子的反演方法流程图。 具体实施方式 0065 下面将结合本发明实施例中的附图, 对本发明实施例中的技术方案进行清楚、 完 整地描述, 显然, 所描述的实施例仅仅是本发明一部分实施例, 而不是全部的实施例。基于 本发明中的实施例, 本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例, 都属于本发明保护的范围。 0066 本发明提出一种振幅补偿因子的反演方法。如图 1 所示。该方法包括 : 0067 步骤 101、 获取垂直地震剖面数据、 地面地震数据, 。
29、并获取每道垂直地震剖面数据 对应的监控子波记录。 0068 步骤 102、 获取每道垂直地震剖面数据的初至时间一, 并通过所述初至时间一反演 垂直地震剖面的层速度。 0069 步骤 103、 根据垂直地震剖面数据生成频率波数谱, 并根据频率波数谱生成一下行 波。 0070 步骤 104、 在下行波中, 根据初至时间一设置一时窗, 并将时窗外的信号清零, 生成 一第二波场。 0071 步骤 105、 根据第二波场生成一振幅谱, 根据振幅谱获取各个深度点和炮点间的平 均 Q 值, 根据平均 Q 值确定相邻深度点之间的层 Q 值。 0072 步骤 106、 根据垂直地震剖面的层速度和层 Q 值, 生。
30、成垂直地震剖面的层速度和层 说 明 书 CN 104502965 A 8 5/7 页 9 Q 值的关系函数。 0073 步骤 107、 对地面地震数据进行噪音压制, 并利用共中心点道集进行速度谱分析, 获取一层速度体。 0074 步骤 108、 根据关系函数和层速度体, 生成一层 Q 体。 0075 步骤 109、 在经过噪音压制的地面地震数据中, 选取目标区域不同信噪比的多个单 炮数据或共中心点道集数据。 0076 步骤 110、 设置不同的振幅补偿增益限, 通过层 Q 体和层速度体, 对单炮数据或共 中心点道集数据进行反 Q 滤波振幅补偿。 0077 步骤 111、 根据经过噪音压制的地面。
31、地震数据, 设置一第二时窗, 并确定一优化后 的振幅补偿增益限。 0078 步骤 112、 根据层 Q 体、 优化后的振幅补偿增益限、 层速度体对多个单炮数据或共 中心点道集数据进行反 Q 滤波振幅补偿。 0079 步骤 113、 根据反 Q 滤波振幅补偿后的多个单炮数据或共中心点道集数据、 第二时 窗, 确定一优化后的层 Q 体。 0080 步骤 114、 根据优化后的振幅补偿增益限、 优化后的层 Q 体、 层速度体, 对第二波场 进行反 Q 滤波振幅补偿, 生成一第三波场。 0081 步骤 115、 在第三波场中, 获取每个深度点的第三波场的相对能量。 0082 步骤 116、 在监控子波。
32、记录中, 获取每道监控子波记录的初至时间二, 并根据初至 时间二设置一第三时窗, 根据第三时窗确定每个深度点的激发子波的能量。 0083 步骤 117、 根据各深度点的第三波场的相对能量和激发子波的能量, 获取得到各深 度点的第三波场的绝对能量。 0084 步骤 118、 根据初至时间一、 各深度点的第三波场的绝对能量确定振幅补偿因子。 0085 本发明实施例振幅补偿因子的反演方法, 能够为了进行地震资料的振幅恢复, 而 反演出一振幅补偿因子 ; 该反演出的振幅补偿因子与实际的地震资料处理完全结合, 能直 接用于地面地震资料保幅处理中, 反演方法具有很好的稳定性, 方便简洁, 运算效率较高。 。
33、0086 具体的, 上述步骤 103 可以通过如下方式实现 : 0087 对垂直地震剖面数据进行傅里叶变换, 生成频率波数谱。根据频率波数谱进行波 场分离, 生成一下行波。 0088 具体的, 上述步骤 105 可以通过如下方式实现 : 0089 对第二波场进行傅里叶变换, 生成一振幅谱。 根据振幅谱进行谱比法反演, 获取得 到各个深度点和炮点间的平均 Q 值。将平均 Q 值进行平滑处理, 反演得到相邻深度点之间 的层 Q 值。 0090 具体的, 上述步骤 106 可以通过如下方式实现 : 0091 以垂直地震剖面层速度为自变量、 层 Q 值为应变量, 利用最小二乘法拟合出垂直 地震剖面层速。
34、度和层 Q 值之间的关系函数 : 0092 0093 其中, Qi-1,i为层 Q 值, V i-1,i为垂直地震剖面层速度, G(Vi-1,i) 为关系函数, an为拟 合系数, N 为自然数。 说 明 书 CN 104502965 A 9 6/7 页 10 0094 具体的, 上述步骤 108 可以通过如下方式实现 : 0095 将关系函数中的垂直地震剖面层速度替换为层速度体, 生成层 Q 体。 0096 具体的, 上述步骤 111 可以通过如下方式实现 : 0097 将第二时窗内的地震信号进行傅里叶变换得到第二振幅谱。根据第二振幅谱, 确 定第二时窗内的地震信号的分辨率。对比反 Q 滤波。
35、振幅补偿前后地震信号的分辨率, 若按 一增益限进行反Q滤波振幅补偿后的第二时窗内的地震信号的分辨率大于等于反Q滤波振 幅补偿之前的地震信号的分辨率, 确定该增益限为优化后的振幅补偿增益限。 0098 该第二时窗长度小于等于 100 毫秒, 且包含主要目的层。 0099 具体的, 该第二振幅谱为 U(), 为频率, 地震子波 u(t) 的分辨率 R 定义为 : 0100 0101 在预设条件 : U(p) maxU(), 1 2下, 则有 : 0102 0103 其中, b 取不同的值, 代表不同的频率范围。 0104 具体的, 上述步骤 112 可以通过如下方式实现 : 0105 根据层 Q 。
36、体和多个不同的比例系数 k(k 1), 确定得到多个不同的层 Q 体。根据 多个不同的层 Q 体、 优化后的振幅补偿增益限、 层速度体对多个单炮数据或共中心点道集 数据进行反 Q 滤波振幅补偿。 0106 具体的, 上述步骤 113 可以通过如下方式实现 : 0107 在第二时窗内的地震信号的分辨率达到最高值时, 确定第二时窗内的地震信号的 分辨率对应的层 Q 体为优化后的层 Q 体。 0108 具体的, 上述步骤 115 可以通过如下方式实现 : 0109 在第三波场中, 对每个深度点处的垂直地震剖面数据求取地震信号振幅的平方 和, 再进行开方运算, 生成每个深度点的第三波场的相对能量。 0。
37、110 具体的, 上述步骤 116 可以通过如下方式实现 : 0111 求取第三时窗的地震信号振幅的平方和, 再进行开方操作, 生成每个深度点的激 发子波的能量。 0112 该第三时窗的长度小于 200 毫秒。 0113 具体的, 上述步骤 117 可以通过如下方式实现 : 0114 将各深度点的第三波场的相对能量除以该深度点处激发子波的能量, 获取得到各 深度点的第三波场的绝对能量。 0115 具体的, 上述步骤 118 可以通过如下方式实现 : 0116 将各深度点处的初至时间一进行自然对数运算, 将所获得的自然对数设置为自变 量。将各深度点处第三波场的绝对能量进行自然对数运算, 将所获得。
38、的自然对数设置为应 变量。将该自变量和应变量通过最小二乘拟合方法进行直线拟合, 获取得到拟合后的直线 斜率值的绝对值。确定拟合后的直线斜率值的绝对值为振幅补偿因子。 说 明 书 CN 104502965 A 10 7/7 页 11 0117 其中, 进行直线拟合的拟合函数为 : 0118 lnA A0-lnt 0119 t 为初至时间一, A 为第三波场的绝对能量, A0为常数, 为振幅补偿因子。 0120 以上所述的具体实施方式, 对本发明的目的、 技术方案和有益效果进行了进一步 详细说明, 所应理解的是, 以上所述仅为本发明的具体实施方式而已, 并不用于限定本发明 的保护范围, 凡在本发明的精神和原则之内, 所做的任何修改、 等同替换、 改进等, 均应包含 在本发明的保护范围之内。 说 明 书 CN 104502965 A 11 1/1 页 12 图 1 说 明 书 附 图 CN 104502965 A 12 。