CN108897043A - 基于弱度参数检测裂缝型储层流体的方法及装置 - Google Patents
基于弱度参数检测裂缝型储层流体的方法及装置 Download PDFInfo
- Publication number
- CN108897043A CN108897043A CN201811020121.7A CN201811020121A CN108897043A CN 108897043 A CN108897043 A CN 108897043A CN 201811020121 A CN201811020121 A CN 201811020121A CN 108897043 A CN108897043 A CN 108897043A
- Authority
- CN
- China
- Prior art keywords
- fluid
- parameter
- weakness
- fluid factor
- factor
- 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.)
- Granted
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 320
- 238000003379 elimination reaction Methods 0.000 title claims abstract description 124
- 230000008030 elimination Effects 0.000 title claims abstract description 118
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims description 25
- 238000004590 computer program Methods 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 abstract description 3
- 235000013399 edible fruits Nutrition 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000001615 p wave Methods 0.000 description 3
- 238000010008 shearing Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
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. for interpretation or for event detection
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于弱度参数检测裂缝型储层流体的方法及装置,涉及石油勘探技术领域。方法包括:获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;根据入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;根据各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;根据裂缝型储层背景介质的拉梅常数对初始流体因子进行修正,获得第二级流体因子;根据裂缝型储层的裂缝纵横比对第二级流体因子进行修正,获得流体因子结果;将流体因子结果与预先设置的阈值区间进行比较,以确定流体因子结果对应的储层流体类别。
Description
技术领域
本发明涉及石油勘探技术领域,尤其涉及一种基于弱度参数检测裂缝型储层流体的方法及装置。
背景技术
目前,裂缝型油气藏的产量占全世界石油天然气总产量的一半以上,在国内的油气勘探中,裂缝型油气藏是一个重要的勘探领域,如塔里木盆地奥陶系碳酸盐岩裂缝型油气藏,四川致密砂岩裂缝型油气藏。裂缝既是油、气、水等地下流体的运移通道,也是流体的储集空间,裂缝能极大地提升储层的渗透率。裂缝型储层描述包括预测裂缝分布特征和检测裂隙充填物,也就是探测开启并充填流体的有效裂缝。
利用地震资料研究岩石裂隙所贮存流体的类型和特性,是根据岩石物理等效介质理论,将与储层流体有关的异常特性表征为流体因子,再依托流体因子实现储层流体的类型判识。对于各向同性介质,现有技术中通常采用弹性参数的代数组合作为流体因子,常见的有LMR(Lumbda-Mu-Rho)法、Russell法、纵横波阻抗的幂次方等。流体检测的质量主要取决于两方面,一是反演计算的弹性参数是否可靠;二是构建的流体因子对流体类型是否敏感且稳定。可见,对于裂缝型储层,应该用不随方位变化的各向异性参数来构建流体因子。目前,如何能够准确可靠的检测裂缝型储层流体的类型成为了一个亟待解决的问题。
发明内容
本发明的实施例提供一种基于弱度参数检测裂缝型储层流体的方法及装置,以实现准确可靠的检测裂缝型储层流体的类型。
为达到上述目的,本发明采用如下技术方案:
一种基于弱度参数检测裂缝型储层流体的方法,包括:
获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;
根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;
根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;
将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
具体的,所述根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子,包括:
获取裂缝介质的等效弹性系数矩阵:
其中,表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵;λ和μ为各向同性背景介质拉梅常数;ΔT用于表示切向弱度参数;ΔN用于表示法向弱度参数;
将裂缝介质的等效弹性系数矩阵中的元素代入到公式:
计算切向弱度参数ΔT和法向弱度参数ΔN;其中,其中,cij表示裂缝介质的等效弹性系数矩阵中的第i行第j列元素;εV、δV和γ表示各向异性参数;
将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子其中,
具体的,所述根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子,包括:
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf;其中,0≈μf<<λf;
根据公式:对所述初始流体因子进行修正,获得第二级流体因子其中,g为采用背景介质的拉梅常数λ和μ计算得到的参数,d表示裂缝型储层的裂缝纵横比。
具体的,所述根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果,包括:
根据公式:以及裂缝型储层的裂缝纵横比d,对所述第二级流体因子进行修正,获得流体因子结果Ifluid;
其中,
具体的,所述预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间;
所述将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别,包括:
将所述流体因子结果与所述流体因子阈值区间进行比较,以确定所述流体因子结果所在的流体因子阈值区间,并确定该流体因子阈值区间对应的储层流体类别。
一种基于弱度参数检测裂缝型储层流体的装置,包括:
三维地震资料处理单元,用于获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;
各向异性参数求取单元,用于根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;
初始流体因子计算单元,用于根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;
第二级流体因子获得单元,用于根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;
流体因子结果获得单元,用于根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;
储层流体类别确定单元,用于将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
另外,所述初始流体因子计算单元,具体用于:
获取裂缝介质的等效弹性系数矩阵:
其中,表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵;λ和μ为各向同性背景介质拉梅常数;ΔT用于表示切向弱度参数;ΔN用于表示法向弱度参数;
将裂缝介质的等效弹性系数矩阵中的元素代入到公式:
计算切向弱度参数ΔT和法向弱度参数ΔN;其中,其中,cij表示裂缝介质的等效弹性系数矩阵中的第i行第j列元素;εV、δV和γ表示各向异性参数;
将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子其中,
另外,所述第二级流体因子获得单元,具体用于:
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf;其中,0≈μf<<λf;
根据公式:对所述初始流体因子进行修正,获得第二级流体因子其中,g为采用背景介质的拉梅常数λ和μ计算得到的参数,d表示裂缝型储层的裂缝纵横比。
另外,所述流体因子结果获得单元,具体用于:
根据公式:以及裂缝型储层的裂缝纵横比d,对所述第二级流体因子进行修正,获得流体因子结果Ifluid;
其中,
另外,所述储层流体类别确定单元中预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间;
所述储层流体类别确定单元,具体用于:
将所述流体因子结果与所述流体因子阈值区间进行比较,以确定所述流体因子结果所在的流体因子阈值区间,并确定该流体因子阈值区间对应的储层流体类别。
一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现以下步骤:
获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;
根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;
根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;
将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
一种计算机设备,包括存储器、处理器及存储在存储上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;
根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;
根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;
将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
本发明实施例提供的一种基于弱度参数检测裂缝型储层流体的方法及装置,首先获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。本发明可以实现准确可靠的检测裂缝型储层流体的类型的目的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图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介质分界面处的纵波反射系数是各向异性参数、入射角和方位角的函数:
上式中:φs为裂缝平面的对称轴方向,即裂缝面法向,(φ-φs)为测线与裂缝面法向的夹角;α,β,Z,G分别是纵波速度、横波速度、纵波阻抗和剪切模量;εV,δV,γ为各向异性参数。分别表示上下两层介质间纵波速度的平均值、横波速度的平均值、纵波阻抗的平均值和剪切模量的平均值,Δα,ΔZ,ΔG,ΔεV,ΔδV,Δγ分别表示上下两层介质间纵波速度的差值、纵波阻抗的差值、剪切模量的差值以及三个各向异性参数的差值。
对上式进行简化,地震纵波的反射振幅R与入射角θ、方位角φ以及各向异性参数之间的关系式如下:
其中A,B,C,D,E,F是入射角θ和方位角φ的组合:
各向异性参数的反演方程组包括下述超定方程组:
Mn×6X6×1=Rn×1
上式中n为覆盖次数,M,X和R分别为:
设定阻尼因子,然后利用奇异值分解法求解所述超定方程组,求解的结果为:
X=M+R
其中M+为矩阵M的广义逆。
通过求解得到的X,可计算第i个地层分界面的各向异性参数差值:
Δγi=X(6)
从上到下累加各向异性参数差值,可得第n个地层的各向异性参数:
对于裂缝型储层,参数εV,δV,γ描述了介质的各向异性程度。
步骤203、获取裂缝介质的等效弹性系数矩阵:
其中,表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵;λ和μ为各向同性背景介质拉梅常数,可通过AVO(振幅随偏移距的变化,Amplitude variationwith offset)三参数反演方法计算而得到;ΔT用于表示切向弱度参数;ΔN用于表示法向弱度参数;
步骤204、将裂缝介质的等效弹性系数矩阵中的元素代入到公式:
则有:
从而可计算切向弱度参数ΔT和法向弱度参数ΔN。其中, 其中,cij表示裂缝介质的等效弹性系数矩阵中的第i行第j列元素;εV、δV和γ表示各向异性参数。
对于计算得到的切向弱度参数ΔT和法向弱度参数ΔN,如图4所示,裂缝型储层充填不同流体时,切向弱度参数ΔT和法向弱度参数ΔN随裂缝密度变化(crack density表示裂缝密度)。
步骤205、将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子其中,
对于初始流体因子如图5所示,初始流体因子随裂缝密度变化很小。
步骤206、获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf。
其中,0≈μf<<λf。
步骤207、根据公式:对所述初始流体因子进行修正,获得第二级流体因子
其中,g为采用背景介质的拉梅常数λ和μ计算得到的参数,d表示裂缝型储层的裂缝纵横比,可通过地下介质裂隙结构反演计算得到。
此处,对于初始流体因子的公式,原理如下:
在地下储层含有裂缝时,其弹性系数矩阵C66可由Schoenberg模型或者Hudson模型求得。
Schoenberg模型也称为线性滑动模型,是指裂缝表面位移不连续,位移与穿过裂缝的应力呈线性关系的模型。设本构坐标系下裂隙的法向与x轴平行,裂隙密度为e,纵横比为d,各向同性背景的拉梅(Lame)常数为λ和μ,裂隙填充物的拉梅常数为λf和μf。引入切向弱度参数ΔT和法向弱度参数ΔN,裂缝介质的等效弹性系数矩阵如公式(1):
其中:
弹性系数矩阵C66为6×6阶的矩阵,上标Schoenberg表示Schoenberg模型。
而Hudson模型中,假定介质包含比地震波长小得多的定向稀疏排列裂隙,裂隙是彼此分离的薄片,裂隙内包含的气体、流体等物质的体积模量比背景介质小。裂缝介质的等效弹性系数矩阵为公式(2):
其中:
其中,U11和U33表达式如公式(2.b):
其中M、k表达式如公式(2.c),与充填流体的拉梅常数相关:
对应相同的裂缝介质,Hudson模型与Schoenberg模型所描述的弹性系数矩阵相等,即
由和得公式(3):
化简后得公式(3.a):
将公式(2.b)的U11,U33代入公式(3.a),得到公式(3.b):
而将切向弱度参数ΔT和法向弱度参数ΔN的比值记为初始流体因子
流体的第二拉梅常数μf接近于零值,相比第一拉梅常数λf小得多,即有:
0≈μf<<λf
依据公式(2.c),有:
0≈M<<κ
此时,公式(4)可变换为公式(4.a):
步骤208、根据公式:以及裂缝型储层的裂缝纵横比d,对所述第二级流体因子进行修正,获得流体因子结果Ifluid。
此处,由于:
0≈μf<<λf;
因此存在
对于流体因子结果中的流体因子,其阈值区间的划分如图6所示,其中,probability是指概率,hydrocarbon indicator是指烃类指示,说明了流体因子阈值区间是如何划分的。在井点处,有纵横波速度和密度、可以计算得到背景介质的拉梅常数(λ、μ);也能解释得到裂缝密度e和裂缝纵横比d;从而计算获流体因子的值;把这些值(很多散点)生成直方图。当各种流体的直方图没有重叠时,能很轻易的划分阈值区间;当直方图存在部分重合时,则在允许的错误识别率之下,根据概率分布图划定阈值区间。
步骤209、将流体因子结果与流体因子阈值区间进行比较,以确定所述流体因子结果所在的流体因子阈值区间,并确定该流体因子阈值区间对应的储层流体类别。
需要说明的是,不同的储层流体类别,例如油、气和水可以对应不同的流体因子阈值区间;则在得到流体因子结果后,可根据其所在的流体因子阈值区间直接确定相应的储层流体类别。
此处,为了研究目标储层的储层流体类别,可从整个三维数据的流体因子结果中提取目标储层的流体因子沿层切片,例如如图7所示,从而便于分析切片上各点的含油气性,便于优选有利勘探的位置。
本发明实施例提供的一种基于弱度参数检测裂缝型储层流体的方法,首先获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。本发明可以实现准确可靠的检测裂缝型储层流体的类型的目的。
对应于上述如图1和图2所示的方法实施例,如图8所示,本发明实施例提供一种基于弱度参数检测裂缝型储层流体的装置,包括:
三维地震资料处理单元31,用于获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集。
各向异性参数求取单元32,用于根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数。
初始流体因子计算单元33,用于根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子。
第二级流体因子获得单元34,用于根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子。
流体因子结果获得单元35,用于根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果。
储层流体类别确定单元36,用于将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
另外,所述初始流体因子计算单元33,具体用于:
获取裂缝介质的等效弹性系数矩阵:
其中,表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵;λ和μ为各向同性背景介质拉梅常数;ΔT用于表示切向弱度参数;ΔN用于表示法向弱度参数;
将裂缝介质的等效弹性系数矩阵中的元素代入到公式:
计算切向弱度参数ΔT和法向弱度参数ΔN;其中,其中,cij表示裂缝介质的等效弹性系数矩阵中的第i行第j列元素;εV、δV和γ表示各向异性参数。
将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子其中,
另外,所述第二级流体因子获得单元34,具体用于:
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf;其中,0≈μf<<λf。
根据公式:对所述初始流体因子进行修正,获得第二级流体因子其中,g为采用背景介质的拉梅常数λ和μ计算得到的参数,d表示裂缝型储层的裂缝纵横比。
另外,所述流体因子结果获得单元35,具体用于:
根据公式:以及裂缝型储层的裂缝纵横比d,对所述第二级流体因子进行修正,获得流体因子结果Ifluid。
其中,
另外,所述储层流体类别确定单元36中预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间。
则储层流体类别确定单元36,具体可以:
将所述流体因子结果与所述流体因子阈值区间进行比较,以确定所述流体因子结果所在的流体因子阈值区间,并确定该流体因子阈值区间对应的储层流体类别。
本发明实施例提供的一种基于弱度参数检测裂缝型储层流体的装置,首先获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。本发明可以实现准确可靠的检测裂缝型储层流体的类型的目的。
另外,本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现以下步骤:
获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集。
根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数。
根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子。
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子。
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果。
将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
另外,本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集。
根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数。
根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子。
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子。
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果。
将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (12)
1.一种基于弱度参数检测裂缝型储层流体的方法,其特征在于,包括:
获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;
根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;
根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;
将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
2.根据权利要求1所述的基于弱度参数检测裂缝型储层流体的方法,其特征在于,所述根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子,包括:
获取裂缝介质的等效弹性系数矩阵:
其中,表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵;λ和μ为各向同性背景介质拉梅常数;ΔT用于表示切向弱度参数;ΔN用于表示法向弱度参数;
将裂缝介质的等效弹性系数矩阵中的元素代入到公式:
计算切向弱度参数ΔT和法向弱度参数ΔN;其中,其中,cij表示裂缝介质的等效弹性系数矩阵中的第i行第j列元素;εV、δV和γ表示各向异性参数;
将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子其中,
3.根据权利要求2所述的基于弱度参数检测裂缝型储层流体的方法,其特征在于,所述根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子,包括:
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf;其中,0≈μf<<λf;
根据公式:对所述初始流体因子进行修正,获得第二级流体因子其中,g为采用背景介质的拉梅常数λ和μ计算得到的参数,d表示裂缝型储层的裂缝纵横比。
4.根据权利要求3所述的基于弱度参数检测裂缝型储层流体的方法,其特征在于,所述根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果,包括:
根据公式:以及裂缝型储层的裂缝纵横比d,对所述第二级流体因子进行修正,获得流体因子结果Ifluid;
其中,
5.根据权利要求4所述的基于弱度参数检测裂缝型储层流体的方法,其特征在于,所述预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间;
所述将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别,包括:
将所述流体因子结果与所述流体因子阈值区间进行比较,以确定所述流体因子结果所在的流体因子阈值区间,并确定该流体因子阈值区间对应的储层流体类别。
6.一种基于弱度参数检测裂缝型储层流体的装置,其特征在于,包括:
三维地震资料处理单元,用于获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;
各向异性参数求取单元,用于根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;
初始流体因子计算单元,用于根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;
第二级流体因子获得单元,用于根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;
流体因子结果获得单元,用于根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;
储层流体类别确定单元,用于将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
7.根据权利要求6所述的基于弱度参数检测裂缝型储层流体的装置,其特征在于,所述初始流体因子计算单元,具体用于:
获取裂缝介质的等效弹性系数矩阵:
其中,表示Schoenberg模型下的6乘6阶裂缝介质的等效弹性系数矩阵;λ和μ为各向同性背景介质拉梅常数;ΔT用于表示切向弱度参数;ΔN用于表示法向弱度参数;
将裂缝介质的等效弹性系数矩阵中的元素代入到公式:
计算切向弱度参数ΔT和法向弱度参数ΔN;其中,其中,cij表示裂缝介质的等效弹性系数矩阵中的第i行第j列元素;εV、δV和γ表示各向异性参数;
将切向弱度参数ΔT和法向弱度参数ΔN的比值确定为初始流体因子其中,
8.根据权利要求7所述的基于弱度参数检测裂缝型储层流体的装置,其特征在于,所述第二级流体因子获得单元,具体用于:
获得裂缝型储层中裂隙填充物的第一拉梅常数λf和第二拉梅常数μf;其中,0≈μf<<λf;
根据公式:对所述初始流体因子进行修正,获得第二级流体因子其中,g为采用背景介质的拉梅常数λ和μ计算得到的参数,d表示裂缝型储层的裂缝纵横比。
9.根据权利要求8所述的基于弱度参数检测裂缝型储层流体的装置,其特征在于,所述流体因子结果获得单元,具体用于:
根据公式:以及裂缝型储层的裂缝纵横比d,对所述第二级流体因子进行修正,获得流体因子结果Ifluid;
其中,
10.根据权利要求9所述的基于弱度参数检测裂缝型储层流体的装置,其特征在于,所述储层流体类别确定单元中预先设置的阈值区间包括通过测井资料求得的储层流体类别分别为油、气和水对应的流体因子阈值区间;
所述储层流体类别确定单元,具体用于:
将所述流体因子结果与所述流体因子阈值区间进行比较,以确定所述流体因子结果所在的流体因子阈值区间,并确定该流体因子阈值区间对应的储层流体类别。
11.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现以下步骤:
获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;
根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;
根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;
将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
12.一种计算机设备,包括存储器、处理器及存储在存储上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现以下步骤:
获得裂缝型储层的三维地震资料,并将三维地震资料处理为分方位的入射角道集;
根据所述入射角道集,采用振幅随方位角变化反演算法计算得到各向异性参数;
根据所述各向异性参数,计算得到裂缝型储层的切向弱度参数和法向弱度参数,并将切向弱度参数和法向弱度参数的比值确定为初始流体因子;
根据裂缝型储层背景介质的拉梅常数对所述初始流体因子进行修正,获得第二级流体因子;
根据裂缝型储层的裂缝纵横比对所述第二级流体因子进行修正,获得流体因子结果;
将所述流体因子结果与预先设置的阈值区间进行比较,以确定所述流体因子结果对应的储层流体类别。
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 true CN108897043A (zh) | 2018-11-27 |
CN108897043B 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) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109709610A (zh) * | 2018-12-27 | 2019-05-03 | 南方科技大学 | 一种岩石裂缝探测方法及系统 |
CN110161562A (zh) * | 2019-04-29 | 2019-08-23 | 中国石油大学(华东) | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
CN113312583A (zh) * | 2021-05-26 | 2021-08-27 | 中南大学 | 裂缝性页岩气储层性能参数预测方法及系统 |
CN113534241A (zh) * | 2021-07-09 | 2021-10-22 | 成都理工大学 | 一种方位频变流体因子识别储层的方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012087796A2 (en) * | 2010-12-21 | 2012-06-28 | Schlumberger Canada Limited | Hydraulic fracture characterization using borehole sonic data |
CN103412327A (zh) * | 2013-08-01 | 2013-11-27 | 中国石油化工股份有限公司胜利油田分公司地质科学研究院 | 一种裂缝性储层的粘弹性参数提取方法 |
CN104007462A (zh) * | 2014-04-16 | 2014-08-27 | 孙赞东 | 基于衰减各向异性的裂缝预测方法 |
CN104018826A (zh) * | 2014-04-16 | 2014-09-03 | 孙赞东 | 基于非线性算法的裂缝及裂缝内流体识别方法 |
CN104375182A (zh) * | 2014-11-10 | 2015-02-25 | 中国石油天然气股份有限公司 | 一种裂缝型储层流体的识别方法及装置 |
CN104614764A (zh) * | 2015-01-28 | 2015-05-13 | 中国石油大学(华东) | 一种基于扰动弹性阻抗反演的裂缝储层流体识别方法 |
WO2015145257A2 (en) * | 2014-03-28 | 2015-10-01 | Cgg Services Sa | Method and apparatus for characterizing azimuthal anisotropy using statistical moments |
-
2018
- 2018-09-03 CN CN201811020121.7A patent/CN108897043B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012087796A2 (en) * | 2010-12-21 | 2012-06-28 | Schlumberger Canada Limited | Hydraulic fracture characterization using borehole sonic data |
CN103412327A (zh) * | 2013-08-01 | 2013-11-27 | 中国石油化工股份有限公司胜利油田分公司地质科学研究院 | 一种裂缝性储层的粘弹性参数提取方法 |
WO2015145257A2 (en) * | 2014-03-28 | 2015-10-01 | Cgg Services Sa | Method and apparatus for characterizing azimuthal anisotropy using statistical moments |
CN104007462A (zh) * | 2014-04-16 | 2014-08-27 | 孙赞东 | 基于衰减各向异性的裂缝预测方法 |
CN104018826A (zh) * | 2014-04-16 | 2014-09-03 | 孙赞东 | 基于非线性算法的裂缝及裂缝内流体识别方法 |
CN104375182A (zh) * | 2014-11-10 | 2015-02-25 | 中国石油天然气股份有限公司 | 一种裂缝型储层流体的识别方法及装置 |
CN104614764A (zh) * | 2015-01-28 | 2015-05-13 | 中国石油大学(华东) | 一种基于扰动弹性阻抗反演的裂缝储层流体识别方法 |
Non-Patent Citations (4)
Title |
---|
罗辑,等: "基于方位弹性阻抗反演的裂缝型储层流体检测方法", 《石油地球物理勘探》 * |
薛姣,等: "基于等效介质模型的裂缝参数AVOA反演", 《石油地球物理勘探》 * |
谢春辉,等: "基于AVAZ反演的裂缝型储层流体识别方法", 《2015年物探技术研讨会》 * |
陈怀震,等: "基于各向异性岩石物理的缝隙流体因子AVAZ反演", 《地球物理学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109709610A (zh) * | 2018-12-27 | 2019-05-03 | 南方科技大学 | 一种岩石裂缝探测方法及系统 |
CN109709610B (zh) * | 2018-12-27 | 2020-07-17 | 南方科技大学 | 一种岩石裂缝探测方法及系统 |
CN110161562A (zh) * | 2019-04-29 | 2019-08-23 | 中国石油大学(华东) | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
CN110161562B (zh) * | 2019-04-29 | 2020-09-18 | 中南大学 | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
CN113312583A (zh) * | 2021-05-26 | 2021-08-27 | 中南大学 | 裂缝性页岩气储层性能参数预测方法及系统 |
CN113534241A (zh) * | 2021-07-09 | 2021-10-22 | 成都理工大学 | 一种方位频变流体因子识别储层的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108897043B (zh) | 2020-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Herwanger | Seismic geomechanics: how to build and calibrate geomechanical models using 3D and 4D seismic data | |
Cipolla et al. | Seismic-to-simulation for unconventional reservoir development | |
Shuey | A simplification of the Zoeppritz equations | |
Rich et al. | Unconventional geophysics for unconventional plays | |
US20200132869A1 (en) | Seismic rock physics inversion method based on large area tight reservoir | |
CN102788994B (zh) | 一种储层裂缝的确定方法 | |
CN103257361B (zh) | 基于Zoeppritz方程近似式的油气预测方法及系统 | |
CN107045143B (zh) | 一种预测裂缝发育的方法及装置 | |
Szydlik et al. | 3D PP/PS prestack depth migration on the Volve field | |
CN105068117A (zh) | 用于裂缝性介质的avo反演方法、装置和设备 | |
CN108897043A (zh) | 基于弱度参数检测裂缝型储层流体的方法及装置 | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
CN104375182A (zh) | 一种裂缝型储层流体的识别方法及装置 | |
CN104570128A (zh) | 一种基于岩石脆性因子的非常规油气储层预测方法 | |
CN104614764A (zh) | 一种基于扰动弹性阻抗反演的裂缝储层流体识别方法 | |
Yasin et al. | Fault and fracture network characterization using seismic data: a study based on neural network models assessment | |
Rodríguez-Pradilla | Microseismic monitoring of a hydraulic-fracturing operation in a CBM reservoir: Case study in the Cerrejón Formation, Cesar-Ranchería Basin, Colombia | |
CN102590863A (zh) | 确定反射系数频散、反射能量及识别气层的方法及装置 | |
Sun et al. | Fractured reservoir modeling by discrete fracture network and seismic modeling in the Tarim Basin, China | |
Sharifi et al. | Full-angle extended elastic impedance | |
Taleghani et al. | Post-treatment assessment of induced fracture networks | |
Gong et al. | Application of multi-level and high-resolution fracture modeling in field-scale reservoir simulation study | |
Amanipoor et al. | The application of géostatistical methods to prepare the 3D petrophysical model of oil reservoir | |
Ali et al. | Seismic modelling of a fractured carbonate reservoir in Abu Dhabi, United Arab Emirates | |
Bejaoui et al. | Characterization and simulation of a complex fractured carbonate field offshore Tunisia |
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 |