唐媛媛，蘇英鋒，關(guān)慶捷，劉迎曦2，孫秀珍

（1.大連醫(yī)科大學(xué)附屬第二醫(yī)院耳鼻咽喉頭頸外科，遼寧 大連 116027；2.大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室，遼寧 大連116021）

鼻腔結(jié)構(gòu)異常者有限元數(shù)值模型建立及氣流場特征分析

唐媛媛1，蘇英鋒1，關(guān)慶捷1，劉迎曦2，孫秀珍1

（1.大連醫(yī)科大學(xué)附屬第二醫(yī)院耳鼻咽喉頭頸外科，遼寧 大連 116027；2.大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室，遼寧 大連116021）

目的通過對健康人和鼻腔結(jié)構(gòu)異常（鼻中隔偏曲）者鼻腔氣道的三維重建，計算機建立有限元數(shù)值模型，分析2者鼻腔氣道氣流場分布特征，進而從生物力學(xué)角度探討鼻腔結(jié)構(gòu)與功能的關(guān)系。方法根據(jù)鼻中隔偏曲者及健康北方人（各20例）的鼻腔CT薄層掃描結(jié)果，應(yīng)用表面重建有限元剖分的方法對其進行三維重建，數(shù)值模擬流經(jīng)鼻腔氣道的氣流特征。結(jié)果鼻中隔偏曲者雙側(cè)鼻腔氣流流量分布以氣道寬敞側(cè)為主，且最大氣流量集中在寬敞側(cè)總鼻道中部；鼻腔氣道壓強在偏曲最為明顯處下降最快，占?xì)獾揽倝簭姴畹?9.65%。健康人雙側(cè)鼻腔氣流量受鼻周期影響，以一側(cè)鼻腔為主，最大氣流量分布在總鼻道中部和下部；鼻腔氣道壓強于鼻閾處下降最快，約占?xì)獾揽倝簭姴畹?8.78%。結(jié)論用計算機模擬的方法建立鼻腔結(jié)構(gòu)的數(shù)值模型，分析其解剖結(jié)構(gòu)異常與氣流場特征之間的關(guān)系，是一種新的研究鼻腔結(jié)構(gòu)、功能與疾病相關(guān)性的可行的科學(xué)手段，可以為以優(yōu)化氣流路徑、改變氣流分布為目的的手術(shù)治療方案及手術(shù)前后的個性化評估提供參考。

鼻中隔偏曲；三維重建；數(shù)值模型；鼻腔結(jié)構(gòu)異常

鼻腔最重要的功能即為呼吸，呼吸氣流主要受到鼻腔結(jié)構(gòu)的影響［1～3］。鼻腔氣流的研究是了解鼻腔結(jié)構(gòu)和功能之間關(guān)系的紐帶。鼻腔結(jié)構(gòu)復(fù)雜，目前常用的鼻腔功能檢查手段（如鼻聲反射和鼻阻力測定）不足以反映鼻腔通氣功能的分區(qū)段變化［4］，隨著相關(guān)臨床檢測設(shè)備和生物數(shù)值模擬研究方法的進步，已有很多學(xué)者采用計算流體動力學(xué)的研究方法對鼻腔氣流場及溫、濕度場等進行研究。研究顯示計算流體力學(xué)在反應(yīng)鼻腔生理功能基本特征的同時，也可用于分析鼻腔結(jié)構(gòu)變化對功能的影響，因此，本研究應(yīng)用計算流體動力學(xué)方法（computational fluid dynamics，CFD）對鼻中隔偏曲患者及健康人的鼻腔氣流流場進行研究，探討鼻腔氣流與鼻腔解剖結(jié)構(gòu)和疾病之間的相互影響，為進一步闡述鼻腔疾病的發(fā)病機制提供參考。

1 材料與方法

1.1 材料

鼻中隔偏曲組（實驗組）：選擇鼻中隔偏曲伴有不同程度鼻塞、頭痛、鼻出血等癥狀的患者20例，年齡20～60歲，其中男14例，女6例，均為北方人，病史10～20年，平均16.7年。經(jīng)前鼻鏡及鼻內(nèi)鏡檢查鼻中隔呈“C”形偏曲（可包含一側(cè)或雙側(cè)鼻甲代償性改變），并除外鼻息肉及其他鼻腔疾病者。健康人組（對照組）：選擇20～60歲的20名健康北方人為對照組，男12例，女8例。經(jīng)詳細(xì)詢問病史，證明既往無慢性上呼吸道疾病病史，無鼻部外傷及手術(shù)史，近3～6個月無上呼吸道急性病史，應(yīng)用前鼻鏡及鼻內(nèi)鏡檢查無鼻科疾病體征或無明顯解剖異常。

1.2 方法

1.2.1 多層螺旋CT掃描：對受試者應(yīng)用螺旋CT（德國西門子公司，型號：Somatom Emotion16）進行鼻部連續(xù)掃描（圖1、圖2），范圍自鼻尖至鼻咽后壁，層厚為3 mm，軟組織窗（窗寬2 000，窗位400）。應(yīng)用CT機自帶的軟件將獲得的DICOM影像學(xué)數(shù)據(jù)轉(zhuǎn)化為BMP格式圖像。

圖1 鼻中隔偏曲者鼻腔冠狀位連續(xù)CT圖像Fig.1 The continuous coronal CT image of nasal septum deviation

圖2 健康人鼻腔冠狀位連續(xù)CT圖像Fig.2 The continuous coronal CT image of healthy nasal cavity

1.2.2 三維重建：應(yīng)用Matlab軟件將上述獲得的BMP格式圖像數(shù)字化，用自編譯的程序識別鼻腔氣道的邊界，并將其導(dǎo)入Ansys10.0軟件中，借助該軟件的前處理功能人工選取每一層面能夠反映氣道結(jié)構(gòu)特征的關(guān)鍵點連成線，進一步構(gòu)建面、體，選用四面體單元對體進行有限元剖對鼻腔氣流通道進行表面三維重建。

1.2.3 計算氣流場：應(yīng)用Ansys10.0后處理功能結(jié)合黏性流體運動方程Navier-Stockes方程對鼻腔氣流場進行計算。邊界條件設(shè)定：（1）鼻腔氣道壁無滑移邊界條件，即Vs=0 m/s；（2）前鼻孔為1個標(biāo)準(zhǔn)大氣壓（P=101 325 Pa）；（3）鼻咽部氣流速度（m/s）根據(jù)國人生理狀態(tài)下呼吸頻率和潮氣量數(shù)值范圍，選擇呼吸周期（T）=3 s，潮氣量為500 mL，即通氣量為10 L/min，假設(shè)呼吸時間相同，呼吸過程模擬為正弦變化模式，由模型計算得到鼻咽部氣道面積，再通過測量1次呼吸量，計算得到鼻咽部氣流平均速度。

根據(jù)三維重建模型，得出鼻咽部的橫截面積、鼻腔前后部的壓差和鼻咽部的氣流流量，根據(jù)計算公式［5］：鼻阻力＝鼻腔前后的壓差/氣體流速，計算鼻阻力值。同時我們利用鼻阻力計對上述受試者進行鼻阻力檢測。

1.3 統(tǒng)計學(xué)處理

應(yīng)用統(tǒng)計分析軟件SPSS 12.0對實驗組和對照組鼻腔氣道不同部位的氣流分布情況分別進行統(tǒng)計分析，實驗數(shù)據(jù)用x±s表示，雙側(cè)對應(yīng)部位均采用兩樣本t檢驗。P＜0.05認(rèn)為差異有統(tǒng)計學(xué)意義。

2 結(jié)果

2.1 氣流流速及氣道面積分布情況

鼻中隔偏曲者鼻腔氣流速度分布云圖見圖3。氣道寬敞側(cè)的氣流流速較氣道狹窄側(cè)快；鼻腔氣道寬敞側(cè)截面積較大。2者均于總鼻道中部最為明顯，見表1。

健康人雙側(cè)鼻腔氣流速度分布云圖見圖4。受鼻周期影響，20例受試者中有16例（80%）以一側(cè)鼻腔通氣為主。雙側(cè)鼻腔氣流流速主要通氣側(cè)較非主要通氣側(cè)快，尤以總鼻道中部及下部為明顯。鼻腔氣道截面積雙側(cè)差異無統(tǒng)計學(xué)意義（P＞0.05）。見表2。

圖3 鼻中隔偏曲者鼻腔氣流分布云圖Fig.3 The airflow distribution in patients with nasal septum deviation

表1 實驗組雙側(cè)氣道對應(yīng)部位氣流量、氣流速度、氣道面積對比Tab.1 The comparison of fluence，airflow rate，area of airway of two airways in the experiment group

圖4 健康人鼻腔氣流分布云圖Fig.4 The airflow distribution of healthy nasal cavity

2.2 鼻腔氣道壓強分布情況

2組受試者鼻腔氣道壓強均從前鼻孔到鼻咽部逐漸下降。實驗組中于氣道狹窄側(cè)鼻中隔偏曲最為明顯處（即鼻腔氣道最狹窄處）氣壓下降最明顯（圖5），約占?xì)獾揽倝簭姴钪档?9.65%，從前鼻孔到鼻咽部氣道壓強差值的變化范圍為（248±12）Pa。對照組中鼻閾附近氣壓下降最快（圖6），約占?xì)獾揽倝簭姴钪档?8.78%，從前鼻孔到鼻咽部的氣道壓強差值變化范圍為（131±91）Pa。

圖5 鼻中隔偏曲者鼻腔壓強分布云圖Fig.5 The nasal cavity pressure of nasal septum deviation

表2 對照組雙側(cè)氣道對應(yīng)部位氣流量、氣流速度、面積對比Tab.2 The comparison of fluence，airflow rate，area of airway of two air passage in the control group

2.3 鼻阻力檢測

據(jù)重建模型計算數(shù)據(jù)，實驗組鼻腔氣道阻力為（0.488±0.27）kPa/（L·s），對照組為（0.164±0.03）kPa/（L·s），兩組差異有統(tǒng)計學(xué)意義（t=2.382，P＜0.05）。通過鼻阻力計測量，得出實驗組鼻腔氣道阻力為（0.479±0.09）kPa/（L·s），對照組為（0.175± 0.11）kPa/（L·s），二者差異有統(tǒng)計學(xué)意義（t=2.327，P＜0.05）。101 194 101 210 101 226 101 242 101 258 101 275 101 291 101 303 101 323 101 339

圖6 健康人鼻腔鼻腔壓強分布云圖Fig.6 The nasal cavity pressure of healthy nasal cavity

3 討論

上氣道的諸結(jié)構(gòu)均與功能有密切的相關(guān)性，鼻腔的主要功能就是對吸入的氣流進行加溫、加濕、過濾灰塵以及嗅覺等。鼻腔結(jié)構(gòu)異常主要表現(xiàn)在鼻中隔和下鼻甲，鼻中隔偏曲的程度和下鼻甲的大小決定鼻腔有效氣流通路的截面積，研究人鼻腔鼻竇的生物力學(xué)數(shù)值模型，是從生物力學(xué)角度來研究上氣道解剖結(jié)構(gòu)與相關(guān)疾病的關(guān)系，從而服務(wù)于臨床。為了驗證計算機數(shù)值模型的可靠性，在實驗過程中，我們對所有受試者根據(jù)數(shù)值模型計算出的鼻阻力數(shù)據(jù)，應(yīng)用鼻阻力計進行鼻阻力測量進行驗證，兩種方法均得出實驗組鼻阻力較對照組增大，而2組經(jīng)兩種方法測得的鼻阻力差異無統(tǒng)計學(xué)意義（P＞0.05）。因此，可以認(rèn)為本研究建立的數(shù)值模型基本反映出鼻腔結(jié)構(gòu)異常者和健康人的鼻腔氣道特征。

本實驗數(shù)值模擬的健康人上氣道鼻腔氣流流場與孫秀珍等［6］研究的正常國人鼻腔氣流場特征相符，同時研究發(fā)現(xiàn)，鼻中隔偏曲者和健康人呼吸時同樣會受到鼻周期的影響，而與鼻周期相比，鼻腔解剖結(jié)構(gòu)異常對流經(jīng)鼻腔的氣流分布影響更為明顯。鼻中隔偏曲者雙側(cè)鼻腔容積不對稱，呼吸時通過氣道狹窄側(cè)鼻腔的氣流流量明顯減少，尤其以該側(cè)總鼻道中部減少更為明顯。隨著氣流量減少，氣流流速也隨之減慢，使該側(cè)中鼻道的鼻竇口正常引流功能受到影響，從而導(dǎo)致該側(cè)鼻竇炎的發(fā)生。為增加其氣流量，該側(cè)下鼻甲多有不同程度的萎縮變小，表現(xiàn)為黏膜變薄、骨質(zhì)吸收并伴隨有向外側(cè)移位等，從而減少該側(cè)鼻腔氣道阻力，我們認(rèn)為這是人體鼻腔的一種自適應(yīng)性變化［7］。對于氣道寬敞側(cè)，由于較大的氣流長期慢性沖擊，該側(cè)鼻腔組織代償性增生，中鼻道黏膜病變的程度和累及范圍甚至超過對側(cè)，同樣導(dǎo)致該側(cè)鼻竇炎發(fā)病率增高。而且研究中我們發(fā)現(xiàn)，鼻中隔高位偏曲者較低位偏曲者在偏曲最為明顯處的總鼻道中部雙側(cè)鼻腔氣流流量及流速差異更大，這一點可能是此前學(xué)者［8］提到的高位偏曲者鼻竇炎發(fā)病率較高的原因之一。

鼻腔氣流場特征研究同樣可以應(yīng)用于患者嗅覺改變的病因分析。鼻中隔偏曲者氣道狹窄側(cè)氣流流量減少、流速減慢，且氣流由于狹窄氣道的分流作用導(dǎo)致局部湍流增多［9］，同時因嗅裂區(qū)本身的解剖特點，使得進入該區(qū)的氣流減少更為明顯。根據(jù)Zhang等［10，11］研究上氣道模型中微小粒子在氣道中的傳送和沉積作用，以及曾敏捷等［12］通過氣體動力學(xué)模型研究微顆粒在人體上呼吸道中的運動沉積，探討顆粒沉積率與氣流流速和顆粒慣性大小的關(guān)系等研究結(jié)果，我們有理由認(rèn)為隨著進入該側(cè)嗅裂區(qū)的氣流分布情況變化，使得進入并沉積于嗅區(qū)的引起嗅覺的顆粒物質(zhì)明顯減少甚至缺失，從而導(dǎo)致嗅覺的減退甚至喪失。

因此，我們認(rèn)為對病態(tài)上氣道結(jié)構(gòu)進行數(shù)值模擬，分析其解剖結(jié)構(gòu)異常與氣流場特征之間的關(guān)系，可用于評估上氣道結(jié)構(gòu)—功能—疾病之間互為因果的發(fā)生機制，進而為以優(yōu)化氣流路徑、改變氣流分布為目的手術(shù)治療提供參考。

［1］Yu S，Liu YX，Sun XZ，et al.Influence of nasal structure on the distribution of airflow in nasal cavity［J］.Rhinology，2008，46（2）：137-143.

［2］Wexler D，Segal R，Kimbell J.Aerodynamic effects of inferior turbinate reduction：computational fluid dynamics simulation［J］.Arch Otolaryngol Head Neck Surg，2005，131（12）：1102-1107.

［3］劉迎曦，于申，孫秀珍.數(shù)值模擬鼻甲的切除對鼻腔內(nèi)氣體流場的影響［J］.生物醫(yī)學(xué)工程學(xué)雜志，2008，25（6）：1315-1318.

［4］曹潔，陳寶元，朱寶玉，等.阻塞性睡眠呼吸暫停綜合征患者鼻阻力變化［J］.中華結(jié)核和呼吸雜志，2000，23（12）：725-726.

［5］黃選兆，汪吉寶.實用耳鼻咽喉科學(xué)［M］.北京：人民衛(wèi)生出版社1998：112-141.

［6］孫秀珍，劉迎曦，蘇英鋒，等.鼻腔氣道三維重建和氣流流場的數(shù)值模擬與分析［J］.臨床耳鼻咽喉頭頸外科雜志，2007，21（23）：1057-1059.

［7］孫秀珍，唐媛媛，劉迎曦，等.鼻中隔偏曲者鼻腔結(jié)構(gòu)自適應(yīng)改變的特征分析［J］.中華耳鼻咽喉頭頸外科雜志，2008，43（5）：351-354.

［8］姜菲菲，閻艾慧，王強，等.鼻中隔偏曲與鼻竇炎關(guān)系的CT影像學(xué)研究［J］.解剖科學(xué)進展，2010，16（2）：135-137.

［9］蘇英鋒，孫秀珍，劉迎曦，等.鼻腔加溫功能特征及其與氣流場關(guān)系的研究［J］.力學(xué)學(xué)報，2012，44（3）：607-613.

［10］Zhang Z，Kleinstreuer C.Airflow structures and nano-particle deposition in a human［J］.J Computational Physics，2004，198（1）：178-210.

［11］Zhang Z.Aerosol transport and deposition in a triple bifurcation bronchial airway model with local tumors［J］.Inhalation Toxicol，2002，14（11）：1111-1133.

［12］曾敏捷，胡桂林，樊建人，等.微顆粒在人體上呼吸道中運動沉積的數(shù)值模擬［J］.浙江大學(xué)學(xué)報（工學(xué)版），2006，40（7）：1164-1167.

（編輯武玉欣）

EstablishmentofFinite Element Modeland Analysisof Airflow Characteristicsin Patientswith NasalStructuralAbnormalities

TANGYuan-yuan1，SUYing-feng1，GUANQing-jie1，LIUYing-xi2，SUNXiu-zhen1
（1.Department of Otolaryngology Head and Neck Surgery，The Second Affiliated Hospital of Dalian Medical University，Dalian 116027，China；2.Department of Engineering Mechanics，Dalian University ofTechnology，Dalian 116023，China）

ObjectiveTo analyze the airflow characteristics and investigate the relationship of the structure and the function of nasal cavity by the three-dimensionalreconstruction ofnasalairway ofpatientswith structuralabnormalities（nasalseptum deviation）and healthy people and the establishmentoffinite elementmodelby computer.MethodsOn the basisofCTimaging ofthe nasalcavity in patients with structuralabnormalities（nasal septum deviation，n=20）and healthy people（n=20），three-dimensional reconstruction of nasal airway was conducted by resurfacing finite element subdivision to simulate the characteristics of airflow in nasal cavity.ResultsThe airflow mainly went through the commodious side of the nose in patients with nasal septum deviation and the maximum fluence appeared in the middle part of meatus nasi communis.The airway pressure decreased most significantly in the most flank-curvature part of nasal septum deviation，accounting approximately 79.65%of the total pressure.In healthy people，the bilateral airflow was affected by nasal cycle and was mainly characterized by one nasal cavity，and the maximum fluence was observed in the middle and the inferior part of meatus nasi communis.The airway pressure decreased most significantly in limen nasi，accounting approximately 58.78%of the total pressure.ConclusionNumerical modeling of nasal cavity can be used to analyze the relationship between the nasal structuralabnormalities and the airflow characteristics，which is a scientific method to analyze the association ofnasalstructure and function with disease and can be used forpre-and post-operative individualevaluation ofoperative therapeutic regimen targeting atoptimizing airway and altering airflow distribution.

nasal septum deviation；three-dimensional reconstruction；numerical model；nasal structural abnormalities

R318.5

A

0258-4646（2015）03-0209-05

國家自然科學(xué)基金（10672036；11072055；10902022；11472074）

唐媛媛（1980-），女，主治醫(yī)師，碩士.

孫秀珍，E-mail：sunxiuzhen001@163.com

2014-11-13

網(wǎng)絡(luò)出版時間：