A Credibility Theory-Based Robust Optimization Model to Hedge Price Uncertainty of DSO with Multiple Transactions

: This paper addresses the deregulated electricity market arising in a distribution system with an electricity transaction. Under such an environment, the distribution system operator (DSO) with a distributed generator faces the challenge of electricity price uncertainty in a spot market. In this context, a credibility theory-based robust optimization model with multiple transactions is established to hedge the uncertain spot price of the DSO. Firstly, on the basis of credibility theory, the spot price is taken as a fuzzy variable and a risk aversion-based fuzzy opportunity constraint is proposed. Then, to exploit the resiliency of multiple transactions on hedging against uncertain spot price, the spot market, option contract and bilateral contract integrating power ﬂow constraints are studied, because it is imperative for DSO to consider the operational constraints of the local network in the electricity market. Finally, the clear equivalence class is adopted to transform the risk aversion constraint into a deterministic robust optimization one. Under the premise of considering the expected cost of the DSO, the optimal electricity transaction strategy that maximizes resistance to uncertain spot price is pursued. The rationality and effectiveness of the model are veriﬁed with a modiﬁed 15-node network. The results show that the introduction of option contracts and bilateral contracts reduces the electricity transaction cost of DSO by USD 28.5. In addition, under the same risk aversion factor, the cost of the proposed model is reduced by USD 195.18 compared with robust optimization, which avoids the over-conservatism of traditional robust optimization.


Introduction
The distribution system operator (DSO) is responsible for maintaining the security of supply and power quality through investment, construction and reconfiguration of the existing distribution system [1][2][3].With the deepening reform of the electricity market environment, as a stakeholder, DSO with distributed generator (DG) plays an important role in the electricity transaction of the distribution system.That is, DSO purchases electricity in the upper wholesale market to meet customer demand as well as maximization of its utility [4][5][6].However, the price in the spot market is characterized by uncertainty due to fluctuations in electricity demand, fuel price and renewable power generation [7,8].Moreover, the forecast error of the spot price is inevitable [9].Thus, in order to obtain the optimal electricity transaction and expected utility, DSO has to capture the uncertain spot price from the perspective of risk aversion [10].In view of this, two main questions need to be answered: how to use a portfolio of electricity purchase transactions to hedge against risk brought by uncertain spot price and how to assess the risk when formulating an optimal electricity transaction strategy under the premise of the expected cost?

Literature Review
For risk decision problems with uncertain electricity price, researchers mainly use three kinds of optimization, including robust optimization [11][12][13][14][15], stochastic optimiza-tion [16][17][18][19][20] and fuzzy optimization [21][22][23][24][25][26].For instance, in [11], the uncertainty of selling/purchasing price in an electricity market is handled by robust optimization with a polyhedral uncertain set.Accordingly, an economical optimal solution is obtained in consideration of the undesired deviation of the market electricity price from the forecasted one.In [12], an adaptive robust optimization is developed to study the uncertain price in a real-time market.In [13], a maximum-minimum-maximum robust optimization model considering the price deviation in the electricity market is proposed, which improves the robustness of system operation against forecast uncertainty.In [14], the uncertainty of market price is dealt with by the upper deviation of forecast price and the robust electricity trading strategies of risk neutrality and risk aversion are compared.In [15], the uncertainty of the electricity market price is modeled based on robust optimization.In this model, instead of the predicted electricity price, the maximum and minimum amounts of the electricity price are considered.However, robust optimization mainly focuses on the worst-case of an uncertain problem and does not fully investigate the risk preference characteristic of a decision maker [27].Stochastic optimization is one of the most commonly used methods for managing uncertain price.In [16], by assuming that the uncertain price follows the normal distribution, a scenario generation-based stochastic framework is developed.In this framework, the risk associated with uncertain electricity price is considered through downside risk constraints.In [17], uncertain electricity price is regarded as a random variable and a stochastic optimization model based on the Monte Carlo sampling method is established in microgrid (MG) optimal operation.In [18], a multi-stage stochastic programming method is developed.The bidding strategy in the spot market is described as a Markov decision problem and solved by approximate dual dynamic programming.In [19], the decision-making problem of a retailer under uncertainty is discussed based on stochastic optimization.In [20], the random scenario method is derived to simulate the uncertain spot price and the conditional value at risk is proposed to evaluate the risk of the electricity trading strategy.However, stochastic optimization requires repeated sampling and the solution efficiency is reduced [28].
In addition, the work in [29] states that besides the random feature, the uncertainty also includes the fuzzy feature.Thus, fuzzy optimization, including fuzzy rough set [21], image fuzzy set [22], neutral particle set [23] and so on, has been studied in the operation of a power system over the past few years.In [24], the authors consider the fuzzy feature of uncertain electricity price, in which the fuzzy feature is approximated with a fuzzy number.In [25], a fuzzy set theory-based MG energy management model is established for price uncertainty.In this model, the uncertain electricity price is characterized by triangular fuzzy numbers.In [26], a risk measurement method based on credibility theory is proposed to evaluate the fuzziness of uncertain wind power.However, few studies have addressed the robust power trading model based on credibility theory in view of the uncertainty of spot price in the electricity market.
On the other hand, the existing research on electricity transaction mainly focuses on the discussion of a return model and the formulation of transaction strategy.For example, In [30], the transaction strategy of the power retailer in the spot market is analyzed and the profit model of the retailer under the background of new power reform is discussed.In [31], a deterministic multi-objective optimization model with the goal of profit maximization and peak demand minimization is established to study the short-term decision-making problem of the retailer.In [32], the optimal bidding strategy of an energy hub in the power market is studied under the protection of energy network information privacy.These studies do not take into account the risk assessment of retailers in electricity trading.In [33], a risk decision-making model of electricity transaction with the goal of profit maximization is established.The model analyzes the impact of different electricity transaction combinations on profits from the retailer's point of view.In [34], a risk management model for DSO portfolio with multiple electricity purchase markets is constructed.The impact of different risk appetites on transaction strategy is studied.In [35], the retailer tries to hedge against uncertainty through three trading platforms.The electricity price uncertainty is modeled with the auto regressive integrated moving average method and the retailer's electricity transaction strategy is determined.In [36], based on the portfolio optimization theory, the optimization model of electricity purchase and sale portfolio is constructed to explore the influence of different factors on the purchase and sale risk in the multi-level electricity market.However, in most cases, DSO is not only responsible for trading electricity but also should consider the operational constraints of the distribution system.Therefore, it is necessary to consider network topology as well as power flow constraints of the distribution system when selecting electricity purchase transactions.
Table 1 reports the majority of the studies presented within the last decade; however, most of the existing electricity purchase strategies do not take into account the topological constraints of the network.In addition, to the best of the authors' knowledge of this paper, there are no studies in the literature addressing the robust electricity trading model based on credibility theory in view of the uncertainty of spot price in the electricity market.

Our Contributions
To address the above issues, this paper develops a risk aversion DSO electricity transaction model based on the credibility theory.The proposed model can help decision makers determine the optimal combination of electricity purchase transactions under an acceptable risk level, considering the uncertain electricity price and power flow constraints when formulating electricity transaction strategy.The main contributions of this paper are as follows:

•
Based on credibility theory, a risk aversion-based fuzzy chance constraint model is proposed.In the model, the uncertain spot price is designed as a fuzzy variable and its credibility distribution is derived to assess the uncertain risk.The proposed model optimizes the credibility that the expected objective is met, from which decision makers can assess the risk of transaction strategy.

•
Multiple transactions, including the spot market, option contract and bilateral contract, are considered to hedge the risk caused by uncertain price, and the impact of different electricity transaction combinations on DSO cost is analyzed while considering power flow constraints.

•
A clear equivalence class method with fuzzy chance constraint is used to transform the proposed model into a deterministic robust optimization model.The effectiveness of the model is verified with a modified 15-node network.

Organization of the Research
The rest of this paper is organized as follows.The credibility function associated with forecast error percentage of spot price is derived in Section 2. The multiple electricity transactions model is established in Section 3. A credibility theory-based robust optimization model to hedge uncertain spot price of DSO with multiple transactions is proposed in Section 4. Case studies and related analysis are introduced in Section 5. Finally, this paper concludes in Section 6.

Problem Formulation
Fuzzy decision-making is a kind of method to solve problems with fuzzy nature, but the traditional fuzzy decision-making has not established a complete axiomatic system.This leads to unconvincing decision-making conclusions until the credibility theory-based uncertainty measurement is established.It makes up for the disadvantage that possibility measure does not have self-duality [37] and provides a new tool for scholars to study fuzzy decision-making problems.In credibility theory, the credibility measure is developed to describe the credibility of fuzzy events [38].It holds that events with credibility 1 must occur and events with credibility 0 do not occur, which avoids the decision-making confusion that may be caused by the traditional calculation of membership degree.
The credibility measure can be expressed by the minimum supremum of variable in a fuzzy event set.For any set A ∈ R, the credibility measure of fuzzy variable ξ ∈ A is defined as [39]: where sup x∈A µ(x) and 1 − sup µ(x) denote the possibility measure and necessity measure of A, respectively.A c represents the complement of the set A, and µ is the membership function of the fuzzy variable.The average value of the possibility measure and the necessity measure in Equation ( 1) is used to ensure the establishment of duality.In addition, the credibility measure satisfies the following four axioms: In this paper, the uncertain spot price is designed as a fuzzy variable and the credibility distribution function including possibility measure and necessity measure is derived to evaluate the uncertain risk.The uncertain risk measurement model based on credibility theory also satisfies these four axioms.

Credibility Distribution Function Associated with Forecast Error Percentage of Spot Price
In electricity transactions, there are inevitable errors in the forecast of spot price [40].Assume that the forecast error percentage of spot price is ε and the mathematical expression is as follows: where λ sm t and λ sm t are the actual spot price and the forecast spot price, respectively.The membership function µ associated with forecast error percentage of spot price can be expressed as the Cauchy distribution [24].The uncertain spot price is taken as the fuzzy variable and its mathematical expression can be described as: where E + and E − , respectively, represent the statistical average of positive and negative error percentages and ω is the weighting factor.
After derivation, we can obtain the credibility function of ε: Proof.According to Equation ( 1), for ε ∈ R, the mathematical expression of the credibility measure is If ε > 0, we have and Combining Equations ( 6) and ( 7), if ε > 0, we have If ε ≤ 0, we have and Combining Equations ( 9) and ( 10), if ε ≤ 0, we have This completes the proof.
The credibility and membership functions associated with forecast error percentage of spot price are shown in Figure 1, where

Objective Function
DSO conducts electricity transaction through the spot market, bilateral contract and option contract.Consider T hour periods and T 1 and T 2 to be, respectively, the peak and non-peak period sets of power demand, satisfying T 1 + T 2 = T.Call option contract is only for T 1 and trading volume does not change over time.In addition, suppose there are N bilateral contracts for DSO to choose and the details of the electricity transaction cost are as follows.
(1).The mathematical expression of DSO's cost function C sm in spot market is as follows: where b ∈ B is the range of network node.λ sm t and p sm t,b are the electricity price and trading volume of DSO in spot market during time t, respectively.
(2).The cost function C bc of DSO in the bilateral contract electricity transaction can be expressed as: where λ bc n is the electricity price with respect to bilateral contract n and p bc t,b,n is the trading volume of nth bilateral contract selected by node b during time t.
(3).The cost function C oc of DSO from the option contract electricity transaction can be calculated as: where p oc t,b is the call option contract trading volume of node b during time t, λ ck and λ 0 are the strike price and premium of the call option, respectively.If λ sm t > λ ck , DSO executes the option contract and its option contract purchases electricity at the fixed price; if λ sm t < λ ck , DSO abandons the exercise option and its option contract purchases electricity at the spot market price.
(4).The power generation cost function C dg of DSO can be expressed as: where p t,b is the active power output of DG at node b during time t and λ dg is the power generation cost price of DG.The mathematical expression of the cost C of the DSO with the spot market, option contract and bilateral contract in an electricity transaction is as follows:

Constraints
In order to ensure that the system operates in a safe and reliable environment, the electricity transaction must meet the following constraints.
(1).Active power output constraint of DG: where P max t,b is the maximum active power output of DG at node b during time t.(2).Node voltage constraint: where v t,b is the square of the voltage at node b during time t.v min t,b and v max t,b are the maximum and minimum values of node voltage at node b during time t, respectively.
(3).Contract volume constraint for bilateral contract: where s b,n is a binary variable.If node b selects contract n, then s b,n = 1; otherwise, s b,n = 0. p bc,min n , p bc,max n are the minimum and maximum contract volumes of bilateral contract n, respectively.
(4).During time t, the total amount of electricity purchased in spot market, bilateral contract and option contract of DSO equals the amount of active power injected from the power grid.Its mathematical expression can be described as: where P grid t,b is the amount of active power injected from the power grid at node b during time t.
(5).The power flow constraints of the distribution network are as follows: where l ∈ L is the range of network line.f p t,l , f q t,l are the active and reactive power flow of line l during time t, respectively.D p t,b and D q t,b are the active load and reactive load of node b during time t, respectively.a t,l is the square of the current of line l during time t.Q grid t,b is the reactive power injected from the power grid at node b during time t.S 2 t,l is the upper limit of the apparent power of line l during time t, and R l , X l , G b and B b are the parameters of resistance, reactance, admittance and conductance of distribution network, respectively.s(l) is the power outflow end of line l and r(l) is the power inflow end of line l.The balance constraints of active and reactive power are shown in Equations ( 22) and (23).Equation (24) relates the line flow to the node voltage.Equation (25) represents the apparent power flow limitation of each line transmitting node and Equation ( 26) is a quadratic curve constraint, which convexes the original non-convex AC OPF problem [41].Under quite unrestricted assumptions, the rationality of this convexity is proved in [42].Equation (27) represents the apparent power flow limitation of each line receiving node.

Robust Optimization Model for DSO Based on Credibility Theory
The forecast error of spot price is inevitable [43].DSO with different risk preferences needs to hedge the risk caused by forecast error while considering operational cost as well as power flow constraints.Given a certain electricity purchase cost, DSO pursues an electricity transaction strategy that maximizes resistance to the uncertain spot price.In view of this, this paper establishes a credibility theory-based robust optimization model to hedge price uncertainty of DSO with multiple transactions.max |ε| (28a) ) where C 0 is the minimum cost of DSO when the spot price equals the forecasted spot price.σ is the risk aversion factor, which indicates the DSO's aversion to the risk due to the uncertain spot price.α is the credibility index and the physical meaning is equivalent to the probability confidence.Equation (28b) is expressed as the credibility that the actual cost of DSO less than the expected cost is not less than α.Equation (28c) represents the expected cost of DSO.When σ is larger, expected cost C e is higher, indicating that DSO has a greater degree of risk aversion.q is the decision variable, which represents the amount of electricity traded by DSO in each market.Generally, when the actual spot price takes the maximum, the DSO's cost is the highest, so Equation (28b) can be expressed as In view of the fact that the above formula belongs to the fuzzy chance constraint and it is difficult to solve directly, one way to solve the fuzzy chance constraint is to convert it into a clear equivalence class and then use the traditional solving process to calculate the clear equivalence model.According to [44], we can obtain the following theorem: Theorem 1. Suppose ξ is degenerated into a one-dimensional fuzzy variable and its membership function is µ.If the function g(x, ξ) has the form g(x, ξ) = h(x) − ξ, then Cr{g(x, ξ) ≤ 0} ≥ α, if and only if h(x) ≤ K α , where x and g are the decision vector and constraint, respectively.Moreover, When ε ≥ 0, α ≥ 1/2, according to the credibility measure function and the above theorem, the robust optimization model shown in Equation ( 28) can be expressed as ) The model considers that the actual spot price fluctuates within a certain range of the predicted spot price.The obtained electricity transaction strategy can ensure that the cost of DSO is less than the expected cost and the credibility is not less than α.The solution process of the proposed method is given by Algorithm 1.
Algorithm 1 Solution process 1: Given system data and forecasted spot price; 2: Considering constraints ( 17)- (27), calculate (16) to obtain the minimum cost C 0 of DSO with predicted spot price; 3: Give DSO risk aversion factor σ or expected cost C e ; 4: Obtain the membership function associated with forecast error percentage of spot price according to the credibility theory and derive its credibility distribution;

Case Analysis
To prove the validity of the proposed model, a modified 15-node distribution network system is selected for numerical study in this paper and the structure of the distribution network system is shown in Figure 2 (the specific parameters are in [45]).Two DGs are set at nodes 1 and 12, respectively.The capacity of each DG is set to 0.15 MW and their power generation cost is USD 30/MWh.The forecast spot price and distribution system load are shown in Figure 3. Suppose that the peak period of electricity consumption is 8 : 00 ∼ 24 : 00 and the rest of the period is non-peak period.The call option strike price λ ck = USD 64.3/MWh and the option premium λ 0 = USD 2.3/MWh.In addition, set the weighting factor ω = 0.33 in the credibility function.In the competitive power market, the cost of DSO depends on its own power generation and power purchase plan.DSO can purchase electricity through different combinations of spot market, option contract and bilateral contract.Assume that DSO has five bilateral contracts to choose from non-peak and peak periods, respectively.The detailed parameters of the bilateral contract are shown in Table 2.In order to verify the effectiveness of the credibility theory-based robust optimization model to hedge price uncertainty, this paper selects different electricity transaction scenarios.Scenario 1, DSO only purchases electricity through the spot market; scenario 2, DSO purchases electricity through the spot market and option contract; scenario 3, DSO purchases electricity through the spot market, bilateral contract and option contract.

Comparison of Transaction Cost under Deterministic Spot Electricity Price
First of all, it is assumed that the actual spot price equals the predicted value.The minimum cost of DSO is obtained by solving the deterministic electricity transaction model.The optimal DSO electricity transactions in different scenarios are shown in Figure 4.In scenario 1, DG output is 7.20 WMh, the spot market purchase is 26.48 MWh and the cost of DSO is USD 1684.0.In scenario 2, DG output is 7.20 WMh, the spot market and option contract purchase are 13.70 MWh and 12.78 MWh, respectively, and the cost of DSO is USD 1682.5.In scenario 3, the option contract purchase is 8.09 MWh, the bilateral contract purchase is 6.42 MWh, the DG output is 7.20 MWh, the spot market purchase is 11.96 MWh and the cost of DSO is USD 1655.6.It can be seen that with the increase of transaction form, the electricity purchase cost of DSO gradually decreases.
For scenario 3, DSO's DG output and electricity transaction in spot market, bilateral contract and option contract are shown in Figure 5.As can be seen from the figure, the nonpeak period electricity transaction market is mainly in the spot market and bilateral contract and peak period electricity transaction market is mainly in the spot market and option contract.Overall, DG output, bilateral contract, spot market and option contract accounted for 21.4%, 19.1%, 35.5% and 24.0% of the total electricity consumption, respectively.The bilateral contract trading volume of each node is shown in Table 3.We can see that bilateral contract transaction is mainly in nodes 2 and 13.Nodes 2 and 13 choose contracts 3, 4, 5, 7, 8, 9, 10 to trade electricity.This is because these two nodes have a high load demand and multiple bilateral contracts can be selected to meet their own demand.The other nodes with low load demand only choose a bilateral contract to trade electricity during peak and non-peak periods.In addition, since the load demand of node 14 is too small, there is no suitable bilateral contract for it to choose, so it meets its own demand through the spot market and option contract.

Robust Optimization Model to Hedge Price Uncertainty of DSO with Multiple Transactions
Assuming that the risk aversion factor σ is 0.1 in scenario 3, the calculated resistible percentage of prediction error is 24.6% and the credibility is 0.83.In this case, if the forecast error percentage of spot price is within the range [0, 24.6%], the cost of DSO is less than or equal to USD 1821.16.If the forecast error percentage of spot price exceeds this range, the actual cost cannot be guaranteed.The electricity transaction strategy of DSO is shown in Figure 6.The credibility associated with the actual cost lower than the expected cost is 0.83, from which the decision maker can assess the risk of the trading strategy.As can be seen in Figure 6, the load demand of DSO at non-peak period is mainly satisfied through DG output and electricity purchase in spot market, while at peak period it is mainly satisfied through option contract and power purchase in spot market.In 24 h, the option contract purchase is 13.77MWh, accounting for 40.90%; the bilateral contract purchase is 9.99 MWh, accounting for 29.67%; the DG output is 6.76 MWh, accounting for 20.08%; and the spot market purchase is 3.15 MWh, accounting for 9.35%.Comparing with Figure 5, it can be seen that considering the uncertainty of spot price, the trading volume of the bilateral contract and the option contract increase, the spot market trading volume decreases and the output of DG decreases slightly.
The bilateral contract trading volume of each node under the uncertainty of spot price is shown in Table 4.We can see that the bilateral contract transaction is still mainly in nodes 2 and 13, but they choose contracts 1-10 to trade electricity and other nodes have also increased the trading volume of bilateral contract.Comparing with the electricity purchase strategy in the deterministic environment, it can be found that in order to reduce the risk caused by the uncertainty of the spot price, DSO increases the trading volume of bilateral contract and option contract.

The Influence of Different Risk Aversion Coefficients of DSO on Electricity Transaction
The curve of electricity transaction strategy with risk aversion factor σ is shown in Figure 7.As the value of σ increases, the purchase volumes of bilateral contract and option contract increase, the purchase volume in the spot market decreases and DG's output remains basically unchanged.The results show that as the expected cost increases, DSO increases the trading volumes of option contract and fixed-price bilateral contract, while reducing volume in spot market with uncertain price.In this way, the robustness of the electricity transaction strategy is increased.
The changes of robustness and credibility with risk aversion factor in different scenarios are shown in Table 5.It can be seen that the credibility increases as the risk aversion factor increases.This shows that the stronger the risk aversion awareness of DSO, the higher the credibility of the expected goal realization.This is because the greater the risk aversion factor, the higher the expected cost.The robustness factor increases with the increasing of expected cost.This shows that the greater the expected cost of the DSO, the lower the acceptance of risk.The more conservative the electricity purchase strategy, the stronger the ability of the resulting electricity transaction strategy to resist risk.In scenario 1, DSO only trades electricity from the spot market.In the event of a bad price that is not conducive to the transaction, there is no electricity purchase plan that can replace or avoid market transactions and it has to accept the market risk caused by price uncertainty.Therefore, the system robustness of scenario 1 is lower than those of other scenarios.
In scenario 2, DSO purchases electricity through bilateral contract and the spot market.The use of fixed-price bilateral contract to purchase electricity avoids to a certain extent the market risk caused by the uncertainty of spot price.
In scenario 3, DSO conducts electricity transaction through spot market, option contract and bilateral contract.It has more means to actively control electricity purchase cost and possible risk losses through reasonable selection of transaction combination and allocation of electricity purchase ratio.In addition, in order to fully demonstrate the effectiveness of the proposed model in dealing with an uncertain problem, this paper compares the proposed robust optimization model with the robust optimization model (RO) [46] and the stochastic optimization model (SO) [47].The optimization results obtained by different optimization methods are shown in Table 6.It can be seen from the table that when the risk aversion factor is 0.05, the cost of the proposed model is reduced by 11.2% compared with robust optimization.This is due to the fact that RO is the worst-case cost of uncertain variables and the resulting electricity trading strategy is too conservative.Compared with random optimization, the solution speed is increased by 80.0%.In addition, compared with RO and SO, the proposed model considers the degree of risk aversion of the decision maker and the decision maker can choose the appropriate risk aversion factor according to their ability to bear the risk.In addition, in order for readers to better understand the proposed model, the credibility theory-based robust optimization model for a user is provided in Appendix A.

Conclusions
Based on the credibility theory, this paper establishes a robust optimization model to hedge price uncertainty of DSO with multiple transactions.This proposed model provides the electricity transaction strategy under different expected cost and the risk-averse DSO achieves the expected goal by rationally allocating the proportion of electricity purchases in different transaction markets.The results of calculation examples show that: (1) Increasing option contract and bilateral contract trading volumes can reduce the electricity transaction cost of DSO by USD 28.5.(2) As the expected cost increases (the degree of risk aversion DSO increases), DSO will increase the purchase of electricity in option contract and bilateral contract, reduce the trading volume in spot market with uncertain price and increase the robustness of electricity transaction strategy.(3) The proposed robust model takes into account the risk aversion of decision maker and obtains the credibility of the expected goal realization.Compared with random optimization, the solution speed is increased by 80.0%.In addition, under the same risk aversion factor, the cost of the proposed model is reduced by USD 195.18 compared with robust optimization and avoids the over-conservatism of traditional robust optimization.This method provides new tools and ideas for electricity transaction decision maker and risk assessment.
This research work only considers the uncertainty of spot electricity price in electricity transaction.In fact, DSO also faces multiple uncertainties brought by renewable energy and demand.In future research work, we will study how to extend the proposed model to measure the multivariate uncertainty and uncertainty coupling.In order to achieve the goal of energy conservation and emission reduction, the impact of green certificates and carbon emissions trading on electricity trading strategy will be studied in the future.

Data Availability Statement:
The data presented in this study are openly available in [45].

Conflicts of Interest:
The authors report no conflict of interest.The authors alone are responsible for the content and writing of the paper.reactive power injected from the power grid, MW a t,l square of the current of line l S t,l upper limit of the apparent power of line l s(l) power outflow end of line l r(l) power inflow end of line l

Figure 1 .
Figure 1.Credibility and membership functions associated with forecast error percentage of spot price.

Figure 3 .
Figure 3. Forecast spot price and distribution system load.

Figure 4 .
Figure 4. DSO electricity transaction strategy in different scenarios.

Figure 5 .
Figure 5. DSO electricity trading volume in each market.

Figure 6 .
Figure 6.Electricity transaction strategy when the risk aversion factor is 0.1.

Figure 7 .
Figure 7.The change curve of electric energy trading volume with σ.

Author
Contributions: L.-P.S.: data curation, formal analysis, writing-original draft preparation, validation.J.-J.C.: resources, conceptualization, visualization, formal analysis writing.L.-W.P.: methodology, visualization.Z.-J.Y.: visualization.All authors have read and agreed to the published version of the manuscript.Funding: This research was supported by Shandong Province Natural Science Foundation of China (No. ZR2022ME047).Institutional Review Board Statement: Not applicable.Informed Consent Statement: Not applicable.

σ
risk aversion factor C e expected cost of DSO, USD R l , X l resistance and reactance of distribution network G b , B b admittance and conductance of distribution network p bc,min n minimum contract volume of bilateral contract, MW p bc,max n maximum contract volume of bilateral contract, MW f active load and reactive load, MW v min t,b , v max t,b maximum and minimum values of node voltage

Table 1 .
Taxonomy of recent research works.

Table 3 .
Bilateral contract trading volume of each node.

Table 4 .
The bilateral contract transaction value of each node under the uncertainty of spot price.

Table 5 .
The changes of robustness coefficient and credibility with risk aversion factors in different scenarios.

Table 6 .
Optimization results under different methods.