CN104215652A - 确定油气饱和度的方法和装置 - Google Patents
确定油气饱和度的方法和装置 Download PDFInfo
- Publication number
- CN104215652A CN104215652A CN201410415544.4A CN201410415544A CN104215652A CN 104215652 A CN104215652 A CN 104215652A CN 201410415544 A CN201410415544 A CN 201410415544A CN 104215652 A CN104215652 A CN 104215652A
- Authority
- CN
- China
- Prior art keywords
- echo
- string
- sounding
- reservoir
- short
- 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
- 238000000034 method Methods 0.000 title claims abstract description 80
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 97
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 64
- 239000012530 fluid Substances 0.000 claims abstract description 49
- 238000005259 measurement Methods 0.000 claims abstract description 42
- 238000009792 diffusion process Methods 0.000 claims abstract description 41
- 238000001228 spectrum Methods 0.000 claims abstract description 32
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims description 43
- 229930195733 hydrocarbon Natural products 0.000 claims description 37
- 150000002430 hydrocarbons Chemical class 0.000 claims description 37
- 239000004215 Carbon black (E152) Substances 0.000 claims description 36
- 239000007789 gas Substances 0.000 claims description 24
- 230000005311 nuclear magnetism Effects 0.000 claims description 9
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 7
- 239000001257 hydrogen Substances 0.000 claims description 7
- 229910052739 hydrogen Inorganic materials 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 abstract 2
- 239000008398 formation water Substances 0.000 abstract 1
- 238000004611 spectroscopical analysis Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000010287 polarization Effects 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000000240 pulsed field-gradient spin echo nuclear magnetic resonance spectroscopy Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种确定油气饱和度的方法和装置,其中,该方法包括:获取长回波间隔测量回波串和短回波间隔测量回波串;根据实际地层信息,在设定的实验温度和实验压力下,确定地层水扩散系数;根据所述地层水扩散系数,求取水扩散弛豫时间;将长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间加入到短回波间隔T2谱中,求取构建回波串,根据长回波间隔测量回波串和所述构建回波串确定储层的流体性质,并计算储层的油气饱和度。本发明解决了现有的核磁共振移谱法需要很强的地区经验的问题,拓宽了核磁共振测井流体识别方法的使用范围。
Description
技术领域
本发明涉及石油测井技术领域,特别涉及一种确定油气饱和度的方法和装置。
背景技术
在石油勘探和测井技术领域,核磁共振是一种重要的储层流体识别方法,尤其是在束缚水含量较高的“低阻”油层存在时,常规电阻率流体识别方法难以准确获取储层的流体信息,核磁共振流体识别技术不受束缚水等的影响,可以准确进行流体识别,具有较好的应用前景。但是,核磁共振测井的横向弛豫时间T2的分布会受到流体性质以及孔隙结构的影响,储层中的油水核磁共振T2在一定程度内重合,增加了油水性质识别的难度。
目前,在利用核磁共振测井进行油水识别时,主要采用差谱法和移谱法进行流体识别,然而,这两种方式都存在一定的问题:
1)对于差谱法来说,采用长、短不同的等待时间(即TW)状态获得的差谱信息进行流体性质识别,该方法往往需要多种约束条件,如储层润湿性为水湿、短等待时间能将水信号完全极化以及测井资料高信噪比等条件。由于孔隙结构复杂,常常遇到短等待时间无法将孔隙中水信号完全极化的情况,此时会出现水层有差谱信号的现象,从而导致水层误解释油层致使核磁测井符合率降低。
2)对于移谱法来说,主要依靠的是经验,理论上仍然相当牵强,在实际应用中难以实现油气水的定量识别。
发明内容
本发明实施例提供了一种确定油气饱和度的方法,以达到拓宽核磁共振测井流体识别方法的使用范围的目的,该方法包括:
采用核磁共振测井仪中带有双回波间隔(Echo-time,TE)信息的采集模式进行测井,获取长回波间隔测量回波串和短回波间隔测量回波串;
对所述长回波间隔测量回波串和短回波间隔测量回波串进行反演,得到长回波间隔下的横向弛豫时间T2谱和短回波间隔下的横向弛豫时间T2谱,并求取核磁共振测井噪声信息的标准方差;
根据实际地层信息,在设定的实验温度和实验压力下,确定地层水扩散系数;
根据所述地层水扩散系数,求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间;
将长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间加入到短回波间隔T2谱中,求取构建回波串,并对所述构建回波串进行反演得到构建T2谱;
将所述长回波间隔测量回波串和所述构建回波串作差值处理,将各个回波对应的差值作求和处理;
根据求和处理的结果和所述标准方差,确定储层的流体性质;
在确定储层的流体性质为含油或者含气的情况下,计算储层的油气饱和度。
在一个实施例中,短回波间隔为0.9ms或者1.2ms,长回波间隔为3.6ms或4.8ms。
在一个实施例中,按照以下公式求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间:
其中,T2l,D表示长回波间隔下的水扩散弛豫时间,T2s,D表示短回波间隔下的水扩散弛豫时间,Dw表示地层水扩散系数,γ表示氢原子的旋磁比,TEl表示长回波间隔,TEs表示短回波间隔,G表示核磁共振测井仪的静磁场梯度。
在一个实施例中,按照以下公式求取构建回波串:
其中,表示构建回波串,Mi表示短回波间隔T2谱对应的幅度值,表示短回波间隔T2谱,t表示采集时间。
在一个实施例中,按照以下公式对将所述长回波间隔测量回波串和所述构建回波串作差值处理:
其中,ECHO(TEl)表示长回波间隔测量回波串,ΔE表示差值处理的结果;
按照以下公式对求得的回波串中各个回波对应的差值作求和处理:
其中,NE表示回波串中的回波个数。
在一个实施例中,根据求和处理的结果和所述标准方差,确定储层的流体性质,包括:
按照以下方式确定储层的流体性质:
当ΔEsum/NE>δ时,确定储层含油,当ΔEsum/NE<-δ时,确定储层含气,当-δ≤ΔEsum/NE≤δ时,确定储层为水层,δ表示核磁共振测井噪声信息的标准方差。
在一个实施例中,按照以下公式计算储层的油气饱和度:
其中,|ΔEmax|是|ΔEi|中的最大幅度值,表示核磁孔隙度,HI表示油或者气的含氢指数。
本发明实施例提供了一种确定油气饱和度的装置,以达到拓宽核磁共振测井流体识别方法的使用范围的目的,该装置包括:
测井模块,用于采用核磁共振测井仪中带有双回波间隔TE信息的采集模式进行测井,获取长回波间隔测量回波串和短回波间隔测量回波串;
反演模块,用于对所述长回波间隔测量回波串和短回波间隔测量回波串进行反演,得到长回波间隔下的横向弛豫时间T2谱和短回波间隔下的横向弛豫时间T2谱,并求取核磁共振测井噪声信息的标准方差;
扩散系数确定模块,用于根据实际地层信息,在设定的实验温度和实验压力下,确定地层水扩散系数;
弛豫时间确定模块,用于根据所述地层水扩散系数,求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间;
构建模块,用于将长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间加入到短回波间隔T2谱中,求取构建回波串,并对所述构建回波串进行反演得到构建T2谱;
差值处理模块,用于将所述长回波间隔测量回波串和所述构建回波串作差值处理,将各个回波对应的差值作求和处理;
流体性质确定模块,用于根据求和处理的结果和所述标准方差,确定储层的流体性质;
油气饱和度计算模块,用于在确定储层的流体性质为含油或者含气的情况下,计算储层的油气饱和度。
在一个实施例中,短回波间隔为0.9ms或者1.2ms,长回波间隔为3.6ms或4.8ms。
在一个实施例中,弛豫信息确定模块具体用于按照以下公式求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间:
其中,T2l,D表示长回波间隔下的水扩散弛豫时间,T2s,D表示短回波间隔下的水扩散弛豫时间,Dw表示地层水扩散系数,γ表示氢原子的旋磁比,TEl表示长回波间隔,TEs表示短回波间隔,G表示核磁共振测井仪的静磁场梯度。
在一个实施例中,所述构建模块具体用于按照以下公式求取构建回波串:
其中,表示构建回波串,Mi表示短回波间隔T2谱对应的幅度值,表示短回波间隔T2谱,t表示采集时间。
在一个实施例中,所述差值处理模块具体用于按照以下公式对将所述长回波间隔测量回波串和所述构建回波串作差值处理:
其中,ECHO(TEl)表示长回波间隔测量回波串,ΔE表示差值处理的结果;
按照以下公式对求得的回波串中各个回波对应的差值作求和处理:
其中,NE表示回波串中的回波个数。
在一个实施例中,所述流体性质确定模块具体用于按照以下方式确定储层的流体性质:
当ΔEsum/NE>δ时,确定储层含油,当ΔEsum/NE<-δ时,确定储层含气,当-δ≤ΔEsum/NE≤δ时,确定储层为水层,δ表示核磁共振测井噪声信息的标准方差。
在一个实施例中,所述油气饱和度计算模块具体用于按照以下公式计算储层的油气饱和度:
其中,|ΔEmax|是|ΔEi|中的最大幅度值,表示核磁孔隙度,HI表示油或者气的含氢指数。
在本发明实施例中,将在实验室中实际测得的地层水扩散系数加入到了短回波间隔T2谱中求取构建回波串,并得到回波串反演求取构建T2谱,避免了现有的核磁共振移谱法需要很强的地区经验的问题,拓宽了核磁共振测井流体识别方法的使用范围。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是本发明实施例的确定油气饱和度的方法流程图;
图2是本发明实施例的油层的核磁共振构建回波串饱和度求取实例示意图;
图3是本发明实施例的气层的核磁共振构建回波串饱和度求取实例示意图;
图4是本发明实施例的水层的核磁共振构建回波串饱和度求取实例示意图;
图5是本发明实施例的确定油气饱和度的装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
在本例中,提供了一种确定油气饱和度的方法,在该方法中,主要是基于核磁共振侧近构建回波串来就确定油气饱和度的,以解决现有技术中核磁共振测井流体识别方法中存在的问题,主要有:差谱法在大孔径储层无法确定饱和度,移谱法定量识别难度大的问题,本例所提供的确定油气饱和度的方法,如图1所示,包括:
步骤101:采用核磁共振测井仪中带有双TE信息的采集模式进行测井,获取长回波间隔测量回波串和短回波间隔测量回波串;
步骤102:对所述长回波间隔测量回波串和短回波间隔测量回波串进行反演,得到长回波间隔T2谱和短回波间隔T2谱,并求取核磁共振测井噪声信息的标准方差;
步骤103:根据实际地层信息,在设定的实验温度和实验压力下,确定地层水扩散系数;
步骤104:根据所述地层水扩散系数,求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间;
步骤105:将所述长回波间隔下的水扩散弛豫时间和所述短回波间隔下的水扩散弛豫时间加入到短回波间隔T2谱中,求取构建回波串,并对所述构建回波串进行反演得到构建T2谱;
步骤106:将所述长回波间隔测量回波串和所述构建回波串作差值处理,将各个回波对应的差值作求和处理;
步骤107:根据求和处理的结果和所述标准方差,确定储层的流体性质;
步骤108:在确定储层的流体性质为含油或者含气的情况下,计算储层的油气饱和度。
在上述实施例中,将在实验室中实际测得的地层水扩散系数加入到了短回波间隔T2谱中求取回波串,并得到回波串反演求取构建T2谱,避免了现有的核磁共振移谱法需要很强的地区经验的问题,拓宽了核磁共振测井流体识别方法的使用范围。进一步的,通过上述方式,解决了现有技术中差谱法在大孔径储层中无法确定饱和度,移谱法定量识别难度大的技术问题,达到了不仅可以有效确定大孔径储层的饱和度,还可以进行简单的定量识别的目的,拓宽了现有核磁共振测井流体识别方法。
具体的,在实施的过程中,核磁共振采集模式的长等待时间通常为10到13秒,以确保所测量的储层流体完全极化,短回波间隔通常0.9ms或者1.2ms,以确保孔隙度等储层信息的准确求取,长回波间隔通常为3.6ms或4.8ms。上述的核磁共振测井仪可是以MRIL-Prime核磁共振测井仪,因此可以采用D9TWE3、D9TWE4、D9TW312、D9TW412等包含双TE的测井模式。
具体的,可以按照以下公式求取标准方差:
其中,Xi表示第i个噪声信号,表示噪声平均值,n表示要处理的噪声的个数。
在上述步骤104中,可以按照以下公式求取所述核磁共振测井仪在长回波间隔下和短回波间隔下的水扩散弛豫信息:
其中,T2l,D表示长回波间隔下的水扩散弛豫时间,T2s,D表示短回波间隔下的水扩散弛豫时间,Dw表示地层水扩散系数,γ表示氢原子的旋磁比,TEl表示长回波间隔,TEs表示短回波间隔,G表示核磁共振测井仪的静磁场梯度。
在上述步骤106中,可以按照以下公式求取构建回波串:
其中,表示构建回波串,Mi表示短回波间隔T2谱对应的幅度值,表示短回波间隔T2谱,t表示采集时间。
具体实施时,可以先按照以下公式对将所述长回波间隔测量回波串和所述构建回波串作差值处理:
其中,ECHO(TEl)表示长回波间隔测量回波串,ΔE表示差值处理的结果;
然后,按照以下公式对求得的回波串中各个回波对应的差值作求和处理:
其中,NE表示回波串中的回波个数。
在上述步骤108中,根据求和处理的结果和所述标准方差,确定储层的流体性质所依据的原理如下:
1)当ΔEsum/NE>δ时,确定储层含油;
2)当ΔEsum/NE<-δ时,确定储层含气;
3)当-δ≤ΔEsum/NE≤δ时,确定储层为水层,不含油气,其中,δ表示核磁共振测井噪声信息的标准方差。
具体的,最终可以按照以下公式计算储层的油气饱和度:
其中,|ΔEmax|是|ΔEi|中的最大幅度值,表示核磁孔隙度,HI表示油或者气的含氢指数。
下面结构一个具体的实施例进行说明,然而值得注意的是该具体实施例仅是为了更好地说明本发明,并不构成对本发明的不当限定。
该储层油气饱和度的确定方法,包括:
步骤1:利用带有双TE信息的采集模式进行测井,获取不同回波间隔(TE)的两组回波串信息,并采用多指数反演算法获取各自的T2谱,并求取核磁测井噪声信息的标准方差δ;
例如:对D9TWE3模式的测井资料数据进行拆分,获取双TE资料(其中,两组回波等待时间均为12.988s,回波间隔分别为0.9ms和3.6ms),将双TE资料进行多指数反演,根据反演的结果求取噪声道的标准方差;
步骤2:将测井对应储层的地层水装入无磁夹持器中,根据实际地层信息,设计相应的实验温度和压力,利用频率为2MHz的核磁共振岩心分析仪,并采用PFGSE脉冲序列确定地层水的扩散系数D,并求取长短回波间隔下核磁共振测井仪器的水扩散弛豫信息;
具体的,将现场采集的地层水装入无磁夹持器的流体容器中并加温加压,利用频率为2MHz的核磁共振岩心分析仪,并采用PFGSE脉冲序列的多组回波数据,利用线性拟合方法确定高温高压下地层水的扩散系数D。
步骤3:将长、短回波时间的扩散弛豫信息以及短回波间隔(TEs)对应到T2谱中,然后求取构建回波串,并利用多指数算法确定构建T2谱。
具体的,将0.9ms和3.6ms回波时间的扩散弛豫信息以及短回波间隔(TEs)对应到T2谱中,采用如下公式求取长回波间隔(TEl)的回波串信息,此时的回波串成为构建回波串,并利用多指数反演算法确定构建T2谱;
步骤4:将长回波间隔3.6ms的测量回波串和构建回波串作差值处理,得到并将每个采集时间的差值求和(其值为ΔEsum),当ΔEsum/NE>δ时(此时回波个数NE为500,δ表示标准方差),认为储层含油信号,并将ΔE进行多指数反演求取油谱信息;当ΔEsum/NE<-δ时,认为储层含气信息,并将-ΔE进行多指数反演求取气谱信息;当-δ≤ΔEsum/NE≤δ,认为储层为水层。
步骤5:利用以下关系求取油气饱和度:
其中,|ΔEmax|为|ΔEi|回波串中的最大幅度值,为核磁孔隙度,其值为用短回波间隔所得的T2谱求和所得,HI为油气含氢指数,利用实验数据确定,油的HI为0.8,气的HI为0.3。
如图2所示是油层的核磁共振构建回波串饱和度求取实例示意图,如图3所示是气层的核磁共振构建回波串饱和度求取实例示意图,如图4所示是水层的核磁共振构建回波串饱和度求取实例示意图,在图2至4中,第一道为岩性曲线,包括自然伽马(GR)、自然电位(SP)和井径曲线,第二道为深度道,单位为m,第三道为电阻率曲线,包括深、中、浅电阻率,第四道为孔隙度曲线道,包括中子、密度和声波资料,第五道为长等待时间(12.988s)短回波间隔(0.9ms)测量横向弛豫时间T2谱,第六道长等待时间(12.988s)长回波间隔(3.6ms)测量横向弛豫时间T2谱,第七道包括构建谱和油谱结果,第八道为油气饱和度和流体剖面,包括束缚水孔隙度、可动水孔隙度和油气孔隙度。
在本例中提出了一种确定油气饱和度的方法,通过该方法将实验室中确定的地层水扩散系数加入到了测井确定的T2谱中,然后根据加入地层水扩散系数的T2谱重构了回波串和T2谱,通过比对构建回波串和长回波间隔的回波串之间的差异来确定储层中的流体信息以及饱和度,避免了现有的核磁共振移谱法需要很强的地区经验的问题,拓宽了现有核磁共振测井流体识别方法。
基于同一发明构思,本发明实施例中还提供了一种确定油气饱和度的装置,如下面的实施例所述。由于确定油气饱和度的装置解决问题的原理与确定油气饱和度的方法相似,因此确定油气饱和度的装置的实施可以参见确定油气饱和度的方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。图5是本发明实施例的确定油气饱和度的装置的一种结构框图,如图5所示,包括:测井模块501、反演模块502、扩散系数确定模块503、弛豫信息确定模块504、构建模块505、差值处理模块506、流体性质确定模块507和油气饱和度计算模块508,下面对该结构进行说明。
测井模块501,用于采用核磁共振测井仪中带有双TE信息的采集模式进行测井,获取长回波间隔测量回波串和短回波间隔测量回波串;
反演模块502,用于对所述长回波间隔测量回波串和短回波间隔测量回波串进行反演,得到长回波间隔T2谱和短回波间隔T2谱,并求取核磁共振测井噪声信息的标准方差;
扩散系数确定模块503,用于根据实际地层信息,在设定的实验温度和实验压力下,确定地层水扩散系数;
弛豫时间确定模块504,用于根据所述地层水扩散系数,求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间;
构建模块505,用于将长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间加入到短回波间隔T2谱中,求取构建回波串,并对所述构建回波串进行反演得到构建T2谱;
差值处理模块506,用于将所述长回波间隔测量回波串和所述构建回波串作差值处理,将各个回波对应的差值作求和处理;
流体性质确定模块507,用于根据求和处理的结果和所述标准方差,确定储层的流体性质;
油气饱和度计算模块508,用于在确定储层的流体性质为含油或者含气的情况下,计算储层的油气饱和度。
在一个实施例中,短回波间隔为0.9ms或者1.2ms,长回波间隔为3.6ms或4.8ms。
在一个实施例中,弛豫时间确定模块504具体可以按照以下公式求取所述核磁共振测井仪在长回波间隔下和短回波间隔下的水扩散弛豫时间:
其中,T2l,D表示长回波间隔下的水扩散弛豫信息,T2s,D表示短回波间隔下的水扩散弛豫信息,Dw表示地层水扩散系数,γ表示氢原子的旋磁比,TEl表示长回波间隔,TEs表示短回波间隔,G表示核磁共振测井仪的静磁场梯度。
在一个实施例中,构建模块206具体可以按照以下公式求取构建回波串:
其中,表示构建回波串,Mi表示短回波间隔T2谱对应的幅度值,表示短回波间隔T2谱,t表示采集时间。
在一个实施例中,所述差值处理模块具体用于按照以下公式对将所述长回波间隔测量回波串和所述构建回波串作差值处理:
其中,ECHO(TEl)表示长回波间隔测量回波串,ΔE表示差值处理的结果;
按照以下公式对求得的回波串中各个回波对应的差值作求和处理:
其中,NE表示回波串中的回波个数。
在一个实施例中,流体性质确定模块507具体用于按照以下方式确定储层的流体性质:
当ΔEsum/NE>δ时,确定储层含油,当ΔEsum/NE<-δ时,确定储层含气,当-δ≤ΔEsum/NE≤δ时,确定储层为水层,δ表示核磁共振测井噪声信息的标准方差。
在一个实施例中,油气饱和度计算模块508具体用于按照以下公式计算储层的油气饱和度:
其中,|ΔEmax|是|ΔEi|中的最大幅度值,表示核磁孔隙度,HI表示油或者气的含氢指数。
在另外一个实施例中,还提供了一种软件,该软件用于执行上述实施例及优选实施方式中描述的技术方案。
在另外一个实施例中,还提供了一种存储介质,该存储介质中存储有上述软件,该存储介质包括但不限于:光盘、软盘、硬盘、可擦写存储器等。
从以上的描述中,可以看出,本发明实施例实现了如下技术效果:在对得到的长回波间隔测量回波串和短回波间隔测量回波串进行反演得到长回波间隔T2谱和短回波间隔T2谱之后,根据实际的地层信息确定地层水扩散系数,然后将地层水扩散系数加入到了短回波间隔T2谱中求取回波串,并得到回波串反演求取构建T2谱,再对长回波间隔测量回波串和构建回波串作差值处理,根据差值处理后的结果确定储层的流体性质,并进行储层的油气饱和度分析。通过上述方式,解决了现有技术中差谱法在大孔径储层中无法确定饱和度,移谱法定量识别难度大的技术问题,达到了不仅可以有效确定大孔径储层的饱和度,还可以进行简单的定量识别的目的,拓宽了现有核磁共振测井流体识别方法。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (14)
1.一种确定油气饱和度的方法,其特征在于,包括:
采用核磁共振测井仪中带有双回波间隔TE信息的采集模式进行测井,获取长回波间隔测量回波串和短回波间隔测量回波串;
对所述长回波间隔测量回波串和短回波间隔测量回波串进行反演,得到长回波间隔下的横向弛豫时间T2谱和短回波间隔下的横向弛豫时间T2谱,并求取核磁共振测井噪声信息的标准方差;
根据实际地层信息,在设定的实验温度和实验压力下,确定地层水扩散系数;
根据所述地层水扩散系数,求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间;
将长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间加入到短回波间隔T2谱中,求取构建回波串,并对所述构建回波串进行反演得到构建T2谱;
将所述长回波间隔测量回波串和所述构建回波串作差值处理,将各个回波对应的差值作求和处理;
根据求和处理的结果和所述标准方差,确定储层的流体性质;
在确定储层的流体性质为含油或者含气的情况下,计算储层的油气饱和度。
2.如权利要求1所述的方法,其特征在于,短回波间隔为0.9ms或者1.2ms,长回波间隔为3.6ms或4.8ms。
3.如权利要求1所述的方法,其特征在于,按照以下公式求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间:
其中,T2l,D表示长回波间隔下的水扩散弛豫时间,T2s,D表示短回波间隔下的水扩散弛豫时间,Dw表示地层水扩散系数,γ表示氢原子的旋磁比,TEl表示长回波间隔,TEs表示短回波间隔,G表示核磁共振测井仪的静磁场梯度。
4.如权利要求1所述的方法,其特征在于,按照以下公式求取构建回波串:
其中,表示构建回波串,Mi表示短回波间隔T2谱对应的幅度值,表示短回波间隔T2谱,t表示采集时间。
5.如权利要求4所述的方法,其特征在于:
按照以下公式对将所述长回波间隔测量回波串和所述构建回波串作差值处理:
其中,ECHO(TEl)表示长回波间隔测量回波串,ΔE表示差值处理的结果;
按照以下公式对求得的回波串中各个回波对应的差值作求和处理:
其中,NE表示回波串中的回波个数。
6.如权利要求5所述的方法,其特征在于,根据求和处理的结果和所述标准方差,确定储层的流体性质,包括:
按照以下方式确定储层的流体性质:
当ΔEsum/NE>δ时,确定储层含油,当ΔEsum/NE<-δ时,确定储层含气,当-δ≤ΔEsum/NE≤δ时,确定储层为水层,δ表示核磁共振测井噪声信息的标准方差。
7.如权利要求5所述的方法,其特征在于,按照以下公式计算储层的油气饱和度:
其中,|ΔEmax|是|ΔEi|中的最大幅度值,表示核磁孔隙度,HI表示油或者气的含氢指数。
8.一种确定油气饱和度的装置,其特征在于,包括:
测井模块,用于采用核磁共振测井仪中带有双回波间隔TE信息的采集模式进行测井,获取长回波间隔测量回波串和短回波间隔测量回波串;
反演模块,用于对所述长回波间隔测量回波串和短回波间隔测量回波串进行反演,得到长回波间隔下的横向弛豫时间T2谱和短回波间隔下的横向弛豫时间T2谱,并求取核磁共振测井噪声信息的标准方差;
扩散系数确定模块,用于根据实际地层信息,在设定的实验温度和实验压力下,确定地层水扩散系数;
弛豫时间确定模块,用于根据所述地层水扩散系数,求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间;
构建模块,用于将长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间加入到短回波间隔T2谱中,求取构建回波串,并对所述构建回波串进行反演得到构建T2谱;
差值处理模块,用于将所述长回波间隔测量回波串和所述构建回波串作差值处理,将各个回波对应的差值作求和处理;
流体性质确定模块,用于根据求和处理的结果和所述标准方差,确定储层的流体性质;
油气饱和度计算模块,用于在确定储层的流体性质为含油或者含气的情况下,计算储层的油气饱和度。
9.如权利要求8所述的装置,其特征在于,短回波间隔为0.9ms或者1.2ms,长回波间隔为3.6ms或4.8ms。
10.如权利要求8所述的装置,其特征在于,弛豫信息确定模块具体用于按照以下公式求取所述核磁共振测井仪在长回波间隔下的水扩散弛豫时间和短回波间隔下的水扩散弛豫时间:
其中,T2l,D表示长回波间隔下的水扩散弛豫时间,T2s,D表示短回波间隔下的水扩散弛豫时间,Dw表示地层水扩散系数,γ表示氢原子的旋磁比,TEl表示长回波间隔,TEs表示短回波间隔,G表示核磁共振测井仪的静磁场梯度。
11.如权利要求8所述的装置,其特征在于,所述构建模块具体用于按照以下公式求取构建回波串:
其中,表示构建回波串,Mi表示短回波间隔T2谱对应的幅度值,表示短回波间隔T2谱,t表示采集时间。
12.如权利要求11所述的装置,其特征在于,所述差值处理模块具体用于按照以下公式对将所述长回波间隔测量回波串和所述构建回波串作差值处理:
其中,ECHO(TEl)表示长回波间隔测量回波串,ΔE表示差值处理的结果;
按照以下公式对求得的回波串中各个回波对应的差值作求和处理:
其中,NE表示回波串中的回波个数。
13.如权利要求12所述的装置,其特征在于,所述流体性质确定模块具体用于按照以下方式确定储层的流体性质:
当ΔEsum/NE>δ时,确定储层含油,当ΔEsum/NE<-δ时,确定储层含气,当-δ≤ΔEsum/NE≤δ时,确定储层为水层,δ表示核磁共振测井噪声信息的标准方差。
14.如权利要求12所述的装置,其特征在于,所述油气饱和度计算模块具体用于按照以下公式计算储层的油气饱和度:
其中,|ΔEmax|是|ΔEi|中的最大幅度值,表示核磁孔隙度,HI表示油或者气的含氢指数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410415544.4A CN104215652B (zh) | 2014-08-21 | 2014-08-21 | 确定油气饱和度的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410415544.4A CN104215652B (zh) | 2014-08-21 | 2014-08-21 | 确定油气饱和度的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104215652A true CN104215652A (zh) | 2014-12-17 |
CN104215652B CN104215652B (zh) | 2016-08-31 |
Family
ID=52097370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410415544.4A Active CN104215652B (zh) | 2014-08-21 | 2014-08-21 | 确定油气饱和度的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104215652B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104990854A (zh) * | 2015-07-06 | 2015-10-21 | 中国石油天然气股份有限公司 | 确定束缚水饱和度的方法及装置 |
CN105464654A (zh) * | 2015-12-15 | 2016-04-06 | 中国石油天然气股份有限公司 | 确定致密砂岩饱和度指数的方法及装置 |
CN105891250A (zh) * | 2016-03-30 | 2016-08-24 | 中国石油天然气股份有限公司 | 致密砂岩储层原始含水饱和度的确定方法 |
CN106285661A (zh) * | 2016-08-24 | 2017-01-04 | 中国石油天然气股份有限公司 | 一种基于判定指数的低阻油层识别方法及装置 |
CN106644879A (zh) * | 2016-11-17 | 2017-05-10 | 中国石油天然气股份有限公司 | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 |
CN108474250A (zh) * | 2016-03-04 | 2018-08-31 | 哈利伯顿能源服务公司 | 用于确定地下地层的孔隙度和孔隙类型的多深度研究核磁共振测井 |
CN109594971A (zh) * | 2018-12-21 | 2019-04-09 | 中国石油天然气集团有限公司 | 基于核磁共振测井增强扩散气层识别因子的流体性质识别方法 |
CN109779619A (zh) * | 2019-01-21 | 2019-05-21 | 中国石油天然气集团有限公司 | 一种通过核磁t2谱计算地层含水饱和度的方法 |
CN109856688A (zh) * | 2019-01-31 | 2019-06-07 | 中国石油天然气集团有限公司 | 基于核磁测井双tw极化增强法的流体性质识别方法 |
CN110133553A (zh) * | 2019-05-10 | 2019-08-16 | 浙江大学 | 一种超短回波时间磁共振指纹弛豫时间测量方法 |
CN113125485A (zh) * | 2019-12-30 | 2021-07-16 | 中国石油天然气股份有限公司 | 柱塞岩样含油饱和度无损测量方法及装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5696448A (en) * | 1995-06-26 | 1997-12-09 | Numar Corporation | NMR system and method for formation evaluation using diffusion and relaxation log measurements |
JPH09203715A (ja) * | 1996-01-24 | 1997-08-05 | Kobe Steel Ltd | 核磁気共鳴緩和時間の測定方法 |
CN101943669B (zh) * | 2010-09-01 | 2012-02-22 | 苏州纽迈电子科技有限公司 | 低场核磁共振测量钻井液含油率的方法 |
CN102998322B (zh) * | 2011-09-14 | 2014-08-06 | 中国石油天然气股份有限公司 | 恒定梯度场核磁共振岩样分析方法及仪器 |
CN102508293B (zh) * | 2011-11-28 | 2013-05-15 | 中国石油大学(北京) | 一种叠前反演的薄层含油气性识别方法 |
CN103675907A (zh) * | 2012-09-20 | 2014-03-26 | 中国石油化工股份有限公司 | 一种基于岩相约束的avo反演油气检测方法 |
CN103760081B (zh) * | 2013-12-31 | 2016-01-06 | 中国石油天然气股份有限公司 | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及*** |
-
2014
- 2014-08-21 CN CN201410415544.4A patent/CN104215652B/zh active Active
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104990854A (zh) * | 2015-07-06 | 2015-10-21 | 中国石油天然气股份有限公司 | 确定束缚水饱和度的方法及装置 |
CN105464654A (zh) * | 2015-12-15 | 2016-04-06 | 中国石油天然气股份有限公司 | 确定致密砂岩饱和度指数的方法及装置 |
CN105464654B (zh) * | 2015-12-15 | 2018-06-01 | 中国石油天然气股份有限公司 | 确定致密砂岩饱和度指数的方法及装置 |
CN108474250A (zh) * | 2016-03-04 | 2018-08-31 | 哈利伯顿能源服务公司 | 用于确定地下地层的孔隙度和孔隙类型的多深度研究核磁共振测井 |
CN105891250A (zh) * | 2016-03-30 | 2016-08-24 | 中国石油天然气股份有限公司 | 致密砂岩储层原始含水饱和度的确定方法 |
CN106285661B (zh) * | 2016-08-24 | 2020-01-07 | 中国石油天然气股份有限公司 | 一种基于判定指数的低阻油层识别方法及装置 |
CN106285661A (zh) * | 2016-08-24 | 2017-01-04 | 中国石油天然气股份有限公司 | 一种基于判定指数的低阻油层识别方法及装置 |
CN106644879A (zh) * | 2016-11-17 | 2017-05-10 | 中国石油天然气股份有限公司 | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 |
CN106644879B (zh) * | 2016-11-17 | 2019-09-10 | 中国石油天然气股份有限公司 | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 |
CN109594971A (zh) * | 2018-12-21 | 2019-04-09 | 中国石油天然气集团有限公司 | 基于核磁共振测井增强扩散气层识别因子的流体性质识别方法 |
CN109779619A (zh) * | 2019-01-21 | 2019-05-21 | 中国石油天然气集团有限公司 | 一种通过核磁t2谱计算地层含水饱和度的方法 |
CN109779619B (zh) * | 2019-01-21 | 2023-02-28 | 中国石油天然气集团有限公司 | 一种通过核磁t2谱计算地层含水饱和度的方法 |
CN109856688A (zh) * | 2019-01-31 | 2019-06-07 | 中国石油天然气集团有限公司 | 基于核磁测井双tw极化增强法的流体性质识别方法 |
CN110133553A (zh) * | 2019-05-10 | 2019-08-16 | 浙江大学 | 一种超短回波时间磁共振指纹弛豫时间测量方法 |
CN110133553B (zh) * | 2019-05-10 | 2020-06-05 | 浙江大学 | 一种超短回波时间磁共振指纹弛豫时间测量方法 |
CN113125485A (zh) * | 2019-12-30 | 2021-07-16 | 中国石油天然气股份有限公司 | 柱塞岩样含油饱和度无损测量方法及装置 |
CN113125485B (zh) * | 2019-12-30 | 2023-10-31 | 中国石油天然气股份有限公司 | 柱塞岩样含油饱和度无损测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104215652B (zh) | 2016-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104215652A (zh) | 确定油气饱和度的方法和装置 | |
Mao et al. | Estimation of permeability by integrating nuclear magnetic resonance (NMR) logs with mercury injection capillary pressure (MICP) data in tight gas sands | |
US11650347B2 (en) | Fast measurement and interpretation of downhole multi-dimensional measurement | |
US11435304B2 (en) | Estimating downhole fluid volumes using multi-dimensional nuclear magnetic resonance measurements | |
CN104453874B (zh) | 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法 | |
CN102042011B (zh) | 利用常规测井资料构建伪核磁t2谱的方法 | |
CN101303417B (zh) | 用于分析具有共用和不同特性的数据的方法 | |
MX2008008999A (es) | Pelicula de capas multiples resistente a la intemperie. | |
US20150177351A1 (en) | Methods of investigating formation samples using nmr data | |
CN103674811A (zh) | 一种核磁共振孔隙度测量的校正方法、装置及*** | |
Heaton et al. | Saturation and viscosity from multidimensional nuclear magnetic resonance logging | |
CN102913240B (zh) | 一种储层流体识别方法 | |
CN103670390A (zh) | 一种水淹层测井评价方法及*** | |
CN106154343A (zh) | 计算致密油储层的含油饱和度的方法 | |
Gu et al. | A novel method for NMR data denoising based on discrete cosine transform and variable length windows | |
CN107808032A (zh) | 页岩气储层含水饱和度计算方法及*** | |
Anand et al. | Predicting effective permeability to oil in sandstone and carbonate reservoirs from well-logging data | |
Chen et al. | New approaches of 3D nuclear magnetic resonance inversion for improving fluid typing | |
Al-Sulami et al. | The unconventional shale reservoirs of jafurah basin: An integrated petrophysical evaluation using cores and advanced well logs | |
Ye et al. | Method of determining unconventional reservoir saturation with NMR logging | |
CN109612900A (zh) | 一种储层岩石渗透率预测方法、装置及存储介质 | |
Wu et al. | A new method for predicting capillary pressure curves based on NMR echo data: Sandstone as an example | |
Feng et al. | A new empirical method for constructing capillary pressure curves from conventional logs in low-permeability sandstones | |
Sun et al. | Apparent hydrogen index and its correlation with heavy oil viscosity | |
Bachman et al. | Practical Down Hole Dielectric and Diffusion-based Magnetic Resonance Workflow for Saturations and Viscosity of Heavy Oil Reservoirs Using a Laboratory Viscosity Calibration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |