基于FFT的重力梯度相关成像研究

《基于FFT的重力梯度相关成像研究》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、2012年第3期王利宏等：基于FFtr的重力梯度相关成像研究15，JB：、y．z(6)关系数图。重力和重力梯度曲线能准确地和理论值1．2基于FF．r滤波的相关系数成像相对应，重力异常相关系数图中能辨别出浅层的上述重力异常和重力梯度异常的相关系数，其A1、A2和深层的C1，对其他异常均无反应。在重力成像结果横向上能分辨出异常体的横向位置，但纵梯度异常相关图中只能看出A1、A2引起的异常，主向上会忽略一些异常体。为了提高成像的纵向分辨要是因为浅层的异常引起的梯度异常比较尖锐，掩率，提出用FFTr对所测得重力异常和重力梯度异常盖了中深层引起的梯度异常，导致其相关系数成像进行空间滤波分析。效果不好。

2、必须通过滤波对资料进行进一步处理。在模型中，浅层小的突变体引起的异常小而尖，从重力异常Ag频谱、重力梯度分量Vzx频谱在频率域反映为高频。深层较大的异常体引起的异和重力梯度分量Vzz频谱中可以看出数据的主频集常大而平缓，在频率域反映为低频。通过低通、带通、中在0～20Hz和2O～60Hz，认为0～2OHz是深层高通滤波器分别将不同深度层异常数据提取出来，大范围异常引起的低频，20～60Hz是中层引起的利用逆变换得到过滤后的某一深度层的异常值。利中频，6OHz以上是浅层小范围引起的高频。用这一深度层密度异常体引起的异常数据进行归一滤波器分别设置低通、高通、带通滤波器。低通化相关计算，再将不同深

3、度层密度异常体引起异常滤波器的高频截止频率选择为20Hz，带通滤波器数据的相关系数进行叠加。最后用相关系数成像，能选择梯形框滤波器，阻带截止频率为20Hz和60提高成像的纵向分辨率。Hz，通带截止频率为25Hz和55Hz，高通滤波器的截止频率为60Hz，滤波器的频率响应函数如图22合成模型参数所示。合成模型由处于3个不同深度层的大小和剩余密度各异的7个直立长方体组成，这些长方体的几何参数和剩余密度如表1所列。表1模型几何参数和剩余密度序号沿长度，m沿z长度／m上顶深度／m)频率／HzA11oo1oo1005(a)低正滤波器A2l0olo010o—5B15o0300700—2．5B250030

4、07002．5B35003007oo2．5C1l50010002Ooo1．5C21O0010oO20oo-1图1显示了直立长方体在模型中的分布情况，其中Cl、C2分布在最深层，B1、B2、B3处于中间层，(b)带通滤波器形体最小的长方体A1、A2处于最浅层。频率／Hzc)高通滤波器图2滤波器频率响应利用低通滤波器得到的深层密度体的重力异常和重力梯度异常计算相关系数C1，利用带通滤波器图1模型示意得到中层密度体的重力异常和重力梯度异常计算相关系数c2，利用高通滤波器得到浅层密度体的重力3模拟结果分析异常和重力梯度异常计算相关系数C3，将C1、C2对模型重力异常、重力梯度异常进行计算得到和C3在

5、空间上进行叠加，得到最终的空间相关系其曲线，直接计算3个参数的相关系数得到地下相数，如图3所示。