位场各阶垂向导数的ISVD算法及其应用

位场垂向导数广泛应用于位场数据处理与解释当中，目前求取位场各阶垂向导数方法大致可以分为两类：①在频率域进行；②在空间域进行。针对这两类方法存在着稳定性较差的问题，这里介绍了计算位场各阶垂向导数的ISVD算法，并通过对单一地质体模型和地质体组合模型产生的布格重力异常进行处理分析，证明了ISVD算法在计算位场各阶垂向导数时，特别是在计算位场高阶垂向导数时，比在频率域及空间域计算结果具有较好的稳定性，ISVD算法能够有效地识别出浅层小尺度地质体边界。为了验证该方法对实际资料的处理效果，对某区块实际布格重力异常数据应用ISVD算法试算各阶垂向导数，结合地质背景等资料划分出断裂体系。结果证明，应用ISVD算法计算的位场各阶垂向导数，对于确定浅层小尺度地质体边界有较好的效果。

Isvd algorithm and its application for calculating each order vertical derivative of potential field

The vertical derivative of potential field is widely used in processing and interpreting of potential field data. The methods of calculating each order vertical derivative of potential field can mainly be divided into two kinds： one is calcu- lated in the frequency domain, the other is in the space domain. Because of the problem of poor stability for these methods, this paper describes the ISVD algorithm for calculating each order vertical derivative of potential field. The result of process- ing and analyzing bouguer gravity anomalies of single geological model and the vertically integrated geological models causing by ISVD algorithm prove that each order vertical derivative of potential field by ISVD algorithm calculating, especially the high order vertical derivative, has better stability than it calculating in the frequency domain or space domain, and it can be used to effectively identify shallow small-scale geological boundaries. In order to verify the effect of ISVD algorithm applying to the actual data, each order vertical derivative of bouguer gravity anomaly data of a block by ISVD algorithm calculating i- dentify shallow small-scale geological boundaries and determine the fault system combining with the geological background information. The results prove that effect of determining shallow small-scale geological boundaries basing on each vertical derivative of potential field calculated by ISVD algorithm is better.