全空间水平层状介质垂直磁偶源的瞬变响应1【精品论文大全】 .doc

精品论文推荐全空间水平层状介质垂直磁偶源的瞬变响应1胡博岳建华杨海燕（中国矿业大学资源与地球科学学院徐州221008）email: 摘要：全空间瞬变电磁场的分布规律是亟待解决的矿井瞬变电磁法的理论基础。本文详细介绍了全空间水平层状介质垂直磁偶源瞬变响应的计算方法，并考虑了巷道存在对场的影响。 建立了全空间层状介质模型；从电磁场矢量位出发，根据边界条件得到一组方程，解之得到 各层的电磁场分量表达式；将表达式离散化，应用汉克尔变换滤波算法求解；并用余弦变换 多项式法将频率域的结果变换到时间域。用频率域的解析解验证了该公式的正确性，并对比 了巷道有无时场的变化情况，说明了巷道的影响是客观存在的，应在 3d 模型空间中进一步 研究。该算法为巷道全空间瞬变电磁场的分布特性提供了理论基础。关键词：全空间 磁偶源 汉克尔变换 巷道 中图分类号：p631. 31、引言矿井物探固有的特殊理论问题是全空间地球物理场问题，巷道对场的影响一直是复杂的 难以解决的难题。目前应用较为广泛的是矿井瞬变电磁法，它借鉴的是地面电磁法的基本理 论1-2，尚未形成自己的理论体系。弄清电磁场的全空间分布规律和巷道的影响是真实认识 井下地电响应特征的基础。真正在全空间条件下研究巷道对电磁场影响的很少，stefi krivochieva 和 michel chouteau3给出了层状介质中半空间和全空间水平放置的方形回线 产生电磁场的数值模拟，采用的是中心回线和分离回线装置，发射线圈在地面上时，是水平 的不接地线圈，正方形或矩形；在井下时，发射线圈在形状和大小上有一定变化，取决于巷 道支护间可以通过的路径；接收线圈是小的单匝或多匝线圈。由于受井下施工环境的限制，采用磁偶极子源作为发射源更为便利，因此研究巷道全空 间条件下磁偶源的响应是首要解决的问题。本文以垂直磁偶源分离回线装置为例，详细介绍 了垂直磁偶源瞬变电磁场分量的计算方法，包括水平层状介质模型的建立，公式推导，汉克 尔变换滤波算法和时频转换等一系列问题，并和频率域均匀半空间的解析解计算的结果相对 比，验证了公式的正确性。对比了巷道存在时和无巷道时场分量的变化情况。2、计算方法1本课题得到国家自然基金项目资助，项目编号：406740748建立全空间水平层状介质模型（n 层），并建立柱坐标系，设各层的参数为s j , hj （电导率和层厚），源位于第 i 和 i + 1 层交界处，如图 1 所示。矢量位 a* 在每一层介质中都满足 亥姆霍兹方程，用下标 j 来表示层位，则有j j j2 a* - k 2 a* = 0z由于 a* 仅有 z 分量 a* ，于是zj z2 a*j- k 2 a* = 0j其通解为a* = iwm0 pm (c eu j z + d e-u j z )j(lr ) d lz j 4p 0 j j 0式中， pm = si 为偶极距， c j , dj 为待求系数，2 2 2u j =l + k j , k j= -iws j m0 , j = 1,l, n 。 边图 1 全空间水平层状介质模型界条件为：电磁场切向分量在各个界面上连续，在各层分界面上，下列等式成立*a = a ，*zjzj +1azjz= azj +1z这里 j = 1, 2,l,n 。 z - 时， d1 = 0 ， z 时， cn = 0 。设从 - 到 各层的界面mm m到 坐 标 原 点 的 垂 直 距 离 分 别 为 h1 , h2 ,l, hn-1 ， 其 中 h1 , hn = 。 设 各 层 的 厚 度 为h2 , h3 ,l, hn-1 。根据边界条件可建立一系列方程，解之得到如下计算步骤：（1） z j = -iwm / u j , j = 1, 2,l, n ； z(1)= z1 ， z(n)= zn ；m = 2,l, i -1， z (m) = zz (m-1) + z th (u h )；m(m-1)zth (u h) + zk = i +1,l, n -1， z (k ) = zm mmkk kz (k +1) + z th (u h )k(k +1)zth (u h) + z（2）求源所在层的系数k kk 。l z (i +1) - zc = i ，iu z (i -1)th (u h) + ziiz + z (i +1)ii -1iiii -1z (i -1) + z th (u h ) l z (i -1) - z z (i +1) - zd = i 1 + i e-2ui hi-1iu z (i -1) + z z (i -1)th (u h) + z ii z + z (i +1) =ii -1i iii -1) iz (i -1) + z th (u h（3）求其它层的系数iic e-ui hi-1 + d eui hi-1 + l e-ui hi-1(i -2)ci -1 =1 + z(i -2)ui- zi -1 e-2ui-1hi-1e ， dui-1hi-1i -1z+ z= z (i -2)- zi -1 -2ui-1hi-2e ci -1i -1z (i -2) + zi -1c e + d e-uk +1h k uk +1h k k +1 k +1 uk h kz - z(k -1) k -2uk hk -1k = i - 2,l, 2 时， ck =z ( k -1) - z1 + k kz (k -1) + ze e-2uk hk l， dk = (k -1)z + zke ckc = 2z(1)e(u1 -u2 )h1 c ，cci + di + u= i ，dz (i +2) + z= i +1 e2ui+1hi+1 c12z + z (1)2i +1z (i + 2) + z1 + i +1 e-2ui+1hi+1z (i + 2) - zi +1z (i +2) - zi +1i +1m = i + 2,l, n -1时， cm =i +1c e + d eum-1h m-1 -um-1h m-1 m-1 m-1 z ( m+1) + ze-um h m-1 ，1 + m mz ( m+1) - ze2um hm z (m+1) + zd = m e2um hm c2z (n)d =e(un-1 +un )hn-1 cmmz (m+1) - zmn z (n) - za；n-1n-1*。2（4）利用快速汉克尔变换计算频率域电磁场。仍以上面推导的源所在水平的 h z 为例，iwm h= -k 2 a* + zjz z z 2 h = pm(u2 - k 2 ) c+ d + l (lr ) d lz 4p 0 i i iu 0 i 根据快速汉克尔变换理论，得到快速汉克尔变换滤波算法如下：g (r ) = 1v f (nd ) h * (m - n ) dr n=-*其中 g ( r ) 是输出函数，f (nd ) 是输入函数，r 是收发距，hv (m - n) d 是快速汉克尔变 换滤波系数，其计算方法及值可参考相关文献3-5，文中计算所用的滤波系数是文献4中给出的 61 点 j 0 系数。 h z 可写为h = pm 1 f + f+ f + f h n zz (i +1) - z (r4p n=-1234 ) 0 令 f =i，则-(i 1) ( )z + z (i +1) zth u hi 1+ zii-iii -1iz (i -1) + z th (u h)l 3l 3 z (i -1) - z l 3f1 =uif ， f = i i i2u z (i -1) + ze-2ui hi-1， f3 = ff2 ， f4 = 。ui（5）将频率域的场变换到时间域。用余弦变换多项式法7-8计算和早期的，( ) 2m0n im h (w j )w j - im h (w j +1 )w j +1 cos ()cos ()p t 2b t =j =0w j +1- w j w j t -w j +1t b(t) 2m0n re h (w j ) - re h (w j +1 )cos ()cos ()p t 2t=j =0w j +1- w j w j t -w j +1t 晚期的用差分算法计算9。同时计算了定义的晚期视电阻率。3、实例和讨论本文以垂直磁偶极子所在水平的频率域的垂向磁场 h z 为例进行了模型试算。 首先，并用频率域均匀半空间水平大地的解析解来验证算法的正确性。图 2 是在相同的条件下（磁偶源置于电导率为 0.0 1sm 的大地表面，观测点距磁偶源 100m）用上述推导的全空间层状介质模型的公式和均匀半空间解析式计算的结果8的对比曲线（横轴为频率（单位 hz），纵轴为 h z （单位 a/m），圈为全空间层状介质公式计算结果，实线为解析解。 由该对比曲线可以看出，两种计算方法一致性很好，证明了推导的公式是正确的。其次，计算时间域均匀全空间（磁偶源置于电导率为 0.0 1sm 的全空间中，观测点 距磁偶源 100m）和半空间（磁偶源置于电导率为 0.0 1sm 的大地表面，观测点距磁偶源100m）和晚期视电阻率的对比曲线，虚线表示半空间的值，圈表示全空间的值，如图3 所示。可见，半空间场的理论不能很好的适用于全空间场。图 4 给出一个七层模型的算例， 虚线表示负值，实线表示正值。 各层参数为r1 = 1000, h1 = 100,r2 = 50, h2 = 200,r3 = 100, h3 = 50,r4 = 1, h4 = 20,r5 = 50, h5 = 100, r6 = 100, h6 = 50, r7 = 1000, h7 = 100 （取标准单位），收发距为 20 米。图 2 h z 实部和虚部对比曲线图 3 均匀全空间和半空间对比曲线最后考虑巷道对瞬变场的影响，将源所在的层面看作是巷道，建立一个上下对称的七层 模型，第四层模拟巷道，巷道存在时 r4 = 10000 ，巷道不存在时 r4 = 100 ，巷道层的高度 为 4m ，其它各层参数为：r1 = 50, h1 = 1000, r2 = 1, h2 = 50, r3 = 50, h3 = 100, r5 = 50, h5 = 100, r6 = 1, h6 = 50,r7 = 50, h7 = 1000 （取标准单位），收发距为 20 米。分别计算有巷道和无巷道的情况进行对比，如图 5 所示，虚线表示巷道存在的结果，圈表示无巷道的结果。说明巷道影响是存在 的，由于研究的是 1d 模型，反映出的影响不明显，若建立 3d 空间模型，可得到更准确的 结果。4、结论图 4 七层模型算例本文推导的公式适用于全空间水平层状介质垂直磁偶源瞬变响应的计算，可用均匀半空间的解析解验证其正确性。根据推导结果可进一步计算任意水平和深度位置的三个场分量， 为获得全空间瞬变电磁场的分布特性奠定了基础。对比了均匀全空间和半空间及巷道有无时 曲线的变化情况，说明了全空间场具有自己的特性，并且巷道对电磁场的影响是存在的。进 一步的研究任务是，结合井下施工技术参数，在一定的约束条件下进行大量计算，得到较为准确的巷道影响因子，为矿井瞬变电磁法的正反演提供理论依据。图 5 巷道有无的对比曲线参考文献：1 李貅瞬变电磁测深的理论与应用西安：陕西科学技术出版社，2002 2 朴化荣电磁测深法原理北京：地质出版社，19903 stefi krivochieva and michel chouteau. whole space modeling of a layered earth in time-domain electromagnetic measurements. journal of applied geophysics, 2002,vol.50:375-391.4 adel a. mohsen, e.a. hashish. the fast hankel transform. geophysical prospecting, 1994, 42, 131-1395 guptasarma, b. singh. new digital linear filters for hankel j0 and j1 transforms. geophysical prospecting,1997, 45, 745-7626 阮百尧均匀水平大地上频率域垂直磁偶源电磁场数值滤波解法桂林工学院学报，2005，25（1）7 李金铭，罗延钟电法勘探新进展北京：地质出版社，19968 殷长春，李吉松任意角度瞬变电磁测深正演问题计算方法物探与化探，1994，18（4）9 米萨克 n.纳比吉安主编,赵经祥等译.电磁法（第一卷 理论）.北京：地质出版社,1992.transient response of vertical magnetic bipolar source in full space plane layer mediumhu bo, yue jian-hua, yang hai-yan(school of resource and geo-science, china university of mining and technology，xuzhou, jiangsu,221008, china)abstractdistributed features of full space transient electromagnetic fields are the unsolved academic base of tem in mine. the paper introduces a method of calculating the transient response of vertical magnetic bipolar source in full space plane layer medium in detail, and the effect of laneway is considered. a full space plane layer medium model wa