论文摘要
黄土高原约占我国国土总面积的6.6%,由于黄土自生结构的特殊性,在外力因素的干扰,尤其是水的作用下,会表现出屈服、软化、强度降低及湿陷等特性,导致地质灾害频频发生,其灾害总数约占据全国总数的三分之一。而对于甘肃黑方台地区,由于长期农业漫灌,地下水位不断抬升,导致平均每年有3-5次滑坡地质灾害发生,严重制约了当地的社会、经济、文化发展。目前对于黑方台黄土滑坡的致灾因素、灾变机理、成灾模式以及演化过程等研究,已有大量的研究成果,而对于黑方台黄土滑坡的滑距预测,却鲜有报道,因此展开对黑方台黄土滑坡的空间滑距预测研究,可以为该地区的滑坡预警及防控提供科学依据。本文围绕着滑距预测这一目的,基于黑方台2015—2017年发生的7起典型黄土滑坡滑前及滑后地形数据,首先通过离散元数值模拟软件对其中6起滑坡进行以参数为变量的系列正交试验,得到各组参数及其对应的模拟形态(滑源区横向宽、滑源区纵向宽、滑坡的壁高、滑坡整体落差、滑距、堆积区横向宽);然后以此作为样本进行神经网络训练,确定出模拟形态与参数之间的映射关系;进而利用此映射关系反演得到这6起滑坡的最优模拟参数,并以这6组最优模拟参数的范围作为该地区黄土滑坡数值模拟参数的取值范围;接着通过对第7个滑坡的正演结果,验证了该参数取值范围的可靠性;最后以该取值范围的均值作为参数对陈家7#潜在滑坡进行了滑距预测。主要得到以下结论:(1)确定了黑方台黄土滑坡离散元数值模拟所应采用的本构模型、计算时间步长、网格搜索大小以及离散元软件等。根据对离散元理论的研究可知,不同的接触模型对应的力学机理不同,根据颗粒运动过程的受力特征,得出黑方台黄土滑坡颗粒与颗粒之间的本构模型应采用JKR含湿颗粒接触模型,颗粒与几何体之间的接触模型应采用RVD滚动接触模型,时间步长的确定方法应采用瑞利波法,网格的搜索边长值应设为1.5dmax,最后通过离散元软件中接触模型库、颗粒接触搜索算法以及时间步长计算方法等三方面的比较,最终确定选用EDEM离散元软件对黑方台黄土滑坡进行模拟研究。(2)实现了EDEM离散元模型的颗粒快速平衡填充。由于目前大多数离散元程序不能直接建模或建模过程十分的繁琐,故本文通过Surfer、Polywork等软件对黑方台的无人机数据进行处理,构建出了滑坡的三维模型,并结合离散元PFC软件以及EDEM软件中API开发的快速填充模块,采用“逐层brick填充法”实现对EDEM离散元模型的颗粒快速平衡填充;然后根据前人力学试验研究、室内休止角试验、以及GEMM材料数据库确定出材料颗粒粒径、泊松比、弹性模量、密度、颗粒与滑源区几何体之间的静态摩擦系数、滚动摩擦系数等6个参数的取值,得到颗粒与颗粒之间、颗粒与非滑源区几何体之间的滚动摩擦系数、静摩擦系数、模型整体的碰撞恢复系数以及JKR模型的能量密度等6个参数的合理取值范围,为正交试验的设计提供依据。(3)利用神经网络方法,反演出6起滑坡对应的最优参数,并以这6组最优参数作为黑方台滑坡的参数取值范围,通过第7起滑坡的正演,验证了该取值范围的准确性。由于离散元中宏观与细观力学参数存在着很高的非线性,为快速、有效的得到能反应实际情况的离散元细观参数,本文首先通过对参数的敏感性分析,确定出JRK模型的表面密度能应为100J/m2,然后根据“敏感性大的参数设置较小变化梯度,敏感性小的参数设置较大变化梯度”的原则进行正交试验设计,通过6起滑坡共1548次的数值计算,得到不同参数下所对应的模拟形态;然后以此作为神经网络训练样本,构建出模拟形态与参数之间的映射关系模型,随后利用该模型实现这6起滑坡的参数反演,并利用这6起滑坡反演得到的最优参数作为黑方台黄土滑坡的参数取值范围;最后通过第七起滑坡的正演,验证了该取值范围的准确性,避免了相同地质条件下模拟时人工标定力学参数所带来的复杂性和随机性。(4)首次提出“限定后缘”有限元强度折减法对黑方台黄土滑坡的潜在滑动面进行搜索。“限定后缘”有限元强度折减法是以裂缝及落水洞的分布作为后缘控制性边界,利用位移等值线作为潜在滑面划分的依据对潜在滑面进行搜索;然后通过对2019.03焦家6#滑坡、2015.01焦家4#滑坡以及2017.10党川9#滑坡滑面的搜索,证明了该方法的可行性、适用性和准确性,为简便、快速、准确的搜索出黑方台黄土滑坡的潜在滑动面提供一种新的思路。(5)实现陈家7#潜在滑坡的滑距及运动过程预测分析。为确定出陈家7#潜在滑坡的滑体形态,本文通过无人机影像及现场的实际调查,以裂缝及落水洞的分布作为控制性边界,圈定出CJ7#滑坡可能发生的后缘边界,利用“限定后缘”强度折减法对其潜在滑面进行搜索,确定出CJ7#滑坡的潜在滑面形态,然后利用ItasCAD三维地质建模软件构建出潜在滑坡体的三维离散元模型,接着以反演出的陈家段已有滑坡离散元细观参数的均值作为参数,通过离散元法得到CJ7#滑坡沿限定后缘边界发生滑动时的滑距约为380m,堆积最大宽度约为97m,滑坡前缘剪出口与滑坡堆积最前缘的高差约为59m,最大堆积厚度约为10m。滑坡发生时,各区域的启动响应先后顺序为:B3>B2>B1>T2>T1,滑坡发生后,各区域的相对前后及堆积层序位置与滑前基本保持一致。
论文目录
文章来源
类型: 硕士论文
作者: 周小棚
导师: 许强
关键词: 离散元,有限元,神经网络,参数反演,滑面搜索,滑距预测
来源: 成都理工大学
年度: 2019
分类: 基础科学,工程科技Ⅱ辑
专业: 地质学,地质学,工业通用技术及设备
单位: 成都理工大学
基金: 国家自然科学基金重点项目—溃散性滑坡成因机理,监测预警与定量风险评价(NO.41630640)
分类号: P642.22
DOI: 10.26986/d.cnki.gcdlc.2019.000865
总页数: 118
文件大小: 14757K
下载量: 152
相关论文文献
- [1].海原地震典型黄土滑坡群形成机制与运动特征研究[J]. 成都大学学报(自然科学版) 2019(04)
- [2].灌溉诱发突发性黄土滑坡机理研究[J]. 山地学报 2020(01)
- [3].陕西黄土地区滑坡的一般勘察方法和工程处理措施[J]. 中国金属通报 2020(03)
- [4].伊犁苏阿苏沟东岸黄土滑坡群发育特征及成因分析[J]. 新疆地质 2020(02)
- [5].黄土滑坡活动的时间分布规律研究[J]. 绿色科技 2019(06)
- [6].新疆伊犁则克台河黄土滑坡分布规律及控制要素分析[J]. 青海交通科技 2018(05)
- [7].黄土滑坡机理分析与防治对策[J]. 河南建材 2019(04)
- [8].陕西泾阳南塬多序次黄土滑坡演化特征及成灾模式[J]. 地质科技情报 2019(05)
- [9].兰州市降雨型黄土滑坡灾害空间分布特征[J]. 测绘科学 2016(12)
- [10].废弃窑洞诱发黄土滑坡机理与回填技术研究[J]. 洛阳理工学院学报(自然科学版) 2017(02)
- [11].延安市宝塔区群发性浅层黄土滑坡的地形条件研究[J]. 甘肃水利水电技术 2017(09)
- [12].山西省吕梁典型黄土滑坡特性分析[J]. 山西建筑 2017(31)
- [13].陕北黄土高原不同地貌类型区黄土滑坡频率分布[J]. 地球科学 2016(02)
- [14].天水市麦积区税湾地震黄土滑坡特征及其形成机制[J]. 地质力学学报 2016(01)
- [15].新疆伊犁典型特大型黄土滑坡群成因机制及变形特征[J]. 水利与建筑工程学报 2016(04)
- [16].黄土滑坡研究现状及其展望[J]. 地下水 2015(01)
- [17].基于非饱和土理论的黄土滑坡降雨预警模型[J]. 水力发电 2020(08)
- [18].基于神经网络的离散元细观参数标定方法研究——以甘肃黑方台黄土滑坡为例[J]. 岩石力学与工程学报 2020(S1)
- [19].论黄土滑坡和崩塌灾害发育特征及其稳定性[J]. 山西建筑 2017(33)
- [20].泾阳县泾河南塬黄土滑坡形成机理分析[J]. 地质灾害与环境保护 2018(01)
- [21].降雨诱发浅表层黄土滑坡机理实验研究[J]. 工程地质学报 2017(05)
- [22].灌溉诱发型黄土滑坡离心模型试验研究[J]. 岩土工程学报 2016(S2)
- [23].泥流型黄土滑坡的特征与成因[J]. 冰川冻土 2015(01)
- [24].灌溉型黄土滑坡变形破坏特征的离心机模型研究[J]. 建筑安全 2015(05)
- [25].彭建兵 黄土高原灾害“追凶”30年[J]. 陕西画报 2020(02)
- [26].“黄土滑坡致灾空间预测研究”项目获国家基金资助[J]. 甘肃科学学报 2013(04)
- [27].黄土滑坡诱发因素及其形成机理研究[J]. 工程地质学报 2011(04)
- [28].多级旋转黄土滑坡基本类型及特征分析[J]. 工程地质学报 2011(05)
- [29].应变能密度在黄土滑坡临滑预报中的应用[J]. 城市道桥与防洪 2010(04)
- [30].对陕北典型的黄土滑坡工程勘察的探讨[J]. 科技风 2009(09)