压力溶解和自由面溶解/沉淀对双重孔隙岩体中热–水–应力耦合影响的有限元分析

 引入裂隙开度的压力溶解和自由面溶解/沉淀模型,针对一个假设的位于饱和双重孔隙–裂隙岩体中的高放废物地质处置库,拟定3种计算工况：(1) 裂隙开度是压力溶解和自由面溶解/沉淀的函数;(2) 裂隙开度仅随压力溶解而变化(这2种工况中基岩的孔隙率亦是应力的函数);(3) 裂隙开度和基岩的孔隙率均为常数,进行热–水–应力耦合的二维有限元分析,考察岩体中的温度、裂隙开度及渗透系数、孔隙水压力、地下水流速和应力的变化、分布情况。结果表明：自由面溶解/沉淀引起的裂隙开张及渗透系数增量的绝对值要明显大于压力溶解产生的裂隙闭合及渗透系数减量的绝对值,而压力溶解对裂隙的开度及其渗透系数的影响较小;同时计入压力溶解和自由面溶解/沉淀相比于仅考虑压力溶解,同时计入压力溶解和自由面溶解/沉淀的裂隙开度及其渗透系数分别约为仅考虑压力溶解时的1.5和7.0倍;在释热温度场的作用下,计算域中的裂隙水压力随时间先上升再下降,但变化幅度不大;3种工况下岩体中的应力量值及分布差别很小。     

基于双重孔隙介质模型的渗流-应力耦合并行数值分析

为模拟天然岩石中孔隙和裂隙的不同渗透性能及应力耦合特性,基于双重孔隙介质模型,建立渗流–应力耦合分析的迭代计算模型。岩石基质渗透系数和孔隙率通过体积应变进行更新。裂隙系统的耦合模型通过立方体单元模型来建立,渗透系数和孔隙率随有效应力的变化进行更新。考虑到耦合分析时步多、计算量大的问题,采用基于element-by-element策略的有限元并行计算方法进行数值模拟。编制相应的耦合分析并行程序CoupledGF,并在分布存储的并行机上实现。对包含一个生产井和一个注入井的封闭区域进行渗流–应力耦合分析,模拟约22 000个时步。计算结果表明,并行耦合分析模型对于渗流应力的耦合分析是很有效的。     

基于双重孔隙介质模型的渗流–应力耦合并行数值分析

为模拟天然岩石中孔隙和裂隙的不同渗透性能及应力耦合特性，基于双重孔隙介质模型，建立渗流–应力耦合分析的迭代计算模型。岩石基质渗透系数和孔隙率通过体积应变进行更新。裂隙系统的耦合模型通过立方体单元模型来建立，渗透系数和孔隙率随有效应力的变化进行更新。考虑到耦合分析时步多、计算量大的问题，采用基于element-by-element策略的有限元并行计算方法进行数值模拟。编制相应的耦合分析并行程序CoupledGF，并在分布存储的并行机上实现。对包含一个生产井和一个注入井的封闭区域进行渗流–应力耦合分析，模拟约     

偏应力往返作用下饱和砂土变形特性的试验研究

在固结排水条件下保持球应力不变,进行了偏应力往返作用下的动三轴试验,分析了饱和砂土的体积变形特性。结果表明:(1) 体应变可视具体条件表现为C1型剪缩(初次加荷)、C2型剪缩(卸荷前未产生剪胀)、C3型剪缩(卸荷前产生剪胀)、剪胀P及剪刚R(再加载无胀缩)等物态有规律的变化;(2) 各物态的体应变均与所作用的偏应力比有良好的线性关系。其参数与综合反映应力特征及土性特征的各物态起始体应变间均有良好的指数或直线关系;(3) 体应变可跟踪偏应力的往返变化,区分加荷、卸荷段的各不同物态,依序考虑各物态起始时的体应变,在对应确定相应参数的基础上进行计算,计算的结果与试验关系间有良好的一致性。     

常压至高压下中砂剪切特性及应力–剪胀关系

为研究高围压范围内砂土相对密实度和围压对土体强度和变形特性的影响,对3种不同相对密实度砂土试样在常至高围压下进行常规三轴固结排水剪切试验,获得偏应力–轴向应变–体应变关系曲线,同时进行颗粒破碎分析。结果表明:在常至中压范围(0.8 MPa≤σ_3≤2 MPa),应力–应变曲线均表现出不同程度的应变软化,其剪胀性随相对密实度增加和围压的降低而增强;当进入高围压范围时(σ_32 MPa),应力–应变曲线逐渐向应变硬化型转变,试样体积逐渐趋于剪缩。颗粒破碎程度随着围压和密实度的增大而增大,在高围压时由于中密和密砂剪切后期出现了明显的颗粒破碎,导致剪切过程中出现了二次相变。不同密实度土体的破坏内摩擦角和对数围压表现良好的线性关系,拟合确定了破坏内摩擦角随对数围压增加的衰减率,同时基于Bolton应力–剪胀关系拟合确定了试验砂土的临界状态内摩擦角,建立了剪胀指标与初始相对密实度及平均有效应力的关系式,为高压情况下砂土地基稳定性分析等提供强度参数。     

用渗透性随埋深变化趋势 反求岩体变形模量的简化模型

 取得大尺度范围岩体的宏观力学参数是地质工程领域的难题。渗透性随埋深衰减的趋势恰恰反映岩体在大尺度范围的渗流场–应力场耦合行为,为解决这个问题提供契机。根据自重应力作用下岩体的应变–渗透耦合机制,以孔隙度和裂隙宽度为中间参数,建立利用渗透性随埋深变化趋势反求岩体变形模量的简化模型,推导了岩块体积模量和结构面法向刚度的反演公式。这种简化模型是在假设岩体物质结构相对均匀、只存在单组裂隙、无构造应力的理想情况下建立的,但为获取大尺度范围岩体的宏观应力–应变参数提供一种新的参考途径。     

基于双参负指数分布的改进的岩体本构关系

在结构面迹长、间距服从双参数负指数分布以及对迹长进行统计区间[α1,α2]积分的基础上,结合室内岩石力学实验,建立改进的岩体统计本构模型。通过引入剩余剪应力比h,建立结构面产状与地应力之间的函数关系,量化主地应力方向、岩体结构面产状对岩体力学性质各向异性的控制性作用。研究结果表明:当地应力方向发生变化时,将改变剩余剪应力比h值,进而引起岩体应力–应变曲线发生变化;在此改进的岩体统计本构模型基础上,计算云南小湾水电站坝基岩体的应力–应变曲线,揭示该坝基岩体在不同应力条件下的各向异性特征。     

基于全应力–应变曲线及起裂应力的岩石脆性特征评价方法

脆性作为岩石重要的力学指标,准确评价岩石脆性特征对岩体工程稳定性具有重要意义。目前,常用的岩石脆性指数评价方法主要基于岩石压缩过程中峰后应力–应变曲线变化规律,而对岩石峰前应力状态关注较少。综合考量峰后应力跌落速率及起裂应力值至峰值应力值之间的应力增长速率,建立一种基于岩体应力–应变全过程的脆性指数计算方法,更加准确、合理地描述岩体脆性特征。试验结果表明:单轴条件下当全程采用轴向位移控制时,脆性较强的岩石难以准确获得其峰后曲线,容易降低基于峰后应力曲的常规脆性评价方法的准确性,而所建立的脆性指标可根据峰前曲线进行判断,能够准确评价岩石脆性程度。在三轴条件下,该脆性指标更为全面地反映了围压对大理岩脆性的抑制作用。建立的脆性指标基于岩石整个压缩过程的应力状态,能够更加显著地体现不同因素对岩石脆性特征的影响。试验结果很好地验证了该指标的准确性与优越性,研究结果对丰富与完善岩体脆性评价体系具有重要的作用。     

温度-应力耦合条件下膏岩渗透性试验研究

为研究深地环境下能源储库膏岩盖层气密性,利用基于四川大学MTS815岩石力学试验系统改进的THM多场耦合三轴渗流测试平台,国内首次展开不同温度、不同围压条件下膏岩加载变形破坏全过程气体渗透性试验。试验结果表明:(1)天然状态下膏岩的渗透率处于10-17～10-15 m2,属于致密低渗岩体介质。(2)相同温度下,随围压增大,膏岩破坏应力–应变曲线具有脆–延性转变特征,且膏岩峰值应力不断提高;相同围压下,随温度升高,膏岩表现出热损伤及热软化效应,峰值应力逐渐降低,塑性变形能力增强。(3)膏岩初始渗透率随围压增大而降低,呈指数型递减关系,5 MPa围压下膏岩初始渗透率为4.13×10-17 m2,15,25,35 MPa围压下初始渗透率分别降低37.23%,65.86%,75.79%。(4)不同温度–应力耦合条件下膏岩渗透性均表现出先降低后升高的演化规律,岩石整体渗透率水平随温度及围压的升高不断降低。     

渗流–应力–化学耦合作用下岩石裂隙渗透特性试验研究

 为研究渗流–应力–化学耦合作用下岩石裂隙渗透特性变化规律,设计3组试验工况,在改变渗透压以及化学溶液的条件下,分别测定每种工况下的渗出水流量、渗出水离子浓度值以及渗出水pH值变化情况,进而得出裂隙渗透特性变化情况。通过处理试验数据,总结分析各因素对裂隙渗透特性的影响,并建立裂隙开度变化率与渗出水中钙离子浓度值之间的关系式。试验结果表明,渗出水流量、裂隙开度总体趋势是随着时间逐渐减小的,最终趋于稳定状态;增大渗透压,稳定状态会被打破,裂隙的流量和开度都会增大,但最终趋于另一个稳定状态;化学溶液对岩体裂隙的侵蚀性大,对岩体渗透性的影响更明显。通过分析和提炼渗出水流量、裂隙开度、渗出水的离子浓度值以及渗出水的pH值等随时间变化的数据,及它们之间的内部关系,在理论上描述岩体裂隙在渗流–应力–化学耦合作用下的渗透特性,进一步揭示渗流–应力–化学耦合作用机制。     

CT尺度砂岩渗流与应力关系试验研究

岩石渗流与应力关系研究是进行岩石渗流场与应力场耦合分析的关键。运用岩石高压三轴加载装置和渗透压加载装置,对砂岩进行了渗流与应力关系试验,同时借助SOMATOMPLUS螺旋CT扫描机进行实时观测。通过试验结果分析,推出了基于CT数的岩石孔隙率公式,在此基础上,分析了岩石应力–应变过程中孔隙率、渗透速度、渗流速度、微孔隙直径、渗透率等的变化规律。结果发现:岩石渗透参数的变化与岩石受力损伤–破裂过程密切相关。在初期的压密阶段,岩石的孔隙率、渗透速度、渗流速度、微孔隙直径、渗透率等随应力的增大而减小;当岩石的内部出现微裂纹后,岩石的孔隙率、渗透速度、渗流速度、微孔隙直径、渗透率等随应力的增大而增大,从宏观应力–应变关系看,从微裂纹出现到宏观破坏出现前,岩石还处于弹性变形阶段;当岩石宏观破坏时,岩石的孔隙率、渗透速度、渗流速度、微孔隙直径、渗透率等达到最大值。同时还发现:在渗透水压力作用下,受压砂岩的微裂纹起裂应力占岩石峰值强度的45%,而同样干岩样中微裂纹起裂应力占岩石峰值强度的55%以上,也就是说,渗透水压力使砂岩样的强度损失10%。     

剪胀对地下工程岩体位移的影响——以加拿大Donkin-Morien隧道为例

根据建立的岩石剪胀角模型,分析岩石峰值内摩擦角和剪胀角的关系,得出岩石在零围压时的峰值剪胀角小于并近似等于峰值内摩擦角,并假设岩石和岩体的剪胀角遵循相似的变化趋势,结合Hoek-Brown强度准则和GSI岩体分级系统,实现剪胀角模型从完整岩石到岩体的转化。采用程序语言在FLAC3D中编写岩体剪胀角模型程序模块。以加拿大Donkin-Morien隧道为工程实例,研究围压和塑性剪切应变依赖的岩体剪胀对隧道渐进开挖过程中围岩位移的影响,论证恒定的剪胀角值不能准确表达隧道开挖边界附近的岩体位移,而考虑围压和塑性剪切应变为影响因素的岩体剪胀角模型能够合理描述围岩的位移分布,模拟结果与实际测量值具有很好的一致性。研究成果可为岩体非线性力学行为的研究和地下工程岩体的稳定性控制提供理论和实践基础。     

基于梯度塑性理论的岩样单轴压缩扩容分析

采用梯度塑性理论，对岩样剪应变局部化引起的扩容进行了理论分析。假设岩石的剪切本构关系为弹性-应变软化双线性，局部化启动于应力峰值强度，利用局部塑性剪应变与局部塑性体积应变的线性关系，得到了局部塑性体积应变、局部塑性体积增量及剪胀引起的剪切带总塑性体积增量的解析式，这体现了该理论在研究剪胀问题时的优越性。另外，还得到了弹性阶段及应变软化阶段的轴向应力-体积应变曲线的理论关系。塑性体积应变是专指由剪切带剪胀而引起的，因而，轴向应力．体积应变不具有尺寸效应，与局部化带的尺寸无关，但扩容角、剪切降模量及泊松比却对该曲线有重要影响。在弹性阶段及应变软化阶段轴向应力-体积应变均呈线性。在相同的应力水平下，扩容角越大则剪胀程度越大；剪切降模量越大，剪胀程度越小。在应变软化阶段，泊松比不影响塑性体积应变。     

考虑应变软化、剪胀和渗流的水工隧洞解析解

基于统一强度理论和弹脆塑性软化模型,考虑不同工况下主应力顺序、岩石应变软化、剪胀和渗流等综合影响,推导了圆形水工压力隧洞应力和洞壁位移解析解。通过工程算例,比较本文方法与传统方法的不同,说明考虑第一主应力的变化和渗流体积力更符合工程实际,并得出统一强度理论参数、软化特性参数和剪胀特性参数对隧洞切向应力和洞壁位移的影响特性。隧洞设计时,应依据围岩材料特性合理地选取统一强度理论参数,并考虑由于软化特性参数和剪胀特性参数的减小对隧洞切向应力和洞壁位移增大的影响。该结果为圆形水工压力隧洞弹塑性分析提供了理论依据,对工程设计有一定的参考价值。     

高孔岩石中局部变形带的理论 和形成条件研究进展

 近几年在野外和实验室都发现高孔隙岩石中局部化变形的压缩带、剪切带和膨胀带,是目前岩石力学、岩石物理和本构模型研究的热点之一。基于局部变形带在水利、石油、核废料处置、垃圾处理、环境污染的治理和地质构造、工程地质和岩土工程中的实际意义,介绍以局部化分岔理论为基础的多孔岩石的压缩带、剪切带的理论、形成条件和判别条件,以及与压缩带、剪切带和膨胀带相应的临界硬化模量。基于帽子模型,介绍椭圆形屈服帽上的压缩和剪切带,以及它们与s-e曲线之间的密切关系。在低围压下的s-e曲线上的第一个零模量点也就是应力平台对应于压缩带,即为帽子模型的体积屈服面;继续加载使s-e曲线继续上升到标本产生硬化,直到s-e曲线上出现第二个零模量点,也就是s-e曲线的峰值对应于剪切破坏,此时应力状态处于帽子模型的剪切屈服面上。同时还给出大量不同颗粒尺寸、不同孔隙度砂岩轴对称压缩试验结果,它们也可用1/4椭圆帽子模型来描述。     

考虑围压影响的深埋圆隧围岩应变软化与剪胀特性分析

 不同围压下岩石应变软化与剪胀特性不同,若在隧洞开挖中考虑围岩塑性区域内变化围压影响,其应力–应变场求解方式将区别于既有文献中的传统方法。根据围压影响下应变软化围岩的临界塑性剪切应变变化特征,给出改进的判断围岩是否进入塑性残余区域的规则;引入考虑围压与临界塑性剪切应变的非线性剪胀模型。基于Hoek-Brown屈服准则,根据一定径向应力增量将围岩塑性软化与残余区域分层,采用有限差分法对围岩应力–应变场进行求解;为分析围压对围岩稳定性的影响,根据临界塑性剪切应变与剪胀系数变化与否,设定4种非线性力学模型,深入分析并比较4种力学模型下临界塑性剪切应变、剪胀系数与围岩变形等在塑性软化与残余区域的分布规律。研究结果表明：地质强度指标GSI较小时,考虑围压影响下的围岩应力–应变场与未考虑时差异明显;此时临界塑性剪切应变的减小对开挖边界的围岩剪胀性具一定抑制作用。     

渗流水压力分级加载岩石蠕变特性研究

采用RLW-2000M微机控制煤岩流变仪,以细粒砂岩为研究对象,对三轴压缩条件下岩石渗流水压力分级加载蠕变试验进行了蠕变特性研究。重点分析了渗流水压力分级加载时蠕变条件下岩石的应变、渗流体积(体积速率)演化曲线,同时对不同渗流水压力分级加载条件下的岩石蠕变与渗透系数演化曲线进行了分析和对比。试验结果表明:分级加载渗流水压力作用下细粒砂岩的蠕变曲线符合蠕变演化三阶段特征;随着渗流水压的逐级加载,瞬时轴向、横向应变和瞬时泊松比呈增大的趋势;随着应变的累积,横向应变量大于轴向应变量呈扩容效应,直至发生蠕变破坏;分级加载渗流水压力作用下渗流体积曲线呈线性演化。研究认为:渗透系数先瞬时减小后增大,初始渗透系数都具有记忆性,即岩石孔隙通道具有记忆特征,孔隙通道经过变形→闭合→冲蚀→形成新通道过程。