瞬变电磁三维正演的体积分方程法.pdfVIP

瞬变电磁三维正演的体积分方程法 唐新功1’2胡文宝2严良俊2 1北京大学地球物理系北京100871 3 2中国石油天然气集团公司物探重点实验室长江大学研究室湖北荆州43402 由于电磁法时间域正演算法的难度和复杂性，只有极少数简单的三维模型可以求出解析解，大部分三维模型 的计算都需要借助于数值方法。体积分方程法的求解虽然会遇到更困难的数学阿题，但它只需要计算电导率异常 区的贡献，因此，体积分方程法计算瞬变电磁场正演响应是比较合适的。我们使用基于张量格林函数的体积分方 程法进行三维模型瞬变电磁响应的正演研究，首先在频率域内计算电磁场分量的频率域响应，然后利用快速数字 滤波技术将计算结果转换到时间域。计算体积分方程时，需要对异常体进行离散化处理，为此，把异常体划分为 Ⅳ个小元素，每一小元素近似取其中心的值。为了在计算中既考虑精度同时又兼顾速度，离散时采用变步长法， 即在测线附近网格要细一些，远离测线时可以粗一些；异常体埋深浅时，网格要细一些，埋深大时可以粗一些。 计算时，根据实际情况需要，长方体小元素还可以进一步细分为若干立方体的小元素，然后沿这些立方体小元素 对张量格林函数进行积分。张量格林函数在长方体上的积分可以看成是立方体小元素积分的累加。 层状地层中三维异常体中的总电场和总磁场满足如下的离散积分方程； E(r)=Ei,(r)+∑(盯。。一仃，)f?(r；r，，’)E。(r，，) (1) =1 H(r)=Hp(r)+∑(口。一口，)f，(r；rJI’)E。(o) (2) 其中r：：和r?为异常体第n个小元素的张量格林函数，仃“和仃i分别为异常体第n个小元素和第7层的电导。 …’ J盯 率。(1)和(2)式等号右边第一项为背景场，或层状地层的场，可以用快速汉克尔变换法计算；第二项为异常体 的散射场或异常场，用格林函数法求解。 然后使用线性数字滤波法把频率域响应转换到时间域，用数字滤波法近似近计算正弦和余弦变换的富氏变换 法更为合适。使用Anderson(1975，1982)的滞后快速数字滤波技术可以计算这些正弦和余弦变换。用正弦和余弦变 换进行数字滤波的精度与滤波系数、滞后褶积截断误差以及频率域响应的采样及精度有关，丽最终的时间域结果 的精确度还是取决于最初的频率域的响应。数字滤波法一般需要在20～40个频率上(每个数量级中取5～8个点) 计算出三维频率域响应的虚分量，才能得出比较精确的值。计算中为了保证足够的计算精度，在频率域响应的计 算中使用了从0．Ol赫兹到1000赫兹的五个数量级中的总共51个频率点进行计算，即每个数量级中取10个点。 在向时间域的转换中，同样每个数量级中取lO个点，共计算了从0．001秒到lO秒的四个数量级中的共41个时间 点的三维体的电磁场瞬态响应。 分别与半空问均匀模型的解析解和lD模型的数值解进行了比较，验证了软件的准确性。通过多个算例表明， 使用体积分方程法计算瞬变电磁法三维正演响应是有效的方法，正演的结果可以帮助我们深入理解电磁波的传播 规律，为指导野外勘探部署、提高电磁法的解释水平以及未来的瞬变电磁法三维反演打下良好的基础。 ·623·