本发明涉及地球物理技术领域,特别是涉及一种全方位观测的emn广域电磁法。
背景技术:
广域电磁法形式多样,目前以电偶源下观测电场的e-ex法应用最为广泛,近年来,e-ex广域电磁法在油气、矿产、地热、水文、工程等方面的应用越来越多。e-ex法在实际施工中,要求测量电极mn与发射电极ab的夹角小于2°,但受地形、地物的影响,有时无法将mn布置得与ab平行,这样观测到的电场并不是ex,而是含有其他方向电场的矢量值。
技术实现要素:
针对上述现有技术的不足,本发明提供了一种全方位观测的emn广域电磁法,在测量电极mn与发射电极ab不平行的时候也可以进行测量。
为了解决上述技术问题,本发明采用了如下的技术方案:
一种全方位观测的emn广域电磁法,包括以下步骤:
s1、测出测量端的测量电极mn与以发射端的发射电极ab为x轴的夹角α;
s2、测量出水平电偶源场中任意方向的电场分量emn1;
s3、对任意方向的电场分量emn1进行静态校正得到电场分量emn2;
s4、用迭代的方式求出静态校正后的电场分量emn2对应的全域视电阻率,并用于反演处理中。
其中,水平电偶源在地表处圆柱坐标系
由公式
即
其中,d=(1+ikr)e-ikr;i为发射端的供电电流;l为偶极矩;ρ为观测点的电阻率;k为波数;r为收发距,即电偶源至观测点的距离;
作为优化,步骤s3中,对实测的电场分量emn1进行静态校正的方法为空间滤波法。
作为优化,s3.1、假设观测点位于远区;设置角度系数q;
即远区的观测值近似为:
s3.2、将公式(6)转换成:
s3.3、对
作为优化,步骤s4中,求emn2的全域视电阻率的具体步骤为:
s4.1、计算出远区的视电阻率ρ(0);
s4.2、将ρ(0)带入公式(4)的右侧,求得下一次迭代视电阻率ρ(1);
s4.3、判断|ρ(1)-ρ(0)|/ρ(0)是否不大于ε,若是,则停止迭代,得到全域视电阻率;若否,则返回s4.2;其中,ε为给定误差限制。
作为优化,步骤s3.1还包括定义
本发明的有益效果是:
本发明能够在测量电极mn与发射电极ab不平行的时候也可以进行测量并进行数据处理,应用更广泛。
附图说明
图1为本发明所述的一种全方位观测的emn广域电磁法的方法流程图。
图2为本发明所述的一种全方位观测的emn广域电磁法的坐标系示意图。
图3为角度系数q的二维等值线图。
具体实施方式
下面结合附图对本发明作进一步的详细说明。
在进行具体实施例的说明之前,先解释一下:
物探的数据处理的方法,一般是基于均匀半空间的前提下进行的,再拓展到水平层状地层型等。
电偶源即发射端,主要是指发射端的尺寸相对于测量端与发射端之间的距离,可以认为很小,所以可以认为发射端就是电偶源。观测点为测量电极mn中点。
而偶极矩则是大小等于发射电极距ab、方向与发射端电流i方向一致的矢量。
本发明中的所有公式符号均通用。
如图1所示,一种全方位观测的emn广域电磁法,包括以下步骤:
s1、测出测量端的测量电极mn与以发射端的发射电极ab为x轴的夹角α。夹角α可以通过测量两个发射电极a、b以及两个测量电极m、n的坐标得到。
如图2所示,s2、测量出水平电偶源中任意方向的电场分量emn1。
均匀半空间表面水平水平电偶源在地表处圆柱坐标系
由坐标变换公式
即
其中,d=(1+ikr)e-ikr;i为发射端的供电电流;l为偶极矩;ρ为观测点的电阻率;k为波数;r为收发距,即电偶源至观测点的距离;
电场分量emn1可由现有的电场强度测试仪测出。
由现有公式k2=ω2με-iωμσ,σ=1/ρ,μ为大地导磁率;ω为发射端的发射电流圆频率可知,d参数中含有电阻率ρ和频率参数ω,实质就是电磁响应的特征。根据公式(3)可知emn1与d参数有关,也就是说,观测emn1也可以进行频率测深。
s3、对任意方向的电场分量emn1进行静态校正的方法为空间滤波法,得到电场分量emn2。具体为:
s3.1、假设观测点位于远区;设置角度系数q,该系数为无量纲;
在远区,d≈0,即公式(3)和公式(5)结合可近似得到远区的电场分量:
这里的远区是指观测点位于以电偶源为中心,收发半径使得d=(1+ikr)e-ikr<0.02的区域;达到远区要求的收发距r,与频率ω、电阻率ρ有关,是一个相对复杂的关系式。
s3.2、将公式(6)转换成:
s3.3、对
s4、用迭代的方式求出静态校正后的电场分量emn2对应的全域视电阻率,并用于反演处理中。
求emn2的全域视电阻率的具体步骤为:
s4.1、计算出远区的视电阻率ρ(0);
s4.2、将ρ(0)带入公式(4)的右侧,求得下一次迭代视电阻率ρ(1);
s4.3、判断|ρ(1)-ρ(0)|/ρ(0)是否不大于ε,若是,则停止迭代,得到第一全域视电阻率;若否,则返回s4.2。ε为给定误差限制。若误差小于规定值,停止迭代。若误差大于规定值,令ρ(0)=ρ(1),返回s4.2,直到满足要求。
通过以上步骤计算出正确的全域视电阻率,再通过成熟的反演方法,来达到勘探的目的。
在观测方位角
为了简化讨论,本实施例假定观测点在远区,得到远区计算式(6),可知电场分量emn1的大小与
将公式(5)定义为q函数,在观测方位角
从图3中可以看出,当
从图3中可以判断出,在其他观测方位角
另外,在e-ex法中通过验证ey的表达式已经可以得出ey与频率无关,与几何尺寸有关,不能通过观测ey来进行频率测深,表明ey不具备频率勘探的可能,因此,应该避开ab与mn夹角大于75°或小于-75°的区域,该区域主要观测的是ey分量。
所以,最终定义q=0的区域为观测禁区,|q|<0.75的区域为观测盲区。
野外观测过程中的角度关系应选择在图3中满足|q|>0.75的区域。
实际勘探中,例如对水平层状地层进行勘探,假设水平层状地层在水平电偶源中,地表观测点处的电场分量er和
其中,
通过用括号+下标0表示贝塞尔函数零阶积分,括号+下标1表示贝塞尔函数一阶积分可得到公式(8)和公式(9);
则水平层状地层中任意方向的电场分量为:
其中,ω为发射端的发射电流圆频率;μ为大地导磁率;m为空间频率;r*和r为空间频率特性函数;h为各层的分层厚度;n为地层的层数;σ各层的为大地电导率;i为发射端的供电电流;l为电偶源极距的长度;i为纯虚数;k为波数。
通过计算公式(10)可以计算出水平层状地层水平电偶源的电场分量emn′,通过步骤s4计算出正确的全域视电阻率,将计算的全域视电阻率与实际测量并静态校正后的电场分量emn2得到的全域电阻率进行拟合反演,求取水平层状地层型的参数即可得到地下的电性模型,达到反演的目的。
最后应说明的是:本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等统计数的范围之内,则本发明也意图包含这些改动和变型。