MCM 2021A 真菌网络问题（范文二）：Murray 锥度机制、规模缩放与有限源共享

一、摘要

在范文一建立"电阻树 + 基尔霍夫流量"模型的基础上，本篇深入剖析** Murray 锥度的物理与经济含义，并系统回答题目第二问：菌丝网络如何随蘑菇规模增长、又如何在有限食物源下竞争。我们首先证明，半径按 $r_e\propto q_e^{1/3}$ 锥度，能在"降低阻力"与"节约材料"之间取得最优平衡——同等送达下，锥度比全边等粗方案节省 $53.5\%$ 的生物量。随后考察规模缩放：随蘑菇数 $N$ 从 $4$ 增至 $28$ ，总菌丝长度以 $L\propto N^{0.53}$ 次线性增长（共享主干被摊销），而总生物量以 $B\propto N^{0.96}$ 近线性增长，使每株材料成本保持有界、呈轻微 U 形，即菌落生长不会带来人均成本爆炸。最后引入有限源通量**约束：当两菌共享同一食物源且源最大通量仅够一菌时，每菌实际获得量被腰斩（惩罚约 $48\%$ ），说明"源通量"才是竞争的真正瓶颈，而非建网本身。

二、问题重述与建模动机

题目第二问要求：在生长约束下优化网络，使养分送达最大化或成本最小化，并讨论规模与共享。范文一已说明 MST 拓扑最优，本篇进一步回答：① 为什么半径必须锥度而非等粗？② 菌落扩大时，建网成本如何随规模演化？③ 多个菌落共享同一食物源会发生什么？这三个子问题分别对应微观结构设计、宏观规模规律、种间资源竞争。

三、Murray 锥度的推导与意义

3.1 权衡的数学形式

对一条承载流量 $q$ 、长度 $L$ 、半径 $r$ 的菌丝边：阻力 $R=L/r^4$ ，材料（生物量）微元 $dV=\pi r^2 L$ 。在固定总送达 $D$ 下，目标是 $\min \sum \pi r_e^2 L_e$ 且各边流量由拓扑决定。由变分可得最优半径满足
$r_e \propto q_e^{1/3}$
即 Murray 定律。其几何直觉是：阻力对半径极敏感（四次方），故高流量主干必须适度加粗以"压住"阻力；而分支流量小，保持纤细即可，避免材料浪费。

3.2 锥度 vs 等粗的定量对比

图1 给出 $r=q^{1/3}$ 曲线（凹），对照线性 $r=q$ ：凹函数意味着把材料从低流量边"转移"到高流量边能同时降低阻力与总材料。基准场景下，若强制全边等粗（取主干 Murray 半径 $r_0=1.303$ ），生物量达 $1112.5$ ；而 Murray 锥度仅 $517.9$ （图2），节省 $53.5\%$ 。图3 显示代价：等粗方案因全面加粗，耗散（ $\sum Rq^2=1580$ ）反而低于 Murray（ $3533$ ）——即以约 2.2 倍的耗散为代价，换来一半的材料。在自然界"建网材料稀缺、养分供给充足"的设定下，这一交换是划算的。

从能量角度再审视这一权衡：等粗方案把全部边加粗到主干级别，使每条支线的阻力都极小，故总耗散最低；但它浪费了大量材料在"本不需要粗"的末梢上。Murray 锥度则精确地把材料投在"刀刃"上——主干承担了绝大部分流量，故值得加粗；末梢流量小，细一点阻力虽高却对全局耗散贡献有限。这种"好钢用在刀刃上"的资源分配思想，不仅在生物网络中存在，也是工程管理、通信网络设计的通用智慧。我们甚至可以给出一个定量判据：当材料成本权重高于能耗权重时，Murray 锥度恒优于等粗方案，这也与真菌"优先节省生物量"的生存策略一致。

四、规模缩放规律

4.1 总长度的次线性增长

令蘑菇在自然地块中随机散布，统计 $N=4,7,\dots,28$ 时的最小生成树总长度。结果 $L$ 从 $117.7$ 增至 $364.7$ ，但对 $N$ 做双对数拟合得 $L\propto N^{0.53}$ （ $R^2=0.95$ ）——显著次线性。这是因为新增蘑菇大多就近接入已有主干，共享一段"干线"，使总延伸被摊销，而非与 $N$ 同步线性外扩。

4.2 生物量与每株成本

总生物量 $B$ 拟合为 $B\propto N^{0.96}$ （ $R^2=0.97$ ），近乎线性（图5）。这看似与"次线性"矛盾，实则源于 Murray 凹性：主干流量随 $N$ 增长使主干半径按 $N^{1/3}$ 增粗，抵消了长度的次线性收益。然而每株成本 $B/N$ 始终保持有界（图6），在大 $N$ 时回落至约 $47$ – $49$ ，呈轻微 U 形——菌落扩大不会造成人均材料成本失控。图7 汇总两指数： $0.53$ （长度）显著小于 $0.96$ （生物量），印证"共享主干摊销"主要压低的是物理延伸而非材料总量。

上述缩放规律具有重要的生态学含义。若生物量随 $N$ 超线性增长（指数大于 1），则大菌落将面临"规模不经济"，生长受阻；而本模型指数 $0.96<1$ ，说明菌丝网络具备天然的规模经济性——尽管不强烈，但足以支撑菌落持续扩张而人均成本不恶化。U 形每株成本曲线进一步提示：存在一个"最优规模"区间（本例约 $N=16$ – $22$ ），此时共享主干的摊销收益与主干加粗的材料代价达到最佳平衡。这对理解真菌菌落为何常以特定密度生长提供了定量线索。

五、有限食物源下的共享竞争

5.1 源通量约束

前文假设食物源供给充足。现实里倒木释放养分的通量有限。设单菌（6 株）需求 $d_1=6.169$ ，令源最大通量 $\Phi=d_1$ （恰好供一菌吃饱）。若两菌共享同一源，有限通量按需求比例分配：
$\text{菌}_i\text{获得} = \Phi\cdot \frac{d_i}{d_1+d_2}$
这一设定把问题从"如何建网最省"切换到"资源总量受限时如何分配"，恰好对应真菌生态中常见的"同一倒木被多菌落争夺"的情景。它与前文的建网优化形成互补：建网是内部效率问题，共享是外部竞争问题，二者共同决定菌落的净收益。

5.2 共享即腰斩

计算得两菌需求 $d_1=6.169, d_2=5.728$ ，共享后分别仅得 $3.199$ 与 $2.970$ ，获取量惩罚约 $48\%$ （图8）。对照之下，两菌各自独立网络的生物量 $B_1=371.7, B_2=271.5$ 表明"建网"本身成本有限——真正的瓶颈是源通量。生态启示：相邻菌落应避免挤占同一倒木，或发展更高效的跨菌落转运，否则将陷入"共同贫穷"。

这一结论还可以用博弈论的语言重新表述：在源通量固定的前提下，两菌落陷入了典型的"公地悲剧"——每个菌落都有动机抢占更多源连接，但 collectively 任意分配都使总获取不超过 $\Phi$ ，且越分散越亏。打破困局的唯一出路是扩大资源基数（寻找新倒木）或提升转运效率（建立菌落间通道以摊薄单源压力）。这与人类社会的公共资源治理逻辑如出一辙，再次彰显了该模型的跨领域解释力。

六、讨论

本篇把"为什么是 Murray 锥度""规模如何缩放""共享如何竞争"三条线索串起：微观上锥度是最优材料分配，宏观上长度次线性使扩展经济，而种间共享的瓶颈在源而非网。这些结论为第三篇的拓扑优化优先级与鲁棒性分析提供了定量支撑。

值得强调的是，三组团规律之间存在内在张力，也蕴含设计智慧。锥度节省了材料，却使耗散略升；规模扩展压低了长度，却因主干增粗使材料近线性增长；共享提升了空间利用率，却因源通量受限而腰斩个体收益。建模的价值，正是把这些彼此牵制的因素同时摆上台面，让决策者看清"省了什么、牺牲了什么"。对真菌而言，其演化显然在"材料最省"与"送达足够"之间找到了稳健的平衡点——而这正是我们模型所复现的均衡态。后续第三篇将进一步检验：当拓扑最优、锥度合理、规模适中时，网络是否还能抵御失效与扰动。

七、结论

我们得出：Murray 锥度比等粗方案省 $53.5\%$ 材料；菌丝长度随规模次线性增长（指数 $0.53$ ）、生物量近线性（指数 $0.96$ ）、每株成本有界；有限源共享使每菌获取量腰斩（惩罚 $48\%$ ）。全部数字由附录 Python 独立复现。

需要说明，本文的规模缩放基于"蘑菇均匀随机散布于固定地块"的设定，其幂指数 $0.53/0.96$ 是这一空间假设下的经验值；若蘑菇呈极度聚集或规则栅格分布，指数会相应变化，但"长度次线性、材料近线性、人均成本有界"的定性结论保持稳健。这种对设定敏感的定量、对机制不敏感的结论，正是物理建模应有的品格。

从更宏观的视角看，本篇揭示的缩放规律其实是一类普适现象：凡是以"树状层级"组织资源输送的系统——真菌菌丝、动物循环系统、电力系统、供水管网——都面临"总长 vs 材料 vs 阻力"的三元权衡，且都倾向于 Murray 型锥度与某种次线性规模律。这使得我们的真菌模型具有了超越单一赛题的方法论价值：它是一把可迁移的"尺子"，用来衡量任何层级化输运网络的效率边界。这也正是数学建模的魅力——从一个具体生物问题出发，最终触及的是普适的结构原理。

参考文献

[1] COMAP. MCM 2021 Problem A: Fungi. 2021.
[2] Murray C D. The physiological principle of minimum work. 1926.
[3] Banavar J R, et al. Supply–demand balance and metabolic scaling. PNAS, 2010.

附录：核心 Python 实现（可独立运行复现上述数字）

import os, sys
_HERE = os.path.dirname(os.path.abspath(__file__))
sys.path.insert(0, os.path.abspath(os.path.join(_HERE, "..", "..", "..", "tools")))
import gen_mcm2021a as G

D = G.gen_mcm2021a()
print("=== MCM2021A Murray锥度与规模（Q2）关键结果 ===")
print("Murray 锥度比等粗半径节省生物量 = %.4f" % D["murray_vs_const_biomass_save"])
print("等粗半径 r0 =", round(D["const"]["r0"], 4), " 等粗生物量 =", D["const"]["B"])
print("\n规模缩放（N, 长度L, 生物量B, 每株B/N, 效率E）：")
for s in D["scale"]:
    print("  N=%2d  L=%.3f  B=%.4f  每株=%.4f  E=%.5f" % (
        s["N"], s["L"], s["B"], s["per"], s["E"]))
print("\n拟合：B(N) ∝ N^%.4f (R2=%.4f) ; L(N) ∝ N^%.4f (R2=%.4f)" % (
    D["scale_p"], D["scale_r2"], D["scale_q"], D["scale_r2_len"]))
sh = D["share"]
print("\n有限源共享：d1=%.4f d2=%.4f Phi=%.4f" % (sh["d1"], sh["d2"], sh["Phi"]))
print("共享后菌1=%.4f 菌2=%.4f  惩罚=%.4f / %.4f" % (
    sh["share1"], sh["share2"], sh["penalty1"], sh["penalty2"]))
print("独立网络生物量 B1=%.4f B2=%.4f（建网非瓶颈）" % (sh["B1"], sh["B2"]))