层状盐岩体非线性蠕变本构模型

 针对由不同岩层交替而组成的层状盐岩体,建立泥岩夹层和盐岩复合体代表单元,根据应变协调原理,从细观力学分析角度通过考虑泥岩夹层弹性性质、盐岩弹性及蠕变力学特性以及两相体积含量建立层状盐岩体宏观各向异性非线性蠕变增量型本构模型,分析层状盐岩在蠕变过程中因保持细观应变协调而产生的应力重分布问题,给出新本构模型ABAQUS有限元二次开发增量迭代算法实现方法,对一层状盐岩体简单试样算例进行初步计算分析验证,分析计算与试验结果吻合较好,表明该本构模型能反映细观应变协调时层状盐岩体宏观蠕变力学特性,该模型将为层状盐岩体内硐室长期稳定性分析提供理论计算基础。     

考虑温度和体积应力的分数阶蠕变损伤Burgers模型

盐腔是石油、天然气和高放射性核废料深地储存的理想屏蔽围岩,考虑到深地盐岩地层中温度和体积应力联合影响,基于细观概率体元强度的威布尔(Weibull)分布构建损伤变量,利用Drucker-Prager强度准则描述微观线弹性脆性破坏,引入考虑塑性损伤的广义胡克定量描述样品尺度损伤,从而建立三轴压缩条件下静态本构模型.进一步...     

基于组合模型法的贯通节理岩体动态损伤本构模型

针对贯通节理岩体动态变形特点并结合已有岩石动态本构模型的相关研究成果,将贯通节理岩体变形过程中的动态应力视为贯通节理岩体静态应力分量与相应动态应力分量的叠加。其中贯通节理岩体静态应力分量采用考虑岩石细观损伤的非线性元件、节理面闭合及剪切变形元件等3个基本元件的串联来模拟,动态应力分量采用黏性元件来模拟,从而建立了贯通节理岩体动态单轴压缩损伤本构模型。其次,根据贯通节理岩体在单轴压缩荷载下往往会沿节理面发生剪切破坏的特点,在前述已建立的损伤本构模型中引人节理剪切破坏准则对该模型进行修正,从而更好地考虑了节理剪切强度对该模型的影响,最终建立了考虑节理剪切强度的贯通节理岩体单轴压缩损伤本构模型。最后利用该模型对贯通节理岩体在压缩荷载作用下的力学特性进行了分析计算,重点讨论了节理倾角对岩体单轴动态压缩峰值强度的影响规律。研究结果表明随着节理倾角的变化,节理岩体将发生岩块张拉或剪切破坏、沿节理面的剪切破坏及上述两种破坏模式的复合破坏,相应地节理岩体的单轴压缩动态峰值强度也随之有较大变化。     

冻融循环作用下砂岩的力学特性及细观损伤本构模型研究

针对寒区岩体工程中岩石的冻融问题,选取砂岩为试样,通过进行室内冻融循环试验、扫描电子显微镜观测和三轴压缩试验对砂岩质量损失、微观结构和力学特性进行了分析。然后基于Lemaitre应变等效假设理论,通过引入能够反映岩石冻融破坏过程中的细观冻融损伤变量和力损伤变量来描述岩石材料的劣化程度及损伤演化规律,并采用连续损伤力学理论,建立了冻融与围压耦合作用下岩石的损伤演化方程及细观损伤本构模型。采用理论推导的方法得出所需的模型参数表达式,最后利用冻融岩石的三轴压缩试验数据对该模型的合理性和准确性进行了验证。将试验曲线的峰值点与模型理论曲线的峰值点进行对比,结果表明两者吻合度较好,该损伤本构模型能够较好地反映岩石三轴压缩过程的应力-应变峰值特性,验证了该模型及模型参数确定方法的合理性与可靠性。该模型拓展了岩石在冻融与围压耦合作用下的损伤模型,进一步的揭示了岩石在冻融与围压耦合作用下的损伤机制和破坏规律。     

层状岩体多节理本构模型与试验验证

 根据含多节理层状岩体典型力学特征,在所提出的考虑一组节理面的层状岩体复合材料模型基础上,建立考虑多组结构面特性的层状岩体多节理本构模型,用来描述其在强度和变形方面的各向异性以及应变硬化–软化特征,进而采用VC++编程语言将该模型嵌入到FLAC3D软件中,实现其非线性数值计算功能。然后,利用所提出的本构模型与三维节理网络模型相结合,初步提出基于连续介质力学理论的合成岩体模型(SRMM)分析方法。基于这些模型理论和分析方法,建立现场真三轴试验岩样的三维合成岩体计算模型,通过与现场试验结果对比,表明所提出的本构模型在描述多节理岩体力学特性方面是可行的、合适的。     

软岩的损伤西原流变模型及FLAC3D实现

流变是工程软岩岩体的一个重要特性。流变过程中,岩体参数随时间发生损伤衰减。作为常用的本构模型之一,西原流变模型能较为综合地反映岩体弹-黏弹-黏塑性流变性质,但其元件各参数通常假定为常数,对加速蠕变阶段的岩体状态描述不够准确,与实际出入较大。基于西原模型,结合由损伤导致的参数衰减,提出了考虑损伤的西原流变模型及结合三轴压缩试验确定损伤阈值的方法,并进行了FLAC^3D自定义本构模型二次开发。用本模型对已有文献中的三轴试验进行数值模拟,得到的结果与原试验对比,结果表明：考虑损伤的西原流变模型与实际情况较为吻合;实验加载曲线中体积应变有明显拐点,之后损伤急剧增大,可确定为损伤阈值。     

非贯通节理岩体单轴压缩动态损伤本构模型

非贯通节理岩体是同时含有节理、裂隙等宏观缺陷及微裂隙、微孔洞等细观缺陷的复合损伤地质材料,基于此提出了在非贯通节理岩体动态损伤本构模型中应同时考虑宏、细观缺陷的观点。首先对基于细观动态断裂机理的经典动态损伤本构模型——TCK模型进行了阐述,其次针对目前节理岩体损伤变量定义中仅考虑节理几何参数而未考虑其强度参数的不足,基于能量原理和断裂力学理论推导得出了同时考虑节理几何及强度参数的宏观损伤变量(张量)的计算公式;第三,基于Lemaitre等效应变假设推导了综合考虑宏、细观缺陷的复合损伤变量(张量);第四,借鉴前人基于复合材料力学的观点,考虑了节理法向及切向刚度等变形参数对岩体动态力学特性的影响,进而建立了基于TCK模型的非贯通节理岩体单轴压缩动态损伤本构模型。并利用该模型讨论了载荷应变率、节理内摩擦角、节理厚度、节理法向及切向刚度和节理倾角等对岩体动态力学特性的影响规律。计算结果与目前的理论及试验研究结果比较吻合,从而说明了该模型的合理性。     

层状岩体结构面特征及其拟合描述模型

为了研究层状岩体的压缩破坏特征,并建立相应的理论描述模型,利用快速拉格朗日元法对三轴压缩情况下的层状岩体进行数值计算,分析了层状压缩强度与结构面倾角的关系,建立了相应的理论模型。利用所建模型,对参考国外文献中的试验结果进行描述,并探讨了加载速率和岩体压缩强度的关系。结果表明：所建拟合模型能够较好描述层状岩体的强度特征,得到了数值试验和室内试验的验证;岩体压缩强度随加速速率增大而增大,二者符合线性关系特征。     

泥岩夹层对盐岩变形和破损特征的影响分析

针对我国大多数盐矿的多层盐岩地质构造特征，对含泥岩夹层盐岩、纯泥岩和纯盐岩3种岩芯试样进行单轴压缩和不同围压下三轴压缩试验研究，对比分析3种试样的变形和破坏特性。试验结果表明：泥岩夹层对盐岩体的变形和破坏特性有明显的影响，强度高于盐岩的泥岩夹层却先于盐岩出现横向拉伸破坏；此外还观察到应力-应变曲线的“应力跌落”现象。针对试验结果，利用Cosserat介质扩展理论对泥岩夹层的影响进行理论分析。分析结果表明，泥岩和盐岩力学特性上的不匹配导致二者界面附近泥岩体等效受到横向拉伸应力作用，这很好地解释试验结果，这一分析结果可对进一步进行层状盐岩体内油（气）储库洞室稳定性分析提供理论基础。     

盐岩的压缩试验研究与损伤模型模拟

基于C.S.Desai提出的扰动态概念,对围压0～25MPa范围内盐岩在三轴压缩下的应力-应变关系进行初步模拟分析。在此模拟过程中,材料的完整态采用ZWT模型进行描述,完全扰动态则采用Mohr-Coulomb摩擦力的表述形式。由此研究模型中相关系数与围压的关系。研究结果表明:此模型能够较好地应用于盐岩三轴压缩过程的模拟分析。盐岩的三轴损伤模型的计算结果与试验结果吻合较好。     

盐岩非线性蠕变损伤本构模型及其工程应用

盐岩以其良好的流变、低渗透率及损伤自我恢复等特性,是目前国内外公认的能源、废弃物贮存和高放射废物地质处置的首选介质。结合金坛储气库盐岩三轴蠕变的研究成果,建立盐岩三维蠕变损伤的本构方程和损伤演化方程,并将建立的本构方程编制成有限元计算程序,模拟金坛储气库在注采过程中的蠕变和损伤演化的影响范围。研究成果对金坛储气库的运行压力设计具有一定的参考意义。     

蠕变作用下层状盐岩界面剪切应力变化规律研究

层状盐岩储气库在长期运行期间,盐岩层与泥岩夹层两种不同岩石材料会因较大的蠕变差异引起变形不协调,在交界面产生剪切应力导致界面滑移破损,使油气沿着界面破损通道泄漏。结合我国在层状盐岩地质构造中建设能源储库的实际情况,通过对湖北云应地区层状盐岩中的盐岩和泥岩夹层的三轴蠕变试验表明,盐岩蠕变率比泥岩蠕变率至少大一个数量级。建立蠕变本构模型,对储库长期运行期内层状盐岩界面剪切应力的分布和变化规律进行数值分析,并讨论了盐岩稳态蠕变率不同时,界面层间剪切应力的变化情况。研究结果表明,界面层间剪切应力集中分布在靠近腔体的小段范围内,其后层间剪切应力迅速降低,且随着蠕变时间的增加,层间剪切应力增加的越来越快;盐岩的蠕变特性越好,层间剪切应力越大,当界面剪切应力达到界面抗剪强度时,就会发生破坏。蠕变试验和数值模拟结果,对于层状盐岩渗漏稳定性和密闭性研究具有重要指导意义。     

层状盐岩储气库物理力学特性与极限运行压力

分析研究层状盐岩中盐岩与夹层的孔隙率与渗透率,指出在未遭受变形破坏条件下,盐岩及其夹层本质致密,不会对储气库安全造成影响.但是,对盐岩与夹层的单轴和三轴力学特性、长期蠕变性进行实验研究与理论分析认为,盐岩与夹层间存在的力学特性差异及变形不协调,会造成2种岩性交界处的剪切破坏,影响储库安全性.在此基础上,提出层状盐岩储气库极限运行压力确定原则,包括顶板稳定、蠕变控制、腔体致密及裸井致密等,并结合实例进行分析研究.所得结果对我国层状盐岩中天然气储库的建造与运行具有一定指导意义及应用价值.     

含夹层盐岩双重介质耦合损伤模型研究

针对中国地下油气储库建设中所出现的含夹层盐岩问题,考虑夹层和盐岩层之间存在地质界面,采用以节点位移和孔隙压力为自由度的界面单元来模拟水力损伤造成的地层界面的开裂、扩展和流体渗漏;并基于多孔介质流–固耦合理论,建立含夹层盐岩双重介质耦合损伤模型。该模型克服了等效连续介质模型不能正确反映地层界面的渗流问题,又克服了双重介质模型不能考虑地层界面开裂问题。在此基础上,采用数值模拟技术,研究高压流体在泥岩夹层与盐岩的界面渗透及其开裂扩展特征,结果表明,高压流体沿腔体围岩渗漏过程中,含夹层盐岩界面呈扇形状张开,沿界面通道流体压力逐步降低。因此,在层状盐岩储库运营过程中,要严格控制腔体压力,避免在含夹层盐岩分层界面上产生油气渗漏,保持腔体的致密性及稳定性。     

层状盐岩中裂纹扩展规律的细观实验研究

 层状盐岩力学行为的研究是地下盐岩储库建设的科学基础。借助可以进行原位加载下实时观测的SEM实验系统,对载荷作用下层状盐岩开裂破坏时的裂纹扩展规律进行研究,并从能量耗散及释放的角度探讨层状盐岩破坏的细观机制。观察发现,层状盐岩中岩盐与泥岩夹杂的细观结构有着较大差异,从而影响着盐岩的宏观力学行为。层状盐岩中裂纹扩展路径一般不是沿着岩盐与夹杂的界面,而是在岩盐或夹杂中扩展。且当裂纹在泥岩夹杂中扩展时,能量耗散明显,形成不规则的裂纹,其分形维数大于1,并有可能出现分叉,此时盐岩承载能力较弱,外载做功还较小时就会导致破坏;当裂纹在岩盐中扩展时,伴随较多的弹性能快速释放,形成较平直的裂纹,其分形维数近似为1,并有可能出现较大张开,此时盐岩承载能力较强,吸收外载做功的能量较高。     

层状盐岩溶腔储气库长期运行稳定性研究

针对我国盐岩矿床地质结构的典型特征:单层厚度薄、软弱夹层多,盐岩夹层对储气库的建造及稳定性具有很大影响,提出了层状盐岩溶腔储气库稳定性研究内容,并运用块裂介质固流耦合数值模拟方法对其影响作了进一步分析.主要结论如下:(1)层状盐岩矿床储气库建造中必然会受到高盐份泥岩夹层的影响.夹层的存在使腔体的稳定性减弱,延缓了腔体内流体的输运、对流扩散过程,增加了建腔的时间,盐岩与盐岩夹层由于蠕变率不同而造成损伤,形成裂隙气体渗漏.(2)高盐份泥岩夹层对储气库的稳定性影响很大,由弹塑性数值模拟可知,运行气压为8-24MPa,形状为椭球腔,且长短轴为7/4时结构最稳定.(3)盐岩夹层对储气库的渗透性具有一定影响.蠕变损伤在盐岩与夹层交接面处产生的裂隙是气体渗透的主通道,具有较强的导气能力,但渗透速度随着时间的延长而减小,最后达到稳定状态,且最外缘气体压力始终低于0.2MPa,从考虑气体渗透的角度,在层状盐岩矿床内建造储气库是可行的.所得结论可为我国薄盐岩矿床内建造储气库提供一定的参考依据.     

岩石破裂过程分析(RFPA2D)系统的细观单元本构关系及验证

从岩石的细观非均匀性特点出发，应用弹性损伤力学，对岩石破裂过程分析数值（RFPA^2D)模拟系统的本构关系进行了描述。用RFPA^2D系统模拟了岩石试样在单轴或双轴载荷作用下的强度特征和破坏形态，通过与实验结果对比验证了该本构关系和数值模拟系统的合理性和有效性。