李思琦，申華,2，朱勇，梁舜，覃開蓉

1.大連理工大學 電信學部生物醫(yī)學工程系，遼寧 大連 116024；2.大連東軟信息學院 電子工程系，遼寧 大連 116023

一種評估肱動脈內(nèi)皮功能的電阻抗信號檢測與分析系統(tǒng)

李思琦1，申華1,2，朱勇1，梁舜1，覃開蓉1

1.大連理工大學 電信學部生物醫(yī)學工程系，遼寧 大連 116024；2.大連東軟信息學院 電子工程系，遼寧 大連 116023

動脈內(nèi)皮功能受損或者功能障礙是動脈硬化早期的主要表現(xiàn)。定量評估肱動脈內(nèi)皮功能有助于及時發(fā)現(xiàn)心血管功能異常，從而為心腦血管疾病的早期診斷和治療提供一定的依據(jù)。目前臨床上主要采用彩色超聲多普勒技術檢測肱動脈血流介導擴張（FMD）來評估肱動脈內(nèi)皮功能。雖然這種檢測方法精確度與分辨率更高，但設備相對昂貴，操作復雜。本文提出一種基于前臂電阻抗信號的肱動脈內(nèi)皮功能檢測與分析系統(tǒng)。該系統(tǒng)通過檢測FMD實施前后前臂的生物電阻抗信號，并運用小波分析技術對前臂生物電阻抗信號的特征參數(shù)進行分析計算，包括特征頻率、絕對能量、平均幅值，以及FMD實施前后這些特征量的相對變化率。由于FMD實施前后前臂電阻抗信號的變化主要由前臂動脈管徑的變化引起，因此，本文提出的前臂電阻抗信號檢測與分析系統(tǒng)為定量評估肱動脈內(nèi)皮功能提供了一種新思路。與目前臨床上使用的內(nèi)皮功能檢測方法相比，本文提出的電阻抗信號檢測與分析系統(tǒng)具有操作簡單、經(jīng)濟實用、快速便捷的特點，不僅能用于家庭或社區(qū)醫(yī)療的監(jiān)測，并且可以實現(xiàn)在線檢測，為進一步研制智能可穿戴電阻抗動脈功能監(jiān)測儀器提供了一定的方法學基礎。

肱動脈內(nèi)皮功能；生物電阻抗信號；小波分析；特征頻率；絕對能量；平均幅值

1 研究背景

近年來，因心腦血管疾病發(fā)病和死亡人數(shù)急劇增多，治療費用也急劇增加，給患者家庭和社會帶來了沉重的負擔。動脈硬化作為一種全身性的病理變化過程，會直接導致重要臟器的損壞，已成為心腦血管疾病的獨立危險因子之一[1-6]。動脈內(nèi)皮功能受損或者功能障礙是動脈硬化早期的主要表現(xiàn)[7]。

動脈內(nèi)皮作為血管流動和血管壁之間的功能性屏障，不僅維持血管結構的完整，還在調(diào)節(jié)血管細胞生長、調(diào)節(jié)抗凝及纖溶系統(tǒng)、介導炎癥與免疫、調(diào)節(jié)白細胞與血小板在血管內(nèi)皮粘附、調(diào)節(jié)脂質(zhì)氧化、調(diào)節(jié)血管通透性等方面扮演十分重要的角色[7-9]。如果動脈內(nèi)皮功能受損或者功能障礙，血液中的脂肪會進入血管壁，與中膜的平滑肌細胞相互作用，最終在動脈壁形成動脈粥樣硬化斑塊[10]。因此，正確評價動脈內(nèi)皮功能、及時發(fā)現(xiàn)動脈內(nèi)皮功能異常，對于心腦血管疾病的早期預警起著至關重要的作用[11-12]。

肱動脈不僅是測量血壓的部位，也是檢測全身動脈功能的一個重要“窗口”。Celermajer等[13]使用超聲多普勒技術無創(chuàng)傷檢測肱動脈血流介導性擴張（Flow-Mediated Dilation，F(xiàn)MD）引起的動脈管徑擴張率評估動脈內(nèi)皮功能。隨后的一系列研究[11,14]表明通過FMD評估的肱動脈內(nèi)皮功能與心臟冠狀動脈的功能緊密相關。目前，利用高空間分辨率的超聲多普勒技術測量肱動脈FMD評估內(nèi)皮功能被美國心臟學會認為是一項“有前途的技術”，已成為臨床診斷的“金標準”[8]。然而，用該方法檢測動脈管徑的變化需要對動脈管徑進行精確定位，不僅操作程序煩瑣、耗時，而且相關的儀器設備笨重、價格昂貴，不便于推廣到社區(qū)或者家庭簡單使用。

生物電阻抗（Electrical Bio-Impedance，EBI）技術[15]是利用生物組織與器官的電特性及其變化規(guī)律提取與人體生理、病理狀況相關的生物醫(yī)學信息的檢測技術，基本測量方式是在體外（表）向檢測對象施加安全的電磁激勵（電流、電壓、電磁場），在體外（表）檢測相應電阻抗的變化獲取相關的信息。其特點是無創(chuàng)、無害，廉價、操作簡單和功能信息豐富，盡管精度較低，但易于被受試者接受。本文利用生物電阻抗技術測量肱動脈FMD實施前后前臂生物電阻抗信號，對檢測的生物電阻抗信號進行小波分析，提出了一種能夠分析肱動脈內(nèi)皮功能的簡便方法?；谠摲椒?，設計出一種具有操作簡單、價格低廉等優(yōu)點的動脈內(nèi)皮功能的電阻抗信號檢測分析裝置。

2 原理與方法

用氣袖無創(chuàng)傷阻斷肱動脈血流3～5 min后釋放，將引起下游前臂動脈的血流動力學環(huán)境特別是作用于動脈內(nèi)皮的壁面剪應力發(fā)生改變，改變的剪應力將誘發(fā)前臂動脈內(nèi)皮釋放血管舒張因子一氧化氮（NO），進一步引起前臂動脈的管徑發(fā)生改變，這一過程稱為FMD。在該過程中前臂動脈管徑的變化引起前臂動脈的總電阻抗發(fā)生改變。通過檢測靜息狀態(tài)和FMD實施后電阻抗信號的變化情況可以反映FMD引起的動脈管徑變化，從而達到間接評估動脈內(nèi)皮功能的目的。

前臂可看成一個電導體，其電阻抗由前臂內(nèi)所有動脈、靜脈、皮膚、肌肉、骨骼和其他組織的電導特性共同決定。由于FMD過程中只有動脈血流的電阻抗特性才發(fā)生明顯變化，因此所測量的前臂電阻抗的改變就是前臂動脈段電阻抗的改變。將前臂內(nèi)的動脈等效于一個縱向長度L不變的圓筒形腔室，其半徑R則在動脈正常搏動和FMD實施過程中發(fā)生改變（圖1），該動脈段的電阻抗Za表達為

式中，σ為血液電導率。

假定除動脈外的組織和器官總阻抗為常量Zp，則前臂總電阻抗Z為Zp和Za的并聯(lián)，表示為

進一步假定外圍組織和器官的阻抗Zp在FMD實施前后保持不變，則對公式(2)兩邊微分可得到

由于Zp＞＞Za，可近似為

即，前臂電阻抗Z的改變就是前臂內(nèi)動脈電阻抗Za的改變。另外，從公式(1)可以看出，前臂動脈電阻抗Za的變化主要取決于動脈內(nèi)半徑R的變化。

圖1 前臂動脈的等效圓筒模型

假定靜息狀態(tài)的前臂電阻抗信號為Z0(t)，F(xiàn)MD實施后的前臂電阻抗信號為ZF(t)，選擇Morlet母小波對連續(xù)的電阻抗波形Z0(t)和ZF(t)進行小波分析，可得到具有五個特征峰的頻譜圖，頻率區(qū)間為0.01～5.0 Hz，每個特征峰對應不同的頻率區(qū)間，反映了五種不同的生理活動[16-17]。其中，頻率區(qū)間0.01～0.1 Hz對應動脈內(nèi)皮的生理活動[16-17]，計算頻譜圖中該頻率區(qū)間的特征頻率fc，絕對能量Ec，平均幅值Ac。分別表達如下[18]

（1）特征頻率fc

則fc為頻率區(qū)間[fc1,fc2]內(nèi)的特征頻率。

（2）絕對能量Ec

（3）平均幅值Ac

其中，Ec(fc1,fc2)表示信號在頻率間隔[fc1,fc2]和時間間隔[t1,t2]內(nèi)的能量。Ac(fc1,fc2)表示在頻率間隔[fc1,fc2]和時間間隔[t1,t2]內(nèi)的平均小波幅值。

這樣，利用上述方法對FMD實施前的電阻抗信號Z0(t)進行分析，可得FMD實施前的3個特征參數(shù)：fc0，Ec0，Ac0，對FMD實施后的電阻抗信號ZF(t)進行類似分析，可得FMD實施后的另外3個特征參數(shù)：fcF，EcF，AcF，進一步，計算特征頻率、絕對能量、平均幅值在FMD實施前后的相對變化百分比，分別表示為

3 系統(tǒng)設計與實現(xiàn)

基于以上所述的原理和方法，本文設計和實現(xiàn)的評估肱動脈內(nèi)皮功能的電阻抗信號檢測與分析裝置，見圖2。該裝置包括檢測系統(tǒng)和分析系統(tǒng)兩部分（圖3），其中檢測系統(tǒng)主要用于檢測前臂生物電阻抗信號；而分析系統(tǒng)主要基于小波分析技術分析計算各項生物電阻抗信號特征參數(shù)[18]。圖3(a)用于檢測前臂生物電阻抗信號的測量裝置包括前臂上雙電極、前臂下雙電極、高頻低幅恒流源發(fā)生器、交流電壓檢測電路、信號放大器、濾波電路、A/D轉(zhuǎn)換電路和計算機。該測量裝置的前臂上雙電極置于待測手臂的上部，前臂下雙電極置于待測手臂的下部見圖3(b)；高頻低幅恒流源發(fā)生器和一對前臂上下電極構成電流回路；另一對前臂上下電極和交流電壓檢測電路構成回路檢測電壓變化。高頻低幅恒流源發(fā)生器給前臂上電極和前臂下電極加一個一定頻率ω的高頻低幅的恒流電源，使前臂通過一個強度恒定的交流電流I，這樣在前臂長度為L段兩端就建立起一定的交流電壓U。交流電壓檢測電路檢測到的交流電壓U經(jīng)放大器、濾波、A/D轉(zhuǎn)換將模擬的電壓信號轉(zhuǎn)換成數(shù)字信號輸入計算機，并根據(jù)如下公式(12)計算電阻抗Z隨時間的變化曲線：

由公式(12)計算的電阻抗信號存入計算機。

圖2 評估肱動脈內(nèi)皮功能的電阻抗信號檢測與分析裝置實物圖

圖3 系統(tǒng)設計與實現(xiàn)

分析系統(tǒng)主要由計算機、肱動脈內(nèi)皮功能參數(shù)計算模塊、顯示屏、硬盤等構成。將檢測系統(tǒng)輸出的電阻抗信號作為原始數(shù)據(jù)，并對其進行小波分解去噪，基于去噪后的前臂電阻抗信號，根據(jù)公式(5)～(11)分別計算描述肱動脈內(nèi)皮功能的信號特征參數(shù)，包括fc0，Ec0，Ac0，fcF，EcF，AcF，Δf，ΔE，ΔA。本系統(tǒng)的軟件部分采用C++實現(xiàn)，可實現(xiàn)病例創(chuàng)建、電阻抗波形顯示、性能選擇（包括量程、增益、采樣率和采樣時間的選擇）、以及分析計算結果的保存和打印。

4 實驗驗證

4.1 實驗方案

隨機選擇7名年齡22～23歲的健康大學生作為受試者，所有受試者血壓、心率正常，沒有心腦血管疾病。受試者在實驗前了解本實驗的詳細流程以及可能存在的風險，并簽署書面知情同意書。

整個實驗檢測過程中，將檢測裝置的電極置于待測前臂的上下兩端，上下電極距離1.5 cm，用常規(guī)水銀血壓計自帶的氣袖作為肱動脈加壓裝置。首先檢測并記錄FMD實施前靜息狀態(tài)前臂電阻抗信號Z0(t)；用血壓計自帶的氣袖給肱動脈加壓，阻斷肱動脈血流3分鐘后釋放，記錄釋放后的前臂電阻抗信號ZF(t)。

4.2 實驗結果與討論

圖4顯示的是利用本文裝置所檢測的受試者1在FMD實施前后一段前臂電阻抗信號的對比。FMD實施后前臂電阻抗信號的幅值小于FMD實施前的幅值，根據(jù)公式(1)可知，F(xiàn)MD實施后前臂動脈管徑增大，這符合血流介導的肱動脈擴張現(xiàn)象。圖4所示的電阻抗信號進行小波變換，得到的受試者1在FMD實施前后電阻抗信號頻譜圖的比較見圖5。

圖4 受試者No. 1在FMD實施前后一段前臂電阻抗信號的對比

圖5 受試者No.1在FMD實施前后電阻抗信號頻譜圖的比較

研究表明[16-17]，人體血流信號的頻譜可劃分出五個頻帶，每個頻帶代表不同的生理活動見表1。觀察圖5可發(fā)現(xiàn)，頻帶II、III和IV的特征峰不夠明顯，說明神經(jīng)活動、肌源性活動和呼吸活動對應的前臂電阻抗信號頻譜特征峰不明顯，造成這種現(xiàn)象的原因還有待于進一步探究；而頻帶I和頻帶V的特征峰明顯，且FMD實施前后特征峰值可產(chǎn)生明顯變化，說明內(nèi)皮活動和心臟活動對應的前臂電阻抗信號頻譜特征峰明顯，且在FMD實施前后變化明顯。鑒于本文系統(tǒng)主要用于評估肱動脈內(nèi)皮功能，因此只對圖5中頻帶I對應的內(nèi)皮細胞活動進行分析。

表1 人體血流信號頻譜的頻帶劃分[16-17]

表2列出了7名受試者在FMD實施前后前電阻抗信號9個特征參量，包括FMD實施前的特征頻率fc0、絕對能量Ec0和平均幅值Ac0，F(xiàn)MD實施后的特征頻率fcF、絕對能量EcF和平均幅值AcF，以及特征頻率、絕對能量和平均幅值在FMD實施前后的相對變化百分比Δf、ΔE和ΔA。從表2中可以看出，在FMD實施前后，特征頻率的相對變化絕對值在3.16%以內(nèi)；FMD實施后絕對能量和平均幅值均減小。但由于個體差異，不同的受試者性別、年齡、身高、體重不同，導致其身體狀況、健康程度、生理功能都不同，因此動脈內(nèi)皮細胞活動對應的絕對能量和平均幅值的相對減小量存在較大的差異。

5 結論與展望

正確評價動脈內(nèi)皮功能心腦血管疾病的預警起著至關重要的作用，本文利用生物電阻抗技術測量肱動脈FMD實施前后前臂生物電阻抗信號，并基于小波分析的方法設計出一種具有操作簡單、價格低廉等優(yōu)點的動脈內(nèi)皮功能的電阻抗信號檢測分析裝置。

針對測量的前臂電阻抗信號進行頻譜分析，可以提取到5個頻段的特征參數(shù)，分別是內(nèi)皮細胞活動、神經(jīng)活動、肌原性活動、呼吸活動和心臟活動，其中與肱動脈內(nèi)皮功能相關的是內(nèi)皮細胞活動。通過對比計算FMD實施前后的特征參數(shù)變化量，發(fā)現(xiàn)FMD實施后特征頻率基本不變，絕對能量、平均幅值均發(fā)生了不同程度的減小。由于FMD實施前后前臂電阻抗信號的變化主要由前臂動脈管徑的變化引起，因此，可以認為前臂電阻抗信號的變化反映了前臂動脈管徑的變化，即內(nèi)皮的功能。本文提出的前臂電阻抗信號檢測與分析系統(tǒng)為定量評估肱動脈內(nèi)皮功能提供了一種新思路。盡管本文初步的研究結果表明，F(xiàn)MD實施前后前臂電阻抗信號的這些特征參數(shù)發(fā)生了不同程度的變化，然而，我們?nèi)孕枰獙Ω嗟臉颖具M行測量與分析，通過對比健康對照人群與有動脈功能受損的疾病患者的前臂生物電阻抗信號，獲得健康對照人群的各個特征參數(shù)變化量和不同患者的特征參數(shù)變化量，以便建立評估規(guī)則，為進一步臨床評估肱動脈內(nèi)皮功能提供參考。

本文提出的檢測與分析系統(tǒng)具有操作簡單、價格低廉等優(yōu)點，可應用于家庭及社區(qū)的醫(yī)療檢測。隨著智能手機與可穿戴醫(yī)療設備的快速發(fā)展和廣泛應用，本文研究也為進一步研制智能可穿戴電阻抗動脈功能監(jiān)測儀器提供了一定的方法學基礎。

表2 FMD實施前后前臂電阻抗信號特征量的比較

[1]Laurent S,Boutouyrie P,Asmar R,et al.Aortic stiffness is an independent predictor of all-cause and cardiovascular mortality in hypertensive patients[J].Hypertension,2001,37(5):1236-1241.

[2]O’rourke MF,Mancia G.Arterial stiffness[J].J Hypertens,1999,17(1):1-4.

[3]Iannuzzi A,Licenziati MR,Acampora C,et al.Carotid artery stiffness in obese children with the metabolic syndrome[J].Am J Cardiol,2006,97(4):528-531.

[4]Favero G,Paganelli C,Buffoli B,et al.Endothelium and its alterations in cardiovascular diseases:life style intervention[J].Biomed Res Int,2014,2014:801896.

[5]Ramos-Casals M.The ill endothelium:how atherosclerosis begins in lupus[J].Rheumatology,2014,53(4):583-585.

[6]呂輝洋,劉璐,劉青,等.頸動脈粥樣硬化在動脈粥樣硬化性心血管疾病的研究新進展[J].心血管病防治知識:學術版,2015,(2):136-139.

[7]Ross R,Glomset JA.Atherosclerosis and the Arterial Smooth Muscle Cell[J].Science,1973,180(4093):1332-1339.

[8]Anderson TJ.Assessment and treatment of endothelial dysfunction in humans[J].J Am Coll Cardiol,1999,34(3):631-638.

[9]王慧群,陳愛華.血管內(nèi)皮與動脈粥樣硬化[J].醫(yī)學綜述,2008,14(4):598-600.

[10]林艾雯,陳竹君.動脈粥樣硬化與內(nèi)皮細胞損傷機制的研究進展[J].嶺南心血管病雜志,2015,21(4):580-582.

[11]Anderson TJ,Uehata A,Gerhard MD,et al.Close relation of endothelial function in the human coronary and peripheral circulations[J].J Am Coll Cardiol,1995,26(5):1235-1241.

[12]葛念嶧,馮新恒,謝道銀,等.冠心病血管內(nèi)皮功能的無創(chuàng)評估[J].臨床心血管病雜志,2000,16(11):491-494.

[13]Celermajer DS,Sorensen KE,Gooch VM,et al.Non-invasive detection of endothelial dysfunction in children and adults at risk of atherosclerosis[J].Lancet(London,England),1992,340 (8828):1111-1115.

[14]Takase B,Uehata A,Akima T,et al.Endothelium-dependent flow-mediated vasodilation in coronary and brachial arteries in suspected coronary artery disease[J].Am J Cardiol,1998,82 (12):1535-9,A7-8.

[15]董秀珍.生物電阻抗技術研究進展[J].中國醫(yī)學物理學雜志,2004,21(6):311-317.

[16]Benedi?i? M,Bernjak A,Stefanovska A,et al.Continuous wavelet transform of laser-Doppler signals from facial microcirculation reveals vasomotion asymmetry[J].Microvasc Res,2007,74(1):45-50.

[17]Stefanovska A,Bracic M,Kvernmo HD.Wavelet analysis of oscillations in the peripheral blood circulation measured by laser Doppler technique[J].IEEE Trans Biomed Eng,1999,46(10):1230-1239.

[18]覃開蓉,李思琦,申華,et al.一種肱動脈內(nèi)皮功能的電阻抗檢測分析裝置及方法[P].中國專利:201610092345.3,2016-02-19. (已公開)

Detecting and Analyzing System for Brachial Artery Endothelial Function

The damage or dysfunction of arterial endothelial function is the primary performance at the early stage of atherosclerosis. Quantitative evaluation of brachial artery endothelial function is helpful to discover the cardiovascular dysfunction in time,thus provides some evidence for the early diagnosis and treatment of cardiovascular and cerebrovascular diseases. To date,brachial artery endothelial function is evaluated at clinic by measuring the brachial fow-mediated vasodilation (FMD) with a Color Doppler Ultrasonography. Although the accuracy and resolution of this method are much higher,its equipment is relatively expensive,and its operation is very complicated. In this study,a detecting and analyzing system is proposed for evaluating brachial artery endothelial function based on forearm biological electrical impedance signal. By detecting forearm biological impedance signal before and after brachial FMD,the characteristic parameters of biological electro-impedance signals,including characteristic frequency,absolute energy,average amplitude and their relative alteration rate,are calculated using wavelet analysis technology. As the changes in the forearm biological electro-impedance before and after FMD are mainly caused by the changes in the brachial arterial diameter,the detecting and analyzing system of forearm biological electro-impedance signal would provide a new idea for evaluating quantitatively brachial artery endothelial function. In comparison with existing methods at clinic,the detecting and analyzing system of forearm biological electro-impedance signal is detection-simpler,cheaper,more practical and convenient,it can not only be used in the family or community medicine,but also realize on-line detection and provide with a methodology base for further development of the intelligent and wearable electrical impedance analyzer for monitoring arterial function.

arterial endothelial function;biological impedance signal;wavelet analysis;characteristic frequency;absolute energy;average amplitude

LI Si-qi1,SHEN Hua1,2,ZHU Yong1,LIANG Shun1,QIN Kai-rong1
1.Department of Biomedical Engineering,Faculty of Electronic Information and Electrical Engineering,Dalian University of Technology,Dalian Liaoning 116024,China;2.Department of Electronic Engineering,Dalian Neusoft University of Information,Dalian Liaoning 116023,China

TP274;R54

A

10.3969/j.issn.1674-1633.2016.09.003

1674-1633(2016)09-0009-05

2016-08-26

國家自然科學基金項目（31370948）；遼寧省公益研究科學基金項目（2014001027）。

覃開蓉，教授。

通訊作者郵箱：krqin@dlut.edu.cn