(57)【要約】【課題】脳梗塞後の血流再建手術について、リモデリングを考慮した脳血管系の過渡血流量を評価し、手術後の出血等のリスクを予測することを可能とすること。【解決手段】脳血管系の複数の系統のそれぞれの血管梗塞発症前の健常時での血液の流量または流動抵抗、血管梗塞発症時においてクロスフローが発生する複数の系統間の接続関係、接続関係にある系統間における流動抵抗、およびクロスフローによって脳血管系におけるリモデリングが達成されるまでの時間を表す時定数の入力を受け付け、入力された各系統の健常時での流量または流動抵抗、複数の系統間の接続関係、系統間における流動抵抗、および時定数を用いて、クロスフローによる血管梗塞発症後の各経過時点での各系統における流量を算出する。【選択図】図１
2015018320.tif 000002
【特許請求の範囲】
【請求項１】
血管梗塞発症後の血流量の予測を行うためにコンピュータ装置が実行する血流量予測方法であって、
脳血管系の複数の系統のそれぞれの血管梗塞発症前の健常時での血液の流量または流動抵抗、前記血管梗塞発症時においてクロスフローが発生する前記複数の系統間の接続関係、前記接続関係にある系統間における流動抵抗、および前記クロスフローによって前記脳血管系におけるリモデリングが達成されるまでの時間を表す時定数の入力を受け付ける入力ステップと、
前記入力ステップにて入力された前記各系統の健常時での流量または流動抵抗、前記接続関係、前記接続関係にある系統間における流動抵抗、および前記時定数を用いて、前記クロスフローによる前記血管梗塞発症後の各単位時間の経過時点での前記各系統における流量を算出する算出ステップと、
を含む血流量予測方法。
【請求項２】
前記入力ステップにて入力された前記健常時での流量または流動抵抗と、前記算出ステップにて算出された前記各系統における流量とから、前記各経過時点での前記各系統における流動抵抗を算出し、この算出された流動抵抗と前記各系統における流量との積を前記各経過時点での前記各系統の圧力損失とし、前記各系統の圧力損失の差から生じる系統間の圧力差によって、前記クロスフローが生じると仮定することを特徴とする請求項１に記載の血流量予測方法。
【請求項３】
前記算出ステップは、
前記健常時での前記各系統における流量の合計値を脳に対する必要流量とし、前記血管梗塞発症後の単位時間のある経過時点における前記各系統における流量の和が前記必要流量に漸近するように、一定の比率を前記ある経過時点での前記各系統への血液の供給量に掛けて次の経過時点の前記各系統の供給量を算出することで、前記血管梗塞発症後の各経過時点での前記各系統における流量を算出することを特徴とする請求項１に記載の血流量予測方法。
【請求項４】
前記算出ステップにて算出された前記各経過時点での前記各系統への供給量と、前記各系統における流量とから算出する、前記健常時での前記各系統における流量に対する前記血管梗塞発症後の血流減少比の時間積分を前記各系統の血管ダメージ量とし、前記健常時での前記各系統における流量に対する前記血管梗塞発症後に行われる血流再建手術後の血流増大比の時間積分を前記各系統の血管への負荷量とし、前記血管ダメージ量と前記血管への負荷量とを比較する再出血リスク評価ステップをさらに含む請求項１に記載の血流量予測方法。
【請求項５】
請求項１から４のいずれか一項に記載の方法をコンピュータ装置に実行させるコンピュータプログラム。
【請求項６】
請求項１から４のいずれか一項に記載の方法をコンピュータ装置に実行させるコンピュータプログラムを記憶した記憶媒体。
【請求項７】
プロセッサと、請求項１から４のいずれか一項に記載の方法を前記プロセッサに実行させるコンピュータプログラムを記憶した記憶媒体と、を備えるコンピュータ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脳梗塞発症患者の血流再建手術後の血流量の予測解析に関し、再出血リスクの評価等への応用が可能な予測解析に関する。
【背景技術】
【０００２】
脳血管系は、図８に示すような構造となっている。すなわち、総頸動脈と椎骨動脈を通じて心臓から血液が供給され、供給された血液は、前大脳動脈、中大脳動脈、後大脳動脈に分かれてそれぞれの灌流域の毛細血管網に配分される。これらの動脈の間の交通動脈は、Willis動脈輪を構成し血管異常発生時に、ある程度血流を融通しあうことができる。さらに下流のそれぞれの灌流域の毛細血管網の間にもクロスフローされうる。
【０００３】
ところで、脳動脈の梗塞や痙攣に起因する動脈の狭窄化は脳組織に深刻なダメージを与え、生命の危機や多くの障害をもたらす。一般にその治療として、ステントやバイパスによる血流再建手術が行われる。
【０００４】
動脈が狭窄化した時に下流側灌流域の虚血による組織ダメージを防止するため、流動抵抗を調整し血流量を再配分するリモデリングが知られている。リモデリングは外的負荷が生じるイベントに対して生体に変化が生じる適応現象であるが、負荷の緩和の代償に別の問題を引き起こすことも少なくない。そのメカニズムや要する時間や血流再配分の達成程度などは、十分明らかにはされているとは言えない。
【０００５】
現在、急速な技術進歩により、核磁気共鳴を用いた血管造影法（MRA）やさらに骨や組織を消去し血管のみを抽出するデジタル・サブトラクション・アンギオグラフィー（DSA）によって、血管系の三次元画像を得ることができる。ただし、比較的大きい血管系や、動脈瘤や狭窄のような病変の検出は可能であるが、血管壁や毛細血管のような微細構造までは測定できない。
【０００６】
動脈瘤や分岐などの三次元流路構造からの計算流体力学(CFD)を用いた血流挙動の解析は、非特許文献１において述べられているように広く行われている。また、非特許文献２に述べられているように、広い範囲の血管系を一次元ネットワークとして脈動を含めて計算する手法や、非特許文献３に述べられているように、体温調節における血流の熱エネルギー輸送の役割の数値解析技術も開発されている。
【先行技術文献】
【非特許文献】
【０００７】
【非特許文献１】大島まり、「脳血管障害における医用画像に基づく血流シミュレーション」、ながれ、一般社団法人日本流体力学会、2002年、第21巻、p.122-128
【非特許文献２】M. S. Olufsen,“Modeling the Arterial System with Reference to an Anesthesia Simulator（血流システムのモデリングと麻酔シミュレータへの適用）”, 1998年, Ph. D. Thesis, Roskilde University, Denmark
【非特許文献３】横山真太郎、「生体内熱移動現象」、北海道大学図書刊行会、1993年
【発明の概要】
【発明が解決しようとする課題】
【０００８】
脳血管系は繊細で複雑な構造を持つため、脳神経外科手術において、術後に血流量が変わることにより却って別の障害を発生させることがないよう配慮することは極めて重要である。従来の知見から脳神経外科医のためのガイドラインが定められているが、さらに判断材料を加えるために、血流シミュレーションによる障害発生のリスク評価が望まれている。しかしながら、上述したような従来の数値解析技術は、固定された（時間変化のない）血管系の条件下における血流にその対象が限られており、動脈狭窄化や血流再建手術等のイベントに対するリモデリングの過程を追う時間変化を取扱うことは困難である。よって、従来の数値解析技術は、術後の過灌流にともなって血管に過大な負荷が加わり障害が生じることのリスク評価には適していない。
【０００９】
また、生体反応を含む複雑な状態変化について純粋に数学モデルを構築することは、個々人の身体条件の相違あるいは確率的なふるまいのため、極めて困難である。メカニズムの知られていない多種多様な対象については、部分的なモデルの精密化よりはむしろ全体的に単純なモデルを用いつつ実例に合うようにパラメータを調整する手法が有用な手段となることがある。
【００１０】
本発明の目的は、上記課題を解決するべく、脳梗塞の発症と血流再建手術について、リモデリングを考慮しつつ脳血管系の過渡血流量を評価し、手術後の出血等のリスクを予測することが可能な予測方法を提供することである。
【課題を解決するための手段】
【００１１】
上記課題を解決するために、本発明の一態様は、血管梗塞発症後の血流量の予測を行うためにコンピュータ装置が実行する血流量予測方法であって、脳血管系の複数の系統のそれぞれの血管梗塞発症前の健常時での血液の流量または流動抵抗、前記血管梗塞発症時においてクロスフローが発生する前記複数の系統間の接続関係、前記接続関係にある系統間における流動抵抗、および前記クロスフローによって前記脳血管系におけるリモデリングが達成されるまでの時間を表す時定数の入力を受け付ける入力ステップと、前記入力ステップにて入力された前記各系統の健常時での流量または流動抵抗、前記接続関係、前記接続関係にある系統間における流動抵抗、および前記時定数を用いて、前記クロスフローによる前記血管梗塞発症後の各単位時間の経過時点での前記各系統における流量を算出する算出ステップと、を含む血流量予測方法である。
【００１２】
また、本発明の他の態様は、上記の方法をコンピュータ装置に実行させるコンピュータプログラムである。
【００１３】
また、本発明の他の態様は、上記の方法をコンピュータ装置に実行させるコンピュータプログラムを記憶した記憶媒体である。
【００１４】
また、本発明の他の態様は、プロセッサと、上記の方法を前記プロセッサに実行させるコンピュータプログラムを記憶した記憶媒体と、を備えるコンピュータ装置である。
【発明の効果】
【００１５】
本発明によれば、脳梗塞の発症と血流再建手術について、リモデリングを考慮しつつ脳血管系の過渡血流量を評価し、手術後の出血等のリスクを予測することが可能である。
【図面の簡単な説明】
【００１６】
【図１】本発明の一実施形態に係る血流量予測プログラムにおける計算の流れの一例を示すフロー図である。
【図２】本実施形態に係る血流量予測プログラムによる血流再建手術後の過渡血流分布の計算例1を示す図である。
【図３】本実施形態に係る血流量予測プログラムによる血流再建手術後の過渡血流分布の計算例２を示す図である。
【図４】本実施形態に係る血流量予測プログラムによる血流再建手術後の過渡血流分布の計算例３を示す図である。
【図５】本実施形態に係る血流量予測プログラムによる血流再建手術後の過渡血流分布の計算例４を示す図である。
【図６】本実施形態に係る血流量予測プログラムによる血流再建手術後の過渡流路抵抗分布の計算例４を示す図である。
【図７】本発明の一の実施形態に係る血流量予測プログラムを実行するコンピュータ装置のハードウェア構成の一例を示す図である。
【図８】脳動脈系統の概略を示す図である。
【発明を実施するための形態】
【００１７】
以下、図面を参照して、本発明の実施形態について説明する。なお、以下の説明において参照する各図では、他の図と同等部分は同一符号によって示される。
【００１８】
（血流量予測方法）
以下、本実施形態に係る血流量予測方法について説明する。なお、以下の説明においては、本実施形態の血流量予測方法は、コンピュータ装置上で動作するシュミレーションプログラム（血流量予測プログラム）によって実現されるものとして説明する。
【００１９】
上述したように、脳血管系はその仕組が複雑であり、個人差も大きい。さらに、リモデリングに要する時間や血流再配分の達成度など多くの不確定要素がある。このため、数学的に厳密なモデルをたてて詳細な数値計算を行っても、これらの不確定要素が存在するために信頼性の高い評価は難しい。上述したように、このような場合は数学モデルを簡単にしつつ、できるだけ多くの実例に合うようにパラメータを調整する方法が適している。また、調整のためのパラメータは多いほど評価結果が精密になると期待される一方、その数が多すぎる場合には、症例に合わせて設定することの難易度が上がるために、かえって評価結果の精度が下がる可能性もある。よって、適度な数のパラメータを選択する必要がある。
【００２０】
本実施形態に係る血流量予測プログラムにおける計算モデルは、次の通りである。
【００２１】
まず、血管系統毎に、流量と流動抵抗の初期値が与えられる（これらの値は、例えばユーザ（医療従事者など）がコンピュータ装置のキーボード等の入力デバイスを操作して入力することによって与えられる）。また、初期状態は血管梗塞発生前の健常状態とする。リモデリングは、ある系統入口が狭隘化したような異常が発生したあと、全系統の合計流量と各系統の流量がともに、脳梗塞発生前の健全状態である初期値に近づくことを目的として機能するものとする。リモデリングが作用する状態に従って変化するパラメータは、各系統の流動抵抗と全系統の合計流量に設定される。
【００２２】
また、リモデリングは、流量の少ない系統は流動抵抗を小さくして他の系統から流量を引き込みやすくし、逆に流量の多い系統は流動抵抗を大きくして自身の流量を削減して他の系統に血液を供給するよう作用するものとする。流動抵抗と流量から各系統の圧力損失が決まり、系統間に生じた圧力差によって系統間のクロスフローが発生する。クロスフロー流量は、系統間の圧力差とクロスフローの流動抵抗とから求められる。
【００２３】
また、全系統の合計流量が初期値より小さいときは、合計流量が初期値に近づくよう各系統の入口流量が増加し、逆に初期値より大きいときは、合計流量が初各系統の入口流量が減少し初期値に漸近する働きがあるものとする。
【００２４】
上記のリモデリングは生体的な応答であるので、ある時定数の時間遅れをもって達成されるものとする。
【００２５】
ここで、血管系統数はN個あるものとし、n番目の系統の初期、すなわち健常時の流量をFLI(n)とすると、脳全体の必要流量FLNEEDは、次式のようにその総和である。
【００２６】
【数１】
2015018320.tif 000003
【００２７】
これらの血流はWillis動脈輪で合流しているため、血液の供給圧力はすべて同一と考えてよい。供給圧力を規格化して“1”とすると、各系統の初期の流動抵抗FRI(n)は流量の逆数で与えられる。
【００２８】
【数２】
2015018320.tif 000004
【００２９】
ここで、系統iと系統jの間のクロスフロー抵抗をFRC(i, j)とする。初期状態では系統間に圧力差がないのでクロスフロー流量FLC(n, i)はゼロである。また、初期状態での各系統への供給量FLIN(n)は健常時の系統流量に等しいものとする。これらは次式で表わされる。
【００３０】
【数３】
2015018320.tif 000005
【００３１】
【数４】
2015018320.tif 000006
【００３２】
脳梗塞が発生するイベントでは、発生する系統への供給流量が瞬時にゼロになり（すなわち、FLIN(n)＝0）、その状態が血流再建手術まで継続すると仮定する。また血流再建手術では、瞬時に初期流量すなわち健常時の系統流量が再び供給される（すなわち、FLIN(n)＝FLI(n)）と仮定する。本実施形態の血流量予測プログラムは、これらの外的条件により血流の非平衡が生じ、リモデリングにより各系統に必要な流量が再配分されるまでの過渡変化を計算することが目的である。
【００３３】
各系統における流量FL(n)は入口流量FLIN(n)とクロスフロー流量FLC(n, i)の和である。クロスフロー流量FLC(n, i)は、当該系統に流入するクロスフローについては正、流出するクロスフローについては負の値をとる。
【００３４】
【数５】
2015018320.tif 000007
【００３５】
各系統の圧力損失DP(n)はそれぞれの流量FL(n)と流動抵抗FR(n)の積で表わされる。
【００３６】
【数６】
2015018320.tif 000008
【００３７】
ここで、次のタイムステップのためのリモデリングを考える。各系統において、現在の流量FL(n)が健常時の流量FLI(n)より小さいときは流動抵抗FR(n)を減少させ、逆に大きいときは流動抵抗FR(n)を増して健常時の流量に近づくよう作用するものとする。生体応答の時間遅れをあらわす時定数τを考慮して、タイムステップΔt後の流動抵抗を次式で計算する。
【００３８】
【数７】
2015018320.tif 000009
【００３９】
ここで、タイムステップ前後の各パラメータは、それぞれ、上付き添え字のnewとoldで示す（以下、同様）。例えば、上記式では、タイムステップ前の流動抵抗はFRold(n)、タイムステップ後の流動抵抗はFRnew(n)である。タイムステップΔtが時定数τより十分小さいとき、指数関数をTaylor展開し高次の項を省略すると次のように近似できる。以下、同様の近似を用いる。
【００４０】
【数８】
2015018320.tif 000010
【００４１】
過渡計算の途中で各系統の流量を合わせた供給流量は次式の通りである。
【００４２】
【数９】
2015018320.tif 000011
【００４３】
全供給流量が脳の全必要量になるように、各系統の流量も同じ時間遅れをもって調整されるものとする。
【００４４】
【数１０】
2015018320.tif 000012
【００４５】
クロスフローは系統n, i間の圧力損失DP(n), DP(i)の差に比例し、クロスフロー流動抵抗FRC(n, i)に反比例するが、同じ時間遅れをもって追随するものと仮定する。
【００４６】
【数１１】
2015018320.tif 000013
【００４７】
このように、本実施形態の血流量予測プログラムでは、これらの過渡変化を所定の時間まで追跡することができる。
【００４８】
また、虚血による系統の血管ダメージ量を、健全時の血流量に対する減少比の時間積分から次式によって推定することができる。
【００４９】
【数１２】
2015018320.tif 000014
【００５０】
一方、過大な灌流による血管への負荷量を、健全時の血流量に対する増大比の時間積分から次式によって推定することができる。
【００５１】
【数１３】
2015018320.tif 000015
【００５２】
上記のように推定される血管のダメージ量Damage(n)と負荷量Load(n)との比較から、再出血リスクを評価することが可能である。
【００５３】
（処理フロー）
図１は、本実施形態に係る血流量予測プログラムの処理の一例を示すフロー図である。ユーザの入力を受け付けることにより（ステップＳ１０２）、各系統の初期流量FLI(n)または流動抵抗FRI(n)、クロスフロー流動抵抗FRC(n, i)およびリモデリングの時定数τを得てtime＝0の初期状態を設定する（ステップＳ１０４）。
【００５４】
次に、過渡計算に入ってタイムステップΔtを順次進め、当該タイムステップにおける血流閉塞や再建のイベントがあれば、イベントを設定する（ステップＳ１０６）。これらのイベントは、例えば、ユーザによって入力されることで設定されてもよい。
【００５５】
次に、各系統において前タイムステップの流量FL(n)から圧力損失DP(n)を算出し、入口流量FLIN(n)とクロスフロー流量FLC(n, i)と流動抵抗FR(n)とから当該タイムステップの各系統の流量FL(n)を算出する（ステップＳ１０８）。
【００５６】
次に、系統毎に初期流量FLI(n)と現流量（当該タイムステップでの流量）FL(n)との差に比例し、リモデリングの時間遅れをもって各系統の流動抵抗を変化させる（ステップＳ１１０）。全流量が必要流量（初期の全流量）に満たないとき、各系統の入口流量FLIN(n)をステップＳ１１０と同じリモデリングの時間遅れをもって増加させる。また、全流量が必要流量より多い場合には、各系統の入口流量FLIN(n)を減少させる（ステップＳ１１２）。
【００５７】
時間はタイムステップを進め（time＝time＋Δt）（ステップＳ１１４）、同じ過渡計算を所定の終了時間まで繰り返す（ステップＳ１１６）。
【００５８】
ステップＳ１０２において、入力で設定されるべきパラメータは、各系統の流量と流動抵抗値の各初期値FLI(n)，FRI(n)、クロスフローの接続関係（いずれの系統間でクロスフローが発生するか）とクロスフロー流動抵抗値FR(n, i)、リモデリングに要する過渡変化の時定数τである。そして、上記処理フローの出力として得られるのは、各系統の流量と流動抵抗値の過渡変化（すなわち、各タイムステップにおける各系統の流量と流動抵抗値）である。これらの計算結果から、上記式（１２）や式（１３）により、過灌流や脳組織ダメージを評価することが可能である。
【００５９】
（具体例）
簡単な系統を扱った本実施形態に係る血流量予測プログラムによる計算結果の具体例を図２〜図６に示す。本具体例では、前大脳動脈、中大脳動脈、後大脳動脈の各系統を単純化し、Ａ系統、Ｂ系統、Ｃ系統の流体力学的に同等な３つの系統に血流が供給されるものとする。ただし、Ａ−Ｂ系統の間、Ｂ−Ｃ系統の間にはクロスフローが生じうるが、Ａ−Ｃ系統の間にはクロスフローは生じないものとする。各系統の流量配分は入口流量と流動抵抗により決まるものとする。健全な初期状態の流量および流動抵抗値は各系統で無次元数の“１”とする。これは初期値で規格化した相対値とみなすことができる。ここでは、クロスフローの抵抗値も簡単のため“１”とする。リモデリングに要する時間遅れの時定数はタイムステップでもある任意時間単位の１００倍である。（言い換えると、リモデリングの時定数の１００分の１を１時間単位と設定したということである。）
以下に説明する４つの計算例により、本実施形態に係る血流量予測プログラムによって、血流再建手術のタイミングによる過灌流と虚血領域の脳組織ダメージについて定性的に目安としうる結果が得られることが分かる。
【００６０】
［計算例１］
計算例１では、初期は健全状態であるとし、５０時間単位の時点でＡ系統において脳梗塞または動脈痙攣で狭窄が発生したケースを考える。狭窄は時間がかかる場合もあるが、ここでは簡単のために瞬時にＡ系統の流量がゼロになると仮定する。また、血流再建手術まではＡ系統はクロスフローによってのみ血流が供給される。５００時間単位後にステントまたはバイパスによる血流再建手術が行われ、その後、再度のリモデリングにより初期の健全状態に漸近する。以上の経過のシミュレーション結果を図２に示す。
【００６１】
発症後、流量が減少したＡ系統において流動抵抗が減少するリモデリングにより、Ｂ系統からのクロスフローによって徐々に血流量が回復する。Ｂ系統はクロスフローによる流出のため流量がやや低下し、Ｃ系統は全流量の低下を補うために流量が増加する。５００時間単位の後にこれらの作用の結果として全系統の流量がほぼ初期流量と等しくなる。この間のＡ系統の虚血状態の継続時間は、血流が遮断されているために酸素や栄養の供給が妨げられることでこの範囲の脳組織に与えるダメージの目安となる。
（５００時間単位において血流再建手術が行われ、その結果Ａ系統入口に健全時と同じ流量の血流が回復するものとする。Ａ系統は入口流量とＢ系統からのクロスフローで過灌流の状態になり、ピーク値は初期流量の約２倍である。全流量も過大になるため、Ｂ系統とＣ系統の流量も低下する。健全状態に近づけるリモデリング効果により、約１０００時間単位を経て初期状態と同じ流量配分および流動抵抗分布が再現される。
【００６２】
［計算例２］
計算例２では、上記計算例１よりも血流再建手術が早期に行われ、１００時間単位において実施されたケースを考える。本例のシミュレーション結果を図３に示す。Ａ系統の血流は一時的にゼロになるがリモデリングによる流動抵抗の変化が進まないうちに血流再建されるため、過灌流も初期状態の１．３倍程度に収まっている。虚血状態の持続時間も短いので、脳組織のダメージも少ないと考えられる。シミュレーション結果のグラフ化により発症後の早期の処置が重要であることを視覚的に認識することができる。
【００６３】
［計算例３］
計算例３では、上記具体例２よりやや遅く発症後の２００時間単位において血流再建手術が行われたケースを考える。本例のシミュレーション結果を図４に示す。手術前にリモデリングが進み、過灌流のピークは初期流量の１．６倍になる。この結果により、上記計算例２と比べて発症後の処置を講じる時間差の影響が表れていることが分かる。
【００６４】
［計算例４］
計算例４では、発症後に血流再建手術が行われなかったケースを考える。本例のシミュレーション結果を図５に示す。時間をかけて各系統の流動抵抗が調整され、Ａ系統の入口流量がゼロのままであるにもかかわらず、約１０００時間単位の後では各系統の流量が初期流量と同じ状態に到達する。しかしながら、この状態が望ましいものでないことは、図６に示される各系統の流動抵抗値の変化から分かる。
【００６５】
図６に示されるように、最終的な流動抵抗値については、Ｂ系統はほぼ初期状態のままであるが、Ａ系統は約０．２倍となる一方、Ｃ系統は約１．６倍になる必要がある。脳血管のリモデリングにおいてこれほどの流動抵抗値の変化が起これば、脳血管系の構造に多大な負担となることが推定できる。
【００６６】
（コンピュータ装置の構成例）
図７は、本実施形態に係る血流量予測プログラムを実行するコンピュータ装置のハードウェア構成の一例を示す図である。
【００６７】
図７に示されるように、本実施形態に係る血流量予測プログラムを実行するコンピュータ装置１は、ＣＰＵ（Central Processing Unit）１１、ＲＯＭ１３、ＲＡＭ１５、ハードディスク１７等の記憶デバイス、ディスプレイが接続されるためのディスプレイインタフェース１９、キーボードやマウス等のユーザの入力操作のためのユーザ入力デバイス２１、リムーバブルメディアが接続されるためのリムーバブルメディアインタフェース２３、プリンタが接続されるためのプリンタインタフェース２５、ネットワークと接続されるためのネットワークインタフェース２７等の一般的なコンピュータの構成と同様の構成により実現されうる。コンピュータ装置１は、上記説明した各処理を実行するためのプログラムを、例えば、ハードディスク１７等の記憶装置に記憶している。そして、ＣＰＵ１１がハードディスク１７等に記憶されているプログラムを読み出して実行することによって、上記説明した各処理が実現されうる。また、各処理に使用される各種のデータは、各装置のハードディスク１７やＲＡＭ１５等に記憶領域に随時記憶されうる。
【００６８】
なお、上述したコンピュータ装置のハードウェア構成はあくまで一例であり、これに限定されるものではない。本実施形態に係る血流量予測プログラムが実行可能な構成であれば、他のハードウェア構成を採用可能である。
【００６９】
（まとめ）
脳梗塞等の脳動脈血管閉塞の発症後、また手術後の過渡血流挙動の概略をあらわすことができる。また、症例適用の蓄積により、流動抵抗やリモデリング時間遅れなどのパラメータの精度を上げることができる。脳梗塞患者の閉塞率と経過時間から、血流再建手術後の過灌流を予測し、過大流量のピークによる荷重に血流欠損時に損傷を受けた血管の耐久性を考慮した再建手術可否の判断ができる。また、再出血リスクを回避する手順（段階的血流再建など）を検討する手段を提供することができる。
【００７０】
ここまで、本発明の一実施形態について説明したが、本発明は上述した実施形態に限定されず、その技術的思想の範囲内において種々異なる形態にて実施されてよいことは言うまでもない。
【００７１】
また、本発明の範囲は、図示され記載された例示的な実施形態に限定されるものではなく、本発明が目的とするものと均等な効果をもたらすすべての実施形態をも含む。さらに、本発明の範囲は、各請求項により画される発明の特徴の組み合わせに限定されるものではなく、すべての開示されたそれぞれの特徴のうち特定の特徴のあらゆる所望する組み合わせによって画されうる。
【符号の説明】
【００７２】
１  コンピュータ装置
１１  ＣＰＵ
１３  ＲＯＭ
１５  ＲＡＭ
１７  ハードディスク
１９  ディスプレイインタフェース
２１  ユーザ入力デバイス
２３  リムーバブルメディアインタフェース
２５  プリンタインタフェース
２７  ネットワークインタフェース
【図１】
2015018320.tif 000016
【図２】
2015018320.tif 000017
【図３】
2015018320.tif 000018
【図４】
2015018320.tif 000019
【図５】
2015018320.tif 000020
【図６】
2015018320.tif 000021
【図７】
2015018320.tif 000022
【図８】
2015018320.tif 000023