CN108734751A - 一种基于随机分形的多孔介质细观结构模型生成方法 - Google Patents
一种基于随机分形的多孔介质细观结构模型生成方法 Download PDFInfo
- Publication number
- CN108734751A CN108734751A CN201810446651.1A CN201810446651A CN108734751A CN 108734751 A CN108734751 A CN 108734751A CN 201810446651 A CN201810446651 A CN 201810446651A CN 108734751 A CN108734751 A CN 108734751A
- Authority
- CN
- China
- Prior art keywords
- pore
- model
- porous media
- random
- boundary
- 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 29
- 239000011148 porous material Substances 0.000 claims abstract description 42
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims abstract description 6
- 230000009466 transformation Effects 0.000 claims abstract description 3
- 239000011159 matrix material Substances 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 15
- 230000003698 anagen phase Effects 0.000 claims description 8
- 230000004069 differentiation Effects 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 230000007773 growth pattern Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 3
- 239000007790 solid phase Substances 0.000 description 7
- 230000037431 insertion Effects 0.000 description 3
- 238000003780 insertion Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000002948 stochastic simulation Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000000151 deposition Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000012010 growth Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/001—Texturing; Colouring; Generation of texture or colour
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明提供一种基于随机分形的多孔介质细观结构模型生成方法,可模拟多孔介质内部细观结构的分布情况。通过编写程序语言模拟孔隙生成,实现孔隙模型的规范化处理,针对多孔介质中孔隙边界的复杂、不确定等形状,通过林顿伊尔系统仿真模式进行模拟,实行字符串序列的等价变换,将字符串序列转化为边界演变规则作为边界生成的依据。同时,基于分形几何的迭代方式,并结合随机分布概率模拟孔隙结构的生成,得出多样化复杂的多孔介质细观结构模型。本发明的仿真效果良好,可有效揭示多孔介质内部细观结构特性,在软件模型建立与三维重构方面可实现广泛应用。
Description
技术领域
本发明属于数字图像与材料工程互相结合应用的领域,涉及一种基于随机分形的多孔介质细观结构模型生成方法,适用于通过编写程序语言模拟多孔介质细观结构的分布并最终生成模型图像。
背景技术
多孔介质内部细观结构的特征对其宏观力学性能有重要的影响,诸如材料的渗透性、抗压性等均与其内部细观结构的分布、排列有着密切的关系。因此,如何更切合实际地模拟多孔介质内部的细观结构并生成图像模型已成为该领域研究的热点问题,通过图像的生成,可以有效揭示多孔介质细、宏观之间的关系。目前在生成多孔介质细观结构图像的方法,主要有离散的减少距离方法生成圆形、椭圆形的两相粒状多孔介质,或采用沉积生成算法生成多孔介质。但这些方法均受限于固定孔隙模型、随机性较差,且算法过于复杂化,成图的理想程度较低,无法实现模型与实体的有效逼真与高吻合度。对此研究了一种利用随机分形迭代模式模拟多孔介质细观结构的生成方法,该方法基于林顿伊尔系统仿真模式模拟孔隙的复杂边界,由分形几何迭代方式生成多样化无规则孔隙结构,并通过分布概率的大小,随机排列于多孔介质固相内部,初始定义多孔介质的目标孔隙率从而决定孔隙的无规则定量分布,最终基于程序汇编语言将数字信息转化为图像并绘制成图。整个模型生成过程逻辑严谨,模型与实体有着较高的匹配度,能有效表征实体多孔介质内部结构特征,在软件建模平台以及多孔介质三维重构的过程中实现广泛的应用。
发明内容
针对现有技术的缺陷,发明一种基于随机分形的多孔介质细观结构模型生成方法,以定义多孔介质目标孔隙率作为控制整个模拟的终结条件,确定绘图区域范围,初始假设区域为多孔介质的固相,令孔隙为生长相,固相为非生长相。由林顿伊尔系统仿真模式模拟孔隙复杂边界,通过分形几何迭代方法随机模拟多样化无规则孔隙,遍历绘图区域建立全局域坐标矩阵,同时引入随机分布概率元素,实现孔隙模型在各坐标的随机分布,建立孔隙初始元大小、孔隙迭代次数与分布概率的关联,最终实现复杂孔隙模型的生成与嵌入,并导出模型图像。该方法生成的孔隙结构随机性强、形状多样化组合,对多孔介质细观结构的模拟更为逼真、有效。
一种基于随机分形的多孔介质细观结构模型生成方法,其发明内容如下:
1.复杂孔隙边界模型的生成规则。
(1)基于林顿伊尔系统仿真模式生成字符串序列,读取并识别字符串序列特征,实现字符串的自动存储;
(2)对字符串进行等价变换,构造模型生成规则,实现孔隙边界复杂模拟。
2.基于分形几何迭代方式生成多样化无规则孔隙。
(1)定义孔隙生长模型,给出几种可能的初始元,对初始元进行编号处理并赋值ID号提供识别功能,为后续绘图实时调用;
(2)引入随机因子,并结合林顿伊尔系统仿真模式随机迭代过程,实现孔隙的演变处理过程;
(3)运用分形几何迭代方式,汇合初始元编号、随机因子形成相互对应关系,实时模拟孔隙演变过程,随机迭代N次后得出复杂孔隙形状;
3.初始元大小与绘制区域坐标对应关系的建立。
(1)遍历绘图区域,给定全局域内随机坐标矩阵;
(2)定义随机分布概率元素与坐标矩阵的对应关系,组装新矩阵,以分布概率元素确定初始元的大小,以坐标元素确定孔隙生成位置,实现孔隙在固相内部的生成与嵌入;
4.基于图形的像素值计算孔隙率的大小,控制整体绘制过程。
(1)整个绘图过程执行黑白二值处理,黑色块为固相,白色块为孔隙,通过计算图像白色区域块像素占全区域的比例,即可得出实时孔隙率的大小;
(2)以实时跟踪的孔隙率大小实现迭代与绘图循环次数的控制,直至达到初始设定的目标孔隙率,完成整个绘图过程,生成多孔介质细观结构模型。
本发明具有以下优点:
(1)可依据林顿伊尔系统的字符串序列控制复杂孔隙边界的生成,随机效果好,符合实际材料结构特征。
(2)基于分形迭代方式构造多样化无规则孔隙模型,可有效模型孔隙的封闭、连通、走向、大小各异等更为全面的性状,仿真效果更好。
(3)建立绘图区域坐标矩阵、分布概率与孔隙模型迭代、嵌入方式的对接关系,由分布概率操控孔隙模型大小,并决定孔隙的分形迭代次数,实现的多重分布,模型结构进一步优化。
(4)建立图像的像素值识别功能,将二值图像数据经过计算反演实时孔隙率大小,通过孔隙率大小智能控制整个绘图流程,以达到目标孔隙率为终止条件,完成模型构建,生成模型图像,有效表征多孔介质内部结构特征。
附图说明
附图1为林顿伊尔系统仿真字符串序列与对应孔隙边界生成模型。
附图2为几种可能的孔隙初始元与对应迭代模型。
附图3为模型结构生成流程图。
附图4为不同孔隙率对应的多孔介质细观结构模型图。
具体实施方式
下面结合附图对本发明做进一步的说明。
通过编写程序语言模拟多孔介质细观结构,主要包括:基于林顿伊尔系统随机模拟孔隙边界;基于分形几何随机迭代过程生成孔隙模型;遍历固相区域并建立全局域坐标矩阵,引入随机分布概率元素与坐标矩阵进行组装,建立组装矩阵与孔隙生成元模型、大小与孔隙迭代次数的对应关系、于非生长相内部生成孔隙结构;识别图像单像素大小,将其反演为孔隙率从而控制整个模型生成流程,最终生成多孔介质细观模型,具体实现步骤如下:
1.选定模型生成区域,确定目标孔隙率控制参数。
(1)给定绘图区域,并对其进行坐标标定,将整片区域初始化,内部以全黑色进行填充,计算图像的总像素R,即假定最初状态全为非生长相(固相);
(2)设定目标孔隙率n,即以n的大小作为目标控制参数,决定着整个模型生成是否的终结条件。
2.基于林顿伊尔系统的仿真模式随机模拟孔隙边界。
(1)建立林顿伊尔系统仿真模式。给定初始字符,定义字符替换规则将初始字符进行逐步替换,形成新字符串序列,其替换基本思路如下:
替换规则:b→a,a→ab
演变规律:b→a→ab→aba→abaab→abaababa→abaababaabaab→…
(2)对新字符串序列实行等价替换,用“F”替换“b”,用“+”或“-”替换“a”,“F”表示孔隙边界生成,“+”或“-”表示边界生成的两个随机方向,“+”为右转60度,“-”为左转60度;
(3)设置孔隙边界生成规则,声明孔隙生成起始点、起始边界长度、边界生长方向以及迭代步长等控制参数,读取字符串序列从而控制孔隙边界生成模式完成边界生成模型过程,如图1所示为几种可能的孔隙边界生成模型;
3.基于分形几何迭代方式随机生成孔隙模型。
(1)给定几种可能的最初孔隙模式,并对初始元赋值ID识别号,并随机进行第一次模型演变,得出孔隙的初始元模型,如图2中所示的初始元;
(2)结合第1步的林顿伊尔系统仿真孔隙边界生成模型,同时加入随机因子,基于分形几何的迭代方式,实时进行边界、孔隙形状的演变,生成孔隙模型,如图2的几种初始元经过迭代后的孔隙模型;
4.建立非生长相内部坐标矩阵与分布概率的关联,实现孔隙大小、形状多样化模拟。
(1)建立非生长相内部随机坐标矩阵A[xi,yi],实现坐标矩阵对绘图区域的完全遍历;
(2)建立分布概率pi,分布概率应满足如下条件:
将各分布概率随机嵌入坐标矩阵A[xi,yi]中,形成新的矩阵B[xi,yi,pi],在矩阵B中,坐标元素xi,yi决定着孔隙生成位置的随机分布,分布概率元素pi决定着孔隙初始元的大小与迭代次数,pi越大,孔隙初始元面积越大,且迭代次数越多,从而实现大、小孔隙与形状多样化孔隙的随机分布,生成孔隙,孔隙模型以白色块嵌入、生长于非生长相内部。
5.计算白色区域图像像素大小r,将白色区域所占比例转换为实时孔隙率大小np,实现模型生成过程循环迭代步数的实时操控。
(1)每完成一步迭代过程,系统进行白色区域的识别与像素值的计算,同时进行实时孔隙率np的计算,即:np=r/R;
(2)循环终结条件的判定,若np<n,则重复2、3、4、5步操作;若np≥n,则终止循环迭代过程,完成孔隙的演变与嵌入,整个绘图流程如图3所示。
6.终止迭代,生成并导出图像,得出多孔介质细观结构模型图。如图4所示为不同目标孔隙率的多孔介质细观结构模型图。
Claims (1)
1.一种基于随机分形的多孔介质细观结构模型生成方法,其特征在于,包括:
(1)基于林顿伊尔系统仿真模式,生成字符串发展序列并自动存储,以字符串序列模式实现孔隙边界模型的随机控制,通过字符串的拓展与演变,系统识别字符串排列特征,并执行字符串的等价变换,建立起孔隙边界模型的构造规则,实现孔隙边界的复杂连续性处理,体现形态各异孔隙边界的组合方式,更好的模拟实际孔隙边界;
(2)定义异形孔隙初始元模型,运用分形几何迭代方式对孔隙初始元进行有限次迭代,模拟孔隙演变过程,结合林顿伊尔系统仿真模式并加入随机因子决定初始元迭代方向,实现孔隙生成过程的多样化、吻合度高与随机性强的无规则模拟以及孔隙模型复杂变化的动态生长方式;
(3)对绘图区进行全局域遍历,建立非生长相内部坐标矩阵,同时给出区域内部随机概率分布元素,形成坐标矩阵、随机概率相统一的组装矩阵,控制孔隙初始元的大小与迭代步数,实现孔隙大小的随机分布,孔隙结构模型进一步优化处理;
(4)以指定色块进行孔隙结构对非生长相的实时嵌入,并跟踪该色块对全区域所占的比例,通过该比例反演孔隙率的实时大小,通过与目标孔隙率的比较,判定孔隙生成条件,自动控制循环迭代过程。以达到循环迭代条件而终止整个模型生成,得出多孔介质细观结构模型二值图像,有效表征实际多孔介质内部结构特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810446651.1A CN108734751B (zh) | 2018-05-08 | 2018-05-08 | 一种基于随机分形的多孔介质细观结构模型生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810446651.1A CN108734751B (zh) | 2018-05-08 | 2018-05-08 | 一种基于随机分形的多孔介质细观结构模型生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108734751A true CN108734751A (zh) | 2018-11-02 |
CN108734751B CN108734751B (zh) | 2022-02-11 |
Family
ID=63938076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810446651.1A Active CN108734751B (zh) | 2018-05-08 | 2018-05-08 | 一种基于随机分形的多孔介质细观结构模型生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108734751B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829975A (zh) * | 2019-03-14 | 2019-05-31 | 河海大学 | 一种变孔隙度多孔介质构建方法 |
CN113298923A (zh) * | 2021-05-31 | 2021-08-24 | 北京大学 | 一种仿真电子画的生成方法、装置、设备及存储介质 |
CN115859402A (zh) * | 2022-12-01 | 2023-03-28 | 哈尔滨师范大学 | 一种基于改进四参数随机生长法的植物纤维地膜结构模拟方法及系统 |
CN116884540A (zh) * | 2023-06-15 | 2023-10-13 | 湖北工业大学 | 多相材料三维模型生成与连通性判断方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6021662A (en) * | 1997-12-15 | 2000-02-08 | Institut Francais Du Petrole | Method for modeling fluid displacements in a porous medium |
US20020173915A1 (en) * | 2000-07-17 | 2002-11-21 | Patrick Egermann | Method for modelling fluid displacements in a porous environment taking into account hysteresis effects |
CN104535475A (zh) * | 2015-01-08 | 2015-04-22 | 中国石油大学(北京) | 碳酸盐岩微观结构的确定方法及装置 |
CN104573198A (zh) * | 2014-12-23 | 2015-04-29 | 长江大学 | 基于随机分形理论的数字岩心及孔隙网络模型重构方法 |
-
2018
- 2018-05-08 CN CN201810446651.1A patent/CN108734751B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6021662A (en) * | 1997-12-15 | 2000-02-08 | Institut Francais Du Petrole | Method for modeling fluid displacements in a porous medium |
US20020173915A1 (en) * | 2000-07-17 | 2002-11-21 | Patrick Egermann | Method for modelling fluid displacements in a porous environment taking into account hysteresis effects |
CN104573198A (zh) * | 2014-12-23 | 2015-04-29 | 长江大学 | 基于随机分形理论的数字岩心及孔隙网络模型重构方法 |
CN104535475A (zh) * | 2015-01-08 | 2015-04-22 | 中国石油大学(北京) | 碳酸盐岩微观结构的确定方法及装置 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829975A (zh) * | 2019-03-14 | 2019-05-31 | 河海大学 | 一种变孔隙度多孔介质构建方法 |
CN113298923A (zh) * | 2021-05-31 | 2021-08-24 | 北京大学 | 一种仿真电子画的生成方法、装置、设备及存储介质 |
CN113298923B (zh) * | 2021-05-31 | 2023-11-07 | 北京大学 | 一种仿真电子画的生成方法、装置、设备及存储介质 |
CN115859402A (zh) * | 2022-12-01 | 2023-03-28 | 哈尔滨师范大学 | 一种基于改进四参数随机生长法的植物纤维地膜结构模拟方法及系统 |
CN115859402B (zh) * | 2022-12-01 | 2023-09-12 | 哈尔滨师范大学 | 一种基于改进四参数随机生长法的植物纤维地膜结构模拟方法及系统 |
CN116884540A (zh) * | 2023-06-15 | 2023-10-13 | 湖北工业大学 | 多相材料三维模型生成与连通性判断方法 |
CN116884540B (zh) * | 2023-06-15 | 2024-05-10 | 湖北工业大学 | 多相材料三维模型生成与连通性判断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108734751B (zh) | 2022-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108734751A (zh) | 一种基于随机分形的多孔介质细观结构模型生成方法 | |
CN109816049A (zh) | 一种基于深度学习的装配监测方法、设备及可读存储介质 | |
CN103745497B (zh) | 植物生长建模方法和系统 | |
CN107239595A (zh) | 基于nx软件pts模块的轮胎三维建模系统及其建模方法 | |
CN105354873B (zh) | 用于多孔介质三维重构的模式密度函数模拟方法 | |
CN112102494B (zh) | 骨架线引导的树状点云表面重建方法及装置 | |
CN105809733A (zh) | 基于SketchUp 的校园三维手绘地图构建方法 | |
CN109146120A (zh) | 一种输变电工程施工任务进度状态的三维模型标识方法 | |
Huijser et al. | Procedural natural systems for game level design | |
CN104616327B (zh) | 一种基于曲面细分的着色器简化方法、装置及图形渲染方法 | |
CN110757806A (zh) | 一种面向3d打印的双类形状密铺方法及系统 | |
CN109785427A (zh) | 利用二维图快速进行三维建模的方法 | |
Krs et al. | PICO: procedural iterative constrained optimizer for geometric modeling | |
Zhang et al. | Synthetic modeling method for large scale terrain based on hydrology | |
CN105678835A (zh) | 数字三维写意水墨山水画的建模绘制与渲染方法 | |
CN113486429B (zh) | 一种基于插值算法的空间汇交结构自动化建造方法 | |
CN106484511A (zh) | 一种谱姿态迁移方法 | |
CN107169218B (zh) | 一种纸张数字化信息植入的模拟方法 | |
CN104794750A (zh) | 一种基于空间殖民算法的树点云三维重建方法 | |
Montazeri et al. | A Practical Ply-Based Appearance Modeling for Knitted Fabrics | |
CN106780722B (zh) | 一种同种异态森林场景构建方法及系统 | |
CN113076571B (zh) | 一种三维衣物实时仿真编辑方法及系统 | |
Sun et al. | Engineering applications of 2D and 3D finite element mesh generation in hydrogeology and water resources | |
CN114445586A (zh) | 一种基于生成对抗网络的三维仿生设计方法及系统 | |
CN113722884B (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 |