Eclipseoffice練習(xí)5：用ECLIPSEOFFICE和SIMOPT進(jìn)行歷史擬合

練習(xí)5：用ECLIPSEOFFICE和SIMOPT進(jìn)行歷史擬合在該練習(xí)中存在一已有的模擬及其相應(yīng)的井歷史含水率數(shù)據(jù)，并將被輸入到ECLIPSEOffice中提供給基本工程，以檢驗(yàn)?zāi)Ｐ团c實(shí)際井生產(chǎn)間的匹配。工程將被輸出到SimOpt中，此處可以使用ECLIPSE100的gradient項(xiàng)來(lái)進(jìn)行匹配分析，并將建立回歸問(wèn)題以改善歷史擬合。擬合模型將從SimOpt中輸出，并同ECLIPSEOffice工程一起形成新方案。最終，模型可用于預(yù)測(cè)?？墒褂肧imOpt中的LinearUncertaintyAnalysis模塊進(jìn)行預(yù)測(cè)分析。該練習(xí)被劃分為多個(gè)階段，每個(gè)階段是獨(dú)立的，并能夠獨(dú)立使用。提示：如果用戶(hù)剛接觸￡?！?石Office，則應(yīng)依次完成所有推薦的階段。目錄：依據(jù)現(xiàn)有的模型建立ECLIPSEOffice工程。在SimOpt中，使用輸出ECLIPSEOffice工程在SimOpt中建立歷史擬合工程。分析當(dāng)前歷史擬合，并設(shè)計(jì)回歸問(wèn)題。輸出SimOpt模型，以輸入ECLIPSEOffice中。在ECLIPSEOffice中的多個(gè)預(yù)測(cè)運(yùn)算。一、依據(jù)現(xiàn)有的模型建立ECLIPSEOffice工程在文件夾office_course\exercises\exercise5中啟動(dòng)ECLIPSEOffice（或在別的文件夾中）。ECLIPSEOffice：File|NewProject調(diào)新工程SNARKSIMoECLIPSEOffice：Case|Importo輸入數(shù)據(jù)集SNARKSIM.DATAoECLIPSEOffice：Module|RunManageroRunManager：Submit|Runs。ECLIPSEOffice:Module|ResultsVieweroResultViewerModule：File|OpenCurrentCase|SUMMARYo在Extract/LoadSummaryVectors面板上，選擇ReadAllSummaryVectors0在Extract/LoadSummaryVectors面板上，選擇ReadAllReports。在Extract/LoadSummaryVectors面板上，單擊Load。13單擊丫65,輸入所有的變量。ResultsViewerModule：View|SummaryVectorsSummaryVectors：New|User在UserTemplates面板的X-Axis列表中，選擇TIME。在UserTemplates面板的X-Axis列表中，選擇WWCT：PROD1。單擊AddToList按鈕。單擊OK（UserTemplates面板上）。比較模擬數(shù)據(jù)和實(shí)際測(cè)數(shù)據(jù)現(xiàn)在，建立ECLIPSEOffice工程和運(yùn)行模擬，其結(jié)果可與實(shí)測(cè)井含水率數(shù)據(jù)進(jìn)行比較。SummaryVectors：File|Observed|ColumnFormato打開(kāi)文件prod1.wwcto單擊OK,使用ColumnFormatUserData面板中的缺省列格式。SummaryVectors：New|Observed在SelectObservedVectorsforplotting面板的X-AxisVectors列表中選擇TIMEo在SelectObservedVectorsforplotting面板的X-AxisVectors列表中選擇WWCT：PROD1o單擊AddtoListo選擇AddtoGrapho單擊OKo提示：如果所作圖線(xiàn)是新圖線(xiàn)，則可略去第7步。討論對(duì)于井PROD1,模型結(jié)果和實(shí)測(cè)結(jié)果間擬合較好，但在生產(chǎn)的早期和晚期，二者產(chǎn)水相差較大。這可能是斷層間水的流度太大，或Z方向的流度偏高。在該練習(xí)的第三階段（在SimOpt中建立歷史擬合工程），將使用SimOpt檢查這些參數(shù)并盡量獲得一較好的擬合結(jié)果。然而，在這樣做之前，必須為SimOpt準(zhǔn)備好模型?！霸赟imOpt中，使用輸出ECLIPSEOffice工程”部分的練習(xí)將檢查從ECLIPSEOffice工程至USimOpt工程的傳遞。二、在SimOpt中，使用輸出ECLIPSEOffice工程如果沒(méi)有進(jìn)行上面的操作，應(yīng)使用SNARKSIM.DATA建立ECLIPSEOffice工程。提示：如果不知道怎樣建立后。乙/08石Office工程，可參見(jiàn)第一步中的相關(guān)內(nèi)容。為SimOpt做模型準(zhǔn)備的第一步是：在ECLIPSEOffice工程中建立新方案。1ECLIPSEOffice：Case|AddCase|Clone所創(chuàng)建的新方案包含了與缺省方案相同的信息。注意：SimOpt將為輸入的數(shù)據(jù)添加有_1”腳標(biāo)的文件名（如BASE.DATA變?yōu)锽ASE_1.DATA）。因此SimOpt版的模型將不會(huì)受到原始數(shù)據(jù)集的影響。它將拒絕文件名“_1”腳標(biāo)的數(shù)據(jù)文件。2在CaseNames樹(shù)中，單擊新方案，以激活它。ECLIPSEOffice：Case|Show改變方案名為：ORIGINAL.單擊OK。ECLIPSEOffice：File|SaveProjecto建立INIT文件SimOpt需要INIT文件和數(shù)據(jù)集。在從ECLIPSEOffice輸出模型ORIGINAL之前，應(yīng)創(chuàng)建INIT文件。ECLIPSEOffice：Module|DataManagerDataManagerModule：Section|GridGridSection：Subsection|GRIDKeywordsGRIDKeywordSection：View|Description在KeywordType列表中選擇OperationalKeywords。GRIDKeywordSection：EditIInserKeyword在KeywordSelection面板中選擇OutputINIT文件。在GRIDKeywordSection面板上，單擊Apply。GridKeywordSection：File|CloseGridSection：File|Close單擊Yes，保存GridSection的當(dāng)前內(nèi)容。DataMangagerModule：File|CloseOffice：File|RunManager.提示：該運(yùn)行僅產(chǎn)生/Na文件（對(duì)模型沒(méi)有起任何改變）。模擬可被關(guān)閉以保存時(shí)間。單擊TurnOffSimulationoRunManagerModule：File|WritData運(yùn)行結(jié)束后，文件必須寫(xiě)入（不帶NOSIM關(guān)鍵字），為SimOpt讀入它做準(zhǔn)備。單擊TurnOffSimulation以便不再檢測(cè)到它。RunManagerModule：File|WriteData單擊OK,覆蓋old.DATA文件。討論SimOpt需要一工作模型和INIT文件來(lái)建立工程。ECLIPSEOffice提供了一套檢查數(shù)據(jù)和編輯關(guān)鍵字的功能來(lái)為SimOpt歷史擬合工程創(chuàng)建文件。在SimOpt中建立歷史擬合工程啟動(dòng)SimOpt,并調(diào)入工程SNARKoSimOpt：Project|Import|SimulationModel輸入文件ORIGINAL_E100.DATASimOpt：Project|Import|Initfile輸入文件ORIGINAL_E100.INITSimOpt：Project|Import|ObservedData|GrafUserDataFormat輸入文件GRAF.WWCTo到現(xiàn)在為止，SimOpt工程的基本數(shù)據(jù)已經(jīng)輸入，應(yīng)保存工程。SimOpt：ProjectISave注意：輸入的井含水率數(shù)據(jù)Wwct）沒(méi)有評(píng)價(jià)其測(cè)量誤差。因此，對(duì)于加。加有必要計(jì)算實(shí)測(cè)數(shù)據(jù)和模擬結(jié)果間的RMS（均方誤差），且必須手工加入。SimOpt：DataIControlList單擊BlockAction（SelectedSet）上的SetError。35在Error文本框中輸入誤差值0.01。單擊ApplyBlockAction0單擊Closeo敏感性模型現(xiàn)在準(zhǔn)備用模型結(jié)果數(shù)據(jù)和實(shí)測(cè)數(shù)據(jù)來(lái)進(jìn)行敏感性分析。在該方案中，模型中的三個(gè)斷層和Z方向的傳導(dǎo)率將用到敏感性模型中。SimOpt：Parameter|Add選擇ParmType列表中的ZTrans,并打擊Applyo選擇Domain按鈕上的Faulto選擇ParmType列表中的FaultTr。在DomainNames列表中單擊FAULTS1,然后按住Shift,單擊FAULTS3,選中所有的三個(gè)斷層。單擊Applyo單擊Closeo45單擊Version1框（ParameterVersions數(shù)據(jù)樹(shù)上）左側(cè)的"+”符號(hào)（在主屏幕的中央?yún)^(qū)域）。定義傳導(dǎo)率乘法器斷層傳導(dǎo)率乘法器沿每個(gè)斷層的整個(gè)長(zhǎng)度是有效的。然而，Z方向的傳導(dǎo)率乘法器對(duì)大多數(shù)網(wǎng)格均有效，因此應(yīng)該定義它的有效區(qū)域。46SimOpt：Parameter|RegionEditor47RegionEditor：Region|Define48在表的第1]行，輸入如下的數(shù)二據(jù)，然后單擊OK：I1(nI:14)I2(nI:14)J1(nJ:14)J2(nJ:14)K1(nK:14)K2(nK:14)11411415單擊OK。RegionEditor：FileIExit。定義開(kāi)始值現(xiàn)在，應(yīng)定義我們所關(guān)注的參數(shù)，設(shè)置它們的初始值和值的范圍。SimOpt：ParameterIControlList提示：斷層的傳導(dǎo)率乘法器的變化范圍是（封閉斷層）到/（無(wú)斷層存在）。如果乘法器沒(méi)有賦初值，可采用以下經(jīng)驗(yàn)準(zhǔn)則：設(shè)置近封閉斷層（敏感性較大）的情況，將對(duì)儲(chǔ)層流動(dòng)產(chǎn)生較大的影響（比起無(wú)斷層存在時(shí)）。因此，其初始值通常取J，并允許該值在0到1間變化。提示：估計(jì)垂向傳導(dǎo)率的誤差系數(shù)為0?；蚋?。因此，Z方向的傳導(dǎo)率值可以在初始值的基礎(chǔ)上增加或降低100倍，即從0.1到100。52在ParameterControl面板的表中輸入如下數(shù)據(jù):TypeDomainModifierModifierMinModifierMaxZTransRegion110.01100FaultTrFAULTS10.011E-51FaultTrFAULTS20.011E-51FaultTrFAULTS30.011E-51單擊ParameterControl面板上的Apply。單擊Closeo討論在SimOpt中的工程包含一模型、INIT文件、實(shí)測(cè)井生產(chǎn)數(shù)據(jù)和完善歷史擬合所需的參數(shù)?，F(xiàn)在可以完成工程，并能夠計(jì)算敏感性。選擇ECLIPSE100中的Gradients項(xiàng)，能夠在單一的模擬運(yùn)行中進(jìn)行敏感性分析。接下來(lái)是檢查敏感性分析，為回歸提供基礎(chǔ)以完善歷史擬合。三、分析當(dāng)前歷史擬合并設(shè)計(jì)回歸問(wèn)題SimOpt：Simulate在RunMode按鈕上選擇Gradients并單擊Start。提示：右擊ObservedData菜單中的井PROD1節(jié)點(diǎn)并選擇Plot，可以觀察模擬過(guò)程中的進(jìn)程。這將做實(shí)測(cè)數(shù)據(jù)的曲線(xiàn)，以與模擬結(jié)果相重疊。3當(dāng)模擬結(jié)束后，主窗體底部的提示信息窗口將給出提示信息。所有實(shí)測(cè)數(shù)據(jù)和模擬結(jié)果數(shù)據(jù)間的均方根吻合度大約為24。注意：FAULTS2和FAULTS3的乘法器被當(dāng)作多余信息。這意味它們對(duì)于數(shù)據(jù)擬合的影響是重要，因此可以將它們關(guān)掉。SimOpt：Parameter|ControlList在表的Activity列的下拉菜單中設(shè)置FAULTS2和FAULTS3為Inactiveo單擊Applyo單擊Closeo觀察其余的參數(shù)在分析中，沒(méi)有涉及到其余的參數(shù)，而這些參數(shù)是可以觀察的。SimOpt：Match|MatchAnalysis單擊UpdateoRMS文件夾顯示：井PROD2的RMS和敏感性為0。因此，該井的擬合不受到參數(shù)選擇變化的影響，所以可從系統(tǒng)中提取實(shí)測(cè)數(shù)據(jù)。SimOpt：Data|ControlListo在表中，設(shè)置PROD2的WWCT數(shù)據(jù)狀態(tài)為Inactove,并單擊Applyo注意：這樣并不會(huì)從工程中刪除數(shù)據(jù)。僅可從鄧計(jì)算中移走數(shù)據(jù)。通過(guò)設(shè)置數(shù)據(jù)狀態(tài)為ctive，可將數(shù)據(jù)添加回RMS計(jì)算中。單擊Closeo啟動(dòng)歷史擬合分析的按鈕13單擊Updateo這將重新進(jìn)行分析計(jì)算（不包括井PROD2的數(shù)據(jù)）。單擊Correlation文件夾。關(guān)系矩陣的非對(duì)角線(xiàn)元素大約是0.33,因此兩個(gè)剩余參數(shù)間的聯(lián)系不足夠緊密。如果參數(shù)間聯(lián)系過(guò)緊，則會(huì)出現(xiàn)數(shù)學(xué)上的ill-posed問(wèn)題（即回歸無(wú)法改善擬合）。在問(wèn)題是well-posed的情況下，是可以使用非線(xiàn)性回歸的方法來(lái)迭代更新模擬值，以達(dá)到改善歷史擬合的目的。注意：應(yīng)記住重要的一點(diǎn)是：分析問(wèn)題是線(xiàn)性的，但在乘法器中改變模型后就變?yōu)榱朔蔷€(xiàn)性問(wèn)題。因此，超過(guò)3或4次迭代回歸的方法是不可取的（因?yàn)槊恳淮蔚鷷?huì)使模型進(jìn)一步遠(yuǎn)離有效值）。SimOpt：Match|Simulate在RunMode單選框上選擇Regresso設(shè)置MaxRunsSimulator為4。單擊Starto提示：如果想觀察RMSfit進(jìn)程，可單擊刊otRms按鈕（每次模擬后，這將更新工程）。同時(shí)，井?dāng)?shù)據(jù)的時(shí)間監(jiān)視將有效。最終的RMSfit是10.3,對(duì)原始數(shù)據(jù)集的較大改動(dòng)將輸入到SimOpt中。19從實(shí)測(cè)數(shù)據(jù)樹(shù)上選擇PROD1WWCT項(xiàng)，觀察實(shí)測(cè)數(shù)據(jù)和模擬結(jié)果數(shù)據(jù)間的擬合。20選擇Data|Plot,觀察改善了的歷史擬合。討論靈敏性分析過(guò)后，是利用SimOpt中的工具來(lái)進(jìn)行遞變率分析。該分析提供了井條件回歸問(wèn)題的功能。SimOpt可用于運(yùn)行ECLIPSE100模擬器，多次更新所選擇的參數(shù)，以完善歷史擬合。四、輸出SimOpt模型，以便輸入到ECLIPSEOffice中根據(jù)歷史擬合所形成的模型可從SimOpt中輸出，并用于ECLIPSEOffice中的預(yù)測(cè)。SimOpt：Project|Export|SimulationModelo在ProjectExportComments面板上單擊OK按鈕。以MATCHED.DATA名保存數(shù)據(jù)集。SimOpt：Project|Exit。5單擊Yes,保存工程。討論擬合模擬（即乘法器的應(yīng)用）可從SimOpt輸出。這將在原來(lái)的模擬中添加乘法器關(guān)鍵字，以便能在以后的應(yīng)用中，恢復(fù)到原始模擬。五、在ECLIPSEOffice中的預(yù)測(cè)1ECLIPSEOffice：CaseIAddCaseINew。2在CaseNames數(shù)據(jù)樹(shù)上，選擇新方案。ECLIPSEOffice：Case|Show。改變方案名為：PREDICTION.Office：Case|Import輸入MATCHED.DATAo添加預(yù)測(cè)時(shí)間段現(xiàn)在，歷史擬合的模型是在ECLIPSEOffice中應(yīng)用，接下來(lái)是延長(zhǎng)模擬時(shí)間段（即添加預(yù)測(cè)時(shí)間段）。ECLIPSEOffice：Module|DataManagerDataManager：Section|Schedule9在Time-Dates列表中，選擇數(shù)據(jù)集中最后的日期（1January2030）oScheduleSection：Time|Insert在NewTimeEntryChoice單選框上選擇TimeStepo12輸入365days作為時(shí)間步。13在NewTime面板上的Num輸入框中輸入5。14單擊OK。這將添加五個(gè)模擬時(shí)間步長(zhǎng)（每個(gè)均為一年）到當(dāng)前的數(shù)據(jù)集中。定義流量控制以下是在預(yù)測(cè)時(shí)間段內(nèi)定義產(chǎn)量控制。ECLIPSEOffice：Event|New在EventTypes列表中選擇WellControls和Limitso在Event列表中選擇ProductionWellControlo18單擊OK。19在ScheduleSection中為新關(guān)鍵字WCONPROD專(zhuān)w入如下信息：WellPROD1Open/ShutflagOpenControlORATOilRate10