CN112796743A - 芯上油藏结构生成方法、系统、计算机设备、终端及应用 - Google Patents
芯上油藏结构生成方法、系统、计算机设备、终端及应用 Download PDFInfo
- Publication number
- CN112796743A CN112796743A CN202110011890.6A CN202110011890A CN112796743A CN 112796743 A CN112796743 A CN 112796743A CN 202110011890 A CN202110011890 A CN 202110011890A CN 112796743 A CN112796743 A CN 112796743A
- Authority
- CN
- China
- Prior art keywords
- throat
- pore
- core
- pores
- connection information
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000009825 accumulation Methods 0.000 title claims abstract description 12
- 238000002910 structure generation Methods 0.000 title claims description 14
- 239000011148 porous material Substances 0.000 claims abstract description 131
- 238000012217 deletion Methods 0.000 claims abstract description 17
- 230000037430 deletion Effects 0.000 claims abstract description 17
- 239000011435 rock Substances 0.000 claims abstract description 16
- 238000012163 sequencing technique Methods 0.000 claims abstract description 5
- 238000002474 experimental method Methods 0.000 claims description 13
- 238000004088 simulation Methods 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 2
- 239000003921 oil Substances 0.000 description 48
- 238000010586 diagram Methods 0.000 description 8
- 230000035699 permeability Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000003429 pore cell Anatomy 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Geometry (AREA)
- Geophysics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于油气田开发技术领域，公开了一种芯上油藏结构生成方法、系统、计算机设备、终端及应用，根据真实岩心结构确定岩石孔隙个数、孔喉比、平均配位数，最大配位数；芯片左右两侧均匀设置一列孔隙；以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的圆形孔隙；选取孔隙圆心为剖分节点，并根据剖分结果提取喉道的连接信息；根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序；遍历每个喉道；定义喉道删除概率；对未删除的连接信息进行喉道填充，做矩形，将矩形内格点填充为孔隙格点；连接完毕后将入口出口打开，完成芯上油藏结构构造。本发明生成的多孔介质具有孔喉随机分布、结构特征符合真实岩心的优点。
Description
技术领域
本发明属于油气田开发技术领域，尤其涉及一种芯上油藏结构生成方法、系统、计算机设备、终端及应用。
背景技术
目前：在油水两相流动的微观渗流实验或模拟过程中，需要设计用于实验的多孔介质芯片或者微观渗流模型即芯上油藏。传统的微观多孔介质模型一种为规则排列的孔隙，孔喉大小分布均匀，所有配位数为固定值，与真实岩心孔喉结构特征差异巨大；另一种随机多孔介质模型则无法获得精确的孔喉结构，对孔喉也无法进行定量表征，同样不符合真实岩心孔喉分布特征。以上两种方案均导致多孔介质拓扑结构信息丢失，使得通过微观渗流模拟或实验得到的渗流现象与渗流参数不准确。
多孔介质的孔喉结构分布特征直接决定流体在其中流动过程时的渗流特性，因此设计微观渗流模型或芯片时应保证孔喉分布与真实岩心孔喉分布。用于评价孔喉结构的参数主要为孔隙大小分布、喉道大小分布、孔喉比及配位数等，为使微观渗流实验及模拟结果更加合理，精确明晰油水两相微观渗流机制，急需提出一种可以定量制定芯上油藏结构的方法。
通过上述分析，现有技术存在的问题及缺陷为：
(1)传统的微观多孔介质模型一种为规则排列的孔隙，孔喉大小分布均匀，所有配位数为固定值，与真实岩心孔喉结构特征差异巨大。
(2)传统的随机多孔介质模型则无法获得精确的孔喉结构，对孔喉也无法进行定量表征，同样不符合真实岩心孔喉分布特征。
解决以上问题及缺陷的难度为：建立拓扑结构与真实岩心一致的多孔介质结构，使得孔喉大小、配位数、孔喉比等结构参数与真实岩心结构严格等价。
解决以上问题及缺陷的意义为：突破了芯上油藏结构无法按照真实岩心结构定制的瓶颈，解决了基于芯上油藏的微观渗流实验及模拟结果渗流特性不准的问题。
发明内容
针对现有技术存在的问题，本发明提供了一种芯上油藏结构生成方法、系统、计算机设备、终端及应用。
本发明是这样实现的，一种芯上油藏结构生成方法，所述芯上油藏结构生成方法包括：
确定芯上油藏区域大小，根据真实岩心结构确定岩石孔隙个数、孔喉比、平均配位数，最大配位数；
芯片左右两侧均匀设置一列孔隙，作为芯片的入口和出口；
以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的圆形孔隙；
选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息；
根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序；
遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系，若线段位于圆内或与圆相交则删除该喉道的连接信息；
定义喉道删除概率，遍历每个喉道，生成随机数；
对未删除的连接信息进行喉道填充，做矩形，将矩形内格点填充为孔隙格点；
连接完毕后将入口出口打开，完成芯上油藏结构构造。
进一步，所述芯上油藏结构生成方法确定芯上油藏区域大小为2500×1000，根据真实岩心结构确定岩石孔隙个数为200、孔喉比α＝3.0、平均配位数Cavg＝3.5，最大配位数Cmax＝8；孔隙半径rp分布符合截断威布尔分布：
其中rpmax＝50、rpmin＝10，δ＝0.3、γ＝3.2为任意数，x为0-1的随机数；
喉道半径分布rt由下式确定：
rt＝(1/α)min(rp1,rp2)；
其中α为孔喉比，rp1、rp2为与喉道相连的两个孔隙的半径；
于芯片左右两侧均匀设置一列半径为(50+10)/2＝30的孔隙，每列孔隙个数为
作为芯片的入口和出口。
进一步，所述芯上油藏结构生成方法以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的200-30＝170个圆形孔隙；
选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息。
进一步，所述芯上油藏结构生成方法根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序，保留长度位于前95％的连接信息，删除长度位于后5％即过长喉道的连接信息。
进一步，所述芯上油藏结构生成方法遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系，若线段位于圆内或与圆相交则删除该喉道的连接信息；
定义喉道删除概率p＝(Cmax-Cavg)/Cmax＝0.5625，遍历每个喉道，每次生成0-1的随机数，如果与喉道相连的两个孔隙的配位数都大于1且随机数小于0.5625，则删除该喉道连接信息。
进一步，所述芯上油藏结构生成方法对未删除的连接信息进行喉道填充，以连接节点分别为(10,26)、(45,45)为例，连接两点并在法向方向做直线，以喉道半径大小6为界限截断，分别得到(x1,y1)、(x2,y2)、(x3,y3)及(x4,y4)四个顶点，以四个顶点做矩形。
本发明的另一目的在于提供一种计算机设备，所述计算机设备包括存储器和处理器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，使得所述处理器执行如下步骤：
确定芯上油藏区域大小，根据真实岩心结构确定岩石孔隙个数、孔喉比、平均配位数，最大配位数；
芯片左右两侧均匀设置一列孔隙，作为芯片的入口和出口；
以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的圆形孔隙；
选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息；
根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序；
遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系，若线段位于圆内或与圆相交则删除该喉道的连接信息；
定义喉道删除概率，遍历每个喉道，生成随机数；
对未删除的连接信息进行喉道填充，做矩形，将矩形内格点填充为孔隙格点；
连接完毕后将入口出口打开，完成芯上油藏结构构造。
本发明的另一目的在于提供一种信息数据处理终端，所述信息数据处理终端用于实现所述的芯上油藏结构生成方法。
本发明的另一目的在于提供一种实施所述芯上油藏结构生成方法的芯上油藏结构生成系统，所述芯上油藏结构生成系统包括：
参数确定模块，用于确定芯上油藏区域大小，根据真实岩心结构确定岩石孔隙个数、孔喉比、平均配位数，最大配位数；芯片左右两侧均匀设置一列孔隙，作为芯片的入口和出口；
圆形孔隙生成模块，用于以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的圆形孔隙；
喉道连接信息提取模块，用于选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息；
喉道排序模块，用于根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序；
位置关系确定模块，用于遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系；
喉道删除概率定义模块，用于定义喉道删除概率，遍历每个喉道，每次生成0-1的随机数；
喉道填充模块，用于对未删除的连接信息进行喉道填充，做矩形，将矩形内格点填充为孔隙格点；
油藏结构构造模块，用于连接完毕后将入口出口打开，完成芯上油藏结构构造。
本发明的另一目的在于提供一种所述芯上油藏结构生成方法在微观渗流实验及模拟中微流控芯片及模型的定制与设计中的用途。
结合上述的所有技术方案，本发明所具备的优点及积极效果为：本发明为保证微观渗流实验多孔介质模型即芯上油藏结构符合真实岩心分布，使得实验过程中渗流参数与流体分布更符合油藏地下情况。本发明生成的多孔介质具有孔喉随机分布、结构特征符合真实岩心的优点。本发明涉及微观渗流实验及模拟中的一种基于三角剖分的芯上油藏结构生成方法，适用于微观渗流实验及模拟中所需的微流控芯片及模型的定制与设计。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对本申请实施例中所需要使用的附图做简单的介绍，显而易见地，下面所描述的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的芯上油藏结构生成方法流程图。
图2是本发明实施例提供的芯上油藏结构生成系统的结构示意图；
图2中：1、参数确定模块；2、圆形孔隙生成模块；3、喉道连接信息提取模块；4、喉道排序模块；5、位置关系确定模块；6、喉道删除概率定义模块；7、喉道填充模块；8、油藏结构构造模块。
图3是本发明实施例提供的芯上油藏结构生成方法的实现流程图。
图4是本发明实施例提供的孔隙随机分布图，白色为孔隙，黑色为岩石。
图5是本发明实施例提供的孔隙三角剖分结果图，灰色表示孔隙，线段表示连接信息示意图。
图6是本发明实施例提供的连接删除示意图；(a)位置表示过长连接，(b)位置表示随机删除连接，(c)位置表示与孔隙重叠的连接。
图7是本发明实施例提供的喉道填充连接示意图；(a)为矩形区域确定示意图，(b)为填充示意图。
图8是本发明实施例提供的芯上油藏定制结构结果图；白色为孔隙，黑色为岩石.
图9是本发明实施例提供的基于芯上油藏定制结构的数值模拟结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
针对现有技术存在的问题，本发明提供了一种芯上油藏结构生成方法、系统、计算机设备、终端及应用，下面结合附图对本发明作详细的描述。
如图1所示，本发明提供的芯上油藏结构生成方法包括以下步骤：
S101：确定芯上油藏区域大小，根据真实岩心结构确定岩石孔隙个数、孔喉比、平均配位数，最大配位数；
S102：芯片左右两侧均匀设置一列孔隙，作为芯片的入口和出口；
S103：以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的圆形孔隙；
S104：选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息；
S105：根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序，保留长度位于前95％的连接信息，删除长度位于后5％即过长喉道的连接信息；
S106：遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系，若线段位于圆内或与圆相交则删除该喉道的连接信息；
S107：定义喉道删除概率，遍历每个喉道，每次生成0-1的随机数，如果与喉道相连的两个孔隙的配位数都大于1且随机数小于0.5625，则删除该喉道连接信息；
S108：对未删除的连接信息进行喉道填充，做矩形，将矩形内格点填充为孔隙格点；
S109：连接完毕后将入口出口打开，完成芯上油藏结构构造。
本发明提供的芯上油藏结构生成方法业内的普通技术人员还可以采用其他的步骤实施，图1的本发明提供的芯上油藏结构生成方法仅仅是一个具体实施例而已。
如图2所示，本发明提供的芯上油藏结构生成系统包括：
参数确定模块1，用于确定芯上油藏区域大小，根据真实岩心结构确定岩石孔隙个数、孔喉比、平均配位数，最大配位数；芯片左右两侧均匀设置一列孔隙，作为芯片的入口和出口；
圆形孔隙生成模块2，用于以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的圆形孔隙；
喉道连接信息提取模块3，用于选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息；
喉道排序模块4，用于根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序；
位置关系确定模块5，用于遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系；
喉道删除概率定义模块6，用于定义喉道删除概率，遍历每个喉道，每次生成0-1的随机数；
喉道填充模块7，用于对未删除的连接信息进行喉道填充，做矩形，将矩形内格点填充为孔隙格点；
油藏结构构造模块8，用于连接完毕后将入口出口打开，完成芯上油藏结构构造。
下面结合附图对本发明的技术方案作进一步的描述。
本发明实施例提供的芯上油藏结构生成方法，如图3所示，包括以下步骤：
(1)确定芯上油藏区域大小为2500×1000，根据真实岩心结构确定岩石孔隙个数为200、孔喉比α＝3.0、平均配位数Cavg＝3.5，最大配位数Cmax＝8；孔隙半径rp分布符合截断威布尔分布：
其中rpmax＝50、rpmin＝10，δ＝0.3、γ＝3.2为任意数，x为0-1的随机数；
喉道半径分布rt由下式确定：
rt＝(1/α)min(rp1,rp2)；
其中α为孔喉比，rp1、rp2为与喉道相连的两个孔隙的半径；
(2)于芯片左右两侧均匀设置一列半径为(50+10)/2＝30的孔隙，每列孔隙个数为
作为芯片的入口和出口。
(3)以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的200-30＝170个圆形孔隙，结果如图4所示。
(4)选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分(Lee D.T.,Schachter B.J.,International Journal of Computer&Information Sciences,1980,9(3),219-242)，如图5所示，并根据剖分结果提取喉道的连接信息即图3中的线段。
(5)根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序，保留长度位于前95％的连接信息，删除长度位于后5％即过长喉道的连接信息，过长的连接信息如图6中(a)位置所示。
(6)遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系，若线段位于圆内或与圆相交则删除该喉道的连接信息，相交连接信息如图6中(c)位置所示。
(7)定义喉道删除概率p＝(Cmax-Cavg)/Cmax＝0.5625，遍历每个喉道，每次生成0-1的随机数，如果与喉道相连的两个孔隙的配位数都大于1且随机数小于0.5625，则删除该喉道连接信息，随机删除的过多连接信息如图6中(b)位置所示。
(8)对未删除的连接信息进行喉道填充，以连接节点分别为(10,26)、(45,45)为例，连接两点并在法向方向做直线，以喉道半径大小6为界限截断，分别得到(x1,y1)、(x2,y2)、(x3,y3)及(x4,y4)四个顶点，以四个顶点做矩形，如图7(a)所示，然后将矩形内格点填充为孔隙格点，填充示意图如图7(b)所示。
(9)连接完毕后将入口出口打开，完成芯上油藏结构构造，如图8所示。
(10)采用该方法生成的结构定制微流控芯片进行微观渗流实验或直接进行微观渗流数值模拟，采用格子玻尔兹曼方法的水驱油微观数值模拟结果如图9所示，其可清晰地捕捉每个孔喉中的局部流动事件，如贾敏效应、盲端剩余油等，为微观渗流机制及提高原油采收率研究提供了有力手段。
本发明基于微观渗流模拟得到芯上油藏的渗透率(Cihan A.,Sukop M.C.,TynerJ.S.,et al.Analytical predictions and lattice Boltzmann simulations ofintrinsic permeability for mass fractal porous media.Vadose Zone Journal,2009,8(1):187-196)，该方案方法构造的多孔介质结构渗透率模拟计算结果与真实岩心室内驱替实验结果进行对比误差仅为5.2％，而传统方法构造的规则多孔介质及随机多孔介质渗透率误差分别高达25.1％及30.8％。
应当注意，本发明的实施方式可以通过硬件、软件或者软件和硬件的结合来实现。硬件部分可以利用专用逻辑来实现；软件部分可以存储在存储器中，由适当的指令执行系统，例如微处理器或者专用设计硬件来执行。本领域的普通技术人员可以理解上述的设备和方法可以使用计算机可执行指令和/或包含在处理器控制代码中来实现，例如在诸如磁盘、CD或DVD-ROM的载体介质、诸如只读存储器(固件)的可编程的存储器或者诸如光学或电子信号载体的数据载体上提供了这样的代码。本发明的设备及其模块可以由诸如超大规模集成电路或门阵列、诸如逻辑芯片、晶体管等的半导体、或者诸如现场可编程门阵列、可编程逻辑设备等的可编程硬件设备的硬件电路实现，也可以用由各种类型的处理器执行的软件实现，也可以由上述硬件电路和软件的结合例如固件来实现。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种芯上油藏结构生成方法，其特征在于，所述芯上油藏结构生成方法包括：
确定芯上油藏区域大小，根据真实岩心结构确定岩石孔隙个数、孔喉比、平均配位数，最大配位数；
芯片左右两侧均匀设置一列孔隙，作为芯片的入口和出口；
以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的圆形孔隙；
选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息；
根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序；
遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系，若线段位于圆内或与圆相交则删除该喉道的连接信息；
定义喉道删除概率，遍历每个喉道，生成随机数；
对未删除的连接信息进行喉道填充，做矩形，将矩形内格点填充为孔隙格点；
连接完毕后将入口出口打开，完成芯上油藏结构构造。
2.如权利要求1所述的芯上油藏结构生成方法，其特征在于，所述芯上油藏结构生成方法确定芯上油藏区域大小为2500×1000，根据真实岩心结构确定岩石孔隙个数为200、孔喉比α＝3.0、平均配位数Cavg＝3.5，最大配位数Cmax＝8；孔隙半径rp分布符合截断威布尔分布：
其中rpmax＝50、rpmin＝10，δ＝0.3、γ＝3.2为任意数，x为0-1的随机数；
喉道半径分布rt由下式确定：
rt＝(1/α)min(rp1,rp2)；
其中α为孔喉比，rp1、rp2为与喉道相连的两个孔隙的半径；
于芯片左右两侧均匀设置一列半径为(50+10)/2＝30的孔隙，每列孔隙个数为
作为芯片的入口和出口。
3.如权利要求1所述的芯上油藏结构生成方法，其特征在于，所述芯上油藏结构生成方法以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的200-30＝170个圆形孔隙；
选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息。
4.如权利要求1所述的芯上油藏结构生成方法，其特征在于，所述芯上油藏结构生成方法根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序，保留长度位于前95％的连接信息，删除长度位于后5％即过长喉道的连接信息。
5.如权利要求1所述的芯上油藏结构生成方法，其特征在于，所述芯上油藏结构生成方法遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系，若线段位于圆内或与圆相交则删除该喉道的连接信息；
定义喉道删除概率p＝(Cmax-Cavg)/Cmax＝0.5625，遍历每个喉道，每次生成0-1的随机数，如果与喉道相连的两个孔隙的配位数都大于1且随机数小于0.5625，则删除该喉道连接信息。
6.如权利要求1所述的芯上油藏结构生成方法，其特征在于，所述芯上油藏结构生成方法对未删除的连接信息进行喉道填充，以连接节点分别为(10,26)、(45,45)为例，连接两点并在法向方向做直线，以喉道半径大小6为界限截断，分别得到(x1,y1)、(x2,y2)、(x3,y3)及(x4,y4)四个顶点，以四个顶点做矩形。
7.一种计算机设备，其特征在于，所述计算机设备包括存储器和处理器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，使得所述处理器执行如下步骤：
确定芯上油藏区域大小，根据真实岩心结构确定岩石孔隙个数、孔喉比、平均配位数，最大配位数；
芯片左右两侧均匀设置一列孔隙，作为芯片的入口和出口；
以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的圆形孔隙；
选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息；
根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序；
遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系，若线段位于圆内或与圆相交则删除该喉道的连接信息；
定义喉道删除概率，遍历每个喉道，生成随机数；
对未删除的连接信息进行喉道填充，做矩形，将矩形内格点填充为孔隙格点；
连接完毕后将入口出口打开，完成芯上油藏结构构造。
8.一种信息数据处理终端，其特征在于，所述信息数据处理终端用于实现权利要求1～6任意一项所述的芯上油藏结构生成方法。
9.一种实施权利要求1～6任意一项所述芯上油藏结构生成方法的芯上油藏结构生成系统，其特征在于，所述芯上油藏结构生成系统包括：
参数确定模块，用于确定芯上油藏区域大小，根据真实岩心结构确定岩石孔隙个数、孔喉比、平均配位数，最大配位数；芯片左右两侧均匀设置一列孔隙，作为芯片的入口和出口；
圆形孔隙生成模块，用于以孔隙大小为半径，按孔隙分布规律随机在芯片区域中生成互不重叠的圆形孔隙；
喉道连接信息提取模块，用于选取孔隙圆心为剖分节点，基于所有节点进行Delaunay三角剖分，并根据剖分结果提取喉道的连接信息；
喉道排序模块，用于根据连接两端孔隙的坐标计算喉道长度，按照喉道长度自短至长排序；
位置关系确定模块，用于遍历每个喉道，查找与该喉道两端节点连接的所有孔隙节点集合U，依次判断该喉道构成的线段与集合U中孔隙构成的圆的位置关系；
喉道删除概率定义模块，用于定义喉道删除概率，遍历每个喉道，每次生成0-1的随机数；
喉道填充模块，用于对未删除的连接信息进行喉道填充，做矩形，将矩形内格点填充为孔隙格点；
油藏结构构造模块，用于连接完毕后将入口出口打开，完成芯上油藏结构构造。
10.一种如权利要求1～6任意一项所述芯上油藏结构生成方法在微观渗流实验及模拟中微流控芯片及模型的定制与设计中的用途。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110011890.6A CN112796743B (zh) | 2021-01-06 | 2021-01-06 | 芯上油藏结构生成方法、系统、计算机设备、终端及应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110011890.6A CN112796743B (zh) | 2021-01-06 | 2021-01-06 | 芯上油藏结构生成方法、系统、计算机设备、终端及应用 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112796743A true CN112796743A (zh) | 2021-05-14 |
| CN112796743B CN112796743B (zh) | 2022-06-28 |
Family
ID=75808478
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110011890.6A Active CN112796743B (zh) | 2021-01-06 | 2021-01-06 | 芯上油藏结构生成方法、系统、计算机设备、终端及应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112796743B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113268791A (zh) * | 2021-06-08 | 2021-08-17 | 清华大学 | 一种桌面油藏及其制备方法 |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB191506718A (en) * | 1915-05-05 | 1916-05-05 | Eric Leopold Hugo Cosby | Improvements in or relating to Well Drilling Tools for Oil, Gas and Water Wells. |
| CN104331579A (zh) * | 2014-11-19 | 2015-02-04 | 中国石油大学(华东) | 一种低渗透储层原油边界层的模拟方法 |
| US20160010440A1 (en) * | 2013-08-02 | 2016-01-14 | Olga Nikolaevna Sizonenko | Oil production intensification device and method |
| CN106021788A (zh) * | 2016-05-31 | 2016-10-12 | 中国石油大学(华东) | 一种基于孔喉结构特征的致密储层分级评价标准划分方法 |
| CN107939387A (zh) * | 2017-11-14 | 2018-04-20 | 中国石油大学(北京) | 一种微观岩石网络模型制作方法 |
| CN108150162A (zh) * | 2018-01-16 | 2018-06-12 | 中国石油大学(北京) | 一种微观油藏模型及其制作方法和使用方法 |
| CN108733902A (zh) * | 2018-05-03 | 2018-11-02 | 中国石油天然气股份有限公司 | 油藏渗透率的确定方法、装置及存储介质 |
| CN109087301A (zh) * | 2018-06-29 | 2018-12-25 | 中国石油大学(华东) | 一种基于居中轴线和表面模型的岩心喉道分割方法 |
| CN109444201A (zh) * | 2018-05-25 | 2019-03-08 | 中国石油大学(华东) | 一种测定致密岩心多孔介质中多相流体流动特征的核磁共振实验装置及方法 |
| CN110302853A (zh) * | 2019-07-03 | 2019-10-08 | 清华大学 | 一种芯片油藏及其制作方法 |
| CN110362842A (zh) * | 2018-04-09 | 2019-10-22 | 长江大学 | 基于多种形状孔喉的随机孔隙网络模型建模方法 |
| CN110702580A (zh) * | 2019-10-21 | 2020-01-17 | 西安石油大学 | 一种基于信息熵的致密砂岩储层孔喉非均质性表征方法 |
| CN111151316A (zh) * | 2020-01-16 | 2020-05-15 | 西安石油大学 | 一种可视化研究微观渗吸与孔喉比关系的微流控芯片、实验装置及方法 |
| CN111151315A (zh) * | 2020-01-16 | 2020-05-15 | 西安石油大学 | 一种研究微米尺度孔喉配位数与渗吸效率的微流控芯片、实验装置及方法 |
| CN112052634A (zh) * | 2020-09-07 | 2020-12-08 | 中国石油大学(华东) | 低渗油藏注水过程中颗粒运移影响的分析方法及装置 |
| CN112098293A (zh) * | 2020-08-03 | 2020-12-18 | 西南石油大学 | 基于孔隙裂缝双重介质气藏非稳态气水两相渗流模拟方法 |
-
2021
- 2021-01-06 CN CN202110011890.6A patent/CN112796743B/zh active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB191506718A (en) * | 1915-05-05 | 1916-05-05 | Eric Leopold Hugo Cosby | Improvements in or relating to Well Drilling Tools for Oil, Gas and Water Wells. |
| US20160010440A1 (en) * | 2013-08-02 | 2016-01-14 | Olga Nikolaevna Sizonenko | Oil production intensification device and method |
| CN104331579A (zh) * | 2014-11-19 | 2015-02-04 | 中国石油大学(华东) | 一种低渗透储层原油边界层的模拟方法 |
| CN106021788A (zh) * | 2016-05-31 | 2016-10-12 | 中国石油大学(华东) | 一种基于孔喉结构特征的致密储层分级评价标准划分方法 |
| CN107939387A (zh) * | 2017-11-14 | 2018-04-20 | 中国石油大学(北京) | 一种微观岩石网络模型制作方法 |
| CN108150162A (zh) * | 2018-01-16 | 2018-06-12 | 中国石油大学(北京) | 一种微观油藏模型及其制作方法和使用方法 |
| CN110362842A (zh) * | 2018-04-09 | 2019-10-22 | 长江大学 | 基于多种形状孔喉的随机孔隙网络模型建模方法 |
| CN108733902A (zh) * | 2018-05-03 | 2018-11-02 | 中国石油天然气股份有限公司 | 油藏渗透率的确定方法、装置及存储介质 |
| CN109444201A (zh) * | 2018-05-25 | 2019-03-08 | 中国石油大学(华东) | 一种测定致密岩心多孔介质中多相流体流动特征的核磁共振实验装置及方法 |
| CN109087301A (zh) * | 2018-06-29 | 2018-12-25 | 中国石油大学(华东) | 一种基于居中轴线和表面模型的岩心喉道分割方法 |
| CN110302853A (zh) * | 2019-07-03 | 2019-10-08 | 清华大学 | 一种芯片油藏及其制作方法 |
| CN110702580A (zh) * | 2019-10-21 | 2020-01-17 | 西安石油大学 | 一种基于信息熵的致密砂岩储层孔喉非均质性表征方法 |
| CN111151316A (zh) * | 2020-01-16 | 2020-05-15 | 西安石油大学 | 一种可视化研究微观渗吸与孔喉比关系的微流控芯片、实验装置及方法 |
| CN111151315A (zh) * | 2020-01-16 | 2020-05-15 | 西安石油大学 | 一种研究微米尺度孔喉配位数与渗吸效率的微流控芯片、实验装置及方法 |
| CN112098293A (zh) * | 2020-08-03 | 2020-12-18 | 西南石油大学 | 基于孔隙裂缝双重介质气藏非稳态气水两相渗流模拟方法 |
| CN112052634A (zh) * | 2020-09-07 | 2020-12-08 | 中国石油大学(华东) | 低渗油藏注水过程中颗粒运移影响的分析方法及装置 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113268791A (zh) * | 2021-06-08 | 2021-08-17 | 清华大学 | 一种桌面油藏及其制备方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112796743B (zh) | 2022-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3710710B2 (ja) | Ｉｃレイアウトにおけるポリゴン表現 | |
| US6816997B2 (en) | System and method for performing design rule check | |
| Lauther | A min-cut placement algorithm for general cell assemblies based on a graph representation | |
| EP0848342A2 (en) | Intergrated circuit design method, database apparatus for designing integrated circuit and integrated circuit design support apparatus | |
| CN110851956B (zh) | 一种建设工程管道施工定位信息自动化计算、标注及出图方法 | |
| CN112796743B (zh) | 芯上油藏结构生成方法、系统、计算机设备、终端及应用 | |
| CN110457843A (zh) | 三维管线模型生成方法、装置、计算机设备和存储介质 | |
| US9817941B2 (en) | Methods, systems, and articles of manufacture for implementing high current carrying interconnects in electronic designs | |
| US11334705B2 (en) | Electrical circuit design using cells with metal lines | |
| CN111597678B (zh) | 一种多洪源遭遇条件洪水淹没特征的解耦方法 | |
| CN112800716B (zh) | 一种集成电路版图布线中线网开路检测方法 | |
| CN103077330A (zh) | 存在垂向水量交换情况下的地下水全局流线可视化方法 | |
| US7174526B2 (en) | Accurate density calculation with density views in layout databases | |
| CN112231877B (zh) | 模拟核电厂第一跨厂房泄洪方案的方法和装置 | |
| CN109711680A (zh) | 基于gis数字化技术的河流水电规划方法 | |
| CN100367284C (zh) | 边界数据的内外判定方法 | |
| US7277837B2 (en) | Method for flow analysis utilizing cell elevation topology | |
| EP3118816A1 (en) | Computational modelling for engineering analysis | |
| López et al. | Postbreakthrough behavior in flow through porous media | |
| CN114882144A (zh) | 河流拓扑结构的生成方法、装置、计算设备及介质 | |
| RU2626974C2 (ru) | Способ определения карстовой области | |
| KR20220161491A (ko) | 가상 제조 환경에서 국부적 CDU (critical dimension uniformity) 모델링 및 제어를 수행하기 위한 시스템 및 방법 | |
| CN113177285A (zh) | 基于排水模型的城市排水系统拓扑关系问题标识与修正法 | |
| EP0644496B1 (en) | Method and system for dividing analyzing region in device simulator | |
| CN109326004B (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 |