CN111811468B - 高心墙堆石坝沉降变形的预测方法 - Google Patents
高心墙堆石坝沉降变形的预测方法 Download PDFInfo
- Publication number
- CN111811468B CN111811468B CN202010922824.XA CN202010922824A CN111811468B CN 111811468 B CN111811468 B CN 111811468B CN 202010922824 A CN202010922824 A CN 202010922824A CN 111811468 B CN111811468 B CN 111811468B
- Authority
- CN
- China
- Prior art keywords
- deformation
- dam
- dam body
- high core
- fill
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种高心墙堆石坝沉降变形的预测方法,其包括:以高心墙堆石坝的建基面作为参考面,依据坝体内已布设的沉降变形监测点及其位置确定各监测点到参考面的距离z i ;判断拟进行变形预测的高心墙堆石坝处于施工期还是运行期;当高心墙堆石坝处于施工期时,依据沉降变形物理成因、组成以及距离z i ,分别计算坝体沉降变形、牵连变形和蠕变变形,并累加得到施工期坝体沉降变形;当高心墙堆石坝处于运行期时,根据监测点与浸润线相对位置和距离z i ,采用对应的总沉降变形计算方法进行总沉降变形的计算。本发明所述沉降变形预测方法更加符合工程实际,能够为高心墙堆石坝的变形调控和安全管控提供更加准确的数据支撑。
Description
技术领域
本发明涉及大坝安全监测技术领域,具体涉及一种高心墙堆石坝沉降变形的预测方法。
背景技术
高心墙堆石坝由于具有对地基条件适应性好、能就地取材和充分利用建筑物开挖料、投资省、施工方法简单以及抗震性能好等优点,在坝工建设中被广泛采用。目前,在我国西部水电开发中,将有一批高心墙堆石坝工程投入建设和运行,如长河坝(坝高240m)、两河口(坝高295m)、古水(坝高305m)、双江口(坝高314m)等,这些高心墙堆石坝变形机理复杂,变形预测难度大,因不协调变形产生的变形破坏风险高,因此,有必要从沉降变形机理与物理成因入手,研究沉降变形的高精度预测方法,为高堆石坝变形调控和安全管控提供技术支撑。
目前心墙堆石坝变形预测一般采用传统的数学模型法,这种方法是通过建立描述效应量与因变量之间的统计学模型,确定它们之间的定量关系。传统的统计回归模型进行测点变形预测,仅仅是采用多元线性回归模型进行系数拟合,未考虑测点变形的物理成因,极易出现拟合精度低、鲁棒性差、适用性低等情况,而且对于时效分量的求解,按照上述形式构建回归方程很有可能造成回归模型起始点偏差较大的现象。
发明内容
针对现有技术中的上述不足,本发明提供的高心墙堆石坝沉降变形的预测方法解决了现有技术未考虑变形监测点变形的物理成因致使预测精度差的问题。
为了达到上述发明目的,本发明采用的技术方案为:
提供一种高心墙堆石坝沉降变形的预测方法,其包括:
S1、选取高心墙堆石坝的建基面作为参考面,根据坝体内已布设的沉降变形监测点及其位置确定各监测点到参考面的距离;
S2、判断高心墙堆石坝为施工期或运行期,若为施工期,则进入步骤S3,若为运行期,则进入步骤S7;
S3、根据坝体材料的邓肯-张模型参数及监测点到参考面的距离,计算施工期坝体压缩产生的沉降变形S di :
其中,H i 为第i个变形监测点到坝表面的垂直距离,单位为m;z i 为第i个变形监测点到参考面的距离,单位为m;n为坝体材料的邓肯-张模型参数;B为坝体变形系数;
S4、根据坝基覆盖层深度、坝基材料的邓肯-张模型参数及变形监测点到参考面的距离,计算施工期坝体压缩产生的牵连变形S li :
其中,H l 为坝基覆盖层深度,单位为m;n l 为坝基材料的邓肯-张模型参数;B l 为坝基变形系数;
式中,K p 为蠕变变形系数;α为坝体堆石料蠕变参数;t j 为坝体开始填筑至变形监测点获得初始测值的时间,单位为天;t c 为坝体开始填筑至进行变形预测的时间,单位为天;t p 为0至t j 中间过程中选定的时间节点,单位为天;p为时间节点变量,1≤p≤j;j为选取的时间节点个数,视精度和预测时间长短而定;e为自然数;
S7、判断变形监测点i到参考面的距离是否小于变形监测点i处对应的浸润线高程H wi ,若是进入步骤S8,否则进入步骤S9;
S8、根据变形监测点i到参考面的距离和变形预测点i到坝体表面的垂直距离,计算运行期坝体任意监测点的总沉降变形S zi :
S9、根据变形预测点i处对应的浸润线高程和变形预测点i到坝体表面的垂直距离,计算运行期坝体任意监测点处的总沉降变形S zi :
本发明的有益效果为:本方案在进行沉降变形预测时,首先获取预测点的相对位置关系,并基于相对位置关系可以得到监测点所在位置的应力应变状态,最终通过应力应变关系得到每个预测点的沉降变形;由于本方案在预测时充分考虑高心墙堆石坝上述的物理成因,相对现有技术中数学模拟模型而言,综合考虑了预测点的位置关系和应力应变状态,使得本方案的沉降变形预测模型能较好地模拟和反映测点变形趋势,显著提升变形预测精度。
附图说明
图1为高心墙堆石坝沉降变形的预测方法的流程图。
图2为坝内任意点与坝高、覆盖层深度的位置关系图。
图3为坝体任意点沉降变形受坝体压缩和坝基压缩的影响。
图4为5-3FLAC3d数值计算算例中沉降变形随坝高的变化。
图5为坝填筑高度与时间的对应关系。
图6为施工期VE1-18测点沉降测值和计算值。
图7为施工期VE2-37测点沉降测值和计算值。
图8为运行期VE1-25 实测值与计算值对比。
图9为LD59 实测值与计算值对比。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
参考图1,图1示出了高心墙堆石坝沉降变形的预测方法的流程图;如图1所示,该方法S包括S1至S9。
在步骤S1中,选取高心墙堆石坝的建基面作为参考面,根据坝体内已布设的沉降变形监测点及其位置确定各监测点到参考面的距离;高心墙堆石坝的结构示意参考图2。
在步骤S2中,判断高心墙堆石坝为施工期或运行期,若为施工期,则进入步骤S3,若为运行期,则进入步骤S7;
在步骤S3中,根据坝体材料的邓肯-张模型参数及监测点到参考面的距离,计算施工期坝体压缩产生的沉降变形S di :
其中,H i 为第i个变形监测点到坝表面的垂直距离,单位为m;z i 为第i个变形监测点到参考面的距离,单位为m;n为坝体材料的邓肯-张模型参数;B为坝体变形系数;
在本发明的一个实施例中,所述施工期坝体压缩产生的变形S di 的构建方法包括:
S31、构建坝体内任意一点的压缩变形的初始模型:
式中,α d 为坝体压应力修正系数;γ为填筑料的容重;E s 为填筑料的压缩模量;z 1,z 2均为积分变量;
S32、采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E s :
S34、将公式(1-3)和公式(1-4)带入公式(1-2),可得:
在步骤S4中,根据坝基覆盖层深度、坝基材料的邓肯-张模型参数及变形监测点到参考面的距离,计算施工期坝体压缩产生的牵连变形S li :
其中,H l 为坝基覆盖层深度,单位为m;n l 为坝基材料的邓肯-张模型参数;B l 为坝基变形系数;
在本发明的一个实施例中,所述施工期坝体压缩产生的牵连变形S li 的构建方法包括:
S41、当不考虑附加应力在坝基内的应力扩散衰减时,有覆盖层变形引起的坝体牵连变形为:
式中,E ls 为坝基压缩模量;a l 为坝基内压应力修正系数;采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E ls :
S43、将公式(2-3)和公式(2-4)带入公式(2-2),可得:
实施时,本方案优选所述坝基应力相关系数B l 的计算公式为:
其中,a l 为坝基压应力修正系数;β l 、η l 分别为坝基内的大主应力及小主应力修正系数;γ为填筑料的容重;k el 、R fl 分别为坝基材料邓肯-张模型参数;pa为大气压强。
式中,K p 为蠕变变形系数;α为坝体堆石料蠕变参数;t j 为坝体开始填筑至变形监测点获得初始测值的时间,单位为天;t c 为坝体开始填筑至进行变形预测的时间,单位为天;t p 为0至t j 中间过程中选定的时间节点,单位为天;p为时间节点变量,1≤p≤j;j为选取的时间节点个数,视精度和预测时间长短而定;e为自然数。
在三轴蠕变试验中,轴向蠕变变形与体积蠕变变形和剪切蠕变变形之间满足如下关系:
轴向蠕变计算公式为:
假设中主应力可近似表示为:
令
那么,式(5-3)中一点的最终轴向蠕变可近似表示如式(5-5)所示的形式:
根据式(5-9),合并同类项,则任意层的沉降蠕变变形量可如式(5-10)计算:
进一步合并式(5-10)中的同类项后可写成如式(5-11)所示的一般形式。
式中,为蠕变变形系数;为坝体堆石料蠕变参数;t j 为坝体开始填筑至变形监测点获得初始测值的时间,单位为天;t c 为坝体开始填筑至进行变形预测的时间,单位为天;t p 为0至t j 中间过程中选定的时间节点,单位为天;p为时间节点变量,;为选取的时间节点个数,视精度和预测时间长短而定;e为自然数;
在步骤S7中,判断变形监测点i到参考面的距离是否小于变形监测点i处对应的浸润线高程H wi ,若是进入步骤S8,否则进入步骤S9;
在步骤S8中,根据变形监测点i到参考面的距离和变形预测点i到坝体表面的垂直距离,计算运行期坝体任意监测点的总沉降变形S zi :
在步骤S9中,根据变形预测点i处对应的浸润线高程和变形预测点i到坝体表面的垂直距离,计算运行期坝体任意监测点处的总沉降变形S zi :
。
由图3显示的变形规律可知,因坝体压缩产生的沉降变形随计算点所在的坝高位置表现为先增后减的规律,而坝基压缩对坝体变形的影响表现为越远离大坝建基面,则因坝基覆盖层压缩产生的坝体牵连变形越小。为进一步验证计算成果的合理性,对上述算例采用FLAC3d进行了对比分析。结果表明,FLAC3d计算的坝体沉降变形随坝高的分布规律(见图4)与图3解析计算的分布规律一致,表明本方案在施工期引入坝体沉降变形物理力学成因,能够提高预测精度。
下面结合结合瀑布沟高心墙堆石坝,对本方案的预测方法的效果进行说明:
选取瀑布沟高心墙堆石坝的坝底作为参考面,选取瀑布沟高心墙堆石坝0+240m监测断面进行验证,在施工期阶段,选取变形监测点VE2-37在大坝施工期2008年02月23日~2009年01月05日期间及变形监测点VE1-18在大坝施工期2009年02月20日~2009年08月25日期间的相对位置、大坝的材料参数及这两个变形监测点的实测数据参考表1至表3:
表1 VE1-18、VE2-37参数表(固定参数)
表2 VE1-18参数表(变化参数
表3 VE2-37参数表(变化参数)
基于变形监测点VE1-18和VE2-37的相对位置及大坝的材料参数,采用本方案的预测方法进行预测,变形监测点VE1-18和VE2-37的预测数据和实测数据分别参见图6和图7。
通过图6和图7可以得出,变形监测点计算值与实测值的平均相对误差均在10%以内,精度较高,可以采用本方案的方法得到的施工期预测数据与实测值历时过程对比示例,本方案从监测点变形序列的发展过程来看,本方案的预测方法能较好地描述测点变形的历时过程,预测效果较好。
在运行期间,选取变形监测点VE1-25在大坝施工期2009年10月11日~2012年10月30日期间及变形监测点LD59在大坝施工期2009年11月06日~2010年10月05日期间的水位高及材料湿化变形参数,变形监测点到坝表面的高度及两个变形监测点的实测数据表4至表6:
表4 VE1-25、LD59参数表(固定参数)表中无单位
表5 VE1-25参数表(变化参数)
表6 LD59参数表(变化参数)
基于获取的数据采用本方案的预测方法进行预测,变形监测点VE1-25和LD59和预测数据参考图8和图9。
通过图8和图9可以得出,采用本方案的预测方法得到的预测值与实测值的相对误差在10%以内,精度较高,满足工程应用需求。
Claims (5)
1.高心墙堆石坝沉降变形的预测方法,其特征在于,包括:
S1、选取高心墙堆石坝的建基面作为参考面,根据坝体内已布设的沉降变形监测点及其位置确定各监测点到参考面的距离;
S2、判断高心墙堆石坝为施工期或运行期,若为施工期,则进入步骤S3,若为运行期,则进入步骤S7;
S3、根据坝体材料的邓肯-张模型参数及监测点到参考面的距离,计算施工期坝体压缩产生的沉降变形S di :
其中,H i 为第i个变形监测点到坝表面的垂直距离,单位为m;z i 为第i个变形监测点到参考面的距离,单位为m;n为坝体材料的邓肯-张模型参数;B为坝体变形系数;
S4、根据坝基覆盖层深度、坝基材料的邓肯-张模型参数及变形监测点到参考面的距离,计算施工期坝体压缩产生的牵连变形S li :
其中,H l 为坝基覆盖层深度,单位为m;n l 为坝基材料的邓肯-张模型参数;B l 为坝基变形系数;
式中,K p 为蠕变变形系数;α为坝体堆石料蠕变参数;t j 为坝体开始填筑至变形监测点获得初始测值的时间,单位为天;t c 为坝体开始填筑至进行变形预测的时间,单位为天;t p 为0至t j 中间过程中选定的时间节点,单位为天;p为时间节点变量,1≤p≤j;j为选取的时间节点个数;e为自然数;
S7、判断变形监测点i到参考面的距离是否小于变形监测点i处对应的浸润线高程H wi ,若是进入步骤S8,否则进入步骤S9;
S8、根据变形监测点i到参考面的距离和变形预测点i到坝体表面的垂直距离,计算运行期坝体任意监测点的总沉降变形S zi :
S9、根据变形预测点i处对应的浸润线高程和变形预测点i到坝体表面的垂直距离,计算运行期坝体任意监测点处的总沉降变形S zi :
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010922824.XA CN111811468B (zh) | 2020-09-04 | 2020-09-04 | 高心墙堆石坝沉降变形的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010922824.XA CN111811468B (zh) | 2020-09-04 | 2020-09-04 | 高心墙堆石坝沉降变形的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111811468A CN111811468A (zh) | 2020-10-23 |
CN111811468B true CN111811468B (zh) | 2020-12-08 |
Family
ID=72860638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010922824.XA Active CN111811468B (zh) | 2020-09-04 | 2020-09-04 | 高心墙堆石坝沉降变形的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111811468B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113449437B (zh) * | 2021-07-26 | 2022-04-12 | 青岛理工大学 | 一种黏土斜墙土石坝渗漏预警方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109917450A (zh) * | 2019-04-10 | 2019-06-21 | 长江水利委员会长江科学院 | 一种全断面网络化布设的土石坝位移监测系统及方法 |
CN110987618A (zh) * | 2019-12-13 | 2020-04-10 | 中国电建集团中南勘测设计研究院有限公司 | 一种堆石料试验装置及其方法 |
CN111210519A (zh) * | 2019-12-25 | 2020-05-29 | 四川大学 | 一种高心墙堆石坝变形参数反演方法 |
CN111475981A (zh) * | 2020-04-22 | 2020-07-31 | 水发规划设计有限公司 | 一种水库土石坝渗透变形和稳定性的数值模拟方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN203249668U (zh) * | 2013-05-07 | 2013-10-23 | 山东大学 | 一种坝体浸润线监测模型试验装置 |
CN104196010A (zh) * | 2014-08-14 | 2014-12-10 | 中国电建集团西北勘测设计研究院有限公司 | 一种面板堆石坝强制预沉降施工方法 |
JP6696290B2 (ja) * | 2016-04-28 | 2020-05-20 | 株式会社大林組 | ロックフィルダムにおけるロックゾーンの品質管理方法 |
-
2020
- 2020-09-04 CN CN202010922824.XA patent/CN111811468B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109917450A (zh) * | 2019-04-10 | 2019-06-21 | 长江水利委员会长江科学院 | 一种全断面网络化布设的土石坝位移监测系统及方法 |
CN110987618A (zh) * | 2019-12-13 | 2020-04-10 | 中国电建集团中南勘测设计研究院有限公司 | 一种堆石料试验装置及其方法 |
CN111210519A (zh) * | 2019-12-25 | 2020-05-29 | 四川大学 | 一种高心墙堆石坝变形参数反演方法 |
CN111475981A (zh) * | 2020-04-22 | 2020-07-31 | 水发规划设计有限公司 | 一种水库土石坝渗透变形和稳定性的数值模拟方法及装置 |
Non-Patent Citations (2)
Title |
---|
"面板堆石坝堆石体湿化变形分析方法研究";方国宝;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20071215(第06期);正文第3.3.3小节,第3.4节 * |
"高心墙堆石坝施工期变形时空预测模型研究";卢祥 等;《工程科学与技术》;20170731;第49卷(第4期);正文第1.2.1小节、图1,2 * |
Also Published As
Publication number | Publication date |
---|---|
CN111811468A (zh) | 2020-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016185B (zh) | 一种冰碛湖溃决洪水峰值流量预测的计算方法 | |
CN108487904B (zh) | 一种基于图版的消除末端效应的相渗曲线校正方法 | |
CN107515025B (zh) | 古城墙病害等级评估方法 | |
Göğüş et al. | Broad-crested weirs with rectangular compound cross sections | |
CN111811468B (zh) | 高心墙堆石坝沉降变形的预测方法 | |
CN105155425A (zh) | 满堂支架流体仿真模拟预压施工方法 | |
Shi et al. | A statistical model of deformation during the construction of a concrete face rockfill dam | |
CN110969370B (zh) | 建筑结构构件的质量风险分析方法 | |
CN104111137B (zh) | 基于三次多项式模型的锚杆承载力预测方法 | |
CN114493702A (zh) | 一种基于有限元分析及地理信息的农村房产损失评估方法 | |
CN111950065A (zh) | 高心墙堆石坝施工期水平变形的预测方法 | |
Wang et al. | Dynamic tensile test of mass concrete with Shapai Dam cores | |
De Vos et al. | The current draft of the application rules of Eurocode 7 in Belgium for the design of pile foundations | |
CN110287634A (zh) | 一种基于体积力施加的坝肩变形模拟方法 | |
Jostad et al. | Cyclic capacity of monopiles in sand under Partially drained conditions: a numerical approach | |
CN110197015A (zh) | 一种坝基预应力锚索有效拉应力测定方法 | |
CN112861207B (zh) | 复合地层的沉降预测的方法、设备及计算机存储介质 | |
CN114662341A (zh) | 岩体临界滑动面极限分析方法 | |
CN111898176B (zh) | 一种成层地基逐渐加荷条件下固结度的计算方法 | |
ERGUN et al. | COMPARISON OF GEOMETRICAL IMPERFECTION DEFINITIONS ON ENCASED COMPOSITE COLUMNS | |
RU2457455C2 (ru) | Способ определения герметичности запорных арматур нефтепровода | |
CN110020792A (zh) | 基于洪水洪峰-洪量联合的高堆石坝施工度汛风险预测方法 | |
Wessels et al. | Flow-gauging structures in South African rivers part 2: calibration | |
Mozer et al. | Time-dependent behavior of concrete beams | |
Balmforth | Flow over side weirs |
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 |