CN109190136A - 面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法 - Google Patents
面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法 Download PDFInfo
- Publication number
- CN109190136A CN109190136A CN201810570435.8A CN201810570435A CN109190136A CN 109190136 A CN109190136 A CN 109190136A CN 201810570435 A CN201810570435 A CN 201810570435A CN 109190136 A CN109190136 A CN 109190136A
- Authority
- CN
- China
- Prior art keywords
- rock mass
- weight
- mechanics parameters
- subsidence
- mining
- 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
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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于加权反演的数值模拟力学参数校准方法。适用于矿区地表沉陷数值模型建立、校准及地表沉陷动态预计研究使用。其步骤为:建立初始数值模型;利用遗传算法分别反演适用于非充分采动和充分采动情况下数值模型静态岩体力学参数;确定不同采动程度下岩体力学参数对评价指标的权重;根据非充分采动和充分采动情况下模型静态岩体力学参数及不同采动程度下岩体力学参数的权重,通过加权反演方法,获得面向地表沉陷动态预计的数值模型岩体力学参数。其步骤简单,兼顾煤矿开采的时空特点,克服了常规反演得到的数值模型岩体力学参数难以满足地表沉陷动态预计的需求,同时解决了常规数值模型岩体力学参数校准的问题,具有广泛的实用性。
Description
技术领域
本发明涉及一种煤矿开采数值模型岩土力学参数反演方法,尤其适用于地表沉陷动态预计的煤矿开采数值模型岩体力学参数加权反演方法。
技术背景
地下煤矿开采引起的地表沉陷在世界各矿区普遍存在。地表沉陷可能会导致地面建筑物损坏。为了保护地面建筑物,必须准确预测地面沉陷,使建筑物的损害程度小于其抗变形能力,从而合理设计地下开采。
考虑到采矿活动是一个时空过程,地表动态移动与变形对地面建筑物的影响可能更为显著。由于采动影响下地表动态移动与变形和工作面地质采矿条件、开采方式、工作面推进速度及开采时间有关,准确预测地表动态移动与变形是一件极其困难的事。部分学者围绕地表动态移动与变形预测开展了探索性研究,提出了多种时间影响函数模型。Knothe提出了用于地表变形动态预计的Knothe函数,这也是目前应用最广泛的地表动态移动与变形预测方法。俄罗斯A.A.Baryakh等建立了考虑时间因素的动态地表沉陷预计方法;刘玉成利用数学分析方法改进了Knothe时间函数,可求出沉陷速度最大的时间;徐洪钟针对Knothe时间函数描述的不足,提出了基于Logistic增长模型的时间函数;但目前大多数地表沉陷动态预测模型处在理论研究阶段,没有较多或较好的工程实例支撑。
随着计算机技术和现代力学理论的不断发展,数值模拟方法因能较为全面考虑矿区地质情况,越来越广泛的被用于地表沉陷预计。例如,Xu利用a 3-D finite differencemethod预测WUTONG mine开采地表沉陷;Vyazmensky利用等效连续体和混合方法实现了对地表沉陷的预测;Guo采用数值模拟方法研究了Wongawilli strip pillar的地表沉陷规律;Ma采用数值模拟方法研究了不同充填率下充填开采地表沉陷规律。但上述研究主要是关于开采稳定后的地表沉陷进行预测,而对于开采过程中地表动态移动与变形预测则鲜有研究。主要原因是由于选取的岩体力学参数难以满足地表沉陷动态预计的需求。因此,目前缺少面向地表沉陷动态预计的数值模型岩体力学参数反演方法。
发明内容
针对上述技术问题,提供了一种解决了现有数值模型用于地表沉陷动态预计精度较差的难题。对于数值模型的建立与校准、地表沉陷动态预计、地面建筑物保护等均具有重要的实际意义和应用价值的面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法。
为实现上述技术目的,本发明的面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法,其特征在于步骤如下:
步骤1:建立研究区域初始数值模型,包括区域地质采矿条件、工作面分布情况、地应力监测结果及室内岩体力学参数,利用FLAC3D建立要研究区域的初始模型;
步骤2:利用遗传算法和初始数值模型,反演得到适用于非充分采动情况下地表沉陷预计的静态岩体力学参数;
步骤3:利用不同采动程度下的地表移动与变形实测数据,并根据不同岩体力学参数设置反演过程中各力学参数的相应权重,从而得到面向地表沉陷动态预计的数值模型岩体力学参数反演方法;
步骤4:结合不同采动程度下岩体力学参数的权重,利用加权反演方法,获取面向地表沉陷动态预计的数值模型岩体力学参数,根据得到的不同采动程度下岩体力学参数值及相应的权重,利用上述计算公式可以计算得到面向地表沉陷动态预计的全局数值模型岩体力学参数。
建立研究区域初始数值模型具体步骤为:
面向要进行地表沉陷动态预计的实验区域,收集研究区域地质采矿条件、工作面分布情况、地应力监测结果及室内岩体力学参数实验结果,根据获取的研究区域所有资料,利用FLAC3D建立要研究区域的初始模型,同时为了减小模型边界效应的影响,在模型边界预留1.4H的宽度,H为采深。
反演非充分采动和充分采动情况下数值模型静态岩体力学参数:
开采沉陷领域通常用下沉、移动、倾斜、水平变形及曲率来评价地表移动与变形,参数反演时候可根据实际情况将地表下沉、移动、倾斜、水平变形及曲率作为评价指标,根据研究区域非充分采动情况下实测的地表下沉、水平移动、地表倾斜、水平变形及曲率等评价指标,利用遗传算法和上述FLAC3D建立的初始数值模型,反演得到适用于非充分采动情况下地表沉陷预计的静态岩体力学参数,岩体力学参数包括弹性模量、泊松比、抗拉强度、内聚力和内摩擦角,其中弹性模量和泊松比称为岩体力学刚度参数,抗拉强度、内聚力和内摩擦角称为岩体力学强度参数,设上述得到弹性模量E1、泊松比P1、抗拉强度S1、内聚力C1和内摩擦角F1,同时根据研究区域充分采动情况下实测的地表下沉、水平移动、地表倾斜、水平变形及曲率的评价指标,利用遗传算法和上述FLACD建立的初始数值模型,反演得到适用于充分采动情况下地表沉陷动态预计的静态岩体力学参数,此时得到弹性模量E2、泊松比P2、抗拉强度S2、内聚力C2和内摩擦角F2。
确定不同采动程度下岩体力学参数对评价指标的权重具体步骤为:
开采沉陷是一个典型的时空过程,不同采动程度下模型岩体力学参数对地表下沉、水平移动、地表倾斜、水平变形及曲率评价指标的敏感性不同,可以充分利用不同采动程度下的地表移动与变形实测数据,并根据不同岩体力学参数设置反演过程中各力学参数的相应权重,从而得到面向地表沉陷动态预计的数值模型岩体力学参数反演方法,
为得到上述目的,基于上述FLAC3D建立的数值模型,就岩体力学参数对地表下沉、水平移动、倾斜、水平变形和曲率的影响进行单因素分析。所谓单因素分析就是改变弹性模量、泊松比、抗拉强度、内聚力、内摩擦角的某一参数时,其他参数取为FLAC3D建立的初始模型中的相应值。考虑到岩体力学参数之间的量纲差异较大,将FLAC3D建立的初始模型评价权重设为1,依据差分法定义相应权重,利用公式:式中,为指标ui对力学参数xj的权重;Dui(x)为指标变化值;Dxj为参数变化值;
根据上述建立的公式和单因素试验结果,可以得到不同采动程度下岩体力学参数对评价指标的权重,包括充分采动和非充分采动情况下岩体力学参数对评价指标的权重,得到非充分采动情况下弹性模量的权重为PE1、泊松比的权重为PP1、抗拉强度的权重为PS1、内聚力的权重为PC1、内摩擦角的权重为PF1;充分采动情况下弹性模量的权重为PE2、泊松比的权重为PP2、抗拉强度的权重为PS2、内聚力的权重为PC2、内摩擦角的权重为PF2。
获取面向地表沉陷动态预计的数值模型岩体力学参数具体步骤为:
反演得到了非充分采动情况下弹性模量E1、泊松比P1、抗拉强度S1、内聚力C1、内摩擦角F1以及充分采动情况下弹性模量E2、泊松比P2、抗拉强度S2、内聚力C2、内摩擦角F2,上述步骤3得到了非充分采动情况下弹性模量权重PE1、泊松比权重PP1、抗拉强度权重PS1、内聚力权重PC1、内摩擦角权重PF1以及充分采动情况下弹性模量权重PE2、泊松比权重PP2、抗拉强度权重PS2、内聚力权重PC2、内摩擦角PF2,结合不同采动程度下岩体力学参数的权重,利用加权反演方法,就可以获取面向地表沉陷动态预计的数值模型岩体力学参数,计算公式为:
X=X1×PX1+X2×PX2
根据得到的不同采动程度下岩体力学参数值及相应的权重,利用上述计算公式可以计算得到面向地表沉陷动态预计的全局数值模型岩体力学参数,式中X1为非充分采动参数取值,PX1为非充分采动参数权重,X2为充分采动参数取值,PX2为充分采动参数权重。
有益效果:本发明兼顾了开采沉陷的时空过程,创造性的利用不同采动程度下岩体力学参数对评价指标的动态敏感性,基于加权的思想,提出了面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法,解决了现在岩体力学参数反演方法应用于地表沉陷动态预计难以满足工程需求的难题,改变了数值模型只能应用于地表沉陷静态分析与预计的现状,提高了数值模型应用于地表沉陷动态预计的精度。本方法对于数值模型的建立与校准、地表沉陷动态预计、地面建筑物保护等均具有重要的实际意义和应用价值。
附图说明
图1为本发明实施的一种面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法流程图;
图2为本发明实施例得到的岩体力学参数权重与塑性区占开采影响区百分比的关系图;
图3为本发明实施例得到的不同岩体力学参数反演方法用于地表沉陷动态预计结果对比图。
具体实施方式
下面将结合图和具体实施过程对本发明做进一步详细说明:
如图1所示,本发明的面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法,包括如下步骤:
步骤1:建立研究区域初始数值模型,包括区域地质采矿条件、工作面分布情况、地应力监测结果及室内岩体力学参数,利用FLAC3D建立要研究区域的初始模型;
步骤2:利用遗传算法和初始数值模型,反演得到适用于非充分采动情况下地表沉陷预计的静态岩体力学参数;
步骤3:利用不同采动程度下的地表移动与变形实测数据,并根据不同岩体力学参数设置反演过程中各力学参数的相应权重,从而得到面向地表沉陷动态预计的数值模型岩体力学参数反演方法;
步骤4:结合不同采动程度下岩体力学参数的权重,利用加权反演方法,获取面向地表沉陷动态预计的数值模型岩体力学参数,根据得到的不同采动程度下岩体力学参数值及相应的权重,利用上述计算公式可以计算得到面向地表沉陷动态预计的全局数值模型岩体力学参数。
具体的:
建立研究区域初始数值模型具体步骤为:
面向要进行地表沉陷动态预计的实验区域,收集研究区域地质采矿条件、工作面分布情况、地应力监测结果及室内岩体力学参数实验结果,根据获取的研究区域所有资料,利用FLAC3D建立要研究区域的初始模型,同时为了减小模型边界效应的影响,在模型边界预留1.4H的宽度,H为采深。
反演非充分采动和充分采动情况下数值模型静态岩体力学参数:
开采沉陷领域通常用下沉、移动、倾斜、水平变形及曲率来评价地表移动与变形,参数反演时候可根据实际情况将地表下沉、移动、倾斜、水平变形及曲率作为评价指标,根据研究区域非充分采动情况下实测的地表下沉、水平移动、地表倾斜、水平变形及曲率等评价指标,利用遗传算法和上述FLAC3D建立的初始数值模型,反演得到适用于非充分采动情况下地表沉陷预计的静态岩体力学参数,岩体力学参数包括弹性模量、泊松比、抗拉强度、内聚力和内摩擦角,其中弹性模量和泊松比称为岩体力学刚度参数,抗拉强度、内聚力和内摩擦角称为岩体力学强度参数,设上述得到弹性模量E1、泊松比P1、抗拉强度S1、内聚力C1和内摩擦角F1,同时根据研究区域充分采动情况下实测的地表下沉、水平移动、地表倾斜、水平变形及曲率的评价指标,利用遗传算法和上述FLACD建立的初始数值模型,反演得到适用于充分采动情况下地表沉陷动态预计的静态岩体力学参数,此时得到弹性模量E2、泊松比P2、抗拉强度S2、内聚力C2和内摩擦角F2。
确定不同采动程度下岩体力学参数对评价指标的权重具体步骤为:
开采沉陷是一个典型的时空过程,不同采动程度下模型岩体力学参数对地表下沉、水平移动、地表倾斜、水平变形及曲率评价指标的敏感性不同,可以充分利用不同采动程度下的地表移动与变形实测数据,并根据不同岩体力学参数设置反演过程中各力学参数的相应权重,从而得到面向地表沉陷动态预计的数值模型岩体力学参数反演方法,
为得到上述目的,基于上述FLAC3D建立的数值模型,就岩体力学参数对地表下沉、水平移动、倾斜、水平变形和曲率的影响进行单因素分析。所谓单因素分析就是改变弹性模量、泊松比、抗拉强度、内聚力、内摩擦角的某一参数时,其他参数取为FLAC3D建立的初始模型中的相应值。考虑到岩体力学参数之间的量纲差异较大,将FLAC3D建立的初始模型评价权重设为1,依据差分法定义相应权重,利用公式:式中,为指标ui对力学参数xj的权重;Dui(x)为指标变化值;Dxj为参数变化值;
根据上述建立的公式和单因素试验结果,可以得到不同采动程度下岩体力学参数对评价指标的权重,包括充分采动和非充分采动情况下岩体力学参数对评价指标的权重,得到非充分采动情况下弹性模量的权重为PE1、泊松比的权重为PP1、抗拉强度的权重为PS1、内聚力的权重为PC1、内摩擦角的权重为PF1;充分采动情况下弹性模量的权重为PE2、泊松比的权重为PP2、抗拉强度的权重为PS2、内聚力的权重为PC2、内摩擦角的权重为PF2。
获取面向地表沉陷动态预计的数值模型岩体力学参数具体步骤为:
反演得到了非充分采动情况下弹性模量E1、泊松比P1、抗拉强度S1、内聚力C1、内摩擦角F1以及充分采动情况下弹性模量E2、泊松比P2、抗拉强度S2、内聚力C2、内摩擦角F2,上述步骤3得到了非充分采动情况下弹性模量权重PE1、泊松比权重PP1、抗拉强度权重PS1、内聚力权重PC1、内摩擦角权重PF1以及充分采动情况下弹性模量权重PE2、泊松比权重PP2、抗拉强度权重PS2、内聚力权重PC2、内摩擦角PF2,结合不同采动程度下岩体力学参数的权重,利用加权反演方法,就可以获取面向地表沉陷动态预计的数值模型岩体力学参数,所述参数共有5个,每个参数均包含非充分采动和充分采动的取值,计算公式为:
X=X1×PX1+X2×PX2
根据得到的不同采动程度下岩体力学参数值及相应的权重,利用上述计算公式可以计算得到面向地表沉陷动态预计的全局数值模型岩体力学参数,式中X1为非充分采动参数取值,PX1为非充分采动参数权重,X2为充分采动参数取值,PX2为充分采动参数权重。
具体实施例:
步骤1:建立研究区域初始数值模型
选择兖州矿区某工作面为例进行本发明方法的详细说明。该工作面煤层为近水平煤层,采厚是5m。岩层从下到上依次为粉砂岩、煤、细砂岩、粉砂岩、砾岩、细砂岩、粉砂岩、细砂岩、中砂岩及表土层。根据该工作面的地质采矿条件、地应力监测结果、工作面分布及室内岩体力学参数实验结果,利用FLAC3D软件建立了一个长1200m,宽1200m,高300m的一个数值模型。模型网格平面尺寸为30m×30m,高度随岩层高度略有不同,共划分12600单元,23064个节点。模型底部边界取u=v=0全约束边界,顶部为自由边界,左右边界为水平位移约束边界。
步骤2:反演非充分采动和充分采动情况下数值模型静态岩体力学参数
选择地表下沉和地表水平移动为评价指标。根据该工作面第6期非充分采动情况下的实测地表下沉和地表水平移动值,利用遗传算法和初始模型,反演得到该工作面非充分采动情况下数值模型弹性模量、泊松比、内聚力、抗拉强度及内摩擦角等岩体力学静态参数;同时根据该工作面第13期充分采动情况下的实测地表下沉和地表水平移动值,利用遗传算法和初始模型,反演得到该工作面充分采动情况下数值模型弹性模量、泊松比、内聚力、抗拉强度及内摩擦角等岩体力学静态参数。
步骤3:确定不同采动程度下岩体力学参数对评价指标的权重
根据上述FLAC3D建立的初始模型,利用单因素实验方法研究不同采动程度下岩体力学参数对地表下沉和水平移动的权重。考虑到岩体力学参数之间的量纲差异较大,将FLAC3D建立的初始模型评价权重设为1,依据差分法定义相应权重,计算方式如下:
式中,为指标ui对力学参数xj的权重;Dui(x)为指标变化值;Dxj为参数变化值
如图2所示,根据上述建立的公式和单因素试验结果,可以得到不同采动程度下岩体力学参数对评价指标的权重,包括充分采动和非充分采动情况下岩体力学参数对评价指标的权重。根据图2可以得到非充分采动和充分采动情况下岩体力学刚度参数和强度参数的权重。岩体力学刚度参数包括弹性模型和泊松比,岩体力学强度参数包括抗拉强度、内聚力和内摩擦角。
步骤4:获取面向地表沉陷动态预计的数值模型岩体力学参数
步骤2中已通过遗传算法反演得到了适用于非充分采动和充分采动情况下地表沉陷预计的静态岩体力学参数,再结合步骤3得到的实验区不同采动程度下岩体力学参数对评价指标的权重,利用加权的思想,就可以得到面向地表下沉动态预计的数值模型岩体力学参数。
本发明方法应用于了兖州矿区某工作面。图3为本发明实施例得到的不同岩体力学参数反演方法用于地表沉陷动态预计结果对比图。图中实测值为该工作面某点实测的地表下沉值;加权反演是指利用本发明中的方法得到的岩体力学参数进行的地表沉陷动态预计结果;常规反演是指基于第13期实测数据利用遗传算法得到的岩体力学参数进行的地表沉陷动态预计结果。6期时,常规反演和实测值相差120mm,加权反演和实测值仅相差17mm;9期时,常规反演和实测值相差184mm,加权反演和实测值相差52mm;11期时,常规反演和实测值相差89mm,加权反演和实测值相差40mm;13期时,加权反演和实测值相差65mm;加权反演和实测值相差31mm;14期时,常规反演和实测值相差265mm,加权反演和实测值相差58mm。因此,本发明能显著提高地表沉陷动态预计精度,兼顾了煤矿开采的时空特点,克服了常规反演得到的数值模型岩体力学参数难以满足地表沉陷动态预计的需求。
Claims (5)
1.一种面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法,其特征在于步骤如下:
步骤1:建立研究区域初始数值模型,包括区域地质采矿条件、工作面分布情况、地应力监测结果及室内岩体力学参数,利用FLAC3D建立要研究区域的初始模型;
步骤2:利用遗传算法和初始数值模型,反演得到适用于非充分采动情况下地表沉陷预计的静态岩体力学参数;
步骤3:利用不同采动程度下的地表移动与变形实测数据,并根据不同岩体力学参数设置反演过程中各力学参数的相应权重,从而得到面向地表沉陷动态预计的数值模型岩体力学参数反演方法;
步骤4:结合不同采动程度下岩体力学参数的权重,利用加权反演方法,获取面向地表沉陷动态预计的数值模型岩体力学参数,根据得到的不同采动程度下岩体力学参数值及相应的权重,利用上述计算公式可以计算得到面向地表沉陷动态预计的全局数值模型岩体力学参数。
2.根据权利要求1所述的面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法,其特征在于建立研究区域初始数值模型具体步骤为:
面向要进行地表沉陷动态预计的实验区域,收集研究区域地质采矿条件、工作面分布情况、地应力监测结果及室内岩体力学参数实验结果,根据获取的研究区域所有资料,利用FLAC3D建立要研究区域的初始模型,同时为了减小模型边界效应的影响,在模型边界预留1.4H的宽度,H为采深。
3.根据权利要求1所述的面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法,其特征在于反演非充分采动和充分采动情况下数值模型静态岩体力学参数:
开采沉陷领域通常用下沉、移动、倾斜、水平变形及曲率来评价地表移动与变形,参数反演时候可根据实际情况将地表下沉、移动、倾斜、水平变形及曲率作为评价指标,根据研究区域非充分采动情况下实测的地表下沉、水平移动、地表倾斜、水平变形及曲率等评价指标,利用遗传算法和上述FLAC3D建立的初始数值模型,反演得到适用于非充分采动情况下地表沉陷预计的静态岩体力学参数,岩体力学参数包括弹性模量、泊松比、抗拉强度、内聚力和内摩擦角,其中弹性模量和泊松比称为岩体力学刚度参数,抗拉强度、内聚力和内摩擦角称为岩体力学强度参数,设上述得到弹性模量E1、泊松比P1、抗拉强度S1、内聚力C1和内摩擦角F1,同时根据研究区域充分采动情况下实测的地表下沉、水平移动、地表倾斜、水平变形及曲率的评价指标,利用遗传算法和上述FLACD建立的初始数值模型,反演得到适用于充分采动情况下地表沉陷动态预计的静态岩体力学参数,此时得到弹性模量E2、泊松比P2、抗拉强度S2、内聚力C2和内摩擦角F2。
4.根据权利要求3所述的面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法,其特征在于确定不同采动程度下岩体力学参数对评价指标的权重具体步骤为:
开采沉陷是一个典型的时空过程,不同采动程度下模型岩体力学参数对地表下沉、水平移动、地表倾斜、水平变形及曲率评价指标的敏感性不同,可以充分利用不同采动程度下的地表移动与变形实测数据,并根据不同岩体力学参数设置反演过程中各力学参数的相应权重,从而得到面向地表沉陷动态预计的数值模型岩体力学参数反演方法,
为得到上述目的,基于上述FLAC3D建立的数值模型,就岩体力学参数对地表下沉、水平移动、倾斜、水平变形和曲率的影响进行单因素分析。所谓单因素分析就是改变弹性模量、泊松比、抗拉强度、内聚力、内摩擦角的某一参数时,其他参数取为FLAC3D建立的初始模型中的相应值。考虑到岩体力学参数之间的量纲差异较大,将FLAC3D建立的初始模型评价权重设为1,依据差分法定义相应权重,利用公式:式中,为指标ui对力学参数xj的权重;Dui(x)为指标变化值;Dxj为参数变化值;
根据上述建立的公式和单因素试验结果,可以得到不同采动程度下岩体力学参数对评价指标的权重,包括充分采动和非充分采动情况下岩体力学参数对评价指标的权重,得到非充分采动情况下弹性模量的权重为PE1、泊松比的权重为PP1、抗拉强度的权重为PS1、内聚力的权重为PC1、内摩擦角的权重为PF1;充分采动情况下弹性模量的权重为PE2、泊松比的权重为PP2、抗拉强度的权重为PS2、内聚力的权重为PC2、内摩擦角的权重为PF2。
5.根据权利要求3所述的面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法,其特征在于获取面向地表沉陷动态预计的数值模型岩体力学参数具体步骤为:
反演得到了非充分采动情况下弹性模量E1、泊松比P1、抗拉强度S1、内聚力C1、内摩擦角F1以及充分采动情况下弹性模量E2、泊松比P2、抗拉强度S2、内聚力C2、内摩擦角F2,上述步骤3得到了非充分采动情况下弹性模量权重PE1、泊松比权重PP1、抗拉强度权重PS1、内聚力权重PC1、内摩擦角权重PF1以及充分采动情况下弹性模量权重PE2、泊松比权重PP2、抗拉强度权重PS2、内聚力权重PC2、内摩擦角PF2,结合不同采动程度下岩体力学参数的权重,利用加权反演方法,就可以获取面向地表沉陷动态预计的数值模型岩体力学参数,计算公式为:
X=X1×PX1+X2×PX2
根据得到的不同采动程度下岩体力学参数值及相应的权重,利用上述计算公式可以计算得到面向地表沉陷动态预计的全局数值模型岩体力学参数,式中X1为非充分采动参数取值,PX1为非充分采动参数权重,X2为充分采动参数取值,PX2为充分采动参数权重。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810570435.8A CN109190136A (zh) | 2018-06-05 | 2018-06-05 | 面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810570435.8A CN109190136A (zh) | 2018-06-05 | 2018-06-05 | 面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109190136A true CN109190136A (zh) | 2019-01-11 |
Family
ID=64948626
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810570435.8A Pending CN109190136A (zh) | 2018-06-05 | 2018-06-05 | 面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109190136A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110991042A (zh) * | 2019-12-03 | 2020-04-10 | 中国矿业大学 | 面向覆岩含水层下地下气化中隔离煤柱与气化炉宽度设计方法 |
CN111737895A (zh) * | 2020-06-12 | 2020-10-02 | 鞍钢集团矿业有限公司 | 一种动态评价露天矿地下采空区顶板稳定性的方法 |
CN112184902A (zh) * | 2020-09-21 | 2021-01-05 | 东华理工大学 | 一种面向越界开采识别的地下开采面反演方法 |
CN113221228A (zh) * | 2021-06-04 | 2021-08-06 | 中国电建集团成都勘测设计研究院有限公司 | 水电站地下洞群围岩力学参数反演方法 |
CN116976751A (zh) * | 2023-08-28 | 2023-10-31 | 山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队) | 一种基于时空数据融合的煤矿采空区勘察*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339528A (zh) * | 2016-08-09 | 2017-01-18 | 鞍钢集团矿业有限公司 | 一种露天铁矿端帮地下开采诱发地表移动范围预测方法 |
-
2018
- 2018-06-05 CN CN201810570435.8A patent/CN109190136A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339528A (zh) * | 2016-08-09 | 2017-01-18 | 鞍钢集团矿业有限公司 | 一种露天铁矿端帮地下开采诱发地表移动范围预测方法 |
Non-Patent Citations (1)
Title |
---|
李怀展: "《无井式煤炭地下气化岩层移动机理与控制研究》", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑(月刊)2018年第01期》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110991042A (zh) * | 2019-12-03 | 2020-04-10 | 中国矿业大学 | 面向覆岩含水层下地下气化中隔离煤柱与气化炉宽度设计方法 |
WO2021109922A1 (zh) * | 2019-12-03 | 2021-06-10 | 中国矿业大学 | 面向覆岩含水层下地下气化中隔离煤柱与气化炉宽度设计方法 |
CN111737895A (zh) * | 2020-06-12 | 2020-10-02 | 鞍钢集团矿业有限公司 | 一种动态评价露天矿地下采空区顶板稳定性的方法 |
CN111737895B (zh) * | 2020-06-12 | 2023-11-14 | 鞍钢集团矿业有限公司 | 一种动态评价露天矿地下采空区顶板稳定性的方法 |
CN112184902A (zh) * | 2020-09-21 | 2021-01-05 | 东华理工大学 | 一种面向越界开采识别的地下开采面反演方法 |
CN113221228A (zh) * | 2021-06-04 | 2021-08-06 | 中国电建集团成都勘测设计研究院有限公司 | 水电站地下洞群围岩力学参数反演方法 |
CN116976751A (zh) * | 2023-08-28 | 2023-10-31 | 山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队) | 一种基于时空数据融合的煤矿采空区勘察*** |
CN116976751B (zh) * | 2023-08-28 | 2024-02-23 | 山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队) | 一种基于时空数据融合的煤矿采空区勘察*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109190136A (zh) | 面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法 | |
CN106326528B (zh) | 一种露天矿端帮地下开采诱发地表裂缝分布规律预测方法 | |
CN106339528B (zh) | 一种露天铁矿端帮地下开采诱发地表移动范围预测方法 | |
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
Cloetingh et al. | Neotectonics and intraplate continental topography of the northern Alpine Foreland | |
CN110135030A (zh) | 一种采空区地表沉降的预测方法 | |
CN108595878A (zh) | 一种边坡失稳机制判定方法 | |
CN108827233A (zh) | 一种两层采空区地面沉降的预测方法 | |
CN108021763A (zh) | 一种老采空区残余沉降预测方法 | |
CN111006987B (zh) | 一种储层裂缝孔渗参数预测方法 | |
Bányai et al. | Recurrent landsliding of a high bank at Dunaszekcső, Hungary: geodetic deformation monitoring and finite element modeling | |
Livio et al. | Growth of bending-moment faults due to progressive folding: Insights from sandbox models and paleoseismological implications | |
White et al. | Extensional fault-propagation folds: Mechanical models and observations from the Modoc Plateau, northeastern California | |
CN108665146B (zh) | 大型露天矿山边坡稳定性评价精度确定方法 | |
CN105242328A (zh) | 古热岩石圈厚度的确定方法及装置 | |
Dennison et al. | Stress modelling and seismicity on the Tanton fault: A case study in a South African gold mine | |
CN104820237B (zh) | 计算地层的方差体的方法 | |
Karim et al. | Volumetric estimation of coal resources in seam VI for require backfill materials of Barapukuria coal mine, Dinajpur, Bangladesh | |
Nuric et al. | Numerical modeling and computer simulation of ground movement above underground mine | |
Unruh et al. | Kinematics of postseismic relaxation from aftershock focal mechanisms of the 1994 Northridge, California, earthquake | |
CN110989034B (zh) | 一种回归-分形插值法反演测井横波时差方法 | |
CN112733242A (zh) | 一种基于物质点法的边坡大变形的确定方法 | |
Guo et al. | Study on surface subsidence law induced by deep mining of large steep metal deposit | |
CN106650036B (zh) | 水平井钻井方法和装置 | |
CN114036785B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190111 |
|
WD01 | Invention patent application deemed withdrawn after publication |