基于离散元方法对走滑拉分盆地演化及次级断裂扩展过程研究

拉分盆地是一种与走滑断裂带密切相关的特殊拉张构造,因其重要的构造意义,及其与火山活动、中小地震群集、特殊的成矿作用间的伴生关系而受到研究者的高度重视。关于拉分盆地的形成演化过程,已有较多的研究成果,但是由于研究手段的限制,缺少对盆地演化中次级断裂扩展过程的研究。基于离散元的数值计算方法是研究断裂扩展方式的理想方法。本文采用基于离散元的颗粒流方法,揭示纯走滑拉分盆地发育过程中的断裂扩展和连接过程,为拉分盆地演化机理和断裂扩展提供新的研究方法。同时,根据主走滑断层与块体运动方向的夹角不同,建立不同的张扭性拉分盆地模型,系统研究张扭性盆地的断裂扩展和演化机理。将上述理论研究结果与死海盆地等经典拉分盆地实例相结合,探讨了死海盆地、土耳其Cinarcik盆地、哥伦比亚El Paraiso盆地等的形成演化机理和断裂扩展方式。      

基于堆积层厚度变化的侏罗系顺向岸坡变形机制转换研究

三峡库区侏罗系顺向岸坡堆积体滑坡众多,其滑动模式存在一定差异。首先统计分析了192个三峡库区侏罗系层位发育的堆积体滑坡滑体及碎石土的工程地质性质和强度参数。在此基础上,运用数值分析软件对堆积体厚度变化引起的滑坡变形机制进行模拟分析。结果表明:堆积体厚度范围在15m及以下滑坡会沿着岩土界线面滑动、15~35m时滑坡会沿着层内剪切面滑动;厚度范围在35m以上时,堆积体滑坡内部可能存在着多层滑带,即滑坡可能沿着层内剪切面滑动或者沿着岩土界线面滑动。堆积体厚度范围在15m及以下时滑坡的治理措施可采用布置抗滑桩;厚度范围15~35m时可采用排水+布置抗滑桩的滑坡治理措施;厚度范围在35m以上时,可采用滑坡前期监测预报+后期根据滑坡发育情况相结合的滑坡防治措施。     

FLAC3D在复杂工程条件下的应用

利用建筑垃圾堆山造景是现代城市的一个新方向。但由于堆山工程占地范围大,荷载大,往往会产生一系列的工程地质问题,在软土地区这个问题更为突出。为使堆山工程能顺利进行,必须对地基进行一定的处理,而常规计算只能选择典型剖面或者地段进行局部分析,不能反映真实的三维影响特征,但数值分析可以有效解决这一问题。利用有限差分软件FLAC3D,对天津某堆山工程进行分析,为优化地基处理方案提供有力依据。     

地表暖化影响下温度示踪地下水流速方法

为构建温度示踪方法测算地下水流速技术体系,并应用于区域地下水资源评价,基于最小二乘法和垂向一维非稳定流水-热运移方程数值解法,提出地表暖化情形下地下水流速计算方法,并对雷州半岛东北部地下水流速进行测算。结果表明:研究区域地下水补给速度为0.796m/a,入渗以西北部降水和运河渗漏为主;地下水排泄速度为0.269m/a,排泄入海主要发生在东海岛、南三岛和硵州岛附近。温度示踪解析区域地下水流动情况与地下水位分布情况基本一致,观测和计算地温数据具有较强相关性(R^2>0.50)和较低均方根误差(均值0.748),表明提出方法率定得到的地下水流速具有较强的可靠性。参数敏感性分析结果表明,地质体热扩散率和地表温度均对地温计算结果产生较明显的影响,参数的准确率定对利用地温计算地下水流速十分重要。     

基于模拟-优化模型的山东羊庄盆地地下水可开采量研究

为了准确量化与优化地下水水源地开采量,在分析山东羊庄盆地水资源分布特点、开发利用现状及存在问题的基础上,以区域经济发展和水资源可持续利用为前提,通过地质、水文地质资料的整理分析,建立起整个盆地地下水流概念模型;采用Visual Modflow模拟软件,并与优化模型软件Ground-Water Management Process(GWM)相耦合,以水源地开采量最大为目标,以水位不低于限定的控制水位为约束条件,评价出盆地最大可开采量为5 730万m^3/a;进而确定出宋屯、后石湾、南庄、魏庄、西石楼、许坡、龙山头、羊庄、东于、段庄以及王杭水源地的优化开采量分别为0.18万、1.40万、2.56万、1.86万、3.80万、1.70万、0.40万、0.40万、0.50万、1.40万和1.50万m^3/d。     

超声波振动下不同应力条件对岩石强度影响的试验

振动有助于碎岩,以往关于振动碎岩机理的研究大多在低频率段展开。为填补超高频率段下振动碎岩机理的空白,采用单轴动静组合加载模式,开展了超声波振动下不同应力条件对岩石强度影响的试验研究,其中超声波振动频率为20kHz,预压范围为100～500N。研究结果表明:当预压小于200N时,岩石内部应力状态无法满足强度准则,岩石强度下降不明显;当预压大于等于200N时,岩石强度随振动时间的增加而逐渐降低且存在最优预压力值(400N)使得岩石强度最低。缩短振动频率与岩石固有频率的差值有利于提高超声波振动碎岩效率。     

缝网间距对高温岩体储留层温度影响规律模拟

为研究高温岩体储留层注水过程中缝网间距对储留层岩体温度的影响规律,建立了考虑热对流和热传导效应的热流固耦合力学模型。以水平井多簇缝网地热开发为例,将缝网及其周围岩体简化为等效的多孔介质,利用ABAQUS二次开发功能实现了流固耦合和热对流耦合。模拟结果表明：储留层岩体温度的变化按注水持续时间分为3个阶段：初温保持阶段,该阶段储留层岩体温度保持原始温度,缝网间距的不同对储留层岩体整体温度无影响;快速降温阶段,该阶段储留层岩体温度下降随注水时间呈现负指数变化规律,缝网间距与储留层岩体温度呈正相关关系;温度缓慢下降平稳阶段,该阶段初期缝网间距与储留层岩体温度的关系和快速降温阶段一致,但后期缝网间距与储留层岩体温度呈负相关关系。空间上,在岩体温度高的区域,缝网间距越小,岩体温度变化梯度越大;在岩体温度较低的区域,缝网间距越小,岩体温度变化梯度越小。     

基于LAHARZ的泥石流堆积范围预测模型的建立及应用

为了探究泥石流的堆积范围,利用LAHARZ软件,对北京市密云县泥石流沟喇嘛栅子南沟进行了数值模拟。结合泥石流沟小流域1：10 000数字高程模型图,模拟了泥石流的堆积范围。首先利用中国部分地区泥石流体积和堆积范围的数据资料,获得了泥石流体体积与其堆积范围的新的统计模型B=11.42V^0.7156;然后通过模拟沟道与实际沟道的对比,确定了最佳沟道阈值为15 000;再结合现场调查统计和降雨历史资料,确定了10年、20年、50年和100年一遇暴雨条件的泥石流体积值,分别为56 500、72 900、94 200和113 100 m³;最后在此基础上对该条泥石流沟的堆积范围进行了预测。结果表明,100年一遇的暴雨条件下泥石流堆积面积为48 729 m²,到达最远距离约为490 m,已影响下游村庄。     

基于NAD算法的声波方程时间四阶差分解法

波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。     

南海夏季风爆发的数值预报试验

分析了1986年南海夏季风爆发的环流演变特征,由经向风速剖面图看出,南风首先在中南半岛迅速加强,然后向南海发展,南海夏季风爆发同孟加拉湾低压的发展密切相关。通过地形和非绝热单因子敏感性数值预报试验表明,地形作用和凝结潜热等非绝热作用对南海夏季风的爆发都是很重要的,包含有这两种作用在内的控制试验成功地预报出了南海夏季风爆发的中期演变过程,单独的地形作用或非绝热作用都不能预报出南海夏季风的爆发。