[发明专利]一种正演模拟蚀变地层随钻多极子声波频散曲线的方法

权利要求书

1.一种正演模拟随钻多极子声波频散曲线的方法，其特征在于，包括：

步骤S1：获取在参考状态下预设角频率ω对应的随钻多极子频散曲线的相速度ν；

步骤S2：设定弹性波的声速随径向变化的函数，根据公式(1)计算地层弹性参数扰动函数

其中，Δλ(r)和Δμ(r)表示参考状态拉梅系数λ和μ的改变量，随径向距离r变化的函数，δ_ij(i,j＝L,γ,M,υ)是克罗尼克函数；

步骤S3：通过摄动积分函数(2)计算扰动后角频率的改变量△ω；

其中，ρ₀表示介质参考状态质量密度，△ρ表示密度的改变量，ω是本征频率，是表示参考状态下波(序号为m)位移场的γ向分量，对柱坐标γ可取径向r，环向θ，和轴向z三个方向，下角标“，”表示对其后面的字母代表的变量求偏导，上标“*”表示共轭，N表示计算模型中介质的层数，r₀＝0，r₁为随钻仪器内径，r₂为随钻仪器外径，r₃为井径，r_N＝∞；

步骤S4：根据公式(3)，计算地层弹性参数扰动后随钻多极子频散曲线相速度改变量Δv；

其中，ω为预设角速度，△ω表示扰动后角频率的改变量，ν散曲线的相速度；

步骤S5：根据公式(4)，计算得到蚀变地层频散曲线相速度V_r；

V_r＝v+Δv (4)

其中，ν为参考状态下频散曲线的相速度，Δv表示扰动后随钻多极子频散曲线相速度改变量Δv；

步骤S6：分别计算在频率区间内每个角频率ω对应的相速度V_r，正演模拟得到所述频率区间内蚀变地层随钻多极子声波测井频散曲线。

2.根据权利要求1所述的方法，其特征在于，所述步骤S1具体为：

获取参考状态下的地层弹性参数，以及井径、随钻仪器弹性、内外径，根据公式(1)的频散方程，计算角频率ω对应的随钻多极子频散曲线相速度ν；

D(k,ω)＝|H_ij|_12×12＝0 (5)

其中，k是波数，ω是角频率，H_ij为12×12矩阵元素。

3.根据权利要求2所述的方法，其特征在于，所述步骤S2具体为：

设定弹性波参数随径向变化的函数，计算得到拉梅系数随径向变化变化函数的函数Δλ(r)和Δμ(r)，根据公式(1)计算地层弹性参数扰动函数

4.根据权利要求1所述的方法，其特征在于，所述步骤S6中所述频率区间为不低于截止频率任意设置的频率区间。