《基于壓力波的人體脈搏波傳播速度無創(chuàng)檢測研究》由會員分享，可在線閱讀，更多相關《基于壓力波的人體脈搏波傳播速度無創(chuàng)檢測研究（13頁珍藏版）》請在裝配圖網上搜索。

1、 基于壓力波的人體脈搏波傳播速度無創(chuàng)檢測研究 26卷3期 2007年6月 中國生物醫(yī)學工程 ChineseJournalofBiomedicalEoain*ering V01.26No.3 June2007 基于壓力波的人體脈搏波傳播速度無創(chuàng)檢測研究 李婷虞鋼 (中國科學院力學研究所,北京1OOO80) 摘要:脈搏波傳播速度的變化往往是血管彈性程度發(fā)生變化的最早體征,因此脈搏波速度的檢測對動脈硬化等 疾病的診斷有較高的參考價值,其檢測精度主要取決于波速參考點的提取算法.在分析實際采集的脈搏壓力波信 號特征的基礎上,對波峰,波足等特征點進行了識別,并提出了一種根據脈搏波

2、前沿內斜率最大值點的個數和位置 確定惟一波速參考點的新算法.臨床實驗證明,新方法能夠有效地提取脈搏波速度,有助于動脈血管彈性程度的 輔助診斷. 關鍵詞:脈搏波速度;波速參考點;波速分析 ResearchofMethodforPulseWaveVelocityAnalysis BasedonPressureWaveDetection LITing,YUGang (InstituteofMechanics,ChineseAcademyofSc/encez,Being100080) Abstract:Pulsewavevelocity(PWV)isaprominentmarkerof

3、arterialstiffnessandisanvaluablefactortodiagnose arteriosclerosis.111eprecisionofthePWVdetectionisdeterminedbythePⅣVreferencepointextractionalgorithm.In thisworkthewavepeaksandfeetweredistinguishedafterthewavecharacterswereinvestigatedindetail.Anew algorithmdeterminingtheuniquePWVreferencepointby

4、thenumberandpositionsofcurrentmaximumgradient pointswaspmpesed.ExperimentalresultsshowedthatthepresentedalgorithmcouldbeusedtocalculatethePⅣV accuratelyanditisprovedtobeagoodassistantmethodtoassessbloodvesselstiffness. Keywords:pulsewavevelocity;PWVreferencepoint;analysis. 中圖分類號R318.04;文獻標識碼A文章編

5、號0258.8021(2007)03.0321.05 引言 在心血管系統(tǒng)中,心臟周期性射血引起的血液 壓力脈動和血管壁振動在動脈血管系中的傳播稱為 脈搏波在動脈中的傳播.人體脈搏波傳播速度 (PWV)的大小在很大程度上依賴于血管的物理特性 和幾何特性,特別是與動脈血管壁的彈性特性密切 相關….研究發(fā)現(xiàn),血管彈性程度越好,Pwv越小, 隨著年齡的增長或其他原因引起動脈硬化時,動脈 管壁發(fā)生的變異將使其彈性模量增大,從而使PWV 增大.,與舒張壓和收縮壓相比,與心血管病 發(fā)病率和死亡率的相關性更大j.因此,在心 收稿日期:2006-07.03,修回日期:2007.04-12

6、. *通訊作者.E-mail:gyu@imeeh.aC.Cll 血管系統(tǒng)評價中得到了廣泛的應用,目前已經成為 評價動脈血管彈性程度的重要標準. 二十世紀2O年代以來,國際上在PWV檢測方 面做了大量的研究,并有相應的產品問世J,國內在 這方面的研究工作還比較少.目前,國際上檢測 PwV的方法主要分為基于壓力波的檢測方法,基于 超聲多普勒的檢測方法和基于核磁共振的檢測方 法.超聲法和核磁共振法的計算量大,儀器成本高; 基于壓力波的檢測方法計算量少,儀器成本相對較 低,其檢測精度取決于波速參考點的選擇和提取. 在基于壓力波的檢測中,一般以每個波的起始 點或斜率最大值點作為

7、波速參考點.事實上,對于 322中國生物醫(yī)學工程26卷 一 些皮膚松弛的老年人,波形起始點往往容易被前 一 舒張期末的噪聲淹沒而不易提取;另外,對應于一 個脈搏波通常有多個斜率最大值點,波速參考點只 需要一個,而相關文獻中并沒有提到如何得到最終 參考點.針對這個問題,本研究提出了一種在脈搏 波前沿升支中由斜率最大值點的個數和位置確定惟 一 波速參考點的算法,實驗證明,該算法抗干擾能力 強,穩(wěn)定性好. 1脈搏信號檢測平臺 在臨床上,通常采集人體某些骨骼附近體表動 脈處的脈搏波,如橈動脈,頸動脈,踝動脈等,主要由 于這些部位容易暴露,組織結構好(表層皮膚和肌肉

8、 都較薄,血管下方為堅硬的骨頭,檢測時血管不易滑 動,因此容易得到精確的脈搏信號. 為實現(xiàn)脈搏信號的檢測,自行研制了脈搏信號 檢測平臺,其結構見圖1. 數據處I1個人信息ll脈搏波 理模塊ll數據庫lI數據庫 人機接口 圖1脈博信號檢測平臺結構圖 Fig.1Theframeworkofthepulsewavedetectionsystem 脈博信號檢測平臺由生理壓力傳感器,信號預 處理電路,A/D轉換電路,人機接口,標準通信接口, 數據處理模塊和數據構成.其工作原理如下:壓力 傳感器同時采集兩路脈搏信號,放大,濾波,模數轉 換后由標準通信接口送往計算機,計算機實現(xiàn)脈搏

9、 信號和患者基本信息的數據庫存儲,并通過鍵盤,鼠 標,顯示設備和打印設備實現(xiàn)人機交互,計算機內配 有相應的脈搏分析軟件系統(tǒng),實現(xiàn)數據的處理.采 集到的脈搏信號見圖2. 2脈搏波傳播速度分析 脈搏波速度Pwv定義為兩脈搏測量點間的體 表距離Ax與脈搏傳播時間△f的比值,即 圖2采集得到的脈搏信號 Fig.2Thepulsewavedetected PiVV=(1) 計算脈搏速度的關鍵在于脈搏波傳播時間△f 的計算.同時測得的相隔一定距離的兩路脈搏波對 應參考點間的時間差即為脈搏波的傳播時間△f. 2.1波速參考點提取 提取波速參考點前,首先對采集到的脈搏信號 進行

10、中值濾波,濾除信號中的脈沖干擾,然后提取每 個脈搏波的波峰點和波足點(如圖2所示),將波速 參考點的搜索范圍限定在波足和波峰之間,即波形 的前沿升支范圍內,最后根據該段信號中斜率最大 值點的個數和位置確定最終波速參考點. 2.1.1中值濾波 信號中的脈沖干擾極易被識別為斜率最大值 點,而中值濾波可以有效地抑止脈沖噪聲.因此,在 識別波形特征點前,首先對信號進行中值濾波.在 每1.5ms采樣一次的情況下,濾波器長度取11,即 將每個數據點的值用其對稱鄰域內11個點中大小 居于中間的值代替. 2.1.2波峰點提取 統(tǒng)計分析采集到的脈搏波信號發(fā)現(xiàn),雖然每個 脈搏波的波峰值

11、大小并不相同,但它們總是在一個 范圍內波動的,波動范圍基本不會超過最大波形高 度的0.3倍,因此可以認為波峰是在整段波形的最 大值附近,大于其領域內所有點的點.通過對大量 脈搏波測量數據的分析,認為把此鄰域長度定義為 90個數據點可以有效地抑制波值附近的小波動,準 確地識別波峰點.假設采集到的脈搏波為數字信號 序列p[j],其中jE[0,N],則實際波峰P[]的判 定條件為 f(max{P[]}一P[])≤△,, 【p[j]≥P[j],∈[J一45,J+45] 式中,△:(max{P[}一min{P[})×0.3,為波 峰點所在位置. 信號序列的最大值max{P[門}

12、和最小值 min{P[]}通過搜索比較整個序列中的所有信號點 獲得. 3期李婷等:基于壓力波的人體脈搏波傳播速度無創(chuàng)檢測研究323 2.1.3波足點提取 從統(tǒng)計分析采集到的脈搏波信號發(fā)現(xiàn),雖然每 個脈搏波的波足值大小并不相同,但它們總是在一 個范圍內波動的,波動范圍基本不會超過最大波形 高度的0.3倍,因此可以認為波足是在整段波形的 最小值附近,小于其領域內所有點的點.通過對大 量脈搏波測量數據的分析,認為把此領域長度定義 為90個數據點可以有效地抑制波足附近的小波動, 準確地識別波足點.假設采集到的脈搏波為數字信 號序列pEj],其中_『∈[0,N],則實際波足

13、pEj]的 判定條件為 f(WEde]一min{PEj]})≤A… 【PEj]≤PEj],J∈[J一45,J+45] 式中,A=(ma】({P[J]}一min{P[J]})×0.3,J為波 足點所在位置. 在具體實現(xiàn)上,可以利用特征點之間的關系,簡 化計算.考慮到同一個脈搏波的波足一定出現(xiàn)在波 峰的前部,而波峰識別先于波足識別,因此在識別波 足時可以利用已經識別所得的波峰,從波峰點開始, 向前搜索,將搜索得到的第一個滿足波足條件的數 據點確定為波足點.利用波峰輔助搜索波足算法的 時問僅為不利用波峰直接搜索波形算法的1/3. 2.1.4波速參考點確定 斜率最大值點的范

14、圍限定在脈搏波的前沿升支 范圍內,因此對于每一個脈搏波只需從波足和波峰 之間的數據點中搜索斜率最大的點即可.在實際判 定中,第i個脈搏波上升沿中斜率最大值點 P[～..]的判定準則為 P[～,]一P[一+1]=Illax{P[一P[_,+1]} ∈[fl')-(4) 式中,_『…為第個脈搏波上升沿中斜率最大值點 的位置,為第i個脈搏波的波足位置,為第i個 脈搏波的波峰位置. 實際提取所得每個脈搏波上升沿中斜率最大值 點一般有幾個到十幾個,而計算脈搏波速度的參考 點只需一個,必須根據已提取的多個斜率最大值點, 最終確定一個惟一的波速參考點.若以～[1], ～[2],…

15、,～[n]表示第i個脈搏波上升沿內 斜率最大值點的位置,…[.j}]表示其中第k個斜 率最大值點的位置,k∈[1,n],確定最終參考點的 具體做法是: 當n≤8時,以數據序列的位置重心～:寺 ∑～. [k]作為最終參考點位置; 當n>8時,采用拉依達判定準則確定最終參 考點位置. (1)計算～[1],～. [2],…,～[n]的平 均值 一一 . :i/∑J～,.[](5),一一.i～,sLJL (2)計算各點與平均位置點偏差的絕對值 P[]:1….一….[(6) (3)計算～[1],～,.[2],…,….[n]序列 的標準偏差 (7) (4)剔除

16、偏差過大的點 若測量數據的誤差符合正態(tài)分布,則數據落入 Pk]≤2a范圍的概率為95.4%,因此可以認為超 出此范圍的點誤差過大,應去除.具體作法為:若 P[k]≥2口去除～[k],回到第一步重新計算, 直到剩下的所有數據點滿足P.[.j}]≤2a為止. (5)將最終_『～值作為參考點的位置 2.2脈搏速度分析 對于同時測得的兩路脈搏波信號,若第一路脈 搏波中完整波形的個數為n.,第二路脈搏中完整波 形的個數為n,兩路波形中每一對完整脈搏波都對 應著一個時間差,因此對應于完整波形對的時間差 個數為n:min(n.,n),實際測得的這n個時間差 并不相同,而真實的脈搏波

17、傳播時間只有一個.這 些測量值只是真實傳播時間的一個反映,它們受到 被測者情緒或者肌肉和血管的當時狀態(tài)的影響而產 生波動.考慮到它們中的大部分只是在一個小范圍 內波動,只是偶爾有一些波動較大的時間差,因此對 這n個時間差進行處理,去掉誤差較大的值,最終 得到一個最接近真實脈搏傳播時間的值. 若第一對脈搏波的時間差為△￡.,第二對脈搏 波的時間差為△￡,…,第n對脈搏波的時間差為 △￡,最終△￡的確定方式如下: (1)求序列△￡.,△￡2,…,△￡的平均值△￡ A￡: 奎(8)n. 中國生物醫(yī)學工程26卷 (2)求序列中每個時間差At對于均值At的偏 差的絕對值

18、 P[f]=I△t一△tI(9) (3)求序列△t.,△t,…,△t的標準偏差 (10) (4)剔除誤差較大的時間差 若IP[f]I>20"則認為偏差過大,剔除△t,回 到第一步重新計算 若IP[i]I≤2a則進入下一步. (5)確定最終的At, 若最終剩余時間差序列為△t.,△t,…,△t,則 ● Atf=1∑Ati—td(11) .'i=1 式中,t為兩路信號通過系統(tǒng)的延遲時間差. 3實驗驗證 目前,醫(yī)學上對于脈搏波速度并沒有一個統(tǒng)一 的定義,不同檢測方法對脈搏波速度的定義和算法 各不相同,使得目前市面上各種儀器測得的脈搏波 速度之間可比性不大.在

19、這種情況下,穩(wěn)定性和測 量的可重復性便成為了評價檢測方法的重要指標. 為驗證本研究提出的檢測方法對檢測的可重 復性,進行了可重復性實驗.實驗對象為166個志 愿者,其中男84人,女82人,在相同的實驗條件下 對每個樣本進行4次PwV測量.實驗證明,即使波 足特征不明顯,通過本方法也可以很好地提取波速 參考點,并不影響irv的計算精度. 定義4次測量值的標準偏差與均值之比為偏差 率,越小說明測量的可重復性越好.對實驗結果的 偏差率進行了統(tǒng)計分析,分為三個年齡組: (1)20～59歲年齡組,共46人,其中男19人,女 27人; (2)60～69歲年齡組,共56人,其中男2

20、9入,女 27人; (3)70～89歲年齡組,共64人,其中男36人,女 28人. 各年齡組Ir,r偏差率的具體分布情況見圖3 ～ 圖5,總體偏差率的分布見圖6. 從檢測結果(圖3～圖5)可以看出,各年齡組 岡Irv檢測的偏差率分布趨勢相同,主要分布在0.03 以下,其次是0.03～0.05;從總體檢測情況來看, 岡Irv測量的總體偏差率主要分布在0.03以下,偏 06o0 0.5o0 040o 0.300 O.2oo 0lo0 0O00 <0030.03~o050.05~oI pwv偏差事 圖3PWV檢測偏差率分布(20—59歲年齡組J Fig

21、.3ThedistributionofPWVdeviationratio(20—59 agegroupJ <003003~o05005-0l 'pwv偏差事 圖4PWV檢測偏差率分布(6o一躑歲年齡組) Fig.4ThedistributionofPWVdeviationratio(6O一69 agegroupJ O7o0 06o0 O5o0 04o0 零 筮O3o0 02o0 O10o 00o0 <003003-005005-0l PWV偏差事 圖5PWV檢測偏差率分布(70—89歲年齡組J Fig.5ThedistributionofPⅥr

22、,,deviationratio(70—89 agegroup) 差率在0.05以下的樣本數約占總樣本數的90%. 由于人體是一個相當復雜的時變系統(tǒng),不同時刻測 得的Irv值不僅與血管壁本身的物理性質有關,還 與被測者當時的情緒以及血管壁周圍組織,神經,肌 肉的狀態(tài)有關,這使得PWV不是一成不變,而是在 咖啪徹徹枷姍瑚咖 OOOOOOOOOO 翠筮 3期李婷等:基于壓力波的人體脈搏波傳播速度無創(chuàng)檢測研究 07o0 06O0 05o0 墨04o0 03o0 02o0 0lo0 0Oo0 <003003~o05005--4)l PWV偏差率 圖6

23、PWV檢測的總體偏差率分布I男84人.女82人) Fig.6ThetotaldistributionofP,)'r,,deviationratio(84 menand82women) 一 個小范圍內波動,因此,本方法對PWv的檢測結 果比較真實地反映了血管的彈性程度,且檢測的可 重復性基本滿足臨床檢測的要求. 4結論 PWv是評價血管彈性程度的重要指標,本研究 提出了根據脈搏波前沿升支內斜率最大值點的個數 和位置確定惟一波速參考點的算法.實驗證明,該 算法在波足特征不明顯的情況下也能很好地確定惟 一 波速參考點,PWV計算結果能夠反映血管的彈性 程度,基本滿足臨床檢

24、測的需要. 參考文獻 柳兆榮,李惜惜.血液動力學原理和方法[M].上海:復旦大 學出版社.1997.179—180. J0hnl(J,LiYingzhu,DoreneO'Hara.Highpulsewavevelocityin hypertensivesareassociatedwitllincreasedvascalarstiffnessand waverdloctions[J].AJH,1998,11(4):161A. WomersleyJ.Mathematicaltheoryofoscillatingflowinallelastic tube[J].JPhysiol,19

25、55,l(2):37—8p. TineWH,JeorgenJ.SusanneR.et.Pulsewavevelocityand cardlovascuL~rdiseaseinageneralpopulation[J].AJH,2005.18 (5):A14. BarnesW,SmithD,kqJ,eta/.Non—involvecharacterizationof arterialsystemfunction[J].^J.II,2OO4.17(5):$56. Pre88ⅡlalIG,NewgardP.Ati"a/18~ucerforthecontinuousexternal ~tn*ementofarterialbloodpre~ure[J].IEEETransBiomedEng, 1963.10(4):73—81. KandaT,Nak~L,nuraE,MoritaniT,eta/.Arterialpulsewavevelocity andriskfactorsforperipheralvasculardisease[J].EurJApp] Physiol,2OOO.82(1—2):1—7. ]]]]]]] 1234567 [[[[[[[