[发明专利]一种多辐射场源瞬变电磁法多分量全域视电阻率定义方法

权利要求书

1.一种多辐射场源瞬变电磁法多分量全域视电阻率定义方法，其特征在于，包括以下步骤：

1)磁场强度的x和z分量

将磁场强度的x和z分量记为B_p(ρ，C，t)(p＝x，z)，其中C表示测点在空中的位置坐标参数。给定初值在的邻域内对B_p(ρ，C，t)(p＝x，z)进行泰勒展开

Bp(ρ,C,t)=Bp(ρτ(0),C,t)+Bp′(ρτ(0),C,t)(ρ-ρτ(0))+Bp′′(ρτ(0),C,t)2!(ρ-ρτ(0))2+...+Bp(n)(ρτ(0)C,t)n!(ρ-ρτ(0))n+Rn(ρ)(p=x,z)---(7)]]>

由于该邻域非常小，在该邻域内用一小段直线对曲线进行近似，即保留(7)式的前两项，有

Bp(ρ,C,t)≈Bp(ρτ(0),C,t)+Bp′(ρτ(0),C,t)(ρ-ρτ(0))(p=x,z)---(8)]]>

对上式进行变换

ρ=Bp(ρ,C,t)-Bp(ρτ(0),C,t)Bp′(ρτ(0),C,t)+ρτ(0)(p=x,z)---(9)]]>

将(9)式写成迭代的形式

ρτ(i+1)≈Δρτ(i)+ρτ(i)(i=0,1,2,...)---(10)]]>

其中

Δρτi=Bp(ρ,C,t)-Bp(ρτ(i-1),C,t)Bp′(ρτ(i-1),C,t)(p=x,z)---(11)]]>

(10)式的迭代终止条件为

|Bp(ρ,C,t)-Bp(ρτ(i),C,t)Bp(ρ,C,t)|<ϵ(p=x,z)---(12)]]>

ε是给定的迭代终止误差限，B_p(ρ，C，t)(p＝x，z)是测得的磁场分量；

2)磁场强度的y分量

ρτ(k)(i)=-16{ρτ(k-1)(i+2)+ρτ(k)(i-2)-4[ρτ(k-1)(i+1)+ρτ(k)(i-1)]}(i=1,...,M)]]>

ρτ(i-1)=2ρτ(i)-ρτ(i+1)(i=1)ρτ(i+1)=2ρτ(i)-ρτ(i-1)(i=M)ρτ(i-2)=ρτ(i+2)-2[ρτ(i)-ρτ(i-1)](i=1)ρτ(i+2)=ρτ(i-2)+2[ρτ(i)-ρτ(i-1)](i=M)---(13)]]>

其中k表示插值的次数，M表示参与插值的总个数，一般在需要插值的区域左右两端各取2个点参与计算。