利用相关系数法确定三种位场垂向导数换算的滤波参数

垂向导数被广泛应用于位场数据的处理和解释中。在波数域中常规垂向导数算子受噪声的影响较大。通常解决方法是在垂向导数算子中引入低通滤波器,来压制高频噪声的干扰。Butterworth低通滤波、Hanning窗滤波以及向上延拓滤波是三种效果较好的低通滤波,但存在滤波参数较难选择的问题。针对这个问题,本文利用相关系数法确定各方法的滤波参数,避免了人为因素对求导结果的影响。模型试验和松辽盆地镇赉地区实际数据处理结果表明,使用相关系数法确定滤波参数后能快速获得较好的垂向导数结果,其中Butterworth低通滤波的导数结果相比其他两种方法的效果要好。     

位场垂向高阶导数的Tikhonov正则化迭代法

在位场数据处理中,垂向导数具有重要的物理意义。其在一定程度上可以划分不同深度和大小异常源产生的叠加异常,且导数的阶次越高,这种分辨能力就越强,但通常认为高阶导数的换算是不稳定的。本文在Tikhonov正则化求位场垂向高阶导数的基础上,结合迭代法进行逐次逼近,提出了位场高阶导数的Tikhonov正则化迭代法,并且得到Tikhonov正则化迭代法的递推公式。通过对该方法的滤波特性分析可以看出,该方法计算的位场垂向高阶导数具有一定的稳定性及保幅性。模型试验和实际数据的处理表明,该方法计算结果较常规FFT求导法有更高的稳定性和实用价值。     

Chebyshev低通滤波器在位场数据波数域导数换算中的应用

导数换算作为常用的位场数据处理手段具有重要的物理意义。由于实际测量误差和噪声的影响,波数域导数换算存在计算过程不稳定、求导结果精度差的问题。为了降低噪声干扰,提高导数计算精度,提出了在波数域常规导数算子基础上附加Chebyshev低通滤波器的任意阶导数换算方法。通过分析该滤波器的滤波特性,并结合位场异常径向平均功率谱曲线特征,确定了Chebyshev低通滤波法的滤波参数,减少了人为因素对求导结果的干扰。二度体和三度体模型试验表明,与其他3种方法(常规FFT(fast Fourier transformation)法、向上延拓法和ISVD(integrated second vertical derivative)法)相比,该方法计算的导数与理论值的均方根误差最小,求导结果精度较高。在虎林盆地布格重力异常数据处理中,利用该方法计算的垂向二阶导数受噪声干扰小,结果可靠性较高,依此划分出的8条大断裂和11条小断裂在以往研究中均得到证实。     

基于导数迭代法的重力归一化总梯度

针对传统重力归一化总梯度方法对油气藏分辨率不够和深度计算精度较低的问题,本文利用导数迭代法的向下延拓结果,计算重力归一化总梯度。将传统方法和本文方法同时用于模型试验,分析谐波数N、埋深和充满系数对G~H场特征的影响。结果表明,相比于传统方法,本文方法受谐波数N的影响较弱,G~H场对油气藏的分辨率更高,反演的油气藏深度更加准确。通过实际应用,证实了本文方法的有效性。