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

摘要

本发明公开了一种正演模拟蚀变地层随钻多极子声波频散曲线的方法，用于正演模拟随钻声波的频散曲线。本发明提供的方法包括：获取在参考状态下预设角频率及对应的相速度；根据公式计算地层弹性参数扰动函数；通过摄动积分函数得到扰动后角频率的改变量；计算地层弹性参数扰动后随钻多极子频散曲线相速度改变量；进而得到蚀变地层频散曲线相速度；分别计算在频率区间内每个角频率对应的相速度，正演模拟得到所述频率区间内蚀变地层随钻多极子声波测井频散曲线。在本发明中，通过摄动积分计算地层参数扰动后的影响，克服了现有正演频散曲线无法模拟井外连续蚀变的不足，并大大提升连续蚀变模型的频散曲线的计算速度。

主权项

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)，计算得到蚀变地层频散曲线相速度Vr；Vr＝v+Δv  (4)其中，ν为参考状态下频散曲线的相速度，Δv表示扰动后随钻多极子频散曲线相速度改变量Δv；步骤S6：分别计算在频率区间内每个角频率ω对应的相速度Vr，正演模拟得到所述频率区间内蚀变地层随钻多极子声波测井频散曲线。