Multiscale modelling of hepatitis B virus at cell level of organization

1.   Introduction

Load frequency control (LFC) is significant for the safe and stable work of the power grids ^[1,2,3]. With the continuous breakthrough of new energy technology, new energy power generation methods are gradually combined with traditional power generation methods and the penetration rate of modern smart grids continues to increase. Wind power and photovoltaic power are the most widely used forms of new energy generation. However, wind energy and photovoltaic energy are characterized by dispersion, randomness, fluctuation and intermittency ^[4], which can affect the power grid, especially the frequency.

Usually, when the frequency of the power system fluctuates, the aggregator load changes its working state according to the scheduling instruction to prevent frequency fluctuation ^[5,6]. Most dynamic demand control models have been with thermostatically controlled devices, with air conditioning loads (ACs) accounting for a huge percentage of them ^[7,8]. However, most literature focused on how ACs can be controlled based on optimized thought, such as the method based on the load frequency modulation model ^[9]. Only a few studies have applied ACs to resolve external perturbations ^[10,11].

With the progress of ultra-high voltage technology and the increasing proportion of new energy in power generation ^[12,13,14], the communication data between power grid systems has become more complex, which also increases the communication pressure between power grids. In recent years, there has been a growing interest in event-triggered mechanisms (ETM) due to their potential in reducing communication and computation resources while maintaining stability and performance in complex systems ^{[15,16,17,18,19]}. ETM is designed to trigger actions and update processes in a system only when certain predefined conditions, or events, are met rather than at fixed time intervals. This smart adaptive strategy is particularly beneficial for systems in which resources are constrained or in dynamic environments where unnecessary updates could be costly. Kazemy et al. ^[15] have addressed synchronization of master-slave neural networks under deception attacks using event-triggered output feedback control to enhance stability and resilience.

While considering the limited network resources, we should also consider the security problems faced by networked systems ^[20]. Especially in recent years, industrial open networks have been attacked frequently. At present, network attacks are mainly divided into denial-of-service (DoS) attacks and spoofing attacks. DoS attacks can cause the loss of transmitted data ^[21], while spoofing attacks are launched by attackers that attempt to disrupt the availability of data ^[22]. Foroush and Martinez ^[23] showed that the general security problem of networked is DoS attack, which attempts to transmit a large amount of invalid data and intentionally interfere with network communication resources. Persis and Tesi ^[24] proved that on the premise of ensuring the asymptotic stability of the system, and the average active time of DoS attacks cannot exceed a certain percentage. In recent years, the research on the security of event-triggered networked systems has made some progress ^[15].

Inspired by the analysis presented above, this paper explores the application of an observer-based ETM for secure LFC in power systems, with the aim of ensuring the safety and stability of the power system. Considering the inherent instability caused by the increasing penetration of renewable energy in the grid, we treat it as a perturbation. In addition, this study proposes the integration of AC demand response as a control strategy to effectively regulate system frequency fluctuations. Additionally, this paper addresses key challenges related to network bandwidth limitations and potential DoS attacks by incorporating the ETM. A novel aspect of this work is the application of ETM to the observer rather than directly to the system state to mitigate the complexity and impracticalities associated with direct state measurement in real-world applications. The contribution of this paper is as follows:

ⅰ) The volatility associated with renewable energy generation is modeled as a perturbation, accounting for its unpredictable nature in the power system. To enhance the LFC, ACs demand response is incorporated as an effective tool to dynamically adjust the system's frequency in response to fluctuating generation conditions, thereby improving overall system resilience.

ⅱ) In contrast to conventional ETM that rely on direct measurement of system states, this paper introduces an observer-based ETM. By using an observer to estimate system states, this approach avoids the challenges and limitations of directly measuring these values, offering a practical solution for real-time control in large-scale, decentralized power networks.

Notations. $U > 0$ means the matrix $U$ is positive definite. $*$ stands for the corresponding symmetric term in the matrix. $diag\{\cdots\}$ stands for diagonal matrix.

2.   Problem formulation

2.1.   Model description

First, considering the random disturbance of wind power and photovoltaic power generation, a multi-area smart grid control model is established with ACs participating in grid frequency regulation, as shown in Figure 1 ^[25]. System parameters and meanings are shown in Table 1, where $k_{aci} = \frac{m_{i}c_{pi}k}{\mathrm{EER}}$.

According to Figure 1, we can obtain

where $D_{aci} = \frac{2\pi d_{aci}}{M_{i}}$, $T_{ii} = \sum_{j = 1, j\neq i}^{n}T_{ij}, i, j = 1, 2, \ldots, n$. Define

Then, we can obtain that

where

From (2.2), the state-space form of multi-area ACs LFC model can be described as

with

2.2.   Full-order state observer

The full-order state observer is constructed as

where $\hat{\mathrm{y}}(\mathfrak{t})$ is the observer output; $\hat{\mathrm{x}}(\mathfrak{t})$ is the observer state and $\mathcal{L}$ is the observer gain to be solved.

2.3.   Observer-based ETM under DoS attacks

In this paper, inspired by ^[26,27], and periodic event-triggered control schemes ^[28], we select the ETM activations:

where the difference between the most recent data sample and the current one is referred to as $\zeta(\mathfrak{t})$, i.e., $\zeta(\mathfrak{t}) = \hat{\mathrm{x}}(\mathfrak{t}_k+lh)-\hat{\mathrm{x}}(\mathfrak{t}_k)$, $lh$ is the intervals between $\mathfrak{t}_{k}$ and $\mathfrak{t}_{k+1}$, $l\in N, i_kh\in (\mathfrak{t}_k, \mathfrak{t}_{k+1}]$, and $\sigma$ is a threshold parameter.

Remark 1. Unlike traditional ETM that rely on state controllers, this paper introduces a novel approach using observer packets for triggering. This approach leverages the observer pattern, which offers significant advantages in terms of modularity and flexibility. Specifically, the observer pattern enables the seamless addition of new observers to the system without disrupting or modifying other components, ensuring a high degree of decoupling. Similarly, the observed object itself is free to evolve or modify its internal implementation without causing any adverse effects on the observers.

In addition, consider the time-varying delay $d_k$ caused by the network between the observer and the event generator ^[29]. Define $d(\mathfrak{t})$ as follows:

where $I_{1} = [\mathfrak{t}_{k}+d_k, \mathfrak{t}_{k}+d_k+h)$, $I_{2}^{l} = [\mathfrak{t}_{k}+d_k+lh, \mathfrak{t}_{k}+d_k+lh+h)$, $I_{2} = \cup_{l = 1}^{l = m-1}I_{2}^{l}$, $I_{3} = [\mathfrak{t}_{k}+d_k+mh, \mathfrak{t}_{k+1}+d_{k+1})$, $[\mathfrak{t}_{k}+d_k, \mathfrak{t}_{k+1}+d_{k+1}) = I_{1}\cup I_{2}\cup I_{3}$. Moreover, $\hat{\mathrm{x}}(\mathfrak{t}_k)$ and $\hat{\mathrm{x}}(\mathfrak{t}_k+\iota h)$ with $\iota = 1, 2, \dots, m$ satisfy (2.5); $d_m$ ($d_M$) denotes the lower (upper) bound on the time delay.

Remark 2. In a multi-area power system, data transmission is routine, but network delays are inevitable. These delays can arise from factors like communication latency or network congestion, and if ignored, they may compromise the accuracy of the measured data. In this study, we define the delay as $d_k$ and focus on the case where $d_k$ is shorter than the sampling interval. This ensures that the integrity of the data sequence is preserved, with each measurement arriving and being processed within its designated sampling period.

For $t\in [\mathfrak{t}_{k}+d_k, \mathfrak{t}_{k+1}+d_{k+1})$, we can define

In what follows, inspired by ^[30], the observer-based ETM $\mathrm{u}(\mathfrak{t})$ can be written as

where $\mathcal{K}$ is the control gain matrix.

Based on (2.6) and (2.7), the observer-based ETM $\mathrm{u}(\mathfrak{t})$ can be rewritten as

For a DOS attack in a network, define $\tau(i_{k}h)$ as the attack state, $\tau(i_{k}h) = 1$ as the attack occurred, and $\tau(i_{k}h) = 0$ as the attack did not occur. In a typical power system, most DoS attacks fail to achieve their purpose, and only a fraction of these attacks cause notable disruptions. We divide DoS attacks into the following four types:

In addition, the energy of general DoS attacks is limited. Facing such attacks, we introduce elastic variable $\zeta_{DoS}(\mathfrak{t})$ on the event triggering mechanism to represent the packet loss and delay caused by DoS attacks on the system ^[31]. We can obtain

where $\Delta_{DoS}$ denotes the longest duration of a DoS attack and $\mathfrak{t}_{k+1}^{DoS}$ refers to an upcoming sampling moment affected by DoS.

2.4.   Modeling augmented systems

Substituting (2.8) into sysetm (2.3) yields

Substituting (2.8) into system (2.4) yields

Define the state error as $\delta(\mathfrak{t}) = \mathrm{x}(\mathfrak{t})-\hat{\mathrm{x}}(\mathfrak{t})$, then

Define the augmented state as $\tilde{\mathrm{x}}^T(\mathfrak{t}) = \left[\begin{matrix}\mathrm{x}^T(\mathfrak{t}) & \delta^T(\mathfrak{t})\end{matrix}\right]$, the augmented system is

where

Definition 1 ^[32,33] For an arbitrary initial condition, given $\gamma > 0$, the augmented system (2.15) is thought to be stable with $H_{\infty}$ performance if the following two conditions are met:

1) for $\omega(\mathfrak{t}) = 0$, the system is asymptotically stable;

2) for $\omega(\mathfrak{t})\in\ell_{2}[0, \infty)$, the following condition is satisfied:

Lemma 1. ^[34] For a given time delays constant $d_m$ and $d_M$ satisfying $d(\mathfrak{t})\in[d_{m}, d_{M}]$, for $d_{M} > d_{m} > 0$, there exist matrices $\mathsf{R}_2$ and $\mathsf{H}$ when

for which the following inequality holds:

3.   Main results

Two theorems are given in this section. Theorem 1 illustrates the stability and $H_{\infty}$ performance of the augmented system, and Theorem 2 gives the design of an observer-based controller.

Theorem 1. For some given positive constants $\sigma$, $d_m$ and $d_M$, the augmented system (2.15) is asymptotically stable with an $H_{\infty}$ performance index $\gamma$, if there exist positive definite matrices $\mathsf{P}$, $\mathsf{Q}_1$, $\mathsf{Q}_2$, $\mathsf{R}_1$, $\mathsf{R}_2$, $\mathsf{H}$ and $\mathsf{V}$ with appropriate dimensions such that

where

Proof. Select Lyapunov function as

From the augmented system (2.15), we have

Using Jessen inequality, Lemma 1 and the Schur complement theorem, at the same time to join the event trigger condition, it can be concluded that

where

Further, we achieve

According to (3.4) and (3.5), the presence of $\lambda > 0$ ensures that $\dot{\mathbb{V}}(\mathfrak{t}) < 0$ is satisfied when $\omega(\mathfrak{t}) = 0$, and we can achieve $lim_{t\to\infty}\mathbb{V}(\mathfrak{t}) = 0$.

By Definition 1, the augmented system (2.15) is asymptotically stable.

When $\omega(\mathfrak{t})\neq0$, both sides of (3.5) integrate from $0$ to $+\infty$,

At zero original condition, we obtain

Since $\Omega < 0$, there exists a proper positive $\varepsilon$ so that

Multiply the upper expression by $e^{\varepsilon \mathfrak{t}}$ and integrate the result. Subsequently,

It is evident that

In the presence of an attack, the system's security performance bounds uniformly, with a degradation no greater than $\Delta$. This finishes the proof.

Theorem 2. For some given positive constants $\sigma$, $d_m$ and $d_M$, the augmented system (2.15) is asymptotical stable with an $H_{\infty}$ performance index $\gamma$, if there exist positive definite matrices $\mathsf{X}_{i}$, $\mathsf{Y}$, $\mathsf{S}$, $\overline{\mathsf{Q}}_{i}$, $\overline{\mathsf{R}}_{i}$, $\overline{\mathsf{H}}$ and $\overline{\mathsf{V}}$ $(i = 1, 2)$ with appropriate dimensions such that

where

The control coefficient matrix of the observer-based controller is shown by: $\mathcal{K} = \mathsf{Y}\mathsf{X_2}^{-1}$, $\mathcal{L} = \mathsf{S}\overline{\mathsf{X}}_2^{-1}$, with $\mathcal{C}\mathsf{X_2} = \overline{\mathsf{X}}_2\mathcal{C}$.

Proof. Set $\mathsf{P} = \begin{bmatrix}\mathsf{P_1} & {\bf 0}\\ * & \mathsf{P_2}\end{bmatrix}$ and define $\mathsf{X} = \mathsf{P}^{-1} = \begin{bmatrix}\mathsf{X_1} & {\bf 0}\\ * & \mathsf{X_2}\end{bmatrix}$, $\overline{\mathsf{Q}}_i = \mathsf{X}^T\mathsf{Q}_i\mathsf{X}$, $\overline{\mathsf{R}}_i = \mathsf{X}^T\mathsf{R}_i\mathsf{X}$ $(i = 1, 2)$, $\overline{\mathsf{H}} = \mathsf{X}^T\mathsf{H}\mathsf{X}$, $\overline{\mathsf{V}} = \mathsf{X_2}^T\mathsf{V}\mathsf{X_2}$, $\mathsf{Y} = \mathcal{K}\mathsf{X_2}$, $\mathsf{S} = \mathcal{L}\overline{\mathsf{X}}_2$. For $\mathsf{X_2} = \mathsf{W}\begin{bmatrix}\mathsf{X_{21}} & {\bf 0}\\ * & \mathsf{X_{22}}\end{bmatrix}\mathsf{W}^T$, on the basis of ^[35], there exist $\overline{\mathsf{X}}_2 = \mathsf{M}\mathsf{N}\mathsf{X_{22}}\mathsf{N}^{-1}\mathsf{M}^T$ such that $\mathcal{C}\mathsf{X_2} = \overline{\mathsf{X}}_2\mathcal{C}$.

Define $\Phi = diag\{\mathsf{X}, \mathsf{X}, \mathsf{X}, \mathsf{X}, \mathsf{X_2}, I, I, I, I, I\}$, then pre- and postmultiply $\Phi$ and its transpose on both sides of (3.1). It is easy to know that $-\mathsf{X}\overline{\mathsf{R}}_{i}^{-1}\mathsf{X}\leq\overline{\mathsf{R}}_{i}-2\mathsf{X}, (i = 1, 2)$, $-\mathsf{X_2}\overline{\mathsf{V}}^{-1}\mathsf{X_2}\leq\overline{\mathsf{V}}-2\mathsf{X_2}$. Finally, we can conclude that (3.1) is a sufficient condition to guarantee (3.6) holds. The proof is accomplished.

Remark 3. The power system under investigation is a complex continuous system, whose stabilization is further complicated by disturbances such as wind power generation, making the calculation of controller parameters a challenging task. Lyapunov theory is well used to facilitate analysis and predict system behavior, especially for complex systems.

4.   Simulation results

In this part, we consider a simulation example for a two-area power system and demonstrate the usefulness of this method. The parameter values of system (2.15) are shown in Table 2 ^[25].

Set $\sigma = 0.01$, $d_{M} = 0.1$, $d_{m} = 0.02$ and the sampling period is $h = 0.05$. The initial condition, in this simulation analysis, is given by $\tilde{x}(0) = \left[\begin{matrix}0.001 & 0.001 & \ldots & 0.001 \end{matrix} \right]$. Next, by solving Theorem 2, the control gain $\mathcal{K}$ and the observer gain $\mathcal{L}$ are given by

Meanwhile, $H_{\infty}$ performance index $\gamma = 15.4285$. Here we set the disturbance $\omega(t)$ occurs at every sampling instant as $\omega(t) = \frac{\sqrt{0.16}}{1+t^2}$, and the probability that the system suffers from an effective DoS is $0.3$.

In this simulation, we analyze the effects of a DoS attack on system dynamics and control responses.

Figure 2 illustrates DoS attack signals from $0-10\; s$ denoted by $\tau(i_k h)$, highlighting the attack's temporal characteristics and intervals of disrupted transmission.

The frequency deviation $\Delta f_i$ in the affected region, as depicted in Figure 3, shows fluctuations due to the DoS attack. However, it gradually stabilizes and tends toward zero, indicating a system tendency to regain frequency balance over time. This response underscores the system's capacity for frequency recovery even under compromised conditions.

Figure 4 presents the output power deviation $\Delta P_{aci}$ of the ACs in the region under DoS attack. Notably, the fluctuations are contained within a narrow range, remaining below 0.1, demonstrating a limited impact on the output power deviation of the ACs. This outcome suggests that the control strategy maintains a relatively stable response in the face of attack-induced disturbances.

The area control error $ACE_i$ in Figure 5 provides insights into area control performance under the DoS attack.

Meanwhile, Figure 6 depicts the error between the observer state and the system state, revealing a similar pattern between them. While exhibiting oscillatory behavior, the rapid reduction of the error, however, is not completely zero due to the presence of perturbations. This suggests that the observer maintains some accuracy despite the attacks.

The timing of event triggers, illustrated in Figure 7, reveals that the maximum trigger interval significantly exceeds the sampling period. This long trigger interval shows the efficacy of the proposed ETM in maintaining system stability and efficiency under DoS conditions. The extended trigger interval also highlights the ETM's capability to manage communication load, triggering transmissions only when essential and thus enhancing safety and stability in the presence of network-based attacks.

5.   Conclusions

In this paper, we have studied the problem of LFC with ACs under DoS attack considering ETM in network environment. When constructing the system model, we have considered the randomness of wind power generation and photovoltaic power generation and have also added the ACs to participate in the frequency regulation. In order to improve the utilization efficiency of network resources while resisting DoS attacks, an ETM has been proposed. During DoS attacks, system stability criteria and control design methods have been derived using $H_{\infty}$ stability theory and LMI techniques for co-designing controller and observer gains. Finally, the proposed control strategy has been applied to the two-area smart grid, and the results have shown that the control method is effective.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

This work was supported by Science and Technology Project of State Grid Anhui Electric Power Co., Ltd. (Grant No. 521205240016).

Conflict of interest

The authors declare there are no conflicts of interest.