王耀東，陳桂波，張燁，韓穎，譚勇，蔡紅星

（長春理工大學 理學院，長春 130022）

基于Davies理論的粗糙面光散射特性研究

王耀東，陳桂波，張燁，韓穎，譚勇，蔡紅星

（長春理工大學 理學院，長春 130022）

基于Davies理論，利用BRDF描述散射特性，系統考察了BRDF隨角度、波長及粗糙度等參數的變化規(guī)律。數值結果表明，粗糙面的BRDF在形態(tài)上有一定的類鏡面散射特性同時在鏡面方向出現“高光”區(qū)域；入射天頂角越大則散射越趨于鏡面反射；波長越小，粗糙面的相對粗糙程度越大，則Davies理論模型非鏡面反射部分比重也越大。結果表明研究可為復雜實體散射分析與應用提供一定的理論依據和指導。

Davies理論；雙向反射分布函數；粗糙面

研究粗糙面散射特性［1］，進而掌握其光散射規(guī)律是對復雜物體散射［2］規(guī)律研究的基礎和前提［3］，而表面散射特性通常用雙向反射分布函數（Bidirectional Reflectance Distribution Function，BRDF）來描述。

BRDF最早由Nicodemus于1970年提出［4］。各國學者在理論上、儀器設備研制的應用［5］上進行了大量的研究。理論模型方面先后有Davies于1954年假定表面高度起伏服從高斯分布并應用基爾霍夫近似提出Davies模型［6］，Torrance和Sparrow于1967年由基于微平面分布理論并假定粗糙面由V型凹槽構［7］成提出Torrance-Sparrow模型［8］，Maxwell和Beard于1973年提出包含表層和內層反射貢獻的漆面材料的Maxwell Beard模型［9］，Phong于1975提出基于經驗模型且主要在計算機圖形學中作為光照模擬的Phong模型［10］，Ove Steinvall于2000年提出能夠很好地描述許多不同類型材質的單站散射特性的四參數單站BRDF模型［11］，吳振森等人于2002年對Torrance-Sparrow模型進行修正，將粗糙表面的高度起伏分布和斜率分布推廣到非高斯分布提出五參數半經驗統計模型［12］，Sun Yinlong于2007年提出基于微面元高斯分布統計概率而建立的，可以很好地匹配光滑表面和粗糙表面的Sun BRDF模型［13］，等等。

其中Davies理論模型基于電磁場Maxwell方程，其物理意義明確且應用范圍相對較廣［14，15］，因此本文采用Davies模型研究粗糙面光散射特性。本文首先介紹了描述散射特性的BRDF及Davies理論模型，然后應用Davies理論模型計算了粗糙面的BRDF進而分析得到了粗糙面的散射特性和規(guī)律。

1 理論模型

1.1 BRDF

BRDF是描述表面散射特性的一種方式，其定義為［16］:

如圖1所示，θi和θr為入射和出射天頂角，φi和φr為入射和出射方位角，公式（1）中dEi表示入射光在目標面元造成輻照度的增量，dLr表示該面元在接收到這個輻照亮度增量后在出射方向亮度的增量。Li表示入射輻射亮度，dωi為入射立體角。即BRDF為對于樣品表面某一待分析的微面元，其沿出射方向（θr,φr）（探測方向）的亮度增量與引起此亮度增量的入射方向（θi,φi）照度增量的比值。

圖1 出射入射光束的幾何示意圖

1.2 Davies理論模型

Davies理論模型主要是通過假設不規(guī)則表面的高度起伏在其平均值附近呈高斯分布，并基于基爾霍夫近似來描述電磁能量從粗糙表面散射的數學公式。其適用于條件下的反射，式中σ為表面粗糙度均方根，a為表面自相關長度［17］。模型的BRDF包含兩個部分，分別是鏡面反射部分和非鏡面反射部分。

其中鏡面反射部分表示為:

非鏡面反射部分表示為:

式中，ρdh,s是平滑表面的方向半球反射率，可由菲涅爾公式決定［14］。

假定φi=0°，θi=30°，λ=500nm，ρdh,s=28.70%，σ=0.07μm，a=2.2μm，計算取定參數［15］后粗糙表面Davies理論模型的BRDF，其散射形態(tài)如圖2所示。

圖2 粗糙表面Davies理論模型的散射形態(tài)圖

圖2在直角坐標系下，待研究散射材料平放于xOy面上。入射光從圖中箭頭方向入射（θi=30°，φi=0°）。設A(x,y,z)為圖形表面上一點，O為原點，可有向量其中θr，φr標示r→（即 OA）向量的方向，同時也為fr(θi,φi;θr,φr)中的出射方向，分別為天頂角和方位角。的大小也同時表示這個出射方向下的 fr值。

由圖2可見，粗糙面Davies理論模型形態(tài)上表現為一定的類鏡面散射特性。由于粗糙面Davies理論模型分為鏡面反射部分和非鏡面反射部分，鏡面反射部分只有鏡面反射方向存在散射光，所以在圖中鏡面反射方向點處 fr函數不連續(xù)，出現極窄的尖峰，也就是出現“高光”區(qū)域［18］。

2 結果與分析

由公式（2）和（3）可知，因變量有θi，φi，θr和φr參變量有λ（波長），ρdh,s（方向半球反射率），σ（表面粗糙度均方根）和a（表面自相關長度）。假定φi=0°，θi=30°，λ=500nm，ρdh,s=28.70%，σ=0.07μm，a=2.2μm，計算取定參數［15］后粗糙表面Davies理論模型的BRDF散射形態(tài)如圖2。

由于研究參數對于BRDF的影響首先得選定其他參數的值。這里波長選取中間值500nm。方向半球反射率、表面粗糙度均方根和表面自相關長度選取Kapon樣品的參數值ρdh,s=28.70%（對應波長500nm下的值）、σ=0.0700μm和a=2.2000μm［15］。

2.1 Davies理論模型中入射角度的影響

假定φi=0°，λ=500nm，ρdh,s=28.70%，σ=0.07μm，a=2.2μm后，計算粗糙表面Davies理論模型的BRDF散射形態(tài)，不同入射角下BRDF散射形態(tài)如圖3，可見在入射天頂角增大時，非鏡面反射部分的形態(tài)保持基本不變的同時整體縮小，而鏡面反射部分則逐漸增強。由于整體形態(tài)基本不變，所以下面討論隨其他參數的影響規(guī)律時取入射天頂角θi=30°這一定值，又由于模型本身假定材料并無方向性，所以入射方位角取任一值均可，為方便起見這里取φi=0°。

圖3 Davies理論模型 fr隨θi變化關系圖（參數：φi=0°；λ=500nm；ρdh,s=28.70%；σ=0.07μm；a=2.2μm）

2.2 Davies理論模型中λ及ρdh,s的影響

波長λ對于粗糙表面Davies理論模型的BRDF的影響，從公式本身來說，可以先探討取定其他參數后 fr隨波長λ的變化關系，假定φi=0°，θi=30°，ρdh,s=28.70%，σ=0.07μm，a=2.2μm后，計算粗糙表面Davies理論模型的BRDF散射形態(tài)，不同入射角下BRDF散射形態(tài)如圖4。

圖4 Davies理論模型 fr隨λ變化關系圖（ρdh,s=28.70%）（參數：φi=0°；θi=30°；ρdh,s=28.70%；σ=0.07μm；a=2.2μm）

由公式（2）和（3）可見，如果其他參數均不變，則隨著波長λ增長，鏡面反射部分逐漸增強；而非鏡面反射部分則分為系數部分和指數部分，系數部分逐漸減弱，指數部分逐漸增強，整體的結果則是逐漸減弱。隨波長的具體變化關系如圖5所示。

圖5 Davies理論模型各項值在鏡面反射方向隨波長的關系（參數：φi=φr-180°=0°；θi=θr=30°；ρdh,s=28.70%；σ=0.07μm；a=2.2μm）

從圖5中可見，Davies理論模型在波長較短時，模型中非鏡面反射項起主要作用。隨著波長的增大，非鏡面反射項的值整體在下降，鏡面反射項的值則緩慢增長，所以鏡面反射項的作用將逐漸增強。這主要是由于隨著波長的增大，粗糙面相對波長的粗糙情況越來越小，從而散射中非鏡面項的比重越來越小了。

從Davies理論模型公式出發(fā)，可以人為討論只改變波長而其他參數均不變的情況下的影響規(guī)律。但是，參數中方向半球反射率ρdh,s并非與波長λ無關的參數，對于某一確定的材料將ρdh,s=ρdh,s(λ)，本文選取Kapon樣品為例討論，其ρdh,s(λ)隨波長λ取值如表1。

表1 Kapon樣品ρdh,s(λ)隨λ的關系［15］

所以實際上針對已知的特定研究材料，ρdh,s不再是獨立參變量，這時假定 φi=0°，θi=30°，σ=0.07μm，a=2.2μm，計算粗糙表面Davies理論模型的BRDF，其散射形態(tài)如圖6。

圖6 Davies理論模型 fr隨λ變化關系圖（ρdh,s=ρdh,s(λ)）（參數：φi=0°；θi=30°；σ=0.07μm；a=2.2μm）

圖6與圖4對比可見，當考慮ρdh,s=ρdh,s(λ)后，散射整體形態(tài)隨波長由遞減轉化成了先增大后減小關系，而這個轉變主要由于ρdh,s在模型中對整個散射起到整體縮放的作用。從而，由于ρdh,s在短波段過渡到中間波段時有快速增大過程而后基本保持不變，使得原本單調遞減的整體散射形態(tài)在短波段到中間波段有了遞增，而后ρdh,s基本不變整體形態(tài)又出現再次遞減。

2.3 Davies理論模型粗糙面參數的影響

2.3.1 Davies理論模型σ的影響

表面粗糙度均方根σ是表面輪廓高度起伏的標準差，其值大體反映了表面輪廓起伏的高低。圖7展示假定φi=0°，θi=30°，λ=500nm，ρdh,s=28.70%，a=2.2μm后，計算粗糙表面Davies理論模型的BRDF在不同表面粗糙度均方根σ下散射形態(tài)，隨著表面粗糙度均方根σ增大，表面輪廓起伏越來越大，表面也越粗糙，所以鏡面反射項逐漸減小，而非鏡面反射項逐漸增大。

圖7 Davies理論模型fr隨σ變化關系圖（參數：φi=0°；θi=30°；λ=500nm；ρdh,s=28.70%；a=2.2μm）

選取鏡面反射方向（θr=30°，φr=180°）的 fr進行更為具體的分析，如圖8可見隨著表面粗糙度均方根σ的增大，fr逐漸由鏡面反射部分為主逐漸轉變?yōu)榉晴R面反射部分為主。

圖8 鏡面反射方向下 fr隨σ變化關系圖

2.3.2 Davies理論模型a的影響

表面自相關長度a，是表面自相關函數R(τ)下降到其最大值的倍的τ［17］。一般用來估計隨機表面的特征波長，如果輪廓上兩點的距離大于自相關長度，則可以認為這兩點是互不相關的，用來說明表面輪廓的隨機性。

假定φi=0°，θi=30°，λ=500nm，ρdh,s=28.70%，σ=0.07μm后，計算粗糙表面Davies理論模型的BRDF在不同表面自相關長度a下散射形態(tài)由圖9可見，隨著表面自相關長度a增大，非鏡面反射項逐漸增大。從公式（2）可見，表面自相關長度a并不影響鏡面反射項。圖10則選取鏡面反射方向（θr=30°，φr=180°）的 fr進行更為具體的分析。

圖9 Davies理論模型 fr隨a變化關系圖（參數：φi=0°；θi=30°；λ=500nm；ρdh,s=28.70%；σ=0.07μm）

圖10 鏡面反射方向下 fr隨a變化關系圖

3 結論

基于Davies理論粗糙面的光散射形態(tài)表現出類散射特性，有較強的方向散射特性，其BRDF包含鏡面反射項和非鏡面反射項，BRDF在鏡面反射方向出現“高光”區(qū)域。

入射方位角不影響B(tài)RDF的整體散射形態(tài)和大小，隨著入射天頂角的增大，鏡面反射部分逐漸增強，非鏡面反射部分逐漸減弱。在不考慮方向半球反射率隨波長變化時，波長增大主要通過影響非鏡面反射部分的系數部分使其逐漸減弱，而鏡面反射部分則緩慢增強。在考慮方向半球反射率隨波長變化時，方向半球反射率相當于對鏡面和非鏡面部分整體進行縮放。

表面粗糙度均方根和自相關長度增大都會使Davies理論模型的非鏡面反射部分增強，而前者增大使鏡面反射部分減弱，后者不影響鏡面反射部分。

［1］ 于昊.散射法表面粗糙度測量的數學模型分析［J］.長春理工大學學報:自然科學版，2006，29（1）:109-112.

［2］ 任智斌，姜會林，付躍剛，等.微球體與微橢球體光散射特性的研究［J］.長春理工大學學報:自然科學版，2006，29（1）:28-31.

［3］ 郭立新，王蕊，吳振森.隨機粗糙面散射的基本理論與方法［M］.北京:科學出版社，2010:1-4.

［4］ Nicodemus F E.Geometrical considerations and nomenclature for reflectance［J］.Applied Optics，1970，9（6）:1474-1475.

［5］ 孫成明，趙飛，袁艷.基于光譜的天基空間點目標特征提取與識別［J］.物理學報，2015，64（3）:034202-1-7.

［6］ Shen Y J，Zhang Z M，Tsai B K.Bidirectional reflectance distribution function of rough silicon wafers［J］.International Journal of Thermophysics，2001，22（4）:1311-1326.

［7］ 劉宏，朱京平，王凱.基于隨機表面微面元理論的二向反射函數幾何衰減因子修正［J］.物理學報，2015，64（18）:184213-1-6.

［8］ Torrance K E，Sparrow E M.Theory for off-specular reflection from roughened surfaces［J］.Journal of the Optical Society of America，1967，57（9）: 1105-1114.

［9］ Maxwell J R，Beard J，Weiner S，et al.Bidirectional reflectance model validation and utilization［R］.Ann Arbor:Environmental Research Institute of Michigan，1973:10-35.

［10］ Bui Tuong Phong.Illumination for computer generated pictures［J］.Communications of ACM，1975，18（6）:311-317.

［11］ Ove Steinvall.Effects of target shape and reflection on laser radar cross sections［J］.Applied Optics，2000，39（24）:4381-4391.

［12］ Wu Zhensen，Xie Donghui，Xie Pinhua，et al.Modeling reflectance function from rough surface and algorithms［J］.Acta Optica Sinica，2002，22（8）: 897-901.

［13］ Sun Yinlong.Statistical ray method for deriving reflection models of rough surfaces［J］.Journal of theOpticalSociety ofAmerica，2007，24（3）: 1105-1114.

［14］ Shen Y J，Zhang Z M，Tsai B K，et al.Bidirectional reflectance distribution function of rough silicon wafers［J］.International Journal of Thermo physics，2001，22（4）:1311-1326.

［15］ Bai Lu，Wu Zhensen，Cao Yunhua，and Huang Xun. Spectral scattering characteristics of space target in near-UV to visible bands［J］.Optics Express，2014，22（7）:8515-8524.

［16］ Dai Yongjiang.Laser Radar Technology［M］.Publish House of Electronics Industry，Beijing:2010: 409-412.

［17］ 袁長良，丁至華，武文堂.表面粗糙度及其測量［M］.北京:機械工業(yè)出版社，1989:16-17.

［18］ 黃冰峰.物體表面反射特性測定方法的研究［D］.武漢:華中科技大學，2007:8-9.

Study on Scattering Characters of the Rough Surface Based on Davies Theory

WANG Yaodong，CHEN Guibo，ZHANG Ye，HAN Ying，TAN Yong，CAI Hongxing
（School of Science，Changchun University of Science and Technology，Changchun 130022）

In this paper，regularity of scattering has been investigated using BRDF based on Davies theory changing with angle，wavelength，roughness，etc.The numerical simulation showed that the shape of rough surface scattering has a certain scattering of mirror-like and has a highlight area in the specular reflection direction.The scattering shape will be more like specular reflection with the zenith angle of incident increase.The off-specular component will be more important with the wavelength decrease，relative roughness increase，in Davies model.This article research results can provide necessary theory and guidance for complex entity scattering.

Davies theory；BRDF；rough surface

O436

A

1672-9870（2016）05-0045-06

2016-04-12

王耀東（1990-），男，碩士研究生，E-mail:gannannu@163.com

陳桂波（1979-），男，副教授，E-mail:guibochen@126.com