一种改进的基于导数迭代的磁场向下延拓方法

高精度的地磁场数据库是地磁匹配导航的基础.在构建地磁数据库时,需要将空中的航磁测量数据向下延拓至某一水平面,而向下延拓在数学上具有不适定性,通常利用迭代方法逐次逼近"适定最优解".本文针对导数迭代法在迭代过程中的噪声累积问题,提出一种改进方法,首先利用稳定的向上延拓将观测面和向上延拓、向下延拓同等高度平面上的磁场值联系...     

位场向下延拓的正则化方法

本文讨论了正则化方法位场向下延拓的四个频率响应公式,分析了它们的特点。这些方法大都能有效地压制由于观测误差或高频干扰所引起的振荡,使下延结果趋于稳定。给出了各种试验的结果。     

位场向下延拓的正则化方法

本文讨论了正则化方法位场向下延拓的四个频率响应公式,分析了它们的特点。这些方法大都能有效地压制由于观测误差或高频干扰所引起的振荡,使下延结果趋于稳定。给出了各种试验的结果。     

位场向下延拓的改进迭代维纳滤波法

根据维纳滤波理论导出的位场向下延拓滤波器为最佳下延滤波器,但因其实现需要已知待求位场和噪声的功率谱而在实际应用中受到限制.针对该问题,本文首先提出一种基于位场径向平均功率谱的位场噪声水平估计方法,进而利用偏差准则求取正则化参数,实现位场正则化向下延拓;然后将位场正则化下延结果的功率谱作为待求位场功率谱的估计初值,采用带修正项的迭代维纳滤波方法来更新对待求位场功率谱的估计,最后提出本文的位场向下延拓改进迭代维纳滤波方法.基于理论重力模型数据及航磁实测数据进行了向下延拓对比试验,结果表明,改进迭代法具有较好的收敛性,且下延精度优于Tikhonov正则化法和递增型维纳滤波法.     

位场向下延拓三种迭代方法之比较

位场向下延拓在重磁资料解释和用于位场导航的基准数据库构建中发挥着重要作用.本文针对第一类Fredholm积分方程的三种空间域迭代解法:迭代Tikhonov正则化法、Landweber正则化迭代法和积分选代法,基于算子理论和不适定问题的正则化处理方法,首先利用傅里叶变换将空间域迭代法变换到波数域,然后由数学归纳法推导得到这三种迭代法对应的波数域位场向下延拓算子;由Landweber迭代法和积分迭代法在迭代形式上的相似性,探讨了它们在位场向下延拓中的异同及各自优势.模型对比分析表明:(1)两种迭代正则化方法在正则化参数选择合适的条件下,其向下延拓的效果要明显优于积分迭代法,且当收敛到相同误差水平时,迭代Tikhonov正则化法在迭代次数上要远远小于Landweber选代法,但迭代Tikhonov正则化方法存在对正则化参数敏感的问题;(2)从实际应用上讲,由于积分迭代法不存在正则化参数的选择问题,所以该迭代法具有较强的实用性,但需考虑其波数域向下延拓算子时噪声的放大效应.     

位场向下延拓的迭代最小二乘法

将位场向下延拓作为一种反问题,结合正则化方法和快速傅里叶算法的优点,本文提出了一种解决位场向下延拓问题的迭代最小二乘法.新方法将位场向下延拓问题转化为最小二秉问题,采用梯度下降法对最小二秉问题进行迭代求解,在迭代过程中,利用快速傅里叶变换方法解决了计算量大的问题.新方法对延拓数据中的噪声具有很好的抑制作用.利用不合噪声...     

正则化等效层重力向下延拓方法

近年来,利用时移微重力技术进行储层开发监测受到国内外学者广泛关注.时移微重力观测数据存在信噪比低,信号弱的问题,难以实现储层内物质运移的定量解释.为压制数据噪声,增强有效弱信号,本文研究了利用Tikhonov正则化方法反演等效层（源）,并由等效源实现重力场向下延拓的方法;在此基础上,本文推导了波数域正则化等效源向下延拓算子.针对向下延拓场幅值衰减问题,提出了正则化等效源迭代补偿算法.通过模拟数据实验研究了不同深度正则化等效源滤波算子及向下延拓算子的波数响应;与波数域Tikhonov正则化向下延拓方法相比,正则化等效源向下延拓方法的延拓精度更高、更稳定.最后,将基于迭代补偿的正则化等效源向下延拓技术应用于实测时移微重力数据证实了该方法能够有效增强局部异常,实现时移微重力数据大深度稳定向下延拓.

     

噪声对磁场向下延拓迭代法的计算误差影响分析

磁场向下延拓可以增强数据的分辨率.在现有的向下延拓方法中,迭代法具有较好的计算效果和较大的向下延拓高度.首先介绍了迭代法的基本原理,分析迭代法的特点;然后利用迭代法对不同噪声水平的模型数据和实测数据进行向下延拓计算,研究噪声对迭代法的计算误差影响.通过延拓结果和实际数据比较分析,客观评价了迭代法的计算精度、可靠性及适应性,对推广应用该方法具有实质性意义.     

迭代法与FFT法位场向下延拓效果的比较

将水平观测面上的实测位场值,垂直投影至下部的延拓水平面上,作为该水平面上的位场初始值. 根据该水平面上的初始值,用快速傅里叶变换(FFT)的方法向上延拓计算观测面上的位场值. 用观测面上的实测值与计算值的差值,对延拓面上的位场值进行校正. 如此反复迭代,直至观测面上的实测值与计算值的差值小到可以忽略. 这种空间域的迭代法原理简单,不用解线性代数方程组,有较高的计算速度和良好的延拓效果. 本文用迭代法对模型数据和实际数据进行向下延拓,对比了迭代法与常规的FFT法在位场向下延拓中的效果,迭代法显著优于FFT法.     

位场向下延拓的波数域迭代法及其收敛性

提出了位场向下延拓的波数域迭代法. 对水平面上的位场观测值进行Fourier变换,得到其波谱. 根据第一类Fredholm积分方程的空间域迭代解法,推导出计算向下延拓水平面上位场波谱的波数域迭代公式. 在波数域中进行迭代,一直进行到相继两次迭代近似解的差值最大绝对值小于给定的精度,或迭代达到给定的最大迭代次数. 对这种迭代近似解进行Fourier逆变换,得到向下延拓的位场. 数值计算结果表明：与空间域迭代法比较,这种波数域迭代法简单、快速,并有同样好的向下延拓效果. 本文还证明了这种迭代法是收敛的,并给出了它的收敛特性和滤波特性.     

地磁正常场的选取与地磁异常场的计算

根据2003年中国地磁观测数据（包括135年地磁测点和35个地磁台）以及我国邻近地区38个IGRF计算点的地磁数据,计算中国地磁异常场的分布。选取两种地磁场模型作为地磁正常场,一是国际参考地磁场的球谐模型,二是中国地磁场泰勒多项式模型。根据各个测点的地磁异常值（观测值减去模型计算值）,用球冠谐分析方法计算地磁异常场的球冠谐模型,并绘制2003年中国地磁异常（△D,△I,△F,△X,△Y,△Z）。分析和讨论了中国地磁异常场。     

补偿向下延拓方法研究及应用

位场向下延拓是重、磁处理和解释的常用方法,但其不稳定性限制了其在资料处理及反演中的应用.本文基于补偿圆滑滤波思想以及空间域向下延拓迭代法,通过逐次补偿的办法实现位场的稳定向下延拓.同时,在频率域空间给出了该下延方法的频率域响应因子,并讨论了其低通滤波特性,理论模型和实际位场资料试验表明该方法向下延拓稳定性具有较高的延拓精度.将其应用于重力密度界面反演中,改进反演的稳定性,实际莫霍界面反演表明下延因子具备实用性.     

补偿向下延拓方法研究及应用

位场向下延拓是重、磁处理和解释的常用方法,但其不稳定性限制了其在资料处理及反演中的应用.本文基于补偿圆滑滤波思想以及空间域向下延拓迭代法,通过逐次补偿的办法实现位场的稳定向下延拓.同时,在频率域空间给出了该下延方法的频率域响应因子,并讨论了其低通滤波特性,理论模型和实际位场资料试验表明该方法向下延拓稳定性具有较高的延拓精度.将其应用于重力密度界面反演中,改进反演的稳定性,实际莫霍界面反演表明下延因子具备实用性.     

约束等效源稳定向下延拓方法研究

在等效源向下延拓方法的理论基础上,研究提出加入观测平面以下空间实测数据作为约束的等效源约束向下延拓法,用于二维位场数据向下延拓.与传统的向下延拓方法和等效源法延拓相比,等效源约束向下延拓法抗噪声性强,能延拓出由于延拓深度过大而衰减至难以探测的高频信号,使延拓结果更接近真实值.理论重力模型数据的延拓结果显示该方法有效,延拓结果优于传统波数域迭代法和等效源法.利用等效源约束延拓法对我国某海域水面实测重力值进行向下延拓并与水下实测值进行对比,结果表明该方法延拓结果最接近实测值.本文方法为我国深水油气近海底重磁高精度探测研究提供一定的技术支撑.

     

非偶极子磁场西向漂移的频散特征

根据第8代国际参考地磁场的资料,本文采用Briggs方法对1900～2000年期间m=1～10次谐波的磁位的全球平均西向漂移速度进行了计算和分析,得出了与Hide的理论预测完全不同的频散特征,即地磁场西向漂移呈现出负频散特征.本文还就Briggs方法得到的结果与采用Malin纬度剖面法的结果进行了比较和讨论,结果表明,正是由于非偶极子磁场的某些谐波分量及其漂移速度的全球分布存在着纬度依赖性,导致了两种方法计算出的个别谐波的漂移速度存在着较大的差异.     

2000年中国地区地磁场长期变化的区域特征

本文用中国地磁场复测点和地磁台站资料,建立了2000年代中国参考地磁场长期变化模型CGRF-SV2000.模型显示,中国地区地磁场变化比较平缓,X、Y、Z、H、D、I、F七个地磁要素的“无符号平均年变率”分别为12.2 nT/a、8.2 nT/a、43.8 nT/a、11.8 nT/a、0.96(′)/a、2.99(′)/a、22.4 nT/a,比国际参考地磁场IGRF给出的全球年变率约小1/3到1/2.各地磁要素的变化显示,地磁北极正在向中国移近,或者说,中国正在向高地磁纬度方向移动,平均移动速度约为3(′)/a.磁偏角变化还显示,中国地区东西部偏角差异继续扩大.作为检验和对比研究,本文利用第8代国际参考地磁场（IGRF）模型,分析了全球地磁场长期变化的时空特征,讨论了全球长期变模型IGRF_SV与中国长期变模型CGRF_SV的异同点.对比分析表明,中国地区地磁场的长期变化与全球长期变化总趋势基本符合,但是,CGRF_SV也表现出一些特有的局部异常特征.