一種直接生成共轉(zhuǎn)換點(diǎn)道集的方法
【專利摘要】本發(fā)明提供一種直接生成共轉(zhuǎn)換點(diǎn)道集的方法,假設(shè)轉(zhuǎn)換波地震數(shù)據(jù)已預(yù)處理好,已獲得準(zhǔn)確的動(dòng)校正參數(shù)和計(jì)算轉(zhuǎn)換點(diǎn)位置的速度等參數(shù),按炮檢距分組,分時(shí)窗疊加獲得時(shí)間和空間變化的共轉(zhuǎn)換點(diǎn)道集。本發(fā)明實(shí)現(xiàn)簡便、適應(yīng)大數(shù)據(jù)量處理。
【專利說明】一種直接生成共轉(zhuǎn)換點(diǎn)道集的方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及轉(zhuǎn)換波地震資料處理【技術(shù)領(lǐng)域】,具體地,涉及一種直接時(shí)空變生成共轉(zhuǎn)換點(diǎn)道集的方法。
【背景技術(shù)】
[0002]在多分量轉(zhuǎn)換波地震資料處理中,共轉(zhuǎn)換點(diǎn)道集的抽取方法一直是一項(xiàng)關(guān)鍵技術(shù)。因?yàn)檗D(zhuǎn)換波的下行波是縱波、上行波是橫波,造成傳播路徑的不對(duì)稱性,轉(zhuǎn)換點(diǎn)位于從中心點(diǎn)偏向接收點(diǎn)的一側(cè),且隨深度及界面上下縱橫波速度的變化而變化。這使得基于常規(guī)共中心點(diǎn)(CMP)道集的速度分析、動(dòng)校正和水平疊加技術(shù)不能適用于轉(zhuǎn)換波,轉(zhuǎn)換波疊加成像必須按共轉(zhuǎn)換點(diǎn)進(jìn)行道集分選。目前,有關(guān)轉(zhuǎn)換點(diǎn)位置的計(jì)算已有很多方法,主要包括漸近線法(ACP)、單層各向同性介質(zhì)迭代解、多層介質(zhì)指定深度各向同性近似解(Thomsen,1999)、多層介質(zhì)指定旅行時(shí)各向同性近似解(Thomsen, 1999)和多層VTI介質(zhì)指定旅行時(shí)近似解(Yuan&Li,2001)等。另外,周竹生等(1993)導(dǎo)出了一種快速、高精度轉(zhuǎn)換點(diǎn)的迭代算法;許士勇等(2002)給出了四次方程的近似表達(dá)式;姚陳等(2005)給出了三維傾角轉(zhuǎn)換波共轉(zhuǎn)換點(diǎn)的解析表達(dá)式;苑春方等(2005,2006,2006)推導(dǎo)出了單層各向同性介質(zhì)與深度有關(guān)的唯一的解析解和含有旅行時(shí)間唯一解析解,傾斜界面上的轉(zhuǎn)換點(diǎn)的精確解;Chang等(2005)利用傾斜界面上轉(zhuǎn)換波射線路徑的三角幾何關(guān)系,推導(dǎo)出傾斜界面上的轉(zhuǎn)換點(diǎn)的近似公式;畢麗飛等(2005)用水平界面推導(dǎo)轉(zhuǎn)換點(diǎn)坐標(biāo)迭代解類似的方法。但目前,共轉(zhuǎn)換點(diǎn)道集的選排方法仍然是以整道選排為主,主要是漸近逼近道分類法、依賴深度的選排方法。漸近逼近道分類法是通過一個(gè)工區(qū)的固定速度比,計(jì)算轉(zhuǎn)換點(diǎn)的漸近線,然后將整個(gè)工區(qū)的轉(zhuǎn)換波數(shù)據(jù)抽取到漸近線道集中,即ACP道集,這種方法能一定程度上能夠滿足深部地層的成像要求,但是淺層的成像效果較差。依賴深度的選排方法只對(duì)某一深度最為正確,都是能在目標(biāo)層聚焦成像,更適合目的層較淺的轉(zhuǎn)換波處理,更適合更新速度分析,并且指定目的層計(jì)算提供了多種轉(zhuǎn)換點(diǎn)位置的計(jì)算方法,但在所選反射層的上面或下面,反射界面的疊后成像可能會(huì)變得模糊。另外,分層選排方法能夠使不同深度的轉(zhuǎn)換波均能達(dá)到較好的聚焦效果,但需要在不同層之間進(jìn)行道集的拼接,特別是在大數(shù)據(jù)量的情況下實(shí)際操作困難。
[0003]上述的共轉(zhuǎn)換點(diǎn)道集的選排方法都有一定的應(yīng)用條件,并不能很好的滿足轉(zhuǎn)換波對(duì)轉(zhuǎn)換點(diǎn)在時(shí)間上和空間上變化的要求。
【發(fā)明內(nèi)容】
[0004]本發(fā)明目的是提供一種能適應(yīng)時(shí)間上和空間上變化的共轉(zhuǎn)換點(diǎn)道集的選排方法,一種直接時(shí)空變生成共轉(zhuǎn)換點(diǎn)道集的方法。
[0005]本發(fā)明通過以下技術(shù)方案實(shí)現(xiàn):
[0006]一種直接生成共轉(zhuǎn)換點(diǎn)道集的方法,具體步驟包括:
[0007]I)采集原始轉(zhuǎn)換波地震數(shù)據(jù),經(jīng)預(yù)處理和速度分析,得到轉(zhuǎn)換波動(dòng)校正和共轉(zhuǎn)換點(diǎn)位置計(jì)算需要的參數(shù),并對(duì)預(yù)處理后轉(zhuǎn)換波地震數(shù)據(jù)進(jìn)行動(dòng)校正,且通過切除消除動(dòng)校拉伸的影響,
[0008]步驟1)所述的預(yù)處理包括加載觀測(cè)系統(tǒng)、靜校正、波場(chǎng)分離、疊前去噪、振幅恢復(fù)、預(yù)置共轉(zhuǎn)換點(diǎn)CCP號(hào);
[0009]2)步驟1)處理后的轉(zhuǎn)換波地震數(shù)據(jù)按炮檢距大小進(jìn)行分組,每一組對(duì)應(yīng)一個(gè)炮檢距范圍;
[0010]3)輸入一組炮檢距范圍內(nèi)的數(shù)據(jù);
[0011]4)按道頭中的共轉(zhuǎn)換點(diǎn)CCP號(hào)讀取轉(zhuǎn)換波用來計(jì)算轉(zhuǎn)換點(diǎn)位置的參數(shù);
[0012]5)對(duì)地震道進(jìn)行分時(shí)窗處理,
[0013]步驟5)所述的分時(shí)窗處理是指將地震道分成多個(gè)時(shí)窗,時(shí)窗長度不變,時(shí)窗重疊,重疊時(shí)窗的重疊長度為時(shí)窗內(nèi)樣點(diǎn)數(shù)的百分比perc, perc取值為0-100 ;
[0014]6)從第一個(gè)時(shí)窗開始,計(jì)算時(shí)窗中點(diǎn)的共轉(zhuǎn)換點(diǎn)CCP號(hào),并由這個(gè)共轉(zhuǎn)換點(diǎn)號(hào)重新獲得該時(shí)窗用來計(jì)算轉(zhuǎn)換點(diǎn)位置的參數(shù),之后重新計(jì)算時(shí)窗中點(diǎn)相對(duì)準(zhǔn)確的共轉(zhuǎn)換點(diǎn)CCP號(hào),直接把時(shí)窗內(nèi)的樣點(diǎn)值與窗函數(shù)相乘疊在共轉(zhuǎn)換點(diǎn)地震道上,即轉(zhuǎn)換波共反射面元的寬度上,相應(yīng)共轉(zhuǎn)換點(diǎn)地震道每個(gè)樣點(diǎn)的覆蓋次數(shù)加上窗函數(shù)的數(shù)值,
[0015]步驟6)所述的計(jì)算轉(zhuǎn)換點(diǎn)位置的參數(shù)是預(yù)先速度分析得到的,是轉(zhuǎn)換波均方根速度ve2、縱波均方根速度各向異性參數(shù)Xrff、縱橫波垂直速度比Y0和縱橫波有效速度
比 reff,
[0016]步驟6)所述的共轉(zhuǎn)換點(diǎn)CCP號(hào),是指按工區(qū)面元?jiǎng)澐值?,與常規(guī)處理中縱波共中心點(diǎn)CMP號(hào)一致,
[0017]步驟6)所述的共轉(zhuǎn)換點(diǎn)CCP號(hào)是由下述方法計(jì)算的:
[0018]xc = xm+D (1)
[0019]式中xM、xc為共中心點(diǎn)CMP和共轉(zhuǎn)換點(diǎn)CCP點(diǎn)的水平坐標(biāo),D為轉(zhuǎn)換點(diǎn)偏離中心點(diǎn)的水平距離,
[0020]D = xp -x/2(2)
[0021]式中Xp為共轉(zhuǎn)換點(diǎn)與炮點(diǎn)之間的水平距離,
[0022]已知共中心點(diǎn)CMP號(hào)As,計(jì)算相應(yīng)的共轉(zhuǎn)換點(diǎn)CCP號(hào)jxc:
[0023]
【權(quán)利要求】
1.一種直接生成共轉(zhuǎn)換點(diǎn)道集的方法,具體步驟包括: 1)采集原始轉(zhuǎn)換波地震數(shù)據(jù),經(jīng)預(yù)處理和速度分析,得到轉(zhuǎn)換波動(dòng)校正和共轉(zhuǎn)換點(diǎn)位置計(jì)算需要的參數(shù),并對(duì)預(yù)處理后轉(zhuǎn)換波地震數(shù)據(jù)進(jìn)行動(dòng)校正,且通過切除消除動(dòng)校拉伸的影響, 步驟I)所述的預(yù)處理包括加載觀測(cè)系統(tǒng)、靜校正、波場(chǎng)分離、疊前去噪、振幅恢復(fù)、預(yù)置共轉(zhuǎn)換點(diǎn)CCP號(hào); 2)步驟I)處理后的轉(zhuǎn)換波地震數(shù)據(jù)按炮檢距大小進(jìn)行分組,每一組對(duì)應(yīng)一個(gè)炮檢距范圍; 3)輸入一組炮檢距范圍內(nèi)的數(shù)據(jù); 4)按道頭中的共轉(zhuǎn)換點(diǎn)CCP號(hào)讀取轉(zhuǎn)換波用來計(jì)算轉(zhuǎn)換點(diǎn)位置的參數(shù); 5)對(duì)地震道進(jìn)行分時(shí)窗處理, 步驟5)所述的分時(shí)窗處理是指將地震道分成多個(gè)時(shí)窗,時(shí)窗長度不變,時(shí)窗重疊,重疊時(shí)窗的重疊長度為時(shí)窗內(nèi)樣點(diǎn)數(shù)的百分比perc, perc取值為0-100 ; 6)從第一個(gè)時(shí)窗開始,計(jì)算時(shí)窗中點(diǎn)的共轉(zhuǎn)換點(diǎn)CCP號(hào),并由這個(gè)共轉(zhuǎn)換點(diǎn)號(hào)重新獲得該時(shí)窗用來計(jì)算轉(zhuǎn)換點(diǎn)位置的參數(shù),之后重新計(jì)算時(shí)窗中點(diǎn)相對(duì)準(zhǔn)確的共轉(zhuǎn)換點(diǎn)CCP號(hào),直接把時(shí)窗內(nèi)的樣點(diǎn)值與窗函數(shù)相乘疊在共轉(zhuǎn)換點(diǎn)地震道上,即轉(zhuǎn)換波共反射面元的寬度上,相應(yīng)共轉(zhuǎn)換點(diǎn)地震道每個(gè)樣點(diǎn)的覆蓋次數(shù)加上窗函數(shù)的數(shù)值, 步驟6)所述的計(jì)算轉(zhuǎn)換點(diǎn)位置的參數(shù)是預(yù)先速度分析得到的,是轉(zhuǎn)換波均方根速度D縱波均方根速度%、各向異性參數(shù)Xeff、縱橫波垂直速度比Ytl和縱橫波有效速度比reff, 步驟6)所述的共轉(zhuǎn)換點(diǎn)CCP號(hào),是指按工區(qū)面元?jiǎng)澐值?,與常規(guī)處理中縱波共中心點(diǎn)CMP號(hào)一致, 步驟6)所述的共轉(zhuǎn)換點(diǎn)CCP號(hào)是由下述方法計(jì)算的:
xC = xm+D 式中xM、x。為共中心點(diǎn)CMP和共轉(zhuǎn)換點(diǎn)CCP點(diǎn)的水平坐標(biāo),D為轉(zhuǎn)換點(diǎn)偏離中心點(diǎn)的水平距離, 式中Xp為共轉(zhuǎn)換點(diǎn)與炮點(diǎn)之間的水平距離, 已知共中心點(diǎn)CMP號(hào)4*,計(jì)算相應(yīng)的共轉(zhuǎn)換點(diǎn)CCP號(hào)Lc'
【文檔編號(hào)】G01V1/28GK103645498SQ201310553398
【公開日】2014年3月19日 申請(qǐng)日期:2013年11月8日 優(yōu)先權(quán)日:2013年11月8日
【發(fā)明者】陳海峰, 李向陽, 狄?guī)妥? 魏建新, 蔡志光, 陳雙全 申請(qǐng)人:中國石油大學(xué)(北京)