From fe5bb447bbc9d325246cf99b980556af62955254 Mon Sep 17 00:00:00 2001
From: lx <ex_lixiang17@cosl.com.cn>
Date: 星期五, 18 七月 2025 17:32:19 +0800
Subject: [PATCH] plot

---
 ErrorAnalysis.UI/FrmMain.cs |  154 ++++++++++++++++++++++++++++++++++++++++++++-------
 1 files changed, 133 insertions(+), 21 deletions(-)

diff --git a/ErrorAnalysis.UI/FrmMain.cs b/ErrorAnalysis.UI/FrmMain.cs
index 7e35121..02cbec4 100644
--- a/ErrorAnalysis.UI/FrmMain.cs
+++ b/ErrorAnalysis.UI/FrmMain.cs
@@ -1,7 +1,15 @@
 using ErrorAnalysis.Repository;
 using ErrorAnalysis.Repository.Entity;
 using ErrorAnalysis.Service;
+using ErrorAnalysis.Service.Model;
+using ErrorAnalysis.UI.Utility;
+using OxyPlot;
+using OxyPlot.Axes;
+using OxyPlot.Legends;
+using OxyPlot.Series;
+using OxyPlot.WindowsForms;
 using System.ComponentModel.DataAnnotations;
+using System.Security.Cryptography;
 using System.Windows.Forms;
 using static System.Windows.Forms.VisualStyles.VisualStyleElement.Menu;
 
@@ -50,12 +58,13 @@
                 txtDataPath.Text = fileDlg.FileName;
                 txtToolType.Text = fileDlg.SafeFileName.Substring(0, fileDlg.SafeFileName.IndexOf("Mcnp"));
                 _connectionString = $"Data Source={fileDlg.FileName};";
-                var modelRepository = new COModelRepository(_connectionString);
+                RepositoryInstance.Instance.ReLoad(_connectionString);
+
                 pbLoadDB.Value = 50;
                 _filters = new Dictionary<string, string>();
                 Task.Run(() =>
                 {
-                    _models = modelRepository.GetCOModels();
+                    _models = RepositoryInstance.Instance.COModelRepository.GetCOModels();
                     BeginInvoke(new System.Action(() =>
                     {
                         pbLoadDB.Value = 80;
@@ -103,9 +112,20 @@
                 return;
             }
 
+            if (!double.TryParse(nudSoPrecisionRequired.Value.ToString(), out double targetErrorRatio) || nudSoPrecisionRequired.Value <= 0)
+            {
+                MessageBox.Show("Please enter a reasonable so precision required");
+                return;
+            }
+
             if (cmbLithology.SelectedIndex == -1)
             {
                 MessageBox.Show("Please select a lithology");
+                return;
+            }
+            if (curModels.Count == 0)
+            {
+                MessageBox.Show("Parameters error!");
                 return;
             }
 
@@ -113,9 +133,8 @@
                 _model = curModels[0];
             _model = curModels.FirstOrDefault(m => m.GravelFillPercent == nudGravelFIllPercent.Value.ToString());
             if (_model == null)
-            {
+                _model = curModels.First();
 
-            }
             if (rdoLockSpeed.Checked)
             {
                 if (!double.TryParse(nudSpeed.Value.ToString(), out double speed) || nudSpeed.Value <= 0)
@@ -123,6 +142,10 @@
                     MessageBox.Show("Please enter a reasonable speed");
                     return;
                 }
+                var calcSpeed = cmbSpeedUnit.Text == "ft/hr" ? UnitConvert.FtHr2MMin(speed) : UnitConvert.MHR2MMin(speed);
+
+                var result = ErrorRatioCalc.GetErrorRatioResult(_model.ModelID, porosity, sw, depth, true, calcSpeed, targetErrorRatio / 100);
+                Plot(result);
 
             }
             else
@@ -171,27 +194,41 @@
             {
                 case "cmbBit":
                     cmbCasingOD.Items.Clear();
+                    _filters.Remove("CasingOD");
                     goto case "cmbCasingOD";
                 case "cmbCasingOD":
                     cmbCasingID.Items.Clear();
+                    _filters.Remove("CasingID");
                     goto case "cmbCasingID";
                 case "cmbCasingID":
                     cmbScreenOD.Items.Clear();
+                    _filters.Remove("ScreenOD");
                     goto case "cmbScreenOD";
                 case "cmbScreenOD":
                     cmbScreenID.Items.Clear();
+                    _filters.Remove("ScreenID");
                     goto case "cmbScreenID";
                 case "cmbScreenID":
                     cmbTubeOD.Items.Clear();
+                    _filters.Remove("TubeOD");
                     goto case "cmbTubeOD";
                 case "cmbTubeOD":
                     cmbTubeID.Items.Clear();
+                    _filters.Remove("TubeID");
                     break;
                 case "cmbTubeID":
                     break;
                 default:
                     break;
             }
+            cmbCasingHoldUp.Items.Clear();
+            _filters.Remove("CasingFluid");
+            cmbScreenHoldUp.Items.Clear();
+            _filters.Remove("ScreenFluid");
+            cmbTubeHoldUp.Items.Clear();
+            _filters.Remove("TubeFluid");
+            cmbLithology.Items.Clear();
+            _filters.Remove("Lithology");
 
             var unSelectedCmbs = allCmb.Where(c => c.SelectedIndex == -1).ToList();
             var curModels = GetFilterResult();
@@ -206,8 +243,10 @@
                     cmb.Enabled = false;
 
                 cmb.Items.AddRange(currentItems);
+
                 if (currentItems.Length == 1)
                     cmb.Text = currentItems[0];
+
             }
 
             LinkageModelProperty();
@@ -226,12 +265,44 @@
             else
                 _filters.Add(property, filterVal);
 
-            var curModels = GetFilterResult();
-
             var allCmb = parentControl.Controls.OfType<Control>().Where(c => c is ComboBox).Select(c => c as ComboBox).ToList();
-
+            switch (control.Name)
+            {
+                case "cmbCasingHoldUp":
+                    cmbScreenHoldUp.Items.Clear();
+                    _filters.Remove("ScreenFluid");
+                    goto case "cmbScreenHoldUp";
+                case "cmbScreenHoldUp":
+                    cmbTubeHoldUp.Items.Clear();
+                    _filters.Remove("TubeFluid");
+                    break;
+                default:
+                    break;
+            }
+            var unSelectedCmbs = allCmb.Where(c => c.SelectedIndex == -1).ToList();
+            var curModels = GetFilterResult();
             var modelType = typeof(COModelTable);
 
+            foreach (var cmb in unSelectedCmbs)
+            {
+                cmb.Items.Clear();
+                cmb.Enabled = true;
+                var cmbProperty = modelType.GetProperty(cmb.Tag.ToString());
+                var currentItems = curModels.Select(c => cmbProperty.GetValue(c).ToString() == "0" ? "water" : cmbProperty.GetValue(c).ToString() == "100" ? "oil" : cmbProperty.GetValue(c).ToString() == "200" ? "gas" : "none").Distinct().ToArray();
+                if (currentItems.Length == 0)
+                    cmb.Enabled = false;
+
+                cmb.Items.AddRange(currentItems);
+                if (_filters.ContainsKey(cmb.Tag.ToString()) && currentItems.Contains(_filters[cmb.Tag.ToString()]))
+                    cmb.Text = _filters[cmb.Tag.ToString()];
+                else
+                {
+                    _filters.Remove(cmb.Tag.ToString());
+                    if (currentItems.Length == 1)
+                        cmb.SelectedIndex = 0;
+                }
+
+            }
             //foreach (var cmb in allCmb)
             //{
             //    if (cmb.Name != control.Name)
@@ -249,17 +320,25 @@
             //}
 
             var lithology = curModels.Select(x => x.Lithology).Distinct().ToArray();
-            var curLithology = cmbLithology.Text;
             cmbLithology.Items.Clear();
             cmbLithology.Items.AddRange(lithology);
-            cmbLithology.Text = curLithology;
 
-            if (lithology.Length == 1)
-                cmbLithology.Text = lithology[0];
+            if (_filters.ContainsKey("Lithology") && lithology.Contains(_filters["Lithology"]))
+                cmbLithology.Text = _filters["Lithology"];
+            else
+                _filters.Remove("Lithology");
         }
 
         private void LinkageModelProperty()
         {
+            cmbCasingHoldUp.Items.Clear();
+            _filters.Remove("CasingFluid");
+            cmbScreenHoldUp.Items.Clear();
+            _filters.Remove("ScreenFluid");
+            cmbTubeHoldUp.Items.Clear();
+            _filters.Remove("TubeFluid");
+            cmbLithology.Items.Clear();
+            _filters.Remove("Lithology");
             var curModels = GetFilterResult();
             if (curModels.Count == 0)
                 return;
@@ -290,9 +369,9 @@
 
 
             var casingFulids = curModels.Select(x => x.CasingFluid == "0" ? "water" : x.CasingFluid == "100" ? "oil" : x.CasingFluid == "200" ? "gas" : "none").Distinct().ToArray();
+
             if (casingFulids.Length == 1 && casingFulids[0] == "999")
             {
-                cmbCasingHoldUp.Items.Clear();
                 cmbCasingHoldUp.Enabled = false;
                 nudCasingHoldUp.Enabled = false;
             }
@@ -301,14 +380,13 @@
                 cmbCasingHoldUp.Enabled = true;
                 nudCasingHoldUp.Enabled = true;
 
-                cmbCasingHoldUp.Items.Clear();
                 cmbCasingHoldUp.Items.AddRange(casingFulids);
             }
 
             var screenFulids = curModels.Select(x => x.ScreenFluid == "0" ? "water" : x.ScreenFluid == "100" ? "oil" : x.ScreenFluid == "200" ? "gas" : "none").Distinct().ToArray();
+
             if (screenFulids.Length == 1 && screenFulids[0] == "none")
             {
-                cmbScreenHoldUp.Items.Clear();
                 cmbScreenHoldUp.Enabled = false;
                 nudScreenHoldUp.Enabled = false;
             }
@@ -316,15 +394,12 @@
             {
                 cmbScreenHoldUp.Enabled = true;
                 nudScreenHoldUp.Enabled = true;
-                cmbScreenHoldUp.Items.Clear();
                 cmbScreenHoldUp.Items.AddRange(screenFulids);
             }
 
             var tubeFluids = curModels.Select(x => x.TubeFluid == "0" ? "water" : x.TubeFluid == "100" ? "oil" : x.TubeFluid == "200" ? "gas" : "none").Distinct().ToArray();
             if (tubeFluids.Length == 1 && tubeFluids[0] == "none")
             {
-
-                cmbTubeHoldUp.Items.Clear();
                 cmbTubeHoldUp.Enabled = false;
                 nudTubeHoldUp.Enabled = false;
             }
@@ -332,19 +407,16 @@
             {
                 cmbTubeHoldUp.Enabled = true;
                 nudTubeHoldUp.Enabled = true;
-
-                cmbTubeHoldUp.Items.Clear();
                 cmbTubeHoldUp.Items.AddRange(tubeFluids);
             }
 
             var lithology = curModels.Select(x => x.Lithology).Distinct().ToArray();
             var curLithology = cmbLithology.Text;
-            cmbLithology.Items.Clear();
             cmbLithology.Items.AddRange(lithology);
             cmbLithology.Text = curLithology;
 
             if (lithology.Length == 1)
-                cmbLithology.Text = lithology[0];
+                cmbLithology.SelectedIndex = 0;
         }
 
         private List<COModelTable> GetFilterResult()
@@ -369,5 +441,45 @@
             else
                 _filters.Add(property, filterVal);
         }
+
+        private void Plot(ErrorRatioResult errorResult)
+        {
+            pnlCOPic.Controls.Clear();
+            var plotModel = new PlotModel();
+
+            var xAxis = new LinearAxis { Position = AxisPosition.Bottom, Title = "Porosity(%)" };
+            var yAxis = new LinearAxis { Position = AxisPosition.Left, Title = "C/O" };
+
+            plotModel.Axes.Add(xAxis);
+            plotModel.Axes.Add(yAxis);
+
+            // 创建曲线系列
+            var oSeries = new LineSeries { Title = "Water Line", Color = OxyColors.MediumBlue };
+            var cSeries = new LineSeries { Title = "Oil Line", Color = OxyColors.DarkRed };
+
+            foreach (var o in errorResult.OWOL)
+                oSeries.Points.Add(new DataPoint(o[0], o[1]));
+            foreach (var c in errorResult.CWOL)
+                cSeries.Points.Add(new DataPoint(c[0], c[1]));
+
+            plotModel.Legends.Add(new Legend { LegendPosition = LegendPosition.TopLeft });
+            // 添加曲线到图表模型
+            plotModel.Series.Add(oSeries);
+            plotModel.Series.Add(cSeries);
+
+            // 可视化代码
+            // 创建 Windows 窗体以显示图表
+            var plotView = new PlotView { Model = plotModel };
+            
+            plotView.Width = 600;
+            plotView.Height = 450;
+            var bitmap = new Bitmap(plotView.Width, plotView.Height);
+            plotView.DrawToBitmap(bitmap, new Rectangle(0, 0, bitmap.Width, bitmap.Height));
+            bitmap.Save(AppDomain.CurrentDomain.BaseDirectory + "co.png", System.Drawing.Imaging.ImageFormat.Png);
+
+            plotView.Dock = DockStyle.Fill;
+            pnlCOPic.Controls.Add(plotView);
+            
+        }
     }
 }

--
Gitblit v1.9.3