Computation hydrodynamics come simulate calculate blood flow reserve fraction system and method
The application be the applying date on December 28th, 2016, Application No. 201611234903.1 and it is entitled " should
The divisional application of the Chinese invention application of the system and method that calculating blood flow reserve fraction is simulated with Fluid Mechanics Computation ".
Technical field
The present invention relates to Human physiology and it is hemodynamic calculate analysis system and method, more particularly to
A kind of computation hydrodynamics (CFD) come simulate calculate blood flow reserve fraction system and for vascular tree calculating fluid force
It learns (CFD) model and corresponding entrance boundary condition and the method for export boundary condition is set.
Background technology
Coronary heart disease Etiological is the coronary stenosis caused by artery sclerosis.Coronary stenosis can cause coronary blood flow and other
Hemodynamic great change causes myocardium opposite and absolute ischemic.Coronary angiography and intravascular ultrasound are regarded as examining
" goldstandard " of disconnected coronary heart disease, but they can only carry out stenosis imaging evaluation, and it is narrow on earth to Distal blood flow production
It is unknown much influences (functional evaluation) have been given birth to.It proposes to calculate coronary blood by pressure measurement from nineteen ninety-five Pijls et al.
The New Set of stream-blood flow reserve fraction (Fractional Flow Reserve, FFR) (Pijls NH, Van Gelder B,
Van der Voort P,Peels K,Bracke FA,Bonnier HJ,El Gamal MI.Fractional flow
reserve a useful index to evaluate the influence of an epicardial coronary
stenosis on myocardial blood flow.Circulation.1995Dec 1;92(11):3183-93 is by drawing
It is for reference with being hereby incorporated by), by long-term basis and clinical research, FFR has become to coronary stenosis Evaluation of Functional
" goldstandard " and the indispensable guiding instrument of coronary revascularization.
FFR is as invasive measurement (i.e. intrusive mood measurement), inevitably with certain risk and may be to machine
Body causes a degree of wound.With the development of science and technology and research level, especially medical image imaging and reconstruction technique, meter
The fields such as fluid operator mechanics, high-performance calculation, numerical simulation calculation FFR (virtual FFR) just become a fast-developing forward position
Direction.Compared to invasive FFR, virtual FFR is non-invasive measurement.It is compared according to clinical test, virtual FFR is careful to individual patient
Each blood vessel rank scale on, from accuracy, sensibility, specificity, positive predictive value, negative predictive value to cardiac muscle lack
The diagnosis of blood has very high ability, more aobvious superior in the comparison with the method by imaging evaluation.
Because the superiority of virtual the noninvasive of FFR, quantitative analysis and diagnosis capability, more and more scientific researches and commercial undertaking pair
The technology carries out basis and clinical research.Since this is an emerging direction, various relevant researchs and improved method, mode are equal
Have and dabble, but totally can be attributed to two general orientation:1. the application of medical image and Model Reconstruction 2. are used as Fluid Mechanics Computation
(CFD) key factor of accuracy:Meet the research and processing of the boundary condition of Human physiology feature.Mainstream at present
Computing system include the FFR based on coronary artery CT computer tomography angiographyCT, made based on rotation blood vessel
The vFFR of shadow, the FFR based on Quantitative Coronary Angiography and thrombolysis in myocardial infarctionQCADeng.
As a cutting edge technology, virtual FFR equally faces various technological challenges and problem.Wherein, it is this to calculate the time
Can technology as a key of clinical expansion and new " goldstandard ".For using accuracy and model complexity to be leading
Method calculates a case time cost and was even up to for more than 100 a hours;Method for quickly calculating, optimal time can be with
It foreshortens within 10 minutes, but experienced analyst is required to be operated for this, it is necessary to the help of invasive measurement and to complexity
The processing capacity of vessel tree model be also restricted.Even if equilibrium is carried out between speed and accuracy, using virtual FFR
Carry out simulation calculate take estimation still surpass 7 it is small when, this causes the clinical expansion of this method very big limitation.In addition, accurate
Property aspect, the perimeter strip for being how rationally to obtain/set CFD model in the virtual FFR another significant challenges calculated using CFD
Part.Boundary condition it is suitable whether, determine to the full extent simulation calculate accuracy.
At present, the blood vessel computing system in nearly all mainstream forward position is both needed to assume when setting the boundary condition of CFD model
Or statistics obtains the flow distribution ratio of human vas tree, this is also existing computational methods for individual patient, even refine
To a bottleneck of the physiological characteristic accuracy of each blood vessel of patient itself.Based on the set perimeter strip of this hypothesis or statistics
Part cannot embody the specificity of individual, affect the accuracy that simulation calculates.In addition, existing blood vessel computing system is in order to automatic
And the fluid parameter of the vessel branch of the clinical concern of prediction exactly, it usually needs (such as coronal to relevant entire vascular tree
The entire vascular tree of artery) CFD modelings and calculating are carried out, it takes long, and consumption calculations resource is big, hinders its entrance therefore and faces
Bed is promoted.
The system for calculating FFR and the stream for vascular tree can be quickly and automatically simulated therefore, it is necessary to a kind of via CFD
Mechanics model sets corresponding entrance boundary condition and the method for export boundary condition, which can be directed to patient
Body, it is independent (independent of analyst experience without invasive measurement intervention), it is automatic, quickly and accurately draw to such as
The fluid parameter of vascular tree coronarius everywhere, and its Evaluation of Functional is obtained accordingly, including FFR and to suspected lesion position
It puts, influence of the lesion to myocardial ischemia, complicated hemadostewnosis situation is (such as single branch, multivessel disease, same blood vessel many places disease
Become, moderate stenosis lesion, continuous lesion, diffusivity lesion, in-stent restenosis etc.) etc. assessment etc..Preferably the system and
Method can also be directed to the inclination of the vessel branch progress computing resource of clinical concern, so as in the limited premise of computing resource
Under, the calculating speed of entire vascular tree and the accuracy of the calculating of the clinical vessel branch paid close attention to can be taken into account.
The content of the invention
For above the technical issues of recognizing, the present invention proposes a kind of computation hydrodynamics to simulate calculating blood
Stream lays in the system of fraction and is the corresponding entrance boundary condition of CFD model setting and export boundary condition of vascular tree
Method, the system and method do not do any hypothesis or system in the setting of boundary condition to the specific flow distribution ratio of vascular tree
Meter, but drawn by calculating (such as being calculated according to the medical image calculating of individual or the CFD of vascular tree) actual prediction, from
And ensure that boundary condition is individual specificity, which thereby enhance the accuracy of result of calculation.Also, system and method is another
One important application is actually to predict the flow that each blood vessel of the clinical concern under complex network exports exactly, such as brain
Blood vessel, and overall computation time is taken into account, overall computation time can be controlled in the applicable scope of clinic.Also, the system
It can be completely noninvasive with method and automatically realizes, analyst can be easy to use without rich experiences, without invasive doctor
The support of image is treated, alleviate the burden of analyst while decreases the pain of patient.
According to the first aspect of the invention, provide a kind of computation hydrodynamics and calculate blood flow reserve fraction to simulate
System, including:
Vessel tree model generation module, for the geometry for obtaining medical image, performing the vascular tree for splitting and rebuilding individual
Model;
Grid generation module is calculated, grid is calculated for being generated for the geometrical model of the vascular tree, it is described so as to establish
The CFD model of vascular tree, including three-dimensional (3D) CFD model;
Boundary condition setup module, the CFD model for the vascular tree to be established set corresponding entrance boundary
Condition and export boundary condition;
Attribute setup module, for setting the physical attribute of blood and flow equation;
Solver, for based on the entrance boundary condition and export boundary condition, set physical attribute and flowing
Equation, the CFD model of the vascular tree to being established solve, to obtain the fluid parameter of vascular tree everywhere;And
Post-processing module, obtained by being solved to the 3D CFD models by vascular tree of the solver to being established
To the fluid parameter post-processed, to obtain the blood flow reserve fraction of vascular tree everywhere, wherein,
The entrance boundary condition and the export boundary condition, at least the three-dimensional computations hydrodynamics mould for vascular tree
Type is individual specificity's to obtain those boundary conditions of the fluid parameter for post processing.
The CFD model of the vascular tree can be used for asking sometimes together with corresponding entrance boundary condition and export boundary condition
Solution obtains the fluid parameter of blood vessel everywhere, but is used only for simulation sometimes and obtains intermediate data, such as determining the 3D
The entrance boundary condition and export boundary condition of CFD model, the entrance and exit boundary condition well embody individual specificity.
No matter which kind of purposes embodies the entrance boundary condition and export boundary condition of individual specificity, it can improve what simulation calculated
Accuracy.
Preferably, for the 3D CFD models so that the entrance boundary condition that solution obtains fluid parameter is each
The flow of inlet, and determined according to individual medical image.
Preferably, the CFD model for calculating the vascular tree that grid generation module is established further includes one-dimensional (1D)
CFD model;The solver also solves the 1D CFD models of the vascular tree, is joined with obtaining the fluid of vascular tree everywhere
Number;The boundary condition setup module will obtain the entrance of fluid parameter for the 3D CFD models of the vascular tree for solving
The entrance boundary condition at place and/or the export boundary condition in exit are arranged to the solver to the vascular tree
1D CFD models solve obtained vascular tree and correspond to the fluid parameter everywhere.The fluid parameter usually can be with
It is any one of flow, speed and pressure of blood, entrance and exit boundary condition can set identical or different kind
Fluid parameter, as long as convergence is preferable when solving.
Preferably, the export boundary condition of fluid parameter is obtained as individual for solving for the 3D CFD models
Vascular tree each exit capilary resistance.Alternatively, the export boundary condition passes through the 1D CFD moulds to vascular tree
Type solves to obtain, and particularly, the CFD model for calculating the vascular tree that grid generation module is established further includes 1D
CFD model;The solver also solves the 1D CFD models of the vascular tree, to obtain the pressure of vascular tree everywhere;
The export boundary condition in the exit of the 3D CFD models of the vascular tree is arranged to by the boundary condition setup module
The solver carries out the 1D CFD models of the vascular tree to solve the pressure of obtained vascular tree correspondence everywhere.
The system can calculate the capilary resistance in each exit of the vascular tree of individual by a variety of methods.Example
Such as, in the system, the boundary condition setup module is also set for the 3DCFD models of the vascular tree under quiescent condition
Entrance boundary condition and export boundary condition;The solver is also based on the entrance boundary condition and outlet border under quiescent condition
Condition solves the 3D CFD models of the vascular tree, to obtain the pressure and flow in each exit of vascular tree;It is described
Boundary condition setup module is carried out according under entrance boundary condition of the solver under quiescent condition and export boundary condition
Solve the pressure in each exit of obtained vascular tree and flow or the entrance according to the solver under quiescent condition
Carry out solving the flow in each exit of obtained vascular tree and the tranquillization shape under boundary condition and export boundary condition
Outlet pressure under state, to calculate the capilary resistance in each exit of the vascular tree of individual.
Preferably, for the vascular tree 3D CFD models set entrance boundary condition and export boundary condition be not limited to
Entrance boundary condition and export boundary condition under quiescent condition can also be according to individual medical image and/or empirical equation
The three-dimensional computations fluid mechanic model of the vascular tree sets initial entrance boundary condition and export boundary condition, the solution
Device is also based on the initial entrance boundary condition and export boundary condition to the three-dimensional computations hydrodynamics mould of the vascular tree
Type is solved, to obtain the pressure and flow in each exit of vascular tree.These initial entrance boundary conditions and outlet side
The fluid parameter that boundary's condition solves cannot directly supply post processing, but accurate for the capilary resistance for calculating exit
Degree is much of that, and the capilary resistance in obtained exit can embody the specificity of individual, and energy as export boundary condition
Enough realize the good convergence of CFD model.The boundary condition setup module is according to the solver on the initial entrance side
Carry out solving the pressure and flow in each exit of obtained vascular tree under boundary's condition and export boundary condition, to calculate individual
Vascular tree each exit capilary resistance
For another example in the system, the CFD model for calculating the vascular tree that grid generation module is established is also
Including 1D CFD models;The solver also solves the 1D CFD models of the vascular tree, to obtain vascular tree everywhere
Flow and pressure;The boundary condition setup module seeks the 1D CFD models of the vascular tree according to the solver
It solves obtained vascular tree and corresponds to flow and pressure everywhere, to calculate the resistance of the capilary in each exit of the vascular tree of individual
Power.
Preferably, for the 3D CFD models for solve obtain the export boundary condition of fluid parameter can not also be according to
Rely in the calculating of the capilary resistance of vascular tree, but calculate to realize by the coupling of 3D-1D.Particularly, the calculating
Grid generation module calculates grid for the entirety of the geometrical model of the vascular tree or branch's generation, so as to establish the vascular tree
Whole 1D CFD models, and establish the entirety of the vascular tree or the 3D CFD models of branch;The boundary condition is set
It is the blood at congestive state lower inlet that module, which is put, by the entrance boundary condition setting of the whole 1D CFD models of the vascular tree
Its export boundary condition is arranged to vein pressure by flow, and when the 3DCFD models for the entirety for establishing the vascular tree are with supply and demand
When solution obtains fluid parameter, corresponding entrance boundary condition setting is congestive state lower inlet by the boundary condition setup module
The fluid parameter at place;And when the 3D CFD models for the branch for establishing the vascular tree obtain fluid parameter for solving, it is described
Corresponding entrance boundary condition setting is whole 1D of the solver to the vascular tree by boundary condition setup module
CFD model carries out solving the fluid parameter of obtained vascular tree correspondence everywhere.
According to the second aspect of the invention, provide a kind of CFD model for vascular tree and corresponding entrance boundary item is set
The method of part and export boundary condition, the CFD model is by obtaining medical image, performing and split and rebuild individual vascular tree
Geometrical model, the then geometrical model generation for the vascular tree calculates grid and establishes, and including 3D CFD models,
In, at least for vascular tree three-dimensional computations fluid mechanic model the entrance boundary condition and the export boundary condition it is equal
For individual specificity's.
The entrance boundary condition can be arranged to individual specificity's in several ways.Preferably, for the 3D
CFD model is the flow of each inlet for the entrance boundary condition that solution obtains fluid parameter, and according to individual
Medical image determines.
It is highly preferred that the CFD model for the vascular tree established further includes 1D CFD models, the method further includes:
The 1D CFD models of the vascular tree are solved, to obtain vascular tree fluid parameter everywhere and the 3D will be used for
CFD model obtains the entrance boundary condition of inlet of fluid parameter and/or the outlet side in exit for solving
Boundary's condition setting be by the 1D CFD models to the vascular tree solve obtained vascular tree correspond to everywhere described in
Fluid parameter.The fluid parameter can be any one of flow, speed and pressure of blood.Entrance boundary condition and go out
The fluid parameter of identical or different species may be employed in mouth boundary condition, as long as ensuring good convergence process.
Optionally, for the 3D CFD models export boundary condition of fluid parameter is obtained for solving
For the capilary resistance in each exit of the vascular tree of individual.The capilary resistance in each exit of the vascular tree of the individual
Power can calculate by various modes.Preferably, the capilary resistance in each exit of individual vascular tree is calculated as below
It obtains:For the entrance boundary condition and export boundary condition under the 3D CFD models setting quiescent condition of the vascular tree;It is based on
Entrance boundary condition and export boundary condition under quiescent condition solve the 3D CFD models of the vascular tree, to obtain
The pressure and flow in each exit of vascular tree;It is carried out according under the entrance boundary condition under quiescent condition and export boundary condition
Solve the flow in each source of obtained vascular tree and the outlet pressure under the quiescent condition or according to quiescent condition
Under entrance boundary condition and export boundary condition under solve the pressure and flow in each source of obtained vascular tree, come
Calculate the capilary resistance in each exit.
Preferably, the capilary resistance in each exit of individual vascular tree can be calculated as below to obtain:According to individual
Medical image and/or empirical equation set initial entrance boundary condition for the three-dimensional computations fluid mechanic model of the vascular tree
And export boundary condition;Based on the initial entrance boundary condition and export boundary condition to the three-dimensional computations of the vascular tree
Fluid mechanic model is solved, to obtain the pressure and flow in each exit of vascular tree;According in the initial entrance
Carry out solving the pressure and flow in each exit of obtained vascular tree under boundary condition and export boundary condition, it is a to calculate
The capilary resistance in each exit of the vascular tree of body.
Furthermore it is preferred that the capilary resistance in each exit of the vascular tree of individual is calculated as below to obtain:Described in foundation
The 1D CFD models of vascular tree;The 1D CFD models of the vascular tree are solved, with obtain the flow of vascular tree everywhere and
Pressure;Solve obtained vascular tree according to the 1D CFD models to the vascular tree and correspond to flow and pressure everywhere,
To calculate the capilary resistance in each exit of the vascular tree of individual.
Preferably, the method can also be set in the following way obtains for the 3D CFD models for solving
The export boundary condition of fluid parameter:The CFD model for the vascular tree established further includes 1D CFD models;To the blood vessel
The 1D CFD models of tree are solved, to obtain the pressure of vascular tree everywhere;It and will be for the 3D CFD moulds of the vascular tree
Type is arranged to for solving the export boundary condition in the exit for obtaining fluid parameter through the 1D to the vascular tree
CFD model carries out solving the pressure of obtained vascular tree correspondence everywhere.
Preferably, the method can also set the entrance of the 3DCFD models in a manner that 3D-1D is coupled and calculated
Boundary condition:The CFD model for the vascular tree established includes the whole 1D CFD models of the vascular tree and the blood
The entirety of Guan Shu or the 3D CFD models of branch;By the entrance boundary condition setting of the whole 1D CFD models of the vascular tree
For the blood flow at congestive state lower inlet, its export boundary condition is arranged to vein pressure;And work as and establish the vascular tree
Whole 3D CFD models when obtaining fluid parameter for solving, by corresponding entrance boundary condition setting under congestive state
The blood flow of inlet;And when the 3D CFD models for the branch for establishing the vascular tree obtain fluid parameter for solving, it will
Corresponding entrance boundary condition setting is to carry out solving obtained blood by the whole 1D CFD models to the vascular tree
The corresponding flow everywhere of pipe tree.
Description of the drawings
Below with reference to the accompanying drawings each embodiment of the present invention is described in detail, wherein
Fig. 1 shows computation hydrodynamics according to an embodiment of the invention and calculates blood flow reserve fraction to simulate
The block diagram of system;
Fig. 2 shows the schematic diagram for the flow that each component of system as shown in Figure 1 sequentially performs;
Fig. 3 (a)-Fig. 3 (c) shows the example of CFD model and boundary condition in accordance with another embodiment of the present invention,
In, Fig. 3 (a) shows the relation between the sectional view of medical image and the geometrical model that segmentation is rebuild, and Fig. 3 (b) shows segmentation weight
The 3-D geometric model of the obtained vascular tree including sustainer, left and right artery and vessel branch is built, Fig. 3 (c) shows left dynamic
The CFD model of the vascular tree of arteries and veins and corresponding boundary condition, including vascular wall and entrance and exit boundary condition.
Fig. 4 shows that the system and method that progress 3D-1D couplings in accordance with another embodiment of the present invention calculate are obtained
The schematic diagram of left artery and its 3D-1D coupling networks of branch.
Fig. 5 shows that the system and method that the progress 3D-1D couplings of another embodiment according to the present invention calculate are counted
The diagram of FFR in the left artery of obtained individual everywhere.
Specific embodiment
As shown in Figure 1, computation hydrodynamics according to an embodiment of the invention calculates blood flow reserve point to simulate
Several systems 1, the system 1 include:Vessel tree model generation module 11, for obtaining medical image, performing and split and rebuild
The geometrical model of the vascular tree of individual, for example, geometrical model such as Fig. 3 that the sectional view of acquired medical image and segmentation are rebuild
(a) shown in;Grid generation module 12 is calculated, grid is calculated for being generated for the geometrical model of the vascular tree, so as to establish
The CFD model of the vascular tree, the CFD model include 3D CFD models;Boundary condition setup module 13, for be established
The CFD model of the vascular tree corresponding entrance boundary condition and export boundary condition are set;Attribute setup module 14, is used for
The physical attribute and flow equation of blood are set;Solver 15, for being based on the entrance boundary condition and outlet border item
Part, set physical attribute and flow equation, the CFD model of the vascular tree to being established solves, to obtain blood
The fluid parameters of Guan Shu everywhere;And post-processing module 16, for by 3D of the solver to the vascular tree established
CFD model, which carries out solving the obtained fluid parameter, to be post-processed, to obtain the blood flow reserve fraction of vascular tree everywhere,
Wherein, 3D CFD models are at least by the boundary condition module and the entrance boundary condition of individual specificity and described is set
Export boundary condition supplies the fluid parameter of the vascular tree post-processed everywhere to calculate.Therefore ensure that individual vascular tree everywhere
The accuracy of fluid parameter, and then the fluid ginseng in vascular tree at each vessel branch of clinical concern can be predicted exactly
Number.
Above-mentioned each component of the system is completed complete FFR and is calculated jointly, is related to multiple important steps, including medicine
The processing of image, the segmentation of vascular tree and three-dimensional reconstruction, geometrical model and CFD calculate the generation of grid, CFD simulations calculate and after
Processing etc..Illustrate the flow that each component of system shown in FIG. 1 sequentially performs by taking Fig. 2 as an example below.
The pre-treatment that FFR is calculated is completed by vessel tree model generation module 11, based on the medical image to individual, such as
Meet the various medical images of DICOM specifications blood vessel analyzed and split as a result, rebuild individual vascular tree geometry
Model.The vessel tree model generation module 11, for example, can utilize include vmtk (The Vascular Modeling
Toolkit various open source softwares), based on that the blood vessel is analyzed and split as a result, calculating and rebuilding generation blood vessel
The vessel centerline and vascular wall of tree, so as to form the geometrical model of the vascular tree.
The CFD calculating sections flow that FFR is calculated by the calculating grid generation module 12, boundary condition setup module 13,
Attribute setup module 14 and the cooperation of solver 15 are completed.
Particularly, the grid generation module 12 that calculates is the blood vessel using the various open source softwares for including vmtk
The geometrical model generation of tree meets the high quality calculating grid that CFD calculates requirement, so as to establish the CFD model of the vascular tree,
In general, the 3D CFD models of the vascular tree are established for calculating the fluid parameter of vascular tree everywhere.If described in simple foundation
Though the 1D CFD models of vascular tree simultaneously calculate the fluid parameter of vascular tree everywhere rapidly accordingly, it can almost be regarded as real
When, but the accuracy for calculating the vascular tree fluid parameter everywhere of gained is still not enough to be directly used in the clinics of vascular conditions and examines
It is disconnected, but be enough as the boundary condition at the corresponding entrance and exit of the 3D CFD models of vascular tree, and fully demonstrate
The specificity of individual can greatly improve the accuracy that the simulation of 3D CFD models calculates.
The boundary condition setup module 13 sets the boundary condition for meeting Human physiology characteristic for individual patient.
CFD calculates grid and is regarded as being made of the point for spatially representing actual geometrical model.The intact vessel train of thought of accurate human body
CFD calculating at least need hundreds thousand of grids, be related to millions of equations, matrix computations, therefore, on engineer application almost
Complete human vas train of thought can not be simulated.It is general in practical application only to consider to the interested subregion of result, i.e., by side
Define justice subregion, hereinafter we be referred to as vascular tree, the visual clinical demand of the vascular tree and it is different, for example,
Can be complete vascular tree coronarius, entrance boundary (being hereinafter also referred to as " entrance ") be defined as left and right sustainer with
The crotch of sustainer, outlet border (being hereinafter also referred to as " exporting ") are defined as the capilary of vascular tree coronarius
Place.By the setting of boundary condition, different vascular trees can also be defined, below only with complete vascular tree coronarius
Exemplified by illustrate that the individual specificity of boundary condition is set, but vascular tree is without being limited thereto.
Fig. 3 (a)-Fig. 3 (c) shows the example of CFD model and boundary condition in accordance with another embodiment of the present invention,
In, Fig. 3 (a) shows that the geometrical model that the sectional view of medical image, segmentation rebuild and relation between the two, Fig. 3 (b) are shown
The 3-D geometric model of the obtained vascular tree including sustainer, left and right sustainer and major vessel branches, Fig. 3 are rebuild in segmentation
(c) CFD model of the vascular tree of aorta sinistra and corresponding boundary condition are shown, including vascular wall and entrance and exit perimeter strip
Part (outlet 1,2 ..., 8).
In the prior art, the boundary condition setup module 13 can only be by counting or assuming the inlet of left and right sustainer
Flow proportional mode, to set the entrance boundary condition of the vascular tree.The prior art also for example sets common access stream
Amount:It is measured and simulated using high-volume clinical case FFR and compared, determine the optimal value to left and right sustainer inlet flow rate, application
The value is in the calculating of all cases.This setup cannot embody the specificity of individual, statistics or the entrance side for assuming gained
Boundary's condition may not be suitable for individual, so as to reduce the accuracy of CFD result of calculations.
The present invention provides the sides that the CFD model for vascular tree sets corresponding entrance boundary condition and export boundary condition
A variety of examples of method realize the individual specificity of entrance boundary condition and the export boundary condition.The example of these methods
It can be used and (hereinafter how be given when describing the example of these methods by system 1 by the above correlation module in system 1
In correlation module come the explanation implemented, but be not intended to limit the module of implementation) or by additional in the system
Module is transferred to the boundary condition and sets using determining the entrance boundary condition of individual specificity and export boundary condition
Module 13 is for using and set.
Attribute setup module 14 is used to set the physical attribute and flow equation of blood, the physical attribute example of the blood
Such as but blood is not limited to as Newtonian fluid and is glued using laminar flow, the density of blood for the individual for meeting Human physiology characteristic and blood flow
Property etc., flowing of the flow equation for example including the blood for UNSTEADY FLOW, based on can not press stream Navier-Stokes
Equation group etc..User can carry out being manually entered what is set either default setting or carry from system as the case may be
It is selected in the physical attribute of blood and flow equation storehouse.
Solver 15 be used for based on set entrance boundary condition and export boundary condition, set physical attribute and
Flow equation, the CFD model of the vascular tree to being established solve, to obtain intermediate result or obtain vascular tree everywhere
Fluid parameter.It is three-dimensional that the CFD model can include a peacekeeping, as described above, in order to ensure the accuracy of result of calculation,
1D CFD models solve the fluid parameter of the vascular tree of gained everywhere generally only as intermediate result, are not used directly for blood vessel
The clinical diagnosis of pathological state.3D CFD models can then solve under corresponding entrance and exit boundary condition and obtain vascular tree
Fluid parameter everywhere is then only used for solving obtaining intermediate result rather than blood for its entrance and exit boundary condition set sometimes
The fluid parameters of Guan Shu everywhere.
The post processing that FFR is calculated is realized using post-processing module 16.Particularly, the post-processing module 16 be used for by
The 3D CFD models of vascular tree of the solver 15 to being established are carried out after solving the obtained fluid parameter progress
Reason, to obtain the blood flow reserve fraction (FFR) of vascular tree everywhere.For example, it is closed according to what is solved to 3D CFD models as a result, utilizing
It is formula FFR=Pd/Pa(PdFor the pressure of blood vessel everywhere, PaFor the average pressure at sustainer), blood vessel can be calculated everywhere
FFR.It is alternatively possible to it is post-processed by the result solved to 3D CFD models, to generate analysis of cases report, preferably
Ground, the FFR values of stenosis distal end can write in analysis of cases report as final calculated value, if FFR<0.75, then show
This is narrow to cause the functional ischemic of blood vessel, it is necessary to carry out interventional treatment.The result solved to 3D CFD models further includes it
The parameter of his clinical concern, as pull back (pull-back) curve, blood vessel blood flow velocity everywhere and pressure distribution, blood vessels of FFR is each
Flow distribution, the shear stress of vascular wall everywhere of branch etc., the parameter of these clinic concerns can also optionally write case
In analysis report, for the reference of doctor's behaviours diagnosis.
The CFD model setting that the present invention is vascular tree is detailed below and obtains the fluid ginseng of blood vessel everywhere for solving
Several corresponding entrance boundary conditions and a variety of examples of the method for export boundary condition, each example realize entrance boundary
The individual specificity of condition and the export boundary condition.At least for vascular tree 3D CFD models entrance boundary condition and
The export boundary condition can be derived by by the medical image of individual or the Fluid Mechanics Computation mould via individual vascular tree
The solution of type accordingly obtains.Especially, in order to solve to obtain the fluid parameter for post processing, for of 3D CFD models
The entrance boundary condition and export boundary condition of body specificity can be derived according to the medical image or historical data of individual
It arrives, can also accordingly be obtained via the solution of 1D CFD models, in initial entrance and can also be gone out by 3D CFD models
Solution under mouthful boundary condition is drawn.
Entrance boundary condition
By taking vascular tree coronarius as an example, this method individually calculates left and right sustainer, and entrance boundary is left and right active
The crotch of arteries and veins and sustainer (branch) sets the respective blood flow of the left and right sustainer of the crotch as entrance boundary condition,
Two ways can be specifically subdivided into.
First way is, according to the medical image of individual, to calculate the myocardial mass of the case, utilize myocardial mass and master
The relation of blood flow volume, obtains aortic blood flow, and left and right sustainer shunts (aorta sinistra shunting 60% and the right side in proportion
Aorta shunting 40%), the blood flow that the left and right sustainer of the infall each shunts so is obtained, and is set to entrance
Boundary condition.Note that this is only to determine the flow of each inlet example according to the medical image of individual, utilize
Correspondence between myocardial mass and aortic blood flow according to the cardiovascular difference set to be modeled, can also use
Other between the blood flow of information and inlet blood vessel of other individual medical images in other medical images are right
It should be related to the blood flow for calculating each blood vessel in inlet, and be set to the entrance of the CFD model of the specific cardiovascular tree
Boundary condition, it is noted that this set method of entrance boundary condition can be adapted for both 1D CFD models and 3D CFD models.
(such as vessel tree model in system 1 can also be utilized to give birth to by establishing the complete 1D CFD models of vascular tree
Into module 11 and calculate grid generation module 12), and (such as utilizing the solver 15 in the system 1) is solved to it to calculate
Go out flow everywhere in complete vascular tree, then corresponding flow everywhere is applied in the complete vascular tree obtained by by the calculating
At the entrance boundary of the 3D CFD models of the vascular tree, the entrance boundary condition for being arranged to 3D CFD models (such as utilizes institute
State the boundary condition setup module 13 in system 1).The second way compensates for the defects of first way well:Sometimes
The blood flow of inlet blood vessel can not be determined according to the medical image of individual.Although the second way by taking flow as an example, also may be used
Using any one of flow, pressure and speed using blood or two kinds as entrance and exit boundary condition.Entrance
It is also possible using fluid parameter of the same race with export boundary condition, as long as ensuring that the convergence of 3D CFD models meets clinical want
It asks.Even without the entrance and exit boundary condition that individual specificity is set for the 1D CFD models, in view of the 1D CFD models
The blood vessel feature of individual itself is reflected, by CFD iterative solutions, 3D is corresponded in the complete vascular tree that when convergence is obtained
Flow everywhere at the entrance boundary of CFD model also well reflects individual specificity, so as to also just be set for 3D CFD models
The entrance boundary condition of individual specificity is put.1D CFD calculating only needs considerably less number of grid, and calculating speed is very fast
Fast (being regarded as calculating in real time) is calculated by 1D CFD promoting to obtain the entrance boundary condition of 3D CFD models as a result,
While result of calculation accuracy, also without the calculating speed for the entire flow for influencing the fluid parameter such as FFR.
Export boundary condition
There is the set-up mode of two class export boundary conditions:
First way is, by export boundary condition be arranged to for individual patient vascular tree coronarius everywhere
Capilary resistance;The second way is, similar to second of set-up mode of entrance boundary condition, to apply 1D CFD result of calculations
In the exit of 3D CFD models, the export boundary condition of 3D CFD models is arranged to the solution of the 1D CFD models of vascular tree
Corresponding pressure everywhere in the blood vessel of gained.
First way can be subdivided into two kinds of submodes:
I. calculated based on 3D CFD
Using the dynamic (dynamical) phenomenon of human bloodstream, i.e., (Resting Condition) vessel segment under quiescent condition
It is the fact that pressure change very small (no matter vessel segment whether there is narrow), the entrance boundary condition setting of 3D CFD models is defeated for the heart
Average discharge Q in the output cyclein, its export boundary condition is arranged to corresponding average pressure pout, by 3D CFD moulds
Type is solved under the entrance and exit boundary condition under as above quiescent condition, to calculate the stream predicted at each outlet border
Measure Qi,out(i is the number of each outlet), utilizes relational expression Pout=Qi,outRiCalculate the capilary resistance at each outlet border
Ri。
It is of course also possible to according to obtained by being solved under the entrance boundary condition under quiescent condition and export boundary condition
The pressure P in each exit of vascular tree arrivedi,outWith flow Qi,out, utilize relational expression Pi,out=Qi,outRiCalculate each outlet
The capilary resistance R of boundaryi。
Although the capilary resistance in each exit is calculated using the entrance and exit side under quiescent condition based on 3D CFD
Boundary's condition is relatively convenient, but also has other clinical feasible methods.For example, the entrance and exit boundary condition can be root
According to the obtained initial entrance boundary condition of individual medical image and/or empirical equation and export boundary condition (herein with " just
Begin " come with being differentiated for the setting of 3D CFD models directly to calculate those boundary conditions of the fluid parameter for post processing).
For example, can average cardiac output be calculated according to the ultrasonic Doppler image of such as individual heart, according on sustainer
Valve blood pump out rear flow how the empirical equation for the ratio distributed in vessel branches at different levels, can obtain corresponding (all everywhere
Such as left and right sustainer, capilary) flow as initial entrance boundary condition and export boundary condition.
It is calculated based on 1D CFD:1D CFD calculating only needs considerably less number of grid, and calculating speed is very fast (can to recognize
To calculate in real time), using the advantage, the 1D CFD models of complete coronary arterial tree are established, according to any of the above-described kind of entrance side
The set-up mode of boundary's condition carrys out inlet porting boundary condition, and export boundary condition can then be all set to vein pressure and (be usually
0).Using calculating of the 1D CFD under the entrance boundary condition and export boundary condition of the setting, obtain in 1D CFD models
The corresponding flow and pressure everywhere of each outlet of corresponding 3D CFD models, utilizes relational expression Pi,out=Qi,outRiIt calculates
The capilary resistance R of each outlet border of 3DCFD modelsi.Fig. 4 is shown as the schematic diagram of 3D-1D coupling networks.This application can pole
The big solution simplified during 3D CFD are calculated to export boundary condition capilary resistance.
Then pressure in complete vascular tree coronarius everywhere is calculated using 1D CFD in the second way, and will
Corresponding pressure applies to the corresponding outlet border of 3D CFD models.Certainly, except pressure, can also use other kinds of
Fluid parameter, including flow, speed etc., the export boundary condition as 3D CFD models.
As above in order to set, the 1D CFD involved by export boundary condition are calculated and 3D CFD calculating can be by system 1
In vessel tree model generation module 11, calculate grid generation module 12, boundary condition setup module 13, attribute setup module 14
It is realized with solver 15 to coordinate.For example, the model life of 1D and 3D can be integrated respectively in each above-mentioned component in system 1
Into the generation of, grid, the function that boundary condition is set, attribute sets and solves, can also be included in each above-mentioned component for 1D
With the discrete sub- component of 3D, realize that 1D CFD are calculated together by each sub- component of 1D, it is real together by each sub- component of 3D
Existing 3D CFD are calculated.
For example, solver 15 can also include 1D CFD solvers and 3D CFD solvers, available for blood vessel is calculated as above
Entrance boundary condition and/or export boundary condition (the 3DCFD models and the corresponding entrance boundary condition of the 3D CFD models of tree
It is used to calculate fluid parameter in vascular tree everywhere with export boundary condition for clinical diagnosis), it can also be used to calculate vascular tree
In fluid parameter everywhere for clinical diagnosis.
Based on high-performance calculation CFD solvers, it can simulate and calculate left and right aorta vessel tree congestive state
(Hyperemia) other relevant blood flow parameters of FFR and clinical concern under.
For example, apply congestive state bottom left, aorta dextra blood flow as entrance in the boundary condition setup module 13
Boundary condition, here, congestive state bottom left, aorta dextra blood flow can be according to the medical images and/or empirical equation of individual
It determines, can also be drawn by 1D CFD model solutions, it can also by 3D CFD models on initial inlet and outlet border
Under the conditions of solution show that example has been presented above, and apply at each outlet border that as above various modes obtain
Capilary resistance as export boundary condition in the case of, the attribute setup module 14 can set vascular fluid be newton
Fluid and for laminar flow, sets the density of blood for meeting Human physiology characteristic and blood flow viscosity, by solver 15 to vascular tree
3D CFD models calculate UNSTEADY FLOW until result restrains, so as to obtain the fluid parameter of vascular tree everywhere.Post-processing module
The step of 16 execution, details are not described herein, and FFR everywhere in vascular tree is calculated accordingly.
Above method calculates FFR, and the calculating dependent on capilary resistance coronarius is as export boundary condition, and most
It models and calculates by the 3D CFD to complete vascular tree eventually to obtain the fluid parameter in vascular tree everywhere.In practice,
Need to be further simplified and accelerate the setting of boundary condition, it is also desirable to further speed up total calculating speed of fluid parameter.This
Outside, doctor require CFD modelings can not only accurate simulation calculate complete coronary artery, it is specific can also individually to simulate calculating
Vessel branch (such as suspected lesion section), and the accuracy calculated the simulation of coronary artery everywhere is distinguishing.For example, doctor
It can wish to obtain complete relatively accurate FFR distributions coronarius by real-time 1D CFD Modeling Calculations come rapid, make
Obtaining therefrom manual or automatic can select the specific vessel branches for needing more clinical concerns, and the specific blood vessel is divided
Zhi Jinhang is relatively time-consuming but the higher 3D CFD Modeling Calculations of accuracy.
The present invention provides a kind of progress 3D-1D couplings calculating, and (inventor is referred to as virtual digit and subtracts shadow angiocardiography
Art) system and boundary condition setting method, with upgrade requirement more than satisfaction, which is combined 3D-1D
The characteristics of CFD and respective advantage greatly simplify the requirement to boundary condition.This method not only equally can relatively accurately mould
Intend calculating complete coronary artery, it more can be more accurate and pointedly individually simulation calculates specific vessel branch (such as doubtful disease
Become section), therefore closer greatly shorten the specific blood vessel for calculating required time and preferentially ensuring more clinical concerns
The accuracy that the simulation of branch calculates.
This method is Newtonian fluid using the blood in the setting blood vessel of attribute setup module 14 and is laminar flow, and setting meets people
Density of blood and the blood flow viscosity of body physiology characteristic, and flow equation is arranged to UNSTEADY FLOW until result is for solving
Device 15 is restrained.Specific requirements depending on solving fluid parameter, 3D CFD grids and model for this specific requirements can be entire
The vascular tree generation of left and right sustainer, or the specific vessel branch generation of certain section;1D CFD grids and model have been then
Whole vascular tree coronarius generation, the step can pass through the vessel tree model generation module 11 in system 1 and calculating net
Lattice generation module 12 is realized.
It is that the 1D CFD models of vascular tree and 3D CFD models are set accordingly using the boundary condition setup module 13
Boundary condition.Wherein, it is congestive state bottom left by the entrance boundary condition setting of 1D CFD models, aorta dextra blood flow, it will
The export boundary condition of 1D CFD models is arranged to vein pressure (being usually 0).
It is the vascular tree of entire left and right sustainer for what 3D CFD models were simulated or is only some specific blood vessel
Branch sets different entrance boundary conditions.Particularly, the vascular tree of entire left and right sustainer is such as simulated, by its entrance side
Boundary's condition setting is application congestive state bottom left, aorta dextra blood flow;As simulated certain specific vessel branch, then by its entrance side
Boundary's condition setting be using 1D CFD models under the entrance boundary condition and export boundary condition (vein pressure) being set as above
In the flow that correspondence position is calculated, other fluid parameters in addition to flow can also be used here.
The export boundary condition of 3D CFD models is arranged to using 1D CFD models in the entrance boundary item being set as above
The pressure value being calculated under part and export boundary condition (vein pressure) in correspondence position can also be used except pressure here
Other fluid parameters in addition.
So realize that the coupling of 3D-1D CFD calculates.Particularly, grid, boundary condition based on foregoing 1D and 3D
Processing, calculated using 1D CFD solvers, by the result of calculation of 1D CFD, according to the foregoing place to 3D boundary conditions
Reason, is applied to the boundary of 3D, is calculated using 3D CFD solvers, is joined with the fluid needed for obtaining in vascular tree everywhere
Number.
During 3D, 1D CFD are calculated, the Navier-Stokes (N-S) that can not press stream is respectively provided with using attribute setup module 14
Equation group (quality, momentum conservation equation):
▽ u=0
U is fluid velocity vectors, and p is pressure, and ρ is fluid density, and υ is the kinematic viscosity of fluid, so as to 1DCFD solvers
With 3D CFD solvers (or solver 15 of the function by incorporating 1D CFD and 3D CFD solutions) based on the N-S that can not press stream
Equation group solves.
CFD calculates post processing:According to the 3D CFD being calculated as a result, utilizing relational expression FFR=Pd/Pa(PdIt is each for blood vessel
The pressure at place, PaFor the average pressure at sustainer), the FFR of blood vessel everywhere can be calculated.It is alternatively possible to by 3D
The result that CFD model solves is post-processed, to generate analysis of cases report, it is preferable that the FFR values of stenosis distal section can be made
It is write for final calculated value in analysis of cases report, if FFR<0.75, then show that the narrow functionality that can cause blood vessel lacks
Blood is, it is necessary to carry out interventional treatment.The parameter of other clinic concerns is further included to the result that 3D CFD models solve, as FFR pulls back
(pull-back) curve, blood vessel blood flow velocity everywhere and pressure distribution, the flow distribution of each branch of blood vessel, vascular wall be everywhere
Shear stress etc., parameters of these clinic concerns can also be optionally write in analysis of cases report, so that doctor's behaviours are examined
Disconnected reference.
By taking the aorta sinistra being intended in individually simulation coronary artery as an example, Fig. 4 shows to be included for complete coronary artery
The 3D CFD grids of the aorta sinistra generation of the 1D CFD grids of distal end Angiogenesis and only clinical concern, by above-mentioned
3D-1D couplings FFR distributions everywhere are calculated in aorta sinistra, as shown in Figure 5, have in practice been verified that the 3D-1D is coupled
Result of calculation provide sufficient accuracy and it is clinically-acceptable take, can control within 1 hour.
The present invention is highly integrated and automates entire calculation process, the reconstruction of processing, computation model including medical image,
Calculating in itself and the post processing of result, realizes the quick and precisely calculating of high-volume case.Such as and it is unrestricted, C/C+ can be used
+ language is developed and optimizes above-mentioned 1D and/or 3D CFD simulation calculation process and it couples calculation process, each in system 1
A module and its realized entire calculation process (including FFR calculate pre-treatment, FFR calculate CFD calculating sections and
Post processing, interpretation of result and report generation that FFR is calculated etc.) automations such as Python, Bash Script can be used to realize,
So as to while analog result accuracy is ensured, improve user-friendliness.
The explanation of above example is only intended to help to understand core concept of the invention.It should be pointed out that for this field
Those of ordinary skill for, without departing from the principle of the present invention, can also to the present invention carry out it is several improvement and repair
Decorations, but these improvement and modification are also fallen into the range of the claims in the present invention are claimed.