面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法技术

一种基于加权反演的数值模拟力学参数校准方法。适用于矿区地表沉陷数值模型建立、校准及地表沉陷动态预计研究使用。其步骤为：建立初始数值模型；利用遗传算法分别反演适用于非充分采动和充分采动情况下数值模型静态岩体力学参数；确定不同采动程度下岩体力学参数对评价指标的权重；根据非充分采动和充分采动情况下模型静态岩体力学参数及不同采动程度下岩体力学参数的权重，通过加权反演方法，获得面向地表沉陷动态预计的数值模型岩体力学参数。其步骤简单，兼顾煤矿开采的时空特点，克服了常规反演得到的数值模型岩体力学参数难以满足地表沉陷动态预计的需求，同时解决了常规数值模型岩体力学参数校准的问题，具有广泛的实用性。

Weighted inversion method of rock mechanics parameters for numerical model for dynamic prediction of surface subsidence

A calibration method of numerical simulation mechanical parameters based on weighted inversion is presented. It is suitable for the establishment, calibration and dynamic prediction of surface subsidence in mining area. The steps are as follows: establishing the initial numerical model; inversing the static rock mechanics parameters of the numerical model under inadequate mining and full mining respectively by genetic algorithm; determining the weight of the rock mechanics parameters to the evaluation index under different mining degrees; establishing the static rock mechanics parameters of the model under inadequate mining and full mining conditions and rock mechanics parameters under different mining degrees. By weighted inversion method, rock mechanics parameters of numerical model for dynamic prediction of surface subsidence are obtained. The procedure is simple, taking into account the temporal and spatial characteristics of coal mining. It overcomes the difficulty of the rock mechanics parameters of the conventional inversion numerical model to meet the needs of dynamic prediction of surface subsidence. At the same time, it solves the problem of calibration of rock mechanics parameters of the conventional numerical model, and has a wide range of practicability.

【技术实现步骤摘要】
面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法
本专利技术涉及一种煤矿开采数值模型岩土力学参数反演方法，尤其适用于地表沉陷动态预计的煤矿开采数值模型岩体力学参数加权反演方法。技术背景地下煤矿开采引起的地表沉陷在世界各矿区普遍存在。地表沉陷可能会导致地面建筑物损坏。为了保护地面建筑物，必须准确预测地面沉陷，使建筑物的损害程度小于其抗变形能力，从而合理设计地下开采。考虑到采矿活动是一个时空过程，地表动态移动与变形对地面建筑物的影响可能更为显著。由于采动影响下地表动态移动与变形和工作面地质采矿条件、开采方式、工作面推进速度及开采时间有关，准确预测地表动态移动与变形是一件极其困难的事。部分学者围绕地表动态移动与变形预测开展了探索性研究，提出了多种时间影响函数模型。Knothe提出了用于地表变形动态预计的Knothe函数，这也是目前应用最广泛的地表动态移动与变形预测方法。俄罗斯A.A.Baryakh等建立了考虑时间因素的动态地表沉陷预计方法；刘玉成利用数学分析方法改进了Knothe时间函数，可求出沉陷速度最大的时间；徐洪钟针对Knothe时间函数描述的不足，提出了基于Logistic增长模型的时间函数；但目前大多数地表沉陷动态预测模型处在理论研究阶段，没有较多或较好的工程实例支撑。随着计算机技术和现代力学理论的不断发展，数值模拟方法因能较为全面考虑矿区地质情况，越来越广泛的被用于地表沉陷预计。例如，Xu利用a3-Dfinitedifferencemethod预测WUTONGmine开采地表沉陷；Vyazmensky利用等效连续体和混合方法实现了对地表沉陷的预测；Guo采用数值模拟方法研究了Wongawillistrippillar的地表沉陷规律；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为非充分采动本文档来自技高网...

【技术保护点】
1.一种面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法，其特征在于步骤如下：步骤1：建立研究区域初始数值模型，包括区域地质采矿条件、工作面分布情况、地应力监测结果及室内岩体力学参数，利用FLAC3D建立要研究区域的初始模型；步骤2：利用遗传算法和初始数值模型，反演得到适用于非充分采动情况下地表沉陷预计的静态岩体力学参数；步骤3：利用不同采动程度下的地表移动与变形实测数据，并根据不同岩体力学参数设置反演过程中各力学参数的相应权重，从而得到面向地表沉陷动态预计的数值模型岩体力学参数反演方法；步骤4：结合不同采动程度下岩体力学参数的权重，利用加权反演方法，获取面向地表沉陷动态预计的数值模型岩体力学参数，根据得到的不同采动程度下岩体力学参数值及相应的权重，利用上述计算公式可以计算得到面向地表沉陷动态预计的全局数值模型岩体力学参数。

【技术特征摘要】
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所述的面向地表沉陷动态预计的数值模型岩体力学参数...

【专利技术属性】
技术研发人员：李怀展，郭广礼，查剑锋，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：江苏,32