位场数据曲化平的迭代法

位场数据曲化平是位场数据处理解释中的重要运算,但是它的计算量和计算的复杂性影响了它在许多处理和解释方法技术中的应用.本文提出一种位场数据曲化平的迭代方法,即通过把位场数据曲化平视为平面位场数据向上延拓的反问题,得到曲化平的线性积分方程,再把曲面上位场数据视为曲面平均高程面上的位场数据,利用向下延拓的波数域广义逆算法把平均高程面上的位场数据向下延拓到设定平面上,再根据曲面和其平均高程面的相对起伏对设定平面上的向下延拓数据进行起伏校正,最后再把所得平面上的位场数据向上延拓得到曲面上的位场数据,并进行迭代.把这种方法用于三维理论模型数据和实际磁场数据的曲化平处理均获得了理想的结果.     

位场曲化平的插值-迭代法

将起伏曲面B上的位场向下延拓至曲面最低点的平面A的插值-迭代法步骤是：1）将曲面B上的场值放置在水平面A上具有相同水平坐标的点上,作为A上的初值;2）用若干水平面切割B,从A的初值,用快速傅里叶变换法(FFT)向上延拓出这些平面的场值,用插值的方法从这些平面的场值计算曲面B的场值;3）根据B上的实测值与计算值的差值,对A上的值进行加权改正;4）重复步骤2）和3）,直到B上的差值小到可以忽略.这种插值-迭代法具有高的计算速度,比通常的FFT法延拓得更深,可以超过10倍点距.文中给出计算实例.     

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

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

位场延拓的积分-迭代法

本文介绍一种新的位场延拓方法——积分-迭代法.将起伏面上的实测位场值,垂直投影至起伏面下部的一个水平面上,作为该水平面上的位场初始值.根据该水平面上的初始值,用积分方法计算起伏面上的位场值.用起伏面上的实测值与计算值的差值,对水平面上的位场值进行校正.如此反复迭代,直至起伏面上的实测值与计算值的差值小到可以忽略.有了水平面上的位场值后,就可以用积分的方法或其他方法计算水平面以上的任意曲面或水平面的位场值.该方法原理简单,不用解线性代数方程组,有较高的计算速度.它特别适用于位场向下延拓,有良好的延拓效果.本文还介绍了积分迭代法的应用实例.     

位场向下延拓的波数域广义逆算法

位场向下延拓是位场数据处理和反演中的重要运算,但是它的不稳定性影响了它在许多处理和反演方法技术中的应用.本文通过把位场向下延拓视为向上延拓的反问题,得到向下延拓的褶积型线性积分方程,再利用Fourier变换矩阵的正交对称特性,并结合矩阵的奇异值分解和广义逆原理,提出了一种稳定的不需要进行求逆运算的位场向下延拓广义逆方法——波数域广义逆算法,解决了位场大深度向下延拓的不稳定性问题.把这种方法用于三维理论模型数据和实际磁场数据的向下延拓获得了理想的结果.     

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

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

空间域位场延拓新方法研究

在空间域进行位场延拓,需要数值求解第一类Fredholm积分方程,由于所得方程组系数矩阵不是稀疏矩阵,求解该方程组需要的计算机内存大,计算量大,导致延拓算法在一般计算机上难以实现,阻碍了对空间域位场延拓方法的研究.在分析系数矩阵结构特征的基础上,本文证明了方程组系数矩阵是对称的分块Toeplitz型矩阵.利用系数矩阵的对称性和分块Toeplitz型矩阵与向量相乘的快速算法,解决了系数矩阵的存储和计算问题,使得空间域位场延拓成为可能,为研究新的位场延拓方法和分析延拓误差提供了一条新的途径.利用模型数据和实测资料,对空间域位场向上延拓、空间域积分迭代法向下延拓进行了检验,结果证实了空间域位场延拓的可行性和正确性.     

利用积分方程的加速迭代算法计算随钻电磁波电阻率测量仪器的响应

采用改进型逐次逼近解法(MSAM)和Aitken加速技术相结合的迭代算法计算二维积分方程.该算法将纵向成层原状地层作为背景地层,将计算区域限制在井眼和侵入带内,具有未知量数目少、收敛速度快、计算精度高的优点.利用该算法对随钻电磁波电阻率测量仪器在轴对称二维地层中的响应进行数值模拟.模拟结果显示,幅度衰减曲线和相位移曲线受井眼、侵入和围岩的影响程度不同,二者径向探测深度和垂向分辨率也有差异,利用补偿后的幅度衰减电阻率曲线和相位移电阻率曲线交叉点的坐标可精确确定地层层界面位置.     

场位解析延拓的稳定化算法

本文讨论场位解析延拓的稳定化算法。首先叙述问题的古典提法,并且为以后作准备,将它化成另一类不适定问题。接着建立问题解的连续依赖性估计。然后根据这个估计式,将问题化成求解一个条件变分问题。最后利用正则化方法解这个变分问题,得到问题解的一个稳定化的计算公式。     

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

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

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

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

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

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

逐步逼近曲化平方法研究

位场资料的常规处理解释方法是建立在平面数据理论之上的,然而实际的位场数据大多为曲面数据,若把曲面数据当成平面数据进行处理,必然导致很大的误差,因此有必要进行“曲化平”处理. 曲化平目前存在的主要问题是计算精度低、计算量大,因此,研究快速、精度高且适合大数据量处理的曲化平方法具有重要的价值. 本文在已有的泰勒级数曲化平方法基础上提出了逐步逼近技术和平均平面技术,使得曲面位场资料处理的精度得到了很大提高. 最后通过理论模型和实际资料的处理验证了该方法的应用效果.     

二维位场在有源空间的解析延拓

位场的向下延拓对重、磁等位场资料的处理与解释具有重要意义,但现有的延拓方法在位场向下接近场源时会产生剧烈震荡,计算过程很不稳定,而且理论上都不能越过场源. 本文在复数域里定义了一个具有位场和几何双重特性的新函数PFGR,提出了一套新的延拓方法,在理论上实现了位场在包含场源的下半空间里的解析延拓,经过模型试算验证,表明该方法可以延拓越过场源,准确确定场源奇点,计算过程稳定,对测量中的随机噪声也具有较好的抗干扰能力.     

位场波数域转换算法误差方程及其应用(英文)

偏移抽样理论是本文作者建立的更普遍的傅立叶变换数值计算的理论. 基于这一理论,作者在现文中导出了位场波数域转换算法误差方程,该方程不仅给出了更灵活的位场波数域转换算法,而且揭示了位场波数域转换中的误差规律.源于该方程的DFT0η η(0.5,0.5)化极技术可以大大提高低纬度(包括磁赤道)化极磁异常的分辨率和精度.该方程所揭示的波数域位场高通转换中边缘振荡的规律性(来源、形成机理和基本性质),从理论上指出了改善现有高通转换位场数据拓边技术效果的途径.     

位场向下延拓的最小曲率方法

针对位场向下延拓的不适定性,我们将位场向下延拓视为向上延拓的反问题,提出以位场最小曲率作为约束条件来求解稳定的下延位场.我们将剖面位场向上延拓表达式用傅里叶矩阵的形式表示,以矩阵乘法形式给出延拓的表达式,同时向待反演的下延位场引入最小曲率约束,得到向下延拓的最小曲率解,并利用正交变换给出了更为简洁的频率域解.随后,利用Kronecker积将上述全部结果拓展至三维位场,给出了三维位场向下延拓的最小曲率解.此外,我们将位场数据的填充、扩充问题与向下延拓问题统筹考虑,提出一种新的向下延拓迭代格式,该算法面向实际资料处理需求、无须预扩充或填补数据.下延迭代时,对原始数据直接向下延拓,而空白部分利用上一次下延位场估计的上延值替代其空白值并对其向下延拓,直至获得最小曲率约束下稳定的向下延拓结果.同时,我们也讨论了利用改进L曲线和广义交叉验证(GCV)计算正则参数最优估计的问题.对理论模型和实际航空重力资料进行了向下延拓检验,处理结果表明位场向下延拓的最小曲率方法解能满足实际位场资料对向下延拓处理的需求,具有较高的下延精度.