基于地层信息的三维洞室可视化仿真技术研究

在大规模的岩土工程建设中隧洞和洞室群是最为常见的地下结构。由于受岩土工程所处复杂地质条件的影响,针对地层和洞室的三维GIS综合建模及可视化分析是信息岩土工程面临的难题。讨论了基于栅格模型的地层三维建模技术,提出了由多层数字等高模型（DEM）生成三维地层的自适应控制算法。在此基础上,系统研究了地层中的三维洞室可视化仿真技术,提出了地层信息和洞室的三维综合建模技术,通过地下洞室的三维实体模型生成、洞室和地裙相减、洞室和地层相交等运算,实现了基于地层模型的洞室开挖与漫游,最终生成了赋含地层信息的三维洞室,并给出了一个典型的工程应用实例。     

统一弹塑性本构模型在FLAC3D中的计算格式

统一强度理论是一个有着坚实理论基础、可表达各种岩土材料强度特性且被广泛应用的新的强度理论。FLAC3D是一个具有强大的计算分析功能且专门针对岩土工程问题开发并被广泛应用的数值分析软件。若能将二者结合起来,无疑会大大促进岩土工程领域相关问题的解决。针对这一问题,研究了统一弹塑性本构模型在有限差分方法中应用的基本理论格式。根据FLAC3D软件中UDM接口计算格式的要求,将统一强度理论引入其中,详细推导了统一弹塑性本构模型在FLAC3D软件中应用的计算公式。由于统一屈服面在应力空间内由12个面组成,在应力角为[0,π/3]的范围内,统一屈服面由两个相交面组成,为计算塑性因子及处理计算中应力超出屈服面的应力调整问题,利用应力角分析了应力空间的分区,推导了区域分界面的公式。最后,编制了相应的UDM接口程序对一圆形隧洞进行了弹塑性分析,对比了计算结果与解析解,结果很好地验证了此计算格式及相应接口程序的正确性。统一强度理论和FLAC3D软件的结合,将使二者的优点得以充分发挥,以解决更多的岩土工程问题。     

考虑温度损伤的盐岩蠕变本构关系研究

根据统计力学原理,以分形岩石力学为桥梁,对盐岩在温度与应力耦合作用下蠕变特性进行了研究,导出了考虑围压效应的损伤变量表达式。综合前人的成果,推导出温度–应力耦合下的盐岩损伤方程。对盐岩在温度与三轴应力共同作用下的蠕变试验数据进行研究,分析了盐岩蠕变过程中衰减蠕变阶段、稳态蠕变阶段和加速蠕变阶段的各自特征。对于广义Bingham 蠕变模型在衰减和稳态蠕变阶段引入一非线性函数,在加速蠕变阶段引入损伤,建立了盐岩考虑温度损伤的蠕变本构关系。利用盐岩的三轴蠕变试验数据对新的蠕变损伤模型参数进行辩识与比较分析。结果表明,该模型能够很好地反映在不同温度作用下盐岩的衰减蠕变阶段、稳态蠕变阶段和加速蠕变阶段的规律,验证了该模型的合理性和正确性。     

河流相砂泥岩薄互层地震反射特征研究

A sedimentary geological model is established in order to study the seismic reflection characteristics of channel sand bodies. Synthetic seismic shot gathers are simulated using the acoustic wave equation and then are prestack time migrated. On the imaged data, the reflection characteristics and instantaneous attributes are analyzed and log-constrained impedance inversion is tested. Because of wave field interference, the experimental results show that seismic events do not definitely correspond to the channel sand bodies and that seismic modes of occurrence do not represent the actual ones. The seismic events formed by wave interference may lead to errors and pitfalls in sand body interpretation. The corresponding relations between instantaneous seismic attributes and sedimentary sands are not well established. Log-constrained impedance inversion improves the resolution of channel sands. However, if the inverted resolution is forced to be too high, artifacts related to the initial model may occur.     

锦屏二级水电站引水隧洞长期稳定性数值分析

针对锦屏II水电站深埋引水隧洞高地应力、强渗透压的复杂赋存环境及硬脆性岩体的特点,采用凝聚力弱化-摩擦强化（CWFS）模型反映岩体的脆塑性变形特点。采用黏弹性模型反映其黏性变形规律,基于现场的位移监测资料,通过遗传规划和粒子群算法的流变模型及参数识别方法得到黏弹脆塑性流变模型及其参数,利用大型岩土工程分析软件FLAC3D,模拟了隧洞群岩体的流变力学行为、喷锚支护型式及其与围岩相互作用机制以及工程结构体在高渗透压下的流固耦合效应。通过模拟结果的分析,采用位移、位移速率、应力、破坏接近度及支护结构的受力状况等状态信息的综合集成分析方法,对此工程围岩及结构体在复杂赋存环境中的长期稳定性状态进行了评价,并论证了现有设计方案的合理性,为指导工程设计和施工起到重要的作用。     

岩土力学参数空间变异性的集合卡尔曼滤波估值

岩土参数具有结构性和随机性的空间变异特征,该特征导致岩土参数具有不确定性。以地质统计学作为岩土参数空间变异性分析的理论基础,将分布于研究区的岩土参数视为区域化变量,变异函数既描述了岩土参数整体的空间结构性变化,又描述了其局部的随机性变化,用变异函数理论模型作为描述岩土参数空间变异规律的数学模型。引入集合卡尔曼滤波（EnKF）分析方法,利用时空分布的观测数据,对岩土参数空间变异性进行估值。数值算例表明,EnKF能够有效地融合观测数据,较好地提供岩土参数空间变异性的估值。     

考虑体积应力效应的大理岩硬化软化力学模型

以锦屏二级水电站大理岩室内三轴压缩试验结果为基础,推导了考虑体积应力效应的塑性内变量公式,并研究了大理岩强度参数及剪胀特性随内变量的变化规律,结果表明,大理岩在峰前阶段凝聚力稍有降低,在峰后阶段迅速减小;内摩擦角随内变量增加而增大,并在内变量为0.8～0.9时达到最大值;不同围压下大理岩剪胀角均随内变量增大,且表现为先增大后减小的变化趋势,但低围压和高围压条件下剪胀角变化趋势不同。在此基础上建立了大理岩考虑体积应力效应的力学模型,并通过FLAC3D模拟了大理岩室内常规三轴压缩试验,模拟结果与室内试验结果吻合很好,表明该模型可较好地反映大理岩的主要力学特性。研究方法和成果可为其他深部脆性岩体的变形破坏分析提供重要的参考价值和借鉴意义。     

不需提取激发脉冲的探地雷达波形反演方法

利用麦克斯韦方程的波形反演方法,能充分利用雷达数据中反映地下介质特性的波形、走时、相位等信息.它要事先已知激发脉冲信号,需从实际资料中提取,但受到诸多因素的影响.本文尝试避免从实际雷达资料提取激发脉冲的方法.首先,借助于一个已知地质模型上的探地雷达系统采集的记录以及给定激发脉冲、给定激发接收天线的正演模拟信号,计算出转换算子.其次,将转换算子作用于实测雷达信号,将其转换成给定激发脉冲激发的、给定激发接收天线采集的信号.最后,用转换后的数据进行波形反演.二维模型数据的转换表明,转换后的数据与正演模拟数据一致,反演结果差异甚小.实验室实测数据的反演结果与实际情况相符.     

西昆仑库地煌斑岩的年代学及地球化学特征

在西昆仑库地发现的煌斑岩脉为钙碱性、钾质闪煌岩,侵位在加里东期的变质增生楔杂岩中,采用角闪石单矿物＾４０Ａｒ／＾３９Ａｒ同位素年代学方法确定粕斑岩的侵位年龄为４０４Ｍａ。煌斑岩的稀土元素总量较高,配分图呈轻稀土富集的右倾型式;微量元素以相对富含固离子亲石元素、亏损高场强元素为特征;结合地球化学、同位素特征和大地构造背景分析,煌斑岩来源于受原特提斯俯冲组分释放富含大离子不相容元素流体交代过的亏损地幔