三维感应仪器系数

1 感应电动势

三维阵列感应测井中，在各向同性介质中单发单收线圈结构中，主分量$xx,yy,zz$接收线圈中的感应电动势分别为

$$V_{zz}=I_{z}\frac{i\omega \mu (A_{TZ}A_{RZ})(N_{TZ}N_{RZ})}{2\pi L^{3}} (1-ikL)e^{ikL}$$

$$V_{xx}=I_{x}\frac{i\omega \mu (A_{TX}A_{RX})(N_{TX}N_{RX})}{4\pi L^{3}}(-1+ikL+k^{2}L^{2})e^{ikL}$$

$$V_{yy}=V_{xx}$$
其中， $A_{TZ},A_{TX},A_{RZ},A_{RX}$分别为Z和X方向发射和接收线圈的有效面积； $N_{TZ},N_{TX},N_{RZ},N_{RX}$分别为Z和X方向发射和接收线圈的匝数； $k$为传播系数，$k=\sqrt{i\omega \mu \sigma}=\frac{1+i}{\delta}$， $\delta=\sqrt{\frac{2}{\omega \mu \sigma}}$为趋肤深度，它表示电磁场在介质中传播时的穿透深度。

2 仪器系数推导

2.1 $zz$分量仪器系数

将$k=\frac{1+i}{\delta}$分别代入感应电动势$V_{zz}$表达式，经过化简整理可得：

$$V_{zz}=I_{z} \frac{\omega \mu (A_{TZ}A_{RZ})(N_{TZ}N_{RZ})}{2\pi L^{3}} \left\{ -e^{-\frac{L}{\delta}} \left[ sin \left( \frac{L}{\delta}\right) \left( 1+\frac{L}{\delta} \right) -\frac{L}{\delta}cos(\frac{L}{\delta})\right] + ie^{-\frac{L}{\delta}} \left[ cos\left( \frac{L}{\delta}\right)\left( 1+\frac{L}{\delta}\right)+\frac{L}{\delta}sin\left( \frac{L}{\delta}\right)\right] \right\}$$
对上式展开为 $\frac{L}{\delta}$的幂级数形式，

$$V_{zz}=I_{z} \frac{\omega \mu (A_{TZ}A_{RZ})(N_{TZ}N_{TR})}{2\pi L^{3}} \left\{ \left[ -\left( \frac{L}{\delta} \right)^{2}+\frac{2}{3} \left( \frac{L}{\delta} \right)^{3} -\frac{2}{15} \left( \frac{L}{\delta} \right)^{5} +... \right]+ i\left[ 1-\frac{2}{3} \left( \frac{L}{\delta} \right)^{3} +\frac{1}{2} \left( \frac{L}{\delta} \right)^{4}+... \right] \right\}$$
从上式幂级数展开式 $V_{zz}$实部展开式第一项可以看出，实部信号与地层电导率呈线性变化，且 $V_{zz}^{R} \propto \left( \frac{1}{L} \right)$； $V_{zz}$虚部展开式第一项为直藕信号且 $V_{zz}^{X} \propto \left( \frac{1}{L^{3}}\right)$，当地层电导率 $\sigma_{f}$较小时，虚部信号基本上全部为直藕信号，可以作为无用信号舍弃，当地层电导率 $\sigma_{f}$比较大时， $V_{zz}$虚部反映地层电导率。

上式提取$\left( \frac{L}{\delta} \right)^{2}=\frac{L^{2}\omega \mu \sigma}{2}$，可得

$$V_{zz}=I_{z} \frac{(\omega \mu)^{2} (A_{TZ}A_{RZ})(N_{TZ}N_{TR})}{4\pi L} \sigma \left[ \left( -1+\frac{2}{3}\frac{L}{\delta}+... \right)+i \frac{\delta ^{2}}{L^{2}} \left( 1-\frac{2}{3}\left( \frac{L}{\delta} \right)^{3}+... \right) \right]$$
视电导率定义为

$$\sigma_{zz}=-\frac{V_{zz}/ I_{zz}}{K_{zz}}$$
则仪器系数为(当地层电导率 $\sigma_{f}$很小时， $\delta$很大，级数展开式后面项可以忽略。仪器所测量的视电导率 $\sigma_{zz}$非常接近地层电导率 $\sigma_{f}$。)

$$K_{zz}=\frac{(\omega \mu)^{2} (A_{TZ}A_{RZ})(N_{TZ}N_{TR})}{4\pi L}$$
由前所述，视电导率 $\sigma_{zz}$的（减去直藕信号）表达式为

$$\sigma_{zz}=\frac{-2i}{\omega \mu L^{2}}[(1-ikL)e^{ikL}-1]$$
视电导率 $\sigma_{zz}$对实部和虚部展开为 $\frac{L}{\delta}$的幂级数形式为

$$\sigma_{zz}^{R}=\sigma \left[ 1-\frac{2}{3} \left( \frac{L}{\delta} \right)+\frac{2}{15} \left( \frac{L}{\delta} \right)^{3}-\frac{1}{18} \left( \frac{L}{\delta} \right)^{4} +... \right]$$
$$\sigma_{zz}^{X}=\sigma \left[ \frac{2}{3} \left( \frac{L}{\delta} \right)-\frac{1}{2} \left( \frac{L}{\delta} \right)^{2} +\frac{2}{15} \left( \frac{L}{\delta} \right)^{3}+... \right]$$

2.2 $xx,yy$分量仪器系数

同理，将$k=\frac{1+i}{\delta}$分别代入感应电动势$V_{xx}$表达式，经过化简整理可得：

$$V_{xx}=I_{x}\frac{\omega \mu (A_{TX}A_{RX})(N_{TX}N_{RX})}{4\pi L^{3}} \left\{ -ie^{ -\frac{L}{\delta} } \left[ (1+\frac{L}{\delta}) cos\left( \frac{L}{\delta} \right) +\left( 2\left( \frac{L}{\delta} \right)^{2} + \frac{L}{\delta}\right) sin\left( \frac{L}{\delta} \right) \right] +e^{ -\frac{L}{\delta} } \left[ sin\left( \frac{L}{\delta} \right) \left( 1+\frac{L}{\delta}\right)-cos\left( \frac{L}{\delta} \right) \left( \frac{L}{\delta} +2 \left( \frac{L}{\delta} \right)^{2} \right) \right]\right\}$$
对上式展开为 $\frac{L}{\delta}$的幂级数形式，

$$V_{xx}=I_{x}\frac{\omega \mu (A_{TX}A_{RX})(N_{TX}N_{RX})}{4\pi L^{3}} \left\{ \left[ -\left( \frac{L}{\delta} \right)^{2}+\frac{4}{3}\left( \frac{L}{\delta} \right)^{3}-\frac{8}{15}\left( \frac{L}{\delta} \right)^{5}+...\right] -i \left[ -1+\frac{4}{3}\left( \frac{L}{\delta} \right)^{3}+\frac{3}{2}\left( \frac{L}{\delta} \right)^{4}+... \right] \right\}$$
上式提取 $\left( \frac{L}{\delta} \right)^{2}=\frac{L^{2}\omega \mu \sigma}{2}$，可得

$$V_{xx}=I_{x}\frac{(\omega \mu)^2 (A_{TX}A_{RX})(N_{TX}N_{RX})}{8\pi L} \sigma \left\{ \left[-1+\frac{4}{3}\left( \frac{L}{\delta} \right)+... \right] -i\frac{\delta^2}{L^2} \left[ -1-\frac{4}{3}\left( \frac{L}{\delta} \right)^{3} +...\right]\right\}$$
由视电导率定义可得 $xx,yy$仪器系数为

$$K_{xx}=\frac{(\omega \mu)^2 (A_{TX}A_{RX})(N_{TX}N_{RX})}{8\pi L}$$
$$K_{yy}=\frac{(\omega \mu)^2 (A_{TY}A_{RY})(N_{TY}N_{RY})}{8\pi L}$$

2.3 与$z$相关交叉分量仪器系数

由于$zx,xz,zy,yz$交叉分量在各向同性地层中为零，无法用前述的方法进行推导。这里给出其仪器系数分别为（给出仪器系数的条件为在低电导率两层地层的边界处$\sigma_{xz},\sigma_{zx}=\pm(\sigma_{1}-\sigma_{2})$）

$$K_{zx}=\frac{(\omega \mu)^2 (A_{TZ}A_{RX})(N_{TZ}N_{RX})}{16\pi L}$$
$$K_{xz}=\frac{(\omega \mu)^2 (A_{TX}A_{RZ})(N_{TX}N_{RZ})}{16\pi L}$$
$$K_{yz}=\frac{(\omega \mu)^2 (A_{TY}A_{RZ})(N_{TY}N_{RZ})}{16\pi L}$$
$$K_{zy}=\frac{(\omega \mu)^2 (A_{TZ}A_{RY})(N_{TZ}N_{RY})}{16\pi L}$$

2.4 $xy,yx$交叉分量仪器系数

类似$xx,yy$分量，这里直接给出$xy,yx$交叉分量仪器系数

$$K_{xy}=\frac{(\omega \mu)^2 (A_{TX}A_{RY})(N_{TX}N_{RY})}{8\pi L}$$

$$K_{yx}=\frac{(\omega \mu)^2 (A_{TY}A_{RX})(N_{TY}N_{RX})}{8\pi L}$$