[發(fā)明專利]一種井下逆時(shí)偏移成像識(shí)別隱蔽儲(chǔ)層的方法及裝置

有效

• 申請(qǐng)?zhí)枺?/td>: 201310601138.2
• 申請(qǐng)日: 2013-11-25
• 公開（公告）號(hào): CN103758511A
• 公開（公告）日: 2014-04-30
• 發(fā)明（設(shè)計(jì)）人: 李寧;郭宏偉;王兵;武宏亮;肖承文;張宮;馮慶付
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油天然氣股份有限公司
• 主分類號(hào): E21B49/00
• 分類號(hào): E21B49/00
• 代理公司: 北京三友知識(shí)產(chǎn)權(quán)代理有限公司 11127
• 代理人: 任默聞

1.一種井下逆時(shí)偏移成像識(shí)別隱蔽儲(chǔ)層的方法，其特征在于，所述方法包括：

獲取遠(yuǎn)探測(cè)聲波測(cè)井資料以及與之相關(guān)的常規(guī)及電成像測(cè)井資料；

對(duì)所述遠(yuǎn)探測(cè)聲波測(cè)井資料進(jìn)行預(yù)處理，并提取縱、橫波時(shí)差信息；

利用多尺度相關(guān)對(duì)比方法在所述遠(yuǎn)探測(cè)聲波測(cè)井資料中提取反射波信息；

根據(jù)所述縱、橫波時(shí)差信息以及反射波信息，利用地震疊前逆時(shí)偏移算法，生成逆時(shí)偏移成像結(jié)果；

根據(jù)所述的逆時(shí)偏移成像結(jié)果，結(jié)合所述常規(guī)及電成像測(cè)井資料，確定所述隱蔽儲(chǔ)層所在層位及深度。

2.根據(jù)權(quán)利要求1所述的井下逆時(shí)偏移成像識(shí)別隱蔽儲(chǔ)層的方法，其特征在于，所述對(duì)所述遠(yuǎn)探測(cè)聲波測(cè)井資料進(jìn)行預(yù)處理，包括：

對(duì)所述遠(yuǎn)探測(cè)聲波測(cè)井資料進(jìn)行增益恢復(fù)、濾除斯通利波及低頻噪聲。

3.根據(jù)權(quán)利要求1所述的井下逆時(shí)偏移成像識(shí)別隱蔽儲(chǔ)層的方法，其特征在于，所述利用多尺度相關(guān)對(duì)比方法在所述遠(yuǎn)探測(cè)聲波測(cè)井資料中提取反射波信息，具體包括：

對(duì)共炮點(diǎn)道集做雙樹復(fù)小波變換，其次再對(duì)不同道上、同一尺度的分解系數(shù)進(jìn)行抽取組成共尺度道集，進(jìn)而實(shí)現(xiàn)對(duì)每一個(gè)共尺度道集做雙樹復(fù)小波反變換，得到共尺度道集的時(shí)域形式，最后再對(duì)每一時(shí)域的共尺度道集做慢度時(shí)間相關(guān)法，得到不同尺度的慢度時(shí)間相關(guān)法圖，從而有效提取出所述反射波信息。

4.根據(jù)權(quán)利要求1所述的井下逆時(shí)偏移成像識(shí)別隱蔽儲(chǔ)層的方法，其特征在于，在根據(jù)所述縱、橫波時(shí)差信息以及反射波信息，根據(jù)疊前逆時(shí)偏移算法，生成儲(chǔ)層成像結(jié)果的步驟中，包括：

觀測(cè)系統(tǒng)的轉(zhuǎn)換：分析遠(yuǎn)探測(cè)聲波系統(tǒng)采集系統(tǒng)與二維單邊地震采集系統(tǒng)的異同，將所述二維單邊地震采集系統(tǒng)做90°轉(zhuǎn)換，同時(shí)將所述二維單邊地震采集系統(tǒng)的參數(shù)按照遠(yuǎn)探測(cè)聲波激發(fā)、接收關(guān)系與所述遠(yuǎn)探測(cè)聲波系統(tǒng)采集系統(tǒng)進(jìn)行一一對(duì)應(yīng)；

尺度轉(zhuǎn)換：采用放大計(jì)算目標(biāo)體N倍的同時(shí)相應(yīng)縮小主頻N倍的方法，將所述遠(yuǎn)探測(cè)聲波測(cè)井資料的主頻降低到地震頻帶范圍；

偏移速度模型重構(gòu)：將速度網(wǎng)格在空間尺度上放大N倍，根據(jù)所選取的有限差分階數(shù)確定網(wǎng)格差分密度，解決空間頻散問題。

5.根據(jù)權(quán)利要求4所述的井下逆時(shí)偏移成像識(shí)別隱蔽儲(chǔ)層的方法，其特征在于，在完成所述觀測(cè)系統(tǒng)轉(zhuǎn)換、尺度轉(zhuǎn)換以及偏移速度模型重構(gòu)后，采用高階有限差分法求解帶PML邊界的二維聲波方程完成波場(chǎng)正演和波場(chǎng)逆推，其計(jì)算公式為：

∂px∂t+d(x)px=v2∂Ax∂x∂pz∂t+d(z)pz=v2∂Az∂z∂Ax∂t+d(x)Ax=∂px∂x+∂pz∂x∂Az∂t+d(z)Az=∂px∂z+∂pz∂z;]]>

其中，聲波波場(chǎng)為p=p(x,z,t)，速度場(chǎng)為v=v(x,z)，將所述聲波波場(chǎng)p=p(x,z,t)分解為x和z兩個(gè)方向：p=p^x+p^z，并引入中間變量A=A^x+A^z，d(x),d(z)為所述PML邊界中的吸收參數(shù)。