1、三維折射波靜校正（譯、校對）1.0技術RefractionMiser射靜校正量的。SicsCompuion 模塊是用于建立折射模型和計算折2.0二維與三維折射靜校正方法折射靜校正算法可以歸納為兩類：RVI 和 RVD。RVI 技術是根據特定的射線模型進行射線追蹤，并將地面上每一點的延遲時計算出來。RVI 技術計算出的延遲時是地表一致性的。例如，互換算法（Hawkins,1961）將射線路徑歸納為圖（1）AGB圖 1 互換算法的射線路徑計算方程為方程（1）：TG=（TAG+TGB-TAB）/2（1）TG 為 G 點的延遲時，TGB、TAG、TAB 分別為縱波沿路程 GB、AG 和 AB 的旅行時

2、。擴展廣義互換法（EGRM）就是一種 RVI 技術，在 EGRM 模塊中有詳細說明。RVI 技術在復雜多變的地區具有良好的空變特性，但是 RVI 技術的射線路徑模型只局限于二維，如圖（1）所示站點A、G、B 必須在一條直線上。RVD 技術是利用最小二乘法，通過迭代計算出 點、檢波點的延遲時和折射層的速度。計算過程中，由于折射層速度的確性可能影響到長的精度。但是由于 RVD 技術能夠將點與檢波點的延遲時分離開，所以它計算出的短波長靜校正精度比較高，更為重要的一點是 RVD 技術不僅僅局限于二維數據。REFRACTION-MISER 模塊中的最小二乘法就是 RVD 算法，其實本模塊就是嚴格按照三維

3、設計的。比較EGRM 與REFR ACTION-MISER（用加載三維網格的二維測線作試驗）， REFR ACTION-MISER 計算出的結果要比EGRM 效果好。1.2折射靜校正計算原理對于只存在一個風化層的情況下，見圖 2，最基本的折射原理可用公式（2）表示：圖 2單個風化層的折射模型TAB=AX/V1+XY/V2+YB/V1（2）AX 為YB 為XY 為波沿AX波沿 YB的路徑的路徑波沿折射界面的路徑V1 風化層速度V2 折射層速度都知道波由一種介質到另一種介質，如果兩種介質的速度不同，就會發生折射，當折射角達到 90 度時，波就稱為折射波，只有下伏地層速度大于上覆地層速度時，才能產生

4、折射波。在這種情況下入射角稱為臨界角，根據Snell 定律，臨界角等于 sin-1（v1/v2），圖 2 中臨界角為AX 與 XX之間的夾角。方程（2）轉化為方程（3）：TAB=ZACOS/V1+AB/V2+ZBCOS/V1（3）ZACOS/V1 AB/V2 ZBCOS/V1ZAA 點延遲時線性校正部分B 點延遲時A 點處風化層厚度ZB ABB 點處風化層厚度A、B 兩點距離臨界角SIN-1（V1/V2）方程（3）是將圖 2 所示的射線路徑分為三部分，如圖 3 所示，第一部分是與 A點有關，是射線以速度V1 沿AX的旅行時減去射線以速度V2 沿折射界面AX的旅行時；第二部分與 B 點有關，計算

5、方法與 A 點相同；第三部分是射線以速度V2 通過AB 兩點所花費的時間；這是一個線性的射線模型。圖 3單個風化層的折射模型 折射波射線路徑被劃分為 檢點延遲時和線性校正三個射線模型方程 3 可以進一步簡化，并擴展到 N 層：TAB= TA+AB/Vn+1+ TB（4）這里B 為多層延遲時的總和，可表示為：ZI 為第I 層的厚度VI為第I 層的速度 為第I 層的臨界角REFRACTION-MISER 利用-方法分離出一系列的點、檢波點延遲時和折射層速度，用最小二乘法的觀點來看，這都與拾取的初至時間相關。1.3初始值的給定程序自動分析每個 CMP 中各道的拾取時間和偏移距的函數，并沿偏移距的方向

6、量取拾取時間的速度，將這個速度作為折射速度的初始值，另外在計算過程中，每一點延遲時的初始值等于這一點的旅行時減去 linear moveout 部分，再除以 2。1.4最小二乘法延遲時的計算REFRACTION-MISER 最大的優點就是給定延遲時、linear moveout 的初始值，就能夠迭代出最終計算結果。在與折射靜校正有關的文獻中，解決這個問題有很多方法：Hson 和Rul（1984）的基于模型的廣義線性反演法，Farrell分解法等。和Euwema（1984）方程（4）所體現的是拾取時間最一般的表達方式。對于一個三維工區拾取時間點多達幾百萬個，用最小二乘法分析，將這些拾取時間分為三

7、部分，對于每一個拾取時間，這三部分必須被計算出來。1、 source delay time點延遲時2、 detector delay time檢波點延遲時3、 cell-consistent refraction velocities折射層速度在迭代過程中，分別計算出這三部分，換句話說，可以先計算 source delaytime，再計算 refraction velocities，最后計算 detector delay time，當然計算次序是可選擇的，由用戶控制（1.5 Iterations）1.4.1站號標識方程 4 還有很多種表達方式，它的解不是唯一的，克服方程的唯一性的一個簡單的方法

8、就是將那些位置相近的點和檢波點進行耦合，這樣耦合在一起的點和檢波點就具有相同的延遲時。通常標識一個單都是通過道頭字SION-NUM-SOURCE 和IDENT-NUM，檢波點用 SION-NUM-DATECT，耦合的目的就是將點與其周圍位置相近的檢波點分為一組，并賦給它們一個相同的組站號。本模塊是將上面講到的一組點、檢波點作為一個耦合站點，一個組站點的延遲時就是指耦合后的檢點組的延遲時。先計算耦合組站點的延遲時，再計算單個點、檢波點的延遲時，這樣有助于提高計算結果的精度，使其更加收斂、準確、可靠。另外，耦合之后，應當確定 refraction velocity，一旦確定就可以計算出與這個速度相

9、一致的點、檢波點靜校正量。本模塊提供了四種耦合類型：source-detector coupling，sourource coupling，detector-detector coupling 和delay-time smoothing。定義前三種耦合類型有兩種方式：manual 和 automatic（人工和自動）。最重要的耦合類型是 source-detector coupling，source-detector coupling 是在SOURCE-DETECTOR 迭代選項中（Orderof Estimation Terms parameter， iteration set）將一個一個的

10、檢點組定義成一個一個的組站點以便計算出這些組站點的延遲時。Sourource coupling 是計算那些位置相近耦合到一起的 點的延遲時，Detector-Detector coupling 是計算那些位置相近耦合到一起的檢波點的延遲時。Delay-time smoothing 在 1。4。 3Automatic Coupling Specification 中詳細說明。有些情況下 、檢站點之間本身就是耦合的，例如對于二維測線來講，當 點與檢波點位置非常接近時，就會出現點與檢波點的站號相同的情況，這就是通常所說的點在檢波點放，這種情況下就沒有必要再去定義網格或耦合站點了，因為數據本身點和檢波

11、點就是耦合的，這種情況下建議采用自動耦合方式：automatic。1.4.2Manual Coupling Specification采 用 人工 偶合 方式 Manual ， 需 要 在本 模 塊之 前再 加一 個模 塊GRID-DEFINE，用來定義一個地面網格，這樣就可以計算出點、檢波點耦合后的組站號。對于一個有代表性的三維工區來講，定義的地面網格應當確保在線與檢波線交叉處的點和檢波點落入同一個面元內。如果網格定義的合適的話，那么在所有線與檢波線交叉處的面元內所包含的點、檢波點的個數和組合方式都應當是相同的。圖 4、圖 5 為兩個例子。這樣落入同一個面元中的所有點和檢波點就被耦一組，并被

12、賦予相同的組站號，通過上圖可以看到，有些面元中只有兩個點，有些面元中只有一個檢波點，而有些面元中同時有兩和一個檢波點，在檢線交叉處的面元中，檢點個數與組合方式完全相同。圖 4 中兩個點與一個檢波點被耦合到一起，且具有相同的組站號，這種思路是通過以一種系統的方式將整個工區內的點和檢波點到一起，這樣就可以提供一種點、檢波點統計的耦合方式。5-6定義的地面網格的面元大小是可變的，一個大的面元可以同時個點和 5-6 個檢波點，最佳的面元尺寸大小為能夠將檢線交叉處的檢點劃入同一個面元內最小面元的大小。在檢線交叉處以外的地方， 點、檢波點分別被分配到不同的面元中去。如圖 4、圖 5 所示。人工耦合（man

13、ual）時，組站號必須放在用戶定義的道頭字中，這個道頭字由用戶在選項Source and Detector Definition parameter set 的Secondary Identifier Literal for Source and Secondary Identifier Literal for Detector 中填寫。GE 庫中必須定義這些道頭字，也就是說 GE 庫中的點表格中必須有 Secondary IdentifierLiteral for Source 中所填寫的道頭字，檢波點表格中必須有 Secondary IdentifierLiteral for Detect

14、or 中所填寫的道頭字。另外拾取時間表格中也必須有這些道頭字（也就是midpo表格）。當選用人工耦合（manual）選項時，SOURCE-DETECT，COUPLED-SOURCE，COUPLED-DETECT 迭代過程中所用的道頭字就是上述的組站號。通常， 選用人工耦合（ manual） 選項時， 流程中還應當包含兩個模塊TRACE-HEADER-MODIFY 和 GRID-DEFINE。將 SION-NUM-SOURCE 和SSION-NUM-DETECT 中的值分別保存到SION-NUM-SOURCE-FIELD 和ION-NUM-DETECT-FIELD 中去，然后打開 GRID-DE

15、FINE 模塊的PrestackRenumbering 選項中Renumber Sion Nnmber at Source and Detector Location 參數中的 UPDATE-PRESERVE 選鍵，然后把field sion 用做secondary literal，這時field sion 中放的就是組站號。1.4.3自動耦合定義自動耦合時，本模塊是以點和檢波點作為初至拾取時間的基本的。也就是說，具有相同道頭字組合 SION-NUM-SOURCE 和 IDENT-NUM，或者具有 SION-NUM-SOURCE 或 IDENT-NUM 和另外一個用戶定義的第二個道頭字組合的所

16、有拾取時間是屬于同一個單的。同樣，具有相同道頭字SION-NUM-DETECT，或具有 SION-NUM-DETECT 和另外一個用戶定義的第二個道頭字的所有拾取時間是屬于同一個檢波點的。用戶是否定義第二個道頭字取決于是否填寫參數：Secondary identifier literal for sour和Secondary identifier literal for detectors（source and detector definition 選項）。在檢點自動耦合時，REFRACTON-MISER 模塊把那些認為位置相近的檢點的初至時間歸為一組并賦給它們一個相同的組站號。如果計算過程

17、中有一或一個檢波點被認為遠離其他或檢波點，那么它也被認為是一組并賦予一個組站號。這些組站號的用途就是確定哪些拾取時間可以被放到一塊分析，以及得到的延遲時便于檢點耦合。盡管這些組站號可以輸出到庫文件中去，它的作用只是自動耦合過程的一個質量。用戶可以查看這個信息以確定耦合是否充分。在分組過程中，REFRACTION-MISER 是按照檢點列表的順序開始檢索的，首先順著列表檢查點和檢波點的X、Y 坐標，找到一，這一要有一個檢波點在定義檢波點半徑范圍內，如果找到這么一，那么這一和它的所有在檢波點半徑范圍內的檢波點就成為第一個檢點組，其他檢波點都將不列入這一組，接下來，檢索這一組內的其他點，能劃入這一組

18、的點要同時滿足兩個條件：一、這些要在選定的那一的點半徑范圍之內，二、這一至少要有一個檢波點在檢波點范圍之內。接下來給這一組的所有點和檢波點分配一個相同的組站號。這些點和檢波點隨之從列表中刪除，接下來在檢表中重復以上過程，一直到再也找不到一滿足以上條件。剩下的點、檢波點就各自被賦予一個單獨的檢組站號。檢波點之間的耦合與檢點之間的耦合過程相似但是是獨立進行的。下面是 4 個自動耦合的例子（iteration 選項中的order of estimation terms數）：參例 1：rs=0，rd=30見圖 6：對 SOURCE-DETECT 選項， 點和其左右兩邊最近的兩個檢波點耦合，并分配到一個

19、相同的組站號。如果定義 COUPLED-SOURCE 選項，作業就會強行中斷，這是因為用戶要求進行 點耦合，但是由于 點耦合半徑為 0，所以就得不到耦合站號。在這種情況下 rs=110 就比較合適。對其他選項來講，站號保持原來的不變。例 2：rs=105，rd=60見圖 7：對 SOURCE-DETECT 選項，點 2、29 檢波點 1、3、17 和 18 耦合到一起并賦予同一個站號 2。同樣點 8、25 檢波點 7、9、21、和 22 耦合到一起并賦予同一個組站號 8。對于 COUPLED-SOURCE 選項，點 2、5 和 29 耦合到一起得到同一個組站號 2，例 3：rs=60,rd=3

20、0點 8、11 和 35 耦合到一起得到同一個組站號 8。見圖 8：對SOURCE-DETECT 選項， 點 4、57 和檢波點 3、5 耦合到一起得到相同的組站號 4。對COUPLED-SOURCE 選項，同的組站號 4。點 4、57 和 18 耦合到一起得到相例 4：rs=30,rd=30見圖 9：對 SOURCE-DETECT 選項， 點 2、3 和檢波點 1、4 耦合到一起得到相同的組站號 2，同樣 點 6、7 和檢波點 5、8 耦合到一起得到相同的組站號 6。對 COUPLED-SOURCE 選項來講， 點 2、3 耦合到一起得到相同的組站號 2，點 6、7 耦合到一起得到相同的組站

21、號 6。1.4.4站耦合與延遲時平滑（略）1.4.5refraction velocity estimation and smoothing通過參數 spatial smoothing，用戶可以定義一個矩形平滑因子來平滑折射層的速度場。這個因子是按二維方式定義的，需要填寫兩個參數：INLINE 方向面元個數和 CROSSLINE 方向面元個數。例如，如果用戶定義 INLINE 方向為 7， CROSSLINE 方向為 11，那么就會產生如圖（10）所示的矩形平滑因子。用于平滑的面元和面元尺寸是在grid 參數中定義的，如果沒有定義此參數，則缺省為觀測系統中的面元大小。折射速度估算包括分割折射層

22、的射線路徑。分割尺寸由 general 選項中的partition sizeforrefractionvelocity 參數控制。最大的和缺省的分割尺寸等于最大偏移距。如果應用缺省值，就不作任何射線分割，同時速度就由檢中心點確定，可以通過速度平滑參數修改。當選擇不作分割，建議平滑參數為 1 倍到 1。5 倍的最大偏移距。當Partition size 給定一個小于最大偏移距的值時，就可以得到較為準確的折射層速度，當然比較費機時。花費大量的機時去作射線分割到底值不值得取決于假設不作分割時折射速度與延遲時估算相比，變化是否大。用模型數據實驗結果表明，partitioning 在一定程度上提高了延遲

23、時的估算質量。所以在實際計算過程中，partitiong 肯定是有好處的。建議 partition size 為檢波線距的一半， smoothing diameter 為 1.5 倍的 partition size 或 3 倍的網格面元大小，這就足夠了。1.4.6refraction velocity constrageneral parameter set 中的refraction velocity constra參數是計算程序中的一個約束參數。這個參數將延遲時與折射速度的非單一解強行計算出單一的延遲時。因為實際折射速度與程序計算時所給定的初始速度差異很大。這樣速度與延遲時分離的結果肯定不好

24、。尤其在拾取的地方（如三維工區的邊緣）延遲時和速度就更為確了。這時如果沿折射層分離射線路徑，就會提高工區邊緣計算結果的精度。當此參數設為 0。0 時，就沒有約束條件。在工區邊緣計算出的延遲時和速度就與初始值相差甚遠，如果這個值大于 0。0，折射速度就被約束到一個與初始值相近的值上去。15 較輕，但是是非常必要的，約束值越大，就越有利于折射速度向初始值靠近，大一點的約束值可以將那些與折射速度不符的拾取時間放到其它項中去（如geometry errors）。在應用 constra應當特別，當工區邊緣折射速度確實是變化十分劇烈（鹽丘存在），即使用一個非常小的約束值都會導致錯誤的結果。在這樣的工區據的

25、經驗，不用constra都比用一個即使非常小的constra值要好，另外這種情況下，partition 也不要用。因此，當工區邊緣折射速度沒有突變時，如果選用了partition 那么建議用一個小的 constra值。如果不用 partition 或地層有突變時，建議不用 constra。另外不建議使用非常大的 constra值，因為一個特別大的值會導致折射速度過度的偏向初始值。不過在拾取時間點相對較少的地方，大一點的 constra于 5）有利于得到較準確的結果。1.5iterations延遲時和折射速度的計算求解過程就叫作迭代。迭代次數由用戶定義，每次值（大迭代的計算結果的準確性都會提高，

26、直到結果收斂。Iteration parameter set 中的order of estimation terms 參數，定義迭代次序，可選的延遲時計算類型有點的 SOURCE、檢波點的 DETECT、 與耦合的COUPLED-SOURCE 、檢與檢耦合的 COUPLED-DETECT 、檢點耦合的SOURCE-DETECT。速度計算類型有 REFRACT-VEL。在迭代計算過程中，每一個計算類型都有應有的次序。下面所列舉的迭代次序是所建議的基本次序。在這個例子中，首先是 REFRACT-VEL 和 SOURCE-DETECT 迭代 15 次。這里在迭代REFRACT-VEL 過程中，建議平

27、滑直徑smoothing diameter 取值為所分析到的最深的折射層的最大偏移距的一半。迭代過程中點和檢波點得到充分點 COUPLED-SOURCE 和檢波點的耦合是十分關鍵的。流程接下來是COUPLED-DETECT 的兩次迭代。經過迭代，就依據拾取時間和現有的折射層以最小二乘法將延遲時和折射層速度計算出來了。1.5.1geometry error estimation本模塊還可以計算出點和檢波點在位置上的水平誤差。基本是初至時間的動校時，偏移距錯誤與否對動校量影響非常大。因此通過動校可以檢查出觀測系統中的錯誤。例如，取一條中間放的二維測線，這條測線是沿 X軸方向施工的，則觀測系統錯誤的

28、旅行時可以表示為：如果檢點坐標在inline 方向上有錯誤，那么旅行時的變化取決于誤差量與折射速度的商。如果檢點在另外方向上發生偏離，就會導致延遲時的計算誤差。因此方程很容易推廣到三維的形式。通常對于一或一個檢波點來講，在一個折射層所在的偏移距范圍內有很多道。任何一個點坐標的誤差都會影響到它的所有檢波點，影響大小取決于檢點的誤差方向、誤差大小和方向本身所帶來的誤差。檢波點受點影響與點受檢波點影響是一樣的，因此如果檢波點坐標出現錯誤，它同樣影響到與它有關的歸入 geometry error點。因此在折射方程與最小二乘法的算法中將這種terms 選項中。能否成功的估算出觀測系統誤差主要取決于 GE

29、 庫中的信息。假設有一的所有檢波點都向一個方向偏離理論位置， 那么情況會是什么樣子呢？一個INLINE 方向的觀測系統錯誤導致某一所有的偏移距都有相同的變化，因此動校時對每一道的初至時間的影響是等效的。如果這種情況發生，觀測系統的錯誤就很難被發現。在觀測系統誤差估算過程中，對三維工區來講，INLINE 方向的錯誤容易被發現，CROSSLINE 方向的錯誤就不容易被發現。因為 INLINE 方向的錯誤對初至時間影響較CROSSLINE 方向的錯誤大的多，因此對 INLINE 方向的誤差估算效果好，而CROSSLINE 方向效果就一般。工區邊緣，基本都是單邊放，這些地區估算效果不是很好。本模塊只是

30、將那些誤差較大的觀測系統錯誤識別出來。在iteration parameter set 選項的order of estimation terms 參數中有兩個迭代選項可以估算觀測系統錯誤SOURCE-GEOMETRY、DETECT-GEOMETRY。觀測系統估算有一個缺陷就是它的自動行為。對于那些符合計算結果的拾取時間的選取的權系數就要高，對于那些不符合計算結果的拾取時間的系數就很低，即使這些拾取時間是由于觀測系統錯誤而導致的不符合計算結果。解決這個問題的辦法就是建立一個參數表，這個參數表在第一次觀測系統誤差分析的迭代過程中就參與計算。不過，一般傾向于在觀測系誤差分析之前就消除那些離散點，所以

31、作了下面一個折中的流程, 這是一個比較典型的迭代流程，其中包括觀測系統誤差估算，具體參數看下面流程：例 6:*ITER在觀測系統估算過程中，要給定一個門檻值，只有那些超過門檻值的拾取時間才被認為是錯誤的。例如，在拾取時間中加入一些噪音，那么超出 10 英寸或者稍微小一些的值就都不可靠了。門檻值在 geometry parameter set 中的 sourcegeometry error tolerance 和detecteometry error tolerance 中填寫。另外每一步估算所計算出的錯誤的拾取時間的個數也是有限制的。由geometry parameter set 中的um p

32、ercentage ofdisplaced sour和um percentage ofdisplaced detectors 控制。估算觀測系統誤差理想的情況就是點和檢波點在一條直線上。這種排列對于計算INLINE 方向誤差較為理想，但是對于 CROSSLINE 方向誤差效果就差的多了。因此系統提供了一個平衡因子參數來提高 檢點在 INLINE、CROSSLINE方向上誤差的估算的質量。由 geometry parameter set 中的minimum in-line balance factor 和minimum crossline balance factor。最大值是 1。0，最小值是

33、 0。0，0。0的意思是估算過程中不考慮排列之間的關系。對于大多數數據來講 0。1 或 0。2就足夠了。迭代結束后，觀測系統誤差就輸出到以下道頭字：DELTA-X。SOURCE DELTA-X。DETECT DELTA-Y。SOURCE DELTA-Y。DETECT這里不建議用這個估算方法去解決觀測系統問題。就象其他統計方法一樣，這些選項只是估算出有問題的拾取時間，但是估算結果受多種的影響里，例如噪音、拾取錯誤、錯誤的偏移距等等，不僅僅只受觀測系統錯誤的影響。因此這種估算只是一種用于分析段。1.5.2weighing and rejecting為了確保延遲時計算結果的可靠性，應當將那些錯誤的拾

34、取時間刪除掉。這里提供了三種方法：第一種類型是，拾取時間質量因子，每個拾取時間有一個值。是用來權衡各個拾取時間點對計算結果的貢獻量大小的。在輸入的 midpotable 中的道頭字 FB-PICK-QUAL-FACTOR 中放著。如果這個值小于等于 0 或大于 1，這個拾取時間將不參與計算。如果表格中沒有這個道頭字，所有的拾取時間自動設為1。第二種類型是在每次迭始的時候根據拾取時間與上次迭代結果所偏離程度對這個拾取時間確定一個系數。這個系數是在 weighting 選項中填寫的。在這個選項中，需要填寫 iteration number，slope of weighting curve 和cut

35、off deviation 參數。這些參數通過內插、外推，就可以為每次迭代都提供一個權系數。 slope of weighting curve 定義斜坡值。cutoff deviation 定義權值為 0。5 時的偏移量。這里建議用戶在所有迭代過程中slope of weighting curve 都置為 5，然后根據實際情況來變化cutoff deviation 的大小。權函數在 0。001 到 0。999 之間變化。不同偏移量的權值系數都會輸出到一個表格中。最后每個拾取時間得到一個權值系數，這個系數等于上面所計算出的權值系數與庫文件中的 FB-PICK-QUAL-FACTOR 的乘積。一般

36、來講 cutoff deviation 的初始值應當稍微大一點（如 100ms）。隨著迭代次數的增加，偏移量是逐漸減小的，減小到定義的最小值時，后面的迭代過程中就保持這個常量。（見例 7）。例 7：*WEIGHT這個例子中，在第 7 次迭代時就開始計算權值，偏移量為 100ms。在第 20次迭代時，偏移量就減少為 20ms。在第 7 次與第 20 次迭代之間偏移量是在 100ms與 20ms 之間內插的。當迭代次數超過 20 次，權值都與第 20 次的值保持相同。每個點、檢波點最后一次迭代計算出的權系數都要在寬行中打印出來。weighing parameter set 選項的 weight t

37、olerance 參數中需要填寫一個最低權系數門檻值，權系數低于這個門檻值的的那部分拾取時間就被打印出來。每個折射層依次在寬行中打印出來，但是只有最深的那個折射層的信息才能輸出到庫文件中去。第三種類型，是設計一個窗，拾取時間值作為 Y 坐標，偏移距作為 X坐標，這個窗的中心就是通過時間-偏移距的一條直線。在general 選項填寫pick acceptance corridor zero-offset reference time 和pick acceptance corridor velocity 參數，分別表示定義窗口的零偏移距時的參考時間和參考速度。pick acceptancecorr

38、idor tolerance 為時窗的寬度。在general 選項中還有一個參數percentage of picks used inysis，填寫這個參數就表示可以任意的排除一定比例的拾取時間。一般來講，如果認為只用一部分拾取時間和用全部的計算出的效果一樣好，就可以填寫這個參數，另外可以節省大量的機時。1.5.3waterlayer如果地表有一個明顯的水界面（不是指小溪或），或者在聯片處理時合并一個沒有水層和一個有水層的工區，這時就碰到水層引起的靜校正問題。在general 選項中將參數 water layer 設為 YESapplication 選項中的 uphole applicatio

39、n type可以解決這個問題。將 uphole定義為 RAYPATH ，并且定義below-refractor uphole velocity 參數。水層被認為是一個特殊的常速層。如果定義了水層，折射層的層數應當加 1（例沒有必要定義水層的偏移距分布范圍。層下面的泥層就是第二層）。本模塊例 8：本模塊一般是從觀測系統中的以下道頭字中抽取水層信息的：DEPTH-SOURCEUPHOLE-TIME-SOURCEWATER-DEPTH-SOURCE DEPTH-DETECT UPHOLE-TIME-DETECT WATER-DEPTH-DETECT當道頭字水深 WATER-DEPTH-SOURCE、

40、WATER-DEPTH-DETECT 為 0，本 模 塊 就 認 為 此 處 沒 有 水 層 ， UPHOLE-TIME-SOURCE 、UPHOLE-TIME-DETECT、DEPTH-SOURCE、DEPTH-DETECT 等道頭字含義正常。如果參數 water layer 設為 YES，但是觀測系統庫文件中沒有水深的道頭字，本模塊會自動在庫文件中這個道頭字并設為 0。這有利于無水區向有水區過渡和工區合并。（1。9water uphole corrections 中有更詳細的介紹）。1.6Diving WavesREFRACTIONN-MISER 提供了兩種方法識別回折波。這兩種方法可以單

41、獨用也可以組合起來用，它們都是基于延遲時的外推公式，這個公式表達的是半空間上覆風化層的速度是連續的即存在一個線性的速度梯度，波的旅行時。速度梯度的存在使得波沒有沿半空間頂部，而是向下到了半空間中去。在這種情況下，旅行時表達式為方程 7：DSDd SDVr RDA點延遲時檢波點延遲時檢距折射速度參考距離回折波系數Zanzi 和 Carlini1991 年，a*（SD-RD）3 這一項表示的是當折射層是連續介質時，初至時間的上凹的曲率，RD 是接收到回折波的最小偏移距，相當于傳統靜校正分析過程中的臨界偏移距。RD 的最大值都要小于或等于用于靜校正分析的最小偏移距。一般情況下，在一個工區這個值是個常

42、數。一些建議這個值應當為 0，但是這樣計算出的延遲時有可能是一個負值。因此，有必要通過某種方式對這個臨界距離進行估算。方程 8 就給出了一個估算公式：Z風化層厚度臨界角在識別回折波的第法中，用戶要給出系數 a 的值。這個系數可以人為給出，也可以用初至拾取時間的回歸算法將其計算出來。回歸算法是在統計分析包中進行的。將一個有代表性的拾取時間樣點輸入到統計分析包中，然后擬和方程 9 所示的曲線。A、B、C 為一般項。接下來，在 general 選項中，給定 c 值refraction velocitygradient factor，給定 RD 值refraction velocity grdient

43、 mininum offset。本模塊在將拾取時間分解為延遲時與折射速度之前就將回折波去除掉。去掉速度梯度影響后的拾取時間就象是波以一個恒速通過一個簡單的折射層一樣。系數 C 也可以通過下列步驟計算出來。首先，將單按偏移距分選，顯示單，這樣偏移距是由小到大排列的，最準確的方法是將數據作一個線性動校，接下來，在回折波開始出現的偏移距處拉一條直線與初至相切，這個偏移距就是參考偏移距 RD。如果數據體真的出現了回折波，初至速度是隨著偏移距的增加而逐漸增大的，另外，隨偏移距的增大初至時間較切線時間早。然后量出直線的時差，就可以用公式 10 計算回折波的系數 a 了：例 10，假設在偏移距 10000f

44、t 處測得時差為 50ms，參考偏移距為 3000ft，那么，回折波系數為：這個例子中，general 選項中，參數refractor velocity gradient mininum offset參數應當填寫 3000，refractor velocity gradient factor 參數應當填寫 0。1458。第二種處理回折波的方法是用最小二乘法求解方程 7，這里要給出 RD 的值，具體參數在例 11 中給出：例 11這個迭代序列中，延遲時與折射速度的初始化的同時，加入了回折波系數的初始化過程。方程 10 中的時差受兩個控制，一、線性速度，二、偏移距差的三次冪，這兩個都受人為影響，因

45、此計算出的時差結果有可能不是很準確，尤其是當折射層只有在遠偏移距才能接收到的情況下。線性時差的誤差有可能導致延遲時的偏差，因此如果選用 diving-wave 迭代選項，用戶在定義偏移距范圍時應當和實際的一樣大。另外本模塊，不允許作多層折射層的回折波迭代，所以用戶可以作一層風化層的回折波迭代（或海洋數據的第一個泥層）。另外 diving-wave項還應當作一個大的平滑，在 spatial smoothing 選項中的in-line diving-wave term smoothing diameter 和crossline diving-wave term smoothing diameter

46、 參數中定義平滑直徑。建議平滑直徑應但取值為三倍的最深折射層所對應的偏移距。如果平滑直徑設為 0，平滑因子就會對整個工區都產生很大的影響，并且系數 a 為一個常數。在實際生產過程中用戶都希望能把回折波處理用于多層折射層的靜校正計算過程中去，只要把多個折射層當作續的折射層，并簡化偏移距范圍的定義和提高延遲時計算次數就可以了。1.7建立折射模型一旦拾取時間被分解為延遲時和折射速度，就必須建立一個折射層的高程模型，以便后續的靜校正量的計算。三維折射層高程的計算和二維的是一樣的（EGRM 中有詳細的介紹）。建立折射層模型要按以下步驟進行：1、將用戶定義的風化層速度作全工區內插。2、如果要用井口信息修改

47、風化層速度，就可以將井口信息轉化為風化層速度（見 1。10），然后將這個速度進行全工區內插并將用戶定義的風化層速度替換掉。3、如果定義高程替代參數（elev override），就會用新的風化層速度替代先前的風化層速度，這個新的速度是由替代高程、延遲時和折射速度三項計算得來的。反之如果不用替代高程，則風化層的高程就由延遲時、折射速度和用戶提供的風化層速度計算得出。4、如果定義了折射層高程平滑參數，那么平滑的是最深的那個折射層的高程，并且風化層的速度也相應的修改了，這樣最終的深度-速度模型就與延遲時和折射速度相一致了。本模塊提供兩種建模方法。第一種用已知的風化層速度計算出風化層高程。第二種用已知

48、的風化層高程計算出風化層速度。這兩種方法都是為了建立風化層的速度和深度模型，使這些模型與延遲時、折射速度、以知的或假設的近地表的速度或結構相一致。第法是將延遲時直接轉化為折射層的高程，所用風化層速度有兩種：一、用戶直接定義風化層速度（velocity specification 選項中的layer velocity 參數）。二、用井口時間修改風化層速度（ uphole application 選項中的 use uphole information 參數）。對于只有一個風化層的近地表模型的點的延遲時可以用公式 11 表示：給定延遲時、速度v1、v2，折射層的深度就可以表示為公式 12：折射層的高

49、程就等于點高程減去折射層厚度，見公式 13：對于多層模型，在計算第二層、第三層的高程的過程和計算第一層的相同，參照公式 5。通常是不知道風化層速度的空間高頻的，尤其是有時風化層速度只是通過有限個點粗略估算出來的。當用井口時間去計算風化層速度時，由于井是稀疏分散的所以在某種程度上來講，得到的速度是確的。由于每個點、檢波點的延遲時的精度很高，用一個確的風化層速度去計算風化層高程就會導致計算結果與實際的風化層位置不相一致。這就直接影響著靜校正的結果。為了解決這個問題，用戶可以選用折射高程替代和平滑選項。如果用戶有一些關于風化層速度的數據，那么對于定義準確的風化層速度替換（refractor elev

50、ation override）是十分有用的，計算出的風化層速度是與 refractor elevation override 相一致的，而不是與某種假設出來的風化層速度相一致的。另一方面，如果井口信息比較準確，那么對于用由井口信息計算出的風化層速度來計算折射層高程是十分有用的，接下來將這個高程作一個輕微的平滑，即相當于將風化層速度作了一個輕微的平，這個平滑后的速度是與延遲時的變化相一致的。平滑的程度直接關系到風化層速度和高程的度。一般來講，近地表反演出不同的模型就得出不同的靜校正量（這些靜校正量之間不是一個簡單的時差關系），靜校正量之間的差別取決于最底部的折射層高程的形態。在refractor

51、 elevation 中定義折射層高程替代參數，在spatial smoothing 中定義折射層高程平滑參數。1.8uphole options為了確保點和檢波點的延遲時相一致，本模塊可以將深井校正到地表。作校正就要打開uphole application 選項。對于震源的井，提供了兩種校正方式： vertical 和 raypath 。 要求井口時間和分別放到道頭字UPHOLE-TIME-SOURCE 和 DEPTH-SOURCE 中去。如果點表格中沒有這些值或有些不正確需要修改，在選項uphole override 中可以填寫正確的井口時間。垂直井口校正（vertical 方式）是簡單的

52、將井口時間加到各自點所對應的初至拾取時間上去。這種方法比較簡單，它是假設地表的點的延遲時就等于井口時間和井 初至時間之和。見圖 11，假設在 B 點可以接受到的初至時間，這個時間就等于井口時間（從 A 到 B），加上真正的初至時間（從 A 到 C）。由圖 11 可以看出這種假設的射線路徑是由 B 到 A 再到 C，只有風化層速度從地表到折射界面的變化梯度是由小到大變化的時，這個射線路徑才近似于真實情況，另外還要保證位置到地表的距離要比到折射界面的距離近。射線路徑校正（raypath），這里提供了一種最好的井口校正方法，就是先從A 點以v2 的速度校到 D 點，再以井口速度從 D 點校到 B 點

53、，這就是射線校正。對于每個井所對應的初至時間，其射線路徑校正量見方程 14：這個方程重新創建了一個射線路徑從 B 到 D 再到 C，在大多數情況下，這種方法的校正量要比垂直校正量大。其實這種方法是將井口的延遲時加到點的延遲時上去了。方程 14 可以進一步簡化為方程 15：射線路徑校正是假設點位于折射界面以上，如果井了折射界面，這種校正就確了。折射界面的點的延遲時就會有偏差，但是，由于REFRACTION-MISER 有分離點、檢波點延遲時的功能，所以在計算點延遲時的過程中這個偏差就得到補償，換句話說，這種校正的確性可以轉化為點延遲時的計算誤差。1.9waterupholecorrections

54、如果在有水的環境中施工，如河流、淺海等，點和檢波點都放置在水下，這種情況類似于井下激發，如圖 12 所示， 檢點都在水中，從圖 12 中不難看出，水層引起部分靜校正問題，本模塊就可以解決這個問題。由于在施工中水層的信息了解的比較多，要解決這部分靜校正就要從一開始就作起。觀測系統中要有點深度 DEPTH-SOURCE 、檢波點深度 DEPTH-DETECT 、點位置的水深WATER-DEPTH-SOURCE、檢波點位置的水深 WATER-DEPTH-DETECT 四個道頭信息。縱波在水中的速度 VW 為 1500m/s，如果能夠估算一個比較準確的折射層速度V3，由S 到D 的折射波的拾取時間就可

55、以變化A 到B 之間折射波的旅行時，換句話說，就是相當于在水底激發、接收。點S 處，校正量見公式 16：轉換成點延遲時的校正量見 公式 17：計算檢波點 D 處的校正量過程與上述點相同。這個校正過程與校正井口時間類似，校正時用減去法。在 General parameter set 中定義 Water Layer 參數，就可以進行水層校正了。用戶需要在 Uphole Application parameter set 中定義 Below-Refractor Uphole Velocity 參數，這個折射速度指的是計算過程中分析到的最深的一個折射層的折射速度。這個速度必須大于水層速度，另外，在 U

56、phole Application parameter set將 Uphole Application Type 定義為RAYPATH。1.10用井口信息修正風化層速度如果 Uphole Application parameter set 中 Use Uphole Information 參數選項為YES 或RESOLVE，本模塊就會用井口時間、深度信息，在計算延遲時和折射層速度時估算新的風化層速度。這個新的風化層速度會在全區內插，這樣在沒有井口信息的點、檢波點位置也會得到對應的速度。首先，本模塊會根據剩余延遲時判斷井是否打穿了折射層，剩余延遲時見公式 18：其中Td 本 點處 點的延遲時 T

57、u 為本 點處的井口時間 V2 折射層速度Vu 本 點處井口速度如果剩余延遲時大于零就位于折射層之上，也就是井沒有打穿折射層。本模塊提供兩種方式計算風化層的速度和厚度，方法的選取取決于速度信息的可靠度。如果認為整個風化層的速度和由井口速度大致相等，那么就選用第法，選擇 RESOLVE 選鍵，井口速度就等于井口深度除以井口時間，對于這種方法來講，整個風化層包括以上和井底到折射界面之間，風化層的速度都用井口速度。選用第二種方法，就應當知道井底到折射界面之間的速度變化規律，例如微測井表明，井底的速度明顯大于井口速度，這時就應當選擇 YES 選項，這種情況下，以上的風化層用井口速度，井底到折射界面之間

58、的風化層用用戶定義的第一層層速度-layer velocity。如果剩余延遲時等于零，就表明井底剛好達到折射界面，選擇 YES 或RESOLVE 計算結果是一樣的。如果剩余延遲時小于零，就意味著井打穿了折射界面，這時就可以表示為兩部分，Z1 風化層厚度，Z2 風化層底到井底之間的深度，見公式 19：同樣井口時間也可以表示為兩部分之和，從地表到折射界面之間的時間和折射界面到井底之間的時間，見公式 20：此處延遲時公式見公式 21公式 18、19、20 包括三個未知量：Z1、Z2、V2，其中 Zu、Tu 已知，Td 和V2 可以迭代出來，因此可以直接計算出 V1。這種情況下 YES 或 RESOL

59、VE 選鍵，計算結果是一樣的。總之，本模塊可以判斷每個點位置處井是否打穿折射界面，可以計算風化層速度，如果風化層速度小于用戶在 Uphole Application parameter set 中定義的最小風化層速度Minnum Weathering Velocity，則計算出的風化層速度不被采用。1.11靜校正量的計算本模塊計算井校正的方法是從地表一個折射層、一個折射層往下剝，各折射層用各自的層速度，到最深的一個折射層時，用戶需要定義一個中間基準面，用最深折射層的速度將靜校正計算到中間基準面上。中間基準面應當非常平滑或是一個平面，并且接近于最深折射層的高程。之后就是用用戶定義的替換速度將靜校

60、正量計算到最終基準面上去。最終基準面和替換速度在 Datum parameter set 中的 Final Datum Elevation 和 Replacement Velocity 中定義，中間基準面在ermediate Datum 或ermediate Datum Bias 或Spatial Smoothing 參數中定義。1.12地表基準面和剩余折射靜校正在不同靜校正方法中，如：野外靜校正、高程靜校正、折射靜校正一直就存在很多爭論。風化層速度、厚度模型不同、基準面和替換速度不同都可以導致計算出的靜校正量不同。折射靜校正的一個目的就是希望消除長波長靜校正和近短波長靜校正，另外一個目的就是