[發明專利]一種基于射孔約束的EnKF微地震事件位置反演方法

有效

• 申請號: 201210313570.7
• 申請日: 2012-08-30
• 公開（公告）號: CN102879801A
• 公開（公告）日: 2013-01-16
• 發明（設計）人: 尹陳;劉鴻;李亞林;何光明;巫芙蓉;陳愛萍;康亮;巫駿
• 申請（專利權）人: 中國石油集團川慶鉆探工程有限公司地球物理勘探公司
• 主分類號: G01V1/00
• 分類號: G01V1/00;G01V1/28
• 代理公司: 北京銘碩知識產權代理有限公司 11286
• 代理人: 譚昌馳

1.一種基于射孔約束的EnKF微地震事件位置反演方法，包括下述步驟：

(1)建立速度場；

(2)選擇第k個事件點；

(3)通過求解程函方程得到第j個檢波器接收到的該事件點的初至旅行時其中，j＝1，...，Nr，Nr表示檢波器的總個數；

(4)應用集合卡爾曼濾波方法將步驟(3)中計算得到的初至旅行時與第j個檢波器實際接收到的初至旅行時進行擬合；

(5)檢查每一個檢波器所接收到的該事件點的旅行時是否全部擬合，如果檢查結果為否，則返回到步驟(3)，并使變量j增1，以對第j+1個檢波器進行擬合，如果檢測結果為是，則完成事件定位。

2.如權利要求1所述的方法，其中，通過計算聲波測井和射孔時間數據來建立速度場。

3.如權利要求1所述的方法，其中，步驟(4)包括下述步驟：

在貝葉斯框架下定義目標函數；

應用集合卡爾曼濾波來進行事件定位，首先隨機生成Ne個事件位置的實現；

對第j個檢波器接收到的微震事件初至旅行時進行擾動，生成Ne個觀測數據向量其中，i＝1...Ne；

通過擬合第j個檢波器接收到的事件的初至旅行時來更新事件位置向量X；

使變量j增1，重復上述過程以對第j+1個檢波器接收到的事件初至旅行時進行擬合。

4.如權利要求3所述的方法，其中，Ne個事件位置的實現為50-100個事件位置的實現。

5.如權利要求3所述的方法，其中，所述目標函數由下式表示：

O(X)=12(X-Xp‾)TCX-1(X-Xp‾)+12(dobs-dcalc)TCD-1(dobs-dcalc),]]>

其中，O(X)表示為確定微地震事件位置所建立的目標函數，X表示事件位置向量，表示上次迭代所更新的數據的平均值，在第一次計算時，為射孔位置參數，C_X是位置參數誤差協方差矩陣，CX=δx20000δy20000δz20000δt2,]]>其中是事件位置參數的標準偏差，d_obs表示觀測數據，d_calc表示正演計算得到的預測數據，C_D表示測量誤差協方差矩陣，其中是觀測數據的標準偏差，其中，k＝1，…，Nobs，Nobs為測量數據維數。