的初至波 與每個第二波形的互相關函數值。
[0036] 圖2示出根據本發(fā)明示例性實施例的基于相似性的微地震信號凸顯方法中獲取 互相關函數值步驟S300的流程圖。
[0037] 參照圖2,在步驟S310,根據第一波形的持續(xù)時間長度,將所述每個微地震事件的 第二波形劃分為一個或多個子波。
[0038] 這里,可根據第一波形的持續(xù)時間長度,按照時間的先后順序,將所述每個微地震 事件的第二波形劃分為一個或多個子波。作為示例,當第一波形的持續(xù)時間長度為n秒時, 可以以n秒為單位,將第二波形按照時間的先后順序劃分為一個或多個子波,這里,按時間 順序劃分的第二波形的最后一個子波的時間長度可小于n秒,n為正數。例如,當第一波形 的持續(xù)時間長度為1秒,第二波形的持續(xù)事件長度為3秒時,可將第二波形的前1秒內的波 形劃分為第1個子波,將第二波形的第1秒至第2秒內的波形劃分為第2個子波,將第二波 形的第2秒至第3秒內的波形劃分為第3個子波,S卩,按照時間的先后順序,將該第二波形 劃分為3個子波。
[0039] 在步驟S320,將第一波形的波函數與所述每個微地震事件的第二波形的每個子波 的波函數進行互相關,以獲取第一波形與所述每個子波的互相關函數值。
[0040] 這里,可將主地震事件的初至波(第一波形)的波函數與劃分的每個微地震事件 的第二波形的每個子波的波函數進行互相關,以獲取初至波與每個子波的互相關函數。這 里,由于根據等式(1)可知同一壓裂段內記錄的微地震事件的波函數具有相似性,因此,如 果劃分的子波的波函數(劃分的子波的d(t))與初至波的波函數(初至波的d(t))的相似 度越高,則說明該子波越接近地震信號,因此,當初至波與該子波進行互相關后,得到的互 相關函數值就越高,反之亦然。
[0041] 圖3示出根據本發(fā)明另一示例性實施例的基于相似性的微地震信號凸顯方法中 獲取互相關函數值步驟S300的流程圖。
[0042] 如圖3所示,除了圖2所示的步驟S310與S320之外,步驟S300還可包括:
[0043] 在步驟S330,將第一波形的波函數進行自相關,獲得自相關函數值。
[0044] 在步驟S340,利用所述自相關函數值,將獲取的互相關函數值進行歸一化。
[0045] 這里,為了后續(xù)計算的方便,可以以第一波形的自相關函數值為單位,將獲取的所 有的互相關函數值進行歸一化。
[0046] 再次參照圖1,在步驟S400,根據獲取的互相關函數值調整所述每個微地震事件 的第二波形。
[0047] 這里,可根據獲取的互相關函數值調整所述每個微地震事件的第二波形的幅度。
[0048] 作為示例,可根據與每個微地震事件的第二波形的每個子波對應的歸一化的互相 關函數值調整每個微地震事件的第二波形的每個子波的幅度。
[0049] 這里,由于當子波與主地震事件的初至波的越接近(子波越接近地震信號)時,所 述子波對應的互相關函數值越大,從而所述子波對應的歸一化的互相關函數值越大;而當 子波與主地震事件的初至波之間的差異越大(子波越接近噪聲)時,所述子波對應的互相 關函數值越小,從而所述子波對應的歸一化的互相關函數值越小,因此,可根據與每個子波 對應的根據歸一化的互相關函數值來調整每個微地震事件的第二波形的每個子波的幅度。
[0050] 作為示例,如果互相關函數值越大,則對與所述互相關函數值對應的子波的調整 程度越大。優(yōu)選地,如果互相關函數值越大,則可對與所述互相關函數值對應的子波的幅度 的凸顯程度越大。從而可對更接近地震信號的子波的幅度進行較大幅度的凸顯,而對更接 近噪聲的子波的幅度進行較小幅度的凸顯,從而可準確、有效地凸顯地震信號。
[0051] 作為示例,可通過下面的等式來調整所述每個微地震事件的第二波形的每個子波 的每個點的幅度:
[0052]D(tJ;i) =d(tJ;i) (l+CorJ;i)2 (2)
[0053] 其中,d(tj;i)為除主地震事件之外的微地震事件之中的第j微地震事件的第i個 子波的與時間對應的一個點的幅度,Coru為與第j微地震事件的第i個子波相對應的 歸一化的互相關函數值,D(ty)為所述點的調整后的幅度,其中,i和j均為大于等于1的 整數。
[0054] 通過上述步驟S100至S400,可對一個道的除主地震事件之外的每個微地震事件 的微地震信號的進行凸顯,而為了對其它道的除主地震事件之外的每個微地震事件的微地 震信號進行凸顯,可針對所述一個道之外的其它至少一個道中的每個道執(zhí)行步驟S100至 步驟S400。
[0055] 在根據本發(fā)明示例性實施例的基于相似性的微地震信號凸顯方法中,可利用預定 壓裂段內記錄的各微地震事件的波形之間的相似性,通過主地震事件的波形來凸顯除主地 震事件之外的微地震事件的微地震信號的強度,從而提高對記錄的微地震事件的信號進行 處理的處理精度。
[0056] 盡管已經參照其示例性實施例具體顯示和描述了本發(fā)明,但是本領域的技術人員 應該理解,在不脫離權利要求所限定的本發(fā)明的精神和范圍的情況下,可以對其進行形式 和細節(jié)上的各種改變。
【主權項】
1. 一種基于相似性的微地震信號凸顯方法,包括: (A) 確定針對預定壓裂段記錄的微地震事件中的主地震事件; (B) 針對一個道,獲取主地震事件在所述道的第一波形,并獲取除主地震事件之外的每 個微地震事件在所述道的第二波形; (C) 獲取第一波形與所述每個微地震事件的第二波形的互相關函數值; (D) 根據獲取的互相關函數值調整所述每個微地震事件的第二波形。
2. 如權利要求1所述的方法,其中,步驟(C)包括: (Cl)根據第一波形的持續(xù)時間長度,將所述每個微地震事件的第二波形劃分為一個或 多個子波; (C2)將第一波形的波函數與所述每個微地震事件的第二波形的每個子波的波函數進 行互相關,以獲取第一波形與所述每個子波的互相關函數值。
3. 如權利要求2所述的方法,其中,步驟(C)還包括: (C3)將第一波形的波函數進行自相關,獲得自相關函數值; (C4)利用所述自相關函數值,將獲取的互相關函數值進行歸一化。
4. 如權利要求3所述的方法,其中,步驟(D)包括:根據與每個微地震事件的第二波形 的每個子波對應的歸一化的互相關函數值調整每個微地震事件的第二波形的每個子波的 幅度。
5. 如權利要求2所述的方法,其中,步驟⑶包括:互相關函數值越大,對與所述互相 關函數值對應的子波的調整程度越大。
6. 如權利要求5所述的方法,其中,步驟(D)包括:通過下面的等式來調整所述每個微 地震事件的第二波形的每個子波的每個點的幅度: D(tJ;i) = d(tJ;i) (l+CorJ; i)2 其中,(Ktli)為除主地震事件之外的微地震事件之中的第j微地震事件的第i個子波 的與時間ty對應的一個點的幅度,Cor u為與第j微地震事件的第i個子波相對應的歸一 化的互相關函數值,D(tp)為所述點的調整后的幅度,其中,i和j均為大于等于1的整數。
7. 如權利要求1所述的方法,其中,步驟(A)包括: (Al)確定針對預定壓裂斷記錄的所有微地震事件之中波形中的地震信號與噪聲的信 噪比最高的預定個微地震事件; (A2)將所述預定個微地震事件中的地震信號能量最高的微地震事件確定為主地震事 件。
8. 如權利要求1所述的方法,其中,第一波形為所述主地震事件的初至波的波形。
9. 如權利要求1所述的方法,還包括:針對所述一個道之外的其它至少一個道中的每 個道執(zhí)行步驟(B)至步驟(D)。
【專利摘要】本發(fā)明提供一種基于相似性的微地震信號凸顯方法,包括:(A)確定針對預定壓裂段記錄的微地震事件中的主地震事件;(B)針對一個道,獲取主地震事件在所述道的第一波形,并獲取除主地震事件之外的每個微地震事件在所述道的第二波形;(C)獲取第一波形與所述每個微地震事件的第二波形的互相關函數值;(D)根據獲取的互相關函數值調整所述每個微地震事件的第二波形。
【IPC分類】G01V1-28
【公開號】CN104749627
【申請?zhí)枴緾N201510127755
【發(fā)明人】康亮, 尹陳, 劉鴻, 袁楓堯, 劉麗婷, 巫芙蓉, 曹立斌, 蔡謙, 王金鋼, 代龍
【申請人】中國石油集團川慶鉆探工程有限公司地球物理勘探公司
【公開日】2015年7月1日
【申請日】2015年3月23日