地应力对采动裂隙岩体渗流特性影响的研究

《地应力对采动裂隙岩体渗流特性影响的研究》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

地应力对采动裂隙岩体渗流特性影响的研究本文由jiefengfei贡献第26卷第9期2009年9月长江科学院院报Vol.26No.9Sep.2009JournalofYangtzeRiverScientificResearchInstitute文章编号:1001-5485(2009)09-0035-04地应力对采动裂隙岩体渗流特性影响的研究孙玉杰1,邬爱清1,张宜虎1,武涛2(1.长江科学院水利部岩土力学重点实验室,武汉430010;2中国地质大学工程学院,武汉.430074)摘要:岩体总是赋存于一定的渗流场和应力场中的,在裂隙岩体中开挖洞室将打破最初的平衡状态,引起渗流场和应力场发生变化,因此研究地应力场对渗流场的影响具有重要的意义。采用UDEC程序,利用数值试验的方法,研究了块体边界尺寸以及初始地应力对在节理裂隙岩体中开挖洞室引起隧洞内涌水量变化的影响。结果表明,随着块体边界尺寸和初始应力比的增加,开挖后隧洞内的涌水量均减小。关键词:地应力;地下洞室;渗流;节理岩体文献标识码:A中图分类号:TU459.41概述张后全、杨天鸿、等应用FRFPA渗流与应力耦合分析系统,对煤层顶板随着开采的逐步进行,采动裂隙逐渐向上发展并最终与含水层连通,进而发生顶板突水的全过程进行了数值模拟,直观地显示了煤层顶板的变形、破坏过程以及渗流场在整个岩体的运移过程及其突水前后渗流场的变化情况,较好地揭示了顶板突水过程。随着国民经济建设的迅猛发展,基本建设规模的不断扩大,人类的活动空间越来越大,使得其开发所依托的岩土体问题也越显突出,裂隙介质地下水运动规律就是其中的突出问题之一。诸如大型水利水电工程中的地下水问题、核废料的深部贮存、垃圾填埋污水下渗、输油管道渗漏造成的地下水污染问题等都与裂隙岩体的渗流有关。而岩体是受结构面的纵横切割、具有一定结构的多裂隙体,同时赋存于一定的地质环境中,地应力场和地下渗流场是岩体物理力学环境中的2个重要组成部分。裂隙岩体渗流受应力的影响,而渗流场的改变又会影响到应力场的改变,两者相互影响、相互作用称为耦合作用。由立方定理可知,裂隙开度的变化对渗流的影响很大,而应力会直接引起裂隙开度的变化肖江、高喜才[4][3]2D率。[1]等研究了采动过程中因顶板冒落、放顶跨落造成顶板透水的灾害机理,并结合具体工程地质条件,利用固-液耦合相似材料模拟方法对首采工作面综放开采过程进行了研究,分析了裂隙带与冒落带导水通道分布特征,总结了含水层潜水渗流规律。黄涛、立杨[5]中在裂隙围岩介质渗透性能等效处理的基础上,提出了渗流与应力耦合环境下裂隙围岩隧道涌水量预测计算的确定性数学模型方法, 包括理论解析法、经验解析法和水文地质数值模拟方法。马立强、张东升[6]。因此,等运用FLAC数值软件的应力渗流耦合系3D研究应力对渗流的影响具有重要的意义。目前,国内外学者多集中研究天然岩体的水力学特征,而对于受扰动后的岩体渗流特征研究较少。殷黎明、杨春和[1]统,结合工程实践,系统分析了采动岩体的渗流矢量、渗流速度、孔隙压力及导水边界线动态发展的普遍规律和分布特征,进行了顶板突水危险区域的识别。杨会军、王梦恕[7]等分析了裂隙岩体渗流运动的基本规律,研应用岩体非连续介质渗流的究了应力对渗流作用影响的机理,并结合甘肃某工程进行了现场应力测量及高压压水试验,得出了该地区渗透系数随地应力值的增大呈负指数递减的变化规律。陈伟、阮怀宁[2]研究成果,结合具体工程进行了最大涌水区段渗流数值模拟,模拟了各种工况条件下地下水的渗流特性以及地表水体的水位变化。隧洞作为地下建筑物,修建过程中将不可避免地穿越不同水文地质、工程地质条件的裂隙岩体。当隧洞通过裂隙岩体的含水区段时,由于隧洞的开利用自己编制的随机连续模型导入有限元程序进行了2场耦合分析,预测了锦屏二级水电站隧洞开挖过程中地下水的最大入渗收稿日期:2008–1112;修回日期:2008–1209––基金项目:国家重点自然科学基金项目(50509003)作者简介:孙玉杰(1982–),男,山东临沂人,硕士,主要从事岩石力学与工程方面的研究工作,(电话)13954637126(电子信箱)sdustsy@j163.com。36长江科学院院报2009年挖人为破坏了原有地下水的渗流条件,使得隧洞自身成为地下水以不同形式(渗出、滴流、股流及大范围突水等)向外排泄的地下廊道,形成涌水灾害。对于大埋深的越岭隧道来说,裂隙岩体处于相对较高的地应力场中,同时穿越含水区段的裂隙岩体中有地下水渗流场的存在,与地应力场构成一种特定的动态平衡体系。因此,研究地应力对隧洞开挖后裂隙岩体中渗流的影响十分必要且很有意义。[9]Indraratnaet(1998)基于UDEC利用数值模型对裂隙岩体中开挖洞室所引起的洞室中涌水量开展了研究。而本文将采用数值模拟的方法讨论初始地应力比(水平地应力与垂直地应力的比值)与典型块体尺寸对裂隙岩体中开挖隧洞所引起的隧洞中总的涌水量的影响。[8]由流动域中中心压力的大小决定。按块体接触条件的不同,裂隙岩体中流体的速率有2种计算方法。(1)点接触。点接触分为角-边接触和角-角接触,设流动域的流体压力为p1,流动域#的流体压力为p2,则由流动域到流动域#流体的流速为v=-kCp,式中kC为接触处的渗透系数。式中: W(1)p=p2-p1+wg(y2-y1),(2)为流体密度;g为重力加速度;y1,y2为2个流动域的中心坐标。(2)边-边接触。首先定义接触长度,图1中lD和lE分别为接触点D和点E的接触长度,然后利用平行板裂隙中的立方定律计算流动速度,即q=-kja3p,l(3)为2UDEC基本算法本文研究的重点是尝试研究裂隙岩体中开挖隧洞时,地下洞室内总的涌水量随着块体尺寸以及地应力比值如何变化。在本文研究中,利用UDEC(IASTCA,1996)水力全耦合模型来模拟裂隙中水的流动。UDEC(UniversalDistinctElementCode)是IASCAT公司开发的针对非连续介质的平面离散元程序,在数学求解方式上采用了与FLAC一致的有限差分法。对于水力全耦合的渗流分析,当水流主要是由裂隙网络控制时,UDEC程序是非常合适的。因此本文中假设岩石基质是不透水的,水流主要是通过水力连通的裂隙网络涌入隧洞内。有些研究者,包括Herbert(1996),Zhang等(1996)和Hencher(1997)利用UDEC程序进行过裂隙岩体介质中的渗流模拟。UDEC中被裂隙所包围的岩块可以被模拟为刚体或者是可变形体。UDEC通过域分析流体在裂隙中的流动。图1中将域顺序标号为~,假定域内充满各向等压流体,域域之间通过接触与邻域发生作用。接触顺序标号为A~F。域、、表示节理,域#!∀表示2个节理的焦点,域为空洞。不计重力时,假设流体压力在流动域中的分布是均匀的。考虑重力时,流体压力则按线性分布的静水压力计算。流体的流动是由相邻流体域的压力差决定的,其中流动域中流体压力的大小图1通过域模拟流体在节理裂隙中的流动Fig1Flowinjointfissure.modeledbyflowbetweendoainsm式中:kj为裂隙的渗透系数(理论值为1/12),度;a为接触的水力开度。a=a0+un,流体的动力粘滞系数;l为2流体域之间的接触长其中:a0为裂隙在无法向应力时的开度,un为裂隙的法向开度(张开为正)。水力开度的最小值为ares,最大值为amax。裂隙宽度随节理法向应力的变化如图2,图中节理法向应力!n为正时表示节理处于受压状态。计算过程中,每计算一个时步,重新生成系统的几何形状,而后计算出所有接触的裂隙宽度以及所有域的体积,之后利用上面的公式计算出各接触处的流量。最后,再考虑流入域的流体流量综合以及由于周围块体的位移增图2UDEC中水力开度a与法向应力!n之间的关系Fig2Relationbetween.hydraulicaperture(a)andjointnoralstress(!n)minUDEC量而产生的域体积的变化,按下式计算出域内的孔隙压力,p=p0+kwQt/v-kwv/vm,(4)式中:p0为前一时步的孔隙压力;Q为通过孔隙周围的所有接触点流入该孔隙的流量之和;kw为流体的体积模量;v=v-v0,Vm=(v+v0)/2其中,v和v0分,别为现在时步和前一时步孔隙的体积;步。计算出域内孔隙压力后,可以计算流体作用在其周围岩块的力。将该力与诸如接触点力和外力荷载等力迭加,施加在块体的节点上。这样得到不透t为计算时第9期孙玉杰等 地应力对采动裂隙岩体渗流特性影响的研究37水岩块的总应力以及节理的有效法向应力。3数值试验3.1数值模型本模型考虑在地表面下100m深处开挖直径为6m的圆形隧洞,并假定地下水位在地表下10m处。模型中考虑节理是2组连续节理,倾角(与水平方向夹角)分别为30∃150∃(如图3所和示)。采用节理面接触模型以及莫尔库伦模型分别用于表征节理与完整岩块的物理响应。3.2模型岩体材料性能以及边界条件节理模型性质见表1。岩块、节理和孔隙流体的材料性能见表2。表1模型所采用节理参数Table1变量节理倾角/(∃)节理间距/m间隙长度/m节理迹长/m图4UDEC模拟中所用的边界条件Fig4Boundaryconditionsappliedto.inUDECmodel假定的边界条件(水力边界和应力边界)如图4所示。流体压力pyy=∀y,w图3含有两组节理的节理模型Fig3Jointmodelcontaining.wtocontinuousjointsets(5)式中:∀是水的重度,y是距离水位面的垂直距离。w最初是没有水涌入开挖区域的,开挖区域内是静水压力状态。一旦隧洞开挖后,在隧洞开挖表面就会变为常量的大气压力。最初的压应力状态是由重力引起的各向同性状态。!yy=式中:gggy,!xx=#!yy,1!zz=#!yy,(6)2为岩块密度;y为水平面下的垂直距离;!xx,!yy,!zz为x,y,z方向的初始应力;#1,#2为初始应力比值因子,并假定在分析中保持为常量。上面的数值试验与简单的平板之间的层流是明显不同的,平板之间的层流节理隙宽被假定为恒定的常量。当前的分析采用离散元法能够模拟出开挖引起围岩的变形对隙宽的影响,反过来,隙宽的变化又影响渗流和应力的重分布。在下一步的迭代过程中将当前的应力用于耦合的渗流变形分析。水力耦合分析分2个步骤进行:第一步是通过施加相应的边界条件模型在初始应力场下达到平衡状态;第二步是假设隧洞的开挖是瞬时完成的,之后模型再次达到新的平衡状态。在本文的研究中,对隧洞内涌水量的以下两个影响因素进行了模拟分析:不同的块体尺寸;#不同的初始应力比(垂直应力和水平应力)。3.0&10102.0&10101003210&10-3.30&10-4.12&1010.10&1010.1003120&10-3.40&10-4.10008&10-42&109流体(水)第2组节理15035.030.Jointparaetersofjointmodelm第1组节理304.005.0表2岩体和水的材料性质Table2Materialpropertiesofrockmassandwater材料变量岩块节理第一组节理第二组节理体积弹性模量/2.7&1010(N%m-2)剪切弹性模量/0.7&1010(N%m-2)岩密度/ 块2500(kg%m-3)粘聚力/672&106.(N%m-2)内摩擦角/(∃)27.0法向刚度/(N%m-2)切向刚度/(N%m-2)裂隙渗透系数/(Pa%s)摩擦角/(∃)初始隙宽/m残余隙宽/m密度/(kg%m-3)流体动力粘滞系数/(Pa%s)体积模量/(N%m-2)33数值试验结果分析.UDEC数值模拟分析的结果见图5和图6所示。渗流分析是采用UDEC的稳定渗流分析模块进行分析的。对于给定的节理形式、块体边界尺寸以及边界条件,隧洞内的涌水量是将与隧洞外围相交的裂隙流量相加得到的。数值分析的结果表明隧洞的涌水量是随着不同的参数变量而改变的。331块体边界尺寸的影响..隧洞内总的涌水量大小与各种不同块体尺寸之间的关系如图5所示,该图表明随着块体尺寸增大,隧洞内总的涌水量在减小。原因是随着块体尺寸节理38长江科学院院报2009年隧洞内涌水量也在减小,这与国内外学者的研究结果是相吻合的,为研究裂隙岩体的渗流规律分析打下了基础。但由于所采用的节理为贯通节理,与实际裂隙岩体中的情况还是有差距的,对渗流应力耦合的机理研究也还不够明晰,有待于进一步的研究。通过本文也证实了UDEC进行裂隙岩体耦合渗流分析是有效的。参考文献:图5块体边界大小对渗流量的影响Fig5Effectsofboundaryblocksizeonflorate.wfordifferentinsitustresses[2][3][4][1]殷黎明,杨春和,王贵宾,等.地应力对裂隙岩体渗流特性影响的研究[J].岩石力学与工程学报,2005,l7(24):3071-3075.陈伟,阮怀宁.锦屏二级水电站裂隙岩体渗流-应力耦合分析[J].水力发电,2008,34(1):25-28.张后全,杨天鸿,赵德深,等.采场工作面顶板突水的渗流场分析[J].矿井地质,2004,17(2):25-27.肖江,高喜才,马岳谭,等.富水覆岩采动裂隙渗流相似模拟研究[J].采矿与安全工程学报,200825(1):,50-53.[5][6]图6初始地应力比对渗流量的影响Fig6Effectofinsituinitialstressratioonflowrate.[7][8][9]黄涛,杨立中.渗流与应力耦合环境下裂隙围岩隧道涌水量的预测研究[J].铁道学报,199921(6):75-80.,马立强,张东升,缪协兴,等.FLAC3D模拟采动岩体渗流规律[J].湖南科技大学学报(自然科学版),2006,21(3):1-5.杨会军,王梦恕.深埋长大隧道渗流数值模拟[J].岩石力学与工程学报,2006,25(3):511-519.张金才,张玉桌.应力对裂隙岩体渗流影响的研究[J].岩土工程学报,199820(2):19-22.,INDRARATNAB,RANJITHPG,AZIZN.NumericalPredictionofInflotoUndergroundCavityUsingaCoupledwHydroGeoechanicalModel[C].InternationalConferencemonGeoechanics/GroundControlinMingandUndermgroundConstruction,DepartmentofCili,MiningandEwilronmentalEngineeringUniversityofWollongon,Australia,1998,(2):863–871.[10]ITSCA.UDECUniversalDistinctEleentCodeUser∋sm,Manual[K].MinneapdisUSA:ITSCAConsultingGroup,,1996.[11]INDRARATNAB,WANGJC.EffectsofStressChangeonWaterInflostoUndergroundExcavation[J].AustralwianGeoechanics200129:99-114.m,,[12]LIAOQH,HENCHERSR.NumericalModellingoftheHydroMechanicalBehaviourofFracturedRockMasses[J].117.(编辑:曾小汉)(下转第43页)IntJRockMech&MinSc.i1997,(3-4):... .增大与隧洞外围相交的裂隙内水压力减小,因此隧洞内的涌水量变小。332初始地应力比的影响..隧洞内总的涌水量与初始地应力比!h/!v之间的关系如图6所示。与估计的相同,随着初始地应力比增加隧洞内总的涌水量减小。随着初始应力比从05到225的变化,隧洞内涌水量的降低值达到..40%~80。这个结果与Indraratna和Wang%[11](1996)所做得的降低70是相一致的。Liao和%Hencher(1997)[12]利用UDEC对初始水平应力比、边界条件以及块体体积对裂隙岩体的渗透特征进行过模拟。他们的结果表明总的渗透系数随着初始水平应力与垂直应力的比值增加在减小。本文的分析结果与其分析是相一致的,不同之处在于本文所采用的模型以及利用涌水量代替了渗透系数。4结论本文在裂隙岩体开挖隧洞中采用数值试验,得到如下结论:随着初始水平应力与垂直应力比值的增大,隧洞内的涌水量在减小;随着块体尺寸的增加第9期崔德山等离子土固化剂加固红黏土的X射线衍射试验43[8]DREVERJIThePreparationofOrientedGlayMineral.SpeciensforXrayDiffractionAnalysisbyaFilterMemmbranePeelTechnique[J].AmericanMineralogist1973,,58:553-554.[10]MOOREDM,REYNOLDSRC.XrayDiffractionandtheldentificationandAnalysisofCayMinerals[M].London:OxfordUniversityPress1997,.[11]丁锐.水玻璃与黏土矿物之间表面反应的实验研究[J].岩土工程学报,1999,21(3):334-337.(编辑:曾小汉)[9]RUANChuandeWARDCR.QuantitativeXraypow,derdiffractionanalysisofclaymineralsinAustraliancoalsusingRietveldmethods[J].2002,(21):227-240.AppliedClayScience,XRayDiffractionAnalysisofRedClayReinforcedbyISSCUIDeshan,XIANGWei,JOACHIMRohn(1FacultyofEngineering,ChinaUniversityofGeosciencesWuhan.,112430074China,;91056Germany),2DepavtentofAppliedGeologyUniversityErlangenNuremberg,Erlangen.m,AbstractTheengineeringproblemsbroughtbythecontradictionbetweenthebadphysicalnatureandgoodme:chanicalpropertiesofredclayaredrawingmoreconcernintheacademicandtheengineeringcirclesMakinguseof.ionicsoilstabilizer(ISS)caneffectivelychangethephysicalpropertiesandmechanicalperformancesofredclay inHanyangdistrictofWuhanCityInordertostudythephasecompositionandthespacingofcrystalplaneofredclay.reinforcedbyISStheXraydiffractionofredclaywasdoneundertheconditionsofMg,2+saturatedK,+saturationglycerolsaturationandheatingto550(bytheSiemensD5000XraydiffractometerTheresultsshowedthat,.thephasecompositionofredclayreinforcedbyISShasnotsignificantvariationandthespacingofcrystalplane,hasasmallchange.Keywordsionicsoilstabilizerredclay;Xraydiffractionphasespacingofcrystalplane:;;;(上接第38页)EffectofGeostressonPereabilityofExcavatedFracturedRockMassmSUNYujie,WUAiqing,ZHANGYihu,WUTao1112(1KeyLaboratoryofGeotechnicalMechanicsandEngineeringoftheMinistryofWaterResourcesDivisionof.GeotechnicalEngineering,YangtzeRiverScientificResearchInstituteWuhan,2FacultyofEngineeringChinaUniversityofGeosciences,Wuhan.,430010China,;430074China),AbstractThegeostressfieldandseepagefieldalaysexistinrockmassthecavityexcavationinthefractured:w,rockmasswilldisturbtheinitialequilibriumandcausethechangeofgeostressandseepage,soitisiportanttomstudytheeffectofgeostressonthepermeabilityofrockmass.Theeffectofrockblockboundarysizeandinitialgeostressonthechangeofblowingdischargeintheexcavatedcavityortunnelisstudiedquantitativelybyanumericaltes.Itisfoundthatwiththeincreaseofblocksizeandtheratioofinitialhorizontalstresstoverticalstressthet,blowingdischargeintheexcavatedtunnelwillbedecreased.Keywordsgeostressundergroundcavityseepagefracturedrockmass:;;;