# Patent

## Title
Detecting unauthorized consumption of electrical energy

## Abstract

Aspects of detecting the unauthorized consumption of electrical energy are described. In some embodiments, a system includes a guided surface waveguide probe that launches a guided surface wave along a surface of a terrestrial medium. The system further includes metering systems that are distributed within a geographical region associated with the guided surface waveguide probe. The system also includes at least one computing device and memory storing computer instructions that cause the at least one computing device to generate an energy flow map using data obtained from the metering systems.

## Background
CROSS REFERENCE TO RELATED APPLICATIONS
This application is related to co-pending U.S. Non-provisional Patent Application entitled “Excitation and Use of Guided Surface Wave Modes on Lossy Media,” which was filed on Mar. 7, 2013 and assigned application Ser. No. 13/789,538, and was published on Sep. 11, 2014 as Publication Number US2014/0252886 A1, and which is incorporated herein by reference in its entirety. This application is also related to co-pending U.S. Non-provisional Patent Application entitled “Excitation and Use of Guided Surface Wave Modes on Lossy Media,” which was filed on Mar. 7, 2013 and assigned application Ser. No. 13/789,525, and was published on Sep. 11, 2014 as Publication Number US2014/0252865 A1, and which is incorporated herein by reference in its entirety. This application is further related to co-pending U.S. Non-provisional Patent Application entitled “Excitation and Use of Guided Surface Wave Modes on Lossy Media,” which was filed on Sep. 10, 2014 and assigned application Ser. No. 14/483,089, and which is incorporated herein by reference in its entirety. This application is further related to co-pending U.S. Non-provisional Patent Application entitled “Excitation and Use of Guided Surface Waves,” which was filed on Jun. 2, 2015 and assigned application Ser. No. 14/728,507, and which is incorporated herein by reference in its entirety. This application is further related to co-pending U.S. Non-provisional Patent Application entitled “Excitation and Use of Guided Surface Waves,” which was filed on Jun. 2, 2015 and assigned application Ser. No. 14/728,492, and which is incorporated herein by reference in its entirety.
BACKGROUND
For over a century, signals transmitted by radio waves involved radiation fields launched using conventional antenna structures. In contrast to radio science, electrical power distribution systems in the last century involved the transmission of energy guided along electrical conductors. This understanding of the distinction between radio frequency (RF) and power transmission has existed since the early 1900's.



BRIEF DESCRIPTION OF THE DRAWINGS
Many aspects of the present disclosure can be better understood with reference to the following drawings. The components in the drawings are not necessarily to scale, emphasis instead being placed upon clearly illustrating the principles of the disclosure. Moreover, in the drawings, like reference numerals designate corresponding parts throughout the several views.
FIG. 1 is a chart that depicts field strength as a function of distance for a guided electromagnetic field and a radiated electromagnetic field.
FIG. 2 is a drawing that illustrates a propagation interface with two regions employed for transmission of a guided surface wave according to various embodiments of the present disclosure.
FIG. 3 is a drawing that illustrates a guided surface waveguide probe disposed with respect to a propagation interface of FIG. 2 according to various embodiments of the present disclosure.
FIG. 4 is a plot of an example of the magnitudes of close-in and far-out asymptotes of first order Hankel functions according to various embodiments of the present disclosure.
FIGS. 5A and 5B are drawings that illustrate a complex angle of incidence of an electric field synthesized by a guided surface waveguide probe according to various embodiments of the present disclosure.
FIG. 6 is a graphical representation illustrating the effect of elevation of a charge terminal on the location where the electric field of FIG. 5A intersects with the lossy conducting medium at a Brewster angle according to various embodiments of the present disclosure.
FIG. 7 is a graphical representation of an example of a guided surface waveguide probe according to various embodiments of the present disclosure.
FIGS. 8A through 8C are graphical representations illustrating examples of equivalent image plane models of the guided surface waveguide probe of FIGS. 3 and 7 according to various embodiments of the present disclosure.
FIGS. 9A and 9B are graphical representations illustrating examples of single-wire transmission line and classic transmission line models of the equivalent image plane models of FIGS. 8B and 8C according to various embodiments of the present disclosure.
FIG. 10 is a flow chart illustrating an example of adjusting a guided surface waveguide probe of FIGS. 3 and 7 to launch a guided surface wave along the surface of a lossy conducting medium according to various embodiments of the present disclosure.
FIG. 11 is a plot illustrating an example of the relationship between a wave tilt angle and the phase delay of a guided surface waveguide probe of FIGS. 3 and 7 according to various embodiments of the present disclosure.
FIG. 12 is a drawing that illustrates an example of a guided surface waveguide probe according to various embodiments of the present disclosure.
FIG. 13 is a graphical representation illustrating the incidence of a synthesized electric field at a complex Brewster angle to match the guided surface waveguide mode at the Hankel crossover distance according to various embodiments of the present disclosure.
FIG. 14 is a graphical representation of an example of a guided surface waveguide probe of FIG. 12 according to various embodiments of the present disclosure.
FIG. 15A includes plots of an example of the imaginary and real parts of a phase delay (ΦU) of a charge terminal T1 of a guided surface waveguide probe according to various embodiments of the present disclosure.
FIG. 15B is a schematic diagram of the guided surface waveguide probe of FIG. 14 according to various embodiments of the present disclosure.
FIG. 16 is a drawing that illustrates an example of a guided surface waveguide probe according to various embodiments of the present disclosure.
FIG. 17 is a graphical representation of an example of a guided surface waveguide probe of FIG. 16 according to various embodiments of the present disclosure.
FIGS. 18A through 18C depict examples of receiving structures that can be employed to receive energy transmitted in the form of a guided surface wave launched by a guided surface waveguide probe according to the various embodiments of the present disclosure.
FIG. 18D is a flow chart illustrating an example of adjusting a receiving structure according to various embodiments of the present disclosure.
FIG. 19 depicts an example of an additional receiving structure that can be employed to receive energy transmitted in the form of a guided surface wave launched by a guided surface waveguide probe according to the various embodiments of the present disclosure.
FIG. 20A depicts a schematic symbol that represents various guided surface waveguide probes according to various embodiments of the present disclosure.
FIG. 20B depicts a schematic symbol that represents various guided surface wave receive structures according to various embodiments of the present disclosure.
FIG. 21 depicts a networked environment according to various embodiments of the present disclosure.
FIG. 22 depicts a scenario in which an unauthorized guided surface wave receive structure is operating in a geographical region according to various embodiments of the present disclosure.
FIG. 23A depicts an example of an expected energy usage map according to various embodiments of the present disclosure.
FIG. 23B depicts an example of an actual energy usage map according to various embodiments of the present disclosure.
FIGS. 24-26 are flow charts illustrating examples of an energy usage monitor generating expected energy usage maps according to various embodiments of the present disclosure.
FIGS. 27A-27B is a flow chart illustrating an example of an energy usage monitor determining whether an unauthorized guided surface wave receive structure is operating in a geographical region.
FIG. 28 is a schematic block diagram that provides one example illustration of a computing environment according to various embodiments of the present disclosure



DETAILED DESCRIPTION
To begin, some terminology shall be established to provide clarity in the discussion of concepts to follow. First, as contemplated herein, a formal distinction is drawn between radiated electromagnetic fields and guided electromagnetic fields.
As contemplated herein, a radiated electromagnetic field comprises electromagnetic energy that is emitted from a source structure in the form of waves that are not bound to a waveguide. For example, a radiated electromagnetic field is generally a field that leaves an electric structure such as an antenna and propagates through the atmosphere or other medium and is not bound to any waveguide structure. Once radiated electromagnetic waves leave an electric structure such as an antenna, they continue to propagate in the medium of propagation (such as air) independent of their source until they dissipate regardless of whether the source continues to operate. Once electromagnetic waves are radiated, they are not recoverable unless intercepted, and, if not intercepted, the energy inherent in the radiated electromagnetic waves is lost forever. Electrical structures such as antennas are designed to radiate electromagnetic fields by maximizing the ratio of the radiation resistance to the structure loss resistance. Radiated energy spreads out in space and is lost regardless of whether a receiver is present. The energy density of the radiated fields is a function of distance due to geometric spreading. Accordingly, the term “radiate” in all its forms as used herein refers to this form of electromagnetic propagation.
A guided electromagnetic field is a propagating electromagnetic wave whose energy is concentrated within or near boundaries between media having different electromagnetic properties. In this sense, a guided electromagnetic field is one that is bound to a waveguide and may be characterized as being conveyed by the current flowing in the waveguide. If there is no load to receive and/or dissipate the energy conveyed in a guided electromagnetic wave, then no energy is lost except for that dissipated in the conductivity of the guiding medium. Stated another way, if there is no load for a guided electromagnetic wave, then no energy is consumed. Thus, a generator or other source generating a guided electromagnetic field does not deliver real power unless a resistive load is present. To this end, such a generator or other source essentially runs idle until a load is presented. This is akin to running a generator to generate a 60 Hertz electromagnetic wave that is transmitted over power lines where there is no electrical load. It should be noted that a guided electromagnetic field or wave is the equivalent to what is termed a “transmission line mode.” This contrasts with radiated electromagnetic waves in which real power is supplied at all times in order to generate radiated waves. Unlike radiated electromagnetic waves, guided electromagnetic energy does not continue to propagate along a finite length waveguide after the energy source is turned off. Accordingly, the term “guide” in all its forms as used herein refers to this transmission mode of electromagnetic propagation.
Referring now to FIG. 1, shown is a graph 100 of field strength in decibels (dB) above an arbitrary reference in volts per meter as a function of distance in kilometers on a log-dB plot to further illustrate the distinction between radiated and guided electromagnetic fields. The graph 100 of FIG. 1 depicts a guided field strength curve 103 that shows the field strength of a guided electromagnetic field as a function of distance. This guided field strength curve 103 is essentially the same as a transmission line mode. Also, the graph 100 of FIG. 1 depicts a radiated field strength curve 106 that shows the field strength of a radiated electromagnetic field as a function of distance.
Of interest are the shapes of the curves 103 and 106 for guided wave and for radiation propagation, respectively. The radiated field strength curve 106 falls off geometrically (1/d, where d is distance), which is depicted as a straight line on the log-log scale. The guided field strength curve 103, on the other hand, has a characteristic exponential decay of e−ad/√{square root over (d)} and exhibits a distinctive knee 109 on the log-log scale. The guided field strength curve 103 and the radiated field strength curve 106 intersect at point 112, which occurs at a crossing distance. At distances less than the crossing distance at intersection point 112, the field strength of a guided electromagnetic field is significantly greater at most locations than the field strength of a radiated electromagnetic field. At distances greater than the crossing distance, the opposite is true. Thus, the guided and radiated field strength curves 103 and 106 further illustrate the fundamental propagation difference between guided and radiated electromagnetic fields. For an informal discussion of the difference between guided and radiated electromagnetic fields, reference is made to Milligan, T., Modern Antenna Design, McGraw-Hill, 1st Edition, 1985, pp. 8-9, which is incorporated herein by reference in its entirety.
The distinction between radiated and guided electromagnetic waves, made above, is readily expressed formally and placed on a rigorous basis. That two such diverse solutions could emerge from one and the same linear partial differential equation, the wave equation, analytically follows from the boundary conditions imposed on the problem. The Green function for the wave equation, itself, contains the distinction between the nature of radiation and guided waves.
In empty space, the wave equation is a differential operator whose eigenfunctions possess a continuous spectrum of eigenvalues on the complex wave-number plane. This transverse electro-magnetic (TEM) field is called the radiation field, and those propagating fields are called “Hertzian waves.” However, in the presence of a conducting boundary, the wave equation plus boundary conditions mathematically lead to a spectral representation of wave-numbers composed of a continuous spectrum plus a sum of discrete spectra. To this end, reference is made to Sommerfeld, A., “Uber die Ausbreitung der Wellen in der Drahtlosen Telegraphie,” Annalen der Physik, Vol. 28, 1909, pp. 665-736. Also see Sommerfeld, A., “Problems of Radio,” published as Chapter 6 in Partial Differential Equations in Physics—Lectures on Theoretical Physics: Volume VI, Academic Press, 1949, pp. 236-289, 295-296; Collin, R. E., “Hertzian Dipole Radiating Over a Lossy Earth or Sea: Some Early and Late 20th Century Controversies,” IEEE Antennas and Propagation Magazine, Vol. 46, No. 2, April 2004, pp. 64-79; and Reich, H. J., Ordnung, P. F, Krauss, H. L., and Skalnik, J. G., Microwave Theory and Techniques, Van Nostrand, 1953, pp. 291-293, each of these references being incorporated herein by reference in its entirety.
The terms “ground wave” and “surface wave” identify two distinctly different physical propagation phenomena. A surface wave arises analytically from a distinct pole yielding a discrete component in the plane wave spectrum. See, e.g., “The Excitation of Plane Surface Waves” by Cullen, A. L., (Proceedings of the IEE (British), Vol. 101, Part IV, August 1954, pp. 225-235). In this context, a surface wave is considered to be a guided surface wave. The surface wave (in the Zenneck-Sommerfeld guided wave sense) is, physically and mathematically, not the same as the ground wave (in the Weyl-Norton-FCC sense) that is now so familiar from radio broadcasting. These two propagation mechanisms arise from the excitation of different types of eigenvalue spectra (continuum or discrete) on the complex plane. The field strength of the guided surface wave decays exponentially with distance as illustrated by curve 103 of FIG. 1 (much like propagation in a lossy waveguide) and resembles propagation in a radial transmission line, as opposed to the classical Hertzian radiation of the ground wave, which propagates spherically, possesses a continuum of eigenvalues, falls off geometrically as illustrated by curve 106 of FIG. 1, and results from branch-cut integrals. As experimentally demonstrated by C. R. Burrows in “The Surface Wave in Radio Propagation over Plane Earth” (Proceedings of the IRE, Vol. 25, No. 2, February, 1937, pp. 219-229) and “The Surface Wave in Radio Transmission” (Bell Laboratories Record, Vol. 15, June 1937, pp. 321-324), vertical antennas radiate ground waves but do not launch guided surface waves.
To summarize the above, first, the continuous part of the wave-number eigenvalue spectrum, corresponding to branch-cut integrals, produces the radiation field, and second, the discrete spectra, and corresponding residue sum arising from the poles enclosed by the contour of integration, result in non-TEM traveling surface waves that are exponentially damped in the direction transverse to the propagation. Such surface waves are guided transmission line modes. For further explanation, reference is made to Friedman, B., Principles and Techniques of Applied Mathematics, Wiley, 1956, pp. pp. 214, 283-286, 290, 298-300.
In free space, antennas excite the continuum eigenvalues of the wave equation, which is a radiation field, where the outwardly propagating RF energy with Ez and Hφ in-phase is lost forever. On the other hand, waveguide probes excite discrete eigenvalues, which results in transmission line propagation. See Collin, R. E., Field Theory of Guided Waves, McGraw-Hill, 1960, pp. 453, 474-477. While such theoretical analyses have held out the hypothetical possibility of launching open surface guided waves over planar or spherical surfaces of lossy, homogeneous media, for more than a century no known structures in the engineering arts have existed for accomplishing this with any practical efficiency. Unfortunately, since it emerged in the early 1900's, the theoretical analysis set forth above has essentially remained a theory and there have been no known structures for practically accomplishing the launching of open surface guided waves over planar or spherical surfaces of lossy, homogeneous media.
According to the various embodiments of the present disclosure, various guided surface waveguide probes are described that are configured to excite electric fields that couple into a guided surface waveguide mode along the surface of a lossy conducting medium. Such guided electromagnetic fields are substantially mode-matched in magnitude and phase to a guided surface wave mode on the surface of the lossy conducting medium. Such a guided surface wave mode can also be termed a Zenneck waveguide mode. By virtue of the fact that the resultant fields excited by the guided surface waveguide probes described herein are substantially mode-matched to a guided surface waveguide mode on the surface of the lossy conducting medium, a guided electromagnetic field in the form of a guided surface wave is launched along the surface of the lossy conducting medium. According to one embodiment, the lossy conducting medium comprises a terrestrial medium such as the Earth.
Referring to FIG. 2, shown is a propagation interface that provides for an examination of the boundary value solutions to Maxwell's equations derived in 1907 by Jonathan Zenneck as set forth in his paper Zenneck, J., “On the Propagation of Plane Electromagnetic Waves Along a Flat Conducting Surface and their Relation to Wireless Telegraphy,” Annalen der Physik, Serial 4, Vol. 23, Sep. 20, 1907, pp. 846-866. FIG. 2 depicts cylindrical coordinates for radially propagating waves along the interface between a lossy conducting medium specified as Region 1 and an insulator specified as Region 2. Region 1 can comprise, for example, any lossy conducting medium. In one example, such a lossy conducting medium can comprise a terrestrial medium such as the Earth or other medium. Region 2 is a second medium that shares a boundary interface with Region 1 and has different constitutive parameters relative to Region 1. Region 2 can comprise, for example, any insulator such as the atmosphere or other medium. The reflection coefficient for such a boundary interface goes to zero only for incidence at a complex Brewster angle. See Stratton, J. A., Electromagnetic Theory, McGraw-Hill, 1941, p. 516.
According to various embodiments, the present disclosure sets forth various guided surface waveguide probes that generate electromagnetic fields that are substantially mode-matched to a guided surface waveguide mode on the surface of the lossy conducting medium comprising Region 1. According to various embodiments, such electromagnetic fields substantially synthesize a wave front incident at a complex Brewster angle of the lossy conducting medium that can result in zero reflection.
To explain further, in Region 2, where an ejωt field variation is assumed and where ρ≠0 and z≧0 (with z being the vertical coordinate normal to the surface of Region 1, and ρ being the radial dimension in cylindrical coordinates), Zenneck's closed-form exact solution of Maxwell's equations satisfying the boundary conditions along the interface are expressed by the following electric field and magnetic field components:




  
    
      
        
          
            H
            
              2
              ⁢
              φ
            
          
          =
          
            A
            ⁢
            
              
            
            ⁢
            
              e
              
                
                  -
                  
                    u
                    2
                  
                
                ⁢
                z
              
            
            ⁢
            
              
                H
                1
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γρ
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        
          
            E
            
              2
              ⁢
              ρ
            
          
          =
          
            
              A
              ⁡
              
                (
                
                  
                    u
                    2
                  
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    
                      ωɛ
                      o
                    
                  
                
                )
              
            
            ⁢
            
              e
              
                
                  -
                  
                    u
                    2
                  
                
                ⁢
                z
              
            
            ⁢
            
              
                H
                1
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γρ
                
                )
              
            
          
        
        ,
        and
      
    
    
      
        (
        2
        )
      
    
  
  
    
      
        
          E
          
            2
            ⁢
            z
          
        
        =
        
          
            A
            ⁡
            
              (
              
                
                  -
                  γ
                
                
                  ωɛ
                  o
                
              
              )
            
          
          ⁢
          
            e
            
              
                -
                
                  u
                  2
                
              
              ⁢
              z
            
          
          ⁢
          
            
              
                H
                0
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γρ
                
                )
              
            
            .
          
        
      
    
    
      
        (
        3
        )
      
    
  




In Region 1, where the ejωt field variation is assumed and where ρ≠0 and z≧0, Zenneck's closed-form exact solution of Maxwell's equations satisfying the boundary conditions along the interface is expressed by the following electric field and magnetic field components:




  
    
      
        
          
            H
            
              1
              ⁢
              φ
            
          
          =
          
            A
            ⁢
            
              
            
            ⁢
            
              e
              
                
                  u
                  1
                
                ⁢
                z
              
            
            ⁢
            
              
                H
                1
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γρ
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        4
        )
      
    
  
  
    
      
        
          
            E
            
              1
              ⁢
              ρ
            
          
          =
          
            
              A
              ⁡
              
                (
                
                  
                    -
                    
                      u
                      1
                    
                  
                  
                    
                      σ
                      1
                    
                    +
                    
                      j
                      ⁢
                      
                        
                      
                      ⁢
                      
                        ωɛ
                        1
                      
                    
                  
                
                )
              
            
            ⁢
            
              e
              
                
                  u
                  1
                
                ⁢
                z
              
            
            ⁢
            
              
                H
                1
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γρ
                
                )
              
            
          
        
        ,
        and
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          E
          
            1
            ⁢
            z
          
        
        =
        
          
            A
            ⁡
            
              (
              
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γ
                
                
                  
                    σ
                    1
                  
                  +
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    
                      ωɛ
                      1
                    
                  
                
              
              )
            
          
          ⁢
          
            e
            
              
                u
                1
              
              ⁢
              z
            
          
          ⁢
          
            
              
                H
                0
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γρ
                
                )
              
            
            .
          
        
      
    
    
      
        (
        6
        )
      
    
  




In these expressions, z is the vertical coordinate normal to the surface of Region 1 and p is the radial coordinate, Hn(2)(−jγφ is a complex argument Hankel function of the second kind and order n, u1 is the propagation constant in the positive vertical (z) direction in Region 1, u2 is the propagation constant in the vertical (z) direction in Region 2, σ1 is the conductivity of Region 1, ω is equal to 2πf, where f is a frequency of excitation, ∈o is the permittivity of free space, ∈1 is the permittivity of Region 1, A is a source constant imposed by the source, and γ is a surface wave radial propagation constant.
The propagation constants in the ±z directions are determined by separating the wave equation above and below the interface between Regions 1 and 2, and imposing the boundary conditions. This exercise gives, in Region 2,




  
    
      
        
          u
          2
        
        =
        
          
            
              -
              j
            
            ⁢
            
              
            
            ⁢
            
              k
              o
            
          
          
            
              1
              +
              
                (
                
                  
                    ɛ
                    r
                  
                  -
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    x
                  
                
                )
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




and gives, in Region 1,

u1=−u2(∈r−jx).  (8)

The radial propagation constant γ is given by





  
    
      
        
          γ
          =
          
            
              j
              ⁢
              
                
                  
                    k
                    o
                    2
                  
                  +
                  
                    u
                    2
                    2
                  
                
              
            
            =
            
              j
              ⁢
              
                
                  
                    k
                    o
                  
                  ⁢
                  n
                
                
                  
                    1
                    +
                    
                      n
                      2
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        9
        )
      
    
  




which is a complex expression where n is the complex index of refraction given by

n=√{square root over (∈r−jx)}.  (10)

In all of the above Equations,





  
    
      
        
          x
          =
          
            
              σ
              1
            
            
              ωɛ
              o
            
          
        
        ,
        and
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        
          
            k
            o
          
          =
          
            
              ω
              ⁢
              
                
                  
                    μ
                    o
                  
                  ⁢
                  
                    ɛ
                    o
                  
                
              
            
            =
            
              
                λ
                o
              
              
                2
                ⁢
                π
              
            
          
        
        ,
      
    
    
      
        (
        12
        )
      
    
  




where ∈r comprises the relative permittivity of Region 1, σ1 is the conductivity of Region 1, ∈o is the permittivity of free space, and μo comprises the permeability of free space. Thus, the generated surface wave propagates parallel to the interface and exponentially decays vertical to it. This is known as evanescence.

Thus, Equations (1)-(3) can be considered to be a cylindrically-symmetric, radially-propagating waveguide mode. See Barlow, H. M., and Brown, J., Radio Surface Waves, Oxford University Press, 1962, pp. 10-12, 29-33. The present disclosure details structures that excite this “open boundary” waveguide mode. Specifically, according to various embodiments, a guided surface waveguide probe is provided with a charge terminal of appropriate size that is fed with voltage and/or current and is positioned relative to the boundary interface between Region 2 and Region 1. This may be better understood with reference to FIG. 3, which shows an example of a guided surface waveguide probe 200a that includes a charge terminal T1 elevated above a lossy conducting medium 203 (e.g., the Earth) along a vertical axis z that is normal to a plane presented by the lossy conducting medium 203. The lossy conducting medium 203 makes up Region 1, and a second medium 206 makes up Region 2 and shares a boundary interface with the lossy conducting medium 203.
According to one embodiment, the lossy conducting medium 203 can comprise a terrestrial medium such as the planet Earth. To this end, such a terrestrial medium comprises all structures or formations included thereon whether natural or man-made. For example, such a terrestrial medium can comprise natural elements such as rock, soil, sand, fresh water, sea water, trees, vegetation, and all other natural elements that make up our planet. In addition, such a terrestrial medium can comprise man-made elements such as concrete, asphalt, building materials, and other man-made materials. In other embodiments, the lossy conducting medium 203 can comprise some medium other than the Earth, whether naturally occurring or man-made. In other embodiments, the lossy conducting medium 203 can comprise other media such as man-made surfaces and structures such as automobiles, aircraft, man-made materials (such as plywood, plastic sheeting, or other materials) or other media.
In the case where the lossy conducting medium 203 comprises a terrestrial medium or Earth, the second medium 206 can comprise the atmosphere above the ground. As such, the atmosphere can be termed an “atmospheric medium” that comprises air and other elements that make up the atmosphere of the Earth. In addition, it is possible that the second medium 206 can comprise other media relative to the lossy conducting medium 203.
The guided surface waveguide probe 200a includes a feed network 209 that couples an excitation source 212 to the charge terminal T1 via, e.g., a vertical feed line conductor. According to various embodiments, a charge Q1 is imposed on the charge terminal T1 to synthesize an electric field based upon the voltage applied to terminal T1 at any given instant. Depending on the angle of incidence (θi) of the electric field (E), it is possible to substantially mode-match the electric field to a guided surface waveguide mode on the surface of the lossy conducting medium 203 comprising Region 1.
By considering the Zenneck closed-form solutions of Equations (1)-(6), the Leontovich impedance boundary condition between Region 1 and Region 2 can be stated as

{circumflex over (z)}×{right arrow over (H)}2(ρ,φ,0)={right arrow over (J)}S,  (13)

where {circumflex over (z)} is a unit normal in the positive vertical (+z) direction and {right arrow over (H)}2 is the magnetic field strength in Region 2 expressed by Equation (1) above. Equation (13) implies that the electric and magnetic fields specified in Equations (1)-(3) may result in a radial surface current density along the boundary interface, where the radial surface current density can be specified by

Jρ=(ρ′)=−A H1(2)(−jγρ′)  (14)

where {circumflex over (z)} is a constant. Further, it should be noted that close-in to the guided surface waveguide probe 200 (for ρ<<λ), Equation (14) above has the behavior





  
    
      
        
          
            J
            close
          
          ⁡
          
            (
            
              ρ
              ′
            
            )
          
        
        =
        
          
            
              -
              
                A
                ⁡
                
                  (
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                  )
                
              
            
            
              π
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    γρ
                    ′
                  
                
                )
              
            
          
          =
          
            
              -
              
                H
                φ
              
            
            =
            
              -
              
                
                  
                    I
                    o
                  
                  
                    2
                    ⁢
                    
                      πρ
                      ′
                    
                  
                
                .
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




The negative sign means that when source current (I0) flows vertically upward as illustrated in FIG. 3, the “close-in” ground current flows radially inward. By field matching on Hφ “close-in,” it can be determined that





  
    
      
        A
        =
        
          
            -
            
              
                
                  I
                  o
                
                ⁢
                γ
              
              4
            
          
          =
          
            -
            
              
                ω
                ⁢
                
                  
                
                ⁢
                
                  q
                  1
                
                ⁢
                γ
              
              4
            
          
        
      
    
    
      
        (
        16
        )
      
    
  




where q1=C1V1, in Equations (1)-(6) and (14). Therefore, the radial surface current density of Equation (14) can be restated as





  
    
      
        
          
            J
            ρ
          
          ⁡
          
            (
            
              ρ
              ′
            
            )
          
        
        =
        
          
            
              
                I
                o
              
              ⁢
              γ
            
            4
          
          ⁢
          
            
              
                H
                1
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  -
                  
                    jγρ
                    ′
                  
                
                )
              
            
            .
          
        
      
    
    
      
        (
        17
        )
      
    
  




The fields expressed by Equations (1)-(6) and (17) have the nature of a transmission line mode bound to a lossy interface, not radiation fields that are associated with groundwave propagation. See Barlow, H. M. and Brown, J., Radio Surface Waves, Oxford University Press, 1962, pp. 1-5.

At this point, a review of the nature of the Hankel functions used in Equations (1)-(6) and (17) is provided for these solutions of the wave equation. One might observe that the Hankel functions of the first and second kind and order n are defined as complex combinations of the standard Bessel functions of the first and second kinds

Hn(1)(x)=Jn(x)+jNn(x), and  (18)

Hn(2)(x)=Jn(x)−jNn(x),  (19)

These functions represent cylindrical waves propagating radially inward (Hn(1)) and outward (Hn(2)), respectively. The definition is analogous to the relationship e+jx=cos x±j sin x. See, for example, Harrington, R. F., Time-Harmonic Fields, McGraw-Hill, 1961, pp. 460-463.

That Hn(2)(kρρ) is an outgoing wave can be recognized from its large argument asymptotic behavior that is obtained directly from the series definitions of Jn(x) and Nn(x). Far-out from the guided surface waveguide probe:




  
    
      
        
          
            
              
                H
                n
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                x
                )
              
            
            ⁢
            
              →
              
                x
                →
                ∞
              
            
            ⁢
            
              
                
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    j
                  
                  
                    π
                    ⁢
                    
                      
                    
                    ⁢
                    x
                  
                
              
              ⁢
              
                j
                n
              
              ⁢
              
                e
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  x
                
              
            
          
          =
          
            
              
                2
                
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  x
                
              
            
            ⁢
            
              j
              n
            
            ⁢
            
              e
              
                -
                
                  j
                  ⁡
                  
                    (
                    
                      x
                      -
                      
                        π
                        4
                      
                    
                    )
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        
          20
          ⁢
          a
        
        )
      
    
  




which, when multiplied by ejωt, is an outward propagating cylindrical wave of the form ej(ωt−kρ) with a 1/√{square root over (ρ)} spatial variation. The first order (n=1) solution can be determined from Equation (20a) to be





  
    
      
        
          
            
              H
              1
              
                (
                2
                )
              
            
            ⁡
            
              (
              x
              )
            
          
          ⁢
          
            →
            
              x
              →
              ∞
            
          
          ⁢
          
            
              
                
                  2
                  ⁢
                  
                    
                  
                  ⁢
                  j
                
                
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  x
                
              
            
            ⁢
            
              e
              
                
                  -
                  j
                
                ⁢
                
                  
                
                ⁢
                x
              
            
          
        
        =
        
          
            
              2
              
                π
                ⁢
                
                  
                
                ⁢
                x
              
            
          
          ⁢
          
            
              e
              
                -
                
                  j
                  ⁡
                  
                    (
                    
                      x
                      -
                      
                        π
                        2
                      
                      -
                      
                        π
                        4
                      
                    
                    )
                  
                
              
            
            .
          
        
      
    
    
      
        (
        
          20
          ⁢
          b
        
        )
      
    
  




Close-in to the guided surface waveguide probe (for ρ<<λ), the Hankel function of first order and the second kind behaves as





  
    
      
        
          
            H
            1
            
              (
              2
              )
            
          
          ⁡
          
            (
            x
            )
          
        
        ⁢
        
          →
          
            x
            →
            0
          
        
        ⁢
        
          
            
              2
              ⁢
              
                
              
              ⁢
              j
            
            
              π
              ⁢
              
                
              
              ⁢
              x
            
          
          .
        
      
    
    
      
        (
        21
        )
      
    
  




Note that these asymptotic expressions are complex quantities. When x is a real quantity, Equations (20b) and (21) differ in phase by √{square root over (j)}, which corresponds to an extra phase advance or “phase boost” of 45° or, equivalently, λ/8. The close-in and far-out asymptotes of the first order Hankel function of the second kind have a Hankel “crossover” or transition point where they are of equal magnitude at a distance of ρ=Rx.

Thus, beyond the Hankel crossover point the “far out” representation predominates over the “close-in” representation of the Hankel function. The distance to the Hankel crossover point (or Hankel crossover distance) can be found by equating Equations (20b) and (21) for −jγp, and solving for Rx. With x=σ/ω∈o, it can be seen that the far-out and close-in Hankel function asymptotes are frequency dependent, with the Hankel crossover point moving out as the frequency is lowered. It should also be noted that the Hankel function asymptotes may also vary as the conductivity (σ) of the lossy conducting medium changes. For example, the conductivity of the soil can vary with changes in weather conditions.
Referring to FIG. 4, shown is an example of a plot of the magnitudes of the first order Hankel functions of Equations (20b) and (21) for a Region 1 conductivity of σ=0.010 mhos/m and relative permittivity ∈r=15, at an operating frequency of 1850 kHz. Curve 115 is the magnitude of the far-out asymptote of Equation (20b) and curve 118 is the magnitude of the close-in asymptote of Equation (21), with the Hankel crossover point 121 occurring at a distance of Rx=54 feet. While the magnitudes are equal, a phase offset exists between the two asymptotes at the Hankel crossover point 121. It can also be seen that the Hankel crossover distance is much less than a wavelength of the operation frequency.
Considering the electric field components given by Equations (2) and (3) of the Zenneck closed-form solution in Region 2, it can be seen that the ratio of Ez and Eρ asymptotically passes to




  
    
      
        
          
            
              E
              z
            
            
              E
              ρ
            
          
          =
          
            
              
                
                  (
                  
                    
                      
                        -
                        j
                      
                      ⁢
                      
                        
                      
                      ⁢
                      γ
                    
                    
                      u
                      2
                    
                  
                  )
                
                ⁢
                
                  
                    
                      H
                      0
                      
                        (
                        2
                        )
                      
                    
                    ⁡
                    
                      (
                      
                        
                          -
                          j
                        
                        ⁢
                        
                          
                        
                        ⁢
                        γρ
                      
                      )
                    
                  
                  
                    
                      H
                      1
                      
                        (
                        2
                        )
                      
                    
                    ⁡
                    
                      (
                      
                        
                          -
                          j
                        
                        ⁢
                        
                          
                        
                        ⁢
                        γρ
                      
                      )
                    
                  
                
              
              ⁢
              
                →
                
                  ρ
                  →
                  ∞
                
              
              ⁢
              
                
                  
                    ɛ
                    r
                  
                  -
                  
                    j
                    ⁢
                    
                      σ
                      
                        ωɛ
                        0
                      
                    
                  
                
              
            
            =
            
              n
              =
              
                tan
                ⁢
                
                  
                
                ⁢
                
                  θ
                  i
                
              
            
          
        
        ,
      
    
    
      
        (
        22
        )
      
    
  




where n is the complex index of refraction of Equation (10) and θi is the angle of incidence of the electric field. In addition, the vertical component of the mode-matched electric field of Equation (3) asymptotically passes to





  
    
      
        
          
            E
            
              2
              ⁢
              
                
              
              ⁢
              z
            
          
          ⁢
          
            →
            
              ρ
              →
              ∞
            
          
          ⁢
          
            
              (
              
                
                  q
                  free
                
                
                  ɛ
                  0
                
              
              )
            
            ⁢
            
              
                
                  γ
                  3
                
                
                  8
                  ⁢
                  π
                
              
            
            ⁢
            
              e
              
                
                  -
                  
                    u
                    2
                  
                
                ⁢
                z
              
            
            ⁢
            
              
                e
                
                  -
                  
                    j
                    ⁡
                    
                      (
                      
                        γρ
                        -
                        
                          π
                          /
                          4
                        
                      
                      )
                    
                  
                
              
              
                ρ
              
            
          
        
        ,
      
    
    
      
        (
        23
        )
      
    
  




which is linearly proportional to free charge on the isolated component of the elevated charge terminal's capacitance at the terminal voltage, qfree=Cfree×VT.

For example, the height H1 of the elevated charge terminal T1 in FIG. 3 affects the amount of free charge on the charge terminal T1. When the charge terminal T1 is near the ground plane of Region 1, most of the charge Q1 on the terminal is “bound.” As the charge terminal T1 is elevated, the bound charge is lessened until the charge terminal T1 reaches a height at which substantially all of the isolated charge is free.
The advantage of an increased capacitive elevation for the charge terminal T1 is that the charge on the elevated charge terminal T1 is further removed from the ground plane, resulting in an increased amount of free charge qfree to couple energy into the guided surface waveguide mode. As the charge terminal T1 is moved away from the ground plane, the charge distribution becomes more uniformly distributed about the surface of the terminal. The amount of free charge is related to the self-capacitance of the charge terminal T1.
For example, the capacitance of a spherical terminal can be expressed as a function of physical height above the ground plane. The capacitance of a sphere at a physical height of h above a perfect ground is given by

Celevated sphere=4π∈oa(1+M+M2+M3+2M4+3M5+ . . . ),  (24)

where the diameter of the sphere is 2a, and where M=a/2h with h being the height of the spherical terminal. As can be seen, an increase in the terminal height h reduces the capacitance C of the charge terminal. It can be shown that for elevations of the charge terminal T1 that are at a height of about four times the diameter (4D=8a) or greater, the charge distribution is approximately uniform about the spherical terminal, which can improve the coupling into the guided surface waveguide mode.

In the case of a sufficiently isolated terminal, the self-capacitance of a conductive sphere can be approximated by C=4π∈oa, where a is the radius of the sphere in meters, and the self-capacitance of a disk can be approximated by C=8∈oa, where a is the radius of the disk in meters. The charge terminal T1 can include any shape such as a sphere, a disk, a cylinder, a cone, a torus, a hood, one or more rings, or any other randomized shape or combination of shapes. An equivalent spherical diameter can be determined and used for positioning of the charge terminal T1.
This may be further understood with reference to the example of FIG. 3, where the charge terminal T1 is elevated at a physical height of hp=H1 above the lossy conducting medium 203. To reduce the effects of the “bound” charge, the charge terminal T1 can be positioned at a physical height that is at least four times the spherical diameter (or equivalent spherical diameter) of the charge terminal T1 to reduce the bounded charge effects.
Referring next to FIG. 5A, shown is a ray optics interpretation of the electric field produced by the elevated charge Q1 on charge terminal T1 of FIG. 3. As in optics, minimizing the reflection of the incident electric field can improve and/or maximize the energy coupled into the guided surface waveguide mode of the lossy conducting medium 203. For an electric field (∈∥) that is polarized parallel to the plane of incidence (not the boundary interface), the amount of reflection of the incident electric field may be determined using the Fresnel reflection coefficient, which can be expressed as




  
    
      
        
          
            
              Γ
              ||
            
            ⁡
            
              (
              
                θ
                i
              
              )
            
          
          =
          
            
              
                E
                
                  ||
                  
                    ,
                    R
                  
                
              
              
                E
                
                  ||
                  
                    ,
                    i
                  
                
              
            
            =
            
              
                
                  
                    
                      (
                      
                        
                          ɛ
                          r
                        
                        -
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          x
                        
                      
                      )
                    
                    -
                    
                      
                        sin
                        2
                      
                      ⁢
                      
                        θ
                        i
                      
                    
                  
                
                -
                
                  
                    (
                    
                      
                        ɛ
                        r
                      
                      -
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        x
                      
                    
                    )
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    θ
                    i
                  
                
              
              
                
                  
                    
                      (
                      
                        
                          ɛ
                          r
                        
                        -
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          x
                        
                      
                      )
                    
                    -
                    
                      
                        sin
                        2
                      
                      ⁢
                      
                        θ
                        i
                      
                    
                  
                
                +
                
                  
                    (
                    
                      
                        ɛ
                        r
                      
                      -
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        x
                      
                    
                    )
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    θ
                    i
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        25
        )
      
    
  




where θi is the conventional angle of incidence measured with respect to the surface normal.

In the example of FIG. 5A, the ray optic interpretation shows the incident field polarized parallel to the plane of incidence having an angle of incidence of θi, which is measured with respect to the surface normal ({circumflex over (z)}). There will be no reflection of the incident electric field when Γ∥(θi)=0 and thus the incident electric field will be completely coupled into a guided surface waveguide mode along the surface of the lossy conducting medium 203. It can be seen that the numerator of Equation (25) goes to zero when the angle of incidence is

θi=arctan(√{square root over (∈r−jx)})=θi,B,  (26)

where x=σ/ω∈o. This complex angle of incidence (θi,B) is referred to as the Brewster angle. Referring back to Equation (22), it can be seen that the same complex Brewster angle (θi,B) relationship is present in both Equations (22) and (26).

As illustrated in FIG. 5A, the electric field vector E can be depicted as an incoming non-uniform plane wave, polarized parallel to the plane of incidence. The electric field vector E can be created from independent horizontal and vertical components as

{right arrow over (E)}(θi)=Eρ{circumflex over (ρ)}+Ez{circumflex over (z)}.  (27)

Geometrically, the illustration in FIG. 5A suggests that the electric field vector E can be given by





  
    
      
        
          
            
              E
              ρ
            
            ⁡
            
              (
              
                ρ
                ,
                z
              
              )
            
          
          =
          
            
              E
              ⁡
              
                (
                
                  ρ
                  ,
                  z
                
                )
              
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              θ
              i
            
          
        
        ,
        and
      
    
    
      
        (
        
          28
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            
              E
              z
            
            ⁡
            
              (
              
                ρ
                ,
                z
              
              )
            
          
          =
          
            
              
                E
                ⁡
                
                  (
                  
                    ρ
                    ,
                    z
                  
                  )
                
              
              ⁢
              
                cos
                ⁡
                
                  (
                  
                    
                      π
                      2
                    
                    -
                    
                      θ
                      i
                    
                  
                  )
                
              
            
            =
            
              
                E
                ⁡
                
                  (
                  
                    ρ
                    ,
                    z
                  
                  )
                
              
              ⁢
              sin
              ⁢
              
                
              
              ⁢
              
                θ
                i
              
            
          
        
        ,
      
    
    
      
        (
        
          28
          ⁢
          b
        
        )
      
    
  




which means that the field ratio is





  
    
      
        
          
            E
            ρ
          
          
            E
            z
          
        
        =
        
          
            1
            
              tan
              ⁢
              
                
              
              ⁢
              
                θ
                i
              
            
          
          =
          
            tan
            ⁢
            
              
            
            ⁢
            
              
                ψ
                i
              
              .
            
          
        
      
    
    
      
        (
        29
        )
      
    
  




A generalized parameter W, called “wave tilt,” is noted herein as the ratio of the horizontal electric field component to the vertical electric field component given by




  
    
      
        
          W
          =
          
            
              
                E
                ρ
              
              
                E
                z
              
            
            =
            
              
                
                W
                
              
              ⁢
              
                e
                
                  j
                  ⁢
                  
                    
                  
                  ⁢
                  Ψ
                
              
            
          
        
        ,
        or
      
    
    
      
        (
        
          30
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            1
            W
          
          =
          
            
              
                E
                ρ
              
              
                E
                z
              
            
            =
            
              
                tan
                ⁢
                
                  
                
                ⁢
                
                  θ
                  i
                
              
              =
              
                
                  1
                  
                    
                    W
                    
                  
                
                ⁢
                
                  e
                  
                    
                      -
                      j
                    
                    ⁢
                    
                      
                    
                    ⁢
                    Ψ
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        
          30
          ⁢
          b
        
        )
      
    
  




which is complex and has both magnitude and phase. For an electromagnetic wave in Region 2, the wave tilt angle (Ψ) is equal to the angle between the normal of the wave-front at the boundary interface with Region 1 and the tangent to the boundary interface. This may be easier to see in FIG. 5B, which illustrates equi-phase surfaces of an electromagnetic wave and their normals for a radial cylindrical guided surface wave. At the boundary interface (z=0) with a perfect conductor, the wave-front normal is parallel to the tangent of the boundary interface, resulting in W=0. However, in the case of a lossy dielectric, a wave tilt W exists because the wave-front normal is not parallel with the tangent of the boundary interface at z=0.

Applying Equation (30b) to a guided surface wave gives




  
    
      
        
          tan
          ⁢
          
            
          
          ⁢
          
            θ
            
              i
              ,
              B
            
          
        
        =
        
          
            
              E
              z
            
            
              E
              ρ
            
          
          =
          
            
              
                u
                2
              
              γ
            
            =
            
              
                
                  
                    ɛ
                    r
                  
                  -
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    x
                  
                
              
              =
              
                n
                =
                
                  
                    1
                    W
                  
                  =
                  
                    
                      1
                      
                        
                        W
                        
                      
                    
                    ⁢
                    
                      
                        e
                        
                          
                            -
                            j
                          
                          ⁢
                          
                            
                          
                          ⁢
                          Ψ
                        
                      
                      .
                    
                  
                
              
            
          
        
      
    
    
      
        (
        31
        )
      
    
  




With the angle of incidence equal to the complex Brewster angle (θi,B), the Fresnel reflection coefficient of Equation (25) vanishes, as shown by





  
    
      
        
          
            Γ
            ||
          
          ⁡
          
            (
            
              θ
              
                i
                ,
                B
              
            
            )
          
        
        =
        
          
            
              
                
                  
                    
                      (
                      
                        
                          ɛ
                          r
                        
                        -
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          x
                        
                      
                      )
                    
                    -
                    
                      
                        sin
                        2
                      
                      ⁢
                      
                        θ
                        i
                      
                    
                  
                
                -
                
                  
                    (
                    
                      
                        ɛ
                        r
                      
                      -
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        x
                      
                    
                    )
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    θ
                    i
                  
                
              
              
                
                  
                    
                      (
                      
                        
                          ɛ
                          r
                        
                        -
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          x
                        
                      
                      )
                    
                    -
                    
                      
                        sin
                        2
                      
                      ⁢
                      
                        θ
                        i
                      
                    
                  
                
                +
                
                  
                    (
                    
                      
                        ɛ
                        r
                      
                      -
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        x
                      
                    
                    )
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    θ
                    i
                  
                
              
            
            ⁢
            
              |
              
                
                  θ
                  i
                
                =
                
                  θ
                  
                    i
                    ,
                    B
                  
                
              
            
          
          =
          0.
        
      
    
    
      
        (
        32
        )
      
    
  




By adjusting the complex field ratio of Equation (22), an incident field can be synthesized to be incident at a complex angle at which the reflection is reduced or eliminated. Establishing this ratio as n=√{square root over (∈r−jx)} results in the synthesized electric field being incident at the complex Brewster angle, making the reflections vanish.

The concept of an electrical effective height can provide further insight into synthesizing an electric field with a complex angle of incidence with a guided surface waveguide probe 200. The electrical effective height (heff) has been defined as




  
    
      
        
          h
          eff
        
        =
        
          
            1
            
              I
              0
            
          
          ⁢
          
            
              ∫
              0
              
                h
                p
              
            
            ⁢
            
              
                I
                ⁡
                
                  (
                  z
                  )
                
              
              ⁢
              
                
              
              ⁢
              d
              ⁢
              
                
              
              ⁢
              z
            
          
        
      
    
    
      
        (
        33
        )
      
    
  




for a monopole with a physical height (or length) of hp. Since the expression depends upon the magnitude and phase of the source distribution along the structure, the effective height (or length) is complex in general. The integration of the distributed current I(z) of the structure is performed over the physical height of the structure (hp), and normalized to the ground current (I0) flowing upward through the base (or input) of the structure. The distributed current along the structure can be expressed by

I(z)=IC cos(β0z),  (34)

where β0 is the propagation factor for current propagating on the structure. In the example of FIG. 3, IC is the current that is distributed along the vertical structure of the guided surface waveguide probe 200a. 

For example, consider a feed network 209 that includes a low loss coil (e.g., a helical coil) at the bottom of the structure and a vertical feed line conductor connected between the coil and the charge terminal T1. The phase delay due to the coil (or helical delay line) is θc=βpIC, with a physical length of IC and a propagation factor of




  
    
      
        
          
            β
            p
          
          =
          
            
              
                2
                ⁢
                π
              
              
                λ
                p
              
            
            =
            
              
                2
                ⁢
                π
              
              
                
                  V
                  f
                
                ⁢
                
                  λ
                  0
                
              
            
          
        
        ,
      
    
    
      
        (
        35
        )
      
    
  




where Vf is the velocity factor on the structure, λ0 is the wavelength at the supplied frequency, and λp is the propagation wavelength resulting from the velocity factor Vf. The phase delay is measured relative to the ground (stake) current I0.

In addition, the spatial phase delay along the length lw of the vertical feed line conductor can be given by θy=βwlw where βw is the propagation phase constant for the vertical feed line conductor. In some implementations, the spatial phase delay may be approximated by θy=βwhp, since the difference between the physical height hp of the guided surface waveguide probe 200a and the vertical feed line conductor length lw is much less than a wavelength at the supplied frequency (λ0). As a result, the total phase delay through the coil and vertical feed line conductor is Φ=θc+θy, and the current fed to the top of the coil from the bottom of the physical structure is

IC(θc+θy)=I0ejΦ,  (36)

with the total phase delay Φ measured relative to the ground (stake) current I0. Consequently, the electrical effective height of a guided surface waveguide probe 200 can be approximated by

heff=1/I0∫0hpI0ejΦ cos(β0z)dz≅hpejΦ,  (37)

for the case where the physical height hp<<λ0. The complex effective height of a monopole, heff=hp at an angle (or phase shift) of Φ, may be adjusted to cause the source fields to match a guided surface waveguide mode and cause a guided surface wave to be launched on the lossy conducting medium 203.

In the example of FIG. 5A, ray optics are used to illustrate the complex angle trigonometry of the incident electric field (E) having a complex Brewster angle of incidence (θi,B) at the Hankel crossover distance (Rx) 121. Recall from Equation (26) that, for a lossy conducting medium, the Brewster angle is complex and specified by




  
    
      
        
          tan
          ⁢
          
            
          
          ⁢
          
            θ
            
              i
              ,
              B
            
          
        
        =
        
          
            
              
                ɛ
                r
              
              -
              
                j
                ⁢
                
                  σ
                  
                    ωɛ
                    o
                  
                
              
            
          
          =
          
            n
            .
          
        
      
    
    
      
        (
        38
        )
      
    
  




Electrically, the geometric parameters are related by the electrical effective height (heff) of the charge terminal T1 by

Rx tan ψi,B=Rx×W=heff=hpejΦ,  (39)

where ψi,B=(π/2)−θi,B is the Brewster angle measured from the surface of the lossy conducting medium. To couple into the guided surface waveguide mode, the wave tilt of the electric field at the Hankel crossover distance can be expressed as the ratio of the electrical effective height and the Hankel crossover distance





  
    
      
        
          
            h
            eff
          
          
            R
            x
          
        
        =
        
          
            tan
            ⁢
            
              
            
            ⁢
            
              ψ
              
                i
                ,
                B
              
            
          
          =
          
            
              W
              
                R
                ⁢
                
                  
                
                ⁢
                x
              
            
            .
          
        
      
    
    
      
        (
        40
        )
      
    
  




Since both the physical height (hp) and the Hankel crossover distance (Rx) are real quantities, the angle (Ψ) of the desired guided surface wave tilt at the Hankel crossover distance (Rx) is equal to the phase (Φ) of the complex effective height (heff). This implies that by varying the phase at the supply point of the coil, and thus the phase shift in Equation (37), the phase, Φ, of the complex effective height can be manipulated to match the angle of the wave tilt, ‘1’, of the guided surface waveguide mode at the Hankel crossover point 121: Φ=Ψ.

In FIG. 5A, a right triangle is depicted having an adjacent side of length R, along the lossy conducting medium surface and a complex Brewster angle ψi,B measured between a ray 124 extending between the Hankel crossover point 121 at Rx and the center of the charge terminal T1, and the lossy conducting medium surface 127 between the Hankel crossover point 121 and the charge terminal T1. With the charge terminal T1 positioned at physical height hp and excited with a charge having the appropriate phase delay Φ, the resulting electric field is incident with the lossy conducting medium boundary interface at the Hankel crossover distance Rx, and at the Brewster angle. Under these conditions, the guided surface waveguide mode can be excited without reflection or substantially negligible reflection.
If the physical height of the charge terminal T1 is decreased without changing the phase shift Φ of the effective height (heff), the resulting electric field intersects the lossy conducting medium 203 at the Brewster angle at a reduced distance from the guided surface waveguide probe 200. FIG. 6 graphically illustrates the effect of decreasing the physical height of the charge terminal T1 on the distance where the electric field is incident at the Brewster angle. As the height is decreased from h3 through h2 to h1, the point where the electric field intersects with the lossy conducting medium (e.g., the Earth) at the Brewster angle moves closer to the charge terminal position. However, as Equation (39) indicates, the height H1 (FIG. 3) of the charge terminal T1 should be at or higher than the physical height (hp) in order to excite the far-out component of the Hankel function. With the charge terminal T1 positioned at or above the effective height (heff), the lossy conducting medium 203 can be illuminated at the Brewster angle of incidence (ψi,B=(π/2)−θi,B) at or beyond the Hankel crossover distance (Rx) 121 as illustrated in FIG. 5A. To reduce or minimize the bound charge on the charge terminal T1, the height should be at least four times the spherical diameter (or equivalent spherical diameter) of the charge terminal T1 as mentioned above.
A guided surface waveguide probe 200 can be configured to establish an electric field having a wave tilt that corresponds to a wave illuminating the surface of the lossy conducting medium 203 at a complex Brewster angle, thereby exciting radial surface currents by substantially mode-matching to a guided surface wave mode at (or beyond) the Hankel crossover point 121 at Rx.
Referring to FIG. 7, shown is a graphical representation of an example of a guided surface waveguide probe 200b that includes a charge terminal T1. An AC source 212 acts as the excitation source for the charge terminal T1, which is coupled to the guided surface waveguide probe 200b through a feed network 209 (FIG. 3) comprising a coil 215 such as, e.g., a helical coil. In other implementations, the AC source 212 can be inductively coupled to the coil 215 through a primary coil. In some embodiments, an impedance matching network may be included to improve and/or maximize coupling of the AC source 212 to the coil 215.
As shown in FIG. 7, the guided surface waveguide probe 200b can include the upper charge terminal T1 (e.g., a sphere at height hp) that is positioned along a vertical axis z that is substantially normal to the plane presented by the lossy conducting medium 203. A second medium 206 is located above the lossy conducting medium 203. The charge terminal T1 has a self-capacitance CT. During operation, charge Q1 is imposed on the terminal T1 depending on the voltage applied to the terminal T1 at any given instant.
In the example of FIG. 7, the coil 215 is coupled to a ground stake 218 at a first end and to the charge terminal T1 via a vertical feed line conductor 221. In some implementations, the coil connection to the charge terminal T1 can be adjusted using a tap 224 of the coil 215 as shown in FIG. 7. The coil 215 can be energized at an operating frequency by the AC source 212 through a tap 227 at a lower portion of the coil 215. In other implementations, the AC source 212 can be inductively coupled to the coil 215 through a primary coil.
The construction and adjustment of the guided surface waveguide probe 200 is based upon various operating conditions, such as the transmission frequency, conditions of the lossy conducting medium (e.g., soil conductivity σ and relative permittivity ∈r), and size of the charge terminal T1. The index of refraction can be calculated from Equations (10) and (11) as

n=√{square root over (∈r−jx)},  (41)

where x=σ/ω∈o with ω=2πf. The conductivity σ and relative permittivity ∈r can be determined through test measurements of the lossy conducting medium 203. The complex Brewster angle (θi,B) measured from the surface normal can also be determined from Equation (26) as

θi,B=arctan(√{square root over (∈r−jx)}),  (42)

or measured from the surface as shown in FIG. 5A as





  
    
      
        
          ψ
          
            i
            ,
            B
          
        
        =
        
          
            π
            2
          
          -
          
            
              θ
              
                i
                ,
                B
              
            
            .
          
        
      
    
    
      
        (
        43
        )
      
    
  




The wave tilt at the Hankel crossover distance (WRx) can also be found using Equation (40).

The Hankel crossover distance can also be found by equating the magnitudes of Equations (20b) and (21) for −jγρ, and solving for Rx as illustrated by FIG. 4. The electrical effective height can then be determined from Equation (39) using the Hankel crossover distance and the complex Brewster angle as

heff=hpejΦ=Rx tan ψi,B.  (44)

As can be seen from Equation (44), the complex effective height (heff) includes a magnitude that is associated with the physical height (hp) of the charge terminal T1 and a phase delay (Φ) that is to be associated with the angle (Ψ) of the wave tilt at the Hankel crossover distance (Rx). With these variables and the selected charge terminal T1 configuration, it is possible to determine the configuration of a guided surface waveguide probe 200.

With the charge terminal T1 positioned at or above the physical height (hp), the feed network 209 (FIG. 3) and/or the vertical feed line connecting the feed network to the charge terminal T1 can be adjusted to match the phase (Φ) of the charge Q1 on the charge terminal T1 to the angle (Ψ) of the wave tilt (W). The size of the charge terminal T1 can be chosen to provide a sufficiently large surface for the charge Q1 imposed on the terminals. In general, it is desirable to make the charge terminal T1 as large as practical. The size of the charge terminal T1 should be large enough to avoid ionization of the surrounding air, which can result in electrical discharge or sparking around the charge terminal.
The phase delay θc of a helically-wound coil can be determined from Maxwell's equations as has been discussed by Corum, K. L. and J. F. Corum, “RF Coils, Helical Resonators and Voltage Magnification by Coherent Spatial Modes,” Microwave Review, Vol. 7, No. 2, September 2001, pp. 36-45, which is incorporated herein by reference in its entirety. For a helical coil with H/D>1, the ratio of the velocity of propagation (v) of a wave along the coil's longitudinal axis to the speed of light (c), or the “velocity factor,” is given by




  
    
      
        
          
            V
            f
          
          =
          
            
              v
              c
            
            =
            
              1
              
                
                  1
                  +
                  
                    20
                    ⁢
                    
                      
                        (
                        
                          D
                          s
                        
                        )
                      
                      2.5
                    
                    ⁢
                    
                      
                        (
                        
                          D
                          
                            λ
                            o
                          
                        
                        )
                      
                      0.5
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        45
        )
      
    
  




where H is the axial length of the solenoidal helix, D is the coil diameter, N is the number of turns of the coil, s=H/N is the turn-to-turn spacing (or helix pitch) of the coil, and λo is the free-space wavelength. Based upon this relationship, the electrical length, or phase delay, of the helical coil is given by





  
    
      
        
          θ
          c
        
        =
        
          
            
              β
              p
            
            ⁢
            H
          
          =
          
            
              
                
                  2
                  ⁢
                  π
                
                
                  λ
                  p
                
              
              ⁢
              H
            
            =
            
              
                
                  2
                  ⁢
                  π
                
                
                  
                    V
                    f
                  
                  ⁢
                  
                    λ
                    0
                  
                
              
              ⁢
              
                H
                .
              
            
          
        
      
    
    
      
        (
        46
        )
      
    
  




The principle is the same if the helix is wound spirally or is short and fat, but Vf and θc are easier to obtain by experimental measurement. The expression for the characteristic (wave) impedance of a helical transmission line has also been derived as





  
    
      
        
          Z
          c
        
        =
        
          
            
              60
              
                V
                f
              
            
            ⁡
            
              [
              
                
                  ln
                  ⁡
                  
                    (
                    
                      
                        
                          V
                          f
                        
                        ⁢
                        
                          λ
                          0
                        
                      
                      D
                    
                    )
                  
                
                -
                1.027
              
              ]
            
          
          .
        
      
    
    
      
        (
        47
        )
      
    
  




The spatial phase delay θy of the structure can be determined using the traveling wave phase delay of the vertical feed line conductor 221 (FIG. 7). The capacitance of a cylindrical vertical conductor above a prefect ground plane can be expressed as




  
    
      
        
          
            C
            A
          
          =
          
            
              
                2
                ⁢
                
                  πɛ
                  o
                
                ⁢
                
                  h
                  w
                
              
              
                
                  ln
                  ⁡
                  
                    (
                    
                      h
                      a
                    
                    )
                  
                
                -
                1
              
            
            ⁢
            Farads
          
        
        ,
      
    
    
      
        (
        48
        )
      
    
  




where hw is the vertical length (or height) of the conductor and a is the radius (in mks units). As with the helical coil, the traveling wave phase delay of the vertical feed line conductor can be given by





  
    
      
        
          
            θ
            y
          
          =
          
            
              
                β
                w
              
              ⁢
              
                h
                w
              
            
            =
            
              
                
                  
                    2
                    ⁢
                    π
                  
                  
                    λ
                    w
                  
                
                ⁢
                
                  h
                  w
                
              
              =
              
                
                  
                    2
                    ⁢
                    π
                  
                  
                    
                      V
                      w
                    
                    ⁢
                    
                      λ
                      0
                    
                  
                
                ⁢
                
                  h
                  w
                
              
            
          
        
        ,
      
    
    
      
        (
        49
        )
      
    
  




where βw is the propagation phase constant for the vertical feed line conductor, hw is the vertical length (or height) of the vertical feed line conductor, VW is the velocity factor on the wire, λ0 is the wavelength at the supplied frequency, and λw is the propagation wavelength resulting from the velocity factor Vw. For a uniform cylindrical conductor, the velocity factor is a constant with Vw≈0.94, or in a range from about 0.93 to about 0.98. If the mast is considered to be a uniform transmission line, its average characteristic impedance can be approximated by





  
    
      
        
          
            Z
            w
          
          =
          
            
              60
              
                V
                
                  w
                  ⁢
                  
                    
                  
                
              
            
            ⁡
            
              [
              
                
                  ln
                  ⁡
                  
                    (
                    
                      
                        h
                        w
                      
                      a
                    
                    )
                  
                
                -
                1
              
              ]
            
          
        
        ,
      
    
    
      
        (
        50
        )
      
    
  




where Vw≈0.94 for a uniform cylindrical conductor and a is the radius of the conductor. An alternative expression that has been employed in amateur radio literature for the characteristic impedance of a single-wire feed line can be given by





  
    
      
        
          Z
          w
        
        =
        
          138
          ⁢
          
            
          
          ⁢
          
            
              log
              ⁡
              
                (
                
                  
                    1.123
                    ⁢
                    
                      
                    
                    ⁢
                    
                      V
                      w
                    
                    ⁢
                    
                      λ
                      0
                    
                  
                  
                    2
                    ⁢
                    π
                    ⁢
                    
                      
                    
                    ⁢
                    a
                  
                
                )
              
            
            .
          
        
      
    
    
      
        (
        51
        )
      
    
  




Equation (51) implies that Zw for a single-wire feeder varies with frequency. The phase delay can be determined based upon the capacitance and characteristic impedance.

With a charge terminal T1 positioned over the lossy conducting medium 203 as shown in FIG. 3, the feed network 209 can be adjusted to excite the charge terminal T1 with the phase shift (Φ) of the complex effective height (heff) equal to the angle (Ψ) of the wave tilt at the Hankel crossover distance, or Φ=Ψ. When this condition is met, the electric field produced by the charge oscillating Q1 on the charge terminal T1 is coupled into a guided surface waveguide mode traveling along the surface of a lossy conducting medium 203. For example, if the Brewster angle (θi,B), the phase delay (θy) associated with the vertical feed line conductor 221 (FIG. 7), and the configuration of the coil 215 (FIG. 7) are known, then the position of the tap 224 (FIG. 7) can be determined and adjusted to impose an oscillating charge Q1 on the charge terminal T1 with phase Φ=Ψ. The position of the tap 224 may be adjusted to maximize coupling the traveling surface waves into the guided surface waveguide mode. Excess coil length beyond the position of the tap 224 can be removed to reduce the capacitive effects. The vertical wire height and/or the geometrical parameters of the helical coil may also be varied.
The coupling to the guided surface waveguide mode on the surface of the lossy conducting medium 203 can be improved and/or optimized by tuning the guided surface waveguide probe 200 for standing wave resonance with respect to a complex image plane associated with the charge Q1 on the charge terminal T1. By doing this, the performance of the guided surface waveguide probe 200 can be adjusted for increased and/or maximum voltage (and thus charge Q1) on the charge terminal T1. Referring back to FIG. 3, the effect of the lossy conducting medium 203 in Region 1 can be examined using image theory analysis.
Physically, an elevated charge Q1 placed over a perfectly conducting plane attracts the free charge on the perfectly conducting plane, which then “piles up” in the region under the elevated charge Q1. The resulting distribution of “bound” electricity on the perfectly conducting plane is similar to a bell-shaped curve. The superposition of the potential of the elevated charge Q1, plus the potential of the induced “piled up” charge beneath it, forces a zero equipotential surface for the perfectly conducting plane. The boundary value problem solution that describes the fields in the region above the perfectly conducting plane may be obtained using the classical notion of image charges, where the field from the elevated charge is superimposed with the field from a corresponding “image” charge below the perfectly conducting plane.
This analysis may also be used with respect to a lossy conducting medium 203 by assuming the presence of an effective image charge Q1′ beneath the guided surface waveguide probe 200. The effective image charge Q1′ coincides with the charge Q1 on the charge terminal T1 about a conducting image ground plane 130, as illustrated in FIG. 3. However, the image charge Q1′ is not merely located at some real depth and 180° out of phase with the primary source charge Q1 on the charge terminal T1, as they would be in the case of a perfect conductor. Rather, the lossy conducting medium 203 (e.g., a terrestrial medium) presents a phase shifted image. That is to say, the image charge Q1′ is at a complex depth below the surface (or physical boundary) of the lossy conducting medium 203. For a discussion of complex image depth, reference is made to Wait, J. R., “Complex Image Theory—Revisited,” IEEE Antennas and Propagation Magazine, Vol. 33, No. 4, August 1991, pp. 27-29, which is incorporated herein by reference in its entirety.
Instead of the image charge Q1′ being at a depth that is equal to the physical height (H1) of the charge Q1, the conducting image ground plane 130 (representing a perfect conductor) is located at a complex depth of z=−d/2 and the image charge Q1′ appears at a complex depth (i.e., the “depth” has both magnitude and phase), given by −D1=−(d/2+d/2+H1)≠H1. For vertically polarized sources over the Earth,




  
    
      
        
          d
          =
          
            
              
                
                  
                    
                      γ
                      e
                      2
                    
                    +
                    
                      k
                      0
                      2
                    
                  
                  2
                
                
                  γ
                  e
                  2
                
              
              ≈
              
                2
                
                  γ
                  e
                
              
            
            =
            
              
                
                  d
                  r
                
                +
                
                  j
                  ⁢
                  
                    
                  
                  ⁢
                  
                    d
                    i
                  
                
              
              =
              
                
                  
                  d
                  
                
                ⁢
                ∠ζ
              
            
          
        
        ,
        
          

        
        ⁢
        where
      
    
    
      
        (
        52
        )
      
    
  
  
    
      
        
          
            γ
            e
            2
          
          =
          
            
              j
              ⁢
              
                
              
              ⁢
              
                ωμ
                1
              
              ⁢
              
                σ
                1
              
            
            -
            
              
                ω
                2
              
              ⁢
              
                μ
                1
              
              ⁢
              
                ɛ
                1
              
            
          
        
        ,
        and
      
    
    
      
        (
        53
        )
      
    
  
  
    
      
        
          
            k
            o
          
          =
          
            ω
            ⁢
            
              
                
                  μ
                  o
                
                ⁢
                
                  ɛ
                  o
                
              
            
          
        
        ,
      
    
    
      
        (
        54
        )
      
    
  




as indicated in Equation (12). The complex spacing of the image charge, in turn, implies that the external field will experience extra phase shifts not encountered when the interface is either a dielectric or a perfect conductor. In the lossy conducting medium, the wave front normal is parallel to the tangent of the conducting image ground plane 130 at z=−d/2, and not at the boundary interface between Regions 1 and 2.

Consider the case illustrated in FIG. 8A where the lossy conducting medium 203 is a finitely conducting Earth 133 with a physical boundary 136. The finitely conducting Earth 133 may be replaced by a perfectly conducting image ground plane 139 as shown in FIG. 8B, which is located at a complex depth z1 below the physical boundary 136. This equivalent representation exhibits the same impedance when looking down into the interface at the physical boundary 136. The equivalent representation of FIG. 8B can be modeled as an equivalent transmission line, as shown in FIG. 8C. The cross-section of the equivalent structure is represented as a (z-directed) end-loaded transmission line, with the impedance of the perfectly conducting image plane being a short circuit (zS=0). The depth z1 can be determined by equating the TEM wave impedance looking down at the Earth to an image ground plane impedance zin seen looking into the transmission line of FIG. 8C.
In the case of FIG. 8A, the propagation constant and wave intrinsic impedance in the upper region (air) 142 are




  
    
      
        
          
            γ
            o
          
          =
          
            
              j
              ⁢
              
                
              
              ⁢
              ω
              ⁢
              
                
                  
                    μ
                    o
                  
                  ⁢
                  
                    ɛ
                    o
                  
                
              
            
            =
            
              0
              +
              
                j
                ⁢
                
                  
                
                ⁢
                β
              
            
          
        
        ,
        and
      
    
    
      
        (
        55
        )
      
    
  
  
    
      
        
          
            z
            o
          
          =
          
            
              
                j
                ⁢
                
                  
                
                ⁢
                
                  ωμ
                  o
                
              
              
                γ
                o
              
            
            =
            
              
                
                  μ
                  o
                
                
                  ɛ
                  o
                
              
            
          
        
        ,
      
    
    
      
        (
        56
        )
      
    
  




In the lossy Earth 133, the propagation constant and wave intrinsic impedance are





  
    
      
        
          
            γ
            e
          
          =
          
            
              j
              ⁢
              
                
              
              ⁢
              
                
                  ωμ
                  1
                
                ⁡
                
                  (
                  
                    
                      σ
                      1
                    
                    +
                    
                      j
                      ⁢
                      
                        
                      
                      ⁢
                      
                        ωɛ
                        1
                      
                    
                  
                  )
                
              
            
          
        
        ,
        and
      
    
    
      
        (
        57
        )
      
    
  
  
    
      
        
          
            Z
            e
          
          =
          
            
              j
              ⁢
              
                
              
              ⁢
              
                ωμ
                1
              
            
            
              γ
              
                e
                ⁢
                
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        58
        )
      
    
  




For normal incidence, the equivalent representation of FIG. 8B is equivalent to a TEM transmission line whose characteristic impedance is that of air (zo), with propagation constant of γo, and whose length is z1. As such, the image ground plane impedance Zin seen at the interface for the shorted transmission line of FIG. 8C is given by

Zin=Z0 tan h(γ0z1).  (59)

Equating the image ground plane impedance Zin associated with the equivalent model of FIG. 8C to the normal incidence wave impedance of FIG. 8A and solving for z1 gives the distance to a short circuit (the perfectly conducting image ground plane 139) as





  
    
      
        
          
            z
            1
          
          =
          
            
              
                1
                
                  γ
                  o
                
              
              ⁢
              
                
                  tanh
                  
                    -
                    1
                  
                
                ⁡
                
                  (
                  
                    
                      z
                      e
                    
                    
                      z
                      o
                    
                  
                  )
                
              
            
            =
            
              
                
                  1
                  
                    γ
                    o
                  
                
                ⁢
                
                  
                    tanh
                    
                      -
                      1
                    
                  
                  ⁡
                  
                    (
                    
                      
                        γ
                        o
                      
                      
                        γ
                        e
                      
                    
                    )
                  
                
              
              ≈
              
                1
                
                  γ
                  e
                
              
            
          
        
        ,
      
    
    
      
        (
        60
        )
      
    
  




where only the first term of the series expansion for the inverse hyperbolic tangent is considered for this approximation. Note that in the air region 142, the propagation constant is γo=jβo, so Zin=jZo tan βoz1 (which is a purely imaginary quantity for a real z1), but ze is a complex value if σ≠0. Therefore, Zin=Ze only when z1 is a complex distance.

Since the equivalent representation of FIG. 8B includes a perfectly conducting image ground plane 139, the image depth for a charge or current lying at the surface of the Earth (physical boundary 136) is equal to distance z1 on the other side of the image ground plane 139, or d=2×z1 beneath the Earth's surface (which is located at z=0). Thus, the distance to the perfectly conducting image ground plane 139 can be approximated by




  
    
      
        d
        =
        
          
            2
            ⁢
            
              z
              1
            
          
          ≈
          
            
              2
              
                γ
                e
              
            
            .
          
        
      
    
    
      
        (
        61
        )
      
    
  




Additionally, the “image charge” will be “equal and opposite” to the real charge, so the potential of the perfectly conducting image ground plane 139 at depth z1=−d/2 will be zero.

If a charge Q1 is elevated a distance H1 above the surface of the Earth as illustrated in FIG. 3, then the image charge Q1′ resides at a complex distance of D1=d+H1 below the surface, or a complex distance of d/2+H1 below the image ground plane 130. The guided surface waveguide probe 200b of FIG. 7 can be modeled as an equivalent single-wire transmission line image plane model that can be based upon the perfectly conducting image ground plane 139 of FIG. 8B. FIG. 9A shows an example of the equivalent single-wire transmission line image plane model, and FIG. 9B illustrates an example of the equivalent classic transmission line model, including the shorted transmission line of FIG. 8C.
In the equivalent image plane models of FIGS. 9A and 9B, Φ=θy+θc is the traveling wave phase delay of the guided surface waveguide probe 200 referenced to Earth 133 (or the lossy conducting medium 203), θc=βpH is the electrical length of the coil 215 (FIG. 7), of physical length H, expressed in degrees, θy=βwhw is the electrical length of the vertical feed line conductor 221 (FIG. 7), of physical length hw, expressed in degrees, and θd=βo d/2 is the phase shift between the image ground plane 139 and the physical boundary 136 of the Earth 133 (or lossy conducting medium 203). In the example of FIGS. 9A and 9B, Zw is the characteristic impedance of the elevated vertical feed line conductor 221 in ohms, Zc is the characteristic impedance of the coil 215 in ohms, and ZO is the characteristic impedance of free space.
At the base of the guided surface waveguide probe 200, the impedance seen “looking up” into the structure is Z↑=Zbase. With a load impedance of:




  
    
      
        
          
            Z
            L
          
          =
          
            1
            
              j
              ⁢
              
                
              
              ⁢
              ω
              ⁢
              
                
              
              ⁢
              
                C
                T
              
            
          
        
        ,
      
    
    
      
        (
        62
        )
      
    
  




where CT is the self-capacitance of the charge terminal T1, the impedance seen “looking up” into the vertical feed line conductor 221 (FIG. 7) is given by:





  
    
      
        
          
            Z
            2
          
          =
          
            
              
                Z
                W
              
              ⁢
              
                
                  
                    Z
                    L
                  
                  +
                  
                    
                      Z
                      W
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            w
                          
                          ⁢
                          
                            h
                            w
                          
                        
                        )
                      
                    
                  
                
                
                  
                    Z
                    W
                  
                  +
                  
                    
                      Z
                      L
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            w
                          
                          ⁢
                          
                            h
                            w
                          
                        
                        )
                      
                    
                  
                
              
            
            =
            
              
                Z
                W
              
              ⁢
              
                
                  
                    Z
                    L
                  
                  +
                  
                    
                      Z
                      W
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            θ
                            y
                          
                        
                        )
                      
                    
                  
                
                
                  
                    Z
                    W
                  
                  +
                  
                    
                      Z
                      L
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            θ
                            y
                          
                        
                        )
                      
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        63
        )
      
    
  




and the impedance seen “looking up” into the coil 215 (FIG. 7) is given by:





  
    
      
        
          Z
          base
        
        =
        
          
            
              Z
              c
            
            ⁢
            
              
                
                  Z
                  2
                
                +
                
                  
                    Z
                    c
                  
                  ⁢
                  
                    tanh
                    ⁡
                    
                      (
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        
                          β
                          p
                        
                        ⁢
                        H
                      
                      )
                    
                  
                
              
              
                
                  Z
                  c
                
                +
                
                  
                    Z
                    2
                  
                  ⁢
                  
                    tanh
                    ⁡
                    
                      (
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        
                          β
                          p
                        
                        ⁢
                        H
                      
                      )
                    
                  
                
              
            
          
          =
          
            
              Z
              c
            
            ⁢
            
              
                
                  
                    Z
                    2
                  
                  +
                  
                    
                      Z
                      c
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            θ
                            c
                          
                        
                        )
                      
                    
                  
                
                
                  
                    Z
                    c
                  
                  +
                  
                    
                      Z
                      2
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            θ
                            c
                          
                        
                        )
                      
                    
                  
                
              
              .
            
          
        
      
    
    
      
        (
        64
        )
      
    
  




At the base of the guided surface waveguide probe 200, the impedance seen “looking down” into the lossy conducting medium 203 is Z↓=Zin, which is given by:





  
    
      
        
          
            Z
            in
          
          =
          
            
              
                Z
                o
              
              ⁢
              
                
                  
                    Z
                    s
                  
                  +
                  
                    
                      Z
                      o
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        [
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            
                              β
                              o
                            
                            ⁡
                            
                              (
                              
                                d
                                ⁢
                                
                                  /
                                
                                ⁢
                                2
                              
                              )
                            
                          
                        
                        ]
                      
                    
                  
                
                
                  
                    Z
                    o
                  
                  +
                  
                    
                      Z
                      s
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        [
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            
                              β
                              o
                            
                            ⁡
                            
                              (
                              
                                d
                                ⁢
                                
                                  /
                                
                                ⁢
                                2
                              
                              )
                            
                          
                        
                        ]
                      
                    
                  
                
              
            
            =
            
              
                Z
                o
              
              ⁢
              
                tanh
                ⁡
                
                  (
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    
                      θ
                      
                        d
                        ⁢
                        
                      
                    
                  
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        65
        )
      
    
  




where ZS=0.

Neglecting losses, the equivalent image plane model can be tuned to resonance when Z↓+Z↑=0 at the physical boundary 136. Or, in the low loss case, X↓+X↑=0 at the physical boundary 136, where X is the corresponding reactive component. Thus, the impedance at the physical boundary 136 “looking up” into the guided surface waveguide probe 200 is the conjugate of the impedance at the physical boundary 136 “looking down” into the lossy conducting medium 203. By adjusting the load impedance ZL of the charge terminal T1 while maintaining the traveling wave phase delay Φ equal to the angle of the media's wave tilt Ψ, so that Φ=Ψ, which improves and/or maximizes coupling of the probe's electric field to a guided surface waveguide mode along the surface of the lossy conducting medium 203 (e.g., Earth), the equivalent image plane models of FIGS. 9A and 9B can be tuned to resonance with respect to the image ground plane 139. In this way, the impedance of the equivalent complex image plane model is purely resistive, which maintains a superposed standing wave on the probe structure that maximizes the voltage and elevated charge on terminal T1, and by equations (1)-(3) and (16) maximizes the propagating surface wave.
It follows from the Hankel solutions, that the guided surface wave excited by the guided surface waveguide probe 200 is an outward propagating traveling wave. The source distribution along the feed network 209 between the charge terminal T1 and the ground stake 218 of the guided surface waveguide probe 200 (FIGS. 3 and 7) is actually composed of a superposition of a traveling wave plus a standing wave on the structure. With the charge terminal T1 positioned at or above the physical height hp, the phase delay of the traveling wave moving through the feed network 209 is matched to the angle of the wave tilt associated with the lossy conducting medium 203. This mode-matching allows the traveling wave to be launched along the lossy conducting medium 203. Once the phase delay has been established for the traveling wave, the load impedance ZL of the charge terminal T1 is adjusted to bring the probe structure into standing wave resonance with respect to the image ground plane (130 of FIG. 3 or 139 of FIG. 8), which is at a complex depth of −d/2. In that case, the impedance seen from the image ground plane has zero reactance and the charge on the charge terminal T1 is maximized.
The distinction between the traveling wave phenomenon and standing wave phenomena is that (1) the phase delay of traveling waves (θ=βd) on a section of transmission line of length d (sometimes called a “delay line”) is due to propagation time delays; whereas (2) the position-dependent phase of standing waves (which are composed of forward and backward propagating waves) depends on both the line length propagation time delay and impedance transitions at interfaces between line sections of different characteristic impedances. In addition to the phase delay that arises due to the physical length of a section of transmission line operating in sinusoidal steady-state, there is an extra reflection coefficient phase at impedance discontinuities that is due to the ratio of Zoa/Zob, where Zoa and Zob are the characteristic impedances of two sections of a transmission line such as, e.g., a helical coil section of characteristic impedance Zoa=Zc (FIG. 9B) and a straight section of vertical feed line conductor of characteristic impedance Zob=Zw (FIG. 9B).
As a result of this phenomenon, two relatively short transmission line sections of widely differing characteristic impedance may be used to provide a very large phase shift. For example, a probe structure composed of two sections of transmission line, one of low impedance and one of high impedance, together totaling a physical length of, say, 0.05λ, may be fabricated to provide a phase shift of 90° which is equivalent to a 0.25λ resonance. This is due to the large jump in characteristic impedances. In this way, a physically short probe structure can be electrically longer than the two physical lengths combined. This is illustrated in FIGS. 9A and 9B, where the discontinuities in the impedance ratios provide large jumps in phase. The impedance discontinuity provides a substantial phase shift where the sections are joined together.
Referring to FIG. 10, shown is a flow chart 150 illustrating an example of adjusting a guided surface waveguide probe 200 (FIGS. 3 and 7) to substantially mode-match to a guided surface waveguide mode on the surface of the lossy conducting medium, which launches a guided surface traveling wave along the surface of a lossy conducting medium 203 (FIG. 3). Beginning with 153, the charge terminal T1 of the guided surface waveguide probe 200 is positioned at a defined height above a lossy conducting medium 203. Utilizing the characteristics of the lossy conducting medium 203 and the operating frequency of the guided surface waveguide probe 200, the Hankel crossover distance can also be found by equating the magnitudes of Equations (20b) and (21) for −jγρ, and solving for Rx as illustrated by FIG. 4. The complex index of refraction (n) can be determined using Equation (41), and the complex Brewster angle (θi,B) can then be determined from Equation (42). The physical height (hp) of the charge terminal T1 can then be determined from Equation (44). The charge terminal T1 should be at or higher than the physical height (hp) in order to excite the far-out component of the Hankel function. This height relationship is initially considered when launching surface waves. To reduce or minimize the bound charge on the charge terminal T1, the height should be at least four times the spherical diameter (or equivalent spherical diameter) of the charge terminal T1.
At 156, the electrical phase delay Φ of the elevated charge Q1 on the charge terminal T1 is matched to the complex wave tilt angle Ψ. The phase delay (θc) of the helical coil and/or the phase delay (θy) of the vertical feed line conductor can be adjusted to make Φ equal to the angle (Ψ) of the wave tilt (W). Based on Equation (31), the angle (Ψ) of the wave tilt can be determined from:




  
    
      
        W
        =
        
          
            
              E
              ρ
            
            
              E
              z
            
          
          =
          
            
              1
              
                tan
                ⁢
                
                  
                
                ⁢
                
                  θ
                  
                    i
                    ,
                    B
                  
                
              
            
            =
            
              
                1
                n
              
              =
              
                
                  
                  W
                  
                
                ⁢
                
                  
                    e
                    
                      j
                      ⁢
                      
                        
                      
                      ⁢
                      Ψ
                    
                  
                  .
                
              
            
          
        
      
    
    
      
        (
        66
        )
      
    
  




The electrical phase Φ can then be matched to the angle of the wave tilt. This angular (or phase) relationship is next considered when launching surface waves. For example, the electrical phase delay Φ=θc+θy can be adjusted by varying the geometrical parameters of the coil 215 (FIG. 7) and/or the length (or height) of the vertical feed line conductor 221 (FIG. 7). By matching Φ=Ψ, an electric field can be established at or beyond the Hankel crossover distance (Rx) with a complex Brewster angle at the boundary interface to excite the surface waveguide mode and launch a traveling wave along the lossy conducting medium 203.

Next at 159, the load impedance of the charge terminal T1 is tuned to resonate the equivalent image plane model of the guided surface waveguide probe 200. The depth (d/2) of the conducting image ground plane 139 of FIGS. 9A and 9B (or 130 of FIG. 3) can be determined using Equations (52), (53) and (54) and the values of the lossy conducting medium 203 (e.g., the Earth), which can be measured. Using that depth, the phase shift (θd) between the image ground plane 139 and the physical boundary 136 of the lossy conducting medium 203 can be determined using θd=βo d/2. The impedance (Zin) as seen “looking down” into the lossy conducting medium 203 can then be determined using Equation (65). This resonance relationship can be considered to maximize the launched surface waves.
Based upon the adjusted parameters of the coil 215 and the length of the vertical feed line conductor 221, the velocity factor, phase delay, and impedance of the coil 215 and vertical feed line conductor 221 can be determined using Equations (45) through (51). In addition, the self-capacitance (CT) of the charge terminal T1 can be determined using, e.g., Equation (24). The propagation factor (βp) of the coil 215 can be determined using Equation (35) and the propagation phase constant (βp) for the vertical feed line conductor 221 can be determined using Equation (49). Using the self-capacitance and the determined values of the coil 215 and vertical feed line conductor 221, the impedance (Zbase) of the guided surface waveguide probe 200 as seen “looking up” into the coil 215 can be determined using Equations (62), (63) and (64).
The equivalent image plane model of the guided surface waveguide probe 200 can be tuned to resonance by adjusting the load impedance ZL such that the reactance component Xbase of Zbase cancels out the reactance component Xin of Zin, or Xbase+Xin=0. Thus, the impedance at the physical boundary 136 “looking up” into the guided surface waveguide probe 200 is the conjugate of the impedance at the physical boundary 136 “looking down” into the lossy conducting medium 203. The load impedance ZL can be adjusted by varying the capacitance (CT) of the charge terminal T1 without changing the electrical phase delay Φ=θc+θy of the charge terminal T1. An iterative approach may be taken to tune the load impedance ZL for resonance of the equivalent image plane model with respect to the conducting image ground plane 139 (or 130). In this way, the coupling of the electric field to a guided surface waveguide mode along the surface of the lossy conducting medium 203 (e.g., Earth) can be improved and/or maximized.
This may be better understood by illustrating the situation with a numerical example. Consider a guided surface waveguide probe 200 comprising a top-loaded vertical stub of physical height hp with a charge terminal T1 at the top, where the charge terminal T1 is excited through a helical coil and vertical feed line conductor at an operational frequency (fo) of 1.85 MHz. With a height (H1) of 16 feet and the lossy conducting medium 203 (e.g., Earth) having a relative permittivity of ∈r=15 and a conductivity of σ1=0.010 mhos/m, several surface wave propagation parameters can be calculated for fo=1.850 MHz. Under these conditions, the Hankel crossover distance can be found to be Rx=54.5 feet with a physical height of hp=5.5 feet, which is well below the actual height of the charge terminal T1. While a charge terminal height of H1=5.5 feet could have been used, the taller probe structure reduced the bound capacitance, permitting a greater percentage of free charge on the charge terminal T1 providing greater field strength and excitation of the traveling wave.
The wave length can be determined as:




  
    
      
        
          
            λ
            o
          
          =
          
            
              c
              
                f
                o
              
            
            =
            
              162.162
              ⁢
              
                
              
              ⁢
              meters
            
          
        
        ,
      
    
    
      
        (
        67
        )
      
    
  




where c is the speed of light. The complex index of refraction is:

n=√{square root over (∈r−jx)}=7.529−j 6.546,  (68)

from Equation (41), where x=σ1/ω∈o with ω=2πfo, and the complex Brewster angle is:

θi,B=arctan(√{square root over (∈r−jx)})=85.6−j 3.744°.  (69)

from Equation (42). Using Equation (66), the wave tilt values can be determined to be:





  
    
      
        W
        =
        
          
            1
            
              tan
              ⁢
              
                
              
              ⁢
              
                θ
                
                  i
                  ,
                  B
                
              
            
          
          =
          
            
              1
              n
            
            =
            
              
                
                  
                  W
                  
                
                ⁢
                
                  e
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    Ψ
                  
                
              
              =
              
                0.101
                ⁢
                
                  
                
                ⁢
                
                  e
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        40.614
                        °
                      
                      .
                    
                  
                
              
            
          
        
      
    
    
      
        (
        70
        )
      
    
  




Thus, the helical coil can be adjusted to match Φ=Ψ=40.614°

The velocity factor of the vertical feed line conductor (approximated as a uniform cylindrical conductor with a diameter of 0.27 inches) can be given as Vw≈0.93. Since hp<<λo, the propagation phase constant for the vertical feed line conductor can be approximated as:




  
    
      
        
          
            β
            w
          
          =
          
            
              
                2
                ⁢
                π
              
              
                λ
                w
              
            
            =
            
              
                
                  2
                  ⁢
                  π
                
                
                  
                    V
                    w
                  
                  ⁢
                  
                    λ
                    0
                  
                
              
              =
              
                0.042
                ⁢
                
                  
                
                ⁢
                
                  
                    m
                    
                      -
                      1
                    
                  
                  .
                
              
            
          
        
        ⁢
        
          
        
      
    
    
      
        (
        71
        )
      
    
  




From Equation (49) the phase delay of the vertical feed line conductor is:

θy=βwhw≈βwhp=11.640°.  (72)

By adjusting the phase delay of the helical coil so that θc=28.974°=40.614°−11.640°, Φ will equal Ψ to match the guided surface waveguide mode. To illustrate the relationship between Φ and Ψ, FIG. 11 shows a plot of both over a range of frequencies. As both Φ and Ψ are frequency dependent, it can be seen that their respective curves cross over each other at approximately 1.85 MHz.

For a helical coil having a conductor diameter of 0.0881 inches, a coil diameter (D) of 30 inches and a turn-to-turn spacing (s) of 4 inches, the velocity factor for the coil can be determined using Equation (45) as:




  
    
      
        
          
            V
            f
          
          =
          
            
              1
              
                
                  1
                  +
                  
                    20
                    ⁢
                    
                      
                        (
                        
                          D
                          s
                        
                        )
                      
                      2.5
                    
                    ⁢
                    
                      
                        (
                        
                          D
                          
                            λ
                            o
                          
                        
                        )
                      
                      0.5
                    
                  
                
              
            
            =
            0.069
          
        
        ,
      
    
    
      
        (
        73
        )
      
    
  




and the propagation factor from Equation (35) is:





  
    
      
        
          
            β
            p
          
          =
          
            
              
                2
                ⁢
                π
              
              
                
                  V
                  f
                
                ⁢
                
                  λ
                  o
                
              
            
            =
            
              0.564
              ⁢
              
                
              
              ⁢
              
                
                  m
                  
                    -
                    1
                  
                
                .
              
            
          
        
        ⁢
        
          
        
      
    
    
      
        (
        74
        )
      
    
  




With θc=28.974°, the axial length of the solenoidal helix (H) can be determined using Equation (46) such that:





  
    
      
        H
        =
        
          
            
              θ
              c
            
            
              β
              p
            
          
          =
          
            35.2732
            ⁢
            
              
            
            ⁢
            
              inches
              .
            
          
        
      
    
    
      
        (
        75
        )
      
    
  




This height determines the location on the helical coil where the vertical feed line conductor is connected, resulting in a coil with 8.818 turns (N=H/s).

With the traveling wave phase delay of the coil and vertical feed line conductor adjusted to match the wave tilt angle (Φ=θc+θy=Ψ), the load impedance (ZL) of the charge terminal T1 can be adjusted for standing wave resonance of the equivalent image plane model of the guided surface wave probe 200. From the measured permittivity, conductivity and permeability of the Earth, the radial propagation constant can be determined using Equation (57)

γe=√{square root over (jωu1(σ1+jω∈1))}=0.25+j 0.292 m−1,  (76)

And the complex depth of the conducting image ground plane can be approximated from Equation (52) as:





  
    
      
        
          
            d
            ≈
            
              2
              
                γ
                e
              
            
          
          =
          
            3.364
            +
            
              j
              ⁢
              
                
              
              ⁢
              3.963
              ⁢
              
                
              
              ⁢
              meters
            
          
        
        ,
      
    
    
      
        (
        77
        )
      
    
  




with a corresponding phase shift between the conducting image ground plane and the physical boundary of the Earth given by:

θd=βo(d/2)=4.015−j 4.73°.  (78)

Using Equation (65), the impedance seen “looking down” into the lossy conducting medium 203 (i.e., Earth) can be determined as:

Zin=Zo tan h(jθd)=Rin+jXin=31.191+j 26.27 ohms.  (79)

By matching the reactive component (Xin) seen “looking down” into the lossy conducting medium 203 with the reactive component (Xbase) seen “looking up” into the guided surface wave probe 200, the coupling into the guided surface waveguide mode may be maximized. This can be accomplished by adjusting the capacitance of the charge terminal T1 without changing the traveling wave phase delays of the coil and vertical feed line conductor. For example, by adjusting the charge terminal capacitance (CT) to 61.8126 pF, the load impedance from Equation (62) is:




  
    
      
        
          
            Z
            L
          
          =
          
            
              1
              
                j
                ⁢
                
                  
                
                ⁢
                ω
                ⁢
                
                  
                
                ⁢
                
                  C
                  T
                
              
            
            =
            
              
                -
                j
              
              ⁢
              
                
              
              ⁢
              1392
              ⁢
              
                
              
              ⁢
              ohms
            
          
        
        ,
      
    
    
      
        (
        80
        )
      
    
  




and the reactive components at the boundary are matched.

Using Equation (51), the impedance of the vertical feed line conductor (having a diameter (2a) of 0.27 inches) is given as




  
    
      
        
          
            Z
            w
          
          =
          
            
              138
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    
                      1.123
                      ⁢
                      
                        
                      
                      ⁢
                      
                        V
                        w
                      
                      ⁢
                      
                        λ
                        0
                      
                    
                    
                      2
                      ⁢
                      π
                      ⁢
                      
                        
                      
                      ⁢
                      a
                    
                  
                  )
                
              
            
            =
            
              537.534
              ⁢
              
                
              
              ⁢
              ohms
            
          
        
        ,
      
    
    
      
        (
        81
        )
      
    
  




and the impedance seen “looking up” into the vertical feed line conductor is given by Equation (63) as:





  
    
      
        
          Z
          2
        
        =
        
          
            
              Z
              W
            
            ⁢
            
              
                
                  Z
                  L
                
                +
                
                  
                    Z
                    w
                  
                  ⁢
                  
                    tanh
                    ⁡
                    
                      (
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        
                          θ
                          y
                        
                      
                      )
                    
                  
                
              
              
                
                  Z
                  w
                
                +
                
                  
                    Z
                    L
                  
                  ⁢
                  
                    tanh
                    ⁡
                    
                      (
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        
                          θ
                          y
                        
                      
                      )
                    
                  
                
              
            
          
          =
          
            
              -
              j
            
            ⁢
            
              
            
            ⁢
            835.438
            ⁢
            
              
            
            ⁢
            
              ohms
              .
            
          
        
      
    
    
      
        (
        82
        )
      
    
  




Using Equation (47), the characteristic impedance of the helical coil is given as





  
    
      
        
          
            Z
            c
          
          =
          
            
              
                60
                
                  V
                  f
                
              
              ⁡
              
                [
                
                  
                    ln
                    ⁡
                    
                      (
                      
                        
                          
                            V
                            f
                          
                          ⁢
                          
                            λ
                            0
                          
                        
                        D
                      
                      )
                    
                  
                  -
                  1.027
                
                ]
              
            
            =
            
              1446
              ⁢
              
                
              
              ⁢
              ohms
            
          
        
        ,
      
    
    
      
        (
        83
        )
      
    
  




and the impedance seen “looking up” into the coil at the base is given by Equation (64) as:





  
    
      
        
          Z
          base
        
        =
        
          
            
              Z
              c
            
            ⁢
            
              
                
                  Z
                  2
                
                +
                
                  
                    Z
                    c
                  
                  ⁢
                  
                    tanh
                    ⁡
                    
                      (
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        
                          θ
                          c
                        
                      
                      )
                    
                  
                
              
              
                
                  Z
                  c
                
                +
                
                  
                    Z
                    2
                  
                  ⁢
                  
                    tanh
                    ⁡
                    
                      (
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        
                          θ
                          c
                        
                      
                      )
                    
                  
                
              
            
          
          =
          
            
              -
              j
            
            ⁢
            
              
            
            ⁢
            26.271
            ⁢
            
              
            
            ⁢
            
              ohms
              .
            
          
        
      
    
    
      
        (
        84
        )
      
    
  




When compared to the solution of Equation (79), it can be seen that the reactive components are opposite and approximately equal, and thus are conjugates of each other. Thus, the impedance (Zip) seen “looking up” into the equivalent image plane model of FIGS. 9A and 9B from the perfectly conducting image ground plane is only resistive or Zip=R+j0.

When the electric fields produced by a guided surface waveguide probe 200 (FIG. 3) are established by matching the traveling wave phase delay of the feed network to the wave tilt angle and the probe structure is resonated with respect to the perfectly conducting image ground plane at complex depth z=−d/2, the fields are substantially mode-matched to a guided surface waveguide mode on the surface of the lossy conducting medium, a guided surface traveling wave is launched along the surface of the lossy conducting medium. As illustrated in FIG. 1, the guided field strength curve 103 of the guided electromagnetic field has a characteristic exponential decay of e−ad/√{square root over (d)} and exhibits a distinctive knee 109 on the log-log scale.
In summary, both analytically and experimentally, the traveling wave component on the structure of the guided surface waveguide probe 200 has a phase delay (Φ) at its upper terminal that matches the angle (Ψ) of the wave tilt of the surface traveling wave (Φ=Ψ)=Under this condition, the surface waveguide may be considered to be “mode-matched”. Furthermore, the resonant standing wave component on the structure of the guided surface waveguide probe 200 has a VMAX at the charge terminal T1 and a VMIN down at the image plane 139 (FIG. 8B) where Zip=Rip+j 0 at a complex depth of z=−d/2, not at the connection at the physical boundary 136 of the lossy conducting medium 203 (FIG. 8B). Lastly, the charge terminal T1 is of sufficient height H1 of FIG. 3 (h≧Rx tan ψi,B) so that electromagnetic waves incident onto the lossy conducting medium 203 at the complex Brewster angle do so out at a distance (≧Rx) where the 1/√{square root over (r)} term is predominant. Receive circuits can be utilized with one or more guided surface waveguide probes to facilitate wireless transmission and/or power delivery systems.
Referring back to FIG. 3, operation of a guided surface waveguide probe 200 may be controlled to adjust for variations in operational conditions associated with the guided surface waveguide probe 200. For example, an adaptive probe control system 230 can be used to control the feed network 209 and/or the charge terminal T1 to control the operation of the guided surface waveguide probe 200. Operational conditions can include, but are not limited to, variations in the characteristics of the lossy conducting medium 203 (e.g., conductivity σ and relative permittivity ∈r), variations in field strength and/or variations in loading of the guided surface waveguide probe 200. As can be seen from Equations (31), (41) and (42), the index of refraction (n), the complex Brewster angle (θi,B), and the wave tilt (|W|ejΨ) can be affected by changes in soil conductivity and permittivity resulting from, e.g., weather conditions.
Equipment such as, e.g., conductivity measurement probes, permittivity sensors, ground parameter meters, field meters, current monitors and/or load receivers can be used to monitor for changes in the operational conditions and provide information about current operational conditions to the adaptive probe control system 230. The probe control system 230 can then make one or more adjustments to the guided surface waveguide probe 200 to maintain specified operational conditions for the guided surface waveguide probe 200. For instance, as the moisture and temperature vary, the conductivity of the soil will also vary. Conductivity measurement probes and/or permittivity sensors may be located at multiple locations around the guided surface waveguide probe 200. Generally, it would be desirable to monitor the conductivity and/or permittivity at or about the Hankel crossover distance Rx for the operational frequency. Conductivity measurement probes and/or permittivity sensors may be located at multiple locations (e.g., in each quadrant) around the guided surface waveguide probe 200.
The conductivity measurement probes and/or permittivity sensors can be configured to evaluate the conductivity and/or permittivity on a periodic basis and communicate the information to the probe control system 230. The information may be communicated to the probe control system 230 through a network such as, but not limited to, a LAN, WLAN, cellular network, or other appropriate wired or wireless communication network. Based upon the monitored conductivity and/or permittivity, the probe control system 230 may evaluate the variation in the index of refraction (n), the complex Brewster angle (θi,B), and/or the wave tilt (|W|ejΨ) and adjust the guided surface waveguide probe 200 to maintain the phase delay (Φ) of the feed network 209 equal to the wave tilt angle (Ψ) and/or maintain resonance of the equivalent image plane model of the guided surface waveguide probe 200. This can be accomplished by adjusting, e.g., θy, θc and/or CT. For instance, the probe control system 230 can adjust the self-capacitance of the charge terminal T1 and/or the phase delay (θy, θc) applied to the charge terminal T1 to maintain the electrical launching efficiency of the guided surface wave at or near its maximum. For example, the self-capacitance of the charge terminal T1 can be varied by changing the size of the terminal. The charge distribution can also be improved by increasing the size of the charge terminal T1, which can reduce the chance of an electrical discharge from the charge terminal T1. In other embodiments, the charge terminal T1 can include a variable inductance that can be adjusted to change the load impedance ZL. The phase applied to the charge terminal T1 can be adjusted by varying the tap position on the coil 215 (FIG. 7), and/or by including a plurality of predefined taps along the coil 215 and switching between the different predefined tap locations to maximize the launching efficiency.
Field or field strength (FS) meters may also be distributed about the guided surface waveguide probe 200 to measure field strength of fields associated with the guided surface wave. The field or FS meters can be configured to detect the field strength and/or changes in the field strength (e.g., electric field strength) and communicate that information to the probe control system 230. The information may be communicated to the probe control system 230 through a network such as, but not limited to, a LAN, WLAN, cellular network, or other appropriate communication network. As the load and/or environmental conditions change or vary during operation, the guided surface waveguide probe 200 may be adjusted to maintain specified field strength(s) at the FS meter locations to ensure appropriate power transmission to the receivers and the loads they supply.
For example, the phase delay (Φ=θy+θc) applied to the charge terminal T1 can be adjusted to match the wave tilt angle (Ψ). By adjusting one or both phase delays, the guided surface waveguide probe 200 can be adjusted to ensure the wave tilt corresponds to the complex Brewster angle. This can be accomplished by adjusting a tap position on the coil 215 (FIG. 7) to change the phase delay supplied to the charge terminal T1. The voltage level supplied to the charge terminal T1 can also be increased or decreased to adjust the electric field strength. This may be accomplished by adjusting the output voltage of the excitation source 212 or by adjusting or reconfiguring the feed network 209. For instance, the position of the tap 227 (FIG. 7) for the AC source 212 can be adjusted to increase the voltage seen by the charge terminal T1. Maintaining field strength levels within predefined ranges can improve coupling by the receivers, reduce ground current losses, and avoid interference with transmissions from other guided surface waveguide probes 200.
The probe control system 230 can be implemented with hardware, firmware, software executed by hardware, or a combination thereof. For example, the probe control system 230 can include processing circuitry including a processor and a memory, both of which can be coupled to a local interface such as, for example, a data bus with an accompanying control/address bus as can be appreciated by those with ordinary skill in the art. A probe control application may be executed by the processor to adjust the operation of the guided surface waveguide probe 200 based upon monitored conditions. The probe control system 230 can also include one or more network interfaces for communicating with the various monitoring devices. Communications can be through a network such as, but not limited to, a LAN, WLAN, cellular network, or other appropriate communication network. The probe control system 230 may comprise, for example, a computer system such as a server, desktop computer, laptop, or other system with like capability.
Referring back to the example of FIG. 5A, the complex angle trigonometry is shown for the ray optic interpretation of the incident electric field (E) of the charge terminal T1 with a complex Brewster angle (θi,B) at the Hankel crossover distance (Rx). Recall that, for a lossy conducting medium, the Brewster angle is complex and specified by equation (38). Electrically, the geometric parameters are related by the electrical effective height (heff) of the charge terminal T1 by equation (39). Since both the physical height (hp) and the Hankel crossover distance (Rx) are real quantities, the angle of the desired guided surface wave tilt at the Hankel crossover distance (WRx) is equal to the phase (Φ) of the complex effective height (heff). With the charge terminal T1 positioned at the physical height hp and excited with a charge having the appropriate phase Φ, the resulting electric field is incident with the lossy conducting medium boundary interface at the Hankel crossover distance Rx, and at the Brewster angle. Under these conditions, the guided surface waveguide mode can be excited without reflection or substantially negligible reflection.
However, Equation (39) means that the physical height of the guided surface waveguide probe 200 can be relatively small. While this will excite the guided surface waveguide mode, this can result in an unduly large bound charge with little free charge. To compensate, the charge terminal T1 can be raised to an appropriate elevation to increase the amount of free charge. As one example rule of thumb, the charge terminal T1 can be positioned at an elevation of about 4-5 times (or more) the effective diameter of the charge terminal T1. FIG. 6 illustrates the effect of raising the charge terminal T1 above the physical height (hp) shown in FIG. 5A. The increased elevation causes the distance at which the wave tilt is incident with the lossy conductive medium to move beyond the Hankel crossover point 121 (FIG. 5A). To improve coupling in the guided surface waveguide mode, and thus provide for a greater launching efficiency of the guided surface wave, a lower compensation terminal T2 can be used to adjust the total effective height (hTE) of the charge terminal T1 such that the wave tilt at the Hankel crossover distance is at the Brewster angle.
Referring to FIG. 12, shown is an example of a guided surface waveguide probe 200c that includes an elevated charge terminal T1 and a lower compensation terminal T2 that are arranged along a vertical axis z that is normal to a plane presented by the lossy conducting medium 203. In this respect, the charge terminal T1 is placed directly above the compensation terminal T2 although it is possible that some other arrangement of two or more charge and/or compensation terminals TN can be used. The guided surface waveguide probe 200c is disposed above a lossy conducting medium 203 according to an embodiment of the present disclosure. The lossy conducting medium 203 makes up Region 1 with a second medium 206 that makes up Region 2 sharing a boundary interface with the lossy conducting medium 203.
The guided surface waveguide probe 200c includes a coupling circuit 209 that couples an excitation source 212 to the charge terminal T1 and the compensation terminal T2. According to various embodiments, charges Q1 and Q2 can be imposed on the respective charge and compensation terminals T1 and T2, depending on the voltages applied to terminals T1 and T2 at any given instant. I1 is the conduction current feeding the charge Q1 on the charge terminal T1 via the terminal lead, and I2 is the conduction current feeding the charge Q2 on the compensation terminal T2 via the terminal lead.
According to the embodiment of FIG. 12, the charge terminal T1 is positioned over the lossy conducting medium 203 at a physical height H1, and the compensation terminal T2 is positioned directly below T1 along the vertical axis z at a physical height H2, where H2 is less than H1. The height h of the transmission structure may be calculated as h=H1−H2. The charge terminal T1 has an isolated (or self) capacitance C1, and the compensation terminal T2 has an isolated (or self) capacitance C2. A mutual capacitance CM can also exist between the terminals T1 and T2 depending on the distance therebetween. During operation, charges Q1 and Q2 are imposed on the charge terminal T1 and the compensation terminal T2, respectively, depending on the voltages applied to the charge terminal T1 and the compensation terminal T2 at any given instant.
Referring next to FIG. 13, shown is a ray optics interpretation of the effects produced by the elevated charge Q1 on charge terminal T1 and compensation terminal T2 of FIG. 12. With the charge terminal T1 elevated to a height where the ray intersects with the lossy conductive medium at the Brewster angle at a distance greater than the Hankel crossover point 121 as illustrated by line 163, the compensation terminal T2 can be used to adjust hTE by compensating for the increased height. The effect of the compensation terminal T2 is to reduce the electrical effective height of the guided surface waveguide probe (or effectively raise the lossy medium interface) such that the wave tilt at the Hankel crossover distance is at the Brewster angle as illustrated by line 166.
The total effective height can be written as the superposition of an upper effective height (hUE) associated with the charge terminal T1 and a lower effective height (hLE) associated with the compensation terminal T2 such that

hTE=hUE+hLE=hpej(βhp+ΦU)+hdej(βhd+ΦL)=Rx×W,  (85)

where ΦU is the phase delay applied to the upper charge terminal T1, ΦL is the phase delay applied to the lower compensation terminal T2, β=2π/λp is the propagation factor from Equation (35), hp is the physical height of the charge terminal T1 and hd is the physical height of the compensation terminal T2. If extra lead lengths are taken into consideration, they can be accounted for by adding the charge terminal lead length z to the physical height hp of the charge terminal T1 and the compensation terminal lead length y to the physical height hd of the compensation terminal T2 as shown in

hTE=(hp+z)ej(β(hp+z)+φU)+(hd+y)ej(β(hd+y)+ΦL)=Rx×W.  (86)

The lower effective height can be used to adjust the total effective height (hTE) to equal the complex effective height (heff) of FIG. 5A.

Equations (85) or (86) can be used to determine the physical height of the lower disk of the compensation terminal T2 and the phase angles to feed the terminals in order to obtain the desired wave tilt at the Hankel crossover distance. For example, Equation (86) can be rewritten as the phase shift applied to the charge terminal T1 as a function of the compensation terminal height (hd) to give




  
    
      
        
          
            Φ
            U
          
          ⁡
          
            (
            
              h
              d
            
            )
          
        
        =
        
          
            -
            
              β
              ⁡
              
                (
                
                  
                    h
                    p
                  
                  +
                  z
                
                )
              
            
          
          -
          
            j
            ⁢
            
              
            
            ⁢
            
              
                ln
                ⁡
                
                  (
                  
                    
                      
                        
                          R
                          x
                        
                        ×
                        w
                      
                      -
                      
                        
                          (
                          
                            
                              h
                              d
                            
                            +
                            y
                          
                          )
                        
                        ⁢
                        
                          e
                          
                            j
                            ⁢
                            
                              
                            
                            ⁢
                            
                              (
                              
                                
                                  β
                                  ⁢
                                  
                                    
                                  
                                  ⁢
                                  
                                    h
                                    d
                                  
                                
                                +
                                
                                  β
                                  ⁢
                                  
                                    
                                  
                                  ⁢
                                  y
                                
                                +
                                
                                  Φ
                                  L
                                
                              
                              )
                            
                          
                        
                      
                    
                    
                      (
                      
                        
                          h
                          p
                        
                        +
                        z
                      
                      )
                    
                  
                  )
                
              
              .
            
          
        
      
    
    
      
        (
        87
        )
      
    
  




To determine the positioning of the compensation terminal T2, the relationships discussed above can be utilized. First, the total effective height (hTE) is the superposition of the complex effective height (hUE) of the upper charge terminal T1 and the complex effective height (hLE) of the lower compensation terminal T2 as expressed in Equation (86). Next, the tangent of the angle of incidence can be expressed geometrically as




  
    
      
        
          
            tan
            ⁢
            
              
            
            ⁢
            
              ψ
              E
            
          
          =
          
            
              h
              
                T
                ⁢
                
                  
                
                ⁢
                E
              
            
            
              R
              x
            
          
        
        ,
      
    
    
      
        (
        88
        )
      
    
  




which is equal to the definition of the wave tilt, W. Finally, given the desired Hankel crossover distance Rx, the hTE can be adjusted to make the wave tilt of the incident ray match the complex Brewster angle at the Hankel crossover point 121. This can be accomplished by adjusting hp, ΦU, and/or hd.

These concepts may be better understood when discussed in the context of an example of a guided surface waveguide probe. Referring to FIG. 14, shown is a graphical representation of an example of a guided surface waveguide probe 200d including an upper charge terminal T1 (e.g., a sphere at height hT) and a lower compensation terminal T2 (e.g., a disk at height hd) that are positioned along a vertical axis z that is substantially normal to the plane presented by the lossy conducting medium 203. During operation, charges Q1 and Q2 are imposed on the charge and compensation terminals T1 and T2, respectively, depending on the voltages applied to the terminals T1 and T2 at any given instant.
An AC source 212 acts as the excitation source for the charge terminal T1, which is coupled to the guided surface waveguide probe 200d through a coupling circuit 209 comprising a coil 215 such as, e.g., a helical coil. The AC source 212 can be connected across a lower portion of the coil 215 through a tap 227, as shown in FIG. 14, or can be inductively coupled to the coil 215 by way of a primary coil. The coil 215 can be coupled to a ground stake 218 at a first end and the charge terminal T1 at a second end. In some implementations, the connection to the charge terminal T1 can be adjusted using a tap 224 at the second end of the coil 215. The compensation terminal T2 is positioned above and substantially parallel with the lossy conducting medium 203 (e.g., the ground or Earth), and energized through a tap 233 coupled to the coil 215. An ammeter 236 located between the coil 215 and ground stake 218 can be used to provide an indication of the magnitude of the current flow (I0) at the base of the guided surface waveguide probe. Alternatively, a current clamp may be used around the conductor coupled to the ground stake 218 to obtain an indication of the magnitude of the current flow (I0).
In the example of FIG. 14, the coil 215 is coupled to a ground stake 218 at a first end and the charge terminal T1 at a second end via a vertical feed line conductor 221. In some implementations, the connection to the charge terminal T1 can be adjusted using a tap 224 at the second end of the coil 215 as shown in FIG. 14. The coil 215 can be energized at an operating frequency by the AC source 212 through a tap 227 at a lower portion of the coil 215. In other implementations, the AC source 212 can be inductively coupled to the coil 215 through a primary coil. The compensation terminal T2 is energized through a tap 233 coupled to the coil 215. An ammeter 236 located between the coil 215 and ground stake 218 can be used to provide an indication of the magnitude of the current flow at the base of the guided surface waveguide probe 200d. Alternatively, a current clamp may be used around the conductor coupled to the ground stake 218 to obtain an indication of the magnitude of the current flow. The compensation terminal T2 is positioned above and substantially parallel with the lossy conducting medium 203 (e.g., the ground).
In the example of FIG. 14, the connection to the charge terminal T1 located on the coil 215 above the connection point of tap 233 for the compensation terminal T2. Such an adjustment allows an increased voltage (and thus a higher charge Q1) to be applied to the upper charge terminal T1. In other embodiments, the connection points for the charge terminal T1 and the compensation terminal T2 can be reversed. It is possible to adjust the total effective height (hTE) of the guided surface waveguide probe 200d to excite an electric field having a guided surface wave tilt at the Hankel crossover distance Rx. The Hankel crossover distance can also be found by equating the magnitudes of equations (20b) and (21) for −jγρ, and solving for Rx as illustrated by FIG. 4. The index of refraction (n), the complex Brewster angle (θi,B and ψi,B), the wave tilt (|W|ejΨ) and the complex effective height (heff=hpejΦ) can be determined as described with respect to Equations (41)-(44) above.
With the selected charge terminal T1 configuration, a spherical diameter (or the effective spherical diameter) can be determined. For example, if the charge terminal T1 is not configured as a sphere, then the terminal configuration may be modeled as a spherical capacitance having an effective spherical diameter. The size of the charge terminal T1 can be chosen to provide a sufficiently large surface for the charge Q1 imposed on the terminals. In general, it is desirable to make the charge terminal T1 as large as practical. The size of the charge terminal T1 should be large enough to avoid ionization of the surrounding air, which can result in electrical discharge or sparking around the charge terminal. To reduce the amount of bound charge on the charge terminal T1, the desired elevation to provide free charge on the charge terminal T1 for launching a guided surface wave should be at least 4-5 times the effective spherical diameter above the lossy conductive medium (e.g., the Earth). The compensation terminal T2 can be used to adjust the total effective height (hTE) of the guided surface waveguide probe 200d to excite an electric field having a guided surface wave tilt at Rx. The compensation terminal T2 can be positioned below the charge terminal T1 at hd=hT−hp, where hT is the total physical height of the charge terminal T1. With the position of the compensation terminal T2 fixed and the phase delay ΦU applied to the upper charge terminal T1, the phase delay ΦL applied to the lower compensation terminal T2 can be determined using the relationships of Equation (86), such that:




  
    
      
        
          
            Φ
            U
          
          ⁡
          
            (
            
              h
              d
            
            )
          
        
        =
        
          
            -
            
              β
              ⁡
              
                (
                
                  
                    h
                    p
                  
                  +
                  y
                
                )
              
            
          
          -
          
            j
            ⁢
            
              
            
            ⁢
            
              
                ln
                ⁡
                
                  (
                  
                    
                      
                        
                          R
                          x
                        
                        ×
                        w
                      
                      -
                      
                        
                          (
                          
                            
                              h
                              p
                            
                            +
                            z
                          
                          )
                        
                        ⁢
                        
                          e
                          
                            j
                            ⁢
                            
                              
                            
                            ⁢
                            
                              (
                              
                                
                                  β
                                  ⁢
                                  
                                    
                                  
                                  ⁢
                                  
                                    h
                                    p
                                  
                                
                                +
                                
                                  β
                                  ⁢
                                  
                                    
                                  
                                  ⁢
                                  z
                                
                                +
                                
                                  Φ
                                  L
                                
                              
                              )
                            
                          
                        
                      
                    
                    
                      (
                      
                        
                          h
                          d
                        
                        +
                        y
                      
                      )
                    
                  
                  )
                
              
              .
            
          
        
      
    
    
      
        (
        89
        )
      
    
  




In alternative embodiments, the compensation terminal T2 can be positioned at a height hd where Im{ΦL}=0. This is graphically illustrated in FIG. 15A, which shows plots 172 and 175 of the imaginary and real parts of ΦU, respectively. The compensation terminal T2 is positioned at a height hd where Im{ΦU}=0, as graphically illustrated in plot 172. At this fixed height, the coil phase ΦU can be determined from Re{ΦU}, as graphically illustrated in plot 175.

With the AC source 212 coupled to the coil 215 (e.g., at the 50Ω point to maximize coupling), the position of tap 233 may be adjusted for parallel resonance of the compensation terminal T2 with at least a portion of the coil at the frequency of operation. FIG. 15B shows a schematic diagram of the general electrical hookup of FIG. 14 in which V1 is the voltage applied to the lower portion of the coil 215 from the AC source 212 through tap 227, V2 is the voltage at tap 224 that is supplied to the upper charge terminal T1, and V3 is the voltage applied to the lower compensation terminal T2 through tap 233. The resistances Rp and Rd represent the ground return resistances of the charge terminal T1 and compensation terminal T2, respectively. The charge and compensation terminals T1 and T2 may be configured as spheres, cylinders, toroids, rings, hoods, or any other combination of capacitive structures. The size of the charge and compensation terminals T1 and T2 can be chosen to provide a sufficiently large surface for the charges Q1 and Q2 imposed on the terminals. In general, it is desirable to make the charge terminal T1 as large as practical. The size of the charge terminal T1 should be large enough to avoid ionization of the surrounding air, which can result in electrical discharge or sparking around the charge terminal. The self-capacitance Cp and Cd of the charge and compensation terminals T1 and T2 respectively, can be determined using, for example, equation (24).
As can be seen in FIG. 15B, a resonant circuit is formed by at least a portion of the inductance of the coil 215, the self-capacitance Cd of the compensation terminal T2, and the ground return resistance Rd associated with the compensation terminal T2. The parallel resonance can be established by adjusting the voltage V3 applied to the compensation terminal T2 (e.g., by adjusting a tap 233 position on the coil 215) or by adjusting the height and/or size of the compensation terminal T2 to adjust Cd. The position of the coil tap 233 can be adjusted for parallel resonance, which will result in the ground current through the ground stake 218 and through the ammeter 236 reaching a maximum point. After parallel resonance of the compensation terminal T2 has been established, the position of the tap 227 for the AC source 212 can be adjusted to the 50Ω point on the coil 215.
Voltage V2 from the coil 215 can be applied to the charge terminal T1, and the position of tap 224 can be adjusted such that the phase (0) of the total effective height (hTE) approximately equals the angle of the guided surface wave tilt (WRx) at the Hankel crossover distance (Rx). The position of the coil tap 224 can be adjusted until this operating point is reached, which results in the ground current through the ammeter 236 increasing to a maximum. At this point, the resultant fields excited by the guided surface waveguide probe 200d are substantially mode-matched to a guided surface waveguide mode on the surface of the lossy conducting medium 203, resulting in the launching of a guided surface wave along the surface of the lossy conducting medium 203. This can be verified by measuring field strength along a radial extending from the guided surface waveguide probe 200.
Resonance of the circuit including the compensation terminal T2 may change with the attachment of the charge terminal T1 and/or with adjustment of the voltage applied to the charge terminal T1 through tap 224. While adjusting the compensation terminal circuit for resonance aids the subsequent adjustment of the charge terminal connection, it is not necessary to establish the guided surface wave tilt (WRx) at the Hankel crossover distance (Rx). The system may be further adjusted to improve coupling by iteratively adjusting the position of the tap 227 for the AC source 212 to be at the 50Ω point on the coil 215 and adjusting the position of tap 233 to maximize the ground current through the ammeter 236. Resonance of the circuit including the compensation terminal T2 may drift as the positions of taps 227 and 233 are adjusted, or when other components are attached to the coil 215.
In other implementations, the voltage V2 from the coil 215 can be applied to the charge terminal T1, and the position of tap 233 can be adjusted such that the phase (0) of the total effective height (hTE) approximately equals the angle (Ψ) of the guided surface wave tilt at Rx. The position of the coil tap 224 can be adjusted until the operating point is reached, resulting in the ground current through the ammeter 236 substantially reaching a maximum. The resultant fields are substantially mode-matched to a guided surface waveguide mode on the surface of the lossy conducting medium 203, and a guided surface wave is launched along the surface of the lossy conducting medium 203. This can be verified by measuring field strength along a radial extending from the guided surface waveguide probe 200. The system may be further adjusted to improve coupling by iteratively adjusting the position of the tap 227 for the AC source 212 to be at the 50Ω point on the coil 215 and adjusting the position of tap 224 and/or 233 to maximize the ground current through the ammeter 236.
Referring back to FIG. 12, operation of a guided surface waveguide probe 200 may be controlled to adjust for variations in operational conditions associated with the guided surface waveguide probe 200. For example, a probe control system 230 can be used to control the coupling circuit 209 and/or positioning of the charge terminal T1 and/or compensation terminal T2 to control the operation of the guided surface waveguide probe 200. Operational conditions can include, but are not limited to, variations in the characteristics of the lossy conducting medium 203 (e.g., conductivity σ and relative permittivity ∈r), variations in field strength and/or variations in loading of the guided surface waveguide probe 200. As can be seen from Equations (41)-(44), the index of refraction (n), the complex Brewster angle (θi,B and ψi,B), the wave tilt (|W|ejΨ) and the complex effective height (heff=hpejΦ) can be affected by changes in soil conductivity and permittivity resulting from, e.g., weather conditions.
Equipment such as, e.g., conductivity measurement probes, permittivity sensors, ground parameter meters, field meters, current monitors and/or load receivers can be used to monitor for changes in the operational conditions and provide information about current operational conditions to the probe control system 230. The probe control system 230 can then make one or more adjustments to the guided surface waveguide probe 200 to maintain specified operational conditions for the guided surface waveguide probe 200. For instance, as the moisture and temperature vary, the conductivity of the soil will also vary. Conductivity measurement probes and/or permittivity sensors may be located at multiple locations around the guided surface waveguide probe 200. Generally, it would be desirable to monitor the conductivity and/or permittivity at or about the Hankel crossover distance R, for the operational frequency. Conductivity measurement probes and/or permittivity sensors may be located at multiple locations (e.g., in each quadrant) around the guided surface waveguide probe 200.
With reference then to FIG. 16, shown is an example of a guided surface waveguide probe 200e that includes a charge terminal T1 and a charge terminal T2 that are arranged along a vertical axis z. The guided surface waveguide probe 200e is disposed above a lossy conducting medium 203, which makes up Region 1. In addition, a second medium 206 shares a boundary interface with the lossy conducting medium 203 and makes up Region 2. The charge terminals T1 and T2 are positioned over the lossy conducting medium 203. The charge terminal T1 is positioned at height H1, and the charge terminal T2 is positioned directly below T1 along the vertical axis z at height H2, where H2 is less than H1. The height h of the transmission structure presented by the guided surface waveguide probe 200e is h=H1−H2. The guided surface waveguide probe 200e includes a probe coupling circuit 209 that couples an excitation source 212 to the charge terminals T1 and T2.
The charge terminals T1 and/or T2 include a conductive mass that can hold an electrical charge, which may be sized to hold as much charge as practically possible. The charge terminal T1 has a self-capacitance C1, and the charge terminal T2 has a self-capacitance C2, which can be determined using, for example, equation (24). By virtue of the placement of the charge terminal T1 directly above the charge terminal T2, a mutual capacitance CM is created between the charge terminals T1 and T2. Note that the charge terminals T1 and T2 need not be identical, but each can have a separate size and shape, and can include different conducting materials. Ultimately, the field strength of a guided surface wave launched by a guided surface waveguide probe 200e is directly proportional to the quantity of charge on the terminal T1. The charge Q1 is, in turn, proportional to the self-capacitance C1 associated with the charge terminal T1 since Q1=C1V, where V is the voltage imposed on the charge terminal T1.
When properly adjusted to operate at a predefined operating frequency, the guided surface waveguide probe 200e generates a guided surface wave along the surface of the lossy conducting medium 203. The excitation source 212 can generate electrical energy at the predefined frequency that is applied to the guided surface waveguide probe 200e to excite the structure. When the electromagnetic fields generated by the guided surface waveguide probe 200e are substantially mode-matched with the lossy conducting medium 203, the electromagnetic fields substantially synthesize a wave front incident at a complex Brewster angle that results in little or no reflection. Thus, the surface waveguide probe 200e does not produce a radiated wave, but launches a guided surface traveling wave along the surface of a lossy conducting medium 203. The energy from the excitation source 212 can be transmitted as Zenneck surface currents to one or more receivers that are located within an effective transmission range of the guided surface waveguide probe 200e. 
One can determine asymptotes of the radial Zenneck surface current Jρ(ρ) on the surface of the lossy conducting medium 203 to be J1(ρ) close-in and J2(ρ) far-out, where




  
    
      
        
          Close
          ⁢
          
            -
          
          ⁢
          in
          ⁢
          
            
          
          ⁢
          
            (
            
              ρ
              <
              
                λ
                /
                8
              
            
            )
          
          ⁢
          
            :
          
        
        ⁢
        
          

        
        ⁢
        
          
        
        ⁢
        
          
            
              
                
                  J
                  ρ
                
                ⁡
                
                  (
                  ρ
                  )
                
              
              ∼
              
                J
                1
              
            
            =
            
              
                
                  
                    I
                    1
                  
                  +
                  
                    I
                    2
                  
                
                
                  2
                  ⁢
                  πρ
                
              
              +
              
                
                  
                    
                      E
                      ρ
                      
                        Q
                        ⁢
                        
                          
                        
                        ⁢
                        S
                      
                    
                    ⁡
                    
                      (
                      
                        Q
                        1
                      
                      )
                    
                  
                  +
                  
                    
                      E
                      ρ
                      
                        Q
                        ⁢
                        
                          
                        
                        ⁢
                        S
                      
                    
                    ⁡
                    
                      (
                      
                        Q
                        2
                      
                      )
                    
                  
                
                
                  Z
                  ρ
                
              
            
          
          ,
          and
        
      
    
    
      
        (
        90
        )
      
    
  
  
    
      
        
          Far
          ⁢
          
            -
          
          ⁢
          out
          ⁢
          
            
          
          ⁢
          
            (
            
              ρ
              >>
              λ/8
            
            )
          
          ⁢
          
            :
          
        
        ⁢
        
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              
                J
                ρ
              
              ⁡
              
                (
                ρ
                )
              
            
            ⁢
            
              ∼
            
            ⁢
            
              J
              2
            
          
          =
          
            
              
                j
                ⁢
                
                  
                
                ⁢
                γω
                ⁢
                
                  
                
                ⁢
                
                  Q
                  1
                
              
              4
            
            ×
            
              
                
                  2
                  ⁢
                  γ
                
                π
              
            
            ×
            
              
                
                  e
                  
                    
                      -
                      
                        (
                        
                          α
                          +
                          
                            j
                            ⁢
                            
                              
                            
                            ⁢
                            β
                          
                        
                        )
                      
                    
                    ⁢
                    ρ
                  
                
                
                  ρ
                
              
              .
            
          
        
      
    
    
      
        (
        91
        )
      
    
  




where I1 is the conduction current feeding the charge Q1 on the first charge terminal T1, and I2 is the conduction current feeding the charge Q2 on the second charge terminal T2. The charge Q1 on the upper charge terminal T1 is determined by Q1=C1V1, where C1 is the isolated capacitance of the charge terminal T1. Note that there is a third component to J1 set forth above given by (EρQ1)/Zρ, which follows from the Leontovich boundary condition and is the radial current contribution in the lossy conducting medium 203 pumped by the quasi-static field of the elevated oscillating charge on the first charge terminal Q1. The quantity Zρ=jωμo/γe is the radial impedance of the lossy conducting medium, where γe(jωμ1σ1−ω2μ1∈1)1/2.

The asymptotes representing the radial current close-in and far-out as set forth by equations (90) and (91) are complex quantities. According to various embodiments, a physical surface current J(ρ), is synthesized to match as close as possible the current asymptotes in magnitude and phase. That is to say close-in, |J(ρ)| is to be tangent to |J1|, and far-out |J(β)| is to be tangent to |J2|. Also, according to the various embodiments, the phase of J(ρ) should transition from the phase of J1 close-in to the phase of J2 far-out.
In order to match the guided surface wave mode at the site of transmission to launch a guided surface wave, the phase of the surface current |J2| far-out should differ from the phase of the surface current |J1| close-in by the propagation phase corresponding to e−jβ(ρ2−ρ1) plus a constant of approximately 45 degrees or 225 degrees. This is because there are two roots for √{square root over (γ)}, one near π/4 and one near 5π/4. The properly adjusted synthetic radial surface current is




  
    
      
        
          
            J
            ρ
          
          ⁡
          
            (
            
              ρ
              ,
              φ
              ,
              0
            
            )
          
        
        =
        
          
            
              
                I
                o
              
              ⁢
              γ
            
            4
          
          ⁢
          
            
              
                H
                1
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γρ
                
                )
              
            
            .
          
        
      
    
    
      
        (
        92
        )
      
    
  




Note that this is consistent with equation (17). By Maxwell's equations, such a J(ρ) surface current automatically creates fields that conform to





  
    
      
        
          
            H
            φ
          
          =
          
            
              
                
                  -
                  γ
                
                ⁢
                
                  
                
                ⁢
                
                  I
                  o
                
              
              4
            
            ⁢
            
              e
              
                
                  -
                  
                    u
                    2
                  
                
                ⁢
                z
              
            
            ⁢
            
              
                H
                1
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                –jγρ
                )
              
            
          
        
        ,
      
    
    
      
        (
        93
        )
      
    
  
  
    
      
        
          
            E
            ρ
          
          =
          
            
              
                
                  -
                  γ
                
                ⁢
                
                  
                
                ⁢
                
                  I
                  o
                
              
              4
            
            ⁢
            
              (
              
                
                  u
                  2
                
                
                  j
                  ⁢
                  
                    
                  
                  ⁢
                  
                    ωɛ
                    o
                  
                
              
              )
            
            ⁢
            
              e
              
                
                  -
                  
                    u
                    2
                  
                
                ⁢
                z
              
            
            ⁢
            
              
                H
                1
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γρ
                
                )
              
            
          
        
        ,
        
          

        
        ⁢
        and
      
    
    
      
        (
        94
        )
      
    
  
  
    
      
        
          E
          z
        
        =
        
          
            
              
                -
                γ
              
              ⁢
              
                
              
              ⁢
              
                I
                o
              
            
            4
          
          ⁢
          
            (
            
              
                -
                γ
              
              
                ωɛ
                o
              
            
            )
          
          ⁢
          
            e
            
              
                -
                
                  u
                  2
                
              
              ⁢
              z
            
          
          ⁢
          
            
              
                H
                0
                
                  (
                  2
                  )
                
              
              ⁡
              
                (
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  γρ
                
                )
              
            
            .
          
        
      
    
    
      
        (
        95
        )
      
    
  




Thus, the difference in phase between the surface current |J2| far-out and the surface current |J1| close-in for the guided surface wave mode that is to be matched is due to the characteristics of the Hankel functions in equations (93)-(95), which are consistent with equations (1)-(3). It is of significance to recognize that the fields expressed by equations (1)-(6) and (17) and equations (92)-(95) have the nature of a transmission line mode bound to a lossy interface, not radiation fields that are associated with groundwave propagation.

In order to obtain the appropriate voltage magnitudes and phases for a given design of a guided surface waveguide probe 200e at a given location, an iterative approach may be used. Specifically, analysis may be performed of a given excitation and configuration of a guided surface waveguide probe 200e taking into account the feed currents to the terminals T1 and T2, the charges on the charge terminals T1 and T2, and their images in the lossy conducting medium 203 in order to determine the radial surface current density generated. This process may be performed iteratively until an optimal configuration and excitation for a given guided surface waveguide probe 200e is determined based on desired parameters. To aid in determining whether a given guided surface waveguide probe 200e is operating at an optimal level, a guided field strength curve 103 (FIG. 1) may be generated using equations (1)-(12) based on values for the conductivity of Region 1 (σ1) and the permittivity of Region 1 (∈1) at the location of the guided surface waveguide probe 200e. Such a guided field strength curve 103 can provide a benchmark for operation such that measured field strengths can be compared with the magnitudes indicated by the guided field strength curve 103 to determine if optimal transmission has been achieved.
In order to arrive at an optimized condition, various parameters associated with the guided surface waveguide probe 200e may be adjusted. One parameter that may be varied to adjust the guided surface waveguide probe 200e is the height of one or both of the charge terminals T1 and/or T2 relative to the surface of the lossy conducting medium 203. In addition, the distance or spacing between the charge terminals T1 and T2 may also be adjusted. In doing so, one may minimize or otherwise alter the mutual capacitance CM or any bound capacitances between the charge terminals T1 and T2 and the lossy conducting medium 203 as can be appreciated. The size of the respective charge terminals T1 and/or T2 can also be adjusted. By changing the size of the charge terminals T1 and/or T2, one will alter the respective self-capacitances C1 and/or C2, and the mutual capacitance CM as can be appreciated.
Still further, another parameter that can be adjusted is the probe coupling circuit 209 associated with the guided surface waveguide probe 200e. This may be accomplished by adjusting the size of the inductive and/or capacitive reactances that make up the probe coupling circuit 209. For example, where such inductive reactances comprise coils, the number of turns on such coils may be adjusted. Ultimately, the adjustments to the probe coupling circuit 209 can be made to alter the electrical length of the probe coupling circuit 209, thereby affecting the voltage magnitudes and phases on the charge terminals T1 and T2.
Note that the iterations of transmission performed by making the various adjustments may be implemented by using computer models or by adjusting physical structures as can be appreciated. By making the above adjustments, one can create corresponding “close-in” surface current J1 and “far-out” surface current J2 that approximate the same currents J(ρ) of the guided surface wave mode specified in Equations (90) and (91) set forth above. In doing so, the resulting electromagnetic fields would be substantially or approximately mode-matched to a guided surface wave mode on the surface of the lossy conducting medium 203.
While not shown in the example of FIG. 16, operation of the guided surface waveguide probe 200e may be controlled to adjust for variations in operational conditions associated with the guided surface waveguide probe 200. For example, a probe control system 230 shown in FIG. 12 can be used to control the coupling circuit 209 and/or positioning and/or size of the charge terminals T1 and/or T2 to control the operation of the guided surface waveguide probe 200e. Operational conditions can include, but are not limited to, variations in the characteristics of the lossy conducting medium 203 (e.g., conductivity σ and relative permittivity ∈r), variations in field strength and/or variations in loading of the guided surface waveguide probe 200e. 
Referring now to FIG. 17, shown is an example of the guided surface waveguide probe 200e of FIG. 16, denoted herein as guided surface waveguide probe 200f. The guided surface waveguide probe 200f includes the charge terminals T1 and T2 that are positioned along a vertical axis z that is substantially normal to the plane presented by the lossy conducting medium 203 (e.g., the Earth). The second medium 206 is above the lossy conducting medium 203. The charge terminal T1 has a self-capacitance C1, and the charge terminal T2 has a self-capacitance C2. During operation, charges Q1 and Q2 are imposed on the charge terminals T1 and T2, respectively, depending on the voltages applied to the charge terminals T1 and T2 at any given instant. A mutual capacitance CM may exist between the charge terminals T1 and T2 depending on the distance there between. In addition, bound capacitances may exist between the respective charge terminals T1 and T2 and the lossy conducting medium 203 depending on the heights of the respective charge terminals T1 and T2 with respect to the lossy conducting medium 203.
The guided surface waveguide probe 200f includes a probe coupling circuit 209 that comprises an inductive impedance comprising a coil L1a having a pair of leads that are coupled to respective ones of the charge terminals T1 and T2. In one embodiment, the coil L1a is specified to have an electrical length that is one-half (1⁄2) of the wavelength at the operating frequency of the guided surface waveguide probe 200f. 
While the electrical length of the coil L1a is specified as approximately one-half (1⁄2) the wavelength at the operating frequency, it is understood that the coil L1a may be specified with an electrical length at other values. According to one embodiment, the fact that the coil L1a has an electrical length of approximately one-half the wavelength at the operating frequency provides for an advantage in that a maximum voltage differential is created on the charge terminals T1 and T2. Nonetheless, the length or diameter of the coil L1a may be increased or decreased when adjusting the guided surface waveguide probe 200f to obtain optimal excitation of a guided surface wave mode. Adjustment of the coil length may be provided by taps located at one or both ends of the coil. In other embodiments, it may be the case that the inductive impedance is specified to have an electrical length that is significantly less than or greater than 1⁄2 the wavelength at the operating frequency of the guided surface waveguide probe 200f. 
The excitation source 212 can be coupled to the probe coupling circuit 209 by way of magnetic coupling. Specifically, the excitation source 212 is coupled to a coil LP that is inductively coupled to the coil L1a. This may be done by link coupling, a tapped coil, a variable reactance, or other coupling approach as can be appreciated. To this end, the coil LP acts as a primary, and the coil L1a acts as a secondary as can be appreciated.
In order to adjust the guided surface waveguide probe 200f for the transmission of a desired guided surface wave, the heights of the respective charge terminals T1 and T2 may be altered with respect to the lossy conducting medium 203 and with respect to each other. Also, the sizes of the charge terminals T1 and T2 may be altered. In addition, the size of the coil L1a may be altered by adding or eliminating turns or by changing some other dimension of the coil L1a. The coil L1a can also include one or more taps for adjusting the electrical length as shown in FIG. 17. The position of a tap connected to either charge terminal T1 or T2 can also be adjusted.
Referring next to FIGS. 18A, 18B, 18C and 19, shown are examples of generalized receive circuits for using the surface-guided waves in wireless power delivery systems. FIGS. 18A and 18B-18C include a linear probe 303 and a tuned resonator 306, respectively. FIG. 19 is a magnetic coil 309 according to various embodiments of the present disclosure. According to various embodiments, each one of the linear probe 303, the tuned resonator 306, and the magnetic coil 309 may be employed to receive power transmitted in the form of a guided surface wave on the surface of a lossy conducting medium 203 according to various embodiments. As mentioned above, in one embodiment the lossy conducting medium 203 comprises a terrestrial medium (or Earth).
With specific reference to FIG. 18A, the open-circuit terminal voltage at the output terminals 312 of the linear probe 303 depends upon the effective height of the linear probe 303. To this end, the terminal point voltage may be calculated as

VT=∫0heEinc·dl,  (96)

where ∈inc is the strength of the incident electric field induced on the linear probe 303 in Volts per meter, dl is an element of integration along the direction of the linear probe 303, and he is the effective height of the linear probe 303. An electrical load 315 is coupled to the output terminals 312 through an impedance matching network 318.

When the linear probe 303 is subjected to a guided surface wave as described above, a voltage is developed across the output terminals 312 that may be applied to the electrical load 315 through a conjugate impedance matching network 318 as the case may be. In order to facilitate the flow of power to the electrical load 315, the electrical load 315 should be substantially impedance matched to the linear probe 303 as will be described below.
Referring to FIG. 18B, a ground current excited coil 306a possessing a phase shift equal to the wave tilt of the guided surface wave includes a charge terminal TR that is elevated (or suspended) above the lossy conducting medium 203. The charge terminal TR has a self-capacitance CR. In addition, there may also be a bound capacitance (not shown) between the charge terminal TR and the lossy conducting medium 203 depending on the height of the charge terminal TR above the lossy conducting medium 203. The bound capacitance should preferably be minimized as much as is practicable, although this may not be entirely necessary in every instance.
The tuned resonator 306a also includes a receiver network comprising a coil LR having a phase shift Φ. One end of the coil LR is coupled to the charge terminal TR, and the other end of the coil LR is coupled to the lossy conducting medium 203. The receiver network can include a vertical supply line conductor that couples the coil LR to the charge terminal TR. To this end, the coil LR (which may also be referred to as tuned resonator LR-CR) comprises a series-adjusted resonator as the charge terminal CR and the coil LR are situated in series. The phase delay of the coil LR can be adjusted by changing the size and/or height of the charge terminal TR, and/or adjusting the size of the coil LR so that the phase Φ of the structure is made substantially equal to the angle of the wave tilt Ψ. The phase delay of the vertical supply line can also be adjusted by, e.g., changing length of the conductor.
For example, the reactance presented by the self-capacitance CR is calculated as 1/jωCR. Note that the total capacitance of the structure 306a may also include capacitance between the charge terminal TR and the lossy conducting medium 203, where the total capacitance of the structure 306a may be calculated from both the self-capacitance CR and any bound capacitance as can be appreciated. According to one embodiment, the charge terminal TR may be raised to a height so as to substantially reduce or eliminate any bound capacitance. The existence of a bound capacitance may be determined from capacitance measurements between the charge terminal TR and the lossy conducting medium 203 as previously discussed.
The inductive reactance presented by a discrete-element coil LR may be calculated as jωL, where L is the lumped-element inductance of the coil LR. If the coil LR is a distributed element, its equivalent terminal-point inductive reactance may be determined by conventional approaches. To tune the structure 306a, one would make adjustments so that the phase delay is equal to the wave tilt for the purpose of mode-matching to the surface waveguide at the frequency of operation. Under this condition, the receiving structure may be considered to be “mode-matched” with the surface waveguide. A transformer link around the structure and/or an impedance matching network 324 may be inserted between the probe and the electrical load 327 in order to couple power to the load. Inserting the impedance matching network 324 between the probe terminals 321 and the electrical load 327 can effect a conjugate-match condition for maximum power transfer to the electrical load 327.
When placed in the presence of surface currents at the operating frequencies power will be delivered from the surface guided wave to the electrical load 327. To this end, an electrical load 327 may be coupled to the structure 306a by way of magnetic coupling, capacitive coupling, or conductive (direct tap) coupling. The elements of the coupling network may be lumped components or distributed elements as can be appreciated.
In the embodiment shown in FIG. 18B, magnetic coupling is employed where a coil LS is positioned as a secondary relative to the coil LR that acts as a transformer primary. The coil LS may be link-coupled to the coil LR by geometrically winding it around the same core structure and adjusting the coupled magnetic flux as can be appreciated. In addition, while the receiving structure 306a comprises a series-tuned resonator, a parallel-tuned resonator or even a distributed-element resonator of the appropriate phase delay may also be used.
While a receiving structure immersed in an electromagnetic field may couple energy from the field, it can be appreciated that polarization-matched structures work best by maximizing the coupling, and conventional rules for probe-coupling to waveguide modes should be observed. For example, a TE20 (transverse electric mode) waveguide probe may be optimal for extracting energy from a conventional waveguide excited in the TE20 mode. Similarly, in these cases, a mode-matched and phase-matched receiving structure can be optimized for coupling power from a surface-guided wave. The guided surface wave excited by a guided surface waveguide probe 200 on the surface of the lossy conducting medium 203 can be considered a waveguide mode of an open waveguide. Excluding waveguide losses, the source energy can be completely recovered. Useful receiving structures may be E-field coupled, H-field coupled, or surface-current excited.
The receiving structure can be adjusted to increase or maximize coupling with the guided surface wave based upon the local characteristics of the lossy conducting medium 203 in the vicinity of the receiving structure. To accomplish this, the phase delay (Φ) of the receiving structure can be adjusted to match the angle (Ψ) of the wave tilt of the surface traveling wave at the receiving structure. If configured appropriately, the receiving structure may then be tuned for resonance with respect to the perfectly conducting image ground plane at complex depth Z=−d/2.
For example, consider a receiving structure comprising the tuned resonator 306a of FIG. 18B, including a coil LR and a vertical supply line connected between the coil LR and a charge terminal TR. With the charge terminal TR positioned at a defined height above the lossy conducting medium 203, the total phase shift Φ of the coil LR and vertical supply line can be matched with the angle (Ψ) of the wave tilt at the location of the tuned resonator 306a. From Equation (22), it can be seen that the wave tilt asymptotically passes to




  
    
      
        
          W
          =
          
            
              
                
                W
                
              
              ⁢
              
                e
                
                  j
                  ⁢
                  
                    
                  
                  ⁢
                  Ψ
                
              
            
            =
            
              
                
                  E
                  ρ
                
                
                  E
                  z
                
              
              ⁢
              
                ⟶
                
                  ρ
                  →
                  ∞
                
              
              ⁢
              
                1
                
                  
                    
                      ɛ
                      r
                    
                    -
                    
                      j
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                          σ
                          1
                        
                        
                          ωɛ
                          o
                        
                      
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        97
        )
      
    
  




where ∈r comprises the relative permittivity and σ1 is the conductivity of the lossy conducting medium 203 at the location of the receiving structure, ∈o is the permittivity of free space, and ψ=2πf, where f is the frequency of excitation. Thus, the wave tilt angle (Ψ) can be determined from Equation (97).

The total phase shift (Φ=θc+θy) of the tuned resonator 306a includes both the phase delay (θc) through the coil LR and the phase delay of the vertical supply line (θy). The spatial phase delay along the conductor length lw of the vertical supply line can be given by θy=βwlw, where βw is the propagation phase constant for the vertical supply line conductor. The phase delay due to the coil (or helical delay line) is θc=βplC, with a physical length of lC and a propagation factor of




  
    
      
        
          
            β
            p
          
          =
          
            
              
                2
                ⁢
                π
              
              
                λ
                p
              
            
            =
            
              
                2
                ⁢
                π
              
              
                
                  V
                  f
                
                ⁢
                
                  λ
                  0
                
              
            
          
        
        ,
      
    
    
      
        (
        98
        )
      
    
  




where Vf is the velocity factor on the structure, λ0 is the wavelength at the supplied frequency, and λp is the propagation wavelength resulting from the velocity factor Vf. One or both of the phase delays (θc+θy) can be adjusted to match the phase shift Φ to the angle (Ψ) of the wave tilt. For example, a tap position may be adjusted on the coil LR of FIG. 18B to adjust the coil phase delay (θc) to match the total phase shift to the wave tilt angle (Φ=Ψ). For example, a portion of the coil can be bypassed by the tap connection as illustrated in FIG. 18B. The vertical supply line conductor can also be connected to the coil LR via a tap, whose position on the coil may be adjusted to match the total phase shift to the angle of the wave tilt.

Once the phase delay (Ψ) of the tuned resonator 306a has been adjusted, the impedance of the charge terminal TR can then be adjusted to tune to resonance with respect to the perfectly conducting image ground plane at complex depth z=−d/2. This can be accomplished by adjusting the capacitance of the charge terminal T1 without changing the traveling wave phase delays of the coil LR and vertical supply line. The adjustments are similar to those described with respect to FIGS. 9A and 9B.
The impedance seen “looking down” into the lossy conducting medium 203 to the complex image plane is given by:

Zin=Rin+jXin=Zo tan h(jβo(d/2)),  (99)

where βo=ω√{square root over (μo∈o)}. For vertically polarized sources over the Earth, the depth of the complex image plane can be given by:

d/2≈1/√{square root over (jωμ1σ1−ω2μ1∈1)},  (100)

where μ1 is the permeability of the lossy conducting medium 203 and ∈1=∈r∈o.

At the base of the tuned resonator 306a, the impedance seen “looking up” into the receiving structure is Z↑=Zbase as illustrated in FIG. 9A. With a terminal impedance of:




  
    
      
        
          
            Z
            R
          
          =
          
            1
            
              j
              ⁢
              
                
              
              ⁢
              ω
              ⁢
              
                
              
              ⁢
              
                C
                
                  R
                  ⁢
                  
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        101
        )
      
    
  




where CR is the self-capacitance of the charge terminal TR, the impedance seen “looking up” into the vertical supply line conductor of the tuned resonator 306a is given by:





  
    
      
        
          
            Z
            2
          
          =
          
            
              
                Z
                W
              
              ⁢
              
                
                  
                    Z
                    R
                  
                  +
                  
                    
                      Z
                      w
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          
                            jβ
                            w
                          
                          ⁢
                          
                            h
                            w
                          
                        
                        )
                      
                    
                  
                
                
                  
                    Z
                    w
                  
                  +
                  
                    
                      Z
                      R
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          
                            jβ
                            w
                          
                          ⁢
                          
                            h
                            w
                          
                        
                        )
                      
                    
                  
                
              
            
            =
            
              
                Z
                W
              
              ⁢
              
                
                  
                    Z
                    R
                  
                  +
                  
                    
                      Z
                      w
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          jθ
                          y
                        
                        )
                      
                    
                  
                
                
                  
                    Z
                    w
                  
                  +
                  
                    
                      Z
                      R
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          jθ
                          y
                        
                        )
                      
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        102
        )
      
    
  




and the impedance seen “looking up” into the coil LR of the tuned resonator 306a is given by:





  
    
      
        
          Z
          base
        
        =
        
          
            
              R
              base
            
            +
            
              j
              ⁢
              
                
              
              ⁢
              
                X
                base
              
            
          
          =
          
            
              
                Z
                R
              
              ⁢
              
                
                  
                    Z
                    2
                  
                  +
                  
                    
                      Z
                      R
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            p
                          
                          ⁢
                          H
                        
                        )
                      
                    
                  
                
                
                  
                    Z
                    R
                  
                  +
                  
                    
                      Z
                      2
                    
                    ⁢
                    
                      tanh
                      ⁡
                      
                        (
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            p
                          
                          ⁢
                          H
                        
                        )
                      
                    
                  
                
              
            
            =
            
              
                Z
                c
              
              ⁢
              
                
                  
                    
                      Z
                      2
                    
                    +
                    
                      
                        Z
                        R
                      
                      ⁢
                      
                        tanh
                        ⁡
                        
                          (
                          
                            j
                            ⁢
                            
                              
                            
                            ⁢
                            
                              θ
                              c
                            
                          
                          )
                        
                      
                    
                  
                  
                    
                      Z
                      R
                    
                    +
                    
                      
                        Z
                        2
                      
                      ⁢
                      
                        tanh
                        ⁡
                        
                          (
                          
                            j
                            ⁢
                            
                              
                            
                            ⁢
                            
                              θ
                              c
                            
                          
                          )
                        
                      
                    
                  
                
                .
              
            
          
        
      
    
    
      
        (
        103
        )
      
    
  




By matching the reactive component (Xin) seen “looking down” into the lossy conducting medium 203 with the reactive component (Xbase) seen “looking up” into the tuned resonator 306a, the coupling into the guided surface waveguide mode may be maximized.

Referring next to FIG. 18C, shown is an example of a tuned resonator 306b that does not include a charge terminal TR at the top of the receiving structure. In this embodiment, the tuned resonator 306b does not include a vertical supply line coupled between the coil LR and the charge terminal TR. Thus, the total phase shift (Φ) of the tuned resonator 306b includes only the phase delay (θc) through the coil LR. As with the tuned resonator 306a of FIG. 18B, the coil phase delay θc can be adjusted to match the angle (Ψ) of the wave tilt determined from Equation (97), which results in Φ=Ψ. While power extraction is possible with the receiving structure coupled into the surface waveguide mode, it is difficult to adjust the receiving structure to maximize coupling with the guided surface wave without the variable reactive load provided by the charge terminal TR.
Referring to FIG. 18D, shown is a flow chart 180 illustrating an example of adjusting a receiving structure to substantially mode-match to a guided surface waveguide mode on the surface of the lossy conducting medium 203. Beginning with 181, if the receiving structure includes a charge terminal TR (e.g., of the tuned resonator 306a of FIG. 18B), then the charge terminal TR is positioned at a defined height above a lossy conducting medium 203 at 184. As the surface guided wave has been established by a guided surface waveguide probe 200, the physical height (hp) of the charge terminal TR may be below that of the effective height. The physical height may be selected to reduce or minimize the bound charge on the charge terminal TR (e.g., four times the spherical diameter of the charge terminal). If the receiving structure does not include a charge terminal TR (e.g., of the tuned resonator 306b of FIG. 18C), then the flow proceeds to 187.
At 187, the electrical phase delay Φ of the receiving structure is matched to the complex wave tilt angle Ψ defined by the local characteristics of the lossy conducting medium 203. The phase delay (θc) of the helical coil and/or the phase delay (θy) of the vertical supply line can be adjusted to make Φ equal to the angle (Ψ) of the wave tilt (W). The angle (Ψ) of the wave tilt can be determined from Equation (86). The electrical phase Φ can then be matched to the angle of the wave tilt. For example, the electrical phase delay Φ=θc+θy can be adjusted by varying the geometrical parameters of the coil LR and/or the length (or height) of the vertical supply line conductor.
Next at 190, the load impedance of the charge terminal TR can be tuned to resonate the equivalent image plane model of the tuned resonator 306a. The depth (d/2) of the conducting image ground plane 139 (FIG. 9A) below the receiving structure can be determined using Equation (100) and the values of the lossy conducting medium 203 (e.g., the Earth) at the receiving structure, which can be locally measured. Using that complex depth, the phase shift (θd) between the image ground plane 139 and the physical boundary 136 (FIG. 9A) of the lossy conducting medium 203 can be determined using θd=βo d/2. The impedance (Z1) as seen “looking down” into the lossy conducting medium 203 can then be determined using Equation (99). This resonance relationship can be considered to maximize coupling with the guided surface waves.
Based upon the adjusted parameters of the coil LR and the length of the vertical supply line conductor, the velocity factor, phase delay, and impedance of the coil LR and vertical supply line can be determined. In addition, the self-capacitance (CR) of the charge terminal TR can be determined using, e.g., Equation (24). The propagation factor (βp) of the coil LR can be determined using Equation (98), and the propagation phase constant (βw) for the vertical supply line can be determined using Equation (49). Using the self-capacitance and the determined values of the coil LR and vertical supply line, the impedance (Zbase) of the tuned resonator 306a as seen “looking up” into the coil LR can be determined using Equations (101), (102), and (103).
The equivalent image plane model of FIG. 9A also applies to the tuned resonator 306a of FIG. 18B. The tuned resonator 306a can be tuned to resonance with respect to the complex image plane by adjusting the load impedance ZR of the charge terminal TR such that the reactance component Xbase of Zbase cancels out the reactance component of Xin of Zin, or Xbase+Xin=0. Thus, the impedance at the physical boundary 136 (FIG. 9A) “looking up” into the coil of the tuned resonator 306a is the conjugate of the impedance at the physical boundary 136 “looking down” into the lossy conducting medium 203. The load impedance ZR can be adjusted by varying the capacitance (CR) of the charge terminal TR without changing the electrical phase delay Φ=θc+θy seen by the charge terminal TR. An iterative approach may be taken to tune the load impedance ZR for resonance of the equivalent image plane model with respect to the conducting image ground plane 139. In this way, the coupling of the electric field to a guided surface waveguide mode along the surface of the lossy conducting medium 203 (e.g., Earth) can be improved and/or maximized.
Referring to FIG. 19, the magnetic coil 309 comprises a receive circuit that is coupled through an impedance matching network 333 to an electrical load 336. In order to facilitate reception and/or extraction of electrical power from a guided surface wave, the magnetic coil 309 may be positioned so that the magnetic flux of the guided surface wave, Hφ, passes through the magnetic coil 309, thereby inducing a current in the magnetic coil 309 and producing a terminal point voltage at its output terminals 330. The magnetic flux of the guided surface wave coupled to a single turn coil is expressed by

=∫∫ACSμrμo{right arrow over (H)}·{circumflex over (n)}dA  (104)

where  is the coupled magnetic flux, μr is the effective relative permeability of the core of the magnetic coil 309, μo is the permeability of free space, {right arrow over (H)} is the incident magnetic field strength vector, {circumflex over (n)} is a unit vector normal to the cross-sectional area of the turns, and ACS is the area enclosed by each loop. For an N-turn magnetic coil 309 oriented for maximum coupling to an incident magnetic field that is uniform over the cross-sectional area of the magnetic coil 309, the open-circuit induced voltage appearing at the output terminals 330 of the magnetic coil 309 is





  
    
      
        
          V
          =
          
            
              
                -
                N
              
              ⁢
              
                
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  F
                
                
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
              
            
            ≈
            
              
                -
                j
              
              ⁢
              
                
              
              ⁢
              
                ωμ
                r
              
              ⁢
              
                μ
                0
              
              ⁢
              
                NHA
                CS
              
            
          
        
        ,
      
    
    
      
        (
        105
        )
      
    
  




where the variables are defined above. The magnetic coil 309 may be tuned to the guided surface wave frequency either as a distributed resonator or with an external capacitor across its output terminals 330, as the case may be, and then impedance-matched to an external electrical load 336 through a conjugate impedance matching network 333.

Assuming that the resulting circuit presented by the magnetic coil 309 and the electrical load 336 are properly adjusted and conjugate impedance matched, via impedance matching network 333, then the current induced in the magnetic coil 309 may be employed to optimally power the electrical load 336. The receive circuit presented by the magnetic coil 309 provides an advantage in that it does not have to be physically connected to the ground.
With reference to FIGS. 18A, 18B, 18C and 19, the receive circuits presented by the linear probe 303, the mode-matched structure 306, and the magnetic coil 309 each facilitate receiving electrical power transmitted from any one of the embodiments of guided surface waveguide probes 200 described above. To this end, the energy received may be used to supply power to an electrical load 315/327/336 via a conjugate matching network as can be appreciated. This contrasts with the signals that may be received in a receiver that were transmitted in the form of a radiated electromagnetic field. Such signals have very low available power, and receivers of such signals do not load the transmitters.
It is also characteristic of the present guided surface waves generated using the guided surface waveguide probes 200 described above that the receive circuits presented by the linear probe 303, the mode-matched structure 306, and the magnetic coil 309 will load the excitation source 212 (e.g., FIGS. 3, 12 and 16) that is applied to the guided surface waveguide probe 200, thereby generating the guided surface wave to which such receive circuits are subjected. This reflects the fact that the guided surface wave generated by a given guided surface waveguide probe 200 described above comprises a transmission line mode. By way of contrast, a power source that drives a radiating antenna that generates a radiated electromagnetic wave is not loaded by the receivers, regardless of the number of receivers employed.
Thus, together one or more guided surface waveguide probes 200 and one or more receive circuits in the form of the linear probe 303, the tuned mode-matched structure 306, and/or the magnetic coil 309 can make up a wireless distribution system. Given that the distance of transmission of a guided surface wave using a guided surface waveguide probe 200 as set forth above depends upon the frequency, it is possible that wireless power distribution can be achieved across wide areas and even globally.
The conventional wireless-power transmission/distribution systems extensively investigated today include “energy harvesting” from radiation fields and also sensor coupling to inductive or reactive near-fields. In contrast, the present wireless-power system does not waste power in the form of radiation which, if not intercepted, is lost forever. Nor is the presently disclosed wireless-power system limited to extremely short ranges as with conventional mutual-reactance coupled near-field systems. The wireless-power system disclosed herein probe-couples to the novel surface-guided transmission line mode, which is equivalent to delivering power to a load by a wave-guide or a load directly wired to the distant power generator. Not counting the power required to maintain transmission field strength plus that dissipated in the surface waveguide, which at extremely low frequencies is insignificant relative to the transmission losses in conventional high-tension power lines at 60 Hz, all of the generator power goes only to the desired electrical load. When the electrical load demand is terminated, the source power generation is relatively idle.
Referring next to FIGS. 20A-20B, shown are examples of various schematic symbols that are used with reference to the discussion that follows. With specific reference to FIG. 20A, shown is a symbol that represents any one of the guided surface waveguide probes 200a, 200b, 200c, 200e, 200d, 200f, or any variations thereof. In the following drawings and discussion, a depiction of this symbol will be referred to as a guided surface waveguide probe P. For the sake of simplicity in the following discussion, any reference to the guided surface waveguide probe P is a reference to any one of the guided surface waveguide probes 200a, 200b, 200c, 200e, 200d, 200f, or variations thereof.
Similarly, with reference to FIG. 20B, shown is a symbol that represents a guided surface wave receive structure that may comprise any one of the linear probe 303 (FIG. 18A), the tuned resonator 306 (FIGS. 18B-18C), or the magnetic coil 309 (FIG. 19). In the following drawings and discussion, a depiction of this symbol will be referred to as a guided surface wave receive structure R. For the sake of simplicity in the following discussion, any reference to the guided surface wave receive structure R is a reference to any one of the linear probe 303, the tuned resonator 306, the magnetic coil 309, or variations thereof.
With reference to the remaining figures, examples of systems and methods for detecting the unauthorized consumption of electrical energy are described. As described above, a guided surface waveguide probe P can be used to transmit electrical energy via guided surface waves. As will be described in further detail below, maps indicating the flow of electrical energy in a geographic region associated with the guided surface waveguide probe P can be constructed. The map can be used to detect and locate guided surface wave receive structures R and their associated electrical loads that may be obtaining electrical energy without authorization from the entity that operates the guided surface waveguide probe P.
As a non-limiting example, metering systems can be distributed throughout the geographical region to measure the electric and magnetic fields present at various locations in the geographical region. Using the electric and magnetic field measurements, the energy flow vectors or Poynting vectors Ŝ that represent the directional energy flux density of the electromagnetic fields of the guided surfaces waves transmitted can be calculated. In addition, the energy flow vectors can be plotted on a map to represent the direction and magnitude of the flow of energy at the various locations in the geographical region. The energy flow map can then be compared to a map that represents the flow of energy that would be expected if there were no unauthorized guided surface waveguide receive structures R operating in the geographic region. If there is a discrepancy between the map representing the determined flow of energy and the map representing the expected flow of energy, it can be determined that an unauthorized guided surface wave receive structure R is operating in the geographic region. Additionally, the energy flow map can be used to locate the unauthorized guided surface wave receive structure R.
With reference to FIG. 21, shown is a networked environment 2100 according to various embodiments. The networked environment 2100 includes a computing environment 2103, one or more probe controllers 2104, one or more receive structure controllers 2105, and one or more metering systems 2106 that can be in data communication with each other via a network 2109. The network 2190 can include, for example, the Internet, one or more intranets, extranets, wide area networks (WANs), local area networks (LANs), other suitable networks, or any combination of two or more such networks. The networked environment 2103 can include wired networks, wireless networks, satellite networks, cable networks, Ethernet networks, and other types of networks.
The computing environment 2103 includes, for example, a server computer or any other system providing computing capability. Alternatively, the computing environment 2103 can employ a plurality of computing devices that can be arranged, for example, in one or more server banks or computer banks or other arrangements. Such computing devices can be located in a single installation or can be distributed among many different geographical locations. For example, some computing devices can be located in the metering systems 2106. In addition, the computing environment 2103 can include a plurality of computing devices that together may comprise a hosted computing resource, a grid computing resource and/or any other distributed computing arrangement. In some cases, the computing environment 2103 can correspond to an elastic computing resource where the allotted capacity of processing, network, storage, or other computing-related resources can vary over time.
Various applications and systems can be executed in the computing environment 2103. Also, various data can be stored in a data store 2113 that can be accessible to the computing environment 2103. The data store 2103 illustrated in FIG. 21 can represent one or more data stores 2113. The data stored in the data store 2113 can be associated with the operation of the various components described below.
The components executed in the computing environment 103 can include an energy usage monitor 2116 and other applications, services, processes, systems, engines, or functionality not discussed in detail herein. The energy usage monitor 2116 monitors the operation of the associated guided surface waveguide probes P and guided surface wave receive structures R. In addition, the energy usage monitor 2116 detects whether another guided surface wave receive structure R is obtaining electrical energy via guided surface waves without authorization from the entity that generates the guided surface waves. In other words, the energy usage monitor 2116 detects whether an unauthorized guided surface wave receive structure R is obtaining electrical energy from a guided surface waveguide probe P. In addition, the energy usage monitor 2116 determines the location of the unauthorized guided surface wave receive structure R, as will be described in further detail below.
The data stored in the data store 2113 includes expected energy flow maps 2119, actual energy flow maps 2123, and potentially other data. An expected energy flow map 2119 represents the flow of electrical energy that would be expected to occur given a particular configuration of one or more guided surface waveguide probes P and guided surface wave receive structures R. By contrast, an actual energy flow map 2123 represents the flow of electrical energy that has been determined to have actually occurred at particular locations within a geographical region. In some examples, an expected energy flow map 2119 and an actual energy flow map 2123 are embodied as graphical maps that graphically represent energy flow vectors at various geographical locations. As will be described in further detail below, the energy usage monitor 2116 generates the expected energy flow maps 2119 and the actual energy flow maps 2123.
The probe controller 2104 is representative of a plurality of probe controllers 2104 that can be in communication with other components through the network 2109. The probe controller 2104 can collect information regarding the operation of a corresponding guided surface waveguide probe P. For example, the probe controller 2104 can obtain information regarding an instantaneous current, voltage, power and/or other information associated with the guided surface waveguide probe P. The probe controller 2104 can also include a positioning system, such as a global positioning system (GPS), that can generate data representing the location of the corresponding guided surface waveguide probe P. In addition, the probe controller 2104 can include a transceiver that can communicate with the energy usage monitor 2116.
The receive structure controller 2105 is representative of a plurality of receive structure controllers 2105 that can be in communication with other components through the network 2109. The receive structure controller 2105 can collect information regarding the operation of a corresponding guided surface wave receive structure R. For example, the receive structure controller 2105 can obtain information regarding an instantaneous current, voltage, power and/or other information associated with the guided surface wave receive structure R. The receive structure controller 2105 can also include a positioning system, such as a global positioning system (GPS), that can generate data representing the location of the corresponding guided surface wave receive structure R. In addition, the receive structure controller 2105 can include a transceiver that can communicate with the energy usage monitor 2116.
The metering system 2106 is representative of a plurality of metering systems 2106 that can be in communication with other components through the network 2109. The metering system 2106 measures electric fields and magnetic fields at the particular location where the metering system 2106 is located in a geographical region. To this end, the metering system 2106 can include electric field meters 2106 and magnetic field meters 2109. The electric field meters 2126 and the magnetic field meters 2129 can sample and measure the magnitudes of the orthogonal components of the electric and magnetic fields, respectively. In some examples, the metering system 2106 includes three electric field meters 2106, and each electric field meter 2126 measures the magnitude of the electric field in one of the {circumflex over (x)}, ŷ, and {circumflex over (z)} directions. Similarly, the metering system 2126 in some examples includes three magnetic field meters 2129, and each magnetic field meter 2109 measures the magnitude of the magnetic field in one of the {circumflex over (x)}, ŷ, and {circumflex over (z)} directions. After measuring the orthogonal components of the electric and magnetic fields, the metering system 2106 can transmit this data to the energy usage monitor 2116. The measured electric fields can be represented as an electric field vector Ê, and the measured magnetic fields can be represented as a magnetic field vector Ĥ. Alternatively, the metering system 2106 may calculate a Poynting vector Ŝ corresponding to the measured electric and magnetic fields and send data indicating the magnitude and direction of the poynting vector Ŝ to the energy usage monitor 2116.
Several metering systems 2106 are distributed throughout the geographical region in which the guided surface waveguide probes P and/or guided surface wave receive structures R operate. In some examples, the metering systems 2106 are arranged in an array, thereby forming a grid in the geographical region in which the guided surface waveguide probes P and/or guided surface wave receive structures R operate. The metering systems 2106 can sample and measure the magnitudes of the orthogonal components of the electric and magnetic fields at their respective locations and transmit this data to the computing environment 2103 for further processing, as will be described in further detail below.
Next, a general description of the operation of the networked environment 2100 is provided. To begin, the energy usage monitor 2116 can generate one or more expected energy flow maps 2119. As discussed above, an expected energy flow map 2119 can represent the energy flow that would be expected to occur given a particular configuration of one or more guided surface waveguide probes P and guided surface wave receive structures R. For example, one expected energy flow map 2119 can represent the expected flow of energy if all of the guided surface wave receive structures R were obtaining the maximum amount of energy for which they were rated. Other expected energy flow maps 2119 can represent the expected flow of energy if the guided surface wave receive structures R were obtaining various amounts of energy. In addition, expected energy flow maps 2119 can be created for particular frequencies and/or frequency ranges at which the guided surface waveguide probes P and/or guided surface wave receive structures R operate.
The energy usage monitor 2116 can generate expected energy flow maps 2119 in various ways. In one approach, an expected energy flow map 2119 is based on one or maps depicting actual energy flow that were created in the past or that were created using historical data. For example, the energy usage monitor 2116 can generate such a map using historical data from times when it was known to a high degree of certainty that no unauthorized guided surface wave receive structures R were operating in the geographic region. In other words, the resulting expected energy flow map 2119 can represent conditions that existed when the geographical region was in a known operating state.
To generate an expected energy flow map 2119 based on historical data, the energy usage monitor 2116 can obtain the electric field vectors Ê and the magnetic field vectors Ĥ as determined by the metering systems 2106 at the corresponding historical time. Using this data, the energy usage monitor 2116 can then calculate energy flow vectors at the location of each metering system 2106 based on the electric field vectors Ê and the magnetic field vectors Ĥ. Each energy flow vector can, for example, indicate the magnitude and direction of the energy flux density at the location of the field meter M. Accordingly, in some examples, an energy flow vector is a Poynting vector. The Poynting vector can be calculated as the cross product of the electric field vector and the magnetic field vector, which can be represented as:

Ŝ=Ê×Ĥ  (106)

where Ŝ is the Poynting vector, Ê is the electric field vector, and Ĥ is the magnetic field vector. Thus, by using Equation 106, for example, the energy usage monitor 2116 can generate an energy flow vector that indicates the magnitude and direction of the flow of energy at the corresponding location.

Once the energy usage monitor 2116 has determined the energy flow vectors corresponding to the metering systems 2106 at the historic time, the energy usage monitor 2116 can plot the energy flow vectors on a map using the known locations of the metering systems 2106. The result is an energy flow map based on historical data from when the conditions of the geographical area were in a known operating state. The energy usage monitor 2116 can then store this map in the data store 2113 as an expected energy flow map 2119, because the expected energy flow in the geographical region can be expected to be substantially similar if similar operating conditions exist at a future time. In some examples, the process of creating an energy flow map described above can be performed multiple times, and the energy usage monitor 2116 can average the results of the multiple maps to create the expected energy flow map 2119.
In another approach of creating an expected energy flow map 2119, the energy usage monitor 2116 generates an expected energy flow map 2119 based on data that represents and/or corresponds to the amount of energy being consumed by each guided surface wave receive structure R. For example, the energy usage monitor 2116 can communicate with the receive structure controllers 2105 to obtain data that represents the location of the guided surface wave receive structures R and the instantaneous power being delivered to the guided surface wave receive structures R. Because the power being delivered to a guided surface wave receive structure R corresponds to the amount of energy being obtained by the guided surface wave receive structure R, the energy usage monitor 2116 can apply this information to plot the relative direction and magnitude of energy flow at various locations in the geographical region.
For instance, if there are two guided surface wave receive structures R in the geographical region, and if a much greater amount of power is being delivered to the first guided surface wave receive structure R power than the second guided surface wave receive structure R, much more energy would be flowing to the first guided surface wave receive structure R relative to the amount of energy that would be flowing to the second guided surface wave receive structure R. Accordingly, the energy usage monitor 2116 can construct an expected energy flow map 2119 indicating that a relatively large amount of energy is expected to flow to the first guided surface wave receive structure R, while a relatively small amount of energy is expect to flow to the second guided surface wave receive structure R. The energy usage monitor 2116 can plot corresponding information for all of the guided surface wave receive structures R of which the energy usage monitor 2116 has information. Thus, the energy usage monitor 2116 can generate an expected energy flow map 2119 based on energy consumption data associated with one or more guided surface wave receive structures R.
In another approach of creating an expected energy flow map 2119, the energy usage monitor 2116 generates an expected energy flow map 211 by running simulations of the guided surface waveguide probes P and/or guided surface wave receive structures R in a geographical region. To this end, the energy usage monitor 2116 can include a simulator system that simulates operation of the guided surface waveguide probes P and/or guided surface wave receive structures R based on the known locations of guided surface waveguide probes P and guided surface wave receive structures R. The simulator system can determine expected electric fields, magnetic fields, energy flow vectors, and/or other information based on various input parameters. Examples of input parameters for the simulator system include the number of end users that are obtaining power from the respective guided surface wave receive structures R, the amount of current being obtained by the respective guided surface waveguide receive structures R, voltages associated with various components of the respective guided surface wave receive structures R, and/or other parameters. In some examples, the receive structure controllers 2105 can provide this information to the energy usage monitor 2116.
After the energy usage monitor 2116 has simulated the operation of the guided surface waveguide probes P and/or the guided surface wave receive structures R to determine expected electric fields, magnetic fields, and/or energy flow vectors, the energy usage monitor 2116 can plot the expected energy flow vectors on a map. The result is an expected energy flow map 2119 based on a simulation of the operation of the guided surface waveguide probes P and/or the guided surface wave receive structures R. In some examples, the energy usage monitor 2116 can perform simulations based on real time data to generate an expected energy flow map 2116 that corresponds to the present operating conditions of the guided surface waveguide probes P and/or guided surface wave receive structures R.
After one or more expected energy flow maps 2119 have been generated for a geographical region associated with the energy usage monitor 2116, the energy usage monitor 2116 can begin the process of determining whether an unauthorized guided surface wave receive structure R is operating in the geographic region. To this end, the energy usage monitor 2116 can obtain the actual electric field vectors Ê and the magnetic field vectors Ĥ as determined by the metering systems 2106. Using this data, the energy usage monitor 2116 can then calculate actual energy flow vectors at the location of each metering system 2106 based on the electric field vectors Ê and the magnetic field vectors Ĥ. For example, the energy usage monitor 2116 can use Equation 106 to calculate actual energy flow vectors for the location of each metering system 2106. As discussed above, each actual energy flow vector can indicate the magnitude and direction of the flow of energy at the corresponding location.
Once the energy usage monitor 2116 has determined the actual energy flow vectors, the energy usage monitor 2116 can plot the actual energy flow vectors on the corresponding locations on map and then store the map as an actual energy flow map 2123. As discussed above, the actual energy flow map 2123 can represent the flow of electrical energy that has been determined to have actually occurred at particular locations within a geographical region. In some examples, the actual energy flow map 2123 is embodied as a graphical map in which the energy flow vectors are graphically depicted on a map.
After the energy usage monitor 2116 has created an actual energy flow map 2123, the energy usage monitor 2116 can compare the actual energy flow map 2123 to an expected energy flow map 2119 that was created using one or more of the techniques described above. For embodiments in which multiple expected energy flow maps 2119 corresponding to various scenarios are stored in the data store 2113, the energy usage monitor 2116 can select the particular expected energy flow map 2119 that corresponds to the expected operating conditions that existed for the actual energy flow map 2123. For example, if it is believed that all guided surface wave receive structures R were obtaining the maximum amount of energy for which the guided surface wave receive structures R were rated, the energy usage monitor 2116 can select the expected energy flow map 2119 that corresponds to such a scenario.
The energy usage monitor 2116 can compare the actual energy flow map 2123 to the expected energy flow map 2119 to determine whether there is a discrepancy between the actual energy flow map 2123 and the expected energy flow map 2119. The energy usage monitor 2116 can determine that a discrepancy exists, for example, if the difference between one or more energy flow vectors represented in the actual energy flow map 2123 and one or more energy flow vectors represented in the expected energy flow map 2119 is greater than a particular threshold. For example, if the difference in magnitude and/or direction between an energy flow vector for a location in an expected energy flow map 2119 and an energy flow vector for the same location in an actual energy flow map 2123 is exceeds a particular threshold, the energy usage monitor 2116 can determine that a discrepancy exists between the expected energy flow map 2119 and the actual energy flow map 2123.
When a discrepancy between the expected energy flow map 2119 and the actual energy flow map 2123 is detected, the energy usage monitor 2116 can conclude that an unauthorized guided surface wave receive structure R is operating in the geographical region associated with the energy usage monitor. In addition, because the energy flow vectors in the actual energy flow map 2123 can indicate the direction of energy flow at various locations, the actual energy flow map 2123 can be used to identify the location where the unauthorized guided surface wave receive structure R is located. For example, in examples in which the metering systems 2106 form a grid, the actual energy flow map 2123 can be used to identify the particular grid segment in which the unauthorized guided surface wave receive structure R is located, as will be described in further detail below.
The energy usage monitor 2116 can initiate an action in response to detecting that an unauthorized guided surface wave receive structure R is operating. For example, the energy usage monitor 2116 can command with the probe controllers 2104 to command one or more guided surface waveguide probes P to stop generating guided surface waves. As another example, the energy usage monitor 2116 can transmit a message to notify law enforcement or other individuals that an unauthorized guided surface wave receive structure R has been detected. Furthermore, the energy usage monitor 2116 can request one or more metering systems 2106 to send updated data regarding their operation so that the energy usage monitor 2116 can generate an updated actual energy flow map 2123 in order to verify whether the previous actual energy flow map 2123 was accurate.
With reference to FIG. 22, shown is a scenario depicting a geographical region 2203 in which a guided surface waveguide probe P and an authorized guided surface wave receive structure R, referred to herein as the authorized guided surface wave receive structure RAuth, is operating. In addition, metering systems 2106, referred to herein as the metering systems 2106a-2106i, are located in the geographical region 2203. In the example shown, the metering systems 2106a-2106i are arranged to form a grid in the geographical region 2203.
In the scenario depicted in FIG. 22, the authorized guided surface wave receive structure RAuth is obtaining electrical energy from guided surface waves launched by the guided surface waveguide probe P. The energy usage monitor 2116 is aware that the authorized guided surface wave receive structure RAuth is operating and expects the authorized guided surface wave receive structure RAuth to be affect the depicted flow of energy an actual energy flow map 2123 that may be generated for the situation depicted in FIG. 22.
However, as illustrated in FIG. 22, an unauthorized guided surface wave receive structure R, referred to herein as the unauthorized guided surface wave receive structure RUnauth, is also operating in the geographical region 2203 and obtaining electrical energy from guided surface waves launched by the guided surface waveguide probe P. The unauthorized guided surface wave receive structure RUnauth is obtaining electrical energy from the guided surface waveguide probe P without authorization of the entity that is operating the guided surface waveguide probe P. For the scenario depicted in FIG. 22, the energy usage monitor 2116 is not aware that the unauthorized guided surface wave receive structure RUnauth is operating in the geographical region 2203. Accordingly, the energy usage monitor 2116 would not expect for energy flow vectors in the actual energy flow map 2123 to indicate that the unauthorized guided surface wave receive structure RUnauth exists.
With reference to FIG. 23A, shown is a rendering of an example of an expected energy flow map 2119 that corresponds to the scenario in which the authorized guided surface wave receive structure RAuth is the only guided surface wave receive structure R operating in the geographical region 2203. It is noted that in other examples in which there are a larger quantity of guided surface wave receive structures R, the complexity of the expected energy flow map 2119 may be more complex. In addition, the quantity of metering systems 2106 can be increased to provide an expected energy flow map 2119 depicting a higher resolution of the expected energy flow.
As shown in FIG. 23A, the expected energy flow map 2119 includes arrows 2303a-2303i that represent expected energy flow vectors for the geographical region 2203. For the example illustrated in FIG. 23A, the lengths of the arrows 2303a-2303i correspond to the relative magnitudes of the expected energy flow vectors, and the direction of the arrows indicate the direction of the expected energy flow vectors. In addition, to aid the present discussion, the location of the authorized guided surface wave receive structure RAuth is plotted on the expected energy flow map 2119.
The expected energy flow map 2119 shown in FIG. 23A corresponds to the situation in which the authorized guided surface wave receive structure RAuth is the only guided surface wave receive structure R operating in the geographical region 2203. As such, all of the arrows 2303a-2303i point to the location of the authorized guided surface wave receive structure RAuth. In addition, the arrows 2303a, 2303b, 2303d, and 2303e, which are the arrows 2303a-2303i that are closest to the authorized guided surface wave receive structure RAuth, are longer than the other arrows 2303c, 2303f, and 2303g-2303i, indicating that the magnitudes of the expected energy flow vectors near the authorized guided surface wave receive structure RAuth are larger than at other locations.
With reference to FIG. 23B, shown is a rendering of an example of an actual energy flow map 2123 of the geographical region 2203 that corresponds to the scenario illustrated in FIG. 22. For the example shown in FIG. 23B, the actual energy flow map 2123 includes arrows 2306a-2306i that represent actual energy flow vectors determined for the geographical region 2203. The lengths of the arrows 2306a-2306i correspond to the relative magnitudes of the actual energy flow vectors, and the direction of the arrows indicate the direction of the actual energy flow vectors. In addition, to aid the present discussion, the locations of the authorized guided surface wave receive structure RAuth and the unauthorized guided surface wave receive structure RUnauth are plotted on the act energy flow map 2123.
As described above, the energy usage monitor 2116 can compare the expected energy flow vector map 2119 that is rendered in FIG. 23A with the actual energy flow vector map 2123 that is rendered in FIG. 23B to determine whether there is a discrepancy between the expected energy flow map 2119 and the actual energy flow map 2123. For example, the energy usage monitor 2116 can compare individual ones of the expected energy flow vectors, which are represented by the arrows 2303a-2303i in FIG. 23A, to the corresponding actual energy flow vectors, which are represented by the arrows 2306a-2306i in FIG. 23B.
If the difference in magnitude or direction between one or more of the expected energy flow vectors and the corresponding actual energy flow vectors exceeds a particular amount, the energy usage monitor 2116 can conclude that a guided surface waveguide receive structure R is operating without authorization in the geographical region 2203. For example, because the magnitude of the actual energy flow vector that is represented by the arrow 2306i in FIG. 23B is much greater than the expected energy flow vector that is represented by the arrow 2303i in FIG. 23A, the energy usage monitor 2116 concludes that there is a guided surface wave receive structure R operating in the geographical region 2203 without authorization. In addition, because the direction of the actual energy flow vector represented by the arrow 2306h in FIG. 23B is much different from the direction of the actual energy flow vector represented by the arrow 2303h in FIG. 23A, the energy usage monitor 2116 concludes that there is a guided surface wave receive structure R operating in the geographical region 2203 without authorization.
As described above, the energy usage monitor 2116 can also use the actual energy flow map 2306d to identify the location of the unauthorized guided surface wave receive structure RUnauth. For example, in the scenario depicted in FIGS. 23A-23B, the arrows 2306f, 2306h, and 2306i in FIG. 23B, which represent actual energy flow vectors, point towards the location of the unauthorized guided surface wave receive structure RAuth. By pointing towards the unauthorized guided surface wave receive structure RUnauth, the arrows 2306f, 2306h, and 2306i indicate that a relatively large amount of energy is flowing towards the unauthorized guided surface wave receive structure RAuth. In addition, the length of the arrow 2306e in FIG. 23B, which represents the actual energy flow, is much smaller than the length of the arrow 2303e in FIG. 23A, which represents the expected energy flow at that location. Accordingly, arrows 2306e, 2306f, 2036h, and 2306i for the actual energy flow map 2123 indicate that a relatively large amount of energy is flowing to the grid segment defined by the metering systems 2106e, 2106f, 2106h, and 2106i. As a result, the energy usage monitor 2116 concludes that the unauthorized guided surface wave receive structure RUnauth is operating in this grid segment.
As described above, the energy usage monitor 2116 can initiate an action in response to detecting that the unauthorized guided surface wave receive structure RUnauth. For example, the energy usage monitor 2116 can communicate with a probe controller 2104 command the guided surface waveguide probe P operating in the geographical region 2203 to stop generating guided surface waves. As another example, the energy usage monitor 2116 can transmit a message to notify law enforcement or other individuals that the unauthorized guided surface wave receive structure RUnauth has been detected. The message can also specify the location of the unauthorized guided surface wave receive structure RUnauth that was detected.
With reference to FIG. 24, shown is a flowchart that depicts an example of the operation of a portion of the energy usage monitor 2116 according to various embodiments. In particular, the flowchart of FIG. 24 depicts an example of the energy usage monitor 2116 generating an expected energy flow map 2119 based on actual energy consumption data associated with guided surface wave receive structures R. The flowchart of FIG. 24 provides an example of the many types of functional arrangements that can be employed to implement the operation of the energy usage monitor 2116 as described herein. As an alternative, the flowchart of FIG. 24 may be viewed as depicting an example of elements of a method implemented by the computing environment 2103.
Beginning at box 2403, the energy usage monitor 2116 obtains geographic location data for the guided surface waveguide receive structures R operating in the geographical region being monitored by the energy usage monitor 2116. The geographic location data can be provided to the energy usage monitor 2116 by the receive structure controllers 2105. In some examples, the geographic location data specifies the locations of the guided surface wave receive structures R.
At box 2406, the energy usage monitor 2116 obtains actual energy consumption data that is associated with the guided surface wave receive structures R. The actual energy consumption data can be provided to the energy usage monitor 2116 by the receive structure controllers 2105. The actual energy consumption data can represent or indicate the amount of energy being obtained by individual ones of the guided surface wave receive structure R. For example, the actual energy consumption data for a guided surface wave receive structure R can represent the amount of energy, power, and/or current being obtained by the guided surface wave receive structure R.
As shown at box 2409, the energy usage monitor 2116 then generates expected energy flow vectors based on the actual energy consumption data. For instance, if the actual energy consumption data indicates that a relatively large amount of power is being delivered to a particular guided surface wave receive structure R, the energy usage monitor 2116 can generate expected energy flow vectors depicting that a relatively large amount of energy is expected to be flowing to the guided surface wave receive structure R. By contrast, if the actual energy consumption data indicates that a relatively small amount of power is being delivery to another guided surface wave receive structure R, the energy usage monitor 2116 can generate expected energy flow vectors depicting that a relatively small amount of energy is expected to be flowing to that guided surface wave receive structure R.
As indicated at box 2413, the energy usage monitor 2116 then plots the generated energy flow vectors on a map based on the geographic location data for the guided surface wave receive structures R. For example, if the energy usage monitor 2116 determines an expected energy flow vector at the location of a guided surface wave receive structure R, the energy usage monitor 2116 can plot the expected energy flow vector where the obtained geographic location data indicates the guided surface wave receive structure R is located.
The energy usage monitor 2116 then stores the generated map as an expected energy flow map 2119 in the data store 2113 (FIG. 21) at box 2416. Thereafter, the process ends.
With reference to FIG. 25, shown is a flowchart that depicts another example of the operation of a portion of the energy usage monitor 2116 according to various embodiments. In particular, the flowchart of FIG. 25 depicts an example of the energy usage monitor 2116 generating an expected energy flow map 2119 based on simulating the operation of one or more guided surface waveguide probes P and one or more guided surface wave receive structures R in the region monitored by the energy usage monitor 2116. The flowchart of FIG. 25 provides an example of the many types of functional arrangements that can be employed to implement the operation of the energy usage monitor 2116 as described herein. As an alternative, the flowchart of FIG. 25 may be viewed as depicting an example of elements of a method implemented by the computing environment 2103.
Beginning at box 2503, the energy usage monitor 2116 simulates the operation of the one or more guided surface waveguide probes P and/or one or more guided surface wave receive structures that are present in the region monitored by the energy usage monitor 2116. To this end, the energy usage monitor 2116 can include a simulator system that simulates operation of the guided surface waveguide probes P and/or guided surface wave receive structures R. The simulator system can determine expected electric fields, magnetic fields, energy flow vectors, and/or other information based on various input parameters, as described above.
As shown at box 2506, the energy usage monitor 2116 then calculates simulated electric field vectors for locations in the geographical region monitored by the energy usage monitor 2116. In some examples, the energy usage monitor 2116 calculates simulated electric field vectors for the locations of the metering systems 2106.
As shown at box 2509, the energy usage monitor 2116 calculates simulated magnetic field vectors for locations in the region monitored by the energy usage monitor 2116. In some examples, the energy usage monitor 2116 calculates simulated magnetic field vectors for the locations of the metering systems 2106. The energy usage monitor 2116 then moves to box 2513 and calculates expected energy flow vectors based on the simulated electric and magnetic field vectors. For example, the energy usage monitor 2116 can use Equation 106 to determine the expected energy flow vectors based on the simulated electric and magnetic field vectors.
As shown at box 2516, the energy usage monitor 2116 then plots the generated energy flow vectors on a map. Thereafter, the energy usage monitor 2116 stores the generated map as an expected energy flow map 2119 in the data store 2113 (FIG. 21), as shown at box 2519. Then, the process ends.
With reference to FIG. 26, shown is a flowchart that depicts another example of the operation of a portion of the energy usage monitor 2116 according to various embodiments. In particular, the flowchart of FIG. 26 depicts an example of the energy usage monitor 2116 generating an expected energy flow map 2119 based on electric and magnetic field readings obtained from metering systems 2106 when it is known to a high degree of certainty that only authorized guided surface wave receive structures R are operating in the geographic region monitored by the energy usage monitor 2116. The flowchart of FIG. 26 provides an example of the many types of functional arrangements that can be employed to implement the operation of the energy usage monitor 2116 as described herein. As an alternative, the flowchart of FIG. 26 may be viewed as depicting an example of elements of a method implemented by the computing environment 2103.
Beginning at box 2603, the energy usage monitor 2116 communicates with the probe controllers 2104 to configure the guided surface waveguide probes P to launch guided surface waves using one or more of the techniques described above. At box 2606, the energy usage monitor 2116 also communicates with the receive structure controllers 2105 to configure the guided surface wave receive structures R to obtain electrical energy from the launched guided surface waves. At this point, it is assumed that there are no unauthorized guided surface wave receive structures R operating in the geographic region.
As shown at box 2609, the energy usage monitor 2116 then obtains electric field readings from each metering system 2106. For example, for each metering system 2106, the energy usage monitor 2116 can obtain data representing the magnitude of the electric field in each of the {circumflex over (x)}, ŷ, and {circumflex over (z)} directions. Using this data, the energy usage monitor 2116 determines an electric field vector for each metering system 2106, as shown at box 2613. Each electric field vector can indicate the magnitude and direction of the electric field at the location of the corresponding metering system 2106.
At box 2609, the energy usage meter 2116 obtains magnetic field readings from each metering system 2106. For example, for each metering system 2106, the energy usage monitor 2116 can obtain data representing the magnitude of the electric field in each of the {circumflex over (x)}, ŷ, and {circumflex over (z)} directions. Using this data, the energy usage monitor 2116 determines a magnetic field vector for each metering system 2106, as shown at box 2619. Each magnetic field vector can indicate the magnitude and direction of the magnetic field at the location of the corresponding metering system 2106.
As indicated at box 2623, the energy usage monitor 2116 then calculates an energy flow vector for each metering system 2106 based on the determined electric and magnetic field vectors. For example, the energy usage monitor 2116 can use Equation 106 to determine a Poynting vector for each metering system 2106.
The energy usage monitor 2116 then plots each energy flow vector at the corresponding location on a map, as shown at box 2626, thereby generating an expected energy flow map 2119. At box 2629, the energy usage monitor can then store the expected energy flow map 2119 in the data store 2113. Thereafter, the process ends.
With reference to FIGS. 27A-27B, shown is a flowchart that depicts another example of the operation of a portion of the energy usage monitor 2116 according to various embodiments. In particular, the flowchart of FIGS. 27A-27B depicts an example of the energy usage monitor 2116 detecting whether an unauthorized guided surface wave receive structure R. The flowchart of FIG. 26 provides an example of the many types of functional arrangements that can be employed to implement the operation of the energy usage monitor 2116 as described herein. As an alternative, the flowchart of FIG. 26 may be viewed as depicting an example of elements of a method implemented by the computing environment 2103.
Beginning with box 2703, the energy usage monitor 2116 obtains electric field readings from each metering system 2106. For example, for each metering system 2106, the energy usage monitor 2116 can obtain data representing the magnitude of the electric field in each of the {circumflex over (x)}, ŷ, and {circumflex over (z)} directions. Using this data, the energy usage monitor 2116 determines an electric field vector for each metering system 2106, as shown at box 2706. Each electric field vector can indicate the magnitude and direction of the electric field at the location of the corresponding metering system 2106.
At box 2709, the energy usage meter 2116 obtains magnetic field readings from each metering system 2106. For example, for each metering system 2106, the energy usage monitor 2116 can obtain data representing the magnitude of the electric field in each of the {circumflex over (x)}, ŷ, and {circumflex over (z)} directions. Using this data, the energy usage monitor 2116 determines a magnetic field vector for each metering system 2106, as shown at box 2713. Each magnetic field vector can indicate the magnitude and direction of the magnetic field at the location of the corresponding metering system 2106.
As indicated at box 2716, the energy usage monitor 2116 then calculates an energy flow vector for each metering system 2106 based on the determined electric and magnetic field vectors. For example, the energy usage monitor 2116 can use Equation 106 to determine a Poynting vector Ŝ for each metering system 2106. The energy usage monitor 2116 can then plot each energy flow vector at the corresponding location on a map to generate an actual energy flow map 2123, as indicated at box 2719.
The energy usage monitor 2116 can then proceed to box 2723, which is shown on FIG. 27B. At box 2723, the energy usage monitor 2116 can retrieve an expected energy flow map 2119 stored in the data store 2113. In some examples, multiple expected energy flow maps 2119 are stored in the data store 2113, and the energy usage monitor 2116 can select a particular expected energy flow map 2119 that corresponds to the expected operating conditions for the geographical region, as discussed above.
At box 2726, the energy usage monitor 2116 compares the actual energy flow map 2123 to the retrieved expected energy flow map 2119. For example, the energy usage monitor 2116 can compare individual ones of the expected energy flow vectors in the expected energy flow map 2119 to corresponding actual energy flow vectors in the actual energy flow map 2123.
As shown at box 2729, the energy flow monitor 2116 can detect whether there is a discrepancy between the expected energy flow map 2119 and the actual energy flow map 2123. For example, if the difference in magnitude and/or direction between an energy flow vector for a location in an expected energy flow map 2119 and an energy flow vector for the same location in an actual energy flow map 2123 exceeds a particular threshold, the energy usage monitor 2116 can determine that a discrepancy exists between the expected energy flow map 2119 and the actual energy flow map 2123.
If a discrepancy is not detected between the expected energy flow map 2119 and the actual energy flow map 2123, the process ends as shown. Otherwise, if a discrepancy is detected, the energy usage monitor 2116 determines that an unauthorized guided surface wave receive structure is operating, as shown at box 2733.
As shown at box 2736, the energy usage monitor 2116 can also initiate an action in response to detecting a discrepancy. For example, the energy usage monitor 2116 can communicate with one or more probe controllers 2104 to command one or more guided surface waveguide probes P to stop generating guided surface waves. As another example, the energy usage monitor 2116 can transmit a message to notify law enforcement or other individuals that an unauthorized guided surface wave receive structure R has been detected. Furthermore, the energy usage monitor 2116 can request one or more metering systems 2106 to send updated data regarding their operation so that the energy usage monitor 2116 can generate an updated actual energy flow map 2123 in order to verify whether the previous actual energy flow map 2123 was accurate. After box 2736, the process ends.
With reference to FIG. 28, shown is a schematic block diagram of the computing environment 2103 according to an embodiment of the present disclosure. The computing environment 103 includes one or more computing devices 2803. Each computing device 2803 includes at least one processor circuit, for example, having a processor 2806 and a memory 2809, both of which are coupled to a local interface 2813. To this end, each computing device 2803 may comprise, for example, at least one server computer or like device. The local interface 2813 may comprise, for example, a data bus with an accompanying address/control bus or other bus structure as can be appreciated.
Stored in the memory 2809 are both data and several components that are executable by the processor 2806. In particular, stored in the memory 2809 and executable by the processor 2806 is the energy usage monitor 2116 and potentially other components. Also stored in the memory 2809 may be the data store 2113 and other data. In addition, an operating system may be stored in the memory 2809 and be executable by the processor 2806.
It is understood that there may be other components that are stored in the memory 2809 and are executable by the processor 2806 as can be appreciated. Where any component discussed herein is implemented in the form of software, any one of a number of programming languages may be employed such as, for example, C, C++, C≠, Objective C, Java®, JavaScript®, Pen, PHP, Visual Basic®, Python®, Ruby, Flash®, or other programming languages.
A number of software components are stored in the memory 2809 and are executable by the processor 2806. In this respect, the term “executable” means a program file that is in a form that can ultimately be run by the processor 2806. Examples of executable programs may be, for example, a compiled program that can be translated into machine code in a format that can be loaded into a random access portion of the memory 2809 and run by the processor 2806, source code that may be expressed in proper format such as object code that is capable of being loaded into a random access portion of the memory 2809 and executed by the processor 2806, or source code that may be interpreted by another executable program to generate instructions in a random access portion of the memory 2809 to be executed by the processor 2806, etc. An executable program may be stored in any portion or component of the memory 2809 including, for example, random access memory (RAM), read-only memory (ROM), hard drive, solid-state drive, USB flash drive, memory card, optical disc such as compact disc (CD) or digital versatile disc (DVD), floppy disk, magnetic tape, or other memory components.
The memory 2809 is defined herein as including both volatile and nonvolatile memory and data storage components. Volatile components are those that do not retain data values upon loss of power. Nonvolatile components are those that retain data upon a loss of power. Thus, the memory 2809 may comprise, for example, random access memory (RAM), read-only memory (ROM), hard disk drives, solid-state drives, USB flash drives, memory cards accessed via a memory card reader, floppy disks accessed via an associated floppy disk drive, optical discs accessed via an optical disc drive, magnetic tapes accessed via an appropriate tape drive, and/or other memory components, or a combination of any two or more of these memory components. In addition, the RAM may comprise, for example, static random access memory (SRAM), dynamic random access memory (DRAM), or magnetic random access memory (MRAM) and other such devices. The ROM may comprise, for example, a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or other like memory device.
Also, the processor 2806 may represent multiple processors 2806 and/or multiple processor cores and the memory 2809 may represent multiple memories 2809 that operate in parallel processing circuits, respectively. In such a case, the local interface 2813 may be an appropriate network that facilitates communication between any two of the multiple processors 2806, between any processor 2806 and any of the memories 2809, or between any two of the memories 2809, etc. The local interface 2813 may comprise additional systems designed to coordinate this communication, including, for example, performing load balancing. The processor 2806 may be of electrical or of some other available construction.
Although the energy usage monitor 2116, and other various systems described herein may be embodied in software or code executed by hardware as discussed above, as an alternative the same may also be embodied in dedicated hardware or a combination of software/hardware and dedicated hardware. If embodied in dedicated hardware, each can be implemented as a circuit or state machine that employs any one of or a combination of a number of technologies. These technologies may include, but are not limited to, discrete logic circuits having logic gates for implementing various logic functions upon an application of one or more data signals, application specific integrated circuits (ASICs) having appropriate logic gates, field-programmable gate arrays (FPGAs), or other components, etc. Such technologies are generally well known by those skilled in the art and, consequently, are not described in detail herein.
The flowcharts of FIGS. 24-26 and 27A-27B show the functionality and operation of an implementation of portions of the energy usage monitor 2116. If embodied in software, each block may represent a module, segment, or portion of code that comprises program instructions to implement the specified logical function(s). The program instructions may be embodied in the form of source code that comprises human-readable statements written in a programming language or machine code that comprises numerical instructions recognizable by a suitable execution system such as a processor 2806 in a computer system or other system. The machine code may be converted from the source code, etc. If embodied in hardware, each block may represent a circuit or a number of interconnected circuits to implement the specified logical function(s).
Although the flowcharts of FIGS. 24-26 and 27A-27B a specific order of execution, it is understood that the order of execution may differ from that which is depicted. For example, the order of execution of two or more blocks may be switched relative to the order shown. Also, two or more blocks shown in succession in FIGS. 24-26 and 27A-27B may be executed concurrently or with partial concurrence. Further, in some embodiments, one or more of the blocks shown in FIGS. 24-26 and 27A-27B may be skipped or omitted. In addition, any number of counters, state variables, warning semaphores, or messages might be added to the logical flow described herein, for purposes of enhanced utility, accounting, performance measurement, or providing troubleshooting aids, etc. It is understood that all such variations are within the scope of the present disclosure.
Also, any logic or application described herein, including the energy usage monitor 2116, that comprises software or code can be embodied in any non-transitory computer-readable medium for use by or in connection with an instruction execution system such as, for example, a processor 2806 in a computer system or other system. In this sense, the logic may comprise, for example, statements including instructions and declarations that can be fetched from the computer-readable medium and executed by the instruction execution system. In the context of the present disclosure, a “computer-readable medium” can be any medium that can contain, store, or maintain the logic or application described herein for use by or in connection with the instruction execution system.
The computer-readable medium can comprise any one of many physical media such as, for example, magnetic, optical, or semiconductor media. More specific examples of a suitable computer-readable medium would include, but are not limited to, magnetic tapes, magnetic floppy diskettes, magnetic hard drives, memory cards, solid-state drives, USB flash drives, or optical discs. Also, the computer-readable medium may be a random access memory (RAM) including, for example, static random access memory (SRAM) and dynamic random access memory (DRAM), or magnetic random access memory (MRAM). In addition, the computer-readable medium may be a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or other type of memory device.
Further, any logic or component described herein, including the energy usage monitor 2116, may be implemented and structured in a variety of ways. For example, one or more components described may be implemented as modules or components of a single application. Further, one or more components described herein may be executed in shared or separate computing devices or a combination thereof. For example, a plurality of the components described herein may execute in the same computing device 2803 or in multiple computing devices 2803 in the same computing environment 2103. Additionally, it is understood that terms such as “application,” “service,” “system,” “engine,” “module,” and so on may be interchangeable and are not intended to be limiting.
Disjunctive language such as the phrase “at least one of X, Y, or Z,” unless specifically stated otherwise, is otherwise understood with the context as used in general to present that an item, term, etc., may be either X, Y, or Z, or any combination thereof (e.g., X, Y, and/or Z). Thus, such disjunctive language is not generally intended to, and should not, imply that certain embodiments require at least one of X, at least one of Y, or at least one of Z to each be present.
The above-described embodiments of the present disclosure are merely possible examples of implementations set forth for a clear understanding of the principles of the disclosure. Many variations and modifications may be made to the above-described embodiment(s) without departing substantially from the spirit and principles of the disclosure. All such modifications and variations are intended to be included herein within the scope of this disclosure and protected by the following claims. In addition, all optional and preferred features and modifications of the described embodiments and dependent claims are usable in all aspects of the disclosure taught herein. Furthermore, the individual features of the dependent claims, as well as all optional and preferred features and modifications of the described embodiments are combinable and interchangeable with one another.

## Claims

1. A system, comprising:
a guided surface waveguide probe configured to launch a guided surface wave along a surface of a terrestrial medium;
a plurality of metering systems distributed within a geographical region associated with the guided surface waveguide probe;
at least one computing device; and
memory storing a plurality of computer instructions that, when executed, cause the at least one computing device to at least generate an energy flow map using at least data obtained from the plurality of metering systems.

2. The system of claim 1, wherein the plurality of computer instructions further cause the at least one computing device to at least:
compare the energy flow map to an expected energy flow map; and
detect that a comparison of the energy flow map to the expected energy flow map indicates that an unexpected guided surface wave receive structure exists in the geographical region.

3. The system of claim 2, wherein the plurality of computer instructions further cause the at least one computing device to at least generate the expected energy flow map based on energy consumption data associated with a guided surface wave receive structure.

4. The system of claim 2, wherein the plurality of computer instructions further cause the at least one computing device to at least:
calculate at least one expected electric field vector;
calculate at least one expected magnetic field vector;
calculate at least one energy flow vector based on the at least one expected electric field vector and the at least one expected magnetic field vector; and
generate the expected energy flow map based on the at least one expected magnetic field vector.

5. The system of claim 2, wherein the plurality of computer instructions further cause the at least one computing device to at least generate the expected energy flow map based on historical energy flow data for the geographical region.

6. The system of claim 2, wherein the plurality of computer instructions further cause the at least one computing device to at least:
calculate an energy flow vector associated with at least one of the plurality of metering systems, wherein the energy flow vector represents an energy magnitude and an energy direction associated with the at least one of the plurality of metering systems; and
generate the energy flow map based on the energy flow vector.

7. The system of claim 6, wherein the energy flow vector is a Poynting vector.

8. The system of claim 6, wherein the plurality of computer instructions further cause the at least one computing device to at least calculate a cross product of an electric field vector and a magnetic field vector associated with the at least one of the plurality of metering systems to calculate the energy flow vector.

9. A method, comprising:
launching, using a guided surface waveguide probe, a guided surface wave along a surface of a terrestrial medium;
obtaining, by at least one computing device, data from a plurality of metering systems distributed within a geographical region associated with the guided surface waveguide probe; and
generating, by the at least one computing device, an energy flow map based on the data from the plurality of metering systems.

10. The method of claim 9, further comprising:
comparing, by the at least one computing device, the energy flow map to an expected energy flow map; and
detecting, by the at least one computing device, that a comparison of the energy flow map to the expected energy flow map indicates that an unexpected guided surface wave receive structure exists in the geographical region.

11. The method of claim 10, further comprising generating, by the at least one computing device, the expected energy flow map based on energy consumption data associated with the geographical region.

12. The method of claim 10, further comprising:
calculating, by the at least one computing device, at least one expected electric field vector;
calculating, by the at least one computing device, at least one expected magnetic field vector;
calculating, by the at least one computing device, at least one energy flow vector based on the at least one expected electric field vector and the at least one expected magnetic field vector; and
generating, by the at least one computing device, the expected energy flow map based on the at least one expected energy flow vector.

13. The method of claim 10, further comprising generating, by the at least one computing device, the expected energy flow map based on historical data for the geographical region.

14. The method of claim 9, further comprising:
calculating, by the at least one computing device, an energy flow vector associated at least one of the plurality of metering systems; and
generating, by the at least one computing device, the energy flow map based on the energy flow vector.

15. The method of claim 14, wherein the energy flow vector is a Poynting vector.

16. The method of claim 14, further comprising calculating, by the at least one computing device, the energy flow vector by at least calculating a cross product of an electric field vector and a magnetic field vector associated with the at least one of the plurality of metering systems.

17. A system, comprising:
a guided surface waveguide probe configured to launch a guided surface wave along a surface of a terrestrial medium;
a plurality of meters;
at least one computing device; and
memory storing a plurality of computer instructions that, when executed, cause the at least one computing device to at least generate an energy flow map based on electric field data and magnetic field data obtained from the plurality of meters.

18. The system of claim 17, wherein the plurality of computer instructions further cause the computing device to at least:
compare the energy flow map to an expected energy flow map; and
detect that a comparison of the energy flow map to the expected energy flow map indicates that an unexpected guided surface wave receive structure exists in a geographical region in which the guided surface waveguide probe operates.

19. The system of claim 18, wherein the plurality of computer instructions further cause the at least one computing device to at least generate the expected energy flow map based on energy consumption data associated with a guided surface wave receive structure.

20. The system of claim 18, wherein the plurality of computer instructions further cause the at least one computing device to at least generate the expected energy flow map based on historical data for a geographical region in which the guided surface waveguide probe operates.