位场向下延拓新方法在海磁数据融合中的应用.pdfVIP

位场向下延拓新方法在海磁数据融合中的应用 肖鹏飞1，徐世浙2，孟令顺1 (1．吉林大学地探学院，长春，130026；2．浙江大学地球科学系，杭州，310027) 1．引言 航磁数据是在高空测定的，且往往测量面起伏，数据量巨大。要把它与海平面测量的海磁数据融合，建立 三维空间的磁场数据库，就必须找到一种计算速度快、精度高的向下延拓方法。对于平面延拓，最常用的方法 是快速傅立叶变换(FFT)，但FFT主要用于向上延拓。在向下延拓时，当延拓深度较大，FFT会出现强烈振荡， 延拓效果极差。航磁和海磁的融合中，海磁向上延拓出三维空间磁场数据库是属于平面向上延拓问题，可以用 FFT较好地解决。但海磁的测量范围小，航磁的测量范围大，要建立三维空间磁场数据库，主要依靠航磁。航磁 向上、向下延拓出三维空间磁场数据库是属于曲面向上、向下延拓问题，不能用FFT计算。目前，国内外最常 用的曲面向上延拓方法是等效源法(候重初，1984，Ivan，1994，等等)，此方法需解线形代数方程组，要求大 量的计算时间。至于曲面向下延拓问题，当延拓深度较大时(几百米)，等效源也是很不稳定的。所以说，到目 前为止，还没有一种适合于海磁和航磁融合、建立兰维空间磁场数据库的高效，又有一定精度的资料处理方法。 我们采用的方法是将起伏面上的位场观测值放置在要向下延拓的水平面的垂直投影点上，把它作为此水平面上 的位场初值。然后用积分的方法，从水平面上的位场初值计算出观测面上的位值。用观测值与计算值之差，对 水平面上的位值进行改正。如此迭代进行，获得水平面上的位值。 2．技术方法 首先将磁场观测面上f点的磁场甜。，，放置在欲向下延拓磁场的水平面的垂直投影点f上，作为水平面上的 磁场的初值。第二步，当观测面与水平蘧间没有场源时，磁场满足拉普拉斯方程。此时，可用下面的积分公式， 由水平面上的磁场的初值“：’计算观测面上f点的磁场．甜_}1)ln：f掣罢掣扩，其中’r是观测面 f —n 么刀，’ 上的f点至水平面上的点的矢量，n是水平面的法向，COS(?I，n)是r与n夹角的余弦．由此式计算的观测面 上的磁场，用“g’表示。然后，用已知的观测面上的磁场与”：’的差值，对“：’进行校正，得到新的甜竿’。重 “。≈”：”“．。一般迭代20次即可。在延拓出水平线上的磁场后，可用水平线磁场向上的延拓公式，延拓出高 于该水平线的其它水平线上的磁异常。可以看出，这一新方法的原理并不复杂，但用它实现位场的快速、准确 的延拓，需要采取许多技术措旋，如异常范围的选取，延拓点的距离,3D延拓计算速度的提高等等，特别是如何 编制实用化、可视化的程序，还要做很大的努力。 3．应用 在某地由航磁、海磁获得的资料，我们正在建立密度为lOm*lOm*Sm的三维背景磁场的数据库。本方法所用 延拓的精度是用航磁向下延拓至海平瓶与海磁测量结果的对比，或海磁向上延拓至航磁高度与航磁测量结果的 对比来确定的。在航磁和海磁位置准确的条件下，精度达到磁异常强度的7-10％。在计算速度上用二维位场延 拓方法，对剖面数据进行延拓时，用!oo，000数据点延拓出100，000数据点的时间小于60s：用三维位场延 拓方法，对平面数据进行延拓时，据点延拓据点的时间小于5小时。 参考文献 1． 徐世浙等，起伏地形位场快速延拓新方法，高校地质学报，VoI．8，No．3，Sept．2002， 2． 侯重初，蔡宗熹，刘奎俊，剖面曲线上的位场转换系统，地球物理学报，27(6)：573．581，1984 3． Ivan continuationof fieldfro,na M．Upward potential polyhedral 1994，42(3)：391—404 · ·． 331 位场向下延拓新方法在海磁数据融合中的应用 作者： 肖鹏飞， 徐世浙， 孟令顺 作者单位： 肖鹏飞,孟令顺(吉林大学地探学院,长春,130026)， 徐世浙(浙江大学地球科学系,杭州