[發(fā)明專利]一種變偏移距VSP帶限角度積分疊前時(shí)間偏移方法和裝置

有效

• 申請?zhí)枺?/td>: 202010643464.X
• 申請日: 2020-07-07
• 公開（公告）號: CN111751875B
• 公開（公告）日: 2022-05-20
• 發(fā)明（設(shè)計(jì)）人: 李建國
• 申請（專利權(quán)）人: 中油奧博（成都）科技有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/34;G01V1/36
• 代理公司: 成都方圓聿聯(lián)專利代理事務(wù)所(普通合伙) 51241
• 代理人: 李鵬

本發(fā)明公開了一種變偏移距VSP帶限角度積分疊前時(shí)間偏移方法和裝置，所述方法包括：包括以下步驟：S1.輸入2維深度速度模型、變偏移距VSP共檢波點(diǎn)道集、變偏移距VSP共炮點(diǎn)走時(shí)表、變偏移距VSP共檢波點(diǎn)走時(shí)表、偏移參數(shù)；S2.建立炮點(diǎn)、檢波點(diǎn)在2維地質(zhì)模型中的坐標(biāo)；S3.設(shè)計(jì)帶限角度濾波器的帶限角度加權(quán)系數(shù)；S4.實(shí)現(xiàn)變偏移距VSP共檢波點(diǎn)帶限角度積分疊前時(shí)間偏移成像；S5.將得到的帶限角度積分疊前時(shí)間偏移重排成共成像道集，共成像道集疊加，得到變偏移VSP帶限角度積分疊前時(shí)間偏移疊加成像。本發(fā)明能夠直接生成變偏移距VSP帶限角度的積分疊前時(shí)間偏移成像，無需生成疊前角度道集，有效提高了數(shù)據(jù)處理效率。

技術(shù)領(lǐng)域

本發(fā)明涉及地震數(shù)據(jù)偏移成像方法，特別是涉及一種變偏移距VSP帶限角度積分疊前時(shí)間偏移方法和裝置。

背景技術(shù)

變偏移距VSP(Walkaway VSP)是指多個(gè)檢波器組成的陣列陳放在觀測井中，地表布設(shè)一條人工激發(fā)的多個(gè)炮點(diǎn)組成的炮線，是一種VSP的二維觀測方式，在炮點(diǎn)人工激發(fā)地震數(shù)據(jù)的條件下，實(shí)現(xiàn)變偏移距VSP井中數(shù)據(jù)采集，采集系統(tǒng)立體圖如圖1所示；變偏移距VSP采集，其實(shí)就是在每一個(gè)炮點(diǎn)激發(fā)時(shí)，由各個(gè)檢波點(diǎn)的檢波器進(jìn)行數(shù)據(jù)采集；變偏移距VSP可以觀測到反射縱波、反射轉(zhuǎn)換橫波，利用偏移成像研究井旁構(gòu)造、儲層預(yù)測。

王楠等研究了《表驅(qū)Kirchhoff疊前時(shí)間偏移角度域成像方法》，文獻(xiàn)利用兩點(diǎn)射線追蹤計(jì)算走時(shí)和反射張角，實(shí)現(xiàn)地面地震數(shù)據(jù)Kirchhoff疊前時(shí)間偏移角度域成像。王珺等研究了《克希霍夫法VSP多波聯(lián)合成像》，文獻(xiàn)實(shí)現(xiàn)了偏移距VSP反射P波、反射S波、透射P波、透射S波克希霍夫法聯(lián)合成像。

現(xiàn)有的積分疊前時(shí)間偏移方法多數(shù)是針對地面地震觀測方式的，變偏移距VSP的積分疊前時(shí)間偏移文獻(xiàn)很少，數(shù)據(jù)處理較為不便。

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

本發(fā)明的目的在于克服現(xiàn)有技術(shù)的不足，提供一種變偏移距VSP帶限角度積分疊前時(shí)間偏移方法和裝置，能夠直接生成變偏移距VSP帶限角度的積分疊前時(shí)間偏移成像，無需生成疊前角度道集，有效提高了數(shù)據(jù)處理效率。

本發(fā)明的目的是通過以下技術(shù)方案來實(shí)現(xiàn)的：一種變偏移距VSP帶限角度積分疊前時(shí)間偏移方法，包括以下步驟：

S1.輸入2維深度速度模型、變偏移距VSP共檢波點(diǎn)道集、變偏移距VSP共炮點(diǎn)走時(shí)表、變偏移距VSP共檢波點(diǎn)走時(shí)表、偏移參數(shù)；

S2.從輸入的數(shù)據(jù)中獲取炮點(diǎn)坐標(biāo)、炮點(diǎn)深度、檢波點(diǎn)坐標(biāo)、檢波點(diǎn)深度、井口大地坐標(biāo)、井口模型坐標(biāo)，建立炮點(diǎn)、檢波點(diǎn)在2維地質(zhì)模型中的坐標(biāo)；

S3.在步驟S2的2維地質(zhì)模型坐標(biāo)系中，采用射線追蹤計(jì)算炮點(diǎn)、檢波點(diǎn)在共散射點(diǎn)處的射線角度，并從步驟S1輸入的數(shù)據(jù)中獲取角度下限、角度上限，設(shè)計(jì)帶限角度濾波器的帶限角度加權(quán)系數(shù)；

S4.在步驟S2的2維地質(zhì)模型坐標(biāo)系中，利用步驟S1的模型計(jì)算時(shí)深表，利用步驟S1中的走時(shí)表計(jì)算共散射走時(shí)，用步驟S3得到的帶限角度加權(quán)系數(shù)孔徑內(nèi)積分求和，實(shí)現(xiàn)變偏移距VSP共檢波點(diǎn)帶限角度積分疊前時(shí)間偏移成像；

S5.將步驟S4得到的帶限角度積分疊前時(shí)間偏移重排成共成像道集，共成像道集疊加，得到變偏移VSP帶限角度積分疊前時(shí)間偏移疊加成像。

進(jìn)一步地，所述步驟S1包括：

S101.輸入變偏移距VSP共檢波點(diǎn)波場數(shù)據(jù){VSPData}，讀取炮點(diǎn)坐標(biāo){SX,SY}、炮點(diǎn)深度SZ、檢波點(diǎn)坐標(biāo){RX,RY}、檢波點(diǎn)深度RZ、井口大地坐標(biāo){WMX,WMY}；

S102.輸入2維網(wǎng)格深度速度模型{ModVel}，讀取網(wǎng)格節(jié)點(diǎn)數(shù)MM×MN、網(wǎng)格大小dg、井口模型坐標(biāo)MWMX信息；

S103.輸入變偏移距VSP共炮點(diǎn)走時(shí)表{RayTS}；