# Patent

## Title
Broadband vibrational energy harvesting

## Abstract

A system that converts environmental vibrational energy into electrical energy includes a transducer that undergoes oscillating movement in response to the vibrational energy in order to produce an oscillating electrical signal. Power electronics process the oscillating electrical signal. A control system (including at least one control element of the power electronics, at least one sensor and control electronics) carries out a control scheme that dynamically varies the dampening of the oscillating movement of the transducer over time. The control scheme is based upon a predetermined parametric relation involving a plurality of variables derived from the properties measured by the at least one sensor. In several embodiments, the plurality of variables includes a first variable representing excitation frequency of the transducer. In another embodiment, the predetermined parametric relation represents relative phase between two variables derived from the properties measured by the at least one sensor.

## Background
CROSS-REFERENCE TO RELATED APPLICATION
The present application claims the benefit of U.S. Provisional Patent Application Ser. No. 61/587,221 filed on Jan. 17, 2012, which is incorporated herein by reference in its entirety.
FIELD
The present application relates to electromechanical systems for harvesting useful electrical energy from environmental vibrational energy. More specifically, the present application relates to electromagnetic and piezoelectric generators that convert environmental vibrational energy to electrical energy.
BACKGROUND
Harvesting energy from the environmental vibration sources has received a renewed attention in recent years in part due to their applications in wireless networks and portable electronics.
A number of studies in the literature have investigated utilizing discrete electrical components (DuToit et al., 2005; Adhikari et al., 2009; Renno et al., 2009; Cammarano et al., 2010; Trimble, 2011) or linear control systems (Kong et al., 2010; Cassidy et al., 2011; Trimble, 2011).
In another class of studies, switching electronics are utilized to increase harvested power (Guyomar et al., 2005; Lefeuvre et al., 2006; Makihara et al., 2006). These methods were applied to piezoelectric harvesters to generate a favorable electrical load on the active element and increase the voltage output of the harvester which were also shown to increase the power output.
Switching electrical loads to improve transducer performance has been investigated in various vibration control applications (Jalili, 2002; Liu et al., 2005; Wagg and Neild, 2011). Control of shunt resistance (Clark, 2000) and capacitance (Davis and Lesieutre, 2000) electrically connected across a piezoelectric transducer was shown to improve vibration damping performance of these systems.
SUMMARY
This summary is provided to introduce a selection of concepts that are further described below in the detailed description. This summary is not intended to identify key or essential features of the claimed subject matter, nor is it intended to be used as an aid in limiting the scope of the claimed subject matter.
Embodiments are provided for a system that converts environmental vibrational energy into electrical energy for storage in at least one electrical power storage device. The system includes a transducer that undergoes oscillating movement in response to the environmental vibrational energy. The transducer produces an oscillating electrical signal in response to such oscillating movement. Power electronics is operably coupled between the transducer and the at least one electrical power storage device. The power electronics process the oscillating electrical signal produced by the transducer. The power electronics include at least one control element having a configuration that provides variable dampening of the oscillating movement of the transducer in response to at least one control signal supplied thereto. The system further includes at least one sensor for measuring properties of the oscillating movement of the transducer over time. Control electronics is operably coupled to the at least one control element of the power electronics and to the at least one sensor. The control electronics carry out a control scheme that generates and supplies the at least one control signal to the at least one control element over time in a manner that controls the at least one control element to dynamically vary the dampening of the oscillating movement of the transducer over time. The control scheme is based upon a predetermined parametric relation involving a plurality of variables derived from the properties measured by the at least one sensor. In several embodiments, the plurality of variables includes a first variable representing excitation frequency of the transducer. The plurality of variables can also include a second variable representing relative position of the transducer, a third variable representing relative velocity of the transducer, and a fourth variable representing relative acceleration of the transducer.
In one embodiment, the parametric relation includes a term representing a damping factor for controlling the damping of the oscillating movement of the transducer, where the term is equivalent to the sum of a mechanical damper with a constant damping coefficient and a spring with a predetermined spring constant. The predetermined spring constant can be configured such that the system resonance matches a particular off-resonant excitation frequency. The term can be constrained by a lower bound that is greater than or equal to zero such that the damping factor term is positive. Alternatively, the term can be constrained by a lower bound that is less than zero such that the damping factor term can be both negative and positive.
In another embodiment, the parametric relation includes a term representing a damping factor for controlling the damping of the oscillating movement of the transducer, wherein the term is equivalent to a mechanical damper with a constant damping coefficient and a simulated mass. The simulated mass can be configured such that the system resonance matches a particular off-resonant excitation frequency. The term can be constrained by a lower bound that is greater than or equal to zero such that the damping factor term is positive. Alternatively, the term can be constrained by a lower bound that is less than zero such that the damping factor term can be both negative and positive.
In yet another embodiment, the control scheme carried out by the control electronics can be based upon a predetermined parametric relation representing relative phase between two variables derived from the properties measured by the at least one sensor. The two variables can have a phase relationship of 180 degrees out of phase in resonant vibrations, and the control scheme can be configured to bring the two variables into a target phase relationship of 180 degrees out of phase for cases when the phase relationship of the two variables varies from the target phase relationship. The control scheme need not require measurement of the excitation frequency of the transducer over time.
In the various embodiments, the parametric relation of the control scheme can be configured to extend the bandwidth of oscillatory movements of the transducer that produce oscillating electrical signals by the transducer in order to improve the off-frequency energy harvesting performance of the system.
The transducer can include at least one coil and magnet that move relative to one another, wherein the at least one coil produces the oscillating electrical signal. The at least one control element of the power electronics can be realized by a variable load circuit operably coupled to the output of the at least one coil.
Alternatively, the transducer can include at least one piezoelectric element that produces the oscillating electrical signal. The at least one control element of the power electronics can be realized by a variable load circuit operably coupled to the output of the at least one piezoelectric element.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a functional block diagram of an exemplary system that converts environmental vibrational energy into electrical energy for storage in at least one electrical power storage device in accordance with embodiments described herein,
FIG. 2 is a schematic representation of a mechanical model of the electromagnetic vibrational energy harvester of FIG. 1;
FIG. 3 is a schematic representation of an electrical model of the electromagnetic vibrational energy harvester of FIG. 1;
FIG. 4 is a graph of the simulation results showing the power harvested by the system of FIG. 1 as a function of excitation frequency, where the system employs a control scheme having a constant electrical damping factor; the different curves have different constant electrical damping factors as shown;
FIGS. 5(A), 5(B) and 5(C) are graphs of the simulation results showing the transducer excitation, transducer response (relative velocity), and the time varying electrical damping factor of the system of FIG. 1 in the time domain, where the system employs a control scheme referred to herein as a simulated spring control scheme; FIG. 5(A) is a graph of the simulation results of the system for a subcritical excitation. FIG. 5(B) is a graph of the simulation results of the system for resonant excitation. FIG. 5(C) is a graph of the simulation results of the system for a supercritical excitation;
FIGS. 6(A), 6(B) and 6(C) are graphs of the simulation results showing the transducer excitation, transducer response (relative velocity), and the time varying electrical damping factor of the system of FIG. 1 in the time domain, where the system employs a control scheme referred to herein as a simulated mass control scheme; FIG. 6(A) is a graph of the simulation results of the system for a subcritical excitation. FIG. 6(B) is a graph of the simulation results of the system for resonant excitation. FIG. 6(C) is a graph of the simulation results of the system for a supercritical excitation;
FIGS. 7(A), 7(B) and 7(C) are graphs of the simulation results showing the transducer excitation, transducer response (relative velocity), and the time varying electrical damping factor of the system of FIG. 1 in the time domain, where the system employs a control scheme referred to herein as a phase-matching control scheme; FIG. 7(A) is a graph of the simulation results of the system for a subcritical excitation. FIG. 7(B) is a graph of the simulation results of the system for resonant excitation. FIG. 7(C) is a graph of the simulation results of the system for a supercritical excitation;
FIG. 8 is a graph of the simulation results showing the power harvested by the system of FIG. 1 as a function of excitation frequency, where the system employs a number of different semi-active control schemes as described in the present application; the different curves correspond to different control schemes as shown;
FIG. 9 is a graph of the simulation results showing the power harvested by the system of FIG. 1 as a function of excitation frequency, where the system employs a phase-matching control scheme as described in the present application; the different curves correspond to different values for a parameter of the phase-matching control scheme;
FIG. 10 is a graph of the simulation results showing the power harvested by the system of FIG. 1 as a function of modulated excitation frequency, where the system employs a phase-matching control scheme as described in the present application; the different curves correspond to different values for a parameter of the phase-matching control scheme as well as a curve (labeled “Resonant Opt.”) corresponding to control scheme having a constant electrical damping factor;
FIG. 11 is a graph of the simulation results showing the power harvested by the system of FIG. 1 as a function of excitation frequency, where the system employs a number of different bipolar control schemes as described in the present application; the different curves correspond to different control schemes as shown;
FIG. 12 is a schematic diagram of a hydrocarbon producing well according to an embodiment of the present application; and
FIG. 13 is a schematic diagram of a subsurface drilling string according to an embodiment of the present application.



DETAILED DESCRIPTION
FIG. 1 is a schematic diagram of a system in accordance with the present application. The system includes an electromagnetic energy harvester 102 that cooperates with power electronics 104 to convert environmental vibrational energy 106 into electrical energy that is stored in one or more electrical power storage devices 108. For example, the electrical power storage device(s) 108 can be realized by one or more capacitors (or batteries) that store electrical energy and output the stored electrical energy as DC voltage (current). The electrical energy stored by the electrical power storage device(s) 108 can be supplied to one or more electrical loads 110 for powering the electrical loads. Such supply can involve conversion of the electrical energy output from the electrical power storage device(s) 108 into a form suitable for consumption by the electrical loads 110. For example, DC voltage (current) output by the electrical power storage device(s) 108 can be converted into other DC voltage (current) signals (i.e., up-converted or down-converted) and/or possibly converted to AC voltage (current) power signals (such as standard 110V AC or 220V AC power signals).
The electromagnetic energy harvester 102 includes a transducer of inertial mass m attached by a structural spring of a static spring constant k to a vibrating structure secured to a reference frame as shown in the schematic diagram of FIG. 2. For example, the transducer can be realized by one or more coils attached to the inertial mass, and one or more magnets attached to the vibrating structure. Alternatively, the transducer can be realized by one or more magnets attached to the inertial mass, and one or more coils attached to the vibrating structure. Environmental vibration from a vibration source causes the transducer to oscillate and thus causes the coil(s) to cut through and interact with flux formed by the magnet to generate current oscillations in the coil(s). The present application treats this electromagnetic transducer as a variable rate mechanical damper. In order to understand this approximation, one must understand the details of the electromagnetic transduction carried out by the electromagnetic transducer.
More specifically, the electromotive force, V, generated by the electromagnetic transducer is proportional to the speed, {dot over (x)}, of the transducer for a large range of transducers (including the prototypical moving coil transducer) as follows:

V=T{dot over (x)},  (1)

where T is the constant of proportionality of the electromagnetic transducer.
The constant of proportionality T also linearly relates the coil current i to the transducer force FE as follows:

FE=−Ti  (2)

The relationship between the electromotive force, V, and the coil current i is determined by the power electronics 104 that are used in conjunction with the harvester 102.

In the illustrative embodiment, the electrical model of the harvester 102 and the power electronics 104 includes a resistance RC, of the coil(s) of the electromagnetic transducer of the harvester 102 and a variable load of the power electronics 104, which is illustrated by the variable resistance RL as shown in FIG. 3. The variable load resistance RL is a simplified model of rectification, load matching or other more complicated circuitry that is normally utilized for the power electronics 104. Non-limiting examples of variable load circuits included in the power electronics are DC-DC buck converter, buck-boost converter (Kong et al., 2010) and H-bridge switching network (Cassidy et al., 2011). The inductance L of the coil of the electromagnetic transducer, which provides a negligible effect in the low frequency range of the system (where the frequency ω is given as (RC+RL)/L), is also neglected. Using this model, the coil current i can be represented as:




  
    
      
        i
        =
        
          V
          
            
              R
              C
            
            +
            
              R
              L
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




Using Equations (1)-(3), the transducer force can be expressed as:





  
    
      
        
          F
          E
        
        =
        
          
            
              -
              
                
                  T
                  2
                
                
                  
                    R
                    C
                  
                  +
                  
                    R
                    L
                  
                
              
            
            ⁢
            
              x
              .
            
          
          =
          
            
              -
              
                b
                e
              
            
            ⁢
            
              x
              .
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




where be is the electrical damping coefficient.
With this equation, one can see that the coil(s) of the electromagnetic transducer behaves analogous to a mechanical damper with a variable damping rate. The electrical damping coefficient be can be equated to




  
    
      T
      2
    
    
      (
      
        
          R
          C
        
        +
        
          R
          L
        
      
      )
    
  
  ,




and can be controlled by varying the load resistance RL. In this manner, the variable load resistance RL of the power electronics 104 acts as a control element that can be varied to control the damping rate of the electromagnetic transducer. The term electrical damping is utilized in this application to refer to electromechanical transduction. Therefore the mechanical energy is converted into the electrical domain in the process and not dissipated.

For a pure AC, zero-inductance power source, the transfer of power to the load is maximized if the load resistance is matched to the coil resistance. Variability in the input impedance of the electrical loads 110 that will be powered by the harvester 102 and associated power storage devices 108 uses control electronics 112 that achieve such load matching. Applicability of such a simple load matching control strategy is limited since this approach results in optimal power harvesting in devices which output a pure, constant-frequency sinusoid. Another implicit assumption is that the harvester force has no effect on the mechanical system. Considering that vibrational harvesters operate near their resonance, where the stiffness forces are balanced with inertial forces, a small amount of change in damping can result in a pronounced difference in performance. It is thus important to analyze the mechanical system of the energy harvester 102 to determine the optimal control strategy carried out by the control electronics 112.
In the illustrative embodiment, the energy harvester 102 and associated power electronics 104 are modeled by a mass-spring resonator as shown in FIG. 2. The mass-spring resonator has two dampers:

    
    
        a fixed damper with a damping coefficient bm to model the intrinsic mechanical damping of the system; and
        a variable damper with a damping coefficient be to model the electromagnetic transducer.
    
    


The control electronics 112 of the energy harvester 102 can be configured to carry out a control scheme that selects an appropriate value of the damping coefficient be to enhance the system bandwidth. The control scheme can be a simulated spring control scheme, a simulated mass control scheme or a phase-matching control scheme as summarized below.
In the simulated spring control scheme, the damping coefficient be is varied such that the electromagnetic transducer behaves similar to a spring and a fixed-rate damper whenever possible. This is particularly useful if the background vibration is near-sinusoidal with a background excitation frequency not matching the resonant frequency. The spring term for this control scheme is adjusted to match what is required to match these two frequencies. This results in a damping coefficient be of the from:




  
    
      
        
          b
          e
        
        =
        
          
            b
            opt
          
          +
          
            
              (
              
                
                  
                    ω
                    2
                  
                  ⁢
                  m
                
                -
                k
              
              )
            
            ⁢
            
              y
              
                y
                .
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




where ω is the background excitation frequency, bopt is the fixed rate optimal damper simulated by the electromagnetic transducer, and





  
    (
    
      
        
          ω
          2
        
        ⁢
        m
      
      -
      k
    
    )
  
  ⁢
  
    y
    
      y
      .
    
  




is the spring term that is adjusted by the control scheme.

The simulated spring term [i.e.,




  
    
      (
      
        
          
            ω
            2
          
          ⁢
          m
        
        -
        k
      
      )
    
    ⁢
    
      y
      
        y
        .
      
    
  
  ]




of Eqn. 5 is based on the displacement y, the relative speed {dot over (y)}, the mass m, the spring constant k, and the base excitation frequency ω. The displacement y and the relative speed {dot over (y)} can be derived from the output of the sensor(s) 114. The base excitation frequency ω can be derived from analog or digital processing of a sensor signal. The parameter m is dictated by the inertial mass of the transducer, and the parameter k is dictated by the spring constant of the transducer. The damping coefficient be can be limited by predetermined bounds for both semi-active (uni-polar) schemes where the electrical damping is limited to positive values (be>0) as well as bi-polar schemes.

In the simulated mass control scheme, the damping coefficient be is varied such that the electromagnetic transducer behaves with a simulated mass that is adjusted such that resonance frequency of the system matches the background excitation frequency. The damping coefficient be in this case is:




  
    
      
        
          b
          e
        
        =
        
          
            b
            opt
          
          +
          
            
              k
              ⁡
              
                (
                
                  
                    1
                    
                      ω
                      2
                    
                  
                  -
                  
                    1
                    
                      ω
                      0
                      2
                    
                  
                
                )
              
            
            ⁢
            
              
                y
                 ̈
              
              
                y
                .
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




where bopt is the fixed rate optimal damper simulated by the electromagnetic transducer, and





  
    k
    ⁡
    
      (
      
        
          1
          
            ω
            2
          
        
        -
        
          1
          
            ω
            0
            2
          
        
      
      )
    
  
  ⁢
  
    
      y
       ̈
    
    
      y
      .
    
  




is the mass term that is adjusted by the control scheme.

The simulated mass term [i.e.,





  
    
      k
      ⁡
      
        (
        
          
            1
            
              ω
              2
            
          
          -
          
            1
            
              ω
              0
              2
            
          
        
        )
      
    
    ⁢
    
      
        y
         ̈
      
      
        y
        .
      
    
  
  ]




of Eqn. 6 is based on the relative speed {dot over (y)}, the relative acceleration ÿ, spring constant k, the base excitation frequency ω, and the resonance frequency ω0. The relative speed {dot over (y)} and the relative acceleration ÿ can be derived from the output of the sensor(s) 114. The parameter k is dictated by the spring constant of the transducer. The parameter ω0 is the natural frequency of the transducer given by ω0=√{square root over (k/m)} where m is the inertial mass of the transducer. The base excitation frequency ω can be derived from analog or digital processing of a sensor signal. The damping coefficient be can be limited by predetermined bounds for both semi-active (uni-polar) schemes where the electrical damping is limited to positive values (be>0) as well as bi-polar schemes.

In resonant vibration of the system, there are fixed phase relationships that exist between system variables. For example the excitation acceleration {umlaut over (x)} is 180 degrees out-of-phase with the relative velocity {dot over (y)}. For non-resonant excitations however, these phase relationships vary from the resonant case. In the phase-matching control scheme these relationships are used to adjust the damping coefficient be to move a non-resonance vibration excitation towards resonance. One such control law is of the form:




  
    
      
        
          b
          e
        
        =
        
          {
          
            
              
                
                  b
                  e
                  
                    m
                    ⁢
                    
                      
                    
                    ⁢
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                
              
              
                
                  
                    
                      x
                       ̈
                    
                    ⁢
                    
                      y
                      .
                    
                  
                  ≤
                  0
                
              
            
            
              
                
                  b
                  e
                  
                    m
                    ⁢
                    
                      
                    
                    ⁢
                    a
                    ⁢
                    
                      
                    
                    ⁢
                    x
                  
                
              
              
                
                  
                    
                      x
                       ̈
                    
                    ⁢
                    
                      y
                      .
                    
                  
                  >
                  0
                
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




where bemin and bemax are the lower and upper bounds on the damping coefficient be.

The relative speed {dot over (y)} and the excitation acceleration {umlaut over (x)} can be derived from the output of the sensor(s) 114. The bounds bemin and bemax of Eqn. (7) can be dictated by the capabilities of the electromagnetic harvester or imposed by the designer to improve system performance for both semi-active (uni-polar) schemes where the electrical damping is limited to positive values (be>0) as well as bi-polar schemes.

In dimensionless form, the equation of motion of the transducer of FIG. 2 is given as:

ÿ+2(ζe+ζm){dot over (y)}+y=−{umlaut over (x)},  (8)

with a base excitation of {umlaut over (x)}=sin(ωt).
For brevity all variables in the remainder of this application will be switched with their dimensionless form. Their dimensional forms will be denoted with tilde ({tilde over (●)}). The damping factor ζe=be/2ω0m and ωm=bm/2ω0m can be derived from the natural frequency ω0 of the system given by ω0=√{square root over (k/m)}. The dimensionless quantities in the equations can be defined in terms of their dimensional counterparts as x={tilde over (x)}ω02/a0, y={tilde over (y)}ω02, t={tilde over (t)}ω0, and ω={tilde over (ω)}/ω0, where a0 is the input acceleration amplitude.
Given this model, the control electronics 112 is configured to select an appropriate value of the damping coefficient be (and thus the corresponding electrical damping factor ζe) in order to maximize the output power. For resonant harvesting where ω=1, this is accomplished by setting be equal to bm (such that ζe=ζm). However, when the excitation source is not resonant (i.e., ω≠1), this is no longer the case.
Power harvested in constant electrical damping systems as a function of excitation frequency is presented in FIG. 4 with ζm=0.01. The values in this plot are calculated using a time-difference numerical simulation in the time domain. The system is excited until the initial disturbances are damped out. A time-average of the electrically damped power, P=2ζe{dot over (y)}2, is calculated over a period much longer than the natural and excitation period of the system. At resonance where ω=1, the power output to the electromagnetic transducer is maximized by setting be equal to bm (such that ζe=ζm). The bandwidth of the system can be increased by increasing the electrical damping factor; however, this comes in the cost of decreasing peak (resonant) power. The power output from the generalized optimal damping factor ζGO is the most power that can be converted to electrical energy at a given input frequency.
FIG. 4 demonstrates the traditional trade-off made in constant electrical damping systems. The resonant power harvested is maximized when the electrical damping is matched with the mechanical damping; however, if the system is not properly tuned to the excitation frequency, the off-resonant power output drops off rapidly. The controller 112 can address this concern by increasing electrical damping and thus the bandwidth of the system; however, this also reduces the resonant peak. Therefore the trade-off that is often required in these systems is between the bandwidth and peek resonant power.
Another way to think about this is that, at an off-resonant input frequency, the optimal damping coefficient is higher than that at resonance. This optimal damping factor ζGO can be analytically calculated in the frequency domain based upon the complex relative speed of the harvester 102 as follows:




  
    
      
        
          
            Y
            .
          
          ⁡
          
            (
            ω
            )
          
        
        =
        
          
            
              iω
              3
            
            
              
                (
                
                  1
                  -
                  
                    ω
                    2
                  
                
                )
              
              +
              
                2
                ⁢
                
                  iω
                  ⁡
                  
                    (
                    
                      
                        ζ
                        e
                      
                      +
                      
                        ζ
                        m
                      
                    
                    )
                  
                
              
            
          
          .
        
      
    
    
      
        (
        9
        )
      
    
  




Using the complex amplitude theorem (Eq. 2.15.14 of Melcher J R (1981) Continuum Electromechanics. MIT Press), the average power harvested can be given as:

P=1⁄2(2ζe){dot over (Y)}(ω){dot over (Y)}*(ω),  (10)

where {dot over (Y)}*(ω) is the complex conjugate of relative speed.
As with other variables, the power is in dimensionless form and can be related to its dimensional form by P={tilde over (P)}ω0/ma02. The generalized optimal damping factor ζGO is determined by setting dP/dζe=0, and picking the positive root (as the negative root is not physically significant) as follows:




  
    
      
        
          ζ
          GO
        
        =
        
          
            
              (
              
                
                  
                    (
                    
                      
                        1
                        -
                        
                          ω
                          2
                        
                      
                      
                        2
                        ⁢
                        ω
                      
                    
                    )
                  
                  2
                
                +
                
                  ζ
                  m
                  2
                
              
              )
            
            
              1
              /
              2
            
          
          .
        
      
    
    
      
        (
        11
        )
      
    
  




The electrical damping factor ζe can then be set to the generalized optimal damping factor ζGO. As illustrated in FIG. 4, the generalized optimal damping ζGO broadens the power peak without sacrificing the resonant performance. However, in order to implement this control strategy, partial state feedback is necessary to determine the excitation frequency ω.

In an illustrate embodiment, the control electronics 112 is configured to vary the variable load resistance RL (control element) of the power electronics 104 in order to control the electrical damping factor ζe (or the corresponding electrical damping coefficient be) of the electromagnetic transducer in a manner that improves the off-resonant performance of the energy harvester 102. The control electronics 112 utilizes one or more sensors 114 that measure the base excitation x of the vibrating structure and the displacement y of the transducer relative to the vibrating structure and provide such measurements as inputs to the control electronics 112. The sensors 114 can also possibly measure the relative speed, {dot over (y)}, through the voltage output of the electromagnetic transducer.
For certain applications, the control electronics 112 can operate in a semi-active mode where the electrical damping is limited to positive values (ζe>0). In this semi-active mode, the electrical damping factor ζe can never be less than zero. As described above, the harvester 102 provides its peak power output at resonant frequency since the mechanical motion is best amplified at this critical frequency.
In one control scheme suitable for the semi-active mode, which is referred to herein as the simulated spring control scheme, the load resistances RL is modulated to change the electrical damping factor ζe (or the corresponding electrical damping coefficient be) such that the damping factor ζe satisfies the following parametric relations:




  
    
      
        
          
            ζ
            e
          
          =
          
            max
            ⁡
            
              [
              
                
                  ζ
                  e
                  min
                
                ,
                
                  
                    ζ
                    e
                    sim
                  
                  +
                  
                    
                      y
                      
                        2
                        ⁢
                        
                          y
                          .
                        
                      
                    
                    ⁢
                    
                      (
                      
                        
                          ω
                          2
                        
                        -
                        1
                      
                      )
                    
                  
                
              
              ]
            
          
        
        ⁢
        
          

        
        ⁢
        
          for
          ⁢
          
            
          
          ⁢
          the
          ⁢
          
            
          
          ⁢
          case
          ⁢
          
            
          
          ⁢
          where
        
        ⁢
        
          

        
        ⁢
        
          
            
              max
              ⁡
              
                [
                
                  
                    ζ
                    e
                    min
                  
                  ,
                  
                    
                      ζ
                      e
                      sim
                    
                    +
                    
                      
                        y
                        
                          2
                          ⁢
                          
                            y
                            .
                          
                        
                      
                      ⁢
                      
                        (
                        
                          
                            ω
                            2
                          
                          -
                          1
                        
                        )
                      
                    
                  
                
                ]
              
            
            <
            
              ζ
              e
              max
            
          
          ;
        
        ⁢
        
          

        
        ⁢
        or
      
    
    
      
        (
        
          12
          ⁢
          A
        
        )
      
    
  
  
    
      
        
          
            ζ
            e
          
          =
          
            ζ
            e
            max
          
        
        ⁢
        
          

        
        ⁢
        
          for
          ⁢
          
            
          
          ⁢
          the
          ⁢
          
            
          
          ⁢
          case
          ⁢
          
            
          
          ⁢
          where
        
        ⁢
        
          

        
        ⁢
        
          
            max
            ⁡
            
              [
              
                
                  ζ
                  e
                  min
                
                ,
                
                  
                    ζ
                    e
                    sim
                  
                  +
                  
                    
                      y
                      
                        2
                        ⁢
                        
                          y
                          .
                        
                      
                    
                    ⁢
                    
                      (
                      
                        
                          ω
                          2
                        
                        -
                        1
                      
                      )
                    
                  
                
              
              ]
            
          
          ≥
          
            
              ζ
              e
              max
            
            .
          
        
      
    
    
      
        (
        
          12
          ⁢
          B
        
        )
      
    
  




The parametric relation of Eqn. (12A) sets the electrical damping factor ζe to the maximum of two values for the case where such maximum is less than a predetermined upper bound ζemax. Otherwise, the electrical damping factor ζe is set to predetermined upper bound ζemax as provided by the parametric relation of Eqn. (12B). The upper bound ζemax may be dictated by the limitations of the electromagnetic transducer.
The first value of the parametric relation of Eqn. (12A) represents a predetermined lower bound ζemin. For example, ζemin can be set to 0 to ensure that the mechanical system does not receive power from the electrical system. The second value of the parametric relation of Eqn. (12A) represents the output equivalent to the sum of the outputs of a mechanical damper with a constant damping coefficient ζesim and a spring with a predetermined spring constant. As adding a spring to the harvester 102 would shift the resonant frequency, this strategy seeks to take an off-resonant excitation frequency and add a spring to the system such that the system resonance matches this frequency. Note that the simulated spring term [i.e.,




  
    
      y
      
        2
        ⁢
        
          y
          .
        
      
    
    ⁢
    
      (
      
        
          ω
          2
        
        -
        1
      
      )
    
  
  ]




of Eqn. (12A) can be negative, and thus the strategy is applicable to both subcritical and supercritical excitations. The simulated spring term [i.e.,





  
    
      y
      
        2
        ⁢
        
          y
          .
        
      
    
    ⁢
    
      (
      
        
          ω
          2
        
        -
        1
      
      )
    
  
  ]




of the parametric relation of Eqn. (12A) is based on the displacement y, the relative speed {dot over (y)} and the base excitation frequency ω. The displacement y and the relative speed {dot over (y)} are derived from the output of the sensor(s) 114. The base excitation frequency ω can be derived from analog or digital processing of a sensor signal. At steady-state all system states (e.g. position, current) oscillate with a frequency equivalent to the excitation frequency. Therefore, by processing the time history of a sensor output it is possible to determine the excitation frequency. A common algorithm utilized for determining the frequency content of a signal is the Fourier transform.

In the simulated spring control scheme, for critical excitations where ω=1, the electrical damping coefficient is constant and equal to the mechanical damping coefficient (ζe=ζm).
The behavior of the simulated spring control scheme is shown in time domain simulations results of FIGS. 5(A), 5(B) and 5(C). Note the simulated spring control scheme behaves identical to an optimally damped harvester in the resonant case where ω=1 as shown in 5(B). The electrical damping coefficient is constant and equal to the mechanical damping coefficient (ζe=ζm). For subcritical excitations where ω<1, the simulated spring control scheme results in a high damping factor while the relative speed {dot over (y)} is leading excitation position x as shown in FIG. 5(A). For supercritical excitations where ω>1, the simulated spring control scheme generates a high damping factor at times when the relative speed is lagging. For both the subcritical and supercritical excitations (where ω≠1), the system seeks to simulate the required positive or negative spring constant. The system achieves this when the signs of the displacement y and the relative speed {dot over (y)} allow for a simulated spring load and when the transducer is capable of applying the load.
In the simulated spring control scheme, the damping of the electromagnetic transducer behaves similar to a spring and a fixed-rate damper whenever possible. This is particularly useful in maximizing the harvested power in the event that the background vibration is near-sinusoidal with a frequency not matching the resonant frequency. The added spring is adjusted to match what is required such that resonant frequency of the system matches the excitation frequency.
In another control scheme suitable for the semi-active mode, which is referred to herein as the simulated mass control scheme, the load resistance RL is modulated to change the electrical damping factor ζe (or the corresponding electrical damping coefficient be) such that the damping factor ζe satisfies the following parametric relations:




  
    
      
        
          
            
              ζ
              e
            
            =
            
              max
              ⁡
              
                [
                
                  
                    ζ
                    e
                    min
                  
                  ,
                  
                    
                      ζ
                      e
                      sim
                    
                    +
                    
                      
                        
                          y
                           ̈
                        
                        
                          2
                          ⁢
                          
                            y
                            .
                          
                        
                      
                      ⁢
                      
                        (
                        
                          
                            1
                            
                              ω
                              2
                            
                          
                          -
                          1
                        
                        )
                      
                    
                  
                
                ]
              
            
          
          ,
          
            

          
          ⁢
          
            for
            ⁢
            
              
            
            ⁢
            the
            ⁢
            
              
            
            ⁢
            case
            ⁢
            
              
            
            ⁢
            where
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              max
              ⁡
              
                [
                
                  
                    ζ
                    e
                    min
                  
                  ,
                  
                    
                      ζ
                      e
                      sim
                    
                    +
                    
                      
                        
                          y
                           ̈
                        
                        
                          2
                          ⁢
                          
                            y
                            .
                          
                        
                      
                      ⁢
                      
                        (
                        
                          
                            1
                            
                              ω
                              2
                            
                          
                          -
                          1
                        
                        )
                      
                    
                  
                
                ]
              
            
            <
            
              ζ
              e
              max
            
          
          ;
        
        ⁢
        
          

        
        ⁢
        or
      
    
    
      
        (
        
          13
          ⁢
          A
        
        )
      
    
  
  
    
      
        
          
            
              ζ
              e
            
            =
            
              ζ
              e
              max
            
          
          ,
          
            

          
          ⁢
          
            for
            ⁢
            
              
            
            ⁢
            the
            ⁢
            
              
            
            ⁢
            case
            ⁢
            
              
            
            ⁢
            where
          
        
        ⁢
        
          

        
        ⁢
        
          
            max
            ⁡
            
              [
              
                
                  ζ
                  e
                  min
                
                ,
                
                  
                    ζ
                    e
                    sim
                  
                  +
                  
                    
                      
                        y
                         ̈
                      
                      
                        2
                        ⁢
                        
                          y
                          .
                        
                      
                    
                    ⁢
                    
                      (
                      
                        
                          1
                          
                            ω
                            2
                          
                        
                        -
                        1
                      
                      )
                    
                  
                
              
              ]
            
          
          ≥
          
            
              ζ
              e
              max
            
            .
          
        
      
    
    
      
        (
        
          13
          ⁢
          B
        
        )
      
    
  




The parametric relation of Eqn. (13A) sets the electrical damping factor ζe to the maximum of two values for the case where such maximum is less than a predetermined upper bound ζemax. Otherwise, the electrical damping factor ζe is set to predetermined upper bound ζemax according to the parametric relation of Eqn. (13B). The upper bound ζemax may be dictated by the limitations of the electromagnetic transducer.
The first value of the parametric relation of Eqn. (13A) represents a predetermined lower bound ζemin. For example, ζemin can be set to 0 to ensure that the mechanical system does not receive power from the electrical system. The second value of the parametric relation of Eqn. (13A) seeks to simulate a system with a resonance matching the excitation frequency by modifying the mass of the transducer. The simulated mass term [i.e.,




  
    
      
        y
         ̈
      
      
        2
        ⁢
        
          y
          .
        
      
    
    ⁢
    
      (
      
        
          1
          
            ω
            2
          
        
        -
        1
      
      )
    
  
  ]




of the parametric relation of Eqn. (13A) is based on the relative speed {dot over (y)}, the relative acceleration ÿ and the base excitation frequency ω. The relative speed {dot over (y)} and the relative acceleration ÿ are derived from the output of the sensor(s) 114. The base excitation frequency ω can be derived from analog or digital processing of a sensor signal. At steady-state all system states (e.g. position, current) oscillate with a frequency equivalent to the excitation frequency. Therefore, by processing the time history of a sensor output it is possible to determine the excitation frequency. A common algorithm utilized for determining the frequency content of a signal is the Fourier transform.

In the simulated mass control scheme, for critical excitations where ω=1, the electrical damping coefficient is constant and equal to the mechanical damping coefficient (ζe=ζm).
The behavior of the simulated mass control scheme is shown in time domain simulations results of FIGS. 6(A), 6(B) and 6(C). Note the simulated mass control scheme behaves identical to an optimally damped harvester in the resonant case where ω=1 as shown in 6(B). The electrical damping coefficient is constant and equal to the mechanical damping coefficient (ζe=ζm). FIG. 6(A) shows the simulation results of the simulated mass control scheme for subcritical excitations where ω<1. FIG. 6(C) shows the simulation results of the simulated mass control scheme for supercritical excitations where ω>1. The simulation results are similar to the results of the simulated spring control scheme. One difference is the oscillating nature of the damping factor for subcritical excitations as shown in FIG. 6(A). This is mainly caused by the fact that the control system is generating a negative simulated mass which leads to instability in the control function.
In the simulated mass control scheme, for off-resonance excitation, the damping of the electromagnetic transducer is controlled by a simulated mass that is adjusted to match what is required such that resonant frequency of the system matches the off-resonant excitation frequency. This is particularly useful in maximizing the harvested power from the off-resonant vibrational energy.
In resonant vibrations the excitation acceleration {umlaut over (x)} and the relative velocity {dot over (y)} are 180° out-of-phase. In yet another control scheme suitable for the semi-active mode, which is referred herein to as the phase-matching control scheme, the load resistance RL is modulated to change the electrical damping factor ζe (and the corresponding electrical damping coefficient be) such that the excitation acceleration {umlaut over (x)} and the relative velocity {dot over (y)} are brought into the 180° out-of-phase relationship for the cases in which they are not. This can be accomplished by the following parametric relations:

ζe=ζemax for the case that {umlaut over (x)}{dot over (y)}≦0, or  (14A)

ζe=ζemax for the case that {umlaut over (x)}{dot over (y)}>0.  (14B)

ζemax represents a predetermined upper bound, which may be dictated by the limitations of the electromagnetic transducer. ζemin represents a predetermined lower bound, which can be set to ζm in a typical case. The parametric relations of Eqns. (14A) and (14B) are based on the sign of the product of the excitation acceleration {umlaut over (x)} and the relative velocity {dot over (y)}, which provides a measure of relative phase offset between the excitation acceleration {umlaut over (x)} and the relative velocity {dot over (y)}. The excitation acceleration {umlaut over (x)} and the relative velocity {dot over (y)} can be derived from the output of the sensor(s) 114. In an alternate embodiment, the control system can utilize the phase relationship between the excitation position x and relative velocity {dot over (y)} since these variables are in-phase during resonant excitations. Any deviation from this phase relationship caused by an off-resonant excitation can be detected by a parametric relation similar to the one given in Eqns. (14A) and (14B).

Although the phase-matching control scheme appears to be a more indirect approach, it is important to note that this phase-matching control scheme does not require the determination of excitation frequency ω, and thus may be implemented with less electronic complexity. It is also possible to generate similar control actions utilizing resonant phase relationships between other variables (such as displacement variables, velocity variables, acceleration variables, current variables and voltage variables of the transducer as well as current variables and voltage variables of the variable load circuitry of the power electronics 104) as measured by the output of the sensor(s) 114 of the system.
The behavior of the phase-matching control scheme is shown in time domain simulations results of FIGS. 7(A), 7(B) and 7(C). Note the phase-matching control scheme behaves identical to an optimally damped harvester in the resonant case where ω=1 as shown in 7(B). FIG. 7(A) shows the simulation results of the phase-matching control scheme for subcritical excitations where ω<1. FIG. 7(C) shows the simulation results of the phase-matching control scheme for supercritical excitations where ω>1. The simulation results are similar to the results of the simulated spring control scheme.
Alternate phase-matching control schemes can employ another variable pair that has a phase relationship of in-phase or 180 degrees out of phase in resonant vibrations. In such configuration, the load resistance RL can be modulated to change the electrical damping factor ζe (and the corresponding electrical damping coefficient be) to bring the variable pair into a target phase relationship of in-phase or 180 degrees out of phase for cases, respectively, when the phase relationship of the variable pair varies from the target phase relationship.
FIG. 8 shows the average power harvested by the control schemes described above for comparison of the power generation capabilities of these schemes. FIG. 8 also shows the power generated with the generalized optimal damping factor of Eqn. (8). Since the generalized optimal damping factor provides the maximum average power output that is possible at a given excitation frequency with a constant damping factor, it serves as a good baseline comparison to the other control strategies. The plots show the average power harvested under harmonic loading conditions. In order to provide a direct comparison, the electronic damping factor is bounded to ζemax=10 for all of the generated plots. The lower bound ζemin is set to zero for all of the control schemes except for the phase-matching control, for which ζemin=ζm=0.01. The generalized optimal electrical damping provides a larger bandwidth than a system optimized only for the resonant peak. All of the other control schemes utilizing the time-varying damping factor provide a broader peak without sacrificing resonant performance, with phase-matching approach providing the best overall response characteristics.
In the simulated mass and spring control schemes, the lower bound ζemin of the electrical damping factor is not an important parameter as long as it can be set to be as low as the mechanical damping factor. However, in the phase-matching control scheme, the lower bound ζemin of the electrical damping factor can affect the harvesting performance as demonstrated in FIG. 9. The system performance can be improved slightly for off-resonant excitations with smaller damping factors; however, this results in a penalty in power harvested near resonance. In one embodiment, the parameter ζemin can be set to ζm in order to give overall performance as evident from FIG. 9.
An additional advantage of the phase-matching control scheme is that it allows the system to seek resonant operation quicker under a time-varying harmonic excitation. In order to demonstrate this, the system was simulated with a phase-modulated input signal as follow:

{dot over (x)}=sin(ωt+π sin(ωpt)),  (15)

where ωp is the modulation frequency of the base excitation.
The power output of three phase-matching control schemes of varying ζemin for ζm=0.01 under resonant excitations (ω=1) of varying modulation frequency ωp is plotted in FIG. 10. The power output for the optimal constant damping factor is also plotted. All of the systems result in the same power output for low levels of the modulation frequency ωp since under this small perturbation from a resonant excitation, the system response gets minimally affected. However, as ωp increases, the phase-matching control output becomes higher as this system allows for the inertial mass to seek the resonant output quicker.
For other applications, the control electronics 112 can operate in a bipolar mode where the electrical damping can have negative values. This bipolar mode allows the system to convert electrical energy back into the mechanical domain to enhance harvesting and thus requires capacitive or inductive energy stored on the electronics. For this bipolar mode, the simulated spring control scheme, the simulated mass control scheme and the phase-matching control scheme as described above can be adapted by allowing the lower bound ζemin of the electrical damping factor to be set to negative values. FIG. 11 is a plot of exemplary simulation results from these control schemes where the lower bound ζemin of the electrical damping factor is set to −10. The average power output of the mechanical system is shown over an interval much longer than the excitation frequency (20 cycles), thus the parameter observed <P> is net harvested energy in FIG. 11. Note that the simulated spring control scheme appears to be more effective than the simulated mass control for subcritical excitations. For supercritical excitations, the opposite is the case. Compared with results obtained with ζemin=0, the peak is broader and thus allows for a more robust harvester to be implemented. Simulated mass control shows more drawbacks for excitations far from resonant frequency.
Advantageously, the various control schemes for the electromagnetic energy harvesting system as described herein show potential in maximizing power output. Specifically, the control schemes enhance the off-resonant peak power, effectively increasing the harvester bandwidth without reducing resonant performance. Such benefits can reduce the constraints for system tuning or in effect realize a more robust harvester system.
It is also contemplated that other circuit models of the harvester 102 and power electronics 104 can be used. In this case, one or more different control elements can be controlled by the control electronics 112 in order to control the damping rate of the coil(s) of the electromagnetic transducer of the harvester 102 for broadband response.
It is also contemplated that the various control schemes as described herein can be applied to vibrational energy harvesters that utilize a piezoelectric transducer. Although the piezoelectric materials provide a non-analogous electromagnetic transduction mechanism to the electromagnetic transducers described above, their force profiles are governed by the electrical load that is placed on them, and hence the harvesting electronics can be utilized to implement the various control schemes as described herein. For example, the phase-matching control system may be implemented to vary the resistance of a shunt resistor Rs connected across the piezoelectric element as follows:

Rs=R1 for the case that {umlaut over (x)}{dot over (y)}≦0, or  (16A)

Rs=R2 for the case that {umlaut over (x)}{dot over (y)}>0.  (16B)

where R1 and R2 are different shunt resistance levels.
The off-resonant response of the energy harvester system as described herein is particularly important in their performance under broadband vibrations. The ability of the energy harvesting system to either adapt to an non-resonant input (tuning) or have a broad frequency response to accept a range of inputs allow for improved efficiency in cases where a non-harmonic background vibrations are present. Utilizing adaptable electronics to improve harvesting efficiency has advantages as there is no added complexity to the mechanical system.
As discussed earlier, a certain amount of electronics are commonly utilized with vibrational energy harvesters to provide a desired load on the system to increase the output power. These electronics also often utilize feedback control schemes that operate at switching frequencies far exceeding the harvester mechanical resonance. It is thus possible to implement a wide range of control schemes without adding complexity to existing systems.
The technological advances of the embodiments described above can be utilized in subsurface tools for oil & gas exploration and production as these embodiments can aid the designers of such subsurface tools to implement simple methods of collocating devices and power sources that derive electrical energy from harvesting of environmental vibrational energy.
FIG. 12 shows an embodiment of a subsurface tool 201 that includes a tubular string 203 (a production string, for example) that extends into a wellbore 205. The tubular string 203 may include a central passageway that communicates a flow 209 from a subterranean formation zone 211 (or to a formation zone in the case of an injection well). The zone 211 represents one out of many possible zones traversed by the wellbore 205. The zone 211 may be defined (i.e., isolated from other zones) by one or more packers 213 (one being depicted in FIG. 12).
The flow 209 is a primary source of vibrational energy downhole, and this vibrational energy can be captured by a vibrational energy harvesting system 215 housed in a power generation tool body 217 for purposes of converting the vibrational energy into downhole electrical power. This electrical power, in turn, may be used to power one or more downhole power-consuming components, such as sleeve valves, ball valves, motors, actuators, sensors, sound sources, electromagnetic signaling sources, or equipment to fire “smart bullets” into a well casing, perforating gun firing heads, controllers, microprocessors, Micro Electrical Mechanical Sensors (MEMS), telemetry systems (transmitters or receivers), etc., depending on the particular application. The vibrational harvesting system 215 of the power generation tool body 217 can be based on the embodiments described above.
The tubular string 203 can include one or more vibration enhancement mechanisms (one shown as 219) that enhance the generation of vibrational energy for capture by the vibrational energy harvesting system 215. For example, the flow 209 can enter the vibration enhancement mechanism 219, which operates to produce a locally more turbulent flow that flows uphole. The creation of this more turbulent flow, in turn, amplifies the vibrational energy, thereby leading to the increased production of downhole power. The vibrational harvesting system 215 can be located in proximity to (within ten feet, for example) to the vibration enhancing mechanism 219.
The vibration enhancing mechanism 219 can include a structural member (such as a flow path obstructer or diverter or Venturi-type orifice or flexible moving member) that is configured to create turbulence in the flow 209 in order to generate environmental vibrations downhole in response to the flow 209. The vibration enhancing mechanism 219 can also be realized by varying the wall thickness of a portion of the tubular string 203 in order to enhance the fundamental vibration mode of the tubular string portion. Other suitable vibration enhancing mechanisms can also be used.
In other embodiments, vibrational energy harvesting systems as described herein can be integrated as part of the subsurface production tool to harvest vibrational energy generated by artificial lift devices (such as gas lift valves and electrical submersible pumps). In yet other embodiments, the vibrational energy harvesting systems as described herein can be integrated as part of the surface-located production tools as well as subsea-located tree structures that monitor and control the operation of subsea wells in order to harvest vibration from the flow of produced fluid.
FIG. 13 depicts a drill string 301 that includes a mud motor 303 and a drill bit 305. The drill string 301 may include a tool body 307 with sensors that are used for purposes of monitoring operation of the drill string 301 and monitoring general operation of the drilling. The sensors are located close to the drill bit 305. A particular challenge with this arrangement is that the sensor tool body 307 may be located away from a power source and thus, electrical cables may have to span across the mud motor 303 for purposes of delivering power to the sensors of the tool body 307. However, in an embodiment of the present application, the sensor tool body 307 may be in close proximity to a vibrational energy harvesting system 309 housed in a power generation tool body 311 for purposes of converting downhole vibrational energy of the drill string 301 (which occurs during operation of the mud motor 303 and interaction between the drill bit 305 and the formation rock). This electrical power, in turn, may be used to power one or more downhole power-consuming components, such as the sensors, telemetry system (transmitter and receiver) and/or other electrical components of the sensor tool body 307. The vibrational harvesting system 309 of the power generation tool body 311 can be based on the embodiments described above. Due to this arrangement, cabling does not have to be extended across the mud motor 303 for purposes of delivering power to the electrical components of the sensor tool body 307. Furthermore, due to the generation of electrical power downhole, various control lines and electrical cables do not need to be extended from the surface of the well.
There have been described and illustrated herein several embodiments of a system that converts environmental vibrational energy into electrical energy for storage in at least one electrical power storage device and control schemes and methodologies employed therein. While particular embodiments of the invention have been described, it is not intended that the invention be limited thereto, as it is intended that the invention be as broad in scope as the art will allow and that the specification be read likewise. Thus, while particular mechanical models and electrical models have been disclosed for a vibrational energy harvester, it will be appreciated that other mechanical models and electrical models can be used as well. In addition, while particular variables and constraints have been disclosed for the control schemes described herein, it will be understood that the control schemes can be adapted to utilize different variables and constraints. Moreover, while particular system configurations have been disclosed, it will be appreciated that other system configurations could be used as well. It will therefore be appreciated by those skilled in the art that yet other modifications could be made to the provided invention without deviating from its spirit and scope as claimed.

## Claims

1. A system for converting environmental vibrational energy into electrical energy for storage in at least one electrical power storage device, the system comprising:
a transducer that undergoes oscillating movement in response to the environmental vibrational energy, the transducer producing an oscillating electrical signal in response to the oscillating movement;
power electronics, operably coupled between the transducer and the at least one electrical power storage device, the power electronics processing the oscillating electrical signal produced by the transducer, wherein the power electronics includes at least one control element having a configuration that provides variable dampening of the oscillating movement of the transducer in response to at least one control signal supplied thereto;
at least one sensor for measuring properties of the oscillations of the transducer over time; and
control electronics operably coupled to the at least one control element of the power electronics and to the at least one sensor, the control electronics carrying out a control scheme that generates and supplies the at least one control signal to the at least one control element over time in a manner that controls the at least one control element to dynamically vary the dampening of the oscillating movement of the transducer over time, wherein the control scheme is based upon a predetermined parametric relation involving a plurality of variables derived from the properties measured by the at least one sensor, wherein the plurality of variables include a first variable representing excitation frequency of the transducer.

2. A system according to claim 1, wherein:
the plurality of variables further includes a second variable representing relative position of the transducer.

3. A system according to claim 1, wherein:
the plurality of variables further includes a third variable representing relative velocity of the transducer.

4. A system according to claim 1, wherein:
the plurality of variables further includes a fourth variable representing relative acceleration of the transducer.

5. A system according to claim 1, wherein:
the plurality of variables further includes a fifth variable representing electrical current of the transducer.

6. A system according to claim 1, wherein:
the plurality of variables further includes a sixth variable representing electrical voltage of the transducer.

7. A system according to claim 1, wherein:
the predetermined parametric relation includes a term representing a damping factor for controlling the damping of the oscillating movement of the transducer, which results in the transducer force to be equivalent to the sum of a spring with a predetermined spring constant and a mechanical damper with a constant damping coefficient.

8. A system according to claim 7, wherein:
the predetermined spring constant is configured such that the theoretical system resonance under the influence of a mechanical spring of the predetermined spring constant matches a particular off-resonant excitation frequency.

9. A system according to claim 7, wherein:
the part of the term representing a spring with a predetermined spring constant is based on the plurality of variables, which include the first variable representing excitation frequency of the transducer, a second variable representing relative position of the transducer, and a third variable representing relative velocity of the transducer.

10. A system according to claim 7, wherein:
the damping factor term is constrained by a lower bound, the lower bound greater than or equal to zero such that the damping factor term is always positive.

11. A system according to claim 8, wherein:
the damping factor term is constrained by a lower bound, the lower bound less than zero such that the damping factor term can be both negative and positive.

12. A system according to claim 1, wherein:
the predetermined parametric relation includes a term representing a damping factor for controlling the damping of the oscillating movement of the transducer, which results in the transducer force to be equivalent to the sum of a mechanical damper with a constant damping coefficient and a simulated mass.

13. A system according to claim 12, wherein:
the simulated mass is configured such that the theoretical system resonance under the influence of an equivalent physical mass matches a particular off-resonant excitation frequency.

14. A system according to claim 12, wherein:
the part of the term representing the simulated mass is based on the plurality of variables, which include the first variable representing excitation frequency of the transducer, a second variable representing relative velocity of the transducer, and a third variable representing relative acceleration of the transducer.

15. A system according to claim 12, wherein:
the damping factor term is constrained by a lower bound, the lower bound greater than or equal to zero such that the damping factor term is always positive.

16. A system according to claim 12, wherein:
the damping factor term is constrained by a lower bound, the lower bound less than zero such that the damping factor term can be both negative and positive.

17. A system according to claim 1, wherein:
the parametric relation of the control scheme is configured to extend the bandwidth of oscillatory movements of the transducer that produce oscillating electrical signals by the transducer.

18. A system according to claim 1, wherein:
the transducer comprises at least one coil and magnet that move relative to one another, wherein the at least one coil produces the oscillating electrical signal.

19. A system according to claim 18, wherein:
the at least one control element of the power electronics comprises variable load circuitry operably coupled to the output of the at least one coil.

20. A system according to claim 1, wherein:
the transducer comprises at least one piezoelectric element that produces the oscillating electrical signal.

21. A system according to claim 20, wherein:
the at least one control element of the power electronics comprises variable load circuitry operably coupled to the output of the at least one piezoelectric element.

22. A system according to claim 1, wherein:
the transducer and power electronics are configured to convert environmental vibrational energy into electrical energy, wherein the environmental vibrational energy is caused by fluid flow through a subterranean well or by drilling operations in a subterranean well.

23. A system for converting environmental vibrational energy into electrical energy for storage in at least one electrical power storage device, the system comprising:
a transducer that undergoes oscillating movement in response to the environmental vibrational energy, the transducer producing an oscillating electrical signal in response to the oscillating movement of the transducer;
power electronics, operably coupled between the transducer and the at least one electrical power storage device, the power electronics processing the oscillating electrical signal produced by the transducer, wherein the power electronics includes at least one control element having a configuration that provides variable dampening of the oscillating movement of the transducer in response to at least one control signal supplied thereto;
at least one sensor for measuring properties of the oscillations of the transducer over time; and
control electronics operably coupled to the at least one control element of the power electronics and to the at least one sensor, the control electronics carrying out a control scheme that generates and supplies the at least one control signal to the at least one control element over time in a manner that controls the at least one control element to dynamically vary the dampening of oscillations of the transducer over time, wherein the control scheme is based upon a predetermined parametric relation representing relative phase between two variables derived from the properties measured by the at least one sensor.

24. A system according to claim 23, wherein:
the two variables have a phase relationship of in-phase or 180 degrees out of phase in resonant vibrations, and the control scheme is configured to bring the two variables into a target phase relationship of in-phase or 180 degrees out of phase for cases, respectively, when the phase relationship of the two variables varies from the target phase relationship.

25. A system according to claim 23, wherein:
at least one of the variables represents a time-varying characteristic of the transducer selected from the group consisting of:
i) excitation acceleration,
ii) excitation velocity,
iii) excitation position,
iv) relative acceleration,
v) relative velocity,
vi) relative position,
vii) electrical current, and
viii) electrical voltage.

26. A system according to claim 23, wherein:
the control scheme includes first and second terms each representing constant damping factors for controlling the damping of oscillations of the transducer, wherein the first term represents a constant damping factor less than the constant damping factor represented by the second term.

27. A system according to claim 26, wherein:
the control scheme selects the constant damping factor of one of the first and second terms based upon a measure of the phase relationship of the two variables.

28. A system according to claim 27, wherein:
the measure of the phase relationship of the two variables is derived by multiplying the two variables.

29. A system according to claim 23, wherein:
the parametric relation of the control scheme is configured to extend the bandwidth of oscillations of the transducer that produce oscillating electrical signals by the transducer.

30. A system according to claim 23, wherein:
the transducer comprises at least one coil and magnet that move relative to one another, wherein the at least one coil produces the oscillating electrical signal.

31. A system according to claim 30, wherein:
the at least one control element of the power electronics comprises variable load circuitry operably coupled to the output of the at least one coil.

32. A system according to claim 23, wherein:
the transducer comprises at least one piezoelectric element that produces the oscillating electrical signal.

33. A system according to claim 32, wherein:
the at least one control element of the power electronics comprises variable load circuitry operably coupled to the output of the at least one piezoelectric element.

34. A system according to claim 21, wherein:
the transducer and power electronics are configured to convert environmental vibrational energy into electrical energy, wherein the environmental vibrational energy is caused by fluid flow through a subterranean well or by drilling operations in a subterranean well.