[发明专利]一种工作面覆岩断裂特征定量化分析方法

权利要求书

1.一种工作面覆岩断裂特征定量化分析方法，其特征在于，包括以下步骤：

获取煤矿的微震事件时空数据集，其中，所述微震事件时空数据集包括若干微震事件时空数据样本；

对所述微震事件时空数据集进行聚类分析，根据评价指标确定最优聚类簇个数，基于所述最优聚类簇个数获取事件簇；

获取事件簇质心的时空数据，分析在时间序列下，覆岩在工作面开采过程中的断裂特征；

获取事件簇质心的能量特征，基于所述事件簇质心的时空数据和所述能量特征，分析在空间关系下，采空区最终形成的覆岩断裂特征；

基于所述工作面开采过程中的断裂特征和所述采空区最终形成的覆岩断裂特征，对断裂线联通过程进行分析，获得覆岩断裂线扩展过程；

基于所述工作面开采过程中的断裂特征和所述采空区最终形成的覆岩断裂特征，获得垂直于工作面推进方向上的断裂线分布，结合钻孔柱状图覆岩岩性信息，得到工作面覆岩断裂线特征示意图；

基于所述工作面覆岩断裂线特征示意图和所述覆岩断裂线扩展过程，根据所述事件簇质心的时空数据，得到定量化的覆岩断裂特征。

2.根据权利要求1所述的工作面覆岩断裂特征定量化分析方法，其特征在于，所述微震事件时空数据样本的特征包括时间t和空间三维坐标x、y、z。

3.根据权利要求1所述的工作面覆岩断裂特征定量化分析方法，其特征在于，对所述微震事件数据集进行聚类分析的过程中采用k-means聚类算法。

4.根据权利要求1所述的工作面覆岩断裂特征定量化分析方法，其特征在于，所述评价指标包括：

指标1：轮廓系数平均值；

指标2：轮廓系数标准差；

指标3：大于轮廓系数平均值样本个数占总样本数比例；

指标4：计算每个簇中大于轮廓系数均值占这个簇样本数的比例，所有簇的比例的平均值；

指标5：计算每个簇中大于轮廓系数均值占这个簇样本数的比例，所有簇的比例的标准差。

5.根据权利要求4所述的工作面覆岩断裂特征定量化分析方法，其特征在于，对所述微震事件时空数据集进行聚类分析，根据评价指标确定最优聚类簇个数包括：

对5项所述评价指标进行归一化处理，基于所述微震事件时空数据集和归一化后的5项所述评价指标，构建评价指标随聚类簇个数变化曲线图，提取所述评价指标随聚类簇个数变化曲线图中5项所述评价指标的若干交汇点，比较所述交汇点，确定最优交汇点，其中，所述最优交汇点为多项所述评价指标趋于平稳的起始点，完成确定最优聚类簇个数。

6.根据权利要求5所述的工作面覆岩断裂特征定量化分析方法，其特征在于，获取所述事件簇质心的时空数据，分析在时间序列下，覆岩在工作面开采过程中的断裂特征包括：

获取所述事件簇质心的时空数据，按事件簇质心时间由前到后的顺序，将事件簇中微震事件依次绘制于三维图中，分析在时间序列下，所述覆岩在工作面开采过程中的断裂特征，其中，所述覆岩在工作面开采过程中的断裂特征为沿工作面推进方向的多个拱形分布事件簇。

7.根据权利要求6所述的工作面覆岩断裂特征定量化分析方法，其特征在于，基于所述事件簇质心的时空数据和所述能量特征，分析在空间关系下，采空区最终形成的覆岩断裂特征包括：

按所述事件簇中微震事件能量和由高到低的顺序依次绘制于三维图中，分析在空间关系下，采空区最终形成的覆岩断裂特征，其中，所述采空区最终形成的覆岩断裂特征根据破裂程度强弱分为多个层级。

8.根据权利要求7所述的工作面覆岩断裂特征定量化分析方法，其特征在于，基于所述工作面开采过程中的断裂特征和所述采空区最终形成的覆岩断裂特征，对断裂线联通过程进行分析，得到覆岩断裂线扩展过程包括：

将各所述拱形分布事件簇依次添加至所述覆岩断裂特征的不同层级分布图中，将既被标记为所述拱形分布事件簇，又被标记为当前层级事件簇的事件用不同颜色显示，获得工作面开采过程中形成的断裂与所述采空区最终形成的覆岩断裂联通情况，获得覆岩断裂线扩展过程。