[发明专利]一种基于射孔约束的EnKF微地震事件位置反演方法

说明书

技术领域

本发明属于石油天然气地震勘探领域，主要应用于石油地震勘探的微地震监测，具体地说，本发明涉及一种能够更精确地确定微地震事件位置的反演方法。

背景技术

在微地震压裂中，我们可以得到每一个检波器接收到的每一个事件的P波和S波的初至旅行时。通过三分量检波器，事件的方位角可以从矢量图中得到。P波和S波旅行时包含了从事件位置到检波器之间的距离信息，比如事件发生时间。由矢量图得到的方位角给出了事件位置的方位信息。这里我们用向量d_obs来代表测量数据向量，包括每一检波器接收到的每一事件的纵横波旅行时和方位角。

在微震监测中需要确定的参数包括每个事件的位置(x，y，z)和发生时间(t₀)。假设有Nr个检波器和N个事件点，则参数向量X的维数为4×N，如X＝[x_i，y_i，z_i，t_0i]^T，其中i＝1，2，…，N。如果P波和S波的旅行时以及方位角都是可用的，那么测量数据向量最大维数为3×Nr，如d＝[P_j，S_j，A_j]^T，其中j＝1，2，…，Nr。

以前的文献中，地震定位是基于位置参数和不同接收点的初至旅行时之间是线性关系的，其表达形式如下面的式1：

A^kΔX^k+1+S^k+1＝ΔT^k…………………………………………式1

其中T是剩余旅行时，即对每一个接收器计算得到的旅行时与观测到的旅行时之间的差异，可能包括P波或者S波旅行时之一，或者两者都包括。向量X包括所有事件的位置参数(x，y，z)和事件的起始时间t₀。矩阵A是敏感性矩阵，由第k次迭代预测得到的事件位置参数计算的旅行时得到。向量S包含由错误的速度结构引起的校正参数，向量S可以用不同的算法来处理得到，比如静态位置关系法(ST)，特定震源位置关系法(SSST)。当向量S是用静态位置关系法得到时，每一个检波器接收的信号s值与事件位置无关。当向量S是用特定震源位置关系法得到时，校正关系是接收器位置与事件位置间的一个方程。在两种相关定位方法中，都至少需要用迭代方法来求解X_K+1和S_K+1中的一个。每一次迭代，事件位置参数都用下面加入了校正关系S的式2进行更新：

X_k+1＝X^k+ΔX^k+1………………………………………………式2

用更新后的新的位置向量X，我们进行正演模型得到旅行时，重新计算T。S_k+1通过求解方程1得到。

上述的校正关系可以减少由不精确的速度模型反演得到的微地震事件位置的不确定性范围。然而，其也对更多参数进行优化，因此增加了反演问题的不适定性，例如，其使最优化问题趋于拥有多解性。因此上述方法可能会使计算得到的微地震事件位置远离真实的事件位置，即异常值。

发明内容

本发明的一方面提供了一种能够减少传统定位方法的不确定性，以更精确地确定微地震事件位置的基于射孔约束的EnKF微地震事件位置反演方法。

根据本发明的一方面，一种基于射孔约束的EnKF微地震事件位置反演方法包括下述步骤：(1)建立速度场；(2)选择第k个事件点；(3)通过求解程函方程得到第j个检波器(其中，j＝1，…，Nr，Nr表示检波器的总个数)接收到的该事件点的初至旅行时；(4)应用集合卡尔曼滤波方法将步骤(3)中计算得到的初至旅行时与第j个检波器实际接收到的初至旅行时进行拟合；(5)检查每一个检波器所接收到的该事件点的旅行时是否全部拟合，如果检查结果为否，则返回到步骤(3)，并使变量j增1，以对第j+1个检波器进行拟合，如果检测结果为是，则完成事件定位。

根据本发明的一方面，可通过计算声波测井和射孔时间数据来建立速度场。

根据本发明的一方面，步骤(4)可包括下述步骤：在贝叶斯框架下定义目标函数；应用集合卡尔曼滤波来进行事件定位，首先随机生成Ne个事件位置的实现；对第j个检波器接收到的微震事件初至旅行时进行扰动，生成Ne个观测数据向量i＝1...Ne；通过拟合第j个检波器接收到的事件的初至旅行时来更新事件位置向量X；使变量j增1，重复上述过程以对第j+1个检波器接收到的事件初至旅行时进行拟合。