1.一种探测射孔金属套管外地层电阻率的井中时域电磁法，其特征在于，包括如下步骤：S1，根据完井和射孔资料，在射孔层段，给井下装置的发射线圈通以脉冲电流，产生瞬时大功率的脉冲磁场激发，在距发射线圈一定距离处用接收线圈记录套管井中磁场垂直分量Hz0产生的感生电动势ε0(t)随时间变化的数据，对发射的电流波形和接收的感生电动势信号ε0(t)进行全时段数字记录；

S2，根据S1中井的资料，建立径向分为井中流体、套管、套外地层的三层介质模型，并采用与S1中一致的脉冲电流模拟激发脉冲磁场；推导柱坐标系下所述三层介质模型的频域方程组，利用快速余弦变换算法计算频域解，再用快速正弦变换算法将频域解转换为时域解，该时域解的磁场垂直分量Hzb产生的感生电动势εb(t)，与S1中记录的感生电动势ε0(t)的差值，即为异常场观测信号Δε0(t)：Δε0(t)＝ε0(t)‑εb(t)      (1)；

S3，在S2中所述的三层介质模型中添加S1中井的射孔和套外电阻率异常体，构成三维地层模型，并采用与S1中一致的脉冲电流模拟激发脉冲磁场；以S2中所述三层介质模型为背景，射孔和套管外异常体可视作异常场源，根据异常场分离技术，将所述三维地层模型的磁场分解成背景场和异常场之和，通过进一步代换得到异常场的磁场双旋度方程，将该方程按照柱坐标系分解成各个方向的微分方程组；

S4，对S3中所述的三维地层模型进行交错网格剖分，将S3中得到的微分方程组，转换为差分方程组，引入稳定性条件和边界条件，利用时域有限差分法计算三维地层模型的异常场；所述三维地层模型的异常场的磁场垂直分量Hza产生感生电动势εa(t)，与S2中计算的异常场观测信号Δε0(t)进行比较，若计算值与观测值的剩余值达到预置精度，则所述三维地层模型即为反演结果；若未达到预置精度，返回至步骤S3，根据剩余值修改三维地层模型，重新计算比较，直至达到预置精度，得到反演结果；

S5，重复步骤S1～S4，得到所有射孔层段的反演电阻率值，根据所有射孔层段的反演电阻率值获得套管外储层电阻率分布的二维图像；

S6，结合测井完井资料与采油资料，由套管外储层电阻率的分布可以确定储层中剩余油的分布状况。

2.根据权利要求1所述的一种探测射孔金属套管外地层电阻率的井中时域电磁法，其特征在于：套管外地层的电阻率测量深度范围为覆盖有效储层及上下围岩地层一定厚度的范围。

3.根据权利要求1所述的一种探测射孔金属套管外地层电阻率的井中时域电磁法，其特征在于：S2中所述的εb(t)是三层介质模型的垂直磁场分量产生的感生电动势，由数值计算所得。

4.根据权利要求1所述的一种探测射孔金属套管外地层电阻率的井中时域电磁法，其特征在于：S3中所述的异常场分离技术是指，为了便于在数值模拟中实现对源的处理，将总t t场e和h表示为背景场和异常场之和，即：

t s p

e＝e+e        (2)

t s p

h＝h+h        (3)

式中，h为磁场，e为电场，上标t，s和p分别表示总场、异常场和背景场。

5.根据权利要求1所述的一种探测射孔金属套管外地层电阻率的井中时域电磁法，其特征在于：S3中所述的磁场双旋度方程为：式中，地层电阻率和磁导率参数(ρ，μ)分别是背景参数(ρb，μb)和异常参数(ρa，μa)之和；

s p p

h、h分别表示异常磁场和背景磁场，e表示背景电场。

6.根据权利要求1所述的一种探测射孔金属套管外地层电阻率的井中时域电磁法，其特征在于：S4中，引入的稳定性条件为时间步长与空间步长必须满足：式中，Δt是时间步长，σ和μmin分别是整个模型中最小电导率和最小磁导率，Δmin为最小空间步长，t为传播持续的时间。

7.根据权利要求1所述的一种探测射孔金属套管外地层电阻率的井中时域电磁法，其特征在于：S3中的磁场双旋度方程里不含电场的异常场，其磁场由该场点前后时刻的磁场以及同时刻的周边磁场加权合成，在时间轴上逐步迭代，最后推出整个时空的磁场值。

8.根据权利要求7所述的一种探测射孔金属套管外地层电阻率的井中时域电磁法，其特征在于：所述磁场双旋度方程的边界处理方法是：在r和z轴的左边界采用后向差分，右边界采用前向差分，在θ坐标轴上，边界是首尾衔接的。