滑动窗口混合加权法磁异常保真去噪研究

实测磁异常通常含有明显的噪声干扰,常用的滑动窗口平均滤波法虽然有明显的去噪效果,但对有意义的异常信号也会造成明显的幅值和宽度上的失真.S-G平滑滤波法具有异常形态保真的优点,但在平稳区的去噪效果欠佳.本文提出异常评价识别具体方法,并据此将平均法与S-G法动态加权融合,从而实现了既能对异常区进行信号的保护,又能在全区达到有效去噪的目的,为磁异常的有效去噪形成简单而有效的方法.     

基于非均匀快速傅里叶变换的最小二乘反演地震数据重建

不规则采样地震数据的重建是地震数据分析处理的重要问题.本文给出了一种基于非均匀快速傅里叶变换的最小二乘反演地震数据重建的方法,在最小二乘反演插值方程中,引入正则化功率谱约束项,通过非均匀快速傅里叶变换和修改周期图的方式,自适应迭代修改约束项,使待插值数据的频谱越来越接近真实的频谱,采用预条件共轭梯度法迭代求解,保证了解的稳定性和收敛速度.理论模型和实际地震数据插值试验证明了本文方法能够去除空间假频,速度快、插值效果好,具有实用价值.     

南海重磁异常特征及火成岩分布(英文)

南海火成岩油气藏具有广阔的勘探前景,综合利用地球物理方法圈划与识别火成岩体、研究火成岩分布是火成岩油气藏研究的基础。针对南海重磁场特征,采用低纬度、变倾角化极技术进行了磁异常化极处理,利用优选延拓方法实现重磁异常分离并提取南海海域浅部火成岩重磁异常信息,利用磁异常三维相关成像给出南海火成岩的三维空间等效分布,在重磁梯度突出局部异常边界信息的基础上,通过梯度加权的重磁相关分析勾画不同类型火成岩的平面展布,火成岩的分布特征显示出受地壳深部结构及断裂构造的控制与影响。     

克里金法在GPS数据内插中的应用

介绍了克里金法在GPS数据内插中的应用、变差函数的自动拟合,以及在进行克里金内插编程时一些值得注意的问题等,通过实验结果看出内插过程取得了很好的效果。     

基于位场分离方法的极地冰盖重力异常剥离：以模型试验为例

为了提高极地重力数据处理的质量,降低巨厚冰层对重力数据解释的影响,将极地大陆冰盖引起的重力异常从叠加异常中剥离出来成为极地重力数据处理的重要步骤之一。针对该问题,本文建立两个极地大陆冰层构造模型,对不同模型进行正演获得对应布格重力异常,再尝试利用三种位场分离方法,包括趋势分析法、向上延拓法和优化滤波法,对冰层引起的重力异常进行剥离试验,基于模型试验结果对比分析不同位场分离方法对极地冰层重力效应剥离的应用效果。模型试验结果显示,趋势分析法在冰层界面起伏较平缓情况下效果较好,但当冰层界面起伏复杂时难以用数学表达式精确拟合区域场,因此易在结果中引入虚假异常;向上延拓法在局部场源较浅的异常分离中效果较好,但由于极地模型冰盖深度较深及延拓高度选取不合理,在试验中出现了分离不足的情况;优化滤波法在相对复杂模型中分离效果相对好,但由于频段选取的主观性,使得该方法也易出现分离不足或过度情况。综合上述模型试验表明,在剥离极地冰盖引起的重力异常时,由于极地冰盖厚度起伏等地质情况复杂多变,应结合实际地质情况及其他可获得的地球物理资料进行综合分析,选择适用于具体情况的分离方法。     

重力和重力梯度数据三维相关成像

本文提出了重力异常三维相关成像方法和重力梯度数据三维相关成像方法,并提出了基于异常分离的三维相关成像方法来提高成像分辨率.通过合成Y型岩脉模型和合成多个直立长方体组合模型的重力异常和重力梯度数据试验分析,验证了本文三维相关成像方法可显示出异常地质体的空间赋存状态和等效剩余质量分布,具有良好的纵向和横向分辨率.     

近地表密度估计的重力贝叶斯分析方法及在云南地区的应用

基于布格重力异常相对于地形起伏光滑分布的约束条件,从一维自由空气重力异常数据出发,采用贝叶斯方法估算近地表岩石密度,同时采用三次B样条函数拟合布格重力异常,获取光滑分布的布格重力异常.数据拟合和光滑约束之间的权重采用Akaike贝叶斯准则（ABIC准则）自动确定.均匀剖分模型和不均匀剖分模型数据试验都验证了该方法的有效性.相关参数评价表明,足够多的样条系数可以提高估计结果的准确性,样条系数的个数接近测点数时可获得较稳定的估计结果.增大异常的噪声水平时,ABIC准则可有效地自动增大先验光滑约束的权重.云南地区两条重力剖面应用结果表明,剖面沿线的近地表密度值起伏变化明显（达2.45~2.8 g·cm^-3）,前寒武纪和古生代地层密度相对较高（主要为2.53~2.75 g·cm^-3）,而中生代密度较低（2.45~2.73 g·cm^-3）;本文估计的近地表密度结果与区域物性资料及地表地质特征较吻合;估计的剖面布格重力异常具有光滑性;红河断裂两侧近地表密度差异较大,可达0.4 g·cm^-3.本文获得的两条剖面近地表密度结构和布格重力异常为该区深部结构与构造研究提供更可靠的重力基础数据.

     

优化滤波方法及其在中国大陆布格重力异常数据处理中的应用

在优选延拓法的理论基础上,研究提出基于格林等效层概念和维纳滤波器的优化滤波法,用于对重力异常数据进行去噪和分离.与传统向上延拓法和优选延拓法相比,优化滤波法分离异常与延拓高度无关,不需要已知延拓高度,具有一定的优势.理论重力模型数据的去噪和异常分离试验表明优化滤波法有效,异常分离效果优于传统向上延拓法和带通滤波法.利用优化滤波法对中国大陆重力异常数据去噪和异常分离,得到有效的布格重力异常和区域重力异常.以中国大陆深地震探测推断的莫霍面深度信息为约束,对区域重力异常数据进行密度界面约束反演,得到中国大陆莫霍面深度分布.本文方法为中国大陆深部探测和区域构造研究提供一定的技术支撑.     

基于GPU并行的重力、重力梯度三维正演快速计算及反演策略

利用NVIDIA CUDA编程平台,实现了基于GPU并行的重力、重力梯度三维快速正演计算方法.采用当前在重力数据约束反演或联合反演中流行的物性模型(密度大小不同、规则排列的长方体单元)作为地下剖分单元,对任意三维复杂模型体均可用很多物性模型进行组合近似,利用解析方法计算出所有物性模型在计算点的异常值并累加求和,得到整个模型体在某一计算点引起的重力(或重力梯度)值.针对精细的复杂模型体产生的问题,采用GPU并行计算技术,主要包括线程有效索引与优化的并行归约技术进行高效计算.在显卡型号为NVIDIA Quadro 2000相对于单线程CPU程序,重力和重力梯度Uxx、Uxy正演计算可以分别达到60与50倍的加速.本文还讨论了GPU并行计算在两种反演方法中的策略,为快速三维反演技术提供了借鉴.     

位场数据网格化的反插值法

位场不规则分布数据的网格化是位场数据分析处理的首要问题.本文借鉴反插值法的原理,提出利用基于预条件共轭梯度的反插值法实现位场数据的网格化.其中,插值算子采用高斯权系数,滤波算子采用Laplacian算子,预条件算子采用滤波算子的逆.通过理论模型和实际航磁数据的网格化试验分析,验证了本文的反插值法适合地球物理位场特征,网格化速度快,精度高,效果好.