关于《tissue‑growth‑based synthetic tree generation and perfusion simulation》的学术研究报告
一、 研究作者、机构与发表信息
本研究的主要作者为Hyun Jin Kim, Hans Christian Rundfeldt, Inpyo Lee, 以及 Seungmin Lee。其中，Hyun Jin Kim为通讯作者。研究团队来自韩国科学技术院（Korea Advanced Institute of Science and Technology, KAIST）机械工程系，第二作者亦与德国卡尔斯鲁厄理工学院（Karlsruhe Institute of Technology）有关联。该研究成果以原创论文（original paper）形式发表于期刊《biomechanics and modeling in mechanobiology》2023年第22卷，论文于2023年2月10日被接受，并于2023年3月4日在线发表。
二、 研究学术背景
本研究的核心科学领域是计算生物力学与生物医学工程，具体聚焦于血管网络生成与多尺度血流灌注模拟。
研究动机与背景：基于医学影像的血流模拟技术在临床应用中前景广阔，但面临重大挑战。其中一个关键难题是如何在影像分辨率之外，对患者特异性的血管网络进行建模。这些微观血管网络对于组织获得氧气和养分至关重要，它们与组织之间存在着“供需关系”。为了模拟组织灌注并预测缺血风险，需要生成既符合生理学原理（如最小能量消耗、满足组织代谢需求）又能在合理时间内完成计算的血管网络。现有的主流方法，如约束性构造优化（Constrained Constructive Optimization， CCO）算法，虽能生成接近最小能量消耗的血管树，但存在计算成本高、难以并行化、且在多个源血管竞争时可能产生不现实的血管网络等局限。
研究目标：本研究旨在开发一种新的合成血管树生成算法。该算法应能：1）基于组织生长与血管生长的互动，生成解剖学上真实、符合生理学原理（如供应附近组织、满足代谢需求、遵循最小功原理）的血管网络；2）大幅降低计算耗时，并能有效并行化；3）与灌注模拟框架结合，用于量化组织灌注并预测患者特定几何结构中的缺血风险区域。
三、 研究详细工作流程
本研究主要包含三大核心流程：1）基于组织生长的合成血管树生成算法开发；2）多尺度血流与灌注模拟框架构建；3）患者特异性几何模型的构建与算法验证。
流程一：基于组织生长的合成血管树生成算法 该方法模拟了组织生长引发血管新生的过程，其工作流程可概括为“初始化 -> 迭代生长与优化”。 1. 初始化：从医学影像中分割出主要动脉作为初始血管树 T_init。将其所有末端（出口）定义为潜在的“源点”。将这些源点的半径缩放至用户定义的最小血管半径 r_min。利用Murray定律（或幂律）估算每个源点的初始血流 q_i ∝ r_min^k。随后，根据代谢需求反推该初始血流所能供应的初始组织体积 ω_init，即从离血管出口最近的组织开始“占据”，直到被占据组织的总代谢需求与该出口血流匹配。 2. 迭代生长与优化：这是一个循环过程，直到组织生长至目标组织 ω_t。 * 组织生长：在每次迭代中，当前已灌注的组织 ω_i 会按照预设的生长速率向量向目标组织 ω_t 生长。生长基于Voronoi图划分，每个当前血管树 T_i 的末端负责其附近一个组织分区的生长。 * 代谢需求计算：生长后的新组织 ω_{i+1} 的代谢需求总量 q_{ω_{i+1}} 被计算出来。与上一迭代组织需求的差值，决定了需要新增多少末端血管段来满足新增的代谢需求。新增末端数量 n_{i+1} = round( (q{ω{i+1}} - q_{ω_i}) / q_min )，其中 q_min 是单个末端段的最小流量。 * 血管生长与优化：为每个新的末端位置，在现有血管树的邻近区段中寻找最佳分叉位置。通过构建候选血管树 T_{c,i}，并计算其损失函数 L，来评估和选择最优结构。损失函数定义为：L = Σ ( 8μ l_p / π r_p^4 * q_p^2 + β_p π r_p^2 l_p )，其中第一项代表血流粘性功，第二项代表血管代谢维持功。优化目标是最小化此损失函数。在此过程中，算法强制遵循血管分叉的几何约束（如幂律 r0^k = r1^k + r2^k）和血流-半径关系（qp ∝ rp^k）。源点处的血管半径会随着迭代逐渐增长，直到达到影像分割得到的原始半径值。 * 并行化：该算法的关键优势在于其固有的可并行性。由于每个组织分区的生长和新末端血管的生成优化可以独立进行，因此能显著减少计算时间。
流程二：多尺度血流与灌注模拟框架 生成的血管树被用于模拟组织中的血液灌注，采用显式耦合的多尺度方法。 1. 一维血管网络血流模拟：在生成的血管树（包括分割的主干和合成的分支）中，使用简化的一维Navier-Stokes方程模拟血流和压力。假设血管壁不可变形，并应用质量守恒和线性动量平衡方程。在血管树入口施加流入边界条件，在末端施加阻力或Windkessel模型边界条件。 2. 三维组织灌注模拟：将组织视为多孔介质，采用Darcy定律模拟小动脉和毛细血管水平的血液灌注。通过求解Darcy速度向量和毛细血管压力的方程，并设置源项和汇项来分别代表血液从血管末端进入组织以及从组织静脉端流出的过程。 3. 显式耦合：在每一个时间步内，进行迭代耦合。将一维模型计算得到的血管末端压力作为Darcy模型的源压力输入；同时，将Darcy模型计算得到的每个组织单元流回血管网络的流量，作为一维模型末端的流量边界条件反馈回去。如此迭代直至两者收敛。
流程三：患者特异性几何构建与验证流程 研究使用ITK-Snap、3D Slicer、SimVascular等开源软件从医学影像数据中分割和重建患者特定的血管腔和组织几何模型，并生成用于模拟的三维和二维离散化网格。 验证分为两个层面： 1. 理想模型验证：在一个球体组织中，设置两个对称的源点。 * 不同分辨率：生成不同终端数量的血管树（567 到 24，493个末端），验证算法在不同树密度下的稳定性和灌注均匀性。 * 非均匀代谢需求：设置径向递减和径向递增两种代谢需求分布，验证算法能根据组织代谢需求自适应调整血管结构。 * 与CCO算法对比：在相同条件下，与传统的约束性构造优化（CCO）算法在计算时间和生成的血管总体积上进行直接比较。 2. 患者特异性模型演示： * 脑血管模型：应用算法于一组脑血管动脉和大脑半球，假设脑组织代谢需求均匀。 * 冠状动脉模型：应用算法于一组冠状动脉和左右心室，为左心室（高代谢需求）和右心室（低代谢需求）设置不同的均匀代谢率，以演示算法对不同生理条件的响应。
四、 研究主要结果
理想模型结果： 1. 不同分辨率下的表现：生成的血管树对称，两个源点的灌注体积和流量基本相等，符合预期。组织血流总体均匀，仅在与源点距离较远处因压力损失略有下降。随着末端段数量增加，合成的血管树总体积相应增加。 2. 非均匀代谢需求下的适应：算法成功生成了适应不同代谢模式的血管树。对于径向递减的代谢需求，算法在组织中心生成了管径较大的血管结构，而外围血管相对萎缩。对于径向递增的代谢需求，则呈现相反趋势，中心血管较细，向外血管增粗（但由于末端位置依赖于网格，在最外围又变细）。尽管代谢需求分布不同，两个对称源点的总流量和灌注区域仍保持平衡，证明了算法的鲁棒性。 3. 与CCO算法的对比：这是本研究的关键结果。在生成相同数量末端血管时，生长基优化算法（GBO）的计算时间较CCO算法大幅减少（例如，生成5003个末端时，GBO耗时17.3分钟，CCO耗时457.4分钟，加速比达26倍）。两者生成的血管树总体积非常接近，差异小于5%，表明GBO在极大提升效率的同时，并未牺牲血管网络的“经济性”（即最小功原理）。此外，GBO算法生成的灌注区域在不同分辨率下更稳定，且由于其“就近生长”的原则，避免了CCO算法有时会出现的血管“舍近求远”供应远处组织的不现实情况。
患者特异性模型结果： 1. 脑血管模型：成功为58个出口生成了包含9713个末端的合成血管树。灌注模拟显示，在一个健康脑中，组织血流分布非常均匀，符合生理预期。 2. 冠状动脉模型：成功为左右心室生成了不同结构的血管树。左心室因代谢需求高、体积大，生成了更密集、管径相对更大的血管网络（血管体积11.8 mL）；右心室则生成了较稀疏、管径更小的网络（血管体积0.7 mL）。生成的灌注区域图显示，每条分割的冠状动脉主要供应其邻近的心室组织，这与已知的生理学和临床观察相符。
五、 研究结论与价值
本研究成功开发并验证了一种基于组织生长的合成血管树生成算法及配套的多尺度灌注模拟框架。
结论：该算法能够依据组织代谢需求，快速生成符合生理学原理（供应附近组织、满足代谢需求、遵循最小功原理）的血管网络。与广泛使用的CCO算法相比，在生成质量相近血管网络的前提下，计算效率获得数量级提升，且具备良好的并行扩展能力。将生成的血管网络与多尺度灌注模型耦合，能够在患者特异性几何结构上实现从大动脉到组织毛细血管的血流模拟。
价值： * 科学价值：提出了一种新的血管生成建模范式，将组织生长动力学与血管优化生长紧密结合，为理解血管网络形态发生提供了新的计算模型。 * 应用价值：该框架具有直接的转化医学潜力。可用于构建患者特异性、全心脏或全脑的精细血流动力学模型，辅助评估冠状动脉疾病或脑血管疾病中特定血管阻塞对下游组织灌注的影响，量化缺血风险，为个性化医疗决策（如手术规划、介入治疗评估）提供定量工具。
六、 研究亮点
七、 其他有价值的讨论
作者在讨论部分也坦诚指出了当前研究的局限性，为未来工作指明了方向： 1. 算法假设分割的血管是组织生长的唯一起点，而现实中可能存在其他“不可见”的源点（如穿支动脉）。 2. 目前末端血管段的位置依赖于组织体积网格的质量，可能导致血管树密度分布不均。尽管血管直径由代谢需求决定，在一定程度上抵消了此影响，但仍需改进。 3. 采用的单室达西灌注模型存在局限性，未来需要引入多室模型以更好地表征不同尺度血管的异质性。 4. 需要更多利用离体解剖数据或在体高分辨率影像数据进行验证。
这些讨论体现了研究的严谨性，并清晰地规划了未来的发展路径。