[发明专利]一种用于数值模拟的等效载荷施加方法

权利要求书

1.一种用于数值模拟的等效载荷施加方法，其特征在于，包括如下步骤：

(1)根据矿井的地质条件和开采情况确定所研究的区域，针对其煤系特征建立数值模型并完成模型网格的划分；根据实验室内测得的煤岩体物理力学性质确定模型的材料参数；

(2)根据步骤(1)中划分的模型网格，获得模型顶部网格的数量及各网格的尺寸；同时，提取模型顶部各网格中心点的坐标；

(3)采集位于模型顶部上方的相应地表的高程数据，并将所述高程数据转化为三维坐标的形式；

(4)将步骤(3)中得到的高程三维坐标导入具有插值功能的软件中进行网格数据化，得到地表网格；设置地表网格的几何参数，包括其分别在X、Y方向的最大值和最小值、网格的数量和各网格的边长、节点个数，保证所述几何参数与模型顶部网格相一致，实现与模型顶部网格节点的一一对应；

(5)将地表网格中所划分的各网格中心的z’_ij坐标减去与其相对应的模型顶部网格中心的z_ij坐标，得到模型顶部各网格与其上方地表之间的地层厚度h_ij，继而得到模型顶部各网格上覆岩层产生的自重应力σ_ij；其中，i表示第i行网格，j表示第j列网格；

(6)通过模拟软件内置的函数语言，将步骤(5)中求解的自重应力σ_ij等效施加至模型顶部各相应的网格上。

2.根据权利要求1所述的一种用于数值模拟的等效载荷施加方法，其特征在于，所述步骤(1)中，煤岩体物理力学性质包括煤岩体密度、剪切模量、体积模量、内聚力、内摩擦角和抗拉强度。

3.根据权利要求1所述的一种用于数值模拟的等效载荷施加方法，其特征在于，所述步骤(1)中，模型网格的划分方式为等尺寸网格划分，即所划分的各网格为等边长的正方形网格。

4.根据权利要求1所述的一种用于数值模拟的等效载荷施加方法，其特征在于，所述步骤(1)中，模型网格的划分方式为不等尺寸网格划分，所述的不等尺寸网格具体是指模型顶部各网格依然为等尺寸的正方形网格；除模型顶部各网格之外的模型其他部分的网格则根据距离重点研究区域的远近，其网格尺寸发生变化，即靠近重点研究区域的网格尺寸小于远离重点研究领域内的网格尺寸。

5.根据权利要求1所述的一种用于数值模拟的等效载荷施加方法，其特征在于，所述步骤(3)中，采集地表高程数据的方法有两种：一种是根据三维井上、井下对照图，利用CAD中的list命令提取所研究区域对应的地表等高线坐标；另一种是选取所研究区域附近的多个钻孔柱状图，利用图中记录的地表高程坐标，采用插值的方式获得区域内的地表高程数据。

6.根据权利要求1所述的一种用于数值模拟的等效载荷施加方法，其特征在于，所述步骤(5)中，自重应力σ_ij的计算公式为σ_ij＝ρgh_ij，式中，式中ρ为上覆岩层的平均密度，g为重力加速度。

7.根据权利要求1所述的一种用于数值模拟的等效载荷施加方法，其特征在于，数值模型的模拟方法为FLAC^3D、FLAC^2D、UDEC、3DEC、PFC有限元或者离散元分析方法中的一种。