大型未破裂顱內(nèi)動(dòng)脈瘤的血流動(dòng)力學(xué)特點(diǎn)分析

李 良，張 瑩，楊新健，陳家亮，王盛章，李海云，于紅玉，汪 陽(yáng)，李傳輝

(1北京大學(xué)第一醫(yī)院，北京 100050;2北京市神經(jīng)外科研究所;3復(fù)旦大學(xué)力學(xué)與工程科學(xué)系;4首都醫(yī)科大學(xué)生物醫(yī)學(xué)工程學(xué)院)

研究發(fā)現(xiàn)，動(dòng)脈瘤的直徑是動(dòng)脈瘤破裂最重要的獨(dú)立危險(xiǎn)因子，即大動(dòng)脈瘤存在更高的自然破裂風(fēng)險(xiǎn)［1～5］。血流動(dòng)力學(xué)作用力被認(rèn)為在動(dòng)脈瘤的形成、發(fā)展和破裂的自然病程中具有極其重要的作用［6，7］。2010 年7 月 ～2011 年6月，我們對(duì)12 例直徑≥15 mm 的未破裂顱內(nèi)動(dòng)脈瘤進(jìn)行了血流動(dòng)力學(xué)數(shù)值模擬，對(duì)其流速、流動(dòng)模式、壁面切應(yīng)力(WSS)、切應(yīng)力震蕩因子(OSI)等進(jìn)行了分析，并將動(dòng)脈瘤和其載瘤動(dòng)脈的參數(shù)比較。現(xiàn)報(bào)告如下。

1 資料與方法

1.1 臨床資料 12例未發(fā)生破裂出血的直徑≥15 mm顱內(nèi)囊性動(dòng)脈瘤患者，其中男3例，女9例;年齡39～67歲。首發(fā)臨床表現(xiàn)為頭痛伴眼瞼下垂和(或)復(fù)視、視物模糊4例，單側(cè)肢體無(wú)力伴言語(yǔ)不清4例，頭痛頭暈伴有耳鳴2例，頭痛伴單眼失明及單純頭痛頭暈各1例。動(dòng)脈瘤部位:頸內(nèi)動(dòng)脈后交通段4例，前交通段3例，頸內(nèi)動(dòng)脈眼動(dòng)脈段2例，頸內(nèi)動(dòng)脈C2段、海綿竇段及大腦中動(dòng)脈分叉各1例。形狀分型:側(cè)方型和頂端(分叉處)型各6例。

1.2 血流動(dòng)力學(xué)數(shù)值模擬方法 本研究所采用的模擬方法詳見(jiàn)我們的前期研究［8，9］。

1.2.1 數(shù)值模擬的硬件及軟件 本研究所應(yīng)用的計(jì)算機(jī)為特別配置的惠普工作站，主頻為Intel Xeon 3.0 G，內(nèi)存為8.0 G，硬盤為160 G。主要軟件有:本課題組自編程序及Matlab 7.0、Ansys ICEM，Ansys CFX，Photoshop7.0 等。

1.2.2 數(shù)值模擬和流體分析 對(duì)12例患者調(diào)整重建后的三維DSA圖像，去掉不重要的分支，重點(diǎn)分析動(dòng)脈瘤及其載瘤動(dòng)脈。采用本課題組自編寫的程序，及Matlab7.0計(jì)算，合成STL格式的三維圖像。采用ANSYS ICEM軟件進(jìn)行三維體網(wǎng)格劃分。為兼顧計(jì)算的精度和計(jì)效率，本研究劃分的網(wǎng)格數(shù)量為1500000～1800000個(gè)。采用ANSYS CFX 11.0軟件行流體分析計(jì)算。血管入口條件設(shè)定為通過(guò)經(jīng)顱多普勒超聲獲得的搏動(dòng)性血流，出口條件設(shè)定為自由流動(dòng)出口。每個(gè)心動(dòng)周期被劃分為800步，每步時(shí)間為0.001 s。計(jì)算兩個(gè)心動(dòng)周期后取第二個(gè)心動(dòng)周期計(jì)算所得的數(shù)值進(jìn)行分析。在每個(gè)心動(dòng)周期內(nèi)部，動(dòng)脈瘤的血流動(dòng)力學(xué)參數(shù)的分布一致，只是隨著血壓的波動(dòng)，發(fā)生了數(shù)值上的變化，最大的時(shí)刻出現(xiàn)在血壓最高峰的收縮期末。

1.2.3 圖像顯示和數(shù)據(jù)分析 采用ANSYS CFX 11.0軟件結(jié)合自編程序行圖像顯示和數(shù)據(jù)采集。取動(dòng)脈瘤上的數(shù)據(jù)和其載瘤動(dòng)脈(取動(dòng)脈瘤頸部為中心1 cm長(zhǎng)度的載瘤動(dòng)脈)上的數(shù)據(jù)相對(duì)比。用SPSS11.5統(tǒng)計(jì)學(xué)專用軟件包進(jìn)行統(tǒng)計(jì)學(xué)分析。所得計(jì)量數(shù)據(jù)以表示，比較采用配對(duì)t檢驗(yàn)。P≤0.05為差異有統(tǒng)計(jì)學(xué)意義。

2 結(jié)果

2.1 血流流動(dòng)模式及流速 側(cè)方動(dòng)脈瘤血流流動(dòng)模式為入瘤血流從瘤頸口下游側(cè)流向動(dòng)脈瘤內(nèi)，并沿著瘤周壁回流至載瘤動(dòng)脈中，再?gòu)牧鲱i口上游流出，形成一個(gè)主要的渦流。頂端(分叉處)動(dòng)脈瘤血流流動(dòng)模式為血流從入流動(dòng)脈直接沖擊入動(dòng)脈瘤，后沿瘤周壁回流至兩出流動(dòng)脈，形成兩側(cè)兩個(gè)主要渦流。所有的大型動(dòng)脈瘤皆發(fā)現(xiàn)瘤內(nèi)部渦流形態(tài)復(fù)雜，在主要渦流內(nèi)部有大面積紊亂的低速流動(dòng)。動(dòng)脈瘤瘤頸附近速度最大，瘤頂附近流速最低，頂部存在一個(gè)面積較大的低速區(qū)，最小速度接近于0。

2.2 顱內(nèi)動(dòng)脈瘤與其載瘤動(dòng)脈內(nèi)的血流動(dòng)力學(xué)參數(shù) 見(jiàn)表1。瘤壁平均低WSS區(qū)域占動(dòng)脈瘤總表面積的百分?jǐn)?shù)平均為24.74%。

表1 大型動(dòng)脈瘤與其載瘤動(dòng)脈血流動(dòng)力學(xué)參數(shù)比較()

表1 大型動(dòng)脈瘤與其載瘤動(dòng)脈血流動(dòng)力學(xué)參數(shù)比較()

注:與載瘤動(dòng)脈比較，*P ＜0.01

動(dòng)脈 WSS(Pa)OSI動(dòng)脈瘤 4.990 ±1.565* 0.0690 ±0.0274*載瘤動(dòng)脈13.850 ±4.092 0.0040 ±0.0022

3 討論

臨床上大約75%的巨大動(dòng)脈瘤發(fā)現(xiàn)時(shí)均未破裂，因此臨床存在大動(dòng)脈瘤不容易破裂的誤解。本研究12例首發(fā)癥狀均為占位壓迫和缺血表現(xiàn)等動(dòng)脈瘤破裂前的先兆癥狀，亦印證了上述臨床事實(shí)。大量研究結(jié)果表明，動(dòng)脈瘤直徑是動(dòng)脈瘤破裂的最重要的獨(dú)立危險(xiǎn)因子。2003年國(guó)際未破裂動(dòng)脈瘤研究協(xié)作組發(fā)表大樣本多中心研究結(jié)果，認(rèn)為動(dòng)脈瘤瘤體越大，發(fā)生破裂引起SAH機(jī)會(huì)越大［1］。

為了驗(yàn)證動(dòng)脈瘤越大越容易破裂的內(nèi)在原因，本研究應(yīng)用計(jì)算機(jī)流體動(dòng)力學(xué)方法對(duì)大型未破裂動(dòng)脈瘤的內(nèi)在血流動(dòng)力學(xué)因素進(jìn)行了研究，并與其各自的載瘤動(dòng)脈的血流動(dòng)力學(xué)參數(shù)作對(duì)比。研究結(jié)果發(fā)現(xiàn)，所有的大型動(dòng)脈瘤皆發(fā)現(xiàn)瘤內(nèi)部渦流形態(tài)復(fù)雜，在主要渦流內(nèi)部有大面積紊亂的低速流動(dòng)。大型動(dòng)脈瘤壁與其載瘤動(dòng)脈相比具有明顯的低WSS和高 OSI。

Cebral等［6］對(duì)62例破裂和未破裂動(dòng)脈瘤內(nèi)血流流動(dòng)模式的研究證明，未破裂動(dòng)脈瘤大部分具有簡(jiǎn)單穩(wěn)定的血流模式，破裂動(dòng)脈瘤大都具有復(fù)雜紊亂的血流模式。提示復(fù)雜紊亂的血流模式與動(dòng)脈瘤破裂高風(fēng)險(xiǎn)相關(guān)。本研究中大型動(dòng)脈瘤皆發(fā)現(xiàn)瘤內(nèi)部渦流紊亂，提示破裂風(fēng)險(xiǎn)，且瘤頂部存在大面積低速血流，流速甚至接近0。局部血流停滯會(huì)導(dǎo)致紅細(xì)胞集聚、血小板和白細(xì)胞在瘤壁內(nèi)側(cè)的集合和黏附，導(dǎo)致附壁小血栓的形成以及血小板和血漿纖維蛋白單體進(jìn)一步浸潤(rùn)動(dòng)脈瘤壁，增加動(dòng)脈瘤破裂風(fēng)險(xiǎn)［10］。

目前的觀點(diǎn)認(rèn)為動(dòng)脈瘤的破裂與WSS不足密切相關(guān)［7，11～13］，維持正常動(dòng)脈瘤壁形態(tài)需要的 WSS約為2 Pa，WSS低于1.5 Pa時(shí)會(huì)導(dǎo)致內(nèi)皮細(xì)胞的變性與凋亡，促進(jìn)動(dòng)脈瘤的生長(zhǎng)以及最后的破裂［14］。WSS不足直接導(dǎo)致血管內(nèi)皮細(xì)胞出現(xiàn)形變，內(nèi)皮細(xì)胞間隙加大，屏障功能喪失，血液內(nèi)的有害物質(zhì)侵入，使瘤壁彈性纖維及膠原纖維斷裂、缺失，從而形成了動(dòng)脈瘤壁的特有病理學(xué)形態(tài)。隨著動(dòng)脈瘤生長(zhǎng)，WSS進(jìn)一步降低，形成惡性循環(huán)直至破裂。本研究發(fā)現(xiàn)大型動(dòng)脈瘤壁上的WSS與載瘤動(dòng)脈相比顯著偏低，且低WSS區(qū)域平均占24.74%，此提示大型動(dòng)脈瘤壁的內(nèi)皮細(xì)胞更容易變性、凋亡而引起破裂。OSI作為評(píng)價(jià)脈動(dòng)血流的一個(gè)重要參數(shù)，可反映收縮期和舒張期動(dòng)脈瘤內(nèi)血流方向快速變化引起WSS的振蕩大小，大的OSI表示W(wǎng)SS的方向變化較大。變化著方向的WSS不斷作用于瘤壁內(nèi)膜上，達(dá)到一定閾值時(shí)造成內(nèi)膜上內(nèi)皮細(xì)胞損傷、脫落，膠原纖維暴露，血小板在頸部、頂部受損的內(nèi)膜表面聚集，附著于膠原纖維上，導(dǎo)致內(nèi)膜增厚。反復(fù)內(nèi)膜損傷造成內(nèi)膜進(jìn)行性增厚，血栓形成，最終由纖維蛋白溶解反應(yīng)導(dǎo)致瘤的破裂。本研究發(fā)現(xiàn)與載瘤動(dòng)脈相比，大型動(dòng)脈瘤具有顯著偏高的OSI值，說(shuō)明大型動(dòng)脈瘤壁WSS震蕩明顯大于載瘤動(dòng)脈。因?yàn)閯?dòng)脈瘤壁比載瘤動(dòng)脈更薄，故能耐受的WSS震蕩更小、破裂風(fēng)險(xiǎn)更高。

本研究對(duì)動(dòng)脈瘤模擬分析的結(jié)果很好的詮釋了大動(dòng)脈瘤破裂風(fēng)險(xiǎn)更高的深層次血流動(dòng)力學(xué)原因。鑒于大動(dòng)脈瘤較之小動(dòng)脈瘤更容易因?yàn)轭^痛、占位壓迫、缺血性發(fā)作等破裂前先兆癥狀而被提早發(fā)現(xiàn)，以及大動(dòng)脈瘤自身存在的更大的治療難度和破裂風(fēng)險(xiǎn)，我們認(rèn)為，針對(duì)未破裂大型動(dòng)脈瘤的血流動(dòng)力學(xué)特點(diǎn)采取相關(guān)治療策略可避免動(dòng)脈瘤破裂出血后導(dǎo)致更高的治療風(fēng)險(xiǎn)和更差的預(yù)后。

［1］Wiebers DO，Whisnant JP，Huston J，et al.Unruptured intracranial aneurysms:Natural history，clinical outcome，and risks of surgical and endovascular treatment［J］.Lancet，2003，362(9378):103-110.

［2］Matsumoto KM，Oota S，Aoki M，et al.Natural history of unruptured cerebral aneurysms of the unoperated and observed cases［J］.No Shinkei Geka，2005，33(1):35-41.

［3］Ojini FI.Natural history of cerebral saccular aneurysms［J］.West Afr J Med，2005，24(1):70-74.

［4］Rinkel GJ.Natural history，epidemiology and screening of unruptured intracranial aneurysms［J］.J Neuroradiol，2008，35(2):99-103.

［5］Mocco J，Komotar RJ，Lavine SD，et al.The natural history of unruptured intracranial aneurysms［J］.Neurosurg Focus，2004，17(5):E3.

［6］Cebral JR，Castro MA，Burgess JE，et al.Characterization of cerebral aneurysms for assessing risk of rupture by using patient-specific computational hemodynamics models［J］.AJNR Am J Neuroradiol，2005，26(10):2550-2559.

［7］Shojima M，Oshima M，Takagi K，et al.Magnitude and role of wall shear stress on cerebral aneurysm:Computational fluid dynamic study of 20 middle cerebral artery aneurysms［J］.Stroke，2004，35(11):2500-2505.

［8］Luo B，Yang X，Wang S，et al.High shear stress and flow velocity in partially occluded aneurysms prone to recanalization［J］.Stroke，2011，42(3):745-753.

［9］張義森，楊新健，喬愛(ài)科，等.支架置入對(duì)顱內(nèi)動(dòng)脈分叉部動(dòng)脈瘤血流動(dòng)力學(xué)的影響［J］.中華神經(jīng)外科雜志，2010，26(3):215-219.

［10］Valencia A，Solis F.Blood flow dynamics and arterial wall interaction in a saccular aneurysm model of the basilar artery［J］.Comp Struct，2006，84(21):1326-1337.

［11］Malek AM，Alper SL，Izumo S.Hemodynamic shear stress and its role in atherosclerosis［J］.JAMA，1999，282(21):2035-2042.

［12］Alvaro V，Hernan M，Rodrigo R，et al.Blood flow dynamics in patient-specific cerebral aneurysm models:The relationship between wall shear stress and aneurysm area index［J］.Med Eng Phys，2008，30(3):329-340.

［13］Brent U，Jenn Stroud R.Numerical simulation of saccular aneurysm hemodynamics:Influence of morphology on rupture risk［J］.J Biomech，2007，40(12):2716-2722.

［14］He X，Ku DN.Pulsatile flow in the human left coronary artery bifurcation:Average conditions［J］.J Biomech Eng，1996，118(1):74-82.