超聲粒子圖像測速技術(shù)評估左心室渦旋特性的初步研究

張耀楠，周偉，牛麗麗

1.西安思源學(xué)院 電子信息工程學(xué)院，陜西 西安 710038；2.東北大學(xué) 中荷生物醫(yī)學(xué)與信息工程學(xué)院， 遼寧 沈陽110169；3.中國科學(xué)院深圳先進(jìn)技術(shù)研究院 生物醫(yī)學(xué)與健康工程研究所， 廣東 深圳 518055

超聲粒子圖像測速技術(shù)評估左心室渦旋特性的初步研究

張耀楠1，2，周偉2,3，牛麗麗3

1.西安思源學(xué)院 電子信息工程學(xué)院，陜西 西安 710038；2.東北大學(xué) 中荷生物醫(yī)學(xué)與信息工程學(xué)院， 遼寧 沈陽110169；3.中國科學(xué)院深圳先進(jìn)技術(shù)研究院 生物醫(yī)學(xué)與健康工程研究所， 廣東 深圳 518055

左心室渦旋特性可被認(rèn)為是心臟健康潛在的指示器，通過超聲粒子圖像測速技術(shù)計算心臟內(nèi)渦旋特性是一個較新的研究方向，關(guān)鍵是要證明其具有足夠的可靠性和精確性，以便用于臨床應(yīng)用。為此，本文構(gòu)建了一個超聲粒子圖像測速技術(shù)系統(tǒng)，搭建了一個人體血液循環(huán)仿真系統(tǒng)，通過聚乙烯醇（PVA）制作了左心室仿體，以便在受控條件下測試超聲粒子圖像測速技術(shù)計算心臟內(nèi)渦旋特性的可行性。系統(tǒng)可以對仿體進(jìn)行連續(xù)成像，通過調(diào)整掃描線密度、圖像寬度、圖像深度等參數(shù)，最快可以實(shí)現(xiàn)每秒鐘上千幀的成像速度。實(shí)驗過程中，在液體水池中打入超聲造影微泡，開啟脈動泵循環(huán)多個周期，對仿體進(jìn)行實(shí)驗采集。通過超聲粒子圖像測速技術(shù)可得到心臟內(nèi)血液流動的速度分布，進(jìn)一步可計算出心臟內(nèi)渦旋的相關(guān)參數(shù)，包括能量損失（DI）、能量波動（EF）、渦旋相對強(qiáng)度（RS）、渦旋的波動（W）等。結(jié)果證明所構(gòu)建的系統(tǒng)、仿體、相關(guān)算法是可行的，可為下一步工作打好基礎(chǔ)。

超聲；粒子圖像測速；左心室；渦旋

0 引言

研究表明，左心室血液的流動動態(tài)可以被認(rèn)為是心臟健康一個潛在的指示器[1-2]。研究者在研究過程中發(fā)現(xiàn)，在心室舒張充盈的過程中，血流會在心室腔內(nèi)形成渦旋。在心室等容收縮期，心室的收縮會加大渦旋的強(qiáng)度，血液在心室內(nèi)以渦旋的形式流動，可以減少血液之間和血流與心室壁之間的碰撞。在心室的射血期，渦旋形式的血流有利于射血時血流方向的改變，從而提高心臟射血的效率。因此，渦旋可以最小化血液流動過程中能量的損失，并減少心臟射血時所需要的能量[3-9]。心臟相當(dāng)于一個大血泵，可為全身供血，在心臟收縮期時，心臟通過逆時針的扭轉(zhuǎn)和擠壓將血液泵至全身各處；在心臟舒張期時，通過順時針方向?qū)⒀撼榛匦呐K。當(dāng)心臟遭受心肌梗塞、心力衰竭、肥厚型心肌癥等疾病危害時，心臟的功能會發(fā)生紊亂，心臟內(nèi)的血流動力學(xué)參數(shù)會發(fā)生變化。

目前用于計算心臟內(nèi)渦旋參數(shù)的主要方法有：相位對比的心臟核磁共振成像技術(shù)和基于彩色多普勒的血流向量圖。1986年Nayler GL等人提出利用相位對比技術(shù)對血流速度進(jìn)行測量，相位對比的核磁共振成像（Phase-contrast Cardiac Magnetic Resonance）技術(shù)應(yīng)用大小相等、方向相反的兩葉組成的雙極梯度磁場，第一葉使靜止和流動組織的質(zhì)子群進(jìn)行相位積聚，隨之應(yīng)用第二葉，靜止質(zhì)子喪失相位，總相位為零，而流動的血液在兩個梯度之間，經(jīng)過不同的正負(fù)梯度，累加后產(chǎn)生相位位移。此時相位位移與血液在流速編碼方向的流動成正比，從而通過公式可計算得到心臟內(nèi)血流的速度場。近些年由于科學(xué)技術(shù)的發(fā)展，磁共振血流成像在三維成像的基礎(chǔ)上添加了時間上分辨率，成為四維核磁共振（4D MRI）血流成像[10]。

基于彩色多普勒的血流向量圖將血流分解為進(jìn)入觀測平面的基本流和平面內(nèi)封閉的渦流，可測量心動周期不同時刻心腔內(nèi)任意一點(diǎn)真實(shí)的血流速度。首先通過彩色多普勒成像，掃描心室內(nèi)的血流區(qū)域，得到心室內(nèi)各個位置沿掃描波束方向的血液流動狀況。再根據(jù)流函數(shù)和流距離函數(shù)得到速度場在垂直于波束方向的速度分量，從而得到心室內(nèi)血流速度場二維分布。

核磁共振的成像方法在三維空間具有很高的空間分辨率，但是該方法在時間分辨率方面受到了極大的限制，計算非常消耗時間。彩色多普勒具有廉價、耗時短、時間分辨率高等優(yōu)點(diǎn)，但是該方法尚未被體外模型驗證，且對垂直于聲束方向的渦旋測量不準(zhǔn)確。

粒子圖像測速（Particle Image Velocimetry，PIV）技術(shù)是用于流體顯像的一種光學(xué)方法，被廣泛應(yīng)用于流體力學(xué)實(shí)驗室，能測量流體的瞬時速度相關(guān)特性。其原理為在流體中充分植入具有良好流體動力學(xué)特征的粒子示蹤劑，使用脈沖激光照射所測流場的切面區(qū)域，通過成像記錄系統(tǒng)獲取兩次或多次曝光的粒子圖像，形成兩幅或多幅PIV實(shí)驗圖像，再利用圖像互相關(guān)方法進(jìn)行分析，從而得出每一小區(qū)域中粒子圖像的平均位移，由此確定流場切面上整個區(qū)域的二維速度。超聲粒子圖像測速技術(shù)（Echo-PIV）采用超聲波束作為顯像源，微泡作為粒子示蹤劑，能夠顯示心腔內(nèi)血流動力。2000年，echo-PIV初次用于實(shí)驗室顯像挾沙水流中的高嶺土粒子，隨后該技術(shù)成功應(yīng)用于心腔內(nèi)血流實(shí)驗和臨床研究[11-13]。

通過超聲粒子圖像測速技術(shù)來計算心臟內(nèi)渦旋特性還是一個較新的研究方向，關(guān)鍵是要證明其具有足夠的可靠性和精確性，以便用于臨床應(yīng)用。為此，本文構(gòu)建了一個超聲粒子圖像測速技術(shù)系統(tǒng)，搭建了一個人體血液循環(huán)仿真系統(tǒng)，通過聚乙烯醇（PVA）制作了左心室仿體，以便在受控條件下測試超聲粒子圖像測速技術(shù)計算心臟內(nèi)渦旋特性的可行性，報道如下。

1 實(shí)驗系統(tǒng)

1.1 超聲粒子圖像測速技術(shù)總體流程

超聲粒子圖像測速技術(shù)總體示意圖，見圖1。本文采用超聲線陣探頭進(jìn)行圖像數(shù)據(jù)的采集，通過超聲探頭對仿體進(jìn)行連續(xù)的B-mode超聲成像，然后將前后兩幀B-mode超聲圖像劃分為有限個小的窗口，對前后兩張圖像對應(yīng)位置的小窗口進(jìn)行互相關(guān)分析，互相關(guān)平面極值所對應(yīng)的位置即為最佳匹配點(diǎn)即小窗口運(yùn)動到的位置，從而可以計算粒子的位移，由于B-mode超聲成像的成像頻率已知，所以可以求得整個流場的速度矢量圖。

圖1 超聲粒子圖像測速技術(shù)總體示意圖

1.2 人體血液循環(huán)仿真系統(tǒng)的搭建

人體血液循環(huán)仿真系統(tǒng)，見圖2。脈動泵采用Harvard Apparatus公司生產(chǎn)的大動物血泵，可以通過大動物血泵調(diào)整脈動頻率、射血分?jǐn)?shù)、射血體積等參數(shù)。超聲采集設(shè)備使用的是Sonix RP系統(tǒng)，通過該系統(tǒng)可以對仿體進(jìn)行連續(xù)成像，通過調(diào)整掃描線密度、圖像寬度、圖像深度等參數(shù)，最快可以實(shí)現(xiàn)每秒鐘上千幀的成像速度。流量計可以實(shí)時提供流入仿體液體的流量信息。試驗中采用的仿體分為兩種，一種是自制的聚乙烯醇（PVA）材料仿體，一種是購買的硅橡膠材料仿體，實(shí)驗的過程中在液體水池中打入超聲造影微泡，開啟脈動泵循環(huán)多個周期，待超聲造影微泡均勻分布在液體中時，對仿體進(jìn)行實(shí)驗采集。在仿體的后方放置吸收超聲的材料，以減少超聲數(shù)據(jù)采集時的噪聲干擾。

圖2 人體血液循環(huán)仿真系統(tǒng)

1.3 用于算法驗證的左心室仿體

左心室硅膠仿體，見圖3。此仿體購買自Shelley Medical Imaging Techologies公司。

圖3 左心室硅膠仿體

該仿體擁有左心室、主動脈弓、主動脈、冠狀動脈，上方為大腦和上肢供給血液的血管。實(shí)驗時，將此仿體連接到血液循環(huán)系統(tǒng)中，打入超聲造影微泡，采集數(shù)據(jù)，應(yīng)用Echo-PIV算法對圖像進(jìn)行計算。由于此仿體是由硅膠制成的，對于超聲波的反射很強(qiáng)，在成像的過程中表現(xiàn)出較多噪聲，因此主要是采用這個仿體進(jìn)行算法的驗證。

1.3.1 聚乙烯醇（PVA）水凝膠的制備

聚乙烯醇PVA水溶液在室溫下可以通過鏈段之間的氫鍵逐漸形成水凝膠，但是這種水凝膠的力學(xué)性能較差，使用價值較低。為了獲得高強(qiáng)度、高含水的PVA水凝膠，有研究者通過不同的交聯(lián)方法和反應(yīng)條件來直接影響高分子的網(wǎng)絡(luò)結(jié)構(gòu)從而改善其性能。

PVA水凝膠的制備按照交聯(lián)方法可分為化學(xué)交聯(lián)和物理交聯(lián)。化學(xué)交聯(lián)又分為輻射交聯(lián)和化學(xué)試劑交聯(lián)。輻射交聯(lián)主要是利用電子束、紫外線等直接輻射PVA溶液，使PVA分子間通過自由基而交聯(lián)在一起。化學(xué)試劑交聯(lián)則是采用化學(xué)交聯(lián)劑使PVA分子間發(fā)生化學(xué)反應(yīng)而交聯(lián)形成凝膠，常用的化學(xué)交聯(lián)劑有醛類、硼酸、環(huán)氧氯丙烷以及可與PVA通過配位絡(luò)合形成凝膠的重金屬鹽等。物理交聯(lián)主要有反復(fù)凍結(jié)法和凍結(jié)部分脫水法。

水凝膠中PVA濃度為10wt%、散射子（Sigmacell Cellulose）濃度為3wt%，制作過程如下：

（1）準(zhǔn)確稱取PVA聚合物加入燒杯中，加入定量的純凈水，在燒杯中加入磁轉(zhuǎn)子，并放入恒溫加熱磁力攪拌器中進(jìn)行加熱。30 ℃下攪拌1～2 h，使PVA聚合物充分溶脹，以利于后面加熱時PVA聚合物的溶解。

（2）調(diào)整溫度到95 ℃，加熱攪拌1 h，觀察溶液中沒有懸浮的顆粒物，基本變?yōu)榫煌该鲿r即可，在加熱的過程中為了防止水分的流失，可以在容器上面覆蓋一層錫箔紙膜。

（3）在PVA聚合物完全溶解后，加入已經(jīng)稱好的Sigmacell Cellulose，繼續(xù)加熱攪拌15 min，使Sigmacell Cellulose均勻分布于PVA水溶液中，這時溶液呈現(xiàn)為乳白色。

（4）停止加熱攪拌后，由于粘度較大，在溶液中可能會有氣泡，把容器放入超聲清潔器，震蕩15 min，除去溶液內(nèi)部的氣泡，然后靜置冷卻至30 ℃。

1.3.2 左心室仿體的制備

左心室仿體模具，見圖4。制作左心室仿體主要以圖4模型為基礎(chǔ)，以PVA水凝膠為材料。制作仿體需要一大一小兩套模具（圖4A），首先使用小模具（圖4A-1）制作出一個實(shí)心模型，然后將此模型放入大模具（圖4A-2）中，在大模具中形成空余位置，并在該位置加入PVA。具體流程如下：

（1）配置5.5%wt的瓊脂糖溶液，將溶液放入微波爐中進(jìn)行加熱直至完全溶解。

（2）將溶液倒入小模具（圖4A-1）中，靜置一段時間待液體完全凝固。打開模型（圖4C-1），即可得到左心室瓊脂糖模型（圖4B），將瓊脂糖模型放入大模具中（圖4C-2），可以看出大模具中瓊脂糖模型周圍有一部分空余的位置。

（3）從大模具上方小孔處注射PVA水凝膠直至溢出，然后將大模具密封。將密封好的大模具放入-20 ℃的冰箱中凍融12 h。最后將凍融好的仿體取出，通過擠壓碾碎心室仿體內(nèi)的瓊脂糖，用水沖出瓊脂糖，得到PVA心室模型（圖4D）。

圖4 左心室仿體模具

2 左心室渦旋參數(shù)的計算

在心室舒張充盈的過程中，血流會在心室腔內(nèi)形成渦旋，渦旋可以最小化血液流動過程中能量的損失，并減少心臟射血時所需要的能量，所以渦旋是心臟內(nèi)血液流動的重要參數(shù)。通過Echo-PIV算法可以得到心臟內(nèi)血液流動的速度分布，進(jìn)一步可計算出心臟內(nèi)渦旋的相關(guān)參數(shù)。

渦量是描述旋渦運(yùn)動最重要的物理量之一，定義為流體速度矢量v的旋度，渦量的單位是秒分之一（s-1）。渦量的計算公式：

渦量是1個周期性的變量，通過對多幅圖像相同位置渦量值組成的序列進(jìn)行傅里葉變化，得到圖像各個位置的零階諧波值和一階諧波值，即公式(2)中的分量ω0(x，y)，ω1(x，y)：

渦旋的相對強(qiáng)度是評判心臟功能的重要參數(shù)，當(dāng)發(fā)生心力衰竭等疾病時，心臟的收縮舒張功能會減弱，從而使得心臟左心室內(nèi)渦旋的相對強(qiáng)度降低。渦旋的相對強(qiáng)度（Relative Strength，RS）是一階諧波強(qiáng)度值與零階諧波值的比，可通過公式(3)計算得到：

其中：

渦旋的波動是平方歸一化的渦量變化：

公式(5)中，T為一次心跳持續(xù)時間，LV代表左心室。液體流動時相互之間會產(chǎn)生摩擦，從而導(dǎo)致能量損失。能量損失也是評價心臟功能的一個重要參數(shù)，當(dāng)心臟發(fā)生病變時，左心室內(nèi)的血流不能繼續(xù)保持良好的渦旋結(jié)構(gòu)，從而會加大能量的損失，其計算公式如下：

其中，ρ代表液體的密度，μ代表液體的粘度，DI代表能量的損失。

能量的波動是相對于平均能量的能量變化，其計算公式如下：

上式中下角標(biāo)“0”代表平均速度。

3 實(shí)驗

3.1 PVA仿體實(shí)驗

當(dāng)心臟發(fā)生心肌梗死、肥厚性心肌癥時，心肌的彈性會發(fā)生變化，同時心臟左心室內(nèi)的血流動力學(xué)參數(shù)也會發(fā)生變化。本文采用不同硬度的仿體測試對血流動力學(xué)參數(shù)的影響。在制作PVA仿體時，會經(jīng)歷凍融階段，不同凍融周期制作出的仿體會出現(xiàn)不同的硬度。通過實(shí)驗，可測量1～8個凍融周期的仿體彈性模量，本研究使用了3周期、7周期的仿體進(jìn)行了實(shí)驗。

3.2 實(shí)驗數(shù)據(jù)的處理

對采集的超聲圖像進(jìn)行處理，步驟如下：

（1）輸入連續(xù)的含有超聲造影微泡的心臟左心室長軸切面圖像，總幀數(shù)N應(yīng)涵蓋至少一個心動周期內(nèi)超聲系統(tǒng)所采集的圖像。

（2）在第一幀圖像上選擇一個感興趣區(qū)域（Region of Interest，ROI）。

（3）對第n幀和第n+1幀圖像的ROI進(jìn)行匹配計算。將ROI劃分為多個分析窗口（次窗口），兩幅圖中對應(yīng)的兩個次窗口進(jìn)行二維互相關(guān)運(yùn)算，獲得該次窗口所代表的血流的位移矢量，然后利用三點(diǎn)高斯峰擬合算法進(jìn)行亞像素分析，利用全局和局部中值濾波器去除錯誤矢量，并采用雙線性內(nèi)插算法對錯誤矢量進(jìn)行替換。通過迭代提高計算的精度，最后減小窗口大小，提高空間分辨率。依次對ROI每一對次窗口進(jìn)行該運(yùn)算，得到第n幀圖像ROI內(nèi)的二維位移矢量分布圖。

（4）n=n+1，判斷n是否大于或等于N。若“否”，則返回第（3）步。若“是”，則進(jìn)行第（5）運(yùn)算。通過對第（3）步的循環(huán)運(yùn)算，得到輸入圖像ROI的二維位移矢量分布圖。換言之，得到了若干個心動周期內(nèi)每一個心跳時刻心臟內(nèi)血液流動的位移矢量圖。

（5）通過步驟（4）得到的位移矢量圖和超聲獲取圖像的時間間隔，計算得到心臟內(nèi)血液速度分布圖。

（6）通過左心室內(nèi)流體的速度，計算出心臟左心室內(nèi)血流動力學(xué)相關(guān)參數(shù)。

4 實(shí)驗結(jié)果與討論

4.1 硅膠仿體實(shí)驗結(jié)果與討論

在硅膠仿體實(shí)驗中，實(shí)驗參數(shù)如下：脈動泵頻率35 Hz、脈動泵每搏射血體積35 mL、超聲探頭12 MHz（線陣探頭），線密度256。硅膠仿體實(shí)驗估計渦旋參數(shù)，見表1。

表1 硅膠仿體實(shí)驗估計渦旋參數(shù)

在人體的整個心動周期內(nèi)，心臟血流并非都處于渦旋狀態(tài)，但在心臟填充期，左心室內(nèi)血流處于渦旋狀態(tài)。對于整個實(shí)驗系統(tǒng)來說，脈動泵相當(dāng)于左心房，可為心室供血，本研究選擇脈動泵供血時期獲得的B超圖像進(jìn)行研究（圖5）。

圖5 硅膠仿體實(shí)驗中供血時期的B超圖像

本研究同時計算了速度（Velocity）、剪切力（Wall Shear Stress，WSS）、渦量（Vorticity）等信息。速度、剪切力、渦量的彩色編碼矢量圖，見圖6～8。從圖中可知，本研究采用的算法可精確地計算出心臟左心室仿體中的血液流動。

4.2 PVA實(shí)驗仿體結(jié)果與討論

本研究采用3周期彈性模量為168.52 kPa和7周期彈性模量為296.04 kPa的仿體進(jìn)行實(shí)驗。每個周期的仿體分別進(jìn)行3組實(shí)驗。實(shí)驗參數(shù)如下：脈動泵頻率15 Hz、脈動泵每搏射血體積15 mL、超聲探頭12 MHz（線陣探頭），線密度128。PVA仿體實(shí)驗結(jié)果，見表2。

由實(shí)驗結(jié)果可知，兩組間DI、EF、RS、W等參數(shù)無統(tǒng)計學(xué)差異。彈性模量大的仿體其形狀比彈性模量小的仿體小，且差異有統(tǒng)計學(xué)意義（P＜0.05）。

本研究對仿體內(nèi)部的血流速度、剪切力和渦量進(jìn)行了比較，并對這3個參數(shù)同樣的尺度進(jìn)行了彩色編碼。彩色編碼對比圖像，見圖9。

本研究中，兩組之間彩色編碼的度量相同，圖像的右端貼有不同顏色彩色編碼所對應(yīng)的值，紅色代表數(shù)值大，藍(lán)色代表數(shù)值小。通過對圖像的觀察，我們發(fā)現(xiàn)7周期的速度、剪切應(yīng)力、渦量均比3周期大。PVA仿體實(shí)驗結(jié)果，見表3。

表2 PVA仿體實(shí)驗結(jié)果

表3 PVA仿體實(shí)驗結(jié)果

圖6 速度彩色編碼矢量圖

圖7 剪切力彩色編碼矢量圖

圖8 渦量彩色編碼矢量圖

圖9 彩色編碼對比圖像

5 結(jié)束語

左心室血液的流動動態(tài)可被認(rèn)為是心臟健康一個潛在的指示器，通過超聲粒子圖像測速技術(shù)來計算心臟內(nèi)渦旋特性還是一個較新的研究方向。本文構(gòu)建了一個超聲粒子圖像測速技術(shù)系統(tǒng)，搭建了一個人體血液循環(huán)仿真系統(tǒng)，通過聚乙烯醇（PVA）制作了左心室仿體，完成了一些初步的實(shí)驗，結(jié)果證明所構(gòu)建的系統(tǒng)、仿體、相關(guān)算法是可行的，為下一步工作打好了基礎(chǔ)。

本研究由中科院深圳先進(jìn)技術(shù)研究院生物醫(yī)學(xué)與健康工程研究所完成，仍有以下需要改進(jìn)的地方：

（1）實(shí)驗期間，由于缺少采集大鼠心臟信號的探頭，未能完成動物實(shí)驗。若能購買到實(shí)驗設(shè)備進(jìn)行動物實(shí)驗，可使實(shí)驗結(jié)果更有說服力。

（2）可與醫(yī)院進(jìn)行合作，采集病人信息，有利于進(jìn)一步探索超聲粒子圖像測速技術(shù)對心臟功能的評估。

（3）改善實(shí)驗算法，使算法能夠更加適用于心室這種非規(guī)則形狀內(nèi)流場的計算，能夠更加準(zhǔn)確的計算出心室內(nèi)的流場。

（4）改善體外實(shí)驗系統(tǒng)，使之能夠更加接近于人體真實(shí)的情況，能夠使采集到的數(shù)據(jù)更接近于真實(shí)。

[1] Sengupta PP,Pedrizzetti G,Kilner PJ,et al.Emerging trendsin CV flow visualization[J].JACC Cardiovasc Imaging,2012,5(3):305-316.

[2] Kanski M,Arvidsson PM,T?ger J,et al.Left ventricular fluid kinetic energy time curves in heart failure from cardiovascular magnetic resonance 4D flow data[J].J Cardiovasc Magn Reson,2015,17:111.

[3] Khalafvand SS,Hung TK,Ng EY,et al.Kinematic, Dynamic, and Energy Characteristics of Diastolic Flow in the Left Ventricle[J].Comput Math Methods Med,2015:701945.

[4] Kheradvar A,Gharib M.Influence of ventricular pressure drop on mitral annulus dynamics through the process of vortex ring formation[J].Ann Biomed Eng,2007,35(12):2050-2064.

[5] Jang J,Ahn CY,Jeon K,et al.A Reconstruction Method of Blood Flow Velocity in Left Ventricle Using Color Flow Ultrasound[J].Comput Math Methods Med,2015:108274.

[6] Choi YJ,Constantino J,Vedula V,et al.A New MRI-Based Model of Heart Function with Coupled Hemodynamics and Application to Normal and Diseased Canine Left Ventricles[J].Front Bioeng Biotechnol,2015,3:140

[7] Elbaz MS,Calkoen EE,Westenberg JJ,et al.Vortex flow during early and late left ventricular filling in normal subjects: quantitative characterization using retrospectively-gated 4D flow cardiovascular magnetic resonance and three-dimensional vortex core analysis[J].J Cardiovasc Magn Reson,2014,16:78.

[8] Subramanian A,Mu H,Kadambi JR,et al.Particle image velocimetry investigation of intravalvular flow fields of a bileaflet mechanical heart valve in a pulsatile flow[J].J Heart Valve Dis,2000,9(5):721-731.

[9] Taylor TW,Yamaguchi T.Flow patterns in three-dimensional left ventricular systolic and diastolic flows determined from computational fluid dynamics[J].Biorheology,1995,32(1):61-71.

[10] Nilsson A,Bloch KM,T?ger J,et al.Accuracy of four-dimensional phase-contrast velocity mapping for blood flow visualizations: a phantom study[J].Acta Radiol,2013,54(6):663-671.

[11] Hong GR,Pedrizzetti G,Tonti G,et al.Characterization and quantification of vortex flow in the human left ventricle by contrast echocardiography using vector particle image velocimetry[J].JACC Cardiovasc Imaging,2008,1(6):705-717.

[12] Agati L,Cimino S,Tonti G,et al.Quantitative analysis of intraventricular blood flow dynamics by echocardiographic particle image velocimetry in patients with acute myocardial infarction at different stages of left ventricular dysfunction[J].Eur Heart J Cardiovasc Imaging,2014,15(11):1203-1212.

[13] 張波,許婧,智光.超聲粒子圖像測速技術(shù)檢測心腔內(nèi)渦流的研究進(jìn)展[J].中國醫(yī)藥導(dǎo)報,2015,(11):30-33.

Preliminary Study of Application of the Ultrasonic Particle Image Velocimetry Technique in Evaluation of Left Ventricular Vortex Characteristics

ZHANG Yao-nan1,2, ZHOU Wei2,3, NIU Li-li3
1. College of Electronics and Information Engineering, Xi’an Siyuan University, Xi’an Shaanxi 710038, China; 2. Sino-Dutch Biomedical and Information Engineering School, Northeastern University, Shenyang Liaoning 110169, China; 3. Institute of Biomedical and Health Engineering, Shenzhen Institutes of Advanced Technology, Chinese Academy of Science, Shenzhen Guangdong 518055, China

The vortex properties of left ventricles can be considered as a potential indicator for heart’s health status. Using the ultrasonic particle image velocimetry (PIV) technique to calculate the characteristics of vortex in the heart is a relatively new research direction, and the key is to prove its adequate reliability and accuracy for clinical applications. For this purpose, a ultrasonic particle image velocimetry system was constructed to build a human blood circulation simulation system. And the polyvinyl alcohol (PVA) was used to make the left ventricular phantoms, so that the ultrasonic PIV technique could be tested in controlled conditions for the calculation of the vortex characteristics in the heart. The system achieved continuous imaging of the body, by adjusting the scanning line density, image width, image depth and other parameters, with the fastest speed of thousands of frames per second. Cardiac blood flow velocity profiles were obtained through the algorithm of acoustic PIV and further calculation of relevant parameters for the vortex in the heart, including energy loss (DI) and energy fluctuations (EF), vortex relative strength (RS) and eddy fluctuations (W). The results show that the constructed system, the simulated body and the related algorithms are feasible, which is the basis for the next step of work.

ultrasound; particle image velocimetry; left ventricles; vortex

R445.1

A

10.3969/j.issn.1674-1633.2016.06.007

1674-1633(2016)06-0036-06

2016-01-27

遼寧省自然科學(xué)基金資助項目（201202071），西安思源學(xué)院校級重大科研項目（XASY-B1601）。

本文作者：張耀楠，教授，主要研究方向為醫(yī)學(xué)影像分析。

作者郵箱：johanzyn@qq.com