2022年电磁法勘探--可控源音频大地电磁测深法 .pdf
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《2022年电磁法勘探--可控源音频大地电磁测深法 .pdf》由会员分享,可在线阅读,更多相关《2022年电磁法勘探--可控源音频大地电磁测深法 .pdf(24页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、2-216 2.1 电磁法勘探-可控源音频大地电磁测深法(CSAMT)由于天然场源的随机性和信号微弱,MT 法需要花费巨大努力来记录和分析野外数据。为克服MT法的这个缺点,加拿大多伦多大学教授D.W.Strangway和他的学生Myron Goldstein 提出了利用人工(可控)场源的音频大地电磁法(CSAMT)。这种方法使用接地导线或不接地回线为场源,在波区测量相互正交的电、磁场切向分量,并计算卡尼亚电阻率,以保留 AMT 法的一些数据解释方法。自 20 世纪 70 年代中期,CSAMT 法得到实际应用,一些公司相继生产用于CSAMT 法测量的仪器和应用解释软件。进入 80 年代后,该方法
2、的理论和仪器得到很大发展,应用领域也扩展到普查、勘探石油、天然气、地热、金属矿产、水文、工程、环境保护等各个方面,从而成为受人重视的一种地球物理方法。虽然 CSAMT 法属于一种人工源的频率电磁深测,但和通常的频率域电磁测深不同。这主要因为CSAMT 法测量两个相互垂直的电磁场切向分量计算卡尼亚电阻率,因而具有较强的抗干扰能力,且更容易获得对地电变化较灵敏的相位差信息;又由于波区电磁场十分接近平面波,因而其资料处理、解释也较为简便,可以保留AMT法中的许多解释方法。CSAMT 和 AMT 或 MT 亦有不同,根本原因是CSAMT 法使用了人工场源,因而极化方向明显,信噪比高,易于观测。但是,由
3、于使用了人工场源,CSAMT法必然受场源效应影响,这主要包括非平面波效应、场源附加效应、阴影效应和测深通道的弯曲。2.2.1 CSAMT 基本理论CSAMT 有 2 种常用的场源水平电偶极子和垂直磁偶极子,此处注重讨论其场的特征和快速计算方法。2.2.1.1 水平层状半空间上水平如图 2.2.1 所示,N 层水平层状介质中第n 层的电阻率和层厚度分别记为n 和 hn。水平电偶极子(接地导线)位于层状介质表面,偶极矩为P=IdL(I 为谐变电流)。选取公共坐标原点位于偶极子中心的柱坐标系和直角坐标系,使 x 轴指向偶极矩方向(即=0 的方向),z 轴垂直向下,则地表面的电磁场分布可通过直接求解场
4、所满足的非齐次亥姆霍兹方程或通过求电型 Lorentz 势所满足的方程和边界条件得到。此处直接写出准静态极限下柱坐标系中地表电磁场表达式:Er=2IdLcosir0*11mm RJ1(mr)dm-101mmRJ0(mr)dm+1r01mRJ1(mr)dm(2.2.1a)E=2IdLsin1r01mRJ1(mr)dm-i0*1mmmRJ0(mr)dm+ir0*11mmRJ1(mr)dm(2.2.1b)名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页,共 24 页 -2-217 图 2.2.1 水平层状介质上的水平电偶极子Hr=-2IdLsin0*1mmmRJ1(mr)dm+r01*mR*
5、1mmmRJ0(mr)dm(2.2.1c)H=2IdLrcos0*1mmmRJ1(mr)dm(2.2.1d)Hz=2IdLsin02*1mmm RJ1(mr)dm(2.2.1e)式中:R*=cothm1h1+coth-112mmcoth(m2h2+,+coth-11NNmm)R=cothm1h1+coth-11122mmcoth(m2h2+,+coth-111NNNNmm)mj=22m-k,k2j=i/j,cothx=(ex+e-x)/(ex-e-x)特别地,当 N=1 时,可得到均匀半空间表面电磁场各分量的表达式:Er=32IdLrcos1+eikr(1-ikr)(2.2.2a)E=32Id
6、Lrsin2-eikr(1-ikr)(2.2.2b)名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页,共 24 页 -2-218 Hr=-232IdLrsin I1(2ikr)K1(2ikr)+6ikrI1(2ikr)K0(2ikr)-I0(2ikr)K1(2ikr)(2.2.2c)H=22IdLrcosI1(2ikr)K1(2ikr)(2.2.2d)Hz=-2432IdLk rsin1-eikr(1-ikr-13k2r2)(2.2.2e)以上各式中,为磁导率,为谐变电流的圆频率,I0(2ikr)、I1(2ikr)和 K0(2ikr)、K1(2ikr)分别为以2ikr为宗量的第一和第二
7、类虚宗量贝塞尔函数,下标“0”或“1”表示贝塞尔函数的阶数,k 为波数。在准静态条件下:k=(1+i)2可以看出,综合参数ikr 决定了电磁场随收发距、频率、电导率的变化规律。准静态极限下:ikr=(1+i)r/,p=r/称为电距离或感应数,它实质上是以趋肤深度为单位来表示的观察点到场源的距离。借助感应数p,可以将电磁场的变化区域分为:p1 即感应数很大的区域,称为“远区”或“平面波场区”;p1 即介于前两者之间的区域,称为“过渡带”或“中间区”。不难理解,感应数 p 的“大”或“小”,也就是场区的“远”或“近”,不但与观察点到场源的距离有关,而且与大地的电导率和使用的频率有关。在不改变发收距
8、和大地电导率 的条件下,频率 不同,趋肤深度不同,电磁场的特性也不同。当介质为低阻时,不大的发收距r 就可以满足“远区”条件,反之,在高阻介质中,即使使用高频,也不得不使 r 相当大才能满足“远区”的条件。(1)近区响应(p 1)当 r即 p1)当 r 即 p1 时,带入幂级数和贝塞尔函数的渐进展开式,可得到(2.2.2)各式的渐近表达式:Er32IdLrcos(2.2.9a)E3IdLrsin(2.2.9b)Hr43siniIdLer(2.2.9c)H-43cos2iIdLer(2.2.9d)Hz-243sin2iIdLer(2.2.9e)4ieHEZ(2.2.9f)比较(2.2.9)各式可
9、以发现,远区电磁场有如下特点:(1)以弛豫频率为极限(极化速度远小于场变化频率),电场与频率无关,但与电导率成反比,由动态平衡的“静电荷”决定;(2)磁场与电导率和频率平方根反比,对电阻率没有电场灵敏,即磁场相对均匀,是感应场;(3)垂直磁场与电导率成反比,即比水平磁场对电导率更灵敏,但衰减也更快;(4)水平电场与磁场随收发距变化规律相同,相位差是/4 或 785.4 mrad;(5)波阻抗与电阻率和频率平方根正比,与场源参数无关,与平面电磁波的波阻抗完全相同。可以定义卡尼亚电阻率:221HEc(2.2.10)当采用直角坐标时,Ex=Ercos-Esin=32IdLr(1-3sin2)(2.2
10、.12a)名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页,共 24 页 -2-220 Hy=Hrsin+Hcos=-432iIdLer(1-3sin2)(2.2.12b)4iyxxyeHEZ(2.2.13)也可以得到同样的视电阻率定义=22221|1|2|ExExHyfHy(2.2.14)实际工作中,多用MKS 制单位。此时E 以 mV/km 为单位,H 以(1=10-2/4A/m)为单位,磁导率设为真空中的值,则有:=221|5|ExfHy(2 m(2.2.16)(2.2.16)式就是以卡尼亚(Cagniard)命名的计算视电阻率的公式。卡尼亚电阻率对远区,也就是说在满足平面波的条
11、件下是有效的,并且是在 MT 法和满足远区条件的CSAMT 法中常用的关系式。(3)过渡带响应(p1)此时 r,p1,电磁场各分量由(2.2.2)式严格地描述,电磁场没有简单的初等函数表达式;电场与电导率、频率相关,受“静态”积累电荷和变化的感应磁场影响;磁场与电导率和频率相关,是感应场,受变化的电场影响;电场与磁场随收发距变化规律不同,由近区特点向远区过渡,常形成“低谷”或“凹陷”(参见图 2.2.3)。此带的波阻抗与电导率,频率及场源参数方位角都有关系,没有简单的表达式。2.2.1.2设水平层状介质分布和图2.2.1 相同,移去水平电偶极子,放置一个中心在坐标原点、磁矩为 M 的垂直磁偶极
12、子。同样选择有公共原点的一个柱坐标系和一个直角坐标系,使 z轴垂直向下,即与磁偶极子磁矩方向相同。在准静态极限下,可写出电磁场各分量表达式为:E=2Mi02*1mmmRJ1(mr)dm(2.2.17a)Hr=2M01*mR2*1mmmRJ1(mr)dm(2.2.17b)Hz=2M03*1mmmRJ0(mr)dm(2.2.17c)特别地,当 N=1 时,可得到均匀半空间表面的电磁场为:E=-42Mr3-eikr(3-3ikr-k2r2)(2.2.18a)Hr=-24MkrI1(2ikr)K1(2ikr)-I2(2ikr)K2(2ikr)(2.2.18b)名师资料总结-精品资料欢迎下载-名师精心整
13、理-第 5 页,共 24 页 -2-221 Hz=2252Mkk r9-eikr(9-9ikr-4k2r2+ik3r3)(2.2.18c)与水平电偶极子的电磁场类似,综合参数ikr 决定了电磁场随收发距、频率、电导率的变化规律。(1)近区响应(p 1)当 r即 p1)当 r即 p1 时,带入幂级数和贝塞尔函数的渐进展开式,可得到(2.2.18)各式的渐近表达式:E-432Mr(2.2.23a)Hr 432Mr(2.2.23b)Hz592Mr(2.2.23c)4ieHEZ(2.2.23d)可以看出,磁偶极子远区电磁场有如下特点:(1)频率很高时(场变化频率远快于介质极化速度),电场与频率无关,但
14、与电导率成反比;(2)磁场与电导率和频率平方根反比;(3)垂直磁场与电导率成反比,即比水平磁场对电导率更灵敏,但衰减也更快;(4)水平电场与磁场随收发距变化规律相同,相位差是/4 或 785.4 mrad。但场随 r 的增大而衰减的速度比电偶极的快;(5)波阻抗与电阻率和频率平方根正比,与场源参数无关,与水平电偶极子的完全相同,可以定义同样的卡尼亚电阻率,即:在“大”的“电距离”上,用磁偶极子做场名师资料总结-精品资料欢迎下载-名师精心整理-第 6 页,共 24 页 -2-222 源所测的电阻率与电偶极子所得的完全相同。(3)过渡带响应(p1,r)如同水平电偶极子的情形一样,磁偶极子的过渡带响
15、应由(2.2.18)式严格地描述,场及波阻抗的性质都比较复杂,与电导率、频率及场源参数有关。2.2.1.3 水平分层介质表面谐变偶极子场的数值计算方法水平层状介质表面人工源激励的电磁场表达式中都含有零阶和一阶贝塞尔函数的无穷积分,这类积分实质上都是汉克尔变换式。自20 世纪 70 年代初发展起来的快速汉克尔变换方法是计算这类积分最有效的工具之一。快速汉克尔变换考虑含 Bessel(贝塞尔)函数在(0,)区间上的积分:g(r)=0f()Jv(r)d(2.2.24)其中 Jv 是 v 阶第一类 Bessel 函数,实数 v-1,引入如下变换式:=e-x/r0 x,y(-,)(2.2.25)r=r0
16、ey式中 x,y 是快速 Hankel 变换中的新变量,r0是选定常数,定义如下新函数:F(x)=f(),G(y)=g(r)r则(2.2.24)式可重写为:G(y)=F(x)Jv(ey-x)ey-xdx=F(x)Hv(y-x)dx=F*Hv(2.2.26)即 G 是函数 F和 Hv 的褶积,这里 Hv(x)=Jv(ex)ex。根据抽样定理,对F(x)进行抽样,并利用抽样定理写出:()F x=nF(n)P(x-n)其中 为抽样间隔,且:P(x)=sin(x)/x将()F x代入(2.2.26)式,则整理得:()G y=nF(n)H*v(y-n)再对 y 离散化,上式可写为()G m=nF(n)H
17、*v(m-n)=nf(01ren)01ren H*v(m-n)(2.2.27)名师资料总结-精品资料欢迎下载-名师精心整理-第 7 页,共 24 页 -2-223 即:g(r)=1rnf(1ren)1ren H*v(m-n)(2.2.28)式中:H*v(y)=P(x)Hv(y-x)dx(2.2.29)称为汉克尔滤波系数,与抽样间隔有关。目前已有不同抽样间隔的不同长度的汉克尔滤波系数,可根据要求选用。再根据(2.2.27)式的快速变换,可以计算含贝塞尔函数的无穷积分。2.2.2 CSAMT野外工作方法和技术2.2.2.1如图 2.2.6 所示,CSAMT 包括张量、矢量和标量三种方式,取决于测量
18、分量的数量和使用场源的数量。(1)张量 CSAMT张量 CSAMT 利用两个方向的场源来测定五个分量Ex、Ey、Hx、Hy、Hz)与大地电磁场不同,CSAMT场源不是全方位的,所以需要两个场源。为了完全确定阻抗张量,总共需要测量10 个分量。图 2.2.6 张量、矢量和标量CSAMT法的定义(a),(b)用两个场源,且每个场源做5 次的全张量测量(c)每个场源观测3 个分量的部分张量测量(d)单一场源 5 次观测的矢量测量(e)单一场源2次观测的标量测量(f)用一个场源在一个方向上只做电场观测的标量CSAET 法,在少量点上测磁场用以把电场数据转换为近似的电阻率。张量测量最好用于构造很复杂的地
19、区和测深点距比地质构造尺寸大很多的地区。当测深名师资料总结-精品资料欢迎下载-名师精心整理-第 8 页,共 24 页 -2-224 点布得很近,地质构造的面貌可以靠高的测量密度来直接填图,因而张量测量的量(如倾斜度)就变得不重要了。即使是在复杂的地区,高密度的标量或矢量数据确定的构造与低密度张量结果相近。但具有强烈的局部非各向同性地区例外,在这些地区用全张量解出的阻抗可能要好些。(2)矢量 CSAMT矢量 CSAMT 利用单一个场源来测量4 个或 5个分量(Ex、Ey、Hx、Hy,有时加测Hz)。矢量 CSAMT 数据提供了关于地下二维或三维构造的信息,但比张量测量的信息少。由于少收集和处理5
20、0%的数据,矢量测量比张量便宜。矢量CSAMT 在非各向同性不强的地区确定复杂地质构造是有用的,出于经济的考虑,在二维地区常常宁可做矢量而不做张量 CSAMT。(3)标量 CSAMT标量 CSAMT 利用一个场源测量两个分量(Ex 和 Hy 或者 Ey 和 Hx),标量测量对于一维的层状条件或者走向已知的二维条件是足够了,在更复杂的条件下,则取决于数据的密度。图 2.2.7 是一个典型的标量CSAMT 示意图。在复杂地质条件下采用单一测线的标量CSAMT 是非常冒险的。例如,当偶极方向恰好垂直于断层(TM 方式)时,用标量数据可以很容易地确定线性的陡倾斜断层;可是如果偶极平行走向(TE 方式)
21、,断层的解释和定位就很困难。因此,在二维和三维地区做标量测量通常都用密的测网,密的测网能部分地克服缺少多分量数据(像矢量和张量数据中)的缺陷。主要的例外是区域性非各向同性很强的地方,在这种情况下张量或者矢量数据也许好些。标量 CSAMT 测量的主要吸引力在于它的成本相当低,这就是为什么至今取得的CSAMT 资料绝大部分是标量结果。图2.2.7 典型的标量CSAMT 测量布置示意图(4)标量 CSAET最简单的、也是目前所有商业仪器及野外所用的CSAMT 形式可称为可控源音频大地电流法(CSAET),它系统地测量电场,只在个别点测量磁场,从而把电场的测量值转换为近似的卡尼亚电阻率。CSAET 通
22、常在发射偶极中垂线15的扇形区域内测量,它也得不到有意义的相位资料。在磁场完全均匀,地质情况不复杂的层状地区用来普查填图是令人满意的。然而当这些假设条件不成立时则会引起麻烦。图2.2.8 是地堑断层的上升和下降断块上的两个一维模拟结果,图 2.2.8(a)表明,位于下降断块上的6 号点的磁场较4 号点的高,这个较大的差平方后计算卡尼亚电阻率可以引起大的电阻率误差。图2.2.8(b)显示名师资料总结-精品资料欢迎下载-名师精心整理-第 9 页,共 24 页 -2-225 在 6 号点测到的 CSAMT 和 CSAET 的响应。CSAMT 响应正确地描绘了高/低/高的地质断面,但是在这同一个点,用
23、在 CSAET 排列中点所得较低的磁场值来归一化,则得出错误的结果。层序好像是高/低/高/低/高,其最大的误差几乎达到100%。CSAET 资料在冲积盖层下“生”出了一层假高阻层。反过来,错误地选择磁场测点,其结果也可能使真实的层“消失”。如果地质条件更加复杂,可能使结果更坏。因此,除非是在极简单的地质环境中普查,在野外应避免使用CSAET。图 2.2.8 CSAET 引起的视电阻率误差(a)在 6 个偶极 CSAET 排列中,在两个测点上的磁场测量值(b)6 号测深点 CSAMT(E、H 均在 6 号点测)和CSAET(E 在 6 号测而 H 在排列中心2.2.2.2 CSAMT 场源与测量
24、模式CSAMT通常使用两种场源:接地导线(水平电偶极子)和不接地回线(垂直磁偶极子)。在均匀大地中,两种场源的远区卡尼亚电阻率与相位是相同的,但信号强度不相同。在大多数实际的野外条件下,电偶极的信噪比优于回线源,而且布置也相对方便,此外,电偶极还可以控制场的方向,所以在CSAMT 或者长距离的TEM 勘探中很少使用回线源。但在平坦的戈壁荒漠、冻土带等地区,电极接地非常困难,常使用不接地回线。不接地回线的另一个优点是场衰减快,因此,使用较小的收发距即可满足远区的要求。在使用接地导线场源时,根据发射偶极AB、接收偶极 MN 和测线布设方向相对于地质构造走向的关系,CSAMT 有 TM(Transv
25、erse Magnetic)和 TE(Transverse Electric)两种测量模式。(1)TM 模式发射偶极 AB、接收偶极MN 及测线方向垂直于地质构造走向布设。TM 模式横向分辨名师资料总结-精品资料欢迎下载-名师精心整理-第 10 页,共 24 页 -2-226 能力较强,观测的电场受静态影响、地形影响较严重。(2)TE 模式发射偶极AB、接收偶极MN 及测线方向平行于地质构造走向布设。TE 模式垂向分辨能力较高,观测的电场受静态影响、地形影响较小。2.2.2.3 CSAMT 测量覆盖范围CSAMT 测量中,有限场源的使用对在平面上允许采集数据的范围提出了一些限制。确定允许勘探范
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年电磁法勘探-可控源音频大地电磁测深法 2022 电磁 勘探 可控 音频 大地 测深
![提示](https://www.taowenge.com/images/bang_tan.gif)
限制150内