Next Article in Journal
A Hybrid Adaptive Controller Applied for Oscillating System
Previous Article in Journal
A Review of Reliability and Fault Analysis Methods for Heavy Equipment and Their Components Used in Mining
Previous Article in Special Issue
Bioconversion of the Brown Tunisian Seaweed Halopteris scoparia: Application to Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Health Monitoring of Pressure Regulating Stations in Gas Distribution Networks Using Mathematical Models

by
Shivom Sharma
1,
Iftekhar A. Karimi
1,*,
Shamsuzzaman Farooq
1,*,
Lakshminarayanan Samavedham
1 and
Rajagopalan Srinivasan
2
1
Department of Chemical & Biomolecular Engineering, National University of Singapore, Singapore 117585, Singapore
2
Department of Chemical Engineering, Indian Institute of Technology Madras, Chennai 600036, India
*
Authors to whom correspondence should be addressed.
Energies 2022, 15(17), 6264; https://doi.org/10.3390/en15176264
Submission received: 21 July 2022 / Revised: 21 August 2022 / Accepted: 24 August 2022 / Published: 28 August 2022
(This article belongs to the Special Issue Design, Integration and Optimization of Biorefineries)

Abstract

:
Many cities have extensive distribution networks that supply natural or town gas to domestic, industrial, and power plant consumers. A typical network may have hundreds of pressure regulating stations that are of different types and capacities, but most legacy networks are sparsely instrumented. The reliability of these stations is the first priority for ensuring uninterrupted gas supplies; hence, condition monitoring and prescriptive maintenance are critical. In this study, mathematical models were developed for two types of commonly used regulators: spring-loaded and lever-type regulators. We also considered three faults that are typically of interest: filter choking, valve seat damage, and diaphragm deterioration. The proposed methodologies used the available measured data and mathematical models to diagnose faults, track prognoses, and estimate the remaining useful life of the regulators. The applicability of our proposed methodologies was demonstrated using real data from an existing distribution network. To facilitate industrial use, the methodologies were packaged into a user-friendly dashboard that could act as an interface with the operational database and display the health status of the regulators.

1. Introduction

Huge networks of pipelines, called gas distribution networks (GDNs), that have hundreds of pressure regulating stations (PRSs) are used to supply natural or town gas to industrial, commercial, and residential users. A PRS steps down its inlet gas pressure and maintains its outlet gas pressure at a desired set point, regardless of the gas demand profile. Gas distribution companies carefully monitor PRSs for both gradual condition deterioration and abrupt faults. Any problems need to be identified and rectified quickly to avoid disruptions in gas supplies to consumers as the companies may face regulatory penalties for such disruptions and/or lose consumer confidence.
A PRS has two components: a gas filter and a pressure regulator (PR) (Figure 1). PRs can be divided into two categories: spring-loaded (or self-operated or direct acting) regulators and pilot regulators. A spring-loaded PR has three parts: a restriction device (orifice), a sensing device (diaphragm), and a loading device (spring). Pilot regulators are used for better control accuracy and the spring is replaced by a more constant diaphragm loading force (e.g., gas pressure). PRSs can have single, dual or multiple stages. Multiple stages offer steadier outlet pressures.
Most studies on PRSs have focused on modeling and stability issues. In the literature, three approaches have been used to model PRs: an equation-based first principles approach [1], a data-driven approach [2], and a bond graph simulation technique [3,4]. Nabi et al. [5] developed a mathematical model for dome-loaded PRs and explored the impact of several geometrical and operational parameters on their performance. Rami et al. [6] developed mathematical models for direct acting and pilot-controlled PRs and identified the parameters that were responsible for their instability. Zafer and Luecke [7] developed a dynamic mathematical model for a spring-loaded PRs, studied the possible causes of vibrations, and suggested design modifications. Wang et al. [8] developed a dynamic model and self-tuning method for pneumatic pressure regulating stations. Several other studies [9,10,11,12] have also addressed the dynamic modeling of PRs. Recently, several studies have focused on the development of data-driven methodologies for the monitoring and fault diagnosis of gas pressure regulators [13,14,15,16,17].
Leo et al. [1] simulated a GDN to study the health monitoring of PRSs. They studied four faults: spring decay, filter choking, diaphragm malfunction, and valve seat damage. Three scenarios were studied for filter choking and spring decay: Pout (outlet pressure) and flow data at the same sampling rates, Pin (intermediate pressure) and Pout data at the same sampling rates, and Pout and cumulative flow data at different sampling rates. To monitor the valve seat condition, the weekly average of the daily maximum Pout values was used, whereas the hourly mean of Pout was used to monitor the diaphragm condition.
Based on our industry interactions, PRSs can develop five types of health degradations and/or faults: the overloading of the downstream network, spring decay, filter choking, valve seat damage, and diaphragm malfunctions. Spring decay, filter choking, and the overloading of downstream networks cause downstream under-pressure, whereas valve seat damage results in over-pressure. PRSs can experience slow (e.g., spring decay and diaphragm malfunctions) and/or fast (e.g., filter choking and valve seat damage) degradation. Over time, the springs lose their stiffness and the diaphragms lose their flexibility of movement due to fatigue (cyclic stress). The filters that are placed upstream of PRs to remove particles can choke with time. In the case of valve seat damage, the flows through regulators can stop completely.
Over the years, the maintenance paradigm for physical systems has improved from reactive maintenance (fixing faults after they occur) to preventive (performing routine repairs), then predictive (predicting failures and then repairing them), and finally prescriptive maintenance (predicting failures and recommending action). PRSs are usually monitored for gradual health deterioration and abrupt faults. Hence, real-time monitoring and automated prescriptive maintenance are crucial for the smooth operation of GDNs. As seen in Figure 1, pressure can be measured at three locations: upstream of the gas filter, between the filter and regulator, and downstream of the PR. Most existing PRSs are sparsely instrumented. The intermediate pressure (Pin in Figure 1) is not measured in many PRSs. Most often, the data that are collected from PRSs are not automatically transferred to the respective control rooms. Depending on the availability of pressure and flow data, either model-based or data-based methods can be applied for PRS health monitoring.
Currently, gas distribution companies perform routine preventive maintenance on PRSs. Maintenance teams visit PRSs every six months for inspection and repair. In the worst-case scenario, the sudden failure of a PRS can lead to unexpected disturbances in the gas supply. This study developed methodologies for the prescriptive maintenance of filters, valve seats, and diaphragms. Spring decay was not considered as it occurs rarely and very slowly (the average lifetime of springs is more than 10 years). Firstly, first principles models were developed for spring-loaded and lever-type PRs. Note that lever-type PRs also use springs. Data that were provided by the PRS vendors were used to calibrate the developed mathematical models. The first principles health monitoring models could be computationally expensive. Hence, as an alternative, the first principles models were simplified to estimate flows more quickly with comparable accuracy. Most gas distribution companies only measure flows at the end-use points (households, industry, power plants, etc.) and the huge number of end users makes it difficult to estimate the total gas consumption within a geographical area. The developed first principles models could provide estimates of the gas flows through PRSs. Gas companies could use this information for advanced planning to avoid the overloading of downstream networks.
The main contributions of the present study are summarized below:
  • Mathematical models for two different regulator types: spring-loaded and level-type regulators (the spring-loaded regulators had a significantly different valve seat geometry than those that were studied by Leo et al. [1]);
  • Precise PR dimensions and calibration data that were obtained from the PRS vendors to determine the mathematical model parameters (in contrast, Leo et al. [1] assumed some parameters);
  • In contrast to Leo et al. [1], this work used real data from a large existing GDN comprising both industrial and residential end-users (the GDN did not measure the flows through the PRs and rarely measured Pin);
  • Significantly modified and advanced monitoring algorithms to address the following challenges for more robust and improved inferences: (a) higher noise levels in industrial data compared to the white noise that has been added to synthetic data, (b) the absence of some data from the GDN, and (c) pressure sensor errors in real data;
  • A filter monitoring methodology that used Pin (or Po when Pin was not available) and Pout, whereas valve seat and diaphragm monitoring only required Pout (these advances are detailed in the Fault Monitoring Methodologies section).
In the prescriptive maintenance of PRSs, the estimation of remaining useful life (RUL) is crucial for uninterrupted gas supplies. The proposed health monitoring methodologies could predict RUL values for the different parts of PRSs. A dashboard (pressure regulator health monitor, PRHM) was developed to track all of the PRSs from one place. It had a user-friendly GUI (graphical user interface) and it could be integrated with online operational data. In addition to supervising various faults, the PRHM also monitored sudden abnormalities in the PRSs.

2. Mathematical Model for Spring-Loaded PRs

Figure 2a presents a schematic of a gas PR series 280 (HON 280, [18]). This spring-loaded PR is available in three sizes: 2″ (DN 50), 3″ (DN 80), and 4″ (DN 100). Leo et al. [1] presented a first principles model for the HON 280 that considered the force balance on the valve seat and the diaphragm assembly, as well as the mass and energy balances across the flow chamber. Since PRs have much faster dynamics than downstream gas networks, they reach steady state instantaneously. Hence, it was appropriate to develop a steady-state model for PRs.
This work considered a different geometry of PR valve seats from that in Leo et al. [1] (Figure 2b); However, the model assumptions and equations for the HON 280 PR were the same as those that were used by Leo et al. [1], except for the force balance. The PR model is represented by Equations (1)–(6). The diaphragm and valve seat moved up and down together. Equation (1) presents the force balance on the diaphragm and valve seat assemblies. The chambers above and below the diaphragm were at atmospheric pressure (Patm) and Pout, respectively. The spring was pre-compressed with an initial force of F0. When the spring was extended (x > 0), then the spring pulled the diaphragm upward. In Equations (1) and (2), Pin is the inlet pressure, K is the spring constant, AD is the diaphragm area, AO is the plunger area, Ap is the peripheral flow area, and d is the diameter of the orifice. In Figure 2b, D, do and dt represent the diameters of the diaphragm, plunger, and shaft, respectively.
The mass and energy balances across the flow chamber that were used by Leo et al. [1] are given by Equations (3) and (4). In Equations (3)–(6), Tin is the inlet temperature, Tout is the outlet temperature, mout is the outlet mass flow, Cd is the orifice discharge coefficient, R is the gas constant, ρ is the density of the gas, γ is the gas heat capacity ratio, μ is the Joule–Thomson expansion coefficient, and N is as defined by Andersen [19].
Kx = F 0 + P atm A D P out A D + P in A O P out A O
A p =   π dx
T out = T in μ ( P in P out )
m out = C d A p NP in γ [ 2 / ( γ + 1 ) ] ( γ + 1 ) / 2 ( γ 1 ) / RT in
N = { 1 if   P out / P in ( 2 γ + 1 ) γ / ( γ 1 )   ( P out P in ) 2 / γ ( P out P in ) ( γ + 1 ) / γ ( γ 1 2 ) ( 2 γ 1 ) ( γ 1 ) / ( γ + 1 ) otherwise
Q = m out × 3600 / ρ
In the PR model, x, Tout, μ, γ, and mout were unknowns. Equation (1) could be solved to compute the spring extension x. Next, Ap could be calculated using Equation (2), in which μ, γ, and ρ depended on the temperature, pressure, and gas composition (X). Following the work of Leo et al. [1], neural network (NN) models were used to estimate μ, γ, and ρ (Equations (7)–(9)) to avoid the use of the computationally expensive Peng–Robinson equation of state. The ANN models were two-layer feedforward networks with one hidden layer and an output layer. In total, 10000 data points were generated for the temperature, pressure, and composition using Latin hypercube sampling. Interested readers are referred to Leo et al. [1] for more details on the ANNs. Equation (3) could be solved iteratively to obtain Tout and μ. Then, mout was computed using Equations (4) and (6) was used to calculate the hourly volumetric flow through the PR.
μ =   f NN ( P out ,   T out , X )
γ =   f NN ( P out ,   T out , X )
ρ =   f NN ( P out ,   T out , X )
In the mathematical models for spring-loaded (Equations (1)–(6)) and lever-type (Equations (3)–(6) and (10)–(15)) PRs, Equation (3) required most of the computation time. In this equation, there were two unknowns (Tout and μ). The initial value of Tout was assumed to estimate μ using the neural network model (Equation (7)). After that, the μ value was used to compute Tout by solving Equation (3). When the assumed Tout value (that was used to estimate μ in the neural network model) and the computed Tout value (that was obtained after solving Equation (3)) were within the threshold limits, they were used as the final values of Tout and μ.
In this study, the gas mixture that flowed through the PR was about 60% hydrogen, 12% methane, 1% ethane, 4.2% pentane, 4% CO, 17% CO2, 3.8% N2, and 1% O2. It was found that Tout was very close to Tin. Hence, the first principles (PR) models were simplified by assuming Tout = Tin. With hundreds of PRs and hundreds of thousands of data samples for each regulator, this simplification to the mathematical models could significantly reduce computation time.

3. Mathematical Model for Lever-Type PRs

Figure 3a presents a schematic of a gas PR series 277 (HON 277). This lever-type PR is available in a 2″ size and is used as a service regulator. The working mechanisms of the HON 277 and HON 280 PRs are very different. The HON 277 PR has a lever between the diaphragm and plunger.
Figure 3b shows the torque balance for this PR. F1 is the force on the plunger that is due to the inlet gas pressure, whereas F2 is the net force on the diaphragm. Equation (10) presents the torque balance around the stationary point in Figure 3b. This equation could be solved to obtain the value of the spring extension (x). For a given value of x, the plunger displacement (y) could be calculated using Equations (11)–(14). Here, Equation (11) provided the initial angle (θ1) of the lever when x = 0. The final angle (θ2) of the lever for x was determined using Equation (12). Additionally, θd represents the angular movement of the lever that corresponded to the spring extension x (Equation (13) and Figure 3c). The displacement of the plunger could be computed using Equation (14). Finally, the peripheral flow area (Ap) is given by Equation (15).
The mass and energy balances across the flow chamber were the same for both regulators. The complete mathematical model for the lever-type PR contained Equations (3)–(6) and (10)–(15). The solution procedure was also the same for both models.
( Kx + P out A D F 0 P atm A D ) L 2 = P in A o L 1
θ 1 = tan 1 ( M / L 2 )
θ 2 = tan 1 ( [ M x ] / L 2 )
θ d =   θ 1   θ 2
y =   L 1 tan ( θ d )
A p =   π dy

4. Model Parameters

Table 1 provides the dimensions of the regulators (dt, d, do, D, L1, L2, and M). In this work, a Cd (coefficient of discharge) of 0.87 was used. The initial force on the spring (F0) was computed using Equation (16). It required set points for P out and the diaphragm area. To estimate the spring constant (K), the vendor data were used. Table S1 in the Supplementary Materials presents the vendor data (Pin, Pout, and flow) for all of the PRs. Figure 4 shows the algorithm of the optimization-based methodology that was used to estimate K. The optimizer (fmincon in MATLAB) provided the initial K for the first principles (PR) model. The model used Pin and flow data from the vendors and other parameters (dimensions, AD, Cd, and F0, as shown in Table 1) to predict Pout. An objective function (OF) was calculated using the Pout value that was predicted by the model and the Pout value that was predicted by the vendor, as shown in Figure 4. The optimizer decided on the next K value using this objective function value. The optimization search was terminated after converging to a final K value that minimized the objective function.
F 0 =   P out ,   SP × A D

5. Fault Monitoring Methodologies

This section discusses the methodologies that were developed for fault monitoring in PRSs. Firstly, data preprocessing is discussed before addressing the potential faults in PRSs, namely filter choking, valve seat damage, and diaphragm deterioration. In general, gas distributors do not measure flows through PRSs. Furthermore, most gas companies only measure P 0 and P out . Figure 5 summarizes the methodologies that were developed to monitor the three faults (filter choking, valve seat damage, and diaphragm deterioration).

5.1. Data Preprocessing

  • Data were imported from Excel into MATLAB. Each regulator had an Excel file with the date, time, inlet pressure (P0), intermediate pressure (Pin), and outlet pressure (Pout) data. Some did not have Pin data.
  • When one pressure value was missing for a data point, then all of its pressure values were replaced by those for the previous data point. A pressure value was considered an outlier when it was more than three standard deviations away from the mean value. Outlier pressure values were replaced by the previous pressure value.
  • Data were checked for valid pressure drops across the filter (ΔPf). The number (SN-ΔPf) of data samples with negative ΔPf values was computed. When SN-ΔPf was lower than 0.1% of the number of all samples, then those faulty pressure values were replaced by the mean values. Otherwise, New Pin = Old Pin − max (ΔPf).
  • Important information was extracted, such as the start and end dates of the data, the number of days of data, and the number of samples for each day.
  • Finally, the pressure data were visualized by plotting P0, Pin, Pout, and ΔPf against time.

5.2. Filter Choking

As discussed earlier, both spring decay and filter choking cause low pressure in downstream networks. In general, spring decay is much slower than filter choking. For relatively short periods, the spring constant can be assumed to be constant. Hence, filter choking can be considered as the main reason for pressure decrease in downstream networks. An ideal regulator maintains the outlet pressure at Pout,SP for any flow and Pin values. However, Pout decreases with the flow and Pin values (droop curve, [20]) for real regulators.
For a constant P0, Pin decreases as the filter chokes. When the weekly cumulative gas flow (WC-Q) within a network remains nearly the same, then Pout also decreases to compensate for the lower Pin in order to maintain the same WC-Q. Hence, filter choking impacts Pout. However, the model with no filter choking predicted a higher WC-Q value for the measured pressures. Thus, an increase in WC-Q could indicate filter choking. When Pin values were not available, then P0 values could be used as proxies, and the flow could be estimated using our mathematical model. As shown in Figure S1 (Supplementary Materials), the impact of Pout was more significant than that of Pin on the calculated flow. Hence, the inferences that were drawn were not compromised. Thus, WC-Q was estimated from the data using our simplified mathematical model. The filter RUL was estimated by fitting a line through the estimated WC-Q values. Filter maintenance is recommended when the RUL values exceed an agreed threshold value. The monitoring steps that were involved in our model were as follows:
  • Input threshold WC-Q values for the PRs (4″, 3″, 2″, and 2″ SR);
  • Calculate gas flows (Q) for all data points using the simplified mathematical model (as discussed in Section 2 and Section 3) Using the model parameters for the 4″, 3″, 2″, and 2″ SR PRs that are presented in Table 1;
  • In the absence of Pin data, QM,2P = f (PR type, model parameters, P0, and Pout);
  • When Pin is available, QM,3P = f (PR type, model parameters, Pin, and Pout);
  • PR type: spring-loaded PR or lever-type PR;
  • Calculate WC-Q;
  • Plot WC-Q against time and obtain the best linear fit;
  • From the linear fit, estimate the number of weeks for the current WC-Q to reach the threshold WC-Q. That was the estimate for the filter RUL.

5.3. Valve Seat Damage

The valve seats make repeated contact with the orifices during the operation of PRs. Contaminant particles can also be deposited from various sources, including pipe wear and tear. Damaged valve seats can allow tiny leaks, even when the valves are completely shut. This is known to produce a creep effect. This creep effect manifests as a downstream over-pressure in the network during low demand periods (e.g., around midnight).
After analyzing the available data, we confirmed that low flows could exist in the PRs between 12 am and 3 am. Following the work of Leo et al. [1], the available Pout data for the gas network during this period were checked carefully and were found to be too noisy to reliably determine the daily maximum Pout value. Hence, to minimize the impacts of abrupt fluctuations and sensor noise, the average Pout value was computed for that period. The average Pout value seemed more reliable than the daily Pout value for monitoring valve seat condition.
  • Input threshold Pout value for creep effect (valve seat damage);
  • For each day, compute the average Pout value from 12 am to 3 am;
  • Plot the average Pout against time and obtain the best linear fit;
  • Using the linear fit, estimate the number of days needed for the current Pout value to reach the threshold Pout value. That was the estimate for the valve seat RUL.

5.4. Diaphragm Malfunction

The diaphragm and valve seat move together to maintain Pout near the set point. When the diaphragm senses a deviation in Pout from the set point, the valve seat moves up or down accordingly. The diaphragm is thus subjected to cyclic stress, which can eventually lead to its fatigue failure. The lifetime of diaphragms is measured using up-down cycles (200,000; [21]). The up-down movement of diaphragm is measured using the spring extension, but it is very difficult to measure the spring extensions in hundreds of regulators.
After closely inspecting the real operational data, we noticed that Pout experienced more fluctuations in small flows around midnight. There can be multiple up-down cycles during each hour of PR operation. Hence, a methodology to count all of the up-down cycles in Pout was developed.
  • Detrend Pout data for each hour (Figure S2 in the Supplementary Materials). For this, the linear best fit was subtracted from the raw data to produce a zero mean for the detrended data;
  • Find the local minimum and maximum values (Figure S2 in the Supplementary Materials). The local minimum and maximum values were used to identify the data points at which the diaphragm movement changed from upward to downward or from downward to upward;
  • Construct up-down cycles using the local minimum and maximum values (Figure 6a);
  • Determine the number of cycles per day using (maximum and minimum values − 1)/3 (Figure 6b);
  • Compute percentage deterioration in diaphragm condition over a year using Equation (17).
Annual   Use   of   a   Diaphragm = Number   of   cycles   in   a   year Maximum   number   of   cycles %

6. Results and Discussion

6.1. Pressure Regulator Health Monitor (PRHM) Dashboard

A user-friendly dashboard (the PRHM) was developed to monitor the health of PRs within gas networks. The hundreds of PRs within a network can be divided into groups based on their geographical locations for easier analysis and comparison. Through the PRHM, PR groups can be updated by adding or removing PRs.
The PRHM has two tabs: “Group Analysis” and “Performance Comparison”. The first tab is used to evaluate some or all PRs within a PR group. The second tab is used to analyze and compare the performances of the PRs in different groups. For each PR, data files are updated automatically as needed. The user can choose the start and end dates for the analysis. The GUI of the PRHM is shown in Figure S3 in the Supplementary Materials. The left-hand side of the PRHM (“Working Directory”, “Add New Data”, “Arrange/Edit”, “Analyze Group”, and “Plots Display”) indicates its capabilities. When a particular regulator group is chosen for analysis (e.g., the PR Group 2 in Figure S3 in the Supplementary Materials), then the right-hand side table in the PRHM presents all of the essential information for that group. The data in the table can be modified as needed before starting the analysis. When a regulator is chosen for analysis, then the color of the entire row changes (e.g., in Figure S3 in the Supplementary Materials, in which a color change to yellow is shown).
In this study, the PR Group 2 contained 15 PRs: RG2_1, RG2_2, …, RG2_15. Each PR had data for different dates, as shown in Figure S3 in the Supplementary Materials. In total, 10 of the 15 PRs were chosen for analysis and then the “Analyze Group” button was clicked. In total, there were 404,694 data points for the 10 PRs. The PRHM dashboard showed the progress of calculations and error messages (when applicable). In the case of error messages, it provided suggestions for how to resolve the errors. The PRHM took about 7 minutes to analyze the 10 PRs using a Windows workstation with 16 GB of RAM and a i7 7700 CPU at 3.60 GHz. After completing the analysis, the PRHM showed three tables for the filter RUL, the creep effect (valve seat) RUL, and the annual diaphragm usage (Figure 7). Users could choose a PR from the drop-down list under the “Plots Display” menu to visualize the pressures and flows. In Figure 7, the inlet pressure, intermediate pressure, outlet pressure, and gas flow data are shown for RG2_5. The sliders below the plots could be used to adjust their display time windows.
All of the results from the analyses were automatically saved in different folders. Old results were removed after one month. Each regulator group had a separate folder for their results (e.g., “PR Group 2″ was the name of the folder for PR Group 2). Each group’s folder had five subfolders for the various plots: (1) the daily number of data points for each PR; (2) the inlet, intermediate, and outlet pressure data and the pressure difference across each filter; (3) the daily average flow through each PR; (4) the weekly cumulative flow through each PR; and (5) the daily average Pout from 12 am to 3 am for each PR.

6.2. Filter Choking

As explained previously, the weekly cumulative flow through a PR is used to monitor filter choking. Figure 8 presents the plots for the weekly cumulative flow against time for eight selected PRs (RG2_1, RG2_3, RG2_5, RG2_7, RG2_8, RG2_9, RG2_10, and RG2_13). For these PRs, the weekly cumulative flow increased monotonically.
The linear regression line is shown in each plot, along with the slope and intercept. Table S2 in the Supplementary Materials shows the confidence intervals for the intercept and slope of the fitted lines. The confidence intervals were calculated for fitted lines with positive slopes and their parameters were statistically significant. The threshold value of the weekly cumulative flow that was used to calculate the RUL was 20% above the average weekly flow that was predicted by the model for the first week. The resultant RUL values are shown in Figure 8 and Table 2. For RG2_2 and RG2_6, the linear regression lines had negative slopes due to the poor quality of the pressure data. In these cases, the dashboard did not display any filter RUL estimations.

6.3. Valve Seat Damage (Creep Effect)

Figure 9 presents the average Pout values (12 am to 3 am) for several days. The average daily Pout value only increased for RG2_7, RG2_8, RG2_9, and RG2_10. The Pout values are expected to increase for damaged valve seats; hence, it could be concluded that the other PRs had no issues with their valve seats. The fitted lines for various PRs are shown in Figure 9. Table S3 in the Supplementary Materials shows the confidence intervals for the intercept and slope of the lines. The confidence intervals were only computed for lines with positive slopes and the line parameters were statistically significant (except RG2_10, which had a lower confidence interval because the slope was negative). The threshold value of the average daily Pout was 2.5 kPa. The RUL for each valve seat can be seen in Figure 9 and Table 2.

6.4. Diaphragm Malfunction

For the selected PRs in PR Group 2, Figure S5 in the Supplementary Materials shows the daily diaphragm cycles. In general, diaphragms have a lifetime of about 10 years or 200,000 cycles [21]. The annual diaphragm usage (i.e., the number of diaphragm cycles over one year expressed as a percentage of the number of cycles over the lifetime of the diaphragm) of the PRs are presented in Table 2. The 10 diaphragms completed around 11.8% of their lifetime cycles over one year. This correlated well with the expected lifetime of 10 years, but earlier replacement could be required (depending on when the diaphragm was last changed).

7. Conclusions

In this work, mathematical models for spring-loaded PRs and lever-type PRs were developed. As gas distribution networks can have hundreds of PRs, the mathematical models were simplified for faster computation. Methodologies for monitoring filter choking, valve seat damage, and diaphragm deterioration were developed for PRs, which demonstrated high reliability and accuracy. The methodologies combined the mathematical models with the available experimental pressure data that were measured upstream and downstream of regulators to infer the health of various components within PRSs. It was demonstrated that filter choking, which occurs relatively rapidly, could be monitored by tracking increases in estimated weekly cumulative flows through the PRs. Similarly, the development of valve seat damage in the regulators, which is a slow process, could be assessed over a long period of time by observing creeping increases in average daily Pout values for a period in the day when no flow was expected. Diaphragm condition could be monitored by the number of up-down cycles using the Pout values. The methodologies for monitoring these faults were extended to estimate the RUL of each component. For industrial users, a user-friendly PRHM dashboard was created to provide a one-stop platform for monitoring the health of PRs within gas networks and planning the timely maintenance of deteriorating components. This work could help toward increasing the system reliability of legacy systems and ensuring uninterrupted gas supplies for customers.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/en15176264/s1, Figure S1: Impact of filter choking on the flow calculated by the mathematical model; Figure S2: (a) detrend Pout, (b) local maximum ●, (c) local minimum ■; Figure S3: Input data table to analyze a group (e.g., regulator group 2) of pressure regulators; Figure S4: Daily number of data samples for RG2_1, RG2_2, RG2_3, RG2_9, RG2_10 and RG2_13 pressure regulators; Figure S5: Daily number of diaphragm cycles for different pressure regulators; Table S1: Vendor data (inlet and outlet pressures, and flow) for all pressure regulators. The capacities shown in the tables are given in terms of natural gas specific gravity (SG = 0.6). For all other gases, multiply flow by the correction factor (= 0.6/SG of gas)0.5; Table S2: Confidence intervals of estimated coefficients (intercept and slope) of linear regression model for filter (RUL). The confidence intervals were calculated for fitted line with positive slope; Table S3: Confidence intervals of estimated coefficients (intercept and slope) of linear regression model for valve seat (RUL). The confidence intervals were calculated for fitted line with positive slope.

Author Contributions

Conceptualization, S.S., I.A.K., S.F., L.S. and R.S.; Formal analysis, S.S.; Methodology, S.S.; Supervision, I.A.K., S.F., L.S. and R.S.; Writing–original draft, S.S.; Writing–review & editing, I.A.K., S.F., L.S. and R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Energy Innovation Research Program (EIRP) Award NRF2017EWT-EP003-020 from the Energy Market Authority (EMA) of Singapore. The EIRP is a competitive grant call initiative that is driven by the Energy Innovation Program Office and funded by the National Research Foundation (NRF) of Singapore.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

List of Symbols

ADArea of diaphragm (m2)
AOArea of plunger (m2)
ApPeripheral flow area (m2)
CdDischarge coefficient of orifice
dDiameter of orifice (m)
DDiameter of diaphragm (m)
doDiameter of plunger (m)
dtDiameter of shaft (m)
F0Initial force on spring (N)
KSpring constant (N/m)
L1, L2, MLever dimensions (m)
mout Outlet mass flow (kg/s)
PatmAtmospheric pressure (kPa)
P0Inlet pressure (kPa)
PinIntermediate pressure (kPa)
PoutOutlet pressure (kPa)
QFlow (m3/h)
RGas constant (J/(K.mol))
TinInlet temperature (K)
Tout Outlet temperature (K)
WC-QWeekly cumulative flow (m3/h)
xSpring extension (m)
XGas composition
yPlunger displacement (m)
ΔPfPressure difference across filter (kPa)
μJoule–Thomson expansion coefficient
ρDensity of gas (kg/m3)
γHeat capacity ratio of gas
θ1, θ2, θdAngles of lever
GUIGraphical user interface
HONHoneywell
MModel
NNNeural network
OFObjective function
PPressure
PRPressure regulator
PRHMPressure regulator health monitor
PRS Pressure regulating station
RMGRegel and Messtechnik GmbH
RULRemaining useful life
SService regulator
SNNumber of samples
SPSet point
VVendor

References

  1. Leo, M.B.; Dutta, A.; Farooq, S.; Karimi, I.A. Simulation and health monitoring of a pressure regulating station. Comput. Chem. Eng. 2020, 139, 106824. [Google Scholar] [CrossRef]
  2. Rao, H.N.; Xiang, L.R.; Karimi, I.A. Modelling pressure regulators using operational data. In Foundations of Computer Aided Process Operations/Chemical Process Control; San Antonio, Texas. 2017. Available online: https://folk.ntnu.no/skoge/prost/proceedings/focapo-cpc-2017/FOCAPO-CPC%202017%20Contributed%20Papers/106_FOCAPO_Contributed.pdf (accessed on 2 April 2021).
  3. Afshari, H.H.; Zanj, A.; Novinzadeh, A.B. Dynamic analysis of a nonlinear pressure regulator using bondgraph simulation technique. Simul. Model. Pract. Theory 2010, 18, 240–252. [Google Scholar] [CrossRef]
  4. Dasgupta, K.; Karmakar, R. Modelling and dynamics of single-stage pressure relief valve with directional damping. Simul. Model. Pract. Theory 2002, 10, 51–67. [Google Scholar] [CrossRef]
  5. Nabi, A.; Dayan, J. Dynamic model for a dome-loaded pressure regulator. J. Dyn. Syst. Meas. Control. Trans. ASME 2000, 122, 290–297. [Google Scholar] [CrossRef]
  6. Rami, E.G.; Jean-Jacques, B.; Pascal, G.; François, M. Stability study and modelling of a pressure regulating station. Int. J. Press. Vessel. Pip. 2005, 82, 51–60. [Google Scholar] [CrossRef]
  7. Zafer, N.; Luecke, G.R. Stability of gas pressure regulators. Appl. Math. Model. 2008, 32, 61–82. [Google Scholar] [CrossRef]
  8. Wang, X.S.; Cheng, Y.H.; Peng, G.Z. Modeling and self-tuning pressure regulator design for pneumatic pressure load systems. Control Eng. Pract. 2007, 15, 1161–1168. [Google Scholar]
  9. Prescott, S.L.; Ulanicki, B. Dynamic modeling of pressure reducing valves. J. Hydraul. Eng. 2003, 129, 804–812. [Google Scholar]
  10. Shahani, A.R.; Esmaili, H.; Aryaei, A.; Mohammadi, S.; Najar, M. Dynamic simulation of a high-pressure regulator. JCARME 2011, 1, 17–28. [Google Scholar]
  11. Ramzan, M.; Maqsood, A. Dynamic modeling and analysis of a high pressure regulator. Math. Probl. Eng. 2016, 2016, 1307181. [Google Scholar] [CrossRef]
  12. Patil MGand Barjibhe, R.B. Flow analysis of gas pressure regulator by numerical and experimental method. Int. Res. J. Eng. Technol. (IRJET) 2018, 5, 747–752. [Google Scholar]
  13. Yinxing, D.; Huaixiu, W.; Yahui, W.; Fangwen, C. Research on fault diagnosis algorithm of gas pressure regulator based on compressed sensing theory. In Proceedings of the 2017 29th Chinese Control and Decision Conference, Chongqing, China, 28–30 May 2017. [Google Scholar]
  14. Yun, A.; Yahui, W.; Yuexiao, L. Research on gas pressure regulator fault diagnosis based on deep confidence network (DBN) Theory. In Proceedings of the 2017 Chinese Automation Congress, Jinan, China, 20–22 October 2017. [Google Scholar]
  15. Song, Y.; Wang, Y. A fault detection method for gas pressure regulators based on improved dynamic canonical correlation analysis. In Proceedings of the 2021 33rd Chinese Control and Decision Conference, Kunming, China, 22–24 May 2021. [Google Scholar]
  16. Jia, J.; Wang, B.; Ma, R.; Deng, Z.; Fu, M. State monitoring of gas regulator station based on feature selection of improved grey relational analysis. IEEE Internet Things J. 2022. [Google Scholar] [CrossRef]
  17. Wang, B.; Jia, J.; Deng, Z.; Fu, M. A state monitoring method of gas regulator station based on evidence theory driven by time-domain information. IEEE Trans. Ind. Electron. 2022, 69, 694–702. [Google Scholar] [CrossRef]
  18. RMG by Honeywell. Gas Pressure Regulator Technical Datasheet. Available online: https://www.honeywellprocess.com/library/marketing/brochures/HON_280-284%20Brochure.pdf (accessed on 2 April 2021).
  19. Andersen, B.W. The Analysis and Design of Pneumatic Systems; Wiley: New York, NY, USA, 1967. [Google Scholar]
  20. Emerson. The Industry Standard for Pressure Regulators. Available online: https://www.emerson.com/documents/automation/catalog-regulators-mini-catalog-fisher-en-125484.pdf (accessed on 15 October 2019).
  21. Sutar, S.S.; Kunbade, P.R.; Jamadade, S.S. Fatigue life estimation of pressure reducing valve diaphragm. IJEAT 2015, 4, 180–188. [Google Scholar]
Figure 1. An overview of a PRS (gas filter and PR) with a downstream network.
Figure 1. An overview of a PRS (gas filter and PR) with a downstream network.
Energies 15 06264 g001
Figure 2. (a) A simple schematic of the HON 280 gas PR and (b) the force balance on the diaphragm and valve seat assemblies.
Figure 2. (a) A simple schematic of the HON 280 gas PR and (b) the force balance on the diaphragm and valve seat assemblies.
Energies 15 06264 g002
Figure 3. (a) A schematic of the HON 277 gas PR, (b) the torque balance for the lever-type PR (HON 277), and (c) the computation of the displacement of the lever and flow area.
Figure 3. (a) A schematic of the HON 277 gas PR, (b) the torque balance for the lever-type PR (HON 277), and (c) the computation of the displacement of the lever and flow area.
Energies 15 06264 g003
Figure 4. The optimization-based methodology that was used to estimate the K values using vendor data.
Figure 4. The optimization-based methodology that was used to estimate the K values using vendor data.
Energies 15 06264 g004
Figure 5. An overview of the developed PR fault monitoring methodologies.
Figure 5. An overview of the developed PR fault monitoring methodologies.
Energies 15 06264 g005
Figure 6. (a) The number of cycles in one hour (1, 2, …21 are the local minimum and maximum values) and (b) the daily number of cycles.
Figure 6. (a) The number of cycles in one hour (1, 2, …21 are the local minimum and maximum values) and (b) the daily number of cycles.
Energies 15 06264 g006
Figure 7. The GUI for the PR health monitor dashboard showing the estimated RUL values for the filter and valve seat, the annual diaphragm usage, and the pressure and flow plots.
Figure 7. The GUI for the PR health monitor dashboard showing the estimated RUL values for the filter and valve seat, the annual diaphragm usage, and the pressure and flow plots.
Energies 15 06264 g007
Figure 8. The weekly cumulative flow and filter RUL estimation for different PRs.
Figure 8. The weekly cumulative flow and filter RUL estimation for different PRs.
Energies 15 06264 g008
Figure 9. The average daily Pout values (midnight to 3 am) and the estimation of valve seat RUL. No creep effect was observed for the other PRs.
Figure 9. The average daily Pout values (midnight to 3 am) and the estimation of valve seat RUL. No creep effect was observed for the other PRs.
Energies 15 06264 g009aEnergies 15 06264 g009b
Table 1. The dimensions and model parameters for the two PRs.
Table 1. The dimensions and model parameters for the two PRs.
Spring-Loaded PRLever-Type PR
4″3″2″ 2″ SR
dt (m)0.0190.01410.0145L1 (m)0.026
d (m)0.09980.0630.0497L2 (m)0.095
do (m)0.1120.0880.06M (m)0.042
D (m)0.320.2960.272d (m)0.032
D (m)0.239
Cd0.870.870.87Cd0.87
Pout,SP (kPa)151515Pout,SP (kPa)17
F0 (N)120.6103.287.2F0 (N)76.3
K (N/m)51,00038,80026,200K (N/m)3,790
Table 2. The estimated RUL values for the filter and the valve seat and the annual diaphragm usage for the different PRs.
Table 2. The estimated RUL values for the filter and the valve seat and the annual diaphragm usage for the different PRs.
PR NameFilter RULValve Seat RULAnnual Diaphragm Usage
(d)(d)(%)
RG2_1491013.5
RG2_20012.5
RG2_3399013.1
RG2_5413013.2
RG2_6009.4
RG2_7337234412.7
RG2_832561210.8
RG2_9491453612.9
RG2_10500863810.2
RG2_1337909.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sharma, S.; Karimi, I.A.; Farooq, S.; Samavedham, L.; Srinivasan, R. Health Monitoring of Pressure Regulating Stations in Gas Distribution Networks Using Mathematical Models. Energies 2022, 15, 6264. https://doi.org/10.3390/en15176264

AMA Style

Sharma S, Karimi IA, Farooq S, Samavedham L, Srinivasan R. Health Monitoring of Pressure Regulating Stations in Gas Distribution Networks Using Mathematical Models. Energies. 2022; 15(17):6264. https://doi.org/10.3390/en15176264

Chicago/Turabian Style

Sharma, Shivom, Iftekhar A. Karimi, Shamsuzzaman Farooq, Lakshminarayanan Samavedham, and Rajagopalan Srinivasan. 2022. "Health Monitoring of Pressure Regulating Stations in Gas Distribution Networks Using Mathematical Models" Energies 15, no. 17: 6264. https://doi.org/10.3390/en15176264

APA Style

Sharma, S., Karimi, I. A., Farooq, S., Samavedham, L., & Srinivasan, R. (2022). Health Monitoring of Pressure Regulating Stations in Gas Distribution Networks Using Mathematical Models. Energies, 15(17), 6264. https://doi.org/10.3390/en15176264

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop