1. Introduction
The development of electric field measurement systems for various scientific and engineering applications, as well as for research purposes, is widespread and was extensively studied both in the past and in the last few decades. The dramatically increasing prospects of HVDC line utilization for power transmission in modern power systems has brought about the need to assess the electric environment (i.e., the electric field, the ion current density, and the space charge density) under the transmission lines to establish safety guidelines. Consequently, the knowledge of these parameters at and immediately above ground level is of great significance due to the possible effects on health and safety during operation [
1,
2,
3,
4]. Electric field estimation is also important in dusty phenomena close to the earth surface [
5], for aeronautical and astronautical projects and rocket launching [
6], and for airborne measurements to create a profile of the fair weather electric field [
7].
A major area of research associated with electric field measurement is the behavior of the electric field in meteorological and atmospheric phenomena with a special focus set on weather forecasting. In atmospheric electricity, most measured surface quantity is the vertical electric field, also known as the potential gradient [
8]. The local electric field has been measured to provide a better understanding of thundercloud electrification and the energy exchange between lightning and the atmosphere across the globe [
9]; to forecast lightning incidents; and to define an electric field threshold for lightning strikes on aircraft [
10]. In [
11], an investigation of electrical activity in thunderstorms was conducted, whereby the electrostatic field was measured during flights inside electrified clouds, in which lightning strikes occur, to provide insight beyond that of surface measurements. In fair weather conditions, the disturbance of the atmospheric electric field indicates the presence of aerosols [
12], and strong variations of the field before earthquakes and during magnetic storms have been observed [
13,
14].
Early measurements of the near-surface, fair weather atmospheric electric field in the absence of electrified clouds is around 100 V/m, maintained by the global electric circuit (GEC) [
15]; the direction of the electric field in fair weather conditions at the surface is negative. When electrified clouds are present, the field can increase from a few kV/m to tens of kV/m as it is significantly influenced by the cloud’s charge regions and lightning neutralization of the charges [
9].
Various types of DC electric field sensors have been developed for science and engineering projects. In general, electric field meters provide either the magnitude of the electric field in V/m or the potential relative to a ground reference. The primary devices employed for electric field measurement are induction probes, field mills based on charge induction, and optical sensors.
The induction probe measures the potential of a conductive plate or antenna that has been equilibrated to the local ambient field with respect to a reference ground plane. However, the employment of induction probes for long-term electric field measurements is not preferrable, due to their susceptibility to ambient space charges and the need for frequent re-zeroing in shielded conditions [
16]. As the challenging part of DC field measurement is the accumulation of random space charges near the sensors that interfere with the measurement of the electric field, for the construction of most sensors oscillating or vibrating devices are employed to create pseudo-AC fields. The conversion from a DC electric field to an AC signal benefits measurements by minimizing the effects of the space charge buildup and the spurious offsets. These devices operate either by alternately shielding or exposing a conductive plate, periodically exposing it to the local field, or by oscillating the top plate of a parallel-plate capacitor. They are advantageous compared to induction probes in that they provide higher sensitivity and they do not require constant re-zeroing.
The three prevalent sensors that adhere to this principle of operation are the shutter-type electric field mill, the vibrating plate electric field meter and the cylindrical field mill [
5,
17]. The electric field is determined by the measurement of the modulated capacitively induced charges or currents that are sensed by the conducting electrodes. The complexity in the manufacturing process of the existing electric field sensors, such as the cylindrical field mill, as well as their volume and the costly components, is not compatible with mass production. Additionally, although the presence of the space charge complicates electric field measurements in general, it is comparatively less difficult to measure electric fields at ground level than those at points above ground with cylindrical field mills. While optical sensors offer fast response and low noise, they also offer a limited dynamic range and require a complicated and costly manufacturing process [
18,
19,
20].
Therefore, the instrument of choice for the electric field measurement in atmospheric applications is the electric field mill or shutter-type field mill, since its construction is non-complex; it offers a wide dynamic measurement range; it is characterized by stable sensitivity, good signal-to-noise-ratio, and fast response; and it is suitable for mass production.
Extensive research can be found regarding the improvement of various aspects of the field mill structure, its proper installation, and its effective operation monitoring. Nevertheless, limited works address the noise introduced by the sensor interface circuitry and especially the preamplification stage, while the challenges that arise due to the weak signal induced by the sensor and its low frequency (order of Hz) have not been adequately analyzed.
In this paper, an overview of the electric field mill sensor characteristics and its most important sensitivity-related parameters is provided. In addition, the most prevalent field mill sensor interface topologies are presented. A low-power design of the front-end circuitry for a field mill sensor is proposed. The main goal of the presented work was to create a sensor interface that satisfies the specifications defined in the atmospheric electric field measurement applications.
One of the requirements taken into consideration throughout the design process of the field mill sensor and its front-end circuitry is the capability of the sensor to measure electric fields with a range sufficient to include both near-surface fair weather and foul weather electric fields. Consequently, the sensor should provide adequate sensitivity and resolution to measure fields of several orders of magnitude efficiently and accurately. Furthermore, the energy-autonomous operation of the sensor for extended periods of time is essential since the installation is frequently in isolated environments where maintenance is difficult. Furthermore, a cost-effective sensor design should not significantly load the motor of the sensor, and the interface circuit design should be easily realized, utilizing commercially available discrete ICs, while simultaneously being easily integrated using standard IC fabrication processes. Additionally, the sensor and its interface should consume minimal power in order for it to be feasible to supply their operation through energy-harvesting systems. Integrating the sensor interface design allows the customization of the preamplification stage, while simultaneously offering a low-power and compact design solution.
To address these requirements, the integrated version of the proposed sensor interface included a noise-optimized op-amp for the preamplification stage that sets the noise threshold low. For the sensor interface to be designed using discrete components while also providing the feasibility of integration, an alternative simplified approach to its design—i.e., the phase-sensitive detection stage—was adopted. An intermittent style of operation is assumed; this will prolong the energy-autonomous operation of the sensing system and contribute, along with the low-power circuitry, to the avoidance of the need for frequent battery replacement. The overall performance of the sensor interface is optimized to be comparable to the existing designs found in the literature or to outperform them in terms of resolution, sensitivity, and power consumption.
3. Proposed Electric Field Mill Sensor and Front-End Circuitry
3.1. Field Mill Design
The top view of the electric field mill sensor sensing electrodes, whose specifications determined the design of the front-end stage of the sensor interface, is depicted in
Figure 4. Two sets of sensing electrodes, set A and set B, each of which consists of two vanes, are utilized for the differential measurement. The current signal induced by set A will have a phase difference of 180° from the current induced by set B.
The main parameters of the field mill sensor appear in
Table 1.
The dynamic range of the measurement should be adequate to cover both fair-feather electric fields and fields during thunderstorms while providing high accuracy. Therefore, a field mill that can measure accurately ±20 kV/m with a resolution of 10 V/m would be sufficient.
Using the above parameters for the theoretically expected induced current for an electric field mill results in currents of the order of 320 pA for fair weather conditions, where the electric field is in the hundreds of V/m, and of approximately 64.8 nA for foul weather conditions. Therefore, a low voltage and current noise density of the op-amp used for the preamplifier realization is essential.
The application under study employs a motor that operates intermittently, i.e., is idle (standby mode) for most of the time and is only activated each time for a few seconds. The average current consumption of this style of operation is calculated by:
Therefore, for a duration of operation of 5 s and a 55-s idle state, which is an adequate duration both for acquiring a measurement and for detecting changes of the electric field before thunderstorms, the motor would consume approximately 0.65 mA. As a result, the contribution of the sensor interface circuits becomes significant in the total power consumption and should be minimized to increase the life of the energy-autonomous operation of the sensing device.
The proposed sensor interface includes an alternative simplified version of the phase-sensitive detection circuit, instead of the analog multiplier, the synchronous demodulator, or the lock-in amplifier commonly used in most works. The interface can be optionally combined with a nano-power microcontroller unit (MCU) (e.g., MSP430 series by Texas Instruments), which can control the motor activation and contribute to signal processing.
3.2. Front-End Circuits
Figure 5 depicts the proposed sensor interface, which consists of a preamplification stage, a filtering stage, and a phase-sensitive detection stage. An optical encoder, which generates a signal,
VOPTO, that indicates the exposure state of the sensing electrodes, is also employed for the electric field sign extraction.
The preamplification stage is realized by two separate channels, one for each set of sensing electrodes, each of which consists of a transimpedance amplifier that converts the current signal induced by each set into a voltage signal. The relation for the gain the transimpedance amplifier provides is:
where
R is the feedback resistance, and the subscript
X can represent either electrode set A or electrode set B. It should be noted that the transimpedance amplifier does not cause a phase shift to the original sensor signal, which translates to a 0° phase difference between the output voltage signal,
VA(
t), and the current of electrode set A,
IA(
t), and the corresponding signals of channel B.
Bandpass filtering of the signal was considered necessary since the induced signal from the electrodes, which has an amplitude proportional to the field magnitude, is expected to have a frequency equal to the rotation frequency multiplied by the number of vanes:
A multiple feedback narrow bandpass filter for each channel, whose central frequency was around the expected frequency of the signal, which is 25.5 Hz according to (10), was selected to exclude unwanted spurious signals as well as slow linear trends in the data. An additional benefit of the narrow bandpass filter is the exclusion of the odd harmonics of the signal frequency, which are caused by the shape of the stator. All frequencies, except for the frequency derived by (10), are filtered out.
The filtering stage contributes to the amplification of the voltage signal. Additionally, the filter causes a phase shift to the original sensor signal.
A voltage subtractor provides the difference, VD(t), of the signals VFA(t) and VFB(t).
The final stage of the sensor interface is the phase-sensitive detection, which, combined with the processed optical encoder signal, provides an output voltage signal that contains information on both the magnitude and the polarity of the electric field. This configuration uses the same logic as in [
9], where the optical encoder signal is processed in the digital domain in the same manner as the raw sensor signal and used as a reference to extract the field polarity by its phase difference. In this work, the realization of the optical encoder signal processing is conducted in the analog domain.
From
Figure 5, it can be observed that the phase-sensitive detection stage consists of two operational amplifiers with a shutdown feature, i.e., a high signal activates the op-amp operation, while a low signal deactivates it. The first op-amp is in inverse amplifier configuration, while the second is in a unit-gain (buffer) configuration.
As neither the transimpedance amplifier nor the voltage subtractor cause a delay to the raw sensor signal and the only phase shift is introduced by the bandpass filter, filtering the signal generated by the optical sensor in the same manner as the raw sensor signal results in an approximately 0° phase difference between the processed optical encoder signal and the filtered sensor signal for the positive electric fields and an approximately 180° phase difference for the negative electric fields.
The filtered optical encoder signal is passed through a comparator that generates a clean square pulse waveform, en, that controls the activation of the buffer amplifier and its inverted version, en, which is responsible for the activation of the inverting amplifier. When the optical encoder signal is high, the unity gain amplifier is activated by en, and the inverting amplifier is deactivated by en_; therefore, for the respected duration the output signal is the same as the input signal, VD(t). When VOPTO is low, it activates the inverting amplifier, and the output signal is the inverted input signal for the respected duration. A low-pass filter is employed to convert the derived voltage signal, VS(t), into a DC voltage, Vout, which is proportional to the ambient electric field. Observing the phase-sensitive detection configuration, the derived conclusion is that voltage Vout will be higher than the mid-supply level (VDD/2) when the electric field is positive, and vice versa.
4. Experimental Measurements
The experimental measurements to verify the operation of the proposed field mill interface depicted in
Figure 5 were conducted on a PCB containing a discrete implementation of the proposed sensor interface using a prototype 3D-printed electric field mill sensor. The motor had a rotation frequency of 12.75 Hz when operating at 5 V. A calibration setup with similar structure to the one depicted in
Figure 3 was employed to create known electric fields ranging from −6.5 kV/m to +6.5 kV/m.
Table 2 includes the values of all the passive components used.
Discrete ICs that are characterized by low noise and low power consumption were used for the preamplification stage, the bandpass filter, and the phase-sensitive detection stage, whose models appear in
Figure 5.
The oscilloscope measurements appear in
Figure 6a–d, where the top two measurement windows in
Figure 6a,b show the sensor signal and the optical encoder signal before (
V1(
t) and
VOPTO) and after processing (
VD(
t) and en) for the positive and negative electric fields, respectively. The bottom two measurement windows in
Figure 6c,d show the differential signal after subtraction,
VD(
t), along with the processed optical encoder signal, en, and the output waveform of the phase-sensitive detection circuit,
VS(
t), before low-pass filtering. To analyze the amplitude-frequency response, an FFT of the voltage waveform at each processing stage is performed using a digital storage oscilloscope (SIGLENT,
SDS1102CML+). The derived diagrams are depicted in
Figure 6e.
As can be observed in
Figure 6e, before filtering the voltage signal
VA(
t) contains several unwanted frequency components caused by both environmental (external) noise sources and internal circuitry noise contribution. After filtering, the voltage signal
VFA(t) contains only the desired frequency of 25.5 Hz amplified by the bandpass filter and a 50 Hz common-mode noise picked up by both electrodes. After the voltage subtraction, the remaining frequency is the theoretically expected frequency of 25.5 Hz, and the signal is doubled in amplitude. The 50 Hz common-mode noise is greatly eliminated.
Figure 7a depicts the discrete implementation of the proposed sensor interface and the calibration setup. In the diagram of
Figure 7b, the output voltage is plotted against the applied electric field.
5. Optimized IC Sensor Interface
The optimized integrated sensor interface appears in
Figure 8. During simulations the theoretical model’s Equation (4) were used to simulate the current induced by the sensing electrodes.
The implementation of the integrated version of the sensor interface has some slight design differences from the discrete version, in order to improve the silicon area consumption. An alternative differential transimpedance amplifier configuration [
37] was selected for the integrated implementation instead of the classic configuration that consists of one separate amplification channel for each set of electrodes since it provides differential-to-single-ended signal amplification, employing only one operational amplifier and, therefore, decreasing the power consumption, the silicon real estate, and the need for an op-amp subtractor, as well as the need for a second bandpass filter, is avoided. The relation for the gain the transimpedance amplifier provides, is:
Low noise introduction in the preamplification stage is of great significance in applications where accurate measurements are necessary since the minimum electric field that can be measured depends on this noise. Therefore, a custom-designed operational amplifier optimized for the specific requirements of low noise combined with low power consumption was employed for this stage (
Figure 9a) [
38]. As a low-voltage operation combined with low power consumption was required for the integrated version of the proposed sensor interface, a compact amplifier that combines low voltage and power efficiency was needed. This amplifier’s output rail-to-rail configuration is necessary to maximize the dynamic range and, by extension, the signal-to-noise ratio. The floating class AB driver prevents the contribution of the class AB stage to the noise and offset introduced by the amplifier [
39].
The dimensions of the optimized op-amp MOSFETs appear in
Table 4. The Miller compensation capacitor C
M is 12 pF, and the Miller compensation resistor is 1.5 kΩ.
Advanced noise simulations were conducted (
Figure 9b,c), and the final design parameters appear in
Table 4. The signal induced from the sensing electrodes has a low frequency (in the order of tens of Hz). At this frequency band, a flicker noise is the dominant source of noise, and therefore, input transistors with a large
WL area are required to reduce its contribution. A trade-off between the current consumption and the noise introduced by the differential pairs exists, where an increased current results in decreased noise [
40]. Low power is also a requirement; hence, the current at the differential pair of the op-amp was selected as a good balance between the consumption and the input voltage noise.
The fabrication technology transistor models (XFAB 0.18 μm) that were employed for the design of the amplifier are supplied in the BSIM3V3 format. BSIM3V3 contains advanced thermal and flicker noise models [
41].
To calculate the noise output in closed-loop op-amps, the noisy op-amp is modeled as a noiseless op-amp with an equivalent input-referred, mean-square noise voltage. The effect of all the noise sources in the circuit is represented by a single source at the input.
A hands-on simulation-based methodology for the noise optimization of the amplifier utilizing the advanced noise simulation tools provided by Cadence was conducted. Noise simulators provide the possibility to determine not only the equivalent input-referred noise, but also the type of noise sources that constitute the extracted noise and the noise contribution of each component. Using these features combined with the theory on noise sources, the efficient noise optimization of the amplifier was feasible.
At first, during the initial advanced noise analyses of the operational amplifier before the optimization of the design, it was confirmed that the prevailing noise sources were the flicker noise and the thermal noise. This observation is consistent with the operational amplifier noise theory. In addition, according to [
39], the noise of the amplifier is determined mainly by the input transistors and the summing circuit since the floating class-AB control is shifted into the summing circuit. This was confirmed by simulations of the noise contribution of each component.
The dominant noise source in the frequency domain of the induced sensor signal—which is in the order of tens of Hz—is the flicker (or 1/
f) noise. Flicker noise is described by the following equation:
where
K is a process dependent constant,
W and
L are the transistor’s width and length, respectively, and
Cox is its gate capacitance per unit area. This was confirmed by noise simulations of the amplifier, where in low frequencies the noise decreases at a rate of 10 dB/decade, which is in accordance with the flicker noise behavior. Thermal noise is caused by the random thermally excited vibration of the charge carrier in a conductor. It is present in all resistive elements, and it is dependent on the absolute temperature. A the MOSFET level, it is described by the following equation:
where k is the Boltzmann’s constant,
T is the temperature in Kelvin,
rds is the channel resistance, and
gm is the device transconductance [
42]. The input-referred voltage noise density at 25.5 Hz at this stage was 1.25 μV/√Hz.
In the next step, the noise contribution of each transistor at the frequency of interest (25.5 Hz) was extracted to determine the prevailing noise contributors. The differential pair transistors presented the highest noise contribution percentage; therefore, according to the flicker noise equation, increasing their WL area, would decrease their noise contribution. For that reason, parametric simulations sweeping the WL area of the input transistors were conducted to determine the dimensions of each transistor finger. The voltage noise density at 25.5 Hz at this stage was 363.2 nV/√Hz.
To further decrease the noise contributions of the input stage, the width of the transistors was increased. Parametric simulations sweeping the number of fingers (n_f) of the input transistors were conducted to determine the total width of each transistor. The voltage noise density at 25.5 Hz at this stage was 233.3 nV/√Hz.
At this point in the simulations, the noise contribution of the input transistors became comparable to the noise contribution of the transistors of the summing stage of the amplifier; therefore, a similar approach was adopted to decrease the noise introduced by this stage. Parametric simulations sweeping the area of transistors M7-M10 and M11-M14 were conducted. The voltage noise density at 25.5 Hz at this stage was 66.7 nV/√Hz.
At this point, no further optimization of the design in terms of noise was conducted since the noise was adequately low. The final transistor dimensions appear in
Table 4.
Figure 10 sums up all the optimization stages conducted, where
n_f is the number of fingers of the input transistors, a is the constant that is multiplied by both the nominal width and the nominal length of the input transistors, and b is the constant that is multiplied by both the nominal width and the nominal length of the summing circuit transistors M7-M10 and M11-M14.
The nominal dimensions for transistors M1–M2 and M7–M10 were W = 30 μm and L = 2 μm, and for M3–M4 and M11–M14, they were W = 10 μm and L = 2 μm.
Table 5 includes the noise parameters as well as the power consumption of the sensor interface.
To interpret the noise introduced by the preamplification stage, the transient noise of the preamplification stage is provided for the input and output node. In
Figure 11, the transient noise analysis provides the peak-to-peak voltage noise.
The related amplitude, which is about 140 μVp-p at the preamplifier’s input and approximately 1.6 mVp-p at its output, is significant and concerns the preamplifier after optimization, depicting that the noise optimization of the preamplification stage is crucial.
Simulations including the transient noise set by the internal noise of the preamplification stage combined with the filtering stage were conducted—assuming an ideal noiseless environment—and a spectrum analysis to deconstruct the time domain representation of the derived signals into the frequency domain representation was performed (
Figure 12). The spectrum analysis after the filtering stage showed significant elimination of the internal noise that is introduced by the preamplification stage.
The total power consumption of the sensor interface contributes significantly to the consumption of the total field mill sensor system due to the intermittent style of operation of the motor; therefore, to prolong its energy-autonomous operation, the sensing interface should consume as little as possible. From
Table 4, it is derived that the preamplifier is the main contributor to the total power consumption, and therefore, its consumption should be limited.
The integration of the on-chip bandpass filter is a challenging process since dealing with low frequencies generally results in the filters requiring passive components with high values that are difficult to integrate. By selecting the minimum possible values for the capacitors (100 pF) with a high area capacitance (MIM capacitor) included in the filter configuration and high values for the resistors, as well as highly resistive components provided by the fabrication technology (high ohmic N+ poly resistor), the integration challenge is overcome. Low power 3.3 V nMOS and pMOS transistors were used throughout the whole interface design.
Figure 13 depicts the schematic of the op-amp circuit with the enable feature, where the control signals en_p and en_n are externally controlled by signals en and en_ for the buffer and the inverting amplifier configuration, respectively. The same op-amp design was employed for the implementation of the bandpass filter op-amp as well as the comparator, but without the enable feature.
The main waveforms of the phase-sensitive detection stage appear in
Figure 14.
The output voltage after low-pass filtering is plotted versus the electric field in
Figure 15. The maximum achieved sensitivity is 45 mV/kV/m.
The main characteristics and performance parameters of the proposed sensor interface are compared in
Table 6 with the other works in the literature that employ an electric field mill.
From
Table 6, it is apparent that the proposed sensor and its sensing interface consume significantly less power than the field mill sensing systems that are employed in foul weather measurements, even when they operate in continuous mode, whereas they provide high sensitivity combined with a single power supply of 3 V. The low power consumption aspect is crucial when an extended energy-autonomous operation is required. It should be noted that while an ADC was not employed in the design of the processing channel, the theoretical resolution of the interface, should a 16-bit ADC be used, was calculated and included in
Table 6 for comparison reasons.
Using a relatively larger sensing electrode area combined with the low noise preamplification stage, an improved sensitivity and resolution are achieved, and accurate fair weather measurements become feasible even at a low rotation frequency of the motor. Therefore, a carefully designed field mill sensor that follows the guidelines discussed, combined with the proposed sensor interface would be ideal for measurements of both fair and foul weather conditions.