低纬度磁异常化极方法应用效果对比

低纬度磁异常化极的方法很多,包括在空间域化极和在频率域化极,比较起来,频率域方法计算简单快速,目前化极主要在频率域进行.然而大部分低纬度化极方法只压制了噪声并没有针对化极因子做改动,也就没有真正达到压制南北向条带状拉长的目的.本文针对其中4种在频率域对化极因子做处理的方法(双曲正弦法,压制因子法,直接阻尼法,伪倾角法),设计模型进行对比,并用南海某地区的实测数据进行4种方法的化极效果对比,选出一种最符合南海地区实际情况的化极方法.并对原有的直接阻尼法和伪倾角法做了适当的改进.     

改进的阻尼法低纬度磁异常化极方法

中低纬度地区的磁异常资料进行化极处理能减少或消除斜磁化影响,对提高磁测资料的利用程度有重要意义.本文在总结分析原有的低纬度磁异常化极过程的基础上,改进了原有化极方法的滤波过程.改进后的阻尼化极方法,采用将振幅谱和相位谱分开进行滤波的过程,过程中只需要考虑对化极因子的振幅谱进行压制即可.通过公式分析表明该方法更加合理,减少了人为因素的影响;经单个模型和组合模型试验表明,此化极方法准确、稳定、可靠,可以有效地突出低纬度地区的异常体信息.     

Hartley变换化极

本文采用一种新颖的积分变换Hartley变换,系统推导了其频率域化极方法,提出在实数空间利用Hartley变换进行磁异常化极.结合现有的低纬度化极方法,讨论了低纬度Hartley变换化极特性,实现低纬度化极.理论模型计算表明建立的Hartley变换化极方法准确、可靠,实际资料处理表明该化极方法具有实用性.     

低纬度磁异常化极的伪倾角方法改进

基于改造化极因子的低磁纬度频率域化极方法具有计算速度快、控制参数少、操作简单、化极稳定等优点.本文分析压制因子法和伪倾角法的化极因子特征及其控制参数的影响,在此基础上改进伪倾角法的化极因子,即在磁偏角垂直方向及附近采用伪倾角法化极因子,而在其他方向采用常规频率域化极因子.改进后的伪倾角法既能有效压制磁偏角垂直方向及附近化极因子的放大作用,使得化极稳定,又能减少其他方向有效信号的化极特征的损失,提高化极精度.理论模型数据试验表明本文改进方法有效.利用本文改进方法对南海海域磁总场异常数据进行了变磁倾角化极,得到南海海域化极磁异常,这为研究南海大地构造特征和岩浆活动提供重要的参考资料.     

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

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

低纬度磁异常化极方法——压制因子法

磁异常化极转换是磁异常解释的重要基础，针对频率域磁异常化极现状，本文对低纬度化极的困难原因进行了深入分析，并提出了针对性的化极措施——压制因子法. 与已有研究比较，该方法具有简单直观、控制参数少因而便于实际操作的特点，模型计算表明了方法的有效性. 在我国南海低磁纬度地区实测磁异常的计算解释中，应用该化极方法，取得了良好的效果.     

一种改进的低纬度磁场化极方法——变频双向阻尼因子法

本文在考察分析原有的低纬度化极阻尼因子的基础上,提出了对原有阻尼因子的改造方法,即提出变频双向阻尼因子,改造后的阻尼因子不仅考虑了对振幅的压制,并兼顾了对相位谱的改造,同时还增加了频率控制因子.通过模型数据和实际资料检验,表明在低纬度地区,特别是在极低纬度地区,变频双向阻尼方法的化极结果能更有效地突出地质信息.     

磁赤道处化极方法

化向地磁极(化极)是最基本的磁测资料处理方法之一,化极能消除或减少斜磁化影响,提高对磁测资料的认识程度和解释水平,对研究地壳产生的磁异常具有重要意义.但低纬度地区特别是磁赤道处,化极处理很不稳定甚至奇异,一直是位场研究的难点.针对地磁纬度较低特别是磁赤道地区磁异常化极的困难,利用从磁北极处垂直磁化向低纬度地区水平磁化方向转换稳定的特点,提出"狭义化赤"概念,并将其与低纬度磁异常"倒相"解释方法结合,提出专门用于磁赤道处化极的方法.该方法扩展了现有的化极理论,实现了磁赤道处的稳定化极.区别于目前任何方法,专门用于(近)水平磁化条件下的化极计算,具有原理简单,实现方便,收敛速度快等特点.对理论模型和实际资料计算表明这种针对磁赤道地区磁异常的化极处理方法是稳定、可靠的.     

低纬度磁异常的初始模型约束全卷积神经网络化极方法

磁异常化极是磁测数据处理的重要基础工作,低纬度磁异常化极运算不稳定.深度学习是一种有效解决不稳定问题的数据驱动方法.然而,现有基于深度学习的位场数据处理仅以位场响应作为输入,缺乏场景外样本数据的通用性,往往导致深度学习网络预测结果的泛化性能不足.因此,为了克服这一问题以及进一步提高低纬度磁异常化极的精度.受“知识&数据”联合驱动深度学习网络结构的启发,本文提出一种基于初始模型约束的低纬度磁异常化极全卷积神经网络结构.该网络结构的输入包含两部分,一部分为低纬度磁异常,另外一部分为初始模型,网络输出为垂直磁化方向磁异常.训练与测试的磁异常数据采用基于网格点格架函数的空间域快速正演生成,以及对低纬度磁异常采用稳定的频率转换来获取初始模型,另外本文对数据集中10%的随机样本加入了5%的高斯噪声,以此增强全卷积神经网络结构的鲁棒性.理论模型试验验证了本文所提方法的有效性、精确性以及鲁棒性.最后,将本文方法应用于实际磁测数据,取得了良好的效果.     

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

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

添加正则化项的黏声逆时偏移成像方法研究

在实际地球介质中传播的地震波会产生衰减和频散现象,因此其更接近于黏弹性介质,在地震处理中补偿这些黏性影响是十分必要的。基于波动方程的叠前深度偏移中进行吸收衰减补偿更准确,也更有物理意义,直接求解双程波动方程的逆时偏移(RTM)能够成像大倾角复杂构造,具有诸多优势。然而当考虑吸收衰减补偿时通常会产生不稳定现象,大部分研究都是在逆时偏移的波场延拓中进行波数域的低通滤波来解决这个问题。本文采用广义标准线性固体的黏声波动方程进行吸收衰减补偿的Q--RTM方法,通过添加正则化项的方式来稳定延拓过程。添加正则化项本质上是低通滤波,滤波窗口是指数形式,在时空域有明确的形式,可以阻止发生高频不稳定。与直接滤波相比,正则化参数可以是空变的,因此比较适合剧烈变化的区域,我们还发现震源归一化互相关成像条件更适合Q--RTM方法。      

基于反演的稳定高效衰减补偿方法

反Q滤波方法是提高地震数据分辨率的一种有效途径,可以用来补偿振幅和校正相位.常规的反Q滤波方法一般基于波场延拓理论,具有不稳定性或振幅补偿不足的缺点.本文基于波场延拓的正Q滤波方程,借鉴反演的思想以及正则化策略提出了一种新的衰减补偿方法,该方法稳定、精确,利用该方法可最终得到高分辨的地震记录.该方法仅计算有效频带内的频率分量,提高了计算效率.模拟数据以及实际数据处理验证了本文方法的有效性.     

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

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

基于混合时频分析技术的地震数据噪声压制(英文)

针对复杂地质结构、陡倾角相干噪声、空间采样不均匀等情况下F-x域反褶积去噪技术的不足,提出首先应用具有时-频聚集性度量准则的广义S变换将时间-空间域的地震数据变换至时间-频率-空间域(t-f-x)的数据,在t-f-x域中对每一个频率切片应用经验模态分解(EMD),移除噪声占主导地位的本征模态函数以压制相干和随机噪声的滤波方法。模型分析表明第一本征模态函数表征的高频信息以噪声为主,移除第一本证模态函数可以达到压制噪声的目的。经广义S变换后形成t-f-x域中EMD滤波方法等效于具有依赖于空间位置、频率、高波数截断特征的自适应f-k滤波。此滤波方法考虑了数据的局部时-频特征,且具有执行简单的特点。与AR预测滤波方法比较,此法滤除的成分包含较少的低波数的信息,滤除的成分非常的局部化,且获得结果没有表现出过度平滑的特征。实际资料的应用表明在经广义S变换后形成t-f-x域中运用EMD滤波方法能够有效地压制随机和陡倾角相干噪声。     

一种基于正则化策略的稳定衰减补偿逆时偏移

地下介质通常具有黏滞性,地震波在地下介质的传播过程中将不可避免地伴随着振幅衰减和速度频散,进而影响地震成像的准确性和分辨率.基于黏滞介质的衰减补偿逆时偏移能沿地震波的传播路径恢复其所经历的振幅衰减和相位畸变,有效提升成像效果.然而,由于地层对地震波的吸收衰减效应呈指数变化,衰减补偿过程中高低频分量的非同步增长易导致补偿算法数值不稳定.为此,本文提出了一种基于正则化策略的稳定衰减补偿逆时偏移方法.该方法基于解耦的常Q分数阶拉普拉斯算子黏滞声波方程描述地震波在地下介质中的传播效应,将振幅正则化因子引入该方程的时间-波数域的解析解中,以确保在补偿过程中地震波场能稳定延拓.二维、三维合成数据以及实际资料的偏移算例均证实了该方法的可行性和有效性,所提出的方法能有效地处理衰减补偿中的不稳定问题,明显提升地震资料的偏移成像质量.     

地震层析成像的实用正则化

正则化是求解地震层析成像反演问题必须的。一般情况下,地震层析成像方程组非常庞大,合理选择正则化参数非常困难。本文提出一种在线性层析成像反演中实际选择正则化参数的算法。该算法基于大多数层析成像常用的数据统计假设。首先利用Lanczos双对角变换将方程组转换到Krylov子空间。在转换的子空间中,方程组转换成标准阻尼最小二乘正规方程。标准方程组的解可以写成正则化参数的显函数,这使得可以进行计算机方便的选取正则化参数。正则化参数选取的第二个准则是数字化模拟的研究。如果转换空间的维数远小于原始模型空间的维数,此算法的计算效果非常好,特别适用于大型地震层析成像问题。     

直升机航空瞬变电磁自适应正则化一维反演方法研究

本文以中心回线式直升机航空瞬变电磁法为例对时间域航空电磁资料的一维反演算法进行理论研究.由于时间域航空电磁反演问题是典型的不适定问题,多解性严重且稳定性能较差.为此,论文使用Tikhonov正则化反演方法进行一维解释,可以得到用最光滑模型、最平坦模型和最小构造模型约束条件下的反演结果,而正则化因子是根据各次迭代的数据目...     

基于分形径向谱的位场向下延拓截止波数自动确定方法

向下延拓是重磁位场数据处理与解释的一项重要技术,并因其固有的不稳定性而成为研究的热点.为获得稳定向下延拓结果,波数域向下延拓一般通过附加低通滤波器或改造向下延拓因子来完成.由此,滤波器或改造因子的截止波数则是精确向下延拓的关键.本文基于分形修正径向平均功率谱的物理特性,提出一种位场向下延拓截止波数的自动确定方法.基于理论重力模型和航磁实测数据的向下延拓对比实验结果表明：(1)本文所提出的自动确定方法物理意义明确,能快速有效地确定截止波数并进而获得向下延拓正则参数;(2)基于本文方法的正则化向下延拓结果优于改进导数迭代法的向下延拓结果.     

大地电磁自适应正则化反演算法

针对大地电磁正则化反演中正则化因子的选取困难问题提出了自适应正则化反演算法（Adaptive Regularized Inversion Algorithm, ARIA）. 在该算法中， ①提出了一种新的数据方差处理方法：数据方差规范化，使得数据方差的大小只对数据的拟合发生影响，不对数据目标函数和模型约束目标函数的权重产生影响，从而减少了正则化因子取值的影响因素；②提出了粗糙度核矩阵的概念，并给出了由基本结构插值基函数计算粗糙度核矩阵的公式，使得模型目标函数的构建更为简便、直接；③根据数据目标函数、模型约束目标函数和正则化因子之间的关系，提出了两种正则化因子自适应调节方法. 本文详细阐述了最平缓模型约束下的大地电磁一维连续介质反演的ARIA实现，以几个算例的分析比较来说明ARIA的有效性.     

位场各阶垂向导数换算的新正则化方法

位场垂向导数大量应用于位场数据处理与解释中.当前广泛采用的位场各阶垂向导数换算方法为基于Laplace方程并结合波数域和空间域方法的具有递推特性的ISVD(integrated second vertical derivative)算法.本文在位场垂向导数换算的正则化方法和径向平均功率谱的基础上,提出一种位场各阶垂向导数换算的新正则化方法.新正则化方法仅需通过分析位场径向平均功率谱来确定一个截止波数,即可稳定换算位场各阶垂向导数.理论模型和实测数据实验结果表明:(1)新正则化方法物理意义明确、计算简单,且各阶垂向导数换算的稳定性和精度明显优于ISVD算法;(2)在用新正则化方法求得各阶垂向导数的基础上,利用泰勒级数法可以获得大深度、高精度的位场向下延拓结果.