CN107153074B - 一种基于高精度成像的微观剩余油定量评价方法 - Google Patents

一种基于高精度成像的微观剩余油定量评价方法 Download PDF

Info

Publication number
CN107153074B
CN107153074B CN201710168433.1A CN201710168433A CN107153074B CN 107153074 B CN107153074 B CN 107153074B CN 201710168433 A CN201710168433 A CN 201710168433A CN 107153074 B CN107153074 B CN 107153074B
Authority
CN
China
Prior art keywords
remaining oil
microscopic
data
displacement
dimensional
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
Application number
CN201710168433.1A
Other languages
English (en)
Other versions
CN107153074A (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.)
Oceanographic Instrumentation Research Institute Shandong Academy of Sciences
Original Assignee
Oceanographic Instrumentation Research Institute Shandong Academy of Sciences
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 Oceanographic Instrumentation Research Institute Shandong Academy of Sciences filed Critical Oceanographic Instrumentation Research Institute Shandong Academy of Sciences
Priority to CN201710168433.1A priority Critical patent/CN107153074B/zh
Publication of CN107153074A publication Critical patent/CN107153074A/zh
Application granted granted Critical
Publication of CN107153074B publication Critical patent/CN107153074B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Geology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Fluid Mechanics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Pulmonology (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明属于微观剩余油定量表征技术领域,公开了一种基于高精度成像的微观剩余油定量评价方法,主要包括以下步骤:用高精度扫描设备扫描岩心,获取并分析包含岩心结构信息的三维数字图像数据体;根据实际需求设计并开展微观驱替实验;获取不同驱替状态下流体分布数据;获得剩余油三维分布数据体;获得剩余油参数;定量评价微观剩余油分布规律。本发明可定性观测岩心孔隙结构和流体分布变化、定量评价岩心中的微观剩余油分布变化,分析不同开采方式、不同开采阶段下微观剩余油的分布规律,因地制宜提高采收率。

Description

一种基于高精度成像的微观剩余油定量评价方法
技术领域
本发明属于微观剩余油定量表征技术领域,具体涉及一种基于高精度成像的微观剩余油定量评价方法。
背景技术
油田开发的最终目的就是如何从地层孔隙中尽可能多的采出原油,在国内的老油田大部分进入中后期油田开发的情况下,研究油气储层岩石孔隙中的剩余油在宏观层面和微观层面的分布规律有助于我们了解油藏目前的开发情况。因此,从微观、小尺度微纳级别研究剩余油的分布,开展基于微观驱替实验的微观剩余油定量评价,深入到孔隙空间内部展开分析,了解和掌握剩余油的微观形成机理和分布规律,进而建立宏观与微观的关系,为后期提高采收率的措施找到理论上的依据。
在微观层面上,储层岩石孔隙空间结构的区别是带来剩余油的分布状态巨变的主要原因。与常规的仿真模型和二维视图研究方法不同,目前,借助高精度成像设备以及图像处理技术,真实三维再现每一类岩心中,每一种驱替方式对其中剩余油分布的影响,并开展基于拓扑结构的定量评价,是为孔隙尺度上的提高采收率(EOR)方法研究提供实验基础的重要手段。
吼道:连接各个孔隙空间之间的狭窄部分。
发明内容
本发明为克服现有技术的不足,采用下述技术方案予以实现:
一种基于高精度成像的微观剩余油定量评价方法具体包括以下步骤:
S1:用高精度成像设备获取样品的三维数字图像;
S2:设计微观驱替实验,获取不同驱替状态下流体分布数据;
S21:根据需求设计驱替实验,通过不同的实验条件模拟不同的驱替方式。
S22:利用碘化试剂标定驱替液相或者剩余油,;基于扫描设备,开展微观驱替实验,在相应的观测点,扫描获取岩心三维流体分布数据。
S3:定量分析剩余油三维分布数据体;
S31:获得剩余油三维数据体,将S2的岩心三维流体分布数据进行数据二值化,分离目标区域对应的所有剩余油点,标记为目标点,并赋值1,其余背景区域(驱替相、骨架等)标记为背景点或噪声点,赋值0,得到只包含0,1的二值数据体;
S32:获得欧式距离图,利用三维空间欧式距离变换方法,对步骤S51的二值数据体进行三维距离测度处理,每一个目标体素点都将被赋值为该点与最近的边界处的距离,边界点则被赋值为0,得到数据体的欧氏距离图,
S33:基于S32中得到的欧氏距离图确定扫描顺序,利用基于三维空间邻居关系的标签方法,按顺序扫描识别,依次以不同数字标定剩余油数字图像中相互连通的剩余油簇,刻画此分辨率下的剩余油簇三维分布;
S34:基于S33中得到的剩余油簇三维体积分布,定义剩余油簇微观赋存形态,以剩余油簇体积与欧拉数作为判定标准,统计此分辨率下的剩余油赋存形态变化;
S35:分割孔隙空间,基于S32的欧氏距离图,利用分水岭算法,计算距离图像中的“盆地”或是“山峰”,将各个孔隙空间区别开,计算每个孔隙空间及吼道空间的最大内切球半径获得孔隙半径和吼道半径,以及剩余油在每个孔隙空间中的赋存体积。
基于S34的剩余油赋存形态和S35的孔隙空间分割数据,利用三维空间拓扑结构方法相互匹配,计算各剩余油基本参数,包括饱和度、体积、表面积、欧拉数、连通程度、赋存形态、赋存位置等。
S4:定量评价微观剩余油分布规律。
优选地,高精度扫描设备为CT。
选取岩心样品,基于CT设备开展不同条件下的微观驱替实验,扫描驱替过程中岩心内的流体分布图像,定性观测岩心孔隙结构和流体分布变化,通过图像二值化、欧式距离测度、分水岭分割方法以及三维图像形态拓扑匹配方法定量分析岩心中剩余油赋存形态以及分布位置变化,得到剩余油饱和度变化曲线、不同赋存状态的剩余油相对含量、微观剩余油饱和度分布曲线等。
本发明根据不同研究需求,基于CT仪器开展不同实验条件下的微观驱替实验,数字化驱替实验过程以及岩心中的流体分布,依据拓扑方法对数据体开展分割和匹配,最终获取微观剩余油分布数据,定量分析、评价相应的微观剩余油分布规律。本技术方案可定性观测岩心孔隙结构和流体分布变化、定量评价岩心中的微观剩余油分布变化,分析不同开采方式、不同开采阶段下微观剩余油的分布规律,因地制宜提高采收率。
附图说明
图1:本发明一种基于高精度成像的微观剩余油定量评价方法流程示意图;
图2:本发明实施例对砂岩开展水驱后的剩余油三维分布图,分别为饱和油阶段,注入1PV、5PV、15PV、50PV(残余油阶段);
图3:本发明实施例砂岩孔隙空间分割图;
图4:本发明实施例砂岩水驱过程中不同剩余油赋存形态;
图5:本发明实施例砂岩开展水驱后的微观剩余油赋存形态变化,分别为饱和油阶段,注入1PV、5PV、15PV、50PV(残余油阶段);
图6:本发明实施例开展水驱后原始饱和油阶段微观剩余油分布曲线图;
图7:本发明实施例开展水驱后残余油阶段微观剩余油分布曲线图。
具体实施方式
下面结合附图和具体实施方式对本发明的技术方案作进一步详细的说明。
如图1-7所示:选取砂岩为样品,基于CT设备开展不同条件下的微观驱替实验,扫描驱替过程中砂岩内的流体分布图像,定性观测砂岩孔隙结构和流体分布变化,通过图像二值化、欧式距离测度、分水岭分割方法以及三维图像形态拓扑匹配方法定量分析砂岩中剩余油赋存形态以及分布位置变化,得到剩余油饱和度变化曲线、不同赋存状态的剩余油相对含量、微观剩余油饱和度分布曲线等,得到图2-图7。
具体包括以下步骤:
S1:用高精度扫描设备扫描砂岩岩心,获取包含岩心结构信息的三维数字图像数据体,如图3,其中浅灰色为孔隙空间,判断此扫描参数下可获取剩余油在孔隙中的清晰分布图像;
S2:设计并开展微观驱替实验,获取不同驱替状态下流体分布数据;
根据需求设计砂岩微观水驱油实验,模拟常温时,水驱对剩余油分布的影响,分析高含水阶段时的剩余油微观分布规律。利用碘化试剂配制驱替油相与水相,加入18%的碘化试剂标定水相,测试其在扫描设备中的射线吸收情况;基于扫描设备,开展微观驱替实验,在注入体积分别为1PV,5PV,10PV等的观测点,扫描获取岩心三维数据。
S3:获取砂岩中的微观剩余油参数
S31:CT数据二值化,获得剩余油三维分布数据体;
S32:获得欧式距离图,利用基于三维空间邻居关系的欧式距离变换方法,对步骤S1和步骤S2的数据体进行三维距离测度处理,得到欧氏距离图;
S33:基于S32的欧氏距离图,利用基于三维空间邻居关系的标签方法,依次以不同数字标定剩余油数字图像中相互连通的剩余油,刻画此分辨率下的剩余油体积分布;
S34:基于S33的剩余油体积分布,利用欧拉数理论,将剩余油赋存分布形态归纳为六种如图4,并统计分析每种赋存形态在不同驱替时刻下的分布变化规律如图5。
S35:分割孔隙空间,基于S32的欧氏距离图,利用分水岭算法,将各个孔隙空间及孔隙空间中的剩余油分割开,如图3,计算每个孔隙空间及吼道空间的最大内切球半径获得孔隙半径和吼道半径,以及剩余油在每个孔隙中的赋存体积,统计得到剩余油分布位置图如图6-7;
S4:获得剩余油参数,结合剩余油形态、位置分布定量评价微观剩余油分布规律。可获知本发明实施例砂岩原始含油饱和度阶段几乎全部为连片状剩余油,其中连通性最好的剩余油簇几乎占到全部剩余油体积的99%,剩余油主要分布在较大孔隙中,在半径<25μm的小孔隙、小孔喉中分布较少。驱替开始后此大块剩余油油簇很快断裂、转化为其他赋存形态。水驱15PV后,部分剩余油以网络状剩余油仍然滞留在孔隙中,且其剩余油体积没有发生很大变化,仅吼道状剩余油在长期水驱作用下被开采出来。水驱结束后,剩余油以孤立状、吼道状、薄膜状的剩余油形态残留在孔隙半径约35μm的孔隙中,但由于岩心水驱油过程中很快形成水的优势通道,部分孔喉比较高处的大孔隙中(孔隙半径约50-60μm)也有剩余油富集。针对此类剩余油分布规律,可以在开发过程中通过注采调控技术来进一步改善开发效果,例如利用不稳定注水等水动力学方法,形成不稳定压力场,使中高含水期原来由于绕流和阻塞等被捕集在半径约45-60μm的大孔隙中的簇状剩余油,能够重新动用。
实施例仅说明本发明的技术方案,而非对其进行任何限制;尽管参照前述实施例对本发明进行了详细的说明,对于本领域的普通技术人员来说,依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或替换,并不使相应技术方案的本质脱离本发明所要求保护的技术方案的精神和范围。

Claims (4)

1.一种基于高精度成像的微观剩余油定量评价方法,其特征在于,包括以下步骤:
S1:用高精度成像设备获取三维数字图像;
S2:设计微观驱替实验,获取不同驱替状态下流体分布数据;
S3:定量分析剩余油分布数据,获得剩余油参数;
S4:定量评价微观剩余油分布规律;
所述S2根据需求设计驱替实验,通过不同的实验条件模拟不同的驱替方式以及开采阶段;所述S2中的微观驱替实验利用碘化试剂标定驱替液相或者剩余油;基于成像设备,开展微观驱替实验,在相应的观测点,扫描获取岩心流体分布三维数据;
所述S3包括以下步骤:
S31:获得剩余油三维数据体,将S2的三维流体分布图像数据进行数据二值化,分离目标区域对应的所有剩余油点,标记为目标点,并赋值1,其余背景区域标记为背景点,赋值0,得到只包含0,1的二值数据体;
S32:获得欧式距离图,利用三维空间欧式距离变换方法,对步骤S31的二值数据体进行三维距离测度处理,并对步骤S1的数据体进行三维距离测度处理,每一个目标体素点都将被赋值为该点与最近的边界处的距离,边界点则被赋值为0,得到数据体的欧氏距离图;
S33:基于S32中得到的欧氏距离图确定扫描顺序,利用基于三维空间邻居域关系的标签方法,按顺序扫描识别,依次以不同数字标定剩余油数字图像中相互连通的剩余油簇,刻画此分辨率下的剩余油簇三维体积分布;
S34:基于S33中得到的剩余油簇三维体积分布,定义剩余油微观赋存形态,以剩余油簇体积与欧拉数作为判定标准,统计此分辨率下的剩余油赋存形态变化;
S35:分割孔隙空间,基于S32的欧氏距离图,利用分水岭算法,计算距离图像中的“盆地”或是“山峰”,将各个孔隙空间区别开,计算每个孔隙空间及吼道空间的最大内切球半径获得孔隙半径和吼道半径,以及剩余油在每个孔隙空间中的赋存体积。
2.根据权利要求1所述的一种基于高精度成像的微观剩余油定量评价方法,其特征在于:S4是基于S34的剩余油赋存形态和S35的孔隙空间分割数据,利用三维空间拓扑结构方法相互匹配,计算各剩余油基本参数。
3.根据权利要求2所述的一种基于高精度成像的微观剩余油定量评价方法,其特征在于:所述的剩余油基本参数为饱和度、体积、表面积、欧拉数、连通程度、赋存形态、赋存位置。
4.根据权利要求1所述的一种基于高精度成像的微观剩余油定量评价方法,其特征在于:高精度成像设备为CT。
CN201710168433.1A 2017-03-21 2017-03-21 一种基于高精度成像的微观剩余油定量评价方法 Active CN107153074B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710168433.1A CN107153074B (zh) 2017-03-21 2017-03-21 一种基于高精度成像的微观剩余油定量评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710168433.1A CN107153074B (zh) 2017-03-21 2017-03-21 一种基于高精度成像的微观剩余油定量评价方法

Publications (2)

Publication Number Publication Date
CN107153074A CN107153074A (zh) 2017-09-12
CN107153074B true CN107153074B (zh) 2019-10-29

Family

ID=59791764

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710168433.1A Active CN107153074B (zh) 2017-03-21 2017-03-21 一种基于高精度成像的微观剩余油定量评价方法

Country Status (1)

Country Link
CN (1) CN107153074B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107957488B (zh) * 2017-10-27 2021-01-01 中国石油天然气股份有限公司 剩余油量的确定方法和装置
CN107975366B (zh) * 2017-11-24 2020-07-17 北京科技大学 一种堵水剂在多孔介质中封堵调剖可视化评价方法
CN108825222B (zh) * 2018-06-02 2021-08-03 东北石油大学 运用分形理论确定功能型聚合物驱剩余油分布及运移规律的方法
CN108843311B (zh) * 2018-06-02 2021-08-03 东北石油大学 运用分形理论确定水驱剩余油分布及运移规律的方法
CN108678738B (zh) * 2018-06-02 2021-11-12 东北石油大学 运用分形理论确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律的方法
CN108868715B (zh) * 2018-06-02 2020-05-01 东北石油大学 运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法
CN108952696B (zh) * 2018-06-02 2021-08-03 东北石油大学 运用分形理论确定化学驱剩余油分布及运移规律的方法
CN108843310B (zh) * 2018-06-02 2021-11-12 东北石油大学 运用分形理论确定基质-高渗条带水驱剩余油动用规律的方法
CN109934906B (zh) * 2018-08-16 2023-03-14 清能艾科(深圳)能源技术有限公司 原油饱和度的获取方法和装置、计算机设备、存储介质
CN110644979B (zh) * 2019-09-03 2021-08-27 中国石油大学(北京) 孔隙流体初始赋存状态的获取方法及其装置
CN110984971B (zh) * 2019-10-14 2023-06-02 中国石油天然气股份有限公司 一种基于形成机理分析的微观剩余油分类方法
CN112150424B (zh) * 2020-09-16 2023-02-24 中国石油大学(华东) 一种基于荧光薄片的微观剩余油定量分析方法
CN112945829B (zh) * 2021-02-07 2023-05-26 西安石油大学 一种致密砂岩储层水驱剩余油分析方法及系统
CN116152317B (zh) * 2023-04-24 2023-07-07 北京润泽创新科技有限公司 基于数字岩心技术的原位对比分析剩余油赋存特征方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104076046B (zh) * 2013-03-28 2016-08-31 中国石油化工股份有限公司 多孔介质中剩余油微观分布图像采集与定量表征方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104076046B (zh) * 2013-03-28 2016-08-31 中国石油化工股份有限公司 多孔介质中剩余油微观分布图像采集与定量表征方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种新的三维欧式距离变换方法及在数字岩心中的应用;王鑫等;《科学通报》;20161212;第62卷(第15期);第1662-1669页 *

Also Published As

Publication number Publication date
CN107153074A (zh) 2017-09-12

Similar Documents

Publication Publication Date Title
CN107153074B (zh) 一种基于高精度成像的微观剩余油定量评价方法
Yang et al. New pore space characterization method of shale matrix formation by considering organic and inorganic pores
CN101929973B (zh) 裂缝储层含油气饱和度定量计算方法
CN100429528C (zh) 一种沉积相控制预测砂岩油藏储层属性的方法
CN108802073B (zh) 一种基于数字岩心的岩石电性参数获取方法及装置
Wienhöfer et al. Predicting subsurface stormflow response of a forested hillslope–the role of connected flow paths
US9135378B2 (en) Method of developing a reservoir from a technique of selecting the positions of wells to be drilled
US9422800B2 (en) Method of developing a petroleum reservoir from a technique for selecting the positions of the wells to be drilled
Jeannin et al. Speleological investigations
CN102518421B (zh) 物理模拟可视化实验装置及其形成方法
CN105626023A (zh) 低渗透油藏垂直压裂裂缝方位试井确定方法
CN105759310A (zh) 复杂非均质储层介质中地震波衰减及速度频散预测方法
An et al. The microscale analysis of reverse displacement based on digital core
CN104298797A (zh) 一种确定缝洞型油藏高导流通道圈闭剩余油的方法
Hohenbrink et al. Texture-depending performance of an in situ method assessing deep seepage
CN110019594A (zh) 一种基于数值模拟的地质知识库建立方法
Feldmann et al. Near surface roughness estimation: A parameterization derived from artificial rainfall experiments and two-dimensional hydrodynamic modelling for multiple vegetation coverages
CN103345867A (zh) 承压井抽水动态实验仪
CN112903555A (zh) 考虑孔隙各向异性的多孔介质渗透率计算方法及装置
Yao et al. Upscaling of carbonate rocks from micropore scale to core scale
Salah Geostatistical analysis of groundwater levels in the south Al Jabal Al Akhdar area using GIS
Jagelke et al. Conceptualization and implementation of a regional groundwater model for the Neckar catchment in the framework of an integrated regional model
Romanov et al. Karstification of an aquifer along the Birs river, Switzerland—A modeling approach
Lončar et al. Hydraulic‐hydrology analysis of the turbulent seepage flow within karst aquifer of the Golubinka spring catchment
CN113568065A (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
GR01 Patent grant
GR01 Patent grant
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 266200, Qingdao, Shandong, Qingdao, Qingdao, the core of the blue Silicon Valley, blue Silicon Valley business center, phase one, building No. 1.

Applicant after: Inst. of Marine Apparatus & Instruments, Shandong Prov. Academy of Sciences

Address before: 266071 Shandong city of Qingdao province Zhejiang City Road No. 28

Applicant before: Inst. of Marine Apparatus & Instruments, Shandong Prov. Academy of Sciences