基于偶函数地震响应以估算薄层厚度的方法及装置

摘要

本发明提供一种基于偶函数地震响应以估算薄层厚度的方法，包括：对零相位的奇函数地震响应进行90度相位转换以形成偶函数地震响应；对偶函数地震响应进行匹配追踪时频分析，以得到峰值振幅与峰值频率；将峰值振幅与峰值频率联合估算薄层厚度。实施上述方法的装置，包括：预处理模块用于对零相位的奇函数地震响应进行90度相位转换以变为偶函数地震响应；匹配追踪时频分析模块用于对偶函数地震响应进行分析，以得到峰值振幅与峰值频率；薄层厚度获取模块用于将峰值振幅与峰值频率联合估算薄层厚度。本发明验证了90度相位转换与薄层奇偶关系转换的等效性，确定偶函数对薄层厚度估算的优势，通过峰值振幅与峰值频率联合，有利于对薄层厚度进行预测。

主权项

一种基于偶函数地震响应以估算薄层厚度的方法，包括以下步骤：S100、对确定地震资料中符合顶底反射界面反射系数比为‑1的薄层位置，对薄层位置的子波相位进行估算并判断，若薄层位置的子波相位为零相位，对零相位的奇函数地震响应进行90度相位转换以形成偶函数地震响应；S200、利用下式对偶函数地震响应进行匹配追踪时频分析，以与其相对应的得到其峰值振幅和峰值频率：$u\left(t,f\right)=\underset{j=1}{\overset{J}{\Sigma }}{a}_j{W}_R^{\prime }\left(f,f_j\right)g\left(t\right)env\left[w_R^{\prime }\left(t-t_j,f_j\right)\right]\exp \left\{i\left[2\pi f\left(t-t_j\right)+{\phi }_j\right]\right\}$$g\left(t\right)=\left\{\begin{array}{c}\frac{1}{T}{e}^{-{\mathrm{at}}^2},\left|t\right|\le T\\ 0,\left|t\right|>T\end{array}\right.$其中，w_R'为匹配原子时间域表式，W_R'为匹配原子频率域表示，g(t)是高斯窗，控制能量干涉，env表示取匹配原子的时间包络，J是参与匹配的原子个数，u为偶函数地震响应的时频关系；j为匹配原子的第几个；a为匹配原子的振幅；f为频率采样点；t为时间采样点；T为高斯窗的能量持续时间；t_j为匹配原子的延时参数、f_j为匹配原子的主频参数、a_j为匹配原子的主振幅参数、φ_j为匹配原子的相位参数；S300、将峰值振幅与峰值频率联合，利用下式以估算出薄层的厚度：H＝a(As‑Asmin)/(Asmax‑Asmin)+b(fs‑fsmin)/(fsmax‑fsmin)+c；其中，H为当前采样点估计的薄层厚度；a为对应峰值振幅项系数；b为对应峰值频率项的系数；c为常数；As为当前采样点的峰值振幅；Asmin为整个薄层统计的最小峰值振幅；Asmax为整个薄层统计的最大峰值振幅；fs为当前采样点的峰值频率；fsmin为整个薄层统计的最小峰值频率；fsmax为整个薄层统计的最大峰值频率。