地震數(shù)據(jù)處理vista軟件使用手冊

--#-超aJj00500i540?亠0051000-1000-1500H2000超aJj00500i540?亠0051000-1000-1500H2000H2500H?|Trace:239Time:2068.00,2069.66Ampl:-378.881不知該模型數(shù)據(jù)已作過濾波處理還是模擬時的原因，數(shù)據(jù)反射波能量很強，深部的反射二、Marmousi模型模型基本信息：單邊放炮，96道接收，共240炮，道間距25m,炮間距25m，偏移距200m（8個道間距）。注：該模型炮點位于檢波器排列之后。采樣率4000微秒，每道采樣750個點?；镜奶幚矸椒ㄅc上面處理模型的方法是差不多的。2.1數(shù)據(jù)的輸入數(shù)據(jù)的輸入與上面處理模型時是一樣的，就不再重復了，下圖是加入Marmousi模型數(shù)據(jù)：SeismicData:[1]—[1]:KarKousi細魚呱幽犧膨勺筮劉潤麴mainseismicSHOTCHAN也很強，反射信息很豐富，信噪比較高，不需作其他的處理就可以做動校正和水平疊加了。之前仍需給數(shù)據(jù)加觀測系統(tǒng)。2.2建立觀測系統(tǒng)建立觀測系統(tǒng)的方法與前面也是一樣的，也不再重復了，僅就個別參數(shù)作些說明。下圖是炮檢關系和炮點坐標：

口回區(qū)I口回區(qū)I員區(qū)4BZfE?y矽圭至醬園團±lTfd冒兀ShotSupplementalSpread-SheetRECORD#FIELD#SHOT#FIRSTCHAN#LASTCHAN#BEFOREGAPFROMBEFOREGAPTOAFTERGAPFROMAFTERGAPTOSHOTBULKSHOTDEPTHSHOTUPHOLESHOELE▲DEF.INCR11.00001.001.000.000.000.0000.0000.000-I?216[1]216□319.00196216.00311.000.000.000.0000.0000.000217[1]217□320.00196217.00312.000.000.000.0000.0000.000218[1]218口321.00196218.00313.000.000.000.0000.0000.000219[1]219□322.00196219.00314.000.000.000.0000.0000.000220[1]220口323.00196220.00315.000.000.000.0000.0000.000221[1]221□324.00196221.00316.000.000.000.0000.0000.000222[1]222口325.00196222.00317.000.000.000.0000.0000.000223[1]223□326.00196223.00318.000.000.000.0000.0000.000224[1]224□327.00196224.00319.000.000.000.0000.0000.000225[1]225□328.00196225.00320.000.000.000.0000.0000.000226[1]226□329.00196226.00321.000.000.000.0000.0000.000227[1]227口330.00196227.00322.000.000.000.0000.0000.000228[1]228□331.00196228.00323.000.000.000.0000.0000.000229[1]229口332.00196229.00324.000.000.000.0000.0000.000230[1]230□333.00196230.00325.000.000.000.0000.0000.000231[1]231口334.00196231.00326.000.000.000.0000.0000.000232[1]232□335.00196232.00327.000.000.000.0000.0000.000233[1]233□336.00196233.00328.000.000.000.0000.0000.000234[1]234□337.00196234.00329.000.000.000.0000.0000.000235[1]235□338.00196235.00330.000.000.000.0000.0000.000236[1]236口339.00196236.00331.000.000.000.0000.0000.000237[1]237□340.00196237.00332.000.000.000.0000.0000.000238[1]238口341.00196238.00333.000.000.000.0000.0000.000239[1]239□342.00196239.00334.000.000.000.0000.0000.000240[1]240□343.00196240.00335.000.000.000.0000.0000.000?AddfAppend)ManyRows炮點是位于檢波器排列之后的，即第1個檢波器第96個檢波器一炮點因此和上面二層水平介質模型的炮檢關系有一點小的差別，只要正確認請炮檢的位置關系和準確確定炮檢坐標，其實是一樣的。同樣點擊關和吃計算CMP和檢波器的覆蓋次數(shù)，點擊破顯示,見下圖:下圖是檢波器坐標：

同樣點擊關和吃計算CMP和檢波器的覆蓋次數(shù)，點擊破顯示,見下圖:2-DGeoMetry:°[1]:Manousi100200CMP300400600m￡關儼吃啟前知圏蚩緒S2-DGeoMetry:°[1]:Manousi100200CMP300400600m￡關儼吃啟前知圏蚩緒S晟恩-診闘無Sub-SurfaceFoldCMP:7Fold:40[0]Recv:0.00CMP:0"45403530252015100—|0丄Shot:0.00CMP最大覆蓋次數(shù)：N=96*l/(2*1)=48次，因此上面所建的觀測系統(tǒng)是正確的。按住保存按鈕園不放，在出現(xiàn)的按鈕組合中園翳.帑帑逍囪中選擇翳將相關的信息寫入道頭，并還可以選擇瞎將炮檢關系及他們的坐標寫入文本文件，以便以后查看。2.3速度分析選取幾個CMP道集進行速度分析，可以選取CMP200，300，400，500，600這6個CMP道集進行分析，獲得這幾個CMP道集的速度譜，常速度疊加圖和一個最佳炮檢距道集,以便拾取這幾個CMP的速度。工作流如下：VelZoncOutputOutput.Output

VelZoncOutputOutput.Output雙擊VelZone圖標，填入?yún)?shù)，如下圖:右鍵點擊VelZone圖標一DataInputControl—Data/headerSelection并選擇CMP道集，如下圖所示：并設置常速疊加和速度譜掃描范圍，設為100，5000，100，設置完畢后，點擊執(zhí)行工作流。執(zhí)行完畢即可得到這幾個CMP道集的速度譜，常速度疊加圖和一個最佳炮檢距道集。選擇Interactive/VelocityTools/InteractiveVelocityAnalysi進行速度分析，在彈出的對話框中選擇速度譜，一個最佳炮檢距道集和常速度疊加數(shù)據(jù)集，如下圖：點擊OK后即可對各CMP道集進行速度分析，拾取各CMP速度，見下圖:

2-DTelocityPick:ovngy、坦基|蠢1罵創(chuàng)魚|創(chuàng)丄|Velocity:1814[Time:1068.02-DTelocityPick:ovngy、坦基|蠢1罵創(chuàng)魚|創(chuàng)丄|Velocity:1814[Time:1068.0口回岡|Semblance:0.791|ZOOMMODE可以通過左右移動*-J'…來選擇不同的CMP道集,拾取不同的CMP的速度。拾取完后，點擊保存按鈕保存為速度文件，為后面的動校正提供速度信息。因模型中部結構比較復雜，反射信息比較復雜，很難準確拾取到反射層的速度,造成速度不準確或不充分，對后面的處理有一定的影響。2.4動校正

得到速度信息后，可以對按CMP排列的地震數(shù)據(jù)進行動校正。工作流如下:動校正后淺層且大炮檢距的地震道會發(fā)生動校正拉伸現(xiàn)象，導致淺層的同相軸遭到嚴重損害，為處理方便，根據(jù)動校正拉伸后頻率下降，對動校正后的數(shù)據(jù)設置一個帶通濾波器，濾掉因拉伸而導致的畸變。OrmsbylOutputOrmsbylOutput雙擊Input圖標，設置輸入數(shù)據(jù)及其排列方式，如下所示：雙擊Nmo命令圖標，設置動校正速度文件選擇為速度分析中得到的速度文件，可以點擊按鈕WEW查看速度信息，見下圖：羈⑧T|Time:2446.26Vel:2285.2分芳y羈⑧T|Time:2446.26Vel:2285.2分芳y園[slfETD99001112223344455566777681357913579135791357917665432210988765443210033333333332222222222222CDP:294點擊雹查看各CMP的速度信息，見下圖:口回區(qū)50010001500200015002000250030003500Vnns:IntervalVelocityPreviousVelocityZZE?CDP:200|Time:2098.1MsVel:3758帶通濾波器設置如下:設置完畢后點擊執(zhí)行工作流，執(zhí)行完畢后得到動校正后的結果。見下圖:SHOTCHAN5]n>o2108106105152150879295I懾冊MI必Ml■則圳PP則口回區(qū)50010001500200015002000250030003500Vnns:IntervalVelocityPreviousVelocityZZE?CDP:200|Time:2098.1MsVel:3758帶通濾波器設置如下:設置完畢后點擊執(zhí)行工作流，執(zhí)行完畢后得到動校正后的結果。見下圖:SHOTCHAN5]n>o2108106105152150879295I懾冊MI必Ml■則圳PP則爛IM刪0!般0劃1111掘II規(guī)1111GGG-15002000J講顧8^Ampl:-26.954-1500?snn-500|Trace:Time:2668.00/2669.66VelocityDisplayFile:D:\zxB\lar>ousiTest-files\OutputMaraousiTest.vel云?]Tj云?]Tj?|4IVelocityMS2.5水平疊加道，以加強有效信號，消除干擾，提水平疊加將動校正后同一個CMP的各道疊加為道，以加強有效信號，消除干擾，提高信噪比。水平疊加工作流見下圖:Output-Output-El雙擊各命令圖標并設置完參數(shù)后，執(zhí)行工作流即可得到水平疊加剖面，見下圖:雙擊各命令圖標并設置完參數(shù)后，執(zhí)行工作流即可得到水平疊加剖面，見下圖:利用已給的Marmousi模型速度（已知的是z-層速度，應將其轉化為t-Rms）文件，對Marmousi模型做動校正并水平疊加后也得到一個水平疊加剖面，但是由于CMP號對應不上，效果也沒有改善多少，也有可能是方法本身的問題。水平疊加剖面見下圖：SeismicData:[8]-stack2阻妙％謚F團勺帑飾碰劇MAINseismicZ5I創(chuàng)購鰹?121勿旦圍CMP2033034035036031000HIlhUhlli9\[……SX.!：I!B?Z5I創(chuàng)購鰹?121勿旦圍CMP2033034035036031000HIlhUhlli9\[……SX.!：I!B?刪卿聊腳酬熬刪!■儷1觀腳嘶卩\i\AmI：!i?\H蝴呀（“幌務曲卿囁冊蹦昭吧mi（：iI伽刪申屹腫比》湍箱柯般耶即is熾!%腳翩刪熾估;:時羽專）卡犧吧幣:曲觸熾絲鷲輜曲哪剛鹹冊酬怖禰制僦擁州瓣悅陽徹獺沁￡□腳:於懸烈濫般砌勰熾怨刪桝慌總m三斛卿巾怡刖?暢砂覘圳偏腳幀吹豐叫”m代is（咄妙m側際熾唸齡呼營：：：：*；；直①:“』濾濫琛”ISI!怵制牖怖U；慨般瀚;;黑薦士R劇―魔澆空滋弦汕屈刪3it林貓曲跟S裁忻觀：；爲:;；⑴:：曲摘絢淤I瞬癮斟腳啊P艸紂桔釜，“，，/仔叫繪::沁:二I即酬C：篇泅）減懾私陽啊熾陸Ii”淵'?"?乂瓦說小滸叱火Wqm=~:..爍*：:;:/占冊".碎加.腫」、;?:."?）ii.iin冊加血腫.mi（i仞""'"：：?。?同;”?.?；■.：'：.:.、：》?”，”?.沁”:?:護曲瀘訊蕊爐：沽起沁耕私臉?縱靛離翔空洲慚曲帶縊赧眈姍炯黜怖晰倔般鵲&戒淞撕船吧褊觴鄴《洱鍛劇滋詼；怎臨贏跡我誦亦汨軸闕3卩側少；代滬&.:聲獄:泌涉擁訓4丹淵I陳陶瞞備‘網(wǎng)i僦卩加麗削曲型師％必測渝I；沁辭I妙朗世::屈加）3腳牴觸I池汕加申?啊陽仙賤售.、、?；"：.、『?.；?“.?嘟空莎帀"ENFO觀愴臨必綁贈刪隔葩瀾臨:朋鬆;懈供磁唸麵蒯翩輒彌羅也砌矽碎腸懈擁勵夠緞策:;卅匯、?毗；汎:.4藝E?畑J“:尸門沖』.*認”：'黑￡；；：（卅；：?、,：'：.“叫川；，：.“淞'""；山：少"丁、?**.嚴闇;小；.“、臥：心刪粘昭騒H比刃曲p湘則廿加世;期;鵬蹶覽幣涮兩訃吻剁編測克溢兩購I1500H2000H他MS]2500H|:||||||||Trace:557Time:2544.00/2543.12|:||||||||Trace:557Time:2544.00/2543.12Ampl:4.917+卅?川,：；：叫z*.:,.>r>-沛?箱詒煤淋氓inmm11inmumumiiiiiiiiniiiiiiiiiiiiiiiiii11對數(shù)據(jù)做反褶積壓縮地震子波后再做動校正和水平疊加，疊加剖面的局部有所改善，見下圖:SeismicData:[11]=stack5[■SeismicData:[11]=stack5[■因為Marmousi模型結構比較復雜，反射信息很復雜，僅僅通過動校正和水平疊加得到各CMP自激自收的信號，這是很難得到地下反射界面的信息的，成像效果也不可能很好，這是模型本身的結構和方法本身所決定的。

.建工程(project)File—newproject［移動盤中處理慎重］選取適當?shù)奶幚砦募顚懴卤韴D1.1二.把新的數(shù)據(jù)文件加到工程上1.新建一個地震數(shù)據(jù)到上述工程文件中(newseismicdatatoaddtoproject)N—2D/3D圖2.1把需要處理的數(shù)據(jù)(如sgy等)加到1步所建的文件中(addseismicdatatoselecteddataset)“+”fsgy/sg2出現(xiàn)圖2.2在右上框中選擇所需文件到下框中，然后一路確定，最后出現(xiàn)圖2.3，主意圖2.3中的DATA#1國薩題［1］：NEW2其中第一個按鈕顯地震記錄圖，第二個按鈕顯到頭編輯，第三個按鈕顯到頭信息圖2.2圖2.3三.建立觀測系統(tǒng)1.Interactive—geometrywindowdisplay圖3.1出現(xiàn)圖3.2圖3.3選中寫著“[l]:new2-DDATA”的按鈕，再點“0K”，既是建立觀測系統(tǒng)2.出現(xiàn)圖3.4圖3.4點“LOAD”出現(xiàn)圖3.5圖3.5加載外部做好的觀測系統(tǒng)

ReadGeometrySpread-Sheet—Readfromgeometryfilein冋rnfWWN1!區(qū)帀眾1（第二個）圖3.6GeometryParameters—CMPBindegaults（圖3.6中第一個）出現(xiàn)圖3.7點“AUTO-CALCULATE”自動計算出有關參數(shù)后再點“ok”圖3.7計算，即GeometryParameters—CalculateFold/Offset（圖3.6中第一個）把觀測系統(tǒng)加入sgy文件，WriteGeometrySpread-Sheet—WriteBinningInfoToHeaders（圖3.6第三個）一些按鈕的功能如圖3.6抽道集JobFlow—openflowfile+commandwindow在右框（vistaflowcommandwindow）中分別選選“input”和“output”到左框（project）中。雙擊“input”，選擇參數(shù)如圖4.1所示，點擊“ok”

圖4.1右鍵“input”點擊markforexecution（此時左下方顏色由黃變綠）雙擊“output”，給輸出文件名取名為“bbb”，設置參數(shù)如圖4.2所示，點擊“ok”圖4.25.如法炮制，右鍵“output”點擊markforexecution（此時左下方顏色由黃變綠）

6.將輸入輸出兩過程連接起來，點擊CreateFlowCommandLinks(如圖4.3所示第七個按鈕)圖4.3此時先點擊“input”拖動鼠標到“output”，可將二者連起來，即過程是由前頭前到箭頭后,如圖4.4：圖4.4點擊圖圖4.4點擊圖4.5所示按鈕即取消連接狀態(tài)UseMousetoDefineFlowConnectionLinks圖4.5執(zhí)行。即點擊圖4.3最后一個按鈕(ExecuteCurrentFlow)保存圖4.3第二個按鈕的第二項(writegeometryspread-sheet—writebinninginfotoheaders)帶通濾波JobFlow—openflowfile+commandwindow普通帶通濾波設立流程如圖5.1所示|工mput~|Output.El|工mput~|Output.El圖5.13.時變帶通濾波(可設定多個時窗)可將圖5.1中間的“Ormsbyh”替換為圖5.2中的“TVOrmsh”O(jiān)rmsbyITVOrnsl]OrmsbyITVOrnsl]NotchBVJorthFILTERING、'上,11■■■■圖5.24.其余的濾波方法可如法炮制，如陷頻濾波(Notch)、巴特沃斯濾波器(BWorth)、(FILTPan)、(Shape)FK濾波其中1-6為生成濾波文件顯示所要FK的地震記錄圖

圖6.1點擊上圖鼠標所指按鈕，即如圖6.2最后一個按鈕，出現(xiàn)圖6.3□PH兩岡耶阿細圜圖6.2~^\|W]0[H1W]W]I||f-KSpectrumAnalysisWindow圖6.3選擇右面下拉菜單如圖6.3所示，然后選擇要做FK的記錄（其中第八個按鈕為以方框選擇，第九個按鈕為以多半形選擇），然后點第十個按鈕“+”開始生成FK參數(shù)文件，選擇適當?shù)念l率范圍，如圖6.4所示圖6.44.濾波時各按鈕的作用囁國副冋頤同E園図聞酹網(wǎng)圜取1劃匠I圖6.5第五個：當選擇點在-0.0左邊，則濾去選擇點左邊所有的波；當選擇點在-0.0右邊，則濾去選擇點右邊所有的波。第六個：選擇垂直帶狀區(qū)域，左鍵選定一邊，鼠標拖動到另一邊。第七個：選擇水平帶狀區(qū)域，左鍵選定一邊，鼠標拖動到另一邊。第八個：切餅濾波，左鍵選定扇形一邊，鼠標拖動到扇形另一邊。第九個：選擇多邊形區(qū)域，選擇最后一個點后雙擊。第十個：刪除所有選定的區(qū)域。第十一個：部分刪除所選區(qū)域(雙擊要刪除區(qū)即可)第十二個：要做FK的輸入地震記錄第十三個：要做FK的輸入地震記錄的頻譜第十四個:執(zhí)行FK濾波按鈕第十五個：FK后的地震記錄第十六個：FK后的地震記錄的頻譜第十七個：FK前后地震記錄的比較第十八個：把FK寫入讀入的數(shù)據(jù)使用適當?shù)陌粹o，選擇要濾去波的部分，如圖6.6：然后單擊第十四個按鈕(FKfilter)進行FK濾波，再單擊第十八個按鈕(applyFKfiltertoinputdata)，濾波前后對比如圖6.7所示；DefineBoxAreaforFilterArea,SingleClickfAndRelease]toDefine2684.(1266J9849.(48431.：47013.!45595.12760.11342.：924.544177.J8506.77088.9671.21253.4835.6M17.80.0Wave:-0.14eque

n5-0.4DefineBoxAreaforFilterArea,SingleClickfAndRelease]toDefine2684.(1266J9849.(48431.：47013.!45595.12760.11342.：924.544177.J8506.77088.9671.21253.4835.6M17.80.0Wave:-0.14eque

n5-0.4-0.3-0.2WAVENUMBER-0.1-0.00.1圖6.6O區(qū)畫畫QE]司QW]岡刪刪圜團團劃|fe|Input/OutputSeismicDataTRACE[AfterFK]10152025303540455(TRACE[BeforeFK]10152025303540455(TRACE[AfterFK]10152025303540455(TRACE[BeforeFK]10152025303540455(TIMEi300ood.Av*-一星PAWf"套1u1;j<*4^-.=』*」+^..rJ*a-3>￡dr-*n1T-JrflH.it"nrfTrace：14Time:463.0Amp:-0.96圖6.7保存濾波文件。單擊圖6.5第三個按鈕，在工作空間中存一個后綴*.fkl為的濾波文件。七?頻譜分析顯示所要FK的地震記錄圖點擊上圖鼠標所指按鈕，即如圖6.2最后一個按鈕，出現(xiàn)圖6.3。把6.3圖右方的下拉菜單改為如圖7.1所示[5~||百[H|||CopyAnalysisWindow

圖7.1用倒數(shù)第二、三個按鈕選擇要頻譜分析的記錄，單擊最后一個按鈕“+”，如圖7.2H\7W\回團副HmainseismictoolbarSHOTCHANSHOTCHAN200TM300EIM400s)5006000.6-0.2.-0.2.inTrace：14inTrace：140.714286[0.000000]圖7.24.單擊工具欄中最后一個按鈕(seismicanalysiswindow)5?把圖7.1的右邊下拉菜單改為圖7.3所示[W]II|FrequencyAnalysisWindow圖7.36.單擊圖7.3最后一個按鈕“+”，在彈出的如圖7.4所示的下拉菜單中選擇適當?shù)膮?shù)后點擊“ok”，則出現(xiàn)頻譜分析圖如圖7.5所示圖7.4?ITrace：10600-Amplitude200^'/0.8-斗眄600-Amplitude200^'/0.8-斗眄0.4-nn-〒丨?丿」?Freq：26.84HzDegrees：-128.701[15.020]PHASE圖7.5八.動校正1?選取動校正命令框如圖8.1Output-圖8.12.Input參數(shù)設置如圖8.2所示，其中gp_cmp是加觀測系統(tǒng)抽過道集、此處要進行動校正的文件。圖8.23.Nmo參數(shù)設置如圖8.3圖8.3其中在file處選則速度文件，如本例中速度文件在D:\vista練習\vel.vel處StretchMute為動校正切除量5.output參數(shù)如圖8.4設置，gp_nmo是動校正后輸出文件圖8.4九.疊加1.選取命令框如圖9.1所示圖9.12.Input參數(shù)設置如圖9.2所示，其中gp_nmo動校正過的、此處要疊加的文件。

圖9.23.cmpstk參數(shù)設置如圖9.3圖9.35.output參數(shù)如圖9.4設置，gp_stack是疊加后輸出文件圖9.4十.速度分析1.建立流程如圖10.1所示，上面是相似（疊加）速度掃描，下面是常速掃描左面是抽共偏移距道集的輸出圖10.12.Velzone參數(shù)設置如圖10.2所示；其中“inputdatelist”是輸入的、抽過道集的、此處要進行速度分析的文件；“#binsfromcenterbin”后的數(shù)字表從要得出的同向軸兩邊個多出現(xiàn)同向軸的數(shù)量。

圖10.23.Offsrts參數(shù)設置如圖10.3（本例默認）圖10.34.Sembl參數(shù)設置如圖10.4,其中startvelocity為掃描起始速度，endvelocity為掃描終止速度，velocityincr為速度間隔

圖10.45.Output參數(shù)設置如圖10.5所示，其中outputname為相似掃描后速度文件名圖10.56.CVS參數(shù)設置如圖10.6所示，其中velocityfunctionlist中的數(shù)字分別為起始掃描速度，終止掃描速度和掃描速度間隔，這些數(shù)據(jù)可用“NEW”按鈕進行修改圖10.67.（下）output參數(shù)設置如圖10.7，