[发明专利]微震事件时空分形维数计算方法、装置及电子设备

权利要求书

1.一种微震事件时空分形维数计算方法，其特征在于，包括：

获取微震事件的聚类分析结果；

基于所述聚类分析结果确定椭球计算总模型和不同时间段对应的椭球计算分模型；其中，时间段是基于作业工况和爆破时间进行划分的，所述椭球计算分模型的参数比例与所述椭球计算总模型的参数比例相同；所述椭球计算总模型是基于由所述聚类分析结果得到的同类微震事件的分布范围的置信区间确定的；所述椭球计算分模型是以时间为自变量，基于不同时间段内所述同类微震事件的分布范围确定的；

基于所述椭球计算总模型和所述椭球计算分模型确定微震事件时空分形维数；其中，所述微震事件时空分形维数是通过对数据点集合进行线性拟合得到的，所述数据点集合包括所述椭球计算分模型的长轴参数和相关指数，所述相关指数为不同时间段的微震数量与微震事件总数量平方的比值。

2.根据权利要求1所述的方法，其特征在于，所述基于所述聚类分析结果确定椭球计算总模型和不同时间段对应的椭球计算分模型的步骤，包括：

基于所述聚类分析结果得到同类微震事件的分布范围，并确定所述分布范围的置信区间；

基于所述分布范围的置信区间确定椭球计算总模型；

基于所述聚类分析结果得到不同时间段内所述同类微震事件的分布范围，并以时间为自变量确定不同时间段对应的椭球计算分模型。

3.根据权利要求1所述的方法，其特征在于，所述基于所述椭球计算总模型和所述椭球计算分模型确定微震事件时空分形维数的步骤，包括：

计算各个所述椭球计算分模型的相关指数；

对所述椭球计算分模型的长轴参数和所述相关指数进行对数运算得到数据点集合；

基于最小二乘法对所述数据点集合进行线性拟合，得到拟合斜率；

基于所述拟合斜率确定微震事件时空分形维数。

4.根据权利要求3所述的方法，其特征在于，所述计算各个所述椭球计算分模型的相关指数的步骤，包括：

基于所述椭球计算总模型计算微震事件总数量；

基于所述椭球计算分模型计算不同时间段的微震数量；

基于所述微震事件总数量和所述不同时间段的微震数量计算各个所述椭球计算分模型的相关指数。

5.根据权利要求4所述的方法，其特征在于，所述基于所述微震事件总数量和各个时间段的微震数量计算各个所述椭球计算分模型的相关指数的步骤，包括：

按照以下相关指数计算公式计算各个所述椭球计算分模型的相关指数：

其中，表示微震事件总数量，表示不同时间段的微震数量，n表示时间段的数量，，表示第n个时间段对应的椭球计算分模型的相关指数，表示椭球计算分模型的长轴参数。

6.根据权利要求3所述的方法，其特征在于，所述基于最小二乘法对所述数据点集合进行线性拟合，得到拟合斜率的步骤，包括：

按照以下公式对数据点集合进行线性拟合：

其中，表示椭球计算分模型的长轴参数。

7.根据权利要求6所述的方法，其特征在于，所述基于所述拟合斜率确定微震事件时空分形维数的步骤，包括：

当所述拟合斜率满足自相似性时，则将所述拟合斜率确定为微震事件时空分形维数。