文章快速检索    
  同济大学学报(自然科学版)  2018, Vol. 46 Issue (8): 1138-1146.  DOI: 10.11908/j.issn.0253-374x.2018.08.018
0

引用本文  

毛人杰, 尤建新, 段春艳, 石涌江, 尤筱玥. 基于模糊随机层次分析法的供应商可持续评估模型改进[J]. 同济大学学报(自然科学版), 2018, 46(8): 1138-1146. DOI: 10.11908/j.issn.0253-374x.2018.08.018.
MAO Renjie, YOU Jianxin, DUAN Chunyan, SHI Yongjiang, YOU Xiaoyue. An Improved Sustainable Supplier Evaluation Model Based on Stochastic Analytical Hierarchy Process Method and Fuzzy Set Theory[J]. Journal of Tongji University (Natural Science), 2018, 46(8): 1138-1146. DOI: 10.11908/j.issn.0253-374x.2018.08.018

基金项目

国家自然科学基金(71402090,71671125,71701153);“博士后国际交流计划”派出项目(20160087)

第一作者

毛人杰(1976—),男,博士生,主要研究方向为管理理论与工业工程.E-mail: maorjb@gmail.com

通信作者

尤筱玥(1991—),女,博士生,主要研究方向为战略管理与财务管理. E-mail: yxyrachel@sina.com

文章历史

收稿日期:2017-11-20
基于模糊随机层次分析法的供应商可持续评估模型改进
毛人杰 1, 尤建新 1, 段春艳 1,2, 石涌江 3, 尤筱玥 1,3     
1. 同济大学 经济与管理学院,上海 200092;
2. 华盛顿大学 放射肿瘤系,西雅图 98195;
3. 剑桥大学 制造研究院,剑桥 CB3 0FS
摘要:利用第三方供应商可持续性评估平台,跨国企业能够更有效监督全球范围内的供应商在可持续发展中做出的努力与改进.在使用第三方平台评估供应商可持续性表现时,发现用于评价和处理指标权重和供应商指标得分的方法缺乏精确度.为此,提出采用随机层次分析法和三角模糊数理论对评估模型进行优化.为检验改进模型的有效性,采用企业实际案例,对比分析原始模型与改进模型的评估结果,发现改进模型能够将供应商可持续性表现有效差异化,得到区分度较明显的评价排序.
关键词供应商评估    可持续    随机层次分析法    三角模糊数    
An Improved Sustainable Supplier Evaluation Model Based on Stochastic Analytical Hierarchy Process Method and Fuzzy Set Theory
MAO Renjie 1, YOU Jianxin 1, DUAN Chunyan 1,2, SHI Yongjiang 3, YOU Xiaoyue 1,3     
1. College of Economics and Management, Tongji University, Shanghai 200092, China;
2. Department of Radiation Oncology, University of Washington, Seattle 98195, USA;
3. Institute for Manufacturing, University of Cambridge, Cambridge CB3 0FS, UK
Abstract: Based on a third-party platform for evaluating supplier sustainability, transnational companies can effectively monitor how much efforts and improvements their global suppliers have made in sustainability. However, methods on the basis of the platform for evaluating and calculating the weights of indexes and values suffered from a lack of accuracy. A stochastic analytical hierarchy process (SAHP) method and the triangular fuzzy number were adopted to optimize the evaluation model. Besides, a practical case about the sustainability assessment of suppliers was analyzed by the original method and the improved method separately. The comparative results show that the proposed model successfully differentiates the sustainable performance of suppliers and obtains a distinct rank.
Key words: supplier evaluation    sustainability    stochastic analytical hierarchy process(SAHP)    triangular fuzzy number    

由于自然资源的有限性以及对企业社会责任的多种关切,企业的可持续发展压力日益加剧.政府、客户、盈利组织和非营利组织都要求企业在可持续发展方面做出努力,制造可持续性更高的产品[1].1994年,Webb[2]研究了产品对环境的影响,建议通过建立环境准则来选择合适的原材料,同时注重再生利用,并提出了绿色采购的概念.随后,学术界对绿色供应链的研究愈发活跃,并于1996年由国际标准化组织推出ISO14000等一系列指标.而在企业实践方面,低碳经济、可持续发展理念也随着科技发展而逐渐融入到企业文化中,谋求产品质量保证的同时提高绿色环保能力,从而体现企业的社会责任.

产品的可持续性不仅包含了所取原材料和工艺流程的可持续,而且向前向后延展至整个产品生命周期的可持续,其中涉及的除了企业自身的可持续性改善之外,还要求供应商的协同合作.因此,有必要对企业的供应商在可持发展中做出的努力进行有效监督和完善.为此,C公司与同行业部分其他企业合作,建立了针对供应商可持续性评估的第三方平台,即ECO系统,要求同意参与评估的所有供应商根据流程提供相关证明材料和数据,并根据评估结果进行定向改进,从而促使供应商与企业向可持续发展共同迈进.

然而,考虑到ECO系统在业界的使用便捷性,初期的评价指标只有4项,并且其运算方式全部采用加权平均法,严重影响了评估结果的准确度和参考价值.因此,C公司提出对系统的评价流程进行改进,例如丰富或细化评价指标、改善运算模型.目前,针对可持续供应商选择的研究已有若干,包括模型研究[3-4]、评价指标/准则筛选[5]、不同压力对可持续供应商发展的影响等[6],这些成果将对ECO系统的完善途径提供有力的理论支持.

本文从C公司的角度,对ECO系统的供应商可持续性表现评估进行改进.首先,从ECO系统中统计总结系统流程的缺陷及其后果,包括指标内容宽泛、计算方式不准确以及供应商排序结果重复.随后,在参考若干决策模型的基础上提出采用三角模糊数来描述专家评估内容,结合层次分析法来确定评价指标的相对权重,并最终确定若干供应商的可持续性表现优劣排序.与单纯的加权平均算法相比较,改进模型的计算方式能够明显提高评估结果的准确度和可靠度,有效避免排序重复情况.

1 ECO评价系统的不足与改进 1.1 ECO评价系统

为了提高监督工作的有效性与公平性,C公司借助ECO系统完成供应商可持续性评估工作.ECO系统是独立于企业与供应商之外的、针对供应商可持续性评估的第三方平台.其优势体现于:①公平性,能够采用统一的标准对不同企业的供应商可持续性做出客观评价;②多功能性,兼具评价模块与改进模块功能,能够有效监督供应商在可持续发展中做出的努力与改进;③可共享性,平台中的供应商在加入之前需授权其评估结果的信息共享,因此,平台注册成员可以对比自己与全球范围内的供应商在可持续发展中的表现.

目前,ECO系统的评估准则由专家团队根据若干环保公约、人权宣言、以及参考文献研究制定.评估指标共有4条,包括环境(environment,ENV)、劳工和人权(labor practices,LAB)、商业道德(fair business practices,FBP)和可持续采购(sustainable procurement,SUP).在参评过程中,供应商需根据问卷提交相应的证明材料或报告,随后ECO的专家团队将以此作为依据进行完成度评价,并综合得到指标的分数,最后用加权平均的方式得到供应商可持续性评估综合得分.

1.2 不足与改进

虽然ECO系统已被广泛使用于企业可持续性评估中,但在实践时,还存在以下不足之处:

(1) 评价指标权重时未考虑指标间的相互影响关系;

(2) 评价供应商不同指标得分时难以给出精确数值.在实际评估过程中,由于供应商可持续发展意识不足,时常出现缺乏支持文件或随意上传文件的情况.如图 1所示,在随机抽取的300份供应商可持续性发展能力评估报告中,有30%的供应商上传的支持文件小于5份,其中上传0份文件的占到14.33%;

图 1 供应商提供的支持文件数量统计 Fig.1 Supporting documents provided by suppliers statistics

(3) 评估指标较少且简单采用加权平均方法计算供应商综合得分,再经四舍五入之后经常出现同分现象,造成无法判别供应商优劣.

供应商选择是一个多准则决策问题(multi-criteria decision-making,MCDM)[7].研究表明,模糊集理论(fuzzy set theory)能够有效解决MCDM评价信息的模糊性和不完整性问题.Zhao等[8]将直觉模糊集与组合权重应用于改进VIKOR模型(VIseKriterijumska Optimizacija I Kompromisno Resenje),由此提高了供应商选择过程的精确性和客观性.尤筱玥和尤建新[9]则通过区间二元语义与VIKOR模型相结合,规避了排序结果重复问题,同时避免了数据采集和处理过程中的信息缺失情况.高建伟等[10]利用区间直觉模糊集来表达和处理模糊不确定信息,并提出一种基于前景理论的随机直觉模糊决策方法.付畅俭和皮冬洋[11]基于犹豫模糊软集,提出了一种考虑个体偏好和可信度的多准则群决策方法.Gou等[12]将Atanassov直觉模糊集(Atanassov’s intuitionistic fuzzy set,A-IFS)、毕达哥斯拉模糊集(Pythagorean fuzzy set,PFS)和犹豫模糊集(hesitant fuzzy set,HFS)3种模糊数集结合,用于处理多种非精确数值的计算.

Jalao等[13]在处理MCDM问题时,将随机变量与层次分析法(analytic hierarchy process,AHP)结合,提出随机层次分析法(stochastic analytical hierarchy process,SAHP)用于计算指标的权重.Cobuloglu和Büyüktahtakn[14]在其基础上将精确值、区间值和三角数三者相结合来解决评价信息的不确定性.

为解决评价指标少、权重不合理、计算方式缺乏精确度等问题,本文在基于C公司供应商可持续性评估实际工作基础上,通过参考大量供应商可持续性评估资料,将原有的4个评估指标进行拓展,给出供应商在细分指标下的得分.同时,采用SAHP计算4个指标和细化指标的权重.该模型主要有两方面优点:①将评价指标进行分层处理,在计算权重时充分考虑指标与指标之间的相互影响;②在评价指标间的相对重要性时,决策者可以根据个人喜好选择合适的方法,如精确值、三角模糊数(triangular fuzzy number)[15]等,来表示评价值,给予决策者充分的灵活性.针对ECO系统计算方式的改进中,将采用三角模糊数来表示供应商在每个细化指标下的得分,并根据规则运算给出供应商可持续表现排序.

2 供应商可持续性评估模型 2.1 问题描述

假设C公司决定在m个供应商Ai(i=1, 2, …, m)中选出一个可持续发展能力最好的供应商,增加与该供应商的业务合作量作为激励措施.根据选定的评价指标,由l个来自不同部门的专家Ek(k=1, 2, …, l)组成的评价小组,对备选供应商的可持续表现进行客观数据收集和主观经验判断.根据专家的职位、知识及经验赋予各专家权重,用λk(k=1, 2…, l)表示.供应商可持续性评估流程如图 2所示.

图 2 供应商可持续性评估流程 Fig.2 Process for evaluating sustainability of suppliers
2.2 计算评价指标权重

在参考C公司原有的4个评价指标的基础上,结合若干关于供应商可持续性评价指标的研究成果[11, 15-19],选取适用的指标,包括污染控制能力、环境管理体系、企业绿色形象、环境竞争能力;童工使用、工作时间、公平待遇、工资福利;公平竞争、商业诚信、保密与知识产权;明确关键要求、供应商培训、监管力度.评价指标的描述详见表 1.

下载CSV 表 1 供应商可持续性评价指标描述与来源 Tab.1 Description and source of indexes for evaluating supplier sustainable performance

基于AHP模型的具体评价指标权重计算步骤如下:

(1) 建立层次结构模型.

应用AHP方法解决评价指标权重问题时,首先需要将评价问题条理化、层次化,构造出一个含有最高层(目标层)、中间层(包括准则层和指标层)、最底层(方案层)的层次结构模型.令准则层元素用(B1, B2, …, Bp)表示,指标层元素用(C1, C2, …, Cq)表示.

(2) 构建判断矩阵.

采用Saaty[20]提出的1~9比例标度法, 见表 2, 对准则以及同准则下指标进行两两重要性比较,构造出判断矩阵.假设每个评价矩阵共有t个元素,Fi(i=1, 2, …, t)参与两两比较,则判断矩阵R如下:

$ \mathit{\boldsymbol{R = }}{\left( {{r_{ij}}} \right)_{t \times t}} = \begin{array}{*{20}{c}} {}&{\;\;\;{F_1}\;\;\;\;\;{F_2}\;\;\;\;\;\; \cdots \;\;\;\;{F_t}}\\ {\begin{array}{*{20}{c}} {{F_1}}\\ {{F_2}}\\ \cdots \\ {{F_t}} \end{array}}&{\left[ {\begin{array}{*{20}{c}} 1&{{r_{12}}}& \cdots &{{r_{1t}}}\\ {1/{r_{12}}}&1& \cdots &{{r_{2t}}}\\ \cdots&\cdots&\cdots&\cdots \\ {1/{r_{1t}}}&{1/{r_{2t}}}& \cdots &1 \end{array}} \right]} \end{array} $ (1)
下载CSV 表 2 1~9比例标度及描述 Tab.2 1~9 proportional scale description

式中:rij表示元素Fi相对元素Fj的相对重要性评价值.

在给定两个元素的相对重要性评价值时,专家可以选择精确值、区间值或三角模糊数.如:当专家认为FiFj相比,前者“一定”比后者明显重要,则rij=5;当专家认为FiFj相比,前者可能和后者同等重要也可能比后者稍微重要,则rij=[1, 3];当专家认为FiFj相比,前者至少和后者同等重要,至多比后者明显重要或最有可能比后者稍微重要,则rij=[1, 3, 5].

(3) 将判断矩阵R转化为随机判断矩阵R′=(rij)

根据Jalao等[13],假设评价值rij服从概率密度函数为fij(rij|θij)的分布,若rij为精确值,则rij~fij=rij;若为三角模糊数,则rij服从三角分布,即rij~fij(rijLL, rijML, rijUL)=Tij(rijLL, rijML, rijUL);若rij为区间值,则rij服从均匀分布,即rij~fij(rijLL, rijML, rijUL)=Uhj(rijLL, rijML, rijUL).随后采用贝塔函数(beta function)来统一拟合服从这3种分布的评价值.令rij服从贝塔分布B(rij|αij, βij, rijLL, rijUL,其中(αij, βij)为形状参数,rijLL, rijUL)为位置参数,且rijLL < rij < rijULαij, βij≥1,根据矩估计法可以得到参数估计值.

$ E\left( {{r_{ij}}} \right) = r_{ij}^{LL} + \frac{{{\alpha _{ij}}}}{{{\alpha _{ij}} + {\beta _{ij}}}}\left( {r_{ij}^{UL} - r_{ij}^{LL}} \right) $ (2)
$ {\rm{Var}}\left( {{r_{ij}}} \right) = \frac{{{\alpha _{ij}}{\beta _{ij}}}}{{{{\left( {{\alpha _{ij}} + {\beta _{ij}}} \right)}^2}\left( {{\alpha _{ij}} + {\beta _{ij}} + 1} \right)}}\left( {r_{ij}^{UL} - r_{ij}^{LL}} \right) $ (3)

式中:E(rij)为评价值rij的期望;Var(rij)为评价值rij的方差.令式(2)和式(3)计算结果分别等于样本的均值与方差,即E(rij)= ${\bar r_{ij}} $,Var(rij)=Sij2,可以得到形状参数估计值如下:

$ \begin{array}{l} {\alpha _{ij}} = \left( {\frac{{{{\bar r}_{ij}} - r_{ij}^{LL}}}{{r_{ij}^{UL} - r_{ij}^{LL}}}} \right)\\ \;\;\;\;\;\;\;\;\left( {\frac{{\left( {\frac{{{{\bar r}_{ij}} - r_{ij}^{LL}}}{{r_{ij}^{UL} - r_{ij}^{LL}}}} \right)\left( {1 - \left( {\frac{{{{\bar r}_{ij}} - r_{ij}^{LL}}}{{r_{ij}^{UL} - r_{ij}^{LL}}}} \right)} \right)}}{{\frac{{{{\left( {{S_{ij}}} \right)}^2}}}{{{{\left( {r_{ij}^{UL} - r_{ij}^{LL}} \right)}^2}}}}} - 1} \right) \end{array} $ (4)
$ \begin{array}{l} {\beta _{ij}} = \left( {1 - \frac{{{{\bar r}_{ij}} - r_{ij}^{LL}}}{{r_{ij}^{UL} - r_{ij}^{LL}}}} \right)\\ \;\;\;\;\;\;\;\;\left( {\frac{{\left( {\frac{{{{\bar r}_{ij}} - r_{ij}^{LL}}}{{r_{ij}^{UL} - r_{ij}^{LL}}}} \right)\left( {1 - \left( {\frac{{{{\bar r}_{ij}} - r_{ij}^{LL}}}{{r_{ij}^{UL} - r_{ij}^{LL}}}} \right)} \right)}}{{\frac{{{{\left( {{S_{ij}}} \right)}^2}}}{{{{\left( {r_{ij}^{UL} - r_{ij}^{LL}} \right)}^2}}}}} - 1} \right) \end{array} $ (5)

综上所述,应用矩估计法将评价值rij转化为服从贝塔分布的随机变量rij,规则如下:

(1) 当rij为精确值时,rij=rijrij为精确值;

(2) 当rij~Uij(rijLL, rijUL)时,rij~B(αij=1, βij=1, aijLL, aijUL);

(3) 当rij~Tij(rijLL, rijML, rijUL)时,rij~B(αij, βij, aijLL, aijUL),其中αijβij由式(4)、式(5)计算得到,计算时$ \begin{array}{l} {{\bar r}_{ij}} = \left( {r_{ij}^{LL} + r_{ij}^{ML} + r_{ij}^{UL}} \right)/3, {\left( {{S_{ij}}} \right)^2} = \left( {{{\left( {r_{ij}^{LL}} \right)}^2} + } \right. \end{array}$$ \begin{array}{l} \left. {{{\left( {r_{ij}^{ML}} \right)}^2} + {{\left( {r_{ij}^{UL}} \right)}^2} - r_{ij}^{LL}r_{ij}^{ML} - r_{ij}^{LL}r_{ij}^{UL} - r_{ij}^{ML}r_{ij}^{UL}} \right)/18 \end{array}$.

(4) 将服从贝塔分布的随机变量rij转化为精确值bij.

根据Kerman[21],通过封闭式近似得到的贝塔分布中值m(αij, βij)如下:

$ m\left( {{\alpha _{ij}},{\beta _{ij}}} \right) \approx \frac{{{\alpha _{ij}} - 1/3}}{{{\alpha _{ij}} + {\beta _{ij}} - 2/3}} $ (6)

式中:当αijβij时,m(αij, βij)满足$ \frac{{{\alpha _{ij}} - 1}}{{{\alpha _{ij}} + {\beta _{ij}} - 2}} \le m\left( {{\alpha _{ij}}, {\beta _{ij}}} \right) \le \frac{{{\alpha _{ij}}}}{{{\alpha _{ij}} + {\beta _{ij}}}}$.

利用上述贝塔分布的中值m(αij, βij)以及Jensen[22]中的贝塔分布公式,可将rij转化为精确值bij[14]如下:

$ {b_{ij}} = r_{ij}^{LL} + m\left( {{\alpha _{ij}},{\beta _{ij}}} \right)\left( {r_{ij}^{UL} - r_{ij}^{LL}} \right) $ (7)

(5) 一致性检验.

通过计算一致性比率CR来检验判断矩阵是否具有一致性[23].通常认为,当CR < 0.1时,判断矩阵具有满意的一致性,否则需要回到(2),调整判断矩阵赋值.

(6) 算术平均法计算层次单排序.

根据判断矩阵,可以计算相对于上一层元素本层与之有联系的所有元素的权重排序如下:

$ {w_i} = \frac{1}{t}\sum\limits_{j = 1}^t {\frac{{{b_{ij}}}}{{\sum\limits_{k = 1}^t {{b_{kj}}} }}} ,\;\;\;\;i = 1,2, \cdots ,t $ (8)

(7) 计算元素综合权重.

准则层及指标层中各元素的综合权重可以计算如下:

$ {{w'}_i} = \sum\limits_{k = 1}^t {{\lambda _k}{w_i}} $ (9)

准则层和指标层所有元素的权重分别为${w_{{B_j}}}\left( {j = 1, 2, \ldots , p} \right) $${w_{{C_j}}}\left( {j = 1, 2, \ldots , q} \right) $.

2.3 供应商可持续性表现得分排序

根据层次结构指标层的所有指标(C1, C2, …, Cq),专家需要评估所有供应商在每个指标下的表现得分.

(1) 构建评价矩阵

根据专家组成员的评价意见,构建问题的评价矩阵Gk.

$ {\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{ = }}{\left( {\tilde x_{ij}^k} \right)_{m \times q}} = \begin{array}{*{20}{c}} {}&{{C_1}\;\;\;\;{C_2}\;\;\;\;\; \cdots \;\;\;\;{C_t}}\\ {\begin{array}{*{20}{c}} {{A_1}}\\ {{A_2}}\\ \cdots \\ {{A_m}} \end{array}}&{\left[ {\begin{array}{*{20}{c}} {\tilde x_{11}^k}&{\tilde x_{12}^k}& \cdots &{\tilde x_{1q}^k}\\ {\tilde x_{21}^k}&{\tilde x_{22}^k}& \cdots &{\tilde x_{2q}^k}\\ \cdots&\cdots&\cdots&\cdots \\ {\tilde x_{m1}^k}&{\tilde x_{m2}^k}& \cdots &{\tilde x_{mq}^k} \end{array}} \right]} \end{array} $ (10)

式中:${{\tilde x}^k}_{ij} = \left( {x_{ij}^{LL}, x_{ij}^{ML}, x_{ij}^{UL}} \right) $为三角模糊数,表示专家Ek针对备选供应商Ai在评价指标Cj下给出的表现评价值.评价时,使用的评价等级及对应的三角模糊数如表 3所示[19].

下载CSV 表 3 评价等级及对应三角模糊数 Tab.3 Evaluation scale and triangular fuzzy numbers

(2) 构建集体评价矩阵

利用模糊加权平均(fuzzy weighted averaging,FWA)算子聚合l位专家给出的评价值,即将l个评价矩阵Gk聚合成一个集体评价矩阵$ \mathit{\boldsymbol{G}} = {\left( {{{\tilde x}_{ij}}} \right)_{m \times q}}$.

$ \begin{array}{l} \mathit{\boldsymbol{G = }}\sum\limits_{k = 1}^l {{\lambda _k}\tilde x_{ij}^k} = \\ \;\;\;\;\;\;\left( {\sum\limits_{k = 1}^l {{\lambda _k}x_{ij}^{LL}} ,\sum\limits_{k = 1}^l {{\lambda _k}x_{ij}^{ML}} ,\sum\limits_{k = 1}^l {{\lambda _k}x_{ij}^{UL}} } \right) \end{array} $ (11)

(3) 将聚合评价值转化为精确值

利用三角函数去模糊化公式[24]将加权评价值${{\tilde x}_{ij}} = \left( {x_{ij}^{LL}, x_{ij}^{ML}, x_{ij}^{UL}} \right) $转化为精确值zij如下:

$ {z_{ij}} = \frac{{x_{ij}^{LL} + 2x_{ij}^{ML} + x_{ij}^{UL}}}{4} $ (12)

(4) 计算综合得分

根据综合权重以及供应商在各指标下的得分来计算供应商可持续表现综合得分zi如下:

$ {z_i} = \sum\limits_{j = 1}^q {{w_{{C_j}}}{z_{ij}}} $ (13)
$ {{z'}_i} = \sum\limits_{j = 1}^p {{w_{{B_j}}}{z_i}} $ (14)

式中:zi数值越大表示供应商Ai的可持续发展能力越好.最终根据zi进行降序排列,得到所有备选供应商的可持续发展能力排序.

3 案例分析 3.1 模型应用

C公司现有A1、A2和A3三个供应商提供同类型产品.根据供应商在ECO系统上问卷的回答情况及其提供的相关证明材料,需要对三个供应商的可持续表现进行评价.评价小组由三位专家组成,分别来自C公司采购部、ECO专职审计人员和C公司可持续发展部,其权重分别为λ1=0.3,λ2=0.4和λ3=0.3.

3.1.1 计算评价指标权重

由于计算过程涉及数据较多,此处仅列出ENV准则下4个指标综合权重的计算过程作为SAHP方法的算例.

(1) 建立层次结构模型,如图 3所示.

图 3 供应商可持续发展能力评价指标体系 Fig.3 Evaluating system of sustainable supplier performance

(2) 构造ENV准则下4个指标间的重要性判断矩阵RB1-C,如表 4所示.

下载CSV 表 4 RB1-C比较判断矩阵及其转化 Tab.4 Comparison of matrix RB1-C and its transformation

(3) 根据式(2)~式(5)将RB1-C比较判断矩阵中的相对重要性评价值转化为服从贝塔分布的随机变量,如表 4所示.

(4) 根据式(6)和式(7)将用贝塔函数变量表示的相对重要性评价值转为精确值,如表 4所示.

(5) 计算BB1-C精确值判断矩阵的一致性比率CR值为CR=0.027 8 < 0.1,表示该判断矩阵通过一致性检验.

(6) 根据式(8)和式(9)计算B1准则下各指标综合权重,如表 5所示.

下载CSV 表 5 B1准则下指标综合权重 Tab.5 Comprehensive weight of indexes under B1 criterion

以此类推,得到B2、B3和B4准则及其包含指标的综合权重,如表 6所示.

下载CSV 表 6 供应商可持续性评估指标综合权重 Tab.6 Comprehensive weight of indexes for sustainable supplier evaluation

表 6可知,劳工与人权(B2)以及环境(B1)准则的权重最高,分别为40.9%和34.9%,且两者之和远大于50%.可见,在供应商可持续性评价问题中,劳工与人权和环境问题上的表现在供应商可持续发展中扮演非常重要的角色.而可持续采购B4所占比重最小,表示供应商对其下游合作伙伴(供应商的供应商)的可持续发展监管能力对企业来说并不是重点,关键仍在于供应商自身的可持续性表现.同时,根据细分指标的权重结果,环境管理体系的建设是环境准则下的重点;劳工与人权中对童工使用的监督和避免最为重要,综合两者所属准则的权重,可以认为环境管理体系和童工使用是所有14项评价指标中最应当着重监督的对象.

3.1.2 供应商综合得分排序

(1) 根据式(11)可以得到供应商可持续表现的集体评价矩阵,如表 7所示.

下载CSV 表 7 供应商可持续表现聚合评价矩阵 Tab.7 Group matrix of sustainable supplier performance

(2) 根据式(12)去模糊化,并根据式(13)和式(14)计算得到备选供应商在各评价指标下的得分以及其综合分数,如表 8所示.

下载CSV 表 8 供应商可持续表现综合得分 Tab.8 Overall rating of sustainable supplier performance

根据表 8可知,3个供应商在可持续表现上的排序为:A2>A3>A1,表示供应商A2在可持续发展表现中最为突出.而供应商A1则由于在B1和B2两个重要准则上的得分较低,导致其在综合得分上明显低于其余两家供应商.

4 方法比较

ECO系统的原始评价模型主要采用加权平均计算方式,首先根据问卷的完成度评价供应商在各准则方面的表现,随后将4项得分加权平均计算得出供应商的可持续性综合分数.表 9为ECO评价小组3位专家以百分制给出的供应商可持续表现原始得分;同时给出准则权重,以1—5分递增来表示权重由低到高,并通过加权平均得到准则的综合权重, 见表 10;最后再次利用加权平均方法计算供应商可持续表现综合得分,见表 11.在系统设计之初,ECO团队为提高数据的可视化效果,采取分数梯度的方式,对于指标得分,无论是原始数据还是计算结果,均在小数点四舍五入后归整为个位数为0或5的数字,如86分记作85分,87.5分记作90分,以此类推.

下载CSV 表 9 供应商可持续表现原始评分 Tab.9 Original rating of sustainable supplier performance
下载CSV 表 10 初始准则权重评分 Tab.10 Original weight of criteria
下载CSV 表 11 供应商可持续表现原始得分 Tab.11 Original overall rating of sustainable supplier performance

表 6表 10计算结果进行比较,发现采用SAHP方法得到的B1和B2准则权重(分别为0.349和0.409)明显高于B3和B4的权重(分别为0.125和0.118),这主要由于SAHP充分考虑了指标之间的相互影响关系;而采用ECO原始评价模型时,B3的权重(0.205)与B1和B2的权重(分别为0.317和0.342)较为接近.其次,与简单的加权平均方法相比,SAHP能够为决策者提供合适的随机变量来表示评价值,免于受到精确数值的限制.最后,由于采用SAHP方法计算的过程中将4个准则进行了细分,通过评价细化的若干指标以提高评价过程的精确性,从而能够更有效地将供应商的综合得分进行差异化,帮助判别其在可持续方面的表现.

此外,通过比较表 8表 11的供应商综合得分,认为加权平均以及四舍五入的计算方式虽然便捷,但导致了大量有效信息的丢失,使原本能够差异化的供应商评价结果趋同,由此不仅容易产生同分问题而且还会导致排序结果错乱.采用三角模糊数进行计算,能够有效避免此类问题出现,并且得到区分度明显的供应商可持续性表现排序结果.因此,为保证分数梯度并避免同分现象,将模糊集理论与SAHP相结合的改进模型更能够提高ECO系统的有效性,帮助企业寻找到真正具有可持续发展能力的供应商.

5 结论

从供应商可持续性评估实例出发,先后采用SAHP和三角模糊数对评估模型的权重计算方法和供应商指标得分排序方法进行优化.针对指标单一缺乏多样性的问题,通过参考供应商可持续性评估相关文献,对原有评价指标进行拓展,并利用细化后的评价指标对供应商进行评价,得到更为精确的评价值.相较于简单的加权平均方式,SAHP能够为决策者提供合适的随机变量来表示评价值,免于受到精确数值的限制.同时,采用三角模糊数能够避免有效信息的丢失,有效解决同分问题,得到明确排序结果,提高了ECO系统的精确性和可靠性,从而帮助企业寻找到真正具有可持续发展能力的供应商.后续研究可以将SAHP方法拓宽到其他供应商评价体系中以验证该模型的普适性,同时也可以与其他随机多指标方法相结合,例如后悔值理论、前景理论等,进一步提高供应商可持续性评估流程的精确性.

参考文献
[1]
KANNAN D. Role of multiple stakeholders and the critical success factor theory for the sustainable supplier selection process[J]. International Journal of Production Economics, 2017, 195: 391
[2]
WEBB L L. Green purchasing: forging a new link in the supply chain[J]. RESOUR, 1994, 1(6): 14
[3]
FALLAHPOUR A, OLUGU E U, MUSA S N, et al. A decision support model for sustainable supplier selection in sustainable supply chain management[J]. Computers & Industrial Engineering, 2017, 105: 391
[4]
SHABANPOUR H, YOUSEFI S, SAEN R F. Future planning for benchmarking and ranking sustainable suppliers using goal programming and robust double frontiers DEA[J]. Transportation Research Part D: Transport and Environment, 2017, 50: 129 DOI:10.1016/j.trd.2016.10.022
[5]
LUTHRA S, GOVINDAN K, KANNAN D, et al. An integrated framework for sustainable supplier selection and evaluation in supply chains[J]. Journal of Cleaner Production, 2017, 140: 1686 DOI:10.1016/j.jclepro.2016.09.078
[6]
SANCHA C, LONGONI A, GIMÉNEZ C. Sustainable supplier development practices: drivers and enablers in a global context[J]. Journal of Purchasing and Supply Management, 2015, 21(2): 95 DOI:10.1016/j.pursup.2014.12.004
[7]
CHEN W, ZOU Y. An integrated method for supplier selection from the perspective of risk aversion[J]. Applied Soft Computing, 2017, 54: 449 DOI:10.1016/j.asoc.2016.10.036
[8]
ZHAO J S, YOU X Y, LIU H C, et al. An extended VIKOR method using intuitionistic fuzzy sets and combination weights for supplier selection[J]. Symmetry, 2017, 9(9): 169 DOI:10.3390/sym9090169
[9]
尤筱玥, 尤建新. 基于区间二元语义VIKOR的外包供应商选择模型[J]. 同济大学学报(自然科学版), 2017, 45(9): 1407
YOU Xiaoyue, YOU Jianxin. Outsourcing supplier selection by interval 2-tuple linguistic VIKOR method[J]. Journal of Tongji University (Natural Science), 2017, 45(9): 1407 DOI:10.11908/j.issn.0253-374x.2017.09.022
[10]
高建伟, 刘慧晖, 谷云东. 基于前景理论的区间直觉模糊多指标决策方法[J]. 系统工程理论与实践, 2014, 34(12): 3175
GAO Jianwei, LIU Huihui, GU Yundong. Interval-valued intuitionistic fuzzy multi-criteria decision-making method based on prospect theory[J]. Systems Engineering—Theory & Practice, 2014, 34(12): 3175 DOI:10.12011/1000-6788(2014)12-3175
[11]
付畅俭, 皮冬洋. 基于个体偏好和可信度的犹豫模糊软集多指标决策方法[J]. 统计与决策, 2017(8): 37
FU Changjian, PI Dongyang. Hesitant fuzzy soft sets multicriteria decision-making method based on preference and credibility levels[J]. Statistics and Decision, 2017(8): 37
[12]
GOU X, XU Z, LIAO H. Alternative queuing method for multiple criteria decision making with hybrid fuzzy and ranking information[J]. Information Sciences, 2016, 357: 144 DOI:10.1016/j.ins.2016.03.046
[13]
JALAO E R, WU T, SHUNK D. A stochastic AHP decision making methodology for imprecise preferences[J]. Information Sciences, 2014, 270: 192 DOI:10.1016/j.ins.2014.02.077
[14]
COBULOGLU H I, BVYVKTAHTAKIN E. A stochastic multi-criteria decision analysis for sustainable biomass crop selection[J]. Expert Systems with Applications, 2015, 42(15-16): 6065 DOI:10.1016/j.eswa.2015.04.006
[15]
苏世彬, 黄瑞华. 基于三角模糊数的属性层次模型[J]. 系统工程理论与实践, 2006, 26(12): 115
SU Shibin, HUANG Ruihua. Attribute hierarchical mode based on triangular fuzzy number[J]. Systems Engineering-Theory & Practice, 2006, 26(12): 115 DOI:10.3321/j.issn:1000-6788.2006.12.017
[16]
袁宇, 关涛, 闫相斌, 等. 基于混合VIKOR方法的供应商选择决策模型[J]. 控制与决策, 2014, 29(3): 551
YUAN Yu, GUAN Tao, YAN Xiangbin, et al. Based on hybrid VIKOR method decision making model for supplier selection[J]. Control and Decision, 2014, 29(3): 551
[17]
AKMAN G. Evaluating suppliers to include green supplier development programs via fuzzy c-means and VIKOR methods[J]. Computers & Industrial Engineering, 2015, 86: 69
[18]
TSENG M L, CHIU A S F. Evaluating firm's green supply chain management in linguistic preferences[J]. Journal of Cleaner Production, 2013, 40: 22 DOI:10.1016/j.jclepro.2010.08.007
[19]
QIN J, LIU X, PEDRYCZ W. An extended TODIM multi-criteria group decision making method for green supplier selection in interval type-2 fuzzy environment[J]. European Journal of Operational Research, 2017, 258(2): 626 DOI:10.1016/j.ejor.2016.09.059
[20]
SAATY T L. A scaling method for priorities in hierarchical structures[J]. Journal of Mathematical Psychology, 1977, 15(3): 234 DOI:10.1016/0022-2496(77)90033-5
[21]
KERMAN J. A closed-form approximation for the median of the beta distribution[J/OL]. [2017-05-06]https://arxiv.org/pdf/1111.0433v1.pdf.
[22]
JENSEN P A. Continuous Distributions[J/OL]. [2017-05-08] https://www.me.utexas.edu/~jensen/ORMM/computation/unit/rvadd/continuous_dist/beta.html.
[23]
雷星晖, 尤筱玥. 基于层次分析法支持决策的外包服务供应商绩效评价[J]. 同济大学学报(自然科学版), 2014, 42(11): 1770
LEI Xinghui, YOU Xiaoyue. Evaluation of outsourcing service providers' performance on analytic hierarchy process to support decision-making[J]. Journal of Tongji University (Natural Science), 2014, 42(11): 1770 DOI:10.11908/j.issn.0253-374x.2014.11.024
[24]
王欣荣, 樊治平. 一种模糊有序加权(FOWA)算子及其应用[J]. 模糊系统与数学, 2003, 17(4): 67
WANG Xinrong, FAN Zhiping. A fuzzy ordered weighted average (FOWA) operator and its application[J]. Fuzzy Systems and Mathematics, 2003, 17(4): 67 DOI:10.3969/j.issn.1001-7402.2003.04.010