一种多阶段网络数据包络分析模型评估包装绿色度的方法与流程
本发明涉及包装科学与技术领域，尤其是指一种进行商品包装环境影响评估的方法，该方法是利用多阶段网络数据包络分析(multi-stagenetworkdataenvelopmentanalysis,简称mndea)模型评估包装绿色度；属于商品包装环境影响分析技术领域。
背景技术：
包装一直是现代社会经济与生活不可或缺的构成要素。随着包装行业的发展，出现了一些不合理的现象，比如欠包装导致内装物的损坏和腐败，过度包装则消耗太多的资源而且产生太多的包装废弃物，导致资源浪费和环境污染。因此，在满足经济发展以及人们工作生活需求的同时，最大限度的减少包装对环境的影响，促进包装与环境友好可持续发展，成了包装领域的重要课题。因此，对绿色包装的评价成了研究人员广泛关注的问题。为方便计算与分析，在本发明中将包装绿色度定义为包装全生命周期中每个阶段的环境效率值与阶段环境影响因子乘积之和。
随着社会经济可持续化发展，包装绿色度的评价已经成为包装行业的研究热点。目前，国际上公认的包装绿色度评价方法是：生命周期评价方法(lifecycleanalysis,lca)。但该评价方法在科学研究领域中仍然存在一些缺陷，主要包括：(1)在lca评价前需要主观确定权重，权重不同得到的评价结果可能会完全不同，导致评价的标准不统一，各评价体系之间不具有可比性。(2)评估过程复杂，耗时较长，难以及时反映市场和技术的动态变化等。引入一种新型的包装评价模型对促进绿色包装及包装评价的发展具有重要意义。本发明的优点是：(1)无需对数据进行无量纲化处理，无需任何权重假设，以决策单元每阶段的实际数据求得最优权重，结果更客观；(2)分阶段评估商品包装系统的环境影响可以清晰的量化各阶段的环境效率值；(3)根据各阶段的环境效率值，针对各阶段具体的环境影响因子给出具体的改善措施；(4)ndea(网络数据包络分析)模型求解分析途径广泛，求解过程简单。在包装整个生命周期的各个阶段中合理分配资源，对于实现科学发展观及可持续发展战略都有十分重要的现实意义。
dea(数据包络分析)的基本原理：每一个包装系统作为一个决策单元(dmu,decisionmakingunits)，再由多个dmu构成被评价群体，通过对投入和产出比率的综合分析，以dmu的各个投入和产出指标的权重为变量进行评价运算，确定有效生产前沿面，并根据各dmu与有效生产前沿面的距离状况，确定各dmu是否dea有效。数据包络分析方法实际上是一种凸分析方法，被评价的所有有效决策单元(dmu)即效率值为1的dmu构成一个有效前沿包络面，无效率dmu即效率值介于0～1之间的dmu则被包裹在面内。
传统dea的公式：
产出投入比率：
令
minθλ≥0
λ≥0
i＝1，2，…，m；r＝1，2，…，q；j＝1，2，…，j
传统dea模型将整个系统决策过程看作一个“黑箱”，不对其内部的结构进行分析，因此，直接应用dea对产品包装生命周期进行环境影响评估时，其分析结果不全面、不清晰，因此很有必要对此加以完善。
通过专利检索没发现有与本发明相同技术的专利文献报道，与本发明有一定关系的专利主要有以下几个：
1、专利号为cn200910022627.6，名称为“一种建筑部品绿色度的评价方法及评价装置”，申请人为：西安建筑科技大学的发明专利，该专利公开了一种建筑部品绿色度的评价方法及使用该评价方法的评价装置，该方法包括以下步骤：进行市场调查及资料查阅，将常用建筑部品分类；运用wbs方法逐层分解评价指标，建立绿色建筑部品评价指标体系，并收集和计算典型建筑部品的原始数据，建立建筑部品信息库；录入待评价建筑部品的数据信息；调整及修改指标权重，即根据行业标准数据判断是否调整相关内容，确定建筑部品绿色度评价指标的权重；建立及调用指标效用函数；对待评价建筑部品的绿色度进行评价；得出评价结果，并将评价结果进行数据保存。该专利是从建筑部品的全生命周期出发，采用定量分析与定性分析相结合，最终实现建筑部品绿色化评价方法的建立。
2、专利号为cn201220534304.2，名称为“一种绿色度评价系统”，申请人为：北京中建建筑科学研究院有限公司；中国矿业大学的实用新型专利，该专利公开了一种绿色度评价系统，包括：评分标准模块，所述评分标准模块用于对各影响因素进行评分；层次结构模块，所述层次结构模块与评分标准模块连接并且用于对各项影响因素的评分进行权重分配从而得到最后评价分值和相对应的绿色度；隶属度模块，所述隶属度模块与评分标准模块连接并且用于将得到的评价分值的隶属度呈现出来；以及结果输出模块，所述结果输出模块与隶属度模块连接并且用于对评价分值、绿色度和隶属度进行显示。
3、专利号为cn201580060181.x，名称为“用于网络分析和报告的系统和方法”，申请人为：卡特伯德网络股份有限公司的发明专利，该专利公开了一种用于网络分析和报告的系统和方法，与常规系统相比可以从多个源收集并分析资产和网络数据，并且使用此类数据来呈现各种系统和软件应用程序之间的网络连接的更完整且准确的表示及规定网络上的安全控制的操作的策略。包括由计算机系统从多个不同类型的源收集数据，其中，收集的数据包括网络数据和资产数据；由计算机系统基于所收集的数据中的网络数据来识别网络业务事件和与网络业务事件有关的多个网络资产；由计算机系统基于所收集的数据中的资产数据来识别所述多个网络资产之间的连接；由计算机系统生成流程信息图，其描绘所述多个网络资产和所述多个网络资产之间的连接，其中，该流程信息图描绘在网络资产之间被允许的网络业务和在网络资产之间被阻止的网络业务；以及经由与计算机系统通信的用户接口的显示器来呈现流程信息图。
通过对上述这些专利的仔细分析，这些专利虽然都涉及到绿色度评价以及网络分析，也提出了一些改进技术方案，但通过仔细分析，所提出的这些专利都完全没有考虑如何解决传统dea模型将整个系统决策过程看作一个“黑箱”，不对其内部的结构进行分析，因此，直接应用dea对产品包装生命周期进行环境影响评估时，其分析结果不全面、不清晰的问题，所以前面所述的问题依然存在，仍有待进一步加以研究。
技术实现要素：
本发明的目的在于针对现有包装绿色度评价中存在的问题，提供一种简单，可靠，实用的基于多阶段网络数据包络分析模型的包装绿色度评价方法。该方法理引入多阶段网络数据包络分析模型，将“黑箱”打开，同时对内部结构进行包装绿色度评估，从而准确定位导致包装绿色度较低的具体阶段，为优化包装设计提供决策支持。
为了达到这一目的，本发明提供了一种多阶段网络数据包络分析模型评估包装绿色度的方法，首先基于生命周期评价法，建立商品包装系统环境影响的多输入多输出多阶段评价指标体系，再建立商品包装系统环境影响网络数据包络评估模型；然后基于网络数据包络分析法，进行模型求解，获得商品包装系统各阶段的环境影响值；再利用统计分析法计算商品包装绿色度；分析计算结果，得出评价结论。
进一步地，所述的建立商品包装系统环境影响的多输入多输出多阶段评价指标体系是先根据评估目的及需求，选择商品包装系统作为被评估对象(即决策单元)，确定被评估对象的个数及功能单位，再基于生命周期评价法建立商品包装系统环境影响的多输入多输出多阶段评价指标体系。
进一步地，所述的确定被评估对象的个数及功能单位是根据评估目的及需求确定被评估的包装对象(决策单元)，然后给定评估的系统功能单位，在相同的对比环境下比对，这样才有可比性，否则就失去了对比意义；包装对象(决策单元)的个数j(j＝1，2，…，j)为输入输出指标及中间变量总数的2-3倍，以保证统计分析结果的准确性。
进一步地，多阶段网络数据包络分析模型评估包装绿色度的方法包括以下步骤：
1)基于生命周期评价法建立商品包装系统环境影响的多输入多输出多阶段评价指标体系；
2)建立商品包装系统环境影响的mndea评估模型；
3)基于ndea利用计算机编程进行模型求解，获得商品包装系统各阶段的环境影响值；
4)将包装对环境的影响分为两部分，一部分是整个生命过程中资源能源的消耗，另一部分是整个生命过程中的环境排放；因此，以ωp(阶段p＝1，2，…，n)代表各阶段的环境影响因子，利用每个阶段的外部输入与非期望产出(环境排放)与整个系统的外部输入及非期望产出总值的比值计算得到；
5)利用统计分析法计算商品包装绿色度；
6)分析计算结果，得出评价结论。
进一步地，将包装绿色度定义为包装生命周期中每个阶段的环境效率值与阶段环境影响因子乘积的和，记为β，则
进一步地，所述的步骤1)建立商品包装系统环境影响的多输入多输出多阶段评价指标体系是选取各阶段的评估指标选取；包括每个阶段的外部输入指标，期望产出指标(包含中间变量)及非期望产出指标(环境排放)选取：各项指标选取内容包括：
1)各阶段的外部输入指标包括：资源及能源的消耗量(其中包含：原材料获取及加工的工艺设备的消耗，包装生产及装配过程中工艺设备的消耗，包装分销及各环节的运输消耗，包装废弃处理及再生过程的消耗)；
2)期望产出指标包括：将包装生产及装配过程设定为p(p＝1，2，…，n)阶段，在p-1阶段时指的是中间变量(包装系统的过程产品)，在p阶段则是产出指标，包括：包装再利用量，回收率，保质期等；
3)非期望产出指标包括：环境排放量(例如包括co2，so2，no2等物质的排放量)。
进一步地，所述的步骤2)中建立商品包装系统环境影响的mndea评估模型是指选择链式串行mndea模型，并且中间变量作为前一阶段的产出同时也作为后一阶段的输入，公式如下：
其目标函数为：
maxβ＝w1gθ1+w2gθ2+…+wpgθp
式中：ω1为阶段1环境影响因子；ω2为阶段2环境影响因子；ωp为阶段p环境影响因子。
约束条件包括每个阶段的环境效率值及环境影响因子：
环境效率值：
第一阶段的环境效率值：
在这里所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标。第二阶段的环境效率值：
同上，所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标，第一阶段的过程输出指标也作为第二阶段的输入指标，记为中间变量。
第p阶段的环境效率值：
同上，所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标，第p-1阶段的过程输出指标也作为第p阶段的输入指标，记为中间变量。
边界约束条件：
环境影响因子：
第一阶段的环境影响因子：
第一阶段的资源能源消耗量及环境排放量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第一阶段的环境影响因子。
第二阶段的环境影响因子：
第二阶段的资源能源消耗量，环境排放量及第一阶段产生的中间变量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第二阶段的环境影响因子。
第p阶段的环境影响因子：
第p阶段的资源能源消耗量，环境排放量及第p-1阶段产生的中间变量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第p阶段的环境影响因子。
包装绿色度：
将包装全生命周期分为p阶段，每个阶段的环境影响值与阶段环境影响因子乘积之和即为包装绿色度：
将上述公式整理得出，包装绿色度的公式为：
进一步地，所述的步骤3)中基于ndea利用计算机编程进行模型求解是指，各阶段的评估指标选取中的外部输入指标，包括资源能源的消耗，在阶段p-1中其前一个阶段的期望输出也作为后一个阶段的部分输入，即中间变量，在阶段p中的期望输出作为整个系统的总输出，各阶段的非期望输出指的是环境排放；构建商品包装系统环境影响的mndea评估模型的相关变量假设如下：
yrj代表第j(j＝1，2，…，j)个商品包装对象在阶段p(p＝3，4，…，n)的第rp(rp＝1，2，…，q)种期望产出(也代表整个包装系统的产出)；
uj代表阶段p(p＝3，4，…，n)的组合权重系数；
θp代表商品包装在阶段p(p＝1，2，…，n)产生的环境效率值；
ωp代表商品包装在阶段p(p＝1，2，…，n)的环境影响因子；
β为整个商品包装系统的包装绿色度。
进一步地，所述的各阶段的评估指标选取中的外部输入指标包括：
进一步地，所述的步骤4)中，基于ndea利用计算机编程进行模型求解是指采用链式串行mndea模型进行包装绿色度的计算，公式如下：
其目标函数为：
maxβ＝w1gθ1+w2gθ2+…+wpgθp
约束条件包括每个阶段的环境效率值及环境影响因子：
环境效率值：
第一阶段的环境效率值：
在这里所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标。第二阶段的环境效率值：
同上，所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标，第一阶段的过程输出指标也作为第二阶段的输入指标，记为中间变量。
第p阶段的环境效率值：
同上，所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标，第p-1阶段的过程输出指标也作为第p阶段的输入指标，记为中间变量。
边界约束条件：
环境影响因子：
第一阶段的环境影响因子：
第一阶段的资源能源消耗量及环境排放量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第一阶段的环境影响因子。
第二阶段的环境影响因子：
第二阶段的资源能源消耗量，环境排放量及第一阶段产生的中间变量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第二阶段的环境影响因子。
第p阶段的环境影响因子：
第p阶段的资源能源消耗量，环境排放量及第p-1阶段产生的中间变量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第p阶段的环境影响因子。
包装绿色度：
将包装全生命周期分为p阶段，每个阶段的环境影响值与阶段环境影响因子乘积之和：
将上述公式整理得出，包装绿色度的公式为：
即为包装绿色度。
本发明的优点在于：
1、本发明引入多阶段网络数据包络分析模型，将“黑箱”打开，同时对内部结构进行包装绿色度评估，从而准确定位导致包装绿色度较低的具体阶段，为优化包装设计提供决策支持。
2、本发明方法理论体系清晰，相比其他方法更加客观、直观，且易于计算、使用方便，能比较完整、科学地反映整个商品包装系统各个阶段的环境影响值，以便针对环境影响值高的阶段进行进一步的改善。
附图说明
图1是本发明所述的多阶段商品包装系统流程示意图；
图2是本发明所述的mndea构建包装绿色度评价模型的示意图。
具体实施方式
下面结合附图和具体实施例来进一步阐述本发明。
实施例一
通过附图1和2可以看出，本发明涉及一种mndea模型评估包装绿色度的方法，其具体实施步骤如下：
1)根据评估目的及需求，选择商品包装系统作为被评估对象(即决策单元)，确定被评估对象的个数及功能单位；
假设本实施例是对不同形式的包装包装1000l牛奶的环境影响进行评价，检验其包装绿色度的大小。下面选择了三种包装形式：纸塑复合(低密度聚乙烯-纸板-低密度聚乙烯3层复合材料)牛奶包装盒和塑料(低密度聚乙烯材料)牛奶包装盒，纸塑铝(纸塑铝6层复合材料)复合牛奶包装盒为待评估包装单元，前两种包装的容量为250ml，所以分别共需包装4000个，纸塑铝复合包装盒的容量为1l，共需1000个包装，具体如表1所示。为使数据集更加完备充分，确保结果的准确性，采用三次埃尔米特插值法进行数据集扩充，将三个决策单元扩充到五个。若在实际生产应用中具有完备充分的数据集时，则可省略插值计算。
表1三种包装形式的具体信息
2)基于生命周期评价法建立商品包装系统环境影响的多输入多输出多阶段评价指标体系；
其中本实施例选取牛奶包装生命周期的原材料获取，包装生产及包装废弃处理三阶段进行研究，其中各阶段的投入产出指标选取参考与包装相关的lca评价以及其他dea评价如下所示：
1)投入指标：第一阶段的投入指标包括原油消耗量为
2)非期望产出(环境排放)指标：第一阶段的非期望产出指标包括二氧化碳排放量为
3)中间变量指标：第一阶段到第二阶段的中间变量为原材料产量为
4)期望产出指标：在第三阶段中的期望产出指标变量为再生原材料量为y3。
表2三阶段ndea模型投入变量、产出变量相关数据
4)基于网络数据包络分析法利用matlab编程进行模型求解，获得商品包装系统各阶段的环境影响值；
以塑料牛奶包装盒，纸塑包装及纸塑复合材料包装盒为ndea中的决策单元：塑料牛奶包装(dmu1)、纸塑包装(dmu3)、纸塑复合材料牛奶包装(dmu5)，为了保证包装绿色度的准确性利用三次埃尔米特插值补充了两组dmu数据，运用matlab2015b进行编程计算，对三个真实包装对象及两个虚拟包装对象的综合包装绿色度进行评估。
5)将包装对环境的影响分为两部分，一部分是整个生命过程中输入指标(资源能源等)的消耗，另一部分是整个生命过程中的环境排放。因此，ωp代表各阶段的环境影响因子，利用每个阶段的输入量与非期望产出(环境排放)与整个系统的外部输入及非期望产出总值的比值计算得到，见表3：
表3各阶段的环境影响度ωp
6)利用统计分析法计算商品包装绿色度；
首先将表2中的数据输入多阶段网络数据包络分析模型的matlab程序中进行计算，得出各决策单元的包装绿色度值，见表4：
表4ndea包装绿色度评价模型的包装绿色度评估结果及阶段环境效率值
7)分析计算结果，得出评价结论；
由表3可以看出塑料牛奶包装盒(dmu1)的包装绿色度最高，不管是原材料获取还是包装生产及废弃物处理阶段相较其他单元都是有效的，由于该包装形式需要的原材料较少，包装生产中的材料消耗量也较少，而且塑料包装的包装废弃物处理阶段的回收质量也不是最低的，所以导致其它包装单元的包装绿色度较低。纸塑铝复合牛奶包装盒(dmu5)及纸塑复合牛奶包装盒(dmu3)的包装绿色度值都小于1，而且纸塑复合牛奶包装盒(dmu3)在包装废弃物处理阶段的环境效率值最低，仅有0.242519，而且影响因子所占比重较大。纸塑铝复合牛奶包装盒(dmu5)则包装生产阶段的环境效率值最低，仅有0.18146，而且影响因子所占比重也较大，这些都会给包装绿色度带来影响。所以这两种包装的包装绿色度都有待进一步提高。原因有(1)塑料包装及纸塑铝复合包装材料在生产过程中消耗大量的纸板，生产纸板时会产生大量环境排放。(2)塑料包装及纸塑复合包装盒的回收率都有待提高，而且这种材料的再利用依靠成熟的纸塑分离技术，但是我国的该项技术还在发展中，并不是很成熟。因此，评估结果表明该方法在包装评价中具有很好的可行性。
上述所列实施例，只是结合附图对本发明的技术方案进行清楚、完整的描述；显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
通过上述实施例的描述，可以得知，本发明涉及一种多阶段网络数据包络分析模型评估包装绿色度的方法，首先基于生命周期评价法，建立商品包装系统环境影响的多输入多输出多阶段评价指标体系，再建立商品包装系统环境影响网络数据包络评估模型；然后基于网络数据包络分析法，进行模型求解，获得商品包装系统各阶段的环境影响值；再利用统计分析法计算商品包装绿色度；分析计算结果，得出评价结论。
进一步地，所述的建立商品包装系统环境影响的多输入多输出多阶段评价指标体系是先根据评估目的及需求，选择商品包装系统作为被评估对象(即决策单元)，确定被评估对象的个数及功能单位，再基于生命周期评价法建立商品包装系统环境影响的多输入多输出多阶段评价指标体系。
进一步地，所述的确定被评估对象的个数及功能单位是根据评估目的及需求确定被评估的包装对象(决策单元)，然后给定评估的系统功能单位，在相同的对比环境下比对，这样才有可比性，否则就失去了对比意义；包装对象(决策单元)的个数j(j＝1，2，…，j)为输入输出指标及中间变量总数的2-3倍，以保证统计分析结果的准确性。
进一步地，多阶段网络数据包络分析模型评估包装绿色度的方法包括以下步骤：
1)基于生命周期评价法建立商品包装系统环境影响的多输入多输出多阶段评价指标体系；
2)建立商品包装系统环境影响的mndea评估模型；
3)基于ndea利用计算机编程进行模型求解，获得商品包装系统各阶段的环境影响值；
4)将包装对环境的影响分为两部分，一部分是整个生命过程中输入指标(资源能源)的消耗，另一部分是整个生命过程中的环境排放；因此，ωp(阶段p＝1，2，…，n)代表各阶段的环境影响因子，利用每个阶段的输入量与非期望产出(环境排放)与整个系统的外部输入及非期望产出总值的比值计算得到；
5)利用统计分析法计算商品包装绿色度；
6)分析计算结果，得出评价结论。
进一步地，将包装绿色度定义为包装生命周期中每个阶段的环境效率值与阶段环境影响因子乘积的和，记为β，则
进一步地，所述的步骤1)建立商品包装系统环境影响的多输入多输出多阶段评价指标体系是选取各阶段的评估指标选取；包括每个阶段的外部输入指标，期望产出指标(包含中间变量)及非期望产出指标(环境排放)选取：各项指标选取内容包括：
1)各阶段的外部输入指标包括：资源及能源的消耗量(其中包含：原材料获取及加工的工艺设备的消耗，包装生产及装配过程中工艺设备的消耗，包装分销及各环节的运输消耗，包装废弃处理及再生过程的消耗)；
2)期望产出指标包括：将包装生产及装配过程设定为p(p＝1，2，…，n)阶段，在p-1阶段时指的是中间变量(包装系统的过程产品)，在p阶段则是产出指标，包括：包装再利用量，回收率，保质期等；
3)非期望产出指标包括：环境排放量(例如co2，so2，no2等物质的排放量)。
进一步地，所述的步骤2)中建立商品包装系统环境影响的mndea评估模型是指选择链式串行mndea模型，并且中间变量作为前一阶段的产出同时也作为后一阶段的输入，公式如下：
其目标函数为：
maxβ＝w1gθ1+w2gθ2+…+wpgθp
约束条件包括每个阶段的环境效率值及环境影响因子：
环境效率值：
第一阶段的环境效率值：
在这里所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标。第二阶段的环境效率值：
同上，所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标，第一阶段的过程输出指标也作为第二阶段的输入指标，记为中间变量。
第p阶段的环境效率值：
同上，所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标，第p-1阶段的过程输出指标也作为第p阶段的输入指标，记为中间变量。
边界约束条件：
环境影响因子：
第一阶段的环境影响因子：
第一阶段的资源能源消耗量及环境排放量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第一阶段的环境影响因子。
第二阶段的环境影响因子：
第二阶段的资源能源消耗量，环境排放量及第一阶段产生的中间变量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第二阶段的环境影响因子。
第p阶段的环境影响因子：
第p阶段的资源能源消耗量，环境排放量及第p-1阶段产生的中间变量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第p阶段的环境影响因子。
包装绿色度：
将包装全生命周期分为p阶段，每个阶段的环境影响值与阶段环境影响因子乘积之和：
将上述公式整理得出，包装绿色度的公式为：
即为包装绿色度。
进一步地，所述的步骤3)中基于ndea利用计算机编程进行模型求解是指，各阶段的评估指标选取中的外部输入指标，包括资源能源的消耗，在阶段p-1中其前一个阶段的期望输出也作为后一个阶段的部分输入，即中间变量，在阶段p中的期望输出作为整个系统的总输出，各阶段的非期望输出指的是环境排放；构建商品包装系统环境影响的mndea评估模型的相关变量假设如下：
yrj代表第j(j＝1，2，…，j)个商品包装对象在阶段p(p＝3，4，…，n)的第rp(rp＝1，2，…，q)种期望产出(也代表整个包装系统的产出)；
uj代表阶段p(p＝3，4，…，n)的组合权重系数；
θp代表商品包装在阶段p(p＝1，2，…，n)产生的环境效率值；
ωp代表商品包装在阶段p(p＝1，2，…，n)的环境影响因子；
β为整个商品包装系统的包装绿色度。
进一步地，所述的各阶段的评估指标选取中的外部输入指标包括：
进一步地，所述的步骤4)中，基于ndea利用计算机编程进行模型求解是指采用链式加性mndea模型进行包装绿色度的计算，公式如下：
其目标函数为：
maxβ＝w1gθ1+w2gθ2+…+wpgθp
约束条件包括每个阶段的环境效率值及环境影响因子：
环境效率值：
第一阶段的环境效率值：
在这里所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标。
第二阶段的环境效率值：
同上，所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标，第一阶段的过程输出指标也作为第二阶段的输入指标，记为中间变量。
第p阶段的环境效率值：
同上，所有的环境排放等非期望产出约束都被当作另一种形式的消耗转化为输入指标，第p-1阶段的过程输出指标也作为第p阶段的输入指标，记为中间变量。
边界约束条件：
环境影响因子：
第一阶段的环境影响因子：
第一阶段的资源能源消耗量及环境排放量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第一阶段的环境影响因子。
第二阶段的环境影响因子：
第二阶段的资源能源消耗量，环境排放量及第一阶段产生的中间变量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第二阶段的环境影响因子。
第p阶段的环境影响因子：
第p阶段的资源能源消耗量，环境排放量及第p-1阶段产生的中间变量的和与包装全生命周期共p阶段的总资源能源消耗量，总环境排放量及作为输入量的总中间变量之和的比值，即为第p阶段的环境影响因子。
包装绿色度：
将包装全生命周期分为p阶段，每个阶段的环境影响值与阶段环境影响因子乘积之和：
将上述公式整理得出，包装绿色度的公式为：
即为包装绿色度。
本发明的优点在于：
1、本发明引入多阶段网络数据包络分析模型，将“黑箱”打开，同时对内部结构进行包装绿色度评估，从而准确定位导致包装绿色度较低的具体阶段，为优化包装设计提供决策支持。
2、本发明方法理论体系清晰，相比其他方法更加客观、直观，且易于计算、使用方便，能比较完整、科学地反映整个商品包装系统各个阶段的环境影响值，以便针对环境影响值高的阶段进行进一步的改善。
- 还没有人留言评论。精彩留言会获得点赞！