strata中文培訓(xùn)

1、1STRATA 培訓(xùn)2STRATA 培訓(xùn)STRATA 課程大綱課程大綱測井曲線時(shí)深對(duì)比-118練習(xí)3：Erskine 3D 測井曲線時(shí)深對(duì)比-131測井曲線層間插值建模-1563STRATA 培訓(xùn)地震反演基礎(chǔ)地震反演基礎(chǔ)反演反演：從地震數(shù)據(jù)中提取它所包含的潛在地質(zhì)信息的過程.傳統(tǒng)上，反演是在疊后地震數(shù)據(jù)上進(jìn)行,目的是提取聲波阻抗體聲波阻抗體.近年來，反演已經(jīng)被擴(kuò)展到疊前數(shù)據(jù)體,目的是既提取聲波阻抗又聲波阻抗又提取橫波阻抗體提取橫波阻抗體，從而可以推得孔隙流體。另外，最新的發(fā)展認(rèn)為可以利用反演結(jié)果直接預(yù)測巖性參數(shù)如孔隙度和含水飽和度。4STRATA 培訓(xùn)幾種不同反演方法幾種不同反演方法遞歸法遞歸

2、法: 傳統(tǒng)帶限反演傳統(tǒng)帶限反演顏色反演顏色反演: 遞歸反演的改進(jìn)遞歸反演的改進(jìn)稀疏脈沖法稀疏脈沖法(兩個(gè)不同算法兩個(gè)不同算法): 約束以獲得盡可能少的同相軸約束以獲得盡可能少的同相軸基于模型法基于模型法: 迭代更新層狀初始模型迭代更新層狀初始模型彈性波阻抗彈性波阻抗: 對(duì)對(duì)AVO數(shù)據(jù)的改進(jìn)數(shù)據(jù)的改進(jìn)LMR: 對(duì)對(duì)AVO數(shù)據(jù)的改進(jìn)數(shù)據(jù)的改進(jìn)聯(lián)合反演：聯(lián)合反演： 對(duì)對(duì)AVO數(shù)據(jù)的改進(jìn)數(shù)據(jù)的改進(jìn)在在STRATA軟件中提供以下這些反演方法軟件中提供以下這些反演方法(六種方法七個(gè)模塊六種方法七個(gè)模塊):5STRATA 培訓(xùn) 阻阻 抗抗縱波縱波橫波橫波彈性波彈性波反射系數(shù)反射系數(shù)子波子波地震道地震道所有反

3、演方法中共同的正演模型所有反演方法中共同的正演模型:反演中的一般正演模型反演中的一般正演模型6STRATA 培訓(xùn)波阻抗波阻抗反射系數(shù)反射系數(shù)反演反演子波子波地震道地震道反演試圖反推得到這個(gè)正演模型反演試圖反推得到這個(gè)正演模型:反演模型反演模型縱波、橫波、彈性阻抗縱波、橫波、彈性阻抗7STRATA 培訓(xùn)縱波阻抗縱波阻抗或或橫波阻抗橫波阻抗或或 彈性阻抗彈性阻抗11iiiiiIIRIII 阻抗阻抗 反射系數(shù)反射系數(shù)縱縱 波波 阻阻 抗抗 =橫橫 波波 阻阻 抗抗 =彈彈 性性 阻阻 抗抗 = PVSV復(fù)雜公式復(fù)雜公式 (后面介紹后面介紹)8STRATA 培訓(xùn)反射系數(shù)反射系數(shù) 地震道地震道*SWRN

4、oise地震道地震道 = 子波與反射系數(shù)的褶積子波與反射系數(shù)的褶積+噪音噪音.注釋注釋不存在模型化的多次波不存在模型化的多次波.沒考慮傳播損失和幾何擴(kuò)散沒考慮傳播損失和幾何擴(kuò)散.沒考慮頻率吸收沒考慮頻率吸收.(1) 子波可能是時(shí)變的子波可能是時(shí)變的.9STRATA 培訓(xùn)反射系數(shù)與子波褶積反射系數(shù)與子波褶積的結(jié)果是移除了大量的結(jié)果是移除了大量高頻細(xì)節(jié)高頻細(xì)節(jié):10STRATA 培訓(xùn)在時(shí)間域中的褶積就是頻率在時(shí)間域中的褶積就是頻率域中的乘積域中的乘積.從這些圖中可以看出從這些圖中可以看出,子波的子波的作用是將地震頻譜中高頻和作用是將地震頻譜中高頻和低頻都消除了低頻都消除了.理論上講理論上講,反演就

5、是試圖將這反演就是試圖將這些失去的頻率區(qū)域進(jìn)行恢復(fù)些失去的頻率區(qū)域進(jìn)行恢復(fù).11STRATA 培訓(xùn)反演中反演中“多解性多解性”問題問題所有反演算法都有多解性問題所有反演算法都有多解性問題.存在多于一種地質(zhì)模型可以與地震數(shù)據(jù)存在多于一種地質(zhì)模型可以與地震數(shù)據(jù)相一致，要在這些可能的模型中決定一相一致，要在這些可能的模型中決定一個(gè)模型，就需要去選擇地震數(shù)據(jù)以外的個(gè)模型，就需要去選擇地震數(shù)據(jù)以外的的一些信息。的一些信息。通常用以下兩種方法來使用地震以外的通常用以下兩種方法來使用地震以外的信息：信息： 初始猜測模型的建立初始猜測模型的建立 對(duì)最終結(jié)果與初始猜測模型背離對(duì)最終結(jié)果與初始猜測模型背離幅度進(jìn)行

6、約束幅度進(jìn)行約束最后反演結(jié)果既依賴于最后反演結(jié)果既依賴于 “其它信息其它信息”也取也取決于地震數(shù)據(jù)決于地震數(shù)據(jù).12STRATA 培訓(xùn)帶限帶限 (遞歸遞歸) 反演反演遞歸反演遞歸反演, 也稱為帶限反演是最簡單和最早一種反演方法也稱為帶限反演是最簡單和最早一種反演方法.11-iiiiiZZrZZ1*11-iiiirZZr1*11-inirZZr從第一層開始從第一層開始, 可以用這個(gè)公式遞推后續(xù)可以用這個(gè)公式遞推后續(xù)每層的阻抗每層的阻抗. :ith +1 層的阻抗可以從上一層層的阻抗可以從上一層 ith 確定確定:從定義反射系數(shù)開始從定義反射系數(shù)開始:13STRATA 培訓(xùn)在這個(gè)簡單的實(shí)例中可知在

7、這個(gè)簡單的實(shí)例中可知:(a) 告訴我們?nèi)绻o出一個(gè)脈沖，告訴我們?nèi)绻o出一個(gè)脈沖，那么我們就可以恢復(fù)阻抗的真實(shí)值那么我們就可以恢復(fù)阻抗的真實(shí)值, 但是但是(b) 中看到如果我們用一個(gè)子波與中看到如果我們用一個(gè)子波與脈沖進(jìn)行褶積，那就無法恢復(fù)阻抗脈沖進(jìn)行褶積，那就無法恢復(fù)阻抗的真實(shí)值。的真實(shí)值。Z2 = 818Z2 = 1500Z1 = 1000Z3 = 1227Z4 = 1004Z1 = 1000帶限反演帶限反演14STRATA 培訓(xùn)輸入地震數(shù)據(jù)輸入地震數(shù)據(jù)遞歸反演數(shù)據(jù)遞歸反演數(shù)據(jù)遞歸反演遞歸反演獲得獲得的結(jié)果與輸入的結(jié)果與輸入的地震數(shù)據(jù)頻的地震數(shù)據(jù)頻帶一樣帶一樣.從圖中可以看從圖中可以看出出

8、,與測井比較與測井比較反演結(jié)果損失反演結(jié)果損失了高頻細(xì)節(jié)了高頻細(xì)節(jié).15STRATA 培訓(xùn)步驟步驟1:遞歸反演遞歸反演 的初始背景模型是通過對(duì)井阻抗濾的初始背景模型是通過對(duì)井阻抗濾波獲得波獲得:10-Hz 高截濾波高截濾波16STRATA 培訓(xùn)步驟步驟 2:對(duì)地震道運(yùn)用遞歸運(yùn)算對(duì)地震道運(yùn)用遞歸運(yùn)算. (注釋注釋: 這幾乎等同于在相位域進(jìn)這幾乎等同于在相位域進(jìn)行行-90 度轉(zhuǎn)換度轉(zhuǎn)換):1*11-iiiirZZr17STRATA 培訓(xùn)步驟步驟 3:濾波后的模型加上濾波后的模型加上比例化比例化反演道集就得到了最終結(jié)果反演道集就得到了最終結(jié)果:+=18STRATA 培訓(xùn)遞歸反演中幾個(gè)要點(diǎn)遞歸反演中幾

9、個(gè)要點(diǎn):不考慮子波不考慮子波. 這意味輸入的地震數(shù)據(jù)必須是零相位這意味輸入的地震數(shù)據(jù)必須是零相位. 如果有如果有提取子波可行提取子波可行Strata 軟件可以自動(dòng)地對(duì)地震數(shù)據(jù)軟件可以自動(dòng)地對(duì)地震數(shù)據(jù) “去相位去相位”.即使地震數(shù)據(jù)是零相位即使地震數(shù)據(jù)是零相位, 實(shí)際子波的旁瓣實(shí)際子波的旁瓣在遞歸算法中也會(huì)在遞歸算法中也會(huì)被解釋為巖性的變化被解釋為巖性的變化.反演的結(jié)果其頻帶與地震數(shù)據(jù)反演的結(jié)果其頻帶與地震數(shù)據(jù)頻帶一樣頻帶一樣.(1) 對(duì)對(duì)地震數(shù)據(jù)的加權(quán)值地震數(shù)據(jù)的加權(quán)值(與反射系數(shù)匹配與反射系數(shù)匹配)決定了遞歸反演的結(jié)決定了遞歸反演的結(jié)果阻抗值范圍是否合適果阻抗值范圍是否合適.19STRATA

10、 培訓(xùn)基于模型的反演基于模型的反演*SWRNoise基于模型反演基于模型反演 是從褶積模型方程式開始的是從褶積模型方程式開始的:假定地震道假定地震道S,和子波和子波 W, 是已知的是已知的.假定噪音是隨機(jī)的并與地震信號(hào)不相關(guān)假定噪音是隨機(jī)的并與地震信號(hào)不相關(guān).求解反射系數(shù)求解反射系數(shù), R, 以滿足這個(gè)等式以滿足這個(gè)等式. 這是一個(gè)非線性問題這是一個(gè)非線性問題,所所以求解過程是迭代進(jìn)行的以求解過程是迭代進(jìn)行的.20STRATA 培訓(xùn)步驟步驟 1:基于模型反演基于模型反演 的初始背景模型是通過對(duì)井阻抗分段形的初始背景模型是通過對(duì)井阻抗分段形成的成的:用戶用毫秒用戶用毫秒(MS)來定義層分塊大小來

11、定義層分塊大小.所有的層分段開始時(shí)都被設(shè)成一樣所有的層分段開始時(shí)都被設(shè)成一樣大小大小(用毫秒用毫秒).21STRATA 培訓(xùn)步驟步驟 2:利用層狀模型和利用層狀模型和已知已知的子波計(jì)算人工合成地震道的子波計(jì)算人工合成地震道.這是與實(shí)際地震道的比較這是與實(shí)際地震道的比較.通過分析人工合成地震記錄通過分析人工合成地震記錄與實(shí)際地震道的誤差，每層與實(shí)際地震道的誤差，每層(分段分段)的厚度和震幅值都進(jìn)行的厚度和震幅值都進(jìn)行修改以減少誤差修改以減少誤差.這個(gè)步驟通過一系列遞歸不這個(gè)步驟通過一系列遞歸不斷重復(fù)斷重復(fù).SyntheticSeismic22STRATA 培訓(xùn)輸入地震數(shù)據(jù)輸入地震數(shù)據(jù)基于模型反演

12、基于模型反演基于模型反演基于模型反演 獲得了寬頻的結(jié)獲得了寬頻的結(jié)果果.問題是高頻成分問題是高頻成分可能來自于初始可能來自于初始猜測模型而不是猜測模型而不是來自地震數(shù)據(jù)來自地震數(shù)據(jù).23STRATA 培訓(xùn)遞歸反演遞歸反演基于模型反演基于模型反演這是遞歸反演和基這是遞歸反演和基于模型反演的比較于模型反演的比較.通常基于模型反演通?；谀Ｐ头囱菘梢垣@得更細(xì)致信可以獲得更細(xì)致信息息,但這個(gè)結(jié)果實(shí)際但這個(gè)結(jié)果實(shí)際上相當(dāng)類似上相當(dāng)類似.24STRATA 培訓(xùn)基于模型反演的幾個(gè)要點(diǎn)基于模型反演的幾個(gè)要點(diǎn):由于已經(jīng)知道子波由于已經(jīng)知道子波,在計(jì)算過程中它的影響從地震數(shù)據(jù)中排除了在計(jì)算過程中它的影響從地震數(shù)據(jù)

13、中排除了.也就是也就是說地震數(shù)據(jù)沒有必要是零相位說地震數(shù)據(jù)沒有必要是零相位,而只要而只要子波與地震數(shù)據(jù)等相位子波與地震數(shù)據(jù)等相位就可以就可以.估算的估算的子波若有誤差子波若有誤差,將導(dǎo)致反演結(jié)果出差錯(cuò)將導(dǎo)致反演結(jié)果出差錯(cuò).地震有效地震有效分辨率分辨率得到提高得到提高.反演結(jié)果可能在很大程度上依賴原始猜測初始模型反演結(jié)果可能在很大程度上依賴原始猜測初始模型.解決的方法是對(duì)初始解決的方法是對(duì)初始模型模型進(jìn)行濾波進(jìn)行濾波.(1) 與其它反演方法一樣與其它反演方法一樣,存在存在多解性多解性問題問題.25STRATA 培訓(xùn)稀疏脈沖反演稀疏脈沖反演稀疏脈沖反演稀疏脈沖反演 假定實(shí)際假定實(shí)際反射可以認(rèn)為是由

14、一系反射可以認(rèn)為是由一系列大脈沖里夾雜有小脈列大脈沖里夾雜有小脈沖背景沖背景.稀疏脈沖反演稀疏脈沖反演假定假定 只有只有大脈沖有意義大脈沖有意義.該方法通該方法通過檢查地震道來尋找大脈過檢查地震道來尋找大脈沖的位置沖的位置.26STRATA 培訓(xùn)稀疏脈沖反演稀疏脈沖反演 每次建每次建立一個(gè)脈沖反射序列，立一個(gè)脈沖反射序列， 通過增加脈沖直到地震通過增加脈沖直到地震道被足夠準(zhǔn)確地反演完道被足夠準(zhǔn)確地反演完成。成。阻抗塊的振幅值是由阻抗塊的振幅值是由基基于模型反演于模型反演算法來確定算法來確定27STRATA 培訓(xùn)輸入地震輸入地震稀疏脈沖反演稀疏脈沖反演稀疏脈沖反演稀疏脈沖反演生成一個(gè)寬帶生成一個(gè)

15、寬帶高頻反演結(jié)果高頻反演結(jié)果.28STRATA 培訓(xùn)基于模型反演基于模型反演稀疏脈沖反演稀疏脈沖反演稀疏脈沖反演稀疏脈沖反演結(jié)果與結(jié)果與基于模基于模型反演型反演結(jié)果類結(jié)果類似似.主要區(qū)別是缺主要區(qū)別是缺少非常薄層的少非常薄層的細(xì)節(jié)細(xì)節(jié).29STRATA 培訓(xùn)稀疏脈沖反演幾個(gè)要點(diǎn)稀疏脈沖反演幾個(gè)要點(diǎn):只有當(dāng)?shù)卣饠?shù)據(jù)只有當(dāng)?shù)卣饠?shù)據(jù)有脈沖有脈沖，稀疏脈沖反演才得到同相軸；，稀疏脈沖反演才得到同相軸；它試圖利用地震數(shù)據(jù)獲得它試圖利用地震數(shù)據(jù)獲得最簡單的可能模型最簡單的可能模型；通常得到的結(jié)果比地質(zhì)本身實(shí)際的通常得到的結(jié)果比地質(zhì)本身實(shí)際的同相軸少同相軸少；它比它比基于模型反演基于模型反演更更少的依賴于初

16、始猜測模型少的依賴于初始猜測模型。30STRATA 培訓(xùn)顏色反演顏色反演顏色反演顏色反演 是對(duì)是對(duì)遞歸反演遞歸反演的修改的修改, 最早是由最早是由BP公司的公司的 Lancaster and Whitcombe 在在2000 年年SEG年會(huì)上提出來的年會(huì)上提出來的.這種反演過程這種反演過程, 尋找一個(gè)簡單操作因子尋找一個(gè)簡單操作因子 O, 對(duì)地震道進(jìn)行直接轉(zhuǎn)換求得反演對(duì)地震道進(jìn)行直接轉(zhuǎn)換求得反演結(jié)果結(jié)果:*IO S作者在頻率域確定簡單操作因子作者在頻率域確定簡單操作因子, O . 通過比較地震數(shù)據(jù)和實(shí)際反演結(jié)果，認(rèn)為簡單操作因子的相位是通過比較地震數(shù)據(jù)和實(shí)際反演結(jié)果，認(rèn)為簡單操作因子的相位是-

17、90度。度。31STRATA 培訓(xùn)簡單操作因子震幅譜用這種方式來求取簡單操作因子震幅譜用這種方式來求取:利用工區(qū)的一組井利用工區(qū)的一組井, 所所有井的聲阻抗的振幅譜有井的聲阻抗的振幅譜做交匯圖做交匯圖.通過理論預(yù)測通過理論預(yù)測,我們可以我們可以擬合一條直線來擬合一條直線來 代表代表“理想理想” 輸出阻抗譜輸出阻抗譜.Log(Frequency)Log(Impedance)聲阻抗震幅譜聲阻抗震幅譜32STRATA 培訓(xùn)然后然后, 利用井旁一組地震道利用井旁一組地震道, 求求取取平均地震譜平均地震譜.地震譜地震譜Frequency (Hz)操作因子譜操作因子譜從以上兩個(gè)頻譜從以上兩個(gè)頻譜,求取操作

18、因子譜求取操作因子譜. 使得地震譜形狀在地震頻帶內(nèi)轉(zhuǎn)換使得地震譜形狀在地震頻帶內(nèi)轉(zhuǎn)換為為阻抗譜阻抗譜.33STRATA 培訓(xùn)顏色反演操作因子顏色反演操作因子Time (ms)將獲得的震幅譜做將獲得的震幅譜做-90 度相移就生成了度相移就生成了顏色反顏色反演操作因子演操作因子.將它與所有地震道褶積將它與所有地震道褶積.34STRATA 培訓(xùn)輸入地震輸入地震顏色反演顏色反演顏色反演顏色反演 結(jié)果非常類似結(jié)果非常類似于遞歸反演于遞歸反演.其中一個(gè)區(qū)別是顏色反其中一個(gè)區(qū)別是顏色反演的結(jié)果是相對(duì)聲阻抗演的結(jié)果是相對(duì)聲阻抗,既有正值也有負(fù)值既有正值也有負(fù)值.-30000+300035STRATA 培訓(xùn)遞歸

19、反演遞歸反演顏色反演顏色反演-30000+30004600830012000相對(duì)相對(duì) AI絕對(duì)絕對(duì) AI36STRATA 培訓(xùn)顏色反演幾點(diǎn)小結(jié)顏色反演幾點(diǎn)小結(jié)除了定義整體阻抗趨勢外除了定義整體阻抗趨勢外,與初始模型基本無關(guān)與初始模型基本無關(guān).運(yùn)算非常快運(yùn)算非?？?用戶所需要定義的參數(shù)非常少用戶所需要定義的參數(shù)非常少.假定地震數(shù)據(jù)是零相位的假定地震數(shù)據(jù)是零相位的.在最初實(shí)施中，該方法可以獲得在最初實(shí)施中，該方法可以獲得相對(duì)相對(duì)阻抗體，雖然我們可阻抗體，雖然我們可以選擇添加低頻阻抗背景趨勢。以選擇添加低頻阻抗背景趨勢。37STRATA 培訓(xùn)AVO 反演反演基本褶積模型反演假定的是零偏移距地震數(shù)據(jù)基

20、本褶積模型反演假定的是零偏移距地震數(shù)據(jù).褶積反演法不能運(yùn)用到帶有褶積反演法不能運(yùn)用到帶有AVO響應(yīng)的地震數(shù)據(jù)響應(yīng)的地震數(shù)據(jù), 因?yàn)轳薹e反演法沒有明因?yàn)轳薹e反演法沒有明確地考慮確地考慮 VP/VS 的變化的變化.將反演方法延伸到疊前將反演方法延伸到疊前AVO 地震數(shù)據(jù)的處理地震數(shù)據(jù)的處理, 目前這些算法主要運(yùn)用于目前這些算法主要運(yùn)用于:彈性阻抗反演彈性阻抗反演Lambda-Mu-Rho (LMR)反演反演聯(lián)合反演聯(lián)合反演38STRATA 培訓(xùn)彈性阻抗彈性阻抗 222tansinsinCBAR彈性阻抗概念最初由彈性阻抗概念最初由Connolly 先生提出先生提出(The Leading Edge,

21、 18, no. 4, 438-452 (1999).他利用他利用Aki-Richards 方程方程,將反射震幅與入射角聯(lián)系起來將反射震幅與入射角聯(lián)系起來:,21:PPVVAwhere222421PSSSPSPPVVVVVVVVB.21:PPVVCand注釋注釋:在常規(guī)反演理論中假定在常規(guī)反演理論中假定 B=C=0, 而沒而沒有考慮有考慮VP/VS 的變化的變化.39STRATA 培訓(xùn) 11022PPVAIRAVAI注意注意, 對(duì)于零偏移距對(duì)于零偏移距:通過類推通過類推, Connolly 定義了一種新阻抗類型定義了一種新阻抗類型: EIEIEIRln2121222(1 tan)( 8sin)

22、(1 4sin)( )KKPSEIVV22PSVVK 其中通過數(shù)學(xué)處理通過數(shù)學(xué)處理,最終最終:40STRATA 培訓(xùn)這張圖重疊顯示了同一口井位置彈性阻抗與聲波阻抗這張圖重疊顯示了同一口井位置彈性阻抗與聲波阻抗.在碳在碳?xì)渚奂瘏^(qū)域彈性阻抗值顯示異常低值氫聚集區(qū)域彈性阻抗值顯示異常低值.41STRATA 培訓(xùn)道集AVO 分析近角度疊加 stack at 1遠(yuǎn)角度疊加 stack at 2反演成彈性阻抗 EI(1)反演成彈性阻抗 EI(2)這種反演的工作流程是從迭前道集開始，首先這種反演的工作流程是從迭前道集開始，首先產(chǎn)生兩組角度疊加剖面產(chǎn)生兩組角度疊加剖面, 分別對(duì)其進(jìn)行反演分別對(duì)其進(jìn)行反演.42

23、STRATA 培訓(xùn)這就獲得兩種反演結(jié)果這就獲得兩種反演結(jié)果:遠(yuǎn)角度反演遠(yuǎn)角度反演近角度反演近角度反演43STRATA 培訓(xùn)近角度反演與遠(yuǎn)角度反近角度反演與遠(yuǎn)角度反演進(jìn)行交繪顯示演進(jìn)行交繪顯示.遠(yuǎn)角度反演異常低值區(qū)遠(yuǎn)角度反演異常低值區(qū)可被解釋為潛在含有碳可被解釋為潛在含有碳?xì)渚酆衔餁渚酆衔?近角度反演近角度反演遠(yuǎn)角度反演遠(yuǎn)角度反演44STRATA 培訓(xùn)將交繪圖中的異常低將交繪圖中的異常低值點(diǎn)從交繪圖中成像值點(diǎn)從交繪圖中成像到原始地震數(shù)據(jù)剖面到原始地震數(shù)據(jù)剖面上以進(jìn)行解釋上以進(jìn)行解釋.45STRATA 培訓(xùn)Lambda-Mu-Rho (LMR)2S2P2P2P2S2SSPZ2Z:so)2()V(Z

24、:and)V(Z:thereforeVand2VLMR方法運(yùn)用了方法運(yùn)用了VP, VS, 之間的關(guān)系式及拉梅常數(shù)之間的關(guān)系式及拉梅常數(shù) 和和 :注意：最終結(jié)果表達(dá)式是注意：最終結(jié)果表達(dá)式是 和和 ，可以通過縱波，可以通過縱波、橫波阻抗、橫波阻抗 ZP和和ZS的關(guān)系的關(guān)系式來表示。式來表示。LMR 方法最初是由方法最初是由 Goodway et al (SEG Expanded Abstracts, 1997)提出來的提出來的。類似于彈性阻抗法，這種方法對(duì)處理含。類似于彈性阻抗法，這種方法對(duì)處理含AVO效應(yīng)的數(shù)據(jù)比常規(guī)反演更有利。效應(yīng)的數(shù)據(jù)比常規(guī)反演更有利。46STRATA 培訓(xùn)GathersA

25、VO AnalysisRP Estimate RS EstimateCross-plotInvert to ZPInvert to ZSTransform to and LMR反演流程首先從疊前數(shù)據(jù)分反演流程首先從疊前數(shù)據(jù)分析得到縱橫波反射系數(shù)體析得到縱橫波反射系數(shù)體RP 和和 RS對(duì)對(duì)RP 和和 RS分別反演會(huì)生成縱橫波分別反演會(huì)生成縱橫波阻抗體阻抗體ZP 和和 ZS .這些數(shù)據(jù)體會(huì)通過類似這些數(shù)據(jù)體會(huì)通過類似Goodway等關(guān)系式或交匯圖進(jìn)行轉(zhuǎn)換。等關(guān)系式或交匯圖進(jìn)行轉(zhuǎn)換。LMR 反演流程反演流程47STRATA 培訓(xùn)如圖展示了對(duì)如圖展示了對(duì)Albert地區(qū)地區(qū)典型含氣砂巖運(yùn)用典型含氣砂巖

26、運(yùn)用LMR法法得到的結(jié)果，圖中橢圓內(nèi)得到的結(jié)果，圖中橢圓內(nèi)指示為含氣砂巖。指示為含氣砂巖。 上剖面為上剖面為體體, 下剖面為下剖面為體體 LMR 反演例子反演例子48STRATA 培訓(xùn)左圖為左圖為與與交匯圖，紅色區(qū)域指示氣層交匯圖，紅色區(qū)域指示氣層 (低低) 藍(lán)色區(qū)域指示泥巖和水砂巖。藍(lán)色區(qū)域指示泥巖和水砂巖。這些區(qū)域投到剖面上如下圖所示，可以指示含這些區(qū)域投到剖面上如下圖所示，可以指示含氣砂巖區(qū)域（紅色部分）氣砂巖區(qū)域（紅色部分）lambda-rhomu-rhoLMR 反演實(shí)例反演實(shí)例49STRATA 培訓(xùn)聯(lián)合反演聯(lián)合反演聯(lián)合反演聯(lián)合反演可以對(duì)輸入的疊前角道集數(shù)據(jù)反演得到可以對(duì)輸入的疊前角道

27、集數(shù)據(jù)反演得到ZP, ZS, 或者可能密度體。這種方或者可能密度體。這種方法的優(yōu)點(diǎn)在于它可以將這些參數(shù)進(jìn)行約束，這樣就可以使結(jié)果更穩(wěn)定并且盡量避免法的優(yōu)點(diǎn)在于它可以將這些參數(shù)進(jìn)行約束，這樣就可以使結(jié)果更穩(wěn)定并且盡量避免多解性問題。多解性問題。123( )PPPSDRc Rc Rc R2122222231tan8sin1tan2sin2SPcccVV 1212.PPPSSSDVRVVRVRwhere:我們從 Fattis version of the Aki-Richards 方程開始。這個(gè)公式用角度作為變量來表示反射系數(shù)：50STRATA 培訓(xùn)聯(lián)合反演聯(lián)合反演Ln(Zs)Ln(Zp)Ln()L

28、n(Zp)SLDLln()ln()ln( )ln()SPcSPcDZkZkLmZmL 假定水砂巖可以用一個(gè)線性關(guān)系式來表示。假定水砂巖可以用一個(gè)線性關(guān)系式來表示。聯(lián)合反演聯(lián)合反演尋找背景趨勢的偏差尋找背景趨勢的偏差 :51STRATA 培訓(xùn)聯(lián)合反演聯(lián)合反演聯(lián)合反演聯(lián)合反演 生成生成Zp、Zs、Density體體,以及以及推導(dǎo)得到的各種推導(dǎo)得到的各種組合。組合。52STRATA 培訓(xùn)聯(lián)合反演聯(lián)合反演聯(lián)合反演的解釋與其它聯(lián)合反演的解釋與其它AVO反演結(jié)果類似：反演結(jié)果類似：53STRATA 培訓(xùn)它是地質(zhì)統(tǒng)計(jì)反演的一種形式，可以通過對(duì)輸入地震數(shù)據(jù)體運(yùn)算產(chǎn)生大量的它是地質(zhì)統(tǒng)計(jì)反演的一種形式，可以通過對(duì)

29、輸入地震數(shù)據(jù)體運(yùn)算產(chǎn)生大量的反演結(jié)果以此來解決多解性的問題。每一個(gè)運(yùn)算結(jié)果要與地震數(shù)據(jù)一致，并反演結(jié)果以此來解決多解性的問題。每一個(gè)運(yùn)算結(jié)果要與地震數(shù)據(jù)一致，并且保持如方差圖中所包含的所期望的連續(xù)條件。通過分析這些結(jié)果以估算得且保持如方差圖中所包含的所期望的連續(xù)條件。通過分析這些結(jié)果以估算得到一個(gè)非確定性結(jié)果，以及最大可能結(jié)果。到一個(gè)非確定性結(jié)果，以及最大可能結(jié)果。隨機(jī)反演隨機(jī)反演Stochastic InversionThis is a form of geostatistical inversion which explicitly addresses the non-uniqueness

30、 problem by producing a large range of inversion results for a given input seismic volume. Each of the results is consistent with the seismic data, and honors the expected continuity conditions, as contained in the variograms. These results are analyzed to give an estimate of the uncertainty in the

31、result, along with the most probable result.54STRATA 培訓(xùn)一般地震反演流程一般地震反演流程(1) 建模型建模型選擇測井選擇測井每口井進(jìn)行相關(guān)每口井進(jìn)行相關(guān)提取子波提取子波讀取讀取 / 拾取地震反演拾取地震反演(2) 完成反演完成反演選擇反演方法和參數(shù)選擇反演方法和參數(shù)對(duì)反演結(jié)果進(jìn)行對(duì)反演結(jié)果進(jìn)行QC(3) 解釋反演結(jié)果解釋反演結(jié)果做數(shù)據(jù)切片做數(shù)據(jù)切片做交繪圖做交繪圖輸出到輸出到EMERGE 軟件軟件55STRATA 培訓(xùn)反演結(jié)果反演結(jié)果QC兩種方法兩種方法:誤差圖誤差圖(1) 交互驗(yàn)證交互驗(yàn)證我們?nèi)绾沃婪囱萁Y(jié)果是我們?nèi)绾沃婪囱萁Y(jié)果是否有效

32、否有效?輸入地震輸入地震反演結(jié)果反演結(jié)果56STRATA 培訓(xùn)從阻抗道利用已知的子波從阻抗道利用已知的子波可以計(jì)算人工合成道可以計(jì)算人工合成道.理想情況理想情況, 這個(gè)人工合成這個(gè)人工合成道應(yīng)該與輸入地震道非常道應(yīng)該與輸入地震道非常相似相似.輸入地震輸入地震人工合成道人工合成道57STRATA 培訓(xùn)從輸入地震道中減去人工合從輸入地震道中減去人工合成道成道,就得到就得到反演誤差圖反演誤差圖.如果反演效果不錯(cuò)如果反演效果不錯(cuò),誤差圖上誤差圖上應(yīng)該只有非常小的震幅應(yīng)該只有非常小的震幅,并沒并沒有聚集在某處有聚集在某處.由于存在多解性，即使微小由于存在多解性，即使微小的誤差圖也不能保證反演結(jié)的誤差圖也

33、不能保證反演結(jié)果就是正確的果就是正確的.輸入地震輸入地震反演誤差反演誤差58STRATA 培訓(xùn)第二種質(zhì)量控制方法是第二種質(zhì)量控制方法是交互驗(yàn)證交互驗(yàn)證. 在這個(gè)驗(yàn)證中在這個(gè)驗(yàn)證中,我們?cè)诔跷覀冊(cè)诔跏寄Ｐ椭袆h除一口井始模型中刪除一口井,在其位置做反演在其位置做反演,然后將反演結(jié)果與隱然后將反演結(jié)果與隱藏井進(jìn)行比較藏井進(jìn)行比較.隱藏井隱藏井反演結(jié)果反演結(jié)果反演誤差反演誤差59STRATA 培訓(xùn)每口井分析誤差每口井分析誤差,以便識(shí)別有問題的井以便識(shí)別有問題的井.60STRATA 培訓(xùn)反演作為反演作為 EMERGE一種屬性一種屬性最新使用反演結(jié)果是將其做為最新使用反演結(jié)果是將其做為EMERGE軟件的軟

34、件的輸入輸入, EMERGE可以直接預(yù)測孔隙度和其它巖性可以直接預(yù)測孔隙度和其它巖性體體.反演反演EMERGE孔隙度體孔隙度體61STRATA 培訓(xùn)練習(xí)練習(xí) 1: 楔形模型練習(xí)楔形模型練習(xí)第一個(gè)練習(xí)對(duì)楔型模型做反演第一個(gè)練習(xí)對(duì)楔型模型做反演. 目的是在簡單例子上學(xué)習(xí)反演的基本步驟目的是在簡單例子上學(xué)習(xí)反演的基本步驟.在在WINDOWS(UNIX)上鍵入上鍵入geoview或或從從Windows啟動(dòng)啟動(dòng)GEOVIEW 軟件軟件.我們將創(chuàng)建一個(gè)新的我們將創(chuàng)建一個(gè)新的GEOVIEW 數(shù)據(jù)庫數(shù)據(jù)庫. 選擇選擇New, 如顯示如顯示, 然后在顯示的第一個(gè)然后在顯示的第一個(gè)菜單上點(diǎn)擊菜單上點(diǎn)擊 OK（如上

35、圖）（如上圖）點(diǎn)擊點(diǎn)擊Well Explorer出現(xiàn)如下圖所示界面出現(xiàn)如下圖所示界面。62STRATA 培訓(xùn)現(xiàn)在我們?yōu)檫@個(gè)工區(qū)創(chuàng)建一個(gè)新的數(shù)據(jù)庫現(xiàn)在我們?yōu)檫@個(gè)工區(qū)創(chuàng)建一個(gè)新的數(shù)據(jù)庫如果如果Geoview工區(qū)已打開可工區(qū)已打開可以點(diǎn)以點(diǎn) Database / New創(chuàng)建創(chuàng)建如果這是如果這是GEOVIEW第一次在第一次在這臺(tái)機(jī)器上運(yùn)行這臺(tái)機(jī)器上運(yùn)行 ，你可以看，你可以看到如右圖窗口，選擇到如右圖窗口，選擇New并并點(diǎn)擊點(diǎn)擊OK.63STRATA 培訓(xùn)將數(shù)據(jù)庫命名將數(shù)據(jù)庫命名 wedge_database Ok:現(xiàn)在，可以看到現(xiàn)在，可以看到 Well Explorer 中沒有井顯示：中沒有井顯示：6

36、4STRATA 培訓(xùn)現(xiàn)在，為這個(gè)數(shù)據(jù)庫添加井曲線，如下步驟：現(xiàn)在，為這個(gè)數(shù)據(jù)庫添加井曲線，如下步驟：點(diǎn)點(diǎn) Import Data:選擇選擇 Wedge_well_logs.las 文文件并點(diǎn)件并點(diǎn)Next65STRATA 培訓(xùn)在接下來的菜單中，選擇默認(rèn)的即可：在接下來的菜單中，選擇默認(rèn)的即可：Note that we are creating a new well called Wedge_well:最后，我們可以讀取聲波曲線等如下：最后，我們可以讀取聲波曲線等如下：在這個(gè)菜單在這個(gè)菜單, 選擇默認(rèn)的參數(shù)即可，如選擇默認(rèn)的參數(shù)即可，如XY坐標(biāo)：坐標(biāo)：66STRATA 培訓(xùn)接下來點(diǎn)接下來點(diǎn)OK，

37、出現(xiàn)如下提示窗口，問你選擇哪種單位，出現(xiàn)如下提示窗口，問你選擇哪種單位， 點(diǎn)擊點(diǎn)擊 Ok 即可：即可：67STRATA 培訓(xùn)當(dāng)曲線加載完成后，出現(xiàn)如下窗口：當(dāng)曲線加載完成后，出現(xiàn)如下窗口：點(diǎn)點(diǎn) “1” 選擇這口井選擇這口井, 然后點(diǎn)然后點(diǎn)Display well. 曲線就會(huì)曲線就會(huì)顯示出來。顯示出來。68STRATA 培訓(xùn)接下來啟動(dòng)接下來啟動(dòng)Strata，選擇，選擇 GEOVIEW 主窗口上的主窗口上的 STRATA 按鈕按鈕.選擇選擇 Start New Project:給新工區(qū)命名：給新工區(qū)命名： “wedge_project”: 69STRATA 培訓(xùn)在在STRATA 窗口上窗口上,通過

38、通過 Data Manager / Import Data / Open Seismic / From SEG-Y File 加載加載地震數(shù)據(jù)地震數(shù)據(jù):在在File Selection 菜單菜單, 選擇文件選擇文件 “wedge.sgy”. 點(diǎn)擊點(diǎn)擊 Next :70STRATA 培訓(xùn)這菜單有幾頁來定義文件格式和數(shù)據(jù)的幾何形狀這菜單有幾頁來定義文件格式和數(shù)據(jù)的幾何形狀(觀測系統(tǒng)觀測系統(tǒng)). 下一頁下一頁,盡管數(shù)據(jù)是一條盡管數(shù)據(jù)是一條2D測線測線,我們可以建立我們可以建立3D觀測系統(tǒng)觀測系統(tǒng).接下來接下來, 需要告訴程序在道頭中沒有需要告訴程序在道頭中沒有Inline & Xline 數(shù)

39、數(shù) 也沒有也沒有 X & Y坐標(biāo)坐標(biāo).71STRATA 培訓(xùn)點(diǎn)擊點(diǎn)擊Next 兩次兩次. 地地震數(shù)據(jù)掃描后震數(shù)據(jù)掃描后,出現(xiàn)這出現(xiàn)這樣界面樣界面:電擊電擊Ok 利用這些缺省利用這些缺省值值. 后面的練習(xí)將更后面的練習(xí)將更仔細(xì)地觀測分析地震仔細(xì)地觀測分析地震數(shù)據(jù)的加載數(shù)據(jù)的加載.72STRATA 培訓(xùn)現(xiàn)在可以看見這個(gè)楔型模型。現(xiàn)在可以看見這個(gè)楔型模型。另外，還會(huì)有如下圖所示窗口另外，還會(huì)有如下圖所示窗口顯示出來，定義井在地震上的顯示出來，定義井在地震上的位置，可以在位置，可以在Xline處填處填45，然后點(diǎn)擊然后點(diǎn)擊OK即可。即可。73STRATA 培訓(xùn)這口這口井就會(huì)顯示在這個(gè)模型上井就

40、會(huì)顯示在這個(gè)模型上 然而，然而，井曲線并不一定對(duì)應(yīng)準(zhǔn)確，所以需要井曲線并不一定對(duì)應(yīng)準(zhǔn)確，所以需要做校正，點(diǎn)做校正，點(diǎn)Well / Edit/Correlate Well:讓你選擇井來進(jìn)行校正，選擇這口讓你選擇井來進(jìn)行校正，選擇這口井，然后點(diǎn)井，然后點(diǎn)Edit.74STRATA 培訓(xùn)eLog 窗口就會(huì)出現(xiàn)窗口就會(huì)出現(xiàn)，右側(cè)顯示的是模，右側(cè)顯示的是模型數(shù)據(jù)型數(shù)據(jù)75STRATA 培訓(xùn)eLog 軟件窗口彈出軟件窗口彈出(顯示井曲線顯示井曲線 wedge log).完成井曲線對(duì)比完成井曲線對(duì)比,在井口位置求一個(gè)平均在井口位置求一個(gè)平均地震復(fù)合道地震復(fù)合道. 缺省值是在井口缺省值是在井口+/-1周圍道周

41、圍道求平均地震道求平均地震道.點(diǎn)擊點(diǎn)擊Ok 使用缺省值使用缺省值.進(jìn)行井曲線對(duì)比，點(diǎn)擊進(jìn)行井曲線對(duì)比，點(diǎn)擊Correlate 按鈕按鈕:76STRATA 培訓(xùn)eLog 軟件窗口顯示如下軟件窗口顯示如下: 藍(lán)色顯示人工合成記錄藍(lán)色顯示人工合成記錄道道. 紅色顯示復(fù)合地震道紅色顯示復(fù)合地震道. 原始地震道用黑色顯示原始地震道用黑色顯示.在對(duì)井曲線進(jìn)行對(duì)比前需在對(duì)井曲線進(jìn)行對(duì)比前需要一個(gè)更準(zhǔn)確的子波要一個(gè)更準(zhǔn)確的子波.77STRATA 培訓(xùn)用于計(jì)算合成記錄的子波是缺省的稱為用于計(jì)算合成記錄的子波是缺省的稱為“current wavelet”. 查看缺省子波查看缺省子波, 點(diǎn)擊點(diǎn)擊 Wavelet /

42、 Display Current Wavelet:缺省子波是帶通濾波缺省子波是帶通濾波(頻率頻率: 5/10 - 50/60. ) 點(diǎn)擊點(diǎn)擊Frequency 標(biāo)簽可在頻標(biāo)簽可在頻率域顯示子波率域顯示子波.78STRATA 培訓(xùn)求取求取 “統(tǒng)計(jì)統(tǒng)計(jì)” 子波（不使用井信息），生成零相位子波子波（不使用井信息），生成零相位子波,其其頻率域與地震數(shù)據(jù)的響應(yīng)一致，點(diǎn)擊頻率域與地震數(shù)據(jù)的響應(yīng)一致，點(diǎn)擊 Wavelet / Extract Wavelet / Statistical:在子波提取的界面在子波提取的界面,第一頁是設(shè)置分第一頁是設(shè)置分析窗口析窗口. 使用缺省值使用缺省值. 點(diǎn)擊點(diǎn)擊Next 兩次

43、兩次.79STRATA 培訓(xùn)使用缺省值使用缺省值. 這些參數(shù)后面再討論這些參數(shù)后面再討論.點(diǎn)擊點(diǎn)擊OK 提取統(tǒng)計(jì)子波提取統(tǒng)計(jì)子波.新子波如左圖所示新子波如左圖所示. 檢查它的頻率響應(yīng)檢查它的頻率響應(yīng),點(diǎn)點(diǎn)擊擊Frequency 標(biāo)簽標(biāo)簽.80STRATA 培訓(xùn)在窗口的底部子波名字也變了在窗口的底部子波名字也變了.井對(duì)比的過程是將藍(lán)色的合成記錄道上的井對(duì)比的過程是將藍(lán)色的合成記錄道上的點(diǎn)對(duì)應(yīng)于紅色地震復(fù)合道上的點(diǎn)點(diǎn)對(duì)應(yīng)于紅色地震復(fù)合道上的點(diǎn). 用鼠標(biāo)用鼠標(biāo)左鍵點(diǎn)中，將四個(gè)同相軸連接左鍵點(diǎn)中，將四個(gè)同相軸連接.如右圖所示點(diǎn)擊如右圖所示點(diǎn)擊Stretch 按鈕按鈕:新子波自動(dòng)成為新子波自動(dòng)成為 “c

44、urrent wavelet”.這意味它被用于在井這意味它被用于在井對(duì)比窗口中合成記錄的制作對(duì)比窗口中合成記錄的制作. 注意到藍(lán)色合成記錄道有微小注意到藍(lán)色合成記錄道有微小的變化的變化.81STRATA 培訓(xùn)點(diǎn)擊點(diǎn)擊 Ok (在在Check Shot Analysis窗口窗口).井對(duì)比類似于井對(duì)比類似于check shot 校正校正. 故出現(xiàn)兩個(gè)新窗口故出現(xiàn)兩個(gè)新窗口, 允許設(shè)允許設(shè)置一些參數(shù)置一些參數(shù). 后面更仔細(xì)談?wù)撨@些參數(shù)后面更仔細(xì)談?wù)撨@些參數(shù).82STRATA 培訓(xùn)點(diǎn)擊點(diǎn)擊Wavelet / Extract Wavelet / Use Well:選擇選擇 wedge_well 點(diǎn)擊點(diǎn)擊

45、Next 使用缺省使用缺省值點(diǎn)擊值點(diǎn)擊 OK (在最后一頁在最后一頁).這時(shí)你會(huì)發(fā)現(xiàn)合成記錄拉伸后這時(shí)你會(huì)發(fā)現(xiàn)合成記錄拉伸后 與地震數(shù)據(jù)擬合的更好與地震數(shù)據(jù)擬合的更好. 現(xiàn)在可以利用井與地震之間的相關(guān)求取更好的子波現(xiàn)在可以利用井與地震之間的相關(guān)求取更好的子波.83STRATA 培訓(xùn)接下來都選擇默認(rèn)即可接下來都選擇默認(rèn)即可這一頁定義合成道在井這一頁定義合成道在井位置處進(jìn)行合成位置處進(jìn)行合成選擇時(shí)窗：選擇時(shí)窗：84STRATA 培訓(xùn)這一頁定義子波提取方法這一頁定義子波提取方法以及基本的子波參數(shù)：以及基本的子波參數(shù)：接下來的課程討論一下這接下來的課程討論一下這些參數(shù)些參數(shù)85STRATA 培訓(xùn)新子波

46、看起來與上一個(gè)很相象新子波看起來與上一個(gè)很相象:但是點(diǎn)擊但是點(diǎn)擊Frequency 標(biāo)簽顯示有標(biāo)簽顯示有一個(gè)非常小的相位校正一個(gè)非常小的相位校正.86STRATA 培訓(xùn)求取了子波，點(diǎn)擊求取了子波，點(diǎn)擊Ok 查看新查看新的對(duì)比結(jié)果的對(duì)比結(jié)果.菜單告知新對(duì)比后的聲波曲線將用名字菜單告知新對(duì)比后的聲波曲線將用名字P-wave_corr. 點(diǎn)擊點(diǎn)擊Ok :點(diǎn)擊點(diǎn)擊File / Exit Window eLog 窗口窗口):87STRATA 培訓(xùn)現(xiàn)在加載地震層位用于初始猜測現(xiàn)在加載地震層位用于初始猜測模型的內(nèi)插。模型的內(nèi)插。在在Horizon Selection, 使用缺使用缺省值名字點(diǎn)擊省值名字點(diǎn)擊O

47、k. 點(diǎn)擊點(diǎn)擊No (在在 dialog 窗口詢問是否顯示層位的平窗口詢問是否顯示層位的平面圖面圖).點(diǎn)擊點(diǎn)擊 Horizon / Pick Horizons:88STRATA 培訓(xùn)在在 STRATA 窗口底部顯示層位拾取的一些選項(xiàng)窗口底部顯示層位拾取的一些選項(xiàng):89STRATA 培訓(xùn)如下圖所示：如下圖所示：注意：一個(gè)選擇參數(shù)注意：一個(gè)選擇參數(shù)“Mode”，我們選擇，我們選擇 “Rubber Band”. 如果我們選擇這項(xiàng)后，必須將如果我們選擇這項(xiàng)后，必須將剖面上目的層點(diǎn)亮，點(diǎn)下左鍵剖面上目的層點(diǎn)亮，點(diǎn)下左鍵，釋放鼠標(biāo)后，釋放鼠標(biāo)后， 層位會(huì)被拾層位會(huì)被拾取出來。取出來。也可以選擇也可以選擇

48、“Right and Left Repeat”. 我們選擇用這種方我們選擇用這種方式。式。90STRATA 培訓(xùn)拾取第一層位拾取第一層位, 置鼠標(biāo)在波置鼠標(biāo)在波峰峰150 ms周圍周圍:點(diǎn)擊左鍵點(diǎn)擊左鍵. 同相軸自動(dòng)左右同相軸自動(dòng)左右拾取拾取:T開始拾取第二層位開始拾取第二層位, 點(diǎn)擊點(diǎn)擊 (在在 Horizon 下拉菜單下拉菜單).設(shè)置拾取模式設(shè)置拾取模式“Left & Right Repeat”.91STRATA 培訓(xùn)使用缺省名字使用缺省名字, Horizon 2. 然后改變?nèi)缓蟾淖?抓獲抓獲Snap 選項(xiàng)為波谷選項(xiàng)為波谷 Trough:點(diǎn)擊楔型頂部的波谷點(diǎn)擊楔型頂部的波谷.用類似

49、方式用類似方式, 拾取楔型底部拾取楔型底部以及在以及在450 ms 周圍的波谷波谷如圖示如圖示.完成后點(diǎn)擊完成后點(diǎn)擊Ok存儲(chǔ)拾取的存儲(chǔ)拾取的層位層位.92STRATA 培訓(xùn)現(xiàn)在可以建立初始猜測模型現(xiàn)在可以建立初始猜測模型, 利用井曲線和四個(gè)拾取的層位利用井曲線和四個(gè)拾取的層位. 點(diǎn)擊點(diǎn)擊 Model / Build/Rebuild a Model (在在STRATA窗口窗口):第一頁設(shè)置新模型的名字第一頁設(shè)置新模型的名字. 利用缺利用缺省參數(shù)省參數(shù),參數(shù)含義后面課程再詳細(xì)參數(shù)含義后面課程再詳細(xì)解釋解釋. 連續(xù)點(diǎn)擊連續(xù)點(diǎn)擊Next 最后點(diǎn)擊最后點(diǎn)擊OK 建立模型建立模型.93STRATA 培訓(xùn)初

50、始猜測模型生成，初始猜測模型生成，。在在wedge 地震數(shù)地震數(shù)據(jù)后面初始猜測阻抗據(jù)后面初始猜測阻抗模型用彩色顯示模型用彩色顯示.94STRATA 培訓(xùn)完成第一種反演完成第一種反演, 點(diǎn)擊點(diǎn)擊 Invert / Bandlimited:在反演菜單上在反演菜單上, 將輸出名字將輸出名字Output Volume File改為改為 “bandlimited ,其它用缺省其它用缺省值值:95STRATA 培訓(xùn)點(diǎn)擊點(diǎn)擊OK 得到反演結(jié)果得到反演結(jié)果:96STRATA 培訓(xùn)現(xiàn)在進(jìn)行第二種反演現(xiàn)在進(jìn)行第二種反演(在初始模型窗口在初始模型窗口). 點(diǎn)擊點(diǎn)擊Invert / Model Based:改變輸出名

51、字改變輸出名字 Output Volume File“model_based”:使用其它缺省獲得基于模型的使用其它缺省獲得基于模型的反演結(jié)果反演結(jié)果Model Based:97STRATA 培訓(xùn)鑒別反演結(jié)果好不好的一個(gè)重鑒別反演結(jié)果好不好的一個(gè)重要標(biāo)準(zhǔn)是要標(biāo)準(zhǔn)是 “error display”,點(diǎn)擊點(diǎn)擊眼睛可以從反演窗口查看此項(xiàng)眼睛可以從反演窗口查看此項(xiàng),如右圖如右圖在參數(shù)菜單選擇在參數(shù)菜單選擇“model_based derived Synthetic Error”, 如右圖所示如右圖所示,點(diǎn)擊點(diǎn)擊Ok.如右圖如右圖,可以看到原始楔形模型與可以看到原始楔形模型與反演結(jié)果生成合成記錄之間的差反

52、演結(jié)果生成合成記錄之間的差異異, 這個(gè)值接近于這個(gè)值接近于0,指示反演結(jié)果指示反演結(jié)果是成功的。是成功的。98STRATA 培訓(xùn)在在Open Window Data Menu 顯示了剛生成的所有數(shù)據(jù)體顯示了剛生成的所有數(shù)據(jù)體:該練習(xí)結(jié)束該練習(xí)結(jié)束. 點(diǎn)擊點(diǎn)擊File / Exit Project (在在 STRATA 軟件任何窗口軟件任何窗口)點(diǎn)擊點(diǎn)擊Yes (在警告診斷窗口在警告診斷窗口)使使STRATA 關(guān)閉所有窗口并保存工區(qū)關(guān)閉所有窗口并保存工區(qū).99STRATA 培訓(xùn)練習(xí)練習(xí) 2, 對(duì)加拿大西部的一塊對(duì)加拿大西部的一塊3D數(shù)據(jù)做反演數(shù)據(jù)做反演. 第一步是讀一口井曲線到已經(jīng)定義好第一步是

53、讀一口井曲線到已經(jīng)定義好的的GEOVIEW 數(shù)據(jù)庫中數(shù)據(jù)庫中.在做這個(gè)之前，先創(chuàng)建一個(gè)數(shù)據(jù)庫在做這個(gè)之前，先創(chuàng)建一個(gè)數(shù)據(jù)庫點(diǎn)擊點(diǎn)擊 Database / New.定義這個(gè)數(shù)據(jù)庫名稱定義這個(gè)數(shù)據(jù)庫名稱 erskine_database.練習(xí)練習(xí) 2: Erskine 3D 初始化初始化100STRATA 培訓(xùn)在在 GEOVIEW Well Explorer 窗口窗口, 點(diǎn)擊點(diǎn)擊Import Data / Logs, Check Shots, Tops, Deviated Geometry from Files:101STRATA 培訓(xùn)改變改變 Destination Well Name 域?yàn)橛驗(yàn)?/p>

54、 “erskine_well”點(diǎn)擊點(diǎn)擊Next :在在 File Import 頁面頁面, 選擇文件選擇文件e“erskine_log.las” 點(diǎn)擊點(diǎn)擊Next .102STRATA 培訓(xùn)顯示有一口井聲波曲線顯示有一口井聲波曲線. 點(diǎn)擊點(diǎn)擊Ok 讀取這根曲線讀取這根曲線.點(diǎn)擊點(diǎn)擊Next (利用缺省位置信息利用缺省位置信息)選取默認(rèn)單位選取默認(rèn)單位103STRATA 培訓(xùn)在井曲線讀進(jìn)在井曲線讀進(jìn)GEOVIEW后后, 點(diǎn)擊點(diǎn)擊erskine_well 然后點(diǎn)擊然后點(diǎn)擊Display Well:erskine 聲波顯示如右圖聲波顯示如右圖.104STRATA 培訓(xùn)在在GEOVIEW 主窗口上點(diǎn)擊

55、主窗口上點(diǎn)擊STRATA啟動(dòng)啟動(dòng) STRATA 軟軟件件.定義工區(qū)為定義工區(qū)為 “erskine_project”選擇選擇 Start New Project:在在STRATA 建立新工區(qū)對(duì)建立新工區(qū)對(duì)erskine 地震數(shù)據(jù)進(jìn)行反演地震數(shù)據(jù)進(jìn)行反演.105STRATA 培訓(xùn)第一步將第一步將3D 地震數(shù)據(jù)讀進(jìn)地震數(shù)據(jù)讀進(jìn)STRATA軟件軟件.點(diǎn)擊點(diǎn)擊Data Manager / Import Data / Open Seismic / From SEG-Y File:選擇文件選擇文件 “ersk3d.sgy”點(diǎn)擊點(diǎn)擊Next :106STRATA 培訓(xùn)同樣做為三維地震數(shù)據(jù)體加載地震數(shù)據(jù)同樣做為

56、三維地震數(shù)據(jù)體加載地震數(shù)據(jù).這個(gè)文件在道頭中沒有這個(gè)文件在道頭中沒有Inline & Xline 號(hào)或號(hào)或 X & Y 坐標(biāo)坐標(biāo). 設(shè)置設(shè)置 primary key to “CDP (Rectangular)”. 這是指數(shù)據(jù)體是由一系列測線組成這是指數(shù)據(jù)體是由一系列測線組成,每條線有相同的每條線有相同的道數(shù)道數(shù).107STRATA 培訓(xùn)改變觀測系統(tǒng)鍵入改變觀測系統(tǒng)鍵入 “155” 如圖顯示如圖顯示. INLINE線將由程序計(jì)算線將由程序計(jì)算.點(diǎn)擊點(diǎn)擊Ok 加載地震數(shù)據(jù)體加載地震數(shù)據(jù)體.點(diǎn)擊點(diǎn)擊Next 兩次選擇缺省直到最后一頁兩次選擇缺省直到最后一頁. 軟件最初假定只軟件最初假

57、定只有單一一條有單一一條INLINE線線.108STRATA 培訓(xùn)在在 Well To Seismic Map 菜單菜單, 如下所示將井設(shè)置在地震如下所示將井設(shè)置在地震數(shù)據(jù)上數(shù)據(jù)上. 在在Plot box 中中wedge_well 井沒有被選用因?yàn)樗疀]有被選用因?yàn)樗粚儆谶@個(gè)數(shù)據(jù)體不屬于這個(gè)數(shù)據(jù)體:在地震顯示中在地震顯示中,鍵入鍵入 “24” 作為作為Inline 號(hào)顯示的號(hào)顯示的結(jié)果如右圖所示結(jié)果如右圖所示:109STRATA 培訓(xùn)對(duì)這個(gè)數(shù)據(jù)體要建立初始猜測模型對(duì)這個(gè)數(shù)據(jù)體要建立初始猜測模型,需要一組層位文件需要一組層位文件.首首先拾取一個(gè)層位然后讀進(jìn)以前拾取的一組層位先拾取一個(gè)層位然后讀

58、進(jìn)以前拾取的一組層位.點(diǎn)擊點(diǎn)擊Horizon / Pick Horizons:使用缺省名字使用缺省名字 “Horizon 1”.點(diǎn)擊點(diǎn)擊Ok.點(diǎn)擊點(diǎn)擊“Yes” (診斷窗口詢問是否診斷窗口詢問是否顯示層位平面圖顯示層位平面圖).110STRATA 培訓(xùn)向以前一樣向以前一樣, 將鼠標(biāo)點(diǎn)在靠近層位處拾取單個(gè)層位將鼠標(biāo)點(diǎn)在靠近層位處拾取單個(gè)層位:平面圖窗口顯示這條平面圖窗口顯示這條INLINE線上拾取的時(shí)線上拾取的時(shí)間間.111STRATA 培訓(xùn)軟件利用這個(gè)人工拾取層作為引導(dǎo)可以自動(dòng)拾取整個(gè)軟件利用這個(gè)人工拾取層作為引導(dǎo)可以自動(dòng)拾取整個(gè)3D數(shù)數(shù)據(jù)體層位據(jù)體層位.點(diǎn)擊點(diǎn)擊Options / Autom

59、atic Picking:點(diǎn)擊點(diǎn)擊Ok (在在Automatic Picking.在層位的平面圖上觀測在層位的平面圖上觀測,最左邊上的幾條最左邊上的幾條INLINE線上可能拾取的層位位置有問題線上可能拾取的層位位置有問題.112STRATA 培訓(xùn)顯示顯示inline 2 測線測線(在在STRATA 窗口窗口).用手動(dòng)方式在同相軸上周圍重新拾取層位用手動(dòng)方式在同相軸上周圍重新拾取層位:我們發(fā)現(xiàn)自動(dòng)拾取的結(jié)果我們發(fā)現(xiàn)自動(dòng)拾取的結(jié)果在一個(gè)小區(qū)域上層位跳出在一個(gè)小區(qū)域上層位跳出一段一段.113STRATA 培訓(xùn)現(xiàn)在再重新自動(dòng)拾取現(xiàn)在再重新自動(dòng)拾取.在自動(dòng)拾取菜單上發(fā)現(xiàn)缺省參數(shù)是清除以前自動(dòng)拾在自動(dòng)拾取

60、菜單上發(fā)現(xiàn)缺省參數(shù)是清除以前自動(dòng)拾取的層位取的層位,而保留人工拾取的層位作為新引導(dǎo)而保留人工拾取的層位作為新引導(dǎo):點(diǎn)擊點(diǎn)擊Ok 得到新層位結(jié)果得到新層位結(jié)果.人工拾取第二條人工拾取第二條INLINE線后模型得到了改進(jìn)去引導(dǎo)自動(dòng)拾取線后模型得到了改進(jìn)去引導(dǎo)自動(dòng)拾取:114STRATA 培訓(xùn)下面我們讀進(jìn)層位下面我們讀進(jìn)層位.首先首先刪除剛拾取的層位刪除剛拾取的層位:層位層位Ok 刪除拾取的層位刪除拾取的層位.點(diǎn)擊點(diǎn)擊Horizon / Import Horizons / From File:選擇五個(gè)文件從選擇五個(gè)文件從 erskine1.pik 到到 erskine5.pik.115STRATA 培訓(xùn)文件類型是缺省文件類型是缺省Ge