航空瞬变电磁一维正演的连分式算法的研究.pdfVIP

第２８卷 第４期 湖南科技大学学报（自然科学版） Ｖｏｌ．２８Ｎｏ．４ 　　 　 ２０１３年　 １２月 ＪｏｕｒｎａｌｏｆＨｕｎａｎＵｎｉｖｅｒｓｉｔｙｏｆＳｃｉｅｎｃｅ＆Ｔｅｃｈｎｏｌｏｇｙ（ＮａｔｕｒａｌＳｃｉｅｎｃｅＥｄｉｔｉｏｎ） Ｄｅｃ．　２０１３ ① 航空瞬变电磁一维正演的连分式算法研究 龙剑波，强建科 （中南大学 地球科学与信息物理学院，有色金属成矿预测教育部重点实验室，湖南 长沙４１００８３） 摘　要：针对快速Ｈａｎｋｅｌ变换精度不高的问题，改进了连分式算法，使之能够计算余弦变换，并与快速汉克尔变换算法作 了比较，结果表明：无论是计算积分收敛或核函数快速衰减的余弦变换，还是计算核函数震荡增加的发散型Ｈａｎｋｅｌ积分时， 连分式算法都具有精度高、计算稳定的特点，而滤波法计算核函数震荡增加的发散型Ｈａｎｋｅｌ积分的误差较大．最后，把连分式 算法应用于航空瞬变电磁一维正演模拟计算，得到了满意的瞬变响应，其计算精度、速度和稳定性很好，为瞬变电磁模拟计算 提供了一种新的计算方法． 关键词：余弦变换；连分式；数值计算；瞬变电磁场 中图分类号：Ｐ６３１．３２５Ｐ６３１．３２６　　　文献标识码：Ａ　　　文章编号：１６７２－９１０２（２０１３）０４－００８６－０６ ［５］ ［６］ 　　在瞬变电磁场理论计算中，磁场值是与观测参数 法 以及数字滤波法 等．Ｆｉｌｏｎ方法是一种较早的 感应电动势密切相关的物理量，磁场分量的高精度数 直接数值求积方法，一般分段越多，精度越高，在达 值求解是瞬变电磁法正演研究和反演解释的重要一 到相同精度的计算时，计算量随着ｋ的增大会明显 环，其在一维理论计算式中往往表现为Ｈａｎｋｅｌ积分 增加，目前已经较少直接应用；折线逼近法利用分段 或变换形式［１］．以中心回线瞬变电磁２．５维有限元数 插值函数近似原函数，然后在每一分段内进行数值 值模拟为例，矩形线源在背景模型上产生的垂直磁场 求积，理论上逼近精度可以足够高，但逼近步长很难 分量Ｈ求解中，目前应用较多的一个方法是将其化 控制，须视核函数的实际情况处理．现在地球物理电 ｚ 为拉氏傅氏域下的一个余弦变换，对该余弦变换进行 磁法中常用的是基于线性卷积理论的数字滤波法 求解之后，需进行反傅氏变换（可转换为余弦变换）， （ＤｉｇｉｔａｌＦｉｌｔｅｒ），又称为快速汉克尔变换 （Ｆａｓｔ 完成之后可用数值微分求解感应电动势．在以上过程 ＨａｎｋｅｌＴｒａｎｓｆｏｒｍ），该法首先由Ｇｈｏｓｈ于１９７１年引 ［７］ 中，计算的积累误差会对最终结果产生严重影响，特 入地球物理电法勘探正演计算 ，之后得到了许多 别是２次余弦变换的误差不容忽视． 改进和应用，并已经成为地球物理中汉克尔变换计 余弦变换是Ｈａｎｋｅｌ变换的一种特例，属于高震