New digital linear filters for Hankel J0 and J1 transforms

The numerical evaluation of certain integral transforms is required for the interpretation of some geophysical exploration data. Digital linear filter operators are widely used for carrying out such numerical integration. It is known that the method of Wiener–Hopf minimization of the error can be used to design very efficient, short digital linear filter operators for this purpose. We have found that, with appropriate modifications, this method can also be used to design longer filters. Two filters for the Hankel J0 transform (61-point and 120-point operators), and two for the Hankel J1 transform (47-point and 140-point operators) have been designed. For these transforms, the new filters give much lower errors compared to all other known filters of comparable, or somewhat longer, size. The new filter operators and some results of comparative performance tests with known integral transforms are presented. These filters would find widespread application in many numerical evaluation problems in geophysics.

[1]  D. Guptasarma.  OPTIMIZATION OF SHORT DIGITAL LINEAR FILTERS FOR INCREASED ACCURACY , 1982 .

[2]  O. Koefoed.  ERROR PROPAGATION AND UNCERTAINTY IN THE INTERPRETATION OF RESISTIVITY SOUNDING DATA , 1976 .

[3]  H. Johansen,et al.  FAST HANKEL TRANSFORMS , 1979 .

[4]  D. T. Biewinga,et al.  TRANSFORMATION OF DIPOLE RESISTIVITY SOUNDING MEASUREMENTS OVER LAYERED EARTH BY LINEAR DIGITAL FILTERING , 1974 .

[5]  O. Koefoed.  A note on the linear filter method of interpreting resistivity sounding data , 1972 .

[6]  E. Hashish,et al.  The fast Hankel transform1 , 1994 .

[7]  W. L. Anderson.  COMMENT ON “OPTIMIZED FAST HANKEL TRANSFORM FILTERS” BY NIELS BØIE CHRISTENSEN1 , 1991 .

[8]  Walter L. Anderson,et al.  Algorithm 588: Fast Hankel Transforms Using Related and Lagged Convolutions , 1982, TOMS.

[9]  U. Das,et al.  TRANSFORMATION OF DIPOLE TO SCHLUMBERGER SOUNDING CURVES BY MEANS OF DIGITAL LINEAR FILTERS , 1977 .

[10]  Walter L. Anderson,et al.  A hybrid fast Hankel transform algorithm for electromagnetic modeling , 1989 .

[11]  Luiz Rijo,et al.  Comment on ‘The fast Hankel transform’ by A.A. Mohsen and E.A. Hashish1 , 1996 .

[12]  K. Sørensen,et al.  The fields from a finite electrical dipole—A new computational approach , 1994 .

[13]  O. Koefoed,et al.  Geosounding Principles: Resistivity Sounding Measurements , 1980 .

[14]  D. Ghosh,et al.  COMPUTATION OF TYPE CURVES FOR ELECTROMAGNETIC DEPTH SOUNDING WITH A HORIZONTAL TRANSMITTING COIL BY MEANS OF A DIGITAL LINEAR FILTER , 1972 .

[15]  Walter L. Anderson,et al.  Improved digital filters for evaluating Fourier and Hankel transform integrals , 1975 .

[16]  O. Koefoed,et al.  Determination of resistivity sounding filters by the Wiener-Hopf least-squares method , 1979 .

[17]  D. Ghosh.  THE APPLICATION OF LINEAR FILTER THEORY TO THE ' DIRECT INTERPRETATION OF GEOELECTRICAL RESISTIVITY SOUNDING MEASUREMENTS * , 1971 .

[18]  J. Nissen,et al.  AN OPTIMIZED DIGITAL FILTER FOR THE FOURIER TRANSFORM , 1986 .

[19]  Niels B. Christensen,et al.  Optimized fast Hankel transform filters , 1990 .