Next Article in Journal
CoreBug: Improving Effort-Aware Bug Prediction in Software Systems Using Generalized k-Core Decomposition in Class Dependency Networks
Next Article in Special Issue
Co-Training Method Based on Semi-Decoupling Features for MOOC Learner Behavior Prediction
Previous Article in Journal
Existence and Uniqueness Theorem for a Solution to a Class of a Third-Order Nonlinear Differential Equation in the Domain of Analyticity
Previous Article in Special Issue
Selection of Cold Chain Logistics Service Providers Based on a Grey AHP and Grey COPRAS Framework: A Case Study in Vietnam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A MULTIMOORA-Based Risk Evaluation Approach for CCUS Projects by Utilizing D Numbers Theory

1
Heilongjiang Province Key Laboratory of Geographical Environment Monitoring and Spatial Information Service in Cold Regions, Harbin Normal University, Harbin 150025, China
2
School of Management, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(5), 204; https://doi.org/10.3390/axioms11050204
Submission received: 6 April 2022 / Revised: 21 April 2022 / Accepted: 22 April 2022 / Published: 26 April 2022

Abstract

:
As the global climate warms, carbon emissions must be reduced in order to alleviate the human climate crisis. Carbon capture, utilization and storage (CCUS) is an emerging technology that can reduce carbon emissions. However, most of the CCUS projects have ended in failure. The reason can be attributed to insufficient risk assessment. To this end, the purpose of this study is to construct a comprehensive risk assessment model for CCUS projects. The main body of this research is divided into two parts. First, in order to evaluate the CCUS project, a risk indicator system is constructed. In what follows, a decision-making framework for risk assessment under the D numbers environment is proposed, including two stages of decision-making preparation and decision-making process. The main task of the preparation stage is to gather evaluation experts and collect decision-making information. In the decision-making stage, this paper takes the D numbers theory as the core (acting on the effective expression and fusion of subjective evaluation information), respectively, proposes the method of determining the weight of risk evaluators, the fusion method of decision-making information from different experts, and the comprehensive decision model based on the MULTIMOORA method. In order to verify the effectiveness of the constructed model, the case of CCUS project site selection in Shengli power plant is analyzed, and the results showed that the third site is the best option. This study finds the importance of a comprehensive and timely risk assessment for the successful implementation of CCUS projects, and suggests that stakeholders carry out a risk assessment of CCUS projects prior to implementation based on the method presented in this paper, so as to improve the success rate.

1. Introduction

As the global climate warms, the concentration of carbon must be controlled, and there are two effective ways to do this. The first is to reduce carbon emissions, and the second is to deploy negative emissions technologies (NETs) [1] to remove carbon from the atmosphere and sequester it reliably. Carbon capture and storage (CCS) is the process of separating CO 2 from industrial or related energy sources and transporting it to a storage site, where it is isolated from the atmosphere for a long time. CCS technologies mainly fall into four categories [2,3,4,5,6,7]: pre-combustion capture, oxygen-containing fuel combustion [8], chemical chain combustion [9], and post-combustion capture. CCS technology is of great significance to achieve global climate goals [10].
In recent years, the global climate change situation has become increasingly serious [11]. According to the International Energy Agency (IEA), the higher the temperature rise, the more CO 2 reductions CCS will have to contribute [12]. The Intergovernmental Panel on Climate Change (IPCC) stressed the importance of achieving net zero CO 2 emissions by the middle of this century, and proposed that the large-scale application of CCS technology is the key to achieving zero carbon emissions, so CCS technology is of great significance to achieve global climate goals [10].
CCS is an important approach to deep CO 2 emission reduction, but the potential of capture and storage, implementation difficulty, and socio-economic benefits of different methods vary greatly [13]. With the development of CCS technology and the deepening of understanding, China first proposed Carbon Capture, Utilization and Storage (CCUS) technology at the Xiangshan Conference in Beijing in 2006, and introduced the technology of CO 2 resource utilization [14]. CCUS technology, after purifying captured CO 2 , puts it into a new production process for recycling and reuse, and recycles CO 2 . It can not only achieve carbon emission reduction, but also produce economic benefits, so it is more practical and operational. After years of development, CCUS technology has been accepted and used worldwide [15].
CCUS is an important carbon emission reduction technology, and more and more CCUS projects will be invested in because of its effectiveness in addressing climate change [16]. With the operation of CCUS project, more and more environmental risks as well as economic and social problems are being discovered [17], which can be attributed to the risk management lagging behind the development of CCUS technology [18,19]. The main obstacles facing CCUS projects in the commercialization stage are weak awareness of environmental risks, lack of legal constraints, and small financial support [20]. In addition to defects in basic research, supporting technology, and management policies, other factors that hinder the large-scale deployment of CCUS projects include immaturity of technology, law, economy, environment, and society [21]. The failure experience of many CCUS projects shows that environmental impact and risk are one of the key factors to determine whether CCUS projects can be successfully implemented, which is also the focus of public concern [22]. There is no practical process for assessing the risk of CCUS projects, from either the perspective of government or industry, particularly in critical areas such as the selection of carbon storage sites.
Risk assessment is a complete process composed of risk identification, risk analysis, and risk assessment. How the process is carried out depends not only on the background of the risk management process, but also on the methods and techniques used to execute the risk assessment work [23]. For CCUS projects, various risk assessment tools were proposed. A series of risk management processes were proposed by ISO 31000: 2018 [24], including determining assessment content → risk identification → analysis → evaluation → treatment → risk monitoring → review → communication → negotiation. OXAND conducts risk assessments for carbon storage projects based on the ISO 31000 risk management framework [25]. IFPEN, SINTEF, and TNO jointly developed the integrated carbon risk assessment tool [26]. IACRAS defined risk assessment as four modules: scenario definition, scenario analysis, impact assessment, and uncertainty assessment [27]. The CCS environmental risk assessment proposed by the UK Environment Agency uses a typical source → path → receptor model, and evaluates the possibility and consequences of exposure at four levels: high, medium, low, and very low [28]. Recently, Liu et al. [29] evaluated the risk of CCUS projects based on decision modeling, constructed risk indicators from the perspective of sustainable development, and proposed the risk assessment methods under intuitionistic fuzzy and linguistic environments.
Through the analysis of the existing literature on CCUS project risk assessment, it is found that there is not a complete evaluation index system, and especially few studies on measuring the implementation risk of CCUS project from the dimension of sustainable development. In addition, evaluation information generally comes from human subjective judgment, while existing evaluation methods do not provide flexible and effective information expression. For this problem, the theory of D numbers is a good choice [30,31], which has the ability to express uncertain information and provides broader conditions than the Dempster–Shafer theory. However, the fusion rules in the traditional D numbers theory have been questioned because they can not meet the characteristics of associative law, so this study proposes a new method to solve this problem. Furthermore, the current algorithm is more single-dimension assessment, which weakens the reliability of the assessment. Therefore, this paper will provide a comprehensive risk assessment method based on the MULTIMOORA method to provide integrated judgment for CCUS project risks.
This paper contributes to the risk evaluation for CCUS projects through three aspects of innovation. First, it introduces the D numbers to express the subjective judgment of the evaluators, and further proposes a new generation of fusion rule for D numbers to satisfy the associative law. Secondly, a comprehensive risk indicator system for evaluating CCUS projects is developed to guide the conduct of risk assessment. Thirdly, a MULTIMOORA-based decision-making model is proposed to characterize the degree of risk of the CCUS projects. The constructed model addresses the limitations of existing risk assessment approaches by optimizing evaluation indicators and improving evaluation algorithms.
The rest of this paper is organized as follows. Some preliminaries, such as the D numbers theory and MULTIMOORA method, are introduced in Section 2. Section 3 gives a complete description of the research methods, including the risk assessment index system and the comprehensive decision-making model. Section 4 shows the usage and effectiveness of the proposed decision-making method through an example of CCUS project risk assessment. The conclusion of this paper is given in Section 5, followed by some future research plans.

2. Preliminaries

In this section, the basic theory of D numbers and the MULTIMOORA method are introduced, respectively.

2.1. D Numbers

The D numbers theory is the expansion of the Dempster–Shafer theory [32,33]. It has more advantages to represent and handle uncertain information because the elements do not need to be mutually exclusive in D numbers and the completeness constraint is released in D numbers. As a generalization of evidence theory, the D numbers theory has a wide application, especially in linguistic assessment [34,35,36,37]. It is also been applied in other fields, such as healthcare waste management [38,39], risk assessment [40], and fuzzy decision-making [41]. The basic knowledge of D numbers will be introduced as follows.
Definition 1. 
Let Ω be a finite nonempty set, D numbers is a mapping formulated by:
D : Ω [ 0 , 1 ]
with:
B Ω D ( B ) 1 a n d D ( ) = 0
whereis an empty set and B is a subset of Ω.
An illustrative example is given below to show the D numbers.
Example 1. 
Suppose a project is assessed, and the assessment score is represented by an interval [ 0 , 100 ] . If an expert gives their assessment result by D numbers, it could be:
D ( { b 1 } ) = 0.2 D ( { b 3 } ) = 0.6 D ( { b 1 , b 2 , b 3 } ) = 0.1
where b 1 = [ 0 , 45 ] , b 2 = [ 38 , 73 ] , and b 3 = [ 61 , 100 ] . Note that the set of { b 1 , b 2 , b 3 } is not a frame of discernment actually, because the elements in the set of { b 1 , b 2 , b 3 } are not mutually exclusive. Due to D ( { b 1 } ) + D ( { b 3 } ) + D ( { b 1 , b 2 , b 3 } ) = 0.9 , the information is incomplete.
For a discrete set Ω = { b 1 , b 2 , , b i , , b n } , where b i R and b i b j if i j , a special form of D numbers can be expressed by:
D ( { b 1 } ) = v 1 D ( { b 2 } ) = v 2 D ( { b i } ) = v i D ( { b n } ) = v n
or simply denoted as D = { ( b 1 , v 1 ) , ( b 2 , v 2 ) , , ( b i , v i ) , , ( b n , v n ) } , where v i > 0 and i = 1 n v i 1 .
Some properties of D numbers are introduced as follows.
Definition 2. 
Permutation invariability. If there are two D numbers that:
D 1 = { ( b 1 , v 1 ) , , ( b i , v i ) , , ( b n , v n ) }
and:
D 2 = { ( b n , v n ) , , ( b i , v i ) , , ( b 1 , v 1 ) } ,
then D 1 D 2 .
Example 2. 
If there are two D numbers:
D 1 = { ( 0 , 0.7 ) , ( 1 , 0.3 ) } a n d D 2 = { ( 1 , 0.3 ) , ( 0 , 0.7 ) }
Then:
D 1 D 2
Definition 3. 
For D = { ( b 1 , v 1 ) , ( b 2 , v 2 ) , , ( b i , v i ) , , ( b n , v n ) } , the integration representation of D is defined as:
I ( D ) = i = 1 n b i v i
where b i R , v i > 0 and i = 1 n v i 1 .
Example 3. 
Let D = { ( 1 , 0.2 ) , ( 2 , 0.1 ) , ( 3 , 0.3 ) , ( 4 , 0.3 ) , ( 5 , 0.1 ) } , then:
I ( D ) = 1 × 0.2 + 2 × 0.1 + 3 × 0.3 + 4 × 0.3 + 5 × 0.1 = 3.0
Next, a combination rule is proposed to combine two D numbers as below.
Definition 4. 
Let D 1 and D 2 be two D numbers, that:
D 1 = { ( b 1 1 , v 1 1 ) , , ( b i 1 , v i 1 ) , , ( b n 1 , v n 1 ) }
D 2 = { ( b 1 2 , v 1 2 ) , , ( b j 2 , v j 2 ) , , ( b m 2 , v m 2 ) }
The combination of D 1 and D 2 , indicated by D = D 1 D 2 , is defined as:
D ( b ) = v
with:
b = b i 1 + b j 2 2
v = v i 1 + v j 2 2 v i 1 + v j 2 2 C C
C = j = 1 m i = 1 n ( v i 1 + v j 2 2 ) , i = 1 n v i 1 = 1 a n d j = 1 m v j 2 = 1 ; j = 1 m i = 1 n ( v i 1 + v j 2 2 ) + j = 1 m ( v c 1 + v j 2 2 ) , i = 1 n v i 1 < 1 a n d j = 1 m v j 2 = 1 ; j = 1 m i = 1 n ( v i 1 + v j 2 2 ) + i = 1 n ( v i 1 + v c 2 2 ) , i = 1 n v i 1 = 1 a n d j = 1 m v j 2 < 1 ; j = 1 m i = 1 n ( v i 1 + v j 2 2 ) + j = 1 m ( v c 1 + v j 2 2 ) + i = 1 n ( v i 1 + v c 2 2 ) + v c 1 + v c 2 2 , i = 1 n v i 1 < 1 a n d j = 1 m v j 2 < 1 .
where v c 1 = 1 i = 1 n v i 1 and v c 2 = 1 j = 1 m v j 2 . Note that superscript in the above equations is not exponent.
Example 4. 
If two D numbers:
D 1 = { ( 0 , 0.7 ) , ( 1 , 0.3 ) } a n d D 2 = { ( 0 , 0.6 ) , ( 1 , 0.4 ) }
the combination of D 1 and D 2 using Equations (5)–(8) is:
D = { ( 0.0 , 0.325 ) , ( 0.5 , 0.500 ) , ( 1.0 , 0.175 ) }

2.2. MULTIMOORA Method

In 2006, Brauersand and Zavadskas [42] first proposed the MOORA (Multi-Objective Optimization on the basis of a Ratio Analysis) method, which includes two submethods, the ratio system method (RSM) and the reference point method (RPM). Subsequently, they extended MOORA and proposed the MULTIMOORA (Multi-Objective Optimizaion on the basis of a Ratio Analysis plus the full MULTIplicative form) method [43], which added a new submethod called the full-multiplicative form (FMF) on the basis of MOORA. The MULTIMOORA method is widely used in decision making [44,45,46] and evaluation [47,48,49]. Assuming X is the initial decision matrix, expressed as X = ( x i j ) m × n , where x i j is the evaluation value of alternative A i under attribute C j , i = {1,2,…,m}, j = {1,2,…,n}. In order to facilitate comparison, the initial decision matrix needs to be standardized, and the standardized decision matrix X * = ( x i j * ) m × n is obtained, and x i j * is the standardized form of x i j [43]:
x i j * = x i j / i = 1 m x i j 2

2.2.1. The Ratio System Method (RSM)

The evaluation value of alternative A i under the ratio system method is:
y i = j = 1 g x i j * j = g + 1 n x i j *
where g and n g represent the number of benefit-type and cost-type attributes, respectively. The larger the value of y i , the better the corresponding alternative. Therefore, the optimal scheme under the ratio system method is:
A R S M * = { A i | m a x i y i ) }

2.2.2. The Reference Point Method (RPM)

First, the optimal reference point r j for each attribute is determined as:
r j = m a x j ( x i j * ) , j g m i n j ( x i j * ) , j > g
The deviation of the attribute value x i j * from the corresponding reference point r j can be represented as r j x i j * . Thus, the maximum bias of each alternative, i.e., the evaluation value of each alternative under the reference point method, can be expressed as:
z i = m a x j | r j x i j * |
The smaller the value of z i , the better the corresponding alternative. Therefore, the optimal alternative under the reference point method is:
A R P M * = { A i | m i n i z i ) }

2.2.3. The Full-Multiplicative Form (FMF)

The evaluation value of each alternative under full-multiplicative form can be expressed as:
u i = j = 1 g x i j * j = g + 1 n x i j *
where j = 1 g x i j * is the product of standardized evaluation value of benefit-type attributes in alternative A i and j = g + 1 n x i j * is the product of cost-type attributes in alternative A i . The larger the value of the u i , the better the corresponding alternative. Therefore, the optimal alternative under full-multiplicative form is:
A F M F * = { A i | m a x i u i ) }

3. Research Methodology

Figure 1 shows the structure of the proposed risk evaluation model. The proposed model consists of two parts, namely the decision-making preparation stage and the decision-making process stage of risk assessment. The former includes forming an expert team, constructing an evaluation index system, collecting evaluation information (represented by D numbers), and transforming it into several decision matrices. The latter mainly covers the determination of the weight of risk evaluators, the aggregation of multiple decision matrices (using the defined fusion rule), the determination of attribute weights, and the proposal of a decision model based on the weighted MULTIMOORA method. The specific methods are discussed as follows.

3.1. A Comprehensive Risk Indicator System for Evaluating CCUS Projects

In this paper, an indicator system is constructed to evaluate the risk of CCUS projects from five sustainability dimensions: economic, social, environmental, governance, and technology [50]. The sources of risk evaluation indicators mainly refer to the indicator system of Liu et al. [29], other references, and the opinions of relevant experts. In this paper, a two-layer evaluation indicator system is used. The specific content of the evaluation indicates are shown in Table 1, and the structure, explanation, and source are also given.

3.2. An Integrated Decision-Making Model

3.2.1. Problem Description of Risk Evaluation for CCUS Projects

The development alternatives of CCUS projects are represented as set O = { o 1 , o 2 , , o m } . In order to accurately measure the risk of each element and select the best alternative, experts in the field (represented as E = { e 1 , e 2 , , e q } , corresponding weight as W E = { ω e 1 , ω e 2 , , ω e q } , ω e k 0 , ω e k = 1 ) are invited to evaluate the alternatives against each risk indicator (denoted as C = { C 1 , C 2 , , C n } , corresponding weight as W C = { ω C 1 , ω C 2 , , ω C n } , ω C j 0 , ω C j = 1 ). The subjective judgment of expert e k on alternative A i against index c j is expressed as a D number ( D i j k ) m × n , in which the evaluation level is defined in Table 2. That is, the frame of discernment is Ω = { C + , B + , A + , N , A , B , C } .
The ultimate goal of this paper is to complete the selection of optimal alternatives for CCUS project by calculating risks; to this end, the decision-making information from different experts needs to be aggregated first, then the MULTIMOORA method is used to sort all alternatives, and finally the optimal alternative can be selected.

3.2.2. The New Combination Rule for D Numbers

In Definition 4, the combination rule for two D numbers is given, but it must be pointed out that the combination operation defined in Definition 4 does not preserve the associative property. It is clear that ( D 1 D 2 ) D 3 D 1 ( D 2 D 3 ) ( D 1 D 3 ) D 2 . To determine the order of combination, μ j is used as the order variable for each D j in [54], but how to obtain μ j has not been mentioned. To ensure that multiple D numbers can be combined correctly and efficiently, in this paper, a combination operation for multiple D numbers is developed as follows.
Definition 5. 
Let D 1 , D 2 , ⋯ , D C are c D numbers, which are shown as follows:
D 1 = { ( b 1 1 , v 1 1 ) , , ( b i 1 , v i 1 ) , , ( b n 1 , v n 1 ) }
D 2 = { ( b 1 2 , v 1 2 ) , , ( b j 2 , v j 2 ) , , ( b m 2 , v m 2 ) }
D C = { ( b 1 C , v 1 C ) , , ( b k C , v k C ) , , ( b p C , v p C ) }
The first step is to obtain a D numbers by averaging all the subset B of the finite nonempty set Ω:
D A V G = { ( b 1 A V G , v 1 A V G ) , , ( b t A V G , v t A V G ) , , ( b q A V G , v q A V G ) }
where b t A V G = b i 1 = b j 2 = = b k C i f t h e y a r e t h e s a m e , v t A V G = v i 1 + v j 2 + + v k C C .
The second step is to combine the averages of c D numbers c 1 times using Definition 4L
D 1 D 2 D C = D A V G D A V G , , D A V G c
An example illustrates in detail the use of the proposed new fusion rule of D numbers, and more detailed steps are shown in Appendix A. It is clear that the new combination rule has no limit to the number of D numbers and it satisfies the associative property.

3.2.3. Define the Weights of the Risk Evaluators

In order to determine the weight of each evaluator, this study measures the entropy of all evaluation information made by each evaluator. Since the evaluation information is expressed by D numbers, the entropy function of D number proposed in the literature [55] is used in this paper, so the entropy of evaluator e k is:
E ( e k ) = 1 m n j = 1 n i = 1 m E ( D i j k ) ,
where E ( D i j k ) represents the entropy of D number D i j k , which is defined as follows.
Let D = { ( b 1 , v 1 ) , ( b 2 , v 2 ) , , ( b i , v i ) , , ( b n , v n ) } be a D number that contains n mutually exclusive and exhaustive hypotheses. When the information of D is incomplete, an operation is needed to complete it. Otherwise, skip this step.
v ˜ i = v i + 1 n ( 1 j = 1 n v j )
where v i is the belief of proposition b i . Thus, the complete D number can be denoted as D ˜ = { ( b 1 , v ˜ 1 ) , ( b 2 , v ˜ 2 ) , , ( b i , v ˜ i ) , , ( b n , v ˜ n ) } , where v ˜ i > 0 and i = 1 n v ˜ i = 1 . Then, the entropy function of D ˜ can be defined as
E ( D ˜ ) = i = 1 n v ˜ i · I ( v ˜ i )
where I ( v ˜ i ) indicates the information content defined as:
I ( v ˜ i ) = l o g ( 1 v ˜ i )
According to the theory of belief entropy [55], entropy represents the magnitude of uncertainty. The greater the entropy provided by the evaluator, the greater the uncertainty, i.e., the less information, and the smaller the weight, and vice versa. Thus, the weight of evaluator can be defined as:
ω e k = 1 E ( e k ) q k = 1 q E ( e k )

3.2.4. Calculate the Aggregated Risk Evaluation Matrix

In the above two subsections, the aggregation method of D numbers and the determination method of the evaluator’s weight have been proposed, respectively. In what follows, the evaluation information of evaluators on different alternatives against each attribute needs to be aggregated, so a weighted average aggregation method should be proposed in this paper.
Definition 6. 
Let the evaluation information for alternative A i against attribute C j by all the evaluators be represented as a set of D numbers { D i j 1 , D i j 2 , , D i j k } on the frame of discernment Ω = { E + , D + , C + , B + , A + , N , A , B , C , D , E } . The weight vector of evaluators is W E = { ω e 1 , ω e 2 , , ω e q } . The weighted average value of each proposition in D number is calculated as follows:
v ˜ i j = k = 1 q ω e k v i j k
so the weighted average D number is D ˜ i j = ( b i j , v ˜ i j ) . The fusion rule of D numbers is used to perform q 1 self-fusion for D ˜ i j , and the D number after aggregation of each evaluator can be obtained:
D i j = D ˜ i j D ˜ i j , , D ˜ i j q 1
In order to lay a foundation for the use of the MULTIMOORA method [42,43], it is necessary to transform the decision matrix based on D number expression into numerical form, so the element I i j in the decision matrix representing the risk evaluation of alternative A i against attribute C j is calculated as:
I i j = θ Ω b i j θ v i j θ
where b i j θ indicates the numerical form of the grade (the second column in Table 2) for proposition θ and v i j θ represents its reliability. The decision matrix in numerical form has been obtained.

3.2.5. The MULTIMOORA-Based Risk Evaluation Approach

In order to evaluate CCUS projects, multiple risk indicators can be involved, and different indicators have different degrees of importance, so the weight of indicators needs to be determined. This subsection considers from two perspectives, namely, the uncertainty of evaluation information and the differentiation degree of indicators for alternatives.
Firstly, the weight of risk indicators is defined from the perspective of uncertainty (refer to the method in Section 3.2.3), and the entropy of indicator C j is defined as:
E ( C j ) = 1 m i = 1 m E ( D i j ) ,
the weight of indicator C j can be determined as:
ω ( C j ) = 1 E ( C j ) n j = 1 n E ( C j )
Then, the weight of indicators can be defined from the perspective of discrimination degree based on entropy weight method.
(1) The decision matrix should be normalized as:
p i j = I i j i = 1 m I i j
where I i j is the numerical form of the decision matrix obtained by Equation (26).
(2) The entropy value can be calculated as:
e j = 1 l n ( m ) j = 1 n p i j l n ( p i j )
(3) The degree of divergence of the risk indicator C j can be calculated as:
d i v j = 1 e j
where d i v j represents the inherent contrast strength of risk indicator C j , and the higher d i v j is, the more important C j is. Thus, the weight of indicator C j can be determined as:
ω ( C j ) = d i v j j = 1 n d i v j
By combining the weights of the above two aspects, the final weight of indicator C j can be obtained as:
ω ( C j ) = ω ( C j ) ω ( C j ) j = 1 n ω ( C j ) ω ( C j )
Furthermore, combined with the MULTIMOORA method and the obtained weights, a weighted version of the MULTIMOORA method can be defined for risk evaluation. The decision matrix aggregated by multiple evaluators can be expressed as I = ( I i j ) m × n , where I i j is the numerical evaluation value of alternative A i against attribute C j , i = {1,2,…,m}, j = {1,2,…,n}. The standardized decision matrix I * = ( I i j * ) m × n of I can be obtained as:
I i j * = I i j / i = 1 m I i j 2
The evaluation value of alternative A i under the ratio system method can be defined as:
Y i = j = 1 g ω ( C j ) I i j * j = g + 1 n ω ( C j ) I i j *
where g and n g represent the number of benefit-type and cost-type attributes, respectively. The larger the value of y i , the better the corresponding alternative. Therefore, the optimal scheme under the ratio system method is A R S M * = { A i | m a x i Y i ) } .
The optimal reference point R j for attribute C j is determined as:
R j = max j ( I i j * ) , j g m i n j ( I i j * ) , j > g
The evaluation value of each alternative under the reference point method can be expressed as:
Z i = m a x j { ω j | R j I i j * | }
The smaller the value of z i , the better the corresponding alternative. Therefore, the optimal alternative under the reference point method is A R P M * = { A i | m i n i Z i ) } .
The evaluation value of each alternative under full-multiplicative form can be expressed as:
U i = j = 1 g ( I i j * ) ω j j = g + 1 n ( I i j * ) ω j
where j = 1 g x i j * is the product of standardized evaluation value of benefit-type attributes in alternative A i and j = g + 1 n x i j * is the product of cost-type attributes in alternative A i . The larger the value of the U i , the better the corresponding alternative. Therefore, the optimal alternative under full-multiplicative form is A F M F * = { A i | m a x i U i ) } .
Finally, the decision results obtained by RSM, RPM, and FMF models are combined for comprehensive decision making. In this paper, the geometric mean method [56] is used for calculation:
S i = Y i × 1 Z i × U i 3
where S i represents the score index of each alternative, and the larger it is, the higher the alternative ranking will be.

4. Case Study

4.1. Case Description

This section selects the Shengli power plant as a case to illustrate the effectiveness of the proposed decision-making model in CCUS projects risk assessment. The site selection is of great significance for the smooth implementation of the CCUS project [57]. Due to the difficulty and high requirements of carbon storage, strict risk assessment is required for the selection of storage site. Based on the literature [29], considering six dimensions (i.e., geographical and climatic, geological conditions, drainage conditions, distance from cities, ecological protection standards, and land-use costs), four alternatives are selected for the CCUS project of the Shengli power plant. In what follows, the risk assessment model proposed in this paper will be utilized to determine the best one from the four alternative sites based on the constructed risk indicator system.

4.2. Decision-Making Process

The background of the risk assessment is introduced above, and the decision-making process is executed below.
Decision preparation: three experts in this field are invited to form an expert group, denoted as E = { e 1 , e 2 , e 3 } . Experts judge the performance of the four alternative sites under various risk indicators based on their own experience, and give the evaluation information in the form of D numbers, as shown in Table 3, i.e., three decision matrices are obtained.
Decision process: firstly, determine the weight of the experts. According to the method proposed in Section 3.2.3, the weights of three experts are ω e 1 = 0.0347 , ω e 2 = 0.5473 , and ω e 3 = 0.4180 , respectively; secondly, the novel fusion rule for D numbers proposed in Section 3.2.2 is employed to aggregate the risk assessment results of multiple experts, and the results shown in Table 4 are obtained; thirdly, the MULTIMOORA-based risk evaluation approach is used to compare and sort multiple alternative sites. The last three columns of Table 5 give the weight calculation results of the indicators from the perspectives of uncertainty and differentiation degree. In addition, Table 5 also records the standardized matrix of D number. The weighted version of the MULTIMOORA method is enabled, and the results of R S M , P R M , and F M F are obtained, respectively, as shown in Table 6. For the final decision, the results of these three methods are aggregated through Equation (39) to obtain the scores shown in Table 6.

4.3. Results and Analysis

The risk assessment results for the CCUS project site selection problem are shown in Table 6. The table indicates that the three sub-methods of the MULTIMOORA-based D numbers method (i.e., RSM, PRM, AND FMF) come to a consistent conclusion: Option 3 Option 4 Option 1 Option 2. In addition, the fusion result of the three methods is the same conclusion. In order to achieve further visualization of the risk assessment results, two of the results of the three methods (i.e., RSM, PRM, AND FMF), are, respectively, taken to establish two-dimensional coordinates, and the results shown in Figure 2, Figure 3 and Figure 4 are obtained. In order to clearly compare different positions, the background of the coordinate system is set to a color that changes uniformly from bottom left to top right; the redder the color, the higher the risk, while the bluer the color, the lower the risk. From a geometric point of view, the closer to the origin, the higher the risk. In other words, from the origin to the upper right, the risk decreases. According to the above principles, the four alternative sites can be easily compared in terms of risk level from Figure 2 and Figure 3. However, in Figure 4, Options 1, 3, and 4 are relatively close, indicating that the discrimination of the four alternatives is not obvious in the dimensions of PRM and FMF, but it still does not affect the acquisition of the final conclusion.
Based on the above analysis, the practical significance of this paper lies in: (1) the information expression and processing method of CCUS project risk assessment based on the D numbers theory can effectively represent the subjective judgment of experts and improve the flexibility of the assessment process; (2) the constructed assessment index system can effectively guide the risk assessment of CCUS projects and promote the systematization and standardization of the assessment; (3) the risk assessment method based on MULTIMOORA can analyze the results of risk assessment from different dimensions, which is conducive to the visualization of risks and novel discoveries.

4.4. Sensitivity Analysis

In this section, a sensitivity analysis is performed to verify the impact of changes in certain parameters on CCUS project risk assessment. The weight of the indicator is crucial to the result of the evaluation method. The weight of the indicator in Equation (33) comes from two aspects, namely uncertainty and differentiation. A new weight synthesis method is defined as: ω ( C j ) = ω ( C j ) β + ω ( C j ) ( 1 β ) . Parameter β is used to adjust the proportion of the two aspects. In this experiment, β is set to change from 0 to 1 at an interval of 0.1 to observe the final evaluation result, as shown in Figure 5.
Figure 5 shows that in R S M and F M F , the ranking results of the four options are very stable, while in P R M , there are some fluctuations. However, the best option is always in a stable state, and the final results of the MULTIMOORA method show that the ranking of the four options is stable, indicating that the evaluation method in this paper is robust under different indication weights.

4.5. Comparison and Analysis

In order to further demonstrate the effectiveness of the assessment method in this study, in this subsection, our method is compared and analyzed from both qualitative and quantitative aspects. The selected comparison methods are six multi-attribute decision-making methods based on the D numbers theory.
In terms of qualitative analysis, comparisons are made from the perspectives of the main contributions, applications, and aggregation method, and the analysis of the results are shown in Table 7.
Based on the results of comparative experiments, only our study essentially improves the fusion rules of D numbers, make its basic theory more perfect, and only our study has made a prominent contribution in terms of application, i.e., it constructs the indicator system of CCUS project risk evaluation. In addition, our study also obtains decision results from more perspectives, which is more robust. The above aspects can illustrate the effectiveness of the method in this paper.
In addition, we also made a quantitative comparison of the above methods, and the results are shown in Table 7. All methods identified Option 3 as the best, followed by Option 4. There are only differences in the ranking of the last two options, but this does not affect the selection of the final optimal option. The results of the quantitative analysis illustrate the availability and effectiveness of the proposed method.

5. Conclusions

Based on the D numbers theory and the MULTIMOORA method, this study puts forward a comprehensive decision-making model for the risk assessment of CCUS projects. This paper applies the theory of uncertain multi-attribute decision making (MADM) to the risk assessment of CCUS projects, which has the following novelties:
(1) Using the existing literature and relevant domain knowledge, this paper constructs a risk indicator system of CCUS project risk assessment, which lays the attribute foundation for the risk assessment based on the MADM method;
(2) The D numbers theory is used to express the evaluator’s information, and a novel fusion rule is proposed, which not only solves the problem that the traditional fusion method does not meet the combination law, but also realizes the effective information expression and fusion of CCUS project risk assessment;
(3) From the perspective of information entropy, the weights of evaluators and risk factors are determined to capture valuable information and make the decision results more accurate;
(4) The traditional MULTIMOORA method is extended by using the D numbers theory, which not only effectively expresses the uncertainty of evaluators’ assessment, but also makes the evaluation results more robust due to the use of three different functions for decision making.
In terms of practical significance, the CCUS project risk assessment model developed in this study examined a large number of indicators in the fields of economy, society, environment, governance, and technology, providing sufficient reference information and decision-making basis for stakeholders. Furthermore, the D numbers used in this study and the improved fusion rules can help decision-makers flexibly express their subjective judgments. In addition, the evaluation model based on the MULTIMOORA method can analyze risks from different dimensions, which is conducive to risk visualization and novel discovery.
This study also has some limitations. On the one hand, the constructed risk assessment index system is still relatively incomplete. On the other hand, using the constructed data to demonstrate and verify the method, its effectiveness needs to be further improved. This leads to the scope of future research: (1) more diversified expression methods of evaluation information are considered in order to more accurately express the subjective judgment of evaluators, (2) the indicator system of CCUS project risk assessment need to be further improved, (3) more extensive decision-making methods other than the MULTIMOORA method, such as ORESTE and PROMETHEE approaches, are employed to evaluate the risk of CCUS projects, and (4) experts could be invited to participate in the risk assessment of CCUS projects.

Author Contributions

Conceptualization, P.L. and L.W.; methodology, P.L.; software, L.C.; validation, L.F., P.L. and L.W.; formal analysis, L.W.; investigation, L.W.; resources, L.F.; data curation, L.C.; writing—original draft preparation, P.L.; writing—review and editing, L.W.; visualization, L.F.; supervision, L.W.; project administration, L.C.; funding acquisition, L.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 42071079, 41671100).

Acknowledgments

The authors greatly appreciate the reviewers’ suggestions and the editor’s encouragement.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Let D 1 , D 2 , D 3 be three D numbers that:
D 1 = { ( 0 , 0.6 ) , ( 1 , 0.4 ) } , D 2 = { ( 0 , 0.5 ) , ( 1 , 0.5 ) } , D 3 = { ( 0 , 0.1 ) , ( 1 , 0.9 ) }
The average D numbers of D 1 , D 2 , and D 3 can be obtained by Equation (17):
D A V G = { ( 0 , 0.6 + 0.5 + 0.1 3 ) , ( 1 , 0.4 + 0.5 + 0.9 3 ) } = { ( 0 , 0.4 ) , ( 1 , 0.6 ) }
Then, the combination of D 1 , D 2 and D 3 can be calculated by Equation (18) as follows:
D 1 D 2 D 3 = D A V G D A V G D A V G = { ( 0 , 0.4 ) , ( 1 , 0.6 ) } { ( 0 , 0.4 ) , ( 1 , 0.6 ) } { ( 0 , 0.4 ) , ( 1 , 0.6 ) }
First, we calculate { ( 0 , 0.4 ) , ( 1 , 0.6 ) } { ( 0 , 0.4 ) , ( 1 , 0.6 ) } . Based on Equation (8):
C = 0.4 + 0.4 2 + 0.4 + 0.6 2 + 0.6 + 0.4 2 + 0.6 + 0.6 2 = 2
Based on Equations (6) and (7), b 11 = ( 0 + 0 ) / 2 = 0 , b 12 = ( 0 + 1 ) / 2 = 0.5 , b 21 = ( 1 + 0 ) / 2 = 0.5 , b 22 = ( 1 + 1 ) / 2 = 1 . v 11 = ( ( 0.4 + 0.4 ) / 2 ) / 2 = 0.2 , v 12 = ( ( 0.4 + 0.6 ) / 2 ) / 2 = 0.25 , v 21 = ( ( 0.6 + 0.4 ) / 2 ) / 2 = 0.25 , v 22 = ( ( 0.6 + 0.6 ) / 2 ) / 2 = 0.3 . So { ( 0 , 0.4 ) , ( 1 , 0.6 ) } { ( 0 , 0.4 ) , ( 1 , 0.6 ) } = { ( 0 , 0.2 ) , ( 0.5 , 0.5 ) , ( 1 , 0.3 ) } .
So, D 1 D 2 D 3 = { ( 0 , 0.2 ) , ( 0.5 , 0.5 ) , ( 1 , 0.3 ) } { ( 0 , 0.4 ) , ( 1 , 0.6 ) } .
Based on Equation (8):
C = 0.2 + 0.4 2 + 0.2 + 0.6 2 + 0.5 + 0.4 2 + 0.5 + 0.6 2 + 0.3 + 0.4 2 + 0.3 + 0.6 2 = 2.5
Based on Equations (6) and (7), b 11 = ( 0 + 0 ) / 2 = 0 , b 12 = ( 0 + 1 ) / 2 = 0.5 , b 21 = ( 0.5 + 0 ) / 2 = 0.25 , b 22 = ( 0.5 + 1 ) / 2 = 0.75 , b 31 = ( 1 + 0 ) / 2 = 0.5 , b 32 = ( 1 + 1 ) / 2 = 1 . v 11 = ( ( 0.2 + 0.4 ) / 2 ) / 2.5 = 0.12 , v 12 = ( ( 0.2 + 0.6 ) / 2 ) / 2.5 = 0.16 , v 21 = ( ( 0.5 + 0.4 ) / 2 ) / 2.5 = 0.18 , v 22 = ( ( 0.5 + 0.6 ) / 2 ) / 2.5 = 0.22 , v 31 = ( ( 0.3 + 0.4 ) / 2 ) / 2.5 = 0.14 , v 32 = ( ( 0.3 + 0.6 ) / 2 ) / 2.5 = 0.18 . So D 1 D 2 D 3 = { ( 0 , 0.2 ) , ( 0.5 , 0.5 ) , ( 1 , 0.3 ) } { ( 0 , 0.4 ) , ( 1 , 0.6 ) } = { ( 0 , 0.12 ) , ( 0.25 , 0.18 ) , ( 0.50 , 0.30 ) , ( 0.75 , 0.22 ) , ( 1 , 0.18 ) } .

References

  1. National Academies of Sciences, Engineering, and Medicine. Negative Emissions Technologies and Reliable Sequestration: A Research Agenda; National Academies Press: Washington, DC, USA, 2019. [Google Scholar]
  2. Abuelgasim, S.; Wang, W.; Abdalazeez, A. A brief review for chemical looping combustion as a promising CO2 capture technology: Fundamentals and progress. Sci. Total Environ. 2021, 764, 142892. [Google Scholar] [CrossRef] [PubMed]
  3. Alalwan, H.A.; Alminshid, A.H. CO2 capturing methods: Chemical looping combustion (CLC) as a promising technique. Sci. Total Environ. 2021, 788, 147850. [Google Scholar] [CrossRef] [PubMed]
  4. Krzywanski, J.; Żyłka, A.; Czakiert, T.; Kulicki, K.; Jankowska, S.; Nowak, W. A 1.5 D model of a complex geometry laboratory scale fuidized bed clc equipment. Powder Technol. 2017, 316, 592–598. [Google Scholar] [CrossRef]
  5. Álvarez Criado, Y.; Arias Rozada, B.; Abanades García, J.C. Calcium looping CO2 capture system for back-up power plants. Energy Environ. Sci. 2017, 10, 1994–2004. [Google Scholar] [CrossRef] [Green Version]
  6. Zylka, A.; Krzywanski, J.; Czakiert, T.; Idziak, K.; Sosnowski, M.; Grabowska, K.; Prauzner, T.; Nowak, W. The 4th Generation of CeSFaMB in numerical simulations for CuO-based oxygen carrier in CLC system. Fuel 2019, 255, 115776. [Google Scholar] [CrossRef]
  7. Krzywanski, J.; Czakiert, T.; Nowak, W.; Shimizu, T.; Żyłka, A.; Idziak, K.; Sosnowski, M.; Grabowska, K. Gaseous emissions from advanced CLC and oxyfuel fluidized bed combustion of coal and biomass in a complex geometry facility: A comprehensive model. Energy 2022, 251, 123896. [Google Scholar] [CrossRef]
  8. Adanez, J.; Abad, A.; Garcia-Labiano, F.; Gayan, P.; Luis, F. Progress in chemical-looping combustion and reforming technologies. Prog. Energy Combust. Sci. 2012, 38, 215–282. [Google Scholar] [CrossRef] [Green Version]
  9. Tsai, W.H. Modeling and simulation of carbon emission-related issues. Energies 2019, 12, 2531. [Google Scholar] [CrossRef] [Green Version]
  10. Change, P.C. Global Warming of 1.5 C; World Meteorological Organization: Geneva, Switzerland, 2018. [Google Scholar]
  11. IEA. Energy Technology Perspectives 2017; IEA: Paris, France, 2017; Available online: https://www.iea.org/reports/energy-technology-perspectives-2017 (accessed on 17 February 2022).
  12. Li, Q.; Liu, G.; Cai, B.; Leamon, G.; Liu, L.C.; Chen, Z.A.; Li, X. Public awareness of the environmental impact and management of carbon dioxide capture, utilization and storage technology: The views of educated people in China. Clean Technol. Environ. Policy 2017, 19, 2041–2056. [Google Scholar] [CrossRef]
  13. Lee, S.Y.; Lee, I.B.; Han, J. Design under uncertainty of carbon capture, utilization and storage infrastructure considering profit, environmental impact, and risk preference. Appl. Energy 2019, 238, 34–44. [Google Scholar] [CrossRef]
  14. Qin, J.; Li Yl, W.; Wu, D. CCUS global progress and China’s policy suggestions. Pet. Geol. Recovery Effic. 2020, 27, 20–28. [Google Scholar]
  15. Yao, X.; Zhong, P.; Zhang, X.; Zhu, L. Business model design for the carbon capture utilization and storage (CCUS) project in China. Energy Policy 2018, 121, 519–533. [Google Scholar] [CrossRef]
  16. Li, Q.; Chen, Z.; Zhang, J.T.; Liu, L.C.; Li, X.; Jia, L. Positioning and revision of CCUS technology development in China. Int. J. Greenh. Gas Control. 2016, 46, 282–293. [Google Scholar] [CrossRef]
  17. Minchener, A.J. Gasification based CCS challenges and opportunities for China. Fuel 2014, 116, 904–909. [Google Scholar] [CrossRef]
  18. He, M.; Luis, S.; Rita, S.; Ana, G.; Euripedes, V., Jr.; Zhang, N. Risk assessment of CO2 injection processes and storage in carboniferous formations: A review. J. Rock Mech. Geotech. Eng. 2011, 3, 39–56. [Google Scholar] [CrossRef] [Green Version]
  19. Liu, L.C.; Li, Q.; Zhang, J.T.; Cao, D. Toward a framework of environmental risk management for CO2 geological storage in China: Gaps and suggestions for future regulations. Mitig. Adapt. Strateg. Glob. Chang. 2016, 21, 191–207. [Google Scholar] [CrossRef]
  20. Jiang, K.; Ashworth, P.; Zhang, S.; Liang, X.; Sun, Y.; Angus, D. China’s carbon capture, utilization and storage (CCUS) policy: A critical review. Renew. Sustain. Energy Rev. 2020, 119, 109601. [Google Scholar] [CrossRef]
  21. Viebahn, P.; Vallentin, D.; Höller, S. Prospects of carbon capture and storage (CCS) in China’s power sector—An integrated assessment. Appl. Energy 2015, 157, 229–244. [Google Scholar] [CrossRef] [Green Version]
  22. China, M. Technical Guideline on Environmental Risk Assessment for Carbon Dioxide Capture, Utilization and Storage (on Trial); Ministry of Environmental Protection of the People’s Republic of China: Beijing, China, 2016.
  23. Gerstenberger, M.C.; Christophersen, A.; Buxton, R.; Allinson, G.; Hou, W.; Leamon, G.; Nicol, A. Integrated risk assessment for CCS. Energy Procedia 2013, 37, 2775–2782. [Google Scholar] [CrossRef] [Green Version]
  24. International Organization for Standardization. Risk Management–Principles and Guidelines; International Organization for Standardization: Geneva, Switzerland, 2009; pp. 15–21. [Google Scholar]
  25. Dias, S.; Le Guen, Y.; Poupard, O.; Shtivelman, V. Risk assessment of MUSTANG project experimental site–Methodological development. Energy Procedia 2011, 4, 4109–4116. [Google Scholar] [CrossRef] [Green Version]
  26. Wollenweber, J.; Busby, D.; Wessel-Berg, D.; Nepveu, M.; Codreanu, D.B.; Grimstad, A.A.; Sijacic, D.; Maurand, N.; Lothe, A.; Wahl, F.; et al. Integrated carbon risk assessment (ICARAS). Energy Procedia 2013, 37, 4825–4832. [Google Scholar] [CrossRef] [Green Version]
  27. Hnottavange-Telleen, K. Common themes in risk evaluation among eight geosequestration projects. Energy Procedia 2013, 37, 2794–2801. [Google Scholar] [CrossRef] [Green Version]
  28. Agency, E. Environmental Risk Assessment for Carbon Capture and Storage 2011; Technical Report 40; Environment Agency: Bristol, UK, 2011.
  29. Liu, B.; Liu, S.; Xue, B.; Lu, S.; Yang, Y. Formalizing an integrated decision-making model for the risk assessment of carbon capture, utilization, and storage projects: From a sustainability perspective. Appl. Energy 2021, 303, 117624. [Google Scholar] [CrossRef]
  30. Chakraborty, S.; Chattopadhyay, R.; Chakraborty, S. An integrated D-MARCOS method for supplier selection in an iron and steel industry. Decis. Mak. Appl. Manag. Eng. 2020, 3, 49–69. [Google Scholar]
  31. Božanić, D.; Milić, A.; Tešić, D.; Salabun, W.; Pamučar, D. D numbers–FUCOM–fuzzy RAFSI model for selecting the Group of construction machines for enabling mobility. Facta Univ. Ser. Mech. Eng. 2021, 19, 447–471. [Google Scholar] [CrossRef]
  32. Dempster, A. Upper and lower probabilities induced by a multivalued mapping. Ann. Math. Stat. 1967, 38, 325–339. [Google Scholar] [CrossRef]
  33. Deng, X.; Jiang, W. A total uncertainty measure for D numbers based on belief intervals. Int. J. Intell. Syst. 2019, 34, 3302–3316. [Google Scholar] [CrossRef] [Green Version]
  34. Lai, H.; Liao, H. A multi-criteria decision making method based on DNMA and CRITIC with linguistic D numbers for blockchain platform evaluation. Eng. Appl. Artif. Intell. 2021, 101, 104200. [Google Scholar] [CrossRef]
  35. Seiti, H.; Hafezalkotob, A.; Herrera-Viedma, E. A novel linguistic approach for multi-granular information fusion and decision-making using risk-based linguistic D numbers. Inf. Sci. 2020, 530, 43–65. [Google Scholar] [CrossRef]
  36. Mi, X.; Tian, Y.; Kang, B. MADA problem: A new scheme based on D numbers and aggregation functions. J. Intell. Fuzzy Syst. 2021, 40, 11231–11255. [Google Scholar] [CrossRef]
  37. Liu, P.; Zhang, X. A novel approach to multi-criteria group decision-making problems based on linguistic D numbers. Comput. Appl. Math. 2020, 39, 1–27. [Google Scholar] [CrossRef]
  38. Mi, X.; Tian, Y.; Kang, B. A hybrid multi-criteria decision making approach for assessing health-care waste management technologies based on soft likelihood function and D-numbers. Appl. Intell. 2021, 51, 6708–6727. [Google Scholar] [CrossRef]
  39. Pamučar, D.; Puška, A.; Željko, S.; Ćirović, G. A new intelligent MCDM model for HCW management: The integrated BWM–MABAC model based on D numbers. Expert Syst. Appl. 2021, 175, 114862. [Google Scholar] [CrossRef]
  40. Mo, H. A SWOT method to evaluate safety risks in life cycle of wind turbine extended by D number theory. J. Intell. Fuzzy Syst. 2021, 40, 4439–4452. [Google Scholar] [CrossRef]
  41. Tian, Y.; Mi, X.; Liu, L.; Kang, B. A new soft likelihood function based on d numbers in handling uncertain information. Int. J. Fuzzy Syst. 2020, 22, 2333–2349. [Google Scholar] [CrossRef]
  42. Brauers, W.K.; Zavadskas, E.K. The MOORA method and its application to privatization in a transition economy. Control. Cybern. 2006, 35, 445–469. [Google Scholar]
  43. Brauers, W.K.M.; Zavadskas, E.K. Project management by MULTIMOORA as an instrument for transition economies. Technol. Econ. Dev. Econ. 2010, 16, 5–24. [Google Scholar] [CrossRef]
  44. Mi, X.; Liao, H.; Liao, Y.; Lin, Q.; Lev, B.; Al-Barakati, A. Green suppler selection by an integrated method with stochastic acceptability analysis and MULTIMOORA. Technol. Econ. Dev. Econ. 2020, 26, 549–572. [Google Scholar] [CrossRef]
  45. Wang, J.; Ma, Q.; Liu, H.C. A meta-evaluation model on science and technology project review experts using IVIF-BWM and MULTIMOORA. Expert Syst. Appl. 2021, 168, 114236. [Google Scholar] [CrossRef]
  46. Zavadskas, E.K.; Bausys, R.; Lescauskiene, I.; Usovaite, A. MULTIMOORA under interval-valued neutrosophic sets as the basis for the quantitative heuristic evaluation methodology HEBIN. Mathematics 2021, 9, 66. [Google Scholar] [CrossRef]
  47. Fedajev, A.; Stanujkic, D.; Karabašević, D.; Brauers, W.K.; Zavadskas, E.K. Assessment of progress towards “Europe 2020” strategy targets by using the MULTIMOORA method and the Shannon Entropy Index. J. Clean. Prod. 2020, 244, 118895. [Google Scholar] [CrossRef]
  48. Hafezalkotob, A.; Hafezalkotob, A.; Liao, H.; Herrera, F. Interval MULTIMOORA method integrating interval Borda rule and interval best–worst-method-based weighting model: Case study on hybrid vehicle engine selection. IEEE Trans. Cybern. 2019, 50, 1157–1169. [Google Scholar] [CrossRef] [PubMed]
  49. Rani, P.; Mishra, A.R. Fermatean fuzzy Einstein aggregation operators-based MULTIMOORA method for electric vehicle charging station selection. Expert Syst. Appl. 2021, 182, 115267. [Google Scholar] [CrossRef]
  50. Griggs, D.; Stafford-Smith, M.; Gaffney, O.; Rockström, J.; Öhman, M.C.; Shyamsundar, P.; Steffen, W.; Glaser, G.; Kanie, N.; Noble, I. Sustainable development goals for people and planet. Nature 2013, 495, 305–307. [Google Scholar] [CrossRef]
  51. Roussanaly, S.; Anantharaman, R.; Lindqvist, K.; Zhai, H.; Rubin, E. Membrane properties required for post-combustion CO2 capture at coal-fired power plants. J. Membr. Sci. 2016, 511, 250–264. [Google Scholar] [CrossRef]
  52. Ramírez, A.; Hoogwijk, M.; Hendriks, C.; Faaij, A. Using a participatory approach to develop a sustainability framework for carbon capture and storage systems in The Netherlands. Int. J. Greenh. Gas Control 2008, 2, 136–154. [Google Scholar] [CrossRef] [Green Version]
  53. Rosa, L.; Reimer, J.A.; Went, M.S.; D’Odorico, P. Hydrological limits to carbon capture and storage. Nat. Sustain. 2020, 3, 658–666. [Google Scholar] [CrossRef]
  54. Deng, X.; Hu, Y.; Deng, Y.; Mahadevan, S. Environmental impact assessment based on D numbers. Expert Syst. Appl. 2014, 41, 635–643. [Google Scholar] [CrossRef]
  55. Xia, J.; Feng, Y.; Liu, L.; Liu, D.; Fei, L. On entropy function and reliability indicator for D numbers. Appl. Intell. 2019, 49, 3248–3266. [Google Scholar] [CrossRef]
  56. Karaca, C.; Ulutaş, A. Supplier performance evaluation by using SWARA and MULTIMOORA. Econ. Manag. Econom. 2017, 137–151. [Google Scholar]
  57. Bachu, S. CO2 storage in geological media: Role, means, status and barriers to deployment. Prog. Energy Combust. Sci. 2008, 34, 254–273. [Google Scholar] [CrossRef]
Figure 1. Process of the proposed risk assessment approach.
Figure 1. Process of the proposed risk assessment approach.
Axioms 11 00204 g001
Figure 2. The risk level of each option in RSM and PRM coordinates.
Figure 2. The risk level of each option in RSM and PRM coordinates.
Axioms 11 00204 g002
Figure 3. The risk level of each option in RSM and FMF coordinates.
Figure 3. The risk level of each option in RSM and FMF coordinates.
Axioms 11 00204 g003
Figure 4. The risk level of each option in PRM and FMF coordinates.
Figure 4. The risk level of each option in PRM and FMF coordinates.
Axioms 11 00204 g004
Figure 5. MULTIMOORA values under different weights, in which figure (a) corresponds to RSM, figure (b) corresponds to PRM, figure (c) corresponds to FMF and figure (d) is their geometric mean.
Figure 5. MULTIMOORA values under different weights, in which figure (a) corresponds to RSM, figure (b) corresponds to PRM, figure (c) corresponds to FMF and figure (d) is their geometric mean.
Axioms 11 00204 g005
Table 1. The comprehensive risk indicator system for evaluating CCUS projects.
Table 1. The comprehensive risk indicator system for evaluating CCUS projects.
DimensionIndicatorMeaningReferences
Economy
Cost C 1 The economic cost of the CCUS project, including investment cost, operation and maintenance
costs, unit capture costs, payback period uncertainty, and affordability.
[51]
Market C 2 The economic system in CCUS market allocation that plays a fundamental role includes
market obstacle, market uncertainty, market competition, and market maturity.
[52]
Industrial
development
C 3 Technology, market size, and value of CO2 capture. It includes the development of the capture industry,
transportation industry, storage industry, utilization industry, and monitoring industry.
[15,23]
Society
Social
acceptance
C 4 The degree to which CCUS projects are accepted by the public, stakeholders, governments, and NGOs,
including fairness, equality, access, and so on.
[23]
Justice C 5 Including equal management between regions and generations, equal availability for regions and generations,
equal accessibility for regions and generations, and equal negative impacts for regions and generations
[52]
Social
benefits
C 6 The contribution of the CCUS project to society.Proposed
by authors
Environment
Climate
pollution
C 7 The CCUS process causes certain substances to enter the atmosphere and, thus, endanger human
health, including air pollution, soil pollution, water pollution, and so on.
[12]
Resources C 8 Resources to meet human needs in CCUS projects, including the energy consumed during the life cycle of a CCUS
and the diversity of energy sources.
[52,53]
Health and
security
C 9 Including capture safety, transportation safety, storage safety, utilization safety, monitoring safety,
alarm and management systems, impact on human health, catastrophic events, etc.
[52]
Governance
Management C 10 Participation degree of government, enterprises, research institutions, and the public in management.[16,29]
Policy and
regulation
C 11 National laws and local government support for CCUS programs, including domestic policies
and regulations and international policies and regulations.
[29]
Demonstration C 12 Demonstration to other industries in carbon capture, transportation, utilization, storage, and monitoring.[16]
Technology
Technological
advancement
C 13 The maturity, flexibility, complexity, and reliability of CCUS technology. It also includes the adoption of
technical processes and the use of equipment manufacturer process methods in CCUS.
[16]
Technological
potential
C 14 The development of carbon transport technologies includes the application and expansion of technologies,
knowledge created, and innovative breakthroughs.
[20]
Technology
management
C 15 The process of efficiently collecting, storing, processing, and applying CCUS data using computer hardware
and software technologies.
[29]
Table 2. The assessment standard for D numbers.
Table 2. The assessment standard for D numbers.
Assessment GradeNumerical RatingDescription
C + +3Extremely good
B + +2Good
A + +1Somewhat good
N0Medium
A −1Somewhat poor
B −2Poor
C −3Extremely poor
Table 3. The assessment information of C O 2 storage site selection represented by D numbers.
Table 3. The assessment information of C O 2 storage site selection represented by D numbers.
DMs/IndicatesOption 1Option 2Option 3Option 4
e 1
C 1 {(2, 0.3),(3, 0.5)}{(1, 0.4),(2, 0.6)}{(0, 0.4),(1, 0.6)}{(−1, 0.5),(2, 0.5)}
C 2 {(1, 0.4),(2, 0.6)}{(0, 0.3),(1, 0.7)}{(1, 0.6)}{(0, 0.8),(1, 0.2)}
C 15 {(0, 0.5),(1, 0.5)}{(1, 0.3),(2, 0.6)}{(1, 0.5),(2, 0.5)}{(1, 0.8),(2, 0.2)}
e 2
C 1 {(2, 0.4),(3, 0.5)}{(1, 0.2),(2, 0.8)}{(0, 0.9),(1, 0.1)}{(−1, 1.0)}
C 2 {(1, 0.7),(2, 0.3)}{(0, 0.7),(1, 0.3)}{(1, 0.7),(2, 0.2)}{(0, 0.8),(1, 0.2)}
C 15 {(0, 1.0)}{(1, 0.1),(2, 0.9)}{(1, 0.6),(2, 0.4)}{(1, 0.6),(2, 0.4)}
e 3
C 1 {(2, 0.6),(3, 0.3)}{(1, 0.7),(2, 0.3)}{(1, 1.0)}{(−1, 0.7),(2, 0.3)}
C 2 {(1, 0.5),(2, 0.5)}{(0, 0.8),(1, 0.2)}{(1, 0.3),(2, 0.6)}{(0, 1.0)}
C 15 {(0, 0.8),(1, 0.2)}{(1, 0.6),(2, 0.3)}{(1, 0.3),(2, 0.7)}{(1, 0.6),(2, 0.4)}
Table 4. The combined assessment information of C O 2 storage site selection.
Table 4. The combined assessment information of C O 2 storage site selection.
IndicatesOption 1Option 2Option 3Option 4
C 1 {(2, 0.09),(2.25, 0.26)
(2.5, 0.33),(2.75, 0.24)
(3, 0.08)}
{(1, 0.07),(1.25, 0.24)
(1.5, 0.33),(1.75, 0.26)
(2, 0.1)}
{(0, 0.15),(0.25, 0.32)
(0.5, 0.33),(0.75, 0.19)
(1, 0.01)}
{(−1, 0.14),(−0.125, 0.31)
(0.5, 0.33),(1.25, 0.21)
(2, 0.02)}
C 2 {(1, 0.1),(1.25, 0.27)
(1.5, 0.33),(1.75, 0.23)
(2, 0.07)}
{(0, 0.1),(0.25, 0.27)
(0.5, 0.33),(0.75, 0.23)
(1, 0.07)}
{(1, 0.1),(1.25, 0.26)
(1.5, 0.33),(1.75, 0.24)
(2, 0.07)}
{(0, 0.15),(0.25, 0.31)
(0.5, 0.33),(0.75, 0.19)
(1, 0.02)}
C 15 {(0, 0.15),(0.25, 0.32)
(0.5, 0.33),(0.75, 0.19)
(1, 0.01)}
{(1, 0.06),(1.25, 0.22)
(1.5, 0.33),(1.75, 0.28)
(2, 0.11)}
{(1, 0.08),(1.25, 0.25)
(1.5, 0.33),(1.75, 0.25)
(2, 0.09)}
{(1, 0.1),(1.25, 0.27)
(1.5, 0.33),(1.75, 0.23)
(2, 0.07)}
Table 5. Intermediate results of the MULTIMOORA-based risk evaluation approach.
Table 5. Intermediate results of the MULTIMOORA-based risk evaluation approach.
IndicatesOption 1Option 2Option 3Option 4 ω ω ω
C 1 2.491.520.400.290.06780.10280.1028
C 2 1.480.481.480.410.05890.05250.0456
C 3 2.361.313.281.830.05590.06650.0548
C 4 0.221.323.462.540.06610.10290.1004
C 5 0.540.793.921.920.07860.04360.0506
C 6 2.522.433.703.740.04810.07170.0509
C 7 1.414.124.244.530.03640.12370.0664
C 8 2.422.512.831.840.09280.03240.0444
C 9 1.720.383.611.540.05150.00800.0061
C 10 2.970.422.810.310.09880.08570.1249
C 11 0.840.522.283.520.09430.10280.1430
C 12 2.634.053.532.540.03690.02720.0148
C 13 2.613.383.633.500.02370.06680.0234
C 14 2.970.461.361.790.11290.08100.1349
C 15 0.391.541.501.480.07730.03240.0370
Table 6. The results from the weighted version of the MULTIMOORA method.
Table 6. The results from the weighted version of the MULTIMOORA method.
OptionsRSMPRMFMFGeometric Mean Method
Y i Ranking Z i Ranking U i Ranking S i Ranking
Option 10.261230.088930.435831.08573
Option 20.107440.099640.287640.67714
Option 30.452810.057910.821011.85901
Option 40.337620.080620.587221.34972
Table 7. Qualitative and quantitative comparison results.
Table 7. Qualitative and quantitative comparison results.
MethodMain ContributionApplicationAggregation MethodRanking
[33]Uncertainty quantification
of D numbers
Feature evaluation for
classification
- O p t i o n 3 O p t i o n 4 O p t i o n 2 O p t i o n 1
[34]LDNs, DNMA, and CRITICBlockchain platform evaluationDouble normalization-
based multiple aggregation
O p t i o n 3 O p t i o n 4 O p t i o n 1 O p t i o n 2
[36]Consider the attitudinal
feature of decision makers
Car performance assessmentPower ordered weighted
averaging operator
O p t i o n 3 O p t i o n 4 O p t i o n 2 O p t i o n 1
[38]Consider multi-granular
linguistic terminology
Health-care waste management
technologies assessment
Soft likelihood function O p t i o n 3 O p t i o n 4 O p t i o n 1 O p t i o n 2
[40]Strengths-weaknesses-
opportunities-threats
Assessment of safety risks in
life cycle of wind turbine
Traditional fusion rules
of D numbers
O p t i o n 3 O p t i o n 4 O p t i o n 2 O p t i o n 1
[41]Deal with uncertain
information of D numbers
Investigation of the
criminal case
Soft likelihood function O p t i o n 3 O p t i o n 4 O p t i o n 1 O p t i o n 2
Our methodNovel fusion rule
A index system
D_MULTIMOORA method
CCUS project risk assessmentNovel fusion rule O p t i o n 3 O p t i o n 4 O p t i o n 1 O p t i o n 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, P.; Che, L.; Wan, L.; Fei, L. A MULTIMOORA-Based Risk Evaluation Approach for CCUS Projects by Utilizing D Numbers Theory. Axioms 2022, 11, 204. https://doi.org/10.3390/axioms11050204

AMA Style

Li P, Che L, Wan L, Fei L. A MULTIMOORA-Based Risk Evaluation Approach for CCUS Projects by Utilizing D Numbers Theory. Axioms. 2022; 11(5):204. https://doi.org/10.3390/axioms11050204

Chicago/Turabian Style

Li, Peilin, Lina Che, Luhe Wan, and Liguo Fei. 2022. "A MULTIMOORA-Based Risk Evaluation Approach for CCUS Projects by Utilizing D Numbers Theory" Axioms 11, no. 5: 204. https://doi.org/10.3390/axioms11050204

APA Style

Li, P., Che, L., Wan, L., & Fei, L. (2022). A MULTIMOORA-Based Risk Evaluation Approach for CCUS Projects by Utilizing D Numbers Theory. Axioms, 11(5), 204. https://doi.org/10.3390/axioms11050204

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop