地震資料數(shù)字處理實(shí)驗(yàn)部分

地震資料數(shù)字處理實(shí)驗(yàn)部分二、實(shí)驗(yàn)?zāi)康?、了解地震記錄格式（SEG-Y）2、掌握Matlab中用于地震數(shù)據(jù)分析的有關(guān)函數(shù)難點(diǎn)：程序設(shè)計(jì)重點(diǎn)：3D觀測(cè)系統(tǒng)CDP面元生成和Inline、crossline道頭生成3、掌握3D觀測(cè)系統(tǒng)CDP面元生成和Inline、crossline道頭生成三、地震記錄格式地震數(shù)據(jù)文件格式,由勘探地球物理學(xué)家協(xié)會(huì)SEG(SocietyofExplorationGeophysicists)制定。(全世界統(tǒng)一的地震數(shù)據(jù)格式)1、野外(時(shí)序)記錄格式:SEGA；SEGB；SEGC；SEGD2、室內(nèi)(道序)記錄格式:SEGY(解編后的炮集記錄、CMP道集、水平疊加、疊前/疊后偏移均為此格式)SEG_Y文件結(jié)構(gòu)：三、SEG-Y格式......EndoffileEOF三、SEG-Y格式2)IBM實(shí)數(shù)格式：1)IEEE實(shí)數(shù)格式：Where:sisthesignbit=0(+)/1(-)eistheexponent=(p+bias)fisthefraction.bias=127(IEEE)or64(IBM)（1）數(shù)據(jù)格式三、SEG-Y格式（2）（text

header）40行說明信息（實(shí)例）：石東1井西區(qū)準(zhǔn)格爾盆地腹部石東1井西處理類型:偏移處理系統(tǒng):FOCUS4.3客戶：新疆石油管理局采樣間隔2ms,記錄長(zhǎng)度6000ms空間位置三維三、SEG-Y格式(3)（ReelHeader)卷頭信息

(400bytes):三、SEG-Y格式三、SEG-Y格式（4）（TraceHeader）道頭信息

(240bytes)線內(nèi)道順序號(hào)

卷內(nèi)道順序號(hào)

野外記錄號(hào)（inline）原始野外記錄中的道號(hào)炮點(diǎn)號(hào)

CDP號(hào)(xline)

CDP道集內(nèi)的道序號(hào)道識(shí)別碼l=地震數(shù)據(jù)；2=死道；3=無效道(空道）;4=爆炸信號(hào)；5=井口道；6＝掃描道；7＝計(jì)時(shí)信號(hào)；三、SEG-Y格式垂直疊加道數(shù)

水平疊加道數(shù)

數(shù)據(jù)類型：1=生產(chǎn)；2=試驗(yàn)炮點(diǎn)到接收點(diǎn)的距離(如果排列與激發(fā)前進(jìn)方向相反取負(fù)值)(分米)接收點(diǎn)的地面高程。高于海平面為正，低于海平面為負(fù)(cm)炮點(diǎn)的地面高程(cm)炮井深度(正數(shù)，cm)接收點(diǎn)基準(zhǔn)面高程(cm)炮點(diǎn)基準(zhǔn)面高程(cm)炮點(diǎn)的水深(cm)接收點(diǎn)的水深(cm)三、SEG-Y格式為（41一68字節(jié)中）高程和深度的比例因子=l，土10，土100，土1000或者土10000。如果為正，乘以因子；如果為負(fù)，則除以因子為（73—88字節(jié)中）坐標(biāo)的比例因子=1，土10，土100，土1000如果為正，乘以因子；如果為負(fù)，除以因子。三、SEG-Y格式三、SEG-Y格式三、SEG-Y格式三、SEG-Y格式每個(gè)樣值可用4字節(jié)實(shí)數(shù)表示三、SEG-Y格式（5）小結(jié){`ILINENO',189,4,`n/a',`In-linenumber'},...{`XLINENO',193,4,`n/a',`Cross-linenumber'}},四、S4m中地震數(shù)據(jù)分析常用函數(shù)（1）讀SEGY格式地震數(shù)據(jù)及顯示地震剖面

f1=‘f:/seismic_data/js02_L126_stack.segy’;%文件路徑及文件名%地震數(shù)據(jù)以單精度存放在sa.traces(樣點(diǎn)，道號(hào))中sa=single(read_segy_file(f1,{‘times’,100,2200},{‘traces’,1:150}));%讀地震數(shù)據(jù)至結(jié)構(gòu)體變量sa中

s_wplot(sa,{‘a(chǎn)nnotation’,‘cdp’});%顯示地震剖面s_ispectrum(sa,{'plottype','wiggle'},{'annotation','cdp'},{'frequencies',0,80});（2）Interactivespectrumcomputation(交互譜計(jì)算)（2）Interactivespectrumcomputation(交互譜計(jì)算)fsa=s_filter(sa,{‘ormsby’,5,10,20,40});%濾波s_compare(sa,fsa,{‘a(chǎn)nnotation’,‘cdp’});%比較顯示（3）帶通濾波及比較顯示（4）顯示卷頭、道頭、拉平（5）AGC、振幅譜（6）地震包函數(shù)(S4M)（6）地震包函數(shù)（6）地震包函數(shù)（6）地震包函數(shù)（6）地震包函數(shù)（6）地震包函數(shù)（6）測(cè)井包函數(shù)（6）測(cè)井包函數(shù)（6）地震包函數(shù)（7）讀測(cè)井?dāng)?shù)據(jù)（Las文件）、顯示測(cè)井曲線1、了解TAUP變換的原理和方法2、掌握Matlab中用于TAUP變換的有關(guān)函數(shù)難點(diǎn)：程序設(shè)計(jì)TAUP變換3、TAUP變換前后地震數(shù)據(jù)的對(duì)比分析實(shí)驗(yàn)?zāi)康腡AUP變換上圖中，當(dāng)波到達(dá)A點(diǎn)時(shí)，在地面上所觀察的早已經(jīng)通過了B點(diǎn)，和地下C點(diǎn)了。很明顯，地面和地下所觀察到的速度比波本身的速度要快，如果我們知道從震源O出發(fā)的波面到達(dá)B點(diǎn)和C點(diǎn)的時(shí)間，那我們能知道從O點(diǎn)到A點(diǎn)的時(shí)間嗎？直線映射為點(diǎn)雙曲線映射為橢圓nTau=(Tau_max-Tau_min)/DTau+1;%%%%%%%nP=(P_max-P_min)/DP+1;%%%%%%%TP(1:nTau,1:nP)=0;forindexP=1:nPP=P_min+DP*(indexP-1);forindexTau=1:nTauTau=Tau_min+DTau*(indexTau-1);forindextrace=1:trace_numt=Tau+P*(indextrace-1)*Dx;if(t<=interest_tmax&&t>=interest_tmin)indext=fix(t/Dt)+1;TP(indexTau,indexP)=TP(indexTau,indexP)+Rec(indext,indextrace);

時(shí)間域代碼實(shí)現(xiàn)示意工業(yè)實(shí)現(xiàn)一般用頻率域，效率高，但實(shí)現(xiàn)較復(fù)雜輸入地震數(shù)據(jù)

TAUP變換后一、實(shí)驗(yàn)?zāi)康淖寣W(xué)生熟悉掌握地震數(shù)據(jù)處理中的速度分析、動(dòng)校正、剩余靜校正疊加的方法的使用速度掃描實(shí)驗(yàn)二、實(shí)驗(yàn)原理地震數(shù)據(jù)處理流程圖三、主要儀器設(shè)備SUN計(jì)算機(jī)、Linux操作系統(tǒng)、地震處理軟件系統(tǒng)四、實(shí)驗(yàn)內(nèi)容與要求實(shí)驗(yàn)內(nèi)容:1、地震數(shù)據(jù)的速度分析。2、CMP道集動(dòng)校正的應(yīng)用。3、剩余靜校正量的計(jì)算與應(yīng)用。4、水平疊加處理實(shí)驗(yàn)要求:

運(yùn)用地震數(shù)據(jù)處理軟件對(duì)所給地震數(shù)據(jù)進(jìn)行抽道、速度分析、動(dòng)校正、剩余靜校正以及疊加處理。實(shí)驗(yàn)完成后，要求學(xué)生嚴(yán)格按照實(shí)驗(yàn)報(bào)告有關(guān)編寫格式進(jìn)行實(shí)驗(yàn)報(bào)告編寫，提供紙質(zhì)和電子版報(bào)告。五、處理模塊及參數(shù)描述

經(jīng)去噪和反褶積處理后的二維地震數(shù)據(jù)：測(cè)線一條，炮數(shù)99炮，每炮120道，中間放炮。

速度分析：為動(dòng)校正、偏移、時(shí)深轉(zhuǎn)換等處理提供速度參數(shù)；為巖性分析、物性分析、油氣預(yù)測(cè)等提供速度信息；為構(gòu)造解釋提供速度參數(shù)。動(dòng)校正：反射波時(shí)距曲線上非零炮檢距上的波至?xí)r間與零炮檢距道上的時(shí)間之差稱為改道的正常時(shí)差，為校正正常時(shí)差所用的速度稱為動(dòng)校正速度。動(dòng)校正是用來消除炮檢距對(duì)反射波旅行時(shí)的影響。剩余靜校正：野外靜校正后，在地震數(shù)據(jù)中仍然殘留有各種剩余靜態(tài)時(shí)移，通常在疊前必需估計(jì)出這類剩余靜態(tài)時(shí)移值，并在CMP道集中加以校正，這種校正稱為剩余靜校正。水平疊加：水平疊加是將CMP道集記錄經(jīng)NMO校正后疊加起來，目的是壓制隨機(jī)噪音，提高地震信噪比。六、處理效果分析1、地震數(shù)據(jù)的速度分析2、CMP道集動(dòng)校正的應(yīng)用3、剩余靜校正

4、水平疊加處理總結(jié)

通過本次地震數(shù)據(jù)處理課程設(shè)計(jì)，對(duì)地震數(shù)據(jù)處理的流程，包括速度分析、動(dòng)校正、剩余靜校正及疊加的方法的使用有了更深一步的