基于裂隙网络模拟技术的结构面分布分维数计算
- 格式:pdf
- 大小:181.03 KB
- 文档页数:5
裂隙网络中流体的运移的模拟裂隙网络是一种丰富多样的地下水储层,它由许多微小的裂隙和孔隙组成。
裂隙网络在地下水资源的储存和运移中起到了重要作用。
裂隙网络中流体的运移模拟对于地下水资源的合理利用和地质灾害的预测具有重要意义。
本文将介绍裂隙网络中流体的运移模拟方法和相关应用。
裂隙网络中的流体运移模拟是基于流体力学和岩石力学等学科的研究成果,结合地下水流动和裂隙网络结构特点而发展起来的。
裂隙网络的结构特点使得流体在其中的运移过程十分复杂,需要多学科的知识和模型的综合应用。
裂隙网络中的流体运移模拟主要包括流体介质的数值模拟和物理模拟两种方法。
数值模拟是通过建立数学模型,采用数值计算的方法来模拟流体在裂隙网络中的运移过程。
数值模拟的基本思想是将裂隙网络划分为一个个小单元,建立各个单元间的质量和能量守恒方程,并采用数值方法求解这些方程。
数值模拟方法可以较好地模拟裂隙网络中的流体运移过程,但需要较高的计算能力和较长的计算时间。
物理模拟是通过实验室的物理实验来模拟裂隙网络中的流体运移过程。
物理模拟主要采用模拟岩石样品或者模拟模型来代替真实的裂隙网络进行实验。
物理模拟可以较真实地再现流体在裂隙网络中的运移过程,但需要耗费大量的时间和成本,并且受到模拟样品尺度和模型参数的限制。
裂隙网络中流体运移模拟的应用主要包括地下水资源的评价和管理、地下工程的设计和施工、地质灾害的预测和防范等方面。
通过流体运移模拟,可以评估地下水资源的储量和可持续利用量,为地下水资源的合理开发和利用提供科学依据。
在地下工程方面,流体运移模拟可以预测地下水对工程的影响,合理设计和施工地下工程,减少地下水对工程的破坏。
在地质灾害预测方面,流体运移模拟可以模拟裂隙网络中的流体运移过程,预测地下水位的变化和地下水的响应,为地下水灾害的预防和防范提供依据。
裂隙网络中流体的运移模拟方法和应用是一个复杂而重要的研究领域。
通过数值模拟和物理模拟,可以较好地模拟流体在裂隙网络中的运移过程,并应用于地下水资源评价、地下工程设计和地质灾害预测等方面。
页岩裂缝网络的几何特征二维表征及连通性分析第一章:绪论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数字模拟,能够真实、忠实地反映出页岩储层的内部结构和裂缝网络,为几何特征和流体特性方面的研究提供了途径。
4.3.结构面迹线测试方法结构面迹线长是描述结构面延伸程度、分布规模的重要参数,也是岩体工程稳定性分析的重要参数之一。
由于各类结构面迹线往往是随机分布的,通常的做法是,通过抽样测试和统计分析,确定其分布模型和特征参数,再模拟出工程开挖面的迹线分布网络,可为岩体工程的选址、开挖、支护和稳定性分析提供重要依据。
目前国内外常用的迹线抽样测试方法有测线法和统计窗法。
下面分别介绍这两种方法任期统计分析结果。
4.3.1. 测线法(Prist 1981)在岩体临空面上(工程围岩表面、岩体露头等)布置测线(取样线),逐条量测那些与测线相交的迹线长度和产状。
由于抽样测试面的限制,测出的可能是全迹线长度,也可能是半迹线的长度,还可能是未见迹线两端点的部分迹线长度,如图4.7所示。
为便于统计分析,首先明确各各迹线的定义及其测取范围。
即:图4.7 结构面迹线的测线测试法(1)结构面迹线长度l :结构面与临空而交线的实际长度.即那些两端点均位于临空面内的迹线长度;(2)平均迹线长度l :同一组结构面实际迹线长度的平均值;(3)结构面迹线端点:迹线的消失或中止点;(4)迹线端点密度μ:平均迹线长度的倒数,1lμ=; (5)半迹线长度h l :测线与迹线交点至迹线一端点的距离;(6)删截半迹线长度x l :测线与迹线交点至删截线与迹线交点之间的距离。
其中,删截线为在测试面上与测线平行的辅线。
通过统计分析得到各种类迹线平均长度与迹线长度期望值的关系,由测到的平均值就可以由下列公式估计样本所在区域迹线的平均值。
迹线长度倒数样本均值μ和方差σ与测试平均值g μ的关系:21g μσμμ=+ (4.8)删截半迹线长度倒数测试平均值1xl 与迹线长度(倒数)样本均值μ和方差σ的关系: 1ln()x n r l nμ-=- (当迹线长度为负指数分布时) (4.9)x μ= (当迹线长度为正态分布时) (4.10)式中,n 和r 分别为半迹线长度和删截半迹线长度的样本数。
天然裂缝影响下的复杂压裂裂缝网络模拟一、本文概述随着石油工业的发展,复杂压裂技术已成为提高油气采收率的重要手段。
然而,在实际压裂过程中,天然裂缝的存在往往会对裂缝扩展产生显著影响,使得裂缝网络的形成变得极为复杂。
因此,对天然裂缝影响下的复杂压裂裂缝网络进行模拟研究,对于优化压裂设计、提高油气采收率具有重要的理论和实践意义。
本文旨在通过数值模拟方法,深入研究天然裂缝对复杂压裂裂缝网络的影响。
我们将对天然裂缝的几何特征和分布规律进行详细分析,建立符合实际地质条件的天然裂缝模型。
在此基础上,我们将利用先进的数值计算方法,模拟复杂压裂过程中的裂缝扩展、交汇和融合等现象,揭示天然裂缝对裂缝网络形成的影响机理。
我们还将分析不同压裂参数下裂缝网络的演化规律,为优化压裂设计提供理论依据。
本文的研究内容将涉及地质建模、数值计算、数据分析等多个方面,综合运用了多种学科的知识和方法。
通过本文的研究,我们期望能够为复杂压裂技术的进一步发展和应用提供有力支持,为石油工业的可持续发展做出贡献。
二、天然裂缝对压裂裂缝网络的影响天然裂缝是地质构造中普遍存在的特征,它们对压裂裂缝网络的形成和扩展具有显著影响。
天然裂缝的存在,一方面可以作为压裂液流动的潜在通道,促进裂缝的扩展;另一方面,也可能成为压裂裂缝的屏障,限制裂缝的进一步延伸。
因此,在压裂模拟中,准确考虑天然裂缝的影响至关重要。
在压裂过程中,天然裂缝与压裂裂缝的交互作用会导致裂缝网络的复杂性增加。
当压裂裂缝遇到天然裂缝时,可能产生多种情况。
一种情况是压裂裂缝可能沿天然裂缝延伸,形成更长的裂缝,从而提高储层的连通性。
另一种情况是压裂裂缝在天然裂缝处发生转向,形成分支裂缝,增加裂缝网络的密度。
天然裂缝的开启程度和方位角也会对压裂裂缝的扩展产生影响。
为了准确模拟天然裂缝对压裂裂缝网络的影响,需要采用先进的数值模拟方法。
这些方法能够考虑天然裂缝的几何特征、力学性质和流体流动特性,从而更准确地预测压裂裂缝网络的形成和扩展。
裂隙网络中流体的运移的模拟裂隙网络是指由许多连通裂隙组成的岩石或土壤体系,这些裂隙构成了复杂的网络结构。
在这种介质中,流体的流动受到裂隙结构的限制,流体在不同大小、形状、方向的裂隙中反复流动和扩散,形成多尺度、多种流态的输运过程。
裂隙介质中流体的运移模拟对于石油开采、地下水资源利用、污染物运移等领域有着重要的应用。
裂隙介质的物理特性裂隙介质普遍存在于地球内部,而且广泛应用于石油勘探、水文地质工程、核废料储存等方面。
裂隙介质的物理特性包括裂隙几何特征、孔隙度、渗透率、孔隙分布等。
其中,渗透率指的是单位时间内单位面积上液体通过岩土体系的流量,其大小取决于裂隙的数量、大小、连通性以及液体在裂隙中的流动特性。
此外,裂隙介质中的流体运移行为受到多种因素的影响,如地下水压力、温度、溶质浓度、流体黏度等。
裂隙介质中流体的运移是一个多尺度、多相、多孔介质的过程,因此需要采用适当的模拟方法进行。
常用的模拟方法有数值模拟和实验研究两种。
数值模拟方法数值模拟是通过建立裂隙介质的数学模型,利用计算机进行数值计算,模拟裂隙介质中流体的运移过程。
数值模拟方法主要包括:有限元法、有限差分法、边界元法、离散元法等。
其中,有限元法是最常用的数值模拟方法之一,在此方法中,裂隙介质的实际形态被离散化为很多小的元素,再利用模拟程序模拟流体在裂隙中的运动。
数值模拟方法的优点是可以模拟不同尺度的裂隙,而且可以进行多尺度、多相现象的模拟,能够提供裂隙介质中流体的动态信息和运行规律。
但这种方法需要考虑到裂隙介质的各种特性,建立准确的数学模型,并在计算中充分考虑非线性、非平衡、非稳态等要素。
实验研究方法实验研究是通过建立实验方法,通过实验来研究裂隙介质中流体的运动规律。
常见的实验方法有压汞法、稳态方法、非稳态方法、核磁共振等。
其中,压汞法是一种广泛应用的实验方法,可以通过测量岩石或土壤样品中渗透率来研究裂隙介质中的流体运移规律。
实验研究方法的优点是设备和技术相对简单,经济实用,能够模拟裂隙介质中的流体运移规律,但需要考虑到实验条件对实验结果的影响,并克服试验现象不够复杂、实验过程难以模拟真实情况等限制。
岩石表面裂纹的计盒维数的matlab算法分形几何的产生源于两个数学问题,一是函数的可微性问题.另一个是维数问题,它的基本概念是分数维数,数字特征是幂律关系,几何或物理特征是自相似性。
它适用于描述自然界中一切不规则的、杂乱无章的现象,并提出了一种定量化的描述方法。
分形几何是一门以不规则几何形态为研究对象的几何学。
自然界中普遍存在着不规则现象,所以分形几何学又可以被称为描述大自然的几何学。
分形几何学建立以后,引起了各个学科领域的关注。
分形几何不仅在理论上,而且在实用上都具有重要价值。
分形理论是法国数学家Mamlelbort在20世纪70年代创建的一门研究学科,主要以不规则、复杂无序而又存在某种内在联系的几何形态为研究对象,与传统几何中的整形概念相对,是在传统几何空间基础上的一个重大提升。
它可以更加准确地描述自然界中物体、现象的特点、本质。
Mamlelbort提出:自然界中大多数几何对象都以分形的形式存在,不存在绝对意义的整形。
自相似性是分形理论的重要特性,系统的自相似性是:根据不同的时间尺度和空间尺度,自然界中结构、过程的特征都存在相似性。
具体表现在,某系统(某结构)的局部性质与整体性质相似,另外,部分与部分间也会存在自相似。
自相似的表现形式相对复杂,并不是对局部进行一定倍数的放大后和整体全部重合,这改变的只是系统的外部表现形式,系统(结构)的形态、复杂性、不规则程度,不会因为放大或者缩小等行为而产生变化。
分形可以依照自相似性的程度,来划分为:有规分形和无规分形,有规分形是指描述对象的自相似性程度非常高,能够通过数学模型来描述;无规分形是具有统计学意义上的分形。
维数是几何中重要的参量,用来表示空间中某一点所处位置所用到的独立坐标数目。
传统的欧式几何中一般采用整数来表示维数,例如:直线用一维表示,平面用二维表示,空间用三维表示。
分形几何中维数不一定用整数来表示,可以是分数。
虽然分形对象变化多样,但它们有着共同的特征存在,也就是可以通过分形维数来表征研究对象的不规则程度。
基于二维裂隙网络模拟的岩块搜索与岩层追踪方法李爱华;朱江【摘要】脉动压力在岩石河床裂隙中的传播是导致大坝下游基岩冲刷破坏的主要原因之一.岩体裂隙网络作为水体的赋存空间,其结构特性决定了水流脉动压力波的传播规律.应用蒙特卡洛模拟方法,编译Fortran程序语言,模拟得出与天然岩体裂隙网络在统计上完全等效的、具有相似结构特性的仿真裂隙网络,实现管网线元和节点的自动编号存储.借助矩阵理论定义了衔接矩阵和回路矩阵,从数学意义上精确描述裂隙网络.提出了基岩冲刷破坏过程中,孤立岩块的自动搜索和逐层追踪的新方法.文中以10 m× 10 m范围内的两组裂隙为例,依据裂隙网络的统计分布规律,通过蒙特卡洛模拟,给出由两组裂隙组成裂隙网络的模拟结构图以及块体逐层自动追踪的结果.研究结果为高坝下游基岩冲刷过程的模拟研究奠定基础.【期刊名称】《水利水运工程学报》【年(卷),期】2014(000)006【总页数】6页(P65-70)【关键词】裂隙网络;蒙特卡洛模拟方法;块体搜索;逐层追踪【作者】李爱华;朱江【作者单位】中国石油大学(华东),山东青岛266580;中国石油大学(华东),山东青岛266580【正文语种】中文【中图分类】TU45岩体裂隙中脉动压力的传播是导致高坝下游冲刷破坏的主要原因之一。
自1973年的国际大坝会议以来,各国学者对脉动压力在缝隙中传播并导致水力冲刷破坏进行了大量的试验研究和数值模拟[1-7],但这些研究大都局限在简单一维线裂隙或者二维面裂隙内,并没有考虑实际岩体复杂的裂隙结构特性。
事实上,岩体是由岩块和分割岩块的裂隙网络所组成的结构体。
岩体裂隙网络作为水体的赋存空间,其特殊的结构特性—如结构面密度、空间分布及其粗糙度、开度、填充情况等,从很大程度上决定了岩体裂隙介质中水流脉动压力传播的特殊性和复杂性。
国内潘别桐、陈建平等在统计学和概率论基础上提出了结构面的二维、三维网络模拟技术—即蒙特卡洛模拟方法,现已被广泛应用于裂隙渗流研究领域。
第23卷 第20期 岩石力学与工程学报 23(20):3465~3469
2004年10月 Chinese Journal of Rock Mechanics and Engineering Oct.,2004
2003年5月19日收到初稿,2003年7月6日收到修改稿。 作者 荣 冠 简介:男,32岁,1992年毕业于长春地质学院地质系,现任讲师、博士研究生,主要从事工程地质和岩体稳定性方面的教学与研究工作。E-mail:rg_mail@yeah.net。
基于裂隙网络模拟技术的结构面分布分维数计算 荣 冠 周创兵 (武汉大学水利水电学院 武汉 430072)
摘要 介绍了分形基本概念及工程岩体分形特性,结合裂隙网络模拟技术,分析了计算结构面分布分维数的方法,并给出了工程应用实例。讨论了计算裂隙分维数有效测度范围问题,对裂隙间距、迹长、交角和岩体RQD等影响因素与分维数的关系进行了初步分析,结果表明:间距对分维数影响最大,迹长和交角则相应次之。岩体裂隙分布分维数较RQD更能全面反映岩体的结构特征,用于对工程岩体进行质量评价是可行的,但对于测度尺寸和分维数与岩体工程性质关系等问题有待深入研究。 关键词 工程地质,分形维数,裂隙,网络模拟 分类号 TU 452 文献标识码 A 文章编号 1000-6915(2004)20-3465-05
STUDY ON FRACTAL DIMENSION OF DISCONTINUITY DISTRIBUTY BASED ON SIMULATION TECHNIQUE OF DISCONTINUITY NETWORK
Rong Guan,Zhou Chuangbing (School of Water Resources and Hydropower Engineering,Wuhan University, Wuhan 430072 China)
Abstract The fractal concept and the fractal characteristics of rock masses are discussed in general. The box dimension method is studied based on the technique of network simulation,and a calculation sample of fractal dimension for the engineering rock masses in Three Gorges Reservoir region is given. The problem of proper range of box size is discussed. Various factors,such as fracture space,total length of discontinuities,and inclination between joints,are analyzed. The results show that the space of discontinuities affects fractal dimension markedly. It is certain that the index of fractal dimension is better than the index,RQD,on characterizing jointed rock masses and classifying engineering rock masses. At last,some problems on the range of box size and the relationship between fractal dimension and engineering properties of rock masses are pointed out for further study. Key words engineering geology,fractal dimension,discontinuity,network simulation
1 分形概论与岩体分形特征 1.1 分形概念 1973年B. Mandelbrot首先提出分形(fractal)概念,1975年其专著《分形,机遇和维数》的出版是分形理论诞生的标志[1]。分形理论的诞生为揭示呈现纷繁复杂现象但隐藏着精细及相似结构的事物和对其进行定量研究提供了有效方法。分形可分为线性分形(自相似性分形)与非线性分形。目前研究最为广泛的为自相似性分形。自相似性是指局部与整体在形态、功能、信息、时间及空间等方面具有统计意义的相似性,在一定的观测尺度范围内分形对象的结构具有相似性。实际上自然界中的分形不存在具有严格数学意义上的自相似性,往往只具有统计意义的自相似性,而且这种相似性并不是在任何• 3466 • 岩石力学与工程学报 2004年 尺度上都能成立,通常只是在某些特定的尺度范围内才成立。 1.2 工程岩体的分形特征 20世纪80年代以来,分形几何开始应用于研究工程岩体的有关问题,这一现代非线性科学理论在解决如岩体类的复杂事物时,在一定程度上发挥了经典理论所不能发挥的作用,且应用越来越广泛深入[2]。实际上工程岩体的裂隙网络并不严格地遵循自相似性,而只是一种统计意义的自相似或局部一定尺度上的分形结构。判断系统是否为分形结构首先可根据Mandelbrot对分形定义来判断:分形是Hausdorff-Besicovitch维严格大于拓扑维的集合。如果系统的分维数与拓扑维相一致,则这个系统只是欧几里德系统或非分形系统;如果分维数大于拓扑维则为分形系统。其次可由分形的直观定义判断:分形是某种程度上,部分与整体相似的形态。分析系统是否具有自相似性或尺度不变性,如果系统具自相似性,则为分形系统,否则为非分形系统。 工程岩体中的结构面是地质历史的产物,其形成与地质历史环境紧密相关。岩体结构面的不少性质往往具有模糊性、不确定性和随机性等特征,但通过研究发现,岩体结构面的长度分布、粗糙性、张开度、间距分布和空间分布等特征在一定的尺度范围内表现出具有自相似性的规律,即反映岩体结构面存在分形特征。例如统计裂隙间距一般服从负指数或对数正态分布,根据文[3]研究指出间距服从以上分布规律只是一定程度、层次和发展阶段的表现,而在根本上服从分形分布。 分形理论在研究岩体裂隙、岩石破碎、岩体结构、地下水渗流等复杂问题中显示出了较强的能力,然而分形理论中关键参数——分维数在工程岩体中进行实测往往存在一定难度。本文主要结合裂隙网络模拟技术研究岩体裂隙分布分维数的计算方法,分析有关因素对分维数的影响和探讨分维数在工程领域的应用。 2 裂隙岩体网络模拟技术 从上世纪80年代初开始研究裂隙网络模拟技术以来,至目前基于Monte-Carlo原理的模拟技术已逐渐成熟,成果已广泛应用于解决岩体工程各种实际问题。裂隙网络模拟研究过程一般包括以下4个环节: (1) 野外采样:现场实测岩体中裂隙的几何参数。一般需要统计的参数有:裂隙组数、产状、间距、迹长、张开情况和结构面粗糙度等。常见的确定裂隙几何参数的采样方法有测线法和统计窗法。 (2) 分组和统计:在现场采样的基础上对裂隙样本进行分组统计研究,包括对样本进行分组和各组样本随机变量(如走向、倾向、倾角、间距、迹长等)的统计。 (3) 对样本分布形式进行拟合检验,判断各随机变量的统计分布形式及分布参数。 (4) 根据裂隙各随机变量的统计分布模型,生成符合裂隙分布规律的随机数,并以此生成裂隙网络图。 常规的裂隙网络统计模拟方法已为大家所熟悉,对具体技术在此不再介绍。但常规模拟技术仍存在不少值得改进和探讨的地方。如对结构面分组可采用聚类方法[4],它具有客观具体的特点。另外在建立结构面概率模型时,通常假定结构面各要素的概率分布函数为正态分布、对数正态分布、负指数分布及均匀分布,这与实际情况是有出入的。在实际应用中,当结构面符合以上4种分布规律时可按Monte-Carlo方法产生随机数,而不符合上述分布规律时可采用直接法[5]产生随机数。 应该指出Monte-Carlo随机抽样方法解决了已知概率分布的随机变量抽样的数学问题,在分布模拟中占有十分重要的地位。而对于未知概率分布的随机变量抽样采用直接法可得到较高精度的结果。 应用计算机生成岩体裂隙网络与实际裂隙分布之间的差异主要取决于现场测量、概率模型选择及计算机模拟方法等过程。只要正确地把握每一步,最终得出的模拟结果是能够近似反映岩体裂隙分布真实状态的。
3 分形维数及计算方法 分形几何的最重要概念是分形维数(fractal dimension),简称分维,最早由Hausdorff于1919年提出。实际研究中主要采用容量维数、信息维数、相关维数3种分维数。其中容量维数(也叫盒维数)的定义为:设F是平面上的有界点集,总可以找到一个矩形使F包含于其中。将矩形分割为若干个边长为ε 的小方格,计被F占领的小方格的数目为)(εN,则可定义F的容量维数为
)/1ln()(lnlim)(0cεεεNFd→= (1) 第23卷 第20期 荣 冠等. 基于裂隙网络模拟技术的结构面分布分维数计算 • 3467 • 容量维数的定义不仅适用于平面点集,同样适用于直线和空间点集。式(1)提供了近似计算点集容量维数的方法。计算F集容量维数的具体算法为:先构造一边长为a的盒子覆盖F集,然后变换不同边长值ε对应形成若干个小盒子,计算包含有F点集的小盒子数)(εN,经过多次变换可得到一系列ε-)(εN数据。再作)/1ln(ε与)(lnεN关系的散点图,采用最小二乘法求其斜率,即为容量维数。这就是通用的计算维数的覆盖法。在研究岩体结构面分布分维时,若考虑对象为沿某一测线方向结构面迹线交点的分布集,此时即为线裂分维。若考虑对象为某一测面上的结构面迹线分布情况,此时即为面裂分维数。同理,还存在体裂分维。 4 工程应用实例 选择的工程实例为三峡水库某库段香溪河谷两岸岩体。三峡水库蓄水后该库岸将被淹没,原城镇将就地后靠重建。为给该库段地质灾害防治工程设计和施工提供依据,进行了此库段岩体裂隙调查及网络模拟研究。图1和表1为研究区内岩体裂隙统计模拟结果。 图1 研究区岩体裂隙网络模拟图 Fig.1 Simulated network of discontinuities in studied area 在裂隙网络模拟图的基础上可以较方便地进行结构面分维数的研究。本文主要研究了库岸工程岩体裂隙的线裂分维和面裂分维,采用的计算方法为上面介绍的确定容量维数的覆盖法。研究线裂分维时,首先生成40 m×40 m大小露头的裂隙网络模拟图;选定0°,15°,30°,…,150°和165°等12个计算线裂分维数的平面方位(参见图1),计算出每个方位上裂隙交点数和位置,这样便可求出在平面内不同方位的线裂分维数。计算面裂分维数则根据模拟出的20 m×20 m区域网络图按覆盖法计算。实际计算分维数只要对原统计模拟程序稍作修改,加上一段在不同条件下生成)/1ln(ε与)(lnεN数对和拟合其直线方程的程序段即可。线裂分维的程序算法十分简单,而确定面裂分维的程序算法中某一小盒子是否被裂隙占领的方法为:首先判断在盒子内是否存在裂隙的端点,如存在,则此盒子被占领,否则继续判断裂隙是否与盒子的4条边相交,如有一边与裂隙相交则盒子被占领,如均未相交则盒子未被占领。根据上面方法研究区岩体裂隙分布分维数的计算结果为:在32/1(~L)024 1/1(L为研究范围内裂隙平均迹长)测度范围,岩体线裂分维数沿平面不同方位为0.284 1~0.376 5范围内变化。在(1/1 024~1/65 536)A (A为裂隙模拟范围的面积)测度范围内岩体面裂分维数为1.492 5。