WO2014208343A1 - Method and apparatus for assessment of depth of anesthesia, and use - Google Patents

Abstract

The present invention provides a method for assessment of the depth of anesthesia, the method including a step in which the major-axis-direction sharpness of a Poincare plot distribution is quantified and a step in which randomness contained in the electroencephalogram is evaluated on the basis of the sharpness to assess the depth of anesthesia.

Description

麻酔深度測定法、麻酔深度測定装置及び使用Anesthesia depth measurement method, anesthesia depth measurement device and use

　本発明は、麻酔深度測定法、麻酔深度測定装置及び麻酔深度判定のための使用に関する。 The present invention relates to an anesthesia depth measurement method, an anesthesia depth measurement device, and use for anesthesia depth determination.

　麻酔法の進歩に伴い術中覚醒が社会的問題となり、BISモニター、Cerebral State Index、AEPモニターなどの脳波由来の麻酔鎮静モニターがガイドライン等で推奨され世界で使用されている。また、的確な麻酔深度を維持することは、術後の認知機能低下を防ぐ上でも有用であり、麻酔深度モニタリングはこのような術後の高次認知機能障害を予防する立場からも推奨される（非特許文献１，２）。 に With the progress of anesthesia, intraoperative awakening has become a social problem, and brain wave-derived anesthesia and sedation monitors such as BIS monitors, Cerebral State Index, and AEP monitors are recommended by guidelines and are used worldwide. Maintaining an appropriate depth of anesthesia is also useful for preventing postoperative cognitive decline, and anesthesia depth monitoring is recommended from the standpoint of preventing such postoperative higher cognitive impairment. (Non-Patent Documents 1 and 2).

　一方で、これらのモニターは術中覚醒を防ぐには十分でなく、また有用性に懐疑的な報告が多くある(非特許文献３～８）。いずれのモニターも麻酔深度推定に長時間の脳波信号を要し解析に時間がかかるため遅延が生じ、リアルタイムに麻酔深度を感知できないことが主要な原因の一つである(非特許文献９～１０)。例えばBISモニターは最短でも60秒以上の遅延が生じる。 On the other hand, these monitors are not enough to prevent intraoperative awakening, and there are many reports skeptical about their usefulness (Non-Patent Documents 3 to 8). Each monitor requires a long time electroencephalogram signal to estimate the depth of anesthesia and takes a long time to analyze, causing a delay, and one of the main causes is that the depth of anesthesia cannot be sensed in real time (Non-Patent Documents 9 to 10). ). For example, a BIS monitor has a delay of 60 seconds or more at the shortest.

　臨床では、麻酔薬濃度が一定に維持されているにもかかわらず、術操作に伴う侵害刺激が加わることで、突如、麻酔深度が相対的に浅くなる状況に多々遭遇する。このような術侵襲に伴う急激な麻酔深度の変化には、現存の麻酔モニターでは直ちに追随できず、これが術中覚醒の主要な原因の一つである。一方で術中覚醒がおきても30秒以上持続しなければ記憶に残らない（非特許文献４）ことが知られる。従って麻酔深度変化を遅延なく感知できれば、このような術中覚醒の問題の多くは回避できるはずである。今後の麻酔深度モニターは、深度変化に俊敏に遅延なく反応できる性能が求められる。 In clinical practice, despite the fact that the anesthetic concentration is maintained at a constant level, there are many situations where the depth of anesthesia suddenly becomes relatively shallow due to the addition of noxious stimuli associated with the operation. Such a rapid change in anesthesia depth due to surgical invasion cannot be immediately followed by an existing anesthesia monitor, and this is one of the main causes of intraoperative arousal. On the other hand, even if awakening occurs during surgery, it is known that it will not remain in memory unless it lasts for 30 seconds or more (Non-patent Document 4). Therefore, if anesthesia depth change can be sensed without delay, many of these problems of intraoperative awakening should be avoided. Future anesthesia depth monitors will need to be able to react quickly to depth changes without delay.

　Poincare plot（ポアンカレプロット）は、低次元の2次元状態空間での簡易な非線形解析法であり、従来心電図のＲＲ間隔を用いた心拍変動解析、自律神経機能解析での研究には用いられている(特許文献１) Poincare plot (Poincare plot) is a simple non-linear analysis method in a low-dimensional two-dimensional state space, and has been used for research in heart rate variability analysis and autonomic function analysis using the RR interval of ECG. (Patent Document 1)

特許第4416249 号Patent No. 4416249

本発明は、従来のバイスペクトラル、スペクトラル解析や、聴覚誘発電位の解析とは異なる新しい方法で、非定常な状態や、短時間脳波信号からでも俊敏に算定できる麻酔深度測定法及び麻酔深度測定装置を開発することを目的とする。 The present invention is a new method different from conventional bispectral, spectral analysis, and auditory evoked potential analysis, and an anesthesia depth measurement method and anesthesia depth measurement device capable of agile calculation even from a non-stationary state or a short-time electroencephalogram signal. The purpose is to develop.

本発明は、以下の麻酔深度測定法、麻酔深度測定装置及び麻酔深度判定のための使用を提供するものである。
項１．　ポワンカレプロット分布の長軸方向への先鋭度を定量化する工程、前記先鋭度に基づき脳波に含まれる乱雑性を評価して麻酔深度を判定する工程を含む、麻酔深度測定法。
項２．　以下の工程を含む、麻酔深度測定法
（１）脳波計により全身麻酔患者の脳波を測定する工程
（２）前記脳波をフィルタリング処理とアナログ／デジタル変換して、脳波時系列信号を出力する工程（但し、アナログ／デジタル変換及びフィルタリング処理の順序は問わない）
（３）前記脳波時系列信号と前記脳波時系列信号の遅延時系列信号を、ある時間での脳波電圧Axと前記ある時間から所定時間の遅延後の脳波電圧Ayのペア（Ax, Ay）でｘｙ平面上にプロットする工程、
（４）前記ｘｙ平面上にプロットされた前記脳波電圧Ayの前記ペア（Ax, Ay）を２つの直交する軸からの変位で定量化し、前記２つの直交する軸に対する前記ペア（Ax, Ay）のそれぞれの変位L1、変位L2に基づいて、前記ペア（Ax, Ay）の前記２つの直交する軸に対する距離に関係づけられた量であるちらばり度V1、ちらばり度V2をそれぞれ求めるちらばり度算出工程、
（５）前記ちらばり度V1と前記ちらばり度V2との比に基づき、麻酔深度を判定する工程。
項３．　前記直交する軸は、前記ｘｙ平面上のそれぞれ長軸（式ｙ＝ｘで表される軸）と短軸（式ｙ＝ｘで表される軸に直交する軸）とである、項２に記載の麻酔深度測定法。
項４．　前記変位L1を前記短軸方向の変位とし、前記変位L2を前記長軸方向の変位とし、前記ちらばり度算出工程は、前記変位L1と前記変位L2を以下の式に従い算出するL1・L2算出工程を含み、L1とL2の比（L1/L2もしくはL2/L1）に基づいて麻酔深度を判定する、項２又は３に記載の麻酔深度測定法 The present invention provides the following anesthesia depth measurement methods, anesthesia depth measurement devices, and uses for anesthetic depth determination.
Item 1. An anesthesia depth measurement method including the step of quantifying the sharpness in the major axis direction of the Poincare plot distribution, and the step of evaluating the randomness included in the electroencephalogram based on the sharpness and determining the depth of anesthesia.
Item 2. Anesthesia depth measurement method including the following steps: (1) A step of measuring the electroencephalogram of a general anesthesia patient by an electroencephalograph (2) A step of outputting the electroencephalogram time-series signal by subjecting the electroencephalogram to filtering processing and analog / digital conversion ( (However, the order of analog / digital conversion and filtering processing does not matter)
(3) The brain wave time series signal and the delayed time series signal of the brain wave time series signal are expressed as a pair (Ax, Ay) of a brain wave voltage Ax at a certain time and a brain wave voltage Ay after a certain time from the certain time. plotting on the xy plane;
(4) The pair (Ax, Ay) of the electroencephalogram voltage Ay plotted on the xy plane is quantified by displacement from two orthogonal axes, and the pair (Ax, Ay) with respect to the two orthogonal axes Based on the respective displacements L1 and L2, the dispersion degree V1 and the dispersion degree V2, which are quantities related to the distance of the pair (Ax, Ay) with respect to the two orthogonal axes, are obtained. Degree calculation process,
(5) A step of determining the depth of anesthesia based on the ratio between the degree of dispersion V1 and the degree of dispersion V2.
Item 3. In the item 2, the orthogonal axes are a long axis (axis represented by the expression y = x) and a short axis (axis orthogonal to the axis represented by the expression y = x), respectively, on the xy plane. Described anesthesia depth measurement method.
Item 4. The displacement L1 is the displacement in the minor axis direction, the displacement L2 is the displacement in the major axis direction, and the dispersion degree calculating step calculates the displacement L1 and the displacement L2 according to the following formulas: L1 / L2 calculation Item 4. Anesthesia depth measurement method according to item 2 or 3, wherein the anesthesia depth is determined based on the ratio of L1 and L2 (L1 / L2 or L2 / L1).

項５．　前記変位L1を前記短軸方向の変位とし、前記変位L2を前記長軸方向の変位とし、前記ちらばり度算出工程は、前記変位L1と前記変位L2を以下の式に従い算出するL1・L2算出工程と Item 5. The displacement L1 is the displacement in the minor axis direction, the displacement L2 is the displacement in the major axis direction, and the dispersion degree calculating step calculates the displacement L1 and the displacement L2 according to the following formulas: L1 / L2 calculation Process and

前記変位L1の標準偏差SD1と前記変位L2の標準偏差SD2とを求めるSD1・SD2算出工程とを含み、
前記ちらばり度V1、前記ちらばり度V2は、それぞれ前記SD1または前記SD2であり、
SD1/SD2またはSD2/SD1に基づいて、麻酔深度を判定する、項２又は３に記載の麻酔深度測定法。
項６．　前記SD1/SD2または前記SD2/SD1と他の麻酔深度の解析手段を組み合わせて前記麻酔深度を判定する、項５に記載の麻酔深度測定法。
項７．　前記他の麻酔深度の解析手段が、前記脳波の周波数解析、バイスペクトラル解析、スペクトラル解析、エントロピー解析、聴覚誘発電位による解析手段、又はこれらの解析手段の２種以上の組み合わせである、項６に記載の麻酔深度測定法。
項８．　全身麻酔患者の脳波を測定する脳波計と、
前記脳波のアナログ／デジタル(A/D)変換手段と、
A/D変換後、またはA/D変換前の脳波をフィルタリング処理するフィルタリング手段と、
前記脳波をフィルタリング処理して得られる脳波時系列信号と前記脳波時系列信号の遅延時系列信号を、ある時間での脳波電圧Axと前記ある時間から所定時間の遅延後の脳波電圧Ayのペア（Ax, Ay）でｘｙ平面上にプロットし、
前記ｘｙ平面上にプロットされた前記脳波電圧Ayの前記ペア（Ax, Ay）を２つの直交する軸からの変位で定量化し、
前記２つの直交する軸に対する前記ペア（Ax, Ay）のそれぞれの変位V1、変位V2に基づいて、前記ペア（Ax, Ay）の前記２つの直交する軸に対する距離に関係づけられた量であるちらばり度V1、ちらばり度V2をそれぞれ求めてちらばり度を算出する演算手段と、
前記ちらばり度V1と前記ちらばり度V2との比に基づき、麻酔深度を判定する麻酔深度判定手段と、
前記判定結果を表示する表示手段と、
を備えた麻酔深度測定装置。
項９．　前記直交する軸は、前記ｘｙ平面上のそれぞれ長軸（式ｙ＝ｘで表される軸）と短軸（式ｙ＝ｘで表される軸に直交する軸）とであって、
前記変位L1を前記短軸方向の変位とし、前記変位L2を前記長軸方向の変位とし、
前記変位L1と前記変位L2を以下の式に従い算出し、 SD1 and SD2 calculation step for obtaining the standard deviation SD1 of the displacement L1 and the standard deviation SD2 of the displacement L2,
The degree of dispersion V1 and the degree of dispersion V2 are the SD1 or the SD2, respectively.
Item 4. Anesthetic depth measurement method according to Item 2 or 3, wherein the depth of anesthesia is determined based on SD1 / SD2 or SD2 / SD1.
Item 6. Item 6. The method for measuring the depth of anesthesia according to item 5, wherein the depth of anesthesia is determined by combining the SD1 / SD2 or SD2 / SD1 and another means for analyzing anesthesia depth.
Item 7. Item 6. The other anesthetic depth analysis means is the electroencephalogram frequency analysis, bispectral analysis, spectral analysis, entropy analysis, auditory evoked potential analysis means, or a combination of two or more of these analysis means. Described anesthesia depth measurement method.
Item 8. An electroencephalograph to measure the electroencephalogram of a general anesthesia patient;
An analog / digital (A / D) conversion means for the electroencephalogram;
Filtering means for filtering brain waves after A / D conversion or before A / D conversion;
An electroencephalogram time-series signal obtained by filtering the electroencephalogram and a delayed time-series signal of the electroencephalogram time-series signal are represented by a pair of electroencephalogram voltage Ax at a certain time and electroencephalogram voltage Ay after a certain time delay from the certain time ( Plot on the xy plane with (Ax, Ay)
Quantifying the pair (Ax, Ay) of the electroencephalogram voltage Ay plotted on the xy plane by displacement from two orthogonal axes;
Based on the respective displacements V1 and V2 of the pair (Ax, Ay) relative to the two orthogonal axes, the amount is related to the distance of the pair (Ax, Ay) relative to the two orthogonal axes Calculation means for calculating the degree of dispersion by obtaining the degree of dispersion V1 and the degree of dispersion V2, respectively;
Based on the ratio of the degree of dispersion V1 and the degree of dispersion V2, the anesthetic depth determination means for determining the depth of anesthesia,
Display means for displaying the determination result;
An anesthesia depth measuring device with
Item 9. The orthogonal axes are a long axis (axis represented by the expression y = x) and a short axis (axis orthogonal to the axis represented by the expression y = x), respectively, on the xy plane,
The displacement L1 is the displacement in the minor axis direction, the displacement L2 is the displacement in the major axis direction,
Calculate the displacement L1 and the displacement L2 according to the following formula,

前記変位L1の標準偏差SD1と前記変位L2の標準偏差SD2とを求め、
前記ちらばり度V1、前記ちらばり度V2は、それぞれ前記SD1または前記SD2であり、
SD1/SD2またはSD2/SD1に基づいて、麻酔深度を判定する麻酔深度判定手段と、
を備えた項８に記載の麻酔深度測定装置。
項１０．　前記麻酔深度判定手段が、前記SD1/SD2または前記SD2/SD1と他の麻酔深度の解析手段を組み合わせて前記麻酔深度を判定する手段である、項９に記載の麻酔深度測定装置。
項１１．　前記他の麻酔深度の解析手段が、前記脳波の周波数解析、バイスペクトラル解析、スペクトラル解析、エントロピー解析、聴覚誘発電位による解析手段、又はこれらの解析手段の２種以上の組み合わせである、項１０に記載の麻酔深度測定装置。
項１２．　ちらばり度V1と前記ちらばり度V2との比の麻酔深度判定のための使用であって、脳波計により全身麻酔患者の脳波を測定し、
前記脳波をアナログ／デジタル変換後、またはアナログ／デジタル変換前に、フィルタリング処理して脳波時系列信号を出力し、
前記脳波時系列信号と前記脳波時系列信号の遅延時系列信号を、ある時間での脳波電圧Axと前記ある時間から所定時間の遅延後の脳波電圧Ayのペア（Ax, Ay）でｘｙ平面上にプロットし、
前記ｘｙ平面上にプロットされた前記脳波電圧Ayの前記ペア（Ax, Ay）を２つの直交する軸からの変位で定量化し、
前記２つの直交する軸に対する前記ペア（Ax, Ay）のそれぞれの変位L1、変位L2に基づいて、前記ペア（Ax, Ay）の前記２つの直交する軸に対する距離に関係づけられた量であるちらばり度V1、ちらばり度V2をそれぞれ求めてちらばり度V1と前記ちらばり度V2との比を算出することを含む、使用。
項１３．　前記直交する軸は、前記ｘｙ平面上のそれぞれ長軸（式ｙ＝ｘで表される軸）と短軸（式ｙ＝ｘで表される軸に直交する軸とであって、
前記変位L1を前記短軸方向の変位とし、前記変位L2を前記長軸方向の変位とし、
前記変位L1と前記変位L2を以下の式に従い算出し、 Obtain the standard deviation SD1 of the displacement L1 and the standard deviation SD2 of the displacement L2,
The degree of dispersion V1 and the degree of dispersion V2 are the SD1 or the SD2, respectively.
Based on SD1 / SD2 or SD2 / SD1, anesthetic depth determination means for determining the depth of anesthesia,
The anesthesia depth measuring apparatus of claim | item 8 provided with.
Item 10. Item 10. The anesthetic depth measurement device according to Item 9, wherein the depth of anesthesia determination unit is a unit that determines the depth of anesthesia by combining the SD1 / SD2 or SD2 / SD1 and another anesthesia depth analysis unit.
Item 11. Item 10. The other anesthetic depth analysis means is the electroencephalogram frequency analysis, bispectral analysis, spectral analysis, entropy analysis, auditory evoked potential analysis means, or a combination of two or more of these analysis means. The anesthesia depth measuring device as described.
Item 12. The use of the ratio of the degree of dispersion V1 and the degree of dispersion V2 for determining the depth of anesthesia, measuring the electroencephalogram of a general anesthesia patient with an electroencephalograph,
After analog / digital conversion of the electroencephalogram or before analog / digital conversion, filtering process is performed to output an electroencephalogram time-series signal,
The brain wave time series signal and the delayed time series signal of the brain wave time series signal are represented on the xy plane by a pair (Ax, Ay) of a brain wave voltage Ax at a certain time and a brain wave voltage Ay after a certain time from the certain time. Plot to
Quantifying the pair (Ax, Ay) of the electroencephalogram voltage Ay plotted on the xy plane by displacement from two orthogonal axes;
An amount related to the distance of the pair (Ax, Ay) to the two orthogonal axes based on the respective displacements L1, L2 of the pair (Ax, Ay) relative to the two orthogonal axes Use including calculating the degree of dispersion V1 and the degree of dispersion V2, respectively, and calculating the ratio of the degree of dispersion V1 and the degree of dispersion V2.
Item 13. The orthogonal axes are respectively a long axis (axis represented by the expression y = x) and a short axis (axis orthogonal to the axis represented by the expression y = x) on the xy plane,
The displacement L1 is the displacement in the minor axis direction, the displacement L2 is the displacement in the major axis direction,
Calculate the displacement L1 and the displacement L2 according to the following formula,

前記変位L1の標準偏差SD1と前記変位L2の標準偏差SD2とを求め、
前記ちらばり度V1、前記ちらばり度V2は、それぞれ前記SD1または前記SD2である、
SD1/SD2またはSD2/SD1の麻酔深度判定のための項１２に記載の使用。
項１４．　前記SD1/SD2または前記SD2/SD1を、脳波の周波数解析、バイスペクトラル解析、スペクトラル解析、エントロピー解析、聴覚誘発電位による解析手段からなる群から選ばれる少なくとも1種と組み合わせて麻酔深度の判定を行う、項１３に記載の使用。
項１５．　前記SD1/SD2または前記SD2/SD1を麻酔深度の測定のためのアルゴリズムと組み合わせて麻酔深度の判定を行う、項１３に記載の使用。 Obtain the standard deviation SD1 of the displacement L1 and the standard deviation SD2 of the displacement L2,
The dispersion degree V1 and the dispersion degree V2 are the SD1 or the SD2, respectively.
Item 15. The use according to Item 12, for determining the depth of anesthesia for SD1 / SD2 or SD2 / SD1.
Item 14. The SD1 / SD2 or SD2 / SD1 is combined with at least one selected from the group consisting of electroencephalogram frequency analysis, bispectral analysis, spectral analysis, entropy analysis, and auditory evoked potential analysis means to determine the depth of anesthesia Item 14. Use according to Item 13.
Item 15. Item 14. The use according to Item 13, wherein the SD1 / SD2 or SD2 / SD1 is combined with an algorithm for measuring the depth of anesthesia to determine the depth of anesthesia.

脳波信号（4 sec. 及び10 sec.）のポワンカレプロットの描出パターンをちらばり度V1とちらばり度V2の比で定量化したところ、いずれも吸入麻酔薬濃度に応じた変化を示すことがわかった。更に、ちらばり度V1とちらばり度V2の比は、脳波のパワースペクトル解析におけるδ帯域低周波成分とβ帯域高周波成分の比率(δ/β)とも高い相関をもつことも明らかになった。麻酔深度が深まると、脳波上δ帯域低周波成分が増大しβ帯域高周波成分が減少することが知られている。これらより、前頭誘導脳波のポワンカレプロット定量化指標（V1とV2の比）が麻酔深度測定に有用であることがわかる。本法は、秒単位の短時間の時系列信号から容易に実現できる方法であり、リアルタイム麻酔深度モニタリングとしての臨床応用が期待できる。 When the pattern of the Poincare plot of the electroencephalogram signal (4 sec. And 10 sec.) Was quantified by the ratio of the dispersity V1 and dispersity V2, it was found that both showed changes depending on the concentration of the inhaled anesthetic. It was. Furthermore, the ratio between the degree of dispersion V1 and the degree of dispersion V2 was also found to have a high correlation with the ratio (δ / β) of the δ band low frequency component and the β band high frequency component in the power spectrum analysis of the electroencephalogram. It is known that as the anesthetic depth increases, the δ band low frequency component on the electroencephalogram increases and the β band high frequency component decreases. From these, it can be seen that the Poincare plot quantification index (ratio of V1 and V2) of the frontal lead electroencephalogram is useful for measuring the depth of anesthesia. This method can be easily realized from a time-series signal of a short time in seconds, and clinical application as real-time anesthesia depth monitoring can be expected.

　本発明は、新しい非線形法を脳波深度解析に用いることで、脳活動の急な変化にも対応でき、また解析に必要な脳波信号も短く秒単位である。従って、迅速に反応する。 The present invention can cope with a sudden change in brain activity by using a new non-linear method for the electroencephalogram depth analysis, and the electroencephalogram signal necessary for the analysis is short and in seconds. Therefore, it reacts quickly.

ランダム信号とランダム信号に線形関係成分を含有した信号のポアンカレプロットのシミュレーション。Simulation of a Poincare plot of a random signal and a signal containing a linear relationship component in the random signal. 脳波のポワンカレプロットとその定量化Evan's Poincare plot and its quantification 脳波収集の概要図Overview of EEG collection ポアンカレプロット定量化のための説明図Illustration for Poincare plot quantification 実施例１のポアンカレプロットを示す。The Poincare plot of Example 1 is shown. 実施例２のSD1/SD2の結果を示す。The result of SD1 / SD2 of Example 2 is shown. 脳波ポワンカレ定量化指標（SD1/SD2）とパワースペクトルの関係を示す。The relationship between EEG Poincare quantification index (SD1 / SD2) and power spectrum is shown. 遅延時間がポワンカレプロット（SD1/SD2指標）に及ぼす影響を示す。The influence of delay time on the Poincare plot (SD1 / SD2 index) is shown. （A）深麻酔時（セボフルラン吸入麻酔濃度3％）の脳波解析例、（B）浅麻酔時（セボフルラン吸入麻酔濃度3％から0.5%に変更約6分後）の脳波解析例(A) EEG analysis example during deep anesthesia (sevoflurane inhalation anesthesia concentration 3%), (B) EEG analysis example under shallow anesthesia (sevoflurane inhalation anesthesia concentration changed from 3% to 0.5% after about 6 minutes) ポワンカレプロット定量値(SD1/SD2)と麻酔深度指標（BIS, SEF95）との比較Comparison between Poincare plot quantitative value (SD1 / SD2) and anesthesia depth index (BIS, SEF95) SD1/SD2と30秒後のSEF95、BISとの関係（Point=158）Relationship between SD1 / SD2 and SEF95 and BIS after 30 seconds (Point = 158) 吸入セボフルラン濃度を0.5%から3%に上げて、再度0.5%に戻した時のポワンカレプロット定量値(SD1/SD2)と麻酔深度指標（BIS, SEF95）との比較Comparison of Poincare plot quantitative value (SD1 / SD2) and depth of anesthesia index (BIS, SEF95) when inhalation sevoflurane concentration is increased from 0.5% to 3% and back to 0.5% 麻酔をプロポフォール静脈麻酔薬で維持した時の麻酔覚醒前60分間の推移Transition of 60 minutes before anesthesia awakening when anesthesia was maintained with propofol intravenous anesthetic 本発明の麻酔深度測定法のフローチャートFlow chart of anesthesia depth measurement method of the present invention 10秒毎に同時採取したpoincare plotのSD1/SD2値とSEF95値 (51 cases, 32351 points)　青丸:プロポフォール(13 cases, 9059 points), 緑:デスフルラン(22 cases, 14805 points), オレンジ:セボフルラン(16cases, 8487 oiunts), 黒線: 全データへの回帰直線、赤線: 混合モデルの回帰直線SD / SD2 and SEF95 values of poincare samples collected every 10 seconds (51 cases, 32351 points) Blue circle: Propofol (13 cases, 9059 points), Green: Desflurane (22 cases, 14805 points), Orange: Sevoflurane (16cases, 8487 oiunts), black line: 回 帰 regression line to all data, red line: regression line of mixed model 全データ(51 cases, 32,351 points)のSD1/SD2とSEF95の関係図とそのBland-Altman plotRelationship diagram of SD1 / SD2 and SEF95 of all data (51 cases, 32,351 points) and its Bland-Altman plot propofol(A), desflurane(B) sevoflurane (C)それぞれの麻酔におけるSD1/SD2とSEF95の回帰分析とBland-Altman plot.Proofol (A), desflurane (B) sevoflurane (C) SD1 / SD2 and SEF95 regression analysis and Bland-Altman plot. 35歳女性のセボフルラン麻酔中の経過例A 35-year-old woman undergoing sevoflurane anesthesia 56歳男性のセボフルラン麻酔中の経過例A 56-year-old man undergoing sevoflurane anesthesia

　現在、さまざまな手法を取り入れた麻酔深度計が開発され臨床使用されているが、術中覚醒を予防し適度な麻酔深度を維持する上で充分ではない。脳は複雑系、非線形システムであるので、その挙動をモニターする麻酔深度測定には、ダイナミックな信号を解析する非線形的多面的アプローチが有用である。ポワンカレプロットは、前提やモデルを必要とせずに、現時間信号と遅延時間信号とを対でプロットすることで、ダイナミックな信号を過去の履歴の枠組みの中で分析する非線形手法である。 Currently, anesthesia depth meters incorporating various techniques have been developed and are in clinical use, but are not sufficient to prevent intraoperative awakening and maintain an appropriate depth of anesthesia. Since the brain is a complex system and a nonlinear system, a nonlinear multi-faceted approach to analyze a dynamic signal is useful for measuring the depth of anesthesia to monitor its behavior. The Poincare plot is a non-linear method that analyzes a dynamic signal in the framework of a past history by plotting a current time signal and a delay time signal in pairs without requiring a premise or a model.

　図１に示すように、ホワイトノイズなどのランダムな複雑な信号のポワンカレプロットは一般に卵丸状の分布を示すが、信号の乱雑度が低下するにつれて、細長い分布になる。脳波は覚醒時には無秩序でランダムであるが、麻酔を深めるに従って同期化して規則性を増す。このような麻酔深度を増すにつれて観察される脳波の乱雑度の低下は、脳波ポワンカレプロットパターンの先鋭度の増加として現れる。 As shown in FIG. 1, the Poincare plot of a random complex signal such as white noise generally shows an egg-shaped distribution, but becomes an elongated distribution as the signal randomness decreases. EEGs are disordered and random when awake, but synchronize and become more regular as anesthesia deepens. The decrease in the degree of electroencephalogram observed as the depth of anesthesia increases as the sharpness of the electroencephalogram Poincare plot pattern increases.

　本発明の好ましい実施形態においては、ポワンカレプロット分布の長軸（ｙ＝ｘ軸）方向への先鋭度を定量化することで脳波に含まれる乱雑性を評価し、麻酔深度を測定する。好ましい実施形態において、定量化は、ポワンカレプロット分布の長軸（ｙ＝ｘ軸）方向のちらばり度をｖ２、短軸方向のちらばり度をｖ１、並びにその比（Ｖ１／Ｖ２）を用いる。ちらばり度は、長軸、短軸方向の標準偏差（SD2,SD1）、並びにその比（SD1/SD2）を用いるのが典型的であるが、ポワンカレプロット分布の中心から長軸、短軸方向への変位距離（L2,　L1）、並びにその比（L1/L2）を用いることもできる。また、パターン分析の軸は、長軸、短軸方向に限らず、長軸、短軸方向のベクトル成分を含有するものなら定量化が可能である。例えば、X軸方向、Y軸方向のベクトルを用いて、(短軸, X軸方向), (長軸, X軸方向), (短軸, Y軸方向), (長軸, Y軸方向), ( X軸方向, Y軸方向) 等の変位の組み合わせ、並びにそれらの標準偏差の組み合わせ等を用いることもできる。 In a preferred embodiment of the present invention, the degree of anesthesia included in the electroencephalogram is evaluated by quantifying the sharpness of the Poincare plot distribution in the long axis (y = x axis) direction, and the anesthetic depth is measured. In a preferred embodiment, quantification uses v2 for the degree of dispersion in the major axis (y = x axis) direction of the Poincare plot distribution, v1 for the degree of dispersion in the minor axis direction, and the ratio (V1 / V2). For the degree of dispersion, the standard deviation (SD2, SD1) and ratio (SD1 / SD2) in the major axis and minor axis direction are typically used, but from the center of the Poincare plot distribution to the major axis and minor axis direction. Displacement distance (L2, L1) and its ratio (L1 / L2) can also be used. In addition, the axis of pattern analysis is not limited to the major axis and minor axis directions, and can be quantified as long as it contains vector components in the major axis and minor axis directions. For example, using vectors in the X-axis direction and Y-axis direction, (short axis, X-axis direction), (long-axis, X-axis direction), (short-axis, Y-axis direction), (long-axis, Y-axis direction) , (X-axis direction, Y-axis direction), etc., as well as combinations of their standard deviations.

　このようにして、脳波ポワンカレプロット分布の長軸（ｙ＝ｘ軸）方向への先鋭度をちらばり度として定量化することで、脳波信号の乱雑度の低下と線形度を計測し、麻酔深度を測定する。 In this way, the degree of sharpness in the major axis (y = x axis) direction of the EEG Poincare plot distribution is quantified as the degree of dispersion, thereby measuring the degree of randomness and linearity of the electroencephalogram signal and determining the depth of anesthesia. Measure.

　本発明において、脳波計は特に限定されず、通常使用されている脳波計を広く使用することができる。 In the present invention, the electroencephalograph is not particularly limited, and a commonly used electroencephalograph can be widely used.

　アナログ／デジタル変換(A/D変換)手段は特に限定されず、通常使用されているA/D変換を広く使用することができる。 Analog / digital conversion (A / D conversion) means is not particularly limited, and A / D conversion that is normally used can be widely used.

　脳波のフィルタリングは、低周波と高周波の両方を遮断するのが好ましい。 EEG filtering is preferably to block both low and high frequencies.

　脳波のフィルタリング手段としては、高周波数通過フィルター、低周波数通過フィルターを使用することができる。 As a means for filtering the electroencephalogram, a high frequency pass filter or a low frequency pass filter can be used.

　A/D変換及びフィルタリング後の一定の周波数範囲の脳波電圧の信号は、パーソナルコンピュータ等の演算手段に出力される。演算手段は、2次元状態空間における脳波時系列信号と前記脳波時系列信号の遅延時系列信号を、ある時間での脳波電圧Axと前記ある時間から所定時間の遅延後の脳波電圧Ayのペア（Ax, Ay）でｘｙ平面上にプロットし（ポアンカレプロット）、前記ｘｙ平面上にプロットされた前記脳波電圧Ayの前記ペア（Ax, Ay）を２つの直交する軸からの変位で定量化し、前記２つの直交する軸に対する前記ペア（Ax, Ay）のそれぞれの変位V1、変位V2に基づいて、前記ペア（Ax, Ay）の前記２つの直交する軸に対する距離に関係づけられた量であるちらばり度V1、ちらばり度V2をそれぞれ求めてちらばり度を算出する。ちらばり度V1、V2は、例えば以下のSD1、SD2として求めることができる。すなわち、前記ｘｙ平面上にプロットされた前記脳波電圧Ayの前記ペア（Ax, Ay）の２つの直交する軸は、前記ｘｙ平面上のそれぞれ長軸（式ｙ＝ｘで表される軸）と短軸（式ｙ＝ｘで表される軸に直交する軸）とであって、前記変位L1を前記短軸方向の変位とし、前記変位L2を前記長軸方向の変位とし、前記変位L1と前記変位L2を以下の式に従い算出し、 The EEG voltage signal in a certain frequency range after A / D conversion and filtering is output to a computing means such as a personal computer. The computing means converts a brain wave time series signal in a two-dimensional state space and a delayed time series signal of the brain wave time series signal into a pair of an electroencephalogram voltage Ax at a certain time and an electroencephalogram voltage Ay after a certain time from the certain time ( Ax, Ay) is plotted on the xy plane (Poincare plot), and the pair (Ax, Ay) of the electroencephalogram voltage Ay plotted on the xy plane is quantified by displacement from two orthogonal axes, Based on the respective displacements V1 and V2 of the pair (Ax, Ay) with respect to two orthogonal axes, an amount related to the distance of the pair (Ax, Ay) with respect to the two orthogonal axes. The spread degree is calculated by obtaining the spread degree V1 and the spread degree V2, respectively. The scatter degrees V1 and V2 can be obtained, for example, as SD1 and SD2 below. That is, two orthogonal axes of the pair (Ax, Ay) of the electroencephalogram voltage Ay plotted on the xy plane are respectively a long axis (axis expressed by the equation y = x) on the xy plane. A short axis (axis perpendicular to the axis represented by the expression y = x), wherein the displacement L1 is a displacement in the short axis direction, the displacement L2 is a displacement in the long axis direction, and the displacement L1 The displacement L2 is calculated according to the following formula:

前記変位L1の標準偏差SD1と前記変位L2の標準偏差SD2とを求め、前記ちらばり度V1、前記ちらばり度V2は、それぞれ前記SD1または前記SD2である。 A standard deviation SD1 of the displacement L1 and a standard deviation SD2 of the displacement L2 are obtained, and the scatter degree V1 and the scatter degree V2 are the SD1 or the SD2, respectively.

　さらに前記演算手段は、L1の標準偏差SD1とL2の標準偏差SD2、並びに、SD1/SD2比またはSD2/SD1比を算出する。SD1とSD2の比は、SD1/SD2とSD2/SD1のいずれを使用してもよいが、SD1/SD2を用いるとその対数がpower spectrumのβ／δ比と高い相関関係をもつので、解釈が便利である。演算手段としては、ＣＰＵが挙げられるが、これに限定されず任意の演算手段を用いることができる。 Further, the calculation means calculates the standard deviation SD1 of L1 and the standard deviation SD2 of L2, and the SD1 / SD2 ratio or SD2 / SD1 ratio. The ratio of SD1 and SD2 can be either SD1 / SD2 or SD2 / SD1, but when SD1 / SD2 is used, the logarithm has a high correlation with the power / spectrum β / δ ratio. Convenient. An example of the calculation means is a CPU, but is not limited to this, and any calculation means can be used.

　V1とV2の比、例えばSD1とSD2の比の信号は麻酔深度判定手段に送られて、麻酔の深度を判定する。麻酔の深度は覚醒(例えば０)から最も深い麻酔(例えば１００)までを連続的な数値で表してもよく、覚醒（１）、浅い麻酔２(覚醒の可能性有り)、適切な麻酔３、深い麻酔４(麻酔による副作用もしくは事故の可能性有り)などのように段階的に判定してもよい。判定結果の信号は表示手段に出力されて適切に麻酔深度を表示する。 ¡A signal of the ratio of V1 and V2, for example, the ratio of SD1 and SD2, is sent to the anesthesia depth determination means to determine the depth of anesthesia. The depth of anesthesia may be expressed as a continuous number from awakening (eg, 0) to deepest anesthesia (eg, 100): awakening (1), shallow anesthesia 2 (possible arousal), appropriate anesthesia 3, The judgment may be made step by step such as deep anesthesia 4 (there is a side effect or accident due to anesthesia). The determination result signal is output to the display means to appropriately display the depth of anesthesia.

　表示手段は麻酔の深度を視覚的（文字、色、光など）、聴覚的（音声、警告音など）に表示する。視覚的表示手段としてはディスプレイが挙げられ、聴覚的表示手段としてはスピーカー、音源チップ(例えばＣＰＵなどの音源)などが挙げられる。 The display means displays the depth of anesthesia visually (character, color, light, etc.) and auditory (sound, warning sound, etc.). The visual display means includes a display, and the auditory display means includes a speaker, a sound source chip (for example, a sound source such as a CPU), and the like.

　本発明の演算手段は、Poincare plot（ポアンカレプロット）の分布分散パターンを2方向で定量化する。ポアンカレプロットは、低次元の2次元状態空間での簡易な非線形解析法であり、従来心電図のＲＲ間隔を用いた心拍変動解析、自律神経機能解析での研究には用いられているが、これを脳波の時系列信号そのものに応用しその分布変動を定量化解析することで麻酔深度を測定する発想は現在にいたるまでない。 The calculation means of the present invention quantifies the distribution distribution pattern of Poincare plot (Poincare plot) in two directions. The Poincare plot is a simple non-linear analysis method in a low-dimensional two-dimensional state space, which has been used for research in heart rate variability analysis and autonomic function analysis using the RR interval of the electrocardiogram. There is no idea to measure the depth of anesthesia by applying it to the EEG time series signal itself and quantifying and analyzing the distribution variation.

　本発明で麻酔深度が判定可能な麻酔薬は特に限定されず、全ての麻酔薬が適用対象である。例えば現在使用されている麻酔薬としては、セボフルラン、デスフルラン、イソフルラン、プロポフォール、ミダゾラム、ジアゼパム、デクスメデトミジンなどが挙げられる。これらの他に、今後開発される麻酔薬は全て本発明の麻酔深度の測定対象となる。 The anesthetic that can determine the depth of anesthesia in the present invention is not particularly limited, and all anesthetics are applicable. For example, currently used anesthetics include sevoflurane, desflurane, isoflurane, propofol, midazolam, diazepam, dexmedetomidine and the like. In addition to these, all anesthetics that will be developed in the future are subject to anesthetic depth measurement according to the present invention.

　本発明は、V1とV2の比、好ましくはSD1とSD2の比により麻酔深度を判定する。V1とV2の比、特にSD1とSD2の比は麻酔薬にかかわらず麻酔深度と相関するが、麻酔薬ごとに判定基準を決定しておくこともできる。 In the present invention, the depth of anesthesia is determined by the ratio of V1 and V2, preferably the ratio of SD1 and SD2. The ratio between V1 and V2, especially the ratio between SD1 and SD2, correlates with the depth of anesthesia regardless of the anesthetic, but the criterion can be determined for each anesthetic.

　全身麻酔中の脳内電気的活動は、状態空間内で周期性・準周期性の軌道をとることが予想され、より低次元空間での解析が有用と本発明者は考えた。そこで、ポワンカレプロットを用いて２次元状態空間で脳波時系列信号を解析し、麻酔深度測定を行う。即ち、脳波時系列信号をその遅延時系列信号とで、２次元状態空間にポワンカレプロットとして描出して、その分布形態を定量化することで、麻酔深度を測定する。 The brain electrical activity during general anesthesia is expected to take periodic or quasi-periodic trajectories in the state space, and the present inventor thought that analysis in a lower dimensional space is useful. Therefore, an anesthesia depth measurement is performed by analyzing a brain wave time-series signal in a two-dimensional state space using a Poincare plot. That is, the anesthesia depth is measured by rendering the electroencephalogram time-series signal as a delayed time-series signal as a Poincare plot in a two-dimensional state space and quantifying the distribution form.

　図14は、本発明の麻酔深度測定法のフローチャートを示す図である。
以下、ちらばり度V1、ちらばり度V2としてSD1とSD2を使用して本発明をより詳細に説明する。 FIG. 14 is a view showing a flowchart of the anesthetic depth measurement method of the present invention.
Hereinafter, the present invention will be described in more detail using SD1 and SD2 as the dispersion degree V1 and the dispersion degree V2.

　ポワンカレプロット分布パターンは、プロット集団の短軸、長軸の２次元２方向のちらばりの比率を用いて定量的に解析する。即ち、短軸方向の変位L1の標準偏差(SD1:Standard Deviation)と長軸方向の変位L2の標準偏差SD2の比率を、SD1/SD2指標として定量化する(図２)。 The Poincare plot distribution pattern is quantitatively analyzed using the ratio of the two-dimensional and two-dimensional dispersion of the short axis and long axis of the plot group. That is, the ratio between the standard deviation (SD1: Standard Deviation) of the displacement L1 in the short axis direction and the standard deviation SD2 of the displacement L2 in the long axis direction is quantified as an SD1 / SD2 index (FIG. 2).

　以下に、前頭部の脳波を測定して麻酔深度を予測する本発明の１つの好ましい実施形態について説明するが、本発明は後頭部、頭頂部、左右の側頭部など、任意の脳波を測定して実施することもできる。 The following describes one preferred embodiment of the present invention for measuring the frontal brain waves to predict the depth of anesthesia. The present invention measures any brain waves such as the occipital region, the top of the head, and the left and right temporal regions. It can also be implemented.

　前頭部に脳波電極を装着して前頭誘導脳波を測定し、アナログ/デジタル変換、及びフィルタリング処置する。その脳波時系列信号をポアンカレプロットして、分布のSD1/SD2を算出することで定量化する。詳細を以下に記載する。 脳 Electroencephalogram electrodes are attached to the frontal region, frontal induced electroencephalograms are measured, and analog / digital conversion and filtering are performed. The EEG time-series signal is Poincare plotted and quantified by calculating SD1 / SD2 of the distribution. Details are described below.

１．脳波時系列信号の計測
前頭導出脳波を用いる場合の解析例を以下に述べる。脳波測定の電極配置の基準である国際10-20法に基づき、前頭誘導脳波を導出する。ここで、(Fp1：frontal montage）, Cz：reference）とし、更に neutral electrode をFpzに装着する。脳波計により導出した前頭誘導脳波は、256Hzサンプリングしてアナログデジタル変換（A/D変換）後、0.5Hz以下の低周波成分と30Hz以上の高周波成分を高周波数通過フィルター、及び低周波数通過フィルターをもちいて除去し、0.5Hz-30Hz成分を抽出する。この脳波時系列信号を、ポワンカレ定量化解析に用いる(図３)。 1. An example of analysis when using an electroencephalogram derived from the frontal measurement of EEG time-series signals is described below. Based on the international 10-20 method, which is the standard for electrode placement for electroencephalogram measurement, the frontal lead electroencephalogram is derived. Here, (Fp1: frontal montage), Cz: reference) and a neutral electrode is attached to Fpz. The frontal brain waves derived from the electroencephalograph are sampled at 256 Hz and converted from analog to digital (A / D conversion). Remove using 0.5Hz-30Hz component. This EEG time-series signal is used for the Poincare quantitative analysis (FIG. 3).

２．脳波出力信号のポワンカレ定量化解析
出力された脳波時系列信号を、その遅延時系列信号とで、2次元平面にポアンカレプロットする。即ち、
ある時系列信号：　　 ( a,  b,  c,  d,  e,  f,  g,  h,  i,  j……….. )　と、
その遅延時系列信号： (    a,　b,　c,  d,  e,  f,  g,  h,  i,  j…….)  との対で、
( a, b )、( b, c )、( c, d )、( d, e )、( e, f )、( f, g )….を順次、x-y平面上にプロットする。
遅延に用いる時間は、サンプリング間隔の1倍から数倍を選択する。2段階遅延させて解析する場合には、
( a,  b,  c,  d,  e,  f,  g,  h,  i,  j……….. )　の時系列信号に対して
(        a,  b,  c,  d,  e,  f,  g,  h,  i,  j……….. )
の遅延時系列信号が得られ、この場合は ( a, c )、( b, d )、( c, e )、( d, f )、( e, g )、( f, h )….のペアを順次x-y平面上にプロットすることになる。 2. The electroencephalogram time series signal output by the Poincare quantification analysis of the electroencephalogram output signal is Poincare plotted on a two-dimensional plane with the delayed time series signal. That is,
A time-series signal: (a, b, c, d, e, f, g, h, i, j ……… ..)
In combination with the delayed time series signal: (a, b, c, d, e, f, g, h, i, j …….)
(a, b), (b, c), (c, d), (d, e), (e, f), (f, g)... are sequentially plotted on the xy plane.
The time used for the delay is selected from 1 to several times the sampling interval. When analyzing with a two-stage delay,
(a, b, c, d, e, f, g, h, i, j ……….)
(a, b, c, d, e, f, g, h, i, j ……… ..)
In this case, (a, c), (b, d), (c, e), (d, f), (e, g), (f, h)…. Pairs will be plotted sequentially on the xy plane.

　このようにして、得られたポワンカレプロット分布パターンを定量的に解析するため、2方向でポワンカレ分布のちらばりを定量化する。定量化は、ｘ軸とｙ軸との双方に４５°の角度の直線（y=x）方向のちらばり（SD2）とそれに直角に交わる方向(SD1)の2方向で計測する。 In this way, in order to quantitatively analyze the obtained Poincare plot distribution pattern, the dispersion of the Poincare distribution is quantified in two directions. For quantification, measurement is performed in two directions, a scatter (SD2) in a straight line (y = x) direction at an angle of 45 ° on both the x-axis and the y-axis, and a direction (SD1) perpendicular thereto.

　ポアンカレプロットは時間遅れをサンプリング間隔として設定した場合、ある時間での脳波電圧がAxであるなら、そのサンプリングの1点先（遅延後）の脳波電圧Ay とのペアで（Ax, Ay）でｘｙ平面上にプロットされることになる。このようにして得た集団のx軸の平均（すなわち平均脳波電圧）をMxとすれば、ある点A(Ax, Ay)の集団内のちらばり具合は、長軸（式ｙ＝ｘで表される軸）、及び点M (Mx, Mx）を通り短軸（式ｙ＝ｘで表される軸に直交する軸）の2方向の変位として、距離L2及びL1で表される。 In the Poincare plot, when the time delay is set as the sampling interval, if the electroencephalogram voltage at a certain time is Ax, it is xy in the pair (Ax, Ay) with the electroencephalogram voltage Ay one point ahead of the sampling (after delay) It will be plotted on the plane. If the average of the x-axis of the group obtained in this way (that is, the average electroencephalogram voltage) is Mx, the dispersity in the group at a certain point A (Ax, Ay) is expressed by the long axis (formula y = x). ), And the displacement in two directions of the short axis (axis perpendicular to the axis represented by the equation y = x) passing through the point M (Mx, Mx), is represented by distances L2 and L1.

　L2、L1は、図４に表したように、正方形とその対角線の図形的性質から、以下の算出式で導かれる。 L2 and L1 are derived by the following calculation formula from the graphic properties of the square and its diagonal, as shown in FIG.

　ポワンカレプロット分布パターンは、このちらばりの短軸方向の変位L1の標準偏差(SD1:Standard Deviation)と長軸方向の変位L2の標準偏差(SD2)の比率として、SD1/SD2で定量化できる。 The Poincare plot distribution pattern can be quantified by SD1 / SD2 as the ratio of the standard deviation (SD1: Standard Deviation) of the displacement L1 in the short axis direction and the standard deviation (SD2) of the displacement L2 in the long axis direction.

　本発明は、SD1とSD2の比のみを用いて麻酔深度判定を行ってもよく、脳波の周波数解析、バイスペクトラル解析、スペクトラル解析、エントロピー解析、聴覚誘発電位による解析手段、これらの解析手段の２種以上をSD1/SD2比とさらに組み合わせて麻酔深度の判定を行ってもよい。例えばBISモニターは現在最も普及しているモニターであり、ＢＩＳモニターとSD1/SD2比に基づく本発明の解析手段を組み合わせることで、または、BIS等のアルゴリズムに組み込むことで、より正確に麻酔深度を測定又は判定することができる。 In the present invention, the depth of anesthesia may be determined using only the ratio of SD1 and SD2. EEG frequency analysis, bispectral analysis, spectral analysis, entropy analysis, auditory evoked potential analysis means, 2 of these analysis means The anesthesia depth may be determined by further combining species or more with the SD1 / SD2 ratio. For example, the BIS monitor is the most popular monitor at present, and by combining the BIS monitor and the analysis means of the present invention based on the SD1 / SD2 ratio, or incorporating it into an algorithm such as BIS, the anesthesia depth can be more accurately determined. Can be measured or determined.

　以下、本発明を図面を用いてより詳細に説明するが、本発明がこれら実施例に限定されないことは言うまでもない。 Hereinafter, although the present invention will be described in more detail with reference to the drawings, it goes without saying that the present invention is not limited to these examples.

実施例１：脳波ポアンカレプロット定量化による麻酔深度評価の代表例
全身麻酔中の37歳女性の異なる4段階の麻酔深度（セボフルラン吸入麻酔0.5%, 1%, 2%, 3%）における前頭誘導脳波(10 sec.) とそのポワンカレプロット例を表示する。麻酔深度を増すにつれて、プロットは卵円型から長く伸びる形状に変化し、SD1/SD2は、0.133 (0.5%)、 0.110(1%)、0.077(2%)、0.048(3%)と、段階的に減少した。結果を図５に示す。 Example 1: Representative example of evaluation of depth of anesthesia by quantification of electroencephalogram Poincare plot Frontal lead electroencephalogram at different four stages of anesthesia (sevoflurane inhalation anesthesia 0.5%, 1%, 2%, 3%) during general anesthesia (10 sec.) And its Poincare plot example are displayed. As the depth of anesthesia increases, the plot changes from an oval shape to a shape that extends longer, with SD1 / SD2 being 0.133 (0.5%), 0.110 (1%), 0.077 (2%), and 0.048 (3%). Decreased. The results are shown in FIG.

実施例２：脳波ポワンカレ定量化指標（SD1/SD2）とセボフルラン吸入麻酔深度の関係
20名において、セボフルラン吸入麻酔濃度（0.5%, 1%, 2%, 3%）変化による4段階の麻酔深度を設定し、そのポワンカレプロットをSD1/SD2指標により定量化した。 Example 2: Relationship between EEG Poincare quantification index (SD1 / SD2) and sevoflurane inhalation anesthesia depth
In 20 patients, four stages of anesthesia depth were set according to changes in sevoflurane inhalation anesthesia concentration (0.5%, 1%, 2%, 3%), and the Poincare plot was quantified using the SD1 / SD2 index.

　Repeated measures analysis of variance (ANOVA) , 及び Tukey multiple-comparison testにより、各麻酔深度変化に追従したSD1/SD2指標の有意な変化が認められている（図６、表１）。 Significant changes in the SD1 / SD2 index following each anesthesia depth change were observed by Repeated measurements (analysis) (variance) (ANOVA) (8), and Tukey (multiple-comparison) tests (Fig. 6, Table 1).

実施例３：脳波ポワンカレ定量化指標（SD1/SD2）とパワースペクトルの関係
20名の4段階麻酔深度(セボフルラン0.5%, 1%, 2%, 3%)において、SD1/SD2とパワースペクトルを同時に解析し、そのβ領域(13-3Hz)の周波数成分とδ領域(0.5-4Hz)の周波数成分のパワー比率を算出し、SD1/SD2値との関係を調べたところ、r=0.92の高い相関が認められ、SD1/SD2の麻酔深度測定法としての有用性が支持された(図７)。 Example 3: Relationship between EEG Poincare quantification index (SD1 / SD2) and power spectrum
At the four stages of anesthesia depth (sevoflurane 0.5%, 1%, 2%, 3%) of 20 people, SD1 / SD2 and power spectrum were analyzed at the same time, and the frequency component of the β region (13-3Hz) and the δ region (0.5 -4Hz) frequency component power ratio was calculated and the relationship with the SD1 / SD2 value was examined.A high correlation of r = 0.92 was observed, supporting the usefulness of SD1 / SD2 as an anesthetic depth measurement method. (FIG. 7).

実施例４：遅延時間がポワンカレプロット（SD1/SD2指標）に及ぼす影響について
各麻酔レベルにおいて、ポワンカレプロットにおける時間遅れ座標の遅延時間をサンプリング時間の１～5倍まで変化させると、SD1が上昇し卵円状になる。しかし、遅延時間を5段階に変化させた時、どの麻酔レベルにおいてもSD1/SD2はほぼ比例的に変化した。従って、遅延時間が麻酔深度に依存するSD1/SD2の変化率に及ぼす影響が小さいことがわかる。このことから、時間遅れの設定はサンプリング時間の1-5倍の範囲ではポワンカレプロットの定量的評価にほとんど影響しないことが示唆された。図８は遅延がSD1/SD2値に及ぼす影響を示す。
（実施例1～4をまとめた論文: Kazuko Hayashi, Nobuhiro Mukai, Teiji Sawa. Poincare analysis of the electroencephalogram during sevoflurane anesthesia. Clinical NeuroPhysiology, 2014, in press、DOI information: 10.1016/j.clinph.2014.04.019.) Example 4: Effect of delay time on Poincare plot (SD1 / SD2 index) When the delay time of the time delay coordinate in the Poincare plot is changed to 1 to 5 times the sampling time at each anesthesia level, SD1 increases. Oval shape. However, when the delay time was changed to 5 levels, SD1 / SD2 changed almost proportionally at any anesthesia level. Therefore, it can be seen that the influence of the delay time on the change rate of SD1 / SD2 depending on the depth of anesthesia is small. This suggests that the time delay setting has little effect on the quantitative evaluation of the Poincare plot in the range of 1-5 times the sampling time. FIG. 8 shows the effect of delay on the SD1 / SD2 value.
(Paper summarizing Examples 1-4: Kazuko Hayashi, Nobuhiro Mukai, Teiji Sawa. Poincare analysis of the electroencephalogram during sevoflurane anesthesia. Clinical NeuroPhysiology, 2014, in press, DOI information: 10.1016 / j.clinph.2014.04.019. )

実施例５
本実施例では、BIS専用の脳波電極を前額に装着してBISモニタリングを行うと同時に、導出される脳波を128Hzで収集後、ポアンカレプロット定量解析を施行し、BISモニターから得られる脳波指標（BIS,　SEF95）との比較検討を施行した。 Example 5
In this example, BIS monitoring is performed with a BIS dedicated electroencephalogram electrode attached to the forehead. At the same time, the derived electroencephalogram is collected at 128 Hz, then the Poincare plot quantitative analysis is performed, and the electroencephalogram index obtained from the BIS monitor ( A comparative study with BIS, SEF95) was conducted.

　BIS値は、脳波の周波数解析、バイスペクトル解析等を基に算出した値を、意識レベルに相関する臨床値として0-100までの数値に経験的統計的アルゴリズムにより人為的に当てはめた指標である。覚醒時には100に、深麻酔時には0に近づき、適度な深度の麻酔状態では、40-60とされる。しかし、BIS値は使用する麻酔薬の種類により鎮静度との解離がみられ（Br J Anaesth 2005; 94: 336-40）、また麻酔薬以外のさまざまな因子にも影響を受けることが知られる。 The BIS value is an index artificially applied to values from 0-100 as clinical values correlated with the level of consciousness, based on frequency analysis, bispectrum analysis, etc. . It approaches 100 when awakened, close to 0 during deep anesthesia, and 40-60 under moderately deep anesthesia. However, BIS values are dissociated from sedation depending on the type of anesthetic used (Br J Anaesth 2005; 94: 336-40) and are also known to be affected by various factors other than anesthetics. .

　SEF95は、全体のPower Spectrumの内、95% がその周波数以下に存在する周波数を指す。一般に麻酔が深くなると、SEF95値は小さくなり、麻酔深度を反映する生理的パラメータとして重用される。(Anesth Analg 2004; 98: 1336-40)
　比較検討の結果、解析例より下記の点が示唆された。
1. ポアンカレプロット定量値（SD1/SD2）は、麻酔薬の種別に依らず、spectral edge frequency 95% (SEF95) と強い相関をもつ。
2. ポアンカレプロット定量値（SD1/SD2）は、麻酔深度変化に際して、BIS値やSEF95値より、数十秒、より早く反応する。
3. ポアンカレプロット定量値（SD1/SD2）は、BIS値が真の麻酔深度を反映していないと思える状況においても、安定して推移した。
以下に詳細を記述する。 SEF95 refers to the frequency where 95% of the entire Power Spectrum exists below that frequency. In general, when anesthesia is deepened, the SEF95 value decreases and is used as a physiological parameter reflecting the depth of anesthesia. (Anesth Analg 2004; 98: 1336-40)
As a result of comparative study, the following points were suggested from the analysis example.
1. Poincare plot quantitative value (SD1 / SD2) has a strong correlation with spectral edge frequency 95% (SEF95) regardless of the type of anesthetic.
2. Poincare plot quantitative value (SD1 / SD2) responds faster than BIS and SEF95 values by several tens of seconds when changing anesthesia depth.
3. Poincare plot quantitative values (SD1 / SD2) remained stable even in situations where the BIS value did not seem to reflect the true depth of anesthesia.
Details are described below.

（方法）
現在汎用される麻酔深度計（BISモニター）を用いて全身麻酔中にBISモニタリングと同時に脳波を収集しポワンカレプロット解析を行い、ポワンカレプロット定量値(SD1/SD2)とBIS値等の脳波深度指標との関係を調べた。 (Method)
Using general-purpose anesthesia meters (BIS monitors), brain waves are collected at the same time as BIS monitoring during general anesthesia, and Poincare plot analysis is performed. Echo depth indicators such as Poincare plot quantitative values (SD1 / SD2) and BIS values I investigated the relationship.

　BISモニターより導出した前頭誘導脳波を128Hzサンプリングでアナログデジタル変換（A/D変換）後、Artifactのチェックとbase line driftに対する補正を行い、この前処理脳波時系列信号を、ポワンカレ定量化解析に用いた。これらのポワンカレ定量化解析は、科学技術計算ソフトであるMATLAB ( The Language of Technical Computing, ver. 7.8.0.347 (R2009a) , The MathWorks ^TM, USA )を用い、収集脳波から自動的に施行した。 The frontal brain waves derived from the BIS monitor are converted from analog to digital (A / D conversion) at 128 Hz sampling, then Artifact is checked and base line drift is corrected, and this preprocessed brain wave time series signal is used for Poincare quantitative analysis. It was. These Poincaré quantification analyzes were automatically performed from collected brain waves using MATLAB (The Language of Technical Computing, ver. 7.8.0.347 (R2009a), The MathWorks ^™ , USA), which is scientific and technical calculation software.

　実施例５のポワンカレプロット定量解析値SD1/SD2は、ポワンカレプロットの設定遅延時間と、導出脳波の構成が実施例１～４と異なるため、同じ麻酔深度でも実施例１～４のSD1/SD2計測値と異なる値を示す。今回はBISモニター用の脳波電極を用いており、脳波導出の基準電極の位置が異なるため、導出脳波自身の性状が異なる。また、今回BISモニターより収集した脳波のサンプリング周波数は128Hzであり、ポワンカレプロット解析における最少遅延時間は1/128秒となり、実施例１～４の設定遅延時間の1/256秒の2倍である。このため、ポワンカレプロットパターンは全体的により膨らんだ卵型に近い性状に近づき、同じ麻酔深度でもSD1/SD2計測値より大きい値となっている。 The Poincare plot quantitative analysis value SD1 / SD2 of Example 5 is different from Examples 1 to 4 in the delay time set in the Poincare plot and the configuration of the derived electroencephalogram, so the SD1 / SD2 measurement of Examples 1 to 4 is the same even at the same anesthesia depth. Indicates a value different from the value. This time, an electroencephalogram electrode for a BIS monitor is used, and the position of the reference electrode for deriving the electroencephalogram is different, so the properties of the derived electroencephalogram itself are different. Moreover, the sampling frequency of the EEG collected from the BIS monitor this time is 128Hz, and the minimum delay time in the Poincare plot analysis is 1/128 seconds, which is twice the set delay time of Examples 1 to 4, 1/256 seconds. . For this reason, the Poincare plot pattern is close to the swelled egg shape as a whole, and is larger than the SD1 / SD2 measurement value even at the same anesthetic depth.

1．自動解析の提示
　以下にポワンカレプロットの自動解析例を示す。
提示例は、セボフルラン吸入麻酔濃度を3％から0.5％へ低下させた時の脳波変化のコンピュータ上の自動解析画面である。図９（A）は、セボフルラン吸入麻酔濃度が3％時の、図９（B）はセボフルラン吸入麻酔濃度の設定を3%から0.5％に低下後、約360秒(6分) 後の脳波解析を示した。それぞれ、上段に4秒の脳波信号が、左下にその4秒間の脳波信号のポワンカレプロットパターンと4秒間の脳波信号のポワンカレプロット定量値：SD1/SD2が描出されている。右下には、そのポワンカレプロット定量値（SD1/SD2）の現在値までの推移が示される。深麻酔時（3％、図９（A））には、脳波信号はゆるやかな波を示し、そのポワンカレプロットは、細長い線状パターンとなる。浅麻酔時（0.5％、図９（B））には、脳波は細かな波で形成され、そのポワンカレプロットは、より膨らんだ楕円状のパターンに変化している。ポワンカレプロット定量値（SD1/SD2）は、深麻酔時（3％）の0.12前後から浅麻酔時（0.5% ）には0.3近くまで上昇する。その経過が、パソコンの画面上に秒単位で観察される。 1． Presentation of automatic analysis An example of automatic analysis of the Poincare plot is shown below.
The presented example is an automatic analysis screen on the computer of EEG changes when the sevoflurane inhalation anesthetic concentration is reduced from 3% to 0.5%. Fig. 9 (A) shows an electroencephalogram analysis after about 360 seconds (6 minutes) after the sevoflurane inhalation anesthesia concentration is 3%, and Fig. 9 (B) shows that the sevoflurane inhalation anesthesia concentration is reduced from 3% to 0.5%. showed that. In the upper part, the electroencephalogram signal of 4 seconds is drawn in the upper part, and the Poincare plot pattern of the electroencephalogram signal for 4 seconds and the Poincare plot quantitative value: SD1 / SD2 of the electroencephalogram signal for 4 seconds are drawn in the lower left. In the lower right, the change to the current value of the Poincare plot quantitative value (SD1 / SD2) is shown. During deep anesthesia (3%, FIG. 9 (A)), the electroencephalogram signal shows a gentle wave, and its Poincare plot has an elongated linear pattern. At the time of superanesthetic (0.5%, FIG. 9B), the electroencephalogram is formed as a fine wave, and its Poincare plot changes to a more swollen elliptical pattern. The Poincare plot quantitative value (SD1 / SD2) increases from around 0.12 during deep anesthesia (3%) to nearly 0.3 during shallow anesthesia (0.5%). The progress is observed in seconds on the computer screen.

2. ポワンカレプロット定量値(SD1/SD2)と麻酔深度指標（BIS, SEF95）との比較
＜Sevoflurane吸入麻酔＞　セボフルラン全身麻酔時のBIS値とspectral edge frequency 95% (SEF95) の経時間的推移をポワンカレプロット定量値(SD1/SD2)の推移と共に2症例において提示する。SD1/SD2値は5秒毎の計測値の2点の平均として過去10秒の脳波から10秒毎に算出し、BISモニターにおいて10秒間隔で更新されるBIS及びSEF95値と同じタイミングで比較した。 2. Comparison of Poincare plot quantitative value (SD1 / SD2) and anesthesia depth index (BIS, SEF95) <Sevoflurane inhalation anesthesia> BIS value during sevoflurane general anesthesia and spectral edge frequency 95% (SEF95) over time Presented in 2 cases along with the change of Poincare plot quantitative value (SD1 / SD2). The SD1 / SD2 value is calculated every 10 seconds from the brain wave of the past 10 seconds as the average of 2 points measured every 5 seconds, and compared with the BIS and SEF95 values updated at 10-second intervals on the BIS monitor. .

（CASE1）吸入セボフルラン濃度を3%(深麻酔）から0.5%(浅麻酔)に低下させ、再度3%に戻した時の約27分間の推移を示す。時間軸の一目盛は2分間、計測値はともに5点の移動平均として示した。SD1/SD2は、BISやSSEF95よりも30-40秒程度早く推移し、麻酔深度変化に対して迅速に反応する傾向が明らかである(図１０)。 (CASE1) Shows the transition for about 27 minutes when the inhaled sevoflurane concentration is reduced from 3% (deep anesthesia) to 0.5% (shallow anesthesia) and back to 3%. The scale of the time axis is 2 minutes, and the measured values are shown as a moving average of 5 points. SD1 / SD2 is about 30-40 seconds faster than BIS and SSEF95, and it is clear that there is a tendency to respond quickly to changes in the depth of anesthesia (FIG. 10).

　同症例におけるSD1/SD2値と30秒後のSEF95値、BIS値との関係を、この27分間の10秒毎のデータで示す(図１１)。ポワンカレプロット定量値(SD1/SD2)は、BIS及びSEF95値と良好な相関が認められる。特にSD1/SD2は、30秒先のSEF95と強い相関を示し、SEF95の推移を先行することが示唆される。 The relationship between the SD1 / SD2 value, the SEF95 value after 30 seconds and the BIS value in the same case is shown as data every 10 seconds for 27 minutes (FIG. 11). The Poincare plot quantitative value (SD1 / SD2) has a good correlation with BIS and SEF95 values. In particular, SD1 / SD2 shows a strong correlation with SEF95 30 seconds ahead, suggesting that the transition of SEF95 precedes.

(CASE2）吸入セボフルラン濃度を0.5%から3%に上げて、再度0.5%に戻した時の約20分間の推移を図１２に示す。ここでもSD1/SD2がBISやSEF95よりも数十秒早く反応している。 (CASE2) FIG. 12 shows the transition of about 20 minutes when the inhaled sevoflurane concentration is increased from 0.5% to 3% and then returned to 0.5% again. Again, SD1 / SD2 reacts tens of seconds faster than BIS and SEF95.

＜Propofol静脈麻酔＞　
(Case3)　麻酔をプロポフォール静脈麻酔薬で維持した時の麻酔覚醒前60分間の推移を図１３に示す。覚醒の約50分前と約20分前にBIS値が一過性に上昇する期間が認められるが、SD1/SD2値は安定して推移した。一方、SD1/SD2はSEF95とは強い相関関係を認め、またより早期に反応する傾向を認めた(図１３)。本症例は、SD1/SD2がBISよりも真の麻酔深度を反映する上で有用であることを証明するものではないが、脳波由来の生理的パラメータであるSEF95とは連動して推移しており、BISの欠陥を補うものとしての有用性も期待できる。 <Propofol intravenous anesthesia>
(Case 3) The transition for 60 minutes before anesthesia awakening when anesthesia was maintained with propofol intravenous anesthetic is shown in FIG. SD1 / SD2 values remained stable, although there was a period in which BIS values rose transiently about 50 minutes before and 20 minutes before waking. On the other hand, SD1 / SD2 showed a strong correlation with SEF95, and tended to react earlier (FIG. 13). This case does not prove that SD1 / SD2 is more useful in reflecting the true depth of anesthesia than BIS, but it is linked to SEF95, a physiological parameter derived from brain waves. Also, it can be expected to be useful as a supplement to BIS defects.

実施例６：本発明の方法を用いたセボフルラン、デスフルラン、プロポフォール麻酔深度解析所見
　本発明が特定の麻酔薬に限らず、吸入麻酔薬、静脈麻酔薬を含む複数の麻酔薬において適応できることを確認するため、3種の麻酔薬（計51症例、32351ポイント、吸入麻酔薬：セボフルラン:n=16、デスフルラン:n=22、及び静脈麻酔薬：プロポフォール:n=13）に関して、麻酔中のSD1/SD2と麻酔深度指標の一つであるSEF95(BISモニター由来のSEF95)との関係を回帰分析した。その結果、手術麻酔レベルでは浅麻酔から深麻酔に渡る広い範囲でR²=0.904の強い相関を認めた(図１５)。（線形回帰式：eSEF95=61.9*SD1/SD2+1.3, R²=0.904, RMSE=0.899, p<0.0001）得られた線形回帰式を用いてSD1/SD2から推定した推定SEF95値(eSEF95)とSEF95のBland-Altman plot は、バイアスのない均一なパターンで、系統誤差を認めなかった(図１６)。また各麻酔薬におけるSD1/SD2とSEF95の関係は、セボフルラン：R²=0.916,　デスフルラン：R²=0.903, プロポフォール：R²=0.899といずれも高い相関を示した(図１７、表２)。これらから、本法が特定の麻酔薬に限らず複数の麻酔薬において、手術麻酔深度に応じた脳波の周波数変化を反映することが明らかとなった。より汎用性を有する麻酔深度モニタリングが可能であることを示唆する所見である。 Example 6: Sevoflurane, Desflurane, Propofol Anesthesia Depth Analysis Findings Using the Method of the Present Invention Confirm that the present invention is applicable in multiple anesthetics including but not limited to specific anesthetics, inhalation anesthetics, intravenous anesthetics Therefore, SD1 / SD2 during anesthesia for 3 types of anesthetics (total 51 cases, 32351 points, inhalation anesthetic: sevoflurane: n = 16, desflurane: n = 22, and intravenous anesthetic: propofol: n = 13) And SEF95 (SEF95 derived from BIS monitor), one of the anesthetic depth indices, was regression analyzed. As a result, at the level of surgical anesthesia, a strong correlation of R ² = 0.904 was recognized in a wide range from shallow anesthesia to deep anesthesia (FIG. 15). (Linear regression: eSEF95 = 61.9 * SD1 / SD2 + 1.3, R ² = 0.904, RMSE = 0.899, p <0.0001) Estimated SEF95 value (eSEF95) estimated from SD1 / SD2 using the obtained linear regression The Bland-Altman plot of SEF95 was a uniform pattern with no bias, and no systematic error was observed (FIG. 16). The relationship between SD1 / SD2 and SEF95 in each anesthetic was highly correlated with sevoflurane: R ² = 0.916, desflurane: R ² = 0.903, and propofol: R ² = 0.899 (FIG. 17, Table 2). From these, it became clear that this method reflects the change in the frequency of the electroencephalogram according to the surgical anesthesia depth in a plurality of anesthetics as well as a specific anesthetic. This finding suggests that anesthesia depth monitoring with more versatility is possible.

実施例７：本法解析値の安定性と変動に関しての検討
　本発明の方法を用いて、5秒間の脳波から麻酔深度を測定する際の測定値の変動を、2症例に関して変動係数から検討した。 Example 7: Examination on stability and fluctuation of analysis value of this method Using the method of the present invention, the fluctuation of the measurement value when measuring the depth of anesthesia from the electroencephalogram for 5 seconds was examined from the coefficient of variation for two cases. .

　麻酔中、5秒毎に5秒脳波から連続的にSD1/SD2（青色）とSEF95（茶色）を算出し、それらの30秒間の変動係数（CV：標準偏差／平均）を比較した。脳波は、0.5Hz未満、30Hz以上をフィルタリング除去して前処理したものを両者の検討に用いた。SEF95は、Hammingウインドウを適用したWelch 法によるクロスパワースペクトル密度から求めた(図１８，図１９)。 During anesthesia, SD1 / SD2 (blue) and SEF95 (brown) were calculated continuously from the 5 second brain wave every 5 seconds, and their coefficient of variation (CV: standard deviation / average) was compared for 30 seconds. EEG was pre-processed by filtering out less than 0.5Hz and above 30Hz, and used for both studies. SEF95 was obtained from the cross power spectral density by the Welch method using a Hamming window (FIGS. 18 and 19).

　各図の上段は、5秒毎のSD1/SD2（薄青色）、SEF95（薄茶色）である。10点（1分）の移動平均をそれぞれ濃青色（SD1/SD2）、濃茶色（SEF95）で示した。下段に、SD1/SD2（青色）とSEF95（薄茶）の6点（30秒）の変動係数を示した。全体的にSD1/SD2の変動係数はSEF95の変動係数より小さく、SD1/SD2の算出がスペクトラム解析に基づくSEF95の算出にくらべて安定しており、短時間脳波からでも可能であることを示唆する所見である。 The upper part of each figure is SD1 / SD2 (light blue) and SEF95 (light brown) every 5 seconds. The moving average of 10 points (1 minute) is shown in dark blue (SD1 / SD2) and dark brown (SEF95), respectively. The bottom row shows the coefficient of variation of 6 points (30 seconds) for SD1 / SD2 (blue) and SEF95 (light brown). Overall, the coefficient of variation of SD1 / SD2 is smaller than that of SEF95, suggesting that the calculation of SD1 / SD2 is more stable than the calculation of SEF95 based on spectrum analysis and is possible even from short-term EEGs. It is a finding.

Claims (15)

1. ポワンカレプロット分布の長軸方向への先鋭度を定量化する工程、前記先鋭度に基づき脳波に含まれる乱雑性を評価して麻酔深度を判定する工程を含む、麻酔深度測定法。 An anesthesia depth measurement method including the step of quantifying the sharpness in the major axis direction of the Poincare plot distribution, and the step of evaluating the randomness included in the electroencephalogram based on the sharpness and determining the depth of anesthesia.
2. 以下の工程を含む、麻酔深度測定法
   （１）脳波計により全身麻酔患者の脳波を測定する工程
   （２）前記脳波をフィルタリング処理とアナログ／デジタル変換して、脳波時系列信号を出力する工程（但し、アナログ／デジタル変換及びフィルタリング処理の順序は問わない）
   （３）前記脳波時系列信号と前記脳波時系列信号の遅延時系列信号を、ある時間での脳波電圧Axと前記ある時間から所定時間の遅延後の脳波電圧Ayのペア（Ax, Ay）でｘｙ平面上にプロットする工程、
   （４）前記ｘｙ平面上にプロットされた前記脳波電圧Ayの前記ペア（Ax, Ay）を２つの直交する軸からの変位で定量化し、前記２つの直交する軸に対する前記ペア（Ax, Ay）のそれぞれの変位L1、変位L2に基づいて、前記ペア（Ax, Ay）の前記２つの直交する軸に対する距離に関係づけられた量であるちらばり度V1、ちらばり度V2をそれぞれ求めるちらばり度算出工程、
   （５）前記ちらばり度V1と前記ちらばり度V2との比に基づき、麻酔深度を判定する工程。 Anesthesia depth measurement method including the following steps: (1) A step of measuring the electroencephalogram of a general anesthesia patient by an electroencephalograph (2) A step of outputting the electroencephalogram time-series signal by subjecting the electroencephalogram to filtering processing and analog / digital conversion ( (However, the order of analog / digital conversion and filtering processing does not matter)
   (3) The brain wave time series signal and the delayed time series signal of the brain wave time series signal are expressed as a pair (Ax, Ay) of a brain wave voltage Ax at a certain time and a brain wave voltage Ay after a certain time from the certain time. plotting on the xy plane;
   (4) The pair (Ax, Ay) of the electroencephalogram voltage Ay plotted on the xy plane is quantified by displacement from two orthogonal axes, and the pair (Ax, Ay) with respect to the two orthogonal axes Based on the respective displacements L1 and L2, the dispersion degree V1 and the dispersion degree V2, which are quantities related to the distance of the pair (Ax, Ay) with respect to the two orthogonal axes, are obtained. Degree calculation process,
   (5) A step of determining the depth of anesthesia based on the ratio between the degree of dispersion V1 and the degree of dispersion V2.
3. 前記直交する軸は、前記ｘｙ平面上のそれぞれ長軸（式ｙ＝ｘで表される軸）と短軸（式ｙ＝ｘで表される軸に直交する軸）とである、請求項２に記載の麻酔深度測定法。 The orthogonal axes are a major axis (axis represented by the expression y = x) and a minor axis (axis orthogonal to the axis represented by the expression y = x), respectively, on the xy plane. Anesthesia depth measurement method described in 1.
4. 前記変位L1を前記短軸方向の変位とし、前記変位L2を前記長軸方向の変位とし、前記ちらばり度算出工程は、前記変位L1と前記変位L2を以下の式に従い算出するL1・L2算出工程を含み、L1とL2の比（L1/L2もしくはL2/L1）に基づいて麻酔深度を判定する、請求項２又は３に記載の麻酔深度測定法
   

   

   The displacement L1 is the displacement in the minor axis direction, the displacement L2 is the displacement in the major axis direction, and the dispersion degree calculating step calculates the displacement L1 and the displacement L2 according to the following formulas: L1 / L2 calculation The anesthesia depth measuring method according to claim 2 or 3, comprising a step of determining the depth of anesthesia based on the ratio of L1 and L2 (L1 / L2 or L2 / L1)
   

   
5. 前記変位L1を前記短軸方向の変位とし、前記変位L2を前記長軸方向の変位とし、前記ちらばり度算出工程は、前記変位L1と前記変位L2を以下の式に従い算出するL1・L2算出工程と
   

   

   前記変位L1の標準偏差SD1と前記変位L2の標準偏差SD2とを求めるSD1・SD2算出工程とを含み、
   前記ちらばり度V1、前記ちらばり度V2は、それぞれ前記SD1または前記SD2であり、
   SD1/SD2またはSD2/SD1に基づいて、麻酔深度を判定する、請求項２又は３に記載の麻酔深度測定法。 The displacement L1 is the displacement in the minor axis direction, the displacement L2 is the displacement in the major axis direction, and the dispersion degree calculating step calculates the displacement L1 and the displacement L2 according to the following formulas: L1 / L2 calculation Process and
   

   

   SD1 and SD2 calculation step for obtaining the standard deviation SD1 of the displacement L1 and the standard deviation SD2 of the displacement L2,
   The degree of dispersion V1 and the degree of dispersion V2 are the SD1 or the SD2, respectively.
   The anesthetic depth measurement method according to claim 2 or 3, wherein the depth of anesthesia is determined based on SD1 / SD2 or SD2 / SD1.
6. 前記SD1/SD2または前記SD2/SD1と他の麻酔深度の解析手段を組み合わせて前記麻酔深度を判定する、請求項５に記載の麻酔深度測定法。 The anesthetic depth measurement method according to claim 5, wherein the anesthetic depth is determined by combining the SD1 / SD2 or SD2 / SD1 and another anesthetic depth analysis means.
7. 前記他の麻酔深度の解析手段が、前記脳波の周波数解析、バイスペクトラル解析、スペクトラル解析、エントロピー解析、聴覚誘発電位による解析手段、又はこれらの解析手段の２種以上の組み合わせである、請求項６に記載の麻酔深度測定法。 The other anesthetic depth analysis means is the electroencephalogram frequency analysis, bispectral analysis, spectral analysis, entropy analysis, auditory evoked potential analysis means, or a combination of two or more of these analysis means. Anesthesia depth measurement method described in 1.
8. 全身麻酔患者の脳波を測定する脳波計と、
   前記脳波のアナログ／デジタル(A/D)変換手段と、
   A/D変換後、またはA/D変換前の脳波をフィルタリング処理するフィルタリング手段と、
   前記脳波をフィルタリング処理して得られる脳波時系列信号と前記脳波時系列信号の遅延時系列信号を、ある時間での脳波電圧Axと前記ある時間から所定時間の遅延後の脳波電圧Ayのペア（Ax, Ay）でｘｙ平面上にプロットし、
   前記ｘｙ平面上にプロットされた前記脳波電圧Ayの前記ペア（Ax, Ay）を２つの直交する軸からの変位で定量化し、
   前記２つの直交する軸に対する前記ペア（Ax, Ay）のそれぞれの変位V1、変位V2に基づいて、前記ペア（Ax, Ay）の前記２つの直交する軸に対する距離に関係づけられた量であるちらばり度V1、ちらばり度V2をそれぞれ求めてちらばり度を算出する演算手段と、
   前記ちらばり度V1と前記ちらばり度V2との比に基づき、麻酔深度を判定する麻酔深度判定手段と、
   前記判定結果を表示する表示手段と、
   を備えた麻酔深度測定装置。 An electroencephalograph to measure the electroencephalogram of a general anesthesia patient;
   An analog / digital (A / D) conversion means for the electroencephalogram;
   Filtering means for filtering brain waves after A / D conversion or before A / D conversion;
   An electroencephalogram time-series signal obtained by filtering the electroencephalogram and a delayed time-series signal of the electroencephalogram time-series signal are represented by a pair of electroencephalogram voltage Ax at a certain time and electroencephalogram voltage Ay after a certain time delay from the certain time ( Plot on the xy plane with (Ax, Ay)
   Quantifying the pair (Ax, Ay) of the electroencephalogram voltage Ay plotted on the xy plane by displacement from two orthogonal axes;
   Based on the respective displacements V1 and V2 of the pair (Ax, Ay) relative to the two orthogonal axes, the amount is related to the distance of the pair (Ax, Ay) relative to the two orthogonal axes Calculation means for calculating the degree of dispersion by obtaining the degree of dispersion V1 and the degree of dispersion V2, respectively;
   Based on the ratio of the degree of dispersion V1 and the degree of dispersion V2, the anesthetic depth determination means for determining the depth of anesthesia,
   Display means for displaying the determination result;
   An anesthesia depth measuring device with
9. 前記直交する軸は、前記ｘｙ平面上のそれぞれ長軸（式ｙ＝ｘで表される軸）と短軸（式ｙ＝ｘで表される軸に直交する軸）とであって、
   前記変位L1を前記短軸方向の変位とし、前記変位L2を前記長軸方向の変位とし、
   前記変位L1と前記変位L2を以下の式に従い算出し、
   

   

   前記変位L1の標準偏差SD1と前記変位L2の標準偏差SD2とを求め、
   前記ちらばり度V1、前記ちらばり度V2は、それぞれ前記SD1または前記SD2であり、
   SD1/SD2またはSD2/SD1に基づいて、麻酔深度を判定する麻酔深度判定手段と、
   を備えた請求項８に記載の麻酔深度測定装置。 The orthogonal axes are a long axis (axis represented by the expression y = x) and a short axis (axis orthogonal to the axis represented by the expression y = x), respectively, on the xy plane,
   The displacement L1 is the displacement in the minor axis direction, the displacement L2 is the displacement in the major axis direction,
   Calculate the displacement L1 and the displacement L2 according to the following formula,
   

   

   Obtain the standard deviation SD1 of the displacement L1 and the standard deviation SD2 of the displacement L2,
   The degree of dispersion V1 and the degree of dispersion V2 are the SD1 or the SD2, respectively.
   Based on SD1 / SD2 or SD2 / SD1, anesthetic depth determination means for determining the depth of anesthesia,
   The anesthesia depth measuring device according to claim 8, comprising:
10. 前記麻酔深度判定手段が、前記SD1/SD2または前記SD2/SD1と他の麻酔深度の解析手段を組み合わせて前記麻酔深度を判定する手段である、請求項９に記載の麻酔深度測定装置。 The anesthesia depth measurement device according to claim 9, wherein the depth of anesthesia determination means is means for determining the depth of anesthesia by combining the SD1 / SD2 or SD2 / SD1 and other anesthesia depth analysis means.
11. 前記他の麻酔深度の解析手段が、前記脳波の周波数解析、バイスペクトラル解析、スペクトラル解析、エントロピー解析、聴覚誘発電位による解析手段、又はこれらの解析手段の２種以上の組み合わせである、請求項１０に記載の麻酔深度測定装置。 11. The other anesthetic depth analysis means is the electroencephalogram frequency analysis, bispectral analysis, spectral analysis, entropy analysis, auditory evoked potential analysis means, or a combination of two or more of these analysis means. The anesthesia depth measuring device described in 1.
12. ちらばり度V1と前記ちらばり度V2との比の麻酔深度判定のための使用であって、脳波計により全身麻酔患者の脳波を測定し、
    前記脳波をアナログ／デジタル変換後フィルタリング処理して脳波時系列信号を出力し、
    前記脳波時系列信号と前記脳波時系列信号の遅延時系列信号を、ある時間での脳波電圧Axと前記ある時間から所定時間の遅延後の脳波電圧Ayのペア（Ax, Ay）でｘｙ平面上にプロットし、
    前記ｘｙ平面上にプロットされた前記脳波電圧Ayの前記ペア（Ax, Ay）を２つの直交する軸からの変位で定量化し、
    前記２つの直交する軸に対応する前記ペア（Ax, Ay）のそれぞれの変位L1、変位L2に基づいて、前記ペア（Ax, Ay）の前記２つの直交する軸に対する距離に関係づけられた量であるちらばり度V1、ちらばり度V2をそれぞれ求めてちらばり度V1と前記ちらばり度V2との比を算出することを含む、使用。 The use of the ratio of the degree of dispersion V1 and the degree of dispersion V2 for determining the depth of anesthesia, measuring the electroencephalogram of a general anesthesia patient with an electroencephalograph,
    EEG time-series signal is output by performing filtering after analog / digital conversion of the EEG,
    The brain wave time series signal and the delayed time series signal of the brain wave time series signal are represented on the xy plane by a pair (Ax, Ay) of a brain wave voltage Ax at a certain time and a brain wave voltage Ay after a certain time from the certain time. Plot to
    Quantifying the pair (Ax, Ay) of the electroencephalogram voltage Ay plotted on the xy plane by displacement from two orthogonal axes;
    A quantity related to the distance of the pair (Ax, Ay) to the two orthogonal axes based on the respective displacements L1, L2 of the pair (Ax, Ay) corresponding to the two orthogonal axes And calculating the ratio of the degree of dispersion V1 and the degree of dispersion V2 by obtaining the degree of dispersion V1 and the degree of dispersion V2, respectively.
13. 前記直交する軸は、前記ｘｙ平面上のそれぞれ長軸（式ｙ＝ｘで表される軸）と短軸（式ｙ＝ｘで表される軸に直交する軸とであって、
    前記変位L1を前記短軸方向の変位とし、前記変位L2を前記長軸方向の変位とし、
    前記変位L1と前記変位L2を以下の式に従い算出し、
    

    

    前記変位L1の標準偏差SD1と前記変位L2の標準偏差SD2とを求め、
    前記ちらばり度V1、前記ちらばり度V2は、それぞれ前記SD1または前記SD2である、
    SD1/SD2またはSD2/SD1の麻酔深度判定のための請求項１２に記載の使用。 The orthogonal axes are respectively a long axis (axis represented by the expression y = x) and a short axis (axis orthogonal to the axis represented by the expression y = x) on the xy plane,
    The displacement L1 is the displacement in the minor axis direction, the displacement L2 is the displacement in the major axis direction,
    Calculate the displacement L1 and the displacement L2 according to the following formula,
    

    

    Obtain the standard deviation SD1 of the displacement L1 and the standard deviation SD2 of the displacement L2,
    The dispersion degree V1 and the dispersion degree V2 are the SD1 or the SD2, respectively.
    Use according to claim 12 for determining the depth of anesthesia of SD1 / SD2 or SD2 / SD1.
14. 前記SD1/SD2または前記SD2/SD1を、脳波の周波数解析、バイスペクトラル解析、スペクトラル解析、エントロピー解析、聴覚誘発電位による解析手段からなる群から選ばれる少なくとも1種と組み合わせて麻酔深度の判定を行う、請求項１３に記載の使用。 The SD1 / SD2 or SD2 / SD1 is combined with at least one selected from the group consisting of electroencephalogram frequency analysis, bispectral analysis, spectral analysis, entropy analysis, and auditory evoked potential analysis means to determine the depth of anesthesia 14. Use according to claim 13.
15. 前記SD1/SD2または前記SD2/SD1を麻酔深度の測定のためのアルゴリズムと組み合わせて麻酔深度の判定を行う、請求項１３に記載の使用。 14. Use according to claim 13, wherein the SD1 / SD2 or SD2 / SD1 is combined with an algorithm for measuring depth of anesthesia to determine the depth of anesthesia.

Images