CN104975853B - 一种获取地层孔隙度的方法及装置 - Google Patents
一种获取地层孔隙度的方法及装置 Download PDFInfo
- Publication number
- CN104975853B CN104975853B CN201510368565.XA CN201510368565A CN104975853B CN 104975853 B CN104975853 B CN 104975853B CN 201510368565 A CN201510368565 A CN 201510368565A CN 104975853 B CN104975853 B CN 104975853B
- Authority
- CN
- China
- Prior art keywords
- time difference
- formation
- mrow
- msub
- slowness
- 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
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 130
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000002474 experimental method Methods 0.000 claims description 65
- 238000004364 calculation method Methods 0.000 claims description 20
- 239000012530 fluid Substances 0.000 claims description 11
- 210000004209 hair Anatomy 0.000 claims description 8
- 238000005056 compaction Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 239000011435 rock Substances 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 2
- 230000006870 function Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000003860 storage Methods 0.000 description 4
- 239000004568 cement Substances 0.000 description 3
- 230000005251 gamma ray Effects 0.000 description 3
- 238000003780 insertion Methods 0.000 description 3
- 230000037431 insertion Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例提供了一种获取地层孔隙度的方法及装置,所述方法包括:获取目标地层的声波测井资料、井斜角曲线资料以及层位时差各向异性比;根据所述目标地层的井斜角曲线资料以及层位时差各向异性比,按照拟合方程计算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差,所述拟合方程是根据已知的实验地层训练得到的;根据所述目标地层声波测井资料中的时差以及所述目标地层的慢度差计算所述目标地层的水平时差;根据所述目标地层的水平时差计算所述目标地层的孔隙度。采用本发明实施例提供的方法可以计算得到目标地层的水平时差,由水平时差计算得到的地层孔隙度更加接近于地层的实际孔隙度。
Description
技术领域
本发明属于石油勘探资料处理领域,特别涉及获取地层孔隙度的方法及装置。
背景技术
地层孔隙度是储层评价的重要参数之一,地层孔隙度与测井电位、钻井时间、油气储量以及含油饱和度等石油勘探参数紧密相连,因此,获取地层孔隙度是石油勘探的一个重要工作步骤。现有技术中,通常根据电缆声波测得的时差计算得到所述地层孔隙度。
在实现本申请过程中,发明人发现现有技术中至少存在如下问题:在很多地层中,尤其在大斜度井以及水平井中,存在时差各向异性现象。这时电缆声波测得的时差是水平时差和垂直时差的综合值。而各向异性现象多数是由于地层内夹杂着具有不同时差的其他地层而产生的。计算所述孔隙度所需要的时差一般是水平时差,因此使用电缆声波测得的时差计算得到的地层孔隙度偏离真实地层孔隙度,进而根据所述地层孔隙度计算得到的含油饱和度等参数值误差也相对较大。
发明内容
针对现有技术中根据声波测井仪器测得的时差计算得到的孔隙度误差较大,本发明提出一种获取地层孔隙度的方法及装置,所述获取地层孔隙度的方法及装置具体是这样实现的:
一种获取地层孔隙度的方法,所述方法包括:
获取目标地层的声波测井资料、井斜角曲线资料以及层位时差各向异性比;
根据所述目标地层的井斜角曲线资料以及层位时差各向异性比,按照拟合方程计算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差,所述拟合方程根据已知的实验地层训练得到;
根据所述目标地层声波测井资料中的时差以及所述目标地层的慢度差计算所述目标地层的水平时差;
根据所述目标地层的水平时差计算所述目标地层的孔隙度。
可选的,在本发明一实施例中,所述根据已知的实验地层训练得到拟合方程,包括:
获取所述实验地层的水平时差,计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差;
根据所述实验地层的水平时差,计算所述实验地层的慢度差,所述慢度差为所述时差与所述实验地层的水平时差的差值;
建立所述井斜角以及所述慢度差的坐标系,在所述坐标系中插入所述不同井斜角、不同时差各向异性比条件下的所述慢度差的离散点;
将所述离散点拟合成曲线,并获取所述曲线的拟合方程。
可选的,在本发明一实施例中,所述计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差,包括:
利用声波测井仪采用有限差分法计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差,其中,所述井斜角从0到90°范围内按15°步长递增,所述时差各向异性比从1.0到1.2范围内按0.02步长递增。
可选的,在本发明一实施例中,所述声波测井仪包括双发双收声波测井仪。
可选的,在本发明一实施例中,所述建立所述井斜角以及所述慢度差的坐标系,具体包括:
建立所述井斜角正弦值以及所述慢度差的直角坐标系。
可选的,在本发明一实施例中,所述拟合方程包括一元二次方程。
可选的,在本发明一实施例中,所述按照所述拟合方程,计算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差,具体包括:
根据所述实验地层井斜角以及慢度差的坐标系中的曲线,按照所述拟合方程,计算与所述实验地层相同的时差各向异性比条件下,所述目标地层在不同井斜角条件下的慢度差,当所述井斜角是步长15°的倍数时,按照所述拟合方程求解所述目标地层的慢度差;当所述井斜角不是15°的倍数时,采用插值法计算所述目标地层的慢度差。
可选的,在本发明一实施例中,根据所述目标地层的水平时差计算所述目标地层的孔隙度,其中所述孔隙度的计算公式为:
式中,Φt为孔隙度,Δtma、Δtf、Δtsh分别为岩石骨架声波水平方向时差、地层流体声波水平方向时差、地层流体声波水平方向时差,Vsh为地层泥质含量,Cp为声波压实校正系数,Δt为所述目的地层的水平时差。
一种获取地层孔隙度的装置,包括:
资料获取模块,用于获取目标地层的声波测井资料、井斜角曲线资料以及层位时差各向异性比;
目标慢度差计算模块,用于根据所述目标地层的井斜角曲线资料以及层位时差各向异性比,按照拟合方程计算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差,所述拟合方程根据拟合装置计算得到;
水平时差计算模块,用于根据所述目标地层声波测井资料中的时差以及所述目标地层的慢度差计算所述目标地层的水平时差;
孔隙度计算模块,用于根据所述目标地层的水平时差计算所述目标地层的孔隙度。
可选的,在本发明一实施例中,所述目标慢度差计算模块还包括拟合装置,所述拟合装置包括:
实验时差计算模块,用于获取实验地层的水平时差,计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差;
慢度差计算模块,用于根据所述实验地层的水平时差,计算所述实验地层的慢度差,所述慢度差为所述时差与所述实验地层的水平时差的差值;
坐标系建立模块,用于建立所述井斜角以及所述慢度差的坐标系,在所述坐标系中插入所述不同井斜角、不同时差各向异性比条件下的所述慢度差的离散点;
拟合模块,用于将所述离散点拟合成曲线,并获取所述曲线的拟合方程。
可选的,在本发明一实施例中,所述实验时差计算模块,还包括:
时差计算子模块,用于利用声波测井仪采用有限差分法计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差,其中,所述井斜角从0到90°范围内按15°步长递增,所述时差各向异性比从1.0到1.2范围内按0.02步长递增。
可选的,在本发明一实施例中,所述声波测井仪包括双发双收声波测井仪。
可选的,在本发明一实施例中,所述坐标系建立模块,还包括:
坐标系建立子模块,用于建立所述井斜角正弦值以及所述慢度差的直角坐标系。
可选的,在本发明一实施例中,所述拟合方程包括一元二次方程。
可选的,在本发明一实施例中,所述目标慢度差计算模块,还包括:
目标慢度差计算子模块,用于根据所述实验地层井斜角以及慢度差的坐标系中的曲线,按照所述拟合方程,计算与所述实验地层相同的时差各向异性比条件下,所述目标地层在不同井斜角条件下的慢度差,当所述井斜角是步长15°的倍数时,按照所述拟合方程求解所述目标地层的慢度差;当所述井斜角不是15°的倍数时,采用插值法计算所述目标地层的慢度差。
可选的,在本发明一实施例中,所述孔隙度计算模块中孔隙度的计算公式为:
式中,Φt为孔隙度,Δtma、Δtf、Δtsh别为岩石骨架声波水平方向时差、地层流体声波水平方向时差、地层流体声波水平方向时差,Vsh为地层泥质含量,Cp为声波压实校正系数,Δt为所述目的地层的水平时差。
本发明提出的获取地层孔隙度方法及装置,通过数值模拟法建立实验地层时差各向异性比条件下,声波测井响应与井斜角的关系,再经过拟合和计算获得目标地层的水平时差,通过水平时差计算得到的地层孔隙度的孔隙度。由水平时差计算得到的地层孔隙度更加接近于地层的实际孔隙度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的获取地层孔隙度的方法流程图;
图2为本发明实施例提供的拟合方程拟合过程的方法流程图;
图3为本发明实施例提供的某工区的井斜角及慢度差的坐标示意图;
图4为本发明实施例提供的获取地层孔隙度的装置结构示意图;
图5为本发明实施例提供的拟合装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
地层形成是逐层沉积的,因此水平方向的时差最能反映地层的真实情况。因此,针对现有技术中根据声波测井仪器直接的测得的时差计算得到的孔隙度误差较大的问题,本发明提出的获取地层孔隙度方法,通过数值模拟法建立实验地层时差各向异性比条件下,声波测井响应与井斜角的关系,在经过拟合和计算获得目标地层的水平时差,通过水平时差计算得到的地层孔隙度更接近真实情况的孔隙度。
如图1的方法流程图所示,获取地层孔隙度的方法步骤如下:
步骤S101:获取目标地层的声波测井资料、井斜角曲线资料以及层位时差各向异性比。
声波测井资料包括声波在测井中的传播速度,而传播速度的倒数就是声波在测井中的传播时差。所述声波测井资料可以通过声波测井仪器获得。所述井斜角是指油井中某点的铅垂线和地球铅垂线的夹角,所述井斜角曲线资料反映了目标地层的井眼轨迹的斜度。所述时差各向异性比是指目的地层垂直方向时差与水平方向时差之比值。
步骤S102:根据所述目标地层的井斜角曲线资料以及层位时差各向异性比,按照拟合方程计算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差,所述拟合方程是根据已知的实验地层训练得到的。
所述拟合方程根据已知的实验地层训练得到的,所述拟合方程表示了在不同层位时差各向异性比条件下的井斜角和慢度差的关系,所以在已知所述目标地层的井斜角曲线资料以及层位时差各向异性比的条件下,可以获取所述目标地层的慢度差。
下面介绍所述拟合方程的拟合过程,如图2所示,所述拟合过程的流程图如下所示:
步骤S201:获取实验地层的水平时差,计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差。
首先选定一个已知的实验地层,并获取所述实验地层的水平时差(即水平方向时差)。计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差。具体地,利用声波测井仪测量并采用有限差分法计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差,其中,所述不同井斜角可以从0到90°范围内按15°的步长递增,所述不同时差各向异性比可以从1.0到1.2范围内按0.02的步长递增。
所述声波测井仪包括双发双收声波测井仪,相对于单发双收声波测井仪来说,所述双发双收声波测井仪可以消除深度误差,并具有井眼补偿的作用,使得测得的时差的精度更高。
步骤S202:可以根据所述实验地层的水平时差,计算所述实验地层的慢度差,所述慢度差可以为所述时差与所述实验地层的水平时差的差值。
步骤S203:建立所述井斜角以及所述慢度差的坐标系,在所述坐标系中插入所述不同井斜角、不同时差各向异性比条件下的所述慢度差的离散点。
建立所述井斜角以及所述慢度差的坐标系。具体地,为了方便后续的计算,可以对所述井斜角进行归一化,因此,以井斜角正弦值作为直角坐标系的横坐标,以所述慢度差作为所述直角坐标系的纵坐标。在所述坐标系中插入所述不同井斜角、不同时差各向异性比条件下的所述慢度差的离散点。
举例说明,图3是某工区的井斜角及慢度差的坐标系,其中,所述井斜角正弦值是所述坐标系的横坐标,所述慢度差是纵坐标,单位是微秒每米(μs/m),曲线上的点是在按照步长计算的不同井斜角正弦值、不同时差各向异性比条件下的慢度差值。
步骤S204:将所述离散点拟合成曲线,并获取所述曲线的拟合方程。
将插入的慢度差离散点拟合成曲线,并获取所述曲线的拟合方程,所述拟合方程可以包括一元二次方程。
如图3所示,将各离散点拟合成曲线,并获取拟合方程,获取图3中各拟合曲线的拟合方程如表1所示,其中,各向异性比为R,慢度差为y,井斜角正弦值为x:
表1 图3所示各曲线的拟合方程
时差各向异性比 | 拟合方程 |
1.04 | y=-9.8559x2+1.2641x+12.078 |
1.06 | y=-10.321x2+2.0125x+15.74 |
1.08 | y=-12.429x2+3.6068x+19.828 |
1.1 | y=-24.991x2+5.5619x+22.022 |
1.12 | y=-31.794x2+8.4162x+25.954 |
1.14 | y=-24.753x2-2.5181x+30.031 |
1.16 | y=-34.727x2+4.9642x+31.74 |
1.18 | y=-37.78x2+3.374x+35.886 |
1.2 | y=-39.303x2+3.1946x+37.953 |
步骤S204中,获取了实验地层在不同的时差各向异性比条件下,慢度差以及井斜角正弦值之间的关系,即所述拟合方程。根据步骤S101获取的目标地层的资料,按照所述拟合方程,可以计算得到所述目标地层的慢度差。具体地,,计算与所述实验地层相同的时差各向异性比条件下,所述目标地层在不同井斜角条件下的慢度差,当所述井斜角是步长15°的倍数时,直接按照所述方程求解所述目标地层的慢度差;当所述井斜角不是15°的倍数时,采用插值法计算所述目标地层的慢度差。例如,设井斜角0°时直角坐标系上的值为(x0,y0),井斜角15°时坐标系上的值为(x1,y1),则井斜角13°时值y为:
步骤S103:根据所述目标地层声波测井资料中的时差以及所述目标地层的慢度差计算所述目标地层的水平时差。
与步骤S202相同,所述水平时差为所述目标地层声波测井资料中的时差何所述目标地层的慢度差的差值。
步骤S104:根据所述目标地层的水平时差计算所述目标地层的孔隙度。
所述目标地层的孔隙度的计算公式为:
式中,Φt为孔隙度,Δtma、Δtf、Δtsh别为岩石骨架声波水平方向时差、地层流体声波水平方向时差、地层流体声波水平方向时差,Vsh为地层泥质含量,Cp为声波压实校正系数,Δt为所述目的地层的水平时差。
上式中,地层泥质含量Vsh的计算公式如下:
式中,ΔGR为自然伽马相对值,GR为当前深度点自然伽马测井值,GR砂为纯砂岩地层的自然伽马测井值,GR泥为纯泥岩的自然伽马测井值,GCUR为与地层年代有关的经验系数,老地层的GCUR=2,新地层的GCUR=3.7。
下面介绍与上述实施例方法相对应的获取地层孔隙度的装置,如图4所示,所述装置40包括:资料获取模块41,目标慢度差计算模块42,水平时差计算模块43以及孔隙度计算模块44,其中,
资料获取模块41,用于获取目标地层的声波测井资料、井斜角曲线资料以及层位时差各向异性比。
目标慢度差计算模块42,用于根据所述目标地层的井斜角曲线资料以及层位时差各向异性比,按照拟合方程计算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差,所述拟合方程是根据已知的实验地层训练得到的。
如图4所示,目标慢度差计算模块42,还包括:
目标慢度差计算子模块421,用于根据所述实验地层井斜角以及慢度差的坐标系中的曲线,按照所述拟合方程,计算与所述实验地层相同的时差各向异性比条件下,所述目标地层在不同井斜角条件下的慢度差,当所述井斜角是步长15°的倍数时,直接按照所述方程求解所述目标地层的慢度差;当所述井斜角不是15°的倍数时,采用插值法计算所述目标地层的慢度差。
如图5所示,所述目标慢度差计算模块42还包括拟合装置50,所述拟合装置50,包括:实验时差计算模块51,慢度差计算模块52,坐标系建立模块53,拟合模块54,其中,
实验时差计算模块51,用于获取实验地层的水平时差,计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差。
如图5所示,所述实验时差计算模块51还包括:
时差计算子模块511,用于利用声波测井仪采用有限差分法计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差,其中,所述不同井斜角在0到90°范围内按15°步长递增,所述不同时差各向异性比在1.0到1.2范围内按0.02步长递增。所述声波测井仪包括双发双收声波测井仪。
慢度差计算模块52,用于根据所述实验地层的水平时差,计算所述实验地层的慢度差,所述慢度差为所述时差与所述实验地层的水平时差的差值。
坐标系建立模块53,用于建立所述井斜角以及所述慢度差的坐标系,在所述坐标系中插入所述不同井斜角、不同时差各向异性比条件下的所述慢度差的离散点。
如图5所示,所述坐标系建立模块53,还包括:
坐标系建立子模块531,用于建立所述井斜角正弦值以及所述慢度差的直角坐标系。所述拟合方程可以包括一元二次方程。
拟合模块54,用于将所述离散点拟合成曲线,并获取所述曲线的拟合方程。
水平时差计算模块43,用于根据所述目标地层声波测井资料中的时差以及所述目标地层的慢度差计算所述目标地层的水平时差。
孔隙度计算模块44,用于根据所述目标地层的水平时差计算所述目标地层的孔隙度。
所述目标地层的孔隙度的计算公式为:
式中,Φt为孔隙度,Δtma、Δtf、Δtsh别为岩石骨架声波水平方向时差、地层流体声波水平方向时差、地层流体声波水平方向时差,Vsh为地层泥质含量,Cp为声波压实校正系数,Δt为所述目的地层的水平时差。
本发明实施例介绍的上述技术方案具有如下有益效果:本发明提出的获取地层孔隙度方法及装置,通过数值模拟法建立实验地层时差各向异性比条件下,声波测井响应与井斜角的关系,再经过拟合和计算获得目标地层的水平时差,通过水平时差计算得到的地层孔隙度的孔隙度。由水平时差计算得到的地层孔隙度更加接近于地层的实际孔隙度。
本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrative logical block),单元,和步骤可以通过电子硬件、电脑软件,或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability),上述的各种说明性部件(illustrative components),单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明实施例中所描述的各种说明性的逻辑块,或单元都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列(FPGA)或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。
本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息。可选地,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中。可选地,处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中,本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (14)
1.一种获取地层孔隙度的方法,其特征在于,所述方法包括:
获取目标地层的声波测井资料、井斜角曲线资料以及层位时差各向异性比;
根据所述目标地层的井斜角曲线资料以及层位时差各向异性比,按照拟合方程计算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差,所述拟合方程根据已知的实验地层按照下述方式训练得到:
获取所述实验地层的水平时差,计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差;
根据所述实验地层的水平时差,计算所述实验地层的慢度差,所述慢度差为所述时差与所述实验地层的水平时差的差值;
建立所述井斜角以及所述慢度差的坐标系,在所述坐标系中插入所述不同井斜角、不同时差各向异性比条件下的所述慢度差的离散点;
将所述离散点拟合成曲线,并获取所述曲线的拟合方程;根据所述目标地层声波测井资料中的时差以及所述目标地层的慢度差计算所述目标地层的水平时差;
根据所述目标地层的水平时差计算所述目标地层的孔隙度。
2.根据权利要求1所述的获取地层孔隙度的方法,其特征在于,所述计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差,包括:
利用声波测井仪采用有限差分法计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差,其中,所述井斜角从0到90°范围内按15°步长递增,所述时差各向异性比从1.0到1.2范围内按0.02步长递增。
3.根据权利要求2所述的获取地层孔隙度的方法,其特征在于,所述声波测井仪包括双发双收声波测井仪。
4.根据权利要求1所述的获取地层孔隙度的方法,其特征在于,所述建立所述井斜角以及所述慢度差的坐标系,具体包括:
建立所述井斜角正弦值以及所述慢度差的直角坐标系。
5.根据权利要求1所述的获取地层孔隙度的方法,其特征在于,所述拟合方程包括一元二次方程。
6.根据权利要求1所述的获取地层孔隙度的方法,其特征在于,所述按照所述拟合方程,计算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差,具体包括:
根据所述实验地层的井斜角以及慢度差的坐标系中的曲线,按照所述拟合方程,计算与所述实验地层相同的时差各向异性比条件下,所述目标地层在不同井斜角条件下的慢度差,当所述井斜角是步长15°的倍数时,按照所述拟合方程求解所述目标地层的慢度差;当所述井斜角不是15°的倍数时,采用插值法计算所述目标地层的慢度差。
7.根据权利要求1所述的获取地层孔隙度的方法,其特征在于,根据所述目标地层的水平时差计算所述目标地层的孔隙度,其中所述孔隙度的计算公式为:
<mrow>
<msub>
<mi>&Phi;</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>*</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mfrac>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
式中,Φt为孔隙度,Δtma、Δtf、Δtsh分别为岩石骨架声波水平方向时差、地层流体声波水平方向时差、地层流体声波水平方向时差,Vsh为地层泥质含量,Cp为声波压实校正系数,Δt为所述目标地层的水平时差。
8.一种获取地层孔隙度的装置,其特征在于,包括:
资料获取模块,用于获取目标地层的声波测井资料、井斜角曲线资料以及层位时差各向异性比;
目标慢度差计算模块,用于根据所述目标地层的井斜角曲线资料以及层位时差各向异性比,按照拟合方程计算所述目标地层在不同井斜角、不同层位时差各向异性比条件下的慢度差,所述拟合方程根据已知的实验地层按照下述方式训练得到:
获取所述实验地层的水平时差,计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差;
根据所述实验地层的水平时差,计算所述实验地层的慢度差,所述慢度差为所述时差与所述实验地层的水平时差的差值;
建立所述井斜角以及所述慢度差的坐标系,在所述坐标系中插入所述不同井斜角、不同时差各向异性比条件下的所述慢度差的离散点;
将所述离散点拟合成曲线,并获取所述曲线的拟合方程;水平时差计算模块,用于根据所述目标地层声波测井资料中的时差以及所述目标地层的慢度差计算所述目标地层的水平时差;
孔隙度计算模块,用于根据所述目标地层的水平时差计算所述目标地层的孔隙度。
9.根据权利要求8所述的获取地层孔隙度的装置,其特征在于,所述实验时差计算模块,还包括:
时差计算子模块,用于利用声波测井仪采用有限差分法计算所述实验地层在不同井斜角、不同时差各向异性比条件下的时差,其中,所述井斜角从0到90°范围内按15°步长递增,所述时差各向异性比从1.0到1.2范围内按0.02步长递增。
10.根据权利要求9所述的获取地层孔隙度的装置,其特征在于,所述声波测井仪包括双发双收声波测井仪。
11.根据权利要求8所述的获取地层孔隙度的装置,其特征在于,所述坐标系建立模块,还包括:
坐标系建立子模块,用于建立所述井斜角正弦值以及所述慢度差的直角坐标系。
12.根据权利要求8所述的获取地层孔隙度的装置,其特征在于,所述拟合方程包括一元二次方程。
13.根据权利要求8所述的获取地层孔隙度的装置,其特征在于,所述目标慢度差计算模块,还包括:
目标慢度差计算子模块,用于根据所述实验地层的井斜角以及慢度差的坐标系中的曲线,按照所述拟合方程,计算与所述实验地层相同的时差各向异性比条件下,所述目标地层在不同井斜角条件下的慢度差,当所述井斜角是步长15°的倍数时,按照所述拟合方程求解所述目标地层的慢度差;当所述井斜角不是15°的倍数时,采用插值法计算所述目标地层的慢度差。
14.根据权利要求8所述的获取地层孔隙度的装置,其特征在于,
所述孔隙度计算模块中孔隙度的计算公式为:
<mrow>
<msub>
<mi>&Phi;</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>*</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mfrac>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
式中,Φt为孔隙度,Δtma、Δtf、Δtsh别为岩石骨架声波水平方向时差、地层流体声波水平方向时差、地层流体声波水平方向时差,Vsh为地层泥质含量,Cp为声波压实校正系数,Δt为所述目标地层的水平时差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510368565.XA CN104975853B (zh) | 2015-06-29 | 2015-06-29 | 一种获取地层孔隙度的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510368565.XA CN104975853B (zh) | 2015-06-29 | 2015-06-29 | 一种获取地层孔隙度的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104975853A CN104975853A (zh) | 2015-10-14 |
CN104975853B true CN104975853B (zh) | 2018-01-05 |
Family
ID=54272769
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510368565.XA Active CN104975853B (zh) | 2015-06-29 | 2015-06-29 | 一种获取地层孔隙度的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104975853B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108150160B (zh) * | 2017-12-21 | 2021-07-06 | 西安石油大学 | 一种地层中欠压实作用超压的求取方法 |
CN110244385B (zh) * | 2018-03-07 | 2020-10-23 | 中国石油化工股份有限公司 | 一种大斜度井测井曲线的校正方法 |
CN109162696B (zh) * | 2018-10-08 | 2022-06-03 | 陕西延长石油(集团)有限责任公司研究院 | 一种利用补偿声波计算地层含油饱和度的方法 |
CN109339778B (zh) * | 2018-11-12 | 2021-11-16 | 中国石油大学(华东) | 一种定量评价射孔穿透深度的声学测井方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6925031B2 (en) * | 2001-12-13 | 2005-08-02 | Baker Hughes Incorporated | Method of using electrical and acoustic anisotropy measurements for fracture identification |
US7529150B2 (en) * | 2006-02-06 | 2009-05-05 | Precision Energy Services, Ltd. | Borehole apparatus and methods for simultaneous multimode excitation and reception to determine elastic wave velocities, elastic modulii, degree of anisotropy and elastic symmetry configurations |
CA2878470C (en) * | 2012-07-10 | 2021-04-27 | Statoil Petroleum As | Anisotropy parameter estimation |
CN103984022B (zh) * | 2014-04-16 | 2015-12-02 | 刘致水 | 大斜度井纵、横波测井速度各向异性同时校正方法 |
CN104632202B (zh) * | 2014-12-25 | 2017-06-13 | 中国石油天然气股份有限公司 | 确定干粘土三孔隙度测井参数值的方法及装置 |
-
2015
- 2015-06-29 CN CN201510368565.XA patent/CN104975853B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104975853A (zh) | 2015-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105317431B (zh) | 一种用于解释和评价水平井测井参数的方法 | |
CA2903451C (en) | Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data | |
CN104975853B (zh) | 一种获取地层孔隙度的方法及装置 | |
CN106154322B (zh) | 测井曲线校正方法和装置 | |
CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
BRPI0611629B1 (pt) | Método para determinar a permeabilidade de uma formação de sub-superfície de dados de registro sônicos e dados de registro de poço e para produzir hidrocarbonetos de uma formação de sub-superfície | |
CN103912268A (zh) | 基于toc的页岩储层含气饱和度确定方法 | |
CN105759312A (zh) | 一种储层特征法井震标定方法 | |
CN106154343B (zh) | 计算致密油储层的含油饱和度的方法 | |
CN102156297B (zh) | 基于砂岩油藏叠后地震数据的流体替换方法 | |
SA520411976B1 (ar) | تعيين الآثار الطباقية الكيميائية لخزان بخواص فيزيائية للصخر وعكس زلزالي | |
CN105089615A (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
CN106054254B (zh) | 一种测井曲线校正方法及装置 | |
CN104375183A (zh) | 一种获得断层平面封堵性的方法和装置 | |
US10495781B2 (en) | Dual porosity methodology for mineral volume calculations in source rock | |
CN105201491B (zh) | 一种确定水淹强度的方法及装置 | |
CN112147688B (zh) | 一种储层油气识别方法和装置 | |
CN107238859B (zh) | 计算虚拟井合成地震记录的方法及系统 | |
CN112115592A (zh) | 一种基于岩石物理实验的大斜度/水平井电阻率校正方法 | |
CN106960264A (zh) | 低渗透砂砾岩地层压力预测方法及装置 | |
CN110160927A (zh) | 一种确定束缚水饱和度的方法及装置 | |
Engler | A new approach to gas material balance in tight gas reservoirs | |
Walker et al. | Stochastic inversion for facies: A case study on the Schiehallion field | |
Craig et al. | Using Maps of Microseismic Events to Define Reservoir Discontinuities | |
Hadibeik et al. | Petrophysical properties of unconventional low-mobility reservoirs (shale gas and heavy oil) by using newly developed adaptive testing approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |