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

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

口回區(qū)I口回區(qū)I員區(qū)4BZfE?y矽圭至醬園團(tuán)±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個檢波器一炮點因此和上面二層水平介質(zhì)模型的炮檢關(guān)系有一點小的差別，只要正確認(rèn)請炮檢的位置關(guān)系和準(zhǔn)確確定炮檢坐標(biāo)，其實是一樣的。同樣點擊關(guān)和吃計算CMP和檢波器的覆蓋次數(shù)，點擊破顯示,見下圖:下圖是檢波器坐標(biāo)：

同樣點擊關(guān)和吃計算CMP和檢波器的覆蓋次數(shù)，點擊破顯示,見下圖:2-DGeoMetry:°[1]:Manousi100200CMP300400600m￡關(guān)儼吃啟前知圏蚩緒S2-DGeoMetry:°[1]:Manousi100200CMP300400600m￡關(guān)儼吃啟前知圏蚩緒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)的按鈕組合中園翳.帑帑逍囪中選擇翳將相關(guān)的信息寫入道頭，并還可以選擇瞎將炮檢關(guān)系及他們的坐標(biāo)寫入文本文件，以便以后查看。2.3速度分析選取幾個CMP道集進(jìn)行速度分析，可以選取CMP200，300，400，500，600這6個CMP道集進(jìn)行分析，獲得這幾個CMP道集的速度譜，常速度疊加圖和一個最佳炮檢距道集,以便拾取這幾個CMP的速度。工作流如下：VelZoncOutputOutput.Output

VelZoncOutputOutput.Output雙擊VelZone圖標(biāo)，填入?yún)?shù)，如下圖:右鍵點擊VelZone圖標(biāo)一DataInputControl—Data/headerSelection并選擇CMP道集，如下圖所示：并設(shè)置常速疊加和速度譜掃描范圍，設(shè)為100，5000，100，設(shè)置完畢后，點擊執(zhí)行工作流。執(zhí)行完畢即可得到這幾個CMP道集的速度譜，常速度疊加圖和一個最佳炮檢距道集。選擇Interactive/VelocityTools/InteractiveVelocityAnalysi進(jìn)行速度分析，在彈出的對話框中選擇速度譜，一個最佳炮檢距道集和常速度疊加數(shù)據(jù)集，如下圖：點擊OK后即可對各CMP道集進(jìn)行速度分析，拾取各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的速度。拾取完后，點擊保存按鈕保存為速度文件，為后面的動校正提供速度信息。因模型中部結(jié)構(gòu)比較復(fù)雜，反射信息比較復(fù)雜，很難準(zhǔn)確拾取到反射層的速度,造成速度不準(zhǔn)確或不充分，對后面的處理有一定的影響。2.4動校正

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

.建工程(project)File—newproject［移動盤中處理慎重］選取適當(dāng)?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”自動計算出有關(guān)參數(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è)置參數(shù)如圖4.2所示，點擊“ok”圖4.25.如法炮制，右鍵“output”點擊markforexecution（此時左下方顏色由黃變綠）

6.將輸入輸出兩過程連接起來，點擊CreateFlowCommandLinks(如圖4.3所示第七個按鈕)圖4.3此時先點擊“input”拖動鼠標(biāo)到“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普通帶通濾波設(shè)立流程如圖5.1所示|工mput~|Output.El|工mput~|Output.El圖5.13.時變帶通濾波(可設(shè)定多個時窗)可將圖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點擊上圖鼠標(biāo)所指按鈕，即如圖6.2最后一個按鈕，出現(xiàn)圖6.3□PH兩岡耶阿細(xì)圜圖6.2~^\|W]0[H1W]W]I||f-KSpectrumAnalysisWindow圖6.3選擇右面下拉菜單如圖6.3所示，然后選擇要做FK的記錄（其中第八個按鈕為以方框選擇，第九個按鈕為以多半形選擇），然后點第十個按鈕“+”開始生成FK參數(shù)文件，選擇適當(dāng)?shù)念l率范圍，如圖6.4所示圖6.44.濾波時各按鈕的作用囁國副冋頤同E園図聞酹網(wǎng)圜取1劃匠I圖6.5第五個：當(dāng)選擇點在-0.0左邊，則濾去選擇點左邊所有的波；當(dāng)選擇點在-0.0右邊，則濾去選擇點右邊所有的波。第六個：選擇垂直帶狀區(qū)域，左鍵選定一邊，鼠標(biāo)拖動到另一邊。第七個：選擇水平帶狀區(qū)域，左鍵選定一邊，鼠標(biāo)拖動到另一邊。第八個：切餅濾波，左鍵選定扇形一邊，鼠標(biāo)拖動到扇形另一邊。第九個：選擇多邊形區(qū)域，選擇最后一個點后雙擊。第十個：刪除所有選定的區(qū)域。第十一個：部分刪除所選區(qū)域(雙擊要刪除區(qū)即可)第十二個：要做FK的輸入地震記錄第十三個：要做FK的輸入地震記錄的頻譜第十四個:執(zhí)行FK濾波按鈕第十五個：FK后的地震記錄第十六個：FK后的地震記錄的頻譜第十七個：FK前后地震記錄的比較第十八個：把FK寫入讀入的數(shù)據(jù)使用適當(dāng)?shù)陌粹o，選擇要濾去波的部分，如圖6.6：然后單擊第十四個按鈕(FKfilter)進(jìn)行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]岡刪刪圜團(tuán)團(tuán)劃|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的地震記錄圖點擊上圖鼠標(biāo)所指按鈕，即如圖6.2最后一個按鈕，出現(xiàn)圖6.3。把6.3圖右方的下拉菜單改為如圖7.1所示[5~||百[H|||CopyAnalysisWindow

圖7.1用倒數(shù)第二、三個按鈕選擇要頻譜分析的記錄，單擊最后一個按鈕“+”，如圖7.2H\7W\回團(tuán)副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所示的下拉菜單中選擇適當(dāng)?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ù)設(shè)置如圖8.2所示，其中g(shù)p_cmp是加觀測系統(tǒng)抽過道集、此處要進(jìn)行動校正的文件。圖8.23.Nmo參數(shù)設(shè)置如圖8.3圖8.3其中在file處選則速度文件，如本例中速度文件在D:\vista練習(xí)\vel.vel處StretchMute為動校正切除量5.output參數(shù)如圖8.4設(shè)置，gp_nmo是動校正后輸出文件圖8.4九.疊加1.選取命令框如圖9.1所示圖9.12.Input參數(shù)設(shè)置如圖9.2所示，其中g(shù)p_nmo動校正過的、此處要疊加的文件。

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

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

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