CN114757028A - 一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法 - Google Patents

一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法 Download PDF

Info

Publication number
CN114757028A
CN114757028A CN202210394632.5A CN202210394632A CN114757028A CN 114757028 A CN114757028 A CN 114757028A CN 202210394632 A CN202210394632 A CN 202210394632A CN 114757028 A CN114757028 A CN 114757028A
Authority
CN
China
Prior art keywords
slope
angle
length
interlayer
section
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
CN202210394632.5A
Other languages
English (en)
Other versions
CN114757028B (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.)
Liaoning Technical University
Original Assignee
Liaoning Technical University
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 Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN202210394632.5A priority Critical patent/CN114757028B/zh
Publication of CN114757028A publication Critical patent/CN114757028A/zh
Application granted granted Critical
Publication of CN114757028B publication Critical patent/CN114757028B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/23—Dune restoration or creation; Cliff stabilisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Geometry (AREA)
  • Algebra (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
  • Tires In General (AREA)

Abstract

本发明公开一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法，方法首先垂直边坡走向方向和垂直软弱夹层走向方向分别做剖面，测量两剖面的相关参数，同时获取边坡岩土力学参数；然后作辅助线，选取剖面上的单个条块由极限平衡法得到边坡稳定性系数；接着将剖面旋转一定的角度形成新的剖面，分析新的剖面中对应的条块的边坡稳定性系数；最后将两稳定性系数的差值定义为边坡稳定性差函数，绘制边坡稳定性差函数与旋转角度之间的关系曲线，找到极值点，将极值位置对应的剖面作为边坡主滑方向。本发明的方法能够定量分析含顺倾软弱夹层边坡主滑方向，能够为边坡稳定性分析与形态优化提供理论基础，同时可为露天煤矿安全高效生产提供强有力保障。

Description

一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法
技术领域
本发明涉及露天煤矿开采领域，尤其涉及一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法。
背景技术
随着我国露天煤矿产量规模的不断增大和数量的不断增加，边坡稳定性问题愈来愈突出，露天煤矿含顺倾软弱夹层边坡滑坡的现象时常发生，且滑体规模大，极难防治，造成的损失极为惨重，严重限制着露天煤矿安全、高效生产。在以往边坡稳定性研究中，当软弱夹层走向方向与边坡走向方向平行时，边坡的主滑方向为垂直边坡走向；当软弱夹层走向方向与边坡走向方向不平行时，边坡主滑方向的确定方法尚未建立。
在以往边坡稳定性研究中，当软弱夹层为水平或近水平时，选取垂直于边坡面作为研究问题的剖面位置；但当软弱夹层为倾斜时，选取计算剖面位置则是改为软弱夹层。在垂直边坡切制剖面中，边坡角较大，软弱夹层倾角较缓，而在垂直软弱夹层切制剖面中边坡角较小，软弱夹层倾角较大，所以在两剖面位置中间必定存在一剖面位置，边坡稳定系数最小，此种位置则可作为露天煤矿含顺倾软弱夹层边坡主滑方向。
因此，需要提出一种新的针对含软弱夹层边坡稳定性主滑方向确定的方法，用于解决此种条件下边坡主滑方向的确定，为露天煤矿含顺倾软弱夹层边坡稳定性分析提供理论基础。
发明内容
针对上述现有技术的不足，本发明提供一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法。
为解决上述技术问题，本发明所采取的技术方案是：一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法，包括如下步骤：
步骤1：垂直于边坡走向方向切割剖面odb，垂直于软弱夹层走向方向切割剖面od'b'，两剖面与水平面的交线ob和ob'相交于点o，测量∠bob'的度数；
步骤2：作剖面odb与软弱夹层交线为oc，剖面od'b'与软弱夹层交线为oc'；
步骤3：测量剖面dob处的参数，包括边坡高度H，即线段db的长度；坡顶到软弱夹层的高度H1，即线段dc的长度；边坡角A的度数，即∠dob的度数作为已知量；
步骤4：测量剖面d'ob'处软弱夹层倾角D的度数，即∠c'ob'的度数作为已知量；
步骤5：获取边坡岩土力学参数；
进一步的，所述边坡岩土力学参数包括：容重γ、软弱夹层内摩擦角
Figure BDA0003598336840000011
和软弱夹层内聚力c。
步骤6：将bd沿着bo方向移动单位长度，与ob、oc和od分别交于点a、h和f，通过点a作线段bb'的平行线与线段ob'交于a'点，通过点a'作d'b'的平行线，与oc'和od'分别交于点h'和f'；
步骤7：以剖面odb中第一个条块dfhc为研究对象，由极限平衡法稳定系数表达公式得到边坡稳定性系数Fs1；
步骤7.1：已知边坡高度H，根据Δdob三角函数关系求得线段ob的长度，记为L1；
步骤7.2：由Δcob三角函数关系得到(H-H1)/L1＝tanB，通过反三角函数，得出剖面odb处软弱夹层倾角的角度B；
步骤7.3：分别过点c和d作平行于ab的辅助线，与af和af的延长线分别交于点g和点e，由Δghc和Δdef三角函数关系、L1的长度、角度B和边坡角A，求得条块dfhc中ef和gh的高度，分别记为H3和H4；
步骤7.4：计算条块左侧fg的高度H2＝H1+H3-H4；
步骤7.5：条块dfhc为不规则平行四边形，其面积通过海伦公式进行求解，设中间系数为Z1，求得条块dfhc的面积
Figure BDA0003598336840000021
其中，L3、L4分别为条块dfhc的顶宽fd和底宽hc的长度；
步骤7.6：由极限平衡法稳定系数表达公式计算边坡稳定性数Fs1，具体公式如下：
Figure BDA0003598336840000022
式中：
Figure BDA0003598336840000023
Figure BDA0003598336840000024
Figure BDA0003598336840000031
步骤8：将剖面odb旋转角度ωi后与面dbb'd'交于dibi，形成新的剖面odibi，与软弱夹层交线为oci；以剖面odibi中第一个条块difihici为研究对象，由极限平衡法稳定系数表达公式得到边坡稳定性系数Fsi1，过程如下：
步骤8.1：在Δbobi中，由三角函数关系和线段ob的长度L1，以及旋转角度ωi后边坡角度Ai求得Li1长度；
步骤8.2：在Δciobi中，已知obi的长度Li1，db＝dibi，由三角函数关系，求得旋转角度ωi后剖面odibi处的软弱夹层倾角Bi的值；
步骤8.3：在Δbobi中，已知ob的长度L1，由三角函数关系，求得旋转角度ωi后条块difihici的宽度aibi，记为Li5；
步骤8.4：分别过点ci和di作平行于aibi的辅助线，与aifi和aifi的延长线分别交于点gi和点ei，由Δgihici和Δdieifi三角函数关系，Li1长度，旋转ωi后边坡角Ai和软弱夹层倾角角度Bi，求得条块difihici中的eifi和gihi高度，分别记为Hi3和Hi4；
步骤8.5：旋转后条块面积同样通过海伦公式进行求解，设中间系数为Zi＝(Hi1+Li3+Hi2+Li4)/2，可求得条块面积
Figure BDA0003598336840000032
Figure BDA0003598336840000033
式中：
Figure BDA0003598336840000034
Figure BDA0003598336840000041
Figure BDA0003598336840000042
Figure BDA0003598336840000043
其中，Li3、Li4分别为条块difihici的顶宽fidi和底宽hici的长度，L2为ob'的长度；
步骤8.6：由极限平衡法稳定系数表达公式可得边坡稳定性数Fsi1：
Figure BDA0003598336840000044
式中：
Figure BDA0003598336840000045
Figure BDA0003598336840000046
Figure BDA0003598336840000047
Figure BDA0003598336840000048
Figure BDA0003598336840000049
Figure BDA0003598336840000051
Figure BDA0003598336840000052
Figure BDA0003598336840000053
Figure BDA0003598336840000054
ξ10＝tan2 A cos2ωi
Figure BDA0003598336840000055
Figure BDA0003598336840000056
ξ13＝H1-H+ξ15
ξ14＝L2 tan Atanω
Figure BDA0003598336840000057
进一步的，所述旋转ωi后边坡角Ai和软弱夹层倾角角度Bi的计算过程如下：
由Δbob'三角函数关系，求得线段bb'的长度
Figure BDA0003598336840000058
由Δc'ob'三角函数关系，求出d'c'的高度
Figure BDA0003598336840000059
将cc'与bb'之间的夹角记为θ，则
Figure BDA00035983368400000510
计算cibi的长度
Figure BDA00035983368400000511
在Δbobi中，由三角函数关系、L1的长度和ωi角度，可求得obi的长度Li1：
Figure BDA00035983368400000512
在Δdiobi中，由三角函数关系，并结合所求得obi的长度Lii，得到旋转后边坡角度Ai：
Ai＝arctan(tan Acosωi)
在Δciobi中，由三角函数关系，并结合cibi的长度Z和obi的长度Lii，得到旋转后软弱夹层倾角Bi：
Figure BDA0003598336840000061
步骤9：将Fs1i与Fs1之差定义为边坡稳定性差函数ΔFs；
所述边坡稳定性差函数ΔFs的具体公式如下：
Figure BDA0003598336840000062
式中：
Figure BDA0003598336840000063
Figure BDA0003598336840000064
Figure BDA0003598336840000065
Figure BDA0003598336840000066
Figure BDA0003598336840000067
Figure BDA0003598336840000071
Figure BDA0003598336840000072
Figure BDA0003598336840000073
Figure BDA0003598336840000074
λ10＝H-H1+λ6。
步骤10：绘制边坡稳定性差函数ΔFs与ωi的曲线，曲线的极值处边坡的稳定系数最小，将此处位置对应的剖面为边坡主滑方向。
采用上述技术方案所产生的有益效果在于：本发明提供的方法考虑到顺倾软弱夹层对不同剖面位置边坡稳定性的影响，提出边坡稳定性差函数，并对公式基于Matlab进行求解和化简，弥补以往在相似工程背景边坡稳定性计算、分析和边坡形态设计时，仅考虑垂直边坡或垂直软弱夹层的边坡稳定性分析，造成计算结果没能真实反映边坡的稳定性状态，影响边坡形态设计或治理方案制定等工作。本发明提出的方法为能够为边坡稳定性分析与形态优化提供理论基础，同时可为露天煤矿安全高效生产提供强有力保障。
附图说明
图1为本发明实施例中确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法的流程图；
图2为本发明实施例中两剖面选取位置的示意图；
图3为本发明实施例中基于两剖面绘制的几何形态图；
图4为本发明实施例中图3是左视图；
图5为本发明实施例中边坡稳定性差函数ΔFs的曲线图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
本实施例以位于内蒙古内的某露天煤矿为例，该露天煤矿边坡岩体同样为软岩倾斜软弱夹层，垂直边坡的剖面中，边坡高度H＝120m，边坡角A为16°，坡顶到软弱夹层的高度H1为60m；垂直软弱夹层的剖面中，软弱夹层倾角D为13°，坡顶到软弱夹层的高度H5为30m；两剖面线在水平面的夹角α为37°。边坡岩体力学参数为：容重γ＝19kN/m3，软弱夹层内摩擦角
Figure BDA0003598336840000081
软弱夹层内聚力c＝9.18kpa。
如图1所示，本实施例中一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法如下所述。
步骤1：垂直于边坡走向方向切割剖面odb，垂直于软弱夹层走向方向切割剖面od'b'，两剖面与水平面的交线ob和ob'相交于点o，测量∠bob'的度数；
本实施例中，两剖面选取位置如图2所示。
步骤2：作剖面odb与软弱夹层交线为oc，剖面od'b'与软弱夹层交线为oc'；
步骤3：测量剖面dob处的参数，包括边坡高度H，即线段db的长度；坡顶到软弱夹层的高度H1，即线段dc的长度；边坡角A的度数，即∠dob的度数作为已知量；
步骤4：测量剖面d'ob'处软弱夹层倾角D的度数，即∠c'ob'的度数作为已知量；
本实施例中，基于两剖面绘制的几何形态图如图3所示，图3的左视图如图4所示。
步骤5：获取边坡岩土力学参数；
进一步的，所述边坡岩土力学参数包括：容重γ、软弱夹层内摩擦角
Figure BDA0003598336840000082
和软弱夹层内聚力c。
步骤6：将bd沿着bo方向移动单位长度，与ob、oc和od分别交于点a、h和f，通过点a作线段bb'的平行线与线段ob'交于a'点，通过点a'作d'b'的平行线，与oc'和od'分别交于点h'和f'；如图3、4所示。
步骤7：以剖面odb中第一个条块dfhc为研究对象，由极限平衡法稳定系数表达公式得到边坡稳定性系数Fs1；
步骤7.1：已知边坡高度H，根据Δdob三角函数关系求得线段ob的长度，记为L1；
步骤7.2：由Δcob三角函数关系得到(H-H1)/L1＝tanB，通过反三角函数，得出剖面odb处软弱夹层倾角的角度B；
步骤7.3：分别过点c和d作平行于ab的辅助线，与af和af的延长线分别交于点g和点e，由Δghc和Δdef三角函数关系、L1的长度、角度B和边坡角A，求得条块dfhc中ef和gh的高度，分别记为H3和H4；辅助线的做法如图3、4所示。
步骤7.4：计算条块左侧fg的高度H2＝H1+H3-H4；
步骤7.5：条块dfhc为不规则平行四边形，其面积通过海伦公式进行求解，设中间系数为Z1，求得条块dfhc的面积
Figure BDA0003598336840000091
其中，L3、L4分别为条块dfhc的顶宽fd和底宽hc的长度；
步骤7.6：由极限平衡法稳定系数表达公式计算边坡稳定性数Fs1，具体公式如下：
Figure BDA0003598336840000092
式中：
Figure BDA0003598336840000093
Figure BDA0003598336840000094
Figure BDA0003598336840000095
步骤8：将剖面odb旋转角度ωi后与面dbb'd'交于dibi，形成新的剖面odibi，与软弱夹层交线为oci，旋转后剖面位置如图3、4所示；以剖面odibi中第一个条块difihici为研究对象，由极限平衡法稳定系数表达公式得到边坡稳定性系数Fsi1，过程如下：
步骤8.1：在Δbobi中，由三角函数关系和线段ob的长度L1，以及旋转角度ωi后边坡角度Ai求得Li1长度；
步骤8.2：在Δciobi中，已知obi的长度Li1，db＝dibi，由三角函数关系，求得旋转角度ωi后剖面odibi处的软弱夹层倾角Bi的值；
步骤8.3：在Δbobi中，已知ob的长度L1，由三角函数关系，求得旋转角度ωi后条块difihici的宽度aibi，记为Li5；
步骤8.4：分别过点ci和di作平行于aibi的辅助线，与aifi和aifi的延长线分别交于点gi和点ei，由Δgihici和Δdieifi三角函数关系，Li1长度，旋转ωi后边坡角Ai和软弱夹层倾角角度Bi，求得条块difihici中的eifi和gihi高度，分别记为Hi3和Hi4；
步骤8.5：旋转后条块面积同样通过海伦公式进行求解，设中间系数为Zi＝(Hi1+Li3+Hi2+Li4)/2，可求得条块面积
Figure BDA0003598336840000101
Figure BDA0003598336840000102
式中：
Figure BDA0003598336840000103
Figure BDA0003598336840000104
Figure BDA0003598336840000105
Figure BDA0003598336840000106
其中，Li3、Li4分别为条块difihici的顶宽fidi和底宽hici的长度，L2为ob'的长度；
步骤8.6：由极限平衡法稳定系数表达公式可得边坡稳定性数Fsi1：
Figure BDA0003598336840000107
式中：
Figure BDA0003598336840000111
Figure BDA0003598336840000112
Figure BDA0003598336840000113
Figure BDA0003598336840000114
Figure BDA0003598336840000115
Figure BDA0003598336840000116
Figure BDA0003598336840000117
Figure BDA0003598336840000118
Figure BDA0003598336840000119
ξ10＝tan2 A cos2ωi
Figure BDA00035983368400001110
Figure BDA00035983368400001111
ξ13＝H1-H+ξ15
ξ14＝L2 tan Atanω
Figure BDA00035983368400001112
进一步的，所述旋转ωi后边坡角Ai和软弱夹层倾角角度Bi的计算过程如下：
由Δbob'三角函数关系，求得线段bb'的长度
Figure BDA0003598336840000121
由Δc'ob'三角函数关系，求出d'c'的高度
Figure BDA0003598336840000122
将cc'与bb'之间的夹角记为θ，则
Figure BDA0003598336840000123
计算cibi的长度
Figure BDA0003598336840000124
在Δbobi中，由三角函数关系、L1的长度和ωi角度，可求得obi的长度Li1：
Figure BDA0003598336840000125
在Δdiobi中，由三角函数关系，并结合所求得obi的长度Lii，得到旋转后边坡角度Ai：
Ai＝arctan(tan Acosωi)
在Δciobi中，由三角函数关系，并结合cibi的长度Z和obi的长度Lii，得到旋转后软弱夹层倾角Bi：
Figure BDA0003598336840000126
步骤9：将Fs1i与Fs1之差定义为边坡稳定性差函数ΔFs；
所述边坡稳定性差函数ΔFs的具体公式如下：
Figure BDA0003598336840000127
式中：
Figure BDA0003598336840000131
Figure BDA0003598336840000132
Figure BDA0003598336840000133
Figure BDA0003598336840000134
Figure BDA0003598336840000135
Figure BDA0003598336840000136
Figure BDA0003598336840000137
Figure BDA0003598336840000138
Figure BDA0003598336840000139
λ10＝H-H1+λ6。
步骤10：绘制边坡稳定性差函数ΔFs与ωi的曲线，曲线的极值处边坡的稳定系数最小，将此处位置对应的剖面为边坡主滑方向。
本实施例中，将数据代入所求得边坡稳定性差函数ΔFs如下所示：
Figure BDA00035983368400001310
式中：
Figure BDA0003598336840000141
Figure BDA0003598336840000142
Figure BDA0003598336840000143
Figure BDA0003598336840000144
Figure BDA0003598336840000145
Figure BDA0003598336840000146
Figure BDA0003598336840000147
Figure BDA0003598336840000148
Figure BDA0003598336840000149
Figure BDA00035983368400001410
绘制边坡稳定性差函数ΔFs的曲线如图5所示，该曲线在点(29,0.1949)处出现极值，说明当剖面odb旋转角度ωi＝29°位置时，边坡稳定系数最小，因此将此处位置剖面作为研究此边坡稳定性主滑方向，为边坡稳定性评价与设计提供基础。

Claims (6)

1.一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法，其特征在于，包括如下步骤：
步骤1：垂直于边坡走向方向切割剖面odb，垂直于软弱夹层走向方向切割剖面od'b'，两剖面与水平面的交线ob和ob'相交于点o，测量∠bob'的度数；
步骤2：作剖面odb与软弱夹层交线为oc，剖面od'b'与软弱夹层交线为oc'；
步骤3：测量剖面dob处的参数，包括边坡高度H，即线段db的长度；坡顶到软弱夹层的高度H1，即线段dc的长度；边坡角A的度数，即∠dob的度数作为已知量；
步骤4：测量剖面d'ob'处软弱夹层倾角D的度数，即∠c'ob'的度数作为已知量；
步骤5：获取边坡岩土力学参数；
步骤6：将bd沿着bo方向移动单位长度，与ob、oc和od分别交于点a、h和f，通过点a作线段bb'的平行线与线段ob'交于a'点，通过点a'作d'b'的平行线，与oc'和od'分别交于点h'和f'；
步骤7：以剖面odb中第一个条块dfhc为研究对象，由极限平衡法稳定系数表达公式得到边坡稳定性系数Fs1；
步骤8：将剖面odb旋转角度ωi后与面dbb'd'交于dibi，形成新的剖面odibi，与软弱夹层交线为oci；以剖面odibi中第一个条块difihici为研究对象，由极限平衡法稳定系数表达公式得到边坡稳定性系数Fsi1；
步骤9：将Fs1i与Fs1之差定义为边坡稳定性差函数ΔFs；
步骤10：绘制边坡稳定性差函数ΔFs与ωi的曲线，曲线的极值处边坡的稳定系数最小，将此处位置对应的剖面为边坡主滑方向。
2.根据权利要求1所述的一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法，其特征在于，所述边坡岩土力学参数包括：容重γ、软弱夹层内摩擦角
Figure FDA0003598336830000011
和软弱夹层内聚力c。
3.根据权利要求1所述的一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法，其特征在于，所述步骤7的具体过程如下：
步骤7.1：已知边坡高度H，根据Δdob三角函数关系求得线段ob的长度，记为L1；
步骤7.2：由Δcob三角函数关系得到(H-H1)/L1＝tanB，通过反三角函数，得出剖面odb处软弱夹层倾角的角度B；
步骤7.3：分别过点c和d作平行于ab的辅助线，与af和af的延长线分别交于点g和点e，由Δghc和Δdef三角函数关系、L1的长度、角度B和边坡角A，求得条块dfhc中ef和gh的高度，分别记为H3和H4；
步骤7.4：计算条块左侧fg的高度H2＝H1+H3-H4；
步骤7.5：条块dfhc为不规则平行四边形，其面积通过海伦公式进行求解，设中间系数为Z1，求得条块dfhc的面积
Figure FDA0003598336830000021
其中，L3、L4分别为条块dfhc的顶宽fd和底宽hc的长度；
步骤7.6：由极限平衡法稳定系数表达公式计算边坡稳定性数Fs1，具体公式如下：
Figure FDA0003598336830000022
式中：
Figure FDA0003598336830000023
Figure FDA0003598336830000024
Figure FDA0003598336830000025
4.根据权利要求1所述的一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法，其特征在于，所述步骤8的具体过程如下：
步骤8.1：在Δbobi中，由三角函数关系和线段ob的长度L1，以及旋转角度ωi后边坡角度Ai求得Li1长度；
步骤8.2：在Δciobi中，已知obi的长度Li1，db＝dibi，由三角函数关系，求得旋转角度ωi后剖面odibi处的软弱夹层倾角Bi的值；
步骤8.3：在Δbobi中，已知ob的长度L1，由三角函数关系，求得旋转角度ωi后条块difihici的宽度aibi，记为Li5；
步骤8.4：分别过点ci和di作平行于aibi的辅助线，与aifi和aifi的延长线分别交于点gi和点ei，由Δgihici和Δdieifi三角函数关系，Li1长度，旋转ωi后边坡角Ai和软弱夹层倾角角度Bi，求得条块difihici中的eifi和gihi高度，分别记为Hi3和Hi4；
步骤8.5：旋转后条块面积同样通过海伦公式进行求解，设中间系数为Zi＝(Hi1+Li3+Hi2+Li4)/2，可求得条块面积
Figure FDA0003598336830000031
Figure FDA0003598336830000032
式中：
Figure FDA0003598336830000033
Figure FDA0003598336830000034
Figure FDA0003598336830000035
Figure FDA0003598336830000036
其中，Li3、Li4分别为条块difihici的顶宽fidi和底宽hici的长度，L2为ob'的长度；
步骤8.6：由极限平衡法稳定系数表达公式可得边坡稳定性数Fsi1：
Figure FDA0003598336830000041
式中：
Figure FDA0003598336830000042
Figure FDA0003598336830000043
Figure FDA0003598336830000044
Figure FDA0003598336830000045
Figure FDA0003598336830000046
Figure FDA0003598336830000047
Figure FDA0003598336830000048
Figure FDA0003598336830000049
Figure FDA00035983368300000410
ξ10＝tan2 Acos2ωi
Figure FDA00035983368300000411
Figure FDA0003598336830000051
ξ13＝H1-H+ξ15
ξ14＝L2 tanAtanω
Figure FDA0003598336830000052
5.根据权利要求6所述的一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法，其特征在于，所述旋转ωi后边坡角Ai和软弱夹层倾角角度Bi的计算过程如下：
由Δbob'三角函数关系，求得线段bb'的长度
Figure FDA0003598336830000053
由Δc'ob'三角函数关系，求出d'c'的高度
Figure FDA0003598336830000054
将cc'与bb'之间的夹角记为θ，则
Figure FDA0003598336830000055
计算cibi的长度
Figure FDA0003598336830000056
在Δbobi中，由三角函数关系、L1的长度和ωi角度，可求得obi的长度Li1：
Figure FDA0003598336830000057
在Δdiobi中，由三角函数关系，并结合所求得obi的长度Lii，得到旋转后边坡角度Ai：
Ai＝arctan(tanAcosωi)
在Δciobi中，由三角函数关系，并结合cibi的长度Z和obi的长度Lii，得到旋转后软弱夹层倾角Bi：
Figure FDA0003598336830000061
6.根据权利要求1所述的一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法，其特征在于，所述边坡稳定性差函数ΔFs的具体公式如下：
Figure FDA0003598336830000062
式中：
Figure FDA0003598336830000063
Figure FDA0003598336830000064
Figure FDA0003598336830000065
Figure FDA0003598336830000066
Figure FDA0003598336830000067
Figure FDA0003598336830000068
Figure FDA0003598336830000069
Figure FDA00035983368300000610
Figure FDA0003598336830000071
λ10＝H-H1+λ6。
CN202210394632.5A 2022-04-15 2022-04-15 一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法 Active CN114757028B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210394632.5A CN114757028B (zh) 2022-04-15 2022-04-15 一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210394632.5A CN114757028B (zh) 2022-04-15 2022-04-15 一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法

Publications (2)

Publication Number Publication Date
CN114757028A true CN114757028A (zh) 2022-07-15
CN114757028B CN114757028B (zh) 2023-12-19

Family

ID=82330694

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210394632.5A Active CN114757028B (zh) 2022-04-15 2022-04-15 一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法

Country Status (1)

Country Link
CN (1) CN114757028B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150198513A1 (en) * 2014-01-13 2015-07-16 Hubei University Of Technology Method of critical displacement forecast based on the deformation failure mechanism of slope
CN105178336A (zh) * 2015-08-18 2015-12-23 辽宁工程技术大学 一种软弱基底排土场极限承载力的确定方法
CN105787176A (zh) * 2016-02-25 2016-07-20 辽宁工程技术大学 一种露天煤矿含顺倾弱层边坡的三维稳定性计算方法
CN106677151A (zh) * 2015-11-11 2017-05-17 尤世元 一种滑移面的测定方法
CN109583135A (zh) * 2018-12-24 2019-04-05 重庆大学 针对s型沟槽填方场地边坡稳定性的极限平衡分析方法
CN109914379A (zh) * 2019-04-03 2019-06-21 昆明理工大学 一种基于条间法向力分布特征的边坡稳定性极限平衡计算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150198513A1 (en) * 2014-01-13 2015-07-16 Hubei University Of Technology Method of critical displacement forecast based on the deformation failure mechanism of slope
CN105178336A (zh) * 2015-08-18 2015-12-23 辽宁工程技术大学 一种软弱基底排土场极限承载力的确定方法
CN106677151A (zh) * 2015-11-11 2017-05-17 尤世元 一种滑移面的测定方法
CN105787176A (zh) * 2016-02-25 2016-07-20 辽宁工程技术大学 一种露天煤矿含顺倾弱层边坡的三维稳定性计算方法
CN109583135A (zh) * 2018-12-24 2019-04-05 重庆大学 针对s型沟槽填方场地边坡稳定性的极限平衡分析方法
CN109914379A (zh) * 2019-04-03 2019-06-21 昆明理工大学 一种基于条间法向力分布特征的边坡稳定性极限平衡计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
丁春健: "扎尼河露天煤矿西帮顺倾软岩边坡稳定性研究", 《中国优秀博硕士学位论文全文数据库（硕士）工程科技I辑》第（2023）02期, pages 19 - 31 *

Also Published As

Publication number Publication date
CN114757028B (zh) 2023-12-19

Similar Documents

Publication Publication Date Title
CN111750822B (zh) 一种采煤诱发的覆岩与地表沉陷协同动态预测方法
Barton TBM tunnelling in jointed and faulted rock
Krabbendam et al. Glacial erosion and bedrock properties in NW Scotland: abrasion and plucking, hardness and joint spacing
Hight et al. Some characteristics of London clay
WO2018201714A1 (zh) 一种煤矿瓦斯的精准抽采方法
CN110362955A (zh) 岩质高边坡稳定性分析三维地质力学模型试验方法及应用
Pan et al. Full-scale linear cutting test in Chongqing Sandstone and the comparison with field TBM excavation performance
Hunt Predicting the ground movements above twin tunnels constructed in London Clay
CN110424945B (zh) 页岩油水平井部署方法及装置
CN113343329B (zh) 一种参数化铁路隧道模型的建立方法
CN111101958A (zh) 一种隧道盾构区间下穿绕城施工工法
CN104809266A (zh) 一种基于样条曲线的工作面煤层赋存情况精确预测方法
CN112364474A (zh) 一种基于露天矿分区开采工艺的开采方案优化方法
WO2023000837A1 (zh) 深埋采场覆岩类型评价标准及薄基岩加厚改造设计方法
CN112417551B (zh) 3d打印采空区相似模拟建筑物的方法
CN110362899A (zh) 用于超大直径盾构过基岩突起段的刀盘扭矩计算方法
CN114757028A (zh) 一种确定露天煤矿含顺倾软弱夹层边坡主滑方向的方法
CN111767593A (zh) 一种基于地应力场类型确定井下巷道合理布置方位的方法
Cardu et al. Analysis of the tunnel boring machine advancement on the Bologna-Florence railway link
CN114961751B (zh) 一种土岩复合地层中盾构掘进引起的土体位移预测方法
Warren et al. Chalk engineering geology–channel tunnel rail link and North Downs tunnel
CN106652033A (zh) 一种地质剖面自然网格的剖分方法
CN113323658B (zh) 一种井筒漏失压力计算方法
CN116070416A (zh) 一种确定采场覆岩竖向三带发育高度的综合测试方法
CN113050161B (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