CN114089421B - 一种油气储层非均质性分析方法 - Google Patents

一种油气储层非均质性分析方法 Download PDF

Info

Publication number
CN114089421B
CN114089421B CN202111478167.5A CN202111478167A CN114089421B CN 114089421 B CN114089421 B CN 114089421B CN 202111478167 A CN202111478167 A CN 202111478167A CN 114089421 B CN114089421 B CN 114089421B
Authority
CN
China
Prior art keywords
heterogeneity
size
core
full
cores
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
CN202111478167.5A
Other languages
English (en)
Other versions
CN114089421A (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.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
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 Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202111478167.5A priority Critical patent/CN114089421B/zh
Publication of CN114089421A publication Critical patent/CN114089421A/zh
Application granted granted Critical
Publication of CN114089421B publication Critical patent/CN114089421B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48—Processing data
    • G01V1/50—Analysing data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Sampling And Sample Adjustment (AREA)

Abstract

本发明提供了一种油气储层非均质性分析方法，适用于油气地质领域，该方法包括：采集待评价区内目标油气储层的全尺寸岩心，对比成像测井与岩心特征确定岩心方位；轴向平分岩心为n块扇形体，并制取其轴切面的铸体薄片；获取每个铸体薄片的平均孔隙度Poi和平均孔隙半径roi，构建基于非均质性评价参数K1i＝Poi/Po1和K2i＝roi/ro1的非均质性指标函数，计算多深度点非均质性指标函数平均值；在按上述方法获取多口钻井非均质性指标函数的基础上，利用插值方法绘制其平面等值线图，以此评价油气储层非均质性。相对于现有技术，该方法步骤简单、实用且更为经济，可快速分析储层非均质性。

Description

一种油气储层非均质性分析方法
技术领域
本发明涉及一种油气储层非均质性分析方法，尤其适用于油气地质领域。
背景技术
储层非均质性是低孔隙度、低渗透率非常规油气藏开发时面临的重要问题。储层非均质性是指储层各种性质在空间变化的属性，主要表现为岩石物质组成的非均质性和孔隙空间的非均质性。当前，油气储层非均质性可以利用综合地质分析法、储层地质建模法和实验室分析法开展分析。
申请公布号CN104375204A的发明专利提出了一种分析储层非均质性的方法和装置，基于岩心刻度电成像测井资料和核磁共振测井资料，对储层进行宏观和微观非均质性分析；申请公布号CN107219159A的发明专利提出一种确定储层非均质性的方法和装置，依据待测储层各个深度组中所有地层内容，在测试渗透率基础上计算获取纵向非均质性指数，进而确定非均质性的评价结果；申请公布号CN112228050A的发明专利提出了一种定量评价致密油储层宏观非均质性的方法及其应用，通过获取目标致密油储层的夹层密度、砂地比、渗砂层表征参数、砂体结构指数、孔隙度及含油指数，建立储层非均质性综合评价模型；申请公布号CN112836393A的发明专利提出了一种基于多尺度熵分析储层非均质性的方法，其主要以测井资料为基础，以归一化和粗粒化的处理方式，将测井信号根据不同尺度得到多个不同信息，构建多尺度熵曲线，实现储层非均质性评价。上述专利涉及的方法步骤较为繁琐，实用性、经济性有待提升。
发明内容
针对现有技术的不足之处，提供一种步骤简单、操作性强、实用性和经济性并存的油气储层非均质性分析方法。
为实现上述技术目的，本发明的一种油气储层非均质性分析方法，在待评价区获取多口钻井目标油气储层段的成像测井资料，以及每口钻井中与成像测井深度对应的多块全尺寸岩心，对比成像测井与岩心特征确定岩心的方位；将所有全尺寸岩心轴向平分岩心为n块扇形体，并制取其轴切面的铸体薄片；通过观测获取每个全尺寸岩心制作的铸体薄片，获取每个铸体薄片的平均孔隙度和平均孔隙半径信息，再利用每个铸体薄片的平均孔隙度和平均孔隙半径信息计算出非均质性评价参数，进而获得非均质性指标函数，计算每口钻井中目标油气储层段全尺寸岩心的非均质性指标函数平均值，利用插值方法绘制其平面等值线图，以此评价油气储层非均质性。
具体步骤：
步骤1，在待评价区内选择一口钻井，获取该钻井目标油气储层段的成像测井资料和与其对应深度的多块全尺寸岩心，由于取出的全尺寸岩心无法体现岩心方向，因此通过比对并匹配岩心与成像测井表现出的层理、裂缝产状特征，利用成像测井扫描图中的方向标注确定岩心方位，为非均质性分析提供岩心的空间方位基础；
步骤2，对确定了岩心方位的该钻井目标油气储层段全尺寸岩心沿其轴向平均切成n块大小形态一致的扇形体，以方位区间[0，360/n)的扇形体记为第1块，顺时针块号依次增大到n，取每块扇形体同侧轴切面制作成矩形铸体薄片，从而获得n片矩形铸体薄片；
步骤3，通过测量计算平均数的方法分别获取每片矩形铸体薄片的平均孔隙度Poi和平均孔隙半径roi信息，然后利用每片矩形铸体薄片的平均孔隙度Poi和平均孔隙半径roi信息计算非均质性评价参数K1i＝Poi/Po1和K2i＝roi/ro1，i＝1,2,···,n；
步骤4，利用非均质性评价参数计算非均质性指标函数：Zi＝K1i×K2i，i＝1,2,···,n，通过不同编号的矩形铸体薄片中记载的平均孔隙度Poi和平均孔隙半径roi判断该全尺寸岩心采集处不同方位区间非均质性指标函数值，重复步骤2和步骤3，获得k块全尺寸岩心的非均质性指标函数，求k块不同深度采集的全尺寸岩心的非均质性指标函数平均值
Figure BDA0003394398420000021
步骤5，按照上述步骤1～步骤4的分析方法，重复计算待评价区域不少于10口钻井的不同方位区间非均质性指标函数平均值，利用插值方法将这些非均质性指标函数平均值Zai绘制平面等值线图，通过该平面等值线图评价区域储层非均质性。
步骤2中已获取了岩心方位的全尺寸岩心沿其轴向平均切成n＝8块，岩心方位区间依次为0°～45°、45°～90°、90°～135°、135°～180°、180°～225°、225°～270°、270°～315°和315°～360°，块号依次为1-8，铸体薄片规格为矩形常规尺寸，大小为a×2cm，a为岩心半径。
步骤5中所用的插值方法优选Kriging法。
有益效果：
目前的储层非均质性评价方法步骤繁琐，实用性和经济性欠佳。鉴于此，本发明的目的在于提供一种步骤简单、操作性强、实用性和经济性并存的油气储层非均质性分析方法，该方法首先将确定了方向的全尺寸岩心沿其轴向平均切成n块大小和形态一致的扇形体，并取其同侧轴切面制作铸体薄片，计算不同铸体薄片的平均孔隙度和平均孔隙半径，计算非均质性评价参数，获取非均质性评价指标函数，计算每口钻井中目标油气储层段全尺寸岩心的非均质性指标函数平均值，利用插值方法绘制其平面等值线图，以此评价储层非均质性。
附图说明
图1为本发明一种油气储层非均质性分析方法的流程示意图；
图2为本发明中全尺寸岩心扇形体切割示意图。
具体实施方式
下面结合附图对本发明的实施例做进一步说明
如图1所示，本发明的一种油气储层非均质性分析方法，首先在待评价区获取多口钻井目标油气储层段的成像测井资料，以及每口钻井中与成像测井深度对应的多块全尺寸岩心，对比成像测井与岩心特征确定岩心的方位；将所有全尺寸岩心轴向平分岩心为n块扇形体，并制取其轴切面的铸体薄片；通过观测获取每个全尺寸岩心制作的铸体薄片，获取每个铸体薄片的平均孔隙度和平均孔隙半径信息，再利用每个铸体薄片的平均孔隙度和平均孔隙半径信息计算出非均质性评价参数，进而获得非均质性指标函数，计算每口钻井中目标油气储层段全尺寸岩心的非均质性指标函数平均值，利用插值方法绘制其平面等值线图，以此评价油气储层非均质性。
具体步骤：
步骤1，在待评价区内选择一口钻井，获取该钻井目标油气储层段的成像测井资料和与其对应深度的多块全尺寸岩心，由于取出的全尺寸岩心无法体现岩心方向，因此通过比对并匹配岩心与成像测井表现出的层理、裂缝产状特征，利用成像测井扫描图中的方向标注确定岩心方位，为非均质性分析提供岩心的空间方位基础；
步骤2，对确定了岩心方位的该钻井目标油气储层段全尺寸岩心沿其轴向平均切成n块大小形态一致的扇形体，以方位区间[0，360/n)的扇形体记为第1块，顺时针块号依次增大到n，取每块扇形体同侧轴切面制作成矩形铸体薄片，从而获得n片矩形铸体薄片；
步骤3，通过测量计算平均数的方法分别获取每片矩形铸体薄片的平均孔隙度Poi和平均孔隙半径roi信息，然后利用每片矩形铸体薄片的平均孔隙度Poi和平均孔隙半径roi信息计算非均质性评价参数K1i＝Poi/Po1和K2i＝roi/ro1，i＝1,2,···,n；
步骤4，利用非均质性评价参数计算非均质性指标函数：Zi＝K1i×K2i，i＝1,2,···,n，通过不同编号的矩形铸体薄片中记载的平均孔隙度Poi和平均孔隙半径roi判断该全尺寸岩心采集处不同方位区间非均质性指标函数值，重复步骤2和步骤3，获得k块全尺寸岩心的非均质性指标函数，求k块不同深度采集的全尺寸岩心的非均质性指标函数平均值
Figure BDA0003394398420000031
步骤5，按照上述步骤1～步骤4的分析方法，重复计算待评价区域不少于10口钻井的不同方位区间非均质性指标函数平均值，利用Kriging插值方法将这些非均质性指标函数平均值Zai绘制平面等值线图，通过该平面等值线图评价区域储层非均质性。
实施例一
步骤2中已获取了岩心方位的全尺寸岩心沿其轴向平均切成n＝8块，岩心方位区间依次为0°～45°、45°～90°、90°～135°、135°～180°、180°～225°、225°～270°、270°～315°和315°～360°，块号依次为1-8，铸体薄片规格为矩形常规尺寸，大小为a×2cm，a为岩心半径。

Claims (3)

1.一种油气储层非均质性分析方法，其特征在于：在待评价区获取多口钻井目标油气储层段的成像测井资料，以及每口钻井中与成像测井深度对应的多块全尺寸岩心，对比成像测井与岩心特征确定岩心的方位；将所有全尺寸岩心轴向平分岩心为n块扇形体，并制取其轴切面的铸体薄片；通过观测获取每个全尺寸岩心制作的铸体薄片，获取每个铸体薄片的平均孔隙度和平均孔隙半径信息，再利用每个铸体薄片的平均孔隙度和平均孔隙半径信息计算出非均质性评价参数，进而获得非均质性指标函数，计算每口钻井中目标油气储层段全尺寸岩心的非均质性指标函数平均值，利用插值方法绘制其平面等值线图，以此评价油气储层非均质性；
具体步骤为：
步骤1，在待评价区内选择一口钻井，获取该钻井目标油气储层段的成像测井资料和与其对应深度的多块全尺寸岩心，由于取出的全尺寸岩心无法体现岩心方向，因此通过比对并匹配岩心与成像测井表现出的层理、裂缝产状特征，利用成像测井扫描图中的方向标注确定岩心方位，为非均质性分析提供岩心的空间方位基础；
步骤2，对确定了岩心方位的该钻井目标油气储层段全尺寸岩心沿其轴向平均切成n块大小形态一致的扇形体，以方位区间[0，360/n)的扇形体记为第1块，顺时针块号依次增大到n，取每块扇形体同侧轴切面制作成矩形铸体薄片，从而获得n片矩形铸体薄片；
步骤3，通过测量计算平均数的方法分别获取每片矩形铸体薄片的平均孔隙度Poi和平均孔隙半径roi信息，然后利用每片矩形铸体薄片的平均孔隙度Poi和平均孔隙半径roi信息计算非均质性评价参数K1i＝Poi/Po1和K2i＝roi/ro1，i＝1,2,···,n；
步骤4，利用非均质性评价参数计算非均质性指标函数：Zi＝K1i×K2i，i＝1,2,···,n，通过不同编号的矩形铸体薄片中记载的平均孔隙度Poi和平均孔隙半径roi判断该全尺寸岩心采集处不同方位区间非均质性指标函数值，重复步骤2和步骤3，获得k块全尺寸岩心的非均质性指标函数，求k块不同深度采集的全尺寸岩心的非均质性指标函数平均值
Figure FDA0004053738740000011
步骤5，按照上述步骤1～步骤4的分析方法，重复计算待评价区域不少于10口钻井的不同方位区间非均质性指标函数平均值，利用插值方法将这些非均质性指标函数平均值Zai绘制平面等值线图，通过该平面等值线图评价区域储层非均质性。
2.按照权利要求1所述的一种油气储层非均质性分析方法的步骤，其特征在于：步骤2中已获取了岩心方位的全尺寸岩心沿其轴向平均切成n＝8块，岩心方位区间依次为0°～45°、45°～90°、90°～135°、135°～180°、180°～225°、225°～270°、270°～315°和315°～360°，块号依次为1-8，铸体薄片规格为矩形常规尺寸，大小为a×2cm，a为岩心半径。
3.按照权利要求1所述的一种油气储层非均质性分析方法的步骤，其特征在于：步骤5中所用的插值方法为Kriging法。
CN202111478167.5A 2021-12-06 2021-12-06 一种油气储层非均质性分析方法 Active CN114089421B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111478167.5A CN114089421B (zh) 2021-12-06 2021-12-06 一种油气储层非均质性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111478167.5A CN114089421B (zh) 2021-12-06 2021-12-06 一种油气储层非均质性分析方法

Publications (2)

Publication Number Publication Date
CN114089421A CN114089421A (zh) 2022-02-25
CN114089421B true CN114089421B (zh) 2023-03-21

Family

ID=80306641

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111478167.5A Active CN114089421B (zh) 2021-12-06 2021-12-06 一种油气储层非均质性分析方法

Country Status (1)

Country Link
CN (1) CN114089421B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112965114A (zh) * 2021-02-07 2021-06-15 中国石油大学(华东) 一种海上深层天然气储层甜点评价方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103278436B (zh) * 2013-02-01 2015-02-04 西安石油大学 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法
CN103353462B (zh) * 2013-06-17 2015-12-09 中国石油大学(华东) 一种基于核磁共振成像的岩石非均质性定量评价方法
CN105447298B (zh) * 2014-09-29 2018-09-11 中国石油化工股份有限公司 一种储层分析的方法及其应用
CN104375204A (zh) * 2014-11-21 2015-02-25 中国海洋石油总公司 一种分析储层非均质性的方法和装置
CN106093350B (zh) * 2016-08-10 2018-05-04 中国石油天然气股份有限公司 确定非均质碳酸盐岩储层饱和度指数的方法
CN107219159B (zh) * 2017-04-19 2019-07-09 中国石油天然气股份有限公司 一种确定储层非均质性的方法和装置
CN107807407B (zh) * 2017-09-30 2019-10-11 中国石油天然气股份有限公司 一种油气区带有效性评价方法和装置
CN107727679B (zh) * 2017-11-03 2018-08-14 中国科学院地质与地球物理研究所 一种表征深层碳酸盐岩岩石物理学特征方法
CN109697752B (zh) * 2018-11-20 2023-07-25 中国石油天然气集团有限公司 基于岩心ct图像孔隙信息提取定量表征岩心非均质性方法
FR3094742B1 (fr) * 2019-04-04 2021-04-16 Ifp Energies Now Procédé pour exploiter un réservoir pétrolier fracturé ayant une taille de pores hétérogène
CN110702580B (zh) * 2019-10-21 2022-02-11 西安石油大学 一种基于信息熵的致密砂岩储层孔喉非均质性表征方法
CN112228050B (zh) * 2020-10-22 2023-12-26 中国石油天然气股份有限公司 一种定量评价致密油储层宏观非均质性的方法及其应用
CN113435066B (zh) * 2021-08-26 2021-11-12 北京润泽创新科技有限公司 一种基于数字岩心技术的测井解释储层评价方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112965114A (zh) * 2021-02-07 2021-06-15 中国石油大学(华东) 一种海上深层天然气储层甜点评价方法

Also Published As

Publication number Publication date
CN114089421A (zh) 2022-02-25

Similar Documents

Publication Publication Date Title
AU2011258594B2 (en) Method for obtaining consistent and integrated physical properties of porous media
CN105115874B (zh) 多源信息融合的多组分三维数字岩心构建方法
CN110644980B (zh) 一种超低渗透油藏储层综合分类评价方法
AU2009316427B2 (en) Method for determining rock physics relationships using computer tomographic images thereof
EP2359311B1 (en) Method for determining elastic-wave attenuation of rock formations using computer tomograpic images thereof
CN108694264B (zh) 一种确定页岩气储层渗透率的方法
CN108444881B (zh) 一种适用于陆相泥页岩微纳米尺度储集空间的表征方法
CN109375283B (zh) 一种砂岩储层3d渗透率演化史的分析方法
CN112145165B (zh) 一种微裂缝-孔隙型储层动静态渗透率转换方法
CN111963159A (zh) 一种砾岩致密油储层流体性质的识别方法
CN115030707A (zh) 一种油页岩“甜点”快速评价方法
CN106777649B (zh) 一种裂缝型储集层孔隙结构定量评价方法
CN113216945B (zh) 一种致密砂岩储层渗透性定量评价方法
CN114089421B (zh) 一种油气储层非均质性分析方法
CN108802192B (zh) 一种钙屑砂岩储层孔隙类型识别方法
CN112012727B (zh) 获得气相有效渗透率的方法及储层产能的预测方法
WO2022204297A1 (en) The combined discrete gaussian analysis of micp and nmr t2 distributions of multi-modal carbonate rocks
CN113530533B (zh) 基于三品质参数的甜点区预测方法及系统
CN103675945A (zh) 一种测定孔洞型储层的饱和度的方法及设备
CN109580678B (zh) 利用数字岩心技术快速识别评价低阻油气层方法
CN114991745B (zh) 一种页岩油储层甜点识别方法和装置
CN116502557B (zh) 基于网格分析法的储层压裂后可动流体饱和度计算方法
CN116956643B (zh) 一种烃源岩纵向/横向非均质特征定量表征方法
Beg Multiscale Pore Network Modeling of Hierarchical Media with Applications to Improved Oil and Gas Recovery
CN114441402B (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