一种虚拟肝静脉压力梯度的测量方法
【专利摘要】本发明属于早期无创性诊断领域，涉及一种虚拟肝静脉压力梯度的测量方法。本发明所述的一种虚拟肝静脉压力梯度的测量方法包括：肝静脉?门静脉系统三维建模；有限元网格划分数学模型；以及流体力学仿真计算虚拟肝静脉压力梯度(vHVPG)。本发明优化并完善肝静脉?门静脉系统三维建模、有限元网格划分和流体力学仿真计算，构建并验证一种更具诊断优势的vHVPG检测新技术，为门脉高压患者的早期无创诊断提供一种安全无创、准确量化的新途径。
【专利说明】
-种虚拟肝静脉压力梯度的测量方法
技术领域
[0001] 本发明属于早期无创性诊断领域，设及一种虚拟肝静脉压力梯度的测量方法。
【背景技术】
[0002] 在全球范围内，肝脏疾病(病毒性感染、酒精性肝病、非酒精性脂肪肝及相关的肝 硬化、肝细胞肝癌）是导致死亡的主要原因之一，仅在中国肝脏疾病就威胁了近3亿人的健 康和生命，严重加剧了全球疾病的负担。W乙型肝炎病毒感染引起的慢性肝病为例，据世界 卫生组织报道，全球约20亿人曾感染乙肝病毒，其中2.4亿人为慢性乙肝病毒感染者，每年 约有65万人死于乙肝病毒感染所致的肝硬化、肝细胞肝癌等。在全球肝硬化患者中，由乙肝 病毒感染引起的为30%。而在我国，由于感染乙肝病毒而导致的肝硬化比例高达60%。因 此，肝硬化的高发病率和较差的临床结局已造成了国内外严重的社会公共健康问题。
[0003] 口脉高压是肝硬化失代偿期的重要表现之一，由于其本身并无症状，临床上口脉 高压的发现和诊断常由于出现晚期严重的并发症，如：急性静脉曲张破裂出血、顽固性腹 水、肝性脑病、Π 脉高压性胃肠病、肝肾综合症、肝肺综合症和继发感染等。因此，Π 脉高压 及其并发症严重影响了肝硬化患者的生活质量和长期预后，对于口脉高压(尤其是早期无 症状性口脉高压）的诊断和监测是肝硬化等终末期肝病治疗链中最重要的环节之一。
[0004] 在最新的2015版Baveno VI共识（Π 脉高压的风险分层和个体化管理）中，有创性 肝静脉压力梯度化epatic venous pressure gradient,HVPG)的测量继续被推荐为诊断临 床显著性口脉高压的金标准，即HVPG>10mmHg可确诊（de Franchis R#*,Baveno VI Faculty.Expanding consensus in portal hypertension:Report of the Baveno VI Consensus 胖orkshop:Str曰tifying risk 曰nd individu曰lizing c曰re for port曰1 hypertension. J Hepatol. 2015 ;63(3): 743-752.)。该侵袭性方法是通过颈内静脉穿刺置 管，依次经过颈内静脉、上腔静脉、右屯、房、下腔静脉进入肝静脉，分别测量肝静脉自由压 (free hepatic venous pressure,FHVP)和肝静脉模入压（wedged hepatic venous pressure,WHVP)，计算两者差值得到HVPG[Bloom S#,Kemp W,Lubel J*.Po;rtal hypertension:pathophysiology.diagnosis and management. Intern Med J.2015:45 (1): 16-26.]。此外，HVPG对于治疗应答、并发症风险和长期预后也可提供重要信息，其数值 的变化被推荐为研究预后的替代指标。此外，HVPG测量还被鼓励用于研究口脉高压创新诊 断和治疗方法的临床试验中。然而，测量HVPG也存在一些问题，首先其检测手段为有创性， 在口脉高压早期尚无严重并发症时很难被患者所接受;其次，测量过程需由经过专业化培 训的介入医生在肝静脉造影的辅助下进行，运不但增加了受试者的射线暴露和造影剂过敏 可能，还存在有一定的技术操作风险;加之其高昂的诊断费用，HVPG测定在我国受到很大限 审IJ，目前仅有部分Ξ级甲等医院开展。除此之外，Π 脉高压的有创性诊断技术还包括:超声 引导下口脉穿刺测压和开腹手术中直接口脉测压[de Franchis R#*,Dell'Era A. Invasive and noninvasive methods to diagnose portal hypertension and eso地ageal varices.Clin Liver Dis.2014;18(2):293-302.]。然而，直接测量口脉压力 的风险更大、对操作者的技术水平要求也更高，在临床实践中很难被医生和患者所采纳，目 前仍多用于动物实验性研究。
[0005] 口脉压力的无创性评估是目前该领域研究的重要内容，主要包括W下Ξ个方面： (1)检测肝内阻力:肝内阻力增加是导致口脉高压形成的重要因素。内源性缩血管物质的释 放和肝组织结构的素乱可促使肝内阻力的显著升高。有研究提示，血清内皮素-1表达水平 与HVPG存在相关性，但该结论尚缺乏大样本量临床试验的支持。另一方面，肝内机械结构的 改变也会导致肝内阻力增加，临床上多采用肝瞬时弹性成像(門broScan)技术，通过对肝脏 硬度值的测定来评估肝硬化程度，该方法操作简便、重复性好。但也有报道指出，門broScan 的诊断准确性受到肥胖、肋间隙狭窄、炎症等因素干扰，对口脉高压并发症(如高风险食管 静脉曲张)的预测价值不高，需联合其他无创指标进行综合判断。（2)测定口脉循环血量：Π 脉高压的直接因素为循环血量的积聚升高超过了机体的代偿，当肝硬化伴发重度口脉高压 时，屯、输出量也会明显升高且与HVPG存在高度相关性。有研究利用彩色多普勒超声(color Doppler ulhasound ,CDUS)对肝硬化患者的口脉直径和血流速度进行测定，在此基础上计 算口脉系统的循环血量。然而，CDUS测量技术的主观性较强且易受到肥胖等因素干扰，尚无 准确的临床证据表明经血流速度和口脉直径计算的口脉压力与HVPG具有较高的诊断一致 性。此外，CT血管成像(CT angiogra地y，CTA)和磁共振静脉显像也可W评估口脉系统的循 环血量。然而，其测定原理仍是基于口脉直径和侧枝循环来间接评价口脉压力。（3)仿真计 算口脉压力:虚拟仿真技术是计算机图形学、医学图像处理、软件工程计算等多学科交叉融 合的结晶。W计算机仿真辅助医学影像获得的虚拟人体数据为疾病的介入诊断和治疗提供 了全新的思路。目前已有多项国际前瞻性临床试验表明，基于CTA和流体力学计算获得的血 流储备分数对于冠状动脉的功能性狭窄具有高度诊断价值，相关仿真计算软件也已获得美 国食品药品监督管理局的上市批准。
[0006] 口脉压力的各种测定技术目前仍受到各自干扰因素的影响，或者有创风险大、操 作难度高，或者干扰因素多、数值变异大。此外，国内外该领域研究的主要问题还表现在:评 估范围局限于口脉主干，未考虑肝静脉-Π 静脉系统的整体影响;仅考虑到口脉血液动力学 因素，忽略了肝硬化患者的血液流变学改变，如血液黏度的变化和红细胞变形能力的降低 都会促使口脉阻力的升高。
[0007] 目前国内外相关报道中，口脉高压的诊断方法，特别是早期无创性诊断技术，仍有 待进一步研究完善。
【发明内容】
[0008] 本发明的目的在于提供一种虚拟肝静脉压力梯度的测量方法，W此为基础可W构 建一种更具诊断优势的虚拟肝静脉压力梯度检测新技术，为口脉高压患者的早期无创诊断 提供一种新的途径。
[0009] 本发明所述的一种虚拟肝静脉压力梯度的测量方法，包括W下步骤：
[0010] A.从标本的肘正中静脉注射造影剂，进行CT血管成像(CTA)，获取包括肝静脉期在 内的CTA图层序列，导出图层序列，格式为dicom,层厚1.25mm，图像分辨率512 X 512像素；
[0011] B.将所获取的CTA图层序列导入医学影像控制软件MIMICS,选择肝静脉期图层序 列，设置图像序列的方位，MIMICS软件自动识别图像序列，生成肝静脉期CTA图像序列的冠 状位、矢状位和水平位图像；
[001^ C.寻找图像中的肝静脉-Π 静脉系统（目标），利用MIMICS软件的Τ'虹esholding(阔 值算法)工具，W尽量包含目标的CT值、尽量排除临近目标的周围肝脏等软组织的CT值为原 贝1J，设定阔值范围，W提取目标;利用MIMICS软件的Region growing(区域增长)工具，选定 目标，W提取只与目标在空间结构上有连接的结构;利用MIMICS软件的化Iculate 3D from mask(3维建模)工具，选择quality(质量)为medium(中等精度），建立初步的肝静脉-口静 脉Ξ维模型；
[0013] D.利用MIMICS软件的Crop mask(裁剪蒙板)工具进一步提取目标结构，剔除部分 非目标结构;再利用MIMICS软件的Edit masks in 3D(S维编辑蒙版)工具，剔除剩余的非 目标结构，只保留肝静脉-Π 静脉系统;反复利用MIMICS软件的Edit masks in 3D(S维切 害d)工具和Edit mask(二维编辑蒙板)工具，进一步对肝静脉-Π 静脉系统进行选择性填充、 剔除噪点像素，从而重建内腔封闭的肝静脉-Π 静脉系统实屯、Ξ维模型；
[0014] Ε.利用MIMICS软件的Smoothing(光滑)操作，将Ξ维模型进行表面光滑处理;将光 滑处理后的实屯、Ξ维模型的几何模型数据选择Ansys area elementUnsys软件面文件)格 式（.inp)导出；
[0015] F.在ANSYS经典模式下，导入所述面文件(格式为.inp)，将长度单位统一为国际单 位m;W面为基础建立肝静脉-Π 静脉系统模型的实屯、模型体；
[0016] G.通过布尔操作对肝静脉-Π 静脉系统模型的血流入口及出口作垂直切面，得到 模拟肝静脉自由压(vFHVP)的开放的几何模型;完成上述操作后，将文件W后缀为.IGS格式 导出备用（IGS是一种Ξ维数值模型文件格式，ANSYS WoriAench模块可读取）；
[0017] H.建立ANSYS Workbench有限元计算平台，包括GeometiT几何模型模块、Fluent流 体计算模块和Results模块（即C抑-POST后处理模块）。通过Geomet巧模块将IGS文件导入， 在Mesh(划分网格）单元中，划分对象为导入的数值模型，网格划分方法设置为 Tetrahe化ons(四面体型），在Physics Preference(物理设置）中选C抑（计算流体力学分 析），在Solver Preference(求解设置）中选Fluent(使用Fluent求解流场）；考虑到运算精 确度W及计算机运行速度，对划分网格尺寸进行限定，max face size(最大面尺寸)设置为 1.5mm，max size(最大尺寸）设为4mm;完成W上设置后，通过Generate Mesh(产生网格)完 成网格划分；
[0018] I.在流体力学计算模块Fluent的solution(解决方案）中设置材料参数(血液密 度、血液粘度、血管壁密度），使模型的物理属性接近人体生物学属性，提高仿真的准确度； 求解控制参数(计算步长、迭代次数、最大循环次数)和边界条件(命名血流入口面，并赋予 速度值，命名出口面后赋予压力值，未命名的血管壁设置为wall), 口静脉血流雷诺数Re< 2000，故仿真流体设置为层流;运算初始化设置为从入口面开始;完成上述参数设定后，模 拟血管壁与血液的流体-固体禪合，计算获得仿真Ξ维血管模型的压力分布和血流分布；
[0019] J.在results后处理模块中，对结果进行读取，通过contour(轮廓)操作显示肝-Π 静脉模型压力分布图；利用软件自带的calculators选项卡，计算获取虚拟肝静脉自由压 (vHWP)数值；
[0020] K.创建一个直径大于或等于被截断血管的圆柱体来模拟阻断球囊，通过布尔运算 对肝右静脉进行阻断，得到模拟肝静脉模入压(vWHVP)的开放的几何模型，并WIGS文件导 出备用（参考步骤F);将IGS文件导入到ANSYS wor化ench中，在产生的两个截面分别赋予 Om/s的速度模拟血流停滞，其余材料参数、边界条件和求解控制参数不变，计算得到虚拟肝 静脉模入压(vWHVP)，在results模块中读取虚拟肝静脉模入压(vWHVP);
[0021] L.求出所述vWHVP与所述vFHVP的差值，即为虚拟肝静脉压力梯度(vHVPG)。
[0022] 根据本发明所述的虚拟肝静脉压力梯度的测量方法的进一步特征，所述步骤I中， 所述材料参数被设置为:血液密度= 1050kg/m3、血液粘度= 0.00化g/m · S、血管壁密度= 1150kg/m\
[0023] 如前所述，口脉高压作为肝硬化失代偿期的重要特征，其临床明确诊断往往滞后 于晚期严重的并发症发生。然而，目前国内外的口脉高压诊断技术或者有创风险大、操作难 度高，或者干扰因素多、数值变异大。本发明首先提出了虚拟肝静脉压力梯度（virtual hepatic venous pressure gradient,vHVPG)的概念，并在其方法学构建中取得重要数据， 有利于探索并建立一种安全无创、准确量化的口脉高压诊断技术。本发明优化并完善肝静 脉-Π 静脉系统Ξ维建模、有限元网格划分和流体力学仿真计算，构建并验证一种更具诊断 优势的vHVPG检测新技术，为口脉高压患者的早期无创诊断提供一种新的途径。
【附图说明】
[0024] 图1为从CT血管成像(CTA)导出包括肝静脉期在内的CTA图层序列。
[0025] 图2为选择导入的肝静脉期CTA图层序列，层厚1.25mm。
[0026] 图3为设置图层序列的方位。
[0027] 图4为自动生成的冠状位、矢状位和水平位图像。
[0028] 图5为化resholding (阔值算法)工具设定阔值范围。
[00巧]图6为Region growing(区域增长)工具选定目标。
[0030] 图7为Cal州late 3D打om mask(S维建模），选择中等精度。
[0031] 图8为建立初步的肝-Π 静脉系统。
[0032] 图9为化op mask(裁剪蒙板)工具进一步提取目标结构，剔除部分非目标结构。
[0033] 图10为化op mask(裁剪蒙板)工具操作后保留的Ξ维模型。
[0034] 图11为Edit masks in 3D(S维编辑蒙版)工具，剔除剩余的非目标结构。
[00巧]图12为Edit masks in 3D(S维编辑蒙版)工具操作后保留的肝静脉-Π 静脉系 统。
[0036] 图13为Edit mask(二维编辑蒙板)工具在冠状位、矢状位和水平位编辑蒙板。
[0037] 图14为反复利用MIMICS的Edit masks in 3D(S维编辑蒙版）工具和Edit mask (二维编辑蒙板)工具操作后的肝静脉-Π 静脉实体模型。
[003引图15为对肝静脉-Π 静脉系统进行Smoothing(光滑)操作。
[0039] 图16为将光滑处理后的模型选择.inp格式导出。
[0040] 图17为W面为基础建立肝静脉-Π 静脉系统模型的实屯、模型体。
[0041 ]图18为模拟肝静脉自由压(FHVP)的开放的几何模型。
[0042] 图19为ANSYS Workbench有限元分析计算平台。
[0043] 图20为solution(解决)单元界面。
[0044] 图21为设定求解控制参数及迭代数目。
[0045] 图22为求解获得仿真Ξ维血管模型的压力分布和血流分布。
[0046] 图23为虚拟肝静脉自由压vFHVP模拟及压力分布图。
[0047] 图24为读取虚拟肝静脉自由压vFHVP数值。
[004引图25为虚拟肝静脉模入压vWHVP数值模型。
[0049] 图26为虚拟肝静脉模入压vWHVP结果及压力分布图。
[0050] 图27为流体力学仿真计算vHVPG与有创测得的压力的比较结果图。
【具体实施方式】
[0051 ]下面结合附图对本发明进行详细说明。
[0052] 本发明所述的一种虚拟肝静脉压力梯度的测量方法包括:肝静脉-Π 静脉系统Ξ 维建模;有限元网格划分数学模型；W及流体力学仿真计算vHVPG。
[0053] 1.从标本的肘正中静脉注射造影剂，进行CT血管成像(CTA)，获取包括肝静脉期在 内的CTA图层序列，导出图层序列，格式为dicom,层厚1.25mm，图像分辨率512 X 512像素。如 图1所示。
[0054] 2.将所获取的CTA图层序列导入医学影像控制软件MIMICS,选择肝静脉期图层序 列（如图2所示），设置图像序列的方位(如图3所示），MIMICS软件自动识别图像序列，生成肝 静脉期CTA图像序列的冠状位、矢状位和水平位图像(如图4所示）。
[0055] C.寻找图像中的肝静脉-Π 静脉系统（目标），利用MIMICS软件的Τ'虹esholding(阔 值算法)工具，W尽量包含目标的CT值、尽量排除临近目标的周围肝脏等软组织的CT值为原 则，设定阔值范围，W提取目标（如图5所示）；利用MIMICS软件的Region growing!；区域增 长)工具，选定目标，W提取只与目标在空间结构上有连接的结构(如图6所示）；利用MIMICS 软件的Calculate 3D from mask(S维建模）工具，选择quality(质量）为medium(中等精 度），建立初步的肝静脉-Π 静脉Ξ维模型(如图7和图8所示）。
[0056] 4.利用MIMICS软件的Crop mask(裁剪蒙板)工具进一步提取目标结构，剔除部分 非目标结构(如图9和图10所示）；再利用MIMICS软件的Edit masks in 3D(S维编辑蒙版） 工具，剔除剩余的非目标结构，只保留肝静脉-Π 静脉系统（如图11和图12所示）；反复利用 MIMICS软件的Edit masks in 3D(S维切割）工具和Edit mask(二维编辑蒙板)工具，进一 步对肝静脉-Π 静脉系统进行选择性填充、剔除噪点像素，从而重建内腔封闭的肝静脉-Π 静脉系统实屯、Ξ维模型。如图13至图15所示。
[0化7] 5.利用MIMICS软件的Smoothing(光滑)操作，将Ξ维模型进行表面光滑处理;将光 滑处理后的实屯、Ξ维模型的几何模型数据选择Ansys area elementUnsys软件面文件)格 式（.inp)导出。如图16所示。
[0化引6.在ANSYS经典模式下，导入所述面文件(格式为.inp)，将长度单位统一为国际单 位m;W面为基础建立肝静脉-Π 静脉系统模型的实屯、模型体(如图17所示）。
[0059] 7.通过布尔操作对肝静脉-Π 静脉系统模型的血流入口及出口作垂直切面，得到 模拟肝静脉自由压(vFHVP)的开放的几何模型(如图18所示）；完成上述操作后，将文件W后 缀为.IGS格式导出备用（IGS是一种Ξ维数值模型文件格式，ANSYS Workbench模块可读 取)。
[0060] 8.建立ANSYS Wor化ench有限元计算平台（如图19所示），包括Geomet巧几何模型 模块、Fluent流体计算模块和Results模块（即C抑-POST后处理模块）。通过Geomet巧模块将 IGS文件导入，在Mesh(划分网格)单元中，划分对象为导入的数值模型，网格划分方法设置 为Tetr址e化ons(四面体型），在Physics Preference(物理设置）中选C抑（计算流体力学分 析），在Solver Preference(求解设置）中选Fluent(使用Fluent求解流场）；考虑到运算精 确度W及计算机运行速度，对划分网格尺寸进行限定，max face size(最大面尺寸)设置为 1.5mm，max size(最大尺寸）设为4mm;完成W上设置后，通过Generate Mesh(产生网格)完 成网格划分；
[0061 ] 9.在流体力学计算模块Fluent的solution(解决）中设置材料参数(血液密度、血 液粘度、血管壁密度），使模型的物理属性接近人体生物学属性，提高仿真的准确度(如图20 所示）。求解控制参数(计算步长、迭代次数、最大循环次数)和边界条件(命名血流入口面， 并赋予速度值，命名出口面后赋予压力值，未命名的血管壁设置为wall), 口静脉血流雷诺 数Re<2000，故仿真流体设置为层流;运算初始化设置为从入口面开始;完成上述参数设定 后，模拟血管壁与血液的流体-固体禪合，计算获得仿真Ξ维血管模型的压力分布和血流分 布。如图21和图22所示。
[0062] 10.在results后处理模块中，对结果进行读取，通过contour(轮廓）操作显示肝- 口静脉模型压力分布图；利用软件自带的calculators选项卡，计算获取虚拟肝静脉自由压 (vFHVP)数值。如图23和图24所示。
[0063] 11.创建一个直径大于或等于被截断血管的圆柱体来模拟阻断球囊，通过布尔运 算对肝右静脉进行阻断，得到模拟肝静脉模入压(vWHVP)的开放的几何模型，并WIGS文件 导出备用（参考步骤F);将IGS文件导入到ANSYS woricbench中，在产生的两个截面分别赋予 Om/s的速度模拟血流停滞，其余材料参数、边界条件和求解控制参数不变，计算得到虚拟肝 静脉模入压(vWHVP)，在results模块中读取虚拟肝静脉模入压(vWHVP)。如图25和图26所 /J、- 〇
[0064] 12.求出所述vWHVP与所述vFHVP的差值，即为虚拟肝静脉压力梯度(vHVPG)。
[0065] 实施例:流体力学仿真计算vHVPG
[0066] 在流体力学仿真模块Fluent中，设置材料参数:血液密度= 1050kg/m3、血液粘度 = 0.00化g/m · S、血管壁密度=1150kg/m3。求解控制参数(计算步长、迭代次数、最大循环 次数)和边界条件(血流速度），模拟血管壁与血液的"流体-固体禪合"，仿真FHVP和WHVP的 血液流线分布(图27的A图和B图）。
[0067] 通过后处理模块分析计算，获得虚拟肝静脉自由压（virtual free hepatic venous pressure, vFHVP)和虚拟肝静脉模入压（virtual wedged hepatic venous pressure,vWHVP)分别为730化和3226化（图27的C图和D图），计算两者差值得到vHVPG = 2496化。
[0068] 通过该病例在经颈静脉肝内口体分流术中有创测得的压力可知，FHVP和WHVP分别 为606化和3139化（图27的E图和F图），计算可得HVPG = 2533化。
[0069] 比较结果显示，本发明对于虚拟肝静脉压力梯度(vHVPG)的测定准确性非常接近 于在临床诊断试验水平，有创HVPG已经纳入国际指南推荐的金标准，而本发明的无创测量 方法能克服前期参考标准不理想且易受全麻开腹干扰的问题，因此，本发明提出的vHVPG无 创检测新技术切实可行，全面评估vHVPG的诊断价值，有望为肝硬化口脉高压的无创诊断提 供新的思路，为改善口脉高压患者的生活质量、减轻家庭和社会的疾病负担起到积极推动 作用。
[0070] 计算流体力学是解决流体计算问题行之有效的途径，其主流方法为有限元分析 法:将一个连续的求解域任意分成适当形状的许多微小单元并于各小单元分片构造插值函 数，根据极值原理(变分或加权余量法)将问题的控制方程转化为所有单元上的有限元方 程，把总体的极值作为各单元极值之和，将局部单元总体合成，形成嵌入了指定边界条件的 代数方程组，求解该方程组即可得到各节点待求的函数值。流体计算的基本方程是化Vier- Stokes方程，其在直角坐标中可写成：
[0071]
[0072] 式中，Δ是拉普拉斯算子，P是流体密度，P是压力，u、v、w是流体t时刻在点（x，y，z) 处的速度分量，Χ、Υ、Ζ是外力的分量，常数μ是动力粘性系数。肝静脉-Π 静脉系统的血管壁 实际变形量较小，血管中血液黏度较低且可视为恒定的黏度(流体的内摩擦剪切力与单位 距离上的两层流体间的相对速度近似成线性关系），血液密度近似为常量，故肝静脉-Π 静 脉系统中的血液可近似为不可压缩的牛顿流体。判断流体运动状态是层流、素流或端流的 无量纲参量是雷诺数(Reynolds number,Re)，其定义Re = pVDAi(P是流体密度、V是流体速 度、D是内直径、μ是流体粘度）。当Re<2000时流体为层流状态，当Re〉3000时为端流。计算可 知肝静脉-Π 静脉系统中的血液做层流运动。
【主权项】
1. 一种虚拟肝静脉压力梯度的测量方法，其特征在于，包括以下步骤： A. 从标本的肘正中静脉注射造影剂，进行CT血管成像(CTA)，获取包括肝静脉期在内的 CTA图层序列，导出图层序列，格式为dicom，层厚1.25mm，图像分辨率512 X 512像素； B. 将所获取的CTA图层序列导入医学影像控制软件MMICS，选择肝静脉期图层序列，设 置图像序列的方位，MMICS软件自动识别图像序列，生成肝静脉期CTA图像序列的冠状位、 矢状位和水平位图像； C. 寻找图像中的肝静脉-门静脉系统（目标），利用MMICS软件的ThreshoIding(阈值算 法)工具，以尽量包含目标的CT值、尽量排除临近目标的周围肝脏等软组织的CT值为原则， 设定阈值范围，以提取目标;利用MIMICS软件的Regiongrowing(区域增长)工具，选定目标， 以提取只与目标在空间结构上有连接的结构;利用MHOCS软件的Calculate 3D from mask (三维建模)工具，选择quality(质量)为medium(中等精度），建立初步的肝静脉-门静脉三 维模型； D. 利用MMICS软件的Crop mask(裁剪蒙板)工具进一步提取目标结构，剔除部分非目 标结构;再利用MMICS软件的Edit masks in 3D(三维编辑蒙版)工具，剔除剩余的非目标 结构，只保留肝静脉-门静脉系统;反复利用MMICS软件的Edit masks in 3D(三维切割)工 具和Edit mask(二维编辑蒙板)工具，进一步对肝静脉-门静脉系统进行选择性填充、剔除 噪点像素，从而重建内腔封闭的肝静脉-门静脉系统实心三维模型； E. 利用MMICS软件的Smoothing(光滑)操作，将三维模型进行表面光滑处理;将光滑处 理后的实心三维模型的几何模型数据选择Ansys area element(Ansys软件面文件）格式 (.inp)导出； F. 在ANSYS经典模式下，导入所述面文件(格式为.inp)，将长度单位统一为国际单位m; 以面为基础建立肝静脉-门静脉系统模型的实心模型体； G. 通过布尔操作对肝静脉-门静脉系统模型的血流入口及出口作垂直切面，得到模拟 肝静脉自由压(vFHVP)的开放的几何模型;完成上述操作后，将文件以后缀为.IGS格式导出 备用（IGS是一种三维数值模型文件格式，ANSYS Workbench模块可读取）； H. 建立ANSYS Workbench有限元计算平台，包括Geometry几何模型模块、Fluent流体计 算模块和Results模块（即CFD-P0ST后处理模块）。通过Geometry模块将IGS文件导入，在 Mesh(划分网格)单元中，划分对象为导入的数值模型，网格划分方法设置为Tetrahedrons (四面体型），在Physics Preference(物理设置）中选CFD(计算流体力学分析），在Solver Preference (求解设置）中选Fluent (使用Fluent求解流场）；考虑到运算精确度以及计算机 运行速度，对划分网格尺寸进行限定，max face size(最大面尺寸)设置为1.5mm，max size (最大尺寸)设为4mm;完成以上设置后，通过Generate Mesh(产生网格)完成网格划分； I. 在流体力学计算模块Fluent的solution(解决方案）中设置材料参数(血液密度、血 液粘度、血管壁密度），使模型的物理属性接近人体生物学属性，提高仿真的准确度;求解控 制参数(计算步长、迭代次数、最大循环次数)和边界条件(命名血流入口面，并赋予速度值， 命名出口面后赋予压力值，未命名的血管壁设置为wal 1 )，门静脉血流雷诺数Re〈2000，故仿 真流体设置为层流;运算初始化设置为从入口面开始;完成上述参数设定后，模拟血管壁与 血液的流体-固体耦合，计算获得仿真三维血管模型的压力分布和血流分布； J. 在results后处理模块中，对结果进行读取，通过contour(轮廓)操作显示肝-门静脉 模型压力分布图；利用软件自带的calculators选项卡，计算获取虚拟肝静脉自由压 (vHWP)数值； K. 创建一个直径大于或等于被截断血管的圆柱体来模拟阻断球囊，通过布尔运算对肝 右静脉进行阻断，得到模拟肝静脉楔入压(vWHVP)的开放的几何模型，并以IGS文件导出备 用(参考步骤F);将IGS文件导入到ANSYS workbench中，在产生的两个截面分别赋予Om/s的 速度模拟血流停滞，其余材料参数、边界条件和求解控制参数不变，计算得到虚拟肝静脉楔 入压(vWHVP)，在results模块中读取虚拟肝静脉楔入压(vWHVP); L. 求出所述vWHVP与所述vFHVP的差值，即为虚拟肝静脉压力梯度(vHVPG)。2.根据权利要求1所述的虚拟肝静脉压力梯度的测量方法，其特征在于:所述步骤I中， 所述材料参数被设置为:血液密度= l〇50kg/m3、血液粘度= 0.005kg/m · s、血管壁密度= 1150kg/m3。
【文档编号】G06T19/00GK105825070SQ201610211335
【公开日】2016年8月3日
【申请日】2016年4月6日
【发明人】祁小龙, 李国新, 黄家乐, 朱炎杰, 赵永昭, 刘燕娜
【申请人】南方医科大学南方医院