高銘澤 陳廣新 陳曉琴 胡明成 成志國(guó) 林 強(qiáng)▲

1.牡丹江醫(yī)學(xué)院醫(yī)學(xué)影像學(xué)院，黑龍江牡丹江 157011；2.牡丹江醫(yī)學(xué)院附屬紅旗醫(yī)院，黑龍江牡丹江 157011

頸動(dòng)脈為顱內(nèi)組織供血，頸動(dòng)脈斑塊是引起腦卒中的重要原因之一[1]。頸動(dòng)脈狹窄的形成原因比較復(fù)雜，血管內(nèi)膜損傷、血管壁蛋白纖維作用減弱導(dǎo)致的血管壁硬化、變厚等均促成頸動(dòng)脈粥樣硬化[2]。已有研究證實(shí)，血流動(dòng)力學(xué)因素是頸動(dòng)脈粥樣硬化和血栓形成的重要因素[3]。目前，現(xiàn)有的醫(yī)學(xué)檢查手段仍無(wú)法對(duì)體內(nèi)的血流動(dòng)力學(xué)分布特征進(jìn)行有效測(cè)量。近些年計(jì)算流體力學(xué)（computational fluid dynamics，CFD）快速發(fā)展，將CFD技術(shù)應(yīng)用于血流動(dòng)力學(xué)仿真分析可實(shí)現(xiàn)頸動(dòng)脈內(nèi)的血流動(dòng)力學(xué)分布特征參數(shù)，對(duì)頸動(dòng)脈狹窄的發(fā)生發(fā)展、病理機(jī)制研究有重要的意義[4]。本研究基于MIMICS軟件對(duì)個(gè)體化CTA真實(shí)患者數(shù)據(jù)構(gòu)建輕中度三種狹窄程度的幾何模型、利用CFD仿真軟件SimVascular進(jìn)行血流動(dòng)力學(xué)仿真計(jì)算，分析不同狹窄程度對(duì)頸動(dòng)脈的血流動(dòng)力學(xué)指標(biāo)的特征分布規(guī)律，并探討?yīng)M窄率大小與頸動(dòng)脈斑塊的形成關(guān)系。

1 模型構(gòu)建

1.1 數(shù)據(jù)采集與幾何模型構(gòu)建

采集牡丹江醫(yī)學(xué)院附屬紅旗醫(yī)院影像科1例頸動(dòng)脈狹窄患者的CTA影像數(shù)據(jù)作為研究對(duì)象，經(jīng)牡丹江醫(yī)學(xué)院附屬紅旗醫(yī)院醫(yī)學(xué)倫理委員會(huì)批準(zhǔn)。影像采集設(shè)備為GE 64排螺旋CT，電壓120 kV，電流300 mA，層距0.5 mm，層厚0.3 mm，螺距1.75 mm，矩陣512×512，共獲得120張圖像。使用MIMICS軟件（Materialise公司，比利時(shí)）對(duì)CTA文件進(jìn)行閾值分割、蒙版編輯、區(qū)域生長(zhǎng)、三維重建等操作生成幾何模型（STL格式），利用3-matic（MIMICS軟件附帶）對(duì)頸動(dòng)脈初步模型進(jìn)行優(yōu)化、光滑生成頸動(dòng)脈狹窄模型，并在此基礎(chǔ)上構(gòu)建正常頸動(dòng)脈、輕度狹窄、重度狹窄幾何模型，見(jiàn)圖1。

1.2 CFD仿真計(jì)算模型構(gòu)建

采用SimVascular（開(kāi)源血流動(dòng)力學(xué)仿真軟件）劃分網(wǎng)格，網(wǎng)格邊界層設(shè)為5層，增長(zhǎng)率為1.2，入口采用速度入口，入口速度曲線見(jiàn)圖2。出口采用壓力出口，血液生理壓力為13 339 Pa[5]。計(jì)算模型采用非牛頓血液黏度模型Carreau-Yasuda，血液密度為1060 kg/m3[6]。在Carreau-Yasuda模型中，剪切率與動(dòng)力黏度關(guān)系見(jiàn)圖3。血液為各向同性、不可壓縮、層流，其控制方程為Navier-Stokes方程。計(jì)算時(shí)間步間隔0.01 s，周期為0.8 s，總計(jì)算時(shí)間為2.4 s即3個(gè)周期。迭代精度設(shè)為10-4。

圖2 入口速度曲線

圖3 非牛頓血液黏度模型

2 結(jié)果

2.1 血流速度分布特征分析

由頸動(dòng)脈內(nèi)血流速度（峰值時(shí)刻）流線圖（圖4）可知，頸動(dòng)脈狹窄處附近形成了血流速度的最大值區(qū)，血液在流經(jīng)狹窄區(qū)域之前，頸動(dòng)脈血流在頸總動(dòng)脈分布均勻，血液到達(dá)分叉部位狹窄區(qū)域血流發(fā)生了較大的流動(dòng)分離，血液的穩(wěn)定層被破壞，呈現(xiàn)復(fù)雜的流態(tài)分布。隨著狹窄率增高，血管壁對(duì)血液的流動(dòng)速度產(chǎn)生了影響，在狹窄率較高的區(qū)域，血流速度更快，在狹窄部位出現(xiàn)局部的渦流、回流、二次流現(xiàn)象，頸動(dòng)脈狹窄部位的血流速度、方向的震蕩對(duì)動(dòng)脈粥樣硬化斑塊的穩(wěn)定性、形成、發(fā)展以及血管壁的重塑、修復(fù)具有重要的影響。

圖4 血液流線圖

2.2 壁面切應(yīng)力分布特征分析

壁面切應(yīng)力分布圖（圖5）結(jié)果表明，在頸動(dòng)脈狹窄處呈現(xiàn)壁面切應(yīng)力較高的情況，尤其是狹窄率越高，壁面切應(yīng)力的值越高、高壁面切應(yīng)力的區(qū)域面積越大。隨著頸動(dòng)脈狹窄率增高，狹窄區(qū)域過(guò)高的壁面切應(yīng)力易導(dǎo)致血管斷裂，也易導(dǎo)致斑塊的脫落，這也構(gòu)成了缺血性腦卒中發(fā)生的可能。高壁面切應(yīng)力使紅細(xì)胞聚集，增加斑塊破裂、脫落堵塞血管的風(fēng)險(xiǎn)。

圖5 壁面切應(yīng)力分布特征

2.3 壁面壓力分析

頸動(dòng)脈心動(dòng)周期壓力分布曲線（圖6）結(jié)果表明，在心動(dòng)周期內(nèi)動(dòng)脈狹窄部位存在負(fù)壓現(xiàn)象。對(duì)比正常、輕度狹窄、重度狹窄的頸動(dòng)脈狹窄部位的壁面壓力，隨著狹窄率越大，狹窄部位壁面壓力的波動(dòng)范圍較大，壓力值逐漸減小。三種頸動(dòng)脈壁面壓力值差異較大，頸動(dòng)脈狹窄對(duì)遠(yuǎn)端血管的血液流動(dòng)產(chǎn)生影響，狹窄部位的低壓力導(dǎo)致遠(yuǎn)端的血流動(dòng)力不足，以至于出現(xiàn)血液回流的現(xiàn)象，甚至導(dǎo)致腦部供血不足，進(jìn)而導(dǎo)致腦卒中發(fā)生。

圖6 頸動(dòng)脈狹窄部位壓力曲線圖

3 討論

頸動(dòng)脈是為大腦供血的重要通道，頸動(dòng)脈斑塊導(dǎo)致頸動(dòng)脈狹窄，嚴(yán)重威脅著大腦的血流運(yùn)輸，腦供血不足、腦卒中的發(fā)生也與頸動(dòng)脈狹窄有著密切聯(lián)系[7]。動(dòng)脈粥樣硬化的形成發(fā)展最重要的因素是血流動(dòng)力學(xué)因素。因此，不同狹窄率頸動(dòng)脈的血流動(dòng)力學(xué)參數(shù)如壁切應(yīng)力、壁面壓力、剪切應(yīng)力等特征分布情況，可為頸動(dòng)脈狹窄的形成機(jī)制提供理論指導(dǎo)，也可為頸動(dòng)脈狹窄的疾病診療提供理論指導(dǎo)[7-11]。本研究利用個(gè)體化真實(shí)病例的影像CTA數(shù)據(jù)，使用MIMICS軟件進(jìn)行三維重建獲得個(gè)體化頸動(dòng)脈的初步三維模型，在此基礎(chǔ)上利用3-matic軟件構(gòu)建不同狹窄率的頸動(dòng)脈幾何模型。這種構(gòu)建模型的技術(shù)可以方便、快捷、有針對(duì)性地創(chuàng)建目標(biāo)模型。

頸動(dòng)脈狹窄是由于動(dòng)脈粥樣硬化引起的血管局部結(jié)構(gòu)異常、腦卒中的發(fā)生。頸動(dòng)脈分叉部位復(fù)雜的血流動(dòng)力學(xué)特征分布如壁切應(yīng)力、壁壓力、血流速度等。頸動(dòng)脈血液在頸動(dòng)脈狹窄區(qū)域流動(dòng)時(shí)的異?，F(xiàn)象如血流回流、低速流動(dòng)、渦流等刺激血細(xì)胞的滯留、沉積，并刺激平滑肌細(xì)胞和內(nèi)皮細(xì)胞，導(dǎo)致增值平滑肌細(xì)胞，上述因素也是導(dǎo)致動(dòng)脈粥樣硬化的基礎(chǔ)[12]。頸動(dòng)脈血液的流動(dòng)跟血管壁的壓力密切相關(guān)，頸動(dòng)脈的血流速度和供血量不足會(huì)促進(jìn)該紅細(xì)胞的沉積和滯留，而壓力梯度較小會(huì)導(dǎo)致頸動(dòng)脈血流速度變慢、供血量不足[13]。一些研究表明[7-15]，動(dòng)脈硬化與血栓的形成與壁切應(yīng)力密切相關(guān)，低壁切應(yīng)力對(duì)近血管壁的邊界層區(qū)域血液正常的層流狀態(tài)產(chǎn)生破壞作用，低壁切應(yīng)力導(dǎo)致血流出現(xiàn)渦流、回流，血細(xì)胞在近壁面流動(dòng)時(shí)，不斷與血管內(nèi)壁層內(nèi)皮細(xì)胞和平滑肌細(xì)胞摩擦，促進(jìn)血小板活性的增加，使血小板在損傷部位黏附聚集，從而導(dǎo)致血管內(nèi)膜增生，血管發(fā)生重構(gòu)，進(jìn)而產(chǎn)生動(dòng)脈粥樣硬化并引起血管狹窄。臨床上對(duì)頸動(dòng)脈狹窄的治療可采用支架介入等方式，支架的大小、形狀、體積等設(shè)計(jì)對(duì)狹窄頸動(dòng)脈的血液流場(chǎng)及應(yīng)力變化起到調(diào)節(jié)作用。

本研究對(duì)個(gè)體化不同狹窄率的頸動(dòng)脈非牛頓、層流的血流動(dòng)力學(xué)仿真計(jì)算研究，比較血流速度、壓力、壁面切應(yīng)力的分布特征，研究結(jié)果表明：①不同狹窄率的頸動(dòng)脈在狹窄部位血流分布存在較大差異，即血流速度與頸動(dòng)脈狹窄率呈正相關(guān)，狹窄中心血流速度較大，靠近壁面區(qū)域的血流速度較小，狹窄部位出現(xiàn)渦流并且動(dòng)脈越狹窄，渦流越明顯加強(qiáng)；②壁面切應(yīng)力隨著動(dòng)脈狹窄率增大，在狹窄部位出現(xiàn)峰值，并且頸動(dòng)脈無(wú)狹窄區(qū)域也受到影響。壁面切應(yīng)力過(guò)大對(duì)血管硬化和斷裂有積極作用，頸動(dòng)脈斑塊脫落的危險(xiǎn)因素增加，缺血性卒中風(fēng)險(xiǎn)增大；③在一個(gè)心動(dòng)周期內(nèi)，頸動(dòng)脈狹窄的壓力出現(xiàn)負(fù)壓現(xiàn)象，即產(chǎn)生負(fù)壓效應(yīng)。負(fù)壓效應(yīng)導(dǎo)致在頸動(dòng)脈狹窄部位血細(xì)胞沉積，進(jìn)而形成頸動(dòng)脈粥樣硬化斑塊，進(jìn)一步造成動(dòng)脈狹窄。

綜上所述，本研究采用考慮血液的非牛頓效應(yīng)的單相流數(shù)值模擬分析，這種建模方法雖然能夠比較準(zhǔn)確地考察頸動(dòng)脈的血流動(dòng)力學(xué)指標(biāo)，但是忽略了頸動(dòng)脈彈性血管的影響。因此，本研究存在一定的不足，在進(jìn)一步的研究中，筆者將考慮采用非牛頓流體的流固耦合分析，以提高模擬的準(zhǔn)確性。