基于相似性的微地震信号凸显方法
【技术领域】
[0001] 本发明涉及对地震信号进行处理的领域，更具体地讲，涉及一种基于相似性的微 地震信号凸显方法。
【背景技术】
[0002] 微地震压裂检测技术作为一项重要的非常规油气藏勘探技术，在油气藏勘探开发 (主要包括：储层压裂检测、油藏动态检测等）方面具有广泛的应用，是储层压裂检测、油藏 动态检测过程中最精确、最及时、信息最丰富的检测手段。而在利用微地震压裂检测技术进 行油气藏勘探开发时，对记录的微地震事件的信号（地震信号和噪声信号）进行处理的处 理精度，直接影响油气藏勘探开发的准确度。
[0003] 微地震事件的震源与常规地震的震源不同，其能量比较弱，约相当于几克到几十 克炸药的能量。因此，通常地面监测到的微地震事件的信噪比很低，因而，凸显记录的微地 震事件的微地震信号（即，提高信噪比）是对记录的微地震事件的信号进行处理的关键。
[0004] 目前，在对记录的微地震事件的信号进行处理时，通常通过压制噪声（即，降低噪 声波的幅度）来增大信噪比。然而，由于记录的微地震事件的信号强度很弱，在对噪声进行 压制的过程中，能量很弱的微地震信号也可能会被压制，从而导致处理后的信号的精确度 不高（即，不能准确凸显地震信号）而影响储层检测、油藏动态检测的准确度。
[0005] 因此，现有的对记录的微地震事件的信号进行处理的方法不能获得精度较高的微 地震信号。
【发明内容】
[0006] 本发明的示例性实施例在于提供一种基于相似性的微地震信号凸显方法。所述方 法能够克服现有技术中对记录的微地震事件的信号进行处理时处理精度不高的缺陷。
[0007] 根据本发明示例性实施例，提供一种基于相似性的微地震信号凸显方法，包 括：（A)确定针对预定压裂段记录的微地震事件中的主地震事件；(B)针对一个道，获取主 地震事件在所述道的第一波形，并获取除主地震事件之外的每个微地震事件在所述道的第 二波形；(C)获取第一波形与所述每个微地震事件的第二波形的互相关函数值；(D)根据获 取的互相关函数值调整所述每个微地震事件的第二波形。
[0008] 可选地，步骤（C)包括：(C1)根据第一波形的持续时间长度，将所述每个微地震事 件的第二波形划分为一个或多个子波；(C2)将第一波形的波函数与所述每个微地震事件 的第二波形的每个子波的波函数进行互相关，以获取第一波形与所述每个子波的互相关函 数值。
[0009] 可选地，步骤（C)还包括：（C3)将第一波形的波函数进行自相关，获得自相关函数 值；(C4)利用所述自相关函数值，将获取的互相关函数值进行归一化。
[0010] 可选地，步骤（D)包括：根据与每个微地震事件的第二波形的每个子波对应的归 一化的互相关函数值调整每个微地震事件的第二波形的每个子波的幅度。
[0011] 可选地，步骤（D)包括：互相关函数值越大，对与所述互相关函数值对应的子波的 调整程度越大。
[0012] 可选地，步骤（D)包括：通过下面的等式来调整所述每个微地震事件的第二波形 的每个子波的每个点的幅度：
[0013] D(tJ；i) =d(tJ；i) (l+CorJ；i)2
[0014] 其中，d(tj;i)为除主地震事件之外的微地震事件之中的第j微地震事件的第i个 子波的与时间ty对应的一个点的幅度，Cor^为与第j微地震事件的第i个子波相对应的 归一化的互相关函数值，D(ty)为所述点的调整后的幅度，其中，i和j均为大于等于1的 整数。
[0015] 可选地，步骤（A)包括：(A1)确定针对预定压裂断记录的所有微地震事件之中波 形中的地震信号与噪声的信噪比最高的预定个微地震事件；(A2)将所述预定个微地震事 件中的地震信号能量最高的微地震事件确定为主地震事件。
[0016] 可选地，第一波形为所述主地震事件的初至波的波形。
[0017] 可选地，所述方法还包括：针对所述一个道之外的其它至少一个道中的每个道执 行步骤⑶至步骤（D)。
[0018] 在根据本发明示例性实施例的基于相似性的微地震信号凸显方法中，可利用预定 压裂段内记录的各微地震事件的波形之间的相似性，通过主地震事件的波形来凸显除主地 震事件之外的微地震事件的微地震信号的强度，从而提高对记录的微地震事件的信号进行 处理的处理精度。
【附图说明】
[0019] 通过下面结合示例性地示出实施例的附图进行的描述，本发明示例性实施例的上 述和其他目的和特点将会变得更加清楚，其中：
[0020] 图1示出根据本发明示例性实施例的基于相似性的微地震信号凸显方法的流程 图。
[0021] 图2示出根据本发明示例性实施例的基于相似性的微地震信号凸显方法中获取 互相关函数值步骤的流程图。
[0022] 图3示出根据本发明另一示例性实施例的基于相似性的微地震信号凸显方法中 获取互相关函数值步骤的流程图。
【具体实施方式】
[0023] 以下，将参照附图更充分地描述本发明的示例性实施例，示例性实施例在附图中 示出。然而，可以以许多不同的形式实施示例性实施例，并且不应被解释为局限于在此阐述 的示例性实施例。相反，提供这些实施例从而本公开将会彻底和完整，并将完全地将示例性 实施例的范围传达给本领域的技术人员。
[0024] 根据本发明的示例性实施例的基于相似性的微地震信号凸显方法可以由相应的 电子设备来实施，也可通过计算机程序来实施。例如，所述方法可通过用于执行对微地震信 号进行凸显的专用设备或专用程序来执行。
[0025] 图1示出根据本发明示例性实施例的基于相似性的微地震信号凸显方法的流程 图。
[0026] 如图1所示，在步骤S100,确定针对预定压裂段记录的微地震事件中的主地震事 件。
[0027] 在任一检测点（任一道）记录的任一微地震事件的波形可用下面的等式来表示：
[0028] d (1)
[0029] 其中，d(t)表示记录的微地震事件的波函数（即，记录的微地震事件的信号）， S(t)表示微地震事件的震源函数，G(t)为格林函数，表示微地震的震源点与接收点（例如， 检波器）之间的介质结构，R(t)为接收函数，表示接收点（例如，检波器）所处位置的介质 结构，t表示时间。
[0030] 对于同一压裂段，当发生微地震时，岩石以相同的方式破裂，因此，所述预定压裂 段内发生的各微地震事件的震源函数（即，s(t))相似；而对于同一检测点，检测点与位于 同一压裂段的各震源点之间的地质结构是相似的，因此，对于所述预定压裂段内发生的各 微地震事件，G(t)和R(t)均相似，因此，根据等式（1)，所述预定压裂段内发生的各微地震 事件的微地震记录d(t)相似，S卩，波形相似。从而，可在所述预定压裂段的所有微地震事件 中选择一个微地震事件作为主地震事件，进而利用记录的主地震事件的波函数（主地震事 件的d(t))与发生在相同压裂段内的其它微地震事件（除主地震事件之外的微地震事件） 的波函数（其它微地震事件的d(t))的相似性，来对记录的其它微地震事件的波形进行处 理。
[0031] 在一实施例中，在步骤S100,可先确定针对预定压裂断记录的所有微地震事件之 中波形中的地震信号与噪声的信噪比最高的预定个微地震事件。这里，根据实际需要的不 同，所述预定个微地震事件可以是多个微地震事件，也可以是一个微地震事件。
[0032] 然后，可将所述预定个微地震事件中的地震信号能量最高的微地震事件确定为主 地震事件。这里，当确定的所述预定个微地震事件为一个微地震事件时，可直接将该微地震 事件确定为主地震事件。
[0033] 在步骤S200,针对一个道，获取主地震事件在所述道的波形（以下，称为：第一波 形），并获取除主地震事件之外的每个微地震事件在所述道的波形（以下，称为：第二波 形）。这里，根据等式（1)，同一道记录的所有微地震事件的波形（即，等式（1)中的d(t)) 具有相似性，因此，为了使用在步骤S100确定的主地震事件的波形来对其它微地震事件的 波形进行处理，需要获取针对一个道的第一波形和所述每个微地震事件的第二波形。
[0034] 在步骤S300,获取第一波形与所述每个微地震事件的第二波形的互相关函数值。
[0035] 作为示例，第一波形可为所述主地震事件的初至波的波形。这里，可将主地震事件 的初至波分别与所述每个微地震事件的第二波形进行互相关，以获取主地震事件