CN112051292A - 一种孔喉结构动态表征方法及装置 - Google Patents

一种孔喉结构动态表征方法及装置 Download PDF

Info

Publication number
CN112051292A
CN112051292A CN202010862095.3A CN202010862095A CN112051292A CN 112051292 A CN112051292 A CN 112051292A CN 202010862095 A CN202010862095 A CN 202010862095A CN 112051292 A CN112051292 A CN 112051292A
Authority
CN
China
Prior art keywords
core sample
static
throat
curve
dynamic
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.)
Pending
Application number
CN202010862095.3A
Other languages
English (en)
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.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Publication of CN112051292A publication Critical patent/CN112051292A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/28Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/44Sample treatment involving radiation, e.g. heat
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N5/00Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
    • G01N5/04Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by removing a component, e.g. by evaporation, and weighing the remainder

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本说明书实施例提供一种孔喉结构动态表征方法及装置。所述方法包括:在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线;基于所述静态核磁共振T2谱曲线和静态喉道分布曲线计算转换系数;在至少一个测量压力下测量对应于所述岩心样本的动态核磁共振T2谱曲线;根据所述动态核磁共振T2谱曲线和所述转换系数确定动态喉道分布曲线;所述动态喉道分布曲线,用于表征所述岩心样本在压力动态变化下的孔喉结构。通过上述方法,在求取得到核磁共振T2谱曲线和喉道分布曲线之间的转换系数之后,能够将对应于岩心样本的不同压力下的核磁共振T2谱曲线转化为相应的喉道分布曲线,进而对所述岩心样本的动态孔喉结构进行表征。

Description

一种孔喉结构动态表征方法及装置
技术领域
本说明书实施例涉及地质勘探开发技术领域,特别涉及一种孔喉结构动态表征方法及装置。
背景技术
孔喉结构指地层中的孔隙喉道所对应的结构,例如可以是地层中连接不同储集区域的狭窄通道。孔喉结构是影响储层储存能力和渗流能力的关键因素,掌握储层的孔喉结构对于储层中油气资源的勘探开发具有重要的指导意义。
在实际生产中,往往会在对储层进行开采之前探明储层所对应的孔喉结构。但随着地层中生产活动的进行,地层中的压力也在不断变化,进而导致地层的孔喉结构也会随之变化。开采之前所探明的孔喉结构可能会与实际生产过程中的孔喉结构具有一定偏差,进而对生产过程中的决策造成影响。因此,目前亟需一种能够获取压力变动情况下地层所对应的孔喉结构的方法。
发明内容
本说明书实施例的目的是提供一种孔喉结构动态表征方法及装置,以解决如何获取不同压力下地层所对应的孔喉结构的问题。
为了解决上述技术问题,本说明书实施例提出一种孔喉结构动态表征方法,包括:
在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线;
基于所述静态核磁共振T2谱曲线和静态喉道分布曲线计算转换系数;
在至少一个测量压力下测量对应于所述岩心样本的动态核磁共振T2谱曲线;
根据所述动态核磁共振T2谱曲线和所述转换系数确定动态喉道分布曲线;所述动态喉道分布曲线,用于表征所述岩心样本在压力动态变化下的孔喉结构。
本说明书实施例还提出一种孔喉结构动态表征装置,包括:
静态曲线获取模块,用于在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道半径曲线;
系数计算模块,用于基于所述静态核磁共振T2谱曲线和静态喉道半径曲线计算转换系数;
T2谱曲线获取模块,用于在至少一个测量压力下测量对应于所述岩心样本的动态核磁共振T2谱曲线;
喉道分布曲线获取模块,用于根据所述动态核磁共振T2谱曲线和所述转换系数确定动态喉道分布曲线;所述动态喉道分布曲线,用于表征所述岩心样本在压力动态变化下的孔喉结构。
由以上本说明书实施例提供的技术方案可见,本说明书实施例在固定压力下获取到对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线,计算得到能够将这两种曲线相互进行转换的转换系数,在利用核磁共振获取在不同的测量压力下对应于岩心样本的动态核磁共振T2谱曲线,进而可以利用所述转换系数将所述动态核磁共振T2谱曲线转换为相应的动态喉道分布曲线,从而获取到所述岩心样本在压力动态变化的情况下孔喉结构的动态变化情况。上述方法利用核磁共振的方法对不同压力下的岩心样本进行测量,避免了测量本身对岩心样本所造成的影响,保证了测量的准确性,在对所述动态核磁共振T2谱曲线进行转换后,所得到的动态喉道分布曲线也能够准确地对岩心样本的孔喉结构进行表征。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本说明书实施例一种孔喉结构动态表征方法的流程图;
图2为本说明书实施例一种孔喉结构动态表征方法的示意图;
图3为本说明书实施例一种孔喉结构动态表征装置的模块图。
具体实施方式
下面将结合本说明书实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本说明书一部分实施例,而不是全部的实施例。基于本说明书中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本说明书保护的范围。
为了解决上述技术问题,本说明书实施例提出了一种孔喉结构动态表征方法。如图1所示,所述孔喉结构动态表征方法包括以下步骤。
S110:在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线。
固定压力指具有某一固定压力值的压力。在维持压力固定不变的情况下,防止压力变化对获取静态核磁共振T2谱曲线和静态喉道分布曲线所造成的干扰,提高测量结果的准确性。在一些实施方式中,所述固定压力可以指对应于所述岩心样本的围压和孔隙压力都分别维持一个特定的压力值。
岩心样本是针对待勘探或待开发的地层进行钻取所得到的样本。所述岩心样本的孔喉结构可以进一步地反映对应地层的孔喉结构,进而对实际应用做出指导。
所述静态核磁共振T2谱曲线是所述岩心样本在所述固定压力下由核磁共振实验下所得到的实验数据所得到的曲线。基于核磁共振理论,所述静态核磁共振T2谱曲线可以用于反映所述岩心样本的孔隙度,间接地体现出所述岩心样本的孔喉结构。具体的采取核磁共振实验获取所述静态核磁共振T2谱曲线的过程可以基于实际情况由操作人员进行实施,在此不做赘述。
所述静态喉道分布曲线是所述岩心样本在所述固定压力下通过高压压汞实验所得到的曲线。由于汞对于绝大多数岩石都是非湿润相,在向岩石中注汞时,若向汞所施加的压力大于或等于孔隙喉道的毛管压力,汞就能克服毛管阻力进入孔隙中。根据进汞的孔隙体积分数和对应的压力,就能够得到相应的高压压汞曲线,进而可以将所述高压压汞曲线转化为喉道分布曲线。具体的采取高压压汞实验获取所述静态喉道分布曲线的过程可以基于实际情况由操作人员进行实施,在此不做赘述。
但是,在进行高压压汞实验时,由于进汞饱和度一般会小于100%,高压压汞实验所得到的喉道半径曲线并不能反映岩心样本中的全部孔隙结构,因此利用高压压汞实验虽然能够直接得到喉道分布情况,但并不能完全反映所述岩心样本的孔喉结构。此外,高压压汞实验会对岩心样本的结构造成一定的影响,在针对同一岩心样本进行多次高压压汞实验时并不能保证后续实验结果的准确性。因此,可以结合核磁共振实验和高压压汞实验来确定所述岩心样本的孔喉结构。
在一些实施方式中,由于获取静态喉道分布曲线时所采取的高压压汞实验会对所述岩心样本的孔喉结构造成一定的影响,因此,为了保证所述静态核磁共振T2谱曲线和静态喉道分布曲线能够相对应,可以先获取对应于所述岩心样本的静态核磁共振T2谱曲线,再获取所述静态喉道分布曲线,从而保证实验结果的准确性。
在一些实施方式中,可以在获取对应于所述岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线之前,对所述岩心样本进行预处理。预处理的目的是为了在后续过程中能够更好的对所述岩心样本进行测量。
具体的,对所述岩心样本进行预处理可以包括:切割所述岩心样本至预设处理尺寸;烘干所述岩心样本至恒重;将所述岩心样本抽真空后注入模拟地层水至所述岩心样本饱和。
在切割所述岩心样本时,可以按照预先所设定的形状及尺寸对所述岩心样本进行切割,这样一来,在重复实施所述孔喉结构动态表征方法时,不会因为所述岩心样本的形状或尺寸不同而对实验结果产生影响。此外,在切割所述岩心样本时,为了克服所述岩心样本表面不均匀所带来的影响,还可以磨平所述岩心样本的端面。
在切割完成后,还可以对所述岩心样本的尺寸进行测量,利用一个具体的示例进行说明,例如可以将所述岩心样本切割为直径为25mm,长度为70mm的圆柱。在将所述岩心样本切割为圆柱的情况下时,可以对所述岩心样本的直径和长度进行测量,测量得到的数据可以用于后续步骤中的计算。具体的,可以利用游标卡尺对所述岩心样本进行测量,以确保测量结果精确至0.01mm。
在对所述岩心样本切割完成后,可以对所述岩心样本进行烘干处理。具体的,可以是将所述岩心样本烘干至重量不再变化的程度,例如,相邻的两次称量所相差的差值小于预设差值阈值时,可以确定所述岩心样本的烘干操作完成。在一个具体的示例中,可以将所述岩心样本放置于104℃的环境中并保持烘干时间在12小时以上,在烘干12小时之后每间隔一小时进行一次称量,当相邻的两次称量的差值小于10mg时完成所述烘干操作。烘干操作能够便于后续步骤中在所述岩心样本中注入模拟地层水并饱和的操作。
为了防止向所述岩心样本注水时由水敏性造成影响,需要配置对应于所述岩心样本的模拟地层水,例如,具有特定比例的生理盐水。在对所述岩心样本进行抽真空后,可以将所述模拟地层水注入所述岩心样本中,以使所述岩心样本具有饱和的模拟地层水。在一个具体的示例中,对所述岩心样本抽真空的时间可以在十二小时以上,饱和模拟地层水的时间可以在十二小时以上,且饱和的压力在20MPa以上。
在一些实施方式中,获取所述静态核磁共振T2谱曲线和静态喉道分布曲线之前,还可以将所述岩心样本划分为第一岩心样本和第二岩心样本,并在固定压力下获取对应于所述第一岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线,在后续步骤中获取对应于所述第二岩心样本的动态核磁共振T2谱曲线。由于在获取静态喉道分布曲线时所采用的高压压汞实验会对岩心样本的结构造成影响,若在后续步骤中仍然获取该岩心样本在不同的测量压力下的动态核磁共振T2谱曲线,无法保证所述岩心样本的孔喉结构是完全由于压力的变化而引起改变的,进而对实验结果造成一定的影响。因此,将所述岩心样本划分为第一岩心样本和第二岩心样本并分别进行实验能够保证实验的精确性,从而实现准确地对岩心样本的动态孔喉结构进行表征。
S120:基于所述静态核磁共振T2谱曲线和静态喉道分布曲线计算转换系数。
由于所述静态核磁共振T2谱曲线和静态喉道分布曲线都可以用于对岩心样本的孔隙结构进行一定程度上的反映,因此所述静态核磁共振T2谱曲线和静态喉道分布曲线之间具有一定的关联性。可以通过所述静态核磁共振T2谱曲线和静态喉道分布曲线之间的相似程度计算出相应的转换系数。
所述转换系数可以用于表示静态核磁共振T2谱曲线中的核磁共振T2值与静态喉道分布曲线中对应的喉道半径值之间的大小比例关系。由于静态核磁共振T2谱曲线和静态喉道分布曲线具有相似的变化趋势,使得两种曲线在各处具有近似的幅值转换比例系数,因此,利用所述转换系数可以实现所述静态核磁共振T2谱曲线和静态喉道分布曲线之间的相互转化。
所述转换系数与所述静态核磁共振T2谱曲线和静态喉道分布曲线之间例如可以满足如公式T2=Cr所对应的关系,式中,T2为静态核磁共振T2谱曲线中的T2值,单位为ms,C为转换系数,单位为ms/μm,r为静态喉道分布曲线中的喉道半径值,单位为μm。
具体的,计算所述转换系数的过程可以是先分别获取所述静态核磁共振T2谱曲线中的T2值和喉道分布曲线中的孔道半径值,再利用最小二乘法计算所述T2值和所述孔道半径值之间的误差,在结合所述误差的基础上拟合所述T2值和孔道半径值得到转换系数。
由于实际应用中核磁共振实验与高压压汞实验的精确度之间存在一定的差异,因此所述静态核磁共振T2谱曲线中的T2值与所述喉道分布曲线中的孔道半径值可能并不能完全对应。因此,在计算转换系数之前还需要在所述核磁共振T2谱曲线和喉道分布曲线中选取恰当的区间来选取T2值和孔道半径值进行计算,以提高转换系数的准确程度。具体的,例如可以将所述静态喉道分布曲线和所述静态核磁共振T2谱曲线按照由大孔向小孔进行累积的方式,得到累积曲线,在所述累积曲线上选取合理的区间并将该区间上的喉道半径与相同饱和度下的T2值进行比较。由于两种曲线上的数据并不能完全一一对应,因此可以利用等间距纵向插值法计算两种方法相同累积频率下的T2值和孔道半径值,再利用最小二乘法计算所述T2值和所述孔道半径值之间的误差,基于所述误差对所述T2值和孔道半径值进行拟合进而可以得到所述转换系数。在考虑所述核磁共振T2谱曲线和喉道分布曲线的对应关系以及彼此之间的误差的情况下对转换系数的计算能够进一步地提高所述计算结果的准确性,从而能够在后续步骤中更为精确地获取所述岩心样本的动态孔喉结构。
S130:在至少一个测量压力下测量对应于所述岩心样本的动态核磁共振T2谱曲线。
为了测量所述岩心样本的动态核磁共振T2谱曲线,可以改变所述岩心样本对应的压力,并使得所述压力符合岩心样本对应的地层在实际开采过程中的压力变化情况,进而在后续步骤中确定所述岩心样本在压力动态变化的情况下的孔喉结构变动情况。因此,所述至少一个测量压力可以是对应于实际生产过程中的压力变动情况的具有不同压力值的压力。
由于在实际生产中主要是地层的孔隙压力变动而引起对应的孔喉结构变动,因此,在一些实施方式中,为了使实验结果更为符合孔喉结构动态变化时的压力变化特点,不同的测量压力可以指对应于岩心样本的不同的孔隙压力。
由于核磁共振实验对岩心样本基本不会造成结构上的影响,因此,利用核磁共振实验对不同测量压力下的岩心样本进行测量时可以保证测量结果的准确性,避免由于实验本身对岩心样本的结构造成影响。
当测量所述动态核磁共振T2谱曲线时,可以在某一测量压力下对所述岩心样本执行核磁共振实验,并测量在该测量压力下的核磁共振T2谱曲线;接着将对应于所述岩心样本的压力调整至另一测量压力对应的压力值,再次执行核磁共振实验以获取调整后的测量压力所对应的核磁共振T2谱曲线。重复执行上述步骤直至测量得到对应于所有测量压力的核磁共振T2谱曲线,作为动态核磁共振T2谱曲线。
具体的利用核磁共振实验获取所述岩心样本的动态核磁共振T2谱曲线的操作可以基于实际情况由操作人员进行实施,在此不做赘述。
在一些实施方式中,可以在保证所述岩心样本的围压不变的情况下,改变对应于所述岩心样本的孔隙压力,进而获取对应于所述岩心样本的动态核磁共振T2谱曲线。由于在实际生产开发中地层中的孔隙结构变化主要是由于孔隙压力降低而造成了相应的孔隙结构的变化,因此在保证岩心样本的围压不变的情况下对孔隙压力进行调整更符合实际情况中的变化,从而使结果更为明显、合理。
具体的,在调整所述岩心样本的围压的过程中,可以保证围压和孔隙压力同步增加,并保持围压和孔隙压力之间的压力差不大于预设压力差阈值,从而维持岩心样本的结构的稳定性。通过上述操作将所述岩心样本的围压和孔隙压力均升至模拟地层压力后,在维持所述围压为模拟地层压力的情况下,逐步调整所述孔隙压力,例如可以是在每次调整孔隙压力之后将所述岩心样本静止一段固定的时间,以使所述岩心样本的孔喉结构变化完全,再对所述岩心样本的核磁共振T2谱曲线进行测量。重复上述步骤并统计所有不同的孔隙压力下对应于所述岩心样本的核磁共振T2谱曲线作为动态核磁共振T2谱曲线。
下面结合附图2对获取所述动态核磁共振T2谱曲线的场景示例进行介绍。基于附图2中的设备,首先将核磁共振设备开机并等待12小时以上,以使岩石夹持器内的温度保持稳定,之后再将所述岩心样本放置于核磁共振设备机箱的岩心夹持器内,利用真空泵将所述岩心样本抽真空后,利用注射泵将模拟地层水注入岩心后,使岩心完全饱和,具体的,所述模拟地层水可以为与所述岩心样本所对应的地层内的地层水相对应的合成盐水。接着可以利用手摇泵将所述岩心样本的围压逐步升至模拟地层压力,在此过程中,利用注射泵将岩心夹持器内的孔压也逐步升为模拟地层压力,并保证所述围压和所述孔隙压力的差值始终不大于5MPa。待所述围压和孔隙压力稳定12小时之后再测量在所述模拟地层压力下的核磁共振T2谱曲线。测量完毕后,通过调节注射泵及回压阀逐步降低岩心孔隙压力,降压完毕后等待至少8小时以使得所述岩心样本的孔喉结构变化完全,且岩心内的流速稳定,之后再测量所述岩心样本对应于降低后的孔隙压力的核磁共振T2谱曲线。重复上述步骤以获取所述动态核磁共振T2谱曲线,并传输至数据收集系统中以在后续步骤中获取动态喉道分布曲线。
需要说明的是,在上述场景示例中用于对所述岩心样本施加围压的氟油不含氢信号,使得所述氟油不会对核磁共振实验造成影响。此外,为了保证实验的顺利进行,所述岩心夹持器所能承受的最大工作压力需要高于所述岩心样本所对应的研究区内的地层压力。
S140:根据所述动态核磁共振T2谱曲线和所述转换系数确定动态喉道分布曲线;所述动态喉道分布曲线,用于表征所述岩心样本在压力动态变化下的孔喉结构。
动态喉道分布曲线即为所述岩心样本分别对应于多个测量压力的喉道分布曲线。由于喉道分布曲线可以对所述岩心样本的孔喉结构进行定量直观的描述,因此,通过所述动态喉道分布曲线可以表征所述岩心样本在压力动态变化下的孔喉结构,具体的,可以用于模拟目标地层在围压不变,孔隙压力变化下的孔喉结构变动情况,进行对实际勘探开发提供参考价值。
在得到所述动态核磁共振T2谱曲线之后,基于步骤S120中求取得到的转换系数,可以将所述动态核磁共振T2谱曲线转换为动态喉道分布曲线。
具体的,可以是利用公式r=T2/C将所述动态核磁共振T2谱曲线中的T2值转换为喉道半径,式中,r为喉道半径,T2为T2值,C为转换系数。在获取到所述喉道半径后,根据计算得到的喉道半径可以构建相应的喉道分布曲线。
在一些实施方式中,所述孔喉结构动态表征方法还可以包括获取所述岩心样本的渗透率和孔隙度。渗透率可以表征岩石中允许流体通过的能力,孔隙度可以表征岩石中孔隙体积占岩石体积的百分比。因此,所述渗透率和孔隙度也可以在结合所述动态喉道分布曲线的情况下对所述岩心样本在压力动态变化下的孔喉结构。
具体的,获取所述孔隙度可以是在获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线之前气测所述岩心样本的孔隙度。具体的气测孔隙度的操作步骤可以基于实际情况由操作人员进行实施,在此不再赘述。
所述岩心样本的渗透率可以分为气测渗透率和液测渗透率。对于所述气测渗透率的获取可以是在获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线之前气测所述岩心样本的渗透率。所述气测渗透率的具体操作方法可以参照获取气测孔隙度的具体步骤,在此不再赘述。
所述液测渗透率可以是在获取动态核磁共振T2谱曲线时,利用在不同压力下岩心内的流体流速对所述岩心样本的渗透率进行测量。具体的,例如可以是利用公式
Figure BDA0002648480230000081
计算所述岩心样本的液测渗透率,式中,K为液测渗透率,Q为岩心样本内的模拟地层水的流量,μ为模拟地层水粘度,L为岩心样本的长度,A为岩心样本的横截面积,ΔP为岩心样本入口压力。
综上,能够综合所获取到的孔隙度、渗透率和动态喉道半径分布曲线对所述岩心样本的孔喉结构进行较为准确的动态表征。
通过本说明书实施例中所介绍的孔喉结构动态表征方法,基于核磁共振实验对岩心样本不会造成损害的特点,获取在各个测量压力下对应于岩心样本的核磁共振T2谱曲线,再利用所述核磁共振T2谱曲线与喉道分布曲线之间的转换系数将所述核磁共振T2谱曲线转换为动态喉道分布曲线,进而定量且直接地对所述岩心样本在不同的测量压力下的孔喉结构进行表征,避免了测量本身对岩心样本所造成的影响,保证了测量的准确性。基于所得到的孔喉结构的动态变化情况,对实际生产中的储层变动情况提供参考价值,进而对实际生产提供帮助。
基于上述孔喉结构动态表征方法,本说明书还提出一种孔喉结构动态表征装置的实施例。如图3所示,所述孔喉结构动态表征装置具体包括:
静态曲线获取模块310,用于在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道半径曲线;
系数计算模块320,用于基于所述静态核磁共振T2谱曲线和静态喉道半径曲线计算转换系数;
T2谱曲线获取模块330,用于在至少一个测量压力下测量对应于所述岩心样本的动态核磁共振T2谱曲线;
喉道分布曲线获取模块340,用于根据所述动态核磁共振T2谱曲线和所述转换系数确定动态喉道分布曲线;所述动态喉道分布曲线,用于表征所述岩心样本在压力动态变化下的孔喉结构。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的,计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本说明书可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本说明书的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本说明书各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本说明书可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本说明书,本领域普通技术人员知道,本说明书有许多变形和变化而不脱离本说明书的精神,希望所附的权利要求包括这些变形和变化而不脱离本说明书的精神。

Claims (10)

1.一种孔喉结构动态表征方法,其特征在于,包括:
在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线;
基于所述静态核磁共振T2谱曲线和静态喉道分布曲线计算转换系数;
在至少一个测量压力下测量对应于所述岩心样本的动态核磁共振T2谱曲线;
根据所述动态核磁共振T2谱曲线和所述转换系数确定动态喉道分布曲线;所述动态喉道分布曲线,用于表征所述岩心样本在压力动态变化下的孔喉结构。
2.如权利要求1所述的方法,其特征在于,所述在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线之前,还包括:
对所述岩心样本进行预处理;其中,包括:切割所述岩心样本至预设处理尺寸;烘干所述岩心样本至恒重;将所述岩心样本抽真空后注入模拟地层水至所述岩心样本饱和。
3.如权利要求1所述的方法,其特征在于,所述在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线,包括:
绘制对应于岩心样本的高压压汞曲线;
利用所述高压压汞曲线获取静态喉道分布曲线。
4.如权利要求1所述的方法,其特征在于,所述在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线之前,还包括:
将所述岩心样本划分为第一岩心样本和第二岩心样本;
相应的,所述在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线,包括:
在固定压力下获取对应于所述第一岩心样本的静态核磁共振T2谱曲线和静态喉道分布曲线;
所述在至少一个测量压力下测量对应于所述岩心样本的动态核磁共振T2谱曲线,包括:
在至少一个测量压力下测量对应于所述第二岩心样本的动态核磁共振T2谱曲线。
5.如权利要求1所述的方法,其特征在于,所述基于所述静态核磁共振T2谱曲线和静态喉道分布曲线计算转换系数,包括:
分别获取所述静态核磁共振T2谱曲线中的T2值和静态喉道分布曲线中的孔道半径值;
利用最小二乘法计算所述T2值和所述孔道半径值之间的误差;
结合所述误差拟合所述T2值和孔道半径值得到转换系数。
6.如权利要求1所述的方法,其特征在于,所述在至少一个测量压力下测量对应于所述岩心样本的动态核磁共振T2谱曲线,包括:
在预设固定围压下测量所述岩心样本对应于至少一个孔隙压力的动态核磁共振T2谱曲线。
7.如权利要求1所述的方法,其特征在于,所述根据所述动态核磁共振T2谱曲线和所述转换系数确定动态喉道分布曲线,包括:
利用公式r=T2/C将所述动态核磁共振T2谱曲线中的T2值转换为喉道半径,式中,r为喉道半径,T2为T2值,C为转换系数;
根据所述喉道半径构建对应于所述至少一个测量压力的喉道分布曲线。
8.如权利要求1所述的方法,其特征在于,所述方法还包括:
获取所述岩心样本的渗透率和孔隙度;所述渗透率和孔隙度用于结合所述动态喉道分布曲线表征所述岩心样本在压力动态变化下的孔喉结构。
9.如权利要求8所述的方法,其特征在于,所述渗透率包括气测渗透率和液测渗透率;
相应的,所述基于所述静态核磁共振T2谱曲线和静态喉道分布曲线计算转换系数之前,还包括:
气测所述岩心样本的气测渗透率;
所述根据所述动态核磁共振T2谱曲线和所述转换系数确定动态喉道分布曲线之前,还包括:
利用公式
Figure FDA0002648480220000021
计算所述岩心样本的液测渗透率,式中,K为液测渗透率,Q为岩心样本内的模拟地层水的流量,μ为模拟地层水粘度,L为岩心样本的长度,A为岩心样本的横截面积,ΔP为岩心样本入口压力。
10.一种孔喉结构动态表征装置,其特征在于,包括:
静态曲线获取模块,用于在固定压力下获取对应于岩心样本的静态核磁共振T2谱曲线和静态喉道半径曲线;
系数计算模块,用于基于所述静态核磁共振T2谱曲线和静态喉道半径曲线计算转换系数;
T2谱曲线获取模块,用于在至少一个测量压力下测量对应于所述岩心样本的动态核磁共振T2谱曲线;
喉道分布曲线获取模块,用于根据所述动态核磁共振T2谱曲线和所述转换系数确定动态喉道分布曲线;所述动态喉道分布曲线,用于表征所述岩心样本在压力动态变化下的孔喉结构。
CN202010862095.3A 2020-07-27 2020-08-25 一种孔喉结构动态表征方法及装置 Pending CN112051292A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2020107304107 2020-07-27
CN202010730410 2020-07-27

Publications (1)

Publication Number Publication Date
CN112051292A true CN112051292A (zh) 2020-12-08

Family

ID=73600132

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010862095.3A Pending CN112051292A (zh) 2020-07-27 2020-08-25 一种孔喉结构动态表征方法及装置

Country Status (1)

Country Link
CN (1) CN112051292A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112816516A (zh) * 2021-02-07 2021-05-18 中国石油天然气集团有限公司 一种饱含甲烷岩样高温高压核磁共振t2谱实验室测量方法
CN112924356A (zh) * 2021-01-28 2021-06-08 中国石油大学(北京) 一种储层动态渗流特征获取方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103926267A (zh) * 2014-04-28 2014-07-16 西安石油大学 一种定量评价应力敏感过程中孔喉变化程度的方法
CN104634718A (zh) * 2015-03-05 2015-05-20 中国石油大学(华东) 应用核磁共振表征致密砂岩孔径分布的标定方法
CN106153662A (zh) * 2016-06-17 2016-11-23 北京大学 岩心应力敏感性的测量方法
US10161891B1 (en) * 2017-11-03 2018-12-25 Institute Of Geology And Geophysics Chinese Academy Of Sciences Method for characterizing rock physical characteristics of deeply buried carbonate rocks

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103926267A (zh) * 2014-04-28 2014-07-16 西安石油大学 一种定量评价应力敏感过程中孔喉变化程度的方法
CN104634718A (zh) * 2015-03-05 2015-05-20 中国石油大学(华东) 应用核磁共振表征致密砂岩孔径分布的标定方法
CN106153662A (zh) * 2016-06-17 2016-11-23 北京大学 岩心应力敏感性的测量方法
US10161891B1 (en) * 2017-11-03 2018-12-25 Institute Of Geology And Geophysics Chinese Academy Of Sciences Method for characterizing rock physical characteristics of deeply buried carbonate rocks

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
GUANGFENG LIU ET AL.: "Experimental determination of dynamic pore-throat structure characteristics in a tight gas sandstone formation with consideration of effective stress", 《MARINE AND PETROLEUM GEOLOGY》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112924356A (zh) * 2021-01-28 2021-06-08 中国石油大学(北京) 一种储层动态渗流特征获取方法及装置
CN112816516A (zh) * 2021-02-07 2021-05-18 中国石油天然气集团有限公司 一种饱含甲烷岩样高温高压核磁共振t2谱实验室测量方法

Similar Documents

Publication Publication Date Title
CN105891089B (zh) 储层渗透率确定方法及装置
CN112051292A (zh) 一种孔喉结构动态表征方法及装置
US10359542B2 (en) Generating dynamically calibrated geo-models in green fields
CN107688037B (zh) 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法
CN109932301A (zh) 一种计算致密储层自发渗吸两相流体相对渗透率的方法
CN106355571A (zh) 一种白云岩储层质量的确定方法及装置
CN110470581B (zh) 确定储层应力敏感程度的方法、装置及存储介质
CN109597964B (zh) 理论油藏渗流参数值确定方法、系统、设备、可读介质
CN112147689B (zh) 一种碳酸盐岩储层解释方法及解释量版建立方法
CN111208566B (zh) 基于sca模型的孔缝参数反演方法、装置及存储介质
McPhee et al. Capillary pressure
CN111707695A (zh) 基于核磁共振数据确定储层渗透率的方法、装置及设备
CN112464587B (zh) 一种预测储层毛管压力曲线的方法、装置及设备
CN106837318B (zh) 岩层稠油含量获取方法和装置
CN113914853B (zh) 深水沉积体系砂泥岩薄互层砂岩储层含水饱和度确定方法
CN111485874A (zh) 确定水平井储层含水饱和度的方法及装置
CN117744302A (zh) 渗透率预测方法、装置及计算机可读存储介质
CN115271367A (zh) 高原隧道围岩分级方法、装置、设备以及存储介质
CN114441401B (zh) 覆压孔隙度的确定方法、装置、计算机设备及存储介质
CN111597671B (zh) 一种基于概率分布的裂缝网络复杂度确定方法及系统
CN110805425B (zh) 原油压缩系数计算方法、装置、设备及可读存储介质
CN112415598A (zh) 近地表地层品质因子模型的建立方法和装置
CN112924356A (zh) 一种储层动态渗流特征获取方法及装置
CN111751879A (zh) 储层孔隙度预测方法及装置
CN111856571B (zh) 测井频段速度向地震频段速度转换方法、装置及设备

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20201208