Optimal Scheduling Strategy for Urban Distribution Grid Resilience Enhancement Considering Renewable-to-Ammonia Coordination

Abstract

The integration of numerous distributed energy sources into the power system offers exciting opportunities to enhance the resilience of distribution networks. It is worth noting that the renewable-to-ammonia system has the potential to alleviate the multi-temporal and spatial imbalance of the power system. Therefore, this paper proposes a mathematical model for a renewable-to-ammonia system, taking into account the material balance and power balance of each unit. Based on this, this paper further explores the optimization scheduling method for flexible ammonia loads in distribution networks. A relaxation method for branch flow models in distribution networks based on second-order cone programming is proposed. An optimization scheduling model for flexible ammonia loads in distribution networks is constructed to minimize network loss. Moreover, considering the environmental advantages of the renewable-to-ammonia system, this paper compares the changes in hydrogen production technologies under different carbon emission constraints. Finally, a case study of the IEEE 33-node system is adopted to verify the effectiveness of the proposed model and method. It indicates that the renewable-to-ammonia system has environmental benefits and can reduce network loss to a certain extent.

1. Introduction

In recent years, with the accelerated development of renewable energy in China, the inherent spatiotemporal imbalance between energy resources and demands has led to increasingly severe challenges with regard to consumption [1,2,3]. For example, there is a significant mismatch between the distribution of energy resources and energy consumption in China, which is particularly evident in the disparity between wind and solar resources in the Western and Three-North regions and the energy demand in these areas. Specifically, these regions are endowed with abundant renewable energy resources, such as wind and solar power, that are capable of generating substantial amounts of clean electricity. However, due to their relatively low levels of economic development, weak industrial infrastructure, and low population density and energy consumption levels, the demand for electricity in these regions is not high. This lack of demand leads to the underutilization of energy resources, preventing the full potential of renewable energy production from being realized. In addition, its uncertainty also poses threats to the safe and stable operation of the power system [4,5]. Therefore, it is necessary to explore new, adjustable resources on the demand side. Among many adjustment methods, producing hydrogen by water electrolysis has been widely regarded as an efficient means to mitigate renewable energy fluctuation, alleviate the multi-temporal and spatial imbalance of the power system, and achieve more extensive optimization of flexible resources [6,7,8].

Unlike electricity, hydrogen is also a raw material in the transportation, chemical, and construction industries. A series of processes, such as hydrogen storage, transmission, and consumption, still need to be considered after production. Therefore, it is essential to first consider the downstream utilization of hydrogen and then explore the flexibility of water electrolysis in the power system [9,10,11]. It is worth noting that ammonia production, as the largest hydrogen consumption area in China, has a concentrated demand for hydrogen and significant economies of scale, possessing the potential for large-scale consumption of renewable energy. Therefore, these distributed energy resources, such as renewable-to-ammonia systems, provide exciting opportunities to enhance the resilience of urban distribution networks [12,13].

Much work so far has focused on urban distribution network scheduling methods that include flexibility. The flexible resources that can be scheduled are mainly distributed across four aspects: source, network, load, and storage. On the node side, there are various distributed power supplies [14], nodes connected to the upper main network, a variety of controllable load [15,16,17] and energy storage devices [18], and an intelligent soft switch (SOP) on the grid side [19].

From the point of view of the node side, the day-ahead scheduling model of the urban distribution network is constructed to promote the consumption of renewable energy, which effectively reduces the amount of wind and light abandonment and network loss [20]. By integrating the distributed multi-energy storage system into the urban distribution network, the voltage regulation capacity of the urban distribution network is enhanced by the interaction of multi-energy loads. The results show that the proposed strategy can effectively alleviate the expected voltage overlimit problem [21]. A two-stage urban distribution network management method is proposed. The upper layer uses a dynamic optimal power flow model to reduce network loss, and the lower layer optimizes the power injection of the energy storage system to eliminate the imbalance at the common coupling points in the urban distribution network [22]. Considering multiple types of resources on the distribution side, including shared electric vehicle charging stations, distributed photovoltaic power supplies, and energy storage systems, an urban distribution network expansion planning model is proposed [23]. A temperature control characteristic model of air conditioning load is constructed, and a method of reducing pressure and regulating the temperature of the air conditioning load based on a multi-agent system is proposed [24]. In this paper, a load forecasting method for large-scale electric vehicle charging is proposed, and the influence of disordered charging on the network loss and voltage of the urban distribution network is analyzed [25]. In addition, the existing literature has studied the optimal scheduling strategies of various flexible load access urban distribution networks, such as electrolytic hydrogen production [26] and smart buildings [27].

From the point of view of the grid, aiming at the lowest possible transmission loss and reliability cost of the system network, the access scheme of SOP in an AC-DC hybrid distribution system is optimized and determined. The results show that the proposed method can effectively improve the reliability and economy of the system’s operation [28]. Considering SOP, building heat storage, distributed power supply, and other resources, the steady-state SOP model is established, and its access scheme is planned. The example shows that SOP plays a complementary and coordinated role in the thermoelectric system [29]. Aiming at the annual comprehensive cost of the urban distribution network, the planning model of SOP is established to maximize the investment benefit [30]. By combining the energy storage system with SOP, the urban distribution network’s robust optimization models of day-ahead hour level, day–day minute level, and real-time level are constructed to improve the voltage stability of the system [31]. A two-layer cooperative optimization method for the urban distribution network and SOP based on reinforcement learning is proposed, in which topological structure is dynamically selected at the upper layer and topological selection is modified at the lower layer, which effectively improves the flexibility of the urban distribution network operation and enables it to adapt to the uncertain scenario of the source load [32,33,34].

However, the above literature rarely considers the optimal scheduling methods of hydrogen downstream load access in the urban distribution network, ignoring its potential downstream adjustable ability [35,36,37]. Therefore, this paper pays further attention to the flexible operation mode of a renewable-to-ammonia system, which is beneficial for achieving a “grid-friendly” chemical industry cluster. The major contributions of this paper are as follows: First of all, we propose a mathematical model for a renewable-to-ammonia system, taking into account the material balance and power balance of the hydrogen production, storage, and utilization process. Secondly, based on the mathematical model and second-order cone programming, we propose an optimal dispatching model for the urban distribution network, considering the flexibility of the green ammonia industry with the lowest possible network loss as the objective function. Lastly, considering the environmental advantages of the renewable-to-ammonia system, this paper compares the changes in hydrogen production routes under low-hydrocarbon and clean-hydrogen constraints.

2. Process Design and System Modeling of Ammonia Production from Renewable Energy

This section mainly introduces the process flow and system modeling of the renewable-to-ammonia (RE2A) process. The RE2A system is shown in Figure 1.

2.1. Renewable-to-Ammonia Technological Process

The RE2A process mainly consists of an air separation unit, a renewable energy electrolytic hydrogen production system, a hydrogen storage tank, and a synthetic ammonia unit. By harnessing wind and solar energy, water undergoes electrolysis to produce H₂ and O₂, as shown in Equation (1). As a stable gas flow rate is crucial for the synthesis of ammonia, the generated H₂ is directed into the hydrogen storage tank to adjust the gas flow rate within an appropriate range. Subsequently, the H₂ from the storage tank is proportionally mixed with N₂ obtained from the air separation unit, and then the hydrogen-to-nitrogen ratio of the synthesis gas is adjusted to 3:1 through nitrogen addition units. The compressed synthesis gas is introduced into the ammonia synthesis tower, where NH₃ is produced via the H-B process at temperatures of 350–550 °C and pressures of 10–46 MPa, as depicted in Equation (2). Finally, liquid ammonia is separated from the gas mixture during the cooling stage.

$${\mathrm{H}}_2\mathrm{O}\to {\mathrm{H}}_2+{\displaystyle \frac{1}{2}}{\mathrm{O}}_2,\Delta {\mathrm{H}}_{298\mathrm{K}}^{\mathrm{o}}=285.8\mathrm{kJ}/\mathrm{mol}$$

(1)

$${\mathrm{N}}_2+{\mathrm{H}}_2\rightleftharpoons 2{\mathrm{NH}}_3,\Delta {\mathrm{H}}_{298\mathrm{K}}^{\mathrm{o}}=-92.4\mathrm{kJ}/\mathrm{mol}$$

(2)

Compared with the traditional coal ammonia production process, the renewable energy green ammonia production process has differences in raw gas production, power sources, and other aspects. For example, the air separation and compression cooling in the renewable energy green ammonia production process uses electric energy instead of coal to provide power. Apart from the above differences, the process flow of the two is basically the same, and the green ammonia production process from renewable energy is also known as the “green H-B process”.

2.2. Renewable-to-Ammonia System Model

2.2.1. Constraint on Material Flow Balance

According to Equation (2) mentioned earlier, the ratio of hydrogen, nitrogen, and ammonia in the ammonia synthesis reactor is 3:1:2.

$${\mathrm{m}}_t^{{\mathrm{H}}_2,\mathrm{out}}:{\mathrm{m}}_t^{{\mathrm{N}}_2}:{\mathrm{m}}_t^{{\mathrm{NH}}_3}=3:1:2,\forall t\in \mathrm{T},$$

(3)

where $m_t^{{\mathrm{H}}_2,\mathrm{out}}$ is the rate of H₂ consumption in the RE2A process; $m_t^{{\mathrm{N}}_2}$ is the rate of N₂ consumption in the RE2A process; and $m_t^{\mathrm{N}{\mathrm{H}}_3}$ is the rate of NH₃ production in the RE2A process.

2.2.2. Power Balance Constraint

Based on the renewable energy-based ammonia production process, the power generation side primarily includes solar power, wind power, and power systems. On the consumption side, the main components are hydrogen electrolysis, air separation for nitrogen, user loads, and ammonia synthesis. Among these, ammonia synthesis, hydrogen electrolysis, and air separation for nitrogen are flexible loads that can participate in the power balance adjustment on the distribution network side. The formula for this process is as follows:

$$P_t^{\mathrm{WP}}+P_t^{\mathrm{PV}}+P_t^{\mathrm{RE}2\mathrm{A},\mathrm{grid}}=P_t^{\mathrm{EW}}+P_t^{\mathrm{NS}}+P_t^{\mathrm{AS}}, \forall t\in \mathrm{T},$$

(4)

$$P_t^{\mathrm{NS}}=\varphi m_t^{{\mathrm{N}}_2}={\displaystyle \frac{1}{3}}\varphi m_t^{{\mathrm{H}}_2,\mathrm{out}},\forall t\in \mathrm{T},$$

(5)

$$P_t^{\mathrm{AS}}=\psi m_t^{{\mathrm{NH}}_3}={\displaystyle \frac{2}{3}}\psi m_t^{{\mathrm{H}}_2,\mathrm{out}},\forall t\in \mathrm{T},$$

(6)

where $P_t^{\mathrm{WP}}$ is the generation of wind power; $P_t^{\mathrm{PV}}$ is the generation of photovoltaic power; $P_t^{\mathrm{EW}}$ is the power consumption of hydrogen production in the RE2A process; $P_t^{\mathrm{NS}}$ is the power consumption of air separation in RE2A process; $P_t^{\mathrm{AS}}$ is the power consumption of ammonia synthesis in RE2A process; $P_t^{\mathrm{RE}2\mathrm{A},\mathrm{grid}}$ is the purchased electricity quantity for the RE2A process; $\varphi$ is the energy consumption coefficient of air separation in the RE2A process; and $\psi$ is the energy consumption coefficient of ammonia synthesis in the C2A process.

Constraint (4) represents the limit for the power balance at each moment. Constraints (5) and (6) represent the power consumption of nitrogen separation and ammonia synthesis, respectively.

2.2.3. Constraints of the Electrolytic Cells

$$C_{t+1}^{\mathrm{HS}}=m_t^{{\mathrm{H}}_2,\mathrm{in}}={\displaystyle \frac{{\eta }^{\mathrm{EW}}{P}_t^{\mathrm{EW}}·3600}{\mathrm{HHV}}},\forall t\in \mathrm{T},$$

(7)

$${\mu }_t^{\mathrm{EW}}{k}^{\mathrm{EW},\min }{C}^{\mathrm{EW}}\le P_t^{\mathrm{EW}}\le {\mu }_t^{\mathrm{EW}}{k}^{\mathrm{EW},\max }{C}^{\mathrm{EW}},\forall t\in \mathrm{T},$$

(8)

where $m_t^{{\mathrm{H}}_2,\mathrm{in}}$ is the rate of electrolytic hydrogen production or the rate of hydrogen flowing into the hydrogen storage tank; ${\eta }^{\mathrm{EW}}$ is the efficiency of electrolytic hydrogen production; $\mathrm{HHV}$ is the high calorific value of hydrogen; ${\mu }_t^{\mathrm{E}\mathrm{W}}$ is the start–stop variable of the electrolytic cell; $k^{\mathrm{E}\mathrm{W},\mathrm{m}\mathrm{a}\mathrm{x}}$ is the upper limit of the electrolytic cell load; and $k^{\mathrm{E}\mathrm{W},\mathrm{m}\mathrm{i}\mathrm{n}}$ is the lower limit of the electrolytic cell load.

2.2.4. Constraints of the Energy Storage Unit

$$C_{t+1}^{\mathrm{HS}}=C_t^{\mathrm{HS}}+\left({\eta }^{\mathrm{HS}}{m}_t^{{\mathrm{H}}_2,\mathrm{in}}-m_t^{{\mathrm{H}}_2,\mathrm{out}}/{\eta }^{\mathrm{HS}}\right),\forall t\in \mathrm{T},$$

(9)

$$C_0^{\mathrm{HS}}=C_T^{\mathrm{HS}},$$

(10)

$$0\le C_t^{\mathrm{HS}}\le C^{\mathrm{HS}},\forall t\in \mathrm{T},$$

(11)

where $C_0^{\mathrm{HS}}$ is the hydrogen storage capacity of the hydrogen storage tank at the initial time; $C_t^{\mathrm{HS}}$ is the hydrogen storage capacity of the hydrogen storage tank at time t; $C_T^{\mathrm{HS}}$ is the hydrogen storage capacity of the hydrogen storage tank at the final moment; and ${\eta }^{\mathrm{HS}}$ is the inflow and outflow efficiency of hydrogen in the hydrogen storage tank.

2.2.5. Constraints of Ammonia Synthesis Unit

$${\mu }_t^{\mathrm{AS}}{k}^{\mathrm{AS},\min }{m}_t^{{\mathrm{NH}}_3}\le m_t^{{\mathrm{H}}_2,\mathrm{out}}\le {\mu }_t^{\mathrm{AS}}{k}^{\mathrm{AS},\max }{m}_t^{{\mathrm{NH}}_3},\forall t\in \mathrm{T},$$

(12)

$$-r^{\mathrm{down}}{m}_t^{{\mathrm{NH}}_3}\le m_{t+1}^{{\mathrm{H}}_2,\mathrm{out}}-m_t^{{\mathrm{H}}_2,\mathrm{out}}\le r^{\mathrm{up}}{m}_t^{{\mathrm{NH}}_3},\forall t\in \mathrm{T},$$

(13)

where ${\mu }_t^{\mathrm{AS}}$ is the start–stop variable of the synthetic ammonia reactor in the RE2A process; $r^{\mathrm{down}}$ is the landslide rate of hydrogen flowing into the ammonia reactor; and $r^{\mathrm{up}}$ is the climb rate of hydrogen flowing into the ammonia synthesis reactor.

Constraint (12) represents the lower and upper limits of the hydrogen outflow rate. Constraint (13) represents the ramp-up and ramp-down limits for ammonia synthesis reactors.

2.2.6. Renewable Energy Output Constraints

$$0\le P_t^{\mathrm{WP}}\le \sum P_{t,\max }^{\mathrm{WP}},\forall t\in \mathrm{T},$$

(14)

$$0\le P_t^{\mathrm{PV}}\le \sum P_{t,\max }^{\mathrm{PV}},\forall t\in \mathrm{T},$$

(15)

where $P_{s,t,max}^{\mathrm{WP}}$ is the maximum amount of electricity generated by wind power; $P_{s,t,max}^{\mathrm{PV}}$ is the maximum amount of electricity generated by photovoltaics.

3. Distribution Network Scheduling Model

3.1. Objective Function

Electrolytic hydrogen production and synthetic ammonia can be used as flexible resources that can be dispatched in the distribution network to promote the consumption of renewable energy and adjust the power balance. This section considers the minimum network transmission loss as the objective function to optimize the operation mode of the distribution network. The objective function can be expressed as follows:

$$\min {\displaystyle \sum _{t\in \mathrm{T}}{P}_t^{\mathrm{Loss}}}$$

(16)

$$P_t^{\mathrm{Loss}}={\displaystyle \sum _{ij\in E}{I}_{ij,t}^2}{r}_{ij}, \forall t\in \mathrm{T}$$

(17)

where ${ P}_t^{\mathrm{Loss}}$ is the sum of the active losses; $I_{ij}$ is the tributary current; $r_{ij}$ is the branch resistance.

3.2. Operational Constraints

The operational constraints include renewable-to-ammonia system constraints, branch power flow constraints of distribution networks, power balance constraints, and carbon emission constraints, which are introduced in sequence below.

3.2.1. Renewable-to-Ammonia System Constraints

Constraints related to ammonia production from renewable energy sources have been elaborated on in the previous chapters. They mainly include the material flow balance constraint Formula (3), the electrolytic cell constraint Formulas (7) and (8), the power consumption constraint Formulas (5) and (6), the hydrogen storage tank constraint Formulas (9)–(11), the renewable energy output constraint Formulas (14) and (15), the ammonia production upper and lower limit constraint Formula (12), and the synthetic ammonia climbing constraint Formula (13).

3.2.2. Distribution Network Branch Power Flow Constraints

The power flow model for distribution network branches can be described as follows:

$$P_{j,t}^{\mathrm{in}}={\displaystyle \sum _{l\in \mathsf{\delta }(j)}{P}_{jl,t}}-{\displaystyle \sum _{i\in \mathsf{\pi }(j)}\left(P_{ij,t}-I_{ij,t}^2{r}_{ij}\right)}, \forall j\in \mathrm{B}$$

(18)

$$Q_{j,t}^{\mathrm{in}}={\displaystyle \sum _{l\in \mathsf{\delta }(j)}{Q}_{jl,t}}-{\displaystyle \sum _{i\in \mathsf{\pi }(j)}\left(Q_{ij,t}-I_{ij,t}^2{x}_{ij}\right)}, \forall j\in \mathrm{B},$$

(19)

$$V_{j,t}^2=V_{i,t}^2-2\left(P_{ij,t}{r}_{ij}+Q_{ij,t}{x}_{ij}\right)+I_{ij,t}^2\left(r_{ij}^2+x_{ij}^2\right), \forall ij\in \mathrm{E}$$

(20)

$$I_{ij,t}^2={\displaystyle \frac{P_{ij,t}^2+Q_{ij,t}^2}{V_{i,t}^2}}, \forall ij\in \mathrm{E},$$

(21)

$${\underset{¯}{I}}_{ij}\le I_{ij,t}\le {\overline{I}}_{ij}, \forall ij\in \mathrm{E}$$

(22)

$${\underset{¯}{V}}_j\le V_{j,t}\le {\overline{V}}_j, \forall j\in \mathrm{B}$$

(23)

where $P_{j,t}^{\mathrm{in}}$ is node j, which injects active power; $P_{ij,t}$, $P_{jl,t}$ is the active power at the head of the branch ij, jl; $Q_{j,t}^{\mathrm{in}}$ is node j, which injects active power; $Q_{ij,t}$, $Q_{jl,t}$ is the reactive power at the head end of branch ij, jl; $x_{ij}$ is branch reactance; I_ij is the lower branch current limit; ${\overline{I}}_{ij}$ is the upper limit of the branch current; $V_{i,t}$, $V_{j,t}$ is the node voltage; V_j is the node voltage lower limit; ${\overline{V}}_j$ is the node voltage cap; $B$ is a collection of all nodes; $E$ is a collection of all the slip roads; $\delta \left(j\right)$ is a set of branch end nodes with j as the first node; and $\pi \left(j\right)$ is a collection of nodes at the head of a branch with j as the end node.

The above constraints can be divided into equality constraints and inequality constraints, wherein the equality constraints are the active power balance constraint (18), the reactive power balance constraint (19), the voltage and current relationship constraint (20), and the power and voltage and current relationship constraint (21). The inequality constraints are the branch current safety constraint Formula (22) and the node voltage safety constraint Formula (23).

It should be noted that BFM type (21) is a quadratic constraint, resulting in a nonlinear programming model. Here, let $\tilde{I}$_ij,t = $I_{ij,t}^2$ and $\tilde{V}$_ij,t = $V_{ij,t}^2$ and relax the equality constraint to an inequality constraint, which can be written as follows:

$$P_{ij,t}^2+Q_{ij,t}^2\le {\tilde{I}}_{ij,t}{\tilde{V}}_{j,t}, \forall ij\in \mathrm{E}$$

(24)

Second-order conic relaxation (SOCR) is performed on Equation (24) by introducing the second-order conic constraint.

$${‖\begin{array}{c}2P_{ij,t}\\ 2Q_{ij,t}\\ {\tilde{I}}_{ij,t}-{\tilde{V}}_{j,t}\end{array}‖}_2\le {\tilde{I}}_{ij,t}+{\tilde{V}}_{j,t}, \forall ij\in \mathrm{E}$$

(25)

The optimal power flow problem has completed its transformation from a nonlinear programming problem to a convex optimization problem, and the difficulty of solving it has significantly decreased.

3.2.3. Power Balance Constraints

The power balance constraints can be described as follows:

$$P_{i,t}^{\mathrm{in}}=P_{i,t}^{\mathrm{WP}}+P_{i,t}^{\mathrm{PV}}-P_{i,t}^{\mathrm{EW}}-P_{i,t}^{\mathrm{NS}}-P_{i,t}^{\mathrm{AS}}-P_{i,t}^{\mathrm{L}}, \forall i\in \mathrm{B},t\in \mathrm{T}$$

(26)

$$Q_{i,t}^{\mathrm{in}}=-Q_{i,t}^{\mathrm{L}}, \forall i\in \mathrm{B},t\in \mathrm{T}$$

(27)

where $P_{i,t}^L$ is the active power required for the load at node i; $Q_{i,t}^L$ is the reactive power required for the load at node i.

3.2.4. Carbon Emission Constraints

Compared with the traditional coal ammonia production process, the advantage of renewable energy ammonia production is that it is clean with no carbon emission, and that it has good environmental benefits. Therefore, the carbon emission constraint of hydrogen production is set here, and the change in hydrogen production methods from ammonia from coal to ammonia from renewable energy sources is compared under the constraints of the no-carbon-emission constraint, the low-hydrocarbon constraint, and the clean hydrogen and renewable hydrogen constraint.

$$f^{\mathrm{coal}}{\displaystyle \sum _{i\in \mathrm{B}}{\displaystyle \sum _{t\in \mathrm{T}}{\dot{n}}_{i,t}^{{\mathrm{H}}_2}}}\le \mathrm{M}{\displaystyle \sum _{i\in \mathrm{B}}{\displaystyle \sum _{t\in \mathrm{T}}({\dot{n}}_{i,t}^{{\mathrm{H}}_2}+{\dot{n}}_{i,t}^{{\mathrm{H}}_2,\mathrm{out}})}},$$

(28)

where $f^{\mathrm{coal}}$ is the carbon emission coefficient of coal-to-hydrogen production, 25 kgCO₂/kgH₂; M is the carbon emission threshold for hydrogen production.

According to the green hydrogen certification standard issued by the China Hydrogen Energy Alliance, when considering the low-carbon hydrogen constraint, the value of M is 14.51 kgCO₂/kgH₂. When considering the constraints of clean hydrogen and renewable hydrogen, the value of M is 4.9 kgCO₂/kgH₂.

4. Case Study

4.1. Data Description

In this example, the IEEE33-node system is used to further verify the optimal scheduling method of the distribution network with a flexible load of ammonia production from renewable energy sources. In the IEEE33-node system, two renewable energy ammonia plants, RE2A1 and RE2A2, are connected at nodes 33 and 18, respectively, where renewable power is provided by wind power plants and photovoltaic power plants. The thermal power plant is connected to node 19, and the photovoltaic power plant is connected to node 25. In addition, a coal-to-ammonia plant is connected at node 33 to compare the change in the hydrogen production technology path under the consideration of carbon emission constraints, as shown in Figure 2. In this example, the optimal scheduling period is 24 h, and the step is set to 1 h. It can be seen from the above that this problem is a mixed-integer second-order cone programming problem, which can be solved by CPLEX on the MATLAB-YALMIP R2020b platform.

In this example, three scenarios are set up to compare the impact of the renewable energy ammonia production system on the distribution network from different angles, as shown below:

• Scenario 1: RE2A1 and RE2A2 provide renewable energy generation but are not used for ammonia production; C2A does not operate.
• Scenario 2: RE2A1 and RE2A2 provide renewable energy generation and are used for ammonia production; C2A does not operate.
• Scenario 3: RE2A1 provides renewable energy generation and is used for ammonia production; RE2A2 provides renewable energy generation but is not used for ammonia production; C2A operates ammonia production.

In the above scenario, scenario 1 and scenario 2 compare the influence of the synthetic ammonia flexible load on the safety state of the distribution network operation and show the optimal scheduling results for hydrogen ammonia. Scenario 3 compares the scheduling of RE2A and C2A in the distribution network under different carbon emission constraints for hydrogen production.

4.2. Optimizing Scheduling Results

This section mainly analyzes the results of the example based on three aspects: distribution network security, hydrogen and ammonia dispatching results, and environmental benefits.

4.2.1. Distribution Network Safety Analysis

This paper mainly reflects the safety state of lines and nodes from the network loss, power flow, and voltage distribution.

Table 1. Comparison of network losses in different scenarios.

	Scenario 1	Scenario 2
Network loss (kWh)	971.19	956.70

shows the network loss comparison between scenario 1 and scenario 2.

In scenario 1, wind power and photovoltaic power are not used for ammonia production, and the load at each node in the distribution network is constant. Therefore, the network loss of the distribution network in scenario 1 is relatively high. In scenario 2, the flexible load of electrolytic hydrogen production and synthetic ammonia can absorb part of the electricity when the wind or light resources are abundant, which can reduce the network loss to a certain extent.

In addition, this example analyzes the branch power flow distribution and node voltage distribution of the distribution network, considering the industrial flexibility of green ammonia, as shown in Figure 3 and Figure 4.

As can be seen from the figures, there is no blocking of the lines in the distribution network within one day, and the voltage value of each node is 1 p.u. It fluctuates smoothly in the range of 1.05 p.u., and no voltage is too low or over the limit. Therefore, flexible and controllable load and power flow optimization can improve the stability of the distribution network operation.

4.2.2. Analysis of Hydrogen and Ammonia Dispatch Results

The production of hydrogen ammonia is shown in Figure 5. It can be found that from the perspective of total output, there is surplus renewable electricity for hydrogen production in the night gale period or in the noon sufficient light period, while in the evening peak period of electricity consumption, hydrogen production decreases significantly. The production of ammonia in the first few periods is relatively small, mainly due to the high power consumption of hydrogen production and the ability to produce ammonia after the completion of hydrogen production, while the production trends of ammonia and hydrogen in the latter period are similar.

From the perspective of RE2A systems connected at different nodes, the hydrogen production trends of RE2A1 and RE2A2 are similar. In contrast, the ammonia production of RE2A1 is mainly concentrated in the later period, while the ammonia production of RE2A2 is mainly concentrated in the earlier period. This is mainly because hydrogen production is closely related to the output trend of renewable energy, and there are hydrogen storage tanks in the process from the completion of hydrogen production to ammonia synthesis, and the optimal time period for ammonia production at different nodes can be determined according to the network loss of the line.

4.2.3. Environmental Benefit Analysis

Compared with the C2A system, the biggest advantage of the RE2A system is that it is clean and low-carbon. Figure 6 and Figure 7, respectively, show the production of hydrogen under the low-hydrocarbon constraint and under the clean and renewable hydrogen constraint.

When there is no carbon emission constraint, all hydrogen is produced by C2A. The hydrogen production of RE2A and C2A is equivalent after considering the low-hydrocarbon constraint. After considering the constraints of clean hydrogen and renewable hydrogen, the hydrogen production of RE2A increased significantly. This indicates that in the context of the low-carbon energy transition, RE2A will have more significant advantages.

In summary, considering the optimal scheduling of the green ammonia industry flexibility in the distribution network can reduce network loss to a certain extent and ensures that the power flow and voltage in the distribution network are in a safe state. In addition, the technical route of ammonia production from renewable energy can promote the transformation of the chemical industry to a low-carbon energy economy.

The choice of ammonia production pathways primarily aims to highlight the significant environmental advantages of using renewable energy for ammonia synthesis compared to traditional coal-based methods. Utilizing renewable energy for ammonia synthesis can effectively reduce greenhouse gas emissions and lower environmental pollution, aligning with global sustainability goals and contributing to climate change mitigation. In contrast, traditional coal-based ammonia production consumes substantial amounts of fossil fuels and emits significant quantities of carbon dioxide and other harmful substances, placing a considerable burden on the environment.

In this context, the focus of this discussion is to showcase the environmental benefits of renewable energy-based ammonia production, emphasizing its positive role in reducing carbon footprints and promoting the development of green industry. Consequently, although economic factors such as production costs, investment requirements, and operational expenses are important considerations for practical implementation, this discussion does not delve into a detailed comparison of these costs. The primary objective is to underline the environmental advantages of renewable energy technologies, rather than to compare the economic costs of various production pathways.

This paper uses the IEEE33-node system for case study analysis primarily to demonstrate the feasibility of the renewable energy-based ammonia production system in reducing network losses and enhancing system security within the distribution network. The analysis indicates that the renewable energy ammonia production system has significant regulatory potential and can also be applied to other large-scale distribution network systems, such as IEEE141, due to their similar topological structures. The focus of this paper is to showcase the effectiveness of the proposed method in distribution networks, and therefore, validation of the case study was not performed across multiple systems.

5. Conclusions

In summary, the conclusions drawn from the paper are as follows: The renewable energy-based ammonia production system, as a new type of flexible resource on the demand side, can reduce network losses to some extent and can enhance the safety and resilience of urban power grids through the optimal scheduling of hydrogen and ammonia production. Additionally, compared to traditional coal-based ammonia production, the renewable energy ammonia production system offers significant environmental benefits.

This paper addresses the scheduling issues of green hydrogen chemical processes integrated into the power system, with ammonia as the final product, and provides optimization strategies for scheduling. Given the limited research on green hydrogen chemical systems in the power sector, future research could focus on the following aspects:

(1)

This paper focuses on ammonia as a downstream product of hydrogen, and the findings are specific to this product. However, other chemical products, such as methanol, which have similar production processes to ammonia, could also be studied as downstream chemical products in further research.

(2)

This paper only considers optimization scheduling strategies for green hydrogen chemical processes aimed at enhancing the safety of active distribution networks. As a flexible resource on the demand side, green hydrogen chemicals can provide ancillary services to the power system, such as peak shaving, frequency regulation, and reserves, due to their flexible control technologies. Exploring these market opportunities could enhance the market competitiveness of green hydrogen chemicals. Therefore, further research on the role of green hydrogen chemical systems in power balance regulation is needed.

(3)

This paper considers hydrogen energy solely as an electrification carrier for the chemical industry. However, green hydrogen chemical systems involve multiple sectors, including power generation, hydrogen production, and chemical processes, with transactions across various markets such as electricity, hydrogen, carbon, and ammonia. Different sectors may be managed by different market investors. Therefore, a key focus for future research will be on exploring how to reasonably coordinate the distribution of benefits while ensuring overall efficiency.