【doc】 声场分离技术及其在近场声全息中的应用.doc

声场分离技术及其在近场声全息中的应用第54卷第2期2005年2月1000329o/20o5/54(o2)/078909物理ACTAPHYSICASINICAVo1.54,No.2,February,20052005Chin.Phys.Soc.声场分离技术及其在近场声全息中的应用于飞陈剑李卫兵陈心昭(合肥工业大学动态测试中心,合肥230009)(2003年8月4日收到;2004年5月24日收到修改稿)提出空间声场分离技术,突破了近场声全息(NAH)的应用局限.它们的局限在于全息面一侧的声场必须是自由声场,即要求所有的声源必须位于另一侧.利用波数域内的波场外推理论及声压的标量叠加原理,建立起声场分离技术的双全息面实现方法,利用波数域内的Euler公式及粒子振速的矢量叠加原理,建立起该技术的单全息面实现方法.该技术的一个突出优点是在具有背景噪声的全息测量情况下,可以消除背景噪声对全息变换结果的影响.理论的推导表明该技术方法的正确性,而仿真算例和实验则显示该技术的可行性和有效性.关键词:声全息,波数域,声场分离,背景噪声PACC:4320,43631.引言关于声学逆问题的研究近来非常活跃,一类是利用测得的散射场信息进行物形反演1,另一类利用测得的二维声压信息重建三维空间声场.近场声全息技术的诞生,打破了传统全息技术成像分辨率的波长限制,可以精确地识别噪声源和可视化空间声场.j.通过记录两维全息面上的声压数据来重建整个三维空间内的声压场,粒子振速矢量场及声强矢量场,全息技术的这一优良特性受到众多研究人员的关注.为了降低全息面声压需同时测量的条件,Loyau和Pascal提出了基于声强测量的宽带声全息技术(BAHIM),该技术利用相位梯度场同声能量场之间的关系,通过二维切向声强和均方声压的测量来重建全息面上的复声压.为了使声全息技术能够应用于重建任意形状物体的表面振速,声压等物理量,可以通过Helmhohz积分理论建立声源表面声学参量同场中声压参量之间的关系h.这种声全息技术被称作为基于边界元法(BEM)的声全息,其本身也存在着一些不可避免的缺点,如表面积分方程中奇异积分的处理,解的非唯一性的处理等问题.同时全息重建的实体对象必须是封闭的实体,需要在声场中布置大量的测量点,虽然通过迭代的方法可以进行测量点的优化选取,但这种过程也是非常耗时的.最近wu等提出了Helmhohz方程最小二乘法的声全息技术,该技术通过将辐射声压在一系列基函数上展开,展开系数由假设解和实际测量的声压相匹配确定,采用最小二乘法来最小化此过程中所引入的误差,并将此方法应用于重建内,外声辐射问题”.同时,还有许多的研究人员在如何减小重建误差提高全息重建的精度方面做出了很多研究.但在过去的声全息变换技术研究中,全息测量的背景噪声干扰问题一直没有得到很好的解决,由此造成的结果是全息测量只能消声室环境下或将所有的背景干扰移除之后进行,如此就大大限制了其在工程测量中的推广应用.本文从近场声全息原理出发,利用声波不同方向传播的特点,针对平面,柱面全息测量面建立波数域内的声场分离公式.利用声场分离技术分离后的全息声压,重建源面上和声场中的声学参量得到了令人满意的结果.由于分离声场技术是在全息面上进行的,它可以应用于全息面是平面或柱面的各种类型全息技术中.2.近场声全息原理在开始讨论声场分离技术之前,首先简单地回国家自然科学基金(批准号:50275044)及教育部科学技术研究重点项目(批准号:02126)资助的课题E-mail:yufei6666sina.eom.ca790物理54卷顾一下平面和柱面声全息技术的原理,由此来定义相关的符号和引出全息重建的基本公式.2.1.平面声全息技术由理想流体媒质中小振幅声波的波动方程,可以得到不依赖于时间的时谐声波场的Helmhohz方程(,Y,z)+P(,Y,z)=0,(1)式中P(,Y,z)为空间点的复声压,是直角坐标,Y,z的函数;=叫/c=2rr/2L为声波数,c为声速,为声波长.对于z>0的空间为自由声场的情况,即所有声源均位于z=0平面以下,由格林公式可以得到方程(1)的解,即任意平面z(为大于零的常数)上的声压同边界平面=0上的声压,粒子法向振速在波数域内的关系为P(,z)=P(,)exp(i,z),(2)P(,z)=P0ckV(,)exp(ikz)/,(3)式中lD.为声介质的密度,P(,z)和P(,)为面z和面z=0上声压的二维连续Fourier变换,也就是波数域内的声压,(,)为在源平面z=0上波数域内的粒子法向振速.(2),(3)式中的和分别对应直角坐标和Y的空间频率(也称作波数),而式中的:为当+时,=一(:+).(4)当:+>时,=i(+:)一.(5)由(2),(3)式,对于任意的两平面z=z和z=z.(z>z.>0)可以建立更一般的关系:P(,z)=P(,zs)Xexpik(z一z.),(6)P(,zH):lDockV(,zs)Xexpik(z一zs)/:.(7)由(6),(7)式可得,若已知z=z.平面的声压或微粒法向振速便可以预测出更远处z=z的声压情况,再通过Euler公式可以很方便的推出z=z的微粒法向振速,还可以进一步导出二阶声场参量声强及远场的指向性等;同时可得近场声全息重建的基本公式P(,z)=P(,z)exp一ik(z一z.),(8)(,s)=P(,z)exp一ik(z一zs)/p.ck.(9)若已知z=z(全息面)的声压数据也可以反演更近表面z=z.(重建面)的声压,粒子法向振速等参量.2.2.柱面声全息技术对于方程(1),令=reos0,Y=rsin0,可以将直角坐标系转换为柱面坐标系.在柱面坐标下采用分离变量法求解方程(1),设方程解的形式为P(r,0,)=P,(r)p(O)p(z),(10)将(3)式代人柱面坐标下的Helmhohz方程可以得到1dP1dp,1dP+孥+z:0,(11)Pu由于含有P的项仅与坐标z有关,故可设该项等于常数一,结果如下所示:+:0.(12)dDz.一/同理,也可以将仅与角度0有关的项表示为+z:0,(13)PdO.若将(12)和(13)式代人(11)式,(11)式可以重新记为+÷+一):0P.(4)+了+Ir一7?(14当时,后=一,(15)当>时,=i:一.(16)通过以上分析得到的(12)和(13)式均为二阶的常微分方程,在取值为(15)式时,(14)式为Bessel方程,在取值为(16)式时,(14)式为修正的Bessel方程.由方程(12)一(14)可以解出方程(11)的通解为p(删量e瑚D:)日(n1)()e+D(n2()日(n2(r)ed:,(17)式中日,日(n2分别表示n阶的第一类和第二类Hankel函数,并分别对应向外发散传播的声波和向内收敛的声波.若将Fourier级数也看作为广义的Fourier变换,2期于飞等:声场分离技术及其在近场声全息中的应用791那么D(:)H?(r)+D(:)H(r)可看作是P(r,0,)的二维Fourier变换.对于外辐射问题,则柱面r=r上的声压和粒子径向振速同柱面r=r(r<r)上声压之间的关系可以表示为吣,rs)=吣),)z,r去豁吣(19)式中上标”表示求导数.(18)和(19)式便为柱面声全息重建外辐射问题的基本公式.对于内辐射问题,只需将(18)和(19)式中的上标(1)换成(2)便可以用来重建内辐射问题了.3.声场分离技术原理3.1.平面分离声场3.1.1.双平面全息实现方法图1表明全息面与声源面之间的空间位置关系:共有两个声源面和两个全息面,其中全息面1与坐标平面z=0重合,其余三个面均与该面平行.两个声源面上的声源可以是真正的平面声源或是点声源,也可以是虚构的声源.,D,:用来表示各平面之间的距离.声源面1全息面1全息面2声源面2声源1图1双平面分离声场方法的测量示意图全息面1上的任意点(.,Y.)的复声压记为P.(.,Y.);声源面1上的声源在该点引起的声压记为P.(.,Y.);声源面2上的声源在该点引起的声压记为P:.(.,Y.).因为声压为标量,所以可以得到同样,对于全息面2上点(:,y:)的复声压P:(,Y)也表示为两源面上声源在该面上产生的声压之和P2(2,Y2)=Pl2(2,Y2)+P22(2,Y2),(20)其中,P.:(:,Y:),P(:,Y:)为源面1,2上的声理论,便得到了在波数域内的双平面全息声场分离公式.对(27)式取二维Fourier逆变换,得到分离后的空间域内的复声压.对于近场声全息,可以将波数域的声压直接代人(8),(9)式来重建源面(或重建面)上的声压和粒子振速等声学参量.3.1.2.单平面全息实现方法全息面与声源面之间的空间位置关系如图2所示,共有两个声源面和一个全息面,其中全息面与坐标平面=0重合,全息面与声源面1,以及声源面2之间的距离分别为.,Z全息面上的任意点(,Y)的复声压记为P(,Y);声源面1上的声源在该点引起的声压记为P.(,Y);声源面2上的声源在该点引起的声压记为P:(,Y).因为声压为标量,所以可以得到P.(-,Y.)=P(.,Y.)+P2l(,Y).(19)点(,P(,Y)=P.(,Y)+P:(,Y).(28)Y)的粒子法向振速记为(,Y),声源面1792物理54卷声源1i啄L且J1/.,.y,/./.声A/AZ/,/./2图2单平面分离声场方法的测量示意图声源2上的声源在该点引起的微粒法向振速,(,Y)和声源面2上的声源引起的微粒法向振速,(,Y)之间的方向相反,故可以设(,Y)和,(,Y)的方向为正,则可以得到(,Y)=.(,Y)一:(,Y).(29)对(28),(29)式两边分别取二维连续Fourier变换,可以得到波数域内该全息平面上的各声压和各法向振速之间的关系P(k,k)=P.(k,k)+P:(k,k),(30)V(k,k,)=V(k,k)一(k,k).(31)波数域内的Euler公式,可以通过对时间域内的Euler公式作变换得到,即要先对时间作Fourier变换得到频域内的Euler公式,再对它作空间平面(x,Y)的二维连续Fourier变换.该公式也可以由(6),(7)式联立求解得到,即可以在波数域内建立单一侧声源在全息面上引起的声压和粒子法向振速之间关系.(k,k,)=P(k,k):/p0ck,(32)(k,k)=P:(k,k):/p0ck.(33)将(33),(32)式分别代入(30),(31)式,可以解得源面1上的声源在全息面上单独产生的声压为1P(k,k,)=寺(P(k,k,)+V(k,k)l0.ck/k:,(34)Pll(k:,r1)=该式为波数域内的单平面全息声场分离公式3.2.柱全息面分离声场原理3.2.1.双柱面全息实现方法图3表明了双柱面全息分离声场方法中两全息面同两个声源之间的关系,其中全息面1和2的半径分别为r.,r:.声源sl在全息面1内,声源s2在全息面2外,同时坐标轴垂直面(,Y)向外.与双平面全息分离声场方法的原理类似,全息面1上总的复声压为P,(r,.),可以表示为Pl(rl,1)=Pll(rl,1)+P2l(rl,1),(35)式中由声源Sl引起的声压为P.(r,),声源S2引起的声压为P:(r,0).面1图3双柱面分离声场方法的测量示意图对(35)式两端做二维Fourier变换,可以得到Pl(k,r1):Pll(k:,r1)+P2l.(k:,r1),(36)同理对全息柱面2,也可以得到P2(k,r2)=Pl2(k:,r1)+P22(k:,r1).(37)由柱面声全息原理中的不同全息面之间声波在波数域之间的关系,可以得到P1ln(r1)=r2),(38)r1)=r2).(39)联立(36)一(39)式可以解得声源S1在全息面1上产生的声压为!,r.)”(krr.)H(2(r:)一P:(k:,r:)(krr.)H(r.)“(krr)(krr:)一日(j_2)J_该式为双柱面全息声场分离公式.3.2.2.单柱面全息实现方法图4表明了单柱面全息分离声场方法中全息面(40)同两个声源之间的关系,全息面的半径为r.由声压标量和粒子振速矢量的叠加原理,可以得知全息面上的波数域声压和径向振速可以表示成分别由声源2期于飞等:声场分离技术及其在近场声全息中的应用793S1和S2引起的声压和振速之和声源2面图4单柱面分离声场方法的测量示意图P.(k:,r)=P(k:,r)=R.(k;,r)+P(k:,r),(41)I,n(k;,r)=V.(k:,r)一Vz(k,r),(42)由波数域的Euler方程,同一柱面上的声压同径向振速之间的关系为r)=r),(43)r)=啪:r).(44)将(41)一(44)式联立可以解得声源S1在全息面上产生的声压为P(:,r)日(,r)日(r)+ip.ckVo(k:,r)日”(krr)1t?(krr)/k,该式为单柱面全息声场分离公式4.数值计算与讨论1t?(krr)1t?¨(,r)+日(,r)日(r)前面的理论推导,均是建立在无限大全息面上的连续Fourier变换基础上的,在理论上和近场声全息的理论推导一样是严谨的.但是在实际的测量中只能在有限大小测量面上的离散点上进行,因此声场分离技术的原理还需进一步的检验和验证.下面通过在对较为典型的两个数值模拟来说明其可行性和有效性.仿真1对脉动球振动的声场模拟,用来验证双柱面全息声场分离技术.在图3所示的坐标系中,在点(0.00,0.15,一0.30)处放置一个半径为0.03m的脉动球声源1,球的振动频率为400Hz,振动幅值为1m/s,在点(0.00,一0.40,0.30)处放置同样的一个脉动球声源2.全息面1的半径r.为0.25m,全息面2的半径r,为0.30m.全息面的轴向长度为2m,采样间隔为0.05m,周向的采样间隔为9.在全息面1和2上模拟测量声压数据,利用提出的双柱面分离声场公式(40),分离出声源1单独作用在全息面1上的声压.将全息柱面上一条母线(:0,Y=一0.25)上的分离后声压同理论值,模拟测量值进行比较,结果如图5所示.从该图可以看出:无论是实部或虚部,分离后的声压同声源1单独作用产生的理论声压值都非常接近,而直接模拟测量的声压值同理论值之间的差距比较大.(45)氇崮救图5分离后的声压同理论值,测量值之间的比较+分离后的声压值,声源1上声源单独作用时的理论值,一一一在全息面上的模拟测量值利用柱面声全息技术,由分离后声压重建更靠近声源1的柱面上的声压,其幅值分布柱面图展开后如图6(a)所示,其中重建面的半径为0.20m.声源1在该重建面上产生的理论声压幅值分布图,展开后如图6(b)所示,而直接采用全息面1上的模拟测量声压来重建该面上的声压幅值图,展开后如图6(c)所示.很容易看出:通过空间声压分离后重建的794物理54卷结果图6(a),去除了背景噪声的影响,和该面上的理论幅值分布(图6(b)符合得非常好,可以很精确地判断出声源的位置和声压幅值的大小,而直接采用全息面1上的模拟测量声压重建的结果(图6(c)因受到背景干扰的影响,在点(0.00,一0.20,0.30)附近重建出一个虚假声源,在展开图上对应的点也就是坐标z=0.30,=一/2的位置.这样由于背景噪声的影响,就会造成对声源的错误识别,过去解决此种影响的办法是将背景噪声移去,或是将被测目标放置在消声室内,然后再进行测量和全息重建声场以识别噪声源.这样就严重约束了全息技术在工程实践中的应用,因为在实际的工业测量中移去所有的声源是不很容易实现的,而利用消声室测量的办法,费用是非常昂贵的.为了更好地比较声压分离后重建声场的效果,在图6的三个子图中抽取最典型的一行作比较,抽取角度=一/2的行,比较结果如图7所示.从该图中可以看出,在点z=一0.3附近,分离声场前,后重建的结果相差不大,和理论值非常接近;在坐标Z大于一0.15的重建点上,分离前,后重建的结果具有显着的差异,而分离后重建的结果和理论值比较符合,误差也非常小.仿真2简支板(3,3)阶模态振动产生的声场模拟,用来验证单平面全息声场分离技术.在图2所示的坐标系中,在声源面1上布置一个长和宽都为1m的简支撑薄板声源,该薄板的中心在z轴上,板四周为无限大障板.设该板的(3,3)模态的共振频率为200Hz,振速的最大值为1m/s.声源面1与全息面之间的距离=0.1m,全息面的尺寸为2m×2m,采样间隔在和Y轴方向上均为5cm.在源面2上放置同样的薄板声源,仅振速的最大值为0.5m/s其余条件与源面1上的板源均相同.声源面2与全息面之间的距离z=0.1m.两板源的振动方向相反,并且相位相同,同以(3,3)模态振动.在全息面上测量全息数据,分别直接利用该全息数据重建源面1的法向振速和采用分离声场技术分离的声压数据重建法向振速,并将重建结果同理论值相比较,重建面上各点的振速重建误差如图8所示.图8(a)由于所采用的全息数据为两声源面上的声源引起的声压之和,故重建结果与理论值之差便为源面2上的振速分布情况.图8(b)为采用提出的分离声场技术分离结果重建的误差分布,从该图可以发现较大的误差主要集中在板声源的边界处,6040馨20n016040馨爱20016040j罾馨2001图6声场分离前,后重建声压幅值同理论值的比较(a)声场分离后重建的声压幅值图;(b)声源1单独作用的理论声压幅值图;(c)声场分离前重建的声压幅值图这是由于此处的振速的导数不连续造成的.如果定义声全息重建总的误差E.为,(46)式中Ei为重建面上各点的重建误差;i为重建面上各点法向振速的理论值.利用这种定义的总误差,可以计算出直接利用测量数据重建的总误差E=49.96ff,利用分离后2期于飞等:声场分离技术及其在近场声垒息中的应用795馨出丑圈7分离静后的重建声压幅值同理论值的单行比较+分离后的声压重建结果,亩源l单独作用的理论值,<圳制氍i丹离前的声压重宦结果圈8分离声压前,看重建法向振速误羲廿布图(a)分离声压前重建结果;(b)分离声压后重建结果的声压数据重建的总误差E.=6.36%其中,直接利用测量数据重建的总误差主要是由声源面2上的干扰声源引起的.而利用测量分离后的数据重建总误差主要是来自于NAH过程带来的误差.5.实验验证全息测量实验在半消声室中进行,实验室本底噪声小于16dBA,实验装置在消声室中的实物照片如图9所示两个全息面的大小均为1.0mx0.7m的平面,沿水平方(轴)和竖直方向(Y轴)的全息测量间隔均为5cm.目标声源的音箱纸盆中心坐标为f一0.15m,0.Om,一0.125m).干扰音箱的纸盆中心坐标位于(0,15m,0.0m,0.225m).全息面1与坐标平面重合,其中心点为坐标原点,与垒息面2之问的距离D为7.5cm.围9双平全息面分离声场技术实验装置宴钫圈利用分离后的全息声压重建娠面上的有功声强,结果如图10(b)所示.和直接利用全息面l上的带有背景干扰的声压数据重建l碌面上的有功声强(如图10(a)所示)相比:图10(b)很好地去除了背景源对全息重建结果的影响.同时,还可以从图10(a)发现,较为发散的伪声源位置与干扰声源的位置是相对应的.最后.将背景声源移开后再测量全息面1上的复声压,重建目标声源面1上的有功声强结果如图10(c)所示.比较图10(b)与(c)可以知道,采用分离声场技术后的全息声压重建声源结果几乎可以达到直接移去背景声潭后重建的效果.关于实验参数选取和误差敏感性方面的详细研究将另文介绍.通过以上两个仿真算例和一个实验,考核了空间声场分离技术在离散化后的可行性和有救性,及其在近场声全息中的应用.以上两例全息面上的声压均为直接测量得到的声压值,当然也可以采用BAHIM中通过测量二维切向声强的方式来获得全息面上的复声压.由于分离声场技术可以分离来自不同方向的声波,因此,在以平面全息面或柱面全796物理54卷翟.圣1翟.x/m专.图1O采用不同方法重建的目标源面有功声强分布比较(a)未利用分离技术重建的结果;(b)采用分离技术后重建的结果;(c)移去背景干扰后直接重建的结果息面作为全息测量面的BEM或HELS声全息技术中,关于背景噪声的处理问题也可以采用本文提出的声场分离技术来实现.6.结论通过理论分析,利用声压标量和粒子振速矢量的叠加原理以及波场的外推和波数域的Euler公式,建立了空间声场分离公式.利用该公式可以分离出全息面两侧声源各自在全息面上引起的声压,进而通过分离后的声压来重建声场中的声学参量,就能够排除另一侧声源带来的干扰.针对理论推导的公式,进行了仿真模拟和实验验证,结果表明:对在离散化和有限尺寸的全息面上测得的数据来说,声场分离技术依然是可行的,有效的,即便用于较复杂的声场的分析;并和未采用该技术而产生的伪声源情况做出了对比.本文提出了平,柱全息面的声场分离技术,将NAH和BANIM的应用范围进一步拓宽:在声全息测量的过程中没有必要再去移除全息面某一侧的声源,这就进一步降低对工业测量的要求,对声全息技术进一步地推向工程实践具有非常积极的意义.11JYouYX,MiaoGP2002ActaPhys.Sin.512038(inChinese)6尤云祥,缪国平物理2002512038【2JYouYXetal2001ActaPhys.Sin.501103(inChinese)尤云7祥等物理2001501103l3JYouYX,MiaoGP2002ActaPhys.Sin.51270(inChinese)8尤云祥,缪国平2002物理512709l4JWilliamsEGandMaynardJD1980Phys.Rev.Letts.455541Ol5JWilliamsEGeta/1980J.Acoust.Soc.Am.68340MaynardJD,WilliamsEGandLeeY1985J.Acoust.Soc.Am.781395VeronesiWAandMaynardJD1987J.Acoust.Soc.Am.8113o7LoyauTandPascalJC1988J.Acoust.Soc.Am.841744MannIIIJAandPascalJC1992NoisCon.Eng.J.393VeronesiWAandMaynardJD1989J.Acoust.Soc.Am.8513o72期于飞等:声场分离技术及其在近场声全息中的应用79711GardnerBK,BerharRJ1988Trans.ASME,J.Vib.Aoust舶flab.D.11084KimGT,LeeBT199OJBaiMR1992J.Acoust.KimBK,IhJG1996J.KimBK.IhJG2000J.WangZ,WuSF1997J.Sound.V/b.136245SOc.Am.92533Acoust.SOc.Am.1003003Acoust.Soc.Am.1o73289Acoust.Soc.Am.1022020202122WuSF2000J.Acoust.Soc.Am.1072511RayessNE,WuSF2000J.Acoust.SOc.Am.1072955WuSF,RayessNE,ZhaoX2001J.Acoust.SOc.Am.1092771NamKU,KimYH1999J.Acoust.SOc.Am.1061655SaijyouK,YoshikawaS2001J.Acoust.SOc.Am.1102007KwonHS,KimYH1995J.Acot.SOc.Am.982104Soundfieldseparationtechniqueanditsapplicationsinnear?fieldacousticholographyYuFeiChenJianLiWei?-BingChenXin?-Zhao(DynamicMeuringCenter,HefeiUniversityofTechnology,Hefei230009,Ch/na)(Received4August2003;reusedmanuscriptreceived24May2004)AbstractSoundfieldseparationtechniquehasbeenproposed,whichovercomesthelimitationonapplicationsofnear-fieldacousticholography(NAH).Thelimitationisthatthesoundfieldononesideofholographicplanemustbefree,thatistosay,allthesoundsourcesmustbeconfinedtotheotherside.Soundfieldseparationtechniquewithdoubleholographicplaneshasbeenbuiltaccordingtothewavefieldextrapolationtheoremandthesuperpositiontheoremofsoundpressurescalarinwavenumberdomain,andthetechniquewithsingleholographicplaneshasalsobeenbuiltaccordingtoEulerSformulationandthevectorsuperpositiontheoremofparticlevelocityinwavenumberdomain.Aparticularfeatureofthetechniqueisthattheinfluenceofbackgroundnoiseonreconstructedresultscanberemovedwhileretainingthebackgroundnoiseintheholographicmeasurements.Thederivationoftheprincipleverifiesthetechniquetheoretically,andthecomputersimulationsandexperimentdemonstrateitsfeasibilityandeffectiveness.Keywords:acousticholography,wavenumberdomain,soundfieldseparation,backgroundnoisePACC:4320,4363ProjectsupportedbyNationalNaturalScienceFoundationofChina(GrantNo.50275044)andbyKeyProgramofChineseMinistryofEducation(GrantNo.02126).