[发明专利]基于微地震事件的储层地质力学参数反演方法

权利要求书

1.基于微地震事件的储层地质力学参数反演方法，其特征在于，该基于微地震事件的储层地质力学参数反演方法包括： 

步骤1，将微地震事件云进行高斯分布函数连续化处理； 

步骤2，选择预测模型，使用测井实测结果作为初始参数，在实际施工压力下进行分时间段的预测，得到预测的微地震事件和当前的储层地质力学参数； 

步骤3，将预测的微地震事件结果进行步骤1中的高斯分布函数连续化处理，对比实际微地震事件云，获得其预测误差分布，并计算出卡尔曼滤波因子；步骤4，使用卡尔曼滤波因子对预测结果进行修正，完成储层的地质力学参更新；以及 

步骤5，按照步骤2、步骤3、步骤4对所有时间段数据进行更新，实现初始地质力学参数的更新，获得实际储层地质力学参数。 

2.根据权利要求1所述的基于微地震事件的储层地质力学参数反演方法，其特征在于，该基于微地震事件的储层地质力学参数反演方法还包括，在步骤1之前，对实际微地震数据进行震源定位处理，获得微地震事件云。 

3.根据权利要求1所述的基于微地震事件的储层地质力学参数反演方法，其特征在于，在步骤1中，将微地震事件云转换为微地震事件密度分布函数， 

并分时间段进行处理，微地震事件信息中包含了其震动时间，此震动与地下某时刻裂缝的生成相关。 

4.根据权利要求3所述的基于微地震事件的储层地质力学参数反演方法，其特征在于，在步骤1中，将整个空间中所有的已知微地震事件经过如式（1）所示的函数转换后，获得空间任意一个位置的密度分布函数值，式（1）为： 

其中，为高斯核函数，n_MEQ为某一时刻微地震事件个数，u_MEQ表示微地震事件空间坐标，u表示当前计算微地震事件密度分布的位置空间坐标，∑表示高斯函数的协方差矩阵，s(u)表示经过转化后任意位 置的微地震事件密度分布函数值。 

5.根据权利要求1所述的基于微地震事件的储层地质力学参数反演方法，其特征在于，在步骤2中，选择渗流-应力-损伤耦合模型作为预测模型。 

6.根据权利要求1所述的基于微地震事件的储层地质力学参数反演方法，其特征在于，在步骤3中，采用改进的集合卡尔曼滤波器，通过微地震事件反演出储层的地质力学参数，集合卡尔曼滤波方程采用(3)式时间更新预测方程和（4）式状态更新校正程表示： 

其中，X表示系统观测和预测的状态，Z_t-1表示已知的和时间相关的控制条件，表示预测模型参数误差的一个随机数，f_t表示从时间t-1到时间t的状态预测方程，在此处f_t表示在一定压力条件下进行正演模拟时储层参数和微地震事件的变化过程，K表示卡尔曼增益矩阵，表示状态协方差矩阵，C_d表示结果协方差矩阵，H表示测量矩阵，表示状态更新值，表示状态预测值，y^j为第j个参数模型实际观测值的扰动量，m^j表示储层的地质力学参数，，s^j表示当前获得的微地震实际密度分布函数，N_e为需要反演的地质力学参数的个数，如仅需要反演孔隙度和弹性模量，则N_e=2；为方程（3）（4）赋初值，然后通过递推的方法来求解方程在不同时刻的预测和更新值，获得储层的地质力学参数值。 

7.根据权利要求1所述的基于微地震事件的储层地质力学参数反演方法，其特征在于，该基于微地震事件的储层地质力学参数反演方法还包括，在步骤5之后，将实际储层地质力学参数输出，作为进行压裂条件下地下储层裂缝分布预测的基础数据；掌握了实际储层地质力学参数后，根据实际地质力学参数进行压裂数值模拟，预测地下裂缝发育及微地震信号的信息。 

8.根据权利要求1所述的基于微地震事件的储层地质力学参数反演方法，其特征在于，在将实际储层地质力学参数输出的步骤中，输出结果包含了在压裂模拟中不变的地质力学参数和压裂后发生变化的地质力学参数两部分。