(57)【要約】      （修正有）【課題】個人の胸腔等の流体源から除去される流体の量を把握することが所望される治療の機器および方法を提供する。【解決手段】流体収集キャニスタと関連付けられた排液管の端部における少なくとも１つの圧力値と、収集キャニスタから離れた場所における排液管内の少なくとも１つの圧力値とを使用することによって、患者の胸腔から、排液管を通して排液される流体の体積および／または流率を含む、患者パラメータを判定するデバイス、コンピュータ可読プログラム、および方法。圧力値は、非線形ソルバを用いて処理され、患者パラメータを判定する。少なくとも１つの第１の圧力値は、複数の圧力値を含む。【選択図】図３
2018140200.tif 000002
【特許請求の範囲】
【請求項１】
本明細書に記載の発明。
【発明の詳細な説明】
【技術分野】
【０００１】
（関連出願）
本願は、２０１３年３月１５日に出願された、米国仮特許出願第６１／７９４，７１２号、および、２０１４年３月１４日に出願された、米国非仮特許出願第１４／２１２，３３６号に対する優先権を主張するものであり、これらの全体の内容は、参照により本明細書中に援用される。
【０００２】
（連邦政府の補助金に関する記載）
本発明の実施形態は、連邦政府の補助金または助成金を得て着想または実践化されていない。
【０００３】
本発明の実施形態は、個人の胸腔等の流体源から除去される流体の量を把握することが所望される治療の機器および方法を対象とする。
【背景技術】
【０００４】
罹患個人の胸腔から流体を除去することが望ましい、いくつかの疾患状態および傷害が存在する。臨床医は、治癒プロセスの進行、疾患の進行、およびそのような排液をもたらすために使用されるデバイスの機能を判定するために、除去される流体の量を監視する。典型的には、患者は、胸腔の中に挿入される排液管を有し、管は、胸部排液ユニット（ＣＤＵ）として知られる特殊な閉鎖されたコンテナの中に排液させる。密閉されたコンテナ内の液体を測定するための種々の手段が存在し、それぞれ、観察または測定を行う臨床医の時間、あるいは測定を行い、測定データを臨床医または病院情報システムに知らせるためのセンサおよび信号伝達デバイスを提供する際の機器コストのいずれかにおけるコストと関連付けられる。
【０００５】
例えば、限定ではないが、いくつかのデバイスは、コンテナの重量または重量変化を測定するための１つまたはそれを上回る歪みゲージを特徴とする。歪みゲージの正確度が増加するにつれて、ゲージは、より高価となり、より損傷を受けやすくなる。歪みゲージは、特に、容易に損傷され、よりロバストであるように設計される場合、重量変化に対してより低い分解能を有する。
【０００６】
ＣＤＵに関して患者を監視する臨床スタッフは、患者の流体排液体積レベルを記録することが要求されるが、注意深く対処しなければならず、時間がかかり、したがって、コストがかかるため、そうでなければ望ましくあり得るほど頻繁に、測定またはチェックを行わない。排液デバイスは、通常、床に置かれ、臨床医が、流体レベルを読み取るために、手を伸ばし、膝をつき、かつ例えば、ＣＤＵが、使用の間、偶発的にひっくり返される場合、誤差を調節することを忘れないようにすることを要求する。
【発明の概要】
【課題を解決するための手段】
【０００７】
本発明の実施形態は、第１の端部および第２の端部を有する排液管を有する患者において、少なくとも１つの患者パラメータを判定するための方法、コンピュータプログラム、およびデバイスを対象にし、第１の端部は、患者の胸腔の中に嵌合され、第２の端部は、収集容器内に受容される。本明細書で使用されるように、用語「患者パラメータ」は、一例として、限定ではないが、胸腔から除去される流体（ガスおよび／または液体）の体積、除去される流体の流率、除去される流体の液体部分を保持するコンテナの現在の体積、濾過デバイスおよび管の現在の機能、ならびに患者呼吸ステータスを含む、排液管と関連付けられた患者パラメータを指す。
【０００８】
本方法、コンピュータプログラム、およびデバイスの実施形態は、収集容器および排液管のうちの少なくとも１つが、（ｉ）収集容器の近位の場所における排液管および（ｉｉ）収集容器のうちの少なくとも１つ内の圧力を感知する第１の圧力センサと、収集容器の遠位の場所における該排液管内の圧力を感知する、第２の圧力センサとを有することを特徴とする。すなわち、第１の圧力センサは、収集容器内部または収集容器に近接する排液管内のいずれかに位置し、収集容器と近接して関連付けられる。第２の圧力センサは、収集容器から離れた、患者により近位の排液管内に位置付けられる。
【０００９】
本方法は、第１の圧力センサからの少なくとも第１の圧力値および第２の圧力センサからの少なくとも第１の圧力値を得るステップを含む。本方法はさらに、非線形ソルバを通して、圧力値によって少なくとも１つの患者パラメータを計算するステップを含む。
【００１０】
本方法の一実施形態は、各圧力値と時間を関連付けるステップと、複数の連続圧力値を経時的に得るステップとを特徴とする。連続圧力値は、好ましくは、タイミングが図られた間隔において、時間的に分離される、すなわち、周期的である。
【００１１】
本方法の実施形態は、コンピュータ手段による実施のために非常に好適である。本明細書で使用されるように、用語「コンピュータ手段」は、メインフレームコンピュータ、サーバ、従来のデスクトップおよびラップトップコンピュータ、マイクロコントローラ、Ｒａｓｐｂｅｒｒｙ  Ｐｉ等の車載コンピュータ、ならびに携帯電話およびタブレットデバイス等のハンドヘルド通信デバイス、インターネットおよび／またはワイヤレス通信デバイスを通してリンクされるオフサイトコンピュータ、および同等物上に担持されるＣＰＵ等、デバイスの内部にある、またはデバイスと別個かつ別々に担持される、コンピュータまたは中央処理ユニット（ＣＰＵ）を指す。好ましくは、第１の圧力値および連続圧力値は、コンピュータメモリ内に記憶される。本明細書で使用されるように、用語「コンピュータメモリ」は、ＣＰＵと信号通信し、そのようなＣＰＵによって処理されることができる、チップ、ディスク、メモリドライブ、ハードドライブ、フラッシュドライブ、および同等物等のデータ記憶を指す。例えば、ＣＰＵは、最初に、メモリからの圧力値および連続圧力値にアクセスし、そのようなデータを編成するステップ、少なくとも１つの患者パラメータを計算するステップ、および１つの患者パラメータを表示する出力を生成するステップを行う。
【００１２】
本明細書で使用されるように、用語「非線形ソルバ」は、仮説／モデル数学的コンピュータプログラムの意味における機械学習プロセスを指す。非線形ソルバは、既知の入力および出力を用いて訓練され、訓練された非線形ソルバは、第１の圧力値および第２の圧力値ならびに全ての他の連続圧力値を受信し、仮説／モデルを通してそのような値を処理し、患者パラメータを得るようにエクスポートまたはプログラムされる。
実施形態は、非線形ソルバが、サポートベクトルマシンアプローチ、人工神経ネットワーク、遺伝子アルゴリズム、および遺伝子プログラミングから成る群から選択される、１つまたはそれを上回るアプローチにおいて実装されることを特徴とする。詳細な議論は、線形回帰仮説および機械学習プロセスの形態における非線形方程式を説明することになるであろう。
【００１３】
例えば、１つの非線形方程式は、以下の線形回帰仮説である。
【数１】
2018140200.tif 000003
前述で使用されるように、
Ｔ＝代表的現在の時間（ｔ＝０）、すなわち、データ後の固定時間点であり、
ｔ＝整数時間サンプルであり、０は、最新であり、Ｎ−は、最も古く、
ｆ（ｘ，ｙ）＝ｘおよびｙに依存する未知の関数であり、
ｈ（Ｔ）＝流体キャニスタ内の流体の体積であり、
Ｐｖ（ｔ）＝通気時の時間の関数としての圧力であり、
Ｐｃ（ｔ）＝時間の関数として患者のＣＤＵ側における圧力であり、
Ｎ＝計算のために要求される時間サンプル数であり、
Ｔｈｅｔａ＝以下に列挙される全Ｔｈｅｔａのセットであり、
Ｔｈｅｔａｂｉａｓ＝仮説のための任意のＤＣバイアス定数であり、
Ｔｈｅｔａｃ，ｔ＝時間ｔ時のＣＤＵ圧力に関する任意のパラメータ係数であり、
Ｔｈｅｔａｖ，ｔ＝時間ｔ時の通気圧力に関する任意のパラメータ係数である。
【００１４】
線形回帰方程式は、以下
【数２】
2018140200.tif 000004
のように表される、「ｍ」値の訓練セットを用いて表される、コスト関数Ｊを有する。
【００１５】
未知の「Ｔｈｅｔａ」値は、以下
【数３】
2018140200.tif 000005
のように、「ｍ」個のアイテムの訓練セットを用いて判定される。
【００１６】
例えば、限定ではないが、非線形方程式およびその導関数は、勾配降下法、粒子群最適化、および人工蜂コロニーから成る群から選定される、１つまたはそれを上回る最適化アルゴリズムによって解かれる。
１つの方法は、以下
【数４】
2018140200.tif 000006
のように導出されるＴｈｅｔａを特徴とする。
Ｘｘを伴うＪ（Ｔｈｅｔａ）に関する偏導関数は、以下
【数５】
2018140200.tif 000007
のように表される、Ｔｈｅｔａｘと関連付けられたＢｉａｓ、Ｐｖ、またはＰｃパラメータである。
【００１７】
本発明の特徴を具現化する、さらなる非線形方程式は、以下
【数６】
2018140200.tif 000008
のように表される。
式中、
Ｔ＝代表的現在の時間（ｔ＝０）、すなわち、データ後の固定時間点であり、
ｔ＝整数時間サンプルであり、０は、最新であり、Ｎ−は、最も古く、
ｆ（ｘ，ｙ）＝ｘおよびｙに依存する未知の関数であり、
ｈ（Ｔ）＝流体キャニスタ内の流体の体積であり、
Ｐｖ（ｔ）＝通気時の時間の関数としての圧力であり、
Ｐｃ（ｔ）＝時間の関数として患者のＣＤＵ側における圧力であり、
Ｎ＝計算のために要求される時間サンプル数であり、
Ｔｈｅｔａ＝以下に列挙される全Ｔｈｅｔａのセットであり、
Ｔｈｅｔａｂｉａｓ＝仮説のための任意のＤＣバイアス定数であり、
Ｔｈｅｔａｃ，ｔ＝時間ｔ時のＣＤＵ圧力に関する任意のパラメータ係数であり、
Ｔｈｅｔａｖ，ｔ＝時間ｔ時の通気圧力に関する任意のパラメータ係数である。
【００１８】
本発明のさらなる実施形態は、少なくとも１つの患者パラメータを判定するためのコンピュータ可読プログラムを対象とする。コンピュータ可読プログラムは、第１の圧力センサからの少なくとも１つの第１の圧力値および第２の圧力センサからの少なくとも１つの第１の圧力値を得るために、メモリを有するコンピュータ処理ユニットのための命令を備える。コンピュータ可読プログラムはさらに、コンピュータ手段に、非線形ソルバを通して、圧力値によって少なくとも１つの患者パラメータを計算するように命令する。
【００１９】
コンピュータ可読プログラムは、コンピュータ手段を動作させる、またはその上で起動し、例えば、タイミングが図られた間隔において、経時的に、少なくとも１つの第１の圧力値および任意の連続圧力値を受信し、そのようなデータをメモリ内に記憶する。データは、表示のための１つまたはそれを上回る患者パラメータを計算するために編成および使用される。これは、前述の非線形方程式および方法に従って処理される。
【００２０】
本発明のさらなる実施形態は、少なくとも１つの患者パラメータを判定するためのデバイスを対象とする。本デバイスは、（ｉ）収集容器の近位の場所における排液管または（ｉｉ）収集容器自体のうちの少なくとも１つ内の圧力を感知するように構築および配列される、第１の圧力センサと、収集容器の遠位の場所における該排液管内の圧力を感知するように構築および配列される、第２の圧力センサとを備える。第１の圧力センサおよび第２の圧力センサは、第１の圧力センサからの第１の圧力値および第２の圧力センサからの第１の圧力値を得るために、メモリを有するコンピューティング手段と信号通信するように設定されるために構築および配列される。コンピュータ手段は、非線形ソルバを通して、圧力値によって少なくとも１つの患者パラメータを計算する。
【００２１】
好ましいデバイスは、コンピュータ手段を備える。コンピュータ手段は、好ましくは、例えば、限定ではないが、タイミングが図られた間隔において、経時的に、複数の第１の圧力値および任意の連続圧力値を受信し、これは、メモリ内に記憶され、１つの患者パラメータを計算するために編成および処理される。
【００２２】
データは、表示のために、１つまたはそれを上回る患者パラメータを計算するように編成および使用される。これは、前述の非線形方程式および方法に従って処理され、以下の発明を実施するための形態により完全に説明される。
【００２３】
本発明の方法、コンピュータ可読プログラム、およびデバイスは、患者の胸腔から除去される流体の監視を可能にする。そのような流体は、ガス、水、血清、血液、血漿、リンパ液、膿、治療用注入流体、およびそれらの混合物等の任意の生物学的流体を含む。用語「ガス」は、酸素、窒素、ヘリウム、および酸化窒素等の治療用目的のために投与される空気およびガスを指すために使用される。本方法、コンピュータ可読プログラム、およびデバイスの実施形態は、ガス体積および液体体積ならびに流体の対応する流率の計算を可能にする。
【００２４】
これらならびに他の特徴および利点は、以下に概略される図を閲覧し、以下の発明を実施するための形態を検証することによって、当業者に明白となるであろう。
本発明は、例えば、以下を提供する。
（項目１）
第１の端部および第２の端部を有する排液管を有する患者において、少なくとも１つの患者パラメータを判定するための方法であって、前記第１の端部は、前記患者の胸腔の中に嵌合され、第２の端部は、収集容器と流体連通し、前記収集容器および前記排液管のうちの少なくとも１つは、（ｉ）前記収集容器の近位の場所における前記排液管および（ｉｉ）前記収集容器のうちの少なくとも１つ内の圧力を感知する、第１の圧力センサと、前記収集容器の遠位の場所における前記排液管内の圧力を感知する、第２の圧力センサとを有し、
前記方法は、
ａ．前記第１の圧力センサからの少なくとも第１の圧力値および前記第２の圧力センサからの少なくとも１つの第１の圧力値を得るステップと、
ｂ．非線形ソルバを通して、前記圧力値によって少なくとも１つの患者パラメータを計算するステップと、
を含む、方法。
（項目２）
前記少なくとも１つの第１の圧力値は、複数の圧力値を含む、項目１に記載の方法。
（項目３）
前記複数の圧力値を含む圧力値は、タイミングが図られた間隔におけるものである、項目２に記載の方法。
（項目４）
前記圧力値は、弁ステータス、時間、患者吸入ステータス、および管ステータスから成る群から選択される、１つまたはそれを上回る事象、ステータス、または値と関連付けられる、項目２−３のいずれかに記載の方法。
（項目５）
前記第１の圧力値および連続圧力値は、コンピュータメモリ内に記憶される、項目２−４のいずれかに記載の方法。
（項目６）
前記少なくとも１つの患者パラメータを計算するステップは、コンピュータによって行われる、項目１−５のいずれかに記載の方法。
（項目７）
前記少なくとも１つの患者パラメータを計算するステップは、コンピュータ手段によって行われる、項目１−６のいずれかに記載の方法。
（項目８）
前記コンピュータ手段は、前記少なくとも１つの患者パラメータを表示する出力を生成する、項目７のいずれかに記載の方法。
（項目９）
前記コンピュータ手段は、非線形方程式によって、前記少なくとも１つの患者パラメータを計算する、項目７−８のいずれかに記載の方法。
（項目１０）
前記非線形方程式は、サポートベクトルマシンアプローチ、人工神経ネットワーク、遺伝子アルゴリズム、および遺伝子プログラミングから成る群から選択される、１つまたはそれを上回るアプローチにおいて実装される、項目９に記載の方法。
（項目１１）
前記非線形方程式は、以下の線形回帰仮説、
【化１】
2018140200.tif 000009
であり、
Ｔ＝代表的現在の時間（ｔ＝０）、すなわち、データ後の固定時間点であり、
ｔ＝整数時間サンプルであり、０は、最新であり、Ｎ−は、最も古く、
ｆ（ｘ，ｙ）＝ｘおよびｙに依存する未知の関数であり、
ｈ（Ｔ）＝流体キャニスタ内の流体の体積であり、
Ｐｖ（ｔ）＝通気時の時間の関数としての圧力であり、
Ｐｃ（ｔ）＝時間の関数として患者のＣＤＵ側における圧力であり、
Ｎ＝計算のために要求される時間サンプル数であり、
Ｔｈｅｔａ＝以下に列挙される全Ｔｈｅｔａのセットであり、
Ｔｈｅｔａｂｉａｓ＝仮説のための任意のＤＣバイアス定数であり、
Ｔｈｅｔａｃ，ｔ＝時間ｔ時のＣＤＵ圧力に関する任意のパラメータ係数であり、
Ｔｈｅｔａｖ，ｔ＝時間ｔ時の通気圧力に関する任意のパラメータ係数である、
項目９に記載の方法。
（項目１２）
前記線形回帰方程式は、以下
【化２】
2018140200.tif 000010
のように表される、「ｍ」値の訓練セットを用いて表される、コスト関数Ｊを有する、項目１１に記載の方法。
（項目１３）
未知の「Ｔｈｅｔａ」値は、以下
【化３】
2018140200.tif 000011
のように、「ｍ」個のアイテムの訓練セットを用いて判定される、項目１２に記載の方法。
（項目１４）
前記非線形方程式およびその導関数は、勾配降下法、粒子群最適化、および人工蜂コロニーから成る群から選定される、１つまたはそれを上回る最適化アルゴリズムによって解かれる、項目１３に記載の方法。
（項目１５）
Ｔｈｅｔａは、以下
【化４】
2018140200.tif 000012
のように導出され、Ｘｘを伴うＪ（Ｔｈｅｔａ）に関する偏導関数は、以下
【化５】
2018140200.tif 000013
のように表される、Ｔｈｅｔａｘと関連付けられたＢｉａｓ、Ｐｖ、またはＰｃパラメータである、項目１２に記載の方法。
（項目１６）
前記非線形方程式は、以下
【化６】
2018140200.tif 000014
のように表され、
式中、
Ｔ＝代表的現在の時間（ｔ＝０）、すなわち、データ後の固定時間点であり、
ｔ＝整数時間サンプルであり、０は、最新であり、Ｎ−は、最も古く、
ｆ（ｘ，ｙ）＝ｘおよびｙに依存する未知の関数であり、
ｈ（Ｔ）＝流体キャニスタ内の流体の体積であり、
Ｐｖ（ｔ）＝通気時の時間の関数としての圧力であり、
Ｐｃ（ｔ）＝時間の関数として患者のＣＤＵ側における圧力であり、
Ｎ＝計算のために要求される時間サンプル数であり、
Ｔｈｅｔａ＝以下に列挙される全Ｔｈｅｔａのセットであり、
Ｔｈｅｔａｂｉａｓ＝仮説のための任意のＤＣバイアス定数であり、
Ｔｈｅｔａｃ，ｔ＝時間ｔ時のＣＤＵ圧力に関する任意のパラメータ係数であり、
Ｔｈｅｔａｖ，ｔ＝時間ｔ時の通気圧力に関する任意のパラメータ係数である、
項目９に記載の方法。
（項目１７）
前記流体は、ガス、水、血清、血液、血漿、リンパ液、膿、治療用注入流体、およびそれらの混合物から成る群のうちの任意の１つを含む、生物学的流体である、項目１−１６のいずれかに記載の方法。
（項目１８）
第１の端部および第２の端部を有する排液管を有する患者において、少なくとも１つの患者パラメータを判定するためのコンピュータ可読プログラムであって、前記第１の端部は、患者の胸腔の中に嵌合される、および／またはそれと流体連通を確立し、第２の端部は、収集容器と流体連通を確立し、前記収集容器および前記排液管のうちの少なくとも１つは、（ｉ）前記収集容器の近位の場所における前記排液管および（ｉｉ）前記収集容器うちの少なくとも１つ内の圧力を感知する第１の圧力センサと、前記収集容器の遠位の場所における前記排液管内の圧力を感知する、第２の圧力センサとを有し、前記コンピュータ可読プログラムは、
ａ．前記第１の圧力センサからの第１の圧力値および前記第２の圧力センサからの第１の圧力値を得ることと、
ｂ．非線形ソルバを通して、前記圧力値によって少なくとも１つの患者パラメータを計算することと
を行うために、メモリを有するコンピュータ手段のための命令を備える、コンピュータ可読プログラム。
（項目１９）
前記少なくとも１つの第１の圧力値は、経時的な複数の圧力値を含む、項目１８に記載のコンピュータ可読プログラム。
（項目２０）
前記複数の圧力値は、タイミングが図られた間隔で作成される、項目１９に記載のコンピュータ可読プログラム。
（項目２１）
前記複数の圧力値は、コンピュータメモリ内に記憶される、項目１９に記載のコンピュータ可読プログラム。
（項目２２）
前記少なくとも１つの患者パラメータを計算するステップは、コンピュータ手段によって行われる、項目１８−２１のいずれかに記載のコンピュータ可読プログラム。
（項目２３）
前記コンピュータ手段は、前記少なくとも１つの患者パラメータを表示する出力を生成する、項目２２のいずれかに記載のコンピュータ可読プログラム。
（項目２４）
前記コンピュータ手段は、非線形方程式によって、前記少なくとも１つの患者パラメータを計算する、項目２３のいずれかに記載のコンピュータ可読プログラム。
（項目２５）
前記非線形方程式は、サポートベクトルマシンアプローチ、人工神経ネットワーク、遺伝子アルゴリズム、および遺伝子プログラミングから成る群から選択される、１つまたはそれを上回るアプローチにおいて実装される、項目２４に記載のコンピュータ可読プログラム。
（項目２６）
前記非線形方程式は、以下の線形回帰仮説
【化７】
2018140200.tif 000015
であり、
Ｔ＝代表的現在の時間（ｔ＝０）、すなわち、データ後の固定時間点であり、
ｔ＝整数時間サンプルであり、０は、最新であり、Ｎ−は、最も古く、
ｆ（ｘ，ｙ）＝ｘおよびｙに依存する未知の関数であり、
ｈ（Ｔ）＝流体キャニスタ内の流体の体積であり、
Ｐｖ（ｔ）＝通気時の時間の関数としての圧力であり、
Ｐｃ（ｔ）＝時間の関数として患者のＣＤＵ側における圧力であり、
Ｎ＝計算のために要求される時間サンプル数であり、
Ｔｈｅｔａ＝以下に列挙される全Ｔｈｅｔａのセットであり、
Ｔｈｅｔａｂｉａｓ＝仮説のための任意のＤＣバイアス定数であり、
Ｔｈｅｔａｃ，ｔ＝時間ｔ時のＣＤＵ圧力に関する任意のパラメータ係数であり、
Ｔｈｅｔａｖ，ｔ＝時間ｔ時の通気圧力に関する任意のパラメータ係数である、
項目２４に記載のコンピュータ可読プログラム。
（項目２７）
前記線形回帰方程式は、以下
【化８】
2018140200.tif 000016
のように表される、「ｍ」値の訓練セットを用いて表される、コスト関数Ｊを有する、項目２６に記載のコンピュータ可読プログラム。
（項目２８）
未知の「Ｔｈｅｔａ」値は、以下
【化９】
2018140200.tif 000017
のように、「ｍ」個のアイテムの訓練セットを用いて判定される、項目２７に記載のコンピュータ可読プログラム。
（項目２９）
前記非線形方程式およびその導関数は、勾配降下法、粒子群最適化、および人工蜂コロニーアプローチによって解かれる、項目２６に記載のコンピュータ可読プログラム。
（項目３０）
Ｔｈｅｔａは、以下
【化１０】
2018140200.tif 000018
のように導出され、
Ｘｘを伴うＪ（Ｔｈｅｔａ）に関する偏導関数は、以下
【化１１】
2018140200.tif 000019
のように表される、Ｔｈｅｔａｘと関連付けられたＢｉａｓ、Ｐｖ、またはＰｃパラメータである、項目２７に記載のコンピュータ可読プログラム。
（項目３１）
前記非線形方程式は、以下
【化１２】
2018140200.tif 000020
のように表され、式中、
Ｔ＝代表的現在の時間（ｔ＝０）、すなわち、データ後の固定時間点であり、
ｔ＝整数時間サンプルであり、０は、最新であり、Ｎ−は、最も古く、
ｆ（ｘ，ｙ）＝ｘおよびｙに依存する未知の関数であり、
ｈ（Ｔ）＝流体キャニスタ内の流体の体積であり、
Ｐｖ（ｔ）＝通気時の時間の関数としての圧力であり、
Ｐｃ（ｔ）＝時間の関数として患者のＣＤＵ側における圧力であり、
Ｎ＝計算のために要求される時間サンプル数であり、
Ｔｈｅｔａ＝以下に列挙される全Ｔｈｅｔａのセットであり、
Ｔｈｅｔａｂｉａｓ＝仮説のための任意のＤＣバイアス定数であり、
Ｔｈｅｔａｃ，ｔ＝時間ｔ時のＣＤＵ圧力に関する任意のパラメータ係数であり、
Ｔｈｅｔａｖ，ｔ＝時間ｔ時の通気圧力に関する任意のパラメータ係数である、
項目２４に記載のコンピュータ可読プログラム。
（項目３２）
前記流体は、ガス、水、血清、血液、血漿、リンパ液、膿、治療用注入流体、およびそれらの混合物から成る群のうちの任意の１つを含む、生物学的流体である、項目１８−３１のいずれかに記載のコンピュータ可読プログラム。
（項目３３）
第１の端部および第２の端部を有する排液管を有する患者において、少なくとも１つの患者パラメータを判定するためのデバイスであって、前記第１の端部は、前記患者の胸腔の中に嵌合される、および／またはそれと流体連通を確立し、第２の端部は、収集容器内に受容され、前記デバイスは、
ａ．（ｉ）前記収集容器の近位の場所における排液管および（ｉｉ）前記収集容器のうちの少なくとも１つ内の圧力を感知するように構築および配列される、第１の圧力センサと、
ｂ．前記収集容器の遠位の場所における前記排液管内の圧力を感知するように構築および配列される、第２の圧力センサと、
を備え、
ｃ．前記第１の圧力センサおよび前記第２の圧力センサは、
ｉ．前記第１の圧力センサからの第１の圧力値および前記第２の圧力センサからの第１の圧力値を得ることと、
ｉｉ．非線形ソルバを通して、前記圧力値によって少なくとも１つの患者パラメータを計算することと
を行うために、メモリを有するコンピューティング手段と信号通信するために構築および配列される、デバイス。
（項目３４）
前記少なくとも１つの第１の圧力値は、複数の圧力値を含む、項目３３に記載のデバイス。
（項目３５）
前記複数の圧力値は、タイミングが図られた間隔で作成される、項目３４に記載のデバイス。
（項目３６）
メモリを有するコンピュータ手段をさらに備える、項目３３−３５のいずれかに記載のデバイス。
（項目３７）
前記複数の圧力値は、コンピュータメモリ内に記憶される、項目３４−３６のいずれかに記載のデバイス。
（項目３８）
前記コンピュータ手段は、前記少なくとも１つの患者パラメータを表示する出力を生成する、項目３６−３７のいずれかに記載のデバイス。
（項目３９）
前記コンピュータ手段は、非線形方程式によって、前記少なくとも１つの患者パラメータを計算する、項目３６−３８のいずれかに記載のデバイス。
（項目４０）
前記非線形方程式は、サポートベクトルマシンアプローチ、人工神経ネットワーク、遺伝子アルゴリズム、および遺伝子プログラミングから成る群から選択される、１つまたはそれを上回るアプローチにおいて実装される、項目３９に記載のデバイス。
（項目４１）
前記非線形方程式は、以下の線形回帰仮説
【化１３】
2018140200.tif 000021
であり、
Ｔ＝代表的現在の時間（ｔ＝０）、すなわち、データ後の固定時間点であり、
ｔ＝整数時間サンプルであり、０は、最新であり、Ｎ−は、最も古く、
ｆ（ｘ，ｙ）＝ｘおよびｙに依存する未知の関数であり、
ｈ（Ｔ）＝流体キャニスタ内の流体の体積であり、
Ｐｖ（ｔ）＝通気時の時間の関数としての圧力であり、
Ｐｃ（ｔ）＝時間の関数として患者のＣＤＵ側における圧力であり、
Ｎ＝計算のために要求される時間サンプル数であり、
Ｔｈｅｔａ＝以下に列挙される全Ｔｈｅｔａのセットであり、
Ｔｈｅｔａｂｉａｓ＝仮説のための任意のＤＣバイアス定数であり、
Ｔｈｅｔａｃ，ｔ＝時間ｔ時のＣＤＵ圧力に関する任意のパラメータ係数であり、
Ｔｈｅｔａｖ，ｔ＝時間ｔ時の通気圧力に関する任意のパラメータ係数である、
項目３９に記載のデバイス。
（項目４２）
前記線形回帰方程式は、以下
【化１４】
2018140200.tif 000022
のように表される、「ｍ」値の訓練セットを用いて表される、コスト関数Ｊを有する、項目４１に記載のデバイス。
（項目４３）
未知の「Ｔｈｅｔａ」値は、以下
【化１５】
2018140200.tif 000023
のように、「ｍ」個のアイテムの訓練セットを用いて判定される、項目４２に記載のデバイス。
（項目４４）
前記非線形方程式およびその導関数は、勾配降下法、粒子群最適化、および人工蜂コロニーアプローチによって解かれる、項目４１に記載のデバイス。
（項目４５）
Ｔｈｅｔａは、以下
【化１６】
2018140200.tif 000024
のように導出され、Ｘｘを伴うＪ（Ｔｈｅｔａ）に関する偏導関数は、以下
【化１７】
2018140200.tif 000025
のように表される、Ｔｈｅｔａｘと関連付けられたＢｉａｓ、Ｐｖ、またはＰｃパラメータである、項目４３に記載のデバイス。
（項目４６）
前記非線形方程式は、以下
【化１８】
2018140200.tif 000026
のように表され、
式中、
Ｔ＝代表的現在の時間（ｔ＝０）、すなわち、データ後の固定時間点であり、
ｔ＝整数時間サンプルであり、０は、最新であり、Ｎ−は、最も古く、
ｆ（ｘ，ｙ）＝ｘおよびｙに依存する未知の関数であり、
ｈ（Ｔ）＝流体キャニスタ内の流体の体積であり、
Ｐｖ（ｔ）＝通気時の時間の関数としての圧力であり、
Ｐｃ（ｔ）＝時間の関数として患者のＣＤＵ側における圧力であり、
Ｎ＝計算のために要求される時間サンプル数であり、
Ｔｈｅｔａ＝以下に列挙される全Ｔｈｅｔａのセットであり、
Ｔｈｅｔａｂｉａｓ＝仮説のための任意のＤＣバイアス定数であり、
Ｔｈｅｔａｃ，ｔ＝時間ｔ時のＣＤＵ圧力に関する任意のパラメータ係数であり、
Ｔｈｅｔａｖ，ｔ＝時間ｔ時の通気圧力に関する任意のパラメータ係数である、
項目３９に記載のデバイス。
（項目４７）
前記流体は、ガス、水、血清、血液、血漿、リンパ液、膿、治療用注入流体、およびそれらの混合物から成る群のうちの任意の１つを含む、生物学的流体である、項目３３−４６のいずれかに記載のデバイス。
（項目４８）
収集容器をさらに備える、項目３３−４７のいずれかに記載のデバイス。
（項目４９）
排液管をさらに備える、項目３３−４７のいずれかに記載のデバイス。
（項目５０）
前記コンピュータ手段は、体積データ中の誤差、流動データ中の誤差、フィルタ詰まり、呼吸ステータス、排液管詰まり、ライン内の液体、および過剰流体移動から成る群から選択される、少なくとも１つの付加的患者パラメータを監視する、項目３６−４９のいずれかに記載のデバイス。
（項目５１）
前記患者パラメータは、収集容器内に収集される液体流体の推定値である、項目１−１７のいずれかに記載の方法。
（項目５２）
前記患者パラメータは、収集容器内に収集される液体流体の推定値である、項目１８−３３のいずれかに記載のコンピュータ可読プログラム。
（項目５３）
前記患者パラメータは、収集容器内に収集される液体流体の推定値である、項目３３−５０のいずれかに記載のデバイス。
（項目５４）
前記非線形ソルバは、訓練される、項目１−１７および５１のいずれかに記載の方法。
（項目５５）
前記非線形ソルバは、訓練される、項目１８−３２および５２のいずれかに記載のコンピュータ可読プログラム。
（項目５６）
前記非線形ソルバは、訓練される、項目３３−５０および５３のいずれかに記載のデバイス。
【図面の簡単な説明】
【００２５】
【図１】図１は、例示的実施形態による、患者の胸腔から排液される流体の体積および／または流動を判定するための一例示的流体分析器を示す。
【００２６】
【図２】図２は、例示的実施形態による、さらに例示的詳細において図１の電子制御モジュールを示す、概略図である。
【００２７】
【図３】図３は、例示的実施形態による、時系列圧力値を収集および記憶するための一例示的方法を図示する、フロー図である。
【００２８】
【図４】図４は、例示的実施形態による、胸部排液流体の測定のために、図３の記憶された時系列圧力値を処理するための一例示的方法を図示する、フロー図である。
【００２９】
【図５】図５は、例示的実施形態による、図１の排液デバイスをモデル化する、一例示的電子回路を図示する、概略図である。
【００３０】
【図６】図６は、例示的実施形態による、機械学習プロセスを描写する。
【発明を実施するための形態】
【００３１】
本発明の実施形態は、患者の胸腔等の外部源からの流体体積を測定するためのデバイス、システム、および方法に関して詳細に説明されるであろう。ここで、図１を参照すると、時として、本明細書では、流体分析器と称され、概して、番号１００によって指定される、本発明の特徴を具現化するデバイスが、描写される。
【００３２】
図１は、患者の胸腔ｌ１０等のソースから排液される流体の体積および／または流動を判定するための一例示的流体分析器１００を示す。図１の実施例では、流体分析器１００は、例えば、二重管腔カテーテル１０９等のカテーテルまたは管類を使用して、患者の胸腔１１０から排液された流体を収集するように動作する、排液デバイス１２０内で実装される。図１の実施例では、排液デバイス１２０は、真空源１１２に結合され、患者１１０からの流体排液を促進するように示される。
【００３３】
排液デバイス１２０は、カテーテル１０９を通して流体キャニスタ１０８の中への流体の排液を制御するように動作する、電子制御モジュール（ＥＣＭ）１０２を含むように示される。描写されるように、ＥＣＭ１０２は、デバイス１２０のユーザ（例えば、臨床スタッフ）と相互作用するためのユーザインターフェース１０６を含む。流体キャニスタ１０８は、例えば、患者の胸腔から排液される流体を収集する、胸部排液ユニット（ＣＤＵ）である。ＥＣＭ１０２は、（ａ）カテーテル１０９の通気経路の患者接続部近傍の圧力を感知する、通気圧力センサ１２６と、（ｂ）同一患者点から生じる、排液管および収集システムの端部における圧力を感知する、ＣＤＵ圧力センサ１２８とから圧力値を受信する。通気圧力センサ１２６は、特に、カテーテル１０９が二重管腔タイプではない場合、カテーテル１０９の排液パイプ上に位置してもよい。排液デバイス１２０および／またはカテーテル１０９とともに位置付けられる、付加的圧力センサも、本明細書の範囲から逸脱することなく、システム１００と併用されてもよい。
【００３４】
排液デバイス１２０は、本明細書の範囲から逸脱することなく、他の構成要素（例えば、弁、フィルタ、センサ、コンバータ、インターフェース等）を含んでもよい。排液デバイス１２０は、圧力センサ１２６および１２８からの時系列の圧力値を処理し、患者１１０から排液される流体体積および／または流体の流率を判定する。ＥＣＭ１０２は、限定ではないが、ＣＤＵフィルタ詰まり、患者呼吸、患者吸入、患者呼気、呼吸数／分、呼吸困難患者（例えば、咳、窒息、笑っている、泣いている等）、管詰まりおよび閉塞または詰まり度、ライン内の流体、ならびに患者からの流体ダンピング症状等の他の指標を判定してもよい。
【００３５】
図２は、さらなる例示的詳細においてＥＣＭ１０２を示す、概略図である。図１および２は、以下の説明とともに最も良く概観される。ＥＣＭ１０２は、ユーザインターフェース１０６、プロセッサ２０６、メモリ２０８、およびセンサインターフェース２１０を含むように示される。メモリ２０８は、プロセッサ２０６によって実行されると、患者１１０から排液される流体の体積および／または流率を判定するために、排液デバイス１２０内で機能性を実装する、機械可読命令を含む、ソフトウェア１０４を記憶するように例証的に示される。センサインターフェース２１０は、通気圧力センサ１２６およびＣＤＵ圧力センサ１２８の両方と結合し、アナログ／デジタルコンバータを含み、感知された圧力値をソフトウェア１０４に提供する。ソフトウェア１０４は、データ取得器２１２および計算機２１４を含む。データ取得器２１２は、圧力値を圧力センサ１２６および１２８から経時的に収集するように動作し、これらの値を時系列圧力データ２２０としてメモリ２０８内に記憶する。時系列圧力データ２２０は、センサインターフェース２１０を介して、圧力センサ１２６および１２８から周期的に読み取られる圧力値を含有する。計算機２１４は、時系列圧力データ２２０を処理し、キャニスタ１０８内の流体の体積を判定し、随意に、キャニスタ１０８の中への流体の流率を判定するように動作する。計算機２１４は、センサインターフェース２１０から読み取られた圧力値および定義された体積に基づいて、これらの圧力値と流体体積および／または流率との間の関係を自動的に学習する、学習アルゴリズムを含んでもよい。計算機２１４は、各圧力センサ１２６およびｌ２８において識別された圧力変動に基づいて、圧力特徴データ２２２を生成する。圧力特徴データ２２２は、以下、通気圧力特徴および／またはＣＤＵ圧力特徴と称され得る。計算機２１４はまた、キャニスタ１０８内の流体体積およびキャニスタ１０８の中への流体流率の一方または両方を含む、流体測定値２２４を生成する。データ取得器２１２および計算機２１４は、並行して動作してもよい。
【００３６】
図３は、時系列圧力値を読み取り、記憶するための一例示的方法３００を図示する、フロー図である。図４は、胸部排液流体の測定のための記憶された時系列圧力データを処理するための一例示的方法を図示する、フロー図である。方法３００および４００は、例えば、ＥＣＭ１０２のソフトウェア１０４内に実装されるコンピュータ上で起動する、例えば、別個の処理スレッドである。図１から４は、以下の説明とともに最も良く概観される。
【００３７】
ステップ３０２では、方法３００は、第１の圧力値を通気圧力センサから読み取る。ステップ３０２の一実施例では、センサインターフェース２１０を使用して、プロセッサ２０６によって実行される、ソフトウェア１０４は、第１の圧力値を通気圧力センサ１２６から読み取る。ステップ３０４では、方法３００は、第２の圧力値をＣＤＵ圧力センサから読み取る。ステップ３０４の一実施例では、プロセッサ２０６によって実行される、ソフトウェア１０４は、センサインターフェース２１０を使用して、第２の圧力値をＣＤＵ圧力センサ１２８から読み取る。
【００３８】
ステップ３０６では、方法３００は、ステップ３０２および３０４の圧力値を時系列圧力データとして記憶する。ステップ３０６の一実施例では、データ取得器２１２は、ステップ３０２および３０４の圧力値をメモリ２０８に時系列圧力データ２２０として記憶する。
【００３９】
方法３００は、周期的に繰り返し、圧力値を収集および記憶し、第１および第２の圧力値に基づいて、流体パラメータを計算する。
【００４０】
ステップ４０２では、方法４００は、時系列圧力データをメモリから読み取る。ステップ４０２の一実施例では、計算機２１４は、時系列圧力データ２２０をメモリ２０８から読み取る。ステップ４０４では、方法４００は、特徴を時系列圧力データから抽出する。ステップ４０４の一実施例では、計算機２１４は、圧力特徴を時系列圧力データ２２０から抽出し、これらの圧力特徴を圧力特徴データ２２２としてメモリ２０８内に記憶する。
【００４１】
ステップ４０６は、決定である。ステップ４０６では、方法４００は、ステップ４０４において判定された特徴が有効であるかどうかを判定し、これらの特徴が有効である場合、方法４００は、ステップ４０８に進み、そうでなければ、方法４００は、ステップ４０２に進む。
【００４２】
ステップ４０８では、方法４００は、非線形ソルバを通して、特徴を処理する。ステップ４０８の一実施例では、計算機２１４は、非線形方程式を通して、圧力特徴データ２２２を処理し、流体体積および流体流率の一方または両方を含む、流体測定値２２４を生成する。ステップ４０８の別の実施例では、計算機２１４は、事前に訓練されたサポートベクトルマシンを通して、圧力特徴データ２２２を処理する。ステップ４０８の別の実施例では、計算機２１４は、人工神経ネットワーク、遺伝子アルゴリズム、および遺伝子プログラミングのうちの１つまたはそれを上回るものを使用して、圧力特徴データ２２２を処理する。例えば、計算機２１４は、収集される訓練、妥当性確認、および試験データセットのうちの１つまたはそれを上回るものに基づいて、データを利用してもよい。
【００４３】
ステップ４１０は、決定である。ステップ４１０では、方法４００は、ステップ４０８において生成されたデータが有効であるかどうかを判定し、本データが有効である場合、方法４００は、ステップ４１２に進み、そうでなければ、方法４００は、ステップ４０２に進む。
【００４４】
ステップ４１２では、方法４００は、判定された体積および他の値を報告する。ステップ４１２の一実施例では、計算機２１４は、流体測定２２４をユーザインターフェース１０６に出力する。
【００４５】
方法４００は、周期的に繰り返し、時系列圧力データ２２０に基づいて、流体測定値２２４を判定および出力する。
【００４６】
図５は、排液デバイス１２０の異なる区分の抵抗、インダクタンス、および静電容量を表す、１つの電子アナログ５００を示す。したがって、例示的実施形態によると、電子アナログ５００は、排液デバイス１２０をモデル化する。しかしながら、排液デバイス１２０は、代替として、本明細書の範囲から逸脱することなく、他の様式でモデル化され得る。電子アナログ５００は、可変キャパシタによってモデル化される、流体キャニスタ１０８のみを除き、両ブランチに対して固定される。各圧力センサ１２６、１２８によって「被られる」患者および「電気ネットワーク」差によって生成される動的圧力信号を前提として、本タイミングは、流体キャニスタ１０８内の流体体積に応じて、非線形に変動するであろう。流体キャニスタ１０８内の流体体積の増加は、圧力タイミングに影響を及ぼす回路の支流内の空気容量を減少させるであろう。フィルタ（図５のＲｆ）は、システムの抵抗が変動（経時的に増加）し得る。本抵抗変動は、線形影響であって、経時的に徐々に生じ得る。
【００４７】
図６は、システムまたはデバイス１００のための機械学習プロセスを描写する。ステップ６００、６１０、６２０、および６４０は、訓練として特徴付けられる。ステップ６００は、訓練のための入力および出力を収集する。これらの入力および出力は、制御された設定において実験的に得られ、訓練ステップ６１０において、入力および出力のマトリクスを生成するであろう。訓練ステップ６１０は、承認または適切性に関して評価される、１つまたはそれを上回る出力を生成する。遺伝子アルゴリズムおよびプログラムは、ステップ６４０において、モデルへの随意の変更を伴う１つの値にコンパイルされる、訓練の間の所望の出力および全予期される出力を基にする適切性に基づいて、入力および出力を評価する。モデルが容認可能結果を生成する場合、モデルは、ステップ６３０において、デバイス１００にエクスポートされる。
【００４８】
モデルは、ステップ６５０において、デバイス１００の中に統合される。また、モデルは、ステップ６６０において、モデルのための入力パラメータを計算するために使用される。モデルは、ステップ６７０およびステップ６８０において、入力パラメータを処理するために使用され、デバイス１００は、出力結果に応じて作用する。
【００４９】
（例示的アプローチ）
ある実施形態では、非線形方程式仮説が、流体キャニスタ１０８内の流体の体積を計算するために使用される。固定モデルまたは動的および自己補正モデルのいずれかが、使用されてもよい。固定モデルは、典型的には、定数（重量または質量）が変化せず、再訓練を伴わずに、ソフトウェア１０４内でそのまま使用され得るように所望される。仮説は、固定モデルまたは動的および自己保定モデルのいずれに対しても同一であろうが、方程式は、モデルのために最適化される必要があり得る。
（仮説）
ＬＥＴ：
・Ｔ＝代表的現在の時間（ｔ＝０）、すなわち、データ後の固定時間点
・ｔ＝整数時間サンプルであり、０は、最新であり、Ｎ−は、最も古い
・ｆ（ｘ，ｙ）＝ｘおよびｙに依存する未知の関数
・Ｊ（Ｔｈｅｔａ）＝Ｔｈｅｔａに関する仮説対実際のコスト
・Ｖ（Ｔ）＝流体キャニスタ１０８内の流体の体積（実際）
・ｈ（Ｔ）＝流体キャニスタ１０８内の流体の体積（仮説）
・Ｐｖ（ｔ）＝例えば、圧力センサ１２６から判定される、通気時の時間の関数としての圧力（または、通気圧力特徴、および他のパラメータ）
・Ｐｃ（ｔ）＝例えば、圧力センサ１２８から判定される、時間の関数としての患者のＣＤＵ側における圧力（または、ＣＤＵ圧力特徴、および他のパラメータ）
・ｍ＝訓練セット内のサンプルの数
・Ｎ＝計算のために要求される時間サンプル数
・Ｔｈｅｔａ＝以下に列挙される全Ｔｈｅｔａのセット
・Ｔｈｅｔａｂｉａｓ＝仮説のための任意のＤＣバイアス定数
・Ｔｈｅｔａｃ，ｔ＝時間ｔ時のＣＤＵ圧力に関する任意のパラメータ係数
・Ｔｈｅｔａｖ，ｔ＝時間ｔ時の通気圧力に関する任意のパラメータ係数
【００５０】
Ｖ（Ｔ）を概算するために、以下の線形回帰仮説を立てる。
【数７】
2018140200.tif 000027
【００５１】
コスト関数Ｊは、以下
【数８】
2018140200.tif 000028
のように、「ｍ」個のアイテムの訓練セットとともに、示されるように表されることができる。
【００５２】
未知の「Ｔｈｅｔａ」値は、以下
【数９】
2018140200.tif 000029
のように、「ｍ」個のアイテムの訓練セットとともに、見出され得る。
【００５３】
本方程式が解法または修正され得る、いくつかの方法がある。方程式を解法または修正する可能性として考えられるアプローチの１つは、「勾配降下法」アプローチ／アルゴリズムの使用を通したものである。勾配降下法アプローチは、方程式のための解の良好な近似を提供する。代替方法として、粒子群最適化システム等のより最新かつより高速の最適化系が挙げられる。方程式は、人工蜂コロニー（ＡＢＣ）を使用して、低オーバーヘッドを伴って、高速で解かれることが可能である。別のアプローチは、「Ｔｈｅｔａ」を解くために、ＭａｔＬａｂ、Ｒ、またはＯｃｔａｖｅ等の外部数学パッケージを使用することである。これらの解法アプローチは、当技術分野において公知である。Ｔｈｅｔａのこれらの値は、例えば、ソフトウェア１０４の形態等において、ＥＣＭ１０２の中に入力または別様に組み込まれる。本アプローチは、典型的には、任意の静的「Ｔｈｅｔａ」群に対して十分である。
【００５４】
勾配降下法は、以下
【数１０】
2018140200.tif 000030
のように、「Ｔｈｅｔａ」を下方移動させるために、「Ｔｈｅｔａ」の観点からのコスト関数のための偏導関数を要求する。
【００５５】
Ｔｈｅｔａｘと関連付けられたＢｉａｓ、Ｐｖ、またはＰｃパラメータである、Ｘｘを伴うＪ（Ｔｈｅｔａ）に関する偏導関数は、以下
【数１１】
2018140200.tif 000031
のように表されることができる。
【００５６】
前述の方程式仮説が収束しない場合、２乗項、および潜在的に、さらに３乗項が、収束達成を支援するために追加されてもよい。付加的累乗項のそのような組み込みは、解が高バイアス条件（適合不足）を呈する場合、要求されるであろう。反対の問題は、あまりに多くのパラメータが存在する場合に生じ得、これは、高分散（過剰適合）と称される。過剰適合は、規則化およびより多くのサンプルの追加等、いくつかの技法を通して解決されることができる。
【００５７】
以下は、例えば、一次モデルが満足のゆくように収束しない場合に使用される、より高次の項を含む、二次方程式仮説である。
【数１２】
2018140200.tif 000032
【００５８】
前述の数学的モデルの正確度は、導出されるパラメータの使用によって改良され得る。例えば、最大値および２つの波形間の時間異が、圧力波形から導出されることができる。大規模な時系列にわたる分析ではなく、導出されるパラメータデータセットの使用は、コスト計算における「誤差」項を低減させるであろう。未加工データのパラメータ抽出を伴わない場合、より多くのサンプルが、誤差を排除するために処理される必要があるであろう。
【００５９】
方程式のさらなる複素解は、その固有モード関数またはＩＭＦへの時間信号の分解を通して得られる。固有モード関数（ＩＭＦ）は、複雑性／計算の増加を犠牲において、正確度を改善するために使用されることができる。最近の研究として、ヒルベルト・ホアン変換（ＨＨＴ）の形態において、Ｎｏｒｄｅｎ  Ｅ．  Ｈｕａｎｇによって行われたものがある。本アプローチは、経験モード分解（ＥＭＤ）およびヒルベルトスペクトル分析（ＨＡＳ）を使用して、正確度を改善するが、計算の増加を犠牲にするという短所を伴う。特定の実施形態では、前述の数学的方程式仮説が、試験され、ＨＨＴ技法は、数学的仮説に満足がいかない場合のみ、実装される。
【００６０】
排液デバイス１２０は、少なくとも部分的に、圧力センサ１２６、１２８からの圧力データを使用して、流体キャニスタ１０８内の流体の体積に加え、情報を判定可能である。加えて、ある実施形態はさらに、１つまたはそれを上回る弁センサ（図示せず）または弁のプログラムされたタイミングを含み、排液デバイス１２０内の１つまたはそれを上回る弁の状態を検出する。これらの実施形態では、ＥＣＭ１０２は、例えば、圧力センサ１２６、１２８からのデータとともに、弁センサからの弁データを使用し、付加的情報を判定する。ＥＣＭ１０２は、例えば、流体キャニスタ１０８の流体体積を判定するための前述のものと同様の技法等の線形回帰技法を使用して、本付加的情報を判定する。別の実施例として、いくつかの実施形態では、ＥＣＭ１０２は、ロジスティック回帰（分類、確率出力／クラス提示とも呼ばれる）を使用して、付加的情報を判定する。
【００６１】
前述のように、排液デバイス１２０は、流体キャニスタ内１０８の流体の概算体積を推定可能である。そのような推定の正確度は、流体キャニスタ１０８内のシステムの「空気容量」の変化に起因する、流体体積に反比例する。加えて、ある実施形態は、経時的にフィルタ変化を特性評価する情報を使用して、わずかな体積変化も検出可能である。排液デバイス１２０は、例えば、フィルタのダンピングに起因する、可変レジスタの抵抗変化を検出することによって、フィルタ変化を判定する。さらに、いくつかの実施形態は、流体キャニスタ１０８の体積がどのように経時的に変化し得るかに基づいて、破損縫目または縫合糸等、重度のダンピング症状患者１１０を検出可能である。
【００６２】
いくつかの実施形態は、排液デバイス１２０内における、患者圧力緩和弁１３０の動作からの弁データおよび圧力センサ１２８からのＣＤＵ圧力データを含む、線形回帰モデルを介して、流体キャニスタ１０８を通る空気流動（また、流体と見なされる）の判定が可能である。正確度は、例えば、システムからのフィードバックおよび現在の空気通過における変化を特性評価する情報を使用することによって増加される。さらに、いくつかの実施形態は、１回換気量を流体キャニスタ１０８に流入させ、解放前に、気圧に近似するように弁を保持し、さらに正確度を助長することが可能である。これらの実施形態では、排液デバイス１２０は、典型的には、療法に有意に影響を及ぼすことなく、より優れた臨床データを達成するために、一瞬だけ、圧力を保持する。安全限界が、このような場合、特徴に対して設定されるであろう。圧力時間減衰データは、通過された体積／ボーラスの概算を支援し、ある時間周期にわたって合計され、空気流率を報告するために使用されることができる。
【００６３】
呼吸数は、類似アプローチを使用して見出され得るが、大量の時間データセット（または、前述のデータセットからのサンプル）を要求し得る。
【００６４】
排液デバイス１２０のある実施形態は、１つまたはそれを上回るロジスティック回帰モデルを使用して、以下の事象のうちの１つまたはそれを上回るものの存在または存在確率を判定可能である。
・線形回帰からの体積データが有効である
・例を訓練するために十分に類似しない圧力信号
・報告される体積における可能性として考えられる誤差
・流体キャニスタ１０８のフィルタ詰まり（ＣＤＵを交換してください）
・患者が呼吸している
・患者が吸入している
・患者が呼気している
・前述の事象に基づいて計算される呼吸数／分
・呼吸困難患者等の患者状態（例えば、咳、窒息、笑っている、泣いている等）
・管詰まり（例えば、完全閉塞）
・管閉塞度
・ライン内の流体
・患者ダンピング流体
【００６５】
本明細書に開示されるシステムおよび方法のある実施形態は、以下の利点のうちの１つまたはそれを上回るものを達成する。
・より良好な患者回復時間。
・管ライン一掃デバイス（ＴＬＣ）が、排液デバイス１２０の一部として埋め込まれ、必要に応じて、排液管を一掃し、従属ループを検出するために、モニタ流体収集および率、空気流動および率を監視することができる。
・ＴＬＣは、従属ループおよび流体収集が存在するという高信頼度があるとき、アクティブ化されることができる。これは、不必要なＴＬＣ起動による療法の中断が少なく、バッテリ寿命を延長させるであろう。詰まりがなく、ＴＬＣをアクティブ化するとき、患者へのリスクが少ない。
・患者安全性は、より優れた監視を通して改良され、例えば、システムは、患者がシステムに適切に接続されているかどうか、および／または患者が呼吸しているかどうかを検出することができる。臨床上の支援は、このような極端な場合には、招集されることができる。
・システムは、排液管が閉塞されているか、または経路を遮断する流体を有するかどうかを判定し、流体妨害物に応じて、ＴＬＣを起動し、臨床スタッフにねじれた管および同等物をアラートすることができる。
・システムは、患者記録に記録するために、または研究において使用するために、または中央ナースステーションに送信するために、データ内のより多くの値を提供する。
・中央ナースステーションの中にネットワーク化されたシステムは、看護スタッフに、以下に関する情報を提供することができる。
・呼吸数（患者管の切断、睡眠または非睡眠関連無呼吸）
・患者の苦痛
・流体の体積
・空気流率
・管開存
・患者からの能動的ダンピング流体
【００６６】
ある実施形態では、本発明の方法、デバイス、およびソフトウェアは、データの訓練セット、次いで、ＥＣＭから記録された試験セットを捕捉後、実装される。データの訓練セットは、例えば、１０ｍｌ体積インクリメントにおける多くのサンプルを含有する。計算における高分散（補間されたサンプルが、予期される誤差より高くあり得る）は、付加的訓練サンプルならびにコスト関数および／またはパラメータセットへの修正を用いて対処される。
【００６７】
種々の例示的実施形態では、サンプルデータセットは、異なる時間／体積における、「不良」データを含有してもよい。「不良」データは、システムモデルが違反され、もはや十分な予測確実性を伴って作用しない測定と見なされるであろう。「不良」データ実施例の１つは、流体が排液ラインまたはライン内に取り付けられたクランプ内にトラップされるときに取得されたデータである。付加的「不良」実施例は、咳および無呼吸等、呼吸の急変を含み得る。別の実施例は、異なる体積における空気漏出を含み得る。
【００６８】
当業者は、データを分析し、「ｔｈｅｔａ」を見出し、次いで、要求に応じて、体積および他のパラメータを解くために使用される、Ｍａｔｌａｂ等の種々のツールおよびＯｃｔａｖｅと呼ばれるフリーツールの可用性を認識するであろう。
【００６９】
（ロジスティック回帰）
ある実施形態では、ロジスティック回帰が、入力の分類を生成するシステムを開発するために適用される（例えば、詰まりステータスに関するはい／いいえ、データ有効／データ無効、患者呼吸等）。ロジスティック回帰は、前述の線形回帰と非常に類似し、方程式仮説は、その出力が０〜１である、シグモイド関数によってラップされる。
【数１３】
2018140200.tif 000033
【００７０】
ロジスティック仮説によって表されるように、方程式は、以下の形態をとり得る。
【数１４】
2018140200.tif 000034
【００７１】
コスト関数は、極小値の回避を支援するために、凸解を保つように変更される。コスト関数の新しい形態では、ｙは、クラスに属する場合、１であって、そうでなければ、０である。
【数１５】
2018140200.tif 000035
【００７２】
偏導関数は、分類子ｙ（０または１）の項においてのみ類似する。
【数１６】
2018140200.tif 000036
【００７３】
分類は、例えば、Ｏｎｅ  ｖｓ  Ａｌｌ構造において行われる。すなわち、分類毎に、神経ネットワークが、訓練されるクラスのみ「クラス内」または１であるように訓練されるであろう。全ての他のサンプルは、「クラス外」または別様に０に設定されるであろう。これは、所望される分類毎に異なる「ｔｈｅｔａ」セットが存在するであろうことを意味する。
【００７４】
（サポートベクトルマシンアプローチ）
サポートベクトルマシンは、分類に非常に好適であって、かつ特殊である。全前述のロジスティック回帰を訓練されるサポートベクトルマシン（ＳＶＭ）に分割することによって、クラスが入力データ内に存在する確率を生成する、計算上効率的システムを得る。
【００７５】
ＳＶＭは、ＳＶＭが、正値例をその超平面（ｎ−次元空間）内に「記憶」するという点において、ロジスティック回帰と異なる。次いで、サンプルが学習された例のうちの１つにどれだけ近いかを報告する。カーネル選択は、有用であって、超平面内の既知の正値に関連して、サンプルの確率を定義する。「ｌｉｂＳＶＭ」と呼ばれるオープンソースかつ最適化されたＣ＋＋ライブラリが、種々の例示的実施形態に従って、クラス予測のために、プロトタイプ化するために使用される、および／または排液デバイス１２０の中に統合されてもよい。
【００７６】
教示されていない条件に対する一般化のための機械学習システムのために、特徴自体が、一般的である必要がある。例えば、圧力センサ１２６、１２８の両方の絶対的圧力の大きさの使用は、良好に一般化され得ない。患者１１０に対しておよび患者によって印加される圧力は、（ｉ）処方負圧および患者の呼吸系によって印加される力と、（ｉｉ）患者の胸膜腔内の死体積の空気量とに関連するであろう。
【００７７】
圧力センサ１２６、１２８から受信されたデータは、流体キャニスタ１０８の流体体積に依存し、より一般的特徴が、形成されることができる。以下は、ある実施形態において使用され得る、特徴の部分的リストである。
・ピーク遅延（例えば、マイクロ秒またはミリ秒単位の時間）−最遅圧力センサ波形（流体キャニスタ１０８の最も可能性の高い圧力）を使用して、そのピーク圧力から通気圧力ピークのピーク圧力に戻る時間を報告する。
・通気に関する周期最大−最小（圧力）
・通気圧力ピーク間の距離（ミリ秒等の時間の単位）
・流体キャニスタ１０８に関する周期最大−最小（圧力）
・流体キャニスタ１０８の圧力ピーク間の距離（ミリ秒等の時間の単位）
・通気圧力ピーク時の流体キャニスタ１０８と通気との間の圧力デルタ
・ＣＤＵ圧力ピーク時の流体キャニスタ１０８と通気との間の圧力デルタ
・通気圧力の傾き（ピークにつながる）
・流体キャニスタ１０８圧力の傾き（ピークにつながる）
・圧力曲線下面積
・弁タイミング（開放／閉鎖時間）
【００７８】
本発明の一実施形態は、遺伝子プログラミング（ＧＰ）を通した機械学習を特徴とする。仮定された解を最適化する、遺伝子アルゴリズムおよび機械学習アプローチと異なり、ＧＰは、問題を解くための「プログラム」を生成する。学習プロセスの結果は、適用される問題を解く、「公式／ソースコード」である。
【００７９】
以下を含む、遺伝子プログラミングに対するいくつかの異なるアプローチがある。
−命令の関数木に基づく、ＧＰ−標準的遺伝子プログラミング
−ＬＧＰ−線形遺伝子プログラミング、線形に処理される固定命令長（高速）
−ＣＧＰ−デカルト遺伝子プログラミング、ランダムに配線されたランダム命令点の固定グリッド
ソフトウェア１０４内に容易に埋め込まれるＣコードを生成することができる、市販のＬＧＰアプリケーションの実施例は、Ｄｉｓｃｉｐｕｌｕｓ  ５  Ｇｅｎｅｔｉｃ  Ｐｒｏｇｒａｍｍｉｇ  Ｐｒｅｄｉｃｔｉｖｅ  Ｍｏｄｅｌｉｎｇ（ＲＭＬ  Ｔｅｃｈｎｏｌｏｇｉｅｓ，  Ｉｎｃ．（Ｌｉｔｔｌｅｔｏｎ，  Ｃｏｌｏｒａｄｏ）ｈｔｔｐ：／／ｗｗｗ．ｒｍｌｔｅｃｈ．ｃｏｍ／）である。
【００８０】
ある実施形態では、通気圧力センサ１２６は、省略されてもよい。これは、例えば、外部から印加される刺激の特性が既知である場合、生じてもよい。体積および他の特徴は、依然として、既知の患者ソースに基づいて、流体キャニスタ１０８内のＣＤＵ圧力センサ１２８を使用して判定されてもよい。例えば、患者に心肺機械または人工呼吸器が取り付けられている場合、そのような外部刺激は、既知である。
【００８１】
変更が、本明細書の範囲から逸脱することなく、前述の方法およびシステムに行われ得る。したがって、前述の説明に含有される、または付随の図面に図示される事項は、例証として解釈されるべきであって、限定的意味ではないことに留意されたい。以下の請求項は、本明細書に説明される全ての一般的および具体的特徴、ならびに文言上、本明細書の範囲内にあると言え得る、本方法およびシステムの範囲の全ての記載を網羅することが意図される。
【図１】
2018140200.tif 000037
【図２】
2018140200.tif 000038
【図３】
2018140200.tif 000039
【図４】
2018140200.tif 000040
【図５】
2018140200.tif 000041
【図６】
2018140200.tif 000042
【外国語明細書】
2018140200000001.pdf