CN105825070A - 一种虚拟肝静脉压力梯度的测量方法 - Google Patents

一种虚拟肝静脉压力梯度的测量方法 Download PDF

Info

Publication number
CN105825070A
CN105825070A CN201610211335.7A CN201610211335A CN105825070A CN 105825070 A CN105825070 A CN 105825070A CN 201610211335 A CN201610211335 A CN 201610211335A CN 105825070 A CN105825070 A CN 105825070A
Authority
CN
China
Prior art keywords
vein
hepatic vein
blood
portal
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610211335.7A
Other languages
English (en)
Inventor
祁小龙
李国新
黄家乐
朱炎杰
赵永昭
刘燕娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southern Hospital Southern Medical University
Original Assignee
Southern Hospital Southern Medical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southern Hospital Southern Medical University filed Critical Southern Hospital Southern Medical University
Priority to CN201610211335.7A priority Critical patent/CN105825070A/zh
Publication of CN105825070A publication Critical patent/CN105825070A/zh
Priority to PCT/CN2016/107208 priority patent/WO2017173834A1/zh
Priority to US16/090,289 priority patent/US10827935B2/en
Priority to ZA2018/06596A priority patent/ZA201806596B/en
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/021—Measuring pressure in heart or blood vessels
    • A61B5/0215—Measuring pressure in heart or blood vessels by means inserted into the body
    • A61B5/02152—Measuring pressure in heart or blood vessels by means inserted into the body specially adapted for venous pressure
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/481—Diagnostic techniques involving the use of contrast agents
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00—Manipulating 3D models or images for computer graphics
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00—Manipulating 3D models or images for computer graphics
    • G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/41—Medical
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Cardiology (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Vascular Medicine (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Physiology (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Architecture (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明属于早期无创性诊断领域，涉及一种虚拟肝静脉压力梯度的测量方法。本发明所述的一种虚拟肝静脉压力梯度的测量方法包括：肝静脉‑门静脉系统三维建模；有限元网格划分数学模型；以及流体力学仿真计算虚拟肝静脉压力梯度(vHVPG)。本发明优化并完善肝静脉‑门静脉系统三维建模、有限元网格划分和流体力学仿真计算，构建并验证一种更具诊断优势的vHVPG检测新技术，为门脉高压患者的早期无创诊断提供一种安全无创、准确量化的新途径。

Description

一种虚拟肝静脉压力梯度的测量方法
技术领域
本发明属于早期无创性诊断领域，涉及一种虚拟肝静脉压力梯度的测量方法。
背景技术
在全球范围内，肝脏疾病(病毒性感染、酒精性肝病、非酒精性脂肪肝及相关的肝硬化、肝细胞肝癌)是导致死亡的主要原因之一，仅在中国肝脏疾病就威胁了近3亿人的健康和生命，严重加剧了全球疾病的负担。以乙型肝炎病毒感染引起的慢性肝病为例，据世界卫生组织报道，全球约20亿人曾感染乙肝病毒，其中2.4亿人为慢性乙肝病毒感染者，每年约有65万人死于乙肝病毒感染所致的肝硬化、肝细胞肝癌等。在全球肝硬化患者中，由乙肝病毒感染引起的为30％。而在我国，由于感染乙肝病毒而导致的肝硬化比例高达60％。因此，肝硬化的高发病率和较差的临床结局已造成了国内外严重的社会公共健康问题。
门脉高压是肝硬化失代偿期的重要表现之一，由于其本身并无症状，临床上门脉高压的发现和诊断常由于出现晚期严重的并发症，如：急性静脉曲张破裂出血、顽固性腹水、肝性脑病、门脉高压性胃肠病、肝肾综合症、肝肺综合症和继发感染等。因此，门脉高压及其并发症严重影响了肝硬化患者的生活质量和长期预后，对于门脉高压(尤其是早期无症状性门脉高压)的诊断和监测是肝硬化等终末期肝病治疗链中最重要的环节之一。
在最新的2015版BavenoVI共识(门脉高压的风险分层和个体化管理)中，有创性肝静脉压力梯度(hepaticvenouspressuregradient,HVPG)的测量继续被推荐为诊断临床显著性门脉高压的金标准，即HVPG≥10mmHg可确诊(deFranchisR#*,BavenoVIFaculty.Expandingconsensusinportalhypertension:ReportoftheBavenoVIConsensusWorkshop:Stratifyingriskandindividualizingcareforportalhypertension.JHepatol.2015；63(3):743-752.)。该侵袭性方法是通过颈内静脉穿刺置管，依次经过颈内静脉、上腔静脉、右心房、下腔静脉进入肝静脉，分别测量肝静脉自由压(freehepaticvenouspressure,FHVP)和肝静脉楔入压(wedgedhepaticvenouspressure,WHVP)，计算两者差值得到HVPG[BloomS#,KempW,LubelJ*.Portalhypertension:pathophysiology,diagnosisandmanagement.InternMedJ.2015；45(1):16-26.]。此外，HVPG对于治疗应答、并发症风险和长期预后也可提供重要信息，其数值的变化被推荐为研究预后的替代指标。此外，HVPG测量还被鼓励用于研究门脉高压创新诊断和治疗方法的临床试验中。然而，测量HVPG也存在一些问题，首先其检测手段为有创性，在门脉高压早期尚无严重并发症时很难被患者所接受；其次，测量过程需由经过专业化培训的介入医生在肝静脉造影的辅助下进行，这不但增加了受试者的射线暴露和造影剂过敏可能，还存在有一定的技术操作风险；加之其高昂的诊断费用，HVPG测定在我国受到很大限制，目前仅有部分三级甲等医院开展。除此之外，门脉高压的有创性诊断技术还包括：超声引导下门脉穿刺测压和开腹手术中直接门脉测压[deFranchisR#*,Dell'EraA.Invasiveandnoninvasivemethodstodiagnoseportalhypertensionandesophagealvarices.ClinLiverDis.2014；18(2):293-302.]。然而，直接测量门脉压力的风险更大、对操作者的技术水平要求也更高，在临床实践中很难被医生和患者所采纳，目前仍多用于动物实验性研究。
门脉压力的无创性评估是目前该领域研究的重要内容，主要包括以下三个方面：(1)检测肝内阻力：肝内阻力增加是导致门脉高压形成的重要因素。内源性缩血管物质的释放和肝组织结构的紊乱可促使肝内阻力的显著升高。有研究提示，血清内皮素-1表达水平与HVPG存在相关性，但该结论尚缺乏大样本量临床试验的支持。另一方面，肝内机械结构的改变也会导致肝内阻力增加，临床上多采用肝瞬时弹性成像(FibroScan)技术，通过对肝脏硬度值的测定来评估肝硬化程度，该方法操作简便、重复性好。但也有报道指出，FibroScan的诊断准确性受到肥胖、肋间隙狭窄、炎症等因素干扰，对门脉高压并发症(如高风险食管静脉曲张)的预测价值不高，需联合其他无创指标进行综合判断。(2)测定门脉循环血量：门脉高压的直接因素为循环血量的积聚升高超过了机体的代偿，当肝硬化伴发重度门脉高压时，心输出量也会明显升高且与HVPG存在高度相关性。有研究利用彩色多普勒超声(colorDopplerultrasound,CDUS)对肝硬化患者的门脉直径和血流速度进行测定，在此基础上计算门脉系统的循环血量。然而，CDUS测量技术的主观性较强且易受到肥胖等因素干扰，尚无准确的临床证据表明经血流速度和门脉直径计算的门脉压力与HVPG具有较高的诊断一致性。此外，CT血管成像(CTangiography,CTA)和磁共振静脉显像也可以评估门脉系统的循环血量。然而，其测定原理仍是基于门脉直径和侧枝循环来间接评价门脉压力。(3)仿真计算门脉压力：虚拟仿真技术是计算机图形学、医学图像处理、软件工程计算等多学科交叉融合的结晶。以计算机仿真辅助医学影像获得的虚拟人体数据为疾病的介入诊断和治疗提供了全新的思路。目前已有多项国际前瞻性临床试验表明，基于CTA和流体力学计算获得的血流储备分数对于冠状动脉的功能性狭窄具有高度诊断价值，相关仿真计算软件也已获得美国食品药品监督管理局的上市批准。
门脉压力的各种测定技术目前仍受到各自干扰因素的影响，或者有创风险大、操作难度高，或者干扰因素多、数值变异大。此外，国内外该领域研究的主要问题还表现在：评估范围局限于门脉主干，未考虑肝静脉-门静脉系统的整体影响；仅考虑到门脉血液动力学因素，忽略了肝硬化患者的血液流变学改变，如血液黏度的变化和红细胞变形能力的降低都会促使门脉阻力的升高。
目前国内外相关报道中，门脉高压的诊断方法，特别是早期无创性诊断技术，仍有待进一步研究完善。
发明内容
本发明的目的在于提供一种虚拟肝静脉压力梯度的测量方法，以此为基础可以构建一种更具诊断优势的虚拟肝静脉压力梯度检测新技术，为门脉高压患者的早期无创诊断提供一种新的途径。
本发明所述的一种虚拟肝静脉压力梯度的测量方法，包括以下步骤：
A.从标本的肘正中静脉注射造影剂，进行CT血管成像(CTA)，获取包括肝静脉期在内的CTA图层序列，导出图层序列，格式为dicom，层厚1.25mm，图像分辨率512×512像素；
B.将所获取的CTA图层序列导入医学影像控制软件MIMICS，选择肝静脉期图层序列，设置图像序列的方位，MIMICS软件自动识别图像序列，生成肝静脉期CTA图像序列的冠状位、矢状位和水平位图像；
C.寻找图像中的肝静脉-门静脉系统(目标)，利用MIMICS软件的Thresholding(阈值算法)工具，以尽量包含目标的CT值、尽量排除临近目标的周围肝脏等软组织的CT值为原则，设定阈值范围，以提取目标；利用MIMICS软件的Regiongrowing(区域增长)工具，选定目标，以提取只与目标在空间结构上有连接的结构；利用MIMICS软件的Calculate3Dfrommask(三维建模)工具，选择quality(质量)为medium(中等精度)，建立初步的肝静脉-门静脉三维模型；
D.利用MIMICS软件的Cropmask(裁剪蒙板)工具进一步提取目标结构，剔除部分非目标结构；再利用MIMICS软件的Editmasksin3D(三维编辑蒙版)工具，剔除剩余的非目标结构，只保留肝静脉-门静脉系统；反复利用MIMICS软件的Editmasksin3D(三维切割)工具和Editmask(二维编辑蒙板)工具，进一步对肝静脉-门静脉系统进行选择性填充、剔除噪点像素，从而重建内腔封闭的肝静脉-门静脉系统实心三维模型；
E.利用MIMICS软件的Smoothing(光滑)操作，将三维模型进行表面光滑处理；将光滑处理后的实心三维模型的几何模型数据选择Ansysareaelement(Ansys软件面文件)格式(.inp)导出；
F.在ANSYS经典模式下，导入所述面文件(格式为.inp)，将长度单位统一为国际单位m；以面为基础建立肝静脉-门静脉系统模型的实心模型体；
G.通过布尔操作对肝静脉-门静脉系统模型的血流入口及出口作垂直切面，得到模拟肝静脉自由压(vFHVP)的开放的几何模型；完成上述操作后，将文件以后缀为.IGS格式导出备用(IGS是一种三维数值模型文件格式，ANSYSWorkbench模块可读取)；
H.建立ANSYSWorkbench有限元计算平台，包括Geometry几何模型模块、Fluent流体计算模块和Results模块(即CFD-POST后处理模块)。通过Geometry模块将IGS文件导入，在Mesh(划分网格)单元中，划分对象为导入的数值模型，网格划分方法设置为Tetrahedrons(四面体型)，在PhysicsPreference(物理设置)中选CFD(计算流体力学分析)，在SolverPreference(求解设置)中选Fluent(使用Fluent求解流场)；考虑到运算精确度以及计算机运行速度，对划分网格尺寸进行限定，maxfacesize(最大面尺寸)设置为1.5mm，maxsize(最大尺寸)设为4mm；完成以上设置后，通过GenerateMesh(产生网格)完成网格划分；
I.在流体力学计算模块Fluent的solution(解决方案)中设置材料参数(血液密度、血液粘度、血管壁密度)，使模型的物理属性接近人体生物学属性，提高仿真的准确度；求解控制参数(计算步长、迭代次数、最大循环次数)和边界条件(命名血流入口面，并赋予速度值，命名出口面后赋予压力值，未命名的血管壁设置为wall)，门静脉血流雷诺数Rε<2000，故仿真流体设置为层流；运算初始化设置为从入口面开始；完成上述参数设定后，模拟血管壁与血液的流体-固体耦合，计算获得仿真三维血管模型的压力分布和血流分布；
J.在results后处理模块中，对结果进行读取，通过contour(轮廓)操作显示肝-门静脉模型压力分布图；利用软件自带的calculators选项卡，计算获取虚拟肝静脉自由压(vFHVP)数值；
K.创建一个直径大于或等于被截断血管的圆柱体来模拟阻断球囊，通过布尔运算对肝右静脉进行阻断，得到模拟肝静脉楔入压(vWHVP)的开放的几何模型，并以IGS文件导出备用(参考步骤F)；将IGS文件导入到ANSYSworkbench中，在产生的两个截面分别赋予0m/s的速度模拟血流停滞，其余材料参数、边界条件和求解控制参数不变，计算得到虚拟肝静脉楔入压(vWHVP)，在results模块中读取虚拟肝静脉楔入压(vWHVP)；
L.求出所述vWHVP与所述vFHVP的差值，即为虚拟肝静脉压力梯度(vHVPG)。
根据本发明所述的虚拟肝静脉压力梯度的测量方法的进一步特征，所述步骤I中，所述材料参数被设置为：血液密度＝1050kg/m3、血液粘度＝0.005kg/m·s、血管壁密度＝1150kg/m3。
如前所述，门脉高压作为肝硬化失代偿期的重要特征，其临床明确诊断往往滞后于晚期严重的并发症发生。然而，目前国内外的门脉高压诊断技术或者有创风险大、操作难度高，或者干扰因素多、数值变异大。本发明首先提出了虚拟肝静脉压力梯度(virtualhepaticvenouspressuregradient,vHVPG)的概念，并在其方法学构建中取得重要数据，有利于探索并建立一种安全无创、准确量化的门脉高压诊断技术。本发明优化并完善肝静脉-门静脉系统三维建模、有限元网格划分和流体力学仿真计算，构建并验证一种更具诊断优势的vHVPG检测新技术，为门脉高压患者的早期无创诊断提供一种新的途径。
附图说明
图1为从CT血管成像(CTA)导出包括肝静脉期在内的CTA图层序列。
图2为选择导入的肝静脉期CTA图层序列，层厚1.25mm。
图3为设置图层序列的方位。
图4为自动生成的冠状位、矢状位和水平位图像。
图5为Thresholding(阈值算法)工具设定阈值范围。
图6为Regiongrowing(区域增长)工具选定目标。
图7为Calculate3Dfrommask(三维建模)，选择中等精度。
图8为建立初步的肝-门静脉系统。
图9为Cropmask(裁剪蒙板)工具进一步提取目标结构，剔除部分非目标结构。
图10为Cropmask(裁剪蒙板)工具操作后保留的三维模型。
图11为Editmasksin3D(三维编辑蒙版)工具，剔除剩余的非目标结构。
图12为Editmasksin3D(三维编辑蒙版)工具操作后保留的肝静脉-门静脉系统。
图13为Editmask(二维编辑蒙板)工具在冠状位、矢状位和水平位编辑蒙板。
图14为反复利用MIMICS的Editmasksin3D(三维编辑蒙版)工具和Editmask(二维编辑蒙板)工具操作后的肝静脉-门静脉实体模型。
图15为对肝静脉-门静脉系统进行Smoothing(光滑)操作。
图16为将光滑处理后的模型选择.inp格式导出。
图17为以面为基础建立肝静脉-门静脉系统模型的实心模型体。
图18为模拟肝静脉自由压(FHVP)的开放的几何模型。
图19为ANSYSWorkbench有限元分析计算平台。
图20为solution(解决)单元界面。
图21为设定求解控制参数及迭代数目。
图22为求解获得仿真三维血管模型的压力分布和血流分布。
图23为虚拟肝静脉自由压vFHVP模拟及压力分布图。
图24为读取虚拟肝静脉自由压vFHVP数值。
图25为虚拟肝静脉楔入压vWHVP数值模型。
图26为虚拟肝静脉楔入压vWHVP结果及压力分布图。
图27为流体力学仿真计算vHVPG与有创测得的压力的比较结果图。
具体实施方式
下面结合附图对本发明进行详细说明。
本发明所述的一种虚拟肝静脉压力梯度的测量方法包括：肝静脉-门静脉系统三维建模；有限元网格划分数学模型；以及流体力学仿真计算vHVPG。
1.从标本的肘正中静脉注射造影剂，进行CT血管成像(CTA)，获取包括肝静脉期在内的CTA图层序列，导出图层序列，格式为dicom，层厚1.25mm，图像分辨率512×512像素。如图1所示。
2.将所获取的CTA图层序列导入医学影像控制软件MIMICS，选择肝静脉期图层序列(如图2所示)，设置图像序列的方位(如图3所示)，MIMICS软件自动识别图像序列，生成肝静脉期CTA图像序列的冠状位、矢状位和水平位图像(如图4所示)。
C.寻找图像中的肝静脉-门静脉系统(目标)，利用MIMICS软件的Thresholding(阈值算法)工具，以尽量包含目标的CT值、尽量排除临近目标的周围肝脏等软组织的CT值为原则，设定阈值范围，以提取目标(如图5所示)；利用MIMICS软件的Regiongrowing(区域增长)工具，选定目标，以提取只与目标在空间结构上有连接的结构(如图6所示)；利用MIMICS软件的Calculate3Dfrommask(三维建模)工具，选择quality(质量)为medium(中等精度)，建立初步的肝静脉-门静脉三维模型(如图7和图8所示)。
4.利用MIMICS软件的Cropmask(裁剪蒙板)工具进一步提取目标结构，剔除部分非目标结构(如图9和图10所示)；再利用MIMICS软件的Editmasksin3D(三维编辑蒙版)工具，剔除剩余的非目标结构，只保留肝静脉-门静脉系统(如图11和图12所示)；反复利用MIMICS软件的Editmasksin3D(三维切割)工具和Editmask(二维编辑蒙板)工具，进一步对肝静脉-门静脉系统进行选择性填充、剔除噪点像素，从而重建内腔封闭的肝静脉-门静脉系统实心三维模型。如图13至图15所示。
5.利用MIMICS软件的Smoothing(光滑)操作，将三维模型进行表面光滑处理；将光滑处理后的实心三维模型的几何模型数据选择Ansysareaelement(Ansys软件面文件)格式(.inp)导出。如图16所示。
6.在ANSYS经典模式下，导入所述面文件(格式为.inp)，将长度单位统一为国际单位m；以面为基础建立肝静脉-门静脉系统模型的实心模型体(如图17所示)。
7.通过布尔操作对肝静脉-门静脉系统模型的血流入口及出口作垂直切面，得到模拟肝静脉自由压(vFHVP)的开放的几何模型(如图18所示)；完成上述操作后，将文件以后缀为.IGS格式导出备用(IGS是一种三维数值模型文件格式，ANSYSWorkbench模块可读取)。
8.建立ANSYSWorkbench有限元计算平台(如图19所示)，包括Geometry几何模型模块、Fluent流体计算模块和Results模块(即CFD-POST后处理模块)。通过Geometry模块将IGS文件导入，在Mesh(划分网格)单元中，划分对象为导入的数值模型，网格划分方法设置为Tetrahedrons(四面体型)，在PhysicsPreference(物理设置)中选CFD(计算流体力学分析)，在SolverPreference(求解设置)中选Fluent(使用Fluent求解流场)；考虑到运算精确度以及计算机运行速度，对划分网格尺寸进行限定，maxfacesize(最大面尺寸)设置为1.5mm，maxsize(最大尺寸)设为4mm；完成以上设置后，通过GenerateMesh(产生网格)完成网格划分；
9.在流体力学计算模块Fluent的solution(解决)中设置材料参数(血液密度、血液粘度、血管壁密度)，使模型的物理属性接近人体生物学属性，提高仿真的准确度(如图20所示)。求解控制参数(计算步长、迭代次数、最大循环次数)和边界条件(命名血流入口面，并赋予速度值，命名出口面后赋予压力值，未命名的血管壁设置为wall)，门静脉血流雷诺数Rε<2000，故仿真流体设置为层流；运算初始化设置为从入口面开始；完成上述参数设定后，模拟血管壁与血液的流体-固体耦合，计算获得仿真三维血管模型的压力分布和血流分布。如图21和图22所示。
10.在results后处理模块中，对结果进行读取，通过contour(轮廓)操作显示肝-门静脉模型压力分布图；利用软件自带的calculators选项卡，计算获取虚拟肝静脉自由压(vFHVP)数值。如图23和图24所示。
11.创建一个直径大于或等于被截断血管的圆柱体来模拟阻断球囊，通过布尔运算对肝右静脉进行阻断，得到模拟肝静脉楔入压(vWHVP)的开放的几何模型，并以IGS文件导出备用(参考步骤F)；将IGS文件导入到ANSYSworkbench中，在产生的两个截面分别赋予0m/s的速度模拟血流停滞，其余材料参数、边界条件和求解控制参数不变，计算得到虚拟肝静脉楔入压(vWHVP)，在results模块中读取虚拟肝静脉楔入压(vWHVP)。如图25和图26所示。
12.求出所述vWHVP与所述vFHVP的差值，即为虚拟肝静脉压力梯度(vHVPG)。
实施例：流体力学仿真计算vHVPG
在流体力学仿真模块Fluent中，设置材料参数：血液密度＝1050kg/m3、血液粘度＝0.005kg/m·s、血管壁密度＝1150kg/m3。求解控制参数(计算步长、迭代次数、最大循环次数)和边界条件(血流速度)，模拟血管壁与血液的“流体-固体耦合”，仿真FHVP和WHVP的血液流线分布(图27的A图和B图)。
通过后处理模块分析计算，获得虚拟肝静脉自由压(virtualfreehepaticvenouspressure,vFHVP)和虚拟肝静脉楔入压(virtualwedgedhepaticvenouspressure,vWHVP)分别为730Pa和3226Pa(图27的C图和D图)，计算两者差值得到vHVPG＝2496Pa。
通过该病例在经颈静脉肝内门体分流术中有创测得的压力可知，FHVP和WHVP分别为606Pa和3139Pa(图27的E图和F图)，计算可得HVPG＝2533Pa。
比较结果显示，本发明对于虚拟肝静脉压力梯度(vHVPG)的测定准确性非常接近于在临床诊断试验水平，有创HVPG已经纳入国际指南推荐的金标准，而本发明的无创测量方法能克服前期参考标准不理想且易受全麻开腹干扰的问题，因此，本发明提出的vHVPG无创检测新技术切实可行，全面评估vHVPG的诊断价值，有望为肝硬化门脉高压的无创诊断提供新的思路，为改善门脉高压患者的生活质量、减轻家庭和社会的疾病负担起到积极推动作用。
计算流体力学是解决流体计算问题行之有效的途径，其主流方法为有限元分析法：将一个连续的求解域任意分成适当形状的许多微小单元并于各小单元分片构造插值函数，根据极值原理(变分或加权余量法)将问题的控制方程转化为所有单元上的有限元方程，把总体的极值作为各单元极值之和，将局部单元总体合成，形成嵌入了指定边界条件的代数方程组，求解该方程组即可得到各节点待求的函数值。流体计算的基本方程是Navier-Stokes方程，其在直角坐标中可写成：
&rho; d u d t = - &part; p &part; x + &rho; X + &mu; &Delta; u &rho; d v d t = - &part; p &part; y + &rho; Y + &mu; &Delta; v &rho; d w d t = - &part; p &part; z + &rho; Z + &mu; &Delta; w ,
式中，Δ是拉普拉斯算子，ρ是流体密度，p是压力，u、v、w是流体t时刻在点(x，y，z)处的速度分量，X、Y、Z是外力的分量，常数μ是动力粘性系数。肝静脉-门静脉系统的血管壁实际变形量较小，血管中血液黏度较低且可视为恒定的黏度(流体的内摩擦剪切力与单位距离上的两层流体间的相对速度近似成线性关系)，血液密度近似为常量，故肝静脉-门静脉系统中的血液可近似为不可压缩的牛顿流体。判断流体运动状态是层流、紊流或湍流的无量纲参量是雷诺数(Reynoldsnumber,Re)，其定义Re＝ρVD/μ(ρ是流体密度、V是流体速度、D是内直径、μ是流体粘度)。当Re<2000时流体为层流状态，当Re>3000时为湍流。计算可知肝静脉-门静脉系统中的血液做层流运动。

Claims (2)

1.一种虚拟肝静脉压力梯度的测量方法，其特征在于，包括以下步骤：
A.从标本的肘正中静脉注射造影剂，进行CT血管成像(CTA)，获取包括肝静脉期在内的CTA图层序列，导出图层序列，格式为dicom，层厚1.25mm，图像分辨率512×512像素；
B.将所获取的CTA图层序列导入医学影像控制软件MIMICS，选择肝静脉期图层序列，设置图像序列的方位，MIMICS软件自动识别图像序列，生成肝静脉期CTA图像序列的冠状位、矢状位和水平位图像；
C.寻找图像中的肝静脉-门静脉系统(目标)，利用MIMICS软件的Thresholding(阈值算法)工具，以尽量包含目标的CT值、尽量排除临近目标的周围肝脏等软组织的CT值为原则，设定阈值范围，以提取目标；利用MIMICS软件的Regiongrowing(区域增长)工具，选定目标，以提取只与目标在空间结构上有连接的结构；利用MIMICS软件的Calculate3Dfrommask(三维建模)工具，选择quality(质量)为medium(中等精度)，建立初步的肝静脉-门静脉三维模型；
D.利用MIMICS软件的Cropmask(裁剪蒙板)工具进一步提取目标结构，剔除部分非目标结构；再利用MIMICS软件的Editmasksin3D(三维编辑蒙版)工具，剔除剩余的非目标结构，只保留肝静脉-门静脉系统；反复利用MIMICS软件的Editmasksin3D(三维切割)工具和Editmask(二维编辑蒙板)工具，进一步对肝静脉-门静脉系统进行选择性填充、剔除噪点像素，从而重建内腔封闭的肝静脉-门静脉系统实心三维模型；
E.利用MIMICS软件的Smoothing(光滑)操作，将三维模型进行表面光滑处理；将光滑处理后的实心三维模型的几何模型数据选择Ansysareaelement(Ansys软件面文件)格式(.inp)导出；
F.在ANSYS经典模式下，导入所述面文件(格式为.inp)，将长度单位统一为国际单位m；以面为基础建立肝静脉-门静脉系统模型的实心模型体；
G.通过布尔操作对肝静脉-门静脉系统模型的血流入口及出口作垂直切面，得到模拟肝静脉自由压(vFHVP)的开放的几何模型；完成上述操作后，将文件以后缀为.IGS格式导出备用(IGS是一种三维数值模型文件格式，ANSYSWorkbench模块可读取)；
H.建立ANSYSWorkbench有限元计算平台，包括Geometry几何模型模块、Fluent流体计算模块和Results模块(即CFD-POST后处理模块)。通过Geometry模块将IGS文件导入，在Mesh(划分网格)单元中，划分对象为导入的数值模型，网格划分方法设置为Tetrahedrons(四面体型)，在PhysicsPreference(物理设置)中选CFD(计算流体力学分析)，在SolverPreference(求解设置)中选Fluent(使用Fluent求解流场)；考虑到运算精确度以及计算机运行速度，对划分网格尺寸进行限定，maxfacesize(最大面尺寸)设置为1.5mm，maxsize(最大尺寸)设为4mm；完成以上设置后，通过GenerateMesh(产生网格)完成网格划分；
I.在流体力学计算模块Fluent的solution(解决方案)中设置材料参数(血液密度、血液粘度、血管壁密度)，使模型的物理属性接近人体生物学属性，提高仿真的准确度；求解控制参数(计算步长、迭代次数、最大循环次数)和边界条件(命名血流入口面，并赋予速度值，命名出口面后赋予压力值，未命名的血管壁设置为wall)，门静脉血流雷诺数Rε<2000，故仿真流体设置为层流；运算初始化设置为从入口面开始；完成上述参数设定后，模拟血管壁与血液的流体-固体耦合，计算获得仿真三维血管模型的压力分布和血流分布；
J.在results后处理模块中，对结果进行读取，通过contour(轮廓)操作显示肝-门静脉模型压力分布图；利用软件自带的calculators选项卡，计算获取虚拟肝静脉自由压(vFHVP)数值；
K.创建一个直径大于或等于被截断血管的圆柱体来模拟阻断球囊，通过布尔运算对肝右静脉进行阻断，得到模拟肝静脉楔入压(vWHVP)的开放的几何模型，并以IGS文件导出备用(参考步骤F)；将IGS文件导入到ANSYSworkbench中，在产生的两个截面分别赋予0m/s的速度模拟血流停滞，其余材料参数、边界条件和求解控制参数不变，计算得到虚拟肝静脉楔入压(vWHVP)，在results模块中读取虚拟肝静脉楔入压(vWHVP)；
L.求出所述vWHVP与所述vFHVP的差值，即为虚拟肝静脉压力梯度(vHVPG)。
2.根据权利要求1所述的虚拟肝静脉压力梯度的测量方法，其特征在于：所述步骤I中，所述材料参数被设置为：血液密度＝1050kg/m3、血液粘度＝0.005kg/m·s、血管壁密度＝1150kg/m3。
CN201610211335.7A 2016-04-06 2016-04-06 一种虚拟肝静脉压力梯度的测量方法 Pending CN105825070A (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201610211335.7A CN105825070A (zh) 2016-04-06 2016-04-06 一种虚拟肝静脉压力梯度的测量方法
PCT/CN2016/107208 WO2017173834A1 (zh) 2016-04-06 2016-11-25 一种虚拟肝静脉压力梯度的测量方法
US16/090,289 US10827935B2 (en) 2016-04-06 2016-11-25 Method of determining virtual hepatic venous pressure gradient
ZA2018/06596A ZA201806596B (en) 2016-04-06 2018-10-04 Method of determining virtual hepatic venous pressure gradient

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610211335.7A CN105825070A (zh) 2016-04-06 2016-04-06 一种虚拟肝静脉压力梯度的测量方法

Publications (1)

Publication Number Publication Date
CN105825070A true CN105825070A (zh) 2016-08-03

Family

ID=56525779

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610211335.7A Pending CN105825070A (zh) 2016-04-06 2016-04-06 一种虚拟肝静脉压力梯度的测量方法

Country Status (4)

Country Link
US (1) US10827935B2 (zh)
CN (1) CN105825070A (zh)
WO (1) WO2017173834A1 (zh)
ZA (1) ZA201806596B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017173834A1 (zh) * 2016-04-06 2017-10-12 祁小龙 一种虚拟肝静脉压力梯度的测量方法
CN107480675A (zh) * 2017-07-28 2017-12-15 祁小龙 一种基于放射组学的肝静脉压力梯度计算模型的构建方法
CN107945878A (zh) * 2017-11-16 2018-04-20 祁小龙 一种基于放射组学的肝静脉压力梯度计算模型的构建方法
CN107928655A (zh) * 2017-12-14 2018-04-20 江苏理工学院 一种基于人体动脉3d血管网络模型的中心动脉血压测定方法
CN109758136A (zh) * 2019-02-28 2019-05-17 四川大学华西医院 基于门脉血管特征的肝静脉压力梯度测量方法
CN109924997A (zh) * 2018-05-02 2019-06-25 上海大学 一种基于医学图像的无创性门静脉压力梯度测量方法
CN110074809A (zh) * 2019-06-17 2019-08-02 祁小龙 Ct图像的肝静脉压力梯度分类方法以及计算机设备
CN110111429A (zh) * 2019-03-16 2019-08-09 哈尔滨理工大学 一种检测单像素血管的方法
CN110216872A (zh) * 2019-06-17 2019-09-10 薛启煌 一种基于cta数据的肝脏三维重建方法及其3d打印模型
CN110270004A (zh) * 2019-07-25 2019-09-24 苏州恒瑞宏远医疗科技有限公司 Hvpg球囊导管及其制作方法
CN110769741A (zh) * 2017-06-30 2020-02-07 佐治亚州立大学研究基金会 用于检测肝纤维化的非侵入性方法
CN112906323A (zh) * 2021-04-08 2021-06-04 兰州大学 一种虚拟肝静脉压力梯度测量方法
CN113033121A (zh) * 2021-04-14 2021-06-25 兰州大学 门静脉高压经颈静脉肝内门体分流术支架直径的选择方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110705161B (zh) * 2019-09-27 2022-11-25 湘潭大学 考虑摩擦热分布的机械密封不对中状态热仿真分析方法
US20210278936A1 (en) * 2020-03-09 2021-09-09 Biosense Webster (Israel) Ltd. Electrophysiological user interface
CN111462217A (zh) * 2020-03-13 2020-07-28 司法鉴定科学研究院 一种颅内血肿体积的测量方法及系统
CN114121227A (zh) * 2021-02-24 2022-03-01 首都医科大学附属北京佑安医院 一种肝硬化患者上腹部增强ct影像组学特征处理方法及其应用
CN117473906B (zh) * 2023-12-26 2024-04-19 浙江远算科技有限公司 一种基于流体动力学仿真的风电机舱后处理方法和介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103796579A (zh) * 2011-09-01 2014-05-14 微创医学科技有限公司 检测门静脉和/或肝静脉压力的方法和门静脉高压监控系统
CN103976720A (zh) * 2013-12-17 2014-08-13 上海交通大学医学院附属仁济医院 血管内压力的无创测量方法及其装置
CN104107039A (zh) * 2013-04-17 2014-10-22 上海市同济医院 一种无创性门静脉血流动力学参数测定方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BRPI0819434A2 (pt) * 2007-11-08 2019-09-24 Univ Utah Res Found método de tratamento ou prevenção de doenças hepáticas e uso de antagonistas á angiogênese em condições de proliferação venosa anormal
ES2368404B1 (es) * 2008-12-01 2012-09-24 Fundacio Clinic Per A La Recerca Biomedica Cateter percutáneo y método para medir presión en vasos sanguíneos.
US9846765B2 (en) * 2012-11-09 2017-12-19 Siemens Healthcare Gmbh System and method for patient specific modeling of liver tumor ablation
WO2014190170A1 (en) * 2013-05-22 2014-11-27 The Regents Of The University Of California Biomarker composite test for hepatic vein pressure gradient and cirrhosis treatment
CN105825070A (zh) 2016-04-06 2016-08-03 南方医科大学南方医院 一种虚拟肝静脉压力梯度的测量方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103796579A (zh) * 2011-09-01 2014-05-14 微创医学科技有限公司 检测门静脉和/或肝静脉压力的方法和门静脉高压监控系统
CN104107039A (zh) * 2013-04-17 2014-10-22 上海市同济医院 一种无创性门静脉血流动力学参数测定方法
CN103976720A (zh) * 2013-12-17 2014-08-13 上海交通大学医学院附属仁济医院 血管内压力的无创测量方法及其装置

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
XIAOLONG QI: "A novel noninvasive assessment of hepatic venous", 《ARCH MED SCI》 *
XIAOLONG QI: "Virtual portal pressure gradient from anatomic CT angiography", 《JOURNAL OF HEPATOLOPY》 *
朱炎杰: "门静脉压力的无创性评估现状及进展", 《临床肝胆病杂志》 *
祁小龙: "《门静脉压力测定方法及进展》", 《中华实用诊断与治疗杂志》 *
祁小龙: "三维模型有限元分析和流体力学计算在门静脉压力无创检测中的应用", 《中华医学会第三次全国肝纤维化，肝硬化学术会议》 *
祁小龙: "基于门静脉三维模型的新型无创性门静脉压力测定", 《2012中国消化系疾病学术大会论文汇编》 *
马婉蓉: "肝静脉压力梯度研究进展", 《中华实用诊断与治疗杂志》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10827935B2 (en) 2016-04-06 2020-11-10 Xiaolong Qi Method of determining virtual hepatic venous pressure gradient
WO2017173834A1 (zh) * 2016-04-06 2017-10-12 祁小龙 一种虚拟肝静脉压力梯度的测量方法
CN110769741A (zh) * 2017-06-30 2020-02-07 佐治亚州立大学研究基金会 用于检测肝纤维化的非侵入性方法
US11576609B2 (en) 2017-06-30 2023-02-14 Georgia State University Research Foundation, Inc. Noninvasive methods for detecting liver fibrosis
CN107480675A (zh) * 2017-07-28 2017-12-15 祁小龙 一种基于放射组学的肝静脉压力梯度计算模型的构建方法
CN107480675B (zh) * 2017-07-28 2020-12-04 祁小龙 一种基于放射组学的肝静脉压力梯度计算模型的构建方法
CN107945878A (zh) * 2017-11-16 2018-04-20 祁小龙 一种基于放射组学的肝静脉压力梯度计算模型的构建方法
WO2019095516A1 (zh) * 2017-11-16 2019-05-23 祁小龙 一种基于放射组学的肝静脉压力梯度计算模型的构建方法
CN107928655A (zh) * 2017-12-14 2018-04-20 江苏理工学院 一种基于人体动脉3d血管网络模型的中心动脉血压测定方法
CN107928655B (zh) * 2017-12-14 2023-07-11 江苏理工学院 一种基于人体动脉3d血管网络模型的中心动脉血压测定方法
CN109924997A (zh) * 2018-05-02 2019-06-25 上海大学 一种基于医学图像的无创性门静脉压力梯度测量方法
CN109758136A (zh) * 2019-02-28 2019-05-17 四川大学华西医院 基于门脉血管特征的肝静脉压力梯度测量方法
CN110111429A (zh) * 2019-03-16 2019-08-09 哈尔滨理工大学 一种检测单像素血管的方法
CN110111429B (zh) * 2019-03-16 2022-11-18 哈尔滨理工大学 一种检测单像素血管的方法
CN110216872A (zh) * 2019-06-17 2019-09-10 薛启煌 一种基于cta数据的肝脏三维重建方法及其3d打印模型
CN110074809A (zh) * 2019-06-17 2019-08-02 祁小龙 Ct图像的肝静脉压力梯度分类方法以及计算机设备
CN110074809B (zh) * 2019-06-17 2023-06-20 祁小龙 Ct图像的肝静脉压力梯度分类方法以及计算机设备
CN110270004A (zh) * 2019-07-25 2019-09-24 苏州恒瑞宏远医疗科技有限公司 Hvpg球囊导管及其制作方法
CN112906323A (zh) * 2021-04-08 2021-06-04 兰州大学 一种虚拟肝静脉压力梯度测量方法
CN113033121B (zh) * 2021-04-14 2023-03-21 兰州大学 门静脉高压经颈静脉肝内门体分流术支架直径的选择方法
CN113033121A (zh) * 2021-04-14 2021-06-25 兰州大学 门静脉高压经颈静脉肝内门体分流术支架直径的选择方法

Also Published As

Publication number Publication date
ZA201806596B (en) 2019-06-26
US20190117094A1 (en) 2019-04-25
WO2017173834A1 (zh) 2017-10-12
US10827935B2 (en) 2020-11-10

Similar Documents

Publication Publication Date Title
CN105825070A (zh) 一种虚拟肝静脉压力梯度的测量方法
Xiao et al. Multi-scale computational model of three-dimensional hemodynamics within a deformable full-body arterial network
CN107491636B (zh) 一种基于计算流体力学的脑血管储备力仿真系统和方法
Xu et al. Assessment of boundary conditions for CFD simulation in human carotid artery
Birchall et al. Analysis of haemodynamic disturbance in the atherosclerotic carotid artery using computational fluid dynamics
Pedrizzetti et al. On estimating intraventricular hemodynamic forces from endocardial dynamics: a comparative study with 4D flow MRI
WO2022166281A1 (zh) 血流动力学指标数据的处理方法和系统
Dueñas-Pamplona et al. A comprehensive comparison of various patient-specific CFD models of the left atrium for atrial fibrillation patients
Wang et al. Aortic arch morphogenesis and flow modeling in the chick embryo
CN107945878B (zh) 一种基于放射组学的肝静脉压力梯度计算模型的构建方法
Hodneland et al. A new framework for assessing subject-specific whole brain circulation and perfusion using MRI-based measurements and a multi-scale continuous flow model
CN103976720B (zh) 利用仿真技术建立血管模型的方法
CN107480675B (zh) 一种基于放射组学的肝静脉压力梯度计算模型的构建方法
Lai et al. Fluid mechanics of blood flow in human fetal left ventricles based on patient-specific 4D ultrasound scans
Wang et al. Spatial distribution of wall shear stress in common carotid artery by color Doppler flow imaging
Bartolo et al. Numerical predictions of shear stress and cyclic stretch in pulmonary hypertension due to left heart failure
Sampath et al. Optimized time-resolved echo particle image velocimetry–particle tracking velocimetry measurements elucidate blood flow in patients with left ventricular thrombus
Moroseos et al. Comparison of Simpson's method and three-dimensional reconstruction for measurement of right ventricular volume in patients with complete or corrected transposition of the great arteries
Müller et al. Sensitivity of flow simulations in a cerebral aneurysm
Antonova et al. Numerical analysis of 3D blood flow and common carotid artery hemodynamics in the carotid artery bifurcation with stenosis
Long Ko et al. Pulsatile hemodynamics in patient-specific thoracic aortic dissection models constructed from computed tomography angiography
Calò et al. Impact of wall displacements on the large-scale flow coherence in ascending aorta
Chen et al. Non-invasive assessment of intracranial wall shear stress using high-resolution magnetic resonance imaging in combination with computational fluid dynamics technique
de Vecchi et al. A novel methodology for personalized simulations of ventricular hemodynamics from noninvasive imaging data
Canè et al. Mixed impact of torsion on LV hemodynamics: A CFD study based on the Chimera technique

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160803

WD01 Invention patent application deemed withdrawn after publication