顱內動脈瘤伴開窗畸形血流動力學數值模擬分析的初步研究

段玉霞 陳偉建 劉一駿 劉瑾瑾 李瑞

顱內動脈瘤伴開窗畸形血流動力學數值模擬分析的初步研究

段玉霞 陳偉建 劉一駿 劉瑾瑾 李瑞

目的探討構建顱內動脈瘤（IA）伴開窗畸形三維數字化模型的可行性，并建立三維有限元模型。方法選取1例前交通動脈瘤伴開窗畸形的急診患者為研究對象，利用GE Lightspeed 64排螺旋CT獲得患者顱腦CT原始數據。應用Mimics軟件連續讀入DICOM格式的顱腦CT原始數據，采用閾值分割與手動分割相結合的圖像處理方法得到最感興趣部位，并選擇3D重建按鈕進行三維重建。再應用Geomagic studio軟件以“抽殼法”對三維重建模型的表面三角面片形態網格和血管壁截面作進一步優化處理，再現血管壁的真實形態。優化后的三維數字化模型以標準模板庫（STL）格式保存并輸出。應用CFX軟件進行三維剛性模型的數值仿真，并分析與動脈瘤生長及破裂相關的血流動力學參數。結果聯合應用Mimics軟件和Geomagic studio軟件首次成功建立了IA伴開窗畸形的三維數字化模型。模型進一步生成三維有限元模型，并作了相關血流動力學分析，形象地顯示了瘤中各區域動脈瘤壁的壁剪切力分布情況。在血流沖擊瘤頸處，壁剪切力最大，瘤體次之；在瘤頂處剪切力最小；在開窗處，壁剪切力明顯增大。結論借助可靠的CT成像技術和優化的建模方法，聯合應用Mimics軟件和Geomagic studio軟件可以成功建立IA伴開窗畸形的三維數字化模型，并分析動脈瘤各個區域和開窗處的血流動力學參數。

顱內動脈瘤開窗畸形三維數字化模型壁剪切力

顱內動脈瘤（IA）是由動脈血管壁病理性局限性擴張產生腦血管瘤樣突起的一種腦血管疾病。IA嚴重威脅著人類健康，一旦破裂，致殘率、致死率均較高。顱內動脈開窗畸形（FIA）是由于胚胎發育期間原始血管融合不全所致的一種少見的先天性動脈血管變異[1-2]，通常情況下FIA本身診斷價值有限，但是FIA合并其他腦血管疾病（尤其是IA）的影像學表現對診斷和治療的意義較大。目前，CT血管造影（CTA）能清晰地顯示顱內血管的分布，能明確診斷IA、FIA或其他血管異常性疾病。隨著現代醫學影像技術的發展，數字化成為影像學的必然趨勢[3]。目前，IA血流動力學研究主要應用工程學中的計算流體力學（CFD）方法實現，應用CFD方法建模實用性強、功能強大；建立準確、有效的三維有限元模型是關鍵所在[4]。國內外對IA三維有限元模型的研究較多，但關于IA伴開窗畸形三維有限元模型的研究較少，且相關領域的基礎研究尚處于空白階段。故本研究以顱腦螺旋CT原始數據為基本數據進行IA伴開窗畸形的三維重建與數字化模型建立，并對動脈瘤模型進行三維數值模擬，測算其血流動力學參數，為IA發生、生長及破裂原因和機制的研究提供參考。

1 材料和方法

1.1 建模環境在Windows平臺下建立腦動脈瘤伴開窗畸形的三維數字化模型。電腦配置：惠普P7370 Intel（R）Core（TM）2 Duo CPU 2.00GHz，內存2.00GB，操作系統Windows 7。軟件：醫學三維圖像生成及編輯處理軟件Mimicsl 10.0（比利時Materialise公司），逆向工程軟件Geomagic studio12.0（美國Raindrop公司）。

1.2 原始數據采集在溫州醫科大學附屬第一醫院選取1例前交通動脈瘤伴開窗畸形患者作為研究對象，研究方案取得患者知情同意，并通過本院醫學倫理學委員會審查。患者男性，48歲。使用GE Lightspeed 64排螺旋CT掃描儀獲得患者的顱腦CT原始數據。掃描范圍：自枕大孔下緣至頂結節水平，掃描線與顱底平行。基本掃描參數設置：層厚0.625mm，層間距0.5mm，螺距0.984：1，掃描速度0.5s/r，電壓100kV，電流500mA。使用高壓注射器（Medrad Stellant D CT Injector，Warrendale，USA）經肘靜脈注入碘克沙醇注射劑（注射速度4.0ml/s）。數據以醫學數字影像通信標準（DICOM）格式存儲。

1.3 動脈瘤伴開窗畸形三維數字化模型的構建應用Mimics軟件連續讀入DICOM格式的顱腦CT原始數據，正確標記圖像的方位，調整顯示閾值，使血管影像盡可能地顯示清楚。采用閾值自動分割與手動分割相結合的方法對各個層面圖像進行細致處理，得到最感興趣部位，主要包括動脈瘤、開窗畸形和載瘤動脈區域，同時去除周圍干擾的細小分支血管以及骨質。點擊3D重建按鈕，采用三維計算功能進行三維重建，三維重建模型以標準模板庫（STL）格式保存并輸出。應用Geomagic studio軟件直接讀取STL格式的三維重建模型，并采用三角面片的形式進行網格化處理，去除不必要的部分。選擇簡化多邊形，再細化多邊形，調用編輯多邊形工具，利用翻轉邊、分裂邊及收縮邊的功能進一步對模型表面三角形面片網格進行優化處理，特別是對開窗畸形處的網格進行精細的優化處理。為了更好地再現血管壁的真實形態，采用一種非常巧妙的方法-“抽殼法”對血管壁截面進行優化（圖1），首先在距離血管壁末端0.5cm處對環繞一周的三角面片進行優化處理，然后使用填充工具對截面進行填充，最后采取抽殼法復原長度為0.5cm的血管管徑，使血管截面形態更加符合實際情況，這對血流動力學的定量分析很重要。優化后的三維數字化模型的每個三角面片盡量保持等邊或直角三角形形態，且保證三角形面片角度的勻稱和邊長的合適。最后邀請2位經驗豐富的生物力學專家對優化后的模型進行評估分析（雙盲法），在達成一致意見且確定符合要求后，將其以STL格式保存并輸出。

圖1 在Geomagic studio軟件中運用“抽殼法”處理IA伴開窗畸形三維重建模型的血管壁截面

1.4 流體力學分析（1）控制方程：假定血液是層流且為黏性、不可壓縮的牛頓流體，其流體力學控制方程包括連續方程、動量方程和能量方程。本研究不考慮能量（如熱量）的傳遞，因此不考慮能量方程。同時假設不考慮重力因素，流體的基本方程是不可壓縮Navier-Strokes方程來控制。設定血液密度為1 050kg/m3，黏性系數為0.003 5Pa/s，入口血流速度為0.5m/s。（2）邊界條件：由于顱內動脈搏動不明顯，故本研究將動脈壁設定為剛性，假設血流和管壁之間沒有滑移和物質交換。（3）數值算法：Navier-Strokes方程組采用質量及動量守恒形式的有限元法來離散，網格單元為四面體型網格。通過SIMPLE算法來計算質量流量，定常計算采用時間迭代法，直至達到指定的收斂標準。

2 結果

CTA顯示前交通動脈瘤伴開窗畸形，大小8.6mm×4.1mm。本研究聯合應用Mimics軟件、Geomagic studio軟件使圖像的識別、重建及優化更加準確、有效，首次成功重建了IA伴開窗畸形的三維數字化模型（圖2a-b）。重建的IA伴開窗畸形三維數字化模型是數字解剖有限元前處理模型，不但具有良好的生物學形態和結構相似性，而且可從不同角度清晰完整地重現IA、開窗畸形和載瘤動脈區域以及三者之間的關系（圖3a-b）。根據剪切力的計算結果和相關參照文獻，將剪切力分為4個區域：高剪切力（＞3Pa）、中剪切力區（＞1.5～3Pa）、低剪切力區（＞0.5～1.5Pa）、極低剪切力區（≤0.5Pa）。本例患者IA旁邊的開窗處可見旋渦狀高切力區，瘤頸處為中切力區，瘤頂區為低切力區（圖4a-b，插頁）。

圖2 某例患者的VR重建圖和三維數字化重建模擬圖像[男，48歲，主訴為“突發意識不清5h”急診入院，臨床診斷為“蛛網膜下腔出血”；a：前交通動脈瘤伴開窗畸形的CTA重建VR圖像；b：前交通動脈瘤伴開窗畸形的三維數字化重建模型圖像（聯合應用Mimics和Geomagic studio軟件處理）；b的圖像質量明顯高于a]

圖3 前交通動脈瘤伴開窗畸形三維數字化重建模型圖像[a：頂視圖（自上而下），b：俯視圖（自下而上）]

3 討論

IA是動脈管壁病理性局限性擴張產生的腦血管瘤樣突起病，與先天性平滑肌發育不良導致局部管壁缺陷、后天性腔內壓力增高導致血管內彈力層破壞等有關。FIA是一種非常少見的先天性動脈血管變異，與原始胚胎血管融合不全有關。FIA的起點正常，位置正常，但行程部分分叉成雙支走行的表現在臨床上很少見[1-2]，且大部分FIA本身不會引起明顯臨床癥狀，通常在尸檢或CTA時偶然發現[5]。因此，FIA的臨床意義主要是其合并囊狀動脈瘤、動靜脈畸形等顱內其他血管異常時對血管介入治療路徑、神經外科治療手術入路的影響[6]。約90%的IA起自頸內動脈系統，其中起自前交通動脈區占30%～35%，大多數動脈瘤以蒂（或稱瘤頸）與載瘤動脈相連，本例患者就屬于這種情況。FIA的發生與原始胚胎血管的發育異常相關，常見于前交通動脈區，其次為椎基底動脈系統和大腦中動脈區[7-8]。根據顯微外科研究結果，前交通動脈開窗畸形的發生率約為40%[9]。本例患者開窗畸形位于前交通動脈，前交通動脈開窗畸形可能和胚胎時期大腦前、中動脈之間的原始血管網吻合支殘留及血流動力學等因素有關[5,10]。

FIA常伴有IA，相關文獻認為FIA與IA形成具有相關性[11-13]。然而，Aditya等[7]發現合并與不合并動脈瘤的患者開窗畸形發生率比較差異無統計學意義；孫曉劼等[8]亦認為FIA與IA之間不存在明確關系。FIA的結構特點是血管腔部分重復變成2個分開的內皮襯里的通道，且在開窗近遠端分叉處伴有部分中膜缺損[14]。有一種假說認為渦流造成了開窗畸形近端和遠端部分血管壁中層內膜的缺陷，導致了動脈瘤的發生[5]。根據杭曉杰等[15]對FIA與動脈瘤解剖關系的分型，本例患者動脈瘤與開窗畸形緊鄰但不位于窗內。相關研究表明，該類型的FIA有明顯的血管內渦流，會使血管內膜明顯變薄，從而使腦動脈瘤的形成機會明顯增大[16]。此外，有研究表明FIA患者有較高的動脈瘤發生率，時因為血管壁中層的缺陷同時伴隨著血流動力學壓力的增高[14]。本例患者開窗處顯示旋渦狀高切力區，表示此處血流動力學壓力增高，證實了以上結論。

本研究在充分考慮IA伴開窗畸形影像學特征和臨床研究價值的基礎上，借助可靠的CT成像技術和優化的建模方法，聯合應用Mimics軟件和Geomagic studio軟件成功建立了IA伴開窗畸形的三維數字化模型，并分析了與動脈瘤生長、破裂相關的血流動力學參數，為相關臨床血流動力學研究和手術治療措施提供了新思路。

[1] Dey M,Awad I A.Fenestration of supraclinoid internal carotid Artery andassociatedaneurysm:embryogenesis,recognition,and management[J].World Neurosurg,2011,76(6):592-595.

[2] Bentura J E,Figueiredo E G,de Monaco B A,et al.Vertebrobasilar artery junction aneurysm associated with fenestration[J].Arq Neuropsiquiatr,2010,68(2):312-314.

[3] 蘇秀云,劉蜀彬.Mimics軟件臨床應用：計算機輔助外科入門技術[J].北京:人民軍醫出版社,2011:1.

[4] Steinman D A,Milner J S,Norley C J,et al.Image-based computational simulation of flow dynamics in a giant intracranial a-neurysm[J].Am J Neuroradiol,2003,24:559-566.

[5] Sanders W P,Sorek P A,Mehta B A,et al.Fenestration of intracranial arteries with special attention to associated aneurysm and other anomalies[J].AJNR,1993,14(3):675-680.

[6] 蔣飆,何偉良,許曉英,等.有孔型椎基底動脈的MR血管成像診斷[J].中華放射學雜志,2004,38(5):525-528.

[7] Aditya B,Ritchard I A,Jeremy W,et al.Intracranial arterial fenestrations:frequency on CT angiography and association with other vascular lesions[J].Surg RadiolAnat,2008,30(5):397-401.

[8] 孫曉劼,程敬亮,金焱,等.CTA對顱內動脈開窗畸形的診斷價值[J].臨床放射學雜志2010,29(6):738-741.

[9] Gast AN,van RooijWJ,SluzewskiM.Fenestrations of the anterior communicating artery:incidence on 3D angiography and relationship to aneurysms[J].AJNR,2008,29(2):296-298.

[10] 趙兵,楊華,劉健,等.大腦前交通動脈成窗現象的解剖學特征及臨床意義[J].中國腦血管雜志,2007,4(4):168-171.

[11] 張瑩,楊新建,李海云,等.椎基底動脈開窗畸形處動脈瘤的血流動力學數值模擬分析研究[J].中國臨床神經外科雜志,2008,13(2):72-75.

[12] Peluso J P,van Rooij W J,Sluzewski M,et al.Aneurysms of the vertebrobasilar junction:incidence,clinical presentation,and outcome of endovascular treatment[J].AJNR,2007,28(9):1747-1751.

[13] Sun Z K,Li M,Li M H,et al.Fenestrations accompanied by intracranial aneurysms assessed with magnetic resonance angiography[J].NeurolIndia,2012,60(1):45-49.

[14] Hacein B L,MuszynskiC A,Varelas P N.Saccular aneurysmassociated with posterior cerebral artery fenestration manifesting as a subarachnoid hemorrhage in a child[J].AJNR,2002,23:1291-1294.

[15] 杭曉杰,程鳴,楊亞旭,等.顱內動脈開窗畸形伴腦動脈瘤的MRA診斷價值[J].中國醫學計算機成像雜志,2013,19(4):355-358.

[16] Tasker A D,Byme J V.Basilar artery fenestration in association with aneurysms of the posterior cerebral circulation[J].Neuroradiology,1997,39:185-189.

Three dimensional digital model of hemodynamics in in tracranial aneurysms with fenestration

DUAN Yuxia,CHEN Weijian,LIU Yijun,et al.
Department of Radiology，the First Affiliated Hospital of Wenzhou Medical University,Wenzhou 325000，China

ObjectiveTo construction a three-dimensional digital model of hemodynamics in intracranial aneurysm with fenestration.MethodsAn emergent patient of anterior communicating artery aneurysm with fenestration admitted in the First Affiliated Hospital of Wenzhou Medical University was selected as the study object.The patient underwent brain CT scan with GE Lightspeed 64-channel multidetector CT scanner.Mimics software was used to read consistent CT raw data of DICOM format directly and the three-dimensional function was used to reconstruct the model by threshold segmentation combined with manual segmentation method and to obtain the most important region of interest.Then Geomagic studio software was used to further optimize surface meshes with triangle patch shape and sections of vessel wall of the three-dimensional reconstruction model;in particular,the'Shell Method'was utilized for optimizing sections of vessel wall in order to reflect the real nature of blood vessel.The optimized three-dimensional digital model was stored and sent out by STL format,and then the hemodynamic parameters related to the growth and rupture of aneurysms were analyzed.ResultsThe three-dimensional digital model of intracranial aneurysm with fenestration was reconstructed successfully.The further three-dimensional finite element model and related hemodynamics were analyzed.The wall shear stress was given in the distribution of various parts of aneurysm.In the window areas,the wall shear stress increased considerably.ConclusionThe three-dimensional digital model of intracranial aneurysm with fenestration has been reconstructed successfully,with this model the hemodynamics of each region and the window of the aneurysm can be analyzed.

Intracranial aneurysm Fenestration Three-dimensional digital model The wall shear stress

10.12056/j.issn.1006-2785.2017.39.23.2016-1519

溫州市科技計劃項目（Y20140027）

325000溫州醫科大學附屬第一醫院放射科（段玉霞、陳偉健、劉瑾瑾、李瑞）；浙江省人民醫院放射科（劉一駿）

陳偉建，E-mail：wyyycwj@163.com

2016-09-27）

（本文編輯：陳丹）