地震剖面线条化的一种简单算法

地震剖面线条化是应用模式识别方法自动追踪剖面中的同相轴, 并自动绘制线条图.首先探测地震道中的周波, 并描述周波属性, 然后根据相邻三道周波相关的目标函数最小原则追踪同相轴, 描述同相轴属性, 并记录同相轴属性表.根据同相轴属性表, 可对地震剖面进行基于同相轴属性的统计分析和滤波以及定量解释.本文对现有的地震剖面线条化方法进行了改进简化, 实现地震剖面同相轴快速自动追踪和自动绘制线条图.合成数据和实际数据试验表明该方法计算简单、快速、有效.      

基于剩余异常相关成像的重磁物性反演方法

将场源区剖分成长方体单元,通过采集的重磁数据反演出这些单元的密度或者磁化率变化,勾画出场源的分布图像,这种方式是重磁三维反演的重要方向.重磁相关成像通过计算测量的重磁异常与地下各点在测区上的重磁异常的归一化相关,显示出异常地质体的空间赋存状态和等效剩余重磁物性.该方法计算速度快,方法简单、稳定,但是反演的结果只是在-1到+1之间的等效物性,不能够直接反演剩余密度或者磁化率,并且无法引入已知的地质约束.本文通过对物性模型的正演和实测结果的残差进行相关成像,迭代更新物性模型实现对物性参数的反演过程.模型实验证明该方法相对相关成像不仅能提高分辨率,还能够得到真正的物性参数.     

危机矿山重磁资料精细处理与解释:以湖北省大冶铁矿为例

针对湖北大冶矿区复杂的地质问题, 基于高精度重磁实测数据, 采用三维重磁数据概率分析和物性反演技术, 实现了大冶矿区重磁资料的精细化处理与解释, 为大冶危机矿山深部及外围找矿提供了重要的参考.在无约束条件下对高磁异常源进行了三维概率成像反演, 并结合地质先验信息完成了矿区航磁异常的带地形人机交互三维反演, 预测了矿区沉积岩类与闪长岩体接触带的空间延展形态, 定量反演了矿区的三维物性特征, 反演结果显示在矿区深部(1000m以下) 及接触带弯折部位为有利的找矿区域, 目前通过实钻已得到证实.      

剩磁影响下磁性体磁化方向估计的一种简单算法

磁测资料处理解释通常需要已知场源的磁化方向信息.在剩磁影响下,总磁化方向与感应磁化方向并不一致,这使得处理解释变得复杂.本文研究提出从近等轴状磁总场异常估计磁化方向的偶极源互相关法,即用磁偶极源的理论磁异常与实测磁异常作互相关,用一系列递变的偶极源位置和磁化方向做试错,取互相关最大化对应的磁化方向作为估计结果.理论数据试验和中国西北某金属矿勘探区航磁资料试验表明本文方法能从近等轴状磁异常有效估计出总磁化方向.     

基于优选延拓的重力异常分离方法及其应用(英文)

本文在Pawlowski优选延拓思想的基础上,提出了在深源长波长信息和浅源短波长信息互不相关、延拓高度足够大的前提下,基于优选向上延拓算子的重力异常分离方法与步骤,通过分析理论模型重力异常分离算子的频率响应随不同延拓高度的变化规律,给出了估计最佳延拓高度的办法.利用本文方法对某铁矿区实际重力异常进行处理,有效地分离出了区域异常和局部异常,异常分离效果好.     

重力异常曲化平的改进插值—迭代法

重力异常曲化平即将起伏观测面上的重力异常延拓到平面上,从而为频率域处理和反演提供平坦观测面的重力异常数据.本文在常规插值—迭代法基础上,给出重力异常曲化平的改进插值—迭代法,即在异常迭代修正过程中引入起伏观测面修正因子,加快曲化平迭代收敛,促进曲化平效果提升.理论模型试验表明本文方法适用于观测面起伏较大、延拓跨度大的复...     

磁异常和梯度的频率域三维成像方法

三维反演在磁数据定量解释中具有重要作用.常用的空间域三维反演方法通常需要大量的正演和反演计算,因此对大规模数据的反演效率较低.三维成像是另一种定性和定量解释磁数据的重要方法.文中给出了一种磁异常与梯度三维成像的频率域迭代方法,该方法可以提高成像效率,适用于大规模数据的快速成像.笔者推导了磁总场异常和异常梯度频率域正演公式和成像公式,并将一种深度尺度因子引入成像公式中,提高了深度精度;笔者采用了迭代优化算法,减小了拟合误差,进一步提高了成像精度.通过理论模型数据试验和中国新疆某金属矿床实测数据,验证了本文方法的有效性、可行性.     

估计磁场源位置的倾斜角梯度算法改进

倾斜角梯度法在未知场源物性先验知识和构造指数情况下,可快速估计出场源的边界位置和深度分布,对于大面积实测磁异常的快速解释具有应用前景。笔者在Salem的倾斜角梯度法基础上,提出了一些算法改进。在导数换算上采用空间域有限差分法换算磁异常的水平导数,采用ISVD算法换算磁异常的高阶垂直导数。在剔除坏解上,采用磁异常广义倾斜角梯度的总水平梯度滤波法和场源深度有效性筛选法相结合。理论模型数据试验表明本方法简单快速,效果明显。     

反插值法实现地球物理数据快速网格化

本文借鉴反演理论，采用反插值法实现地球物理数据的快速网格化.首先由已知点与未知网格点的反距离拓扑关系建立反演方程，用已知点值直接计算出其所在网格的未知网格点值，再利用Laplacian算子滤波，使得模型光滑，且能量最小化.利用预条件共轭梯度法求解网格化方程组，并结合螺旋坐标系思想和Wilson\|Burg谱分解法，将二维数据的滤波处理转换到一维空间进行处理，实现二维数据快速滤波.另外，引入了阻尼约束，保证求解稳定、迭代收敛.最后，应用该方法对合成数据和实际资料进行了试验.     

青藏高原东北缘重力异常多尺度横向构造分析

本文研究了青藏高原东北缘地区布格重力异常特征,采用优化滤波法和归一化总水平梯度垂直导数法对研究区重力异常进行多尺度分离和横向构造分析.分离出的多尺度重力异常特征表明:1) 青藏高原东北缘地区大致以东经106°线为界,有一条醒目的重力异常梯级带,即贺兰山—六盘山—川滇南北构造带的北段,其东西两侧布格重力异常特征在形态和走向上截然不同,意味着两侧密度结构和构造特征存在明显差异. 2) 鄂尔多斯地块内部定边以北,重力异常高带走向由北东向转为近南北向,推测定边附近存在一个密度或构造界面,其两侧物质组成和构造特征具有差异,对比大尺度重力异常和中尺度重力异常,表明异常特征的这种差异主要是由上地幔深部结构引起的. 3) 青藏高原东北部各块体深部边界位置与地表构造分布不同,反映出该区构造复杂,深浅构造差异大. 4) 由于印度—欧亚板块碰撞及随后印度板块持续向北的挤压作用,造成青藏高原东北缘中、下地壳物质在巨大的北东向推挤力和鄂尔多斯刚性块体阻挡的共同作用下,沿着相对软弱的秦岭造山带方向蠕动.依据多尺度重力异常及其横向构造特征,综合推断出研究区内五条断裂带,即秦岭地轴北缘断裂带、海原—六盘山断裂带、香山—天景山断裂带、烟筒山断裂带和青铜峡—固原断裂带,并分析了它们在地壳深部的可能展布特征.