1. College of Geo-Exploration Science and Technology, Jilin University, Changchun 130026, China;
2. Department of Earth Sciences, Nanjing University, Nanjing 210093, China
Having investigated the design of classical FIR digital filter, the authors put forward a method of FIR filter design based on Chebyshev best uniform approximation which can be used to separate potential anomalies. In an experiment on the theoretical model, the maximum error of regional anomaly separation by low pass filter of Hanning window is 6.266×10-6 m/s2 , the mean square error is 2.115×10-6 m/s2, and the maximum percentage error is 22.2%. The errors are always higher than 10.1% when the computing points are beyond ±9 km. On the contrary, regional and residual anomalies separated by the best uniform approximation theory fit fairly well with the theoretical anomaly, with the curves in fairly good agreement. The maximum error of regional anomaly separation is 3.101×10-6 m/s2 , the mean square error is 0.989×10-6 m/s2 , and the maximum percentage error is 7.76% for several data on the boundary and is lower than 4.1% for the other points. The method was used to separate Bouguer gravity anomalies along the Sunwu-Jiayin profile with satisfactory results obtained. It is proved that the anomaly separation is quite complete.
[1] Pawlowski R S. Preferential continuation for potential-field anomaly enhancement[J].Geophysics,1995,60(2): 390.
[2] Pedersen L B. Relations between potential fields and some equivalent sources[J].Geophysics,1991,56: 961.
[3] Syberg F J R. A Fourier method for the regional residual problem of potential fields[J].Geophysical Prospecting, 1972,20:47.
[4] 张凤旭,张凤琴,刘财,等. 基于余弦变换的匹配滤波方法分离重磁异常[J].石油地球物理勘探,2006, 41(2), 216.
[5] 穆石铭,申宁华,孙运生. 区域地球物理数据处理方法及其应用[M].长春:吉林科学技术出版社,1990.
[6] 胡广书. 数字信号处理:理论、算法与实现[M]. 北京:清华大学出版社,1997.
[7] McClellan J H, Parks T W. A unified approach to the design of optimum FIR linear phase digital f-ilters[J]. IEEE Trans. Circuit Theory, 1973, 20:697.
[8] McClellan J H, Parks T W, Rabiner L R. A computer program for designing optimum FIR linear phase digital filters[J].IEEE Trans. Audio Electroacoustics, 1973, 21:506.
[9] Rabiner L R, McClellan J H, Parks T W. FIR digital filter design techniques using weighted Chebyshev approximation[J].Proc IEEE, 1975, 63:595.
[10] 大庆油田石油地质志编写组. 中国石油地质志(卷2):上册[M].北京:石油工业出版社,1993.