差异风化型危岩力学模型及破坏机制研究

差异风化型危岩多形成于由砂岩和泥岩等软、硬岩性组成的互层陡倾岩质边坡中,其崩塌破坏属地质灾害中的常见形态。以普洱渡危岩高边坡为例,推导出差异风化型危岩岩腔后壁泥岩压应力随岩腔深度增加而增大的计算公式。据普洱渡危岩高边坡地形地貌和地质构造特征建立了两类4种差异风化型危岩力学模型：第1类为主控结构面贯通率等于1时的泥岩基座压裂破坏型和转动破坏型危岩;第2类为主控结构面贯通率小于1时的坠落型和倾倒型危岩。利用岩石强度理论推导出研究区砂岩的二次抛物线型摩尔强度包络方程,以极限平衡理论和摩尔强度理论推导了4种危岩体在自重、地震力和裂隙水压力共同作用下的崩塌破坏判别表达式,并由此反演出危岩体的临界崩塌边界方程,得出危岩体崩塌边界（岩腔深度、危岩体厚度和高度、主控裂隙深度）之间的关系,为现场预测和判别危岩体的稳定性和崩落时序提供简便可靠的依据。     

页岩各向异性特征的试验研究

为研究彭水页岩气区块储层的各向异性特征,开展了石柱县龙马溪组页岩的单轴和三轴压缩试验,分析了其力学特性、强度特征和破裂模式的各向异性,并揭示了其破坏机制的各向异性。结果表明：(1) 龙马溪组页岩具有明显的各向异性特征,弹性模量在平行层理方向最大,垂直层理方向最小,且围压的增加使其增加速率不断减小;0°、30°和60°、90°页岩的泊松比随围压的增加呈现出了相反的变化规律,这可能与页岩层理间孔隙和微裂缝的良好发育有关。(2) 相同围压下,0°试样强度最高,90°次之,30°最低,总体上呈现出两边高、中间低的U型变化规律,而不同角度的Hoek-Brown强度准则能较好地反映其强度的各向异性特征。(3) 页岩破裂模式的各向异性是由破坏机制的各向异性引起的,而强度的各向异性是由破坏机制的各向异性控制的。单轴压缩时,0°页岩为沿层理的张拉劈裂破坏,30°为沿层理的剪切滑移破坏,60°为贯穿层理和沿层理的复合剪切破坏,90°为贯穿层理的张拉破坏。三轴压缩时,0°为贯穿层理的共轭剪切破坏,30°为沿层理的剪切滑移破坏,60°和90°为贯穿层理的剪切破坏;页岩地层的层状沉积结构和层理间的弱胶结作用是破坏机制各向异性的根源。研究结果为水平井井壁的稳定性分析和水力压裂施工设计等提供了技术参考。     

基于多重分形特征的岩体结构面剪切强度研究

Barton剪切强度模型是目前工程实践中普遍采用的强度公式,而实践应用时,该模型中结构面粗糙性系数（JRC）的评估存在主观性和片面性的缺点。鉴于此,应用多重分形理论,提出了采用多重分形参数准确量化JRC的方法。首先应用数字图像处理技术获取了结构面三维形貌数据信息;然后采用投影覆盖法进行了结构面的分形维数计算,重点对结构面的多重分形特征进行了分析,通过对比分析了构造的15个结构面的多重分形特征值。结果表明：结构面越粗糙,多重分形特征参数 和 值越大, 或 能够很好地描述结构面的形貌特征。最后对9组石膏试件进行了剪切试验,通过不同正应力下对应的剪应力数据分析,结合结构面形貌多重分形参数 或 ,以Barton剪切强度公式为基础,应用最小二乘法原理,给出了JRC与 及JRC与 关系。这样在工程实际中,可以用参数 或 对JRC进行估算,克服了JRC人为估值主观性及采用二维标准轮廓线评估片面性的缺点,为准确评估结构面粗糙度提供了一条新的途径。在此基础上,应用Barton模型可准确估算岩体结构面的剪切强度。     

热-水-应力-迁移耦合条件下双重孔隙-裂隙介质的抗剪强度及有限元分析

考虑裂隙的连通率、间距、孔隙基质和裂隙材料在表征单元（REV）中的体积分数,并假定双重孔隙-裂隙介质的等效内摩擦角保持常数,而等效黏聚力是固有黏聚力、等效塑性应变、基质吸力、溶质浓度及温度的函数,提出了一种在热-水-应力-迁移耦合条件下确定表征单元内任一平面上等效的黏聚力及内摩擦角的方法。针对一个假定的位于非饱和双重孔隙-裂隙岩体中的高放废物地质处置模型进行了数值模拟及分析。结果表明,基质吸力对于等效黏聚力的增强作用大于等效塑性应变和溶质浓度的减弱作用,使得等效黏聚力得到了提高,故减少了围岩中的塑性区;由此岩体应力、孔（裂）隙水压力及流速、孔（裂）隙溶质浓度的分布及量值也发生相应的改变。     

弹性纵波垂直通过非线性结构面的傅里叶级数分析法

利用傅里叶级数计算任意函数形式的弹性纵波垂直通过非线性结构面的透射波和反射波的速度波形。采用双曲模型描述结构面的变形特性,基于位移不连续模型,结合波振面处动量守恒定律,推导了弹性纵波垂直通过非线性结构面的基本方程。假设应力波在含结构面岩体中传播时,结构面的存在不改变应力波波形函数的最小正周期,运用傅里叶级数理论和周期延拓方法,得到了任意函数形式的弹性纵波垂直入射时透射波和反射波速度波形的傅里叶级数解,并验证了傅里叶级数解是合理的。利用傅里叶级数解,分析了单一频率的正弦谐波入射至结构面时,透射波中各阶谐波的振幅和相位与谐波阶数的关系。研究结果表明,各阶谐波的振幅与阶数呈负指数关系衰减,前7阶谐波振幅的衰减指数为谐波阶数的二次函数,当谐波阶数大于7时,衰减指数为谐波阶数的一次函数;各阶谐波的相位与谐波阶数呈线性关系。     

基于MPGA的复杂应力状态边坡稳定性分析

极限平衡法在边坡稳定性分析中处于主导地位,该算法通常建立在一定的假设之上,没有考虑边坡土体的非线性应力-应变关系,不能用于未达到极限状态以及加固后边坡的稳定性分析。为了解决以上问题,将有限元计算与多种群遗传算法（MPGA）相结合,建立一种基于MPGA的复杂应力状态边坡稳定性分析通用模型,通过数值应力场求解安全系数,并为多种群遗传算法构建适应度函数;再利用多种群遗传算法为安全系数的计算提供滑移面。为了保证分析的高效、合理,根据滑移面发展趋势,动态产生初始滑移面,并增加一个滑移面约束条件。最后,通过均质边坡和软弱夹层边坡两个典型算例分析验证了该方法的合理性;通过分析土钉加固的软弱夹层边坡,证明了该方法可用于加固等复杂应力状态边坡的稳定性分析。     

南北地震带壳幔结构成像研究

<正>对地壳和上地幔的速度结构进行成像是研究地球深部结构及动力学问题的一个有效途径。南北地震带地区是青藏高原与相邻块体的过渡带,这里的地质构造十分复杂,地震活动频繁。在南北地震带地区进行壳幔结构成像研究,有利于提高对青藏高原及其与相邻块体相     

底板导水裂隙带中含承压水裂隙的力学特性

华北煤田煤系地层的底部为含水丰富的奥陶系灰岩,该层灰岩与地表水系和煤系地层含水层存在着密切的水力联系。许多大型突水是采掘空间与含水层之间形成突水通道引起的。因此,分析开采工过程中底板岩体中含承压水裂隙的断裂力学特性具有重要意义。本文采用作者提出的数值方法分析了含高压水裂隙在开采过程中的断裂力学特性。假设裂隙面上作用着均匀水压力,考虑地应力的影响,采用断裂力学的叠加原理,分析了采煤工作面推进过程中含承压水裂隙的断裂力学特性,讨论了开采过程中底板承压水导升。结果表明:在工作面推进过程中,位于底板岩层且与含承压水岩层连通的裂隙,会在承压水水压力和扰动应力的共同作用下,产生破坏,导水裂隙带的高度会增加;这就增加了底板突水的危险性。     

2008-2015年北京高速公路道面结冰特征分析

道面结冰是北京地区高速公路交通安全的主要危害因素。利用北京地区2008-2015年28个ROSA交通气象站的结冰数据分析了北京地区高速公路道面结冰的基本特征。结果表明：（1）北京市高速公路道面结冰空间差异显著,不同高速公路结冰特征显著不同,同一高速公路不同路段结冰特征也明显不同,这与高速公路各路段局地气候存在差异有关。（2）道面结冰主要发生在北部和东部高速公路所在路段,以小汤山西桥站、顾家庄桥站、六道口桥站和丁各庄桥站道面结冰灾害最为严重。冬季道面结冰月变化显著,11月和3月,结冰次数较少,结冰持续时间短;12月、1月和2月结冰次数多,结冰持续时间长。（3）持续时间越长的结冰过程结冰次数越少。不同持续时间的结冰过程在各时段的结冰次数和累积结冰时长均存在明显差异,两者在各时刻的演变特征基本相同,峰值均出现在22：00,且在00：00之后呈明显减少趋势。（4）大部分站点慢车道比快车道更容易结冰,少部分站点快、慢车道结冰率近似,有的站点甚至快车道结冰率远高于慢车道。快、慢车道水或覆盖物厚度与交通站点所处位置的局地气象条件密切相关。（5）气温、道面温度和大地温度与水或覆盖物的厚度变化呈反相位,即温度越低,水或覆盖物的厚度越大。这可为北京地区道面结冰预报预警方法的开展提供着眼点和依据。     

陇中黄土高原LAS在非均匀下垫面的观测优势及其在区域陆面模式验证中的初步应用

目前通用的通量观测技术涡动相关仪（EC）在区域陆面模式验证中存在能量不闭合和空间代表性有限的问题,寻求改进EC观测热通量的新技术是提高陆面模式验证效果的关键环节,大孔径闪烁仪（LAS）的出现有效改善了这一现状。基于黄土高原定西站2010年1月和6月的同步综合观测资料以及目前比较有代表性的陆面过程模式CLM的模拟数据,分析研究了LAS对EC观测地表能量不平衡问题的改进以及LAS对EC在区域陆面模式验证中的提高,结果表明：LAS可有效解决EC观测存在的地表能量不平衡问题,提高EC的地表能量闭合度,在非均匀下垫面LAS观测优势突出;利用LAS观测的感热通量进行区域陆面模式的验证,能够很大程度地避免EC能量不闭合和空间尺度不匹配在验证中造成的偏差,LAS观测更适合于大尺度模拟的验证,验证效果更好。