CN111814321A - 一种再现堰塞坝堆积特性的数值模拟方法 - Google Patents
一种再现堰塞坝堆积特性的数值模拟方法 Download PDFInfo
- Publication number
- CN111814321A CN111814321A CN202010586131.8A CN202010586131A CN111814321A CN 111814321 A CN111814321 A CN 111814321A CN 202010586131 A CN202010586131 A CN 202010586131A CN 111814321 A CN111814321 A CN 111814321A
- Authority
- CN
- China
- Prior art keywords
- dam
- numerical simulation
- accumulation
- damming
- reproducing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种再现堰塞坝堆积特性的数值模拟方法，包括如下步骤：对堰塞坝坝体材料进行3D扫描，获取块石点云数据库；根据块石点云数据库获取材料级配曲线特征；根据级配曲线特征设置球形颗粒级配滑体或者块石颗粒级配滑体；根据现场滑坡堰塞坝特点，在颗粒流软件中设置计算参数信息；通过对坝体材料参数标定，选取数值模拟合适的计算参数；堰塞坝堆积完成后，测量相关堰塞坝坝体堆积信息；对堰塞坝坝体堆积信息进行回归分析，给出预测模型。本发明将数值模拟技术应用到堰塞坝堆积演化过程当中，通过数值模拟的方法再现堰塞坝堆积特性，构建坝体形态预测模型，有效避免了物理模型试验耗时耗力的缺点，可以方便、准确的认知堰塞坝堆积特性。
Description
技术领域
本发明属于边坡工程与滑坡治理领域，具体涉及一种再现堰塞坝堆积特性的数值模 拟方法。
背景技术
地震或降雨引起的滑坡，特别是在高山峡谷地区，危害性极大，容易形成堵塞河流的堰塞坝，威胁上下游群众的生命财产安全。
目前，对堰塞坝堆积形态的研究主要集中在物理模型试验研究，但是试验研究耗时 耗力，且不易总结出较好的规律。
所以，需要一个新的技术方案来解决这个问题。
发明内容
发明目的：为了克服现有技术中存在的无法准确快速描述堰塞坝堆积演化过程的不 足，提供一种再现堰塞坝堆积特性的数值模拟方法，能够准确快速地描述堰塞坝堆积特性。
技术方案：为实现上述目的，本发明提供一种再现堰塞坝堆积特性的数值模拟方法， 包括如下步骤：
S1：对堰塞坝坝体材料进行3D扫描，获取块石点云数据库；
S2：根据块石点云数据库获取材料级配曲线特征；
S3：根据级配曲线特征设置球形颗粒级配滑体或者块石颗粒级配滑体；
S4：根据现场滑坡堰塞坝特点，在颗粒流软件中设置计算参数信息；
S5：通过对坝体材料参数标定，选取数值模拟合适的计算参数；
S6：堰塞坝堆积完成后，测量相关堰塞坝坝体堆积信息；
S7：对堰塞坝坝体堆积信息进行回归分析，给出预测模型。
进一步的，所述步骤S4中计算参数信息包括滑体、滑面和河床。
进一步的，所述颗粒流软件中滑体为球形颗粒或者块石颗粒，滑面和河床为刚性墙。
进一步的，所述步骤S5中通过自由堆积试验来确定数值模拟合适的计算参数。
进一步的，所述步骤S6中堰塞坝坝体堆积信息包括坝高、坝长、坝宽、上下游绝 对倾角、坝体前倾角。
进一步的，所述步骤S1中块石点云数据库的格式为stl。
进一步的，坝体堆积形态预测模型主要通过多元回归来实现，参数分为输入组与输 出组。
数值模拟技术的发展为研究滑坡运动和堰塞坝堆积特性提供了新的途径，可以实现 启动到堆积的全过程模拟与监测，凭借其能模拟大变形、大位移的优势，以及实现滑体颗粒的平动、转动、分离和流动等特点，可以在滑坡堆积运动种发挥其巨大作用。因此， 采用数值模拟的方法再现堰塞坝堆积演化过程是一种较好的手段，能够反映出堆积的实 际情况。
有益效果：本发明与现有技术相比，将数值模拟技术应用到堰塞坝堆积演化过程当 中，通过数值模拟的方法再现堰塞坝堆积特性，构建坝体形态预测模型，有效避免了物理模型试验耗时耗力的缺点，可以方便、准确的认知堰塞坝堆积特性。
附图说明
图1为本发明的流程图；
图2为3D扫描块石的点云数据；
图3为堰塞坝坝体材料级配曲线设置；
图4为PFC3D中搭建的数值模型；
图5为堰塞坝坝体材料自然堆积角标定过程；
图6为堰塞坝堆积形态图。
具体实施方式
下面结合附图和具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本 发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图1所示，本发明提供一种再现堰塞坝堆积特性的数值模拟方法，包括如下步骤：
(1)对现场堰塞坝坝体材料进行3D扫描，获取stl格式的块石点云数据库；堰塞 坝坝体堆积材料不宜过小或过大，尺寸需满足3D扫描的要求。
(2)获取材料级配曲线特征，根据级配曲线特征设置球形颗粒级配滑体或者块石颗粒级配滑体；滑体设置可以根据计算机计算效率、块石所占比例、细颗粒土所占比例 等因素将滑体设置为球形颗粒或者是块石颗粒，但都需满足堰塞坝坝体材料级配曲线。
(3)根据现场滑坡堰塞坝特点，在颗粒流软件PFC3D中设置滑体、滑面、河床等 计算参数信息；PFC3D中的滑体为球形颗粒或者块石颗粒，滑面和河床为刚性墙，主要 的考虑因素包括滑体体积含量、滑面长度、滑面角度、滑面出口宽度、河床形状和河谷 角度。
(4)通过对坝体材料参数标定，选取数值模拟合适的计算参数；主要是通过选择合适的模型，通过自由堆积试验来确定。
(5)堰塞坝堆积完成后，测量相关坝体信息，包括坝高、坝长、坝宽、上下游绝 对倾角、坝体前倾角。
(6)对堰塞坝坝体信息进行回归分析，给出预测模型；堆积形态预测模型主要通过多元回归来实现，参数分为输入组与输出组。
本实施例中将上述方法实际应用于现场堰塞坝坝体材料，其具体步骤如下：
1、对现场堰塞坝坝体材料进行3D扫描，构建了可以导入PFC3D的stl格式的块石点云数据库，块石点云数据库具体如图2所示。
2、根据现场级配曲线调查，分别绘制了堰塞坝三条级配曲线，分别为平均线、上包线和下包线，考虑到计算效率，本次实施例中选用球形颗粒，具体颗粒级配曲线如图 3中所示。
3、设置滑体、滑面、河床等信息，本实施例中滑体体积设置为0.4m*0.4m*0.4m， 滑面宽度1m，滑出口宽度1m，河谷形状设置为三角形，河床角度0°，坡面角度30°， 具体如图4所示。
4、本实施例中堰塞坝坝体材料主要为砾石，其自然堆积角一般为25°-40°，对自然堆积角的坝体材料进行标定，采用PFC3D中的线性阻抗模型，标定过程如图5所示， 具体参数参照表1。
表1数值模型计算参数
| 最小颗粒半径 | 1.0cm | 摩擦系数 | 0.4 |
| 颗粒粒径比 | 1.66 | 阻抗摩擦系数 | 0.6 |
| 颗粒密度 | 2650kg/m<sup>3</sup> | 法向阻尼 | 0.2 |
| 孔隙率 | 0.5 | 切向阻尼 | 0.2 |
| 接触法向刚度 | 5.0e7N/m | 接触切向刚度 | 5.0e7N/m |
5、堰塞坝堆积完成后，堰塞坝堆积形态如图6所示，包括纵河向剖面图、横河向 剖面图和俯视形态图。
6、分别统计不同工况下的坝高、坝长、上下游倾角和坝体前倾角，以滑动距离l、出口宽度a、坡面角度正切值tanη和河谷角度正切值作为tanθ作为输入数据，以坝体高 度h、长度L、上游倾角正切值tanβu、下游倾角正切值tanβd和前倾角正切值
为输 出数据，对其变化规律进行多项式回归分析，回归公式如下所示；
h＝-0.0509+0.06961a+0.0124tanη+0.0497tanθ+0.0164l2-0.2766a2-0.1066la
L＝1.0214+0.07728a+1.0129tanθ+0.1506tanη2
tanβu＝0.04629l+0.7277a-0.7569tanη+0.2793tanη2+2.5468tanθ2-0.479la
tanβd＝0.8366-0.6602tanη-0.7561tanθ+0.0241l2+0.2534tanη2-0.1516la
本实施例中所获得的预测模型就是上面的回归公式，可根据实际试验条件对坝体形 态进行预测。预测模型包含输入值和输出值，预测模型是通过将输入值代入回归公式，得到输出值，则输出值就是根据模型得到的预测值。
Claims (5)
1.一种再现堰塞坝堆积特性的数值模拟方法，其特征在于：包括如下步骤：
S1：对堰塞坝坝体材料进行3D扫描，获取块石点云数据库；
S2：根据块石点云数据库获取材料级配曲线特征；
S3：根据级配曲线特征设置球形颗粒级配滑体或者块石颗粒级配滑体；
S4：根据现场滑坡堰塞坝特点，在颗粒流软件中设置计算参数信息；
S5：通过对坝体材料参数标定，选取数值模拟合适的计算参数；
S6：堰塞坝堆积完成后，测量相关堰塞坝坝体堆积信息；
S7：对堰塞坝坝体堆积信息进行回归分析，给出预测模型。
2.根据权利要求1所述的一种再现堰塞坝堆积特性的数值模拟方法，其特征在于：所述步骤S4中计算参数信息包括滑体、滑面和河床。
3.根据权利要求2所述的一种再现堰塞坝堆积特性的数值模拟方法，其特征在于：所述颗粒流软件中滑体为球形颗粒或者块石颗粒，滑面和河床为刚性墙。
4.根据权利要求1所述的一种再现堰塞坝堆积特性的数值模拟方法，其特征在于：所述步骤S5中通过自由堆积试验来确定数值模拟合适的计算参数。
5.根据权利要求1所述的一种再现堰塞坝堆积特性的数值模拟方法，其特征在于：所述步骤S6中堰塞坝坝体堆积信息包括坝高、坝长、坝宽、上下游绝对倾角、坝体前倾角。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010586131.8A CN111814321B (zh) | 2020-06-24 | 2020-06-24 | 一种再现堰塞坝堆积特性的数值模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010586131.8A CN111814321B (zh) | 2020-06-24 | 2020-06-24 | 一种再现堰塞坝堆积特性的数值模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111814321A true CN111814321A (zh) | 2020-10-23 |
| CN111814321B CN111814321B (zh) | 2023-05-23 |
Family
ID=72844859
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010586131.8A Active CN111814321B (zh) | 2020-06-24 | 2020-06-24 | 一种再现堰塞坝堆积特性的数值模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111814321B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113380119A (zh) * | 2021-06-01 | 2021-09-10 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种测定堰塞体材料斜坡运动和谷底堆积的试验装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108090670A (zh) * | 2017-12-15 | 2018-05-29 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种基于水流侵蚀能力的堰塞坝形成评估方法 |
| CN110532682A (zh) * | 2019-08-28 | 2019-12-03 | 四川建筑职业技术学院 | 一种泥石流堰塞坝溃决时流深预测方法 |
| CN110532684A (zh) * | 2019-08-28 | 2019-12-03 | 四川建筑职业技术学院 | 一种泥石流堰塞坝溃决后的冲出量预测方法 |
-
2020
- 2020-06-24 CN CN202010586131.8A patent/CN111814321B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108090670A (zh) * | 2017-12-15 | 2018-05-29 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种基于水流侵蚀能力的堰塞坝形成评估方法 |
| CN110532682A (zh) * | 2019-08-28 | 2019-12-03 | 四川建筑职业技术学院 | 一种泥石流堰塞坝溃决时流深预测方法 |
| CN110532684A (zh) * | 2019-08-28 | 2019-12-03 | 四川建筑职业技术学院 | 一种泥石流堰塞坝溃决后的冲出量预测方法 |
Non-Patent Citations (2)
| Title |
|---|
| HAOWU 等: "Laboratory-scale investigation of the 3-D geometry of landslide dams in a U-shaped valley", 《ENGINEERING GEOLOGY》 * |
| YUANYUAN ZHOU 等: "3D DEM investigation on the morphology and structure of landslide dams formed by dry granular flows", 《ENGINEERING GEOLOGY》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113380119A (zh) * | 2021-06-01 | 2021-09-10 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种测定堰塞体材料斜坡运动和谷底堆积的试验装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111814321B (zh) | 2023-05-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Quantitative analysis of microscopic structure and gas seepage characteristics of low-rank coal based on CT three-dimensional reconstruction of CT images and fractal theory | |
| Chinnarasri et al. | Embankment dam breach and its outflow characteristics | |
| CN107091911A (zh) | 一种河流水质预测方法 | |
| CN114218840B (zh) | 河口航道水沙运动及其地形演变整体建模及可视化系统 | |
| CN111814321A (zh) | 一种再现堰塞坝堆积特性的数值模拟方法 | |
| CN113128102B (zh) | 铁路道床道砟侵入物多尺度分析模型快速构建的方法 | |
| Imre et al. | Case studies and benchmark examples for the use of grading entropy in geotechnics | |
| CN113585161A (zh) | 一种冲积河流深泓横向迁移预测模型的构建方法 | |
| CN114357903B (zh) | 基于河床形态参数的冲积河流阻力计算、模拟方法及装置 | |
| CN116151421A (zh) | 抽水蓄能电站选址方法、数据处理终端及可读存储介质 | |
| CN115481467A (zh) | 复杂地质沉管隧道基槽分土质超欠挖分析方法及装置 | |
| CN113916916A (zh) | 一种用于页岩数字岩心三维渗流-颗粒流动耦合的模拟方法 | |
| Abdel-Kawy et al. | Effects of core characteristics on seepage through earth dams | |
| CN111104746A (zh) | 一种基于小波分析的河漫滩高程确定方法 | |
| Rotnicka et al. | Vertical profiles of aeolian mass flux above different sand surfaces and sand surfaces covered with pebbles | |
| Zhao et al. | Seepage–fractal model of embankment soil and its application | |
| Reitano et al. | Erosional response of granular material in landscape models | |
| Deshpande et al. | Advent of sheet flow in suction affected alluvial channels | |
| CN114611431A (zh) | 一种资料贫乏地区泥沙长时间序列重构方法 | |
| CN114547922A (zh) | 基于多边形网格的城市内涝分析方法、装置和存储介质 | |
| Lv et al. | Simulation Study on hydrological process of soil cracks in open-pit coal mine dump | |
| Luo et al. | Study on the Evolution of Tailings Dam Break Disaster under Complex Environment | |
| Barzegari et al. | Numerical investigation of cavitation on spillways. A case study: Aydoghmush dam | |
| Kökpinar et al. | Physical and numerical modeling of shoreline evaluation of the Kizilirmak river mouth, Turkey | |
| Nurafza et al. | Facies connectivity modelling: analysis and field study |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |