CN111861107A - 一种基于强度折减法的滑坡动态稳定性评价方法 - Google Patents
一种基于强度折减法的滑坡动态稳定性评价方法 Download PDFInfo
- Publication number
- CN111861107A CN111861107A CN202010513680.2A CN202010513680A CN111861107A CN 111861107 A CN111861107 A CN 111861107A CN 202010513680 A CN202010513680 A CN 202010513680A CN 111861107 A CN111861107 A CN 111861107A
- Authority
- CN
- China
- Prior art keywords
- landslide
- creep
- reduction
- model
- soil
- 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
- 238000000034 method Methods 0.000 title claims abstract description 67
- 230000009467 reduction Effects 0.000 title claims abstract description 67
- 238000011156 evaluation Methods 0.000 title claims abstract description 20
- 239000002689 soil Substances 0.000 claims abstract description 62
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 49
- 230000007774 longterm Effects 0.000 claims abstract description 46
- 238000004458 analytical method Methods 0.000 claims abstract description 40
- 230000008878 coupling Effects 0.000 claims abstract description 21
- 238000010168 coupling process Methods 0.000 claims abstract description 21
- 238000005859 coupling reaction Methods 0.000 claims abstract description 21
- 239000012530 fluid Substances 0.000 claims abstract description 20
- 239000007787 solid Substances 0.000 claims abstract description 15
- 238000004088 simulation Methods 0.000 claims abstract description 12
- 230000000737 periodic effect Effects 0.000 claims abstract description 5
- 238000012544 monitoring process Methods 0.000 claims description 41
- 238000006073 displacement reaction Methods 0.000 claims description 30
- 238000004364 calculation method Methods 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 18
- 238000012360 testing method Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 11
- 239000011435 rock Substances 0.000 claims description 6
- 230000009471 action Effects 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 4
- 238000007596 consolidation process Methods 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 230000000630 rising effect Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 abstract description 5
- 230000002265 prevention Effects 0.000 abstract description 4
- 235000015220 hamburgers Nutrition 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000035772 mutation Effects 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 230000005476 size effect Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Images
Classifications
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Fluid Mechanics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Primary Health Care (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明提供一种基于强度折减法的滑坡动态稳定性评价方法，步骤如下：根据滑坡所处的库区水位周期调度情况建立概化模型，利用滑坡岩土体的参数数据计算力学分析时标和流体分析时标，比较力学分析时标和流体分析时标的大小，确定流固耦合模拟方法；建立滑带土的蠕变本构模型，识别滑带土的蠕变参数；计算滑带土的长期抗剪强度；构建滑坡模型，计算滑坡的应力场、应变场和渗流场；计算滑坡稳定系数。本发明提供的方法揭示了水动力条件下的滑坡变形特征与演化机理，实现了对库区蠕变型滑坡稳定性的动态评价，为同类型滑坡的预测预警与防治提供了基础。
Description
技术领域
本发明涉及地质灾害防治领域，尤其涉及一种基于强度折减法的滑坡动态稳定性评价方法。
背景技术
自一系列国家水电发展规划实施以来，中国西南地区兴建了一大批水利水电工程，这些工程的建设和运行显著改变了库区水动力条件，造成了大量水库滑坡的复活和产生，严重威胁到人类的生命财产安全。目前评价滑坡稳定性系数最常用的规范法中，难以考虑到水位实时变化与滑坡岩土体本构关系对其的影响，这种静态的评价方法对滑坡当前与长期安全状态难以作出准确判断。
随着计算机技术的发展，自R.W.Clough在1960年首次提出“有限元”以来，数值模拟技术在岩土工程分析中取得了巨大发展，其中由Itasca公司研发的FLAC/FLAC3D被广泛应用于土木、水利、地质等领域。FLAC/FLAC3D采用强度折减法求解滑坡稳定性系数，与传统刚体极限平衡法相比，它不需要事先指定滑面，也可以对具有复杂三维地形的滑坡进行求解，同时可以将岩土体本构关系与滑坡稳定性联系起来。然而，目前大多数分析均采用了M-C模型，没有考虑到滑坡岩土体材料蠕变特性导致的其临界应变状态改变和强度劣化；同时由于之前版本计算效率较低，大多数研究采用将GEOSLOPE软件得到的渗流场以“watertable”命令导入来实现对滑坡地下水位的模拟，没有实现真耦合。这两个不足导致了目前对水动力蠕变型滑坡的演化机理及稳定性评价等研究存在偏差；同时由于规范法中求解滑坡稳定性系数的公式中不包含时间指标，因此难以评价在在不同内外动力作用下滑坡的动态稳定性。
发明内容
有鉴于此，本发明提供了一种基于强度折减法的滑坡动态稳定性评价方法，该方法揭示了水动力条件下的滑坡变形特征与演化机理，实现了对库区蠕变型滑坡稳定性的动态评价，为同类型滑坡的预测预警与防治提供了基础。
本发明提供一种基于强度折减法的滑坡动态稳定性评价方法，包括以下步骤：
步骤S1，根据滑坡所处的库区水位周期调度情况建立概化模型，利用滑坡岩土体的参数数据计算力学分析时标和流体分析时标，比较力学分析时标和流体分析时标的大小，确定流固耦合模拟方法；
步骤S2，建立滑带土的蠕变本构模型，识别滑带土的蠕变参数；
步骤S3，计算滑带土的长期抗剪强度；
步骤S4，构建滑坡模型，在滑坡涉水部位施加步骤S1的概化模型，利用步骤S2的蠕变参数和步骤S3的长期抗剪强度设置滑坡岩土体本构模型，根据步骤S1的流固耦合模拟方法设置耦合方式，在滑坡模型上设置与现场GNSS监测对应的监测点，计算滑坡的应力场、应变场和渗流场；
步骤S5，对滑带土长期抗剪强度指标进行线性增加折减，每折减一次，重复计算滑坡的应力场、应变场和渗流场，并记录当前监测点的位移值，绘制位移-折减系数曲线，直至曲线出现突变，停止折减，记录最后一次折减系数Fmax，对最后一次折减系数Fmax进行线性减小折减，每折减一次，重复计算滑坡的应力场、应变场和渗流场，记录当前监测点的位移值，将该当前监测点的位移值插入到位移-折减系数曲线上，直至插值点前未发生位移突变，记录最后一次插值点对应的折减系数Fs，该折减系数Fs即为滑坡当前状态下的稳定系数。
进一步地，步骤S1中，根据库水位上升、高水位保持、库水位下降、低水位保持四个阶段对滑坡进行概化，建立水位与时间的函数表达式，记为H＝f(t)，H为水位高程，t为天数；若低水位保持时间较短，可简化为三阶段。
进一步地，步骤S1中，力学分析时标的计算公式为：
式中，tm为力学分析时标；G为剪切模量；Lc m为模型特征长度；Ku为土体不排水体积模量，Ku＝K+α2M，其中，K为固体体积模量；α为比奥系数，α的值为1；M为比奥模量，当α＝1时，M＝Kf/n，Kf为流体体积模量，n为土体孔隙率；
流体分析时标的计算公式为：
式中，tf为流体分析时标；Lc f为渗流特征长度；co为广义固结系数，
其中，k为渗透系数。
进一步地，步骤S2中，采用主动筛选法建立蠕变本构模型，具体过程为：整理滑坡滑带土的三轴蠕变试验数据，根据玻尔兹曼线性叠加原理分别得到不同围压下轴向应变-时间关系曲线；根据轴向应变-时间关系曲线的特点，确定蠕变本构模型应该包含的元件类型，如弹性元件、粘性元件等；从元件模型中选择出包含这些元件的模型，再将选择出的模型的典型应变-时间曲线与得到的轴向应变-时间关系曲线进行对照，确定最合理且参数最少的模型为蠕变本构模型。应注意，在选择的蠕变本构模型中必须包含强度参数。
进一步地，蠕变参数的识别过程为：根据蠕变本构模型确定所需要识别的参数，利用蠕变试验数据通过理论公式求解和进行最小二乘法拟合，确定不同围压、不同应力水平条件下蠕变参数的多组数值，并得到由拟合值所描述的曲线与试验曲线的拟合度R2(R2越大，表明曲线拟合效果越好)；然后基于现场监测数据，在数值模型设立对应监测点，选择某一时段作为反演区间，求解相同条件下不同蠕变参数组合滑坡的位移，通过将数值计算结果与监测数据进行比较，与监测数据最接近的数值组合即选择为合理的滑带土的蠕变参数。
进一步地，步骤S3中，滑带土的长期抗剪强度的计算过程为：
S301，计算滑带土长期强度，对于3阶段蠕变曲线，在某级偏应力作用下，若应变-时间曲线上有稳态蠕变阶段向加速蠕变阶段过渡的拐点，则该应力为长期强度；对于2阶段的蠕变曲线，采用等时曲线法确定长期强度，即绘制不同时刻轴向应变-应力等时曲线，确定曲线上的拐点，其对应的偏应力即为长期强度；
S302，在正应力-剪应力平面下，绘制长期强度与其对应围压的摩尔应力圆，绘制各摩尔应力圆的公切线并求得其方程，然后根据摩尔-库伦准则即可计算得出土体的长期抗剪强度。
进一步地，步骤S5中，对长期抗剪强度指标按照公式
和
进行线性增加折减或线性减小折减，其中，F为折减系数；c为初始粘聚力；c′为折减后粘聚力；
为初始内摩擦角；
为折减后内摩擦角；线性增加折减的过程为：从1.1开始，以0.1为间隔线性增加；线性减小折减的过程为：从最后一次折减系数Fmax开始，以0.01为间隔线性减少。
进一步地，步骤S5中，根据GBT32864-2016《滑坡防治工程勘查规范》对滑坡的动态稳定性进行评价，即滑坡稳定系数Fs＜1.00时，滑坡状态为不稳定；1.00≤滑坡稳定系数Fs＜1.05时，滑坡状态为欠稳定；1.05≤滑坡稳定系数Fs＜1.15，滑坡状态为基本稳定；滑坡稳定系数Fs≥1.15，滑坡状态为稳定。
本发明提供的技术方案带来的有益效果是：本发明提供的方法基于流固耦合理论，在FLAC3D软件中施加概化后的滑坡真实水动力条件进行渗流求解，采用滑带土蠕变本构进行力学求解，通过将两种作用相结合，揭示水动力条件下滑坡变形特征与演化机理，实现了对库区蠕变型滑坡稳定性的动态评价，为同类型滑坡的预测预警与防治提供了基础。
附图说明
图1为本发明的一种基于强度折减法的滑坡动态稳定性评价方法的流程示意图。
图2为本发明的滑坡水动力条件概化图。
图3为本发明的滑带土蠕变试验轴向应变-时间曲线(围压＝200kPa)。
图4为Burgers模型典型曲线。
图5为Cvisc元件模型示意图。
图6为蠕变参数反演位移曲线。
图7为滑带土长期抗剪强度。
图8为滑坡数值模型。
图9为滑坡多场特征云图。
图10为滑坡监测位移曲线。
图11为蠕变位移监测点与折减系数曲线。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述。
请参考图1，本发明的实施例提供了一种基于强度折减法的滑坡动态稳定性评价方法，包括以下步骤：
步骤S1，根据滑坡所处的库区水位周期调度情况建立概化模型，利用滑坡岩土体的参数数据计算力学分析时标和流体分析时标，比较力学分析时标和流体分析时标的大小，确定流固耦合模拟方法。
具体地，步骤S1中，根据库水位上升、高水位保持、库水位下降、低水位保持四个阶段对滑坡进行概化，建立水位与时间的函数表达式，记为H＝f(t)，H为水位高程，t为天数；若低水位保持时间较短，可简化为三阶段。
在FLAC3D中，由于完全耦合分析计算量十分庞大，在很多情况下，可以将所模拟问题简化为强耦合、弱耦合和不耦合三种方式中的一种，简化依据主要有时标比值和扰动属性，可以根据表1来进行选择。
表1流固耦合模拟方法
表1中，tm为力学分析时标，其计算公式为：
式中，ρ为密度；G为剪切模量；Lc m为模型特征长度；Ku为土体不排水体积模量，Ku＝K+α2M，其中，K为固体体积模量；α为比奥系数，默认为1；M为比奥模量，当α＝1时，M＝Kf/n，Kf为流体体积模量，n为土体孔隙率。
tf为流体分析时标，其计算公式为：
式中，Lc f为渗流特征长度；co为广义固结系数，
其中，k为渗透系数。
Rk为流固刚度比，其计算公式为：
在渗流计算时，为了保持真实扩散率以及系统特征时标，流体体积模量Kf需要调整为：
通过整理滑坡在以上计算中所涉及的参数，对滑坡数值分析的力学分析时标tm与流体分析时标tf进行计算比较，同时结合扰动属性，选择出最优的流固耦合模拟方法。
步骤S2，建立滑带土的蠕变本构模型，识别滑带土的蠕变参数。
常用的蠕变本构模型有经验模型与元件模型。由于经验模型缺乏理论依据，元件模型建立的本构模型更能真实地反映岩土体本身力学性质，同时具有简单易用等特点，使得后者在实际工程中使用最为广泛。常见的元件模型有Kelvin模型、Poynting-Thomson模型、Bingham模型、Burgers模型、Nishihara模型和Cvisc模型等。
具体地，步骤S2中，采用主动筛选法建立蠕变本构模型，具体过程为：整理滑坡滑带土的三轴蠕变试验数据，根据玻尔兹曼线性叠加原理分别得到不同围压下轴向应变-时间关系曲线；根据轴向应变-时间关系曲线的特点，确定蠕变本构模型应该包含的元件类型，如弹性元件、粘性元件等；从元件模型中选择出包含这些元件的模型，再将选择出的模型的典型应变-时间曲线与得到的轴向应变-时间关系曲线进行对照，确定最合理且参数最少的模型为蠕变本构模型。应注意，在选择的蠕变本构模型中必须包含强度参数。
岩土蠕变参数识别方法分为正分析与反分析，本实施例将两种方法进行了结合，具体地，蠕变参数的识别过程为：根据蠕变本构模型确定所需要识别的参数，利用蠕变试验数据通过理论公式求解和进行最小二乘法拟合，确定不同围压、不同应力水平条件下蠕变参数的多组数值，并得到由拟合值所描述的曲线与试验曲线的拟合度R2(R2越大，表明曲线拟合效果越好)；然后基于现场监测数据，在数值模型设立对应监测点，选择某一时段作为反演区间，求解相同条件下不同蠕变参数组合滑坡的位移，通过将数值计算结果与监测数据进行比较，与监测数据最接近的数值组合即选择为合理的滑带土的蠕变参数。
步骤S3，计算滑带土的长期抗剪强度。
具体地，步骤S3中，滑带土的长期抗剪强度的计算过程为：
S301，计算滑带土长期强度，对于3阶段蠕变曲线，在某级偏应力作用下，若应变-时间曲线上有稳态蠕变阶段向加速蠕变阶段过渡的拐点，则该应力为长期强度；对于2阶段的蠕变曲线，采用等时曲线法确定长期强度，即绘制不同时刻轴向应变-应力等时曲线，确定曲线上的拐点，其对应的偏应力即为长期强度。
S302，在正应力-剪应力平面下，绘制长期强度与其对应围压的摩尔应力圆，绘制各摩尔应力圆的公切线并求得其方程，然后根据摩尔-库伦准则即可计算得出土体的长期抗剪强度。
步骤S4，在FLAC3D中构建滑坡模型；设置位移与渗流边界条件，在滑坡涉水部位通过FISH语言施加步骤S1的概化模型，其余位置设置为不透水边界；利用步骤S2的蠕变参数和步骤S3的长期抗剪强度设置滑坡岩土体本构模型，根据步骤S1的流固耦合模拟方法设置耦合方式，在滑坡模型上设置与现场GNSS监测对应的监测点，设置监测指标，如孔隙水压力、位移等；计算滑坡的应力场、应变场和渗流场，获得滑坡应力场、应变场与渗流场分布规律。
步骤S5，对滑带土长期抗剪强度指标进行线性增加折减，每折减一次，重复计算滑坡的应力场、应变场和渗流场，并记录当前监测点的位移值，绘制位移-折减系数曲线，直至曲线出现突变，停止折减，记录最后一次折减系数Fmax，对最后一次折减系数Fmax进行线性减小折减，每折减一次，重复计算滑坡的应力场、应变场和渗流场，记录当前监测点的位移值，将该当前监测点的位移值插入到位移-折减系数曲线上，直至插值点前未发生位移突变，记录最后一次插值点对应的折减系数Fs，该折减系数Fs即为滑坡当前状态下的稳定性系数。
具体地，步骤S5中，对长期抗剪强度指标按照公式
和
进行线性增加折减或线性减小折减，其中，F为折减系数；c为初始粘聚力；c′为折减后粘聚力；
为初始内摩擦角；
为折减后内摩擦角；线性增加折减的过程为：从1.1开始，以0.1为间隔线性增加；线性减小折减的过程为：从最后一次折减系数Fmax开始，以0.01为间隔线性减少。
根据GBT32864-2016《滑坡防治工程勘查规范》(见表2)对滑坡的动态稳定性进行评价。
表2滑坡稳定状态划分
本发明以呷爬滑坡为例实施评价方法，呷爬滑坡位于锦屏一级水电站坝址上游约11.5km的雅砻江右岸，滑坡纵长约880m，宽约320～400m，面积约0.28km2，总方量约1300万m3，具体过程为：
(1)统计呷爬滑坡所处的雅砻江库区水位周期调度情况，按照库水位上升、1880m高水位保持、库水位下降、1800m保持四个阶段对其进行概化，结果如图2所示，以2017年6月19日～2018年6月14日一个完整水文年为例，将概化模型用函数关系表示为：
式中，H为水位高程；t为天数。
(2)根据地勘资料获得滑坡岩土体的相关参数，见表3(表3中各参数均取国际单位)，将各参数代入力学分析时标和流体分析时标的计算公式，计算得到tm＝0.25，tf＝1.125，tm在tf范围内且未满足tm＜＜tf，同时对于水库滑坡，扰动显然是由于库水位波动所引起的，因此本实施例采用两步耦合的方式来实现流固耦合。
表3呷爬滑坡岩土体相关参数取值
(3)根据针对滑坡滑带土已开展的蠕变试验，收集试验数据，根据玻尔兹曼线性叠加原理分别得到不同围压下轴向应变-时间关系曲线，如图3为围压＝200kPa时的曲线。曲线包含两个明显的特点：在施加荷载后短时间内，土体产生明显的瞬时弹性变形，因此模型应包含弹性元件；随着试验进行，试样产生减速蠕变，一定时间后蠕变速率保持不变，表现为粘弹性特征，因此模型应包含粘性元件。经过筛选，发现Burgers模型的应变-时间曲线(图4)最与其相似，同时也包含了这两种元件，因此将土体的蠕变本构模型初步定为Burgers模型。然而，由于在后续分析中需要使用强度折减法来研究滑坡的稳定性，在FLAC3D中，对不具有强度参数的模型无法进行折减，因此需要在Burgers元件模型的基础上在串联一个M-C元件，以产生塑性变形，进而将模型升级为Cvisc模型(见图5)。
(4)Cvisc模型共包含6个参数，EM、ηM、EK、ηK、c、
分别为弹性模量、马克斯韦尔粘性系数、粘弹性模量、开尔文粘性系数、粘聚力和内摩擦角。后两个参数为抗剪强度指标，常为已知，因此还需确定前四个参数。在土体未进入塑性之前(图3曲线没有出现加速蠕变)，土体表现出弹黏性特性，则Cvisc模型的应力-时间关系可简化为：
式中，σ为应力；ε为应变。
根据试验应变-时间曲线(图3)，可以获得模型的初始应变ε0和稳定蠕变时的曲线斜率m。当t＝0，由上式可得ε0＝σ/EM；当t较大时，
趋于0，则m＝1/ηM；通过最小二乘法可以求出另外两个参数EM、ηM。滑带土最终的蠕变拟合参数如表4。
表4滑带土蠕变参数
从表4中可知共获得了16组滑带土蠕变参数取值，由于室内试验具有尺寸效应等原因，需要进一步通过位移反分析法来确定用于后续稳定性系数计算的参数。本实例选取2017.4.10～2017.6.19低水位维持这一阶段为例，分别将16组参数代入计算，计算过程见步骤(7)-(12)(这里不需要考虑强度衰减)，将16组参数计算得到的相同位置位移监测值与现场真实GNSS监测数据进行比较，最终发现当EM＝19978.6661MPa、ηM＝1.30E+07MPa·h、EK＝41221.29717MPa、ηK＝254029.3004MPa·h时，模拟值与监测值最为接近(图6)。因此选取该组参数作为滑带土的蠕变参数。
(5)计算滑带土长期强度。本案例土体蠕变曲线为2阶段的蠕变曲线，因此采用采用等时曲线法。即绘制不同时刻轴向应变-应力等时曲线，确定曲线上的拐点，其对应的偏应力即为长期强度。最终得出围压为100kpa、200kpa、300kpa、400kpa时的长期强度分别为224.2kpa、293.15kpa、425.94kpa、520.8kpa。
(6)计算滑带土长期抗剪强度。参考图7，在正应力-剪应力平面下，绘制长期强度与其对应围压的摩尔应力圆，绘制各摩尔应力圆的公切线并求得其方程τ＝38.36+1.49σ，其中，τ是剪应力，σ是正应力，然后根据摩尔-库伦准则即可计算得出土体的长期抗剪强度指标(见图7)，其中粘聚力c＝38.36kpa，内摩擦角
(7)在FLAC3D中构建滑坡模型(见图8)。
(8)设置位移与渗流边界条件，其中模型两侧、底部和Y方向采用单向位移约束条件，顶部为自由边界条件；滑坡左侧涉水部分的水力边界施加由FISH语言所编写的步骤(1)的概化模型(本实例采用2017年6月19日～2018年6月14日)，根据地下水信息右侧采用1950m固定水头边界，其余位置为不透水边界。
(9)设置滑坡岩土体本构模型，其中滑带土选用蠕变模型，由步骤(4)得蠕变参数，由步骤(6)得长期抗剪强度指标。由于在计算时渗流时间步与蠕变时间步都特别小，为了提高计算效率本实例仅考虑深层滑带的影响。
(10)设置耦合方式，由步骤(2)知本实例宜采用两步法，因此在求解时确保渗流分析时长与力学分析时长一致。在使用“set fluid off mech on”命令后，将流体模量设置为0，防止因单元体积应变变化引起孔隙水压力的改变。
(11)如图8设置与现场GNSS监测对应的监测点，并记录两点处位移值。
(12)求解，分时间段保存滑坡应力场、应变场与渗流场云图。图9为库水位下降结束后滑坡的多场信息图，图10为滑坡监测曲线，数值解与现场值规律一致，模拟方法效是可靠的，鉴于此，结合滑坡的多场信息可以对滑坡的变形特征和演化机理更深入一步的研究，如：库水位下降阶段滑坡变形速率增大，这表明水位下降造成的高额动水压力致使滑坡变形加剧，应加强在库水位下降阶段滑坡变形的监测工作。
(13)以上述库水位变化三阶段中的水为下降阶段为例，在滑带关键部位设置2个监测点(图8中C1和C2)。对滑带土长期抗剪强度指标进行折减，从1.1开始，以0.1为间隔线性增加。每折减一次，重复计算滑坡应力场、应变场与渗流场，并记录当前监测点位移值，绘制位移-折减系数曲线(图11)，直到曲线出现突变，停止折减，记最后一次折减系数为Fmax。
(14)将Fmax以0.01为间隔线性减少，每折减一次，重复计算滑坡应力场、应变场与渗流场，并记录当前监测点位移值，并将当前位移插入到步骤(13)中所获得的位移-折减系数曲线，直到插值点前未发生位移突变，记录最后一次插值点所对应的折减系数为Fs，则为滑坡当前状态下的稳定系数。折减系数与监测点位移值的关系曲线如图9，由此确定滑坡在当前工况下的稳定系数为1.36。
(15)根据GBT32864-2016《滑坡防治工程勘查规范》，对滑坡的动态稳定性进行评价。由1.36>1.15，表明滑坡在当前工况下处于稳定状态，但是考虑到滑坡的持续变形，应继续加强监测。为了进行比较，本实例还进行了另外两个工况的计算，一个是不考虑蠕变本构，一个是不考虑长期抗剪强度，计算结果如表5所示。
表5不同工况的稳定系数
从表5可以看出，库水位下降对滑坡稳定性影响较大，同时是否考虑滑带土的蠕变特性以及在蠕变本构中是否考虑长期抗剪强度对稳定性的变化也有显著影响。三种工况下滑坡稳定性系数分别下降了7.1％、9.1％和11.7％，这表明对于蠕变型滑坡，若在稳定性评价时不能衡量到蠕变对其的定量影响，将会导致滑坡稳定状态评价偏危险。目前常规评价滑坡稳定性的方法均是静态的，而本发明所提出的方法既考虑考虑到滑坡水动力条件的变化，也考虑到滑坡岩土体自身特性的影响，将滑坡的稳定性系数与时间挂钩，实现了对滑坡稳定状态的动态评价。
在本文中，所涉及的前、后、上、下等方位词是以附图中零部件位于图中以及零部件相互之间的位置来定义的，只是为了表达技术方案的清楚及方便。应当理解，所述方位词的使用不应限制本申请请求保护的范围。
在不冲突的情况下，本文中上述实施例及实施例中的特征可以相互结合。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于强度折减法的滑坡动态稳定性评价方法，其特征在于，包括以下步骤：
S1，根据滑坡所处的库区水位周期调度情况建立概化模型，利用滑坡岩土体的参数数据计算力学分析时标和流体分析时标，比较力学分析时标和流体分析时标的大小，确定流固耦合模拟方法；
S2，建立滑带土的蠕变本构模型，识别滑带土的蠕变参数；
S3，计算滑带土的长期抗剪强度；
S4，构建滑坡模型，在滑坡涉水部位施加步骤S1的概化模型，利用步骤S2的蠕变参数和步骤S3的长期抗剪强度设置滑坡岩土体本构模型，根据步骤S1的流固耦合模拟方法设置耦合方式，在滑坡模型上设置与现场GNSS监测对应的监测点，计算滑坡的应力场、应变场和渗流场；
S5，对滑带土长期抗剪强度指标进行线性增加折减，每折减一次，重复计算滑坡的应力场、应变场和渗流场，并记录当前监测点的位移值，绘制位移-折减系数曲线，直至曲线出现突变，停止折减，记录最后一次折减系数，对最后一次折减系数进行线性减小折减，每折减一次，重复计算滑坡的应力场、应变场和渗流场，记录当前监测点的位移值，将该当前监测点的位移值插入到位移-折减系数曲线上，直至插值点前未发生位移突变，记录最后一次插值点对应的折减系数，该折减系数即为滑坡当前状态下的稳定系数。
2.根据权利要求1所述的基于强度折减法的滑坡动态稳定性评价方法，其特征在于，步骤S1中，根据库水位上升、高水位保持、库水位下降、低水位保持四个阶段对滑坡进行概化，建立水位与时间的函数表达式，记为H＝f(t)，H为水位高程，t为天数。
3.根据权利要求1所述的基于强度折减法的滑坡动态稳定性评价方法，其特征在于，步骤S1中，力学分析时标的计算公式为：
式中，tm为力学分析时标；G为剪切模量；Lc m为模型特征长度；Ku为土体不排水体积模量，Ku＝K+α2M，其中，K为固体体积模量；α为比奥系数，α的值为1；M为比奥模量，当α＝1时，M＝Kf/n，Kf为流体体积模量，n为土体孔隙率；
流体分析时标的计算公式为：
式中，tf为流体分析时标；Lc f为渗流特征长度；co为广义固结系数，
其中，k为渗透系数。
4.根据权利要求1所述的基于强度折减法的滑坡动态稳定性评价方法，其特征在于，步骤S2中，采用主动筛选法建立蠕变本构模型，具体过程为：整理滑坡滑带土的三轴蠕变试验数据，根据玻尔兹曼线性叠加原理分别得到不同围压下轴向应变-时间关系曲线；根据轴向应变-时间关系曲线的特点，确定蠕变本构模型应该包含的元件类型；从元件模型中选择出包含这些元件的模型，再将选择出的模型的典型应变-时间曲线与得到的轴向应变-时间关系曲线进行对照，确定最合理且参数最少的模型为蠕变本构模型。
5.根据权利要求1所述的基于强度折减法的滑坡动态稳定性评价方法，其特征在于，蠕变参数的识别过程为：根据蠕变本构模型确定所需要识别的参数，利用蠕变试验数据确定不同围压、不同应力水平条件下蠕变参数的多组数值；然后基于现场监测数据，在数值模型上设立对应监测点，选择某一时段作为反演区间，求解相同条件下不同蠕变参数组合滑坡的位移，将数值计算结果与监测数据进行比较，与监测数据最接近的数值组合即选择为滑带土的蠕变参数。
6.根据权利要求1所述的基于强度折减法的滑坡动态稳定性评价方法，其特征在于，步骤S3中，滑带土的长期抗剪强度的计算过程为：
S301，计算滑带土长期强度，对于3阶段蠕变曲线，在某级偏应力作用下，若应变-时间曲线上有稳态蠕变阶段向加速蠕变阶段过渡的拐点，则该应力为长期强度；对于2阶段的蠕变曲线，采用等时曲线法确定长期强度，即绘制不同时刻的轴向应变-应力等时曲线，确定曲线上的拐点，其对应的偏应力即为长期强度；
S302，在正应力-剪应力平面下，绘制长期强度与其对应围压的摩尔应力圆，绘制各摩尔应力圆的公切线并求得其方程，然后根据摩尔-库伦准则即可计算得出土体的长期抗剪强度。
7.根据权利要求1所述的基于强度折减法的滑坡动态稳定性评价方法，其特征在于，步骤S5中，对长期抗剪强度指标按照公式
和
进行折减，其中，F为折减系数；c为初始粘聚力；c′为折减后粘聚力；
为初始内摩擦角；
为折减后内摩擦角；线性增加折减的过程为：从1.1开始，以0.1为间隔线性增加；线性减小折减的过程为：从最后一次折减系数开始，以0.01为间隔线性减少。
8.根据权利要求1所述的基于强度折减法的滑坡动态稳定性评价方法，其特征在于，步骤S5中，根据GBT32864-2016《滑坡防治工程勘查规范》对滑坡的动态稳定性进行评价，即滑坡稳定系数Fs＜1.00时，滑坡状态为不稳定；1.00≤滑坡稳定系数Fs＜1.05时，滑坡状态为欠稳定；1.05≤滑坡稳定系数Fs＜1.15，滑坡状态为基本稳定；滑坡稳定系数Fs≥1.15，滑坡状态为稳定。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010513680.2A CN111861107B (zh) | 2020-06-08 | 2020-06-08 | 一种基于强度折减法的滑坡动态稳定性评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010513680.2A CN111861107B (zh) | 2020-06-08 | 2020-06-08 | 一种基于强度折减法的滑坡动态稳定性评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111861107A true CN111861107A (zh) | 2020-10-30 |
| CN111861107B CN111861107B (zh) | 2021-03-30 |
Family
ID=72986341
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010513680.2A Active CN111861107B (zh) | 2020-06-08 | 2020-06-08 | 一种基于强度折减法的滑坡动态稳定性评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111861107B (zh) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112632829A (zh) * | 2020-12-28 | 2021-04-09 | 河海大学 | 基于abaqus-python脚本的滑坡安全系数区间分析方法 |
| CN112733412A (zh) * | 2020-12-28 | 2021-04-30 | 河海大学 | 一种水动力作用滑坡运动机制研究的速度等效表征方法 |
| CN112861409A (zh) * | 2021-02-26 | 2021-05-28 | 山东大学 | 单桩基础承载能力计算方法、系统、存储介质及设备 |
| CN113128060A (zh) * | 2021-04-23 | 2021-07-16 | 四川贡嘎雪新材料有限责任公司 | 一种矿厂老采区潜在滑坡稳定性分析方法 |
| CN113686646A (zh) * | 2021-08-05 | 2021-11-23 | 中国科学院武汉岩土力学研究所 | 一种确定岩石长期强度的方法 |
| CN114018727A (zh) * | 2021-11-05 | 2022-02-08 | 中国地质大学（武汉） | 一种大变形全过程滑带土抗剪强度确定方法 |
| CN114036794A (zh) * | 2021-11-12 | 2022-02-11 | 中国海洋大学 | 一种海底斜坡震后稳定性分析方法 |
| CN114065504A (zh) * | 2021-11-12 | 2022-02-18 | 浙江华东建设工程有限公司 | 一种基于原位监测的海底斜坡震后稳定性分析方法 |
| CN114216831A (zh) * | 2021-12-16 | 2022-03-22 | 长沙学院 | 一种联合中性点计算及有限元模拟的滑坡稳定性辨识方法 |
| CN116244813A (zh) * | 2023-05-11 | 2023-06-09 | 中国铁路设计集团有限公司 | 一种基于触探技术的土体强度参数反演方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105069248A (zh) * | 2015-08-21 | 2015-11-18 | 成都理工大学 | 巨型滑坡多块滑动计算方法 |
| CN108538026A (zh) * | 2018-02-23 | 2018-09-14 | 青岛理工大学 | 一种边坡稳定性的多指标综合评价与预警方法 |
| CN108920754A (zh) * | 2018-05-25 | 2018-11-30 | 重庆交通大学 | 基于应变软化剪切本构模型和位移变化的结构面控制边坡稳定性动态评价方法 |
| US20190113496A1 (en) * | 2018-11-26 | 2019-04-18 | Huiming Tang | Landslide experimental device and experimental method for simulating constant seepage flow |
| CN110210084A (zh) * | 2019-05-18 | 2019-09-06 | 西南交通大学 | 基于时间融合的力学变化判断滑坡可能性的方法及系统 |
| CN110261573A (zh) * | 2019-05-16 | 2019-09-20 | 同济大学 | 一种高位岩质滑坡稳定性动态数值评价方法 |
| JP2020060078A (ja) * | 2018-10-12 | 2020-04-16 | 日本信号株式会社 | 土砂崩落の発生時刻を推定するシステム |
-
2020
- 2020-06-08 CN CN202010513680.2A patent/CN111861107B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105069248A (zh) * | 2015-08-21 | 2015-11-18 | 成都理工大学 | 巨型滑坡多块滑动计算方法 |
| CN108538026A (zh) * | 2018-02-23 | 2018-09-14 | 青岛理工大学 | 一种边坡稳定性的多指标综合评价与预警方法 |
| CN108920754A (zh) * | 2018-05-25 | 2018-11-30 | 重庆交通大学 | 基于应变软化剪切本构模型和位移变化的结构面控制边坡稳定性动态评价方法 |
| JP2020060078A (ja) * | 2018-10-12 | 2020-04-16 | 日本信号株式会社 | 土砂崩落の発生時刻を推定するシステム |
| US20190113496A1 (en) * | 2018-11-26 | 2019-04-18 | Huiming Tang | Landslide experimental device and experimental method for simulating constant seepage flow |
| CN110261573A (zh) * | 2019-05-16 | 2019-09-20 | 同济大学 | 一种高位岩质滑坡稳定性动态数值评价方法 |
| CN110210084A (zh) * | 2019-05-18 | 2019-09-06 | 西南交通大学 | 基于时间融合的力学变化判断滑坡可能性的方法及系统 |
Non-Patent Citations (4)
| Title |
|---|
| LI WANG等: "Construction of 3D Creep Model of Landslide Slip-Surface Soil and Secondary Development Based on FLAC3D", 《ADVANCES IN CIVIL ENGINEERING》 * |
| XINLI HU等: "Model Test and Numerical Analysis on the Deformation and Stability of a Landslide Subjected to Reservoir Filling", 《GEOFLUIDS》 * |
| 孙淼军等: "蠕动型滑坡滑带土蠕变特性研究", 《岩土力学》 * |
| 汤罗圣: "三峡库区堆积层滑坡稳定性与预测预报研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112733412A (zh) * | 2020-12-28 | 2021-04-30 | 河海大学 | 一种水动力作用滑坡运动机制研究的速度等效表征方法 |
| CN112632829A (zh) * | 2020-12-28 | 2021-04-09 | 河海大学 | 基于abaqus-python脚本的滑坡安全系数区间分析方法 |
| CN112733412B (zh) * | 2020-12-28 | 2023-03-24 | 河海大学 | 一种水动力作用滑坡运动机制研究的速度等效表征方法 |
| CN112861409A (zh) * | 2021-02-26 | 2021-05-28 | 山东大学 | 单桩基础承载能力计算方法、系统、存储介质及设备 |
| CN113128060A (zh) * | 2021-04-23 | 2021-07-16 | 四川贡嘎雪新材料有限责任公司 | 一种矿厂老采区潜在滑坡稳定性分析方法 |
| CN113686646A (zh) * | 2021-08-05 | 2021-11-23 | 中国科学院武汉岩土力学研究所 | 一种确定岩石长期强度的方法 |
| CN113686646B (zh) * | 2021-08-05 | 2022-07-12 | 中国科学院武汉岩土力学研究所 | 一种确定岩石长期强度的方法 |
| CN114018727B (zh) * | 2021-11-05 | 2023-09-12 | 中国地质大学（武汉） | 一种大变形全过程滑带土抗剪强度确定方法 |
| CN114018727A (zh) * | 2021-11-05 | 2022-02-08 | 中国地质大学（武汉） | 一种大变形全过程滑带土抗剪强度确定方法 |
| CN114036794A (zh) * | 2021-11-12 | 2022-02-11 | 中国海洋大学 | 一种海底斜坡震后稳定性分析方法 |
| CN114065504A (zh) * | 2021-11-12 | 2022-02-18 | 浙江华东建设工程有限公司 | 一种基于原位监测的海底斜坡震后稳定性分析方法 |
| CN114036794B (zh) * | 2021-11-12 | 2022-07-22 | 中国海洋大学 | 一种海底斜坡震后稳定性分析方法 |
| CN114216831A (zh) * | 2021-12-16 | 2022-03-22 | 长沙学院 | 一种联合中性点计算及有限元模拟的滑坡稳定性辨识方法 |
| CN114216831B (zh) * | 2021-12-16 | 2023-08-01 | 长沙学院 | 一种联合中性点计算及有限元模拟的滑坡稳定性辨识方法 |
| CN116244813B (zh) * | 2023-05-11 | 2023-07-18 | 中国铁路设计集团有限公司 | 一种基于触探技术的土体强度参数反演方法 |
| CN116244813A (zh) * | 2023-05-11 | 2023-06-09 | 中国铁路设计集团有限公司 | 一种基于触探技术的土体强度参数反演方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111861107B (zh) | 2021-03-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111861107B (zh) | 一种基于强度折减法的滑坡动态稳定性评价方法 | |
| Surarak et al. | Stiffness and strength parameters for hardening soil model of soft and stiff Bangkok clays | |
| Preisig et al. | Hydromechanical rock mass fatigue in deep-seated landslides accompanying seasonal variations in pore pressures | |
| Pollen et al. | Estimating the mechanical effects of riparian vegetation on stream bank stability using a fiber bundle model | |
| CN104406623B (zh) | 基于地下水位与位移监测的边坡动态稳定性系数测定方法 | |
| Wang et al. | Landsliding triggered by reservoir operation: a general conceptual model with a case study at Three Gorges Reservoir | |
| Morita et al. | A quick method to determine subsidence, reservoir compaction, and in-situ stress induced by reservoir depletion | |
| CN104252009B (zh) | 一种恢复石油与天然气盖层脆塑性演化史的方法 | |
| McCracken et al. | Origin of the active drumlin field at Múlajökull, Iceland: new insights from till shear and consolidation patterns | |
| Pham et al. | Stability of slope and seepage analysis in earth dam using numerical finite element model | |
| CN107145671A (zh) | 一种油藏数值模拟方法及系统 | |
| Hewitt et al. | The influence of a poroelastic till on rapid subglacial flooding and cavity formation | |
| Wang et al. | Analysis of the abutment movements of high arch dams due to reservoir impoundment | |
| Komasi et al. | Seepage and stability analysis of the Eyvashan Earth Dam under drawdown conditions | |
| Yang et al. | Numerical simulation of the free surface and water inflow of a slope, considering the nonlinear flow properties of gravel layers: a case study | |
| Wright et al. | An examination of slope stability computation procedures for sudden drawdown | |
| Liu et al. | Hysteresis of dam slope safety factor under water level fluctuations based on the LEM coupled with FEM method | |
| Wu et al. | Research on the deformation mechanisms of a tunnel-landslide system based on the point safety factor of the interface | |
| Terbouche et al. | Analysis of pore water pressures in an earth dam under operating conditions (case of Taksebt Dam, Algeria) | |
| Karunawardena | Consolidation analysis of Sri Lankan peaty clay using elasto-viscoplastic theory | |
| Fokker et al. | Estimating the distribution of salt cavern squeeze using subsidence measurements | |
| Wei | Quantifying the Effects of Water Management Decisions on Streambank Stability | |
| Bedayat et al. | Drainage of poroelastic fractures and its implications on the performance of naturally fractured reservoirs | |
| Foroutan et al. | Coupling Effect of Suction Variation on Riverbank Stability | |
| Halleux | Enhancing soft organic clay strength assessment via CPT: a numerical study supported by field data |
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 |