1. Introduction
In recent years, the network has undergone a significant revolution. On one hand, there have been notable changes in electric assets. For instance, photovoltaic plants are becoming increasingly common on the roofs of industrial, commercial, and private buildings [
1,
2]. The number of electric vehicles has rocketed in the last few years thanks to commercial maneuvers and country-based incentives [
3,
4]. Electric scooters have also become a popular means of private transportation, requiring a power supply for charging [
5,
6]. On the other hand, monitoring techniques and associated software have been evolving. The advancement and widespread adoption of artificial intelligence (AI) and machine learning (ML) techniques have brought about significant changes in power system monitoring techniques [
7,
8]. The common thread among all these developments is instrumentation. The new assets can be monitored through sensors and distributed measurement systems deployed throughout the power network. These sensors provide measurements for the monitoring system and the learning algorithms used in ML. Therefore, it can be concluded that sensors and instrumentation play a significant role.
Voltage and current measurements are the primary activities conducted at network nodes and branches. However, numerous other quantities are also measured, such as humidity, temperature, power, energy, pressure, etc. [
9,
10,
11]. All these electrical and environmental parameters are measured to obtain the most accurate health status of the grid. As a result, studies and standards dedicated to sensors and instrumentation are published on a daily basis. Regarding standards, the IEC 61869 series comprises fifteen main documents that pertain to instrument transformers (ITs). Readers can refer to documents IEC 61869-1 and -6 [
12,
13] or the generic specifications applicable to ITs and low-power instrument transformers (LPITs), respectively. In recent literature, a calibration technique associated with partial discharge measurement using a high-frequency current transformer (HFCT) is described in [
14]. ITs with digital output are calibrated in [
15]. The authors in [
16] characterize a measurement chain that includes LPIT and the phasor measurement unit (PMU). The aim in [
17] is to accurately measure distorted signals, while [
18] deals with design solutions for ITs, and [
19] focuses on failure investigation techniques to be applied to ITs.
In this paper, the IT considered is the Rogowski coil. Due to its features, the Rogowski coil is widely adopted in various applications beyond power systems. Consequently, Rogowski coil modeling has been a research topic for the past few decades, and numerous works can be found in the literature. For example, an electro-thermal model written in VHDL-AMS language and based on finite element analysis (FEA) is described in [
20]. In [
21], the authors design a Rogowski coil from scratch for three-phase current measurement in electric motors. A highly detailed model is presented in [
22], relying on equations that require the complete geometry of the device and the space in which it is installed. A four-layer PCB-based Rogowski coil is designed in [
23]. A simple model, although based on several preliminary measurements, is detailed in [
24]. The same authors introduced a new testing signal, known as the sinc response, to be measured by the Rogowski coil to obtain its input-output response [
25]. The high-frequency behavior of the Rogowski coil is studied in [
26] using the classical lumped model, while [
27] explores another high-frequency modeling approach based on a black-box conceptualization of the Rogowski derived from preliminary measurements on the device.
From the literature review, it becomes apparent that there is a lack of very simple models that do not require extensive preliminary measurements. Furthermore, obtaining information on the internal structure and design of the device is not always feasible. This makes the modeling even more complicated and not generalizable for every commercial Rogowski coil. To address this issue, this paper introduces a new and straightforward method to obtain the frequency model of a commercial Rogowski coil. In practice, it is common to purchase a device for various applications rather than designing a custom Rogowski coil from scratch. The proposed model is based on a few inputs obtained from the device’s datasheets, and a set of equations is developed to derive a transfer function dependent on the cross-sectional shape. The lack of knowledge about commercial devices’ internal structure is an obstacle to many existing models. Subsequently, the obtained model is implemented in the Matlab 2023a environment for preliminary efficacy testing. Finally, the model is experimentally validated through measurements performed on commercial Rogowski coils. The results clearly demonstrate the validity and accuracy of the frequency model, confirming its potential for implementation in simulation environments such as real-time power network simulators.
The remainder of this paper is structured as follows:
Section 2 includes the theoretical concepts and motivation behind the work.
Section 3 is dedicated entirely to the model description and contextualization. Model validation is presented in
Section 4. Finally,
Section 5 concludes this paper and provides suggestions for future research.
3. The Modelling Procedure
The proposed modeling procedure starts with the geometrical measurements on the RC and ends with the transfer function (TF) that correlates the frequency relation between the output voltage and the input current. The procedure is summarized in the flowchart depicted in
Figure 2.
The first step is to measure the inner diameter 2
a and the diameter of the cross-section
d. The frequency modeling procedure is developed for three cross-section geometries: circular, square/rectangular, and oval. In the latter two geometries, an extra measurement is needed, as explained in
Figure 3. As for the oval case, commercial devices do not present ideal oval shapes. Therefore, the geometry is often schematized, as shown in
Figure 3.
The second step consists of extracting some information from the datasheet. The two key parameters needed for the modeling are the accuracy class and the transformation ratio (TR). Other typical parameters, given in the RC datasheet, are listed in
Table 1. Note that the TF also exploits the information about the rated burden. However, if not provided by the manufacturer, the standard value of 2 MΩ given in [
13] can be used.
The third step allows us to obtain the mutual induction,
starting from (1). Substituting the derivative over time with the angular speed
:
where
is the generic frequency and
is the rated transformation ratio (at 50 Hz). Once
is found, it is possible to obtain the
TR at each frequency of interest. For this assumption, the parasitic parameters will be neglected.
The fourth step regards geometry. Using the quantities obtained in the first step, the perimeter and the cross-section area of the RC can be calculated.
At this point, using consolidated expressions of the mutual inductance [
29,
30], the number of turns of the RC can be obtained from the reversed formula. Note that, from a practical perspective, the lack of knowledge of the number of turns is one of the limiting factors of RC modeling. Hence:
where
is the overall permeability, and
,
, and
are the number of turns of the rectangular, circular, and oval cross-section, respectively. They can be used to obtain, in addition to the geometrical parameters, the length
of the single coil, section
, and radius
of the wire used to wound the RC.
The last step is the calculation of the parasitic parameters of the RC:
where
is the overall permittivity (to be fixed depending on the used material),
and
are the parasitic capacitance and inductance, respectively. The suffix
is used to switch among the cross-section parameters (
). As for the winding resistance,
, it can be easily obtained with the second Ohm’s law, measured, or extracted from the datasheet. Of course, there will be datasheets with extra information that can be used instead of some parameter expressions. However, the proposed method has been generalized to avoid any exceptional cases of missing (extra) inputs. A small note on permeability and permittivity is necessary. The general notation has been used to highlight that, depending on the RC adopted, the user may need to insert or not insert the relative permeability/permittivity. Further tests on these specific parameters demonstrated an almost negligible effect of the permittivity/permeability value on the accuracy parameters.
With the results from (2)–(9), the RC transfer function in the “
s” domain can be obtained as:
where
is the rated burden, or better, the impedance of the measuring device connected in cascade to the RC. The
TF in (10) must be customized with the parameters associated with the studied cross-section. To better contextualize the
TF in (10), a schematic with the equivalent circuit is depicted in
Figure 4. All the symbols in the picture were previously described.