(57)【要約】【課題】測定された血管内圧の時間変化の波形から一心拍毎に前負荷の指標となる値を算出する循環系指標算出プログラム、循環系指標算出装置、循環系指標算出システム及び循環系指標算出方法を提供する。【解決手段】循環系指標算出装置１は、患者４の血管内圧の時間変化のうち、心周期の拡張期の波形に対する回帰曲線を一心拍毎に求め、当該回帰曲線の収束値から患者４の循環系の指標を算出する循環系指標算出手段１０２と、算出された循環系の指標の時間変化を表示処理する表示処理手段１０３とを有する。【選択図】図２
2019062990.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータを、
患者の血管内圧の時間変化のうち、心周期の拡張期の波形に対する回帰曲線を一心拍毎に求め、当該回帰曲線の収束値から当該心拍における前記患者の循環系の指標を算出する算出手段として機能させるための循環系指標算出プログラム。
【請求項２】
前記算出手段は、前記血管内圧として動脈圧の時間変化に対する回帰曲線を求め、循環系としての血管の流路抵抗と血管のコンプライアンスに対する一心拍の圧力応答を表す減衰関数を前記回帰曲線とする請求項１に記載の循環系指標算出プログラム。
【請求項３】
前記算出手段は、複数心拍の前記収束値の平均から前記循環系の指標を算出する請求項１又は２に記載の循環系指標算出プログラム。
【請求項４】
前記算出手段は、前記循環系の指標として平均循環充満圧を算出する請求項１〜３のいずれか１項に記載の循環系指標算出プログラム。
【請求項５】
前記算出手段は、前記減衰関数の前記収束値を算出するとともに、さらに前記減衰関数の時定数を算出する請求項２に記載の循環系指標算出プログラム。
【請求項６】
前記コンピュータを、
前記算出手段が算出した前記循環系の指標の時間変化を表示処理する表示処理手段としてさらに機能させる請求項１〜５のいずれか１項に記載の循環系指標算出プログラム。
【請求項７】
前記表示処理手段は、前記循環系の指標の時間変化と、前記患者の生体情報の時間変化とを同期して表示する請求項６に記載の循環系指標算出プログラム。
【請求項８】
患者の血管内圧の時間変化のうち、心周期の拡張期の波形に対する回帰曲線を一心拍毎に求め、当該回帰曲線の収束値から当該心拍における前記患者の循環系の指標を算出する算出手段を有する循環系指標算出装置。
【請求項９】
患者の血管内圧を測定する測定器と、
前記測定器によって測定された前記患者の血管内圧の時間変化のうち、心周期の拡張期の波形に対する回帰曲線を一心拍毎に求め、当該回帰曲線の収束値から当該心拍における前記患者の循環系の指標を算出する算出手段を有する循環系指標算出装置と、
前記循環系指標算出装置が算出した前記循環系の指標の時間変化を表示処理する表示装置とを備える循環系指標算出システム。
【請求項１０】
患者の血管内圧の時間変化のうち、心周期の拡張期の範囲を定めるステップと、
心周期の拡張期の範囲の波形に対する回帰曲線を一心拍毎に求めるステップと、
当該回帰曲線の収束値から当該心拍における前記患者の循環系の指標を算出するステップとを有する循環系指標算出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、循環系指標算出プログラム、循環系指標算出装置、循環系指標算出システム及び循環系指標算出方法に関する。
【背景技術】
【０００２】
医学の分野において、心臓、動脈、静脈及び末梢血管を含む循環系の状態を把握するために静脈圧、動脈圧、心拍出量、血液量等の循環系の指標となる値を知ることは、術中の患者の状態を管理するために非常に重要である。
【０００３】
従来の技術として、循環系の状態を把握するための指標の一つとして心拍出量を測定する循環系指標算出装置が提案されている（例えば、特許文献１参照）。
【０００４】
特許文献１に開示された循環系指標算出装置は、熱希釈法を用いて心拍出量を求める際に用いられるものであって、患者の血液温度情報と心電図とを取得し、患者の心臓内に冷水を注入したときの血液温度情報に基づいて熱希釈曲線と心拍出量を算出し、患者の血液温度情報と同時に測定された心電図と熱希釈曲線とを同一画面内に同期させて表示するとともに、心電図のうち期外収縮を含む第一心電図と期外収縮を含まない第二心電図とを異なる表示態様で画面に表示する。医師は、表示された画面を確認し、測定中における期外収縮の有無を判断して、測定結果における期外収縮の影響を考察する。
【０００５】
上記した特許文献１の循環系指標算出装置は、測定中における期外収縮の有無を判断するための情報を表示するため、当該期外収縮の有無を判断するための情報がない場合に比べて再計測の回数を抑制し、患者の負担を抑制するものの、熱希釈法は患者の心臓内に冷水を注入する必要があり、比較的大きな侵襲を伴うものである。
【０００６】
そこで、熱希釈法に比べて侵襲の程度を抑制したものとして、測定された動脈圧波形から心拍出量を算出する循環系指標算出装置が提案されている（例えば、非特許文献１参照）。
【０００７】
非特許文献１に開示された循環系指標算出装置は、血管コンプライアンスを考慮して測定された動脈圧波形を動脈血管量波形へ変換し、変換した動脈血液量波形についてデータを評価するために自己相関処理を行って補正前の心拍出量を算出し、ノモグラムを用いて校正係数を決定して補正前の心拍出量から補正後の心拍出量を算出して表示部に表示する。動脈圧波形から動脈血管量波形への変換は、患者の身長、体重、年齢等の固有情報により定まる係数を用いて予め定められた数式により行われる。また、自己相関処理では、心周期ごとの波形を求めることで心拍数が算出されるとともに、補正前の一回拍出量が自己相関結果から算出され、心拍数と一回拍出量との積から補正前の心拍出量が求められる。なお、ノモグラムの算出方法は非公開であるが、患者の身長、体重、年齢等の固有情報に関するファクターを変数として用いて算出されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２０１６‐１８７５３８号公報
【非特許文献】
【０００９】
【非特許文献１】佐藤慎他、“機器紹介  低侵襲連続的心拍出量モニタシステム「LiDCOrapid」”、循環制御、日本、２０１４年、第３５巻、第２号、１６６‐１７６頁
【発明の概要】
【発明が解決しようとする課題】
【００１０】
上記した非特許文献１の循環系指標算出装置は、測定された動脈圧波形から心拍出量を算出するものであるが、心拍出量は前負荷及び後負荷の影響を受けるものであり、いずれか一方のみに起因する指標を測定するものではない。
【００１１】
従って、本発明の目的は、測定された血管内圧の時間変化の波形から一心拍毎に前負荷の指標となる値を算出する循環系指標算出プログラム、循環系指標算出装置、循環系指標算出システム及び循環系指標算出方法を提供することにある。
【課題を解決するための手段】
【００１２】
本発明の一態様は、上記目的を達成するため、以下の循環系指標算出プログラム、循環系指標算出装置、循環系指標算出システム及び循環系指標算出方法を提供する。
【００１３】
［１］コンピュータを、
患者の血管内圧の時間変化のうち、心周期の拡張期の波形に対する回帰曲線を一心拍毎に求め、当該回帰曲線の収束値から当該心拍における前記患者の循環系の指標を算出する算出手段として機能させるための循環系指標算出プログラム。
［２］前記算出手段は、前記血管内圧として動脈圧の時間変化に対する回帰曲線を求め、循環系としての血管の流路抵抗と血管のコンプライアンスに対する一心拍の圧力応答を表す減衰関数を前記回帰曲線とする前記［１］に記載の循環系指標算出プログラム。
［３］前記算出手段は、複数心拍の前記収束値の平均から前記循環系の指標を算出する前記［１］又は［２］に記載の循環系指標算出プログラム。
［４］前記算出手段は、前記循環系の指標として平均循環充満圧を算出する前記［１］〜［３］のいずれかに記載の循環系指標算出プログラム。
［５］前記算出手段は、前記減衰関数の前記収束値を算出するとともに、さらに前記減衰関数の時定数を算出する前記［２］に記載の循環系指標算出プログラム
［６］前記コンピュータを、
前記算出手段が算出した前記循環系の指標の時間変化を表示処理する表示処理手段としてさらに機能させる前記［１］〜［５］のいずれかに記載の循環系指標算出プログラム。
［７］前記表示処理手段は、前記循環系の指標の時間変化と、前記患者の生体情報の時間変化とを同期して表示する前記［６］に記載の循環系指標算出プログラム。
［８］患者の血管内圧の時間変化のうち、心周期の拡張期の波形に対する回帰曲線を一心拍毎に求め、当該回帰曲線の収束値から当該心拍における前記患者の循環系の指標を算出する算出手段を有する循環系指標算出装置。
［９］患者の血管内圧を測定する測定器と、
前記測定器によって測定された前記患者の血管内圧の時間変化のうち、心周期の拡張期の波形に対する回帰曲線を一心拍毎に求め、当該回帰曲線の収束値から当該心拍における前記患者の循環系の指標を算出する算出手段を有する循環系指標算出装置と、
前記循環系指標算出装置が算出した前記循環系の指標の時間変化を表示処理する表示装置とを備える循環系指標算出システム。
［１０］患者の血管内圧の時間変化のうち、心周期の拡張期の範囲を定めるステップと、
心周期の拡張期の範囲の波形に対する回帰曲線を一心拍毎に求めるステップと、
当該回帰曲線の収束値から当該心拍における前記患者の循環系の指標を算出するステップとを有する循環系指標算出方法。
【発明の効果】
【００１４】
請求項１、８〜１０に係る発明によれば、測定された血管内圧の時間変化の波形から一心拍毎に前負荷の指標となる値を算出することができる。
請求項２に係る発明によれば、血管内圧のうち動脈圧の時間変化に対する回帰曲線を求め、循環系としての血管の流路抵抗と血管のコンプライアンスに対する一心拍の圧力応答を表す減衰関数を回帰曲線として用いることができる。
請求項３に係る発明によれば、複数心拍の収束値の平均から循環系の指標を算出することができる。
請求項４に係る発明によれば、循環系の指標として平均循環充満圧を算出することができる。
請求項５に係る発明によれば、減衰関数の収束値を算出するとともに、さらに減衰関数の時定数を算出することができる。
請求項６に係る発明によれば、算出手段が算出した循環系の指標の時間変化を表示処理することができる。
請求項７に係る発明によれば、循環系の指標の時間変化と、患者の生体情報の時間変化とを同期して表示することができる。
【図面の簡単な説明】
【００１５】
【図１】図１は、実施の形態に係る循環系指標算出システムの構成の一例を示す概略図である。
【図２】図２は、実施の形態に係る循環系指標算出装置の構成例を示すブロック図である。
【図３】図３は、動脈圧情報の構成の一例を示す概略図である。
【図４】図４は、生体情報の構成の一例を示す概略図である。
【図５】図５は、循環系指標情報の構成の一例を示す概略図である。
【図６】図６は、循環系指標算出手段の動作を説明するための動脈圧波形及びその微分値を示すグラフ図である。
【図７】図７は、循環系指標算出手段の動作を説明するための一心拍の動脈圧波形を示すグラフ図である。
【図８】図８は、循環系指標算出手段のフィッティング動作を説明するためのグラフ図である。
【図９】図９は、表示処理手段の表示動作の一例を説明するための概略図である。
【図１０】図１０は、循環系指標算出システムの動作を説明するためのフローチャートである。
【図１１】図１１は、心拍出量及び静脈還流量と平均右心房圧との関係を示すグラフ図である。
【図１２】図１２は、動脈圧の減衰を表すモデルとしての回路図である。
【発明を実施するための形態】
【００１６】
［実施の形態］
（循環系指標算出システムの構成）
図１は、実施の形態に係る循環系指標算出システムの構成の一例を示す概略図である。
【００１７】
この循環系指標算出システム６は、患者４の手術を行う場合等に、測定された患者４の血管内圧から循環系の指標となる値を算出し、医師５に患者４の状態を知らせるために、循環系の指標となる値を表示するためのもの、又は循環系の指標となる値を生体情報と同期して表示するためのものである。ここで、循環系指標算出システム６は、血管内圧の一例として動脈圧を測定し、循環系の指標となる値として平均循環充満圧を算出するものであり、平均循環充満圧は一心拍毎に算出するものとする。
【００１８】
循環系指標算出システム６は、表示部１２と操作部１３とを備えた専用に設計された機器又はＰＣ（Ｐｅｒｓｏｎａｌ  Ｃｏｍｐｕｔｅｒ）やタブレット端末等の情報処理装置であって情報を処理する循環系指標算出装置１と、例えばＡラインにより患者４の動脈圧を測定する動脈圧センサ２と、患者４の心電図や血圧、脈拍、脳波、呼吸等の生体情報を測定する生体センサ３と、動脈圧センサ２及び生体センサ３で測定された情報並びに循環系指標算出装置１で算出された情報を表示する生体情報モニタ３Ａとを有する。循環系指標算出装置１、動脈圧センサ２、生体センサ３、生体情報モニタ３Ａは、医師５によって操作される。
【００１９】
循環系指標算出装置１は、動脈圧センサ２及び生体センサ３から得られる情報に基づいて、患者４の状態を示す指標として平均循環充満圧を算出して生体情報とともに経時的に表示部１２及び／又は生体情報モニタ３Ａに表示するものであって、本体内に情報を処理するための機能を有するＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）やＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）、フラッシュメモリ等の電子部品を備える。
【００２０】
なお、循環系指標算出装置１は、表示部１２及び操作部１３を有しないサーバ装置として構成してもよく、その場合は端末装置の要求に応じて動作し、サーバ装置として構成する場合は遠隔地に配置されるものであってもよい。また、循環系指標算出装置１は、患者４の手術を行う手術室に配置されるものであってもよいし、複数の手術室を監視する部屋に配置されて複数の生体情報モニタ３Ａに接続されるものであってもよい。また、循環系指標算出装置１は、算出した指標に基づいて、例えば、患者４に輸液、輸血するように図示しない注入装置を制御するものであってもよい。
【００２１】
動脈圧センサ２は、患者４の動脈に挿入される針と、針と一端が接続された動脈ラインと、動脈ラインの他端に設けられて動脈圧に比例した電圧を出力する圧トランスデューサとを有する。圧トランスデューサの出力信号は、生体情報モニタ３Ａを介して循環系指標算出装置１に入力される。
【００２２】
各機器は、専用線により相互に通信可能に接続されるが、有線又は無線の通信ネットワークにより接続されてもよいし、イントラネットやＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）等の通信網を用いて接続されるものであってもよい。
【００２３】
循環系指標算出装置１は、上記構成において、生体情報モニタ３Ａを介して動脈圧センサ２及び生体センサ３から出力される情報を受信し、受信した動脈圧の情報に基づいて平均循環充満圧を算出するとともに、生体センサ３から得られた生体の情報を適宜組み合わせて表示部１２に表示する。医師５は、表示部１２又は生体情報モニタ３Ａに表示された情報を確認し、患者４の状態を判断して輸液、輸血等の処置を行う。
【００２４】
なお、循環系指標算出装置１、動脈圧センサ２、生体センサ３及び生体情報モニタ３Ａの機能の全部又は一部を一体に構成してもよいし、各装置の機能の一部又は全部を他の装置に含めてもよい。特に、循環系指標算出装置１の機能を生体情報モニタ３Ａに内蔵するように構成してもよい。また、循環系指標算出装置１、動脈圧センサ２及び生体センサ３の機能の全部又は一部を遠隔地に配置された機器で動作させるように構成してもよい。また、循環系指標算出装置１に対して複数の生体情報モニタ３Ａを対応させ、循環系指標算出装置１に同時に複数の患者の循環系の指標を算出させるものであってもよい。また、循環系指標算出装置１は、動脈圧センサ２及び生体センサ３から直接、動脈圧情報及び生体情報を受信してもよいし、図１に図示しているように生体情報モニタ３Ａを介して受信してもよい。
【００２５】
（循環系指標算出装置の構成）
図２は、実施の形態に係る循環系指標算出装置１の構成例を示すブロック図である。
【００２６】
循環系指標算出装置１は、ＣＰＵ等から構成され、各部を制御するとともに、各種のプログラムを実行する制御部１０と、ＨＤＤやフラッシュメモリ等の記憶媒体から構成され情報を記憶する記憶部１１と、画像及び文字により情報を表示する表示部１２と、操作内容に応じて制御部１０に対して操作信号を出力する操作部１３と、外部装置と通信する通信部１４とを備える。
【００２７】
制御部１０は、後述する循環系指標算出プログラム１１０を実行することで、動脈圧情報取得手段１００、生体情報取得手段１０１、循環系指標算出手段１０２、表示処理手段１０３及び情報出力手段１０４等として機能する。
【００２８】
動脈圧情報取得手段１００は、通信部１４及び生体情報モニタ３Ａを介して動脈圧センサ２から動脈圧に比例した信号を取得して動脈圧の数値に変換し、時刻情報とともに動脈圧情報１１１として記憶部１１に記録する。
【００２９】
生体情報取得手段１０１は、通信部１４及び生体情報モニタ３Ａを介して生体センサ３から心電図や血圧計、パルスオキシメータ等の情報を取得し、時刻情報とともに生体情報１１２として記憶部１１に記録する。
【００３０】
循環系指標算出手段１０２は、動脈圧情報１１１から循環系の指標となる値として平均循環充満圧を算出し、時刻情報とともに循環系指標情報１１３として記憶部１１に記録する。
【００３１】
循環系指標算出手段１０２は、動脈圧情報１１１を、時刻を横軸に動脈圧を縦軸に描いた波形に対し、心周期の拡張期の波形に対する回帰曲線を一心拍毎に求め、当該回帰曲線の収束値から平均循環充満圧を求める。回帰曲線は、毛細血管を含む動脈と静脈間の全血管の流路抵抗と、血管自体の拡張収縮に伴う管内容量の変化とを考慮したモデルとし、指数関数的に減衰するものを採用し、当該回帰曲線を定義するパラメータから収束値と時定数が算出される。詳しくは後述する。
【００３２】
表示処理手段１０３は、動脈圧情報１１１、生体情報１１２及び循環系指標情報１１３の全て又は適宜選択した一部を表示部１２に表示処理する。表示方法は、リアルタイム表示、履歴表示、複数拍動の平均値表示、予測値の表示を適宜選択することができ、複数の情報を表示する場合には各情報を時刻で同期して表示するものとする。表示方法は数値によるもの、グラフによるもの、色によるもの等、その方法は限定されない。
【００３３】
情報出力手段１０４は、動脈圧情報１１１及び循環系指標情報１１３の全て又は適宜選択した一部を、生体情報モニタ３Ａに表示された生体情報１１２とともに表示するように生体情報モニタ３Ａに対して出力するものである。特に、情報出力手段１０４は、生体情報モニタ３Ａにおいて動脈圧情報１１１、生体情報１１２及び循環系指標情報１１３が同期して表示されるように情報を出力することが好ましい。
【００３４】
記憶部１１は、制御部１０を上述した各手段１００‐１０４として動作させる循環系指標算出プログラム１１０、動脈圧情報１１１、生体情報１１２及び循環系指標情報１１３等を記憶する。
【００３５】
図３は、動脈圧情報１１１の構成の一例を示す概略図である。
【００３６】
動脈圧情報１１１は、動脈圧情報取得手段１００が動脈圧センサ２から取得した情報であって、値を取得した時刻ｔと、動脈圧ｐとを有する。
【００３７】
図４は、生体情報１１２の構成の一例を示す概略図である。
【００３８】
生体情報１１２は、生体情報取得手段１０１が生体センサ３から取得した情報であって、値を取得した時刻ｔと、心電図Ｖｈ、心拍数Ｈｒ、血圧Ｐｂ、体温Ｔｂ…等の生体情報とを有する。なお、図に示した生体情報は例示であってこれらに限定されるものではない。
【００３９】
図５は、循環系指標情報１１３の構成の一例を示す概略図である。
【００４０】
循環系指標情報１１３は、循環系指標算出手段１０２が動脈圧情報１１１から循環系の指標となる値として算出した情報であって、値を算出した時刻ｔと、循環系の指標である平均循環充満圧Ｐｓｆとを有する。
【００４１】
（循環系指標算出システムの動作）
次に、本実施の形態の作用を上記に説明した構成を前提とし、図１〜図１０、図１２を参照しつつ、（１）基本動作、（２）平均循環充満圧算出動作及び（３）表示処理動作に分けて説明する。
【００４２】
（１）基本動作
図１０は、循環系指標算出システム６の動作を説明するためのフローチャートである。
【００４３】
まず、医師５は、患者４が手術室に入室した後、患者４に対し動脈圧センサ２の一部としての針を動脈に刺入し、生体センサ３として心電図、非観血的血圧計、パルスオキシメータ、脳波計等を装着する。なお、動脈圧センサ２の一部としての針は動脈ラインを介して圧トランスデューサに接続されて動脈圧が計測される。
【００４４】
次に、必要に応じて、麻酔を行うため患者４の静脈に静脈留置針の刺入を行い、麻酔用の鎮静薬ポンプ及び鎮痛薬ポンプをそれぞれ接続する。また、人工呼吸器を患者４に取り付ける。
【００４５】
その後、循環系指標算出装置１の動脈圧情報取得手段１００は、通信部１４及び生体情報モニタ３Ａを介して動脈圧センサ２から動脈圧に比例した信号を取得して（Ｓ１）動脈圧の数値に変換し、図３に示したように、時刻情報とともに動脈圧情報１１１として記憶部１１に記録する。
【００４６】
また、循環系指標算出装置１の生体情報取得手段１０１は、通信部１４及び生体情報モニタ３Ａを介して生体センサ３から心電図や血圧計、パルスオキシメータ等の情報を取得し（Ｓ２）、図４に示したように、時刻情報とともに生体情報１１２として記憶部１１に記録する。
【００４７】
次に、循環系指標算出手段１０２は、動脈圧情報１１１から循環系の指標となる値として平均循環充満圧を算出する。算出方法の具体的手順を以下に説明する（Ｓ３〜Ｓ７）。
【００４８】
（２）平均循環充満圧算出動作
図６は、循環系指標算出手段１０２の動作を説明するための動脈圧波形及びその微分値を示すグラフ図である。
【００４９】
まず、循環系指標算出手段１０２は、図６上段に示すように、横軸を時刻、縦軸を動脈圧として動脈圧情報１１１を描いた波形とするとともに、図６下段に示すように、当該動脈圧波形を時間で微分した波形とする（Ｓ３）。
【００５０】
動脈圧波形のうち上行脚を時間で微分したものから正のピークが得られ（ｔ＝ｔｐ１、ｔｐ２…）、循環系指標算出手段１０２は、当該ピークの立ち上がりの始まり（ｔ＝ｔＢ０、ｔＢ１…）を一心拍の始まりとする。つまり、循環系指標算出手段１０２は、範囲（ｔＢ０〜ｔＢ１）、範囲（ｔＢ１〜ｔＢ２）…をそれぞれ一心拍の範囲として特定する（Ｓ４）。
【００５１】
次に、循環系指標算出手段１０２は、一心拍の範囲内において拡張期の始点を特定する（Ｓ５）。特定方法の具体的手順を以下に説明する。
【００５２】
図７は、循環系指標算出手段１０２の動作を説明するための一心拍の動脈圧波形を示すグラフ図である。
【００５３】
一心拍の動脈圧波形は、一般的に図７に示すような形状を有しており、重複切痕であるノッチの図中左の谷Ｐｎを境として収縮期（ｔＢｉ‐１〜ｔａｉ）と拡張期（ｔａｉ〜ｔＢｉ）に分けられる。循環系指標算出手段１０２は、ノッチの左の谷Ｐｎを特定するために衝撃波であるピークＰｐと重拍波であるピークＰｄ付近の双方に接する接線ｌｔを引き、当該接線ｌｔと距離ｄｍａｘが最大となる点を計算によって求めてこの点をノッチの左の谷Ｐｎ（ｔ＝ｔａｉ）とする。
【００５４】
次に、循環系指標算出手段１０２は、時刻（ｔ＝ｔａｉ）から予め定めた期間ｔｎｆだけ（例えば、ｔｎｆ＝１００ｍｓｅｃ）後の時刻（ｔ＝ｔＡｉ）を、回帰曲線を求める始点とする。つまり、循環系指標算出手段１０２は、範囲（ｔＡｉ〜ｔＢｉ）の下行脚の波形に対する回帰曲線を求める。なお、上述した下行脚の範囲を求める方法は一例であり、始点及び終点を定められれば方法は限定されるものではない。なお、ｔｎｆは定数であってもよいし、心拍数が高い場合には短く、心拍数が低い場合には長くして、脈拍数に合わせて変更するようにしてもよい。
【００５５】
循環系指標算出手段１０２は、拡張期の減衰波形に、一例として、以下の数式を用いてフィッティングする（Ｓ６）ことで、範囲（ｔＡｉ〜ｔＢｉ）の波形に対する回帰曲線を求める。
【数１】
2019062990.tif 000003
【００５６】
ここで、Ｐｅｓは、ｔ＝ｔＡｉにおける動脈圧であり、Ｐ０は収束値を示している。回帰曲線を求める手法の具体例として、最小二乗法が用いられ、最小二乗法アルゴリズムの一例としてレーベンバーグ・マーカート法を用いてフィッティングを行うことができる。
【００５７】
上記した数式（数１）は、血管の流路抵抗、血管のコンプライアンスを考慮し、血液循環を一時的に停止した場合に動脈圧が減衰して、血管内に均一に血液が分布したときに血管内圧がＰ０に収束するモデルから導出したものである。つまり、Ｐ０は血液循環を一時的に停止した場合の血管内圧の定常値である平均循環充満圧Ｐｓｆに対応する値である。
【００５８】
上記した数式（数１）は、一例として、以下の電気回路を用いたモデルから類推導出することができる。
【００５９】
図１２は、動脈圧の減衰を表すモデルとしての回路図である。
【００６０】
心臓からの拍出による圧力を電圧Ｖに、流量を電流Ｉに、血管のコンプライアンスをキャパシタの静電容量Ｃに、血管の流路抵抗を抵抗の抵抗値Ｒに、平均循環充満圧を定電圧Ｖ０の電源に置き換え、図１２に示すように、キャパシタと抵抗を並列に接続し、これらと電池を直列に接続した回路を疑似モデルとして想定すると、当該減衰回路の電圧Ｖはパルス応答の後に指数関数的に減衰して、値はＶ０に収束するものとなる。なお、当該関数が、循環系の圧力の特性を示すものとして妥当であることを示す根拠として、次に説明する論文で説明された結果との相関が挙げられる。
【００６１】
つまり、上記した数式（数１）を一心拍毎にフィッティングすることにより得られる収束値Ｐ０が以下に示す文献の平均循環充満圧の値と良好に相関することを発明者らは確認した（Ｊａｃｉｎｔａ  Ｊ．  Ｍａａｓ，  ｅｔ  ａｌ．、“Ｂｅｄｓｉｄｅ  Ａｓｓｅｓｓｍｅｎｔ  ｏｆ  Ｔｏｔａｌ  Ｓｙｓｔｅｍｉｃ  Ｖａｓｃｕｌａｒ  Ｃｏｍｐｌｉａｎｃｅ，  Ｓｔｒｅｓｓｅｄ  Ｖｏｌｕｍｅ，  ａｎｄ  Ｃａｒｄｉａｃ  Ｆｕｎｃｔｉｏｎ  Ｃｕｒｖｅｓ  ｉｎ  Ｉｎｔｅｎｓｉｖｅ  Ｃａｒｅ  Ｕｎｉｔ  Ｐａｔｉｅｎｔｓ”、Ａｎｅｓｔｈｅｓｉａ  ＆  Ａｎａｌｇｅｓｉａ、２０１２年、１０月、Ｖｏｌｕｍｅ  １１５、Ｉｓｓｕｅ  ４、ｐ８８０‐８８７参照。）。当該文献によると、上腕を駆血して橈骨動脈圧と中心静脈圧とを測定すると、橈骨動脈圧が減衰するとともに中心静脈圧が上昇する。橈骨動脈圧の減衰の結果、そして中心静脈圧の上昇の結果、両者は同じ値に収束し、当該同値に収束した値が平均循環充満圧Ｐｓｆである。当該平均循環充満圧Ｐｓｆと一心拍毎にフィッティングにより得られる収束値Ｐ０は良好に相関する。
【００６２】
なお、上記近似及び導出された数式は例示であって、血管の流路抵抗、血管のコンプライアンスを考慮したモデルから導出された数式であれば上記の数式（数１）に限られるものではなく、少なくとも収束値が求まればよく、適宜定数を乗じてもよいし、異なる数式を用いてもよい。
【００６３】
図８は、循環系指標算出手段１０２のフィッティング動作を説明するためのグラフ図である。
【００６４】
図８において、破線は連続する拍動の動脈圧波形を示しており、実線は範囲（ｔＡｉ〜ｔＢｉ）に対して上記数式をフィッティングしたものである。実線は、血管内圧の定常値Ｐ０まで伸びており、定常値Ｐ０は血液循環を一時停止した場合の平均循環充満圧Ｐｓｆに対応する。
【００６５】
循環系指標算出手段１０２は、上記数式を破線の動脈圧波形にフィッティングすることでパラメータ（Ｐｅｓ，Ｐ０，τ）を定め、Ｐ０から平均循環充満圧Ｐｓｆを算出するとともに、時定数τを算出する（Ｓ７）。循環系指標算出手段１０２は、図５に示したように、算出した平均循環充満圧Ｐｓｆを、対応する動脈圧波形の時刻とともに循環系指標情報１１３として記憶部１１に格納する。なお、対応する時刻は、一心拍の動脈圧波形の始点、終点又は始点と終点の間の任意の時刻を選択してよい。
【００６６】
次に、循環系指標算出手段１０２は、一心拍毎に得られた平均循環充満圧Ｐｓｆ及び時定数τを、当該心拍に対応する時刻の生体情報１１２と対応付ける（Ｓ８）。つまり、図４及び図５に示した、同時刻の生体情報１１２と循環系指標情報１１３とを対応付ける。
【００６７】
（３）表示処理動作
次に、表示処理手段１０３又は情報出力手段１０４は、平均循環充満圧Ｐｓｆ及びこれに対応付けられた生体情報１１２を表示部１２に表示する（Ｓ９）。
【００６８】
図９は、表示処理手段１０３又は情報出力手段１０４の表示動作の一例を説明するための概略図である。
【００６９】
表示画面１０３Ａは、表示処理手段１０３によって表示部１２に表示された表示内容又は情報出力手段１０４によって生体情報モニタ３Ａに表示された表示内容であり、メニュー、患者情報、時刻情報を表示するヘッダ表示領域１０３０と、生体情報１１２と生体情報１１２を経時的にグラフ表示した生体情報表示領域１０３１と、循環系指標情報１１３と循環系指標情報１１３を経時的にグラフ表示した循環系指標情報表示領域１０３２とを有し、生体情報表示領域１０３１のグラフ表示と循環系指標情報表示領域１０３２のグラフ表示とは同期して表示される。
【００７０】
医師５は、表示部１２又は生体情報モニタ３Ａに表示された上記表示画面１０３Ａの表示内容を確認し、患者４の循環系の状態を判断して輸液、輸血等の処置を適宜行う。
【００７１】
（実施の形態の効果）
上記した実施の形態によれば、患者４の動脈圧の時間変化のうち、心周期の拡張期の波形に対する回帰曲線を一心拍毎に求め、当該回帰曲線の収束値から平均循環充満圧Ｐｓｆを算出したため、測定された血管内圧の時間変化の波形から一心拍毎に前負荷の指標となる値を算出することができる。また、平均循環充満圧Ｐｓｆは、血管容量と血液量とから定まる指標であるため、術前の数値と術中又は術後の数値とを比較することで、患者の血管容量及び血液量、並びに患者の血管容量及び血液量の経過を監視することができる。
【００７２】
また、従来の平均循環充満圧を測定する方法は駆血が必要で、かつ、値が求まるまで少なくとも数十秒を要するため大きな侵襲を伴うものであったが、本願発明は少なくとも駆血を必要としないため上記侵襲を低減することができる。また、従来の平均循環充満圧は、測定までに数十秒を要するものであるから複数心拍における平均値と考えられるが、本願発明は一心拍毎に求めることができる。また、本願発明により求められた複数心拍の平均循環充満圧の平均をとることで精度を向上することができる。
【００７３】
また、従来は輸液の反応性を評価するために一回拍出量変動（ＳＶＶ：Ｓｔｒｏｋｅ  Ｖｏｌｕｍｅ  Ｖａｒｉａｔｉｏｎ）、脈圧変動（ＰＰＶ：Ｐｕｌｓｅ  Ｐｒｅｓｓｕｒｅ  Ｖａｒｉａｔｉｏｎ）を測定していた。ＳＶＶは動脈圧ラインから得られる圧波形に基づいて心拍出量を算定し、その一回心拍出量の呼吸性変動を変化率で表すものである。ＰＰＶは、動脈圧ラインから得られる脈圧の呼吸性変動の変化率を算出するものである。これらの値はどのくらいの輸液を行ってよいか評価するためのものであり、ＰＰＶ及びＳＶＶのいずれも輸液により反応のある数値であるが、呼気時、吸気時により変動するものであり一心拍毎の数値をそれぞれ比較することができないため、患者の血管容量及び血液量、並びに患者の血管容量及び血液量の経過を監視するためには用いることができなかったが、本願で算出する平均循環充満圧Ｐｓｆに対応するＰ０は、一心拍毎に数値が求まるものであるため、術前の数値と術中又は術後の数値とを比較することで、患者の適正な血液量を管理し、患者の血管容量及び血液量の経過を監視することができる。
【００７４】
また、以下に説明するように、平均循環充満圧Ｐｓｆからさらに静脈還流抵抗の解析が可能となる。
【００７５】
図１１は、心拍出量及び静脈還流量と平均右心房圧との関係を示すグラフ図である。
【００７６】
心拍出量曲線及び整脈還流量曲線は、個体内変動及び個体間変動はあるものの概ね図１１に示すような形状となる。静脈還流曲線の横軸との交点が平均循環充満圧Ｐｓｆであり、心拍出量曲線と整脈還流量曲線との交点が循環平衡点Ｐｃｅであり、中心静脈圧Ｃｖｐをその圧力とするものである。中心静脈圧Ｃｖｐは、右心房内圧であり中心静脈ライン等から測定される。
【００７７】
図１１において、静脈還流量曲線の右下がりの傾きの逆数は静脈還流抵抗を示すため、ＰｓｆとＰｃｅの２点が定まれば静脈還流抵抗が求められる。従来、リアルタイムで平均循環充満圧Ｐｓｆが測定できなかったため、臨床では静脈還流抵抗を求められなかったが、本発明によりＰｓｆが算出可能となるため、臨床において静脈還流抵抗が算出可能となる。また、逐次、静脈還流抵抗を算出することで静脈還流抵抗の時間変化を測定することができる。
【００７８】
［他の実施の形態］
なお、本発明は、上記実施の形態に限定されず、本発明の趣旨を逸脱しない範囲で種々な変形が可能である。
【００７９】
例えば、上記実施の形態では大動脈の動脈圧、つまり、体循環の血管内圧に対応する動脈圧を測定したが、肺循環の血管内圧に対応する動脈圧、つまり、肺動脈圧を測定してもよく、数値の絶対値は動脈圧に比べて小さいものの、肺動脈圧からも平均循環充満圧を求めることができる。また、全体循環の血管内圧として大動脈圧及び肺動脈圧を総合して測定してもよく、大動脈圧及び肺動脈圧から平均循環充満圧を求めてもよい。
【００８０】
また、上記実施の形態では、動脈圧の下行脚に対して回帰曲線を求めたが、静脈圧の増加曲線に対して回帰曲線を求めて平均循環充満圧を求めてもよい。静脈圧の変動幅が小さいためノイズ対策等は必要となるものの、理論的には静脈圧からも平均循環充満圧が算出される。
【００８１】
また、上記実施の形態では、動脈圧波形の微分値から一心拍の範囲を定義したが、心電図の波形等の生体情報をはじめとする他の情報から一心拍の範囲を定義してもよい。
【００８２】
また、時定数τのリアルタイム表示、履歴表示、複数拍動の平均値表示、予測値の表示を表示部１２及び／又は生体情報モニタ３Ａにさらに表示するようにしてよいのはもちろん、時定数τから算出される他の生体情報のリアルタイム表示、履歴表示、複数拍動の平均値表示、予測値の表示を表示部１２及び／又は生体情報モニタ３Ａにさらに表示するようにしてもよい。
【００８３】
また、表示部１２、操作部１３、通信部１４は、循環系指標算出装置１の必須の構成でないことはもちろんであり、省略してもよいし、別装置としてもよい。同様に生体センサ３は、循環系指標算出システム６に必須の構成ではなく、省略してもよいし、別システムとして用意してもよい。
【００８４】
上記実施の形態では制御部１０の各手段１００〜１０４の機能をプログラムで実現したが、各手段の全て又は一部をＡＳＩＣ等のハードウエアによって実現してもよい。また、上記実施の形態で用いたプログラムをＣＤ‐ＲＯＭ等の記録媒体に記憶して提供してもよいし、インターネットを介して配信することで提供することもできる。また、クラウド上で動作するプログラムであってもよい。また、上記実施の形態で説明した上記動作の順序の入れ替え、削除、追加等は本発明の要旨を変更しない範囲内で可能である。
【符号の説明】
【００８５】
１        ：循環系指標算出装置
２        ：動脈圧センサ
３        ：生体センサ
３Ａ      ：生体情報モニタ
４        ：患者
５        ：医師
６        ：循環系指標算出システム
１０      ：制御部
１１      ：記憶部
１２      ：表示部
１３      ：操作部
１４      ：通信部
１００    ：動脈圧情報取得手段
１０１    ：生体情報取得手段
１０２    ：循環系指標算出手段
１０３    ：表示処理手段
１０３Ａ  ：表示画面
１０４    ：情報出力手段
１１０    ：循環系指標算出プログラム
１１１    ：動脈圧情報
１１２    ：生体情報
１１３    ：循環系指標情報
１０３０  ：ヘッダ表示領域
１０３１  ：生体情報表示領域
１０３２  ：循環系指標情報表示領域
【図１】
2019062990.tif 000004
【図２】
2019062990.tif 000005
【図３】
2019062990.tif 000006
【図４】
2019062990.tif 000007
【図５】
2019062990.tif 000008
【図６】
2019062990.tif 000009
【図７】
2019062990.tif 000010
【図８】
2019062990.tif 000011
【図９】
2019062990.tif 000012
【図１０】
2019062990.tif 000013
【図１１】
2019062990.tif 000014
【図１２】
2019062990.tif 000015