CMG軟件STARS模塊操作手冊_第1頁
CMG軟件STARS模塊操作手冊_第2頁
CMG軟件STARS模塊操作手冊_第3頁
CMG軟件STARS模塊操作手冊_第4頁
CMG軟件STARS模塊操作手冊_第5頁
已閱讀5頁,還剩106頁未讀 繼續免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

關鍵字:Outprn,wprn,outsrf,gridcart,wsrf,wprn,vamod,vatypeDualpor,difrac,con,mod,djfrac,equalsi,por,matrix,permicpor,prpor,rockcp,thconr,thconw,thcono,thcong,hlosspropcmm,pcrit,tcrit,avg,bvg,Cartesian:直角坐標;defaulted:缺省;文件后綴:dat,lst一、前言概述STARS是三維、四相、多組分、熱采、蒸汽添加劑的模擬器。網格系統:直角坐標、變深度/變厚度。分散組分一種組分穩定地分散在各相內,如聚合物、凝膠、微粒、乳化、泡沫。這些概念與靈活的組分物性軟件包配合使用,如吸附、堵塞、非線性粘度、分散、非平衡的變量傳遞。特別是,泡沫有兩種模擬方法:A)機理的,可模擬泡沫的產生、傳播、破裂。B)經驗的,但十分實用,如歷史擬合。天然裂縫油藏,有四個模擬天然裂縫油內的流動:雙孔隙的DP,雙滲的DK,三維多層(從外向內)的MINC,垂向細分多層的VR。過程/機理不同,使用不同的模型。基本近似:裂縫與巖塊二種介質。裂縫PV和儲量小,是主要的流動通道,巖塊恰好相反。自適應隱式STARS有全隱和自適應式兩種。對大部分問題,僅有少數網格塊要求全隱,絕大多數網格塊要求顯式。自適應隱式適合于錐進問題和非均質嚴重的層狀油藏及高滲薄層。使用的△與全隱式的接近。因而節省CPU時間1/3~1/2,自適應的選取隨時間、空間變化,因而是動態。全隱式的井用全隱式的方法求解PW和井塊的變量。對于高度的層狀油藏的多層完井,這種方法消除了收斂性的困難。井的控制、約束功能十分強,定速率/流壓、經濟限制、修井等。矩陣的解法解法稱為AIMSOL,不完全的高斯分解作預處理,GMRES加速。AIMSOL是針對自適雅可比矩陣。通常使用解法的缺省參數即可。局部的CYZM網格加密有兩種用途:近井加密,如錐進和水平井,靜態裂縫的模擬。用戶指明要細分網格的區域。塊間連接、傳導使用解法的缺省參數即可。靈活的網格系統有直角坐標、圓柱坐標、變深度/變厚度坐標,無論二維和三維均可以使用這三種坐標系。9、水層模型兩種:小水層,且形狀邊界清楚,用水區塊模擬水層。大的,甚至無限大的水層用半解析的水層模擬。單位三種:公制、油田、實驗單位。支持用FORTRAN寫成,可在以下機器上成功運行:IBM、CDC、CRAY、HONEYWELL、DEC、PRIME、OULD、UNISYS、APILLO、SUN、HP、FRS。繪圖軟件名為RESULTS,使用SF2文件系統繪圖輸出。它也可以作部分前處理,如形成網格以及φKSHD等值圖的數字化輸入。使用一個轉換程序可將SR2文件轉換成SRF文件。離散井筒當其水平井筒較長時,由于瞬變性、粘性壓力降、不同的多相流方式,因而產生了沿井筒入、產出的非均質性。STARS采用離散井筒,且井筒塊與油藏聯立求解可解決這些問題。井筒流動為顯式,可選取多種流動關系式。地質力學模型包括:塑性變形、剪切膨脹、吞吐的壓實驅動、高壓注入產生裂縫、近井地層的破裂與出砂。由三個子模型構造的地質力學模型可以處理以上的問題。地質力學模型與油藏模型的連接為顯式的。可控體積有限元網格(CVFE)將求解區域離散成許多三角形的單元。細節看:“RESERVOIRSIMULATIONWITHACONTROL—VOLUMEFINITE—ELEMENTMETHOD”,SPERE(AUG1922)349—35F。包含CVEF網格的面可以是IJ、IK或JK面。還可以命名用混合網格,近進為錐進式網格,遠井為粗大的CVRF網格。二、指導前言1、初始運算:輸入數據文件為*.dat,生成的輸出文件有:打印輸出文件*.out,SR2的指標文件*irf和SR·的主要結果文件*mrf。重啟運行:輸入文件有*.dat、*.*irf、*mrf,輸出文件有*.out、*irf和*mrf。輸入數據的分組有7組輸入數據。按可下順序輸入:I/O控制,油藏描述,組分性質,初始條件,數值方法模型,井數據與循環數據。個別的有例外,如I/O的關鍵字也可以出現在循環數據內。個別組內的關鍵字出現也有順序限制。圓括號圓括號內的項是可選的。斜線表示任選其一,如*RANGEHCK(*ON/*OFF)。花括號內裝一個數表,如{SWKRWKROW(PCOW)}飽和度函數表。網格數組如ARRAA:*POR或*POR{GRID}組分物性如*KEYWORD{NCOMP}*KEYWORD{NUMY}*KEYWORD{NUMX}冒號表示數據的范圍,如il(:i2),1≤il<i≤NX。組分的指定‘comp-des’組分號,或組分名。標題卡*TITLE1(至多40字符)’SimulationRun=1-1989-01-23’*TITLE2(至多80字符)’DualporosityproblemusingtheMINCoption’*TITLE3(至多80字符)‘ThisisA12X12X10Cartesiangridsystem’標題出現在輸出的打印或繪圖文件內。2、重啟運行重啟文件二進制文件,包括初始數據與基本變量,按一定的頻率寫入,寫重啟文件是可選的。重啟的必要性敏感性計算,或歷史擬合;改變井的說明;在作大運算之前,稱作一個短的運算看正確性;從一個運算的尾部分出多個不同的運算;重啟的含義不需每次從頭開始計算,僅需從中間某點開始就可以了。如何作重啟用初始運算重啟文件,使用*WRST,*WRST的后面跟寫重啟的時間步數,如*WRST10。重啟運算以*RESTARTn開始,如想從上次運算的最后一個時步開始重啟,那么置n為空,即*RESTAT,否則用*RESTARTn,n表示時間步數。在重啟運算時允許修初始文件的數據!把控制輸出打印文件的內容輸入數據文件內,不在I/O節內使用*WPRN和*OUTPRN,如果想改變二者的參數,它們還可以再次出現循環數據節內。*WPRN表示寫網格數據,井數據,解法數據的頻率,如果不寫井與網格的數據,那么表為:*WPRN*WELL0*WPRN*GRID0這樣可以大大減小文件的長度,節省存儲單元。如果以上兩行不出現,那么每個時步均要輸出井、網格塊的數據,文件迅速增大。*OUTPUT后跟輸出的項目,如哪些井、網格、油藏的數據要輸出。如要輸出全部井的分井層指標應使用*OUTPUT*WELL*ALL如只輸出全部井的一行摘要,應使用*OUTPUT*WELL*BRIEF。控制模擬結果文件的內容使用*OUTSRF控制模擬結果文件SR2的內容。*OUTSRF出現在I/O節,因參數的改變也可再次出現在循環節內。如果網格與井的數據不出現在輸出打印文件內,那么頻率數為0。*OUTSRF規定了哪些井數、網格數據、油藏數據要打印。可以對井數據、網格數據分別使用變量表。3、描述網格系統必須的:*GRID、*DI、*DJ、*DK。可選的:*DEPTH、*DTOP、*DIP。這些關鍵字位于油藏描述節內,在*POR之前出現。三種類型的網格系統,如*GRID*CART10106*GRID*CVRF10106*GRID*RADIAL101015網格尺寸用數組形式輸入,如:*GRID*CART101012*DI*CON100.0*DJ*CON100.0*DK*CON100.0DKKVAR25.0250.0340.075.03*402*50如果*KDIR*UP,那DK的數據是從底部往上部排列。4、給定空網格塊表示空塊有兩種方法:*NULL和*VAMOD,二者應出現在油藏描述節。在*NULL后跟一個數組,如果該塊定空塊則數為0,活塊的數為1。如*NULL*NULL1:101:101:311:41:31:30如果*NULL不出現,那么全部塊為活塊。5、細分網格用*REFINE定義細分的區域,它出現在油藏描述節,且在*NULL和*POR之前。它的后面還要跟細分的要求。如網格10×10×3的基本塊(1、1、3)細分為X方向2塊,Y方向3塊,Z方向2塊,表為:*REFINE113into232**每個方向細分最大塊數為4對雙孔問題,不允許網格的細分。變深/變厚的網格可以細分,但細分的塊與基本塊的厚度相同,如果要求細分塊的屬性(如φKSHD)與基本塊的物性不同,那么應使用*RG。6、時間步長太小用*OUTPRN*ITER*NEWTON輸出內迭代(*ITER)、外迭代(*NEWTON)的信息,對分析太小的△T+分益收斂失敗原因有:內收斂迭代失敗;牛頓迭代失敗,導致切小△T;物質守恒誤差過大。如果在輸出文件內多處出現“iterationroutinefailstoconverge”RGE”,那么試著采取如下的補救措施:縮小△T,辦法是縮小*DTMAX或減小每步要求的變化*NORM*PRESS*NORM*SATUR。增加迭代次數,即放大*ITERMAX用*SDEGREE增大因式分解的級,但會增大存儲量。如果對第一或第二次牛頓迭代的內迭代收斂失敗,但至少在最后一次又收斂,那么問題并不嚴重。牛頓迭代收斂失敗導致△T致的切小,失敗的原因是基本變量的最大變化超過了要求量的2倍。可能會發生非物理的值,如負的壓力或負的飽和度,或者超過了給定的最大的牛頓迭代次數。如果問題是由最大變化造成的,如不頻繁發生問題就不大。如果出現頻繁切小△T,那么應采取如下的補救辦法:檢查KR和PVT曲線的非線性,曲線應光滑。檢查網格步長,以及網格塊的物性數據的合理性。檢查井的約束,生產井給MINPw,注入井給MAXPw。牛頓迭代次數大造成不收斂,那么增加*NEWTONCYC。如果解振蕩,如Sg一會兒〈0,一會〉0,好么采取1和4的措施。置全油藏或油藏的某個部分全隱式解法。尤其變化大的區域。這種情況的例子有:1)沒有使用重向平衡的初使化。這可能導致初期有大的P、S變化,即使全部井是關閉的。這種情況發生應用全隱式。2)有氣頂。在氣頂的底層設置隱式。3)有高K的網格或網格處。造成P變化小,S變化大。置這此塊為全隱式。壓力收斂客限置為0.1Kpa。如果收劍限比每步變量的實際變化大的多,那么會造成物質守恒誤差。容限比變化小一個數量級。使用*MATBLTOL改變模型的靈敏度。在大部分情況下,要求變化(*NORM*PRESS和*SATUR)容限(*CONVERGE*PRESS和*COMVERGE*SATUR)的缺省值是合適的。但對水平井與裂縫問題,應用小值代替缺省值,錐進問題同樣如此。7、定義井輸入順序為:*WELL(必須)PRODCER/INJECJOR/*SHUTIN/*OPEN(必須)*INCOMP(注油/報時必須)OPERATE(至少需要一個OPERATE)MONITOR(可選,不能與非權的注入井連用,后接*PERFV*GEO)PERF/PERFV/PERRG(必須)[全部處于井數據節內。]定義井類型有四類:*PRODUCER/*INJECTOR/*SAUTIN/*OPEN如WELL1’MEDRIVERP1’VERT11;WELL2’MEDRIVERP2’VERT1515;WELL3’MEDRIVERP1’VERT55;WELL4’MEDRIVERP21’VERT1010*PRODRCER1:2*INJECOR3:4*MOBWEIGHT*PERFV1:21:31.65*PERFV3:42:31.878、怎樣關井與開井關井先用*OPERATE、*MONITOR、*PERFV定義一口井;而后在任何時間*TIME或*DATE就可在井,甚至剛一定義完就可以關井。開井;關閉的井可以在任何時候再一次打開。臂如,將一口吞吐井定義成一注、采兩口井,井別的開關轉換如下:wellI’Ingjector1’injectormobweight1’operatebhp1000tinjw450qual7perf1**Ijkwi11188well2’probucer1’probucer2operateLiquid1000perf2**Ijkwi11188shutin2time10shutin1time17open2time40stop9、操作和監控約束OPERATE和MONITOR表示對井的約束,至少需要一個*OPERATE,*MONITOR可選。每口井的第一個*OPERATE,稱為基本的操作約束。如果違返了監控的操作約束,且帶有*CONT,那么監控約束就變為操作約束。如果有多個操作約違反,那么采取的行動是:最嚴歷的是*STOP或*SHUTIN,輕微的是*CONT,(1)生產井生產井的操作:定速率約束(基本操作約束),或定最小的井底流壓。如果是油井就定油率或涂率,氣井定氣率。往后的約束可以是最小產率。如井1有兩個操作約束:*OPERATE*MAX*OIL12000.0*CONT*OPERATE*MIN*NHP1500.0*CONT(2)注入井注入井的操作:定最大的注入速率(基本的操作約束),或定最大的井底流壓約束。注氣井、注水井分別選氣率、水率作為約束。如*INJECTOR2*OPERATE*MAX*WATER1000.0*STOP*OPERATE*MAX*BHP2250.0*STOP(3)監控約束我們極力主張監控生產井的GOR或含水,如:*OPERATE*MAX*OIL1200.0*CONT*OPERATE*MIN*BHP2500.0*CONT*MONIOR*GOR15000.0*STOP10、給定井的指數為了給定井的指數,要使用*GEOMTRY/*PERF/*PERFV/*PERFRG。*PERF/*PERFV/*PERFRG必要,僅對*QEOMETRY可選,僅流度加權的生產井或注入井才用它。為了內部計算井的指數,*GEMOETRY必須后跟多個參數。完井關健字*PERF/*PERFV/*ERFRG必須后跟完井的位置和已外部算好的井指數。如果*GEOMETRY使用,那么必須后跟一個完井的關健字。*GEO與完井關鍵字一同使用。*PERF最適合用于水平井或斜井,當然也可用于垂直井。它的格式為*WELL1’12—09—18—56’*PERF1**ifjfkfwisetn112:41.241*WELL1’12—09—18—56’**radgeofacwfracskin*GEOMETRY*K37524881.00.0*PERF*GEO1**ifjfkfffsetn1·:411如果*VERT位于*WELL的后面,那么就應使用*PERFV。僅在*PERFV內給定完井段K1:K2。如:WELL2*VERT22*PERFV1**kfwisetn2:41.561GEOMETRY應與*PERFGEO/PERFV*GEO/*PERFRG*GEO聯合使用!基本塊(1,1,3)細分為(3,3,2)個小塊,井處于基本塊(1,1,3)和細分塊(2,2,1:2)內,那么*REFINE332*RANGE113*WELL1*PERFRG1**ifjfkfirjrkrwisetn11、水平井模擬水平井有兩種方法:方法1:注入井作為線源,生產井作為線匯。略去了井筒的摩阻壓力降和流體的舉升效應,(holdepeffects)。對有回流(backflow)的問題常給錯誤的結果,少量輕微的回流沿且不嚴重。K的平面差異可引起水平井筒的回流。*WELL*INJECTORor*PRODUCER*OPERATE*GEOMETRY*PERFGEOGEOMETRY和*PERF相當于輸入了井的生產能力。運行以上的數據組,觀察計算的產率,如果結果不對,那么就去掉*GEO,用*PERF直接輸入用戶計算的井指數。方法2:使用井筒的離散模型。這個模型既可用于水平井,又可用于垂井。它最適合摩阻壓力降、舉升效應顯著的那些問題。要引用*WELLBORE,井筒作為塊的次生孔隙度,因而必須指定相應的壓縮系數、巖石類型、Kr表。停止一個模擬運算*STOP出現在數據文件的最后,只能出現一次,停放的時間或日期與它連用,如*DATE19980908*STOP12、建井數據的概要數據順序如下:*TIME/*DATE是必要的用*DTWELL定義階段的初始△t井別*PROUCER/*INJECTOR*OPERATE(MONITOR)*GEOMETRPERF/PERFV/*PERFRG*SHUTIN*OPEN可以使用不同的關鍵字以下是井約束隨時間的改變定義新井:在調整現有井的參數之前,使用步(1)(3)(4)(5)(6)。可以用*ATER改變任何基本操作約束,且*ALTER與*ALTER與*TIME/*DATE同時使用。可以調動I/O控制的參數,以及傳導率乘子。在以后井的改變中,也可使用*DTWELL和*DTMAX關鍵字數據輸入系統13、關鍵字數據輸入系統關鍵字系統前言格式:關鍵字、數據什。如*MAXERROR10,如果關鍵字不出現,就用缺省值。字符集注釋的字符、引號內的字符不檢查正誤,只原樣輸出打印。CGM的字符集由大小寫的字母、數字0——9、關鍵字指示符,算術運算符=、+、—、/構成。關鍵字指示符以下四種情況均有效:*POR/POR/*por/por。**字符串,表示注釋行。可以用COMMENT重新定義新定義新的注釋符。關鍵字順序關鍵字分為若干組,不同組的關鍵字按要求順序輸入。某組關鍵字內的關鍵字可以按任意順序輸入,除非有特別的說明。LIST、MOLIST、*INCLUDE、*COMENT、*TRADNSLATE、*RANGECHCK可以出殃在數據文件的任何位置。某些如*WRST可以出現我有I/O節,也可出現在井數據節內。關鍵字分為為必須、可選的、條件必須的。可選的具有缺省值,如關鍵字不出現就用缺值。串連關鍵字如必要,一行內可以包含多個關鍵字。大小寫關鍵字、字母數字串可以大寫、小寫或任意的混合,文件名按操作系統的要求書寫,如IBM要求大寫。行的長度最大長度130字符,以后的字符被略去。如需變更,可以在源程序文件KWCOML。INC內設變參數MDLINE,并重新編譯源程序。分隔符關鍵字、數據、字符串用空格、逗號或空行分隔開。字串字串放入引號內,如‘5-35-48-W5’。以下情況是允許的‘THISISTHE“RIGHT”WAY。自己定義關鍵字如需要,用戶可以用*TRANSLATE定義新的關鍵字的對應,等價關系。數值以自由格式輸入數據,實數可以不帶小數點、指數用E、e、D、d表示,一個數值內不能有空格。如果要求輸入整數,而實際輸入帶有小數部分,那么打印出錯誤信息,程序停止運行。以下的實數均為有效:25.040,-3,1.23E+02,0.02D-4,34.e02,+2.3,+.3,-..3以下為無效:34.E0234.E0234.E02有空格,指數帶小數。重復計數。23,。23,。23,。41,。27與3*23,。41,。27等價有效整數的范圍如若干整數是相鄰連接變化的,那么可以用,代替TO,如123456101112等價1:4610:12表格**Swkrwkow(pcow)0.20.01.045.00.28990.0220.676919.030.37780.0180.415310.070.46670.0610.21784.090.55560.1430.08351.800.67820.2890.0123.500.75610.4500.0.101.00001.0000.0.0使用表格的輸入選擇*INT,可將來源不同的表格自動合并。錯誤與警告信息數據文件回應打印輸出時,帶有警告或錯誤的信息。如果沒有使用*NOLIST,那么即使有錯誤也不打印錯誤信息。因此,僅當數據文件完全正確以后,才建議使用*NOLIST。*COMMENT注釋(可選)格式:*COMMENTTEXT如:*MAXERROR14**Changemaximumnuberoferrors空行(可選)為使輸入的數據文件便于閱讀,可以加入空行。*RANGECHECK數據的范圍檢查(可選)格式*RANGECHECK(*ON/*OFF)ON要檢查,OFF不檢查;無RANGECHECK時要檢查。*RANGCHCK*OFF表示非常重要的數據不作范圍檢查,但重要的數據仍要檢查。如要跳過壓力變化限定的檢查;*RANGECHECK*OFF*NORM*PRESS500*RANGHECK*ON每個數據的變化范圍在CGM的內部自動設置。插入文件(可選)格式:*INCLUDE‘FILENAME’*LIST數據文件的列表(可選)LIST/NOLESTLIST與NOLIST可在一個輸入數據文件內多次交替使用,缺省要列表。*COMMENT改變注釋的定義符(可選)格式:*SOMMENT‘ab’缺省時‘ab’注釋可以單占一行,或放在數據行的后面。*TRANSLATE定義等價的關鍵字(可選)格式:*TRANSLATE‘YOUR—KEYWORD’‘CMG—KEYWORD’目的是定義用戶便于記憶新關鍵字,如:*TRANSLATE‘KX’‘PERMI’*TRANSLATE‘NX’‘NDIVIX’*GRID—ARRAY輸入網格性質數組格式:*GRID—ARRAY(ARRY—QUALIFIER)READ—OTIONDATE(ARRAY—MODIFIER)ARRAY—QUALIFIER稱為數組的限定符,包括:*MAERIX*FRACTURE*RGBLOCK—ADDRESS—RANGE*WELLBOREBLOCK—ADDREESS—RANGE*ANNULUSBLOCK—ADDDRESS—RANGE*TUBINGBLOCK—ADDRESS—RANGE如果ARRAY—QUALIFIER不出現,那么就是*ALLELEM。網格地址范圍(BLOCK—ADDRESS—RANGE)指粗網格塊的,形為:il(:i2)ji(:j2)kl(:k2)讀入選擇為*CON/*IVAR/*JVAR/*ALL/*IJK/*EQUALSIARRAY—MODIFIFR稱為數組的修正,指*MOD。離散井筒處理成基本塊的局部細分網格,基本塊IJK,井筒塊地址為ijk/11,如果井筒循環流動,那么油管地址是Ijk/111,環空地址是Ijk/21。如果給塊IJK賦值,缺省的這個值也賦給它的全部細分塊。*GRID—ARRAY輸入細分網格的物性格式:*GRID—ARRAY*RGil(:i2)jl(:j2)kl(:k2)如像:*POR*RG111*ALL(細分為ni=3,nj=2nk=1)08、079、078、081、080760這種方式可為每個細分塊輸入不同的值,是與基本塊的值不同。應與油藏描述節的*REINE聯合使用(細分的定義)il(:i2)jl(:j2)kl(:k2)是被細分的基本塊的區域。任何的讀入選擇對*RG有效。數據個數應與細分塊的總數相同。*WELLBORE輸入井筒的網格數據格式:*WELLBORE(:i2)jl(:j2)kl(:k2)*ANNULUS(:i2)jl(:j2)kl(:k2)*TUBIBNG(:i2)jl(:j2)kl(:k2)*WELLBORE表示將數據分給網塊的井筒單元。對循環井就分給環空與油管。*ANNULUS(:i2)jl(:j2)kl(:k2)的環空單元。如果無油管也分給井筒。*TUBING(:i2)jl(:j2)kl(:k2)的漬管單元。僅當*WELLBORE/*CIRCWELL出現,以上三種數組限定符才有效。例如:*SO*MATRIX*CON0。7SOANNULUS1:625*CON0.05SOTUBING1:625*CON0例如:*SO*MATRIX*0.70SOWELLBORE1:625*CON0表示環空和油管內SO=0.0。有兩種用法用法1:僅有*WELLBORE,無*ANNULUS與*TUBING。用法2:有*ANNULUS和*TUBING,無*WELLBORE。*CON常數值的數組格式:*CONValue例如:*POR*CON0.16*PERMI*CON100*IJK數組的分區域輸入法格式:*JIK{il(:i2)jl(:j2)kl(:k2)value}{}表示一個網格區域給一個常數,可出現多個{}(表),*IJK讀入法在循環數據節內為有用,或復寫缺省的數據,或復寫前給定的數據,如*KRTYPE*IJK1:51:51:14例如,網格為10×10×3,有三種輸入法:方法1:*POR*IJK1:101:101:30.2461:51:510.17正確方法2:*POR*CON0.246*MOD1:51:51=0.17正確方法3*POR*CON0.246*POR*IJK1:51:510.17錯誤*IVARI方向變化的數組輸入法格式:*IVARvalue(1)。。。value(ni)例如,nx=3DIIVAR3000.04000.05000.0*JVARJ方向變化的數組輸入法格式:JVARvalue(1)。。。value(nj)*KVARK方向變化的數組輸入法]格式:KVARvalue(1)。。。value(nk)例如:nz=3*PERMIKVAR2000.050.0500.0*PERMJ*KVAR2000.050.0500.0PERMKKVAR20.040.060.0*ALL全部網格塊的數組輸入法格式:*ALLvalue(1)。。。value(ni*nj*nk)數的輸入順序按先x后y最后z。例如:ni=10nj=3nk=2*POR*ALL.08.080.81.09.12.1509.097.087.011.15.13408.087.157.145.12.12.135.18.029.047.12.12.154.167.187.121.122.08.08.095.13.12.157.17.18.184.122.084.09.11.12.134.157.157.157.18.18.098.09.0908.09.144.143.123.16.165.102.10.10*Y\EPUAKSI數組的全等輸入法(方賂性)格式:*EQUALSI([*常用于傳導率、滲透率、彌散系數。例如:*PERMIMATRIX*CON100.0*PERMJMATRIX*EQUALSI*RERMKMATRIX*EQUALST*2*PERMI*FRACTURE*CON10000*PERMJ*FEACTURE*EQUALSI*PERMK*FRACTURE*EQUALSI*2*MOD數組的修正(條件)格式:*MODil:I2jl:j2kl:k2(+)value(-)(*)(/)(=)或*MOD(*)value(-)(+)(/)(=)緊接修正的數組后面。如果網格區域不出現,則是整個網格。例如,網格*POR*ALL.08.08.081.09.12.15.09.097.087.011.15.134.08.087.157.145.12.135.18.092..047.12.12.154.167.187.187.121.122.08.08.095.13.12.157.17.18.184.122.084.09.11.12.134.157.157.18.18.098.09.09.08.09.144.143.123.16.165.102.10.10*MOD1:31:41+.01521=.13或*MOD*0.95*INT內插表格數據(可選0如果元素所在位置為*INT,那么用相鄰元素的線性插值生成。例如:*SWT**SwKrwKrow0.20.01.00.30.05*INT0.4*INT0.70.50.40*INT0.6*INT0.50.70.8*INT0.81.00.01.01.00.0*INT只能出現在函數列,崦不能出現在自變量*SW列內。如果整個列僅有一個非0的數,那么其它元素也具有相同的值。輸入輸出控制摘要確定輸入、輸出行動的參數,如文件名、單位、標題、選擇、寫的頻率等。選擇表各種輸入和輸出的文件名只需要輸入數據文件名重啟運算僅需要一個輸入的重啟文件名。Irg全部相容文件名缺省。各文件名的給定或缺省是獨立的。缺省依賴于輸入和輸出文件名的根據名。容易作重啟運算的引導。輸入、輸出單位的選取有SI、油田和實驗三種單位。個別量的單位可以改變。輸入、輸出的單位可以不同。基于mol的量,可用質量代替mol。輸出文件有如下的寫選擇:可以用井、網格和數值解的數據。每個量的頻率、;輸出個數可變。對整個網格提供了長長的變量表。可以控制數組輸出的方向排列。某些量可使用特殊的單位,如ppm,pH等。SR2文件具有如下的寫選擇:可以用井、網格、歷史數據。每個量寫出的頻率及數量(個數)可變。對整個網格供了變量表。某些量可使用特殊的單位,如ppm,PH等。可以用二進制(無格式)或XDR格式寫SR2文件。重啟有如下的選擇:寫的頻率。重繞頻率。以上一個時步或特定時步,讀重啟。中斷有如下的選擇:在清除和關閉文件后立即結束運算。在完成當前時步和寫重啟記錄后結束運行。用指令提示用戶,要求回答(交互式)。必要的數據本節無必要的關鍵字,均為可選的。每個關鍵字均有缺省值可供使用。關鍵字順序如果有*FILENAES,應該是第一個關鍵字,如果有*MASSBASIS,*PARTCLSIEE,*PARTCLWT,必須在*OUTPUT和*OUTSRF之前出現。其它節的使用*MAXTEPS可以出現本節和數值控制節內。本節的關鍵字可以或不可以出現在循環數據節,有下表:可出殃于數值節不能出現在數值節*MAXERROR*TIELE1*RESTART*SRFASCII*TITLE2*MASSBASIS*WRST*TITLE3*PARTCLSIXE*REWIND*CHECKONLY*MAXSTEPS*OUTSOLVR*INUNIT*OUTSRF*WEL*OUTPRN*OUTUNIT*OUTSRF*SPECIAL*WPRN*PRNTORIEN*XDR*OTSRFGRIS*DIM*WSRF*SR2PREC模擬結果文件系統SR·文件系統SR2包含三個文件:指標結果文件。Irf,主要結果文件。Mrf,重繞結果文件。Rrf。后處理的繪圖軟件RESULTS和編寫報告的軟件要使用。Irf和mrf文件。作重啟運算也需要三個文件,如果用*rewind寫重啟記錄,那么對重啟運算rrf也是必須的。只有上次重繞烈軍屬入rrf的信息才可以得出(沖掉了老的重啟信息)。SR2文件系統與以前的SRF文件系統不兼容,為了從SR2文件系統產生一個SRF文件系統,應該使用工具轉換程序SR2RF。內部生成缺省的文件名(文件名包括路徑名)有三種類型的文件:輸入數據文件。Dat,輸出文件。Out,輸入重啟文件。Irf。通過提示或命令行輸入數據文件名,輸出文件名為缺省。通過提示,*FIENAME或命令行輸入‘輸入重啟文件名’。使用缺省功能,僅需改變每個運算的*INDEX—IN文件名,就可以完成一系列自展式的重啟運算。例1、輸入數據文件‘cycle。Dat’,*output的缺省是‘cycle。Out’*INDEX—OUT的缺省是‘cycle。Irf’,*‘runl。Mrf’。例2每一個運算部分用命令“stars。Ese—frunldat”完成,缺省的輸出文件是‘runl。Out’,‘runl。Irf’。如果沒有命令行的選擇—frunl。Dat,程序會提無法用戶回答。對第二個運算部分,將‘runl,dat’拷入‘runl2,dat’。在‘runl2,dat’。在‘runl2,dat’中加入*RESTART。執行命令“stars,ese—frun2,dat-frun2,dat—rrunl”。輸出文件是run2,out,run2,irf和run2,mrf。事實上,runl,irf為run2,ir的你親,當你在用RESULTS和報告編寫軟件(REPORTWRITER)查看run2時,井與特殊歷史的第一部分數據就自動在SR2系統的r—unl內找到。對較小的數據文件,適合于如上所述的每個重啟運算的單獨COPY,但是對較大的數據文件,這就變得不實際了,可以注釋掉無用的數據行,還可以將不改動的數據部分*INCLDE的形成插入這樣文件就可變得較小了。運行模型STARS可以交互式或批處理式的運行。如果你采取了命令行無自變量的交互方式,那么你就鍵入提示要輸入的文件名。每一次提示是對輸入數據文件名,每二次提示是對輸入重啟文憑名(如果是重啟運行的話)。僅當你同時用了*FILENAME和*PROMPT,才提示用戶輸入其它文件名。如果使用了命令一行自變量,就不會產生提示。執行批處理方式有兩種方法,特別適合于多個作業的后臺遠行:方法1、把全部的鍵盤提示改為作業運行原本,原本包含了正確的響應。通過命令行一自變量和/或*FILENAME給定全部需要的文件名,這樣就不需要提示了。方法2、放*SCREEN在輸入數據文件的頭部,這樣數據就可以直接讀入標準的輸入裝置,如象“stars,ese,〈run,dat”。這種情況下,必須給定*OUTPT的文件名,因為程序不知道輸入數據文件名,且無缺省值。Stars。ese命令行一自變量(可選)格式:Stars。Ese(-finput-data)(-rinput-restart)(-checkonly)(-dimsum)(-onestep)Stars。Ese是STARS的請求命令,一般是一個可執行文件的名稱。在UNIX內可以是一個局部的文件,一個同文件的聯系,或僅通過你的尋找路線的一個進入。Input_data數據的文件名。Input_restart用前一次運行生成的輸入重啟文件名irf,缺省:如果Input_data不出現,那么將提示用戶輸入;如果這是一個重啟運行,輸入重啟文件名INPUT—RETART不出現,將提示輸入。*FILENAME輸入/輸出文件名(可選)僅用于覆蓋缺省的文件名,或當其未使用命令行一自變量時給定輸入的重啟文件名。輸入數據文件名必須通過標準的輸入裝置給定(鍵盤/作業運行Script)或命令行一自變量“_f”。格式:*FILENAME(S){file_type(name_option)}文件類型有:*OUTPUT*INDEX—OUT*MAIN—RESULTS—OUT*REWIND—OUT*INDEX—IN*MAIN—RESULTS—IN*MAIN-RESULTS—IN*REWIND—IN*GEOMECHOUTname_option是:‘filename’*PROMPT如果文件類型是*OUTPUT,那么name

_option還可允許*SCREEN*OUTPUT是有格式的模擬結果輸出文件(out)*INDEX—OUT是無格式的指標結果文件(irf)當其開始重啟運行時,將將SR2輸入文件的其些信息拷入SR2的輸出文件。井的歷史數據不拷貝。不抹去輸入的SR2文件,直到它包括的信息不再有用為止。*MAIN—RESULTS—OUT是無格式(二進制)的主要結果文件(mrf)*REWIND—OUT以重繞的方式寫入的重啟數據,稱為可重繞的結果文件(irf)。*INDEX—IN是指標結果文件,從它讀入模擬結果和重啟記錄,此文件僅對重啟運行是必要的。*MAIN—RESULTS—IN是主要結果文件,從它讀入模擬結果和重啟錄(二進制),僅對重啟運算它才是必要的。**REWIND—IN是可重繞的結果文憑,從它讀入重繞、重啟記錄(二進制),僅對重啟運算它才是必要的。GEOMECHOUT地質力學模型生成的有格式輸出文件。如果GEOMECHOUT不存在,或者是PROMPT使用了,那么就把這個輸出寫到主要輸出文件(用*OUTPUT)。‘’空串,將使用一個內部生成的缺省文件名。‘filenam’文件名,可接受的文件名依賴于操作系統。*PROMPT提示用戶通過鍵盤打入文件名。除去INDEX—IN外,全部文件類型都有一個內部生成的缺省文件名,為了使用它,在提示下打入一個空白的反應(大概是立即回車)。*PROMPT提示用戶通過鍵盤打入文件名。除去INDEX—IN外,全部文件類型都有一個內部生成的缺省文件名,為了使用它,在提示下打入一個空白的反應(大概立即回車)。*SCREEN對這種文件類型,數據被送到標準的輸出裝置(屏幕/作業日記文件)。缺省:輸入數據文件無缺省的文件名。用戶必須提供,或是通過標準的輸入裝置(鍵盤/作業運行Script),或命令行一自變量“_f”提供。如果任何其它的必要的文件名來通過*FILENAME(包括*PROMPT)給定,那么就命名用內部生成的缺省文件名。條件:如果*FILENAMES存在,那么它必須是數據文件的第一個關鍵字,如不是這樣,就提示用戶輸入所輸入所需的文件名,而任何以后的*FILENAME將不識別。打開全部帶有‘UNKNOWN’狀態的輸出文件,因而不能避免重復的寫入(OVERWITING)。依賴使用的操作系統,一個老的文件將重復寫入或在后面附加。打開全部帶有‘OLD’狀態的輸入文件,因此它們必須存在。*CHECKONLY檢查錯誤的掃描方式(可選)格式:*CHECKONLY如果*CHECKONLY存在,作業只檢查數據文件的錯誤,包括語法的和超范圍的錯誤。如不存在,則要作時步的計算。*TITLE1*TITLE2*TITLE3方案和情況的標識(可選)格式:*TITLE1‘字串’max40*TITLE2‘字串’max80*TITLE3‘字串’max80*CASEID‘字串’max8*CSDATE‘字串’max10TITLE1最多40字符,TITLE2和*TITLE3最多80字符,它將出現在打印輸出文件和SR2文件內。*CASEID最多8字符用于SR2文件標識曲線。*CSDATE最多10字符,模擬運算的日期,如缺省操作系統提供當天的日期。它產出現在I/O控制節內。例如:*TITLE1’DUALPOSITY/DUALPERMEABILTITYRUNNO。1’*TITLE2’RUNBYABSTARFF,DEC16,1988,C、DCO*TITLE3’4200gridblocks;var,thickeness’*CASEID’NOGAS’*INUNIT輸入/輸出數據的單位(可選)格式:INUNT(SI|*FIELD|*LAB){*EXCEPTqnty_nounit_no}OUTUNIT(SI|*FIELD|*LB){*EXCEPTqnty_nounit_no}定義:*EXCEPT、例外,它允許些量不命名用這種單位。qnty_no下表內的量編號。Unit_no下表內的單位騙號。缺省:如果*INUNIT不存在,則為*SI。如果OUTUNIT不存在,則與輸入單位相同。本手冊各關鍵字的缺省值而無單位,這個單位是*SI。解釋:SR2文件內的數據使用的是STARS的內部單位,與輸入單位無關。但是,指標結果文件。Irf數據使用了輸出單位。例:INUNITSI*EXCEPT22**用0F,not單位表量名稱*SI*FIELD*LAB1、時間daysdaysminutes2、溫度degCDegFdegC3、壓力kpaPsiKpa4、長度mFtcm5、體積m3Ft3Cm36、滲透率mdMdMd7、質量kgLbKg8、mol質量(質量基礎)Gmole(kg)Lbmole(lb)Gmole(kg)9、粘度cpCpCp10、能量JoulesBtuJoules1、流量M3BblCm312、表面張力Dyne/cmByne/cmDyne/cm改變單位的選取名稱量號0123時間1daysHrMinYr溫度2degKDegCDegFDegR壓力3kpaPsiAtmBar長度4mFtCm體積5M3Ft3BblCm3滲透率6M2DarcyMicro-m2Md質量7kgLbMol質量8GmoleLbmol粘度9Kpa-dayKpa-hrCp能量10jBtu流量11M3Ft3BblCm3表面張力12Kpa-mN/MDyne/cm*MASSBASIS質量基礎指示符(可選)格式:*MASSBASIS解釋:某些化學驅過程,要求用質量分數去代替mol分數,缺省單位為mol。例如,組分驟合物的分子量極大,對應的mol分數小,基于mol分數加權的混合規則可能不再合適。*MASSBAIS的作用:輸入的mol數據解釋成質量數據。組分的性質是單位質量,k值定義為相的質量分數,而不是mol分數比。相的組成為質量分數。本手冊輸入每個組分的數據進,包含了它的單位的一般定義。例如,密度單是molarmass/volume,molarmass解釋為mole/massbasis。對SI單位,密度單位為gmol/m#,WJGJVKET*MASSBASIS,那么它為kg/m#。唯一的例外是分子量的定義,對SI單位是mass/mol/mole,即KG/GMOL。有組分汽化的地方不要使用*MASSBASIS選擇,這是因為k=汽/液,必須用mole分數定義。對蒸汽驅過程特別是這樣。*MAXERROR錯誤信息的最大個數(可選)格式:*MAXERRORnnmnum=1~100,缺省為20。數據檢查錯誤,如果錯誤個數超過這兒的num,那么模擬停止運行。*RESTART開始和停算的時間步驟格式:*RESTART(nstart)*MAXSTEPSnstop定義:nstart重啟運算開始計算的時間步數。Nstop是停止計算的時間步數。缺省:如果*RESTART不存在,不需要讀重啟良悄,那么第一個時間步數為1。如果它存在,但無nstart,就命名用*INDEX—IN內的最后一次重啟錄。如果*MAXSTEPS不存在,那么nstop=9999。解釋:為了調試輸入的數據文件,可以采用*MAXSTEPS1。*WRST寫重啟記錄(可選)格式:*WRST(*TIME/freg)*REWIND(num)目的:在輸出重啟文件時,二者控制著寫記錄和重繞記錄的頻率。定義:fweg為整數,fweg是頻率,表示時間步數。*TIME或為*DATE。臺果freg=0不產生重啟記錄。TIME表示在每個TIME或*DATE時刻,要把網格的數組長寫入重啟文件。Num重繞重啟文件的時步數。每隔num=0不重繞。如果num=1,那么僅保留最后一次重啟記錄,記一次重啟又重繞一次。缺省:如果*WRST不存在,就不產生重啟文件。如果后面無*TIME或freg,就認為有*EINE。如果*REWIND不存在,不后果繞重啟文件。如果*REWIND存在但無num,就認為num=1。條件:以上二關鍵字可以出現I/O節與循環數據節內。*WPRN打印輸出的頻率(可選)格式:*WPRN(*GRID/*ITER)(*TIME/freg)目的:控制寫入打印輸出文件的頻率,打例子文件名用*OUTPUT設置。定義:*GRID輸出井和網格的數據.*ITER輸出井墑要、迭代、守恒、數據。Freg每freg步寫一次,為0不寫。*TIME在*TIME或*TATE寫也結果。條件:*WPRN可出現在I/O節或循環數據節。*OUTPRN打印輸出文件的項目(可選)格式:*OUTPRN*GRID{*ALL|*NONE|*EXCEPT|*REMOVE|(special_untit)grid_var}*OUTPRN*WELL{*ALL|*BRIEF|NONE|well_var}*OUTPRN*ITER{BRIEF|NEWTON|UNCONV|TSS|NONE}*PARTCLSIXEvol*PARTCLMLWTcmm目的:給定寫入輸出的打印文件的項目。定義:*GRID指定風格變量寫出。網格變量表為:PRES油相壓力SW水飽和度(也可以WATSAT)SO油飽和度(也可以OILSAT)SG氣飽和度(也可以GASSAT)TEMP溫度Y:氣相的組分組成X:油相的組分組成W:水相的組分組成Z、所有相的組分組成SOLCONC固體組分的濃度HEATLOS(已不使用)OBHLOSS熱損失率VPOROSVOID孔隙度FPOROD流體孔隙度POREVOL孔隙體積VISO水粘度VISG油粘度VISG氣粘度KRW水相相對滲透率KRO油相相對滲透率KRG氣相相對滲透率COW水油摶管力PCOG油氣毛管力MOLDENW水相mol密度MOLDENO油相mol密度MOLDENG氣相mol密度MASDENW水相質量密度MASDENO油相質量量密度MASDENG氣相質量密度RFW水相阻力因子RFD油相阻力因子RFG氣相阻力因子FRCFLOW相的分相流ADSORP組分吸附KRINTERkr的插值值IFT局部表面張力CAPN局部PC數LOGIFT表面張力的自然對數LOGCAPNCAPN的自然對數FLUIDH流體的自然對數FLUIDH流體的熱焓WATERHEAD到等價水柱項的深度(看*DTOP)AQWATCUM水層的凈水浸量AQWATCUE水層的水浸速率AQWATCUM水層的凈熱浸量AQWATCUE水層的熱浸速率以下的量受Special_unit的控制:VLKVCMP關鍵組分的成組,用于計算k=汽/液。LLKVCMP關鍵組分的成組,用于計算k=液/液。VISCCMP關鍵組分的成組,用于油水非線性混合粘度。ADSPCMP關鍵組分的成組,用于kr的計算。PLPMCMP關鍵組分的成組,用于kr的計算。通常y、x、w、x的單位為mol分數,如用*MASSBASIS則為質量分數。Special_unit物殊組分性質,如VLKVCMP,LLKVCMP,VIDCOMP,ADSPCMP,RLPMCMP在報告中用特殊單位。MOLFRmoL分數MASFR質量分數ppm百分之幾VOLFR體積分數MOLAR體積克分子的濃度PHPH=14+LOG(體積克分子濃度)NUM單位相體積內的質點數用*SOLCONC或ADSORRP標識的固體吸附濃度單位:MOLEmoles/孔隙體積MASS質量。孔隙體積VOL固體體積/孔隙體積NUM質點/孔隙體積每種物性可以使用一種特殊的單位,例如:*OUTPUT*GRID*PPM*ADSPCMP*RLPMCMP造成*ADSPCMP組分為PPM,*RLPMCMP組分為number密度*WELL指定的井變量寫入輸出文件,井變量表:LAYPWF層號和BHPLAYPHASE層相率、累積LAYGAS層氣率、累積LAYGAQ層油/水率、累積LAYGAS層氣率、累積WELLCOMP井組分/相摘要TYPECONP井組分/相摘要TYPECOMP注、采組分/相的摘要*ITER打印以下迭代結果BRIEF收斂統計NEWTON牛頓迭代摘要TSSNEWTON+時步選擇UNCONVTSS+不收斂變量的細節(僅用于調試)*ALL選取此表的全部項目*NONE無項目。*EXCEPT除去以后的例外變量表,其余項目全選。*REMOVE從表中去掉以后的項目VOL固相,吸附或重要的拖帶組分的一個質點體積,(m3/ft3/cm3),用于計算*NUM的數值密度。缺省為1。E-11cm3,對應為半徑=1。33E-4cm的地球。

CMM固相、吸附或捕集物質一個mol的質量(kg/gmole)此量僅對(MASS命名用,無缺省值。缺省:可選的關鍵字,如果它不存在,那么缺省為:*OUTPRN(*GRID/*WELL)*NONE*OUTPN*ITER*BRIEF如果*PARTCLSIXE不存在,vol=1,E-11cm3如果*MASSBASIS不存在,濃度的缺省單位為mol。如果*MASSBASIS存在,濃度的缺省單位為mass。條件:此關鍵字出殃在I/O、循環數據節。僅當*MASS出現,*PARTCLMLWT才是必須的。例如:*OUTPRN*GRID*SW*SO*SG*PRES*TEMP*PARTCLMLWT13**COKEmw=131b/lbmole*OUTPRN*GRIDMASSSOLCONE*WSRFSR2的輸入頻率(可選)格式:*WSRF(*GRID/*WELL)(*TIME/freg)目的:控制*OUTSRF內項目寫入SR2的頻率。定義:*GRID控制井數據網格數據的寫出頻率*WELL控制井數數據的寫出頻率。Freg頻率,時間步數*WSRF*GRID0不寫出*WSRF*WELL0不允許*TIME在每一個TIME或DATE時刻寫出缺省:如果無*WSRF*GRID,就缺省為*WSRF*GRID*TIME如果無*WSRF*WELL,就缺省為*WSRF*WELL1。條件:*WSRF可出現在I/O,循環數據節內。*XDR模擬結果文件的項目(可選)格式:*OUTSRF*WELL{comp-untit|*DOWNHOLE|*BLOCKP|*COMPONENT(*ONE|*ALL|comp_list)|*LAYER(*NONE|*ALL|WELL_LIST)}OUTSRF*GRID{(*ALL*NONE|*EXCEPT|*REMOVE|(special_his)*SR2PREC(*SINGLE|*DOUBLE)*SRFASCII目的:標記要寫入模擬結果文件的項目.定義:*XDR可以用外數據的表示,或以用戶平臺的二進制格式寫入二進制數據文件。使用XDR后,就允許在一個平臺上寫成SR2二進文件,又可在另一個平臺上讀出SR2二進文件。例如,在UNIX工作站的服務器上生成SR2文件,然后在一臺PC機上用RESULTS或PEPORTWRITER讀出。如果SR2是XDR格式,那么*XDR將出現在指標文件,irf的頂部附近。XDR是一個非標準的選擇*WELL全部井的數據寫入SR2文件。Com_unit除去體積之外,在質量和/或mol單位存儲井的動態,comp_unit可以是以下之一或之二:*MASS用質量單位存儲*MOLE用mole單位存儲*DOWNHOLE除去地面條件外,還以井底條件寫入井的動態,使SR2文件變大。*BLOCKP如與*LAYER同用,輸出每個井塊層的P。*COMPONENT以組分的形成寫出井的動態,缺省為*NONE,*ALL表示全部組分,增大SR2文件。*LAYER寫出井組層的動態,缺省為*NONE,*ALL表示全部井,也可輸入一個井名表。*GRID網格數據寫入SR2文件。Srf—varsrf—var表為:由*OUTPRN提供的grid—var表全部(但要去掉FRCFLOW,POREVOL,POREVOL,增加WATFRFL,OILFRFL,GASFRFL),還有以下的變量:KVALXW組分的k=氣/水(y/w)KVALYX組分的k=氣/油(y/x)KVALXW組分的k=油/水(x/w)KVALWX組分的k=水/油(w/x)CMPDENW水相的組分質量密度CMPDENO油相的級玢質量密度CMPVISO水相的組分粘度CMPVISO油水相的組分粘度CCHLOSS熱損失PERMIX方向的kPERMJY方向的KPERMKX方向的KWATMOBI方向的水相流度OILMOBI方向的油相流度GASMOBI方向相流度TOTMOBI方向的總流度WATFRFL水相的分相流OILFRFL油相的分相流GASFRFL氣相的分相流VISOVMP油粘度的非線性混合中,關鍵組分的組成。VISWCMP水粘度的非線性混合中,關鍵組分的組成。IFTCMP不使用KRSETNKr數據的組分STRESIX方向的有效應力STRESJY方向的有效應力STRESKZ方向的有效應用力STRESSHy—Z或r—Z內的剪切應力STRESMYP最大主應力STRESMNP最小主應力STRNEPL有效塑性變形Special—unit特殊單位與*OUTPRN描述的相同。*ALL選取表內的全部項目。*NONE不選擇取表內的項目。*EXCEPT除去后面列出的項目,表內其余項目均選定。*REMOVE從以前的表內去掉后繼項目。*Special—his。。。。缺省:如無*OUTSRF*GRID等價于*OUTSRF*GRID*NONE如無*OUTSRF*WELL,結果為:僅寫出體積僅寫出地面條件僅寫出相的僅寫出井的如有*OUTSRF*WELL,結果為:如無*MOLE和*MASS,僅寫體積。如無*DOWNHOLE,僅寫地面條件。如無*COMPONENT,不寫組會的。如無*;AYER,不寫層的。如無*OUTSRF*SPECIAL,不寫特殊歷史。條件:*OUTSRF*GRID可出現在I/O和循環數據節內,對*GRID的缺省頻率是*TIME。*WELL和*SPECIAL只能出現在I/O節。解釋:例如*OUTSRF*GRID*OILSAT*GASSAT*WATSAT*PRES*PARTVCMLWT13**COKEMWIS13IB/IBMOLE*OUTSRF*GRID*MASS*SLXONC*OUTSRF*SPECIAL*BLKVAR*MASS*SOLCONC620*OUTSRF*SPECIAL*BLKVAR*MKLE*SKLCONC620*PRNTORIEN網格數組打印的方向(可選)格式:*PRNTORIENirotatijkord*PRINT_REF(*ON|*OFF)目的:提供打印的方向,提供細分網格的打印。定義:irotat打印的軸旋轉標志旋轉標志行號列號平面0(最緊湊的打印)1IJK2IKJ3KIJ4KJI5JIK6JKIIjkord軸的反向標志反向標志行號列號平面0底層處于頁的底部1正常正常正常2正常正常反向3正常反向正常4正常反向反向5反向正常正常6反向正常反向7反向反向正常8反向反向反向*PRINT—REF用戶可以或不可以用*OUTPRN*GRID打印細網格和基本網格的值。*ON可以打印,*OFF不能打印。缺省:如無*PRNTORIRN則為*PRNTORIEN00,給出最緊溱的打印方式,油藏底層在頁底.如無*PRINT—REF,則為*PRINT—REF*ON。如果有*PRINT—REF,但后面頁*ON或*OFF,則為*PRINT—ON。解釋:缺省的打印方向:在打印紙上從左向右為I,從上往下為j,但后面頁*ON或*OFF,則油層,從上一個xoy平面到下一個xoy平面上為k,油層的底層處于頁底層處于頁底。從左至右I增大,從上往下j增大,從上平面往下平面k增大,某些情況下可以用irotat改變打印的方向。通常,在一行內從左至右指數增大,一列內以上往下指增大,平面與列的增大方向同,可以用ijkrd改變指數增大的方向。提示:*PRNTORIEN造成網格數組以標的順序寫向輸出文件,與網格數組輸入選擇*ALL對應。為使一個運算的輸出文件作為另一個運算的輸入文件,應使用*PRNTORIEN,11,拷貝和粘巾要求的數據到新文件去,而且自動抹去注釋“k=、“j=”、“I=”。*OUTSOLVER矩陣解法的打印出(可選)目的:控制解法軟件包AIMSOL結果的打印輸出。格式:*OUTSOLVER(*ON/OFF)缺省:*OUTSOLVER*OFF解釋:許多有關的量(職塊間的連接)的個數與維值可以打印。這些個數用于建立公共存儲,公共存儲優化了可用的存儲空間。還找印出了內迭代殘量減小的細節。*INTERRUPT收集UNIX信號2中斷(可選)目的:當其檢測出中斷信號后,給定應采取的行動。格式:*INTERRUPT(*INTRRACTIVE/*STOP/*RESTART—STOP)定義:*INTERACTIVE提示用戶以交互方式打入指令,用戶可選取*STOP或*RESTART—STOP。*STOP立即結束運算,當前時步未算完,但要關閉輸出文件,為的是防止文件的破壞。*RESTART—STOP算完當前時步,寫出用*OUTPRN和*OUTSRF給定的輸出,定一個重啟記錄,停止運算。缺省:如無*INTERRRUP,或有但二級子關鍵字無,則為*INTERRUPT*INTERACTIVE。解釋:可以用兩種方式傳送一個UNIX中婁(信號2):手動同時按下Ctrl—c,可中斷當前的交互作業用帶有ID“pid”的命令“kill—2pid”可中斷后臺作業。中斷用于失敗一個作時業,文件可以正常的關閉,基些平臺不清洗輸出文件緩沖區,無中數處理,因些會丟失輸出文件。*DATAECHO詳細的數據回應控制(可選)目的:控制輸入數據文件詳細回應格式:*DATAECHO(*ON/*OFF)缺省:如無*DATAECHO,或有但無二級子關鍵字,則為*DATAECHO*ON解釋:回應的數據構成主要輸出打印文件的一部分分,包括網格定義,組分性質,巖石—流體性質,初始條件、數值控制等。回應打印增大了文件尺寸,待文件正確后可以去掉回應。油藏描述摘要本節數據的分類:模擬網格與網格的細分天然裂縫油藏井筒離散巖石性質蓋層熱損失井筒熱損失熱力水層關鍵字順序*GRID是第一個關鍵字,上述(1)——(3)的關鍵字在*POR之前出現。網格的選擇STARS支持如下兩類網格:有限差分網格(FD)笛卡爾的、徑向、變深/厚度。可控體積有限元網格(CVFE)對FD網格的選擇,*GRID、*DI、*DJ、*DK是必須的,*NINEPOINT、*DIP、*REFINE、*VAMOD、*NULL是可選的。對變深網格,*DTOP可以與初始化的*DWOC、*DGOC和*REFDEPTH同時使用。對CVFE網格選擇,*GRID、*ELEMTONODE、*COORDSI、*COORDSJ、*COORDSK是必須的,*DIP、*CVRFHYBRID是可選的。在使用CVFE選擇時,以下的選擇目前不能使用:。裂縫油藏。離散井筒。熱力水層建議用CMG的交互式的網格生成器建立CVFE網格數據。裂縫油藏選擇裂縫油藏選擇為:*DUALPOR/*DUALPERM/*MINC/*SUBDOMAIC/*VERTNOSEG每個方向的裂縫間距為*DIFRAC、*DJFRAC、*DKFRAC當前不能使用CVFE網格選擇。離散井筒的選擇將井筒離散成若干網格塊,這些網格塊的流動方程與油藏流動方程聯之。井筒離散用*WELLBORE實現,可后加子關鍵字*CIRCWELL表示環空與油管和循環流動,當然環空與油管也對直井部隊發,井微離散可從油層通過上蓋直到地面,非油層部分僅附加很少的風格塊。從井筒到蓋層熱失可用*AQUIFER*WELLBORE模擬。當前井筒離散還不能與CVFE網格同時使。混合網格的離散井筒為了有效地模擬單井的SAGD(SteamAssititedGravityDrsinage),井筒必須直接與它上、下的網格塊相連,不的是允許蒸汽上浮和液體的下沉。將離散的井筒加在混合的細分網格內即可實現。巖石性質孔隙度、滲透率和傳導率的修正用*POR、*PERMI、*PERMJ、*PERMW、TRANSI、*TRACSJ、*TRANSK、*TRANSIJ+、*TRANSIJ-、*TRANSIK+、*TRANSIK-。要使用限空間*MATRIX和*FRATURE區分巖塊和裂縫。當其細分網格用了*REFINE、*MINC、*SUBDOMAIN、*VERTNOSEG、*WELLBORE選擇時,要加限空間*RG。巖石類型*ROCKTYRE、*THTYPE是巖石類型數組。巖石類型可按物性分組:壓縮笥:*CPOR、*CTPOR、*CPORPD、*PORMAX熱物性:*ROCKCP、*THCONP、*THCONW、*THCONO、*THCONG、*CPD蓋層熱損:*HLOSSPROP、*HLOSST、*HLOSSTDIFF蓋層熱選擇對半解析的無限厚蓋層熱損模型來說,熱損的方向與蓋層的熱物性用*HLOSSPROP、*HLOSST、*HLOSSTDIFF給定。對CVFR網格,要用*BOUNDELEM、*BOUNDNODE。井筒熱損失選擇指注蒸汽時井筒的熱損失,須要最少的關鍵字為*RJI、*PTO、*RCI、*EH、*DEPYH,且滿足rti〈rto〈rin熱力水層選擇水層或熱力水層選擇為*AQUIFER。水層的位置、方向、性質給定如下:AAQUIFER(*BOTTOM/*BOUNDARY/*REGION/*WELLBORE)*AQGEOM(*RADIAL/*RECTANG)*AQH*AQPOR*AQVISC*AQRCAP*AQRCND*AQCOMP*AQPERM如是裂縫油藏,僅裂縫與水層連結。當前水層選擇不能與CVFE選擇時使用。井至地表的選擇在模擬從完井段到地表的井筒流動時,具有如下的油藏與井的組合選擇:用*GRID和*DK定義地表到定義地表完井段之間的蓋層網格快。臂如1600m厚的蓋層,用16個100塊模擬,只有兩列井筒塊參加時步的數值計算。對于*KDIR*UP、蓋層塊的編號在油層之后:對于*KDIR*DOWN在油層之前。除去包含井筒的兩縱列網格塊之外,蓋層的其它塊通通作為死塊用*NULL標識。用*POR為其它蓋層(非1、2列)列指定0孔隙度除去油藏的完井段外,在網格的定義部分用*WELLBORE定義穿過蓋層塊的離散井筒。典型地,它是一個斜井筒:油怪部分為水平,且垂直到地表。命名用*AQUIFFER選擇,為蓋層網格塊列內的每塊外面的區域指定解析熱損。每個層的導體處理成徑向塊,這些徑向塊的中心即是井筒,按等體積算出各塊的半徑。熱量從井筒向四周各塊傳遞,使用僅熱損失的水層模型,熱量成徑向地流向上蓋層。在地表井的源匯項與離散井筒聯系在一起(源匯項只在井口處)*GRID網格類型(必須)目的:*GRID標志油藏條件數據輸出的開始。格式:*GRID*CARTninjnk((*KDIR)*UP|*DOWN)GRID*RADIALninjnk(*RWrw)((*KDIR)*UP|*DOWN)*GRID*CVFEneletnnodtncvlyr(*ONIJ/*ONIK/*ONJK)(*KDIR)*UP/*DOWN)定義:*CART長方形的笛卡網格*RADIAL徑向圓柱網格*CVFE可控積有限元網格ninjnk三個方向的網格塊數neletCVFEW全部層的三角形(單)元的總數mmodtCVFE全部層可控體積節點的總數ncvlyrCVFE網格的層數*RWRW最里面一個網格塊摧半徑(m),缺塊省值不能為0。為了精確地復制

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論