【文章內容簡介】 此 ： 相反，到站點G5的旅行時差和偏移距差是一個負值，也就是：在圖2和圖3所給的模型中，由于站點G3到A、B等距，旅行時差和偏移距差這兩者都接近于零。在上面例子中，可以繪出偏移距差旅行時差圖，如圖4所示。圖4 在圖2和圖3中所舉的例子中G1G5這5個站點得到的速度圖對圖4中5個點進行最小平方線性擬合后，其斜率的一半就是點A和B之間的速度Vn+1的估算值。這就是EGRM中基本的折射層速度計算方法。在Hybrid EGRM/迭代算法中，另一個方法也可進行速度計算。 EGRM 延遲時計算—單風化層模型圖5 這個折射模型顯示了在EGRM法中采用的折射路徑EGRM法需要三個折射旅行時間去計算一個特定站點的延遲時，這三個旅行時來自于兩個炮點位置。在圖5中為了計算G點的延遲時，畫出了三個傳播路徑。在圖5中，傳播路徑為：從炮點A到接收點B的傳播時間 ，或者反方向傳播。從炮點A到接收點G的傳播時間 ，或者反方向傳播。從點G到接收點B的傳播時間 ，或者反方向傳播。 EGRM法計算G點延遲時的公式為： 把等式（2）代入（7）可得，化簡得：而AB=AG+GB （10）因此，化簡得 這幾乎和截距法中定義的截距時是一致的（除了因子2）。 EGRM延遲時計算—兩個或更多風化層情況圖6 兩層折射模型圖6中兩層的情況下，根據(jù)折射原理，折射傳播時間TAB可描述如下：這里 DA1=地表A點和折射面1上W點之間的射線距離 DA2=折射面1上W點和折射面2上X點之間的射線距離 DB1=地表B點和折射面1上Z點之間的射線距離 DB2=折射面1上Z點和折射面2上Y點之間的射線距離 DR=X點和Y點之間的射線距離 V1=地表和折射面1之間的平均速度 V2=折射面1和2之間的平均速度 V3=折射面2下伏地層的速度等式12可簡化為：這里 ZA1= A點下方層1的厚度 ZA2= A點下方層2的厚度 ZB1= B點下方層1的厚度 ZB2= B點下方層2的厚度 AB= 地表A和B點之間的距離 a= 穿過折射層1 的角度，定義為arcsin（v1/v3） 216。= 穿過折射層2 的由SHELL定理定義的臨界角 arcsin（v2/v3）圖7 兩個風化層的折射模型象單風化層情況一樣，也需要三個折射旅行時來計算延遲時，并且使用同樣的EGRM公式。把等式13代入等式7 ，化減后得：在規(guī)則采集的觀測系統(tǒng)下，偏移距的剩余項接近于零，等式14變成了一個簡單的表達式，實際上，對于任意多層，這個表達式可擴展為一個簡單的公式。這里： ZGj=G點下方層J的厚度Vj=層j的速度?/j= 由Snell定理所定義的arcsin（Vj/Vn+1）注意到當j是最深一層的風化層時（或j=n），臨界角216。和角度?/j是一致的。 GaussSeidel 迭代旅行時分解等式2 表示折射波傳播時間是炮點延遲時、檢波點延遲時和偏移距速度項之和，或對于n層的情況下：這里：TA=炮點A的延遲時TB=檢波點B的延遲時根據(jù)前面折射層速度分析結果，可以立即計算出偏移距—速度項，這個等式中僅剩下炮點項和檢波點項沒有解決。其實，GuassSeidel算法非常簡單，重新整理等式17，炮點項可寫成：同理，檢波點項可寫成：然而，由于TB未知，TA也不能求出（反之亦然），因此，需要得到TA和TB的初始估算值。EGRM模塊中是通過采用截距法的一個變換來獲得這些初始估算值的。對于每一個地面點，EGRM根據(jù)給定的折射分層，抽取出所有以該點為中心點的折射旅行時。例如，在120站點，一個這樣的旅行時可以是從炮點110到檢波點130的旅行時，因為站點120是這個炮、檢點對的中心點或中心站點，其它的旅行時如從炮點134到檢波點106，從炮點101到檢波點139等等。如圖8所示。圖8 這個模型顯示了以120點為中心點的折射路徑因為這些旅行時隨不同的炮、檢點偏移距而變，所以他們有不同的線形校正量。如圖9所示。圖9 以站點120為中心點的旅行時間—偏移距變化圖，黑線表示不同偏移距上的折射旅行時間從每一個這樣的旅行時中減去偏移距—速度項，就等于對他們應用了線形校正或者進行了拉平。把這個中心點上所有被拉平的時間（也就是每一個減掉了偏移距—速度項的剩余時間）進行平均，然后分成一半（因為截距時間是雙程的延遲時）。見圖10。對測線上每一點都這樣做，在這些中心站點上得到的延遲時的粗解稱做初始的結構解。 平均截距時間圖10 應用線性校正后的旅行時和偏移距對應圖。穿過實線的虛線表示120點的平均截距時間初始的結構解被分配給測線上每一個炮點和檢波點。在上面的例子中，如果在120站點有一個炮點，該炮點的初始延遲時就是截距時的一半，對于檢波點也同樣。初始的結構解表示了測線上一個粗略的延遲時估算值?，F(xiàn)在，在GaussSeidel迭代算法中要求提供的炮、檢點的初始值是已知的了。首先，根據(jù)等式18，EGRM模塊重新估算測線上每一炮的延遲時。例如，EGRM抽取出所有屬于炮點101的并且是我們所需要的折射旅行時。因為折射層速度Vn+1已經(jīng)估算出，偏移距—速度項是已知的，根據(jù)初始的結構解，檢波點的延遲時值TB也是已知的，我們就可以重新計算炮點延遲時TA，對測線上的每一炮都這樣做。所有的炮點延遲時重新計算后，用等式19，把所有的檢波點延遲時也用同樣的方式重新計算一遍。這樣一個過程叫做一次迭代。EGRM實際上在炮點域和檢波點域進行了多次迭代。用戶可以控制迭代次數(shù)。 Hybrid EGRM/迭代算法EGRM和GuassSeidel迭代算法代表了兩種不同的延遲時計算方法。每一項技術都有本身的長處和弱項。例如，迭代算法對折射層的速度Vn+1的錯誤非常敏感，這是因為在等式18和19中，偏移距/ Vn+1這一項是一個主要的分量。在偏移距/ Vn+1這一項中產(chǎn)生的所有誤差都被分配到炮點延遲時項和檢波點延遲時項里。偏移距越大，Vn+1的誤差引起的影響就越大。而EGRM采用的是合并、化減折射旅行時的方法，它排除了偏移距這一項，而僅僅在偏移距剩余項中需要Vn+1，而偏移距剩余項通常較小且常常等于零。因此，EGRM法對折射層速度Vn+1的誤差相對來說不敏感。但是，EGRM由于它合并、化減旅行時這種方法本身的特點，它只能允許每個站點有一個延遲時而不能計算出位于同一站點的不同的炮、檢點延遲時。很自然的，迭代算法能夠直接計算炮、檢點各自的延遲時。在彎線或炮點垂向有偏離時，雖然炮點和檢波點有相同的站號，但是炮點和檢波點的延遲時是不同的。EGRM另一個嚴重的問題出現(xiàn)在測線的端點。因為EGRM方法對任一點的延遲時估算它都需要三個折射波旅行時，在測線的端點或者測線中間炮點跳度比較大的時候，這種算法并不總能得到足夠的折射波旅行時信息來估算延遲時。 用非互易估算法估算延遲時。（ TimeDeph Estimation using Nonreciprocal Estimation）。迭代算法通常不存在這些問題。只要在測線端點存在一個炮點，迭代算法就能估算它的延遲時（當然，假定存在所需要的折射層的旅行時，并且取決于折射層速度的精度）。同理，也可估算出位于測線端點的檢波點的延遲時。Hybrid方法綜合了EGRM法和迭代法優(yōu)點。Hybrid方法很容易理解，它首先用EGRM方法得出折射層的速度和炮、檢點延遲時，把這些延遲時作為迭代法的初始值進行迭代計算，而不用上面所描述的用截距法求得的粗略的延遲時作為迭代的初始值。Hybrid算法遵循等式17。對于任意的拾取時間TAB，當給出一系列TA和TB的估算值（從EGRM的結果中得到）后，可以重新計算出折射層速度Vn+1，這個速度和EGRM法所估算出的速度趨勢是一致的。這個等式很簡單：根據(jù)每一個初至拾取時間TAB，延遲時估算值TA和TB可計算出一個新的速度Vn+1。把這些新的速度值列成一個表，可得到測線上一個新的折射層速度場。通常情況下，這種方法在承認EGRM算法得到的結果基礎上，對后續(xù)的迭代算法能有效的約束。Hybrid算法可由EGRM算法求出的長波長的結果，用迭代法估算炮、檢點各自的延遲時。同時，迭代算法還可以解決測線的端點問題。 用非互易估算法估算延遲時在計算一個站點的延遲時，EGRM算法需要三個折射波旅行時間，如圖5和圖7所示。但是，在測線的端點和炮點跳的比較大的區(qū)域，并不能得到這樣三個旅行時。同樣，如果在測線上部分地方由于拾取錯誤或折射傳播的非地表一致性，會使初至拾取時間表現(xiàn)出非地表一致性。對于這些站點，EGRM模塊會認為計算出的延遲時不可靠。在這種情況下，EGRM模塊必須通過其他方法估算延遲時。 Hybrid EGRM/迭代算法中介紹的方法。在EGRM模塊中還提供了另外三種非互易延遲時估算法：1. 延遲時預測法2. 折射面高程預測法3. 延遲時的覆蓋外推法當使用非互易估算法Nonreciprocal Estimation Type for Layer這個選件時，這三種方法都可選用。延遲時預測法在外推測線端點處的延遲時的時候，對與外推區(qū)域相鄰的站點的延遲時做最小平方擬合。用于最小平方擬合的站點數(shù)缺省值為25，也可在Minimum number of station to overlap 參數(shù)中指定。對不在測線端點處的那些延遲時，延遲時的估算使用的是平滑技術。平滑技術用來估算那些變化量超過一定值的個別延遲時。變化量的缺省值為16，也可在延遲時變化量門檻值（timedepth variance threshold）參數(shù)中指定。平滑算法是這樣來做的：用比較可靠的相鄰站點的延遲時（也就是，較小的變化量）來計算一個加權平均系數(shù)。相鄰這些站點上的延遲時的加權系數(shù)和他們的距離成反比，和他們的變化量成反比。在對目標站點的延遲時平滑計算時，離目標站點近的站點比離目標站點遠的站點有更大的影響，同理，有更可靠延遲時的站點也具有更大的影響。計算平滑后的延遲時值Ｚｋ的等式是：　這里Zk＋j＝站點K＋J的延遲時（或折射面高程）Vk＋j＝站點K＋J的延遲時變化量K＝需要估算延遲時的站點M和N大于5，并且還要足夠大，以便在j0或J0時，至少存在一項。另外，對那些延遲時估算值是無效值的那些站點，EGRM給延遲時和變化量都賦值為1，對這些站點，需要用平滑算法來估算延遲時，實際上，也就是一種內插算法。如果認為該地區(qū)（根據(jù)地質信息或觀察）折射面基本上隨地表而起伏，推薦使用該選項。這是因為平滑和線形外推所得到的延遲時是比較圓滑的，在建立近地表模型時，轉換成的折射層厚度（不是折射面高程）也是圓滑的，這樣，當?shù)乇砥鸱鼤r，折射面高程也就隨之起伏。折射面高程預測法這個方法和延遲時預測法一樣，只不過在估算前，延遲時先被轉換成折射面高程。然后，對這些高程用最小平方擬合的方法進行外推，或者用加權平滑法進行內插。如果測線上某一段地表變化比較大，并且認為該段的初至拾取值不太可靠，我們推薦使用該選項。如果地形起伏較大，但是認為初至時間比較可靠，可以采用下面所介紹的延遲時的覆蓋外推估算法。延遲時的覆蓋外推法通過使用一種幻影技術，在整個測線上可以把拾取初至時間擴展成一個連續(xù)的初至時間序列。這些擴展出來的拾取時間叫做重疊時間。對要求進行延遲時估算的區(qū)域，在最小平方的意義下，我們把這些重疊時間喻為和這個區(qū)域相鄰的站點的延遲時。然后，采用最小平方擬合法，把這些重疊時間值擴展到那個區(qū)域，作為那個區(qū)域的延遲時。這種技術對于外推和內插是一樣的。重疊時間估算法使用單個方向的拾取時間值去估算延遲時。位于估算區(qū)域內的站點的拾取時間可免于錯誤的拾取。在許多情況下，由于炮點不連續(xù)，重疊時間不能有效的估算整個區(qū)域的延遲時。在這種情況下，用這些重疊時間盡可能去估算這個區(qū)域必要的一些延遲時。剩下的估算既可用延遲時投影法也可用折射面高程投影法。如果，在要求進行非互易估算的區(qū)域里，初至時間比較可靠，并且，不想采用Hybrid EGR迭代算法，推薦使用重疊型的外推法。 分析員提供延遲時EGRM模塊允許由分析員來提供測線上任意段、任意或所有折射層的延遲時。如果分析員是通過截距法來提供延遲時的話，對任何解釋出的截距時都要除以2，因為截距時實際上是雙程的延遲時。如果提供了某一層的延遲時，并且分析員不希望用EGRM去估算這個折射層，那么將不應再提供這個折射面的偏移距范圍。 由分析員提供折射面高程，轉換成延遲時分析員可以提供折射面高程來轉換成延遲時。如果在第一次運算時，分析員觀察到EGRM把延遲時外推到一個區(qū)段時，得到的結果不合理，導致這個區(qū)段有一個看起來很奇怪的折射模型。那么，分析員需要重新運行EGRM，實驗不同的外推方案，以得到一個更可靠的結果。分析員可以使用TimeDepth Override 和 Refractor Elevation Override 兩者之一來指定外推法。 由井口數(shù)據(jù)提供延遲時可用井口時間來調整延遲時和速度。在這種方式下，埋在地下的炮點的海拔高程用來定義成一個偽折射面。后續(xù)處理中，當把實際的折射面的延遲時轉換成高程時，偽折射面被作為一個實際的折射面來考慮。如果由井口信息來定義層位可行，那么由井定義的偽折射層得到的靜校正結果，比那種只是把炮點移到地表得到的結果更準確。根據(jù)井口信息對層1的延遲時的調整，可根據(jù)實際情況采用不同的方式。（也就是，在模塊中的