Circularity of Thermodynamical Material Networks: Indicators, Examples, and Algorithms

本文提出利用热力学材料网络（TMNs）将循环经济视为材料网络设计问题，通过图论形式开发了循环性指标，并结合流体与固体材料的数值算例验证了该方法在刻画高动态物料流方面的优势。

要点

这篇论文就像是在给“循环经济”这个宏大的概念，换上了一副精密的“工程眼镜”。

作者 Federico Zocco 觉得，以前大家谈论循环经济（比如回收、再利用），往往像是在看一张静态的“快照”或者统计报表，只知道大概有多少东西进来了、多少东西出去了。但这就像只看一张全家福，却不知道家里每个人下一秒在做什么。

为了更精准地设计循环经济，作者提出了一套叫**“热力学材料网络”（TMN）的新方法，并发明了一套“循环度指标”**。

我们可以用几个生动的比喻来理解这篇论文的核心内容：

1. 从“看照片”到“拍电影”：动态 vs. 静态

• 旧方法（MFA）： 就像是在看一张全家福照片。你知道家里有多少人，每个人手里拿着什么，但你看不到他们下一秒是去厨房拿水，还是去客厅看电视。这种方法数据量巨大，而且很难捕捉到那些“转瞬即逝”的变化（比如一分钟内发生的流动）。
• 新方法（TMN）： 就像是在拍一部高清电影。作者把供应链看作是一个由“房间”（热力学隔间）和“走廊”（管道）组成的复杂建筑。
  • 房间里存放着材料（库存）。
  • 走廊里流动着材料（物流）。
  • 作者用微分方程（数学里的动态语言）来描述这些材料是如何像水流一样，或者像快递包裹一样，在房间里进进出出的。哪怕是一分钟内的变化，都能被精准捕捉。

2. 什么是“循环”？——寻找“闭环”

在论文里，“循环”不仅仅是指“把东西扔了再捡回来”，而是指**“物质流形成了一个完美的闭环”**。

• 比喻： 想象一个游乐场。
  • 线性经济：你买票进去（输入），玩一圈，然后被保安赶出来扔进垃圾桶（输出）。
  • 循环经济：你买票进去，玩完一圈，自动滑回起点，再玩第二圈，永不停歇。
• 作者的工具（图论）： 作者用**“有向图”**（一种数学网络图）来画这个游乐场。
  • 节点是房间。
  • 箭头是流动的方向。
  • 循环度指标就是用来计算：在这个网络里，有多少箭头构成了**“死循环”**（即箭头转了一圈又回到了原点，没有漏出去）？

3. 六大“体检指标”：给网络做 CT

作者发明了一套指标，就像给这个复杂的材料网络做CT 扫描，看看它“循环”得有多好：

• 几何/算术/调和平均值（Geometric/Arithmetic/Harmonic Means）： 这就像是在计算**“循环的强度”**。如果循环里的每一段路（流量）都很顺畅，那循环度就高；如果有一段路堵死了，整个循环的“平均速度”就会受影响。
• 缩放指标（Scaled）vs. 总量指标（Total）：
  • 缩放指标：就像看百分比。比如，“在这个系统里，有多少比例的材料是在完美循环的？”（0% 到 100%）。
  • 总量指标：就像看总吨数。比如，“到底有多少吨材料在转圈圈？”
• 连通性（Connectivity）： 就像看交通网有多发达。每个房间有多少条路进出？路越多，网络越灵活。
• 方向性（Directionality）： 就像看单行道。如果大部分箭头都指向同一个方向（比如从 A 到 B 到 C），那它可能更像一条流水线，而不是一个循环。

4. 两个生动的案例

案例一：流体（像水一样流动）

• 场景： 想象一个自来水厂的循环系统。
• 问题： 作者先展示了一个“错误”的模型，就像水管漏水了，但水表却显示水量不变（违反了质量守恒定律）。这时候算出来的“循环度”是假的，就像在漏水的桶里谈“水循环”。
• 修正： 作者修正了模型，让水量的变化符合物理定律（流进去多少，存多少，流出来多少）。修正后，那些“循环度指标”才变得真实可信。
• 启示： 如果数学模型不遵守物理定律（比如物质凭空消失或产生），算出来的环保指标就是“数字游戏”。

案例二：固体（像快递包裹一样运输）

• 场景： 想象一个塑料瓶的回收之旅。
  • 工厂生产瓶子 -> 卡车运到餐厅 -> 餐厅用完 -> 卡车运到回收厂 -> 回收厂处理 -> 卡车运回工厂。
• 难点： 固体不像水，水是连续流动的，而固体是**“一批一批”**（Batch）运输的。
  • 卡车 A 运走瓶子时，卡车 B 可能还没回来。
  • 这就导致在同一时刻，可能没有形成“闭环”。
• 发现： 作者发现，在这个固体案例中，因为卡车是分时运行的（去的时候没回，回的时候没去），所以在任何瞬间，网络里都没有形成“闭环”。
  • 结果：所有的“循环度指标”瞬间变成了0！
• 深层含义： 这揭示了现实中的一个痛点。如果物流安排得太死板（一辆车跑完一圈再跑下一圈），循环经济在数学上看起来就是“不循环”的。
• 解决方案： 要想提高循环度，必须让多辆车同时运行，或者让物流更频繁，确保在任何时刻，都有材料在“转圈圈”。

5. 总结：为什么这很重要？

这篇论文的核心贡献在于：

1. 拒绝“拍脑袋”： 以前衡量循环经济可能靠经验或简单的统计，现在有了严格的数学公式和物理基础。
2. 看见“微观动态”： 它能捕捉到1 分钟内发生的材料流动变化，这是传统方法做不到的。
3. 指导设计： 它不仅仅是在“评价”一个系统好不好，而是告诉工程师怎么设计才能更好。比如，如果你发现循环度是 0，那就去调整物流计划，让卡车同时跑起来，而不是等一辆车回来再发下一辆。

一句话总结：
作者用**“数学 + 物理”**的方法，把循环经济从一张模糊的“概念图”，变成了一台可以精确计算、实时监测的“精密仪器”，告诉我们如何把地球的“资源大循环”转得更顺畅、更真实。

技术摘要

热力学材料网络（TMN）的循环性：指标、示例与算法

技术总结

1. 研究背景与问题陈述

传统的“开采 - 制造 - 废弃”线性经济模式在资源有限性和环境污染方面不可持续。循环经济（Circular Economy, CE）旨在通过增加重复使用、维修、翻新和回收来闭环材料流。然而，现有的循环经济指标存在两大主要局限：

1. 缺乏数学和物理基础：许多现有指标仅基于代数方程或定性描述，缺乏清晰的物理定义和可重复性，难以进行系统比较。
2. 忽视生命周期闭环程度：许多指标仅关注单一生命周期阶段（如回收阶段）或经济价值，未能从系统层面衡量材料流在整体网络中的闭环程度。
3. 静态模型的局限性：传统的物质流分析（MFA）通常基于静态（恒定）的库存和流量数据，难以捕捉现实中发生频率极高（如小于 1 分钟）的动态变化。

核心问题：如何建立一种基于物理原理、能够描述高度动态物质流、并基于图论系统定义循环经济指标的方法？

2. 方法论：热力学材料网络 (TMN)

本文提出了一种基于**热力学材料网络（Thermodynamical Material Networks, TMNs）**的新框架，将供应链视为一组可连接、修改或移除的热力学隔室（Compartments）。

• 核心概念：

  • TMN 定义：由相互连接的热力学隔室组成的集合，用于传输、存储、使用和转化目标材料。每个隔室通过控制面建模，遵循质量平衡和热力学定律。
  • 图论形式化：将 TMN 映射为加权有向图（Digraph）。
    • 节点（Nodes）：代表存储、转化或使用材料的隔室（$R$ 子集），权重为物质库存（Mass Stock）。
    • 弧（Arcs）：代表在隔室间移动材料的隔室（$T$ 子集），权重为质量流量（Mass Flow Rate）。
  • 动态建模：
    • 流体材料：使用连续时间微分方程（基于质量守恒原理 $\frac{dm}{dt} = \sum \dot{m}_{in} - \sum \dot{m}_{out}$）。
    • 固体材料：使用混合动力学方程（Hybrid Dynamical Equations），结合连续时间和离散时间项，以模拟批次运输（Batch Transport）和离散事件。

• 循环性指标体系：
  基于图论中的**有向循环（Directed Cycles）**概念，作者定义了三类核心循环性指标和四类辅助指标：

  1. 基于循环的指标（核心）：
     • 几何/调和/算术平均缩放循环性 ($\lambda_{GS}, \lambda_{HS}, \lambda_{AS}$)：归一化指标（范围 [0,1]），衡量网络中闭环流动相对于总流动的强度。
     • 几何/调和/算术平均总循环性 ($\lambda_{GT}, \lambda_{HT}, \lambda_{AT}$)：非归一化指标，衡量循环流动的绝对强度。
     • 循环共享 ($\lambda_S$)：衡量被多个循环共享的流量。
     • 循环性 ($\lambda_Y$)：网络中有向循环的数量。
  2. 拓扑辅助指标：
     • 平均连通性 ($\lambda_C$)：节点的平均度数。
     • 方向性 ($\lambda_D$)：衡量流动方向与节点索引顺序的一致性。
  3. 物质分布辅助指标：
     • 总库存 ($\theta_S$)、总流量 ($\theta_F$)、库存分布 ($\theta_D$)、累积 - 耗尽向量 ($\theta_A$)。

3. 主要贡献

1. 物理基础的形式化：首次利用图论和热力学原理，为材料循环性提供了严格的数学定义和物理基础，解决了现有指标缺乏可重复性和物理意义的问题。
2. 动态指标开发：提出了能够处理静态、连续时间动态和混合动态（离散 + 连续）库存与流量的指标计算算法。
3. 算法实现：开发了基于 MATLAB 的算法（公开代码），能够自动计算上述所有指标，包括循环检测、权重计算和动态更新。
4. 多场景验证：通过流体和固体材料两个不同领域的数值案例，验证了方法的适用性。

4. 关键结果与案例研究

案例一：流体材料（动态连续时间）

• 场景：模拟了一个包含 4 个节点和动态流量（正弦/余弦函数）的网络。
• 发现：
  • 如果直接设定静态库存而忽略质量守恒（即 $\frac{dm}{dt} \neq 0$），会导致非物理的循环性指标结果。
  • 修正：通过强制满足质量守恒方程（方程 63），动态更新主对角线（库存）的数值，使得模型物理上合理。
  • 结果：在修正后的模型中，当存在同时非零的循环流动时，缩放循环性指标（$\lambda_{GS}$ 等）达到 1（完全闭环）；当网络拓扑随时间变化（流量开关）时，指标能实时反映循环性的动态变化。

案例二：固体材料（混合动力学，塑料回收）

• 场景：模拟了一个涉及 PET、HDPE 和 PP 塑料的“组织 - 制造商 - 回收中心”网络。包含订单、运输（卡车）、回收和再配送过程。
• 建模特点：使用了混合动力学方程（方程 74-86），模拟了批次运输、交通拥堵导致的延迟以及污染导致的不可回收损耗（30%）。
• 关键发现：
  • 循环性为零的悖论：尽管网络拓扑上存在回路（组织 $\to$ 回收 $\to$ 组织），但由于固体运输是批次化的，去程（$m_{1,2}$）和回程（$m_{2,1}$）在时间上不重叠（即不存在 $t$ 使得 $\dot{m}_{1,2}(t) \neq 0$ 且 $\dot{m}_{2,1}(t) \neq 0$）。
  • 指标结果：所有基于循环的指标（$\lambda_Y, \lambda_{GS}$ 等）均为 0。
  • 物理意义：这揭示了该指标体系的一个关键特性——它严格衡量同时存在的闭环流动。在现实中，如果运输频率低（如一天一班），系统在这些指标下表现为“非循环”，尽管物质最终回到了起点。
  • 数据效率：TMN 模型仅需 34 个参数 即可精确模拟 290 分钟内的所有动态事件。相比之下，若使用传统 MFA 以 1 分钟分辨率捕捉相同事件，需要 2,601 个历史数据点（3 种材料 $\times$ 3 个节点 $\times$ 289 个时间点）。TMN 的数据需求减少了约 76 倍。

5. 意义与结论

• 理论突破：将循环经济从定性的“概念”转化为可量化、可计算的“网络设计问题”，类比于热力学循环或电路网络的设计。
• 动态捕捉能力：TMN 能够捕捉亚分钟级的动态变化（如卡车到达、批次装载），这是传统静态 MFA 无法做到的。
• 设计指导：指标结果（如固体案例中循环性为 0）直接指出了系统设计的改进方向：为了真正提高循环性，必须增加运输频率或优化调度，使得去程和回程流量在时间上重叠，从而形成物理上的“同时闭环”。
• 未来展望：
  • 扩展模型规模以应用于更复杂的工业网络。
  • 改进指标定义，使其不仅依赖于“同时”流动，而是考虑固定时间窗口（如一天）内的循环流动，以更贴合实际运营场景。

总结：本文提出了一种基于物理和图论的强有力框架，通过热力学材料网络（TMN）量化循环经济。它不仅解决了现有指标缺乏物理基础的问题，还通过动态建模揭示了传统静态分析无法发现的系统特性（如时间错位导致的“伪线性”），为设计真正高效的闭环材料系统提供了数学工具和算法支持。