CN107280664A - 一种基于疼痛脑响应的术前及术中麻醉深度评估的方法 - Google Patents

一种基于疼痛脑响应的术前及术中麻醉深度评估的方法 Download PDF

Info

Publication number
CN107280664A
CN107280664A CN201710614670.6A CN201710614670A CN107280664A CN 107280664 A CN107280664 A CN 107280664A CN 201710614670 A CN201710614670 A CN 201710614670A CN 107280664 A CN107280664 A CN 107280664A
Authority
CN
China
Prior art keywords
pain
patient
anesthesia
brain
laser
Prior art date
Legal status (The legal status 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 status listed.)
Pending
Application number
CN201710614670.6A
Other languages
English (en)
Inventor
胡理
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Psychology of CAS
Original Assignee
Institute of Psychology of CAS
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 Institute of Psychology of CAS filed Critical Institute of Psychology of CAS
Priority to CN201710614670.6A priority Critical patent/CN107280664A/zh
Publication of CN107280664A publication Critical patent/CN107280664A/zh
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/369—Electroencephalography [EEG]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • 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/4821—Determining level or depth of anaesthesia
    • 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/4824—Touch or pain perception evaluation
    • A61B5/4827—Touch or pain perception evaluation assessing touch sensitivity, e.g. for evaluation of pain threshold
    • 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • 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
    • 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)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Pain & Pain Management (AREA)
  • Anesthesiology (AREA)
  • Psychology (AREA)
  • Hospice & Palliative Care (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

本发明公开了一种基于疼痛脑响应的术前及术中麻醉深度评估的方法，麻醉前对病人通过激光诱发非手术刺激的疼痛；记录病人每次疼痛刺激所产生的脑电响应信号、口头评价疼痛强度及对应的最强疼痛刺激能量E0；对脑电信号进行信噪比增强处理，并提取N2、P2脑响应的潜伏期变异和振幅变异的参数；建立疼痛强度与N2、P2脑响应参数之间的线性模型；对病人实行麻醉处理，并将麻醉前所得最强激光刺激能量为E0的非手术诱发的激光刺激作用到麻醉病人身体上，对激光刺激所产生的脑电信号执行信噪比增强处理，并提取N2、P2脑响应的潜伏期变异和振幅变异的参数，通过计算得到预测疼痛强度，并通过此疼痛强度，对病人术前及术中的麻醉深度进行评估。本发明将测算时间控制在2秒以内，能够很好地应用于局麻以及全麻的术前及术中检测。

Description

一种基于疼痛脑响应的术前及术中麻醉深度评估的方法
技术领域
本发明属于术前麻醉技术领域，具体涉及一种基于疼痛脑响应的术前及术中麻醉深度评估的方法，用于指导实时调整镇痛方案。
背景技术
术前及术中麻醉深度的精确评估和判断已成为临床麻醉医师高度关注和亟待解决的一项难题。全身麻醉是一种特殊而复杂的状态，包括镇静催眠、记忆缺失、镇痛、应激抑制和肌肉松弛等多方面因素，并历经全麻诱导、维持及苏醒的全过程，在整个过程中均通过使用麻醉药物来控制。
不同手术方式、不同人群及不同程度的术中刺激对麻醉药物反应性亦不同。麻醉药物用量不足造成的麻醉深度过浅，容易导致血流动力学波动、术中知晓和体动等不良反应，使手术无法顺利进行；麻醉药物过量造成的麻醉过深则对呼吸系统、循环系统在内的生命器官功能出现严重抑制，缺血缺氧，严重的可能导致患者脑功能损伤，危及生命。因此，术前及术中麻醉深度评估既有利于控制麻醉质量，又可利用最少的麻醉药物达到最佳的麻醉效果，缩短术后苏醒时间。
目前临床应用的术前麻醉，特别是局麻检测手段还仅限于行为测试或者病人的口头报告，并没有通过任何生理指标来反应麻醉实施的效果。而全麻的术前及术中检测应用的是术中监测技术，常见的有两种监测技术，分别为听觉诱发电位指数(A-line(OR)ARX-index,AAI)和脑电双频指数(bispectral index,BIS)。
BIS反应的是整个大脑皮层的抑制程度，是一种综合性计算结果；但是该技术并没有具体的解剖学定位，而且存在一定的延迟时间，在移动时间平均数模式下，通常的延迟时间为30-60秒。AAI反应的是从内侧膝状体和初级听觉皮层产生的中潜伏期听觉诱发电位波幅与潜伏期的变化，是大脑对声音刺激的主动反应，波形有明确的解剖学定位；然而该技术的局限性在于存在延迟时间(外因输入自动回归模式下)通常为2-6秒。
发明内容
因此，本发明为了解决现有AAI和BIS技术在全麻的术前及术中检测应用中所存在的监测延迟时间长以及不能直接对伤害性系统功能进行评估的技术问题，为此，本发明提供了一种基于疼痛脑响应的术前及术中麻醉深度评估的方法。
所采用的技术方案如下：
一种基于疼痛脑响应的术前及术中麻醉深度评估的方法，所述方法包括：
步骤1、麻醉前对病人通过激光诱发非手术刺激的疼痛；
步骤2、记录病人每次疼痛刺激所产生的脑电信号、口头评价疼痛强度及对应的最强疼痛刺激能量E0；
步骤3、对步骤2中的脑电信号进行信噪比增强处理，并提取N2、P2脑响应的潜伏期变异和振幅变异的参数；
步骤4、建立疼痛强度与N2、P2脑响应参数之间的线性模型；
步骤5、对病人实行麻醉处理，并将麻醉前所得到的最强刺激能量为E0的非手术诱发的激光刺激作用到麻醉病人身体上，对激光刺激所产生的脑电信号执行步骤3，获得该激光刺激相应的N2、P2的脑响应参数；
步骤6、将步骤5得到的N2、P2脑响应参数代入步骤4的线性模型中，得到预测的疼痛强度，从而对病人术前及术中的麻醉深度进行评估。
所述步骤1中采用激光束短暂辐射热脉冲刺激器刺激病人的手背或足背，诱发病人的疼痛。
所述步骤2中的记录病人麻醉前所产生疼痛强度及对应的最强疼痛刺激能量E0，其具体方法是：以相同步长逐渐增加激光刺激强度，直到病人的口头评价疼痛强度的评分出现6-8分为止，依次记录病人每次的脑电响应及口头评价疼痛强度I，并记录最强激光刺激的刺激能量E0。
每次激光刺激持续时间为3-6ms，起始刺激能量为1-2J，激光刺激强度增加的步长为0.2-0.5J，刺激间隔为10-15s。
激光刺激能量的递增序列重复3至5次，将每次所测得的最强激光能量E0取平均值为Et直接作用到麻醉病人身体上。
每次激光刺激后，激光束随机移动，以避免伤害感受器疲劳或致敏。
所述步骤2中对病人进行口头评价疼痛强度的测取方法是：将激光刺激能量直接作用到病人身上，病人根据承受的疼痛对疼痛强度作出口头评分，口头评价疼痛强度的评分标准为0-10分，其中0分为完全不痛，4分为刚开始有疼痛感觉，10分为最痛。
所述步骤3中的脑电信号进行信噪比增强处理，并提取N2、P2脑响应的潜伏期变异和振幅变异的参数，其具体方法是：
步骤3.1，对所采集的脑电信号进行带通滤波BPF以消除不相关的噪声；
步骤3.2，采用独立成分分析ICA去除眼动伪迹；
步骤3.3，采用公共空间模式滤波CSP去除无关伪迹；
步骤3.4，采用多元线性回归MLR提取N2、P2脑响应的潜伏期变异和振幅变异等相关参数:
f(t)＝β1yN(t)+β2y'N(t)+β3yp(t)+β4y'p(t)+β5 (1)
其中：f(t)是每次激光诱发的疼痛的脑响应与时间t之间的变化关系；
β1,β2,β3,β4,β5是每次诱发刺激响应中N2、P2成分的潜伏期变异和振幅变异的常量参数；
yN(t)和yP(t)是N2、P2波形的模板；
y'N(t)和y'P(t)分别是N2、P2模板关于时间的导数。
所述的带通滤波BPF采用1-30Hz进行；所述的独立成分分析ICA中，去除具有大EOG通道贡献和额叶头皮分布的独立成分；所述的CSP滤波时，选择对应于Xpost的最大特征值和Xpre的最小特征值的3个特征向量来重构所有通道每次的脑电信号响应。
所述步骤4中建立疼痛强度与N2、P2脑响应参数之间的线性模型的具体方法：
根据每次的口头评价疼痛强度Ik及在步骤3中所得到的对应试次的参数β1k,β2k,β3k,β4k,β5k建立线性回归模型：
通过最小二乘法获得线性回归模型系数α1,α2,α3,α4,α5和c值。
所述步骤6中获得预测疼痛强度的具体方法是：对麻醉后施加给病人身上的最强激光诱发的疼痛刺激所产生的脑电信号通过步骤3处理后，提取与N2、P2脑响应参数相对应的β1,β2,β3,β4,β5系数值；将β1,β2,β3,β4,β5以及步骤4中所得到的α1,α2,α3,α4,α5和c代入等式(2)中，计算出预测疼痛强度I0。
所述步骤6对病人术前及术中的麻醉深度进行评估的方法是，通过等式(2)中所得到的预测疼痛强度I0值越低，则麻醉程度越高；若I0≤4，则说明麻醉后的病人已暂时失去痛觉。
本发明技术方案，具有如下优点：
A.本发明通过激光刺激诱发疼痛，通过对麻醉前和麻醉后的病人进行激光刺激诱发疼痛，采集病人在两种情况下所产生的脑电信号，根据对麻醉前和麻醉后的病人的脑电信号进行信噪比增强处理及多元线性回归分析，提取出病人的疼痛信号特征，同时结合麻醉前病人的口头评价疼痛强度生成疼痛强度的线性模型，根据疼痛强度的线性模型可以一次性测得病人在麻醉状态下的疼痛强度，从而可以准确对麻醉后激光刺激诱发的疼痛进行实时评估，完成麻醉深度的检测，是对伤害性系统直接进行评估的麻醉深度检测方式。一般情况下本发明可以将延迟时间控制在2秒以内，甚至更快，具有时间延迟小和操作简便的特点，能够很好地应用于局麻以及全麻的术前及术中检测。
B.本发明通过带通滤波、独立成分分析去伪迹、共空间模式滤波相结合的方式有效地提高了激光诱发电位(Laser-evoked potentials，LEP)信号的信噪比，有利于对麻醉后病人的麻醉深度作出更准确的评估；同时，本发明同样适用于手背和脚背，因此操作便携性强、应用范围广，特别是能够应用于局麻的术前及术中检测，实现了对局麻以及全麻手术的麻醉监测。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明所提供的术前及术中麻醉深度评估方法流程示意图；
图2为图1所示麻醉深度评估数据处理过程示意图；
图3为图1中步骤3对所采集的脑电信号进行信噪处理流程示意图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1和图2所示，本发明提供了一种基于疼痛脑响应的术前及术中麻醉深度评估的方法，包括：
【S1】麻醉前对病人通过激光诱发非手术刺激的疼痛。
在病人未麻醉的状态下，采用激光束短暂辐射热脉冲刺激器刺激病人的手背或足背，诱发病人的疼痛。比如将激光脉冲指向病人左手手背。
【S2】记录病人每次激光刺激所产生的脑电信号、口头评价疼痛强度及对应的最强疼痛刺激能量E0。
采用相同步长逐渐增加激光刺激强度，直到病人的疼痛评分出现6-8分为止，依次记录病人每次的脑电响应及口头评价疼痛强度I，并记录病人评分最高时所对应的最强激光刺激的刺激能量E0。本发明采用每次激光刺激持续时间为3-6ms，起始刺激能量为1-2J，激光刺激强度增加的步长0.2-0.5J，刺激间隔为10-15s，直到病人的疼痛评分出现6-8分为止。每次激光刺激后，激光束随机移动，以避免伤害感受器疲劳或致敏。激光刺激能量的递增序列重复3至5次，将每次所测得的最强激光能量E0取平均值为Et，Et也是作为对麻醉后病人进行麻醉深度评价的关键参数。
比如，设置激光刺激持续时间为4ms，起始激光刺激能量为2.5J，逐渐增加激光刺激强度，每次激光刺激后，激光束在5×5cm的范围内随机移动约1厘米，激光刺激间期间隔(ISI)介于10和15s之间。每次激光刺激后让病人对疼痛强度进行口头打分(0-10分)，其中0分为完全不痛，4分为刚开始有疼痛感觉，10分为最痛。
本发明中优选的脑电设备使用64通道Brain Products系统(Brain ProductsGmbH，Munich，Germany；通带：0.01-100Hz；采样速率：1000Hz)，鼻尖用作参考通道，所有通道阻抗保持低于10kΩ。为了监测眼睛运动和眼睛眨眼，同时从4个表面电极记录电眼信号：一对放置在上眼睑和下眼睑上，另一对放置在左右轨道外角1厘米处。
【S3】对【S2】中的脑电信号进行信噪比增强处理，并提取N2、P2脑响应的潜伏期变异和振幅变异等相关参数，如图3所示。
步骤3.1，对所采集的脑电信号进行带通滤波(bandpass filtering,BPF)以消除不相关的噪声；将采集到的脑电信号经过1和30Hz的带通滤波，使用激光刺激前500ms到刺激后1000ms的分析窗口提取脑电波(Electroencephalogram,EEG)，并使用预刺激时间间隔进行基线校正。
步骤3.2，采用独立成分分析(independent component analysis,ICA)去除伪迹；
步骤3.3，采用共空间模式滤波(common spatial pattern,CSP)去除无关伪迹；选择对应于Xpost的最大特征值和Xpre的最小特征值的3个特征向量来重构所有通道每次的脑电信号响应。
步骤3.4、采用多元线性回归(multiple linear regression,MLR)将滤波后的脑电信号曲线生成疼痛强度与脑响应之间的线性模型，拟合出每次LEP响应中N2、P2成分的潜伏期变异和振幅变异的参数：
f(t)＝β1yN(t)+β2y'N(t)+β3yp(t)+β4y'p(t)+β5 (1)
其中：f(t)是每次激光诱发的疼痛的脑响应与时间t之间的变化关系；
β1,β2,β3,β4,β5是每次诱发刺激响应中N2、P2成分的潜伏期变异和振幅变异的常量参数；
这里的N2、P2为人类疼痛脑电成分中的一对复合波，反应的是Aδ纤维介导的对激光刺激的皮层响应，该对复合波在颅顶电极处波幅最大。
yN(t)和yP(t)是N2、P2波形的模板；
y'N(t)和y'P(t)分别是N2、P2模板关于时间的导数。
对于激光刺激后时间间隔(0至500ms)内的Cz电极记录的LEP信号，可以将MLR模型写成如下等式：
f(t)＝aNyN(t+lN)+apyp(t+lp) (3)
其中f(t)是每次激光刺激的激光诱发电位LEP波形随时间t变化的函数，yN(t)和yP(t)是N2、P2波形的模板(每个病人的模板可通过平均该病人所有LEP的波形得到)，aN和aP是N2、P2模板的权重，lN和lP分别是N2、P2模板的潜伏期偏移值。
使用泰勒展开，等式(3)可以写成：
f(t)≈aNyN(t)+lNaNy'N(t)+apyp(t)+lpapy'p(t)+ε (4)
然后将等式(4)转化为上面的等式(1)。
下表1中是对麻醉前的病人进行39试次的试验，通过【S3】用多元线性回归的系数拟合出每次激光刺激诱发电位响应中N2、P2成分的潜伏期变异和振幅变异的参数β1,β2,β3,β4,β5。
表1每个试次的变异参数以及疼痛评分
上表中所给定试次的激光刺激能量越大，激光刺激对象所给出的口头评分也越高。
【S4】建立疼痛强度与N2、P2脑响应参数之间的线性模型。
根据每次的口头评价疼痛强度Ik及在【S3】中所得到的对应试次的参数β1k,β2k,β3k,β4k,β5k建立线性模型：
通过最小二乘法获得线性回归模型系数α1,α2,α3,α4,α5和c值。
比如，通过对表1中的39个试次进行计算，得到α1＝0.9111,α2＝0.0107,α3＝0.0247,α4＝-0.0061,α5＝0.0765，以及c＝5.0957。
【S5】对病人实行麻醉处理，并将麻醉前所得到的最强刺激能量E0的非手术诱发的激光刺激作用到麻醉病人身体上，对激光刺激所产生的脑电信号执行【S3】，获得该激光刺激相应的N2、P2的脑响应参数。
对麻醉后施加给病人身上的最强激光诱发的疼痛刺激所产生的脑电信号通过【S3】处理后，提取与N2、P2脑响应参数相对应的β1,β2,β3,β4,β5系数值；将β1,β2,β3,β4,β5以及【S4】中所得到的α1,α2,α3,α4,α5和c代入等式(2)中，计算出预测疼痛强度I0。
比如，将上表中评分为8的所对应的激光刺激能量啊直接作用到麻醉后的病人身体上，执行【S3】后，得到β1＝-1.7194,β2＝6.8976,β3＝0.1719,β4＝5.9563,β5＝-4.7442，并带入得到I0的值为3.2078。
【S6】将【S5】得到的N2、P2脑响应参数代入【S4】的线性模型中，得到预测的疼痛强度，从而对病人术前及术中的麻醉深度进行评估。
具体地，通过等式(2)中所得到的预测疼痛强度I0值越低，则麻醉程度越高；若I0≤4，则说明麻醉后的病人已暂时失去痛觉。在【S5】中通过计算得到麻醉后病人的预测疼痛强度I0为3.2078，可以依据此值对麻醉后的病人作出已暂时失去痛觉的评估，病人适合进行手术处理。
显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (12)

1.一种基于疼痛脑响应的术前及术中麻醉深度评估的方法，其特征在于，所述方法包括：
步骤1、麻醉前对病人通过激光诱发非手术刺激的疼痛；
步骤2、记录病人每次疼痛刺激所产生的脑电信号、口头评价疼痛强度及对应的最强疼痛刺激能量E0；
步骤3、对步骤2中的脑电信号进行信噪比增强处理，并提取N2、P2脑响应的潜伏期变异和振幅变异的参数；
步骤4、建立疼痛强度与N2、P2脑响应参数之间的线性模型；
步骤5、对病人实行麻醉处理，并将麻醉前所得到的最强刺激能量E0的非手术诱发的激光刺激作用到麻醉病人身体上，对激光刺激所产生的脑电信号执行步骤3，获得该激光刺激相应的N2、P2的脑响应参数；
步骤6、将步骤5得到的N2、P2脑响应参数代入步骤4的线性模型中，得到预测的疼痛强度，从而对病人术前及术中的麻醉深度进行评估。
2.根据权利要求1所述的方法，其特征在于，所述步骤1中采用激光束短暂辐射热脉冲刺激器刺激病人的手背或足背，诱发病人的疼痛。
3.根据权利要求1所述的方法，其特征在于，所述步骤2中的记录病人麻醉前所产生疼痛强度及对应的最强疼痛刺激能量E0，其具体方法是：以相同步长逐渐增加激光刺激强度，直到病人的口头评价疼痛强度的评分出现6-8分为止，依次记录病人每次的脑电响应及口头评价疼痛强度I，并记录最强激光刺激的刺激能量E0。
4.根据权利要求3所述的方法，其特征在于，每次激光刺激持续时间为3-6ms，起始刺激能量为1-2J，激光刺激强度增加的步长为0.2-0.5J，刺激间隔为10-15s。
5.根据权利要求4所述的方法，其特征在于，激光刺激能量的递增序列重复3至5次，将每次所测得的最强激光能量E0取平均值为Et直接作用到麻醉病人身体上。
6.根据权利要求5所述的方法，其特征在于，每次激光刺激后，激光束随机移动，以避免伤害感受器疲劳或致敏。
7.根据权利要求1所述的方法，其特征在于，所述步骤2中对病人进行口头评价疼痛强度的测取方法是：将激光刺激能量直接作用到病人身上，病人根据承受的疼痛对疼痛强度作出口头评分，口头评价疼痛强度的评分标准为0-10分，其中0分为完全不痛，4分为刚开始有疼痛感觉，10分为最痛。
8.根据权利要求1所述的方法，其特征在于，所述步骤3中的脑电信号进行信噪比增强处理，并提取N2、P2脑响应的潜伏期变异和振幅变异的参数，其具体方法是：
步骤3.1，对所采集的脑电信号进行带通滤波BPF以消除不相关的噪声；
步骤3.2，采用独立成分分析ICA去除眼动伪迹；
步骤3.3，采用共空间模式滤波CSP去除无关伪迹；
步骤3.4，采用多元线性回归MLR提取N2、P2脑响应的潜伏期变异和振幅变异等相关参数:
f(t)＝β1yN(t)+β2y'N(t)+β3yp(t)+β4y'p(t)+β5 (1)
其中：f(t)是每次激光诱发的疼痛的脑响应与时间t之间的变化关系；
β1,β2,β3,β4,β5是每次诱发刺激响应中N2、P2成分的潜伏期变异和振幅变异的常量参数；
yN(t)和yP(t)是N2、P2波形的模板；
y'N(t)和y'P(t)分别是N2、P2模板关于时间的导数。
9.根据权利要求8所述的方法，其特征在于，所述的带通滤波BPF采用1-30Hz进行；所述的独立成分分析ICA中，去除具有大EOG通道贡献和额叶头皮分布的独立成分；所述的CSP滤波时，选择对应于Xpost的最大特征值和Xpre的最小特征值的3个特征向量来重构所有通道每次的脑电信号响应。
10.根据权利要求9所述的方法，其特征在于，所述步骤4中建立疼痛强度与N2、P2脑响应参数之间的线性模型的具体方法：
根据每次的口头评价疼痛强度Ik及在步骤3中所得到的对应试次的参数β1k,β2k,β3k,β4k,β5k建立线性回归模型：
<mrow> <msub> <mi>I</mi> <mi>k</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <msub> <mi>&amp;alpha;</mi> <mi>i</mi> </msub> <msub> <mi>&amp;beta;</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> <mo>+</mo> <mi>c</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
通过最小二乘法获得线性回归模型系数α1,α2,α3,α4,α5和c值。
11.根据权利要求10所述的方法，其特征在于，所述步骤6中获得预测疼痛强度的具体方法是：对麻醉后施加给病人身上的最强激光诱发的疼痛刺激所产生的脑电信号通过步骤3处理后，提取与N2、P2脑响应参数相对应的β1,β2,β3,β4,β5系数值；将β1,β2,β3,β4,β5以及步骤4中所得到的α1,α2,α3,α4,α5和c代入等式(2)中，计算出预测疼痛强度I0。
12.根据权利要求11所述的方法，其特征在于，所述步骤6对病人术前及术中的麻醉深度进行评估的方法是，通过等式(2)中所得到的预测疼痛强度I0值越低，则麻醉程度越高；若I0≤4，则说明麻醉后的病人已暂时失去痛觉。
CN201710614670.6A 2017-07-26 2017-07-26 一种基于疼痛脑响应的术前及术中麻醉深度评估的方法 Pending CN107280664A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710614670.6A CN107280664A (zh) 2017-07-26 2017-07-26 一种基于疼痛脑响应的术前及术中麻醉深度评估的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710614670.6A CN107280664A (zh) 2017-07-26 2017-07-26 一种基于疼痛脑响应的术前及术中麻醉深度评估的方法

Publications (1)

Publication Number Publication Date
CN107280664A true CN107280664A (zh) 2017-10-24

Family

ID=60103052

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710614670.6A Pending CN107280664A (zh) 2017-07-26 2017-07-26 一种基于疼痛脑响应的术前及术中麻醉深度评估的方法

Country Status (1)

Country Link
CN (1) CN107280664A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110267590A (zh) * 2017-12-29 2019-09-20 深圳迈瑞生物医疗电子股份有限公司 基于脑电的麻醉深度监测方法和装置
CN113143208A (zh) * 2021-03-12 2021-07-23 深圳大学 基于多维度测量的疼痛敏感性评估系统和方法
CN114209283A (zh) * 2021-12-27 2022-03-22 燕山大学 一种基于脑机空间复杂度的术中疼痛检测方法
CN114391810A (zh) * 2022-01-06 2022-04-26 深圳元海恒鑫医疗科技有限公司 一种基于脑电信号监测儿童麻醉靶控智能化方法及系统

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110267590A (zh) * 2017-12-29 2019-09-20 深圳迈瑞生物医疗电子股份有限公司 基于脑电的麻醉深度监测方法和装置
CN113143208A (zh) * 2021-03-12 2021-07-23 深圳大学 基于多维度测量的疼痛敏感性评估系统和方法
CN113143208B (zh) * 2021-03-12 2023-07-25 深圳大学 基于多维度测量的疼痛敏感性评估系统和方法
CN114209283A (zh) * 2021-12-27 2022-03-22 燕山大学 一种基于脑机空间复杂度的术中疼痛检测方法
CN114391810A (zh) * 2022-01-06 2022-04-26 深圳元海恒鑫医疗科技有限公司 一种基于脑电信号监测儿童麻醉靶控智能化方法及系统

Similar Documents

Publication Publication Date Title
Redinbaugh et al. Thalamus modulates consciousness via layer-specific control of cortex
CN104825159B (zh) 脑电同频检测方法、脑电同频刺激方法与装置
Walter Critical review: The technique and application of electro-encephalography
CN107280664A (zh) 一种基于疼痛脑响应的术前及术中麻醉深度评估的方法
DE112009002463T5 (de) Verfahren und Vorrichtung zur Messung und Behandlung von Schüttelfrost während einer therapeutischen Tempratursteuerung
CN107041747B (zh) 一种用于脊柱整复系统的呼吸指示装置
Mulheren et al. Vibration over the larynx increases swallowing and cortical activation for swallowing
CN106943141B (zh) 一种用于脊柱整复系统的肌电检测装置
CN107041746B (zh) 一种用于脊柱整复系统的磁疗装置
CN107041748B (zh) 一种用于脊柱整复系统心电触发磁疗装置
Jovanovic et al. KITE-BCI: A brain-computer interface system for functional electrical stimulation therapy
CN107041995B (zh) 一种用于脊柱整复系统的磁疗装置
CN107041750B (zh) 一种用于脊柱整复系统的经颅磁刺激装置
Glassman et al. Transcranial magnetic motor-evoked potentials in scoliosis surgery
CN106901730B (zh) 一种脊柱整复系统
CN107041749B (zh) 一种用于脊柱整复系统的皮温变色涂料
CN106901729B (zh) 一种用于脊柱整复系统的药物涂层
Johnson et al. Reduction of muscle artifact in electroencephalographic recording
CN107041745B (zh) 一种用于脊柱整复系统的红外检测装置
Halldórsdóttir Network neuroscience evaluation of the relation between electroencephalography signals during electrophysiological and movement test
Redinbaugh et al. Central thalamus modulates consciousness by controlling layer-specific cortical interactions
Zhang et al. Intraoperative neurological monitoring system for robot assisted minimally invasive spine surgery using electromyography
McAfee Assessing Neuronal Synchrony and Brain Function Through Local Field Potential and Spike Analysis
Walter et al. BCCI-a bidirectional cortical communication interface
Sălăgean et al. Local Field Potential Microstate Analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20171024

RJ01 Rejection of invention patent application after publication