CSAMT简介一

可控源音频大地电磁方法简介

1、概述

可控源音频大地电磁法（英文名称：Controlled-Sourse Audio-frequency Magnetotelluric，简称CSAMT）是一种控制接收区位于远离发射源的远场的可控源电磁测量方法。它是基于音频大地电磁法（Audio Magnetotelluric，AMT）和大地电磁测深法（Magnetotelluric，MT）的基础发展起来的，通过扫描特定的频率范围，通过后续反演处理，可以生成地面的深度剖面电阻率模型，大地电磁法的主要优点是可以具有相当显著深度穿透。
CSAMT方法测量的在测点位置上的正交且水平的电场和磁场（如$E_x,H_y$），通过它们的比值来反演地下电阻率的特性分布情况，该方法普遍用于矿产资源、油气、地热即工程勘查等许多领域之中。

2、特点(zonge.com，geophysics.inc)

1、相较于AMT和MT方法测量天然源的电场和磁场在较宽频带范围内的波动，CSAMT使用人工发射源取代自然源，这三种方法在相互比较时各有优缺点。CSAMT方法的主要优势是信号更强，因此更具有一致性。这使得处理和解释更加容易。CSAMT方法的数据采集速度比AMT和MT测量快得多，可以降低总体测量成本。
2、CSAMT方法的主要缺点与“近场”效应有关，它会扭曲数据，不易后续的反演解释。为了克服这个问题，CSAMT发射机和接收器的分离、定向和间隔要通过检查可用的TDEM或电阻率测深数据仔细确定。
3、可控源音频大地电磁法工作频率范围为 $10^{-1}\sim 10^4$Ｈｚ，在此频段可忽略位移电流的作用，地下介质和空气的磁导率相接近。
4、CSEM / CSAMT在绘制20 ~ 1000米深度范围时非常有用。垂直分辨率是深度的5%到20%，横向分辨率由测点站距决定。
5、静位移效应（现象），所谓静态效应是指当近地表存在局部导电性不均匀体时，电流流过不均匀体表面而在其上形成“积累电荷”，由此产生一个与外电流场成正比(比例系数不随频率变化）的附加电场。其影响的结果是，各电磁场分量对测深频率按出现比例出现放缩，从而使绘于双对数坐标系中的频率测深曲线，沿视电阻率轴（即纵轴）发生上下平移。当局部不均匀体为低阻体时，测深曲线向下平移;而若为高阻体，则向上平移。

3、正演方法

1、正演计算方法多使用有限元法，因为正演的问题是一个准静态问题，而且有场源分布、起伏地表的影响。使用的麦克斯韦方程组形式多为积分形式，由于场源附近的电场和磁场变化非常迅速，如何处理场源奇异性的策略通常有两种：（１）采取将电磁场总场分离成背景场和二次场的策略；（２）采用伪 δ函数来代替场源直接计算总场的策略。
2、通常，主场是计算存在解析公式的一维结构的。这种经典的主场和次场分离对于复杂结构是无效的，因为确定合适的背景结构是困难的。Mitsuhata在考虑地形的CSAMT二维数值模拟的研究中所述：当地下介质复杂的情况下，很难找到一个简单的背景场构造，而采用伪 δ函数来代替场源直接计算总场的策略不存在这个问题，因此更适合于模拟复杂地质结构。在受控源电磁数据建模中，通常将电磁场划分为主（背景）场和二次（散射）场，以避免源奇异性，只对异常物体引起的二次场进行数值计算。然而，这种传统的方案对于缺乏简单背景结构的复杂结构并无效(Yuji Mitsuhata)。

3.1直接计算总场(钟苏美等)

理论公式

在电偶极子或者水平线圈激发情况下，总电场强度矢量 $\mathbf{E}$和总磁场强度矢量 $\mathbf{H}$相互作用的关系满足麦克斯韦方程：
$$\{\begin{array}{l}\oint \mathbf{H}\cdot \mathrm{d}\mathbf{l}=\iint \left(\sigma \mathbf{E}+{\mathbf{J}}^s\right)\cdot \mathrm{d}\mathbf{S}\\ \oint \mathbf{E}\cdot \mathrm{d}\mathbf{l}=i{\mu }_0\omega \iint \mathbf{H}\cdot \mathrm{d}\mathbf{S}\end{array}$$式中：ω表示角频率；${\mu }_0$为真空磁导率；σ为介质的电导率；Ｅ表示电场强度；Ｈ表示磁场强度；${\mathbf{J}}^s$为源电流密度。

场源奇异性处理(Yuji Mitsuhata)

使用伪 δ函数代替场源直接模拟总场以解决场源的奇异性，伪$\delta$函数表达式如下：
$$\begin{array}{c}\delta (x-x_0)=\frac{1}{2\tau }\{\begin{array}{ll}0& (x-x_0\le -2\tau )\\ \frac{(x-x_0+2\tau {)}^2}{2{\tau }^2}& -2\tau <(x-x_0)\le -\tau \\ -\frac{(x-x_0+2\tau {)}^2}{2{\tau }^2}+\frac{2(x-x_0+2\tau )}{{\tau }^2}-1& -\tau <(x-x_0)\le \tau \\ \frac{(x-x_0+2\tau {)}^2}{2{\tau }^2}-\frac{4(x-x_0+2\tau )}{{\tau }^2}+8& \tau <(x-x_0)\le 2\tau \\ 0& 2\tau <(x-x_0)\end{array}\end{array}$$
)

为了可以直接求解总电磁场，虚将原有的研究边界处二次场值为零的边界条件修改为均匀半空间中电偶极子在远区的一次场值的边界条件，并通过公式推导，获得计算均匀半空间中电偶极子在远区一次场值的解析表达式。

3.2 分离一次场和二次场（邓居智等）

该种方法可以降低直接计算电磁场总场的难度，将总场分解为背景场（一次场）和感应场（二次场），背景场（一次场）利用快速汉克尔变换求取，感应场（二次场）采用数值计算求解。将总场表示为一次场和二次场的和：
$$\{\begin{array}{l}\mathbf{E}={\mathbf{E}}^a+{\mathbf{E}}^b\\ \mathbf{H}={\mathbf{H}}^a+{\mathbf{H}}^b\end{array}$$式中：分别$\mathbf{E}、\mathbf{H}$为总场的电场强度和磁场强度，${\mathbf{E}}^b$和${\mathbf{E}}^a$分别为背景场和二次场的电场强度，${\mathbf{H}}^b$和${\mathbf{H}}^a$分别为背景场和二次场的磁场强度。
经过总场分离后，背景场为均匀半空间（水平层状介质）有限长电偶源所激发的电磁场，其满足的麦克斯韦方程积分形式为：
$$\{\begin{array}{l}\oint {\mathbf{H}}^b\cdot \mathrm{d}\mathbf{l}=\iint \left({\sigma }^b{\mathbf{E}}^b+{\mathbf{J}}^e\right)\cdot \mathrm{d}\mathbf{S}\\ \oint {\mathbf{E}}^b\cdot \mathrm{d}\mathbf{l}=i{\mu }_0\omega \iint {\mathbf{H}}^b\cdot \mathrm{d}\mathbf{S}\end{array}$$将总场所满足的麦克斯韦方程减去背景场所满足的麦克斯韦方程，即得到二次场满足麦克斯韦方程的积分形式为：
$$\{\begin{array}{l}\oint {\mathbf{H}}^a\cdot \mathrm{d}\mathbf{l}=\iint \mathbf{J}\cdot \mathrm{d}\mathbf{S}\\ \oint {\mathbf{E}}^a\cdot \mathrm{d}\mathbf{l}=i{\mu }_0\omega \iint {\mathbf{H}}^a\cdot \mathrm{d}\mathbf{S}\end{array}$$式中： $\mathbf{J}$为电流密度，它与背景电场二次电场和电导率的关系为
$$\mathbf{J}=\sigma {\mathbf{E}}^a+\Delta \sigma {\mathbf{E}}^b$$其中 $\Delta \sigma$为剩余电导率：
$$\Delta \sigma =\sigma -{\sigma }^b$$在笛卡儿右手坐标系中，用交错采样剖分网格在研究区域内对剩余电阻率和二次场满足的麦克斯韦积分方程进行离散化，可获得关于地下各网格单元采样点处二次磁场的正演方程:
$$\begin{array}{c}{\mathbf{KH}}_2=\mathbf{s}\end{array}$$其中，$\mathbf{K}$为对称的大型稀疏系数矩阵，${\mathbf{H}}_2$为待求解各网格单元采样点处的二次磁场三分量组成的列向量；$\mathbf{s}$为与一次场及边界场值有关的列向量。求解线性方程组从而获得二次磁场值列向量${\mathbf{H}}_2$。
经过上述变换后，CSAMT 三维数值模拟问题即转化为背景场和二次场的求解，背景场可以通过快速汉克尔变换求取，二次场采用三维交错采样有限差分法进行数值计算，将求解得到的背景场值加上二次场值即为有限长电偶源激发下CSAMT三维电磁场分布。

边界条件

1、求解二次磁场时可以采用边界二次磁场值为零的边界条件，
2、直接求总场使用的边界上的总磁场是均匀半空间中电偶极子在远区产生的一次磁场响应。

4、控制方程及边界条件

求解的微分方程如下： $$\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 \left(\mathbf{A}-\frac{1}{k^2}\nabla \nabla \cdot \mathbf{A}\right)\end{array}\end{array}$$其中$\mathbf{A}$为磁矢势，$\mathbf{E}$为电场强度，$k$为波数
磁矢量$\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}$$

5、磁场与磁矢量势的关系

给出Biot-Savart定律： $$\begin{array}{c}\mathbf{B}(p)=\frac{{\mu }_0}{4\pi }{\int }_V\frac{\mathbf{j}(q)\times \mathbf{L}qp}{L^{3}qp}dV\end{array}$$ 利用方向导数（梯度）作以下变换：
$$\begin{array}{c}\frac{\mathbf{L}qp}{L^{3}qp}={\nabla }^q\frac{1}{L_{qp}}=-{\nabla }^p\frac{1}{L_{qp}}\end{array}$$其中上标字母表示梯度的方向， 利用叉乘的交换性质，可以得到： $$\begin{array}{c}\mathbf{B}(p)=\frac{{\mu }_0}{4\pi }{\int }_V\mathbf{j}(q)\times {\nabla }^q\frac{1}{L_{qp}}dV=\frac{{\mu }_0}{4\pi }{\int }_V{\nabla }^p\frac{1}{L_{qp}}\times \mathbf{j}(q)dV\end{array}$$ 再根据下列的矢量计算恒等式：
$$\begin{array}{c}\nabla \times \left(\mathbf{ab}\right)=\nabla \mathbf{a}\times \mathbf{b}+\mathbf{a}\nabla \times \mathbf{b}\end{array}$$则
$$\begin{array}{c}\begin{array}{rl}& {\nabla }^p\times \frac{\mathbf{j}}{L_{qp}}={\nabla }^p\frac{1}{L_{qp}}\times \mathbf{j}+\frac{{\nabla }^p\times \mathbf{j}}{L_{qp}}\\ & {\nabla }^p\frac{1}{L_{qp}}\times \mathbf{j}={\nabla }^p\times \frac{\mathbf{j}}{L_{qp}}-\frac{{\nabla }^p\times \mathbf{j}}{L_{qp}}\end{array}\end{array}$$代入到$(4)$式中，可得：，
$$\begin{array}{c}\begin{array}{rl}\mathbf{B}(p)& =\frac{{\mu }_0}{4\pi }{\int }_V{\nabla }^p\times \frac{\mathbf{j}}{L_{qp}}-\frac{{\nabla }^p\times \mathbf{j}}{L_{qp}}dV\\ & =\frac{{\mu }_0}{4\pi }{\int }_V{\nabla }^p\times \frac{\mathbf{j}}{L_{qp}}dV-\frac{{\mu }_0}{4\pi }{\int }_V\frac{{\nabla }^p\times \mathbf{j}}{L_{qp}}dV\end{array}\end{array}$$由于矢量电流密度$\mathbf{j}$的分布只与$q$的分布有关，则上式第二个积分项为零，则：
$$\begin{array}{c}\begin{array}{rl}\mathbf{B}(p)& =\frac{{\mu }_0}{4\pi }{\int }_V{\nabla }^p\times \frac{\mathbf{j}}{L_{qp}}dV\\ & ={\nabla }^p\times \frac{{\mu }_0}{4\pi }{\int }_V\frac{\mathbf{j}}{L_{qp}}dV\end{array}\end{array}$$又$\mathbf{B}(p)=\nabla \times \mathbf{A}$，可得：
$$\begin{array}{c}\mathbf{A}(p)=\frac{{\mu }_0}{4\pi }{\int }_V\frac{\mathbf{j}(q)}{L_{qp}}dV\end{array}$$

6、环路电流中轴线上的矢量势A和磁场B

在电流环形分布的情况下，磁矢量势有以下形式：
$$\begin{array}{c}\mathbf{A}=\frac{{\mu }_{0}I}{4\pi }{\oint }_L\frac{d\mathbf{l}}{L_{qp}}\end{array}$$
)

电流环轴线上的磁场

在轴线上每个点的距离$L_{qp}$都相等，所以有： $$\begin{array}{c}\mathbf{A}=\frac{{\mu }_{0}I}{4\pi L_{qp}}{\oint }_{L}d\mathbf{l}\end{array}$$再根据沿任意闭合回路的径向积分为零，得到矢量$\mathbf{A}$在$z$轴上的值为零。同时在$z$轴上的磁感应强度$\mathbf{B}$有对称水平分量，导致在$z$轴上，只与一个垂直分量$B_z$,根据Biot-Savart定理，可知： $$\begin{array}{c}d{B}_z=|d\mathbf{B}|\frac{a}{L_{qp}}=\frac{{\mu }_{0}I}{4\pi }\frac{dl}{L_{qp}^2}\frac{a}{L_{qp}}=\frac{{\mu }_{0}Ia}{4\pi }\frac{dl}{L_{qp}^3}\end{array}$$在闭合环路上进行积分，我们可以得到：
$$\begin{array}{c}{B}_z={\oint }_{L}d{B}_{z}dl={\oint }_L\frac{{\mu }_{0}Ia}{4\pi }\frac{dl}{L_{qp}^3}dl=\frac{{\mu }_{0}Ia}{4\pi }\frac{2\pi a}{L_{qp}^3}=\frac{{\mu }_{0}M}{2\pi {\left(a^2+z^2\right)}^{3/2}}\end{array}$$其中，$M=I\pi a^2=IS$，为平面载流线圈的磁矩。 考虑当距离$z$的值远大于环形电流的半径$a$时，上式中的分母$a$可以忽略，得到：
$$\begin{array}{c}{B}_z=\frac{{\mu }_{0}M}{2\pi z^3},当z\gg a时\end{array}$$
)

在任意点上的电流回路的磁场

在这个环形电流中，矢量势有对称性，导致与变量$\varphi$无关，且每一组对称的线元电流在$p$点产生与$d\mathbf{l}$相同方向的矢势，经过矢量合成，得到的与$x-z$平面相垂直的矢势，因此矢势$\mathbf{A}$只有$\varphi$分量$A_{\varphi }$，再利用三角形余弦定理和勾股定理计算矢径的大小，得： $$\begin{array}{c}{A}_{\varphi }=\frac{{\mu }_{0}I}{4\pi }\oint \frac{d{l}_{\varphi }}{L_{qp}}=\frac{{\mu }_{0}I}{2\pi }{\oint }_0^{\pi }\frac{a\cos \varphi d\varphi }{{\left(a^2+r^2-2ar\cos \varphi +z^2\right)}^{1/2}}\end{array}$$其中，$d{l}_{\varphi }$是$d\mathbf{l}$沿着$\varphi$方向的分量，且
$$\begin{array}{r}d{l}_{\varphi }=dl\cos \varphi =(ad\varphi )\cos \varphi =a\cos \varphi \\ d\varphi L_{qp}={\left(a^2+r^2-2ar\cos \varphi +z^2\right)}^{1/2}\end{array}$$
由于后续的推导无需使用这个积分的结果，这里不进行推到，直接给出结果： $$\begin{array}{c}{A}_{\varphi }=\frac{I{\mu }_0}{\pi k}{\left(\frac{a}{r}\right)}^{1/2}\left[\left(1-\frac{k^2}{2}\right)K-E\right]\end{array}$$其中，其中，$K$和$E$是第一类和第二类的完全椭圆积分。

7、磁偶极子的磁场

假设从载流环的中心到观测点R的距离明显大于环的半径，即 $$\begin{array}{c}R={\left(r^2+z^2\right)}^{1/2}\gg a\end{array}$$则式$(15)$可以化简为：
$$\begin{array}{c}\begin{array}{rl}{A}_{\varphi }& =\frac{{\mu }_{0}Ia}{2\pi }{\int }_0^{\pi }\frac{\cos \varphi d\varphi }{{\left(R^2-2ar\cos \varphi \right)}^{1/2}}\\ & =\frac{{\mu }_{0}Ia}{2\pi R}{\int }_0^{\pi }\frac{\cos \varphi d\varphi }{{\left[1-\left(\frac{2ar}{R^2}\right)\cos \varphi \right]}^{1/2}}\\ & \stackrel{(1-x{)}^{1/2}=1-\frac{1}{2}x}{========}\frac{{\mu }_{0}Ia}{2\pi R}{\int }_0^{\pi }\left(1+\frac{ar}{R^2}\cos \varphi \right)\cos \varphi d\varphi \\ & =\frac{{\mu }_{0}Ia}{2\pi R}\left[{\int }_0^{\pi }\cos \varphi d\varphi +\frac{ar}{R^2}{\int }_0^{\pi }{\cos }^2\varphi d\varphi \right]\\ & =\frac{{\mu }_{0}Ia}{2\pi R}\left[0+\frac{ar}{R^2}{\int }_0^{\pi }\frac{\cos 2\varphi +1}{2}d\varphi \right]\\ & =\frac{{\mu }_{0}Ia}{2\pi R}\frac{ar}{R^2}{\left(\frac{\frac{1}{2}\sin 2\varphi +\varphi }{2}\right)|}_0^{\pi }\\ & =\frac{{\mu }_{0}I{a}^{2}r}{4R^3}\\ & \stackrel{\sin \theta =\frac{r}{R},M=IS}{=======}\frac{{\mu }_{0}M}{4R^2}\sin \theta \end{array}\end{array}$$将矢势写成矢量的形式： $$\begin{array}{c}\begin{array}{r}\mathbf{A}=\frac{{\mu }_0\mathbf{M}\times \mathbf{R}}{4\pi R^3}\end{array}\end{array}$$其中，$\mathbf{M}=IS\mathbf{l}z=M\mathbf{l}z$，$M=IS$为磁矩大小，$\mathbf{l}z$为方向，又因为矢势只与$\varphi$分量有关，则
$$\begin{array}{c}{A}_R=A_{\theta }=0\end{array}$$根据$\mathbf{B}=\nabla \times \mathbf{A}$,再根据矢量微分算符$\nabla$在球坐标下的性质： $$\begin{array}{c}\nabla \times \mathbf{A}=\frac{1}{r\sin \theta }\left[\frac{\partial \left(\sin \theta A\theta \right)}{\partial \theta }\right]{\mathbf{e}}_r+\frac{1}{r}\left[\frac{1}{\sin \theta }\frac{\partial A_r}{\partial \varphi }-\frac{\partial (rA\varphi )}{\partial r}\right]{\mathbf{e}}_{\theta }+\frac{1}{r}\left[\frac{\partial (e{A}_{\theta })}{\partial r}-\frac{\partial A_r}{\partial \theta }\right]{\mathbf{e}}_{\varphi }\end{array}$$结合$(18)$式，得：
$$B_R=\frac{1}{R\sin \theta }\frac{\partial \left(\sin \theta A_{\varphi }\right)}{\partial \theta },B_{\theta }=-\frac{1}{R}\frac{\partial \left(R{A}_{\varphi }\right)}{\partial R},B_{\varphi }=0$$进而， $$\begin{array}{c}{B}_R=\frac{2{\mu }_{0}M}{4\pi R^3}\cos \theta ,B_{\theta }=\frac{{\mu }_{0}M}{4\pi R^3}\sin \theta ,B_{\varphi }=0\end{array}$$

8、非导电介质中磁偶极子的准平稳场

考虑一个磁偶极子，其力矩$M(t)$沿$z$轴定向，并且位于一个球形坐标系的原点，结合$(22)$式，
)

非导电介质中磁偶极子的准平稳场

根据法拉第电磁感应定律，有： $$\begin{array}{c}{E}_{\varphi }=-\frac{1}{2\pi r}\frac{\partial \mathbf{\Phi }}{\partial t}\end{array}$$$$\begin{array}{c}\Phi ={\int }_S\mathbf{B}\cdot \mathrm{d}\mathbf{S}=2\pi {\int }_0^{r}r{B}_z\mathrm{d}r\end{array}$$由几何关系可知， $$\begin{array}{c}\begin{array}{r}{B}_z(t)=B_R\cos \theta -B_{\theta }\sin \theta ,B_z(t)=\frac{{\mu }_{0}M(t)}{4\pi R^3}(3{\cos }^2\theta -1)\end{array}\end{array}$$代入到$(24)$式中进行积分运算得：
$$\begin{array}{c}\begin{array}{rl}& \frac{\partial \mathbf{\Phi }}{\partial t}=\frac{{\mu }_0\dot{M}(t)}{2R^3}{r}^2\&R=(r^2+z^2{)}^{1/2},\dot{M}(t)=dM/dt\\ & E_{\varphi }(t)=-{\mu }_0\frac{\dot{M}(t)}{4\pi R^2}\sin \theta \\ & B_R(t)=\frac{2{\mu }_{0}M(t)}{4\pi R^3}\cos \theta ,B_{\theta }(t)=\frac{{\mu }_{0}M(t)}{4\pi R^3}\sin \theta ,E_{\varphi }(t)=-{\mu }_0\frac{\dot{M}(t)}{4\pi R^2}\sin \theta \end{array}\end{array}$$

9、电场和磁场的复振幅的表达式

)

球坐标中的磁偶极子和矢量势

本节从${\mathbf{E}}^{\ast }=\nabla \times {\mathbf{A}}^{\ast }$出发，考虑电场只具有$\varphi$分量，偶极子为$\mathbf{M}=\mathrm{Re}{M}_0{e}^{-i\omega t}\mathbf{l}z$，在球坐标系下，结合公式(21)，可得：
$$\begin{array}{c}{E}_{\varphi }^{\ast }=\frac{1}{R}\left[\frac{\partial }{\partial R}(R{A}_{\theta }^{\ast })-\frac{\partial A^R}{\partial \theta }\right]\end{array}$$代入以下条件和(35)式，
$$\begin{array}{c}\begin{array}{rl}& A_R^{\ast }=A_z^{\ast }\cos \theta ,A_{\theta }^{\ast }=-A_z^{\ast }\sin \theta \\ & E_{\varphi }^{\ast }=\frac{C}{R^2}(1-ikR)e^{ikR}\sin \theta \\ & E_{\varphi }^{\ast }=\frac{C}{R^2}\sin \theta ,靠近偶极子时,R⟶0\end{array}\end{array}$$

10、(一次场)均匀全空间球坐标系下的亥姆霍兹方程（Helmholtz Equation）

因为空间上的对称性，考虑矢量势只具有单一的$z$分量，在球坐标系下只与变量$R$有关，因此可以表示为：
$$\begin{array}{c}\mathbf{A}={\mathbf{A}}_z(k,R)\end{array}$$ 拉普拉斯算子在球坐标系下表示为：
$$\begin{array}{c}{\nabla }^2=\frac{1}{R^2}\frac{\partial }{\partial R}\left(R^2\frac{\partial }{\partial R}\right)+\frac{1}{R^2\sin \theta }\frac{\partial }{\partial \theta }\left(\sin \theta \frac{\partial }{\partial \theta }\right)+\frac{1}{R^2{\sin }^2\theta }\frac{{\partial }^2}{\partial {\varphi }^2}\end{array}$$忽略掉角项$(\theta ,\varphi )$部分： $$\begin{array}{c}{\nabla }^2=\frac{1}{R^2}\frac{\partial }{\partial R}\left(R^2\frac{\partial }{\partial R}\right)\end{array}$$代入到$(1)$的亥姆霍兹方程中去，由于微分项减少到一个，偏微分就转变未来常微分，可以得到： $$\begin{array}{c}\frac{1}{R^2}\frac{d}{dR}\left(R^2\frac{d{A}_z}{dR}\right)+k^2{A}_z=0\end{array}$$为了求解这个方程需要定义一个新的变量$W=A_{z}R$来帮助我们求解，通过该变量和$(5)$方程中第一项的求导性质： $$\begin{array}{c}\begin{array}{rl}{A}_z& =\frac{W}{R}\\ \frac{d{A}_z}{dR}& =-R^{-2}W+R^{-1}\frac{dW}{dR}\\ R^2\frac{d{A}_z}{dR}& =-W+R\frac{dW}{dR}\\ \frac{d}{dR}\left(R^2\frac{d{A}_z}{dR}\right)& =R\frac{d^{2}W}{d{R}^2}\end{array}\end{array}$$代入到$(32)$式中去，得： $$\begin{array}{c}\begin{array}{r}\frac{1}{R}\frac{d^{2}W}{d{R}^2}+\frac{1}{R}{k}^{2}W=0\\ \frac{d^{2}W}{d{R}^2}+k^{2}W=0\end{array}\end{array}$$上述$(34)$式是普通的二阶线性其次常微分方程，它有指数形式的通解： $$\begin{array}{c}W=\{\begin{array}{l}{e}^{ikR},\\ e^{-ikR}\end{array}\end{array}$$ 最终得到原亥姆霍兹方程的通解为： $$\begin{array}{c}{A}_z=C\frac{e^{ikR}}{R}+D\frac{e^{-ikR}}{R}\end{array}$$考虑$k$的形式为$k=a+ib$, $$-ikR=-i(a+ib)R=-iaR+bR$$当波传到无穷远时，矢势为无穷大，所以舍弃掉第二项，最终得到：
$$\begin{array}{c}{A}_z^{\ast }(k,R)=C\frac{e^{ikR}}{R}\end{array}$$ 在靠近偶极子时，场的主要由偶极子源引起，即当$R⟶0$时，$E_{\varphi }^{\ast }$应该为偶极子产生场的大小，考虑这个边界条件，由（26）式可以得到：
$$\begin{array}{c}\begin{array}{rl}& E_{\varphi }(t)=-{\mu }_0\frac{\dot{M}(t)}{4\pi R^2}\sin \theta (26)\\ & E_{\varphi }^{\ast }(\omega )=\frac{i\omega {\mu }_0{M}_0}{4\pi R^2}\sin \theta \\ & E_{\varphi }^{\ast }=\frac{C}{R^2}\sin \theta =\frac{i\omega {\mu }_0{M}_0}{4\pi R^2}\sin \theta \\ & C=\frac{i\omega {\mu }_0{M}_0}{4\pi }\end{array}\end{array}$$因此，矢量势的复振幅为 $$\begin{array}{c}\begin{array}{r}{\mathbf{A}}^{=}\frac{i\omega {\mu }_0{M}_0}{4\pi }{e}^{ikR}\mathbf{l}z\end{array}\end{array}$$因此，我们得到均匀全空间中电偶极子的矢量势的分量为：
$$\begin{array}{c}\begin{array}{r}{A}_z^{\ast }(k,R)=\frac{i\omega {\mu }_0{M}_0}{4\pi }\frac{e^{ikR}}{R}=\frac{i\omega {\mu }_0{M}_0}{4\pi }{\int }_0\frac{m}{m_2}{e}^{(-m2|z-h|)}{J}_0(mr)dm\end{array}\end{array}$$其中，后续等式运用了索末菲恒等式。 考虑（28）式，对于电场的复振幅，我们有：
$$\begin{array}{c}\begin{array}{r}{E}_{\varphi }^{\ast }=\frac{i\omega {\mu }_0{M}_0}{4\pi }(1-ikR)e^{ikR}\sin \theta \end{array}\end{array}$$

11、(二次场)亥姆霍兹方程的通解

在柱坐标系下，亥姆霍兹方程具有以下形式：
$$\begin{array}{c}\begin{array}{rl}& \frac{{\partial }^2{A}_z^{\ast }}{\partial r^2}+\frac{1}{r}\frac{\partial A_Z^{\ast }}{\partial r}+\frac{\partial A_z^{\ast }}{\partial z^2}{k}^2{A}_z^{\ast }=0\\ & 其中，\frac{\partial A_z^{\ast }}{\partial \varphi }=0\end{array}\end{array}$$利用分离变量法来求解这个偏微分方程，假定：
$$\begin{array}{c}{A}_z^{\ast }=U(r)V(z)\end{array}$$代入到原式中，进行分离变量后，可以得到两个常微分方程如下：
$$\begin{array}{c}\begin{array}{r}\frac{{\partial }^{2}U}{\partial r^2}+\frac{1}{r}\frac{\partial U}{\partial r}+m^{2}U=0\\ \frac{{\partial }^{2}V}{\partial z^2}-(m^2-k^2)V=0\end{array}\end{array}$$其中，$m$为分离的变量。第一个方程的解是第一类和第二类的贝塞尔函数$J_0(mr),Y_0(mr)$，考虑到$Y_0(mr)$在$z$轴上有无穷大的极值，因此把它给舍弃掉。第二个方程式二阶线性常微分方程，使用特征根法求解得到：
$$\begin{array}{c}V=C{e}^{-(m^2-k^2{)}^{1/2}z}+D{e}^{(m^2-k^2{)}^{1/2}z}\end{array}$$考虑到电磁场的有限性，亥姆霍兹方程通解的最终形式为：
$$\begin{array}{c}{A}_z^{\ast }=\frac{i\omega {\mu }_0{M}_0}{4\pi }{\int }_0^{\infty }{C}_m{e}^{-m_{i}z}{J}_0(mr)dm\end{array}$$

12、在均匀半空间表面上的水平电偶极子的谐波场

假设一个水平电偶极子位于如图所示的下半部分空间，在地球表面下的距离h处。

存在水平边界时的水平电偶极子

矢量势的复振幅在任意位置处服从亥姆霍兹方程：
$$\begin{array}{c}\begin{array}{rl}{\nabla }^2{\mathbf{A}}^{\ast }+k& {}^2{\mathbf{A}}^{\ast }=0\\ {\mathbf{B}}^{\ast }=\nabla \times {\mathbf{A}}^{\ast },{\mathbf{E}}^{\ast }& =i\omega {\mathbf{A}}^{\ast }+\frac{1}{\gamma {\mu }_0}\nabla \nabla \cdot {\mathbf{A}}^{\ast }\end{array}\end{array}$$我们假设矢量势有两个分量：一个分量指向偶极矩，另一个分量垂直于界面：
$$\begin{array}{c}\begin{array}{r}{\mathbf{A}}^{\ast }=\left(A_x^{\ast },0,A_z^{\ast }\right)\end{array}\end{array}$$在笛卡尔坐标系中，边界条件可以写成
$$\begin{array}{c}\begin{array}{r}{E}_{1x}^{\ast }=E_{2x}^{\ast },E_{1y}^{\ast }=E_{2y}^{\ast }\\ B_{1x}^{\ast }=B_{2x}^{\ast },B_{1y}^{\ast }=B_{2y}^{\ast }\end{array}\end{array}$$为了满足矢量势的这些条件，在表面位置上存在以下等式：

$$\begin{array}{c}\begin{array}{rl}i\omega A_{1x}^{\ast }+\frac{1}{{\gamma }_1{\mu }_0}\frac{\partial }{\partial x}\nabla \cdot {\mathbf{A}}_1^{\ast }& =i\omega A_{2x}^{\ast }+\frac{1}{{\gamma }_2{\mu }_0}\frac{\partial }{\partial x}\nabla \cdot {\mathbf{A}}_2^{\ast }\\ \frac{1}{{\gamma }_1{\mu }_0}\frac{\partial }{\partial y}\nabla \cdot {\mathbf{A}}_1^{\ast }& =\frac{1}{{\gamma }_2{\mu }_0}\frac{\partial }{\partial y}\nabla \cdot {\mathbf{A}}_2^{\ast }\end{array}\end{array}$$

$$\begin{array}{c}\begin{array}{rl}{\mathbf{B}}^{\ast }& =\nabla \times {\mathbf{A}}^{\ast }\\ & =\left|\begin{array}{ccc}\mathbf{i}& \mathbf{j}& \mathbf{k}\\ \frac{\partial }{\partial x}& \frac{\partial }{\partial y}& \frac{\partial }{\partial z}\\ A_x^{\ast }& 0& A_z^{\ast }\end{array}\right|\end{array}\end{array}$$

$$\begin{array}{c}\begin{array}{rl}\frac{\partial A_{1z}^{\ast }}{\partial y}& =\frac{\partial A_{2z}^{\ast }}{\partial y}\\ \frac{\partial A_{1x}^{\ast }}{\partial z}-\frac{\partial A_{1z}^{\ast }}{\partial x}& =\frac{\partial A_{2x}^{\ast }}{\partial z}-\frac{\partial A_{2z}^{\ast }}{\partial x}\end{array}\end{array}$$通过化简，这些条件等价于两组非常简单的条件：
$$\begin{array}{c}\begin{array}{rl}{A}_{1x}^{\ast }=A_{2x}^{\ast }& ,\frac{A_{1x}^{\ast }}{\partial z}=\frac{A_{2x}^{\ast }}{\partial z}\\ A_{1z}^{\ast }=A_{2z}^{\ast }& ,\frac{1}{{\gamma }_1}\nabla \cdot {\mathbf{A}}_1^{\ast }=\frac{1}{{\gamma }_2}\nabla \cdot {\mathbf{A}}_2^{\ast }\end{array}\end{array}$$

12.1、$A_x^{\ast }$和 $A_z^{\ast }$的积分表示

接下来确定了向量势的垂直分量$A_x^{\ast }$的一个表达式，它也满足亥姆霍兹方程：
$$\begin{array}{c}\{\begin{array}{l}{\nabla }^2{A}_{1x}^{\ast }+k_1^2{A}_{1x}^{\ast }=0,当z<0\\ {\nabla }^2{A}_{2x}^{\ast }+k_2^2{A}_{2x}^{\ast }=0,当z>0\end{array}\end{array}$$同时满足边界条件，
$$\begin{array}{c}{A}_{1x}^{\ast }=A_{2x}^{\ast },\frac{\partial A_{1x}^{\ast }}{\partial z}=\frac{A_{2x}^{\ast }}{\partial z},当z=0\end{array}$$在下半部分空间中，分量$A_{2x}^{\ast }$可以写成一个和如下：
$$\begin{array}{c}\begin{array}{r}{A}_{2x}^{\ast }=A_{2x}^{{(0)}^{\ast }}+A_{2x}^{{(s)}^{\ast }}\end{array}\end{array}$$其中，$A_{2x}^{{(0)}^{\ast }}$为均匀全空间中电偶极子的矢量势的分量，$A_{2x}^{{(s)}^{\ast }}$为表示二次场的矢量势项。利用之前的证明：
$$\begin{array}{c}\begin{array}{r}{A}_{2x}^{{(0)}^{\ast }}=\frac{{\mu }_{0}Idl}{4\pi }\frac{e^{i{k}_{2}R}}{R}=\frac{{\mu }_{0}Idl}{4\pi }{\int }_0^{\infty }\frac{m}{m_2}{e}^{(-m2|z-h|)}{J}_0(mr)dm\end{array}\end{array}$$这里有，
$$\begin{array}{c}\begin{array}{r}R={\left[r^2+(z-h{)}^2\right]}^{1/2},m_2={\left(m^2-k_2^2\right)}^{1/2}\end{array}\end{array}$$再利用之前求解亥姆霍兹方程的结果，可知：
$$\begin{array}{c}\begin{array}{rl}& A_{1x}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }{\int }_0^{\infty }{C}_1{e}^{(m_{1}z}{J}_0(mr)dm\\ & A_{2x}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }{\int }_0^{\infty }\left[\frac{m}{m_2}{e}^{-m_2|z-h|)}+D_1{e}^{-m_{2}z}\right]J_0(mr)dm\end{array}\end{array}$$代入边界条件，我们得到了一个确定$C_1$和$D_1$的方程组:
$$\begin{array}{c}\begin{array}{rl}& C_1=\frac{m}{m_2}{e}^{-m_{2}h}+D1\\ & m_1{C}_1=m{e}^{-m_{2}h}-m_2{D}_1\end{array}\end{array}$$解这个二元方程组，得
$$\begin{array}{c}\begin{array}{r}{C}_1=\frac{2m}{m_1+m_2}{e}^{-m_{2}h},D_1=\frac{m}{m_2}\frac{m_2-n_1}{m_2+m_1}{e}^{-m_{2}h},m_1=(m^2-k_1^2{)}^{1/2}\end{array}\end{array}$$因此，我们得到了水平分量的以下表达式:
$$\begin{array}{c}\begin{array}{rl}& A_{1x}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }{\int }_0^{\infty }\frac{2m}{m_1+m_2}{e}^{-m_{2}h}{e}^{m_{1}z}{J}_0(mr)dm\\ & A_{2x}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }{\int }_0^{\infty }\frac{m}{m_2}\left\{e^{-m_2|z-h|}+\frac{m_2-m_1}{m_2+m_1}{e}^{-m_2(z+h)}\right\}{J}_0(mr)dm\end{array}\end{array}$$
利用贝塞尔函数的微分性质，同理可得
$$\begin{array}{c}\begin{array}{r}{A}_{1z}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }({\rho }_1-{\rho }_2)\cos \varphi {\int }_0^{\infty }\frac{2m^2{e}^{-m_{2}h}{e}^{m_{1}z}}{(m_1+m_2)({\rho }_1{m}_1+{\rho }_2{m}_2)}{J}_1(mr)dm\\ A_{2z}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }({\rho }_1-{\rho }_2)\cos \varphi {\int }_0^{\infty }\frac{2m^2{e}^{-m_{2}h}{e}^{-m_{2}z}}{(m_1+m_2)({\rho }_1{m}_1+{\rho }_2{m}_2)}{J}_1(mr)dm\end{array}\end{array}$$
讨论偶极子位于地表时，上半部分空间是一个绝缘体，即$h=0$且$m_1=m,m_2=m_1$在这个条件下，矢量势的分量为：
$$\begin{array}{c}\begin{array}{rl}& \{\begin{array}{ll}& A_{1x}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }{\int }_0^{\infty }\frac{2m}{m_1+m_2}{e}^{m_{1}z}{J}_0(mr)dm\\ & A_{1z}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }\cos \varphi {\int }_0^{\infty }\frac{2m}{m+m_1}{e}^{mz}{J}_1(mr)dm\end{array},当z<0\\ & \{\begin{array}{ll}& A_{2x}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }{\int }_0^{\infty }\frac{2m}{m+m_1}{e}^{-m_{1}z}{J}_0(mr)dm\\ & A_{2z}^{\ast }=\frac{{\mu }_{0}Idl}{4\pi }\cos \varphi {\int }_0^{\infty }\frac{2m}{m+m_1}{e}^{-m_{1}z}{J}_1(mr)dm,当z>0\end{array}\end{array}\end{array}$$

12.2、地球表面的电磁场

电偶极源电磁场的分布特点

12.3、近区

近区（A） (近场， 感应场微弱)靠近供电偶极，电场水平分量正比于地下电阻率，且与频率无关。近区视电阻率是接-发距r的函数。近区电场E按1/3衰减，磁场H按1/2衰减。近区测量结果与直流电阻率测深相类似

12.4、远区

远区 ©与大地电磁法一样，场是电阻率和频率的函数。
远区的卡尼亚电阻率和探测深度分别由下两式表示：

$$\begin{array}{c}\rho =\frac{1}{5}f\times \left|\frac{E_x}{E_y}\right|,h\approx 356\times \sqrt{\frac{\rho }{f}}\end{array}$$