自編JASON使用說明

JASON操作流程及實(shí)用初級(jí)手冊(cè)

（invertrace部分-試用版）JASON(invertrace)操作流程、八?4前言一：JASON軟件工作的基礎(chǔ)與關(guān)鍵：井?dāng)?shù)據(jù)編輯、子波提取、合成記錄制作與時(shí)深關(guān)系的調(diào)整。二：工作思路：1：所有的井?dāng)?shù)據(jù)應(yīng)能較好地反映地下地質(zhì)情況的變化。2：首要的是要作好井?dāng)?shù)據(jù)的編輯工作；時(shí)深轉(zhuǎn)換要準(zhǔn)確；合成記錄要盡可能的與地震記錄接近；子波提取要合理；極性與地震數(shù)據(jù)相匹配。3：先做InverTrace,了解油藏的大致分布，并在AI體上重新解釋儲(chǔ)層的頂、底面（因?yàn)槭艿卣鸱直媛实南拗?，地震?shù)據(jù)的波峰或波谷不一定代表油藏的反射）。4：通過交會(huì)圖（crossplot）尋找與反映油藏關(guān)系密切的有關(guān)屬性（如：GR、RES、Impedance、Porosity…?.5:在此基礎(chǔ)上進(jìn)行IverMod的屬性反演，但要有一定數(shù)量的井，InverMod反演才可靠。6：在作好InverTrace、InverMod的基礎(chǔ)上，再應(yīng)用StatMod模塊。第一章數(shù)據(jù)連接與加載一：首先要建立一個(gè)JASON的工作目錄。（例sn4jason）加載的數(shù)據(jù)、中間過程數(shù)據(jù)及最終成果將保存在此目錄下，相當(dāng)于一個(gè)SeisProject。二：數(shù)據(jù)連接與加載。在所建的工作目錄下啟動(dòng)JASON,出現(xiàn)JASON工作主菜單。圖1-1JASON工作主菜單（一）地震工區(qū)相關(guān)數(shù)據(jù)連接主菜單fDatalinksfLandmarkfLandmark98（97,2000）出現(xiàn)次級(jí)菜單。圖1-2地震工區(qū)相關(guān)數(shù)據(jù)連接菜單

l：FilefSeisWorksproject（選擇要研究的地震工區(qū)）

2：Select—import—Seismic/propertydata（選擇要加載的地震數(shù)據(jù)體）

import—Horizons（選擇要研究的地震層位）注：若LandMark地震工區(qū)有現(xiàn)成可用井曲線，可import—wells直接加載。3：Select—Tracegate（選擇要研究的地震工區(qū)范圍）4:Select—Timegate（選擇要研究的地震工區(qū)時(shí)間范圍）5：0ptions—DesireJGWformat~Asfile（8bitinteger 可供選擇的數(shù)據(jù)加載格式6：0ptions—Existingfiles—overwrite（Append 對(duì)已存在文件的處理方式）7：Transport—Import（加載數(shù)據(jù)執(zhí)行）注：在整個(gè)數(shù)據(jù)加載過程中，有些選項(xiàng)是可選的，注意選擇。（二）井曲線數(shù)據(jù)加載井曲線數(shù)據(jù)編輯和加載是非常關(guān)鍵的工作，也是做好反演的前提和基礎(chǔ)，首先要對(duì)收集來的井曲線數(shù)據(jù)進(jìn)行編輯和優(yōu)選，剔除無關(guān)和異常數(shù)據(jù)。一般需要測量深度（Measureddepth）、聲波（P-sonic）、自然伽瑪（Gammaray）以及密度（Density）四種曲線數(shù)據(jù)。井曲線數(shù)據(jù)加載過程如下：主菜單一Datalinks—Wells—Welllogimport調(diào)出井曲線數(shù)據(jù)加載主菜單。圖1-3井曲線數(shù)據(jù)加載主菜單1：Input—Settemplatefile（設(shè)置模板格式文件，選擇一口井即可）。注：在之前進(jìn)行井曲線數(shù)據(jù)編輯時(shí)，最好統(tǒng)一格式，便于批量加載。2：Input—Selectfiles（選擇編輯好的井曲線數(shù)據(jù)文件）。3：Edit—Edittemplate調(diào)出井曲線格式編輯加載工作菜單。（1） ：Startdepth—Changeunit（m）（2） ：Genericlog根據(jù)井曲線數(shù)據(jù)選擇列（Colum）,然后編輯每列的單位（Editheader）.例如：Columl：JGWlogtype選MeasureddepthLogunit選mColum2：JGWlogtype選P-sonicLogunit選us/ftColum3：JGWlogtype選GammarayLogunit選gAPIColum4：JGWlogtype選DensityLogunit選g/cm34:Transport—Importselectedwells（1） Edit井坐標(biāo)（2） Importwells第二章合成記錄制作合成記錄制作是否細(xì)致、精確是關(guān)乎反演結(jié)果是否正確合理的關(guān)鍵因素，也是整個(gè)反演過程中的核心工作。一：要注意的問題（一）子波提取1：子波長度：太長不好，一般為100ms左右為宜：2：制作子波的時(shí)窗（timegate）（1）時(shí)窗不能太小，至少是子波長度的三倍以上。（2）當(dāng)目的層較深，地震資料的信噪比較低時(shí)，提取子波的時(shí)窗最好用淺層部分為好，因?yàn)闇\層資料頻率損失較少；而深層資料頻率損失嚴(yán)重，主頻太低，影響地震反演結(jié)果。（3） 提取子波的時(shí)窗，也可以以目的層為中心，其上取0.05ms,其下取0.05ms（4） 井旁地震道掃描，至少為三道；井向測線的投影方向：規(guī)定為自北順時(shí)針方向轉(zhuǎn)過的角度為測線方位。（5） 第一次的合成記錄用雷克子波，主要目的是把地震解釋的層位與井分層對(duì)齊。后邊的合成記錄用從地震數(shù)據(jù)中提取的子波，不斷調(diào)整合成記錄的波形、振幅與時(shí)間。（6） 子波的比例——至關(guān)重要：因?yàn)樵谒械姆囱菟惴ㄖ?，均以合成記錄代替地震?shù)據(jù)進(jìn)行計(jì)算，因此必須保證合成記錄與地震數(shù)據(jù)的最大匹配，使剩余值更少些。（二）波形調(diào)整：盡量保持其原貌進(jìn)行整體調(diào)整，不宜拉伸或壓縮波形。二：合成記錄制作（一）子波提取主菜單fModelingfWavelets調(diào)出子波提取菜單圖2-1子波提取菜單1：提取雷克子波（Ricker）：Edit—Createsyntheticwavelet調(diào)出子波提取界面。首先Wavelettype選Ricker；frequency[Hz]選25（注意其它參數(shù)的選取）然后Calculate出雷克子波，可通過WaveleteditingQC觀察其波形，然后通過調(diào)整相關(guān)參數(shù)調(diào)整其波形。2：提取其它子波：根據(jù)需要通過Estimate提取其它子波。注：因?yàn)樵谡麄€(gè)反演過程中要生成許多目標(biāo)名稱，因此起名字一定要規(guī)范，便于識(shí)別目標(biāo)。（二）合成記錄制作1：主菜單fAnalysis—WelllogeditingandSeismictie調(diào)出合成記錄制作工作界面。2：Inputfwell（輸入井*.wll）。3：InputfHorizons（輸入層*.hor）。4：Input—Seismicdata（輸入層*.mod）。5：Input—Wavelet（輸入子波*.mtr）。6：Input—Tracegate—Definetracegate（選擇道數(shù)，*可選）。7：主菜單fUtilities—Projectmanagenent—Projectparameters（統(tǒng)一參數(shù)，用Change改，忌直接改，P-sonic（40-150）*可選）。注：以上過程無嚴(yán)格順序。8根據(jù)井曲線數(shù)據(jù)計(jì)算波阻抗曲線（CalculateAI）。注：在合成記錄制作工作界面右側(cè)編輯框（Edit）Calculate選項(xiàng)中根據(jù)實(shí)際情況CalculateAI。9：在合成記錄制作工作界面右側(cè)編輯框（Edit）CreateTDRelation選項(xiàng)中用Sonic給一個(gè)時(shí)深關(guān)系。10：Edit—Scalewavelet…得出合成記錄。注：若正常極性及正反射系數(shù)對(duì)應(yīng)波谷時(shí)，Scalefactor應(yīng)是負(fù)值。11：調(diào)整合成記錄使其盡可能的與地震記錄接近。

12：用右鍵調(diào)出編輯對(duì)話框進(jìn)行加井分層等編輯工作。注：合成記錄制作是整個(gè)反演過程中“人為干預(yù)”最多和耗時(shí)最長的過程，也是影響反演結(jié)果合理正確與否的關(guān)鍵環(huán)節(jié)，其后基本是工作站在進(jìn)行各種運(yùn)算，因此一定要認(rèn)真、仔細(xì)地做好合成記錄1:Seis K.M]EditwellSonic=Wellfogeditingandseismicfie//5.1g(5m0e20)」0jEditFileInputEditOutputDisplayNindouPrintHelpjasoncjp,.s/'.ti'iWI---ISI1.5 BE3sn457 57 56221222223i i i1400;15007I6OO71700t1800-1900^200071:Seis K.M]EditwellSonic=Wellfogeditingandseismicfie//5.1g(5m0e20)」0jEditFileInputEditOutputDisplayNindouPrintHelpjasoncjp,.s/'.ti'iWI---ISI1.5 BE3sn457 57 56221222223i i i1400;15007I6OO71700t1800-1900^20007210072200723002400jZ:Seis,Synt 區(qū)MDensity[g/cnA3JDepth/Tine：1600?1800L2000?2200r2400r2600^■2800MD2365.755nTVD2365-7551VT1866.12■i—EditCunnLCancelUndoRedoCalculateUseVelo->SonicSonic->VeloPatchDensitywithGardnerCalculateAlAutonaticTD-correctSonic(ENTDC)Helllogs---Pseudologs,..Inportlogs...Addlog,*.Renovelog...|-SyntheticsUpdatefiutonaticReady.Nell：sn9.ull<TDfunctionpresent)JWavelet：ricker25_wavelet.ntrJSeisnic：snpcl6.nod3traces$Synthetics!3traces$ResiduHB2=Zoon；<SHIFT>+NB1=Scroll?<CNTRLXSH1FT>+MB1=SelectVerticalaxis.SaveResiduals 11SaveUellSaveHellRsSave圖2-2合成記錄制作工作界面第三章計(jì)算波阻抗體波阻抗體計(jì)算基本上需要四步。注：在這四步中，中間生成的路徑名和文件名一定要起規(guī)范，以免混淆。一：建立模型主菜單fModelingfEarthModel—Modelbuilder（withoutTDC（時(shí)深校正））調(diào)出建立模型的工作菜單。圖3-1建立模型工作菜單l:Input—Time/Depthmode（選時(shí)間和深度模式，選time）2:Input—Horizons（輸入層文件,*.hor）3:Edit—Editframework（編輯模型層格架，從底到頂，生成*.frw文件）注：編輯模型層格架時(shí),要通過QC檢査一下格架的合理性。4：Input—Framework（上步生成的*.frw文件）5：Input—SelectdataforEarthModel—Wells（輸入井,*.wll）6:Input—SelectdataforEarthModel—Wells（輸入井曲線）7：Input—Tracegate（輸入模型的范圍）8：Output—Generate（建立模型）二:模型生成主菜單—Modeling—EarthModel—Modelgenerator調(diào)出模型生成的工作菜單。圖3-2模型生成工作菜單1:Input—Solidmodel（第一步生成的*.mod）2:Input—Dataforseismicmodeling—Seismicdata&Wavelets（Seismicdata選從地震工區(qū)輸入的*.mod,子波要適當(dāng)選?。?:Input—Dataforseismicmodeling—Timegate（用層選取時(shí)間范圍）4:Input—Tracegate（輸入模型生成的范圍）5:Edit—Modelingparameters（編輯模型參數(shù)）注：Timesampleinterval-時(shí)間采樣間隔要和地震工區(qū)輸入的*?mod的時(shí)間采樣間隔保持一致。6:Output—Generate（模型生成）在Selectoutputfiles對(duì)話框中要選PTmpedanceintime三:稀疏約束井反演主菜單—Modeling—Invertrace—Constrainedsparsespike調(diào)出稀疏約束井反演工圖3-3稀疏約束井反演工作菜單1：Input—SeismicdataandWavelets（選輸入地震數(shù)據(jù)體和子波）注：可用通過主菜單—Modeling—Wavelets—Estimate—Estimatewaveletamplitudespecturm提取的振幅譜子波。2：Input—Tracegate（選輸入范圍）3：Input—Timegate（選輸入時(shí)間范圍）在“Usehorizons”對(duì)話框中的Sethorizonsfile選項(xiàng)中要選第一步生成的*_Tinterface.hor文件，然后再選控制層。4：Input—Wells/Userlocations?Wells（選輸入井）5：Input—SelectQCtraces（可選）6：Edit—Edittrend（編輯單口井包絡(luò)函數(shù)）7:Edit—Disableallconstraints（關(guān)掉井約束）8：Edit—Advanced—Waveletscaling（子波估算）若Estimate的“Scalefactor”接近于正負(fù)1,則輸入子波可用，否則要重新提取輸入子波以免影響結(jié)果的正確性。9：Output—Generateresults在Generatethefollowingoutput選項(xiàng)中選P_Impedance，其它根據(jù)實(shí)際需要可選。四:道合并（高頻地震和低頻井約束模型道合并）主菜單—Modeling—Invertrace—Tracemerging調(diào)出道合并工作菜單。圖3-4道合并工作菜單1：Input—Lowfreq.data（選第二步生成的*_Timpedance.mod）2：Input—Bandlimiteddata（選第三步生成的*_impedance.mod）3：Input—Wavelet（輸入子波）4：Input—Tracegate（輸入范圍）5：Edit—Designfilter黃線是井的低頻（〈10Hz）,蘭線是地震的高頻（可給大一些），另外Filterlength要給1024。6：Edit—Edittrend7：Output—Generate選CompleteP_Impedancedata,其它根據(jù)實(shí)際需要可選。至此，稀疏約束井反演（Invertrace）已基本完成，得到了我們要反演的阻抗體。但這并不是想要的最終結(jié)果，我們還需要所反演的目的層（砂層組或砂層）剖面和平面波阻抗特征，因此還需要在得到的反演阻抗體上追出目標(biāo)（目的層），分析其波阻抗屬性得到我們想要的最終結(jié)果。一種方法是在地震體上解釋目的層（砂層組或砂層），然后用解釋層去控制波阻抗體提取波阻抗屬性特征，該種方法難度大，精確度??；另一種方法是在波阻抗體上解釋目的層（砂層組或砂層），然后用解釋層去控制波阻抗體提取波阻抗屬性，該種方法雖然精確度高，但解釋難度大且耗時(shí)長。下面介紹一種簡單可行的目標(biāo)追蹤方法（該種方法適用前提：目的層平行或近似平行輸入層文件的某一層）。第四章目標(biāo)追蹤目標(biāo)追蹤即目的層（砂層組或砂層）自動(dòng)化追蹤解釋基本上分三步來完成。一：三維可視化追蹤主菜單—Analysis—Volumeviewandbodychecking調(diào)出三維可視化追蹤工作界面。

FileInputDisplayBodycheckingPrintHelp恃^^^滋DolljjFileInputDisplayBodycheckingPrintHelp恃^^^滋DolljjCurrentsurvey：sn4Currenthorizon：AllhorizonsNothingpickedReady*圖4-1三維可視化追蹤工作界面1：Input—Selectdata（選擇第三章第四步生成的*.mod）2:Input—Tracegate（輸入范圍）3:Input—Browsetracegate（輸入瀏覽范圍）4:Input—Verticalunit—Time5:Input—Verticalgate（選擇頂?shù)椎目刂茖樱?:Bodychecking—Bodychecking（1） :Type-inrange（根據(jù)巖性與波阻抗關(guān)系確定取值范圍）（2） :Connectivity（調(diào)整相關(guān)系數(shù)的大小來確定分辨程度）（3） :Capture（生成層）（4） :Savehorzions（給頂（Top）、底（Base）和厚度（Thickness））二：建立模型主菜單—Modeling—EarthModel—Modelbuilder（withoutTDC（時(shí)深校正））調(diào)出建立模型的工作菜單。

圖4-2建立模型工作菜單1:Input—Horizons（輸入上步生成的層文件,*.hor）2:Edit—Editframework（編輯模型層格架，從底到頂，生成*.frw文件）3：Input—Framework（上步生成的*.frw文件）4：Input—Tracegate（輸入模型的范圍）5：Output—Generate（建立模型）三：屬性提取主菜單—Analysis—Attributeextraction—Layerattributes調(diào)出層屬性提取菜單。圖4-3層屬性提取菜單1：Input—Solidmodel（輸入上步生成的*.mod文件）2：Input—Seismic/propertydata（選第三章第四步生成的波阻抗體*.mod文件）3：Input—Tracegate4：Edit—（選一個(gè)Top）5：Edit—Attributes（選擇要提取的屬性）6：Output—Generate（生成最終結(jié)果）給“Horzionoutputfilename"—個(gè)名字。第五章目標(biāo)編輯在完成上述工作后，我們可利用反演結(jié)果（波阻抗體）對(duì)