【声传播】——球面波的反射

概述

本文整理了球面波的反射问题，核心思路使用角谱法将球面波展开成平面波，利用已知的平面波反射系数，合成反射的球面波。
本文在后面讨论了反射系数与入射角无关的情况下球面反射波的形式以及远程条件下的球面反射波形式。

角谱法展开球面波

利用二维傅里叶变换对，可以将球面波展开成平面波的组合，即角谱法。
注：这里用$\frac{e^{jkR}}{R}$的形式来表示球面波，因为水声里一般用$e^{j(kx-wt)}$来表示平面波。
$$\frac{e^{jkR}}{R}=\iint A(k_x,k_y)e^{j(k_{x}x+k_{y}y+k_{z}z)}\mathrm{d}{\mathrm{k}}_{\mathrm{x}}\mathrm{d}{\mathrm{k}}_{\mathrm{y}}$$
为了求$A(k_x,k_y)$，取$z=0$平面，此时$R=r=\sqrt{x^2+y^2}$，则可将上式化简得到：
$$\frac{e^{jkr}}{r}=\iint A(k_x,k_y)e^{j(k_{x}x+k_{y}y)}\mathrm{d}{\mathrm{k}}_{\mathrm{x}}\mathrm{d}{\mathrm{k}}_{\mathrm{y}}$$
则$A(k_x,k_y)$可以表示为：
$$(2\pi {)}^{2}A(k_x,k_y)=\iint \frac{e^{jkr}}{r}{e}^{-j(k_{x}x+k_{y}y)}\mathrm{d}\mathrm{x}\mathrm{d}\mathrm{y}$$
利用坐标变换式，将直角坐标变换成极坐标形式，波矢域与空间域都变成极坐标形式，积分变量在替换时需要乘上一个雅克比行列式，下图是转换关系：

这么变换的目的就是把上式右边的积分给算出来
$$(2\pi {)}^{2}A(k_x,k_y)={\int }_0^{2\pi }\mathrm{d}\phi {\int }_0^{\infty }{e}^{jr[k-\xi cos(\psi -\phi )]}\mathrm{d}\mathrm{r}$$
先对$r$积分，再对$\phi$做积分，通过查积分表，可以得到结果：
$$A(k_x,k_y)=\frac{j}{2\pi \sqrt{k^2-k_x^2-k_y^2}}=\frac{j}{2\pi k_z}$$
综上，球面波展开成平面波的组合：
$$\frac{e^{jkR}}{R}=\iint \frac{j}{2\pi k_z}{e}^{j(k_{x}x+k_{y}y\pm k_{z}z)}\mathrm{d}{\mathrm{k}}_{\mathrm{x}}\mathrm{d}{\mathrm{k}}_{\mathrm{y}}$$
这里的平面可以是普通的平面波，也可以是非均匀平面波（凋落波），即由于$\xi >k$，$k_z$为纯虚数，如下式，但该波仅在$z>0$时才有物理意义
$$e^{-k_{z}z}{e}^{j(k_{x}x+k_{y}y)}$$

根据平面波的反射系数求解球面反射波的表达式

示意图如下，注意这里实际上是一个二维的场景，已将$x,y$缩为$r$，取的是三维空间中的一个径向的面来分析问题

根据上面得到的结论，可以写出球面反射波的积分形式：
$$p_r=\iint \frac{jV(k_z)}{2\pi k_z}{e}^{j(k_{x}x+k_{y}y+k_z(z+z_0))}\mathrm{d}{\mathrm{k}}_{\mathrm{x}}\mathrm{d}{\mathrm{k}}_{\mathrm{y}}$$
将波矢域的直角坐标形式转换成极坐标形式，同样是为了求出上述积分式：
$$p_r={\int }_0^{\infty }\frac{jV(k_z)}{2\pi k_z}{e}^{j(k_z(z+z_0))}\xi \mathrm{d}\xi {\int }_0^{2\pi }{e}^{j\xi cos(\psi -\phi )}\mathrm{d}\psi$$
后面那个对$\psi$的积分的结果实际上是0阶贝塞尔函数，0阶贝塞尔函数又可以转换为第一类0阶汉克尔函数和第二类0阶汉克尔函数的组合，化简得：
$$p_r=\frac{j}{2}{\int }_{-\infty }^{\infty }\frac{V(k_z)}{k_z}{H}_0^{(1)}(\xi r)e^{j[k_z(z+z_0)]}\xi \mathrm{d}\xi$$

反射系数与入射角无关时的球面反射波表达式

如果反射系数与入射角无关时(即分界面为刚性分界面)，则可以将反射系数$V(k_z)$移到积分号外面，得到：
$$p_r=\frac{jV(k_z)}{2}{\int }_{-\infty }^{\infty }{H}_0^{(1)}(\xi r)\frac{e^{j[k_z(z+z_0)]}}{k_z}\xi \mathrm{d}\xi$$
可以发现，这个积分就是柱面波与球面波之间的变换式，故上式可以转换为球面波形式，其中$R_1=\sqrt{r^2+(z+z_0{)}^2}$：
$$p_r=\frac{V{e}^{jk{R}_1}}{R_1}$$
上式的物理意义就是$O^{{}^{\prime }}$处虚源产生的声场，点源强度乘上了一个反射系数。

一般情况下的远场球面反射波表达式

远场条件$\xi r\gg 1$，此时对0阶第一类汉克尔函数作近似：
$$H_0^{(1)}(\xi r)\approx \sqrt{\frac{2}{\pi \xi r}}{e}^{j(\xi r-\frac{\pi }{4})}$$
代入前面的积分式，可以得到：
$$p_r=\frac{e^{\frac{j\pi }{4}}}{\sqrt{2\pi r}}{\int }_{-\infty }^{\infty }\frac{\sqrt{\xi }V(k_z)}{k_z}{e}^{jw(\xi )}\mathrm{d}\xi$$
$$w(\xi )\equiv \xi r+k_z(z+z_0)$$
因为$w(\xi )$变化剧烈，而$\frac{\sqrt{\xi }V(k_z)}{k_z}$变化相对来说缓慢，故可以用稳相法求积分，即先求出稳相点，然后在稳相点附近展开到二阶项，然后代入积分项，直流项可直接提出来，二阶项用到了换元技巧以及解析函数的围线积分为0的性质，对积分路径进行了变换，最终化简得到了远程条件下的球面反射波的表达式：
$$p_r=\frac{V({\theta }_0)e^{jk{R}_1}}{R_1}$$
其中${\theta }_0$对应的是图上的角度，其实也是稳相点求出的角度，恰好是声源以声线入射再反射到达目标点对应的角度。

Reference：
L.M. Brekhovskikh,&Yu.P. Lysanov. Fundamentals of Ocean Acoustic(Chapter 4)