Introduction

The Intergovernmental Panel on Climate Change’s (IPCC) fifth assessment report unequivocally underscores the intensification of global warming, predominantly caused by human activities. Consequently, addressing the impact of human actions on the environment through greenhouse gas emission reduction has become a paramount concern for humanity. A cost-efficient approach to combat global warming is to establish a carbon market, wherein the pivotal challenge lies in allocating carbon emission quotas. The allocation of carbon quotas typically involves negotiations between the government and diverse enterprises, constituting a consensus problem [1].

So, the allocation of carbon quotas can be seen as the group decision-making [2-4]. Group decision-making is to give full play to the collective wisdom. Numerous participate in the overall process of decision analysis and decision-making. When making decisions on alternatives, because the decision makers often belong to different organizations and have different identities and knowledge backgrounds, the decision makers often have different opinions during the decision-making. At this time, a moderator with rich negotiation and communication experience and personality charm is needed to coordinate all decision makers, so that they can change their views and finally reach an agreement. However, in the coordination process, the moderator always needs to spend some resources to coordinate the decision-makers. In this process, the moderator always wants to spend the least resources and finally reach a consensus [3, 5-9].

Among many research problems of group decision-making, how to minimize the adjustment cost of individual opinions in what we call the consensus reaching process, is a classical problem. Ben-Arieh et al. [10] supposed that the cost of adjusting opinions is linear and studied the issue of minimum cost consensus without budgetary constraints. Ben-Arieh et al. [11] also studied the group consensus costs, and proposed three methods of reaching consensus. In order to study the problem of adjusting opinions between the moderator and the individual, Gong et al. [12] proposed two new types of consensus reaching models, which are involved in minimum cost and maximum return. Ji et al. [13] proposed Risk Maximum Expert Consensus Model (RMECM) based on mean-variance (MV) theory to account for risk factors. The Robust Risk Maximum Expert Consensus Model (R-RMECM) is developed to address uncertainties arising from estimation errors in the mean and covariance matrix of unit adjustment cost. Dong et al. [14] introduced the consensus opinion operator into the decision-making problems first, and the consensus opinion operator adopted the ordered weighted averaging operator and measured the deviation. According to previous studies, Zhang et al. [15] added an aggregation operator into the minimum cost consensus (MCC) model and introduced the idea of soft consensus into the model. Cheng et al. [16] studied the asymmetry in the unit up-adjustment and down-adjustment costs, and three new asymmetry models are proposed: MCCM-DC model, the ε-MCCM-DC model and TB-MCCM-DC model. Based on the above research, Li et al. [17] studied the uncertain scenarios and extended three models, uncertain unit adjustment cost, and initial opinion, which are based on different scenarios discussed. The two-stage stochastic programming is the main method to analyze uncertain parameters.

Existing group-decision-making (GDM) studies always deem that the unit adjustment cost ci, the initial opinion oi in MCC model are known, but the real world is complex and changeable, ci and oi are generally uncertain. So, combining optimization theory with cost consensus problems for group decision-making in uncertain environments is an important branch of research. Stochastic Programming (SP) and Robust Optimization (RO) are two research ideas that exist in existing research approaches to dispose of uncertainties in cost consensus problems [18-21].

In the stochastic programming method, Gong et al. [22] pioneered to combine group decision-making with uncertainty theory to establish a bridge between deterministic group decision and uncertainty theory, and the preferences of decision individuals are fitted using belief degree and uncertainty distribution to construct an opportunity constrained minimum cost consensus model. Ji et al. [23] studied the strategy weight manipulation problem with uncertainty, which applies the uncertainty theory based on confidence degree, and assumed that the uncertain attributes’ values obey uncertainty distribution of linearity, and then constructed a series of hybrid 0-1 planning models, set the strategy weight vector for the required ranking of specific alternatives. The numerical case results of COVID-19 vaccine can reflect the effectiveness of the model well. Liu et al. [24] considered the real-life management problem of “nail households” in urban demolition and eviction, which is a frequent scenario in urban evolution and renewal, and the problem of “nail households” often plagues governments, property developers, and homeowners. As a classic case study, the minimum cost consensus model, multivariate planning, stochastic chance constrained planning and interval number are applied jointly to build a multi-objective minimum cost consensus model based on interval number constraints. But the method relies heavily on the decision maker’s portrayal of the information of the true probability distribution.

In the robust optimization method, Han et al. [25] considered the uncertainty in the input data and established the minimum cost consensus model based on the RO method with four different uncertainty sets. Wei et al. [26] proposed three robust consensus models which used three different aggregation operators, and set novel Box ambiguous set, Ellipsoid ambiguous set, and Polyhedron ambiguous set to study. Jin et al. [27] used the RO method to study the MADM problem and established mixed 0–1 robust optimization model. But in RO models, these worst-case scenarios do not necessarily occur, so consensus cost problems studied by robust optimization-based approaches [1] inevitably lead to over-conservative results.

Although SP and RO are relatively mature methods to study uncertainty problems, the disadvantages of these two methods are also obvious: the results of SP are vulnerable to the accuracy of unknown parameter estimation and the results of RO are too conservative. Therefore, a more reasonable method is needed to deal with uncertainty to ensure that the results are more in line with social reality. In recent years, the Distributionally Robust Optimization (DRO) approach, which has the advantages of SP and RO, is proposed and has been well studied [18, 28-31]. DRO can effectively deal with the uncertainty caused by complex environment, DRO can effectively solve some problems of RO and SP, combining statistical learning and optimization theory to obtain a sufficiently good solution by assuming that the parameters follow certain possible distributions. DRO has found an ambiguity set that contains the true probability distribution of the problem’s random variables. We can fully utilize historical data to construct this ambiguity set and solve the problem. In this study, we adopt the DRO method to avoid the model being too conservative, while also making good use of historical data, ensuring that the model can have good applicability in uncertain environments. But MCCM problems in an uncertain environment will be accompanied by risks, which cannot be completely ignored. The Conditional Value at Risk (CVaR) is a conditional risk measure used to assess financial investment risk. Currently, in the field of group decision-making, more and more researchers are also using CVaR when considering risks and its advantages include:

Tail Risk Consideration: CVaR focuses on the average loss of a portfolio when losses exceed a certain threshold, providing a more comprehensive view of tail risk and better reflecting the extreme events.

Convexity: CVaR is a convex function, making it easier to handle in risk optimization models. Convex optimization problems are generally easier to solve, and CVaR has an advantage in this regard.

Interpretability: The definition of CVaR is relatively simple – it is the average loss when losses exceed a certain threshold. This simplicity makes CVaR easy to interpret and understand for investors and decision-makers.

Due to the above advantages, using the CVaR risk criterion can offset the adverse effects of random variability in the models.

Based on the above findings, this paper uses DRO and CVaR approximation in the TB-MCCM problem, so as to deal with the uncertainty and risk respectively.

The main contributions of this paper are summarized as follows.

(1) Three new distributionally robust two-stage minimum asymmetric adjustment cost consensus models based on CVaR are proposed.

(2) Considering that it is very difficult to obtain historical decision-making data with risk in the carbon market, novel box ambiguous set and polyhedron ambiguous set are constructed, respectively.

(3) The risk expectation cost of the allocation of carbon emission quotas problem under the worst-case condition is measured and three computable linear equivalent forms of the proposed models are derived in order to facilitate calculation.

The rest of this paper is organized as follows.

Section 2 introduces the background and prior knowledge of the minimum cost consensus model, the distributionally robust optimization theory, a coherent risk measure and presents three new distributionally robust two-stage minimum asymmetric adjustment cost consensus models and two ambiguous sets. Section 3 includes numerical examples based on carbon emission quotas and some discussions. Section 5 concludes and presents future work.

Material and Methods

Material

The material used in this paper is some prior knowledge of the minimum cost consensus model, the distributionally robust optimization theory and a coherent risk measure.

In the minimum cost consensus problem, there are n decision makers participating in the discussion. Let oi represent the initial opinion, õ represent the consensus opinion and ci represent the cost adjusting the initial opinion to the consensus opinion. Finally, the model can be described as follows [12]:

(1)
minϕ=∑i=1nci|o~−oi|s.t.o~∈O.

Cheng et al. [16] analyzed the model and found that ci is symmetrical. According to the actual situation, three new minimum cost consensus models with costs which are different in different adjustment directions are proposed. Based on Cheng’s models, the following two-stage stochastic minimum cost consensus models are given by Li et al. [17].

Based on Cheng’s MCCM-DC model, the two-stage stochastic programming method is used and the two-stage stochastic MCCM-DC can be given:

(2)
(TSCM - DC)min𝖤ξ{Q1[o~,ξ(ϖ)]},s.t.omin≤o~≤omax

where ϖ is random event, Q1[o~,ξ(ϖ)] is the optimal value for the second stage recourse problem, õ is the consensus opinion, ξ(ϖ)={c(ϖ),o(ϖ)} and the second stage problem is as follows:

(3)
(Second Stage)minδc(ϖ)Tδs.t.Wδ=o(ϖ)−o~,δ≥0

where c(ϖ)=[c1D(ϖ),c1U(ϖ),…,cnD(ϖ),cnU(ϖ)]T, o(ϖ)=[o1(ϖ),…,on(ϖ)]T, δ=[δ1+,δ1−,…,δn+,δn−]T, o~=[o~,…,o~]T, and

(4)
W=[1−100⋯0000001−1⋯0000⋮⋮⋮⋮⋯⋮⋮⋮⋮0000⋯1−1000000⋯001−1]n×2n

Based on Cheng’s ε-MCCM-DC model, the two-stage stochastic programming method is used and the two-stage stochastic ε-MCCM-DC can be given:

(5)
(ε−TSCM−DC)min𝖤ξ{Q2[o~,ξ(ϖ)]},s.t.omin≤o~≤omax

where ϖ is random event, Q2[o~,ξ(ϖ)] is the optimal value for the second stage recourse problem, õ is the consensus opinion, ξ(ϖ)={c(ϖ),o(ϖ)} and the second stage problem is as follows:

(6)
(Second Stage)minδc(ϖ)Tδs.t.Wδ=o−o~,0≤δ≤ε(ϖ)

where c(ϖ)=[c1D(ϖ),c1U(ϖ),…,cnD(ϖ),cnU(ϖ)]T, ε(ϖ)=[ε1(ϖ),…,εn(ϖ)]T, o=[o1,…,on]T, δ=[δ1+,δ1−,…,δn+,δn−]T, o~=[o~,…,o~]T, and

(7)
W=[1−100⋯0000001−1⋯0000⋮⋮⋮⋮⋯⋮⋮⋮⋮0000⋯1−1000000⋯001−1]n×2n

Based on Cheng’s TB-MCCM-DC model, the two-stage stochastic programming method is used and the two-stage stochastic TB-MCCM-DC can be given:

(8)
(TB - TSCM - DC)min𝖤ξ{Q3[o~,ξ(ϖ)]},s.t.omin≤o~≤omax

where ϖ is random event, Q3[o~,ξ(ϖ)] is the optimal value for the second stage recourse problem, õ is the consensus opinion, ξ(ϖ)={c(ϖ),o(ϖ)} and the second stage problem is as follows:

(9)
(Second Stage)minδc(ϖ)Tδs.t.Wδ=o+θ(ϖ)−o~,δ≥0

where c(ϖ)=[0,c1D(ϖ),c1U(ϖ),0,…,0,ciD(ϖ),ciU(ϖ),0,…,0,cnD(ϖ),cnU(ϖ),0]T, θ(ϖ)=[θ1(ϖ),θ1(ϖ),…,θi(ϖ),θi(ϖ),…,θn(ϖ),θn(ϖ)]T, o=[o1,o1,…,oi,oi,…,on,on]T, o~=[o~,o~,…,o~,o~,…,o~,o~]T∈R2n, δ=[u1+,u1−,v1+,v1−,…,ui+,ui−,vi+,vi−,…,un+,un−,vn+,vn−]T, and

(10)
W=[1−100⋯0000001−1⋯0000⋮⋮⋮⋮⋯⋮⋮⋮⋮0000⋯1−1000000⋯001−1]2n×4n.

Distributionally Robust Optimization Method

The distributed robust optimization model combines the advantages of SP and RO, and the expression is as follows:

(11)
minx∈XsupP∈P𝖤ξ~P[c(ξ)Tx],

where x is the decision variable, X is the decision space of x, c(ξ)Tx is the objective function with random variables, 𝖤ξ~P[·] is the expectation operator, ξ is the column vector of random variables, P is a certain probability distribution of ξ, sup is the abbreviation of supremum and P is the ambiguous set which contains all possible probability distributions.

A Coherent Risk Measure

Suppose μ(x,ξ):Rn×Ω→R is a measurable loss function. The probability of μ(x, ξ) not exceeding a certain threshold ζ can be expressed as Ψ(x,ζ):=∫μ(x,ξ)≤ζp(ξ)dξ. Value at risk (VaR) is defined as follows.

Definition 2.1. [32] Assume that β∈(0,1) is a given confidence level and the function Ψ(x, ξ) is continuous everywhere about ζ. The measure of risk called VaR can be expressed as:

VaRβ(x):=min{ζ∈R:∫μ(x,ξ)≤ζp(ξ)dξ},

where ζ is the threshold value and p(∙) is the probability density function of the random variable.

Artzner et al. [33] proposed the concept of consistent risk measures, i.e., risk measures should satisfy positivity, chi-squaredness, and subadditivity. Since VaR does not satisfy subadditivity, it is not a consistent risk measure. Rockafellar et al. propose a consistent risk measure, CVaR, which is defined as follows.

Definition 2.2. [34] Suppose β∈(0,1) is a given confidence level and the function Ψ(x, ξ) is continuous everywhere about ζ. The mean of the β truncated tail distribution of the loss function μ(x, ξ) and it can be defined as:

CVaRβ(x):=11−β∫μ(x,ξ)≥VaRβ(x)μ(x,ξ)p(ξ)dξ,

where p(·) is the probability density function of the random variable.

Obviously, CVaR indicates that the risk loss is not less than the expected value of VaR. Since the defining equation of CVaR is not easy to calculate under continuous random variables, Rockafellar et al. [34] give an approximate calculation method that CVaR can be equivalently transformed into

CVaRβ(x)=minζ∈RFβ(x,ζ),

where

Fβ(x,ζ):=ζ+11−β∫ξ∈Rn[μ(x,ξ)−ζ]+p(ξ)dξ .

The calculation of Fβ(x, ζ) involves multivariate, non-smooth functions with difficult integration problems. The CVaR approximation under discrete random variables is given by Rockafellar et al. as follows:

Fβ(x,ζ):=ζ+11−β∑i=1npi[μ(x,ξ[i])−ζ]+,

where n indicates the number of scenarios and pi indicates the probability of occurrence of the i scenario.

Lemma 2.1. [35] Assume that β∈(0,1) and pi denote the probability of occurrence of scenario i. The CVaR minimization problem on x∈X is equivalent to the minimization of Fβ(x, ζ) on (x, ζ)∈X×R, i.e.

(12)
minx∈XCVaRβ(x)=min(x,ζ)∈X×Rζ+11−β∑i=1mpi[μ(x,ξ[i])−ζ]+.

Since [∙]+ is a convex function, Fβ(x, ζ) is also convex with respect to (x, ζ)∈X×R.

Methods

Model Construction Considering Risk Aversion with Different Opinion Adjustment Directions

Let oi(ϖ) denote the uncertain initial opinion of DMi. There are two uncertain unit adjustment costs on opinion modification in different directions ciD(ϖ) and ciU(ϖ). Assuming ξ(ϖ)T={ciD(ϖ)T,ciU(ϖ)T,oi(ϖ)T}, λ as a trade-off parameter, is an important parameter used in the formula to balance the expected value and the conditional risk function value (CVaR). The value range of is λ [0,1]. When λ = 0, it indicates that the formula does not consider possible risks at all, while when λ = 1, it indicates that the formula is a complete risk preference type, and the setting of this parameter is influenced by the decision-maker’s risk preference. And the following risk-adjusted cost consensus model (RDRO-DC) with different opinion adjustment directions is given:

(13)
mino~supP∈P(1−λ)𝖤ξ~P[Q1(o~,ξ(ϖ))]+λCVaRα,ξ~P[Q1(o~,ξ(ϖ))]s.t.omin≤o~≤omax,

where Q1(o~,ξ(ϖ)) is the optimal value of the second stage recourse problem, λ∈[0,1] represents the risk aversion measure and the second stage problem is as follows:

(14)
(Second Stage)minδc(ϖ)Tδs.t.Wδ=o(ϖ)−o~ ,δ≥0.
(15)
W=[1−100⋯0000001−1⋯0000⋮⋮⋮⋮⋯⋮⋮⋮⋮0000⋯1−1000000⋯001−1]n×2n,

where c(ϖ)=[c1D(ϖ),c1U(ϖ),…,cnD(ϖ),cnU(ϖ)]T, o(ϖ)=[o1(ϖ),…,on(ϖ)]T, δ=[δ1+,δ1−,…,δn+,δn−]T, o~=[o~,…,o~]T∈Rn, and supP∈P{·} is the worst case of the probability distribution P, i.e., the minimum-maximum robustness criterion.

Model Construction Considering Risk Aversion with Compromise Limits

Similarly, the following worst-case CVaR distribution-based robust model is proposed for the cost consensus modeling problem in conjunction with the finite compromise constraint εi. Let oi(ϖ) denote the uncertain initial opinion of DMi. There are two uncertain unit adjustment costs on opinion modification in different directions ciD(ϖ) and ciU(ϖ). Assuming ξ(ϖ)T={ciD(ϖ)T,ciU(ϖ)T,εi(ϖ)T}, the following risk-adjusted cost consensus model with compromise limit in different opinion adjustment directions (ε-RDRO-DC) is given:

(16)
mino~supP∈P(1−λ)𝖤ξ~P[Q2(o~,ξ(ϖ))]+λCVaRα,ξ~P[Q2(o~,ξ(ϖ))]s.t.omin≤o~≤omax,

where Q2(o~,ξ(ϖ)) is the optimal value of the second stage recourse problem, λ∈[0,1] represents the risk aversion measure and the second stage problem is as follows:

(17)
(Second Stage)minδc(ϖ)Tδs.t.Wδ=o(ϖ)−o~,0≤δ≤ε(ϖ)
(18)
W=[1−100⋯0000001−1⋯0000⋮⋮⋮⋮⋯⋮⋮⋮⋮0000⋯1−1000000⋯001−1]n×2n,

where c(ϖ)=[c1D(ϖ),c1U(ϖ),…,cnD(ϖ),cnU(ϖ)]T, ε(ϖ)=[ε1(ϖ),…,εn(ϖ)]T, o=[o1,…,on]T, δ=[δ1+,δ1−,…,δn+,δn−]T, o~=[o~,…,o~]T∈Rn, and supP∈P{·} is the worst case of the probability distribution P, i.e., the minimum-maximum robustness criterion.

Model Construction Considering Risk-Averse Cost-Free Thresholds

Assuming uncertainty ξ(ϖ)T={ciD(ϖ)T,ciU(ϖ)T,θi(ϖ)T}, the range of opinion adjustment for each DM is denoted by θi(ϖ), which is uncertain and costless. The uncertainty tolerance of DMi is assumed to be θi(ϖ), ui−=[o~−θi(ϖ)−oi]+, ui+=[oi−o~+θi(ϖ)]+, vi−=[o~+θi(ϖ)−oi]+, vi+=[oi−o~−θi(ϖ)]+, where ui−∈[0,o~−θi(ϖ)], ui+∈[o~−θi(ϖ),o~], vi−∈[o~,o~+θi(ϖ)], vi+∈[o~+θi(ϖ),+∞), ui+·ui−=0, vi+·vi−=0 TB-RDRO-DC is expressed in the following form.

(19)
mino~supP∈P(1−λ)𝖤ξ~P[Q3(o~,ξ(ϖ))]+λCVaRα,ξ~P[Q3(o~,ξ(ϖ))]s.t.omin≤o~≤omax,

where Q3(o~,ξ(ϖ)) is the optimal value of the second stage recourse problem, λ∈[0,1] represents the risk aversion measure and the second stage problem is as follows:

(20)
(Second Stage)minδc(ϖ)Tδs.t.Wδ=o+θ(ϖ)−o~,δ≥0
(21)
W=[1−100⋯0000001−1⋯0000⋮⋮⋮⋮⋯⋮⋮⋮⋮0000⋯1−1000000⋯001−1]2n×4n,

where c(ϖ)=[0,c1D(ϖ),c1U(ϖ),0,…,0,ciD(ϖ),ciU(ϖ),0,…,0,cnD(ϖ),cnU(ϖ),0]T, θ(ϖ)=[θ1(ϖ),θ1(ϖ),…,θi(ϖ),θi(ϖ),…,θn(ϖ),θn(ϖ)]T, o=[o1,o1,…,oi,oi,…,on,on]T, o~=[o~,o~,…,o~,o~,…,o~,o~]T∈R2n, δ=[u1+,u1−,v1+,v1−,…,ui+,ui−,vi+,vi−,…,un+,un−,vn+,vn−]T, and supP∈P{·} is the worst case of the probability distribution P, i.e., the minimum-maximum robustness criterion.

Construction of Ambiguous Sets

Two kinds of ambiguous sets, box ambiguous set P1 and polyhedral ambiguous set P2 are constructed [36]:

(22)
P1={p=p0+π|eTπ=0,‖π‖∞≤Ψ},
(23)
P2={p=p0+P1π|eTP1π=0,p0+P1π≥0,‖π‖1≤1},

where p0 is the nominal distribution of discrete probabilities, i.e., the distribution with the highest probability; e is the unit vector; π∈RS is the perturbation vector; Ψ∈[0,1] is the upper limit of fluctuations; and P1∈RS×S is the known perturbation matrix. In polyhedral ambiguous set, the condition eTP1π=0 and the non-negative constraint p0+P1π≥0 ensure that p conforms to the non-negative properties of the probability distribution.

To specify the simplest ambiguous set, it is sensible to consider box-type set and polyhedral ambiguous set, and the resulting problem can be formulated in a computationally tractable way.

Linear Equivalent Reconstruction of the Model

Model Reconstruction Considering Risk Aversion with Different Opinion Adjustment Directions

Through simplification, CVaR is transformed into a computable equation, it is given by:

(24)
maxP∈P𝖤P[Q(x,ξ(ϖ))]=maxP∈PQ(x,ξ(ϖ))TP,

where the probability P has finite support Ω, and P=(P1,P2,⋯,P|Ω|)T, Pϖ>0 is the probability of random event ϖ and ∑ϖ∈ΩPϖ=1.

According to Definition 2.2, CVaR[Q(x,ξ(ϖ))] can be transformed into

(25)
CVaRα,P[Q(x,ξ(ϖ))]=minη∈R+{η+11−αmaxP∈P𝖤P[(Q(x,ξ(ϖ))−η)+]},

where α represents the degree of risk. When α = 0, CVaRα,P[Q(x,ξ(ϖ))] will express as neutral. In contrast, when α → 1, CVaRα,P[Q(x,ξ(ϖ))] will express highly risk-averse, η represents the maximum loss suffered by the function Q(x,ξ(ϖ)) for a given α.

The auxiliary vector t=(t1,t2⋯,t|Ω|)T are introduced into Eq. (25), the Eq. (25) is converted to the following form:

(26)
minη∈R+η+11−αmaxP∈PtTPs.t.Q(x,ξ(ϖ))−eη≤t.t≥0.

Based on the Eq. (25), The objective function of the model (13) can be transformed into

(27)
mino~supP∈Pminη∈R+(1−λ)𝖤ξ~P[Q1(o~,ξ(ϖ))]+λ(η+11−α𝖤ξ~P[(Q1(o~,ξ(ϖ))−η)+]).

The order of the operators supP∈P and minη∈R+ can be changed according to the strongly maximal-minimal nature of Q1(o~,ξ(ϖ)):

(28)
supP∈Pminη∈R+(1−λ)𝖤ξ~P[Q1(o~,ξ(ϖ))]+λ(η+11−α𝖤ξ~P[(Q1(o~,ξ(ϖ))−η)+]).=minη∈R+supP∈P(1−λ)𝖤ξ~P[Q1(o~,ξ(ϖ))]+λ(η+11−α𝖤ξ~P[(Q1(o~,ξ(ϖ))−η)+])=minη∈R+λη+(1−λ)maxP∈P𝖤ξ~P[Q1(o~,ξ(ϖ))]+λ1−αmaxP∈P𝖤ξ~P[(Q1(o~,ξ(ϖ))−η)+].

Therefore, the proposed RDRO-DC model (13) can be equivalently expressed as

(29)
mino~,ηλη+(1−λ)maxP∈P𝖤ξ~P[Q1(o~,ξ(ϖ))]+λ1−αmaxP∈P𝖤ξ~P[(Q1(o~,ξ(ϖ))−η)+],s.t.omin≤o~≤omax.

where maxP∈P𝖤ξ~P[Q1(o~,ξ(ϖ))] and maxP∈P𝖤ξ~P[(Q1(o~,ξ(ϖ))−η)+] depend on the ambiguous set properties of the discrete probability distribution.

Theorem 2.1. Under the box ambiguous set P1, the model (29) can be rewritten as follows:

(30)
mino~,ηλη+(1−λ)[Q1(o~,ξ(ϖ))Tp0+ΨTβ+ΨTγ]+λ1−α[tTp0+ΨTβ′+ΨTγ′]s.t.eμ−β+γ=Q1(o~,ξ(ϖ)),eμ′−β′+γ′=t,Q1(o~,ξ(ϖ))−eη≤t,t≥0,β≥0,γ≥0,β′≥0,γ′≥0,omin≤o~≤omax.Constraints (3.2) - (3.3)

where μ,β,γ,μ′,β′,γ′∈R×R|P1|×R|P1|×R×R|P1|×R|P1| are the auxiliary variables.

Proof: The Eq. (24) can be deduced that

(31)
maxP∈P1𝖤ξ~P[Q1(o~,ξ(ϖ))]=maxP∈P1Q1(o~,ξ(ϖ))TP,

Therefore, the equivalent form of maxP∈P1𝖤ξ~P[(Q1(o~,ξ(ϖ))−η)+] can be expressed as

(32)
maxP∈P1tTPs.t.Q1(o~,ξ(ϖ))−eη≤t.t≥0

Under the box ambiguous set, the model (31) and (32) are linear programming problems and the model (31) is given by the following equation

(33)
maxP∈P1Q1(o~,ξ(ϖ))TP=Q1(o~,ξ(ϖ))Tp0+maxπ{Q1(o~,ξ(ϖ))Tπ|eTπ=0,‖π‖∞≤Ψ},

among them, ‖π‖∞≤max1≤ϖ≤|P1||πϖ|.

In turn, maxπ{Q1(o~,ξ(ϖ))Tπ|eTπ=0,‖π‖∞≤Ψ} can be rewritten as

(34)
maxπQ1(o~,ξ(ϖ))Tπs.t.eTπ=0,−π≤Ψ,π≤Ψ,

among them, Ψ=ψe.

According to the strong duality theory, we can get:

(35)
minμ.β,γΨTβ+ΨTγs.t.eμ−β+γ=Q1(o~,ξ(ϖ)).β≥0,γ≥0

Therefore, the pairwise form of the (31) is expressed as follows:

(36)
minμ.β,γQ1(o~,ξ(ϖ))Tp0+ΨTβ+ΨTγs.t.eμ−β+γ=Q1(o~,ξ(ϖ)).β≥0,γ≥0

In this way, the model (32) of the pairwise planning can be derived in the following form

(37)
minμ′,β′,γ′tTp0+ΨTβ′+ΨTγ′s.t.eμ′−β′+γ′=tQ1(o~,ξ(ϖ))−eη≤t.t≥0β′≥0,γ′≥0

In summary, the model (29) is given by

(38)
mino―,ηλη+(1−λ)[Q1(o~,ξ(ϖ))Tp0+ΨTβ+ΨTγ]+λ1−α[tTp0+ΨTβ′+ΨTγ′]s.t.eμ−β+γ=Q1(o~,ξ(ϖ))eμ′−β′+γ′=tQ1(o~,ξ(ϖ))−eη≤t.t≥0β≥0,γ≥0,β′≥0,γ′≥0omin≤o~≤omaxConstraints (3.2) - (3.3)

With the above proof, it is possible to obtain the model (29) computable reconstruction under the box ambiguous set.

Theorem 2.2. Under the polyhedral ambiguous set P2, the model (29) can be rewritten as follows:

(39)
mino~,η,θ,ζ,v,θ′,ζ′,v′λη+(1−λ)[Q1(o~,ξ(ϖ))Tp0+p0Tθ+v]+λ1−α[tTp0+p0Tθ′+v′]s.t.‖P1TQ1(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v‖P1Tt+P1Tθ′+P1Teζ′‖∞≤v′Q1(o~,ξ(ϖ))−eη≤t,t≥0θ≥0,ζ≥0,θ′≥0,ζ′≥0omin≤o~≤omaxConstraints (3.2) - (3.3)

where θ,ζ,v,θ′,ζ′,v′∈R|P2|×R×R×R|P2|×R×R are the auxiliary variables.

Proof: Under the polyhedral ambiguous set, it can be deduced that

(40)
maxP∈P2Q1(o~,ξ(ϖ))TP=Q1(o~,ξ(ϖ))Tp0+maxπ{Q1(o~,ξ(ϖ))TP1π|eTP1π=0,p0+P1π≥0,‖π‖1≤1},=Q1(o~,ξ(ϖ))Tp0+Γ*(Q1(o~,ξ(ϖ)))

where ‖π‖1≤max1≤ϖ≤|P2||πϖ| and Γ*(Q1(o~,ξ(ϖ))) are the optimal values for the following convex problems:

(41)
maxπ{Q1(o~,ξ(ϖ))TP1π|eTP1π=0,p0+P1π≥0,‖π‖1≤1}.

The Lagrangian function of model (41) is described as follows:

(42)
L(π,θ,ζ,v)=Q1(o~,ξ(ϖ))TP1π+θT(p0+P1π)+ζeTP1π+v(1−‖π‖1).

The Lagrangian pairwise function of model (41) is described as follows:

(43)
g(θ,ζ,v)=maxπL(π,θ,ζ,v)=p0Tθ+v+maxπ{(P1TQ1(o~,ξ(ϖ))+P1Tθ+P1Teζ)π−v‖π‖1},=p0Tθ+v+f*(P1TQ1(o~,ξ(ϖ))+P1Tθ+P1Teζ)

where

(44)
f*(P1TQ1(o~,ξ(ϖ))+P1Tθ+P1Teζ)={0‖P1TQ1(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v∞otherwise.

Therefore, the dual of the (41) is the following problem:

(45)
minθ,ζ,v{p0Tθ+v|‖P1TQ1(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v,θ≥0,ζ≥0}.

On this basis, the model (31) of the equivalence planning can be expressed as follows:

(46)
minμ.β,γQ1(o~,ξ(ϖ))Tp0+p0Tθ+vs.t.‖P1TQ1(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v.θ≥0,ζ≥0,

Similarly, the model (32) of the equivalent form can be rewritten as

(47)
minθ′,ζ′,v′tTp0+p0Tθ′+v′s.t.‖P1Tt+P1Tθ′+P1Teζ′‖∞≤v′Q1(o~,ξ(ϖ))−eη≤t.t≥0θ′≥0,ζ′≥0

In summary, the model (29) is given by

(48)
mino~,η,θ,ζ,v,θ′,ζ′,v′λη+(1−λ)[Q1(o~,ξ(ϖ))Tp0+p0Tθ+v]+λ1−α[tTp0+p0Tθ′+v′]s.t.‖P1TQ1(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v‖P1Tt+P1Tθ′+P1Teζ′‖∞≤v′Q1(o~,ξ(ϖ))−eη≤t.t≥0θ≥0,ζ≥0,θ′≥0,ζ′≥0Constraints (3.2) - (3.3)

With the above proof, it is possible to obtain the model (29) computably reconfigurable under the polyhedral ambiguous set.

Considering Risk Aversion with Compromise Limits for Model Refactoring

Based on the Eq. (25), the objective function of the model (16) can be transformed into

(49)
mino~supP∈Pminη∈R+(1−λ)𝖤ξ~P[Q2(o~,ξ(ϖ))]+λ(η+11−α𝖤ξ~P[(Q2(o~,ξ(ϖ))−η)+]).

The order of the operators supP∈P and minη∈R+ can be changed according to the strongly maximal-minimal nature of Q2(o~,ξ(ϖ)):

(50)
supP∈Pminη∈R+(1−λ)𝖤ξ~P[Q2(o~,ξ(ϖ))]+λ(η+11−α𝖤ξ~P[(Q2(o~,ξ(ϖ))−η)+])=minη∈R+supP∈P(1−λ)𝖤ξ~P[Q2(o~,ξ(ϖ))]+λ(η+11−α𝖤ξ~P[(Q2(o~,ξ(ϖ))−η)+]).=minη∈R+λη+(1−λ)maxP∈P𝖤ξ~P[Q2(o~,ξ(ϖ))]+λ1−αmaxP∈P𝖤ξ~P[(Q2(o~,ξ(ϖ))−η)+]

Therefore, the proposed ε-RDRO-DC model (16) can be equivalently expressed as

(51)
mino~,ηλη+(1−λ)maxP∈P𝖤ξ~P[Q2(o~,ξ(ϖ))]+λ1−αmaxP∈P𝖤ξ~P[(Q2(o~,ξ(ϖ))−η)+],s.t.omin≤o~≤omax

where maxP∈P𝖤ξ~P[Q2(o~,ξ(ϖ))] and maxP∈P𝖤ξ~P[(Q2(o~,ξ(ϖ))−η)+] depend on the ambiguous set properties of the discrete probability distribution.

Theorem 2.3. Under the box ambiguous set P1, the model (51) can be rewritten as follows:

(52)
mino~,ηλη+(1−λ)[Q2(o~,ξ(ϖ))Tp0+ΨTβ+ΨTγ]a+λ1−α[tTp0+ΨTβ′+ΨTγ′]s.t.eμ−β+γ=Q2(o~,ξ(ϖ))eμ′−β′+γ′=tQ2(o~,ξ(ϖ))−eη≤t,t≥0β≥0,γ≥0,β′≥0,γ′≥0omin≤o~≤omaxConstraints (3.5) - (3.6)

where μ,β,γ,μ′,β′,γ′∈R×R|P1|×R|P1|×R×R|P1|×R|P1| are the auxiliary variables.

Proof: The formula (24) can be deduced that

(53)
maxP∈P1𝖤ξ~P[Q2(o~,ξ(ϖ))]=maxP∈P1Q2(o~,ξ(ϖ))TP.

Therefore, the equivalent form of maxP∈P1𝖤ξ~P[(Q2(o~,ξ(ϖ))−η)+] can be expressed as

(54)
maxP∈P1tTPs.t.Q2(o~,ξ(ϖ))−eη≤t.t≥0

Under the box ambiguous set, the model (53) and (54) are linear programming problems and the model (53) is given by the following equation:

(55)
maxP∈P1Q2(o~,ξ(ϖ))TP=Q2(o~,ξ(ϖ))Tp0+maxπ{Q2(o~,ξ(ϖ))Tπ|eTπ=0,‖π‖∞≤Ψ},

among them, ‖π‖∞≤max1≤ϖ≤|P1||πϖ|.

In turn, maxπ{Q2(o~,ξ(ϖ))Tπ|eTπ=0,‖π‖∞≤Ψ} can be rewritten as

(56)
maxπQ2(o~,ξ(ϖ))Tπs.t.eTπ=0−π≤Ψ,π≤Ψ

among them, Ψ=ψe.

According to the strong duality theory, we can get:

(57)
minμ.β,γΨTβ+ΨTγs.t.eμ−β+γ=Q2(o~,ξ(ϖ)).β≥0,γ≥0

Therefore, the pairwise form of the (53) is expressed as follows:

(58)
minμ.β,γQ2(o~,ξ(ϖ))Tp0+ΨTβ+ΨTγs.t.eμ−β+γ=Q2(o~,ξ(ϖ)).β≥0,γ≥0

In this way, the model (54) of the pairwise planning can be derived in the following form:

(59)
minμ′,β′,γ′tTp0+ΨTβ′+ΨTγ′s.t.eμ′−β′+γ′=tQ2(o~,ξ(ϖ))−eη≤t.t≥0β′≥0,γ′≥0

In summary, the model (51) is given by

(60)
mino―,ηλη+(1−λ)[Q2(o~,ξ(ϖ))Tp0+ΨTβ+ΨTγ]+λ1−α[tTp0+ΨTβ′+ΨTγ′]s.t.eμ−β+γ=Q2(o~,ξ(ϖ))eμ′−β′+γ′=tQ2(o~,ξ(ϖ))−eη≤t.t≥0β≥0,γ≥0,β′≥0,γ′≥0omin≤o~≤omaxConstraints (3.5) - (3.6)

With the above proof, it is possible to obtain the model (51) computable reconstruction under the box ambiguous set.

Theorem 2.4. Under the polyhedral ambiguous set P2, the model can be rewritten as follows:

(61)
mino~,η,θ,ζ,v,θ′,ζ′,v′λη+(1−λ)[Q2(o~,ξ(ϖ))Tp0+p0Tθ+v]+λ1−α[tTp0+p0Tθ′+v′]s.t.‖P1TQ2(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v‖P1Tt+P1Tθ′+P1Teζ′‖∞≤v′Q2(o~,ξ(ϖ))−eη≤t,t≥0θ≥0,ζ≥0,θ′≥0,ζ′≥0omin≤o~≤omaxConstraints (3.5) - (3.6)

where θ,ζ,v,θ′,ζ′,v′∈R|P2|×R×R×R|P2|×R×R are the auxiliary variables.

Proof: Under the polyhedral ambiguous set, it can be deduced that

(62)
maxP∈P2Q2(o~,ξ(ϖ))TP=Q2(o~,ξ(ϖ))Tp0+maxπ{Q2(o~,ξ(ϖ))TP1π|eTP1π=0,p0+P1π≥0,‖π‖1≤1},=Q2(o~,ξ(ϖ))Tp0+Γ*(Q2(o~,ξ(ϖ)))

where ‖π‖1≤max1≤ϖ≤|P2||πϖ| and Γ*(Q2(o~,ξ(ϖ))) are the optimal values for the following convex problems:

(63)
maxπ{Q2(o~,ξ(ϖ))TP1π|eTP1π=0,p0+P1π≥0,‖π‖1≤1}.

The Lagrangian function of model (63) is described as follows:

(64)
L(π,θ,ζ,v)=Q2(o~,ξ(ϖ))TP1π+θT(p0+P1π)+ζeTP1π+v(1−‖π‖1).

The Lagrangian pairwise function of model (63) is described as follows:

(65)
g(θ,ζ,v)=maxπL(π,θ,ζ,v)=p0Tθ+v+maxπ{(P1TQ2(o~,ξ(ϖ))+P1Tθ+P1Teζ)π−v‖π‖1},=p0Tθ+v+f*(P1TQ2(o~,ξ(ϖ))+P1Tθ+P1Teζ)

where

(66)
f*(P1TQ2(o~,ξ(ϖ))+P1Tθ+P1Teζ)={0‖P1TQ2(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v∞otherwise.

Therefore, the dual of the (63) is the following problem:

(67)
minθ,ζ,v{p0Tθ+v|‖P1TQ2(o―,ξ(ω))+P1Tθ+P1Teζ‖∞≤v,θ≥0,ζ≥0}.

On this basis, the model (63) of the equivalence planning can be expressed as follows:

(68)
minμ.β,γQ2(o―,ξ(ω))Tp0+p0Tθ+vs.t.‖P1TQ2(o―,ξ(ω))+P1Tθ+P1Teζ‖∞≤v.θ≥0,ζ≥0

Similarly, the model (64) of the equivalent form can be rewritten as

(69)
minθ′,ζ′,v′tTp0+p0Tθ′+v′s.t.‖P1Tt+P1Tθ′+P1Teζ′‖∞≤v′Q2(o―,ξ(ω))−eη≤t.t≥0θ′≥0,ζ′≥0

In summary, the model (51) is equivalent to the following linear programming form

(70)
mino―,η,θ,ζ,v,θ′,ζ′,v′λη+(1−λ)[Q2(o―,ξ(ω))Tp0+p0Tθ+v]+λ1−α[tTp0+p0Tθ′+v′]s.t.‖P1TQ2(o―,ξ(ω))+P1Tθ+P1Teζ‖∞≤v‖P1Tt+P1Tθ′+P1Teζ′‖∞≤v′Q2(o―,ξ(ω))−eη≤t.t≥0θ≥0,ζ≥0,θ′≥0,ζ′≥0Constraints (3.5) - (3.6)

With the above proof, it is possible to obtain the model (51) computably reconfigurable under the polyhedral ambiguous set.

Model Reconstruction Considering Risk Aversion No-Cost Threshold

Based on the Eq. (25), The objective function of the model (19) can be transformed into

(71)
mino~supP∈Pminη∈R+(1−λ)𝖤ξ~P[Q3(o~,ξ(ϖ))]+λ(η+11−α𝖤ξ~P[(Q3(o~,ξ(ϖ))−η)+]).

The order of the operators supP∈P and minη∈R+ can be changed according to the strongly maximal-minimal nature of Q3(o~,ξ(ϖ)):

(72)
supP∈Pminη∈R+(1−λ)𝖤ξ~P[Q3(o~,ξ(ϖ))]+λ(η+11−α𝖤ξ~P[(Q3(o~,ξ(ϖ))−η)+])=minη∈R+supP∈P(1−λ)𝖤ξ~P[Q3(o~,ξ(ϖ))]+λ(η+11−α𝖤ξ~P[(Q3(o~,ξ(ϖ))−η)+]).=minη∈R+λη+(1−λ)maxP∈P𝖤ξ~P[Q3(o~,ξ(ϖ))]+λ1−αmaxP∈P𝖤ξ~P[(Q3(o~,ξ(ϖ))−η)+]

Therefore, the proposed TB-RDRO-DC model (19) can be equivalently expressed as

(73)
mino~,ηλη+(1−λ)maxP∈P𝖤ξ~P[Q3(o~,ξ(ϖ))]+λ1−αmaxP∈P𝖤ξ~P[(Q3(o~,ξ(ϖ))−η)+],s.t.omin≤o~≤omax

where maxP∈P𝖤ξ~P[Q3(o~,ξ(ϖ))] and maxP∈P𝖤ξ~P[(Q3(o~,ξ(ϖ))−η)+] depend on the ambiguous set properties of the discrete probability distribution.

Theorem 2.5. Under the box ambiguous set P1, the model (73) can be rewritten as follows:

(74)
mino~,ηλη+(1−λ)[Q3(o~,ξ(ϖ))Tp0+ΨTβ+ΨTγ]+λ1−α[tTp0+ΨTβ′+ΨTγ′]s.t.eμ−β+γ=Q3(o~,ξ(ϖ))eμ′−β′+γ′=tQ3(o~,ξ(ϖ))−eη≤t,t≥0β≥0,γ≥0,β′≥0,γ′≥0omin≤o~≤omaxConstraints (3.8) - (3.9)

where μ,β,γ,μ′,β′,γ′∈R×R|P1|×R|P1|×R×R|P1|×R|P1| are the auxiliary variables.

Proof: The Eq. (24) can be deduced that

(75)
maxP∈P1𝖤ξ~P[Q3(o~,ξ(ϖ))]=maxP∈P1Q3(o~,ξ(ϖ))TP.

Therefore, the equivalent form of maxP∈P1𝖤ξ~P[(Q3(o~,ξ(ϖ))−η)+] can be expressed as

(76)
maxP∈P1tTPs.t.Q3(o~,ξ(ϖ))−eη≤t.t≥0

Under the box ambiguous set, the model (75) and (76) are linear programming problems. In this case, the model (75) is given by the following equation:

(77)
maxP∈P1Q3(o~,ξ(ϖ))TP=Q3(o~,ξ(ϖ))Tp0+maxπ{Q3(o~,ξ(ϖ))Tπ|eTπ=0,‖π‖∞≤Ψ},

among them, ‖π‖∞≤max1≤ϖ≤|P1||πϖ|.

In turn, maxπ{Q3(o~,ξ(ϖ))Tπ|eTπ=0,‖π‖∞≤Ψ} can be rewritten as

(78)
maxπQ3(o~,ξ(ϖ))Tπs.t.eTπ=0−π≤Ψ,π≤Ψ

among them, Ψ=ψe.

According to the strong duality theory, we can get:

(79)
minμ.β,γΨTβ+ΨTγs.t.eμ−β+γ=Q3(o~,ξ(ϖ)).β≥0,γ≥0

Therefore, the pairwise form of the (75) is expressed as follows:

(80)
minμ.β,γQ3(o―,ξ(ω))Tp0+ΨTβ+ΨTγs.t.eμ−β+γ=Q3(o―,ξ(ω)),β≥0,γ≥0,

In this way, the model (76) of the pairwise planning can be derived in the following form:

(81)
minμ′,β′,γ′tTp0+ΨTβ′+ΨTγ′s.t.eμ′−β′+γ′=tQ3(o~,ξ(ϖ))−eη≤t.t≥0β′≥0,γ′≥0

In summary, the model (73) is equivalent to the following linear programming form

(82)
mino―,ηλη+(1−λ)[Q3(o~,ξ(ϖ))Tp0+ΨTβ+ΨTγ]+λ1−α[tTp0+ΨTβ′+ΨTγ′]s.t.eμ−β+γ=Q3(o~,ξ(ϖ))eμ′−β′+γ′=tQ3(o~,ξ(ϖ))−eη≤t.t≥0β≥0,γ≥0,β′≥0,γ′≥0omin≤o~≤omaxConstraints (3.8) - (3.9)

With the above proof, it is possible to obtain the model (73) computable reconstruction under the box ambiguous set.

Theorem 2.6. Under the polyhedral ambiguous set P2, the model (73) can be rewritten as follows:

(83)
mino~,η,θ,ζ,v,θ′,ζ′,v′λη+(1−λ)[Q3(o~,ξ(ϖ))Tp0+p0Tθ+v]+λ1−α[tTp0+p0Tθ′+v′]s.t.‖P1TQ3(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v‖P1Tt+P1Tθ′+P1Teζ′‖∞≤v′Q3(o~,ξ(ϖ))−eη≤t,t≥0θ≥0,ζ≥0,θ′≥0,ζ′≥0omin≤o~≤omaxConstraints (3.8) - (3.9)

where θ,ζ,v,θ′,ζ′,v′∈R|P2|×R×R×R|P2|×R×R are the auxiliary variables.

Proof: Under the polyhedral ambiguous set, it can be deduced that

(84)
maxP∈P2Q3(o~,ξ(ϖ))TP=Q3(o~,ξ(ϖ))Tp0+maxπ{Q3(o~,ξ(ϖ))TP1π|eTP1π=0,p0+P1π≥0,‖π‖1≤1},=Q3(o~,ξ(ϖ))Tp0+Γ*(Q3(o~,ξ(ϖ)))

where ‖π‖1≤max1≤ϖ≤|P2||πϖ| and Γ*(Q3(o~,ξ(ϖ))) are the optimal values for the following convex problems:

(85)
maxπ{Q3(o~,ξ(ϖ))TP1π|eTP1π=0,p0+P1π≥0,‖π‖1≤1}.

The Lagrangian function of model (85) is as follows:

(86)
L(π,θ,ζ,v)=Q3(o~,ξ(ϖ))TP1π+θT(p0+P1π)+ζeTP1π+v(1−‖π‖1).

The Lagrangian pairwise function of model (85) is as follows:

(87)
g(θ,ζ,v)=maxπL(π,θ,ζ,v)=p0Tθ+v+maxπ{(P1TQ3(o~,ξ(ϖ))+P1Tθ+P1Teζ)π−v‖π‖1},=p0Tθ+v+f*(P1TQ3(o~,ξ(ϖ))+P1Tθ+P1Teζ)

where

(88)
f*(P1TQ3(o~,ξ(ϖ))+P1Tθ+P1Teζ)={0‖P1TQ3(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v∞otherwise.

Therefore, the dual of the (85) is the following problem:

(89)
minθ,ζ,v{p0Tθ+v|‖P1TQ3(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v,θ≥0,ζ≥0}.

On this basis, the model (75) of the equivalence planning can be expressed as follows:

(90)
minμ.β,γQ3(o~,ξ(ϖ))Tp0+p0Tθ+vs.t.‖P1TQ3(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v.θ≥0,ζ≥0

Similarly, the model (76) of the equivalent form can be rewritten as

(91)
minθ′,ζ′,v′tTp0+p0Tθ′+v′s.t.‖P1Tt+P1Tθ′+P1Teζ′‖∞≤v′Q3(o~,ξ(ϖ))−eη≤t.t≥0θ′≥0,ζ′≥0

In summary, the model (73) is given by

(92)
mino~,η,θ,ζ,v,θ′,ζ′,v′λη+(1−λ)[Q3(o~,ξ(ϖ))Tp0+p0Tθ+v]+λ1−α[tTp0+p0Tθ′+v′]s.t.‖P1TQ3(o~,ξ(ϖ))+P1Tθ+P1Teζ‖∞≤v‖P1Tt+P1Tθ′+P1Teζ′‖∞≤v′Q3(o~,ξ(ϖ))−eη≤t.t≥0θ≥0,ζ≥0,θ′≥0,ζ′≥0Constraints (3.8) - (3.9)

With the above proof, it is possible to obtain the model (73) computably reconfigurable under the polyhedral ambiguous set.

Results and Discussion

Due to the large amount of carbon dioxide generated by burning fossil fuels such as oil and coal, these greenhouse gasses can strongly absorb infrared radiation from the ground, leading to an increase in Earth’s temperature. Global warming not only endangers the balance of natural ecosystems, but also affects human health and even threatens human survival.

Carbon emissions trading is a market mechanism adopted to promote global greenhouse gas emissions reduction and reduce global carbon dioxide emissions. The United Nations Intergovernmental Panel on Climate Change passed the United Nations Framework Convention on Climate Change on May 9, 1992 through difficult negotiations. The first additional agreement to the Convention was passed in Kyoto, Japan in December 1997. The Protocol regards market mechanisms as a new path to solve the problem of greenhouse gas emissions reduction represented by carbon dioxide, which regards carbon dioxide emission rights as a commodity, thus forming the trading of carbon dioxide emission rights, abbreviated as carbon trading. A complete carbon market has been established in China, and the key to it is how to allocate carbon emission quotas. Carbon quotas refer to the legal amount obtained by each enterprise through bargaining with the government over a certain period of time.

In this context, the negotiation of carbon quotas between enterprises and the government has actually formed a consensus process, and the numerical simulation cases in this article are all based on the above issues. Assuming that the government is negotiating carbon quotas with 20 companies, and the adjustment costs for each company in different directions are uncertain, we assume that each company has its own initial opinions (ideas) on carbon quotas. All numerical calculations were done with a laptop computer (Intel i7-7300HQ CPU, 32G RAM), using CPLEX 12.10 to run the solution under JDK11.

Calculation Results

The following example explains the proposed distributionally robust cost consensus model considering risk factors. Suppose P0=(0.3,0.6,0.1)T, Ψ=(0.02,0.04,0.06,0.08,0.1), P1=Γ*I=Ψ|P2|*I, where I is unit vector and |P2|=3, α=0.5, and λ=(0.1,0.3,0.5,0.7,0.9).

Example 1. (Based on RDRO-DC)

The proposed model’s initial information data is given by Table 1.

The RDRO-DC model is solved using the CPLEX solver. According to the experimental results in Table 2, the RDRO-DC model’s consensus opinion is 6.13 and the minimum consensus cost is 381.12.

Example 2. (Based on ε-RDRO-DC)

The proposed model’s initial information data is given by Table 3.

The ε-RDRO-DC model is solved using the CPLEX solver. According to the experimental results in Table 4., the ε-RDRO-DC model’s consensus opinion is 5.87 and the minimum consensus cost is 375.6.

Example 3. (Based on TB-RDRO-DC)

The proposed model’s initial information data is given by Table 5.

The TB-RDRO-DC model is solved using the CPLEX solver. According to the experimental results in Table 6., the TB-RDRO-DC model’s consensus opinion is 6.06 and the minimum consensus cost is 378.5.

Table 1.

RDRO-DC initial information data.

ioiciDciU
14.21523
26.71026
36.21122
45.31224
56.11325
65.51131
73.81621
84.91223
94.61524
105.01625
115.71231
126.31126
134.91432
145.31323
152.82139
163.21830
174.51424
186.81125
197.1918
203.61930
Table 2.

RDRO-DC Consensus Cost Results.

ΨλBox Ambiguous SetPolyhedral Ambiguous SetΨλBox Ambiguous SetPolyhedral Ambiguous Set
0.020.1349.6381.20.040.1367.5416.2
0.3362.8403.60.3384.1433.9
0.5374.8427.40.5398.7457.6
0.7391.2459.00.7418.3480.2
0.9417.3492.70.9437.7501.3
0.060.1377.8431.70.080.1403.8462.4
0.3401.1454.00.3427.6483.1
0.5423.6472.80.5459.3506.3
0.7447.2498.30.7481.2534.2
0.9472.8524.70.9506.9557.9
0.10.1415.2486.2
0.3437.9504.9
0.5461.7526.7
0.7482.5551.8
0.9513.6578.3

Comparison and Analysis of Results

Different combinations of Ψ and λ are used in the numerical case. As can be seen in Table 2., Table 4., and Table 6., the optimal target value of the model increases as the parameter Ψ increases; because when the parameter Ψ is reduced, the accurate probabilistic information will be obtained and the model will be more accurate and more convenient to solve. It can also be observed that when the coefficient λ decreases, the optimal target value of the model decreases, as CVaR becomes a smaller proportion of the objective function, and there is no need to consider more risks to increase costs.

As can be seen from Table 7., the optimal target values of the models increase when the level of risk increases. This is because when the value of α tends to 1, CVaR will express highly risk-averse, and this high level of risk will lead to the overly conservative solution, which means the optimal target value increases. By observing the three uncertainty models proposed in this paper, it is observed that the optimal target value under the box ambiguous set is always smaller than the optimal target value under the polyhedral ambiguous set. This is because the probability distribution of the polyhedral ambiguous set is more perturbed than that of the box ambiguous set. In other words, the box ambiguous set is simpler than the polyhedral ambiguous set, so its probability distribution information is more comprehensive than that of the polyhedral ambiguous set.

Table 3.

ε-RDRO-DC initial information data.

ioiciDciUεi
14.2152324
26.7102622
36.2112218
45.3122421
56.1132515
65.5113116
73.8162125
84.9122318
94.6152420
105.0162522
115.7123114
126.3112618
134.9143216
145.3132324
152.8213916
163.2183020
174.5142422
186.8112518
197.191820
203.6193024
Table 4.

ε-RDRO-DC Consensus Cost Results.

ΨλBox Ambiguous SetPolyhedral Ambiguous SetΨλBox Ambiguous SetPolyhedral Ambiguous Set
0.020.1357.2375.60.040.1373.4388.7
0.3371.8393.40.3389.2401.3
0.5384.6417.50.5403.1429.2
0.7402.5438.60.7422.7442
0.9419.8457.20.9446.8458.6
0.060.1384.6398.60.080.1401.5409.7
0.3407402.10.3426.8433.3
0.5423.5433.90.5448.2453.8
0.7451.6444.80.7469.7467.6
0.9478.2468.20.9493.1488.4
0.10.1418.4412.8
0.3439.6443.5
0.5463.9461.6
0.7484.5474.1
0.9517.3495.9

When Ψ = 0, the distributionally robust models will degenerate to a traditional stochastic programming model. The difference between the distributionally robust models and the traditional stochastic programming model will be discussed in the following paragraph. As can be seen from Table 8., under the same risk level, the optimal target values of the uncertainty models are greater than that of the certainty model, this shows that the uncertainty models are more conservative, this conservatism is caused by the robustness of the model. In the traditional stochastic programming model, the uncertainty of the data was not considered, although the uncertainty model increases the cost, the model can avoid the failure caused by the data uncertainty to a certain extent. Therefore, the distributionally robust models are important to the consensus reaching in MCCM.

Compared with the study of Ji et al. [37], According to the results shown in Table 9., it can be found that the distributionally robust models have a significant advantage over the stochastic programming model. Whether based on the box ambiguous set or higher cost polyhedral ambiguous sets, the total consensus cost to be paid is better than the stochastic programming model in most cases. Therefore, when faced with a consensus decision problem that considers the decision maker’s risk preferences, the model is able to obtain consensus opinions at a lower consensus cost.

Figs 1-3. show the comparison of consensus costs of the three models in this paper, the model in Ji et al. [37] and the no-risk model under different α. It can be seen that the addition of risk function will obviously increase the consensus cost, which means that the model is risk sensitive. The total consensus cost of distributionally robust model with risk aversion is higher than that without risk aversion. In other words, it is more difficult to reach consensus when the risk is considered. Meanwhile, the two-stage stochastic minimum cost consensus model with risk aversion is more conservative than that without risk aversion, so the results are worse. However, a comparison with the results under existing stochastic programming-based methods with risk factors shows that the consensus cost of the model proposed in this paper would be better than the existing studies.

Table 5.

TB-RDRO-DC initial information data.

ioiciDciUθi
14.215233
26.710266
36.211225
45.312242
56.113254
65.511317
73.816213
84.912238
94.615246
105.016254
115.712315
126.311263
134.914326
145.313234
152.821398
163.218302
174.514245
186.811256
197.19183
203.619302
Table 6.

TB-RDRO-DC Consensus Cost Results.

ΨλBox Ambiguous SetPolyhedral Ambiguous SetΨλBox Ambiguous SetPolyhedral Ambiguous Set
0.020.1358.9378.50.040.1371.6401.9
0.3367.2391.60.3382.8438.3
0.5381.4413.40.5397.6462.8
0.7393.5436.80.7416.7473.5
0.9411.6453.20.9432.5496.1
0.060.1382.4413.80.080.1397.6418.5
0.3406.1442.70.3416.7432.7
0.5421.3468.90.5438.5458.2
0.7443.5491.30.7467.3481.6
0.9467.3503.40.9491.7504.3
0.10.1427.8422.6
0.3451.1439.8
0.5480.2467.2
0.7493.5483.9
0.9509.3507.5
Table 7.

Optimal results of the model with different α.

αUncertainty modelBox Ambiguous SetPolyhedral Ambiguous SetDeterministic modelOptimal results
0.1RDRO-DC217.8236.7TSMCCM-DC184.6
ε-RDRO-DC218.4232.5ε-TSMCCM-DC189.3
TB-RDRO-DC223.6235.6TB-TSMCCM-DC181.5
0.3RDRO-DC249.8268.3TSMCCM-DC184.6
ε-RDRO-DC242.7268.9ε-TSMCCM-DC189.3
TB-RDRO-DC245.1264.3TB-TSMCCM-DC181.5
0.5RDRO-DC282.6309.4TSMCCM-DC184.6
ε-RDRO-DC284.5312.8ε-TSMCCM-DC189.3
TB-RDRO-DC288.4306.7TB-TSMCCM-DC181.5
0.7RDRO-DC327.1344.8TSMCCM-DC184.6
ε-RDRO-DC325.8341.9ε-TSMCCM-DC189.3
TB-RDRO-DC319.3340.2TB-TSMCCM-DC181.5
0.9RDRO-DC349.6381.2TSMCCM-DC184.6
ε-RDRO-DC357.2375.6ε-TSMCCM-DC189.3
TB-RDRO-DC358.9378.5TB-TSMCCM-DC181.5
Table 8.

Comparison of the results of different models.

αUncertainty modelBox Ambiguous SetPolyhedral Ambiguous SetΨ = 0
0.1RDRO-DC217.8236.7196.2
ε-RDRO-DC218.4232.5197.3
TB-RDRO-DC223.6235.6203.5
0.3RDRO-DC249.8268.3217.8
ε-RDRO-DC242.7268.9219.4
TB-RDRO-DC245.1264.3221.3
0.5RDRO-DC282.6309.4257.9
ε-RDRO-DC284.5312.8253.8
TB-RDRO-DC288.4306.7258.6
0.7RDRO-DC327.1344.8278.1
ε-RDRO-DC325.8341.9279.5
TB-RDRO-DC319.3340.2283.4
0.9RDRO-DC349.6381.2301.2
ε-RDRO-DC357.2375.6306.5
TB-RDRO-DC358.9378.5307.5
Table 9.

Comparison of the results of the Distributionally Robust and Stochastic Programming models.

αDistributionally Robust ModelBox Ambiguous SetPolyhedral Ambiguous SetStochastic Planning ModelOptimal results
0.1RDRO-DC217.8236.7MRMCCM-DC242.7
ε-RDRO-DC218.4232.5ε-MRMCCM-DC246.5
TB-RDRO-DC223.6235.6TB-MRMCCM-DC245.9
0.3RDRO-DC249.8268.3MRMCCM-DC281.9
ε-RDRO-DC242.7268.9ε-MRMCCM-DC273.4
TB-RDRO-DC245.1264.3TB-MRMCCM-DC266.5
0.5RDRO-DC282.6309.4MRMCCM-DC296.1
ε-RDRO-DC284.5312.8ε-MRMCCM-DC316.2
TB-RDRO-DC288.4306.7εB-MRMCCM-DC311.9
0.7RDRO-DC327.1344.8MRMCCM-DC362.8
ε-RDRO-DC325.8341.9ε-MRMCCM-DC371.3
TB-RDRO-DC319.3340.2TB-MRMCCM-DC373.5
0.9RDRO-DC349.6381.2MRMCCM-DC400.2
ε-RDRO-DC357.2375.6ε-MRMCCM-DC398.7
TB-RDRO-DC358.9378.5TB-MRMCCM-DC392.1
Fig. 1.

Comparison of the results of RDRO-DC, stochastic programming method and risk-free model.

https://www.pjoes.com/f/fulltexts/182899/PolJEnvironStud-2024-33-5-5065-5085_Fig_1_min.jpg
Fig. 2.

Comparison of the results of ε-RDRO-DC, stochastic programming approach and risk-free model.

https://www.pjoes.com/f/fulltexts/182899/PolJEnvironStud-2024-33-5-5065-5085_Fig_2_min.jpg
Fig. 3.

Comparison of the results of TB-RDRO-DC, stochastic programming method and risk-free model.

https://www.pjoes.com/f/fulltexts/182899/PolJEnvironStud-2024-33-5-5065-5085_Fig_3_min.jpg

Conclusions and Future Research

Against the backdrop of global warming and increasingly mature carbon emissions trading, the problem of distributional robust consensus decision-making with decision maker risk preferences is considered. The three distributionally robust two-stage minimum asymmetric cost consensus models with risk aversion based on carbon emission quotas are investigated, using a novel box ambiguous set and polyhedral ambiguous set. Through the transformation, the relevant formulas that are convenient for calculation are obtained. Finally, a number of meaningful conclusions can be obtained through numerical examples.

(1) Different values of Ψ, λ and α will lead to the changes in models’ optimal target values The parameter Ψ indicates the accurate probabilistic information, the optimal target value increases as the parameter Ψ increases. The parameter λ represents the proportion of the objective function, the λ decreases, the coefficient λ decreases, the optimal target value decreases. The parameter α shows the level of risk, the optimal target value increases when the value of α tends to 1, which means the model is highly risk-averse.

(2) The optimal target value under the box ambiguous set is always smaller than the optimal target value under the polyhedral ambiguous set. The main reason is: probability distribution information of the box ambiguous set is more comprehensive than that of the polyhedral ambiguous set.

(3) By comparing the distribution robust model with the stochastic programming model, the distribution robust model can use the distribution information in a more robust way than the stochastic programming model. Whether based on the box ambiguous set or polyhedral ambiguous set, the total consensus cost of the distribution robust model to be paid is better than the stochastic programming method in most cases.

(4) The consensus model is sensitive to the risk factor. The total consensus cost of distributionally robust model with risk aversion is higher than that without risk aversion, and it is more difficult to reach consensus when the risk is considered.

In short, it will make the model better avoid uncertainty and risk that introducing distributionally robust method and CVaR into consensus problem. The consensus cost of the models in this paper outperforms the corresponding results under stochastic programming methods and brings new solutions to decision problems in the allocation of carbon emission quotas.

One limitation inherent in this study pertains to its exclusive focus on small-scale group consensus models, involving fewer than 20 decision-makers. In reality, however, large-scale group decision-making scenarios involve the participation of more than 20 individuals, and these have not been comprehensively explored. As a prospective avenue for further research, we intend to delve into the realm of large-scale group decision-making problems in subsequent studies. This will involve leveraging social network methodologies and other innovative approaches to extend the applicability of the findings from this study to larger decision-making groups.

Additionally, another limitation stems from the exclusive consideration of box and polyhedral ambiguity sets. In future research endeavors, we aim to broaden our exploration by incorporating multiple types of ambiguity sets, such as ellipsoid or a combination of box and polyhedral ambiguity sets. This approach will enable us to investigate and understand better the impact of different ambiguity set structures on managing uncertainty.