一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法及系统技术方案

本发明专利技术提供一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法及系统，所述系统在应力条件下计算随钻多极子声波测井地层扰动函数

A method and system for rapidly calculating dispersion curves of multipole acoustic logging while drilling under ground stress condition

The invention provides a method and system for fast calculating the dispersion curve of multipole acoustic logging while drilling under the condition of in-situ stress. The system calculates the formation disturbance function of multipole acoustic logging while drilling under the condition of stress.

【技术实现步骤摘要】
一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法及系统
本专利技术涉及随钻声波测井
，特别涉及一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法及系统。
技术介绍
随钻声波测井已成为深海、深地探测的必备技术之一，其目的是在钻井的过程中确定地层的信息，特别是纵横波速度，为油气田的勘探和开发提供重要的信息。井周异常地应力的存在会给随钻测井造成极大的安全隐患，快速计算不同地应力条件下的随钻多极子声波测井频散曲线对于异常地应力识别、反演地应力参数等至关重要。目前正演模拟地应力条件下多极子声波测井频散曲线响应特征的方法主要有两种：一种是利用有限差分或者有限元等数值解法计算地应力条件下多极子声波测井的频散曲线响应特征，该方法可用于电缆测井和随钻多极子声波测井的频散曲线数值模拟，但由于采用数值解法，计算量巨大，不能满足随钻声波测井的实时计算的需要；另一种是基于微扰理论的摄动积分方法，该方法是一种半数值半解析的模拟方法，可实现地应力条件下多极子声波测井频散曲线的快速计算，目前该方法已成功应用于地应力条件下电缆测井频散曲线快速数值模拟，但对于随钻多极子声波测井的快速数值模拟尚无成功应用的报道。
技术实现思路
针对现有技术的不足，本专利技术提供一种利用摄动积分快速计算地应力条件下随钻多极子声波测井频散曲线的方法，由于摄动积分是一种半解析版数值的方法及系统，它可以实现异常地应力随钻多极子声波测井频散曲线的快速数值模拟。本专利技术公开一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法，所述方法包括以下步骤：S1、根据均匀地层弹性参数及井径、随钻仪器弹性参数、内外径参数，计算角频率ω对应的参考状态以及随钻多极子频散曲线相速度v；S2、利用均匀地层弹性参数计算位移梯度和应变分量，给定异常地应力地层三阶弹性常数和地应力大小，计算随钻多极子声波测井的扰动函数S3、根据摄动积分计算地层地应力条件下角频率的改变量Δω；S4、根据计算随钻多极子声波测井频散曲线相速度改变量；S5、根据步骤S1和步骤S4得到的参考状态频散曲线相速度v和应力条件下频散曲线相速度改变量Δv，计算异常地应力地层频散曲线相速度Vstress；S6、重复步骤S1-S5，使角频率ω遍历所求频率区间，获取频率区间范围内应力地层的随钻多极子声波测井频散曲线。在上述技术方案中，所述步骤S1中弹性参数包括纵波速度、横波速度以及密度，所述的参考状态是均匀地层无异常地应力时的状态。在上述技术方案中，所述步骤S2中扰动函数为：其中，TLM是应力，cLγKM是二阶弹性常数，wν,K是位移梯度，EAB是应变分量，δγν是克罗尼克函数，是三阶弹性常数，且以上运算均满足爱因斯坦求和约定。在上述技术方案中，所述S3中摄动积分计算公式为：其中，ρ0是介质参考状态质量密度；Δρ表示密度的改变量；ω是本征频率；是表示某一模式波位移场的γ向分量，对柱坐标γ可取径向r，环向θ，和轴向z三个方向，下角标“，”表示对其后面的字母代表的变量求偏导。在上述技术方案中，所述步骤S4中计算得出的积分表达式如下：其中，Δv是应力条件下频散曲线相速度改变量，v是S1中计算得到的参考状态频散曲线相速度，r0＝0，r4＝∞，其中忽略了流体非线性效应，分子上包括上标Steel表示仪器和上标Formation表示地层两部分的积分；所述仪器主体由钢材料做成，相比于地层的孔隙岩石，钢的三阶弹性常数很小，因而公式3分子中的仪器积分部分若进一步省略，公式可以进一步的化简为：在上述技术方案中，所述步骤S5中异常地应力地层频散曲线相速度Vstress为Vstress＝v+Δv。本专利技术还公开一种快速计算地应力条件下随钻多极子声波测井频散曲线的系统，所述系统包括参考状态计算模块、扰动函数计算模块、角频率改变量计算模块、相速度改变量计算模块、频散曲线相速度计算模块、应力状态频散曲线绘制模块，其中：参考状态计算模块，根据均匀地层弹性参数及井径、随钻仪器弹性参数、内外径参数，计算角频率ω对应的参考状态以及随钻多极子频散曲线相速度v；扰动函数计算模块，利用均匀地层弹性参数计算位移梯度和应变分量，给定异常地应力地层三阶弹性常数和地应力大小，计算随钻多极子声波测井的扰动函数角频率改变量计算模块，利用摄动积分计算地层地应力条件下角频率的改变量Δω；相速度改变量计算模块，根据计算随钻多极子声波测井频散曲线相速度改变量；频散曲线相速度计算模块，根据得到的参考状态频散曲线相速度v和应力条件下频散曲线相速度改变量Δv，最后计算异常地应力地层频散曲线相速度Vstress；应力状态频散曲线绘制模块，使角频率ω遍历所求频率区间，获取频率区间范围内应力地层的随钻多极子声波测井频散曲线。在上述技术方案中，所述参考状态计算模块中，弹性参数包括纵波速度、横波速度以及密度，所述的参考状态是均匀地层无异常地应力时的状态。本专利技术一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法及系统，具有以下有益效果：本专利技术采用基于微扰理论的摄动积分实现异常地应力条件下随钻多极子声波测井频散曲线的快速数值模拟，可实现两方面的效益：其一是与已有的技术相比，该方法计算速度快，可满足随钻声波测井实时计算的需要；其二是为异常地应力条件下随钻多极子声波测井频散曲线响应特征快速计算为异常地应力识别及反演提供了可靠的正演模型。附图说明图1为本专利技术一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法流程图；图2本专利技术的利用摄动积分计算的应力条件下随钻四极子频散曲线图图3为本专利技术一种快速计算地应力条件下随钻多极子声波测井频散曲线的系统模块图；具体实施方式下面结合附图对本专利技术作进一步详细描述本专利技术公开一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法，所述方法包括以下步骤，如图1所示：S1、根据均匀地层弹性参数及井径、随钻仪器弹性参数、内外径参数，计算角频率ω对应的参考状态以及随钻多极子频散曲线相速度v；其中，所述弹性参数包括纵波速度、横波速度以及密度，所述的参考状态是均匀地层无异常地应力时的状态。S2、利用均匀地层弹性参数计算位移梯度和应变分量，给定异常地应力地层三阶弹性常数和地应力大小，计算随钻多极子声波测井的扰动函数其中所述扰动函数为：其中，TLM是应力，cLγKM是二阶弹性常数，wν,K是位移梯度，EAB是应变分量，δγν是克罗尼克函数，是三阶弹性常数，且以上运算均满足爱因斯坦求和约定。S3、根据摄动积分计算地层地应力条件下角频率的改变量Δω；其中，所述S3中摄动积分计算公式为：其中，ρ0是介质参考状态质量密度；Δρ表示密度的改变量；ω是本征频率；是表示某一模式波位移场的γ向分量，对柱坐标γ可取径向r，环向θ，和轴向z三个方向，下角标“，”表示对其后面的字母代表的变量求偏导。S4、根据计算随钻多极子声波测井频散曲线相速度改变量；其中计算得出的积分表达式如下：其中，Δv是应力条件下频散曲线相速度改变量，v是S1中计算得到的参考状态频散曲线相速度，r0＝0，r4＝∞，其中忽略了流体非线性效应，分子上包括上标Steel表示仪器和上标Formation表示地层两部分的积分；所述仪器主体由钢材料做成，相比于地层的孔隙岩石，钢的三阶弹性常数很小，其中所述仪器优选随钻仪器，因而公式3分本文档来自技高网...

【技术保护点】
1.一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法，其特征在于，所述方法包括以下步骤：S1、根据均匀地层弹性参数及井径、随钻仪器弹性参数、内外径参数，计算角频率ω对应的参考状态以及随钻多极子频散曲线相速度v；S2、利用均匀地层弹性参数计算位移梯度和应变分量，给定异常地应力地层三阶弹性常数和地应力大小，计算随钻多极子声波测井的扰动函数

【技术特征摘要】
1.一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法，其特征在于，所述方法包括以下步骤：S1、根据均匀地层弹性参数及井径、随钻仪器弹性参数、内外径参数，计算角频率ω对应的参考状态以及随钻多极子频散曲线相速度v；S2、利用均匀地层弹性参数计算位移梯度和应变分量，给定异常地应力地层三阶弹性常数和地应力大小，计算随钻多极子声波测井的扰动函数S3、根据摄动积分计算地层地应力条件下角频率的改变量Δω；S4、根据计算随钻多极子声波测井频散曲线相速度改变量；S5、根据步骤S1和步骤S4得到的参考状态频散曲线相速度v和应力条件下频散曲线相速度改变量Δv，计算异常地应力地层频散曲线相速度Vstress；S6、重复步骤S1-S5，使角频率ω遍历所求频率区间，获取频率区间范围内应力地层的随钻多极子声波测井频散曲线。2.根据权利要求1所述一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法，其特征在于，所述步骤S1中弹性参数包括纵波速度、横波速度以及密度，所述的参考状态是均匀地层无异常地应力时的状态。3.根据权利要求1所述一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法，其特征在于，所述步骤S2中扰动函数为：其中，TLM是应力，cLγKM是二阶弹性常数，wν,K是位移梯度，EAB是应变分量，δγν是克罗尼克函数，是三阶弹性常数，且以上运算均满足爱因斯坦求和约定。4.根据权利要求1所述一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法，其特征在于，所述步骤S3中摄动积分计算公式为：其中，ρ0是介质参考状态质量密度；Δρ表示密度的改变量；ω是本征频率；是表示某一模式波位移场的γ向分量，对柱坐标γ可取径向r，环向θ，和轴向z三个方向，下角标“，”表示对其后面的字母代表的变量求偏导。5.根据权利要求1所述一种快速计算地应力条件下随钻多极子声波测井频散曲线的方法，其特征在于，所述步骤S4中计算得出的积分表达式...

【专利技术属性】
技术研发人员：黑创，苏远大，唐晓明，庄春喜，李修权，罗明璋，
申请(专利权)人：长江大学，
类型：发明
国别省市：湖北,42