1.增加压裂井裂缝半长拟合参数 2.参数物理边界统一用框架中的

Bin/Config/Lang/cn/nmNum_cn.ts

@@ -3490,6 +3490,10 @@ Supported types: Vertical, Vertical Fractured, and Horizontal Multi-Fractured We
         <source>Fracture conductivity</source>
         <translation>裂缝导流能力</translation>
     </message>
+    <message>
+        <source>Fracture half length</source>
+        <translation>裂缝半长</translation>
+    </message>
     <message>
         <source>The minimum value of %1 must be greater than zero for automatic fitting.</source>
         <translation>自动拟合时，%1 的最小值和初始值必须大于零。</translation>

Include/nmNum/nmCalculation/nmCalculationAutoFitPSO.h

@@ -372,12 +372,13 @@ private:
 	// 参数索引约定：
 	// 0 k 渗透率；1 skin 表皮系数；2 wellboreC 井筒储集；
 	// 3 phi 孔隙度；4 h 储层厚度；5 Ct 综合压缩系数；
-	// 6 Cf 岩石压缩系数；7 Swi 初始含水饱和度。
+	// 6 Cf 岩石压缩系数；7 Swi 初始含水饱和度；
+	// 8 Dfc 裂缝导流能力；9 fractureHalfLength 裂缝半长。
 	// m_enabledParamIndices 保存被用户勾选的参数索引，粒子的 position 维度与它一致。
-	QVector<bool> m_parameterSelected;      // 完整 9 个参数是否被用户勾选参与拟合。
-	QVector<double> m_parameterLower;       // 完整 9 个参数的搜索下界。
-	QVector<double> m_parameterUpper;       // 完整 9 个参数的搜索上界。
-	QVector<int> m_enabledParamIndices;     // 被勾选参数在完整 9 维体系中的索引。
+	QVector<bool> m_parameterSelected;      // 完整 10 个参数是否被用户勾选参与拟合。
+	QVector<double> m_parameterLower;       // 完整 10 个参数的搜索下界。
+	QVector<double> m_parameterUpper;       // 完整 10 个参数的搜索上界。
+	QVector<int> m_enabledParamIndices;     // 被勾选参数在完整 10 维体系中的索引。
 	QVector<QVector<double> > m_targetLogLogData; // 目标井 history log-log 曲线：time/pressure/derivative。
 	QString m_targetWellName;               // 目标井名称；读写井参数和读取模拟曲线都依赖它。
 

Include/nmNum/nmData/nmDataAutomaticFitting.h

@@ -82,6 +82,13 @@ public:
 	nmDataAttribute& getFractureConductivityMin();
 	void setFractureConductivityMin(const nmDataAttribute& fractureConductivityMin);
 
+	// Getter and Setter for fractureHalfLengthMax
+	nmDataAttribute& getFractureHalfLengthMax();
+	void setFractureHalfLengthMax(const nmDataAttribute& fractureHalfLengthMax);
+	// Getter and Setter for fractureHalfLengthMin
+	nmDataAttribute& getFractureHalfLengthMin();
+	void setFractureHalfLengthMin(const nmDataAttribute& fractureHalfLengthMin);
+
 	// Getter and Setter for iteration count
 	nmDataAttribute& getIterationCount();
 	void setIterationCount(const nmDataAttribute& iterationCount);
@@ -124,6 +131,9 @@ public:
 	bool getFractureConductivitySelected() const;
 	void setFractureConductivitySelected(bool selected);
 
+	bool getFractureHalfLengthSelected() const;
+	void setFractureHalfLengthSelected(bool selected);
+
 private:
 	// 参数最大值
 	nmDataAttribute m_permeabilityMax;
@@ -135,6 +145,7 @@ private:
 	nmDataAttribute m_cfMax;
 	nmDataAttribute m_swiMax;
 	nmDataAttribute m_fractureConductivityMax;
+	nmDataAttribute m_fractureHalfLengthMax;
 
 	// 参数最小值
 	nmDataAttribute m_permeabilityMin;
@@ -146,6 +157,7 @@ private:
 	nmDataAttribute m_cfMin;
 	nmDataAttribute m_swiMin;
 	nmDataAttribute m_fractureConductivityMin;
+	nmDataAttribute m_fractureHalfLengthMin;
 
 	// 迭代参数
 	nmDataAttribute m_iterationCount;      // 迭代步数
@@ -163,6 +175,7 @@ private:
 	bool m_cfSelected;                     // 是否选择岩石压缩系数进行拟合
     bool m_swiSelected;					   // 是否选择初始含水饱和度进行拟合
 	bool m_fractureConductivitySelected; // 是否选择裂缝导流能力进行拟合
+	bool m_fractureHalfLengthSelected;  // 是否选择裂缝半长进行拟合
 };
 
 #endif // NMDATAAUTOMATICFITTING_H

Include/nmNum/nmSubWxs/nmWxAutomaticFitting.h

@@ -92,6 +92,7 @@ private:
 	QCheckBox* m_cfCheckBox;       // 岩石压缩系数
 	QCheckBox* m_swiCheckBox;	   // 初始含水饱和度
 	QCheckBox* m_dfcCheckBox;      // 裂缝导流能力
+	QCheckBox* m_fractureHalfLengthCheckBox; // 裂缝半长
 
 	// 按钮
 	QPushButton* m_reverseBtn;

Src/nmNum/nmCalculation/nmCalculationAutoFitPSO.cpp

@@ -432,7 +432,7 @@ static QString findExecutableInPath(const QString& executableName)
 static QStringList traceParameterNames()
 {
     // trace 和 trace meta 使用的完整参数名顺序。
-    // 这个顺序必须与 buildTraceParameterVector() 和 m_parameterSelected 的 0-8 索引一致。
+    // 这个顺序必须与 buildTraceParameterVector() 和 m_parameterSelected 的 0-9 索引一致。
     QStringList names;
     names << "k"
           << "skin"
@@ -442,7 +442,8 @@ static QStringList traceParameterNames()
           << "Ct"
           << "Cf"
           << "Swi"
-          << "Dfc";
+          << "Dfc"
+          << "fractureHalfLength";
     return names;
 }
 
@@ -1013,7 +1014,7 @@ void nmCalculationAutoFitPSO::writeTraceHeader()
 {
     // trace CSV 字段说明：
     // - 代理模式保持原 k/skin/wellboreC/phi/h/Ct/Cf 契约；
-    // - 非代理模式额外记录 Swi/Dfc，便于复盘信赖域对两项参数的调整；
+    // - 非代理模式额外记录 Swi/Dfc/裂缝半长，便于复盘信赖域调整；
     // - solver_objective 是真实求解器误差；
     // - surrogate_objective 是 Python 代理评分；
     // - screening_decision 说明该粒子为什么跑/不跑真实求解器；
@@ -1037,7 +1038,8 @@ void nmCalculationAutoFitPSO::writeTraceHeader()
          << "Cf";
     if(!isSurrogateScreeningEnabled()) {
         cols << "Swi"
-             << "Dfc";
+             << "Dfc"
+             << "fractureHalfLength";
     }
     cols << "solver_objective"
          << "solver_success"
@@ -1054,7 +1056,8 @@ void nmCalculationAutoFitPSO::writeTraceHeader()
          << "pbest_Cf";
     if(!isSurrogateScreeningEnabled()) {
         cols << "pbest_Swi"
-             << "pbest_Dfc";
+             << "pbest_Dfc"
+             << "pbest_fractureHalfLength";
     }
     cols << "gbest_objective"
          << "gbest_k"
@@ -1066,7 +1069,8 @@ void nmCalculationAutoFitPSO::writeTraceHeader()
          << "gbest_Cf";
     if(!isSurrogateScreeningEnabled()) {
         cols << "gbest_Swi"
-             << "gbest_Dfc";
+             << "gbest_Dfc"
+             << "gbest_fractureHalfLength";
     }
     cols << "enabled_param_indices"
          << "pressure_loss"
@@ -1132,7 +1136,7 @@ void nmCalculationAutoFitPSO::writeTraceMetaFile()
     // 非代理 v5 增加带符号诊断列；代理 PSO 保留原 v3 字段和目标，避免改变
     // 已有模型的训练和回放契约。
     out << "  \"schema_version\": "
-        << (isSurrogateScreeningEnabled() ? 3 : 6) << ",\n";
+        << (isSurrogateScreeningEnabled() ? 3 : 7) << ",\n";
     out << "  \"trace_type\": "
         << jsonEscape(isSurrogateScreeningEnabled()
                       ? "pso_baseline_replay_meta"
@@ -1193,10 +1197,10 @@ void nmCalculationAutoFitPSO::writeTraceMetaFile()
 
 QVector<double> nmCalculationAutoFitPSO::buildTraceParameterVector(const QVector<double>& selectedParameters) const
 {
-    // 将粒子内部使用的“启用参数向量”还原成完整 9 维参数向量。
+    // 将粒子内部使用的“启用参数向量”还原成完整 10 维参数向量。
     // 未启用的参数从当前 DataManager 读取，启用的参数用 selectedParameters 覆盖。
     // trace CSV、候选 CSV、代理训练域检查都需要这个完整向量。
-    QVector<double> fullParams(9, 0.0);
+    QVector<double> fullParams(10, 0.0);
 
     nmDataAnalyzeManager* dataManager = nmDataAnalyzeManager::getCurrentInstance();
 
@@ -1224,12 +1228,14 @@ QVector<double> nmCalculationAutoFitPSO::buildTraceParameterVector(const QVector
                     dynamic_cast<nmDataVerticalFracturedWell*>(pTargetWell);
                 if(fracturedWell) {
                     fullParams[8] = fracturedWell->getDfc().getValue().toDouble();
+                    fullParams[9] = fracturedWell->getFractureHalfLength().getValue().toDouble();
                 }
             } else if(pTargetWell->getWellType() == NM_WELL_MODEL::Horizontal_Fractured_Well) {
                 nmDataHorizontalFracturedWell* fracturedWell =
                     dynamic_cast<nmDataHorizontalFracturedWell*>(pTargetWell);
                 if(fracturedWell) {
                     fullParams[8] = fracturedWell->getDfc().getValue().toDouble();
+                    fullParams[9] = fracturedWell->getFractureHalfLength().getValue().toDouble();
                 }
             }
         }
@@ -1290,7 +1296,8 @@ void nmCalculationAutoFitPSO::writeTraceRow(int generation,
          << traceParamAt(currentParams, 6);
     if(!isSurrogateScreeningEnabled()) {
         cols << traceParamAt(currentParams, 7)
-             << traceParamAt(currentParams, 8);
+             << traceParamAt(currentParams, 8)
+             << traceParamAt(currentParams, 9);
     }
     cols << traceNumber(solverObjective)
          << QString::number(solverSuccess ? 1 : 0)
@@ -1307,7 +1314,8 @@ void nmCalculationAutoFitPSO::writeTraceRow(int generation,
          << traceParamAt(pbestParams, 6);
     if(!isSurrogateScreeningEnabled()) {
         cols << traceParamAt(pbestParams, 7)
-             << traceParamAt(pbestParams, 8);
+             << traceParamAt(pbestParams, 8)
+             << traceParamAt(pbestParams, 9);
     }
     cols << traceNumber(m_globalBestFitness)
          << traceParamAt(gbestParams, 0)
@@ -1319,7 +1327,8 @@ void nmCalculationAutoFitPSO::writeTraceRow(int generation,
          << traceParamAt(gbestParams, 6);
     if(!isSurrogateScreeningEnabled()) {
         cols << traceParamAt(gbestParams, 7)
-             << traceParamAt(gbestParams, 8);
+             << traceParamAt(gbestParams, 8)
+             << traceParamAt(gbestParams, 9);
     }
     cols << csvEscape(enabledIndices.join(";"));
 
@@ -2917,14 +2926,14 @@ void nmCalculationAutoFitPSO::loadParameterBounds()
     // 读取用户勾选的拟合参数及上下界。
     //
     // 这里构建三个核心数组：
-    // - m_parameterSelected[9]：完整参数体系中每个参数是否参与拟合；
-    // - m_parameterLower/Upper[9]：完整参数体系的搜索上下界；
+    // - m_parameterSelected[10]：完整参数体系中每个参数是否参与拟合；
+    // - m_parameterLower/Upper[10]：完整参数体系的搜索上下界；
     // - m_enabledParamIndices：把粒子内部紧凑向量映射回完整参数索引。
     nmDataAnalyzeManager* dataManager = nmDataAnalyzeManager::getCurrentInstance();
     nmDataAutomaticFitting fittingData = dataManager->getAutomaticFittingDataCopy();
 
     // 获取参数选择状态
-    m_parameterSelected.resize(9);
+    m_parameterSelected.resize(10);
     m_parameterSelected[0] = fittingData.getPermeabilitySelected();
     m_parameterSelected[1] = fittingData.getSkinSelected();
     m_parameterSelected[2] = fittingData.getWellboreStorageSelected();
@@ -2934,10 +2943,11 @@ void nmCalculationAutoFitPSO::loadParameterBounds()
     m_parameterSelected[6] = fittingData.getCfSelected();
     m_parameterSelected[7] = fittingData.getSwiSelected();
     m_parameterSelected[8] = fittingData.getFractureConductivitySelected();
+    m_parameterSelected[9] = fittingData.getFractureHalfLengthSelected();
 
     // 获取参数边界
-    m_parameterLower.resize(9);
-    m_parameterUpper.resize(9);
+    m_parameterLower.resize(10);
+    m_parameterUpper.resize(10);
 
     m_parameterLower[0] = fittingData.getPermeabilityMin().getValue().toDouble();
     m_parameterUpper[0] = fittingData.getPermeabilityMax().getValue().toDouble();
@@ -2966,6 +2976,9 @@ void nmCalculationAutoFitPSO::loadParameterBounds()
     m_parameterLower[8] = fittingData.getFractureConductivityMin().getValue().toDouble();
     m_parameterUpper[8] = fittingData.getFractureConductivityMax().getValue().toDouble();
 
+    m_parameterLower[9] = fittingData.getFractureHalfLengthMin().getValue().toDouble();
+    m_parameterUpper[9] = fittingData.getFractureHalfLengthMax().getValue().toDouble();
+
     // 更新启用参数索引
     m_enabledParamIndices.clear();
 
@@ -3474,11 +3487,14 @@ bool nmCalculationAutoFitPSO::startAutoFitting()
             applyParametersToDataManager(m_globalBestPosition);
 
 			// 即使用户此时停止、不再执行最终完整计算，也要把 PEBI 缓存恢复为
-			// 最终已接受的 Dfc，避免缓存仍停留在最后一个被拒绝的候选值。
-			if(m_parameterSelected.size() > 8 && m_parameterSelected[8]) {
+			// 最终已接受的裂缝参数，避免缓存仍停留在最后一个被拒绝的候选值。
+			const bool fractureGridParameterSelected =
+				(m_parameterSelected.size() > 8 && m_parameterSelected[8]) ||
+				(m_parameterSelected.size() > 9 && m_parameterSelected[9]);
+			if(fractureGridParameterSelected) {
 				nmCalculationPebiGrid* pebiGrid = nmCalculationPebiGrid::getInstance();
 				if(!pebiGrid || !pebiGrid->generateOutputPara()) {
-					throw std::runtime_error("Failed to refresh final fracture conductivity");
+					throw std::runtime_error("Failed to refresh final fracture parameters");
 				}
 			}
 
@@ -3711,6 +3727,22 @@ void nmCalculationAutoFitPSO::extractUserInitialValues()
                     }
                 }
                 break;
+
+            case 9: // 裂缝半长
+                if(pTargetWell && pTargetWell->getWellType() == NM_WELL_MODEL::Vertical_Fractured_Well) {
+                    nmDataVerticalFracturedWell* fracturedWell =
+                        dynamic_cast<nmDataVerticalFracturedWell*>(pTargetWell);
+                    if(fracturedWell) {
+                        initialValue = fracturedWell->getFractureHalfLength().getValue().toDouble();
+                    }
+                } else if(pTargetWell && pTargetWell->getWellType() == NM_WELL_MODEL::Horizontal_Fractured_Well) {
+                    nmDataHorizontalFracturedWell* fracturedWell =
+                        dynamic_cast<nmDataHorizontalFracturedWell*>(pTargetWell);
+                    if(fracturedWell) {
+                        initialValue = fracturedWell->getFractureHalfLength().getValue().toDouble();
+                    }
+                }
+                break;
         }
 
         m_initialValues.append(initialValue);
@@ -5534,12 +5566,15 @@ double nmCalculationAutoFitPSO::evaluateFitness(const QVector<double>& parameter
             return 1e10;
         }
 
-		// Dfc 位于 PEBI 裂缝数组 crack[5]，不是每次求解都会重新组装的 Base/CS 参数。
-		// 因此只有勾选 Dfc 时才刷新一次网格输出参数，保证本次真实试算使用刚写入井对象的值。
-		if(m_parameterSelected.size() > 8 && m_parameterSelected[8]) {
+		// Dfc 和裂缝半长都通过 PEBI 裂缝数组传入，不属于每次求解都会重新组装的 Base/CS 参数。
+		// 勾选任一裂缝参数时刷新网格输出，保证本次真实试算使用新的导流能力和端点坐标。
+		const bool fractureGridParameterSelected =
+			(m_parameterSelected.size() > 8 && m_parameterSelected[8]) ||
+			(m_parameterSelected.size() > 9 && m_parameterSelected[9]);
+		if(fractureGridParameterSelected) {
 			nmCalculationPebiGrid* pebiGrid = nmCalculationPebiGrid::getInstance();
 			if(!pebiGrid || !pebiGrid->generateOutputPara()) {
-				DEBUG_OUT(QString("%1: Call #%2 - Failed to refresh PEBI fracture conductivity")
+				DEBUG_OUT(QString("%1: Call #%2 - Failed to refresh PEBI fracture parameters")
 				          .arg(funcName).arg(callCount));
 				return 1e10;
 			}
@@ -5707,7 +5742,7 @@ void nmCalculationAutoFitPSO::updateReservoirParameters(const QVector<double>& p
     nmDataAnalyzeManager* dataManager = nmDataAnalyzeManager::getCurrentInstance();
     nmDataReservoir reservoirData = dataManager->getReservoirDataCopy();
 
-    // paramIndex 是粒子 position 中的索引；i 是完整 9 个参数体系中的索引。
+    // paramIndex 是粒子 position 中的索引；i 是完整 10 个参数体系中的索引。
     // 只有 m_parameterSelected[i] 为 true 时，才从 parameters 中消费一个值。
     int paramIndex = 0;
 
@@ -5754,7 +5789,7 @@ void nmCalculationAutoFitPSO::updateWellParameters(const QVector<double>& parame
     // 更新目标井上的拟合参数。目前井级可拟合参数主要是：
     // - skin：写入第一个 perforation；
     // - wellboreC：写入井筒储集系数；
-    // - Dfc：只写入垂直压裂井或多段压裂水平井的裂缝导流能力。
+    // - Dfc/裂缝半长：只写入垂直压裂井或多段压裂水平井。
     // 如果目标井不存在或没有射孔数据，这里只记录 debug，不抛异常。
     nmDataAnalyzeManager* dataManager = nmDataAnalyzeManager::getCurrentInstance();
 
@@ -5811,6 +5846,24 @@ void nmCalculationAutoFitPSO::updateWellParameters(const QVector<double>& parame
                         }
                     }
                     break;
+
+                case 9: { // 裂缝半长
+                        // 直接修改井对象中的属性，复用已有信号重算裂缝端点。
+                        if(pWell->getWellType() == NM_WELL_MODEL::Vertical_Fractured_Well) {
+                            nmDataVerticalFracturedWell* fracturedWell =
+                                dynamic_cast<nmDataVerticalFracturedWell*>(pWell);
+                            if(fracturedWell) {
+                                fracturedWell->getFractureHalfLength().setValue(value);
+                            }
+                        } else if(pWell->getWellType() == NM_WELL_MODEL::Horizontal_Fractured_Well) {
+                            nmDataHorizontalFracturedWell* fracturedWell =
+                                dynamic_cast<nmDataHorizontalFracturedWell*>(pWell);
+                            if(fracturedWell) {
+                                fracturedWell->getFractureHalfLength().setValue(value);
+                            }
+                        }
+                    }
+                    break;
             }
 
             paramIndex++;

Src/nmNum/nmData/nmDataAutomaticFitting.cpp

@@ -12,6 +12,7 @@ nmDataAutomaticFitting::nmDataAutomaticFitting()
 	m_cfSelected = false;                  // 默认不选中
     m_swiSelected = false;				  // 默认不选中
 	m_fractureConductivitySelected = false; // 仅压裂井可用，默认不选中
+	m_fractureHalfLengthSelected = false;  // 仅压裂井可用，默认不选中
 
 	// 拟合上下界不再使用固定默认值，由自动拟合窗口按数据对象初值和物理边界生成。
 	m_permeabilityMax = nmDataAttribute("Permeability Max", QVariant(), "Darcy");
@@ -23,6 +24,7 @@ nmDataAutomaticFitting::nmDataAutomaticFitting()
 	m_cfMax = nmDataAttribute("Cf Max", QVariant(), "");
 	m_swiMax = nmDataAttribute("Swi Max", QVariant(), "");
 	m_fractureConductivityMax = nmDataAttribute("Fracture Conductivity Max", QVariant(), "md.m");
+	m_fractureHalfLengthMax = nmDataAttribute("Fracture Half Length Max", QVariant(), "m");
 
 	m_permeabilityMin = nmDataAttribute("Permeability Min", QVariant(), "Darcy");
 	m_skinMin = nmDataAttribute("Skin Min", QVariant(), "");
@@ -33,6 +35,7 @@ nmDataAutomaticFitting::nmDataAutomaticFitting()
 	m_cfMin = nmDataAttribute("Cf Min", QVariant(), "");
 	m_swiMin = nmDataAttribute("Swi Min", QVariant(), "");
 	m_fractureConductivityMin = nmDataAttribute("Fracture Conductivity Min", QVariant(), "md.m");
+	m_fractureHalfLengthMin = nmDataAttribute("Fracture Half Length Min", QVariant(), "m");
 
 	// 初始化迭代参数
 	m_iterationCount = nmDataAttribute("Iteration Count", 20, "");
@@ -63,6 +66,7 @@ nmDataAutomaticFitting& nmDataAutomaticFitting::operator=(const nmDataAutomaticF
 		m_cfSelected = other.m_cfSelected;
 		m_swiSelected = other.m_swiSelected;
 		m_fractureConductivitySelected = other.m_fractureConductivitySelected;
+		m_fractureHalfLengthSelected = other.m_fractureHalfLengthSelected;
 
 		// 复制参数最大值
 		m_permeabilityMax = other.m_permeabilityMax;
@@ -74,6 +78,7 @@ nmDataAutomaticFitting& nmDataAutomaticFitting::operator=(const nmDataAutomaticF
 		m_cfMax = other.m_cfMax;
 		m_swiMax = other.m_swiMax;
 		m_fractureConductivityMax = other.m_fractureConductivityMax;
+		m_fractureHalfLengthMax = other.m_fractureHalfLengthMax;
 
 		// 复制参数最小值
 		m_permeabilityMin = other.m_permeabilityMin;
@@ -85,6 +90,7 @@ nmDataAutomaticFitting& nmDataAutomaticFitting::operator=(const nmDataAutomaticF
 		m_cfMin = other.m_cfMin;
 		m_swiMin = other.m_swiMin;
 		m_fractureConductivityMin = other.m_fractureConductivityMin;
+		m_fractureHalfLengthMin = other.m_fractureHalfLengthMin;
 
 		// 复制迭代参数
 		m_iterationCount = other.m_iterationCount;
@@ -109,6 +115,7 @@ rapidjson::Value nmDataAutomaticFitting::ToJsonValue(rapidjson::Document::Alloca
 	fittingObject.AddMember("CfSelected", m_cfSelected, allocator);
     fittingObject.AddMember("SwiSelected", m_swiSelected, allocator);
 	fittingObject.AddMember("FractureConductivitySelected", m_fractureConductivitySelected, allocator);
+	fittingObject.AddMember("FractureHalfLengthSelected", m_fractureHalfLengthSelected, allocator);
 
 	// 序列化参数最大值
 	fittingObject.AddMember("PermeabilityMax", m_permeabilityMax.ToJsonValue(allocator), allocator);
@@ -120,6 +127,7 @@ rapidjson::Value nmDataAutomaticFitting::ToJsonValue(rapidjson::Document::Alloca
 	fittingObject.AddMember("CfMax", m_cfMax.ToJsonValue(allocator), allocator);
 	fittingObject.AddMember("SwiMax", m_swiMax.ToJsonValue(allocator), allocator);
 	fittingObject.AddMember("FractureConductivityMax", m_fractureConductivityMax.ToJsonValue(allocator), allocator);
+	fittingObject.AddMember("FractureHalfLengthMax", m_fractureHalfLengthMax.ToJsonValue(allocator), allocator);
 
 	// 序列化参数最小值
 	fittingObject.AddMember("PermeabilityMin", m_permeabilityMin.ToJsonValue(allocator), allocator);
@@ -131,6 +139,7 @@ rapidjson::Value nmDataAutomaticFitting::ToJsonValue(rapidjson::Document::Alloca
 	fittingObject.AddMember("CfMin", m_cfMin.ToJsonValue(allocator), allocator);
 	fittingObject.AddMember("SwiMin", m_swiMin.ToJsonValue(allocator), allocator);
 	fittingObject.AddMember("FractureConductivityMin", m_fractureConductivityMin.ToJsonValue(allocator), allocator);
+	fittingObject.AddMember("FractureHalfLengthMin", m_fractureHalfLengthMin.ToJsonValue(allocator), allocator);
 
 	// 序列化迭代参数
 	fittingObject.AddMember("IterationCount", m_iterationCount.ToJsonValue(allocator), allocator);
@@ -171,6 +180,9 @@ void nmDataAutomaticFitting::FromJsonValue(const rapidjson::Value& jsonValue)
 	if (jsonValue.HasMember("FractureConductivitySelected") && jsonValue["FractureConductivitySelected"].IsBool()) {
 		m_fractureConductivitySelected = jsonValue["FractureConductivitySelected"].GetBool();
 	}
+	if (jsonValue.HasMember("FractureHalfLengthSelected") && jsonValue["FractureHalfLengthSelected"].IsBool()) {
+		m_fractureHalfLengthSelected = jsonValue["FractureHalfLengthSelected"].GetBool();
+	}
 
 	// 反序列化参数最大值
 	if (jsonValue.HasMember("PermeabilityMax") && jsonValue["PermeabilityMax"].IsObject()) {
@@ -200,6 +212,9 @@ void nmDataAutomaticFitting::FromJsonValue(const rapidjson::Value& jsonValue)
 	if (jsonValue.HasMember("FractureConductivityMax") && jsonValue["FractureConductivityMax"].IsObject()) {
 		m_fractureConductivityMax.FromJsonValue(jsonValue["FractureConductivityMax"]);
 	}
+	if (jsonValue.HasMember("FractureHalfLengthMax") && jsonValue["FractureHalfLengthMax"].IsObject()) {
+		m_fractureHalfLengthMax.FromJsonValue(jsonValue["FractureHalfLengthMax"]);
+	}
 
 	// 反序列化参数最小值
 	if (jsonValue.HasMember("PermeabilityMin") && jsonValue["PermeabilityMin"].IsObject()) {
@@ -229,6 +244,9 @@ void nmDataAutomaticFitting::FromJsonValue(const rapidjson::Value& jsonValue)
 	if (jsonValue.HasMember("FractureConductivityMin") && jsonValue["FractureConductivityMin"].IsObject()) {
 		m_fractureConductivityMin.FromJsonValue(jsonValue["FractureConductivityMin"]);
 	}
+	if (jsonValue.HasMember("FractureHalfLengthMin") && jsonValue["FractureHalfLengthMin"].IsObject()) {
+		m_fractureHalfLengthMin.FromJsonValue(jsonValue["FractureHalfLengthMin"]);
+	}
 
 	// 反序列化迭代参数
 	if (jsonValue.HasMember("IterationCount") && jsonValue["IterationCount"].IsObject()) {
@@ -274,6 +292,9 @@ void nmDataAutomaticFitting::setSwiSelected(bool selected) { m_swiSelected = sel
 bool nmDataAutomaticFitting::getFractureConductivitySelected() const { return m_fractureConductivitySelected; }
 void nmDataAutomaticFitting::setFractureConductivitySelected(bool selected) { m_fractureConductivitySelected = selected; }
 
+bool nmDataAutomaticFitting::getFractureHalfLengthSelected() const { return m_fractureHalfLengthSelected; }
+void nmDataAutomaticFitting::setFractureHalfLengthSelected(bool selected) { m_fractureHalfLengthSelected = selected; }
+
 
 // Getter and Setter implementations for Max values
 nmDataAttribute& nmDataAutomaticFitting::getPermeabilityMax() { return m_permeabilityMax; }
@@ -303,6 +324,9 @@ void nmDataAutomaticFitting::setSwiMax(const nmDataAttribute& swiMax) { m_swiMax
 nmDataAttribute& nmDataAutomaticFitting::getFractureConductivityMax() { return m_fractureConductivityMax; }
 void nmDataAutomaticFitting::setFractureConductivityMax(const nmDataAttribute& fractureConductivityMax) { m_fractureConductivityMax = fractureConductivityMax; }
 
+nmDataAttribute& nmDataAutomaticFitting::getFractureHalfLengthMax() { return m_fractureHalfLengthMax; }
+void nmDataAutomaticFitting::setFractureHalfLengthMax(const nmDataAttribute& fractureHalfLengthMax) { m_fractureHalfLengthMax = fractureHalfLengthMax; }
+
 // Getter and Setter implementations for Min values
 nmDataAttribute& nmDataAutomaticFitting::getPermeabilityMin() { return m_permeabilityMin; }
 void nmDataAutomaticFitting::setPermeabilityMin(const nmDataAttribute& permeabilityMin) { m_permeabilityMin = permeabilityMin; }
@@ -331,6 +355,9 @@ void nmDataAutomaticFitting::setSwiMin(const nmDataAttribute& swiMin) { m_swiMin
 nmDataAttribute& nmDataAutomaticFitting::getFractureConductivityMin() { return m_fractureConductivityMin; }
 void nmDataAutomaticFitting::setFractureConductivityMin(const nmDataAttribute& fractureConductivityMin) { m_fractureConductivityMin = fractureConductivityMin; }
 
+nmDataAttribute& nmDataAutomaticFitting::getFractureHalfLengthMin() { return m_fractureHalfLengthMin; }
+void nmDataAutomaticFitting::setFractureHalfLengthMin(const nmDataAttribute& fractureHalfLengthMin) { m_fractureHalfLengthMin = fractureHalfLengthMin; }
+
 // Getter and Setter implementations for iteration parameters
 nmDataAttribute& nmDataAutomaticFitting::getIterationCount() { return m_iterationCount; }
 void nmDataAutomaticFitting::setIterationCount(const nmDataAttribute& iterationCount) { m_iterationCount = iterationCount; }

Src/nmNum/nmSubWxs/nmWxAutomaticFitting.cpp

@@ -43,13 +43,10 @@ bool nmAutoFitUiNearlyEqual(double left, double right)
 	return std::fabs(left - right) <= DBL_EPSILON * 8.0 * scale;
 }
 
-// 从系统参数表读取物理边界，读取失败时保留调用方提供的兜底边界。
+// 从系统参数表读取物理边界，所有自动拟合参数统一以参数表配置为准。
 bool nmAutoFitReadPhysicalRange(const char* parameterName,
-	double fallbackMin, double fallbackMax, double& minValue, double& maxValue)
+	double& minValue, double& maxValue)
 {
-	minValue = fallbackMin;
-	maxValue = fallbackMax;
-
 	iSysParaHelper* paraHelper = _paraHelper;
 	if(!paraHelper) {
 		return false;
@@ -245,6 +242,7 @@ void nmWxAutomaticFitting::updateParameterVisibility(QTableWidget* table, NM_SOL
 		}
 	}
 	setParameterRowVisible(table, 8, showFractureConductivity); // Dfc
+	setParameterRowVisible(table, 9, showFractureConductivity); // 裂缝半长
 
 	renumberVisibleParameterRows(table);
 }
@@ -255,29 +253,15 @@ bool nmWxAutomaticFitting::getPhysicalParameterRange(int parameterIndex,
 {
 	static const char* parameterNames[] = {
 		"Result_K", "Result_W_Skin", "Result_W_C", "Result_phi",
-		"Result_h", "Result_Cti", "Result_Cf", "Result_Swi", "Result_W_Dfc"
-	};
-	// KAPPA 的边界使用 md、ft、bbl/psi；自动拟合界面使用 Darcy、m、m^3/MPa，
-	// 这里统一换算到界面和 PSO 实际使用的单位：K 除以 1000，h 由 ft 换成 m，
-	// 井筒储集系数的 4.33667154546306e34 bbl/psi 对应约 1e36 m^3/MPa。
-	// Ct/Cf/Swi/Dfc 沿用模型参数表边界。
-	static const double physicalMin[] = {
-		1.01325027383089e-18, -5.0, 0.0, 1.0e-4, 1.0e-5, 1.0e-30, 1.0e-30, 0.0, 0.0
-	};
-	static const double physicalMax[] = {
-		1.01325027383089e42, 5000.0, 1.0e36, 0.9999, 1.0e9, 10.0, 10.0, 1.0, 1.0e30
+		"Result_h", "Result_Cti", "Result_Cf", "Result_Swi", "Result_W_Dfc",
+		"W_FractureHalfLength"
 	};
-
-	if(parameterIndex < 0 || parameterIndex >= 9) {
+	if(parameterIndex < 0 || parameterIndex >= 10) {
 		return false;
 	}
 
-	minValue = physicalMin[parameterIndex];
-	maxValue = physicalMax[parameterIndex];
-	bool rangeRead = true;
-	if(parameterIndex >= 5) {
-		rangeRead = nmAutoFitReadPhysicalRange(parameterNames[parameterIndex],
-			physicalMin[parameterIndex], physicalMax[parameterIndex], minValue, maxValue);
+	if(!nmAutoFitReadPhysicalRange(parameterNames[parameterIndex], minValue, maxValue)) {
+		return false;
 	}
 
 	if(parameterIndex == 7) {
@@ -294,7 +278,7 @@ bool nmWxAutomaticFitting::getPhysicalParameterRange(int parameterIndex,
 		}
 	}
 
-	return rangeRead;
+	return true;
 }
 
 // 将一组上下界同步到表格和自动拟合数据，保证 PSO 读取到同一份配置。
@@ -357,6 +341,10 @@ void nmWxAutomaticFitting::setParameterRange(int parameterIndex,
 		automaticFittingData.getFractureConductivityMin().setValue(minValue);
 		automaticFittingData.getFractureConductivityMax().setValue(maxValue);
 		break;
+	case 9:
+		automaticFittingData.getFractureHalfLengthMin().setValue(minValue);
+		automaticFittingData.getFractureHalfLengthMax().setValue(maxValue);
+		break;
 	default:
 		break;
 	}
@@ -369,14 +357,16 @@ void nmWxAutomaticFitting::setParameterRange(int parameterIndex,
 void nmWxAutomaticFitting::updateRangeForParameter(int parameterIndex,
 	double centerValue)
 {
-	if(!m_parameterTable || parameterIndex < 0 || parameterIndex >= 9
+	if(!m_parameterTable || parameterIndex < 0 || parameterIndex >= 10
 	   || !nmAutoFitUiIsFinite(centerValue)) {
 		return;
 	}
 
 	double physicalMin = 0.0;
 	double physicalMax = 0.0;
-	getPhysicalParameterRange(parameterIndex, physicalMin, physicalMax);
+	if(!getPhysicalParameterRange(parameterIndex, physicalMin, physicalMax)) {
+		return;
+	}
 
 	double reference = centerValue;
 	const bool positiveParameter = parameterIndex != 1;
@@ -437,7 +427,7 @@ void nmWxAutomaticFitting::initializeSuggestedParameterRanges()
 		return;
 	}
 
-	for(int parameterIndex = 0; parameterIndex < 9; ++parameterIndex) {
+	for(int parameterIndex = 0; parameterIndex < 10; ++parameterIndex) {
 		QTableWidgetItem* initialItem = m_parameterTable->item(parameterIndex, 3);
 		if(initialItem) {
 			bool initialOk = false;
@@ -448,8 +438,8 @@ void nmWxAutomaticFitting::initializeSuggestedParameterRanges()
 				// 数据对象没有提供该初值时使用完整物理区间，不回退到旧的默认范围。
 				double physicalMin = 0.0;
 				double physicalMax = 0.0;
-				getPhysicalParameterRange(parameterIndex, physicalMin, physicalMax);
-				if(physicalMax >= physicalMin) {
+				if(getPhysicalParameterRange(parameterIndex, physicalMin, physicalMax)
+				   && physicalMax >= physicalMin) {
 					setParameterRange(parameterIndex, physicalMin, physicalMax);
 				}
 			}
@@ -464,7 +454,7 @@ void nmWxAutomaticFitting::normalizeSavedParameterRanges()
 		return;
 	}
 
-	for(int parameterIndex = 0; parameterIndex < 9; ++parameterIndex) {
+	for(int parameterIndex = 0; parameterIndex < 10; ++parameterIndex) {
 		QTableWidgetItem* minItem = m_parameterTable->item(parameterIndex, 2);
 		QTableWidgetItem* maxItem = m_parameterTable->item(parameterIndex, 4);
 		QTableWidgetItem* initialItem = m_parameterTable->item(parameterIndex, 3);
@@ -474,7 +464,9 @@ void nmWxAutomaticFitting::normalizeSavedParameterRanges()
 
 		double physicalMin = 0.0;
 		double physicalMax = 0.0;
-		getPhysicalParameterRange(parameterIndex, physicalMin, physicalMax);
+		if(!getPhysicalParameterRange(parameterIndex, physicalMin, physicalMax)) {
+			continue;
+		}
 		bool savedMinOk = false;
 		bool savedMaxOk = false;
 		const double savedMin = minItem->text().toDouble(&savedMinOk);
@@ -513,18 +505,19 @@ bool nmWxAutomaticFitting::validateParameterTable(QString& errorMessage, int par
 		errorMessage = tr("The parameter table is unavailable.");
 		return false;
 	}
-	if(parameterIndex < -1 || parameterIndex >= 9) {
+	if(parameterIndex < -1 || parameterIndex >= 10) {
 		errorMessage = tr("The parameter row is invalid.");
 		return false;
 	}
 
 	static const char* parameterNames[] = {
 		"Permeability", "Skin", "Wellbore storage", "Porosity",
-		"Thickness", "Ct", "Cf", "Swi", "Fracture conductivity"
+		"Thickness", "Ct", "Cf", "Swi", "Fracture conductivity",
+		"Fracture half length"
 	};
 
 	const int firstParameterIndex = parameterIndex < 0 ? 0 : parameterIndex;
-	const int lastParameterIndex = parameterIndex < 0 ? 9 : parameterIndex + 1;
+	const int lastParameterIndex = parameterIndex < 0 ? 10 : parameterIndex + 1;
 	for(int currentParameterIndex = firstParameterIndex;
 		currentParameterIndex < lastParameterIndex; ++currentParameterIndex) {
 		// 隐藏参数不参与当前模型拟合，不用它们的历史值阻塞当前设置。
@@ -558,8 +551,8 @@ bool nmWxAutomaticFitting::validateParameterTable(QString& errorMessage, int par
 
 		double physicalMin = 0.0;
 		double physicalMax = 0.0;
-		getPhysicalParameterRange(currentParameterIndex, physicalMin, physicalMax);
-		if(!nmAutoFitUiIsFinite(physicalMin) || !nmAutoFitUiIsFinite(physicalMax)
+		if(!getPhysicalParameterRange(currentParameterIndex, physicalMin, physicalMax)
+		   || !nmAutoFitUiIsFinite(physicalMin) || !nmAutoFitUiIsFinite(physicalMax)
 		   || physicalMax < physicalMin) {
 			errorMessage = tr("The physical range of %1 is invalid.")
 				.arg(tr(parameterNames[currentParameterIndex]));
@@ -717,7 +710,7 @@ void nmWxAutomaticFitting::setupUI()
 void nmWxAutomaticFitting::setupParameterTable()
 {
 	// 创建表格
-	m_parameterTable = new QTableWidget(9, 6, this);
+	m_parameterTable = new QTableWidget(10, 6, this);
 
 	// 设置表头
 	QStringList headers;
@@ -837,6 +830,16 @@ void nmWxAutomaticFitting::setupParameterTable()
 	m_parameterTable->setItem(8, 4, new QTableWidgetItem(QString::number(automaticFittingData.getFractureConductivityMax().getValue().toDouble())));
 	m_parameterTable->setItem(8, 5, new QTableWidgetItem(tr("md.m")));
 
+	// 裂缝半长。该行与 Dfc 一样只对压裂井显示，初值从当前目标井读取。
+	m_parameterTable->setItem(9, 0, new QTableWidgetItem("10"));
+	m_fractureHalfLengthCheckBox = new QCheckBox(tr("Fracture half length"));
+	m_fractureHalfLengthCheckBox->setChecked(automaticFittingData.getFractureHalfLengthSelected());
+	m_parameterTable->setCellWidget(9, 1, m_fractureHalfLengthCheckBox);
+	m_parameterTable->setItem(9, 2, new QTableWidgetItem(QString::number(automaticFittingData.getFractureHalfLengthMin().getValue().toDouble())));
+	m_parameterTable->setItem(9, 3, new QTableWidgetItem());
+	m_parameterTable->setItem(9, 4, new QTableWidgetItem(QString::number(automaticFittingData.getFractureHalfLengthMax().getValue().toDouble())));
+	m_parameterTable->setItem(9, 5, new QTableWidgetItem(tr("m")));
+
 	// 设置表格行为
 	for(int i = 0; i < m_parameterTable->rowCount(); ++i) {
 		for(int j = 0; j < 6; ++j) {
@@ -1037,6 +1040,7 @@ void nmWxAutomaticFitting::onReverseSelection()
 	if(!m_parameterTable->isRowHidden(6)) m_cfCheckBox->setChecked(!m_cfCheckBox->isChecked());
 	if(!m_parameterTable->isRowHidden(7)) m_swiCheckBox->setChecked(!m_swiCheckBox->isChecked());
 	if(!m_parameterTable->isRowHidden(8)) m_dfcCheckBox->setChecked(!m_dfcCheckBox->isChecked());
+	if(!m_parameterTable->isRowHidden(9)) m_fractureHalfLengthCheckBox->setChecked(!m_fractureHalfLengthCheckBox->isChecked());
 }
 
 void nmWxAutomaticFitting::onParameterTableItemChanged(QTableWidgetItem* item)
@@ -1155,7 +1159,8 @@ void nmWxAutomaticFitting::onAccept()
 		m_cCheckBox->isChecked() || m_phiCheckBox->isChecked() ||
 		m_hCheckBox->isChecked() ||
 		m_ctCheckBox->isChecked() || m_cfCheckBox->isChecked() ||
-		m_swiCheckBox->isChecked() || m_dfcCheckBox->isChecked();
+		m_swiCheckBox->isChecked() || m_dfcCheckBox->isChecked() ||
+		m_fractureHalfLengthCheckBox->isChecked();
 
 	if(!hasSelectedParams) {
 		QMessageBox::warning(this, tr("Warning"), tr("Please select at least one parameter for optimization!"));
@@ -1174,6 +1179,7 @@ void nmWxAutomaticFitting::onAccept()
 	if(m_cfCheckBox->isChecked()) selectedParameterNames << tr("Cf");
 	if(m_swiCheckBox->isChecked()) selectedParameterNames << tr("Swi");
 	if(m_dfcCheckBox->isChecked()) selectedParameterNames << tr("Fracture conductivity");
+	if(m_fractureHalfLengthCheckBox->isChecked()) selectedParameterNames << tr("Fracture half length");
 
 	// 启动自动拟合 - 传递双对数历史数据
 	startAutoFitting(targetLogLogData, selectedParameterNames, selectedWellName);
@@ -1194,6 +1200,7 @@ void nmWxAutomaticFitting::onWellSelected(int index)
 	double skinValue = 0.0;
 	double wellboreStorageValue = 0.0;
 	double fractureConductivityValue = 0.0;
+	double fractureHalfLengthValue = 0.0;
 
 	// 查找垂直井
 	for(int i = 0; i < m_verticalWells.size(); ++i) {
@@ -1224,6 +1231,7 @@ void nmWxAutomaticFitting::onWellSelected(int index)
 				skinValue = m_verticalFracturedWells[i].getPerforation(0)->getSkin().getValue().toDouble();
 				wellboreStorageValue = m_verticalFracturedWells[i].getWellboreStorage().getValue().toDouble();
 				fractureConductivityValue = m_verticalFracturedWells[i].getDfc().getValue().toDouble();
+				fractureHalfLengthValue = m_verticalFracturedWells[i].getFractureHalfLength().getValue().toDouble();
 				fracturedWell = true;
 				found = true;
 				break;
@@ -1238,6 +1246,7 @@ void nmWxAutomaticFitting::onWellSelected(int index)
 				skinValue = m_horizontalFracturedWells[i].getPerforation(0)->getSkin().getValue().toDouble();
 				wellboreStorageValue = m_horizontalFracturedWells[i].getWellboreStorage().getValue().toDouble();
 				fractureConductivityValue = m_horizontalFracturedWells[i].getDfc().getValue().toDouble();
+				fractureHalfLengthValue = m_horizontalFracturedWells[i].getFractureHalfLength().getValue().toDouble();
 				fracturedWell = true;
 				found = true;
 				break;
@@ -1268,12 +1277,16 @@ void nmWxAutomaticFitting::onWellSelected(int index)
 		if(fracturedWell && m_parameterTable->item(8, 3)) {
 			m_parameterTable->item(8, 3)->setText(QString::number(fractureConductivityValue));
 		}
+		if(fracturedWell && m_parameterTable->item(9, 3)) {
+			m_parameterTable->item(9, 3)->setText(QString::number(fractureHalfLengthValue));
+		}
 
 		if(m_autoParameterRanges) {
 			updateRangeForParameter(1, skinValue);
 			updateRangeForParameter(2, wellboreStorageValue);
 			if(fracturedWell) {
 				updateRangeForParameter(8, fractureConductivityValue);
+				updateRangeForParameter(9, fractureHalfLengthValue);
 			}
 		}
 	}
@@ -1291,6 +1304,7 @@ void nmWxAutomaticFitting::setAutomaticFittingValue()
 	automaticFittingData.setCfSelected(m_cfCheckBox->isChecked());
 	automaticFittingData.setSwiSelected(m_swiCheckBox->isChecked());
 	automaticFittingData.setFractureConductivitySelected(m_dfcCheckBox->isChecked());
+	automaticFittingData.setFractureHalfLengthSelected(m_fractureHalfLengthCheckBox->isChecked());
 	automaticFittingData.setSurrogateScreeningEnabled(m_surrogateCombo && m_surrogateCombo->currentIndex() == 1);
 
 	// 保存渗透率的最小值和最大值
@@ -1329,6 +1343,10 @@ void nmWxAutomaticFitting::setAutomaticFittingValue()
 	automaticFittingData.getFractureConductivityMin().setValue(m_parameterTable->item(8, 2)->text().toDouble());
 	automaticFittingData.getFractureConductivityMax().setValue(m_parameterTable->item(8, 4)->text().toDouble());
 
+	// 保存裂缝半长的最小值和最大值
+	automaticFittingData.getFractureHalfLengthMin().setValue(m_parameterTable->item(9, 2)->text().toDouble());
+	automaticFittingData.getFractureHalfLengthMax().setValue(m_parameterTable->item(9, 4)->text().toDouble());
+
 	// 保存迭代参数
 	automaticFittingData.getIterationCount().setValue(m_iterationEdit->text().toInt());
 	automaticFittingData.getErrorTolerance().setValue(m_errorLimitEdit->text().toDouble());
@@ -1353,6 +1371,7 @@ void nmWxAutomaticFitting::setAutomaticFittingValue()
 		double newSkinValue = m_parameterTable->item(1, 3)->text().toDouble();
 		double newWellboreStorageValue = m_parameterTable->item(2, 3)->text().toDouble();
 		double newFractureConductivityValue = m_parameterTable->item(8, 3)->text().toDouble();
+		double newFractureHalfLengthValue = m_parameterTable->item(9, 3)->text().toDouble();
 
 		// 直接从数据管理器获取目标井
 		nmDataAnalyzeManager* manager = nmDataAnalyzeManager::getCurrentInstance();
@@ -1382,6 +1401,7 @@ void nmWxAutomaticFitting::setAutomaticFittingValue()
 					nmDataAttribute dfc = fracturedWell->getDfc();
 					dfc.setValue(newFractureConductivityValue);
 					fracturedWell->setDfc(dfc);
+					fracturedWell->getFractureHalfLength().setValue(newFractureHalfLengthValue);
 				}
 			} else if(wellType == NM_WELL_MODEL::Horizontal_Fractured_Well) {
 				nmDataHorizontalFracturedWell* fracturedWell = dynamic_cast<nmDataHorizontalFracturedWell*>(pTargetWell);
@@ -1389,6 +1409,7 @@ void nmWxAutomaticFitting::setAutomaticFittingValue()
 					nmDataAttribute dfc = fracturedWell->getDfc();
 					dfc.setValue(newFractureConductivityValue);
 					fracturedWell->setDfc(dfc);
+					fracturedWell->getFractureHalfLength().setValue(newFractureHalfLengthValue);
 				}
 			}
 
@@ -1622,6 +1643,7 @@ void nmWxAutomaticFitting::updateBestParametersToTable()
 	if(m_cfCheckBox->isChecked()) enabledParams.append(6);     // 岩石压缩系数
 	if(m_swiCheckBox->isChecked()) enabledParams.append(7);    // 初始含水饱和度
 	if(m_dfcCheckBox->isChecked()) enabledParams.append(8);    // 裂缝导流能力
+	if(m_fractureHalfLengthCheckBox->isChecked()) enabledParams.append(9); // 裂缝半长
 
 	// 更新参数值和范围
 	for (int i = 0; i < bestSolution.size() && i < enabledParams.size(); ++i) {