基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法
【專利摘要】本發(fā)明提供了一種基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法,屬于地球物理地震勘探領(lǐng)域。所述方法包括:步驟1:X、Y、Z三分量數(shù)據(jù)輸入;步驟2:縱、橫波道頭加載、道編輯;步驟3:水平分量偏振旋轉(zhuǎn);步驟4:三分量數(shù)據(jù)分選;步驟5:縱、橫波近地表靜校正;步驟6:縱、橫波信號處理;步驟7:縱、橫波CDP道集分選;步驟8:縱、橫波速度分析;步驟9:縱、橫波剩余靜校正;步驟10:縱、橫波疊前去噪;步驟11:縱、橫波PSTM偏移;步驟12:縱、橫波疊加和修飾處理;步驟13:縱、橫波成果數(shù)據(jù)輸出。
【專利說明】基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于地球物理地震勘探領(lǐng)域,具體涉及一種基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法。
【背景技術(shù)】
[0002]多分量勘探的資料采集技術(shù)需要三分量檢波器、縱波震源和橫波震源。在縱波震源激發(fā)情況下,三分量檢波器主要接收P-P、P-SV、P-SH波,在橫波震源激發(fā)情況下,又可以分SV波激發(fā)和SH波激發(fā),并可以分別接收到SV-P、SV-SV和SV-SH波以及SH_P、SH-SV和SH-SH波。在地層為各向同性介質(zhì)條件下,一種震源激發(fā)的地震波場中,一般SV波和SH波不能共存,但在各向異性介質(zhì)條件下,可以同時(shí)記錄到SV、SH波。
[0003]縱波(P)的質(zhì)點(diǎn)振動方向平行于波的傳播方向。橫波(S)質(zhì)點(diǎn)振動方向垂直于波的傳播方向。SV是射線平面內(nèi)振動的橫波(垂直偏振橫波)。SH波是水平面內(nèi)振動的橫波(水平偏振橫波)。轉(zhuǎn)換波(PS、SP)是由不同振型的波場轉(zhuǎn)換為另外一種振型的波場。
[0004]多波勘探可研究工區(qū)內(nèi)地震屬性與含油氣性地層的關(guān)系,總結(jié)多波地震屬性在含油氣地層中的響應(yīng);可綜合利用縱波和橫波地震信息,提取物性參數(shù)、巖性參數(shù)信息,對研究區(qū)主要目的層進(jìn)行多波地震屬性含氣性特征響應(yīng)分析。多波勘探因在巖性、裂隙和流體識別等方面獲得成功,而備受關(guān)注。
[0005]由于橫波勘探的成本高,生產(chǎn)應(yīng)用不普遍,當(dāng)前多數(shù)地區(qū)多波勘探更多的仍然是利用轉(zhuǎn)換波信息進(jìn)行勘探。即采用炸藥震源激發(fā),三分量檢波器接收獲得縱波和轉(zhuǎn)換波的地震勘探,資料處理和解釋以縱波和轉(zhuǎn)換波為主要研究對象。
[0006]但在某些探區(qū),其表層各向異性顯著、速度變化劇烈,在這特殊的地震地質(zhì)條件下,炸藥震源激發(fā)三分量檢波器接收,能夠采集到震源橫波,而且其能量強(qiáng)于傳統(tǒng)的轉(zhuǎn)換波,這一發(fā)現(xiàn)為在這些探區(qū)進(jìn)行震源橫波地震勘探成為可能,但是目前特殊條件下炸藥震源激發(fā)三分量接收的橫波勘探方法還沒有被開發(fā)利用于生產(chǎn),沒有形成炸藥震源橫波地震勘探的配套處理方法。
[0007]震源橫波的產(chǎn)生機(jī)理復(fù)雜,其資料處理主要存在4大難題:①震源橫波的振動方向不穩(wěn)定,炮、檢點(diǎn)間記錄的一致性差,對資料的一致性影響較大。②震源橫波能量總體較弱,干擾波嚴(yán)重,信噪比偏低。③復(fù)雜地表?xiàng)l件由于橫波靜校正受整個(gè)風(fēng)化層的影響,且速度較縱波低,其靜校正量往往變化劇烈,導(dǎo)致靜校正問題突出。④波場復(fù)雜,難以獲得準(zhǔn)確的速度。
【發(fā)明內(nèi)容】
[0008]本發(fā)明的目的在于解決上述現(xiàn)有技術(shù)中存在的難題,提供一種基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法,獲得較轉(zhuǎn)換波剖面品質(zhì)更好的震源橫波剖面,進(jìn)而提高處理成果解釋的可靠性。
[0009]本發(fā)明是通過以下技術(shù)方案實(shí)現(xiàn)的:
[0010]一種基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法,包括:
[0011]步驟1:X、Y、Z三分量數(shù)據(jù)輸入:將野外采集的原始三分量炮集數(shù)據(jù)輸入到處理系統(tǒng);
[0012]步驟2:縱、橫波道頭加載、道編輯:根據(jù)炮點(diǎn)和檢波點(diǎn)信息產(chǎn)生空間屬性文件,將空間屬性文件加載到道頭里;人工進(jìn)行道記錄刪除或修改;生成起始時(shí)間;應(yīng)用網(wǎng)格定義將資料處理的范圍定義處理;輸出加載道頭的炮集原始地震數(shù)據(jù);
[0013]步驟3:水平分量偏振旋轉(zhuǎn):利用橫波的初至波對步驟2輸出的X分量和Y分量數(shù)據(jù)進(jìn)行偏振旋轉(zhuǎn),使得能量得以歸位;水平分量偏振旋轉(zhuǎn)是炸藥震源激發(fā)三分量接收的橫波處理的關(guān)鍵,只有利用橫波的初至波進(jìn)行水平分量偏振旋轉(zhuǎn),使得能量得以歸位后,資料的處理效果才有保證;
[0014]步驟4:三分量數(shù)據(jù)分選;
[0015]步驟5:縱、橫波近地表靜校正;
[0016]步驟6:縱、橫波信號處理:應(yīng)用地表一致性振幅補(bǔ)償、炮域噪聲衰減、地表一致性反褶積以及隨機(jī)噪聲衰減進(jìn)行縱、橫波信號處理;
[0017]步驟7:縱、橫波⑶P道集分選:將步驟6輸出的炮集記錄數(shù)據(jù)分選成按⑶P道集輸出;
[0018]步驟8:縱、橫波速度分析:產(chǎn)生速度譜(包含速度函數(shù)的中間處理成果),然后對速度譜進(jìn)行拾取速度得到速度文件,再從拾取的速度文件生成插值速度場;
[0019]步驟9:縱、橫波剩余靜校正:應(yīng)用剩余靜校正,生成用于目標(biāo)線偏移的速度;
[0020]步驟10:縱、橫波疊前去噪:對步驟9輸出的X分量和Z分量地震數(shù)據(jù)運(yùn)行疊前去噪、QC疊加,并從所得的疊加數(shù)據(jù)體產(chǎn)生切除時(shí)間;縱、橫波疊前去噪就是針對資料信噪比低的情況,利用多域波場分離技術(shù)實(shí)現(xiàn)震源橫波凈化去噪,尤其是重點(diǎn)提高震源橫波的信噪比;
[0021]步驟11:縱、橫波PSTM偏移:進(jìn)行目標(biāo)線的偏移;拾取、生成疊前時(shí)間體偏移的速度場;用PSTM速度體運(yùn)行整塊數(shù)據(jù)體的偏移;拾取、生成疊前時(shí)間偏移后的疊加速度場;
[0022]步驟12:縱、橫波疊加和修飾處理:用疊前時(shí)間偏移后的疊加速度場,產(chǎn)生一個(gè)疊前時(shí)間偏移后的疊加;應(yīng)用疊后剩余振幅補(bǔ)償以便平衡振幅;應(yīng)用濾波、去噪處理得到最終信號增強(qiáng)的疊加;
[0023]步驟13:縱、橫波成果數(shù)據(jù)輸出:將步驟12得到的成果數(shù)據(jù)以SEGY格式輸出。
[0024]所述步驟I中,如果輸入的原始三分量炮集數(shù)據(jù)是按多路方式記錄的格式進(jìn)行記錄的,則首先對數(shù)據(jù)進(jìn)行解編,即對一個(gè)大的矩陣進(jìn)行變換,使變換后的矩陣的行能按地震道讀數(shù)據(jù)。
[0025]所述步驟2中道記錄刪除或修改是指:刪除不正常的道記錄和修改反極性的道記錄,以確保后續(xù)水平分量偏振旋轉(zhuǎn)處理的效果。
[0026]所述步驟4是這樣實(shí)現(xiàn)的:將步驟3輸出的混在一起的三分量地震記錄數(shù)據(jù)分選成X、Y、Z三個(gè)分量,分別輸出;后續(xù)處理只需對能量歸位好的X分量(震源橫波)和Z分量(縱波)進(jìn)行處理。
[0027]所述步驟5是這樣實(shí)現(xiàn)的:
[0028]從步驟4輸出的三個(gè)分量記錄上自動拾取縱、橫波的初至?xí)r間;對拾取時(shí)間表中自動拾取的初至?xí)r間進(jìn)行復(fù)查,并通過文本文件方式加載到道頭里;生成層析靜校正量;創(chuàng)建包含光滑CDP靜校正量的空間屬性文件。由于橫波的近地表靜校正的復(fù)雜性,有時(shí)需采用特殊手段進(jìn)行處理(如可以采用名稱為“震源橫波勘探的近地表靜校正處理方法”的專利中介紹的方法)。
[0029]與現(xiàn)有技術(shù)相比,本發(fā)明的有益效果是:通過對波場性質(zhì)的認(rèn)真細(xì)致地分析和研究,在正確的波場認(rèn)識指導(dǎo)下,該發(fā)明應(yīng)用于實(shí)際地震資料處理,采用成熟的處理技術(shù),獲得了較高信噪比的縱橫波剖面,且橫波剖面品質(zhì)遠(yuǎn)比該區(qū)的轉(zhuǎn)換波剖面好,橫波剖面與縱波剖面面貌整體一致性較好,反射波組特征明顯,可用于地質(zhì)解釋。
【專利附圖】
【附圖說明】
[0030]圖1a是原始三分量單炮記錄中的X分量。
[0031]圖1b是原始三分量單炮記錄中的Y分量。
[0032]圖1c是原始三分量單炮記錄中的Z分量。
[0033]圖2a是利用VSP測到的速度比進(jìn)行拉伸(加靜校正)的Z分量(縱波)。
[0034]圖2b是(加靜校正)的X分量(橫波)。
[0035](注:只需對縱波進(jìn)行拉伸,拉伸后便于兩者對比)
[0036]圖3是非零井源距VSP的Z分量記錄。
[0037]圖4a是X分量疊加剖面。
[0038]圖4b是Z分量疊加拉伸剖面。
[0039]圖5是本發(fā)明方法的步驟框圖。
【具體實(shí)施方式】
[0040]下面結(jié)合附圖對本發(fā)明作進(jìn)一步詳細(xì)描述:
[0041]通過對特定的地震地質(zhì)條件下的多波性質(zhì)分析研究,發(fā)現(xiàn)在近地表存在致密的高速巖層、各向異性等特定的地震地質(zhì)條件下,炸藥在近地表各向異性介質(zhì)中激發(fā)產(chǎn)生的橫波較轉(zhuǎn)換波更為發(fā)育,通過三分量接收,可獲得較轉(zhuǎn)換波品質(zhì)更好的震源橫波。
[0042]通過采用與縱波處理類同的處理流程(由于炸藥震源橫波與縱波的射線路徑都是遵循相同的反射規(guī)律的,除個(gè)別橫波處理特有的模塊外,用于縱波處理的模塊多數(shù)可以用于橫波處理;不像轉(zhuǎn)換波的射線路徑,縱波經(jīng)轉(zhuǎn)換后變成轉(zhuǎn)換波的反射路徑是不對稱的,故處理采取的手段也不同,更為復(fù)雜),可獲得較高品質(zhì)的震源橫波剖面,提高了其可解釋性,
[0043]由于橫波的資料處理技術(shù)(縱波的資料處理技術(shù)均可用于橫波的資料處理,如CDP道集分選等)較轉(zhuǎn)換波處理技術(shù)(如CCP道集分選)成熟,在同等原始資料品質(zhì)條件下,更容易獲得高品質(zhì)的橫波剖面。
[0044]因此,為充分發(fā)揮震源橫波在特殊探區(qū)地震勘探中的作用,在不增加采集成本的情況下,本發(fā)明給出了一套利用震源橫波信息進(jìn)行多波勘探的處理方法,以獲得較轉(zhuǎn)換波剖面品質(zhì)更好的震源橫波剖面,進(jìn)而提高處理成果解釋的可靠性。開辟了特殊探區(qū)多波地震勘探的新途徑。本發(fā)明為類似探區(qū)進(jìn)行多波勘探提供了新的方法,具有很好的應(yīng)用前景。
[0045](I)炸藥震源地震勘探中震源橫波的認(rèn)識
[0046]①通過對采集到的X、Y分量接收到的地震波場性質(zhì)進(jìn)行了研究分析,認(rèn)為在特定的地表非均質(zhì)條件下,在各向異性介質(zhì)中,炸藥震源能夠同時(shí)產(chǎn)生縱波波場和橫波波場(如圖1a至圖lc、圖2a、圖2b和圖3所示。)。
[0047]②研究的X、Y分量記錄中所接收到的優(yōu)勢波的性質(zhì)屬震源橫波,而記錄中的轉(zhuǎn)換波能量偏弱;震源橫波的能量強(qiáng)弱、振動方向與激發(fā)點(diǎn)的地震地質(zhì)條件有關(guān),同時(shí)也與檢波測線的方位角等其它因素有關(guān),不同方位角測線接收的震源橫波在X、Y分量記錄中的能量強(qiáng)弱規(guī)律不同。
[0048]③在多波勘探中,首先要深入了解多波的性質(zhì),只有對多波的性質(zhì)有了正確的認(rèn)識,才能充分挖掘出資料的潛力。震源橫波的正確認(rèn)識為采取正確的處理手段、獲得較高質(zhì)量的處理成果打下了良好的基礎(chǔ)。
[0049](2)震源橫波資料處理方法
[0050]如圖5所示,本發(fā)明方法的具體實(shí)施步驟如下:
[0051 ] 步驟1:X/Y/Z三分量數(shù)據(jù)輸入。
[0052]將野外采集的原始三分量炮集數(shù)據(jù)輸入到處理系統(tǒng)(例如使用生產(chǎn)中常用的Omega、CGGVeritas、Focus或Promax等等地震資料處理系統(tǒng)),如果輸入的野外數(shù)據(jù)是按多路方式記錄的格式(如SEGB、SEGD格式),則需解編將數(shù)據(jù)重排,在數(shù)學(xué)上就是對一個(gè)大的矩陣進(jìn)行變換,使變換后的矩陣的行能按地震道讀數(shù)據(jù)。
[0053]步驟2:縱/橫波道頭加載/道編輯。
[0054]根據(jù)炮點(diǎn)和檢波點(diǎn)信息產(chǎn)生空間屬性文件,將空間屬性文件加載到道頭里(具體實(shí)施時(shí),是運(yùn)行Omega地震資料處理系統(tǒng)的Geometry_update模塊);人工進(jìn)行道記錄刪除或修改(刪除不正常道記錄,修改反極性道記錄);生成起始時(shí)間;應(yīng)用網(wǎng)格定義(運(yùn)行Omega地震資料處理系統(tǒng)的Grid_define模塊,將資料處理的范圍定義出來);輸出加載道頭的炮集原始地震數(shù)據(jù)。
[0055]步驟3:水平分量偏振旋轉(zhuǎn)(步驟2的輸出的X分量、Y分量數(shù)據(jù))。
[0056]利用橫波的初至波進(jìn)行水平分量偏振旋轉(zhuǎn)(與運(yùn)行Omega地震資料處理系統(tǒng)的Rotate_2_comp模塊類似,根據(jù)X分量、Y分量偏振旋轉(zhuǎn)后X分量能量最強(qiáng)、Y分量能量最弱為準(zhǔn),主要得到含震源橫波的X分量),使得能量得以歸位,有利于解決波組的一致性和保幅問題。
[0057]步驟4:三分量數(shù)據(jù)分選。
[0058]將混在一起的三分量地震記錄數(shù)據(jù)(步驟3的輸出),分選成X、Y、Z三個(gè)分量(運(yùn)行Omega地震資料處理系統(tǒng)的TRACE-SELECT模塊),分別輸出以便后續(xù)處理。
[0059]步驟5:縱/橫波近地表靜校正。
[0060]從炮集記錄(指步驟4得到的分選后的數(shù)據(jù))上自動拾取初至?xí)r間(運(yùn)行Omega地震資料處理系統(tǒng)的Shot_Correl_FB_pick模塊),以便用于靜校正估算;對拾取時(shí)間表中自動拾取的初至?xí)r間進(jìn)行復(fù)查,并通過文本文件方式加載到道頭里;生成層析靜校正量(運(yùn)行Omega地震資料處理系統(tǒng)的Refract1n-Tomo模塊);創(chuàng)建包含光滑CDP靜校正量(浮動處理基準(zhǔn)面)的空間屬性文件。
[0061]注:橫波近地表靜校正(可參考名稱為“震源橫波勘探的近地表靜校正處理方法”的專利)是決定資料處理剖面品質(zhì)的基礎(chǔ),只有在橫波近地表靜校正和地表一致性剩余靜校正問題解決后,去噪問題才能得到較好的解決,進(jìn)而才能提聞速度精度。
[0062]步驟6:縱/橫波信號處理。
[0063]應(yīng)用地表一致性振幅補(bǔ)償;應(yīng)用炮域噪聲衰減;應(yīng)用地表一致性反褶積;應(yīng)用隨機(jī)噪聲衰減。其中,地表一致性振幅補(bǔ)償:運(yùn)行Omega地震資料處理系統(tǒng)的SC_Amplitude_Compensat1n模塊,使得振幅一致性提高;炮域噪聲衰減:運(yùn)行Omega地震資料處理系統(tǒng)的Noise_Atten_Shot模塊,壓制相干噪聲、壓制高頻噪聲、抑制低頻噪聲;地表一致性反裙積:運(yùn)行Omega地震資料處理系統(tǒng)的SC_Decon模塊,起到相位校正、提高分辨率作用;隨機(jī)噪聲衰減運(yùn)行Omega地震資料處理系統(tǒng)的Residual_Noise_Atten模塊,壓制隨機(jī)噪聲。
[0064]步驟7:縱/橫波⑶P道集分選。
[0065]將炮集記錄數(shù)據(jù)(步驟6的輸出),分選成按⑶P道集(運(yùn)行Omega地震資料處理系統(tǒng)的Multrway-Sort模塊)輸出以便后續(xù)處理。
[0066]步驟8:縱/橫波速度分析。
[0067]產(chǎn)生速度譜(運(yùn)行Omega地震資料處理系統(tǒng)的Velocity_Field_Create模塊),從拾取的速度文件生成插值速度場(與剩余靜校正結(jié)合可進(jìn)行多次迭代處理(速度應(yīng)用到剩余靜校正后,一般需再次進(jìn)行速度分析,以提高速度精度。迭代處理這是資料處理中常用的手段,即循環(huán)處理的意思)。
[0068]步驟9:縱/橫波剩余靜校正。
[0069]應(yīng)用剩余靜校正(與速度分析結(jié)合可進(jìn)行多次迭代處理(剩余靜校正應(yīng)用到速度分析后,一般需再次進(jìn)行剩余靜校正,以提高剩余靜校正的精度));生成用于目標(biāo)線偏移的速度(運(yùn)行Omega地震資料處理系統(tǒng)的TL_KPSTM_Velocity_Field_Create模塊)。
[0070]步驟10 -M /橫波疊前去噪。
[0071]對步驟9輸出的X分量、Z分量地震數(shù)據(jù)運(yùn)行疊前去噪、QC疊加(Q代表質(zhì)量,C代表控制,QC疊加是指質(zhì)量控制的疊加)(運(yùn)行Omega地震資料處理系統(tǒng)的Area_Group_3DRNA、Post_3DRNA_Stk_Qc模塊),并從該疊加數(shù)據(jù)體產(chǎn)生切除時(shí)間,以便用于疊后數(shù)據(jù)。
[0072]步驟11:縱/橫波PSTM偏移。
[0073]進(jìn)行目標(biāo)線的偏移(運(yùn)行Omega地震資料處理系統(tǒng)的TL_KPSTM模塊);拾取、生成疊前時(shí)間體偏移的速度場(運(yùn)行Omega地震資料處理系統(tǒng)的Volume_KPSTM_Velocity_Field_Create模塊);用PSTM(疊前時(shí)間偏移)速度體運(yùn)行整塊數(shù)據(jù)體的偏移(運(yùn)行Omega地震資料處理系統(tǒng)的KPSTM模塊);拾取、生成疊前時(shí)間偏移后的疊加速度場(運(yùn)行Omega地震資料處理系統(tǒng)的 Post_KPSTM_Velocity_Field_Create 模塊)。
[0074]步驟12:縱/橫波疊加和修飾處理。
[0075]用疊前時(shí)間偏移后的疊加速度場,產(chǎn)生一個(gè)疊前時(shí)間偏移后的疊加(運(yùn)行Omega地震資料處理系統(tǒng)的KPSTM_Stack模塊);應(yīng)用疊后剩余振幅補(bǔ)償以便平衡振幅(運(yùn)行Omega地震資料處理系統(tǒng)的KPSTM_Stack_TVflt_SRAC模塊));應(yīng)用濾波、去噪處理得到最終信號增強(qiáng)的疊加。
[0076]步驟13:縱/橫波成果數(shù)據(jù)輸出。
[0077]輸出縱/橫波處理成果數(shù)據(jù)體(指處理結(jié)束后用于解釋的疊后縱、橫波數(shù)據(jù),三維的數(shù)據(jù)通常稱為體。是將步驟12處理的成果數(shù)據(jù)以SEGY格式輸出,以便解釋人員使用)。
[0078]圖1a至圖1c為一炮原始三分量炮集記錄,其中圖1a是原始三分量單炮記錄中的X分量,圖1b是原始三分量單炮記錄中的Y分量,圖1c是原始三分量單炮記錄中的Z分量。在Z分量的炮集記錄中,箭頭所指的幾組主要反射波組與X分量中的箭頭所指的幾組主要反射波組相對應(yīng)。其中:
[0079]第一對的反射旅行時(shí)在Z分量記錄中的右邊排列的最遠(yuǎn)道大約為2100毫秒,在X分量記錄中的右邊排列的最遠(yuǎn)道大約為3900毫秒,大約是Z分量反射波的1.86倍。
[0080]第二對的反射旅行時(shí)在Z分量記錄中的右邊排列的最遠(yuǎn)道大約為2350毫秒,在X分量記錄中的右邊排列的最遠(yuǎn)道大約為4350毫秒,大約是Z分量反射波的1.85倍。
[0081]第三對的反射旅行時(shí)在Z分量記錄中的右邊排列的最遠(yuǎn)道大約為2700毫秒,在X分量記錄中的右邊排列的最遠(yuǎn)道大約為5100毫秒,大約是Z分量反射波的1.89倍。
[0082]再計(jì)算層間旅行時(shí)的差異(用于波場性質(zhì)分析,說明X分量所觀測到的地震波是橫波,而不是多波勘探中常用的轉(zhuǎn)換波。),Tp2-Tpi = 250,Ts2-Tsi = 450,層間速度比為1.8 ;Tp3-Tp2 = 350,Ts3-Ts2 = 750,層間速度比為2.14。以上計(jì)算的縱橫波速度比均基本滿足從VSP測井資料得到的縱橫波速度比關(guān)系,如表I所示。
[0083]
[0084]
rD 電纜深校正深縱波層速橫波層速I縱橫波層速I
度m 度m度m/S度m/S度比
12420 24126081__3162__1.9230.315
22430~ 2422 —6257 —3189~ 1.9620.325
3— 2740~2732623332071.944~0J20~
4" 2760~2752641931592.032~0340~
5" 2840—2832593833111.793~0274~
6" 3280—3272630630642.058~0345~
7" 3470—3462726334612.099~0353~
8" 3570—3562654931262.095~0352~
9" 3650—3642635430632.074~0349~
10" 4000—3992599433831.772~0266~
114180 4172538229271.8390.290
[0085]表I
[0086]圖2a和圖2b為Z分量與X分量對比圖,其中,圖2a是利用VSP測到的速度比進(jìn)行拉伸(加靜校正)的Z分量(縱波),圖2b是(加靜校正)的X分量(橫波),其中Z分量利用VSP測到的速度比(淺層為2,深層為1.83)進(jìn)行了拉伸(用于波場性質(zhì)分析,說明X分量所觀測到的地震波是橫波,而不是多波勘探中常用的轉(zhuǎn)換波。),可以看到箭頭對應(yīng)的四組強(qiáng)的反射波組有很好的一致性。另外,X分量在零炮檢距附近,反射能量沒有消失,基本上在零入射角附近反射能量最強(qiáng);前后兩邊排列的反射沒有發(fā)生相位反轉(zhuǎn)的現(xiàn)象,不滿足轉(zhuǎn)換波的波場特征。因此,可以確定這四組波的不是轉(zhuǎn)換波,而是橫波。
[0087]圖3為非零井源距VSP的Z分量記錄??梢栽谟涗浿邪l(fā)現(xiàn)很強(qiáng)的橫波,S波-定義為炸藥震源激發(fā),在自由界面發(fā)生波型轉(zhuǎn)換,產(chǎn)生S波并向下傳播。
[0088]圖4a和圖4b為橫波疊加剖面與縱波疊加拉伸剖面對比,其中,圖4a是X分量疊加剖面,圖4b是Z分量疊加拉伸剖面,對比兩圖可見兩剖面的有效波組具有很好的對應(yīng)關(guān)系O
[0089]本發(fā)明屬地球物理地震勘探領(lǐng)域,具體涉及地震資料的橫波信息的利用及處理,發(fā)明可應(yīng)用于特殊條件下炸藥震源激發(fā)三分量接收的震源橫波的多波地震勘探。
[0090]上述技術(shù)方案只是本發(fā)明的一種實(shí)施方式,對于本領(lǐng)域內(nèi)的技術(shù)人員而言,在本發(fā)明公開了應(yīng)用方法和原理的基礎(chǔ)上,很容易做出各種類型的改進(jìn)或變形,而不僅限于本發(fā)明上述【具體實(shí)施方式】所描述的方法,因此前面描述的方式只是優(yōu)選的,而并不具有限制性的意義。
【權(quán)利要求】
1.一種基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法,其特征在于:所述方法包括: 步驟1:x、Y、Z三分量數(shù)據(jù)輸入:將野外采集的原始三分量炮集數(shù)據(jù)輸入到處理系統(tǒng);步驟2:縱、橫波道頭加載、道編輯:根據(jù)炮點(diǎn)和檢波點(diǎn)信息產(chǎn)生空間屬性文件,將空間屬性文件加載到道頭里;人工進(jìn)行道記錄刪除或修改;生成起始時(shí)間;應(yīng)用網(wǎng)格定義將資料處理的范圍定義處理;輸出加載道頭的炮集原始地震數(shù)據(jù); 步驟3:水平分量偏振旋轉(zhuǎn):利用橫波的初至波對步驟2輸出的X分量和Y分量數(shù)據(jù)進(jìn)行偏振旋轉(zhuǎn); 步驟4:三分量數(shù)據(jù)分選; 步驟5:縱、橫波近地表靜校正; 步驟6:縱、橫波信號處理:應(yīng)用地表一致性振幅補(bǔ)償、炮域噪聲衰減、地表一致性反褶積以及隨機(jī)噪聲衰減進(jìn)行縱、橫波信號處理; 步驟7:縱、橫波CDP道集分選:將步驟6輸出的炮集記錄數(shù)據(jù)分選成按CDP道集輸出;步驟8:縱、橫波速度分析:產(chǎn)生速度譜(包含速度函數(shù)的中間處理成果),然后對速度譜進(jìn)行拾取速度得到速度文件,再從拾取的速度文件生成插值速度場; 步驟9:縱、橫波剩余靜校正:應(yīng)用剩余靜校正,生成用于目標(biāo)線偏移的速度; 步驟10:縱、橫波疊前去噪:對步驟9輸出的X分量和Z分量地震數(shù)據(jù)運(yùn)行疊前去噪、QC疊加,并從所得的疊加數(shù)據(jù)體產(chǎn)生切除時(shí)間; 步驟11:縱、橫波PSTM偏移:進(jìn)行目標(biāo)線的偏移;拾取、生成疊前時(shí)間體偏移的速度場;用PSTM速度體運(yùn)行整塊數(shù)據(jù)體的偏移;拾取、生成疊前時(shí)間偏移后的疊加速度場; 步驟12:縱、橫波疊加和修飾處理:用疊前時(shí)間偏移后的疊加速度場,產(chǎn)生一個(gè)疊前時(shí)間偏移后的疊加;應(yīng)用疊后剩余振幅補(bǔ)償以便平衡振幅;應(yīng)用濾波、去噪處理得到最終信號增強(qiáng)的疊加; 步驟13:縱、橫波成果數(shù)據(jù)輸出:將步驟12得到的成果數(shù)據(jù)以SEGY格式輸出。
2.根據(jù)權(quán)利要求1所述的基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法,其特征在于:所述步驟I中,如果輸入的原始三分量炮集數(shù)據(jù)是按多路方式記錄的格式進(jìn)行記錄的,則首先對數(shù)據(jù)進(jìn)行解編,即對一個(gè)大的矩陣進(jìn)行變換,使變換后的矩陣的行能按地震道讀數(shù)據(jù)。
3.根據(jù)權(quán)利要求1所述的基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法,其特征在于:所述步驟2中道記錄刪除或修改是指:刪除不正常的道記錄和修改反極性的道記錄。
4.根據(jù)權(quán)利要求1所述的基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法,其特征在于:所述步驟4是這樣實(shí)現(xiàn)的:將步驟3輸出的混在一起的三分量地震記錄數(shù)據(jù)分選成X、Y、Z三個(gè)分量,分別輸出。
5.根據(jù)權(quán)利要求1所述的基于特殊條件下炸藥震源激發(fā)三分量接收的橫波處理方法,其特征在于:所述步驟5是這樣實(shí)現(xiàn)的:從步驟4輸出的三個(gè)分量記錄上自動拾取縱、橫波的初至?xí)r間;對拾取時(shí)間表中自動拾取的初至?xí)r間進(jìn)行復(fù)查,并通過文本文件方式加載到道頭里;生成層析靜校正量;創(chuàng)建包含光滑CDP靜校正量的空間屬性文件。
【文檔編號】G01V1/30GK104166157SQ201310181853
【公開日】2014年11月26日 申請日期:2013年5月16日 優(yōu)先權(quán)日:2013年5月16日
【發(fā)明者】鄔達(dá)理, 黃忠玉 申請人:中國石油化工股份有限公司, 中國石油化工股份有限公司石油物探技術(shù)研究院