地下采矿诱发斜坡移动变形分析

地下采矿与滑坡的关系问题是近期社会关注的热点问题。基于对此问题已有研究的不足,以及随机介质理论在矿区地表移动与变形预测中的广泛应用,探索采用随机介质理论方法对地下采矿条件下坡体的移动变形规律进行研究。以重庆武隆鸡尾山崩滑体为研究对象,建立了坡体内两个地下采矿剖面的计算模型,并通过计算机编程辅助计算,对地下采矿条件下坡体的移动变形规律进行了分析,并对不同地下采矿工况分别进行了详细研究。计算表明,采矿诱发的坡体不均匀沉降加剧了鸡尾山坡体的移动下滑趋势,认为随机介质理论适用于地下采矿诱发的坡体移动变形分析。     

岩溶区多重介质水流模型研究进展

数值模拟技术广泛应用于地下水资源评价工作中,为水资源合理开发利用提供科学依据。本文针对岩溶区的数值模拟技术展开研究,综述了岩溶区地下水流建模技术的主要方法和研究进展,介绍了SWMM模型、UGRFLOW模型、CAVE模型、CFP模型等多重介质模型的建模方式和运行原理。岩溶区地下水赋存介质复杂,管道、裂隙和空隙多重介质并存,水流特征多样,等效连续介质模型不能很好地刻画岩溶区复杂的地下水流特征。随着对岩溶水流系统研究的深入,二重和三重介质模型逐渐应用于岩溶区地下水流模拟工作中,取得较好成果,但仍存在如忽略中宽裂隙的导水作用、岩溶管道设置单一、模型适用性不强等问题。今后应从解析岩溶含水介质的结构特征入手,关注中宽裂隙在地下水流系统中的导水作用和水流特征,改进管道流模块以更好地刻画复杂的地下河管道,改进不同介质间水流交换量算法,提高模型精确性和适用性。      

裂隙(EDA)介质P波地震记录数值模拟及波场特征分析

通过对潜山面和内幕EDA介质地面地震记录数值模拟及其波场特征分析,在了解地震波场的动力学特征对裂缝介质的敏感性基础上,确定主要以振幅,其次是频率的地震波动力学特征为裂缝检测被检参数,以均方振幅与中心频率作为检验检测结果正确与否的评价参数,为建立裂缝检测方法奠定了理论基础.     

非均匀介质中交错网格高阶有限差分数值模拟

地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。     

不连续面在双重介质热-水-力三维耦合分析中的有限元数值实现

在建立双重介质热-水-力耦合微分控制方程的基础上,提出了裂隙岩体热-水-力耦合的三维力学模型,对不同介质分别建立以节点位移、水压力和温度为求解量的三维有限元格式,开发了双重介质热-水-力耦合分析的的三维有限元计算程序,在有限元数值分析中不连续面应力计算采用等厚度空间8节点节理单元进行离散,而不连续面渗流和热能计算时采用平面4节点等参单元进行离散,这样保证了不同介质之间的水量、热量交换和两类模型接触处节点水头、温度和位移相等。通过高温岩体地热开发算例,揭示了在热-水-力耦合作用下不连续面处于低应力区,其张开度随运行时间的延长呈非线性增加,非稳定渗流阶段不连续面显著地控制着渗流场的整体分布,它的水头远高于拟连续岩体介质的水头,而进入稳定渗流阶段不连续面的控渗作用不明显,由于高温岩体地热开发系统中存在大规模的热量补给,不连续面对岩体温度场分布的影响并不显著。     

湘黔地区下寒武统黑色岩系中镍-钼矿床黄铁矿的成因

重点研究了贵州遵义中南村和湖南张家界三岔镍-钼多金属矿床的黄铁矿,矿床成因为热水喷流沉积型.黄铁矿样品S/Fe比值的平均值为2.0093-2.048,成分特征均属于铁亏损型,说明其形成温度低.Co/Ni变化范围为0.077~4.5,均值为2.197.据Co/Ni的比值可以判断本区的黄铁矿主要为热液成因,矿化物质来源为热水喷流体系内的热液.成分图解标型显示,黑色岩系中高ω(Co+Ni)/ω(Fe)的黄铁矿和高ω(As+Se+Te)/ω(S)的黄铁矿均大量出现,可能与热水喷流沉积的温差变化范围较大有关.而并非正常沉积的产物.S/Se比值为950.8-1059.9,说明生成环境均一程度较高,而且海水温度比较高.由中南村镍-钼矿床的Se/Te值也能判明其为热液成因.关于黄铁矿的研究对本区镍-钼矿床的寻找有重要意义.     

通风块石层非达西效应的试验研究

块石通风路基是青藏铁路建设中保护多年冻土的主要措施,准确地描述空气在块石中的流动是块碎石路基冷却效应数值分析的首要问题.通过风洞试验研究了块石层内部渗流速度和压力梯度的关系,试验中块石的平均粒径为12.5cm,孔隙率为0.487,风速和压力差数据是在不同的风机转速下测得.结果表明:块石层内空气渗流速度随压力梯度增大而增大,压力梯度与渗流速度之间存在非线性关系.基于Forchheimer方程,利用试验数据获得了块石层的渗透率和惯性阻力系数.这些结果可为各种块石结构路基降温效果的数值分析提供理论基础.     

非均质介质中地下水流动与溶质运移模拟——问题与挑战

对大空间尺度和长时间跨度的地下水流动及污染物质运移进行预测的需求, 使水文地质研究面临异乎寻常的挑战.这些需求来自于对核废料地质储放方法的安全性评价、地下水污染状况评价及其治理方案的选择.流动系统的非均质性是地下水流动及物质运移模拟中最主要的困难之一, 这种困难来自对非均质系统进行特征描述(通过原位观测实现)、概念化及模拟.评述了非均质介质中流动运移模拟的一些重要问题与挑战, 讨论了解决的途径.讨论的主题包括: 动力流动的沟道化, 示踪剂穿透曲线, 裂隙岩石中流体流动的多尺度, 观测的不同尺度, 模拟、预测与非均质性以及系统特征描述和预测性模拟的分析.      

岩土介质三维快速拉格朗日数值分析方法研究

研究了FLAC3D（fast lagrangian analysis of continua in 3 dimensions）的特点,并与有限单元法作了比较。FLAC3D方法以结点运动方程为支配方程,追踪了介质从受荷到达到平衡状态的过程,而有限元法是根据介质力学平衡方程直接求解,这是二者主要区别。FLAC3D没有采用介质真实的阻尼特性和结点质量,给出的不是介质所经历的真实过程,不能正确反映过程的影响,因此给出的介质应力和变形计算结果的物理意义是不甚明确的。求解过程中的介质振动,是一种噪音,可引起弹塑性介质计算结果误差,而弹性介质的计算结果几乎不受影响。研究指出了FLAC3D方法的优缺点。还通过算例作了四个方面的研究：(1) Drucker-Prager屈服准则与Mohr-Coulomb屈服准则比较;(2) 膨胀角取值对计算结果的影响;(3) 大变形与小变形对计算结果的影响;(4) 精度设置对计算结果的影响。研究表明,Drucker-Prager准则与Mohr-Coulomb准则结果差异颇大;膨胀角取值对结果的影响是敏感和显著的;一般情况下,取小变形模式是合适的,计算精度取10-5是足够的。     

岩体结构与岩体水力耦合计算模型

水力耦合作用是岩体中力学过程与渗流过程相互作用的物理过程。水力耦合机理的理解是水力耦合分析的关键问题,其耦合机理是由岩体结构特性决定的。在分析岩体水力耦合过程基础上,根据岩体的基本结构及代表性单元体（REV）是否存在提出了建立水力耦合模型的方法。当裂隙岩体中不存在代表性单元体（REV）时,提出了裂隙岩体多重介质流固耦合分析的全耦合数学模型,给出流固耦合模型数值方法求解的数学模型及有限元计算表达式。程序编制和验证工作正在进行中。