CN107942380A - 一种考虑垮落带采空区的计算机数值模拟方法 - Google Patents
一种考虑垮落带采空区的计算机数值模拟方法 Download PDFInfo
- Publication number
- CN107942380A CN107942380A CN201711008161.5A CN201711008161A CN107942380A CN 107942380 A CN107942380 A CN 107942380A CN 201711008161 A CN201711008161 A CN 201711008161A CN 107942380 A CN107942380 A CN 107942380A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- goaf
- mfrac
- rock
- 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 44
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000011435 rock Substances 0.000 claims abstract description 75
- 239000000463 material Substances 0.000 claims abstract description 31
- 239000003245 coal Substances 0.000 claims abstract description 17
- 230000006835 compression Effects 0.000 claims abstract description 12
- 238000007906 compression Methods 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 230000033001 locomotion Effects 0.000 claims description 7
- 238000010008 shearing Methods 0.000 claims description 7
- 238000002474 experimental method Methods 0.000 claims description 6
- 238000012669 compression test Methods 0.000 claims description 4
- 230000006378 damage Effects 0.000 claims description 4
- 238000005553 drilling Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000004880 explosion Methods 0.000 claims description 3
- 238000009530 blood pressure measurement Methods 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 230000004044 response Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 abstract description 4
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 238000005065 mining Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000012512 characterization method Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000009533 lab test Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011236 particulate material Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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/05—Geographic models
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Environmental & Geological Engineering (AREA)
- Computer Hardware Design (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Geophysics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geology (AREA)
- Acoustics & Sound (AREA)
- Computer Graphics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于煤炭采出后采空区受上覆岩层压缩响应的计算机数值模拟技术领域，具体是一种考虑垮落带采空区的计算机数值模拟方法。包括以下步骤：第一步：建模；第二步：确定除采空区以外的岩体力学参数；第三步：将采空区材料赋予双屈服本构模型；第四步：确定采空区材料相关参数；第五步：将所有岩体力学和采空区参数代入模型进行运算，便可得到考虑采空区的数值模拟结果。本发明将垮落角考虑入数值模拟，更加符合实际。采用双屈服模型模拟采空区,其参数基于迭代试错匹配与Salamon应力‑应变曲线拟合确定，实验室试验数据与计算机模拟数据的匹配确保模拟结果的可靠性。
Description
技术领域
本发明属于煤炭采出后采空区受上覆岩层压缩响应的计算机数值模拟技术领域，具体是一种考虑垮落带采空区的计算机数值模拟方法。
背景技术
煤炭从地下采出引起上覆岩层移动和应力重新分布。长壁开采为我国主要的采煤方法，开采后形成采空区，采空区是采场煤炭采出的产物，对地形地貌影响巨大，采空区响应对覆岩运移影响显著。
岩层移动的关键层理论，通过试验、理论研究和现场证据证明了长壁开采覆岩形成三带，即垮落带、裂隙带和弯曲下沉带。图1示为长壁采场覆岩运动规律。垮落带岩体呈散体和碎裂结构，岩块间空隙大、连通性好；裂隙带岩层断裂失去连续性，层状特性仍存在；弯曲下沉带只弯曲不断裂，保持连续性和原始特征。为垮落带采空区发育特征提供了基础依据。
许多地面工程需要确定采空区响应，除了地上工程，许多地下工程位于采空区或其附近，同样涉及采空区响应，如分层开采上分层采空区下布置下分层巷道、沿空留巷、沿空掘巷以及错层位开采技术等。
采空区除了影响地形地貌外，还对采动应力分布具有决定性的重要作用，鉴于采空区的重要性，国内学者对采空区材料及其响应进行了研究。
大量对采空区破碎矸石压实特性的实验研究表明，采空区矸石压缩应力-应变曲线呈指数或近似指数的增长趋势。其中，作为该类试验的权威代表，Pappas等通过在不同煤矿不同长壁工作面取矸石样品进行多组实验室压缩试验，对目前提出的多个矸石压缩应力-应变公式分析，得出Salamon总结出的矸石压缩的应力-应变曲线最符合矸石的压实特性。
然而，采矿科学在煤矿设计中常需要进行相关数值模拟以研究设计的合理性和科学性。如，数值模拟技术对巷道应力的分布、煤柱的设计如今已应用非常广泛和流行。科学的数值模拟方法将大幅提升模拟的准确性，从而提升模拟结论的现实指导性。尽管目前垮落带采空区的实验室试验已较多，获得较多的数据和结论，但目前流行的长壁开采数值模拟均对采空区不做任何处理，及在模型中挖空后便进行运算，采空区的支撑作用没有任何体现，即覆岩压力没有传递至采空区，也没有通过采空区传递至底板，模拟结果必然与实际大相径庭。垮落带采空区的数值模拟技术的，没有一种科学的模拟方法，限制了计算机技术在采矿工程领域的指导作用。目前计算机的运算速度和容量以迅雷不及掩耳之势发展，计算机更新和换代令人应接不暇，利用好计算机，提出科学合理的采空区数值模拟方法必能大幅提高计算机在采矿领域的应用价值和指导意义。
发明内容
本发明为了解决上述问题，提供一种考虑垮落带采空区的计算机数值模拟方法。
本发明采取以下技术方案：一种考虑垮落带采空区的计算机数值模拟方法，所述方法包括以下步骤：
第一步：FLAC3D建模；
第二步：确定除采空区以外的岩体力学参数；
第三步：将采空区材料赋予双屈服本构模型；
第四步：确定采空区材料相关参数；
第五步：将所有岩体力学和采空区参数代入FLAC3D模型中进行运算，便可得到考虑采空区的数值模拟结果。
在所述的第一步中，建模尺寸数据由钻孔地质资料和工作面尺寸确定，其中，FLAC3D模型中煤层岩层的高度及岩性与钻孔获得的岩芯相对应，模型中工作面尺寸按照实际中工作面的尺寸划定；通过地应力测量、模型埋深以及质地构造，确定模型的边界条件，边界条件包括模型四周的固定情况、测压系数(水平应力与竖直应力的比值)以及顶部加载情况。处理方式为：顶部边界自由，底面限制竖直方向的运动，前后左右四个面限制水平运动，顶部可自由运动；顶部加载根据上方未建模型的高度确定，若未建高度为h，覆岩容重为γ，则顶部施加竖直向下的力为γh；侧压系数则根据实测水平应力与竖直应力的比值在模型中赋值。
在所述的第二步中，对岩样进行实验室岩石力学试验，获得RocData软件需要的实验室数据来估算σci,mi.和a并得到广义H-B强度参数mb,s和a，并得到最适M-C准则等价摩擦角和粘聚力c，这些参数用于模型中各岩层的赋值和运算，选取除采空区以外的岩体材料的本构模型，确定各岩层的物理力学参数，本构模型表达式为：
式中，σ1为破坏时的最大主应力；σ3为破坏时的最小主应力或是三轴试验中的围压；mb为岩体的H-B常数；s为与岩体特性有关的材料常数，反映岩体破碎程度，其取值范围0～1；a为表征节理岩体的常数；σci为完整岩石的单轴抗压强度UCS；而：
式中，D为节理岩体弱化因子，主要考虑爆破破坏和应力松弛对节理岩体的扰动程度，取值为0～1；GSI为地质强度指标，GSI的值由工程岩体的岩体结构、结构面特征因素综合确定；
H-B与Mohr-Coulomb准则，下文简写为H-C准则，拟合确定最适M-C准则等价摩擦角和粘聚力c的表达式为：
式中，σ3n＝σ3max/σci，σ3max为最小主应力上限值，
抗压强度和抗拉强度分别为:
σc＝σcisa
要得到岩体的广义H-B强度参数mb,s和a,需要各岩层的以下参数:完整岩石的单轴抗压强度σci；完整岩石参数mi；地质强度指标GSI；扰动系数D，通过实验室数据拟合或通过软件内置典型岩层参数范围表估算σci,mi.和a并得到广义H-B强度参数mb,s和a以及最适M-C准则等价摩擦角和粘聚力c，其中σci、mi、GSI和D参数可通过实验室数据拟合，缺少实验室数据的情况下可根据软件内置典型岩石经验参数得到σci,mi和a，该内置参数通过世界范围内的大量岩石工程收集所得。
在所述的第三步中，采空区材料赋予双屈服本构模型，该本构面模型中应力-应变呈类似指数式增长的方式。该本构模型适于模拟压应力引起体积永久缩减的低胶结粒状散体岩土材料。模型中考虑了由于各向同性压力引起的永久性体积变化，还考虑了FLAC程序中应变软化/硬化模型中剪切破坏和拉破坏判断准则以及体积屈服面。
在所述的第四步中，双屈服模型中所需的参数有：体积模量、剪切模量、剪胀角和摩擦角参数。
参数通过在FLAC3D中进行采空区材料压缩试验与Salamon经验公式迭代匹配得到；
采空区矸石应力-应变曲线符合Salamon经验公式：
式中，σ为加载于材料上的轴向压力，[MPa]；ε为在该压力下的应变，m/m；E0为初始切线模量[MPa]；εm为最大可能应变；且：
其中，σc为垮落区矸石单轴抗压强度，εm只与初始碎胀系数b有关：
初始碎胀系数为：
其中，Hc为煤层上方垮落带高度；m为煤厚。
在所述的第四步中，在FLAC3D中建立简单块体模型，底部固支，顶部不约束，其它四个面约束横向位移，通过在其表面施加竖向速率模拟加载过程，通过迭代试错法来确定体积模量、剪切模量、剪胀角和摩擦角参数，经过反复试验，可最终确定参数，使参数满足数值模拟所得的应力-应变曲线与salamon公式曲线的匹配效果理想。
与现有技术相比，本发明提供了一种考虑垮落带采空区的计算机数值模拟方法。该方法为了改善目前数值模拟中不考虑采空区造成模拟结果不符合实际情况的现状。不考虑采空区，采空区不发挥任何支撑作用，即覆岩压力没有传递至采空区，也没有通过采空区传递至底板，则会造成支承压力偏大，覆岩和底板破环和运移范围偏大，本发明与现广泛应用的技术相比具有以下优点：模型单元随机生成，而非常规矩形，避免人为设定网格造成的影响；煤岩体参数则基于H-B强度准则和M-C准则拟合得出，考虑了岩体结构的弱面，模拟数据更加可靠；将垮落角考虑入数值模拟，更加符合实际。采用双屈服模型模拟采空区,其参数基于迭代试错匹配与Salamon应力-应变曲线拟合确定，实验室试验数据与计算机模拟数据的匹配确保模拟结果的可靠性。
附图说明
图1为覆岩三带发育示意图；
图2为FLAC3D数值模型示例；
图3为RocData软件H-B准则与M-C准则换算实例；
图4为双屈服本构模型表征的力学响应图像；
图5为相似模拟所得垮落角示意图；
图6为数值模拟与Salamon经验公式曲线匹配结果。
具体实施方式
一种考虑垮落带采空区的计算机数值模拟方法，包括以下步骤：
第一步：建模。建模对计算的影响显著，建议采用外接建模接口软件，如ansys等进行建模，可建立出FLAC3D自身无法建立的随机分布的网格模型，计算更加可靠。根据钻孔数据、工作面尺寸等逐层建立数值模型。通过地应力测试，确定模型的边界条件。图2所示为一FLAC3D模型。
第二步：选取除采空区以外的岩体材料的本构模型，确定各岩层的物理力学参数。准确的本构模型是得出可靠模拟结果的基础。Hoek-Brown经验强度准则(下文简称H-B准则)是由Hoek E和Brown ET，在大量矿山岩体力学测试、室内实验基础并根据大量工程实践经验得出，该本构模型考虑了岩体的不连续性，已在世界范围内得到广泛应用和实践检验。该本构模型表达式为：
式中，σ1为破坏时的最大主应力；σ3为破坏时的最小主应力(或是三轴试验中的围压)；mb为岩体的H-B常数；s为与岩体特性有关的材料常数，反映岩体破碎程度，其取值范围0～1；a为表征节理岩体的常数；σci为完整岩石的单轴抗压强度(UCS)。而：
式中，D为节理岩体弱化因子，主要考虑爆破破坏和应力松弛对节理岩体的扰动程度，取值为0～1；GSI为地质强度指标，GSI的值由工程岩体的岩体结构、结构面特征等因素综合确定。
H-B与Mohr-Coulomb准则(下文简写为H-C准则)拟合确定最适M-C准则等价摩擦角和粘聚力c的表达式为：
式中，σ3n＝σ3max/σci，σ3max为最小主应力上限值，与特定的工程类型有关。
抗压强度和抗拉强度分别为:
σc＝σcisa
要得到岩体的广义H-B强度参数(mb,s和a),需要各岩层的以下参数:完整岩石的单轴抗压强度σci；完整岩石参数mi；地质强度指标GSI；扰动系数D。RocData能够通过实验室数据拟合或通过软件内置典型岩层参数范围表估算σci,mi.和a并得到广义H-B强度参数mb,s和a以及最适M-C准则等价摩擦角和粘聚力c。其中σci、mi、GSI和D参数可通过实验室数据拟合，缺少实验室数据的情况下可根据软件内置典型岩石经验参数得到σci,mi和a，该内置参数通过世界范围内的大量岩石工程收集所得。这些参数将用于模型中各岩层的赋值和运算。图3所示为一换算过程例子。
第三步：确定采空区材料的本构模型。合理处理采空区才能得到合理的模拟结论。采空区矸石在逐渐压实过程中强度和模量逐渐增大，双屈服本构模型适用于模拟压应力引起体积永久缩减的低胶结粒状散体岩土材料，符合采空区材料的力学行为。除了将剪切和拉伸破坏包络线纳入FLAC外，通过将体积屈服面即应变软化/硬化模型(或称为“帽盖”)纳入FLAC3D，双屈服模型将各向同性的压力引起的永久体积改变考虑在内。与采空区具有类似性质的材料在压缩时除了发生剪切屈服外，还会发生一种和剪切屈服明显相互作用的屈服模式，例如液压回填材料或者由混凝土胶结形成的小颗粒材料，而双屈服本构模型能够表征这种行为。采空区材料赋予双屈服本构模型的目的是使采空区材料满足在逐渐压实过程中强度和模量呈Salamon经验公式的指数式增长，符合试验现象，与实际接近。图4所示为双屈服破坏模型的力学响应。
第四步：确定采空区材料相关参数。双屈服模型中所需的参数有：体积模量、剪切模量、剪胀角和摩擦角参数。这些参数则通过在FLAC3D中进行采空区材料压缩试验与Salamon经验公式迭代匹配得到。
采空区矸石应力-应变曲线符合Salamon经验公式：
式中，σ为加载于材料上的轴向压力，[MPa]；ε为在该压力下的应变，m/m；E0为初始切线模量[MPa]；εm为最大可能应变。且：
其中，σc为垮落区矸石单轴抗压强度。εm只与初始碎胀系数b有关：
初始碎胀系数为：
其中，Hc为煤层上方垮落带高度，；m为煤厚。
由此可见，确定采空材料参数，需要确定垮落带采空区范围。有若干种确定采空区范围的方式，包括经验公式、相似模拟实验、穿层钻孔、地面钻孔或其它地球物理勘探手段，以及上述两种或几中方法的结合。图5所示为相似模拟所得垮落角示意图。注意，相似模拟得出的垮落带范围的准确性不如现场实测得出的可靠，视工程精度要求适当选择采空区范围测定方法。
确定采空区范围(包括垮落带高度、长度、宽度以及垮落角)后，便可通过以上公式计算出Salamon经验公式的参数E0和εm，Salamon经验公式便可确定。
在FLAC3D中建立简单块体模型，底部固支，顶部不约束，其它四个面约束横向位移。通过在其表面施加竖向速率模拟加载过程。通过迭代试错法来确定体积模量、剪切模量、剪胀角和摩擦角参数。经过反复试验，可最终确定参数，使参数满足数值模拟所得的应力-应变曲线与salamon公式曲线的匹配效果理想。图6所示为匹配结果。
第五步：将所有岩体力学和采空区参数代入模型进行运算，得到考虑采空区的数值模拟结果。
Claims (7)
1.一种考虑垮落带采空区的计算机数值模拟方法，其特征在于，所述方法包括以下步骤：
第一步：FLAC3D建模；
第二步：确定除采空区以外的岩体力学参数；
第三步：将采空区材料赋予双屈服本构模型；
第四步：确定采空区材料相关参数；
第五步：将所有岩体力学和采空区参数代入FLAC3D模型中进行运算，便可得到考虑采空区的数值模拟结果。
2.根据权利要求1所述的考虑垮落带采空区的计算机数值模拟方法，其特征在于：在所述的第一步中，建模尺寸数据由钻孔地质资料和工作面尺寸确定，其中，FLAC3D模型中煤层岩层的高度及岩性与钻孔获得的岩芯相对应，模型中工作面尺寸按照实际中工作面的尺寸划定；通过地应力测量、模型埋深以及质地构造，确定模型的边界条件，边界条件包括模型四周的固定情况、测压系数即水平应力与竖直应力的比值以及顶部加载情况，处理方式为：顶部边界自由，底面限制竖直方向的运动，前后左右四个面限制水平运动，顶部可自由运动；顶部加载根据上方未建模型的高度确定，若未建高度为h，覆岩容重为γ，则顶部施加竖直向下的力为γh；侧压系数则根据实测水平应力与竖直应力的比值在模型中赋值。
3.根据权利要求1或2所述的考虑垮落带采空区的计算机数值模拟方法，其特征在于：在所述的第二步中，过岩样进行实验室岩石力学试验，获得RocData软件需要的实验室数据来估算σci,mi.和a并得到广义H-B强度参数mb,s和a，并得到最适M-C准则等价摩擦角和粘聚力c，这些参数用于模型中各岩层的赋值和运算，选取除采空区以外的岩体材料的本构模型，确定各岩层的物理力学参数，本构模型表达式为：
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mi>b</mi>
</msub>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>a</mi>
</msup>
</mrow>
式中，σ1为破坏时的最大主应力；σ3为破坏时的最小主应力或是三轴试验中的围压；mb为岩体的H-B常数；s为与岩体特性有关的材料常数，反映岩体破碎程度，其取值范围0～1；a为表征节理岩体的常数；σci为完整岩石的单轴抗压强度UCS；而：
<mrow>
<msub>
<mi>m</mi>
<mi>b</mi>
</msub>
<mo>=</mo>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>G</mi>
<mi>S</mi>
<mi>I</mi>
<mo>-</mo>
<mn>100</mn>
</mrow>
<mrow>
<mn>28</mn>
<mo>-</mo>
<mn>14</mn>
<mi>D</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>G</mi>
<mi>S</mi>
<mi>I</mi>
<mo>-</mo>
<mn>100</mn>
</mrow>
<mrow>
<mn>9</mn>
<mo>-</mo>
<mn>3</mn>
<mi>D</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>a</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>G</mi>
<mi>S</mi>
<mi>I</mi>
<mo>/</mo>
<mn>15</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mn>20</mn>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
式中，D为节理岩体弱化因子，主要考虑爆破破坏和应力松弛对节理岩体的扰动程度，取值为0～1；GSI为地质强度指标，GSI的值由工程岩体的岩体结构、结构面特征因素综合确定；
H-B与Mohr-Coulomb准则，下文简写为H-C准则，拟合确定最适M-C准则等价摩擦角和粘聚力c的表达式为：
<mrow>
<mi>c</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>m</mi>
<mi>b</mi>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mn>3</mn>
<mi>n</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>m</mi>
<mi>b</mi>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mn>3</mn>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>a</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>a</mi>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mo>+</mo>
<mi>a</mi>
<mo>)</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>am</mi>
<mi>b</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>m</mi>
<mi>b</mi>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mn>3</mn>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>a</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>+</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
式中，σ3n＝σ3max/σci，σ3max为最小主应力上限值，
抗压强度和抗拉强度分别为:
σc＝σcisa
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>s&sigma;</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>b</mi>
</msub>
</mfrac>
</mrow>
要得到岩体的广义H-B强度参数mb,s和a,需要各岩层的以下参数:完整岩石的单轴抗压强度σci；完整岩石参数mi；地质强度指标GSI；扰动系数D，通过实验室数据拟合或通过软件内置典型岩层参数范围表估算σci,mi.和a并得到广义H-B强度参数mb,s和a以及最适M-C准则等价摩擦角和粘聚力c，其中σci、mi、GSI和D参数可通过实验室数据拟合，缺少实验室数据的情况下可根据软件内置典型岩石经验参数得到σci,mi和a，该内置参数通过世界范围内的大量岩石工程收集所得。
4.根据权利要求3所述的考虑垮落带采空区的计算机数值模拟方法，其特征在于：在所述的第三步中，采空区材料赋予双屈服本构模型。
5.根据权利要求4所述的考虑垮落带采空区的计算机数值模拟方法，其特征在于：在所述的第四步中，双屈服模型中所需的参数有：体积模量、剪切模量、剪胀角和摩擦角参数。
6.根据权利要求5所述的考虑垮落带采空区的计算机数值模拟方法，其特征在于：在所述的第四步中，参数通过在FLAC3D中进行采空区材料压缩试验与Salamon经验公式迭代匹配得到；
采空区矸石应力-应变曲线符合Salamon经验公式：
<mrow>
<mi>&sigma;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>E</mi>
<mn>0</mn>
</msub>
<mi>&epsiv;</mi>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&epsiv;</mi>
<mo>/</mo>
<msub>
<mi>&epsiv;</mi>
<mi>m</mi>
</msub>
</mrow>
</mfrac>
</mrow>
式中，σ为加载于材料上的轴向压力，[MPa]；ε为在该压力下的应变，m/m；E0为初始切线模量[MPa]；εm为最大可能应变；且：
<mrow>
<msub>
<mi>E</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>10.39</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>c</mi>
<mn>1.042</mn>
</msubsup>
</mrow>
<msup>
<mi>b</mi>
<mn>7.7</mn>
</msup>
</mfrac>
</mrow>
其中，σc为垮落区矸石单轴抗压强度，εm只与初始碎胀系数b有关：
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>b</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>b</mi>
</mfrac>
</mrow>
初始碎胀系数为：
<mrow>
<mi>b</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>H</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<mi>m</mi>
</mrow>
<msub>
<mi>H</mi>
<mi>c</mi>
</msub>
</mfrac>
</mrow>
其中，Hc为煤层上方垮落带高度；m为煤厚。
7.根据权利要求6所述的考虑垮落带采空区的计算机数值模拟方法，其特征在于：在所述的第四步中，在FLAC3D中建立简单块体模型，底部固支，顶部不约束，其它四个面约束横向位移，通过在其表面施加竖向速率模拟加载过程，通过迭代试错法来确定体积模量、剪切模量、剪胀角和摩擦角参数，经过反复试验，可最终确定参数，使参数满足数值模拟所得的应力-应变曲线与salamon公式曲线的匹配效果理想。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711008161.5A CN107942380A (zh) | 2017-10-23 | 2017-10-23 | 一种考虑垮落带采空区的计算机数值模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711008161.5A CN107942380A (zh) | 2017-10-23 | 2017-10-23 | 一种考虑垮落带采空区的计算机数值模拟方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107942380A true CN107942380A (zh) | 2018-04-20 |
Family
ID=61936498
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711008161.5A Pending CN107942380A (zh) | 2017-10-23 | 2017-10-23 | 一种考虑垮落带采空区的计算机数值模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107942380A (zh) |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108763725A (zh) * | 2018-05-24 | 2018-11-06 | 西安科技大学 | 基于采空区压实效应的工作面开采全过程数值模拟方法 |
| CN108827233A (zh) * | 2018-09-17 | 2018-11-16 | 中国地质大学（北京） | 一种两层采空区地面沉降的预测方法 |
| CN108959799A (zh) * | 2018-07-18 | 2018-12-07 | 南京维思科汽车科技有限公司 | 一种汽车空调冷凝器的建模方法 |
| CN108986208A (zh) * | 2018-07-11 | 2018-12-11 | 辽宁工程技术大学 | 一种煤矿采空区冒落形态的重构方法 |
| CN109443943A (zh) * | 2018-11-26 | 2019-03-08 | 长江大学 | 岩石样本强度准则材料参数的构建方法及设备 |
| CN109812276A (zh) * | 2019-01-22 | 2019-05-28 | 河南理工大学 | 一种基于动载应力场确定邻空巷道合理位置的方法 |
| CN110032801A (zh) * | 2019-04-15 | 2019-07-19 | 太原理工大学 | 一种含有t形割理网络节理煤体模型的重构方法 |
| CN110044719A (zh) * | 2019-05-30 | 2019-07-23 | 辽宁工程技术大学 | 一种散体矸石真三轴压缩试验装置及方法 |
| CN110489826A (zh) * | 2019-07-30 | 2019-11-22 | 石家庄铁道大学 | 岩体边坡安全系数计算方法、装置及终端设备 |
| CN110609335A (zh) * | 2019-09-17 | 2019-12-24 | 太原理工大学 | 一种基于多手段的残采区复杂条件探测方法 |
| CN110727026A (zh) * | 2019-10-10 | 2020-01-24 | 太原理工大学 | 一种基于地震-海啸连续作用的非弹性位移比谱模型 |
| CN110823729A (zh) * | 2019-11-15 | 2020-02-21 | 中国地质环境监测院 | 一种库区涉水厚层危岩体溃屈失稳的评价方法 |
| CN111472777A (zh) * | 2020-04-14 | 2020-07-31 | 西安科技大学 | 一种边坡下近水平煤层工作面煤柱宽度确定方法 |
| CN112765853A (zh) * | 2021-01-19 | 2021-05-07 | 中国科学院武汉岩土力学研究所 | 一种爆破振动可靠度设计方法 |
| CN113236365A (zh) * | 2021-05-10 | 2021-08-10 | 晋能控股煤业集团同忻煤矿山西有限公司 | 一种特厚煤层顶煤跨落角现场实测方法 |
| RU2774053C1 (ru) * | 2021-08-30 | 2022-06-14 | Аньхуэй Юниверсити оф Сайенс энд Текнолоджи | Многопольная система мониторинга и анализа для тестирования интеллектуального многомерного имитационного моделирования нагрузки |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101806054A (zh) * | 2010-04-02 | 2010-08-18 | 中铁二局股份有限公司 | 高速公路下伏房柱式采空区稳定性评价方法 |
| CN103700141A (zh) * | 2013-11-19 | 2014-04-02 | 中国矿业大学(北京) | 煤矿底板突水通道的多级尺度地质建模分析方法 |
| CN103984807A (zh) * | 2014-04-29 | 2014-08-13 | 天津大学 | 耦合精细地质信息和监控信息的大坝灌浆数值模拟方法 |
| CN104181611A (zh) * | 2014-08-28 | 2014-12-03 | 山东科技大学 | 一种矿井工作面顶底板采动破坏裂隙发育动态监测方法 |
| CN104331533A (zh) * | 2014-09-18 | 2015-02-04 | 山东科技大学 | 一种封闭不良钻孔采动响应数值模拟方法 |
| CN104504461A (zh) * | 2014-12-09 | 2015-04-08 | 鞍钢集团矿业公司 | 一种露天转地下开采诱发矿坑边坡变形破坏预测方法 |
| CN105388265A (zh) * | 2015-10-29 | 2016-03-09 | 重庆大学 | 采空区瓦斯富集及运移的三维模拟方法 |
| CN106285782A (zh) * | 2016-08-30 | 2017-01-04 | 中国矿业大学(北京) | 一种复杂地质赋存环境下的冲击地压预警方法及系统 |
| CN106884677A (zh) * | 2017-04-10 | 2017-06-23 | 大同煤矿集团有限责任公司 | 坚硬顶板特厚煤层开采强矿压预测预控方法 |
-
2017
- 2017-10-23 CN CN201711008161.5A patent/CN107942380A/zh active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101806054A (zh) * | 2010-04-02 | 2010-08-18 | 中铁二局股份有限公司 | 高速公路下伏房柱式采空区稳定性评价方法 |
| CN103700141A (zh) * | 2013-11-19 | 2014-04-02 | 中国矿业大学(北京) | 煤矿底板突水通道的多级尺度地质建模分析方法 |
| CN103984807A (zh) * | 2014-04-29 | 2014-08-13 | 天津大学 | 耦合精细地质信息和监控信息的大坝灌浆数值模拟方法 |
| CN104181611A (zh) * | 2014-08-28 | 2014-12-03 | 山东科技大学 | 一种矿井工作面顶底板采动破坏裂隙发育动态监测方法 |
| CN104331533A (zh) * | 2014-09-18 | 2015-02-04 | 山东科技大学 | 一种封闭不良钻孔采动响应数值模拟方法 |
| CN104504461A (zh) * | 2014-12-09 | 2015-04-08 | 鞍钢集团矿业公司 | 一种露天转地下开采诱发矿坑边坡变形破坏预测方法 |
| CN105388265A (zh) * | 2015-10-29 | 2016-03-09 | 重庆大学 | 采空区瓦斯富集及运移的三维模拟方法 |
| CN106285782A (zh) * | 2016-08-30 | 2017-01-04 | 中国矿业大学(北京) | 一种复杂地质赋存环境下的冲击地压预警方法及系统 |
| CN106884677A (zh) * | 2017-04-10 | 2017-06-23 | 大同煤矿集团有限责任公司 | 坚硬顶板特厚煤层开采强矿压预测预控方法 |
Non-Patent Citations (3)
| Title |
|---|
| SALAMON M D G: ""Mechanism of caving in longwall coal mining"", 《ROCK MECHANICS CONTRIBUTIONS AND CHALLENGES: PROCEEDINGS OF THE 31ST US SYMPOSIUM》 * |
| 王朋飞 等: ""非充分采动采空区域煤岩柱（体）耦合作用机制及应用"", 《岩石力学与工程学报》 * |
| 赵毅鑫 等: ""基于Hoek-Brown参数确定方法的多煤层开采工作面矿压显现规律模拟研究"", 《煤炭学报》 * |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108763725A (zh) * | 2018-05-24 | 2018-11-06 | 西安科技大学 | 基于采空区压实效应的工作面开采全过程数值模拟方法 |
| CN108986208A (zh) * | 2018-07-11 | 2018-12-11 | 辽宁工程技术大学 | 一种煤矿采空区冒落形态的重构方法 |
| CN108986208B (zh) * | 2018-07-11 | 2023-04-07 | 辽宁工程技术大学 | 一种煤矿采空区冒落形态的重构方法 |
| CN108959799A (zh) * | 2018-07-18 | 2018-12-07 | 南京维思科汽车科技有限公司 | 一种汽车空调冷凝器的建模方法 |
| CN108827233A (zh) * | 2018-09-17 | 2018-11-16 | 中国地质大学（北京） | 一种两层采空区地面沉降的预测方法 |
| CN109443943A (zh) * | 2018-11-26 | 2019-03-08 | 长江大学 | 岩石样本强度准则材料参数的构建方法及设备 |
| CN109443943B (zh) * | 2018-11-26 | 2021-01-15 | 长江大学 | 岩石样本强度准则材料参数的构建方法及设备 |
| CN109812276A (zh) * | 2019-01-22 | 2019-05-28 | 河南理工大学 | 一种基于动载应力场确定邻空巷道合理位置的方法 |
| CN109812276B (zh) * | 2019-01-22 | 2020-06-02 | 河南理工大学 | 一种基于动载应力场确定邻空巷道合理位置的方法 |
| CN110032801A (zh) * | 2019-04-15 | 2019-07-19 | 太原理工大学 | 一种含有t形割理网络节理煤体模型的重构方法 |
| CN110044719B (zh) * | 2019-05-30 | 2023-11-17 | 辽宁工程技术大学 | 一种散体矸石真三轴压缩试验装置及方法 |
| CN110044719A (zh) * | 2019-05-30 | 2019-07-23 | 辽宁工程技术大学 | 一种散体矸石真三轴压缩试验装置及方法 |
| CN110489826A (zh) * | 2019-07-30 | 2019-11-22 | 石家庄铁道大学 | 岩体边坡安全系数计算方法、装置及终端设备 |
| CN110489826B (zh) * | 2019-07-30 | 2022-12-20 | 石家庄铁道大学 | 岩体边坡安全系数计算方法、装置及终端设备 |
| CN110609335A (zh) * | 2019-09-17 | 2019-12-24 | 太原理工大学 | 一种基于多手段的残采区复杂条件探测方法 |
| CN110727026A (zh) * | 2019-10-10 | 2020-01-24 | 太原理工大学 | 一种基于地震-海啸连续作用的非弹性位移比谱模型 |
| CN110727026B (zh) * | 2019-10-10 | 2021-05-28 | 太原理工大学 | 一种基于地震-海啸连续作用的非弹性位移比谱模型 |
| CN110823729A (zh) * | 2019-11-15 | 2020-02-21 | 中国地质环境监测院 | 一种库区涉水厚层危岩体溃屈失稳的评价方法 |
| CN111472777A (zh) * | 2020-04-14 | 2020-07-31 | 西安科技大学 | 一种边坡下近水平煤层工作面煤柱宽度确定方法 |
| CN112765853A (zh) * | 2021-01-19 | 2021-05-07 | 中国科学院武汉岩土力学研究所 | 一种爆破振动可靠度设计方法 |
| CN113236365A (zh) * | 2021-05-10 | 2021-08-10 | 晋能控股煤业集团同忻煤矿山西有限公司 | 一种特厚煤层顶煤跨落角现场实测方法 |
| CN113236365B (zh) * | 2021-05-10 | 2024-02-06 | 晋能控股煤业集团同忻煤矿山西有限公司 | 一种特厚煤层顶煤跨落角现场实测方法 |
| RU2774053C1 (ru) * | 2021-08-30 | 2022-06-14 | Аньхуэй Юниверсити оф Сайенс энд Текнолоджи | Многопольная система мониторинга и анализа для тестирования интеллектуального многомерного имитационного моделирования нагрузки |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107942380A (zh) | 一种考虑垮落带采空区的计算机数值模拟方法 | |
| Liu et al. | Stability assessment of the Three-Gorges Dam foundation, China, using physical and numerical modeling—Part I: physical model tests | |
| Das et al. | Assessment of the strength of inclined coal pillars through numerical modelling based on the ubiquitous joint model | |
| US9164194B2 (en) | Method for modeling deformation in subsurface strata | |
| CN106485015B (zh) | 一种矿井断层影响范围的确定方法 | |
| Guo et al. | Simulation of three-dimensional tectonic stress fields and quantitative prediction of tectonic fracture within the Damintun Depression, Liaohe Basin, northeast China | |
| Tutluoglu et al. | Two and three dimensional analysis of a slope failure in a lignite mine | |
| Li et al. | Numerical modelling study of the load sharing law of anti-sliding piles based on the soil arching effect for Erliban landslide, China | |
| Chen et al. | Experimental and numerical study of granite blocks containing two side flaws and a tunnel-shaped opening | |
| Chen et al. | Numerical study on spalling failure of rock surrounding deep buried tunnel based on DEM | |
| Fan et al. | Damage zones induced by in situ stress unloading during excavation of diversion tunnels for the Jinping II hydropower project | |
| Shabanimashcool et al. | Discontinuous modelling of stratum cave-in in a longwall coal mine in the arctic area | |
| Zhang et al. | Seismic performance of loess-mudstone slope by centrifuge tests | |
| Zhang et al. | A physical and numerical model-based research on the subsidence features of overlying strata caused by coal mining in Henan, China | |
| Wang et al. | Stability and airtightness of a deep anhydrite cavern group used as an underground storage space: A case study | |
| CN115324560A (zh) | 利用地应力场模拟确定压裂诱发油气套管变形位置的方法 | |
| Chen et al. | A case study on the height of a water-flow fracture zone above undersea mining: Sanshandao Gold Mine, China | |
| Arasteh et al. | Discontinuous modeling of roof strata caving in a mechanized longwall mine in tabas coal mine | |
| Moss et al. | Physical, analytical, and numerical modeling of reverse‐fault displacement through near‐surface soils | |
| CN109918697A (zh) | 一种矿井冲击地压强度的预测方法及装置 | |
| Ren et al. | Geomechanical models for the quantitatively prediction of multi-scale fracture distribution in carbonate reservoirs | |
| Yan et al. | Earthquake loading response of a slope with an inclined weak intercalated layer using centrifuge modeling | |
| Isik et al. | Deformation modulus of heavily jointed–sheared and blocky greywackes by pressuremeter tests: numerical, experimental and empirical assessments | |
| Tan et al. | Influence of complicated faults on the differentiation and accumulation of in-situ stress in deep rock mass | |
| Yan et al. | A multiscale analysis of adjacent fault dislocation mechanism induced by tunnel excavation based on continuous-discrete coupling 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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180420 |