CN108897043B - 基于弱度参数检测裂缝型储层流体的方法及装置 - Google Patents

基于弱度参数检测裂缝型储层流体的方法及装置 Download PDF

Info

Publication number
CN108897043B
CN108897043B CN201811020121.7A CN201811020121A CN108897043B CN 108897043 B CN108897043 B CN 108897043B CN 201811020121 A CN201811020121 A CN 201811020121A CN 108897043 B CN108897043 B CN 108897043B
Authority
CN
China
Prior art keywords
fluid factor
fluid
parameter
reservoir
weakness parameter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811020121.7A
Other languages
English (en)
Other versions
CN108897043A (zh
Inventor
谢春辉
杨午阳
王恩利
杜炳毅
闫国亮
李海山
杨庆
鄢高韩
何润
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201811020121.7A priority Critical patent/CN108897043B/zh
Publication of CN108897043A publication Critical patent/CN108897043A/zh
Application granted granted Critical
Publication of CN108897043B publication Critical patent/CN108897043B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters

Abstract

本发明提供了一种基于弱度参数检测裂缝型储层流体的方法及装置，涉及石油勘探技术领域。方法包括：获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；根据入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；根据各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；根据裂缝型储层背景介质的拉梅常数对初始流体因子进行修正，获得第二级流体因子；根据裂缝型储层的裂缝纵横比对第二级流体因子进行修正，获得流体因子结果；将流体因子结果与预先设置的阈值区间进行比较，以确定流体因子结果对应的储层流体类别。

Description

基于弱度参数检测裂缝型储层流体的方法及装置
技术领域
本发明涉及石油勘探技术领域，尤其涉及一种基于弱度参数检测裂缝型储层流体的方法及装置。
背景技术
目前，裂缝型油气藏的产量占全世界石油天然气总产量的一半以上，在国内的油气勘探中，裂缝型油气藏是一个重要的勘探领域，如塔里木盆地奥陶系碳酸盐岩裂缝型油气藏，四川致密砂岩裂缝型油气藏。裂缝既是油、气、水等地下流体的运移通道，也是流体的储集空间，裂缝能极大地提升储层的渗透率。裂缝型储层描述包括预测裂缝分布特征和检测裂隙充填物，也就是探测开启并充填流体的有效裂缝。
利用地震资料研究岩石裂隙所贮存流体的类型和特性，是根据岩石物理等效介质理论，将与储层流体有关的异常特性表征为流体因子，再依托流体因子实现储层流体的类型判识。对于各向同性介质，现有技术中通常采用弹性参数的代数组合作为流体因子，常见的有LMR(Lumbda-Mu-Rho)法、Russell法、纵横波阻抗的幂次方等。流体检测的质量主要取决于两方面，一是反演计算的弹性参数是否可靠；二是构建的流体因子对流体类型是否敏感且稳定。可见，对于裂缝型储层，应该用不随方位变化的各向异性参数来构建流体因子。目前，如何能够准确可靠的检测裂缝型储层流体的类型成为了一个亟待解决的问题。
发明内容
本发明的实施例提供一种基于弱度参数检测裂缝型储层流体的方法及装置，以实现准确可靠的检测裂缝型储层流体的类型。
为达到上述目的，本发明采用如下技术方案：
一种基于弱度参数检测裂缝型储层流体的方法，包括：
获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；
根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；
根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；
将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。
具体的，所述根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子，包括：
获取裂缝介质的等效弹性系数矩阵：
Figure BDA0001787079660000021
其中，
Figure BDA0001787079660000022
表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵；λ和μ为各向同性背景介质拉梅常数；ΔT用于表示切向弱度参数；ΔN用于表示法向弱度参数；
Figure BDA0001787079660000023
将裂缝介质的等效弹性系数矩阵
Figure BDA0001787079660000024
中的元素代入到公式：
Figure BDA0001787079660000025
计算切向弱度参数ΔT和法向弱度参数ΔN；其中，
Figure BDA0001787079660000026
其中，cij表示裂缝介质的等效弹性系数矩阵
Figure BDA0001787079660000031
中的第i行第j列元素；εV、δV和γ表示各向异性参数；
将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子
Figure BDA0001787079660000032
其中，
Figure BDA0001787079660000033
具体的，所述根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子，包括：
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf；其中，0≈μf<<λf；
根据公式：
Figure BDA0001787079660000034
对所述初始流体因子
Figure BDA0001787079660000035
进行修正，获得第二级流体因子
Figure BDA0001787079660000036
其中，
Figure BDA0001787079660000037
g为采用背景介质的拉梅常数λ和μ计算得到的参数，
Figure BDA0001787079660000038
d表示裂缝型储层的裂缝纵横比。
具体的，所述根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果，包括：
根据公式：
Figure BDA0001787079660000039
以及裂缝型储层的裂缝纵横比d，对所述第二级流体因子
Figure BDA00017870796600000310
进行修正，获得流体因子结果Ifluid；
其中，
Figure BDA00017870796600000311
具体的，所述预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间；
所述将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别，包括：
将所述流体因子结果与所述流体因子阈值区间进行比较，以确定所述流体因子结果所在的流体因子阈值区间，并确定该流体因子阈值区间对应的储层流体类别。
一种基于弱度参数检测裂缝型储层流体的装置，包括：
三维地震资料处理单元，用于获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；
各向异性参数求取单元，用于根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；
初始流体因子计算单元，用于根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；
第二级流体因子获得单元，用于根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；
流体因子结果获得单元，用于根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；
储层流体类别确定单元，用于将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。
另外，所述初始流体因子计算单元，具体用于：
获取裂缝介质的等效弹性系数矩阵：
Figure BDA0001787079660000041
其中，
Figure BDA0001787079660000042
表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵；λ和μ为各向同性背景介质拉梅常数；ΔT用于表示切向弱度参数；ΔN用于表示法向弱度参数；
Figure BDA0001787079660000043
将裂缝介质的等效弹性系数矩阵
Figure BDA0001787079660000044
中的元素代入到公式：
Figure BDA0001787079660000045
计算切向弱度参数ΔT和法向弱度参数ΔN；其中，
Figure BDA0001787079660000046
其中，cij表示裂缝介质的等效弹性系数矩阵
Figure BDA0001787079660000051
中的第i行第j列元素；εV、δV和γ表示各向异性参数；
将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子
Figure BDA0001787079660000052
其中，
Figure BDA0001787079660000053
另外，所述第二级流体因子获得单元，具体用于：
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf；其中，0≈μf<<λf；
根据公式：
Figure BDA0001787079660000054
对所述初始流体因子
Figure BDA0001787079660000055
进行修正，获得第二级流体因子
Figure BDA0001787079660000056
其中，
Figure BDA0001787079660000057
g为采用背景介质的拉梅常数λ和μ计算得到的参数，
Figure BDA0001787079660000058
d表示裂缝型储层的裂缝纵横比。
另外，所述流体因子结果获得单元，具体用于：
根据公式：
Figure BDA0001787079660000059
以及裂缝型储层的裂缝纵横比d，对所述第二级流体因子
Figure BDA00017870796600000510
进行修正，获得流体因子结果Ifluid；
其中，
Figure BDA00017870796600000511
另外，所述储层流体类别确定单元中预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间；
所述储层流体类别确定单元，具体用于：
将所述流体因子结果与所述流体因子阈值区间进行比较，以确定所述流体因子结果所在的流体因子阈值区间，并确定该流体因子阈值区间对应的储层流体类别。
一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现以下步骤：
获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；
根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；
根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；
将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。
一种计算机设备，包括存储器、处理器及存储在存储上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现以下步骤：
获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；
根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；
根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；
将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。
本发明实施例提供的一种基于弱度参数检测裂缝型储层流体的方法及装置，首先获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。本发明可以实现准确可靠的检测裂缝型储层流体的类型的目的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施提供的一种基于弱度参数检测裂缝型储层流体的方法的流程图一；
图2为本发明实施提供的一种基于弱度参数检测裂缝型储层流体的方法的流程图二；
图3为裂缝型储层充填不同流体时的方位角道集特征示意图；
图4为裂缝型储层充填不同流体时，切向弱度参数ΔT和法向弱度参数ΔN随裂缝密度变化的关系图；
图5为初始流体因子随裂缝密度变化的关系图；
图6为流体因子的阈值区间划分图；
图7为目标储层的流体因子沿层切片示意图；
图8为本发明实施提供的一种基于弱度参数检测裂缝型储层流体的装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明实施提供一种基于弱度参数检测裂缝型储层流体的方法，包括：
步骤101、获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集。
步骤102、根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数。
步骤103、根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子。
步骤104、根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子。
步骤105、根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果。
步骤106、将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。
为了使本领域的技术人员更好的了解本发明，下面列举一个更为详细的实施例，如图2所示，本发明实施例提供一种基于弱度参数检测裂缝型储层流体的方法，包括：
步骤201、获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集。
例如，如图3所示，裂缝型储层充填不同流体时的方位角道集特征不同。
步骤202、根据所述入射角道集，采用振幅随方位角变化反演算法(AmplitudeVersus Azimuth，简称AVAZ)计算得到各向异性参数。
对于垂直排列的平行裂缝，可以将其等效为HTI介质(Horizontal TransverseIsotropy，即具有水平对称轴的横向各向同性介质)模型来研究。若入射角为θ，测线方位角为φ，则HTI介质分界面处的纵波反射系数是各向异性参数、入射角和方位角的函数：
Figure BDA0001787079660000081
上式中：φs为裂缝平面的对称轴方向，即裂缝面法向，(φ-φs)为测线与裂缝面法向的夹角；α,β,Z,G分别是纵波速度、横波速度、纵波阻抗和剪切模量；εV,δV,γ为各向异性参数。
Figure BDA0001787079660000082
分别表示上下两层介质间纵波速度的平均值、横波速度的平均值、纵波阻抗的平均值和剪切模量的平均值，Δα,ΔZ,ΔG,ΔεV,ΔδV,Δγ分别表示上下两层介质间纵波速度的差值、纵波阻抗的差值、剪切模量的差值以及三个各向异性参数的差值。
对上式进行简化，地震纵波的反射振幅R与入射角θ、方位角φ以及各向异性参数之间的关系式如下：
Figure BDA0001787079660000091
其中A，B，C，D，E，F是入射角θ和方位角φ的组合：
Figure BDA0001787079660000092
各向异性参数的反演方程组包括下述超定方程组：
Mn×6X6×1＝Rn×1
上式中n为覆盖次数，M，X和R分别为：
Figure BDA0001787079660000093
设定阻尼因子，然后利用奇异值分解法求解所述超定方程组，求解的结果为：
X＝M+R
其中M+为矩阵M的广义逆。
通过求解得到的X，可计算第i个地层分界面的各向异性参数差值：
Figure BDA0001787079660000109
Figure BDA00017870796600001010
Δγi＝X(6)
从上到下累加各向异性参数差值，可得第n个地层的各向异性参数：
Figure BDA0001787079660000101
Figure BDA0001787079660000102
Figure BDA0001787079660000103
对于裂缝型储层，参数εV,δV,γ描述了介质的各向异性程度。
步骤203、获取裂缝介质的等效弹性系数矩阵：
Figure BDA0001787079660000104
其中，
Figure BDA0001787079660000105
表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵；λ和μ为各向同性背景介质拉梅常数，可通过AVO(振幅随偏移距的变化，Amplitude variationwith offset)三参数反演方法计算而得到；ΔT用于表示切向弱度参数；ΔN用于表示法向弱度参数；
Figure BDA0001787079660000106
步骤204、将裂缝介质的等效弹性系数矩阵
Figure BDA0001787079660000107
中的元素代入到公式：
Figure BDA0001787079660000108
则有：
Figure BDA0001787079660000111
从而可计算切向弱度参数ΔT和法向弱度参数ΔN。其中，
Figure BDA0001787079660000112
Figure BDA0001787079660000113
其中，cij表示裂缝介质的等效弹性系数矩阵
Figure BDA0001787079660000114
中的第i行第j列元素；εV、δV和γ表示各向异性参数。
对于计算得到的切向弱度参数ΔT和法向弱度参数ΔN，如图4所示，裂缝型储层充填不同流体时，切向弱度参数ΔT和法向弱度参数ΔN随裂缝密度变化(crack density表示裂缝密度)。
步骤205、将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子
Figure BDA0001787079660000115
其中，
Figure BDA0001787079660000116
对于初始流体因子
Figure BDA0001787079660000117
如图5所示，初始流体因子随裂缝密度变化很小。
步骤206、获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf。
其中，0≈μf<<λf。
步骤207、根据公式：
Figure BDA0001787079660000118
对所述初始流体因子
Figure BDA0001787079660000119
进行修正，获得第二级流体因子
Figure BDA00017870796600001110
其中，
Figure BDA00017870796600001111
g为采用背景介质的拉梅常数λ和μ计算得到的参数，
Figure BDA00017870796600001112
d表示裂缝型储层的裂缝纵横比，可通过地下介质裂隙结构反演计算得到。
此处，对于初始流体因子
Figure BDA00017870796600001113
的公式，原理如下：
在地下储层含有裂缝时，其弹性系数矩阵C66可由Schoenberg模型或者Hudson模型求得。
Schoenberg模型也称为线性滑动模型，是指裂缝表面位移不连续，位移与穿过裂缝的应力呈线性关系的模型。设本构坐标系下裂隙的法向与x轴平行，裂隙密度为e，纵横比为d，各向同性背景的拉梅(Lame)常数为λ和μ，裂隙填充物的拉梅常数为λf和μf。引入切向弱度参数ΔT和法向弱度参数ΔN，裂缝介质的等效弹性系数矩阵如公式(1)：
Figure BDA0001787079660000121
其中：
Figure BDA0001787079660000122
弹性系数矩阵C66为6×6阶的矩阵，上标Schoenberg表示Schoenberg模型。
而Hudson模型中，假定介质包含比地震波长小得多的定向稀疏排列裂隙，裂隙是彼此分离的薄片，裂隙内包含的气体、流体等物质的体积模量比背景介质小。裂缝介质的等效弹性系数矩阵为公式(2)：
Figure BDA0001787079660000123
其中：
Figure BDA0001787079660000124
其中，U11和U33表达式如公式(2.b)：
Figure BDA0001787079660000131
其中M、k表达式如公式(2.c)，与充填流体的拉梅常数相关：
Figure BDA0001787079660000132
对应相同的裂缝介质，Hudson模型与Schoenberg模型所描述的弹性系数矩阵相等，即
Figure BDA0001787079660000133
由
Figure BDA0001787079660000134
和
Figure BDA0001787079660000135
得公式(3)：
Figure BDA0001787079660000136
化简后得公式(3.a)：
Figure BDA0001787079660000137
将公式(2.b)的U11，U33代入公式(3.a)，得到公式(3.b)：
Figure BDA0001787079660000138
而将切向弱度参数ΔT和法向弱度参数ΔN的比值记为初始流体因子
Figure BDA0001787079660000139
Figure BDA00017870796600001310
流体的第二拉梅常数μf接近于零值，相比第一拉梅常数λf小得多，即有：
0≈μf<<λf
依据公式(2.c)，有：
0≈M<<κ
此时，公式(4)可变换为公式(4.a)：
Figure BDA0001787079660000141
步骤208、根据公式：
Figure BDA0001787079660000142
以及裂缝型储层的裂缝纵横比d，对所述第二级流体因子
Figure BDA0001787079660000143
进行修正，获得流体因子结果Ifluid。
此处，由于：
0≈μf<<λf；
因此存在
Figure BDA0001787079660000144
对于流体因子结果中的流体因子，其阈值区间的划分如图6所示，其中，probability是指概率，hydrocarbon indicator是指烃类指示，说明了流体因子阈值区间是如何划分的。在井点处，有纵横波速度和密度、可以计算得到背景介质的拉梅常数(λ、μ)；也能解释得到裂缝密度e和裂缝纵横比d；从而计算获流体因子的值；把这些值(很多散点)生成直方图。当各种流体的直方图没有重叠时，能很轻易的划分阈值区间；当直方图存在部分重合时，则在允许的错误识别率之下，根据概率分布图划定阈值区间。
步骤209、将流体因子结果与流体因子阈值区间进行比较，以确定所述流体因子结果所在的流体因子阈值区间，并确定该流体因子阈值区间对应的储层流体类别。
需要说明的是，不同的储层流体类别，例如油、气和水可以对应不同的流体因子阈值区间；则在得到流体因子结果后，可根据其所在的流体因子阈值区间直接确定相应的储层流体类别。
此处，为了研究目标储层的储层流体类别，可从整个三维数据的流体因子结果中提取目标储层的流体因子沿层切片，例如如图7所示，从而便于分析切片上各点的含油气性，便于优选有利勘探的位置。
本发明实施例提供的一种基于弱度参数检测裂缝型储层流体的方法，首先获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。本发明可以实现准确可靠的检测裂缝型储层流体的类型的目的。
对应于上述如图1和图2所示的方法实施例，如图8所示，本发明实施例提供一种基于弱度参数检测裂缝型储层流体的装置，包括：
三维地震资料处理单元31，用于获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集。
各向异性参数求取单元32，用于根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数。
初始流体因子计算单元33，用于根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子。
第二级流体因子获得单元34，用于根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子。
流体因子结果获得单元35，用于根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果。
储层流体类别确定单元36，用于将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。
另外，所述初始流体因子计算单元33，具体用于：
获取裂缝介质的等效弹性系数矩阵：
Figure BDA0001787079660000151
其中，
Figure BDA0001787079660000152
表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵；λ和μ为各向同性背景介质拉梅常数；ΔT用于表示切向弱度参数；ΔN用于表示法向弱度参数；
Figure BDA0001787079660000161
将裂缝介质的等效弹性系数矩阵
Figure BDA0001787079660000162
中的元素代入到公式：
Figure BDA0001787079660000163
计算切向弱度参数ΔT和法向弱度参数ΔN；其中，
Figure BDA0001787079660000164
其中，cij表示裂缝介质的等效弹性系数矩阵
Figure BDA0001787079660000165
中的第i行第j列元素；εV、δV和γ表示各向异性参数。
将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子
Figure BDA0001787079660000166
其中，
Figure BDA0001787079660000167
另外，所述第二级流体因子获得单元34，具体用于：
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf；其中，0≈μf<<λf。
根据公式：
Figure BDA0001787079660000168
对所述初始流体因子
Figure BDA0001787079660000169
进行修正，获得第二级流体因子
Figure BDA00017870796600001610
其中，
Figure BDA00017870796600001611
g为采用背景介质的拉梅常数λ和μ计算得到的参数，
Figure BDA00017870796600001612
d表示裂缝型储层的裂缝纵横比。
另外，所述流体因子结果获得单元35，具体用于：
根据公式：
Figure BDA00017870796600001613
以及裂缝型储层的裂缝纵横比d，对所述第二级流体因子
Figure BDA00017870796600001614
进行修正，获得流体因子结果Ifluid。
其中，
Figure BDA00017870796600001615
另外，所述储层流体类别确定单元36中预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间。
则储层流体类别确定单元36，具体可以：
将所述流体因子结果与所述流体因子阈值区间进行比较，以确定所述流体因子结果所在的流体因子阈值区间，并确定该流体因子阈值区间对应的储层流体类别。
本发明实施例提供的一种基于弱度参数检测裂缝型储层流体的装置，首先获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。本发明可以实现准确可靠的检测裂缝型储层流体的类型的目的。
另外，本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现以下步骤：
获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集。
根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数。
根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子。
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子。
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果。
将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。
另外，本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现以下步骤：
获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集。
根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数。
根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子。
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子。
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果。
将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (6)

1.一种基于弱度参数检测裂缝型储层流体的方法，其特征在于，包括：
获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；
根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；
根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；
将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别；
所述预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间；
所述将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别，包括：
将所述流体因子结果与所述流体因子阈值区间进行比较，以确定所述流体因子结果所在的流体因子阈值区间，并确定该流体因子阈值区间对应的储层流体类别；
所述根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子，包括：
获取裂缝介质的等效弹性系数矩阵：
Figure FDA0002479870600000011
其中，
Figure FDA0002479870600000012
表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵；λ和μ为各向同性背景介质拉梅常数；△T用于表示切向弱度参数；△N用于表示法向弱度参数；
Figure FDA0002479870600000021
将裂缝介质的等效弹性系数矩阵
Figure FDA0002479870600000022
中的元素代入到公式：
Figure FDA0002479870600000023
计算切向弱度参数△T和法向弱度参数△N；其中，
Figure FDA0002479870600000024
其中，cij表示裂缝介质的等效弹性系数矩阵
Figure FDA0002479870600000025
中的第i行第j列元素；εV、δV和γ表示各向异性参数；
将切向弱度参数△T和法向弱度参数△N的比值确定为初始流体因子
Figure FDA0002479870600000026
其中，
Figure FDA0002479870600000027
所述根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子，包括：
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf；其中，0≈μf<<λf；
根据公式：
Figure FDA0002479870600000028
对所述初始流体因子
Figure FDA0002479870600000029
进行修正，获得第二级流体因子
Figure FDA00024798706000000210
其中，
Figure FDA00024798706000000211
g为采用背景介质的拉梅常数λ和μ计算得到的参数，
Figure FDA00024798706000000212
d表示裂缝型储层的裂缝纵横比。
2.根据权利要求1所述的基于弱度参数检测裂缝型储层流体的方法，其特征在于，所述根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果，包括：
根据公式：
Figure FDA00024798706000000213
以及裂缝型储层的裂缝纵横比d，对所述第二级流体因子
Figure FDA00024798706000000214
进行修正，获得流体因子结果Ifluid；
其中，
Figure FDA0002479870600000031
3.一种基于弱度参数检测裂缝型储层流体的装置，其特征在于，包括：
三维地震资料处理单元，用于获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；
各向异性参数求取单元，用于根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；
初始流体因子计算单元，用于根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；
第二级流体因子获得单元，用于根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；
流体因子结果获得单元，用于根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；
储层流体类别确定单元，用于将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别；
所述储层流体类别确定单元中预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间；
所述储层流体类别确定单元，具体用于：
将所述流体因子结果与所述流体因子阈值区间进行比较，以确定所述流体因子结果所在的流体因子阈值区间，并确定该流体因子阈值区间对应的储层流体类别；
所述初始流体因子计算单元，具体用于：
获取裂缝介质的等效弹性系数矩阵：
Figure FDA0002479870600000032
其中，
Figure FDA0002479870600000033
表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵；λ和μ为各向同性背景介质拉梅常数；△T用于表示切向弱度参数；△N用于表示法向弱度参数；
Figure FDA0002479870600000041
将裂缝介质的等效弹性系数矩阵
Figure FDA0002479870600000042
中的元素代入到公式：
Figure FDA0002479870600000043
计算切向弱度参数△T和法向弱度参数△N；其中，
Figure FDA0002479870600000044
其中，cij表示裂缝介质的等效弹性系数矩阵
Figure FDA0002479870600000045
中的第i行第j列元素；εV、δV和γ表示各向异性参数；
将切向弱度参数△T和法向弱度参数△N的比值确定为初始流体因子
Figure FDA0002479870600000046
其中，
Figure FDA0002479870600000047
所述第二级流体因子获得单元，具体用于：
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf；其中，0≈μf<<λf；
根据公式：
Figure FDA0002479870600000048
对所述初始流体因子
Figure FDA0002479870600000049
进行修正，获得第二级流体因子
Figure FDA00024798706000000410
其中，
Figure FDA00024798706000000411
g为采用背景介质的拉梅常数λ和μ计算得到的参数，
Figure FDA00024798706000000412
d表示裂缝型储层的裂缝纵横比。
4.根据权利要求3所述的基于弱度参数检测裂缝型储层流体的装置，其特征在于，所述流体因子结果获得单元，具体用于：
根据公式：
Figure FDA00024798706000000413
以及裂缝型储层的裂缝纵横比d，对所述第二级流体因子
Figure FDA00024798706000000414
进行修正，获得流体因子结果Ifluid；
其中，
Figure FDA00024798706000000415
5.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现以下步骤：
获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；
根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；
根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；
将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别；
所述预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间；
所述将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别，包括：
将所述流体因子结果与所述流体因子阈值区间进行比较，以确定所述流体因子结果所在的流体因子阈值区间，并确定该流体因子阈值区间对应的储层流体类别；
所述根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子，包括：
获取裂缝介质的等效弹性系数矩阵：
Figure FDA0002479870600000051
其中，
Figure FDA0002479870600000052
表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵；λ和μ为各向同性背景介质拉梅常数；△T用于表示切向弱度参数；△N用于表示法向弱度参数；
Figure FDA0002479870600000053
将裂缝介质的等效弹性系数矩阵
Figure FDA0002479870600000061
中的元素代入到公式：
Figure FDA0002479870600000062
计算切向弱度参数△T和法向弱度参数△N；其中，
Figure FDA0002479870600000063
其中，cij表示裂缝介质的等效弹性系数矩阵
Figure FDA0002479870600000064
中的第i行第j列元素；εV、δV和γ表示各向异性参数；
将切向弱度参数△T和法向弱度参数△N的比值确定为初始流体因子
Figure FDA0002479870600000065
其中，
Figure FDA0002479870600000066
所述根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子，包括：
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf；其中，0≈μf<<λf；
根据公式：
Figure FDA0002479870600000067
对所述初始流体因子
Figure FDA0002479870600000068
进行修正，获得第二级流体因子
Figure FDA0002479870600000069
其中，
Figure FDA00024798706000000610
g为采用背景介质的拉梅常数λ和μ计算得到的参数，
Figure FDA00024798706000000611
d表示裂缝型储层的裂缝纵横比。
6.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：
获得裂缝型储层的三维地震资料，并将三维地震资料处理为分方位的入射角道集；
根据所述入射角道集，采用振幅随方位角变化反演算法计算得到各向异性参数；
根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子；
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子；
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正，获得流体因子结果；
将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别；
所述预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间；
所述将所述流体因子结果与预先设置的阈值区间进行比较，以确定所述流体因子结果对应的储层流体类别，包括：
将所述流体因子结果与所述流体因子阈值区间进行比较，以确定所述流体因子结果所在的流体因子阈值区间，并确定该流体因子阈值区间对应的储层流体类别；
所述根据所述各向异性参数，计算得到裂缝型储层的切向弱度参数和法向弱度参数，并将切向弱度参数和法向弱度参数的比值确定为初始流体因子，包括：
获取裂缝介质的等效弹性系数矩阵：
Figure FDA0002479870600000071
其中，
Figure FDA0002479870600000072
表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵；λ和μ为各向同性背景介质拉梅常数；△T用于表示切向弱度参数；△N用于表示法向弱度参数；
Figure FDA0002479870600000073
将裂缝介质的等效弹性系数矩阵
Figure FDA0002479870600000074
中的元素代入到公式：
Figure FDA0002479870600000075
计算切向弱度参数△T和法向弱度参数△N；其中，
Figure FDA0002479870600000081
其中，cij表示裂缝介质的等效弹性系数矩阵
Figure FDA0002479870600000082
中的第i行第j列元素；εV、δV和γ表示各向异性参数；
将切向弱度参数△T和法向弱度参数△N的比值确定为初始流体因子
Figure FDA0002479870600000083
其中，
Figure FDA0002479870600000084
所述根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正，获得第二级流体因子，包括：
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf；其中，0≈μf<<λf；
根据公式：
Figure FDA0002479870600000085
对所述初始流体因子
Figure FDA0002479870600000086
进行修正，获得第二级流体因子
Figure FDA0002479870600000087
其中，
Figure FDA0002479870600000088
g为采用背景介质的拉梅常数λ和μ计算得到的参数，
Figure FDA0002479870600000089
d表示裂缝型储层的裂缝纵横比。
CN201811020121.7A 2018-09-03 2018-09-03 基于弱度参数检测裂缝型储层流体的方法及装置 Active CN108897043B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811020121.7A CN108897043B (zh) 2018-09-03 2018-09-03 基于弱度参数检测裂缝型储层流体的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811020121.7A CN108897043B (zh) 2018-09-03 2018-09-03 基于弱度参数检测裂缝型储层流体的方法及装置

Publications (2)

Publication Number Publication Date
CN108897043A CN108897043A (zh) 2018-11-27
CN108897043B true CN108897043B (zh) 2020-09-04

Family

ID=64359224

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811020121.7A Active CN108897043B (zh) 2018-09-03 2018-09-03 基于弱度参数检测裂缝型储层流体的方法及装置

Country Status (1)

Country Link
CN (1) CN108897043B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109709610B (zh) * 2018-12-27 2020-07-17 南方科技大学 一种岩石裂缝探测方法及系统
CN110161562B (zh) * 2019-04-29 2020-09-18 中南大学 倾斜横向各向同性介质中的裂缝弱度反演方法及系统
CN113312583A (zh) * 2021-05-26 2021-08-27 中南大学 裂缝性页岩气储层性能参数预测方法及系统
CN113534241B (zh) * 2021-07-09 2022-04-29 成都理工大学 一种方位频变流体因子识别储层的方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9476998B2 (en) * 2010-12-21 2016-10-25 Schlumberger Technology Corporation Hydraulic fracture characterization using borehole sonic data
CN103412327B (zh) * 2013-08-01 2016-04-06 中国石油化工股份有限公司胜利油田分公司地质科学研究院 一种裂缝性储层的粘弹性参数提取方法
WO2015145257A2 (en) * 2014-03-28 2015-10-01 Cgg Services Sa Method and apparatus for characterizing azimuthal anisotropy using statistical moments
CN104007462B (zh) * 2014-04-16 2016-09-28 彭玲丽 基于衰减各向异性的裂缝预测方法
CN104018826B (zh) * 2014-04-16 2016-05-25 彭玲丽 基于非线性算法的裂缝及裂缝内流体识别方法
CN104375182B (zh) * 2014-11-10 2017-03-08 中国石油天然气股份有限公司 一种裂缝型储层流体的识别方法及装置
CN104614764B (zh) * 2015-01-28 2018-03-30 中国石油大学(华东) 一种基于扰动弹性阻抗反演的裂缝储层流体识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
珠江口盆地白云深水地区轻质原油综合预测技术;刘杰，等;《物探化探计算技术》;20180531;第313页 *

Also Published As

Publication number Publication date
CN108897043A (zh) 2018-11-27

Similar Documents

Publication Publication Date Title
CN108897043B (zh) 基于弱度参数检测裂缝型储层流体的方法及装置
US7679993B2 (en) Method of characterizing a fractured reservoir using seismic reflection amplitudes
CN102096107B (zh) 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法
US9244188B2 (en) System and method for estimating a nuclear magnetic resonance relaxation time cutoff
CN111123354B (zh) 基于频变反射振幅衰减预测致密气层的方法及设备
CN110632652B (zh) 断块油气藏的油气运移方式确定方法及装置
CN104005760A (zh) 基于方位各向异性弹性阻抗的裂缝检测方法
Falcon-Suarez et al. Experimental assessment of pore fluid distribution and geomechanical changes in saline sandstone reservoirs during and after CO2 injection
CN113312583A (zh) 裂缝性页岩气储层性能参数预测方法及系统
CN106323835A (zh) 确定非均质碳酸盐岩储层胶结指数的方法
CN107831540A (zh) 储层物性参数直接提取新方法
CN111381292A (zh) 一种预测砂岩含烃储层的测井解释方法与装置
CN112305602A (zh) 基于叠前多属性与古地貌融合技术的碳酸岩储层预测方法
US11530972B2 (en) Analyzing fractured rock samples
US11187078B2 (en) Method for evaluating connectivity between a first well and a second well in a hydrocarbon production field and related system
CN111241652B (zh) 一种确定地层原油粘度的方法及装置
CN106285655B (zh) 砂岩地层测井评估方法
CN108089228B (zh) 一种确定地层岩石特性的解释数据方法和装置
CN111159917A (zh) 基于岩石物理模型的各向异性流体因子构建方法
Mukherjee et al. Estimation of in-situ stress fields from P-wave seismic data
CN114876454A (zh) 一种大斜度井复杂岩性水淹层识别方法及系统
US20110196613A1 (en) Method of characterizing a co2 plume in a geological storage aquifer
CN115327630A (zh) 声波测井校正方法、装置、设备、存储介质及程序产品
CN111273354A (zh) 地下开采引起的地下水变化分析方法和系统
CN113704950B (zh) 一种气测全烃标准化评价储层含气性的方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant