CN105654525B - 反映真实岩心连通性的2d模型构建方法 - Google Patents

反映真实岩心连通性的2d模型构建方法 Download PDF

Info

Publication number
CN105654525B
CN105654525B CN201410643887.6A CN201410643887A CN105654525B CN 105654525 B CN105654525 B CN 105654525B CN 201410643887 A CN201410643887 A CN 201410643887A CN 105654525 B CN105654525 B CN 105654525B
Authority
CN
China
Prior art keywords
auto
image
correlation function
model
pore
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
CN201410643887.6A
Other languages
English (en)
Other versions
CN105654525A (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 Petroleum and Chemical Corp
Sinopec Shengli Geological Scientific Reserch Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Shengli Geological Scientific Reserch Institute
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 Petroleum and Chemical Corp, Sinopec Shengli Geological Scientific Reserch Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201410643887.6A priority Critical patent/CN105654525B/zh
Publication of CN105654525A publication Critical patent/CN105654525A/zh
Application granted granted Critical
Publication of CN105654525B publication Critical patent/CN105654525B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种反映真实岩心连通性的2D模型构建方法，该方法包括：步骤1，进行岩心微观孔隙结构图像预处理；步骤2，建立自相关函数；步骤3，根据预处理后的岩心微观孔隙结构图像和自相关函数，构建图像自相关函数与距离曲线；步骤4，从平面上应用自相关函数确定具有岩心孔隙结构特征代表性的模型尺度标准，从而确定二维模型的平面范围；步骤5，确定纵向孔隙最大贯穿深度；步骤6，进行纵向孔隙结构映射；以及步骤7，去除孤立孔隙，建立反映真实岩心孔隙结构特征及连通性的模拟模型。该反映真实岩心连通性的2D模型构建方法保证了对岩心内部微观流动模拟的可靠性和准确性，为真实反映微观油水渗流机理提供合理的参考依据。

Description

反映真实岩心连通性的2D模型构建方法
技术领域
本发明涉及油田开发技术领域，特别是涉及到一种反映真实岩心连通性的2D模型构建方法。
背景技术
三维孔隙尺度微观数值模拟模型是基于真实岩心扫描图像，识别岩心内部复杂的孔喉空间，进行精细地网格剖分后构建的。由于孔喉空间边界形态复杂，网格剖分工作量大，网格模型规模大，后期开展数值模拟耗时长，无法满足常规微观数值模拟计算的需求。而且，三维孔隙模型的孔喉空间平面及纵向交叉连通，流体在孔隙中的剥蚀、分离、聚并等过程较多，很难直观的追踪流体的流动形态，分析剩余油微观分布的作用机制。
构建孔隙尺度二维模型，一方面可以直观地展示流体微观流动形态，快速地开展剩余油分布特征的机理分析；另一方面，剖分后网格模型规模小，计算量小，省时省力。因此，二维模型能否反映真实岩心孔隙结构特征及连通性是模型构建的关键技术。为此我们发明了一种新的反映真实岩心连通性的2D模型构建方法，解决了以上技术问题。
发明内容
本发明的目的是提供一种构建的岩心微观模型与真实的岩心内部孔隙空间相似度更高的反映真实岩心连通性的2D模型构建方法。
本发明的目的可通过如下技术措施来实现：反映真实岩心连通性的2D模型构建方法，该反映真实岩心连通性的2D模型构建方法包括：步骤1，进行岩心微观孔隙结构图像预处理；步骤2，建立自相关函数；步骤3，根据预处理后的岩心微观孔隙结构图像和自相关函数，构建图像自相关函数与距离曲线；步骤4，从平面上应用自相关函数确定具有岩心孔隙结构特征代表性的模型尺度标准，从而确定二维模型的平面范围；步骤5，确定纵向孔隙最大贯穿深度；步骤6，进行纵向孔隙结构映射；以及步骤7，去除孤立孔隙，建立反映真实岩心孔隙结构特征及连通性的数模模型。
本发明的目的还可通过如下技术措施来实现：
在步骤1中，选取一张扫描图像作为模型构建对象，构建区域由众多岩石骨架及孔隙空间组成；根据图像灰度值，对孔隙和骨架进行分割处理，处理后图像变为(0,1)的二值图像。
在步骤2中，首先定义指示函数f与模拟区域中任意位置r的关系，用指示函数f(r)来表示：
其中，f(r)为指示函数，无量纲；r表示二值图像的空间位置；
由统计平均得到，孔隙度用指示函数表示构建两点自相关函数该函数用指示函数表示根据研究需要初步确定微观模型构建原点及范围。
在步骤3中，在图像模拟区域内选取任意一点作为原点，由步骤2得到，原点处自相关函数值统计为φ，无限远处，自相关函数值为φ2，做自相关函数与距离的函数曲线。
在步骤4中，由自相关函数与距离曲线可以得到，原点处，自相关函数为φ；距原点距离较远时，自相关函数值为φ2；距原点N个像素，是自相关函数变化的拐点，大于N个像素的距离时，自相关函数变化较平缓，说明N个像素内，图像的非均质程度已经得到最大程度地体现，根据此尺度标准，选取微观模型大于N个像素的图像区域。
在步骤5中，应用孔隙半径分布曲线，选取累积分布频率为95％对应的孔隙大小，该孔隙直径作为纵向孔隙的最大贯穿深度，该深度除以纵向扫描间隔即为可参考的纵向映射图像数。
在步骤6中，在纵向上，将多层扫描图像纵向叠加，层层映射，得到反映底层孔隙连通性的用二维图像。
在步骤7中，去除对流体流动不起作用的孤立孔隙，得到孔隙尺度数模模型的平面图像，根据该图像建立反映真实岩心孔隙结构特征及连通性的数模模型。
本发明中的反映真实岩心连通性的2D模型构建方法，涉及孔隙尺度数值模拟应用及油气田开发提高油藏采收率领域。具体地说是针对岩心扫描图像，通过自相关函数及孔隙大小分布曲线从平面及纵向选取能代表真实岩心孔隙结构特征及连通性的模型范围，建立孔隙尺度模拟模型，为开展岩心微观流动模拟提供基础。应用该方法可以从平面上确定构建微观数值模拟模型的尺度标准，纵向上，兼顾孔隙的纵向深度，将底层孔隙连通性进行映射，修正顶层平面孔隙连通性，从而得到孔隙尺度数模模型的图像基础，根据该图像即可建立反映真实岩心孔隙结构特征及连通性的数模模型。
本发明结合数学统计自相关函数，提出了构建微观模型的最小平面尺度范围，避免了选取最大范围扫描图像建模而带来的数值模拟的较大工作量，同时，与常规的依据单一扫描图像建立微观模型相比，由于本方法考虑了孔隙在纵向上的贯穿深度，平面连通性更好，构建的岩心微观模型与真实的岩心内部孔隙空间相似度更高，从而保证了对岩心内部微观流动模拟的可靠性和准确性，为真实反映微观油水渗流机理提供合理的参考依据。
附图说明
图1为本发明的反映真实岩心连通性的2D模型构建方法的一具体实施例的流程图；
图2是孔隙骨架分割后的二值扫描图像；
图3是图像自相关函数曲线；
图4是模拟模型平面图像；
图5是孔隙半径分布频率曲线；
图6是局部区域孔隙连通性修正前后对比图；
图7是连通性修正后模拟模型平面图像；
图8是去除孤立孔隙后模拟模型平面图像。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合所附图式，作详细说明如下。
如图1所示，图1为本发明的反映真实岩心连通性的2D模型构建方法的流程图。
在步骤101，岩心微观孔隙结构图像预处理。首先选取一张扫描图像作为模型构建对象，构建区域由众多岩石骨架及孔隙空间组成。根据图像灰度值，对孔隙和骨架进行分割处理，处理后图像变为(0,1)的二值图像。如图2，白色区域为岩石骨架，黑色区域为孔隙。流程进入到步骤102。
在步骤102，建立自相关函数。定义指示函数f与模拟区域中任意位置r的关系，用指示函数f(r)来表示。
其中，f(r)为指示函数，无量纲；r表示二值图像的空间位置。
由统计平均得到，孔隙度可以用指示函数表示构建两点自相关函数该函数可以用指示函数表示根据研究需要初步确定微观模型构建原点及范围。
研究区域流程进入到步骤103。
在步骤103，构建图像自相关函数与距离曲线。在图像模拟区域内选取任意一点作为原点，由步骤102得到，原点处自相关函数值统计为φ，无限远处，自相关函数值为φ2。做自相关函数与距离的函数曲线，如图3所示。流程进入到步骤104。
在步骤104，确定平面模型构建最小尺度标准。从平面上应用自相关函数确定具有岩心孔隙结构特征代表性的模型尺度标准，从而确定二维模型的平面范围。在一实施例中，由自相关曲线图3可以得到，原点处，自相关函数为φ＝0.31；距原点距离较远时，自相关函数值在φ2＝0.096上下波动；距原点约20个像素，是自相关函数变化的拐点，大于20个像素的距离时，自相关函数变化较平缓，说明20个像素内，图像的非均质程度已经得到最大程度地体现。根据此尺度标准，选取微观模型大于20个像素的图像区域，如图4所示。对该图像进行连通性分析可发现，无论从哪个方向注入流体，均没有一条连通性通路使得流体从对面出口端流出，模型不具有渗流能力。流程进入到步骤105。
在步骤105，确定纵向孔隙最大贯穿深度。应用孔隙半径分布曲线，选取累积分布频率为95％对应的孔隙大小，该孔隙直径作为纵向孔隙的最大贯穿深度。该深度除以纵向扫描间隔即为可参考的纵向映射图像数。在一实施例中，根据图4对孔隙空间进行特征分析，得到平面孔隙半径分布曲线，如图5所示，对应累积分布频率95％的孔隙半径为50μm。该孔隙直径作为纵向孔隙的最大贯穿深度，该深度除以纵向扫描间隔即为可参考的纵向映射图像数。由于本岩心扫描图像的分辨率为8.5μm，纵向扫描间隔与平面分辨率一致，因此，为反映孔隙的纵向贯穿性，需要选择50×2÷8.5≈12张扫描图像进行映射，修正平面孔隙的连通性。流程进入到步骤106。
在步骤106，纵向孔隙结构映射。在纵向上，将多层扫描图像纵向叠加，层层映射，则可得到反映底层孔隙连通性的用二维图像。在一实施例中，如图6所示，分别为模型局部区域的顶面第1张及第10张孔隙分布图，可以看出，第1张圈内岩石颗粒有胶结，而第10张胶结消失，实际孔隙是连通的。建立二维模型时，则需要将此类孔隙的连通性同时映射到顶面第一张图像，对第1张图像进行连通性修正。将多层扫描图像(层数由步骤105获得)纵向叠加，层层映射，则可得到反映底层孔隙连通性的用二维图像。模型平面图像连通性修正后如图7所示，修正后的模型具备了连通通路，有一定的渗流能力。流程进入到步骤107。
在步骤107，去除孤立孔隙。去除对流体流动不起作用的孤立孔隙，得到孔隙尺度数模模型的平面图像，根据该图像即可建立反映真实岩心孔隙结构特征及连通性的数模模型。如图8所示，依据该二值图像即可建立孔隙模型，模型具有原图像的孔隙结构特征，同时具有更好的连通性。流程结束。

Claims (7)

1.反映真实岩心连通性的2D模型构建方法，其特征在于，该反映真实岩心连通性的2D模型构建方法包括：
步骤1，进行岩心微观孔隙结构图像预处理；
步骤2，建立自相关函数；
步骤3，根据预处理后的岩心微观孔隙结构图像和自相关函数，构建图像自相关函数与距离曲线；
步骤4，从平面上应用图像自相关函数确定具有岩心孔隙结构特征代表性的模型尺度标准，从而确定二维模型的平面范围；
步骤5，确定纵向孔隙最大贯穿深度；
步骤6，进行纵向孔隙结构映射；以及
步骤7，去除孤立孔隙，建立反映真实岩心孔隙结构特征及连通性的模拟模型；
在步骤5中，应用孔隙半径分布曲线，选取累积分布概率为95％对应的孔隙半径大小，孔隙直径作为纵向孔隙的最大贯穿深度，该深度除以纵向扫描间隔即为可参考的纵向映射图像数。
2.根据权利要求1所述的反映真实岩心连通性的2D模型构建方法，其特征在于，在步骤1中，选取一张扫描图像作为模型构建对象，构建区域由众多岩石骨架及孔隙空间组成；根据图像灰度值，对孔隙和骨架进行分割处理，处理后图像变为(0,1)的二值图像。
3.根据权利要求1所述的反映真实岩心连通性的2D模型构建方法，其特征在于，在步骤2中，首先定义指示函数f与模拟区域中任意位置r的关系，用指示函数f(r)来表示：
其中，f(r)为指示函数，无量纲；r表示二值图像的空间位置；
由统计平均得到，孔隙度用指示函数表示构建两点自相关函数该函数用指示函数表示其中h为自相关函数A两点之间的距离；根据研究需要初步确定微观模型构建原点及范围。
4.根据权利要求1所述的反映真实岩心连通性的2D模型构建方法，其特征在于，在步骤3中，在图像模拟区域内选取任意一点作为原点，由步骤2得到，原点处自相关函数值统计为φ，无限远处，自相关函数值为φ2，做自相关函数与距离的函数曲线，其中h为自相关函数A两点之间的距离。
5.根据权利要求1所述的反映真实岩心连通性的2D模型构建方法，其特征在于，在步骤4中，由自相关函数与距离曲线可以得到，原点处，自相关函数为φ；距原点距离较远时，自相关函数值为φ2；其中h为自相关函数A两点之间的距离；距原点N个像素，是自相关函数变化的拐点，大于N个像素的距离时，自相关函数变化较平缓，说明N个像素内，图像的非均质程度已经得到最大程度地体现，根据此尺度标准，选取微观模型大于N个像素的图像区域。
6.根据权利要求1所述的反映真实岩心连通性的2D模型构建方法，其特征在于，在步骤6中，在纵向上，将多层扫描图像纵向叠加，层层映射，得到反映底层孔隙连通性的用二维图像。
7.根据权利要求1所述的反映真实岩心连通性的2D模型构建方法，其特征在于，在步骤7中，去除对流体流动不起作用的孤立孔隙，得到孔隙尺度数模模型的平面图像，根据该图像建立反映真实岩心孔隙结构特征及连通性的数模模型。
CN201410643887.6A 2014-11-10 2014-11-10 反映真实岩心连通性的2d模型构建方法 Active CN105654525B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410643887.6A CN105654525B (zh) 2014-11-10 2014-11-10 反映真实岩心连通性的2d模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410643887.6A CN105654525B (zh) 2014-11-10 2014-11-10 反映真实岩心连通性的2d模型构建方法

Publications (2)

Publication Number Publication Date
CN105654525A CN105654525A (zh) 2016-06-08
CN105654525B true CN105654525B (zh) 2018-06-29

Family

ID=56479571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410643887.6A Active CN105654525B (zh) 2014-11-10 2014-11-10 反映真实岩心连通性的2d模型构建方法

Country Status (1)

Country Link
CN (1) CN105654525B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106227943B (zh) * 2016-07-26 2018-10-23 中国石油大学(华东) 一种二维不同表面粗糙度的随机孔隙生成方法
CN108038903B (zh) * 2017-12-07 2022-01-11 中国石油化工股份有限公司 用于构建岩心模型的三维数字模型生成方法
CN113740224B (zh) * 2021-08-23 2022-12-09 西安交通大学 一种针对多孔介质渗流过程中微观动用孔隙的识别方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101556703A (zh) * 2009-05-16 2009-10-14 中国石油大学(华东) 基于连续切片图像的网络模型建立方法
CN102129548A (zh) * 2011-03-15 2011-07-20 中国石油大学(华东) 多孔介质中流动泡沫结构图像实时采集与识别方法
CN102867302A (zh) * 2012-08-30 2013-01-09 四川大学 基于三维图像信息处理的岩心裂缝识别方法
CN103325117A (zh) * 2013-06-17 2013-09-25 中国石油天然气集团公司 一种基于matlab的岩心图像处理方法及系统
CN103700102A (zh) * 2013-12-20 2014-04-02 电子科技大学 基于ct图像的岩心目标提取方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130325408A1 (en) * 2011-01-28 2013-12-05 Schlumberger Technology Corporation Two dimensional nmr of diffusion and relaxation for material characterization

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101556703A (zh) * 2009-05-16 2009-10-14 中国石油大学(华东) 基于连续切片图像的网络模型建立方法
CN102129548A (zh) * 2011-03-15 2011-07-20 中国石油大学(华东) 多孔介质中流动泡沫结构图像实时采集与识别方法
CN102867302A (zh) * 2012-08-30 2013-01-09 四川大学 基于三维图像信息处理的岩心裂缝识别方法
CN103325117A (zh) * 2013-06-17 2013-09-25 中国石油天然气集团公司 一种基于matlab的岩心图像处理方法及系统
CN103700102A (zh) * 2013-12-20 2014-04-02 电子科技大学 基于ct图像的岩心目标提取方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
图像处理技术在数字岩心建模中的应用;朱益华;《石油天然气学报(江汉石油学院学报)》;20071015;第29卷(第5期);第54-55页第1.1、1.2节，第56页第1.3节 *

Also Published As

Publication number Publication date
CN105654525A (zh) 2016-06-08

Similar Documents

Publication Publication Date Title
CN101556703B (zh) 基于连续切片图像的网络模型建立方法
Lafarge et al. Creating large-scale city models from 3D-point clouds: a robust approach with hybrid representation
CN104850682B (zh) 基于位置的多点地质统计学建模方法
CN105261068B (zh) 基于Micro-CT技术的储层岩心三维实体模型重构方法
US9189862B2 (en) Outline approximation for point cloud of building
CN106803253B (zh) 一种三维岩石图像裂缝识别方法
CN102254354B (zh) 一种实现土的微观结构表面形态三维可视化的方法
CN104573705A (zh) 一种建筑物激光扫描点云数据的聚类方法
CN106126816B (zh) 重复建筑自动感知下的大规模als建筑点云建模方法
CN105279790A (zh) 裂缝网络三维数字岩心建模方法
CN109164026A (zh) 岩石渗流能力评价方法及装置
CN106326517B (zh) 层状裂缝-基质杂交网格建模方法和装置
CN104392486A (zh) 一种点云场景重建方法
CN106324002A (zh) 一种基于岩石分类和多尺度数字岩心的碳酸盐岩孔隙结构表征方法
CN104331699A (zh) 一种三维点云平面化快速搜索比对的方法
CN106157309A (zh) 一种基于虚拟种子点的机载LiDAR地面点云滤波方法
CN106875481B (zh) 一种三维可视化遥感影像地表分类模型的制作方法
CN105654525B (zh) 反映真实岩心连通性的2d模型构建方法
CN104036538A (zh) 土石混合体三维细观结构重建与分析方法及系统
CN107861129A (zh) 一种山地地貌遥感提取方法
CN111754623A (zh) 一种基于多点地质统计学的离散裂缝建模方法
CN105447452A (zh) 一种基于地物空间分布特征的遥感亚像元制图方法
CN115668287A (zh) 源自二维图像的材料性质
CN117115063A (zh) 一种多源数据融合应用方法
CN105653881A (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
GR01 Patent grant
GR01 Patent grant