CN112835099A - 一种边坡地震稳定性判定方法 - Google Patents

一种边坡地震稳定性判定方法 Download PDF

Info

Publication number
CN112835099A
CN112835099A CN202110012336.XA CN202110012336A CN112835099A CN 112835099 A CN112835099 A CN 112835099A CN 202110012336 A CN202110012336 A CN 202110012336A CN 112835099 A CN112835099 A CN 112835099A
Authority
CN
China
Prior art keywords
point
sliding body
stress
sliding
slope
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.)
Pending
Application number
CN202110012336.XA
Other languages
English (en)
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.)
Nuclear Industry Southwest Survey & Design Institute Co ltd
Original Assignee
Nuclear Industry Southwest Survey & Design Institute Co ltd
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 Nuclear Industry Southwest Survey & Design Institute Co ltd filed Critical Nuclear Industry Southwest Survey & Design Institute Co ltd
Priority to CN202110012336.XA priority Critical patent/CN112835099A/zh
Publication of CN112835099A publication Critical patent/CN112835099A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种边坡地震稳定性判定方法,其包括以下步骤:建立岩质基覆边坡的概化模型、选取研究点、计算地震波到达研究点的应力σni、计算下滑剪应力τsi、计算发生剪切破坏的总的下滑力F和抗滑力F、对稳定性进行判定。本方案合理的考虑了滑体自重应力场,控制性结构面的内摩擦角,内聚力以及地震波的透射、反射系数,基岩、覆盖层的物理力学参数,地震波纵波的入射和位移时程等影响因素,来对地震波对滑体的稳定性影响进行判定。有效提升了地震对边坡稳定性影响判定的准确性和可靠性,可为地震滑坡、泥石流的预防和治理工作提供有针对性的依据,并可以实现对边坡稳定状态的实时掌握。

Description

一种边坡地震稳定性判定方法
技术领域
本发明涉及地质活动研究技术领域,具体涉及一种边坡地震稳定性判定方法。
背景技术
大量地震灾害显示,地震中死亡的人没有一个是“振”死的。地震造成人员伤亡的根本原因是地震动对工程结构的破坏以及引发的次生灾害,边坡失稳就是其中最为常见的地震次生灾害之一,在山区或丘陵地区,其危害程度有时甚至超越了地震动引起的结构破坏。因此,针对边坡地震稳定性的研究工作历来是岩土工程与工程地震领域重点课题之一。
特殊的地理环境,活跃的地质构造以及飞速发展的城市化建设又使得边坡地震稳定性问题在我国尤为突出。因此,开展边坡地震稳定性研究,对于我国边坡灾害的风险评估、抗震加固设计以及减轻滑坡灾害造成的损失具有十分重要的理论意义和工程应用价值。
发明内容
针对现有技术的上述不足,本发明提供了一种能准确判定边坡地震稳定性的判定方法。
为达到上述发明目的,本发明所采用的技术方案为:
S1:建立岩质基覆边坡的概化模型;
S2:将概化模型中的滑体结构面均匀网格化;
S3:选取网格化后的滑体结构面上的某一网格点i作为研究点;
S4:利用地震波到达点i的入射纵波、反射纵波和反射横波计算地震波到达点i时在滑体结构面上的应力σni
S5:利用应力σni计算点i上的下滑剪应力τsi
S6:重复步骤S3-S5,计算滑体结构面上其余点的应力和下滑剪应力;
S7:计算滑体结构面发生剪切破坏的总的下滑力F和抗滑力F
Figure BDA0002885475460000021
其中,n为在滑体结构面上选取的研究点个数,i为选取的研究点,τsi为第i个研究点的下滑剪应力,Ai=0.5(Li-1+Li+1)×1,表示第i个研究点所代表的面积,其中Li-1和Li+1分别为第i-l个网格点和i+1个网格点的边长;σ0i表示点i处的法向自重力;
Figure BDA0002885475460000022
和Ci分别表示滑体结构面上i点处的内摩擦角和粘聚力。
S8:若K=F/F≤K0,则滑体危险;若K=F/F>K0,则滑体安全,其中K0为边坡的安全性系数;
S9:计算滑体结构面上发生拉伸破坏的总的拉应力F和抗拉力F′
Figure BDA0002885475460000023
其中,σni为i个研究点的拉应力,σtensioni表示第i个研究点允许拉应力;
S10:若K=F/F≤K0,则滑体危险;若K=F/F>K0,则滑体安全。
本发明的有益效果为:本方案合理的考虑了滑体自重应力场,控制性结构面的内摩擦角,内聚力以及地震波的透射、反射系数,基岩、覆盖层的物理力学参数,地震波纵波的入射和位移时程等影响因素,来对地震波对滑体的稳定性影响进行判定。有效提升了地震对边坡稳定性影响判定的准确性和可靠性,可为地震滑坡、泥石流的预防和治理工作提供有针对性的依据,并可以实现对边坡稳定状态的实时掌握。
附图说明
图1为边坡地震稳定性判定方法的流程图。
图2为岩质基覆边坡的概化模型图。
图3为地震波的反射和透射分析图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,本方案的边坡地震稳定性判定方法包括以下步骤:
S1:建立岩质基覆边坡的概化模型;
大量的地震震害调查结果显示,地震灾区内最为普遍、失稳数量最多的斜坡坡体结构类型为基覆型高陡边坡,而该类型的边坡均存在一个共同点:坡体内部均存在明显的结构面,其两侧材料的性质存在一定差异,并且结构面是影响边坡滑塌的主控因素。基于上述思想,本发明将概化基覆型高陡边坡的几何计算模型,如图2所示。
S2:将概化模型中的滑体结构面均匀网格化;
S3:选取网格化后的滑体结构面上的某一网格点i作为研究点;
地震波传至滑体结构面时,地震波在滑体结构面处将产生大量的透射和反射,并发生波形转换,最终形成了各种类型的地震波如纵波、横波、瑞雷波、洛夫波等,其中以纵波和横波对高陡边坡造成的破坏性最大。因滑体结构面处的作用机理,如图3所示。
S4:利用地震波到达点i的入射纵波、反射纵波和反射横波计算地震波到达点i时在滑体结构面上的应力σni
Figure BDA0002885475460000041
其中,λ1为拉梅常数,
Figure BDA0002885475460000042
为入射纵波,
Figure BDA0002885475460000043
为反射纵波,
Figure BDA0002885475460000044
为反射横波;
Figure BDA0002885475460000045
Figure BDA0002885475460000046
均为波矢,且
Figure BDA0002885475460000047
α为纵波的反射角,α′为横波的反射角;μ1为滑体的水平位移,且
Figure BDA0002885475460000048
Figure BDA0002885475460000049
均为波矢,且
Figure BDA00028854754600000410
σ0为初始自重应力场在点i的法向应力。
S5:利用应力σni计算点i上的下滑剪应力τsi
Figure BDA00028854754600000411
S6:重复步骤S3-S5,计算滑体结构面上其余点的应力和下滑剪应力;
S7:计算滑体结构面发生剪切破坏的总的下滑力F和抗滑力F
Figure BDA00028854754600000412
其中,n为在滑体结构面上选取的研究点个数,i为选取的研究点,τsi为第i个研究点的下滑剪应力,Ai=0.5(Li-1+Li+1)×1,表示第i个研究点所代表的面积,其中Li-1和Li+1分别为第i-l个网格点和i+1个网格点的边长;σ0i表示点i处的法向自重力;
Figure BDA00028854754600000413
和Ci分别表示滑体结构面上i点处的内摩擦角和粘聚力。
S8:若K=F/F≤K0,则滑体危险;若K=F/F>K0,则滑体安全,其中K0为边坡的安全性系数,利用设计所需的安全性系数K0进行边坡安全性的判断;
S9:计算滑体结构面上发生拉伸破坏的总的拉应力F和抗拉力F′
Figure BDA00028854754600000414
其中,σni为i个研究点的拉应力,σtensioni表示第i个研究点允许拉应力;
S10:若K=F/F≤K0,则滑体危险;若K=F/F>K0,则滑体安全。
本方案合理的考虑了滑体自重应力场,控制性结构面的内摩擦角,内聚力以及地震波的透射、反射系数,基岩、覆盖层的物理力学参数,地震波纵波的入射角和位移时程等影响因素,来对地震波对滑体的稳定性影响进行判定。有效提升了地震对边坡稳定性影响判定的准确性和可靠性,可为地震滑坡、泥石流的预防和治理工作提供有针对性的依据,并可以实现对边坡稳定状态的实时掌握。

Claims (4)

1.一种边坡地震稳定性判定方法,其特征在于,包括以下步骤:
S1:建立岩质基覆边坡的概化模型;
S2:将概化模型中的滑体结构面均匀网格化;
S3:选取网格化后的滑体结构面上的某一网格点i作为研究点;
S4:利用地震波到达点i的入射纵波、反射纵波和反射横波计算地震波到达点i时在滑体结构面上的应力σni
S5:利用应力σni计算点i上的下滑剪应力τsi
S6:重复步骤S3-S5,计算滑体结构面上其余点的应力和下滑剪应力;
S7:计算滑体结构面发生剪切破坏的总的下滑力F和抗滑力F
Figure FDA0002885475450000011
其中,n为在滑体结构面上选取的研究点个数,i为选取的研究点,τsi为第i个研究点的下滑剪应力,Ai=0.5(Li-1+Li+1)×1,表示第i个研究点所代表的面积,其中Li-1和Li+1分别为第i-l个网格点和i+1个网格点的边长;σ0i表示点i处的法向自重力;
Figure FDA0002885475450000012
和Ci分别表示滑体结构面上i点处的内摩擦角和粘聚力。
S8:若K=F/F≤K0,则滑体危险;若K=F/F>K0,则滑体安全,其中K0为边坡的安全性系数;
S9:计算滑体结构面上发生拉伸破坏的总的拉应力F和抗拉力F′
Figure FDA0002885475450000013
其中,σni为i个研究点的拉应力,σtensioni表示第i个研究点允许拉应力;
S10:若K=F/F≤K0,则滑体危险;若K=F/F>K0,则滑体安全。
2.根据权利要求1所述的边坡地震稳定性判定方法,其特征在于,所述应力σni为:
Figure FDA0002885475450000021
其中,λ1为拉梅常数,
Figure FDA0002885475450000022
为入射纵波,
Figure FDA0002885475450000023
为反射纵波,
Figure FDA0002885475450000024
为反射横波,
Figure FDA0002885475450000025
Figure FDA0002885475450000026
均为波矢,且
Figure FDA0002885475450000027
α为纵波的反射角,α′为偏振横波的反射角,μ1为滑体的水平位移,
Figure FDA0002885475450000028
均为波矢,且
Figure FDA0002885475450000029
σ0为初始地应力场在研究点的法向应力。
3.根据权利要求2所述的边坡地震稳定性判定方法,其特征在于,所述水平位移μ1为:
Figure FDA00028854754500000210
4.根据权利要求1所述的边坡地震稳定性判定方法,其特征在于,所述下滑剪应力τsi为:
Figure FDA00028854754500000211
CN202110012336.XA 2021-01-06 2021-01-06 一种边坡地震稳定性判定方法 Pending CN112835099A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110012336.XA CN112835099A (zh) 2021-01-06 2021-01-06 一种边坡地震稳定性判定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110012336.XA CN112835099A (zh) 2021-01-06 2021-01-06 一种边坡地震稳定性判定方法

Publications (1)

Publication Number Publication Date
CN112835099A true CN112835099A (zh) 2021-05-25

Family

ID=75926095

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110012336.XA Pending CN112835099A (zh) 2021-01-06 2021-01-06 一种边坡地震稳定性判定方法

Country Status (1)

Country Link
CN (1) CN112835099A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115754225A (zh) * 2022-10-21 2023-03-07 西南交通大学 基于加速度判断上下坡体相互作用的地震稳定性分析方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108873069A (zh) * 2018-05-08 2018-11-23 四川大学 一种边坡地震稳定性的时频分析方法
US20190250291A1 (en) * 2018-02-09 2019-08-15 China University Of Geosciences, Beijing Method and system for acquiring probability of slope failure and destabilization caused by earthquake

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190250291A1 (en) * 2018-02-09 2019-08-15 China University Of Geosciences, Beijing Method and system for acquiring probability of slope failure and destabilization caused by earthquake
CN108873069A (zh) * 2018-05-08 2018-11-23 四川大学 一种边坡地震稳定性的时频分析方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
杨长卫: "岩质边坡的地震动特性及基覆型边坡的滑坡成因机理、稳定性判识、危害范围评价体系旳研究", 中国博士学位论文全文数据库基础科学辑, no. 1, pages 99 - 106 *
杨长卫;张建经;周德培;: "SV波作用下岩质边坡地震稳定性的时频分析方法研究", 岩石力学与工程学报, vol. 32, no. 03, pages 483 - 491 *
耿鹏超: "地震作用下边坡背坡面效应的机理研究", 中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑, no. 2, pages 10 - 16 *
范刚: "含软弱夹层层状岩质边坡地震响应及稳定性判识时频方法研究", 中国博士学位论文全文数据库 工程科技Ⅱ辑, no. 4 *
黄钊: "地震波作用下缓坡失稳破坏的机理研究", 中国优秀硕士学位论文全文数据库 基础科学辑, no. 8, pages 10 - 39 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115754225A (zh) * 2022-10-21 2023-03-07 西南交通大学 基于加速度判断上下坡体相互作用的地震稳定性分析方法
CN115754225B (zh) * 2022-10-21 2024-03-15 西南交通大学 基于加速度判断上下坡体相互作用的地震稳定性分析方法

Similar Documents

Publication Publication Date Title
Ghanaat Failure modes approach to safety evaluation of dams
Tocher Earthquake energy and ground breakage
Seed et al. Dynamic analysis of the slide in the Lower San Fernando Dam during the earthquake of February 9, 1971
Zheng et al. Theoretical and numerical study on the block-flexure toppling failure of rock slopes
Ponter et al. The force exerted by a moving ice sheet on an offshore structure: Part 1. The creep mode
CN112835099A (zh) 一种边坡地震稳定性判定方法
Zou et al. Dislocation of face-slabs of Zipingpu concrete face rockfill dam during Wenchuan earthquake
CN109214054B (zh) 基于位移的边坡动力稳定性分析方法
Krawinkler et al. Effects of soft soils on design spectra
Du et al. Quantitative vulnerability estimation of structures for individual landslide: application to the Metropolitan Area of San Salvador, El Salvador
Wakai et al. Numerical modeling of an earthquake-induced landslide considering the strain-softening characteristics at the bedding plane
Danay et al. Seismic-induced slip of concrete gravity dams
Zhao et al. Pseudo-static stability analysis of wedges based on the nonlinear Barton-Bandis failure criterion
Huang et al. One simplified method for seismic stability analysis of an unsaturated slope considering seismic amplification effect
Azizan et al. Incremental dynamic analysis of koyna dam under repeated ground motions
Győri et al. Methods and uncertainties in liquefaction hazard assessment for nuclear power plants
Kenney A review of recent research on strength and consolidation of soft sensitive clays
Chieffo et al. Seismic vulnerability assessment of a Romanian historical building under near-source earthquake
Partono et al. Development of seismic microzonation maps of Semarang, Indonesia
Moftakhar et al. Comparison of stability criteria for concrete dams in different approximate methods based on finite element analysis
Liao et al. Buckling Analysis of Bedded Rock Slope Subject to Seismic Force–a Case Study in Karakoram Highway, Pakistan
Widyastuti et al. Earthquake vulnerability of port structures in Indonesia
Farrokhzad Depth reduction factor assessment for evaluation of cyclic stress ratio based on site response analysis
Mandal et al. Slope stability analysis by static and dynamic method
Frydman et al. Analysis of seismically triggered slides off Israel

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