CN117371068A - 一种对bif铁矿顺层边坡变形破坏的数值模拟分析方法 - Google Patents
一种对bif铁矿顺层边坡变形破坏的数值模拟分析方法 Download PDFInfo
- Publication number
- CN117371068A CN117371068A CN202311484963.9A CN202311484963A CN117371068A CN 117371068 A CN117371068 A CN 117371068A CN 202311484963 A CN202311484963 A CN 202311484963A CN 117371068 A CN117371068 A CN 117371068A
- Authority
- CN
- China
- Prior art keywords
- slope
- side slope
- numerical simulation
- layer
- state
- 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
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 36
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 title claims abstract description 32
- 238000004458 analytical method Methods 0.000 title claims abstract description 28
- 229910052742 iron Inorganic materials 0.000 title claims abstract description 16
- 239000011435 rock Substances 0.000 claims abstract description 62
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 46
- 238000000034 method Methods 0.000 claims abstract description 22
- 239000011044 quartzite Substances 0.000 claims abstract description 19
- 239000007787 solid Substances 0.000 claims abstract description 13
- SZVJSHCCFOBDDC-UHFFFAOYSA-N iron(II,III) oxide Inorganic materials O=[Fe]O[Fe]O[Fe]=O SZVJSHCCFOBDDC-UHFFFAOYSA-N 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000006073 displacement reaction Methods 0.000 claims description 36
- 238000009826 distribution Methods 0.000 claims description 24
- 230000009471 action Effects 0.000 claims description 22
- 229920006395 saturated elastomer Polymers 0.000 claims description 20
- 239000010438 granite Substances 0.000 claims description 19
- 238000012360 testing method Methods 0.000 claims description 13
- 239000010453 quartz Substances 0.000 claims description 9
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 230000009467 reduction Effects 0.000 claims description 6
- 238000012669 compression test Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000007670 refining Methods 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 11
- 239000003673 groundwater Substances 0.000 abstract description 2
- 239000010410 layer Substances 0.000 description 41
- 238000005065 mining Methods 0.000 description 15
- 238000010586 diagram Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 11
- 238000004519 manufacturing process Methods 0.000 description 8
- 230000002829 reductive effect Effects 0.000 description 8
- 238000010008 shearing Methods 0.000 description 8
- 239000002689 soil Substances 0.000 description 7
- 239000002734 clay mineral Substances 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- KWYUFKZDYYNOTN-UHFFFAOYSA-M Potassium hydroxide Chemical compound [OH-].[K+] KWYUFKZDYYNOTN-UHFFFAOYSA-M 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000009933 burial Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000002301 combined effect Effects 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- YGANSGVIUGARFR-UHFFFAOYSA-N dipotassium dioxosilane oxo(oxoalumanyloxy)alumane oxygen(2-) Chemical compound [O--].[K+].[K+].O=[Si]=O.O=[Al]O[Al]=O YGANSGVIUGARFR-UHFFFAOYSA-N 0.000 description 2
- 239000010433 feldspar Substances 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000010445 mica Substances 0.000 description 2
- 229910052618 mica group Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 235000010755 mineral Nutrition 0.000 description 2
- 229910052627 muscovite Inorganic materials 0.000 description 2
- 229940072033 potash Drugs 0.000 description 2
- BWHMMNNQKKPAPP-UHFFFAOYSA-L potassium carbonate Substances [K+].[K+].[O-]C([O-])=O BWHMMNNQKKPAPP-UHFFFAOYSA-L 0.000 description 2
- 235000015320 potassium carbonate Nutrition 0.000 description 2
- 238000005482 strain hardening Methods 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- XOJVVFBFDXDTEG-UHFFFAOYSA-N Norphytane Natural products CC(C)CCCC(C)CCCC(C)CCCC(C)C XOJVVFBFDXDTEG-UHFFFAOYSA-N 0.000 description 1
- 206010035148 Plague Diseases 0.000 description 1
- 241000923606 Schistes Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- DLHONNLASJQAHX-UHFFFAOYSA-N aluminum;potassium;oxygen(2-);silicon(4+) Chemical compound [O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[Al+3].[Si+4].[Si+4].[Si+4].[K+] DLHONNLASJQAHX-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052626 biotite Inorganic materials 0.000 description 1
- 229910001919 chlorite Inorganic materials 0.000 description 1
- 229910052619 chlorite group Inorganic materials 0.000 description 1
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 239000004567 concrete Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Optimization (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Mathematical Analysis (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明公开了一种对BIF铁矿顺层边坡变形破坏的数值模拟分析方法，采用Rhino三维建模软件建立边坡实体模型，Griddle细化网格模型，并采用FLAC3D有限差分软件进行数值模拟计算，以司家营铁矿东帮N26勘探线‑126m运输平台边坡为研究背景，分析了岩性、水及风化作用等条件影响下边坡的变形破坏规律。结合司家营铁矿区的工程地质条件，定性分析了磁铁石英岩及黑云变粒岩岩石特性、地下水及风化作用等因素对边坡变形破坏的影响。根据有限差分软件FLAC3D数值模拟，定量分析不同因素对边坡岩体变形的影响。本发明的方法对进一步丰富边坡稳定性研究、预防边坡滑坡破坏具有一定的指导意义。
Description
技术领域
本发明涉及矿山安全生产技术领域，尤其涉及一种对BIF铁矿顺层边坡变形破坏的数值模拟分析方法。
背景技术
矿山边坡作为露天矿山采场的主要构成之一，其稳定性问题是矿山安全生产的关键技术问题，同时也是保证矿山正常生产的重要条件。随着露天矿山规模的逐渐扩大，其开挖深度和面积也在不断增加，这些对边坡的稳定性无疑是一个巨大的考验。因此，为了矿区的正常运营，针对边坡的稳定性问题应该尽早、尽快的深入分析研究，提前发现存在的问题和潜在的安全隐患，尽早制定相应的防范对策。
冀东司家营BIF铁矿矿区采场内曾多次发生边坡滑坡、失稳等问题。2012年，司家营铁矿东北部边坡体受风化作用的影响发生严重的滑坡破坏，其岩体主要为角闪岩，受风化作用影响后的角闪岩在长期卸荷作用下，岩体沿着不同方向的解理开裂崩塌成块，造成边坡失稳。此外，司家营铁矿区采场边坡分别在2013年、2016年及2017年也都发生过不同类型和不同程度的失稳破坏，严重威胁着矿区边坡下平台、人员及设备安全，并对矿区正常生产造成了巨大的影响。同时，在处理边坡的滑坡破坏中也耗费了大量的人力与财力。所以，进行矿山边坡稳定性的科学评价对采矿作业的安全开展有着至关重要的意义，直接关系到矿山的安全，也关系着采矿的效益。
司家营铁矿区所处大地构造位置为燕山褶皱带山海关隆起之昌黎凸起的西南边缘地带，褶皱构造发育，边坡岩体位于司家营向斜的东翼、新河背斜的西翼，区内发育少量褶皱和小型弯曲。采区内中元古界长城系大红峪组石英砂岩作为沉积盖层呈角度不整合覆盖在新太古界单塔子岩群白庙子组黑云变粒岩、混合岩化黑云变粒岩、钾长石化云母片岩等变质岩系及BIF(太古代沉积变质铁矿)之上。采场上盘边坡(西帮)属于反倾边坡，上部由第四系粉土、黏土和粉细砂组成(正逐步剥离)，下部为石英砂岩。岩层倾向近西，倾角40°～55°。采场下盘边坡(东帮)属于顺层岩质边坡，第四系和石英砂岩已经全部被剥离，因此以钾长石化白云母片岩、黑云变粒岩夹BIF为主，边坡总体产状为倾向270°，倾角50°～60°。其中位于下部层位的钾长石化白云母片岩经风化后呈土黄色，极易破碎，遇水后呈泥状，是易造成边坡失稳滑塌的岩石类型。根据现有地质钻孔资料，初步推测强风化层底部埋深约100m，中风化层底部埋深约140m，以下为弱风化层。
目前针对司家营矿区围岩存在的工程问题，主要为对复杂地质条件影响下的各类岩性及其工程特征缺乏有效的深入了解。同时，以往对司家营的围岩稳定性分析主要在结构面分析及数值模拟方面，对司家营矿区围岩的岩性方面研究薄弱。本发明通过分析研究磁铁石英岩、黑云变粒岩等围岩的地质工程特性对其稳定性进行评价，可以为矿区的安全生产提供理论依据和科学指导，研究成果具有一定的理论价值和实际意义。
FLAC3D有限差分软件是一款功能强大的数值模拟软件，其在数值模拟分析中地质模型的后处理方面较为擅长，而对于数值模拟分析前期的地质模型建立，其内置建模板块如Extrusion建模等在内的模型建立功能则相对比较薄弱，只能建立一些较为简单的假三维模型，而假三维模型在进行边坡稳定性的分析时约束比较多，并且也不能真实的反映地层等因素对于边坡稳定性的影响。为了最大程度上还原边坡的实际地层特征，真实的反映在水等不同因素影响下的边坡变形破坏特征，本发明使用专业的三维建模软件建立真三维模型，并通过接口软件将其处理导入FLAC3D数值分析软件中，从三维的角度分析水等不同因素对顺层边坡的稳定性影响。
发明内容
本发明的目的是提供一种对BIF铁矿顺层边坡变形破坏的数值模拟分析方法，对进一步丰富边坡稳定性研究、预防边坡滑坡破坏具有一定的指导意义。
为了实现上述目的，本发明提供如下技术方案：
一种对BIF铁矿顺层边坡变形破坏的数值模拟分析方法，包括以下步骤：
S1、通过Rhino建模软件建立边坡实体模型，选取摩尔-库伦本构模型分析边坡各地层岩性工程力学性能，将边坡各地层岩性自上而下划分三个工程地质单元：磁铁石英岩、黑云变粒岩和云母石英岩，其中片麻岩为下伏基岩；通过Rhino软件画出边坡网格，再由Griddle软件的Gsurf模块细化网格模型，最终通过Griddle的Gvol模块转换生成*.f3grid格式文件并导入FLAC3D软件中运行数值模拟计算；
S2、通过FLAC3D数值模拟，定量分析天然原始状态下、饱水状态下、风化状态下以及水和风化共同作用下的顺层边坡的水平及竖向位移分布情况及塑性区的分布情况，采用强度折减法分别对天然原始状态下、饱水状态下、风化状态下以及水和风化共同作用下的边坡整体安全系数进行计算，获取边坡的潜在滑动剪切面。
进一步地，步骤S1通过Rhino建模软件建立边坡实体模型，具体的操作过程为：使用Rhino建模软件的曲线工具将所研究边坡Z坐标轴上剖面的线条轮廓画出，此时边坡的平面线条轮廓已有，再使用Rhino建模软件的曲面工具在线的基础上以线成面，此时边坡的平面轮廓已有，然后使用Rhino建模软件的实体工具在面的基础上以面挤出实体，此时边坡的实体模型呈现出来。
进一步地，步骤S2针对天然原始状态下的边坡通过FLAC3D及强度折减法进行数值模拟分析，数值模拟参数赋值依据边坡在不同含水条件及不同风化程度下岩石样品单轴抗压试验及抗剪强度试验等力学试验和岩石的物理试验，赋天然状态下第一层磁铁石英岩体的物理力学参数为：粘聚力411N、内摩擦角39°，第二层黑云变粒岩的物理力学参数为粘聚力220N、内摩擦角37°；赋饱和状态下第一层磁铁石英岩体的物理力学参数为：、粘聚力390N、内摩擦角35°，第二层黑云变粒岩的物理力学参数为粘聚力202N、内摩擦角34°；赋风化状态下第一层磁铁石英岩的物理力学参数为：粘聚力12N、内摩擦角34°，第二层黑云变粒岩的物理力学参数为粘聚力120N、内摩擦角34°；赋饱水风化状态下第一层磁铁石英岩的物理力学参数为：粘聚力8N、内摩擦角31°，第二层黑云变粒岩的物理力学参数为粘聚力106、内摩擦角33°。
与现有技术相比，本发明的有益效果为：
本发明提出的对BIF铁矿顺层边坡变形破坏的数值模拟分析方法，采用Rhino三维建模软件建立边坡实体模型，Griddle细化网格模型，并采用FLAC3D有限差分软件进行数值模拟计算，以司家营铁矿东帮N26勘探线-126m运输平台边坡为研究背景，分析了岩性、水及风化作用等条件影响下边坡的变形破坏规律。结合司家营铁矿区的工程地质条件，定性分析了磁铁石英岩及黑云变粒岩岩石特性、地下水及风化作用等因素对边坡变形破坏的影响。根据有限差分软件FLAC3D数值模拟，定量分析不同因素对边坡岩体变形的影响。天然原始状态、饱水状态、风化状态及共同作用下边坡岩体的竖向应力分布方面有着相同的分布情况，竖向应力分布从上至下不断增大，有着明显的分层及递进规律。在应力值方面，风化作用影响下边坡浅层的应力值最小。本发明的方法对进一步丰富边坡稳定性研究、预防边坡滑坡破坏具有一定的指导意义。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的边坡模型图；
图2为本发明实施例提供的天然原始状态边坡最大剪切应变增量图。
图3(a)为本发明实施例提供的饱和风化边坡水平位移分布图。
图3(b)为本发明实施例提供的饱和风化边坡竖向位移分布图。
图4为本发明实施例提供的饱和风化边坡总位移云图。
图5为本发明实施例提供的饱和风化边坡塑性区分布图。
图6为本发明实施例提供的饱和风化边坡最大剪切应变增量图。
具体实施方式
为了使本领域的技术人员更好地理解本发明的技术方案，下面将结合附图和实施例对本发明作进一步的详细介绍。
本发明以冀东司家营铁矿N26勘探线-126m运输平台边坡工程为研究对象，定量分析边坡的变形破坏特征，对其在水、风化等不同因素影响下的变形破坏规律进行数值模拟分析，评价其稳定性。
本发明基于对司家营铁矿磁铁石英岩及黑云变粒岩微观、物理、水理及力学特性的研究，可知磁铁石英岩及黑云变粒岩受风化作用内部矿物蚀变并生成绿泥石等黏土矿物，且内部孔、裂隙的发育，使得岩石表现出强度降低、软化性强等显著特征，进一步对边坡的稳定性造成不利影响，对矿区的安全生产极为不利。因此，研究受风化及水影响下边坡岩体的变形破坏特征对矿区的边坡防治及安全生产具有重要意义。
本发明提出的对BIF铁矿顺层边坡变形破坏的数值模拟分析方法，包括以下步骤：
S1、通过Rhino建模软件建立边坡实体模型，具体的操作过程为：根据实际，使用Rhino建模软件的曲线工具将所研究边坡Z坐标轴上剖面的线条轮廓画出，此时边坡的平面线条轮廓已有，再使用Rhino建模软件的曲面工具在线的基础上以线成面，此时边坡的平面轮廓已有，然后使用Rhino建模软件的实体工具在面的基础上以面挤出实体，此时边坡的实体模型呈现出来。由于在后续通过FLAC3D有限元软件进行稳定性分析时，边坡模型是以网格的形式呈现，故需要对实体模型进行网格的划分，具体操作过程为：选中研究边坡的实体模型，使用Rhino软件的网格工具初步划分出边坡实体的网格模型，此时边坡实体的网格模型较为粗略，然后使用外部嵌入在Rhino里Griddle软件的Gsurf模块对边坡网格模型进行更为细致的划分，从而生成细致网格模型；最终通过Griddle的Gvol模块转换生成*.f3grid格式文件并导入FLAC3D软件中进行数值模拟计算。
具体地，本实施例边坡的模拟选取司家营铁矿东帮N26勘探线-126m运输平台旁边坡为研究对象。根据现场揭露，边坡体坡高约16m，坡脚至右边界距离约24m，为顺倾边坡，斜坡坡度较陡，约24°～30°。边坡各地层岩性自上而下按照工程力学性能并结合工程特性共划分三个工程地质单元：磁铁石英岩、黑云变粒岩、云母石英岩，其中片麻岩为下伏基岩。
在划分网格时，根据坡体的岩层分布情况，在尽可能还原坡体真实分层的同时，又要尽可能对重点分析区域内的网格进行细致的划分，对需要重点分析的区域及地层，在其相对应的网格模型区域通过Rhino的“属性”栏设置较小的数值来达到对该区域网格进行更小更密的划分，针对原始边坡所建立的模型及分组情况如图1所示。
本构模型是对岩土材料力学性质特性的经验性描述，代表的是外载条件下岩体、土体的应力—应变关系。对于本构模型的选择，本发明考虑了工程材料的已知力学特性及本构模型的适用范围。在FLAC3D有限差分软件中，内置有12种岩土本构模型以适应各种工程分析的需要，它们是：开挖模型(空模型)；弹性模型(各向同性、横观各向同性和正交各向同性弹性模型)；塑性模型(德鲁克-普拉格模型、摩尔-库伦模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型、双屈服模型、霍克-布朗模型和修正的剑桥模型)。其中，摩尔库伦模型是边坡稳定、地下开挖等工程中最通用的岩土本构模型，它适用于那些剪应力下屈服，但剪应力只取决于最大、最小主应力，而第二主应力对屈服不产生影响的材料，代表的材料类型有：松散或胶结的材料，如土体、岩石、混凝土等，基于此，本发明模拟选取摩尔-库伦本构模型。
S2、通过FLAC3D数值模拟，定量分析天然原始状态下、饱水状态下、风化状态下和水及风化共同作用下的边坡的水平及竖向位移分布情况及塑性区的分布情况，采用强度折减法分别对天然原始状态下、饱水状态下、风化状态下和水及风化共同作用下的边坡整体安全系数进行计算，获取边坡的潜在滑动剪切面。
为了更好的体现后期风化作用及水等因素对边坡体稳定性的影响，并与其形成对比使得后续分析结果更直观、更科学，针对天然原始状态下的边坡通过FLAC3D及强度折减法进行数值模拟分析，数值模拟参数赋值依据边坡在不同含水条件及不同风化程度下岩石样品单轴抗压试验及抗剪强度试验等力学试验和岩石的物理试验，赋天然状态下第一层磁铁石英岩体的物理力学参数为：粘聚力411N、内摩擦角39°，第二层黑云变粒岩的物理力学参数为粘聚力220N、内摩擦角37°；赋饱和状态下第一层磁铁石英岩体的物理力学参数为：、粘聚力390N、内摩擦角35°，第二层黑云变粒岩的物理力学参数为粘聚力202N、内摩擦角34°；赋风化状态下第一层磁铁石英岩的物理力学参数为：粘聚力12N、内摩擦角34°，第二层黑云变粒岩的物理力学参数为粘聚力120N、内摩擦角34°；赋饱水风化状态下第一层磁铁石英岩的物理力学参数为：粘聚力8N、内摩擦角31°，第二层黑云变粒岩的物理力学参数为粘聚力106、内摩擦角33°。可以据此分析天然原始状态下边坡的最大剪切应变增量分布情况并判断其潜在的剪切滑动面，分析结果如图2所示。
通过观察天然原始状态下边坡的最大剪切应变增量分布图，可以看出，天然原始边坡的最大剪切应变增量峰值集中在上覆地层与下伏基岩接触面的位置，说明整个坡体的剪切滑动面有极大可能出现在上覆地层与基岩接触面位置，为边坡潜在的滑动带。在发明针对天然原始状态下的边坡模拟中，由FLAC3D配合其“model fos”模块对边坡整体安全系数进行计算，即通过FLAC3D的“强度折减法”来求解边坡的安全系数，从而分析边坡的稳定性。“强度折减法”的求解过程为将边坡岩土体抗剪强度参数(粘聚力、内摩擦角)不断折减降低，直至将其降低至边坡体达到破坏状态为止，并以岩土体实际抗剪强度与此时临近破坏时的折减后抗剪强度的比值作为安全系数值，当安全系数大于1时，边坡处于稳定状态，当安全系数等于1时边坡处于破坏极限平衡状态，当边坡小于1时边坡处于破坏状态。通过FLAC3D强度折减法分析天然原始状态下的边坡安全系数为4，其整体处于稳定的状态，观察FLAC3D数值模拟后的init.f3sav边坡初始状态文件里的“Zone State By Average”图、“Zone X Displacement”图、“Zone ZDisplacement”图，即区域平均状态图，X方向位移增量图、Z方向位移增量图，可看出边坡区域内部的塑性区域及相关的位移量增量，可知其无明显塑性区及位移量。
根据边坡在不同含水条件及不同风化程度下岩石样品单轴抗压试验及抗剪强度试验等力学试验和岩石的物理试验获取的边坡岩石单轴抗压强度等力学参数，将其代入到等效Mohr-Coulomb强度参数公式中求解Mohr-Coulomb强度参数粘聚力和内摩擦角等抗剪强度参数，并对风化边坡体赋予饱和状态下岩体的Mohr-Coulomb强度参数。通过数值模拟分析在重力作用下水及风化共同作用边坡体的变形破坏特征(包括边坡位移及塑性区、边坡安全系数)。其中边坡首地层(磁铁石英岩)为强风化岩层，第二地层(黑云变粒岩)未出露为弱风化岩层。
(1)边坡位移及塑性区分析
通过Rhino建立边坡网格模型，并由Griddle软件细分生成网格单元后导入FLAC3D，设置边坡模型只在X及Z方向的边界条件及处于静力状态的初始条件等条件并运用强度折减法来分析边坡的稳定性，生成结果文件后，可获得水及风化作用影响下边坡的水平及竖向位移分布情况，如图3(a)和图3(b)所示。由水及风化共同作用影响下边坡的位移分布图可以看出，相较于在饱水及风化作用影响下边坡的水平及竖向位移分布，共同作用影响下的边坡体向下方运输平台水平方向发展的区域主要集中于上覆岩层第一层，且水平方向的位移峰值集中在覆岩层第一层的中下部，且向垂直方向发展的区域主要集中在上覆岩层第一层，位移峰值集中在上覆岩层第一层的中下部，在实际工程中此位移峰值的分布特点对边坡下方的人员、设备会造成极大地威胁。共同作用下边坡在水平及垂直方向的位移量相较于边坡各自在饱水及风化状态下的位移量相比有着非常显著的增大。
受水及风化共同作用影响边坡的总位移分布情况见图4。通过观察此时边坡的总位移分布情况图可以看到，此时边坡的总位移峰值增长至59.2mm，边坡出现比较大的位移。这是由于边坡岩体表层的磁铁石英岩受到较强的风化作用导致，结合对风化磁铁石英岩微观组构的分析，风化后磁铁石英岩内部黑云母等矿物蚀变为黏土矿物，导致岩石吸水性及软化性大大增强，黏土矿物摩氏硬度较低，其遇水后强度呈指数级别下降，进一步造成风化岩石强度较大程度的降低。
受水及风化共同作用扰动下的边坡塑性区分布情况见图5。通过观察在水及风化共同作用扰动下边坡体塑性区的分布情况图，可以看出，与天然原始状态、饱水状态及受风化作用状态下边坡体的塑性区相比较，受到风化及水共同作用下的边坡体塑性区有着非常显著的变化，共同作用下边坡体的塑性区沿浅层风化带不断展布并向内部延伸，塑性区面积大大增加且深度加深，贯穿连通上覆岩层第一层岩体，这表明边坡体上覆岩层第一层已经发生剪切破坏，进而产生上覆岩层首层坡体向运输平台方向延伸的相对滑移。
(2)边坡安全系数
通过FLAC3D对水及风化作用共同影响下边坡模型的强度指标进行强度折减法计算，获取在水及风化作用共同影响下边坡的潜在滑动剪切面及安全系数，如图6所示。
天然原始状态下的边坡体由新鲜岩体组成，边坡整体安全系数较高，岩体比较稳定，没有出现明显的位移及塑性区域，剪切滑动带处于上覆岩体与基岩接触面的位置；受水影响下的边坡体，由于第一二层岩体饱和，主要塑性区在第二层与第三层地层截面处，出现了局部的小塑性区，边坡岩体整体的安全系数有明显的下降，岩体向运输平台方向(X轴负方向)有较小的位移量，但仍处于比较稳定的状态，潜在剪切滑动带没有明显变化；受风化作用影响下的边坡体，剪切滑动带上移至第一层岩体中下位置，由于岩体第一层受到较强的风化作用，岩体浅层发生剪切破坏，出现发生在岩体表面的塑性区，此时受风化影响岩体的位移量比饱水时岩体的位移量要大，然而总体位移量依然不明显，边坡整体安全系数有进一步的下降，但仍然处于比较稳定的状态；受水及风化共同作用下的边坡体，剪切滑动带位置与受风化作用边坡体相比没有明显改变，滑体内出现较大面积的塑性区，且沿滑动带向运输平台方向出现较大的位移，位移量达59.20mm，此时边坡整体安全系数下降至1.164，表明此时边坡体临近破坏状态，如不进行及时有效的防护措施，边坡体对下方运输平台人员及设备安全将造成严重威胁。
与天然原始状态、饱水状态及风化状态下的边坡模型最大剪切应变增量分布情况相对比，可以看出，水及风化作用共同作用下边坡体的最大剪切应变增量的峰值集中在上覆地层第一层的中下区域，其与风化作用影响下边坡体的最大剪切应变增量峰值的集中位置相似，这也有力的说明此时整个边坡的坡体剪切滑动面处于上覆地层第一层的中下区域。该模拟结果表明：在水及风化作用的单独影响下，边坡内部及浅层出现局部的小塑性区，且产生较小的位移，在水及风化共同作用下，边坡体浅层塑性区贯穿连通，出现较大剪切破坏塑性区，并且产生较大的位移量，根据强度折减法计算此时边坡的整体安全系数为1.164，接近1，安全系数较低，此时边坡整体处于危险状态。
综上所述，通过数值模拟过程中边坡模型的总位移云图及塑性区图可以看出，相较于天然原始状态、饱水状态及风化作用下的边坡，受水及风化共同作用下的边坡体上覆地层第一层出现了大范围且贯穿连通的由剪切破坏造成的塑性区，且在坡中沿滑动带出现向运输平台方向延伸的滑动位移，滑动位移达59.20mm，通过软件强度折减法求解的边坡整体安全系数临近1，即边坡体临近破坏状态，可知受水及风化共同作用后的边坡安全系数相较于其他三种状态下的边坡，安全系数有更进一步的降低，且边坡此时整体已经处于危险的状态，如不进行及时有效的包括疏水、排水等在内的治理措施，边坡在外部因素下干扰下极易发生崩塌、滑坡的危险。
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，但这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (3)
1.一种对BIF铁矿顺层边坡变形破坏的数值模拟分析方法，其特征在于，包括以下步骤：
S1、通过Rhino建模软件建立边坡实体模型，选取摩尔-库伦本构模型分析边坡各地层岩性工程力学性能，将边坡各地层岩性自上而下划分三个工程地质单元：磁铁石英岩、黑云变粒岩和云母石英岩，其中片麻岩为下伏基岩；通过Rhino软件画出边坡网格，再由Griddle软件的Gsurf模块细化网格模型，最终通过Griddle的Gvol模块转换生成*.f3grid格式文件并导入FLAC3D软件中运行数值模拟计算；
S2、通过FLAC3D数值模拟，定量分析天然原始状态下、饱水状态下、风化状态下以及水和风化共同作用下的顺层边坡的水平及竖向位移分布情况及塑性区的分布情况，采用强度折减法分别对天然原始状态下、饱水状态下、风化状态下以及水和风化共同作用下的边坡整体安全系数进行计算，获取边坡的潜在滑动剪切面。
2.根据权利要求1所述的对BIF铁矿顺层边坡变形破坏的数值模拟分析方法，其特征在于，步骤S1通过Rhino建模软件建立边坡实体模型，具体的操作过程为：使用Rhino建模软件的曲线工具将所研究边坡Z坐标轴上剖面的线条轮廓画出，此时边坡的平面线条轮廓已有，再使用Rhino建模软件的曲面工具在线的基础上以线成面，此时边坡的平面轮廓已有，然后使用Rhino建模软件的实体工具在面的基础上以面挤出实体，此时边坡的实体模型呈现出来。
3.根据权利要求1所述的对BIF铁矿顺层边坡变形破坏的数值模拟分析方法，其特征在于，步骤S2针对天然原始状态下的边坡通过FLAC3D及强度折减法进行数值模拟分析，数值模拟参数赋值依据边坡在不同含水条件及不同风化程度下岩石样品单轴抗压试验及抗剪强度试验等力学试验和岩石的物理试验，赋天然状态下第一层磁铁石英岩体的物理力学参数为：粘聚力411N、内摩擦角39°，第二层黑云变粒岩的物理力学参数为粘聚力220N、内摩擦角37°；赋饱和状态下第一层磁铁石英岩体的物理力学参数为：、粘聚力390N、内摩擦角35°，第二层黑云变粒岩的物理力学参数为粘聚力202N、内摩擦角34°；赋风化状态下第一层磁铁石英岩的物理力学参数为：粘聚力12N、内摩擦角34°，第二层黑云变粒岩的物理力学参数为粘聚力120N、内摩擦角34°；赋饱水风化状态下第一层磁铁石英岩的物理力学参数为：粘聚力8N、内摩擦角31°，第二层黑云变粒岩的物理力学参数为粘聚力106、内摩擦角33°。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311484963.9A CN117371068A (zh) | 2023-11-09 | 2023-11-09 | 一种对bif铁矿顺层边坡变形破坏的数值模拟分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311484963.9A CN117371068A (zh) | 2023-11-09 | 2023-11-09 | 一种对bif铁矿顺层边坡变形破坏的数值模拟分析方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117371068A true CN117371068A (zh) | 2024-01-09 |
Family
ID=89400317
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311484963.9A Pending CN117371068A (zh) | 2023-11-09 | 2023-11-09 | 一种对bif铁矿顺层边坡变形破坏的数值模拟分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117371068A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117556644A (zh) * | 2024-01-12 | 2024-02-13 | 长江勘测规划设计研究有限责任公司 | 一种陡倾顺层岩质高边坡的顺向坡定量判别方法 |
-
2023
- 2023-11-09 CN CN202311484963.9A patent/CN117371068A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117556644A (zh) * | 2024-01-12 | 2024-02-13 | 长江勘测规划设计研究有限责任公司 | 一种陡倾顺层岩质高边坡的顺向坡定量判别方法 |
| CN117556644B (zh) * | 2024-01-12 | 2024-05-10 | 长江勘测规划设计研究有限责任公司 | 一种陡倾顺层岩质高边坡的顺向坡定量判别方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xu et al. | Discrete element modeling of strata and surface movement induced by mining under open-pit final slope | |
| Liu et al. | Stability assessment of the Three-Gorges Dam foundation, China, using physical and numerical modeling—Part I: physical model tests | |
| Rybak et al. | Geomechanical substantiation of parameters of technology for mining salt deposits with a backfill | |
| Tutluoglu et al. | Two and three dimensional analysis of a slope failure in a lignite mine | |
| Xue et al. | A risk prediction method for water or mud inrush from water-bearing faults in subsea tunnel based on cusp catastrophe model | |
| CN117371068A (zh) | 一种对bif铁矿顺层边坡变形破坏的数值模拟分析方法 | |
| Li et al. | Fracture development at laminated floor layers under longwall face in deep coal mining | |
| Wei et al. | Experiment and numerical simulation of overburden and surface damage law in shallow coal seam mining under the gully | |
| CN103437318A (zh) | 悬浮式海堤的工后沉降计算预测方法 | |
| Su et al. | Stability prediction and optimal angle of high slope in open-pit mine based on two-dimension limit equilibrium method and three-dimension numerical simulation | |
| Coşar | Application of rock mass classification systems for future support design of the Dim Tunnel near Alanya | |
| Oyanguren et al. | Stability analysis of Llerin Rockfill Dam: An in situ direct shear test | |
| Yang et al. | Deformation patterns and failure mechanism of high and steep stratified rock slopes with upper steep and lower gentle style induced by step-by-step excavations | |
| Wang et al. | A new calculation method for tunneling-caused stratum settlement | |
| Huang et al. | Three-dimensional stability analysis of slurry trench based on Mohr-Coulomb nonlinear failure criterion | |
| Tai et al. | Reasonable gangue section length for disposing gangue pollutants in the new green mixed workface | |
| Meng et al. | 3D nonlinear analysis of stilling basin in complex fractured dam foundation | |
| Gou et al. | Evaluation of underground blast-induced ground motions through near-surface low-velocity geological layers | |
| Lu et al. | Monitoring and analysis of ground subsidence and backfill stress distribution in Jinchuan Mine, China | |
| Zhu et al. | Determining the minimum thickness of a crown pillar and optimizing the cross-sectional size of point pillars for the safe exploitation of a subsea gold mine | |
| Wang et al. | Numerical analysis for mining-induced stress and plastic evolution involving influencing factors: high in situ stress, excavation rate and multilayered heterogeneity | |
| Hu et al. | Physical Experiments on the Deformation of Strata with Different Properties Induced by Underground Mining. | |
| Mangushev et al. | ANALYTICAL AND NUMERICAL METHODS FOR DETERMINING THE CARRYING CAPACITY OF A PILE BARETT ON WEAK SOILS IN DEEP PITS | |
| Argimbaev et al. | Determining the safe working height of an power shovel for the development of iron-bearing tailing dumps | |
| Al Heib et al. | Analysis of Faults' Effect on the Stability of Surface Lignite Mining Areas Using the Distinct Element Method |
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 |