数学建模竞赛中工程优化问题的实战解析：以巷道支护为例
1. 从“煤矿巷道支护”看2025五一杯C题的实战价值
每年五一杯数学建模竞赛的C题，总是能精准地戳中一个交叉学科的“痒点”。2025年的这道题，直接把场景锚定在了“煤矿巷道支护”上。如果你乍一看觉得这离自己很遥远，那可能就错过了这道题最精妙的设计。这道题的价值，远不止于解决一个工程问题，它本质上是一次对“如何将复杂现实抽象为可计算模型”的完整演练。对于所有理工科，尤其是自动化、机械、土木、采矿乃至数据科学方向的同学来说，这是一次绝佳的练兵机会。它考察的不是你对煤矿有多了解，而是你面对一个充满不确定性、约束条件复杂、目标函数模糊的真实世界问题时，如何运用数学工具进行建模、求解和优化的核心能力。
我参加过也指导过不少数学建模比赛，一个深刻的体会是：最能拉开差距的，往往不是最后那个炫酷的算法，而是最初的问题分析和模型建立阶段。五一杯C题这次选的“巷道支护”，就是一个典型的“麻雀虽小，五脏俱全”的案例。它涉及力学分析（围岩压力）、材料科学（支护材料属性）、优化理论（成本与安全平衡），甚至还有一点随机性（地质条件的不确定性）。题目没有给你一个现成的公式往里套，它给你的是一个需要你自己去定义、去量化的系统。这正是数学建模从“做题”转向“解决问题”的关键一步。
2. 题目核心：在安全与成本的钢丝上寻找最优解
我们先来拆解一下这道题的核心矛盾，这也是所有工程优化问题的母题：如何在有限的资源（成本）下，最大化系统的可靠性（安全）？具体到巷道支护，这个矛盾被具象化为几个关键变量：
决策变量：这通常就是我们要设计的支护方案参数。可能包括：
- 支护材料的类型（如锚杆、钢拱架、混凝土衬砌的组合）。
- 关键尺寸参数（如锚杆的长度、直径、间距；拱架的截面尺寸、间距）。
- 支护的布置方式（如全断面支护还是局部加强支护）。
目标函数：我们需要优化什么？题目很可能会给出一个多目标甚至需要你自行构造的目标。常见的有：
- 经济性目标：总成本最小化。成本包括材料成本、安装人工成本、后期维护成本。
- 安全性目标：巷道稳定性最大化，或变形量最小化。这通常需要通过一个力学模型来量化，例如计算支护后围岩的应力状态、位移场，判断是否超过安全阈值。
- 很多时候，这两个目标是冲突的。更安全通常意味着更昂贵。因此，问题很可能被构建为一个多目标优化问题，你需要找到一系列“帕累托最优解”，即在不损害一个目标的情况下，无法再改进另一个目标的解集。
约束条件：这是模型的边界，也是现实世界的物理和工程限制。
- 力学约束：支护结构本身的强度必须大于其所受的荷载（围岩压力）。围岩的变形必须在允许范围内（如不影响设备通行、不引发塌方）。
- 几何约束：支护结构的尺寸受巷道断面形状和尺寸的限制。
- 施工约束：某些支护方案可能因为施工工艺太复杂或工期太长而被排除。
- 规范约束：必须符合国家或行业的矿山安全设计规范。
为什么说这个设定很巧妙？因为它迫使你思考模型的“粒度”。你是把围岩当作均匀的弹性体，还是考虑其节理、裂隙？你是把支护结构当作刚性体，还是考虑其弹塑性变形？你对成本的建模是只算材料费，还是包含了折现后的全生命周期成本？不同的假设，会导向完全不同的模型复杂度和求解路径。这里没有标准答案，只有“在给定时间和计算资源下，最合理的近似”。这正是数学建模的艺术所在。
3. 模型构建的三种典型路径与选择逻辑
面对这样的问题，队伍通常会走向几种不同的建模路径。选择哪一条，直接决定了论文的深度和获奖潜力。
3.1 路径一：基于力学原理的解析或数值模型
这是最“硬核”的路径，适合队伍中有力学或土木背景成员的情况。
- 核心思路：建立巷道围岩-支护结构相互作用的力学模型。常见的方法有：
- 弹性力学解析解：对于圆形巷道，在均匀地应力场下，可以套用经典的“厚壁圆筒”或“芬纳公式”等解析解，计算围岩的应力重分布和塑性区范围。支护结构被视为对围岩提供均匀的径向支护力。
- 数值模拟方法：当巷道形状不规则、地质条件复杂时，解析解无能为力。这时需要引入有限元法或有限差分法，使用软件如ANSYS、ABAQUS、FLAC3D或COMSOL进行建模。你可以详细定义围岩的材料本构模型（弹性、弹塑性、应变软化等）、支护结构的单元类型，并进行仿真计算。
- 优势：物理意义清晰，结果可信度高，容易体现理论深度。如果能将数值模拟结果与简化解析解进行对比验证，会是非常大的亮点。
- 挑战：对专业知识要求高，软件学习成本大，计算可能耗时。在论文中需要花大量篇幅解释模型假设和参数取值的依据。
- 实操心得：如果走这条路，千万不要只做一个“黑箱”模拟。论文的重点不是展示软件操作，而是解释“为什么这样设置模型”。例如，你选择Drucker-Prager准则而不是Mohr-Coulomb准则来模拟岩体，理由是什么？支护结构与围岩的接触面是如何处理的？这些细节的讨论，才是评委看重的地方。
3.2 路径二：基于响应面与代理模型的优化框架
这是目前数学建模竞赛中非常主流且高效的思路，特别适合优化问题。
- 核心思路：承认完整的物理模型（如复杂的有限元模型）计算一次成本太高，不适合直接嵌入优化循环。因此，我们采用“两步走”策略：
- 实验设计：在决策变量的设计空间内，选取一批有代表性的样本点（如使用拉丁超立方抽样）。
- 构建代理模型：对每一个样本点，调用一次“高保真模型”（可以是简化力学模型，甚至是一个已知的函数）计算其响应（如最大位移、安全系数、总成本）。然后，用这些样本数据，训练一个计算快速的“代理模型”，比如多项式响应面、克里金模型或神经网络。
- 优化求解：将代理模型作为目标函数和约束，调用成熟的优化算法（如遗传算法、粒子群算法、序列二次规划）进行快速寻优。
- 优势：将耗时的仿真与快速的优化解耦，大大提高了求解效率。思路清晰，通用性强，易于在论文中呈现完整的流程图。
- 挑战：代理模型的精度直接影响优化结果的可靠性。需要谨慎选择实验设计方法和代理模型类型，并进行充分的模型验证（如留出检验集计算误差）。
- 实操心得：样本点的数量和质量是关键。样本太少，代理模型不准；样本太多，失去了“代理”的意义。一个实用的技巧是采用“自适应采样”：先在全局粗采样构建初始代理模型，然后在优化过程中，在最优解附近区域进行加密采样，迭代更新代理模型，逐步逼近真实最优解。
3.3 路径三：基于系统工程的综合评价与决策模型
当问题中的不确定性非常突出，或者目标难以用单一函数量化时，可以转向这条路径。
- 核心思路：将支护方案选择看作一个多属性决策问题。除了成本和变形量，你可能还需要考虑“施工便捷性”、“对环境的影响”、“技术成熟度”等定性或半定量指标。
- 评价指标体系构建：建立一套层次化的评价指标。
- 指标量化与标准化：将不同量纲的指标（如成本“元”、位移“毫米”、施工难度“等级”）通过极差法等方法标准化。
- 权重确定：使用层次分析法、熵权法或组合赋权法，确定各指标的相对重要性。
- 方案排序：使用TOPSIS、灰色关联分析、VIKOR等方法对有限个预设方案进行排序选优。
- 优势：能处理定性指标，更贴近实际的工程决策流程，体现系统性思维。
- 挑战：容易陷入主观性（特别是AHP中判断矩阵的构造）。如果最终只是对几个预设方案排序，可能显得创新性不足。
- 实操心得：将定性与定量结合。例如，可以先用量化模型（如路径一或二）生成一批“帕累托最优”的方案集，然后再用多属性决策方法，引入专家经验等定性因素，从这个优质方案集中选出最终的“推荐方案”。这样既保证了方案的科学性，又体现了决策的综合性。
注意：在实际比赛中，高水平队伍往往会混合使用这些路径。例如，用简化力学模型或数值模拟作为“高保真模型”，构建代理模型进行优化，得到帕累托前沿，最后用决策分析方法选择一个折中方案。这种“组合拳”最能体现建模的综合能力。
4. 求解算法选型：没有最好，只有最合适
模型建好了，用什么算法来求解？这道题很可能引出以下两类算法需求：
4.1 多目标优化算法
当你的模型有“成本最低”和“安全最高”两个冲突目标时，你需要的是多目标优化算法。它的目的不是找到一个唯一解，而是找到那条“帕累托前沿”。
- 经典算法：NSGA-II（带精英策略的非支配排序遗传算法）是绝对的主力。它概念清晰，实现成熟，有大量开源代码（如Python的pymoo库、DEAP库）。其核心是通过“非支配排序”和“拥挤度计算”来维持解集的多样性和收敛性。
- 为什么是它？对于这种可能非凸、非线性的问题，基于种群的进化算法具有全局搜索能力，不易陷入局部最优，且能一次性生成一组解，非常适合探索帕累托前沿的形状。
- 参数调优陷阱：NSGA-II有交叉概率、变异概率、种群大小等参数。新手常犯的错误是直接使用默认值。你必须进行参数敏感性分析。例如，在其他条件不变时，分别增大和减小种群规模，观察帕累托前沿的变化。如果变化不大，说明你的参数设置是稳健的；如果差异显著，则需谨慎选择。在论文中展示这个分析过程，是严谨性的体现。
4.2 混合整数规划/非线性规划求解器
如果你的决策变量中包含离散选择（例如，是否采用某种支护类型，是0-1变量），那么问题可能是一个混合整数非线性规划问题。
- 求解器选择：对于中小规模问题，可以尝试LINGO、GAMS等商业软件，或SCIP、Couenne等开源求解器。对于复杂问题，常采用启发式算法或将其松弛后求解。
- 简化策略：在竞赛有限时间内，一个讨巧的做法是分步优化。先离散地选择几种典型的支护类型组合，对每一种组合，将其尺寸参数作为连续变量进行优化，最后比较各组合下的最优结果。这实际上是将一个MINLP问题分解为多个NLP问题和一个枚举比较问题。
算法部分的核心建议：不要追求算法的复杂度，而要追求应用的恰当性和分析的深度。用一个经典的NSGA-II把问题分析透彻，远比生搬硬套一个不理解的“最新智能算法”要好。论文中必须包含算法关键参数的设置理由、收敛性分析（如迭代过程中种群适应度的变化曲线）以及对结果的可视化（二维或三维的帕累托前沿图）。
5. 论文写作与结果呈现：如何让评委眼前一亮
数学建模竞赛，归根结底是“作文”竞赛。模型再精妙，算法再高级，如果无法清晰、有说服力地表达出来，一切归零。针对这道题，论文写作有几个需要特别注意的点：
5.1 摘要：用一页纸讲一个完整的故事
摘要必须独立成篇，即使不读正文也能了解你的全部工作。结构可以遵循“问题简述 -> 建模思路 -> 求解方法 -> 主要结果 -> 结论特色”的逻辑。
- 模板句式要避免：切忌“本文首先……然后……最后……”。要用叙述性语言串联。例如：“针对煤矿巷道支护方案优化问题，本文构建了一个以经济成本与巷道顶板最大位移为双目标的优化模型。其中，支护抗力与围岩变形关系通过简化弹塑性力学模型描述，并采用响应面法构建了其代理模型以提升优化效率。运用NSGA-II算法求解得到帕累托最优解集，分析发现成本与安全存在明显的权衡关系。进一步引入熵权-TOPSIS决策法，综合考虑施工因素，推荐了某特定地质条件下的最优支护参数。该方案在成本增加不超过15%的前提下，将安全系数提升了约30%。”
- 关键结果必须量化：把最重要的数字写进去，比如“最优方案成本为XX元/米，顶板位移预计为YY毫米”。
5.2 模型假设与符号说明：严谨性的基石
这是很多队伍敷衍了事的部分，恰恰是评委判断你思维严谨性的第一关。
- 假设要合理且必要：例如，“假设围岩为均质、各向同性的连续介质”、“忽略地下水的影响”、“支护结构在服役期内材料性能不发生退化”。每一条假设都要说明其合理性（基于题目背景或一般工程常识）以及它如何简化了问题。
- 符号说明要用表格：三线表，列包括：符号、含义、单位。单位尤其重要，它体现了你对物理量的理解。如果全文单位统一，例如全部使用国际单位制，会是一个加分项。
5.3 灵敏度分析与模型检验：证明你的模型“靠谱”
这是将论文从“完成解答”提升到“优秀论文”的关键环节。
- 灵敏度分析：改变模型中的关键输入参数（如岩体的弹性模量、内聚力，或成本单价），观察输出结果（如最优成本、最优位移）的变化程度。这回答了“如果我的参数估计有误差，结论会多大程度上被颠覆？”这个问题。可以用图表展示，例如绘制某个参数变化时，帕累托前沿的移动情况。
- 模型检验：如果你的模型有解析解或可比较的已知案例，一定要进行对比。例如，将你的简化力学模型在特定条件下退化为经典公式，看结果是否一致。或者，将你的优化结果与工程经验中的常规方案进行对比，讨论差异及其原因。这证明了你的模型的有效性。
5.4 可视化：一图胜千言
- 帕累托前沿图：多目标优化的标准输出。二维图用散点图，三点图可以展示两个目标与一个关键决策变量的关系。用不同的颜色或形状区分不同的方案簇。
- 决策变量影响图：例如，用箱线图或趋势图展示当锚杆间距变化时，成本和位移的分布如何变化。这能直观揭示设计规律。
- 模型示意图：手绘或用绘图软件（如Visio, draw.io）画一个清晰的巷道支护剖面图，标注你的决策变量，能让评委迅速理解你的物理模型。
6. 备赛与实战中的避坑指南
结合这道具体的题目和多年观赛经验，我总结几个最容易“踩坑”的地方：
坑一：盲目追求模型复杂度。看到“巷道支护”，有的队伍第一反应就是去搞一个复杂的有限元模型。结果三天时间，光建模型、调试参数、等计算就花了两天半，最后根本没时间做优化和深入分析。原则是：用尽可能简单的模型抓住问题的本质。一个考虑均匀地应力和圆形巷道的简化解析模型，如果能清晰地揭示成本与安全的权衡关系，并完成完整的优化-决策流程，其价值远胜于一个设置不当、结果无法解释的复杂仿真。
坑二：忽略单位与量纲。这是低级错误但屡见不鲜。力学公式中应力是Pa，你代入的数据是MPa；成本计算中，长度单位一会儿是米一会儿是厘米。这会导致结果数量级完全错误。从建模第一刻起，就统一使用国际单位制，并在所有计算式中保持量纲一致。
坑三：优化结果缺乏工程解释。算法跑出了一组“最优解”，就直接写在论文里。评委可能会问：为什么锚杆间距最优值是1.2米，而不是1米或1.5米？你需要从力学原理上解释：间距太小，成本剧增但安全增益有限（边际效益递减）；间距太大，支护不足以形成连续承载拱。将数学最优解与工程常识关联起来，是画龙点睛之笔。
坑四：论文写成软件说明书或代码附录。论文的核心是逻辑、思想和结果，不是操作步骤。不要在正文里大段贴代码或软件操作截图。代码可以放在附录，正文中只需说明你用了什么方法、什么工具、关键参数如何设置，以及最终得到的结果。
坑五：团队分工与时间管理失控。三个人一上来都扎进建模里，最后没人写论文。理想的分工是：一人主攻模型建立与算法实现（编程能力强），一人主攻模型求解与结果分析（数学和逻辑好），一人主攻论文写作与图表绘制（文字和视觉表达好）。但分工不能僵化，必须保持高频沟通。第一天下午必须确定基本模型框架，第二天下午必须开始论文写作，最后一天全天用于打磨论文和摘要。
2025年五一杯C题的“煤矿巷道支护”，是一道非常经典的工程优化建模题。它像一块试金石，能有效区分出哪些队伍只会套用模型，哪些队伍真正具备了分析问题、简化问题、解决问题的能力。它的价值不在于让你成为采矿专家，而在于训练你那种将模糊现实转化为清晰数学语言的思维框架。这种能力，无论是在学术研究还是未来的工程职业生涯中，都是无价的。