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

权利要求书

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，最终获得深度区间散射衰减曲线。

2.根据权利要求1所述的一种多参数联合反演计算岩石散射衰减的方法，其特征在于：所述步骤4具体为：

根据步骤3计算每个接收器满足最小二乘求解目标函数△E的Q_I和Q_S的全局极小值，将每个深度点的N个接收器计算得到的散射衰减品质因子Q_S的均值作为该深度点的散射衰减品质因子，取倒数得到该深度点的散射衰减值Atten：