Eclipseoffice练习5：用ECLIPSEOFFICE和SIMOPT进行历史拟合

练习5：用ECLIPSEOFFICE和SIMOPT进行历史拟合在该练习中存在一已有的模拟及其相应的井历史含水率数据，并将被输入到ECLIPSEOffice中提供给基本工程，以检验模型与实际井生产间的匹配。工程将被输出到SimOpt中，此处可以使用ECLIPSE100的gradient项来进行匹配分析，并将建立回归问题以改善历史拟合。拟合模型将从SimOpt中输出，并同ECLIPSEOffice工程一起形成新方案。最终，模型可用于预测。可使用SimOpt中的LinearUncertaintyAnalysis模块进行预测分析。该练习被划分为多个阶段，每个阶段是独立的，并能够独立使用。提示：如果用户刚接触£。〃&石Office，则应依次完成所有推荐的阶段。目录：依据现有的模型建立ECLIPSEOffice工程。在SimOpt中，使用输出ECLIPSEOffice工程在SimOpt中建立历史拟合工程。分析当前历史拟合，并设计回归问题。输出SimOpt模型，以输入ECLIPSEOffice中。在ECLIPSEOffice中的多个预测运算。一、依据现有的模型建立ECLIPSEOffice工程在文件夹office_course\exercises\exercise5中启动ECLIPSEOffice（或在别的文件夹中）。ECLIPSEOffice：File|NewProject调新工程SNARKSIMoECLIPSEOffice：Case|Importo输入数据集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面板上）。比较模拟数据和实际测数据现在，建立ECLIPSEOffice工程和运行模拟，其结果可与实测井含水率数据进行比较。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,模型结果和实测结果间拟合较好，但在生产的早期和晚期，二者产水相差较大。这可能是断层间水的流度太大，或Z方向的流度偏高。在该练习的第三阶段（在SimOpt中建立历史拟合工程），将使用SimOpt检查这些参数并尽量获得一较好的拟合结果。然而，在这样做之前，必须为SimOpt准备好模型。“在SimOpt中，使用输出ECLIPSEOffice工程”部分的练习将检查从ECLIPSEOffice工程至USimOpt工程的传递。二、在SimOpt中，使用输出ECLIPSEOffice工程如果没有进行上面的操作，应使用SNARKSIM.DATA建立ECLIPSEOffice工程。提示：如果不知道怎样建立后。乙/08石Office工程，可参见第一步中的相关内容。为SimOpt做模型准备的第一步是：在ECLIPSEOffice工程中建立新方案。1ECLIPSEOffice：Case|AddCase|Clone所创建的新方案包含了与缺省方案相同的信息。注意：SimOpt将为输入的数据添加有_1”脚标的文件名（如BASE.DATA变为BASE_1.DATA）。因此SimOpt版的模型将不会受到原始数据集的影响。它将拒绝文件名“_1”脚标的数据文件。2在CaseNames树中，单击新方案，以激活它。ECLIPSEOffice：Case|Show改变方案名为：ORIGINAL.单击OK。ECLIPSEOffice：File|SaveProjecto建立INIT文件SimOpt需要INIT文件和数据集。在从ECLIPSEOffice输出模型ORIGINAL之前，应创建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的当前内容。DataMangagerModule：File|CloseOffice：File|RunManager.提示：该运行仅产生/Na文件（对模型没有起任何改变）。模拟可被关闭以保存时间。单击TurnOffSimulationoRunManagerModule：File|WritData运行结束后，文件必须写入（不带NOSIM关键字），为SimOpt读入它做准备。单击TurnOffSimulation以便不再检测到它。RunManagerModule：File|WriteData单击OK,覆盖old.DATA文件。讨论SimOpt需要一工作模型和INIT文件来建立工程。ECLIPSEOffice提供了一套检查数据和编辑关键字的功能来为SimOpt历史拟合工程创建文件。在SimOpt中建立历史拟合工程启动SimOpt,并调入工程SNARKoSimOpt：Project|Import|SimulationModel输入文件ORIGINAL_E100.DATASimOpt：Project|Import|Initfile输入文件ORIGINAL_E100.INITSimOpt：Project|Import|ObservedData|GrafUserDataFormat输入文件GRAF.WWCTo到现在为止，SimOpt工程的基本数据已经输入，应保存工程。SimOpt：ProjectISave注意：输入的井含水率数据Wwct）没有评价其测量误差。因此，对于加。加有必要计算实测数据和模拟结果间的RMS（均方误差），且必须手工加入。SimOpt：DataIControlList单击BlockAction（SelectedSet）上的SetError。35在Error文本框中输入误差值0.01。单击ApplyBlockAction0单击Closeo敏感性模型现在准备用模型结果数据和实测数据来进行敏感性分析。在该方案中，模型中的三个断层和Z方向的传导率将用到敏感性模型中。SimOpt：Parameter|Add选择ParmType列表中的ZTrans,并打击Applyo选择Domain按钮上的Faulto选择ParmType列表中的FaultTr。在DomainNames列表中单击FAULTS1,然后按住Shift,单击FAULTS3,选中所有的三个断层。单击Applyo单击Closeo45单击Version1框（ParameterVersions数据树上）左侧的"+”符号（在主屏幕的中央区域）。定义传导率乘法器断层传导率乘法器沿每个断层的整个长度是有效的。然而，Z方向的传导率乘法器对大多数网格均有效，因此应该定义它的有效区域。46SimOpt：Parameter|RegionEditor47RegionEditor：Region|Define48在表的第1]行，输入如下的数二据，然后单击OK：I1(nI:14)I2(nI:14)J1(nJ:14)J2(nJ:14)K1(nK:14)K2(nK:14)11411415单击OK。RegionEditor：FileIExit。定义开始值现在，应定义我们所关注的参数，设置它们的初始值和值的范围。SimOpt：ParameterIControlList提示：断层的传导率乘法器的变化范围是（封闭断层）到/（无断层存在）。如果乘法器没有赋初值，可采用以下经验准则：设置近封闭断层（敏感性较大）的情况，将对储层流动产生较大的影响（比起无断层存在时）。因此，其初始值通常取J，并允许该值在0到1间变化。提示：估计垂向传导率的误差系数为0。或更大。因此，Z方向的传导率值可以在初始值的基础上增加或降低100倍，即从0.1到100。52在ParameterControl面板的表中输入如下数据:TypeDomainModifierModifierMinModifierMaxZTransRegion110.01100FaultTrFAULTS10.011E-51FaultTrFAULTS20.011E-51FaultTrFAULTS30.011E-51单击ParameterControl面板上的Apply。单击Closeo讨论在SimOpt中的工程包含一模型、INIT文件、实测井生产数据和完善历史拟合所需的参数。现在可以完成工程，并能够计算敏感性。选择ECLIPSE100中的Gradients项，能够在单一的模拟运行中进行敏感性分析。接下来是检查敏感性分析，为回归提供基础以完善历史拟合。三、分析当前历史拟合并设计回归问题SimOpt：Simulate在RunMode按钮上选择Gradients并单击Start。提示：右击ObservedData菜单中的井PROD1节点并选择Plot，可以观察模拟过程中的进程。这将做实测数据的曲线，以与模拟结果相重叠。3当模拟结束后，主窗体底部的提示信息窗口将给出提示信息。所有实测数据和模拟结果数据间的均方根吻合度大约为24。注意：FAULTS2和FAULTS3的乘法器被当作多余信息。这意味它们对于数据拟合的影响是重要，因此可以将它们关掉。SimOpt：Parameter|ControlList在表的Activity列的下拉菜单中设置FAULTS2和FAULTS3为Inactiveo单击Applyo单击Closeo观察其余的参数在分析中，没有涉及到其余的参数，而这些参数是可以观察的。SimOpt：Match|MatchAnalysis单击UpdateoRMS文件夹显示：井PROD2的RMS和敏感性为0。因此，该井的拟合不受到参数选择变化的影响，所以可从系统中提取实测数据。SimOpt：Data|ControlListo在表中，设置PROD2的WWCT数据状态为Inactove,并单击Applyo注意：这样并不会从工程中删除数据。仅可从邓计算中移走数据。通过设置数据状态为ctive，可将数据添加回RMS计算中。单击Closeo启动历史拟合分析的按钮13单击Updateo这将重新进行分析计算（不包括井PROD2的数据）。单击Correlation文件夹。关系矩阵的非对角线元素大约是0.33,因此两个剩余参数间的联系不足够紧密。如果参数间联系过紧，则会出现数学上的ill-posed问题（即回归无法改善拟合）。在问题是well-posed的情况下，是可以使用非线性回归的方法来迭代更新模拟值，以达到改善历史拟合的目的。注意：应记住重要的一点是：分析问题是线性的，但在乘法器中改变模型后就变为了非线性问题。因此，超过3或4次迭代回归的方法是不可取的（因为每一次迭代会使模型进一步远离有效值）。SimOpt：Match|Simulate在RunMode单选框上选择Regresso设置MaxRunsSimulator为4。单击Starto提示：如果想观察RMSfit进程，可单击刊otRms按钮（每次模拟后，这将更新工程）。同时，井数据的时间监视将有效。最终的RMSfit是10.3,对原始数据集的较大改动将输入到SimOpt中。19从实测数据树上选择PROD1WWCT项，观察实测数据和模拟结果数据间的拟合。20选择Data|Plot,观察改善了的历史拟合。讨论灵敏性分析过后，是利用SimOpt中的工具来进行递变率分析。该分析提供了井条件回归问题的功能。SimOpt可用于运行ECLIPSE100模拟器，多次更新所选择的参数，以完善历史拟合。四、输出SimOpt模型，以便输入到ECLIPSEOffice中根据历史拟合所形成的模型可从SimOpt中输出，并用于ECLIPSEOffice中的预测。SimOpt：Project|Export|SimulationModelo在ProjectExportComments面板上单击OK按钮。以MATCHED.DATA名保存数据集。SimOpt：Project|Exit。5单击Yes,保存工程。讨论拟合模拟（即乘法器的应用）可从SimOpt输出。这将在原来的模拟中添加乘法器关键字，以便能在以后的应用中，恢复到原始模拟。五、在ECLIPSEOffice中的预测1ECLIPSEOffice：CaseIAddCaseINew。2在CaseNames数据树上，选择新方案。ECLIPSEOffice：Case|Show。改变方案名为：PREDICTION.Office：Case|Import输入MATCHED.DATAo添加预测时间段现在，历史拟合的模型是在ECLIPSEOffice中应用，接下来是延长模拟时间段（即添加预测时间段）。ECLIPSEOffice：Module|DataManagerDataManager：Section|Schedule9在Time-Dates列表中，选择数据集中最后的日期（1January2030）oScheduleSection：Time|Insert在NewTimeEntryChoice单选框上选择TimeStepo12输入365days作为时间步。13在NewTime面板上的Num输入框中输入5。14单击OK。这将添加五个模拟时间步长（每个均为一年）到当前的数据集中。定义流量控制以下是在预测时间段内定义产量控制。ECLIPSEOffice：Event|New在EventTypes列表中选择WellControls和Limitso在Event列表中选择ProductionWellControlo18单击OK。19在ScheduleSection中为新关键字WCONPROD专w入如下信息：WellPROD1Open/ShutflagOpenControlORATOilRate10