论文部分内容阅读
随着重磁勘探仪器探测精度与效率的提高,一些传统的位场数据处理方法因计算精度低、分析条件复杂等缺陷而难以满足当今位场勘探精细解释的需求。鉴于此,改进和开发高精度、高效率的位场数据处理方法具有重要意义。本文在位场数据转换与分离以及边界检测的反演工作方面提出了一些改进措施,取得了较好的效果。位场异常的转换与分离是最常规的异常处理解析方法,异常转换的精度以及异常分离的是否彻底直接影响一些正反演方法(如欧拉反褶积、界面异常反演等)的准确性和分辨率,本文以迭代补偿原理为基础,提出了用于导数换算和向下延拓的近似算子迭代法以及位场分离的波数域迭代法。近似算子迭代法是针对位场导数换算和向下延拓过程中存在的不稳定性问题提出的。该方法利用位场的泰勒级数、Laplace方程和不同阶次导数的转换关系,推导了导数换算和向下延拓的近似算子,将上述近似算子置入波数域迭代中,实现了位场导数换算和向下延拓的稳定计算。在分析近似算子迭代法滤波特性和收敛性的过程中,笔者发现可以利用邻次迭代结果的相关度确定最佳的迭代次数区间。含噪声模型试验表明,近似算子迭代法相比于其他方法具有较强的计算稳定性和较高的计算精度;相关系数与迭代次数的关系曲线可以快速提供最佳的迭代次数,保证了近似算子迭代法较高的运算效率和精度。同样以迭代补偿思想为前提,本文设计了一种幅值较低的区域场分离算子,给出了位场分离的波数域迭代法。该法能够将局部异常频谱中属于区域场的成分剥离出来,返还给区域异常,进而弥补区域异常的信号损失,最终达到改善异常分离效果的目的。模型试验表明,相比于常规向上延拓的异常分离方法,本文方法提取的目标场与真实值更加接近。利用位场资料反演地质体的边界位置可以为构造解释提供重、磁场划分依据。现有的边界检测方法大都与位场异常的梯度变化相关,因此这些方法常常受高频干扰影响,致使识别出的边界信息存在着模糊、缺失或虚假等缺点。针对于此,本文提出了基于差分求导的伪总梯度余弦值法、加权小子域滤波和结合增强型均衡滤波器的欧拉反褶积法三种边界检测技术。这三种方法有各自的特点,它们的处理结果可以相互弥补、验证,从而提高边界识别的精度和分辨率。基于差分求导的伪总梯度余弦值法是考虑了边界点异常的方向信息后,采用子域均值和差分方式计算梯度,以垂向梯度与伪总梯度(水平梯度的幅值被压缩)矢量夹角的余弦值作为滤波输出,结果中的极大值对应边界位置。该方法在计算过程中对边界位置的异常特征进行了两次紧缩,突出了边界与非边界异常的差异性。模型试验结果表明该法具有较好的计算稳定性和边界位置识别效果。加权小子域滤波技术是以直立岩脉重力异常、数值差分及小子域滤波为基础提出的,滤波结果以平面等值线中被紧缩的梯级带作为地质体边界的识别标志,以正、负圈闭异常反映局部地质体剩余密度的分析模式,在空间域内实现了对三方向小子域滤波的加权整合。模型试验表明,加权小子域滤波不但对地质体边界有较高的检测精度,而且能够清晰地反映出区域与局部场的圈闭异常关系(地质体的剩余密度关系)。针对常规Euler方程分散解较多,解的精度受构造指数影响较大的问题,笔者以一种新的增强型均衡滤波器为基础,推导了新的Euler方程。新方程不需要预先给定构造指数和背景场就可以计算出边界的平面位置和深度。模型试验表明,相比于常规Euler方程,新方程对边界的平面位置和深度的反演精度更高。以上边界检测方法中,基于差分求导的伪总梯度余弦值法能够清晰的检测出梯级带特征的构造边界;加权小子域滤波不仅弥补了前者检测非梯级带特征边界能力的不足,更能够反映出边界两侧地质体的异常强弱关系;结合增强型均衡滤波器的欧拉法反褶积法可以很好的填补前两者无法反演边界深度信息的空白。这三种方法配套使用可以更好的反映构造边界的三维分布特征。将以上方法应用到松辽盆地镇赉地区重磁场数据处理中,取得了一些成果。依据位场转换与分离方法的处理结果,结合收集的地质、物性、钻井资料,分析了研究区重、磁场、剩余异常、垂向导数及中、浅部地层的视物性分布特征;根据位场边界检测方法的处理结果,在研究区内划分出11条主要断裂,并将研究区划分为3个二级构造单元和9个次级构造单元;在此基础上,以岩石物性、钻井资料和地震剖面为约束条件,反演了研究区中、新生界底界面深度,预测了4处凹陷为油气勘探的有利远景区。