CN111695177A - 基于位移场监测数据的大坝力学参数随机反演方法及系统 - Google Patents

基于位移场监测数据的大坝力学参数随机反演方法及系统 Download PDF

Info

Publication number
CN111695177A
CN111695177A CN202010407591.XA CN202010407591A CN111695177A CN 111695177 A CN111695177 A CN 111695177A CN 202010407591 A CN202010407591 A CN 202010407591A CN 111695177 A CN111695177 A CN 111695177A
Authority
CN
China
Prior art keywords
dam
module
displacement
monitoring data
inversion
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
Application number
CN202010407591.XA
Other languages
English (en)
Other versions
CN111695177B (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 CN202010407591.XA priority Critical patent/CN111695177B/zh
Publication of CN111695177A publication Critical patent/CN111695177A/zh
Application granted granted Critical
Publication of CN111695177B publication Critical patent/CN111695177B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Civil Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Architecture (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本发明公开一种基于位移场监测数据的大坝力学参数随机反演方法及系统,本发明的随机反演方法通过联合空间多测点的位移监测数据,给出了一阶和二阶摄动意义下的参数随机反演准则,能够对大坝真实的位移量和观测误差进行分离;本发明的随机反演系统包括数据库模块,结构信息模块,随机反演模块和信息展示模块:数据库模块用于采集、汇总大坝基础信息,地理空间信息,监测信息;结构信息模块用于存放工程的几何模型信息和有限元模型信息;随机反演模块用于对大坝力学参数的均值和标准差进行识别;信息展示模块具有工程模型展示功能、监控信息浏览功能和反演评估结果输出功能。本发明能够对大坝的真实运行状态和观测误差进行识别。

Description

基于位移场监测数据的大坝力学参数随机反演方法及系统
技术领域
本发明属于大坝安全领域,具体是一种基于位移场监测数据的大坝力学参数随机反演方法及系统。
背景技术
当前,世界范围内已建成大量的高坝工程,如中国的锦屏一级混凝土拱坝,三峡混凝土重力坝,瑞士的大迪克桑斯重力坝,以及美国的胡佛重力拱坝等。这些高坝工程在防洪、发电、航运、灌溉和供水等方面发挥着至关重要的作用,产生了非常显著的社会和经济效益。大坝在产生巨大工程效益的同时,其工作性态也会随着服役时间的增长而不断变化,会出现一些隐患。部分工程中出现的坝基失稳、坝体开裂、渗透破坏等问题,给相关国家带来了惨重的灾害和巨大的损失。奥地利科恩布莱恩、美国德沃夏克、前苏联萨扬舒申斯克等高坝都发生过严重开裂漏水,修补加固费用巨大;法国马尔帕塞拱坝因坝肩失稳发生溃决,给生命、财产带来巨大损失;中国也有高混凝土坝发生严重裂缝、高压水劈裂等影响安全的实例。因此,对大坝进行合理有效的检查、监控、评价、维护,对确保大坝安全和公共安全具有十分重要的意义。
大坝安全监测数据中含有大量丰富的信息,是大坝结构状态的综合表现,受多种因素的叠加影响。依据实测资料对大坝原型进行反分析,可对大坝前期的设计施工方案进行检验,对当下的工作性态进行评价,对未来的发展趋势进行预测;从而为及时评估和发现大坝的异常迹象提供充分依据,在此基础上制定恰当的水库控制运行计划和大坝的维护管理措施来确保大坝的运行安全,在发生险情时及时发布警报以避免事故的发生。
当前,对大坝结构的反演分析一般是基于各种确定性的数学模型,反演的力学参数某种程度上只能说是较好地反映实际工程情况的“等效参数”。作为挡水建筑物,大坝从建设到运行期间经历复杂内部外部条件变化,本质是一个受多种因素影响的不确定性系统。大坝监测信息的不确定性,主要由本身材料力学参数引起的响应量的摄动和观测误差这两部分组成,而目前的相关分析,没有对观测误差进行识别和剔除。因此,做好监测资料的分析工作,从观测信息中分离出有用的信息,找出影响坝体结构性态的主要因素,才能准确掌握大坝的安全运行状态。
发明内容
本发明的目的是针对现有技术存在的问题,提供一种基于位移场监测数据的大坝力学参数随机反演方法及系统,能对混凝土坝自动化监测信息进行实时处理,对大坝的真实安全状态进行反馈
为实现上述目的,本发明采用的技术方案是:
一种基于位移场监测数据的大坝力学参数随机反演方法,包括以下步骤:
获取大坝各测点的监测数据,对所述监测数据进行预处理;
根据预处理后的监测数据,建立各测点大坝变形的统计模型;
利用所述统计模型,分离大坝位移监测数据中的水压分量;
调用大坝的有限元模型计算对应水荷载下的大坝位移场;
采用一阶摄动随机反演方法或二阶摄动随机反演方法,并调用智能优化算法计算大坝力学参数的平均值和方差;
输出反演结果。
具体地,分离大坝位移监测数据中的水压分量的方法为:在统计模型中将大坝径向位移δ按照成因分为水压分量、温度分量和时效分量,在位移δ中扣除温度分量和时效分量后,得到水压分量如下式所示:
δH=δ-δTθ (1)
其中,δH为大坝径向位移δ中的水压分量;δT为大坝径向位移δ中的温度分量;δθ为大坝径向位移δ中的时效分量。
具体地,采用一阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
Figure BDA0002491966660000021
式中E为数学期望,将式(6)代入式(4)中,可得,
Figure BDA0002491966660000022
定义目标函数为:
Figure BDA0002491966660000023
一阶摄动随机反演的准则为目标函数取得最小值,即在方差最小的意义下,结合优化算法对参数均值
Figure BDA0002491966660000024
进行识别,
Figure BDA0002491966660000025
为X取均值下的位移计算值;
所述随机变量X的方差反演方法为:
对于某一观测时刻k,k=1,2,···,M,令
Figure BDA0002491966660000031
式中,
Figure BDA0002491966660000032
为N×p矩阵,则式(5)可写为:
Figure BDA0002491966660000033
式中,δ(X,Yk)为N×1列阵,{εXk}为p×1列阵,
Figure BDA0002491966660000034
式(10)变换得,
Figure BDA0002491966660000035
将上式两端分别乘以
Figure BDA0002491966660000036
得,
Figure BDA0002491966660000037
变换得,
Figure BDA0002491966660000038
随机变量的方差为:
Figure BDA0002491966660000039
将式(12)的计算结果代入式(10),可得对于某一观测时刻k,由随机变量引起的位移摄动为:
Figure BDA00024919666600000310
式中,{εδXk}为N×1列阵;对于某时刻k,观测误差为:
fk}={εδk}-{εδXk} (16)
联立上述公式(2)至(16),即可计算各测点由随机变量引起的位移方差var(δX)统计值,以及由观测工作产生的位移方差var(δf)统计值。
可选地,采用二阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
将δ(X,Y)在均值处进行包含二次项的Taylor展开,在不考虑随机变量相关条件下,有:
Figure BDA0002491966660000041
对上式取均值,有:
Figure BDA0002491966660000042
对某一观测时刻k,k=1,2,···,M,令:
Figure BDA0002491966660000043
则式(17)可写为:
Figure BDA0002491966660000044
此时,目标函数为:
Figure BDA0002491966660000045
式中
Figure BDA0002491966660000046
Ajk为矩阵Ak的第j行,j=1,2,···,N,即Ajk为1×p列阵;
根据目标函数,即在方差最小的意义下,结合优化算法对二阶意义下的参数均值
Figure BDA0002491966660000047
进行反演;此时
Figure BDA0002491966660000048
是未知的,采用一阶意义下的估计作为估计初值,进行均值
Figure BDA0002491966660000049
计算;
所述随机变量X的方差反演方法为:
式(20)可改写为:
Figure BDA00024919666600000410
式中
Figure BDA00024919666600000411
为N×1列阵;
定义目标函数Lk为:
Figure BDA00024919666600000412
式中Bjk为矩阵Bk的第j行,j=1,2,···,N,即Bjk为1×p列阵;Cjk为矩阵Ck的第j行,j=1,2,···,N,即Cjk为1×p列阵;
在目标函数Lk最小的意义下,对{εXk}进行求解,此时{εXk}满足方程组:
Figure BDA0002491966660000051
上式方程组展开后得:
Figure BDA0002491966660000052
式中Ajik是Ajk的第i个数,i=1,2,···,p;
上式可简化为:
Figure BDA0002491966660000053
此时式中Ak,Bk,Ck为N×1列阵;
上式(26)展开得:
Figure BDA0002491966660000054
根据卡尔丹公式求得一元三次方程的三个根,选目标函数最小的实根作为最终解;
此外对于多参数反演问题,调用优化函数直接对4次多项式目标函数进行寻优求解{εXk};
进一步可得随机变量的方差为:
Figure BDA0002491966660000055
根据上式求出
Figure BDA0002491966660000056
后代入式(21),反复迭代,最终达到
Figure BDA0002491966660000057
Figure BDA0002491966660000058
时结束迭代;此时求出
Figure BDA0002491966660000059
Figure BDA00024919666600000510
即为反演的参数和方差。
与上述反演方法相对应的,本发明还提供了一种基于位移场监测数据的大坝力学参数随机反演系统,包括数据库模块,结构信息模块,随机反演模块和信息展示模块;所述数据库模块用于存储大坝的基础信息数据和各类监测数据;所述结构信息模块用于存放大坝工程的结构模型信息;所述随机反演模块用于对大坝力学参数的均值和方差进行计算;所述信息展示模块用于展示大坝的模型、监测数据和反演结果。
具体地,所述数据库模块包括大坝基础信息模块、地理空间信息模块和监测信息模块;
所述大坝基础信息模块用于记载有工程概况、水库大坝设计、施工、运行相关资料;
所述地理空间信息模块用于记载水文和地质情况统计资料和工程地质勘察试验资料;
所述监测信息模块用于记载各类监测设备的原始测值,主要包括:环境量监测数据、变形监测数据、渗流监测数据和压力(应力)监测数据。
具体地,所述结构模型信息包括工程几何模型信息和有限元模型信息。
具体地,所述随机反演模块包括数据前处理模块,智能优化算法模块,一阶摄动随机反演模块和二阶摄动随机反演模块;
所述监测数据前处理模块,用于建立测点变形统计模型,分离得到测点变形的水压分量;
所述智能优化算法模块包括粒子群算法、灰狼算法、鲸鱼优化算法、蚁群算法、人工鱼群算法等智能寻优算法;可通过GUI及按钮触发群智能寻优算法选择框进行选择;
所述一阶摄动随机反演模块用于对一阶摄动意义下力学参数的均值和方差进行识别;
所述二阶摄动随机反演模块用于进一步对二阶摄动意义下力学参数的均值和方差进行识别。
具体地,所述信息展示模块包括模型轻量化展示、监控信息浏览和评估结果输出功能;其中,借助于BIMFACE轻量化模型转换功能,将所建立的大坝及监测测点三维信息模型在Web端进行展示,同时,可通过浏览器查看大坝监测测点信息和实时评估的结果信息。
与现有技术相比,本发明的有益效果是:
(1)本发明提出的随机反演方法考虑了大坝结构力学性质的不确定性以及观测误差的不确定性,联合空间多测点位移,能够对大坝结构力学参数进行准确识别;
(2)本发明通过随机反演方法能够将大坝的真实位移量和观测误差进行分离,提高数据的精确度和评估结果的可信度;
(3)本发明提出的随机反演系统能对大坝的自动化监测数据进行全链条处理,对大坝的真实安全状态进行实时反馈,从而协助技术人员高效地完成大坝安全监测及管理工作,具有重要的社会经济效益。
附图说明
图1为本发明一种基于位移场监测数据的大坝力学参数随机反演方法的流程图;
图2为本发明一种基于位移场监测数据的大坝力学参数随机反演系统的结构框图;
图3为本发明实施例中某碾压混凝土坝11#坝段有限元模型图;
图4为本发明实施例中大坝各测点变形的水压分量分离结果示意图;
图5为本发明实施例中垂线PL11-1测点位移反馈结果示意图。
具体实施方式
下面将结合本发明中的附图,对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动条件下所获得的所有其它实施例,都属于本发明保护的范围。
如图1所示,本实施例提供了一种基于位移场监测数据的大坝力学参数随机反演方法,包括以下步骤:
获取大坝各测点的监测数据,对所述监测数据进行预处理;
根据预处理后的监测数据,建立各测点大坝变形的统计模型;
利用所述统计模型,分离大坝位移监测数据中的水压分量;
调用大坝的有限元模型计算对应水荷载下的大坝位移场;
采用一阶摄动随机反演方法或二阶摄动随机反演方法,并调用智能优化算法计算大坝力学参数的平均值和方差;
输出反演结果。
具体地,分离大坝位移监测数据中的水压分量的方法为:在统计模型中将大坝径向位移δ按照成因分为水压分量、温度分量和时效分量,在位移δ中扣除温度分量和时效分量后,得到水压分量如下式所示:
δH=δ-δTθ (1)
其中,δH为大坝径向位移δ中的水压分量;δT为大坝径向位移δ中的温度分量;δθ为大坝径向位移δ中的时效分量。
具体地,采用一阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
Figure BDA0002491966660000071
式中E为数学期望,将式(6)代入式(4)中,可得,
Figure BDA0002491966660000081
定义目标函数为:
Figure BDA0002491966660000082
一阶摄动随机反演的准则为目标函数取得最小值,即在方差最小的意义下,结合优化算法对参数均值
Figure BDA0002491966660000083
进行识别,
Figure BDA0002491966660000084
为X取均值下的位移计算值;
所述随机变量X的方差反演方法为:
对于某一观测时刻k,k=1,2,···,M,令
Figure BDA0002491966660000085
式中,
Figure BDA0002491966660000086
为N×p矩阵,则式(5)可写为:
Figure BDA0002491966660000087
式中,δ(X,Yk)为N×1列阵,{εXk}为p×1列阵,
Figure BDA0002491966660000088
式(10)变换得,
Figure BDA0002491966660000089
将上式两端分别乘以
Figure BDA00024919666600000810
得,
Figure BDA00024919666600000811
变换得,
Figure BDA00024919666600000812
随机变量的方差为:
Figure BDA00024919666600000813
将式(12)的计算结果代入式(10),可得对于某一观测时刻k,由随机变量引起的位移摄动为:
Figure BDA00024919666600000814
式中,{εδXk}为N×1列阵;对于某时刻k,观测误差为:
fk}={εδk}-{εδXk} (16)
联立上述公式(2)至(16),即可计算各测点由随机变量引起的位移方差var(δX)统计值,以及由观测工作产生的位移方差var(δf)统计值。
可选地,采用二阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
将δ(X,Y)在均值处进行包含二次项的Taylor展开,在不考虑随机变量相关条件下,有:
Figure BDA0002491966660000091
对上式取均值,有:
Figure BDA0002491966660000092
对某一观测时刻k,k=1,2,···,M,令:
Figure BDA0002491966660000093
则式(17)可写为:
Figure BDA0002491966660000094
此时,目标函数为:
Figure BDA0002491966660000095
式中
Figure BDA0002491966660000096
Ajk为矩阵Ak的第j行,j=1,2,···,N,即Ajk为1×p列阵;
根据目标函数,即在方差最小的意义下,结合优化算法对二阶意义下的参数均值
Figure BDA0002491966660000097
进行反演;此时
Figure BDA0002491966660000098
是未知的,采用一阶意义下的估计作为估计初值,进行均值
Figure BDA0002491966660000099
计算;
所述随机变量X的方差反演方法为:
式(20)可改写为:
Figure BDA00024919666600000910
式中
Figure BDA0002491966660000101
为N×1列阵;
定义目标函数Lk为:
Figure BDA0002491966660000102
式中Bjk为矩阵Bk的第j行,j=1,2,···,N,即Bjk为1×p列阵;Cjk为矩阵Ck的第j行,j=1,2,···,N,即Cjk为1×p列阵;
在目标函数Lk最小的意义下,对{εXk}进行求解,此时{εXk}满足方程组:
Figure BDA0002491966660000103
上式方程组展开后得:
Figure BDA0002491966660000104
式中Ajik是Ajk的第i个数,i=1,2,···,p;
上式可简化为:
Figure BDA0002491966660000105
此时式中Ak,Bk,Ck为N×1列阵;
上式(26)展开得:
Figure BDA0002491966660000106
根据卡尔丹公式求得一元三次方程的三个根,选目标函数最小的实根作为最终解;
此外对于多参数反演问题,调用优化函数直接对4次多项式目标函数进行寻优求解{εXk};
进一步可得随机变量的方差为:
Figure BDA0002491966660000111
根据上式求出
Figure BDA0002491966660000112
后代入式(21),反复迭代,最终达到
Figure BDA0002491966660000113
Figure BDA0002491966660000114
时结束迭代;此时求出
Figure BDA0002491966660000115
Figure BDA0002491966660000116
即为反演的参数和方差。
与上述反演方法相对应的,如图2所示,本实施例还提供了一种基于位移场监测数据的大坝力学参数随机反演系统,包括数据库模块,结构信息模块,随机反演模块和信息展示模块;
所述数据库模块包括大坝基础信息模块、地理空间信息模块和监测信息模块;所述大坝基础信息模块用于记载有工程概况、水库大坝设计、施工、运行相关资料;所述地理空间信息模块用于记载水文和地质情况统计资料和工程地质勘察试验资料;所述监测信息模块用于记载各类监测设备的原始测值,主要包括:环境量监测数据、变形监测数据、渗流监测数据和压力(应力)监测数据。
所述结构模型信息包括工程几何模型信息和有限元模型信息。
所述随机反演模块包括数据前处理模块,智能优化算法模块,一阶摄动随机反演模块和二阶摄动随机反演模块;所述监测数据前处理模块,用于建立测点变形统计模型,分离得到测点变形的水压分量;所述智能优化算法模块包括粒子群算法、灰狼算法、鲸鱼优化算法、蚁群算法、人工鱼群算法等智能寻优算法;可通过GUI及按钮触发群智能寻优算法选择框进行选择;所述一阶摄动随机反演模块用于对一阶摄动意义下力学参数的均值和方差进行识别;所述二阶摄动随机反演模块用于进一步对二阶摄动意义下力学参数的均值和方差进行识别。
所述信息展示模块包括模型轻量化展示、监控信息浏览和评估结果输出功能;其中,借助于BIMFACE轻量化模型转换功能,将所建立的大坝及监测测点三维信息模型在Web端进行展示,同时,可通过浏览器查看大坝监测测点信息和实时评估的结果信息。
下面以一个工程实例来说明本发明:
在中国的西南地区建有一碾压混凝土坝工程。该工程前期建设中,于2007年年底浇筑至坝顶382m高程,其中建基面高程216.43m,坝高165.57m,坝顶宽度14m。坝体上游面折坡点位于270m高程,以下坡度约为1:0.25,下游面坡度约为1:0.70。
所述碾压混凝土重力坝工程的各类工程数据存放于数据库模块,可对相关数据执行查询、输出、处理、管理等工作。工程几何模型信息和有限元模型信息存放于结构信息模块。
如想了解11#坝段坝基的运行情况,操作步骤如下:
步骤(1),输入指令,通过结构信息模块查看并调用该坝段的结构信息、监测仪器布置信息和有限元模型信息;通过查询可知在11#坝段埋设有4个垂线用于监测水位移,坝段有限元模型图及对应测点位置示意图如图3所示;其中在151.5m高程埋设了1个倒垂线(IP-11),在379.2,342,270m高程埋设了3个正垂线(PL11-1,PL11-2,PL11-3,其中PL11-2是一线多垂,包括PL11-2-1和PL11-2-2)。
输入指令,从数据库中选取调用11#坝段对应得到监测资料,包括水位等环境量,各测点对应的实测值序列;本实施案例中选取观测序列从2010年6月28日至2013年4月15日,共59组测值。
步骤(2),数据前处理模块中对11坝段5个测点的监测数据进行预处理,建立各测点大坝变形的统计模型,并分离得到水压分量,其结果如图4所示。
步骤(3),调用结构信息中的有限元计算模型。
步骤(4),点击一阶或二阶摄动随机反演模块,本实施案例选择一阶摄动随机反演模块执行一阶摄动随机反演分析。
步骤(5),点击智能优化算法模块,本实施案例选择灰狼算法;输入坝体及基岩待反演参数取值范围,设置为坝体弹性模量Ec∈[35GPa,55GPa],坝基弹性模量Er∈[20GPa,45GPa];终止条件设置为迭代次数为100次。
步骤(6),再次点击一阶摄动随机反演模块中的触发按钮进行计算,程序进入迭代过程。
步骤(7),点击信息展示模块评估结果输出选择框,查看计算结果;计算结束后输出优化反演结果为:坝体弹性模量均值为45.53GPa,标准差为5.87GPa,变异系数为12.9%;坝基弹性模量均值为31.53GPa,标准差为6.12GPa,变异系数为19.4%。各测点处由随机力学参数引起的位移摄动std(δ′X)和观测误差std(δf′)的分离统计结果如表1所示。
表1各测点处位移摄动量的分离统计结果
测点 IP11 PL11-3 PL11-2-2 PL11-2-1 PL11-1
std(δ′<sub>Y</sub>)(mm) 013 028 047 061 081
std(δ′<sub>f</sub>)(mm) 008 009 007 009 007
同时,如图5所示,本实施例为位移反馈结果提供一个概率输出,即大坝测点位移的置信区间。如95%的置信区间可定义为
Figure BDA0002491966660000131
Figure BDA0002491966660000132
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。

Claims (8)

1.一种基于位移场监测数据的大坝力学参数随机反演方法,其特征在于,包括以下步骤:
获取大坝各测点的监测数据,对所述监测数据进行预处理;
根据预处理后的监测数据,建立各测点大坝变形的统计模型;
利用所述统计模型,分离大坝位移监测数据中的水压分量;
调用大坝的有限元模型计算对应水荷载下的大坝位移场;
采用一阶摄动随机反演方法或二阶摄动随机反演方法,并调用智能优化算法计算大坝力学参数的平均值和方差;
输出反演结果。
2.根据权利要求1所述的一种基于位移场监测数据的大坝力学参数随机反演方法,其特征在于,分离大坝位移监测数据中的水压分量的方法为:在统计模型中将大坝径向位移δ按照成因分为水压分量、温度分量和时效分量,在位移δ中扣除温度分量和时效分量后,得到水压分量如下式所示:
δH=δ-δTθ (1)
其中,δH为大坝径向位移δ中的水压分量;δT为大坝径向位移δ中的温度分量;δθ为大坝径向位移δ中的时效分量。
3.根据权利要求1所述的一种基于位移场监测数据的大坝力学参数随机反演方法,其特征在于,采用一阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
假设X为具有固定均值的随机变量,则表示为:
Figure FDA0002491966650000011
式中,X={x1,x2,…xp}T,p为随机变量的个数;
Figure FDA0002491966650000012
为X的均值,
Figure FDA0002491966650000013
εX为随机变量摄动值;
位移量δ(X,Y)是关于X,Y的函数,Y为确定性参数,Y={y1,y2,…yq}T,q为确定性参数的个数,若观测位移δ为非平稳随机过程变量,具有确定的趋势项,则δ表示为:
Figure FDA0002491966650000014
式中,
Figure FDA0002491966650000015
为δ的均值;εδ为位移摄动量,由随机变量摄动引起的位移摄动εδX和观测误差εf两部分组成;假定观测误差εf相互独立,数学期望为0且服从正态分布;
设有N个位移观测点,每个观测点有M个测值,则第j个测点的位移方差统计表达式为:
Figure FDA0002491966650000021
式中,var为方差算子;δ(X,Yjk)为第j个测点第k次观测量,j=1,2,···,N,k=1,2,···,M;令δjk=δ(X,Yjk),
Figure FDA0002491966650000022
Figure FDA0002491966650000023
为δjk的均值;
对于一阶摄动随机反演,将δ(X,Y)在均值处进行Taylor展开,并将高于一次的项舍去,则有:
Figure FDA0002491966650000024
对上式取均值,有:
Figure FDA0002491966650000025
式中E为数学期望,将式(6)代入式(4)中,可得,
Figure FDA0002491966650000026
定义目标函数为:
Figure FDA0002491966650000027
一阶摄动随机反演的准则为目标函数取得最小值,即在方差最小的意义下,结合优化算法对参数均值
Figure FDA0002491966650000028
进行识别,
Figure FDA0002491966650000029
为X取均值下的位移计算值;
所述随机变量X的方差反演方法为:
对于某一观测时刻k,k=1,2,···,M,令
Figure FDA00024919666500000210
式中,
Figure FDA00024919666500000211
为N×p矩阵,则式(5)可写为:
Figure FDA00024919666500000212
式中,δ(X,Yk)为N×1列阵,{εXk}为p×1列阵,
Figure FDA00024919666500000213
式(10)变换得,
Figure FDA00024919666500000214
将上式两端分别乘以
Figure FDA0002491966650000031
得,
Figure FDA0002491966650000032
变换得,
Figure FDA0002491966650000033
随机变量的方差为:
Figure FDA0002491966650000034
将式(12)的计算结果代入式(10),可得对于某一观测时刻k,由随机变量引起的位移摄动为:
Figure FDA0002491966650000035
式中,{εδXk}为N×1列阵;对于某时刻k,观测误差为:
fk}={εδk}-{εδXk} (16)
联立上述公式(2)至(16),即可计算各测点由随机变量引起的位移方差var(δX)统计值,以及由观测工作产生的位移方差var(δf)统计值。
4.根据权利要求3所述的一种基于位移场监测数据的大坝力学参数随机反演方法,其特征在于,采用二阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
将δ(X,Y)在均值处进行包含二次项的Taylor展开,在不考虑随机变量相关条件下,有:
Figure FDA0002491966650000036
对上式取均值,有:
Figure FDA0002491966650000037
对某一观测时刻k,k=1,2,···,M,令:
Figure FDA0002491966650000041
则式(17)可写为:
Figure FDA0002491966650000042
此时,目标函数为:
Figure FDA0002491966650000043
式中
Figure FDA0002491966650000044
Ajk为矩阵Ak的第j行,j=1,2,···,N,即Ajk为1×p列阵;
根据目标函数,即在方差最小的意义下,结合优化算法对二阶意义下的参数均值
Figure FDA0002491966650000045
进行反演;此时
Figure FDA0002491966650000046
是未知的,采用一阶意义下的估计作为估计初值,进行均值
Figure FDA0002491966650000047
计算;
所述随机变量X的方差反演方法为:
式(20)可改写为:
Figure FDA0002491966650000048
式中
Figure FDA0002491966650000049
为N×1列阵;
定义目标函数Lk为:
Figure FDA00024919666500000410
式中Bjk为矩阵Bk的第j行,j=1,2,···,N,即Bjk为1×p列阵;Cjk为矩阵Ck的第j行,j=1,2,···,N,即Cjk为1×p列阵;
在目标函数Lk最小的意义下,对{εXk}进行求解,此时{εXk}满足方程组:
Figure FDA00024919666500000411
上式方程组展开后得:
Figure FDA0002491966650000051
式中Ajik是Ajk的第i个数,i=1,2,···,p;
上式可简化为:
Figure FDA0002491966650000052
此时式中Ak,Bk,Ck为N×1列阵;
上式(26)展开得:
Figure FDA0002491966650000053
根据卡尔丹公式求得一元三次方程的三个根,选目标函数最小的实根作为最终解;
此外对于多参数反演问题,调用优化函数直接对4次多项式目标函数进行寻优求解{εXk};
进一步可得随机变量的方差为:
Figure FDA0002491966650000054
根据上式求出
Figure FDA0002491966650000055
后代入式(21),反复迭代,最终达到
Figure FDA0002491966650000056
Figure FDA0002491966650000057
时结束迭代;此时求出
Figure FDA0002491966650000058
Figure FDA0002491966650000059
即为反演的参数和方差。
5.一种基于位移场监测数据的大坝力学参数随机反演系统,其特征在于,包括数据库模块,结构信息模块,随机反演模块和信息展示模块;所述数据库模块用于存储大坝的基础信息数据和各类监测数据;所述结构信息模块用于存放大坝工程的结构模型信息;所述随机反演模块用于对大坝力学参数的均值和方差进行计算;所述信息展示模块用于展示大坝的模型、监测数据和反演结果。
6.根据权利要求5所述的一种基于位移场监测数据的大坝力学参数随机反演系统,其特征在于,所述数据库模块包括大坝基础信息模块、地理空间信息模块和监测信息模块;
所述大坝基础信息模块用于记载有工程概况、水库大坝设计、施工、运行相关资料;
所述地理空间信息模块用于记载水文和地质情况统计资料和工程地质勘察试验资料;
所述监测信息模块用于记载各类监测设备的原始测值,主要包括:环境量监测数据、变形监测数据、渗流监测数据和压力监测数据。
7.根据权利要求5所述的一种基于位移场监测数据的大坝力学参数随机反演系统,其特征在于,所述结构模型信息包括工程几何模型信息和有限元模型信息。
8.根据权利要求5所述的一种基于位移场监测数据的大坝力学参数随机反演系统,其特征在于,所述随机反演模块包括数据前处理模块,智能优化算法模块,一阶摄动随机反演模块和二阶摄动随机反演模块;
所述监测数据前处理模块,用于建立测点变形统计模型,分离得到测点变形的水压分量;
所述智能优化算法模块包括粒子群算法、灰狼算法、鲸鱼优化算法、蚁群算法、人工鱼群算法中的一个或多个智能寻优算法;
所述一阶摄动随机反演模块用于对一阶摄动意义下力学参数的均值和方差进行识别;
所述二阶摄动随机反演模块用于进一步对二阶摄动意义下力学参数的均值和方差进行识别。
CN202010407591.XA 2020-05-14 2020-05-14 基于位移场监测数据的大坝力学参数随机反演方法及系统 Active CN111695177B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010407591.XA CN111695177B (zh) 2020-05-14 2020-05-14 基于位移场监测数据的大坝力学参数随机反演方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010407591.XA CN111695177B (zh) 2020-05-14 2020-05-14 基于位移场监测数据的大坝力学参数随机反演方法及系统

Publications (2)

Publication Number Publication Date
CN111695177A true CN111695177A (zh) 2020-09-22
CN111695177B CN111695177B (zh) 2022-10-14

Family

ID=72477662

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010407591.XA Active CN111695177B (zh) 2020-05-14 2020-05-14 基于位移场监测数据的大坝力学参数随机反演方法及系统

Country Status (1)

Country Link
CN (1) CN111695177B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112287608A (zh) * 2020-11-19 2021-01-29 中国水利水电科学研究院 一种基于模糊综合评价的大坝安全预测系统及方法
CN112507591A (zh) * 2020-12-11 2021-03-16 四川大学 一种基于分离式混合建模技术的重力坝坝顶位移监控方法
CN113486420A (zh) * 2021-06-08 2021-10-08 广东省水利水电科学研究院 水库大坝岸坡稳定性监测方法、装置、设备及介质
CN113742814A (zh) * 2021-07-27 2021-12-03 南方电网调峰调频发电有限公司 大坝安全预警方法、装置、计算机设备和存储介质
CN114547951A (zh) * 2022-04-24 2022-05-27 浙江远算科技有限公司 一种基于数据同化的大坝状态预测方法及系统
CN114992808A (zh) * 2022-06-15 2022-09-02 山东大学 基于组合智能算法的热泵空调热管理控制方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107092759A (zh) * 2017-05-05 2017-08-25 河海大学 基于重力坝坝基参数反演的坝体位移监测点优化布置方法
CN107301282A (zh) * 2017-06-12 2017-10-27 天津大学 基于多源监测时序数据的混凝土坝力学参数反演方法
CN109858064A (zh) * 2018-11-22 2019-06-07 南昌大学 一种基于原型振动响应的大坝及地基弹模动力反演方法
CN110990915A (zh) * 2019-11-20 2020-04-10 河海大学 一种混凝土坝安全监测智能反馈系统及方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107092759A (zh) * 2017-05-05 2017-08-25 河海大学 基于重力坝坝基参数反演的坝体位移监测点优化布置方法
CN107301282A (zh) * 2017-06-12 2017-10-27 天津大学 基于多源监测时序数据的混凝土坝力学参数反演方法
CN109858064A (zh) * 2018-11-22 2019-06-07 南昌大学 一种基于原型振动响应的大坝及地基弹模动力反演方法
CN110990915A (zh) * 2019-11-20 2020-04-10 河海大学 一种混凝土坝安全监测智能反馈系统及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
姚磊华等: "渗流问题单参数摄动随机反演方法", 《计算力学学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112287608A (zh) * 2020-11-19 2021-01-29 中国水利水电科学研究院 一种基于模糊综合评价的大坝安全预测系统及方法
CN112287608B (zh) * 2020-11-19 2023-11-24 中国水利水电科学研究院 一种基于模糊综合评价的大坝安全预测系统及方法
CN112507591A (zh) * 2020-12-11 2021-03-16 四川大学 一种基于分离式混合建模技术的重力坝坝顶位移监控方法
CN112507591B (zh) * 2020-12-11 2022-11-18 四川大学 一种基于分离式混合建模技术的重力坝坝顶位移监控方法
CN113486420A (zh) * 2021-06-08 2021-10-08 广东省水利水电科学研究院 水库大坝岸坡稳定性监测方法、装置、设备及介质
CN113486420B (zh) * 2021-06-08 2022-07-01 广东省水利水电科学研究院 水库大坝岸坡稳定性监测方法、装置、设备及介质
CN113742814A (zh) * 2021-07-27 2021-12-03 南方电网调峰调频发电有限公司 大坝安全预警方法、装置、计算机设备和存储介质
CN113742814B (zh) * 2021-07-27 2023-12-19 南方电网调峰调频发电有限公司 大坝安全预警方法、装置、计算机设备和存储介质
CN114547951A (zh) * 2022-04-24 2022-05-27 浙江远算科技有限公司 一种基于数据同化的大坝状态预测方法及系统
CN114992808A (zh) * 2022-06-15 2022-09-02 山东大学 基于组合智能算法的热泵空调热管理控制方法及系统

Also Published As

Publication number Publication date
CN111695177B (zh) 2022-10-14

Similar Documents

Publication Publication Date Title
CN111695177B (zh) 基于位移场监测数据的大坝力学参数随机反演方法及系统
CN112101660B (zh) 基于阶段性注意力机制的降雨型滑坡位移预测模型及方法
CN109145464B (zh) 融合多目标蚁狮优化与迹稀疏正则化的结构损伤识别方法
CN112528365B (zh) 一种地下基础设施结构健康演化趋势预测方法
CN110866974A (zh) 一种基于三维展示的水工监测系统
CN109035730A (zh) 一种考虑服役环境影响的混凝土坝损伤动力预警方法
CN114547951B (zh) 一种基于数据同化的大坝状态预测方法及系统
CN115759378A (zh) 一种基于数字孪生的大坝安全分析预警系统及方法
CN112904266B (zh) 一种电能表寿命预测方法及装置
CN106677763A (zh) 一种基于动态集成建模的油井动液面预测方法
Li et al. Prediction of the failure point settlement in rockfill dams based on spatial-temporal data and multiple-monitoring-point models
CN112417700B (zh) Eh油站基于状态评价的故障诊断系统
CN112307536B (zh) 一种大坝渗流参数反演方法
CN115659729A (zh) 一种基于结构仿真计算的大坝安全监测分析方法及系统
CN116432123A (zh) 一种基于cart决策树算法的电能表故障预警方法
Wu et al. An approach for deformation modulus mechanism of super-high arch dams
CN102156641A (zh) 一种软件成本置信区间预测方法及系统
CN115577587A (zh) 一种历史建筑健康状态监测方法及系统
CN115238573A (zh) 考虑工况参数的水电机组性能劣化趋势预测方法和系统
CN117150925A (zh) 一种水电工程高陡边坡岩体力学参数反分析方法
CN105488343A (zh) 一种电力二次设备故障概率计算方法
CN116720249A (zh) 一种混凝土坝应力场构建方法
CN116167240A (zh) 一种大坝结构损伤的多测点联合监控方法
Liang et al. Damage identification of arch dams based on post-earthquake residual displacement gradients
Rana et al. Evaluating soil parameters considering probabilistic back analysis for slope failures

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