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

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

Info

Publication number
CN111695177B
CN111695177B CN202010407591.XA CN202010407591A CN111695177B CN 111695177 B CN111695177 B CN 111695177B CN 202010407591 A CN202010407591 A CN 202010407591A CN 111695177 B CN111695177 B CN 111695177B
Authority
CN
China
Prior art keywords
dam
module
inversion
monitoring data
displacement
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.)
Active
Application number
CN202010407591.XA
Other languages
English (en)
Other versions
CN111695177A (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

Abstract

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

Description

基于位移场监测数据的大坝力学参数随机反演方法及系统
技术领域
本发明属于大坝安全领域,具体是一种基于位移场监测数据的大坝力学参数随机反演方法及系统。
背景技术
当前,世界范围内已建成大量的高坝工程,如中国的锦屏一级混凝土拱坝,三峡混凝土重力坝,瑞士的大迪克桑斯重力坝,以及美国的胡佛重力拱坝等。这些高坝工程在防洪、发电、航运、灌溉和供水等方面发挥着至关重要的作用,产生了非常显著的社会和经济效益。大坝在产生巨大工程效益的同时,其工作性态也会随着服役时间的增长而不断变化,会出现一些隐患。部分工程中出现的坝基失稳、坝体开裂、渗透破坏等问题,给相关国家带来了惨重的灾害和巨大的损失。奥地利科恩布莱恩、美国德沃夏克等高坝都发生过严重开裂漏水,修补加固费用巨大;法国马尔帕塞拱坝因坝肩失稳发生溃决,给生命、财产带来巨大损失;中国也有高混凝土坝发生严重裂缝、高压水劈裂等影响安全的实例。因此,对大坝进行合理有效的检查、监控、评价、维护,对确保大坝安全和公共安全具有十分重要的意义。
大坝安全监测数据中含有大量丰富的信息,是大坝结构状态的综合表现,受多种因素的叠加影响。依据实测资料对大坝原型进行反分析,可对大坝前期的设计施工方案进行检验,对当下的工作性态进行评价,对未来的发展趋势进行预测;从而为及时评估和发现大坝的异常迹象提供充分依据,在此基础上制定恰当的水库控制运行计划和大坝的维护管理措施来确保大坝的运行安全,在发生险情时及时发布警报以避免事故的发生。
当前,对大坝结构的反演分析一般是基于各种确定性的数学模型,反演的力学参数某种程度上只能说是较好地反映实际工程情况的“等效参数”。作为挡水建筑物,大坝从建设到运行期间经历复杂内部外部条件变化,本质是一个受多种因素影响的不确定性系统。大坝监测信息的不确定性,主要由本身材料力学参数引起的响应量的摄动和观测误差这两部分组成,而目前的相关分析,没有对观测误差进行识别和剔除。因此,做好监测资料的分析工作,从观测信息中分离出有用的信息,找出影响坝体结构性态的主要因素,才能准确掌握大坝的安全运行状态。
发明内容
本发明的目的是针对现有技术存在的问题,提供一种基于位移场监测数据的大坝力学参数随机反演方法及系统,能对混凝土坝自动化监测信息进行实时处理,对大坝的真实安全状态进行反馈
为实现上述目的,本发明采用的技术方案是:
一种基于位移场监测数据的大坝力学参数随机反演方法,包括以下步骤:
获取大坝各测点的监测数据,对所述监测数据进行预处理;
根据预处理后的监测数据,建立各测点大坝变形的统计模型;
利用所述统计模型,分离大坝位移监测数据中的水压分量;
调用大坝的有限元模型计算对应水荷载下的大坝位移场;
采用一阶摄动随机反演方法或二阶摄动随机反演方法,并调用智能优化算法计算大坝力学参数的平均值和方差;
输出反演结果。
具体地,分离大坝位移监测数据中的水压分量的方法为:在统计模型中将大坝径向位移δ按照成因分为水压分量、温度分量和时效分量,在位移δ中扣除温度分量和时效分量后,得到水压分量如下式所示:
δH=δ-δTθ (1)
其中,δH为大坝径向位移δ中的水压分量;δT为大坝径向位移δ中的温度分量;δθ为大坝径向位移δ中的时效分量。
具体地,采用一阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
Figure GDA0003794705450000021
式中E为数学期望,将式(6)代入式(4)中,可得,
Figure GDA0003794705450000022
定义目标函数为:
Figure GDA0003794705450000023
一阶摄动随机反演的准则为目标函数取得最小值,即在方差最小的意义下,结合优化算法对参数均值
Figure GDA0003794705450000024
进行识别,
Figure GDA0003794705450000025
为X取均值下的位移计算值;
所述随机变量X的方差反演方法为:
对于某一观测时刻k,k=1,2,…,M,令
Figure GDA0003794705450000031
式中,
Figure GDA0003794705450000032
为N×p矩阵,则式(5)可写为:
Figure GDA0003794705450000033
式中,δ(X,Yk)为N×1列阵,{εXk}为p×1列阵,
Figure GDA0003794705450000034
式(10)变换得,
Figure GDA0003794705450000035
将上式两端分别乘以
Figure GDA0003794705450000036
得,
Figure GDA0003794705450000037
变换得,
Figure GDA0003794705450000038
随机变量的方差为:
Figure GDA0003794705450000039
将式(12)的计算结果代入式(10),可得对于某一观测时刻k,由随机变量引起的位移摄动为:
Figure GDA00037947054500000310
式中,{εδXk}为N×1列阵;对于某时刻k,观测误差为:
fk}={εδk}-{εδXk} (16)
联立上述公式(2)至(16),即可计算各测点由随机变量引起的位移方差var(δX)统计值,以及由观测工作产生的位移方差var(δf)统计值。
可选地,采用二阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
将δ(X,Y)在均值处进行包含二次项的Taylor展开,在不考虑随机变量相关条件下,有:
Figure GDA0003794705450000041
对上式取均值,有:
Figure GDA0003794705450000042
对某一观测时刻k,k=1,2,…,M,令:
Figure GDA0003794705450000043
则式(17)可写为:
Figure GDA0003794705450000044
此时,目标函数为:
Figure GDA0003794705450000045
式中
Figure GDA0003794705450000046
Ajk为矩阵Ak的第j行,j=1,2,…,N,即Ajk为1×p列阵;
根据目标函数,即在方差最小的意义下,结合优化算法对二阶意义下的参数均值
Figure GDA0003794705450000047
进行反演;此时
Figure GDA0003794705450000048
是未知的,采用一阶意义下的估计作为估计初值,进行均值
Figure GDA0003794705450000049
计算;
所述随机变量X的方差反演方法为:
式(20)可改写为:
Figure GDA00037947054500000410
式中
Figure GDA00037947054500000411
Figure GDA00037947054500000412
为N×1列阵;
定义目标函数Lk为:
Figure GDA00037947054500000413
式中Bjk为矩阵Bk的第j行,j=1,2,…,N,即Bjk为1×p列阵;Cjk为矩阵Ck的第j行,j=1,2,…,N,即Cjk为1×p列阵;
在目标函数Lk最小的意义下,对{εXk}进行求解,此时{εXk}满足方程组:
Figure GDA0003794705450000051
上式方程组展开后得:
Figure GDA0003794705450000052
式中Ajik是Ajk的第i个数,i=1,2,…,p;
上式可简化为:
Figure GDA0003794705450000053
此时式中Ak,Bk,Ck为N×1列阵;
上式(26)展开得:
Figure GDA0003794705450000054
根据卡尔丹公式求得一元三次方程的三个根,选目标函数最小的实根作为最终解;
此外对于多参数反演问题,调用优化函数直接对4次多项式目标函数进行寻优求解{εXk};
进一步可得随机变量的方差为:
Figure GDA0003794705450000055
根据上式求出
Figure GDA00037947054500000510
后代入式(21),反复迭代,最终达到
Figure GDA0003794705450000056
Figure GDA0003794705450000057
时结束迭代;此时求出
Figure GDA0003794705450000058
Figure GDA0003794705450000059
即为反演的参数和方差。
与上述反演方法相对应的,本发明还提供了一种基于位移场监测数据的大坝力学参数随机反演系统,包括数据库模块,结构信息模块,随机反演模块和信息展示模块;所述数据库模块用于存储大坝的基础信息数据和各类监测数据;所述结构信息模块用于存放大坝工程的结构模型信息;所述随机反演模块用于对大坝力学参数的均值和方差进行计算;所述信息展示模块用于展示大坝的模型、监测数据和反演结果。
具体地,所述数据库模块包括大坝基础信息模块、地理空间信息模块和监测信息模块;
所述大坝基础信息模块用于记载有工程概况、水库大坝设计、施工、运行相关资料;
所述地理空间信息模块用于记载水文和地质情况统计资料和工程地质勘察试验资料;
所述监测信息模块用于记载各类监测设备的原始测值,主要包括:环境量监测数据、变形监测数据、渗流监测数据和压力(应力)监测数据。
具体地,所述结构模型信息包括工程几何模型信息和有限元模型信息。
具体地,所述随机反演模块包括数据前处理模块,智能优化算法模块,一阶摄动随机反演模块和二阶摄动随机反演模块;
所述监测数据前处理模块,用于建立测点变形统计模型,分离得到测点变形的水压分量;
所述智能优化算法模块包括粒子群算法、灰狼算法、鲸鱼优化算法、蚁群算法、人工鱼群算法等智能寻优算法;可通过GUI及按钮触发群智能寻优算法选择框进行选择;
所述一阶摄动随机反演模块用于对一阶摄动意义下力学参数的均值和方差进行识别;
所述二阶摄动随机反演模块用于进一步对二阶摄动意义下力学参数的均值和方差进行识别。
具体地,所述信息展示模块包括模型轻量化展示、监控信息浏览和评估结果输出功能;其中,借助于BIMFACE轻量化模型转换功能,将所建立的大坝及监测测点三维信息模型在Web端进行展示,同时,可通过浏览器查看大坝监测测点信息和实时评估的结果信息。
与现有技术相比,本发明的有益效果是:
(1)本发明提出的随机反演方法考虑了大坝结构力学性质的不确定性以及观测误差的不确定性,联合空间多测点位移,能够对大坝结构力学参数进行准确识别;
(2)本发明通过随机反演方法能够将大坝的真实位移量和观测误差进行分离,提高数据的精确度和评估结果的可信度;
(3)本发明提出的随机反演系统能对大坝的自动化监测数据进行全链条处理,对大坝的真实安全状态进行实时反馈,从而协助技术人员高效地完成大坝安全监测及管理工作,具有重要的社会经济效益。
附图说明
图1为本发明一种基于位移场监测数据的大坝力学参数随机反演方法的流程图;
图2为本发明一种基于位移场监测数据的大坝力学参数随机反演系统的结构框图;
图3为本发明实施例中某碾压混凝土坝11#坝段有限元模型图;
图4为本发明实施例中大坝各测点变形的水压分量分离结果示意图;
图5为本发明实施例中垂线PL11-1测点位移反馈结果示意图。
具体实施方式
下面将结合本发明中的附图,对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动条件下所获得的所有其它实施例,都属于本发明保护的范围。
如图1所示,本实施例提供了一种基于位移场监测数据的大坝力学参数随机反演方法,包括以下步骤:
获取大坝各测点的监测数据,对所述监测数据进行预处理;
根据预处理后的监测数据,建立各测点大坝变形的统计模型;
利用所述统计模型,分离大坝位移监测数据中的水压分量;
调用大坝的有限元模型计算对应水荷载下的大坝位移场;
采用一阶摄动随机反演方法或二阶摄动随机反演方法,并调用智能优化算法计算大坝力学参数的平均值和方差;
输出反演结果。
具体地,分离大坝位移监测数据中的水压分量的方法为:在统计模型中将大坝径向位移δ按照成因分为水压分量、温度分量和时效分量,在位移δ中扣除温度分量和时效分量后,得到水压分量如下式所示:
δH=δ-δTθ (1)
其中,δH为大坝径向位移δ中的水压分量;δT为大坝径向位移δ中的温度分量;δθ为大坝径向位移δ中的时效分量。
具体地,采用一阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
Figure GDA0003794705450000071
式中E为数学期望,将式(6)代入式(4)中,可得,
Figure GDA0003794705450000081
定义目标函数为:
Figure GDA0003794705450000082
一阶摄动随机反演的准则为目标函数取得最小值,即在方差最小的意义下,结合优化算法对参数均值
Figure GDA00037947054500000814
进行识别,
Figure GDA0003794705450000083
为X取均值下的位移计算值;
所述随机变量X的方差反演方法为:
对于某一观测时刻k,k=1,2,…,M,令
Figure GDA0003794705450000084
式中,
Figure GDA0003794705450000085
为N×p矩阵,则式(5)可写为:
Figure GDA0003794705450000086
式中,δ(X,Yk)为N×1列阵,{εXk}为p×1列阵,
Figure GDA0003794705450000087
式(10)变换得,
Figure GDA0003794705450000088
将上式两端分别乘以
Figure GDA0003794705450000089
得,
Figure GDA00037947054500000810
变换得,
Figure GDA00037947054500000811
随机变量的方差为:
Figure GDA00037947054500000812
将式(12)的计算结果代入式(10),可得对于某一观测时刻k,由随机变量引起的位移摄动为:
Figure GDA00037947054500000813
式中,{εδXk}为N×1列阵;对于某时刻k,观测误差为:
fk}={εδk}-{εδXk} (16)
联立上述公式(2)至(16),即可计算各测点由随机变量引起的位移方差var(δX)统计值,以及由观测工作产生的位移方差var(δf)统计值。
可选地,采用二阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
将δ(X,Y)在均值处进行包含二次项的Taylor展开,在不考虑随机变量相关条件下,有:
Figure GDA0003794705450000091
对上式取均值,有:
Figure GDA0003794705450000092
对某一观测时刻k,k=1,2,…,M,令:
Figure GDA0003794705450000093
则式(17)可写为:
Figure GDA0003794705450000094
此时,目标函数为:
Figure GDA0003794705450000095
式中
Figure GDA0003794705450000096
Ajk为矩阵Ak的第j行,j=1,2,…,N,即Ajk为1×p列阵;
根据目标函数,即在方差最小的意义下,结合优化算法对二阶意义下的参数均值
Figure GDA0003794705450000097
进行反演;此时
Figure GDA0003794705450000098
是未知的,采用一阶意义下的估计作为估计初值,进行均值
Figure GDA0003794705450000099
计算;
所述随机变量X的方差反演方法为:
式(20)可改写为:
Figure GDA00037947054500000910
式中
Figure GDA0003794705450000101
Figure GDA0003794705450000102
为N×1列阵;
定义目标函数Lk为:
Figure GDA0003794705450000103
式中Bjk为矩阵Bk的第j行,j=1,2,…,N,即Bjk为1×p列阵;Cjk为矩阵Ck的第j行,j=1,2,…,N,即Cjk为1×p列阵;
在目标函数Lk最小的意义下,对{εXk}进行求解,此时{εXk}满足方程组:
Figure GDA0003794705450000104
上式方程组展开后得:
Figure GDA0003794705450000105
式中Ajik是Ajk的第i个数,i=1,2,…,p;
上式可简化为:
Figure GDA0003794705450000106
此时式中Ak,Bk,Ck为N×1列阵;
上式(26)展开得:
Figure GDA0003794705450000107
根据卡尔丹公式求得一元三次方程的三个根,选目标函数最小的实根作为最终解;
此外对于多参数反演问题,调用优化函数直接对4次多项式目标函数进行寻优求解{εXk};
进一步可得随机变量的方差为:
Figure GDA0003794705450000111
根据上式求出
Figure GDA0003794705450000112
后代入式(21),反复迭代,最终达到
Figure GDA0003794705450000113
Figure GDA0003794705450000114
时结束迭代;此时求出
Figure GDA0003794705450000115
Figure GDA0003794705450000116
即为反演的参数和方差。
与上述反演方法相对应的,如图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>x</sub>)(mm) 0.13 0.28 0.47 0.61 0.81
std(δ′<sub>f</sub>)(mm) 0.08 0.09 0.07 0.09 0.07
同时,如图5所示,本实施例为位移反馈结果提供一个概率输出,即大坝测点位移的置信区间。如95%的置信区间可定义为
Figure GDA0003794705450000131
Figure GDA0003794705450000132
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。

Claims (7)

1.一种基于位移场监测数据的大坝力学参数随机反演方法,其特征在于,包括以下步骤:
获取大坝各测点的监测数据,对所述监测数据进行预处理;
根据预处理后的监测数据,建立各测点大坝变形的统计模型;
利用所述统计模型,分离大坝位移监测数据中的水压分量;
调用大坝的有限元模型计算对应水荷载下的大坝位移场;
采用一阶摄动随机反演方法或二阶摄动随机反演方法,并调用智能优化算法计算大坝力学参数的平均值和方差;
输出反演结果;
采用一阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
假设X为具有固定均值的随机变量,则表示为:
Figure FDA0003763466240000011
式中,X={x1,x2,…xp}T,p为随机变量的个数;
Figure FDA0003763466240000012
为X的均值,
Figure FDA0003763466240000013
εX为随机变量摄动值;
位移量δ(X,Y)是关于X,Y的函数,Y为确定性参数,Y={y1,y2,…yq}T,q为确定性参数的个数,若观测位移δ为非平稳随机过程变量,具有确定的趋势项,则δ表示为:
Figure FDA0003763466240000014
式中,
Figure FDA0003763466240000015
为δ的均值;εδ为位移摄动量,由随机变量摄动引起的位移摄动εδX和观测误差εf两部分组成;假定观测误差εf相互独立,数学期望为0且服从正态分布;
设有N个位移观测点,每个观测点有M个测值,则第j个测点的位移方差统计表达式为:
Figure FDA0003763466240000016
式中,var为方差算子;δ(X,Yjk)为第j个测点第k次观测量,j=1,2,···,N,k=1,2,···,M;令δjk=δ(X,Yjk),
Figure FDA0003763466240000017
Figure FDA0003763466240000018
为δjk的均值;
对于一阶摄动随机反演,将δ(X,Y)在均值处进行Taylor展开,并将高于一次的项舍去,则有:
Figure FDA0003763466240000019
对上式取均值,有:
Figure FDA0003763466240000021
式中E为数学期望,将式(6)代入式(4)中,可得,
Figure FDA0003763466240000022
定义目标函数为:
Figure FDA0003763466240000023
一阶摄动随机反演的准则为目标函数取得最小值,即在方差最小的意义下,结合优化算法对参数均值
Figure FDA0003763466240000024
进行识别,
Figure FDA0003763466240000025
为X取均值下的位移计算值;
所述随机变量X的方差反演方法为:
对于某一观测时刻k,k=1,2,···,M,令
Figure FDA0003763466240000026
式中,
Figure FDA0003763466240000027
为N×p矩阵,则式(5)可写为:
Figure FDA0003763466240000028
式中,δ(X,Yk)为N×1列阵,{εXk}为p×1列阵,
Figure FDA0003763466240000029
式(10)变换得,
Figure FDA00037634662400000210
将上式两端分别乘以
Figure FDA00037634662400000211
得,
Figure FDA00037634662400000212
变换得,
Figure FDA00037634662400000213
随机变量的方差为:
Figure FDA00037634662400000214
将式(12)的计算结果代入式(10),可得对于某一观测时刻k,由随机变量引起的位移摄动为:
Figure FDA0003763466240000031
式中,{εδXk}为N×1列阵;对于某时刻k,观测误差为:
fk}={εδk}-{εδXk} (16)
联立上述公式(2)至(16),即可计算各测点由随机变量引起的位移方差var(δX)统计值,以及由观测工作产生的位移方差var(δf)统计值。
2.根据权利要求1所述的一种基于位移场监测数据的大坝力学参数随机反演方法,其特征在于,分离大坝位移监测数据中的水压分量的方法为:在统计模型中将大坝径向位移δ按照成因分为水压分量、温度分量和时效分量,在位移δ中扣除温度分量和时效分量后,得到水压分量如下式所示:
δH=δ-δTθ (1)
其中,δH为大坝径向位移δ中的水压分量;δT为大坝径向位移δ中的温度分量;δθ为大坝径向位移δ中的时效分量。
3.根据权利要求1所述的一种基于位移场监测数据的大坝力学参数随机反演方法,其特征在于,采用二阶摄动随机反演方法计算大坝力学参数的平均值和方差的方法为:
将δ(X,Y)在均值处进行包含二次项的Taylor展开,在不考虑随机变量相关条件下,有:
Figure FDA0003763466240000032
对上式取均值,有:
Figure FDA0003763466240000033
对某一观测时刻k,k=1,2,···,M,令:
Figure FDA0003763466240000034
则式(17)可写为:
Figure FDA0003763466240000041
此时,目标函数为:
Figure FDA0003763466240000042
式中
Figure FDA0003763466240000043
Ajk为矩阵Ak的第j行,j=1,2,···,N,即Ajk为1×p列阵;
根据目标函数,即在方差最小的意义下,结合优化算法对二阶意义下的参数均值
Figure FDA0003763466240000044
进行反演;此时
Figure FDA0003763466240000045
是未知的,采用一阶意义下的估计作为估计初值,进行均值
Figure FDA0003763466240000046
计算;
所述随机变量X的方差反演方法为:
式(20)可改写为:
Figure FDA0003763466240000047
式中
Figure FDA0003763466240000048
为N×1列阵;
定义目标函数Lk为:
Figure FDA0003763466240000049
式中Bjk为矩阵Bk的第j行,j=1,2,···,N,即Bjk为1×p列阵;Cjk为矩阵Ck的第j行,j=1,2,···,N,即Cjk为1×p列阵;
在目标函数Lk最小的意义下,对{εXk}进行求解,此时{εXk}满足方程组:
Figure FDA00037634662400000410
上式方程组展开后得:
Figure FDA00037634662400000411
式中Ajik是Ajk的第i个数,i=1,2,···,p;
上式可简化为:
Figure FDA0003763466240000051
此时式中Ak,Bk,Ck为N×1列阵;
上式(26)展开得:
Figure FDA0003763466240000052
根据卡尔丹公式求得一元三次方程的三个根,选目标函数最小的实根作为最终解;
对于多参数反演问题,调用优化函数直接对4次多项式目标函数进行寻优求解{εXk};
可得随机变量的方差为:
Figure FDA0003763466240000053
根据上式求出
Figure FDA0003763466240000054
后代入式(21),反复迭代,最终达到
Figure FDA0003763466240000055
Figure FDA0003763466240000056
时结束迭代;此时求出
Figure FDA0003763466240000057
Figure FDA0003763466240000058
即为反演的参数和方差。
4.一种基于位移场监测数据的大坝力学参数随机反演系统,基于权利要求1至3任一项所述的大坝力学参数随机反演方法,其特征在于,包括数据库模块,结构信息模块,随机反演模块和信息展示模块;所述数据库模块用于存储大坝的基础信息数据和各类监测数据;所述结构信息模块用于存放大坝工程的结构模型信息;所述随机反演模块用于对大坝力学参数的均值和方差进行计算;所述信息展示模块用于展示大坝的模型、监测数据和反演结果。
5.根据权利要求4所述的一种基于位移场监测数据的大坝力学参数随机反演系统,其特征在于,所述数据库模块包括大坝基础信息模块、地理空间信息模块和监测信息模块;
所述大坝基础信息模块用于记载有工程概况、水库大坝设计、施工、运行相关资料;
所述地理空间信息模块用于记载水文和地质情况统计资料和工程地质勘察试验资料;
所述监测信息模块用于记载各类监测设备的原始测值,主要包括:环境量监测数据、变形监测数据、渗流监测数据和压力监测数据。
6.根据权利要求4所述的一种基于位移场监测数据的大坝力学参数随机反演系统,其特征在于,所述结构模型信息包括工程几何模型信息和有限元模型信息。
7.根据权利要求4所述的一种基于位移场监测数据的大坝力学参数随机反演系统,其特征在于,所述随机反演模块包括数据前处理模块,智能优化算法模块,一阶摄动随机反演模块和二阶摄动随机反演模块;
所述监测数据前处理模块,用于建立测点变形统计模型,分离得到测点变形的水压分量;
所述智能优化算法模块包括粒子群算法、灰狼算法、鲸鱼优化算法、蚁群算法、人工鱼群算法中的一个或多个智能寻优算法;
所述一阶摄动随机反演模块用于对一阶摄动意义下力学参数的均值和方差进行识别;
所述二阶摄动随机反演模块用于对二阶摄动意义下力学参数的均值和方差进行识别。
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 CN111695177A (zh) 2020-09-22
CN111695177B true 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)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112287608B (zh) * 2020-11-19 2023-11-24 中国水利水电科学研究院 一种基于模糊综合评价的大坝安全预测系统及方法
CN112507591B (zh) * 2020-12-11 2022-11-18 四川大学 一种基于分离式混合建模技术的重力坝坝顶位移监控方法
CN113486420B (zh) * 2021-06-08 2022-07-01 广东省水利水电科学研究院 水库大坝岸坡稳定性监测方法、装置、设备及介质
CN113742814B (zh) * 2021-07-27 2023-12-19 南方电网调峰调频发电有限公司 大坝安全预警方法、装置、计算机设备和存储介质
CN114547951B (zh) * 2022-04-24 2022-07-22 浙江远算科技有限公司 一种基于数据同化的大坝状态预测方法及系统
CN114992808B (zh) * 2022-06-15 2023-06-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
渗流问题单参数摄动随机反演方法;姚磊华等;《计算力学学报》;20081231;第25卷(第6期);全文 *

Also Published As

Publication number Publication date
CN111695177A (zh) 2020-09-22

Similar Documents

Publication Publication Date Title
CN111695177B (zh) 基于位移场监测数据的大坝力学参数随机反演方法及系统
CN108984873B (zh) 供水管网实时漏损检测方法、设备、系统及存储介质
CN109145464B (zh) 融合多目标蚁狮优化与迹稀疏正则化的结构损伤识别方法
CN110866974A (zh) 一种基于三维展示的水工监测系统
CN115759378A (zh) 一种基于数字孪生的大坝安全分析预警系统及方法
CN112528365B (zh) 一种地下基础设施结构健康演化趋势预测方法
CN114547951B (zh) 一种基于数据同化的大坝状态预测方法及系统
CN109035730A (zh) 一种考虑服役环境影响的混凝土坝损伤动力预警方法
CN112904266B (zh) 一种电能表寿命预测方法及装置
CN109359770B (zh) 一种基于机器学习预测中暑发生的模型及方法
CN112503400B (zh) 一种供水管网测压点多目标优化布置方法
CN108229726A (zh) 一种碎块石边坡灾变险情预警预报方法
CN104134013B (zh) 一种风力机叶片模态分析方法
CN115659729A (zh) 一种基于结构仿真计算的大坝安全监测分析方法及系统
CN111241660B (zh) 一种基于云计算技术的土木工程实验检测系统及方法
CN116432123A (zh) 一种基于cart决策树算法的电能表故障预警方法
Wu et al. An approach for deformation modulus mechanism of super-high arch dams
CN116167240A (zh) 一种大坝结构损伤的多测点联合监控方法
CN115713312A (zh) 一种水利水电工程环境量不利工况组合识别方法及装置
CN113418146A (zh) 一种用于供水管网的漏损辅助定位控制方法
CN111767649A (zh) 一种软基地区变电站地质变形安全评估方法
CN206627073U (zh) 基于3d激光扫描的地下工程施工快速监测预测装置
CN117709488B (zh) 基于RUN-XGBoost的大坝渗流预测方法
CN113569444B (zh) 一种堤身材料渗透系数随机序贯反演方法
CN103955953A (zh) 一种面向数字土壤制图的地形协同变量选取方法

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