Study on Dynamic Pricing Strategy for Industrial Power Users Considering Demand Response Differences in Master–Slave Game

: With the deepening of power market reform, further study on power trading mechanisms has become the core issue of power market study. The development stage of the industrial electricity market requires efﬁcient and ﬂexible pricing mechanisms. Currently available pricing strategies are inadequate for demand response management. Therefore, this paper provides an in-depth study of the pricing mechanism in the industrial electricity market in the context of electricity market reform. It proposes a demand–response-based dynamic pricing strategy for industrial parks. The method proposes a dynamic pricing strategy for demand-side response in industrial parks based on master–slave game by establishing an exogenous model of demand-side response and incentives. Compared with the existing strategies, the strategy is more efﬁcient and ﬂexible, and effectively improves the economic efﬁciency of power trading and load regulation. Firstly, an exogenous model of demand-side response and incentive is built to characterize the demand-side response cost. The method focuses more on describing the exogenous characteristics of user incentives and response quantities. It only needs to analyze the exogenous indicators and random errors in various typical scenarios. The description of user demand-side response is more efﬁcient. Secondly, a master–slave-game-based dynamic pricing strategy for industrial parks with demand-side response is proposed. The strategy is composed of a two-stage optimization. The primary regulation of customers is achieved by day-ahead time-of-use tariffs. The secondary regulation of customers is achieved by means of the same-day regulation of demand and purchase regarding clean electricity. The proposed two-stage price formation mechanism is more economical, more effective in load regulation, and improves the ﬂexibility of industrial pricing. Finally, a case study is conducted on an industrial power user in a park in Liaoning Province. The results show that the proposed method is signiﬁcantly better than existing methods in terms of improving the economic efﬁciency and load control effectiveness of the pricing strategy.


Introduction
As an important segment of local economic development, industrial parks contain many forms of energy use. There are many types and large quantities of energy used in industrial parks, and there is a wide scope for energy saving. Demand-side response resources are used for multiple types of loads in industrial parks combined with game theory to form a flexible electricity pricing strategy and a suitable pricing strategy. It is of great significance for flexible power consumption, efficient resource allocation, energy saving and carbon reductions in the park.
Demand-side response mechanisms improve the operating efficiency of the electricity market, and enhance the resource allocation efficiency and elasticity of electricity users. Therefore, a series of studies on demand-side response have been conducted by domestic from Game Theory expectations. Using an online questionnaire, the study discovers that the choice patterns in this context also deviate from Game Theory predictions. In Ref. [19], a non-cooperative game model with one master and many slaves is constructed. The model combines game theory in the transaction model and a mathematical model of the integrated source system. It greatly improves the utilization of system and user consumption surplus. In Ref. [20], a two-tier optimization model is developed based on demand-side response combined with coalition games. The same distribute algorithm is used to solve the system, which makes the system much more economical and resilient. In Ref. [21], the operational characteristics of the integrated energy sub-regional and hierarchical synergy are addressed. The information interaction between regions and between users and energy suppliers is considered. The integrated energy system is collaboratively optimized. Ref. [22] presents an equilibrium model for electricity retailers, using mean-variance utility theory to account for their risk preferences. Through theoretical analysis and the Levenberg-Marquardt algorithm, the paper investigates the effects of wholesale price uncertainty and risk preferences on bidding strategies and Nash equilibrium outcomes. Game theory has made progress in electricity market research, but still needs to address issues such as deviations from actual behavior, lack of a comprehensive framework, and failure to consider market dynamics and uncertainty.
However, there are several problems with these studies. Most of them propose tariff decisions and pricing schemes by targeting the internal structure of customer loads and using price elasticity matrices, lacking specific studies on the external characteristics of demand-side response. The existing static pricing strategy applied to industrial park pricing is not flexible enough. There has been insufficient study different types of tariff strategies. Their effectiveness, feasibility and impact in the power system have not been fully studied.
This paper proposes three innovative contributions: (a) This paper fills a gap in the references by proposing an exogenous mode to consider random errors. The model obtains the demand-side response coefficient index of users in each case based on the identification and analysis of historical data. This index is used to build an exogenous model of demand-side response and incentive to characterize the demand-side response cost. It provides a clearer and more effective characterization of demand-side response costs and the relevance of traditional industrial pricing. (b) A strategy is proposed for the first time in this paper considering a two-stage masterslave game of user demand-side response. The first stage of the strategy is to form a time-of-day tariff for the day ahead and to regulate customers once. The second stage is an intraday correction of the tariff based on a model of the external characteristics of the customer's demand-side response and the clean consumption incentive. Traditional pricing methods usually lack a specific consideration of user-demand response, and pricing strategies are too rigid to meet the flexibility needs of industrial users. Compared with traditional pricing methods, this innovative approach effectively meets the flexibility needs of industrial users and improves load regulation efficiency and cost-effectiveness. It improves the problem of traditional pricing strategies that failed to meet the flexibility needs of industrial users. The proposed pricing strategy is more economical and more effective in regulating the load. (c) This paper proposes introducing clean consumption subsidies and environmental cost factors into the incentive model. Through the adjustment of incentive indexes in various scenarios in the model, dynamic pricing is realized in the game to promote the consumption of renewable energy. This enables industrial park users to actively consume clean electricity and makes the effect of peak shaving and valley filling more significant. Traditional methods fail to adequately consider incentives for clean energy consumption. This leads to a low utilization rate for renewable energy. This paper incorporates clean consumption incentives and environmental cost factors into the incentive model. This encourages industrial park users to actively consume clean The remainder of this paper is organized as follows. Section 2 builds an exogenous model of demand-side response and incentive to characterize the demand-side response cost. The method focuses more on describing the exogenous characteristics of user incentives and response quantities. Only the exogenous indicators and random errors in various typical scenarios need to be analyzed. The description of user demand-side response is more efficient. Section 3 presents a master-slave game dynamic pricing strategy that considers user demand response and is tailored to industrial parks. Section 4 analyzes a case study of industrial power users in a park in Liaoning. Finally, the paper is concluded in Section 5.

Demand-Response Strategy for Industrial Parks
Industrial parks, as an important sector of local economic development, have a high demand for electricity and are more sensitive to electricity prices. Therefore, they tend to shift their loads based on electricity prices in order to reduce electricity costs without affecting production activities.
This paper uses demand-side response metrics to characterize the relationship between users, demand-side response volumes and costs. The external characteristics of industrial parks in various time periods and situations are modeled without analyzing the complex load composition inside users and building demand-response costs for various loads. Therefore, this model has the characteristics of high efficiency and clear characterization and reflects the allocation of demand-side resources. Through price incentives, a more efficient and flexible demand-side management can be established. The implementation framework of a demand-response strategy for industrial parks is shown in Figure 1.
tives and response quantities. Only the exogenous indicators and random errors in ous typical scenarios need to be analyzed. The description of user demand-side resp is more efficient. Section 3 presents a master-slave game dynamic pricing strategy considers user demand response and is tailored to industrial parks. Section 4 analy case study of industrial power users in a park in Liaoning. Finally, the paper is concl in Section 5.

Demand-Response Strategy for Industrial Parks
Industrial parks, as an important sector of local economic development, have a demand for electricity and are more sensitive to electricity prices. Therefore, they te shift their loads based on electricity prices in order to reduce electricity costs witho fecting production activities.
This paper uses demand-side response metrics to characterize the relationshi tween users, demand-side response volumes and costs. The external characteristics dustrial parks in various time periods and situations are modeled without analyzin complex load composition inside users and building demand-response costs for va loads. Therefore, this model has the characteristics of high efficiency and clear chara zation and reflects the allocation of demand-side resources. Through price incentiv more efficient and flexible demand-side management can be established. The imple tation framework of a demand-response strategy for industrial parks is shown in F 1.

Analysis of Demand-Response External Characteristics of Users
From the analysis, it is clear that the participation of users in regulation tasks dustrial parks is accompanied by certain losses. The size of the response of the laun user is proportional to the response cost.  From the analysis, it is clear that the participation of users in regulation tasks in industrial parks is accompanied by certain losses. The size of the response of the launched user is proportional to the response cost.
As users increase the response volume, the cost of participating in demand-side response grows faster. The relationship is approximated as a quadratic functional relationship and is approximately expressed by the following Formula (2): where S e is the user's response volume, C(S e ) is the cost of participating in the demandside response, and α cos t and β cos t are the user's demand-side response parameters. In general, the amount of response expected from a user is S e , then the incentive given to user I needs to be greater than or equal to its response cost. This means that users will spontaneously participate in a response only if the loss of participation is comparable to the reward received.
Considering the fluctuations and deviations in user behavior, as well as errors in the critical point due to constraints, there are two main sources of deviation in user incentives and responses: the first type of error comes from the variability of user behavior; the second type of error causes a deviation in the theoretical critical point. Taking into account the response error and constraints of users, Formula (1) is adjusted as follows: where C re is the user's response cost, L is the user's current load, ν is the user's response error. α re and β re are the user's demand response parameters, and their values depend on the user's response characteristics. As shown in Equation (4), the elasticity coefficient is related to the user's demand response and incentive. Parameters α re and β re are used as demand-side response indicators to characterize user demand-side response elasticity.
In Formula (1), when the reward obtained by the user equals the loss incurred, the user's objective is to participate in demand response with the lowest cost and maximum economic benefit. As the effect of lower order errors is relatively negligible, lower-order errors are omitted from the analysis of the user's response elasticity. Substituting Formula (3) into Formula (4), the user's response elasticity is obtained as follows: where ε S e is the response elasticity of the user. Different users have different response elasticity, as shown in Figure 2.  The response behavior of users results in the loss of their own costs. As the degree o user participation in response increases, their losses increase at a faster rate. Therefore, a cost function for user participation in response is proposed as an approximation based on these characteristics: S is the response of user i to the re ceived stimulus in time period t . i  and i  are the user's demand-side response co efficients. The incentive factor is determined by fitting historical data for the user. In this model, once the demand-side response coefficient of the user for that time period is ob tained, its response cost, which is the minimum incentive required, is obtained. The for mation of this cost is no greater than the average electricity price in the market av q , which gives rise to the following constraint:  From Formula (3), it is evident that parameters α re and β re reflect the differences in the elasticity of demand response for different users. They describe the elasticity of the demand response on the demand side. Since the response elasticity of the same user varies at different times, users have different values of α re and β re at different times of day. Identifying the demand response parameters during typical time periods suffices to obtain the demand response situation, as shown in Formula (6), for the load response.
The response behavior of users results in the loss of their own costs. As the degree of user participation in response increases, their losses increase at a faster rate. Therefore, a cost function for user participation in response is proposed as an approximation based on these characteristics: where C t i is the user's response cost function. S t ei is the response of user i to the received stimulus in time period t. α i and β i are the user's demand-side response coefficients. The incentive factor is determined by fitting historical data for the user. In this model, once the demand-side response coefficient of the user for that time period is obtained, its response cost, which is the minimum incentive required, is obtained. The formation of this cost is no greater than the average electricity price in the market q av , which gives rise to the following constraint:

Demand-Side Response Extrinsic Model Considering Random Errors
In Formula (8), due to the impact of external factors on users' behavior, the error term ν t i is added. The uncertainty of response behavior leads to the error term ν t i , which characterizes the noise caused by i users in time period t. This deviation is result from temperature, weather, and typical day types. The error parameters formed by applying various indicators make the proposed model more accurate in characterizing users. These parameters are found in Table 1, which is formed referring to the fuzzy formula in reference [23] to enhance the objectivity and accuracy of the indicators. The input of the fuzzy controller is the error e between the historical and current actual indicators and the rate of change in error e c . e c = e(t) − e(t − 1) The output is the correction factor α ∈ [0, 1] for the actual indicators, and the correction function is obtained as follows: Based on the fuzzy formula, the following indicators are obtained: The specific indicators for demand-side response errors are listed in Table 1. Therefore, the noise term ν t i is represented as: Moreover, since parameters α i and β i are also affected by various environmental factors, they also need to be adjusted accordingly.
The extrinsic model for user demand-side response and incentives is obtained as follows: where λ i is the error value due to other factors. The total benefit of user demand-side response in time period t is finally built as follows: where Y t i is the user's benefit, I 0 is the base incentive, γ is the incentive at the given time. S t ei−max is the maximum response quantity of the user during time period t.

Pricing Strategy Based on Master-Slave Game Theory
This paper proposes a dynamic pricing strategy for industrial parks based on masterslave game theory. The strategy aims to build a pricing mechanism for industrial parks under demand-side response incentives and environmental cost guidance that differ from the traditional master-slave game decision pricing. The proposed pricing strategy's objective function to achieve dynamic pricing in different scenarios through a change in controllable quantities is the innovation of this chapter. The strategy consists of two stages: day-ahead and intra-day price adjustments. In the day-ahead stage, the power retailer develops electricity price information for each time period of the next day based on historical data, load forecasting, and daily load curves. The price includes the following components: where χ t is the purchasing price of electricity, χ u is the cost of using the distribution network, and χ p t is the profit obtained by the power retailer. In the first stage, users adjust their electricity usage arrangements based on the dayahead price information. The day-ahead tariff is corrected by controlling the intra-day incentive change. The second-stage objective function is built to maximize the revenue of electricity sellers and minimize the cost of electricity for customers. Price incentive constraints and park load characteristic constraints are set, and variable intra-day electricity sales prices are solved to achieve dynamic pricing. The process is illustrated in Figure 3.
develops electricity price information for each time period of the next day based on his torical data, load forecasting, and daily load curves. The price includes the following com ponents: where t  is the purchasing price of electricity, u  is the cost of using the distribution network, and p t  is the profit obtained by the power retailer.
In the first stage, users adjust their electricity usage arrangements based on the day ahead price information. The day-ahead tariff is corrected by controlling the intra-day in centive change. The second-stage objective function is built to maximize the revenue o electricity sellers and minimize the cost of electricity for customers. Price incentive con straints and park load characteristic constraints are set, and variable intra-day electricity sales prices are solved to achieve dynamic pricing. The process is illustrated in Figure 3.   When there is renewable energy access, to ensure the park users actively consume clean electricity, a clean consumption subsidy will be developed. The environmental cos factor is introduced on the basis of the traditional renewable energy consumption allow ance. The carbon emissions of users are counted and an environmental cost factor i formed, and the renewable energy consumption allowance received by users is bounded When there is renewable energy access, to ensure the park users actively consume clean electricity, a clean consumption subsidy will be developed. The environmental cost factor is introduced on the basis of the traditional renewable energy consumption allowance. The carbon emissions of users are counted and an environmental cost factor is formed, and the renewable energy consumption allowance received by users is bounded by the environmental cost factor. This means customers consume renewable energy with lower carbon emissions.
For industrial power users, their environmental costs are calculated based on their "carbon footprint". The original environmental cost calculation included industrial waste gas, wastewater, and solid waste. However, this strategy only considers the following two forms of pollution: (a) Industrial waste gas emissions The pollution in the original industrial waste gas is calculated based on the emissions of sulfur dioxide and dust. To account for carbon emissions in decision-making, carbon dioxide emissions are also included in the calculation. The carbon footprint is calculated as follows: where r 1 is the balance factor of forestland in the area where the industrial park is located to absorb waste gas. c 1 , c 2 , and c 3 are the contents of sulfur dioxide, carbon dioxide and dust in the waste gas, respectively. a 1 is the absorption capacity of vegetation for sulfur dioxide, carbon dioxide, and dust. The pollution in industrial wastewater is mainly calculated based on sulfur dioxide and chemical oxygen demand. Carbon emissions are also taken into account in the calculation. The carbon footprint is calculated as follows: where r 2 is the balance factor of wetlands in the area in which the industrial park is located for absorbing wastewater. c 4 , c 5 , and c 6 are the contents of sulfur dioxide, chemical oxygen demand, carbon in the wastewater, respectively. a 2 is the absorption capacity of wetlands for sulfur dioxide, chemical oxygen demand, and carbon.

Demand-Response Incentive
The demand-response incentive formulated by the electricity supplier needs to ensure the required response amount and minimize its own cost. Therefore, according to the analysis of user demand-response cost and incentive in Section 4. The demand-response cost for the electricity supplier is the minimum incentive required for user response in that time period. The incentive required for user i in time period t is given by Formula (20). The formation of this model is based on the external characteristics of demand-response and incentives under different scenarios. The incentive varies with time and application scenarios.

Incentive Model for Electricity Retailers
The incentive model of the power retailer is based on different scenarios, divided into two situations: with or without the integration of clean energy. Different incentives are set according to the response demand and the integration of clean energy. The incentive model for the power retailer when purchasing clean energy is as follows: The incentive model for the power retailer when there is no clean energy access is:

Objective Function
The revenue objective function of the power retailer is: , there is no clean energy access where q t , q c , and q e are the selling price, purchasing price, and clean energy price, respectively. ω is the proportion of clean energy purchases. L(t) is the load at time t, π t i is the cost of incentives for user response and clean energy consumption, which varies with different application times and scenarios as a controllable variable. f C is the risk cost of the electricity retailer. L 1 (t) is the electricity consumption before regulation at time t. ∆L(t) is the electricity consumption adjusted for the load at time t.
Since the clean energy purchase cost in the revenue function of the power supplier is uncertain, the revenue function of the power supplier is also a random variable. Therefore, according to probability theory, the mean and standard deviation of the revenue function of the power supplier is obtained as follows: The upper-level model is built as shown in Formula (28):

Constraint Conditions
The constraints of the upper-level function are the upper and lower bounds of incentive prices. π min t is the minimum required incentive for user participation in demand-response at time t. π max t is the maximum required incentive for user participation in demandresponse and the consumption of renewable energy, which is lower than the average electricity price level.
For pricing, it is necessary to avoid the situation where the seller sets too high a price to ensure their profit. Therefore, the average dynamic electricity price of the seller needs to be within a certain limit. Its average value needs to be less than or equal to the market-average electricity price.

User Lower-Level Benefit Function Objective Function
The objective function is to minimize the user's electricity cost under the given conditions, where Y t i is the benefit that user i obtains from participating in demand-side response: The constraint for user participation in demand response is: where τ t i is the operating state of user i at time t, with τ t i = 1 indicating that the user i starts working at time t; otherwise, it stops working. DT i is the working cycle of user i, and k t i are the response status of user i at time t.
The power balance constraint is defined as: where L base t and L re t are the base load and the responsive load.

Game Theory Model for Power Retailers and Industrial Power Users
The proposed master-slave game model regards the power company as the leader and the users as the followers. Different incentive models are developed for different scenarios in the model, which is divided into two cases: with and without access to clean energy. Different incentives are formulated based on response demands and the use of clean energy, enabling dynamic pricing. Moreover, the incentive model includes clean energy consumption subsidies and environmental cost factors. It encourages users to reduce carbon emissions through the adoption of new production technologies and improved management. As followers, user electricity consumption behavior depends on their electricity cost model and usage limits.
The power-seller's decision set is determined by the user's electricity consumption.
where L t i is the electricity consumption of user i at time t, w t i is the revenue of the sales company for the i-th user at time t. The objective functions are (29) and (30), respectively.
The master-slave game model is as follows: where I is the set of users and θ i is the set of user decisions, namely electricity consumption. θ s is the set of electricity seller decisions, namely electricity prices.
Assuming that the solutions for the upper-and lower-level models and their equilibrium solutions are a, b, and a * , b * , respectively, the game reaches equilibrium when the following conditions are met: The specific game process is as follows: Firstly, the electricity seller, acting as the leader, sends rewards and demand responses to users, initiating the game. Secondly, users select the best countermeasures through the objective function. The electricity seller calculates the optimal price based on user responses and the objective function and sends this price back to users. Finally, this process is repeated until the final result is obtained, and the resulting electricity price is used to update the day-ahead price.

Case Study Analysis
In this paper, data for an industrial customer in a park in Liaoning province are selected. The demand-side response parameters are solved and modeled according to the load excitation-response curves of typical days under various scenarios of historical response data. The dynamic tariffs of two scenarios without and with clean electricity are calculated and analyzed.

Simulation Analysis of User Demand-Response Indicators
As shown in Table 2, historical load incentive-response consumption data provide a simplified exploration and understanding of the response characteristics of industrial loads under different scenarios. They are also effective for the construction of corresponding demand-side response models and dynamic tariff analysis. Users actively participate in demand-side response when the incentive they receive is greater than or equal to their own cost of participating in demand-side response. Therefore, it is assumed that the incentive received by the load in the historical response data is its response cost. The analytical solution for the demand-side response index is shown in Figure 4. Users actively participate in demand-side response when the incentive they receive is greater than or equal to their own cost of participating in demand-side response. Therefore, it is assumed that the incentive received by the load in the historical response data is its response cost. The analytical solution for the demand-side response index is shown in Figure 4.  . Therefore, the demand-side cost model of user 1 for that time period is obtained as: By introducing a set of random noise simulations to the user response data, the modified demand-side response external characteristic model is expressed as: The incentive results are shown in Table 3. From Figure 3, it can be concluded that the average error of this model is only 8.1%, and it performs well in these applications.  The demand-side response indicator of user 1 for that time period is obtained by the data-fitting calculation α t 1 = 2.70, β t 1 = 3.08. Therefore, the demand-side cost model of user 1 for that time period is obtained as: By introducing a set of random noise simulations to the user response data, the modified demand-side response external characteristic model is expressed as: The incentive results are shown in Table 3. From Figure 3, it can be concluded that the average error of this model is only 8.1%, and it performs well in these applications. As shown in Table 4, presents the incentive and response consumption data in Scenario 1 cases. In this scenario, the demand-response coefficient is α t 1 = 3.78, β t 1 = 2.96. The results obtained when applying this parameter to build the model and performing analysis and calculation are shown in Figure 5. As shown in Table 4,presents the incentive and response consumption data in Scenario 1 cases. . The results obtained when applying this parameter to build the model and performing analysis and calculation are shown in Figure 5. As shown in Table 5,presents the incentive and response consumption data in Scenario 2 cases.  As shown in Table 5, presents the incentive and response consumption data in Scenario 2 cases. In this scenario, the demand-response coefficient is α t 1 = 3.56, β t 1 = 2.38. Applying this parameter to build the model and performing analysis and calculation, the results are shown in Figure 6. In this scenario, the demand-response coefficient is Applying this parameter to build the model and performing analysis and calculation, the results are shown in Figure 6. Through the above data analysis, the errors obtained under the two scenarios are shown in Figure 7. These errors perform well in calculations.  Through the above data analysis, the errors obtained under the two scenarios are shown in Figure 7. These errors perform well in calculations. In this scenario, the demand-response coefficient is ing this parameter to build the model and performing analysis and calculation, the res are shown in Figure 6. Through the above data analysis, the errors obtained under the two scenarios shown in Figure 7. These errors perform well in calculations. The traditional model does not take into account noise effects from various facto user demand-side response, nor does it consider random errors from users. To com the average model errors, 30 test points were chosen under the same scenario, conside and not considering random errors. The results are presented in Figure 8. The traditional model does not take into account noise effects from various factors in user demand-side response, nor does it consider random errors from users. To compare the average model errors, 30 test points were chosen under the same scenario, considering and not considering random errors. The results are presented in Figure 8. The results show that the average error of the model that considers noise facto reduced. Therefore, the proposed model in this paper has a clearer and more efficien scription of the load's response to the user's demand side. It is suitable for building mand-side response models for multiple types of users in industrial parks.

Analysis of Dynamic Electricity Prices under Two Scenarios
By analyzing the typical historical daily load incentive response volume curve customers, the demand-side response volume can be obtained by the second stag lower optimization. The required incentives are obtained as follows. Table 6 shows second-stage optimization response volume and required excitation.  The results show that the average error of the model that considers noise factors is reduced. Therefore, the proposed model in this paper has a clearer and more efficient description of the load's response to the user's demand side. It is suitable for building demand-side response models for multiple types of users in industrial parks.

Analysis of Dynamic Electricity Prices under Two Scenarios
By analyzing the typical historical daily load incentive response volume curves for customers, the demand-side response volume can be obtained by the second stage of lower optimization. The required incentives are obtained as follows. Table 6 shows the second-stage optimization response volume and required excitation. It is possible to analyze the tariff in two scenarios: when there is renewable energy access and when there is no renewable energy access. The resulting two-stage dynamic tariff is as follows: (a) Scenario 1: No clean electricity The load curve optimized by this strategy is shown in the Figure 9. As shown in Figure 9, through the two-stage electricity price guidance, the electricity seller can formulate prices that are closer to the needs of users and their own needs. As shown in Figure 10, the two-stage pricing strategy has a more significant effect on load transfer, peak shaving and valley filling.
It is possible to analyze the tariff in two scenarios: when there is renewable en access and when there is no renewable energy access. The resulting two-stage dyna tariff is as follows: (a) Scenario 1: No clean electricity The load curve optimized by this strategy is shown in the Figure 9. As shown in ure 9, through the two-stage electricity price guidance, the electricity seller can formu prices that are closer to the needs of users and their own needs. As shown in Figur the two-stage pricing strategy has a more significant effect on load transfer, peak sha and valley filling. When clean energy (wind and solar power, for example) is integrated, the purcha cost of the electricity seller decreases. Corresponding consumption subsidies are prov to users. The optimized electricity prices formed are as follows: In this scenario, the load curve of the user is shown in Figure 11. This shows that integrating clean energy, the purchasing cost of the electricity seller decreases and the receives a clean consumption subsidy. Therefore, the user is more actively involved in response and consumption of clean energy, resulting in a more significant effect of p shaving and valley filling. It is possible to analyze the tariff in two scenarios: when there is renewable ene access and when there is no renewable energy access. The resulting two-stage dyna tariff is as follows: (a) Scenario 1: No clean electricity The load curve optimized by this strategy is shown in the Figure 9. As shown in ure 9, through the two-stage electricity price guidance, the electricity seller can formu prices that are closer to the needs of users and their own needs. As shown in Figure the two-stage pricing strategy has a more significant effect on load transfer, peak shav and valley filling. When clean energy (wind and solar power, for example) is integrated, the purchas cost of the electricity seller decreases. Corresponding consumption subsidies are provi to users. The optimized electricity prices formed are as follows: In this scenario, the load curve of the user is shown in Figure 11. This shows that a integrating clean energy, the purchasing cost of the electricity seller decreases and the u receives a clean consumption subsidy. Therefore, the user is more actively involved in response and consumption of clean energy, resulting in a more significant effect of p shaving and valley filling. In this scenario, the load curve of the user is shown in Figure 11. This shows that after integrating clean energy, the purchasing cost of the electricity seller decreases and the user receives a clean consumption subsidy. Therefore, the user is more actively involved in the response and consumption of clean energy, resulting in a more significant effect of peak shaving and valley filling.

Analysis and Comparison of Dynamic Electricity Prices under Two Scenarios
The response of the load under the two-stage electricity price optimization strat with and without clean energy integration, was analyzed. As shown in Figures 12 an in the scenario with clean electricity access, due to the significant reduction in electr costs for customers and the incentive of clean consumption subsidies, users are mor tively involved in the response. This further strengthens the peak-shaving and valley ing situation.  Table 7 compares the economic benefits of the first-stage electricity price forma and the two-stage electricity price formation under the two scenarios.

Analysis and Comparison of Dynamic Electricity Prices under Two Scenarios
The response of the load under the two-stage electricity price optimization strategy, with and without clean energy integration, was analyzed. As shown in Figures 12 and 13, in the scenario with clean electricity access, due to the significant reduction in electricity costs for customers and the incentive of clean consumption subsidies, users are more actively involved in the response. This further strengthens the peak-shaving and valley-filling situation.  Figure 11. Electricity price of an industrial power user in a certain park with clean energy integ tion.

Analysis and Comparison of Dynamic Electricity Prices under Two Scenarios
The response of the load under the two-stage electricity price optimization strateg with and without clean energy integration, was analyzed. As shown in Figures 12 and in the scenario with clean electricity access, due to the significant reduction in electric costs for customers and the incentive of clean consumption subsidies, users are more tively involved in the response. This further strengthens the peak-shaving and valley-fi ing situation.  Table 7 compares the economic benefits of the first-stage electricity price formati and the two-stage electricity price formation under the two scenarios.    Figure 11. Electricity price of an industrial power user in a certain park with clean energy integ tion.

Analysis and Comparison of Dynamic Electricity Prices under Two Scenarios
The response of the load under the two-stage electricity price optimization strate with and without clean energy integration, was analyzed. As shown in Figures 12 and in the scenario with clean electricity access, due to the significant reduction in electric costs for customers and the incentive of clean consumption subsidies, users are more tively involved in the response. This further strengthens the peak-shaving and valley-fi ing situation.  Table 7 compares the economic benefits of the first-stage electricity price format and the two-stage electricity price formation under the two scenarios.   Table 7 compares the economic benefits of the first-stage electricity price formation and the two-stage electricity price formation under the two scenarios. As shown in Table 7, the two-stage optimization significantly reduces the purchasing risk cost for users and increases the demand response cost. After clean energy integration, the purchasing risk for clean energy increases, but the purchasing cost decreases. Compared to the scenario without clean energy integration, the profit of the electricity seller increases by 21.46% and the user cost decreases by 8.47% in this scenario. Furthermore, the peakto-valley ratio under the two-stage electricity price is reduced by 24.5% compared to the first-stage electricity price. The peak-to-valley ratio under the scenario with clean energy integration is reduced by 26.83% compared to the scenario without clean energy integration. It can be seen that the two-stage electricity price-formation strategy has a more significant effect on peak-shaving and valley-filling. Therefore, the proposed pricing strategy is more economically efficient and more effective in terms of load control.

Conclusions
This paper proposes a demand-response-based master-slave game dynamic pricing strategy for industrial parks. The following conclusions are drawn from the simulation analysis: (a) A demand-side response externality model is built. It accurately describes the incentive response characteristics of industrial power users with an average error of about 17%, which is an 8.71% reduction compared to models without considering random errors. (b) The proposed model is efficient in describing demand-side response for various types of users in industrial parks. This approach offers a clearer representation of load and user demand-side response characteristics. (c) The proposed demand-response-based master-slave game dynamic pricing strategy for industrial parks shows enhanced economic efficiency. There is a 24.5% reduction in peak-to-peak difference using the two-stage tariff compared to traditional time-of-use pricing. (d) The load curve analysis reveals a 26.83% reduction in peak-to-peak difference when clean power is purchased. This indicates that clean power purchasing contributes to further peak-and-valley reductions, enhancing the effectiveness of the proposed pricing strategy for load regulation.
Future work is planned to take more error factors into account. The constraint relationship between the two will be further explored to bring the model closer to the displayed situation.  Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.
Data Availability Statement: Data will be made available on request.

Conflicts of Interest:
We declare that we do not have any commercial or associative interests that represent a conflict of interest in connection with the submitted work. Profit obtained by the power retailer r 1 Balance factor of forestland in the area c 1 Contents of sulfur dioxide c 2 Carbon dioxide c 3 Dust in the waste gas a 1 Absorption capacity of vegetation for sulfur dioxide, carbon dioxide, and dust r 2 Balance factor of wetlands in the area c 4 Contents of sulfur dioxide c 5 Chemical oxygen demand c 6 Carbon in the wastewater, respectively a 2 Absorption capacity of wetlands for sulfur dioxide, chemical oxygen demand, and carbon q t Selling price q c Purchasing price q e Clean energy price, respectively.

Abbreviations
ω Proportion of clean energy purchasing L(t) Load at time t The cost of incentives f C Risk cost of the electricity retailer L 1 (t) Electricity consumption before regulation at time t ∆L(t) Electricity consumption adjusted for the load at time t Lower-level models a * Upper-level equilibrium solutions b * Lower-level equilibrium solutions