An Optimal Operation Strategy of Regenerative Electric Heating Considering the Difference in User Thermal Comfort

: Regenerative electric heating has gradually become one of the main forms of winter heating with the promotion of “coal to electricity” project. By fully exploiting its regulating capacity, it can effectively achieve a win–win situation of “peak shaving and valley ﬁlling” on the grid side and “demand response” on the customer side. In order to meet the different heating demands of users, a regenerative electric heating optimization and control strategy is proposed, taking into account the difference in users’ thermal comfort. Firstly, the reasons for the difference in user thermal comfort are analyzed, and the differentiated preference factors are calculated based on the maximum likelihood estimation method to design differentiated heating schemes. Then, a dynamic optimization and control model for regenerative electric heating with comfort and economic evaluation indicators is established and solved by using quantum genetic algorithm. Finally, a numerical example is used for simulation analysis. The research results show that the strategy proposed in this paper can take into account the comfort of customers and the economy of peaking and low load shifting, so that the operation of regenerative electric heating can respond to the different needs of different customer groups, and realize ﬂexible adjustment at any time of the day.


Introduction
Accelerating the electric energy substitution, promoting the electrification of urban and rural residents' consumption, and completely getting rid of the dependence on fossil energy are the key factors to achieve the goal of "carbon peak and carbon neutral" [1].In order to optimize the energy structure of winter heating, replace the bulk coal used for heating, and reduce the emission of air pollutants in heating areas, the strategy of electrical energy substitution for energy supply was proposed by most countries to promote the wide application of electric heating [2].Currently, the mainstream electric heating systems are mainly regenerative electric heating and air-source heat pumps.Compared to air source heat pumps, regenerative electric heating systems operate in a more flexible manner, facilitate interaction with the power grid, and have broad application prospects [3].
The construction of electric heating facilities with heat storage function is an important part of the change in residents' lifestyle.However, due to the lack of reasonable and effective operational control measures, and the lack of scientific coordination between electricity consumption and heat release, the existing operating methods cannot fully utilize the economic advantages of heat storage during low load periods, and it is also difficult to achieve stable and comfortable indoor heating temperatures for users.How to achieve a balance and coordination between economy and comfort, and obtain the optimal scheduling and operation plan with the best revenue indicators is an urgent problem to be solved after the large-scale integration of regenerative electric heating [4,5].genetic algorithm to solve the reactive power optimization configuration of AC/DC hybrid power grids.The author from [26] utilizes the advantages of quantum genetic algorithm to solve the optimization calculation problem of power grid fault diagnosis.Quantum genetic algorithm combines the advantages of quantum theory and classical genetic algorithm, and has better population diversity, a faster rate of convergence, and better global optimization ability than ordinary genetic algorithms.
To sum up, most of the existing optimal operation strategies for electric heating considering comfort and economy were single evaluative analysis, without classifying and refining heat users and deeply exploring the differentiated demands for heating between different groups.In order to meet the different heating demands of users, a regenerative electric heating optimization and control strategy is proposed taking into account the difference in users' thermal comfort.
The purpose of this paper is to classify and refine the thermal users, and deeply explore the differentiated demand for heating among different groups, so as to meet the heating demand of users with different characteristics.Regenerative heating optimization for different users can achieve personalized thermal comfort requirements, energy saving optimization, and flexible response according to the differences in thermal comfort of users.This can improve user satisfaction, reduce energy consumption and operating costs, and provide a better heating experience for different users.
In order to meet the differentiated heating needs of different users, reduce energy consumption and operating costs, and improve user satisfaction, this paper innovatively proposes a dynamic optimization regulation model of regenerative electric heating with comfort and economy as evaluation indexes.Firstly, the working principle of regenerative electric heating is described, and the mathematical model of the typical working mode of regenerative electric heating system is established.Secondly, the reasons of thermal comfort difference are analyzed, the differential preference factor is calculated based on the maximum likelihood estimation method, and the differentiated heating scheme is designed.Then, a dynamic optimal control model of regenerative electric heating with evaluation indexes of comfort and economy is established and solved by quantum genetic algorithm.The results show that the model can make the operation of regenerative electric heating respond to the different needs of different user groups and realize the flexible adjustment of day and night.

Structure and Working Principle of Regenerative Electric Heating Equipment
A typical regenerative electric heating structure is shown in Figure 1, and its working principle is that the user sets the storage power, peak and valley hours, and heating temperature through the controller, taking into account their own thermal comfort and Timeof-use Electricity Price Mechanism.The controller controls the direct heating equipment of the regenerative electric heating system in the low load periods, and controls the heat release speed by adjusting the heat dissipation port and circulating fan to provide direct heating for the user.At the same time, the heating is stored until the thermal storage maximum.During peak load periods, the heat is released by the thermal storage unit to provide heating for users.If the residual heat of the heat storage body cannot maintain the room temperature demand, the direct heating equipment begins to assist in heating, thereby maintaining the indoor temperature within the required range.

Direct Heating Equipment
Two flow directions of the energized heating capacity are created by the controller's modification of the operating power of the direct heating equipment: one directs heat to the user, and the other directs heat to the regenerator's heat storage to help control operation.The electro-thermal conversion relationship can be expressed by the following formula: where: is the sum of the electric heat of the direct heating equipment in the t period; ( ) D H Q t is the heat produced by the direct heating equipment to the user in the t period; . ( ) t is the heat discharged by the accumulator for the time period; ( ) is the operating power of direct heating equipment in the t period; and DH η is the elec- tro-thermal conversion efficiency of direct heating equipment.

Heat Storage Equipment
The heat storage equipment is stored by the direct heat storage equipment during the low power consumption period and released on demand during the peak power consumption period, which plays a buffering role in the indoor temperature regulation.In this paper, the temperature change in the regenerator in the heat storage equipment is converted into heat change.The operating characteristics of the heat storage equipment can be expressed as the relationship between the total heat storage, heat storage/release, and heat loss.The dynamic mathematical model is constrained as follows: where: is the total heat storage capacity of the heat storage equipment in the t period; . ( ) is the heat storage capacity of the regenerator in the t period; . ( ) t is the heat released by the regenerator in the t period; .

Modeling of Regenerative Electric Heating System 2.2.1. Direct Heating Equipment
Two flow directions of the energized heating capacity are created by the controller's modification of the operating power of the direct heating equipment: one directs heat to the user, and the other directs heat to the regenerator's heat storage to help control operation.The electro-thermal conversion relationship can be expressed by the following formula: where: ∑ Q(t) is the sum of the electric heat of the direct heating equipment in the t period; Q DH (t) is the heat produced by the direct heating equipment to the user in the t period; Q REH.C (t) is the heat discharged by the accumulator for the time period; P DH (t) is the operating power of direct heating equipment in the t period; and η DH is the electro-thermal conversion efficiency of direct heating equipment.

Heat Storage Equipment
The heat storage equipment is stored by the direct heat storage equipment during the low power consumption period and released on demand during the peak power consumption period, which plays a buffering role in the indoor temperature regulation.In this paper, the temperature change in the regenerator in the heat storage equipment is converted into heat change.The operating characteristics of the heat storage equipment can be expressed as the relationship between the total heat storage, heat storage/release, and heat loss.The dynamic mathematical model is constrained as follows: where: S REH (t) is the total heat storage capacity of the heat storage equipment in the t period; Q REH.C (t) is the heat storage capacity of the regenerator in the t period; Q REH.D (t) is the heat released by the regenerator in the t period; η REH.C and η REH.D are the heat storage and release efficiency of the heat storage equipment, respectively; and ω REH is the energy coefficient of energy dissipation or self-loss of the heat storage equipment itself to the environment.

Regenerative Electric Heating System Model
The regenerative electric heating system purchases electricity from the power grid converts the electric energy into heat energy through the direct heat equipment and then provides heat energy to the user under two working conditions through the optimal scheduling of the heat storage equipment.
where: Q all (t) is the total heat load demand of users in the t period.
The above content provides a complete description of the operating mode of the regenerative electric heating system, which can not only meet the heating needs of users, but also provides the possibility for the power grid to optimize and regulate it.

Analysis of Thermal Comfort Differences among Users
The thermal demand of users in winter directly affects the formulation and implementation of heating control strategies.There is a significant difference in thermal comfort between different user groups, and age is the leading factor.Older groups usually have a higher demand for a warmer environment due to their poor body regulation ability, while younger groups prefer a cooler room temperature environment due to their stronger metabolic and thermoregulatory abilities.Health conditions also have an impact on thermal comfort: susceptible groups with diseases need to maintain a warm environment to alleviate symptoms, while people in good health are usually adapted to carry out various activities at lower ambient temperatures.Economic level is an important constraint, and better-off groups may be more inclined to provide higher heating temperatures in pursuit of greater comfort, while the less well-off groups may pay more attention to energy saving and cost control, and thus choose lower room temperatures within the acceptable range.In order to satisfy the requirements for thermal comfort of various user groups, the start-and-stop of electric heating equipment should be flexibly adjusted in accordance with individual variances and preferences of users.
Here, this paper divides users into two groups based on their age and physical health status: the young group and the old group, and then uses thermal comfort evaluation index to calculate the optimal thermal comfort for different groups, which was based on the predicted mean vote (PMV): where: I PMV is the PMV index value; T a is the indoor ambient temperature; P v is the relative air humidity; A, B and C are known parameters, which are related to the individual characteristics of users.Through a thermal comfort questionnaire survey of 235 users of different ages in a rural area of northeast China, the data of clothing thermal resistance, work nature, and thermal comfort temperature demand are integrated and analyzed to fit the thermal comfort model parameters for each type of user, and the results are shown in Table 1.According to the PMV value range of [−0.5, 0.5] recommended by ISO7730 as the human comfort temperature range, under the condition that the relative humidity of indoor air is 60%, the human metabolic rate is 1.2 and the air flow velocity was va < 0.2 m/s, the optimal thermal comfort temperature of users under different age groups and clothing thermal resistance is calculated, as shown in Table 2 for different users' own characteristic parameters and corresponding optimal thermal comfort temperature.
Economic factors are a key consideration in the assessment and analysis of users' heating needs.The division by age and physical health status alone is not sufficient to fully reflect the actual situation.In order to create a regenerative electric heating optimal regulation strategy that is suitable for customers based on economic sensitivity, we must further segment these two categories.This will allow us to create the best electricity consumption plan based on the relative weighting of customers' preferences for comfort and economy.In this section, history data about the heating behavior of different types of users was collected, including daily heating electricity bills and indoor temperature, in order to classify the persons being surveyed into high-income group (>150,000 CNY/year), middleincome group (50,000-150,000 CNY/year) and low-income group (<50,000 CNY/year) according to their annual household income.For the youth group, most of them are working people, so they can be divided into high-income group, middle-income group, and low-income group according to their work nature and labor ability; considering that most of the elderly group are retired at home or have lost their labor ability, they are divided into middle-income group and low-income group.

Solution of Thermal Comfort Demand Preference with Multiple Features
The differentiated behavior of users' participation in demand response is mainly characterized by the weighting factors of economic and comfort preferences when deciding heating regulation schemes, i.e., for economically oriented users, they prefer to sacrifice comfort for the lowest cost during the whole winter heating period, and vice versa.In this section, based on the classification by class in Section 3.1, the historical data of user heating behavior information is classified, and the maximum likelihood estimation method is used to calculate the differentiated preference factors of users in order to realize the refined simulation of users' differentiated response behavior.The specific solution process is as follows: The discrete set of heating behavior historical data, X, satisfies the probability distribution and can be described as The set of samples contains daily heating electricity bills and indoor temperature and can be described as x = x 1 , x 2 , • • • , x n .The set of demand preference coefficients can be described as and the three together describe the emergence probability of x sample point: The probability of occurrence for the event (X The likelihood function obtains its maximum value when the set of demand preference coefficients takes the value of θ = θ1 , θ2 , • • • , θk , indicating that the selected parameters can reasonably respond to the statistical model, and the set of demand preference coefficients θ at this time is called the maximum likelihood estimate of θ. If the likelihood function is differentiable with respect to the set of demand preference coefficients, the maximum point is obtained by taking the partial derivative of the likelihood function according to the differentiation property of the multivariate function, which is the demand preference coefficients for economy and comfort of the user, as shown in Table 3, from which point it is satisfied that: On this basis, the preference factors of users with multiple characteristics are calculated.The weighted sum representation of the correlation coefficients between each feature parameter (household structure and income status) and thermal comfort temperature are calculated based on Spearman coefficient correlation analysis, which is performed according to the preference coefficients of each category feature and its corresponding user group historical data of heating behavior, as shown in mathematical Formula (5): where ω E , ω C are the user's preference factors for economy and comfort requirements, respectively.

User Thermal Comfort Design Scheme
According to the analysis performed, the thermal comfort design for users has a more reasonable design of their heating temperature range, which significantly improves the economy while satisfying the heating comfort.The thermal comfort design scheme for various user groups according to their thermal comfort demand preferences is shown in Table 4.The vertical comparison shows that a lower temperature range is set during the daytime when nobody is at home, the outdoor temperature is high and the electricity price is more expensive during the peak load period, and a higher temperature range is set at night when people live together, the outdoor temperature is lower and the electricity price is lower during the low load period.The horizontal comparison shows that a higher temperature range can be set for the higher economic level or the aging group, and a lower temperature range can be set for the lower economic level or the younger group.The design of thermal comfort is based on the difference in users' own behavior and avoids the unreasonableness brought by setting the same heating temperature for all groups in 24 h throughout the day, thus achieving the balance of the contradictory factors of heating comfort and economy.

Comfort Indicators
Electricity effect functions of living appliances are divided into two categories: temperature utility functions and time utility functions, which are evaluation indicators for analyzing the quantitative relationship between the utility obtained by electricity users in the process of consuming electric energy and the output of electric energy [27].The temperature utility function is oriented to the temperature-controlled electrical equipment of the energy storage and describes the degree of satisfaction of the user.According to the relationship between the expected temperature and the actual temperature of the user, the temperature utility function is shown in a mathematical Formula (6).(10) where, in the formula: U(t) is the temperature utility function of the t period; ε is the parameter representing the rate of decrease in the utility of the reaction temperature, with a value of (0, 1]; T low and T up are the lower and upper limits of the user's differentiated thermal comfort temperature interval, respectively, taken as in Table 5, in which thermal comfort temperature range was fitted by multiple users features.mean value of temperature utility and temperature effect minimum value was adopted as the user's comfort index, as shown in mathematical Formula (11).(11) where: J 1 is the total indicator of user comfort; U mean is one of the sub-indicators reflecting the mean value of the temperature utility function of user comfort within a regulation cycle; U min is the second sub-indicator reflecting the minimum value of the temperature utility function within a regulation cycle, characterizing the fluctuation in indoor temperature and the difference in the mean value; k 1 and k 2 are the weight distribution coefficients.

Economic Indicators
Peak load Time-of-use Electricity Price can stimulate customers' demand response through price signals, and customers can adjust their electricity consumption behavior to shift more load to the low load periods to reduce heating electricity costs and improve grid operation efficiency.In order to evaluate the economic efficiency of the controlled users under the proposed regulation strategy, the load transfer rate as shown in mathematical Formula ( 8) is proposed as an economic indicator: where: t p and t v are the total hours of peak load and low load periods; Q p (t) is the heat load demand in t hours during the peak load periods; Q v (t) is the heat load demand in t hours during the low load periods; P Rated DH is the rated electrical power of the selected regenerative electric heating equipment; and P DH (t) is the operating power of the equipment during the scheduling periods.

Dynamic Optimal Regulation Model a
Objective function Considering that the main contradiction of users in the winter heating season is "body feeling comfort" and "heating cost under different temperature", this section takes the balance decision of heating electricity cost and comfort of users in a scheduling cycle as the objective function in order to build a dynamic optimization regulation model which takes into account the difference in users' thermal comfort and gives consideration to comfort and economic indicators.It is expressed as: where where: α is the differentiated preference factor for different requirements of heating comfort and economy; the larger the value of α is, the more the objective of dynamic optimization tends to be the user's comfort; the smaller the value of α is, the more the objective of dynamic optimization tends to be the user's economy; M is the preference adjustment control coefficient, and its value range is [0, 1] and its coefficients β reflect the preference factor α on the choice of optimization bias.
b Thermal energy supply and demand balance constraint: See mathematical Formula (4) for details where: T set (t) is the thermal comfort temperature set by different user groups according to their differences in each time period; T out (t) is the local outdoor temperature data of users in the t period; R in−out is the equivalent thermal resistance between indoor-outdoor; and ∆t is the simulation step.
c Direct heating equipment operation constraint: where: P Rated DH is the rated electric power of the selected regenerative electric heating equipment.

d
Operating constraints of heat storage equipment Energy constraint for heat storing and discharging of heat equipment is as follows: See mathematical Formula (3) for details 0 ≤ S REH (t + 1) ≤ S Rated REH (20) where: S Rated REH is the rated capacity of the storage equipment.The heat limit of the heat storing and discharging process constraint is as follows: Heat storage equipment constraint in the regulation process is as follows: where: S REH (1) is the amount of heat stored in the heat accumulator during the initial dispatching period; S REH (T) is the amount of heat remaining in the heat accumulator at the end of the dispatching period.

Quantum Genetic Algorithm
The Quantum Genetic Algorithm (QGA) is an example of the successful application of the quantum computing (QC) concept in the field of the genetic algorithm (GA) [25].To achieve population evolution guided by optimal individual information, the optimization solution of the goal was achieved by utilizing the characteristic that the superposition state of each qubit will shrink and collapse to a certain form during the iteration process.Due to its unique encoding and updating methods, the QGA outperforms the traditional genetic algorithm in all aspects [24], and has the characteristics of small population size, fast rate of convergence, and strong global optimization ability.

Dynamic Optimal Regulation Model Solving Process
It is a nonlinear planning problem to develop a regenerative electric heating optimal regulation strategy that takes into account the differences in the thermal comfort of the users.An algorithm is used to solve the dynamic optimal regulation model developed in this section.Combined with the model built in this chapter, the fitness function is a nonlinear function of optimal comfort or optimal economy directed by the differential preference factor under the influence of preference regulation control coefficients.Quantum bits and quantum logic gates are used to update the chromosome to encode and update the operating state of the direct heat equipment.The solution flow of the dynamic optimal regulation strategy based on QGA is shown in Figure 2, and the specific solution method is as follows:  Step1: Input the initial value or default value of building equivalent parameters; input the equipment parameters; input the user set value.
Step2: Use the temperature forecast data of the next 24 h and calculate the heat demand by mathematical formula (8).
Step3: Apply the QGA according to mathematical formulas ( 9) to ( 17) to solve the dynamic optimal values of regenerative electric heating operating power for 96 time periods in the next 24 h and record the room temperature forecast values corresponding to the optimal solutions.
Step4: Compare the predicted room temperature with the user's thermal comfort temperature interval for each time period: if it meets the user's demand, go to Step5, if not, go back to Step3.Correct the thermal comfort temperature interval for that time period and solve it again.
Step5: Input the optimal values obtained from Step3 into the operation log of regenerative electric heating, and supply heat to customers according to the schedule in each time period based on this optimization result.
The specific steps for solving the optimal value of the input power using QGA are: Step 3-1: Input the original parameters.The input includes the current population evolution number 0 t = and the maximum population iteration number T ; the initial population 0 ( ) Q t with individuals M is randomly generated.
Step 3-2: Initial chromosome coding.Genetic algorithms often use binary encoding to take values of chromosomes, the QGA is based on the introduction of quantum programming using quantum bit probability amplitude to represent the encoding of chromosomes and a quantum bit will be in multiple quantum superposition states with amplitude 1 at the same time, increasing the variation of chromosome values, as shown in the fol- Step1: Input the initial value or default value of building equivalent parameters; input the equipment parameters; input the user set value.
Step2: Use the temperature forecast data of the next 24 h and calculate the heat demand by mathematical Formula (8).
Step3: Apply the QGA according to mathematical Formulas ( 9) to ( 17) to solve the dynamic optimal values of regenerative electric heating operating power {P DH (t 1 ), P DH (t 2 ), • • • , P DH (t 96 )} for 96 time periods in the next 24 h and record the room temperature forecast values corresponding to the optimal solutions.
Step4: Compare the predicted room temperature with the user's thermal comfort temperature interval for each time period: if it meets the user's demand, go to Step5, if not, go back to Step3.Correct the thermal comfort temperature interval for that time period and solve it again.
Step5: Input the optimal values {P DH (t 1 ), P DH (t 2 ), • • • , P DH (t 96 )} obtained from Step3 into the operation log of regenerative electric heating, and supply heat to customers according to the schedule in each time period based on this optimization result.
The specific steps for solving the optimal value of the input power using QGA are: Step 3-1: Input the original parameters.The input includes the current population evolution number t = 0 and the maximum population iteration number T; the initial population Q(t 0 ) with individuals M is randomly generated.
Step 3-2: Initial chromosome coding.Genetic algorithms often use binary encoding to take values of chromosomes, the QGA is based on the introduction of quantum programming using quantum bit probability amplitude to represent the encoding of chromosomes and a quantum bit will be in multiple quantum superposition states with amplitude 1 at the same time, increasing the variation of chromosome values, as shown in the following mathematical formula: The decision variable in the dynamic optimal regulation model established in this paper is the operating power of 96 time periods in the next 24 h.The chromosome structure as described in this paper is divided into 96 subparts, and each subpart represents the individual genes constituted by each chromosome, i.e., the operating power of regenerative electric heating operation power in that time period, so that the i-th gene at k iterations after employing quantum bit encoding is denoted as To make the initial state probabilities of spatial variables equal, α = β = 1 √ 2 is further set.
Step 3-3: Individual evaluation.Measure the number of M individuals in the population Q(t 0 ) and the binary code of the deterministic solution P(t 0 ) can be obtained so as to find the individual fitness value of the population to the corresponding decimal number.
Step 3-4: Record the optimal individual in the population and the corresponding fitness value and take the current optimal individual (the current optimal energy equipment capacity allocation result) as the evolutionary target.
Step 3-5: Quantum revolving gate operation updates the chromosome coding.The QGA updates the chromosome encoding by changing the quantum angle of all chromosome quantum bits in the population and interfering with the quantum state by rotating the plural amplitude, so that the chromosome value converges to the better chromosome, and the quantum revolving gate is shown as follows: Step 3-6: Generate the next generation population.The new population Q(t + 1) is further obtained based on the result of updating the chromosome encoding by using the quantum rotation gate.

Example Analysis 6.1. Example Data
In this section, we analyze an example of the "coal to electricity" conversion project being performed in a rural area of a city in the Liaoning Province of China.The average building area of the town is 82 m 2 and the rated power of regenerative electric heating is 20 kW, taking into account the heat load demand of 100 W/m 2 and the heat storage time of 8 ∼ 10 h per day with a certain margin.The typical daily outdoor temperature data used in this example is the result of the fuzzy C-mean algorithm clustering of the data sampled from the weather observatory near the town from 1 November 2022 to 1 April 2023, as shown in Figure 3.The user set temperature is the thermal comfort temperature interval for each type of user in Table 6.ε is taken as 0.25, k 1 and k 2 are taken as 0.5 when calculating the comfort index, and β is taken as 0.4 in the optimization objective expression.The other main equipment parameters are:  Using the optimal control strategy that takes into account the difference in thermal comfort of users proposed in this paper, the operating power and start-stop time of electric heating equipment are simulated and analyzed, and the optimal control scheme and user thermal comfort satisfaction are compared horizontally and vertically to verify the effectiveness of the method.

Calculation Analysis
By analyzing the objective function, it can be seen that by setting the value of α , users can adjust the operation mode of regenerative electric heating according to their own thermal comfort difference, i.e., adjust the preference of operation optimization results between economy and comfort.In this paper, the evolution algebra is set to 100, and the population number is set to 50 to obtain the preference factor α of the thermal com- fort temperature demand of each user group in Table 3 for three typical days.At the same time, the iterative convergence of the QGA in calculating the preference factor α is ana- lyzed.Figure 4 shows the change in fitness value with the number of iterations.It can be seen from the figure that the population changes with the increase in the number of iterations.At the beginning of the iteration, the fitness value changed greatly, and then most of the fitness values tended to be stable and convergent after 37 generations.The maximum number of iterations was 78 generations.
The response of each optimization index solved by the QGA to the preference factor is shown in Figure 5.The graphs show that the results are in accordance with the predefined preference factor α : the larger the α the more users prefer good heating experi- ence, and the smaller the α the more users prefer good economy at the expense of a certain heating experience.The corresponding preference factor α can be obtained in Figure 5 according to the comfort and economy preference coefficients for the five groups classified in Table 3, which provide parameters to support the subsequent calculation.The preference factor α for the high-income youth group, middle-income youth group, low- income youth group, middle-income elderly group, and low-income elderly group are 2.7, 1.3, 0.5, 1.8, and 1.0, respectively.Using the optimal control strategy that takes into account the difference in thermal comfort of users proposed in this paper, the operating power and start-stop time of electric heating equipment are simulated and analyzed, and the optimal control scheme and user thermal comfort satisfaction are compared horizontally and vertically to verify the effectiveness of the method.

Calculation Analysis
By analyzing the objective function, it can be seen that by setting the value of α, users can adjust the operation mode of regenerative electric heating according to their own thermal comfort difference, i.e., adjust the preference of operation optimization results between economy and comfort.In this paper, the evolution algebra is set to 100, and the population number is set to 50 to obtain the preference factor α of the thermal comfort temperature demand of each user group in Table 3 for three typical days.At the same time, the iterative convergence of the QGA in calculating the preference factor α is analyzed.Figure 4 shows the change in fitness value with the number of iterations.It can be seen from the figure that the population changes with the increase in the number of iterations.At the beginning of the iteration, the fitness value changed greatly, and then most of the fitness values tended to be stable and convergent after 37 generations.The maximum number of iterations was 78 generations.Considering the user thermal comfort design scheme, a rolling optimization decision is made for regenerative electric heating regulation strategy formulation within a scheduling cycle, i.e., the user is allowed to adjust the heating temperature and preference factor in the middle of the day, and it is assumed that the typical temperature data used to formulate the optimal scheduling strategy and the actual temperature obey a normal distribution N(0,0.5)error.The operation of regenerative electric heating is corrected on a The response of each optimization index solved by the QGA to the preference factor is shown in Figure 5.The graphs show that the results are in accordance with the predefined preference factor α: the larger the α the more users prefer good heating experience, and the smaller the α the more users prefer good economy at the expense of a certain heating experience.The corresponding preference factor α can be obtained in Figure 5 according to the comfort and economy preference coefficients for the five groups classified in Table 3, which provide parameters to support the subsequent calculation.The preference factor α for the high-income youth group, middle-income youth group, low-income youth group, middle-income elderly group, and low-income elderly group are 2.7, 1.3, 0.5, 1.8, and 1.0, respectively.To analyze the effect of preference factors on heating comfort and economy, assuming that the ambient temperatures set for the youth group are all 24 C °, and with preference factors α of: 2.7, 1.3, and 0.5, respectively, the operation of regenerative electric heating and indoor ambient temperature changes throughout the day are obtained as shown in Figure 6.The MAE (mean absolute error) and RMSE (root mean square error) between its indoor temperature fluctuations and preset values are: 0. ⋅ .The indicators of various operating parameters show that the larger the value of preference factor α is taken to be, the more inclined to comfort with less room temperature fluctuation, but at the same time, it needs to transfer part of the valley section load to the peak section and lose the economy.
Considering the user thermal comfort design scheme, a rolling optimization decision is made for regenerative electric heating regulation strategy formulation within a sched- To analyze the effect of preference factors on heating comfort and economy, assuming that the ambient temperatures set for the youth group are all 24 • C, and with preference factors α of: 2.7, 1.3, and 0.5, respectively, the operation of regenerative electric heating and indoor ambient temperature changes throughout the day are obtained as shown in Figure 6.The MAE (mean absolute error) and RMSE (root mean square error) between its indoor temperature fluctuations and preset values are: 0.7742, 0.6342, and 0.5086; and 0.9248, 0.7498, and 05034, respectively; and the average temperatures are 24.2774,24.0669, 24.0351; and the electricity consumption in peak load and low periods are: 111.1566 kW • h, 94.1093 kW • h, 75.7083 kW • h and 0 kW • h, 16.9910 kW • h, and 35.43957 kW • h.The indicators of various operating parameters show that the larger the value of preference factor α is taken to be, the more inclined to comfort with less room temperature fluctuation, but at the same time, it needs to transfer part of the valley section load to the peak section and lose the economy.
rolling basis according to the real-time thermal load demand, keeping in mind the combined effect of the temperature error and the variability of thermal comfort temperatures at different times throughout the day, in order to improve the transfer rate of peak load to valley load and to improve the economic efficiency of customers while satisfying comfort.Using the middle-class youth group as an example, rolling optimization decisions are made for three typical days in one scheduling cycle, with users having the option to change the heating temperature and preference monitor at any time.The heating temperature is lowered from 08:00 to 19:00.The optimization results and actual operation results for some time periods are shown in Figures 7-9, and the statistical results of relevant operation parameters and indexes are shown in Table 7.   Considering the user thermal comfort design scheme, a rolling optimization decision is made for regenerative electric heating regulation strategy formulation within a scheduling cycle, i.e., the user is allowed to adjust the heating temperature and preference factor in the middle of the day, and it is assumed that the typical temperature data used to formulate the optimal scheduling strategy and the actual temperature obey a normal distribution N(0,0.5)error.The operation of regenerative electric heating is corrected on a rolling basis according to the real-time thermal load demand, keeping in mind the combined effect of the temperature error and the variability of thermal comfort temperatures at different times throughout the day, in order to improve the transfer rate of peak load to valley load and to improve the economic efficiency of customers while satisfying comfort.Using the middle-class youth group as an example, rolling optimization decisions are made for three typical days in one scheduling cycle, with users having the option to change the heating temperature and preference monitor at any time.The heating temperature is lowered from 08:00 to 19:00.The optimization results and actual operation results for some time periods are shown in Figures 7-9, and the statistical results of relevant operation parameters and indexes are shown in Table 7.Without considering the user thermal comfort design scheme, the optimized results of the operation of the middle-income youth group at Day2 with the same thermal comfort temperature set throughout the day are shown in Figure 10.The actual all-day 24 h operation results differ from the optimized results at the 00:00 starting hour in terms of input power and temperature profiles.According to the optimized results of the operation at the starting time of 00:00, the power consumptions in peak load and low load periods are 41.7961 kW • h, 126.0462 kW • h, and the comfort index J 1 is 0.6172; the actual operation results are 41.7961 kW • h, 127.7921 kW • h, and the comfort index J 1 is 0.6106.With the operating indexes of the design scheme considering the thermal comfort of users, it can be seen that the electricity consumption in peak and low load periods is effectively reduced under the condition of ensuring comfort, which makes the economy and comfort basically consistent with the set expectation.At the same time, the comparison between the design scheme without considering the user's thermal comfort (10) and the design scheme considering the user's thermal comfort (Figure 8) shows that the design scheme considering the user's thermal comfort can effectively reduce the electricity consumption in peak and valley segments, improve the user's comfort index, and reduce temperature fluctuations.
basically consistent with the set expectation.At the same time, the comparison between the design scheme without considering the user's thermal comfort (10) and the design scheme considering the user's thermal comfort (Figure 8) shows that the design scheme considering the user's thermal comfort can effectively reduce the electricity consumption in peak and valley segments, improve the user's comfort index, and reduce temperature fluctuations.

Conclusions
The application of regenerative electric heating is a key part of the popularization of the "coal to electricity" project, and the current typical operation methods cannot meet the differentiated heating needs of different user groups.In this paper, an optimal operation strategy of regenerative electric heating considering the difference in user thermal comfort was proposed, and the specific conclusions are as follows: 1.In order to refine the differentiated behavior of user participation in demand response, the differentiated preference factors are calculated based on the maximum likelihood estimation method to design differentiated heating schemes.Additionally, the thermal comfort scheme was designed based on multiple characteristics of users' thermal comfort preferences.2. The regenerative electric heating optimal control index was intended to achieve a balanced decision between user comfort and economy during the heating duration.

Conclusions
The application of regenerative electric heating is a key part of the popularization of the "coal to electricity" project, and the current typical operation methods cannot meet the differentiated heating needs of different user groups.In this paper, an optimal operation strategy of regenerative electric heating considering the difference in user thermal comfort was proposed, and the specific conclusions are as follows: 1.
In order to refine the differentiated behavior of user participation in demand response, the differentiated preference factors are calculated based on the maximum likelihood estimation method to design differentiated heating schemes.Additionally, the thermal comfort scheme was designed based on multiple characteristics of users' thermal comfort preferences.

2.
The regenerative electric heating optimal control index was intended to achieve a balanced decision between user comfort and economy during the heating duration.Preference criteria are utilized to define different user groups.The QGA was used to solve the dynamic optimal regulation model with the different thermal comfort, and the optimal power of regenerative electric heating was solved for each time period in the scheduling cycle.

3.
The research results show that the strategy proposed in this paper can take into account the comfort of customers and the economy of peak and low load shifting compared with the traditional operation method, so that the operation of regenerative electric heating can respond to the different needs of different customer groups and achieve flexible regulation at all times of the day.

Figure 1 .
Figure 1.Structure diagram of regenerative electric heating.

η
are the heat storage and release efficiency of the heat storage equipment, respectively; and R EH ω is the energy coefficient of energy dissipation or self-loss of the heat storage equipment itself to the environment.

Figure 1 .
Figure 1.Structure diagram of regenerative electric heating.

Figure 2 .
Figure 2. Flow of regenerative electric heating optimal regulation strategy solution with user-difference in thermal comfort.

Figure 2 .
Figure 2. Flow of regenerative electric heating optimal regulation strategy solution with userdifference in thermal comfort.

Figure 4 .
Figure 4. Variation in fitness values with the number of iterations.Different colored lines indicate different comfort levels, varying with the number of iterations.

Figure 5 .
Figure 5. Response of the optimization index to the preference factor (results after fitting for 3 typical days).To analyze the effect of preference factors on heating comfort and economy, assuming that the ambient temperatures set for the youth group are all 24 C °, and with preference factors α of: 2.7, 1.3, and 0.5, respectively, the operation of regenerative electric heating and indoor ambient temperature changes throughout the day are obtained as shown in Figure 6.The MAE (mean absolute error) and RMSE (root mean square error) between its indoor temperature fluctuations and preset values are: 0.7742, 0.6342, and 0.5086; and 0.9248, 0.7498, and 05034, respectively; and the average temperatures are 24.2774,24.0669, 24.0351; and the electricity consumption in peak load and low periods are: 111.1566kW h ⋅ , 94.1093kW h ⋅ , 75.7083kW h ⋅

Figure 4 .
Figure 4. Variation in fitness values with the number of iterations.Different colored lines indicate different comfort levels, varying with the number of iterations.

Figure 4 .
Figure 4. Variation in fitness values with the number of iterations.Different colored lines indicate different comfort levels, varying with the number of iterations.

Figure 5 .
Figure 5. Response of the optimization index to the preference factor (results after fitting for 3 typical days).

Figure 5 .
Figure 5. Response of the optimization index to the preference factor (results after fitting for 3 typical days).

Figure 6 .
Figure 6.Optimized regulation results of regenerative electric heating for different user preferences (Day = 1).

Figure 7 .
Figure 7. Optimized operation results of Day1 for the design scheme of thermal comfort for users ( 1 .3α =

Figure 6 .
Figure 6.Optimized regulation results of regenerative electric heating for different user preferences (Day = 1).

Figure 6 .
Optimized regulation results of regenerative electric heating for different user preferences (Day = 1).

Figure 7 .
Figure 7. Optimized operation results of Day1 for the design scheme of thermal comfort for users ( 1 .3α =

Figure 8 .
Figure 8. Optimized operation results of Day2 for the design scheme of thermal comfort for users (α = 1.3).

Figure 8 .
Figure 8. Optimized operation results of Day2 for the design scheme of thermal comfort for users ( 1.3 α = ).

Figure 9 .
Figure 9. Optimized operation results of Day3 for the design scheme of thermal comfort for users ( 1.3 α = ).

Figure 9 .
Figure 9. Optimized operation results of Day3 for the design scheme of thermal comfort for users (α = 1.3).

Figure 10 .
Figure 10.Optimized operation results of Day2 without considering the user thermal comfort design scheme( 1.3 α = ).

Figure 10 .
Figure 10.Optimized operation results of Day2 without considering the user thermal comfort design scheme (α = 1.3).

Table 1 .
The questionnaire survey results.

Table 2 .
Parameters of thermal comfort model for different users.

Table 3 .
Demand preference coefficients for single feature users.

Table 4 .
Thermal comfort demand preferences of users with multiple characteristics.

Table 5 .
Thermal comfort design scheme based on users' difference in thermal comfort.

Table 7 .
Statistical results of relevant operating parameters and indicators.

Table 7 .
Statistical results of relevant operating parameters and indicators.