0 引言

　　由于不同物質對于不同波長光波的吸收和散射的響應不同，光譜技術被廣泛應用于各種物質屬性的分析與檢測[1]。然而不同物質的非均勻性與高度散射特性是造成光譜檢測與分析結果不準確的主要原因。散射導致穿過介質的光子光程大于光源與光檢測器的距離，同時光譜分析的基本線性定律，即朗伯比爾定律(Lambert-Beer law，A=μa d)不再成立。如果能夠克服光子在介質中散射的問題，則可以極大地提高光譜檢測的精度和應用范圍。所以，在光譜檢測分析中大量使用了降低散射的方法用于提取出弱散射的光子，達到最終提高光譜檢測精度的目的[2-4]。然而提取弱散射光子的同時必然會降低檢測到的光信號的強度，即影響檢測的信噪比。所以在應用降低散射方法提高精度與保證檢測信號的信噪比之間要做權衡。

　　本文應用Monte Carlo模型[5]對3種在光譜技術中常用的降低散射的方法進行了比較分析，即偏振差法[2]、附加吸收劑法[3]和空間濾波法[4]。這3種方法都用于提取出弱散射光子，其基本思路是弱散射光子的光程等于或接近光源與光檢測器距離，從而基于朗伯比爾定律可以得到光強衰減與光子吸收系數之間的近似線性關系。本文對于線性化的效果通過比較光程的方差來判斷，對于不同方法在獲得相似的光程方差的同時的檢測光強也進行了比較。

1 原理

　　1.1 光子吸收和散射與光強衰減的關系

　　由于介質對光的散射和吸收，光波穿過介質后強度發生衰減。吸收導致透射光強度的減少，散射造成透射光在空域和時域的發散。由于光檢測器尺寸和檢測角度的有限性，散射也導致了檢測光強的衰減。標志介質對光的吸收和散射程度的參數為吸收系數（μa）和散射系數（μs），二者均與光波波長相關。

　　朗伯比爾定律是光譜分析的基本定律[1]，在不存在散射的情況下，透射光光強的衰減（Attenuation，A）與介質的吸收系數之間存在線性關系，即：

　　其中I0是入射光強，I是穿過介質的透射光強，二者的關系為：，d是光子在介質中經過的光程長度，本文以毫米為單位。對于透射光譜檢測技術，光源與光檢測器通常共軸放置，所以在沒有散射的情況下，透射光程為常數，即等于介質的厚度，A與μa成線性關系，如式(1)所示。此外，光子光程的方差為零。當存在散射的情況下，散射導致光子光程為不定長的變量。所以，散射導致A與μa關系成非線性。

　　散射導致光子光程超過光源與光檢測器間距離，從而導致光檢測器檢測到的光子在時間上展寬。通常用出射光的時間響應曲線(Temporal Point Spread Function，TPSF)反映光子在介質中光程。圖1顯示當光子穿過充滿散射介質時，透射光的TPSF時間響應曲線展寬(圖1(b))，光強衰減與吸收系數為非線性關系(圖1(c))，光程方差不再為零(圖1(d))。光程方差值越小，意味著光子散射程度越低，光強衰減與吸收系數間更加線性化。所以本文以光子光程的方差作為評價散射降低方法效果的依據。

　　1.2 降低光譜散射方法

　　偏振差法是一種有效的去除散射光的方法，廣泛應用于光譜檢測[2]、光學成像[6]、多層介質特性研究[7]等領域。偏振差法的基本思想是基于弱散射光的偏振保持特性。散射導致光子的偏振狀態呈現隨機化，因此多次散射光呈現隨機偏振狀態。所謂偏振差法，即是在光檢測端加入一個與入射光偏振方向一致的偏振片檢測出保持原有偏振態的部分光，再應用一個與入射光偏振方向正交的偏振片檢測出多次散射的部分光，將二者相減提取出弱散射光[2]。

　　在光譜測量中使用附加吸收劑法[3]是通過加入一定可控數量的吸收劑，使多次散射光進一步衰減，達到A與μa進一步線性化的目標，最終提高光譜測量的精度。本文應用附加吸收劑法濾除光程較長的多次散射光。

　　對于傳輸光譜，應用空間濾波法是基于弱散射光更靠近光軸，而強散射光更加遠離光軸的基本思想。因此共軸放置光檢測器與光源，可以獲得弱散射光。Leith[8]等人證明在傳輸模式下，通過控制光檢測器的檢測尺寸可以得到弱散射光信號，即包含無散射光子和弱散射光子的光信號。

　　1.3 Monte Carlo模型

　　由于Monte Carlo方法可以使仿真光子在不同介質中的傳輸過程更加靈活，因此得到了廣泛的應用。本文采用Monte Carlo方法[5]仿真透射光穿過充滿散射介質的試管的過程，并且對每個光子在傳輸過程中的散射、吸收以及偏振態的改變進行獨立計算，文中共仿真5×107個光子。充滿散射與吸收介質的試管尺寸設定為(40×10×10)mm3，光子從側面（40×10 mm2）中心正入射進入試管。假定以生物研究為光譜分析與檢測的對象[9]，介質的散射和吸收系數范圍設定為μs=2～9 mm-1與μa=0～0.3 mm-1，其各向異性系數（散射角余弦的平均值）g=〈cosθ〉=0.9。仿真過程中每個光子每次散射的步長、散射角度、位置、偏振態都會被計算并且記錄下來。

　　光子偏振態的改變基于Muller矩陣與Stokes矢量進行計算。對于空間濾波法的仿真，只有所有符合光檢測器檢測尺寸與檢測角度的光子會被用于計算出射光強。本文假定光檢測器為圓形，尺寸為2 mm，檢測角度設定為5°。對于附加吸收劑法的仿真，是通過計算所有出射光子的總光程，再應用朗伯比爾定律計算光強，附加吸收劑的吸收系數表示為μaadded。

2 仿真結果與分析

　　2.1 弱散射系數情況下的3種方法比較

　　圖2顯示了在散射系數μs=2 mm-1情況下的TPSF時間響應曲線(圖2(a))、吸收系數μa與光強衰減A(圖2(b))以及與方差(圖2(c))的關系曲線。每幅子圖中包含4條曲線，即未應用任何降低散射方法的透射光的響應曲線（實線），以及應用圓偏振差法（短劃線）、附加吸收劑法（點劃線）、空間濾波法（虛線）曲線。從圖2可以看出，加入降低散射方法可以得到更窄的TPSF曲線(圖2(a))，使光強衰減與吸收系數關系曲線更加線性化(圖2(b))，以及得到更低的光程方差(圖2(c))。但是，從圖2(a)和圖2(b)中更高的衰減值可以看出，降低散射法的光強衰減也更為嚴重。

　　附加吸收劑法是通過添加可控劑量的吸收劑，使多次散射的長光程光子被進一步衰減，從而實現降低光子光程總長度的目標。為了更好地與偏振差法相比較，附加吸收劑劑量控制在0.272 mm-1，即可以獲得與偏振差法相近的光程方差與線性化程度。然而，附加吸收劑法檢測到的光強更弱，光強衰減度更高，意味著信噪比更低。值得注意的是附加吸收劑法的效果取決于添加的附加吸收劑的多少，繼續增大附加吸收劑會使線性化程度更好，但是也會影響輸出的光強值。

　　空間濾波法濾除了多次散射的光子，使檢測到的光子光程近似等于光源到光子檢測器距離。所以，應用空間濾波法可以線性化光強衰減度與吸收系數關系曲線。與偏振差法相比，空間濾波法獲得的透射光TPSF盡管更窄，但是由于仿真噪音的影響，光子分散程度更高，所以方差值更大。3種方法中空間濾波法的光強衰減最為嚴重。同樣，空間濾波法的效果取決于光檢測器的尺寸。增大光檢測器尺寸，可以提高檢測光強，但是會影響提取弱散射光子的效果。

　　3種方法中，偏振差法獲得的檢測光子數目最高(圖2(a))，光強衰減度最低(圖2(b))，光程方差最低，即獲得光子發散程度最低(圖2(c))。所以，偏振差法相比其他兩種方法效果更好。

　　2.2 提高散射系數情況下3種方法比較

　　在散射系數比較低的情況下(μs=2 mm-1)，相比附加吸收劑法和空間濾波法，偏振差法可以檢測到光強更強、分散程度更低的弱散射光子。以下對提高散射系數情況下的3種方法的效果進行比較分析。圖3、圖4、圖5顯示了在散射系數μs=4、7、9 mm-1情況下，3種方法的響應。

　　隨著散射系數的增加，偏振差法仍舊可以提取出弱散射光子，但是偏振態保持特性在降低。偏振差法提取出的光子數目也隨之減少，衰減增加。圖3(a)中TPSF尾部的噪音主要來自于隨機噪聲。在高散射情況下，散射導致光子光程拉長，TPSF尾部拉長，曲線平滑度降低(圖3(a))，光程方差隨μs增加而增加(圖3(c))。同樣，散射導致光強衰減增加，如圖3(b)所示。當μs>9 mm-1之后，無法獲得滿意的偏振保持曲線，所以文中仿真到μ=9 mm-1。

　　圖4顯示了附加吸收劑法在不同散射系數下的響應。為了便于與偏振差法相比較，附加吸收劑的取值使附加吸收劑法與偏振差法取得一致的光程方差值(圖3(c)和圖4(c))。則相應附加吸收劑的值分別為：μaadded=0.199 (實線)、0.158(虛線)和0.133 (短劃線) mm-1。比較圖3(b)和圖4(b)中的光強衰減值可以看出，與偏振差法相比，附加吸收劑法提取出的光強更弱。

　　圖5顯示了空間濾波法在不同散射系數取值下的響應。同樣，隨著μs的增加，TPSF噪聲增大(圖5(a))，方差增加(圖5(c))，光強衰減度也隨之增加(圖5(b))。相比于偏振差法，空間濾波法所獲得的光強更弱。這與選取的檢測器的尺寸大小相關，為保持一致，此處仍舊取光檢測器尺寸為2 mm-1。

　　2.3 3種方法的檢測光強的比較

　　提取出弱散射光子，不可避免地降低了檢測光的強度。在應用降低散射方法的同時，還要權衡對于光檢測強度的影響。所以文中對3種方法在獲得相同光程方差值處的光檢測強度進行了比較。

　　圖6給出了3種方法在散射系數μs=2 mm-1時的光程方差與光子數目的關系曲線。光子數目對應光檢測器檢測到的光強。相同方差的情況下，光子數目值越高，說明該方法的性能越好。可以看出在當前散射系數條件下，偏振差法（實線）給出了最佳的性能，其次是附加吸收劑法（短劃線）和空間濾波法（虛線）。

3 結論

　　本文以提高光譜檢測精度為目標，以光譜分析的基本線性定律——朗伯比爾定律為理論基礎，對3種常用的降低光譜散射的方法進行比較分析。以線性化光強衰減度與吸收系數的關系曲線，縮短光程方差作為評價各種方法的標準。根據本文的仿真，相比于附加吸收劑法和空間濾波法，偏振差法給出了更好的線性化效果，更低的光程方差，可以更加有效地提高光譜檢測的精度。此外應用降低散射的方法不可避免地導致檢測光強的降低。文中比較了3種方法在相同光程方差情況下的光強，仿真結果表明，偏振差法給出了最佳的線性化光譜分析的效果，同時可以獲得最高的檢測光子數目。

　　參考文獻

　　[1] PARSON W.Modern optical spectroscopy[M].New York：Springer，2007.

　　[2] LU B，WU J W.Quantifying spectroscopic concentration ratio by polarization subtraction technique[J].Chinese Optics Letters，2013，11(S2)：S23001.1-4.

　　[3] REICH O，SCHAEL F.Investigation of the optical diffusion coefficient in turbid media by the added-absorber method[C].European Conference on Biomedical Optics，2001，4431：299-305.

　　[4] EVERALL N.Temporal and spatial resolution in transmission raman spectroscopy[J].Applied Spectroscopy，2010，64(1)：52-60.

　　[5] CHANG P C Y.Polarization discrimination for active imaging in scattering media[J].Optics Communications，1999，159(1-3)：1-6.

　　[6] ZHU Q.Experimental and theoretical evaluation of rotating orthogonal polarization imaging[J].Journal of Biomedical Optics，2009，14(14)：034006.1-10.

　　[7] STOCKFORD I M.Analysis of the spatial distribution of polarized light backscattered from layered media[J].Journal of Biomedical Optics，2002，7(3)：313-320.

　　[8] LEITH E N.Realization of time gating by use of spatial filtering[J].Applied Optics，1999，38(8)：1370-1376.

　　[9] SANDELL J L，ZHU C T.A review of in-vivo optical properties of human tissues and its impact on PDT[J].Journal of Biophotonics，2013，4(11-12)：773-787.