海南保亭县甘什岭滑坡成因分析及稳定性评价

简要介绍了保亭县甘什岭滑坡的形态、地层岩性、水文地质环境条件及地质构造等情况,分析了该滑坡的形成原因,并对其进行稳定性评价和危险性分析,稳定计算方法选用剩余推力传递法对滑坡稳定性加以试算,其计算结果与滑坡实际情况基本吻合。通过以上分析,可得以下结论:当滑坡在天然状态下时,稳定系数K的范围在1.17~1.26之间,此时滑坡处于稳定状态;当滑坡在饱水状态下时,稳定系数K的范围在0.99~1.03之间,此时滑坡处于滑移状态。     

基于CT扫描的致密砂岩渗流特征及应力敏感性研究

认识低渗透储层的渗流特征对油气开采和储层改造具有重要意义。为此,研究利用微CT扫描技术对致密砂岩岩样进行扫描,据此建立了能够精细刻画岩样的微观模型,运用COMSOL模拟了流体在岩石孔隙中的渗流特征,研究了致密砂岩的渗流特征及应力敏感性。研究结果表明:流体入口和出口间压差固定时,岩石的渗透率保持不变,与入口出口压力的具体数值无关;不同方向的岩石模型计算渗透率处在同一数量级但有微小差异;在侧向压力作用下,渗流路径变窄,通过渗流路径的整体速度下降,渗透率下降,但在孔隙相对较大的地方,由于路径变窄,流体速度较未加压力前略有上升。      

应用远震有限频率层析成像反演首都圈上地幔速度结构

本文选用首都圈数字地震台网2003年9月~2005年12月记录的300多个远震事件的波形资料,采用分频带多道互相关方法得到三个不同频段的P波相对走时数据共18499个,计算了每个频段的走时灵敏度核,应用有限频率层析成像反演得到首都圈地区的上地幔三维P波速度结构模型.利用检测板估计了反演结果的分辨率,并与射线层析成像方法的结果进行了比较,说明了反演结果的可靠性.研究结果表明,各构造单元具有明显不同的速度结构特征,其差异可到150 km深:燕山隆起区表现高速;太行山隆起区整体以低速为主并存在小范围高速块体;华北盆地、渤海湾下浅层上地幔中存在大范围的强低速异常,其顶面在50~70 km,可视为软流圈顶面的埋深,这一结果说明华北盆地、渤海湾下岩石圈明显减薄;张家口-蓬莱断裂带是上地幔浅部速度结构的变异带,也是岩石圈减薄的边界带,区内大部分强震都发生在该构造带上,由此看来该带上强震的发生不仅与地壳结构的不均匀性有关,还可能有较深的构造背景.     

基于刚体平动运动单元的上限有限元研究

刚性块体极限分析上限法常应用于岩土工程稳定性研究,然而应用时需假定刚性块体破坏模式并递推繁琐的几何关系。为此,提出一种适应性更广的基于非线性规划模型的刚体平动运动单元上限有限元法,并解决了其优化模型初始值的确定问题。通过引入有限单元思想,将计算区域离散成刚体单元,同时以单元速度和节点坐标作为决策变量,由上限定理建立非线性规划模型获得上限解。利用编制的上限有限元程序进行边坡和浅埋隧道稳定性算例验证,表明运动单元上限有限元法能调整速度间断线至较优方位,所得破坏模式特征鲜明,上限解精度高,可广泛应用于边坡、隧道等稳定性分析研究。     

用变分坡恩迭代方法重建二维非均匀介质结构

提出了用于二维轴对称非均匀介质结构的反演和成像的一种新的反演迭代方法－变分坡恩迭代方法（ＶＢＩＭ）,首先利用玻恩近似将非线民生积分方程线性化,然后应用变分方法导出用于反演的电场积分方程,正演数据则利用高效的数值模式匹配方法获得,数值结果表明,ＶＢＩＭ与ＢＩＭ相比,其收敛速度,成像质量等均得到较大的改善。     

海南省澄迈县福山镇高铁站前道路土体滑坡的GEO5数值模拟研究

通过对海南省澄迈县福山高铁站前道路土体滑坡的形态特征、地层岩性、水文地质环境条件及地质构造等情况的分析,结合地形测量、钻探揭示的地层结构及分类,运用GEO5模拟软件中的地层模块建立了滑坡区三维可视化模型,并在此模型基础上,采用边坡稳定性分析模块,运用萨玛法(Sarma),对该滑体边坡在天然状态和饱水状态2种工况下进行数值模拟研究,得出边坡稳定性情况和滑动面位置。模拟结果与野外调查结果相吻合,为今后制定合理的滑坡防治措施提供了理论依据。     

时域瞬变电磁法三维有限差分正演技术研究

瞬变电磁法应用广泛,三维数值模拟是研究复杂地质模型异常响应规律的重要技术手段之一,也是反演的基础.目前瞬变电磁数值模拟的不足主要有两个方面:第一,场源是在地表水平、浅层介质均匀的条件下计算的,限制了应用范围;第二,地下边界采用Dirichlet边界条件,导致计算空间很大,耗时较长.针对上述问题,在三维正演时,场源采用有限长细导线模型,在Maxwell有源差分方程中直接加入电流密度进行计算.在地表面加入空气层,避免了复杂的向上延拓计算,也可以对地形影响下的响应规律进行分析.在空气边界和地下边界均采用CPML吸收边界条件,并改进了CPML的参数分布,能够吸收空气介质和大地介质中的低频电磁波而反射误差极小,在满足计算精度的条件下可以有效减小节点数量.对循环迭代方法进行优化,将计算域、CPML区域和场源的空间循环统一转化为矩阵方式,加快了计算速度,但是空间消耗增大了约4~5倍.采用三维有限差分正演算法对均匀半空间模型、层状模型和地形模型进行了计算,并与解析解进行了对比验证.     

应用远震有限频率层析成像反演中国东北地区上地幔P波三维速度结构

选取黑龙江、 吉林、 辽宁、 内蒙古区域地震台网, 以及NECESSArray流动台阵记录的223个远震事件的波形资料, 采用多道互相关方法得到了22569个P波相对走时数据, 并计算了相应的走时灵敏度核, 应用有限频率层析成像反演得到中国东北地区上地幔600 km以上的P波三维速度结构模型, 利用检测板评估了反演结果的分辨率。 结果表明, 松辽盆地下方80~200 km的深度上呈主体的低速异常, 与这一地区上地幔浅部的高地温值和低密度的特征相互对应, 可能暗示了部分熔融的地幔。 南北重力梯度带两侧的速度结构明显不同, 这一差异可以延伸到200 km以下, 表明在中国东北地区南北重力梯度带有可能是一条上地幔内部结构的变化带, 或是深部结构的分界线。 长白山火山区下呈大范围的低速异常, 并可从上地幔浅部延伸到地幔转换带中, 推测此低速异常可能反映了地幔转换带内上涌的热物质, 上涌的原因则主要是受到太平洋板块俯冲运动的作用。     

北天山中段微震活动检测及其活动性和孕震构造分析

北天山地震带地处中国大陆强震高发区, 孕震构造复杂, 近年来陆续发生了2016年呼图壁M_S6.2地震和2017年精河M_S6.6地震。 由于测震台网相对比较稀疏, 该区域微震监测能力较弱。 本文主要采用波形模板扫描法对北天山中段(43.5°N~44.5°N; 85°E~87.5°E)进行微震事件检测, 并反演精细的一维速度结构, 重新定位地震; 深入分析该区域的地震活动性和孕震构造特征。 经过微震检测, 得到该地区2014年1月至2018年9月期间57902个地震事件, 是原地震目录的10倍, 完备震级从1.2降至0.5。 结果显示, 北天山中段地震十分活跃, 主要分布在北天山山前霍尔果斯—玛纳斯—吐谷鲁背斜带南翼的浅部和南玛纳斯—齐古背斜带深部, 呼图壁地震震后地震活动性有增强的趋势。 研究期间沿背斜构造带走向地震分布不均, 霍尔果斯—玛纳斯—吐谷鲁背斜带西段地震活动多于东段, 南玛纳斯—齐古背斜带东段地震活动显著强于西段。 经过重定位, 发现研究区的地震事件主要发生在褶皱内部的“盲断层”上, 这些隐伏断裂与区域活动断裂和背斜构造共同组成的断层系孕育了北天山山前活跃的地震活动, 并可能成为未来强震的发震构造。     

川南地区龙马溪页岩气体滑脱效应实验研究

针对页岩储层气体滑脱效应特征及其影响机制不清问题，选取四川盆地长宁地区志留系龙马溪组页岩样品，开展了低温氮气吸附孔隙结构表征实验，并利用非稳态脉冲衰竭方法测量了不同围压下氦气、氮气在页岩岩心上的气体渗透率，分析了平均孔隙压力、气体类型、围压对滑脱效应的影响，建立了滑脱因子的预测关系式。结果表明:压力低于2.5 MPa时，页岩气体滑脱效应不能忽略。由于“分子筛效应”的影响，页岩克氏渗透率与测试流体介质类型有关，以氦气为流动介质测试得到的克氏渗透率大于以氮气为流动介质的测试结果。滑脱效应与气体类型有关，龙马溪组页岩的氦气滑脱因子约为氮气滑脱因子的1.7倍。利用滑脱因子计算得到围压为10~40 MPa时，氦气在页岩上的有效渗流孔径为113~166 nm，氮气的有效渗流孔径为66~99 nm，均远大于液氮吸附法测试的平均孔径。建立了龙马溪组页岩气体滑脱因子与克氏渗透率的幂函数关系，为页岩气流动模型的建立提供了基础。