CN105930660B - 一种基于盲信息的堆积体滑坡安全评价方法 - Google Patents

一种基于盲信息的堆积体滑坡安全评价方法 Download PDF

Info

Publication number
CN105930660B
CN105930660B CN201610256635.7A CN201610256635A CN105930660B CN 105930660 B CN105930660 B CN 105930660B CN 201610256635 A CN201610256635 A CN 201610256635A CN 105930660 B CN105930660 B CN 105930660B
Authority
CN
China
Prior art keywords
blind information
blind
information
safety
expression formula
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.)
Expired - Fee Related
Application number
CN201610256635.7A
Other languages
English (en)
Other versions
CN105930660A (zh
Inventor
向志鹏
王环玲
徐卫亚
孟庆祥
杨兰兰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201610256635.7A priority Critical patent/CN105930660B/zh
Publication of CN105930660A publication Critical patent/CN105930660A/zh
Application granted granted Critical
Publication of CN105930660B publication Critical patent/CN105930660B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

本发明公开一种基于盲信息的堆积体滑坡安全评价方法,步骤包括:将粘聚力和内摩擦角定义为堆积体滑坡安全评价的盲信息参数;选取滑带土体样本进行直剪实验,获取各样本的盲信息数据,通过数据处理,分别得到粘聚力盲信息A和内摩擦角盲信息B的初步表达式;以梯形模糊数模拟、表征盲信息,利用α截集技术对梯形模糊数进行定量处理,分别得到粘聚力盲信息A、内摩擦角盲信息B的最终表达式;通过刚体极限平衡方法,构建滑坡安全评价的盲信息模型,得到安全系数的盲信息C的表达式C=A*B;根据安全系数盲信息C的表达式对堆积体滑坡的安全性进行评价。本发明将盲信息理论与刚体极限平衡方法相结合,可为工程技术人员提供更为全面的滑坡失稳信息。

Description

一种基于盲信息的堆积体滑坡安全评价方法
技术领域
本发明涉及一种堆积体滑坡安全评价方法,具体涉及一种基于盲信息的堆积体滑坡安全评价方法。
背景技术
在我国西南高山峡谷以及世界范围内分布着大量的堆积体滑坡,这类滑坡在人类或者自然活动中,灾害频发,给经济建设和人民正常生活带来严重危害,堆积体滑坡的安全性评价具有重要意义。
堆积体作为一种多元介质,一般由作为骨料的碎石、块石和黏土、砂土等充填成分组成,其力学特性复杂,结构具有强烈的不均匀性和随机性。相对于均质边坡,堆积体滑坡的力学参数的不确定性对滑坡的安全性评价影响特别突出。目前工程上对于滑坡的安全性评价主要采用确定性方法,确定性方法最终只能提供一个确定的安全系数,在实际的工程中,有可能出现在计算时被判定为稳定的边坡,实际中发生了失稳破坏的情况。
针对现有技术中存在的上述问题,本发明提出一种新的堆积体滑坡安全评价方法,该方法能够弥补传统确定性方法的缺陷,为工程技术人员提供更为全面的滑坡失稳信息。
发明内容
发明目的:本发明的目的在于提供一种基于盲信息的堆积体滑坡安全评价方法,该方法充分考虑了计算参数的不确定性,能够为工程技术人员提供更为全面的滑坡失稳信息。
技术方案:本发明所述的一种基于盲信息的堆积体滑坡安全评价方法,包括如下步骤:
(1)将粘聚力和内摩擦角定义为堆积体滑坡安全评价的盲信息参数;
(2)选取滑带土体样本进行直剪实验,获取各样本的盲信息数据,通过数据处理,分别得到粘聚力盲信息A和内摩擦角盲信息B的初步表达式;
(3)以梯形模糊数模拟、表征盲信息,利用α截集技术对梯形模糊数进行定量处理,分别得到粘聚力盲信息A、内摩擦角盲信息B的最终表达式;
(4)通过刚体极限平衡方法,构建滑坡安全评价的盲信息模型,得到安全系数的盲信息C的表达式C=A*B;
(5)根据安全系数盲信息C的表达式对堆积体滑坡的安全性进行评价。
具体的,步骤(2)中,数据处理的方法为:剔除错误数据;采用层次聚类分析法,对剩余数据进行数据区段的划分,以各数据区段的数据量占剩余数据总量的百分比作为该区段的可信度。
上述步骤(3)中,以梯形模糊数模拟、表征盲信息的方法为:计算各数据区段内数据的均值E(x),方差δ(x),取a=E(x)-k1δ(x)、d=E(x)+k1δ(x)分别作为梯形模糊数的上下限,取b=E(x)-k2δ(x)、c=E(x)+k2δ(x)分别作为梯形模糊数的最可能值。
上述步骤(4)中,安全系数盲信息C的表达式为:
hα(z)=fα(x).gα(y);
其中,为盲信息A经α截集技术处理后与α值相对应的区间数,fα(x)为盲信息A的可信度;为盲信息B经过α截集技术处理后与α值相对应的区间数,gα(y)为盲信息B的可信度;hα(z)为计算得到的盲信息C的可信度;*为盲信息的加减乘除运算。
上述步骤(5)中,根据边坡安全规范,得出天然工况下堆积体滑坡允许的安全系数,结合安全系数盲信息C的表达式,计算滑坡低于允许安全系数的概率,该概率表征滑坡的失稳安全性。
有益效果:与现有技术相比,本发明的优点在于:本发明通过将盲信息理论与刚体极限平衡方法相结合,应用于岩土工程中滑坡安全评价数值计算领域,一方面综合考虑了堆积体这种岩土体材料结构强烈的不均匀性和随机性的特点,能够充分反映堆积体这种多元介质材料参数的不确定性;另一方面弥补了采用传统确定性方法进行安全性评价时产生的安全系数单一的缺点,能够为工程技术人员提供更为全面的滑坡失稳信息,为滑坡安全性评价提供了新的思路,具有较强的工程意义。
附图说明
图1为本发明基于盲信息的堆积体滑坡安全评价方法的简要流程图;
图2为实施例中堆积体滑坡的典型剖面图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
如图1,本发明的一种基于盲信息的堆积体滑坡安全评价方法,包括如下步骤:
(1)将粘聚力和内摩擦角定义为堆积体滑坡安全评价的盲信息参数;
粘聚力和内摩擦角是影响堆积体滑坡安全性的两个主要指标。
(2)选取滑带土体样本进行直剪实验,获取各样本的盲信息数据,通过数据处理,分别得到粘聚力盲信息A和内摩擦角盲信息B的初步表达式;
数据处理包括:剔除错误数据;然后采用层次聚类分析法,对剩余数据进行数据区段的划分,以各数据区段的数据量占剩余数据总量的百分比作为该区段的可信度;
(3)以梯形模糊数模拟、表征盲信息,利用α截集技术对梯形模糊数进行定量处理,分别得到粘聚力盲信息A、内摩擦角盲信息B的最终表达式;
其中,以梯形模糊数模拟、表征盲信息的方法为:计算各数据区段内数据的均值E(x),方差δ(x),然后取a=E(x)-k1δ(x)、d=E(x)+k1δ(x)分别作为梯形模糊数的上下限,取b=E(x)-k2δ(x)、c=E(x)+k2δ(x)分别作为梯形模糊数的最可能值。
(4)通过刚体极限平衡方法,构建滑坡安全评价的盲信息模型,得到安全系数的盲信息C的表达式C=A*B;
安全系数盲信息C的表达式为:
hα(z)=fα(x).gα(y);
其中,表示盲信息A经α截集技术处理后与α值相对应的区间数,fα(x)表示盲信息A的可信度;表示盲信息B经过α截集技术处理后与α值相对应的区间数,gα(y)表示盲信息B的可信度;hα(z)表示计算后得到的盲信息C的可信度;*代表盲信息的加减乘除运算。
(5)根据安全系数盲信息C的表达式对堆积体滑坡的安全性进行评价。
具体的,根据边坡安全规范,得出天然工况下堆积体滑坡允许的安全系数,结合安全系数盲信息C的表达式,计算滑坡低于允许安全系数的概率,该概率表征滑坡的失稳安全性。
本发明通过将盲信息理论与刚体极限平衡方法相结合,应用于岩土工程中滑坡安全评价数值计算领域,一方面综合考虑了堆积体这种岩土体材料结构强烈的不均匀性和随机性的特点,能够充分反映堆积体这种多元介质材料参数的不确定性;另一方面弥补了采用传统确定性方法进行安全性评价时产生的安全系数单一的缺点,能够为工程技术人员提供更为全面的滑坡失稳信息,为滑坡安全性评价提供了新的思路,具有较强的工程意义。
实施例
下面以某水电站坝前堆积体滑坡为例,对本发明的基于盲信息的堆积体滑坡安全评价方法进行说明。
某水电站坝前堆积体滑坡由砂土构成,充填碎石,块石,充填物为变质砂岩,板岩。地质勘探表明,堆积体稳定性控制层为滑带土层,失稳破坏模式主要为沿着底滑面的整体滑移;滑坡典型剖面图如图2。
依据本发明方法,该滑坡的安全性评价步骤为:
(1)滑坡安全评价盲信息参数的定义:将影响堆积体滑坡安全性的两个主要指标粘聚力和内摩擦角定义为盲信息。
(2)堆积体力学参数实验数据的处理:
在现场钻孔选取不同部位的滑带土体样本,进行室内直剪实验,记录各样本对应的粘聚力和内摩擦角数据,本次实验总共进行了35组滑带土的剪切实验,剔除实验数据有错误的5组,剩余的30组实验数据为有效数据;采用层次聚类分析法,分别对两个盲信息的有效数据进行数据区段的划分,以各数据区段内数据量占总数据量的百分比作为该区段的可信度,得到盲信息的初步表达式,如下:
①内摩擦角的盲信息初步表达式为:
可以看出,内摩擦角在14.8~18.2kPa之间的可信度为0.15,在19.6~27.2kPa之间的可信度为0.65,在28.0~31.2kPa之间的可信度为0.20。
②粘聚力的盲信息初步表达式:
可以看出,粘聚力在20.4°~25.3°之间的可信度为0.20,在26.3°~35.2°之间的可信度为0.55,在33.0°~36.5°之间的可信度为0.25。
(3)梯形模糊参数的计算:以梯形模糊数分别模拟、表征内摩擦角和粘聚力,利用α截集技术对梯形模糊数进行定量处理,得到盲信息的最终表达式。结合工程实际情况,本算例中取α=0.7的截集水平进行计算,得到内摩擦角和粘聚力的盲信息最终表达式如下:
①内摩擦角的盲信息最终表达式:
②粘聚力的盲信息最终表达式:
(4)盲信息的运算及滑坡安全系数盲信息表达式的求解:按照盲信息的运算规则计算安全系数的盲信息表达式。为了能够尽可能的减少计算量,对计算结果进行适当的合并,得到安全系数的盲信息表达式。对于本实施例的整体滑移模式,安全系数的盲信息表达式为:
由上式可知,安全系数在0.82-1.09之间的可能性为12%,在1.10-1.20之间的可能性为10%,其余依次类推。
(5)滑坡失稳安全性评价:根据边坡安全规范,在天然工况下边坡允许的安全系数为1.25,因此可以计算此实施例中安全系数低于1.25的概率为:34%,存在一定的失稳风险。
该结果不仅充分考虑了堆积体材料的不均匀性,而且这样安全系数的表达式比较全面地反映了在各个不同的取值区间的可能性,避免了由于参数取值单一而造成的计算安全系数时的偏差,能够为工程提供全面的滑坡失稳参考信息,对后续滑坡预警以及治理具有重要意义。

Claims (1)

1.一种基于盲信息的堆积体滑坡安全评价方法,其特征在于,包括如下步骤:
(1)将粘聚力和内摩擦角定义为堆积体滑坡安全评价的盲信息参数;
(2)选取滑带土体样本进行直剪实验,获取各样本的盲信息数据,通过数据处理,分别得到粘聚力盲信息A和内摩擦角盲信息B的初步表达式;其中,所述数据处理的方法为:剔除错误数据;采用层次聚类分析法,对剩余数据进行数据区段的划分,以各数据区段的数据量占剩余数据总量的百分比作为该区段的可信度;
(3)以梯形模糊数模拟、表征盲信息,利用α截集技术对梯形模糊数进行定量处理,分别得到粘聚力盲信息A、内摩擦角盲信息B的最终表达式;其中,所述以梯形模糊数模拟、表征盲信息的方法为:计算各数据区段内数据的均值E(x),方差δ(x),取a=E(x)-k1δ(x)、d=E(x)+k1δ(x)分别作为梯形模糊数的上下限,取b=E(x)-k2δ(x)、c=E(x)+k2δ(x)分别作为梯形模糊数的最可能值;
(4)通过刚体极限平衡方法,构建滑坡安全评价的盲信息模型,得到安全系数的盲信息C的表达式C=A*B;
hα(z)=fα(x).gα(y);
上式中,为盲信息A经α截集技术处理后与α值相对应的区间数,fα(x)为盲信息A的可信度;为盲信息B经过α截集技术处理后与α值相对应的区间数,gα(y)为盲信息B的可信度;hα(z)为计算得到的盲信息C的可信度;*为盲信息的加减乘除运算;
(5)根据安全系数盲信息C的表达式对堆积体滑坡的安全性进行评价,根据边坡安全规范,得出天然工况下堆积体滑坡允许的安全系数,结合安全系数盲信息C的表达式,计算滑坡低于允许安全系数的概率,该概率表征滑坡的失稳安全性。
CN201610256635.7A 2016-04-22 2016-04-22 一种基于盲信息的堆积体滑坡安全评价方法 Expired - Fee Related CN105930660B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610256635.7A CN105930660B (zh) 2016-04-22 2016-04-22 一种基于盲信息的堆积体滑坡安全评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610256635.7A CN105930660B (zh) 2016-04-22 2016-04-22 一种基于盲信息的堆积体滑坡安全评价方法

Publications (2)

Publication Number Publication Date
CN105930660A CN105930660A (zh) 2016-09-07
CN105930660B true CN105930660B (zh) 2018-07-24

Family

ID=56838904

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610256635.7A Expired - Fee Related CN105930660B (zh) 2016-04-22 2016-04-22 一种基于盲信息的堆积体滑坡安全评价方法

Country Status (1)

Country Link
CN (1) CN105930660B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108133115B (zh) * 2018-01-12 2019-11-08 河北工业大学 基于数值模拟及极限平衡计算的滑坡危险性评价方法
CN111191393B (zh) * 2019-12-20 2022-02-11 河海大学 基于离散余弦变换的水动力型滑坡非确定性评价方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103077306A (zh) * 2012-12-31 2013-05-01 河海大学 一种基于赫斯特指数的边坡安全评价方法
CN103268420A (zh) * 2013-05-24 2013-08-28 河海大学 一种岩石高边坡的危险性评价方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103077306A (zh) * 2012-12-31 2013-05-01 河海大学 一种基于赫斯特指数的边坡安全评价方法
CN103268420A (zh) * 2013-05-24 2013-08-28 河海大学 一种岩石高边坡的危险性评价方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于盲数理论和响应面法的边坡稳定性分析;邱金伟等;《计算力学学报》;20141231;第31卷;第98-100页 *
基于盲数理论的边坡安全稳定分析研究;赵志峰等;《岩土力学》;20071130;第28卷(第11期);摘要、第2403-2404页 *
考虑计算参数不确定性的边坡稳定分析;赵志峰等;《第一届中国水利水电岩土力学与工程学术讨论会论文集》;20061231;摘要、第520页 *

Also Published As

Publication number Publication date
CN105930660A (zh) 2016-09-07

Similar Documents

Publication Publication Date Title
Sun et al. Dynamic load prediction of tunnel boring machine (TBM) based on heterogeneous in-situ data
WO2019042483A2 (zh) 一种tbm在掘岩体状态实时感知系统和方法
AU2015345707B2 (en) Method for determining earth surface interpenetrated crack distribution and air leakage characteristics in shallow burial coal mining
CN104458309B (zh) 一种用于物理模拟实验中的相似材料配比确定方法
CN101806224B (zh) 井下煤层气抽采产能预测方法
CN104863577B (zh) 利用地震纵波传播时间预测地层孔隙压力的方法
CN104408323A (zh) 一种基于多源信息融合的采场顶板离层水害超前预报方法
CN106124263A (zh) 一种地下矿山采空区相似材料3d模型的制作方法
CN110298107B (zh) 一种基于增量叠加的工作面冲击危险性评价方法
Azocar Investigating the mesh dependency and upscaling of 3D grain-based models for the simulation of brittle fracture processes in low-porosity crystalline rock
Zhimin et al. Predicting the height of water-flow fractured zone during coal mining under the Xiaolangdi Reservoir
Karaman et al. A comparative assessment of rock mass deformation modulus
CN105046080B (zh) 一种岩体质量评价方法
CN104636532A (zh) 一种煤矿瓦斯抽采钻孔封孔深度和长度的确定方法
He et al. A model to estimate the height of the water-conducting fracture zone for longwall panels in western China
CN106295214A (zh) 一种矿井瓦斯浓度预警方法
CN105930660B (zh) 一种基于盲信息的堆积体滑坡安全评价方法
Grenon et al. Analysis of a large rock slope failure on the east wall of the LAB chrysotile mine in Canada: back analysis, impact of water infilling and mining activity
Boah et al. Decline curve analysis and production forecast studies for oil well performance prediction: a case study of reservoir X
CN111666624A (zh) 一种无支护条件下巨跨地下洞室围岩稳定性评价方法
CN108333331A (zh) 小窑采空区浅埋巷道上覆及侧壁岩土稳定性评价方法
Yang et al. Probabilistic model of disc-cutter wear in TBM construction: A case study of Chaoer to Xiliao water conveyance tunnel in China
CN112184033B (zh) 一种碳酸盐岩岩层断层封闭性评价方法
CN114091260A (zh) 一种考虑断裂保护层隔水性的覆岩离层突水风险评价方法
CN106021717A (zh) 一种基于神经网络的地铁开挖引起的地表下沉分析方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180724

CF01 Termination of patent right due to non-payment of annual fee