Optimization Algorithms for Urban Energy Balance: The Application of Thermodynamic Methods in Low-Carbon Ecological City Design

Optimization Algorithms for Urban Energy Balance: The Application of Thermodynamic Methods in Low-Carbon Ecological City Design

Xuefeng Wang Yugui Jia* Mingyue Shao Zongyue Chang Pan Zhang

Hebei Key Laboratory of Water Quality Engineering and Comprehensive Utilization of Water Resources, Hebei University of Architecture, Zhangjiakou 075000, China

Hebei Renewable Energy Heating Engineering Research Center, Hebei University of Architecture, Zhangjiakou 075000, China

Corresponding Author Email: 
wxf1213@hebiace.edu.cn
Page: 
949-959
|
DOI: 
https://doi.org/10.18280/ijht.420323
Received: 
8 January 2024
|
Revised: 
26 April 2024
|
Accepted: 
2 June 2024
|
Available online: 
27 June 2024
| Citation

© 2024 The authors. This article is published by IIETA and is licensed under the CC BY 4.0 license (http://creativecommons.org/licenses/by/4.0/).

OPEN ACCESS

Abstract: 

As the severity of global climate change and the energy crisis intensifies, constructing sustainable urban energy systems has become a focal point of global attention. The optimization of urban energy balance is a crucial step in achieving low-carbon ecological city design. Its objectives are to reduce energy consumption and carbon emissions while ensuring the stability and economy of urban energy supply. This paper focuses on the application of thermodynamic methods in the design of low-carbon ecological cities, aiming to enhance the efficiency and sustainability of urban energy systems through optimization algorithms. Current research is concentrated on exploring efficient management of urban energy flows, yet faces challenges in computational efficiency and model adaptability concerning thermodynamic system modeling and energy allocation strategy optimization. This study first establishes an urban energy system model based on comprehensive thermodynamic principles, considering the diversity and complexity of urban energy consumption. Subsequently, it proposes an optimization algorithm based on multi-energy flow balance to improve urban energy allocation strategies. This algorithm optimizes for system dynamics and non-linearity issues, and enhances computational efficiency for processing large-scale data. The findings of this paper not only provide new theoretical perspectives and computational tools for optimizing urban energy balance but also offer scientific decision-making support for practical low-carbon ecological city planning and management, demonstrating significant theoretical innovation and practical application value.

Keywords: 

urban energy systems, thermodynamics, low-carbon ecological city, energy balance optimization, multi-energy flow balance, computational efficiency

1. Introduction

In today's era, faced with global climate change and an energy crisis, establishing sustainable energy systems has become a core issue in urban development planning [1-9]. Cities, as the main venues for energy consumption, have their energy structures and consumption patterns directly related to the level of carbon emissions [10-12]. Therefore, in the field of ecological city design, the research and implementation of optimizing urban energy balance become particularly critical. This involves not only the effective supply and consumption of energy but also the application of thermodynamic principles in urban energy management and the interaction between these principles and urban planning strategies.

The significance of related research lies in its ability to provide urban planners with a powerful tool to simulate, analyze, and optimize energy flows in a scientific and systematic manner [13, 14]. Through a deep understanding of urban energy systems, more efficient and low-carbon energy allocation schemes can be constructed, thus promoting cities to develop in a greener and more sustainable direction [15-17]. Such research not only helps reduce greenhouse gas emissions and improve the urban ecological environment but also enhances the efficiency of energy use, providing urban residents with more reliable and economical energy supply.

However, existing research methods have some flaws and deficiencies. For instance, in thermodynamic system modeling, many models have not fully considered the challenges brought by urban scale and complexity [18-21]; in terms of energy allocation strategy optimization, the dynamic characteristics and non-linear issues of multi-energy flow balance often have not been effectively addressed [22-25]. Moreover, the computational efficiency of existing methods in processing large-scale data is often not high, which limits their application in actual large urban energy systems.

The main content and research value of this paper are as follows: first, we have constructed a thermodynamic system model for low-carbon ecological city design, which comprehensively considers the diversity and complexity of urban energy systems and fully utilizes thermodynamic principles, offering a new perspective for urban energy balance. Second, we propose an optimization algorithm for urban energy allocation strategy based on multi-energy flow balance, which not only considers the dynamism and non-linearity of the energy system but also improves computational efficiency, making the algorithm applicable to large-scale urban energy planning. Through these two parts of the research, this paper not only promotes the development of low-carbon ecological city design theory but also provides feasible solutions for actual urban planning and management, having significant theoretical and practical implications.

2. Thermodynamic System Modeling for Low-Carbon Ecological City Design

Figure 1 presents the structure diagram of the low-carbon urban energy system. This paper makes significant adjustments to traditional thermodynamic system modeling methods to suit the needs of low-carbon ecological city design. Traditional thermodynamic systems mainly include heat sources, primary networks, heat exchange stations, and heat loads, focusing on the production, transmission, exchange, and final use of thermal energy. Figure 2 shows the relationships between the basic components of the thermodynamic system. However, in the design of low-carbon ecological cities, thermodynamic system modeling needs to more comprehensively consider energy efficiency, environmental impact, and the integration of renewable energies. Firstly, this paper focuses on the primary network, i.e., the main transmission system in the thermodynamic system, and analyzes the heat exchange stations as heat loads. Considering that the secondary network can be neglected due to its shorter length and relatively lower energy consumption, this model emphasizes optimizing the efficiency of the primary network and reducing its energy loss. In the low-carbon modeling approach, this paper not only focuses on the recycling of thermal energy but also enhances the system's capacity to absorb renewable energies, such as solar or geothermal energy, thereby reducing dependence on fossil fuels. Moreover, thermodynamic system modeling in low-carbon ecological city design also emphasizes system operation optimization, including the use of advanced thermodynamic principles and algorithms to predict and manage the flow of thermal energy, ensuring maximum energy efficiency and minimal environmental impact of the heat network. This approach, while considering the complexity of urban energy systems, will also introduce more data and real-time monitoring to achieve dynamic optimization and higher operational flexibility of the thermodynamic system.

Figure 1. Structure diagram of low-carbon urban energy system

Figure 2. Basic structure of the thermodynamic system

In the design of low-carbon ecological cities, the heat source node needs to consider the integration and efficiency of renewable energies. When calculating the heat supply in a fixed period, it is necessary not only to evaluate the output of traditional heat sources but also to consider the contribution and reliability of renewable heat sources such as solar and geothermal energy. The algorithm needs to be able to dynamically adjust the heat source output based on real-time data of renewable energies to reduce reliance on non-renewable energies. Load node modeling should take into account the energy efficiency and heat retention characteristics of buildings, such as insulation performance, automated energy management systems, etc. The calculation method should comprehensively consider external climatic conditions, usage patterns of buildings, and residents' behavioral habits. Assuming the specific heat capacity of water is represented by $z$, the mass flow rates passing through the heat source and load are represented by $W_{u, s}^{G T}$ asd $W_{u, s}^M$, respectively, the supply water temperature by $S_{u, s}^T$, and the return water temperature by $S_{u, s}^e$. The following equation gives the formula for calculating the heat supplied by source node $u$ in period $s$ :

$g_{u, s}^{G T}=z W_{u, s}^{G T}\left(S_{u, s}^t-S_{u, s}^e\right)$     (1)

The following formula calculates the heat consumed by load node $u$ in period $s$ :

$g_{u, s}^M=z W_{u, s}^M\left(S_{u, s}^t-S_{u, s}^e\right)$    (2)

To increase system flexibility, this paper introduces a thermal power relaxation term, allowing for some flexibility in the amount of heat supplied within a certain range. This allows the system to preheat or store thermal energy during non-peak periods, balancing the total heat supply within the period. By relaxing, the system can absorb the instability of renewable energies and optimize system operation through energy storage and demand-side management technologies. Assuming the allowable fluctuation limits of thermal load at time $s$ are represented by $g^{M, M A X}{ }_{u, s}$ and $g^{M, M I N}{ }_{u, s}$, with the thermal power relaxation coefficients represented by $j_{U P}$ and $j_{D O}$, respectively:

$g_{u, s}^{M, U V} \leq g_{\mu, s}^M \leq g_{u, s}^{M, M U X}$     (3)

$\left\{\begin{array}{l}g_{u, s}^{M, M N}=\left(1-j_{D O}\right) g_{u, s}^M \\ g_{u, s}^{M, M A X}=\left(1+j_{D O}\right) g_{u, s}^M\end{array}\right.$    (4)

Assuming the indoor temperature limits are represented by $s_{X M A X}$ and $s_{x M I N}$, with the actual indoor temperature in period $s$ represented by $S_{I N, s}$. The calculation formulas for $j_{U P}$ and $j_{D O}$ are:

$\left\{\begin{array}{l}j_{D O}=\frac{\left|S_{x M I N}-S_{I N}^f\right|}{S_{I N}^f} \\ j_{U P}=\frac{\left|S_{x M I N}-S_{I N}^f\right|}{S_{I N}^f}\end{array}\right.$    (5)

To avoid always being at the lower limit of indoor temperature during optimization, it is necessary to constrain the average temperature. This measure ensures that comfort in living and working is not compromised while pursuing energy efficiency. The constraint can be implemented by setting a lower limit of temperature or introducing a resident satisfaction model. Assuming the ideal indoor temperature is represented by $S_{I V}$, then the calculation formula is:

$\sum_{s=1}^S S_{I N, s} / S=S_{I N}^s$     (6)

Overall, thermodynamic system modeling for low-carbon ecological city design emphasizes the system's overall efficiency and environmental impact. Compared to traditional thermodynamic system modeling, it focuses more on the integration of renewable energies, enhancement of thermal energy storage, and demand-side response capabilities, as well as the impact of user behavior on the performance of the energy system. Moreover, this modeling approach should also include the ability to monitor and adjust the system in real time to adapt to dynamically changing environmental conditions and user demands.

3. Optimization of Low-Carbon Urban Energy Allocation Strategy Based on Multi-Energy Flow Balance

In the design and optimization process of low-carbon urban energy systems, economic viability is the core determinant of a project's feasibility. Compared with traditional energy systems, low-carbon urban energy systems place more emphasis on integrating various energy flows to ensure efficient operation and minimize environmental impact. For such systems, the premise of allocation optimization is to meet the known total load demand, with the optimization goal being to minimize the overall annualized costs, including initial investment and operational costs. This encompasses not only the expenses of conventional energy facilities but also the investment in renewable energy technologies and potential environmental taxes. In this way, low-carbon urban energy systems aim to achieve long-term economic and environmental benefits simultaneously. Assuming the total annual cost is represented by $Q$, the annualized investment construction cost of equipment by ${ }_\nu Q_{I N K}$ and operational expenses by $Q_{O P E}$, the expression is:

$M I N(Q)=M I N\left(Q_{I N V}+Q_{O P E}\right)$     (7)

When constructing the energy system of a low-carbon ecological city, maintaining energy flow balance and ensuring that each facility has adequate production capacity are two basic operational constraints of the system. These constraints aim to ensure that the city's heating, cooling, and electricity loads can be reliably met, meaning that at any moment, there is a balance between energy supply and demand, and the design capacity of various energy devices can cope with the maximum demand under extreme conditions. Assuming the upper limit of the number of devices is represented by $V_{M A X}$, and the upper and lower limits of device capacity are represented by $Z^{M I N}{ }_u$ and $Z^{M A X}{ }_u$, respectively, thus we have the inequalities:

$0 \leq V_u \leq V_{M A X}$    (8)

$Z_u^{M I N} \leq Z_u \leq Z_u^{M A X}$    (9)

This paper proposes an optimization method for low-carbon urban energy allocation strategy based on multi-energy flow balance. This method combines the Particle Swarm Optimization (PSO) algorithm with the method of equal incremental. In the first phase of optimization, the PSO algorithm is used to find the optimal energy allocation scheme to minimize construction investment costs while meeting energy supply and demand balance and low-carbon constraints. The second phase employs the method of equal incremental to search for the optimal configuration where marginal costs are equal under the premise of energy flow balance. The method of equal incremental is used here to refine the solution obtained in the first phase, to minimize the system's operational costs, including fuel costs, operation and maintenance costs, carbon emission costs, etc.

In the optimization problem of this paper, each particle represents a potential energy configuration scheme, including but not limited to the size and layout of renewable energy facilities, traditional energy facilities, energy storage systems, and energy conversion units. Assuming there are $V$ individuals in a $v$-dimensional space, the $k$-th individual is represented by $V_k$, with its position attribute represented by $M P_k$, and its velocity attribute by $N R_k$, the formula is:

$V_k=\left[M P_k ; N R_k\right]$    (10)

The position of a particle can be mapped to the parameters of an energy configuration scheme, such as the capacity, number, and location of various energy facilities. The specific values of these parameters constitute the particle's coordinates in the solution space. The velocity of a particle represents the speed and direction of changes in the parameters of the energy configuration scheme. Each component of velocity corresponds to a dimension of the particle's position, i.e., the rate of change of an energy configuration parameter. The following expressions give the formulas for position and velocity attributes:

$M P_k=\left[m p_{k 1}, m p_{k 2}, m p_{k 3} \cdots \cdots m p_{k v}\right]$      (11)

$N R_k=\left[n r_{k 1}, n r_{k 2}, n r_{k 3} \cdots \cdots n r_{k v}\right]$   (12)

The update of particle positions is based on their current velocity and the individual and collective historical best positions. The update of velocity takes into account both the cognitive component of the individual and the social imitation part. The formula for updating the positions and velocities of the particle swarm is as follows:

$d\left(M P_u\right)=d\left(m p_{u 1}, m p_{u 2}, m p_{u 3} \cdots \cdots m p_{u v}\right)$    (13)

The velocity found in the previous search represents the particle's velocity in the previous iteration, affecting the velocity update for the current iteration and providing "inertia" to avoid abrupt changes during the search process. The velocity of the current search is the particle's velocity in the current iteration, updated based on the individual and population's historical best positions as well as the previous velocity. Assuming the duration of motion is π, the formula for the previous velocity is:

$N R_k^{u-1}=\left[n r_{k 1}^{u-1}, n r_{k 2}^{u-1}, n r_{k 3}^{u-1} \cdots \cdots \cdots n r_{k n}^{u-1}\right]$   (14)

Each particle's best solution found during the iteration process, i.e., the energy configuration scheme that minimizes construction investment costs up to that point. The formula for calculating an individual's historical best position is:

$M P_k^{B E}=\left[m p_{k 1}^{B E}, m p_{k 2}^{B E}, m p_{k 3}^{B E} \cdots \cdots m p_{k n}^{B E}\right]$    (15)

The best solution found by all particles in the entire swarm, i.e., the best energy configuration scheme among all schemes that minimizes construction investment costs. The formula for calculating the swarm's historical best position is:

$M P_{A L}^{B E}=\left[m p_{A L 1}^{B E}, m p_{A L 2}^{B E}, m p_{A L 3}^{B E} \cdots \cdots m p_{A L n}^{B E}\right]$     (16)

The velocity component from the current position to an individual's historical best position encourages the particle to move closer to its personal historical best position, reflecting the ability to learn from its own experience. The formula for calculating the velocity from the current position to the individual's historical best position is:

$N R_{A L ; B E}^u=\left(M P_k^{B E}-M P_k^u\right) / \pi$     (17)

The velocity component from the current position to the swarm's historical best position encourages the particle to move closer to the entire group's historical best position, reflecting the ability to learn through social interaction and information sharing. The formula for calculating the velocity from the current position to the swarm's historical best position is:

$N R_{A L ; B E}^u=\left(M P_{A L}^{B E}-M P_k^u\right) / \pi$     (18)

Assuming the weight coefficients are represented by $Y_1, Y_2$, $Y_3$, the current velocity can be calculated using the following formula:

$N R_k^u=Y_1 * N R_k^{u-1}+Y_2 * N R_{k, B E}^u+Y_3 * N R_{A L, B E}^u$     (19)

To prevent particles from skipping the optimal solution due to excessive velocity, the magnitude of velocity is usually limited to a range, which is specifically set based on the physical and economic limitations of the energy configuration parameters. That is, when $N R^u{ }_k>N R_{M A X}$, then $N R^u{ }_k=N R_{M A X}$; when $N R^u{ }_k<N R_{M I N}$, then $N R^u{ }_k=N R_{M I N}$. Assuming the current velocity is represented by $N R^u{ }_k$, the calculation for the next position $M P^{u+1}{ }_k$ can be conducted using the following formula:

$M P_k^{u+1}=M P_k^u+\pi \cdot M P_k^u$      (20)

The PSO algorithm searches for the optimal energy configuration scheme through several cycles, with the final position reached by the swarm representing the optimal solution, thus achieving economic, reliable, and low-carbon energy supply.

The equal incremental method is a traditional optimization technique used for economic dispatch in power systems, aiming to find the generator operation combination with the lowest generation cost. In low-carbon urban energy systems, this method can be used to optimize the operation of various energy devices to minimize operational costs while meeting energy demands. In low-carbon urban energy systems, the multi-energy flow balance configuration involves the integration and optimization of various energy sources (such as wind power, solar energy, biomass energy, and traditional fossil fuel energies). Compared to traditional urban energy configurations, the design of low-carbon ecological cities needs to pay more attention to the proportion of renewable energies, the comprehensive utilization efficiency of energy, control of carbon emissions, and the reliability of energy supply.

In the energy balance optimization algorithm for low-carbon ecological cities, considering the dynamic adjustment and optimization of the heating system is crucial for achieving multi-energy flow balance. Unlike the static design of traditional urban energy configurations, the energy system of low-carbon ecological cities needs to have high flexibility and adaptability to respond to real-time changes in load and the variability of renewable energies. In this context, the load allocation for each period is not an independent event starting from zero but is based on a continuous process based on the load allocation of the previous moment. This means that when adjusting the load, the goal is to balance the load difference between the current moment and the previous one, rather than reallocating the entire system's load. Specifically, assuming the electricity and heat load demands at time u are represented by $W_{r r}^u$ and $W_{g r}^u$, with the initial state $W^0=0$, then:

$\left\{\begin{array}{l}\Delta W_{r r}=W_{r r}^u-W_{r r}^{u-1} \\ \Delta W_{g r^{\prime}}=W_{g r^{\prime}}^u-W_{g r}^{u-1}\end{array}\right.$      (21)

For the management of thermal loads, if there is an increased demand, the first step is to assess whether the heating equipment with the lowest operational cost has the capacity to increase the load. If the load of these devices cannot meet the increased demand, it is necessary to consider the next lowest-cost equipment. This process requires real-time monitoring and evaluation of the operational costs and available capacity of the equipment. Similarly, if the thermal load demand decreases, priority should be given to reducing the load of the highest operational cost equipment to maintain the economic efficiency of the energy system. If the reduction capacity of the highest-cost equipment is insufficient to meet the demand decrease, adjustments need to be made to the next highest-cost equipment. Figure 3 shows the energy balance calculation process when thermal load demand increases.

Figure 3. Energy balance calculation process for increased thermal load demand

Figure 4. Principle of thermal load distribution affecting electric load distribution

Unlike traditional urban energy configurations, in the energy system of low-carbon ecological cities, energy stations are usually located close to users to minimize energy transmission distances and losses, especially for thermal and cooling energy. Electric energy, due to its relatively high transmission efficiency, can be balanced through electricity grid purchases and sales. However, because of the lower transmission efficiency of thermal energy, it usually needs to be distributed and adjusted at nearby energy stations. During the process of balancing electric loads, the heating loads of combined heat and power units and heat pumps have been determined, meaning their power supply capabilities and electricity consumption are also fixed. Therefore, the main variables to be determined are the charging and discharging states of battery storage systems and the electricity purchase and sale conditions of the grid. The core of this process is to optimize the use of electricity, ensuring supply and demand balance while enhancing system energy efficiency and reducing environmental impact. If the current power supply load cannot meet the electricity demand and there is enough stored electricity, priority should be given to using battery discharge to balance the electric load. When the battery level is insufficient, consider purchasing electricity from the grid. Conversely, if the power supply load exceeds the electricity demand, priority should be given to storing the excess electricity in batteries, and if the batteries are fully charged, consider selling the excess electricity back to the grid. Figure 4 shows the principle of how thermal load distribution affects electric load distribution.

In the multi-energy flow balance configuration for lowcarbon ecological city design, the equal incremental method is particularly important, being an effective computational method for dealing with nonlinear problems. In the low-carbon urban energy system, especially when involving combined heat and power systems, the system's efficiency changes are highly nonlinear, meaning that the generation efficiency and waste heat utilization efficiency are not constant but vary with the change in equipment workload. Assuming the efficiency of the $k$-th unit at time $u$ is represented by $\lambda_k^u$, and the output of the $k$-th unit at time $u$ is represented by $W^u{ }_k$, the calculation formula is:

$\lambda_k^u=d\left(W_k^u\right)$     (22)

The equal incremental method allows us to precisely predict the efficiency of equipment at the current moment based on the working state of the equipment in the previous moment. This method considers the efficiency changes of different units under different load segments, making the management of the energy system more refined and the optimization results closer to actual working conditions. To maintain the nonlinearity of the equipment model, assuming the output of the $k$-th unit at time $u$-1 is represented by $W^{u-1}{ }_k$, with the initial state $\lambda^0{ }_k=0$, this paper uses the method of equal division, and there are:

$\left\{\begin{array}{l}\lambda_k^u=d\left(W_k^{u-1}\right) \\ \lambda_k^0=0\end{array}\right.$     (23)

In the application of low-carbon ecological city design, using the equal incremental method to update and adjust equipment efficiency not only allows for a more precise matching of supply and demand but also achieves optimization of energy use, thereby reducing energy waste and lowering carbon emissions. Additionally, this method is specifically designed to handle nonlinear problems, capable of adapting to the characteristics of efficiency changes with load in complex systems like combined heat and power. Therefore, this method helps achieve more efficient and greener energy management strategies in the balance configuration of multi-energy flows, which is the core goal pursued by low-carbon ecological city design.

4. Experimental Results and Analysis

The research area of this study is a newly constructed campus (as shown in Figure 5), covering an area of 960 mu. Within the campus boundary, there is a drainage ditch running from north to south, which undertakes the task of draining and flood discharge for part of the eastern bridge area during the rainy season. Normally, the riverbed is exposed, and its water conservation function needs to be enhanced. The climate is characterized by a temperate continental monsoon climate, with an average annual precipitation of 409mm, indicating a severe lack of water resources in this arid and semi-arid region. The research area is divided into two geomorphological units: the plateau above the dam and the area below the dam. The plateau above the dam has an elevation between 1300 and 1600 meters, boasting rich scenic resources, with the potential for development in the tens of millions of kilowatts.

The project implements a "wind energy - electrical energy - thermal energy" conversion and local consumption engineering, with the system design plan as shown in Figure 6. The main construction content of the project includes: two 20MW and one 24MW electrode boilers, dismantling the original coal-fired boiler facilities, constructing new electric boiler rooms and test rooms, building a supporting 110kv substation, and constructing about 4km of 110kv lines, with a total investment of 0.68 million yuan, all of which are corporate investments. Upon completion of the project, the campus can fully achieve heating through wind power, saving 7,818.59 tons of standard coal annually.

Table 1 presents the calculation results for the nodes in the thermal system network of the research area. Based on the data from Table 1, it can be observed that the relative errors in supply temperature are generally very small, all below 0.5%, demonstrating the model's high precision in temperature control. The supply temperatures from node 1 to node 5 remain within a stable range, consistent with the minor changes expected under ideal conditions, indicating the model's effective simulation of thermodynamic behaviors in the real system. Especially noteworthy is node 6, where the supply temperature is accurate to an integer, and the relative error is 0. This could be due to its special role in the system or an optimization result in the calculation process. Regarding thermal power, nodes 2 to 5 all have a small amount of heat output, with their relative errors still maintained at a low level, showing that the model's prediction for heat output is similarly accurate. The thermal power for node 6 is negative, which in thermal systems usually indicates the recovery of heat energy or its transport to other parts of the system. Its relative error is still kept within a reasonable range, further verifying the model's precise capability in simulating energy flow and system efficiency.

Table 2 provides the calculation results for the thermal network pipelines of the research area, including the pipeline numbers, the head and tail node numbers for the pipelines, the flow rates through the pipelines, and the relative errors in flow rate calculations. The data shows that the relative errors in flow rate calculations for all pipelines are kept below 0.5%, indicating the model's high accuracy in simulating the dynamics of water flow. Pipeline 1 has the largest flow rate at 1.2154 kg/s, with a relative error of 0.48%, while pipeline 5 has the smallest flow rate at 0.2635 kg/s, with the smallest relative error of 0.24% among all the pipelines. This distribution of flow rates reflects the flow driven by temperature differences between pipelines or differences in flow demand designed into the system. The precise calculation of flow rates and low levels of error across the pipelines demonstrate that the model can accurately predict the flow of the heat carrier between different pipelines, which is crucial for ensuring effective heat distribution and overall energy efficiency of the system.

Figure 5. Location of the research area

Figure 6. Local consumption wind power heating system

Table 1. Calculation results for the nodes in the thermal system network

Node Number

Supply Temperature (℃)

Relative Error (%)

Thermal Power (MW )

Relative Error (%)

1

128.2658

0.08

0

/

2

126.3569

0.16

0.114

/

3

127.5245

0.11

0.139

/

4

127.7841

0.22

0.125

/

5

123.1215

0.33

0.108

/

6

126

0

-0.52

0.12

 

Table 2. Calculation results for thermal network pipelines

Pipeline Number

Head Node

Tail Node

Pipeline Flow (m/(kg.s-1)

Relative Error (%)

1

6

1

1.2154

0.48

2

1

2

0.2563

0.47

3

1

3

0.3458

0.4

4

1

4

0.5789

0.52

5

4

5

0.2635

0.24

Figure 7 displays the electricity generation efficiency of cogeneration units in a low-carbon city energy system of the research area over a 24-hour period, with two sets of data representing efficiency predictions from linear and nonlinear models. In the linear model, the generation efficiency remains constant at 0.3, representing the use of traditional methods or the assumption that efficiency does not change throughout the day. However, the nonlinear model shows variations in efficiency over time, particularly during the high-demand period from 8:00 to 20:00, where efficiency significantly drops to a minimum of 0.12. This variation reflects the impact of various factors on the efficiency of cogeneration units in actual operation, such as load fluctuations, fuel supply changes, or equipment performance. The efficiency changes in the nonlinear model more realistically represent the performance of cogeneration systems in actual operation, demonstrating the ability to respond to changes in energy demand at different times.

Figure 7. Electricity generation efficiency of cogeneration units in low-carbon urban energy system

Figure 8 shows two sets of data representing the waste heat utilization efficiency of cogeneration units of the research area at different times of the day, from both linear and nonlinear model predictions. The linear model predicts a constant waste heat utilization efficiency of 0.92 throughout the day, a simplification that ignores the complex variations that occur in actual operation. In contrast, the nonlinear model reveals slight fluctuations in waste heat utilization efficiency, especially during the evening and morning periods (16:00-8:00), with efficiency slightly decreasing to a minimum of 0.87. These fluctuations correspond to daily changes in energy demand, such as the reduction in thermal load during the night leading to a decrease in waste heat utilization efficiency. The dynamic predictions of the nonlinear model offer a more refined perspective that reflects efficiency changes due to load variations in actual operation.

Figure 8. Waste heat utilization efficiency of cogeneration units in low-carbon urban energy system

In summary, the multi-energy flow balance urban energy configuration strategy optimization algorithm proposed in this paper demonstrates the effectiveness of capturing the dynamic nonlinear characteristics of waste heat utilization efficiency in cogeneration systems. The nonlinear model, compared to the linear model, can more accurately simulate minor changes in efficiency throughout the day, which are crucial for optimizing energy use and enhancing system energy efficiency. Particularly during low-load periods at night, the nonlinear model correctly reflects the phenomenon of reduced efficiency, aiding in the precise adjustment of energy configuration strategies to reduce energy wastage.

Table 3. Flexibility indicators of energy configuration in a low-carbon urban energy system under four thermal energy supply methods

Thermal Energy Supply Method

Upward Flexibility

Downward Flexibility

Total Flexibility

Electric Output Units

Thermal Output Units

Electric Output Units

Thermal Output Units

1

73265.4

55689

158795

55689

343438.4

2

74512

62315.2

156398

62548.2

355773.4

3

77895.2

62785

189654

61741

392075.2

4

78954.2

68956

185325

63215.2

396450.4

Table 3 compares the flexibility of energy configuration in a low-carbon urban energy system under four different thermal energy supply methods, through three indices: upward flexibility, downward flexibility, and total flexibility. These indices reflect the system's capability to respond to changes in energy demand. Methods 1 and 2 show lower upward and downward flexibility, especially the flexibility of thermal output units compared to electrical output units. With the introduction of new energy systems in Methods 3 and 4, a significant improvement in both upward and downward flexibility is observed, particularly the increase in upward flexibility for electrical output units, indicating that the integration of new energy systems enhances the system's responsiveness to demand fluctuations. Methods 2 and 4 exhibit higher total flexibility compared to Methods 1 and 3, due to the consideration of thermal dynamic balance and the introduction of thermal power relaxation items, suggesting that the application of thermal dynamic balance and relaxation items can further improve the system's adaptability and resilience. The data analysis concludes that the multi-energy flow balance urban energy configuration strategy optimization algorithm proposed in this paper effectively enhances the flexibility of the low-carbon urban energy system. Especially with the introduction of new energy systems and consideration of thermal dynamic balance and thermal power relaxation items, the system's flexibility is significantly improved, indicating that the algorithm can adapt to diversified energy supply methods and complex changes in energy demand.

Figure 9. Comparison of net heat power supply and thermal load in low-carbon urban energy system

Figure 9 illustrates the comparison between the net heat power supply and the thermal load of the research area using Methods 2 and 4, showing good consistency with the changes in thermal load, indicating the optimization algorithm's effectiveness in tracking and meeting thermal load demands. At various times, the net heat power supply in Method 4 is slightly higher than in Method 2, especially during periods of decreasing demand (e.g., 20:00 to 24:00), where Method 4's net heat power supply is closer to the thermal load. This is attributed to the addition of new energy systems and the use of thermal power relaxation items, allowing the heating system to adjust the heat supply more precisely, reducing over or under heating. Moreover, during peak hours (e.g., 4:00), Methods 2 and 4 can provide net heat power supply close to or slightly above the thermal load, ensuring the reliability of heat supply. This analysis concludes that the multi-energy flow balance urban energy configuration strategy optimization algorithm proposed in this paper is efficient and reliable in practice. The algorithm, especially Method 4, which considers new energy systems and thermal power relaxation items, allows the heat supply system to flexibly respond to changes in thermal load, ensuring efficient matching of energy supply and the economic operation of the system.

The electricity load distribution graph of the research area given in Figure 10 shows that cogeneration contributes no electricity from 0:00 to 7:00 and 25:00 to 24:00 but provides stable electricity from 8:00 to 24:00, especially during peak hours (12:00 to 20:00) to meet increased electricity load demands. Additionally, heat pumps consume electricity during the morning peak period (8:00 to 11:00), corresponding to their heat supply needs at the same time, and are inactive during off-peak periods. The grid adjusts its electricity contribution dynamically throughout the day, especially providing basic electricity needs when cogeneration is not operating (0:00 to 7:00 and 25:00 to 24:00) and rapidly increasing supply during peak periods (12:00 to 20:00), particularly from 12:00 to 16:00. Batteries perform minor electricity adjustments at 8:00 and 9:00, showcasing their role as an energy buffer. Overall, a dynamic balance between electricity load and supply is achieved through the reasonable distribution of different energy sources, especially during demand peaks, ensuring electricity load demands are met.

(a) Electricity load distribution

(b) Heat load distribution

Figure 10. Results of urban energy configuration based on multi-energy flow balance in a low-carbon city

In the heat load distribution graph of the research area shown in Figure 10, cogeneration, as the main source of thermal energy, provides stable output from 8:00 to 24:00, especially during peak hours (8:00 to 16:00) to meet the highest heat load demands. Heat pumps supply heat only from 8:00 to 12:00, corresponding to their electricity consumption times, indicating the use of electricity to drive heat pumps to meet thermal load demands. This configuration allows for the reasonable use of heat pumps when electricity load is low, and heat load is high, optimizing energy use efficiency. Overall, the satisfaction of heat load also shows the flexible adjustment capability of the heating system, with strong heat supply capacity during peak periods and reduced output on-demand during off-peak periods, avoiding unnecessary energy wastage.

The comprehensive data analysis concludes that the multi-energy flow balance urban energy configuration strategy optimization algorithm effectively achieves a dynamic balance between electricity and thermal loads. The algorithm considers the dynamism and non-linearity of the energy system, coordinating the output of various energy forms (e.g., cogeneration, heat pumps, grid, and batteries) not only to meet electricity and thermal load demands but also to enhance energy use efficiency and reduce energy wastage. The collaborative operation of different energy forms during peak periods is particularly evident, ensuring the reliability of energy supply and the economic operation of the system.

5. Conclusion

This paper focuses on designing an efficient thermodynamic system model for low-carbon ecological cities and proposes a new urban energy configuration strategy optimization algorithm. The core objective of the research is to achieve optimal configuration of the energy system to meet the city's electricity and heat energy demands while enhancing the system's low-carbon operational efficiency. An integrated energy system model was developed, based on thermodynamic principles and considering the diversity and complexity of urban energy systems. The model aims to provide a new method for urban energy balance analysis. An optimization algorithm based on multi-energy flow balance was introduced, addressing the system's dynamism and non-linearity while also improving computational efficiency, making it suitable for large-scale urban energy planning.

The model's precision and reliability in simulating urban thermodynamic systems were validated through the analysis of network nodes and thermal network pipelines. The electricity generation and waste heat utilization efficiency of cogeneration units in the low-carbon urban energy system were evaluated to determine their contribution to the system's overall efficiency. The flexibility of the low-carbon urban energy system under four different thermal energy supply methods was assessed, and the match between actual net heat power supply and heat load was evaluated. A comprehensive analysis and optimization of the city's energy system's electricity and heat load configurations were performed based on multi-energy flow balance.

The thermodynamic system model and optimization algorithm proposed in this paper provide an efficient and feasible solution for the design of low-carbon ecological city energy systems. With the development of new technologies, such as artificial intelligence and the Internet of Things, future research could further integrate these technologies to enhance the model's predictive accuracy and the optimization algorithm's computational speed. Considering the importance of renewable energy resources in achieving low-carbon city goals, future research could explore more deeply how to better integrate renewable energy resources into urban energy systems.

  References

[1] Mishra, P., Singh, G. (2023). Energy management systems in sustainable smart cities based on the internet of energy: A technical review. Energies, 16(19): 6903. https://doi.org/10.3390/en16196903

[2] Wyrwicka, M.K., Więcek-Janka, E., Brzeziński, Ł. (2023). Transition to sustainable energy system for smart cities—Literature review. Energies, 16(21): 7224. https://doi.org/10.3390/en16217224

[3] Isanbaev, V., Baños, R., Gil, C., Gil, M.M., Martínez, F., Alcayde, A. (2023). Towards energy efficiency in microgrids for smart sustainable cities. In Sustainable Smart Cities and Territories International Conference, Manizales, Colombia, pp. 254-265. https://doi.org/10.1007/978-3-031-36957-5_22

[4] Wang, X., Tan, Y., Zhou, G., Jing, G., Francis, E.J. (2024). A framework for analyzing energy consumption in urban built-up areas based on single photonic radar and spatial big data. Energy, 290: 130202. https://doi.org/10.1016/j.energy.2023.130202

[5] Hasfita, F., Sembiring, E., Damanhuri, E. (2022). The potential of city park waste for biofuel feedstock production: A case study in Bandung city, Indonesia. International Journal of Sustainable Development and Planning, 17(7): 2255-2262. https://doi.org/10.18280/ijsdp.170726

[6] Prasetyo, S.D., Regannanta, F.J., Mauludin, M.S., Arifin, Z. (2023). Economic feasibility of solar-powered electric vehicle charging stations: A case study in Ngawi, Indonesia. Mechatronics and Intelligent Transportation Systems, 2(4): 201-210. https://doi.org/10.56578/mits020402

[7] Ragheb, R.A. (2022). Towards resilience: Energy efficiency in urban communities - Case study of New Borg El Arab city in Alexandria, Egypt. International Journal of Sustainable Development and Planning, 17(3): 795-811. https://doi.org/10.18280/ijsdp.170310

[8] de Oliveira, J.A.P. (2019). Sustainability challenges in an urban century: Can we change urbanization paths to make cities the solutions for rather than the drivers of global problems? Challenges in Sustainability, 7(1): 1-4. https://doi.org/10.12924/cis2019.07010001

[9] Ragheb, G., El-Wahab, M.A., Ragheb, R.A. (2022). Sustainable indicators framework for strategic urban development: A case study of Abu Teeg city in Assiut, Egypt. International Journal of Sustainable Development and Planning, 17(1): 91-107. https://doi.org/10.18280/ijsdp.170109

[10] Nie, Y., Zhang, G., Zhong, L., Su, B., Xi, X. (2024). Urban‒rural disparities in household energy and electricity consumption under the influence of electricity price reform policies. Energy Policy, 184: 113868. https://doi.org/10.1016/j.enpol.2023.113868

[11] Shakarami, K., Rahnama, M.R. (2023). Spatial analysis of the impacts of the urban form on the energy consumption of Karaj over the Covid-19 era (2019-2022). Energy and Buildings, 298: 113568. https://doi.org/10.1016/j.enbuild.2023.113568

[12] Tsirantonakis, D., Chrysoulakis, N. (2022). Earth observation data exploitation in urban surface modelling: The urban energy balance response to a suburban park development. Remote Sensing, 14(6): 1473. https://doi.org/10.3390/rs14061473

[13] Bernard, É., Munck, C.D., Lemonsu, A. (2022). Detailed mapping and modeling of urban vegetation: what are the benefits for microclimatic simulations with town energy balance (TEB) at neighborhood scale? Journal of Applied Meteorology and Climatology, 61(9): 1159-1178. https://doi.org/10.1175/JAMC-D-21-0134.1

[14] Saxena, S., Yaghoobian, N. (2022). Diurnal surface heating and roof material effects on urban pollution dispersion: a coupled large-eddy simulation and surface energy balance analysis. Boundary-Layer Meteorology, 184(1): 143-171. https://doi.org/10.1007/s10546-022-00699-5

[15] Zeng, X., Zhou, Z., Liu, C. (2022). Chinese urban energy and carbon congestion effects: A data envelopment analysis and materials balance approach. Journal of Cleaner Production, 341: 130817. https://doi.org/10.1016/j.jclepro.2022.130817

[16] Pan, T., He, S., Liu, Z., Jiang, L., Zhao, Q., Hamdi, R. (2023). Analyzing Changes in urban green spaces and their effect on land temperature from the perspective of surface radiation energy balance in Rizhao city, the central coast of China. Remote Sensing, 15(19): 4785. https://doi.org/10.3390/rs15194785

[17] Yawale, S.K., Hanaoka, T., Kapshe, M. (2021). Development of energy balance table for rural and urban households and evaluation of energy consumption in Indian states. Renewable and Sustainable Energy Reviews, 136: 110392. https://doi.org/10.1016/j.rser.2020.110392

[18] Oliveira, A., Lopes, A., Niza, S., Soares, A. (2022). An urban energy balance-guided machine learning approach for synthetic nocturnal surface Urban Heat Island prediction: A heatwave event in Naples. Science of the Total Environment, 805: 150130. https://doi.org/10.1016/j.scitotenv.2021.150130

[19] Ali-Toudert, F. (2021). Exploration of the thermal behaviour and energy balance of urban canyons in relation to their geometrical and constructive properties. Building and Environment, 188: 107466. https://doi.org/10.1016/j.buildenv.2020.107466

[20] Fei, F., Wang, D. (2023). Application of neural network feedforward in fuzzy PI controller for electric vehicle thermal management system: Modeling and simulation studies. Energies, 17(1): 9. https://doi.org/10.3390/en17010009

[21] Bicer, E.A., Schirmer, P.A., Schreivogel, P., Schrag, G. (2023). Electric vehicle thermal management system modeling with informed neural networks. In 2023 25th European Conference on Power Electronics and Applications (EPE'23 ECCE Europe), Aalborg, Denmark, pp. 1-8. https://doi.org/10.23919/EPE23ECCEEurope58414.2023.10264482

[22] Özdemir, T.H., Yaşar, C.F. (2023). Proton exchange membrane fuel cell thermal system modeling and control. In 2023 10th International Conference on Electrical and Electronics Engineering (ICEEE), Istanbul, Turkiye, pp. 491-499. https://doi.org/10.1109/ICEEE59925.2023.00095

[23] Fu, Y., Zhang, S., Zhang, X., Yan, G. (2023). Research on modeling and simulation of environmental control system for multipurpose thermal utilization of fighter aircraft. In China Aeronautical Science and Technology Conference, Wuzhen, China, pp. 476-484. https://doi.org/10.1007/978-981-99-8864-8_45

[24] Deng, J., Su, Y., Peng, C., Qiang, W., Cai, W., Wei, Q., Zhang, H. (2023). How to improve the energy performance of mid-deep geothermal heat pump systems: Optimization of heat pump, system configuration and control strategy. Energy, 285: 129537. https://doi.org/10.1016/j.energy.2023.129537

[25] Zhou, J., Wu, Y., Zhong, Z., Xu, C., Ke, Y., Gao, J. (2021). Modeling and configuration optimization of the natural gas-wind-photovoltaic-hydrogen integrated energy system: A novel deviation satisfaction strategy. Energy Conversion and Management, 243: 114340. https://doi.org/10.1016/j.enconman.2021.114340