專(zhuān)利名稱(chēng)：測(cè)井曲線(xiàn)與地震記錄全波匹配方法及裝置的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及石油地球物理勘探地震模型正演模擬領(lǐng)域，特別涉及一種測(cè)井曲線(xiàn)與 地震記錄全波匹配方法及裝置。
背景技術(shù)：
在油氣勘探開(kāi)發(fā)過(guò)程中，地球物理測(cè)井?dāng)?shù)據(jù)常用于約束地震數(shù)據(jù)進(jìn)行油氣儲(chǔ)層預(yù) 測(cè)，地球物理測(cè)井?dāng)?shù)據(jù)與地震數(shù)據(jù)匹配的技術(shù)(以下簡(jiǎn)稱(chēng)井地匹配技術(shù))精度直接影響油 氣儲(chǔ)層預(yù)測(cè)精度?，F(xiàn)有井地匹配技術(shù)屬于褶積模型匹配技術(shù)，即將聲波測(cè)井?dāng)?shù)據(jù)轉(zhuǎn)換成聲波速度， 結(jié)合密度測(cè)井?dāng)?shù)據(jù)生成反射系數(shù)序列，反射系數(shù)序列與子波的褶積形成井旁地震記錄。這種由褶積方法形成的地震記錄僅僅是地層一次反射波的結(jié)果。主要存在以下 問(wèn)題(1)沒(méi)有考慮地震波傳播方向并不垂直于地層，事實(shí)上，得到的疊后地震記錄是由多 角度反射地震波疊加而成的；( 沒(méi)有考慮層間多次反射波的存在，散射衰減效應(yīng)被忽略； (3)沒(méi)有考慮轉(zhuǎn)換波的存在，即地震波在介質(zhì)分界面上存在縱波與橫波能量的相互轉(zhuǎn)化。因此在實(shí)現(xiàn)本發(fā)明的過(guò)程中，發(fā)明人發(fā)現(xiàn)現(xiàn)有技術(shù)的缺陷在于技術(shù)精度不高，即 合成地震記錄與實(shí)際井旁地震記錄相似度較低(一般低于90% )。

發(fā)明內(nèi)容
本發(fā)明實(shí)施例提供一種測(cè)井曲線(xiàn)與地震記錄全波匹配方法及裝置，目的在于提升 井地?cái)?shù)據(jù)匹配精度。為達(dá)到上述目的，本發(fā)明實(shí)施例提供一種測(cè)井曲線(xiàn)與地震記錄全波匹配的方法， 該方法包括獲取地震記錄；利用反射率法對(duì)地震記錄進(jìn)行處理，獲得井旁道地震記錄；將 井旁道地震記錄與測(cè)井曲線(xiàn)進(jìn)行匹配。本發(fā)明實(shí)施例還提供一種測(cè)井曲線(xiàn)與地震記錄全波匹配的裝置，該裝置包括獲取單元，獲取地震記錄；處理單元，利用反射率法對(duì)地震記錄進(jìn)行處理，獲得井旁道地震記錄；匹配單元，將井旁道地震記錄與測(cè)井曲線(xiàn)進(jìn)行匹配。本發(fā)明實(shí)施例的有益效果在于，通過(guò)在進(jìn)行測(cè)井?dāng)?shù)據(jù)與地震數(shù)據(jù)匹配時(shí)，同時(shí)考 慮多角度反射地震波、層間多次波及界面轉(zhuǎn)換波，可以提升井地?cái)?shù)據(jù)匹配精度。


此處所說(shuō)明的附圖用來(lái)提供對(duì)本發(fā)明的進(jìn)一步理解，構(gòu)成本申請(qǐng)的一部分，并不 構(gòu)成對(duì)本發(fā)明的限定。在附圖中圖1是本發(fā)明實(shí)施例的測(cè)井曲線(xiàn)與地震記錄全波匹配方法的流程圖；圖2是本發(fā)明實(shí)施例的層間多次波的示意圖；圖3是本發(fā)明實(shí)施例的任一區(qū)域的反射和透射情況的示意圖4是本發(fā)明實(shí)施例的單層區(qū)域的含多次波的演示示意圖；圖5是本發(fā)明實(shí)施例的雙層區(qū)域的傳播演示示意圖；圖6是本發(fā)明實(shí)施例的地震波場(chǎng)的合成圖；圖7是本發(fā)明實(shí)施例的頻率漢明窗和慢度漢明窗示意圖；圖8是本發(fā)明實(shí)施例的Ricker子波時(shí)間域波形示意圖；圖9是本發(fā)明實(shí)施例的Ricker子波頻率域波形示意圖；圖10是本發(fā)明實(shí)施例的由模擬測(cè)井?dāng)?shù)據(jù)提取的反射系數(shù)、褶積記錄、反射率法正 演結(jié)果的示意圖；圖11是本發(fā)明實(shí)施例的三層模型一次波、多次波和轉(zhuǎn)換波記錄的示意圖；圖12是本發(fā)明實(shí)施例的np5-4井速度和密度測(cè)井資料及井旁道地震記錄的示意 圖；圖13是本發(fā)明實(shí)施例的對(duì)np5-4井測(cè)井?dāng)?shù)據(jù)分層及用反射率法模擬記錄的示意 圖；圖14是本發(fā)明實(shí)施例的測(cè)井曲線(xiàn)與地震記錄全波匹配裝置的構(gòu)成圖；圖15是本發(fā)明實(shí)施例的處理單元的構(gòu)成圖。
具體實(shí)施例方式為使本發(fā)明的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，下面結(jié)合附圖對(duì)本發(fā)明實(shí)施 例作進(jìn)一步詳細(xì)說(shuō)明。在此，本發(fā)明的示意性實(shí)施例及其說(shuō)明用于解釋本發(fā)明，但并不作為 對(duì)本發(fā)明的限定。本發(fā)明實(shí)施例提供一種測(cè)井曲線(xiàn)與地震記錄全波匹配的方法，如圖1所示，所述 方法包括步驟101，獲取地震記錄；步驟102，利用反射率法對(duì)地震記錄進(jìn)行處理，獲得井旁道地震記錄；步驟103，將井旁道地震記錄與測(cè)井曲線(xiàn)進(jìn)行匹配。在地震勘探資料中，除了一次反射波、隨機(jī)干擾以外，還存在著多次波和轉(zhuǎn)換波。 層間多次波又稱(chēng)為內(nèi)部多次波，是指所有下行反射發(fā)生在除自由界面以外的其它反射界面 的多次波。圖2是層間多次波的示意圖，如圖2所示，內(nèi)部多次波的反射界面通常很難精確確 定，因此內(nèi)部多次波很難預(yù)測(cè)，也較難進(jìn)行有效地消除。此外內(nèi)部多次波具有與有效波相似 甚至更高的速度，這使得內(nèi)部多次波的壓制更為困難。當(dāng)假設(shè)條件不能滿(mǎn)足時(shí)，多次波就不 但不能被消除，壓制措施甚至可能損傷有效波。在實(shí)際資料的處理和解釋中，多次波的類(lèi)型一般根據(jù)對(duì)時(shí)間、速度的目測(cè)和經(jīng)驗(yàn) 來(lái)確定。但在有些多次波與有效波沒(méi)有明顯的差異，常常難以確定。層間多次波同相軸通常 以不同的正常時(shí)差迭加在一次波同相軸之上，和一次有效反射波相干涉，就往往給處理和 解釋帶來(lái)很多麻煩。如果淺、中層存在良好的反射界面并產(chǎn)生多次波，就有可能掩蓋了中、 深層的一次反射波。在剖面上多次波較強(qiáng)時(shí)，如果在解釋中不能正確地把多次波識(shí)別出來(lái)， 就會(huì)造成錯(cuò)誤的地質(zhì)解釋。在本實(shí)施例中，采用反射率法正演模擬地震記錄，然后做試驗(yàn)分析。反射率法是應(yīng)
4用最廣泛的疊前地震記錄計(jì)算方法，它能夠提供波場(chǎng)的完全解，包括地震波在地層內(nèi)傳播 時(shí)形成的一次波、多次波和轉(zhuǎn)換波。以下首先對(duì)反射率法的原理進(jìn)行簡(jiǎn)單的介紹?；镜姆瓷渎仕惴M的層狀介質(zhì)，其上邊界是半空間或自由表面，下邊界為半 空間。在本實(shí)施例中，討論局限于圓柱坐標(biāo)系統(tǒng)，χ軸(徑向)和ζ軸(縱向)?？臻g為 均勻介質(zhì)，而且設(shè)χ為射線(xiàn)參數(shù)，Z為深度參數(shù)。激發(fā)和接收點(diǎn)可以在模型的任意處。反射率方法和其它的計(jì)算合成記錄方法一樣，是一種對(duì)彈性動(dòng)力學(xué)方程- PCO1U = Nξ + f(L D和關(guān)系式
‘ (1.2)求積分的方法。這里的ξ和ζ是無(wú)限小的壓力和張力。P是密度，U是彈性模量。而f是體力 密度。3是個(gè)四階的彈性張量，在柱空間的均勻介質(zhì)最通常的形式看作是一個(gè)關(guān)于垂直軸線(xiàn) 對(duì)稱(chēng)的橫向同性體。3因此有五個(gè)獨(dú)立的參數(shù)。一個(gè)經(jīng)常被忽略的方面是式(1. 2)根本上是個(gè)頻率域的關(guān)系式，既然3的組成通 常是復(fù)雜的頻率依靠的量。這樣假設(shè)以時(shí)間為變量的一個(gè)暫時(shí)的傅立葉變換
權(quán)利要求
1.一種測(cè)井曲線(xiàn)與地震記錄全波匹配的方法，其特征在于，所述方法包括 獲取地震記錄；利用反射率法對(duì)所述地震記錄進(jìn)行處理，獲得井旁道地震記錄； 將所述井旁道地震記錄與測(cè)井曲線(xiàn)進(jìn)行匹配。
2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述利用反射率法對(duì)所述地震記錄進(jìn)行 處理，獲得井旁道地震記錄，具體包括在頻率慢度域求解反射透射系數(shù)響應(yīng)；根據(jù)所述反射透射系數(shù)響應(yīng)進(jìn)行慢度積分和頻率積分，得到所述井旁地震記錄。
3.根據(jù)權(quán)利要求2所述的方法，其特征在于，所述在頻率慢度域求解反射透射系數(shù)響 應(yīng)，具體包括分別求取自由界面、固體-固體界面和任意區(qū)域的反射透射系數(shù)； 將所述反射透射系數(shù)合成為總的反射透射系數(shù)響應(yīng)。
4.根據(jù)權(quán)利要求2所述的方法，其特征在于，根據(jù)所述反射透射系數(shù)響應(yīng)進(jìn)行慢度積 分和頻率積分，具體包括在頻率域?qū)⑿盘?hào)與漢明窗相乘；在慢度域采用漢明窗與反射率函數(shù)相乘來(lái)消除截?cái)嘈?yīng)。
5.一種測(cè)井曲線(xiàn)與地震記錄全波匹配的裝置，其特征在于，所述裝置包括 獲取單元，獲取地震記錄；處理單元，利用反射率法對(duì)所述地震記錄進(jìn)行處理，獲得井旁道地震記錄； 匹配單元，將所述井旁道地震記錄與測(cè)井曲線(xiàn)進(jìn)行匹配。
6.根據(jù)權(quán)利要求5所述的裝置，其特征在于，所述處理單元具體包括 求解單元，在頻率慢度域求解反射透射系數(shù)響應(yīng)；積分單元，根據(jù)所述反射透射系數(shù)響應(yīng)進(jìn)行慢度積分和頻率積分，得到所述井旁地震記錄。
7.根據(jù)權(quán)利要求6所述的裝置，其特征在于，所述求解單元具體用于分別求取自由界 面、固體-固體界面和任意區(qū)域的反射透射系數(shù)；將所述反射透射系數(shù)合成為總的反射透 射系數(shù)響應(yīng)。
8.根據(jù)權(quán)利要求6所述的裝置，其特征在于，所述積分單元具體用于在頻率域?qū)⑿盘?hào) 與漢明窗相乘；在慢度域采用漢明窗與反射率函數(shù)相乘來(lái)消除截?cái)嘈?yīng)。
全文摘要
本發(fā)明實(shí)施例提供一種測(cè)井曲線(xiàn)與地震記錄全波匹配的方法及裝置，該方法包括獲取地震記錄；利用反射率法對(duì)地震記錄進(jìn)行處理，獲得井旁道地震記錄；將井旁道地震記錄與測(cè)井曲線(xiàn)進(jìn)行匹配。通過(guò)本發(fā)明實(shí)施例，可以同時(shí)考慮多角度反射地震波、層間多次波及界面轉(zhuǎn)換波，提升井地?cái)?shù)據(jù)匹配精度。
文檔編號(hào)G01V1/28GK102096100SQ20101055921
公開(kāi)日2011年6月15日 申請(qǐng)日期2010年11月24日 優(yōu)先權(quán)日2010年11月24日
發(fā)明者王尚旭 申請(qǐng)人:中國(guó)石油大學(xué)(北京), 中國(guó)石油天然氣集團(tuán)公司