一種高效實時脈搏血氧監(jiān)測系統(tǒng)的研究

張根選，石 波，劉勝洋，陳建方，TSAU Young

（1.蚌埠醫(yī)學(xué)院醫(yī)學(xué)影像學(xué)系，安徽蚌埠233030；2.深圳迪美泰數(shù)字醫(yī)學(xué)技術(shù)有限公司，廣東深圳518067）

一種高效實時脈搏血氧監(jiān)測系統(tǒng)的研究

張根選1，石 波1，劉勝洋1，陳建方1，TSAU Young2

（1.蚌埠醫(yī)學(xué)院醫(yī)學(xué)影像學(xué)系，安徽蚌埠233030；2.深圳迪美泰數(shù)字醫(yī)學(xué)技術(shù)有限公司，廣東深圳518067）

針對傳統(tǒng)的脈搏血氧監(jiān)測存在光電檢測電路復(fù)雜、實時性差等缺點(diǎn)，本文提出了一種高效實時的監(jiān)測方法。以光電容積脈搏波描記法（PPG）為基本原理，構(gòu)建了以光頻轉(zhuǎn)換技術(shù)為核心的全數(shù)字硬件系統(tǒng)，采用零相移數(shù)字濾波和小波分析技術(shù)實現(xiàn)脈搏血氧信號的預(yù)處理，運(yùn)用線性回歸分析完成血氧飽和度的計算。采用Fluke的Index 2血氧模擬儀產(chǎn)生的標(biāo)準(zhǔn)信號進(jìn)行測試，在70%～100%范圍，誤差小于2%。結(jié)果表明，該方法切實有效，具有電路簡單、數(shù)據(jù)利用率高等優(yōu)點(diǎn)，可以實現(xiàn)血氧值的秒級更新。

脈搏血氧飽和度；光電容積脈搏波描記法；光頻轉(zhuǎn)換；線性回歸分析

1 引 言

脈搏血氧計是一種無創(chuàng)連續(xù)監(jiān)測血氧飽和度（SpO2）的設(shè)備，現(xiàn)已廣泛應(yīng)用于手術(shù)室、恢復(fù)室、急救病房、危重病房等。現(xiàn)有的脈搏血氧計大多是以光電容積脈搏波描記法（Photoplethysmography，PPG）為基本原理［1－3］。在測量中，SpO2的計算是以特征值R的準(zhǔn)確提取為基礎(chǔ)的。傳統(tǒng)的R值提取方法大多采用脈搏波的峰谷值法［4－5］，不能算是一種嚴(yán)格意義上的實時監(jiān)測。此外，在這些設(shè)備中，脈搏血氧信號的硬件采集系統(tǒng)大多都是采用模擬電路［6－7］，主要由光電探測器、電流—電壓轉(zhuǎn)換電路、模擬放大器、模擬濾波器等部分組成。這種設(shè)計存在電路復(fù)雜、抗干擾能力差、內(nèi)部噪聲大、體積大、功耗高等各項不足。由此可見，脈搏血氧飽和度監(jiān)測作為發(fā)展中的一項技術(shù)，在數(shù)據(jù)采集和血氧信號處理等方面還需要進(jìn)一步改進(jìn)和完善。本文以PPG為基本原理，構(gòu)建了以光頻轉(zhuǎn)換技術(shù)為核心的全數(shù)字硬件系統(tǒng)，采用光頻轉(zhuǎn)換技術(shù)、零相位數(shù)字濾波和小波分析技術(shù)實現(xiàn)脈搏血氧信號的預(yù)處理以及線性回歸算法的使用實現(xiàn)了血氧的高效實時監(jiān)測。

2 方 法

2.1 測量原理

PPG是血氧飽和度的一種無創(chuàng)檢測方法，它是基于Lambert-Beer定律，利用氧合血紅蛋白（HbO2）和脫氧血紅蛋白（Hb）在紅光和紅外光區(qū)的光吸收系數(shù)差異來進(jìn)行的。基于PPG原理所測得的血氧飽和度一般用SpO2表示，計算公式可表示為：

其中，CHbO2和CHb分別表示動脈血液中HbO2和Hb的濃度。

假設(shè)波長為λ，光強(qiáng)為Io的單色光垂直入射指端時，透射光強(qiáng)可分為兩部分：一部分是直流成分，主要反映指端組織中各種非脈動成分，如肌肉、骨骼、色素、脂肪、水和靜脈血等對光的吸收；另一部分為交流成分，是由血管隨心臟搏動而產(chǎn)生的收縮和舒張引起，主要反映動脈血中HbO2和Hb對光的吸收。如果忽略散射、反射等因素造成的衰減，根據(jù)Lambert-Beer定律，透光強(qiáng)度I可以表示為：

其中，εo，Co和L為組織內(nèi)各種非脈動成分的總吸光系數(shù)、光吸收物質(zhì)的濃度和光路徑長度；εHbO2，εHb分別為動脈血液中HbO2和Hb的吸光系數(shù)。

當(dāng)動脈搏動時，入射光穿過指端的光路徑長度會隨之發(fā)生變化，假設(shè)變化為L＋ΔL，對應(yīng)的透射光強(qiáng)度變化為I＋ΔI，這時式（1）可以寫成：

由于透射光中交流成分占直流成分的百分比很小，即ΔI／I＜＜1，這時，將式（2）和式（1）相除，并對兩邊取對數(shù)，可寫成：

如果采用660 nm的紅光和940 nm紅外光作為光源，則用這兩種特定波長的光照射手指時，可以得到：

將（5）和（6）兩式相除可得：

式中，A和B為經(jīng)驗常數(shù)，可以通過實驗定標(biāo)得到；特征值R本文采用線性回歸算法來提取。

2.2 光頻全數(shù)字硬件系統(tǒng)的設(shè)計

硬件系統(tǒng)采用光頻轉(zhuǎn)換技術(shù)［8－9］，結(jié)構(gòu)如圖1所示。整個系統(tǒng)由電源模塊、微處理器模塊、電光模塊、光頻轉(zhuǎn)換模塊（TSL230RD）以及通信接口模塊等組成。

圖1 硬件系統(tǒng)框圖Fig.1 Hardware system block diagram

為了保證運(yùn)算速度，微處理器選擇微芯公司的PIC32MX4XXL。微處理器的I／O接口輸出脈沖控制光源驅(qū)動電路，其中光源為660 nm和940 nm的雙波長發(fā)光二極管。光頻轉(zhuǎn)換器選擇德州儀器的TSL230RD，將穿過人體組織吸收衰減后的透射光轉(zhuǎn)換為頻率信號。微處理器將脈沖頻率捕捉并進(jìn)行數(shù)字濾波預(yù)處理，然后完成血氧信號的計算。

理論上，紅光和紅外光強(qiáng)度信號必須是相同時刻采集才能準(zhǔn)確計算出血氧飽和度的值，但是技術(shù)上由于兩種光同時照射相互干擾而難以實現(xiàn)，實際上對兩種光強(qiáng)度信號的采樣為交替式而不是完全同步，即使發(fā)光二極管按照一定的時序交替發(fā)光。對這兩種波長的光強(qiáng)度信號采樣通常采用恒定采樣間隔，其時間差應(yīng)當(dāng)越小越好。

2.3 信號預(yù)處理

在血氧信號采集過程中，光頻轉(zhuǎn)換器接收到的光信號除了包含脈搏信息里的投射光信號，還包含測量環(huán)境下的背景光信號、運(yùn)動偽差、工頻和其他儀器如高頻電刀等產(chǎn)生的干擾等。傳統(tǒng)的設(shè)計都是采用模擬濾波或者是模擬濾波和數(shù)字濾波并用的方法。本文采用零相位數(shù)字濾波和小波分析去除噪聲信號。

2.3.1 零相位濾波去除血氧信號中的基線漂移

以紅光為例信號進(jìn)行分析（紅外光分析方法相同），紅光信號的采樣點(diǎn)數(shù)設(shè)為M＝3000，紅光原始信號如圖2和圖3所示，信號中除了有嚴(yán)重的基線漂移以外，還嚴(yán)重存在各種頻率的干擾信號。經(jīng)過零相位數(shù)字濾波去除基線漂移后信號描如圖4和圖5所示。

圖2 原始信號的整體波形Fig.2 overall signal waveform of the original signal

圖3 原始信號的局部波形Fig.3 local signal waveform of the original signal

圖4 去除基線漂移后的整體信號波形Fig.4 the whole waveform after remove baseline drift

圖5 去除基線漂移后的局部整體信號波形Fig.5 the localwaveform after remove baseline drift

2.3.2 小波分析去噪

根據(jù)小波變換的理論可知，由于小波變換具有多分辨率的特性，因此可先將其分解，然后重構(gòu)有用信號，經(jīng)過實驗數(shù)據(jù)分析，將a5與d5兩層信號重構(gòu)得到的信號效果最好。圖6和圖7為采用小波去噪方法對去除基線漂移后的血氧信號處理后的結(jié)果。通過比較小波消噪前后的波形圖，可見噪聲得到了明顯的去除，尤其是高頻噪聲。

圖6 去噪后的整體信號波形Fig.6 The overall signal after rejecting noise

2.4 血氧飽和度的實時計算

圖7 小波去噪后的局部信號波形Fig.7 The local signal after rejecting noise

線性回歸是利用數(shù)理統(tǒng)計中的回歸分析，來確定兩種或者兩種以上變量間相互依賴的定量關(guān)系的一種統(tǒng)計分析方法，應(yīng)用十分廣泛。對于生物組織這種強(qiáng)散射性、各向異性的復(fù)雜介質(zhì)，不完全符合Lambert-Beer定律，散射效應(yīng)一方面會改變一部分光的出射方向，另一方面使得光在介質(zhì)中的傳播路徑變長，被吸收的幾率增加。這些都會導(dǎo)致生物組織對入射光衰減的改變。如果兩束透射光強(qiáng)度之間的具有相關(guān)性，就可以使用連續(xù)采樣過程產(chǎn)生的大量樣本數(shù)據(jù)進(jìn)行線性回歸提取R值，從而大幅度提高精度［10］。

假設(shè)紅光和紅外光的透射光強(qiáng)信號序列分別為x＝｛x1，x2，x3，……，xN｝和y＝｛y1，y2，y3，……，yN｝。如果證明兩個信號序列之間存在相關(guān)性，則可以建立一元線性回歸模型，即：

其中，a，b為回歸系數(shù)，利用最小二乘法可以求出：

其中，N為參與運(yùn)算的樣本點(diǎn)的個數(shù)。經(jīng)過推倒可以得到R的表達(dá)式為：

由此可見，R值的提取僅取決于參與運(yùn)算的采樣點(diǎn)的個數(shù)而不依賴于峰谷值，只需要少量的普通采樣點(diǎn)就可以完成血氧值的計算，大大提高了數(shù)據(jù)使用效率，降低了硬件系統(tǒng)的負(fù)擔(dān)，可以實現(xiàn)了真正意義上的血氧飽和度實時測量。

3 結(jié) 果

對公式（9）中的兩個常數(shù)進(jìn)行定標(biāo)后，用線性回歸算法計算的R值計算出血氧飽和度。實驗采用Fluke的Index 2血氧模擬儀產(chǎn)生不同血氧飽和度結(jié)果與標(biāo)準(zhǔn)信號間的誤差如表1所示。

表1 血氧飽和度測試結(jié)果Tab.1 oxygen saturation of blood test results

在70%～100%范圍，誤差小于2%。

4 結(jié) 論

本文提出并驗證了一種高效實時的脈搏血氧監(jiān)測系統(tǒng)，該系統(tǒng)是運(yùn)用線性回歸算法實現(xiàn)數(shù)據(jù)利用率的高效性和血氧值更新的實時性；同時光頻轉(zhuǎn)換技術(shù)和數(shù)字濾波技術(shù)的應(yīng)用具有模擬濾波技術(shù)無與倫比的優(yōu)點(diǎn)。該系統(tǒng)可用于心電監(jiān)護(hù)儀中的血氧監(jiān)測模塊以及穿戴式血氧設(shè)備。

［1］ John A Crowe，Damianos Damianou.The wavelength dependence of the photoplethysmogram and its implication to pulse oximetry［C］.14thAnnual International Conference of the IEEE Engineering in Medicine and Biology Society，1992，6：2423－2424.

［2］ Webster JG.Design of pulse oximeters［M］.British：Institute of Physics Publishing，1997.

［3］ Luo Zhichang，Zhang Song，Yang Yimn.Engineering analysis for pulse wave and its application in clinical practice［M］. Beijing：Science Press，2006：165－179.（in Chinese）羅志昌，張松，楊益民.脈搏波的工程分析與臨床應(yīng)用［M］.北京：科學(xué)出版社，2006：165－179.

［4］ Reichelt S，F(xiàn)iala J，Werber A，et al.Development of an implantable pulse oximeter［J］.IEEE Transactions on Biomedical Engineering，2008，55（2）：581－588.

［5］ Hong Lei，Chang Jinyi，Guan Yujing.A-trous and immune genetic algorithm based resrarch on blood oxygen signal［J］.Journal of Changshu Institute Technology：Natural Sciences，2009，23（2）：106－110.（in Chinese）洪蕾，常晉義，關(guān)玉景.動態(tài)血氧飽和度監(jiān)測算法的研究［J］.常熟理工學(xué)院學(xué)報：自然科學(xué)版，2009，23（2）：106－110.

［6］ Wu Xiaoling，Cai Guiyan.The development of a wearable pulse oximeter sensor and study of the calibration method［J］.Journal of Biomedical Engineering，2009，26（4）：731－738.（in Chinese）吳小玲，蔡桂艷.可穿戴式血氧傳感器的設(shè)計和定標(biāo)方法研究［J］.生物醫(yī)學(xué)工程學(xué)雜志，2009，26（4）：731－738.

［7］ Xu Xiaofeng，Li Zhe，Ling Zhenbao.Design of blood oxygen saturation measure system based on reflection type［J］.Journal of Jinlin University：Information Science Edition，2013，31（3）：260－265.（in Chinese）許曉峰，歷哲，凌振寶.反射式血氧飽和度測量系統(tǒng)設(shè)計［J］.吉林大學(xué)學(xué)報：信息科學(xué)版，2013，31（3）：260－265.

［8］ Yang Tao，Lin Wanhua，Weng Yujie，etal.Design of digital pulse oximeter with high resolution based on opticalfrequency converter［J］.Beijing Biomedical Engineering，2011，30（1）：73－77.（in Chinese）楊濤，林宛華，翁羽潔，等.基于光－頻轉(zhuǎn)換器的高精度數(shù)字脈搏血氧儀的設(shè)計［J］.北京生物醫(yī)學(xué)工程，2011，30（1）：73－77.

［9］ Gong Yushun，Wu Baoming，Gao Dandan，et al.Design of an anti-noise wearable pulse oximetry［J］.Chinese Journal of Sensors and Actuators，2012，25（1）：6－10.（in Chinese）龔渝順，吳寶明，高丹丹，等.一種抗干擾穿戴式血氧飽和度監(jiān)測儀的研制［J］.傳感技術(shù)學(xué)報，2012，25（1）：6－10.

［10］Shi Bo，Liu Shengyang，Chen Jianfang，et al.Features extraction from pulse oximetry signals through linear regression methods：a preliminary experimental study［J］.Chinese Journal of Medical Physics，2013，30（1）：3913－3916.（in Chinese）石波，劉勝洋，陳建方，等.線性回歸算法提取脈搏血氧特征值初步驗證［J］.中國醫(yī)學(xué)物理學(xué)雜志，2013，30（1）：3913－3916.

Efficient real-time pulse oximetry monitoring system

ZHANG Gen-xuan1，SHIBo1，LIU Sheng-yang1，CHEN Jian-fang1，TSAU Young2
（1.Department of Medical Imaging，Bengbu Medical College，Bengbu 233030，China；2.Dimetek Digital Medical Technologies Ltd，Shenzhen 518067，China）

To solve the problems of the circuit complexity and poor real-time performance existed in traditional pulse oximetrymonitoring system，an efficient real-time method is proposed.Based on the photoplethysmograph（PPG），a digital hardware system is constructed with the light to frequency conversion technology.This system adopts zero phase digital filter and wavelet analysis technology to realize pulse oximetry preprocessing，and employs linear regression algorithm to complete the pulse oxygen saturation（SpO2）calculation.By using thismethod，the standard signals in the range of70%～100%derived from Fluke Index 2 XL oxygen simulator have acquired error rate less than 2%.The results indicate that，this efficientmethod has the advantages of low complexity of circuit and high utilization rate of data，and can refresh the SpO2value at the level of second.

pulse oxygen；saturation（SpO2）；PPG；optical frequency conversion；linear regression analysis

TN29

A

10.3969／j.issn.1001-5078.2014.02.0

1001-5078（2014）02-0187-05

國家級大學(xué)生創(chuàng)新訓(xùn)練資助項目（No.201210367013）資助。

張根選（1976－），男，碩士，講師，主要從事生物醫(yī)學(xué)信號檢測、生物醫(yī)學(xué)儀器研究。E-mail：cyzgx123＠163.com

2013-07-04；

2013-09-12