可控源音频大地电磁理论基础

可控源音频大地电磁

理论基础

1、控制方程

电磁场满足麦克斯韦方程组：
$$\begin{array}{c}\begin{array}{rlrl}& \nabla \times \mathbf{H}=\sigma \mathbf{E}+\frac{\partial \mathbf{D}}{\partial t}& & (安培定律)\\ & \nabla \times \mathbf{E}=-\frac{\partial \mathbf{B}}{\partial t}& & (法拉第定律)\\ & \nabla \cdot \mathbf{B}=0& & (磁通量连续定律)\\ & \nabla \cdot \mathbf{D}=\mathbf{q}& & (库伦定律)\end{array}\end{array}$$
其中$\mathbf{E}(V/m)$为电场强度，$\mathbf{H}(A/m)$为磁场强度，$\mathbf{B}(Wb/m^2)$为磁感应强度，$\mathbf{D}(C/m^2)$为电位移矢量，$\mathbf{q}(C/m^3)$为电荷密度，这四个物理量是电磁场的基本量，它们通过介质电性参数$ϵ,\mu$进行联系，在各向同性介质中，介质电性参数为标量。本构方程如下：
$$\begin{array}{c}\begin{array}{r}\mathbf{D}=ϵ\mathbf{E}\\ \mathbf{H}=\frac{\mathbf{B}}{\mu }\\ \mathbf{J}=\sigma \mathbf{E}\end{array}\end{array}$$其中，$\mathbf{J}$为电流密度，$\sigma$为介质导电率，$ϵ$为介质介电常数，$\mu$为磁导率。
对安培定律和法拉第定律两边取旋度，并且考虑在场源外（$\nabla \cdot \mathbf{D}=0$），有：
$$\begin{array}{c}\begin{array}{rl}\nabla \times \nabla \times \mathbf{H}& =\sigma \nabla \times \mathbf{E}+\frac{\partial }{\partial t}\nabla \times \mathbf{D}\\ \nabla (\nabla \cdot \mathbf{H})-{\nabla }^2\mathbf{H}& =\sigma (-\frac{\partial B}{\partial t})+\frac{\partial }{\partial t}ϵ(-\frac{\partial B}{\partial t})\\ -{\nabla }^2\mathbf{H}& =ϵ\mu \frac{{\partial }^2\mathbf{H}}{\partial t^2}+\sigma \mu \frac{\partial \mathbf{H}}{\partial t}\end{array}\end{array}$$
$$\begin{array}{c}\begin{array}{rl}\nabla \times \nabla \times \mathbf{E}& =\frac{\partial }{\partial t}\nabla \times \mathbf{B}\\ \nabla (\nabla \cdot \mathbf{E})-{\nabla }^2\mathbf{E}& =\frac{\partial }{\partial t}\mu (\sigma \mathbf{E}+\frac{ϵ\partial \mathbf{E}}{\partial t})\\ -{\nabla }^2\mathbf{E}& =ϵ\mu \frac{{\partial }^2\mathbf{E}}{\partial t^2}+\sigma \mu \frac{\partial \mathbf{E}}{\partial t}\end{array}\end{array}$$上述分别为$\mathbf{H}$和$\mathbf{E}$满足的微分方程，称为电报方程
在CSAMT方法中，所激发的源是谐变场，即：
$$\begin{array}{c}\mathbf{H}={\mathbf{H}}_0{e}^{-i\omega t},\mathbf{E}={\mathbf{E}}_0{e}^{-i\omega t}\end{array}$$
代入到微分方程组中，可以得到：

$$\begin{array}{c}\begin{array}{rl}\nabla \times \mathbf{H}& =\sigma \mathbf{E}-i\omega ϵ\mathbf{E}=-i\omega {ϵ}^{\ast }\mathbf{E},{ϵ}^{\ast }=ϵ+i\frac{\sigma }{\omega }\\ \nabla \times \mathbf{E}& =i\omega \mu \mathbf{H}\\ \nabla \cdot \mathbf{H}& =0\\ \nabla \cdot \mathbf{E}& =0\end{array}\end{array}$$式中${ϵ}^{\ast }=ϵ+i\frac{\sigma }{\omega }$为复介电常数
同样的在方程两端取旋度，可以得到：
$$\begin{array}{c}\begin{array}{rl}\nabla \times \nabla \times \mathbf{H}& =-i\omega {ϵ}^{\ast }\nabla \times \mathbf{E}\\ \nabla (\nabla \cdot \mathbf{H})-{\nabla }^2\mathbf{H}& =-i\omega {ϵ}^{\ast }\cdot i\omega \mu \mathbf{H}\\ -{\nabla }^2\mathbf{H}& ={\omega }^2\mu {ϵ}^{\ast }\mathbf{H}\end{array}\end{array}$$令$k^2={\omega }^2\mu {ϵ}^{\ast }={\omega }^2\mu ϵ+i\omega \mu \sigma$，$k$称为复波数，或者叫做传播系数,则可以得到：
$$\begin{array}{c}{\nabla }^2\mathbf{H}+k^2\mathbf{H}=0\end{array}$$
同理可以得到：
$$\begin{array}{c}{\nabla }^2\mathbf{E}+k^2\mathbf{E}=0\end{array}$$
上述两式被称为谐变电磁场的基本微分方程——亥姆霍兹（Helmholtz equation）齐次方程。
在求解电磁场边值问题时，使用$\mathbf{H},\mathbf{E}$是不方便的，在这里引入一个矢量位$\mathbf{A}$（磁矢量位），以此来减少求解过程中未知数，方便推导。
$$\begin{array}{c}\mathbf{H}=\nabla \times \mathbf{A}\end{array}$$代入到$(6)$式中，可以得到：
$$\begin{array}{c}\begin{array}{rl}& \nabla \times \mathbf{E}=i\omega \mu \nabla \times \mathbf{A}\\ & \nabla \times \left(\mathbf{E}-i\omega \mu \mathbf{A}\right)=0\end{array}\end{array}$$根据矢量恒等式$\nabla \times (\nabla U)=0$，可得：
$$\begin{array}{c}\begin{array}{rl}& \mathbf{E}-i\omega \mu \mathbf{A}=\nabla U\\ & \mathbf{E}=i\omega \mu \mathbf{A}+\nabla U\end{array}\end{array}$$其中$U$为电磁场的标量位，在直流稳恒电场中，满足$\mathbf{E}=-\nabla U$，再考虑$(6)$式，得：
$$\begin{array}{c}\begin{array}{rl}\nabla \times \nabla \times \mathbf{A}& =-i\omega {ϵ}^{\ast }\left(i\omega \mu \mathbf{A}-\nabla U\right)\\ \nabla \nabla \cdot \mathbf{A}-{\nabla }^2\mathbf{A}& ={\omega }^2\mu {ϵ}^{\ast }\mathbf{A}+i\omega {ϵ}^{\ast }\nabla U\\ \nabla \left(\nabla \cdot \mathbf{A}-i\omega {ϵ}^{\ast }U\right)& ={\nabla }^2\mathbf{A}+k^2\mathbf{A}\end{array}\end{array}$$令上述等式的左端等于零，可以得到以下结论：
$$\begin{array}{c}\begin{array}{rl}\nabla \left(\nabla \cdot \mathbf{A}-i\omega {ϵ}^{\ast }U\right)& =0\\ U& =\frac{\nabla \cdot \mathbf{A}}{i\omega {ϵ}^{\ast }}\\ {\nabla }^2\mathbf{A}+k^2\mathbf{A}& =0\end{array}\end{array}$$上式即为磁矢量位的亥姆霍兹方程，结合式$(11)$和$(13)$，可以得到电场强度$\mathbf{E}$的表达式：
$$\begin{array}{c}\begin{array}{rl}\mathbf{E}& =i\omega \mu \mathbf{A}+\nabla U\\ & =i\omega \mu \mathbf{A}+\nabla \frac{\nabla \cdot \mathbf{A}}{i\omega {ϵ}^{\ast }}\\ & =i\omega \mu \left(\mathbf{A}+\frac{\nabla \nabla \cdot \mathbf{A}}{i\omega {ϵ}^{\ast }\cdot i\omega \mu }\right)\\ & =i\omega \mu \left(\mathbf{A}-\frac{\nabla \nabla \cdot \mathbf{A}}{k^2}\right)\end{array}\end{array}$$至此问题简化为以下微分方程的求解：
$$\begin{array}{c}\begin{array}{rl}& {\nabla }^2\mathbf{A}+k^2\mathbf{A}=0\\ & \mathbf{H}=\nabla \times \mathbf{A}\\ & \mathbf{E}=i\omega \mu \left(\mathbf{A}-\frac{\nabla \nabla \cdot \mathbf{A}}{k^2}\right)\end{array}\end{array}$$

2、边界条件

磁矢量$\mathbf{A}$满足以下边界条件：
$$\begin{array}{c}\begin{array}{rlrl}{\mu }_1{\mathbf{A}}_{x1}& ={\mu }_2{\mathbf{A}}_{x2}& & \frac{\partial {\mathbf{A}}_{x1}}{\partial z}=\frac{\partial {\mathbf{A}}_{x2}}{\partial z}\\ {\mu }_1{\mathbf{A}}_{y1}& ={\mu }_2{\mathbf{A}}_{y2}& & \frac{\partial {\mathbf{A}}_{y1}}{\partial z}=\frac{\partial {\mathbf{A}}_{y2}}{\partial z}\\ {\mathbf{A}}_{z1}& ={\mathbf{A}}_{z2}& & \frac{{\mu }_1}{k_1^2}\nabla \cdot {\mathbf{A}}_1=\frac{{\mu }_2}{k_2^2}\nabla \cdot {\mathbf{A}}_2\end{array}\end{array}$$