測井解釋軟件Jewel Suite中文操作手冊3

測井解釋軟件JewelSuite中文操作手冊3在第2步選擇porosity做為屬性的名。在這里也可以自定義一個新的名字。第2步的文件夾(folder)的名字填寫tutorial(可以用戶自定義)第3步選擇內(nèi)插的算法，這里選擇distanceweighted方法，缺省的冪指數(shù)為2.上圖中，插值的類型有:內(nèi)插、廣義函數(shù)、模擬、過濾、相等模擬方法，這里僅練習內(nèi)插和模擬兩種方法。第4步，首先選擇參與模擬的井，然后選擇模擬的目標曲線。1.模擬算法內(nèi)插算法:,距離權重內(nèi)插,普通克里金內(nèi)插,協(xié)同克里金內(nèi)插,最臨近域內(nèi)插;模擬算法，這些算法使用3D模擬技術模擬3D網(wǎng)格屬性，需要輸入井曲線數(shù)1據(jù):,順序高斯模擬。,協(xié)同高斯模擬。相模擬，使用相模擬的方法對3D網(wǎng)格的相屬性進行模擬，使用井曲線數(shù)據(jù)進行基于模型的模擬:,基于目標的相模擬;,順序指示模擬，使用井數(shù)據(jù)產(chǎn)生網(wǎng)格屬性體，其中井的數(shù)據(jù)為不連續(xù)的數(shù)據(jù)。離散分布模擬:不需要井曲線,常量分布模擬;,一致性模擬分布;,正態(tài)分布模擬;,曲線正態(tài)分布模擬;,三角剖分分布模擬;,直方圖分布模擬;,基于巖性的屬性充填。濾波器:采用6點濾波器發(fā)對屬性進行濾波，確保濾掉偏離結果的屬性。在第6步，點擊execute按鈕，執(zhí)行屬性模擬處理。打開網(wǎng)格模型，選擇顯示網(wǎng)格的孔隙度屬性，顯示圖例:porosity屬性，打開view-propertyexplorer，在colorsetmin填0.05，在選擇colorsetmax處填0.25，這樣就把孔隙度顏色比例進行了重新設置:回到屬性模擬流程窗口，改變第3步的反距離內(nèi)插的指數(shù)為4，點擊第6步的執(zhí)行，查看不同參數(shù)對內(nèi)插孔隙度值的影響，上圖變?yōu)橄聢D:接下來使用不同的插值方法插值出密度屬性，在第2步改變屬性的名稱為density。在第3步選擇ordinarykriging，type=isotropic，缺省lateraltrend為0.5在第4步，選擇進行模擬的井曲線為密度曲線。點擊第6步的執(zhí)行。把前邊進行密度模擬時的參數(shù)lateraltrend修改為1.5，上圖就變?yōu)橄聢D:2.變差分析前邊使用了插值算法進行了屬性插值。接下來，使用慣續(xù)高斯算法進行孔隙度屬性模擬。首先需要創(chuàng)建變差模型，使用變差模型作為慣續(xù)高斯模擬算法的一項輸入?；氐搅鞒坦芾泶翱冢x擇modelvalidation-modelstatistics菜單在第1步，選擇井類型。在第2步，選擇所有的井(可以點擊鼠標右鍵，選擇“全選”)在第3步選擇孔隙度曲線(porosity)。在第4步，可以根據(jù)需要選擇濾波器，針對需要的井段、數(shù)值區(qū)域進行分選，在這次練習中不選濾波器，對全井段進行變差分析。在第5步，選擇variogram按鈕，彈出一個窗口，可以在窗口里定義變差分析的初始參數(shù)。初始參數(shù)定義如下:點擊create按鈕，產(chǎn)生變差結果。在右側的窗口的總部，缺省的變差類型是lateralomni-directional，可以根據(jù)需要選擇不同的變差類型，下圖:由于這個練習僅有5口井，井數(shù)太少，對于平面的變差分析不太適合。因此采用垂向變差分析。上圖。在上圖的右上角，定義變差類型，這里選擇spherical。可以看到下圖是通過計算數(shù)據(jù)點得到的球形變差。這就是基于數(shù)據(jù)點的計算建立變差函數(shù)的過程。在下圖中使用鼠標拖動曲線，就可以通過手工修改變差曲線?？梢愿鶕?jù)需要調(diào)整有2個控制點或者整條曲線與數(shù)據(jù)點吻合好。如下圖:接下來選擇側向(lateral)檢查側向變差函數(shù)。由于井數(shù)較少的限制，計算出來的側向變差可以能不太合理。這時，就需要使用其他的信息來源建立側向變差函數(shù)，例如:相同區(qū)域、或地質(zhì)條件類似的其他油藏的信息。為2000米，可以看到變差函數(shù)定義右上角對話框的水平變差量圖像隨之改變。注意:如果在選擇，可以選擇3個變差:垂向、主要方向(對平)、次要方向(水平，與主方向垂直)變差。在對話框中可以詳細定義每個方向的參數(shù)。在下方的窗口，可以定義每個方向變差函數(shù)的搜索半徑:得到想要的變差函數(shù)后，把它保存起來，點擊，在彈出的窗口中給個自定義的名字，把變差分析的結果儲存起來，以備進行屬性模擬時使用。3.建立屬性模型回到流程的主菜單，進入ModelStatistics流程窗口，打開統(tǒng)計管理，可以點擊快捷鍵。選擇variograms，在左下角選擇前邊進行變差分析之后存儲的變差函數(shù)的名字。在這里，這個變差函數(shù)圖像是只讀的，不可以編輯。保存工區(qū)時，這個變差函數(shù)也隨之保存，以備以后使用。點擊上邊窗口下方的按鈕，選擇慣續(xù)高斯模擬算法。Statismanager窗口就關閉了，自動打開屬性模擬流程的窗口，可以看到，變差函數(shù)的設置都自動填寫到該界面上了(見第3步)。如下圖。橢球和克里金兩項內(nèi)容也被缺省的填上了，這些的填寫是基于橫向變差的數(shù)據(jù)。在第2步，把property的名字改為porositysimulated，檢查其他的設置沒有問題，點擊execute按鈕，使用慣續(xù)高斯算法進行模擬。在網(wǎng)格中，可以看到這個新模擬出來的屬性。點擊這個屬性，把它顯示在網(wǎng)格中，如果想不顯示等值線，在網(wǎng)格的切片上點擊右鍵，點擊contours就可以不顯示等值線了:點擊切片，按鍵盤的?和?鍵，可以逐層顯示切片。屬性模擬可以產(chǎn)生多個實現(xiàn)，在第3步就是定義進行多少個實現(xiàn)的模擬，種子點的可以輸入，也可以點擊右側的圖標自動產(chǎn)生。上圖是做了3個實現(xiàn)的模擬，在查看其他模擬結果時可以和看切片一樣滾動瀏覽:建議:,在模擬時，創(chuàng)建多個實現(xiàn)，在3d窗口和統(tǒng)計直方圖窗口中進行分析;,使用同一的和隨機的方法創(chuàng)建屬性，注意這兩種方法不要用井曲線;,使用不同的屬性模擬方法模擬指定的模型的layers或zones，在第5步使用濾波器。第二節(jié)創(chuàng)建屬性與編輯使用2個數(shù)據(jù):Tutorial2009-GraphicalPropertyEditorBasics.jewelTutorial2009–GraphicalPropertyEditorCreatingaChannel.jewel1.打開Tutorial2009-GraphicalPropertyEditorBasics.jewel工區(qū):使用快捷按鈕調(diào)整顯示比例和角度:打開網(wǎng)格的屬性，顯示網(wǎng)格的地質(zhì)屬性打開graphicalpropertyeditor流程圖框確認在第2步中已經(jīng)選擇了jewelgrid，第3步中選擇了porosity。當打開propertyeditor時，相關的工具條就顯示出來:點擊工具條上的快捷鍵，在第4步的色標上選擇藍色，在網(wǎng)格上點一些藍點:如果要在網(wǎng)格上畫上未定義數(shù)值的點，如下圖選擇:然后左鍵點擊，就可以在網(wǎng)格上畫未定義數(shù)值的點了。選擇快捷鍵，左鍵點擊編輯定義值，右鍵點擊編輯定義背景:選擇快捷鍵可以選擇與網(wǎng)格中相同的值，用于其他編輯工具的編輯。擦除畫空心矩形畫實心矩形畫線充填相同的鄰區(qū)域充填到邊界:鎖定為當前的編輯狀態(tài)。在第5步可以定義編輯層厚度的方式、參數(shù)。改變northing和easting，觀察編輯的變化:2.編輯河道相打開工區(qū):Tutorial2009-GraphicalPropertyEditorCreatingaChannel.jewel調(diào)整好顯示比例:打開網(wǎng)格的mylithology屬性，右鍵點擊mylithology，選擇editclasses，增加channel巖性，定義好顏色。這個channel可以在graphicalpropertyeditor中進行編輯了。打開graphicalpropertyeditor流程，選擇網(wǎng)格的channel屬性。在處選擇northing，選擇allgridlayers在網(wǎng)格的側面畫一個河道的形態(tài)，在網(wǎng)格屬性里shale處不打勾:同一可以用這種方法定義出多種相類型、形狀:第二章繪制屬性圖練習數(shù)據(jù):\Tutorials2009\PropertyMapping\Tutorial2009-PropertyMapping.jewel第一節(jié)從地質(zhì)到模擬網(wǎng)格的粗化打開練習工區(qū)通過下邊的工具調(diào)整橫、縱比例，放大、縮小、旋轉調(diào)整好顯示角度。打開網(wǎng)格的地質(zhì)屬性-Simulation-GeologicalProperties在流程圖中選擇propertymapping流程:在第1步選擇jewelgrid，在第2步選擇simulation。在第2步時，確認在選擇繪圖方式時選擇upscaling。在第3步選擇孔隙度點擊execute執(zhí)行操作，把孔隙度屬性從地質(zhì)網(wǎng)格中粗化到模擬網(wǎng)格中?？梢钥吹竭@simulation網(wǎng)格中產(chǎn)生了一個新的地質(zhì)屬性:mapped_porosity。選擇simulation網(wǎng)格的mapped_porosity顯示如下:第二節(jié)從模擬到地震的網(wǎng)格細化打開顯示HFBasesmallProperties屬性，這個屬性對是含氣飽和度的25個時間步長的模擬。選擇這個屬性，把時間拉到10。選擇地震數(shù)據(jù):模擬結果和地震數(shù)據(jù)都顯示出來:到propertymapping流程窗口，在第1步選擇simulation，第2步選擇haddofield.segy，確認第2步還選擇了downscaling，在第3步只選擇hfbasesmall_sg屬性。點擊execute，執(zhí)行網(wǎng)格細化操作。這樣就把含氣飽和度屬性細化到地震網(wǎng)格上了。顯示模擬的網(wǎng)格，在這里每一個結果代表一個時間點含氣飽和度的情況:與代表相同的意思。不顯示模擬網(wǎng)格，僅顯示地震體?，F(xiàn)在，已經(jīng)把含氣飽和度數(shù)據(jù)存入地震數(shù)據(jù)體中。在地震體中顯示含氣飽和度屬性，在seismicamplitude處點擊右鍵，彈出菜單，選擇SetAsOverlay這樣就把含氣飽和度屬性疊加顯示在地震數(shù)據(jù)的振幅上了。屬性的網(wǎng)格細化過程需要占用大量的硬件內(nèi)存，在內(nèi)存不夠時，JewelSuite會給出提示。出現(xiàn)這種情況時可以采用下面2種方法解決:,重新裝載小的地震數(shù)據(jù)體，可以對地震數(shù)據(jù)進行重采樣，包括針對地震采樣率進行重采樣、主測線/橫測線道間距重采樣等方法減小地震數(shù)據(jù)體。,調(diào)整流程第3步的時間步長。第三節(jié)從地震到地質(zhì)的網(wǎng)格粗化把地震、地質(zhì)網(wǎng)格同時都顯示出來。在propertymapping流程中，第1步選擇HaddoField.segy，第2步選擇jewelgrid，確認在第3步只選擇了seismicamplitude，點擊execute，執(zhí)行操作。在地震網(wǎng)格中顯示出粗化進來的地震振幅屬性:右鍵在左側圖框中點擊地質(zhì)網(wǎng)格-solidandmesh—solid，顯示如下:右鍵在3維窗口中點擊網(wǎng)格，選擇nodisplayfilter-linear第四節(jié)從網(wǎng)格屬性創(chuàng)建偽測井曲線練習數(shù)據(jù):2009\Tutorials2009\PropertyMapping\Tutorial2009-Pseudologs.jewel打開流程:點擊newview，打開wellcorrelationview選擇地質(zhì)網(wǎng)格、井、網(wǎng)格中的孔隙度，點擊即可完成。第三章網(wǎng)格粗化第一節(jié)規(guī)則網(wǎng)格的粗化使用的數(shù)據(jù):2009\Tutorials2009\ModelUpscaling\Tutorial2009-UpscalingRegularGrid.jew打開工區(qū):打開網(wǎng)格粗化流程在第1步選擇regulargrid，在第二步輸入simulationmodel，點擊showmodeldefinitionview打開網(wǎng)格定義窗口。在這里可以定義layer組的數(shù)量。每一個組有一個或多個layer構成，并且與單個模擬的層相當。這些組將用于垂直網(wǎng)格粗化。在模型定義窗口(modeldefinitionview)點擊reservoir_1，然后點擊modelupscaling窗口中的group定義單個層為模擬的層。按住shifte鍵點擊Reservoir_2和Reservoir_4，選擇這2個層，在這2個層的任意一個上點擊鼠標右鍵，選擇grouplayers。同樣把Reservoir_5創(chuàng)建一個組，Reservoir_6到Reservoir_9創(chuàng)建成一個組，Reservoir_10創(chuàng)建組。點擊第6步的execute，執(zhí)行該操作。在3維圖上顯示simulationgrids，如下圖右側，右鍵點擊網(wǎng)格-屬性-layerid:這樣就把原來的細網(wǎng)格粗化成上面的網(wǎng)格圖例。在流程的第3步，可以定義橫向網(wǎng)格的大小，在north方向給3，在east方向給2，如下圖:點擊執(zhí)行，網(wǎng)格會立刻更新，下圖為更新后同樣顯示網(wǎng)格id屬性圖:做下圖左側圖框里的選項，可以看到組的切片:如果把地質(zhì)網(wǎng)格也顯示出來，可以看到地質(zhì)網(wǎng)格比模擬網(wǎng)格細:回到流程管理的主菜單，選擇simulationpreparation:在第一步選擇simulationmodel。注意，在第4步可以選Transmissibilities也可以選delNorth/East/Thickness，因為在這個工區(qū)里沒有斷層，這就意味著這兩種操作都可應用在這里。選擇North/East/Thickness，在第7步點擊execute。在網(wǎng)格上點擊鼠標右鍵Properties,Thickness,SimPreparation_CellThickness，顯示圖例，可以看到下圖:通過這種方法還可以實現(xiàn)以下目的:,可以看到，通過改變油藏的定義的組來改變網(wǎng)格的垂向比例。,可以使用不同的網(wǎng)格粗化方法創(chuàng)建額外的地質(zhì)屬性;,創(chuàng)建非滲透層。第二節(jié)對Jewel網(wǎng)格的粗化使用的數(shù)據(jù):\Tutorials2009\Upscaling\Tutorial2009–UpscalingJewelGrid.jewel打開工區(qū)進入modelupscaling流程:第1步選擇geologicalmodel，這個網(wǎng)格是64×64×37的。第2步輸入名字simulationmodel。第3步定義橫向網(wǎng)格，這個練習既要做橫向的網(wǎng)格粗化，也要做垂向的網(wǎng)格粗化，首先定義橫向網(wǎng)格粗化。North方向網(wǎng)格給2，east方向網(wǎng)格給3。接下來定義垂向網(wǎng)格。第4步點擊showmodeldefinitionview與地質(zhì)網(wǎng)格相關定義一些組(group)這些組就用于垂向網(wǎng)格粗化。按照下圖，在油藏上點擊鼠標右鍵，選擇group，即可建立一些組，如下圖:下一步，使泥巖層在模擬網(wǎng)格中為不活的層。Ctrl+鼠標左鍵，選中5個泥巖層，右鍵點擊，在菜單中選擇impenetrable使泥巖層為不活的層?；氐饺S窗口。在流程的第5步選擇porosity，選擇缺省的Arithmeticmeanweighted網(wǎng)格粗化方法，后邊是CellVolume屬性。還選擇3個滲透率屬性:kx，ky，kz，同樣適用缺省算法進行網(wǎng)格粗化。在第6步，點擊execute，得到了粗化的網(wǎng)格:32x32x13調(diào)整網(wǎng)格的邊界:打開view-PropertyExplorer，在SolutionExplorer窗口內(nèi)點擊Simulationmodel網(wǎng)格的FullVolume，在PropertyExplorer窗口內(nèi)會出現(xiàn)Boundary項，在的Boundary前邊帶+號，點擊+，打開Boundary如下圖:在StartNorth處把0改成5，在StartEast處把0改成5。在改變StartEast的值后必須按TAB鍵后才能更新。其它方便的方法更新邊界，就是用AnimationControl對話框。按F8鍵，會彈出對話框，選擇3dmodel，可以很方便的調(diào)整邊界。這個方法也可以調(diào)整地質(zhì)網(wǎng)格的邊界和地震體的邊界。調(diào)整之后的網(wǎng)格和地質(zhì)網(wǎng)格的十字交叉線在一起顯示:右鍵點擊Simulationmodel選擇Resetboundaries完成網(wǎng)格邊界的重新設置，恢復到完整狀態(tài)的網(wǎng)格。第三節(jié)模擬前的準備使用上一節(jié)的結果數(shù)據(jù)。打開modelupscaling-simulationpreparation流程。在進行數(shù)值模擬以前，要在這個“模擬前的準備”流程里計算網(wǎng)格的一些屬性數(shù)據(jù)包括:,網(wǎng)格的可傳遞的屬性(TX,TY,TZ);,無相鄰連接關系(網(wǎng)格間有斷層的情況);,網(wǎng)格轉換的適宜性說明。在第1步選擇模擬網(wǎng)格。在流程的第2步，輸入預備模擬的名字:simpreparation1(這個例子里會做2個預備模擬)。由于要進行新的預備模擬，不選第3步。在第4步選擇，因為模擬網(wǎng)格有斷層。選擇算法，這種算法可以盡可能的把小的網(wǎng)格合并成有效網(wǎng)格，所以在第5步的合并小網(wǎng)格的操作可以跳過。在第7步，點擊execute，執(zhí)行該操作。這樣，模擬網(wǎng)格就準備好了，這就意味著所有的可傳遞的屬性以及和相鄰網(wǎng)格不連通性都已經(jīng)計算出來。所有這些屬性都是數(shù)值模擬需要用到的。這些屬性都被分開放在不同的文件夾下。在三維窗口中，右鍵點擊模擬網(wǎng)格模型，Properties-Transmissibility-SimPreparation1_TX，或者直接如下圖左側圖框中選擇:可以把一個模擬網(wǎng)格創(chuàng)建成多個預備模擬網(wǎng)格，在創(chuàng)建時考慮到不同的因素:,不同的滲透率屬性;,不同的IJK簡化算法;,不同的小網(wǎng)格合并方法。接下來要使用不同的IJK簡化算法和小網(wǎng)格合并方法來計算。在第2步選擇新的操作，在這里填上新的名字:“SimPreparation2”在第4步，改變IJK算法為CompactIJK，在第5步選擇Mergesmallcells參數(shù)為10%。在第7步，點擊execute，執(zhí)行操作。可以看到這SolutionExplorer窗口中出現(xiàn)了一個新的屬性文件夾:“SimPreparation2”選擇MergedCells屬性，不選noneadjustedcell項，如下圖:從上圖就可以看到所有的mergedsmallcells網(wǎng)格，這些網(wǎng)格的小網(wǎng)格被合并了。這個模擬模型現(xiàn)在為數(shù)值模擬準備了2個預備網(wǎng)格。對與這個模型還可以進行一下幾方面的練習:,在模型定義窗口中改變組，以便改變垂向網(wǎng)格粗化參數(shù);,改變平面網(wǎng)格尺度，以改變平面網(wǎng)格粗化參數(shù);,從地質(zhì)網(wǎng)格到模擬網(wǎng)格粗化更多屬性，使用不同的粗化方法進行粗化;,改變小網(wǎng)格合并方法，在三維窗口中查看小網(wǎng)格合并情況。第四章鉆井設計使用的數(shù)據(jù):\Tutorials2009\WellPlanning\Tutorial2009-WellPlanning.jewel第一節(jié)直井設計打開工區(qū)進入到井位設計流程:點擊快捷鍵Createverticalwell在網(wǎng)格上的任意一點點擊，就會出來一口新井。在流程的第2步，改變井名為INJECDT1，井類型injector;狀態(tài):水井;條件:油藏;顏色:淺藍。在第3步選擇偏移量。把設計井巷北移動300米，向東移動300米，點擊movewholewelldesign。井就會移動。在第3步，選擇Movecopyofselectedwelldesign就會復制一口新井移動。改變它的名字、類型、狀態(tài)、條件、顏色，這個設計井的這些屬性離開會改變。第二節(jié)斜井設計點擊井位設計快捷鍵工具欄的斜井設計快捷鍵，在網(wǎng)格的任意一點點擊，這樣就創(chuàng)建了一口井。接下來可以加這口井的其他拐點。改變設計井的名字、井類型、狀態(tài)、條件、顏色。打開jewelgrid，只顯示切片。選擇快捷鍵，在切片上點擊任意一點，就在該處增加了設計井的拐點。這樣，就可以參考地質(zhì)屬性網(wǎng)格高孔隙度的層設計井位。選擇任意一個方向的剖面，在剖面上點擊，這樣就在剖面上增加了井的拐點。在第4步選擇設計的拐點，在第5步就顯示出該拐點的深度、坐標，修改這里的深度、坐標后。修改后，在流程的最下邊點擊，安裝修改的點重新調(diào)整井軌跡。使用快捷鍵:,移動設計井拐點:鼠標左鍵點擊設計井拐點，在想把拐點移動到的目標位置的網(wǎng)格上再點擊鼠標左鍵，設計井的拐點就移動過去了。,刪除設計井拐點:點擊右鍵-刪除，即可。第五章計算體積計算的體積報告可以是整個模型的、油藏(zones)、網(wǎng)格斷塊(segments)、油藏斷塊(Compartments)。使用數(shù)據(jù):2009\Tutorials2009\Volumetrics\Tutorial2009–VolumetricsReporting.jewel第一節(jié)創(chuàng)建獨立的區(qū)帶(zone、segments、Compartments)這一節(jié)練習地質(zhì)網(wǎng)格模型或模擬網(wǎng)格模型的油藏(zones)、網(wǎng)格斷塊(segments)、油藏斷塊(Compartments)。對于這個指定的區(qū)域，可以計算它的體積，并且可以在三維窗口中獨立的顯示出來。打開練習工區(qū):打開網(wǎng)格的屬性，可以看到所有用于練習計算體積的的地質(zhì)屬性選擇ZoneId_level_1，基于這個屬性的區(qū)帶已經(jīng)創(chuàng)建。在SolutionExplorer窗口中，右鍵點擊JewelGrid，在create里選擇zones這樣就創(chuàng)建了zone(油藏)，在網(wǎng)格文件夾下可以看到zone項。下圖是顯示2個zone。創(chuàng)建網(wǎng)格斷塊(segments):同樣在SolutionExplorer窗口中，右鍵點擊JewelGrid，在create里選擇segments，就創(chuàng)建了網(wǎng)格斷塊:網(wǎng)格斷塊是單獨的斷塊，與其他油藏是完全分開的。當創(chuàng)建了油藏斷塊(segments)，同時也創(chuàng)建了SegmentId，ZoneId_level_1和SegmentId在屬性計算里非常有用。SolutionExplorer窗口中，右鍵點擊Jewel創(chuàng)建油藏斷塊(Compartments)，在在create里選擇Compartments，就創(chuàng)建了油藏斷塊:Grid，油藏斷塊(Compartments)屬性值屬性計算里會經(jīng)常用到。油藏斷塊(Compartments)的創(chuàng)建是在油藏(zones)和網(wǎng)格斷塊(segments)的基礎上創(chuàng)建出來的。第二節(jié)創(chuàng)建體積報告使用上一節(jié)的結果數(shù)據(jù)進行練習。進行體積計算時，使用下面幾種屬性:Ntg,Porosity,OilSaturation,GasSaturation,andCellVolume.打開VolumetricsReporting流程在第一步選擇網(wǎng)格，在第二步選擇zones。在第二步出現(xiàn)一個窗口，選擇Zone_Level_1屬性。下邊出現(xiàn)很多zones，在這里可以選擇任意一個或多個zones。如果網(wǎng)格內(nèi)沒有zone,segment,orcompartment，體積計算功能會自動創(chuàng)建它們。在第三步選擇網(wǎng)格屬性:CellVolume,NettoGross,Porosity,OilSaturation,GasSaturation。在第五步，需要定義計算體積的內(nèi)容。只有當Nettogross被定義或缺省設置時，才計算NetRockVolume。只有當Porosity被定義或缺省設置時，才計算NetPoreVolume。在第六步定義報告中使用的單位類型。缺省時就選用工區(qū)的單位設置。這里設定油:百萬桶(mmbbl)，氣:百萬立方米(MMm3)在第7步選擇generateprobalility。點擊generate。這樣就計算出了體積報告的結果?，F(xiàn)在，在第二步選擇Segments，或者Compartments創(chuàng)建報告。在彈出的窗口中點擊yes就替換了原來的報告，點擊no就重新創(chuàng)建一個報告。右鍵點擊報告，選擇ExporttoMSWord、excel、剪切板:進行這項操作首先在電腦上要安裝微軟的office第三節(jié)分析多個實現(xiàn)的屬性報JewelSuite提供了創(chuàng)建多個實現(xiàn)的體積報告和在直方圖中分析結果的第功能?？梢允褂枚鄠€實現(xiàn)輸入統(tǒng)計報告，例如孔隙度、滲透率屬性。在這個練習里，就是要用多個實現(xiàn)創(chuàng)建屬性計算報告。首先在三維窗口里查看有30個實現(xiàn)的孔隙度結果。這30個實現(xiàn)是用屬性模擬工具流程創(chuàng)建的，使用的是慣續(xù)高斯算法。SolutionExplorer窗口中右擊點擊孔隙度實現(xiàn)的屬性，選擇showhistogram，在下圖是選擇了多個實現(xiàn)顯示在直方圖中:選擇所有實現(xiàn)的孔隙度屬性都顯示在直方圖中，關掉Showrealizationsseparately。右鍵點擊直方圖，選擇CumulativeDistribution,Normal這就顯示了孔隙度的分布曲線，給出了選擇區(qū)域的10%、50%、90%的分布特征:第四節(jié)使用多個實現(xiàn)屬性創(chuàng)建10%，50%，90%的分布報告接著上一節(jié)，在體積報告中使用多個實現(xiàn)的孔隙度，產(chǎn)生10%，50%，90%的報告。打開三維顯示窗口，在第二步選擇zones，在第3步選擇Porosity處選擇PorosityRealizations。其他的設置不變。在第七步選擇Normal分布類型。點擊Generate創(chuàng)建體積報告。產(chǎn)生的體積報告包含多個實現(xiàn)的屬性，自動識別這些屬性，用于產(chǎn)生10%，50%，90%的報告。所選擇的分布類型將被用于計算10%，50%，90%。假定沒有定義輸入數(shù)據(jù)的分布類型，缺省使用empiric分布類型。生成概率屬性計算后，在JewelSuite的網(wǎng)格中檢查體積的計算結果。這些都在SolutionExplorer中的單獨的文件夾內(nèi)。在SolutionExplorer中，打開VolumetricsOutputProperties文件夾，這里包括每一個實現(xiàn)的體積計算結果。這些結果可以用于以后的體積結果/報告分析中使用。例如直方圖分析。接下來分析幾個zones的NPV值:在中，右鍵點擊GridReservoir_1NPVProbability屬性，選擇ShowHistogram?？梢钥吹揭粋€值，NPV值是基于網(wǎng)格zone1的第一個孔隙度的實現(xiàn)。不選，右鍵點擊實現(xiàn)的窗口，選擇所有的實現(xiàn)。右鍵點擊直方圖柱，選擇CumulativeDistribution,Normal注意:,更新網(wǎng)格屬性后一定要更新體積報告。,還可以創(chuàng)建油氣擴張網(wǎng)格由于生成體積報告中附錄1參數(shù)和計量單位缺省單位設置如下:種類縮寫名量綱缺省單位1-WayTravelTimeOWTTimes2-WayTravelTimeTWTTimesAccelerationAcAcft/s2AccelerationSpeedAcSAcSft/s3AINU-AcousticImpedanceAcSlowAcSlows/mAcousticSlownessAmountOfSubstanceAOSAOSmolAreaAAft2AzimuthAzmtAngledegCompressibilityCompComp1/psiCPUTimeCPUTimeTimemsCurvatureCurvDLSdeg/ftDensityDensDenslb/ft3Depth&ElevationDELengthftDiameterDiamLengthInDipDipLengthdegDistanceDistLengthftEffectiveKhKHEffKHEffM3ElectricalPotentialEPEPVoltElectricalResistanceEREROhmFractionparticlesFracFracppmGammaRayGRGRgAPIGDGDkg/m3GasDensityGasFormationVolumeFactorBgBgrb/scfGasGravityGGG-GasOilRatioGORGORscf/stbGeometricalDistanceGDistLengthftInclinationInclAngledeg種類縮寫名量綱缺省單位InterfacialTensionITensITensdynes/cmInversePressureIPIP1/psiKFactorKFKF1/msKFactor,TimeFKTFKTm/s/msLateralVelocityLVelVelft/sLinearTimeLTLTdaysMagneticFieldMFMFμTMassMMkgODODkg/m3OilDensityOilExpansionFactorOEPNU-OilFormationVolumeFactorBoBoRb/stbOilGasRatioOGROGRstb/mmcfOilGravityOGGravity-PercentagePercRatio%PermeabilityPermPermmDPressurePPpsiPressureGradientPGPGpa/mProductivityIndexPIPIm3/s/PaRatioRatioRatioRatioReflectivityReflNU-ResistivityResResohm.mSaltConcentrationWsDenskg/m3ScScPxScreenSeismicAmplitudeSANU-SeismicTimeSeismicTimeTimemsSeismicVelocitySeismicVelVelm/sSimTimeTimedaySimulationTimeSonicSoSos/ftStrikeStrikeAngledegTemperatureTTFTensionTensTensdynes/cm3TimeTimeTimedayTransmissibilityTransTransrb-cp/day-psi種類縮寫名量綱缺省單位VelocityVelVelft/sViscosityVisViscPVolumeVVbblVolumeExpansionFactorVEFVEFscf/rcfVolumeFluidCumulativeVFCVFCmbblVolumeFluidFlowVFFVFFbbl/dayVolumeGasCumulativeVGCVGCmmcfVolumeGasFlowVGFVGFmcf/dayWDWDkg/m3WaterDensityWaterFormationVolumeFactorBwBwrb/stb附錄2屬性計算公式及變量在這里給出了JewelSuite中屬性計算模塊的比較詳細的公式變量。第一節(jié)基本運算模塊名字創(chuàng)建表達式描述注釋表達式例子Exsinsin()Sinefunction.ArgumentisDegreesorRadianssin($Dip$)(optionincalculatoruserinterface).coscos()Cosinefunction.ArgumentisDegreesorRadianscos($Dip$)(optionincalculatoruserinterface).tantan()Tangentfunction.ArgumentisDegreesorRadianstan($Dip$)(optionincalculatoruserinterface).powpow()Powerfunction.pow($Depth$,3)expexp()Exponentialfunction.exp(2.0)loglog()Commonlogarithm.log(10.0)1modmod()Residueofdivision.Alternativesyntaxispossible.mod($LayerId$,10)$LayerId$mod10$LayerId$%10asinasin()Inversesinefunction.ResultisDegreesorRadiansasin($Dip$)(optionincalculatoruserinterface).acosacos()Inversecosinefunction.ResultisDegreesorRadiansacos($Dip$)(optionincalculatoruserinterface).atanatan()Inversetangentfunction.ResultisDegreesorRadiansatan($Dip$)(optionincalculatoruserinterface).sqrtsqrt()Squarerootfunction.sqrt($Dip$)例如:if($North$mod10=0or$East$mod10=0)then$Porosity$elsenull第二節(jié)距離運算模塊名字創(chuàng)建表達式描述注釋表達式例子ExDistanceDistance(objecttype,DistancetonearestnodeofDistancetoseveral2-6Distance(wellBore)"Objectname")specifiedobject(absoluteorobjectscanbeDistance(objecttype,usingnodesprojectedoncalculated.Names"Objectname",x,y,z)plane).mustbetypedDistance(wellBore,separatedbya"W0,W1,W2")comma.String“all”oremptystringmeansallobjectsofcorrespondingtype.Distance(wellBore,Alsodifferentobject"W0,W1,W2",typescanbeusedtriMesh,"S0")(syntax-Distance(objecttype,"Objectname",objecttype,"ObjectDistance(markerSet,name",…))."M",0,0,1)Possibleobjecttypes:,seismic,grid,simGrid,markerSet,grid2D,triMesh,polyLineSet,wellBoreDistanceSDistanceS(objecttype,ProjectionofvectorSeeDistancefunctionDistanceS(triMesh,7,8"Objectname",x,y,z)betweencurrentlocationcomments"S",0,0,1)andnearestnodeofspecifiedobjectoncorrespondingdirection.Nearestnodeissearchedusingprojectiononplane.DistanceFilteredDistanceFiltered(objectDistancetonearestobjectSeeDistancefunctionDistanceFiltered(grid,type,"Objectname",node,whichhavecomments"JOAGrid","Filterproperty","Filtercorrespondingproperty“LayerId”,“K=12”)propertyClass");value.AbsolutedistanceordistanceusingprojectiononDistanceFiltered(objectplanecanbecalculated.type,"Objectname","Filterproperty","FilterpropertyClass",x,y,z)DistanceSFilteredDistanceSFilteredProjectionofvectorSeeDistancefunction9DistanceSFiltered((objecttype,"Objectbetweencurrentlocationcommentsgrid,"JOAGrid",name","Filterproperty",andnearestnodeof“LayerId”,“K=12”)"FilterpropertyClass",specifiedobjectonx,y,z)correspondingdirection.Onlynodes,whichhavecorrespondingpropertyvalue,areused.Nearestnodeissearchedusingprojectiononplane.ShapeProjectionShapeProjection(“ShapeGeneratepropertywithShapeProjection(“P“,10name”,x,y,z)values,whichshowsif0,0,1)propertylocationliesinsideoroutsideofcorrespondingcylindricalsurface.第三節(jié)濾波運算模塊名字創(chuàng)建表達式描述注釋表達式例子Exsmoothsmooth(property,skipnullPropertysmoothing.smooth($Kx$,1,11,12Skipnullflag:flag,center,neighbors)0.125,0.125)1–true,smooth(property,skipnull0–false.flag,center,top,base)smooth(property,skipnullflag,center,top,base,north,south,east,west)filter2Dfilter2D(property,filtertype,Thisfunctionisfilter2D($Depth$,northdimension,eastmeantonlyfor2DIDWFilter,3,3)dimension)grid.Possiblefiltertypes:,MeanFilter,GaussianFilter,MedianFilter,IDWFilter第四節(jié)內(nèi)插運算模塊名字創(chuàng)建表達式描述注釋表達式例子Exlineline(,x1,y1,x2,y2,…)InterpolatedvalueTablecanbepreparedline($Ahd$,0,0,0,13,14usinglinearinStatisticManager.1000)line(,“Tablename”)interpolation.Whenapreparedtableline(,,“Tablename”)isused(byname),theline($Ahd$,“Table1D”)dimension,1Dor2D,willbechosenaccordingtothetabledimension.line($N$,$E$,“Table2D”)splinespline(,x1,y1,x2,y2,…)InterpolatedvalueSeelinefunctionspline($Ahd$,0,0,0,13,14usingcubicsplinecomments1000)spline(,“Tablename”)interpolation.spline(,,“Tablename”)spline($Ahd$,“Table1D”)spline($N$,$E$,“Table2D”)splineNsplineN(,x1,y1,x2,y2,…)InterpolatedvalueSeelinefunctionsplineN($Ahd$,0,1000,13,14usingcubicsplinecomments500,0,1000,1000)splineN(,“Tablename”)interpolation(derivativeinnodepointsaresplineN(,,“Tablename”)equaltozero).splineN($Ahd$,“Table1D”)splineN($N$,$E$,“Table2D”)line2Dline2D(,,numX,numY,Interpolatedvalueline2D($X$,$Y$,2,2,x1,…,y1,….,z1,…)usinglinear0,1,0,1,10,20,20,40)interpolation(2Dcase).spline2Dspline2D(,,numX,numY,Interpolatedvaluespline2D($X$,$Y$,2,x1,…,y1,…,z1,…)usingcubicspline2,0,1,0,1,10,20,20,interpolation(2Dcase).40)splineN2DsplineN2D(,,numX,numY,InterpolatedvaluesplineN2D($X$,x1,…,y1,…,z1,…)usingcubicspline$Y$,2,2,0,1,0,1,10,interpolation(derivative20,20,40)innodepointsareequaltozero,2Dcase).第五節(jié)數(shù)學計算模塊名字創(chuàng)建表達式描述注釋表達式例子Exsinhsinh(expression)Hyperbolicsine.sinh($P$)coshcosh(expression)Hyperboliccosine.cosh($P$)tanhtanh(expression)Hyperbolictangent.tanh($P$)absabs(expression)Absolutevalue.abs($P$)clipclip(expression,min,max)Clipvalueusingdefinedclip($P$,0,1)15boundaries.minmin(expression1,Minimumvalue.min($A$,$B$,$C$,0.5)16expression2,...)maxmax(expression1,Maximumvalue.max($A$,$B$,$C$)expression2,...)meanmean(expression1,Meanvalue.mean($A$,$B$,$C$)17expression2,...)medianmedian(expression1,Medianvalue.median($A$[0],$A$[1])expression2,...)stddevstddev(expression1,Standarddeviation.stddev($A$,$B$,$C$)expression2,...)第六節(jié)屬性提取模塊名字創(chuàng)建表達式描述注釋表達式例子ExGridGrid("Gridname","GridpropertyGridSeveralobjects18Grid("MyGrid","Porosity")name")propertycanbeused.valueinGrid("Gridname","GridpropertynodeGrid("MyGrid","Porosity",0,0,name",x,y,z)nearestto1)currentlocation(absoluteorusingprojectiononplane).SeismicSeismic("Seismicgridname",SeismicSeveralobjectsSeismic("MySeismic",19"Propertyname")volumecanbeused."Amplitude")propertySeismic("Seismicgridname",valuein"Propertyname",x,y,z)nodeSeismic("MySeismic",nearestto"Amplitude",0,0,1)currentlocation(absoluteorusingprojectiononplane).SimGridSimGrid("Simulationgridname",SimulationSeveralobjectsSimGrid("MyGrid","Porosity")"Propertyname")pertyvalueSimGrid("Simulationgridname",SimGrid("MyGrid","Porosity","Propertyname",x,y,z)nearestto0,0,1)currentlocation(absoluteorusingprojectiononplane).Grid2DGrid2D("Gridname","Property2DgridSeveralobjectsGrid2D("MyGrid2D","Depth")name")propertycanbeused.Grid2D("MyGrid2D","Depth",value0,0,1)Grid2D("Gridname","Propertynearesttoname",x,y,z)currentlocation(absoluteorusingprojectiononplane).TriMeshTriMesh("Tri-meshname",Tri-meshSeveralobjectsTriMesh("S1,S2","Depth")"Propertyname")propertycanbeused.valueTriMesh("Tri-meshname",nearesttoTriMesh("S","Depth",0,0,1)"Propertyname",x,y,z)currentlocation(absoluteorusingprojectiononplane).PolyLineSetPolyLineSet("PolyLineSetPolylinesetSeveralobjectsPolyLineSet("P","Depth")name","Propertyname")propertycanbeused.valuePolyLineSet("PolyLineSetnearesttoPolyLineSet("all","Depth",0,name","Propertyname",x,y,z)current0,1)location(absoluteorusingprojectiononplane).WellBoreWellBore("Wellborename","LogLogsvalueSeveralobjectsWellBore("all”,"Porosity")name")nearesttocanbeused.currentlocationWellBore("Wellborename","LogWellBore("W1,W2","Porosity")name",x,y,z)(absoluteorusingprojectiononplane).MarkerSetMarkerSet("Markersetname",MarkersetSeveralobjects20MarkerSet("M","Depth")"Propertyname")propertycanbeused.valueMarkerSet("Markersetname",nearesttoMarkerSet("M","Depth",0,0,"Propertyname",x,y,z)current1)location(absoluteorusingprojectiononplane).PropertyExtractionPropertyExtraction(objecttype,GeneralSeveralobjectsPropertyExtraction(wellBore,"Objectname","Propertyname")function.anddifferent"W1,W2","Porosity")SpecifiedobjecttypescanPropertyExtraction(objecttype,objectbeused.Seealso"Objectname","Propertyname",propertycommentstoPropertyExtraction(surface,x,y,z)valueDistancefunction."S","Depth")nearesttocurrentlocation(absoluteorPropertyExtraction(markerSet,using"M","Depth”,0,0,1)projectiononplane).PropertyExtractionNotNullPropertyExtractionNotNull(objectSpecifiedSimilarasPropertyExtraction21type,"Objectname","PropertyobjectPropertyExtractionNotNull(grid2D,"S","Test")name")propertyfunction,butonlyvaluedefined(notnull)PropertyExtractionNotNull(objectnearesttopropertyvaluestype,"Objectname","Propertycurrentaretakenintoname",x,y,z)location.account.Onlydefinedvalues(notnull)aretakenintoaccount.PropertyExtractionFilteredPropertyExtractionFiltered(objectSpecifiedSimilarasPropertyExtractionFiltered(grid,22type,"Objectname","PropertyobjectPropertyExtraction"Grid","CellCenterDepth",name","Filterproperty","Filterpropertyfunction,butonly"LayerId","K=0",0,0,1)propertyclass")valuepointsthatnearesttocorrespondstoPropertyExtractionFilteredcurrentspecifiedfilter(objecttype,"Objectname",location.class,aretaken"Propertyname","FilterOperty","Filterpropertyclass",definedx,y,z)values(notnull)aretakenintoaccount.Onlypointsspecifiedbyfilterinformationaretakenintoaccount.附表:可加載的數(shù)據(jù)類型表數(shù)據(jù)種類格式描述輸入輸出位圖JOAImageFiles(*.joaimage)可可Imagefiles(*.bmp,*.jpg,*.jpeg,*.gif,*.png)可否2D網(wǎng)格JOA2DGridFiles(*.joa2dgrid)可可ASCII2DGridFiles(*.grid2d.ascii)可可CharismaHorizonGridFiles(*.*)可否CPS‐3GridFiles(*.*)可可Geogr