## 3 Research Methods
### Coding Method of the Solution  
The solution in the algorithm is the energy supply of the gas-turbine generators, heat-storage equipment, power-storage equipment, absorption refrigeration equipment, compression refrigeration equipment, gas-boiler heating, and exchange power of the main grid in time periods of \(T\), in which the value in each field represents the amount of energy supplied during period \(t\) (\(t\) = 1, 2,..., \(T\)). If the energy supplies of the heat-storage equipment and power-storage equipment are negative, it means that energy storage is used during \(t\); if they are positive, it means the energy is discharged. Additionally, when the power exchanged with the main grid is positive, it means that CCHP purchases electricity from the main grid; if it is negative, the opposite is true. An example of a solution for \(T\) = 2 is shown in Table 2, where its superscript is the \(G\)th generation, the first subscript is the \(i\)th group, the second subscript is the time period \(t\), and the third subscript \(j\) is the variable number of each energy supply unit that ranges from 1 to 7. Thus, \(x_{i,j}^{G}\) is \(P_{j}(t)\) of the \(i\)th group of the \(G\)th generation, in which \(j\) is the absorption refrigeration equipment, compression refrigeration equipment, power-storage equipment, gas-turbine generator, main grid exchange power, heat-storage equipment, and gas boiler in order. The sequenceis the same as the decision variable vector \(X(t)=(C_{AC}(t),\ C_{EC}(t),\ P_{ES}(t),\ P_{GT}(t),\ P_{Grid}(t),\ Q_{HS}(t),\) and \(Q_{GB}(t))\) mentioned in Section 3.1, for example, \(x_{i,1,1}^{G}\) is the \(C_{AC}(1)\) of the \(i\)th group of the \(G\)th generation, \(x_{i,1,2}^{G}\) is the \(C_{EC}(1)\) of the \(i\)th group of the \(G\)th generation, and \(x_{i,2,1}^{G}\) is the \(C_{AC}(2)\) of the \(i\)th group of the \(G\)th generation. When starting the update process, the fields in the solution variables were modified to obtain a better performance of the fitness function. According to the constraint formula of energy balance, i.e., Equations (10)-(12), it can be known that the energy output must be equal to the energy input in each time period. Because the cooling energy loads converts absorption refrigeration \(C_{AC}(t)\) into thermal energy load \(Q_{AC}(t)\) according to Equation (14) and concerts compression refrigeration \(C_{EC}(t)\) into electric energy load \(P_{EC}(t)\) according to Equation (15), the value of absorption refrigeration \(C_{AC}(t)\) and compression refrigeration \(C_{EC}(t)\) must be determined first in the updating process of the algorithm in order to confirm the total demand of the electric energy load and thermal energy load. According to Equation (12), the sum of absorption refrigeration \(C_{AC}(t)\) and compression refrigeration \(C_{EC}(t)\) should be equal to the cooling energy load \(C_{load}(t)\), which determines one of them and the other is also determined \begin{table}
\begin{tabular}{c c} \hline
**Parameter** & **Values** \\ \hline \(K_{OM,GT}\) (NTD/kWh) & 0.1805 \\ \(K_{OM,GB}\) (NTD/kWh) & 0.0459 \\ \(K_{OM,WT}\) (NTD/kWh) & 0.1332 \\ \(K_{OM,PV}\) (NTD/kWh) & 0.0432 \\ \(\mu_{Grid}\) (kg/kWh) & 0.637 \\ \(u_{NC}\) (kg/m\({}^{3}\)) & 2.09 \\ \(P_{NC}\) (NTD/m\({}^{3}\)) & 12.2581 \\ \(\eta_{NC}\) (kWh/m\({}^{3}\)) & 9.87 \\ \(\eta_{GT}\) & 0.35 \\ \(\eta_{recry}\) & 0.5 \\ \(\eta_{GB}\) & 0.89 \\ \(\eta_{ES,c}\) & 0.95 \\ \(\eta_{ES,d}\) & 0.95 \\ \(\eta_{HS,c}\) & 0.98 \\ \(\eta_{HS,d}\) & 0.98 \\ \(\sigma_{ES}\) & 0.04 \\ \(\sigma_{HS}\) & 0.02 \\ \(COP_{AC}\) & 1.38 \\ \(COP_{EC}\) & 3 \\ \(\gamma_{ES,c}\) & 0.2 \\ \(\gamma_{ES,d}\) & 0.4 \\ \(\gamma_{HS,c}\) & 0.2 \\ \(\gamma_{HS,d}\) & 0.4 \\ \hline Time period & Electricity price (NTD/kWh) \\ \hline \multicolumn{2}{c}{7:00–22:00} & 3.33 \\ Others & 1.39 \\ \hline \end{tabular}
\end{table} 
Table 1: The parameter values of the system and the electricity price of the time period.