脉搏波形态学研究及其在子痫前期的应用陈婉琳;江锋;陈新忠;封英;缪家骏;焦翠翠;陈杭【摘要】子痫前期是妊娠20周后出现的特有疾病,会导致血液系统、消化系统、神经系统等的异常改变,严重危害孕产妇和围生儿的健康.针对实际临床中存在的子痫前期检测方法复杂、不便于日常监护的问题,提出了基于光电容积脉搏波提取相关特征参数的方法.根据脉搏波下降支的形态特征及生理意义,提取出了新的特征参数——光电容积斜率指数(PSI),通过逐段分析下降支斜率来体现外周血管阻力、血管壁弹性等生理信息的细节变化.基于临床采集的50例孕妇(23名为子痫前期患者,27名为正常孕妇)的光电容积脉搏波数据,对PSI指数区分子痫前期患者与正常孕妇的有效性进行了验证.结果表明,真阳性20例,假阳性1例,真阴性26例,假阴性3例,灵敏性和特异性分别达到87.0％和96.3％,准确性达到92.0％.PSI指数为子痫前期的无创诊断研究提供了理论依据和新的研究方法,对于实现子痫前期的早期筛查、保障孕妇和胎儿的健康有着重要意义.【期刊名称】《天津大学学报》【年(卷),期】2019(052)008【总页数】5页(P857-861)【关键词】光电容积脉搏波;子痫前期;光电容积斜率指数;形态学【作者】陈婉琳;江锋;陈新忠;封英;缪家骏;焦翠翠;陈杭【作者单位】浙江大学生物医学工程与仪器科学学院,杭州 310027;生物医学工程教育部重点实验室,杭州 310027;浙江省心脑血管检测技术与药效评价重点实验室,杭州 310027;浙江大学生物医学工程与仪器科学学院,杭州 310027;生物医学工程教育部重点实验室,杭州 310027;浙江省心脑血管检测技术与药效评价重点实验室,杭州 310027;浙江省心脑血管检测技术与药效评价重点实验室,杭州 310027;浙江大学医学院附属妇产科医院,杭州 310027;浙江大学医学院附属妇产科医院,杭州310027;浙江大学生物医学工程与仪器科学学院,杭州 310027;生物医学工程教育部重点实验室,杭州 310027;浙江省心脑血管检测技术与药效评价重点实验室,杭州310027;浙江大学医学院附属妇产科医院,杭州 310027;浙江大学生物医学工程与仪器科学学院,杭州 310027;生物医学工程教育部重点实验室,杭州 310027;浙江省心脑血管检测技术与药效评价重点实验室,杭州 310027【正文语种】中文【中图分类】TK448.21子痫前期是一种在妊娠20周后出现的疾病，通常表现为高血压、尿蛋白，并可能伴有水肿、头痛和腹部疼痛等[1]，是导致孕产妇和围生儿死亡的主要原因[2]．在临床中，子痫前期的诊断主要通过血压和尿蛋白的检测，该方法检测操作繁琐、无法实现实时检测，不利于孕妇的日常监护．子痫前期会导致全身小动脉痉挛、外周血管阻力增加、血管顺应性下降、心输出量减少等，这些病理生理变化会改变脉搏波的波形高度、波形斜率等形态特征．因此，从形态学角度，光电容积脉搏波与子痫前期之间存在着相关性．目前，已有一些光电容积脉搏波形态学与子痫前期相关的研究．Euliano等[3]研究发现，脉搏波波峰-重搏波波峰时长等参数在子痫前期患者和正常孕妇中差异明显；文献[4-5]研究发现，通过脉搏波传导时间变异性、波峰幅值等参数可预测高血压．但由于研究中存在数据量不足、适用范围有限等问题，这些参数仍未用在临床中．本课题组在光电容积脉搏波形态学研究方面已有一定的基础，提取出了一些新的特征参数，如脉搏波传导时间变异性、下降支衰减常数(QRC)、下降支面积比(ARD)、上升支面积比(ARA)、全周期面积比(ART)、搏动间期(PBI)等已应用于低血压预测[6]、心血管系统状态评估[7-8]、伤害感受监测[9-11]、麻醉深度监测[12-13]等临床研究中，取得了一定的成果．本文就光电容积脉搏波的形态学进行深入研究和探讨，提取出新的特征参数——光电容积斜率指数，该参数可有效区分子痫前期患者和正常孕妇．光电容积脉搏波(photoplethysmography，PPG)波形由上升支和下降支两部分组成，如图1所示．上升支表示心室快速射血期血液流过指端微血管床的流量变化情况．相比于上升支，下降支则相对更加复杂，PB段表示心室舒张，血液容积下降．BCU段称为重搏波，表示“血液返流”．脉搏波的形态特征反映了外周血管阻力大小、心输出量大小等状态．点的斜率表示血容量变化的快慢，外周阻力增大时，下降支前段的斜率较小，下降支后段的斜率较大；点P、U连接线的斜率表示下降支整体血流量变化快慢的平均水平．本研究将以下降支作为主要研究对象，进行特征参数的提取．实验数据由浙江大学医学院附属妇产科医院提供，医护人员通过GE监护仪B650对50名孕妇进行了光电容积脉搏波的采集，其中23名患有子痫前期的孕妇设为实验组，27名正常孕妇设为对照组．处理的原始信号是从监护仪中导出的csv文件数据，采样频率为100Hz．如图2为GE监护仪B650导出的其中一例原始脉搏波数据，横坐标为采样时间，纵坐标为脉搏波相对高度．由于存在高频扰动，实验通过smooth函数对原始脉搏波信号y进行平滑处理，返回与原信号等长的信号，记为y1，其中，xi表示每个点的采样时间，i表示数据点序号，，n为数据点长度．数学原理为波峰P和波谷U是提取的两个主要特征点．波峰P是一个脉搏波周期中的最高点．对于最大值的检索分为两个部分，首先根据一阶导数找到极大值，其次判断该极大值在某个范围内是最大值．实验中通过“搜索窗”函数实现对某个范围内最大值的检索．“搜索窗”是以极大值点为中点的具有一定宽度的“窗”，窗函数返回的是窗内最大值的横坐标．由于单个脉搏波周期为0.5～1.0s，为保证每个极大值能与其附近至少半个周期范围内的数据点进行比较，实验将搜索窗的宽度定义为0.6s，即检索极大值前后各0.3s范围内的最大值，并返回该最大的横坐标．当最大值的横坐标等于极大值的横坐标时，则说明该极大值为该周期脉搏波的波峰，否则该极大值为重搏波点．完成所有波峰的检测后，波谷U即为相邻两个波峰之间的最小值．如图3所示，以一个心动周期为例，将波谷U和波谷U′连接，记连接两个波谷的线段的斜率为k′，记该周期中各点修正后的值为y2，y2和k′满足公式将去除基线漂移的脉搏波起始点归零，得到处理后的脉搏波波形，记为Y，用于之后新的特征参数提取．光电容积斜率指数(PSI)旨在通过逐段分析脉搏波下降支的斜率，体现脉搏波形态变化更多的细节．将光电容积脉搏波下降支分为横坐标相等的10段．由于脉搏波采样率有限，若等分处的横坐标不存在数据点，则向下取最近的点作为代替．这样的取点方法会对特征参数计算产生一定的误差，实验前期已进行评估分析，发现该误差对于最终子痫前期与正常孕妇的分类结果并无影响，如图4所示．光电容积斜率指数PSIij是指脉搏波下降支曲线的斜率，如PSI26即为点P2、P6连线的斜率．计算公式为式中：k为点P0、P10连线的斜率；x0和Y0、x10和Y10分别为点P1、P10的横纵坐标；xi、xj为点Pi、Pj的横坐标，s；Yi、Yj为点Pi、Pj的纵坐标．PSIij 通过除以k值对脉搏波斜率进行了归一化处理．图5是斜率曲线，绿色曲线表示根据传统斜率公式计算得到的斜率，红色曲线表示归一化后的斜率PSI．从图中可以明显看出归一化后的斜率值PSI更加稳定，降低了个体差异带来的误差．通过光电容积斜率指数，进行区分子痫前期患者与正常孕妇的研究．基于先前的分析，本文选择了PSI1(7)进行了相关研究，如图6是实验组和对照组的PSI1(7)分布曲线，实验组的PSI1(7)分布在1.2～1.4之间，对照组的PSI1(7)分布在1.0～1.3之间，可明显观察到实验组的平均值大于对照组的平均值．表1为PSI1(7)的分析参数的统计结果．PSI1(7)的Mann-Whitney U检验P值小于0.05，说明实验组和对照组之间的差异有统计学意义．PSI1(7)与子痫前期的相关系数rs为0.780，即两者存在较高的相关性．对PSI1(7)的均值和标准差做进一步分析，如图7所示．图中实验组的平均值明显大于对照组，与前面分析结果相符．同时，相比于子痫前期，正常孕妇数据的标准差相对更大，分布更加离散．将PSI1(7)的分析结果与医院诊断结果进行对照．50例数据中，真阳性20例，假阳性1例，真阴性26例，假阴性3例，即灵敏性为87.0%，特异性为96.3%，准确性为92.0%．可见，通过光电容积斜率指数PSI1(7)可有效地区分子痫前期患者与正常孕妇．(1) 通过光电容积脉搏波的形态学分析，根据下降支的斜率特征提取出了新的特征参数——光电容积斜率指数，用以表征血管壁弹性和外周血管阻力等生理状况．(2) 实验结果表明，通过光电容积斜率指数可有效地区分子痫前期患者和正常孕妇，为子痫前期的无创诊断研究提供了有力的理论依据和新的研究方法．(3) 由于子痫前期脉搏波缺乏标准数据库，文中所涉及的临床数据均是本次研究期间在合作医院实际采集得到，实验样本数据量较少．随着研究的推进，数据将不断补充，继而增强实验结果的可靠性，提高子痫前期与正常孕妇分类的正确率．【相关文献】［1］ American College of Gynecologists. ACOG practice bulletin No.33：Diagnosis and management of preeclampsia and eclampsia[J]. Obstetrics and Gynacology，2002，99(1)：159-167.［2］ Fabry I G，Richart T，Cheng X，et al. Diagnosis and treatment of hypertensive disorders during pregnancy[J]. Acta Clinica Belgica，2010，65(4)：229-236.［3］ Euliano T Y，Michalopoulos K，Singh S，et al. Photoplethysmography and heart rate variability for the diagnosis of preeclampsia[J]. Obstetric Anesthesiology，2018，126(3)：913-919.［4］ Lan Kunchan，Raknim P，Kao Weifong，et al. Toward hypertension prediction based on PPG-derived HRV signals：A feasibility study[J]. Journal of Medical System，2018，42(6)：103.［5］ Hangsik Shin，Se Dong Min. Feasibility study for the non-invasive blood pressure estimation based on PPG morphology：Normotensive subject study[J]. Biomedical Engineering Online，2017，16：10.［6］孙磊. 基于脉搏波传导时间变异性的低血压预测研究[D]. 杭州：浙江大学生物医学工程与仪器科学学院，2014.Sun Lei. Study of Hypertension Prediction with Pulse Transit Time Variability[D]. Hangzhou：School of Biomedical Engineering and Instrument Science，Zhejiang University，2014(in Chinese).［7］侯凌霄. 基于光电容积脉搏波的舒张期衰减时间参数的研究[D]. 杭州：浙江大学生物医学工程与仪器科学学院，2012.Hou Lingxiao. Study of Diastolic Decay Time Constant from PhotoplethysmographicPluse[D]. Hangzhou：School of Biomedical Engineering and Instrument Science，Zhejiang University，2012(in Chinese).［8］ Hou Lingxiao，Wang Xuan，Chen Xinzhong，et al. A novel non-iterative shape method for estimating the decay time constant of the finger photoplethysmographic pulse[J]. Journal of Zhejiang University：Science A，2011，12(6)：438-445.［9］王选. 基于光电容积脉搏波形态的伤害感受指数的研究[D]. 杭州：浙江大学生物医学工程与仪器科学学院，2012.Wang Xuan. The Surgical Pleth Modal Index for Monitoring Nociception[D]. Hangzhou：School of Biomedical Engineering and Instrument Science，Zhejiang University，2012(in Chinese).［10］ Feng Y，Chen X，Wang X，et al. Area difference ratio for assessing nociceptive balance during laryngoscopy and intubation under intravenous anaethesia：Preliminary investigation of a novel photoplethysmographic variable[J]. European Journal of Anaesthesiology，2015，32(1)：58-59.［11］ Wang X，Chen X，Ye S，et al. Changes in the photoplethysmogram with trachealintubation and remifentanil concentration[J]. Anaesthesia，2012，67(12)：1332-1336. ［12］刘静. 麻醉状态下心血管多参数监护系统的软件开发[D]. 杭州：浙江大学生物医学工程与仪器科学学院，2016.Liu Jing. Software Design for Cardiovascular Multi-parameter Monitoring System During Anesthsia[D]. Hangzhou：School of Biomedical Engineering and Instrument Science，Zhejiang University，2016(in Chinese).［13］付秀泉. 基于脉搏波分析的麻醉手术中无创生理参数检测系统[D]. 杭州：浙江大学生物医学工程与仪器科学学院，2012.Fu Xiuquan. A Non-invasive Physiological Parameter Monitoring System Based on Photoplethysmography in Anesthsia[D]. Hangzhou：School of Biomedical Engineering and Instrument Science，Zhejiang University，2012(in Chinese).。
浙江大学教学、科研用房管理办法（征求意见稿）第一条根据《浙江大学公用房管理条例》制订本实施办法。
第二条本实施办法所涉及的各项定编及收费标准是根据目前情况确定的，今后将根据国家及学校的具体情况，由校公用房工作领导小组研究后进行适当调整。
第三条各院（系）级单位可根据自身实际制订实施细则。
第四条本实施办法中“定额面积”为房屋建筑面积（含厕所、走廊等共用部位面积），单位为平方米（㎡）。
第五条本实施办法将学院（系）的“定额面积”分为“定编面积”和“核算面积”。
其中“定编面积”由实验室实习场所及附属用房、教师办公用房、行政办公用房、教学科研辅助用房组成，主要参照教育部“92招标”核定。
“核算面积”由高层次人才与重点研究基地补贴用房和科研业绩补贴用房组成，根据学校实际情况核定。
第六条学院（系）实验室实习场所及附属用房定编实验室实习场所及附属用房包括三年级以上本科生、研究生的专业基础课、专业课、科研项目以及博士后研究所需的各种实验室、实习工厂农场牧场林场、实验室附属用房（准备室、天平室、仪器室、标本室、模型室、陈列室、动物室、空调室、更衣室、实验人员办公室等）。
学院（系）实验室实习场所及附属用房定编计算公式为：学院（系）实验室实习场所及附属用房定编面积 = （三年级以上本科生人数+研究生人数）×M +研究生人数×B+博士后人数×H其中M为按照学科规模和类别额定的每个三年级以上本科生、研究生的用房面积, B 为每个研究生补助指标，H为每个博士后研究人员核定指标。
各学院M、B、H值如下表：注：（1）、外语科的指标包括语言实验室。
（2）艺术科的指标包括音乐和美术专业的各种琴房、演奏厅、画室、雕塑室、展览厅等。
（参照92指标核定。
按照3年级以上人数我校学科自然规模：工科大约15000人，理科大约3000人、农科大约3000人、医科大约3300人，人文、社会学科大约6500人）第七条学院（系）教师教学与科研的办公用房定编每位全职岗教师原则上只能配置一处办公用房；短期（6个月以下）交流教师原则上不配置专门办公用房，确有需要由学院（系）内部调整公用办公空间解决；项目聘用人员由聘用单位根据工作性质在项目用房中统筹考虑。
