專利名稱::二維疊加速度及均方根速度場(chǎng)閉合差校正方法
技術(shù)領(lǐng)域:
:本發(fā)明涉及一種對(duì)低幅度構(gòu)造圈閉石油地震勘探二維地震工區(qū)速度建場(chǎng)精度的疊加速度及均方根速度場(chǎng)校正方法。
背景技術(shù):
:在二維地震工區(qū)的資料解釋、變速成圖工作中,測(cè)線交點(diǎn)普遍存在反射層時(shí)間和速度場(chǎng)閉合差問(wèn)題。目前,針對(duì)測(cè)線交點(diǎn)反射層時(shí)間閉合差問(wèn)題的研究比較深入,但對(duì)速度場(chǎng)閉合差問(wèn)題的研究還不多見(jiàn)。比較有代表性的有左龍憑等人在2000年提出的對(duì)不同測(cè)線計(jì)算校正算子,對(duì)交點(diǎn)閉合差作時(shí)變和空變時(shí)間差校正及波形校正。另外,陳遵德等人提出采用多段折線最優(yōu)逼近平均速度的方法,把速度曲線分成數(shù)段,每段對(duì)應(yīng)一套地層,在每段內(nèi)求出相應(yīng)的時(shí)間閉合差,采用數(shù)學(xué)方法進(jìn)行校正,從而解決相交測(cè)線速度閉合差問(wèn)題。上述兩種方法都在不同程度上消除了交叉測(cè)線的閉合差問(wèn)題,提高了二維工區(qū)的速度建場(chǎng)精度。但都沒(méi)有從閉合差產(chǎn)生的根本原因入手來(lái)解決問(wèn)題,不可避免地存在著局限性①對(duì)各種原因引起的速度閉合差實(shí)行一刀切,做一次性校正,精度必然受到制約;②均采用數(shù)學(xué)曲線擬合的方法,前者用一個(gè)擬合函數(shù),后者是多段擬合,降低了疊加(平均)速度的精度,同時(shí)也降低了閉合差求取的精度;③具體方法存在個(gè)性缺陷,前者針對(duì)疊加速度做校正,只適合于小地層傾角;后者在把速度縱向分段擬合時(shí),方法本身就可能帶來(lái)速度閉合問(wèn)題,并且該方法只能針對(duì)大套地層,精度受到限制。交叉測(cè)線速度不閉合的問(wèn)題是制約二維地震工區(qū),尤其是低幅度構(gòu)造二維地震工區(qū)速度建場(chǎng)精度的關(guān)鍵因素之一,因此探究交叉測(cè)線速度閉合差產(chǎn)生的根本原因,研發(fā)科學(xué)有效的校正技術(shù)成為當(dāng)前二維地震工區(qū)速度研究的主要內(nèi)容。前人的研究沒(méi)有從速度閉合差產(chǎn)生的根本原因入手,試圖采用純粹的數(shù)學(xué)方法徹底解決地球物理問(wèn)題,不可避免地存在各種缺陷。
發(fā)明內(nèi)容本發(fā)明的目的是根據(jù)形成二維速度場(chǎng)閉合差的根本原因,即地震資料主測(cè)線和聯(lián)絡(luò)測(cè)線不閉合,資料解釋時(shí)的T。閉合差導(dǎo)致疊加速度場(chǎng)主測(cè)線、聯(lián)絡(luò)測(cè)線的閉合差以及在將疊加速度轉(zhuǎn)換為均方根速度時(shí),由于主測(cè)線、聯(lián)絡(luò)測(cè)線視傾角的差別,均方根速度本身存在閉合差這兩個(gè)方面采用兩步法分別對(duì)這兩種閉合差進(jìn)行校正,最后得到精確的二維地震速度場(chǎng)。1、解釋T。閉合差引起的疊加速度閉合差校正技術(shù)。由于T。時(shí)間和疊加速度譜是一一對(duì)應(yīng)的關(guān)系,而在二維剖面解釋過(guò)程中,往往存在時(shí)間閉合差,在編制T。圖時(shí)要對(duì)T。時(shí)間閉合差進(jìn)行校正,同樣這個(gè)時(shí)間閉合差也會(huì)引起速度譜的閉合差,因此要對(duì)這個(gè)速度閉合差進(jìn)行校正。校正時(shí)根據(jù)解釋提供的每條線的T。時(shí)間和T。平面圖,計(jì)算每條線T。時(shí)間與T。平面圖的時(shí)間差A(yù)T(見(jiàn)附圖1),把這個(gè)時(shí)間差作為速度譜的校正量,根據(jù)AT正負(fù)(上抬為負(fù),下移為正)的不同,用不同的速度和校正公式對(duì)速度譜做校正當(dāng)AT>0時(shí),校正速度用地表充填速度v。第一步用Dix公式由速度譜上的疊加速度計(jì)算第n層的層速度<formula>formulaseeoriginaldocumentpage5</formula>第二步由層速度反算疊加速度當(dāng)AT<0時(shí),校正速度用速度譜上的速度。第一步用Dix公式由速度譜上的疊加速度計(jì)算第n層的層速度。第二步由下面公式由層速度反算疊加速度<formula>formulaseeoriginaldocumentpage5</formula>2、主測(cè)線、聯(lián)絡(luò)線視傾角不同引起的均方根速度閉合差校正技術(shù)。對(duì)由于主測(cè)線和聯(lián)絡(luò)測(cè)線視傾角的不同形成的閉合差我們按如下流程進(jìn)行校正(見(jiàn)附圖2):①?gòu)闹鳒y(cè)線、聯(lián)絡(luò)測(cè)線疊加剖面上,分別拾取出剖面的時(shí)間斜率4、d2,②利用地震測(cè)線分出主測(cè)線、聯(lián)絡(luò)測(cè)線的疊加速度Val、Va2;③根據(jù)如下公式計(jì)算主測(cè)線、聯(lián)絡(luò)測(cè)線的視傾角9i和^、均方根速度V、Vk<formula>formulaseeoriginaldocumentpage5</formula>④計(jì)算主測(cè)線與聯(lián)絡(luò)測(cè)線交點(diǎn)處的加權(quán)系數(shù)Cl(i,j)、c2(i,j)<formula>formulaseeoriginaldocumentpage5</formula>⑤對(duì)交點(diǎn)處的加權(quán)系數(shù)q(i,j)、c2(i,j)各自進(jìn)行插值,得到各自加權(quán)系數(shù)散點(diǎn)Ci(x,y)c2(x,y)。⑥歸一化處理最終各點(diǎn)的加權(quán)系數(shù)為c2(x,>Q__qO,力+c20,力qo,力+c2o,力⑦對(duì)主測(cè)線和聯(lián)絡(luò)線均方根速度插值得到各自均方根速度散點(diǎn)vK1(x,y)VK2(x,y)⑧加權(quán)疊加得到最終的均方根速度場(chǎng)<formula>formulaseeoriginaldocumentpage6</formula>在校正主測(cè)線和聯(lián)絡(luò)線的速度閉合差時(shí),關(guān)鍵的難點(diǎn)是剖面時(shí)間斜率的求取,在實(shí)際研究中,利用以下兩種方式來(lái)獲得時(shí)間斜率一是從segy時(shí)間剖面獲得。通過(guò)對(duì)segy時(shí)間剖面的同相軸自動(dòng)追蹤,對(duì)追蹤的曲線求導(dǎo),獲得各剖面沿測(cè)線方向的時(shí)間斜率,然后與測(cè)量成果進(jìn)行坐標(biāo)匹配,獲得空間位置上的時(shí)間斜率。二是通過(guò)T。時(shí)間圖獲得。首先把測(cè)線分成主測(cè)線和聯(lián)絡(luò)測(cè)線兩個(gè)方向,然后沿測(cè)線分別到各層T。圖上求時(shí)間斜率,并對(duì)時(shí)間斜率做空間網(wǎng)格化,插值獲得空間每個(gè)點(diǎn)的時(shí)間斜率。該技術(shù)實(shí)施過(guò)程中要首先進(jìn)行第一步校正,然后再進(jìn)行第二步校正,順序固定,不能顛倒,經(jīng)過(guò)兩步速度閉合差的校正,就可以得到相對(duì)精確的速度場(chǎng)。發(fā)明效果本次研究從地球物理學(xué)的角度出發(fā),探究二維地震速度場(chǎng)閉合差產(chǎn)生的根本原因,對(duì)各種原因引起的速度閉合差進(jìn)行深入的分析研究,結(jié)合構(gòu)造地質(zhì)因素,利用精確的數(shù)學(xué)公式對(duì)不同的速度閉合差進(jìn)行計(jì)算及分步校正,從而科學(xué)有效地消除交叉測(cè)線速度閉合差,切實(shí)提高二維地震工區(qū)速度場(chǎng)精度。該技術(shù)與前人研究成果對(duì)比有以下幾個(gè)優(yōu)點(diǎn)①對(duì)速度閉合差進(jìn)行深入分析研究,根據(jù)其形成的原因分別進(jìn)行計(jì)算并分步進(jìn)行校正,比前人對(duì)速度閉合差實(shí)施"一刀切"的不科學(xué)的校正方法邁進(jìn)了一大步;②用嚴(yán)格的數(shù)學(xué)公式計(jì)算來(lái)求取速度閉合差,這比前人用曲線擬合方式來(lái)求取速度閉合差要精確得多。③可以對(duì)任何目的層速度閉合差進(jìn)行計(jì)算和校正,不受地層厚薄、傾角大小的限制,提高了速度研究精度。2006年以來(lái),利用該技術(shù)成果在塔里木盆地低幅度構(gòu)造區(qū)進(jìn)行了重點(diǎn)研究,在下第三系及志留系共發(fā)現(xiàn)低幅度構(gòu)造圈閉6個(gè),總面積計(jì)55kn^,上鉆6口探井,其中1口井見(jiàn)油氣顯示,5口井獲得工業(yè)油氣流,預(yù)測(cè)資源量天然氣1500億立方,凝析油3500萬(wàn)噸。圖1T。閉合差及其引起的速度閉合差校正示意圖圖2均方根速度閉合差校正流程圖圖3某工區(qū)地震T8反射層平均速度圖(原始)圖4某工區(qū)地震T8反射層平均速度圖(T。閉合差引起的速度閉合差校正)圖5某工區(qū)地震T8反射層平均速度圖(均方根速度閉合差校正)圖6B工區(qū)某反射層平均速度圖(原始)圖7B工區(qū)某反射層平均速度圖(T。閉合差引起的速度閉合差校正)圖8B工區(qū)某反射層平均速度圖(均方根速度閉合差校正)具體實(shí)施方式實(shí)施例1某二維地震工區(qū)勘探目的層埋藏深,構(gòu)造起伏比較平緩,對(duì)地震速度場(chǎng)精度要求比較高。工區(qū)原始地震資料存在較大的速度閉合差(圖3),需要對(duì)其進(jìn)行校正,以提高地震速度場(chǎng)的精度。1、T。閉合差引起的疊加速度閉合差校正①T。閉合差的求取首先把整個(gè)工區(qū)的地震解釋數(shù)據(jù)經(jīng)T。閉合差校正之后進(jìn)行網(wǎng)格插值,形成整個(gè)工.......r。圖的曲面文件減去每條地震測(cè)線沒(méi)有做T。閉合差校正的散點(diǎn)r的at執(zhí)占救坦區(qū)的地震解釋等T。圖,用T。,,數(shù)據(jù),得到沿地震測(cè)線分布的AT散點(diǎn)數(shù)據(jù)。②T。閉合差引起的疊加速度閉合差校正用AT作為疊加速度譜的校正量,根據(jù)AT符號(hào)的正負(fù)(上抬為負(fù),下移為正),用不同的速度及相應(yīng)的校正公式對(duì)速度譜做校正。",當(dāng)AT>0時(shí),校正速度用地表充填速度v。第一步用Dix公式由速度譜上的疊加速度計(jì)算第n層的層速度V—fV/—V,/,第二步由層速度反算疊加速度1v02A/+S《MA/1+2A/,當(dāng)AT<0時(shí),校正速度用速度譜上的速度。第一步用Dix公式由速度譜上的疊加速度計(jì)算第n層的層速度,第二步由下面公式由層速度反算疊加速度i'H《My1+a/n"1^^tij-ti(w)^Ajt工區(qū)速度閉合差基本在70m/s以內(nèi),經(jīng)過(guò)上述對(duì)由T。閉合差引起的疊加速度閉合差的校正,大部分速度閉合差已經(jīng)消除(圖4)。2、主測(cè)線、聯(lián)絡(luò)線視傾角不同引起的均方根速度閉合差校正工區(qū)地層變化相對(duì)平緩,地層傾角不大,因此由主測(cè)線、聯(lián)絡(luò)線視傾角不同引起的均方根速度閉合差相對(duì)較小。①利用T。圖求空間位置上的時(shí)間斜率。工區(qū)面積大、地震測(cè)線多,用地震剖面求取時(shí)間斜率工作量太大,因此用T。圖求空間位置上的時(shí)間斜率。首先由工區(qū)地震解釋成果獲得地震解釋等T。圖,然后沿主測(cè)線和聯(lián)絡(luò)測(cè)線兩個(gè)方向分別到各層T。圖上求時(shí)間斜率,并對(duì)時(shí)間斜率做空間網(wǎng)格插值、平滑獲得空間每個(gè)點(diǎn)的時(shí)間斜率。由于工區(qū)地層傾角比較小,時(shí)間斜率變化平緩,因此該工區(qū)做時(shí)間斜率的空間網(wǎng)格插值時(shí),插值半徑橫向取100m(取值范圍一般在0-1000m)、縱向取50ms(取值范圍一般在0-200ms),平滑系數(shù)取2(取值范圍一般在0-30);②利用工區(qū)地震測(cè)線分出主測(cè)線、聯(lián)絡(luò)測(cè)線的疊加速度Val、Va2;③根據(jù)如下公式計(jì)算主測(cè)線、聯(lián)絡(luò)測(cè)線的視傾角e工和e2、均方根速度vK1、vK2;7|1+-2J&2=ll+lK2「《=咖_1.1+《=cos-④計(jì)算主測(cè)線與聯(lián)絡(luò)測(cè)線交點(diǎn)處的加權(quán)系數(shù)Cl(i,j)、C2(i,j)sin2+cos26>2(/,_/)c2(',_/)_cos2+sin22…2⑤對(duì)交點(diǎn)處的加權(quán)系數(shù)Cl(i,j)、c2(i,j)各自進(jìn)行插值,得到工區(qū)均勻分布的加權(quán)系數(shù)散點(diǎn)q(x,y)c2(x,y)。⑥對(duì)加權(quán)系數(shù)散點(diǎn)做歸一化處理,得到平滑的、均勻變化的加權(quán)系數(shù)(x,力+c2O,力(X,力+c2(JC,力⑦對(duì)主測(cè)線和聯(lián)絡(luò)線均方根速度插值得到各自均方根速度散點(diǎn)vK1(x,y)VK2(x,y)⑧用⑥得到的加權(quán)系數(shù)對(duì)由⑦得到的均方根速度散點(diǎn)進(jìn)行加權(quán)疊加,也就是做均方根速度傾角校正得到最終消除閉合差的均方根速度場(chǎng),、c20,_y)Vff(X,>0=經(jīng)過(guò)主、聯(lián)測(cè)線視傾角不同引起的均方根速度閉合差校正,計(jì)算得出的平均速度閉合差已經(jīng)完全消除(圖5),在經(jīng)過(guò)上述兩步法速度閉合差校正之后,有效提高了二維速度場(chǎng)的精度。用該速度場(chǎng)對(duì)工區(qū)地震解釋成果進(jìn)行時(shí)深轉(zhuǎn)換,新發(fā)現(xiàn)了3個(gè)低幅度構(gòu)造圈閉,其中1號(hào)圈閉已經(jīng)鉆探成功,預(yù)測(cè)資源量天然氣400億立方,凝析油600萬(wàn)噸。實(shí)施例2B工區(qū)地表地下地質(zhì)結(jié)構(gòu)復(fù)雜,二維地震測(cè)網(wǎng)為2X44X4km,存在施工年度不同,處理單位不同等因素,交叉測(cè)線速度閉合問(wèn)題比較嚴(yán)重(圖6),變速成圖的重點(diǎn)是做好閉合差校正工作。采用如下校正步驟1、T。閉合差引起的疊加速度閉合差校正該區(qū)T。閉合差20-45ms,其引起的速度閉合差不容忽視,必須進(jìn)行校正。①T。閉合差的求取首先把整個(gè)工區(qū)的地震解釋數(shù)據(jù)經(jīng)T。閉合差校正之后進(jìn)行網(wǎng)格插值,形成整個(gè)工區(qū)的地震解釋等T。圖,用T。圖的曲面文件減去每條地震測(cè)線沒(méi)有做T。閉合差校正的散點(diǎn)數(shù)據(jù),得到沿地震測(cè)線分布的AT散點(diǎn)數(shù)據(jù)。②T。閉合差引起的疊加速度閉合差校正用AT作為疊加速度譜的校正量,根據(jù)AT符號(hào)的正負(fù)(上抬為負(fù),下移為正),用不同的速度及相應(yīng)的校正公式對(duì)速度譜做校正。當(dāng)AT>0時(shí),校正速度用地表充填速度v。第一步用Dix公式由速度譜上的疊加速度計(jì)算第n層的層速度第二步由層速度反算疊加速度v02A/+Sv》,當(dāng)AT<0時(shí),校正速度用速度譜上的速度。第一步用Dix公式由速度譜上的疊加速度計(jì)算第n層的層速度c第二步由下面公式由層速度反算疊加速度v"=2M從圖6我們可以看出,校正前速度閉合差達(dá)到40m/s,構(gòu)造埋深時(shí)間深度大于3s時(shí),深度誤差將超過(guò)60m,無(wú)法落實(shí)幅度小于100m的圈閉。經(jīng)過(guò)T。閉合差引起的疊加速度閉合差校正后所計(jì)算出的平均速度(圖7)測(cè)線交點(diǎn)閉合差已基本小于15m/s。2、主測(cè)線、聯(lián)絡(luò)線視傾角不同引起的均方根速度閉合差校正①利用地震segy時(shí)間剖面求空間位置上的時(shí)間斜率。工區(qū)面積不大、地震測(cè)線相對(duì)較少,用地震剖面求取時(shí)間斜率以提高時(shí)間斜率的求取精度。由于工區(qū)地層傾角相對(duì)較大,時(shí)間斜率變化較大,因此該工區(qū)做時(shí)間斜率的空間網(wǎng)格插值時(shí),為獲得相對(duì)平滑的趨勢(shì)面,插值半徑橫向取值稍大為300m(取值范圍一般在0-1000m)、縱向取80ms(取值范圍一般在0-200ms),平滑系數(shù)取5(取值范圍一般在0-30);②利用工區(qū)地震測(cè)線分出主測(cè)線、聯(lián)絡(luò)測(cè)線的疊加速度Val、Va2;③根據(jù)如下公式計(jì)算主測(cè)線、聯(lián)絡(luò)測(cè)線的視傾角9i和92、均方根速度|1+&2=-K2力+'M丫…",洲—④計(jì)算主測(cè)線與聯(lián)絡(luò)測(cè)線交點(diǎn)處的加權(quán)系數(shù)Cl(i,j)、c2(i,j)《=cos—'sin2《(/,/)+cos2c2(/,/)=方根速度傾角校正得到最終消除閉合差的均方根速度場(chǎng)cos2《(/,力+sin2P2(/,力2…2⑤對(duì)交點(diǎn)處的加權(quán)系數(shù)Cl(i,j)、c2(i,j)各自進(jìn)行插值,得到工區(qū)均勻分布的加權(quán)系數(shù)散點(diǎn)q(x,y)c2(x,y)。⑥對(duì)加權(quán)系數(shù)散點(diǎn)做歸一化處理,得到平滑的、均勻變化的加權(quán)系數(shù)c!0,力+c20,力c,0,力+c2(x,力⑦對(duì)主測(cè)線和聯(lián)絡(luò)線均方根速度插值得到各自均方根速度散點(diǎn)VK1(x,氛(x,y)⑧用⑥得到的加權(quán)系數(shù)對(duì)由⑦得到的均方根速度散點(diǎn)進(jìn)行加權(quán);t加,也就是做均9<formula>formulaseeoriginaldocumentpage10</formula>經(jīng)過(guò)主、聯(lián)測(cè)線視傾角不同引起的均方根速度閉合差校正,計(jì)算得出的平均速度閉合差已經(jīng)完全消除(圖8),同時(shí)保持了原始速度的整體規(guī)律,速度橫向變化呈現(xiàn)北低南高,中部低速條帶的特征。這與工區(qū)北高(地層尖滅)南低,中部沿?cái)嗔严到y(tǒng)分布低幅度構(gòu)造帶的構(gòu)造地質(zhì)規(guī)律是一致的。權(quán)利要求一種二維疊加速度及均方根速度場(chǎng)校正方法,其特征在于(1)、T0閉合差引起的疊加速度閉合差校正①T0閉合差的求取首先對(duì)整個(gè)二維工區(qū)的地震解釋數(shù)據(jù)做T0閉合差校正,然后網(wǎng)格插值,形成整個(gè)工區(qū)的地震解釋等T0圖,用T0圖的曲面文件減去每條地震測(cè)線原始的沒(méi)有做T0閉合差校正的T0散點(diǎn)數(shù)據(jù),得到沿地震測(cè)線分布的ΔT散點(diǎn)數(shù)據(jù),當(dāng)T0散點(diǎn)數(shù)據(jù)相對(duì)于T0圖上抬時(shí)ΔT為負(fù),反之ΔT為正。②T0閉合差引起的疊加速度閉合差校正用ΔT作為疊加速度譜的校正量,根據(jù)ΔT符號(hào)的正負(fù),用不同的速度及不同的校正公式對(duì)速度譜做校正。當(dāng)ΔT>0時(shí),校正速度用地表充填速度v0第一步用Dix公式由速度譜上的疊加速度計(jì)算第n層的層速度<mrow><msub><mi>v</mi><mi>cn</mi></msub><mo>=</mo><msqrt><mfrac><mrow><msubsup><mi>v</mi><mi>n</mi><mn>2</mn></msubsup><msub><mi>t</mi><mi>n</mi></msub><mo>-</mo><msubsup><mi>v</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></msubsup><msub><mi>t</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow><mrow><msub><mi>t</mi><mi>n</mi></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mfrac></msqrt></mrow>第二步由層速度反算疊加速度<mrow><msub><mi>v</mi><mi>n</mi></msub><mo>=</mo><msqrt><mfrac><mrow><msubsup><mi>v</mi><mn>0</mn><mn>2</mn></msubsup><mi>Δt</mi><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>v</mi><mi>ci</mi><mn>2</mn></msubsup><msub><mi>Δt</mi><mi>i</mi></msub></mrow><mrow><mi>Δt</mi><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>Δt</mi><mi>i</mi></msub></mrow></mfrac></msqrt></mrow>當(dāng)ΔT<0時(shí),校正速度用速度譜上的速度;第一步用Dix公式由速度譜上的疊加速度計(jì)算第n層的層速度。第二步由下面公式由層速度反算疊加速度<mrow><msub><mi>v</mi><mi>n</mi></msub><mo>=</mo><msqrt><mfrac><mrow><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi></mrow><mi>n</mi></munderover><msubsup><mi>v</mi><mi>ci</mi><mn>2</mn></msubsup><msub><mi>Δt</mi><mi>i</mi></msub></mrow><mrow><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>Δt</mi><mi>i</mi></msub><mo>+</mo><mi>Δt</mi></mrow></mfrac></msqrt><mfencedopen=''close=''><mtable><mtr><mtd><msub><mi>t</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow></msub><mo><</mo><mi>Δt</mi><mo><</mo><msub><mi>t</mi><mi>k</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Δt</mi><mi>k</mi></msub><mo>=</mo><msub><mi>t</mi><mi>k</mi></msub><mo>-</mo><mi>Δt</mi></mtd></mtr></mtable></mfenced></mrow>(2)主測(cè)線、聯(lián)絡(luò)線視傾角不同引起的均方根速度閉合差校正①?gòu)闹鳒y(cè)線、聯(lián)絡(luò)測(cè)線疊加剖面上,分別拾取出剖面的時(shí)間斜率d1、d2,②利用地震測(cè)線分出主測(cè)線、聯(lián)絡(luò)測(cè)線的疊加速度Va1、Va2;③根據(jù)如下公式計(jì)算主測(cè)線、聯(lián)絡(luò)測(cè)線的視傾角θ1和θ2、均方根速度VR1、VR2;<mrow><msub><mi>V</mi><mrow><mi>R</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>V</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><msqrt><mn>1</mn><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><msub><mi>d</mi><mn>1</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>,</mo></mrow><mrow><msub><mi>V</mi><mrow><mi>R</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>V</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><msqrt><mn>1</mn><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><msub><mi>d</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac></mrow><mrow><msub><mi>θ</mi><mn>1</mn></msub><mo>=</mo><msup><mi>cos</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mfrac><mn>1</mn><msqrt><mn>1</mn><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><msub><mi>d</mi><mn>1</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>]</mo><mo>,</mo></mrow><mrow><msub><mi>θ</mi><mn>2</mn></msub><mo>=</mo><msup><mi>cos</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mfrac><mn>1</mn><msqrt><mn>1</mn><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><msub><mi>d</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>]</mo></mrow>④計(jì)算主測(cè)線與聯(lián)絡(luò)測(cè)線交點(diǎn)處的加權(quán)系數(shù)c1(i,j)、c2(i,j)<mrow><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>sin</mi><mn>2</mn></msup><msub><mi>θ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>cos</mi><mn>2</mn></msup><msub><mi>θ</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac></mrow><mrow><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>cos</mi><mn>2</mn></msup><msub><mi>θ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>sin</mi><mn>2</mn></msup><msub><mi>θ</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac></mrow>⑤對(duì)交點(diǎn)處的加權(quán)系數(shù)c1(i,j)、c2(i,j)各自進(jìn)行插值,得到各自加權(quán)系數(shù)散點(diǎn)c1(x,y)c2(x,y)。⑥歸一化處理最終各點(diǎn)的加權(quán)系數(shù)為<mfencedopen=''close=''><mtable><mtr><mtd><mfrac><mrow><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac></mtd><mtd><mfrac><mrow><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced>⑦對(duì)主測(cè)線和聯(lián)絡(luò)線均方根速度插值得到各自均方根速度散點(diǎn)vR1(x,y)vR2(x,y)⑧加權(quán)疊加得到最終的均方根速度場(chǎng)<mrow><msub><mi>v</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><msub><mi>v</mi><mrow><mi>R</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><msub><mi>v</mi><mrow><mi>R</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>.</mo></mrow>F2008102394638C0000022.tif全文摘要本發(fā)明涉及一種二維疊加速度及均方根速度場(chǎng)的兩步法校正方法;首先求取T0曲面與T0閉合校正前的每條測(cè)線的時(shí)間差ΔT,利用ΔT對(duì)疊加速度做基準(zhǔn)面校正,ΔT>0時(shí),校正速度用地表充填速度,ΔT<0時(shí),校正速度用速度譜上的速度;針對(duì)在疊加速度轉(zhuǎn)換為均方根速度時(shí)由視傾角引起的均方根速度閉合差,首先獲取主測(cè)線、聯(lián)絡(luò)測(cè)線疊加剖面的時(shí)間斜率,求取不同測(cè)線方向的視傾角和均方根速度,計(jì)算交點(diǎn)處的加權(quán)系數(shù),對(duì)交點(diǎn)處的加權(quán)系數(shù)插值并歸一化處理,利用不同測(cè)線方向的加權(quán)系數(shù)對(duì)主測(cè)線和聯(lián)絡(luò)線均方根速度加權(quán)疊加得到最終的均方根速度場(chǎng);本方法科學(xué)有效地消除了交叉測(cè)線速度閉合差,提高了二維地震工區(qū)速度場(chǎng)精度。文檔編號(hào)G01V1/36GK101750628SQ20081023946公開(kāi)日2010年6月23日申請(qǐng)日期2008年12月11日優(yōu)先權(quán)日2008年12月11日發(fā)明者張耀堂,彭更新,楊舉勇,楊金華,滿益志,王興軍,肖又軍,賴敬容申請(qǐng)人:中國(guó)石油天然氣股份有限公司