基于统计模型的裂隙岩体渗流场的并行数值模拟.pdfVIP

基于统计模型的裂隙岩体渗流场的并行数值模拟 刘耀儒，杨强，覃振朝 摘要：为模拟天然岩石中裂隙对渗流的影响，基于裂隙的空问形态分布及分布的统计特性，应用裂隙网络随机 模拟技术，采用Matlab程序生成二围裂隙网络．由此网络生成包含裂隙的有限元计算网格．考虑到渗流场数值 编制了相应的模拟分析并行程序，并在分布存储的并行机上实现．分别对包含一组裂隙和两组裂隙的岩体进行 了渗流场的数值模拟．计算结果表明。该方法能够较真实的反映裂隙岩体的实际状态，对于模拟岩体中的裂隙 是很有效的；同时，并行计算可有效降低数值模拟的计算时问． 关键词：裂隙岩体；数值分析；渗流；并行计算 1引言 的方程： 在边坡和坝肩的稳定分析、地下洞室开挖和稳定 分析以及油藏模拟方面，考虑渗流场的影响是非常有 必要的．由于天然岩体存在着大量的孔隙和裂隙，这 灿oct唣Nj警NldQ+驱p胆+ 些缺陷不但大大地改变了岩体的力学性质．也严重影 响着岩体的渗透特性。因此，对天然岩体中裂隙的模 管脚釉姜》磅州巧嵩移，，芳矽Q 拟是十分有必要的． 在过去的20a中，裂隙岩体的渗流场模拟得到了 =＼rq，NIdr 很大的发展．许多学者已经做了大量的卓有成效 (2) 的工作，如Louis、Rome、BarenbIatt、oIssonll】、 杜广林12J、王嫒【3l、仵彦卿嘲、王思志14l、张有天 和陈平‘51等． 其中，M和M均为形函数；吼为边界的流量． 皇舌墅套篓鎏誓魁坌錾譬荽蓑苎i三重奎，垦落 2．2描述裂隙参数的方法 委悬萋囊菖蓑羹蜚豢嘉鬈藿％畜蠢蓄孟鍪霪譬 “‘套：蒜翥洋姜署晶：妥至妻特征是岩体中存在着 勰鬈案嫠釜霎彗藩戮船黧晕釜爨 大量淼甏麓鞴志赫茹簇鬈蒹 羹翟薷套眚嬖i釜嚣是冈爨蠹羹桀妻舞譬象嘉喜柔；；i凝翟蒿笔嘉惩荔i；i薹；辫蒿赛薯譬蓄襞 雾燃峰辫鬻髦善燃掣 磊弼善客黑黧黧幕婴嬲蔫掣 龋黧三二嚣粼燃黼雾嬲茧勰嬲竿矽哩一删暇‘ klJ彤饫2裂隰喇彤恹一舣抽砸刀l：自】疋1、殳‘ 娶嚣；篙警善藉。尝嚣擎篇妻嚣蚕蠢；瓣旮尖簟零蕊簇Zi翥；靠占蒜—般描述为固定不变． 法进行数值模拟．编制了相应的并行模拟程序， 度来描述：走向和倾角．裂隙的产状可以用概率 随机数的最大值和最小值；刀用来表征数值积分精 分布来描述，常用的概率分布有均匀分布、 度．一般当门8时就可以满足精度要求：仃和∥分 飚her分布、双变量Fisher分布和Binghanl分 别为统计量的方差和均值。 布等。 (6)开度：统计表明，裂隙的开度的分布符合 3．2迹线在岩体中的描述方法 对数分布或者是指数分布． 迹线在岩体中的分布是随机的．要在岩体中生成 一条迹线，就要先确定它由哪几个随机参数所控制．‘ (7)裂隙的边界：在实际观测中，可以区分两 假定迹线长是一条线段，则描述其在岩体的位置时可 种不同裂隙的边界t一是裂隙终止于与另一裂隙 的相交处：另一种是裂隙终止子完整的岩体． 以用迹线的长度，，迹线中心点坐标(％，％)以及 迹线的倾角方向曰来完全确定．如图l所示．而在实 际中。还要考虑裂隙的张开度．所以迹线在平面上是 2．3裂隙岩体渗流场的模拟步骤 有一定的厚度的．即通常所说的力学隙宽． (1)采用测线量测法或其它方法确定岩体的裂 隙参数： (2)建立各种参数