[发明专利]油藏改造微地震事件散点聚类分析方法及系统

权利要求书

1.油藏改造微地震事件散点聚类分析方法，其特征在于，该油藏改造微地震事件散点聚类分析方法包括：

步骤1、实时获取油藏改造微地震事件散点；

步骤2、计算所有微地震事件的误差椭圆；

步骤3、根据微地震事件的矩震级进行加权运算；

步骤4、计算误差椭圆范围内微地震事件的几何中心；

步骤5、通过不断迭代使得目标函数最小化，进行油藏改造人工裂缝的精确描述和精细刻画；

在步骤2，获取微地震散点后计算误差椭圆，微地震事件散点的误差椭圆计算公式如下：

其中：σ值为残时差方程r_i＝t_i-t₀-T_i(x₀,y₀,z₀)偏微分矩阵的方差，t₀是微地震事件发生的初始时刻，(x₀,y₀,z₀)是微地震事件发生的空间位置，i是检波器的序号，i＝1,2,…,N，N是检波器的总数量，t_i是第i个检波器的观测微地震事件到时，T_i是第i个检波器的理论模拟到时；B值由残时差方程r_i＝t_i-t₀-T_i(x₀,y₀,z₀)偏微分矩阵计算得到；是置信度为α、自由度为β的χ²分布，α是置信度参数，β是自由度参数，即包含独立变量的个数，χ²是卡方分布；

误差椭圆计算公式中B值是由残时差方程r_i＝t_i-t₀-T_i(x₀,y₀,z₀)的偏微分矩阵A来定义，运算公式如下：

B＝A^TA

其中：是由残时差方程r_i＝t_i-t₀-T_i(x₀,y₀,z₀)的偏微分矩阵计算得到，是理论模拟到时函数T的偏微分；A^T为偏微分 矩阵A的转置；

误差椭圆计算公式中σ值由残时差方程r_i＝t_i-t₀-T_i(x₀,y₀,z₀)的偏微分矩阵A的方差来定义，公式如下：

其中：r_i是残时差，i＝1,2,…,N；N是检波器的总数量。

2.根据权利要求1所述的油藏改造微地震事件散点聚类分析方法，其特征在于，在步骤1，微地震事件是通过布设在井中或地面的多个检波器记录波形定位得到的，包括微地震事件的发生时刻、空间位置和震源参数信息。

3.根据权利要求1所述的油藏改造微地震事件散点聚类分析方法，其特征在于，在步骤2，计算所有微地震事件散点的误差椭圆，将其作为散点聚类分析的约束条件，即作为微地震事件可以移动的几何范围。

4.根据权利要求1所述的油藏改造微地震事件散点聚类分析方法，其特征在于，在步骤3，根据微地震事件的矩震级进行加权运算，权重是通过微地震事件的矩震级计算得到的，具体计算公式如下：

M_w＝a·log₁₀(M₀)+b

其中：a是常数因子，取值2/3；b是常数因子，取值-6；M₀是地震矩，由地震波振幅低频成分的大小来决定，它反映了震源点处破裂能量的大小，断层面积越大，激发的长周期地震波的能量越大、周期越长。