二維處理流程的學習資料課件

1、GRISYS 處理系統(tǒng)二維地震資料處理處理流程的建立 1、了解地質任務，工區(qū)所在位置的地理位置、地表地震條件、構造和地層特征，向用戶索取必要的相關資料（近地表信息）。 2、對將要處理的數(shù)據(jù)進行調查分析，通過試驗初步掌握數(shù)據(jù)的信噪比，分析了解野外施工質量及數(shù)據(jù)品質在空間上的變化情況。 3、針對數(shù)據(jù)的實際情況，選定處理系統(tǒng)中對應的處理功能，并對該功能進行試驗，最終選定某一功能選件，即針對性措施。 4、結合以往的處理歷史檔案（新資料首次處理，應結合該地區(qū)相臨資料的處理歷史檔案），征求用戶的意見，建立處理流程，同時由處理責任人審核批準，完成“處理設計”。 一個好的處理流程應具有：科學性、合理性、針對性

2、。二維處理常用模塊 觀測系統(tǒng)定義 DBAGE PLOTGE HLABEL GRISPS IDFGEDB 野外靜校正應用 DBAST HIFIST 道分選 GASORT SORTIN DATINP DATOUT 二維疊前去噪 LINOEL DEGROR HNOIBP COHNMV（2D）RELNOI RNAB2D 振幅補償 AMPANA AMPCOM SCAMPC 地表一致性反褶積 SURSAN SURDEC SCDCON 速度分析 SVEL 二維剩余靜校正 AUREST IRESTAT(交互) 反褶積 TSDCON PRDCON TVDCON MTDCON 疊加 PASTAK GASTAK 疊

3、前時間偏移 HCPSTM DMO PRADMO DIPDMO 疊后修飾 RNATTE POLFIT BLUEFI RAPFIL SMDCOM ZPDCON 偏移 WEMIGR 最終顯示及合成圖 DBAFL AMPEQU DBASI RESECT數(shù)據(jù)帶解編STEP 1GRISYS系統(tǒng)輸入輸出模塊SEGBSEGD ： DBAIFDBAISDEMULTSEGD ： SEGDINSEGY(tape+disk) ： SEGYIN、SGYOUTGRISYS(TAPE) ： TAPEIN、TAPOUTGRISYS(DISK) ： DISKIN、DSKOUT、SORTIN DSKTSO、DATINP、DATO

4、UTCGG ： CGGINP、CGGOUTWGC ： CODEIN、CODOUTSEGDIN SEGD格式磁帶解編 為配合海上數(shù)據(jù)處理，對SEGDIN模塊進行了較大規(guī)模的改進，增加海上SEGD格式數(shù)據(jù)讀取的功能，并能準確地從原始數(shù)據(jù)中讀取相關數(shù)據(jù)并存放在對應的GRISYS數(shù)據(jù)道頭字中。 GRISYS/WS-V8.0 改進內(nèi)容 增加MP參數(shù)：為乘MP因子輸入標志。 這是一個選擇性參數(shù)，填此選件后，程序從數(shù)據(jù)的頭塊中求出因子的實際值與儀器記錄數(shù)值相乘。從而，得出近似儀器增益放大前的數(shù)值。一般數(shù)據(jù)輸入時不用此參數(shù)。僅當做數(shù)值分析輸入時慎用，缺省為不應用因子。同時，取消TRUE參數(shù)。 SEGDIN G

5、RISYS/WS-V8.0 改進內(nèi)容 GRISYS/WS-V8.0 改進內(nèi)容 道頭增加內(nèi)容為： HD80=水、陸檢排列的標識。 HD81=海洋資料采集拖纜的纜號。 HD82=分放1，2，3，4排列。 HD83=接收點站號。 HD84=接收線號。 HD85=炮線號。（SEA參數(shù)起作用） HD86=炮點站號。（SEA參數(shù)起作用） HD88=磁帶卷號。 HD89=通道號。 HD90=唯一道號。 SEGDIN磁盤和磁帶數(shù)據(jù)分析TPUITAPUILISEGYADISGYTAPE 8.0新增TPU、ITAPUIL 磁帶數(shù)據(jù)分析 兩者功能基本相同，TPU為早期的版本，TAPUIL為交互版本。功能主要有： 磁

6、帶數(shù)據(jù)格式自動分析； 磁帶數(shù)據(jù)DUMP； 磁帶數(shù)據(jù)拷貝。 其中磁帶數(shù)據(jù)拷貝又提供多種拷貝方式：帶對帶、帶對盤、盤對帶、合并拷貝以及控制拷貝等等。ISEGYAD 交互SEGY數(shù)據(jù)顯示 本模塊的主要功能是在GRISYS系統(tǒng)中，增加對SEG-Y數(shù)據(jù)磁盤文件的監(jiān)控、對未知數(shù)據(jù)文件的格式分析等。 具體功能為：（1）對數(shù)據(jù)文件打開的同時，對數(shù)據(jù)文件進行格式分析；（2）對SEG-Y數(shù)據(jù)體顯示；（3）對SEG-Y數(shù)據(jù) “C卡”的顯示、修改和轉存；（4）對SEG-Y數(shù)據(jù) “400字節(jié)”的文件頭塊主要內(nèi)容的顯示；（5）對SEG-Y數(shù)據(jù) “240字節(jié)”的道頭主要內(nèi)容的顯示。交互SEGY磁帶數(shù)據(jù)檢查 ISGYTAPE

7、 在GRISYS/WS-V7.0中推出了SEGY磁盤數(shù)據(jù)檢查的交互模塊。此次版本升級我們將SEGY磁帶數(shù)據(jù)檢查補充進來，使得SEGY數(shù)據(jù)檢查功能更加完善。其具有以下幾個特點： 數(shù)據(jù)DUMP； 3200和400字節(jié)頭塊檢查； 顯示參數(shù)可調； 道頭數(shù)據(jù)即時顯示； 任意道頭數(shù)據(jù)曲線繪制。 GRISYS/WS-V8.0 新增內(nèi)容 GRISYS/WS-V8.0 新增內(nèi)容主界面及磁帶機選擇磁帶數(shù)據(jù)顯示即時道頭信息讀取交互SEGY磁帶數(shù)據(jù)檢查 GRISYS/WS-V8.0 新增內(nèi)容交互SEGY磁帶數(shù)據(jù)檢查 其他功能3200 頭塊400 頭塊數(shù)據(jù)DUMP觀測系統(tǒng)定義（專題）STEP 2道編輯STEP 3 IM

8、UTEDI 交互地震道切除與編輯 IMUTEDI 交互地震道切除與編輯(新） IDATVIEW 交互地震數(shù)據(jù)顯示 IDATCOMP 交互數(shù)據(jù)對比顯示 IPRTEST 交互參數(shù)試驗IMUTEDI主 要 功 能顯示地震道道集數(shù)據(jù)及道頭字交互定義切除交互定義整道充零和開時窗充零交互定義地震道極性反轉按道頭字18交互定義整道編輯庫分析顯示指定波形的線性速度IMUTEDI道頭字顯示IMUTEDI交互定義切除imutedi.old交互道編輯線性速度分析imutedi.old退出程序后,程序將用戶定義的所有庫以標準GRISYS編碼格式存入filename.tool_code中將切除庫的二進制數(shù)據(jù)存入muli

9、b.bin中，以備下次使用imutedi.old 隨著地震勘探技術和計算機技術的不斷發(fā)展，地震數(shù)據(jù)體也變得越來越大。為了適應不斷增長得地震數(shù)據(jù)體對處理系統(tǒng)的需求，此次版本升級對于大數(shù)據(jù)量數(shù)據(jù)的處理進行了較大的改進，并新開發(fā)了兩個針對大數(shù)據(jù)體顯示和編輯的模塊。 新版的IMUTEDI； 交互地震數(shù)據(jù)顯示模塊IDATVIEW。 GRISYS/WS-V8.0 新增內(nèi)容 GRISYS/WS-V8.0 新增內(nèi)容交互地震道切除與編輯 IMUTEDI（NEW） 新版的IMUTEDI模塊除了沿襲了原有模塊的主要功能以外，具有使用方便，支持大道集、大數(shù)據(jù)體的輸入和快速定位等特點。 索引文件引入，實現(xiàn)大數(shù)據(jù)體的快速

10、定位檢索； 支持大道集，缺省以2000道為限顯示頁面； 雙滑動桿設計，一個為數(shù)據(jù)體滑動，另一個為道集內(nèi)部滑動； 切除庫和道編輯庫直接讀入。 GRISYS/WS-V8.0 新增內(nèi)容啟動IMUTEDI數(shù)據(jù)加載建立索引文件數(shù)據(jù)操作保存道編輯參數(shù)庫具體操作流程退出 GRISYS/WS-V8.0 新增內(nèi)容IMUTEDI（新）主菜單及子菜單項 GRISYS/WS-V8.0 新增內(nèi)容IMUTEDI（新）支持多編輯庫同時操作，并不互相干擾。 GRISYS/WS-V8.0 新增內(nèi)容IMUTEDI（新）多種道編輯庫的文本文件輸入和輸出，無需二進制文件。 GRISYS/WS-V8.0 新增內(nèi)容IMUTEDI（新）其

11、他各項功能及對話框 GRISYS/WS-V8.0 新增內(nèi)容IDATVIEW 較好地適應大數(shù)據(jù)體輸入； 獨立的數(shù)據(jù)索引方式； 數(shù)據(jù)的快速準確定位； 無極放大功能； 支持雙屏獨立操作； 自動記憶10屏顯示。交互地震數(shù)據(jù)顯示啟動IDATVIEW數(shù)據(jù)加載建立索引文件數(shù)據(jù)操作退出具體操作流程 GRISYS/WS-V8.0 新增內(nèi)容主界面數(shù)據(jù)加載建立索引數(shù)據(jù)操作IDATVIEW無極放大功能IDATVIEW GRISYS/WS-V8.0 新增內(nèi)容 GRISYS/WS-V8.0 新增內(nèi)容IDATVIEW支持雙屏操作；彼此獨立。 GRISYS/WS-V8.0 新增內(nèi)容IDATVIEW數(shù)據(jù)直接定位多種顯示方式；自

12、定義顯示道頭。顯示方式調整 GRISYS/WS-V8.0 新增內(nèi)容交互地震數(shù)據(jù)顯示和對比 IDATCOMP 在GRISYS/WS-V7.0中推出了繪圖文件對比的交互模塊，該模塊在地震資料處理中監(jiān)控處理質量和處理效果起到了很好的作用。但是實際上，地震數(shù)據(jù)的直接對比更能反映地震數(shù)據(jù)真實的變化。為此我們此次推出了GRISYS地震數(shù)據(jù)的對比顯示工具。其主要功能是： 雙數(shù)據(jù)輸入，獨立或同步翻動比較； 顯示參數(shù)統(tǒng)一，并支持多種顯示方式； 鼠標點擊切換數(shù)據(jù)，快捷方便； 支持不同的GRISYS數(shù)據(jù)格式。 GRISYS/WS-V8.0 新增內(nèi)容主界面數(shù)據(jù)類型定義數(shù)據(jù)顯示及對比顯示參數(shù)定義不同類型顯示方式野外靜校

13、正折射波靜校正(專題)STEP 4&5子波整形(可選)STEP 6子波整形 GRISYS處理系統(tǒng)中，子波整形由兩個模塊來完成，一是整形濾波器設計（SHPFIL），二是 整形反褶積（ SHDECN）。SHPFIL：同一工區(qū)采集的資料，由于使用不同震源或者使用不同記錄系統(tǒng)施工，致使信號在工區(qū)內(nèi)不能很好地連續(xù)對比，使用子波整形校正可消除子波之間的差異，從而改善處理效果。 首先輸入?yún)⒖紨?shù)據(jù)道以及需要校正的數(shù)據(jù)道，指定參與計算子波的參考道和需要校正道的范圍和時窗大小，近似計算出參考道的子波和需要校正的道的子波，利用維納濾波器求出這兩種子波的差異，即濾波因子。子波整形：JOB P：DQ, L：595, T

14、：332：DISKIN DFN: GZHY0001， (參考道) DFN: GZHY0000，NTY2 (要校正的道)：SHPFIL MU02, CMP97, TZ4950, CMP100, TZ4950, LW380, ACMP96, ATZ4950, ACMP98, ATZ4950, ALW380, OUT1，SHFN:SHAB1, PSTK, SHOT, MAX4000, MIN100, AMAX4000, AMIN100, L100, TAP50, NOIS1, LIST, FOLD120, HD20, ALL, OUT0, PLOT END子波整形 把整形濾波器設計應用模塊（SHPF

15、IL）求得的整形因子，應用到數(shù)據(jù)上，達到子波一致的目的 。 ：JOB P：DQ, L：595, T：332 ：DBAMU VN01, HD4, (1-500)=M100X50, M1000X3000, TAP10 ：DISKIN DFN: GZHY0000, LTR3000, HD2=(195), NTY2, STOP #：SHDECN (-0.05, 0.015, -0.122, -0.083, -0.123, 0.0026, -0.049, -0.033, -0.072, -0.004, ），LIST ：SHDECN SHFN:SHAB1，LIST ：DSKOUT DFN：GZHY0001

16、，F(xiàn)MT4 ： END參考炮子波整形前單炮GRISYS常用的反褶積模塊參考炮頻譜子波整形前頻譜GRISYS常用的反褶積模塊參考炮子波整形后單炮GRISYS常用的反褶積模塊子波整形后頻譜參考炮頻譜GRISYS常用的反褶積模塊原始CMP道集845Hz濾波后 存在相位差子波整形校正兩種檢波器GRISYS常用的反褶積模塊匹配因子對炸藥震源激發(fā)單炮進行匹配濾波對比GRISYS常用的反褶積模塊可控震源激發(fā)的原始單炮匹配濾波后的單炮GRISYS常用的反褶積模塊子波整形前疊加剖面（INLINE）子波整形后疊加剖面（INLINE）子波整形前疊加剖面（INLINE）子波整形整形后的處理： 子波整形后的數(shù)據(jù)，應進行

17、地表一致性反褶積或多道統(tǒng)計子波反褶積處理，以此來進一步完成子波的統(tǒng)一。因為在預測步長較大時，地表一致性反褶積的主要功能是“子波調整”，壓縮子波的處理由后續(xù)的統(tǒng)計子波反褶積或單道反褶積等模塊來完成，即串聯(lián)反褶積。疊前去噪(專題)STEP 7IPRTEST主 要 功 能可進行頻譜分析、俞氏子波參數(shù)選取、視速度拾取等分析DEGROR 內(nèi)切濾波法壓制面波ATTGRO 自適應衰減面波RELNOI 疊前線性干擾濾除HNOIBP 強能量干擾的分頻壓制TSDCON 兩步法統(tǒng)計子波反褶積PRDCON 預測反褶積參數(shù)試驗DEGRORIPRTEST參數(shù)試驗IPRTESTDEGROR去面波前后單炮參數(shù)試驗HNOIBP

18、IPRTEST參數(shù)試驗HNOIBP-壓制高能干擾前后單炮IPRTEST參數(shù)試驗IPRTESTRELNOI參數(shù)試驗IPRTESTRELNOI-去線性干擾前后單炮參數(shù)試驗IPRTESTPRDCON參數(shù)試驗IPRTESTPRDCON-預測反褶積前后單炮參數(shù)試驗IPRTEST頻譜分析參數(shù)試驗IPRTEST俞氏子波參數(shù)選取參數(shù)試驗 GRISYS/WS-V8.0 改進內(nèi)容IPRTEST新增7個模塊，分別為DBAFL、PVNOEL、LINOEL、FILTER、SFNPRE、POLFIT、SCALER；新增模塊試驗方式與原有模塊基本相同；增加交互拾取擬合參數(shù)功能。振幅補償STEP 8用于分析由于 地表因素 球

19、面擴散 介質吸收 影響而引起的振幅能量衰減情況AMPANA振幅分析AMPANA: JOB P：PPP，L：L302，T：AMPANA: DBAFL VN01，I(10): TAPEIN UNIT01，G302，F(xiàn)1，HD2=(2-11)，STOP: AMPANANX60，NT10，DB，W500-W5000， FL01，(TR1-TR60，I3) : END NXn n為要處理的炮檢距個數(shù) NTm m為對每個炮檢距所規(guī)定的疊加道數(shù) DB 填此標識時，采用對數(shù)比例繪圖 振幅分析AMPANA不同炮檢距平均振幅曲線振幅分析振幅處理 AMPCOM 振幅補償 SCAMPC 地表一致性振幅補償 AMPEQ

20、U 振幅動態(tài)平衡（E選件）AMPEQU空單道平衡選件A整道平衡選件B多時窗區(qū)域平衡C均方根平衡選件D瞬時平衡選件E單炮平衡選件 振幅處理AMPEQU 炮間均衡： DISKIN DFN：GHIFIST1，NTY68（NTY2）： AMPEQU E，NTA60： DSKOUT DFN：GAMPEQU1 E：在用戶定義的單炮記錄道的規(guī)定時窗內(nèi)，求取平均振幅值，再根據(jù)提供的區(qū)域平衡級別系數(shù)計算出相應記錄的平衡系數(shù)，將其與實際振幅值相乘，求其平衡后的振幅值。 振幅處理AMPEQU炮間能量均衡前 炮間能量均衡后振幅處理AMPCOM主 要 功 能消除地震船汽槍能量不同所造成的對地震波振幅的影響 取地震道振幅

21、絕對值 消除球面擴散對振幅的影響消除地層吸收對振幅的影響將地震道振幅規(guī)格化到一給定值 地震道振幅常數(shù)加權 振幅處理AMPCOM： AMPCOM T： AMPCOM T1.5： AMPCOM T2，ABS： AMPCOM ABC，C2： AMPCOM ABC，A3，B4，C2，NR3000： AMPCOM ABS： AMPCOM SDC，TV02，WGT1000 振幅處理AMPCOM振幅統(tǒng)計分析timeAtime振幅處理AMPCOM振幅統(tǒng)計分析timeAtime振幅處理AMPCOM振幅補償前振幅處理振幅補償后AMPCOM振幅處理AMPCOM振幅補償前后單炮振幅處理SCAMPC去除由激發(fā)、接收等因

22、素引起 的接收道的能量差異,包括 幾何擴散和吸收系數(shù)補償 地表一致性振幅道平衡 振幅處理SCAMPC: JOB P：PPP，L：LLL，T：SCAMPC: SCAMPCFP1，LIST1，TI300，TA1200， LW400，NT120，T1500， V2500，MIX300，MAX3250: DSKOUT .: PLOTTRSCL75，STP15，TW50: END 只做幾何擴散和吸收系數(shù)補償 振幅處理SCAMPC: JOB P：PPP，L：LLL，T：SCAMPC: SCAMPCTI300，TA1200，LW400，TCI400，TCA1300， NT120，T1500，V2500，MI

23、X300，MAX3250， STI61，STA152，I1，SHOT87 : DSKOUT .: PLOTTRSCL75，STP15，TW50: END 做幾何擴散和吸收系數(shù)補償及地表一致性道振幅平衡 振幅處理SCAMPC地表一致性振幅補償前單炮振幅處理地表一致性振幅補償后單炮SCAMPC振幅處理一次速度分析剩余靜校正（專題）剩余靜校正見靜校正專題STEP 9 速度分析疊加速度分析 常速掃描速度分析速度分析數(shù)據(jù)解釋速度譜解釋常速掃描解釋 層位拾取 速度處理速度換算： 均方根速度換算層速度，等時速 度內(nèi)插，沿層速度內(nèi)插DBAPV插值交互速度編輯與修改常速掃描疊加的速度與現(xiàn)有速度場的結合功能介紹速

24、度分析所要求的相應處理功能NMO 疊加其它輔助性功能TV 函數(shù)加載TV 函數(shù)輸出Mute 函數(shù)加載Mute 函數(shù)輸出速度處理和速度場質控圖件顯示TV函數(shù)剖面展開顯示疊加速度二維彩色剖面顯示速度等值線圖顯示速度譜解釋窗口常速掃描解釋窗口7.0 對新增和改進部分的說明 作csum 大道集的選項,便于質量控制; 條件約束自動尋找極值點 速度微調 tv函數(shù)展開顯示 譜線 + 變面積顯示（同批量譜顯示） 常速掃描疊加的速度場可選擇修改掃描范圍內(nèi)的時間速度對 Svel 簡單操作過程在%狀態(tài)下或 gscons下啟動svel，彈出主界面按照二級管理（在主界面tool下project manage 項）加載數(shù)據(jù)

25、，同一地區(qū)的數(shù)據(jù)在load中可同時加載多條測線數(shù)據(jù)。主界面File選項，打開本次分析的數(shù)據(jù)。主界面Seisfunc功能項中，選VanalyAnalyse彈出作速度譜的參數(shù)框 。此參數(shù)框中有三個 參數(shù)框必須人工干預 : CMPs for vel. analyse Parameters Muting function definition Parameters Ref velocity curves Parameters 以上三個參數(shù)框填寫好后按OK鍵，其它參數(shù)用缺省值，在速度譜解釋窗口中顯示出速度譜、大道集和變速疊加段 。TV&MU輸出 主界面File窗口選Svel項, 彈出窗口，在Data t

26、ype欄中選擇保存類型。如: TV庫，選TVPair，在Name對應的參數(shù)框填輸出文件名后點 OK，存成內(nèi)部文件。 主界面：ToolTV I/O選項，彈出窗口，點擊OUTPUT，選右邊欄中的文件名，變黑后，在Grisys TV Filename下邊的路徑指向上，填寫輸出文件名，點OK，完成輸出功能。此時可在當前目錄下找到填寫的輸出文件。第一次迭代第二次迭代第四次迭代優(yōu)勢頻帶的速度分析與剩余靜校正的迭代注意事項 TV庫中第一個T切除庫中第一個T時間。 輸入數(shù)據(jù)應經(jīng)過振幅恢復處理的cmp集數(shù)據(jù) 譜圖太窄時，可改變Maxim Offest參數(shù)，使其小于實際偏移距。 偶而出現(xiàn)運行不正常時，一要查看數(shù)入

27、數(shù)據(jù)正確性，二建立新工區(qū)，重新填參數(shù)或進入/grisys/oper/optrun/svel目錄中，清除本工區(qū)的所有中間文件 GRISYS/WS-V8.0 改進內(nèi)容SVEL 8.0改進內(nèi)容支持1ms采樣數(shù)據(jù)；此次修改主要面向大數(shù)據(jù)體；CMP數(shù)/號最大可以支持9999999；改進譜點定義方式，增加起止點及增量定義方式；速度譜一次可做512個；對于沒有解釋的譜點，生成TV文件時，不再輸出。 GRISYS/WS-V8.0 改進內(nèi)容滑動條，CMP個數(shù)2000時自動分屏；滑動條上的數(shù)字為當前頁面的第一個CMP號；拖動滑動桿來定位數(shù)據(jù)。新的速度譜點定義方式:除了保留原有的定義方式之外,新增了起始點、終止點和

28、增量的定義方式。交互速度分析 GRISYS/WS-V8.0 改進內(nèi)容DBATV/DBAVI 速度庫定義增加速度庫的時變空變比例系數(shù)加權功能。具體是：去除原有模塊中的“PE”參數(shù)，新增“TPmPn”參數(shù)，從而實現(xiàn)速度的時空變比例加權。 GRISYS/WS-V8.0 新增內(nèi)容速度平滑編輯 DBAVE DBAVE模塊根據(jù)用戶提供的參數(shù)對速度庫進行如下處理： 速度平滑處理 速度值編輯修改 速度時變比例系數(shù)加權：JOB：DBAFL VN01：DBAMU VN01：DBATV/DBAVI VN01：DBAPV/DBAVS VN01：DISKIN：FANMOV TV01/VI01，MU01：PASTAK L

29、W300 或者 ：GASTAK： DSKOUT： FILTER FL01： AMPEQU L300：PLOTTR SCL100，STP10：END二維疊加作業(yè)地表一致性反褶積（可選）STEP 10地表一致性反褶積（第一步） 地表一致性反褶積可以消除炮點、檢波點隨地表條件和采集因素的變化對子波帶來的影響，可輸出穩(wěn)定和基本一致的地震子波 SURSAN對輸入的地震數(shù)據(jù)按時窗求取其對數(shù)功率譜，將對數(shù)功率譜按同一分量分別進行累加，對共炮點、共檢波點、共偏移距、共中心點等五個分量進行地表一致性相關分析。SURSAN參數(shù)說明 :NINTXa, Wb-Wc a為炮檢距，b、c為時窗起止時間值(ms)，因為只定

30、義一個時窗，其時窗范圍盡量大一些，以便包括較多的有效層。 由道號定義時窗 TRd, We-Wf d為道號，e、f為時窗起止時間值(ms)。 PK(g, h) g相關分析時第一分量所選用的道頭字，h為g所指定的道集在輸入數(shù)據(jù)中的道集數(shù)。 SK(i,j) i為相關分析時第二個分量所選用的道頭字，j為i所指定的道 集在輸入數(shù)據(jù)中的道集數(shù)。 通常情況下取g=2，i=18進行相關分析，所得效果較為理想。 LWl l為相關長度(ms)。 SPFN: Gxxxxxx 為存放PK分量的對數(shù)功率譜的磁盤文件 GPFN: Gxxxxxx 為存放SK分量的對數(shù)功率譜的磁盤文件MAXFRf f為相關分析的最大頻率，缺

31、省為80.0Hz。 NITERSn n為迭代次數(shù)，缺省為10。 NTm m為輸入道數(shù)，缺省為全部輸入道SURSAN使用實例 :第一步：JOB P: LNC, L: SGS, T: DEGROR ：DISKIN DFN: GDBA02, HD2=(1-300),NTY2： FILTER FL01： AMPEQU L300：MUTING BGIN,MU01 ：SURSAN X100，W400W2500， X3050，W370024500，NT30000， MAXFR120，LW200，NITERS20， PK(2, 300)，SK(18, 130)， SPFN：GLF00001，GPFN：GLF0

32、0002 ：END SURSAN地表一致性反褶積（第二步） SURDEC對輸入利用SURSAN模塊算出的地表一致性反褶積的譜分析結果對輸入的地震 ,道先進行地表一致性反褶積算子的求取，然后利用該算子進行地表一致性反褶積。 SURDEC參數(shù)說明 :SPFN: Gxxxxxx 為存放PK分量的對數(shù)功率譜的磁盤文件。 GPFN: Gxxxxxx 為存放SK分量的對數(shù)功率譜的磁盤文件。 Ll l為反褶積算子的長度，MUn n為切除庫號。 SPIKE 脈沖反褶積，為缺省方式。 GAP 預測反褶積。 BPSPK 帶通脈沖反褶積。 BPGAP 帶通預測反褶積。 TWO 填此選件時，既對pk分量方向做反褶積，

33、同時也對sk分量方向做反褶積,缺省為僅做pk分量的反褶積。GAPLNg g為預測步長，缺省為1個采樣間隔。 FLa a為帶通低截頻，缺省為5.0(Hz)。 FHb b為帶通高截頻，缺省為60.0(Hz)。 a, b的范圍應大于有效波的頻率范圍。 CFc c為白噪百分比，缺省為0.1。 CFOUTd d為加入的通放帶外的功率譜的百分比，d缺省為50。此時兼顧了信 噪比與分辨率，d變小時分辨率提高，d變大的信噪比提高SURDEC 使用實例 :第二步 預測反褶積，算子長為120，預測步長為30 （算子的長度不能大于SURSAN中的LW參數(shù)）：JOB P: LNC, L: SGS, T: SURDEC

34、：DISKIN DFN: GDBA02, HD2=(1-300),NTY2: SURDEC TWO, GAP, L120，GAPLN30，F(xiàn)L6, MU01, FH60，CF20，CFOUT40.0， SPFN：GLF00001， GPFN：GLF00002 : ENDSURDEC原始單炮地表一致性反褶積后的單炮二次速度分析剩余靜校正同一次速度分析和剩余靜校正STEP 11反褶積STEP 12反褶積：改造地震記錄子波的一種運算。理想情況為把子波改造成單位脈沖?；炯僭O： 子波是小相位的 反射系數(shù)序列為白噪反褶積的基本概念 地震記錄的褶積模型 為地層反射系數(shù)； 為地震子波。 地震子波包含震源響應

35、、檢波器響應，這些響應是時不變的；現(xiàn)在一般檢波器都是組合接收，還有組合響應（若炮組合也有組合響應），還有地層吸收，而后二者是時變的。 反褶積的目的是想消除地震子波 影響，從而得到地下反射系數(shù) ，克萊伯特形象地將反褶積比喻為顯微鏡，使我們能更細致地觀察到地下物性變化（縱向）。如有剩余子波也希望是零相位的，因為零相位子波分辨率最高。反褶積的基本概念單道反褶積： 脈沖、預測 (TVDCON、PRDCON) 前提均是子波最小相位,并假設反射系數(shù)為白噪,但是實際情況地震子波是混合相位,且頻帶有限,通過實驗得出脈沖反褶積壓縮子波功能最強,預測其次.多道反褶積： TSDCON、SURDEC、SCDCON 前

36、提是假設子波具有地表一致性，子波最小相位和反射系數(shù)為白噪。 SURDEC、SCDCON以子波調整為主，TSDCON以壓縮子波，提高分辨率為主。 反褶積的基本概念頻率分析 SPEANA 頻譜分析 IFANSFL 交互頻譜分析 FBSCAN 頻率域濾波與頻帶掃描 計算并顯示地震信號的振幅譜和相位譜：JOB P：NKK，L：LNK1，T：SPEANA ：DISKIN DFN：G0000321， HD4=(150-152) ：SPEANA HD4， (150-152， W500-W1000)， AMP2,LTWO ：END AMP2，以增益方式繪振幅譜SPEANA頻率分析SPEANA頻率分析 GRIS

37、YS/WS-V8.0 改進內(nèi)容SPEANA 頻譜分析 改進頻譜分析結果的輸出方式，使得其結果和交互頻譜分析的結果保持一致。 GRISYS/WS-V8.0 改進內(nèi)容交互頻譜分析和濾波器設計頻率分析IFSANFL 可實時的對地震記錄道形式的各種數(shù)據(jù)進行Fourier譜分析，包括單道一維譜、多道一維譜及多道的二維譜（即FK譜）； 在一維譜（包括單道和多道）分析的基礎上進行濾波器設計和試驗，濾波器包括低通、帶通、高通、帶陷、雷克子波和俞氏子波； 可進行頻帶掃描分析；可進行時頻譜分析；可進行多時窗的頻譜對比，其中包括不同數(shù)據(jù)集上相同時窗的頻譜對比，相同數(shù)據(jù)集上不同時窗的頻譜對比；可定義并應用濾波器的作用

38、范圍及強度，從而實現(xiàn)時變和空變?yōu)V波；可定義切割濾波；可對得到的濾波器建立相應的GRISYS濾波庫，形成相應批量作業(yè)，且用戶可隨時查看與編輯該作業(yè)文件。 頻率分析主界面數(shù)據(jù)屬性及顯示參數(shù)數(shù)據(jù)顯示IFSANFL頻率分析IFSANFL頻率分析選定分析范圍頻譜單道頻譜濾波器設計器IFSANFL頻率分析選定分析范圍IFSANFL設計濾波器應用效果效果分析數(shù)據(jù)應用頻率分析IFSANFL多時窗對比分析 GRISYS/WS-V8.0 改進內(nèi)容IFSANFL一維譜增加分貝顯示；一維譜加蓋網(wǎng)格，并可修改顯示參數(shù)；改變頻率掃描結果顯示方式，增大顯示面積。 GRISYS/WS-V8.0 改進內(nèi)容交互頻譜分析顯示結果和

39、方式與SPEANA模塊一致 GRISYS/WS-V8.0 改進內(nèi)容交互頻譜分析增大了結果顯示面積，方便對比觀察分頻掃描頻率分析FBSCAN 最常用也是最直觀的原始資料頻率特征分析手段。常用于現(xiàn)場處理的野外施工參數(shù)對比分析以及室內(nèi)處理中原始資料的頻率分析。 實際上就是對輸入數(shù)據(jù)進行多個不同的濾波處理，并將濾波后的數(shù)據(jù)顯示在同一張剖面上進行對比分析。通過剖面上隨著濾波門的變化而產(chǎn)生的數(shù)據(jù)信噪比的變化來了解資料中有效波和干擾波的頻率分布范圍和基本特征。 需要配合PLOTTR模塊使用。FBSCAN：JOB P：NKK，L：LNK1，T：SPEANA ：DISKIN DFN：G0000321， HD2=

40、(150-152) ,STOP：AMPEQU A: FBSCAN BP(7,10,20,28),BP(10,15,30,42), BP(14,20,40,56),BP(21,30,60,84), BP(28,40,80,112),BP(35,50,100,140), BP(42,60,120,168),BP(50,70,140,196), BP(56,80,160,224),M120: PLOTTR SCL100,STP10,GAP2=HD22,HDA2,HDB17：END 頻率分析FBSCAN10-20HZ15-30HZ20-40HZ頻率分析FBSCAN30-60HZ40-80HZ50-10

41、0HZ60-120HZ頻率分析時頻分析 TFSCAN 應用一系列具有不同中心頻率和寬度的三角形濾波器，對地震道進行濾波，將同一個地震道的濾波結果安濾波器中心頻率從小到大順序排列，可研究地震記錄中不同頻率成分在不同時間上的響應，應用于地震地層學分析，時頻補償。頻率分析TFSCAN：JOB P：NKK，L：LNK1，T：SPEANA ：DISKIN DFN：G0000321， HD2=(150-152) ,STOP：AMPEQU A*: TFSCAN：TFSCAN ENVELP : PLOTTR SCL100,STP10,GAP2=HD22,HDA2,HDB17：END ENVELP:指明輸出地震

42、道的包絡，缺省輸出地震道。HD22頻率分析地震道 時頻分析結果 時頻分析結果ENVELP頻率分析地震道 時頻分析結果頻譜分析頻率分析 本模塊用兩種典型子波俞氏子波和雷克子波對記錄進行濾波掃描，處理員可分析掃描結果，用于選擇適當?shù)臑V波參數(shù)和子波反褶積的期望輸出。 當輸入為脈沖函數(shù)時，本模塊輸出為參數(shù)對應的子波；當輸入為地震記錄時，用俞氏子波或雷克子波對記錄做濾波處理，輸出為濾波掃描的結果。頻率分析YUSACN 俞氏及雷克子波掃描YUSCAN: JOB P:2D,L:001,T:T: DATGEN SI1,LTR4000,NT4,RD,D1000,Z4100,1,ZG1000: YUSCAN YU

43、（10，60），YU（25，150），YU（35，300）: AMPEQU A: PLOTTR SCL0600，STP100，LR，HEAD，G0，NR,W0900-W1100: END: DISKIN DFN:HLABEL1,HD2=20,STOP：YUSCAN RICK（10，30，50，80，100，200）*: YUSCAN YU（10，60），YU（25，150），YU（35，300）：AMPEQU L300：PLOTTR SCL0100，STP100，LR，HEAD，GAP3HD22輸入脈沖函數(shù)輸入地震道頻率分析YUSCAN原始單炮RICK（20）RICK（30）RICK（40）頻

44、率分析YUSCANRICK（40）RICK（50）RICK（60）原始單炮頻率分析 反褶積前準備工作 盡可能壓制各種干擾 振幅補償 反Q濾波補償 做好野外靜校、折射靜校正反褶積的基本概念(子波整形 SHPFIL+SHDECN)脈沖反褶積 TVDCON預測反褶積 PRDCON兩步法統(tǒng)計子波反褶積 TSDCON地表一致性反褶積 SURSAN+SUEDEC疊后提高分辨率：SMDCON、ZPDCON GRISYS常用的反褶積模塊 預測反褶積 脈沖反褶積 前提均是子波最小相位，并假設反射系數(shù)為白噪，但是實際情況地震子波是混合相位，且頻帶有限,通過實驗得出脈沖反褶積壓縮子波功能最強，預測其次。 地下地震子

45、波是混合相位用脈沖反褶積； 目標為層間多次和海底全程多次，且信噪比低時，用預測反褶積。 GRISYS常用的反褶積模塊預測反褶積 預測反褶積是通過預測濾波來實現(xiàn)的，預測濾波的基本思想是，設計一個濾波器，使得該濾波器具有某種預測能力，通過它對信號的當前值和過去值的濾波，預測未來某個時刻將要出現(xiàn)的信號成分。預測時刻和當前時刻的距離稱為預測距離。預測濾波器的設計采用了最小平方準則，即要求預測結果和實際信號間誤差的最小平方和最小。 PRDCON是一個預測反褶積模塊，既可以用于消除地震記錄中的長、 短周期多次波，又可以用來提高地震記錄的分辨率。 GRISYS常用的反褶積模塊GRISYS常用的反褶積模塊原始

46、單炮脈沖反褶積預測反褶積 先在炮集記錄上提取統(tǒng)計子波，利用該子波做反褶積處理,消除炮點對子波的影響,第一步期望輸出波形為最小相位； 然后在共接收點道集上提取統(tǒng)計子波，利用該子波進行反褶積處理,消除共接收點對子波的影響,第二步期望輸出為零相位。 TSDCON不僅能校正、壓縮激發(fā)震源子波波形，還可以消除地表接收響應。從這個意義上講，它具有一定的地表一致性反褶積效果。兩步法子波反褶積 TSDCONTSDCONGRISYS常用的反褶積模塊：JOB P：MAQIN，L：301，T：OPTIONSG ：DBAMU VN01，HD2，(1-100)=M100X30， M500X570，M1000X1455，

47、: DISKI N DFN:GG05854，NTY2（NTY68），HD2=(1-350)， : TSDCON SG，NT96，XMIN25，XMAX1455 , MU01， (WN400-WN1100，WF600-WF1300， XN25-XF1455， T100-T1500，L120，WC1010，BR10-BR65， EP0)， (WN1200-WN2300，WF1400-WF2500， XN25-XF1455， T1700-T5000，L120，WC1015，BR10-BR55， EP0) : DSKOUT DF:GTSDCON1，F(xiàn)MT4: ENDTSDCON第一步GRISYS常用的

48、反褶積模塊：JOB P：MAQIN，L：301，T：OPTIONSG ：DBAMU VN01，HD2，(1-100)=M100X30， M500X570，M1000X1455，: DISKIN DFN:GTSDCON2，NTY18: TSDCON ZR，NT48，XMIN25，XMAX1455, MU01， (WN400-WN1100，WF600-WF1300， XN25-XF1455， T100-T1500，L120，WC1010，BR10-BR65， EP0)， (WN1200-WN2300，WF1400-WF2500， XN25-XF1455， T1700-T5000，L120，WC10

49、15，BR10-BR55， EP0) : DSKOUT DF:GTSDCON4: ENDTSDCON第二步GRISYS常用的反褶積模塊原始單炮脈沖反褶積預測反褶積GRISYS常用的反褶積模塊統(tǒng)計子波反褶積 反褶積參數(shù)選?。ㄒ詢刹椒槔?1）計算時窗,應用時窗 2）算子長度 1.52倍子波長度 3）白噪系數(shù) 10051020 4）期望輸出的子波相位 零相位、小相位 5）期望輸出的子波類型 帶通子波 脈沖 FC85Hz 雷克子波 FCpickup啟動拾取狀態(tài)并彈出層控制界面 GRISYS/WS-V8.0 改進內(nèi)容添加層拾取控制點填寫擬合參數(shù)保存參數(shù)交互擬合參數(shù)拾取 GRISYS/WS-V8.0 改進內(nèi)容交互擬合參數(shù)拾取參數(shù)拾取完成后，傳遞給參數(shù)試驗對話框 GRISYS/WS-V8.0 改進內(nèi)容交互擬合參數(shù)拾取擬合處理后對比RAPFIL 徑向預測濾波 是一個多道徑向預測濾波模塊，通過時移求出最大相關值所對