[发明专利]一种多参数联合反演计算岩石散射衰减的方法

摘要

本发明公开了一种多参数联合反演计算岩石散射衰减的方法，首先在深度区间内进行阵列声波测井，得到深度区间内偶极阵列波形数据；计算每个深度点实测波形的功率谱；构建反演目标函数；通过最小二乘获得每个深度点的散射衰减；循环执行步骤2至步骤4，最终获得深度区间散射衰减曲线。从声波测井数据正确分析和提取弹性波的散射衰减，将会提供很有价值的地层信息，与大尺度地震波散射衰减相比，声波测井计算的散射衰减可提供地下岩石小尺度的非均匀性信息。

主权项

1.一种多参数联合反演计算岩石散射衰减的方法，其特征在于：该方法包括以下步骤：步骤1：在深度区间内进行阵列声波测井，得到深度区间内偶极阵列波形数据；步骤2：计算每个深度点实测波形的功率谱；所述步骤2具体为：对步骤1得到的偶极阵列波形数据进行滤波，获得处理深度位置处的偶极子全波列数据v(t)，设定窗函数g(t)，利用短时傅里叶变换(STFT)计算出实测波形不同时刻的功率谱A(t,f)：其中，t为时间，ω为角频率，f为频率，A(t,f)函数给出了以t时为中心的时窗g(t)内信号v(t)的傅里叶振幅谱；步骤3：构建反演目标函数；所述步骤3具体为：构建目标函数△E，它是理论的尾波振幅谱A_C(t,f)和实测的尾波振幅谱A(t,f)的函数，定义为：△E(t,f,Q_I,Q_S)＝|A_C(t,f)‑A(t,f)|²             (2)其中Q_I为固有衰减品质因子，Q_S为散射衰减品质因子；理论的尾波振幅谱A_C(t,f)可由三维非均匀模型的尾波幅度表达式计算得到：其中，t_d为直达波走时，I_D是直达波振幅A_D在时间窗(t₁,t₂)的积分，表示为：步骤4：通过最小二乘获得每个深度点的散射衰减；步骤5：循环执行步骤2至步骤4，最终获得深度区间散射衰减曲线。