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

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