HUIYING

真實(shí)解剖下光傳播與皮層劑量預(yù)測(cè)不準(zhǔn)概述

經(jīng)顱光生物調(diào)制通過近紅外光照射頭皮，試圖調(diào)節(jié)大腦皮層神經(jīng)元活動(dòng)，但光在穿越頭皮、顱骨、腦脊液、灰質(zhì)等多層組織時(shí)，其實(shí)際傳播路徑與到達(dá)皮層的劑量長(zhǎng)期缺乏精確量化?，F(xiàn)有研究多采用簡(jiǎn)化幾何模型，如分層平板、同心球或忽略腦溝回的三維結(jié)構(gòu)，部分模型甚至將腦脊液層合并或省略，導(dǎo)致無法準(zhǔn)確描述光在真實(shí)解剖結(jié)構(gòu)下的散射與分布特征。這種簡(jiǎn)化使得光在腦脊液層可能產(chǎn)生的特殊光學(xué)效應(yīng)被忽視，光在皮層表面的空間分布、穿透深度與組織吸收比例均缺乏可靠預(yù)測(cè)，成為tPBM臨床參數(shù)優(yōu)化與療效評(píng)估的根本瓶頸。

HUIYING

光傳播中的“光暈效應(yīng)”及其形成機(jī)制

在高精度模型中，光的傳播過程展現(xiàn)出以往簡(jiǎn)化模型無法捕捉的重要特征。如圖1A系列所示，當(dāng)光從頭皮進(jìn)入顱骨時(shí)，光吸收區(qū)域基本保持與光源形狀一致的圓形分布；然而當(dāng)光進(jìn)一步傳播至腦脊液層時(shí)，如圖1A3所示，光吸收區(qū)域開始出現(xiàn)明顯的“光暈效應(yīng)”——即光能向周圍擴(kuò)散，形成比光源面積更大的環(huán)形分布。這一效應(yīng)在腦脊液層最為顯著，如圖1A7的局部放大圖所示，并進(jìn)一步延伸至灰質(zhì)表面，如圖1A4與A8所示，使皮層實(shí)際受照范圍大于光源投影區(qū)域。論文通過第3.4節(jié)的參數(shù)影響分析進(jìn)一步驗(yàn)證了這一現(xiàn)象的機(jī)制：如圖5D系列所示，只有當(dāng)腦脊液恢復(fù)其真實(shí)的極低散射系數(shù)（μ_s‘≈0.001 cm?1）時(shí)，光暈效應(yīng)才會(huì)重現(xiàn)；若將腦脊液光學(xué)屬性設(shè)為與灰質(zhì)相同，如圖5A至C系列所示，光暈效應(yīng)隨之消失。這表明腦脊液因其獨(dú)特的弱散射特性，在光傳播中起到了“光擴(kuò)散層”的作用，顯著改變了光在皮層表面的空間分布模式，這是以往簡(jiǎn)化模型所無法揭示的關(guān)鍵機(jī)制。

圖1：800 nm光源在F3腦區(qū)刺激下的光傳播特征分析

圖1系統(tǒng)展示了光從頭皮表面穿透至腦組織的全過程。圖1 A1至A5分別呈現(xiàn)了光在頭皮、顱骨、腦脊液、灰質(zhì)、白質(zhì)表面的吸收輻照度分布，圖中標(biāo)注了10、100、1000 mW/cm2三種輸入條件下各組織的最大吸收值，直觀反映了光能的逐層衰減。圖1 A6至A8對(duì)關(guān)鍵區(qū)域進(jìn)行了局部放大，清晰顯示了腦脊液層出現(xiàn)的“光暈效應(yīng)”及其向灰質(zhì)表面的延伸。圖2B1至B4從冠狀面、軸狀面、矢狀面展示了光吸收的深度分布，圖1 B5至B7則提供了對(duì)應(yīng)的結(jié)構(gòu)MRI切片作為對(duì)照。圖1 C1通過沿激光中心垂直向內(nèi)的切割線繪制了光吸收隨深度變化的曲線，量化了頭皮約65%、顱骨約33.9%的吸收比例。圖1 C2將三種功率密度下的吸收曲線歸一化后對(duì)比，顯示光吸收與輸入功率呈線性關(guān)系。圖1全面揭示了光在真實(shí)解剖結(jié)構(gòu)中的傳播規(guī)律與空間分布特征。

圖5：組織光學(xué)屬性對(duì)光傳播空間分布的影響分析

圖5通過參數(shù)影響分析，驗(yàn)證了腦脊液在光暈效應(yīng)形成中的關(guān)鍵作用。圖5A系列顯示將所有組織設(shè)為灰質(zhì)光學(xué)屬性時(shí)，光吸收呈均勻分布，無光暈效應(yīng)。圖5B系列顯示恢復(fù)頭皮真實(shí)屬性后，光暈效應(yīng)仍未出現(xiàn)。圖5C系列顯示進(jìn)一步恢復(fù)顱骨真實(shí)屬性后，光暈效應(yīng)依然缺失。圖5D系列顯示只有當(dāng)腦脊液恢復(fù)其真實(shí)極低散射系數(shù)后，光暈效應(yīng)才重現(xiàn)，并延伸至下層組織。圖5A6至D6的截面圖進(jìn)一步對(duì)比了四種條件下的光吸收深度分布。圖5系統(tǒng)證明了腦脊液因其弱散射特性，是光暈效應(yīng)的關(guān)鍵來源，揭示了以往簡(jiǎn)化模型無法捕捉的重要光學(xué)機(jī)制。

HUIYING

熱安全邊界的量化需求與現(xiàn)有局限

tPBM在照射過程中，光能被組織吸收后轉(zhuǎn)化為熱量，可能引起頭皮和腦組織溫度升高，但現(xiàn)有臨床研究采用的功率密度差異懸殊，從幾十到上千 mW/cm2不等，缺乏統(tǒng)一的熱安全標(biāo)準(zhǔn)。美國(guó)國(guó)家標(biāo)準(zhǔn)學(xué)會(huì)針對(duì)皮膚表面設(shè)定了330 mW/cm2的激光安全限值，但該限值是否適用于腦組織、不同功率密度下實(shí)際溫升幅度如何、溫度隨刺激時(shí)間的動(dòng)態(tài)變化規(guī)律如何，均缺乏系統(tǒng)性的量化預(yù)測(cè)。現(xiàn)有研究多將光傳播與熱效應(yīng)分開處理，未能在真實(shí)解剖模型中實(shí)現(xiàn)光-熱耦合分析，使得臨床參數(shù)選擇多依賴經(jīng)驗(yàn)或體外實(shí)驗(yàn)，存在潛在熱風(fēng)險(xiǎn)隱患。因此，建立一個(gè)能夠同時(shí)預(yù)測(cè)光分布與溫度變化的耦合模型，對(duì)于界定tPBM的熱安全邊界至關(guān)重要。

HUIYING

高精度多物理場(chǎng)耦合模型的構(gòu)建

為解決上述問題，研究團(tuán)隊(duì)構(gòu)建了一套高分辨率（1mm）的三維有限元仿真模型。模型基于MNI152平均腦模板，精細(xì)分割出頭皮、顱骨、腦脊液、灰質(zhì)、白質(zhì)等組織，完整保留了腦溝回、眼窩等真實(shí)解剖細(xì)節(jié)。如圖2所示，模型能夠清晰展示各組織的三維結(jié)構(gòu)，并在F3腦區(qū)位置設(shè)置圓柱形光源（半徑20 mm，面積約12 cm2）。模型采用輻射傳輸方程的擴(kuò)散近似描述光傳播，同時(shí)耦合Pennes生物傳熱方程計(jì)算溫度變化，實(shí)現(xiàn)了光吸收與熱源之間的物理耦合。組織光學(xué)參數(shù)與熱學(xué)參數(shù)均取自文獻(xiàn)，考慮了血液灌注散熱與皮膚表面對(duì)流散熱。模型支持穩(wěn)態(tài)與瞬態(tài)分析，能夠輸出光吸收分布、功率分配比例、溫度空間分布及其隨刺激時(shí)間的變化曲線，為tPBM的劑量設(shè)計(jì)與安全評(píng)估提供了高保真度的數(shù)值平臺(tái)。

圖2：三維頭模型組織分割與激光光源位置示意圖

圖2展示了基于MNI152模板構(gòu)建的高精度三維頭模型，通過圖像分割技術(shù)將頭部組織精細(xì)劃分為頭皮、顱骨、腦脊液、灰質(zhì)和白質(zhì)等多個(gè)層次。圖2中第一行呈現(xiàn)了各組織的三維重建結(jié)果，清晰顯示了眼窩、顱骨結(jié)構(gòu)以及連續(xù)的腦脊液層。第二行則展示了激光光源在F3腦區(qū)的放置位置，以及光源下方皮層表面的解剖細(xì)節(jié)。圖2為后續(xù)光傳播與熱效應(yīng)仿真提供了真實(shí)的幾何基礎(chǔ)，體現(xiàn)了模型在解剖細(xì)節(jié)上的高保真度。

HUIYING

面向臨床應(yīng)用的仿真實(shí)驗(yàn)設(shè)計(jì)與結(jié)果

本研究雖未開展實(shí)際臨床試驗(yàn)，但完全基于臨床真實(shí)場(chǎng)景構(gòu)建了仿真實(shí)驗(yàn)。研究以F3腦區(qū)為刺激靶點(diǎn)，采用連續(xù)波激光模式，模擬了三種臨床常見功率密度（10、100、1000 mW/cm2）下的光傳播與熱效應(yīng)，并以ANSI 330 mW/cm2皮膚安全限值作為對(duì)照，同時(shí)考察了20分鐘內(nèi)溫度隨刺激時(shí)間的動(dòng)態(tài)變化。

光傳播方面：如圖1所示，模型三維可視化地展示了光從頭皮穿透至灰質(zhì)、白質(zhì)的全過程。圖1 A1至A5分別展示了光在頭皮、顱骨、腦脊液、灰質(zhì)、白質(zhì)表面的吸收輻照度分布，圖中標(biāo)注了三種輸入條件下各組織的最大吸收值。在100 mW/cm2輸入下，頭皮最大吸收為34.9 mW/cm2，顱骨為1.04 mW/cm2，灰質(zhì)表面僅為1.03 mW/cm2，約占輸入值的1%。圖1 B2至B4分別從冠狀面、軸狀面、矢狀面展示了光吸收的深度分布，圖1 C1的線圖分析顯示，光強(qiáng)降至0.01 mW/cm2時(shí)，距灰質(zhì)表面約113 mm。表III進(jìn)一步量化了功率分布：當(dāng)輸入功率為1.7174 W（對(duì)應(yīng)100 mW/cm2）時(shí)，到達(dá)灰質(zhì)的功率為0.0978 W，占比約5.7%。這些數(shù)據(jù)為臨床中“皮層實(shí)際劑量”的設(shè)定提供了直接換算依據(jù)。

熱效應(yīng)方面：如圖3所示，模型預(yù)測(cè)了不同功率密度下各組織的溫度分布。圖3B1至B5顯示，10 mW/cm2刺激下溫度分布與基線幾乎無差異；而圖3D1至D5顯示，1000 mW/cm2刺激下，頭皮局部溫度顯著升高，最大溫升達(dá)3.76°C（圖3E1），灰質(zhì)溫升為0.57°C（圖3E4）。相比之下，100 mW/cm2刺激下頭皮最大溫升為0.38°C，灰質(zhì)溫升僅為0.06°C。圖4A通過沿激光中心垂直向內(nèi)的“切割線”展示了溫度隨深度的變化曲線，顯示溫度在頭皮層快速上升后迅速下降。圖4B匯總了三種功率密度下各組織的溫度范圍，顯示在1000 mW/cm2時(shí)，頭皮最高溫度已達(dá)40.7°C。圖4C與4D繪制了頭皮與腦組織最高溫度隨輸入功率密度的變化曲線，并標(biāo)示了ANSI 330 mW/cm2限值下的對(duì)應(yīng)溫度（頭皮升溫1.24°C，腦組織升溫0.19°C），結(jié)果顯示溫度升高與功率密度并非完全線性，腦組織升溫速率低于頭皮。

圖3：不同功率密度下的組織溫度分布與溫升變化

圖3對(duì)比了無刺激基線狀態(tài)與三種功率密度刺激下的組織溫度分布。圖3A1至A5顯示了無刺激時(shí)各組織的基礎(chǔ)溫度（37°C）。圖3B1至B5、C1至C5、D1至D5分別對(duì)應(yīng)10、100、1000 mW/cm2刺激下的溫度分布，圖中標(biāo)注了各組織的最高溫度值。圖3E1至E5則呈現(xiàn)了相對(duì)于基線的溫度升高值，顯示在1000 mW/cm2條件下，頭皮最大溫升達(dá)3.76°C，灰質(zhì)溫升為0.57°C；而100 mW/cm2條件下頭皮溫升僅0.38°C，灰質(zhì)溫升0.06°C。圖3直觀展示了熱效應(yīng)隨功率密度的變化規(guī)律，并明確了不同劑量下的熱安全邊界。

圖4：溫度隨深度與功率密度的變化規(guī)律

圖4從定量角度分析了溫度變化的空間分布與劑量依賴性。圖4A通過沿激光中心垂直向內(nèi)的切割線繪制了100 mW/cm2條件下溫度隨深度的變化曲線，顯示溫度在頭皮層快速上升后迅速下降，進(jìn)入顱骨后逐漸平緩。圖4B匯總了10、100、1000 mW/cm2三種功率密度下各組織的最低與最高溫度范圍，顯示在1000 mW/cm2時(shí)頭皮最高溫度已達(dá)40.7°C。圖4C與4D分別繪制了頭皮與腦組織最高溫度隨輸入功率密度的變化曲線，并標(biāo)示了ANSI 330 mW/cm2限值下的對(duì)應(yīng)溫度（頭皮升溫1.24°C，腦組織升溫0.19°C），結(jié)果顯示腦組織升溫速率低于頭皮，且溫度升高與功率密度呈非線性關(guān)系。

瞬態(tài)分析結(jié)果：如圖6所示。圖6A展示了在100 mW/cm2刺激下，不同時(shí)間點(diǎn)（20、120、600、1200秒）溫度沿深度的變化，可見頭皮溫度隨刺激時(shí)間迅速上升，而腦組織溫度上升緩慢。圖6B顯示，所有組織的溫度均在約600秒（10分鐘）后趨于平臺(tái)期，此后繼續(xù)延長(zhǎng)刺激時(shí)間不會(huì)產(chǎn)生顯著額外溫升。例如，刺激1200秒（20分鐘）時(shí)，頭皮溫度穩(wěn)定在37.43°C，腦組織溫度為37.09°C，與穩(wěn)態(tài)分析結(jié)果高度一致。

圖6：100 mW/cm2功率密度下溫度隨刺激時(shí)間的瞬態(tài)變化

圖6展示了20分鐘刺激時(shí)間內(nèi)溫度的動(dòng)態(tài)演變規(guī)律。圖6A繪制了20、120、600、1200秒四個(gè)時(shí)間點(diǎn)下溫度沿深度的變化曲線，顯示頭皮溫度隨刺激時(shí)間迅速上升，而腦組織溫度上升緩慢。圖6B呈現(xiàn)了各組織最高溫度隨時(shí)間的演變曲線，顯示所有組織的溫度均在約600秒（10分鐘）后趨于平臺(tái)期，此后繼續(xù)延長(zhǎng)刺激時(shí)間不會(huì)產(chǎn)生顯著額外溫升。例如，刺激1200秒（20分鐘）時(shí)，頭皮溫度穩(wěn)定在37.43°C，腦組織溫度為37.09°C。圖6為臨床刺激時(shí)長(zhǎng)的選擇提供了熱安全方面的量化依據(jù)。

安全性結(jié)論：100 mW/cm2及以下功率密度在常規(guī)刺激時(shí)長(zhǎng)內(nèi)熱效應(yīng)安全可控；1000 mW/cm2導(dǎo)致頭皮溫升3.76°C、灰質(zhì)溫升0.57°C，已接近或超過細(xì)胞興奮性變化閾值（≥0.5°C），且頭皮溫度超過臨床常規(guī)安全上限，故不推薦使用。

HUIYING

總結(jié)

本研究通過構(gòu)建高精度、多物理場(chǎng)耦合的三維有限元模型，系統(tǒng)揭示了tPBM在真實(shí)人腦解剖結(jié)構(gòu)下的光傳播規(guī)律與熱安全性邊界。模型首次在1 mm分辨率下完整保留了腦脊液與腦溝回結(jié)構(gòu)，發(fā)現(xiàn)腦脊液因其極低散射系數(shù)，在光傳播中產(chǎn)生顯著的“光暈效應(yīng)”，改變了光在皮層表面的空間分布，這一機(jī)制在以往簡(jiǎn)化模型中無法捕捉。研究量化了光能在大腦各組織中的實(shí)際分配：僅約1%的輸入光能到達(dá)灰質(zhì)，頭皮吸收約65%的光能。在熱安全方面，100 mW/cm2及以下功率密度在10分鐘刺激時(shí)長(zhǎng)內(nèi)溫升安全可控，而1000 mW/cm2存在明確熱風(fēng)險(xiǎn)。研究結(jié)果為tPBM的臨床參數(shù)選擇、光源設(shè)計(jì)與安全性評(píng)估提供了量化依據(jù)，體現(xiàn)了高保真計(jì)算建模在神經(jīng)調(diào)控器械研發(fā)中的支撐價(jià)值。