/// <summary> 创建初始重启模型 </summary>
public MainFileRestartSimON InitRestartInfoModel(SimONData mainData)
{
SOLVECTRL start = mainData.Key.Find <SOLVECTRL>();
if (start == null)
{
return(null);
}
MainFileRestartSimON restart = new MainFileRestartSimON();
restart.Parent = null;
restart.RestartTime = start.Date;
restart.FileName = Path.GetFileNameWithoutExtension(mainData.FileName);
restart.FilePath = Path.GetDirectoryName(mainData.FilePath);
restart.Index = 0;
restart.ResultFilePath = restart.FilePath + "\\" + Path.GetFileNameWithoutExtension(restart.FileName) + "_rst1.bin";
restart.MainFilePath = restart.FilePath + "\\" + restart.FileName + KeyConfiger.SimONExtend;
// 把主文件的SCH和SOLU部分传进来
restart.Solution = mainData.Key.Find <SOLUTION>();
restart.Schedule = mainData.Key.Find <SCHEDULE>();
restart.Well = mainData.Key.Find <WELL>();
restart.SchPath = restart.FilePath + "\\" + restart.FileName + "_SCH.DAT";
restart.InitPath = restart.FilePath + "\\" + restart.FileName + "_iNIT.DAT";
return(restart);
}
/// <summary> 将Eclipse数模文件转换成SimON数模文件 </summary>
public SimONData ConvertToSimON(EclipseData ecl)
{
// Todo ：Eclipse里面的修改参数没有解析成SimON中修改参数
ecl.RunModify();
RUNSPEC runspec = ecl.Key.Find <RUNSPEC>();
GRID grid = ecl.Key.Find <GRID>();
SOLUTION solution = ecl.Key.Find <SOLUTION>();
SUMMARY summary = ecl.Key.Find <SUMMARY>();
SCHEDULE schedule = ecl.Key.Find <SCHEDULE>();
REGIONS regions = ecl.Key.Find <REGIONS>();
PROPS props = ecl.Key.Find <PROPS>();
SimONData simon = new SimONData();
simon.FileName = ecl.FileName;
simon.FilePath = ecl.FilePath;
simon.MmfDirPath = ecl.MmfDirPath;
simon.InitConstruct();
simon.X = ecl.X;
simon.Y = ecl.Y;
simon.Z = ecl.Z;
// 模型定义
#region - 起始时间 -
SOLVECTRL tuning = new SOLVECTRL("TUNING");
tuning.Date = ecl.Key.Find <START>().StartTime;
simon.Key.Add(tuning);
#endregion
#region - 维数定义 -
RSVSIZE rsvsize = new RSVSIZE("RSVSIZE");
DIMENS dimens = ecl.Key.Find <DIMENS>();
rsvsize.X = dimens.X;
rsvsize.Y = dimens.Y;
rsvsize.Z = dimens.Z;
simon.Key.Add(rsvsize);
#endregion
#region - 单位类型 -
UnitType unitType = UnitType.METRIC;
// 读到METRIC公制单位
METRIC metric = ecl.Key.Find <METRIC>();
if (metric != null)
{
simon.Key.Add(metric);
unitType = UnitType.METRIC;
}
// 单位类型
FIELD field = ecl.Key.Find <FIELD>();
if (field != null)
{
simon.Key.Add(field);
unitType = UnitType.FIELD;
}
#endregion
#region - 流体类型 -
MODELTYPE modeltype = new MODELTYPE("MODELTYPE");
// 流体类型
OIL oil = runspec.Find <OIL>();
WATER water = runspec.Find <WATER>();
GAS gas = runspec.Find <GAS>();
DISGAS disgas = runspec.Find <DISGAS>();
VAPOIL vapoil = runspec.Find <VAPOIL>();
// 黑油
if (oil != null && water != null && gas != null && disgas != null && vapoil == null)
{
modeltype.MetricType = MetricType.BLACKOIL;
}
// 油水
else if (oil != null && water != null && gas == null && disgas == null && vapoil == null)
{
modeltype.MetricType = MetricType.OILWATER;
}
// 气水
else if (oil == null && water != null && gas != null && disgas == null && vapoil == null)
{
modeltype.MetricType = MetricType.GASWATER;
}
// 挥发油
else if (oil != null && water != null && gas != null && disgas != null && vapoil != null)
{
modeltype.MetricType = MetricType.HFOIL;
}
else
{
modeltype.MetricType = MetricType.BLACKOIL;
}
simon.Key.Add(modeltype);
#endregion
#region - 分区维数 -
EQUILREG equilreg = new EQUILREG("EQUILREG");
FIPREG fipreg = new FIPREG("FIPREG");
ROCKREG rockreg = new ROCKREG("ROCKREG");
SATREG satreg = new SATREG("SATREG");
PVTREG pvtreg = new PVTREG("PVTREG");
simon.Key.Add(equilreg);
simon.Key.Add(fipreg);
simon.Key.Add(rockreg);
simon.Key.Add(satreg);
simon.Key.Add(pvtreg);
TABDIMS tabdims = runspec.Find <TABDIMS>();
if (tabdims != null)
{
fipreg.X = tabdims.Fipfqzds4.ToString();
rockreg.X = tabdims.Yslxgs12.ToString();
satreg.X = tabdims.Bhdbs0.ToString();
pvtreg.X = tabdims.Pvtbs1.ToString();
//fipreg.X = "1";
//rockreg.X = "1";
//satreg.X = "1";
//pvtreg.X = "1";
}
EQLDIMS eqldims = runspec.Find <EQLDIMS>();
if (eqldims != null)
{
//equilreg.X = "1";
equilreg.X = eqldims.Phfqs0.ToString();
}
OVERBURD overburd = props.Find <OVERBURD>();
if (overburd != null)
{
//rockreg.X = overburd.Regions.Count.ToString();
}
EQUILMAP equilmap = new EQUILMAP("EQUILMAP");
FIPMAP fipmap = new FIPMAP("FIPMAP");
ROCKMAP rockmap = new ROCKMAP("ROCKMAP");
SATMAP satmap = new SATMAP("SATMAP");
PVTMAP pvtmap = new PVTMAP("PVTMAP");
if (regions != null)
{
EQLNUM eqlnum = regions.Find <EQLNUM>();
if (eqlnum != null)
{
equilmap = eqlnum.ToTableKey <EQUILMAP>();
solution.Add(equilmap);
eqlnum.Delete();
eqlnum.Dispose();
}
// Todo ：非平衡初始化压力需要转换
var pressure = solution.Find <PRESSURE>();
if (pressure != null)
{
POIL poil = pressure.ToTableKey <POIL>();
solution.Add(poil);
pressure.Delete();
pressure.Dispose();
}
if (regions != null)
{
FIPNUM fipnum = regions.Find <FIPNUM>();
if (fipnum != null)
{
fipmap = fipnum.ToTableKey <FIPMAP>();
grid.Add(fipmap);
fipnum.Delete();
fipnum.Dispose();
}
ROCKNUM rocknum = regions.Find <ROCKNUM>();
if (rocknum != null)
{
rockmap = rocknum.TransToTableKeyByName("ROCKMAP", true) as ROCKMAP;
grid.Add(rockmap);
rocknum.Delete();
rocknum.Dispose();
}
SATNUM satnum = regions.Find <SATNUM>();
if (satnum != null)
{
satmap = satnum.ToTableKey <SATMAP>();
grid.Add(satmap);
satnum.Delete();
satnum.Dispose();
}
PVTNUM pvtnum = regions.Find <PVTNUM>();
if (pvtnum != null)
{
pvtmap = pvtnum.ToTableKey <PVTMAP>();
grid.Add(pvtmap);
pvtnum.Delete();
pvtnum.Dispose();
}
}
}
#endregion
#region - 地质模型 -
simon.Key.Add(grid);
#endregion
#region - 断层 -
//var eclFaults = grid.FindAll<OPT.Product.SimalorManager.RegisterKeys.Eclipse.FAULTS>();
//foreach (var v in eclFaults)
//{
//grid.AddRange(this.ConvertToSimON(v));
//v.Delete();
//}
#endregion
#region - 水体 -
//AQUFETP AQUFETP=
// Todo ：Fetkovich水体数据转换
var ct = solution.Find <OPT.Product.SimalorManager.RegisterKeys.Eclipse.AQUCT>();
if (ct != null)
{
var newFetp = this.ConvertToSimON(ct);
solution.Add(newFetp);
ct.Delete();
}
// Todo ：Fetkovich水体数据转换
var fetp = solution.Find <OPT.Product.SimalorManager.RegisterKeys.Eclipse.AQUFETP>();
if (fetp != null)
{
var newFetp = this.ConvertToSimON(fetp);
solution.Add(newFetp);
fetp.Delete();
}
// Todo ：水体连接数据转换
var aquancon = solution.Find <OPT.Product.SimalorManager.RegisterKeys.Eclipse.AQUANCON>();
if (aquancon != null)
{
var newFetp = this.ConvertToSimON(aquancon);
solution.Add(newFetp);
aquancon.Delete();
}
#endregion
#region - 流体模型 岩石模型-
GRAVITY gravity = ecl.Key.Find <GRAVITY>();
if (gravity != null)
{
// Todo ：SimON只解析绝对密度
DENSITY density = this.ConvertTo(gravity, unitType);
gravity.ParentKey.Add(density);
gravity.Delete();
}
List <IRegionInterface> regSoltionKeys = solution.FindAll <IRegionInterface>();
regSoltionKeys.ForEach(l => l.TransToSimONRegion());
simon.Key.Add(solution);
//
List <IRegionInterface> regPropsKeys = props.FindAll <IRegionInterface>();
regPropsKeys.ForEach(l => l.TransToSimONRegion());
//// Todo ：SGWFN 需要特殊转换为 SWGF
//SGWFN sgwfn = props.Find<SGWFN>();
//if (sgwfn != null)
//{
// //props.AddRange(sgwfn.ConvertTo());
// simon.Key.AddRange<SWGF>(sgwfn.ConvertTo());
//}
simon.Key.Add(props);
#endregion
#region - 初始化模型 -
List <EQUIL> equil = solution.FindAll <EQUIL>();
foreach (var item in equil)
{
EQUILPAR equilpar = new EQUILPAR("EQUILPAR");
EQUIL.Item it = item.GetSingleRegion().GetSingleItem();
equilpar.Szstzdhs0 = it.cksd0;
equilpar.Szstljs1 = it.ckyl1;
equilpar.Ctstyxhs2 = it.ysjmsd2;
equilpar.Ctstyxzdhs3 = it.ysjmcmgyl3.ToDefalt("0");
//equilpar.Jxstzds4 = it.yqjmsd4;
equilpar.E100wgzds5 = it.yqjmsd4;
equilpar.E300jxstzds6 = it.yqjmcmgyl5;
item.ParentKey.Add(equilpar);
item.Delete();
}
#endregion
#region - 生产模型 -
WELL well = new WELL("WELL");
// Todo ：添加完井数据 （注意要放到生产模型前面）
simon.Key.Add(well);
// 生产模型
simon.Key.Add(this.ConvertToSimON(schedule, well, ecl.Key.Find <START>().StartTime, simon.HistoryData));
#endregion
// Todo ：转换修正关键字
List <ModifyKey> modifys = ecl.Key.FindAll <ModifyKey>();
grid.AddRange(modifys);
return(simon);
}
/// <summary> 用重启信息生成新重启数据（只包含主文件，初始化文件和生产文件） </summary>
public SimONData ChangeRestartModel(string mainFilePath, RestartInfoModelSimON model)
{
// 不读取INCLUDE部分数据 只读取SCHEDLE中INCLUDE
SimONData data = FileFactoryService.Instance.ThreadLoadFunc <SimONData>(() => new SimONData(mainFilePath, null, l => l.GetParentKey() is SCHEDULE));
// HTodo ：传递输出参数 需要测试
OUTSCHED rptsched = data.Key.Find <OUTSCHED>();
RPTSUM rptsum = data.Key.Find <RPTSUM>();
var incs = data.Key.FindAll <INCLUDE>();
// 设置所有INCLUDE都不生成文件
incs.ForEach(l => l.IsCreateFile = false);
// 保存主文件
SOLUTION sl = data.Key.Find <SOLUTION>();
SCHEDULE sc = data.Key.Find <SCHEDULE>();
WELLCTRL well = data.Key.Find <WELLCTRL>();
// 更改起始时间
SOLVECTRL tuning = data.Key.Find <SOLVECTRL>();
// Todo ：主文件没有在solotion中找
if (tuning == null)
{
tuning = sl.Find <SOLVECTRL>();
}
tuning.Date = model.RestartTime;
model.Solution.Add(tuning);
// 替换数据
sl.ExChangeData(model.Solution);
sc.ExChangeData(model.Schedule);
well.ExChangeData(model.Well);
//
// 设置保存部分数据
List <INCLUDE> slIncludes = sl.FindAll <INCLUDE>();
slIncludes.ForEach(l => l.IsCreateFile = true);
List <INCLUDE> scIncludes = sc.FindAll <INCLUDE>();
scIncludes.ForEach(l => l.IsCreateFile = true);
List <INCLUDE> wellIncludes = well.FindAll <INCLUDE>();
wellIncludes.ForEach(l => l.IsCreateFile = true);
// 保存主文件（目前没用）
//model.MainData = data;
if (rptsched != null)
{
data.Key.Add(rptsched);
}
if (rptsum != null)
{
data.Key.Add(rptsum);
}
//// Todo ：插入标识到第二个关键字
//HeBianGu.Product.SimalorManager.RegisterKeys.SimON.STEPRST restart = new HeBianGu.Product.SimalorManager.RegisterKeys.SimON.STEPRST("STEPRST");
//data.Key.InsertKey(1, restart);
return(data);
}
/// <summary> 清理父节点 </summary>
public void InitParentKey(bool isCoal = false)
{
this.Key.Clear();
this.Key.CreateSingle <SIMSET>("SIMSET");
GRID grid = this.Key.CreateSingle <GRID>("GRID");
if (grid != null)
{
// 增加屏显
ECHO_OFF echo_off = new ECHO_OFF("ECHO_OFF");
grid.Add(echo_off);
INCLUDE include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_REG.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
grid.Add(include);
include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_GPRO.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
grid.Add(include);
include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_GGO.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
grid.Add(include);
include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_GOTH.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
grid.Add(include);
// 关闭屏显
ECHO_ON echo_on = new ECHO_ON("ECHO_ON");
grid.Add(echo_on);
}
this.Key.CreateSingle <MODIFY>("MODIFY");
WELL well = this.Key.CreateSingle <WELL>("WELL");
USE_TF use_tf = new USE_TF("USE_TF");
if (well != null)
{
well.Add(use_tf);
INCLUDE include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_WELL.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
well.Add(include);
}
PROPS props = this.Key.CreateSingle <PROPS>("PROPS");
if (props != null)
{
INCLUDE include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_SAT.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
props.Add(include);
include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_PVT.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
props.Add(include);
}
if (isCoal)
{
ADSORB adsorb = this.Key.CreateSingle <ADSORB>("ADSORB");
INCLUDE include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_ADS.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
adsorb.Add(include);
}
SOLUTION solution = this.Key.CreateSingle <SOLUTION>("SOLUTION");
if (solution != null)
{
INCLUDE include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_INI.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
solution.Add(include);
include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_AQU.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
solution.Add(include);
}
TUNESET tuneset = this.Key.CreateSingle <TUNESET>("TUNESET");
SOLVECTRL solvectrl = new SOLVECTRL("SOLVECTRL");
tuneset.Add(solvectrl);
//if (solvectrl != null)
//{
// INCLUDE include = new INCLUDE("INCLUDE");
// include.FileName = this.FileName.GetFileNameWithoutExtension() + "_TUN.DAT";
// include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
// solvectrl.Add(include);
//}
//SUMMARY summary = this.Key.CreateSingle<SUMMARY>("SUMMARY");
//if (summary != null)
//{
// INCLUDE include = new INCLUDE("INCLUDE");
// include.FileName = this.FileName.GetFileNameWithoutExtension() + "_sum.inc";
// include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
// summary.Add(include);
// //if (include.Keys.Count == 0)
// //{
// // ALL all = new ALL("ALL");
// // include.Add(all);
// //}
//}
SCHEDULE schedule = this.Key.CreateSingle <SCHEDULE>("SCHEDULE");
if (schedule != null)
{
INCLUDE include = new INCLUDE("INCLUDE");
include.FileName = this.FileName.GetFileNameWithoutExtension() + "_SCH.DAT";
include.FilePath = this.FilePath.GetDirectoryName() + "\\" + include.FileName;
schedule.Add(include);
}
//END end = this.Key.CreateSingle<END>("END");
// 输出关键字
OUTSCHED rptsched = new OUTSCHED("RPTSCHED");
this.Key.Add(rptsched);
}
