基于快速余弦变换的低频电磁场数值滤波法正演.pdfVIP

石油天然气学报 (江汉石油学院学报) 2011年 1月 第33卷 第 1期 JournalofOilandGasTechnology(J．JPI) Jan．2011Vo1．33 No．1 基于快速余弦变换的低频 电磁场数值滤波法正演 徐振平(芸 ：。。) 胡文宝 (油气资源与勘查技术教育部重点实验室 (长江大学)，湖北荆州434023) [摘要]利用快速汉克尔变换理论 的推导方法 ，对快速余弦变换 的滤波系数及算法进行 了推导，得 出了快 速余弦变换的滤波系数。使用该系数计算含余弦的高振荡积分 ，降低 了运算的时间，提高 了计算的精度 ； 并在柱状分层 的地层 电磁场 响应的正演分析 中进行 了应用，分析 了快速余弦变换滤波算法在低频 电磁 响 应 的误差 。 [关键词]快速余弦变换 ；快速汉克尔变换 ；滤波系数 ；柱状分层 电磁响应 [中图分类号]P631．322 [文献标识码]A [文章编号]1000—9752 (2011)01—0063—05 在很多应用领域 ，高振荡积分问题是一个非常重要的研究课题n]。当振荡频率远远超过积分结点的 数 目时 ，经典积分方法如梯形公式法、Simpson法 、高斯积分法等对应用广泛的高振荡积分的数值计算 将失去效用 ；傅里叶变换L2或快速傅里叶变换 (FFT)存在积分变量较大时发生畸变的问题 ；Filon方 法、Filon型方法_3j、Levin方法 、Levin型方法 ]、渐 近法 、广义积分方法 以及数值最速下 降法 都 存在处理复杂 、计算精度不够的问题 。当振荡部分为正余弦时，可以采用滤波计算方法 ，利用快速汉克 尔变换系数以及汉克尔变换和余弦变换 的关系得出快速余弦变换系数 ，并用该系数计算高振荡积分口]。 所有这些方法都有一定的应用条件和精度 问题。应用于地球物理频率域 电磁响应计算时 ，由于应用 条件复杂 ，地电模型和电磁场量不同等 ，其响应函数的特性差异很大。为此 ，此次研究基于滤波算法的 数值积分法 ，利用复数的围道积分方法，并引入修正函数以改变系数计算中级数项的收敛性，按照不同 极点分 区段用有限个级数项的和替代无穷求和 ，最后得 出快速余弦变换滤波系数 ，再用该系数计算高振 荡积分 ，收敛快、精度高、运算速度快。 1 柱坐标下低频 电磁响应数学模型 图1是柱坐标下 3层介质 电磁响应的模型，p 和P。分别 为第 1层、第 2层介质的半径 ， 为3层介质的相对磁导率 ， 。 为 3层介质的电导率 ，发射线圈为 T，接收线圈为 R。此 时 ，电场只有 分量 E ，磁场有P分量 H 和 z分量H 。其中 R E 在无源区满足的方程组为： 等a+。鲁a+。吉aP+。(＼k。一fD)，E一≠一0㈩… 式中，k 一 鲫 一i ； 为角频率 ，rad／s； 为介质 的磁导 丁 率 ，H／m；d为 介 质 的 电导率 ，S／m；￡为 介 质 的介 电常 数 ， F／m 。 图 1 三层柱状分层介质模型 在利用边界条件 和分离变量法后加上初始激励 ，并按照 [收稿 日期]2010—09—17 [基金项 目]国家 973 规划项 目 (20(】7cB2096o7)；国家自然科学基金项目 44704073)。 [作者简介]徐振平 (1968一)，男 ，1990年大连理工大学毕业 ，副教授 ，硕士生导师，现主要从事地球物理勘探方面的科研工作。 石油天然气学报(江汉石油学院学报) 2011年 1月 H一一土f+ 1，可以得到如下的积分表达式8[9]： 1“ 、p “p H一一 r厂()COS(1z)d (2) 式 中，M 为发射源的磁偶极矩，Wb ·1TI；_厂()是含修正 的Bessel函数的表达式 ，具体见下式；