考虑渗透系数变化的地下结构温度-渗流耦合分析

地下岩体的温度随着深度的增加呈逐渐上升的趋势，在地下水的作用下，岩体温度-渗流耦合作用明显。分析了地下岩体温度场与渗流场间的相互作用效应，考虑流体渗流和岩体骨架本身的热传导作用以及温度梯度所产生的渗流作用，并计入温度对水体运动粘滞系数及岩体渗透率的影响，给出了完整的基于连续介质模型的地下岩体温度-渗流耦合数学模型。借助...     

裂隙岩体渗流模型研究现状与展望

裂隙岩体渗流对于边坡、地下工程及基础岩土体的承载能力有显著的制约作用。本文简要地介绍了多种裂隙岩体渗流模型研究现状,评述了几类比较有代表性的渗流模型特点以及存在的不足,为选取合理的数学模型用于求解具体的裂隙岩体渗流问题提供了参考依据,并在上述基础上提出了一些需要进一步研究的问题。指出就目前应用最为广泛的等效连续介质模型而言,裂隙岩体几何参数、有效孔隙度以及等效渗透张量的确定仍有待于更深入的研究,而从细观力学结构入手研究渗流耦合模型将具有重大的理论意义和实用价值。     

不同耦合介质爆破爆炸能量传递效率研究

不耦合装药爆破可有效降低孔壁峰值压力,改善爆破效果。针对不同耦合介质爆破时爆炸能量的传递问题,考虑岩体应变率效应,理论分析了爆炸作用下岩体变形及破坏特征,得到不同耦合介质爆破时理论爆炸能量的传递效率,并结合数值模拟研究了岩体性质、炸药类别及不耦合装药系数对不同耦合介质爆破时爆炸能量传递效率差异的影响。研究结果表明,爆炸能量传递效率与耦合介质相关,装药结构相同时,水耦合爆破比空气耦合爆破爆炸能量传递效率高;不同耦合介质爆破爆炸能量传递效率的差异受爆破岩体、炸药性质及不耦合装药系数等因素影响;装药系数相同耦合介质不同的爆破,岩体强度越高,不同耦合介质爆炸能量传递效率差别越大;岩体性质相同时,不同耦合介质爆破间能量传递效率差异随不耦合装药系数的增大而增大,对于乳化炸药在粉砂岩中的爆破,不耦合装药系数由1.28增至3.44时,水耦合爆破传递入周围岩体的能量由空气耦合爆破的1.45倍增至6.52倍。研究结果可对优化爆破设计、改善爆炸能量分布提供参考。     

岩体结构识别的三角多项式函数模式图方法

通过构造工程岩体结构分类的三角多项式函数,基于实际工程的各种岩体结构的典型特征,建立数学模型,并利用其数学特征所展示的韵律进行岩体结构类型的划分和识别。应用这种方法对长江三峡工程坝基块状岩体及清江高坝洲层状坝基岩体典型岩体结构建立三角模式图,作为各种类型岩体结构的标准模式图。岩体结构识别的三角多项式函数模式图方法为工程岩体结构类型的定量划分提供了一种有效的方法。     

岩体变形参数渐变取值模型及工程应用

分析了当前岩体现场试验、变形参数取值和常规计算模式存在的问题。针对河谷、斜坡岩体从表部到深部 ,岩体风化、岩体波速和岩体地应力等彼此均呈现渐变性特征 ,提出了岩体变形参数渐变取值模型 ,并以黄河上游李家峡水电站初期蓄水监测资料 ,论证了岩体变形参数空间上连续性渐变取值的合理性和科学性 ,对今后其它岩体工程的力学参数取值和计算模型的概化提供了指导意义.     

锚杆拉拔力学作用机制的DDA方法模拟研究

为研究锚杆拉拔实验中锚杆与岩体间的力学作用过程和机制,在非连续变形分析(discontinuous deformation analysis, DDA)方法中实现锚杆-岩体界面耦合作用算法,采用多线性黏结滑移界面模型模拟锚杆的滑移脱黏,并针对锚杆拉拔力作用节点和脱黏节点提出专门算法。通过模拟结果与实验结果的对比,对文献中的三线性和四线性黏结滑移界面模型进行参数修正,并验证所提锚杆-岩体界面耦合作用算法的正确性和适用性。通过基于修正参数的多线性黏结滑移界面模型的锚杆拉拔DDA模拟,研究了锚杆的拉拔力-拉拔位移关系,以及不同拉拔力条件下锚杆-岩体界面切应力和锚杆轴力的分布与变化规律,并探讨了锚杆弹性模量对锚杆拉拔力学响应的影响。锚杆-岩体界面耦合作用算法的实现,也为岩体锚固的DDA模拟分析创造了条件。     

页岩基质解吸-扩散-渗流耦合实验及数学模型

采用川南地区龙马溪组页岩样品,设计了页岩基质解吸-扩散-渗流耦合物理模拟实验,揭示了页岩基质气体流动特征以及压力传播规律.推导了页岩气解吸-扩散-渗流耦合数学模型并且利用有限差分法对数学模型进行数值求解,与实验结果相比较表明该数学模型能够很好地描述气体在页岩基质中的流动规律.同时对页岩基质气体流动的影响因素进行了分析,认为页岩基质的渗透率、扩散系数、解吸附常数等因素均能影响页岩基质气体的流量和压力传播规律,在页岩气藏的开发过程中需要考虑这些参数的影响,该数学模型为页岩气井产能计算提供了更准确的计算方法.      

宝泉抽水蓄能电站岩体抗剪参数的选取

本文讨论了岩体抗剪参数C、f值的选取方法, 推导了C、f值的数字特征公式。通过对宝泉抽水蓄能电站工程岩体抗剪参数的选取, 结果表明：由于岩体抗剪参数具有随机性和模糊性, 采用随机-模糊法选取的C、f值要符合实际；文中给出的数字特征公式可对C、f值的波动性及相关性进行检验。     

非连续变形分析方法的改进及应用1）

首先,采用位移收敛判据对非连续变形的岩石断裂分析进行改进并程序实现. 其次,将非连续变形分析方法和裂隙网络渗流模型结合起来,建立了渗流应力耦合分析模型,研究了裂隙岩体变形对渗流的影响以及渗流与应力耦合作用下裂隙岩体的变形破坏特征. 研究发现,主干裂隙控制着渗流场分布,考虑耦合作用时的水头值比不考虑耦合作用时大,说明考虑耦合后坡体更容易变形失稳. 改进前后裂纹扩展形态基本不变,表明程序计算的收敛判据是有效的.     

非连续变形分析方法的改进及应用

首先,采用位移收敛判据对非连续变形的岩石断裂分析进行改进并程序实现.其次,将非连续变形分析方法和裂隙网络渗流模型结合起来,建立了渗流应力耦合分析模型,研究了裂隙岩体变形对渗流的影响以及渗流与应力耦合作用下裂隙岩体的变形破坏特征.研究发现,主干裂隙控制着渗流场分布,考虑耦合作用时的水头值比不考虑耦合作用时大,说明考虑耦合后坡体更容易变形失稳.改进前后裂纹扩展形态基本不变,表明程序计算的收敛判据是有效的.     

SOLUTION AND ANALYSIS OF CIRCULAR TUNNEL FOR THE STRAIN-SOFTENING ROCK MASSES CONSIDERING THE AXIAL IN SITU STRESS AND SEEPAGE FORCE

为了研究轴向应力和渗透力共同作用下软化围岩的应力与位移的变化及分布规律. 基于摩尔-库伦屈服准则及应力-应变软化模型并考虑轴向应力和渗透力的共同作用,将整个塑性区分为有限个同心圆环,以弹塑性交界面处的应力、应变为初始值,并采用微小径向应力增量逐步求出各个圆环上的应力应变及塑性区半径,据此重构了考虑渗透力和轴向力共同作用下软化围岩应力应变特性的逐步求解方法. 利用该方法,推导出软化围岩应力应变的解. 计算结果表明：在考虑轴向应力作用下,塑性区半径和隧道围岩位移都随着渗透力的增加而有所减小;当轴向应力为最小主应力时,渗透力的影响更为显著. 这说明渗透力的存在对于隧道围岩的应力应变分布以及塑性半径和围岩的位移有不可忽略的影响.     

岩体系统渗流场与应力场耦合的广义双重介质模型的应用研究

本文通过实验资料分析, 提出了裂隙岩体应力与渗流之间的分形几何关系式。以研究区岩体系统结构为基础, 提出了岩体系统应力场与渗流场耦合的广义双重介质模型, 并将该模型应用于某水利水电工程的渗流与控制中, 对坝区岩体进行了应力场和渗流场的耦合分析, 取得了比较满意的结果。     

A VARIATIONAL INEQUALITY APPROACH FOR NON-STEADY SEEPAGE FLOW THROUGH THREE-DIMENSIONAL FRACTURE NETWORK

为了求解裂隙岩体有自由面非稳定渗流问题,将Darcy定律延拓至整个研究区域,使得潜在溢出边界条件满足Signorini型边界条件,建立了三维裂隙网络非稳定渗流问题的抛物型变分不等式（parabolic variational inequality,PVI）提法,并证明其与偏微分方程（partial differential equation,PDE）提法的等价性,从而将自由面上的流量条件以及潜在溢出边界上的互补条件转化成自然边界条件,降低该问题求解难度。同时给出了基于PVI提法的有限元数值求解方法,通过与交叉裂隙模型理论解的对比分析,证明了该方法的正确性。最后将该方法对含复杂三维裂隙网络的边坡进行非稳定渗流分析,计算结果表明该方法对于复杂裂隙网络求解具有较强的可靠性和适应性。     

Nonlinear Coupled Mathematical Model for Solid Deformation and Gas Seepage in Fractured Media

Based on detailed investigation into the interactional physical mechanism of solid deformations and gas seepage in rock matrix and fracture, a nonlinear coupled mathematical model of solid deformation and gas seepage is put forward and the FEM model is built up to carry out numerical analysis. The coupled interaction laws between solid deformations and gas seepage in rock matrix and fractures has been emphasized in the model, which is a vital progress for coupled mathematical model of solid deformation and gas seepage of rock mass media. As an example, the methane extraction in fractured coal seam has been numerically simulated. By analyzing the simulation results, the law of methane migration and exchange in rock matrix and fractures is interpreted.     

岩体中饱和渗流应力耦合模型研究进展

在大坝渗流与控制、各种水工隧洞和交通隧道以及地下各种洞室开挖、采矿工程与油藏开采中的水力劈裂、 岩坡和坝基的稳定性研究等领域中经常会遇到岩体的饱和渗流应力耦合问题.本文对这些领域中的渗流应力耦合问题分为6种力学模型: 等效连续介质模型、 裂隙网络模型、 双重介质模型、 断裂力学模型、 连续损伤力学模型和统计模型.前3种主要从经典的黏弹塑性本构着手考虑, 着重对渗流场的处理, 后3种则是从考虑岩体在耦合作用下发生损伤破裂行为方面切入, 着重于岩体结构内部发生质的改变带来的更加复杂的耦合效应.对这些模型进行了详细的介绍, 指出了这些模型的优缺点与适用范围, 给出了一些有代表性的研究成果, 并对该研究领域中未来的研究方向进行了展望.     

考虑轴向力和渗透力时软化围岩隧道解析

为了研究轴向应力和渗透力共同作用下软化围岩的应力与位移的变化及分布规律. 基于摩尔-库伦屈服准则及应力-应变软化模型并考虑轴向应力和渗透力的共同作用,将整个塑性区分为有限个同心圆环,以弹塑性交界面处的应力、应变为初始值,并采用微小径向应力增量逐步求出各个圆环上的应力应变及塑性区半径,据此重构了考虑渗透力和轴向力共同作用下软化围岩应力应变特性的逐步求解方法. 利用该方法,推导出软化围岩应力应变的解. 计算结果表明:在考虑轴向应力作用下,塑性区半径和隧道围岩位移都随着渗透力的增加而有所减小;当轴向应力为最小主应力时,渗透力的影响更为显著. 这说明渗透力的存在对于隧道围岩的应力应变分布以及塑性半径和围岩的位移有不可忽略的影响.      

露天矿边坡饱和渗透系数随机场数值模型研究

为分析强降雨入渗及渗透系数空间变异性对闭坑露天矿边坡渗流场的影响程度，基于非饱和渗流理论和随机场理论，采用非侵入式随机方法，通过FISH语言编写非饱和区单元饱和度、渗透系数与基质吸力的修正函数，建立饱和-非饱和渗流随机场模型，开展强降雨作用的高大陡深岩质边坡渗流特征研究。研究结果表明，修正饱和-非饱和渗流随机场模型能够准确地描述露天矿边坡的降雨入渗过程。且降雨入渗主要影响到露天矿边坡浅层渗流场。随降雨持时变化，在坡面与地下水位线之间形成包围的且逐渐缩小的非饱和区。坡面最早出现暂态饱和区且厚度逐渐增加，但增幅逐渐放缓。同时坡面点孔隙水压力最早达到稳定值零，离坡面越远的点孔隙水压力达到稳定值零的用时会越长。该结论可为闭坑露天矿边坡的地质灾害风险预测提供参考意义。     

孔隙网络模型在渗流力学研究中的应用

孔隙网络模型广泛地用于同多孔介质有关的微观模拟中. 本文简介了利用网络模型研究储层渗流规律的基本思想、主要步骤以及网络模型同其它微观模型相比较的优缺点, 总结了在微观渗流研究中使用的两大类网络模型(准静态网络模型和动态网络模型)的特征及其适用范围, 综述了目前网络模型在研究渗流中的应用现状以及国内外研究比较活跃的几个方面,最后分析了网络模型今后的发展趋向.      

海底隧道涌水量与覆岩厚度关系研究

对半无限空间水文地质模型，分别采用经验解析法和数值模拟法计算海底隧道涌水量，分析海底隧道涌水量与覆岩厚度的关系. 并根据日本最小涌水量法，对上述计算结果进行插值，求出不同海水深度所对应的最小覆岩厚度. 最后，拟合出海水深度与最小覆岩厚度关系的直线公式.