CN105139585A - 一种土质边坡险情智能预警预报方法 - Google Patents
一种土质边坡险情智能预警预报方法 Download PDFInfo
- Publication number
- CN105139585A CN105139585A CN201510455712.7A CN201510455712A CN105139585A CN 105139585 A CN105139585 A CN 105139585A CN 201510455712 A CN201510455712 A CN 201510455712A CN 105139585 A CN105139585 A CN 105139585A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- slope
- soil
- time
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000006073 displacement reaction Methods 0.000 claims abstract description 48
- 238000012544 monitoring process Methods 0.000 claims abstract description 35
- 230000006378 damage Effects 0.000 claims abstract description 21
- 230000035699 permeability Effects 0.000 claims abstract description 21
- 238000009499 grossing Methods 0.000 claims description 20
- 230000008859 change Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000035772 mutation Effects 0.000 claims description 14
- 238000013528 artificial neural network Methods 0.000 claims description 12
- 238000013473 artificial intelligence Methods 0.000 claims description 11
- 238000011156 evaluation Methods 0.000 claims description 11
- 238000013277 forecasting method Methods 0.000 claims description 10
- 230000006870 function Effects 0.000 claims description 10
- 210000002569 neuron Anatomy 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000013139 quantization Methods 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 7
- 230000035508 accumulation Effects 0.000 claims description 6
- 238000009825 accumulation Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 238000006467 substitution reaction Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 4
- 239000002356 single layer Substances 0.000 claims description 4
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 3
- PMGQWSIVQFOFOQ-YKVZVUFRSA-N clemastine fumarate Chemical compound OC(=O)\C=C\C(O)=O.CN1CCC[C@@H]1CCO[C@@](C)(C=1C=CC(Cl)=CC=1)C1=CC=CC=C1 PMGQWSIVQFOFOQ-YKVZVUFRSA-N 0.000 claims description 3
- 239000013307 optical fiber Substances 0.000 claims description 3
- 238000000611 regression analysis Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 230000008595 infiltration Effects 0.000 claims description 2
- 238000001764 infiltration Methods 0.000 claims description 2
- 230000035515 penetration Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000003723 Smelting Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013523 data management Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
- G08B21/02—Alarms for ensuring the safety of persons
- G08B21/10—Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明涉及一种土质边坡险情智能预警预报方法,包括以下步骤:1)现场设置雨量监测站,埋置土壤渗透计和位移计;2)实时采集边坡的降雨工况、土壤渗透性和位移数据并传输到远程客户端;3)数据格式化处理;4)计算降雨工况、土壤渗透性、坡体位移等数据影响边坡安全系数的权重大小;5)确定边坡各时段的安全系数;6)计算各时段的降雨工况、边坡土壤渗透系数及变形、位移与安全系数的关系,预测土质边坡破坏发生时间;7)获取潜在破坏的边坡的临界破坏和失稳时间,预测灾害发生时间;8)综合步骤6)和7)结果,对边坡险情进行早期预警预报。与现有技术相比,本发明更合理定量化、更符合实际地进行土质边坡险情的预警预报。
Description
技术领域
本发明涉及土质边坡工程的灾害预警预报技术领域,尤其是涉及一种土质边坡险情智能预警预报方法。
背景技术
中国是世界上滑坡最严重的国家之一。仅在2006~2011年的6年间全国共发生滑坡157927起。近年来,中国每年因滑坡灾害造成的经济损失都在10亿美元以上。2009年全国仅特大型和大型滑坡灾害共发生16起,直接经济损失达到1.9亿元人民币;2011年死亡、失踪10人以上或直接经济损失1亿元以上的重大滑坡灾害发生7起,直接经济损失5.99亿元人民币。滑坡灾害不仅造成巨大的经济损失,还严重危害人民的生命安全。土质滑坡是滑坡中分布最为广泛、规模大、暴发频率高、突发性强和持续危害性较大的一类致灾体。我国面临的此类型滑坡灾害非常严峻,对滑坡灾害进行科学、合理的评价、预测、管理和控制是十分必要。据不完全统计,我国每年发生滑坡、泥石流、塌陷等各种地质灾害数万起,能够在灾害发生前做出预警的还不到十分之一,如何对这些地质灾害进行预警,切实保护人民的生命和财产安全,是我国当前面临的一个重要课题。
目前,国外内在对于土质边坡险情智能预警预报技术进行了一些相关研究:已建立了公路边坡变形的GPS一机多天线监测系统,实现了公路边坡灾害监测的高精度、自动化和全天候等目标;基于“3S”的大冶铁矿高陡边坡灾害信息处理系统实行全天候自动监测;在边坡灾害预警信息生成方面进行了研究;依据边坡的特性、使用条件和破坏的经济影响,基于评价边坡的预期变形失败和临界读取频率,构建边坡监测系统的可靠性图的方法,评价监控系统预警临界斜坡破坏的有效性和可靠性;运用径向基函数的降雨指标,建立日本泥石流和边坡破坏的早期预警系统;结合遗传算法和最小二乘法,提出预测土质和岩质边坡稳定性的进化多项式回归方法。目前,在灾害预警技术上,基于模糊层次分析原理,建立了公路路堑边坡失稳的预估模型,研制了相应的土质边坡失稳预测风险评估软件,并采用灰色关联原理,分析降雨条件下影响路基稳定性的因素的权重,并提出相应的预防措施。这些方法需要大量的专家经验,一般适用于岩质边坡的预测,而在土质边坡失稳预测还不成熟。不过,国内外这些研究成果还存在精度不高、实施较烦等不足。所以,研发一种高精度、高可靠性的土质边坡险情智能预警预报方法显得尤为紧迫。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种预警预报精度高、可靠性高、定量化和更加合理的土质边坡险情智能预警预报方法。
本发明的目的可以通过以下技术方案来实现:
一种土质边坡险情智能预警预报方法,其特征在于,包括以下步骤:
1)现场设置带有雨量计的雨量监测站,土体中埋置土壤渗透计和位移计;
2)使用雨量计、土壤渗透计和位移计分别实时采集边坡的降雨工况、坡体土壤渗透性和土体位移及变形数据并传输到数据采集站,数据采集站将接收的数据实时自动传输到远程客户端;
3)远程客户端接收数据,并对数据进行格式化处理;
4)采用层次分析灰色关联度方法,计算降雨工况、土壤渗透性、坡体位移、边坡坡度、地下水位和岩土体抗剪强度影响边坡安全系数的权重和灰色关联度;
5)采用弹塑性强度折减法,确定边坡各时段的安全系数,得到以时间为顺序的安全系数数列;
6)依据步骤3)、4)和5)的数据,建立反馈动力学神经网络,并采用人工智能模糊控制方法,计算各时段的降雨工况、坡体土壤渗透系数、土体位移及变形数据与安全系数的关系,预测土质边坡破坏的发生时间,当安全系数为1时,土质边坡发生破坏;
预测土质边坡破坏的发生时间计算式:t=f(Fs,k,u,ρ)。Fs=1时,t为土质边坡发生破坏的时间,式中,Fs为边坡的安全系数;t为土质边坡变形的计算时间;k为土壤渗透系数;u为坡体的变形或位移;ρ为降雨工况。
7)对步骤3)得到的数据采用突变理论预报法进行计算,获取潜在破坏的边坡的临界破坏时间和失稳时间;
8)综合步骤6)和7),在临界破坏时间或失稳时间之前对土质边坡的险情进行早期预警预报。
所述步骤1)中,所述的雨量计、土壤渗透仪和位移计带有自动数据A/D转换功能。
所述步骤2)中,数据采集站采用太阳能供电,所述雨量计、土壤渗透计和位移计与数据采集站之间通过光纤电缆传输数据,数据采集站与远程客户端之间通过物联网的虚拟IP技术以移动信号的方式传输数据。
所述步骤3)中,对边坡的降雨工况、坡体土壤渗透系数、土体位移及变形数据的格式化处理包括等间隔化处理、归一化处理和/或平滑化处理,所述的等间隔化处理具体为:原始监测数据序列x0(ti),各时段的间隔时间为Δti,用下式计算平均时间间隔:
其中ti为监测时间点,n为监测时间点个数,
用下式计算等间隔化后的新时间点t′i:
t′i=(i-1)Δt0+t1,i=1,2,…,n
用下式计算t′i对应的数值x(t′i):
x(t′i)=x0(ti)-Δx0(ti)
其中,
或判断t′i所处的区间,用区间的两端点值的两点插值计算等间隔化后对应的数值x(t′i),由下式计算:
上述等间隔化处理的数据包括降雨工况、坡体土壤渗透系数、土体位移及变形数据。
所述的归一化处理具体为:
在监测数据序列x1,x2....xi中找出最大值xmax和最小值xmin,将监测数据xi替换为(xi-xmin)/(xmax-xmin),
归一化后的所有的监测数据值为0到1之间。
所述的平滑化处理包括邻点中值平滑处理和多点平滑处理,其中邻点中值平滑处理为:
取两个相邻的离散数据之间的中点作为新的离散数据,并包括原来的始点和终点,新的离散数据计算式为:
多点平滑处理分为等权滑动平均和不等权滑动平均两类,其中等权滑动平均处理得到新的离散数据为:
采用五点不等权滑动平均的不等权滑动平均处理得到新的离散数据为:
式中,平滑化处理的数据x包括降雨工况、坡体土壤渗透系数、土体位移及变形数据。
所述步骤4)中,层次分析灰色关联方法具体包括以下步骤:
401)获得一个参考数列X0和n个比较数列X1,X2,…,Xn,数列元素为降雨工况ρ(t)、坡体土壤渗透性k(t)和土体位移u(t),计算n个比较数列各元素的关联系数,其中第i个比较数列第k个元素的关联系数ξi(k)为:
式中:η为分辨系数,0≤η≤1,第1层次最小差为不同k值计算的绝对差|X0(k)-Xi(k)|中的最小值,第2层次最小差为Δ1(min),Δ2(min),…,Δn(min)中最小值,第1层次最大差为按不同k值计算的绝对差|X0(k)-Xi(k)|中的最大值,第2层次最大差为Δ1(max),Δ2(max),…,Δn(max)中最大值;
402)用下式计算第i个比较数列Xi对参考数列X0的关联度ri:
403)用下式计算层次分析的灰色关联度Ri:
式中,w(k)为通过降雨条件下土质边坡稳定性的影响因素层次分析计算,获取的第k个元素的权重。
所述步骤5)具体包括:
501)采用折减系数FS对边坡的岩土体抗剪强度指标C和进行折减:
CF=C/FS
式中,C是岩土体的粘聚力,是岩土体的内摩擦角;
502)用折减后的虚拟抗剪强度指标CF和取代原来的岩土体抗剪强度指标C和使用ABAQUS软件进行有限元稳定性计算,获取边坡安全系数。
所述步骤6)中,反馈动力神经网络的建立具体包括以下步骤:
601)确定单层互连的网络神经元,建立反馈动力神经网络的初始络状态;
602)利用外界输入激发网络进入动态演变过程;
603)待网络稳定,获得网络输出;
604)设定人工智能模糊控制方法的规则来源:包括专家经验和知识、操作员的操作模式和自我学习;
605)采用单输入和单输出的定值控制时间响应,在各阶段分别用状态评估模糊控制法则和目标评估模糊控制法则获得输出,模糊控制器之输出输入关系为:E,CE→CU,
其中,状态评估模糊控制中,前件部变量为输出误差量化值E和误差变化量化值CE,后件部变量为控制器输出变化量CU,其中误差量化值E、误差变化量化值CE分别为:
E=Ke×ek
式中,Ke为误差量化因子,ek为第k次采样得到的误差值,ek=R-Y,R为输出设定值,Y为系统输出,Kee为误差变化量化因子,T为采样周期;
为了权衡系统的快速性和稳定性,系统引入自调节因子α用以自动调整误差以及误差变化量对控制量的权重,模糊控制器输出变化量CU为:
CU=αE+(1-α)EC
式中,α为自调节因子,|E|max取E绝对值的最大值,αmin取α的下限,αmax取α的上限,并且满足0≤αmin≤αmax≤1。
606)选择各阶段的特征点,根据各阶段控制对象反应确定该阶段的控制法则。
所述步骤7)具体包括以下步骤:
701)用位移时间序列(ti,ui)监测经归一化后的数据(xi,yi),采用最小二乘法,求取给定线型和已知数据点的拟合曲线;
702)取一阶导数p′(x),计算出曲线各点的斜率tanα,并计算α=89°时的x值,获得边坡失稳时间预测值:
xi=ti/tn
yi=ui/un
t=tn·x
式中,t为边坡失稳时间预测值,ti为监测时间序列,tn为最后一次边坡位移监测时间,ui为对应ti时刻边坡位移测值,un为对应tn时刻边坡位移测值
703)获得边坡失稳时间预测值上的边坡坡体位移序列X(0)(t),对其进行一次累加得到序列X(1)(t),一次累加计算公式为:
其中,i为原始数据的个数,X(0)(i)为非负的原始数据,将X(1)(t)用Taylor展式截取到5次项:
式中:待定系数A0、…、A5通过回归分析求得,t为累计监测时间,单位为天;
704)对上式求导还原、变量代换得到控制变量的灰色尖点突变模型其中T为状态变量,a、b为状态变量系数;
705)令得到平衡曲面方程:
T3+aT+b=0
706)将平衡曲面的折叠或尖拐点集在a-b平面的投影,得到分叉集方程:
4a3+27b2=0
在分叉集上的(a,b)点为系统的临界状态,对分叉集方程求解,得到相应的边坡系统失稳前后时间差为:
Δt与临界状态的历时之和即为边坡系统破坏时间。
与现有技术相比,本发明具有以下优点:
1)本发明采用物联网的虚拟IP技术远程实时自动采集和传输现场降雨工况和坡体位移的信息,具有实时性的特点。
2)通过后台客户端软件进行人工智能及非平衡动力学灾害预测理论的分析和处理,获取边坡的失稳时间,然后进行边坡险情的早期预警预报,从而提高了土质边坡失稳破坏灾害的预警预报精度和可靠度。
3)本发明从降雨工况(降雨强度和降雨量随时间的变化)、坡体的位移等参数及其与边坡安全系数的关系得出失稳时间,并且进一步对所得失稳时间进行处理,与现有的依据降雨量和滑坡间的关系等其他角度不同,具有动态性,能够获取更加精确的失稳时间,从而进一步提高土质边坡失稳破坏灾害的预警预报精度。
4)突变理论预报方法和人工智能模糊控制分别得到了灾害发生时间,前者可以预测临界失稳时间,后者可以改善长期预测精确性,二者结合使用,进一步提高土质边坡失稳破坏灾害的预警预报精度和可靠性。
附图说明
图1为本发明的流程示意图;
图2为DHNN模型结构;
图3为突变理论预报法中分叉集的示意图;
图4为突变理论预报法中跨越分叉集时状态变量T发生突跳的示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,一种土质边坡险情智能预警预报方法,包括以下步骤:
S1,现场设置带有自动数据转换功能的雨量监测站,埋置带有自动数据转换功能的土壤渗透计及位移计;
S2,雨量计及传感器分别实时采集边坡降雨的工况和坡体变形和位移数据通过光纤电缆传输到数据采集站,数据采集站以移动信号的方式通过物联网的虚拟IP技术将接收的监测数据实时自动传输到远程客户端。
S3,装设有C#语言开发的《实时现场降雨工况及位移现场监测数据管理系统》模块的远程客户端接收监测数据,并进行数据格式化处理,包括等间隔化处理、归一化处理和/或平滑化处理。
等间隔化处理具体为:
设非等间隔边坡的降雨工况、坡体土壤渗透性和土体位移及变形原始监测序列x0(ti),各时段的间隔时间为Δti,且各时段间隔不完全相等,则平均时间间隔为
等间隔化后的新时间点为
t′i=(i-1)Δt0+t1(i=1,2,…,n)
t′i对应的数值x(t′i)可由下式计算
x(t′i)=x0(ti)-Δx0(ti)(i=1,2,…,n)
其中,
或者在编程时,直接判断t′i所处的区间,用区间的两端点值的两点插值计算出等间隔化后相应的数值x(t′i),计算方法为
归一化处理具体为:
在所有监测的数据中找出最大的数xmax和最小的数xmin。然后,把所有监测的数据进行如下计算:
(xi-xmin)/(xmax-xmin)
式中,xi为降雨工况、坡体土壤渗透性和位移的原始数据。
这样,即可实现所有的监测数据都归一化为0到1之间。
平滑化处理具体为:
监测数据的平滑处理可以采用邻点中值平滑处理和多点平滑处理方法。邻点中值平滑处理方法的处理过程是取两个相邻的离散数据之间的中点作为新的离散数据,并包括原来的始点和终点。新的离散数据可按下式计算
多点平滑处理可分为等权和不等权两类,对于等权滑动平均,新的离散数据为
不等权滑动平均,依据窗口大小不同有五点、十一点、十五点、二十一点等不同计算公式,其中五点不等权滑动平均的计算公式为
S4,采用反馈动力学神经网络和人工智能模糊控制方法以及突变理论预报法分别对步骤S3处理过的数据进行计算,预测土质边坡破坏的发生时间及其临界破坏时间。
其中,采用反馈动力学神经网络和人工智能模糊控制方法预测土质边坡破坏的发生时间具体包括以下步骤:
S401,采用弹塑性强度折减法计算边坡安全系数,具体包括:
S4011,采用折减系数FS对边坡的岩土体抗剪强度指标C和进行折减:
CF=C/FS
式中,C是岩土体的粘聚力,是岩土体的内摩擦角;
S4012,用折减后的虚拟抗剪强度指标CF和取代原来的岩土体抗剪强度指标C和进行ABAQUS有限元稳定性计算,获取边坡安全系数。
S402,运用层次分析灰色关联度理论,确定降雨工况、土壤渗透性、坡体位移、边坡坡度、地下水位和岩土体抗剪强度影响边坡安全系数的权重大小,具体包括以下步骤:
S4021,对于一个参考数列X0和若干个比较数列X1,X2,…,Xn,数列元素为降雨工况ρ(t)、坡体土壤渗透性k(t)、位移u(t),计算Xi的第k个元素的关联系数ξi(k),i=1,2….n。
式中:η为分辨系数,在[0,1]取值,通常取为0.5;为2个层次的最小差,第1层次最小差是指在绝对差|X0(k)-Xi(k)|中按不同k值挑选其中最小者,第2层次最小差是在Δ1(min),Δ2(min),…,Δn(min)中挑选其中最小者;同理,定义2个层次的最大差
S4022,根据比较数列Xi对参考数列X0的k个关联系数ξi(k),用下式求出比较数列Xi对参考数列X0的关联度ri:
S4023,层次分析的灰色关联度Ri计算如下:
式中:w(k)为通过降雨条件下土质边坡稳定性的影响因素层次分析计算,获取的第k个元素的权重。
S403,建立反馈动力学神经网络。如图2所示,本发明采用离散型单层反馈神经网络(DHNN网络模型),原理如下:
1)DHNN网中的每个神经元都有相同的功能,其输出xj称为状态。所有神经元的状态的集合构成反馈网络的状态:
X=[x1,x2,…,xn]T
反馈网络的输入就是网络的状态初始值,表示为
X(0)=[x1(0),x2(0),…,xn(0)]T
反馈网络在外界输入激发下,从初始状态进入动态演变过程,变化规律为
xj=f(netj),j=1,2,3,…n
2)DHNN网的转移函数常采用符号函数
式中,净输入为
反馈网络稳定时,每一个神经元的状态都不再改变,此时的稳定状态就是网络的输出,表示为:
3)DHNN神经网络模型由单层互连的神经元ui(i=1,2,3,…n)组成。神经元没有自连接,即:wii=0;神经元与神经元之间的连接是对称的,即wij=wji。权值矩阵为:
4)DHNN反馈神经网络的工作方式分为:
①网络的异步(串行)工作方式:
网络运行时每次只有一个神经元进行状态的调整计算,其它神经元的状态均保持不变,即
②网络的同步(并行)工作方式:
网络的同步工作方式是一种并行方式,所有神经元同时调整状态,即
xj(t+1)=sgn[netj(t)]j=1,2,3…,n
S404,根据反馈动力学神经网络,采用人工智能模糊控制方法,计算各时段的降雨工况、边坡土壤渗透系数及变形、位移与安全系数的关系,预测土质边坡破坏的发生时间,控制规则具体为:
1)规则来源
模糊控制规则的取得方式:
①专家的经验和知识
藉由询问经验丰富的专家,获得系统的知识,并将知识改为if….then的型式,如此便可构成模糊控制规则。除此之外,为了获得最佳的系统性能,常还需要多次使用试误法,以修正模糊控制规则。
②操作员的操作模式
记录操作员将操作模式整理为if….then的型式,可构成一组控制规则。
③学习
通过模糊控制器的自我学习或自我组织,使模糊控制器能够根据设定的目标,增加或修改模糊控制规则。
2)规则型式
模糊控制规则的形式主要可分为二种:
①状态评估模糊控制规则
状态评估模糊控制规则型式如下:
Ri:ifx1isAi1andx2isAi2….andxnisAin
thenyisCi
其中,x1,x2,……,xn及y为语言变量或称为模糊变量,代表系统的态变量和控制变量;Ai1,Ai2,….,Ain及Ci为语言值,代表论域中的模糊集合。该形式还有另一种表示法,是将后件部改为系统状态变量的函数,其形式如下:
Ri:ifx1isAi1andx2isAi2….andxnisAin
theny=f1(x1,x2,…….,xn)
②目标评估模糊控制规则
目标评估模糊控制规则能够评估控制目标,并且预测未来控制信号,其形式如下:
Ri:if(UisCi→(xisA1andyisB1))thenUisCi
3)规则流程
实际应用模糊控制时,最初的问题是控制器的设计,即如何设计模糊控制法则。本发明的模糊控制法则设计概念如下:
采用单输入和单输出的定值控制时间响应,在各阶段分别用状态评估模糊控制法则和目标评估模糊控制法则获得输出,模糊控制器之输出输入关系为:E,CE→CU,
其中,状态评估模糊控制中,前件部变量为输出的误差E和在一个取样周期内E的变化量CE,后件部变量为控制器输出量U的变化量CU。设ek为第k次采样所得到的误差值,E为误差量化值,CE为误差变化量化值,分别为:
ek=R-Y
E=Ke×ek
式中,R表示设定值,Y表示系统输出,Ke为误差量化因子,Kee为误差变化量化因子,T为采样周期;
为了权衡系统的快速性和稳定性,系统引入自调节因子α用以自动调整误差以及误差变化量对控制量的权重,模糊控制器输出变化量CU表示为
CU=αE+(1-α)EC
式中,α为自调节因子,|E|max取E绝对值的最大值,αmin取α的下限,αmax取α的上限,并且满足0≤αmin≤αmax≤1。
由此可知,误差变化量CE是随输出Y的斜率的符号变号,当输出上升时,CE<0,下降时CE>0。
本发明所设计的模糊控制器之输出输入关系为:
E,CE→CU
由于由E与CE推论CU的构造中,CU与E的关系恰巧相当于积分关系U(t)=Ki∫E(t)dt,而CU与CE的关系相当于比例关系U(t)=KpE(t)的缘故,所以又称为FuzzyPI控制。
4)模糊控制
设计模糊控制规则时,是在所设想对控制对象各阶段的反应,记述采取哪一种控制比较好;首先选择各阶段的特征点,记录在模糊控制规则的前件部,然后思考在该点采取的动作,记录在模糊控制规则的后件部。
采用突变理论预报法预测土质边坡的临界破坏时间,具体包括S405和S406:
S405,用位移时间序列(ti,ui)监测经上述归一化变换后的数据(xi,yi),对其作曲线拟合分析,获得最佳指使曲线p(x)。曲线的拟合可采用多种线型,针对边坡的变形特点,本软件系统采用了多项试拟合曲线和指数拟合曲线。为求取给定线型和已知数据点的最佳拟合曲线,可采用最小二乘法。
S406,取一阶导数p′(x),计算出曲线各点的斜率tanα,并计算α=89°时的x值,获得边坡失稳时间预测值:
xi=ti/tn
yi=ui/un
t=tn·x
式中,t为边坡失稳时间预测值,ti为监测时间序列,tn为最后一次边坡位移监测时间,ui为对应ti时刻边坡位移测值,un为对应tn时刻边坡位移测值。
通过以上计算,即可获得对应边坡失稳的时间预测值。
S407,采用突变理论预报法潜在破坏的边坡的临界破坏时间,具体为:
设边坡坡体位移的时间序列为X(0)(t),其一次累加序列为X(1)(t)。一次累加生成的计算公式为
其中,i为原始数据的个数;X(1)(i)为一次累加生成数据,右上角括号中的数字表示累加次数;X(0)(i)为非负的原始数据。
因为任何单变量的连续函数总可由Taylor展式表示成幂级数的形式,所以,X(1)(t)可表示为
式中:A0,A1,A2,…,An为待定系数;t为累计监测时间,单位为d。
当式(403-1)截取到5次项时,可近似表示为
式中:待定系数A0、A1、A2、A3、A4、A5通过回归分析求得。
根据突变理论的分析,只要截取Taylor级数到4次方项就能保证尖点位置的正确性。于是,对式(407-2)求导还原后得
令a0=A1,a1=2A2,a2=3A3,a3=4A4,a4=5A5,则(503-3)式变为
X=a0+a1t+a2t2+a3t3+a4t4(407-4)
为了简化运算,对上式作变量代换,消去三次项,可令
t=Tt-q(407-5)
式中:
将(407-5)式代入式(403-4),得
X=b4Tt 4+b2Tt 2+b1Tt+b0(407-6)
上式中
对(407-6)再作进一步的变量代换,令
这里仅以b4>0的情况为例进行分析,将(403-7)代入(503-6)式得
式中:C对突变分析是无意义的,a、b分别为
式(407-8)即为以T为状态变量(间接的时间状态变量),a、b为控制变量的灰色尖点突变模型。
特别指出,上述变换的目的,除了要构建尖点突变标准形式外,重要的一点还在于找出控制系统状态变化的“重要参量”,而忽略次要变量。
根据突变理论,令得平衡曲面方程
T3+aT+b=0(407-10)
式(407-10)给出了状态变量T和控制变量a、b的关系,详见图3所示的曲面。这一平衡曲面的折叠或尖拐点集称为奇点集,它在a-b平面的投影称为分叉集(见图3)。分叉集方程为
4a3+27b2=0(407-11)
令D=4a3+27b2,由图3可知,分叉集为一半立方抛物线,在(0,0)处有一尖点。在分叉集上的(a,b)点(即a、b满足分叉集方程,D=0)对应于系统的不稳定状态(临界状态),系统可能由一个平衡状态突变到另一个平衡状态。同时分叉集又将控制变量平面分为两个区。在较大的那个区域内(见图3)中AOB以下部分,D<0时,围岩-支护系统也有可能失稳。判断不稳定点的标准是
由以上分析可知,边坡系统稳定性发生突变有以下两种情况
(1)D=0;
(2)D<0且3T2+a<0。
对边坡系统失稳预报而言,跨越分叉集的左支(b<0)的情况,这时对应的点是不稳定状态,T值发生跳跃。
现在来确定对应于左支分叉集上各点的状态变量值。在式(403-11)成立的条件下,当a=0时,方程(407-10)有三重零根,T1=T2=T3=0;当a<0时,方程(403-10)有三个实根,它们是
跨越分叉集时状态变量T发生突跳,如图4所示:
根据前述代换,相应的边坡系统失稳前后时间差为
利用上式便可得到边坡系统在临界状态(D=0)与破坏时的时差Δt,Δt与临界状态的历时之和即为边坡系统破坏时间。对于D<0,3T2+a<0的情况,系统状态在分叉集内,用上述方法确定的破坏时间一般应大于实际破坏时间。
进一步的研究表明,若b4<0,则边坡系统的状态正好与上述讨论的结果相反。
S5,根据上述土质边坡破坏的发生时间及其临界破坏时间,在临界破坏时间或失稳时间之前对土质边坡的险情进行早期预警预报。
Claims (10)
1.一种土质边坡险情智能预警预报方法,其特征在于,包括以下步骤:
1)现场设置带有雨量计的雨量监测站,土体中埋置土壤渗透计和位移计;
2)使用雨量计、土壤渗透计和位移计分别实时采集边坡的降雨工况、坡体土壤渗透性和土体位移及变形数据并传输到数据采集站,数据采集站将接收的数据实时自动传输到远程客户端;
3)远程客户端接收数据,并对数据进行格式化处理;
4)采用层次分析灰色关联度方法,计算降雨工况、土壤渗透性、坡体位移、边坡坡度、地下水位和岩土体抗剪强度影响边坡安全系数的权重和灰色关联度;
5)采用弹塑性强度折减法,确定边坡各时段的安全系数,得到以时间为顺序的安全系数数列;
6)依据步骤3)、4)和5)的数据,建立反馈动力学神经网络,并采用人工智能模糊控制方法,计算各时段的降雨工况、坡体土壤渗透系数、土体位移及变形数据与安全系数的关系,预测土质边坡破坏的发生时间,当安全系数为1时,土质边坡发生破坏;
7)对步骤3)得到的数据采用突变理论预报法进行计算,获取潜在破坏的边坡的临界破坏时间和失稳时间;
8)综合步骤6)和7),在临界破坏时间或失稳时间之前对土质边坡的险情进行早期预警预报。
2.根据权利要求1所述的一种土质边坡险情智能预警预报方法,其特征在于,所述步骤1)中,所述的雨量计、土壤渗透仪和位移计带有自动数据A/D转换功能。
3.根据权利要求1所述的一种土质边坡险情智能预警预报方法,其特征在于,所述步骤2)中,数据采集站采用太阳能供电,所述雨量计、土壤渗透计和位移计与数据采集站之间通过光纤电缆传输数据,数据采集站与远程客户端之间通过物联网的虚拟IP技术以移动信号的方式传输数据。
4.根据权利要求1所述的一种土质边坡险情智能预警预报方法,其特征在于,所述步骤3)中,对边坡的降雨工况、坡体土壤渗透系数、土体位移及变形数据的格式化处理包括等间隔化处理、归一化处理和/或平滑化处理,所述的等间隔化处理具体为:原始监测数据序列x0(ti),各时段的间隔时间为Δti,用下式计算平均时间间隔:
其中ti为监测时间点,n为监测时间点个数,
用下式计算等间隔化后的新时间点t′i:
t′i=(i-1)Δt0+t1,i=1,2,…,n
用下式计算t′i对应的数值x(t′i):
x(t′i)=x0(ti)-Δx0(ti)
其中,
或判断t′i所处的区间,用区间的两端点值的两点插值计算等间隔化后对应的数值x(t′i),由下式计算:
5.根据权利要求4所述的一种土质边坡险情智能预警预报方法,其特征在于,所述的归一化处理具体为:
在监测数据序列x1,x2....xi中找出最大值xmax和最小值xmin,将监测数据xi替换为(xi-xmin)/(xmax-xmin)。
6.根据权利要求4所述的一种土质边坡险情智能预警预报方法,其特征在于,所述的平滑化处理包括邻点中值平滑处理和多点平滑处理,其中邻点中值平滑处理为:
取两个相邻的离散数据之间的中点作为新的离散数据,并包括原来的始点和终点,新的离散数据计算式为:
多点平滑处理分为等权滑动平均和不等权滑动平均两类,其中等权滑动平均处理得到新的离散数据为:
采用五点不等权滑动平均的不等权滑动平均处理得到新的离散数据为:
7.根据权利要求1所述的一种土质边坡险情智能预警预报方法,其特征在于,所述步骤4)中,层次分析灰色关联方法具体包括以下步骤:
401)获得一个参考数列X0和n个比较数列X1,X2,…,Xn,数列元素为降雨工况ρ(t)、坡体土壤渗透性k(t)和土体位移u(t),计算n个比较数列各元素的关联系数,其中第i个比较数列第k个元素的关联系数ξi(k)为:
式中:η为分辨系数,0≤η≤1,第1层次最小差为不同k值计算的绝对差|X0(k)-Xi(k)|中的最小值,第2层次最小差为Δ1(min),Δ2(min),…,Δn(min)中最小值,第1层次最大差为按不同k值计算的绝对差|X0(k)-Xi(k)|中的最大值,第2层次最大差为Δ1(max),Δ2(max),…,Δn(max)中最大值;
402)用下式计算第i个比较数列Xi对参考数列X0的关联度ri:
403)用下式计算层次分析的灰色关联度Ri:
式中,w(k)为通过降雨条件下土质边坡稳定性的影响因素层次分析计算,获取的第k个元素的权重。
8.根据权利要求1所述的一种土质边坡险情智能预警预报方法,其特征在于,所述步骤5)具体包括:
501)采用折减系数FS对边坡的岩土体抗剪强度指标C和进行折减:
CF=C/FS
式中,C是岩土体的粘聚力,是岩土体的内摩擦角;
502)用折减后的虚拟抗剪强度指标CF和取代原来的岩土体抗剪强度指标C和使用ABAQUS软件进行有限元稳定性计算,获取边坡安全系数。
9.根据权利要求1所述的一种土质边坡险情智能预警预报方法,其特征在于,所述步骤6)中,反馈动力神经网络的建立具体包括以下步骤:
601)确定单层互连的网络神经元,建立反馈动力神经网络的初始络状态;
602)利用外界输入激发网络进入动态演变过程;
603)待网络稳定,获得网络输出;
604)设定人工智能模糊控制方法的规则来源:包括专家经验和知识、操作员的操作模式和自我学习;
605)采用单输入和单输出的定值控制时间响应,在各阶段分别用状态评估模糊控制法则和目标评估模糊控制法则获得输出,模糊控制器之输出输入关系为:E,CE→CU,
其中,状态评估模糊控制中,前件部变量为输出误差量化值E和误差变化量化值CE,后件部变量为控制器输出变化量CU,其中误差量化值E、误差变化量化值CE分别为:
E=Ke×ek
式中,Ke为误差量化因子,ek为第k次采样得到的误差值,ek=R-Y,R为输出设定值,Y为系统输出,Kee为误差变化量化因子,T为采样周期;
模糊控制器输出变化量CU为:
CU=αE+(1-α)EC
式中,α为自调节因子,|E|max取E绝对值的最大值,αmin取α的下限,αmax取α的上限,并且满足0≤αmin≤αmax≤1;
606)选择各阶段的特征点,根据各阶段控制对象反应确定该阶段的控制法则。
10.根据权利要求1所述的一种土质边坡险情智能预警预报方法,其特征在于,所述步骤7)具体包括以下步骤:
701)用位移时间序列(ti,ui)监测经归一化后的数据(xi,yi),采用最小二乘法,求取给定线型和已知数据点的拟合曲线;
702)取一阶导数p′(x),计算出曲线各点的斜率tanα,并计算α=89°时的x值,获得边坡失稳时间预测值:
xi=ti/tn
yi=ui/un
t=tn·x
式中,t为边坡失稳时间预测值,ti为监测时间序列,tn为最后一次边坡位移监测时间,ui为对应ti时刻边坡位移测值,un为对应tn时刻边坡位移测值;
703)获得边坡失稳时间预测值上的边坡坡体位移序列X(0)(t),对其进行一次累加得到序列X(1)(t),一次累加计算公式为:
其中,i为原始数据的个数,X(0)(i)为非负的原始数据,
将X(1)(t)用Taylor展式截取到5次项:
式中:待定系数A0、…、A5通过回归分析求得,t为累计监测时间,单位为天;
704)对上式求导还原、变量代换得到控制变量的灰色尖点突变模型其中T为状态变量,a、b为状态变量系数;
705)令得到平衡曲面方程:
T3+aT+b=0
706)将平衡曲面的折叠或尖拐点集在a-b平面的投影,得到分叉集方程:
4a3+27b2=0
在分叉集上的(a,b)点为系统的临界状态,对分叉集方程求解,得到相应的边坡系统失稳前后时间差为:
Δt与临界状态的历时之和即为边坡系统破坏时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510455712.7A CN105139585B (zh) | 2015-07-29 | 2015-07-29 | 一种土质边坡险情智能预警预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510455712.7A CN105139585B (zh) | 2015-07-29 | 2015-07-29 | 一种土质边坡险情智能预警预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105139585A true CN105139585A (zh) | 2015-12-09 |
CN105139585B CN105139585B (zh) | 2017-11-07 |
Family
ID=54724913
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510455712.7A Expired - Fee Related CN105139585B (zh) | 2015-07-29 | 2015-07-29 | 一种土质边坡险情智能预警预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105139585B (zh) |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529197A (zh) * | 2016-12-07 | 2017-03-22 | 中国地质大学(武汉) | 一种滑坡稳定性时变规律分析的方法 |
CN106599372A (zh) * | 2016-11-18 | 2017-04-26 | 中国电建集团成都勘测设计研究院有限公司 | 边坡变形阶段位移时间曲线切线角的确定方法 |
CN106645651A (zh) * | 2017-02-20 | 2017-05-10 | 长沙市玖车测控技术有限公司 | 水土流失监测与预警系统 |
CN106649910A (zh) * | 2016-08-31 | 2017-05-10 | 广西交通科学研究院 | 结合边坡浅层滑动面的包络图判断边坡稳定性的方法 |
CN106710148A (zh) * | 2016-12-28 | 2017-05-24 | 中国科学院深圳先进技术研究院 | 地质灾害预测方法及系统 |
CN107655444A (zh) * | 2017-09-08 | 2018-02-02 | 广西交通科学研究院有限公司 | 结合监测与有限元修正以确定边坡滑动面变动的方法 |
CN107831516A (zh) * | 2017-10-31 | 2018-03-23 | 河海大学 | 融合gnss和地面监测网获取大坝实时高精度位移的方法 |
CN108229726A (zh) * | 2017-12-15 | 2018-06-29 | 同济大学 | 一种碎块石边坡灾变险情预警预报方法 |
CN108256172A (zh) * | 2017-12-26 | 2018-07-06 | 同济大学 | 一种顶管下穿既有箱涵过程中险情预警预报方法 |
CN108921390A (zh) * | 2018-06-08 | 2018-11-30 | 同济大学 | 一种土体大变形灾害简化评估方法 |
CN109034619A (zh) * | 2018-07-26 | 2018-12-18 | 长江勘测规划设计研究有限责任公司 | 一种基于模糊综合层次分析的弃渣场安全评价方法 |
US20190051146A1 (en) * | 2017-08-09 | 2019-02-14 | Institute Of Mountain Hazards And Environment, Chinese Academy Of Sciences | Three-dimensional multi-point multi-index early warning method for risk at power grid tower in landslide section |
CN110018292A (zh) * | 2019-04-19 | 2019-07-16 | 中国地质调查局西安地质调查中心 | 基于水-吸应力-局部安全场的黄土滑坡精细预警技术 |
CN110261573A (zh) * | 2019-05-16 | 2019-09-20 | 同济大学 | 一种高位岩质滑坡稳定性动态数值评价方法 |
CN110309525A (zh) * | 2019-03-22 | 2019-10-08 | 北京北科安地科技发展有限公司 | 一种边坡空间变形及破坏趋势计算方法 |
CN110333336A (zh) * | 2019-07-05 | 2019-10-15 | 东北大学 | 一种降雨条件下监测土坡失稳预警系统及方法 |
CN110457853A (zh) * | 2019-08-20 | 2019-11-15 | 沙洲职业工学院 | 一种含结构面的拱坝坝肩岩体失稳判断方法 |
CN110595535A (zh) * | 2019-08-19 | 2019-12-20 | 湖南强智科技发展有限公司 | 监测方法、装置和存储介质 |
CN110847969A (zh) * | 2019-10-29 | 2020-02-28 | 中国电建集团华东勘测设计研究院有限公司 | 一种适用于岩体条件地下洞室群变形分级预警指标的确定方法 |
CN111144636A (zh) * | 2019-12-22 | 2020-05-12 | 同济大学 | 一种边坡变形预测方法 |
CN111504268A (zh) * | 2020-04-22 | 2020-08-07 | 深圳市地质局 | 一种土质边坡险情智能预警预报方法 |
CN111539568A (zh) * | 2020-04-22 | 2020-08-14 | 深圳市地质局 | 基于无人机和三维建模技术的安全监测系统及方法 |
CN112526104A (zh) * | 2020-11-06 | 2021-03-19 | 马鞍山矿山研究总院股份有限公司 | 一种边坡稳定性监测测预警方法、系统、介质 |
CN113111595A (zh) * | 2021-06-15 | 2021-07-13 | 西南交通大学 | 一种边坡智能预警方法和装置 |
CN113821949A (zh) * | 2021-09-08 | 2021-12-21 | 中电建电力检修工程有限公司 | 一种基于变形监测的岩质边坡稳定安全与可靠度预测方法 |
CN114330144A (zh) * | 2022-02-15 | 2022-04-12 | 西南交通大学 | 一种边坡危石预警方法、装置、设备及可读存储介质 |
CN114485788A (zh) * | 2022-01-12 | 2022-05-13 | 北京科技大学 | 基于倾斜与强振特征的边坡危岩体崩塌预警方法及装置 |
CN114495429A (zh) * | 2021-12-23 | 2022-05-13 | 北京中关村智连安全科学研究院有限公司 | 一种崩塌灾害早期预警方法及系统 |
CN114997548A (zh) * | 2022-08-05 | 2022-09-02 | 武九铁路客运专线湖北有限责任公司 | 一种基于遗传神经网络的动水地质渗流实时预警方法 |
CN116580542A (zh) * | 2023-07-14 | 2023-08-11 | 四川川核地质工程有限公司 | 一种洪水预警方法及系统 |
CN117094178A (zh) * | 2023-10-18 | 2023-11-21 | 矿冶科技集团有限公司 | 边坡稳定性预测方法、装置、电子设备及存储介质 |
CN117609742A (zh) * | 2024-01-24 | 2024-02-27 | 中建安装集团有限公司 | 一种用于实现智能化管理的边坡施工监管系统及方法 |
CN118376287A (zh) * | 2024-06-27 | 2024-07-23 | 朝阳县自然资源事务服务中心 | 一种国土土地资源动态监测系统及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006195650A (ja) * | 2005-01-12 | 2006-07-27 | Chuo Kaihatsu Kk | 斜面崩壊監視予測システム |
CN103198619A (zh) * | 2013-03-08 | 2013-07-10 | 青岛理工大学 | 岩体滑坡的动力监测预警方法 |
CN104179514A (zh) * | 2014-08-18 | 2014-12-03 | 同济大学 | 水下隧道破碎围岩突水预测与渗流控制的方法 |
CN104408899A (zh) * | 2014-11-05 | 2015-03-11 | 同济大学 | 一种山岭高速公路花岗岩残坡积土滑坡远程实时预警方法 |
-
2015
- 2015-07-29 CN CN201510455712.7A patent/CN105139585B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006195650A (ja) * | 2005-01-12 | 2006-07-27 | Chuo Kaihatsu Kk | 斜面崩壊監視予測システム |
CN103198619A (zh) * | 2013-03-08 | 2013-07-10 | 青岛理工大学 | 岩体滑坡的动力监测预警方法 |
CN104179514A (zh) * | 2014-08-18 | 2014-12-03 | 同济大学 | 水下隧道破碎围岩突水预测与渗流控制的方法 |
CN104408899A (zh) * | 2014-11-05 | 2015-03-11 | 同济大学 | 一种山岭高速公路花岗岩残坡积土滑坡远程实时预警方法 |
Non-Patent Citations (3)
Title |
---|
张丽,黄英,任磊,孙宏波: "边坡稳定影响因素及稳定分析方法探讨", 《云南水力发电》 * |
陈新民,罗国煜: "基于经验的边坡稳定性灰色系统分析与评价", 《岩土工程学报》 * |
韩传超,韩艳霞: "边坡稳定性分析方法研究及发展综述", 《能源环境保护》 * |
Cited By (49)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106649910A (zh) * | 2016-08-31 | 2017-05-10 | 广西交通科学研究院 | 结合边坡浅层滑动面的包络图判断边坡稳定性的方法 |
CN106599372A (zh) * | 2016-11-18 | 2017-04-26 | 中国电建集团成都勘测设计研究院有限公司 | 边坡变形阶段位移时间曲线切线角的确定方法 |
CN106529197A (zh) * | 2016-12-07 | 2017-03-22 | 中国地质大学(武汉) | 一种滑坡稳定性时变规律分析的方法 |
CN106529197B (zh) * | 2016-12-07 | 2019-03-22 | 中国地质大学(武汉) | 一种滑坡稳定性时变规律分析的方法 |
CN106710148A (zh) * | 2016-12-28 | 2017-05-24 | 中国科学院深圳先进技术研究院 | 地质灾害预测方法及系统 |
CN106645651A (zh) * | 2017-02-20 | 2017-05-10 | 长沙市玖车测控技术有限公司 | 水土流失监测与预警系统 |
US20190051146A1 (en) * | 2017-08-09 | 2019-02-14 | Institute Of Mountain Hazards And Environment, Chinese Academy Of Sciences | Three-dimensional multi-point multi-index early warning method for risk at power grid tower in landslide section |
CN107655444A (zh) * | 2017-09-08 | 2018-02-02 | 广西交通科学研究院有限公司 | 结合监测与有限元修正以确定边坡滑动面变动的方法 |
CN107831516A (zh) * | 2017-10-31 | 2018-03-23 | 河海大学 | 融合gnss和地面监测网获取大坝实时高精度位移的方法 |
CN107831516B (zh) * | 2017-10-31 | 2020-11-20 | 河海大学 | 融合gnss和地面监测网获取大坝实时高精度位移的方法 |
CN108229726A (zh) * | 2017-12-15 | 2018-06-29 | 同济大学 | 一种碎块石边坡灾变险情预警预报方法 |
CN108256172A (zh) * | 2017-12-26 | 2018-07-06 | 同济大学 | 一种顶管下穿既有箱涵过程中险情预警预报方法 |
CN108921390A (zh) * | 2018-06-08 | 2018-11-30 | 同济大学 | 一种土体大变形灾害简化评估方法 |
CN108921390B (zh) * | 2018-06-08 | 2022-04-05 | 同济大学 | 一种土体大变形灾害简化评估方法 |
CN109034619A (zh) * | 2018-07-26 | 2018-12-18 | 长江勘测规划设计研究有限责任公司 | 一种基于模糊综合层次分析的弃渣场安全评价方法 |
CN110309525A (zh) * | 2019-03-22 | 2019-10-08 | 北京北科安地科技发展有限公司 | 一种边坡空间变形及破坏趋势计算方法 |
CN110018292A (zh) * | 2019-04-19 | 2019-07-16 | 中国地质调查局西安地质调查中心 | 基于水-吸应力-局部安全场的黄土滑坡精细预警技术 |
CN110261573A (zh) * | 2019-05-16 | 2019-09-20 | 同济大学 | 一种高位岩质滑坡稳定性动态数值评价方法 |
CN110261573B (zh) * | 2019-05-16 | 2021-09-03 | 同济大学 | 一种高位岩质滑坡稳定性动态评价方法 |
CN110333336A (zh) * | 2019-07-05 | 2019-10-15 | 东北大学 | 一种降雨条件下监测土坡失稳预警系统及方法 |
CN110595535A (zh) * | 2019-08-19 | 2019-12-20 | 湖南强智科技发展有限公司 | 监测方法、装置和存储介质 |
WO2021031461A1 (zh) * | 2019-08-20 | 2021-02-25 | 沙洲职业工学院 | 一种含结构面的拱坝坝肩岩体失稳判断方法 |
CN110457853A (zh) * | 2019-08-20 | 2019-11-15 | 沙洲职业工学院 | 一种含结构面的拱坝坝肩岩体失稳判断方法 |
CN110457853B (zh) * | 2019-08-20 | 2023-04-18 | 沙洲职业工学院 | 一种含结构面的拱坝坝肩岩体失稳判断方法 |
CN110847969B (zh) * | 2019-10-29 | 2021-06-08 | 中国电建集团华东勘测设计研究院有限公司 | 一种适用于岩体条件地下洞室群变形分级预警指标的确定方法 |
CN110847969A (zh) * | 2019-10-29 | 2020-02-28 | 中国电建集团华东勘测设计研究院有限公司 | 一种适用于岩体条件地下洞室群变形分级预警指标的确定方法 |
CN111144636A (zh) * | 2019-12-22 | 2020-05-12 | 同济大学 | 一种边坡变形预测方法 |
CN111144636B (zh) * | 2019-12-22 | 2022-08-16 | 同济大学 | 一种边坡变形预测方法 |
CN111504268A (zh) * | 2020-04-22 | 2020-08-07 | 深圳市地质局 | 一种土质边坡险情智能预警预报方法 |
CN111539568A (zh) * | 2020-04-22 | 2020-08-14 | 深圳市地质局 | 基于无人机和三维建模技术的安全监测系统及方法 |
CN112526104A (zh) * | 2020-11-06 | 2021-03-19 | 马鞍山矿山研究总院股份有限公司 | 一种边坡稳定性监测测预警方法、系统、介质 |
CN113111595B (zh) * | 2021-06-15 | 2021-08-31 | 西南交通大学 | 一种边坡智能预警方法和装置 |
CN113111595A (zh) * | 2021-06-15 | 2021-07-13 | 西南交通大学 | 一种边坡智能预警方法和装置 |
CN113821949A (zh) * | 2021-09-08 | 2021-12-21 | 中电建电力检修工程有限公司 | 一种基于变形监测的岩质边坡稳定安全与可靠度预测方法 |
CN113821949B (zh) * | 2021-09-08 | 2023-12-26 | 中电建电力检修工程有限公司 | 一种基于变形监测的岩质边坡稳定安全与可靠度预测方法 |
CN114495429A (zh) * | 2021-12-23 | 2022-05-13 | 北京中关村智连安全科学研究院有限公司 | 一种崩塌灾害早期预警方法及系统 |
CN114485788B (zh) * | 2022-01-12 | 2022-10-11 | 北京科技大学 | 基于倾斜与强振特征的边坡危岩体崩塌预警方法及装置 |
CN114485788A (zh) * | 2022-01-12 | 2022-05-13 | 北京科技大学 | 基于倾斜与强振特征的边坡危岩体崩塌预警方法及装置 |
CN114330144A (zh) * | 2022-02-15 | 2022-04-12 | 西南交通大学 | 一种边坡危石预警方法、装置、设备及可读存储介质 |
CN114997548B (zh) * | 2022-08-05 | 2022-10-18 | 武九铁路客运专线湖北有限责任公司 | 一种基于遗传神经网络的动水地质渗流实时预警方法 |
CN114997548A (zh) * | 2022-08-05 | 2022-09-02 | 武九铁路客运专线湖北有限责任公司 | 一种基于遗传神经网络的动水地质渗流实时预警方法 |
CN116580542A (zh) * | 2023-07-14 | 2023-08-11 | 四川川核地质工程有限公司 | 一种洪水预警方法及系统 |
CN116580542B (zh) * | 2023-07-14 | 2023-09-15 | 四川川核地质工程有限公司 | 一种洪水预警方法及系统 |
CN117094178A (zh) * | 2023-10-18 | 2023-11-21 | 矿冶科技集团有限公司 | 边坡稳定性预测方法、装置、电子设备及存储介质 |
CN117094178B (zh) * | 2023-10-18 | 2023-12-29 | 矿冶科技集团有限公司 | 边坡稳定性预测方法、装置、电子设备及存储介质 |
CN117609742A (zh) * | 2024-01-24 | 2024-02-27 | 中建安装集团有限公司 | 一种用于实现智能化管理的边坡施工监管系统及方法 |
CN117609742B (zh) * | 2024-01-24 | 2024-03-26 | 中建安装集团有限公司 | 一种用于实现智能化管理的边坡施工监管系统及方法 |
CN118376287A (zh) * | 2024-06-27 | 2024-07-23 | 朝阳县自然资源事务服务中心 | 一种国土土地资源动态监测系统及方法 |
CN118376287B (zh) * | 2024-06-27 | 2024-08-30 | 朝阳县自然资源事务服务中心 | 一种国土土地资源动态监测系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105139585B (zh) | 2017-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105139585B (zh) | 一种土质边坡险情智能预警预报方法 | |
CN110008301B (zh) | 基于机器学习的区域性地质灾害易发性预测方法及装置 | |
CN103400044B (zh) | 一种改进的水环境安全评价分析方法 | |
CN104408899B (zh) | 一种山岭高速公路花岗岩残坡积土滑坡远程实时预警方法 | |
CN107765347A (zh) | 一种高斯过程回归和粒子滤波的短期风速预测方法 | |
CN115063020B (zh) | 基于风险监测融合的梯级水电站多维安全调度装置及方法 | |
CN113222296A (zh) | 一种基于数字孪生的防洪调度方法 | |
CN103793854A (zh) | 多重组合优化的架空输电线路运行风险信息化评估方法 | |
CN106372799B (zh) | 一种电网安全风险预测方法 | |
CN102346964A (zh) | 道路交通网络区域拥堵实时预测与智能管理系统 | |
CN110288144B (zh) | 一种基于效用理论和d-s证据推理的深基坑安全性评估方法 | |
CN104063612A (zh) | 一种隧道工程风险态势模糊评估方法及评估系统 | |
KR20140068436A (ko) | 시계열 예측 모델을 이용한 이상관측자료 탐지 방법 및 지하수위의 이상관측자료 탐지 방법 | |
CN103530818A (zh) | 一种基于brb系统的供水管网建模方法 | |
CN117113236B (zh) | 一种智慧城市监控系统及数据处理方法 | |
CN117870609A (zh) | 一种基于非完整拱效应的软岩隧道掌子面变形监测方法 | |
CN117191147A (zh) | 一种泄洪大坝水位监测预警方法及系统 | |
CN106759546B (zh) | 基于改进多变量灰色预测模型的基坑变形预测方法及装置 | |
CN115310536A (zh) | 基于神经网络和gcn深度学习模型的水库水位预测预警方法 | |
Zhang et al. | Real-time burst detection based on multiple features of pressure data | |
CN104462855A (zh) | 一种地下构筑物监控数据处理和分析的方法及设备 | |
Dong | Mine gas emission prediction based on Gaussian process model | |
CN101382474B (zh) | 一种桥梁结构安全的安全预警方法 | |
CN107704706A (zh) | 基于灰色理论的泥石流灾害易发性预测方法 | |
CN103353295A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171107 Termination date: 20210729 |