[go: up one dir, main page]

CN114677492B - Coronary microcirculation hemodynamics simulation analysis method and device - Google Patents

Coronary microcirculation hemodynamics simulation analysis method and device Download PDF

Info

Publication number
CN114677492B
CN114677492B CN202210385802.3A CN202210385802A CN114677492B CN 114677492 B CN114677492 B CN 114677492B CN 202210385802 A CN202210385802 A CN 202210385802A CN 114677492 B CN114677492 B CN 114677492B
Authority
CN
China
Prior art keywords
coronary artery
coronary
imr
microcirculation
dimensional geometric
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.)
Active
Application number
CN202210385802.3A
Other languages
Chinese (zh)
Other versions
CN114677492A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202210385802.3A priority Critical patent/CN114677492B/en
Publication of CN114677492A publication Critical patent/CN114677492A/en
Application granted granted Critical
Publication of CN114677492B publication Critical patent/CN114677492B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/02028Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/60ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT 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
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/24Fluid dynamics
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/41Medical

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)
  • Biomedical Technology (AREA)
  • Theoretical Computer Science (AREA)
  • Biophysics (AREA)
  • Primary Health Care (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Epidemiology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Cardiology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • Computational Linguistics (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Physiology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明公开了一种冠状动脉微循环血流动力学仿真分析方法及装置,包括以下步骤:1)获取对象的冠状动脉CT血管造影图像;2)基于对象特定的冠状动脉CT血管造影图像,进行冠状动脉三维重建,得到对象相应的冠状动脉三维几何模型;3)对所述的冠状动脉三维几何模型进行网格划分及边界条件设定;4)实测冠状动脉全部或部分分支IMR值,根据对象实测IMR值对上述边界条件进行修正;5)对所述的冠状动脉三维几何模型进行血流动力学仿真分析,计算和可视化所述的冠状动脉三维几何模型的血流速度分布和压力分布,实现对象整个冠状动脉微循环系统全面展示。本发明将IMR值与计算流体力学相结合,降低了有创手术的医疗风险和手术成本。

The present invention discloses a coronary microcirculation hemodynamics simulation analysis method and device, comprising the following steps: 1) obtaining a coronary CT angiography image of an object; 2) based on the coronary CT angiography image specific to the object, performing a three-dimensional reconstruction of the coronary artery to obtain a corresponding three-dimensional geometric model of the coronary artery of the object; 3) performing grid division and boundary condition setting on the three-dimensional geometric model of the coronary artery; 4) measuring the IMR value of all or part of the branches of the coronary artery, and correcting the above boundary conditions according to the measured IMR value of the object; 5) performing hemodynamics simulation analysis on the three-dimensional geometric model of the coronary artery, calculating and visualizing the blood flow velocity distribution and pressure distribution of the three-dimensional geometric model of the coronary artery, and realizing a comprehensive display of the entire coronary microcirculation system of the object. The present invention combines the IMR value with computational fluid dynamics, reducing the medical risks and surgical costs of invasive surgery.

Description

Simulation analysis method and device for coronary artery microcirculation hemodynamics
Technical Field
The invention relates to the field of hemodynamic simulation analysis, in particular to a method and a system for simulated analysis of coronary microcirculation hemodynamic.
Background
Methods for clinically assessing the state of coronary microcirculation are largely divided into non-invasive and invasive types. The noninvasive evaluation is mainly performed by methods such as transthoracic Doppler echocardiography, single photon emission computer tomography, myocardial magnetic resonance imaging and the like, but the noninvasive evaluation method is greatly influenced by the level of operators, the resolution of equipment and the technical means, and is difficult to accurately evaluate the microcirculation state.
The pressure guide wire is utilized to invasively measure the microcirculation resistance index (Index ofMicrocirculatory Resistance, IMR) so as to quantitatively evaluate the coronary microcirculation resistance, is not easily influenced by the change of the hemodynamics, has the advantages of good repeatability, strong specificity, high reliability and the like, and becomes a gold standard for evaluating the coronary microcirculation diseases clinically.
The pressure guide wire mainly used in clinic at present is provided with a baroreceptor and a temperature-sensitive receptor at a position 3cm away from the tail end, and the guide wire body is provided with another temperature-sensitive receptor. The guidewire was delivered to the distal epicardial coronary artery and baroreceptor readings (Pd) were recorded. By directing the catheter to inject normal saline at room temperature, the time (T) between two temperature-sensitive receptors of the liquid is recorded according to a thermal dilution curve, and the product of Pd and T is the IMR. The IMR measurement operation is complex, only the IMR corresponding to a certain blood vessel is generally measured clinically, and even if the IMR of a plurality of blood vessels is measured, the whole coronary microcirculation system of the subject is difficult to comprehensively evaluate.
Disclosure of Invention
In view of the above, the invention provides a simulation analysis method and device for coronary artery microcirculation hemodynamics, aiming at the problem of difficult overall evaluation of the state of the whole coronary artery microcirculation system in the prior art.
In a first aspect, a method for simulating and analyzing coronary artery microcirculation hemodynamics includes the steps of:
1) Acquiring a coronary artery CT angiography image of a subject;
2) Based on the CT angiography image of the specific coronary artery of the object, carrying out three-dimensional reconstruction of the coronary artery to obtain a three-dimensional geometrical model of the coronary artery corresponding to the object;
3) Performing grid division and boundary condition setting on the three-dimensional geometric model of the coronary artery;
4) The actual measurement of all or part of the branch IMR values of the coronary arteries, and the boundary conditions are corrected according to the actual measurement of the IMR values of the object;
5) And carrying out hemodynamic simulation analysis on the three-dimensional geometric model of the coronary artery, and calculating and visualizing the blood flow velocity distribution and the pressure distribution of the three-dimensional geometric model of the coronary artery to realize comprehensive display of the whole coronary artery microcirculation system of the subject.
Compared with the prior art, the method has the following advantages:
(1) The invention combines the IMR value with computational fluid dynamics, and realizes the hemodynamic simulation of a specific coronary artery microcirculation system of an object through boundary condition correction;
(2) The coronary artery hemodynamic simulation method adopted by the invention can intuitively represent the flow velocity and pressure distribution of the coronary artery, and simultaneously non-invasively measure the hemodynamic parameters such as the flow velocity, the pressure, the coronary artery blood flow reserve, the wall shear stress and the like, thereby reducing the medical risk and the operation cost of invasive operation;
(3) The invention realizes the hemodynamic simulation of the object-specific coronary artery microcirculation system, provides help for explaining the influence of the coronary artery microcirculation change on the hemodynamic parameters such as blood flow, pressure, FFR and the like, and is beneficial to researching the contribution of coronary artery microcirculation disturbance on myocardial ischemia from the aspects of hydrodynamics and coronary artery physiology.
In a refinement, in step 2), the three-dimensional reconstruction of the coronary artery comprises a left coronary artery three-dimensional geometry starting from the left coronary artery inlet and extending to the branching distal end of the left anterior descending branch (LAD) and the left circumflex branch (LCX), and a right coronary artery three-dimensional geometry starting from the right coronary artery inlet and extending to the branching distal end of the Right Coronary Artery (RCA). And identifying and analyzing epicardial stenosis and plaque, and determining the real lumen boundary of the coronary artery to obtain a three-dimensional geometrical model of the coronary artery.
As an improvement, between steps 2) and 3), i.e. after obtaining the object-specific three-dimensional geometric model of the coronary arteries, before meshing, the object-specific three-dimensional geometric model of the coronary arteries is preprocessed by pruning the distal branches with blurry or too small diameters (diameter <1 mm) in view of the resolution of the coronary artery CT angiography images, correcting local error topologies (spikes, hollows, etc.), global smoothing the geometric model, and finally cutting the vessel model perpendicular to the local vessel centerline to generate the inlet and outlet of the coronary artery geometric model as the locations for the inflow and outflow of blood.
As an improvement, in step 3), the meshing and boundary condition setting respectively comprise performing fine meshing, discretizing the model into unstructured tetrahedral units, setting the inlet boundary condition as mean arterial pressure, setting the outlet boundary condition as far-end real microcirculation resistance value, and applying no-slip and solid wall conditions to the vessel wall.
As an improvement, in step 4), the boundary condition correction includes:
a. For simulation of the normal coronary artery microcirculation state of a subject, namely IMR is less than or equal to 25 mmHg.s, the real microcirculation resistance value in the normal microcirculation state is distributed to all branch outlets according to the corrected Mory law, and the corrected Mory law formula is as follows:
wherein, Q 1,Q2,D1,D2 is the flow and diameter of two distal branches at the bifurcation, respectively;
b. Simulation of the subject coronary artery microcirculation disturbance state, i.e. IMR >30 mmHg.s, assigning TMR CMD value under the microcirculation disturbance state to lesion branch outlet, calculation of TMR CMD is based on assumption of IMR and TMR linear relation, and is obtained by multiplying normal TMR value by correction factor lambda, i.e.:
TMRCMD=TMR·λ#
wherein TMR is the real microcirculation resistance value in the normal microcirculation state, lambda is a correction factor, TMR CMD is the microcirculation resistance value in the coronary artery microcirculation disturbance state;
The calculating method of the correction factor lambda is the ratio of the measured IMR value of the object to the normal IMR cut-off value, namely:
Wherein IMR S is the measured IMR value of the subject, IMR n is the normal IMR cut-off value, and in the method, IMR n is defined as 25 mmHg.s.
In a refinement, in step 5), the hemodynamic analysis includes assuming the blood is an incompressible Newtonian fluid, and the computational fluid dynamics (Computational Fluid Dynamics, CFD) modeling is based on an incompressible continuity equation and a Navier-Stokes equation:
And
Wherein, Is the three-dimensional velocity vector of the fluid, ρ is the fluid density, p is the pressure, μ is the blood flow viscosity, and t is the time.
The coronary artery microcirculation hemodynamic simulation analysis device comprises a processor, a memory, a touch display and a power supply, wherein the memory is used for storing instructions, the processor is used for calling the instructions in the memory and executing any one of the first aspect or any one of the possible implementation manners of the first aspect, and the touch display is used for displaying coronary artery CT angiography images, three-dimensional geometric models of the coronary arteries and hemodynamic simulation analysis results and performing interactive operation.
As an improvement, the simulation analysis device for the coronary artery microcirculation hemodynamics further comprises a deep learning module, wherein the deep learning module is used for meshing of the three-dimensional geometric model of the coronary artery, and the rapid and automatic meshing is realized by establishing the deep learning model.
As an improvement, the coronary artery microcirculation hemodynamic simulation analysis device further comprises a 5G communication module, wherein the 5G communication module is used for acquiring a target coronary artery CTA image from the CT imaging equipment in real time and transmitting an analysis result back to the electronic medical record system after the hemodynamic simulation analysis is finished.
Drawings
FIG. 1 is a flow chart of a method of simulated analysis of coronary microcirculation hemodynamics of the present invention;
FIG. 2 is a coronary CT angiography image in an embodiment of the invention;
FIG. 3 is a schematic view of reconstructing a three-dimensional geometric model of a coronary artery in an embodiment of the present invention;
FIG. 4 is a schematic diagram of a three-dimensional geometric preprocessing model of a coronary artery in an embodiment of the invention;
FIG. 5 is a flow velocity profile of a hemodynamic simulation of a microcirculation disturbance in an embodiment of the present invention;
fig. 6 is a pressure distribution of a hemodynamic simulation of a microcirculation disturbance in an embodiment of the present invention.
Detailed Description
The invention will be further described with reference to the drawings and specific examples, but the invention is not limited to these examples only. The invention is intended to cover any alternatives, modifications, equivalents, and variations that fall within the spirit and scope of the invention. In the following description of preferred embodiments of the invention, specific details are set forth in order to provide a thorough understanding of the invention, and the invention will be fully understood to those skilled in the art without such details.
As shown in fig. 1, the hemodynamic simulation method of the coronary artery microcirculation disturbance based on the microcirculation resistance index of the present invention comprises the following steps:
(1) Acquiring a CTA image and an IMR value of an object:
the method comprises the steps of obtaining a plurality of layers of CTA scanning data samples of an object, wherein a 128-layer computed tomography scanner is used for CTA image acquisition, the equipment resolution is 192 x 0.5mm, the scanning matrix is 512 x 512, and a CTA image is stored in a DICOM mode. Measuring the distal coronary pressure P d by using a pressure guide wire with the diameter of 0.014 inch, injecting normal saline into the coronary artery pellets in the maximum hyperemia state, obtaining the average transit time T mn according to a thermal dilution curve, and calculating to obtain an IMR value;
(2) The three-dimensional model reconstruction of the coronary artery of the object CTA image is realized:
And loading CTA data into the chemicals software to reconstruct the three-dimensional geometrical model of the coronary artery of the object. The geometric model includes the ascending aorta, left coronary artery trunk (LMCA), left anterior descending branch (LAD), left circumflex branch (LCX), right Coronary Artery (RCA), and other observable tiny branches. And at the editable three-view interface, using a coronary artery semiautomatic segmentation tool, selecting the position of an aorta at the coronary artery opening as a starting point, selecting the position of the farthest branches of the left and right observable coronary arteries as an ending point, automatically creating a blood vessel path along the trend of the coronary artery blood vessel, and segmenting the lumen of the coronary artery. Based on the difference of CTA image pixel gray values of different tissues, calcified plaques are manually identified and segmented, and finally the model is repaired by operations such as matrix reduction, contour element reduction and the like, so that a real three-dimensional geometrical model of the coronary artery is obtained, as shown in figure 2.
(3) Model preprocessing to obtain a high-quality three-dimensional geometric model of the coronary artery:
The three-dimensional geometrical model of the coronary artery of the object is imported into Geomagic Studio software, the fuzzy or smaller (diameter <1 mm) distal branches are removed, the error geometrical structure such as holes, nails, self-intersection points, triangular faces and the like are modified, and global smoothing processing is carried out. Then, the inlet and outlet of the model are determined, and the blood vessel is cut perpendicular to the central line of the local blood vessel, so that the approximately round coronary inlet and outlet with flat surfaces are obtained.
(4) Grid division and boundary condition setting of a coronary artery model:
The processed coronary artery model was grid-partitioned using ANSYS software, which required to complete the dependency check on grid density, using unstructured tetrahedral grids with a global maximum cell length defined as 3mm and the maximum cell lengths of the model inlet, outlet and vessel walls defined as 0.1mm, and 0.2mm, respectively. Appropriate boundary conditions were assigned at the inlet and outlet of the model, the coronary artery inlet set to an average aortic pressure of 110mmHg and the coronary artery outlet set to the calculated TMR CMD1 value. The no-slip condition applies to the vessel wall.
In one possible embodiment, the clinical measurement of IMR is 78mmhg·s, and the correction factor λ 1 is determined based on the following equation:
Wherein IMR S is a subject-specific clinical IMR measurement, IMR n is a normal IMR cutoff, in this method defined as 25mmhg·s, calculated as λ 1 =3.12.
The true microcirculation resistance value TMR in the normal congestion state was set to 6mmhg·s·ml -1, and TMR CMD1 of this example was calculated as:
TMRCMD1=λ·TMR=3.12×6mmHg·s·mL-1=18.72mmHg·s·mL-1
(5) Numerical method of coronary artery model and CFD simulation:
blood flow simulation was performed using ANSYS software, with blood modeled as an incompressible newtonian fluid, employing an incompressible continuity equation and the Navier-Stokes equation:
wherein, Is the three-dimensional velocity vector of the fluid, ρ is the fluid density, set to 1050kg/m 3, p is the pressure, μ is the blood flow viscosity, set to 0.0035Pa, and t is the time.
And solving a continuity equation and a Navier-Stokes equation through CFD to obtain the flow velocity distribution and the pressure distribution of the coronary microcirculation disturbance model, wherein the flow velocity distribution and the pressure distribution are respectively shown in figures 3 and 4.
The embodiment of the invention provides a simulation analysis device for coronary artery microcirculation hemodynamics, which comprises a processor, a memory, a touch display and a power supply, wherein the memory is used for storing instructions, the processor is used for calling the instructions in the memory and executing any one possible implementation mode of the first aspect or the first aspect, and the touch display is used for displaying coronary artery CT angiography images, a three-dimensional geometric model of the coronary artery and a simulation analysis result of the hemodynamics and performing interactive operation. The coronary artery microcirculation hemodynamic simulation analysis device is communicated with the image archiving and communication system PACS through a network, and can acquire the CTA image of the coronary artery of the subject in real time.
In one possible design, the coronary microcirculation hemodynamic simulation analysis device further comprises a deep learning module, wherein the deep learning module is used for meshing of the three-dimensional geometric model of the coronary artery, and rapid and automatic meshing is realized by establishing the deep learning model.
In one possible design, the coronary artery microcirculation hemodynamic simulation analysis device further comprises a 5G communication module, which is used for acquiring the CTA image of the coronary artery of the subject from the CT imaging device in real time and transmitting the analysis result back to the electronic medical record system after the hemodynamic simulation analysis is finished.
The foregoing is illustrative of the preferred embodiments of the present invention, and is not to be construed as limiting the claims. The present invention is not limited to the above embodiments, and the specific structure thereof is allowed to vary. In general, all changes that come within the meaning and range of equivalency of the claims are to be embraced within their scope.

Claims (3)

1.一种冠状动脉微循环血流动力学仿真分析方法,包括以下步骤:1. A coronary microcirculation hemodynamics simulation analysis method, comprising the following steps: 1)获取对象的冠状动脉CT血管造影图像;1) obtaining a coronary CT angiography image of a subject; 2)基于对象特定的冠状动脉CT血管造影图像,进行冠状动脉三维重建,得到对象相应的冠状动脉三维几何模型;2) based on the subject's specific coronary artery CT angiography image, perform coronary artery 3D reconstruction to obtain a corresponding coronary artery 3D geometric model of the subject; 3)对所述的冠状动脉三维几何模型进行网格划分及边界条件设定;3) Meshing the coronary artery three-dimensional geometric model and setting boundary conditions; 4)实测冠状动脉全部或部分分支IMR值,根据对象实测IMR值对上述边界条件进行修正;4) Measure the IMR values of all or some of the coronary artery branches and modify the above boundary conditions according to the measured IMR values of the subjects; 5)对所述的冠状动脉三维几何模型进行血流动力学仿真分析,计算和可视化所述的冠状动脉三维几何模型的血流速度分布和压力分布,实现对象整个冠状动脉微循环系统全面展示;5) performing hemodynamic simulation analysis on the three-dimensional geometric model of the coronary artery, calculating and visualizing the blood flow velocity distribution and pressure distribution of the three-dimensional geometric model of the coronary artery, and realizing a comprehensive display of the entire coronary microcirculatory system of the subject; 在步骤2)中,所述冠状动脉三维重建包括以下过程:左冠状动脉三维几何结构从左冠状动脉入口开始,延伸至左前降支和左回旋支的分支远端;右冠状动脉三维几何结构从右冠状动脉入口开始,延伸至右冠状动脉的分支远端;识别和分析心外膜狭窄及斑块,确定冠状动脉的真实管腔边界,得到冠状动脉三维几何模型;In step 2), the three-dimensional reconstruction of the coronary arteries includes the following processes: the three-dimensional geometric structure of the left coronary artery starts from the entrance of the left coronary artery and extends to the distal ends of the branches of the left anterior descending artery and the left circumflex artery; the three-dimensional geometric structure of the right coronary artery starts from the entrance of the right coronary artery and extends to the distal ends of the branches of the right coronary artery; the epicardial stenosis and plaque are identified and analyzed, the true lumen boundary of the coronary artery is determined, and the three-dimensional geometric model of the coronary artery is obtained; 在步骤3)中,所述网格划分包括:进行精细网格划分,将模型离散化为非结构四面体单元;所述边界条件设定包括:入口边界条件设置为平均主动脉压,出口边界条件设置为远端真实微循环阻力值,血管壁应用无滑动和实心壁条件;In step 3), the meshing includes: performing fine meshing to discretize the model into unstructured tetrahedral units; the boundary condition setting includes: setting the inlet boundary condition to the mean aortic pressure, setting the outlet boundary condition to the distal real microcirculatory resistance value, and applying no-slip and solid wall conditions to the vascular wall; 在步骤4)中,边界条件修正包括:In step 4), the boundary condition correction includes: a、对于对象正常冠状动脉微循环状态,即IMR≤25mmHg·s,的模拟,根据校正的默里定律将正常微循环状态下的真实微循环阻力值分配到所有分支出口,校正的默里定律公式为:a. For the simulation of the normal coronary microcirculation state of the subject, that is, IMR ≤ 25 mmHg·s, the actual microcirculation resistance value under the normal microcirculation state is distributed to all branch outlets according to the corrected Murray's law. The corrected Murray's law formula is: 其中,Q1,Q2,D1,D2分别为分叉处两个远端分支的流量和直径;Among them, Q 1 , Q 2 , D 1 , and D 2 are the flow and diameter of the two distal branches at the bifurcation, respectively; b、对于对象冠状动脉微循环障碍状态,即IMR>30mmHg·s,的模拟:将微循环障碍状态下的TMRCMD值分配到病变分支出口;TMRCMD的计算基于IMR和TMR线性关系的假设,由正常的TMR值乘以修正因子λ得到,即:b. For the simulation of the coronary microcirculatory disorder state of the subject, that is, IMR>30mmHg·s: the TMR CMD value under the microcirculatory disorder state is assigned to the outlet of the lesion branch; the calculation of TMR CMD is based on the assumption of a linear relationship between IMR and TMR, and is obtained by multiplying the normal TMR value by the correction factor λ, that is: TMRCMD=TMR·λTMR CMD = TMR·λ 其中,TMR为正常微循环状态下的真实微循环阻力值,λ为修正因子,TMRCMD为冠状动脉微循环障碍状态下的微循环阻力值;Among them, TMR is the true microcirculatory resistance value under normal microcirculatory conditions, λ is the correction factor, and TMR CMD is the microcirculatory resistance value under coronary microcirculatory disorder conditions; 修正因子λ的计算方法为对象实测IMR值与正常IMR截止值的比值,即:The correction factor λ is calculated as the ratio of the measured IMR value of the subject to the normal IMR cutoff value, that is: 其中,IMRS为对象实测IMR值,IMRn为正常IMR截止值,本方法中,定义IMRn为25mmHg·s;Among them, IMR S is the measured IMR value of the subject, IMR n is the normal IMR cutoff value, and in this method, IMR n is defined as 25 mmHg·s; 在步骤5)中,血流动力学分析包括:将血液假定为不可压缩的牛顿流体,计算流体动力学建模基于不可压缩的连续性方程和Navier-Stokes方程:In step 5), the hemodynamic analysis includes: assuming that the blood is an incompressible Newtonian fluid, and the computational fluid dynamics modeling is based on the incompressible continuity equation and the Navier-Stokes equation: and 其中,是流体三维速度矢量,ρ是流体密度,p是压力,μ是血流粘性,t是时间。in, is the three-dimensional velocity vector of the fluid, ρ is the fluid density, p is the pressure, μ is the blood flow viscosity, and t is the time. 2.如权利要求1所述的冠状动脉微循环血流动力学仿真分析方法,其特征在于:获得所述的冠状动脉三维几何模型后,网格划分之前,对所述的冠状动脉三维几何模型进行预处理:对所述的冠状动脉三维几何模型局部错误拓扑结构进行校正,再进行全局平滑处理,最后,垂直于局部血管中心线环切割血管模型,生成冠状动脉几何模型的入口和出口,作为血流进入和离开的位置。2. The coronary microcirculation hemodynamics simulation analysis method as described in claim 1 is characterized in that: after obtaining the coronary artery three-dimensional geometric model and before meshing, the coronary artery three-dimensional geometric model is preprocessed: the local erroneous topological structure of the coronary artery three-dimensional geometric model is corrected, and then global smoothing is performed. Finally, the vascular model is cut perpendicular to the local vascular centerline to generate the entrance and exit of the coronary artery geometric model as the location where blood flow enters and leaves. 3.一种冠状动脉微循环血流动力学仿真分析装置,包括:处理器、存储器、触控显示器和电源;其中,存储器用于存储指令;处理器用于调用存储器中的指令,执行上述权利要求1-2任一项所述的方法;触控显示器用于冠状动脉CT血管造影图像、冠状动脉三维几何模型和血流动力学仿真分析结果显示,以及交互操作。3. A coronary microcirculation hemodynamic simulation and analysis device, comprising: a processor, a memory, a touch display and a power supply; wherein the memory is used to store instructions; the processor is used to call the instructions in the memory to execute the method described in any one of claims 1-2; the touch display is used to display coronary CT angiography images, coronary three-dimensional geometric models and hemodynamic simulation analysis results, as well as interactive operations.
CN202210385802.3A 2022-04-13 2022-04-13 Coronary microcirculation hemodynamics simulation analysis method and device Active CN114677492B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210385802.3A CN114677492B (en) 2022-04-13 2022-04-13 Coronary microcirculation hemodynamics simulation analysis method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210385802.3A CN114677492B (en) 2022-04-13 2022-04-13 Coronary microcirculation hemodynamics simulation analysis method and device

Publications (2)

Publication Number Publication Date
CN114677492A CN114677492A (en) 2022-06-28
CN114677492B true CN114677492B (en) 2024-12-06

Family

ID=82077507

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210385802.3A Active CN114677492B (en) 2022-04-13 2022-04-13 Coronary microcirculation hemodynamics simulation analysis method and device

Country Status (1)

Country Link
CN (1) CN114677492B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116167293B (en) * 2022-12-30 2025-10-03 深圳北芯生命科技股份有限公司 Calculation method of fractional flow reserve of main vascular trunks based on computational fluid dynamics
CN116453697B (en) * 2022-12-30 2023-12-19 徐州医科大学 Coronary artery stenosis hemodynamic simulation method and system based on FFR fitting
CN117438055B (en) * 2023-10-08 2025-07-04 中国人民解放军陆军军医大学第一附属医院 A non-invasive MRR calculation method based on computational fluid dynamics
CN117174244A (en) * 2023-11-03 2023-12-05 北京唯迈医疗设备有限公司 Balloon type selection method and system based on intraoperative image simulation

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108550189A (en) * 2018-05-03 2018-09-18 苏州润迈德医疗科技有限公司 Microcirculation drag index computational methods based on contrastographic picture and fluid mechanic model
CN113015497A (en) * 2019-10-15 2021-06-22 Ai医疗股份公司 Method and device for simulating blood flow of blood vessel inherent to object

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8315812B2 (en) * 2010-08-12 2012-11-20 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
WO2019071249A1 (en) * 2017-10-06 2019-04-11 Emory University Methods and systems for determining hemodynamic information for one or more arterial segments
CN115482212A (en) * 2021-12-17 2022-12-16 深圳市阅影科技有限公司 Method and device for determining functional indexes of blood vessels based on blood flow model

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108550189A (en) * 2018-05-03 2018-09-18 苏州润迈德医疗科技有限公司 Microcirculation drag index computational methods based on contrastographic picture and fluid mechanic model
CN113015497A (en) * 2019-10-15 2021-06-22 Ai医疗股份公司 Method and device for simulating blood flow of blood vessel inherent to object

Also Published As

Publication number Publication date
CN114677492A (en) 2022-06-28

Similar Documents

Publication Publication Date Title
JP7483079B2 (en) Method and system for determining treatment by modifying a patient-specific geometric model - Patents.com
JP7304390B2 (en) Method and system for patient-specific blood flow modeling
CN114677492B (en) Coronary microcirculation hemodynamics simulation analysis method and device
US10134129B2 (en) Method and system for hemodynamic computation in coronary arteries
US11589924B2 (en) Non-invasive assessment and therapy guidance for coronary artery disease in diffuse and tandem lesions
US10803995B2 (en) Method and system for non-invasive functional assessment of coronary artery stenosis using flow computations in diseased and hypothetical normal anatomical models
US10872698B2 (en) Method and system for enhancing medical image-based blood flow computations using physiological measurements
JP6378779B2 (en) System and method for determining blood flow characteristics using a flow ratio
JP6215469B2 (en) Method and system for modeling blood flow using boundary conditions for optimized diagnostic capabilities
JP5868554B2 (en) System and method for estimating blood flow characteristics from blood vessel shape and physiology
CN108109698B (en) System for calculating fractional flow reserve and method for setting boundary conditions
KR102130254B1 (en) Method and apparatus for simulating blood flow of a subject-specific blood vessel
WO2019025270A9 (en) Non-invasive assessment and therapy guidance for coronary artery disease in diffuse and tandem lesions
CN112384138B (en) Method, device, system and storage medium for acquiring blood flow of great artery of heart table
CN115985504A (en) Method for reconstructing coronary artery model based on coronary angiography and related equipment

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant