CN104330433A - 一种获得目的储层横向弛豫时间分布的方法及装置 - Google Patents

一种获得目的储层横向弛豫时间分布的方法及装置 Download PDF

Info

Publication number
CN104330433A
CN104330433A CN201410586502.7A CN201410586502A CN104330433A CN 104330433 A CN104330433 A CN 104330433A CN 201410586502 A CN201410586502 A CN 201410586502A CN 104330433 A CN104330433 A CN 104330433A
Authority
CN
China
Prior art keywords
crude oil
intrinsic
distribution
magnetic resonance
object reservoir
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
Application number
CN201410586502.7A
Other languages
English (en)
Other versions
CN104330433B (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 CN201410586502.7A priority Critical patent/CN104330433B/zh
Publication of CN104330433A publication Critical patent/CN104330433A/zh
Application granted granted Critical
Publication of CN104330433B publication Critical patent/CN104330433B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本申请实施例提供一种获得目的储层横向弛豫时间分布的方法及装置。该方法包括:根据目的储层的温度和压力信息,设置相应的实验温度和压力,获取目的储层自由状态地层水的核磁共振本征横向弛豫时间分布、目的储层原油的核磁共振本征横向弛豫时间分布以及目的储层地层水的核磁共振本征横向弛豫时间分布;再计算出原油的含氢指数、地层水的含氢指数、原油的扩散系数和地层水的扩散系数,接着,根据这些数据将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号,再将自旋回波信号进行反演获得目的储层在完全含水状态下的核磁共振横向弛豫时间分布。

Description

一种获得目的储层横向弛豫时间分布的方法及装置
技术领域
本发明涉及石油测井技术领域,尤其涉及一种获得目的储层横向弛豫时间分布的方法及装置。
背景技术
在石油测井技术领域,核磁共振技术得到了广泛应用,核磁共振技术是确定储层油气孔径分布信息,实现复杂储层的定量评价的有效手段。目前利用核磁共振技术获取储层油气孔径分布信息的主要方法是将核磁共振横向弛豫时间T2分布转换为孔径分布,且必须为完全含水状态下的核磁共振横向弛豫时间T2分布才能准确获取孔径分布信息,但直接利用自由水核磁共振横向弛豫时间T2分布,就忽略了实际储层中油气会对T2分布产生的影响,使得获得的核磁共振横向弛豫时间T2分布无法准确的进行储层评价。
发明内容
本申请的目的是提供一种获得目的储层横向弛豫时间分布的方法及装置,能够去除储层中油气对核磁共振横向弛豫时间T2分布的影响。
为了实现上述目的,本发明提供了一种获得目的储层横向弛豫时间分布的方法,该方法包括:
S101:根据目的储层的温度和压力信息,设置相应的实验温度和压力,获取目的储层岩心中目的储层地层水的核磁共振本征横向弛豫时间T2分布,所述目的储层岩心中目的储层地层水是自由状态的地层水;
S102:根据目的储层温度和压力信息,设置相应的实验温度和压力,获取目的储层原油的核磁共振本征横向弛豫时间T2分布和目的储层地层水的核磁共振本征横向弛豫时间T2分布,根据所述目的储层原油的核磁共振本征横向弛豫时间T2分布确定原油的含氢指数,根据所述目的储层的地层水的核磁共振本征横向弛豫时间T2分布确定地层水的含氢指数;
S103:获取与所述目的储层的温度和压力相同的温度和压力下所述原油的扩散系数、以及获取与所述目的储层的温度和压力相同的温度和压力下所述地层水的扩散系数;
S104:利用核磁共振测井获取目的储层下包含原油、地层水和岩心的自旋回波信号,将所述原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号;所述原油的本征弛豫和扩散弛豫综合贡献是通过所述原油的本征横向弛豫时间T2分布、所述原油的含氢指数和所述原油的扩散系数确定的,所述自由状态地层水的本征弛豫和扩散弛豫综合贡献是通过所述自由状态的地层水的本征横向弛豫时间T2分布、所述地层水的含氢指数和所述地层水的扩散系数确定的;
S105:将所述目的储层在完全含水状态下的自旋回波信号进行反演,获取目的储层在完全含水状态下的核磁共振横向弛豫时间T2分布。
在一个优选的实施例中,所述自由状态的地层水的核磁共振本征横向弛豫时间T2分布是采用CPMG脉冲序列得到自旋回波信号,将自旋回波信号反演得到的。
在一个优选的实施例中,所述目的储层原油的核磁共振本征横向弛豫时间T2分布是采用CPMG脉冲序列得到自旋回波信号,将自旋回波信号反演得到的。
在一个优选的实施例中,所述目的储层地层水的核磁共振本征横向弛豫时间T2分布是采用CPMG的脉冲序列得到自旋回波信号,将自旋回波信号反演得到的。
在一个优选的实施例中,所述目的储层原油的含氢指数是利用标准水样将所述原油的核磁共振本征横向弛豫时间T2分布进行刻度,计算的公式如下:
I Ho = Σ i n P i / V 0 Σ i n P i , s / V s
上式中,为求和公式;IHo为所述原油的含氢指数;Pi为所述原油的核磁共振横向弛豫时间T2分布的第i个区间孔隙度分量;V0为所述原油体积;Pi,s为标准水样的核磁共振横向弛豫时间T2的第i个区间孔隙度分量;Vs为标准水样的体积。
在一个优选的实施例中,所述目的储层地层水的含氢指数是利用标准水样将所述地层水的核磁共振本征横向弛豫时间T2分布进行刻度,计算的公式如下:
I Hw = Σ i n P i ′ / V w Σ i n P i , s / V s
上式中,为求和公式;IHw为所述地层水的含氢指数;Pi'为所述地层水的核磁共振横向弛豫时间T2分布的第i个区间孔隙度分量;Vw为所述地层水体积;Pi,s为标准水样的核磁共振横向弛豫时间T2的第i个区间孔隙度分量;Vs为标准水样的体积。
在一个优选的实施例中,所述目的储层原油的扩散系数是在脉冲梯度磁场的条件下,采用PFGSE脉冲序列得到原油自旋回波信号,获取原油回波间隔和原油回波幅度,所述原油回波间隔和所述原油回波幅度之间的斜率是所述原油的扩散系数。
PFGSE脉冲为脉冲梯度场自旋回波脉冲。
在一个优选的实施例中,所述目的储层地层水的扩散系数是在脉冲梯度磁场的条件下,采用PFGSE脉冲序列得到地层水自旋回波信号,获取地层水回波间隔和地层水回波幅度,所述地层水回波间隔和所述地层水回波幅度之间的斜率是所述地层水的扩散系数。
在一个优选的实施例中,所述将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层条件在完全含水状态下的自旋回波信号具体通过公式获取:
ECHOW ( t ) = ECHO ( t ) - φ S o I Ho e - t [ 1 T 2 o + D o ( γG T E ) 2 12 ] + φ S o I Hw e - t [ 1 T 2 W + D w ( γG T E ) 2 12 ]
其中,ECHOW(t)为目的储层在完全含水状态下的自旋回波信号;ECHO(t)为目的储层包含原油、地层水和岩心自旋回波信号;φ为核磁共振孔隙度;So为含油饱和度;γ为氢核的旋磁比,是常数;G为磁场梯度;TE为CPMG脉冲序列的回波间隔;IHo为原油的含氢指数;IHw为地层水的含氢指数;T2o为原油的本征弛豫时间;T2w为目的储层岩心中自由状态的地层水的本征横向弛豫时间;Do为目的储层原油的扩散系数;Dw为目的储层地层水的扩散系数;t为采集时间。
本申请另一方面还提供一种获得目的储层横向弛豫时间分布的的装置,该装置包括:
第一横向弛豫时间分布获取单元,用于获取与所述目的储层的温度和压力相同的温度和压力下的目的储层岩心中目的储层地层水的核磁共振本征横向弛豫时间T2分布,所述目的储层岩心中目的储层地层水是自由状态的地层水;
第二横向弛豫时间分布获取单元,用于获取与所述目的储层的温度和压力相同的温度和压力下的目的储层原油的核磁共振本征横向弛豫时间T2分布,以及获取与所述目的储层的温度和压力相同的温度和压力下的目的储层地层水的核磁共振本征横向弛豫时间T2分布;
含氢指数获取单元,用于根据第二横向弛豫时间分布获取单元得到的所述目的储层原油的核磁共振本征横向弛豫时间T2分布确定原油的含氢指数,用于根据第二横向弛豫时间分布得到的所述目的储层的地层水的核磁共振本征横向弛豫时间T2分布确定地层水的含氢指数;
扩散系数获取单元,用于获取与所述目的储层的温度和压力相同的温度和压力下所述原油的扩散系数,以及获取与所述目的储层的温度和压力相同的温度和压力下所述地层水的扩散系数;
自旋回波信号获取单元,用于利用核磁共振测井获取目的储层下包含原油、地层水和岩心的自旋回波信号,将所述原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号;所述原油的本征弛豫和扩散弛豫综合贡献是通过所述原油的本征横向弛豫时间T2分布、所述原油的含氢指数和所述原油的扩散系数确定的,所述自由状态地层水的本征弛豫和扩散弛豫综合贡献是通过所述自由状态的地层水的本征横向弛豫时间T2分布、所述地层水的含氢指数和所述地层水的扩散系数确定的;
第三横向弛豫时间分布获取单元,用于将自旋回波信号获取单元得到的目的储层在完全含水状态下的自旋回波信号进行反演,获取目的储层在完全含水状态下的核磁共振横向弛豫时间T2分布。
由以上本申请实施例提供的技术方案可见,本申请根据目的储层的温度和压力信息,设置好实验温度和压力,获取目的储层自由状态的地层水的核磁共振横向弛豫时间T2分布、目的储层原油核磁共振横向弛豫时间T2分布、目的储层原油的含氢指数、目的储层地层水的含氢指数、目的储层原油的扩散系数和目的储层地层水的扩散系数,根据这些数据将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号,再将自旋回波信号进行反演获得目的储层在完全含水状态下的核磁共振横向弛豫时间T2分布,与现有技术相比提高了储层评价准确性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种获得目的储层在完全含水状态下核磁共振的横向弛豫时间分布的方法的流程图;
图2是目的储层温度压力下实验所测的岩心中地层水的核磁共振横向弛豫时间分布与将原油的本征弛豫和扩散弛豫综合贡献替换为地层水的本征弛豫和扩散弛豫综合贡献,获取的目的储层条件下完全含水状态的横向弛豫时间分布;
图3是获得目的储层在完全含水状态下核磁共振的横向弛豫时间分布的装置示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
现有技术中为了获取完全含水状态下的核磁共振横向弛豫时间T2分布,进行准确储层评价,可以利用自由水核磁共振横向弛豫时间T2分布,但是忽略了实际储层中油气会对核磁共振横向弛豫时间T2分布产生的影响,无法准确的进行目的储层评价。
为解决上述问题,下面详细说明本申请实施例的具体实现。
如图1所示的是本申请实施例提供的一种获得目的储层在完全含水状态下核磁共振的横向弛豫时间分布的方法的流程图,该方法包括:
S201:根据目的储层的温度和压力信息,设置相应的实验温度和压力,获取目的储层岩心中目的储层地层水的核磁共振本征横向弛豫时间T2分布,目的储层岩心中目的储层地层水是自由状态的地层水。
在实际中,将目的储层的岩心洗油洗盐处理后,配比与实际储层矿化度相同的地层水,并进行抽空加压饱和,也即地层水在岩心的孔隙中,将饱和好的岩心装入核磁夹持器中。根据实际地层信息,设计相应的实验温度和压力(80℃,20Mpa),压力包括围压(为20Mpa)和孔压(为20Mpa);利用频率为2MHz的核磁共振岩心分析仪,并采用CPMG脉冲序列确定目的储层条件下完全含水状态下的自旋回波信号,再将自旋回波信号进行反演得到岩心中地层水的核磁共振横向弛豫时间T2w,也即自由状态下地层水的核磁共振横向弛豫时间T2w;由于没有采用脉冲梯度,此时的横向弛豫时间没有扩散弛豫贡献,属于本征弛豫时间。
CPMG脉冲序列为一个90°脉冲后跟着一系列180°脉冲的脉冲序列,是由Carr,Pucell,Meiboom和Gill的名字命名的。
S202:根据目的储层温度和压力信息,设置相应的实验温度和压力,获取目的储层原油的核磁共振本征横向弛豫时间T2分布和目的储层地层水的核磁共振本征横向弛豫时间T2分布,根据目的储层原油的核磁共振本征横向弛豫时间T2分布确定原油的含氢指数,根据目的储层的地层水的核磁共振本征横向弛豫时间T2分布确定地层水的含氢指数。
在实际中,将目的储层原油样品和地层水样品分别装入核磁夹持器中,根据目的储层温度压力信息,设计相应的实验温度和压力(80℃,20Mpa),利用频率为2MHz的核磁共振岩心分析仪,没有施加脉冲梯度磁场,并采用CPMG脉冲序列分别确定原油的自旋回波信号和地层水的自旋回波信号,再将原油的自旋回波信号进行反演得到原油的核磁共振本征横向弛豫时间T2分布,地层水的自旋回波信号进行反演的到地层水的核磁共振本征横向弛豫时间T2分布,再根据原油的核磁共振本征横向弛豫时间T2分布确定原油的含氢指数IHo,根据地层水的核磁共振本征横向弛豫时间T2分布确定地层水的含氢指数IHw,含氢指数具体是通过下式获取:
I H = Σ i n P i / V Σ i n P i , s / V s
上式中,为求和公式;IH为所测流体的含氢指数;Pi为所测流体的核磁共振横向弛豫时间T2分布的第i个区间孔隙度分量;V为所测流体的体积;Pi,s为标准水样的核磁共振横向弛豫时间T2的第i个区间孔隙度分量,标准水样是指普通的纯净水;Vs为标准水样的体积。
S203:获取与目的储层的温度和压力相同的温度和压力下原油的扩散系数、以及获取与目的储层的温度和压力相同的温度和压力下地层水的扩散系数。
在实际中,将目的储层原油样品和地层水样品分别装入核磁夹持器中,根据目的储层温度压力信息,设计相应的实验温度和压力(80℃,20Mpa),在脉冲梯度磁场的条件下,采用PFGSE脉冲序列分别确定原油的自旋回波信号和地层水的自旋回波信号;根据原油的自旋回波信息确定原油回波间隔和原油回波幅度,上述原油回波间隔(X轴)和原油回波幅度(Y轴)之间的斜率是原油的扩散系数;根据地层水的自旋回波信息确定地层水回波间隔和地层水回波幅度,上述地层水回波间隔(X轴)和地层水回波幅度(Y轴)之间的斜率是地层水的扩散系数。
S204:利用核磁共振测井获取目的储层下包含原油、地层水和岩心的自旋回波信号,将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号。
在实际中,原油的本征弛豫和扩散弛豫综合贡献是通过原油的本征横向弛豫时间T2分布、原油的含氢指数和原油的扩散系数确定的,自由状态地层水的本征弛豫和扩散弛豫综合贡献是通过自由状态的地层水的本征横向弛豫时间T2分布、地层水的含氢指数和地层水的扩散系数确定的。将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号具体通过公式获取:
ECHOW ( t ) = ECHO ( t ) - φ S o I Ho e - t [ 1 T 2 o + D o ( γG T E ) 2 12 ] + φ S o I Hw e - t [ 1 T 2 W + D w ( γG T E ) 2 12 ]
其中,ECHOW(t)为目的储层完全含水状态的自旋回波信号;ECHO(t)为目的储层包含原油、地层水和岩心自旋回波信号,是核磁共振测井获取的;是原油的本征弛豫和扩散弛豫综合贡献;是自由状态地层水的本征弛豫和扩散弛豫综合贡献;φ为核磁共振孔隙度,可以从回波串中求取,具体步骤为将自旋回波串进行反演,求取T2的第i个孔隙度分量Pi,将孔隙度分量Pi求和得到核磁共振孔隙度;So为含油饱和度,可以利用核磁或者电阻率等方法求取,如利用核磁中的差谱法、移谱法或者二维核磁共振技术求取,或者利用电阻率的archie公式求取;γ为氢核的旋磁比,是常数,单位为rad·/(s·T);G为磁场梯度,为常数,单位为G/cm,可通过仪器确定;TE为CPMG脉冲序列的回波间隔,为常数,可以从仪器采集模式获取;IHo为原油的含氢指数;IHw为地层水的含氢指数;T2o为原油的本征弛豫时间;T2w为目的储层岩心中自由状态的地层水的本征横向弛豫时间;Do为目的储层原油的扩散系数;Dw为目的储层地层水的扩散系数;t为采集时间。
S205:将目的储层在完全含水状态下的自旋回波信号进行反演,获取目的储层在完全含水状态下的核磁共振横向弛豫时间T2分布。
在实际中,将自旋回波信号进行反演获取横向弛豫时间T2分布可以通过最小二乘法计算获得,但本领域技术人员应该清楚在实际中还可以是其他反演算法,例如奇异值分解法(即SVD法)和模平滑法,本申请实施例并不以此为限。
如图2所示的分别是目的储层温度压力(80℃,20Mpa)下实验所测的岩心中地层水的核磁共振横向弛豫时间T2分布201,与将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态的地层水的本征弛豫和扩散弛豫综合贡献,进行校正获取的目的储层条件下完全含水状态的横向弛豫时间T2分布202,其中横坐标为核磁共振横向弛豫时间T2,单位为ms,纵坐标为核磁共振横向弛豫时间T2分布的区间孔隙度分量Pi。可见采用校正得到的目的储层在完全含水状态下的核磁共振横向弛豫时间T2分布与目的储层实验所测的岩心中地层水的核磁共振横向弛豫时间T2分布相符合。
由此可见,本申请实施例提供的一种获得目的储层在完全含水状态下核磁共振的横向弛豫时间分布的方法,根据目的储层的温度和压力信息,设置好实验温度和压力,获取目的储层自由状态的地层水的核磁共振横向弛豫时间T2w分布、目的储层原油核磁共振横向弛豫时间T2o分布、目的储层原油的含氢指数IHo、目的储层地层水的含氢指数IHw、目的储层原油的扩散系数Do和目的储层地层水的扩散系数Dw,根据这些数据将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号,再将自旋回波信号进行反演获得目的储层在完全含水状态下的核磁共振横向弛豫时间T2分布,与现有技术相比提高了储层评价准确性。
本申请实施例还相应提供一种获得目的储层在完全含水状态下核磁共振的横向弛豫时间分布的装置300,如图3所示,该装置300包括:第一横向弛豫时间分布获取单元301、第二横向弛豫时间分布获取单元302、含氢指数获取单元303、扩散系数获取单元304、自旋回波信号获取单元305和第三横向弛豫时间分布获取单元306,其中:
第一横向弛豫时间分布获取单元301用于获取与目的储层的温度和压力相同的温度和压力下的目的储层岩心中目的储层地层水的核磁共振本征横向弛豫时间T2分布,上述目的储层岩心中目的储层地层水是自由状态的地层水;
第二横向弛豫时间分布获取单元302用于获取与目的储层的温度和压力相同的温度和压力下的目的储层原油的核磁共振本征横向弛豫时间T2分布,以及获取与目的储层的温度和压力相同的温度和压力下的目的储层地层水的核磁共振本征横向弛豫时间T2分布;
含氢指数获取单元303用于根据第二横向弛豫时间分布获取单元302得到的目的储层原油的核磁共振本征横向弛豫时间T2分布确定原油的含氢指数,用于根据第二横向弛豫时间分布单元302得到的目的储层的地层水的核磁共振本征横向弛豫时间T2分布确定地层水的含氢指数;
扩散系数获取单元304用于获取与目的储层的温度和压力相同的温度和压力下所述原油的扩散系数,以及获取与目的储层的温度和压力相同的温度和压力下地层水的扩散系数;
自旋回波信号获取单元305用于利用核磁共振测井获取目的储层下包含原油、地层水和岩心的自旋回波信号,将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号;原油的本征弛豫和扩散弛豫综合贡献是通过原油的本征横向弛豫时间T2分布、原油的含氢指数和原油的扩散系数确定的,自由状态地层水的本征弛豫和扩散弛豫综合贡献是通过所述自由状态地层水的本征横向弛豫时间T2分布、地层水的含氢指数和地层水的扩散系数确定的;
第三横向弛豫时间分布获取单元306用于将自旋回波信号获取单元305得到的目的储层在完全含水状态下的自旋回波信号进行反演,获取目的储层在完全含水状态下的核磁共振横向弛豫时间T2分布。
由此可见,本申请实施例提供的一种获得目的储层在完全含水状态下核磁共振的横向弛豫时间分布的装置,根据目的储层的温度和压力信息,设置好实验温度和压力,获取目的储层自由状态的地层水的核磁共振横向弛豫时间T2w分布、目的储层原油核磁共振横向弛豫时间T2o分布、目的储层原油的含氢指数IHo、目的储层地层水的含氢指数IHw、目的储层原油的扩散系数Do和目的储层地层水的扩散系数Dw,根据这些数据将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号,再将自旋回波信号进行反演获得目的储层条件下完全含水状态的核磁共振横向弛豫时间T2分布,与现有技术相比提高了储层评价准确性。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (10)

1.一种获得目的储层横向弛豫时间分布的方法,其特征在于,该方法包括:
S1:根据目的储层的温度和压力信息,设置相应的实验温度和压力,获取目的储层岩心中目的储层地层水的核磁共振本征横向弛豫时间T2分布,所述目的储层岩心中目的储层地层水是自由状态的地层水;
S2:根据目的储层温度和压力信息,设置相应的实验温度和压力,获取目的储层原油的核磁共振本征横向弛豫时间T2分布和目的储层地层水的核磁共振本征横向弛豫时间T2分布,根据所述目的储层原油的核磁共振本征横向弛豫时间T2分布确定原油的含氢指数,根据所述目的储层的地层水的核磁共振本征横向弛豫时间T2分布确定地层水的含氢指数;
S3:获取与所述目的储层的温度和压力相同的温度和压力下所述原油的扩散系数,以及获取与所述目的储层的温度和压力相同的温度和压力下所述地层水的扩散系数;
S4:利用核磁共振测井获取目的储层下包含原油、地层水和岩心的自旋回波信号,将所述原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号;所述原油的本征弛豫和扩散弛豫综合贡献是通过所述原油的本征横向弛豫时间T2分布、所述原油的含氢指数和所述原油的扩散系数确定的,所述自由状态地层水的本征弛豫和扩散弛豫综合贡献是通过所述自由状态的地层水的本征横向弛豫时间T2分布、所述地层水的含氢指数和所述地层水的扩散系数确定的;
S5:将所述目的储层在完全含水状态下的自旋回波信号进行反演,获取目的储层在完全含水状态下的核磁共振横向弛豫时间T2分布。
2.根据权利要求1所述的方法,其特征在于,所述自由状态的地层水的核磁共振本征横向弛豫时间T2分布是采用CPMG脉冲序列得到自旋回波信号,将自旋回波信号反演得到的。
3.根据权利要求1所述的方法,其特征在于,所述目的储层原油的核磁共振本征横向弛豫时间T2分布是采用CPMG脉冲的脉冲序列得到自旋回波信号,将自旋回波信号反演得到的。
4.根据权利要求1所述的方法,其特征在于,所述目的储层地层水的核磁共振本征横向弛豫时间T2分布是采用CPMG的脉冲序列得到自旋回波信号,将自旋回波信号反演得到的。
5.根据权利要求1所述的方法,其特征在于,所述目的储层原油的含氢指数是利用标准水样将所述原油的核磁共振本征横向弛豫时间T2分布进行刻度,计算的公式如下:
I Ho = Σ i n P i / V 0 Σ i n P i , s / V s
上式中,为求和公式;IHo为所述原油的含氢指数;Pi为所述原油的核磁共振横向弛豫时间T2分布的第i个区间孔隙度分量;V0为所述原油体积;Pi,s为标准水样的核磁共振横向弛豫时间T2的第i个区间孔隙度分量;Vs为标准水样的体积。
6.根据权利要求1所述的方法,其特征在于,所述目的储层地层水的含氢指数是利用标准水样将所述地层水的核磁共振本征横向弛豫时间T2分布进行刻度,计算的公式如下:
I Hw = Σ i n P i ′ / V w Σ i n P i , s / V s
上式中,为求和公式;IHw为所述地层水的含氢指数;Pi'为所述地层水的核磁共振横向弛豫时间T2分布的第i个区间孔隙度分量;Vw为所述地层水体积;Pi,s为标准水样的核磁共振横向弛豫时间T2的第i个区间孔隙度分量;Vs为标准水样的体积。
7.根据权利要求1所述的方法,其特征在于,所述目的储层原油的扩散系数是在脉冲梯度磁场的条件下,采用脉冲梯度场自旋回波(PFGSE)脉冲序列得到原油自旋回波信号,获取原油回波间隔和原油回波幅度,所述原油回波间隔和所述原油回波幅度之间的斜率是所述原油的扩散系数。
8.根据权利要求1所述的方法,其特征在于,所述目的储层地层水的扩散系数是在脉冲梯度磁场的条件下,采用PFGSE脉冲序列得到地层水自旋回波信号,获取地层水回波间隔和地层水回波幅度,所述地层水回波间隔和所述地层水回波幅度之间的斜率是所述地层水的扩散系数。
9.根据权利要求1所述的方法,其特征在于,所述将原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号具体通过公式获取:
ECHOW ( t ) = ECHO ( t ) - φ S 0 I Ho e - t [ 1 T 2 o + D o ( γG T E ) 2 12 ] + φ S o I Hw e - t [ 1 T 2 W + D w ( γG T E ) 2 12 ]
其中,ECHOW(t)为目的储层在完全含水状态下的自旋回波信号;ECHO(t)为目的储层包含原油、地层水和岩心的自旋回波信号;φ为核磁共振孔隙度;So为含油饱和度;γ为氢核的旋磁比,是常数;G为磁场梯度;TE为CPMG脉冲序列的回波间隔;IHo为原油的含氢指数;IHw为地层水的含氢指数;T2o为原油的本征弛豫时间;T2w为目的储层岩心中自由状态的地层水的本征横向弛豫时间;Do为目的储层原油的扩散系数;Dw为目的储层地层水的扩散系数;t为采集时间。
10.一种获得目的储层横向弛豫时间分布的装置,其特征在于,该装置包括:
第一横向弛豫时间分布获取单元,用于获取与所述目的储层的温度和压力相同的温度和压力下的目的储层岩心中目的储层地层水的核磁共振本征横向弛豫时间T2分布,所述目的储层岩心中目的储层地层水是自由状态的地层水;
第二横向弛豫时间分布获取单元,用于获取与所述目的储层的温度和压力相同的温度和压力下的目的储层原油的核磁共振本征横向弛豫时间T2分布,以及获取与所述目的储层的温度和压力相同的温度和压力下的目的储层地层水的核磁共振本征横向弛豫时间T2分布;
含氢指数获取单元,用于根据第二横向弛豫时间分布获取单元得到的所述目的储层原油的核磁共振本征横向弛豫时间T2分布确定原油的含氢指数,用于根据第二横向弛豫时间分布得到的所述目的储层的地层水的核磁共振本征横向弛豫时间T2分布确定地层水的含氢指数;
扩散系数获取单元,用于获取与所述目的储层的温度和压力相同的温度和压力下所述原油的扩散系数,以及获取与所述目的储层的温度和压力相同的温度和压力下所述地层水的扩散系数;
自旋回波信号获取单元,用于利用核磁共振测井获取目的储层下包含原油、地层水和岩心的自旋回波信号,将所述原油的本征弛豫和扩散弛豫综合贡献替换为自由状态地层水的本征弛豫和扩散弛豫综合贡献,获取目的储层在完全含水状态下的自旋回波信号;所述原油的本征弛豫和扩散弛豫综合贡献是通过所述原油的本征横向弛豫时间T2分布、所述原油的含氢指数和所述原油的扩散系数确定的,所述自由状态地层水的本征弛豫和扩散弛豫综合贡献是通过所述自由状态的地层水的本征横向弛豫时间T2分布、所述地层水的含氢指数和所述地层水的扩散系数确定的;
第三横向弛豫时间分布获取单元,用于将自旋回波信号获取单元得到的目的储层在完全含水状态下的自旋回波信号进行反演,获取目的储层在完全含水状态下的核磁共振横向弛豫时间T2分布。
CN201410586502.7A 2014-10-28 2014-10-28 一种获得目的储层横向弛豫时间分布的方法及装置 Active CN104330433B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410586502.7A CN104330433B (zh) 2014-10-28 2014-10-28 一种获得目的储层横向弛豫时间分布的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410586502.7A CN104330433B (zh) 2014-10-28 2014-10-28 一种获得目的储层横向弛豫时间分布的方法及装置

Publications (2)

Publication Number Publication Date
CN104330433A true CN104330433A (zh) 2015-02-04
CN104330433B CN104330433B (zh) 2016-08-03

Family

ID=52405203

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410586502.7A Active CN104330433B (zh) 2014-10-28 2014-10-28 一种获得目的储层横向弛豫时间分布的方法及装置

Country Status (1)

Country Link
CN (1) CN104330433B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105464654A (zh) * 2015-12-15 2016-04-06 中国石油天然气股份有限公司 确定致密砂岩饱和度指数的方法及装置
CN106066494A (zh) * 2016-05-24 2016-11-02 中国地质大学(北京) 一种火成岩核磁共振孔隙度校正方法和t2分布校正方法
CN107941837A (zh) * 2017-11-14 2018-04-20 山东省分析测试中心 低场核磁共振鉴别阿胶品质的方法
CN108005633A (zh) * 2017-11-27 2018-05-08 中石化石油工程技术服务有限公司 致密气藏的二维核磁共振测井观测模式参数确定方法
CN108458960A (zh) * 2018-03-27 2018-08-28 中国石油大学(华东) 一种富有机质泥页岩含氢组分、孔隙度及孔径的评价方法
CN110244369A (zh) * 2019-06-28 2019-09-17 中国石油大学(北京) 储层束缚和可动流体分布确定方法、装置及系统
CN112710688A (zh) * 2019-10-24 2021-04-27 中国石油天然气股份有限公司 核磁共振纵向弛豫获取方法及系统
US11143607B2 (en) * 2020-03-13 2021-10-12 King Fahd University Of Petroleum And Minerals Method for evaluation of permeability anisotropy using NMR diffusion measurements for oil and gas wells
CN113781452A (zh) * 2021-09-15 2021-12-10 大庆油田有限责任公司 页岩含油饱和度的确定方法及装置、电子设备和存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08313462A (ja) * 1995-05-19 1996-11-29 Japan Magnet Technol Kk 西瓜の熟成度評価方法
CN1231428A (zh) * 1998-04-08 1999-10-13 施卢默格海外有限公司 利用核磁共振和其它测井数据评估岩层构造的方法
WO2001042817A1 (en) * 1999-12-10 2001-06-14 Schlumberger Limited Nuclear magnetic resonance method and logging apparatus
CN1132019C (zh) * 1995-06-26 2003-12-24 纽马公司 采用扩散和弛豫测井测量技术评估地岩的核磁共振系统和方法
US20040027123A1 (en) * 2002-08-09 2004-02-12 Heaton Nicholas J. Method for Detecting Hydrocarbons from NMR Data
CN102368095A (zh) * 2011-09-10 2012-03-07 吉林大学 多指数拟合核磁共振地下水探测信号弛豫时间谱提取方法
CN102944571A (zh) * 2012-10-17 2013-02-27 中国地质大学(北京) 一种测量煤中不同状态水分含量的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08313462A (ja) * 1995-05-19 1996-11-29 Japan Magnet Technol Kk 西瓜の熟成度評価方法
CN1132019C (zh) * 1995-06-26 2003-12-24 纽马公司 采用扩散和弛豫测井测量技术评估地岩的核磁共振系统和方法
CN1231428A (zh) * 1998-04-08 1999-10-13 施卢默格海外有限公司 利用核磁共振和其它测井数据评估岩层构造的方法
WO2001042817A1 (en) * 1999-12-10 2001-06-14 Schlumberger Limited Nuclear magnetic resonance method and logging apparatus
US20040027123A1 (en) * 2002-08-09 2004-02-12 Heaton Nicholas J. Method for Detecting Hydrocarbons from NMR Data
CN102368095A (zh) * 2011-09-10 2012-03-07 吉林大学 多指数拟合核磁共振地下水探测信号弛豫时间谱提取方法
CN102944571A (zh) * 2012-10-17 2013-02-27 中国地质大学(北京) 一种测量煤中不同状态水分含量的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
石玉江 刘天定: "二维核磁共振测井技术在长庆油田的应用", 《测井技术》 *
谭茂金 邹友龙: "(T2,D)二维核磁共振测井混合反演方法与参数影响分析", 《地球物理学报》 *

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105464654B (zh) * 2015-12-15 2018-06-01 中国石油天然气股份有限公司 确定致密砂岩饱和度指数的方法及装置
CN105464654A (zh) * 2015-12-15 2016-04-06 中国石油天然气股份有限公司 确定致密砂岩饱和度指数的方法及装置
CN106066494A (zh) * 2016-05-24 2016-11-02 中国地质大学(北京) 一种火成岩核磁共振孔隙度校正方法和t2分布校正方法
CN106066494B (zh) * 2016-05-24 2018-05-18 中国地质大学(北京) 一种火成岩核磁共振孔隙度校正方法和t2分布校正方法
CN107941837A (zh) * 2017-11-14 2018-04-20 山东省分析测试中心 低场核磁共振鉴别阿胶品质的方法
CN108005633A (zh) * 2017-11-27 2018-05-08 中石化石油工程技术服务有限公司 致密气藏的二维核磁共振测井观测模式参数确定方法
CN108458960A (zh) * 2018-03-27 2018-08-28 中国石油大学(华东) 一种富有机质泥页岩含氢组分、孔隙度及孔径的评价方法
CN108458960B (zh) * 2018-03-27 2019-10-29 中国石油大学(华东) 富有机质泥页岩含氢组分、孔隙度及孔径的评价方法
CN110244369A (zh) * 2019-06-28 2019-09-17 中国石油大学(北京) 储层束缚和可动流体分布确定方法、装置及系统
CN112710688B (zh) * 2019-10-24 2023-08-22 中国石油天然气股份有限公司 核磁共振纵向弛豫获取方法及系统
CN112710688A (zh) * 2019-10-24 2021-04-27 中国石油天然气股份有限公司 核磁共振纵向弛豫获取方法及系统
US11143607B2 (en) * 2020-03-13 2021-10-12 King Fahd University Of Petroleum And Minerals Method for evaluation of permeability anisotropy using NMR diffusion measurements for oil and gas wells
US20210372956A1 (en) * 2020-03-13 2021-12-02 King Fahd University Of Petroleum And Minerals Downhole method for determining geologic permeability
US20210372955A1 (en) * 2020-03-13 2021-12-02 King Fahd University Of Petroleum And Minerals Method for determining permeability in a hydrocarbon formation
US20220003695A1 (en) * 2020-03-13 2022-01-06 King Fahd University Of Petroleum And Minerals Method for downhole determination of permeability anisotropy using nmr
US20220003696A1 (en) * 2020-03-13 2022-01-06 King Fahd University Of Petroleum And Minerals NMR method for determining permeability in geologic formation
US11493462B2 (en) * 2020-03-13 2022-11-08 King Fahd University Of Petroleum And Minerals Measurement method for determining resistivity and permeability in a borehole
US11525792B2 (en) * 2020-03-13 2022-12-13 King Fahd University Of Petroleum And Minerals NMR method for determining permeability in geologic formation
US11543369B2 (en) * 2020-03-13 2023-01-03 King Fahd University Of Petroleum And Minerals Method for determining permeability in a hydrocarbon formation
US11543370B2 (en) * 2020-03-13 2023-01-03 King Fahd University Of Petroleum And Minerals Downhole method for determining geologic permeability
US11555795B2 (en) * 2020-03-13 2023-01-17 King Fahd University Of Petroleum And Minerals Method for downhole determination of permeability anisotropy using NMR
US20210372954A1 (en) * 2020-03-13 2021-12-02 King Fahd University Of Petroleum And Minerals Measurement method for determining resistivity and permeability in a borehole
CN113781452A (zh) * 2021-09-15 2021-12-10 大庆油田有限责任公司 页岩含油饱和度的确定方法及装置、电子设备和存储介质

Also Published As

Publication number Publication date
CN104330433B (zh) 2016-08-03

Similar Documents

Publication Publication Date Title
CN104330433B (zh) 一种获得目的储层横向弛豫时间分布的方法及装置
US10444171B2 (en) Absolute porosity and pore size determination of pore types in media with varying pore sizes
CN104215652B (zh) 确定油气饱和度的方法和装置
Yan et al. A robust NMR method to measure porosity of low porosity rocks
CN104453874A (zh) 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法
WO2014126883A3 (en) Estimating molecular size distributions in formation fluid samples using a downhole nmr fluid analyzer
CN103675722B (zh) 岩石t2-g实验采集参数自动匹配方法
CN103674811A (zh) 一种核磁共振孔隙度测量的校正方法、装置及系统
US10534055B2 (en) NMR method for determining non-oil volume of a rock sample
Liang et al. Two-dimensional nuclear magnetic resonance method for wettability determination of tight sand
CA3107041A1 (en) Assessment of inaccessible pore volume for polymer flooding
AU2017204581A1 (en) Nuclear magnetic resonance tool calibration
CN106290103B (zh) 一种页岩气储层中粘土微孔孔隙度的测定方法
CN106644879A (zh) 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置
Bernard Instruments and field work to measure a magnetic resonance sounding
CN110029990B (zh) 一种核磁共振测井方法和装置
CN106066494B (zh) 一种火成岩核磁共振孔隙度校正方法和t2分布校正方法
Hiller et al. First Measurements of Surface Nuclear Magnetic Resonance Signals Without an Oscillating Excitation Pulse–Exploiting Non‐Adiabatic Prepolarization Switch‐Off
Mitchell et al. A general approach to T2 measurements in the presence of internal gradients
CN106556876B (zh) 一种基于多频偏共振激发的三维核磁共振叠前反演方法
Liu et al. Characterizing fluid presence and transport in rock cores at reservoir-like conditions via spatially resolved NMR relaxation/diffusion maps
XIE et al. A method for multiple echo trains jointing inversion of NMR relaxation measurements
Pan et al. Advantages of the optimum pulse moment in surface nmr and application in groundwater exploration
US20150293195A1 (en) Nuclear magnetic resonance t2 recovery pulse
CN114577836A (zh) 孔隙信息确定方法、装置、设备、介质及程序产品

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