"穿刺針が刺入される穿刺対象臓器の血流量を推定するシステムであって、前記穿刺針の刺入部位に一定の熱量を付与する加熱手段と、前記熱量の付与後における前記加熱手段の近傍の温度変化を測定可能な温度センサと、前記温度センサの計測値に基づき前記血流量を推定する血流量推定手段とを備え、前記血流量推定手段では、予め記憶された数式により、前記温度センサの計測値から時間に対する温度上昇率を算出し、当該温度上昇率に基づいて前記血流量を算出することを特徴とする穿刺対象臓器の血流量推定システム。前記加熱手段と前記温度センサは、前記穿刺針に一体的に設けられていることを特徴とする請求項1記載の血流量推定システム。前記穿刺針は、電磁波による熱エネルギーを前記刺入部位に付与するための電極を有する電極針であり、前記加熱手段による加熱は、前記電極による加熱と別個独立に行われることを特徴とする請求項1又は2記載の血流量推定システム。請求項3記載の血流量推定システムを備えた温度分布推定システムであって、前記血流量に基づき、前記穿刺対象臓器内での温度分布を推定する温度分布推定手段を更に備え、前記温度分布推定手段では、予め記憶された熱伝導方程式に、前記血流量と、前記電極から前記穿刺対象臓器に与えられた熱量とを代入することで、前記温度分布を算出することを特徴とする穿刺対象臓器の温度分布推定システム。穿刺針が刺入された穿刺対象臓器の状態を解析する解析装置において、前記穿刺針の刺入部位に一定の熱量を付与したときのその近傍の温度変化に基づき、前記穿刺対象臓器の血流量を推定する血流量推定手段を備え、前記血流量推定手段は、予め記憶された数式により、時間に対する前記近傍の温度上昇率を算出し、当該温度上昇率に基づいて前記血流量を算出することを特徴とする穿刺対象臓器の解析装置。穿刺針が刺入された穿刺対象臓器の状態を解析する解析装置のコンピュータを機能させるプログラムにおいて、前記穿刺針の刺入部位に一定の熱量を付与したときのその近傍の温度の計測値から、予め記憶された数式により、時間に対する前記近傍の温度上昇率を算出し、当該温度上昇率に基づいて前記血流量を算出する血流量推定手段として前記コンピュータを機能させることを特徴とする穿刺対象臓器の解析装置用プログラム。しかしながら、前記従来の電極針にあっては、臓器の刺入部位の温度管理しかできないため、医師は、術中に、焼灼対象となる腫瘍或いは臓器全体の温度分布及び当該温度分布の経時的な変化を把握することができず、次の不都合を招来する。一般的に、RFAは、超音波画像下で医師の目視によって行われるが、超音波画像はモノクロであるため、腫瘍と焼灼部分の境界が不明瞭で目視しにくい。また、電極針から臓器に所定の熱量を付与した際に、当該臓器の血流量が多いと血流に熱が奪われ易くなり、腫瘍への伝熱率が低下して焼灼範囲が意図した範囲よりも小さくなる。逆に、臓器の血流量が少ないと腫瘍への伝熱率が高くなり、焼灼範囲が意図した範囲よりも大きくなる。血流量は、術中の患者の状態によって変化することから、刺入部位の温度をモニタリングし、電極針から臓器に与えられるジュール熱の熱量を所望の値に維持していても、血流量の変化によって臓器内の焼灼範囲の大きさが変わることになる。このため、電極針からのジュール熱の出力調整は、医師の経験や勘に依存しており、治療成績にばらつきが生じ易い。以上により、臓器の焼灼範囲を正確にコントロールするには、術中、血流量によって変わる臓器の温度分布及びその変化に関する情報を逐次取得することが重要となる。従って、電極針の刺入部位の温度管理をしていても、焼灼対象となる腫瘍或いは臓器全体の温度分布とその変化を術中に正確に把握できなければ、正確な範囲で焼灼することができず、腫瘍の未焼灼領域の残留や正常組織への焼灼過多を招来する虞がある。本発明は、このような不都合に着目して案出されたものであり、その目的は、穿刺対象臓器内の血流量を推定し、当該血流量から穿刺対象臓器の温度分布とその変化を把握することに寄与する穿刺対象臓器の血流量推定システム、温度分布推定システム、解析装置、及び解析装置用プログラムを提供することにある。前記目的を達成するため、本発明は、穿刺対象臓器における穿刺針の刺入部位に一定の熱量を付与したときのその近傍の温度の計測値から、予め記憶された数式により、時間に対する前記近傍の温度上昇率を算出し、当該温度上昇率に基づいて前記血流量を算出する血流量推定手段を備えた点に主たる特徴を有する。本発明によれば、穿刺術中において、穿刺対象臓器内における血流量を経時的に推定することができ、当該血流量に左右する穿刺対象臓器の温度分布の推定を経時的に行うことができる。これにより、この温度分布に関する情報を医師に提示可能になる他、当該情報に基づいて、穿刺対象臓器を加熱する電極からの出力エネルギーを自動的に制御することができ、患者の術中の血流量の変化に拘らず、所望の範囲の焼灼をより正確に行うことができる。以下、本発明の実施形態について図面を参照しながら説明する。図1には、本実施形態に係る温度分布推定システムの概略構成図が示されている。この図において、前記温度分布推定システム10は、穿刺対象臓器の腫瘍部分に穿刺針としての電極針12を刺し、当該電極針12に供給されるラジオ波によるジュール熱で腫瘍部分を焼灼するRFAラジオ波焼灼療法時に、穿刺対象臓器の温度分布を推定するシステムである。この温度分布推定システム10は、図示しないラジオ波発生装置からの高周波電流が通電する電極12Aを有する電極針12にそれぞれ内蔵された加熱手段14及び温度センサ15と、温度センサ15の計測値に基づき、穿刺対象臓器内の温度分布を経時的に求める解析装置17とを備えている。前記加熱手段14は、図示しない電源に接続され、電極12Aとは別の加熱用細線からなり、電極12Aからのラジオ波エネルギーによる刺入部位の加熱とは別個独立して、当該刺入部位の周囲に一定の熱量を継続して供給するようになっている。なお、加熱手段14としては、加熱用細線に限定されず、腫瘍部分の凝固壊死させるための電磁波による熱エネルギーの付与とは別に、刺入部位の周囲に一定の熱量を継続して供給できる限りにおいて種々のものを採用することができる。前記温度センサ15は、加熱手段14の近傍の温度、すなわち、刺入部位の近傍の電極針12の温度を計測可能になっており、当該計測値に対応する電気信号が経時的に逐次、解析装置17に送信される。つまり、術中、加熱手段14から穿刺対象臓器に一定熱量が与えられ、それによって生じる温度上昇が温度センサ15で経時的に計測される。前記解析装置17は、CPU等の演算処理装置及びメモリやハードディスク等の記憶装置等からなるコンピュータによって構成され、当該コンピュータを以下の各手段として機能させるためのプログラムがインストールされている。この解析装置17は、温度センサ15の計測値に基づき、穿刺対象臓器の血流量の推定値を経時的に求める血流量推定手段19と、血流量推定手段19で求めた血流量に基づき、穿刺対象臓器内での各部位の温度の推定値を経時的に求める温度分布推定手段20とを備えている。前記血流量推定手段19では、予め記憶された以下の式1により、温度センサ15の計測値温度Tから、加熱手段14による継続した加熱の開始時からの経過時間tにおける電極針12の温度上昇率Uが算出され、当該温度上昇率Uから血流量Fが算出される。なお、温度上昇率Uと血流量Fの関係は、本発明者らが鋭意、実験研究を行った結果、式1として導出されたものであり、同式1中、a、bは、数値解析に基づくシミュレーションによる実験等の実験研究によって特定された定数となる。また、式1は、加熱手段14からの発熱量のうち、血流量が多い場合、血流に吸収される熱量が多くなって電極針12の温度上昇量が小さくなる一方、血流量が少ない場合、血流に吸収される熱量が少なくなって電極針12の温度上昇量が大きくなる点に着目し、実験研究により導出されたものである。前記温度分布推定手段20では、前記経過時間tにおける血流量Fと、電極12Aから穿刺対象臓器に与えられるジュール熱の熱量QPとが、予め記憶された次式2の熱伝導方程式に代入され、当該方程式を解くことで、穿刺対象臓器の温度分布が算出される。ここで、温度分布は、穿刺対象臓器内の適当な基準位置に設定された座標系から見た任意要素に対する位置ベクトルBと、経過時間tとで表される関数θB,tの形で求められる。また、ジュール熱の熱量QPは、電極12Aへの供給電力量Pと肝臓の電気伝導率によって求まる関数であり、経過時間tにおける前記供給電力量Pの入力によって特定される。上式2において、ρは穿刺対象臓器の密度、cは穿刺対象臓器の比熱、λは穿刺対象臓器の熱伝導率、ρbは血液の密度、cbは血液の比熱、θbは血液の温度であり、これら数値ρ,c,λ,ρb,cb,θbは、定数として予め記憶されている。すなわち、ここでは、経過時間t毎の温度関数θB,tが求められ、穿刺対象臓器の座標系における各位置の温度が経過時間t毎に特定され、同一となる温度部分毎に肝臓を領域分けした温度分布が経時的に得られる。ここで得られた温度分布のデータは、図示しない表示装置に表示することでRFAを行う医師に提示し、或いは、電極12Aからのジュール熱の発熱量の自動制御に用いることができる。なお、前記実施形態では、加熱手段14と温度センサ15を電極針12に内蔵した場合を説明したが、本発明はこれに限らず、これらを電極針12と別に設けることも可能である。また、前記温度分布推定システム10は、穿刺対象臓器の血流量を推定する血流量推定システムを含んで構成されているが、本発明では、前記温度分布推定手段20を設けずに、血流量推定システムのみの構成とすることもできる。更に、本発明は、RFAの際の適用のみならず、マイクロ波凝固療法PMCT等の他の穿刺術を行う際の適用も可能となる。その他、本発明における装置各部の構成は図示構成例に限定されるものではなく、実質的に同様の作用を奏する限りにおいて、種々の変更が可能である。"