【正文】 速度控制輸出（Velocity Diagnostic Printout）因?yàn)樗俣仁钦凵潇o校正計(jì)算中非常重要的方面，EGRM提供了兩級(jí)控制輸出。 速度當(dāng)使用速度繪圖（Plot Velocities）參數(shù)時(shí)，在折射速度分析結(jié)束時(shí)，得到的速度可繪出來(lái)。圖21，22是拾取時(shí)間繪圖的例子。然而，海底以下的那些折射面仍然是比較光滑的，延遲時(shí)的差異是海底到折射面之間的地層速度突變引起的。這些被舍棄的速度值可以要求熟悉該地區(qū)的地質(zhì)家來(lái)提供。 海洋資料處理要考慮的因素標(biāo)準(zhǔn)的海洋采集的觀測(cè)系統(tǒng)和陸上采集有一個(gè)根本的不同：在海洋采集中，由于電纜的漂移，檢波點(diǎn)位置是不固定的。例如，考慮下列已給定了偏移距的這些炮點(diǎn)：這些偏移距范圍必須變成相對(duì)于中心點(diǎn)的偏移距范圍。實(shí)際上，假如井穿過(guò)了折射層，這種射線路徑校正將是錯(cuò)誤的。（EGRM延遲時(shí)算法要求折射旅行時(shí)來(lái)自于地表），缺省的話，EGRM將對(duì)炮點(diǎn)進(jìn)行校正。如果分析員不希望EGRM去計(jì)算某一層的延遲時(shí)和該層下的折射面速度，那么就不應(yīng)該給出該層的偏移距范圍。圖17和圖18列舉了兩種常見情況下中間基準(zhǔn)面的選取。這些延遲值的正確性是未知的，它依賴于V1和V2的準(zhǔn)確估算。下面的例子可以看到確實(shí)需要這樣一個(gè)中間基準(zhǔn)面。分析員應(yīng)該認(rèn)識(shí)到，在本章討論的refraction elevation overrids和折射面高程平滑這兩者實(shí)際上都是一個(gè)對(duì)近地表的認(rèn)識(shí)和解釋過(guò)程。我們建議采用25點(diǎn)的平滑半徑，這樣將只改進(jìn)靜校正的高頻分量而不影響長(zhǎng)波長(zhǎng)分量。但是，我們不提倡這樣做。折射面高程的應(yīng)用有兩種方式。圖15 由圖13中的延遲時(shí)轉(zhuǎn)換成的折射模型，表1：由圖14中模型得到的靜校正 表2：由圖15中模型得到的靜校正在表1和表2 中，使用了相同的延遲時(shí)和折射層速度。假定地表是水平的， EGRM還計(jì)算出的折射層的速度是4000m/s。把延遲時(shí)轉(zhuǎn)換成折射面高程對(duì)于單風(fēng)化層模型，延遲時(shí)公式11和16可被表示成三個(gè)分量：折射面深度（ZG），速度（V1）和速度（V2），Q是arcsin（v1/v2）。不同的風(fēng)化層速度可以得到不同的折射模型和。直達(dá)波是在地表或近地表中傳播、直接到達(dá)檢波器的那部分地震能量，同樣的，沿直達(dá)波斜率畫線，在零偏移距處它們應(yīng)該交到零時(shí)間上。當(dāng)由井口得到的延遲時(shí)超過(guò)EGRM計(jì)算出的延遲時(shí)，井一定打穿了折射面。偽折射層選件中第三個(gè)方法（速度選件）提供了一個(gè)避免這種情況的方法。EGRM用這些井口信息把深井炮向上校正到地表。這些井的測(cè)量信息是否足夠用來(lái)定義一條線上的風(fēng)化層速度剖面，要根據(jù)井的間隔大小來(lái)看。在進(jìn)行這種測(cè)量時(shí)，每口井要打穿折射層，垂向的旅行時(shí)間測(cè)量值通常隨著每口井的深度的不同而變化。這三種方案在下面的章節(jié)討論?？刹捎闷渌椒▉?lái)求臨界距離。V1可以取任意值，只要它不超過(guò)V2就行，否則將不能得到折射波。TD（G）和V2已知，但是風(fēng)化層速度V1和風(fēng)化層厚度ZG還是未知的。當(dāng)井穿過(guò)一個(gè)實(shí)際存在的折射面時(shí)，我們認(rèn)為井太深了。 由井口數(shù)據(jù)提供延遲時(shí)可用井口時(shí)間來(lái)調(diào)整延遲時(shí)和速度。如果，在要求進(jìn)行非互易估算的區(qū)域里，初至?xí)r間比較可靠，并且，不想采用Hybrid EGR迭代算法，推薦使用重疊型的外推法。對(duì)要求進(jìn)行延遲時(shí)估算的區(qū)域，在最小平方的意義下，我們把這些重疊時(shí)間喻為和這個(gè)區(qū)域相鄰的站點(diǎn)的延遲時(shí)。如果認(rèn)為該地區(qū)（根據(jù)地質(zhì)信息或觀察）折射面基本上隨地表而起伏，推薦使用該選項(xiàng)。對(duì)不在測(cè)線端點(diǎn)處的那些延遲時(shí)，延遲時(shí)的估算使用的是平滑技術(shù)。但是，在測(cè)線的端點(diǎn)和炮點(diǎn)跳的比較大的區(qū)域，并不能得到這樣三個(gè)旅行時(shí)。Hybrid算法遵循等式17。因?yàn)镋GRM方法對(duì)任一點(diǎn)的延遲時(shí)估算它都需要三個(gè)折射波旅行時(shí)，在測(cè)線的端點(diǎn)或者測(cè)線中間炮點(diǎn)跳度比較大的時(shí)候，這種算法并不總能得到足夠的折射波旅行時(shí)信息來(lái)估算延遲時(shí)。在偏移距/ Vn+1這一項(xiàng)中產(chǎn)生的所有誤差都被分配到炮點(diǎn)延遲時(shí)項(xiàng)和檢波點(diǎn)延遲時(shí)項(xiàng)里。因?yàn)檎凵鋵铀俣萔n+1已經(jīng)估算出，偏移距—速度項(xiàng)是已知的，根據(jù)初始的結(jié)構(gòu)解，檢波點(diǎn)的延遲時(shí)值TB也是已知的，我們就可以重新計(jì)算炮點(diǎn)延遲時(shí)TA，對(duì)測(cè)線上的每一炮都這樣做。對(duì)測(cè)線上每一點(diǎn)都這樣做，在這些中心站點(diǎn)上得到的延遲時(shí)的粗解稱做初始的結(jié)構(gòu)解。對(duì)于每一個(gè)地面點(diǎn)，EGRM根據(jù)給定的折射分層，抽取出所有以該點(diǎn)為中心點(diǎn)的折射旅行時(shí)。 EGRM延遲時(shí)計(jì)算—兩個(gè)或更多風(fēng)化層情況圖6 兩層折射模型圖6中兩層的情況下，根據(jù)折射原理，折射傳播時(shí)間TAB可描述如下：這里 DA1=地表A點(diǎn)和折射面1上W點(diǎn)之間的射線距離 DA2=折射面1上W點(diǎn)和折射面2上X點(diǎn)之間的射線距離 DB1=地表B點(diǎn)和折射面1上Z點(diǎn)之間的射線距離 DB2=折射面1上Z點(diǎn)和折射面2上Y點(diǎn)之間的射線距離 DR=X點(diǎn)和Y點(diǎn)之間的射線距離 V1=地表和折射面1之間的平均速度 V2=折射面1和2之間的平均速度 V3=折射面2下伏地層的速度等式12可簡(jiǎn)化為：這里 ZA1= A點(diǎn)下方層1的厚度 ZA2= A點(diǎn)下方層2的厚度 ZB1= B點(diǎn)下方層1的厚度 ZB2= B點(diǎn)下方層2的厚度 AB= 地表A和B點(diǎn)之間的距離 a= 穿過(guò)折射層1 的角度，定義為arcsin（v1/v3） 216。這就是EGRM中基本的折射層速度計(jì)算方法。接下來(lái)敘述怎樣去分析折射波旅行時(shí)并計(jì)算出折射波速度和延遲時(shí)。=由斯奈兒定理定義的臨界角=arcsin（v1/v2）。第三步是根據(jù)折射模型計(jì)算靜校正。 折射波原理和折射波靜校正利用折射波求取靜校正稱做折射靜校正。這種算法依賴于速度項(xiàng)的估算精度，所以它對(duì)速度項(xiàng)計(jì)算的誤差更為敏感。在EGRM模塊中提供了三種方法來(lái)提取折射波的信息，即EGRM法、GaussSeidel迭代法和Hybrid EGRM/迭代法。d. 檢波點(diǎn)表格（Station Tables）指的是Primary Stations和Secondary Stations兩個(gè)參數(shù)表。b. 如果沒有替換速度和地表高程，本模塊就不會(huì)進(jìn)行靜校正計(jì)算，因此，如果野外沒有提供靜校正，也沒提供替換速度和地表高程，作業(yè)會(huì)報(bào)錯(cuò)。用戶提供的野外靜校正量在Primary shots或Secondary Shots或Primary Stations或Secondary Stations中輸入，如果上述表格中有有效靜校正量輸入，本模塊不重新計(jì)算靜校正量，如果上述表格中沒有靜校正量輸入，則本模塊自動(dòng)計(jì)算。在Primary Shots，Secondary Shots或ShotPattern Mapping三個(gè)參數(shù)之一中定義了Pattern號(hào)。如果原始觀測(cè)系統(tǒng)中存在的一些信息在流程中定義了，則觀測(cè)系統(tǒng)保持不變。Geometry Database選項(xiàng)中的Database Write Mode參數(shù)，就是控制創(chuàng)建表格信息方式的。5．如果一個(gè)參數(shù)僅在下列一個(gè)表格中定義（Primary Stations,Primary Shots,Secondary Stations,Secondary Shots），那么，這個(gè)參數(shù)值可以按上述規(guī)則進(jìn)行計(jì)算，如果參數(shù)同時(shí)在Primary和Secondary兩個(gè)表格中出現(xiàn)，Primary表格中參數(shù)按以上規(guī)則計(jì)算，而Secondary表格中參數(shù)不再進(jìn)行內(nèi)插或外推計(jì)算，只是Secondary表格中定義的參數(shù)值取代Primary表格中對(duì)應(yīng)位置的參數(shù)值，也就是說(shuō)Secondary表格是用來(lái)描述異常點(diǎn)的。6．Station Number Increment站號(hào)增量7．Trace Number Increment道增量GEOMETRY_LAND模塊用上述地球物理語(yǔ)言文件來(lái)構(gòu)件這些表格，并通過(guò)擴(kuò)展、內(nèi)插、外推的方法產(chǎn)生盡可能多的信息，因此處理員并不需要填寫一些多余的信息。 TOFS = SAMP_INT+TSA其中 TOFS = 第一個(gè)樣點(diǎn)時(shí)間 SAMP_INT = 數(shù)據(jù)采樣間隔 TSA = Time Shift Alignment12．Source Type震源的激發(fā)方式。3．Ident Number炮點(diǎn)標(biāo)識(shí)號(hào)4．XY Coordinates炮點(diǎn)的xy坐標(biāo)5．Surface Elevation炮點(diǎn)的地表高程6．Datum Elevation炮點(diǎn)的基準(zhǔn)面高程（參照海平面）7．Static Correction炮點(diǎn)靜校正量。本模塊不能夠?qū)缀闻帕行畔⒅苯哟娣诺降卣饠?shù)據(jù)道頭中去，這一點(diǎn)需要通過(guò)模塊Geometry_Updata來(lái)實(shí)現(xiàn)。參與此項(xiàng)工作的籌劃與組織者有：龔仁彬、華元生、王宇超、王小衛(wèi)、張建偉、王建華、胡自多、王紅旗等。本次翻譯是根據(jù)（版本：）各模塊的重要性有選擇地翻譯了39個(gè)模塊，約30萬(wàn)字 。 模塊簡(jiǎn)介設(shè)計(jì)GEOMETRY_LAND模塊的目的就是將野外地震隊(duì)提供的觀測(cè)系統(tǒng)的信息編制成規(guī)則的，與其一一對(duì)應(yīng)的表格信息，這樣使用起來(lái)更加方便，本模塊產(chǎn)生的規(guī)則觀測(cè)系統(tǒng)信息包括以下表格信息：1．檢波點(diǎn)表格（Detector Table）：描述地面上各站點(diǎn)的情況；2．炮點(diǎn)表格（Source Tables）：描述每一炮的信息；3．Pattern Tables：描述測(cè)線的炮檢關(guān)系；4．Exception Tables：描述編輯的一些信息。（如果野外沒有提供靜校正量值，那么本模塊根據(jù)炮點(diǎn)的地表高程、基準(zhǔn)面高程、替換速度、井口時(shí)間、風(fēng)化層厚度、風(fēng)化層速度、xy坐標(biāo)和General parameter set中Spatial Parameter Interpolation Type參數(shù)中定義的內(nèi)插類型等自動(dòng)計(jì)算出炮點(diǎn)的靜校正量；如果野外提供了靜校正量，則這個(gè)量作為首選，本模塊就不再計(jì)算新的靜校正量了， Field Static Corrections）。這個(gè)值主要用于計(jì)算野外靜校正時(shí)用的。另外，這些地球物理語(yǔ)言文件包括的參數(shù)要么是描寫地震測(cè)線炮檢點(diǎn)的空間位置關(guān)系，要么是通過(guò)一些已知定義的信息來(lái)計(jì)算出這些空間位置關(guān)系。（ Primary and Secondary Data Tables）Table 1和Table 2定義了炮點(diǎn)表格和檢波點(diǎn)表格中各個(gè)道頭字的求取順序。如果此參數(shù)選用NEW選鍵，就意味著要?jiǎng)?chuàng)建的觀測(cè)系統(tǒng)原先并不存在，模塊中定義的表格信息就會(huì)寫到新建的觀測(cè)系統(tǒng)中去。流程中定義了一些觀測(cè)系統(tǒng)中沒有的信息，則這些信息被寫到觀測(cè)系統(tǒng)中去。如果ShotPattern Mapping中定義了Pattern號(hào)，在Primary Shots或Secondary Shots中也定義了Pattern號(hào)，則ShotPattern Mapping中的Pattern值將Primary Shots或Secondary Shots中的Pattern值覆蓋掉，定義Pattern有兩種方式：relative（相對(duì)的）或absolute（絕對(duì)的）兩種方式。參照table 1和table 2。c. 檢波點(diǎn)表格和炮點(diǎn)表格都可以定義炮點(diǎn)的井口時(shí)間和井深，但不一定在兩個(gè)表格中都定義。e. 炮點(diǎn)表格（Shot Tables）指的是Primary Shots和Secondary Shots兩個(gè)參數(shù)表。這些技術(shù)既可以單獨(dú)使用，也可以組合使用。但是，在外推方面，迭代分解法要優(yōu)于EGRM法，而且迭代分解法可以分別計(jì)算炮、檢點(diǎn)的延遲時(shí)，而EGRM卻不能。折射靜校正的估算分三個(gè)步驟完成。首先，從地表開始，逐層下剝，計(jì)算垂直穿過(guò)該層的時(shí)間延遲，然后，從最深的折射面校正到一個(gè)基準(zhǔn)面或中間基準(zhǔn)面，校正速度為最深的折射層速度。在圖1中是A點(diǎn)到X點(diǎn)之間的射線與垂直于折射面上X點(diǎn)的直線之間的角度。 折射層速度分析圖2 這個(gè)折射模型顯示了從A點(diǎn)到G1，G2，G3，G4，G5點(diǎn)的折射路徑，折射路徑穿過(guò)n個(gè)折射層，在第n層的底達(dá)到臨界角，產(chǎn)生折射。在Hybrid EGRM/迭代算法中，另一個(gè)方法也可進(jìn)行速度計(jì)算。= 穿過(guò)折射層2 的由SHELL定理定義的臨界角 arcsin（v2/v3）圖7 兩個(gè)風(fēng)化層的折射模型象單風(fēng)化層情況一樣，也需要三個(gè)折射旅行時(shí)來(lái)計(jì)算延遲時(shí)，并且使用同樣的EGRM公式。例如，在120站點(diǎn)，一個(gè)這樣的旅行時(shí)可以是從炮點(diǎn)110到檢波點(diǎn)130的旅行時(shí)，因?yàn)檎军c(diǎn)120是這個(gè)炮、檢點(diǎn)對(duì)的中心點(diǎn)或中心站點(diǎn)，其它的旅行時(shí)如從炮點(diǎn)134到檢波點(diǎn)106，從炮點(diǎn)101到檢波點(diǎn)139等等。 平均截距時(shí)間圖10 應(yīng)用線性校正后的旅行時(shí)和偏移距對(duì)應(yīng)圖。所有的炮點(diǎn)延遲時(shí)重新計(jì)算后，用等式19，把所有的檢波點(diǎn)延遲時(shí)也用同樣的方式重新計(jì)算一遍。偏移距越大，Vn+1的誤差引起的影響就越大。 用非互易估算法估算延遲時(shí)。對(duì)于任意的拾取時(shí)間TAB，當(dāng)給出一系列TA和TB的估算值（從EGRM的結(jié)果中得到）后，可以重新計(jì)算出折射層速度Vn+1，這個(gè)速度和EGRM法所估算出的速度趨勢(shì)是一致的。同樣，如果在測(cè)線上部分地方由于拾取錯(cuò)誤或折射傳播的非地表一致性，會(huì)使初至拾取時(shí)間表現(xiàn)出非地表一致性。平滑技術(shù)用來(lái)估算那些變化量超過(guò)一定值的個(gè)別延遲時(shí)。這是因?yàn)槠交途€形外推所得到的延遲時(shí)是比較圓滑的，在建立近地表模型時(shí)，轉(zhuǎn)換成的折射層厚度（不是折射面高程）也是圓滑的，這樣，當(dāng)?shù)乇砥鸱鼤r(shí)，折射面高程也就隨之起伏。然后，采用最小平方擬合法，把這些重疊時(shí)間值擴(kuò)展到那個(gè)區(qū)域，作為那個(gè)區(qū)域的延遲時(shí)。 分析員提供延遲時(shí)EGRM模塊允許由分析員來(lái)提供測(cè)線上任意段、任意或所有折射層的延遲時(shí)。在這種方式下，埋在地下的炮點(diǎn)的海拔高程用來(lái)定義成一個(gè)偽折射面。在這種情況下，程序?qū)φ凵涿?的求解強(qiáng)行終止，并且相當(dāng)多的風(fēng)化層信息被丟失。圖11（來(lái)自于一個(gè)簡(jiǎn)單模型的一系列拾取時(shí)間的圖形）