CN111912953A - 一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法 - Google Patents

一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法 Download PDF

Info

Publication number
CN111912953A
CN111912953A CN202010756699.XA CN202010756699A CN111912953A CN 111912953 A CN111912953 A CN 111912953A CN 202010756699 A CN202010756699 A CN 202010756699A CN 111912953 A CN111912953 A CN 111912953A
Authority
CN
China
Prior art keywords
excavation
slope
mine
monitoring
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
Application number
CN202010756699.XA
Other languages
English (en)
Other versions
CN111912953B (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.)
Qindao University Of Technology
Original Assignee
Qindao University Of Technology
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 Qindao University Of Technology filed Critical Qindao University Of Technology
Priority to CN202010756699.XA priority Critical patent/CN111912953B/zh
Publication of CN111912953A publication Critical patent/CN111912953A/zh
Application granted granted Critical
Publication of CN111912953B publication Critical patent/CN111912953B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/18Status alarms
    • G08B21/182Level alarms, e.g. alarms responsive to variables exceeding a threshold
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/23Dune restoration or creation; Cliff stabilisation

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
  • Medicinal Chemistry (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Emergency Management (AREA)
  • Food Science & Technology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Business, Economics & Management (AREA)

Abstract

本发明公开一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法,首先确定深凹矿山边坡特征参数,并选取开挖位移监测点;实时监测边坡位移,得到边坡开挖方量;确定边坡开挖方量与边坡滑移动力增载量之间的关系;确定矿山边坡开挖动力增载位移响应率;确定矿山边坡开挖动力增载位移响应比;确定基于矿山边坡开挖动力增载位移响应比的失效概率评价指标;根据得到的失效概率实现对矿山边坡稳定性的测定。本方案根据边坡的失稳动力灾变机理,以矿山边坡每次开挖挖方量作为动力增载量,以开挖引起的位移作为对应的响应量,将边坡失稳的动因、机理与位移响应变化进行有机耦合,有效解决传统极限平衡法和位移时序法应用的不足与局限。

Description

一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法
技术领域
本发明属于边坡稳定性评价与监测预警技术领域,具体涉及一种基于开挖量和位移实时监测的深凹矿山开采边坡稳定性测定及评价方法。
背景技术
自20世纪90年代,随着矿产资源的开发和露天采矿业的不断发展,许多矿山已转为深凹露天矿,边坡不断加高加陡,滑坡事故频繁发生,危害日趋严重。因此,如何对该类滑坡的卸荷失稳过程进行有效监测预警并将其造成的损失降低到最低程度,是我国矿山防灾减灾目标的重要组成部分,也是灾害预警与防治领域的当务之急。
在传统矿山边坡稳定性预测方法中,应用最广泛的是极限平衡法和位移时序法。极限平衡法是一种以边坡失稳机理与受力条件清晰明确为基础的力学评价模型,并且具有明确的失稳判据,即稳定性系数Fs=1。但该评价模型是静态的,不能将时间因素考虑在内,并且对于矿山边坡而言,持续的开挖扰动贯穿整个边坡开挖过程,坡体内始终在进行应力重分布,开挖过程中边坡高度和形状不断发生变化,这使得无法考虑土体应力应变关系的极限平衡法在矿山边坡的稳定性评价中存在一定的局限性;边坡位移时序预测法作为边坡稳定性评价方法之一,只能通过单维度的位移量或位移速率的变化预测边坡稳定性,无法从多方面、多维度考虑引起边坡位移发生变化的因素,即不能究其本源对边坡位移变化的力学动因作出解释,所以位移时序预测法的准确性会受到很大影响,无法精准判断出边坡的复杂动力变化规律。
因此,寻求和建立一种能同时克服传统极限平衡法与位移时序预测法二者弊端的方法,对矿山开采边坡稳定性评价与监测预警及防治具有重要的工程应用价值。
发明内容
本发明针对现有矿山开挖边坡稳定性监测方法与预测模型存在的不足,提出一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法,以深凹矿山边坡开挖量作为动力增载量,以开挖引起的位移或位移速率作为其动力响应量,将边坡开挖量与位移或位移速率极限耦合集成,依此确定开采边坡动力增载位移响应比评价参数的稳定性预测判据,实现对边坡稳定性的测定。
本发明是采用以下的技术方案实现的:一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法,包括以下步骤:
步骤一、确定深凹矿山边坡特征参数,并选取开挖位移监测点,所述矿山边坡特征参数包括边坡开挖台阶数量为m,各开挖步高度hi,各开挖坡角θi
步骤二、实时监测边坡位移,得到边坡开挖方量;
步骤三、根据步骤二中得到的边坡开挖方量,确定边坡开挖方量与边坡滑移动力增载量之间的关系;
步骤四、根据步骤二实时监测的边坡位移及步骤三得到的滑移动力增载量确定矿山边坡开挖动力增载位移响应率;
步骤五、根据步骤四得到的结果确定矿山边坡开挖动力增载位移响应比,由此建立矿山边坡动力增载位移响应比评价模型;
步骤六、将矿山边坡开挖动力增载位移响应比作为评价参数应用到边坡稳定性判定中,确定基于矿山边坡开挖动力增载位移响应比的失效概率评价指标;
步骤七、进而根据步骤六得到的失效概率实现对矿山边坡稳定性的测定。
进一步的,根据权利要求1所述的基于挖方量监测的深凹矿山开采边坡稳定性测定方法,其特征在于:所述步骤二具体包括以下步骤:
(1)分析获得矿山边坡开挖过程中各监测点k布设后第i步开挖完成时的位移变化量
Figure BDA0002611819560000021
ΔSi k=Si k-Si-1 k (1)
式中:
Figure BDA0002611819560000022
为监测点k布设后第i-1步开挖完成时的合位移,k=1,2,…,m,
Figure BDA0002611819560000023
监测点k布设后第i步开挖完成时的合位移;
(2)得到监测点k布设后第i开挖步开挖方量:
Figure BDA0002611819560000024
式中:
Figure BDA0002611819560000025
为监测点k布设后第i开挖步上、下断面面积;
Figure BDA0002611819560000026
作为各开挖步的开挖影响范围,
Figure BDA0002611819560000027
θi为第i开挖步的坡角;d为坡体厚度,Li为开挖高度,Li=hi
进一步的,所述步骤三具体包括以下步骤:
(1)确定矿山边坡开挖方量与开挖坡角的关系:
Figure BDA0002611819560000028
(2)确定矿山边坡开挖方量与开挖坡角均值的关系:
Figure BDA0002611819560000031
式中:
Figure BDA0002611819560000032
为监测点k布设后第1-i步开挖完成时,整体开挖方量的上、下断面面积;Li(1-i)为1-i步整体开挖高度;
Figure BDA0002611819560000033
为监测点k布设后第1-i步整体开挖方量的影响范围;a为开挖台阶宽度,
Figure BDA0002611819560000034
为开挖坡角的均值,
Figure BDA0002611819560000035
(3)确定边坡滑移动力增载量与开挖坡角均值的关系,进而得到边坡开挖方量与边坡滑移动力增载量之间的关系;
其中,将监测点k布设后第1-i开挖步的稳定性系数降低变化量ΔFsi k与开挖体积Vi k的乘积作为该开挖步的滑移动力增载量ΔTi k,即:
Figure BDA0002611819560000036
得到开挖坡角θi、开挖坡角均值
Figure BDA0002611819560000037
与滑移动力增载量成正相关关系。
进一步的,所述步骤四中,边坡开挖动力增载位移响应率采用以下方式确定:
将边坡监测点k布设后第i步开挖完成时的下滑动力增量ΔTi k作为该开挖步的滑移动力增载量,边坡开挖引起的对应位移
Figure BDA0002611819560000038
作为该开挖步的位移响应量,进而得到边坡监测点k 布设后第i步开挖完成时的动力增载位移响应率
Figure BDA0002611819560000039
Figure BDA00026118195600000310
进一步的,所述步骤一中,开挖位移监测点的位置采用以下选取方式:
(1)矿山边坡开挖前,在第一个开挖台阶临空面拐角处设置一个初始位移监测点;
(2)此后,在每开挖到一个台阶时则在新的台阶开挖临空面拐角处布设一个位移监测点,并设布设的监测点为k(k=1,2,…,m)。
进一步的,所述步骤五中,将矿山边坡监测点k布设后第i步开挖完成时的动力增载位移响应率
Figure BDA00026118195600000311
与初始动力增载位移响应率
Figure BDA00026118195600000312
之比定义为监测点k的动力增载位移响应比,由此建立的矿山边坡动力增载位移响应比评价模型为:
Figure BDA00026118195600000313
进一步的,所述步骤六在确定基于矿山边坡开挖动力增载位移响应比的失效概率评价指标时,具体采用以下方式:
(1)对矿山边坡监测点k在布设后第i步开挖完成时,分别计算其动力增载位移响应比的均值
Figure BDA0002611819560000041
与标准差
Figure BDA0002611819560000042
Figure BDA0002611819560000043
Figure BDA0002611819560000044
(2)获得矿山边坡监测点k的动力增载位移响应比可靠度指标
Figure BDA0002611819560000045
Figure BDA0002611819560000046
(3)得到矿山边坡监测点k布设后第i步开挖完成时发生破坏的概率,即失效概率
Figure BDA0002611819560000047
Figure BDA0002611819560000048
式中,
Figure BDA0002611819560000049
为标准正态分布函数。
与现有技术相比,本发明的优点和积极效果在于:
本方案根据边坡的失稳动力灾变机理,以矿山边坡每次开挖挖方量作为动力增载量,以开挖引起的位移作为对应的响应量,确定边坡动力增载位移响应比参数和评价模型,并将其作为矿山边坡开挖过程中边坡稳定性预测预报的评价参数;同时,根据数理统计基本原理,建立了动力增载位移响应比与边坡失效概率之间的关系,并以此关系为依据确定出基于边坡动力增载位移响应比评价参数的稳定性预测判据,该方法可将边坡失稳的动因、机理与位移响应变化进行有机耦合,有效解决传统极限平衡法和位移时序法应用的不足与局限。
附图说明
图1为本发明实施例所述边坡稳定性测定方法的流程图;
图2为本发明实施例所述位移监测仪器布设示意图;
图3为本发明实施例边坡开挖示意图。
具体实施方式
为了能够更加清楚地理解本发明的上述目的、特征和优点,下面结合附图及实施例对本发明做进一步说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用不同于在此描述的其他方式来实施,因此,本发明并不限于下面公开的具体实施例。
本方案提出一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法,包括以下步骤:
步骤一、确定深凹矿山边坡特征参数,并选取开挖位移监测点,所述矿山边坡特征参数包括边坡开挖台阶数量为m,各开挖步高度hi,各开挖坡角θi
步骤二、实时监测边坡位移,得到边坡开挖方量;
步骤三、根据步骤二中得到的边坡开挖方量,确定边坡开挖方量与边坡滑移动力增载量之间的关系;
步骤四、根据步骤二实时监测的边坡位移及步骤三得到的滑移动力增载量确定矿山边坡开挖动力增载位移响应率;
步骤五、根据步骤四得到的结果确定矿山边坡开挖动力增载位移响应比,由此建立矿山边坡动力增载位移响应比评价模型;
步骤六、将矿山边坡开挖动力增载位移响应比作为评价参数应用到边坡稳定性判定中,确定基于矿山边坡开挖动力增载位移响应比的失效概率评价指标;
步骤七、进而根据步骤六得到的失效概率实现对矿山边坡稳定性的测定。
为了更好地阐述本发明,下面结合甘肃省内的某石英石矿床,对该矿山边坡开挖的稳定性进行分析评价,确定其失效概率,具体的实施步骤如下:
为实现上述目的,根据莫尔-库伦强度理论(原理1)、弹塑性力学、岩土力学原理(原理2)和统计学原理(原理3),对本发明进行详细说明,其具体步骤如下:
步骤一:深凹矿山边坡特征参数的确定与开挖位移监测点选取
1)深凹矿山边坡特征参数的确定
在深凹矿山边坡动态开挖过程中,随着边坡开挖深度、坡面角度不断增加,边坡将表现出坡体逐渐增高加陡、坡形自上而下朝着背离临空面方向发展且呈凹陷形态的特征。对待评价的深凹矿山边坡进行初步勘探和测绘,并根据深凹矿山边坡的动态开挖特点,分析边坡的范围和尺寸,确定边坡开挖台阶数量为m,各开挖步高度hi,各开挖步坡面角度θi作为深凹矿山边坡特征参数(i取1,2,…m,即坡面角的数量与开挖台阶数量相同),其中,θi>θi-1;本实施例中,开挖台阶数量为7,各开挖步高度12,各开挖步坡面角度θi分别为:68°、70°、 72°、74°、76°、78°、80°。
2)矿山边坡开挖位移监测点选取
矿山边坡开始开挖前,在第一个台阶开挖临空面拐角处设置一个初始位移监测点,此后在开挖过程中每开挖到一个台阶时则在新的台阶开挖临空面拐角处布设一个位移监测点,监测点记为k(k=1,2,…,m),位移监测仪器布设如图2所示。本实施例中,位移监测设备采用高精度测斜仪,通过使用专业的超前钻孔设备并采用搭接法进行仪器埋设孔位的钻孔工作,监测各台阶监测点的合位移。在边坡台阶上打下钻孔后,将测斜仪安置在测斜管内,之后将测斜管放入测斜孔,各测斜孔通过数据线相连,并接入现场的数据采集箱,数据采集箱通过自身的无线发射装置,将实时采集的数据传入服务器,避免人工采集的弊端,达到安全且高效的监测效果。
步骤二:边坡位移实时监测与边坡开挖挖方量的确定
1)矿山开采边坡位移实时监测的确定
根据施工方案,对矿山边坡每次开挖引起的位移响应进行实时监测,同时通过监测设备的无线电采集装置将各监测点k(k=1,2,…,m)的位移监测数据传输到主机,然后对监测数据进行预处理,并将各监测点k(k=1,2,…,m)在布设后第i步开挖完成时的合位移
Figure BDA0002611819560000061
录入Excel 表格。
根据公式(1),计算矿山边坡开挖过程中各监测点k(k=1,2,…,m)布设后第i步开挖完成时的位移变化量
Figure BDA0002611819560000062
计算结果见表1:
ΔSi k=Si k-Si-1 k (1)
式中:
Figure BDA0002611819560000063
为监测点k布设后第i-1步开挖完成时的合位移;
Figure BDA0002611819560000064
监测点k布设后第i步开挖完成时的合位移。
表1监测点k在各开挖步下的位移变化量(单位:m)
Figure BDA0002611819560000065
2)矿山边坡实时开挖挖方量的确定
根据《矿山地质规范》中矿体体积的计算规定,同时考虑到矿山中部位置距离边坡较远,开挖量对边坡影响较小,只将开挖至边坡附近且受开挖深度直接影响的挖方量参与计算,并将各开挖步前3m开挖方量作为初始开挖体积。监测点k布设后第1开挖步开挖体积计算如下(开挖示意图见图3),计算结果见表2:
Figure BDA0002611819560000071
式中:
Figure BDA0002611819560000072
为监测点k布设后第1开挖步上、下断面面积,
Figure BDA0002611819560000073
可作为各开挖步的开挖影响范围,
Figure BDA0002611819560000074
L1为开挖高度,L1=h1;a为开挖台阶宽度,查阅《滑坡防治工程设计与施工技术规范》(DZ/T0219-2006)进行取值,各开挖步的取值相同;d为坡体厚度,θ1为第1开挖步的坡角。
监测点k布设后第i开挖步开挖体积计算如下,计算结果见表2:
Figure BDA0002611819560000075
式中:
Figure BDA0002611819560000076
为监测点k布设后第i开挖步上、下断面面积;θi为第i开挖步的坡角;Li为开挖高度,Li=hi;其他符号意义同式(2)。
表2各监测点布设后第i开挖步体积(单位:m3)
Figure BDA0002611819560000077
步骤三:边坡开挖方量与边坡滑移动力增载量关系的确定
1)深凹矿山边坡开挖方量与开挖坡角关系
根据步骤二中深凹矿山边坡实时开挖方量的计算可知,开挖方量与开挖坡角的关系如下:
Figure BDA0002611819560000081
由式(4)可知,各开挖步的开挖高度hi、坡体厚度d、开挖影响范围
Figure BDA0002611819560000082
为相对不变的物理几何参数,开挖坡角θi为动态变化参数,开挖坡角θi越大,则开挖体积Vi k也越大。
2)深凹矿山边坡开挖方量与开挖坡角均值的关系
为研究深凹矿山边坡整体开挖方量与开挖坡角均值的关系,在明确了深凹矿山边坡各开挖步坡角是逐渐增大的同时,则坡角均值也在逐渐增加。将各开挖步前2m开挖方量作为初始开挖体积;监测点k布设后第1步开挖完成时的整体开挖方量即该步开挖体积本身,见式 (2);监测点k布设后第2步开挖完成时,整体开挖方量与坡角均值的关系如下,计算结果见表3:
Figure BDA0002611819560000083
式中:
Figure BDA0002611819560000084
为监测点k布设后第1-2步开挖完成时,整体开挖方量的上、下断面面积; Li(1-2)为1-2步整体开挖高度;
Figure BDA0002611819560000085
为开挖坡角的均值,
Figure BDA0002611819560000086
监测点k布设后第i步开挖完成时,整体开挖方量与坡角均值的关系如下,计算结果见表3:
Figure BDA0002611819560000087
式中:
Figure BDA0002611819560000088
为监测点k布设后第1-i步开挖完成时,整体开挖方量的上、下断面面积; Li(1-i)为1-i步整体开挖高度;
Figure BDA0002611819560000089
为监测点k布设后第1-i步整体开挖方量的影响范围;
Figure BDA00026118195600000810
为开挖坡角的均值,
Figure BDA00026118195600000811
由式(6)可知,随着坡角θi的增大,坡角均值
Figure BDA00026118195600000812
也在逐渐增大,式中hi、d、
Figure BDA00026118195600000813
均为相对不变的物理几何参数,
Figure BDA0002611819560000091
为动态变化参数,开挖坡角均值
Figure BDA0002611819560000092
越大,则整体开挖体积
Figure BDA0002611819560000093
越大。
表3各监测点布设后整体第i开挖步体积(单位:m3)
Figure BDA0002611819560000094
3)边坡滑移动力增载量与开挖坡角均值的关系
根据步骤二,对于无明显结构面的完整岩体边坡,监测点k布设后第i步开挖完成时边坡的整体稳定性系数Fsi k如下,计算结果见表4:
Figure BDA0002611819560000095
式中:γ为岩(土)体重度(KN/m3);
Figure BDA0002611819560000096
为监测点k布设后第1-i步开挖完成时的整体体积;αf为破裂角,取
Figure BDA0002611819560000097
A为潜在结构面面积(m2),坡体厚度取单位1时, A=1*l’=l’,即结构面长度,对于无明显结构面完整岩体边坡,l’可根据破裂角及《建筑边坡工程技术规范》(GB 50330-2013)中边坡勘探点距离计算取值;
Figure BDA0002611819560000098
为岩(土)体的内摩擦角; c为岩(土)体的粘聚力。
随着矿山边坡开挖坡角均值越来越大,边坡下滑动力增加,进而使边坡的稳定性系数Fsi降低。故将监测点k布设后第1-i开挖步的稳定性系数降低变化量ΔFsi k(见式8)与开挖体积 Vi k的乘积作为该开挖步的滑移动力增载量ΔTi k,根据式(9)计算如下,计算结果见表5:
Figure BDA0002611819560000099
Figure 3
根据式(9)可知,在其他几何、物理参数相对不变的情况下,因开挖坡角逐渐增大,坡角均值也逐渐增大,且
Figure BDA0002611819560000102
故随着开挖坡角θi、开挖坡角均值
Figure BDA0002611819560000103
的增大,深凹矿山边坡的滑移动力增载量也逐渐增大,开挖坡角θi、开挖坡角均值
Figure BDA0002611819560000104
与滑移动力增载量成正相关关系。
表4监测点k布设后第i步开挖完成时整体稳定性系数
Figure BDA0002611819560000105
表5监测点k布设后第i开挖步下滑动力增量变化量
Figure BDA0002611819560000106
步骤四:矿山边坡开挖动力增载位移响应率的确定
以弹塑性力学、岩土力学原理中的边坡受力变形规律为依据,将矿山边坡监测点k布设后第1步开挖完成时的下滑动力增量
Figure BDA0002611819560000107
作为初始滑移动力增载量,初始开挖引起的位移
Figure BDA0002611819560000108
作为初始位移响应量。根据式(10)可计算边坡监测点k布设后第1步开挖完成时的初始动力增载位移响应率
Figure BDA0002611819560000109
Figure BDA0002611819560000111
将边坡监测点k布设后第i步开挖完成时的下滑动力增量ΔTi k作为该开挖步的滑移动力增载量,边坡开挖引起的对应位移
Figure BDA0002611819560000112
作为该开挖步的位移响应量。根据式(11)可计算边坡监测点k布设后第i步开挖完成时的动力增载位移响应率
Figure BDA0002611819560000113
Figure BDA0002611819560000114
式(10)、(11)的计算结果见表6。
表6监测点k布设后各开挖步的动力增载响应率
Figure BDA0002611819560000115
步骤五:矿山边坡开挖动力增载位移响应比的确定
根据步骤四的计算结果,将矿山边坡监测点k布设后第i步开挖完成时的动力增载位移响应率
Figure BDA0002611819560000116
与初始动力增载位移响应率
Figure BDA0002611819560000117
之比定义为监测点k的动力增载位移响应比,由此建立矿山边坡动力增载位移响应比评价模型,即:
Figure BDA0002611819560000118
计算结果见表7。
表7监测点k布设后各开挖步的动力增载响应比
Figure BDA0002611819560000119
Figure BDA0002611819560000121
将矿山边坡开挖动力增载位移响应比作为评价参数应用到边坡稳定性判定中,对边坡进行稳定性评价与监测预警。
步骤六:基于矿山边坡开挖动力增载位移响应比的失效概率评价指标的确定
对矿山边坡监测点k在布设后第i步开挖完成时,分别根据式(13)、(14)计算其动力增载位移响应比的均值
Figure BDA0002611819560000122
与标准差
Figure BDA0002611819560000123
计算结果见表8和表9:
Figure BDA0002611819560000124
Figure BDA0002611819560000125
表8监测点k的动力增载响应比均值
Figure BDA0002611819560000126
表9监测点k的动力增载响应比标准差
Figure BDA0002611819560000127
根据式(15)计算矿山边坡监测点k的动力增载位移响应比可靠度指标
Figure BDA0002611819560000128
计算结果见
表10:
Figure BDA0002611819560000131
表10监测点k的动力增载响应比可靠度
Figure BDA0002611819560000132
根据式(16)计算矿山边坡监测点k布设后第i步开挖完成时可能发生破坏的概率,即失效概率
Figure BDA0002611819560000133
计算结果见表11:
Figure BDA0002611819560000134
式中,
Figure BDA0002611819560000135
为标准正态分布函数,其值可查询标准正态函数分布取值表。
表11各监测点失效概率
Figure BDA0002611819560000136
步骤七:矿山边坡稳定性评价及监测预警
将失效概率
Figure BDA0002611819560000137
作为矿山边坡稳定性的监测预警与评价参数,根据边坡失效概率的大小实时判定矿山边坡开挖过程中的稳定性状态,并对边坡稳定性等级进行划分与监测预警:
当监测点k各开挖层的
Figure BDA0002611819560000138
均小于等于0.3时,判定边坡处于稳定状态,可继续进行下一个台阶的开挖;当监测点k第i开挖层的
Figure BDA0002611819560000139
为0.3~0.6时,判定边坡处于欠稳定状态,此时应发布一级预警,同时对第i层的台阶坡面进行补充勘察与初步加固处理等;当监测点k第i 开挖层的
Figure BDA00026118195600001310
大于0.6时,判定边坡稳定性较差,随时可能发生失稳破坏,此时应发布二级预警,并及时开挖边坡进行系统勘察与加固、支护等,必要时刻及时撤离人员及设备至安全区域。
本实施例中,监测点1~7的失效概率虽然整体呈降低趋势,但
Figure BDA0002611819560000141
均大于0.6,可判定该矿山边坡处于不稳定状态,随时可能发生失稳破坏,此时应发布二级预警,并及时开挖边坡进行系统勘察与加固、支护等,必要时刻及时撤离人员及设备至安全区域。
本发明根据莫尔-库伦强度理论、弹塑性力学和数理统计基本原理,提出了一种基于开挖方量和位移实时监测与分析的深凹矿山边坡稳定性测定新方法。该方法可有效克服传统极限平衡法无法考虑土体应力应变关系等方面的局限,同时又解决了传统位移时序预测方法无法定量描述边坡动力作用机理、无稳定的统一失稳判据等问题,采用的判据可以准确的确定边坡在开挖过程中不同阶段的稳定性状态,可为开挖型边坡的稳定性评价和治理提供有效依据。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作其它形式的限制,任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例应用于其它领域,但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型,仍属于本发明技术方案的保护范围。

Claims (7)

1.一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法,其特征在于,包括以下步骤:
步骤一、确定深凹矿山边坡特征参数,并选取开挖位移监测点,所述矿山边坡特征参数包括边坡开挖台阶数量为m,各开挖步高度hi,各开挖坡角θi
步骤二、实时监测边坡位移,得到边坡开挖方量;
步骤三、根据步骤二中得到的边坡开挖方量,确定边坡开挖方量与边坡滑移动力增载量之间的关系;
步骤四、根据步骤二实时监测的边坡位移及步骤三得到的滑移动力增载量确定矿山边坡开挖动力增载位移响应率;
步骤五、根据步骤四得到的结果确定矿山边坡开挖动力增载位移响应比,由此建立矿山边坡动力增载位移响应比评价模型;
步骤六、将矿山边坡开挖动力增载位移响应比作为评价参数应用到边坡稳定性判定中,确定基于矿山边坡开挖动力增载位移响应比的失效概率评价指标;
步骤七、进而根据步骤六得到的失效概率实现对矿山边坡稳定性的测定。
2.根据权利要求1所述的基于挖方量监测的深凹矿山开采边坡稳定性测定方法,其特征在于:所述步骤二具体包括以下步骤:
(1)分析获得矿山边坡开挖过程中各监测点k布设后第i步开挖完成时的位移变化量
Figure FDA0002611819550000011
ΔSi k=Si k-Si-1 k (1)
式中:
Figure FDA0002611819550000012
为监测点k布设后第i-1步开挖完成时的合位移,k=1,2,…,m,
Figure FDA0002611819550000013
监测点k布设后第i步开挖完成时的合位移;
(2)得到监测点k布设后第i开挖步开挖方量:
Figure FDA0002611819550000014
式中:
Figure FDA0002611819550000015
为监测点k布设后第i开挖步上、下断面面积;
Figure FDA0002611819550000016
作为各开挖步的开挖影响范围,
Figure FDA0002611819550000017
θi为第i开挖步的坡角;d为坡体厚度,Li为开挖高度,Li=hi
3.根据权利要求2所述的基于挖方量监测的深凹矿山开采边坡稳定性测定方法,其特征在于:所述步骤三具体包括以下步骤:
(1)确定矿山边坡开挖方量与开挖坡角的关系:
Figure FDA0002611819550000021
(2)确定矿山边坡开挖方量与开挖坡角均值的关系:
Figure FDA0002611819550000022
式中:
Figure FDA0002611819550000023
为监测点k布设后第1-i步开挖完成时,整体开挖方量的上、下断面面积;Li(1-i)为1-i步整体开挖高度;
Figure FDA0002611819550000024
a为监测点k布设后第1-i步整体开挖方量的影响范围;a为开挖台阶宽度,
Figure FDA0002611819550000025
为开挖坡角的均值,
Figure FDA0002611819550000026
(3)确定边坡滑移动力增载量与开挖坡角均值的关系,进而得到边坡开挖方量与边坡滑移动力增载量之间的关系;
其中,将监测点k布设后第1-i开挖步的稳定性系数降低变化量ΔFsi k与开挖体积Vi k的乘积作为该开挖步的滑移动力增载量ΔTi k,即:
Figure FDA0002611819550000027
得到开挖坡角θi、开挖坡角均值
Figure FDA0002611819550000028
与滑移动力增载量成正相关关系。
4.根据权利要求3所述的基于挖方量监测的深凹矿山开采边坡稳定性测定方法,其特征在于:所述步骤四中,边坡开挖动力增载位移响应率采用以下方式确定:
将边坡监测点k布设后第i步开挖完成时的下滑动力增量ΔTi k作为该开挖步的滑移动力增载量,边坡开挖引起的对应位移
Figure FDA0002611819550000029
作为该开挖步的位移响应量,进而得到边坡监测点k布设后第i步开挖完成时的动力增载位移响应率
Figure FDA00026118195500000210
Figure FDA00026118195500000211
5.根据权利要求1所述的基于挖方量监测的深凹矿山开采边坡稳定性测定方法,其特征在于:所述步骤一中,开挖位移监测点的位置采用以下选取方式:
(1)矿山边坡开挖前,在第一个开挖台阶临空面拐角处设置一个初始位移监测点;
(2)此后,在每开挖到一个台阶时则在新的台阶开挖临空面拐角处布设一个位移监测点,并设布设的监测点为k(k=1,2,…,m)。
6.根据权利要求4所述的基于挖方量监测的深凹矿山开采边坡稳定性测定方法,其特征在于:所述步骤五中,将矿山边坡监测点k布设后第i步开挖完成时的动力增载位移响应率
Figure FDA0002611819550000031
与初始动力增载位移响应率
Figure FDA0002611819550000032
之比定义为监测点k的动力增载位移响应比,由此建立的矿山边坡动力增载位移响应比评价模型为:
Figure FDA0002611819550000033
7.根据权利要求6所述的基于挖方量监测的深凹矿山开采边坡稳定性测定方法,其特征在于:所述步骤六在确定基于矿山边坡开挖动力增载位移响应比的失效概率评价指标时,具体采用以下方式:
(1)对矿山边坡监测点k在布设后第i步开挖完成时,分别计算其动力增载位移响应比的均值
Figure FDA0002611819550000034
与标准差
Figure FDA0002611819550000035
Figure FDA0002611819550000036
Figure FDA0002611819550000037
(2)获得矿山边坡监测点k的动力增载位移响应比可靠度指标
Figure FDA0002611819550000038
Figure FDA0002611819550000039
(3)得到矿山边坡监测点k布设后第i步开挖完成时发生破坏的概率,即失效概率
Figure FDA00026118195500000310
Figure FDA00026118195500000311
式中,
Figure FDA00026118195500000312
为标准正态分布函数。
CN202010756699.XA 2020-07-31 2020-07-31 一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法 Active CN111912953B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010756699.XA CN111912953B (zh) 2020-07-31 2020-07-31 一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010756699.XA CN111912953B (zh) 2020-07-31 2020-07-31 一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法

Publications (2)

Publication Number Publication Date
CN111912953A true CN111912953A (zh) 2020-11-10
CN111912953B CN111912953B (zh) 2022-08-19

Family

ID=73287121

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010756699.XA Active CN111912953B (zh) 2020-07-31 2020-07-31 一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法

Country Status (1)

Country Link
CN (1) CN111912953B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113239510A (zh) * 2021-03-18 2021-08-10 昆明理工大学 一种矿山边坡最优开挖坡率计算方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103103974A (zh) * 2013-02-02 2013-05-15 青岛理工大学 利用监测位移方向性参数检测边坡稳定性的方法
CN104535346A (zh) * 2014-12-17 2015-04-22 青岛理工大学 一种混凝土重力危坝的抗滑稳定性检测方法
JP2016023443A (ja) * 2014-07-17 2016-02-08 首都高速道路株式会社 橋梁の制振構造及び橋梁の制振構造の設定方法
CN105606063A (zh) * 2016-01-27 2016-05-25 青岛理工大学 一种基于正交应变比的土层边坡稳定性测定方法
CN108867666A (zh) * 2018-05-25 2018-11-23 重庆交通大学 基于开挖变形的结构面控制边坡稳定性评价方法
CN109680669A (zh) * 2018-12-18 2019-04-26 青岛理工大学 一种基于地下水动力加载的基坑稳定性监测预警方法
CN109785584A (zh) * 2019-01-29 2019-05-21 青岛理工大学 复合水动力库岸边坡稳定性预测方法
CN109811809A (zh) * 2019-03-20 2019-05-28 四川路桥华东建设有限责任公司 一种高边坡与锚碇基坑开挖相互影响机制的模型试验系统及方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103103974A (zh) * 2013-02-02 2013-05-15 青岛理工大学 利用监测位移方向性参数检测边坡稳定性的方法
JP2016023443A (ja) * 2014-07-17 2016-02-08 首都高速道路株式会社 橋梁の制振構造及び橋梁の制振構造の設定方法
CN104535346A (zh) * 2014-12-17 2015-04-22 青岛理工大学 一种混凝土重力危坝的抗滑稳定性检测方法
CN105606063A (zh) * 2016-01-27 2016-05-25 青岛理工大学 一种基于正交应变比的土层边坡稳定性测定方法
CN108867666A (zh) * 2018-05-25 2018-11-23 重庆交通大学 基于开挖变形的结构面控制边坡稳定性评价方法
CN109680669A (zh) * 2018-12-18 2019-04-26 青岛理工大学 一种基于地下水动力加载的基坑稳定性监测预警方法
CN109785584A (zh) * 2019-01-29 2019-05-21 青岛理工大学 复合水动力库岸边坡稳定性预测方法
CN109811809A (zh) * 2019-03-20 2019-05-28 四川路桥华东建设有限责任公司 一种高边坡与锚碇基坑开挖相互影响机制的模型试验系统及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HE KEQIANG 等: "Unload–load displacement response ratio parameter and its application in prediction of debris landslide induced by rainfall", 《ENVIRONMENTAL EARTH SCIENCE》, vol. 76, no. 1, 3 January 2017 (2017-01-03), pages 1 - 16, XP036203381, DOI: 10.1007/s12665-016-6372-0 *
杨国俊等: "某水电站导流洞入口边坡开挖稳定性分析", 《安全与环境工程》, vol. 19, no. 05, 30 September 2012 (2012-09-30), pages 111 - 114 *
贺可强: "复杂堆积层滑坡的稳定性评价方法分析与展望", 《青岛理工大学学报》, vol. 37, no. 01, 31 December 2016 (2016-12-31), pages 1 - 9 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113239510A (zh) * 2021-03-18 2021-08-10 昆明理工大学 一种矿山边坡最优开挖坡率计算方法
CN113239510B (zh) * 2021-03-18 2022-05-20 昆明理工大学 一种矿山边坡最优开挖坡率计算方法

Also Published As

Publication number Publication date
CN111912953B (zh) 2022-08-19

Similar Documents

Publication Publication Date Title
Zhou et al. Research on theory and technology of floor heave control in semicoal rock roadway: Taking longhu coal mine in Qitaihe mining area as an Example
CN105260575A (zh) 一种基于神经网络的巷道围岩变形预测方法
Madani Esfahani et al. Fault detection in 3D by sequential Gaussian simulation of Rock Quality Designation (RQD) Case study: Gazestan phosphate ore deposit, Central Iran
CN103225509B (zh) 一种大采高综采采场顶板分类与支护阻力确定方法
CN110298107B (zh) 一种基于增量叠加的工作面冲击危险性评价方法
CN110595887B (zh) 一种浅埋硬岩大型地下洞室成洞条件量化评价方法
Hu et al. Study on failure depth of coal seam floor in deep mining
CN105653811B (zh) 高压变电站深回填区旋挖钻进中的入岩深度确定方法
CN106339798A (zh) 一种基于模糊综合评判法的软岩斜井井筒支护体系稳定性评价方法
CN108665146B (zh) 大型露天矿山边坡稳定性评价精度确定方法
CN111912953B (zh) 一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法
Li et al. Prediction model of fracture depth and water inrush risk zoning in deep mining coal seam floor
CN111324960A (zh) 一种含有浅埋隐伏球状空洞地基承载力的量化方法
CN106777629B (zh) 一种基于工程边界条件的膨胀土计算模型及其应用方法
CN112948924A (zh) 基于覆岩结构的近松散层开采导水裂隙带高度确定方法
CN117371068A (zh) 一种对bif铁矿顺层边坡变形破坏的数值模拟分析方法
CN104632207A (zh) 测定钻孔灌注桩孔壁致塌影响因素的方法
Bi et al. Water abundance comprehensive evaluation of coal mine aquifer based on projection pursuit model
CN106407569A (zh) 一种厚松散层薄基岩条件下地表下沉值计算方法
CN111598366A (zh) 一种实时钻井辅助决策方法及系统
CN108416174A (zh) 大型露天矿山边坡设计安全系数取值方法
FANG et al. Intelligent recognition of tunnel stratum based on advanced drilling tests
CN113804853A (zh) 煤矿瓦斯含量检测方法及系统
Chao et al. A prediction of borehole stability while drilling preliminary prospecting wells based on seismic impedance
CN110472273B (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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: He Keqiang

Inventor after: Wu Lili

Inventor after: Bao Fangge

Inventor after: Sun Linna

Inventor after: Gao Ming

Inventor before: He Keqiang

Inventor before: Bao Fangge

Inventor before: Wu Lili

Inventor before: Sun Linna

Inventor before: Gao Ming

GR01 Patent grant
GR01 Patent grant