[发明专利]一种基于地震信号包络峰值的地层介质品质因子估计方法

摘要

本发明涉及一种基于地震信号包络峰值的地层介质品质因子估计方法，其步骤：1）将两相邻检波器的距离为厚度把介质分成若干个小薄板；2）对于第i个小薄板，用常相位子波逼近该小薄板顶部接收的直达波，用基于MSMG的高阶累积量匹配法估计地震子波参数；3）用基于MSMG的高阶累积量匹配法估计小波变换的核函数的参数；4）计算第i个小薄板上下界面处接收信号的瞬时频率；5）计算第i个小薄板内包络峰值瞬时频率的变化；6）计算第i个小薄板的品质因子Q值；7）重复步骤2）～6），依次计算除第i个小薄板外的其它N-1个小薄板的品质因子；8）利用步骤6）和7）得到的N个小薄板的品质因子Q值来预测储层的含油气性，进而得到估计的衰减曲线。

主权项

1.一种基于地震信号包络峰值的地层介质品质因子估计方法，其包括以下步骤：1）以井中两相邻检波器的距离为厚度把介质分成若干个小薄板；2）对于第i个小薄板，用常相位子波逼近该小薄板顶部接收的直达波，用基于MSMG的高阶累积量匹配法估计地震子波参数；其中MSMG为多变量高斯混合模型；用常相位子波逼近小薄板顶部接收直达波的方法为：（1）在水平层状黏弹介质中，设各层的品质因子Q值为常数，仅考虑平面波的单程波传播，则位于地表处的震源子波传到深度z处时其频率表达式为：U^(Δz,ω)=G(Δz)U^(0,ω)exp[-iωΔzc(ω)-ωΔz2c(ω)Q],]]>其中，Δz为震源子波的传输距离，ω为震源子波的角频率，G(Δz)为独立于频率和吸收的因子，为震源子波的频域表达式，exp为以自然对数e为底的指数函数，c(ω)为震源子波的相速度；（2）采用常相位子波去逼近震源子波，该常相位子波的函数表达式为：其中，为震源子波的相位常数，σ为震源子波的尺度参数，控制着子波的宽度，δ为地震子波的频率参数，t为震源子波的传播时间；用基于MSMG的高阶累积量匹配法估计地震子波参数的过程为：（1）利用已有相位旋转方法估计地震子波的相位参数（2）利用累积量匹配方法估计尺度参数σ和频率参数δ；3）用基于MSMG的高阶累积量匹配法估计小波变换的核函数的参数；4）计算第i个小薄板上下界面处接收信号的瞬时频率；5）计算第i个小薄板内包络峰值瞬时频率的变化；6）计算第i个小薄板的品质因子Q值：将步骤3）得到的匹配地震小波的三个参数σ_i,δ_i和τ_i的值与步骤5）得到的第i个小薄板内包络峰值瞬时频率的变化Δf_p(i)代入包络峰值瞬时频率和品质因子Q值之间的解析关系式，得到品质因子Q值的估计值。所述包络峰值瞬时频率和品质因子Q值之间的解析关系式为：Q≈τδ2κ(η)4πΔfp,]]>其中，κ(η)=1-2πηΦ-1(2πη)exp(-2π2η2),]]>Φ(x)=12π∫-∞xexp(-t22)dt]]>为标准正态分布的概率积分函数，t为积分变量，7）重复步骤2）～6），依次计算除第i个小薄板外的其它N-1个小薄板的品质因子；8）利用步骤6）和7）得到的N个小薄板的品质因子Q值来预测储层的含油气性，进而得到估计的衰减曲线。