CN112507591A - 一种基于分离式混合建模技术的重力坝坝顶位移监控方法 - Google Patents
一种基于分离式混合建模技术的重力坝坝顶位移监控方法 Download PDFInfo
- Publication number
- CN112507591A CN112507591A CN202011463759.5A CN202011463759A CN112507591A CN 112507591 A CN112507591 A CN 112507591A CN 202011463759 A CN202011463759 A CN 202011463759A CN 112507591 A CN112507591 A CN 112507591A
- Authority
- CN
- China
- Prior art keywords
- dam
- dam crest
- displacement
- model
- crest 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.)
- Granted
Links
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 164
- 238000012544 monitoring process Methods 0.000 title claims abstract description 103
- 230000005484 gravity Effects 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005516 engineering process Methods 0.000 title claims abstract description 16
- 238000000926 separation method Methods 0.000 title claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 69
- 230000032683 aging Effects 0.000 claims abstract description 26
- 230000000694 effects Effects 0.000 claims abstract description 15
- 230000002277 temperature effect Effects 0.000 claims abstract description 11
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 9
- 230000002159 abnormal effect Effects 0.000 claims abstract description 4
- 230000036760 body temperature Effects 0.000 claims description 16
- 238000013461 design Methods 0.000 claims description 12
- 239000000463 material Substances 0.000 claims description 12
- 230000000737 periodic effect Effects 0.000 claims description 12
- 239000004567 concrete Substances 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 238000011144 upstream manufacturing Methods 0.000 claims description 9
- 239000011435 rock Substances 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000005192 partition Methods 0.000 claims description 5
- 238000009529 body temperature measurement Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 238000000513 principal component analysis Methods 0.000 claims description 4
- 208000031636 Body Temperature Changes Diseases 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000001550 time effect Effects 0.000 description 5
- 238000013179 statistical model Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000011382 roller-compacted concrete Substances 0.000 description 1
Images
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/08—Construction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- General Engineering & Computer Science (AREA)
- Public Health (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了一种基于分离式混合建模技术的重力坝坝顶位移监控方法,属于大坝安全监控方法技术领域,采用经验模态分解提取重力坝坝顶位移监测数据序列的时效趋势项,再根据坝体等温条件分离监测数据序列中的水压效应量,分别构建坝顶位移的时效分量和水压分量模型;利用建立的时效和水压分量模型从监测数据序列中分解出温度效应量,建立坝顶位移的温度分量模型;最后,叠加时效、水压和温度分量得到重力坝坝顶位移监控模型。根据监控模型得到的坝顶位移估计值和模型误差的概率密度函数,构造坝顶位移的置信区间(置信度可以取95%),对坝顶位移是否异常进行监控。当新获取的坝顶位移监测数据位于置信区间之外时,则发出数据异常预警。
Description
技术领域
本发明属于大坝安全监控方法技术领域,具体涉及一种基于分离式混合建模技术的重力坝坝顶位移监控方法。
背景技术
目前重力坝坝顶位移监控常用的数学模型是多元统计回归模型,如HST模型。多元统计回归模型的数学形式简单,易于使用,对历史监测数据的拟合精度较高。但在多元统计回归模型构建过程中,水压、温度和时效分量模型参数是按最小二乘法同时确定的,为了使整体模型对监测数据的拟合精度最高,各分量模型将相互补偿,从而导致整体模型过度拟合,以及各分量模型不能真实反映坝顶位移的水压、温度和时间效应,模型的预测精度较低,从而增大坝顶位移监控的误判率。
发明内容
本发明的目的在于提供一种基于分离式混合建模技术的重力坝坝顶位移监控方法,解决现有技术中为了使整体模型对监测数据的拟合精度最高,各分量模型将相互补偿,从而导致整体模型过度拟合,以及各分量模型不能真实反映坝顶位移的水压、温度和时间效应,模型的预测精度较低,从而增大坝顶位移监控的误判率的技术问题。
本发明提供了一种基于分离式混合建模技术的重力坝坝顶位移监控方法,包括以下步骤:
S1.采用经验模态分解(EMD)技术将重力坝坝顶位移监测数据序列分解为一系列不同频率的固有模态函数分量(IMF),其中频率最低且具有单调性的IMF为重力坝坝顶位移监测数据序列中的时效趋势项,将其余IMF进行叠加重构即可得到重力坝坝顶位移监测数据序列的周期变化项,该周期变化项是由库水位和坝体温度变化引起的;
S2.根据EMD提取的时效趋势项,采用复合指数函数构建重力坝坝顶位移的时效分量模型,即
F3[θ(t)]=h*exp[kθ(t)j]+q (1)
式中,θ(t)=(t-t0)100,t为时间,t0为监测数据序列的起始时间;h、k、j、q为模型参数;
S3.采用EMD分解、重构得到的重力坝坝顶位移监测数据序列的周期变化项是库水位和坝体温度变化共同作用的结果,由于正常运行条件下重力坝一般都在线弹性范围内工作,因此,当两个时刻的坝体温度相等或非常接近时,可以认为这两个时刻的坝顶位移差值仅是库水位变化引起的。所以,根据坝体等温条件(即两个时刻坝体温度测点的测值相等或温度测值之差小于0.2℃),可以从坝顶位移监测数据序列的周期变化项中分离出库水位变化导致的位移变化值,即坝顶位移监测数据序列的水压效应量;
S4.根据坝顶位移监测数据序列的水压效应量,构建坝顶位移的水压分量模型;
S5.利用建立的时效和水压分量模型从重力坝坝顶位移监测数据序列中分解出温度效应量(即坝体温度变化引起的坝顶位移),构建坝顶位移的温度分量模型;
S6.叠加步骤S2、S4和S5构建的时效、水压和温度分量模型,最终得到重力坝坝顶位移的监控模型;
S7.计算坝顶位移实测值和监控模型估计值的差值(即模型误差),采用核函数拟合模型误差的概率密度函数,根据坝顶位移的模型估计值和模型误差的概率密度函数,构造坝顶位移的置信区间(置信度可以取95%),对坝顶位移是否异常进行监控。当新获取的坝顶位移监测数据位于置信区间之外时,则发出数据异常预警。
进一步的,步骤S2中h、k、j、q模型参数,根据重力坝坝顶位移监测数据序列的时效趋势项,采用1stOpt15PRO软件确定。
进一步的,所述步骤S4中构建坝顶位移的水压分量模型包括以下步骤:
1)根据坝体结构轮廓尺寸、材料分区和坝址地形地质条件等资料,建立重力坝坝段的三维有限元模型。根据坝顶位移测点位置在有限元模型中布置相应结点,以便于提取位移计算结果。
2)确定需要反演分析的坝体混凝土和地基岩体弹性模量,以及各弹性模量的取值区间。根据坝与地基材料弹性模量的取值区间和库水位变化范围,确定弹性模量和库水位的正交试验设计水平数和相应数值。采用正交试验设计对弹性模量和库水位的各水平值进行组合。
3)采用第1)步建立的重力坝坝段有限元模型,按照弹性模量和库水位的不同组合,计算相应的坝顶位移值。根据弹性模量和库水位的不同组合,以及相应的坝顶位移计算值,采用1stOpt15PRO软件确定坝顶位移水压分量模型中的系数a、bi、ci、di、ei,见式(2)。
式中,H(t)为某个时刻库水位对应的上游水深;xi(i=1,2,…,n)代表坝体混凝土和地基岩体的弹性模量。
4)根据第(3)步得到的坝顶位移监测数据序列的水压效应量及对应的上游水深变化值,采用1stOpt15PRO软件对式(2)中的xi(i=1,2,...,n)(即坝体混凝土和地基岩体的弹性模量)进行反分析。
5)将第4)步得到的坝体混凝土和地基岩体弹性模量的反分析值代入式(2),即可建立以上游水深为变量的坝顶位移水压分量模型。
进一步的,所述步骤S5中构建坝顶位移的温度分量模型包括以下步骤:
1)计算坝体温度测点的测值与坝顶位移测值的Spearman相关系数,根据Spearman相关系数绝对值大于0.65的条件,筛选出与坝顶位移测值相关性较强的坝体温度测点,用于构建坝顶位移的温度分量模型。
2)对筛选出的坝体温度测点的测值序列进行主成分分析,并根据主成分对原测值方差的累计贡献率大于95%的条件,确定用于构建温度分量模型的主成分个数。
3)采用坝体温度测值的主成分构建坝顶位移温度分量模型,即
式中,Pi(t)(i=1,2,…,m)为坝体温度测值的主成分,m为对原测值方差的累计贡献率大于95%的主成分个数;fi(i=1,2,…,m)和g为待定系数,可以根据从坝顶位移监测数据序列中分解出的温度效应量,采用最小二乘法确定。
进一步的,采用MATLAB软件中的核函数拟合模型误差的概率密度函数。
本发明的有益效果为:
1.目前重力坝坝顶位移监控常用的数学模型是多元统计回归模型,如HST模型。多元统计回归模型存在过度拟合,泛化能力较低,不能真实反映坝顶位移的水压、温度和时间效应,以及坝顶位移监控误判率高等问题。针对多元统计回归模型的问题,本发明提出了一种基于分离式混合建模技术的重力坝坝顶位移监控方法,通过从历史监测数据序列中分解出坝顶位移的水压、温度和时间效应量,再分别构建坝顶位移的水压、温度和时效分量模型,最后,叠加时效、水压和温度分量得到重力坝坝顶位移的监控模型。本发明提出的方法避免了多元统计回归模型的过度拟合问题,能更准确模拟重力坝坝顶位移的水压、温度与时间效应,模型泛化能力较高,降低了重力坝坝顶位移监控的误判率。
附图说明
为了更清楚地说明本发明实施方式的技术方案,下面将对实施方式中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施方式,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明重力坝坝顶位移监控方法流程图;
图2为某重力坝9#坝段监测仪器布置图;
图3为某重力坝9#坝段坝顶LA9测点顺河向位移的EMD分解和重构结果;
图4为某重力坝9#坝段坝顶LA9测点顺河向位移的时效分量模型及拟合效果;
图5为某重力坝9#坝段坝顶LA9测点顺河向位移差(满足坝体等温条件)及对应的库水深度;
图6为某重力坝9#坝段坝体三维有限元模型及材料分区;
图7为某重力坝9#坝段坝与地基三维有限元模型及地基材料分区;
图8为某重力坝9#坝段坝顶LA9测点处顺河向位移的计算值和模型拟合值;
图9为某重力坝9#坝段坝顶LA9测点处顺河向位移计算值和模型拟合值的相关图;
图10为不同库水位下坝顶顺河向位移水压分量模型的位移输出值;
图11为坝体温度计测值主成分PC1和PC2的历时过程线;
图12为某重力坝9#坝段坝顶LA9测点顺河向位移温度分量的实测值与模型拟合值;
图13为某重力坝9#坝段坝顶顺河向位移监控模型误差的概率密度函数;
图14为采用本发明方法构建的模型进行某重力坝9#坝段坝顶位移监控的结果;
图15为采用HST统计模型进行某重力坝9#坝段坝顶位移监控的结果。
具体实施方式
为使本发明实施方式的目的、技术方案和优点更加清楚,下面将结合本发明实施方式中的附图,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式是本发明一部分实施方式,而不是全部的实施方式。通常在此处附图中描述和展示的本发明实施方式的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
需要说明的是,在不冲突的情况下,本发明中的实施方式及实施方式中的特征可以相互组合。
在本发明的描述中,需要说明的是,术语“上”、“下”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,或者是该发明产品使用时惯常摆放的方位或位置关系,或者是本领域技术人员惯常理解的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
在本发明的描述中,还需要说明的是,除非另有明确的规定和限定,术语“设置”、“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1
本发明的实施方式提供了一种基于分离式混合建模技术的重力坝坝顶位移监控方法,包括以下步骤:
某水电站枢纽拦河大坝采用碾压混凝土重力坝,自左岸至右岸依次布置左岸挡水坝段、河床溢流坝段和右岸挡水坝段。重力坝坝顶高程1334.00m,坝顶长516.00m,最低建基高程1166.00m,最大坝高168m,最大坝底宽度153.2m。本实例主要针对该重力坝9#挡水坝段的坝顶位移,根据其历史监测数据序列,采用本发明方法建立坝顶位移监控模型,构造坝顶位移置信区间,对坝顶位移是否异常进行监控,并与HST模型进行对比。9#坝段长度为20m,建基面高程为1186m,坝段高148m。为了监控运行期间的大坝顺河向位移,在坝顶布置了一套真空激光准直位移监测系统,其中9#坝段坝顶位移测点的编号为LA9。此外,在9#坝段坝体内埋设了18支温度计以监测坝体温度变化,反映坝体温度场分布情况,目前有14支温度计正常在测。9#坝段监测仪器埋设位置信息见表1,监测仪器布置见图2。
采用本发明方法建立9#坝段坝顶位移的监控模型,具体步骤如下:
S1.采用经验模态分解(EMD)技术将9#坝段坝顶位移测点LA9的顺河向位移监测数据序列分解为一系列不同频率的固有模态函数分量(IMF),其中频率最低且具有单调性的IMF为位移监测数据序列中的时效趋势项。将其余IMF进行叠加重构即可得到坝顶位移监测数据序列的周期变化项,该周期变化项是由库水位和坝体温度变化引起的。LA9测点顺河向位移的EMD分解及重构结果见图3。
S2.根据EMD提取的时效趋势项,采用复合指数函数构建9#坝段坝顶顺河向位移的时效分量模型,模型中的系数采用1stOpt15PRO软件确定。建立的时效分量模型如式(4),模型拟合效果见图4。
F3[θ(t)]=-1.21*exp[-0.08θ(t)2.14]+4.45(4)
S3.采用EMD分解、重构得到的坝顶位移监测数据序列的周期变化项是库水位和坝体温度变化共同作用的结果。根据坝体等温条件(即两个时刻14支坝体温度计的测值相等或温度测值之差小于0.2℃),从坝顶位移监测数据序列的周期变化项中分离出库水位变化引起的位移变化值,即坝顶位移监测数据序列的水压效应量。在本实例中,共筛选出63组数据,每组数据包含两个时刻的坝顶位移(满足坝体温度相等的条件)及对应的库水深度,如图5所示。
S4.根据分离出的位移监测数据序列的水压效应量,构建坝顶顺河向位移的水压分量模型。具体步骤如下:
1)根据9#坝段结构轮廓尺寸、材料分区和坝基地质条件等资料,采用ANSYS软件建立该坝段三维有限元模型,材料本构模型选用线弹性模型。根据坝顶LA9测点位置在有限元模型中布置相应结点,以方便提取位移计算结果。有限元模型的模拟范围为:自建基面向下延伸2.5倍坝高(约370m),上下游方向也延伸2.5倍坝高(约370m)。模型坐标系统的X轴正方向为上游指向下游,Y轴的正方向为竖直向上,Z轴的正方向为左岸指向右岸。采用8节点六面体SOLID45单元进行模型离散,整个模型共划分为58452个单元和65931个结点。坝与地基材料的设计参数见表2,有限元模型如图6和图7所示。
2)选取坝体混凝土R1、R2、R3、R4、R5、C2,以及II类、III(1)类、III(2)类基岩的弹性模量作为反演参数。根据坝与地基材料弹性模量的取值区间和库水位变化范围,确定弹性模量和库水位的正交试验设计水平数和相应数值。每种材料的弹性模量取7个水平,并在弹性模量设计值上下一定范围选取设计水平值。在确定库水深度的设计水平值时,以库水位监测数据的历史最大、最小值作为上下限,选择129.33m、132.00m、134.67m、136.00m、138.67m、141.00m、144.00m共7个水位。由于下游水深较小,对坝体位移的影响可以忽略不计。采用正交试验设计对弹性模量和库水位的各水平值进行组合,共得到81组试验样本,结果见表3。
3)采用第1)步建立的重力坝坝段有限元模型,按照弹性模量和库水位的不同组合,计算相应的坝顶顺河向位移值(LA9测点处)。根据81组计算结果,采用1stOpt15PRO软件确定坝顶顺河向位移水压分量模型中的系数a、bi、ci、di、ei,见式(5)。模型拟合效果见图8和图9,模型复相关系数为0.97,说明模型拟合效果很好。
4)根据第(3)步得到的坝顶顺河向位移监测数据序列的水压效应量及对应的上游水深变化值,采用1stOpt15PRO软件对式(5)中的xi(i=1,2,...,n)(即坝体混凝土和基岩弹性模量)进行反分析,结果见表4。
5)将第4)步得到的坝体混凝土和基岩弹性模量的反分析值代入式(5),即可建立以库水深为变量的坝顶顺河向位移水压分量模型,如式(6)所示。不同库水位下坝顶顺河向位移水压分量模型的位移输出值见图10。
F1[H(t)]=0.07+48353808.35*[-1.456180*10-11H(t)+1.202051*10-13H(t)2-1.797984*10-16H(t)3]
(6)
S5.利用已经建立的时效和水压分量模型从9#坝段坝顶LA9测点的顺河向位移监测数据序列中分解出温度效应量(即坝体温度变化引起的坝顶位移),构建坝顶顺河向位移的温度分量模型。具体步骤如下:
1)计算坝体温度计测值与坝顶顺河向位移测值的Spearman相关系数,根据Spearman相关系数绝对值大于0.65的条件,筛选出与坝顶位移测值相关性较强的坝体温度计T4、T8、T10、T14,用于构建坝顶顺河向位移的温度分量模型。
2)对筛选出的4支坝体温度计的测值序列进行主成分分析,结果见表5。由于前两个主成分PC1和PC2对原测值方差的累计贡献率达到了99.42%,故采用PC1和PC2构建坝顶顺河向位移的温度分量模型,PC1和PC2的历时过程线见图11。
根据从坝顶顺河向位移监测数据序列中分解出的温度效应量,采用最小二乘法确定温度分量模型中的系数。坝顶顺河向位移温度分量模型如式(7)所示,模型拟合效果见图12。
F2[P(t)]=-0.22386P1(t)+0.16179P2(t)-62.28 (7)
S6.叠加第(2)步、第(4)步和第(5)步构建的时效、水压和温度分量模型,即式(4)、式(6)和式(7),最终得到9#坝段坝顶顺河向位移的监控模型如下:
F(t)=48353808.35*[-1.456180*10-11H(t)+1.202051*10-13H(t)2-1.797984*10-16H(t)3]-0.22386P1(t)+0.16179P2(t)-1.21*exp[-0.08θ(t)2.14]-57.76 (8)
S7.计算坝顶顺河向位移实测值和监控模型估计值的差值(即模型误差)。采用MATLAB软件中的核函数拟合模型误差的概率密度函数,如图13所示。根据坝顶顺河向位移的模型估计值和模型误差的概率密度函数,构造坝顶顺河向位移的置信区间(置信度取95%),当新获取的坝顶顺河向位移监测数据位于置信区间之外时,则发出数据异常预警。
本实例利用2015年5月1日至2017年5月1日期间的664个监测数据构造坝顶顺河向位移监控模型,再将2017年5月2日至2017年8月1日期间的87个监测数据假定为“新获取的坝顶顺河向位移监测数据”,进行坝顶位移监控模拟,以检验监控模型的性能。图14和图15分别为本发明监控模型和HST统计模型的位移监控情况。由图可见,“新获取的位移监测数据”全部位于本发明方法构造的坝顶位移置信区间内,而有约70%的监测数据位于HST统计模型得到的坝顶位移置信区间之外。由于本实例分析的重力坝一直处于正常工作状态,因此,采用HST统计模型进行该重力坝坝顶位移监控,会出现高达70%的误判率,而本发明方法构建的监控模型未出现误判,能正确反映该重力坝的变形性态。
表1某重力坝9#坝段监测仪器埋设位置考证表
表2某重力坝9#坝段坝与地基材料的设计参数
表3坝与地基弹性模量和库水位的正交试验设计结果
表4坝与地基弹性模量的反分析结果
表5坝体温度计测值序列的主成分分析结果
主成分 | 特征值 | 方差贡献率(%) | 累计贡献率(%) |
PC<sub>1</sub> | 39.91 | 93.16 | 93.16 |
PC<sub>2</sub> | 2.68 | 6.26 | 99.42 |
PC<sub>3</sub> | 0.23 | 0.53 | 99.94 |
PC<sub>4</sub> | 0.02 | 0.06 | 100 |
本发明不局限于上述可选实施方式,任何人在本发明的启示下都可得出其他各种形式的产品,但不论在其形状或结构上作任何变化,凡是落入本发明权利要求界定范围内的技术方案,均落在本发明的保护范围之内。
Claims (5)
1.一种基于分离式混合建模技术的重力坝坝顶位移监控方法,其特征在于,包括以下步骤:
S1.采用EMD技术将重力坝坝顶位移监测数据序列分解为一系列不同频率的固有模态函数分量,其中频率最低且具有单调性的固有模态函数分量为重力坝坝顶位移监测数据序列中的时效趋势项,将其余固有模态函数分量进行叠加重构即可得到重力坝坝顶位移监测数据序列的周期变化项;
S2.根据EMD提取的时效趋势项,采用复合指数函数构建重力坝坝顶位移的时效分量模型,即
F3[θ(t)]=h*exp[kθ(t)j]+q (1)
式中,θ(t)=(t-t0)/100,t为时间,t0为监测数据序列的起始时间;h、k、j、q为模型参数;
S3.采用EMD分解、重构得到的重力坝坝顶位移监测数据序列的周期变化项是库水位和坝体温度变化共同作用的结果,根据两个时刻坝体温度测点的测值相等或温度测值之差小于0.2℃,从坝顶位移监测数据序列的周期变化项中分离出库水位变化导致的位移变化值,即坝顶位移监测数据序列的水压效应量;
S4.根据坝顶位移监测数据序列的水压效应量,构建坝顶位移的水压分量模型;
S5.利用建立的时效和水压分量模型从重力坝坝顶位移监测数据序列中分解出温度效应量,构建坝顶位移的温度分量模型;
S6.叠加步骤S2、S4和S5构建的时效、水压和温度分量模型,最终得到重力坝坝顶位移的监控模型;
S7.计算坝顶位移实测值和监控模型估计值的差值,采用核函数拟合模型误差的概率密度函数,根据坝顶位移的模型估计值和模型误差的概率密度函数,构造坝顶位移的置信区间,对坝顶位移是否异常进行监控,当新获取的坝顶位移监测数据位于置信区间之外时,则发出数据异常预警。
2.根据权利要求1所述的一种基于分离式混合建模技术的重力坝坝顶位移监控方法,其特征在于,步骤S2中h、k、j、q模型参数,根据重力坝坝顶位移监测数据序列的时效趋势项,采用1stOpt15PRO软件确定。
3.根据权利要求1所述的一种基于分离式混合建模技术的重力坝坝顶位移监控方法,其特征在于,所述步骤S4中构建坝顶位移的水压分量模型包括以下步骤:
1)根据坝体结构轮廓尺寸、材料分区和坝址地形地质条件等资料,建立重力坝坝段的三维有限元模型,根据坝顶位移测点位置在有限元模型中布置相应结点,以便于提取位移计算结果;
2)确定需要反演分析的坝体混凝土和地基岩体弹性模量,以及各弹性模量的取值区间,根据坝与地基材料弹性模量的取值区间和库水位变化范围,确定弹性模量和库水位的正交试验设计水平数和相应数值,采用正交试验设计对弹性模量和库水位的各水平值进行组合;
3)采用第1)步建立的重力坝坝段有限元模型,按照弹性模量和库水位的不同组合,计算相应的坝顶位移值,根据弹性模量和库水位的不同组合,以及相应的坝顶位移计算值,采用1stOpt15PRO软件确定坝顶位移水压分量模型中的系数a、bi、ci、di、ei,见式(2);
式中,H(t)为某个时刻库水位对应的上游水深;xi(i=1,2,…,n)代表坝体混凝土和地基岩体的弹性模量;
4)根据S3步得到的坝顶位移监测数据序列的水压效应量及对应的上游水深变化值,采用1stOpt15PRO软件对式(2)中的xi(i=1,2,...,n)进行反分析;
5)将第4)步得到的坝体混凝土和地基岩体弹性模量的反分析值代入式(2),即可建立以上游水深为变量的坝顶位移水压分量模型。
4.根据权利要求1所述的一种基于分离式混合建模技术的重力坝坝顶位移监控方法,其特征在于,所述步骤S5中构建坝顶位移的温度分量模型包括以下步骤:
1)计算坝体温度测点的测值与坝顶位移测值的Spearman相关系数,根据Spearman相关系数绝对值大于0.65的条件,筛选出与坝顶位移测值相关性较强的坝体温度测点,用于构建坝顶位移的温度分量模型;
2)对筛选出的坝体温度测点的测值序列进行主成分分析,并根据主成分对原测值方差的累计贡献率大于95%的条件,确定用于构建温度分量模型的主成分个数;
3)采用坝体温度测值的主成分构建坝顶位移温度分量模型,即
式中,Pi(t)(i=1,2,…,m)为坝体温度测值的主成分,m为对原测值方差的累计贡献率大于95%的主成分个数;fi(i=1,2,…,m)和g为待定系数,根据从坝顶位移监测数据序列中分解出的温度效应量,采用最小二乘法确定。
5.根据权利要求1所述的一种基于分离式混合建模技术的重力坝坝顶位移监控方法,其特征在于,步骤S7中采用MATLAB软件中的核函数拟合模型误差的概率密度函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011463759.5A CN112507591B (zh) | 2020-12-11 | 2020-12-11 | 一种基于分离式混合建模技术的重力坝坝顶位移监控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011463759.5A CN112507591B (zh) | 2020-12-11 | 2020-12-11 | 一种基于分离式混合建模技术的重力坝坝顶位移监控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112507591A true CN112507591A (zh) | 2021-03-16 |
CN112507591B CN112507591B (zh) | 2022-11-18 |
Family
ID=74972594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011463759.5A Active CN112507591B (zh) | 2020-12-11 | 2020-12-11 | 一种基于分离式混合建模技术的重力坝坝顶位移监控方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112507591B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113742939A (zh) * | 2021-09-15 | 2021-12-03 | 水利部南京水利水文自动化研究所 | 一种运行期混凝土坝效应量温度分量模型的构建方法 |
CN116842339A (zh) * | 2023-06-20 | 2023-10-03 | 北京讯腾智慧科技股份有限公司 | 一种大坝数据分析方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080082470A1 (en) * | 2006-09-29 | 2008-04-03 | Ehsan Sobhani Tehrani | Infrastructure health monitoring and analysis |
CN106934208A (zh) * | 2017-01-05 | 2017-07-07 | 中国电建集团华东勘测设计研究院有限公司 | 一种大坝异常监测数据自动识别方法 |
US20180275303A1 (en) * | 2017-03-24 | 2018-09-27 | Xin ZHAN | 4D Time Shift and Amplitude Joint Inversion for Obtaining Quantitative Saturation and Pressure Separation |
CN109117540A (zh) * | 2018-08-02 | 2019-01-01 | 三峡大学 | 一种解决大坝混凝土力学参数反演不唯一性的概率统计分析方法 |
CN111191191A (zh) * | 2019-12-26 | 2020-05-22 | 南昌大学 | 一种精准预测混凝土坝变形效应的组合模型的构建方法 |
CN111259590A (zh) * | 2020-01-17 | 2020-06-09 | 南昌大学 | 一种混凝土坝变形安全监控模型的构建方法 |
CN111551147A (zh) * | 2020-06-09 | 2020-08-18 | 福州大学 | 基于gnss与测量机器人融合的拱坝表面变形监测系统 |
CN111695177A (zh) * | 2020-05-14 | 2020-09-22 | 河海大学 | 基于位移场监测数据的大坝力学参数随机反演方法及系统 |
-
2020
- 2020-12-11 CN CN202011463759.5A patent/CN112507591B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080082470A1 (en) * | 2006-09-29 | 2008-04-03 | Ehsan Sobhani Tehrani | Infrastructure health monitoring and analysis |
CN106934208A (zh) * | 2017-01-05 | 2017-07-07 | 中国电建集团华东勘测设计研究院有限公司 | 一种大坝异常监测数据自动识别方法 |
US20180275303A1 (en) * | 2017-03-24 | 2018-09-27 | Xin ZHAN | 4D Time Shift and Amplitude Joint Inversion for Obtaining Quantitative Saturation and Pressure Separation |
CN109117540A (zh) * | 2018-08-02 | 2019-01-01 | 三峡大学 | 一种解决大坝混凝土力学参数反演不唯一性的概率统计分析方法 |
CN111191191A (zh) * | 2019-12-26 | 2020-05-22 | 南昌大学 | 一种精准预测混凝土坝变形效应的组合模型的构建方法 |
CN111259590A (zh) * | 2020-01-17 | 2020-06-09 | 南昌大学 | 一种混凝土坝变形安全监控模型的构建方法 |
CN111695177A (zh) * | 2020-05-14 | 2020-09-22 | 河海大学 | 基于位移场监测数据的大坝力学参数随机反演方法及系统 |
CN111551147A (zh) * | 2020-06-09 | 2020-08-18 | 福州大学 | 基于gnss与测量机器人融合的拱坝表面变形监测系统 |
Non-Patent Citations (7)
Title |
---|
CHUAN YIN 等: "Separate modeling technique for deformation monitoring of concrete dams", 《STRUCTURAL HEALTH MONITORING》 * |
M. TATIN 等: "Statistical modelling of thermal displacements for concrete dams: Influence of water temperature profile and dam thickness profile", 《ENGINEERING STRUCTURES》 * |
姜振翔 等: "基于小波分解和支持向量机的大坝位移监控模型", 《长江科学院院报》 * |
易恒如 等: "基于挠度公式的重力坝水平位移监测模型水压分量的校验", 《大坝与安全》 * |
李端有等: "混凝土拱坝多测点确定性位移监控模型研究", 《水利学报》 * |
田紫圆 等: "基于机器学习的重力坝变形监测统计模型及应用", 《水利规划与设计》 * |
祖安君: "基于统计模型的低水头混凝土重力坝垂直位移影响因素分析", 《水利水电技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113742939A (zh) * | 2021-09-15 | 2021-12-03 | 水利部南京水利水文自动化研究所 | 一种运行期混凝土坝效应量温度分量模型的构建方法 |
CN116842339A (zh) * | 2023-06-20 | 2023-10-03 | 北京讯腾智慧科技股份有限公司 | 一种大坝数据分析方法及装置 |
CN116842339B (zh) * | 2023-06-20 | 2024-03-12 | 北京讯腾智慧科技股份有限公司 | 一种大坝数据分析方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112507591B (zh) | 2022-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shah et al. | Modelling the effects of spatial variability in rainfall on catchment response. 1. Formulation and calibration of a stochastic rainfall field model | |
CN112507591B (zh) | 一种基于分离式混合建模技术的重力坝坝顶位移监控方法 | |
CN106529055A (zh) | 一种基于应变模态振型相关性的模型修正方法 | |
Vergara et al. | Estimating a-priori kinematic wave model parameters based on regionalization for flash flood forecasting in the Conterminous United States | |
CN112711899B (zh) | 一种蒸发波导高度的融合预测方法 | |
CN114547951B (zh) | 一种基于数据同化的大坝状态预测方法及系统 | |
CN104048676B (zh) | 基于改进粒子滤波的mems陀螺随机误差补偿方法 | |
CN109858064B (zh) | 一种基于原型振动响应的大坝及地基弹模动力反演方法 | |
CN112464584B (zh) | 自由表面流的水位和流量推求方法 | |
CN111191191A (zh) | 一种精准预测混凝土坝变形效应的组合模型的构建方法 | |
Ray et al. | Bayesian calibration of the Community Land Model using surrogates | |
CN113532397B (zh) | 一种基于扩展算法的分区域高程异常拟合方法 | |
CN104091065A (zh) | 一种求解浅水问题模拟间断水流数值的方法 | |
CN103353295B (zh) | 一种精确预测大坝坝体垂直变形量的方法 | |
CN109190291A (zh) | 获取动力触探锤击数修正系数的方法 | |
CN105404709B (zh) | 基于复杂网络的堤坝健康监测敏感测点分析方法 | |
CN112100711B (zh) | 一种基于arima和pso-elm的混凝土坝变形组合预报模型构建方法 | |
CN105893329B (zh) | 一种基于月尺度的潮位资料一致性修正方法 | |
Pham et al. | Spatiotemporal variations of extreme sea levels around the South China Sea: assessing the influence of tropical cyclones, monsoons and major climate modes | |
CN115455867A (zh) | 基于回归分析的坝区流态的推求方法 | |
CN111815769B (zh) | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 | |
Xu et al. | Improving the real-time probabilistic channel flood forecasting by incorporating the uncertainty of inflow using the particle filter | |
CN116413187B (zh) | 一种基于毛管压力资料的储层渗透率预测方法及系统 | |
Zhang et al. | A river stage correction approach using Fourier series | |
CN116660974B (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 |