IntroductionA new method for rapid computation of electromagnetic (EM) fields for high-frequency sounding (HFS) over a layered earth is presented in this report. The essence of the new method uses a Q-factor correction for extending a closed-form half-space analytic solution to a layered earth model. Use of the Q-factor in this context was first studied by Wait (1962; p. 53-61). Kraichman (1976) also discusses when the Q-factor method can be used to provide a good approximation to an exact layered earth solution.According to Wait (1962) and Kraichman (1976), the Q-factor approach may be used if the dipole fields can be approximated by a single normally incident plane wave at the receiver in the top layer. Fortunately, for the high-frequency range of interest with the HFS system, and with moderately conductive and dielectric earth models, the Q-factor method (also called the HFS Q-method) is usually a very good approximation to the exact solution for a layered earth. Examples of computing layered earth responses by the HFS Q-method and exact numerical integration method will be given for comparison of several typical models encountered in shallow HFS studies.Because the Q-factor method does not require any numerical integrations, the layered earth HFS Q-method algorithm is very fast (even on PCs). In general, the computations only require elementary complex functions, and in one case, modified Bessel functions of a complex argument. (See equations (l)-(3) below.)
Mathematical Theory for the HFS SystemThe mathematical theory and system development for the HFS system can be found in Stewart's (1993) PhD thesis, and the references contained therein. A case history paper on use of the HFS in field applications is also given by Stewart et al. (1994). In the latter paper, mention of a quick reconnaissance method of HFS data inversion is given by Anderson (1991), where complex image theory and the Q-factor method was first introduced to obtain fast approximate image solutions for layered models.To obtain improved forward and/or inverse solutions for HFS over a layered earth, I begin with an exact half-space formulation (Wait, 1954) for the magnetic H fields and electric E field components in cylindrical coordinates for a vertical magnetic dipole (VMD) loop source on the earth's surface, with exp(+/otf) time dependence, (5), r is the VMD transmitter and receiver separation, where r > 0 m, and m is the magnetic dipole moment (ampere-m2). The terms 7n(z) and Kn(z) in equation (2) are modified Bessel functions of order n and complex argument z defined in equation (5).In general, the propagation terms (70 in air, yl in earth) in equations (4) can have the electrical conductivity (a), magnetic permeability (/x), and dielectric permittivity (e) vary as a function of angular frequency (co). In addition, each of these parameters can also be complex. However, in this report, I will only consider the special case where each of the material parameters (a, p, e) are real and do not vary with frequency in the earth, or in each layer...