粗糙球体和锥体目标激光散射非相干分量比

思黛蓉，王明军,2*,刘永勤，眭晓林

(1.西安理工大学 自动化与信息工程学院，西安 710048; 2.西安理工大学 陕西省智能协同网络军民共建重点实验室，西安 710126; 3.中国电子科技集团有限公司第十一研究所 固体激光重点实验室，北京 100015)

引 言

激光雷达探测技术在军用和民用工程上有广泛的应用，例如激光测距、测速，激光2维、3维成像[1]，激光距离多普勒域成像和合成孔径成像[2]等。所以，对于粗糙体目标激光散射特征研究是其不可缺少的关键技术之一[3-5]。研究人员将入射激光假设为平面波[6-7]来进行激光雷达目标成像的探测和识别，主要是利用粗糙面与激光散射特性及机理去探测目标[8-9]。当激光束从粗糙的目标表面反射时，散射回波包含了物体的物理尺寸、轮廓等重要的信息，对目标的识别、跟踪及粗糙表面光学特性的反演具有重要的意义[10]。

在粗糙面散射中，众多学者利用基尔霍夫理论研究了随机粗糙表面散射特性，其中，激光入射在粗糙凸体目标的散射特征是科研人员长期以来关注的研究热点之一。ISHIMARU等人[4,11]采用双频互相关函数研究随机脉冲波传播和随机粗糙面脉冲波散射的问题。SALAMI等人[12]根据粗糙表面的统计特性，给出了有阴影效应的相干散射强度大于没有阴影效应的散射强度，非相干散射强度的变化与此相反。CHEN等人[5-6]研究脉冲平面波对2维随机粗糙介质表面的散射特性，给出了激光波束散射双频互相关函数的解析表达式。ZHANG等人以粗糙球体和圆柱体[13]为例，推导了任意形状物体在远场中的相干与非相干分量的表达式[14-15]。YANG等人[16]利用目标散射特性计算旋转粗糙凸物体的时间相关函数，给出了旋转圆锥体目标散斑相关函数。LIU等人[7,17]利用平面波研究了空间锥形目标的激光散射特性，计算了圆锥的相干散射强度，并分析了散射强度随不同参量的变化规律。作者所在团队[18-19]利用高斯波束入射在粗糙椭球体目标上，推导了散射场量二阶统计矩及非相干分量比的表达式。ZHA等人在高斯波束[20]的基础上，研究了在各向异性媒质中圆柱的散射特性，获得了各向异性圆柱的内场以及近场的归一化散射强度分布图[21]。目前，对于高斯分布任意粗糙体目标的理论建模仍以平面波为主[15,19]。然而，对于平面波入射下粗糙体目标的非相干分量比研究较少，本文中研究了激光波束入射下空间球体和锥体目标的散射特征。

本文中以标量Helmholtz积分方程和粗糙体目标散射理论研究为基础，计算了2维粗糙体目标散射场非相干分量比。针对粗糙球体和圆锥体散射特性，详细分析了粗糙目标非相干分量比的影响因素，为后期计算更复杂粗糙体目标理论建模及散射特性等问题的研究奠定基础。

1 粗糙物体平面波散射特性

Fig.1 Schematic diagram of scattering of rough objects[6]

根据标量Helmholtz方程[6]，散射场表达式为:

(1)

(2)

观察点P处的散射场表达式为:

(3)

当表面任一点的主曲率半径远大于入射波长，可采用切平面近似。粗糙面上任一点rΣ′上的散射场和法向导数可分别为:

E′(rΣ′)=(1+Ri)Ei(rΣ′)

(4)

(5)

式中，Ri为菲涅耳反射系数。采用远场近似，(1)式可简化为:

(6)

式中，令V=ki-ks，M=ki+ks,则:ki=k(sinθicosφi,sinθisinφi,-cosθi),ks=k(sinθscosφs,sinθssinφs,cosθs)。

假设目标为导体，粗糙目标的远场散射场可简化为:

(7)

利用远场近似，(7)式中距离为:

(8)

(7)式中指数项的分母|rs-r′|可以近似为R，R是观察点P与光滑面S′原点之间的距离，则粗糙目标的远区散射场可表示为:

exp(iV·r′)dS′

(9)

2 粗糙球体和锥体目标的激光散射特性

2.1 球体表面激光散射的非相干分量比

根据(9)式可知，平均散射场可表示为:

exp(iV·r′)dS′

(10)

散射场场量的二阶统计特征即互相关函数表达式[15]为:

exp[iV(r1′-r2′)]×

〈exp[iVn1′ξ(r1′)]〉〈exp[-iVn2′ξ(r2′)]〉

(11)

式中,r1′为目标原点到光滑面S′的矢径，r2′为目标原点到光滑面S″的矢径，S″是粗糙目标散射场共轭后的光滑面，k1为散射场Es的自由空间波数，k2为散射场Es共轭后的自由空间波数，ki,1表示散射场Es的入射波数矢量，n1′为光滑面S′的单位外法向矢量，ki,2表示散射场Es共轭后的入射波数矢量，n2′为光滑面S″的单位外法向矢量，ξ(r1′)是沿光滑面S′外法线方向的随机高度起伏，ξ(r2′)是沿光滑面S″外法线方向的随机高度起伏。

散射场量强度的非相干分量可表示为:

(ki,2·n2′)exp[iV(r1′-r2′)]×

(χf-χ1χ2)

(12)

式中,K=exp[i(k1-k2)R]/(2πR)2，粗糙面起伏函数ξ(r′)服从高斯分布，粗糙面的粗糙度为σ，相关长度为lc，令Vz=V·n′，则高斯分布的1阶、2阶特征函数[14]分别为χ1=exp(-k12σ2Vz2/2),χ2=exp(-k22σ2Vz2/2),χf=exp[-k2σ2Vz2(1-〈ξ1ξ2〉)]。

定义比值系数γ12为:

(13)

由(11)式、(12)式可定义γ12的数学含义，其为高斯粗糙球体非相干分量与相干分量的比值。如图2所示，入射平面波沿着-z轴方向入射在粗糙球体上，球体中心位于坐标系的原点。假设平面波入射波长λ=1.06μm，球体半径为a，σ为粗糙面的粗糙度，lc为相关长度，入射角θi=0°，方位角φi=0°，其轴线与ks的夹角为θs，即ks=(sinθs,0,cosθs)。

Fig.2 Schematic diagram of plane wave incident rough sphere scattering

2.2 圆锥体表面激光散射的非相干分量比

假设平面波入射在圆锥体侧面，取圆锥底面中心为坐标原点，如图3所示，圆锥体的底面半径为b，高为h，半锥角为α，方位角为φ，则锥面方程[7]为:

Fig.3 Schematic diagram of cone light scattering

(14)

则圆锥曲面上的微元ds=(h-z)tanα·secαdφdz。假设入射波在平面xOz内，即方位角为0°。圆锥的散射波矢量为ks=sinθs(cosφs·x+sinφs·y)+cosθs·z，法向矢量为n′=cosαcosφ·x+cosαsinφ·y+sinα·z，入射波矢量为ki=-sinθi·x-cosθi·z，x,y,z分别为单位正交基向量。

(15)

令R(rs,r′)=|rs-r′|，R是观察点P与光滑面S′原点之间的距离，则平均散射场可表示为:

(16)

利用高频近似，则(16)式的积分可简化为:

{h·exp[ikh[tanα[sinθscos(φ-φs)+sinθicosφ]-cosθs-cosθi]]}

(17)

故锥体目标的相干散射分量为:

Ic=|〈Es〉2|

(18)

根据粗糙面散射理论及(15)式、(16)式，粗糙圆锥体对平面波的非相干散射分量为:

exp[i(kR(rs,r″)+ki·r″)]×

{〈exp[ikV·n′ξ(r′)-ikV·n″ξ(r″)]〉-

〈exp[ikV·n′ξ(r′)]〉×

〈exp[ikV·n″ξ(r″)]}dS′dS″

(19)

由图4可知，引入新的积分变量，假设S′面上任一点处的曲率半径近似等于圆锥体底面半径b，且曲率半径远大于波长λ和粗糙面的相关长度lc，并满足条件Vlc2/a≪1，则V·n′≈V·n″，对dR的积分可以远近似用在r′处的切平面S⊥′(r′)内的积分代替，即积分变量dR≈dR⊥。

Fig.4 Tangential plane approximate calculation of the incoherent scattering intensity of the object[7]

当观察点位于散射场远场时，令ρ(rs,r′)≈R(rs,r′)+ki·r′，将(19)式进一步简化为:

[χ2(V·n′;R⊥)-χ2(V·n′)]dR⊥

(20)

式中，χ2(V·n′;R⊥)=〈exp{ikV·n′[ξ(r′)-ξ(r′+R⊥)]}〉,χ(V·n′)=〈exp[ikV·n′ξ(r′)]〉为1维、2维特征函数。

根据ISHIMARU[4]的描述，平面波对粗糙物体的非相干散射强度可以看作是平面波入射场的单位面积粗糙平面非相干射强度的叠加，则:

If=A·If,0

(21)

式中,A=(2L)2为被照射面积，If,0为单位面积粗糙面的非相干散射强度。

对于圆锥的非相干散射有:

(22)

式中,S(θi,θs)为遮蔽函数，且:

S(θi,θs)=S(θi)S(θs)=

(23)

(23)式表示曲面积分只在照射区域进行。If,0为单位面积粗糙平面的非相干散射强度，ns为粗糙物体平均表面外法向单位矢量。将(17)式代入(20)式得:

(24)

由(18)式、(20)式和(24)式可计算出粗糙圆锥体目标的非相干散射分量比。

3 数值计算及分析

计算结果表明，由图5a中给出3种不同高度起伏条件下对应的γ12值可以发现，粗糙度越大，非相干分量越多，相干分量越少，球体表面越粗糙；由图5b可以看出，球体半径a对非相干分量比的影响较大，在散射角为45°时存在峰值，且随着球体半径a的变大，峰值高度变小，相干分量越多，即表面越光滑，该结论同参考文献[19]中的所述相符。

Fig.5 Normalized incoherent component ratio of rough sphere scattering as a function of scattering angle

图6中给出不同粗糙度条件下，粗糙球体散射归一化非相干分量比随相关长度的变化情况。从图中可以看出,当粗糙度增大时，对应曲线上升的速度越快，即非相干分量占总分量的比重越大。由图7可以发现,当粗糙度增加时,球体半径a不同对非相干分量比有影响，且随着球体半径a增加，球体非相干分量比逐渐变小。

Fig.6 Normalized incoherent component ratio of rough sphere scattering varies with correlation length

Fig.7 Normalized incoherent component ratio of rough sphere scattering varies with roughness

假设圆锥的计算波长λ=1.06μm，圆锥体底面半径b=1m，半锥角α=15°，入射光与z轴的夹角θi=30°，抛光铝材料的复折射率为2.43+10.7i，镀金包覆薄膜材料的复折射率为1.2277+10.3768i，观察方位角φs=30°，粗糙度分别为σ=0.05μm,σ=0.1μm和σ=0.2μm。

图8中给出了粗糙圆锥体目标随散射角的变化情况。可以发现，粗糙度对非相干分量比有很大影响，随着粗糙度的不同，非相干分量比的峰值位置也有所不同，但均位于镜反射方向上。图9中给出圆锥体的底面半径b越大，非相干分量比就越大，函数γ12的下降速度就越快。

Fig.8 Variation of normalized incoherent component ratio of cone scattering under different roughness

Fig.9 Variation of normalized incoherent component ratio of cone scattering at different radii

图10中给出了两种不同目标材料随散射角的变化情况。由图可看出，抛光铝材料圆锥体非相干分量比大于镀金包覆薄膜材料。由图11可以看出，不同的粗糙度下，随着圆锥入射角的增加，函数γ12的峰值位置发生变化，与图8中的变化趋势相似。图中结论与参考文献[19]和参考文献[7]中所述结论相符。

Fig.10 Normalized incoherent component ratio of cones of different materials as a function of scattering angle

Fig.11 Normalized incoherent component ratio of cones with different roughness as a function of incident angle

4 结 论

基于物理光学近似，在平面波激光照射下，从理论上研究2维高斯分布粗糙体目标的远区散射场及双频互相关函数。以粗糙球体和圆锥体目标为例，通过数值方法分析了球体、圆锥体散射的非相干分量比随散射角、粗糙度、半径及目标材料等因素的变化情况。

(1)散射角对粗糙球体散射的非相干分量比有影响，当粗糙度越大，目标非相干散射分量在总散射分量中的比重较大，且与参考文献[18]中的结论相符。

(2)粗糙球体半径a变大，非相干分量比变小，相干分量越大。

(3)粗糙度的不同对圆锥体目标散射的非相干分量比的峰值位置有影响，且峰值均位于镜反射方向上。

(4)圆锥体的底面半径b变化与非相干分量比成正比，对于不同的目标材料，可以发现抛光铝材料的非相干分量比大于镀金包覆薄膜材料。

综上所述，研究高斯粗糙面散射特性实质上是采用定量方法讨论各种粗糙体目标的非相干分量比受哪些因素影响。通过研究粗糙体目标的散射特性，对开展更复杂的目标激光散射统计特性、以及激光散斑探测和识别等问题提供了理论依据。