CN106952000A - 一种喀斯特区域滑坡灾害风险动态评估方法 - Google Patents
一种喀斯特区域滑坡灾害风险动态评估方法 Download PDFInfo
- Publication number
- CN106952000A CN106952000A CN201710205591.XA CN201710205591A CN106952000A CN 106952000 A CN106952000 A CN 106952000A CN 201710205591 A CN201710205591 A CN 201710205591A CN 106952000 A CN106952000 A CN 106952000A
- Authority
- CN
- China
- Prior art keywords
- inducing factors
- flood inducing
- rock mass
- threedimensional model
- stability degree
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Physics (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种喀斯特区域滑坡灾害风险动态评估方法,其包括:选取降雨量、位移矢量和沉降量作为致灾因子;对致灾因子、山体稳定度进行量化、归一化处理和权重确定;以大量区域灾害统计及人工降雨侵蚀试验结果为训练样本,采用改进BP网络计算方法对神经网络进行训练直至达到误差要求,获得致灾因子与山体稳定度安全系数的映射关系;使用激光扫描技术对山体的实貌图像进行处理获取山体的三维几何坐标,根据获得的三维几何坐标进行逐层三维建模得到三维模型,再对三维模型进行网格划分;将映射关系和致灾因子输入到三维模型进行数值模拟评估,得到评估结果。该方法能够有效预测发生滑坡的可能性,以便及时对危险区域进行防治。
Description
技术领域
本发明涉及滑坡灾害风险评估领域,特别是涉及一种喀斯特区域滑坡灾害风险动态评估方法。
背景技术
喀斯特作为特殊的地质地貌,其区域内多发的地质灾害,特别是滑坡灾害已成为影响"三农"乃至区域经济社会发展的重要因素。
目前,关于喀斯特区域的滑坡灾害的研究主要方向是建立灾害数据库,以实现动态监测和预测等预防措施。总的来说,关于喀斯特区域的滑坡的研究不多,特别是滑坡灾害风险评估方面的研究更是少之又少。
发明内容
本发明主要解决的技术问题是提供一种喀斯特区域滑坡灾害风险动态评估方法,能够有效预测发生滑坡的可能性,以便及时对危险区域进行防治。
为解决上述技术问题,本发明采用的一个技术方案是:提供一种喀斯特区域滑坡灾害风险动态评估方法,包括以下步骤:S1:选取降雨量、位移矢量和沉降量作为致灾因子;S2:对致灾因子、山体稳定度进行量化、归一化处理和权重确定;S3:以大量区域灾害统计及人工降雨侵蚀试验结果为训练样本,采用BP人工神经网络方法,给定BP网络层数及相应层之间的权重向量,采用改进BP网络计算方法对神经网络进行训练直至达到误差要求,获得致灾因子与山体稳定度安全系数的映射关系;S4:使用激光扫描技术对山体的实貌图像进行处理获取山体的三维几何坐标,根据获得的三维几何坐标进行逐层三维建模得到三维模型,再对三维模型进行网格划分;S5:将所述映射关系和所述致灾因子输入到三维模型进行数值模拟评估,得到评估结果。
其中,所述对三维模型进行网格划分之前,还包括:根据计算机计算速度确定单元网格的尺寸大小。
其中,所述步骤S2具体包括:S21:采用隶属度函数法对致灾因子进行量化得到量化数据;S22:以sigmoid函数为激活函数对所述量化数据进行归一化处理,将所述量化数据转化为[0.1-0.9]区间的数值;S23:将山体稳定度分类为稳定、基本稳定、不稳定三类,并以sigmoid函数为激活函数进行归一化处理,将山体稳定度转化为[0.1-0.9]区间的安全系数;S24:构建层次结构模型,所述层析结构模型包括由下至上的方案层、准则层和目标层,其中,目标层为山体稳定度,方案层为各致灾因子,采用改进层次分析计算方法得到各层次的权重向量。
本发明的有益效果是:区别于现有技术的情况,本发明通过选取降雨量、位移矢量和沉降量作为致灾因子,再利用大量区域灾害统计及人工降雨侵蚀试验结果得到致灾因子与山体稳定度安全系数的映射关系,并建立喀斯特地貌滑坡的三维模型,通过将映射关系和致灾因子输入到三维模型进行数值模拟评估,直观的展现出山体滑坡危险区域的分布,从而能够有效预测发生滑坡的可能性,以便及时对危险区域进行防治。
附图说明
图1是本发明实施例喀斯特区域滑坡灾害风险动态评估方法的流程图。
图2是图1中步骤S2的具体流程图。
图3是山体的实貌图像与三维模型的对比图。
图4是进行数值模拟评估后得到的山体结构面位移矢量及其弹塑性区域分布的模拟结果示意图。
图5是进行数值模拟评估后得到的山体沉降量的模拟结果示意图。
具体实施例
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参阅图1,是本发明实施例喀斯特区域滑坡灾害风险动态评估方法的流程图。本方法包括以下步骤:
S1:选取降雨量、位移矢量和沉降量作为致灾因子。
其中,确定致灾因子是开展滑坡灾害风险评估的第一步。确定致灾因子,建立评估体系要充分考虑喀斯特地貌的地理特征,本发明根据大量区域灾害统计及人工降雨侵蚀试验结果确定降雨量、位移矢量和沉降量的主要致灾因子。
S2:对致灾因子、山体稳定度进行量化、归一化处理和权重确定。
在本实施例中,如图2所示,步骤S2具体包括:
S21:采用隶属度函数法对致灾因子进行量化得到量化数据。致灾因子的定量指标直接根据实际值以衡量,将各定性指标分为3级,通过给各等级赋予边界值,然后通过线性差值的方法确定其隶属函数,完成定性指标的量化。在其他实施例中,定性指标的量化还可以采用德尔菲法、统计分析法和信息量法其中之一进行。
S22:以sigmoid函数为激活函数对量化数据进行归一化处理,将量化数据转化为[0.1-0.9]区间的数值。
S23:将山体稳定度分类为稳定、基本稳定、不稳定三类,并以sigmoid函数为激活函数进行归一化处理,将山体稳定度转化为[0.1-0.9]区间的安全系数。
S24:构建层次结构模型,所述层析结构模型包括由下至上的方案层、准则层和目标层,其中,目标层为山体稳定度,方案层为各致灾因子,采用改进层次分析计算方法得到各层次的权重向量。
S3:以大量区域灾害统计及人工降雨侵蚀试验结果为训练样本,采用BP人工神经网络方法,给定BP网络层数及相应层之间的权重向量,采用改进BP网络计算方法对神经网络进行训练直至达到误差要求,获得致灾因子与山体稳定度安全系数的映射关系。
S4:使用激光扫描技术对山体的实貌图像进行处理获取山体的三维几何坐标,根据获得的三维几何坐标进行逐层三维建模得到三维模型,再对三维模型进行网格划分。
其中,采用Rhino软件进行逐层三维建模,用Kubrix软件进行网格划分。在本实施例中,对三维模型进行网格划分时,还包括:根据计算机计算速度确定单元网格的尺寸大小。如图3所示,是山体的实貌图像与三维模型的对比图。图中,图(a)为山体的实貌图像,图(b)为三维模型,从图中可以看出,三维模型与三体的轮廓基本一致。
S5:将映射关系和致灾因子输入到三维模型进行数值模拟评估,得到评估结果。
其中,使用Itasca公司的数值计算软件3DEC(3 Dimension Distinct ElementCode)进行数值模拟评估,3DEC是一款基于离散单元法作为基本理论以描述离散介质力学行为的计算分析程序,并且可以进行静力学和动力学分析。模拟中需要的岩石的力学参数通过室内单轴压缩、剪切及抗拉试验获取,模拟中考虑岩体的风化,对于岩体强度采用强度折减系数进行折减。数值模拟评估结果可以直观的表现出岩体不同区域的力学特性以及随时间发展岩体的破坏情况。如图4所示,是山体结构面位移矢量的数值模拟结果在左边,山体结构面的弹塑性区域分布数值模拟结果在右边。如图5所示,是山体沉降量数值模拟结果。通过评估结果,可以预测出山体滑坡危险区域,从而有效预测发生滑坡的可能性。
通过上述方式,本发明通过选取降雨量、位移矢量和沉降量作为致灾因子,再利用大量区域灾害统计及人工降雨侵蚀试验结果得到致灾因子与山体稳定度安全系数的映射关系,并建立喀斯特地貌滑坡的三维模型,通过将映射关系和致灾因子输入到三维模型进行数值模拟评估,直观的展现出山体滑坡危险区域的分布,从而能够有效预测发生滑坡的可能性,以便及时对危险区域进行防治。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (3)
1.一种喀斯特区域滑坡灾害风险动态评估方法,其特征在于,包括以下步骤:
S1:选取降雨量、位移矢量和沉降量作为致灾因子;S2:对致灾因子、山体稳定度进行量化、归一化处理和权重确定;S3:以大量区域灾害统计及人工降雨侵蚀试验结果为训练样本,采用BP人工神经网络方法,给定BP网络层数及相应层之间的权重向量,采用改进BP网络计算方法对神经网络进行训练直至达到误差要求,获得致灾因子与山体稳定度安全系数的映射关系;S4:使用激光扫描技术对山体的实貌图像进行处理获取山体的三维几何坐标,根据获得的三维几何坐标进行逐层三维建模得到三维模型,再对三维模型进行网格划分;S5:将所述映射关系和所述致灾因子输入到三维模型进行数值模拟评估,得到评估结果。
2.根据权利要求1所述的方法,其特征在于,所述对三维模型进行网格划分之前,还包括:根据计算机计算速度确定单元网格的尺寸大小。
3.根据权利要求2所述的方法,其特征在于,所述步骤S2具体包括:S21:采用隶属度函数法对致灾因子进行量化得到量化数据;S22:以sigmoid函数为激活函数对所述量化数据进行归一化处理,将所述量化数据转化为[0.1-0.9]区间的数值;S23:将山体稳定度分类为稳定、基本稳定、不稳定三类,并以sigmoid函数为激活函数进行归一化处理,将山体稳定度转化为[0.1-0.9]区间的安全系数;S24:构建层次结构模型,所述层析结构模型包括由下至上的方案层、准则层和目标层,其中,目标层为山体稳定度,方案层为各致灾因子,采用改进层次分析计算方法得到各层次的权重向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710205591.XA CN106952000A (zh) | 2017-03-31 | 2017-03-31 | 一种喀斯特区域滑坡灾害风险动态评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710205591.XA CN106952000A (zh) | 2017-03-31 | 2017-03-31 | 一种喀斯特区域滑坡灾害风险动态评估方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106952000A true CN106952000A (zh) | 2017-07-14 |
Family
ID=59474068
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710205591.XA Pending CN106952000A (zh) | 2017-03-31 | 2017-03-31 | 一种喀斯特区域滑坡灾害风险动态评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106952000A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109584515A (zh) * | 2018-12-21 | 2019-04-05 | 武汉万集信息技术有限公司 | 山体灾害的预警方法、装置及可读存储介质 |
CN113012399A (zh) * | 2021-02-25 | 2021-06-22 | 中原工学院 | 一种降雨型滑坡预警方法及系统 |
CN113095007A (zh) * | 2021-04-07 | 2021-07-09 | 南京金创有色金属科技发展有限公司 | 一种安全阀长周期运行安全评估技术方法 |
CN113295851A (zh) * | 2021-05-25 | 2021-08-24 | 山东省地质矿产勘查开发局第六地质大队(山东省第六地质矿产勘查院) | 一种矿山地质环境预警装置及其预警方法 |
CN115345511A (zh) * | 2022-08-29 | 2022-11-15 | 中咨数据有限公司 | 一种公路走廊滑坡危险性动态评价方法、评价系统及设备 |
CN117058628A (zh) * | 2023-10-13 | 2023-11-14 | 东华理工大学南昌校区 | 一种基于图像识别的滑坡灾害监测方法、系统及存储介质 |
-
2017
- 2017-03-31 CN CN201710205591.XA patent/CN106952000A/zh active Pending
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109584515A (zh) * | 2018-12-21 | 2019-04-05 | 武汉万集信息技术有限公司 | 山体灾害的预警方法、装置及可读存储介质 |
CN109584515B (zh) * | 2018-12-21 | 2021-12-28 | 武汉万集信息技术有限公司 | 山体灾害的预警方法、装置及可读存储介质 |
CN113012399A (zh) * | 2021-02-25 | 2021-06-22 | 中原工学院 | 一种降雨型滑坡预警方法及系统 |
CN113012399B (zh) * | 2021-02-25 | 2022-06-14 | 中原工学院 | 一种降雨型滑坡预警方法及系统 |
CN113095007A (zh) * | 2021-04-07 | 2021-07-09 | 南京金创有色金属科技发展有限公司 | 一种安全阀长周期运行安全评估技术方法 |
CN113095007B (zh) * | 2021-04-07 | 2023-11-07 | 南京金创有色金属科技发展有限公司 | 一种安全阀长周期运行安全评估技术方法 |
CN113295851A (zh) * | 2021-05-25 | 2021-08-24 | 山东省地质矿产勘查开发局第六地质大队(山东省第六地质矿产勘查院) | 一种矿山地质环境预警装置及其预警方法 |
CN113295851B (zh) * | 2021-05-25 | 2023-04-28 | 山东省地质矿产勘查开发局第六地质大队(山东省第六地质矿产勘查院) | 一种矿山地质环境预警装置及其预警方法 |
CN115345511A (zh) * | 2022-08-29 | 2022-11-15 | 中咨数据有限公司 | 一种公路走廊滑坡危险性动态评价方法、评价系统及设备 |
CN115345511B (zh) * | 2022-08-29 | 2023-06-06 | 中咨数据有限公司 | 一种公路走廊滑坡危险性动态评价方法、评价系统及设备 |
CN117058628A (zh) * | 2023-10-13 | 2023-11-14 | 东华理工大学南昌校区 | 一种基于图像识别的滑坡灾害监测方法、系统及存储介质 |
CN117058628B (zh) * | 2023-10-13 | 2024-01-26 | 东华理工大学南昌校区 | 一种基于图像识别的滑坡灾害监测方法、系统及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106952000A (zh) | 一种喀斯特区域滑坡灾害风险动态评估方法 | |
Lin et al. | Large-scale ecological red line planning in urban agglomerations using a semi-automatic intelligent zoning method | |
CN105654414A (zh) | 基于开源体系架构和建筑物空间数据库的城市多灾害风险损失评测系统及方法 | |
Guo et al. | Simulation evaluation of urban low-carbon competitiveness of cities within Wuhan city circle in China | |
Zhang et al. | Risk-based estimate for operational safety in complex projects under uncertainty | |
Wang | Site Selection of Fire Stations in Cities Based on Geographic Information System and Fuzzy Analytic Hierarchy Process. | |
CN106909999A (zh) | 小型土石坝退役综合评价方法 | |
CN107480889A (zh) | 一种湿地生态安全预警等级评定方法 | |
CN105184427A (zh) | 一种对农田生态环境进行预警的方法及装置 | |
CN111968019A (zh) | 一种基于cnn与rnn集成的滑坡灾害易发性预测方法 | |
CN102646164A (zh) | 一种结合空间滤波的土地利用变化建模方法及其系统 | |
CN112163745A (zh) | 一种基于组合模型的区域安全风险分析方法及系统 | |
CN117217058A (zh) | 一种抽蓄电站水工结构安全监测方法、系统、设备及介质 | |
Quan | Smart design for sustainable neighborhood development | |
CN115659729A (zh) | 一种基于结构仿真计算的大坝安全监测分析方法及系统 | |
CN114723283A (zh) | 一种面向城市群的生态承载力遥感评估方法及装置 | |
CN114493052A (zh) | 多模型融合自适应新能源功率预测方法和系统 | |
CN114936751A (zh) | 一种基于gis和多线约束条件的停车场选址规划方法 | |
Qin et al. | Multi-scenario simulation for 2060 and driving factors of the eco-spatial carbon sink in the Beibu Gulf Urban Agglomeration, China | |
Hua-xi et al. | Study on spatial prediction and time forecast of landslide | |
Vukicevic et al. | GIS Based Multicrteria Analysis in Integration of SEA Process Into Planning, Case Study: South West Region, Republic of Ireland | |
Eichrodt | Development of a spatial trafficability evaluation system | |
Wang et al. | Geographic information system and system dynamics combination technique for municipal solid waste treatment station site selection | |
CN1329689C (zh) | 压力容器疲劳寿命安全预测方法 | |
Tidwell et al. | Decision support for integrated water-energy planning. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170714 |
|
WD01 | Invention patent application deemed withdrawn after publication |