RU2721878C2 - Способы и системы для анализа данных о содержании глюкозы, измеренных у больного диабетом - Google Patents

Способы и системы для анализа данных о содержании глюкозы, измеренных у больного диабетом Download PDF

Info

Publication number
RU2721878C2
RU2721878C2 RU2017138101A RU2017138101A RU2721878C2 RU 2721878 C2 RU2721878 C2 RU 2721878C2 RU 2017138101 A RU2017138101 A RU 2017138101A RU 2017138101 A RU2017138101 A RU 2017138101A RU 2721878 C2 RU2721878 C2 RU 2721878C2
Authority
RU
Russia
Prior art keywords
glucose
insulin
bolus
measured
threshold
Prior art date
Application number
RU2017138101A
Other languages
English (en)
Other versions
RU2017138101A3 (ru
RU2017138101A (ru
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 Ф. Хоффманн-Ля Рош Аг
Publication of RU2017138101A publication Critical patent/RU2017138101A/ru
Publication of RU2017138101A3 publication Critical patent/RU2017138101A3/ru
Application granted granted Critical
Publication of RU2721878C2 publication Critical patent/RU2721878C2/ru

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
    • A61B5/0015—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by features of the telemetry system
    • A61B5/0022—Monitoring a patient using a global network, e.g. telephone networks, internet
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/14532—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring glucose, e.g. by tissue impedance measurement
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4836—Diagnosis combined with treatment in closed-loop systems or methods
    • A61B5/4839—Diagnosis combined with treatment in closed-loop systems or methods combined with drug delivery
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7282—Event detection, e.g. detecting unique waveforms indicative of a medical condition
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/7405—Details of notification to user or communication with user or patient ; user input means using sound
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/7455—Details of notification to user or communication with user or patient ; user input means characterised by tactile indication, e.g. vibration or electrical stimulation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/746—Alarms related to a physiological condition, e.g. details of setting alarm thresholds or avoiding false alarms
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/7475—User input or interface means, e.g. keyboard, pointing device, joystick
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/483—Physical analysis of biological material
    • G01N33/487—Physical analysis of biological material of liquid biological material
    • G01N33/49—Blood
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Chemical & Material Sciences (AREA)
  • Optics & Photonics (AREA)
  • Medicinal Chemistry (AREA)
  • Pharmacology & Pharmacy (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Emergency Medicine (AREA)
  • Power Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Business, Economics & Management (AREA)
  • General Business, Economics & Management (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Hematology (AREA)
  • Ecology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Food Science & Technology (AREA)
  • Urology & Nephrology (AREA)

Abstract

Группа изобретений относится к медицине, а именно к анализу уровня глюкозы у больного диабетом. Предложена система для реализации способа, причем система содержит контроллер, содержащий машиночитаемый носитель информации с хранящимися в нем машиноисполняемыми командами, при выполнении которых в контроллере контроллером осуществляется: получение значений глюкозы от сенсора глюкозы, подключенного к больному диабетом, максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента, анализ измеренных значений глюкозы с использованием средства вероятностного анализа для определения граничного значения глюкозы (gρ) при пороге вероятности того, что больному диабетом потребуется заданная доза инсулина, определение порога глюкозы
Figure 00000076
с использованием максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента, и сравнение граничного значения глюкозы (gρ) с порогом глюкозы
Figure 00000076
, причем если граничное значение глюкозы (gρ) выше порога глюкозы

Description

Настоящее изобретение в целом относится к обработке данных о содержании глюкозы, измеряемом у больного диабетом, в частности для отправки больному диабетом предупреждения, если вероятность возникновения потребности в инсулине превысит пороговый уровень.
Как известно, люди страдают диабетом первого или второго типа, когда их организм неспособен должным образом регулировать уровень сахара в крови. Многие из этих людей могут использовать непрерывное мониторирование гликемии (НМГ), таким образом контролируя свой уровень глюкозы на постоянной основе. Для осуществления НМГ человеку можно ввести под кожу сенсор глюкозы, способный измерять уровень глюкозы в его интерстициальной жидкости. Сенсор глюкозы может периодически измерять у человека уровень глюкозы с известным временным интервалом, например, каждую минуту, и передавать результаты измерения глюкозы в инсулиновую помпу, глюкометр, смартфон или другой электронный прибор контроля, также называемый монитором.
Объектами изобретения являются, по меньшей мере, способ и система для анализа уровня глюкозы у больного диабетом.
По меньшей мере в одном варианте осуществления изобретения предлагаемый способ анализа уровня глюкозы у больного диабетом включает:
- получение вычислительным устройством измеренных значений глюкозы от системы непрерывного мониторирования гликемии, подключенной к больному диабетом, максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента;
- анализ измеренных значений глюкозы, выполняемый в вычислительном устройстве с использованием средства вероятностного анализа для определения граничного значения глюкозы (gρ) при пороге вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина;
- выполняемое вычислительным устройством определение порога глюкозы
Figure 00000001
с использованием максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента; и
- выполняемое вычислительным устройством сравнение граничного значения глюкозы (gρ) с порогом глюкозы
Figure 00000001
, причем если граничное значение глюкозы (gρ) превышает порог глюкозы
Figure 00000001
, то вычислительное устройство выдает через пользовательский интерфейс предупреждение.
По меньшей мере в одном варианте осуществления изобретения предлагаемая система для анализа уровня глюкозы у больного диабетом включает в себя:
Система для анализа уровня глюкозы у больного диабетом, включающая в себя:
- контроллер, содержащий машиночитаемый носитель информации с хранящимися в нем машиноисполняемыми командами, при выполнении которых в контроллере контроллером осуществляется:
- получение вычислительным устройством измеренных значений глюкозы от сенсора глюкозы, подключенного к больному диабетом, максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента,
- анализ измеренных значений глюкозы, выполняемый в вычислительном устройстве с использованием средства вероятностного анализа для определения граничного значения глюкозы (gρ) при пороге вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина,
- выполняемое вычислительным устройством определение порога глюкозы
Figure 00000001
с использованием максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента, и
- выполняемое вычислительным устройством сравнение граничного значения глюкозы (gρ) с порогом глюкозы
Figure 00000001
, причем если граничное значение глюкозы (gρ) выше порога глюкозы
Figure 00000001
, то вычислительное устройство выдает предупреждение; и
- связанный с контроллером калькулятор болюсов, способный определять болюс инсулина на основании измеренных значений глюкозы.
По меньшей мере в одном варианте осуществления настоящего изобретения предлагаемый способ включает также определение максимального безопасного болюса инсулина. При необходимости предлагаемый способ может также включать выдачу через пользовательский интерфейс тревожного сообщения, если введенный вручную запрос болюса превышает максимальный безопасный болюс инсулина.
По меньшей мере в одном варианте осуществления настоящего изобретения граничное значение глюкозы (gρ) определяют при помощи уравнения:
Figure 00000002
где:
Figure 00000003
- вариация при измерении глюкозы,
ρ - порог вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина, и
g - измеренное значение глюкозы.
По меньшей мере в одном варианте осуществления настоящего изобретения порог глюкозы
Figure 00000001
определяют при помощи уравнения:
Figure 00000004
где:
gm - максимально допустимая глюкоза,
с - поправочный коэффициент,
I0 - минимальная доза инсулина.
По меньшей мере в одном варианте осуществления настоящего изобретения заданная доза инсулина равна одной единице.
По меньшей мере в одном варианте осуществления настоящего изобретения порог глюкозы зависит от времени суток. Кроме того, по меньшей мере в одном варианте осуществления изобретения сутки разбивают на временные блоки, по меньшей мере один из которых имеет иной порог глюкозы, нежели остальные временные блоки.
По меньшей мере в одном варианте осуществления настоящего изобретения порог вероятности выбирают из группы значений, состоящей из 99%, 98%, 90%, 75% и 50%.
По меньшей мере в одном варианте осуществления настоящего изобретения в период сна заданная доза инсулина выше, чем в период бодрствования.
По меньшей мере в одном варианте осуществления настоящего изобретения измеренные значения глюкозы получают путем периодического измерения глюкозы, выполняемого приблизительно каждую минуту.
По меньшей мере в одном варианте осуществления настоящего изобретения по меньшей мере одно из следующего: максимально допустимую глюкозу, минимальную дозу инсулина и поправочный коэффициент, получают из калькулятора болюсов.
По меньшей мере в одном варианте осуществления настоящего изобретения предупреждение может выдаваться посредством изображения, звука, вибрации или любой комбинации вышеперечисленного. Как возможный вариант, предупреждение также может указывать на необходимость измерения глюкозы.
По меньшей мере в одном варианте осуществления настоящего изобретения сгенерированное предупреждение и/или тревожное сообщение может быть временно отключено в течение некоторого периода времени, следующего за приемом пищи, корректирующим болюсом, или измерением глюкозы крови (ГК). По меньшей мере в одном варианте осуществления изобретения этот период времени может быть равен длительности времени задержки действия инсулина, определенного калькулятором болюсов, например, может составлять 15 минут, 30 минут, 45 минут, 1 час или 2 часа.
По меньшей мере в одном варианте осуществления настоящего изобретения предлагаемый способ может включать введение больному диабетом заданного страховочного болюса инсулина, если не проводится измерение глюкозы и не вводится корректирующая доза. Как возможный вариант, измерение глюкозы может представлять собой измерение глюкозы крови в рамках самоконтроля гликемии или дополнительное измерение, выполняемое системой непрерывного мониторирования гликемии и инициируемое пользователем.
По меньшей мере в одном варианте осуществления настоящего изобретения предлагаемый способ может также включать определение вычислительным устройством болюсной дозы инсулина с использованием результата измерения глюкозы.
По меньшей мере в одном варианте осуществления настоящего изобретения при определении болюсной дозы инсулина используют по меньшей мере одно из следующего: результат определения скорости изменения измеренных значений глюкозы, тенденцию в поведении измеренных значений глюкозы и закономерность в поведении измеренных значений глюкозы.
По меньшей мере в одном варианте осуществления настоящего изобретения при определении болюсной дозы инсулина используют по меньшей мере один результат определения скорости изменения измеренных значений глюкозы и прогнозное значение глюкозы, причем корректирующую дозу определяют по формуле:
Figure 00000005
где:
с - поправочный коэффициент,
g - измеренное значение глюкозы,
Figure 00000006
- скорость изменения глюкозы,
τ - прогнозное время, и
gt - прогнозный уровень глюкозы в момент времени t.
По меньшей мере в одном варианте осуществления настоящего изобретения время t составляет примерно от 20 до 60 минут.
По меньшей мере в одном варианте осуществления настоящего изобретения предлагаемый способ, выполняемый в контроллере, включает также определение максимального безопасного болюса инсулина.
По меньшей мере в одном варианте осуществления настоящего изобретения контроллер обеспечивает возможность выдачи на пользовательский интерфейс тревожного сообщения, если запрос болюса, введенный вручную посредством пользовательского интерфейса, превышает максимальный безопасный болюс инсулина.
По меньшей мере в одном варианте осуществления настоящего изобретения предлагаемый способ, осуществляемый при выполнении машиноисполняемых команд, включает также определение скорости изменения глюкозы на основании измеренных значений глюкозы, причем скорость изменения глюкозы используется калькулятором болюсов для модификации болюса инсулина.
По меньшей мере в одном варианте осуществления настоящего изобретения контроллер содержит также пользовательский интерфейс, предлагаемый способ включает также отображение на пользовательском интерфейсе модифицированного болюса инсулина.
По меньшей мере в одном варианте осуществления настоящего изобретения скорость изменения может использоваться для прогноза ожидаемого в будущем уровня глюкозы, предназначенного для использования калькулятором болюсов.
Осуществление изобретения рассматривается ниже на примерах, которые, являясь иллюстративными, не ограничивают возможностей осуществления изобретения и объема его охраны, определяемой его формулой. Приведенное ниже подробное описание иллюстративных вариантов осуществления изобретения поясняется перечисленными ниже чертежами, на которых однотипные элементы обозначены соответствующими ссылочными обозначениями:
На фиг. 1 приведена блок-схема выполнения способа анализа уровня глюкозы у больного диабетом по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 2 приведена блок-схема выполнения способа анализа уровня глюкозы у больного диабетом по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 3 схематически представлена система для анализа уровня глюкозы у больного диабетом по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 4 схематически представлена система непрерывного мониторирования гликемии (НМГ) по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 5 схематически представлен прибор контроля глюкозы по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 6 приведены графики результатов измерений глюкозы (измеренные значения глюкозы) и фактических уровней глюкозы человека по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 7 схематически представлены средство вероятностного анализа и рекурсивный фильтр по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 8 изображены переходы состояний для скрытой марковской модели по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 9 представлено функционирование скрытой марковской модели и фильтра Калмана по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 10 представлено функционирование прогнозирующего алгоритма по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 11 схематически представлен способ прогнозирования у человека уровня глюкозы с использованием средства вероятностного анализа и рекурсивного фильтра по меньшей мере в одном варианте осуществления настоящего изобретения.
На фиг. 12 приведена блок-схема предлагаемой в изобретении системы рекомендации болюсов инсулина по меньшей мере в одном иллюстративном варианте ее выполнения.
На фиг. 13 приведен график изменения глюкозы крови во времени, иллюстрирующий ряд примеров приема болюсов инсулина в отношении текущего и следующего (ближайшего) за ним временных интервалов.
На фиг. 14 приведена блок-схема по меньшей мере одного варианта программного алгоритма рекомендации болюсов инсулина, выполняемого представленной на фиг. 12 системой для определения и рекомендации количеств болюсного инсулина в соответствии с настоящим изобретением.
На фиг. 15-18 показаны примеры данных НМГ и реакция на них, получаемая в результате осуществления предлагаемых в изобретении способов по меньшей мере в одном варианте их осуществления.
На фиг. 19 показан пример данных НМГ и реакция на них, получаемая в результате осуществления предлагаемых в изобретении способов по меньшей мере в одном варианте их осуществления.
Описываемые варианты осуществления изобретения в целом относятся к способам и системам для обработки данных о содержании глюкозы, измеренных у больного диабетом и, в частности, для предупреждения больного диабетом, если вероятность того, что ему потребуется доза инсулина, превысит заданный порог. В контексте настоящего изобретения под "измеренными значениями глюкозы" понимаются уровни глюкозы человека, измеренные сенсором глюкозы, под "фактическим уровнем глюкозы" понимается действительный уровень глюкозы человека, а под "оценочным уровнем глюкозы" - уровень глюкозы человека, оцененный на основании результатов измерения глюкозы.
На фиг. 1 представлен способ 100 анализа уровня глюкозы у больного диабетом по меньшей мере в одном варианте осуществления изобретения. Рассматриваемый в качестве примера способ 100 включает следующие шаги: получение вычислительным устройством множества измеренных значений глюкозы от системы непрерывного мониторирования гликемии, подключенной к больному диабетом (рассматриваемый в качестве примера шаг 110 получения), анализ множества измеренных значений глюкозы, выполняемый в вычислительном устройстве с использованием средства вероятностного анализа для определения порога глюкозы (
Figure 00000007
) и граничного значения глюкозы (
Figure 00000008
) при пороге вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина (рассматриваемый в качестве примера шаг 120 анализа), и выполняемое вычислительным устройством сравнение граничного значения глюкозы (
Figure 00000008
) с порогом глюкозы (
Figure 00000009
), причем если граничное значение глюкозы (
Figure 00000008
) больше или равно порогу глюкозы (
Figure 00000010
), то вычислительное устройство выдает через пользовательский интерфейс предупреждение (рассматриваемый в качестве примера шаг 130 сравнения). По меньшей мере в одном варианте осуществления способа 100 вычислительное устройство при необходимости может определять максимальный безопасный болюс инсулина (рассматриваемый в качестве примера шаг 125). Если вручную введен запрос болюса, который превысит максимальный безопасный болюс, выдается тревожное сообщение (рассматриваемый в качестве примера шаг 135 выдачи тревожного сообщения). По меньшей мере в одном варианте осуществления шага 135 тревожное сообщение может включать в себя извещение, адресованное больному диабетом, и несколько вариантов реагирования на это извещение.
Вычислительное устройство, используемое при осуществлении способа 100 по меньшей мере в одном его варианте, может включать в себя глюкометр, инсулиновую помпу, микропроцессор, связанный с системой непрерывного мониторирования гликемии, микроконтроллер, сотовый телефон, смартфон, персональный цифровой ассистент (карманный компьютер), персональный компьютер или компьютерный сервер. При рассмотрении в данном описании некоторых вариантов осуществления изобретения вычислительное устройство также может называться прибором контроля глюкозы, или монитором глюкозы.
Система непрерывного мониторирования гликемии, используемая при осуществлении способа 100 по меньшей мере в одном его варианте, физически подключена к больному диабетом и выполнена с возможностью автоматического измерения его уровня глюкозы. Результаты измерения глюкозы могут поступать от системы непрерывного мониторирования гликемии с обычными интервалами, например, каждую минуту, каждые 2 минуты, каждые 3 минуты, каждые 4 минуты, каждые 5 минут, каждые 10 минут или каждые 15 минут.
На рассматриваемом в качестве примера шаге 110 получения в вычислительное устройство могут поступать дополнительные источники данных, используемых при определении порога глюкозы и/или граничного значения глюкозы на шаге 120 анализа. Например, к таким дополнительным данным могут относиться максимально допустимая глюкоза, минимальная доза инсулина и поправочный коэффициент, которые могут использоваться при определении порога глюкозы. По меньшей мере некоторые из дополнительных источников данных можно получать от калькулятора болюсов.
Максимально допустимая глюкоза, по меньшей мере в одном варианте осуществления изобретения, выдается калькулятором болюсов. Этот рассматриваемый в качестве примера показатель является либо средним значением верхней и нижней целевых границ, определяющих целевой диапазон глюкозы, либо верхней границей трапецоида, следующего за каким-либо значительным возмущением, таким как прандиальная, т.е. пищевая (принимаемая на прием пищи) или корректирующая дозы инсулина. После получения высокого измеренного значения глюкозы крови (ГК) корректирующая доза инсулина может скорректировать уровень ГК, приведя его к максимально допустимому значению глюкозы. Корректирующей дозой также может быть доза, после которой дальнейшая коррекция не рекомендуется. По меньшей мере в одном варианте осуществления изобретения калькулятор болюсов также содержит коэффициент чувствительности к инсулину или поправочный коэффициент.
Далее, по меньшей мере в одном варианте осуществления изобретения может использоваться рекурсивный фильтр, предоставляющий вычислительному устройству вектор состояния, включающий в себя значение глюкозы, скорость изменения глюкозы и ускорение глюкозы. Кроме того, на шаге 110 получения рекурсивный фильтр может предоставлять вычислительному устройству ковариационную матрицу. Ковариационная матрица содержит меру точности вектора состояния.
Порог глюкозы (
Figure 00000011
) по меньшей мере в одном варианте осуществления шага 120 анализа может вычисляться при помощи уравнения [1]:
Figure 00000012
где:
gm - максимально допустимая глюкоза,
с - поправочный коэффициент,
I0 - минимальная доза инсулина.
Минимальная доза инсулина может быть задана по меньшей мере в одном варианте осуществления изобретения как определенное количество инсулина (например: 0,5, 1, 2, 3, 4 или 5 единиц) или как процентная доля полной суточной дозы больного диабетом (например 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9% или 10%).
Граничное значение глюкозы (
Figure 00000008
) по меньшей мере в одном варианте осуществления шага 120 анализа может вычисляться при помощи уравнения [2]:
Figure 00000013
где:
Figure 00000014
- вариация при измерении глюкозы,
ρ - порог вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина, и
g - результат определения глюкозы, полученный на основании множества измеренных значений глюкозы.
Порог вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина, может составлять, по меньшей мере в одном варианте осуществления настоящего изобретения, 99%, 98%, 95%, 90%, 75% или 50%.
По меньшей мере в одном варианте осуществления настоящего изобретения заданная доза инсулина может зависеть от времени суток. Например, сутки можно разбить на множество временных блоков, по меньшей мере один из которых имеет иной порог глюкозы, нежели по меньшей мере один из остальных временных блоков. Кроме того, в период сна заданная доза инсулина может быть выше, чем в период бодрствования. Заданная доза инсулина, по меньшей мере в одном варианте осуществления изобретения, может составлять 0,5, 1, 2, 3, 4 или 5 единиц.
По меньшей мере в одном варианте осуществления шага 130 сравнения предупреждение может выдаваться посредством изображения, звука, вибрации, любой комбинации вышеперечисленного или любого иного раздражителя, способного побудить пользователя к выполнению определенного задания. Кроме того, по меньшей мере в одном варианте осуществления изобретения предупреждение может указывать на необходимость измерения глюкозы пользователем. Такое измерение глюкозы может представлять собой измерение глюкозы крови в рамках самоконтроля гликемии или дополнительное измерение, выполняемое системой непрерывного мониторирования гликемии. Дополнительное измерение, выполняемое системой непрерывного мониторирования гликемии, может осуществляться системой автоматически или может инициироваться пользователем вручную (в ручном режиме).
На дозирование инсулина также может оказывать влияние ситуация, называемая "наложением инсулина", когда человек дозирует себе инсулин, не учитывая тот инсулин, который еще остается действующим, или текущий уровень глюкозы крови. По меньшей мере в одном варианте осуществления настоящего изобретения предлагается то, как можно использовать НМГ для определения максимального безопасного болюса инсулина. Если введенный вручную запрос болюса превышает по количеству инсулина максимальный безопасный болюс, то человек уведомляется об этом, и ему предлагается на выбор набор возможных действий. Это мероприятие препятствует практике введения болюсов на основе показателя НМГ и уменьшает риск возникновения наложения инсулина.
По меньшей мере в одном варианте осуществления шага 125 алгоритм, используемый для выдачи тревожного сообщения об опасности наложения инсулина, запускается выбором введения болюса вручную. При выборе количества болюсного инсулина контроллер получает значение максимально допустимой глюкозы, обеспечиваемое калькулятором болюсов. Это значение является либо средним значением верхней и нижней целевых границ, определяющих целевой диапазон глюкозы, либо верхней границей трапецоида, следующего за прандиальной и корректирующей дозами инсулина. После получения высокого измеренного значения глюкозы крови (ГК) корректирующая доза инсулина может скорректировать уровень ГК, приведя его к максимально допустимому значению глюкозы. Калькулятор болюсов также содержит коэффициент чувствительности к инсулину (КЧИ), или поправочный коэффициент.
gm максимально допустимая глюкоза
с поправочный коэффициент или КЧИ (т.е. 30 мг/дл/МЕ)
Выдача тревожного сообщения может определяться, по меньшей мере в одном варианте осуществления шага 135, двумя параметрами. Первым является порог вероятности, который должен быть преодолен для выдачи тревожного сообщения. Это определяется вероятностью того, что измеренное значение будет ниже порога. Вторым является максимальная безопасная доза инсулина. Это порог для выдачи тревожного сообщения, и он вычисляется на основе порога вероятности.
ρL Нижний порог вероятности (т.е. 5%)
ρU Верхний порог вероятности (т.е. 95%)
gL Уровень глюкозы для нижнего порога вероятности
gU Уровень глюкозы для верхнего порога вероятности
IS Максимальная безопасная доза инсулина
IU Верхняя или небезопасная доза инсулина. Это доза, превышающая верхние 95% области неопределенности глюкозы.
Для определения необходимости выдачи тревожного сообщения вычисляется граница для нижней квантили ρ.
Figure 00000015
Figure 00000016
Из этой нижней процентили распределения неопределенности глюкозы вычисляется максимальная безопасная доза инсулина:
Figure 00000017
В случае неравенства (I>IS) заданный вручную болюс находится выше максимальной безопасной дозы инсулина. В этом случае больной диабетом уведомляется и получает на выбор набор возможных мер. К таким мерам может относиться по меньшей мере одна из следующих:
- Провести самоконтроль гликемии (СГК) и задействовать калькулятор болюсов.
- Уменьшить заданный вручную болюс до максимальной безопасной дозы инсулина.
- Отменить болюс.
- Ввести количество углеводов для соотнесения с дозой.
- Пометить болюс как прандиальный (вводимый перед приемом пищи), не вводя количество углеводов.
- Игнорировать тревожное сообщение и ввести болюс.
Небезопасная доза инсулина вычисляется по верхней процентили распределения неопределенности глюкозы:
Figure 00000018
В случае неравенства (I>IU) заданный вручную болюс находится в небезопасной области значений дозы инсулина. В этом случае больной диабетом уведомляется и получает на выбор набор возможных мер. К таким мерам может относиться по меньшей мере одна из следующих:
- Провести СГК и задействовать калькулятор болюсов.
- Уменьшить заданный вручную болюс до максимальной безопасной дозы инсулина.
- Отменить болюс.
- Ввести количество углеводов для соотнесения с дозой.
- Пометить болюс как прандиальный, не вводя количество углеводов.
Способ 100, по меньшей мере в одном варианте его осуществления, может дополнительно включать шаг введения больному диабетом страховочного болюса инсулина, если не проводится измерение глюкозы крови и не вводится корректирующая доза (рассматриваемый в качестве примера шаг 140 введения). В рассматриваемом в качестве примера варианте осуществления изобретения страховочный болюс может определяться вычислительным устройством с использованием граничного значения глюкозы.
Способ 100, по меньшей мере в одном варианте его осуществления, также может включать шаг определения вычислительным устройством болюсной дозы инсулина с использованием результата измерения глюкозы (рассматриваемый в качестве примера шаг 150 определения). По меньшей мере в одном варианте осуществления изобретения вычислительное устройство может содержать калькулятор болюсов инсулина, выполненный, например, в рассматриваемом в данном описании варианте его осуществления. По меньшей мере в одном варианте осуществления изобретения вычислительное устройство может определять препрандиальную дозу инсулина. На шаге 150 определения болюсной дозы инсулина может использоваться по меньшей мере одно из следующего: результат определения скорости изменения множества измеренных значений глюкозы, тенденция в поведении множества измеренных значений глюкозы и закономерность в поведении множества измеренных значений глюкозы.
В рассматриваемом в качестве примера варианте осуществления шага 150 определения дозы инсулина на этом шаге может использоваться по меньшей мере один результат определения скорости изменения множества измеренных значений глюкозы и прогнозное значение глюкозы, причем корректирующую дозу определяют при помощи уравнения [3]:
Figure 00000019
где:
с - поправочный коэффициент,
g - результат определения глюкозы, полученный на основании множества измеренных значений глюкозы,
Figure 00000020
- скорость изменения глюкозы,
τ - прогнозное время, и
gt - прогнозный уровень глюкозы в момент времени t.
По меньшей мере в одном варианте осуществления изобретения время t, используемое на рассматриваемом в качестве примера шаге 150, может составлять примерно 20 минут, примерно 30 минут, примерно 40 минут, примерно 50 минут или примерно 60 минут.
По меньшей мере в одном варианте осуществления шага 150 определения дозы инсулина может использоваться поверхность риска. Этот вариант служит для расчета поправки к максимально допустимой глюкозе. На шаге 150, по меньшей мере в этом варианте его осуществления, поверхность риска используется для задания кривой, где эта поверхность риска является нулевой. В этом случае риск гипогликемии указывает на то, что поверхность риска отрицательная, а риск гипергликемии указывает на то, что поверхность риска положительная. Определяют оптимальное значение глюкозы (g0opt), когда скорость изменения глюкозы является нулевой. Определяют также оптимальное значение глюкозы (g1opt) Для текущей скорости ее изменения. Если g1opt≤g0opt, то поправку к максимально допустимой глюкозе устанавливают равной нулю. Если g1opt>g0opt, к максимально допустимой глюкозе прибавляют разность (g1opt-g0opt). Впоследствии для вычисления порога глюкозы используют уточненную таким образом максимально допустимую глюкозу.
По меньшей мере в одном варианте своего осуществления настоящее изобретение может быть реализовано на машиночитаемом носителе информации с хранящимися в нем машиноисполняемыми командами, предназначенными для осуществления способа анализа уровня глюкозы у человека. Способом анализа уровня глюкозы у человека может быть любой вариант осуществления способа 100, рассматриваемого в данном описании.
На фиг. 2 схематически представлен рассматриваемый в качестве примера вариант осуществления способа 200, осуществляемого машиноисполняемыми командами. Рассматриваемый в качестве примера способ 200 включает следующие шаги: получение вычислительным устройством множества измеренных значений глюкозы от системы непрерывного мониторирования гликемии, подключенной к больному диабетом (рассматриваемый в качестве примера шаг 210 получения); анализ множества измеренных значений глюкозы, выполняемый в вычислительном устройстве с использованием средства вероятностного анализа для определения граничного значения глюкозы (
Figure 00000008
) при пороге вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина, и порога глюкозы (
Figure 00000021
) (рассматриваемый в качестве примера шаг 220 анализа); выполняемое вычислительным устройством сравнение граничного значения глюкозы (
Figure 00000008
) с порогом глюкозы (
Figure 00000022
), причем если граничное значение глюкозы (
Figure 00000008
) превышает порог глюкозы (
Figure 00000023
), то вычислительное устройство выдает через пользовательский интерфейс предупреждение (рассматриваемый в качестве примера шаг 230 сравнения). По меньшей мере в одном варианте осуществления способа 200 вычислительное устройство при необходимости может определять максимальный безопасный болюс инсулина (рассматриваемый в качестве примера шаг 225). Если вручную введен запрос болюса, который превысит максимальный безопасный болюс, выдается тревожное сообщение (рассматриваемый в качестве примера шаг 235 выдачи тревожного сообщения). По меньшей мере в одном варианте осуществления шага 235 тревожное сообщение может включать в себя извещение, адресованное больному диабетом, и несколько вариантов реагирования на это извещение.
По меньшей мере в одном варианте осуществления способа 100 или 200, предупреждение, сгенерированное на шаге 130/230, и/или тревожное сообщение, сгенерированное на шаге 135/235, может быть временно отключено в течение периода времени, следующего за приемом пищи, корректирующим болюсом или измерением ГК. По меньшей мере в одном варианте осуществления изобретения этот период времени может быть равен длительности времени задержки действия инсулина, определенного калькулятором болюсов, например, может составлять 15 минут, 30 минут, 45 минут, 1 час или 2 часа.
На фиг. 3 схематически представлена система 300 для анализа уровня глюкозы у больного диабетом по меньшей мере в одном варианте осуществления настоящего изобретения. Эта рассматриваемая в качестве примера система 300 включает в себя контроллер 310, содержащий машиночитаемый носитель 320 информации с хранящимися в нем машиноисполняемыми командами, предназначенными для осуществления настоящего изобретения. По меньшей мере в одном варианте осуществления изобретения при выполнении машиноисполняемых команд осуществляется способ, включающий: получение вычислительным устройством множества измеренных значений глюкозы от сенсора глюкозы, подключенного к больному диабетом; анализ множества измеренных значений глюкозы, выполняемый в вычислительном устройстве с использованием средства вероятностного анализа для определения вероятности (PID) того, что больному диабетом потребуется по меньшей мере заданная доза инсулина; и выполняемое вычислительным устройством сравнение вероятности (PID) с пороговой вероятностью (PT), причем если вероятность (PID) выше пороговой вероятности (PT), то вычислительное устройство выдает предупреждение.
Рассматриваемая в качестве примера система 300 также может включать в себя связанный с контроллером калькулятор 330 болюсов, способный определять болюс инсулина на основании множества измеренных значений глюкозы.
Способ, выполняемый машиноисполняемыми командами 320, по меньшей мере в одном варианте осуществления настоящего изобретения также включает определение скорости изменения глюкозы на основании множества измеренных значений глюкозы, причем скорость изменения глюкозы используется калькулятором болюсов для модификации болюса инсулина.
По меньшей мере в одном варианте осуществления изобретения контроллер 310 содержит также пользовательский интерфейс 340, связанный с контроллером 310. При этом выполняемый системой 300 способ может также включать шаг отображения на пользовательском интерфейсе модифицированного болюса инсулина. Кроме того, скорость изменения глюкозы, полученная на основании множества результатов определения глюкозы, может использоваться для прогноза ожидаемого в будущем уровня глюкозы, предназначенного для использования калькулятором болюсов.
Средства непрерывного мониторирования гликемии
На фиг. 4 представлен вариант выполнения системы непрерывного мониторирования гликемии (НМГ), используемой по меньшей мере в одном варианте осуществления способов и систем, рассматриваемых в данном описании. Рассматриваемая в качестве примера система 410 может включать в себя сенсор 416 глюкозы, имеющий иглу 418, вводимую под кожу 412 больного диабетом. Кончик иглы 418 может находиться в интерстициальной жидкости 414, вследствие чего результаты измерений, выполняемых сенсором 416 глюкозы, основываются на уровне глюкозы в интерстициальной жидкости 414. Сенсор 416 глюкозы может быть поставлен человеку на живот или на другое подходящее для этого место и может быть зафиксирован клейкой лентой (пластырем) или клеем (не показаны). Кроме того, сенсор 416 глюкозы может периодически калиброваться для повышения его точности. Эта периодическая калибровка может способствовать коррекции дрейфа показаний сенсора, обусловленного ухудшением характеристик сенсора и изменениями в физическом состоянии места его установки. Сенсор 416 глюкозы может содержать и другие компоненты, в том числе (без ограничения) беспроводной передатчик 420 и антенну 422. Хотя на фиг. 4 показано, что сенсор 416 глюкозы имеет прямоугольную форму, предполагается, что он может иметь и другие геометрические формы. Хотя для доступа к крови или иной физиологической жидкости человека в сенсоре 416 глюкозы используется игла 418, для доступа к крови или иной жидкости человека с целью измерения глюкозы в нем могут использоваться и другие подходящие приборы, в том числе еще не разработанные.
После проведения измерения сенсор 416 глюкозы может передавать измеренное значение глюкозы в прибор 426 контроля глюкозы (монитор гликемии) с дисплеем 428 и антенной 430 по линии 424 связи. Линия 424 связи может беспроводной, например, работающей на радиочастотах, и результаты измерения глюкозы передаются по такой линии связи посредством электромагнитных волн. Примером системы беспроводной радиосвязи является система типа Bluetooth®, использующая частоту, приблизительно составляющую 2,4 гигагерца (ГГц). В качестве еще одного примера схемы беспроводной связи можно назвать системы беспроводной связи с использованием инфракрасного света, например системы, поддерживаемые Ассоциацией инфракрасной передачи данных (InfraRed Data Association®, сокр. IrDA®). Могут быть предусмотрены и другие подходящие типы беспроводной связи, включая существующие и еще не разработанные технологии. Линия 424 связи может быть однонаправленной или односторонней (т.е. данные могут передаваться только от сенсора 416 глюкозы в прибор 426 контроля глюкозы) либо двунаправленной или двусторонней (т.е. данные могут передаваться между сенсором 416 глюкозы и прибором 426 контроля глюкозы в обоих направлениях). Кроме того, линия 424 связи может обеспечивать передачу данных между двумя или более устройствами (например, между сенсором глюкозы, прибором контроля глюкозы, инсулиновой помпой и т.д.). Хотя на фиг. 4 показана беспроводная линия 424 связи, в качестве альтернативы ей может быть предусмотрена проводная линия связи, например, канал проводной связи по технологии Ethernet. Могут использоваться и другие общественные или частные проводные или беспроводные линии связи.
На фиг. 5 представлен прибор 526 контроля глюкозы, который по меньшей мере в одном варианте его выполнения может содержать дисплей 528, микроконтроллер 532 и устройство 534 ввода. К неограничивающим примерам реализации прибора 526 контроля глюкозы относятся, глюкометр, инсулиновая помпа, сотовый телефон, смартфон, персональный цифровой ассистент (карманный компьютер), персональный компьютер или компьютерный сервер. Микроконтроллер 532 может быть электрически связан с устройством 534 ввода, которое может быть выполнено с возможностью приема множества измеренных значений глюкозы от подключенного к организму человека сенсора глюкозы. Микроконтроллер 532 может быть выполнен с возможностью приема множества измеренных значений глюкозы из устройства 534 ввода. Прибор 526 контроля глюкозы также может быть выполнен с возможностью хранения в запоминающем устройстве (не показано) множества измеренных значений глюкозы, полученных от сенсора 416 глюкозы, на протяжении определенного периода времени. Микроконтроллер 532 также может быть выполнен с возможностью анализа множества измеренных значений глюкозы при помощи средства вероятностного анализа, обеспечивающего определение вероятности неправильной работы сенсора глюкозы. Кроме того, микроконтроллер 532 может быть выполнен с возможностью оценивания уровня глюкозы человека с использованием рекурсивного фильтра, способного взвешивать множество измеренных значений глюкозы по вероятности точности показаний сенсора глюкозы. Наконец, микроконтроллер 532 может быть электрически связан с дисплеем 528 с возможностью передачи на дисплей 528 информации, относящейся к оценочному уровню глюкозы человека. Как обсуждается в данном описании, отображаемая на дисплее информация может включать в себя оценочный уровень глюкозы у человека и/или уровень глюкозы, прогнозируемый у человека в определенный момент времени в будущем, т.е. прогнозный уровень глюкозы. Далее, отображаемая на дисплее информация также может включать в себя оценку качества или неопределенности оценочного уровня глюкозы. Кроме того, отображаемая на дисплее информация может включать в себя уведомления, предупреждения и т.д. в отношении того, что оцениваемые или прогнозируемые уровни глюкозы человека являются гипогликемическими или станут таковыми в определенный момент времени в будущем. Это может происходить, например, при снижении (или прогнозировании снижения) уровня глюкозы человека ниже заданного порога гипогликемии, составляющего, например, 50 миллиграммов глюкозы на децилитр крови (мг/дл). Далее, отображаемая на дисплее 528 информация может включать в себя тот или иной вариант предупреждения, если вероятность того, что человеку потребуется заданная доза инсулина, превышает заданный порог.
Микроконтроллер 532 может быть 8-разрядным, 16-разрядным или любым иным подходящим устройством и может иметь внутрикристальные периферийные устройства, обеспечивающие его функционирование. Например, микроконтроллер 532 может иметь внутреннюю память для хранения компьютерной программы, внутреннюю память для хранения данных и внутреннюю память для энергонезависимого хранения параметров, используемых микроконтроллером во время его работы. Кроме того, микроконтроллер 532 может иметь таймеры, последовательные коммуникационные порты, контроллеры прерываний и тому подобное. Микроконтроллер 532 также может входить в состав специализированной заказной интегральной схемы (ASIC), которая может содержать другие схемы, обеспечивающие функционирование прибора 426 контроля глюкозы. Устройство 534 ввода может быть беспроводным коммуникационным модулем, выполненным с возможностью беспроводного приема измеряемых значений глюкозы от сенсора глюкозы (как показано на фиг. 4). При этом для повышения надежности беспроводного соединения может использоваться антенна 530. В качестве альтернативы использованию беспроводных технологий, устройство 534 ввода может быть коммуникационным модулем с возможностью приема измеряемых значений глюкозы посредством проводного соединения, например, по стандарту Ethernet или аналогичному протоколу. Дисплей 528 может включать в себя жидкокристаллический (ЖК) экран или иную подходящую технологию. Дисплей 528 также может быть выполнен с возможностью тактильной передачи информации человеку, например, за счет вибрирования.
На фиг. 6 приведен пример двумерного графика измеренных значений 640 глюкозы, полученных от сенсора глюкозы, подключенного к больному диабетом. По горизонтальной ось отложено время в часах, а по вертикальной оси - измеренный и фактический уровни глюкозы человека в миллиграммах глюкозы на децилитр крови (мг/дл). Измерения может выполняться сенсором глюкозы автоматически с заданной частотой, например, каждую минуту. Измеренные значения 640 глюкозы в общем случае могут соответствовать фактическому уровню 642 глюкозы человека. Однако, время от времени сенсор глюкозы может функционировать неправильно, например, на временном отрезке 644, показанном на фиг. 6. Неправильное функционирование может быть обусловлено пропаданием сигнала от сенсора, вызванным физиологическими проблемами вблизи иглы сенсора глюкозы, или электрическими проблемами с самим сенсором глюкозы. На временном отрезке 644 неправильной работы сенсора измеряемые сенсором значения 640 глюкозы могут быть существенно ниже фактического уровня 642 глюкозы человека. В конце временного отрезка 644 неправильной работы измеренные значения 640 глюкозы могут восстановиться и снова соответствовать фактическому уровню 642 глюкозы человека.
Кроме того, время от времени сенсор глюкозы также может создавать шум, например, на временном отрезке 646, показанном на фиг. 6. Такой шум может быть обусловлен физическими перемещениями иглы сенсора глюкозы относительно кожи или электрическим шумом, возникающим в самом сенсоре глюкозы. На временном отрезке 646 присутствия шума измеренные значения 640 глюкозы могут колебаться, в результате чего одни результаты измерения оказываются выше, а другие - ниже фактического уровня 642 глюкозы. При наличии шума в сигнале сенсора результаты измерения могут колебаться даже вблизи фактического уровня 642 глюкозы. В конце временного отрезка 646 присутствия шума измеренные значения 640 глюкозы могут восстановиться и снова прийти в близкое соответствие фактическому уровню 642 глюкозы человека. Хотя временные отрезки неправильной работы сенсора и присутствия шума в его сигнале на графике разнесены во времени, они также могут совпадать во времени, а также следовать непосредственно друг за другом. Кроме того, длительность неправильной работы сенсора и/или присутствия шума в его сигнале может быть меньшей или большей, чем показанная на фиг. 6.
Именно с учетом рассмотренных выше проблем уровня техники в настоящем изобретении предлагается оценивать фактический уровень глюкозы человека при наличии шума в сигнале сенсора и/или при неправильной работе сенсора. На фиг. 7 изображена система 750, выполненная с возможностью оценивания уровня глюкозы человека. Система 750 может получать измеренные значения 640 глюкозы от сенсора 416 глюкозы, который может быть подключен к больному диабетом (не показан). Сенсор 416 глюкозы может быть подключен к организму человека описанным выше или любым иным подходящим образом. Сенсор 416 глюкозы может быть выполнен с возможностью периодического, например, выполняемого каждую минуту, измерения уровня глюкозы человека и передачи измеренных значений 640 глюкозы в систему 750 (например, по описанной выше линии связи).
Система 750 для оценивания уровня глюкозы у больного диабетом может включать в себя средство 754 вероятностного анализа и рекурсивный фильтр 752. Средство 754 вероятностного анализа может быть выполнено с возможностью получения измеренных значений 640 глюкозы и вычисления вероятности PA 758 точности сенсора глюкозы, т.е. вероятности того, что сенсор 416 глюкозы функционирует нормально, или правильно. Вероятность PA 758 точности сенсора глюкозы может быть основана только на наблюдаемых данных, таких как измеренные значения 640 глюкозы и/или их изменения. Соответственно, средство 754 вероятностного анализа может использоваться для проведения различия между шумом сенсора, который может иметь нормальное распределение, и неправильной работой сенсора, при которой результаты измерений распределяться иначе, чем по нормальному закону. С каждым типом неопределенности можно обращаться по-своему ввиду различий в распределениях неопределенности.
Средство 754 вероятностного анализа может включать в себя любое число математических алгоритмов, способных анализировать измеренные значения 640 глюкозы и/или их изменения и вычислять вероятность PA 758 точности сенсора глюкозы. Средство 754 вероятностного анализа также может быть выполнено с возможностью получения других типов данных, на основании которых можно определять вероятность PA 758 точности сенсора глюкозы, например, с возможностью получения информации о том, когда человек принимает пищу, когда человек занимается физическими упражнениями и когда человеку вводится инсулин. Кроме того, вероятность PA 758 может быть основана на одном или нескольких результатах измерений, получаемых от прибора измерения импеданса, подключенного к телу человека и выполненного с возможностью измерения у человека импеданса. Такой прибор измерения импеданса может показывать, правильно ли сенсор 416 глюкозы подключен к организму человеку, а значит, точны ли измеренные значения 640 глюкозы, получаемые от сенсора 416 глюкозы. Могут использоваться и другие типы данных.
Средство 754 вероятностного анализа может иметь ряд различных форм, таких как конечный автомат, байесовские модели или другие алгоритмы. В одном варианте осуществления изобретения средство 754 вероятностного анализа может иметь форму простого конечного автомата, в котором вероятность Ра точности сенсора глюкозы всегда может находиться в наборе значений {0,1} (т.е. вероятность PA 758 равна либо 0%, либо 100% в зависимости от состояния конечного автомата). В данном примере система перейдет в состояние неточности сенсора, TA→I, если ΔCG (т.е. изменение текущего измеренного значения глюкозы относительно предыдущего измеренного значения глюкозы) окажется меньше некоторого отрицательного порога, τ1, и перейдет обратно в состояние точности сенсора, TI→A, если ΔCG окажется больше некоторого положительного порога, τ2, или если полученное от сенсора значение CG (т.е. текущее измеренное значение глюкозы) находится в пределах физиологически возможных значений глюкозы (g0 и gmax), и с момента перехода в состояние неточности сенсора прошло определенное время, ΔTA→I>τ3.
Математически это можно представить следующим образом:
TI→A если ΔCG<τ1
TA→I если ΔCG>τ2 или (g0<CG<gmax и ΔTA→I>τ3)
Если не выполнено ни одно из этих условий перехода, то конечный автомат может остаться в своем текущем состоянии. Это лишь один пример того, как средство 754 вероятностного анализа может быть реализовано в форме конечного автомата. Средство 754 вероятностного анализа также может быть реализовано в иных формах.
В еще одном варианте осуществления изобретения средство 754 вероятностного анализа может включать в себя скрытую марковскую модель, имеющую для сенсора глюкозы два состояния: 1) состояние, в котором сенсор глюкозы является точным, т.е. выдает точные показания, и это состояние обозначено как "SA", и 2) состояние, в котором сенсор глюкозы является неточным, т.е. выдает неточные показания, и это состояние обозначено как "SI" Скрытая марковская модель может обеспечивать функции переходов между состояниями, определяющих вероятность перехода из состояния в состояние SA, например, следующую функцию:
Figure 00000024
,
где "CG" - текущее измеренное значение глюкозы, "ΔCG" - изменение измеренного значения глюкозы от предыдущего до текущего, и α1-α4 - константы, зависящие от характеристик сенсора глюкозы. Диапазон выходных значений этой функции составляет от ноля до единицы, где ноль означает, что вероятность точности сенсора равна 0%, а единица - что вероятность точности сенсора равна 100%. Функция "min" выбирает наименьшее из значения математического выражения и единицы (т.е. вероятности 100%). Эта функция перехода может быть основана на текущем значении CG и значении ΔCG. Кроме того, функция перехода может быть сигмоидой, в которой параметры α1 и α3 определяют положение сигмоидального перехода, а параметры α2 и α4 определяют наклон сигмоиды. Эти параметры могут быть настроены под конкретного пациента и/или партию сенсоров.
Продолжая рассматривать пример скрытой марковской модели, вероятность сохранения состояния SI (когда текущим состоянием является состояние SI) может определяться выражением:
Figure 00000025
,
в котором эта вероятность является функцией лишь значения ΔCG и предыдущей вероятности
Figure 00000026
нахождения в состоянии SI или перехода в это состояние. Диапазон выходных значений этой функции составляет от ноля до единицы, где ноль означает, что вероятность точности равна 0%, а единица - что вероятность точности сенсора равна 100%. Функция "max" выбирает наибольшее из значения математического выражения и ноля (т.е. вероятности 0%). Параметр γ - это ограничивающий параметр, меньший единицы и предназначенный для обеспечения постепенного перехода состояния скрытой марковской модели обратно к состоянию SA, если значения CG и ΔCG не позволяют полагать, что модель остается в состоянии SI. Параметр γ может быть константой и может соотноситься с вероятностью сохранения состояния SI, когда значение ΔCG является относительно нормальным. Например, параметр у может выбираться с таким расчетом, чтобы, когда значение ΔCG является относительно нормальным, скрытая марковская модель оставалась в состоянии SI в течение приблизительно 10 минут. Эта функция вероятности также включает в себя сигмоидальную функцию, обнаруживающую быстрые подъемы сигнала CG, связанные с возвратом в состояние SA. Параметр α5 определяет положение сигмоидального перехода, а параметр α6 определяет наклон сигмоиды. Оба этих параметра могут быть настроены под конкретного пациента и/или партию сенсоров.
Текущую вероятность PI перехода в состояние SI можно обозначить как PA→I или PI→I в зависимости от того, является ли текущим состоянием SA или SI. Текущую вероятность PI того, что сенсор глюкозы неточен (т.е. находится в состоянии SI), можно выразить как (SA×PA→I)+(SI×PI→I). Отметим, что одно или другое состояние (SA или SI) равно "1" при нахождении в этом состоянии и равно "0" в других случаях. Это включает в себя вероятность перехода в SI (PA→I) с учетом вероятности нахождения в SA и вероятность сохранения SI, умноженную на вероятность текущего нахождения в SI. Значение PI→I равно 1-PI→A, а вероятность того, что сенсор точен, просто равна PA=1-PI. Таким образом, для данного примера вероятность точности сенсора глюкозы может быть выражена как PA=1-[(SA×PA→I)+(SI×PI→I)].
На фиг. 8 приведено графическое представление двух функций перехода, PA→I и PI→A (т.е. 1-PI→A, вероятность перехода из SI в SA, когда текущим состоянием является SI), наложенное на гистограмму ΔCG. Эта гистограмма имеет компоненту 868 гауссовской формы (формы нормального распределения) с центром в ноле и двумя хвостами, связанными с входом в состояние неправильной работы сенсора и выходом из этого состояния. Две функции перехода нанесены поверх гистограммы, показывая, что их можно настроить таким образом, чтобы они определяли начало хвостов гистограммы. Компонента 868 гауссовской формы может представлять диапазон значений ΔCG, получаемых при нормальной работе сенсора глюкозы. Значения ΔCG, находящиеся в пределах компоненты 868 гауссовской формы, могут быть обусловлены, например, шумом сенсора. Значения ΔCG, находящиеся за пределами компоненты 868 гауссовской формы и слева от нее, могут быть обусловлены переходом сенсора из состояния SA в состояние SI. Форма этого распределения может использоваться для характеристики партии сенсоров глюкозы после производства и может использоваться для кодирования сенсоров. То есть, функции перехода (PA→I и PI→A) можно корректировать (путем регулировки параметров араб) таким образом, чтобы они соответствовали компоненте 868 гауссовской формы для конкретной партии сенсоров глюкозы. Таким образом, скрытая марковская модель может использоваться для определения вероятности того, что сенсор выдает точные показания, PA, только на основании измеренных значений глюкозы и их изменений.
На фиг. 9 показан график, иллюстрирующий пример функционирования скрытой марковской модели во время неправильной работы сенсора глюкозы и в присутствии шума сенсора глюкозы. На графике представлены измеренные значения 640 глюкозы и вероятность PA 758 точности сенсора глюкозы. В период 970 сенсор глюкозы может начать функционировать неправильно, вследствие чего измеренные значения 640 глюкозы станут неточными; одновременно вероятность PA 758 (определенная посредством скрытой марковской модели) может уменьшиться от приблизительно 100% (до периода 970) до около 0% в период 974. Это уменьшение вероятности может быть обусловлено тем, что скрытая марковская модель обнаружила быстрое снижение измеренных значений 640 глюкозы в начале периода 970 (т.е. при первом проявлении неправильной работы). В конце периода 974 сенсор глюкозы может начать работать нормально (т.е. измеренные значения 640 глюкозы снова становятся точными) и вероятность PA 758 может снова увеличиться до приблизительно 100%. Как и до того, это увеличение может быть обусловлено тем, что скрытая марковская модель обнаружила быстрое повышение измеренных значений 640 глюкозы в конце периода 974 (т.е. при возвращении сенсора глюкозы к нормальной работе). Скорость изменения вероятности PA 758 от около 0% до приблизительно 100% может зависеть от того, как быстро сенсор глюкозы переходит от неправильной (неточной) работы к нормальной (точной) работе. Если этот переход относительно быстр, вероятность PA 758 может быстро перейти от около 100% до приблизительно 0%. Если же сенсор глюкозы переходит от неправильной к нормальной работе медленно, переход вероятности PA 758 от около 0% до приблизительно 100% также может быть медленным. Ограничивающий параметр у (входящий в уравнение PI→I) может обеспечивать постепенный переход вероятности PA 758 обратно к состоянию SA, если значения CG и ΔCG дают слабые основания или вообще не позволяют полагать, что модель остается в состоянии SI.
Кроме того, на фиг. 9 показано, что к уменьшению вероятности PA 758 также может приводить наличие шума в показаниях сенсора глюкозы, имеющее место на графике в период 976, причем уменьшение вероятности может зависеть от интенсивности и уровня шума. Как показано на фиг. 9, в период 976 шум сенсора глюкозы может порождать небольшое уменьшение вероятности PA 758. Разумеется, как неправильная работа сенсора глюкозы, так и присутствие шума в его сигнале могут характеризоваться различными уровнями амплитуды сигнала и/или продолжительностью. Кроме того, неправильная работа сенсора глюкозы и присутствие шума в его сигнале могут перекрываться во времени, частично или полностью. Скрытая марковская модель может быть настроена на определение вероятности PA 758 в любых из этих условий. Как это обсуждается ниже, вероятность PA 758 может использоваться в рекурсивном фильтре для минимизации влияния, оказываемого неправильной работой сенсора глюкозы, что позволяет точно оценивать фактический уровень глюкозы человека в случае неправильной работы сенсора глюкозы и/или присутствия шума в его сигнале.
Возвращаясь к фиг. 7, система 750 для оценивания уровня глюкозы у больного диабетом может включать в себя рекурсивный фильтр 752, который может использоваться для оценивания уровня глюкозы человека на основании множества измеренных значений глюкозы, взвешенных по вероятности PA 758 точности сенсора глюкозы. К примерам используемых рекурсивных фильтров относятся фильтр Калмана и расширенный фильтр Калмана (РФК). Разумеется, могут использоваться и многие другие типы рекурсивных фильтров.
В одном варианте рекурсивным фильтром 752 может быть фильтр Калмана (в дальнейшем понятие "фильтр Калмана" также употребляется применительно к "расширенному фильтру Калмана"), настроенный на обработку измеренных значений 640 глюкозы (т.е. необработанных данных от сенсора глюкозы) в линейной системе второго порядка, описываемой приведенными ниже уравнениями. Фильтр Калмана может включать в себя, среди прочего, вектор состояния, представляющий оценочное состояние оцениваемой переменной, которой в данном примере является глюкоза человека. Фильтр Калмана может включать шаг прогнозирования, на котором прогнозируются априорное состояние и ковариация, а также шаг измерения, на котором обновляются апостериорный коэффициент (Kk) усиления фильтра Калмана, вектор состояния и ковариация. Вектор состояния может обновляться при каждом новом поступлении информации (т.е. рекурсивно). В настоящем изобретении переменные в векторе состояния х могут представлять оценку фактического уровня глюкозы человека, основанную на измеренных значениях 640 глюкозы. Оценочный вектор x уровня глюкозы может представлять оценочный уровень g глюкозы человека, его первую производную
Figure 00000027
и его вторую производную
Figure 00000028
. Вектор z измеренных значений глюкозы может включать в себя значения текущей глюкозы CG и их значения изменений глюкозы ΔCG. Могут использоваться и другие динамические модели. Векторы х и z можно представить как
Figure 00000029
и zk=[CG ΔCG]T, где k означает k-ую выборку. Для оценивания вектора х уровня глюкозы может использоваться следующее уравнение:
Figure 00000030
, где k означает k-ую выборку,
Figure 00000031
, Kk - коэффициент усиления фильтра Калмана и PA 758 - вероятность точности сенсора глюкозы (от средства вероятностного анализа). Таким образом, вероятность PA 758 точности сенсора может использоваться для взвешивания измеренных значений глюкозы, реализованного в матрице zk. Матрицы и сопутствующие уравнения для фильтра Калмана могут быть следующими:
Figure 00000032
,
Figure 00000033
,
Figure 00000034
,
Figure 00000035
,
Figure 00000036
,
Figure 00000037
,
Figure 00000038
,
Figure 00000039
, и
Figure 00000040
.
Параметры β1 и β2 в матрице А могут быть заданы незначительно меньшими единицы, чтобы в случае неправильной работы сенсора оценочный уровень глюкозы затухал. Матрица Q может представлять ковариацию шума процесса, а Kk может представлять коэффициент усиления фильтра Калмана. Первоначальные или предварительные оценки для этих параметров можно определять известным из уровня техники образом.
В расширенном фильтре Калмана (РФК) система может быть представлена нелинейной моделью,
Figure 00000041
, а измерения также представляются нелинейной моделью, xk=h(xk). Эта нелинейная модель может включать информацию, поступающую от других источников, uk, к которой могут относиться сведения о приеме пищи, введении инсулина, физических упражнениях или другие входные сведения, которые могут оказывать влияние на модель поведения глюкозы. Нелинейная модель может быть построена на основе частных физиологических моделей поведения глюкозы. Шаг прогнозирования выполняется путем оценка нелинейной модели, а прогнозная неопределенность вычисляется с использованием якобиана модели, Fk, при помощи вектора состояния. Это линейную модель, локализованную вокруг текущего состояния системы. Расширенным фильтром Калмана могут использоваться следующие уравнения:
Figure 00000042
Figure 00000043
Figure 00000044
Figure 00000045
Figure 00000046
Figure 00000047
Figure 00000048
.
После шага прогнозирования текущий результат измерения глюкозы соответствующим сенсором, CGk, может использоваться на шаге коррекции. На шаге коррекции также может использоваться ранее вычисленная вероятность PA 758 точности сенсора глюкозы (полученная от средства вероятностного анализа). Фильтр Калмана может быть настроен для взвешивания текущего измеренного значения глюкозы по вероятности точности сенсора глюкозы. Например, когда вероятность PA 758 низка, влияние текущего измеренного значения глюкозы на фильтр Калмана может стремиться к нулю; и наоборот, когда вероятность PA 758 высока, влияние текущего измеренного значения глюкозы может быть выше.
Проведение различия между неправильной работой сенсора и шумом сенсора может способствовать оцениванию уровня глюкозы человека, и в этом качестве фильтр Калмана может обращаться с этими ситуации по-разному. Для шума сенсора с нормальным распределением фильтр Калмана может быть настроен на исключение такого шума путем усреднения. Это может быть связано с тем, что шум сенсора может быть охарактеризован для каждого типа и/или для каждой партии сенсоров глюкозы, включая, кроме прочего, частотный диапазон шума и соответствующий размах изменений амплитуды измеренных значений глюкозы. Эти характеристики шума могут быть реализованы в некоторых или во всех параметрах фильтра Калмана (например, в параметрах
Figure 00000049
или
Figure 00000050
) таким образом, чтобы настроить фильтр Калмана на отфильтровывание шума и выдачу относительно точного оценочного уровня глюкозы человека даже в присутствии шума сенсора. С другой стороны, ошибке, порождаемой неправильной работой сенсора, нормальное распределение обычно не свойственно, вследствие чего такая ошибка обрабатывается в рамках фильтра Калмана иначе. В одном варианте реализации фильтра Калмана им может использоваться для взвешивания измеренных значений глюкозы вероятность PA 758 (определенная средством вероятностного анализа), что позволяет в основном не учитывать измеренные значения глюкозы при неправильной работе сенсора.
Ковариационная матрица из фильтра Калмана представляет неопределенность сигнала, но не включает неопределенность, вызванную калибровкой. Обусловленная калибровкой неопределенность должна начинаться с высокого значения, а затем уменьшаться на основании числа измерений глюкозы, а также разности в парах значений глюкозы и откалиброванных результатов измерений в рамках НМГ. Для учета этой неопределенности разность между значением глюкозы и результатом измерений в рамках НМГ должна принимать во внимание ковариацию фильтра Калмана. При такой модификации ошибка в сигнале в большинстве случаев не будет влиять на неопределенность, относимую на счет калибровки. Соответственно, калькулятор болюсов по меньшей мере в одном варианте его реализации можно регулировать, используя скорость изменения значений глюкозы, а также вычисление, основанное на оценке рисков.
Пример работы фильтра Калмана представлен на фиг. 9, где отражены измеренные значения глюкозы 970 и оценочный уровень глюкозы 978 человека. Как было указано выше, на этом графике также показано поведение вероятности точности сенсора глюкозы, PA 758. В нормальных условиях оценочный уровень глюкозы 978 человека может повторять измеренные значения глюкозы 970. Однако в период 974 сенсор может функционировать неправильно; одновременно вероятность PA 758 может уменьшиться до около 0% (определенная путем использования средства вероятностного анализа), указывая на низкую вероятность точности сенсора глюкозы. Соответственно, фильтр Калмана может учесть вероятность PA 758, уменьшив значимость измеренных значений глюкозы при оценивании уровня глюкозы человека в период 974 неправильной работы сенсора.
Далее, как показано на фиг. 9, в период 976 измеренные значения глюкозы 970 могут содержать шум. Фильтр Калмана может отфильтровывать этот шум, чтобы в указанный период 976 выдаваемый им оценочный уровень глюкозы 978 был относительно гладким. Хотя в период 976 измеренные значения глюкозы могут содержать шум, вероятность PA 758 в это время может оставаться относительно высокой (например, может оставаться вблизи 100%), так как средство вероятностного анализа может быть способно проводить различие между шумом сенсора и его неправильной работой. По этой причине фильтр Калмана может продолжать придавать значениям глюкозы, измеренным в период 972, относительно высокую важность (как это следует из относительно высоких значений вероятности PA 758 в период 972).
Неопределенность в результатах измерения глюкозы сенсором, как правило, не является постоянной. Ее текущее значение можно оценивать как функцию последних показаний z сенсора, вероятности Ра точности сенсора глюкозы, максимальной неопределенность
Figure 00000051
результата измерения и нормальной неопределенности
Figure 00000052
, связанной с непрерывными измерениями глюкозы. Значение σmax можно вычислять как максимальный физиологический разброс глюкозы у человека с плохо контролируемым диабетом. Его можно оценивать по выборкам данных НМГ. Аналогично, σCG - это минимальная неопределенность для сенсора глюкозы при его нормальной работе. Она может соответствовать работе сенсора в наилучшем случае и может оцениваться по вариации измеренных значений глюкозы по сравнению с данными, полученными путем взятия пробы крови из пальца, когда сенсор работает идеально. Могут существовать и другие методы оценивания неопределенности измерений, включающие использование повышенной частоты сбора данных от сенсора глюкозы. Это можно интерпретировать как вариацию разности между недавними результатами измерения CG (текущего уровня глюкозы) и оцененным состоянием фильтра Калмана.
Оценочный уровень глюкозы человека, определяемый рекурсивным фильтром, может использоваться для прогноза уровня глюкозы человека в определенный момент времени в будущем. Эти оценки также могут использоваться для анализа поведения человека и закономерностей в поведении глюкозы. Возвращаясь к фиг. 7, для прогнозирования возможности и/или времени наступления у человека гипогликемии может использоваться прогнозирующий алгоритм 762, который может инициировать выдачу тревожных сообщений или предупреждений. На вход прогнозирующего алгоритма 762 может подаваться оценочный уровень глюкозы 760 человека из рекурсивного фильтра 752, а также может подаваться неопределенность этого оценочного уровня глюкозы. Вместе с тем, прогнозирующий алгоритм 762 дополнительно может учитывать другие входные данные, в том числе времена приема пищи, количества углеводов, прием лекарственных средств, физические упражнения, дозы инсулина и так далее. На вход прогнозирующего алгоритма 762 также может поступать информация от других источников данных, например, измеренные значения глюкозы (т.е. необработанные данные от сенсора глюкозы) или обработанные данные от сенсора глюкозы. Прогнозирующий алгоритм 762 может использовать регрессию на основе гауссовского процесса для выучивания специфической для конкретного пациента модели прогноза, представленной на фиг. 7 обучающей моделью 764. Прогнозирующий алгоритм 762 также может оценивать неопределенность прогноза, что может допускать возможность коррекции порогов выдачи тревожных сообщений по чувствительности. Пороги выдачи тревожных сообщений также могут корректироваться на основании текущей деятельности пациента; например, чувствительность может быть повышена, когда пациент спит.
В качестве примера, прогноз гипогликемии может выполняться с использованием системной модели фильтра Калмана или расширенного фильтра Калмана. В данном примере шаг прогнозирования,
Figure 00000053
или
Figure 00000054
, неоднократно повторяют в течение заданного времени прогноза, и прогнозное значение сравнивают с конкретным порогом. Например, если фильтр Калмана обновляется каждую минуту, на шаге прогнозирования фильтр Калмана может задействоваться сорок пять раз для прогнозирования того, как уровень глюкозы человека изменится от настоящего времени до момента времени через сорок пять минут в будущем. Модель прогноза может включать дополнительные прогнозные входные данные, такие как предполагаемые приемы пищи, введение инсулина, физические упражнения или входную информацию о прочих ожидаемых в будущем факторах.
В еще одном примере оценочное значение глюкозы g и скорость g его изменения, оцененные рекурсивным фильтром, используются для задания линейного прогноза, который сравнивается с порогом гипогликемии и/или порогом гипергликемии. Прогноз выполняется при помощи следующего уравнения, предусматривающего умножение производной на требуемое время прогноза tpt, и вычисление прогнозного значения глюкозы
Figure 00000055
:
Figure 00000056
.
В качестве примера, конкретные используемые входные векторы могут включать три выборки оценочного уровня глюкозы (CG), взятых в моменты времени t=0, -15, и -30 минут, текущую производную оценочного уровня глюкозы и производную в момент времени t=-15 минут, а также время, прошедшей после последнего приема пищи. Эта информация о приеме пищи, tmeal, а также информация о введении болюса В является факультативной, при этом могут быть включены и другие данные. Математически это можно выразить следующим образом:
Figure 00000057
Figure 00000058
Регрессия на основе гауссовского процесса может использовать следующее уравнение для прогноза будущих уровней глюкозы человека на основании данных для обучения, представленных (Х, у), и точки измерения (х*, у*):
Figure 00000059
,
где k(х, х) - функция ковариации. Для генерирования результатов может использоваться гауссовская функция ковариации, но также здесь могут использоваться и другие функции. Используемая гауссовская функция ковариации имеет вид:
Figure 00000060
.
На фиг. 10 представлено функционирование прогнозирующего алгоритма. Слева (от момента времени t=-40 до 0) показаны измеренные значения глюкозы 1080 от сенсора глюкозы и оценочный уровень глюкозы 1082 человека (т.е. данные на выходе фильтра Калмана). Текущим временем является t=0. Прогнозирующий алгоритм может определять прогнозный уровень глюкозы 1084 человека в определенный момент времени в будущем (т.е. в любой момент времени, больший t=0). Кроме того, прогнозирующий алгоритм может использоваться для прогнозирования возможности и/или времени наступления у человека гипогликемии или гипергликемии, т.е. для прогнозирования того, станет ли уровень глюкозы человека гипо- или гипергликемическим, и/или того, когда это случится. Для человека может быть установлен порог гипогликемии 1086, причем снижение фактического уровня глюкозы с переходом через этот порог означает наступление у человека гипогликемического состояния. Порог гипогликемии 1086 может определяться как уникальный для каждого человека. Для среднего человека этот порог может составлять около 50 мг/мл. Также для каждого человека порог гипогликемии 1086 может колебаться, например, в зависимости от времени, того или иного события или сочетания таких факторов. В качестве примеров, порог гипогликемии 1086 для человека может зависеть от времени суток, от приема человеком лекарственных средств, от того, находится ли сенсор в состоянии временного неиспользования и/или как долго он находится в этом состоянии, и т.д. На фиг. 10 прогнозирующий алгоритм предсказывает, что человек окажется в состоянии гипогликемии в момент времени t=45 (т.е. через 45 минут от текущего момента). Разумеется, с течением времени прогнозирующий алгоритм может продолжать использование самой последней оценки уровня глюкозы (получаемой от фильтра Калмана) и соответственно корректировать прогнозируемые уровни глюкозы.
Помимо того, что прогнозирующий алгоритм может предсказывать значения уровня глюкозы человека в будущем, он также может быть настроен для определения вероятности того, что такой прогноз является точным. Например, прогнозы лишь на одну минуту или на две минуты вперед могут иметь высокую точность, а прогнозы на 60 или 70 минут вперед могут быть относительно неточными. Разумеется, вероятность того, что прогноз точен, может представлять собой непрерывное множество значений, начинающихся вблизи 100% для ближайшего будущего и постепенно уменьшающихся до около 0% по мере увеличения срока прогноза. Эта информация (о точности прогноза) может использоваться, в сочетании с собственно прогнозом, для создания системы предупреждения человека о риске наступления гипогликемии. Как показано на фиг. 10, такая система предупреждения может не выдавать тревожных сообщений 1088, если прогнозный уровень глюкозы 1084 достаточно высок относительно порога гипогликемии 1086; она может выдать предостережение 1090 при приближении прогнозного уровня глюкозы 1084 к порогу гипогликемии 1086 на заданную величину; и она может сообщать об опасности 1092 при падении прогнозного уровня глюкозы 1084 ниже порога гипогликемии 1086.
Прогнозирующий алгоритм, как уже было указано, может включать функцию обучения, способную обучаться особенностям конкретного человека. Функция обучения может создавать обучающие данные, которые могут использоваться в прогнозирующем алгоритме и могут взвешиваться по влиянию, оказываемому ими на формирование прогноза. Уровень влияния обучающих данных может определяться функцией ковариации k(х,х), используемой внутри независимой переменной (регрессора) в регрессии на основе гауссовского процесса.
Прогнозирующий алгоритм может инициализироваться с помощью генерического множества обучающих примеров или без обучающих примеров. По мере получения результатов измерений в прогнозирующий алгоритм и/или функцию обучения могут включаться новые измеренные данные. Есть много возможных алгоритмов для включения новых данных. Такие алгоритмы предусматривают добавление данных в обучающее множество, когда: 1) истек заданный период времени, 2) прогноз на основе конкретных данных оказался неудачным, 3) входные данные в обучающем множестве не представлены или 4) пациент или лицо, осуществляющее уход за ним, вручную включает данные, в том числе все новые данные, если это уместно.
При добавлении новых данных в обучающее множество новые данные могут включаться в виде нового вектора или путем перевзвешивания существующего обучающего вектора. Второй метод выгоден тем, что позволяет сохранить неизменными требования к объему памяти. После добавления дополнительных данных прогнозирующий алгоритм может быть обновлен на устройстве сразу, на персональном компьютере - ретроспективно или также ретроспективно в клинике.
На фиг. 11 представлен способ 1100 оценивания уровня глюкозы у больного диабетом. Способ может включать ряд действий, которые могут выполняться в любом подходящем порядке. При совершении действия 1102 способ 1100 может предусматривать получение вычислительным устройством множества измеренных значений глюкозы от сенсора глюкозы, подключенного к организму человека. При совершении действия 1104 способ 1100 может предусматривать использование вычислительного устройства для анализа множества измеренных значений глюкозы с использованием средства вероятностного анализа, способного определять вероятность точности сенсора глюкозы на основании указанного множества результатов измерения глюкозы. При совершении действия 1106 способ 1100 предусматривать использование вычислительного устройства для оценивания уровня глюкозы человека с использованием рекурсивного фильтра, способного взвешивать множество измеренных значений глюкозы по вероятности точности сенсора глюкозы. Средство вероятностного анализа и рекурсивный фильтр могут иметь описанное выше исполнение.
Калькулятор болюсов инсулина
На фиг. 12 приведена блок-схема системы 1210 рекомендации болюсов инсулина в одном иллюстративном варианте ее выполнения. В рассматриваемом варианте система 1210 рекомендации болюсов инсулина включает в себя блок 1212 рекомендации болюсов, имеющий по меньшей мере схему 1214 управления, электрически соединенную с устройством 1216 визуального отображения, а также с устройством 1218 ввода информации. Схема 1214 управления может быть, например, обычным микропроцессорным управляющим компьютером, способным выполнять один или несколько программных алгоритмов, хотя в альтернативном исполнении схема 1214 управления может представлять собой любую одиночную электронную схему или набор электронных схем, способных работать, как это описано ниже. В некоторых вариантах осуществления изобретения схема 1214 управления может быть электрически подключена к обычному запоминающему устройству 1220, показанному штриховой линией. Устройство 1216 визуального отображения может представлять собой или включать в себя любой обычный визуальный индикатор, в том числе дисплей на электронно-лучевой трубке (ЭЛТ), жидкокристаллический дисплей, газоразрядный (плазменный) дисплей, одно- или многоцветный монитор, сенсорный экран с возможностью ввода данных прикосновением к нему и т.п. Устройство 1218 ввода информации может представлять собой или включать в себя любое обычное устройство ввода, в том числе клавиатуру, вспомогательную клавиатуру, мышь или аналогичное координатно-указательное устройство, один или несколько кодированных или некодированных, чувствительных к прикосновению коммутаторов (ключей), связанных с устройством 1216 отображения, управляемое голосом устройство ввода данных или аналогичные устройства.
В некоторых вариантах осуществления изобретения система 1210 рекомендации болюсов инсулина также может включать в себя дополнительный блок 1230 рекомендации болюсов, показанный на фиг. 12 штриховой линией. Блок 1230 может содержать схему 1232 управления, электрически соединенную с устройством 1234 визуального отображения, а также с устройством 1236 ввода информации, причем схема 1232 управления, устройство 1234 отображения и устройство 1236 ввода информации могут быть выполнены в любой из форм, описанных выше в отношении блока 1212 рекомендации болюсов. Схема 1232 управления также может быть электрически подключена к обычному запоминающему устройству 1238. В этом варианте осуществления изобретения каждый из блоков 1212, 1230 рекомендации болюсов может быть выполнен с возможностью обмена информацией посредством проводного соединения 1240, включающего одну или несколько линий передачи сигналов, физически соединяющих два блока, посредством линии 1242 беспроводной передачи сигналов, такой как канал передачи радиосигналов или канал сотовой телефонной связи, и/или посредством всемирной сети (WWW) 1244, причем в каждом случае могут использоваться обычные технологии.
Система 1210 рекомендации болюсов инсулина выполнена с возможностью определения одного или нескольких конкретных количеств болюсного инсулина и рекомендации введения этих болюсов в кровоток пользователя системы 1210 согласно протоколу рекомендации болюсов инсулина, реализованному в системе 1210 в виде одного или нескольких исполняемые программные алгоритмы. Физическая структура системы 1210 рекомендации болюсов инсулина для выполнения таких программных алгоритмов и для передачи полезной информации между системой 1210 и пользователем может принимать различные формы. В одном иллюстративном варианте осуществления изобретения система 1210 рекомендации болюсов включает в себя, например, только блок 1212 рекомендации болюсов, реализованный в обычном персональном компьютере (ПК), переносном компьютере или ноутбуке, персональном цифровом ассистенте (карманном компьютере) или в ином аналогичном устройстве, либо выполненный в виде портативного, переносного или настольного специализированного устройство. В любом из этих случаев блок 1212 рекомендации болюсов содержит запоминающее устройство 1220, в котором хранится несколько исполняемых программных алгоритмов, и схему 1214 управления, обеспечивающую выполнение этих программных алгоритмов для определения и рекомендации одной или нескольких инъекций конкретных количеств болюсного инсулина в кровоток пользователя согласно протоколу рекомендации болюсов инсулина, как это подробно описывается ниже. В этом варианте осуществления изобретения устройство 1216 отображения (дисплей) управляется схемой 1214 управления в соответствии с программными алгоритмами для вывода пользователю информации и выдачи пользователю предложения ввести информацию, которую пользователь может ввести посредством устройства 1218 ввода информации.
В еще одном иллюстративном варианте осуществления изобретения система 1210 рекомендации болюсов инсулина включает в себя один блок 1212 рекомендации болюсов и другой блок 1230 рекомендации болюсов. В качестве примера такого варианта осуществления изобретения, блок 1212 рекомендации болюсов может быть карманным компьютером или специализированным блоком рекомендации болюсов, как это упоминалось выше, а блок 1230 рекомендации болюсов может быть персональным компьютером (ПК), переносным компьютером или ноутбуком. В этом варианте осуществления изобретения блок 1212 может быть связан с блоком 1230, т.е. может взаимодействовать с ним, посредством либо беспроводного интерфейса 1242, либо проводного интерфейса 1240, который может быть электрически соединен с карманным компьютером или с подставкой (док-станцией) под специализированный блок рекомендации болюсов, способной удерживать блок 1212 и электрически соединять его с блоком 1230 с возможностью обмена данными. В данном примере в запоминающем устройстве 1220 блока 1212 и в запоминающем устройстве 1238 блока 1230 может храниться ряд программных алгоритмов, и пользователь может использовать блок 1212 рекомендации болюсов как мобильный блок рекомендации болюсов инсулина и/или может использовать блок 1230 рекомендации болюсов как стационарный блок рекомендации болюсов инсулина. В данном случае пользователь будет поддерживать базы данных каждого из блоков 1212 и 1230 в актуальном состоянии, периодически синхронизируя базы данных обоих блоков 1212 и 1230 посредством проводного 1240 или беспроводного 1242 интерфейса.
В еще одном примере варианта системы 1210 рекомендации болюсов инсулина, в котором эта система включает в себя блоки 1212 и 1230 рекомендации болюсов, блок 1212 рекомендации болюсов может быть карманным компьютером, персональным компьютером, переносным компьютером или ноутбуком, сотовым телефоном или любым другим узлом или устройством, имеющим возможность доступа во всемирную сеть 1244. В данном примере в запоминающем устройстве 1220 блока 1212 рекомендации болюсов не обязательно должны храниться программные алгоритмы, и более того, блоку 1212 рекомендации болюсов вообще не нужно иметь запоминающего устройства 1220. Блок 1230 рекомендации болюсов, в данном примере, может быть удаленным компьютером или обычным веб-сервером, также выполненным с возможностью доступа во всемирную сеть 1244 и имеющим запоминающее устройство 1238, в котором хранятся программные алгоритмы. Схема 1232 управления удаленного компьютера или веб-сервера 1230 в данном примере обеспечивает выполнение программных алгоритмов на основании информации, предоставляемой через всемирную сеть 1244 пользователем посредством блока 1212 рекомендации болюсов. В этом частном варианте осуществления изобретения пользователь и/или сотрудник медицинского учреждения может получить доступ к веб-странице или веб-сайту, управляемой(-ому) блоком 1230 рекомендации болюсов и внести в схему 1232 управления исходные рабочие параметры и/или пределы для протокола рекомендации болюсов инсулина. Затем пользователь может войти на веб-страницу или веб-сайт и ввести информацию о текущем уровне глюкозы крови, и тогда схема 1232 управления на основании информации о текущем уровне глюкозы крови может определить и рекомендовать, посредством веб-страницы или веб-сайта, одну или несколько инъекций конкретных количеств болюсного инсулина в кровоток пользователя в соответствии с протоколом рекомендации болюсов инсулина, подробно рассматриваемом ниже.
Таким образом, в этом частном варианте осуществления изобретения программные алгоритмы рекомендации болюсов инсулина находятся в удаленном компьютере или веб-сервере 1230, и поэтому блоку 1212 рекомендации болюсов нужно иметь аппаратные средства, достаточные лишь для того, что обеспечить возможность предоставления веб-странице или вебсайту информации о текущем уровне глюкозы крови и просмотра результатов с рекомендацией, выдаваемых на веб-странице или веб-сайте удаленным компьютером или веб-сервером 1230. Хотя на практике в этом варианте осуществления изобретения может быть желательно снабдить блок 1212 рекомендации болюсов запоминающим устройством 1220 и хранить в этом запоминающем устройстве набор программных алгоритмов рекомендации болюсов, чтобы блок 1212 рекомендации болюсов мог выполнять эти программные алгоритмы независимо, если, например, доступ во всемирную сеть 1244 и/или к соответствующим веб-странице или веб-сайту невозможен или практически неосуществим. Кроме того, в этом варианте осуществления изобретения желательно предусмотреть синхронизацию удаленной и/или доступной через сеть базы данных с базой данных, хранящейся в запоминающем устройстве 1220 блока 1212 рекомендации болюсов.
Понятно, что система 1210 рекомендации болюсов инсулина может быть выполнена с возможностью взаимодействия с глюкометром либо иным устройством для автоматического определения глюкозы крови и/или с инсулиновой помпой либо иным устройством для автоматического дозирования или введения инсулина. В вариантах осуществления изобретения, в которых глюкометр или иное устройство для автоматического определения глюкозы крови входит в состав системы 1210 рекомендации болюсов инсулина, управляющий компьютер 1214 может быть выполнен с возможностью управления таким устройством, с использованием обычных методов, таким образом, чтобы оно автоматически выдавало информацию о текущем уровне глюкозы крови, которую система 1210 может затем использовать, как это подробно описывается ниже, для определения и рекомендации введения одного или нескольких количеств болюсного инсулина. В вариантах осуществления изобретения, в которых инсулиновая помпа или другое устройство для автоматического дозирования инсулина входит в состав системы 1210 рекомендации болюсов инсулина, управляющий компьютер 1214 может быть выполнен с возможностью управления таким устройством, с использованием обычных методов, таким образом, чтобы оно автоматически вводило в организм пользователя рекомендованные количества болюсного инсулина.
Как описано выше, представленная на фиг. 12 система 1210 рекомендации болюсов инсулина обеспечивает выполнение ряда программных алгоритмов для определения и рекомендации введения в кровоток пользователя одного или нескольких конкретных количеств болюсного инсулина согласно протоколу рекомендации болюсов инсулина. Протокол введения болюсов инсулина, в рамках осуществления настоящего изобретения, предусматривает разбиение длительного периода времени, например, одних или нескольких суток, на несколько соседних временных интервалов. Для любого такого временного интервала можно задать определенный целевой уровень глюкозы и определенное значение чувствительности к инсулину. Целевой уровень глюкозы для любого временного интервала соответствует некоторому постоянному значению глюкозы, которое пользователь хотел бы поддерживать в течение этого временного интервала. В качестве примера целевого уровня глюкозы можно указать значение 120 мг/дл, хотя различные целевые уровни глюкозы, задаваемые для некоторых временных интервалов, могут принимать иные значения. Значение чувствительности к инсулину для любого временного интервала соответствует снижению глюкозы, выраженному, например, в мг/дл на единицу вливаемого инсулина. В качестве примера показателя чувствительности к инсулину можно указать значение 30 мг/дл/ед., хотя в других случаях могут использоваться и другие значения чувствительности к инсулину. Подробные сведения, относящиеся к примерам программных алгоритмов, которые могут выполняться системой 1210 для осуществления такого протокола рекомендации болюсов инсулина, описаны в патентной заявке США №10/927614, которая принадлежит заявителю и содержание которой включено в данное описание путем ссылки. Вместе с тем, следует иметь в виду, что для осуществления такого протокола рекомендации болюсов инсулина система 1210, в качестве альтернативы вышеуказанным алгоритмам или в дополнение к ним, также может быть запрограммирована для выполнения других обычных программных алгоритмов.
На фиг. 13 представлен график изменения глюкозы крови во времени, на котором показано два временных интервала, заданных в соответствии с описанным выше протоколом рекомендации болюсов инсулина. Текущий временной интервал начинается в момент времени TBCI (англ. сокр. от "Begin Current Interval" - начало текущего интервала) и заканчивается в последующий момент времени TECI (англ. сокр. от "End Current Interval" - окончание текущего интервала). Целевое значение 1350 глюкозы крови, заданное в текущем интервале, обозначено как BGTCI, а чувствительность к инсулину, заданная в текущем интервале, обозначена как ISCI. Следующий временной интервал, т.е. временной интервал, непосредственно наступающий за текущим, начинается в момент времени TBNI (англ. сокр. от "Begin Next Интервал" - начало следующего интервала), совпадающий с моментом TECI окончания текущего интервала, и заканчивается в последующий момент времени TENI (англ. сокр. от "End Next Interval" - окончание следующего интервала). Целевое значение 1352 глюкозы крови, заданное в следующем соседнем интервале, обозначено как BGTNI, а чувствительность к инсулину, заданная в следующем соседнем интервале, обозначена как ISNI.
Пользователь системы 1210 может в любое время получить результат измерения своего уровня глюкозы, используя один или несколько обычных методов. Если результат измерения текущего уровня глюкозы пользователя превышает целевой уровень глюкозы крови для текущего интервала, BGTCI, при использовании обычных систем рекомендации болюсов, как правило, вычисляется рекомендуемый корректирующий болюс инсулина СВ, для чего используется уравнение:
Figure 00000061
где BGM - уровень глюкозы крови, измеренный в текущий момент времени, TC. В этом случае пользователь вблизи момента времени TC вводит рекомендованный корректирующий болюс инсулина, который действует известным образом, снижая уровень глюкозы в течение периода времени TM, где период TM в контексте данного описания задан как длительность сахароснижающего действия введенного болюса инсулина.
На фиг. 13 описанный выше сценарий представлен в текущем интервале, когда в текущий момент времени TC было измерено значение 1354 глюкозы крови. Благодаря тому, что вблизи момента времени TC введен корректирующий болюс СВ, вычисленный как было описано чуть выше, глюкоза крови в данном примере снижается с течением времени (по кривой 1356) до целевого значения 1350 глюкозы, которое соответствует целевому значению глюкозы крови для текущего интервала, BGTCI, и остается на уровне BGTCI в течение времени TM. В качестве численного примера допустим, что измеренное значение глюкозы крови, BGM, соответствующее точке 1354, равно 200 мг/дл, целевое значение глюкозы крови в текущем интервале, BGTCI, равно 120 мг/дл, а чувствительность к инсулину в текущем интервале, ISCI, равна 30 мг/дл/ед. Подстановка этих значений в приведенное выше уравнение для определения корректирующего болюса обычным методом дает (200 мг/дл - 120 мг/дл)/30 мг/дл/ед.=2,667, или приблизительно 2,7 единицы. Таким образом, для снижения у пользователя уровня глюкозы с 200 мг/дл до целевых 120 мг/дл в течение текущего интервала требуется приблизительно 2,7 единицы инсулина.
Как показано в приведенном выше примере, обычное уравнение для определения корректирующего болюса действует удовлетворительно при условии, что корректирующий болюс вводится в текущем временном интервале достаточно быстро (рано), вследствие чего длительность TM сахароснижающего действия введенного болюса инсулина, ограничивается текущим интервалом. Однако использование обычного уравнения для определения корректирующего болюса как основы для вычисления и введения корректирующего болюса в текущий момент времени TC, более поздний в текущем временном интервале, вследствие чего длительность TM сахароснижающего действия введенного болюса перекрывает текущий и следующий за ним временные интервалы, может дать нежелательные результаты. Рассмотрим, например, случай, когда значение глюкозы крови 1358 измеряется в другой текущий момент времени TC, который находится близко к моменту TECI окончания текущего временного интервала (соответствующего моменту TBNI начала следующего временного интервала). При вычислении корректирующего болюса СВ при помощи описанного выше обычного уравнения и при введении этого болюса вблизи момента времени TC глюкоза крови с течением времени (по кривой 1360) снижается, в том числе, в данном примере, в следующем временном интервале, до уровня глюкозы 1362, который может быть значительно ниже целевого значения глюкозы 1352 в этом следующем временном интервале. В качестве еще одного численного примера допустим, что измеренное значение глюкозы крови, BGM, соответствующее точке 1358, также равно 200 мг/дл, целевое значение глюкозы крови в текущем интервале, BGTCI, также равно 120 мг/дл, чувствительность к инсулину в текущем интервале, ISCI, также равна 30 мг/дл/ед., целевое значение глюкозы крови в следующем временном интервале, BGTNI, равно 150 мг/дл, а чувствительность к инсулину в следующем временном интервале, ISNI, равна 40 мг/дл/ед. Поскольку момент времени TC еще находится в текущем временном интервале, приведенное выше обычное уравнение для определения корректирующего болюса снова дает (200 мг/дл - 20 мг/дл)/30 мг/дл/ед.=2,667, т.е. приблизительно 2,7 единицы. При этом допустим, что в данном примере момент времени TC находится в пределах одной минуты от момента TBNI начала следующего временного интервала. Если бы пользователь подождал еще минуту до проведения измерения глюкозы крови, вследствие чего следующий временной интервал стал бы текущим временным интервалом, приведенное выше обычное уравнение для определения корректирующего болюса дало бы (200 мг/дл - 150 мг/дл)/40 мг/дл/ед.=1,25, т.е. приблизительно 1,3 единицы. Таким образом, в данном примере введение болюса инсулина в количестве 2,7 единицы в момент времени TC или вблизи него приводит к поступлению в организм избыточного инсулина в количестве 1,4 единицы, что при чувствительности к инсулину, ISNI, составляющей 40 мг/дл/ед., приводит к чрезмерному снижению глюкозы крови до уровня 1362, который на (40 мг/дл/ед. × 1,4 ед.)=56 мг/дл ниже целевого уровня 1352 глюкозы крови, что в конечном итоге соответствует снижению глюкозы с 200 мг/дл в момент времени 1358 до 94 мг/дл в момент времени 1362.
Один эффективный метод повышения точности определения корректирующего болюса в случаях, когда период времени TM сахароснижающего действия введенного болюса инсулина перекрывает текущий и следующий за ним временные интервалы, заключается в том, чтобы учесть характер поведения во времени целевого уровня глюкозы и чувствительности к инсулину. Например, если h(τ) - относительная величина активности инсулина, использованной для болюса инсулина (данного типа), введенного в момент времени τ=0. В момент времени τ≥TM имеем h(τ)=1. Тогда если IS(T) - чувствительность к инсулину как функция времени, то в текущий момент времени Т=TC единичный болюс инсулина имеет следующую импульсную характеристику:
Figure 00000062
.
где d(τ) производная h(τ) по времени.
Если D(T) - зависящая от времени функция скорости инфузии болюса инсулина, падение глюкозы в данный момент времени относительно BG-∞=BG(T→-∞) дает:
Figure 00000063
.
Чтобы вычислить в момент времени Т=TBNI падение глюкозы крови ΔBGpb, которое произойдет в будущем под действием уже введенных прошлых болюсов инсулина, используется следующее уравнение:
Figure 00000064
.
Для вычисления количества (величины) корректирующего болюса СВ, подлежащего введению в момент времени Т=TC+С при длительности инфузии болюса, пренебрежимо малой по сравнению с TM, используется следующее уравнение:
Figure 00000065
.
Таким образом, потребное падение глюкозы крови (ГК - англ. сокр. BG от Blood Glucose), ΔBGneeded, нужно вычислять, исходя из целевого значения глюкозы крови в момент времени Т=TBNI+TM, или BGTNI. При объединении результатов выражений (4), (5) и (8) уравнение для вычисления количества корректирующего болюса СВ, когда длительность его действия перекрывает текущий и следующий за ним временные интервалы, имеет вид:
Figure 00000066
.
Решение уравнения (9) в отношении интегралов дает уравнение: CB=(BGM-BGTNI)/{[ISCI*(h(TBNI-TC)-h(0))]+[ISNI*(h(TM)-h(TBNI-TC))]} (10), где количество (h(TBNI-TC)-h(0)) соответствует той (процентной) доле действия болюсного инсулина, введенного в момент времени Т=TC или вблизи него, которая будет израсходована или использована на снижение измеренного значения глюкозы крови, BGM, в текущем временном интервале, а количество (h(TM)-h(TBNI-TC)) соответствует той (процентной) доле действия болюсного инсулина, введенного в момент времени Т=TC или вблизи него, которая будет потрачена или использована на снижение измеренного значения глюкозы крови, BGM, в следующем временном интервале.
Возвращаясь к фиг. 13, рассмотрим еще один численный пример, поясняющий влияние уравнения (10) на измеренное значение 1358 глюкозы крови. Допустим, что измеренное значение глюкозы крови, BGM, соответствующее точке 1358, также равно 200 мг/дл, целевое значение глюкозы крови в текущем интервале, BGTCI, также равно 120 мг/дл, чувствительность к инсулину в текущем интервале, ISCI, также равна 30 мг/дл/ед., целевое значение глюкозы крови в следующем временном интервале, BGTNI, также равно 150 мг/дл, а чувствительность к инсулину в следующем временном интервале, ISNI, также равна 40 мг/дл/ед. Допустим также, что к моменту времени Т=TBNI на снижение уровня глюкозы в текущем временном интервале будет потрачено 40% действия болюсного инсулина, введенного в момент времени Т=TC (соответствующий точке 1358) или вблизи него. Соответственно, на снижение уровня глюкозы в следующем временном интервале остается потратить 1-0,4=60% действия болюсного инсулина, введенного в момент времени Т=TC. Тогда согласно уравнению (10) корректирующий болюс СВ, который нужно ввести в момент времени Т=TC (соответствующий точке 1358) или вблизи него, рассчитывается следующим образом: СВ=(200 мг/дл - 50 мг/дл)/[(30 мг/дл/ед. × 0,4)+(40 мг/дл/ед. × 0,6)]=1,388 ед., или приблизительно 1,4 ед. При чувствительности к инсулину в следующем временном интервале 40 мг/дл/ед. введение 1,4 ед. инсулина в момент времени Т=TC (соответствующий точке 1358) приводит к уменьшению глюкозы крови по кривой 1364, составляющему (40 мг/дл/ед.) × 1,4 ед.=56 мг/дл и соответствующему снижению глюкозы крови с 200 мг/дл в момент времени 1358 до конечного уровня глюкозы 1366, составляющего приблизительно 144 мг/дл, что гораздо ближе к целевому значению глюкозы крови 150 мг/дл в следующем временном интервале, чем значение 94 мг/дл, получающееся при использовании обычного уравнения (4) для определения корректирующего болюса.
На фиг. 4 приведена блок-схема одного иллюстративного варианта программного алгоритма 1400 для определения количеств болюсного инсулина, рекомендуемых в различных условиях, представленных рассмотренными выше примерами. Показанный фиг.14 программный алгоритм 1400 рекомендации болюсов инсулина описывается как реализованный в блоке 1212 рекомендации болюсов и выполняемый схемой 1214 управления, причем блок 1212 рекомендации болюсов инсулина выполнен в виде обычного карманного компьютера или портативного специализированного блока рекомендации болюсов инсулина, хотя специалистам также понятна возможность альтернативной реализации алгоритма 1400 в блоке 1212 рекомендации болюсов и/или блоке 1230 рекомендации болюсов, выполненном(-ых) в одной или нескольких описанных выше физических формах.
В рассматриваемом варианте алгоритм 1400 начинается на шаге 1402, на котором схема 1214 управления приводится в действие для определения того, имеется ли новое измеренное значение глюкозы, BGM. Алгоритм 1400, в рассматриваемом примере осуществления изобретения, предполагает, что результаты измерений глюкозы, выполняемых в любой текущий момент времени, TC, вводятся в систему 1210, или поступают в нее иным образом, в момент времени TC или вблизи него. Пока не будет получено новое значение глюкозы, BGM, схема 1214 управления будет продолжать возвращаться к началу шага 1402. При появлении же нового результата измерения глюкозы, BGM, выполнение алгоритма переходит на шаг 1404, на котором схема управления приводится в действие для получения необходимых параметров, относящихся к текущему и следующему за ним временным интервалам. В одном варианте осуществления изобретения эти необходимые параметры хранятся в запоминающем устройстве 20 и/или в исполняемой памяти схемы 1214 управления, и схема 1214 управления приводится в действие для выполнения шага 1404 путем извлечения этих параметров из запоминающего устройства 1220 и/или из исполняемой памяти схемы 1214 управления. В качестве альтернативы, необходимые параметры могут вводиться в систему 1210 или предоставляться ей иным образом при помощи одного или нескольких компонентов, описанных выше со ссылкой на фиг. 12. В любом случае необходимые параметры, получаемые на шаге 1404 в рассматриваемом примере осуществления изобретения, включают текущее значение времени TC, значение времени начала следующего временного интервала, TBNI, длительность TM сахароснижающего действия болюса инсулина, возможно вводимого в организм в момент времени TC или вблизи него, измеренное значение глюкозы, BGM, целевое значение глюкозы для текущего временного интервала, BGTCI, целевое значение глюкозы для следующего временного интервала, BGTNI, чувствительность к инсулину для текущего временного интервала, ISCI, и чувствительность к инсулину для следующего временного интервала, ISNI.
После шага 1404 схема 1214 управления приводится в действие на шаге 1406 для сравнения измеренного значения глюкозы, BGM, с целевым значением глюкозы для текущего временного интервала, BGTCI. Если BGM не превышает BGTCI, выполнение алгоритма 1400 возвращается к началу шага 1402, в противном случае выполнение алгоритма переходит на шаг 1408. Таким образом, алгоритм 1400 не проходит через шаг 1406 до тех пор, пока не появится новый результат измерения глюкозы, BGM, и пока BGM не превысит BGTCI. Разумеется, что алгоритм 1400 может быть встроен в другой алгоритм рекомендации болюсов инсулина, способный выполнять шаги 1402 и 1406. В таких случаях алгоритм 1400 может быть модифицирован для его включения в такой алгоритм рекомендации болюсов инсулина за счет исключения шагов 1402 и 1406.
На шаге 1408 схема 1214 управления приводится в действие для сравнения суммы TC, TM и постоянного значения времени, TK, с временем начала следующего временного интервала, TBNI. Если эта сумма больше TBNI, выполнение алгоритма переходит на шаг 1410. В противном случае, если эта сумма меньше или равна TBNI, выполнение алгоритма переходит на шаг 1412. В одном варианте осуществления изобретения TK равно нулю, а сумма TC и TM представляет значение времени, по отношению к текущему и следующему за ним временным интервалам, когда действие болюсного инсулина, если таковой введен в текущий момент времени TC, окончится. Таким образом, если сумма TC и TM меньше или равна TBNI, это указывает, что момент времени TC находится в текущем временном интервале достаточно близко к началу этого интервала, чтобы длительность TM сахароснижающего действия болюса инсулина, введенного в момент времени TC, была ограничена текущим временным интервалом. В этом случае вычисление точного значения корректирующего болюса СВ можно выполнять при помощи обычного уравнения (4). Если же сумма TC и TM превышает TBNI, это указывает, что момент времени TC находится в текущем временном интервале достаточно близко к концу этого интервала, чтобы длительность TM сахароснижающего действия болюса инсулина, введенного в момент времени TC, перекрывала текущий и следующий за ним временные интервалы. В этом случае вычисление точного значения корректирующего болюса СВ нужно выполнять при помощи описанного выше уравнения (10) для определения корректирующего болюса. Возможны такие варианты реализации алгоритма 1400, в которых схема 1214 управления приводится в действие для вычисления корректирующего болюса СВ согласно уравнению (10), если только период TM сахароснижающего действия корректирующего болюса, подлежащего введению в момент времени TC, переходит в следующий временной интервал и продолжается в нем заданный отрезок времени. В таких вариантах осуществления изобретения постоянная времени TK не будет нулевой, а будет иметь некоторое положительное значение времени, благодаря которому удовлетворение используемому на шаге 1408 неравенству не будет приводить к переходу на шаг 1410, если период TM не перейдет в следующий временной интервал и не продолжится в нем заданный отрезок времени. Например, этот заданный отрезок времени может составлять 30% от периода TM, хотя, разумеется, заданный отрезок времени может принимать и другие значения.
В любом случае, если схема 1214 управления на шаге 1408 определит, что период TM сахароснижающего действия корректирующего болюса, подлежащего введению в момент времени TC, переходит в следующий временной интервал и продолжается в нем достаточное время, выполнение алгоритма переходит на шаг 1410, на котором схема управления приводится в действие для сравнения абсолютного значения разности между BGTCI и BGTNI с постоянной K1, а также для сравнения абсолютного значения разности между ISCI и ISNI с другой постоянной K2. В одном варианте осуществления изобретения K1 и K2 равны нулю, а абсолютные значения указанных разностей равны нулю, только если целевое значение глюкозы и значение чувствительности к инсулину между текущим и следующим за ним временными интервалами не изменяются. В этом случае вычисление точного значения корректирующего болюса СВ может выполняться с использованием обычного уравнения (4). Если же целевое значение глюкозы и/или значение чувствительности к инсулину между текущим и следующим за ним временными интервалами изменяется, вычисление точного значения корректирующего болюса СВ нужно выполнять при помощи описанного выше уравнения (10) для определения корректирующего болюса. Возможны варианты реализации алгоритма 1400, в которых постоянным K1 и K2 присваивают некоторые положительные постоянные значения, чтобы перед вычислением корректирующего болюса СВ согласно уравнению (10) целевое значение глюкозы и/или значение чувствительности к инсулину изменялись более чем на заданные величины. Например, K1 и K2 могут быть равны 5, хотя, разумеется, могут использоваться и другие значения K1 и K2, причем постоянная K1 может быть равна или неравна постоянной K2.
В любом случае, если период TM сахароснижающего действия болюса инсулина, подлежащий введению в момент времени TC, ограничен текущим временным интервалом или по меньшей мере не продолжается в нем достаточное время, и если ни целевые значения глюкозы, ни значения чувствительности к инсулину значительно не изменяются между текущим и следующим за нем временными интервалами, на шаге 1412 схема 1214 управления приводится в действие для вычисления корректирующего болюса СВ согласно уравнению (4), т.е. по формуле: СВ=(BGM-BGTCI)/ISCI. Если же период TM продолжается в следующем имеющемся временном интервале достаточное время, либо если целевое значение глюкозы и/или значение чувствительности к инсулину значительно изменяются между текущим и следующим имеющимся временными интервалами, схема 1214 управления приводится на шаге 1414 в действие для определения (процентной) доли PU действия болюсного инсулина, которая будет использована или потрачена на снижение уровня глюкозы крови в текущем временном интервале.
В рассматриваемом варианте схема 1214 управления приводится в действие для определения этой (процентной) доли как функции текущего времени TC, времени TBNI начала следующего временного интервала и длительности TM сахароснижающего действия корректирующего болюса, который будет введен в момент времени TC или вблизи него. Эта функция может храниться в запоминающем устройстве 1220 в форме одной или нескольких таблиц, графиков, диаграмм, уравнений и т.п., причем в одном конкретном варианте осуществления изобретения эта функция хранится в запоминающем устройстве 1220 в форме двумерной справочной таблицы (таблицы соответствия). В этом варианте осуществления изобретения по одной оси справочной таблицы в ней содержатся значения, соответствующие разности между значениями TBNI и значениями TC, а по другой оси справочной таблицы в ней содержатся значения длительности TM. В этом варианте осуществления изобретения таблица заполнена процентными значениями, соответствующими значениям процентной доли действия болюсного инсулина, которая будет использована или потрачена на снижение уровня глюкозы в текущем временном интервале, в зависимости от значений TM и значений разности времен TBNI-TC. После этого на шаге 1416 схема управления приводится в действие для определения (процентной) доли PR действия болюсного инсулина, которая будет использована или потрачена на снижение уровня глюкозы в следующем временном интервале, путем вычитания из 100% процентной доли, которая была определена на шаге 1414 в отношении текущего временного интервала.
Специалистам должно быть понятно, что в качестве альтернативы шаги 1414 и 1416 можно модифицировать так, чтобы схема 1214 управления приводилась в действие для вычисления (процентной) доли PR действия болюсного инсулина, которая будет использована или потрачена на снижение уровня глюкозы в следующем временном интервале, как функции TC, TBNI и TM с использованием любого из методов, описанных выше в отношении шага 1414, и последующего определения (процентной) доли PU действия болюсного инсулина, которая будет использована или потрачена на снижение уровня глюкозы в текущем временном интервале, путем вычитания из 100% процентной доли, которая была определена на шаге 1414 в отношении следующего временного интервала.
В любом случае выполнение алгоритма переходит с шага 1416 на шаг 1418, на котором схема 1214 управления приводится в действие для вычисления количества корректирующего болюса СВ согласно уравнению (10), т.е. по формуле: СВ=(BGM-BGTNI)/[(ISCI*PU)+(ISNI*PR)]. С любого из шагов 1412 или 1418 выполнение алгоритма возвращается на начало шага 1402.
В одном варианте калькулятора болюсов он может обеспечивать уменьшение дозы инсулина на основании скорости изменения глюкозы. Если уровень глюкозы уменьшается, в калькуляторе болюсов может использоваться прогнозное значение глюкозы. Прогнозное значение глюкозы может вычисляться с использованием линейного прогнозирования или при помощи другой модели. Такая модель может быть физиологической моделью. Обычно вычисляемая корректирующая доза болюса основывается на расхождении между текущим и целевым значениями глюкозы. В одном варианте осуществления изобретения корректирующая доза (Ic) может определяться при помощи уравнения:
Figure 00000067
.
В подобном уравнении с - поправочный коэффициент, g - уровень глюкозы, a gt - прогнозный уровень глюкозы в момент времени t.
В случае если скорость изменения глюкозы известна и является отрицательной, процесс вычисления можно изменить, чтобы использовать прогнозное значение. В нижеследующем случае для определения прогнозного значения глюкозы используется линейная модель
Figure 00000068
.
В этом случае τ - время прогноза, которое может быть задаваться в зависимости от степени осведомленности больного диабетом о гипогликемии и степени неприятия им гипогликемии. По меньшей мере в одном варианте осуществления изобретения для времени прогноза может использоваться диапазон от примерно 20 минут до примерно 60 минут.
В случае использования метода, учитывающего риски, для определения корректирующей дозы используется следующее уравнение:
Figure 00000069
,
где gadj - поправка на максимально допустимую глюкозу.
Если рекомендованная отрицательная корректирующая доза превышает прандиальный болюс (болюс на прием пищи), в результате чего сумма двух болюсов будет отрицательной, для учета разности система может рекомендовать больному диабетом увеличить количество углеводов в еде. В качестве альтернативы, для учета разности система также может уменьшить базальную скорость.
Примеры
На фиг. 15-18 показаны примеры осуществления настоящего изобретения, предусматривающие оценку данных непрерывного мониторирования гликемии (НМГ) по меньшей мере в одном варианте осуществления предлагаемого в изобретении способа. На фиг. 15 представлена оценка набора данных постпрандиального НМГ, в результате которой выдается предупреждение. По меньшей мере на некоторых из фиг. 15-18 используются следующие обозначения: прием пищи - 1500, кривая данных НМГ - 1501, результат измерения при самоконтроле гликемии (СГК) - 1502, результат СГК, указывающий на необходимость корректирующего болюса - 1502b, постпрандиальный трапецоид максимально допустимой глюкозы - 1503, максимально допустимая глюкоза, повышенная из-за корректирующей дозы - 1503b, первый результат НМГ - 1504, неопределенность, связанная с первым результатом НМГ - 1505, второй результат НМГ - 1506, неопределенность, связанная со вторым результатом НМГ - 1507, верхнее целевое значение - 1508, нижнее целевое значение - 1509, целевое значение - 1510, результат НМГ - 1511, неопределенность, связанная с соответствующим результатом НМГ - 1512, прогноз гипогликемии - 1513, прогнозное значение глюкозы, полученное на основании модели - 1514, кривая данных НМГ с нулевой скоростью изменения глюкозы во время приема пищи - 1515, кривая данных НМГ с отрицательной скоростью изменения глюкозы во время приема пищи - 1516, модифицированный корректирующий болюс, основанный на прогнозном значении глюкозы - 1517 и время прогноза - 1518.
На фиг. 15 после постпрандиального трапецоида, связанного с пищевым подъемом глюкозы, сигнал данных НМГ остается высоким. Если значение для I0 было задано равным 0, тревожное сообщение будет выдано при первом результате НМГ (1504). Если же минимальная доза инсулина была задана равной N, система не будет выдавать тревожного сообщения, пока не поступит второй результат НМГ (1506).
На фиг. 16 постпрандиальные данные НМГ не приводят к выдаче предупреждения, поскольку вероятность необходимости введения корректирующего болюса в случае, если было выполнено измерение при СГК, так и не превысило необходимого порога. На фиг. 15 и 16 система не выдает предупреждения, поскольку вероятность возникновения потребности в корректирующем болюсе, составляющем по меньшей мере I0, не превышает требуемого порога.
На фиг. 17 после корректирующего болюса следует предупреждение. В частности, показанный на фиг. 17 пример идентичен показанному на фиг. 15, за исключением того, что трапецоид является результатом корректирующей дозы, а не приема пищи.
На фиг. 18 в пищевой болюс включено снижение безопасности, когда скорость изменения глюкозы отрицательна.
На фиг. 19 показано три сценария выдачи тревожного сообщения об опасности наложения инсулина. Ниже описываются обозначения, используемые на фиг. 19:
- 1900 - болюс, вычисленный калькулятором болюсов на основании результата СГК (1902b)
- 1901 - кривая данных НМГ
- 1902b - результат измерения при СГК, используемый калькулятором болюсов
- 1903b - максимально допустимая глюкоза крови, определенная калькулятором болюсов
- 1904а - запрос заданного вручную болюса для сценария А
- 1906а - значение глюкозы по результатам НМГ для сценария А
- 1907а - значение gρ для сценария А
- 1908а - для сценария А выдается тревожное сообщение, поскольку значение IS будет отрицательным, так как gρ<gm. Таким образом, без результата измерения при СГК или соответствующих данных об углеводах, этот запрос болюса вероятно приведет к наложению инсулина.
- 1904b - запрос заданного вручную болюса для сценария В
- 1905b - максимальный безопасный болюс IS и значение болюса, выбранное для сценария В
- 1906b - значение глюкозы по результатам НМГ для сценария В
- 1907b - значение gρ для сценария В
- 1908b - для сценария В выдается тревожное сообщение, поскольку значение IS меньше значения I. В этом случае больной диабетом выбирает возможность уменьшения количества болюсного инсулина до максимального безопасного болюса (1905b)
- 1904с - запрос заданного вручную болюса для сценария С
- 1905с - запрошенный болюс вводится, так как он меньше максимального безопасного болюса IS для сценария С
- 1906с - значение глюкозы по результатам НМГ для сценария С
- 1907с - значение gρ для сценария С
- 1908с - тревожное сообщение для сценария С отсутствует, поскольку запрошенное количество болюсного инсулина меньше максимального безопасного болюса.
- 1909 - нижняя граница целевого диапазона
- 1910 - целевое значение
- 1911 - верхняя граница целевого диапазона
Использование по меньшей мере одного варианта реализации алгоритма для устранения наложения инсулина можно пояснить на нескольких отдельных примерах.
Пример 1: На помпе вручную задан болюс, причем имеется контроллер, который связан с помпой. Помпа посылает контроллеру запрос на введение болюса в организм. Контроллер определяет, нужно ли отправить тревожное сообщение. Контроллер посылает обратно в помпу состояние тревоги. Если тревожное сообщение необходимо, помпа предупреждает больного диабетом и предоставляет на выбор ряд возможных мер.
Пример 2: На помпе задан вручную болюс, причем контроллер отсутствует либо присутствует, но не связан с помпой. Помпа не может напрямую взаимодействовать со средствами НМГ. В подобном примере система не может реализовать этот алгоритм.
Возможный вариант 1: Принимать все заданные вручную болюсы.
Возможный вариант 2: Предложить больному диабетом запустить калькулятор болюсов и затем принять заданный вручную болюс.
Пример 3: На помпе вручную задан болюс, причем контроллер отсутствует либо присутствует, но не связан с помпой. Помпа не может напрямую взаимодействовать со средствами НМГ. Помпа запрашивает самый последний результат НМГ и его неопределенность. Помпа вычисляет максимальный безопасный болюс на основании истории болюсов и данных НМГ. Если необходимо тревожное сообщение, помпа предупреждает больного диабетом и предоставляет на выбор ряд возможных мер.
Следует иметь в виду, что описанные выше способы и системы могут использоваться для оценивания уровня глюкозы у больного диабетом и предупреждения больного диабетом при достижении заданной вероятности того, что потребуется определенная доза инсулина. Описанные выше способы могут быть запрограммированы и храниться на машиночитаемом носителе информации, на котором записаны машиноисполняемые (исполняемые на компьютере) команды для осуществления этих способов. К таким машиночитаемым носителям информации могут относиться компакт-диски, жесткие диски, флэш-накопители, запоминающие устройства с произвольной выборкой, динамические запоминающие устройства с произвольной выборкой, электрически перепрограммируемые постоянные запоминающие устройства и т.д.
Следует отметить, что указания на выполнение или настройку того или иного компонента, используемого при осуществлении настоящего изобретения, с возможностью воплощения определенного свойства или качества либо с возможностью выполнения определенной функции относятся к конструкции, а не назначению этого компонента. В частности, описание того, каким образом тот или иной компонент выполнен или настроен, означает существующее физическое состояние этого компонента и в этой связи должно трактоваться как описание конструктивных характеристик компонента.
Хотя изобретение описано выше в отношении его конкретных аспектов и на конкретных примерах его осуществления, осуществление изобретения в том виде, как оно заявлено, также возможно при внесении в него различных изменений. Кроме того, хотя выше были описаны различные аспекты изобретения, такие аспекты не обязательно должны использоваться совместно. Таким образом, любые вышеупомянутые изменения могут подпадать под патентные притязания, изложенные в прилагаемой формуле изобретения.

Claims (46)

1. Способ анализа уровня глюкозы у больного диабетом, включающий:
- получение вычислительным устройством измеренных значений глюкозы от системы непрерывного мониторирования гликемии, подключенной к больному диабетом, максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента;
- анализ измеренных значений глюкозы, выполняемый в вычислительном устройстве с использованием средства вероятностного анализа для определения граничного значения глюкозы (gρ) при пороге вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина;
- выполняемое вычислительным устройством определение порога глюкозы
Figure 00000070
с использованием максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента; и
- выполняемое вычислительным устройством сравнение граничного значения глюкозы (gρ) с порогом глюкозы
Figure 00000070
, причем если граничное значение глюкозы (gρ) превышает порог глюкозы
Figure 00000070
, то вычислительное устройство выдает через пользовательский интерфейс предупреждение.
2. Способ по п. 1, в котором граничное значение глюкозы (gρ) определяют по формуле:
Figure 00000071
где:
Figure 00000072
- вариация при измерении глюкозы,
ρ - порог вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина, и
g - измеренное значение глюкозы.
3. Способ по п. 1, в котором порог глюкозы
Figure 00000070
определяют по формуле:
Figure 00000073
где:
gm - максимально допустимая глюкоза,
с - поправочный коэффициент,
I0 - минимальная доза инсулина,
причем по меньшей мере одно из следующего: максимально допустимую глюкозу, минимальную дозу инсулина и поправочный коэффициент, необязательно получают из калькулятора болюсов.
4. Способ по п. 1, в котором в период сна заданная доза инсулина выше, чем в период бодрствования, или заданная доза инсулина равна одной единице.
5. Способ по п. 1, в котором порог глюкозы зависит от времени суток, и/или сутки разбивают на временные блоки, по меньшей мере один из которых имеет иной порог глюкозы, нежели остальные временные блоки.
6. Способ по п. 1, в котором порог вероятности выбирают из группы значений, состоящей из 99%, 98%, 90%, 75% и 50%, и/или измеренные значения глюкозы получают путем периодического измерения глюкозы, выполняемого каждую минуту.
7. Способ по п. 1, в котором предупреждение выдается посредством изображения, звука, вибрации или любой комбинации вышеперечисленного, и/или предупреждение указывает на необходимость измерения глюкозы.
8. Способ по одному из предыдущих пунктов, включающий также определение максимального безопасного болюса инсулина и, при необходимости, выдачу через пользовательский интерфейс тревожного сообщения, если введенный вручную запрос болюса превышает максимальный безопасный болюс инсулина.
9. Способ по одному из предыдущих пунктов, включающий также введение больному диабетом заданного страховочного болюса инсулина, если не проводится измерение глюкозы и не вводится корректирующая доза.
10. Способ по п. 9, в котором измерение глюкозы представляет собой измерение глюкозы крови в рамках самоконтроля гликемии или дополнительное измерение, выполняемое системой непрерывного мониторирования гликемии и инициируемое пользователем.
11. Способ по одному из предыдущих пунктов, включающий также определение вычислительным устройством болюсной дозы инсулина с использованием результата измерения глюкозы, причем при определении болюсной дозы инсулина используют:
- по меньшей мере одно из следующего: результат определения скорости изменения измеренных значений глюкозы, тенденцию в поведении измеренных значений глюкозы и закономерность в поведении измеренных значений глюкозы; или
- по меньшей мере один результат определения скорости изменения измеренных значений глюкозы и прогнозное значение глюкозы,
причем корректирующую дозу определяют по формуле:
Figure 00000074
где:
с - поправочный коэффициент,
g - измеренное значение глюкозы,
Figure 00000075
- скорость изменения глюкозы,
τ - прогнозное время, и
gt - прогнозный уровень глюкозы в момент времени t.
12. Способ по п. 11, в котором время t составляет от 20 до 60 минут.
13. Система для анализа уровня глюкозы у больного диабетом, включающая в себя:
- контроллер, содержащий машиночитаемый носитель информации с хранящимися в нем машиноисполняемыми командами, при выполнении которых в контроллере контроллером осуществляется:
- получение вычислительным устройством измеренных значений глюкозы от сенсора глюкозы, подключенного к больному диабетом, максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента,
- анализ измеренных значений глюкозы, выполняемый в вычислительном устройстве с использованием средства вероятностного анализа для определения граничного значения глюкозы (gρ) при пороге вероятности того, что больному диабетом потребуется по меньшей мере заданная доза инсулина,
- выполняемое вычислительным устройством определение порога глюкозы
Figure 00000070
с использованием максимально допустимой глюкозы, минимальной дозы инсулина и поправочного коэффициента, и
- выполняемое вычислительным устройством сравнение граничного значения глюкозы (gρ) с порогом глюкозы
Figure 00000070
, причем если граничное значение глюкозы (gρ) выше порога глюкозы
Figure 00000070
, то вычислительное устройство выдает предупреждение; и
- связанный с контроллером калькулятор болюсов, способный определять болюс инсулина на основании измеренных значений глюкозы.
14. Система по п. 13, в которой в машиночитаемом носителе информации также хранятся машиноисполняемые команды, при выполнении которых в контроллере контроллером осуществляется определение максимального безопасного болюса инсулина и/или определение скорости изменения глюкозы на основании измеренных значений глюкозы, причем скорость изменения глюкозы используется калькулятором болюсов для модификации болюса инсулина или используется для прогноза ожидаемого в будущем уровня глюкозы, предназначенного для использования калькулятором болюсов.
15. Система по п. 13, содержащая также пользовательский интерфейс, причем контроллер обеспечивает возможность отображения на пользовательском интерфейсе модифицированного болюса инсулина и/или возможность выдачи на пользовательский интерфейс тревожного сообщения, если запрос болюса, введенный вручную посредством пользовательского интерфейса, превышает максимальный безопасный болюс инсулина.
RU2017138101A 2015-04-02 2016-04-01 Способы и системы для анализа данных о содержании глюкозы, измеренных у больного диабетом RU2721878C2 (ru)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US14/677,148 US10617363B2 (en) 2015-04-02 2015-04-02 Methods and systems for analyzing glucose data measured from a person having diabetes
US14/677,148 2015-04-02
PCT/US2016/025502 WO2016161254A1 (en) 2015-04-02 2016-04-01 Methods and systems for analyzing glucose data measured from a person having diabetes

Publications (3)

Publication Number Publication Date
RU2017138101A RU2017138101A (ru) 2019-05-06
RU2017138101A3 RU2017138101A3 (ru) 2019-09-03
RU2721878C2 true RU2721878C2 (ru) 2020-05-25

Family

ID=55752772

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2017138101A RU2721878C2 (ru) 2015-04-02 2016-04-01 Способы и системы для анализа данных о содержании глюкозы, измеренных у больного диабетом

Country Status (10)

Country Link
US (1) US10617363B2 (ru)
EP (1) EP3277180B1 (ru)
CN (1) CN107809943B (ru)
CA (1) CA2979886C (ru)
ES (1) ES2944540T3 (ru)
FI (1) FI3277180T3 (ru)
HK (1) HK1252509A1 (ru)
HU (1) HUE061731T2 (ru)
RU (1) RU2721878C2 (ru)
WO (1) WO2016161254A1 (ru)

Families Citing this family (46)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7959598B2 (en) 2008-08-20 2011-06-14 Asante Solutions, Inc. Infusion pump systems and methods
US9561324B2 (en) 2013-07-19 2017-02-07 Bigfoot Biomedical, Inc. Infusion pump system and method
GB2523989B (en) 2014-01-30 2020-07-29 Insulet Netherlands B V Therapeutic product delivery system and method of pairing
CN107405446B (zh) 2015-02-18 2020-09-29 英赛罗公司 流体输送和输注装置及其使用方法
DE202015005019U1 (de) * 2015-07-18 2016-10-20 Rudolf King Methode zur Anbindung und Nutzung von externen Serverkapazitäten bei mobilen persönlichen und anderen Sicherheitssystemen und Trennung von Sensoren, Funktionsbereichen und Rechenkapazitäten mittels Alarmverlagerung auf externe Computer und Server
US10275573B2 (en) 2016-01-13 2019-04-30 Bigfoot Biomedical, Inc. User interface for diabetes management system
CA3009351A1 (en) 2016-01-14 2017-07-20 Bigfoot Biomedical, Inc. Adjusting insulin delivery rates
US10610643B2 (en) 2016-01-14 2020-04-07 Bigfoot Biomedical, Inc. Occlusion resolution in medication delivery devices, systems, and methods
US10765807B2 (en) 2016-09-23 2020-09-08 Insulet Corporation Fluid delivery device with sensor
WO2018111928A1 (en) 2016-12-12 2018-06-21 Mazlish Bryan Alarms and alerts for medication delivery devices and related systems and methods
US11033682B2 (en) 2017-01-13 2021-06-15 Bigfoot Biomedical, Inc. Insulin delivery methods, systems and devices
US10881792B2 (en) 2017-01-13 2021-01-05 Bigfoot Biomedical, Inc. System and method for adjusting insulin delivery
WO2018204568A1 (en) 2017-05-05 2018-11-08 Eli Lilly And Company Closed loop control of physiological glucose
EP3758008A1 (en) * 2017-06-15 2020-12-30 Novo Nordisk A/S Basal titration with adaptive target glucose level
DE112017007606T5 (de) 2017-06-30 2020-02-27 Mitsubishi Electric Corporation Instabilitätsdetektionvorrichtung, instabilitätsdetektionssystem und instabilitätsdetektionsverfahren
CN109692010A (zh) * 2017-10-20 2019-04-30 爱科来株式会社 测定装置、计算机可读取的记录介质和测定方法
WO2019123510A1 (ja) * 2017-12-18 2019-06-27 三菱電機株式会社 表示制御装置、表示システム、表示装置、表示方法および表示プログラム
WO2019125932A1 (en) 2017-12-21 2019-06-27 Eli Lilly And Company Closed loop control of physiological glucose
USD928199S1 (en) 2018-04-02 2021-08-17 Bigfoot Biomedical, Inc. Medication delivery device with icons
AU2019263490A1 (en) 2018-05-04 2020-11-26 Insulet Corporation Safety constraints for a control algorithm-based drug delivery system
WO2019246381A1 (en) 2018-06-22 2019-12-26 Eli Lilly And Company Insulin and pramlintide delivery systems, methods, and devices
CN112601934B (zh) 2018-09-03 2023-02-17 三菱电机株式会社 信号显示控制装置以及计算机可读取的记录介质
EP3856285A1 (en) 2018-09-28 2021-08-04 Insulet Corporation Activity mode for artificial pancreas system
CN113226448A (zh) * 2018-10-06 2021-08-06 基耐普提克有限公司 用于可穿戴电子数字治疗设备的方法和装置
WO2020077223A1 (en) 2018-10-11 2020-04-16 Insulet Corporation Event detection for drug delivery system
USD920343S1 (en) 2019-01-09 2021-05-25 Bigfoot Biomedical, Inc. Display screen or portion thereof with graphical user interface associated with insulin delivery
US12114972B2 (en) * 2019-02-01 2024-10-15 Medtronic Minimed, Inc. Methods, systems, and devices for continuous glucose monitoring
US11801344B2 (en) * 2019-09-13 2023-10-31 Insulet Corporation Blood glucose rate of change modulation of meal and correction insulin bolus quantity
US11935637B2 (en) 2019-09-27 2024-03-19 Insulet Corporation Onboarding and total daily insulin adaptivity
WO2021113647A1 (en) 2019-12-06 2021-06-10 Insulet Corporation Techniques and devices providing adaptivity and personalization in diabetes treatment
US11833329B2 (en) 2019-12-20 2023-12-05 Insulet Corporation Techniques for improved automatic drug delivery performance using delivery tendencies from past delivery history and use patterns
WO2021138171A1 (en) 2019-12-30 2021-07-08 Roche Diabetes Care, Inc. Computer-implemented diabetes management methods
WO2021141941A1 (en) 2020-01-06 2021-07-15 Insulet Corporation Prediction of meal and/or exercise events based on persistent residuals
US11551802B2 (en) 2020-02-11 2023-01-10 Insulet Corporation Early meal detection and calorie intake detection
US11986630B2 (en) 2020-02-12 2024-05-21 Insulet Corporation Dual hormone delivery system for reducing impending hypoglycemia and/or hyperglycemia risk
US11547800B2 (en) 2020-02-12 2023-01-10 Insulet Corporation User parameter dependent cost function for personalized reduction of hypoglycemia and/or hyperglycemia in a closed loop artificial pancreas system
US11324889B2 (en) 2020-02-14 2022-05-10 Insulet Corporation Compensation for missing readings from a glucose monitor in an automated insulin delivery system
US11607493B2 (en) 2020-04-06 2023-03-21 Insulet Corporation Initial total daily insulin setting for user onboarding
US12031982B2 (en) 2020-04-19 2024-07-09 John J. Daniels Using exhaled breath condensate for testing for a biomarker of COVID-19
CN112133439B (zh) * 2020-09-07 2022-06-03 北京理工大学 一种基于高斯过程的餐前胰岛素剂量个体化决策系统
EP4221588A1 (en) 2020-09-30 2023-08-09 Insulet Corporation Secure wireless communications between a glucose monitor and other devices
US11904140B2 (en) 2021-03-10 2024-02-20 Insulet Corporation Adaptable asymmetric medicament cost component in a control system for medicament delivery
US11738144B2 (en) 2021-09-27 2023-08-29 Insulet Corporation Techniques enabling adaptation of parameters in aid systems by user input
US11439754B1 (en) 2021-12-01 2022-09-13 Insulet Corporation Optimizing embedded formulations for drug delivery
CN114464291B (zh) * 2021-12-22 2022-11-25 北京理工大学 一种基于贝叶斯优化的mdi剂量建议系统
WO2024147928A1 (en) 2023-01-06 2024-07-11 Insulet Corporation Automatically or manually initiated meal bolus delivery with subsequent automatic safety constraint relaxation

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090105572A1 (en) * 2007-10-18 2009-04-23 Lifescan Scotland, Ltd. Method for predicting a user's future glycemic state
US20120166126A1 (en) * 2010-12-22 2012-06-28 Roche Diagnostics Operations, Inc. Calibration of a handheld diabetes managing device that receives data from a continuous glucose monitor
US20140083867A1 (en) * 2012-09-21 2014-03-27 Lifescan, Inc. Method and system to derive multiple glycemic patterns from glucose measurements during time of the day
EP2762073A1 (en) * 2010-01-26 2014-08-06 Roche Diagnostics GmbH Methods and systems for processing glucose data measured from a person having diabetes
US20140235981A1 (en) * 2008-08-31 2014-08-21 Abbott Diabetes Care Inc. Closed Loop Control With Improved Alarm Functions

Family Cites Families (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1985002546A1 (en) * 1983-12-08 1985-06-20 Disetronic Ag Infusion device
US5420108A (en) * 1992-09-14 1995-05-30 Shohet; Isaac H. Method of controlling diabetes mellitus
US7278983B2 (en) * 2002-07-24 2007-10-09 Medtronic Minimed, Inc. Physiological monitoring device for controlling a medication infusion device
US7291107B2 (en) 2004-08-26 2007-11-06 Roche Diagnostics Operations, Inc. Insulin bolus recommendation system
US8758245B2 (en) * 2007-03-20 2014-06-24 Lifescan, Inc. Systems and methods for pattern recognition in diabetes management
WO2009075925A1 (en) 2007-12-13 2009-06-18 Shaya Steven A Method and apparatus to calculate diabetic sensitivity factors affecting blood glucose
US9220456B2 (en) * 2008-04-04 2015-12-29 Hygieia, Inc. Systems, methods and devices for achieving glycemic balance
GB2466183A (en) 2008-12-09 2010-06-16 Cambridge Entpr Ltd Closed loop diabetes management system
JP2012519018A (ja) * 2009-02-26 2012-08-23 モール・リサーチ・アプリケーションズ・リミテッド 糖尿病に関連した治療を自動モニタリングする方法およびシステム
US8257300B2 (en) 2009-05-22 2012-09-04 Abbott Diabetes Care Inc. Safety features for integrated insulin delivery system
DK3988470T3 (da) * 2009-08-31 2023-08-28 Abbott Diabetes Care Inc Visningsindretninger til en medicinsk indretning
US20110071464A1 (en) 2009-09-23 2011-03-24 Medtronic Minimed, Inc. Semi-closed loop insulin delivery
EP2644088B1 (en) * 2009-09-29 2017-01-18 Lifescan Scotland Limited Analyte testing method and device for diabetes management
US8579879B2 (en) 2010-02-19 2013-11-12 Medtronic Minimed, Inc. Closed-loop glucose control startup
US8532933B2 (en) 2010-06-18 2013-09-10 Roche Diagnostics Operations, Inc. Insulin optimization systems and testing methods with adjusted exit criterion accounting for system noise associated with biomarkers
US20120006100A1 (en) * 2010-07-06 2012-01-12 Medtronic Minimed, Inc. Method and/or system for determining blood glucose reference sample times
US8774889B2 (en) 2010-12-22 2014-07-08 Roche Diagnostics Operations, Inc. Patient monitoring system with efficient pattern matching algorithm
US20140088392A1 (en) 2011-02-11 2014-03-27 Abbott Diabetes Care Inc. Feedback from Cloud or HCP to Payer or Patient via Meter or Cell Phone
CN103907116A (zh) 2011-08-26 2014-07-02 弗吉尼亚大学专利基金会 用于糖尿病的自适应咨询控制的方法、系统以及计算机可读介质
US10881339B2 (en) 2012-06-29 2021-01-05 Dexcom, Inc. Use of sensor redundancy to detect sensor failures
US20140066884A1 (en) 2012-08-30 2014-03-06 Medtronic Minimed, Inc. Sensor model supervisor for a closed-loop insulin infusion system
US9623179B2 (en) * 2012-08-30 2017-04-18 Medtronic Minimed, Inc. Safeguarding techniques for a closed-loop insulin infusion system
US10463282B2 (en) 2012-10-04 2019-11-05 Roche Diabetes Care, Inc. System and method for assessing risk associated with a glucose state
US9119528B2 (en) 2012-10-30 2015-09-01 Dexcom, Inc. Systems and methods for providing sensitive and specific alarms
WO2015073211A1 (en) 2013-11-14 2015-05-21 Regents Of The University Of California Glucose rate increase detector: a meal detection module for the health monitoring system
US10722650B2 (en) 2014-03-28 2020-07-28 Roche Diabetes Care, Inc. System and method for adjusting therapy based on risk associated with a glucose state
US9465917B2 (en) 2014-05-30 2016-10-11 Roche Diabetes Care, Inc. Hazard based assessment patterns
AU2016269688A1 (en) * 2015-05-29 2018-01-18 Insulcloud, S.L. Monitoring device for drug application with a drug pen, with logging, communication and alarms

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090105572A1 (en) * 2007-10-18 2009-04-23 Lifescan Scotland, Ltd. Method for predicting a user's future glycemic state
US20140235981A1 (en) * 2008-08-31 2014-08-21 Abbott Diabetes Care Inc. Closed Loop Control With Improved Alarm Functions
EP2762073A1 (en) * 2010-01-26 2014-08-06 Roche Diagnostics GmbH Methods and systems for processing glucose data measured from a person having diabetes
US20120166126A1 (en) * 2010-12-22 2012-06-28 Roche Diagnostics Operations, Inc. Calibration of a handheld diabetes managing device that receives data from a continuous glucose monitor
US20140083867A1 (en) * 2012-09-21 2014-03-27 Lifescan, Inc. Method and system to derive multiple glycemic patterns from glucose measurements during time of the day

Also Published As

Publication number Publication date
EP3277180A1 (en) 2018-02-07
WO2016161254A1 (en) 2016-10-06
US20160287184A1 (en) 2016-10-06
CA2979886C (en) 2024-05-21
CN107809943B (zh) 2022-02-11
RU2017138101A3 (ru) 2019-09-03
FI3277180T3 (fi) 2023-05-23
US10617363B2 (en) 2020-04-14
RU2017138101A (ru) 2019-05-06
EP3277180B1 (en) 2023-03-22
CN107809943A (zh) 2018-03-16
HUE061731T2 (hu) 2023-08-28
CA2979886A1 (en) 2016-10-06
ES2944540T3 (es) 2023-06-22
HK1252509A1 (zh) 2019-05-31

Similar Documents

Publication Publication Date Title
RU2721878C2 (ru) Способы и системы для анализа данных о содержании глюкозы, измеренных у больного диабетом
US10722650B2 (en) System and method for adjusting therapy based on risk associated with a glucose state
EP3479262B1 (en) Control-to-range failsafes
US11605461B2 (en) Risk-based control-to-range blood glucose diabetes management device
CN109564774B (zh) 控制到范围积极性
RU2737339C2 (ru) Способ и устройство для оценивания уровня глюкозы у больного диабетом и автоматического выбора режима разомкнутого или замкнутого управления устройством для введения лекарственного средства
US20180116589A1 (en) Computer-based diabetes management
EP2260462B1 (en) System for optimizing a patient&#39;s insulin dosage regimen
US20090164190A1 (en) Physiological condition simulation device and method
US20230233112A1 (en) Optimizing medication dosage based on analyte sensor data
KR20230113361A (ko) 자동 약물 전달 시스템을 위한 간단한 식사 알림을위한 디바이스 및 방법
EP4283624A1 (en) Assessment of past insulin delivery outcomes to automatically tune mda systems