一种获取地层孔隙度的方法及装置的制造方法
【技术领域】
[0001] 本发明属于石油勘探资料处理领域，特别涉及获取地层孔隙度的方法及装置。
【背景技术】
[0002] 地层孔隙度是储层评价的重要参数之一，地层孔隙度与测井电位、钻井时间、油气 储量以及含油饱和度等石油勘探参数紧密相连，因此，获取地层孔隙度是石油勘探的一个 重要工作步骤。现有技术中，通常根据电缆声波测得的时差计算得到所述地层孔隙度。
[0003] 在实现本申请过程中，发明人发现现有技术中至少存在如下问题：在很多地层中， 尤其在大斜度井以及水平井中，存在时差各向异性现象。这时电缆声波测得的时差是水平 时差和垂直时差的综合值。而各向异性现象多数是由于地层内夹杂着具有不同时差的其他 地层而产生的。计算所述孔隙度所需要的时差一般是水平时差，因此使用电缆声波测得的 时差计算得到的地层孔隙度偏离真实地层孔隙度，进而根据所述地层孔隙度计算得到的含 油饱和度等参数值误差也相对较大。
【发明内容】
[0004] 针对现有技术中根据声波测井仪器测得的时差计算得到的孔隙度误差较大，本发 明提出一种获取地层孔隙度的方法及装置，所述获取地层孔隙度的方法及装置具体是这样 实现的：
[0005] 一种获取地层孔隙度的方法，所述方法包括：
[0006] 获取目标地层的声波测井资料、井斜角曲线资料以及层位时差各向异性比；
[0007] 根据所述目标地层的井斜角曲线资料以及层位时差各向异性比，按照拟合方程计 算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差，所述拟合方程 根据已知的实验地层训练得到；
[0008] 根据所述目标地层声波测井资料中的时差以及所述目标地层的慢度差计算所述 目标地层的水平时差；
[0009] 根据所述目标地层的水平时差计算所述目标地层的孔隙度。
[0010] 可选的，在本发明一实施例中，所述根据已知的实验地层训练得到拟合方程，包 括：
[0011] 获取所述实验地层的水平时差，计算所述实验地层在不同井斜角、不同时差各向 异性比条件下的时差；
[0012] 根据所述实验地层的水平时差，计算所述实验地层的慢度差，所述慢度差为所述 时差与所述实验地层的水平时差的差值；
[0013] 建立所述井斜角以及所述慢度差的坐标系，在所述坐标系中插入所述不同井斜 角、不同时差各向异性比条件下的所述慢度差的离散点；
[0014] 将所述离散点拟合成曲线，并获取所述曲线的拟合方程。
[0015] 可选的，在本发明一实施例中，所述计算所述实验地层在不同井斜角、不同时差各 向异性比条件下的时差，包括：
[0016] 利用声波测井仪采用有限差分法计算所述实验地层在不同井斜角、不同时差各向 异性比条件下的时差，其中，所述井斜角从0到90°范围内按15°步长递增，所述时差各向 异性比从I. 0到1. 2范围内按0. 02步长递增。
[0017] 可选的，在本发明一实施例中，所述声波测井仪包括双发双收声波测井仪。
[0018] 可选的，在本发明一实施例中，所述建立所述井斜角以及所述慢度差的坐标系，具 体包括：
[0019] 建立所述井斜角正弦值以及所述慢度差的直角坐标系。
[0020] 可选的，在本发明一实施例中，所述拟合方程包括一元二次方程。
[0021] 可选的，在本发明一实施例中，所述按照所述拟合方程，计算所述目标地层在不同 井斜角、不同层位时差各向异性比条件下的慢度差，具体包括：
[0022] 根据所述实验地层井斜角以及慢度差的坐标系中的曲线，按照所述拟合方程，计 算与所述实验地层相同的时差各向异性比条件下，所述目标地层在不同井斜角条件下的慢 度差，当所述井斜角是步长15°的倍数时，按照所述拟合方程求解所述目标地层的慢度差； 当所述井斜角不是15°的倍数时，采用插值法计算所述目标地层的慢度差。
[0023] 可选的，在本发明一实施例中，根据所述目标地层的水平时差计算所述目标地层 的孔隙度，其中所述孔隙度的计算公式为：
[0025] 式中，为孔隙度，Δ tM、Δ tf、Δ tsh分别为岩石骨架声波水平方向时差、地层流 体声波水平方向时差、地层流体声波水平方向时差，V sh为地层泥质含量，CpS声波压实校正 系数，At为所述目的地层的水平时差。
[0026] -种获取地层孔隙度的装置，包括：
[0027] 资料获取模块，用于获取目标地层的声波测井资料、井斜角曲线资料以及层位时 差各向异性比；
[0028] 目标慢度差计算模块，用于根据所述目标地层的井斜角曲线资料以及层位时差各 向异性比，按照拟合方程计算所述目标地层在不同井斜角、不同层位时差各向异性比条件 下的慢度差，所述拟合方程根据拟合装置计算得到；
[0029] 水平时差计算模块，用于根据所述目标地层声波测井资料中的时差以及所述目标 地层的慢度差计算所述目标地层的水平时差；
[0030] 孔隙度计算模块，用于根据所述目标地层的水平时差计算所述目标地层的孔隙 度。
[0031] 可选的，在本发明一实施例中，所述目标慢度差计算模块还包括拟合装置，所述拟 合装置包括：
[0032] 实验时差计算模块，用于获取实验地层的水平时差，计算所述实验地层在不同井 斜角、不同时差各向异性比条件下的时差；
[0033] 慢度差计算模块，用于根据所述实验地层的水平时差，计算所述实验地层的慢度 差，所述慢度差为所述时差与所述实验地层的水平时差的差值；
[0034] 坐标系建立模块，用于建立所述井斜角以及所述慢度差的坐标系，在所述坐标系 中插入所述不同井斜角、不同时差各向异性比条件下的所述慢度差的离散点；
[0035] 拟合模块，用于将所述离散点拟合成曲线，并获取所述曲线的拟合方程。
[0036] 可选的，在本发明一实施例中，所述实验时差计算模块，还包括：
[0037] 时差计算子模块，用于利用声波测井仪采用有限差分法计算所述实验地层在不同 井斜角、不同时差各向异性比条件下的时差，其中，所述井斜角从〇到90°范围内按15°步 长递增，所述时差各向异性比从1. 〇到1. 2范围内按0. 02步长递增。
[0038] 可选的，在本发明一实施例中，所述声波测井仪包括双发双收声波测井仪。
[0039] 可选的，在本发明一实施例中，所述坐标系建立模块，还包括：
[0040] 坐标系建立子模块，用于建立所述井斜角正弦值以及所述慢度差的直角坐标系。
[0041] 可选的，在本发明一实施例中，所述拟合方程包括一元二次方程。
[0042] 可选的，在本发明一实施例中，所述目标慢度差计算模块，还包括：
[0043] 目标慢度差计算子模块，用于根据所述实验地层井斜角以及慢度差的坐标系中的 曲线，按照所述拟合方程，计算与所述实验地层相同的时差各向异性比条件下，所述目标地 层在不同井斜角条件下的慢度差，当所述井斜角是步长15°的倍数时，按照所述拟合方程 求解所述目标地层的慢度差；当所述井斜角不是15°的倍数时，采用插值法计算所述目标 地层的慢度差。
[0044] 可选的，在本发明一实施例中，所述孔隙度