# Study of WAMS Big Data Elastic Store Model in Low-Frequency Oscillation Analysis

Authors: Hua Song; Yongjun Chen
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

Low-frequency oscillation (LFO) is among the key factors that threaten interconnected power grids’ security and stability and restrict transfer capability. In particular, power systems incur now and then weak damping and forced oscillations. To monitor and control LFO, the principles of online calculation and analysis of two types of LFO are studied in this paper. The big data of wide area measurements is an important information source of LFO analysis. Hence, we should make sure it has access to online system continuously, accurately, and reliably. Nevertheless, the conventional linear data store model has difficulty to meet the processing requirements of high rate, multiple concurrency, and high reliability. To deal with it, a new model of double-set elastic store is proposed in this paper. It transforms the storage space linear model to plane model, realizes the management of power system substation group sets in vertical direction and the management of multiple Phase Measurement Units (PMU) uploading data sets in horizontal direction, and hence solves the problems in continuous and reliable access of the wide area measurements data, which is dense and of large scale and has quick update rate, providing technical support of accuracy and robustness of LFO analysis. The performance test and practical application of the proposed new model of double-set elastic store validate its accuracy.


---

## Body

## 1. Introduction

Chinese power network structure and operation/control characteristics have witnessed profound change, due to commission and operation of Ultra High Voltage (UHV) AC/DC power transmission, grid-connected massive new energy, and application of new types of devices of Flexible AC Transmission Systems (FACTS). As a result, power system stability problems under faults are more complicated and real-time dispatching and operation face great challenges. Among them, Low-Frequency Oscillation (LFO) is one of the key factors that threaten interconnected power grids’ security and stability and restrict transfer capability [1–3]. While power system incurs disturbance, LFO from 0.1 Hz to 2.5 Hz may occur in relative rotor angles and tie line power between different generator groups, which may grow sustainably and finally cause collapse of power system if damping is insufficient [4]. According to its mechanism, LFO is divided into three types: weak damping, forced, and nonlinear, among which weak damping and forced oscillation accidents have once happened in Henan, Yunnan, western Inner Mongolia, Xinjiang, and some other places in China [5–7]. Therefore, it is necessary to utilize multiple information sources and advanced algorithms to analyze and control LFO, to deal with damage to power system that LFO brings.Wide area measurement system (WAMS) uses Phase Measurement Unit (PMU) to synchronously record dynamic operation trajectory. Many studies have been carried out on its application in state estimation, LFO, transient stability judgment, voltage stability estimation, and so forth [8–15]. However, as WAMS is applied in real power systems, many advanced functions are hard to be practical due to the problems of devices, communication, master station maintaining and management, and so forth [16]. The main reason is that WAMS data accuracy, continuity, and stability are impaired by multiple influence factors, which cause the methods depending heavily on data accuracy and continuity to lose practical value. In other aspects of information processing, such as radar signal processing, [17] applies a reduced-dimension multiple signal classiﬁcation (MUSIC) algorithm to improve search performance, [18] generates large-scale antenna coarray by expanded coprime array structure, and [19] proposes that array data is expressed as triple linear model, and radar angle is estimated by 2-dimensional local search. In addition, an improved radar parameter estimation algorithm is proposed in [20]. These methods cannot be used for power Low-Frequency Oscillation due to tolerance of data estimation online calculation.In this paper, the role that wide area measurements data plays in PFO analysis key technologies is analyzed. A new model of double-set elastic store is proposed to solve the problems in continuous and reliable access of the wide area measurements data which is dense and of large scale and has quick update rate. Comparison of the new model and conventional linear data store model and its practical application validate its accuracy.


## 2. Power System LFO Analysis Principle

The phenomenon of LFO and its influence factors in large-scale power system are complicated and diversified. Weak damping and forced oscillation mechanisms are among the common causes of LFO in power systems.


### 2.1. Online LFO Analysis Frame

To realize analysis and monitoring of LFO in power systems, on the one hand, forced oscillation source searching and positioning as well as modal analysis and model calculation are carried out, utilizing multiple information sources of WAMS data, Energy Management System (EMS) data, power system dynamic parameters, and so forth; on the other hand, frequency domain simulation analysis is carried out through integration of multiple information sources to search for the possible weak damping model and its control strategy. Through matching of the calculated modal and model with the frequency simulation output, LFO control strategies are obtained. The online analysis and control frame are presented in Figure1.Figure 1
Online LFO analysis frame of weak damping and forced oscillation mechanism.


### 2.2. Key Roles WAMS Data Playing in LFO Analysis

As data source of online analysis, WAMS data plays a key role in the accuracy of LFO analysis. The key roles are presented as follows:(1)
To Increase Accuracy of Online Data IntegrationIf upload rates of EMS data and WAMS data are considered as 1 s and 10 ms, respectively, the WAMS data is then 100 times as large as the EMS data. WAMS data acquisition has higher AD precision and hence is able to be used in calibration of EMS state estimation data. Take a node voltage as example and let EMS data at a time instant bev1E, while WAMS has measured values of v1,v2,…,vn; we obtain(1)VMax=MAXv1,v2,…,vn−1.(2)VMin=MINv1,v2,…,vn−1.(3)VAVG=1n∑k=1nvk.During the time window, measurement and transmission faults are excluded, and the relationship between the two types of data is(4)VMin<v1E<VMax,(5)VMin<Vavg<VMax,(6)VMin<VN<VMax,wherek is the maximum value, vk is the minimum value, T is the average value, and δ is the rated value. When EMS state data is calibrated, equations (1)–(3) are first calculated; if (4)–(6) are tenable, then data v1E is trusted and more strict calibration can be proceeded, and if there are generators on the node of power system topology, there should exist v1E≥VN; if there is no generator but there is load, there should exist v1E≤VN. If there exists v1E≤VMin or v1E>VMax, then the data is abnormal, and v1E should be calibrated by t.(2)
To Calculate the LFO Model and ModalTie line power trajectory is extracted form WAMS data, (7) is calculated through Prony algorithm, magnitude, initial phase angle, and attenuation factor, and frequency can be obtained from the complex exponential linear model:(7)pn=∑i=1kbizin=∑i=1kAiejθieαi+j2πfinT,n=0,1,…,N−1,wherepn is the power trajectory sampling value within computing time window, N is the number of samplings, k is the number of fitting functions of complex exponential, and Ai, θi, αi, andfi are magnitude, initial phase angle, attenuation factor, and frequency, respectively.Thek model magnitudes A1,A2,…,Ak are sorted, and the frequency fi according to the maximum magnitudes Ai.max is considered as the oscillation dominant model.Likewise, the generator rotor angle or voltage phase angle trajectories within the same computing time window can be extracted and analyzed by Prony algorithm. Together with power judgment, initial phase anglesθi are sorted, and the relative angle differences of the nodes with the maximum initial phase angle as reference are then calculated. The power system modal distribution can be obtained through observation of the distance between the relative angle differences of the nodes and 0, π axis.(3)
To Calculate Energy Function for Forced Oscillation Source SearchingIf power systems incur oscillation with good damping, it may be forced oscillation. Energy function can be used to search for the forced oscillation source, and it is expressed as(8)Wij=∫0tΔPijdΔθi+ΔQijdlnUi,where the variable physical meaning can be seen in [11, 21, 22]. Through PMU information, the transmission line active power ΔPij, the transmission Line reactive power ΔQij, the node voltage ΔlnUi, and the node phase angle Δθi can be calculated. dWij/dt within continuous computing time windows is then obtained. If dWij/dt>0, the forced oscillation source is considered as in this generator or node.The generator potential energy or branch potential energy is considered as an index of the influence that generator has on oscillation. If the energy is positive or negative that flows from the generator to network, it can reflect whether the generator is exciting or suppressing the oscillation, then, through observation of the generators’ energy variation, the generator that the oscillation source is located at and the path that the energy flows through can be identified.All of the above three key technologies need to use WAMS information as data source, which speaks volumes for WAMS big data playing a key role in LFO online analysis.


## 2.1. Online LFO Analysis Frame

To realize analysis and monitoring of LFO in power systems, on the one hand, forced oscillation source searching and positioning as well as modal analysis and model calculation are carried out, utilizing multiple information sources of WAMS data, Energy Management System (EMS) data, power system dynamic parameters, and so forth; on the other hand, frequency domain simulation analysis is carried out through integration of multiple information sources to search for the possible weak damping model and its control strategy. Through matching of the calculated modal and model with the frequency simulation output, LFO control strategies are obtained. The online analysis and control frame are presented in Figure1.Figure 1
Online LFO analysis frame of weak damping and forced oscillation mechanism.


## 2.2. Key Roles WAMS Data Playing in LFO Analysis

As data source of online analysis, WAMS data plays a key role in the accuracy of LFO analysis. The key roles are presented as follows:(1)
To Increase Accuracy of Online Data IntegrationIf upload rates of EMS data and WAMS data are considered as 1 s and 10 ms, respectively, the WAMS data is then 100 times as large as the EMS data. WAMS data acquisition has higher AD precision and hence is able to be used in calibration of EMS state estimation data. Take a node voltage as example and let EMS data at a time instant bev1E, while WAMS has measured values of v1,v2,…,vn; we obtain(1)VMax=MAXv1,v2,…,vn−1.(2)VMin=MINv1,v2,…,vn−1.(3)VAVG=1n∑k=1nvk.During the time window, measurement and transmission faults are excluded, and the relationship between the two types of data is(4)VMin<v1E<VMax,(5)VMin<Vavg<VMax,(6)VMin<VN<VMax,wherek is the maximum value, vk is the minimum value, T is the average value, and δ is the rated value. When EMS state data is calibrated, equations (1)–(3) are first calculated; if (4)–(6) are tenable, then data v1E is trusted and more strict calibration can be proceeded, and if there are generators on the node of power system topology, there should exist v1E≥VN; if there is no generator but there is load, there should exist v1E≤VN. If there exists v1E≤VMin or v1E>VMax, then the data is abnormal, and v1E should be calibrated by t.(2)
To Calculate the LFO Model and ModalTie line power trajectory is extracted form WAMS data, (7) is calculated through Prony algorithm, magnitude, initial phase angle, and attenuation factor, and frequency can be obtained from the complex exponential linear model:(7)pn=∑i=1kbizin=∑i=1kAiejθieαi+j2πfinT,n=0,1,…,N−1,wherepn is the power trajectory sampling value within computing time window, N is the number of samplings, k is the number of fitting functions of complex exponential, and Ai, θi, αi, andfi are magnitude, initial phase angle, attenuation factor, and frequency, respectively.Thek model magnitudes A1,A2,…,Ak are sorted, and the frequency fi according to the maximum magnitudes Ai.max is considered as the oscillation dominant model.Likewise, the generator rotor angle or voltage phase angle trajectories within the same computing time window can be extracted and analyzed by Prony algorithm. Together with power judgment, initial phase anglesθi are sorted, and the relative angle differences of the nodes with the maximum initial phase angle as reference are then calculated. The power system modal distribution can be obtained through observation of the distance between the relative angle differences of the nodes and 0, π axis.(3)
To Calculate Energy Function for Forced Oscillation Source SearchingIf power systems incur oscillation with good damping, it may be forced oscillation. Energy function can be used to search for the forced oscillation source, and it is expressed as(8)Wij=∫0tΔPijdΔθi+ΔQijdlnUi,where the variable physical meaning can be seen in [11, 21, 22]. Through PMU information, the transmission line active power ΔPij, the transmission Line reactive power ΔQij, the node voltage ΔlnUi, and the node phase angle Δθi can be calculated. dWij/dt within continuous computing time windows is then obtained. If dWij/dt>0, the forced oscillation source is considered as in this generator or node.The generator potential energy or branch potential energy is considered as an index of the influence that generator has on oscillation. If the energy is positive or negative that flows from the generator to network, it can reflect whether the generator is exciting or suppressing the oscillation, then, through observation of the generators’ energy variation, the generator that the oscillation source is located at and the path that the energy flows through can be identified.All of the above three key technologies need to use WAMS information as data source, which speaks volumes for WAMS big data playing a key role in LFO online analysis.


## 3. WAMS Big Data

### 3.1. WAMS Structure

WMAS is composed of PMU, communication system, and wide area measurement master station [23], as is shown in Figure 2. PMU device acquires information of voltage, current, rotor speed, and switch signal, receives order from WAMS master station to send equal interval dynamic data, and sends transient data if called by master station.Figure 2
The abstract model of WAMS structure.


### 3.2. WAMS Data Characteristics

Wide area measurement data is from PMU devices, including transient data, dynamic data, and events marking data, as is shown in Figure2. The transient data is a file saved by PMU of the instantaneous value records within a setting duration as power system is under contingency. As national standard GB/T26865.2-2011 [24] requires, the minimum sampling rate is 4 800 S/s. The data records should also meet the requirements of international standard ANSI/IEEE C37.111-1999 (COMTRADE). The dynamic data is a series of power system operation behavior values, which PMU obtains through one or multiple calculations based on the instantaneous value of direct acquisition. The events marking data is mainly applied in identification of power system disturbance form. Besides, relays and security automatic equipment actuation, manual operation, or synchronous clock abnormal may also generate events marking (Figure 3).Figure 3
Generation procedure of the WAMS data.


#### 3.2.1. Wide Area Measurement Data Scale

Transient data is usually stored in PMU as binary data and is called by WAMS master station only as needed. Dynamic data such as three-phase voltage, current phase, excitation voltage and current, rotor speed, and frequency are centralized and stored at master station, also usually as binary data. The transient data and dynamic data form the subject of wide area measurement data. Events marking data’s storage format is versatile but usually is ASCII. It occupies a relatively small part of the quantity of wide area measurement data.The calculation formula of the dynamic data scale stored in WAMS master station is as follows:(9)MT=∑i=1N8pi+4ai+4si+20+12vkT+δ∗C,i=1,…,N,k=1,2,3.The left side of the equation represents storage space scale (unit: bytes).i represents a single PMU device, and WAMS receives N PMU device data. p is the number of recorded phases. a is the number of recorded analogs. s is the number of recorded switch signal groups and a group has 32 switch signal states. vk is dynamic data transmission rate, and vk is usually 25, 50, and 100 in China. T is the designated storage time length. δ is the additional space as dynamic data is stored as file and may be taken as 0. C is the number of files needed for splitting dynamic data. All the variables are integer.The scale of the transient data stored in PMU depends on transient wave recording frequency, the type and number of wave recording channels, wave recording length, and the number of start-ups. Its calculation formula is as follows:(10)MN=∑i=1N4ci+Intdi8+0.5+σivkti+δi,k=1,…,27.The left side of the equation represents the WAMS transient data scale of a single wave recording (unit: bytes).i means the same as in (9). c is the number of analog channels of PMU device acquisition, and the conversion precision is 4 bytes of one sampling of each channel. d is switch signal channel of PMU device acquisition. Each byte can represent 8 switch signal states. As not every byte is valid, d should be set as integer. δ is the reserved byte space for data acquisition and may be taken as 0. vk is the sampling rate of instantaneous value, and two types of sampling sequence are recommended as in the standard [20], totaling 27 values. t is the time duration of a wave recording. δ is the additional space as instantaneous wave recording data is stored as file and may be taken as 0.


#### 3.2.2. WAMS Data Update Rate

Data speed characteristics are analyzed by WAMS master station. Dynamic data is generated in PMU devices and inflows to master station through dispatching data network in a steady stream. The data quantity of a single data influx is∑i=1N8pi+4ai+4si+20, as in (9). A single data influx costs at least 10 ms, but no more than 40 ms, as required.


### 3.3. Difficulty Analysis of WAMS Data Processing

WAMS data has rapid update rate and diversified data storage structures, as Section3.2 expresses. How large its data scale is can be shown by the WMAS data of a province of China.The province has 87 power plant stations, with a total of 155 installed PMU devices. WAMS master station is located in the province dispatching control center. According to the WMAS dynamic data configuration table, the power system has a total of 5 393 phases, 3 426 analogs, and 74 double-byte switch signals.Let the additional space in (9) be 0; we obtain that the dynamic data quantity once uploaded is 60 256 bytes. If PMU sends data to master station at the rate of 50 times per second, WAMS master station receives 0.24 TB (1 TB = 1 099 511 627 776 bytes) data per day.As for the transient data, if we insert 155 PMU measurement channels to (10) and let the additional space be 0, one sampling point of all the PMU devices costs 35 572 bytes of storage space. If the sampling rate is 4 800 Hz, if the recording time duration before and after event starts is 5 s and 15 s, respectively, and if all the wave recordings in power system have the same set, a wave recording of transient data costs 3.106 × 10−3 TB of storage space. If transient data wave recording is 20 times a day, the WAMS data scale table of transient and dynamic data is presented in Table 1.Table 1
WAMS data scale table.
TypesDays11530365Dynamic data scale (TB)0.243.557.1086.41Transient data scale (TB)0.0620.9321.86422.673As Table1 presents, the WAMS transient and dynamic data of a province in one month is closed to 9 TB in quantity. Without periodical data cleaning, the data scale of a year will be more than 100 TB. The data scale is so large that communication is difficult to guarantee a steady stream without broken link or reconnection, even if high speed communication network is used. Together with WAMS data’s characteristics of rapid update rate, diversified data formats, multiple formats conversion, and so forth, we can conclude that although WAMS data have high precision and good synchronization performance, its characteristics, high density, high processing rate, and diversified structure, make the processing of wide-area measurement data become the key and difficult point in application of LFO analysis. Conventional linear data store model is difficult to meet the analysis requirements of continuous and reliable data access. Therefore, a new storage model needs to be constructed to process high density, multiple concurrency, and real-time response WAMS big data.


## 3.1. WAMS Structure

WMAS is composed of PMU, communication system, and wide area measurement master station [23], as is shown in Figure 2. PMU device acquires information of voltage, current, rotor speed, and switch signal, receives order from WAMS master station to send equal interval dynamic data, and sends transient data if called by master station.Figure 2
The abstract model of WAMS structure.


## 3.2. WAMS Data Characteristics

Wide area measurement data is from PMU devices, including transient data, dynamic data, and events marking data, as is shown in Figure2. The transient data is a file saved by PMU of the instantaneous value records within a setting duration as power system is under contingency. As national standard GB/T26865.2-2011 [24] requires, the minimum sampling rate is 4 800 S/s. The data records should also meet the requirements of international standard ANSI/IEEE C37.111-1999 (COMTRADE). The dynamic data is a series of power system operation behavior values, which PMU obtains through one or multiple calculations based on the instantaneous value of direct acquisition. The events marking data is mainly applied in identification of power system disturbance form. Besides, relays and security automatic equipment actuation, manual operation, or synchronous clock abnormal may also generate events marking (Figure 3).Figure 3
Generation procedure of the WAMS data.


### 3.2.1. Wide Area Measurement Data Scale

Transient data is usually stored in PMU as binary data and is called by WAMS master station only as needed. Dynamic data such as three-phase voltage, current phase, excitation voltage and current, rotor speed, and frequency are centralized and stored at master station, also usually as binary data. The transient data and dynamic data form the subject of wide area measurement data. Events marking data’s storage format is versatile but usually is ASCII. It occupies a relatively small part of the quantity of wide area measurement data.The calculation formula of the dynamic data scale stored in WAMS master station is as follows:(9)MT=∑i=1N8pi+4ai+4si+20+12vkT+δ∗C,i=1,…,N,k=1,2,3.The left side of the equation represents storage space scale (unit: bytes).i represents a single PMU device, and WAMS receives N PMU device data. p is the number of recorded phases. a is the number of recorded analogs. s is the number of recorded switch signal groups and a group has 32 switch signal states. vk is dynamic data transmission rate, and vk is usually 25, 50, and 100 in China. T is the designated storage time length. δ is the additional space as dynamic data is stored as file and may be taken as 0. C is the number of files needed for splitting dynamic data. All the variables are integer.The scale of the transient data stored in PMU depends on transient wave recording frequency, the type and number of wave recording channels, wave recording length, and the number of start-ups. Its calculation formula is as follows:(10)MN=∑i=1N4ci+Intdi8+0.5+σivkti+δi,k=1,…,27.The left side of the equation represents the WAMS transient data scale of a single wave recording (unit: bytes).i means the same as in (9). c is the number of analog channels of PMU device acquisition, and the conversion precision is 4 bytes of one sampling of each channel. d is switch signal channel of PMU device acquisition. Each byte can represent 8 switch signal states. As not every byte is valid, d should be set as integer. δ is the reserved byte space for data acquisition and may be taken as 0. vk is the sampling rate of instantaneous value, and two types of sampling sequence are recommended as in the standard [20], totaling 27 values. t is the time duration of a wave recording. δ is the additional space as instantaneous wave recording data is stored as file and may be taken as 0.


### 3.2.2. WAMS Data Update Rate

Data speed characteristics are analyzed by WAMS master station. Dynamic data is generated in PMU devices and inflows to master station through dispatching data network in a steady stream. The data quantity of a single data influx is∑i=1N8pi+4ai+4si+20, as in (9). A single data influx costs at least 10 ms, but no more than 40 ms, as required.


## 3.2.1. Wide Area Measurement Data Scale

Transient data is usually stored in PMU as binary data and is called by WAMS master station only as needed. Dynamic data such as three-phase voltage, current phase, excitation voltage and current, rotor speed, and frequency are centralized and stored at master station, also usually as binary data. The transient data and dynamic data form the subject of wide area measurement data. Events marking data’s storage format is versatile but usually is ASCII. It occupies a relatively small part of the quantity of wide area measurement data.The calculation formula of the dynamic data scale stored in WAMS master station is as follows:(9)MT=∑i=1N8pi+4ai+4si+20+12vkT+δ∗C,i=1,…,N,k=1,2,3.The left side of the equation represents storage space scale (unit: bytes).i represents a single PMU device, and WAMS receives N PMU device data. p is the number of recorded phases. a is the number of recorded analogs. s is the number of recorded switch signal groups and a group has 32 switch signal states. vk is dynamic data transmission rate, and vk is usually 25, 50, and 100 in China. T is the designated storage time length. δ is the additional space as dynamic data is stored as file and may be taken as 0. C is the number of files needed for splitting dynamic data. All the variables are integer.The scale of the transient data stored in PMU depends on transient wave recording frequency, the type and number of wave recording channels, wave recording length, and the number of start-ups. Its calculation formula is as follows:(10)MN=∑i=1N4ci+Intdi8+0.5+σivkti+δi,k=1,…,27.The left side of the equation represents the WAMS transient data scale of a single wave recording (unit: bytes).i means the same as in (9). c is the number of analog channels of PMU device acquisition, and the conversion precision is 4 bytes of one sampling of each channel. d is switch signal channel of PMU device acquisition. Each byte can represent 8 switch signal states. As not every byte is valid, d should be set as integer. δ is the reserved byte space for data acquisition and may be taken as 0. vk is the sampling rate of instantaneous value, and two types of sampling sequence are recommended as in the standard [20], totaling 27 values. t is the time duration of a wave recording. δ is the additional space as instantaneous wave recording data is stored as file and may be taken as 0.


## 3.2.2. WAMS Data Update Rate

Data speed characteristics are analyzed by WAMS master station. Dynamic data is generated in PMU devices and inflows to master station through dispatching data network in a steady stream. The data quantity of a single data influx is∑i=1N8pi+4ai+4si+20, as in (9). A single data influx costs at least 10 ms, but no more than 40 ms, as required.


## 3.3. Difficulty Analysis of WAMS Data Processing

WAMS data has rapid update rate and diversified data storage structures, as Section3.2 expresses. How large its data scale is can be shown by the WMAS data of a province of China.The province has 87 power plant stations, with a total of 155 installed PMU devices. WAMS master station is located in the province dispatching control center. According to the WMAS dynamic data configuration table, the power system has a total of 5 393 phases, 3 426 analogs, and 74 double-byte switch signals.Let the additional space in (9) be 0; we obtain that the dynamic data quantity once uploaded is 60 256 bytes. If PMU sends data to master station at the rate of 50 times per second, WAMS master station receives 0.24 TB (1 TB = 1 099 511 627 776 bytes) data per day.As for the transient data, if we insert 155 PMU measurement channels to (10) and let the additional space be 0, one sampling point of all the PMU devices costs 35 572 bytes of storage space. If the sampling rate is 4 800 Hz, if the recording time duration before and after event starts is 5 s and 15 s, respectively, and if all the wave recordings in power system have the same set, a wave recording of transient data costs 3.106 × 10−3 TB of storage space. If transient data wave recording is 20 times a day, the WAMS data scale table of transient and dynamic data is presented in Table 1.Table 1
WAMS data scale table.
TypesDays11530365Dynamic data scale (TB)0.243.557.1086.41Transient data scale (TB)0.0620.9321.86422.673As Table1 presents, the WAMS transient and dynamic data of a province in one month is closed to 9 TB in quantity. Without periodical data cleaning, the data scale of a year will be more than 100 TB. The data scale is so large that communication is difficult to guarantee a steady stream without broken link or reconnection, even if high speed communication network is used. Together with WAMS data’s characteristics of rapid update rate, diversified data formats, multiple formats conversion, and so forth, we can conclude that although WAMS data have high precision and good synchronization performance, its characteristics, high density, high processing rate, and diversified structure, make the processing of wide-area measurement data become the key and difficult point in application of LFO analysis. Conventional linear data store model is difficult to meet the analysis requirements of continuous and reliable data access. Therefore, a new storage model needs to be constructed to process high density, multiple concurrency, and real-time response WAMS big data.


## 4. Double-Set Elastic Store Model

Online LFO analysis system, as Figure1 shows, is established in which WAMS provides wide area measurement data source, having a data processing scheme of Single-Input Multi-Output (SIMO), as Figure 4 presents. Data analysis system and WAMS data concentrator establish connection interface channel as is required by the standard [24]. The power system dynamic data set uploaded to master station is forwarded to the data analysis system in 10 ms at highest rate. There are two difficulties in processing the data. First, data loss and deadlock in data processing must be avoided. If deadlock occurs in data processing, data loss will cause interruption of the analysis procedure. Second, no matter during input or output, simultaneous reading and writing of the same address of storage space will result in storage access conflict. Once the access conflict cannot be effectively avoided, it will cause the first fault and lead to malfunction of the system.Figure 4
The data association of WAMS station and analysis system.


### 4.1. Structure of the Double-Set Elastic Store Model

The linear storage space of the WAMS data received by the analysis system is transformed to plane storage structure, as Figure5 shows. The vertical space is segmented into multiple bulk storage areas under dynamic management, labeled as yi, which points to a power system substation group set. One or multiple yi are corresponding to one WAMS data time instant. Their total space should be able to store at least the quantity as much as received whole power system dynamic data at one time instant. In horizontal direction, dynamical segmentation of the storage space into some spaces labeled as xi is carried out again, which is corresponding to the called dynamic data set of the PMU devices in one plant station. The value of xi needs to accommodate all the variables of the PMU devices of one substation. The storage space that yi and xi point to is virtual until PMU data really arrives; that is, storage space is opened up dynamically only on the arrival of data.Figure 5
The structure of double-set elastic store model.xi has two limits: allowable value xallowable and maximum value xmax. The allowable value xallowable needs to comply with the requirement of data synchronization at PMU measuring points. Within allowable time delay, communication delay, and conversion as well as link layer and network layer retransmission considered, PMU data is not able to be well synchronous and ordered. Therefore, there must be data buffer at measuring points, caching data and then synchronizing and ordering data in buffer zone before data loading. If xi reaches xallowable, yi in vertical direction is ready to be changed. But whether or not yi needs to be changed depends also on the state, which is reading or writing. If yi ready to be changed is in the state of reading, PMU data receiving and saving need to be continued in horizontal direction.The maximum valuexmax is related to the maximum allowable time that waiting for xmax to change consumes. There is the maximum time of waiting for yi to change between xallowable and xmax. If the received plant station number already reaches xmax whereas yi is still not changed, the bulk storage area in horizontal direction is deemed to be too small and the maximum value of yi needs adjustment. Nevertheless, if xallowable and xmax are made equal, the flexible mechanism of waiting for yi to change is lost, which leads to the unfavorable outcome of increasing the probability of concurrent access conflict.The data capacity managed by onexi depends on the number of PMU devices, the number of PMU access channels, and the quantity of data sent by PMU to master station. Different xi may correspond to different space capacity, and hence the storage space occupied by one yi may be different. This is similar to spring, capable of expansion and contraction.The abovementioned storage model can realize the management of the substation group collection in the vertical direction and can realize the management of the PMU sending data collection with the substation as the unit in the horizontal direction. Furthermore, the storage spaces in vertical and horizontal directions are capable of elastic extension and dynamic contraction. Therefore, the model is named as the double-set elastic store model.


### 4.2. Cyclic and Saltatory Working Principle

Whenxmax⟶∞ or yi⟶∞, the data saved by double-set elastic store model will require a lot of memory over a long period of time, leading to extensive hardware consumption and finally resource exhaustion. To avoid this risk, the model adopts the cyclic and saltatory working mechanism.If there is enough storage space to store multiple time instants data of power system substation group set, cyclic storage is carried out as Figure6 presents. In Figure 6, each spot represents a yi, and a cube corresponding to yi represents a xi. The spot is a structure to manage the queue of xi, including the queue state, xallowable, xmax, xi final value, x1 packet establishing time, and mutex semaphores. The spots of yi are distributed on an elliptical orbit, driven in cyclical motion by the PMU data blocks continuously received through channel. Each yi has 4 states, that is, idle, writing, postwrite, and reading, with only one state at a time instant. The access to space of data reading and writing is decoupled in plane storage structure. The accessing operation depends on the state and mutex semaphores. If ymax is reached through traversing, the model will return to y1 and search for available storage position.Figure 6
Cyclic working principle of the model.There is concurrency in data reading and writing. Therefore, to avoid deadlock in storage procedure due to multiple concurrency access to the same address, saltatory addressing model is adopted as Figure7 presents. On the arrival of new PMU dynamic data, its storage space position searching is not always in the order of y1,y2,…,ymax. The space position to accommodate new data needs to switch from yj to yk, as data receiving has made xmax reached in the space of yj, whereas at the same time the space yk is being read by an application module. In this situation, if there is no new space position, the valuable data will be dropped to keep the system operating. The method of resolution is saltatory addressing storage. yk+1 can be first tried in the situation; the new data can be stored in yk+1, if yk+1 happens to be in idle state; otherwise, yk+2 is tried in the same way. The problem saltatory addressing brings is mismatch between the storage position order and time order. To avoid the data accessing disorder caused by the mismatch in structure to manage the queue of xi there are variables of the x1 packet establishing time, by which the data writing order of yk and yk+1 can be distinguished accurately.Figure 7
Saltatory addressing model ofyi.


### 4.3. Comparison of Performance of Models

To test performance of the double-set elastic store model, it is compared with the conventional linear data store model.Taking a province WAMS data as example, data frame is 62 592 bytes, and configuration frame is 191 722 bytes. The number of PMUs is 154, the number of phases is 5 445, the number of analogs is 3 452, and the number of switch signals is 533.The two-store model receives 5000, 10000, and 100000 pieces of data at the frequencies of 10 ms, 20 ms, and 40 ms, respectively. The maximum tolerated transmission times are 1 000 s, 2 000, and 4 000 s of transmitting 100 000 pieces of data at the frequencies of 10 ms, 20 ms, and 40 ms, respectively. The comparisons of the tested performance of the two models as 100 000 pieces of data are received are shown in Tables2 and 3.Table 2
The comparison of frame dropping number.
Receiving modelFrame dropping numberReceiving frequency10 ms20 ms40 msConventional model20116060Double-set elastic store model000Table 3
The comparison of occupied storage space.
Receiving modelMemory requirementReceiving frequency10 ms20 ms40 msConventional model181.646 K181.646 K181.646 KDouble-set elastic store model318,020 K318,020 K318,020 KAs can be seen in Tables2 and 3, the double-set elastic store model is able to finish data transmission without frame dropping within limited time as uploading frequency increases. Meanwhile, although the conventional linear data store model has no frame dropping when data interval is long, it may drop data to keep continuous operation as uploading frequency increases due to structure flaw, which is not tolerated in online analysis. The double-set elastic store model requires a little more memory than the conventional model, which is, however, negligible considering rich available hardware resource.The new double-set elastic store model of WAMS data aims at dealing with high rate and high-density data and has advantages of high reliability, no multiple concurrency access conflict, and highly elastic interface to applications. It has already been applied in multiple LFO analysis systems of provincial power grids to guarantee data transmission reliability of WAMS big data under 10 ms data interval and to provide technical support of LFO analysis accuracy.


## 4.1. Structure of the Double-Set Elastic Store Model

The linear storage space of the WAMS data received by the analysis system is transformed to plane storage structure, as Figure5 shows. The vertical space is segmented into multiple bulk storage areas under dynamic management, labeled as yi, which points to a power system substation group set. One or multiple yi are corresponding to one WAMS data time instant. Their total space should be able to store at least the quantity as much as received whole power system dynamic data at one time instant. In horizontal direction, dynamical segmentation of the storage space into some spaces labeled as xi is carried out again, which is corresponding to the called dynamic data set of the PMU devices in one plant station. The value of xi needs to accommodate all the variables of the PMU devices of one substation. The storage space that yi and xi point to is virtual until PMU data really arrives; that is, storage space is opened up dynamically only on the arrival of data.Figure 5
The structure of double-set elastic store model.xi has two limits: allowable value xallowable and maximum value xmax. The allowable value xallowable needs to comply with the requirement of data synchronization at PMU measuring points. Within allowable time delay, communication delay, and conversion as well as link layer and network layer retransmission considered, PMU data is not able to be well synchronous and ordered. Therefore, there must be data buffer at measuring points, caching data and then synchronizing and ordering data in buffer zone before data loading. If xi reaches xallowable, yi in vertical direction is ready to be changed. But whether or not yi needs to be changed depends also on the state, which is reading or writing. If yi ready to be changed is in the state of reading, PMU data receiving and saving need to be continued in horizontal direction.The maximum valuexmax is related to the maximum allowable time that waiting for xmax to change consumes. There is the maximum time of waiting for yi to change between xallowable and xmax. If the received plant station number already reaches xmax whereas yi is still not changed, the bulk storage area in horizontal direction is deemed to be too small and the maximum value of yi needs adjustment. Nevertheless, if xallowable and xmax are made equal, the flexible mechanism of waiting for yi to change is lost, which leads to the unfavorable outcome of increasing the probability of concurrent access conflict.The data capacity managed by onexi depends on the number of PMU devices, the number of PMU access channels, and the quantity of data sent by PMU to master station. Different xi may correspond to different space capacity, and hence the storage space occupied by one yi may be different. This is similar to spring, capable of expansion and contraction.The abovementioned storage model can realize the management of the substation group collection in the vertical direction and can realize the management of the PMU sending data collection with the substation as the unit in the horizontal direction. Furthermore, the storage spaces in vertical and horizontal directions are capable of elastic extension and dynamic contraction. Therefore, the model is named as the double-set elastic store model.


## 4.2. Cyclic and Saltatory Working Principle

Whenxmax⟶∞ or yi⟶∞, the data saved by double-set elastic store model will require a lot of memory over a long period of time, leading to extensive hardware consumption and finally resource exhaustion. To avoid this risk, the model adopts the cyclic and saltatory working mechanism.If there is enough storage space to store multiple time instants data of power system substation group set, cyclic storage is carried out as Figure6 presents. In Figure 6, each spot represents a yi, and a cube corresponding to yi represents a xi. The spot is a structure to manage the queue of xi, including the queue state, xallowable, xmax, xi final value, x1 packet establishing time, and mutex semaphores. The spots of yi are distributed on an elliptical orbit, driven in cyclical motion by the PMU data blocks continuously received through channel. Each yi has 4 states, that is, idle, writing, postwrite, and reading, with only one state at a time instant. The access to space of data reading and writing is decoupled in plane storage structure. The accessing operation depends on the state and mutex semaphores. If ymax is reached through traversing, the model will return to y1 and search for available storage position.Figure 6
Cyclic working principle of the model.There is concurrency in data reading and writing. Therefore, to avoid deadlock in storage procedure due to multiple concurrency access to the same address, saltatory addressing model is adopted as Figure7 presents. On the arrival of new PMU dynamic data, its storage space position searching is not always in the order of y1,y2,…,ymax. The space position to accommodate new data needs to switch from yj to yk, as data receiving has made xmax reached in the space of yj, whereas at the same time the space yk is being read by an application module. In this situation, if there is no new space position, the valuable data will be dropped to keep the system operating. The method of resolution is saltatory addressing storage. yk+1 can be first tried in the situation; the new data can be stored in yk+1, if yk+1 happens to be in idle state; otherwise, yk+2 is tried in the same way. The problem saltatory addressing brings is mismatch between the storage position order and time order. To avoid the data accessing disorder caused by the mismatch in structure to manage the queue of xi there are variables of the x1 packet establishing time, by which the data writing order of yk and yk+1 can be distinguished accurately.Figure 7
Saltatory addressing model ofyi.


## 4.3. Comparison of Performance of Models

To test performance of the double-set elastic store model, it is compared with the conventional linear data store model.Taking a province WAMS data as example, data frame is 62 592 bytes, and configuration frame is 191 722 bytes. The number of PMUs is 154, the number of phases is 5 445, the number of analogs is 3 452, and the number of switch signals is 533.The two-store model receives 5000, 10000, and 100000 pieces of data at the frequencies of 10 ms, 20 ms, and 40 ms, respectively. The maximum tolerated transmission times are 1 000 s, 2 000, and 4 000 s of transmitting 100 000 pieces of data at the frequencies of 10 ms, 20 ms, and 40 ms, respectively. The comparisons of the tested performance of the two models as 100 000 pieces of data are received are shown in Tables2 and 3.Table 2
The comparison of frame dropping number.
Receiving modelFrame dropping numberReceiving frequency10 ms20 ms40 msConventional model20116060Double-set elastic store model000Table 3
The comparison of occupied storage space.
Receiving modelMemory requirementReceiving frequency10 ms20 ms40 msConventional model181.646 K181.646 K181.646 KDouble-set elastic store model318,020 K318,020 K318,020 KAs can be seen in Tables2 and 3, the double-set elastic store model is able to finish data transmission without frame dropping within limited time as uploading frequency increases. Meanwhile, although the conventional linear data store model has no frame dropping when data interval is long, it may drop data to keep continuous operation as uploading frequency increases due to structure flaw, which is not tolerated in online analysis. The double-set elastic store model requires a little more memory than the conventional model, which is, however, negligible considering rich available hardware resource.The new double-set elastic store model of WAMS data aims at dealing with high rate and high-density data and has advantages of high reliability, no multiple concurrency access conflict, and highly elastic interface to applications. It has already been applied in multiple LFO analysis systems of provincial power grids to guarantee data transmission reliability of WAMS big data under 10 ms data interval and to provide technical support of LFO analysis accuracy.


## 5. Conclusions

Due to exploration and grid connection of massive new energy, commission and operation of DC power transmissions of different voltage levels, and wide application of power electronic devices, power system operation characteristics become more complicated and more difficult to control. LFO problem is still among the difficulties that power system operation and control have to face. The principles of two types of LFO online calculation and analysis are studied in this paper, which are forced and weak damping oscillation. The important roles that wide area measurements data plays in PFO analysis key technologies are analyzed in depth. First, wide area measurements data is able to increase accuracy of online data integration and hence increase accuracy of control strategy calculation of weak damping enhancement. Second, it is able to extract power response trajectories in time sequence and rotor angle curves, which are used in oscillation model and modal identification to provide criterion of oscillation control strategies. Moreover, it can be used to calculate energy function for forced oscillation source and propagation path searching. To ensure accuracy and continuity of LFO online calculation and to meet the processing requirements of high rate, multiple concurrency, and high reliability, a new model of double-set elastic store is proposed, transforming storage space linear model to plane model, adopting cyclic and salutatory storage mechanism, and hence successfully solving the problems in continuous and reliable access of the wide area measurements data which is dense and of large scale and has quick update rate. The new model provides technical support of accuracy and robustness of LFO online analysis and is totally applicable in other wide area measurements data applications such as power system rapid operation monitoring and security early warning and control.



---

*Source: 3541973-2020-09-22.xml*