(57)【要約】      （修正有）【課題】患者別の血管の情報を決定する方法に関する。より詳細には、簡略化された冠循環モデルを適用した患者別の心臓血管情報を決定する方法に関するものである。また、患者別に動脈から発源する血管の分枝に対する血流量を決定する方法に関する。【解決手段】コンピューターを利用して、心血管情報を決定する方法は、大動脈から発源する複数の冠状動脈を含む画像データの提供を受けるステップと、画像データを処理して複数の冠状動脈の３次元形状モデルを生成するステップと、生成された３次元形状モデルの血流をシミュレートするステップと、各冠状動脈の冠血流予備量比（ＦＦＲ）を決定するステップとを含んでいる。複数の冠状動脈の３次元形状モデルの血流のシミュレーションのステップで、計算流体力学モデルを適用して、集中パラメータモデルは、冠状動脈、冠状動脈の毛細血管、および冠静脈を含む簡略化された冠動脈循環モデルを使用する。【選択図】図２
2018187406.tif 000002
【特許請求の範囲】
【請求項１】
コンピューターシステムを利用して、患者の血管情報を決定するための方法であって、
動脈から発源する複数の血管の少なくとも一部を含む画像データの入力を受けるステップと、
前記入力された画像データを処理して、それぞれの血管の3次元モデルを生成するステップと、
前記生成された各血管の3次元モデルでは、それぞれの血管の分枝から末端までの血管の長さを求めるステップと、
前記求められた各血管の長さの割合に応じて、前記各血管に流れる血流量の割合を決定するステップとを含む血管情報を決定する方法。
【請求項２】
前記各血管の長さを求めるステップは、
前記各血管の3次元モデルの血管中心線を求めるステップと、
前記それぞれの血管の3次元モデルの断面積が一定の値以下になる末端点を求めるステップと、
前記それぞれの血管の分枝から末端地点までの中心線の長さを求めるステップとを含む請求項1に記載の血管情報を決定する方法。
【請求項３】
前記複数の血管は、右冠状動脈（RCA）、左前下行枝（LAD）および左回旋枝（LCX）を含み、
前記血流量を決定するステップにおいて、前記各血管を流れる血流量の割合は、以下の比例によって定める請求項1又は2に記載の血管情報を決定する方法。
2018187406.tif 000022
ここで、
QLAD：左前下行枝の血流量、 lLAD：左前下行枝の長さ
QLCX：左回旋枝の血流量、 lLCX：左回旋枝の長さ
QRCA：右冠状動脈の血流量、 lRCA：右冠状動脈の長さ
(lRCA)RV：右冠状動脈の右心室の血液供給の部分の長さ
(lRCA)LV：右冠状動脈の左心室の血液供給の部分の長さ
α：右心室の血液供給血管の補正係数
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、患者別の血管の情報を決定する方法に関するものである。より詳細には、簡略化された冠循環モデルを適用した患者別の心臓血管情報を決定する方法に関するものである。また、患者別の動脈から発源する血管の分枝に対する血流を決定する方法に関するものである。
【背景技術】
【０００２】
一般的に、冠血流予備量比（ＦＦＲ、FRACTIONAL FLOW RESERVE）は、冠動脈狭窄の機能的重症度（FUNCTIONAL SEVERITY）を評価する臨床指標として広く使用される。患者のＣＴ（コンピューター断層撮影、COMPUTED TOMOGRAPHY）データに基づいたコンピューターシミュレーション方法は、ＦＦＲ計算のための非侵襲的な方法である。この方法は、計算流体力学モデルを、心血管システムの集中パラメータモデルと結合して狭窄された冠状動脈（冠動脈）の血流力学について詳細な解析結果を提供することができる。
【０００３】
Ｐｉｊｌｓなどは、冠動脈疾患の臨床指標として完全に膨張した状態の微小血管の狭窄の前後の割合を意味する冠血流予備量比（ＦＦＲ）を紹介した。Ｐｉｊｌｓグループは、ガイドワイヤー技術を利用して、ＦＦＲの測定法を開発し、これが冠動脈狭窄の程度を評価する技術となった。
【０００４】
Ｋｉｍ等は、ＦＦＲ値を評価する上で、ＣＴ画像と患者の情報を利用した非侵襲的なシミュレーション方法を提示した。これは、大動脈と冠状血管の血流力学的な計算のための計算流体力学技術と全体の心血管システムの集中パラメータモデルが統合された方法である。このシミュレーションモデルの妥当性と有用性は、いくつかの研究を通じて検証された（Min et al, 2012; Koo et al。、2011）。しかし、Ｋｉｍ等によって開発されたモデルは、計算流体力学モデルに大動脈を含み、集中パラメータモデルに全体的な心血管システムを含んでおり、計算が複雑であり、多くのパラメータの識別が必要とされて、シミュレーションの不確実性が増加される。
【０００５】
一方、大韓民国特許第10-1524955号（発明の名称、患者別の血流のモデリング方法およびシステム）には、患者の心血管情報を決定するための方法およびシステムが開示されている。前記特許に開示された心血管情報の決定方法は、大動脈から発源する複数の冠状動脈の少なくとも一部を含む患者の解剖学的構造の幾何形態に関する患者別のデータを受信して、患者別のデータに基づいて、複数の冠状動脈の少なくとも一部を含む解剖学的構造の第１の部分を表す３次元モデルを作成し、少なくとも部分的に心筋組織の質量またはボリュームに基づいて解剖学的構造の第１の部分内の血流特性に関する物理学ベースのモデルを作成し、前記３次元モデルと物理学ベースのモデルに基づいて、解剖学的構造の第１の部分内の冠血流予備量比（ＦＦＲ：FRACTIONAL FLOW RESERVE）を決定するステップを含んでいる。
【０００６】
前記特許に開示された方法で、物理学ベースのモデルは、３次元モデルの境界を通じた血流を示す集中パラメータモデル（LUMPED PARAMETER MODEL）を使用する。また、集中パラメータモデルを適用する場合、心室の筋肉のボリュームに基づいて血流量の条件を定める。これは、冠状動脈が血液を供給する領域の心室筋のボリュームが大きい領域ほど血流量が多くなるという仮定に基づいている。また、前記特許に開示された方法は、心筋組織のボリュームを求め、スケーリング則（scaling law）を使用しなければならない。スケーリング則を使用するためには、３次元心室モデルのセグメンテーション（segmentation）が不可欠である。すなわち、前記特許の方法を適用するためには、心臓全体に対するセグメンテーション処理を実行する必要があるため、モデルの不確実性が大きくなる。特に、心室の筋肉の体積の場合、厚さの方向の形状が複雑で、セグメンテーションの精度が落ちる。
【０００７】
また、Ｋｉｍ等によって開発された方法や、前記特許に開示された方法は、計算流体力学モデル（COMPUTATIONAL FLUID DYNAMICS MODEL、以下ＣＦＤモデル）に大動脈が含まれており、集中パラメータモデルは体動脈、体静脈、肺血管、左心、右心などを含む閉回路で構成されている。また、血流力学解析のためにＣＦＤモデルと集中パラメータモデルは、いくつかの標準的な代表値を持つパラメータを使用して、これらのパラメータ（例えば、体動脈、体静脈、肺血管などの抵抗と静電容量値）は、それぞれの個々の患者に適用するには不向きである。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】大韓民国特許第10-1524955号、発明の名称、患者別の血流のモデリング方法およびシステム
【非特許文献】
【０００９】
【非特許文献１】Pijls、NH、VanGelder、B.、VanderVoort、P.、Peels、K.、Bracke、FA、Bonnier、HJ、Gamal、MI、1995 Fractional flow reserve：a useful index to evaluate the influence of an epicardial coronary stenosis on myocardial blood flow。 Circulation 92、3183e3193。
【非特許文献２】Kim、HJ、Vignon-Clementel、IE、Coogan、JS、Figueroa、CA、Jansen、KE、Taylor、CA、2010. Patient-specific modeling of blood flow and pressure in human coronary arteries。 Ann。 Biomed。 Eng。 38（10）、3195e3209。
【非特許文献３】Min、JK、Leipsic、J.、Pencina、MJ、Berman、DS、Koo、BK、van Mieghem、C.、Erglis、A.、Lin、FY、Dunning、AM、Apruzzese、P.、Budoff、MJ、Cole 、JH、Jaffer、FA、Leon、MB、Malpeso、J.、Mancini、GB、Park、SJ、Schwartz、RS、Shaw、LJ、Mauri、L.、2012年Diagnostic accuracy of fractional flow reserve from anatomic CT angiography。 JAMA 308、1237e1245。
【非特許文献４】Koo、BK、Erglis、A.、Doh、JH、Daniels、DV、Jegere、S.、Kim、HS、Dunning、A.、DeFrance、T.、Lansky、A.、Leipsic、J.、Min、JK、 2011 Diagnosis of ischemia-causing coronary stenoses by noninvasive fractional flow reserve computed from coronary computed tomographic angiograms。 Results from the prospective multicenter DISCOVER-FLOW（diagnosis of ischemia-causing Stenoses obtained via noninvasive fractional flow Reserve）study。 J. Am。 Coll。 Cardiol。 58、1989e1997。
【非特許文献５】Shim、EB、Chang、KS、1997 Numerical analysis of three-dimensional Bj?ork-Shiley valvular flow in an aorta。 J. Biomech。 Eng。 119（1）、45e51。
【非特許文献６】Shim、EB、Chang、KS、1994 Three-dimensional vortex flow past a tilting disc valve using a segregated finite element scheme。 Comput。 Fluid Dyn。 J. 3（1）、205e222。
【非特許文献７】Shim、EB、Kamm、RD、Heldt、T.、Mark、RG、2000。Numerical analysis of blood flow through a stenosed artery using a coupled multiscale simulation method。 Comput Cardiol。 27、219e222。
【非特許文献８】Einstein、A.、1906. Eine neue bestimmung der molek dimensionen。 Ann。 Phys。19（2）、289。
【非特許文献９】Schreiner、W.、Neumann、F.、Mohl、W.、1990。Therole of intramyocardial pressure during coronary sinus interventions：a computer model study。 IEEE Trans。 Biomed。 Eng。 37、956e967。
【非特許文献１０】Lim、KM、Kim、IS、Choi、SW、Min、BG、Won、YS、Kim、HY、Shim、EB、2009。Computational analysis of the effect of the type of LVAD flow on coronary perfusion and ventricular afterload。 J. Physiol。 Sci。 59（4）、307e316。
【非特許文献１１】Brown、AG、Shi、Y.、Marzo、A.、Staicu、C.、Valverde、I.、Beerbaum、P.、Lawford、PV、Hose、DR、2012年Accuracy vs. computational time：translating aortic simulations to the clinic。 J. Biomech。 45（3）、516e623。
【非特許文献１２】Taylor、CA、Fonte、TA、Min、JK、2013 Computational fluid dynamics applied to cardiac computed tomography for noninvasive quantification of fractional flow reserve：scientific basis。 J. Am。 Coll。 Cardiol。 61（22）、2233e2241。
【非特許文献１３】West、GB、Brown、JH、Enquist、BJ、1997。A general model for the origin of allometric scaling laws in biology。 Science 276（5309）、122e126。
【発明の概要】
【発明が解決しようとする課題】
【００１０】
本発明は、血管の血流解析のための新しい方法とシステムを提供することを目的とする。
【００１１】
本発明の第１の目的は、身体全体の閉鎖型の集中パラメータモデルを使用せずに、血流を解釈しようとする血管の部分についてのみＣＦＤモデルと集中パラメータモデルを適用して血流を解析するための新しい方法とシステムを提供することである。
【００１２】
本発明の第２の目的は、動脈からの分枝された血管が血液を供給する身体領域の体積や質量を求めず、動脈からの枝血管のそれぞれに流れる血流量の割合を定めるための新しい方法を提供することである。
【課題を解決するための手段】
【００１３】
本発明の一の実施形態によれば、心血管情報のうち、冠状動脈の冠血流予備量比（ＦＦＲ）を決定する方法が提供される。
【００１４】
本発明に係るコンピューターシステムを利用して、心血管情報を決定する方法は、大動脈から発源する複数の冠状動脈を含む画像データの提供を受けるステップと、前記画像データを処理して前記複数の冠状動脈の３次元形状モデルを生成するステップと、前記生成された複数の冠状動脈の３次元形状モデルの血流をシミュレーションするステップと、前記血流シミュレーションの結果として、前記の各冠状動脈の冠血流予備量比（ＦＦＲ）を決定するステップとを含んでいる。前記複数の冠状動脈の３次元形状モデルの血流のシミュレーションのステップで、冠状動脈の３次元形状モデルの計算流体力学モデルを適用して、計算流体力学モデルと結合されている集中パラメータモデルは、冠状動脈、冠状動脈の毛細血管、および冠状静脈を含む簡略化された冠動脈循環モデルを使用する。
【００１５】
本発明の他の実施形態によれば、前記血流シミュレーションのステップで、前記冠状動脈の３次元形状モデルの計算流体力学モデルを適用したときに、入口境界条件として大動脈血圧パターンを使用することができる。
【００１６】
本発明の他の実施形態によれば、前記血流シミュレーションのステップは、前記冠状動脈の３次元形状モデルの中心線の長さを求めるステップをさらに含むことができる。また、前記計算流体力学モデルに前記簡略化された冠動脈循環モデルを結合するとき、それぞれの冠状動脈に流れる血流量の割合で冠状動脈の毛細血管の抵抗値を定め、それぞれの冠状動脈に流れる血流量の割合は、冠状動脈の３次元形状モデルの中心線の長さの比率で定めることができる。
【００１７】
本発明の他の実施形態によれば、血管情報の中で、複数の血管に流れる血流量の割合を決定する方法が提供される。
【００１８】
本発明に係るコンピューターシステムを利用して、患者の血管情報を決定するための方法は、動脈から発源する複数の血管の少なくとも一部を含む画像データの入力を受けるステップと、前記入力された画像データを処理して、それぞれの血管の３次元モデルを生成するステップと、前記生成された各血管の３次元モデルでは、それぞれの血管の分枝から末端までの血管の長さを求めるステップと、前記求められた各血管の長さの比率にしたがって、前記の各血管に流れる血流量の割合を決定するステップとを含む。
【００１９】
前記の各血管の長さを求めるステップは、前記各血管の３次元モデルの血管中心線を求めるステップと、前記それぞれの血管の３次元モデルの断面積が一定の値以下になる末端点を求めるステップと、前記の各血管の分枝から末端地点までの中心線の長さを求めるステップとを含む。
【００２０】
発明の他の実施形態によれば、前記複数の血管は、右冠状動脈（RCA）、左前下行枝（LAD）、および左回旋枝（LCX）であり得る。
【００２１】
複数の血管が冠状動脈の場合に、前記の各冠状動脈に流れる血流量の割合を、以下の比例式によって決定することができる。
【式１】
【００２２】
2018187406.tif 000003
ここでは、
QLAD：左前下行枝の血流量、lLAD：左前下行枝の長さ
QLCX：左回旋枝の血流量、lLCX：左回旋枝の長さ
QRCA：右冠状動脈の血流量、lRCA：右冠状動脈の長さ
（lRCA）RV：右冠状動脈の右心室の血液供給の部分の長さ
（lRCA）LV：右冠状動脈の左心室の血液供給の部分の長さ
α：右心室の血液供給血管の補正係数
【００２３】
それぞれの冠状動脈の長さの割合を持って、それぞれの冠状動脈に流れる血流量の割合を求めると、心血管情報の冠血流予備量比（ＦＦＲ）を決定するために活用することができる。従来は、それぞれの冠状動脈が血液を供給する心筋ボリュームや質量に集中パラメータモデルを適用するための冠状動脈の毛細血管の抵抗値を求めたが、それぞれの冠状動脈の長さを持って簡単に集中パラメータモデルに適用するのための冠状動脈の毛細血管の抵抗値を求めることができる。
【００２４】
本発明の他の実施形態によれば、前記複数の血管は、脳や小脳に血液を供給する動脈であることがある。例えば、複数の血管は、頸部動脈、椎骨動脈、外頸動脈、中大脳動脈等であることができる。
【発明の効果】
【００２５】
本発明に係る心血管情報の中、冠状動脈の冠血流予備量比（ＦＦＲ）を決定する方法によれば、計算流体血流力学モデルと集中パラメータモデルのシミュレーション計算量を削減し、従ってＦＦＲ決定のための計算時間を短縮することができる。本発明に係る方法は、心血管情報を決定する方法は、単純化された冠動脈循環モデルを使用して冠状動脈内の血流に対して局部的にシミュレートするためである。特に、身体全体に対して閉回路を構成する従来の集中パラメータモデルを使用した場合に発生する患者一人一人のためのモデルの不確実性を減らすことができる。また、本発明に係る冠動脈の冠血流予備量比（ＦＦＲ）を決定する方法は、冠状動脈の形状についての簡単な境界条件を提示して、従来の方法と比較して計算量と計算時間において、より効率的である。また本発明に係る方法は電算流体力学モデルに大動脈を含まないため、従来の方法によるモデルと比べより単純なモデルである。
【００２６】
本発明に係る血管情報のうち、複数の血管に流れる血流量の割合を決定する方法が提供される。本発明に係る方法によると、動脈から分枝した枝血管が血液を供給する対応する身体の部分のボリュームや質量を求めずに、それぞれの枝血管に流れる血流量の割合を決定することができる。したがって、それぞれの血管が血液を供給する身体の部分のボリュームや体積を求めるために、対応する身体の部分の３次元モデルを求めるための複雑な計算を省略することができる。また、血管の３次元モデルだけを求めるための演算をして、血管に流れる血流量の割合を求めることができるので、体の部分の３次元モデルを求める過程で発生する不確実性を減らすことができる。
【図面の簡単な説明】
【００２７】
【図１】は、本発明に係る簡略化された冠循環モデルの模式図である。
【図２】は、本発明に係る計算流体力学モデルと簡略化された集中パラメータモデルの結合関係を示す模式図である。
【図３】は、本発明において、心拍数65 beats / minと心臓の収縮期/拡張期血圧120/80 mmHg状態での仮想的なPaoとPLVの例を示すグラフである。
【図４】は、本発明において、狭窄部位の直径比（Astenosis/A0）が0.5であり、直径（D）が5mmである仮想血管を示した図である。
【図５】は、図４に図示された仮想血管の狭窄部位の後の部分での軸方向速度を示すグラフである。
【図６】は、本発明において、狭窄がない仮想血管と狭窄された仮想血管に対する３次元幾何学的モデルである。
【図７】は、様々なモデルの冠状動脈の出口流量を比較して示すグラフである。
【図８】は、本発明において、50％狭窄がある仮想血管と狭窄がない仮想血管のＦＦＲ値を等高線図に示したものである。
【図９】は、本発明において、患者の冠状動脈の３次元形状モデルと簡略化された集中パラメータモデルを図で表したものである。
【図１０】は、本発明において、大動脈を含む冠状動脈モデルの模式図と集中パラメータモデルと大動脈の流動パターンを示したものである。
【図１１】は、本発明において、大動脈が除外されたモデルと大動脈を含むモデルに対して計算されたＦＦＲ値を測定された臨床データと比較した図である。
【図１２】は、本発明において、左前下行枝の最大流量でWSS分布を示す図である。
【図１３】は、心臓の血管の分布を示す図である。
【図１４】は、冠状動脈の主血管から分枝する微小血管の模式図である。
【図１５】は、直径が大きい血管から分枝する微小血管の抵抗を表示する概略図である。
【図１６】は、本発明に係る血管の長さを取得する方法を説明するフローチャートである。
【図１７】は、本発明によって得られた血流量と、従来の方法によって求めた血流量のブランドアルトマンプロット(Bland-Altman plot)である。
【図１８】は、図１７に図示されたグラフの相関分析の結果を示す図である。
【発明を実施するための形態】
【００２８】
本発明のその他の目的、特定な長点と新規の特徴は、添付された図面と関連した以下の詳細な説明と好ましい実施例においてからさらに明らかになる。以下では、本発明に係る冠動脈の心筋冠血流予備量比を計算するための新しい患者オーダーメイドモデルと心筋冠血流予備量比決定方法を開示する。
【００２９】
本発明に係る方法は、計算流体力学モデルのために冠状動脈だけを使用して、冠動脈血管系の集中パラメータモデルを含んでいる。本発明に係る方法は、冠状動脈の形状についての簡単な境界条件を提供し、従来の方法に比べて計算量を削減し、計算時間を短縮することができる。また、この方法は、モデルに大動脈を含まないため、従来の方法に比べて、より単純なモデルである。本発明に係る方法の方法の妥当性を検証するために、集中パラメータモデルと結合された計算流体力学モデルを用いて、３次元直線血管に対してシミュレーションした。本発明に係る方法を用いる計算結果を集中パラメータ方法だけを利用してシミュレーションした結果と比較した。また、臨床から得られた実質的な患者の形状について、本発明に係る方法を適用した。冠状動脈の狭窄がある患者のＣＴ画像データを利用して、患者のカスタムモデルを構成し、本発明に係る方法を適用して計算されたＦＦＲ値を臨床的に測定された結果と比較した。また、患者のＣＴ画像データから大動脈を含むモデルに対してＦＦＲ値を計算し、これを大動脈が含まれていないモデルについて計算されたＦＦＲ値と比較した。
【００３０】
狭窄された冠状動脈へのCFDモデル
冠状動脈の局所的計算流体力学モデルに冠動脈システムのシミュレーションを連結する統合的な接近法を使用した。冠循環シミュレーションのための集中パラメータモデルは、微小冠動脈から右心房まで冠状血管が表現されたものを使用し、冠状動脈内の局所の血流をシミュレートするために有限要素法によって実装された計算流体力学モデルを使用した。
【００３１】
そして、狭窄された冠状動脈の局所的計算流体力学モデルは、、狭窄された冠状動脈を貫通する血流の３次元ナビエ-ストークス方程式(Navier-Stokes equations)を計算するためにPISO(pressure implicit splitting of operators)方式を用いた有限要素解析方法を使用した。詳細な解析アルゴリズムは、研究論文（Shim and Chang、1994、1997）に開示されている。
【００３２】
血流は非圧縮性、層流、粘性流体と仮定し、ニュートン流体で処理した。血液の密度は1,060kg / m3と想定した。血液の粘性係数μは、患者のヘマトクリット値とアインシュタインの方法（Einstein、1906）を利用して獲得した。
【式２】
【００３３】
2018187406.tif 000004
ここでは、「HCT」は、患者のヘマトクリット値であり、μ0は水とほぼ同じプラズマの粘度を示す（μ0= 0.0015kg/ms）。
【００３４】
冠動脈循環の集中パラメータモデル
また、冠循環の集中パラメータモデルは、冠状血管に限定した。冠循環は、図１に示すように、三種類の構成要素、すなわち、冠状動脈（coa）、冠状動脈の毛細血管（coc）、冠状静脈（cov）を含んでいる。冠循環の集中パラメータモデルは、それぞれの抵抗と静電容量などの要素で構成され、電気的アナログ観点から数式化される。血液の循環に関連する血流力学的パラメータを使用して、三つの部分の構成要素それぞれについて、圧力、体積、血流を計算する。他の毛細血管系とは対照的に、心臓の収縮期中、血管周囲の左心室の圧力が非常に強くなるため、冠状動脈の血流量は低くなる。拡張期中の心臓の筋肉が緩みながら左心室の毛細血管に流れる血液が妨げられていないので、冠状動脈に流れる血流の移動はより自由になる。この機序をモデル化するために、冠状動脈の心筋内部圧力（図１のPimp）を使用した。この心筋内部の圧力は、心室の圧力と相関関係を有する。
【式３】
【００３５】
2018187406.tif 000005
この式でPLVは、左心室の圧力であり、γは0.75を使用した。このγの値は、シュライナー（Schreiner）などにより提案された。以下で説明するように、冠状血管の毛細血管における区画した周辺の抵抗（Rcoc、Rcov）を考慮にいれ体積依存形の値を使用した。
【式４】
【００３６】
2018187406.tif 000006
RactualはRcocまたはRcovの実際の値であり、Roriginalは Rcoc または Rcov のデフォルト値である。式４においてβは経験的定数であり、Vcocは冠状毛細血管の体積を示す。βはシュライナー等により提示された0.3に設定した。式４に示すように、冠状血管内の血流の流れは、圧力勾配の信号に応じて、巡行（ie、Pcoa> Pcoc）、または逆行することができる。しかし、毛細血管の体積が0に近づくようになれば、もはや押し出されてくることがないので、逆流が中断される。さらに、毛細血管が圧縮されると、抵抗が増加して、血流の動きが調節される。式４の逆圧力勾配に応じて毛細血管の体積が0に近づくと逆流も0に近く減少する。冠循環の集中パラメータ方法を扱う詳細については、既存の研究（Schreiner et al。、1990; Shim et al。、2000）に記載されている。
【００３７】
計算流体力学モデルと集中パラメータモデルの組み合わせ
また、集中パラメータ方法と計算流体力学モデルの統合は、冠循環の集中パラメータ方法と計算流体力学モデルを結合するために、両方のシステムが相互の影響を与えるという家をした。図２は、計算流体力学モデルと集中パラメータモデルの結合を概略的示す。計算流体力学モデルの出口で計算された流量は、集中パラメータモデルに伝達される。常微分方程式を介して集中パラメータモデルの三つの部分の圧力と流量を計算する（Lim et al。、2009）。LPMで計算された冠動脈圧の値は、次のステップの計算のために計算流体力学モデルの出口境界条件としてPoutletを計算するために使用されるように計算流体力学モデルに伝達される。このように結合された方法は、、微小血管と結合された大きな血管の局所的流体力学の解析に効果的であり、これはBrown等の論文（Brown et al。、2012）に詳細に説明されている。
【００３８】
Ｋｉｍ等により提示された従来の方法とは異なり、本発明に係る方法は、大動脈血圧（図２でのPao）の瞬間的なパターンを計算流体力学モデルの入口境界条件として使用した。これらの圧力パターンは、心臓のサイクルと患者の心臓収縮期/拡張期血圧を合わせるために提案された。左心室の時間変化に伴う静電容量モデルを大動脈血圧の規定された一時的なパターンと結合することにより、PLVを入手し、これは式２のPimpの計算のために使用される。図３は、心拍数65 beats / minと心臓収縮期/拡張期血圧120/80 mmHg状態での仮想PaoとPLVの例を示している。
【００３９】
簡略化された冠状動脈の３次元モデル
本発明に係る方法の検証のために、図４に示すように、少しの狭窄があるチューブ（仮想血管）をシミュレートして得られた結果と、既に数値的研究を介して出てきた結果と比較した。狭窄部位の割合であるAstenosis / A0は0.5であり、直径（D）は5mmである。液体の密度と粘性は、それぞれ755 kg / m3、0.00143 Pa・secであり、既存の数値研究に示されたサインの速度を入口境界条件に適用した。
【式５】
【００４０】
2018187406.tif 000007
式５に示した入口速度Vzは軸方向であり、Dは、チューブの直径である。rとtは、それぞれ半径と時間である。図５には、狭窄部位の後の部分での軸方向速度を示す。サイン（正弦）周期入口速度は、Tpの式で表現した。計算された速度波形は、従来の研究結果と非常によく一致した。これは、本発明に係る計算流体力学の方法の時間的、空間的精度を示す。
【００４１】
冠状動脈の患者  オーダーメイドモデル
また、臨床的な実際の形状を利用して、シミュレーションモデルをテストするために、ＣＴ画像データから冠状動脈の３次元形状モデルを生成して、患者オーダーメイドモデルを開発した。臨床データの出所は、大韓民国の蔚山大学病院で、蔚山大学病院のデータの使用は、IRBの承認を受けた。2014年の病院を訪れた患者のうち、直径2.5mm以上の主要な冠状動脈の50％以上の狭窄を示した患者一人を選択した。
【００４２】
ＣＴ画像のセグメンテーションのための活動輪郭アルゴリズムに基づいたセグメンテーションプログラムを開発した。また、冠動脈形状の中心線抽出は、Antica等により提示されたアルゴリズムを用いた。また、HypermeshTM（http://www.altairhyperworks.com）ソフトウェアを使用して、ＣＴ画像データを整えて３Ｄ解析のための格子システムを作成した。開発された格子システムにドロネー（Delaunay）のアルゴリズムを利用して３次元ボリュームメッシュを生成した。
【００４３】
冠状動脈の簡素化された３次元モデルの心血管情報を決定する方法を確認するために、左前下行枝を模倣したストレート血管モデルの50％狭窄があるモデルと狭窄がないモデルの血行動態シミュレーション結果を集中パラメータモデルの結果と比較した。図６には、３次元の幾何学的モデルに狭窄がない血管と狭窄された血管に対して、それぞれ59,353個と121,292個のグリッドが示されている。図７は、このように構成された２つのモデルの出口流量と集中パラメータのモデルでの出口流量を計算して比較したものである。通常の状態での左前下行枝（LAD）の平均流量は1.3cc / sと仮定した。最大充血状態の場合、平均流量を定常状態流量の４倍5.2cc / sに増加した。図７に示すように、冠状動脈の流れは、心室収縮期中に減少し、拡張期に増加した。これらの結果は、文献に報告された冠状動脈の血流速度の波形と一致する。狭窄がないモデルと集中パラメータモデルは、血流のパターンが似ているが、拡張期での50％狭窄モデルは、通常のモデルに比べて流量は減少した。
【００４４】
図８は、50％狭窄があるモデルと狭窄がないモデルのＦＦＲ値を輪郭図に示したものである。狭窄がないモデルは、ＦＦＲ値がほぼすべての部分で約1.0の値が出てきた、狭窄モデルである場合には、狭窄部位の後ろ測定したＦＦＲ値（約0.8）が減少したことが確認できた。興味深いことに、狭窄があるモデルのＦＦＲ最小値は、狭窄病変のくびれ部位で観察された。
【００４５】
患者オーダーメイド型モデルの観点から、本発明に係る方法を検証するために、左回旋枝（LCX）で狭窄が発生した患者の３次元血流力学をシミュレートした。このモデルの３次元形状の構造と境界条件は、前に説明された。図９は、患者のＣＴ画像データとして生成された計算流体力学モデルの格子システムと集中パラメータモデルを示す。下の表１は、シミュレーションのために必要な患者の生理学的なデータを示する。計算のために、最初に心拍出量と心拍数によって測定される休眠時の心拍出量を求めた。冠状動脈の総流量は、心拍出量の4％に仮定された。最大充血状態のために、冠状動脈の流量を休眠状態の4倍に増加した。それぞれの冠状動脈の血管抵抗値は、冠状動脈の総流量と各種類の断面積で計算された。静電容量値も同様に、従来の方法に応じて決定した。下の表２は、集中パラメータモデルの左前下行枝（LAD）、左回旋枝（LCX）、右冠状動脈（RCA）の総抵抗値と静電容量値を示す。
【００４６】
冠状動脈の入口で、収縮/拡張期血圧で測定された患者の心拍数によって定義された大動脈圧の規定されたパターンは、図３の通りである。
【００４７】
【表１】
2018187406.tif 000008
【００４８】
【表２】
2018187406.tif 000009
【００４９】
また、シミュレーションされたＦＦＲ値の大動脈部分の影響を評価するために、大動脈が含まれるモデルに対してシミュレーションを行った。大動脈が含まれるモデルでは、大動脈のないモデルで冠状動脈の入口境界条件を指定することに備えて、大動脈の入・出口境界条件が指定されなければならない。左心室から大動脈の流動パターンは、図１０に示されたようにQmaxの最大流量を有する正弦曲線であると仮定した。測定された拍出量と心周期に基づいてQmaxを計算することができる。入口流量を使用して入口面からの速度パターンを規定し、これは、図１０に示されている。入口面で空間速度分布は均一なものと仮定される。Ｋｉｍらの研究では、大動脈の出口は全体の心血管系の集中パラメータモデルの構成要素に連結されている。本発明に係る方法では、出口境界条件で大動脈圧の特定のパターンを使用した。
【００５０】
下の表３は、格子の数と心周期に応じた計算時間を示す。大動脈のないモデルは、大動脈があるモデルよりも計算がより効率的である。図１１には、両方のモデルの計算されたＦＦＲ値と測定された臨床データの値が図示されている。臨床ＦＦＲはプレッシャーガイドワイヤー（pressure guide-wire Certus、St. Jude Medical Systems、Uppsala、Sweden）によって測定された。左前下行枝（LAD）、左回旋枝（LCX）、右冠状動脈（RCA）で測定されたＦＦＲ値は、３つのケースについてはほとんど類似する。また、大動脈が含まれていないモデルのＦＦＲ輪郭パターンが大動脈が含まれているモデルのパターンとほぼ同じであった。左前下行枝（LAD）から最大流量のときに計算された壁面せん断応力と速度ベクトルは、図１２に示されている。最も高いせん断応力の部分は、左回旋枝（LCX）狭窄の近くにあった。
【００５１】
【表３】
2018187406.tif 000010
【００５２】
ＦＦＲシミュレーションは、非侵襲的であり、安全な方法であるため、冠動脈狭窄の機能の程度を評価するのに非常に有用である。ＦＦＲシミュレーションの従来の方法は、計算流体力学モデルでの大動脈と全体の心血管システムの集中パラメータ方法を含でおり、これは多くの計算量と計算時間を必要とする。しかし、臨床適用が効率的であるためのＦＦＲシミュレーション方法は、個人のコンピューターでも実行することができるほど簡単で速くなければならない。本発明に係る方法は、狭窄した冠動脈のＦＦＲ値を計算するため、新しく効果的な方法を提供する。
【００５３】
本発明に係る方法は、従来の方法といくつかの明確な違いを持っている。第一に、本発明に係る方法に適用されるモデルは、計算流体力学モデルで大動脈を含んでいない。したがって、それぞれ左心室と体循環システムに連結されている大動脈の入口、出口の境界条件を必要としない。第二に、本発明に係る方法に適用されるモデルは、集中パラメータモデルが閉回路ではない。集中パラメータシステムにおいて、閉回路モデルは、動脈と静脈、肺血管、左心、右心がそれぞれの集中パラメータとして含まれていなければならない。また、冠動脈循環システムのための集中パラメータの構成要素もある。しかし、本発明に係る方法、したがって適用されるモデルの集中パラメータの構成は、唯一冠動脈微血管系（冠状動脈、冠状毛細血管、冠状静脈）のみを含んでいる。計算流体力学モデルで大動脈と集中パラメータ方法の閉回路が使用されている従来のモデルは、いくつかのパラメータの（抵抗値及び静電容量値：体動脈、体静脈、肺血管、静脈、左心室、右心室）を必要とする。
【００５４】
これらの各パラメータは、基準ケースで代表値を持つが、それを患者のカスタムケースで適用する場合、患者オーダーメイドモデルとするのは難しい。本発明に係る方法では、大動脈を含まず、規定された圧力パターンを冠状動脈の計算流体力学の境界条件に適用する。この規定された圧力パターンは、図３のように、臨床で測定された血圧を介して抽出する。したがって、閉回路モデルのための集中パラメータモデルの各変数は、本発明に係る方法を適用するためのモデルには必要とせず、したがって、モデルの不確実性が低減される。
【００５５】
計算を考慮したとき、本発明に係る方法は、いくつかの長点がある。表３に示すように、大動脈を含んでいないモデルの３次元格子の数は、大動脈を含むモデルに比べて少なく、これは私たちのモデルの計算的効率に寄与する。また、従来の方法は、大動脈の計算流体力学モデルを全体循環集中パラメータモデルと結合することをベースにしており、これは、各時間帯のステップで二つのモデルとの間の反復的な相互作用を必要とする。本発明に係る方法では、大動脈の計算流体力学モデルと全体の心血管の集中パラメータモデルを冠状動脈の計算流体力学モデルの規定された境界条件に置き換えし、これは大動脈計算流体力学モデルと集中パラメータモデルとの間の繰り返しされた計算を必要としない。計算の精度の面で、本発明に係る方法は、従来の方法と比較して非常に合理的な結果を示す。直線状の血管から計算流体力学モデルの出口の流量は、集中パラメータモデルで計算された値とほぼ同じである。一方、50％狭窄がある場合、流量とＦＦＲは、通常の場合と比較して減少した。図７に示すように、50％狭窄を介して計算された平均流量は実験的観測と比較して類似しており、これは本発明に係る方法の計算精度を示す。
【００５６】
３Ｄ患者オーダーメイドモデルで冠状動脈の形状は、患者のＣＴデータから抽出した。また、冠状動脈の全抵抗と静電容量を定め、それたちを全体の冠動脈流量と冠状動脈の分枝の断面積の割合に応じて、それぞれの冠状動脈の抵抗と静電容量を再分配した。図１１に示すように、右冠状動脈、左回旋枝、左前下行枝でのＦＦＲ値の計算結果は、臨床的に測定された結果と一致する。興味深いことに、大動脈のモデルの計算結果は、大動脈のないモデルの結果とほぼ同じであり、これは大動脈の部分は、冠状動脈ＦＦＲの計算値に影響を与えないことを示す。冠状動脈の高い壁面せん断応力はプラークの成長やプラークの破裂の可能性が高いことを意味する。したがって、壁面せん断応力分布は、狭窄冠動脈の重要な変数である。図１２は、左前下行枝の最大流量でWSS分配を見せてくれて、これは狭窄病変は、比較的高い値を持つことを示す。
【００５７】
動脈からの分枝された複数の血管に流れる血流量の割合を定める方法
冠状動脈や脳血管の血流解析のためには、動脈からの枝血管への血流の分配比率を求める必要がある。冠状動脈の分枝に流れる血流量の割合は、CT-ＦＦＲを計算する際、各冠動脈分枝に連結されたLPMモデルの抵抗値を設定する基準となる。ところが、複数の枝血管への血流分配は、患者に応じてそれぞれ異なり、そのための患者オーダーメイドモデルの方法が必要である。
【００５８】
以下では、本発明の他の側面に沿って、患者に応じて、冠動脈分枝の血流比の設定方法を提供して、ＣＴ−ＦＦＲ分析などのためのシミュレーション技術の精度と利便性を高めようとする。
【００５９】
冠状動脈への全体流量（flow rate）は、大きく見てLAD、LCX、そしてRCAの３つに分けられる。冠状動脈全体に流れる流量（Qtotal）は、患者ごとに異なる割合でそれぞれの分枝に分けられる。もし患者の冠状動脈に流れる血流量の割合を求めることができれば、患者の心臓の末梢血管による抵抗を設定することができる。患者の冠状動脈に流れる血流量は、以下の式６のように表すことができる。
【式６】
【００６０】
2018187406.tif 000011
【００６１】
式５を利用して、冠状動脈に流れる血流量と血圧と抵抗の関係を下の式７のように表すことができる。
【式７】
【００６２】
2018187406.tif 000012
【００６３】
式７で、ΔPは冠状動脈が大動脈から分離されているところでの圧力と心臓の筋肉に流れる血液が、最終的に戻っていく場所である心房圧力との差を意味する。冠状動脈の心臓分布構造は、図１３に示すように、直径が大きい動脈に沿って末端に行きながら、直径が小さい動脈が分枝する形態で構成されている。図１３の心血管の写真はhttp://www.plastinate.com/leistungen/ausgusspl.htmから援用したものである。図１３に示すように分布する血管の血流抵抗を見てみる。ハーゲン・ポアズイユ流れ（Poiseuille flow）の仮定によれば、血管の抵抗は、直径の４乗に反比例するので、相対的に直径が大きい血管ではなく、小さな血管においての抵抗がはるかに大きい。したがって、直径が大きい冠状動脈での抵抗は無視することができる。なので、直径が小さい血管の抵抗のみを考慮してモデル化しても、実際の抵抗と大きな差がないだろう。
【００６４】
図１４は、図１３に図示された冠状動脈の血管分布を模式的に示したものである。図１４を参照すると、直径が大きい冠状動脈から分枝された直径の小さい血管は並列的にn個連結されていることを前提する。この場合、直径が小さい血管の抵抗だけを考慮すると、冠状動脈および枝血管の総抵抗は、以下の式８のように表すことができる。
【式８】
【００６５】
2018187406.tif 000013
【００６６】
式８でそれぞれの直径が小さい血管の抵抗が同じであると仮定すると、以下の式９と抵抗Ｒを求めることができる。
【式９】
【００６７】
2018187406.tif 000014
【００６８】
図１３を参照すると、冠状動脈には分枝から末端に沿ってほぼ一定の間隔で直径が小さい血管が分枝されているので、直径が小さい血管の分枝数nは、直径が大きい冠状動脈の長さに比例する。したがって、冠状動脈の抵抗は、次の式１０のように表すことができる。
【式１０】
【００６９】
2018187406.tif 000015
ここでは、kは任意の比例定数であり、lは冠状動脈の血管の長さである。比例定数kと直径が小さい血管の抵抗R1は、左心室、心臓の筋肉に血液を供給するLADとLCXの場合、ほぼ同じであると仮定することができる。実際に冠状動脈の左心室の部位の微小血管の解剖学的なパターンは、どこでもほぼ同じであると知られている。
【００７０】
式１０を式７の左心室に血液を供給する冠状動脈に対して適用すると、下記の式１１が求められる。
【式１１】
【００７１】
2018187406.tif 000016
【００７２】
しかし、RCAの場合は少し複雑である。RCA分枝の場合、右心室に血液を供給する血管と左心室に血液を供給する血管が共存する。したがって、RCAの場合、右心室に血液を供給する血管と左心室に血液を供給する血管が連結されていることを前提として、次のような式１２で表すことができる。
【式１２】
【００７３】
2018187406.tif 000017
【００７４】
式１２で、左心室（式ではLVと表示）に血液を供給する血管の場合、nの値や比例定数kの値はLAD、LCXと同じ値を使用することができる。一方、右心室に血液を供給する血管に対しては、同一なnの値や比例定数kの値を使用することができない。左心室の場合、人体全体に血液をポンピングするために大きな圧力を発生させるが、右心室の場合には、肺だけに血液を供給すればよいため、比較的小さな圧力を発生させる。実際に左心室内部での血流の圧力は、約100 mmHgに達するが、右心室の場合、左心室の圧力の３分の１から４分の１の間の値を持つ。したがって、右心室の筋層は、左心室に比べて相対的に薄く、分枝された血管の粗密度も低下する。したがって、RCA血管の右心室に血液を供給する部分は、単位動脈の長さ当、分枝されている直７３８径が小さい血管の数が小さくなる。したがって右心室に血液を供給する血管の部分については、左心室に適用したnの値を１よりも大きい適切な定数（α）に分けて、以下の式１３のようにRCAに対する抵抗を求めることができる。
【式１３】
【００７５】
2018187406.tif 000018
【００７６】
式１３を式７のRCAに対する式に代入すると、以下の式１４のようにRCAに流れる血流量を求めることができる。
【式１４】
【００７７】
2018187406.tif 000019
【００７８】
前記の式１０と式１４で、ΔPは大動脈の圧力と心臓の筋肉に流れた血液が、最終的に戻っていく場所である心房圧力との差を意味するので、冠状動脈のLAD、LCX、RCAの両方に対して同じである。したがって、冠状動脈のそれぞれに流れる流量の比は次の式１５のように表すことができる。
【式１５】
【００７９】
2018187406.tif 000020
【００８０】
LAD、LCX、RCAの血管の長さを求めると、式１５から各分枝に流れる血流の比を知ることができる。そして、それぞれの分枝に流れる血流量の割合通りに冠状動脈に流れる血流量（全体心拍出量の約４％）を配布すると、それぞれの冠状動脈に流れる血流量を推定することができる。また、それぞれの冠状動脈に対して推定された血流量に基づいて、それぞれの冠状動脈への全体的な血管の抵抗（LPMでの抵抗）の値を設定することができる。
【００８１】
冠状動脈の３次元形状モデルから血管の長さを取得する方法
本発明に係る一実施形態に、冠状動脈の３次元形状モデルから血管の長さを求め、それぞれの血管に流れる血流量の割合を定める方法について述べる。本発明に係る方法によれば、冠状動脈のLAD、LCX、RCAがそれぞれの血液を供給する心筋の体積を求めず、それぞれの冠状動脈の分枝の長さを求め、それぞれの冠状動脈の分枝に流れる血流量の割合を推定することができる。したがってＦＦＲのような患者のカスタマイズ心血管情報を決定するために、左心室の心臓の筋肉の体積を求めるためのＣＴ画像データの画像解析手順が必要ないようになる。
【００８２】
本発明に係る方法を用いる場合、単にCFD（Computational Fluid Dynamics）モデルの実装のためにＣＴ画像データから抽出された冠状動脈の３次元モデルを利用して、それぞれの冠状動脈の長さだけを計算して血流を計算する。したがって、冠状動脈の３次元モデル解析のための計算量を削減し、計算時間を短縮することになる。
【００８３】
図１６を参照すると、次のような順序で冠状動脈の３次元モデルから冠状動脈の長さを求める。
１）ＣＴ画像データを処理して、冠状動脈の３次元形状モデル生成する。
２）生成された冠状動脈の３次元形状モデルから冠状動脈の中心線を作成する。
３）生成された中心線の末端の面積を計算して、一定面積以下となるポイントを末端地点に選定する。
４）冠状動脈の分枝から末端地点までの中心線の長さを計算する。
次に、冠状動脈がRCAの場合には、次のような順序で冠状動脈の長さを求める。
５）RCAの血管の長さを測定するが、PLA（postero-lateral artery）、PDA（posterior descending artery）に分かれる分枝までの長さと、分枝点以降の長さを計算する。
６）PLAとPDAの分枝になる前までの長さは、一定の値（1/3 - 1/4の間の値）に分けて長さを補正してくれる。
７）PLAとPDAの長さは、左心室に血液を供給するので、LAD、LCXと同じ長さに計算された長さを使用する。
【００８４】
冠動脈血管の長さを利用した冠状動脈の血管への血流の計算
本発明に係る方法の妥当性の検証のために、心筋の体積を利用して、冠状動脈の血流量を推定した結果と、本発明に係る方法に応じて、冠状動脈に流れる血流量を計算した結果とを比較した。
【００８５】
検証には、大韓民国の蔚山大学病院で提供されたＣＴ画像データを使用した。それぞれの冠状動脈の分枝が血液を供給する心臓の筋肉の量を計算するために市販のimage processing softwareであるTerarecon社Aquariusを使用した。蔚山大学校病院から提供されたＣＴ画像データ（25人の患者-50個の血管）に対してそれぞれの冠状動脈の血管の長さを計算して、計算された血管の長さをもとに、それぞれの冠状動脈に流れる血流量を推定した。表４には、血管の長さをもとに推定した冠状動脈に流れる血流量が記載されている。表４で、心臓の筋肉のボリュームをもとに推定した冠状動脈に流れる血流量はFlowvolumeで表示し、本発明により、血管の長さをもとに推定した冠状動脈に流れる血流量はFlowlengthで表示した。
【００８６】
【表４】
2018187406.tif 000021
【００８７】
また、LCAとLCXについて求めた血流量について、Bland-Altman plotを描いて図１７に示した。図１８は、図１７に図示されたBland-Altman plotの相関分析の結果である。図１８を参照すると、本発明に係る方法で求めた血流量と、従来の方法で求めた血流量は相関係数（correlation coefficient）が0.826で相関関係が非常に高いことが分かる。
【００８８】
Bland-Altman plotは二つの測定方法で求められた値の差があるかを分析するためのグラフである。図１７の中央線は相加平均（mean）値を示し、この値が０に近いほど二つの値の差の誤差が少ないことを意味する。SD値も０に近いほど、両方の値が似ていることを示す。
【００８９】
前記で説明され、図面に示された本発明の実施例は、本発明の技術的思想を限定するものと解釈してはならない。本発明の保護範囲は、特許請求の範囲に記載された事項によってのみ制限され、本発明の技術分野における通常の知識を有する者は、本発明の技術的思想を様々な形で改良変更することが可能である。したがって、これらの改良および変更は、通常の知識を有する者に自明なものである限り、本発明の保護範囲に属するようになるだろう。
【符号の説明】
【００９０】
QLAD：左前下行枝の血流
lLAD：左前下行枝の長さ
QLCX：左回旋枝の血流
lLCX：左回旋枝の長さ
QRCA：右冠状動脈の血流
lRCA：右冠状動脈の長さ
(lRCA)RV：右冠状動脈の右心室の血液供給の部分の長さ
(lRCA)LV：右冠状動脈の左心室の血液供給の部分の長さ
α：右心室の血液供給血管の補正係数
【図１】
2018187406.tif 000023
【図２】
2018187406.tif 000024
【図３】
2018187406.tif 000025
【図４】
2018187406.tif 000026
【図５】
2018187406.tif 000027
【図６】
2018187406.tif 000028
【図７】
2018187406.tif 000029
【図８】
2018187406.tif 000030
【図９】
2018187406.tif 000031
【図１０】
2018187406.tif 000032
【図１１】
2018187406.tif 000033
【図１２】
2018187406.tif 000034
【図１３】
2018187406.tif 000035
【図１４】
2018187406.tif 000036
【図１５】
2018187406.tif 000037
【図１６】
2018187406.tif 000038
【図１７】
2018187406.tif 000039
【図１８】
2018187406.tif 000040