浅谈离散裂缝网络(DFN)模型与裂缝随机建模方法

致密砂岩储层天然裂缝建模方法林承焰; 李辉; 马存飞; 任丽华; 陈仕臻; 李师涛; 梁书义【期刊名称】《《中国石油大学学报（自然科学版）》》【年(卷),期】2019(043)005【总页数】13页(P21-33)【关键词】构造裂缝; 致密砂岩; 应力场模拟; 离散裂缝网络; 密度体; 岩石力学【作者】林承焰; 李辉; 马存飞; 任丽华; 陈仕臻; 李师涛; 梁书义【作者单位】中国石油大学(华东)地球科学与技术学院山东青岛266580; 山东省油藏地质重点实验室山东青岛266580; 中国石油大学(华东)深层油气重点实验室山东青岛266580; 胜利油田东胜精攻石油开发集团股份有限公司山东东营257000【正文语种】中文【中图分类】TE122.2天然裂缝的成因受多种因素的综合影响,其发育状态具有很强的不确定性,分布规律具有很强的随机性[1-4]。

如何更加准确地刻画储层中裂缝的空间展布,仅从构造或者岩相单一手段分析裂缝的发育程度存在明显缺陷[5-6]。

从裂缝发育的岩石来看,不同部位岩石脆性组分的含量、岩石结构、岩石粒度等的差异导致对裂缝控制作用的不同,即沉积相控制下的岩性为主要影响因素[7-8];从裂缝形成机制来看,采用有限元数值模拟构造应力场分析方法,认为应力分布情况是主要影响因素[9-14];从裂缝与周围构造体间关系来看[15-16],裂缝发育位置距断层远近为主要影响因素;从地震资料入手,利用蚂蚁追踪、相干体属性等分析方法对识别裂缝具有很好的效果[17]。

从建模手段来看,逐渐从连续裂缝模型发展到离散裂缝网络模型[18-21]。

针对致密砂岩储层,笔者采用基于示性点过程的离散裂缝网络建模方法,在裂缝发育地质模式的控制下,综合裂缝成因力学机制,利用裂缝密度体建立天然裂缝模型。

1 裂缝建模研究思路传统裂缝建模方法对成像测井、岩心的依赖程度较大,存在成本高,井间误差大等问题。

裂缝密度体的提出可以有效解决缺乏特殊测井和取心井,以及井少、井距大的情况下裂缝预测与建模的问题。

页岩裂缝网络的几何特征二维表征及连通性分析第一章：绪论1.1 研究背景和意义1.2 国内外研究现状和进展1.3 研究内容和目标第二章：页岩裂缝网络二维表征方法2.1 裂缝网络特征参数2.2 图像处理方法与算法2.3 页岩裂缝网络二维表征模型第三章：页岩裂缝网络的几何特征分析3.1 裂缝长度分布规律分析3.2 裂缝面积分布规律分析3.3 裂缝角度分布规律分析3.4 裂缝密度与分形维数分析第四章：页岩裂缝网络连通性分析4.1 连通性参数的定义4.2 网络连通性的数学模型4.3 连通性分析实例及结果解释第五章：总结与展望5.1 涵盖的问题及创新点5.2 研究成果与展望5.3 存在的问题与改进方向附录：1. 图像处理算法介绍及源代码2. 统计分析软件介绍及用法第一章：绪论1.1 研究背景和意义近年来，页岩气作为一种新型的能源资源备受重视。

而在页岩气的开采中，掌握页岩储层的特征和裂缝网络是关键的技术之一。

页岩气储层的生产能力与裂缝网络的数量和规模密切相关，因此对页岩裂缝网络的研究成为了页岩气领域的研究热点。

页岩是一种含油气的泥岩，其孔隙度极低，储层不易形成，因此页岩储层的裂缝网络对气体流动和储集具有不可替代的作用。

裂缝网络的结构特征和空间分布决定了气体的流动路径和有效地渗透面积，在页岩储层的开采过程中，不仅会形成不同性质的裂缝网络，而且裂缝之间还有很强的复杂交互作用。

因此，合理地掌握裂缝网络的几何特征、规律和连通性对于评估页岩储层的开发潜力，并指导页岩气开采的灵活性和效率具有重要的意义。

1.2 国内外研究现状和进展页岩裂缝网络的研究主要集中在其几何特征和流体特性方面。

国内外学者尝试从不同角度深入研究页岩裂缝网络，形成了丰富的研究成果。

其中，基于三维数字化数据的几何表征方法和样品分析技术是研究页岩裂缝网络的重要手段。

例如，基于Ct扫描的3D数字模拟，能够真实、忠实地反映出页岩储层的内部结构和裂缝网络，为几何特征和流体特性方面的研究提供了途径。

202019 年 第 6 期采用分离裂缝模型模拟混凝土裂缝扩展时，有限元网格需要不断地调整和再生成，这将导致新的节点拓扑的产生。

为此，本文探讨了基于虚拟裂缝模型的混凝土裂缝扩展模拟方法，通过非线性弹簧模拟虚拟裂缝的黏聚力，旨在避免裂缝扩展过程中有限元网格的不断剖分，提高计算效率和精度。

三点弯曲梁算例结果表明，本文所探讨的方法是合理和有效的。

自20世纪60年代以来，国内外学者在计算机模拟裂缝扩展方面做了很多工作，采用了[1]很多新方法。

有限元法作为强大的分析工具已经得到学术界的广泛认同。

断裂力学中的虚拟裂缝模型对裂缝扩展模拟采用分离网格节点来实现，在裂缝扩展过程中需对单元网格进行不断剖分，这将导致新的节点拓扑的产生，在计算效率和计算精度上都存在一定的缺陷。

本文在分析混凝土断裂过程区材料的软化特性基础上，运用非线性断裂力学模型描述其本构关系，探讨基于虚拟裂缝模型的裂缝自动扩展模拟方法，旨在避免裂缝扩展过程中有限元网格的不断剖分，提高计算效率和精度。

1 虚拟裂缝模型在混凝土宏观裂缝前端存在一个微裂缝区，随着荷载的增加在自由裂缝和微裂缝之间存在裂缝的亚临界扩展。

裂缝的亚临界扩展和微裂缝区合称为断裂过程区。

断裂过程区的存在减小了材料刚度，削弱了材料本身应力传递能力，这种现象就是混凝土的软化特性。

虚拟裂缝模型把裂缝分解为两部分：真实物理裂缝和虚拟裂缝。

前者表征宏观自由裂缝，裂缝面上不传递任何应力，而后者将带状裂缝区简化为一条分离裂缝，即虚拟裂缝。

虚拟裂缝面上可以传递应力，其上某点应力的大小与该点裂缝面的张开位移 之间存在确定的关系。

虚拟裂缝模型通过非线性断裂过程区在混凝土类准脆性材料的断裂分析中引入了基于应变软化机理的非线性本构模型，合理有效地对裂缝由微裂纹向宏观裂缝发展进行了描述，并且虚拟裂缝还表征了由于裂缝张开而出现的非线性连续位移。

裂缝张开位移与拉应力的软化关系是虚拟裂缝的核心内容。

为了数值计算的方便，一般将软化曲线简化成图1示的单直线模型或双直线[2]模型。

如何运用UDEC创建裂隙岩体水力学模型-水力学论文-水利工程论文-水利论文——文章均为WORD文档，下载后可直接编辑使用亦可打印——1、引言库岸边坡常因库水位变动而失稳，尤其是岩石里面的节理裂隙，由于水在裂隙中流动，影响岩石的各种特性，控制着岩体的变性破坏特征。

本文简单论述一下如何运用UDEC建立斜坡离散裂隙网络数值模拟计算模型，可以更加正确的反映岩石内部的结构和构造。

从而可以帮助我们揭示一些裂隙岩体边坡在各种水的作用下的下变形机制，为一些岩质边坡的防治提供科学的理论依据。

2、离散单元法的基本原理离散单元法是1970年由Cundall首次提出的，于1986年由王永嘉引入我国，是专门针对不连续介质问题提出的数值模拟解决方法，它对于边坡稳定性的研究是将所研究的边坡岩土体划分为一个个小块，通过每一个小块间的相互作用，以及力与位移的相互作用建立方程。

通过一次次的迭代，配合所建立的平衡方程，使每一个小块都达到平衡状态。

由于离散单元法是通过计算块体之间的作用得到的结果，所以这种方法可以分析实际岩块间大位移的情况，而且可以详细的解析出岩体内部应力与应变的分布情况。

它还有一个重要特点，既其求解平衡方程是利用时间差分法。

因此该方法在实际工程中可以弥补有限元法的缺点，进而求解非均质和不连续体的大位移和大变形的问题。

2.1离散元程序UDECUDEC（Universal Distinct Element Code）是一款由ITASCA公司基于离散单元法原理开发并推广应用的二维的大型商用数值模拟软件。

UDEC对模拟节理化岩体材料介质在准静态及准动态荷载条件下的反应过程特别合适，它不但能够实现接触的模拟，而且可以自动侦测并识别新的接触产生，并模拟其力学行为。

UDEC数值分析程序是为一系列工程问题开发的专业求解工具，例如：它可以应用于地下结构、地震、矿山、核废料处理、能源等问题的研究。

2.2裂隙岩体离散裂隙网络介质模型研究岩石中有很多断层、节理、裂隙，统称为结构面，在岩石水力学中都称之为裂隙。

裂缝三维地质建模的难点与对策聂永生;田景春;魏生祥;孙利【摘要】裂缝三维地质建模通常包括裂缝强度建模、离散裂缝网络建模和裂缝属性建模3个主要步骤,在这3个步骤中面临着如何将裂缝信息和成果归一化以及如何三维量化描述裂缝开度和渗透率属性参数等诸多难点.对上述问题进行了初步探讨,并总结出具体的解决方法.首先,利用成像测井或常规测井裂缝解释的数据作为基本数据点,以动态分析及地震检测等资料为约束条件,选用合适的地质统计学方法建立裂缝强度模型；其次,综合利用岩心、测井及动态监测等研究成果,建立离散裂缝网络模型；最后,运用岩心裂缝分析和经验公式相结合的方式获得初步的裂缝属性模型,并根据试井分析等方法获得的裂缝渗透率等参数进行调整,得到最终的裂缝属性模型.运用上述方法结合国际主流的裂缝研究与建模软件,可定量描述裂缝的三维分布,为油藏数值模拟及油藏研究提供更好的支持.%Three-dimensional quantitative analysis of fracture is not only the comprehensive embodiment of fracture research results, but also the objective need of fractural reservoirs study and numerical reservoir simulation. With the advance of logging recognition, seismic detecting, dynamic monitoring, physical modeling and numerical simulation methods, we can get more valuable materials than before which bring even higher challenges in fracture three-dimensional modeling. Three main steps included in the fracture three-dimensional quantitative analysis, i.e., fracture intensity modeling, DFN modeling and fracture network properties modeling. It is the most important and the most difficult thing that how to unify the fracture information and research results, how to calculate fractureattributes such as aperture and permeability. This paper has proposed the related solutions and discussions, introduced a method which uses imaging or conventional logging data, dynamic analysis results, seismic detecting data and other data to build fracture intensity model, and uses core fracture data in combination with empirical formulas to build fracture properties models. Combining this method with international fracture modeling software, we can quantitatively describe the three-dimensional distribution of fracture and provide solid support for reservoir study and numerical simulation.【期刊名称】《油气地质与采收率》【年(卷),期】2013(020)002【总页数】3页(P39-41)【关键词】裂缝;蚂蚁体;裂缝强度模型;离散裂缝网络模型;裂缝属性模型【作者】聂永生;田景春;魏生祥;孙利【作者单位】成都理工大学沉积地质研究院,四川成都610059;成都理工大学油气藏地质及开发工程国家重点实验室,四川成都610059【正文语种】中文【中图分类】TE112.3随着中国的裂缝性油藏和受裂缝影响的油藏在已发现油气藏中所占比例逐渐增大，其储量和产量所占比例也逐渐上升，针对储层裂缝的研究日益受到重视［1-2］。

基于Monte-Carlo方法的三维裂隙网络模拟摘要:目前,人们研究裂隙岩体中流动和溶质运移时越来越多的采用离散介质方法,因此需建立离散裂隙网络。

由于裂隙的分布具有不确定性和交叉性,目前无法用确定性的研究方法进行研究。

文章用概率论与统计学的方法分析裂隙的分布及组合特征,利用Monte-Carlo方法进行三维裂隙网络随机模拟,建立了三维裂隙网络的随机模拟模型。

该研究可为进一步分析裂隙网络的连通性和渗透张量分布提供依据。

关键词:离散裂隙网络;概率统计;Monte-Carlo;三维裂隙网络目前裂隙岩体中流动和运移模拟方法可分为三类:①把岩体看作等效连续介质体,建立等效连续介质渗流数学模型;②不考虑岩块中的渗流,而把裂隙作为非连续网络处理,建立裂隙网络非连续介质渗流数学模型,即离散裂隙网络方法;③考虑岩体内裂隙导水,岩块储水而建立的岩体双重介质渗流数学模型。

其中,连续介质模型运用最为成熟,但这种方法因把裂隙空间平均到渗流单元体上,不能很好地刻画裂隙的特殊导水作用,在实际工程中,当岩体的裂隙分布比较稀疏,且岩体中的渗流主要取决于大的裂隙时,运用连续介质方法处理岩体渗流问题则容易引起大的误差。

离散裂隙网络模型最能反映裂隙发育的规律及裂隙水渗流的规律,但是受目前测量手段和计算量的限制,不可能把大范围的研究区内的每一条裂隙都测出来,目前大多只应用于小范围研究。

现阶段,我国主要用的是广义双重介质模型,即对主干裂隙采用离散裂隙网络模型,对于次要裂隙采用连续介质模型。

双重介质模型,结合了上述两种模型的优点。

国外更多的是研究非饱和带裂隙中水的流动和溶质运移,主要与废物处置、油气运移有关,且研究尺度小,研究精度高,一般采用离散裂隙网络模型或双重介质模型。

但是无论是离散裂隙网络模型还是双重介质模型,离散裂隙网络的构建是其基础。

1离散裂隙网络的构建离散裂隙网络方法的前提假设为:提供流体流动的路径的裂隙网络的几何参数(包括位置,形状,方位,尺度和每条裂隙的水力传导系数等)能够用数学描述。

2012年6月7号裂缝综合预测方法及应用研究目前国内外裂缝预测技术大致有4类[1]:①性分析和生产经验总结的裂缝预测方法;②利用测井手段和地震信息识别和预测裂缝;③非线性理论方法检测和识别裂缝,主要应用分形理论和神经网络等技术方法;④通过构造应力来分析预测裂缝；应用Petrel进行构造建模的研究由于Petrel软件不能自动识别断层的接触关系,所以必须对断层的接触关系进行处理。

在处理断层接触关系时应遵循如下原则:1裂系统发育史要明确,早期断层不能切割晚期断层。

2主、从断层要分清,主断层不允许被剪切。

3断层间的接触关系一定要清晰,断面可以相交但不能互相切割。

基于离散裂缝网络模型裂缝储层的随机模拟在裂缝预测研究方面已有学者做过大量的工作，例如，Ｍｕｒｒａｙ［３］提出了曲率法，通过构造主曲率来预测裂缝储层的孔隙度与渗透率；曾锦光等［４］提出了利用屈曲薄板挠曲理论来模拟地层古应力，对裂缝的分布进行预测：模拟得到的最大（或最小）主应力正值区为拉张区（涉及到强度因子），最易产生拉张破裂缝，因而模拟所得主应力越大，理论上的破裂就越发育；王仁［５］从构造应力方面考虑，结合岩石破裂准则对裂缝分布规律进行了数值模拟；丁中一等［６］提出了“二元法”，即破裂法和能量法的结合运用，破裂值代表发育裂缝的可能性，能量法代表裂缝发育能力的大小，应用于裂缝分布预测；宋慧珍等［７］将构造应力场和地质研究紧密结合，运用有限元数值模拟方法，研究裂缝空间状态；刘卫丽等［８］运用地震资料正反演结合对裂缝进行了预测。

上述学者在研究裂缝储层时，都是通过构造应力分析，结合岩心、露头的观测、各种测井曲线分析、地震资料的利用、生产动态资料和相似区块的开发经验，定性半定量确定裂缝的表征参数，根据这些参数从整体概略预测某个区域的储层裂缝发育程度。

由于Ｐｅｔｒｅｌ建立裂缝模型有局限性，而Ｆｒｅｄ能较好地建立裂缝模型，因此，在本次模拟中采用３步建模：①运用Ｐｅｔｒｅｌ建立基质的物性模型；②利用Ｆｒｅｄ建立裂缝物性的空间分布特征；③将区块的基质模型和裂缝模型利用网格数据叠加方法建立整个裂缝渗流区的地质模型。

裂缝性油藏数值模拟方法姚军（中国石油大学山东东营 257061）摘要：目前对天然裂缝性油藏的数值模拟可以大致分为连续性模型和离散性模型两大类；连续性模型又可以分为双重介质模型和单介质模型，双重介质模型主要是以Barrenblatt 和Warren-Root在20世纪60年代提出的双重孔隙/双重渗透模型为基础，在这类模型中认为油藏中每一点都存在有基岩和裂缝两种介质，基岩被相互平行排列的裂缝分割称为单个的岩块，每种介质存在独立的水动力场，通过两种介质间的窜流的将其联系起来；而对于单介质模型，则是通过一定的方法将裂缝的渗透率和基岩的渗透率进行综合的考虑，得出整个油田的有效渗透率，该有效渗透率考虑了裂缝的密度、方位等的影响，然后将该有效渗透率输入到普通的单一介质模拟器中来对裂缝性油藏进行模拟；由于双重介质模型不能够对不连续且控制着流体流动的大裂缝进行准确的模拟等原因，离散性模型在近段时间逐渐发展起来，而其又可以分为离散裂缝网络模型和离散管网模型；在离散裂缝网络模型中，对地质上描述出来的每个裂缝都进行了离散的显式的表示，同时根据局部裂缝的形状决定基岩的几何形状，由于地质上描述的裂缝数目一般较多，相应的在数值模拟中需要的离散点数目也就十分巨大，对模拟造成了一定的困难，所以目前很多的专家和学者又对该方法进行了进一步的改进，有许多简化的方法存在；离散管网模型则是先对所要模拟的区域进行了网格的划分，进而采用管子连接两个网格块，相应的两个网格块之间的传导率也采用管子的传导率来代替，这种方法的特点是数学上比较简单，灵活性较强，同时由于管子只对其连接的两个网格有影响，所以改变管子的传导率只会影响一个方向的传导性，而不会像常规的模拟器那样要同时影响两边的传导性，但是该方法目前研究较少。

0 前言随着世界碳酸盐岩油气田的大规模开发，系统深入研究这类油气田的渗流模式及其在开发中的应用已成为重要课题。

地质学家通过岩芯分析，确认碳酸盐岩(灰岩、白云岩)具有明显可见的裂缝、孔洞，含有密集的树枝状构造的粗裂缝以及连接的孔洞和孔隙。