專利名稱：從已獲得的壓縮波反射信息估算剪切波反射信息的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種從已獲得的壓縮波反射信息估算剪切波反射信息的方法。
在地球物理學(xué)研究中，壓縮波或P波反射特性的用途是眾所周知的。各種地震剖面一般都是從可提供關(guān)于地層特性有用信息的各種壓縮波地下反射特性產(chǎn)生的。近代，在獲得能更多地提供有關(guān)各種地下地層特性有用信息的各種地震剪切波或S波地震剖面方面，人們的關(guān)注日益增長(zhǎng)，因?yàn)檫@些有用信息能與從壓縮波地震剖面得到的信息一起使用。例如，壓縮波地震剖面能提供關(guān)于地下各地層的可壓縮性方面的有用信息，而剪切波地震剖面能提供關(guān)于地下地層剛度方面的有用信息。剪切波地震反射波的噪聲比壓縮波地震反射波的噪聲大，因而使判斷困難。此外對(duì)剪切波的反射波的檢測(cè)比壓縮波的反射波的檢測(cè)更為困難。剪切波反射波的幅值一般比壓縮波的反射波小得多，因而使檢測(cè)困難;并且為了使剪切波面直接在地下地層中傳播以感生出一種剪切波反射波，同為了獲得壓縮波反射特性的要求相比，則前者需要專用傳感器和更多步驟。這就使剪切波反射波特性的獲得更加困難。成本更高而且耗時(shí)更多。
剪切波特性是難以獲得的，然而它在勘探烴類物質(zhì)方面卻很有用處。產(chǎn)生壓縮波反射波幅度反?，F(xiàn)象的烴類礦層并不產(chǎn)生類似的剪切波反射波幅度異?，F(xiàn)象。出現(xiàn)這樣一種結(jié)果是因?yàn)榧羟胁ú⒉豁憫?yīng)于任何流體，因此，對(duì)氣體、油類和水并不產(chǎn)生不同的幅度異?，F(xiàn)象?？墒?，異常巖性卻會(huì)引起壓縮波幅度的異?，F(xiàn)象，例如鹽、煤和硬巖夾層也通常都具有相同的異常剪切波特性。然而，因?yàn)榇蠖鄶?shù)所研究的幅度異?，F(xiàn)象是出現(xiàn)在近海的，而剪切波地震數(shù)據(jù)只能在岸上錄出，所以剪切波信息的這樣一種用途并沒(méi)得到廣泛利用。
若干年來(lái)一直都在進(jìn)行著要把剪切波信息應(yīng)用于提供有用數(shù)據(jù)的努力。由S.Domenico著并于1984年刊載在“Geophysic(地球物理學(xué))”雜志49卷第8期上的題為<根據(jù)剪切壓縮波速度確定巖性和孔隙度>一文和題為<利用剪切、壓縮波地震特性解釋地層的情況研究>(由M.McCormik等人著、在1984年“Geophysics”雜志49卷第5期上刊載的)一文都集中地以Vp/Vs比值作為巖性和孔隙度的一個(gè)有用特征。但是在許多情況下，Vp和Vs是密切相關(guān)的。因此為了從剪切波特性抽取關(guān)于地層透視方面的有用信息，就有必要在相當(dāng)細(xì)的間隔中測(cè)量Vp/Vs比值的較小變化。實(shí)踐證明，這樣一種可能性用現(xiàn)有的剪切波技術(shù)是難以實(shí)現(xiàn)的。
然而，存在有一種Vs與Vp具有截然不同特性的情況，這種情況出現(xiàn)于因在孔隙中存在有氣體而強(qiáng)烈影響巖石特性的地球下表面部分中。當(dāng)某一多孔巖石的孔內(nèi)流體從海水變?yōu)闅怏w而沿著某一斷層上升時(shí)，則Vp一般是減低的。在高度多孔的弱固結(jié)巖石中，Vp將大為減低，并產(chǎn)生反射比和烴類物質(zhì)指示標(biāo)記的明顯橫向變化。相反，所述孔隙中的流體變化只引起Vs的小量變化，結(jié)果使反射比的橫向變化可忽略不計(jì)。因此，當(dāng)烴類物質(zhì)指示標(biāo)記應(yīng)在壓縮波地震剖面上突出為幅度異常時(shí)，則在相應(yīng)的剪切波地震剖面上即使有幅度異常也極小。如果壓縮波幅度異常是由諸如鹽或煤的巖性異常所引起的話，則在剪切波剖面上也同樣應(yīng)顯現(xiàn)。因此，剪切波的最好的可能用途應(yīng)當(dāng)是用來(lái)辨別真實(shí)的(烴類物質(zhì))和假的(巖性)幅度異常現(xiàn)象。
人們久已認(rèn)識(shí)到壓縮波的反射率是入射角的一個(gè)函數(shù)。Zoeppritz公式描述了這種性能，該公式包含剪切波速度信息。雖然剪切波數(shù)據(jù)可通過(guò)從幅度水平錯(cuò)距數(shù)據(jù)提取這樣的剪切波數(shù)據(jù)而根據(jù)傳統(tǒng)的數(shù)據(jù)估算出來(lái)，但是有許多問(wèn)題使它難于從真實(shí)的地震數(shù)據(jù)獲得剪切波反射率的有效估算值。例如，這種剪切波反射估算值(下文稱為偽剪切波估算)對(duì)于NMO速度方面的小量誤差是極為敏感的。各偽剪切波值都是通過(guò)在每個(gè)層位的某一NMO校正后的CDP集合中測(cè)量P波數(shù)據(jù)的幅度特性而得到的。為了防止引入各種速度誤差而對(duì)時(shí)間或速度所進(jìn)行的校正是使各軌跡之間產(chǎn)生近于完全的時(shí)間校正，對(duì)該CDP集合中的每個(gè)反射量來(lái)說(shuō)是需要的。但是，盡管在進(jìn)行NMO校正時(shí)、對(duì)于堆垛速度的確定給予了注意，常用的NMO校正方法還總趨于遺留某些地震數(shù)據(jù)的殘留標(biāo)準(zhǔn)時(shí)差以致對(duì)后來(lái)進(jìn)行的各剪切估算值引入各種誤差。因此本發(fā)明的一個(gè)目的是要避免或減輕這個(gè)問(wèn)題。
因此，本發(fā)明歸屬于一種從壓縮波反射信息產(chǎn)生剪切波反射信息的方法，它包括下列各步驟a)從若干聲波源產(chǎn)生逐漸下降的聲波，并在若干接收機(jī)處接收和記錄在某一地下界面上產(chǎn)生的逐漸上升的波形;
b)把所述錄下的反射波集合成具有某一公共反射邊界的各不同源和接收機(jī)的水平錯(cuò)距組;
c)對(duì)所述集合的反射波進(jìn)行規(guī)定的標(biāo)準(zhǔn)時(shí)差校正;
d)對(duì)所述標(biāo)準(zhǔn)時(shí)差校正后的波進(jìn)行殘留標(biāo)準(zhǔn)時(shí)差的校正;
e)在給定時(shí)間點(diǎn)上，為不同的源和接收機(jī)間的水平錯(cuò)距值測(cè)量所述集合波的幅度;以及f)直接從所述測(cè)得的幅度水平錯(cuò)距值確定剪切波反射信息，使得所述信息等效于通過(guò)用剪切波直接測(cè)得的剪切波反射信息。
現(xiàn)將參考各附圖對(duì)本發(fā)明進(jìn)行更詳細(xì)地描述，附圖中

圖1說(shuō)明利用一系列源和接收機(jī)的一個(gè)海底地震勘探系統(tǒng);
圖2說(shuō)明一個(gè)壓縮反射波軌跡的CDP集合和表示具有水平錯(cuò)距的軌跡幅度變化曲線圖。
圖3以流程圖的吊式說(shuō)明一種根據(jù)在某一進(jìn)行地震勘探的區(qū)域中所得到的壓縮波反射特性而為該區(qū)域估算剪切波反射特性的方法;
圖4a至d以流程圖的形式說(shuō)明另外的一些用于把由地震勘探中所得到的壓縮波特性同用圖3的方法所估算出的偽剪切波反射特性進(jìn)行對(duì)比以便為引起壓縮波反射和人或偽剪切波反射特性的異?，F(xiàn)象而確定存在有烴類物質(zhì)可靠標(biāo)志的方法;
圖5說(shuō)明通過(guò)圖1和2的方法所得到的一種壓縮波地震剖面。
圖6說(shuō)明一種通過(guò)圖3的方法，同時(shí)利用圖5的數(shù)據(jù)而估算出的偽剪切波地震剖面;
圖7說(shuō)明圖5的壓縮波平方值的地震剖面;以及圖8說(shuō)明通過(guò)把壓縮波反射信息與估算出的偽剪切波反射信息進(jìn)行比較的比較分析方法而計(jì)算出的一種地震剖面。
圖9以流程圖的方式說(shuō)明一種對(duì)地震數(shù)據(jù)的殘留標(biāo)準(zhǔn)化時(shí)差(RNMO)進(jìn)行校正的方法;
圖10a說(shuō)明被構(gòu)成圖9方法的一部分的推演濾波器的帶通特性;
圖10b說(shuō)明圖10a的推演濾波器的頻域構(gòu)成;
圖11說(shuō)明對(duì)地震數(shù)據(jù)進(jìn)行RNMO校正的效果;以及圖12以流程圖的方式說(shuō)明對(duì)所產(chǎn)生的壓縮波進(jìn)行計(jì)算其入射角的另一種方法。
參考各附圖，圖1說(shuō)明一般用于海底地震勘測(cè)的源和接收機(jī)組。應(yīng)當(dāng)清楚，所述機(jī)組同樣也適用于陸地上的地震勘測(cè)，并且本發(fā)明可應(yīng)用于這兩種勘測(cè)所收集的數(shù)據(jù)。但是當(dāng)用于陸地勘測(cè)時(shí)，就要用一附加的靜態(tài)校正(這將在下面詳述)，而這種校正在海上環(huán)境中是不需要的。
在收集數(shù)據(jù)期間，一般是使源S1……Sn順序啟動(dòng)以產(chǎn)生部份地從各地表界面層反射回來(lái)的遞減聲波，同時(shí)產(chǎn)生由接收機(jī)R1……Rn所接收并記錄的遞升反射信號(hào)(軌跡)。在后繼的為產(chǎn)生地震剖面的處理過(guò)程期間，所記錄的信號(hào)經(jīng)常是收集作為共深點(diǎn)(CDP)堆的。在圖2中為離得近和遠(yuǎn)的各配對(duì)的源和接收機(jī)以及介于二者距離間的各配對(duì)的源和接收機(jī)示出一個(gè)典型的共深點(diǎn)集合。在這個(gè)集合中的各信號(hào)都已經(jīng)用某一公知方法對(duì)其進(jìn)行過(guò)標(biāo)準(zhǔn)時(shí)差(NMO)校正的。如果涉及陸地發(fā)射，則在圖2中所記錄的信號(hào)集合也都已預(yù)先經(jīng)過(guò)表面一致性靜態(tài)校正處理過(guò)的。NMO校正和表面一致性靜態(tài)校正二者都是本技術(shù)中所公知的，因而這里不再重復(fù)對(duì)此詳述了。然后，通常對(duì)圖2所述的信號(hào)集合進(jìn)行CDP堆疊以便提高反射信號(hào)的信噪比，于是把該堆疊成的軌跡用來(lái)產(chǎn)生壓縮反射波特性的地震剖面。
本發(fā)明的方法使用了對(duì)水平錯(cuò)距進(jìn)行過(guò)CDP集合的壓縮波的幅度變化值〔在這里表示為Ap(X)騁員閽詬魘奔溲《ǖ閔瞎浪慍黽羥脅ǚ瓷涫藎檬縈摯捎美床嚶τ謁鯟DP堆疊的壓縮波地震剖面的偽剪切波地震剖面。正如在本技術(shù)領(lǐng)域中所公知的，水平錯(cuò)距指的是用以產(chǎn)生某一信號(hào)軌跡的、配對(duì)的源和接收機(jī)間的距離。圖2還示出隨著在兩個(gè)時(shí)間選定點(diǎn)t1、t1+△t上所集合成的軌跡的水平錯(cuò)距而變化的幅度變化值A(chǔ)p(X)。
為了搞清本發(fā)明的方法，就有必要對(duì)于隨水平錯(cuò)距而變的幅度變化產(chǎn)生影響的各種因數(shù)進(jìn)行討論。這些因數(shù)是P波反射系數(shù)Rp(在某一反射界面上的壓縮波的反射系數(shù))，各傳送通道因數(shù)(聲波擴(kuò)散和傳播特性)，源-接收機(jī)組和地表界面層的位置的幾何特征、以及諸如傳感器耦合等的隨機(jī)變化(主要是用于陸地式發(fā)射的小型地震儀)以及淺水衰減。關(guān)于這些因數(shù)，通過(guò)適當(dāng)選擇發(fā)射條件可使大多數(shù)因數(shù)(Rp除外)對(duì)值A(chǔ)p(X)影響最小。但是所有這些因素(Rp除外)都不能徹底被消除，而在本發(fā)明的處理步驟中將得到補(bǔ)償，這些步驟將在下面進(jìn)行更詳細(xì)的描述。在本發(fā)明中把影響Ap(X)值的各種因數(shù)進(jìn)行量化，以便為提供剪切波反射波幅值的信息而可確定因剪切波反射波所造成的Ap(X)值的變化影響。
所述反射系數(shù)Rp是必須用本發(fā)明的方法來(lái)確定的最重要的因數(shù)之一，眾所周知的Zoeppritz公式涉及所述反射系數(shù)Rp同各變量ρ、Vp、Vs、ρ′、V′p、V′s的關(guān)系，其中ρ、Vp和Vs分別是在某一地層界面處、上部介質(zhì)中的密度、壓縮波速度和剪切波速度，而ρ′、V′p和V′s分別是在該界面處、下部介質(zhì)中的密度、壓縮波速度和剪切波速度。
在考慮Zoeppritz公式中的這些因數(shù)的關(guān)系的同時(shí)，可作某些假設(shè)。首先假設(shè)，穿過(guò)某一反射邊界的相對(duì)速度/密度的變化都是小的，其次假設(shè)配對(duì)的源和接收機(jī)間的入射角θ都是小的。一般，為了保證第二假設(shè)是準(zhǔn)確的，在所收集的各軌跡中可使用同源和接收機(jī)間的法線所成的角為22°以下的各種角度。在這些假設(shè)的簡(jiǎn)化條件下，Zoeppritz冪級(jí)數(shù)公式可擴(kuò)展成它的分量，通過(guò)略去其二次項(xiàng)而產(chǎn)生下列簡(jiǎn)式Rp(θ)≈P+Qsin2θ (1)P＝RPO(2)Q＝RPO-8T2RSO+(2T2-1/2)△ρ/ρ (3)其中△ρ＝ρ-ρ′;
θ是同配對(duì)的源和接收機(jī)間的聲波通道的法線所成的角;以及T＝Vs/Vp的平均值 (4)R 是規(guī)定的入射P波反射系數(shù);而Rso是規(guī)定的入射S波反射系數(shù)(Rpo和Rso都代表垂直于某一反射邊界而傳送的P和S波的系數(shù))。
對(duì)于上面討論過(guò)的各一次因數(shù)(例如傳送通道因數(shù)、幾何特征、隨機(jī)變化等)來(lái)說(shuō)，它們都沒(méi)有通過(guò)恰當(dāng)?shù)剡x擇發(fā)射條件和預(yù)處理而消除，相反、它們?nèi)栽谟绊懼S水平錯(cuò)距而變的幅度變化Ap(X)，該Ap(X)可由下式表達(dá)
Fi～1+gisin2θ (5)式中i＝1至N個(gè)影響因數(shù)。
因此，可以把所有剩下的影響因數(shù)總的表達(dá)為Ap(θ)≈(F1F2…Fn)Rp(θ) (6)也可將Ap(θ)值表達(dá)為下面的通式Ap(θ)＝P+Qsin2θ (7)其中P＝Rpo(8)以及Q＝GRpo-8T2RSO+(2T2-1/2) (9)其中G =ΣiNg i]]>(10)公式(10)表示由上面論及的各種影響因數(shù)(例如傳送通道因數(shù)、幾何特征、隨機(jī)變化等，它們都未通過(guò)發(fā)送條件或預(yù)處理而從各集合的軌跡中消除)所引起的對(duì)集合軌跡中的總的幅度影響。
對(duì)于以已知水平錯(cuò)距所記錄的軌跡來(lái)說(shuō)，公式(7)可改寫為Ap(x)～P+Qfx2(11)其中X表示由圖1的源和接收機(jī)組的空間定位所確定的已知水平錯(cuò)距值，Ap(X)是為某一確定的水平錯(cuò)距X而測(cè)得的幅值，而f表示為把水平錯(cuò)距換算成入射角(它是根據(jù)NMO校正而導(dǎo)出的)的換算因數(shù)。
f＝(VINT/tVst) (12)其中VINT是某一估算出的區(qū)間速度(P波通過(guò)研究中介質(zhì)的速度)，而Vst是某一估算出的堆垛速度(把某一堆棧排成一行行軌跡所用的有效速度)，以及t是一種雙向傳送時(shí)間。
公式(11)可對(duì)每個(gè)選定時(shí)間點(diǎn)求解P和Q值，同時(shí)利用一種例如迭代最小均方程序以求出P和Q點(diǎn)的最佳“擬合值”。根據(jù)所算得的P和Q值，用公式(8)可求出各Rpo值。此外，如果T，G，和△ρ/ρ各值都是已知，則可將所確定的Q值代入公式(9)以確定Rso(即規(guī)定的剪切波反射系數(shù))。T值〔公式(4)〕對(duì)Q的關(guān)系可說(shuō)明如下T= (Q)/(VP/VS) Q1.4 GRPO-4RSO+1/2△ρ/ρ2.0 GRPO-2RSO4.0 GRPO-1/2RSO-3/8△ρ/ρ表中指出，如T＝Vp/Vs＝2.0，則△ρ/ρ在計(jì)算Q時(shí)為零即消除了一個(gè)因數(shù)。在許多場(chǎng)合下，T值甚至為2.0，但是現(xiàn)在就可利用壓縮和剪切波鉆孔記錄來(lái)獲得各Vp/Vs的直接測(cè)量值，因而在確定了△ρ/ρ和G值的同時(shí)就可根據(jù)T、Q和Rpo的確定值來(lái)確定Rso值。
△ρ/ρ值經(jīng)常可根據(jù)已知的地表特性而相當(dāng)準(zhǔn)確地預(yù)測(cè)出。一般，△ρ/ρ～1/nRpo，其中n是對(duì)某一特定區(qū)域的觀察值，并且一般等于5。也可直接從某一勘探區(qū)域的實(shí)際測(cè)井記錄中作為測(cè)得值而獲得△ρ/ρ值(因?yàn)椤鳓?ρ＝△Vp/Vp)。
G值是可估算的;因?yàn)楫?dāng)源和接收機(jī)組中使用小型陸地振動(dòng)速度地震儀時(shí)，G通常具有G～0的值。對(duì)于海洋爆破所用的水中地震檢測(cè)源和接收機(jī)來(lái)說(shuō)G通常具有G～1的值。但是也有可能更準(zhǔn)確地估算G值。茲將其計(jì)算技術(shù)描述如下選出公式(2)和(3)P＝Rpo(2)Q＝GRpo-8T2Rso+(2T2-1/2)△ρ/ρ (3)對(duì)于所述CDP軌跡集合中的每個(gè)時(shí)間點(diǎn)而言，P和Q是可用例如最小均方計(jì)算方法來(lái)估算的。通過(guò)假設(shè)△ρ/ρ～1/n(Rpo)，就可使公式(3)中的密度項(xiàng)△ρ/ρ去除。于是公式(3)就化成Q=G+( (2T2-1/2)/(nRPO) )-ST2RSO(13)在簡(jiǎn)化Zoeppritz方程的過(guò)程中，曾經(jīng)假定T和G是隨時(shí)間(因而當(dāng)然是深度)而緩慢地變化的。因此估算它們的值所需要的一切就是一種由它們的各組成部份組成的、平滑的或經(jīng)低通濾波的型式。這可根據(jù)在就近的各鉆孔中所測(cè)得P波和S波鉆探記錄數(shù)據(jù)中取得。因此，當(dāng)P和S波的速度(V′p，V′s)都取自相應(yīng)的鉆探記錄中時(shí)，函數(shù)T就可估算為T≈< (V′S)/(V′P) >(14)其中<>指的是一種低通濾波操作。
其次，如果合成的P和S波的反射波R′po和R′so都是從P和S波鉆探記錄中導(dǎo)出的，則通過(guò)組合公式(8)和(9)以及低通濾波而形成下面的估算函數(shù)GG≈( (<PQ>)/(P2) )+8T2((R′POR′SO))/(R′2PO) - ((2T2-1/2))/(n) (15)一旦G是已知值，則(9)式就能解出Rso。
Rso＝(1/n(G+2T2-1/2)P-Q)/8T2(16)因此，首先利用迭代過(guò)程解方程(7)以確定P(其中P＝Rpo)和Q，然后解方程(9)以確定Rso，所述壓縮波反射系數(shù)Rpo和剪切波反射系數(shù)Rso二者都可根據(jù)測(cè)得的隨水平錯(cuò)距而變化的幅度變化值A(chǔ)p(X)而確定。可以用這些值來(lái)同時(shí)作出壓縮波和剪切波地震剖面以便判斷和分析。
接著轉(zhuǎn)向圖3，現(xiàn)將公開(kāi)本發(fā)明的方法。這一方法是可借助一種能在現(xiàn)有的通用數(shù)字計(jì)算機(jī)上運(yùn)行的典型計(jì)算機(jī)程序而實(shí)施的。
在步驟101開(kāi)始，運(yùn)用圖1的地震數(shù)據(jù)采集設(shè)備得到VRMS、T、S諸值和采集(發(fā)射的)和物理的參數(shù)，并輸入到計(jì)算機(jī)，其中VRMS表示堆垛速度(Vst)，T＝<Vs′/Vp′> (17)該T值是一種從實(shí)際鉆孔記錄中所得到的濾波后平滑的Vs′/Vp′型式;以及S= (<RP′RS′>)/(<RP′2>)該S值也是從實(shí)際鉆孔記錄中所得到的濾波后(平滑的)統(tǒng)計(jì)數(shù)據(jù)。這些采集參數(shù)包括用于數(shù)組參數(shù)中的源和接收機(jī)的類型方面的參數(shù)。所述物理參數(shù)包括彈性衰減參數(shù)q。
接著進(jìn)行到步驟103，將按NMO校正的CDP軌跡集合讀入計(jì)算機(jī)(這些軌跡也都是在需要時(shí)預(yù)先處理過(guò)的，以便作附近表面的靜態(tài)校正，例如對(duì)陸地發(fā)射的校正)。然后在步驟104，按照?qǐng)D9的方法對(duì)NMO校正后的CDP軌跡進(jìn)行殘留的NMO校正。
偽剪切波的估算值對(duì)NMO速度的少量誤差都是很敏感的。標(biāo)準(zhǔn)時(shí)差是一種對(duì)于隨增加著的源和接收機(jī)間的水平錯(cuò)距而延長(zhǎng)反射時(shí)間的有規(guī)則的移位。在讀入計(jì)算機(jī)之前，上面所用校正CDP軌跡集合的一種典型標(biāo)準(zhǔn)時(shí)差的校正方法將使沿著時(shí)間層位的各地震事件(seismicevents)拉平。雖然NMO校正方法多少有點(diǎn)效果，但它卻趨于使所述數(shù)據(jù)遺留某些殘留標(biāo)準(zhǔn)時(shí)差?？梢园迅鳂?biāo)準(zhǔn)時(shí)差的誤差表達(dá)為可顯著降低地震位移測(cè)量結(jié)果質(zhì)量的堆垛速度方面的誤差，并將引起各幅值偏差以致又將掩蓋來(lái)自土壤內(nèi)反射層的真實(shí)幅度-水平錯(cuò)距特性。對(duì)殘留標(biāo)準(zhǔn)時(shí)差的校正是根據(jù)主要反射層的幅度偏移性能的確定對(duì)每個(gè)時(shí)間取樣進(jìn)行估算的。
現(xiàn)在來(lái)描述殘留標(biāo)準(zhǔn)化時(shí)差的校正方法，該方法通過(guò)對(duì)在已經(jīng)完成標(biāo)準(zhǔn)時(shí)差(NMO)校正之后為獲得地震數(shù)據(jù)的原始反射層R(X，T)幅度偏移特性的每個(gè)時(shí)間采樣，估算其殘留標(biāo)準(zhǔn)時(shí)差。準(zhǔn)確到一級(jí)(firstorder)的物理特性變化的所需結(jié)果是R(xTO)=ρ(TO)+[QO- (δV)/(Vold) TO( (VO)/(Vold) )2ρ′] (X2)/(TO2VO2) (18)
其中Vold+δV是準(zhǔn)確的堆垛速度;
TO是在規(guī)定入射角時(shí)的雙向傳送時(shí)間;
X是水平錯(cuò)距;
VO是基準(zhǔn)速度;
P和QO均是Zoeppritz公式的一次展開(kāi)式(first order expansion)對(duì)于P波反射比的系數(shù);而P′是P的一次導(dǎo)數(shù)根據(jù)公式(18)，δV是根據(jù)如下的時(shí)間偏移統(tǒng)計(jì)而估算的δV=P ′QO-δVVo l dTO(VOVo l d)2P ′2P ′2]]>(19)其中<>是在某一以TO為中心的窗口中的時(shí)間平均值與<P′P′>相比，<P′QO>將多般是可忽略不計(jì)的。因此具有時(shí)間量綱的公式(19)就變?yōu)棣膖=- (δV)/(Vold) TO( (VO)/(Vold) )2(20)為使時(shí)間偏移t與殘留NMO發(fā)生關(guān)系，就產(chǎn)生TRES=TO-δt( (X2)/(TO2VO2) )(21)借助圖9，公開(kāi)了本發(fā)明的殘留NMO影響的校正方法。在步驟1030，把可靠的CDP數(shù)據(jù)的各隅角頻率(f1、f2、f3、f4)(如在圖10a中所說(shuō)明、并由眾所周知的CDP數(shù)據(jù)分析方法所確定)，諸如NMO-校正后的CDP集合、軌跡長(zhǎng)度、采樣速率以及水平錯(cuò)距等物理參數(shù)輸入某一通用數(shù)字計(jì)算機(jī)。在本發(fā)明的另一些實(shí)施例中也可將關(guān)于信/噪比(S/Nmin)的下限和堆垛速度曲線圖包括在作為該數(shù)字計(jì)算機(jī)的輸入中。當(dāng)需要選定加有時(shí)間偏移的數(shù)據(jù)時(shí)，就可使用該信/噪比。當(dāng)需要對(duì)堆垛速度進(jìn)行校正時(shí)，就可使用這些堆垛速度曲線圖。該數(shù)字計(jì)算機(jī)或者可以是與執(zhí)行圖3的估算偽剪切波算法的同一數(shù)字計(jì)算機(jī)，或者可以是與通過(guò)傳統(tǒng)方法確定該剪切波估算過(guò)程中所用數(shù)字計(jì)算機(jī)相連接的第二數(shù)字計(jì)算機(jī)。
在步驟1050，計(jì)算帶限推演濾波器，其方法是采用一個(gè)加權(quán)的、“最小二乘方”法，該二乘方適于由隅角頻率f1、f2、f3和f4、采樣速率和軌跡長(zhǎng)度(即濾波器長(zhǎng)度)〔這些用于形成推演濾波器〕轉(zhuǎn)化而來(lái)的某一梯形加權(quán)函數(shù)導(dǎo)數(shù)的“iω”頻率響應(yīng)。參考圖10b，應(yīng)當(dāng)注意從時(shí)間零起以+n△t隔開(kāi)的一對(duì)不對(duì)稱的峰值+An的離散傅里葉變換是-2An＊Sin(2πn△f)。因此，將求解濾波器最佳(2m+1)點(diǎn)的過(guò)程簡(jiǎn)化為使m個(gè)正弦波擬合到某一線性曲線的過(guò)程。
返回到圖9的步驟1070，進(jìn)行對(duì)于該輸入的CDP數(shù)據(jù)執(zhí)行RNMO校正的迭代次數(shù)的選定。已經(jīng)證明凡是堆垛速度誤差呈現(xiàn)為大的殘留NMO誤差，即5-10%誤差者，將在少于10次迭代后就可使諸反射事件獲得近于完善的時(shí)間校正。對(duì)于本發(fā)明的一般使用來(lái)說(shuō)，所述CDP的精心分析將會(huì)使速度誤差限制到小于3%。因此認(rèn)為選定五次迭代可特別滿意地消除所述預(yù)期的速度誤差。在步驟1090來(lái)擬合每個(gè)時(shí)間點(diǎn)上的水平錯(cuò)距特性，從而為每個(gè)時(shí)間點(diǎn)計(jì)算出噪聲的估算值。對(duì)于這樣一種計(jì)算來(lái)說(shuō)，使所述水平錯(cuò)距特性在所述輸入CDP特性的每個(gè)時(shí)間點(diǎn)上適于下面的公式P+Q( (X2)/(VO2TO2) )(22)雖然為了計(jì)算穩(wěn)定而推薦某一任意選定的VO＝6，000英尺/秒，但VO可假設(shè)為任意可變的速度。假設(shè)R(TO，X)為在NMO校正時(shí)間TO時(shí)根據(jù)所述輸入CDP軌跡的水平錯(cuò)距數(shù)據(jù)，則R(TO,X)～P(TO)+Q(TO) (X2)/(VO2TO2) 3)于是就可利用在本技術(shù)領(lǐng)域中公知的方法，從公式(23)解出P(TO)和Q(TO)。在可測(cè)得S/Nmin值的本發(fā)明的實(shí)施例中，在這里也可按照下式計(jì)算出信噪比(ρ2(TO))/(ρ2(TO)+σ2(TO)) (24)其中方差δ2(TO)是按照下式估算的δ2(TO)=( 1/(N-2) )[(R(TO′X)-P(TO)-Q(T0)( (X2)/(VO2TO2) )2] (25)
繼續(xù)進(jìn)行到步驟1110，使在步驟1050確定的推演濾波器與在步驟1090的向前10個(gè)P軌跡的傅里葉變換進(jìn)行卷積以計(jì)算頻域P′軌跡，然后再將其變換回到所述時(shí)域。在步驟1130，根據(jù)下面的關(guān)系式確定V/V<P′Q>/<P′2> (26)其中<P′Q>/<P′2>是在各移動(dòng)的時(shí)間窗口中算出的。
由于已知Q=-( (δV)/(Vold) )TO( (VO)/(Vold) )2ρ′(27)P′是通過(guò)對(duì)某一以TO為中心的時(shí)間窗口進(jìn)行加權(quán)最小二乘方擬合而解出的，于是(δV)/(Vold) TO( (VO)/(Vold) )2= (Σwtpt′Qt)/(Σwtpt′2) ≤δt(28)式中Wt是對(duì)該窗口內(nèi)的Tth采樣的余弦加權(quán);
P′、Q和P′2是通過(guò)與某一自乘的(raised)余弦加權(quán)函數(shù)W進(jìn)行卷積然后通過(guò)除法而形成的以便求取每個(gè)時(shí)間點(diǎn)上的估算值。
在步驟1150，為利用各種信噪比的本發(fā)明的實(shí)施例而確定各可靠的時(shí)間。當(dāng)某一信噪比已經(jīng)提供時(shí)，該輸入就對(duì)于應(yīng)當(dāng)使用多少個(gè)由上面所確定的(δV，時(shí)間)值提供一種截止值。在所述曲線的確定過(guò)程中將在各數(shù)值之間利用這些可用值的線性判讀。
繼續(xù)進(jìn)行到步驟1170，以正好與堆垛初始NMO校正數(shù)據(jù)的相同方法來(lái)形成一種堆垛速度函數(shù)Vst。在步驟1190，按照下面的公式修改這些堆垛速度Vnew=Vold+[ (δV)/(Vold) TO(VO)/(Vold) ] (Vold)/(TO) ( (Vold)/(VO) )2(29)
繼續(xù)進(jìn)行到步驟1210，確定通過(guò)速度校正的CDP數(shù)據(jù)的所需迭代次數(shù)是否已經(jīng)完成。如果需要附加該CDP迭代次數(shù)，則程序就返回到步驟1090以便對(duì)在步驟1190中所計(jì)算的校正后速度CDP數(shù)據(jù)進(jìn)行下一次速度校正的迭代。如果在步驟1210確定所有CDP集合已經(jīng)完成必需次數(shù)的速度校正的迭代，從而已經(jīng)消除了NMO校正后CDP數(shù)據(jù)的RNMO水平錯(cuò)距;則本發(fā)明方法就以步驟1250返回到圖3的步驟105。
通過(guò)參考圖11就可看到一個(gè)從地震數(shù)據(jù)中消除RNMO水平錯(cuò)距的實(shí)例，圖11表明在RNMO校正之前的NMO校正的CDP數(shù)據(jù)以及準(zhǔn)確到5、10、15和20次迭代的RNMO數(shù)據(jù)。在圖3的步驟105，在某一時(shí)間記數(shù)器中置定時(shí)間T＝0。隨后在步驟107，則對(duì)于該時(shí)間記數(shù)器中所置定的時(shí)間可按下式計(jì)算出“深度”＝ZZ＝VRMS·時(shí)間/2 (30)這是一種速度VINT的估算區(qū)間，以及關(guān)于對(duì)于水平錯(cuò)距特性的幅度產(chǎn)生影響的各種因數(shù)的決定性分量Gdet。例如，某些決定性因數(shù)是由某些已知公式所支配。這些因數(shù)包括球面展開(kāi)(spherical spreading)，它可定義為SS＝-(1/2)V(3)/V3INT(31)其中V3＝(1/T)∫TOVINT3·tdt (32)其中T＝單向傳送時(shí)間;源、接收機(jī)組的方向性系數(shù)(SRD)，它可定義為SRD＝-1/6(πsl)2/VINT2(33)其中S表示源信號(hào)的頻率l＝機(jī)組長(zhǎng)度;以及衰減(AT)，它可定義為
AT＝-/q(πsl/VINT2)VRMS(34)利用這些公式就可為這些影響因數(shù)而確定Gdet。其后，在步驟109，使計(jì)數(shù)器置位，此時(shí)用于該集合軌跡的近距水平錯(cuò)距值×具最小的入射角θ，而在步驟111，這一入射角θ是用下式估算的sin2θ=( (VINT)/(VRMS) · (X2)/(X2+4X2) )(35)在后繼步驟113，使水平錯(cuò)距記數(shù)器X遞增，它的下一個(gè)水平錯(cuò)距值為X(向具有最遠(yuǎn)水平錯(cuò)距的軌跡移動(dòng)的集合中的下一軌跡的水平錯(cuò)距)，然后在步驟115，做出是否所有水平錯(cuò)距(所有軌跡)已經(jīng)處理完畢的判定。如果為“否”、則就為在步驟113中所遞增的最新水平錯(cuò)距估算入射角θ。
雖然上述方法可產(chǎn)生關(guān)于入射角θ的滿意估算值，但是也可利用下面的另一方法更精確地計(jì)算入射角θ，現(xiàn)將參考圖12來(lái)描述這種方法。
地震勘探一般是對(duì)具有不同地質(zhì)特性的若干層組成的地下地層而進(jìn)行的。例如參見(jiàn)圖1，圖中示出所要研究的兩個(gè)作為舉例的地下地層。為進(jìn)行下面的推導(dǎo)起見(jiàn)，在這里把Ui定義為層1頂部的速度，把t定義為所產(chǎn)生的地震能量進(jìn)入層1時(shí)的時(shí)間，Vi是層1內(nèi)的速度，而α具有加速度的量綱并且是一種輸入?yún)?shù)。因此，地層1內(nèi)的速度Vi是用如下式的步驟2030來(lái)確定的Vi＝Ui+αi(t-ti) (36)然后按照下面的方法、在步驟2050確定各軌跡的Ui值。從下面的速度均方值的定義開(kāi)始V2i＝(1/ti)∫tioV2dt (37)于是ti+1Vi+12-tiVi2=∫trtr-1V1dt(38)]]>假設(shè)該地層的速度是線性地增加的，則將給出ti+1Vi+12-tiVi2=∫trtr-1[Ui+ai(t-ti)]2dt(39)]]>
該式結(jié)果產(chǎn)生下面Ui的二次方程表達(dá)式(ti+1-ti)U2i+αi(ti+1-ti)2Ui+(αi/3)(ti+1-ti)3-ti+1Vi+12+tiV2i＝0 (40)在用步驟2050算出Ui以后，就在步驟2070通過(guò)對(duì)每個(gè)界面處的兩個(gè)速度估算值進(jìn)行平均，并對(duì)各結(jié)果值按某一三次樣條函數(shù)進(jìn)行內(nèi)插而獲得一種濾清的區(qū)間速度函數(shù)。
然后在步驟2090根據(jù)用步驟2070算得的濾清的區(qū)間速度函數(shù)的水平錯(cuò)距和片刻計(jì)算各入射角。為下面的推導(dǎo)起見(jiàn)，使用下面的符號(hào)t＝單向傳送時(shí)間;
ρ＝密度;
α＝壓縮波速度;
β＝剪切波速度;
θ＝入射角(其中θi將表示第i個(gè)反射層的入射角，而θO將表示在表面處的入射角);
P＝sinθiαi＝sinθOαOV(n)＝(1/t)∫tOαn(t′)dt′Zi＝第i層厚度;
Xi＝為第i層的射線所橫過(guò)的側(cè)距;
X＝總的水平錯(cuò)距于是對(duì)任一地層來(lái)說(shuō)Xj＝Zjtanθj＝Zj[pαj+1/2p3α3j…](41)通過(guò)積分就產(chǎn)生x＝∫zOpαdz+(1/2)∫zOp3α3dz+… (42)代入dZ＝αdt，因而X＝ρV2t+(1/2)ρ3V4t+… (43)通過(guò)變換，并代入P＝(Sinθ/α)，則得
sinθ= (αx)/(V2t) - (α)/2 (V4)/(V2) ( (X)/(V2t) )3+…(44)sin2θ= (α2)/(V2) [- ((2x)2)/((2x)2+4V2t2) )(45)在本發(fā)明的這個(gè)實(shí)施例中，在算出入射角以后，該算法在步驟2110返回到圖3的步驟1107，在哪里為了按照?qǐng)D3的步驟111至115的方法估算各對(duì)應(yīng)的入射角θ而已經(jīng)處理了所有的水平錯(cuò)距以后或者利用圖12的方法算出所有的入射角以后，該處理計(jì)算機(jī)就繼續(xù)進(jìn)行到步驟117，在那里，利用公式(7)和最小均方的遞歸處理技術(shù)計(jì)算P和Q值。接著，計(jì)算機(jī)進(jìn)到步驟119，在那里對(duì)不是由步驟107中算得的剪切波速度分量Gdet所造成的水平錯(cuò)距特性估算其幅度分量。這個(gè)估算值G是按照下面公式而進(jìn)行估算的G≈ (PQ)/(P2) +8T2S- ((2T2-1/2))/(n) -Gdet(46)隨后，用步驟121，根據(jù)下式算出剪切波反射比分量RRS=-[Q-[Gdet+δG+ (2T2-1/2)/(n) ]P]/8T2(47)該式是公式(16)的一種重寫型式。于是，計(jì)算機(jī)進(jìn)到步驟123，在那里使時(shí)間計(jì)數(shù)器遞增。在下一步驟125中，確定被處理的軌跡是否已經(jīng)結(jié)束。如果沒(méi)有結(jié)束，則計(jì)算機(jī)返回到步驟107，利用一個(gè)新的軌跡時(shí)間值，并重復(fù)進(jìn)行步驟107到123。如果，在步驟125中確定已經(jīng)結(jié)束，則隨后計(jì)算機(jī)進(jìn)到步驟127，在哪里確定所有的NMO校正的集合是否都已處理過(guò)。如果已全部處理過(guò)，則程序就結(jié)束;如果不是，則計(jì)算機(jī)進(jìn)行到步驟103，并開(kāi)始處理一個(gè)新的軌跡集合。
在步驟121中，計(jì)算機(jī)算出各剪切波反射比系數(shù)，然后才可將其加到這些集合的軌跡中以產(chǎn)生代表剪切波反射波軌跡的各軌跡。于是這些剪切波軌跡就可以傳統(tǒng)的方式(CDP堆垛)來(lái)產(chǎn)生偽剪切波地震剖面。
在用步驟121完成某一偽剪切波地震剖面計(jì)算之時(shí)，開(kāi)始進(jìn)行壓縮波和剪切波地震剖面的幅度異?，F(xiàn)象的對(duì)比分析，以便判定是否存在有正的烴類物質(zhì)標(biāo)記。通過(guò)分別參考圖6和圖7就可看到某一地層的估算的偽剪切波剖面和壓縮波剖面。在圖6和7中，重點(diǎn)突出了包括可指示出或者存在烴類物質(zhì)或者存在巖性地層的幅度異常現(xiàn)象的各地震剖面部份。參考圖4a至d的流程圖可更精確地描述另一些比較壓縮波地震剖面和相關(guān)的估算偽剪切波地震剖面的方法。
借助于圖4a，本發(fā)明的第一種比較分析方法以步驟151開(kāi)始，在相同的增益條件下輸入各地震剖面Rp和Rs。在步驟153，選定包括指出正在研究的地下地層中可能存在烴類物質(zhì)的幅度異?，F(xiàn)象的壓縮波地震剖面部份。用步驟155選擇對(duì)應(yīng)于該已選Rp剖面部分的Rs剖面部分。在步驟157比較所選定的Rp和Rs剖面。通過(guò)為所選定的異?，F(xiàn)象而確定Rp-Rs的值就可進(jìn)行Rp對(duì)Rs的比較。如果在步驟159，指示該兩種所選異?，F(xiàn)象之間有著高度的相似性(即Rp-Rs＝0)，則在步驟161中就對(duì)所選定的Rp剖面部分判定為存在負(fù)的烴類物質(zhì)。如果在步驟159中指示該兩種所選異常現(xiàn)象之間的相似程度低(即Rp-Rs并不趨于零)，則在步驟163中，對(duì)所選定的Rp剖面部分判定為存在正的烴類物質(zhì)。如果在步驟165確定所研究的最后壓縮波地震剖面幅度異?，F(xiàn)象已經(jīng)處理過(guò)了，則該算法就結(jié)束。如果還有別的壓縮波地震剖面幅度異常現(xiàn)象待研究，則該算法就在步驟167返回到步驟153以便作進(jìn)一步處理。
其次借助圖4b，用步驟176開(kāi)始本發(fā)明的第二個(gè)比較分析方法，在該步驟中以相等的增益輸入各地震剖面Rp和Rs。緊接著到步驟178，計(jì)算壓縮波地震剖面同偽剪切波地震剖面的幅度差值(Rp-Rs)。在步驟180，選定包括可指示出在該所研究的地下地層中可能存在烴類物質(zhì)的幅度異?，F(xiàn)象的壓縮波地震剖面部分。在步驟182選定對(duì)應(yīng)于所選Rp剖面部份的(Rp-Rs)剖面部分。在步驟184，將對(duì)應(yīng)于先前選定的Rp剖面部分的所選Rp-Rs剖面同所選定的Rp剖面進(jìn)行對(duì)比。如果在步驟186，指示出低的相似度(即Rp剖面的所選定部分的幅度異?，F(xiàn)象大，而(Rp-Rs)剖面的所選定部分即使有幅度異常也極小)，則在步驟188，對(duì)該Rp剖面的所選定部分判定為存在負(fù)的烴類物質(zhì)。如果在步驟186指示出高度的相似性(即Rp和(Rp-Rs)剖面的所選定部分存在著相同的幅度異常現(xiàn)象)，則在步驟190，對(duì)該Rp剖面的所選定部分判定為存在正的烴類物質(zhì)。如果在步驟192，判定所研究的最后壓縮波地震剖面幅度異?，F(xiàn)象已處理過(guò)，則該算法就結(jié)束。如果還有別的壓縮波地震剖面幅度異?，F(xiàn)象待研究，則該算法就在步驟194返回到步驟180以便作進(jìn)一步處理。
其次借助圖4C，在步驟201，以相等的增益輸入地震剖面Rp和Rs，從而開(kāi)始本發(fā)明的第三種比較分析方法。在步驟203，對(duì)壓縮波地震剖面R乘方，以便計(jì)算出壓縮波地震剖面的R2p。緊接著，進(jìn)到步驟205，計(jì)算壓縮波和偽剪切波地震剖面的幅度之差(Rp-Rs)。在步驟207，通過(guò)將該幅度差值(Rp-Rs)乘以Rp而將該幅度差值定標(biāo)到與R2p地震剖面相同的數(shù)值。通過(guò)分別參考圖7和8，就可看到對(duì)同一地震數(shù)據(jù)進(jìn)行R2p和Rp(Rp-Rs)地震剖面的比較。在圖7和8中，重點(diǎn)突出了與圖5和6幅度異?，F(xiàn)象有關(guān)的地震剖面部分。
返回到圖4C，在步驟209選定一個(gè)壓縮波地震剖面部分，該剖面部份包括有一種可指示出在調(diào)查研究中的地下地層中可能有烴類物質(zhì)的一種幅度異?，F(xiàn)象。在步驟211，選定對(duì)應(yīng)于所選定的Rp剖面部分的Rp(Rp-Rs)和R2p的剖面部分。步驟213將對(duì)應(yīng)于先前所選定的Rp剖面部分的所選定的Rp(Rp-Rs)剖面和所選R2p剖面進(jìn)行對(duì)比。如果在步驟215，指示出低的相似度(即所選R2p剖面部分具有大的幅度異?，F(xiàn)象而所選Rp(Rp-Rs)剖面部分即使有幅度異?，F(xiàn)象也極小)，則在步驟217，對(duì)所選Rp剖面部分判定為存在負(fù)的烴類物質(zhì)。如果在步驟215指示出高的相似度(即所選R2p和Rp(Rp-Rs)剖面部分出現(xiàn)相似的幅度異?，F(xiàn)象)，則在步驟219，對(duì)所選Rp剖面部分判定為存在正的烴類物質(zhì)。如果在步驟221中判定該所研究的最后壓縮波地震剖面幅度異?，F(xiàn)象已經(jīng)處理過(guò)了，則結(jié)束該算法。如果還有別的壓縮波地震剖面幅度異?，F(xiàn)象待研究，則該算法就在步驟223返回到步驟209以便進(jìn)一步處理。
借助圖4d，在步驟226，以同一增益輸入地震剖面Rp和Rs，從而開(kāi)始本發(fā)明的第四種比較分析方法。在步驟228，使壓縮波地震剖面Rp乘以-1，以對(duì)該壓縮波地震剖面算出帶符號(hào)的SgnRp。緊接著，進(jìn)到步驟230，計(jì)算壓縮波和偽剪切波地震剖面的幅度差值(Rp-Rs)。在步驟232，選定一個(gè)壓縮波地震剖面部分，該剖面部分包括有一種可指示出在所研究的地下地層中存在烴類物質(zhì)的幅度異?，F(xiàn)象。在步驟234，選定SgnRp和對(duì)應(yīng)于所選Rp剖面部分的(Rp-Rs)剖面部分。在步驟236，將所選SgnRp剖面同對(duì)應(yīng)于先前所選定Rp剖面部分的所選(Rp-Rs)剖面進(jìn)行對(duì)比。如果在步驟238指示出低的相似度(即所選定的SgnRp剖面部分具有大的幅度異?，F(xiàn)象而所選定的(Rp-Rs)剖面部分幾乎沒(méi)有幅度異常現(xiàn)象)，則在步驟240，對(duì)所選的Rp剖面部分判定為存在負(fù)的烴類物質(zhì)。如果在步驟238，指示出高的相似度(即對(duì)所選SgnRp和(Rp-Rs)剖面部分呈現(xiàn)大的幅度異?，F(xiàn)象)，則在步驟242，對(duì)所選Rp剖面部分判定為存在正的烴類物質(zhì)。如果在步驟244，斷定所研究的最后的壓縮波地震剖面幅度異?，F(xiàn)象還未經(jīng)處理，則該算法就在步驟246返回到步驟232以便進(jìn)一步處理。
權(quán)利要求
1.一種從壓縮波反射信息產(chǎn)生剪切波反射信息的方法，特征在于該方法包括下列步驟a)從若干聲波源產(chǎn)生逐漸下降的聲波，并在若干接收機(jī)處接收和記錄在某一地下界面上產(chǎn)生的逐漸上升的波形；b)把所述錄下的反射波集合成具有某一公共反射邊界的不同源和接收機(jī)的水平錯(cuò)距組；c)對(duì)所述集合的反射波進(jìn)行標(biāo)準(zhǔn)時(shí)差校正；d)對(duì)所述標(biāo)準(zhǔn)時(shí)差校正后的波進(jìn)行殘留標(biāo)準(zhǔn)時(shí)差校正；e)在給定時(shí)間點(diǎn)上，對(duì)不同的源和接收機(jī)間的水平錯(cuò)距值測(cè)量所述集合波的幅度；以及f)直接從所述測(cè)得的幅度水平錯(cuò)距值確定剪切波反射信息以使所述信息等效于通過(guò)用剪切波直接測(cè)得的剪切波反射信息。
2.權(quán)利要求1的方法，特征在于其中步驟(d)包括校正按照下式校正后的所述標(biāo)準(zhǔn)時(shí)差波形的堆垛速度Vnew=Vold+[ (δV)/(Vold)To(Vo)/(Vold) ] (Vold)/(To) ( (Vold)/(Vo) )2其中Vnew是校正后的堆垛速度，Vold是測(cè)得的堆垛速度δV是速度校正量TO是在標(biāo)準(zhǔn)入射角下的雙向傳送時(shí)間，以及VO是某一基準(zhǔn)速度。
3.權(quán)利要求1或2的方法，特征在于還包括下列步驟-為所述集合后的記錄軌跡確定平滑后的區(qū)間速度的函數(shù);以及-根據(jù)所述水平錯(cuò)距和所述平滑后的區(qū)間速度函數(shù)確定所述集合后軌跡的入射角;從而直接根據(jù)所述測(cè)得的幅度錯(cuò)距值和所述集合后軌跡的所述入射角來(lái)確定所述剪切波反射信息。
4.權(quán)利要求3的方法，特征在于其中根據(jù)所述各水平錯(cuò)距和所述平滑后的區(qū)間速度函數(shù)確定所述集合后軌跡的入射角θ的步驟是按照下式計(jì)算的sin2θ= (α2)/(V2) [ ((2x)2)/((2x)2+4V2t2) ]其中α＝在給定深度下的區(qū)間速度，Y＝堆垛速度，X＝源/接收機(jī)間的水平錯(cuò)距，以及t＝傳向給定深度處的單向傳送時(shí)間。
5.任一前述權(quán)利要求的方法，特征在于其中步驟(f)包括如下步驟-為每個(gè)所述時(shí)間點(diǎn)確定某一值Gdet，該值表示隨水平錯(cuò)距而影響所述測(cè)得的幅度變化的已知因數(shù)的量化值;-為每個(gè)所述時(shí)間點(diǎn)估算隨水平錯(cuò)距不同而異的幅度分量G，該分量G不是歸因于某一剪切波分量，并且不是構(gòu)成Gdet的;以及-為每個(gè)所述時(shí)間點(diǎn)、按照下式確定剪切波反射比分量RsRs＝-[Q-[Gdet+δG+ (2T2-1/2)/(n) 〕p]/8T2其中P＝Rpo以及Q＝GRpo-8T2Rso+(2T2-1/2)△ρ/ρ式中Rpo＝標(biāo)準(zhǔn)的入射壓縮波反射系數(shù);Rso＝標(biāo)準(zhǔn)入射剪切波反射系數(shù);T＝(Vs/Vp′)平均值ρ＝密度G＝由隨水平錯(cuò)距而影焐ǚ缺浠乃鲆閻蚴鸕姆扔跋斕淖芏睢
6.權(quán)利要求 述的方法，特征在于其中值δG是按照下式δG≈ ((PQ))/((P2)) +8T2S- ((2T2-1/2))/(n) -Gdet其中<PQ>表示PQ的濾波后的值，而<P2>表示P2的濾波后值。
全文摘要
從壓縮波反射信息產(chǎn)生剪切波反射信息的方法，步驟為a)從若干聲源產(chǎn)生下降著的聲波，并在若干接收機(jī)處接收和記錄某一地下界面上產(chǎn)生的上升著的波；b)把所述錄下的反射波集合成具有公共反射邊界的不同源接收機(jī)間的錯(cuò)距組；c)對(duì)所述集合的反射波進(jìn)行標(biāo)準(zhǔn)時(shí)差校正；d)為所述殘留標(biāo)準(zhǔn)時(shí)差測(cè)量所述集合后的波幅；e)直接從所述測(cè)量結(jié)果確定剪切波反射信息以使其等效于利用剪切波直接測(cè)量所獲得的剪切波反射信息。
文檔編號(hào)G01V1/28GK1032591SQ88106870
公開(kāi)日1989年4月26日 申請(qǐng)日期1988年9月20日 優(yōu)先權(quán)日1987年9月21日
發(fā)明者托馬斯·杰爾斯特魯普·菲奇, 斯坦利·杰拉爾·拉斯特, 尼爾·羅德尼·戈因斯, 雷·史蒂文·斯普拉特 申請(qǐng)人:無(wú)比石油公司