CN104777181B - 致密油核磁共振t2截止值及流体饱和度确定方法、装置 - Google Patents
致密油核磁共振t2截止值及流体饱和度确定方法、装置 Download PDFInfo
- Publication number
- CN104777181B CN104777181B CN201410855862.2A CN201410855862A CN104777181B CN 104777181 B CN104777181 B CN 104777181B CN 201410855862 A CN201410855862 A CN 201410855862A CN 104777181 B CN104777181 B CN 104777181B
- Authority
- CN
- China
- Prior art keywords
- carbon content
- pore throat
- throat radius
- lower limit
- cutoff
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000012530 fluid Substances 0.000 title claims abstract description 32
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 21
- 239000011148 porous material Substances 0.000 claims abstract description 121
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 103
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 103
- 239000011435 rock Substances 0.000 claims abstract description 28
- 238000012360 testing method Methods 0.000 claims description 68
- 238000005259 measurement Methods 0.000 abstract 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 42
- 238000002474 experimental method Methods 0.000 description 26
- 238000009826 distribution Methods 0.000 description 12
- 238000000280 densification Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000005311 nuclear magnetism Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000005119 centrifugation Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000003556 assay Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000003340 mental effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000005349 anion exchange Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种致密油核磁共振T2截止值及流体饱和度确定方法及装置,方法包括:对待测样品进行碳含量测定,生成碳含量测定结果;根据碳含量测定结果确定待测样品的孔喉半径下限值;根据孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值。本方案利用实际样品测定的致密油可动流体孔喉半径下限值求取T2截止值,由于实际样品测定的致密油可动流体孔喉半径下限值采用由粗变细孔喉测定碳含量逐渐逼近的方法,更加具有真实性。
Description
技术领域
本发明涉及石油勘探技术领域,具体的讲是致密油核磁共振T2截止值及流体饱和度确定方法及装置。
背景技术
核磁共振实验技术是目前监测岩石可动流体饱和度的一项重要技术。核磁共振实验中T2弛豫时间谱包含岩石孔隙结构分布和所含流体的信息,这种信息是核磁共振实验测定岩石可动流体饱和度的基础所在。T2分布与孔隙尺寸相关。较大孔隙中的可动流体与岩石孔隙固体表面作用力弱,T2弛豫时间较大,相比而言,束缚流体与岩石孔隙固体表面作用力强,T2弛豫时间较小。
T2截止值的研究是核磁实验进行可动流体与束缚流体区分的关键,大于T2截止值信号流体为可动流体,相应的为束缚流体。目前确定T2截止值主要有储层系数与T2截止值关系计算、毛管力与T2截止值的关系计算、T2几何平均值求取、阴离子交换容量计算、T2谱峰形态、离心实验等方法。
离心法是现有技术中计算T2截止值最为普遍与准确的方法。尽管如此,离心法测定T2截止值并不是适合于所有的致密油。图1为离心实验准确测定致密油中束缚水饱和度适用性示意图。其中,图1(A):大孔型样品含水饱和度与离心力关系;图1(B):大孔与小孔并存型样品离心样品含水饱和度与离心力关系,图1(C):大孔型样品孔隙分布,图1(D):大孔与小孔并存型样品样品孔隙分布。Sw:含水饱和度;Swi:样品实际束缚水饱和度;Swi’:离心实验测定束缚水饱和度;P:离心实验最大离心力;R:孔喉半径;P1:大孔型样品可动水离心出所需要的最大离心力;P2:大孔与小孔并存型样品小孔中可动水离心出所需要的最小离心力;P3:大孔与小孔并存型样品小孔中可动水离心出所需要的最大离心力;)。
如图1(C)所示为大孔型样品孔隙分布,对于孔隙分布为单峰型的大孔型样品,如图1(A)所示,如果孔隙中可动水离心出所需要的最大离心力P1小于离心实验最大离心力P,在实验过程中,随着离心力的逐渐增加,含水饱和度快速降低,当离心力达到P1并进一步增加至P时,含水饱和度则基本没有较大变化,此时即认为样品的束缚水饱和度为Swi。实验测定的可动水与束缚水含量是准确的。但是对于如图1(D)所示的孔隙分布为双峰型的大孔与小孔并存型样品,如图1(B)所示,由于小孔隙中可动水离心出所需要的最小离心力P1大于离心实验最大离心力P,则实验测定过程中即使达到了最大离心力P,也无法将小孔中的可动水离出。但是实际离心实验过程中,由于含水饱和度的变化在离心力从零至P1逐渐增加的过程汇总快速降低,而P1至P过程中,含水饱和度同样基本没有较大变化,此时实验结论则错误性地将实验测得的束缚水饱和度Swi’认为是实际的束缚水饱和度,因此此时计算得到可动水与束缚水含量是不准确的。只有当离心力增加足够大达到小孔中可动水离心出所需要的最小离心力P2时,随着压力的增加,小孔中的可动水才能够离心出,此时的束缚水饱和度为图1B中的Swi。
由此可见,实际的离心实验在处理致密油特别是孔隙非均质性强的样品时,存在诸多不准确性。离心实验中的离心力需要足够大才能够将小孔中的可动水离心出,但是目前离心实验技术中的离心力往往不能够达到相应要求,最重要的是,很难将孔隙分布与离心力进行定量对应。
发明内容
为更加准确、真实的确定致密油核磁共振信号T2截止值,本发明实施例提供了一种致密油核磁共振T2截止值确定方法,方法包括:
对待测样品进行碳含量测定,生成碳含量测定结果;
根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值;
所述的根据所述碳含量测定结果确定待测样品的孔喉半径下限值包括:
根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
此外,本发明实施例还提供一种致密油核磁共振T2截止值确定装置,包括:
测定模块,用于对待测样品进行碳含量测定,生成碳含量测定结果;
孔喉半径下限值确定模块,用于根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
T2截止值确定模块,用于根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值;
所述的孔喉半径下限值确定模块包括:
孔喉半径确定单元,用于根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
下限值确定单元,用于根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
同时,利用上述获得的T2截止值,利用大于T2截止值信号幅值面积与总信号幅值面积比值即可计算可动流体饱和度,本发明还提供一种可动流体饱和度确定方法,方法包括:
确定待测样品的T2截止值;
根据确定的T2截止值的信号幅值面积与总信号面积比值确定可动流体饱和度,其中确定待测样品的T2截止值包括:
对待测样品进行碳含量测定,生成碳含量测定结果;
根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值;
所述的根据所述碳含量测定结果确定待测样品的孔喉半径下限值包括:
根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
同时,本发明还提供一种可动流体饱和度确定装置,装置包括:
截止值确定模块,用于确定待测样品的T2截止值;
饱和度确定模块,用于根据确定的T2截止值的信号幅值面积与总信号面积比值确定可动流体饱和度,其中截止值确定模块包括:
测定模块,用于对待测样品进行碳含量测定,生成碳含量测定结果;
孔喉半径下限值确定模块,根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
T2截止值确定模块,用于根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值;
所述的孔喉半径下限值确定模块包括:
孔喉半径确定单元,用于根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
下限值确定单元,用于根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
本方案利用实际样品测定的致密油可动流体孔喉半径下限值求取T2截止值,由于实际样品测定的致密油可动流体孔喉半径下限值采用由粗变细孔喉测定碳含量逐渐逼近的方法,其数值真实反映了油在吼道中实际能够充注的最小半径值,即可动孔喉半径下限值,因此本方法相比离心法等其他方法更加具有真实性。此外,现有技术进行离心实验中的离心力需要足够大才能够将小孔中的可动水离心出,但是目前离心实验技术中的离心力往往不能够达到相应要求,最重要的是,在实际离心实验中很难判定是否已经完全将可动水离心出岩心,只能通过判断可动水随着离心力继续增大不再明显减少就停止离心,而这种判断明显不具备科学性。本方案通过准确测定致密油可动流体孔喉半径下限值,带入公式即可准确求取T2截止值,进而求取准确的可动流动饱和度。
为让本发明的上述和其他目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附图式,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为现有技术中进行离心实验准确测定致密油中束缚水饱和度适用性示意图;
图2为本发明提供的一种致密油核磁共振T2截止值确定方法的流程图;
图3为本发明实施例中提供的确定待测样品的孔喉半径下限值的流程图;
图4为本发明提供的一种致密油核磁共振T2截止值确定装置的框图;
图5为本发明实施例中确定致密砂岩储层聚集孔喉半径下限测定方法原理图;
图6为本发明方法计算得到的T2截止值与现有技术计算得到的T2截止值的结果比较图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图2所示,本发明提供一种致密油核磁共振T2截止值确定方法,包括:
步骤S301,对待测样品进行碳含量测定,生成碳含量测定结果;
步骤S302,根据碳含量测定结果确定待测样品的孔喉半径下限值;
步骤S303,根据孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值。
具体的,如图3所示,本发明实施例中根据碳含量测定结果确定待测样品的孔喉半径下限值包括:
步骤S3021,根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
步骤S3022,根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
具体的,本方案具体实施时,根据所述孔喉半径下限值和待测样品的岩心弛豫率利用公式r=2ρ2T2确定对应的T2截止值,其中,
ρ2:核磁共振信号T2的表面弛豫率,μm/s;
r为所述孔隙半径下限值,nm。
同时,本发明还提供一种致密油核磁共振信号T2截止值确定装置,如图4所示,包括:
测定模块301,用于对待测样品进行碳含量测定,生成碳含量测定结果;
孔喉半径下限值确定模块302,用于根据碳含量测定结果确定待测样品的孔喉半径下限值;
T2截止值确定模块303,用于根据孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值。
具体的,本发明的具体实施时,孔喉半径下限值确定模块包括:
孔喉半径确定单元,用于根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
下限值确定单元,用于根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
下面结合本方案的原理及具体实施例对本方案做进一步详细说明:
核磁实验中,孔隙中的流体有三种不同的弛豫机制:自由弛豫、表面弛豫和扩散弛豫。弛豫时间T2可表示为:
当采用短TE且孔隙只含饱和流体时,表面弛豫起主要作用,即T2直接与孔隙尺寸成正比:
由公式(2)可知,T2分布图实际上反映了孔隙尺寸的分布。假设孔隙是一个半径为r的圆柱,则公式(2)可简化为:
则:r=2ρ2T2 (4)
符号注释:
T2:通过脉冲序列(CPMG)采集的孔隙流体的横向弛豫时间,ms;
T2z:在足够大的容器中(大到容器影响可忽略不计)孔隙流体的横向弛豫时间,ms;
T2b:表面弛豫引起的横向弛豫时间,ms;
T2k:磁场梯度下由扩散引起的孔隙流体的横向弛豫时间,ms;
ρ2:T2表面弛豫率,μm/s;
孔隙的比表面积,1/μm。
r:孔隙半径,nm。
测定出致密油可动流体孔喉半径下限值,通过公式(4)计算即可得到对应的核磁T2弛豫时间截止值,之后通过信号幅值面积比值即可计算相应可动流体饱和度。
本发明公开的核磁T2弛豫时间截止值的方案具体步骤如下:
步骤(1):选取一块致密油藏产油层段砂岩样品,通过测定样品中由粗变细厚度碳含量确定致密油可动流体孔喉半径下限值rcutoff。如图4所示,其中rcutoff取r4与r5平均值。
步骤(2):通过孔喉半径与T2信号值的关系即公式(2)求取对应的T2截止值,通过大于T2截止值信号幅值面积与总信号幅值面积比值即可计算可动流体饱和度。
图1离心实验准确测定致密油中束缚水饱和度适用性示意图(A:大孔型样品含水饱和度与离心力关系;B:大孔与小孔并存型样品离心样品含水饱和度与离心力关系,C:大孔型样品孔隙分布,D:大孔与小孔并存型样品样品孔隙分布。Sw:含水饱和度;Swi:样品实际束缚水饱和度;Swi’:离心实验测定束缚水饱和度;P:离心实验最大离心力;R:孔喉半径;P1:大孔型样品可动水离心出所需要的最大离心力;P2:大孔与小孔并存型样品小孔中可动水离心出所需要的最小离心力;P3:大孔与小孔并存型样品小孔中可动水离心出所需要的最大离心力;);
图5为本发明实施例中,确定致密砂岩储层聚集孔喉半径下限测定方法原理图(r为相应孔喉半径)
图6本发明方法计算得到的T2截止值与现有技术计算得到的T2截止值的结果比较(符号注释:T2cutoff本发明方法计算得到的T2截止值;T2’cutoff:前人计算T2截止值)。
发明效果:
真实性:本方法创新性利用实际样品测定的致密油可动流体孔喉半径下限值求取T2截止值,由于实际样品测定的致密油可动流体孔喉半径下限值采用由粗变细孔喉测定碳含量逐渐逼近的方法,其数值真实反映了油在吼道中实际能够充注的最小半径值,即可动孔喉半径下限值,因此本方法相比离心法等其他方法更加具有真实性。
准确性:相比目前进行调研的其他几种求取T2截止值的方法,本方法具有准确性。离心实验中的离心力需要足够大才能够将小孔中的可动水离心出,但是目前离心实验技术中的离心力往往不能够达到相应要求,最重要的是,在实际离心实验中很难判定是否已经完全将可动水离心出岩心,只能通过判断可动水随着离心力继续增大不再明显减少就停止离心,而这种判断明显不具备科学性。本方法通过准确测定致密油可动流体孔喉半径下限值,带入公式即可准确求取T2截止值,进而求取准确的可动流动饱和度。
具体实施方式:
吉林油田扶余油层致密油R59井S121进行实验为例。致密油可动流体孔喉半径下限值实验测定数据如下:1号点:吼道半径110nm,碳含量36%;2号点:吼道半径88nm,碳含量24%;3号点:吼道半径76nm,碳含量15%;4号点:吼道半径42nm,碳含量7%;5号点:吼道半径18nm,碳含量0%。取42nm与18nm平均值30nm作为致密油可动流体孔喉半径下限值。核磁T2谱岩心弛豫率取15μm/s,岩心弛豫率一般采用经验值。砂岩的弛豫率经验值一般取值15μm/s。
孔喉半径下限值为30nm,则相应的T2截止值为1.5ms(图6)。
前人针对T2截止值进行了大量计算,致密储层的T2截止值一般取30ms。本方法计算的得到的1.5ms数据相比前人研究更加准确,特别是针对于致密储层,如果采用30ms,则可动流体百分数为大于30ms的信号幅值面积813与总信号幅值面积4739的比值,为17%;本方法计算的得到的1.5ms对应的可动流体百分数为大于1.5ms的信号幅值面积2902与总信号幅值面积4739的比值为61%,其中确定幅值面积即为信号值的积分。由于该井属于2013年新井,笔者选取了邻近井位的产液情况,计算得到的可动流体饱和度为58%。本方法采用的T2截止值计算的可动流体饱和度更加准确。
如上所述,得到对应的核磁T2弛豫时间截止值,之后通过信号幅值面积比值即可计算相应可动流体饱和度,因此本发明还公开了一种可动流体饱和度确定方法,方法包括:
确定待测样品的T2截止值;
根据确定的T2截止值的信号幅值面积与总信号面积比值确定可动流体饱和度,其中确定待测样品的T2截止值包括:
对待测样品进行碳含量测定,生成碳含量测定结果;
根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值。
以及一种可动流体饱和度确定装置,装置包括:
截止值确定模块,用于确定待测样品的T2截止值;
饱和度确定模块,用于根据确定的T2截止值的信号幅值面积与总信号面积比值确定可动流体饱和度,其中截止值确定模块包括:
测定模块,用于对待测样品进行碳含量测定,生成碳含量测定结果;
孔喉半径下限值确定模块,根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
T2截止值确定模块,用于根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种致密油核磁共振T2截止值确定方法,其特征在于,所述的方法包括:
对待测样品进行碳含量测定,生成碳含量测定结果;
根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值;
所述的根据所述碳含量测定结果确定待测样品的孔喉半径下限值包括:
根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
2.如权利要求1所述的方法,其特征在于,所述的根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值包括:
根据所述孔喉半径下限值和待测样品的岩心弛豫率利用公式r=2ρ2T2确定对应的T2截止值,其中,
ρ2:核磁共振信号T2的表面弛豫率,μm/s。
3.一种致密油核磁共振T2截止值确定装置,其特征在于,所述的装置包括:
测定模块,用于对待测样品进行碳含量测定,生成碳含量测定结果;
孔喉半径下限值确定模块,用于根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
T2截止值确定模块,用于根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值;
所述的孔喉半径下限值确定模块包括:
孔喉半径确定单元,用于根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
下限值确定单元,用于根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
4.如权利要求3所述的装置,其特征在于,所述的T2截止值确定模块根据所述孔喉半径下限值和待测样品的岩心弛豫率利用公式r=2ρ2T2确定对应的T2截止值,其中,
ρ2:核磁共振信号T2的表面弛豫率,μm/s。
5.一种可动流体饱和度确定方法,其特征在于,所述的方法包括:
确定待测样品的T2截止值;
根据确定的T2截止值的信号幅值面积与总信号面积比值确定可动流体饱和度,其中确定待测样品的T2截止值包括:
对待测样品进行碳含量测定,生成碳含量测定结果;
根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值;
所述的根据所述碳含量测定结果确定待测样品的孔喉半径下限值包括:
根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
6.如权利要求5所述的方法,其特征在于,所述的根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值包括:
根据所述孔喉半径下限值和待测样品的岩心弛豫率利用公式r=2ρ2T2确定对应的T2截止值,其中,
ρ2:核磁共振信号T2的表面弛豫率,μm/s。
7.一种可动流体饱和度确定装置,其特征在于,所述的装置包括:
截止值确定模块,用于确定待测样品的T2截止值;
饱和度确定模块,用于根据确定的T2截止值的信号幅值面积与总信号面积比值确定可动流体饱和度,其中截止值确定模块包括:
测定模块,用于对待测样品进行碳含量测定,生成碳含量测定结果;
孔喉半径下限值确定模块,根据所述碳含量测定结果确定待测样品的孔喉半径下限值;
T2截止值确定模块,用于根据所述孔喉半径下限值和待测样品的岩心弛豫率确定对应的T2截止值;
所述的孔喉半径下限值确定模块包括:
孔喉半径确定单元,用于根据碳含量测定结果,确定碳含量为零的测定点和相邻的碳含量非零的测定点分别对应的孔喉半径;
下限值确定单元,用于根据碳含量为零的测定点和相邻的碳含量非零的测定点所对应的两孔喉半径的平均值确定待测样品的孔喉半径下限值。
8.如权利要求7所述的装置,其特征在于,所述的T2截止值确定模块根据所述孔喉半径下限值和待测样品的岩心弛豫率利用公式r=2ρ2T2确定对应的T2截止值,其中,
ρ2:核磁共振信号T2的表面弛豫率,μm/s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410855862.2A CN104777181B (zh) | 2014-12-31 | 2014-12-31 | 致密油核磁共振t2截止值及流体饱和度确定方法、装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410855862.2A CN104777181B (zh) | 2014-12-31 | 2014-12-31 | 致密油核磁共振t2截止值及流体饱和度确定方法、装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104777181A CN104777181A (zh) | 2015-07-15 |
CN104777181B true CN104777181B (zh) | 2017-05-10 |
Family
ID=53618783
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410855862.2A Active CN104777181B (zh) | 2014-12-31 | 2014-12-31 | 致密油核磁共振t2截止值及流体饱和度确定方法、装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104777181B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106290440A (zh) * | 2015-05-26 | 2017-01-04 | 中国石油天然气股份有限公司 | 驱油效果检测方法及装置 |
CN105114064B (zh) * | 2015-08-04 | 2018-03-13 | 中国石油天然气股份有限公司 | 确定致密砂岩储层饱和度的方法 |
CN105093340A (zh) * | 2015-08-05 | 2015-11-25 | 中国石油天然气股份有限公司 | 煤层可动水体积的获取方法及获取装置 |
CN107290259B (zh) * | 2016-03-30 | 2019-09-10 | 中国石油化工股份有限公司 | 低渗砂岩储层有效渗流孔喉半径的计算方法 |
CN105891250B (zh) * | 2016-03-30 | 2017-12-19 | 中国石油天然气股份有限公司 | 致密砂岩储层原始含水饱和度的确定方法 |
CN107843611B (zh) * | 2016-09-20 | 2020-07-14 | 中国石油化工股份有限公司 | 低渗透砂岩储层可动流体饱和度核磁共振参数表征新方法 |
CN107228934A (zh) * | 2017-06-27 | 2017-10-03 | 中国石油大学(华东) | 致密砂岩油气充注孔喉半径下限的确定方法 |
US10338014B2 (en) * | 2017-07-27 | 2019-07-02 | Saudi Arabian Oil Company | Estimating formation properties using saturation profiles |
CN107748126B (zh) * | 2017-10-19 | 2019-10-11 | 北京青檬艾柯科技有限公司 | 一种获取岩石孔隙尺寸和孔隙表面弛豫率的核磁共振方法 |
CN110161071A (zh) * | 2019-04-24 | 2019-08-23 | 西安石油大学 | 一种评价致密储层可动流体最小喉道半径的方法 |
CN110044946B (zh) * | 2019-05-22 | 2022-04-15 | 张善文 | 一种低场核磁共振测量炭黑含油率的方法 |
CN110687612B (zh) | 2019-09-17 | 2020-09-08 | 中国石油天然气股份有限公司 | 吸附油和游离油含量连续表征的页岩油分析方法及装置 |
CN110823780B (zh) * | 2019-11-19 | 2021-05-18 | 中南大学 | 一种饱和冻结岩石未冻水含量的计算方法 |
CN113253354B (zh) * | 2021-05-12 | 2024-06-14 | 中国石油天然气集团有限公司 | 一种核磁t2谱截止值的确定方法、装置、介质及电子设备 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2248561C1 (ru) * | 2003-12-15 | 2005-03-20 | Общество с ограниченной ответственностью "ПермНИПИнефть" | Способ определения смачиваемости поровой поверхности пород-коллекторов |
CA2620669A1 (en) * | 2005-09-21 | 2007-03-29 | Idemitsu Kosan Co., Ltd. | Method of determining cyclic carbon content in sample and method of determining content of cyclic carbon material in sample |
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN102706913A (zh) * | 2012-06-19 | 2012-10-03 | 中国石油天然气股份有限公司 | 一种致密砂岩储层成藏孔喉半径下限的测定方法与装置 |
WO2012123882A3 (en) * | 2011-03-15 | 2012-12-27 | Schlumberger Technology B.V. | Examination of porosity by nmr and intrusion porosimetry |
CN103257151A (zh) * | 2013-04-24 | 2013-08-21 | 西安石油大学 | 一种定量评价油气二次运移过程中孔喉动用规律的方法 |
CN103454408A (zh) * | 2013-08-19 | 2013-12-18 | 中国石油天然气股份有限公司 | 一种测算致密砂岩油聚集量的方法与装置 |
CN103674803A (zh) * | 2013-11-18 | 2014-03-26 | 中国石油天然气集团公司 | 一种计算煤岩割理渗透率的方法及*** |
CN103743891A (zh) * | 2013-12-31 | 2014-04-23 | 中国石油天然气股份有限公司 | 致密油储层微纳米孔喉***中可动油赋存状态的测定方法 |
-
2014
- 2014-12-31 CN CN201410855862.2A patent/CN104777181B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2248561C1 (ru) * | 2003-12-15 | 2005-03-20 | Общество с ограниченной ответственностью "ПермНИПИнефть" | Способ определения смачиваемости поровой поверхности пород-коллекторов |
CA2620669A1 (en) * | 2005-09-21 | 2007-03-29 | Idemitsu Kosan Co., Ltd. | Method of determining cyclic carbon content in sample and method of determining content of cyclic carbon material in sample |
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
WO2012123882A3 (en) * | 2011-03-15 | 2012-12-27 | Schlumberger Technology B.V. | Examination of porosity by nmr and intrusion porosimetry |
CN102706913A (zh) * | 2012-06-19 | 2012-10-03 | 中国石油天然气股份有限公司 | 一种致密砂岩储层成藏孔喉半径下限的测定方法与装置 |
CN103257151A (zh) * | 2013-04-24 | 2013-08-21 | 西安石油大学 | 一种定量评价油气二次运移过程中孔喉动用规律的方法 |
CN103454408A (zh) * | 2013-08-19 | 2013-12-18 | 中国石油天然气股份有限公司 | 一种测算致密砂岩油聚集量的方法与装置 |
CN103674803A (zh) * | 2013-11-18 | 2014-03-26 | 中国石油天然气集团公司 | 一种计算煤岩割理渗透率的方法及*** |
CN103743891A (zh) * | 2013-12-31 | 2014-04-23 | 中国石油天然气股份有限公司 | 致密油储层微纳米孔喉***中可动油赋存状态的测定方法 |
Non-Patent Citations (1)
Title |
---|
储层流动孔喉下限研究;时宇等;《西南石油大学学报(自然科学版)》;20091231;第31卷(第6期);第42页第1.2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN104777181A (zh) | 2015-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104777181B (zh) | 致密油核磁共振t2截止值及流体饱和度确定方法、装置 | |
CN103884633B (zh) | 一种确定岩石渗透率的方法及装置 | |
CN104568694B (zh) | 一种致密岩心气‑水相对渗透率的测试方法 | |
US10444171B2 (en) | Absolute porosity and pore size determination of pore types in media with varying pore sizes | |
CN105464654B (zh) | 确定致密砂岩饱和度指数的方法及装置 | |
CN100349013C (zh) | 核磁共振测井t2谱t2截止值的确定方法 | |
CN106050225B (zh) | 一种核磁共振测井t2谱100%纯水谱的确定方法 | |
CN104237957B (zh) | 一种核磁共振测井t2截止值的确定方法 | |
CN106153662A (zh) | 岩心应力敏感性的测量方法 | |
CN203658217U (zh) | 一种渗流启动压力梯度测试实验装置 | |
CN104729972B (zh) | 确定碳酸盐岩气水相对渗透率的方法及装置 | |
CN105891089B (zh) | 储层渗透率确定方法及装置 | |
CN103257151B (zh) | 一种定量评价油气二次运移过程中孔喉动用规律的方法 | |
CN110296931B (zh) | 一种致密砂岩油水相对渗透率信息的表征方法及*** | |
MX2008008999A (es) | Pelicula de capas multiples resistente a la intemperie. | |
CN104101562A (zh) | 一种确定储层渗透率的方法 | |
CN105804734A (zh) | 一种利用核磁共振测井识别稠油储层的方法 | |
Meng et al. | Rock fabric of tight sandstone and its influence on irreducible water saturation in Eastern Ordos Basin | |
CN104698020B (zh) | 一种松散岩心微观孔隙结构特征参数的采集处理方法 | |
CN107907911A (zh) | 基于核磁共振的致密储层含油量测定方法 | |
MX2007008961A (es) | Metodo para analizar datos que tienen propiedades compartidas y distintas. | |
CN109580689B (zh) | 一种核磁共振测井t2截止值的逐点计算方法 | |
CN104297130B (zh) | 低渗透率储层的固有渗透率解释方法及*** | |
CN106644879B (zh) | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 | |
CN105401937B (zh) | 一种基于孔隙结构的饱和度指数预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |