From e0279fe6db39d1071ef04ad0080b887cb6a0d335 Mon Sep 17 00:00:00 2001
From: lx <ex_lixiang17@cosl.com.cn>
Date: 星期一, 04 八月 2025 17:27:47 +0800
Subject: [PATCH] upload
---
ErrorAnalysis.Service/COMergeCalcService.cs | 191 ++++++++++++++++++++++++++++++++++-------------
1 files changed, 138 insertions(+), 53 deletions(-)
diff --git a/ErrorAnalysis.Service/COMergeCalcService.cs b/ErrorAnalysis.Service/COMergeCalcService.cs
index 06861c2..08c53d8 100644
--- a/ErrorAnalysis.Service/COMergeCalcService.cs
+++ b/ErrorAnalysis.Service/COMergeCalcService.cs
@@ -4,11 +4,15 @@
using System.Linq;
using System.Text;
using System.Threading.Tasks;
+using static System.Runtime.InteropServices.JavaScript.JSType;
namespace ErrorAnalysis.Service
{
public class COMergeCalcService
{
+ //Ncoef=1e8*2.54*2.54*2*6*11*2*77*0.2/0.6*speed;
+ //Fcoef=1e8*2.54*2.54*4*6*2*11*77*0.2/0.6*speed;
+
const double _gg = (9.0 - 0.0) / 255.0;//0.03529
const double _offset = 0.0;
@@ -18,69 +22,150 @@
const int _oWinStartIndex = (int)(4.88 / _gg + _offset); // org
const int _oWinEndIndex = (int)(6.36 / _gg + _offset); // width = 6.36-4.88=1.48
- public static double GetFarMergeCResult(string connectionString, string modelId, int porosity)
+ const double _coef = 1e8 * 2.54 * 2.54 * 4 * 6 * 2 * 11 * 77 * 0.2 / 0.6;
+
+ private static (double, double) GetFarInterplolateResult(string modelId, double porosity, double sw, double speed, double depth)
{
- var repository = new COFarResultRepository(connectionString);
- var cResult = repository.GetCOFarResult(modelId, porosity, 0);
- if (cResult == null)
- throw new InvalidDataException("COFarResult not found");
- var cArr = cResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray();
- if (cArr == null)
- throw new InvalidDataException("COFarResult InelasticSpec is null");
+ var cWinLength = _cWinEndIndex - _cWinStartIndex + 1;
+ var oWinLength = _oWinEndIndex - _oWinStartIndex + 1;
+ var oilLine = RepositoryInstance.Instance.COFarResultRepository?.GetCOFarResult(modelId, (int)porosity, 0);
+ var waterLine = RepositoryInstance.Instance.COFarResultRepository?.GetCOFarResult(modelId, (int)porosity, 100);
+ if (oilLine == null || waterLine == null)
+ throw new InvalidDataException("COFarResult Line is null");
+ var oilLineValArr = oilLine.InelasticSpec.Split(',').Select(v => Convert.ToDouble(v)).ToArray();
+ var waterLineValArr = waterLine.InelasticSpec.Split(',').Select(v => Convert.ToDouble(v)).ToArray();
- var length = _cWinEndIndex - _cWinStartIndex + 1;
- var result = cArr.Skip(_cWinStartIndex).Take(length).Sum();
+ var oilLineC = oilLineValArr.Skip(_cWinStartIndex).Take(cWinLength).Sum() * _coef / speed * depth * 0.07;
+ var oilLineO = oilLineValArr.Skip(_oWinStartIndex).Take(oWinLength).Sum() * _coef / speed * depth * 0.07;
+ var waterLineC = waterLineValArr.Skip(_cWinStartIndex).Take(cWinLength).Sum() * _coef / speed * depth * 0.07;
+ var waterLineO = waterLineValArr.Skip(_oWinStartIndex).Take(oWinLength).Sum() * _coef / speed * depth * 0.07;
+
+ double cRes = 0, oRes = 0;
+ if (sw == 0)
+ {
+ cRes = waterLineC;
+ oRes = waterLineO;
+ }
+ else if (sw == 100)
+ {
+ cRes = waterLineC;
+ oRes = waterLineO;
+ }
+ else if (sw > 0 && sw < 100)
+ {
+ cRes = Utility.Interpolate(sw, 100, waterLineC, 0, oilLineC);
+ oRes = Utility.Interpolate(sw, 0, oilLineO, 100, waterLineO);
+ }
+ return (cRes, oRes);
+ }
+
+ private static (double, double) GetNearInterplolateResult(string modelId, double porosity, double sw, double speed, double depth)
+ {
+ var cWinLength = _cWinEndIndex - _cWinStartIndex + 1;
+ var oWinLength = _oWinEndIndex - _oWinStartIndex + 1;
+ var oilLine = RepositoryInstance.Instance.CONearResultRepository?.GetCONearResult(modelId, (int)porosity, 0);
+ var waterLine = RepositoryInstance.Instance.CONearResultRepository?.GetCONearResult(modelId, (int)porosity, 100);
+ if (oilLine == null || waterLine == null)
+ throw new InvalidDataException("COFarResult Line is null");
+ var oilLineValArr = oilLine.InelasticSpec.Split(',').Select(v => Convert.ToDouble(v)).ToArray();
+ var waterLineValArr = waterLine.InelasticSpec.Split(',').Select(v => Convert.ToDouble(v)).ToArray();
+
+ var oilLineC = oilLineValArr.Skip(_cWinStartIndex).Take(cWinLength).Sum() * _coef / speed * depth * 0.07;
+ var oilLineO = oilLineValArr.Skip(_oWinStartIndex).Take(oWinLength).Sum() * _coef / speed * depth * 0.07;
+
+ var waterLineC = waterLineValArr.Skip(_cWinStartIndex).Take(cWinLength).Sum() * _coef / speed * depth * 0.07;
+ var waterLineO = waterLineValArr.Skip(_oWinStartIndex).Take(oWinLength).Sum() * _coef / speed * depth * 0.07;
+
+ double cRes = 0, oRes = 0;
+ if (sw == 0)
+ {
+ cRes = waterLineC;
+ oRes = waterLineO;
+ }
+ else if (sw == 100)
+ {
+ cRes = waterLineC;
+ oRes = waterLineO;
+ }
+ else if (sw > 0 && sw < 100)
+ {
+ cRes = Utility.Interpolate(sw, 100, waterLineC, 0, oilLineC);
+ oRes = Utility.Interpolate(sw, 0, oilLineO, 100, waterLineO);
+ }
+ return (cRes, oRes);
+ }
+
+ /// <summary>
+ /// 鑾峰彇杩滄帰澶碈鎴朞鍊�
+ /// </summary>
+ /// <param name="connectionString">鏁版嵁搴撹繛鎺ュ瓧绗︿覆</param>
+ /// <param name="modelId">绠℃煴ID</param>
+ /// <param name="porosity">瀛旈殭搴�</param>
+ /// <param name="sw">鍚按楗卞拰搴�</param>
+ /// <param name="readSw">璇诲彇鍚按楗卞拰搴︼紙0涓虹⒊锛�100涓烘哀锛�</param>
+ /// <param name="speed">娴嬮��</param>
+ /// <param name="depth">娣卞害</param>
+ /// <returns>纰冲拰姘у��</returns>
+ /// <exception cref="InvalidDataException">瀛旈殭搴﹁秴杩囪寖鍥�</exception>
+ public static (double, double) GetFarMergeCOResult(string modelId, double porosity, double sw, double speed, double depth)
+ {
+ if (porosity > 40)
+ throw new InvalidDataException("Porosity value out of range!");
+
+ (double, double) result = (0, 0);
+ if (porosity % 5 == 0)
+ {
+ result = GetFarInterplolateResult(modelId, porosity, sw, speed, depth);
+ }
+ else
+ {
+ var ceilingPorosity = Math.Ceiling(porosity / 5) * 5;
+ var ceilingResult = GetFarInterplolateResult(modelId, ceilingPorosity, sw, speed, depth);
+
+ var floorPorosity = Math.Floor(porosity / 5) * 5;
+ var floorResult = GetFarInterplolateResult(modelId, floorPorosity, sw, speed, depth);
+
+ result.Item1 = Utility.Interpolate(porosity, floorPorosity, floorResult.Item1, ceilingPorosity, ceilingResult.Item1);
+ result.Item2 = Utility.Interpolate(porosity, floorPorosity, floorResult.Item2, ceilingPorosity, ceilingResult.Item2);
+ }
return result;
}
- public static double GetFarMergeOResult(string connectionString, string modelId, int porosity)
+ /// <summary>
+ /// 鑾峰彇杩戞帰澶碈鎴朞鍊�
+ /// </summary>
+ /// <param name="connectionString">鏁版嵁搴撹繛鎺ュ瓧绗︿覆</param>
+ /// <param name="modelId">绠℃煴ID</param>
+ /// <param name="porosity">瀛旈殭搴�</param>
+ /// <param name="sw">鍚按楗卞拰搴�</param>
+ /// <param name="readSw">璇诲彇鍚按楗卞拰搴︼紙0涓虹⒊锛�100涓烘哀锛�</param>
+ /// <param name="speed">娴嬮��</param>
+ /// <param name="depth">娣卞害</param>
+ /// <returns>纰虫垨姘у��</returns>
+ /// <exception cref="InvalidDataException">瀛旈殭搴﹁秴杩囪寖鍥�</exception>
+ public static (double, double) GetNearMergeCOResult(string modelId, double porosity, double sw, double speed, double depth)
{
- var repository = new COFarResultRepository(connectionString);
- var cResult = repository.GetCOFarResult(modelId, porosity, 100);
- if (cResult == null)
- throw new InvalidDataException("COFarResult not found");
- var cArr = cResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray();
- if (cArr == null)
- throw new InvalidDataException("COFarResult InelasticSpec is null");
+ if (porosity > 40)
+ throw new InvalidDataException("Porosity value out of range!");
- var length = _oWinEndIndex - _oWinStartIndex + 1;
- var result = cArr.Skip(_oWinStartIndex).Take(length).Sum();
+ (double, double) result = (0, 0);
+ if (porosity % 5 == 0)
+ {
+ result = GetNearInterplolateResult(modelId, porosity, sw, speed, depth);
+ }
+ else
+ {
+ var ceilingPorosity = Math.Ceiling(porosity / 5) * 5;
+ var ceilingResult = GetNearInterplolateResult(modelId, ceilingPorosity, sw, speed, depth);
+ var floorPorosity = Math.Floor(porosity / 5) * 5;
+ var floorResult = GetNearInterplolateResult(modelId, floorPorosity, sw, speed, depth);
+
+ result.Item1 = Utility.Interpolate(porosity, floorPorosity, floorResult.Item1, ceilingPorosity, ceilingResult.Item1);
+ result.Item2 = Utility.Interpolate(porosity, floorPorosity, floorResult.Item2, ceilingPorosity, ceilingResult.Item2);
+ }
return result;
}
-
- public static double GetNearMergeCResult(string connectionString, string modelId, int porosity)
- {
- var repository = new CONearResultRepository(connectionString);
- var cResult = repository.GetCONearResult(modelId, porosity, 0);
- if (cResult == null)
- throw new InvalidDataException("CONearResult not found");
- var cArr = cResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray();
- if (cArr == null)
- throw new InvalidDataException("CONearResult InelasticSpec is null");
-
- var length = _cWinEndIndex - _cWinStartIndex + 1;
- var result = cArr.Skip(_cWinStartIndex).Take(length).Sum();
-
- return result;
- }
-
- public static double GetNearMergeOResult(string connectionString, string modelId, int porosity)
- {
- var repository = new CONearResultRepository(connectionString);
- var cResult = repository.GetCONearResult(modelId, porosity, 100);
- if (cResult == null)
- throw new InvalidDataException("CONearResult not found");
- var cArr = cResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray();
- if (cArr == null)
- throw new InvalidDataException("CONearResult InelasticSpec is null");
-
- var length = _oWinEndIndex - _oWinStartIndex + 1;
- var result = cArr.Skip(_oWinStartIndex).Take(length).Sum();
-
- return result;
- }
-
}
}
--
Gitblit v1.9.3