靜校正(討論版)

1、四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司靜校正處理靜校正處理（黃大云）（黃大云）2008.4四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司內(nèi)容內(nèi)容 靜校正概述靜校正概述 高程靜校正高程靜校正 折射靜校正折射靜校正 層析靜校正層析靜校正 剩余靜校正剩余靜校正四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司地震資料常規(guī)處理一般流程地震資料常規(guī)處理一般流程數(shù)據(jù)輸入數(shù)據(jù)輸入觀測系統(tǒng)加載觀測系統(tǒng)加載壞炮壞道去除壞炮壞道去除噪聲壓制噪聲壓制振幅補償振幅補償靜校正靜校正反褶積反褶積疊加疊加剩余靜校正剩余

2、靜校正速度分析速度分析動校正動校正偏移偏移四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 問題問題:在地震資料處理中為什么要做靜校正在地震資料處理中為什么要做靜校正?靜校正概述靜校正概述 1、 為使地下構造形態(tài)不受地表起伏的影響為使地下構造形態(tài)不受地表起伏的影響; 2、一個探區(qū)需使用統(tǒng)一的基準面、一個探區(qū)需使用統(tǒng)一的基準面,以便測線之間接點以便測線之間接點; 3、地震資料處理的一些重要步驟是在反射波時距曲線為雙、地震資料處理的一些重要步驟是在反射波時距曲線為雙 曲線的前提下進行的（速度分析、動校正等），但反射曲線的前提下進行的（速度分析、動校正等），但反射 波時距曲線為雙曲線的條

3、件是：地表水平、上覆介質(zhì)速度波時距曲線為雙曲線的條件是：地表水平、上覆介質(zhì)速度 為常數(shù)。為了后續(xù)處理的需要，應把反射時間校正到炮點、為常數(shù)。為了后續(xù)處理的需要，應把反射時間校正到炮點、 檢波點均在一個水平面上一樣的情況。檢波點均在一個水平面上一樣的情況。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司22202vxttP PV VS SM MxR R021txt21（雙曲線方程）（雙曲線方程）v地表水平、上覆介質(zhì)速度為常數(shù)時地表水平、上覆介質(zhì)速度為常數(shù)時,地震波傳播路徑如下圖地震波傳播路徑如下圖:四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 在地震勘探中，由于受地形

4、起伏、爆炸井深不一、低在地震勘探中，由于受地形起伏、爆炸井深不一、低降速帶的厚度和速度變化等因素的影響，地震波的反射時降速帶的厚度和速度變化等因素的影響，地震波的反射時距曲線已不是理論上的雙曲線分布。如何對記錄的反射時距曲線已不是理論上的雙曲線分布。如何對記錄的反射時間進行校正，使反射時距曲線還原為雙曲線，這就是靜校間進行校正，使反射時距曲線還原為雙曲線，這就是靜校正要解決的主要問題。靜校正通常包括三個方面：正要解決的主要問題。靜校正通常包括三個方面：1 1、野外高程變化引起的時差校正；、野外高程變化引起的時差校正；2 2、低、降速帶引起的時差校正；、低、降速帶引起的時差校正；3 3、以上兩種

5、校正留下的殘余和其它因素引起的剩余時差校正。、以上兩種校正留下的殘余和其它因素引起的剩余時差校正。 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司目前我們常用的靜校正方法目前我們常用的靜校正方法高程靜校正高程靜校正折射靜校正折射靜校正層析靜校正層析靜校正剩余靜校正剩余靜校正四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 解決地形起伏、爆炸井深不一引起的靜校正問解決地形起伏、爆炸井深不一引起的靜校正問題可用野外高程校正方法解決。該方法利用野外測題可用野外高程校正方法解決。該方法利用野外測量成果和預定的基準面高程及其上覆層速度來計算量成果和預定的基準面高程及其上覆層速度

6、來計算校正量。校正量。高程靜校正高程靜校正四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司P PV VS1S1R1R1M M基準面地面S SR R炮點校正量檢波點校正量高程校正示意圖高程校正示意圖誤差來自于把三角形的一個邊誤差來自于把三角形的一個邊的長度看作另兩邊之和，其誤的長度看作另兩邊之和，其誤差隨基準面離地面的距離和炮差隨基準面離地面的距離和炮檢距的增加而變大。檢距的增加而變大。這種計算校正量的這種計算校正量的方法有時誤差太大方法有時誤差太大四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 通常不直接校正到標準基準面而是先校正到通常不直接校正到標準基準面而是先校正

7、到浮動浮動基準面，以后再校正到標準基準面?；鶞拭妫院笤傩Ｕ綐藴驶鶞拭?。問題問題: :1 1、為什么要使用浮動基準面、為什么要使用浮動基準面? ?2 2、在什么情況下須使用、在什么情況下須使用浮動浮動基準面基準面? ?3 3、數(shù)據(jù)在什么時候校正到固定基準面、數(shù)據(jù)在什么時候校正到固定基準面? ?四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司P PV VS1S1R1R1M M基準面地面S SR R炮點校正量檢波點校正量高程校正示意圖高程校正示意圖浮動基準面S2S2R2R2 浮動基準面與地面浮動基準面與地面高差較小高差較小, ,可降低靜可降低靜校正計算誤差校正計算誤差問題問題1:為什

8、么要使用浮動基準面為什么要使用浮動基準面?四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 浮動基準面在局部范圍內(nèi)可看作是平面，它可能是水平面，也可能是傾浮動基準面在局部范圍內(nèi)可看作是平面，它可能是水平面，也可能是傾斜面。如果是水平面，靜校正的目的已經(jīng)達到。如果是傾斜面，則反射波斜面。如果是水平面，靜校正的目的已經(jīng)達到。如果是傾斜面，則反射波走時圖如左圖。我們將該圖做旋轉(zhuǎn)，使傾斜的浮動基準面變?yōu)樗?，原本走時圖如左圖。我們將該圖做旋轉(zhuǎn)，使傾斜的浮動基準面變?yōu)樗?，原本水平的反射界面就變?yōu)閮A斜。這樣一來，我們就把傾斜基準面的問題變成水平的反射界面就變?yōu)閮A斜。這樣一來，我們就把傾斜基準

9、面的問題變成了傾斜反射面的問題。而了傾斜反射面的問題。而傾斜反射面的時差校正傾斜反射面的時差校正將將由傾斜反射面的動校正去解決。由傾斜反射面的動校正去解決。 如果地下反射界面原本就是傾斜的，旋轉(zhuǎn)后變成什么樣呢？如果地下反射界面原本就是傾斜的，旋轉(zhuǎn)后變成什么樣呢？SRRS 靜校正的目不是要將炮點、檢波點校正到同一水平面上嗎？靜校正的目不是要將炮點、檢波點校正到同一水平面上嗎？只校正到浮動基準面，動校正的條件滿足嗎？只校正到浮動基準面，動校正的條件滿足嗎？四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 通常用炮點、檢波點的平滑面作為浮動基準面，平滑長度通常用炮點、檢波點的平滑面作為浮

10、動基準面，平滑長度一般為排列長度。浮動基準面靜校正只計算炮點、檢波點到浮一般為排列長度。浮動基準面靜校正只計算炮點、檢波點到浮動基準面的校正量，以后再計算浮動基準面到標準基準面的校動基準面的校正量，以后再計算浮動基準面到標準基準面的校正量。炮點、檢波點實際靜校正量等于兩部分相加。正量。炮點、檢波點實際靜校正量等于兩部分相加。 當?shù)乇砥鸱髸r須使用浮動基準面當?shù)乇砥鸱髸r須使用浮動基準面!問題問題3:數(shù)據(jù)在什么時候校正到固定基準面數(shù)據(jù)在什么時候校正到固定基準面?問題問題2:在什么情況下須使用浮動基準面在什么情況下須使用浮動基準面?四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司三

11、種可能的選擇三種可能的選擇: 1、動校前、動校前 2、動校后、動校后 3、疊加后、疊加后1 1、動校前、動校前: :與直接使用固定基準面差別不大與直接使用固定基準面差別不大2 2、動校后、動校后, ,即速度分析和動校正在浮動基準面上進行即速度分析和動校正在浮動基準面上進行: :浮動基準面是平滑的浮動基準面是平滑的, ,在局部范圍內(nèi)可看成直線在局部范圍內(nèi)可看成直線, ,這時的這時的靜靜校校正相當于傾斜地層的動校正正相當于傾斜地層的動校正,傾角為浮動基準面傾角和地層傾角為浮動基準面傾角和地層傾角的綜合傾角的綜合。動校后動校后CMP道集內(nèi)同相軸拉平道集內(nèi)同相軸拉平,此時校正到固此時校正到固定基準面后

12、同相軸會有時差定基準面后同相軸會有時差,但因浮動基準面比地表平滑但因浮動基準面比地表平滑,時時差不及差不及直接使用固定基準面那么大直接使用固定基準面那么大, ,有時有時或許或許可用剩余靜?？捎檬Ｓ囔o校正解決。正解決。3 3、疊加后、疊加后: :理論上疊加成像最好理論上疊加成像最好, ,但這時已沒有炮點但這時已沒有炮點、檢波點的檢波點的 概念概念, ,校正須按校正須按CMPCMP點進行。點進行。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 野外高程校正中的替換速度為常數(shù)，當探野外高程校正中的替換速度為常數(shù)，當探區(qū)內(nèi)無低降速帶變化或變化很小可忽略不計時，區(qū)內(nèi)無低降速帶變化或變化很小

13、可忽略不計時，使用野外高程校正即可解決靜校正問題。但對使用野外高程校正即可解決靜校正問題。但對于復雜地區(qū)而言，情況就不那么簡單，因為在于復雜地區(qū)而言，情況就不那么簡單，因為在這些地區(qū)，低降速帶的厚度和速度通常都存在這些地區(qū)，低降速帶的厚度和速度通常都存在變化，這時就不能使用高程靜校正了。變化，這時就不能使用高程靜校正了。 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司折射靜校正折射靜校正 如果低降速帶存在，則它與高速層的分界面如果低降速帶存在，則它與高速層的分界面是一個很好的反射界面。當入射波以臨界角入射是一個很好的反射界面。當入射波以臨界角入射到該界面時，將產(chǎn)生折射波。折射波傳

14、到地面為到該界面時，將產(chǎn)生折射波。折射波傳到地面為檢波器所接收，根據(jù)折射波的初至時間可估算低檢波器所接收，根據(jù)折射波的初至時間可估算低降速帶的速度和炮點及檢波點的延遲時間，這就降速帶的速度和炮點及檢波點的延遲時間，這就是折射靜校正的基本思路。是折射靜校正的基本思路。 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司A點激發(fā)，點激發(fā)，B點接收的折射波到達時間點接收的折射波到達時間TAB為：為： TAB = AC/VR+CD/VS+DB/VR 或或 TAB = AH/（VR cos）+AB-（AH+BG）tan/VS+BG/（VR cos） 考慮到考慮到sin=VR/VS ,且且VR

15、VS , 上式變?yōu)樯鲜阶優(yōu)?TAB = TA +AB/ VS + TB （折射波旅行時方程）（折射波旅行時方程） 其中，其中，TA = AHcos/VR 為炮點的延遲時間為炮點的延遲時間 (時間深度時間深度)， TB = BGcos/VR 為檢波點的延遲時間為檢波點的延遲時間(時間深度時間深度) 。 對以上對以上折射波旅行時方程可用互換法或迭代法求解。折射波旅行時方程可用互換法或迭代法求解。 折射示意圖折射示意圖 臨界角：sin=VR/ VS 地層速度：VSABHGCD低降速帶速度：VR地面折射面四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司因為臨界角因為臨界角：當?shù)徒邓賻俣犬?/p>

16、低降速帶速度 與地層速度與地層速度 差別不大時差別不大時,臨界角會臨界角會很大很大,在有限的排列內(nèi)可能收不到折射或折射波到達時在有限的排列內(nèi)可能收不到折射或折射波到達時,其其它波已先行到達它波已先行到達, 所以低降速帶速度與地層速度的差別必所以低降速帶速度與地層速度的差別必須足夠大才能在記錄上拾取到折射波初至。須足夠大才能在記錄上拾取到折射波初至。 SRVV /sinRVSV問題問題1:在什么情況下能在記錄上拾取到折射波初至在什么情況下能在記錄上拾取到折射波初至?四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 折射靜校正的兩個必要條件：折射靜校正的兩個必要條件： 1 1、存在穩(wěn)定

17、的折射層；、存在穩(wěn)定的折射層； 2 2、能夠從炮記錄的初至中拾取到折射波初至。、能夠從炮記錄的初至中拾取到折射波初至。 如果不存在穩(wěn)定的折射層或者或不能拾取到折射波初至，如果不存在穩(wěn)定的折射層或者或不能拾取到折射波初至，則無法實現(xiàn)折射靜校正。則無法實現(xiàn)折射靜校正。 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司問題問題2:折射波初至在記錄上有什么特征折射波初至在記錄上有什么特征?為什么為什么? 與直達波初至比較與直達波初至比較, ,相鄰道時差小相鄰道時差小, ,所以顯得平直所以顯得平直, ,兩端與直達波相接形成拐點。兩端與直達波相接形成拐點。 因為相鄰道折射波到達時間之間的差由它

18、們傳因為相鄰道折射波到達時間之間的差由它們傳播路徑不同引起的播路徑不同引起的, ,而不同的這段路徑為高速。而不同的這段路徑為高速。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司Omega系統(tǒng)的折射靜校正模塊系統(tǒng)的折射靜校正模塊 EGRM EGRMEGRM是是Extnded Generalized Reciprocal Method(Extnded Generalized Reciprocal Method(擴展擴展常規(guī)互換法常規(guī)互換法) )的簡稱的簡稱, ,該模塊根據(jù)折射初至拾取時間使用互該模塊根據(jù)折射初至拾取時間使用互換法或換法或( (和和) )高斯高斯- -賽德爾迭代法求解獲

19、得炮點項和檢波點項賽德爾迭代法求解獲得炮點項和檢波點項的折射時間深度和折射速度層速度信息，并由此構建出折的折射時間深度和折射速度層速度信息，并由此構建出折射層的速度深度模型，再并根據(jù)速度深度模型計算出靜校射層的速度深度模型，再并根據(jù)速度深度模型計算出靜校正量。正量。 總體靜校正量總體靜校正量= 地面到中間基準面時間地面到中間基準面時間+ +中間基準中間基準面到標準基準面的時間面到標準基準面的時間以折射層速度計算以折射層速度計算以替換速度計算以替換速度計算四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司三種提取折射信息的方法三種提取折射信息的方法: :1 1、擴展常規(guī)互換法擴展常規(guī)互

20、換法,該技術是常規(guī)互換法技術的擴展該技術是常規(guī)互換法技術的擴展 能獲得準確有效的折射延遲時估算能獲得準確有效的折射延遲時估算,且與折射速度場的誤差相對無關且與折射速度場的誤差相對無關。2、高斯高斯- -賽德爾迭代分解法賽德爾迭代分解法( (最小平方法最小平方法) ) 有很強的外推特性有很強的外推特性, ,因該方法要同時估算炮點因該方法要同時估算炮點、檢波點折射延遲時和射檢波點折射延遲時和射速度速度, ,所以與折射速度場的誤差有關所以與折射速度場的誤差有關。3 3、混合、混合EGRM/EGRM/迭代算法迭代算法 以上兩種算法組合以上兩種算法組合。 三種方法既可單獨使用三種方法既可單獨使用, ,也

21、可按某種組合使用也可按某種組合使用。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司擴展常規(guī)互換法擴展常規(guī)互換法四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司互換法折射速度計算互換法折射速度計算AG1G2G3G4G5V1VnVn+1V1VnBG1G2G3G4G5Vn+1DT1=TAG1-TBG1= 負值負值DO1=AG1-BG1= 負值負值DT3=TAG3-TBG3= 0DO3=AG3-BG3= 0DT5=TAG5-TBG5= 正值正值DO5=AG5-BG5= 正值正值DT= 初至時差初至時差DO= 偏移距差偏移距差A、B：炮點：炮點G： 接收點接收點四川石油管理局地

22、球物理勘探公司四川石油管理局地球物理勘探公司DTDOG1G2G3G4G5+_Vn+1= 擬合直線斜率擬合直線斜率/2DT：時間差：時間差DO：偏移距差：偏移距差四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司互換法延遲時間互換法延遲時間(以時間表示的深度以時間表示的深度)計算計算XYABV1V2DADBDR121VDBVDRVDATAB當當V2比比V1大得多時大得多時,上式可近似地改為上式可近似地改為: （1）TAB=折射波到達時間；折射波到達時間； ZA=A點的折射層厚度；點的折射層厚度；ZA= B點的折射層厚度；點的折射層厚度； = A到到B點的地面距離；點的地面距離； =臨界

23、角。臨界角。這說明折射波到達時間可表示為這說明折射波到達時間可表示為A、B兩點的時間深度之和再加上兩點的時間深度之和再加上AB的距離的距離/V2。121VCOSZBVABVCOSZATABAB四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 某一站點某一站點G的時間深度的時間深度TD(G)需用需用A、B兩個炮點的三個折射旅行時兩個炮點的三個折射旅行時: （2） TAB= A點激發(fā)點激發(fā)B點接收的折射時間；點接收的折射時間； TAG= A點激發(fā)點激發(fā)G點接收的折射時間；點接收的折射時間； TBG= B點激發(fā)點激發(fā)G點接收的折射時間。點接收的折射時間。2)()(TABTGBTAGGTD

24、AGBV1臨界角臨界角V2四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 并非所有點都能使用互換法求取其時間深度，如測線兩端、炮點不規(guī)則并非所有點都能使用互換法求取其時間深度，如測線兩端、炮點不規(guī)則段等。對這些點必須使用其它辦法段等。對這些點必須使用其它辦法: :1 1、Time-depth ProjectTime-depth Project 利用相鄰點的時間深度做最小平方（線性）擬合利用相鄰點的時間深度做最小平方（線性）擬合2 2、Refractor Elevation ProjectRefractor Elevation Project 先將先將相臨點的時間深度轉(zhuǎn)換為折射層高

25、程，再做最小平方（線性）擬合相臨點的時間深度轉(zhuǎn)換為折射層高程，再做最小平方（線性）擬合3 3、Overlay Extrapolation of Time-DepthOverlay Extrapolation of Time-Depth 將拾取的初至時間依次連續(xù)擴展到整條測線，這些擴展時間就叫做重復時將拾取的初至時間依次連續(xù)擴展到整條測線，這些擴展時間就叫做重復時間（間（overlay-times），重復時間與要估算時間深度的區(qū)域的臨近站點的時），重復時間與要估算時間深度的區(qū)域的臨近站點的時間進行最小平方意義下的比較，并按最小平方擬合線擴展到估算區(qū)域。這種間進行最小平方意義下的比較，并按最小平方

26、擬合線擴展到估算區(qū)域。這種方法等同于外插和內(nèi)插地結合。方法等同于外插和內(nèi)插地結合。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司折射模型的建立折射模型的建立 為了建立折射模型，先將時間為了建立折射模型，先將時間-深度表達式做一下變換：深度表達式做一下變換： 將式（將式（1）代入式（）代入式（2），經(jīng)一系列計算并考慮到），經(jīng)一系列計算并考慮到 最后得到：最后得到： （3） 其中，其中，ZG為折射層厚度為折射層厚度,在式（在式（3）中，）中，ZG、V1和臨界角和臨界角 都是未知的，但臨界角可由都是未知的，但臨界角可由V1、V2求得。前求得。前面我們用直線擬合方式求得了面我們用直線擬合

27、方式求得了V2,所以實際上只有所以實際上只有ZG和和V1兩個未知數(shù)，而折射層兩個未知數(shù)，而折射層的速度的速度-深度正是由它們確定。為了得到深度正是由它們確定。為了得到ZG和和V1，我們再來看另一個觀測值：臨，我們再來看另一個觀測值：臨界距離，所謂臨界距離即折射波最先道達那一道的炮檢距。臨界距離界距離，所謂臨界距離即折射波最先道達那一道的炮檢距。臨界距離XC與臨界角與臨界角的關系是的關系是:GBAGAB1cos)(VZGGTD21sin1VV)4(tan2ZGXC四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司式（式（4）說明，在）說明，在V2已知的前提下，臨界距離已知的前提下，臨界

28、距離XC也是也是ZG和和V1的函數(shù)。于是由的函數(shù)。于是由 和和兩個方程可最終求得兩個方程可最終求得ZG和和V1。有了各站點（。有了各站點（G）的厚度和速度）的厚度和速度V1，折射層的速度，折射層的速度深度模型也就確定了。深度模型也就確定了。1cos)(VZGGTDtan2ZGXC四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司高斯高斯-賽德爾迭代旅行時分解賽德爾迭代旅行時分解 該方法就是用迭代分解的辦法計算各站點折射層的時間深度，有了時間深度，即可該方法就是用迭代分解的辦法計算各站點折射層的時間深度，有了時間深度，即可按照前面的方法得到折射層的速度深度模型。式（按照前面的方法得到折

29、射層的速度深度模型。式（1）表示折射旅行時是炮點、檢波點時間深度之和再加上偏移距表示折射旅行時是炮點、檢波點時間深度之和再加上偏移距/V2，為便于敘述，可寫成：，為便于敘述，可寫成：顯然迭代計算公式就是顯然迭代計算公式就是 其中，其中，TA =炮點炮點A的時間深度，的時間深度，TB=檢波點檢波點B的時間深度，的時間深度， 稱為構造項或偏移稱為構造項或偏移距距-速度項速度項,由于由于V2可通過前面的互換法求得，這一項可當作已知，這時迭代只需前兩可通過前面的互換法求得，這一項可當作已知，這時迭代只需前兩個公式即可，若使用三個公式迭代，則用互換法求得的構造項用作迭代初值。個公式即可，若使用三個公式迭

30、代，則用互換法求得的構造項用作迭代初值。121VCOSZBVABVCOSZATABTBVABTATAB2TBVABTABTA2TAVABTABTB2TBTATABVAB22VAB四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司TA和和TB的初值估算的初值估算 為了獲得某個站點的初始時間深度，模塊要把以該點為中點的所有折射旅行時收為了獲得某個站點的初始時間深度，模塊要把以該點為中點的所有折射旅行時收集起來：集起來： 如圖所示，要得到站點如圖所示，要得到站點120的初始時間深度，就把炮點的初始時間深度，就把炮點110到檢波點到檢波點130、炮點、炮點134到檢波點到檢波點106、炮點、

31、炮點101到檢波點到檢波點139等等折射旅行時收集起來，它們顯然都是等等折射旅行時收集起來，它們顯然都是偏移距偏移距-速度項的函數(shù)（左圖）。將偏移速度項的函數(shù)（左圖）。將偏移-速度項減去，再求平均，然后除以速度項減去，再求平均，然后除以2，即，即為該站點的初始時間深度（右圖）。為該站點的初始時間深度（右圖）。120101 106110130 134139TBVABTATAB2四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司混合混合EGRM/迭代算法迭代算法 用互換法得到的各站點的時間深度和（或）偏移用互換法得到的各站點的時間深度和（或）偏移距距- -速度值作為初始值迭代計算出各站點

32、最終的時間速度值作為初始值迭代計算出各站點最終的時間深度和（或）偏移距深度和（或）偏移距- -速度值速度值。 迭代法和混合法的折射層速度深度模型建立與前迭代法和混合法的折射層速度深度模型建立與前述互換法相同。述互換法相同。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司靜校正量計算靜校正量計算 計算靜校正量需用兩個基準面，一個是中間基準面計算靜校正量需用兩個基準面，一個是中間基準面, ,一個是最終基準一個是最終基準面。任一站點的靜校正量等于從該點到中間基準面的時間延遲和從中間面。任一站點的靜校正量等于從該點到中間基準面的時間延遲和從中間基準面到最終基準面時間延遲的代數(shù)和。從站點到中

33、間基準面的時間延基準面到最終基準面時間延遲的代數(shù)和。從站點到中間基準面的時間延遲計算使用前面所建立的速度深度模型的速度，從中間基準面到最終基遲計算使用前面所建立的速度深度模型的速度，從中間基準面到最終基準面時間延遲的計算使用替換速度。準面時間延遲的計算使用替換速度。 EGRMEGRM將中間基準面定為折射層底界的一個平滑面。因此，若折射層底將中間基準面定為折射層底界的一個平滑面。因此，若折射層底界高程變化不大，則中間基準面為水平，若折射層底界高程變化較大，界高程變化不大，則中間基準面為水平，若折射層底界高程變化較大，則中間基準面就不是水平的。另外用戶還可以用參數(shù)對中間基準面進行則中間基準面就不是

34、水平的。另外用戶還可以用參數(shù)對中間基準面進行修改或控制。修改或控制。 折射層的速度深度模型可以是單層也可以是多層，但即使我們只要折射層的速度深度模型可以是單層也可以是多層，但即使我們只要求單層模型，模塊也把井底以上的地表當作一層即所謂的虛擬折射層，求單層模型，模塊也把井底以上的地表當作一層即所謂的虛擬折射層，所以實際上仍是兩層。多層模型的靜值計算公式如下：所以實際上仍是兩層。多層模型的靜值計算公式如下：四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司多層模型的靜值計算公式多層模型的靜值計算公式其中，其中， TS = 總的靜校正量總的靜校正量 THICKi = 第第i 層的厚度層的厚

35、度 Vi = 第第 i層的速度層的速度 ERn = 第第 n層的高程層的高程 EI = 中間基準面高程中間基準面高程 Vn+1 = 第第 n層以下的速度層以下的速度 ED = 最終基準高程最終基準高程 VR = 替換速度替換速度 TUP = 井口時間井口時間UPRDInInniiiTVEEVEERVTHICKTS11四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司主要參數(shù)主要參數(shù)一一 處理參數(shù)組處理參數(shù)組(*PROCESS_OPTIONS)1 1、算法類型、算法類型(SOLUTION_TYPE)(SOLUTION_TYPE)三個選項三個選項: :EGRM:EGRM:所有折射層都用擴

36、展常規(guī)互換算法所有折射層都用擴展常規(guī)互換算法ITERATIVE:ITERATIVE:指定層使用迭代法指定層使用迭代法, ,其余用擴展常規(guī)互換算法其余用擴展常規(guī)互換算法HYBRID:HYBRID:使用擴展常規(guī)互換算法使用擴展常規(guī)互換算法/ /迭代法混合算法迭代法混合算法, ,即用擴展常規(guī)即用擴展常規(guī)互換算法求得的時間深度值作為初值進行迭代互換算法求得的時間深度值作為初值進行迭代 缺省缺省:EGRM:EGRM四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司2、最終基準面高程、最終基準面高程(DATUM_ELEV) 缺省缺省: #DATUM_ELEV3、替換速度、替換速度(REPL_VE

37、L) 缺省缺省: #REPLACE_VEL4、初至拾取時間調(diào)整值、初至拾取時間調(diào)整值(BULK_SHIFT) 所有初至拾取時間都加上該值所有初至拾取時間都加上該值 缺省缺省: 0.05、時間深度偏差門檻值、時間深度偏差門檻值(VARIANCE_THRESHOLDL) 計算偏差超過該門檻值的時間深度將做加權平滑計算偏差超過該門檻值的時間深度將做加權平滑,并在迭并在迭 代計算中不用作初值。代計算中不用作初值。 缺省缺省: 16ms四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司二二 折射層參數(shù)組折射層參數(shù)組(*LAYER_CONTROL) 該參數(shù)組以參數(shù)表的形式給出該參數(shù)組以參數(shù)表的形

38、式給出1 1、重復站點最小個數(shù)、重復站點最小個數(shù)(MIN_OVERLAP)(MIN_OVERLAP) 在測線兩端或炮點缺失區(qū)域不能獲得時間深度或折射高程在測線兩端或炮點缺失區(qū)域不能獲得時間深度或折射高程, ,需用最小平需用最小平方擬合外插獲得。本參數(shù)指定做方擬合外插獲得。本參數(shù)指定做最小平方擬合的站點個數(shù)。最小平方擬合的站點個數(shù)。 缺省缺省:25:252 2、折射模型的層號、折射模型的層號(LAYER)(LAYER) 缺省缺省: NOT_USED(: NOT_USED(無缺省無缺省) ) 3 3、用于該折射層處理的起始站點號、用于該折射層處理的起始站點號 缺省缺省:GENERAL_FIRST_

39、STN:GENERAL_FIRST_STN4 4、用于、用于該折射層處理該折射層處理的最后站點號的最后站點號 缺省缺省:GENERAL_LAST_STN:GENERAL_LAST_STN四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司5 5、非互換估算類型、非互換估算類型(ESTIMATION_TYPE)(ESTIMATION_TYPE) 指定生成測線兩端或其它不能使用互換法區(qū)域的時間深度的估算類型指定生成測線兩端或其它不能使用互換法區(qū)域的時間深度的估算類型 TIME_DEPTH TIME_DEPTH 用相鄰站點的時間深度做最小平方擬合外插用相鄰站點的時間深度做最小平方擬合外插 E

40、LEVATION ELEVATION 用相鄰站點的折射高程做最小平方擬合外插用相鄰站點的折射高程做最小平方擬合外插 EXTRAP_TIME_DEPTH EXTRAP_TIME_DEPTH 用相鄰站點的拾取時間做最小平方擬合外插用相鄰站點的拾取時間做最小平方擬合外插, ,然后然后 計算時間深度。計算時間深度。 EXTRAP_ELEV EXTRAP_ELEV 用相鄰站點的折射高程做最小平方擬合外插用相鄰站點的折射高程做最小平方擬合外插, , 然后計算時然后計算時 間深度。間深度。6 6、中間站點號、中間站點號(MID_STN)(MID_STN) 指定該層內(nèi)的站點號指定該層內(nèi)的站點號, ,它的最小偏

41、移距、它的最小偏移距、最大偏移距和內(nèi)插方式由后面的最大偏移距和內(nèi)插方式由后面的參數(shù)指定。參數(shù)指定。 FIRST FIRST 第一個站點第一個站點 LASTLAST 最后一個站點最后一個站點 其它數(shù)組其它數(shù)組 ( (無缺省無缺省) )四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司7 7、最小偏移距、最小偏移距(MIN_OFFSET)(MIN_OFFSET) 該中間站點號的最小偏移距該中間站點號的最小偏移距 大于等于大于等于0 0 缺省缺省:0.0:0.08 8、最大偏移距、最大偏移距(MIN_OFFSET)(MIN_OFFSET) 該中間站點號的最大偏移距該中間站點號的最大偏移距

42、大于等于最小偏移距大于等于最小偏移距 缺省缺省: :最小偏移距最小偏移距 由最小偏移距和最大偏移距確定該偏移距范圍內(nèi)的初至拾取時間將由最小偏移距和最大偏移距確定該偏移距范圍內(nèi)的初至拾取時間將 用作時間深度和折射速度計算。用作時間深度和折射速度計算。9 9、內(nèi)插偏移距、內(nèi)插偏移距 是否對下一個指定站點的偏移距做內(nèi)插是否對下一個指定站點的偏移距做內(nèi)插 YES YES 內(nèi)插內(nèi)插 NO NO 不內(nèi)插不內(nèi)插 插缺省插缺省: :YES 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司三三 迭代求解參數(shù)組迭代求解參數(shù)組(*ITERATIONS) 當算法類型選擇了當算法類型選擇了ITERATIVE

43、ITERATIVE時時, ,必須使用本參數(shù)組必須使用本參數(shù)組1 1、迭代次數(shù)、迭代次數(shù)(MUN_ITER)(MUN_ITER) 重復估算炮點和檢波點項的次數(shù)。重復估算炮點和檢波點項的次數(shù)。 缺省缺省:6:62 2、分解分量項順序、分解分量項順序 SD SD 炮點項炮點項, ,檢波點項檢波點項 DS DS 檢波點項檢波點項, ,炮點項炮點項 DSC DSC 檢波點項檢波點項, ,炮點項炮點項, ,構造項構造項 SDC SDC 炮點項炮點項, , 檢波點項檢波點項, ,構造項構造項 CDS CDS 構造項構造項, ,檢波點項檢波點項, ,炮點項炮點項 CSD CSD 構造項構造項, ,炮點項炮點項

44、, ,檢波點項檢波點項 缺省缺省: SD : SD 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司3 3、使用迭代法的折射層號、使用迭代法的折射層號 無缺省無缺省四四 混合迭代方法參數(shù)組混合迭代方法參數(shù)組(HYBRID_METHOD)(HYBRID_METHOD)1 1、刪除出現(xiàn)過大值的道、刪除出現(xiàn)過大值的道(DELETE_OUTLIERS)(DELETE_OUTLIERS) 刪除或保留那些在統(tǒng)計中出現(xiàn)過大值的道刪除或保留那些在統(tǒng)計中出現(xiàn)過大值的道 NO NO 不刪除不刪除( (建議建議) ) YES YES 刪除刪除 缺省缺省:NO:NO四川石油管理局地球物理勘探公司四川石油管

45、理局地球物理勘探公司2 2、過大值刪除門檻、過大值刪除門檻(OUTLIER_THRESHOLD)(OUTLIER_THRESHOLD) 任何小于該值的拾取時間都不當作過大值任何小于該值的拾取時間都不當作過大值 缺省缺省:5ms:5ms3 3、每個站點的炮點、檢波點時間深度分別計算、每個站點的炮點、檢波點時間深度分別計算(STN_DEP)(STN_DEP) YES YES 分別計算分別計算 NO NO 同一站點的炮點、檢波點時間深度同為它們的平均值同一站點的炮點、檢波點時間深度同為它們的平均值 缺省缺省:YES:YES4 4、混合算法的構造項、混合算法的構造項(SLACK_TERM)(SLACK

46、_TERM) 構造項是偏移距構造項是偏移距/ /折射速度折射速度, ,而偏移距是確定的。本參數(shù)指定構而偏移距是確定的。本參數(shù)指定構 造項是速度還是中間站點一致的構造項造項是速度還是中間站點一致的構造項 VELOCITYS VELOCITYS 折射速度折射速度 STRUCTURAL STRUCTURAL 構造項構造項 缺省缺省: VELOCITYS: VELOCITYS 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 五五 速度控制參數(shù)組速度控制參數(shù)組(VEL） 每個站點的地表速度必須指定，可用本參數(shù)組指定也可用井口控制參數(shù)每個站點的地表速度必須指定，可用本參數(shù)組指定也可用井口控制

47、參數(shù)組指定。若用本參數(shù)組指定，則以下表的形式：組指定。若用本參數(shù)組指定，則以下表的形式： LAYER SIN VEL INTERPOLATE LAYER 101 1400 NO參數(shù)參數(shù)1 1、層號（、層號（LAYERLAYER） 定義速度控制的層號定義速度控制的層號四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司2 2、站號（、站號（STNSTN） 速度控制點的站點號速度控制點的站點號 FIRST FIRST 第一個站點第一個站點 LAST LAST 最后一個站點最后一個站點 其它數(shù)組其它數(shù)組 ( (無缺省無缺省) )3 3、該點的層速度（、該點的層速度（VELVEL）4 4、速度

48、內(nèi)插（、速度內(nèi)插（INTERPOLATEINTERPOLATE） 從該點到下一個指定站點之間要不要速度內(nèi)插從該點到下一個指定站點之間要不要速度內(nèi)插 NO NO 不內(nèi)插不內(nèi)插 YES YES 要內(nèi)插要內(nèi)插 缺?。喝笔。篘ONO四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 六六 時間時間-深度控制參數(shù)組深度控制參數(shù)組(TIME_DEPTHS） 用戶提供時間用戶提供時間- -深度。深度。若指定本參數(shù)組，則以下表的形式：若指定本參數(shù)組，則以下表的形式： LAYER SIN TIME_DEPTHL INTERPOLATE LAYER 101 301 NO參數(shù)參數(shù)1 1、層號（、層號（LA

49、YERLAYER） 指定時間指定時間- -深度控制的層號深度控制的層號2 2、站號（、站號（STNSTN） 時間時間- -深度控制點的站點號深度控制點的站點號 FIRST FIRST 第一個站點第一個站點 LAST LAST 最后一個站點最后一個站點 其它數(shù)組其它數(shù)組 ( (無缺省無缺省) ) 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司3 3、該點的時間、該點的時間- -深度（深度（TIME_DEPTHTIME_DEPTH）4 4、時間、時間- -深度內(nèi)插（深度內(nèi)插（INTERPOLATEINTERPOLATE） 從該點到下一個指定站點之間要不要時間深度內(nèi)插從該點到下一個指定

50、站點之間要不要時間深度內(nèi)插 NO NO 不內(nèi)插不內(nèi)插 YES YES 要內(nèi)插要內(nèi)插 缺?。喝笔。篘ONO說明說明: :使用本參數(shù)組就意味著由用戶給出指定折射層的時間深度，而不希望使用本參數(shù)組就意味著由用戶給出指定折射層的時間深度，而不希望 模塊去估算。如果要這樣做就不應為該折射層提供偏移距范圍。模塊去估算。如果要這樣做就不應為該折射層提供偏移距范圍。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 七七 折射高程控制參數(shù)組折射高程控制參數(shù)組(ELEV） 用戶提供折射高程。若指定本參數(shù)組，則以下表的形式：用戶提供折射高程。若指定本參數(shù)組，則以下表的形式： LAYER ADJUSTME

51、T SIN ELEV INTERPOLATE LAYER ADJUSTMET 101 111 NO參數(shù)參數(shù)1 1、層號（、層號（LAYERLAYER） 指定折射高程控制的層號指定折射高程控制的層號2 2、更改方法（、更改方法（ADJUSTMETADJUSTMET） TIME_DEPTH TIME_DEPTH 將高程更改為時間將高程更改為時間- -深度（速度無效）深度（速度無效） VELOCITY VELOCITY 將高程更改為速度（時間將高程更改為速度（時間- -深度無效）深度無效） 缺省：缺?。?VELOCITYVELOCITY 該參數(shù)是指在繪圖時將高程值變換為時間深度或速度。該參數(shù)是指在繪

52、圖時將高程值變換為時間深度或速度。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司2 2、站號（、站號（STNSTN） 折射高程控制點的站點號折射高程控制點的站點號 FIRST FIRST 第一個站點第一個站點 LAST LAST 最后一個站點最后一個站點 其它數(shù)組其它數(shù)組 ( (無缺省無缺省) )3 3、該點的折射高程（、該點的折射高程（ELEVELEV）4 4、時間、時間- -深度內(nèi)插（深度內(nèi)插（INTERPOLATEINTERPOLATE） 從該點到下一個指定站點之間的折射高程要不要內(nèi)插從該點到下一個指定站點之間的折射高程要不要內(nèi)插 NO NO 不內(nèi)插不內(nèi)插 YES YES

53、要內(nèi)插要內(nèi)插 缺?。喝笔。篘ONO說明：本參數(shù)組通常是在做了第一輪說明：本參數(shù)組通常是在做了第一輪EGRMEGRM之后，分析人員發(fā)現(xiàn)在那些外拓的之后，分析人員發(fā)現(xiàn)在那些外拓的 區(qū)域折射模型不合常理，這時用戶就可迭代地用時間區(qū)域折射模型不合常理，這時用戶就可迭代地用時間- -深度控制參數(shù)組深度控制參數(shù)組 和折射高程控制參數(shù)組給出一個指定的外拓。和折射高程控制參數(shù)組給出一個指定的外拓。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司八八 速度平滑參數(shù)組（速度平滑參數(shù)組（VEL_SMOOTHINGVEL_SMOOTHING）例：例：LAYER SIN SMOOTH_DIAMETER IN

54、TERPOLATELAYER 100 55 NO參數(shù)參數(shù)1 1、層號（、層號（LAYERLAYER） 指定折射高程控制的層號指定折射高程控制的層號2 2、站號（、站號（STNSTN） 折射高程控制點的站點號折射高程控制點的站點號 FIRST FIRST 第一個站點第一個站點 LAST LAST 最后一個站點最后一個站點 其它數(shù)組其它數(shù)組 ( (無缺省無缺省) )3 3、平滑直徑（、平滑直徑（ SMOOTH_DIAMETERSMOOTH_DIAMETER） 該站點速度的平滑直徑，以站點個數(shù)表示。該站點速度的平滑直徑，以站點個數(shù)表示。 缺省：缺?。? 0 （不平滑）（不平滑）四川石油管理局地球物理

55、勘探公司四川石油管理局地球物理勘探公司4 4、平滑直徑內(nèi)插（、平滑直徑內(nèi)插（INTERPOLATEINTERPOLATE） 從該點到下一個指定站點之間的平滑直徑要不要內(nèi)插從該點到下一個指定站點之間的平滑直徑要不要內(nèi)插 NO NO 不內(nèi)插不內(nèi)插 （缺?。ㄈ笔。?YES YES 要內(nèi)插要內(nèi)插四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司時間時間- -深度平滑參數(shù)組（深度平滑參數(shù)組（TIME_SMOOTHINGTIME_SMOOTHING）十十 中間基準面平滑參數(shù)組（中間基準面平滑參數(shù)組（DATUM_SMOOTHINGDATUM_SMOOTHING）十一十一 折射高程平滑參數(shù)組（折射

56、高程平滑參數(shù)組（ELEV_SMOOTHINGELEV_SMOOTHING） 以上三個平滑參數(shù)組的格式與速度平滑參數(shù)組相同。以上三個平滑參數(shù)組的格式與速度平滑參數(shù)組相同。 時間時間- -深度平滑參數(shù)組平滑折射層的時間深度，由于它相當于平深度平滑參數(shù)組平滑折射層的時間深度，由于它相當于平滑所計滑所計算的靜校正量，所以要小心使用，大的平滑直徑將影響長波長和短算的靜校正量，所以要小心使用，大的平滑直徑將影響長波長和短波長靜波長靜校正量的解。校正量的解。 中間基準面平滑參數(shù)組用以產(chǎn)生中間基準面，各站點最低折射深中間基準面平滑參數(shù)組用以產(chǎn)生中間基準面，各站點最低折射深度的平度的平滑結果就是中間基準面。滑結

57、果就是中間基準面。 當折射高程出現(xiàn)明顯的異常，且這種異常很可能是速度引起的時，當折射高程出現(xiàn)明顯的異常，且這種異常很可能是速度引起的時，可用可用折射高程平滑參數(shù)組平滑所計算的折射高程。折射高程的平滑將使折射高程平滑參數(shù)組平滑所計算的折射高程。折射高程的平滑將使該折射該折射層以上的速度得到修正。（折射高程顯的異常根據(jù)質(zhì)量控制圖判斷）層以上的速度得到修正。（折射高程顯的異常根據(jù)質(zhì)量控制圖判斷）四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司十二十二 浮動中間基準面參數(shù)組（浮動中間基準面參數(shù)組（INTERMEDIATE_DATUM） 如前所述如前所述, , ,EGRMEGRM模塊根據(jù)折射

58、深度自動生成中間基準面。用戶可用本參數(shù)組修模塊根據(jù)折射深度自動生成中間基準面。用戶可用本參數(shù)組修改這個這個基準面或重新指定一個平滑中間基準面。改這個這個基準面或重新指定一個平滑中間基準面。 修改中間基準面就是將該中間基準面做上下平移，所用參數(shù)是：修改中間基準面就是將該中間基準面做上下平移，所用參數(shù)是：平滑中間基準面平移（平滑中間基準面平移（DATUM_SHIFTDATUM_SHIFT）缺省：缺?。?.00.0 重新指定平滑中間基準面則用以下參數(shù)表重新指定平滑中間基準面則用以下參數(shù)表: : STN ELEV 60 11001 1、站點號（、站點號（STNSTN） 折射高程控制點的站點號折射高程控

59、制點的站點號 FIRST FIRST 第一個站點第一個站點 LAST LAST 最后一個站點最后一個站點 其它數(shù)組其它數(shù)組 ( (無缺省無缺省) )2 2、該站點的基準面高程（、該站點的基準面高程（ELEVELEV） 無缺省無缺省四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司井口控制參數(shù)組（井口控制參數(shù)組（UPHOLE） 該參數(shù)組將井口信息用于模塊的計算中。該參數(shù)組將井口信息用于模塊的計算中。 參數(shù)參數(shù)1、井口時間應用方向（、井口時間應用方向（DIRECTION） 本參數(shù)指定井口時間應用方向是沿傾斜路徑或是垂直路徑。本參數(shù)指定井口時間應用方向是沿傾斜路徑或是垂直路徑。 RAYPA

60、TH 沿傾斜路徑沿傾斜路徑 VERTICAL 沿垂直路徑沿垂直路徑 缺?。喝笔。?垂直（井較淺時）垂直（井較淺時）BCAD井底井底四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司2 2、虛擬折射層的應用（、虛擬折射層的應用（PSEUDO_REFRACTOR_APPPSEUDO_REFRACTOR_APP） 本參數(shù)指定如何應用井口信息本參數(shù)指定如何應用井口信息 VELOCITIES VELOCITIES 井口信息用于修改用戶指定的近地表速度場井口信息用于修改用戶指定的近地表速度場 ELEVATIONS ELEVATIONS 將井口深度轉(zhuǎn)換為折射層高程（虛擬），并作為第一折射層將井口深