渗透系数对渗流速度影响的热-水-应力耦合二维有限元分析

使用建立的饱和-非饱和介质中热-水-应力耦合弹塑性模型,取一组真实的缓冲材料、岩体及贮存容器（用加热器代替）的物理力学指标作为基本计算条件,以一个假想的核废料处置库为模拟对象,就高、中、低3种渗透系数进行了二维有限元数值分析,考察了近场的渗流速度大小及方向的变化,认识到：在其他计算参数相同的条件下,缓冲材料的渗透系数越高,近场的渗流速度越大,缓冲层达到饱和状态也越快;当缓冲层达到饱和状态后,如果玻璃固化体继续放热,有可能使得近场的温度梯度的作用强于水力梯度的作用,导致渗流速度方向发生逆转,即水分由内向外迁移;依渗透系数不同,其渗流速度的大小、方向有显著差异,这也势必引起放射性核素随地下水流动扩散的复杂多变。     

CO2注入岩体的热-气-应力耦合二维弹塑性有限元分析

目前国内关于CO2地质埋存课题的研究已经起步,其中与岩石力学有关的工作急需加强。为此,考虑到向含气地层中注入CO2时岩石、气体和温度的相互作用,将CO2视为理想气体,使用Drucker-Prager屈服准则和‘无拉应力’判据,建立了一个热-气-应力耦合模型并研制了相应的二维有限元程序。假定了一个由下部注入层和上覆冠石层组成的CO2埋存地质系统,以此为数值模拟的对象,分析了CO2在不同的注入速率和注入时间条件下岩体的中的位移、应力、受拉与塑性破坏区的变化和分布情况,结果显示,为保证CO2注入岩体的稳定,应优择最佳的注入速率和注入时间。     

压力溶解对颗粒聚集岩体中热-水-应力耦合作用的弹塑性有限元分析

将Taron等提出的颗粒聚集体的压力溶解模型引入笔者所研制的孔隙介质热-水-应力耦合有限元程序中,并使用摩尔-库仑准则,针对一个假设的实验室尺度且位于饱和石英颗粒聚集岩体中的高放废物地质处置库模型,拟定弹性分析和弹塑性分析两种计算工况,进行4 a处置时段的数值模拟,考察了岩体中的温度、颗粒界面水膜及孔隙中的溶质浓度、迁移和沉淀质量、孔隙率及渗透系数、孔隙水压力、地下水流速和应力及塑性区的变化、分布情况。结果主要显示：弹塑性分析中由于应力调整和增大了分子扩散系数,使得塑性区的颗粒介质的溶解、迁移和沉淀有明显的变化,并对渗流场（孔隙水压力及流速）和应力场产生显著的影响。但两种工况弹性区中的颗粒介质的溶解、迁移和沉淀差别较小。     

亥废料处置概念库近场热-水-应力耦合模型及数值分析

特别考虑了缓冲层中温度梯度水分扩散、水蒸汽扩散对水连续性及能量守恒的影响，进一步完善了所建立的分析饱和-非饱和介质中热-水-应力耦合现象的控制方程，并针对试验资料，使用所开发的有限元程序对一个核废料处置概念库近场的热-水-应力耦合过程进行了数值模拟，考察了缓冲层及岩体中若干点的温度、饱和度、孔隙水压力及主应力的变化、分布情况，并将部分结果与国外类似软件的计算数据作了对比，结果显示，在定性和定量上二者的规律有某种一致性，从而得出了一定的认识。     

热-水-应力耦合影响的有限元分析

为了考虑应力拉压和压力（化学）溶解对裂隙开度的综合影响,对所建立的双重孔隙-裂隙介质热-水-应力耦合模型中裂隙开度的计算模型作了改进。通过一个假定的高放废物地质处置库算例,就岩体裂隙开度变化的3种工况,分析了岩体中的温度、孔（裂）隙水压力、地下水流速和主应力的变化、分布情况。结果显示：3种工况的计算域中温度场基本相同;孔（裂）隙水渗流场形态相似,但其量值有一定差别;工况1的裂隙开度在应力和压力（化学）溶解的共同作用下闭合量最大,负孔（裂）隙水压力增值最高;核废物的释热效应明显地改变了岩体自重应力场的水平分量,但对其垂直分量影响较小。     

热-水-应力耦合条件下三项因素对核素迁移影响的有限元分析

目前在放射性核素随地下水迁移的研究中还很少涉及温度场和应力场的影响。针对这一问题,引入渗透迁移方程,将所开发的饱和-非饱和孔隙介质中热-水-应力耦合弹塑性模型及其二维有限元程序进行了扩展和改进,使之可以同步地对温度场、渗流场、应力场和放射性核素浓度场的变化进行解析。对一个假想的核废料处置库算例,通过分别改变影响地下水中核素迁移的3个因素：渗透系数、分配系数和分子扩散系数,考察了在热-水-应力-迁移耦合作用条件下近场放射性核素浓度的分布及变化,从而认识到：当缓冲层的渗透系数小到一定程度后,缓冲层中核素迁移基本不再受渗流场的影响;缓冲层分配系数越小,分子扩散系数越大,核素迁移越快。     

非饱和地质体中水-应力耦合二维有限元方法及对锚杆支护的分析

从建立应力平衡方程、水连续性方程入手,采用Galerkin方法将各控制方程分别在空间域和时间域进行离散,开发出了一个用于分析饱和-非饱和岩土介质中水-应力耦合弹塑性问题的二维有限元程序。然后就一个假定的处于渗流场中的饱和-非饱和土体使用预应力锚杆的支护效果问题进行了模拟计算。结果表明,在有地下水赋存的饱和-非饱和条件下,预应力锚杆的作用主要体现在减少了土体中的塑性区,而对位移的约束依部位而异,其对渗流场的影响不大。     

核废料处置概念库近场热-水-应力耦合模型及数值分析

特别考虑了缓冲层中温度梯度水分扩散、水蒸汽扩散对水连续性及能量守恒的影响,进一步完善了所建立的分析饱和-非饱和介质中热-水-应力耦合现象的控制方程,并针对试验资料,使用所开发的有限元程序对一个核废料处置概念库近场的热-水-应力耦合过程进行了数值模拟,考察了缓冲层及岩体中若干点的温度、饱和度、孔隙水压力及主应力的变化、分布情况,并将部分结果与国外类似软件的计算数据作了对比,结果显示,在定性和定量上二者的规律有某种一致性,从而得出了一定的认识。     

核废料处置概念库开挖及封闭后近场热-水-应力耦合数值分析

考虑了缓冲层中温度梯度水分扩散、水蒸汽扩散对水连续性及能量守恒的影响,进一步完善了饱和-非饱和介质中热-水-应力耦合现象的控制方程。应用缓冲材料的实验资料,对核废料处置概念库开挖及封闭后近场的热-水-应力耦合过程进行了数值模拟,以考察围岩中的位移、塑性区、主应力、孔隙水压力和温度的分布及变化。认为开挖是引起洞室围岩中的变形及应力重新分布的主要原因,但坑道封闭后的热-水-应力耦合过程对变形及应力状态也产生了明显的影响。     

裂隙刚度随应力变化对双重孔隙介质热-水-应力耦合影响的有限元分析

引入并修正了变刚度的连续屈服节理模型,同时考虑应力拉压和压力（化学）溶解对裂隙开度的综合影响,对所建立的双重孔隙-裂隙介质热-水-应力耦合有限元计算程序作了改进。通过一个假定的高放废物地质处置库的数值模拟,就岩体裂隙刚度变化的2种工况,分析了岩体中的温度、裂隙刚度、正应力、孔（裂）隙水压力和地下水流速的变化、分布情况。结果显示：与裂隙刚度是常数时相比,裂隙刚度是法向应力的函数时计算域中温度较低;岩体应力的大小也有一定不同,其分布与裂隙刚度“场”有明显的相似性;并且负孔（裂）隙水压力的绝对值要略小一点,约是常数时的98%。     

铁山隧道和大垭口隧道的平面弹塑性有限元分析

对巴彭公路铁山隧道和开万公路大垭口隧道的两个典型剖面的稳定性进行了平面弹塑性有限元计算,分析了围岩中的应力、变形、塑性区以及初期支护和二次衬砌结构的承载能力,得出了相应的认识：（1）在特定的埋深、岩体类别、支护结构及荷载释放比例的条件下,洞室周围有较强的应力集中现象,一定部位及范围内的围岩进入了塑性状态;（2）喷层、二次衬砌及锚杆均在弹性范围内工作,并且有较大的安全储备。     

龙滩水电站航道座滑边坡平面弹塑性有限元分析

基于弹塑性有限元（FEM）模型和Mohr-Coulomb屈服准则,采用平面有限元法对龙滩水电站航道1+016～ 1+080段座滑后边坡在开挖处理过程中的稳定性进行了数值模拟研究,分析了开挖中岩坡的应力场、位移场、塑性屈服区及支护效应。模拟结果表明,边坡在开挖过程中存在连续变形,局部出现塑性屈服,地下水位升降对边坡影响较大,但边坡整体仍处于稳定状态,开挖过程中有效的支护措施对改善边坡的稳定性至关重要。     

考虑饱和-非饱和渗流、温度和应力耦合的三维有限元程序研制

在探索多场耦合理论基础上,研制出了一个考虑饱和-非饱和渗流场和温度场对应力场作用的三维弹塑性有限元程序,该程序采用了以下简化假设：不考虑力学过程和渗流过程对温度场产生的影响,而只考虑变温对力学过程产生的影响;不考虑力学过程和温度场对渗流场产生的影响,而只考虑渗流场对力学过程产生的影响。经与弹性力学解析解对比,验证了该程序计算结果是合理的。     

黏土岩温度-渗流-应力耦合特性试验与本构模型研究进展

高放废物处置库、垃圾填埋场等工程中常常涉及到温度场（T）、渗流场（H）和应力场（M）的耦合作用的问题。从试验和理论模型两个角度综述国内外黏土岩温度-渗流-应力耦合特性的研究进展,主要包括其传热特性、温度影响下的渗流特性、变形、强度、蠕变特性。在此基础上,重点分析了黏土岩水-热迁移模型以及热-力耦合本构模型的适应性。基于上述认识,通过试验研究了比利时Boom clay在温度作用下的强度、渗透性、蠕变性等特征。结果表明：随着温度升高,Boom clay的强度有所降低,渗透性显著增强,蠕变速率明显加快。提出了适用于Boom clay的THM耦合弹塑性损伤模型,计算结果验证了模型能合理反映温度的影响。最后,探讨了黏土岩THM耦合机理研究的不足和今后的研究方向。