1. PROMETHEE for OMCDA

1. PROMETHEE for OMCDA

The PROMETHEE methods belong to the family of outranking methods, which include PROMETHEE I for partial ranking of the alternatives and PROMETHEE II for complete ranking. Other upgrades were also developed in the past decades, such as PROMETHEE III for ranking problems involved interval alternatives, PROMETHEE IV for both partial and complete rankings, and PROMETHEE V for problems with specific constraints.

Consider a problem where a set of possible actions A={a, b, ... } are evaluated on k criteria, i.e. F={f1, f2, ..., fk}. In consideration of the convenience in calculation and without loss of generality, each of the criteria is supposed to be maximized. The problem can be presented in a n×k evaluation matrix, with the actions and criteria being presented in the line and column, respectively. Firstly, a comparison in single criterion is conducted between actions a and b as follows:

Eq. (A.1)

where dk(a,b) is difference between actions a and b under the kth criterion. It is then used in preference function Pj(x). Several types of preference functions were provided in the original PROMETHEE definition, like linear, step or Gaussian preference functions. Among them, the linear criterion preference function is often used in solving MCDA problems, which can be described as follows:

Eq. (A.2)

where qk represents value of an indifference threshold; pk represents value of a strict preference threshold pk. The value of Pk(x) ranges from 0 to 1: when approaching 0, it is considered as negligible by the decision maker and action a is not strictly preferred to action b; when approaching 1, on the contrary, the difference between two actions is significant and action a is highly preferred to action b (as shown in Fig. A.1). When the difference is between two thresholds, an intermediate value is calculated by the linear method as the preference degree.

After estimating the weight for each criterion, the positive and negative flows can be defined, and a multi-criteria preference degree equation is introduced as follows:

Eq. (A.3)

where wk represents weight of criterion fk, assuming the sum value of wk is 1 when

Eq. (A.4)

In order to explore the performance of each alternative ranking against n-1 other ones, the following two outranking flows are defined:

The Positive flow:

Eq. (A.5)

Negative flow:

Eq. (A.6)

The positive flow illustrates how much one decision is outranking over the other alternative. If φ+(a) scores high, then the preference power of action a is distinct. In comparison, the negative flow expresses how much one decision is outranked by the others, and a smaller score of φ-(a) represents strength over other alternatives.

In PROMETHEE I, a partial ranking is defined as follows:

a is preferred to b iif. φ+(a) φ+(b) and φ-(a) φ-(b);

a is indifferent to b iif. φ+(a) =φ+(b) and φ-(a) =φ-(b) ;

a is incomparable to b iif. φ+(a) φ+(b) and φ-(a) φ-(b),or φ+(a) φ+(b) and φ-(a) φ-(b).

PROMETHEE II is evolved on the basis of partial ranking, which combines both positive and negative flows and calculates the net flow by:

Eq. (A.7)

This leads to:

– a is preferred to b iif. φ+(a) φ+(b);

– a is indifferent to b iif. φ+(a) =φ+(b) ;

In addition,

represents net flow for action a under criterion j. All individual net flows of alternatives under each criterion can be summarized in a net flow matrix as follows:

φ1(a) / φ2(a) / … / φj(a) / … / φk(a)
a1 / φ1(a1) / φ2(a1) / … / φj(a1) / … / φk(a1)
a2 / φ1(a2) / φ2(a2) / … / φj(a2) / … / φk(a2)
… / … / … / … / … / …
ai / φ1(ai) / φ2(ai) / … / φj(ai) / φk(ai)
… / … / … / … / … / …
an / φ1(an) / φ2(an) / … / φj(an) / φk(an)

where ai represents action, and each column stands for the net flow of action under a specific criterion. The matrix could provide information related to the performance of alternatives, and is used in the PROMETHEE GAIA process.

2. GAIA for OMCDA

The GAIA method can provide graphical representation for decision makers in a two dimensional view. The conflict among criteria and the impact of the weight on final outranking results can be easily observed in the GAIA plane. Combing with the principal component analysis (PCA) on matrix φ, a plane is define on which only a few information is lost in projection. Before the PCA, the covariance matrix of net flows of all alternatives C was calculated by:

Eq. (A.8)

Where φ’ represents transported matrix of φ. Two k-dimensional eigenvectors are described as u and υ, respectively. In order to minimize the loss by projection, expression u’Cu+v’Cv related to u and υ, has to be maximized by the following description:

Eq. (A.9)

Where cjj is variance ofφj; cjs is covariance matrix between φj and φs; ||γj|| is projection length on the jth criterion; (γj, γs) is scalar product between γj and γs; λ1 and λ2 are the largest and second largest eigenvalues of C, respectively. During the PCA, a multidimensional plane (or say GAIA plane) is built based on λ1 and λ2, in which all the projection of each alternative has the least loss in information.

All projection of criteria γj would be presented as a specific line with individual length and direction. The analysis on projection γj mainly focuses on the following aspects: (1) the length of γj on the GAIA plane represents dominated power of criterion fj; if all alternatives were strongly distinguished under fj, then variance of φj would be very large and the projection of the action would range in an extensive scale, leading to a long projection in the plane, i.e. the longer γj, the more fj distinguishing the alternatives; (2) according to the definition of covariance, criteria expressing the same preference have a large covariance. When two axes of criteria present the same direction in the plane, the criteria are considered as the similar criterion; while when two criteria were independent, covariance cjs between them would close to zero, i.e. axes of two criteria being orthogonal. When two criteria are conflict, cjs would be negative and axes of criteria show an opposite direction in the GAIA plane.

Besides the projection of criterion j on the plane, the PROMETHEE decision axis is also estimated by calculating weight vectors w in the k-dimensional space, which is described as:

Eq. (A.10)

Additionally, the net flows are described in the form of scalar product between vectors ai and w as follows:

Eq. (A.11)

Vector w is decision axis to judge the performance of ai when the weight is taken into account. Since all parameter vectors are discussed in the plane, the projection of weight vectors π is used as the PROMETHEE decision axis. When π is short, the power of the decision axis is weak and the length of π indicates that, some criteria are conflict and the ideal alternative would appear near the origin. When the length of π is long, on the contrary, the alternatives selected in the direction would be desirable. Furthermore, the decision maker could adjust the weight and preference function in the PROMETHEE process, and then the direction and length of decision axes will change correspondingly. Except for the alternatives easy to be distinguished in the GAIA plane, there are some alternatives that are incomparable in PROMETHEE I as their performance is remarkable under different criteria. For these alternatives, PROMETHEE II is introduced for decision makers to make appropriate choice.

Table and Figure Captions List

Fig. A.1 The linear preference function

Table A.1 Preference parameters for scenario 1

Table A.2 Performances of top five pumping-rate groups in scenario 1

Table A.3 Preference parameters for scenario 2

Table A.4 Performances of top five pumping-rate groups in scenario 2

Table A.5 PROMETHEE ranking of five typical pumping-rate groups under 10 years of remediation

Table A.6 Preference parameters for scenario 3

Table A.7 Performances of top five pumping-rate groups in scenario 3

Table A.8 PROMETHEE ranking of five typical pumping-rate groups under 15 years of remediation

Table A.9 Preference parameters for scenario 4

Table A.10 Performances of top five pumping-rate groups in scenario 4

Table A.11 PROMETHEE ranking of five typical pumping-rate groups under 20 years of remediation

Fig. A.1 The linear preference function

Table A.1 Preference parameters for scenario 1

∑Qi / Cc / C / E
Weight / 0.10 / 0.40 / 0.40 / 0.10
Preference function / Linear / Linear / Linear / Linear
p / 100 / 2000 / 2 / 1
q / 50 / 10000 / 15 / 10

Table A.2 Performances of top five pumping-rate groups in scenario 1

Group / ∑Qi / C / Cc / E
Group 32 / 177.0 / 5323.35 / 1.1658 / 13.43
Group 16 / 194.5 / 4668.66 / 1.1235 / 23.04
Group 45 / 281.2 / 2804.14 / 0.9760 / 29.36
Group 24 / 175.9 / 11785.01 / 2.5648 / 27.29
Group 21 / 308.9 / 5515.17 / 2.1084 / 36.66

Table A.3 Preference parameters for scenario 2

∑Qi / Cc / C / E
Weight / 0.15 / 0.33 / 0.45 / 0.07
Preference function / Linear / Linear / Linear / Linear
p / 100 / 2000 / 5 / 1
q / 50 / 10000 / 40 / 20

Table A.4 Performances of top five pumping-rate groups in scenario 2

Group / ∑Qi / C / Cc / E
Group 49 / 147.5 / 352.24 / 0.1037 / 30.56
Group 36 / 181.1 / 3548.39 / 1.2823 / 19.96
Group 28 / 328.5 / 275.81 / 0.1808 / 30.34
Group 32 / 177.0 / 10663.58 / 3.7647 / 13.43
Group 20 / 224.8 / 8593.86 / 3.8538 / 35.45

Table A.5 PROMETHEE ranking of five typical pumping-rate groups under 10 years of remediation

Action / Φ / Φ+ / Φ-
Action49 / 0.8199 / 0.8565 / 0.0366
Action36 / 0.8007 / 0.8179 / 0.0172
Action28 / 0.6586 / 0.7539 / 0.0953
Action32 / 0.6220 / 0.6899 / 0.0679
Action20 / 0.4899 / 0.6049 / 0.1150

Table A.6 Preference parameters for scenario 3

∑Qi / Cc / C / E
Weight / 0.13 / 0.30 / 0.50 / 0.07
Preference function / Linear / Linear / Linear / Linear
p / 100 / 2000 / 10 / 1
q / 50 / 10000 / 40 / 20

Table A.7 Performances of top five pumping-rate groups in scenario 3

Group / ∑Qi / C / Cc / E
Group 26 / 234.3 / 5981.9 / 3.858 / 21.89
Group 13 / 166.8 / 13031.0 / 5.984 / 28.01
Group 29 / 267.7 / 1573.7 / 1.159 / 33.01
Group 39 / 197.2 / 8367.1 / 4.542 / 32.29
Group 49 / 147.5 / 14423.7 / 5.856 / 30.56

Table A.8 PROMETHEE ranking of five typical pumping-rate groups under 15 years of remediation

Action / Φ / Φ+ / Φ-
Action 26 / 0.8135 / 0.8355 / 0.0221
Action 13 / 0.7515 / 0.8096 / 0.0581
Action 29 / 0.7483 / 0.8151 / 0.0668
Action 39 / 0.7449 / 0.8053 / 0.0605
Action 49 / 0.7296 / 0.8055 / 0.0759

Table A.9 Preference parameters for scenario 4

∑Qi / Cc / C / E
Weight / 0.05 / 0.30 / 0.60 / 0.50
Preference function / Linear / Linear / Linear / Linear
p / 100 / 2000 / 15 / 1
q / 50 / 10000 / 50 / 20

Table A.10 Performances of top five pumping-rate groups in scenario 4

Group / ∑Qi / C / Cc / E
Group 47 / 340.3 / 191.0 / 0.228 / 31.99
Group 13 / 166.8 / 8891.6 / 5.207 / 28.01
Group 37 / 211.2 / 10612.0 / 7.866 / 19.53
Group 49 / 147.5 / 12330.1 / 6.384 / 30.56
Group 1 / 243.5 / 11879.4 / 10.153 / 30.14

Table A.11 PROMETHEE ranking of five typical pumping-rate groups under 20 years of remediation

Action / Φ / Φ+ / Φ-
Action 47 / 0.8614 / 0.9150 / 0.0536
Action 13 / 0.8602 / 0.8949 / 0.0347
Action 37 / 0.8482 / 0.8766 / 0.0284
Action 49 / 0.8163 / 0.8635 / 0.0471
Action 1 / 0.7215 / 0.7932 / 0.0717