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,它是理论的尾波振幅谱AC(t,f)和实测的尾波振幅谱A(t,f)的函数,定义为:2
△E(t,f,QI,QS)=|AC(t,f)-A(t,f)| (2)其中QI为固有衰减品质因子,QS为散射衰减品质因子;理论的尾波振幅谱AC(t,f)可由三维非均匀模型的尾波幅度表达式计算得到:其中,td为直达波走时,ID是直达波振幅AD在时间窗(t1,t2)的积分,表示为:步骤4:通过最小二乘获得每个深度点的散射衰减;
步骤5:循环执行步骤2至步骤4,最终获得深度区间散射衰减曲线。
2.根据权利要求1所述的一种多参数联合反演计算岩石散射衰减的方法,其特征在于:所述步骤4具体为:
根据步骤3计算每个接收器满足最小二乘求解目标函数△E的QI和QS的全局极小值,将每个深度点的N个接收器计算得到的散射衰减品质因子QS的均值 作为该深度点的散射衰减品质因子,取倒数得到该深度点的散射衰减值Atten: