[發(fā)明專利]各向異性縱橫波走時點乘定位方法及系統(tǒng)

有效

• 申請?zhí)枺?/td>: 201710941754.0
• 申請日: 2017-10-11
• 公開（公告）號: CN109655919B
• 公開（公告）日: 2020-04-07
• 發(fā)明（設計）人: 余波
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/50
• 分類號: G01V1/50
• 代理公司: 北京思創(chuàng)畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

一種各向異性縱橫波走時點乘定位方法及系統(tǒng)。該方法可以包括：將各向異性縱波走時目標方程、各向異性橫波走時目標方程、各向異性縱橫波時差目標方程點乘，建立目標方程；通過目標方程建立射孔目標方程，通過射孔目標方程對各向異性參數(shù)進行反演，獲得射孔位置所在地層各向異性參數(shù)；建立微地震各向異性定位目標方程，對微地震事件空間位置進行定位反演，獲取微地震事件的初始徑向距離、初始深度距離；通過目標方程對微地震事件的徑向距離、深度距離、各向異性參數(shù)同時進行反演，進而獲得微地震事件的最終定位結果及修正后的各向異性參數(shù)。本發(fā)明可以實現(xiàn)各向異性介質(zhì)微地震定位處理獲得更高精度事件的徑向、深度分布。

技術領域

本發(fā)明涉及井中微地震信號處理領域，更具體地，涉及一種各向異性縱橫波走時點乘定位方法及系統(tǒng)。

背景技術

井中微地震監(jiān)測是微地震觀測方式之一，特點是井下三分量檢波器接收微地震全波場信號，相對于地面微地震監(jiān)測，井中接收到的數(shù)據(jù)信噪比較高、微地震事件個數(shù)與類型較豐富。但是，由于井中微地震檢波器個數(shù)有限(一般12-32級三分量井中檢波器)，導致監(jiān)測范圍較小，用常規(guī)井中微地震定位方法容易出現(xiàn)不穩(wěn)定、精度不高等微地震定位結果。為了解決這一難點，需開發(fā)定位精度較高的新方法。

目前，井中微地震定位方法主要有兩種思路：一是基于縱波、橫波事件旅行時正演，代表算法有網(wǎng)絡搜索法、模擬退火法、geiger法等，優(yōu)點是容易實現(xiàn)，缺點是由于初至相位信號弱導致縱波、橫波事件旅行時難以準確拾取，影響定位結果；第二種定位思路是基于波動方程褶積，代表算法有干涉法、逆時偏移法、被動源成像法，優(yōu)點是不需要拾取事件初至，缺點是對資料信噪比、速度模型要求高、檢波器個數(shù)要求較多數(shù)量，計算成本高。

以上方法中，基于走時定位方法簡單、實用且基于各向同性均勻介質(zhì)假設，適合井中微地震事件定位處理。但是在非常規(guī)致密砂巖氣、頁巖氣藏儲層壓裂微地震開發(fā)中，地層存在非均質(zhì)性，屬于各向異性介質(zhì)，微地震縱橫波走時及傳播路徑不同于各向同性，用目前定位方法無法滿足定位精度要求。因此，有必要開發(fā)一種各向異性縱橫波走時點乘定位方法及系統(tǒng)。

公開于本發(fā)明背景技術部分的信息僅僅旨在加深對本發(fā)明的一般背景技術的理解，而不應當被視為承認或以任何形式暗示該信息構成已為本領域技術人員所公知的現(xiàn)有技術。

發(fā)明內(nèi)容

本發(fā)明提出了一種各向異性縱橫波走時點乘定位方法及系統(tǒng)，可以實現(xiàn)各向異性介質(zhì)微地震定位處理獲得更高精度的微地震事件的徑向、深度分布。

根據(jù)本發(fā)明的一方面，提出了一種各向異性縱橫波走時點乘定位方法。所述方法可以包括：將各向異性縱波走時目標方程OPJ_P、各向異性橫波走時目標方程OPJ_S對各向異性縱橫波走時目標方程OPJ_PS點乘，建立目標方程OPJ；拾取射孔的縱橫波走時T_shoot,P、T_shoot,S，根據(jù)射孔位置、聲波測井縱橫波速度，通過所述目標方程建立射孔目標方程，通過所述射孔目標方程對各向異性參數(shù)進行反演，獲得各向異性參數(shù)優(yōu)選值；拾取微地震事件縱橫波走時T_event,P、T_event,S，根據(jù)所述聲波測井縱橫波速度、所述各向異性參數(shù)優(yōu)選值，通過所述目標方程建立微地震目標方程，通過所述微地震目標方程對微地震事件進行定位反演，獲取所述微地震事件的徑向距離、深度距離；根據(jù)所述微地震事件的徑向距離、深度距離與所述各向異性參數(shù)優(yōu)選值，通過所述目標方程對所述微地震事件的徑向距離、深度距離進行反演，進而獲得所述微地震事件的最終定位結果及對應的修正后的各向異性參數(shù)。

優(yōu)選地，所述目標方程OPJ為：

其中，T_P、T_S分別為拾取縱橫波實際走時，T_Pi、T_Si分別為反演縱橫波走時。