WO2020067505A1 - 呼吸数計測装置 - Google Patents

呼吸数計測装置 Download PDF

Info

Publication number
WO2020067505A1
WO2020067505A1 PCT/JP2019/038359 JP2019038359W WO2020067505A1 WO 2020067505 A1 WO2020067505 A1 WO 2020067505A1 JP 2019038359 W JP2019038359 W JP 2019038359W WO 2020067505 A1 WO2020067505 A1 WO 2020067505A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
patient
pressure
respiratory
respiratory rate
Prior art date
Application number
PCT/JP2019/038359
Other languages
English (en)
French (fr)
Inventor
佑樹 山浦
裕太 森下
章 吉澤
Original Assignee
帝人ファーマ株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 帝人ファーマ株式会社 filed Critical 帝人ファーマ株式会社
Priority to EP19867937.5A priority Critical patent/EP3858409B1/en
Priority to US17/279,215 priority patent/US20220031986A1/en
Priority to CN201980063693.XA priority patent/CN112739403B/zh
Priority to JP2020549479A priority patent/JP7055217B2/ja
Priority to KR1020217009195A priority patent/KR102526829B1/ko
Publication of WO2020067505A1 publication Critical patent/WO2020067505A1/ja

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/021Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes operated by electrical means
    • A61M16/022Control means therefor
    • A61M16/024Control means therefor including calculation means, e.g. using a processor
    • A61M16/026Control means therefor including calculation means, e.g. using a processor specially adapted for predicting, e.g. for determining an information representative of a flow limitation during a ventilation cycle by using a root square technique or a regression analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/0816Measuring devices for examining respiratory frequency
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/087Measuring breath flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7246Details of waveform analysis using correlation, e.g. template matching or determination of similarity
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/0003Accessories therefor, e.g. sensors, vibrators, negative pressure
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/021Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes operated by electrical means
    • A61M16/022Control means therefor
    • A61M16/024Control means therefor including calculation means, e.g. using a processor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/10Preparation of respiratory gases or vapours
    • A61M16/1005Preparation of respiratory gases or vapours with O2 features or with parameter measurement
    • A61M16/101Preparation of respiratory gases or vapours with O2 features or with parameter measurement using an oxygen concentrator
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01DSEPARATION
    • B01D53/00Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
    • B01D53/02Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols by adsorption, e.g. preparative gas chromatography
    • B01D53/04Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols by adsorption, e.g. preparative gas chromatography with stationary adsorbents
    • B01D53/047Pressure swing adsorption
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2505/00Evaluating, monitoring or diagnosing in the context of a particular type of medical care
    • A61B2505/07Home care
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/0057Pumps therefor
    • A61M16/0063Compressors
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/06Respiratory or anaesthetic masks
    • A61M16/0666Nasal cannulas or tubing
    • A61M16/0672Nasal cannula assemblies for oxygen therapy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/10Preparation of respiratory gases or vapours
    • A61M16/14Preparation of respiratory gases or vapours by mixing different fluids, one of them being in a liquid phase
    • A61M16/16Devices to humidify the respiration air
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/20Valves specially adapted to medical respiratory devices
    • A61M16/201Controlled valves
    • A61M16/202Controlled valves electrically actuated
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/0003Accessories therefor, e.g. sensors, vibrators, negative pressure
    • A61M2016/0015Accessories therefor, e.g. sensors, vibrators, negative pressure inhalation detectors
    • A61M2016/0018Accessories therefor, e.g. sensors, vibrators, negative pressure inhalation detectors electrical
    • A61M2016/0021Accessories therefor, e.g. sensors, vibrators, negative pressure inhalation detectors electrical with a proportional output signal, e.g. from a thermistor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/0003Accessories therefor, e.g. sensors, vibrators, negative pressure
    • A61M2016/0027Accessories therefor, e.g. sensors, vibrators, negative pressure pressure meter
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/0003Accessories therefor, e.g. sensors, vibrators, negative pressure
    • A61M2016/003Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/0003Accessories therefor, e.g. sensors, vibrators, negative pressure
    • A61M2016/003Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter
    • A61M2016/0033Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter electrical
    • A61M2016/0039Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter electrical in the inspiratory circuit
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2202/00Special media to be introduced, removed or treated
    • A61M2202/02Gases
    • A61M2202/0208Oxygen
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2205/00General characteristics of the apparatus
    • A61M2205/33Controlling, regulating or measuring
    • A61M2205/3331Pressure; Flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2205/00General characteristics of the apparatus
    • A61M2205/50General characteristics of the apparatus with microprocessors or computers
    • A61M2205/502User interfaces, e.g. screens or keyboards
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2205/00General characteristics of the apparatus
    • A61M2205/50General characteristics of the apparatus with microprocessors or computers
    • A61M2205/52General characteristics of the apparatus with microprocessors or computers with memories providing a history of measured variating parameters of apparatus or patient
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2230/00Measuring parameters of the user
    • A61M2230/40Respiratory characteristics
    • A61M2230/42Rate
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01DSEPARATION
    • B01D2256/00Main component in the product gas stream after treatment
    • B01D2256/12Oxygen
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01DSEPARATION
    • B01D2259/00Type of treatment
    • B01D2259/40Further details for adsorption processes and devices
    • B01D2259/40007Controlling pressure or temperature swing adsorption
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01DSEPARATION
    • B01D2259/00Type of treatment
    • B01D2259/40Further details for adsorption processes and devices
    • B01D2259/40007Controlling pressure or temperature swing adsorption
    • B01D2259/40009Controlling pressure or temperature swing adsorption using sensors or gas analysers
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01DSEPARATION
    • B01D2259/00Type of treatment
    • B01D2259/45Gas separation or purification devices adapted for specific applications
    • B01D2259/4533Gas separation or purification devices adapted for specific applications for medical purposes

Definitions

  • the present disclosure relates to a respiratory rate measuring device used for a pressure fluctuation adsorption type oxygen concentrator.
  • oxygen therapy has been used as one of patients for respiratory diseases such as asthma and obstructive chronic lung disease. This is a therapy in which a patient inhales oxygen gas or oxygen-enriched gas.
  • home oxygen therapy HET: Home Oxygen Therapy
  • QOL quality of life
  • the oxygen concentrator is a device that condenses and discharges about 21% of oxygen present in the air.
  • a pressure fluctuation adsorption type hereinafter, PSA type: Pressure : Swing Adsorption
  • PSA type Pressure : Swing Adsorption
  • the PSA type oxygen concentrator air is taken into an adsorption column filled with an adsorbent for selectively adsorbing nitrogen gas, and an adsorption step of pressurizing the interior of the adsorption column and adsorbing nitrogen gas to the adsorbent is performed.
  • the desorption step of reducing the pressure and discharging the adsorbed nitrogen gas out of the system is repeated.
  • concentrated oxygen gas is generated, and the oxygen concentrator can continuously provide the patient with high-concentration oxygen gas.
  • the method of reducing the pressure in the adsorption cylinder to an atmospheric pressure or lower may be referred to as VPSA or VSA instead of PSA, but the basic principle is the same as PSA. Unify with the name of.
  • COPD chronic obstructive pulmonary disease
  • COPD COPD
  • oxygen therapy respiratory information of a patient, particularly a change in respiratory rate, is a very useful information source for grasping a patient's condition.
  • the respiratory rate can also be measured by using a belly band, etc., but in order to grasp the long-term changes in the patient's condition, the patient must always or regularly wear the belly band. There is a burden on the patient. On the other hand, in general, patients receiving home oxygen therapy often require regular use for at least several hours a day. If the oxygen concentrator has a built-in or attached respiratory rate measuring device, the respiratory rate measuring device can reliably confirm changes in the condition during use of the oxygen concentrator without imposing an additional burden on the patient. it can. Respiratory rate measuring devices are very useful.
  • a fine differential pressure sensor for measuring respiration is attached between the oxygen concentrator and the cannula worn by the patient.
  • Patent Documents 2 and 4 disclose that a respiratory rate and the like can be calculated from a respiratory pattern, but do not describe a specific method.
  • Patent Literature 3 stores the timing at which the pressure waveform changes from falling to rising, counting the respiratory rate at the interval, and multiplying a predetermined detection level rate from the pressure amplitude to exceed the threshold. A method of determining that breathing is performed only when is shown.
  • Patent Document 5 discloses a method of calculating a respiratory rate by FFT (Fast Fourier Transform) processing or TDS processing.
  • FFT Fast Fourier Transform
  • Patent Document 6 discloses a method in which the pressure fluctuation of the oxygen concentrator itself is measured in advance, stored, and subtracted from the detected pressure waveform.
  • the waveform data measured every 100 milliseconds (ms) is sent as it is, for example, and the data amount becomes enormous, and analysis may take time.
  • the method of sending the calculated respiratory data as it is has a drawback that when a respiratory waveform is disturbed due to noise or body movement, a peak value due to the disturb is detected as respiration.
  • the respiratory waveform varies greatly depending on the condition of the patient using the subject, the state of the patient at the time of activity or at bedtime, and the respiratory waveform varies greatly from patient to patient. It is difficult to determine a breath by uniformly setting a threshold for determining a breath.
  • the FFT requires a large amount of calculation, and if the data section used for the calculation is made short, an accurate respiratory cycle cannot be calculated. If the data section is made too long, a slight difference in the respiratory interval between breaths may occur. This makes it difficult to obtain the peak of the respiratory cycle.
  • the TDS process has problems due to the effects of noise and body movement, changes in breathing due to the state of the patient, and differences in the breathing state of each individual patient.
  • the purpose of the respiratory rate measuring device is to solve the above-mentioned problem, and it is used together with the PSA type oxygen concentrator to improve the respiratory rate measuring accuracy.
  • the respiratory rate measuring device is connected to a patient, and is a pressure-concentration adsorption type oxygen concentrating device that condenses oxygen in the air by periodically repeating pressurization and depressurization.
  • a detection unit that detects the pressure in the pipe and / or the gas flow rate in the pipe that supplies the gas to the patient, and outputs the pressure data and / or the gas flow rate data, and performs the patient respiration based on the pressure data and / or the gas flow rate data.
  • the estimating unit estimates the respiratory rate using a waveform indicating at least 10 seconds or more of patient respiratory information data.
  • the estimating unit uses a peak of the autocorrelation coefficient exceeding a predetermined threshold value for estimating the respiratory rate.
  • the threshold value is 0.3 or more and 0.7 or less.
  • the respiratory rate measuring device outputs, when a peak of an autocorrelation coefficient exceeding a predetermined threshold value cannot be obtained, a message indicating that calculation was not possible.
  • the calculation unit is configured to perform a periodic pressure change and / or a periodic pressure change of the oxygen-enriched gas generated by the operation of the oxygen concentrator based on the pressure data and / or the gas flow rate data.
  • the respiratory rate measuring device further includes a storage unit that stores the fluctuation value data measured in advance, and the calculation unit removes the fluctuation value data from the pressure data and / or the previous gas flow rate data. By doing so, it is preferable to extract patient respiration information data.
  • the respiratory rate measuring device is used together with the PSA type oxygen concentrating device, and can improve the measuring accuracy of the respiratory rate without imposing an additional burden on the patient.
  • PSA oxygen concentrator A configuration of a PSA type oxygen concentrator, which is an example of a pressure fluctuation adsorption type oxygen concentrator used together with the respiratory rate measuring device according to one aspect of the embodiment, will be described.
  • FIG. 1 is a diagram showing an example of a schematic configuration of a PSA-type oxygen concentrator.
  • the PSA-type oxygen concentrator 1 has an oxygen generator 11 that takes in air A from outside the PSA-type oxygen concentrator 1 and generates a concentrated oxygen gas.
  • Air A taken into the oxygen generation unit 11 from outside the PSA type oxygen device is compressed by the compressor 111 and sent to the adsorption column 113 via the first switching valve 112.
  • the first switching valve 112 sends compressed air to the adsorption cylinder 113 by communicating any one of the adsorption cylinders 113 with the compressor 111, and opens the other adsorption cylinders to the atmosphere.
  • the adsorption column 113 is filled with an adsorbent for selectively adsorbing nitrogen gas.
  • the compressed air that has passed through the adsorption cylinder 113 has a reduced nitrogen gas concentration and becomes concentrated oxygen gas.
  • the concentrated oxygen gas is stored in the concentrated oxygen buffer tank 115 via the second switching valve 114.
  • the second switching valve 114 connects or disconnects any one of the plurality of adsorption cylinders 113 with the concentrated oxygen buffer tank 115.
  • the oxygen generation unit 11 allows the compressor 111 to communicate with one of the plurality of adsorption cylinders 113 by the first switching valve 112, and connects the adsorption cylinder 113 and the concentrated oxygen buffer tank 115 which communicate with the compressor 111 by the second switching valve 114. Communicate. Therefore, the compressed oxygen buffer tank 115 is communicated with the compressor 111 and any one of the plurality of adsorption cylinders 113, and the generated concentrated oxygen gas is supplied to the concentrated oxygen buffer tank 115.
  • the adsorption cylinder 113 that is not in communication with the compressor 111 is opened to the atmosphere via the first switching valve 112 in a state where the adsorption cylinder 113 is shut off from the concentrated oxygen buffer tank 115 by the second switching valve 114.
  • the pressure in the adsorption cylinder 113 is reduced, and the nitrogen gas adsorbed by the adsorbent is discharged to the outside of the PSA type oxygen concentrator 1.
  • Opening / closing of the first switching valve 112 and the second switching valve 114 is controlled by, for example, a microcomputer unit in a respiratory rate measurement device (not shown).
  • the microcomputer unit can acquire the timing of switching between pressurization and depressurization in the adsorption cylinder 113.
  • the respiratory rate measuring device may be installed inside the oxygen concentrating device 1 or may be installed outside the oxygen concentrating device 1 separately from the oxygen concentrating device 1.
  • the oxygen concentrator 1 may include an oxygen concentration controller that controls the oxygen concentration process including the opening and closing processes of the first and second switching valves.
  • the microcomputer unit can obtain the timing of switching between pressurization and depressurization in the adsorption cylinder 113 from the oxygen concentration control unit.
  • the PSA type oxygen concentrator 1 may be connected to any two or more of the plurality of adsorption columns 113 other than the above-described basic configuration.
  • the PSA-type oxygen concentrator 1 includes a pressure equalizing step for equalizing the pressure of each adsorption column and a purging step for returning a part of the generated oxygen-enriched gas to one of the adsorption columns 113. It may have additional steps.
  • the compressed air is released from the compressor 111 by the first switching valve 112 and released to the atmosphere.
  • the adsorption cylinder 113 that has been released to the atmosphere is connected to the compressor 111 by the first switching valve 112, and shifts to the process of oxygen compression.
  • the plurality of adsorption cylinders 113 alternately repeat compression and release to the atmosphere by the first switching valve 112, so that the concentrated oxygen gas can be continuously supplied.
  • the oxygen flow rate of the concentrated oxygen gas whose pressure has been adjusted by the oxygen generation unit 11 is controlled by an oxygen flow rate control unit 12 including a control valve 121 and a flow meter 122. Supplied to Either the control valve 121 or the flow meter 122 in the oxygen flow control unit 12 may be provided upstream of the flow path, and the oxygen flow control unit 12 may include other components.
  • the PSA-type oxygen concentrator 1 may have a switchable fixed orifice for switching the flow rate instead of the flow meter 122 and the control valve 121.
  • the PSA-type oxygen concentrator 1 may use a flow meter that can be visually observed such as a rotameter as the flow meter 122, and may use a method of adjusting the flow rate by manual operation using a manual flow rate adjusting valve instead of the control valve 121. Alternatively, other flow control methods may be used.
  • the PSA-type oxygen concentrator 1 may be configured without the humidifier 101.
  • FIG. 2 is a diagram showing an example of a schematic configuration of a respiratory rate measuring device.
  • Oxygen generated by the PSA-type oxygen concentrator 1 is supplied to the patient's nasal cavity via a pipe 2 connected to the PSA-type oxygen concentrator 1 and a nasal cannula 3 connected to the pipe 2.
  • the patient is constantly breathing even during oxygen inhalation, and a pressure change caused by the patient's breathing is transmitted to the nasal cannula 3, the pipe 2, and the PSA-type oxygen concentrator 1.
  • the respiratory rate measuring device 4 is connected to the pipe 2 serving as the oxygen supply path in order to acquire the respiratory pressure of the patient.
  • the pipe indicates the pipe including all the pipes between the humidifier 101 (or the oxygen flow controller 12 when the PSA-type oxygen concentrator 1 does not include the humidifier 101) and the nasal cannula 3; May be connected to any position of the pipe.
  • Part or all of the respiratory rate measuring device 4 may be installed inside the PSA type oxygen concentrating device 1 or may be installed outside.
  • the pressure fluctuation added to the respiratory rate measuring device 4 is accompanied by the pressure increase and decrease during the generation of the concentrated oxygen gas. It was found that pressure fluctuations were included.
  • the amplitude of the pressure fluctuation generated during the generation of the concentrated oxygen gas is larger than the amplitude of the respiratory pressure, and the amplitude of the respiratory pressure also decreases due to the pressure loss caused by passing through the oxygen flow path. It is difficult to directly measure a patient's respiratory pressure waveform.
  • the respiratory rate measuring device 4 has a microcomputer unit 7 connected to the pressure sensor 6 and having a calculating unit and an estimating unit.
  • the pressure sensor 6 is preferably a slight differential pressure sensor.
  • the respiratory rate measuring device 4 includes a display unit 8 that displays the measured respiratory rate, for example, a liquid crystal display.
  • the display unit 8 is connected to the microcomputer unit 7 and is controlled by the microcomputer unit 7.
  • the respiratory rate measuring device 4 is electrically connected to, for example, a pressure sensor 6, preferably a slight differential pressure sensor 6, as a detecting unit that detects and outputs the pressure in the pipe. And a microcomputer unit 7.
  • the respiratory rate measuring device 4 may have a pressure smoothing unit including a volume connected to the minute differential pressure sensor 6 and an orifice 5 connecting between the pipe and the volume. The reason for providing the pressure smoothing unit is as follows.
  • the respiratory pressure of a patient is usually about ⁇ 10 to 100 Pa, it is preferable to use a sensor having a range of about ⁇ 100 Pa as the fine differential pressure sensor 6 in order to obtain the respiratory pressure with the respiratory rate measuring device 4.
  • a supply pressure due to the supply of oxygen is constantly generated.
  • a supply pressure due to the supply of oxygen is about 300 Pa even at 1 LPM (liter per minute). Therefore, if one end of the small differential pressure sensor 6 is connected to the oxygen supply path as described above with the other end of the small differential pressure sensor being open to the atmosphere, the measurement range of the small differential pressure sensor 6 will be exceeded. Therefore, it is preferable to apply the pressure after passing through the orifice 5 to the other end of the small differential pressure sensor 6 to acquire the pressure including the breathing information of the patient within the measurement range of the small differential pressure sensor 6.
  • a pressure that does not exceed the measurement range of the fine differential pressure sensor 6 is applied to the other end of the fine differential pressure sensor 6, and the method is not limited to the example of the present embodiment.
  • the other end of the fine differential pressure sensor 6 may be open to the atmosphere as long as the measurement range of the fine differential pressure sensor 6 is larger than the supply pressure by the oxygen supply and the resolution is at a level capable of detecting the pressure fluctuation accompanying the patient's breathing.
  • a pressure sensor for measuring a gauge pressure or an absolute pressure may be used instead of the fine differential pressure sensor 6.
  • the respiratory rate measuring device 4 may use a flow rate sensor instead of the pressure sensor 6.
  • the respiration rate measuring device 4 may have both the pressure sensor 6 and the flow sensor.
  • the pressure sensor 6 and / or the flow sensor detects the pressure in the pipe and / or the gas flow rate in the pipe including the respiration information of the patient in the pipe 2 and outputs the pressure data and / or the gas flow rate data in the pipe. Department.
  • FIG. 3 is a diagram illustrating an example of a configuration block of the microcomputer unit 7.
  • the microcomputer unit 7 having the calculating unit 722 and the estimating unit 723 is connected to the pressure sensor 6, preferably the minute differential pressure sensor 6.
  • the calculating unit 722 and the estimating unit 723 receive the in-pipe pressure data and / or the in-pipe gas flow rate data including the breathing information of the patient detected by the slight differential pressure sensor 6 serving as the detection unit, and execute the processing described below. To obtain the patient's breathing information.
  • the microcomputer unit 7 may be the same microcomputer as the processing unit that processes the oxygen generation function and the display / user interface function of the oxygen concentrator, or may be separated. If they are separated, the switching timing of the cycle T of the PSA is obtained from the microcomputer that processes the oxygen generation function, and is used for calculation.
  • the microcomputer unit 7 includes a storage unit 71 and a processing unit 72.
  • the storage unit 71 includes one or a plurality of semiconductor memories. For example, it has at least one of a nonvolatile memory such as a RAM, a flash memory, an EPROM, and an EEPROM.
  • the storage unit 71 stores a driver program, an operating system program, an application program, data, and the like used for processing by the processing unit 72.
  • the storage unit 71 stores, as a driver program, a device driver program for controlling the small differential pressure sensor 6 and the like as the detection unit.
  • the computer program may be installed in the storage unit 71 from a computer-readable portable recording medium such as a CD-ROM or a DVD-ROM using a known setup program or the like.
  • the program may be downloaded from a program server or the like and installed.
  • the storage unit 71 may temporarily store temporary data relating to a predetermined process.
  • the storage unit 71 stores a threshold 711, a fluctuation value data file 712, and the like used for estimating a respiratory rate.
  • the processing unit 72 includes one or more processors and their peripheral circuits.
  • the processing unit 72 controls the overall operation of the respiratory rate measuring device 4 as a whole, and is, for example, an MCU (Micro Control Unit).
  • MCU Micro Control Unit
  • the processing unit 72 executes processing based on programs (such as an operating system program, a driver program, and an application program) stored in the storage unit 71.
  • the processing unit 72 may execute a plurality of programs (such as application programs) in parallel.
  • the processing unit 72 includes a detection data acquisition unit 721, a calculation unit 722, an estimation unit 723, a respiration rate output unit 724, and the like.
  • processing unit 72 may be implemented in the microcomputer unit 7 as independent integrated circuits, circuit modules, microprocessors, or firmware.
  • Continuous flow 5 LPM data groups acquired when an extension tube is connected 20 m downstream of the respiratory rate measuring device 4 are shown in FIGS. 4, 5 and 6.
  • the continuous flow is one of the concentrated oxygen gas supply methods, and is a method of continuously supplying a constant flow of concentrated oxygen gas.
  • FIG. 4 is a diagram showing pressure data including respiratory information of a patient and PSA pressure data of the PSA type oxygen concentrator 1 when a continuous flow of 5 LPM and an extension tube are connected by 20 m.
  • FIG. 5 is a view showing PSA pressure data of the PSA type oxygen concentrator 1 when there is no patient respiration when the continuous flow is 5 LPM and the extension tube is connected to 20 m, or after the respiratory component is removed.
  • the PSA pressure by the PSA type oxygen concentrator is a periodic pressure change accompanying the cycle of the adsorption cylinder of the PSA type oxygen concentrator 1 generated when oxygen is generated by the PSA type oxygen concentrator 1 described above.
  • the cycle of the PSA pressure waveform coincides with the adsorption cylinder switching cycle of the PSA oxygen concentrator 1.
  • FIG. 6 is a diagram showing the result of performing difference processing on the components of FIG. 4 from the components of FIG. 4 by software.
  • the difference processing is processing for calculating the difference between the two data at an arbitrary time.
  • the respiratory rate measuring device 4 can remove the PSA pressure component and detect the respiratory pressure of the patient's respiratory model.
  • the respiratory rate measuring apparatus 4 allows the patient to breathe even under the condition that the oxygen is continuously inhaled from the PSA type oxygen concentrator 1 through the nasal cannula and the extension tube 20 m at 5 LPM. Information can be obtained.
  • the method of removing the fluctuation value data such as the PSA pressure component and extracting the patient respiration information data is described in Patent Document 6 by measuring, storing and detecting the pressure fluctuation of the oxygen concentrator itself in advance.
  • a method of subtracting from a pressure waveform or a method of subtracting a measured value in real time may be used.
  • the calculating unit 722 determines that the pressure at a certain time t is equal to the flow rate value Y (t). , Y (t) and Y (t ⁇ T), Y (t ⁇ 2T),..., Y (t ⁇ nT) (n is a predetermined integer) X (t) ) Is calculated to extract respiratory information.
  • the original respiratory waveform from which the PSA pressure component has been removed can be reproduced well.
  • FIG. 7 is a conceptual diagram of the arithmetic processing according to the present embodiment.
  • (I) is a respiratory waveform
  • (II) is a model waveform of a PSA waveform.
  • T is the cycle of the PSA waveform, which coincides with the switching cycle of the suction cylinder.
  • (III) is a waveform obtained by adding (I) and (II), and corresponds to the pressure measured in the piping.
  • (IV) is obtained by dividing (III) for each cycle T and superimposing them. As long as the respiratory cycle and the PSA cycle do not completely match, the respiratory waveform appears randomly in the waveform cut at the cycle T.
  • (V) is the result of averaging the waveforms superimposed in (IV).
  • fluctuation value data X (t) indicating a periodic pressure change and / or flow rate change of the oxygen-enriched gas generated by the operation of the PSA-type oxygen concentrator.
  • V are estimated by superposition and averaging
  • a simple moving average for five periods is used, where T is one period.
  • (VI) is a waveform obtained by subtracting (V) from the last T portion of the waveform of (III). It can be seen that the original respiratory waveform (I) has been accurately reproduced.
  • the method of selecting T may be T from the time when one of the plurality of adsorption cylinders starts communicating with the compressor to the time when another adsorption cylinder switches to the state communicating with the compressor.
  • the method of selecting T is defined as T from the time when the state in which one adsorber is in communication with the compressor starts, the time when the other adsorber is in communication with the compressor, and the time when the first adsorber is switched to the state in which it communicates with the compressor again Is also good. Further, as another example, T may be selected as a multiple of each switching time.
  • X (t) is the value after the averaging process at time t
  • Y (t) is the measured value of the pressure at time t.
  • a i is a weighting coefficient in the weighted average, and an arbitrary real number can be selected.
  • n and a i an averaging process with faster convergence, such as a FIR (Finite Impulse Response) filter, is also possible.
  • equation (1) appears to be simply an equation for numerically calculating the time average of Y (t).
  • T is selected to be a sufficiently small value with respect to the fluctuation period of Y (t), whereas in the present embodiment, it is important to select T based on the adsorption / desorption period of the PSA process. is there.
  • This is not a simple averaging process of the measured value Y (t) but a waveform of Y (t) over the period T as one unit, which is traced back by an integral multiple of T. It means taking the average value.
  • the average value of the waveform it is possible to accurately estimate the PSA pressure having a characteristic that appears periodically at intervals of T.
  • the method of measuring pressure was shown as a respiratory rate measurement unit.However, since the flow rate flowing through the piping changes in accordance with the fluctuation of pressure, a method of measuring flow rate instead of pressure may also be used. Can be. Further, a method in which the pressure value and the flow rate value are combined or switched depending on the operating conditions and environmental conditions of the device can be used.
  • the respiration information obtained by the arithmetic processing is so-called raw data that represents real-time information of respiration such as pressure and flow rate as a waveform.
  • This data may be recorded or transmitted as it is, but the data amount becomes enormous and the analysis takes time. Therefore, it is desirable to automatically calculate, record, and / or transmit the respiratory rate data inside the respiratory rate measuring device 4.
  • the respirable waveform that can be detected has many noise components due to pressure fluctuations due to airflow and swaying of the cannula.
  • the peak detection for calculating the timing at which the pressure changes from falling to rising or the method for detecting the timing at which the pressure becomes equal to or greater than the threshold value cannot properly calculate the respiratory cycle and respiratory rate.
  • the respiratory waveforms to be counted appear as peaks represented by a1, a2, and a3 in the figure and valleys that appear immediately after the peaks, and the other portions are noise components. Even if an attempt is made to detect this peak as the timing at which the pressure changes from rising to falling, it is highly likely that part b in the figure will be detected.
  • the autocorrelation coefficient between the original waveform and the waveform shifted by a time ⁇ t from the original waveform is obtained, and ⁇ t is changed so that the autocorrelation coefficient becomes a peak. It has been found that the respiration rate can be detected with high detection power by obtaining.
  • the inventor has found a respiration rate estimating method for estimating a respiration rate per predetermined time from patient respiration information.
  • the autocorrelation coefficient R can be calculated by the following formula.
  • ⁇ t is the amount of time shift
  • t 0 is the data acquisition interval
  • n is the number of data used for one calculation of the autocorrelation coefficient
  • f (t) is the patient respiration information data at time t.
  • n values are obtained at intervals of t 0 .
  • ⁇ and ⁇ are the average value and the standard deviation of f (t), but in the actual calculation, the average and the standard deviation of the acquired n values of f (t) may be used.
  • the calculation method of the autocorrelation coefficient of the present embodiment is an example, and another calculation method of the autocorrelation coefficient may be used.
  • the autocorrelation coefficient R calculated in this embodiment includes a normalization coefficient , The value is defined in the range of -1.0 to +1.0.
  • a waveform in which it was difficult to determine a respiratory cycle by threshold or peak detection in raw patient respiratory information data can easily obtain a respiratory cycle by obtaining an autocorrelation coefficient. I understand.
  • the data section of f (t) used for calculating the autocorrelation coefficient is determined from the range of the respiratory cycle to be calculated. As a result of intensive studies by the inventor, it has been found that the data section of f (t) needs to be at least twice the minimum respiratory cycle to be calculated.
  • the threshold value of the autocorrelation coefficient for determining the peak of the autocorrelation coefficient needs to be set to at least 0.3 or more, preferably between 0.3 and 0.7. . At lower thresholds, the possibility of erroneous determination of an accidental increase in autocorrelation value due to noise or disturbance of the baseline as a peak increases, and at higher threshold values, the increase in autocorrelation value due to the respiratory cycle can be overlooked. Turned out to be high.
  • the respiration rate is calculated as the reciprocal of the respiration interval of the value of ⁇ t when the autocorrelation coefficient is equal to or more than a certain value and takes a peak value.
  • ⁇ t at the point P1 of the first peak is the respiratory interval
  • dividing one minute by the respiratory interval is the respiratory rate per minute.
  • Depending on the respiratory waveform, there may be a plurality of peaks having a value equal to or more than a certain value as shown in FIG. 8, but the peak P2 on the right side is a double cycle peak that appears when ⁇ t is shifted by a plurality of breaths.
  • the leftmost that is, the peak P1 having the smallest ⁇ t, which is estimated to correspond to the basic cycle of respiration, should be used.
  • the information on the respiratory rate can be automatically and accurately calculated from the respiratory waveform acquired by the respiratory rate measuring unit connected to the concentrator.
  • the estimation unit 723 calculates the respiration rate. The accuracy can be further improved by not performing. The determination on the presence or absence of breathing is not tested on the determination based on the variance of the patient respiration information data f (t), and another determination method may be used.
  • the variance ⁇ 2 of the patient respiration information data f (t) is calculated by the following equation.
  • n is the number of data used for one calculation of the autocorrelation coefficient
  • f (t) is the value of the patient respiration information data at time t.
  • ⁇ and ⁇ 2 are the average value and variance of f (t).
  • the threshold for dispersion determination in the case of a PSA-type oxygen concentrator is: When the respiratory rate is less than 8: ⁇ 2 ⁇ 18 (Pa 2 ) For breathing rate 8-10: ⁇ 2 ⁇ 18 (Pa 2 ) When the respiratory rate is 11 or more: ⁇ 2 ⁇ 7 (Pa 2 )
  • the value may vary depending on the oxygen supply device used. Since there are variations depending on parameters such as the respiratory rate, the threshold value for the dispersion determination needs to be appropriately set depending on the oxygen supply device used.
  • FIG. 9 is a flowchart illustrating an example of a process for extracting patient respiration information data.
  • the process of extracting the patient respiration information data shown in FIG. 9 is executed by the microcomputer unit 7 according to a computer program stored in the storage unit 71 in advance.
  • the detection data acquisition unit 721 acquires pressure data Y (t) from the pressure sensor 6, which is a detection unit (ST101).
  • Arithmetic unit 722 calculates an average value X (t) for n cycles of Y (t) (ST102).
  • Arithmetic unit 722 extracts patient respiration information data f (t) by calculating the difference between Y (t) and X (t) (ST103).
  • the arithmetic unit 722 When using the PSA pressure fluctuation measured and stored in advance, the arithmetic unit 722 does not perform the processing of ST102, and after the processing of ST101, calculates the difference between Y (t) and the PSA pressure fluctuation measured and stored in advance as ST103. Extracts patient respiration information data f (t).
  • the process of estimating the respiratory rate is repeatedly performed. For example, every time n pieces of measurement data of the pressure Y (t) are acquired at intervals of t 0 , the patient respiratory information data f (t ) Is extracted and updated.
  • the patient respiration information data f (t) may be extracted and updated at predetermined intervals.
  • FIG. 10 is a flowchart illustrating an example of a process of estimating a respiratory rate based on patient respiratory information data.
  • the processing for estimating the respiratory rate shown in FIG. 10 is executed by the microcomputer unit 7 according to a computer program stored in the storage unit 71 in advance.
  • the estimating unit 723 calculates an average value ⁇ and a standard deviation ⁇ of the patient respiration information data f (t) (ST201).
  • the estimating unit 723 determines whether or not the variance ⁇ 2, which is the square of the standard deviation ⁇ , is equal to or greater than a predetermined threshold TH D (ST202). If the variance ⁇ 2 is less than the threshold value TH D (ST202: NO), the respiratory rate output unit 724 outputs a message indicating that calculation was not possible (ST213), and ends the process.
  • the estimating unit 723 sets the time ⁇ t to zero (0) as the time shift amount (ST203).
  • the estimating unit 723 increments the data acquisition interval t 0 at time ⁇ t (ST204). Estimating section 723 calculates autocorrelation coefficient R ( ⁇ t) at time ⁇ t of patient respiration information data f (t) (ST205). Estimating section 723 writes calculated autocorrelation coefficient R ( ⁇ t) to storage section 71 (ST206). Estimating section 723 determines whether or not time ⁇ t has reached nt 0 (ST207). When the time ⁇ t has not reached nt 0 , the process returns to ST203, and the estimating unit 723 repeats the process (ST207: NO).
  • estimation section 723 reads out autocorrelation coefficient R ( ⁇ t) stored in storage section 71 (ST208).
  • the estimating unit 723 compares the read autocorrelation coefficients R ( ⁇ t), and determines whether there is a maximum autocorrelation coefficient Max (R ( ⁇ t)), that is, whether there is a peak autocorrelation coefficient R ( ⁇ t). Is determined (ST209). Specifically, it is determined whether Max (R ( ⁇ t)) is equal to or greater than a predetermined threshold TH C.
  • the estimating unit 723 determines the maximum autocorrelation coefficient Max (R ( ⁇ t)).
  • the time ⁇ t is set as a breathing interval (ST210). Further, estimation section 723 estimates the respiration rate from respiration interval ⁇ t (ST211). Respiratory rate output section 724 outputs a respiratory rate signal (ST212), and the process ends. For example, the display unit 8 to which the respiratory rate signal is input displays the respiratory rate.
  • the respiratory rate output unit 724 When the maximum autocorrelation coefficient Max (R ( ⁇ t)) is less than the predetermined threshold TH C (ST209: NO), the respiratory rate output unit 724 outputs a message indicating that calculation was not possible (ST213). , And the process ends. While the respiratory rate measuring device 4 is operating, the process of estimating the respiratory rate is repeatedly performed, and the respiratory rate measuring device 4 can update the respiratory rate and display it on the display unit 8. The display on the display unit 8 may be performed at predetermined intervals.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Pulmonology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Emergency Medicine (AREA)
  • Anesthesiology (AREA)
  • Hematology (AREA)
  • Physics & Mathematics (AREA)
  • Physiology (AREA)
  • Surgery (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Chemical & Material Sciences (AREA)
  • Oil, Petroleum & Natural Gas (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

PSA式酸素濃縮装置に用いられ、計測精度を向上可能な呼吸数計測装置であって、患者と接続され、且つ、加圧・減圧を周期的に繰り返すことで空気中の酸素を濃縮する圧力変動吸着方式の酸素濃縮装置(1)から酸素濃縮ガスを患者に供給する配管(2)内の管内圧力及び/又は管内気体流量を検知して、圧力データ及び/又は気体流量データを出力する検知部(6)と、圧力データ及び/又は気体流量データに基づいて、患者呼吸情報データを抽出する演算部(722)と、患者呼吸情報データに基づいて所定時間当たりの呼吸数を推定する推定部(723)と、を有し、推定部(723)は、時間Δtを変化させながら、患者呼吸情報データと、患者呼吸情報データから時間Δtだけずらしたデータとの自己相関係数を求め、自己相関係数がピークとなる時間Δtを呼吸間隔として、呼吸数を推定する、ことを特徴とする呼吸数計測装置(4)。

Description

呼吸数計測装置
 本開示は、圧力変動吸着式酸素濃縮装置に用いられる呼吸数計測装置に関する。
 従来、喘息、閉塞性慢性肺疾患等の呼吸器疾患患者に対するひとつとして酸素療法が行われている。これは、酸素ガスや酸素濃縮ガスを患者に吸入させる療法である。近年では、患者QOL(QOL:Quality of Life)向上を目的に自宅や施設等で酸素吸入をする在宅酸素療法(HOT:Home Oxygen Therapy)が主流になってきており、酸素供給源としては酸素濃縮装置が主に使用されている。
 酸素濃縮装置とは、空気中に存在する約21%の酸素を濃縮して排出する装置である。酸素濃縮装置の多くは圧力変動吸着式(以下、PSA式:Pressure Swing Adsorption)が一般的に用いられる。
 PSA式酸素濃縮装置では、窒素ガスを選択的に吸着する吸着剤が充填された吸着筒に空気を取り込み、吸着筒内を加圧し窒素ガスを吸着剤に吸着させる吸着工程と、吸着筒内を減圧し吸着された窒素ガスを系外に排出する脱着工程が繰り返えされる。吸着工程と脱着工程が繰り返されることにより濃縮酸素ガスが生成され、酸素濃縮装置は患者に連続的に高濃度酸素ガスを提供することができる。なお、脱着工程に於いて吸着筒内を大気圧以下まで減圧する方法は、PSAではなくVPSAやVSA等と呼称されることがあるが、基本的な原理はPSAと同一であり、ここではPSAの呼称で統一する。
 在宅酸素療法を受ける患者の主疾患は慢性閉塞性肺疾患(以下、COPD:Chronic Obstructive Pulmonary Disease)である。COPDとは気管支の狭窄や肺胞壁の破壊により、咳・痰や労作時呼吸困難の症状を呈する不可逆性疾患である。
 COPDの症状が悪化すると息切れや呼吸数の増加が見られ、COPDの急性増悪と呼ばれる「安定期の治療の変更あるいは追加が必要な状態」まで症状が悪化することもある。COPDの急性増悪が起こると、患者は入院するケースが多く、呼吸不全に陥ることや生命の危機に直面することもある。また、患者が退院できたとしても入院以前より安定期の症状が悪化し、入退院を繰り返すことも珍しくない。
 COPDに於いては、患者の急性増悪の予兆又は初期症状がなるべく早い段階で察知されて、症状が入院する程度まで悪化する前に、患者に早期治療を施すことが非常に重要であり、在宅酸素療法において患者の呼吸情報、特に呼吸数の変化は患者の病態を把握するうえで非常に有益な情報源となる。
 呼吸数は、腹バンド等を使用することによっても測定できるが、患者の病態の長期的な変化を把握するためには、患者は腹バンドを常時、あるいは定期的に忘れずに着用する必要がある等、患者が負担を強いられる。一方で一般的に在宅酸素療法を受ける患者は、少なくとも1日数時間の定期的な使用が求められることが多い。酸素濃縮器に呼吸数計測装置が内蔵、あるいは付属していれば、呼吸数計測装置によって患者に追加の負担を強いることなく、かつ確実に酸素濃縮装置使用中の病態の変化を確認することができる。呼吸数計測装置は非常に有用である。
 酸素濃縮装置に内蔵可能な、患者の呼吸数を取得する方法として、特許文献1のように、酸素濃縮装置と患者が装着しているカニューラまでの間に呼吸計測用の微差圧センサを取り付け酸素吸入中の患者呼吸圧を計測し、患者呼吸圧を記録媒体に記録したり、通信データとして送信する方法がある。
 また、特許文献2、4には、呼吸パターンから呼吸数などが算出できることが示されているが、具体的な方法については述べられていない。
 特許文献3には、圧力波形が下降から上昇に転じたタイミングを記憶しておき、その間隔で呼吸数をカウントする方法や、圧力の振幅から所定の検出レベル率を乗じてその閾値を超えたときのみ呼吸と判定する方法が示されている。
 特許文献5には、FFT(FFT:Fast Fourier Transform)処理やTDS処理により呼吸数を算出する方法が開示されている。
 特許文献6には、酸素濃縮装置自身の圧力変動を予め計測し、記憶しておいて、検出した圧力波形から差し引く方法が開示されている。
特開平6-190045号公報 特開2001-286566号公報 特開平7-96035号公報 特表2015-85191号公報 特表2011-518016号公報 国際公開2018-180392号公報
 計算された呼吸データをそのまま送る方法は、例えば100ミリ秒(ms)毎に測定した波形データをそのまま送信することになり、データ量が膨大となり、解析に時間を要することがある。
 計算された呼吸データをそのまま送る方法は、ノイズや体動などにより呼吸波形に乱れが発生した場合、その乱れによるピーク値なども呼吸として検出してしまう欠点がある。呼吸波形は、使用している患者の病態や、活動時、就寝時などの状態によって大きく変化し、さらには患者ごとに呼吸波形は大きく異なる。呼吸を判定するための閾値を一律に設けて、呼吸を判定することは困難である。
 さらに、FFTは計算量が膨大になると共に、演算に用いるデータの区間を短くとると正確な呼吸周期が算出できなくなり、データの区間を長くとりすぎると、呼吸毎の呼吸間隔のわずかなずれの影響で呼吸周期のピークが得られにくくなる。TDS処理については、ノイズや体動の影響、患者の状態による呼吸の変化や、患者個人ごとの呼吸状態の違いに起因する課題がある。
 呼吸数計測装置の目的は、上述の問題を解決するためになされたものであり、PSA式酸素濃縮装置とともに用いられ、呼吸数の計測精度を向上させることを可能することにある。
 実施形態の一側面に係る呼吸数計測装置は、患者と接続され、且つ、加圧・減圧を周期的に繰り返すことで空気中の酸素を濃縮する圧力変動吸着方式の酸素濃縮装置から酸素濃縮ガスを患者に供給する配管内の管内圧力及び/又は管内気体流量を検知して、圧力データ及び/又は気体流量データを出力する検知部と、圧力データ及び/又は気体流量データに基づいて、患者呼吸情報データを抽出する演算部と、患者呼吸情報データに基づいて所定時間当たりの呼吸数を推定する推定部と、を有し、推定部は、時間Δtを変化させながら、患者呼吸情報データと、患者呼吸情報データから時間Δtだけずらしたデータとの自己相関係数を求め、自己相関係数がピークとなる時間Δtを呼吸間隔として、呼吸数を推定する、ことを特徴とする。
 更に、実施形態の一側面に係る呼吸数計測装置は、推定部が、少なくとも10秒以上の患者呼吸情報データを示す波形を用いて、呼吸数を推定する、ことが好ましい。
 更に、実施形態の一側面に係る呼吸数計測装置は、推定部が、所定の閾値を超えた自己相関係数のピークを呼吸数の推定に利用する、ことが好ましい。
 更に、実施形態の一側面に係る呼吸数計測装置は、閾値が、0.3以上且つ0.7以下の値である、ことが好ましい。
 更に、実施形態の一側面に係る呼吸数計測装置は、所定の閾値を超えた自己相関係数のピークを取得できない場合、計算不能であった旨を出力する、ことが好ましい。
 更に、実施形態の一側面に係る呼吸数計測装置は、演算部が、圧力データ及び/又は気体流量データに基づいて、酸素濃縮装置の動作により発生する酸素濃縮ガスの周期的な圧力変化及び/又は流量変化を示す変動値データを推定し、変動値データを圧力データ及び/又は気体流量データから除去することによって、患者呼吸情報データを抽出する、ことが好ましい。
 更に、実施形態の一側面に係る呼吸数計測装置は、予め測定した変動値データを記憶する記憶部を更に有し、演算部は、変動値データを圧力データ及び/又は前気体流量データから除去することによって、患者呼吸情報データを抽出する、ことが好ましい。
 本実施形態によれば、呼吸数計測装置は、PSA式酸素濃縮装置とともに用いられ、患者に追加の負担を強いることなく、呼吸数の計測精度を向上させることを可能する。
 本発明の目的及び効果は、特に請求項において指摘される構成要素及び組み合わせを用いることによって認識され且つ得られるだろう。前述の一般的な説明及び後述の詳細な説明の両方は、例示的及び説明的なものであり、特許請求の範囲に記載されている本発明を制限するものではない。
PSA式酸素濃縮装置の概略構成の一例を示す図である。 呼吸数計測装置の概略構成の一例を示す図である。 マイコン部の一例を示すブロック図である。 連続流5LPM、延長チューブを20m追加した条件における、患者呼吸情報を含んだ圧データを示す図である。 連続流5LPM、延長チューブを20m追加した条件における、演算部により抽出したPSA圧データを示す図である。 連続流5LPM、延長チューブを20m追加した条件における、差分処理後の患者呼吸情報を示す図である。 患者呼吸による圧力変動とPSA圧を含んだ圧力データから、患者呼吸情報とPSA圧を分離する原理を示した概念図である。 差分処理後の患者呼吸情報から、自己相関係数を求めた結果を表す図である。 患者呼吸情報データを抽出する処理の一例を示すフローチャート図である。 患者呼吸情報データに基づいて、呼吸数を推定する処理の一例を示すフローチャート図である。
 以下、実施形態の一側面に係る呼吸数計測装置について、図を参照しつつ説明する。但し、本開示の技術的範囲はそれらの実施の形態に限定されず、特許請求の範囲に記載された発明とその均等物に及ぶ点に留意されたい。なお、以下の説明及び図において、同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
[PSA式酸素濃縮装置]
 実施形態の一側面に係る呼吸数計測装置と共に使用される、圧力変動吸着方式の酸素濃縮装置の一例である、PSA式酸素濃縮装置の構成を説明する。
 図1は、PSA式酸素濃縮装置の概略構成の一例を示す図である。
 PSA式酸素濃縮装置1は、PSA式酸素濃縮装置1の外部から空気Aを取り込んで濃縮酸素ガスを生成する酸素生成部11を有する。
 PSA式酸素装置外部から酸素生成部11に取り込まれた空気Aは、コンプレッサ111で圧縮され、第1切替弁112を介して吸着筒113に送られる。第1切替弁112は、複数の吸着筒113のいずれかとコンプレッサ111を連通することで吸着筒113に圧縮空気を送り込むとともに、その他の吸着筒を大気に開放する。
 吸着筒113には窒素ガスを選択的に吸着する吸着剤が充填されている。吸着筒113を通過した圧縮空気は窒素ガス濃度が低下し、濃縮酸素ガスとなる。濃縮酸素ガスは第2切替弁114を介して濃縮酸素バッファタンク115に貯蔵される。第2切替弁114は複数の吸着筒113のいずれかと濃縮酸素バッファタンク115とを連通または遮断する。
 酸素生成部11は、第1切替弁112によりコンプレッサ111と複数の吸着筒113のいずれかを連通させるとともに、第2切替弁114により該コンプレッサ111と連通した吸着筒113と濃縮酸素バッファタンク115を連通させる。したがって、コンプレッサ111と複数の吸着筒113のいずれかと濃縮酸素バッファタンク115が連通され、生成された濃縮酸素ガスが濃縮酸素バッファタンク115に供給される。他方、コンプレッサ111と連通しない吸着筒113は、第2切替弁114により濃縮酸素バッファタンク115からと遮断された状態で、第1切替弁112を介して大気開放される。これにより、吸着筒113内が減圧し、吸着剤に吸着された窒素ガスをPSA式酸素濃縮装置1外部に放出する。
 第1切替弁112及び第2切替弁114の開閉は、例えば、不図示の呼吸数計測装置内のマイコン部により制御される。マイコン部は吸着筒113内の加圧と減圧を切換えるタイミングを取得することが可能である。呼吸数計測装置は、酸素濃縮装置1の内部に設置されていてもよく、酸素濃縮装置1とは別体に酸素濃縮装置1の外部に設置されていてもよい。
 又、別の例として、酸素濃縮装置1は、第1及び第2切替弁の開閉処理を含む酸素濃縮処理を制御する酸素濃縮制御部を有してもよい。マイコン部は酸素濃縮制御部から吸着筒113内の加圧と減圧を切換えるタイミングを取得することができる。
 PSA式酸素濃縮装置1は、上述の基本構成以外に複数の吸着筒113のいずれか2つ以上の吸着筒を接続してもよい。PSA式酸素濃縮装置1は、それぞれの吸着筒の圧力を均一化する均圧工程や、生成された酸素濃縮ガスの一部を、複数の吸着筒113の内いずれかに還流するパージ工程などの付加的な工程を有していてもよい。
 通常、吸着筒113に圧縮空気が送られた後は、第1切替弁112によって該コンプレッサ111と遮断された状態になり大気解放される。反対に、大気解放されていた吸着筒113は、第1切替弁112によって該コンプレッサ111に連結され、酸素圧縮の過程に移行する。このように、第1切替弁112によって複数の吸着筒113が交互に圧縮と大気解放を繰り返すことで、連続的に濃縮酸素ガスを供給することができる。
 濃縮酸素生成に伴う吸着筒113内部の圧力変化は非常に大きいため、吸着筒113の切り替えに伴って、吸着筒113下流の酸素ガス流路中の圧力には周期的な圧力変動が生じる。濃縮酸素バッファタンク115に貯蔵された濃縮酸素ガスは調圧弁116により該圧力変動が減衰されるように調整される。
 酸素生成部11で圧力調整された濃縮酸素ガスは、コントロールバルブ121と流量計122から成る酸素流量制御部12で酸素流量が制御され、加湿器101を介して酸素供給口13より酸素濃縮装置外に供給される。酸素流量制御部12においてコントロールバルブ121と流量計122はどちらが流路の上流に備えられていてもよく、また酸素流量制御部12にはその他の構成が含まれていても構わない。
 PSA式酸素濃縮装置1は、流量計122とコントロールバルブ121の代わりに、流量を切り替えるための切換え式の固定オリフィスを有してもよい。PSA式酸素濃縮装置1は、流量計122としてロータメータなどの目視可能な流量計を用いるとともにコントロールバルブ121の代わりに手動式の流量調整弁を用いて手動操作により流量を調整する方式を用いてもよく、その他の流量調整方式を用いてもよい。PSA式酸素濃縮装置1は、加湿器101を備えない構成とすることも可能である。
[呼吸数計測装置]
 本発明に係る呼吸数計測装置の概略構成の一例を説明する。
 図2は、呼吸数計測装置の概略構成の一例を示す図である。
 PSA式酸素濃縮装置1で生成された酸素は、PSA式酸素濃縮装置1に接続された配管2、配管2に接続された鼻カニューラ3を経由して患者の鼻腔に供給される。一方、患者は酸素吸入中も常時呼吸をしており、患者の呼吸によって生じる圧変化が鼻カニューラ3や配管2、PSA式酸素濃縮装置1の方へ伝搬される。
 本実施形態においては、患者の呼吸圧を取得すべく、酸素供給経路となる配管2に呼吸数計測装置4が接続される。配管は加湿器101(PSA式酸素濃縮装置1が加湿器101を備えない場合は酸素流量制御部12)と鼻カニューラ3の間の配管を全て含んだものを指しており、呼吸数計測装置4は配管のどの位置に接続されていてもよい。
 呼吸数計測装置4の一部ないし全部は、PSA式酸素濃縮装置1の内部に設置されていてもよく、外部に設置されていてもよい。
 本願発明者らが鋭意検討を行った結果、調圧弁116によって圧力変動を減衰した後であっても、呼吸数計測装置4に付加される圧力変動には濃縮酸素ガス生成時の加減圧に伴う圧力変動が含まれることが明らかとなった。濃縮酸素ガス生成時に生じる圧力変動の振幅は呼吸圧力の振幅に比べて大きく、また呼吸圧力の振幅も酸素流路を通ることによる圧損で低下するため、呼吸数計測装置4で計測される圧力から直接患者の呼吸圧力波形を計測することは困難である。
 本実施形態においては、呼吸数計測装置4は、圧力センサ6に接続される、演算部及び推定部を有するマイコン部7を有する。圧力センサ6は好ましくは微差圧センサである。さらに、好ましくは、呼吸数計測装置4は、計測された呼吸数を表示する、例えば液晶ディスプレイである表示部8を有する。表示部8は、マイコン部7に接続され、マイコン部7によって制御される。
 本実施形態においては、呼吸数計測装置4は、配管内の圧力を検知して出力する検知部として、例えば、圧力センサ6、好ましくは微差圧センサ6と、検知部に電気的に接続されるマイコン部7を有する。呼吸数計測装置4は、微差圧センサ6に接続される容積部と、配管と容積部の間を接続するオリフィス5から構成される圧力平滑化部を有してもよい。圧力平滑化部を設ける理由は以下の通りである。
 患者の呼吸圧は通常±10~100Pa程度であるため、呼吸数計測装置4にて呼吸圧を取得するためには微差圧センサ6はレンジ±100Pa程度のセンサを用いるのが好ましい。酸素濃縮装置から酸素が供給されている状況では、酸素供給による供給圧が常に生じており、例えば、酸素供給による供給圧は、1LPM(litter per minute:リットル/分)でも300Pa程度存在する。そのため微差圧センサの他端を大気解放した状態で、上述のように微差圧センサ6の一端を酸素供給経路に接続すると、微差圧センサ6の測定レンジをオーバーしてしまう。そのため、オリフィス5経由後の圧を微差圧センサ6の他端に印加することで、患者の呼吸情報を含んだ圧を微差圧センサ6の測定レンジ範囲内で取得することが好ましい。
 微差圧センサ6の他端には微差圧センサ6の測定レンジをオーバーしない程度の圧力が印加されていればよく、その方法は本実施形態の例に限定されない。微差圧センサ6の測定レンジが酸素供給による供給圧よりも大きく、かつ患者の呼吸に伴う圧力変動を検出できるレベルの分解能があれば、微差圧センサ6の他端は大気解放でもよい。微差圧センサ6の代わりにゲージ圧や絶対圧を測定する圧力センサを用いてもよい。
 呼吸に伴う圧力変化は、配管に流れる酸素濃縮ガスの微小な流量変動としても現れるため、呼吸数計測装置4は、圧力センサ6の代わりに流量センサを用いてもよい。呼吸数計測装置4は、圧力センサ6及び流量センサの両方を有してもよい。圧力センサ6及び/又は流量センサは、配管2内の患者の呼吸情報を含んだ配管内の管内圧力及び/又は管内気体流量を検知して、圧力データ及び/又は管内気体流量データを出力する検知部となる。
[マイコン部]
 図3は、マイコン部7の構成ブロックの一例を示す図である。
 圧力センサ6、好ましくは微差圧センサ6には、演算部722および推定部723を有するマイコン部7が接続されている。演算部722および推定部723は、検知部である微差圧センサ6で検知された患者の呼吸情報を含んだ配管内の管内圧力データ及び又は管内気体流量データを受け取って、後述の処理を実行することで患者の呼吸情報を取得する。
 マイコン部7は、酸素濃縮装置の酸素生成機能や表示・ユーザインタフェース機能などを処理する処理部と同一のマイコンであってもよいし、分離されていてもよい。分離されている場合は、酸素生成機能を処理するマイコンから、PSAの周期Tの切り替わりのタイミングを入手して、演算に利用する。
 マイコン部7は、記憶部71と処理部72を有する。記憶部71は、1又は複数の半導体メモリにより構成される。例えば、RAMや、フラッシュメモリ、EPROM、EEPROM等の不揮発性メモリの少なくとも一つを有する。記憶部71は、処理部72による処理に用いられるドライバプログラム、オペレーティングシステムプログラム、アプリケーションプログラム、データ等を記憶する。
 例えば、記憶部71は、ドライバプログラムとして、検知部である微差圧センサ6等を制御するデバイスドライバプログラムを記憶する。コンピュータプログラムは、例えば、CD-ROM、DVD-ROM等のコンピュータ読み取り可能な可搬型記録媒体から、公知のセットアッププログラム等を用いて記憶部71にインストールされてもよい。また、プログラムサーバ等からダウンロードしてインストールしてもよい。
 更に、記憶部71は、所定の処理に係る一時的なデータを一時的に記憶してもよい。記憶部71は、呼吸数の推定に利用するための閾値711、変動値データファイル712等を記憶する。
 処理部72は、一又は複数個のプロセッサ及びその周辺回路を有する。処理部72は、呼吸数計測装置4の全体的な動作を統括的に制御するものであり、例えば、MCU(Micro Control Unit)である。
 処理部72は、記憶部71に記憶されているプログラム(オペレーティングシステムプログラム、ドライバプログラム、アプリケーションプログラム等)に基づいて処理を実行する。処理部72は、複数のプログラム(アプリケーションプログラム等)を並列に実行してもよい。処理部72は、検知データ取得部721、演算部722、推定部723、呼吸数出力部724等を有する。
 処理部72が有するこれらの各部は、独立した集積回路、回路モジュール、マイクロプロセッサ、又はファームウェアとしてマイコン部7に実装されてもよい。
[実施形態で行われる演算処理の原理]
 図2に示す構成を用いて鼻カニューラ3から患者呼吸モデルの呼吸圧を印加し、呼吸情報を取得したデータを用いて、本実施形態で行われる処理の原理について説明する。
 連続流5LPM、呼吸数計測装置4の下流側に延長チューブを20m接続した時に取得しデータ群を図4、図5、図6に示す。連続流とは、濃縮酸素ガス供給方式のひとつで、連続的に一定流量の濃縮酸素ガスを供給する方式である。
 図4は、連続流5LPM、延長チューブを20m接続したときの患者の呼吸情報とPSA式酸素濃縮装置1のPSA圧データとを含んだ圧データを示す図である。
 図5は、連続流5LPM、延長チューブを20m接続したときの患者呼吸が無い場合、あるいは呼吸成分を除去した後のPSA式酸素濃縮装置1のPSA圧データを示す図である。
 PSA式酸素濃縮装置によるPSA圧とは、上述したPSA式酸素濃縮装置1により酸素を生成する際に発生するPSA式酸素濃縮装置1の吸着筒周期に伴う周期的な圧変化である。PSA圧の波形の周期はPSA式酸素濃縮装置1の吸着筒切替周期と一致している。
 図6は、図4の成分から図5の成分をソフトウェアにて差分処理をした結果を示す図である。差分処理とは、任意時刻における両データの差を計算する処理である。差分処理により、呼吸数計測装置4は、PSA圧成分を除去し、患者呼吸モデルの呼吸圧を検出することができている。ソフトによる差分処理を実施することで、呼吸数計測装置4は、PSA式酸素濃縮装置1から鼻カニューラと延長チューブ20mを介して酸素を連続的に5LPM吸入している状況下においても患者の呼吸情報が取得可能である。
 PSA圧成分等の変動値データを除去し、患者呼吸情報データを抽出する方法は、特許文献6に記載される、酸素濃縮装置自身の圧力変動を予め計測し、記憶しておいて、検出した圧力波形から差し引く方法によってもよいし、リアルタイムで計測したものを差し引く方法によってもよい。
 PSA式酸素濃縮装置が、少なくとも総運転時間の一部の期間、単一の周期Tで加圧・減圧を切り替える場合、演算部722はある時間tに於ける圧力は流量値Y(t)と、Y(t)及びY(t-T)、Y(t-2T)、・・・、Y(t-nT)(nは予め決められた任意の整数)を平均化処理した値X(t)との差分を算出して呼吸情報を抽出する。PSA圧成分が除去された元の呼吸波形がよく再現できる。
 図7は、本実施形態の演算処理の概念図である。(I)は呼吸波形、(II)はPSA波形のモデル波形である。TはPSA波形の周期であり、吸着筒の切替周期に一致する。(III)は、(I)と(II)を加算した波形であり、配管部で測定される圧力に相当する。(IV)は、(III)を周期Tごとに分割し重ね合わせたものである。呼吸の周期とPSAの周期が完全に一致しない限り、周期Tで切り取った波形に於いては、呼吸波形はランダムに表れる。(V)は(IV)で重ね合わせた波形を平均化処理したものである。測定される圧力データ及び/又は気体流量データ(III)に基づいて、PSA式酸素濃縮装置の動作により発生する酸素濃縮ガスの周期的な圧力変化及び/又は流量変化を示す変動値データX(t)が、重ね合わせて平均することにより推定される(V)。
 一つの例では、Tを1周期として5周期分の単純移動平均が用いられている。(VI)は(III)の波形の最後のT部分から(V)を差し引いた波形である。元の呼吸波形(I)が精度よく再現できていることが分かる。Tの選び方は、複数の吸着筒の内1つの吸着筒がコンプレッサと連通している状態が始まってから、別の吸着筒がコンプレッサと連通する状態に切り替わるまでの時間をTとしてもよい。Tの選び方は、1つの吸着筒がコンプレッサと連通する状態が始まってから、他の吸着筒がコンプレッサと連通する状態を経て、再び最初の吸着筒がコンプレッサに連通する状態に切り替わるまでをTとしてもよい。さらには、別の例として、Tの選び方は、それぞれの切替時間の倍数をTとしてもよい。
 実施例では、酸素濃縮ガスの周期的な圧力変化及び/又は流量変化と推定される、重ね合わせて平均された変動値データX(t)を示す(V)の計算をするに当たって、5周期分の移動平均が用いられているが、平均の仕方は別の方法でもよい。一般的には、(V)のステップに於いて計算する平均値X(t)は以下の数式で表される。本願全体を通じて、iは虚数を示すものではなく、単に整変数iを示すものであるとする。
Figure JPOXMLDOC01-appb-M000001
 X(t)は時間tに於ける平均化処理後の値であり、Y(t)は時間tに於ける圧力の実測値である。また、aiは加重平均に於ける重み付けの係数であり任意の実数を選択することができる。nは平均を取るYの個数を示す。例えばn=4とし、a0~a4を1とすれば、図6の計算と同じになる。また、ai=e-bi(bは任意の正の実数)とし、n=∞とすれば指数平滑平均とすることができる。
 実際の計算では無限大を扱うことは難しいため、X(t)は、X(t)=(1-e-b)Y(t)+e-bX(t-T)等のように逐次計算により漸近値を得る。n及びaiを適切に選ぶことによって、FIR(Finite Impulse Response)フィルタのような、より収束が早い平均化処理も可能である。
 式(1)は一見、単にY(t)の時間平均を数値計算的に求める式に見える。通常の時間平均ではTとして、Y(t)の変動周期に対して十分小さい値を選択するのに対し、本実施形態では、TとしてPSAプロセスの吸脱着周期を基本に選択することが重要である。これは、単純に測定値Y(t)を時間的に平均化処理するのではなく、周期Tに亘るY(t)の波形を一つの単位として、Tの整数倍の時間だけ遡った波形の平均値をとる事を意味する。波形の平均値を取ることにより、Tの間隔で周期的に表れる特性をもつPSA圧を精度よく推定することが可能となる。
 本例は、呼吸数計測部として、圧力を測定する方法を示したが、実際には圧力の変動に伴い、配管を流れる流量も変化するため、圧力の代わりに流量を測定する方法も用いることができる。また、装置の運転条件や環境条件によって圧力値と流量値を組み合わせて、あるいは切り替えて測定する方法も用いることが出来る。
 演算処理により取得した呼吸情報は、圧力や流量などの呼吸のリアルタイムの情報を波形として表した、いわば生データである。このデータをそのまま記録、あるいは送信してもよいが、データ量が膨大となるうえ、解析にも時間がかかる。そのため、呼吸数計測装置4内部で、呼吸数データを自動的に算出した上で、記録、及び/又は送信することが望ましい。
 例に示したように、検出できる呼吸波形は、気流による圧力変動やカニューラの揺れなどに伴うノイズ成分が多い。圧力が下降から上昇に転じるタイミングを算出するピーク検出や、閾値以上の値になったタイミングを検出する方法では、適切に呼吸周期、呼吸数を算出することができない。
 例えば、図6に於いては、カウントすべき呼吸波形は、図中a1、a2、a3であらわされるピーク及びその直後に現れる谷として現れており、それ以外の部分はノイズ成分である。圧力が上昇から下降に転じるタイミングとしてこのピークを検出しようとしても、図中のb部分を検出してしまう可能性が高い。
[呼吸数推定処理]
 ある一定の閾値以上をピークとして判断する方法であっても、閾値をXとして設定した場合は、a2のような低いピークが検出できなくなる可能性があり、Yに設定した場合は、ノイズ成分bをピークとして検出してしまう可能性がある。本例では3呼吸分程度しか波形を示していないが、実際の波形ではここで示した以上にピークやノイズ成分の高さがばらつく可能性があり、ピークを確実にとらえ、ノイズを確実に除去できる閾値の設定は極めて困難である。
 発明者は鋭意検討した結果、上記の方法ではなく、元の波形と、そこから時間Δtだけずらした波形との自己相関係数を求め、Δtを変化させて自己相関係数がピークとなるΔtを求めることで、呼吸数を高い検出力で検出できる事を見出した。発明者は患者呼吸情報から所定時間当たりの呼吸数を推定する呼吸数推定方法を見出した。
 自己相関係数Rは、以下のような計算式にて算出可能である。
Figure JPOXMLDOC01-appb-M000002
 ここで、Δtは時間のずらし量、t0はデータの取得間隔、nは1回の自己相関係数の計算に用いるデータ数、f(t)は、時刻tに於ける患者呼吸情報データの値で、計算に於いては、t0の間隔でn個取得する。μ、σは、f(t)の平均値及び標準偏差であるが、実際の計算上は、取得したn個のf(t)の値の平均、標準偏差を用いてもよい。また、本実施形態の自己相関係数の算出方法は、一例であり、他の自己相関係数の算出方法でもよい。
 本実施形態で算出する自己相関係数Rには、正規化係数
Figure JPOXMLDOC01-appb-M000003
を乗じることにより、値を-1.0~+1.0の範囲に規定している。
 図8は、正規化係数を乗じた上式を用いて、図6で示した算出した波形と、時間Δtだけ時間軸に沿って原波形をずらした波形との自己相関係数の計算結果をΔtに対してプロットして作成したグラフである。自己相関係数は、Δt=0から、t0ずつ増加させて都度計算する。Δt=0は原波形同士の相関となるため1となる。その後、高い相関を持つ点として点P1及び点P2が求められる。
 図8に示されるように、生の患者呼吸情報データでは閾値やピーク検出によって呼吸周期を判断することが難しかった波形が、自己相関係数を求めることにより容易に呼吸周期を求めることが出来ることが分かる。
 自己相関係数の計算に用いるf(t)のデータ区間は算出すべき呼吸周期の範囲から定められる。発明者が鋭意検討した結果、f(t)のデータ区間は算出すべき最低の呼吸周期に対して2倍以上が必要である事が分かった。一般的な成人の呼吸数は1分間に12~20回程度であり、これは呼吸周期3~5秒程度に相当する。f(t)のデータ区間は5秒×2=10秒は必要であることが分かった。
 自己相関係数のピークを判定するための自己相関係数の閾値としては、少なくとも0.3以上、好ましくは0.3以上0.7以下の間に設定する事が必要であることが判明した。それ以下の閾値ではノイズやベースラインの乱れなどによる偶発的な自己相関値の上昇をピークと誤判定する可能性が増加し、それ以上の閾値では、呼吸周期による自己相関値の上昇を見逃す可能性が高いことが分かった。
 呼吸数は、自己相関係数が一定の値以上かつピーク値を取ったときのΔtの値の呼吸間隔の逆数により算出される。図8において最初のピークの点P1に於けるΔtが呼吸間隔となり、1分間を呼吸間隔で除すると1分当たりの呼吸数となる。1分当たりの呼吸数(bpm)は、以下の式で導かれる。
1分当たりの呼吸数(bpm)=60秒÷(ピーク点のステップ数×サンプリング時間)
 呼吸波形によっては、図8の様に一定の値以上のピークが複数点存在する事もあるが、右側のピークP2はΔtを複数呼吸分ずらした時に現れる倍周期のピークである。呼吸数の推定では呼吸の基本周期に該当すると推測される、最も左、すなわちΔtが最も小さいピークP1を用いるべきである。
 本実施形態に示した計算により、濃縮器に接続された呼吸数計測部にて取得した呼吸波形から、呼吸数の情報が自動的に精度よく算出されることができる。
 呼気吸気に基づく、例えば圧力に関連する情報の変化を表す患者呼吸情報データf(t)は、呼吸以外の小さなノイズしか入っていない場合、偶然に一定周期で自己相関係数が高くなってしまうことがある。そこで、患者呼吸情報データf(t)の分散による判定を加えることにより、呼吸が入っているかもしれないと判断できるレベルの大きさの波がないときは、推定部723は、呼吸数の算定をしないようしてより精度を高めることができる。呼吸の有無に判定は、患者呼吸情報データf(t)の分散による判定に検定されず、他の判定方法を使用してもよい。
 患者呼吸情報データf(t)の分散σ2は次式で計算される。
Figure JPOXMLDOC01-appb-M000004
 ここで、nは1回の自己相関係数の計算に用いるデータ数、f(t)は、時刻tに於ける患者呼吸情報データの値。μ、σ2は、f(t)の平均値及び分散である。
 一例として、PSA式酸素濃縮装置の場合の分散判定の閾値は、
呼吸数8未満の場合 :σ2≧18(Pa2
呼吸数8~10の場合:σ2≧18(Pa2
呼吸数11以上の場合:σ2≧7(Pa2
であるが、使用する酸素供給装置によりに値が異なる可能性がある。呼吸数等のパラメーターによりばらつきがあるので、分散判定の閾値は使用する酸素供給装置により適切に設定する必要がある。
 一例として、呼吸数8未満の場合の患者呼吸情報データf(t)の分散σ2の第1閾値THD1はTHD1=18とする。呼吸数8~10の場合の患者呼吸情報データf(t)の分散σ2の第2閾値THD2はTHD2=18とする。呼吸数11以上の場合の患者呼吸情報データf(t)の分散σ2の第3閾値THD3はTHD3=7とする。患者呼吸情報データf(t)が分散σ2の閾値以上の場合、呼吸の有無の判定が行われる。
 図9は、患者呼吸情報データを抽出する処理の一例を示すフローチャート図である。
 図9に示す患者呼吸情報データを抽出する処理は、記憶部71に予め記憶されたコンピュータプログラムに従ってマイコン部7にて実行される。検知データ取得部721は、検知部である圧力センサ6から、圧力データY(t)を取得する(ST101)。演算部722は、Y(t)のn周期分の平均値X(t)を算出する(ST102)。演算部722は、Y(t)とX(t)の差分をとることにより患者呼吸情報データf(t)を抽出する(ST103)。予め計測記憶したPSA圧力変動を用いる場合は、演算部722は、ST102の処理を行わず、ST101の処理の後、ST103として、Y(t)と予め計測記憶したPSA圧力変動の差分をとることにより患者呼吸情報データf(t)を抽出する。
 呼吸数計測装置4が稼働中、呼吸数を推定する処理は繰り返し行われ、例えば圧力Y(t)の測定データがt0の間隔でn個取得されるごとに、患者呼吸情報データf(t)が抽出され、更新される。患者呼吸情報データf(t)は、所定の間隔で抽出され、更新されてもよい。
 図10は、患者呼吸情報データに基づいて、呼吸数を推定する処理の一例を示すフローチャート図である。
 図10に示す呼吸数を推定する処理は、記憶部71に予め記憶されたコンピュータプログラムに従ってマイコン部7にて実行される。推定部723は、患者呼吸情報データf(t)の平均値μと標準偏差σを算出する(ST201)。まず、呼吸数の有無判定を行うため、推定部723は、標準偏差σの平方数である分散σ2が所定の閾値THD以上であるか否かを判断する(ST202)。分散σ2が閾値THD未満のときは(ST202:NO)、呼吸数出力部724は、計算不能であった旨を出力して(ST213)、処理を終了する。分散σ2が閾値THD以上のときは(ST202:YES)、推定部723は、時間のずらし量として時間Δtをゼロ(0)に設定する(ST203)。
 推定部723は、時間Δtにデータの取得間隔t0をインクリメントする(ST204)。推定部723は、患者呼吸情報データf(t)の時間Δtにおける自己相関係数R(Δt)を算出する(ST205)。推定部723は、算出した自己相関係数R(Δt)を記憶部71に書き出す(ST206)。推定部723は、時間Δtがnt0に達したかを判断する(ST207)。時間Δtがnt0に達していないときはST203に戻り、推定部723は処理を繰り返す(ST207:NO)。
 時間Δtがnt0に達したときは(ST207:YES)、推定部723は、記憶部71に記憶されている自己相関係数R(Δt)を読み出す(ST208)。推定部723は、読み出した自己相関係数R(Δt)を比較し、最大となる自己相関係数Max(R(Δt))、すなわちピークとなる自己相関係数R(Δt)があるか否かを判断する(ST209)。具体的には、Max(R(Δt))が所定の閾値THC以上であるか否かが判断される。
 最大となる自己相関係数Max(R(Δt))が所定の閾値THC以上であるとき(ST209:YES)、推定部723は、最大となる自己相関係数Max(R(Δt))の時間Δtを呼吸間隔に設定する(ST210)。更に、推定部723は、呼吸間隔Δtから呼吸数を推定する(ST211)。呼吸数出力部724は、呼吸数信号を出力し(ST212)、処理は終了する。例えば、呼吸数信号が入力された表示部8は呼吸数を表示する。
 最大となる自己相関係数Max(R(Δt))が所定の閾値THC未満のときは(ST209:NO)、呼吸数出力部724は、計算不能であった旨を出力して(ST213)、処理を終了する。呼吸数計測装置4が稼働中、呼吸数を推定する処理は繰り返し行われ、呼吸数計測装置4は、呼吸数を更新し、表示部8に表示することができる。表示部8への表示は、所定の間隔で行われてもよい。
 当業者は、本発明の精神及び範囲から外れることなく、様々な変更、置換、及び修正をこれに加えることが可能であることを理解されたい。
 1  PSA式酸素濃縮装置
 2  配管
 3  鼻カニューラ
 4  呼吸数計測装置
 5  オリフィス
 6  圧力センサ
 7  マイコン部
 71  記憶部
 72  処理部
 721  検知データ取得部
 722  演算部
 723  推定部
 724  呼吸数出力部
 8  表示部

Claims (7)

  1.  患者と接続され、且つ、加圧・減圧を周期的に繰り返すことで空気中の酸素を濃縮する圧力変動吸着方式の酸素濃縮装置から酸素濃縮ガスを患者に供給する配管内の管内圧力及び/又は管内気体流量を検知して、圧力データ及び/又は気体流量データを出力する検知部と、
     前記圧力データ及び/又は前記気体流量データに基づいて、患者呼吸情報データを抽出する演算部と、
     前記患者呼吸情報データに基づいて所定時間当たりの呼吸数を推定する推定部と、を有し、
     前記推定部は、時間Δtを変化させながら、前記患者呼吸情報データと、前記患者呼吸情報データから前記時間Δtだけずらしたデータとの自己相関係数を求め、前記自己相関係数がピークとなる時間Δtを呼吸間隔として、前記呼吸数を推定する、
    ことを特徴とする呼吸数計測装置。
  2.  前記推定部は、少なくとも10秒以上の前記患者呼吸情報データを示す波形を用いて、前記呼吸数を推定する、請求項1に記載の呼吸数計測装置。
  3.  前記推定部は、所定の閾値を超えた前記自己相関係数のピークを前記呼吸数の推定に利用する、請求項1又は2に記載の呼吸数計測装置。
  4.  前記閾値は、0.3以上且つ0.7以下の値である、請求項3に記載の呼吸数計測装置。
  5.  前記所定の閾値を超えた前記自己相関係数のピークを取得できない場合、計算不能であった旨を出力する、請求項1~4の何れか一項に記載の呼吸数計測装置。
  6.  前記演算部は、前記圧力データ及び/又は前記気体流量データに基づいて、酸素濃縮装置の動作により発生する前記酸素濃縮ガスの周期的な圧力変化及び/又は流量変化を示す変動値データを推定し、前記変動値データを前記圧力データ及び/又は前記気体流量データから除去することによって、前記患者呼吸情報データを抽出する、請求項1~5の何れか一項に記載の呼吸数計測装置。
  7.  予め測定した前記変動値データを記憶する記憶部を更に有し、
     前記演算部は、前記変動値データを前記圧力データ及び/又は前記気体流量データから除去することによって、前記患者呼吸情報データを抽出する、請求項1~5の何れか一項に記載の呼吸数計測装置。
PCT/JP2019/038359 2018-09-28 2019-09-27 呼吸数計測装置 WO2020067505A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
EP19867937.5A EP3858409B1 (en) 2018-09-28 2019-09-27 Respiration rate measurement device
US17/279,215 US20220031986A1 (en) 2018-09-28 2019-09-27 Respiratory rate measurement device
CN201980063693.XA CN112739403B (zh) 2018-09-28 2019-09-27 呼吸次数计测装置
JP2020549479A JP7055217B2 (ja) 2018-09-28 2019-09-27 呼吸数計測装置
KR1020217009195A KR102526829B1 (ko) 2018-09-28 2019-09-27 호흡수 계측 장치

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2018-185785 2018-09-28
JP2018185785 2018-09-28

Publications (1)

Publication Number Publication Date
WO2020067505A1 true WO2020067505A1 (ja) 2020-04-02

Family

ID=69953176

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/038359 WO2020067505A1 (ja) 2018-09-28 2019-09-27 呼吸数計測装置

Country Status (6)

Country Link
US (1) US20220031986A1 (ja)
EP (1) EP3858409B1 (ja)
JP (1) JP7055217B2 (ja)
KR (1) KR102526829B1 (ja)
CN (1) CN112739403B (ja)
WO (1) WO2020067505A1 (ja)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06190045A (ja) 1992-10-13 1994-07-12 Teijin Ltd 呼吸測定器、呼吸用気体供給装置及び呼吸モニタリングシステム
JPH0796035A (ja) 1993-09-30 1995-04-11 Teijin Ltd 呼吸数測定装置、呼吸用気体供給装置及び呼吸モニタリングシステム
JP2001266205A (ja) * 2000-03-22 2001-09-28 Toshiba Corp 磁気情報検出装置
JP2001286566A (ja) 2000-02-14 2001-10-16 Fukuda Sangyo:Kk 医療用の濃縮酸素ガス供給装置
JP2002102187A (ja) * 2000-09-29 2002-04-09 Matsushita Electric Ind Co Ltd 生体検出装置
JP2002336207A (ja) * 2001-05-14 2002-11-26 Matsushita Electric Ind Co Ltd 在床異常モニタ装置
JP2011518016A (ja) 2008-04-21 2011-06-23 レール・リキード−ソシエテ・アノニム・プール・レテュード・エ・レクスプロワタシオン・デ・プロセデ・ジョルジュ・クロード 酸素療法治療の患者による遵守を検出するためのデバイス
JP2014068782A (ja) * 2012-09-28 2014-04-21 Teijin Pharma Ltd 酸素供給装置
JP2015085191A (ja) 2013-09-27 2015-05-07 公立大学法人奈良県立医科大学 呼吸監視装置及び呼吸監視方法
WO2018180392A1 (ja) 2017-03-31 2018-10-04 帝人ファーマ株式会社 呼吸情報取得装置および呼吸情報取得方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2004266540B9 (en) * 2003-08-26 2010-04-01 Teijin Limited Oxygen-concentrating device
US7011693B2 (en) * 2003-11-12 2006-03-14 General Motors Corporation Control of a hydrogen purifying pressure swing adsorption unit in fuel processor module for hydrogen generation
ES2823307T3 (es) * 2009-08-13 2021-05-06 Hidetsugu Asanoi Dispositivo para calcular información de forma de onda respiratoria y dispositivo médico que utiliza la información de forma de onda respiratoria
US10255647B2 (en) * 2009-09-28 2019-04-09 Caire Inc. Controlling and communicating with respiratory care devices
NZ625053A (en) * 2009-11-16 2015-11-27 Resmed Ltd Methods and apparatus for adaptable pressure treatment of sleep disordered breathing
WO2014059405A1 (en) * 2012-10-12 2014-04-17 Inova Labs, Inc. Method and systems for the delivery of oxygen enriched gas
GB2519964A (en) * 2013-11-01 2015-05-13 Linde Aktiengesellshcaft Apparatus and method for detecting health deterioration
CA2932826A1 (en) * 2014-02-20 2015-08-27 Covidien Lp Systems and methods for filtering autocorrelation peaks and detecting harmonics
US11071841B2 (en) * 2015-06-22 2021-07-27 Caire Inc. Wearable oxygen generator and docking station enabling higher oxygen flow capacity
KR102426719B1 (ko) * 2016-05-17 2022-07-28 피셔 앤 페이켈 핼스케어 리미티드 유동 치료 장치를 위한 유로 감지
JP2018180392A (ja) 2017-04-18 2018-11-15 株式会社東芝 画像形成装置

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06190045A (ja) 1992-10-13 1994-07-12 Teijin Ltd 呼吸測定器、呼吸用気体供給装置及び呼吸モニタリングシステム
JPH0796035A (ja) 1993-09-30 1995-04-11 Teijin Ltd 呼吸数測定装置、呼吸用気体供給装置及び呼吸モニタリングシステム
JP2001286566A (ja) 2000-02-14 2001-10-16 Fukuda Sangyo:Kk 医療用の濃縮酸素ガス供給装置
JP2001266205A (ja) * 2000-03-22 2001-09-28 Toshiba Corp 磁気情報検出装置
JP2002102187A (ja) * 2000-09-29 2002-04-09 Matsushita Electric Ind Co Ltd 生体検出装置
JP2002336207A (ja) * 2001-05-14 2002-11-26 Matsushita Electric Ind Co Ltd 在床異常モニタ装置
JP2011518016A (ja) 2008-04-21 2011-06-23 レール・リキード−ソシエテ・アノニム・プール・レテュード・エ・レクスプロワタシオン・デ・プロセデ・ジョルジュ・クロード 酸素療法治療の患者による遵守を検出するためのデバイス
JP2014068782A (ja) * 2012-09-28 2014-04-21 Teijin Pharma Ltd 酸素供給装置
JP2015085191A (ja) 2013-09-27 2015-05-07 公立大学法人奈良県立医科大学 呼吸監視装置及び呼吸監視方法
WO2018180392A1 (ja) 2017-03-31 2018-10-04 帝人ファーマ株式会社 呼吸情報取得装置および呼吸情報取得方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3858409A4

Also Published As

Publication number Publication date
KR102526829B1 (ko) 2023-04-27
EP3858409B1 (en) 2024-05-08
EP3858409A4 (en) 2022-06-29
US20220031986A1 (en) 2022-02-03
CN112739403A (zh) 2021-04-30
EP3858409A1 (en) 2021-08-04
CN112739403B (zh) 2024-04-30
JPWO2020067505A1 (ja) 2021-05-20
JP7055217B2 (ja) 2022-04-15
KR20210047936A (ko) 2021-04-30

Similar Documents

Publication Publication Date Title
US5876352A (en) Process for determining the mechanical properties of the respiratory system of a respirated patient and device for carrying out the process
CA2788781C (en) Event-based delay detection and control of networked systems in medical ventilation
JP2011522621A (ja) 患者の努力に比例した換気のためのシステムおよび方法
US20190192796A1 (en) Ventilator with error detection for flow sensors
WO2012059768A1 (en) Apparatus and methods for breath sampling
WO2020067526A1 (ja) 呼吸モニタリング装置
JP6808020B2 (ja) 呼吸情報取得装置および呼吸情報取得方法
JP6432970B2 (ja) 呼吸監視装置及び呼吸監視方法
WO2020067505A1 (ja) 呼吸数計測装置
JP6952909B2 (ja) 呼吸情報取得装置
JP6893432B2 (ja) 呼吸測定装置
CN103513022A (zh) 一种不需要严格控制呼气流量的一氧化氮测量方法和设备
JPH065645U (ja) 人工呼吸装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19867937

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2020549479

Country of ref document: JP

Kind code of ref document: A

ENP Entry into the national phase

Ref document number: 20217009195

Country of ref document: KR

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2019867937

Country of ref document: EP

Effective date: 20210428