首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 203 毫秒
1.
时间域航空电磁法激电效应对电磁扩散的影响   总被引:1,自引:0,他引:1       下载免费PDF全文
由于激发极化效应的影响,时间域航空电磁晚期道信号经常会出现变号现象.基于电阻率的传统反演方法无法对变号数据进行正确反演,因此通常在数据处理中予以剔除.为深入了解极化介质的电磁扩散特征,认识航空瞬变电磁负响应的产生机理,本文研究时间域航空电磁系统的电磁扩散特征.我们以均匀极化、非极化半空间及层状介质模型为例,通过直接积分的方法求解频率域电场响应,并由欧姆定律得到电流响应,再经过汉克尔变换得到时间域电流响应.通过研究电流随时间在地下极化介质中的传播特征研究电磁扩散过程;通过对比不同激电参数对电磁扩散的影响,研究极化介质中感应电流与极化电流的扩散规律,从而合理地解释极化介质中负响应的产生机理.基于本文研究和分析结果,可加深对时间域航空电磁法中激电效应的认识.  相似文献   

2.
SOTEM响应特性分析与最佳观测区域研究   总被引:4,自引:2,他引:2       下载免费PDF全文
电性源短偏移距瞬变电磁法(SOTEM)是目前研究和应用较为广泛的一种人工源时间域电磁法工作装置,对深部资源地球物理精细探测具有一定的实际意义.为了深入理解方法内涵并更好地进行推广应用,本文基于电性源瞬变电磁一维正演理论,研究了SOTEM地下感应电流扩散、多分量电磁响应平面分布、多偏移距衰减等特性,然后根据上述特性研究了SOTEM的最佳观测区域.研究结果表明:电性源在地下可以产生水平和垂直两个方向的感应电流.其中,水平感应电流又分为上部水平感应电流和下部水平感应电流(又称作返回电流),水平感应电流的极大值主要集中于发射源附近并垂直向下扩散;垂直感应电流极大值沿与地面呈45°角的方向向下、向外扩散,并且具有较低的振幅和较快的扩散速度.电性源激发的六个方向的电磁场分量都具有一定的探测能力,但是考虑到地面观测的方便性和各分量的传播、分布特点,大多数情况仅利用垂直磁场分量Hz(B/t)和水平电场分量Ex.其中,Hz仅对低阻目标体敏感,且敏感区域位于赤道向区域,并集中在发射源附近;Ex既对低阻体敏感也对高阻体敏感,对低阻体的敏感区域位于赤道向区域,而对高阻体的敏感区域位于轴向区域,并且敏感区域距发射源的距离与目标体埋深和围岩电性有关.  相似文献   

3.
针对岩层中介电常数和磁导率随深度变化的情况,给出水平磁偶极子源电磁场的水平或垂直空间分布特征及相应的变化规律。采用Kong给出的汉克尔(Hankel)J0变换线性滤波器(241点)和汉克尔J1变换线性滤波器(241点)算法,选用均匀大地和两层地电模型,讨论介质的介电常数和磁导率变化时的水平磁偶极子电磁场空间分布特征。计算结果表明,磁场的峨分量幅值与磁导率反相关;磁场H2和电场Ey分量随着磁导率的增大而增大;低频电磁测深中介电常数变化时对电磁场基本没有影响。  相似文献   

4.
各向异性介质模型电性结构复杂,如何进行合理的网格剖分成为获得高精度正演结果的关键,为此本文开展时间域航空电磁各向异性大地三维自适应有限元正演算法研究.通过结合非结构时间域有限元算法和自适应网格优化技术,实现各向异性介质条件下三维时间域航空电磁自适应正演.考虑到时间域航空电磁响应随时间的衰减特性,为了综合评价不同时刻的后验误差,本文将时间作为加权因子,调整各个时刻后验误差的相对权重,进而实现对浅部和深部网格的同步优化.通过与一维解析结果进行对比验证了本文算法的可靠性.数值实验结果显示电导率各向异性对自适应网格影响严重,其最大主轴电导率的数值及其分布特征直接决定了网格加密效果.此外,各向异性对时间域航空电磁三分量响应的分布形态和异常幅值也会产生严重影响,利用全域视电阻率极性图,可以很好地识别各向异性主轴方向.  相似文献   

5.
井中磁源瞬变电磁响应特征研究   总被引:6,自引:1,他引:5       下载免费PDF全文
井中瞬变电磁波勘探是一个全空间地球物理场问题.采用Gaver-Stehfest逆拉氏变换方法,正演计算了瞬变信号激励下接收线圈上的电磁场响应.分析了包含井眼泥浆、套管、水泥环和地层的轴对称多层介质模型的电磁场响应特征,考察了各层介质参数对井中瞬变电磁响应的影响.不同电导率井眼泥浆的电磁场响应衰减曲线表明,井眼泥浆电导率...  相似文献   

6.
时间域航空电磁系统的探测深度是一项关键技术参数,在系统研制设计和资料解释中都有重要作用.当目标层与围岩电性差异不大或目标层引起的异常响应不明显时,常规的估算探测深度方法不再适用,为克服这种缺陷,本文提出了一种系统探测深度的估算方法,该方法通过模拟电场在均匀半空间模型和层状模型中的扩散过程,记录电场最大幅值在地下的瞬时位置,将系统的探测深度定义为观测到的响应值等于给定的噪声水平时,对应时刻的地下介质中感应电场幅度的最大值对应的深度.论文以我国自主研制的直升机时间域航空电磁系统CHTEM-I为例,根据对不同条件下的计算结果的分析,给出了探测深度与模型电导率、飞行高度和噪声水平之间的关系曲线,并总结出一些能够提高系统探测深度的建议.结果表明半空间模型电导率在0.000295~0.0422S/m范围内,CHTEM-I系统可以在噪声水平与t-0.5成正比条件下达到300m的探测深度.本文方法的估算结果不受模型层厚薄,或目标层与围岩电性差异小等因素影响,因而具有较高的实用性.文中的结论对时间域航空电磁系统设计有理论指导意义,也可用于其它时间域电磁勘探系统探测深度的估算.  相似文献   

7.
全波形时间域航空电磁响应三维有限差分数值计算   总被引:11,自引:8,他引:3       下载免费PDF全文
全波形时域航空电磁测量具有解释精度高、分辨率高等优点,已成为新一代时域航空电磁探测发展的趋势.为此,本文基于无源有限差分法,将发射电流离散为若干梯形脉冲元,将脉冲元在时间上产生的电磁场转换为空间初始场,在差分迭代过程中分时引入差分方程,逐步将有源有限差分方程转换为无源有限差分方程,最终实现了发射电流为任意波形时,全波形时间域航空电磁响应的三维数值计算.计算结果表明:在梯形波发射时,对于典型地电模型的航空时间域电磁响应,计算稳定时间大于10 ms;与积分方程法相比,发射线圈中心点二次场响应平均误差小于1%.任意源三维全波形数值计算的实现为全波形三维反演和仪器设计奠定了基础.  相似文献   

8.
时间域航空电磁数据往往在测量过程中受到天然和人文噪声的干扰.如果不能很好滤除这些电磁噪声,那么将会降低资料质量、影响反演的精度,甚至获得错误的解释结果.本文提出了一种基于核主成分分析的去噪方法,通过核主成分分析提取叠加后衰减曲线的主成分,然后使用能量占比方法分离反映地下介质的有效信号和噪声,最后使用反映地下介质的特定成分进行重构.本文所推荐的去噪方法不仅能剔除天然噪声,例如天电产生的尖脉冲或者振荡,而且能有效地抑制人文噪声.分别使用基于核主成分分析的去噪方法,以及AeroTEM软件的处理方法对同样的吊舱式时间域直升机航空电磁勘查系统实测数据进行处理,并比较其结果.处理结果表明:所推荐的去噪方法要优于AeroTEM软件.  相似文献   

9.
用扩散电场法估算CHTEM-I系统的探测深度   总被引:3,自引:2,他引:1       下载免费PDF全文
时间域航空电磁系统的探测深度是一项关键技术参数,在系统研制设计和资料解释中都有重要作用.当目标层与围岩电性差异不大或目标层引起的异常响应不明显时,常规的估算探测深度方法不再适用,为克服这种缺陷,本文提出了一种系统探测深度的估算方法,该方法通过模拟电场在均匀半空间模型和层状模型中的扩散过程,记录电场最大幅值在地下的瞬时位置,将系统的探测深度定义为观测到的响应值等于给定的噪声水平时,对应时刻的地下介质中感应电场幅度的最大值对应的深度.论文以我国自主研制的直升机时间域航空电磁系统CHTEM-I为例,根据对不同条件下的计算结果的分析,给出了探测深度与模型电导率、飞行高度和噪声水平之间的关系曲线,并总结出一些能够提高系统探测深度的建议.结果表明半空间模型电导率在0.000295~0.0422 S/m范围内,CHTEM-I系统可以在噪声水平与t-0.5成正比条件下达到300 m的探测深度.本文方法的估算结果不受模型层厚薄,或目标层与围岩电性差异小等因素影响,因而具有较高的实用性.文中的结论对时间域航空电磁系统设计有理论指导意义,也可用于其它时间域电磁勘探系统探测深度的估算.  相似文献   

10.
回线源瞬变电磁成像的理论分析及数值计算   总被引:9,自引:4,他引:5       下载免费PDF全文
进一步提高瞬变电磁法对地探测的解释精度,提出了回线源瞬变电磁成像原理及数值计算方法. 讨论了频率域中水平层状介质中瞬变电磁响应,得到一个以波阻抗为积分核的双重积分式;然后对水平层介质下电磁场的解进行上、下行波分离,得到含有以反射系数序列为未知的线性方程组,并给出了求取波阻抗和反射系数的数值解法:对实测磁场值进行域的变换,以均匀半空间下的等效波阻抗代替积分核函数,经过线性数字滤波后,在频率域求出等效波阻抗;把频率域中的波阻抗转换到时间域,以此为参数,构建方程组,在时间域用线性规划法求出反射系数序列. 最终以反射系数为参数进行成像. 对理论模型的数值计算结果表明,用本文提出的成像方法可以增强瞬变电磁法识别地下电性分界面的能力.  相似文献   

11.
起伏地表频域/时域航空电磁系统三维正演模拟研究   总被引:1,自引:1,他引:0       下载免费PDF全文
由于航空电磁系统具有工作频率低、时间延迟短等特点,地形对航空电磁响应有很大影响,忽略地形影响会给航空电磁数据解释造成很大误差.本文将基于非结构化网格的矢量有限元法应用于模拟起伏地表条件下频域/时域(FD/TD)三维航空电磁系统响应.该方法由于采用非结构网格,与传统的结构化网格电磁正演算法相比,能更好地拟合地形和地下不规则异常体,提高对不规则地形和地下介质航空电磁响应的计算精度.通过将计算结果与半空间模型的半解析解及已发表的结果进行对比,检验了本文算法的精度.通过对典型山峰和山谷地形航空电磁响应分析对比,总结了地形对航空电磁响应的影响特征.研究结果对航空电磁地形效应的识别和校正具有指导意义.  相似文献   

12.
2.5维起伏地表条件下时间域航空电磁正演模拟   总被引:7,自引:5,他引:2       下载免费PDF全文
时间域航空电磁作为一种高效地球物理勘探技术特别适合我国地形复杂地区(沙漠、高山、湖泊、沼泽等)资源勘查.然而,这些地区地形起伏较大,对航空电磁响应有严重影响,忽略地形影响会给航空电磁数据解释造成很大误差.到目前为止人们对航空电磁地形效应特征研究十分有限.本文提出了基于非结构化网格的有限元法模拟带地形时间域航空电磁系统响应.该方法与基于结构化网格的有限差分相比能更好地模拟地形.首先通过傅里叶变换将2.5维问题转化成二维问题,利用伽辽金方法对二维问题进行离散.通过使用MUMPS求解器,得到波数域电磁响应.利用反傅里叶变换将波数域电磁响应变换到空间域,并利用正弦变换将其变换到时间域,得到2.5维时间域航空电磁响应.通过将本文的计算结果与半空间模型解析解及其他已发表的结果进行对比,检验了本文算法的精度.最后,我们系统分析了山峰和山谷地形对航空响应的影响特征.本文研究结果对航空电磁地形效应的识别和校正具有指导意义.  相似文献   

13.
航空电磁勘查技术发展现状及展望   总被引:11,自引:9,他引:2       下载免费PDF全文
航空电磁作为一种高效的地球物理勘查技术手段,其发展在国外(加拿大、澳大利亚等国家)已十分成熟.然而,在我国该项技术仍处于发展当中,在国内目前尚未形成具有实际探测能力的航空电磁系统和解释手段.这一现状严重影响了我国对地形地质条件复杂区域(比如广大西部地区)矿产资源勘查的需求.本文旨在通过系统介绍航空电磁勘查技术中的基础理论、关键技术、仪器系统、数据处理、解释及应用,并对未来我国航空电磁勘查技术的发展提出建议,使读者了解该技术未来发展方向和研究热点,以期该项技术在我国得到快速发展并获得广泛应用.  相似文献   

14.
传统时间域航空电磁全波形正演模拟主要采用间接法(褶积算法)和直接法(时域有限差分方法等),然而褶积算法需要获得精确的电流二阶导数,这给发射电流数据采集工作带来极大挑战;时域有限差分方法受到网格和时间步长的严格限制,缺乏灵活性.为解决这些问题,本文采用时域有限元方法,通过直接改变每个时间道上的瞬时电流强度模拟任意发射波形的电磁响应.由于无需计算电流二阶导数,大大提高了正演结果的精度.利用基于非结构四面体网格的矢量有限元方法和后推欧拉技术对时间域电场扩散方程进行空间和时间离散,实现三维航空电磁时间域全波形的直接正演模拟.由此不仅可以模拟复杂的地电结构,而且基于后推欧拉法的无条件稳定性,可以更加灵活地选取时间步长,提高计算效率.通过与1D数值模拟结果进行对比验证了该方法的准确性.本文对三维柱状体模型上HELITEM MULTIPULSE和VTEM系统实际发射波形电磁响应进行模拟,并与褶积算法的结果进行比较,验证了本文算法模拟实际发射波形电磁响应的优越性.对复杂三维地质体模型上不同发射波形电磁响应进行模拟,验证了时间域有限元算法可有效处理复杂地下地质结构.  相似文献   

15.
常规的三维时间域航空电磁模拟通常采用隐式步长方法进行时间离散,需要几次矩阵分解和上百次右端源项回带,计算效率较低.为了提高正演计算效率,本文提出使用有理Krylov方法求解时间域电场扩散方程.首先使用非结构四面体网格进行空间离散,采用Nédélec矢量基函数近似四面体单元内的电场;然后基于有限元离散给出矩阵指数和矢量乘积表示的电场显式解;最后采用有理Arnoldi算法构造Krylov子空间内的正交基函数并进一步求解矩阵指数与矢量的乘积,直接得到任意时刻的电场解向量,避免步长离散过程.此外,本文还提出一种指数加权偏移参数优化方法,使得有理Arnoldi近似在瞬变衰减晚期具备更高的精度,从而降低Krylov子空间阶数并提高计算效率.通过和层状模型解析解的对比验证了有理Krylov方法的精度.针对三维异常体模型使用全局网格和局部网格剖分并和其他数值方法比较,进一步说明了有理Krylov方法的有效性.  相似文献   

16.
航空瞬变电磁法对地下典型目标体的探测能力研究   总被引:2,自引:2,他引:0       下载免费PDF全文
航空电磁法的探测能力受飞行高度、发射波形、发射磁矩和发射基频等因素的影响,致使不同分量间的勘探能力存在差异.航空电磁如对所有磁场和磁感应分量、on-和off-time数据进行观测和解释,不仅数据量大、耗时长,而且出现大量冗余数据.目前国内针对此问题尚无系统解决方法.本文针对吊舱式直升机航空电磁系统,采用积分方程法求解频率域响应,经汉克尔变换转换到时间域,计算了地下三维目标体的B和dB/dt时间域响应.利用异常体响应与背景场响应作比值,并通过设定响应阀值定义最大勘探深度,进而分析不同发射波形、不同分量以及on-和off-time期间的航空电磁系统的探测能力.基于本文分析手段,可根据实际勘探目标,确定一套探测能力较强的航空电磁最佳参数组合,为野外测量和数据处理提供技术指导,高效完成勘探任务.  相似文献   

17.
航空电磁系统的收发线圈受外界因素影响会发生姿态变化,使测量结果产生一定偏差.传统的数据处理方法考虑姿态变化情况单一,且主要以频率域系统为主,而目前复杂姿态变化对时间域航空电磁系统全时响应的影响尚未有系统研究.本文在前人研究基础上拓宽思路,不仅研究姿态的角度变化,同时改进前人计算方法,研究线圈发生位置变化的影响,并将线圈与大地的耦合感应效应考虑在内.本文通过定义姿态变化前后两种坐标系,确定姿态变化的角度和旋转矩阵,并将一维层状半解析解与姿态角度变化和位置变化进行整合,推导出收发线圈任意姿态变化的表达式.以直升机吊舱分离装置为例,对时间域航空电磁系统收发线圈任意姿态变化进行细致分析.基于重叠偶极子的假设,给出可应用于实际工作中的姿态校正因子,以提高实测数据的处理效率与精度.  相似文献   

18.
瞬变电磁法的探测深度问题   总被引:9,自引:1,他引:8       下载免费PDF全文
用解析分析、时域有限差分、时-频分析的方法,以地面中心回线装置和阶跃脉冲激励源为例,分析讨论了瞬变电磁测深法的勘探深度问题,以便为野外勘探工作设计提供依据,达到预期的探测目的.解析计算证实了瞬变场在地下以有限速度传播,数值模拟表示出了准静态条件下瞬变场的反射.研究结果表明,由于时间域电磁场遵循因果律,瞬变电磁法的探测深度主要由观测时间决定. 瞬变电磁场的初始传播速度与大地电阻率无关,继后在大地色散作用下,阶跃脉冲前沿逐渐变得平缓,各频率分量的传播速度与电阻率有关,在低阻地层中探测同样的深度需要较长的观测时间. 最大探测深度是在给定时间内电磁波往返地下某一深度的单程距离,最小探测深度受仪器性能的限制,但是埋藏较浅的异常体也有可能在晚时段被观测到.从时-频密度谱中可得到瞬变电磁场信号时间与频率的关系.  相似文献   

19.
Helicopter electromagnetic (HEM) systems are commonly used for conductivity mapping and the data are often interpreted using an isotropic horizontally layered earth model. However, in regions with distinct dipping stratification, it is useful to extend the model to a layered earth with general anisotropy by assigning each layer a symmetrical 3 × 3 resistivity tensor. The electromagnetic (EM) field is represented by two scalar potentials, which describe the poloidal and toroidal parts of the magnetic field. Via a 2D Fourier transform, we obtain two coupled ordinary differential equations in the vertical coordinate. To stabilize the numerical calculation, the wavenumber domain is divided into two parts associated with small and large wavenumbers. The EM field for small wavenumbers is continued from layer to layer with the continuity conditions. For large wavenumbers, the EM field behaves like a DC field and therefore cannot be sensed by airborne EM systems. Thus, the contribution from the large wavenumbers is simply ignored. The magnetic fields are calculated for the vertical coaxial (VCX), horizontal coplanar (HCP) and vertical coplanar (VCP) coil configurations for a helicopter EM system. The apparent resistivities defined from the VCX, VCP and HCP coil responses, when plotted in polar coordinates, clearly identify the principal anisotropic axes of an anisotropic earth. The field example from the Edwards Aquifer recharge area in Texas confirms that the polar plots of the apparent resistivities identify the principal anisotropic axes that coincide well with the direction of the underground structures.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号