一种滑坡涌浪计算方法技术

本发明专利技术涉及一种滑坡涌浪计算方法，将库岸边坡滑坡体入水到涌浪传播各阶段作为完整统一的过程进行分析研究，包括如下步骤：步骤1：根据已有地形、水位和滑坡体形状位置资料，建立三维有限元模型；步骤2：利用Navier‑Stokes方程描述水体、空气、滑坡体的相互作用，用改进的守恒式Level Set方法准确捕捉流体自由表面，真实模拟涌浪产生过程，为涌浪传播及灾害预测提供合理的前提条件；步骤3：采用沿水深积分的流体控制方程模拟涌浪传播过程，对涌浪运动轨迹和泛滥区域进行数值预报，评价涌浪灾害。本发明专利技术不仅对库岸边坡失稳破坏及涌浪次生灾害的预测具有重大的理论意义和实用价值，而且对其它类似的滑坡、泥石流、堰塞湖等地质灾害的数值模拟亦具有借鉴意义。

【技术实现步骤摘要】
一种滑坡涌浪计算方法
本专利技术涉及一种滑坡涌浪计算方法，属于水利工程灾害预测领域。
技术介绍
水利水电建设的长足发展为社会发展做出了重要贡献，但也带来一系列地质环境问题。相比于一般的边坡失稳问题，库岸边坡失稳会对下游的人民生命财产安全构成更加巨大的威胁。滑坡体滑入水库后产生的涌浪有可能翻越坝顶，淹没下游地区，甚至由于涌浪的冲击作用，大坝存在溃决的风险。在国内水利水电工程建设中，不乏水库滑坡失稳破坏并产生涌浪灾害的先例。如1961年3月柘溪水库塘岩光滑坡，165万m3的滑坡体冲入库区后产生了高达21m的涌浪，造成重大损失。随着一系列高坝大库在我国西南地区建立，水库库岸边坡稳定问题越来越吸引人们的注意。在此背景下，数值模拟作为灾害预测的手段之一，就显的非常必要和重要了。目前的数值模拟基本是首先用能量守恒原理估算滑坡体入水前的滑速，而后依据经验公式计算初始涌浪，最后用计算流体力学的方法求解涌浪的传播过程，实际上只是求解了已知初始条件和边界条件的计算流体力学问题。滑坡体入水产生涌浪的过程是一个复杂的水动力学问题，过程中涉及到三相介质的相互耦合、浸入、波浪破碎等现象，介质间的界面难于捕捉描述。现有的数值方法，包括基于浅水方程的解法，不能反映三相介质的强烈相互作用。主要的界面描述方法有VOF法、LevelSet法等，涌浪自由面的描述在计算流体力学中本身就是很大的难题，数值计算的准确性可以得到进一步改善。相对于涌浪产生来说，涌浪的传播过程研究成果已较为丰富。但此类计算的前提是已知首浪的各种要素，如浪高、速度等。从某种意义上说，现有方法是用一个粗估的初始涌浪数据作为初始计算条件，即使涌浪传播过程的模拟再精细，也很难得到合理可信的计算结果。尽管国内外学者对此已做了大量的研究工作，取得了丰富的研究成果，但基本上是着眼于某个方面的研究，还没有针对从库岸边坡失稳运动到涌浪产生与传播全过程系统的数值模拟。而实际上此类灾害的发生基本是在几十秒到几分钟极短的时间内，库岸边坡从失稳、滑动到涌浪产生再到涌浪传播应该是一个完整连续的过程，每个阶段的准确描述都将对最后灾害预报的可信度产生较大影响。
技术实现思路
本专利技术所要解决的技术问题是针对
技术介绍
的缺点，将库岸边坡滑坡体入水到涌浪传播各阶段作为完整统一的过程进行分析研究。采用原始的不可压缩粘性流体Navier-Stokes方程描述水体、空气、滑坡体的相互作用。采用改进的守恒式LevelSet方法，准确捕捉流体自由表面，真实模拟涌浪产生过程，为涌浪传播及灾害预测提供合理的前提条件。在此基础上，采用沿水深积分的流体控制方程模拟涌浪传播过程，对涌浪运动轨迹和泛滥区域进行数值预报，评价涌浪灾害。本专利技术是在已经利用动力时程强度折减法进行土质库岸边坡动力稳定分析，确定了库岸可能的失稳滑动范围的前提下进行的。本专利技术为实现上述专利技术目的采用如下技术方案：一种滑坡涌浪计算方法，包括如下步骤：步骤1：根据已有地形、水位和滑坡体形状位置资料，建立三维有限元模型；步骤2：结合Navier-Stokes方程和改进的守恒式LevelSet方法，计算流体各运动参数，捕捉自由面；步骤3：采用沿水深积分的流体控制方程模拟涌浪传播过程，对涌浪运动轨迹和泛滥区域进行数值预报。所述步骤2中对土质或松散介质坡体，基于“流体化”的特点，以流体力学基本方程Navier-Stokes方程配以非牛顿流体或拟流体本构关系描述其物理特性和运动过程，由于将滑坡体视作似流体，滑坡体、水体和空气三相介质可用同一本构模型表示：式中，τy0表示流体屈服应力，φ为内摩擦角，I2D是应变张量D的第二不变量，μ为黏性系数，η为表示流体特性的参量。在利用Navier-Stokes方程计算出的流场信息的基础上，采用改进的守恒式双重LevelSet方法捕捉滑坡体、水体和空气三者的交界面，反映三者间强烈的流固耦合作用，真实模拟滑坡体入水产生涌浪的过程，改进的守恒式双重LevelSet方法包含两个部分：指示函数的输运，和指示函数的重新初始化。在三相流中，各点的物理参数用指示函数H表示，各点的物理参数包括密度或黏性；引入指示函数H1、H2，滑坡体、水体和空气三相介质各自区域分别记作Ω1、Ω2、Ω3，滑坡体与水体、水体与空气、滑坡体与空气界面依次用Γ12、Γ23、Γ13表示，P表示物理参数，则各区域参数用统一的方程表示：P＝P1+(P2-P1)H1+(P3-P2)H2在Ω1内，H1＝H2＝0；在Ω2内，H1＝1，H2＝0；在Ω3内，H1＝H2＝1。交界面的位置确定方法，引入标记距离函数和指示函数H，或H＝0.5表示界面，重新初始化方程中包含界面的法线方向，由于指示函数H本身的局限性，不适宜用于求解远离H＝0.5处的界面法线方向，选取一个临界值δ，当时，利用计算界面法线方向，否则用H计算法线方向，具体过程表述如下：(1)求解两个指示函数的输运方程，得到和H；(2)对远离界面处的标记距离函数进行初始化；(3)计算界面弥散区域各点的法线方向；(4)对指示函数H进行初始化；(5)根据H计算靠近界面处的距离标记函数(6)进行物理特性插值，确定界面位置，进行下一步。交界面的位置确定过程中选取的临界值δ一般可以取0.1Δ～0.2Δ，其中Δ表示网格尺寸，当时表示远离界面，表示靠近界面。以步骤2中首浪产生数值模结果作为前提进行连贯的数值分析，用沿水深积分的Navier-Stokes方程即浅水方程作为控制方程求解涌浪的传播过程，通过步骤3的计算，可以得到各时刻各测点的水深数值，对涌浪的运动轨迹和泛滥区域进行数值预报，评价涌浪灾害。典型滑坡体截面产生的初始浪高计算方法如下：在滑坡体进入水体之后的任意时刻，取指示函数H值分布图，取两点，分别为点A和点B，设HA＜0.5，点A位于空气中，HB＞0.5，点B位于水体中，A点高程hA，B点高程hB，根据线性插值远离，容易得出，H＝0.5的点的高程为初始浪高h＝hO-hStill，其中hStill表示静水位。本专利技术采用上述技术方案具有如下有益效果：本专利技术针对库岸边坡失稳到涌浪产生与传播全过程进行数值模拟采用，将库岸边坡从失稳、滑动到涌浪产生再到涌浪传播看作是一个完整连续的过程，采用Navier-Stokes方程和改进的守恒式LevelSet方法计算出首浪的各参数，在此基础上，采用浅水方程计算涌浪的传播，对涌浪运动轨迹和泛滥区域进行数值预报，进行涌浪灾害评价。附图说明图1是本专利技术三维有限元网格图，图2是本专利技术典型断面密度分布图，图3是本专利技术三相流指示函数H1、H2示意图，图4是本专利技术首浪高度计算示意图，图5是本专利技术河道截面处涌浪水位随时间变化过程(靠近滑坡体入水处)，图6是本专利技术河道截面处涌浪水位随时间变化过程(靠近大坝)。具体实施方式下面结合附图和具体实施方式，进一步阐明本专利技术。应理解下述具体实施方式仅用于说明本专利技术而不用于限制本专利技术的范围。本专利技术采用原始的不可压缩粘性流体Navier-Stokes方程描述水体、空气、滑坡体的相互作用。采用改进的守恒式LevelSet方法，准确捕捉流体自由表面，真实模拟涌浪产生过程。在此基础上，采用沿水深积分的流体控制方程模拟涌浪传播过程，对涌浪运动轨迹和泛滥区域进行数值预报，评价涌浪灾害。一种滑坡涌浪计算方法本文档来自技高网...

【技术保护点】
一种滑坡涌浪计算方法，其特征在于包括如下步骤：步骤1：根据已有地形、水位和滑坡体形状位置资料，建立三维有限元模型；步骤2：结合Navier‑Stokes方程和改进的守恒式Level Set方法，计算流体各运动参数，捕捉自由面；步骤3：采用沿水深积分的流体控制方程模拟涌浪传播过程，对涌浪运动轨迹和泛滥区域进行数值预报。

【技术特征摘要】
1.一种滑坡涌浪计算方法，其特征在于包括如下步骤：步骤1：根据已有地形、水位和滑坡体形状位置资料，建立三维有限元模型；步骤2：结合Navier-Stokes方程和改进的守恒式LevelSet方法，计算流体各运动参数，捕捉自由面；步骤3：采用沿水深积分的流体控制方程模拟涌浪传播过程，对涌浪运动轨迹和泛滥区域进行数值预报。2.根据权利要求1所述的滑坡涌浪计算方法，其特征在于所述步骤2中对土质或松散介质坡体，基于“流体化”的特点，以流体力学基本方程Navier-Stokes方程配以非牛顿流体或拟流体本构关系描述其物理特性和运动过程，由于将滑坡体视作似流体，滑坡体、水体和空气三相介质可用同一本构模型表示：式中，τy0表示流体屈服应力，φ为内摩擦角，I2D是应变张量D的第二不变量，μ为黏性系数，η为表示流体特性的参量。3.根据权利要求2所述的滑坡涌浪计算方法，其特征在于在利用Navier-Stokes方程计算出的流场信息的基础上，采用改进的守恒式双重LevelSet方法捕捉滑坡体、水体和空气三者的交界面，反映三者间强烈的流固耦合作用，真实模拟滑坡体入水产生涌浪的过程，改进的守恒式双重LevelSet方法包含两个部分：指示函数的输运，和指示函数的重新初始化。4.根据权利要求3所述的滑坡涌浪计算方法，其特征在于在三相流中，各点的物理参数用指示函数H表示，各点的物理参数包括密度或黏性；引入指示函数H1、H2，滑坡体、水体和空气三相介质各自区域分别记作Ω1、Ω2、Ω3，滑坡体与水体、水体与空气、滑坡体与空气界面依次用Γ12、Γ23、Γ13表示，P表示物理参数，则各区域参数用统一的方程...

【专利技术属性】
技术研发人员：毛佳，赵兰浩，郭博文，李同春，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32