地震資料數(shù)字處理-4

1、本章主要由以下幾部分組成：本章主要由以下幾部分組成： 4.1 概述概述 4.2 動校正動校正 4.3 速度分析速度分析 4.4 靜校正靜校正 4 4. .5 5 水平疊加水平疊加 疊加疊加是地震處理三大核心技術(shù)之一，其目的是壓制隨機壓制隨機干擾、提高地震信噪比。干擾、提高地震信噪比。 與疊加技術(shù)相關(guān)的研究內(nèi)容與疊加技術(shù)相關(guān)的研究內(nèi)容： 速度分析速度分析為疊加提供最佳疊加速度。 動校正動校正消除炮檢距對反射波旅行時的影響。 靜校正靜校正消除地表起伏和低降速帶的變化對反射波旅行時的影響。 高質(zhì)量的動靜校正是獲取最佳疊加剖面的基礎(chǔ)。高質(zhì)量的動靜校正是獲取最佳疊加剖面的基礎(chǔ)。 速度是疊加的關(guān)鍵參數(shù)。關(guān)

2、于速度的測量方法：聲波測井的速度是疊加的關(guān)鍵參數(shù)。關(guān)于速度的測量方法：聲波測井的直接測量法；地震勘探數(shù)據(jù)的間接測量法。直接測量法；地震勘探數(shù)據(jù)的間接測量法。 地震勘探中有關(guān)速度的概念：層速度、平均速度、均方根地震勘探中有關(guān)速度的概念：層速度、平均速度、均方根（rms）速度、瞬時速度、相速度、群速度、動校正速度、瞬時速度、相速度、群速度、動校正（NMO）速度、疊加速度和偏移速度等。速度、疊加速度和偏移速度等。 本章討論根據(jù)地震數(shù)據(jù)來估測速度的方法。 估測速度需要共本章討論根據(jù)地震數(shù)據(jù)來估測速度的方法。 估測速度需要共中心點中心點（CMP）記錄所提供的非零炮檢距數(shù)據(jù)。利用估測出）記錄所提供的非零炮

3、檢距數(shù)據(jù)。利用估測出來的速度校正非零炮檢距時差，把記錄到的數(shù)據(jù)體（在中心點來的速度校正非零炮檢距時差，把記錄到的數(shù)據(jù)體（在中心點-炮檢距炮檢距-時間坐標(biāo)中）壓縮為疊加剖面。時間坐標(biāo)中）壓縮為疊加剖面。 一動校正速度一動校正速度 定義：為校正正常時差所用的速度稱為動校正速度。定義：為校正正常時差所用的速度稱為動校正速度。 1單個水平反射層：單個水平反射層：NMO 速度等于該反射層上部介質(zhì)的速度。速度等于該反射層上部介質(zhì)的速度。 2單個傾斜反射層：單個傾斜反射層：NMO 速度等于該反射層上部介質(zhì)速度除以反射速度等于該反射層上部介質(zhì)速度除以反射層傾角的余弦。若考慮三維空間傾斜反射層，還需考慮方位角因

4、素。層傾角的余弦。若考慮三維空間傾斜反射層，還需考慮方位角因素。 3多層水平反射層：小炮檢距時，某個水平反射層的多層水平反射層：小炮檢距時，某個水平反射層的 NMO 速度等于該速度等于該反射層上覆介質(zhì)的反射層上覆介質(zhì)的 rms 速度。速度。 4多層任意傾斜反射層：只要傾角不大，分布不廣，仍可用雙曲線近似。多層任意傾斜反射層：只要傾角不大，分布不廣，仍可用雙曲線近似。 NMO 速度與疊加速度的差別：速度與疊加速度的差別：NMO 速度是依據(jù)小排列雙曲線形狀分布速度是依據(jù)小排列雙曲線形狀分布旅行時間（旅行時間（Taner 和和 Koehler,1969; Al-Chalabi,1973） ；而疊加速

5、度則是依照；而疊加速度則是依照與整個排列長度數(shù)據(jù)擬合最好的雙曲線。與整個排列長度數(shù)據(jù)擬合最好的雙曲線。 但但是是在實際在實際應(yīng)用應(yīng)用中中通常認(rèn)為這兩種速度是相等的。通常認(rèn)為這兩種速度是相等的。 二速度分析方法二速度分析方法 建立在雙曲線假設(shè)基礎(chǔ)之上的常規(guī)速度分析方法：建立在雙曲線假設(shè)基礎(chǔ)之上的常規(guī)速度分析方法： 122xt 法法 22xt 平面上的反射波雙曲線時距方程為線性方程。 因此， 從平面上的反射波雙曲線時距方程為線性方程。 因此， 從22xt 坐標(biāo)中的最佳擬合直線可估計出零炮檢距上的反射波時間和該反射的疊加坐標(biāo)中的最佳擬合直線可估計出零炮檢距上的反射波時間和該反射的疊加速度。速度。 2

6、速度掃描法速度掃描法 該方法是應(yīng)用一系列常速度值在該方法是應(yīng)用一系列常速度值在 CMP 道集作動校正，并將結(jié)果并列道集作動校正，并將結(jié)果并列顯示，從中選出能使反射波同相軸拉平程度最高的速度作為顯示，從中選出能使反射波同相軸拉平程度最高的速度作為 NMO 速度。速度。 3常速疊加（常速疊加（CVS）法）法 取測線的一小段，用一系列常速度值作疊加處理，不同的速度疊加成取測線的一小段，用一系列常速度值作疊加處理，不同的速度疊加成不同的疊加圖象，稱為不同的疊加圖象，稱為 CVS 圖象。從圖象。從 CVS 圖象中取出獲得最佳疊加的速圖象中取出獲得最佳疊加的速度為疊加速度。度為疊加速度。 三靜校正三靜校正

7、 地表起伏和地表起伏和/或近地表速度的變化所造成的靜態(tài)時移會使水平層狀介或近地表速度的變化所造成的靜態(tài)時移會使水平層狀介質(zhì)中的反射波時距曲線偏離雙曲線。對這種時移所作的校正稱為靜校正。質(zhì)中的反射波時距曲線偏離雙曲線。對這種時移所作的校正稱為靜校正。 1野外靜校正：在野外對估計出來的風(fēng)化層和高程變化所做的初步野外靜校正：在野外對估計出來的風(fēng)化層和高程變化所做的初步校正稱為野外靜校正。校正稱為野外靜校正。 2剩余靜校正：野外靜校正后，在地震數(shù)據(jù)中仍然殘留有各種剩余剩余靜校正：野外靜校正后，在地震數(shù)據(jù)中仍然殘留有各種剩余靜態(tài)時移， 通常在疊前必須估計出這類剩余靜態(tài)時移值， 并在靜態(tài)時移， 通常在疊前

8、必須估計出這類剩余靜態(tài)時移值， 并在 CMP 道集道集中加以校正。這種校正稱為剩余靜校正。中加以校正。這種校正稱為剩余靜校正。 區(qū)域速度函數(shù)或速度分析區(qū)域速度函數(shù)或速度分析初步初步 NMO 動校動校估計剩余靜校正量，估計剩余靜校正量，進(jìn)行剩余靜校正進(jìn)行剩余靜校正重新進(jìn)行速度分析（提高所拾重新進(jìn)行速度分析（提高所拾取的速度質(zhì)量）取的速度質(zhì)量）疊疊加加剩余靜校正剩余靜校正。這是一個多次迭代過程。這是一個多次迭代過程。 四水平疊加四水平疊加 水平疊加是將水平疊加是將 CMP 道集記錄經(jīng)道集記錄經(jīng) NMO 動校后疊加起來， 目的是壓制動校后疊加起來， 目的是壓制隨機噪音，提高地震信噪比。隨機噪音，提高

9、地震信噪比。 注意：疊加和偏移所要求的速度未必相同。注意：疊加和偏移所要求的速度未必相同。 疊加速度與傾角有關(guān)，而偏移速度與傾角無關(guān)。疊加速度與傾角有關(guān)，而偏移速度與傾角無關(guān)。 單個水平地層的動校正單個水平地層的動校正 共中心點道集共中心點道集地震波旅行時方程地震波旅行時方程為為（Pythagoras 理論理論） ： 2222)0()(vxtxt （4.2.1） 式中，式中，x 是震源與接收點之間的距離（偏移距） ，是震源與接收點之間的距離（偏移距） ，v 是反射界是反射界面以上介質(zhì)的速度，而面以上介質(zhì)的速度，而 t（0）是沿垂直路徑是沿垂直路徑 MD 的雙程旅行的雙程旅行時。時。注意深度點對

10、地面的投影，沿著反射層正交線，與中心注意深度點對地面的投影，沿著反射層正交線，與中心點點 M 重合。重合。這只有當(dāng)反射層是水平的情況才如此。這只有當(dāng)反射層是水平的情況才如此。 右圖右圖2是一個共中心點道集（是一個共中心點道集（CMP），），也代表一個共深度點道集（也代表一個共深度點道集（CDP）（在這個（在這個CMP道集中的所有道包含道集中的所有道包含來自同一深度點的反射。）圖中偏移來自同一深度點的反射。）圖中偏移距范圍為距范圍為03150m，道距，道距50m。反射。反射層以上的介質(zhì)速度為層以上的介質(zhì)速度為2264m/s。在給定偏移距上的雙程旅行時在給定偏移距上的雙程旅行時 t（x）與零偏雙程

11、時之間的差稱作動態(tài)時差與零偏雙程時之間的差稱作動態(tài)時差 NMO。 NMO 速度一旦估算出來，炮檢距對波至?xí)r間的影響就能校正。速度一旦估算出來，炮檢距對波至?xí)r間的影響就能校正。 雙曲線時移校正雙曲線時移校正的的數(shù)值方法數(shù)值方法：根據(jù)：根據(jù)原原始始 CMP 道集中道集中 A 的振幅值找出動校后道集上的振幅值找出動校后道集上A的振幅值。的振幅值。 給定給定 t（0） 、x 和和NMOv值根據(jù)方程（值根據(jù)方程（4.2.1）算出）算出 t（x） 。假定是假定是 1003ms，如果采樣間，如果采樣間隔為隔為 4ms，那么該時間就等于第，那么該時間就等于第 250.75 個采樣點。個采樣點。因此，必須采取相

12、鄰的整數(shù)序樣點因此，必須采取相鄰的整數(shù)序樣點上的振幅值通過內(nèi)插或上的振幅值通過內(nèi)插或抽樣定理來計算該時刻的振幅。抽樣定理來計算該時刻的振幅。 NMO 動校值的計算：動校值的計算： 1)0(1)0()0()(212tvxttxttNMONMO 表表4-1 不同炮檢距不同炮檢距 x的的NMO值和已知速度的零炮檢距雙程時間值和已知速度的零炮檢距雙程時間 stNMO, st),0( smvNMO/, mx,1000 mx,2000 0.25 0.5 1 2 4 2000 2500 3000 3500 4000 0.309 0.140 0.054 0.020 0.008 0.780 0.443 0.20

13、1 0.080 0.031 不同速度對均勻介質(zhì)水平面反射的動校正：不同速度對均勻介質(zhì)水平面反射的動校正： 如果所用速度高于介質(zhì)速度（如果所用速度高于介質(zhì)速度（2264m/s） ，雙曲線不能完全拉平，稱為欠校正。） ，雙曲線不能完全拉平，稱為欠校正。 所用速度低于介質(zhì)速度，雙曲線上翹，稱為過校正。所用速度低于介質(zhì)速度，雙曲線上翹，稱為過校正。 傳統(tǒng)速度分析的基礎(chǔ)傳統(tǒng)速度分析的基礎(chǔ):采用方程（采用方程（4.2.2）對）對 CMP 道集通過一系列常速度進(jìn)行道集通過一系列常速度進(jìn)行動校試驗，使該道集的反射曲線拉得最平的速度就是疊前最佳動校正速度。動校試驗，使該道集的反射曲線拉得最平的速度就是疊前最佳動

14、校正速度。 1)0(1)0()0()(212tvxttxttNMONMO二水平層狀介質(zhì)的動校正二水平層狀介質(zhì)的動校正 對于常速層狀介質(zhì)， 地震射線從震源對于常速層狀介質(zhì)， 地震射線從震源 S 至深度點至深度點 D 然后返回接收點然后返回接收點 R，地面，地面中點在中點在 M, 炮檢距為炮檢距為 x。旅行時方程可表示為（。旅行時方程可表示為（Taner，Koehler，1969） ：） ： 63422102)(xCxCxCCxt （4.2.3） 式中式中322120,1),0(CCvCtCrms 是地層厚度和層速度的復(fù)雜函數(shù)。是地層厚度和層速度的復(fù)雜函數(shù)。深度點深度點D 的均方根速度定義為：的均

15、方根速度定義為：Niiirmstvtv122)0()0(1（4.2.4） 此處此處it為第為第 i 層的雙程旅行時間：層的雙程旅行時間：ikktt1)0(。 若排列近似為小排列（炮檢距小于深度）若排列近似為小排列（炮檢距小于深度） 則方程（則方程（4.2.3）中的級數(shù)可省略為：）中的級數(shù)可省略為： 2222)0()(rmsvxtxt （4.2.5） 比較（比較（4.2.1）和）和（4.2.5） ，可見對于水平層狀介質(zhì)，若小排列近似關(guān)系成立，可見對于水平層狀介質(zhì)，若小排列近似關(guān)系成立，NMO 動校速度等于均方根速度。動校速度等于均方根速度。 淺層同相軸淺層同相軸t（0）=0.8s和和t（0）=1

16、.2s只是在大炮檢距處有少許差別，通過省略只是在大炮檢距處有少許差別，通過省略高階項，可用小排列雙曲線近似求出水平層反射時間。高階項，可用小排列雙曲線近似求出水平層反射時間。三三動動校校拉拉伸伸動動校校正正結(jié)結(jié)果果出出現(xiàn)現(xiàn)頻頻率率畸畸變變，同同相相軸軸移移向向低低頻頻。 主主周周期期為為T的的波波形形經(jīng)經(jīng)NMO動動校校之之后后拉拉伸伸為為T，拉拉伸伸量量為為： )0( ttffNMO （4.2.6） 其其中中f 是是主主頻頻，f是是所所引引起起的的頻頻率率變變化化，NMOt由由方方程程（4.2.2）給給出出。 影影響響：大大炮炮檢檢距距上上波波形形拉拉伸伸將將嚴(yán)嚴(yán)重重?fù)p損害害淺淺層層同同相相軸

17、軸疊疊加加效效果果 解解決決辦辦法法：切切除除 表表 4-2 NMO 拉伸拉伸 ff% st),0( smvNMO/, mx1000 mx2000 0.25 0.5 1 2 4 2000 2500 3000 3500 4000 123 28 5 1 0.2 312 89 20 4 0.8 表表 4-2 列出用表列出用表 4-1 中速度關(guān)系所導(dǎo)致的動校拉伸，用頻率變化百分比表示。中速度關(guān)系所導(dǎo)致的動校拉伸，用頻率變化百分比表示?？梢娎熘饕抻跍\層和大炮檢距，譬可見拉伸主要限于淺層和大炮檢距，譬如一個如一個st25. 0)0(，主頻，主頻 30Hz，炮，炮檢距檢距 2000m 處的同相軸，經(jīng)處的

18、同相軸，經(jīng) NMO 動校后，主頻移到了近動校后，主頻移到了近 10Hz。 切除切除量量：拉伸量拉伸量 50%100%以上以上，在不造成質(zhì)量下降的情況下盡量多保留在不造成質(zhì)量下降的情況下盡量多保留CMP 道參與疊加，通常要兼顧信噪比（道參與疊加，通常要兼顧信噪比（S/N）和切除，采取折中）和切除，采取折中。 未切除時，未切除時，CMP的淺層部分可見一個低頻的拉伸帶的淺層部分可見一個低頻的拉伸帶 最佳切除選擇法：采用逐步疊加，根最佳切除選擇法：采用逐步疊加，根據(jù)波形變化選擇切除的方法。疊加次據(jù)波形變化選擇切除的方法。疊加次數(shù)是向近炮點方向逐步增加。數(shù)是向近炮點方向逐步增加。 切除太多是危險的，因為

19、大炮檢距是切除太多是危險的，因為大炮檢距是有效壓制多次波所必需的數(shù)據(jù)。有效壓制多次波所必需的數(shù)據(jù)。四單個傾斜地層的動校正四單個傾斜地層的動校正對于傾斜層，中點對于傾斜層，中點M不再是不再是深度點深度點D在地表的投影。在地表的投影。CDP道集和道集和CMP道集只有在水平層狀地層時才等價，在地道集只有在水平層狀地層時才等價，在地下界面傾斜或速度橫向變化時，這兩種道下界面傾斜或速度橫向變化時，這兩種道集不相同。集不相同。對于道集中的炮對于道集中的炮-檢對來說，檢對來說，不論界面是否傾斜，中點不論界面是否傾斜，中點M總是共中心點，總是共中心點，但所記錄的傾斜反射層但所記錄的傾斜反射層CMP道集中的每道

20、集中的每一炮一炮-檢對，它所反映的地下深度點檢對，它所反映的地下深度點D則不則不相同。相同。Levin（1971）由圖）由圖4-13的的幾何關(guān)系導(dǎo)出具有幾何關(guān)系導(dǎo)出具有傾角地層的時間方程傾角地層的時間方程如下：如下： （4.2.7）22222cos)0()(vxtxt傾斜地層反射波同相軸作正確疊加所要求的動校速度傾斜地層反射波同相軸作正確疊加所要求的動校速度比其上覆介質(zhì)速度大比其上覆介質(zhì)速度大 。cosvvNMO四單個傾斜地層的動校正四單個傾斜地層的動校正Levin把該結(jié)論推廣到三維（把該結(jié)論推廣到三維（3-D）空間傾斜界面，）空間傾斜界面，NMO速度不僅依速度不僅依賴界面傾角，而且依賴炮賴界

21、面傾角，而且依賴炮-檢布排方位：檢布排方位： （4.2.9）方位角是實際剖面方向與構(gòu)造傾向的夾角（圖方位角是實際剖面方向與構(gòu)造傾向的夾角（圖4-14），視傾角的定義），視傾角的定義為：為：（4.2.10） 由此定義重寫（由此定義重寫（4.2.8）式的）式的NMO速度：速度： （4.2.11）該式與適用二維（該式與適用二維（2-D）界面幾何關(guān)系的（）界面幾何關(guān)系的（4.2.8）式形式相同，但式）式形式相同，但式（4.2.8）中所用的是真傾角，而（）中所用的是真傾角，而（4.2.11）式中所用的是視傾角。）式中所用的是視傾角。 2122)cossin1 ( vvNMOcossinsincosvvN

22、MO當(dāng)傾角不超過當(dāng)傾角不超過16o時，時，速度比幾乎不變；傾角速度比幾乎不變；傾角為為16o時，疊加速度與時，疊加速度與真實速度相差真實速度相差4%。 總之，不論總之，不論2-D或或3-D，傾斜層的，傾斜層的NMO速度速度跟傾角有關(guān)。高速水平地層跟低速傾斜層在跟傾角有關(guān)。高速水平地層跟低速傾斜層在炮檢距時差上可以完全相同炮檢距時差上可以完全相同 l五任意傾斜層狀介質(zhì)的動校正五任意傾斜層狀介質(zhì)的動校正 4.2.4 4.2.54.2.8圖圖示示幾幾個個任任意意傾傾斜斜地地層層的的 2-D 界界面面反反射射路路徑徑幾幾何何關(guān)關(guān)系系: 應(yīng)應(yīng)看看到到 CMP 射射線線是是由由中中點點 M 以以該該傾傾斜

23、斜界界面面的的法法射射入入射射到到D，不不是是D點點。零零炮炮檢檢距距時時間間是是 M 到到D的的雙雙程程時時間間。Hubral 和和 kery（1980）求求出出沿沿 SDG 路路程程旅旅行行時時間間： 高高階階項項 2222)0()(NMOvxtxt （4.2.12） NiikkkiiNMOtvtv1112202)coscos()0(cos)0(1 (4.2.13) 式式中中角角度度由由圖圖 4-18 定定義義，對對單單一一傾傾斜斜層層，方方程程（4.2.13）簡簡化化為為方方程程（4.2.8） 。進(jìn)進(jìn)而而，對對于于水水平平層層狀狀地地層層，方方程程（4.2.13）又又可可簡簡化化為為方方

24、程程（4.2.4） 。只只要要傾傾斜斜平平緩緩，而而且且是是小小排排列列，就就能能利利用用雙雙曲曲線線近近似似表表示示旅旅行行時時間間方方程程方方程程（4.2.5） ，并并且且所所需需的的 NMO 動動校校速速度度近近似似等等于于均均方方根根速速度度方方程程（4.2.4）。 Niiirmstvtv122)0()0(1 cosvvNMO 2222)0()(rmsvxtxt 表表 4-3 各種地層模型的各種地層模型的 NMO 速度速度 模模 型型 NMO 速度速度 單個水平層單個水平層 反射界面以上介質(zhì)的速度反射界面以上介質(zhì)的速度 水平層狀地層水平層狀地層 由小排列給出的均方根速度函數(shù)由小排列給出

25、的均方根速度函數(shù) 單個傾斜層單個傾斜層 界面以上介質(zhì)速度被除以傾角的余弦界面以上介質(zhì)速度被除以傾角的余弦 多個任意傾斜層多個任意傾斜層 由小排列、小傾角定義的均方根速度函由小排列、小傾角定義的均方根速度函數(shù)數(shù) 在小排列和小傾角規(guī)定下，它們的在小排列和小傾角規(guī)定下，它們的 NMO 時間都近時間都近似符合雙曲線，并由下式確定：似符合雙曲線，并由下式確定： 2222)0()(NMOvxtxt （4.2.14） 雙曲線的雙曲線的 NMO 速度應(yīng)同疊加速度不同，疊加速度是速度應(yīng)同疊加速度不同，疊加速度是能對能對 CMP 道集作最佳疊加的速度，利用雙曲線形狀來道集作最佳疊加的速度，利用雙曲線形狀來確定出最

26、佳確定出最佳的疊加軌跡：的疊加軌跡： 2222)0()(stststvxtxt （4.2.15） 這里的這里的stv是該排列范圍上與是該排列范圍上與 CMP 道集時間曲線最接近的某條雙曲線的速度， 該雙道集時間曲線最接近的某條雙曲線的速度， 該雙曲線未必是方程（曲線未必是方程（4.2.14）所定義的小排列雙曲線。）所定義的小排列雙曲線。 這類疊加速度這類疊加速度stv與與 NMO 動校速度動校速度NMOv的差別稱為排列長度偏離（的差別稱為排列長度偏離（Alchalabi,1973; Hubral 和和 Krey,1980） 。從方程（） 。從方程（4.2.14）和（）和（4.2.15）可）可見

27、，排列越短，最佳疊加雙曲見，排列越短，最佳疊加雙曲一22xt 法法 動校時間差是由地震數(shù)據(jù)確定速度的基礎(chǔ)。用所得速度作動校正，使動校時間差是由地震數(shù)據(jù)確定速度的基礎(chǔ)。用所得速度作動校正，使 CMP 道集在疊前對齊。道集在疊前對齊。 方程（方程（4.2.15） 在在22)(xxt平面上描述出一條直線，直線的斜率為平面上描述出一條直線，直線的斜率為vst，截距為，截距為)0( t。 表表 4-4 由圖由圖 4-20 合成模型所估算的疊加速合成模型所估算的疊加速度與實測均方根速度度與實測均方根速度 st),0( 疊加速度疊加速度smxt/),(22 實測均方根實測均方根速度速度,sm / 0.4 0

28、.8 1.2 1.6 2000 2264 2519 2828 2000 2264 2533 2806 212122121122212212221111221222021)()(t 22-4kkvkxtvtvkvtvtkxvktxvktxtkvxvxtkvxtnnnRRRRRRR下層：上層：切線斜率方程：的證明關(guān)于二速度掃描二速度掃描 對對 CMP 道集用常速掃描是另一種速度分析道集用常速掃描是另一種速度分析技術(shù)。技術(shù)。 圖圖 4-23 有一個有一個 CMP 道集， 用一系列常速道集， 用一系列常速度，從度，從 5000ft/s 到到 13600ft/s，重復(fù)對道集作，重復(fù)對道集作 NMO 校正

29、，每校正一次得一張圖象，并把它們校正，每校正一次得一張圖象，并把它們并置在一起。并置在一起?，F(xiàn)在來考察同相軸現(xiàn)在來考察同相軸 A 的的 NMO，看到采用小速度時它獲得過校正，采用大速度，看到采用小速度時它獲得過校正，采用大速度時則獲得欠校正，在采用時則獲得欠校正，在采用 8300ft/s 時，時，NMO 校正使同相軸為水平，因此校正使同相軸為水平，因此 8300ft/s 為同相軸為同相軸A 的疊加速度。的疊加速度。同相軸同相軸 B 在相應(yīng)在相應(yīng) 8900ft/s 的速度時校正到水平。的速度時校正到水平。通過這種方法，我們求出通過這種方法，我們求出適合于該道集的適合于該道集的 NMO 動校速度函

30、數(shù)。動校速度函數(shù)。 三常速疊加法三常速疊加法(CVS) 主要原理：根據(jù)一系列常速度形成的疊加數(shù)據(jù)中疊加同相軸的振幅和連續(xù)性來估測疊加速主要原理：根據(jù)一系列常速度形成的疊加數(shù)據(jù)中疊加同相軸的振幅和連續(xù)性來估測疊加速度。度。 實現(xiàn)方法：選擇實現(xiàn)方法：選擇 24 個個 CMP 道集（典型的范圍是道集（典型的范圍是 2448CMP, 但也有整條測線的） ，給定某但也有整條測線的） ，給定某一疊加速度，得到一疊加速度，得到 24 個經(jīng)個經(jīng) NMO 校正后的疊加道，形成一段常速疊加剖面。改變速度得到校正后的疊加道，形成一段常速疊加剖面。改變速度得到一系列疊加段，排列起來形成常速度疊加圖象。按時間對所研究同

31、相軸產(chǎn)生最佳疊加響應(yīng)的一系列疊加段，排列起來形成常速度疊加圖象。按時間對所研究同相軸產(chǎn)生最佳疊加響應(yīng)的原則選擇速度。原則選擇速度。 注意：深部同相軸速度估計分辨率有所降低，原因是注意：深部同相軸速度估計分辨率有所降低，原因是 NMO 時差隨深度迅速減小。時差隨深度迅速減小。 四速度譜四速度譜 如圖所示： 從某一速度例如如圖所示： 從某一速度例如 2000m/s 到某一速度例如到某一速度例如 4300m/s 的各種速度反復(fù)對的各種速度反復(fù)對 CMP道集進(jìn)行道集進(jìn)行 NMO 校正和疊加，把每一種速度所得的疊加結(jié)果并排顯示在速度校正和疊加，把每一種速度所得的疊加結(jié)果并排顯示在速度-雙程零炮檢距雙程零

32、炮檢距時間平面中，稱此為速度譜（時間平面中，稱此為速度譜（Taner 和和 Koehler，1969） 。） 。 在速度為在速度為3000m/s時獲得了最大時獲得了最大疊加振幅，該速度應(yīng)是輸入疊加振幅，該速度應(yīng)是輸入CMP道集中該同相軸的疊加速道集中該同相軸的疊加速度。速度譜上的那些低振幅水度。速度譜上的那些低振幅水平軌跡是小炮檢距分量的疊加，平軌跡是小炮檢距分量的疊加，高振幅區(qū)域是全部炮檢距分量高振幅區(qū)域是全部炮檢距分量疊加結(jié)果（疊加結(jié)果（Sherwood和和Poe,1972）,因此必須保留長炮因此必須保留長炮檢距數(shù)據(jù)以保證分辨率。檢距數(shù)據(jù)以保證分辨率。 2700、2800和和3000m/s

33、分別對應(yīng)分別對應(yīng)淺、中、深三種同相軸。速度譜淺、中、深三種同相軸。速度譜不僅能提供疊加速度不僅能提供疊加速度,而且可以區(qū)而且可以區(qū)分反射波和多次波。分反射波和多次波。 速度分析的目的是通過雙曲線軌跡對整個速度分析的目的是通過雙曲線軌跡對整個 CMP 道集選出使信號具有最佳相干性的速度。道集選出使信號具有最佳相干性的速度。 疊加相干性除了疊加振幅（前面的例子）表示外，還可用其他不同的相干性測量方法：疊加相干性除了疊加振幅（前面的例子）表示外，還可用其他不同的相干性測量方法： l 疊后振幅：疊后振幅： Miititfs1)(, （4.3.1） )(,itif為第為第 i 道雙程時間為道雙程時間為)

34、(it的振幅，的振幅，M 為為 CMP 道集總道數(shù)，道集總道數(shù)，)(it由旅行時雙曲線確定：由旅行時雙曲線確定： 21222)0()(stivxtit （4.3.2） 定義歸一化疊加振幅：定義歸一化疊加振幅： iititfsNS)(, （4.3.3） NS的范圍是的范圍是10 NS。 l 疊疊后后振振幅幅： l 互互相相關(guān)關(guān)求求和和： 非非歸歸一一化化互互相相關(guān)關(guān)求求和和：沿沿著著試試驗驗疊疊加加雙雙曲曲線線采采用用一一個個時時窗窗對對CMP道道集集非非歸歸一一化化互互相相關(guān)關(guān)求求和和 tMiitittMiitiMiitifsffCC12)(,212)(,21)(,2121 （4.3.4） C

35、C可可解解釋釋為為輸輸出出疊疊加加能能量量和和輸輸入入能能量量之之差差的的一一半半。 歸歸一一化化互互相相關(guān)關(guān)求求和和：CC的的歸歸一一化化形形式式是是速速度度譜譜計計算算中中常常用用的的另另一一屬屬性性： tMkkMitkitkititikitkiitiffffMMNCC1112)(,2)(,)(,)(,21) 1(2 l 歸一化的能量互相關(guān)求和歸一化的能量互相關(guān)求和 歸一化的能量互相關(guān)求和表達(dá)式為：歸一化的能量互相關(guān)求和表達(dá)式為： tMiitifCCMECC12)(,)1(2 （4.3.6） ECC的范圍是的范圍是 111 ECCM。 歸一化的輸入歸一化的輸入/輸出能量比，即相似性為：輸出

36、能量比，即相似性為： tMiitittfsMNE12)(,21 （4.3.7） NE和和ECC的關(guān)系：的關(guān)系：)1(11 NEMMECC （4.3.8） NE的范圍是的范圍是10 NE。 表表 4-5 由式（由式（4.3.9）表示的二次疊加特例的相干屬性測量值）表示的二次疊加特例的相干屬性測量值 屬性屬性 5 . 0a 5 . 0a 疊加振幅疊加振幅式（式（4.3.1） )(5 . 1tf )(5 . 0tf NS 相干性相干性式（式（4.3.3） 1 0.333 非歸一化互相關(guān)求和非歸一化互相關(guān)求和 CC式（式（4.3.4） ttf)(5 . 02 ttf)(5 . 02 歸一化互相關(guān)求和歸

37、一化互相關(guān)求和 NCC式（式（4.3.5） 1 1 能量歸一化互相關(guān)求和能量歸一化互相關(guān)求和 ECC式 （式 （4.3.6） 0.8 -0.8 相似性相似性 NE式（式（4.3.7） 0.9 0.1 疊加振幅對極性很敏感，非歸一化互相關(guān)在速度譜上起到突出強反疊加振幅對極性很敏感，非歸一化互相關(guān)在速度譜上起到突出強反射的作用，而歸一化互相關(guān)或歸一化能量互相關(guān)則起到加強速度譜射的作用，而歸一化互相關(guān)或歸一化能量互相關(guān)則起到加強速度譜中的弱反射。正如式（中的弱反射。正如式（4.3.8）所示，相似性跟偏置的歸一化能量互）所示，相似性跟偏置的歸一化能量互相關(guān)求和相差無幾。相關(guān)求和相差無幾。速度譜通常速度

38、譜通常用用圖圖 4-28 中的兩種顯示形式，由時窗排列圖或等值線圖形式選取速度。中的兩種顯示形式，由時窗排列圖或等值線圖形式選取速度。雖然時窗排列圖也用的較普遍，但多數(shù)人偏愛使用等值線圖。雖然時窗排列圖也用的較普遍，但多數(shù)人偏愛使用等值線圖。 另一種有助于拾取速度的量值是每一時窗的最大相干值，將它以時間函數(shù)形式并排另一種有助于拾取速度的量值是每一時窗的最大相干值，將它以時間函數(shù)形式并排在速度譜邊上，如圖在速度譜邊上，如圖 4-28 所示。所示。 五影響速度估算的因素五影響速度估算的因素 下述因素會限制地震資料速度估算的精度和分辨率下述因素會限制地震資料速度估算的精度和分辨率伊爾馬滋伊爾馬滋（Y

39、ilmaz）,1993： （1）排列長度：缺乏大炮檢距信息意味著缺乏辨別速度所需要的重要時差；但大炮檢）排列長度：缺乏大炮檢距信息意味著缺乏辨別速度所需要的重要時差；但大炮檢距區(qū)域的資料有拉伸問題。距區(qū)域的資料有拉伸問題。 （2）疊加次數(shù)：疊加）疊加次數(shù)：疊加 32 次甚至次甚至 16 次對速度譜沒有影響，但再低使峰值發(fā)生嚴(yán)重位移。次對速度譜沒有影響，但再低使峰值發(fā)生嚴(yán)重位移。 （3）S/N 比：存在高幅隨機噪音時也可識別有效信號，但比：存在高幅隨機噪音時也可識別有效信號，但 S/N 不高時精度要受限制。不高時精度要受限制。 （4）切除：會減少淺層疊加次數(shù)）切除：會減少淺層疊加次數(shù),導(dǎo)致切除帶

40、位置的同相軸振幅減弱，它對速度譜有副導(dǎo)致切除帶位置的同相軸振幅減弱，它對速度譜有副作用；校正方法是用切除帶中有效疊加次數(shù)比例乘疊加振幅來實現(xiàn)作用；校正方法是用切除帶中有效疊加次數(shù)比例乘疊加振幅來實現(xiàn)。 （5）時窗寬度：太小，工作量大；太大，缺乏時間分辨率。一般為信號主周期的一半）時窗寬度：太小，工作量大；太大，缺乏時間分辨率。一般為信號主周期的一半到一倍之間，約為到一倍之間，約為 20 到到 40ms。注意淺層周期短，深層周期長。注意淺層周期短，深層周期長。 （6）速度采樣密度：掃描范圍應(yīng)包含一次反射波速度；速度間隔太大會降低分辨率。）速度采樣密度：掃描范圍應(yīng)包含一次反射波速度；速度間隔太大會

41、降低分辨率。 （7）相干屬性量的選擇；）相干屬性量的選擇； （8）對雙曲性正常時差的偏離度；）對雙曲性正常時差的偏離度； （9）數(shù)據(jù)的頻譜寬度。）數(shù)據(jù)的頻譜寬度。 改改善善速速度度譜譜分分析析的的方方法法： l 數(shù)數(shù)據(jù)據(jù)重重新新采采樣樣（重重新新按按時時間間分分樣樣） 。 l 帶帶通通濾濾波波和和自自動動增增益益控控制制（AGC）改改善善互互相相關(guān)關(guān)處處理理。 l 在在分分析析時時利利用用幾幾個個相相鄰鄰 CMP 道道集集：一一種種是是道道集集相相加加求求和和，再再計計算算速速度度譜譜；另另一一種種是是先先算算出出每每一一個個道道集集的的速速度度譜譜，再再進(jìn)進(jìn)行行譜譜相相加加。前前者者節(jié)節(jié)省省

42、機機時時但但未未考考慮慮構(gòu)構(gòu)造造傾傾角角影影響響。 l 速速度度譜譜矩矩陣陣作作平平滑滑。 l 工工區(qū)區(qū)的的工工作作經(jīng)經(jīng)驗驗有有助助于于選選擇擇正正確確的的一一次次反反射射波波疊疊加加速速度度。 六層速度分析六層速度分析 沿著某個有用層次連續(xù)追蹤分析速度，稱為層速度分析（沿著某個有用層次連續(xù)追蹤分析速度，稱為層速度分析（HVA） 。） 。 基本原理基本原理:沿雙曲線分布的時窗算出相干值， 表示出各沿雙曲線分布的時窗算出相干值， 表示出各 CMP 位置的位置的速度函數(shù)。疊加剖面上有構(gòu)造間斷的地方，則按斷層分開的區(qū)段來作速度函數(shù)。疊加剖面上有構(gòu)造間斷的地方，則按斷層分開的區(qū)段來作HVA。 層速度分

43、析是沿著所選定的關(guān)鍵層逐個提取層速度分析是沿著所選定的關(guān)鍵層逐個提取 CMP 位置上速度信息位置上速度信息的有效方法；的有效方法；能改善疊加剖面的質(zhì)量，以備作疊后深度偏移。能改善疊加剖面的質(zhì)量，以備作疊后深度偏移。 注意：注意：HVA 建立在雙曲線動校原理之上，而需作深度偏移的資料建立在雙曲線動校原理之上，而需作深度偏移的資料往往具有復(fù)雜的動校形態(tài)，因而往往具有復(fù)雜的動校形態(tài)，因而 HVA 受到影響。但無論如何，實際受到影響。但無論如何，實際中中 HVA 確實能提供出標(biāo)志層詳細(xì)的速度橫向變化情況。確實能提供出標(biāo)志層詳細(xì)的速度橫向變化情況。 一、定義：一、定義： 消除近地表不規(guī)則所造成時差， 使

44、反射波時距曲線恢復(fù)為雙曲線的過程稱為靜校正。消除近地表不規(guī)則所造成時差， 使反射波時距曲線恢復(fù)為雙曲線的過程稱為靜校正。 野外靜校往往不包括地形突變、風(fēng)化基底和風(fēng)化速度的影響，這些影響必需通過剩野外靜校往往不包括地形突變、風(fēng)化基底和風(fēng)化速度的影響，這些影響必需通過剩余靜校正加以消除。余靜校正加以消除。 二、剩余靜校正量的影響：二、剩余靜校正量的影響： 嚴(yán)重降低速度譜分析的質(zhì)量。嚴(yán)重降低速度譜分析的質(zhì)量。 導(dǎo)致錯誤的疊加剖面，形成暗點和假構(gòu)造。導(dǎo)致錯誤的疊加剖面，形成暗點和假構(gòu)造。 一般在野外靜校正的基礎(chǔ)上作NMO，再作剩余靜校正，然后重做速度分析，更新速度提取。剩余靜校正通常用于陸上資料，但在

45、特定情況下對海上資料也產(chǎn)生重大改進(jìn)。一自動統(tǒng)計靜校正方法（靜校正思路一）一自動統(tǒng)計靜校正方法（靜校正思路一） 1 1假設(shè)條件和特點假設(shè)條件和特點 1 1）假設(shè)：假設(shè)： （1 1）同一炮點在低速帶中入射的時間與入射角無關(guān)，即可認(rèn)為在低速同一炮點在低速帶中入射的時間與入射角無關(guān)，即可認(rèn)為在低速帶中都是垂直入射的。帶中都是垂直入射的。 （2 2）炮點（或接收點）由于地形起伏及低速帶變化所引起的靜校正量時）炮點（或接收點）由于地形起伏及低速帶變化所引起的靜校正量時差是隨機的，其均值為零。差是隨機的，其均值為零。 以點假設(shè)在一定條件下符合實際以點假設(shè)在一定條件下符合實際 情況的。例如，在低速帶劇烈變化的

46、地區(qū)。情況的。例如，在低速帶劇烈變化的地區(qū)。 采用某種方法一定要考慮到實際條件。采用某種方法一定要考慮到實際條件。 tx基準(zhǔn)面或準(zhǔn)基準(zhǔn)面2 2）幾個特點幾個特點 （1 1）某個記錄道的靜校正量包括激發(fā)點靜校正量和接收點靜校正量。設(shè)）某個記錄道的靜校正量包括激發(fā)點靜校正量和接收點靜校正量。設(shè)某道的靜校正量是某道的靜校正量是，激發(fā)點和接收點靜校正量分別是，激發(fā)點和接收點靜校正量分別是S S和和R R ，則有，則有 = =S S+ +R R （4.4.1） 如果提取靜校正量時是求絕對靜校正量，要分別求出如果提取靜校正量時是求絕對靜校正量，要分別求出S S和和R R 。如果僅提。如果僅提取相對靜校正量

47、，只要求出取相對靜校正量，只要求出即可。即可。 （2）同一點放炮，多道接收，各道炮點靜校正量相同，接收點的靜校正）同一點放炮，多道接收，各道炮點靜校正量相同，接收點的靜校正量不同。量不同。 （3 3） 同） 同一接收點接收多個炮點的記錄， 這些記錄的炮點靜校正量不相同一接收點接收多個炮點的記錄， 這些記錄的炮點靜校正量不相同 ，但接收點靜校正量相同。但接收點靜校正量相同。 將同一炮各道的靜校正量相加平均后得到某炮點S的絕對靜校正量。將屬于同一接收點各道的靜校正量相加平均后得到該地面點R的絕對靜校正量。設(shè)設(shè)某某道道總總校校正正量量用用j j（j=1, 2, 3, 24）表表示示，各各激激發(fā)發(fā)點點

48、和和接接收收點點的的靜靜校校正正量量分分別別用用S S j j （j=1, 2, 3, 24） 和和R R j j （j=1, 2, 3, 24）表表示示，對對于于屬屬于于同同一一炮炮點點 S 的的各各道道記記錄錄則則有有 第第 S 炮炮第第一一道道總總校校正正量量 1 1 = =S S+ +R R 1 1 第第 S 炮炮第第二二道道總總校校正正量量 2 2 = =S S+ +R R 2 2 第第 S 炮炮第第二二十十四四道道總總校校正正量量 2 24 4 = =S S+ +R R 2 24 4 以上各道是在同一炮點放炮，炮點靜校正量相同，都是S。但接收點位置不同，所以接收點的靜校正量分別是R

49、 1 ,R 2 ,R 3 ,R 24 。 將將以以上上 24 道道的的總總校校正正量量相相加加再再平平均均，即即 241241241241241241241241jRjSjRjjSjj （4.4.2） 241241241241jRjSjj （4.4.2） 令令 241241jj （4.4.3） 在低速帶變化劇烈地區(qū)，241241jRj是由一些或正或負(fù)的項相加平均的結(jié)果，根據(jù)上面假設(shè)（根據(jù)上面假設(shè)（2 2） ，241241jRj應(yīng)趨于零 241241jRj0 （4.4.5） 所以最后得所以最后得S （4.4.6） （4.4.6）式表示，將屬于同一炮的各道的時差（或這些道的靜校正量）相加平）式表示

50、，將屬于同一炮的各道的時差（或這些道的靜校正量）相加平均后得到了炮點的靜校正量。這樣就求出了某炮點均后得到了炮點的靜校正量。這樣就求出了某炮點 S 的絕對靜校正量的絕對靜校正量S。 例如， 對例如， 對 6 次覆蓋，次覆蓋， 24 道記錄， 同一地面位置道記錄， 同一地面位置 12 次重復(fù)接收 （圖次重復(fù)接收 （圖 4 4- -4848） 。） 。第一炮第一道靜校正量第一炮第一道靜校正量 1 1 = =S 1S 1+ +R 1R 1 第二炮第三道靜校正量第二炮第三道靜校正量 2 2 = =S 2S 2+ +R 3R 3 第十二炮第二十二道靜校正量第十二炮第二十二道靜校正量 12 12 = =S

51、 12S 12+ +R 23R 23 根據(jù)上面同樣的道理，將根據(jù)上面同樣的道理，將 12 道靜校正量相加，道靜校正量相加， AjSjjjRjSjjj 121121)12(121121121121121121 （4.4.7） AjSjjjRjSjjj 121121)12(121121121121121121 （4.4.7） 令令 121121jj （4.4.8） 根據(jù)假設(shè)（根據(jù)假設(shè)（2）有：）有： 121121jSj 0 （4.4.10） 所以，最后有所以，最后有 A （4.4.11） （4.4.11）式表示屬于同一地面位置點接收的各記錄道的時差（或靜校正）式表示屬于同一地面位置點接收的各記錄道

52、的時差（或靜校正量） 相加量） 相加平均后得到該地面位置點的靜校正量， 這樣就求出了某接收點平均后得到該地面位置點的靜校正量， 這樣就求出了某接收點R的絕對靜校正量。的絕對靜校正量。 1）共深度點道集內(nèi)求取相對靜校正量）共深度點道集內(nèi)求取相對靜校正量 這個方法是在共深度點道集內(nèi)相對于參考道（或叫主道）用互相關(guān)法求取相對這個方法是在共深度點道集內(nèi)相對于參考道（或叫主道）用互相關(guān)法求取相對靜校正量， 并將各道相對參考道進(jìn)行靜校正， 以改善共深度點道集的疊加效果。靜校正量， 并將各道相對參考道進(jìn)行靜校正， 以改善共深度點道集的疊加效果。 具體步驟：具體步驟： （1 1） 形成參考道形成參考道目的是把

53、共深度點道集的各記錄道對齊目的是把共深度點道集的各記錄道對齊 設(shè)設(shè) xj（t）表示每個地震道的振幅，）表示每個地震道的振幅，j 表示共深度點道集內(nèi)各道的序號，疊表示共深度點道集內(nèi)各道的序號，疊加結(jié)果用加結(jié)果用 y（t）表示，則有）表示，則有 )(1)(1txNtyNjj （4.4.12） 此處此處 N 是覆蓋次數(shù)是覆蓋次數(shù),設(shè)設(shè) k 是振幅抽樣序號，是振幅抽樣序號，K 是振幅抽樣個數(shù)。是振幅抽樣個數(shù)。 )(1)(1kxNkyNjj k=1, 2, , K （4.4.13） y（k）就是參考道。）就是參考道。 說明： 這里說明： 這里 N 道地震記錄涉及到道地震記錄涉及到 2N 個地面位置點，

54、其靜校正量有正有負(fù)，個地面位置點， 其靜校正量有正有負(fù)，滿足假設(shè)條件（滿足假設(shè)條件（2） ，疊加后初步消除了一部分靜校正量，可作為參考道。） ，疊加后初步消除了一部分靜校正量，可作為參考道。 （2）互相關(guān)方法求取共深度點道集內(nèi)各道的相對靜校正量）互相關(guān)方法求取共深度點道集內(nèi)各道的相對靜校正量 在選定時窗內(nèi)計算共深度點道集中各道和參考道的互相關(guān)，求出靜校在選定時窗內(nèi)計算共深度點道集中各道和參考道的互相關(guān)，求出靜校正量。正量。 互相關(guān)公式：互相關(guān)公式： iTTiixyyxr21)(，M ，210 （4.4.14） M 為最大靜校正量；為最大靜校正量；T1、T2表示時窗起始和終了時間；表示時窗起始和

55、終了時間；T2-T1為為時窗長度。時窗長度。 用（用（4.4.14）式編制程序可求取相關(guān)函數(shù)曲線。相關(guān)函數(shù)曲線的極）式編制程序可求取相關(guān)函數(shù)曲線。相關(guān)函數(shù)曲線的極大值對應(yīng)的大值對應(yīng)的值值k便是此道的相對靜校正量。對共深度點道集內(nèi)的各道便是此道的相對靜校正量。對共深度點道集內(nèi)的各道均用上述互相關(guān)方法求取均用上述互相關(guān)方法求取k值，各道用對應(yīng)的值，各道用對應(yīng)的k值作靜校正。值作靜校正。 注意：這時求出的靜校正量是炮點和接收點相對靜校正量相加的結(jié)果，注意：這時求出的靜校正量是炮點和接收點相對靜校正量相加的結(jié)果，不是單個的炮點或接收點的靜校正量的絕對值。不是單個的炮點或接收點的靜校正量的絕對值。 （3

56、）參考道加工）參考道加工 l 用以上方法將共深度點道集內(nèi)各記錄道調(diào)齊以后，疊加效果得到提高。用以上方法將共深度點道集內(nèi)各記錄道調(diào)齊以后，疊加效果得到提高。 l 但是，由于各深度點道集的參考道不一樣，在各深度點之間同相軸對不齊。但是，由于各深度點道集的參考道不一樣，在各深度點之間同相軸對不齊。 l 用參考道加工的方法使各參考道一致起來，主要方法有組合、混波、三道用參考道加工的方法使各參考道一致起來，主要方法有組合、混波、三道相位均勻化或者扇形濾波等。相位均勻化或者扇形濾波等。 2）用互相關(guān)方法求炮點和接收點的絕對靜校正量）用互相關(guān)方法求炮點和接收點的絕對靜校正量 （1） 方法原理方法原理 炮點靜

57、校正量的求取：求出炮集中各道的相對靜校正量，再用來求取此炮點的炮點靜校正量的求?。呵蟪雠诩懈鞯赖南鄬o校正量，再用來求取此炮點的絕對靜校正量。絕對靜校正量。 以以 24 道記錄為例，用道記錄為例，用 xj示某炮示某炮 S 的各道記錄（的各道記錄（j=1, 2, , 24） ，對應(yīng)） ，對應(yīng) 24 個共個共深度點， 用深度點， 用jx表示該炮第表示該炮第 j 道所在共深度點道集各道疊加的結(jié)果。 對六次覆蓋有：道所在共深度點道集各道疊加的結(jié)果。 對六次覆蓋有： jx=61njnx （j=1, 2, 3, , 24） （4.4.15） 式中式中 n 表示覆蓋次數(shù)序號（表示覆蓋次數(shù)序號（n=1, 2

58、, , 6） 。將每個共深度點道集六道中與） 。將每個共深度點道集六道中與 S 炮無關(guān)炮無關(guān)的其余五道疊加起來，再與的其余五道疊加起來，再與 S 炮有關(guān)的道作互相關(guān)，互相關(guān)公式如下：炮有關(guān)的道作互相關(guān)，互相關(guān)公式如下： )(jjjjxxxr （j=1, 2, 3, , 24） （4.4.16） 式中，式中，rj表示相關(guān)值，表示相關(guān)值，是相關(guān)符號。是相關(guān)符號。 由此可得由此可得 24 個相關(guān)值個相關(guān)值 r1 , r2 ,r24 。令令 yj=jjxx， 式（式（4.4.16）可寫成：）可寫成： jjjyxr （j=1, 2, 3, , 24） 顯然顯然 yj=jjxx相當(dāng)于前面的參考道。相當(dāng)于前

59、面的參考道。 jjjyxr （j=1, 2, 3, , 24） 相關(guān)函數(shù)相關(guān)函數(shù) rj（）具體）具體表示為：表示為：kjTkkjjyxr,0,)( （4.4.17） 式中式中 T 表示相關(guān)時窗長度，表示相關(guān)時窗長度，k 表示時窗中抽樣序號。再將表示時窗中抽樣序號。再將 24 個相關(guān)函數(shù)相加得到個相關(guān)函數(shù)相加得到總的相關(guān)函數(shù)總的相關(guān)函數(shù) R（） ：） ： kjjTkkjjjyxrR,2410,241)()( （4.4.18） （1620， 或或 3240，） 用（用（4.4.18）式編制程序算出）式編制程序算出 R（）值以后，以）值以后，以為橫坐標(biāo)，為橫坐標(biāo)，R（）為縱坐）為縱坐標(biāo)作圖，得到一條

60、標(biāo)作圖，得到一條 R（）的相關(guān)函數(shù)曲線（圖）的相關(guān)函數(shù)曲線（圖 4-51） 。相關(guān)函數(shù)曲線極大值對應(yīng)） 。相關(guān)函數(shù)曲線極大值對應(yīng)的的時間就是第時間就是第 S 炮的炮點靜校正量炮的炮點靜校正量S S 。 接收點靜校正量：首先求取某接收點接收點靜校正量：首先求取某接收點 R 所記錄各道的相對靜校正量，再求取此所記錄各道的相對靜校正量，再求取此接收點的絕對靜校正量。方法同上。接收點的絕對靜校正量。方法同上。 （2）方法步驟：）方法步驟： 由于靜校正量的不精確，需多次疊代，故整個步驟是一個大的循環(huán)過程：由于靜校正量的不精確，需多次疊代，故整個步驟是一個大的循環(huán)過程： 先求全測線上的各炮點靜校正量。并保