纯净水安全检测论文

纯净水安全监控问题的研究与应用摘要：本文我们根据某地区提供的近年来关于各公司的纯净水检测报告，考虑电导率、菌落总数、大肠困群、霉菌四个危害因素对纯净水的安全风险评价的影响建立了数学模型，并利用此模型对纯净水安全监控问题做进一步探讨。针对问题一，提出简单的线性综合评价模型。将数据按各批次进行不同因素的加权均值处理，再无量纲化后，利用变异系法求得权向量,再根据本文已知建立特征矩阵，得到权向量，采用综合法得到四个危险指标的权向量。从而确定纯净水安全风险度模型.利用此模型对所有批次的纯净水进行安全风险评估。针对问题二，将数据按各公司进行不同因素的加权均值处理，利用模型一对该城区的九家水公司按风险度进行排序评价，；然后再按各公司的不同因素的加权均值，分类排序评价得表六，并分析表六中的排序名次（排在前面的危害因素危害程度更大）得出各公司相对其他公司而言的主要可能危害因素。依据公司主要危害因素、各公司的不同因素的加权均值数据、公司风险度、饮用水安全方面有关知识，指出公司的实际特点。针对问题三，我们先对纯净水生产流通环节即仓库和销售网点两类的不合格数进行统计计算出不合格危害因素在仓库的比例占76.19%，销售网点占3.15%，作直方图，可以直观的看出仓库比销售点存在更大的隐患。利用excel对各批次不合格危害因素进行统计。再用Excel将表一的无量纲化后的数据绘制成图表，得到四个危害指标随着检验批次的变化图。根据图表总结分析，我们得出A公司在电导率方面存在严重问题，B在不断完善中，D公司在各方面危害指标的管理存在严重问题。针对问题四，我们运用了规划模型，运用分层求解法求得了公司的规模划分问题，再将将A公司的各点的危害因素进行加权平均，然后进行无量纲化的处理，运用模型一求出各存放点的安全风险度值进而得出A公司各存放点批次的分布。针对问题五，我们根据以上建立的模型分析出了此城区饮用水的在各项危害指标中都存在隐患，并提出了5条相关性的监控意见。【关键词】纯净水安全风险评价安全监督线性综合评价模型变异系数法层次分析法综合权向量法规划安全形势监控措施一、问题的重述日趋加剧的水污染，已对人类的生存安全构成重大威胁，成为人类健康、经济和社会可持续发展的重大障碍。据世界权威机构调查，在发展中国家，各类疾病有8%是由于饮用了不卫生的水而传播的，每年因饮用不卫生水至少造成全球2000万人死亡，因此，水污染被称作"世界头号杀手"。我国政府对纯净水安全问题十分重视，已将纯净水安全作为一项重要的公共管理目标，采取了一系列措施，强化纯净水安全的监管，并取得了初步成效。但纯净水安全问题的总体形势仍不容乐观，依然存在一系列隐忧，近年来食品安全方面的恶性、突发性事件屡屡发生。2007年07月12日，南通一纯净水厂发生造假事件。2008年3月底，贵阳市发生数百人感染甲肝事件，经卫生部中国疾控中心专家组核查，确认“竹源牌”桶装水是造成疫情爆发的主因。2009年03月25日，某大学B区学生饮用了“清清”牌桶装纯净水后，百余学生先后出现集体腹泻事件。2009年2月26日，湖南师范某寝室在长沙爱高普纯净水有限公司订购的桶装纯净水中出现了黑色虫子事件。生物性和化学性污染对纯净水安全的影响愈来愈严重。本问题主要考虑纯净水的以下危害因素:（按照危害的严重性依次给出）“电导率”:是纯净水的特征性指标，反映的是纯净水的纯净程度，以及生产工艺控制的好坏，“电导率”根本达不到国家卫生标准要求，与自来水无异，根本不能算做纯净水。菌落总数:是指纯净水检样经过处理，在一定条件下培养后所取1ml（g）检样中所含菌落的总数。它可以作为判定纯净水被污染程度的指标之一。大肠菌群：反映纯净水加工过程中对大便污染程度的一个指标。数值越高证明污染越严重。霉菌：食物霉变后产生，直接引起中毒，或产生致癌物质，毒害人体。纯净水的安全危机的爆发，往往是日常的监控机制和管理长期存在漏洞的反映。完整、有效的纯净水安全风险分析监测预控，为政府及有关部门实施控制措施提供决策依据和技术支持，可以有效提高纯净水安全监管效率和管理水平，及时化解可能出现的安全危机。近年来，我国在从国家宏观层面探讨建立纯净水安全预警机制的研究方面，已取得了不少理论成果但由于我国地域辽阔，经济社会发展水平很不平衡，如何构建有效的预警机制并应用到饮用水安全监控过程还处于起步阶段。某城区共有九家生产并销售纯净水的公司，其中A公司和B公司规模较大，其余均为小公司。针对该城区提供的近年的关于各公司的纯净水检测报告（见附件），请你利用数学建模的方法回答以下问题：1、结合本问题所给数据，给出纯净水安全风险分析的科学评价方法，确定评价的标准和评价的规则，对该城区所有批次的纯净水进行评判排序。2、对该城区范围内的监控对象（各公司）按风险度进行排序评价，并对它们分类综合评价，指出各公司产品的主要可能的危害因素，并指出同类公司的实际特点。3.对检测出的不合格的样品成因分析：评价纯净水生产流通环节（归为仓库和销售网点两类）的危害因素以及各个危害指标的分布规律，并通过四类危害指标的分析，讨论A、B、D公司的管理状况。4.国家相关部门每年要面对各种专项检验，对于纯净水专项检验的投入经费有限，已知该城区下一年度投在纯净水方面的检验总批次为100个批次，在现有历史数据的基础上，并考虑各公司的实际运行状况，如何设置各公司检验批次的分布，使得抽检方案的针对性最优（即检出的风险性为最大）。5.结合你的工作，请你给该城区食品安全委员会写一篇短文，阐述你的观点，评价该城区的饮用水安全形势并给出监控对策。二、模型假设(1)各个样品有一个危害因素不合格则视为该样品不合格；(2)采样的地点、样品和数量均是随机的；(3)各公司的规模与被检验的总批次数量成正比；(4)假设题目所给数据真实有效；(5)本文只考虑所检测的四个危害因素，其他因素忽略；(6)各个公司在一年内的生产模式基本保持不变；(7)各个公司放在仓库和销售点的产品相同；(8)本批次检验的产品均是待出售的现货或存货；(9)在一年内没有被检查，则视为该批次在这一年为合格。三、符号说明：在第k时间时第j批次的样品经项目i检验的结果；：在第k时间时第j批次的样品经项目i检验的样品数量；：第j批次的样品经项目i检验的结果；：数据经过无量纲化处理后对应的数值；：i因素的权重；：权重向量；:i地点抽检了n批次；:纯净水的安全风险度；:市场占有率:各公司的风险度:各公司的批次数四、模型的分析与准备文中已给出纯净水检测的数据，对此数据进行分析。考虑到以产品的批次来分类，所以需对原始数据进行简单处理。其中，对霉菌和酵母一项“没有检验出”时看作0来处理。按批次进行平均加权：由于各组数据的单位不同，需对其进行无量纲化处理：得到yij。相关数据见下表：Tab.1各批次不同因素的加权平均加权平均后数据无量钢化后数据批号电导率（≤10）菌落总数（≤20）大肠菌群（≤3）霉菌和酵母（≤0）电导率（≤0.148）菌落总数（≤0.0625）大肠菌群（≤1.5）霉菌和酵母菌（≤0）200704153.910100.05809800.50200704252.61000.0386330.003125002007051213000.0148590.00937500200706091.725200.0255570.01562510200707037.2335100.1074290.1093750.50200707045.433200.50.250.08068410.250.25200708094.0050.5100.059510.0015630.50200708104.7320.500.0702820.006250.250200712163.71.670.33300.0549780.0052190.16650200712184.140.5000.0615160.001563002007121927.21000.404160.00312500200712203.370000.050074000200808264.560000.067756000200808278.245160.250.250.1225110.050.1250.25200809013.7534000.0557210.10625002008090221.345000.3164930.14062500200901091.820000.0270430002009011027.80000.4130760002009011167.37.50110.023438012009011225.6860.400.3815750.018750.20五、模型的建立与求解5.1模型一的建立与求解（问题一）本模型我们考虑最简单的线性综合模型：我们在原题所给数据的基础上，利用变异系数法计算各项指标的权重:其中：是第项指标的变异系数、也称为标准差系数；是第项指标的标准差；是第项指标的平均数。各项指标的权重为：由题目所给数据，利用以上公式和无量纲化后的数据计算各指标的权重系数如下表：Tab.2各指标的权重系数电导率菌落总数大肠杆菌霉菌和酵母菌平均值14.82868.15500.17740.0893标准差0.23660.22190.26640.2308变异系数1.38752.96941.52623.0779权重系数0.15480.33140.17030.3435，此权向量完全是由数据计算而得到，客观性较强。但是，由于确定各项指标的权重需要综合考虑评价对象的属性机理等因素，因而不能仅仅依赖数学计算处理后的数据，必须结合各项因素的社会属性等，给出另一个主观性较强的权重向量。这一点题目中给我们提供了各项指标的重要性程度，因此我们可以利用层次分析法，根据“1-9尺度”原则，建立成对比较矩阵：经一致性检验，CR=CI/RI=0.0433<0.1（程序见附录road1.m）,所以A的特征向量可用作权向量。从而得到另一个权重向量：(程序见附录road2.m)。此权向量由生活中的经验常理等得出，主观性较强，与𝑊1可以互补。最后综合考虑这两方面我们可以得到权重向量𝑊。利用综合法，计算𝑊：𝑊=𝑎𝑊1+（1−𝑎）𝑊2其中0≤a≤1，为调整主客观因素比例最终决定各目标权重系数的因子，它可以根据不同系统和实际需要而调整。此处，我们取a=0.6。得到𝑊=(程序见附录road3.m)。将𝑊代入（3）式可得模型一：Q表示安全风险度，Q值越大，此纯净水的质量越不安全。由此，我们可以对各个批次的纯净水进行评判排序。数据如下表：Tab.3所有批次的纯净水安全风险评判排序名次批号Q值名次批号Q值1200705120.007311200708090.09922200901090.008212200712190.1243200704250.012713200901100.12584200712200.015214200808270.13085200712180.019215200809020.13876200808260.020616200707030.14637200712160.045217200901120.15418200809010.04918200706090.17379200708100.063619200707040.424310200704150.098320200901110.5446从上表可以看出，批次20070512的纯净水的安全风险最小，最安全；批次20090111的纯净水安全风险最大，最不安全。进一步分析：结合表一可以看出，批次20080902的水在电导率、霉菌和酵母菌两项检验上严重超标，水质量严重不合格。批次20090110的水在菌落总数一项上超标及其严重，水质量也严重不合格。考虑到每一批次中，有一项检验不合格，则此水质量即为不合格。以上被检验的所有批次中合格的水有11批20070415，20070425，20070512，20070609，20090109，20070810，20071216，20071218，20071220，20080826，20070809。其他9个批次均不合格。5.2模型二的建立与求解（问题二）5.2.1公司的安全风险综合评价Tab.4各公司不同因素的加权平均各公司不同因素的加权均值无量纲化后公司电导率(≤10)菌落总数（≤20）大肠菌群（≤3）霉菌和酵母（≤0）电导率（≤）菌落总数大肠菌群霉菌和酵母A14.932531.650.250.050.71450.20870.19020.1305B3.00875151.68750.37500.14401.00000.32810.0000C2.09713.28571.142900.10030.02171.00000.0000D18.1333398.16660.50.33330.86760.64720.43751.0000E2.61000.12440.00660.00000.0000F3.910100.18710.00000.87500.0000G13000.04780.01980.00000.0000H16.79289.80.200.80340.59200.17500.0000I20.930.66670.333301.00000.20220.29160.0000将量纲化后的数值代入（5）式可得各个公司的值(程序见附录road4.m)，以此来排名：Tab.5所有公司的纯净水安全风险评判排序排名公司Q值1G0.02052E0.03993F0.19814C0.19835A0.34156B0.39797I0.41248H0.45119D0.7627由上表可以看出，公司G的纯净水安全风险最小；公司D的纯净水安全风险最大，最不安全。统计各个公司各批次在抽查中有不合格的项目次数，得出：Tab.6各个公司各批次不合格项目次数统计不合格次数电导率菌落总数大肠杆菌霉菌和酵母A4301B0100C0000D1202E0000F0000G0000H2200I2100分析此表可得：C、E、F、G公司纯净水质量最好，不合格项目为0记录。公司A，D不合格记录最多，其次是公司H、I，B较少。5.2.2公司分类的评价同样再对各个公司分类进行评价。根据表三的数据，我们分别按这四个危害因素将这九个公司进行排序评价：Tab.7各公司分类分析评价表安全风险度电导率菌落总数大肠杆菌霉菌和酵母1IBCD2DDFA3HHDB（0）4AABC（0）5FIIE（0）6BCAF（0）7EGHG（0）8CEE（0）H（0）9GF（0）G（0）I（0）上表中，排名靠前的公司在此项目上的安全风险度越大。从表中可以看出，各个公司在各项的评价排名上有很大区别。其中显而易见，A、D两家公司多次排名靠前，产品安全风险较大，不安全。而E、F、G三家公司排名较后，产品较安全。对表格七的数据进一步分析：可以找出各公司产品的主要可能的危害因素，如下表：Tab.8各个公司产品主要可能的危害因素公司主要的危害因素A霉菌和酵母B菌落总数C大肠杆菌D霉菌和酵母E电导率F大肠杆菌G菌落总数H电导率+菌落总数I电导率由上可以看出，有些公司有多种危害因素，但只按照其危害程度相对其他公司较大的来计。有些公司虽然没有检验出不合格的产品，但是它的某项指标的数值较大，比较接近标准值，是潜在的危害。5.2.3公司的实际特点分析结合表八，做进一步的分析可以找出同类公司的实际特点。利用文中给出的对电导率菌落总数大肠杆菌及霉菌和酵母菌的定义和解释，结合我们查找的有关饮用水安全等方面的知识，对这九家公司的实际特点分别作出说明。结论见下表Tab.9公司的实际特点公司实际特点A电导率菌落总数霉菌和酵母都超出了标准值，且霉菌和酵母是致癌物质，纯净水纯净程度较低，生产工艺控制的不好，污染程度较高，可能会致癌，所以公司生产的纯净水有比较严重的质量问题。B菌落总数极其高，远远超过了标准值，说明该纯净水被污染程度高，该公司的水源地已遭到破坏或者在生产过程中生产工艺控制的不好。C大肠菌群的含量较高，但还没超标，说明该公司的水源受到了大便的污染，质量稍有问题。D霉菌和酵母量很高，电导率和菌落总数也均超标，说明该纯净水有非常严重的质量问题，而且该公司的水风险度最大，对人体很有害，购买时应慎重考虑。E各危害因素都很低，属于比较安全的纯净水。F大肠杆菌的含量相对于其他公司稍高，但没有超标，水源受到一定的污染，其他方面较好，没有什么太大问题。G各危害因素都很低，该公司的水风险度最小，属于比较安全的纯净水。H电导率和菌落总数均超标，该公司的水的风险度较大，说明该纯净水有较严重的质量问题，对人体害处较大，购买时应慎重考虑。I电导率和菌落总数均超标，大肠杆菌的含量也相对较高，且公司的水的风险度最大，说明该纯净水有较严重的质量问题，对人体很有害。5.3模型三的建立与求解（问题三）5.3.1生产环节的评价以样品为单位，进行不合格产品的成因分析。数据中共有9个公司，20个批次和70个样本被检验。按照纯净水生产流通环节来分可分为仓库和销售网点两类。对其不合格样品数量进行统计分析，可得以下数据：Tab.10不合格率分类统计表仓库销售点不合格样品数不合格率不合格样品数不合格率电导率1535.71%13.57%菌落总数2047.62%00.00%大肠杆菌00.00%00.00%霉菌和酵母511.90%00.00%不合格样品总数3276.19%13.57%样品总数4228由上表可以清楚地看出：位于仓库的水不合格率高达76.19%；相比之下，销售点的产品质量要好些，不合格率只有3.57%；可见销售点水的质量问题比较严重，可能是时间过长。由上表数据将两者进行对比做出下表：Fig.1仓库与销售店的不合格率的比较由上图可以看出：销售网点的水少部门在电导率方面存在问题，其他方面还好；而仓库的纯净水在电导率和菌落总数方面存在着很严重的问题，不合格率很高，而且在霉菌和酵母这一项上有较大问题；不过在大肠杆菌上两个都还好，没有不合格的问题。所以可知，导电率偏高是各个地方普遍存在的问题。5.3.2危害指标的分析我们把各个不合格因素进行分类，按照检验的批次顺序，将各个危害指标中不合格的公司进行汇总。得出的结果如下表：Tab.11各批次不合格危害因素统计表危害因素不合格公司电导率HAAIDHIAA菌落总数BDAHIADHA大肠杆菌无霉菌和酵母DAD由上表可以的到各个危害因素在各个公司的不合格批次产品中的数量分布，但由不合格次数不足以判断出完整的各个危害因素的分布规律，还需结合不合格的程度进行综合衡量才行。对于表一中的数据，将无量纲化后数据绘制成图标，得到四个危害指标随着检验批次及公司的变化图。图中的横坐标是以批次的顺序即以时间的变化的。所以此图可以看出四个危害指标随时间的变化的大体走向。Tab.12各公司检验次数ABCDEFGHI11112411132Fig.2检验指标变化趋势图在上图中，有四种不同颜色的实折线和虚直线，当实折线超过相同颜色的虚直线代表超标，检验不合格。各点即为检测数据。结合表九可以得出以下规律，针对四个危害因素分别给出：（1）电导率的分布规律结合表九，从图中的蓝线可以看出，从不合格样本数量来看，电导率偏高的纯净水样本主要分布在A公司，其次是H和I，D公司最少。其他公司在此项没有检验出不合格。从具体数值来看，D厂只抽检了一批电导率不合格的纯净水，但数值高达84.4，若认为抽到电脑率极高的纯净水是小概率事件，在不考虑特殊原因的情况下，恰恰在D公司发生了，所以可以认为D公司此批次大多数的纯净水电导率都没控制好，都有严重的问题。而D公司其他批次的水没有检测出问题，说明D公司在总体上来看，在此项上问题不是很大。结论：电导率这项危害主要分布在A公司，其次是H、I。2007年此危害较少存在，但自2008年此项危害出现率较高。（2）菌落总数的分布规律结合表九，从图中的红线可以看出菌落总数不合格的情况在A、B、D、H、I公司都存在，而A、B、D、H公司都有严重超标现象。其中，同（1）中的D公司类似，B公司只检验出一批菌落总数超标的纯净水，但该值高达800，说明B公司在此批次的纯净水在菌落总数上没有控制好，导致大都出现非常严重的问题，而B公司总体上在此项没太大问题。结论：菌落总数这项危害因素主要分布在A、H、D公司，其次是I。此危害在2007年有一个极大值，但自2008年来得到较好的控制。（3）大肠杆菌的分布规律结合表九，从图中可以看出，检验中并没有发现此项超标的情况，说明总体上纯净水在此方面没有太大问题。但是除了E、C、G非常干净外，其他的所有是都存在不同程度的轻微污染。此危害指标程度不定，随时间的波动性较大。（4）霉菌和酵母的分布规律结合表九，从图中的紫线可以看出，霉菌和酵母超标只出现在A、D两公司，由由上表数据可以看出该危害因素主要分布在D公司，其次是A。此危害指标随时间出现间歇性超标现象5.3.3A由上面的分析A、B、D厂生产的纯净水分别在电导率、菌落总数及霉菌和酵母这三项危害因素上有严重问题。A公司生产的纯净水在三年的纯净程度均不是很理想，纯净水受到不同程度的污染，该企业的管理状况不容乐观，但在2009年生产的纯净水受污染程度明显减少，说明管理得到了一定的改善。但在导电率方面非常差，生产工艺控制方面的缺陷依旧严重，急需改善。B公司生产的纯净水于2007在生产流程中受到了一定的污染，但2008年以来菌落总数已得到很好的控制，其他方面该公司做得很好，所以企业的管理在不断的完善之中。D公司生产的纯净水在电导率、菌落总数、霉菌和酵母这三项指标上都有严重的问题，而且3年以来这些问题一直存在，说明该企业在管理方面有很大的缺陷，生产工艺控制方面做得很糟，其产品对人的健康是一较大的安全隐患。5.4模型四的建立与求解（问题四）5.4.1按公司规模划分此问是典型的规划问题，目的就是使各批次检出的总风险最大，即抽检出的各个公司的水总风险性最大，我们这里运用分层求解的办法求解。首先把100批次分配到各个公司。题目说9个公司中A公司和B公司规模较大，其余均是小公司，所以在分配抽检批次时，必须考虑此项因素，这里我们做一下粗略的划分，参照前面对各公司的抽检次数将各个公司规模划分如下：在第二问中我们依旧求解出各个公司的风险度（见第二问求解），分配时综合考虑市场占有率和风险度两因素：Tab.13公司市场占有率统计表公司ＡＢＣＤＥＦＧＨＩ市场占有率28.57%31.43%5.71%11.43%2.86%2.86%2.86%8.57%5.71%设市场占有率为，风险度为(数据见表四)，则=*,因此当批次总数为100时，分配到各公司的批次数为:=100结果见下表：Tab.14各公司检验批次的分布公司ＡＢＣＤＥＦＧＨＩ批次分配25323220101065.4.2对公司批次进行划分首先我们根据数据表发现各公司的不同时期风险度是动态变化的，如图所示：Fig.3A公司风险度随着时间的变化先增加后减少，说明该公司在风险度低时会疏忽对生产工艺等方面的监督，而使风险度变大，当风险度逐步增大时因为城区监管部门的检验有逐步减少到极小值，如此周期循环下去，符合我们的一般经验。这里取它的极大值作为它的风险度最大值，极小值为它的风险度最小值。我们根据题目所给数据，将A公司的各点的危害因素进行加权平均，然后进行无量纲化的处理，如下：Tab.15公司各点加权平均值无量纲化后的数据A公司纯净水存放点电导率（<=10)菌落总数(<=20)大肠菌群(<=3)霉菌和酵母(<=0)电导率(<=0.4732)菌落总数(<=0.125)大肠菌群(<=3)霉菌和酵母(<=0)仓库21.13084.00000.16670.08331.00000.02500.16670.0833厂成品库7.230035.00001.00000.00000.34220.21881.00000.0000成品库5.3500160.00000.00000.00000.25321.00000.00000.0000店内3.97000.00000.00000.00000.18790.00000.00000.0000货架3.37000.00000.00000.00000.15950.00000.00000.0000货架上0.00000.00000.00000.00000.00000.00000.00000.0000营业部0.00000.00000.00000.00000.00000.00000.00000.0000再根据安全风险度公式求得各点的安全风险度值，具体数据如下表：Tab.16AA公司纯净水存放点电导率(<=0.4732)菌落总数(<=0.125)大肠菌群(<=3)霉菌和酵母(<=0)Q值仓库1.00000.02500.16670.08330.3583厂成品库0.34220.21881.00000.00000.3313成品库0.25321.00000.00000.00000.3782店内0.18790.00000.00000.00000.0572货架0.15950.00000.00000.00000.0486货架上0.00000.00000.00000.00000营业部内0.00000.00000.00000.00000将上表中计算得到的Q值，进行归一化处理，：从而得到A公司各存放点的市场占有率，进而得到A公司各点的批次分布情况，如下：Tab.17AQ值各点占有率各点的检验批次0.358330.53%80.331328.23%70.378232.23%80.05724.87%10.04864.14%100.00%000.00%0同理，我们根据以上方法，分别求得B、C、D、E、F、G、H、I公司各存放点的批次的分布，如下表：Tab.18各公司各点风险度值最大的批次分配方案公司仓库厂成品库成品库店内货架货架上营业部内总数A878110025B6160700332C00003003D7015000022E00000101F00001001G00000000H503200010I600000065.5给食品安全委员会的一封信尊敬的食品委员会的一封信：您好！最近我们对贵城区的纯净水的安全形势做了研究分析，对纯净水的安全风险评判建立了数学模型，利用此模型对纯净水安全监控问题提出一些建议，并对城区九家纯净水公司的产品卫生安全情况作了分析评价。针对该城区近三年来九个纯净水公司的产品质量检测报告，我们经过建模分析发现纯净水中危害因素电导率超标现象越来越普遍，问题很严重；菌落总数出现间歇性超标，但其总体有下降趋势；菌落和霉菌控制方面也一直控制得不好，问题一直存在。同时我们对城区的九家纯净水公司做了安全风险评判，其结果为，G公司产品最安全，其排名依次为。按照目前趋势如果不采取措施，贵城区的饮用水安全环境会逐渐降低，居民的用水安全将受到严重威胁。为了解决此问题我们小组提出几条纯净水安全监控对策：1、缩短对于葛根素生产纯净水的检测周期，检验各公司生产纯净水的水源的各项指标，然后找到各公司的相对薄弱环节并进行科学的指导与改进，并加强薄弱指标的检查力度；2、在生物性和化学性污染日趋严重的大环境下，提倡各公司强化传统生产水的工艺、替换传统的消毒剂、吸附、膜过滤和生物预处理等深度净水技术3、为了提高居民对自身饮用水的安全意识，对每期纯净水抽样检验结果进行公布，并将各公司进行排名，促进各公司之间的竞争，带动纯净水行业的进步，加大排名靠后的公司的检查力度；4、加强各公司对人员培训和生产过程的监督与管理。安全优质的引用水不仅取决于优质的水源、优质的净水工艺以及安全的供水网络，也取决于制水人员的卫生意识以及对水质变化的应急处理能力。因此需加强水厂制水人员、管理人员及监督人员的专业训练；5、对于多次检查不合格的且不接受科学技术的公司依照法律进行停产整治，保障人民饮水安全的保障上面是我们小组在数学建模中获得的结果与想法，希望能给您们带来工作上的帮助。六、模型的评价与改进6.1模型的评价优点：用合理的方式处理不同类型的数据，使其尽量符合实际情况在权重系数的确定中，运用变异系数和层次分析法相结合的形式，变异系数法仅仅用数据求权重，客观性较大；而层次分析法利用了各个指标的重要性程度，主观性强，两者结合使问题的考虑更全面更符合实际在问题求解中运用了大量的表格和图，使结果明了、清晰在求解最有分配方案时，我们利用了层次分析，巧妙的把复杂的问题转化为单个的线性的规划问题，简化了计算过程。缺点：对每个批号的检验结果取平均值，精确度不高在确定建模一的权重时，虽然运用了综合权重法，但尚有不足之处不是很精确6.1.2模型的改进对于我们给出的权重模型，在考虑权重时，我们可以再利用几种权重分析法，例如主成分分析法，负相关系数法等，最后再综合各权重向量得到一个比较贴近实际的权重分配，这样可以增加模型的可信度。对于纯净水安全问题，它会随着人的主观意志而改变，不会是静态的，也不是很确定，单一模型很难满足实际情况，故我们可以建立一个动态的或多种机制综合来衡量。七、参考文献[1]姜启源,谢金星,叶俊,数学模型（第三版）.高等教育出版社2003.8[2]赵静,但琦,数学建模与数学实验（第三版）.高等教育出版社2008.1[3]陈杰,MATLAB宝典.电子工业出版社2007.1八、附录程序road1.m%层次分析法求权重clear;clca=[13571/31350.21/3131/70.21/31];[x,y]=eig(a)eigenvalue=diag(y);lamda=eigenvalue(1)ci1=(lamda-4)/3;cr1=ci1/0.9w1=x(:,1)/sum(x(:,1))程序road2.m%迭代求权重clc;clear;A=[13571/31350.51/3131/70.51/31];[V,B]=eig(A)e0=[0.250.250.250.25]';b=A*e0;x=b(1,:)+b(2,:)+b(3,:)+b(4,:);e1=b/xm=A*e1;x=0;forn=1:4x=x+m(n,:);ende2=m/xm=A*e2;x=0;forn=1:4x=x+m(n,:);ende3=m/xm=A*e3;x=0;forn=1:4x=x+m(n,:);ende4=m/xm=A*e4;x=0;forn=1:4x=x+m(n,:);ende5=m/xm=A*e5;x=0;forn=1:4x=x+m(n,:);ende6=m/xm=A*e6;x=0;forn=1:4x=x+m(n,:);ende7=m/xm=A*e7;x=0;forn=1:4x=x+m(n,:);ende8=m/x程序road3.m%安全风险度计算clear;clc;W1=[0.15480.33140.17030.3435]';W2=[0.52890.25580.14770.0676]';w=0.6*W1+0.4*W2y=[0.05809800.500.0386330.003125000.0148590.009375000.0255570.015625100.1074290.1093750.500.08068410.250.250.059510.00156250.500.0702820.006250.2500.0549780.005218750.166500.0615160.0015625000.404160.003125000.0500740000.0677560000.1225110.050.1250.250.0557210.10625000.3164930.140625000.0270430000.41307600010.0234375010.3815750.018750.20];Q=y*w程序road4.m%公司安全风险值clear;clc;W1=[0.15480.33140.17030.3435]';W2=[0.52890.25580.14770.0676]';w=0.6*W1+0.4*W2y=[0.71450.20870.19020.13050.14401.00000.32810.00000.10030.02171.00000.00000.86760.64720.43751.00000.12440.00660.00000.00000.18710.00000.87500.00000.04780.01980.00000.00000.80340.59200.17500.00001.00000.20220.29160.0000];Q=y*w

附录资料：不需要的可以在删除安全生产设备管理方案第一章生产设备的安全机器是现代生产的核心。然而，随着生产率的提高，机器又给生产现场带来危险。传统认为，对机器危险做适当控制既花钱又影响生产效率。另外，工人在操作机器或在机器周围工作时可能拆掉防护装置或拒绝使用个人防护用具，再有，许多企业根本不重视机器的安全。本章将阐述对机器安全方面的许多消极态度的不合理性。消除机器危险，同时增加产量完全是有可能的。即便是要安装防护装置，也不一定就意味着花钱多或降低生产率。在某些特别的位置确实存在着机器危险情况，比如：作业点、机器动力传输点或移动加工件周围。本章重点阐述作业点的危险，因为这个位置的防护措施通常比较困难，但又与生产效率有密切关系。其他危险当然也不能忽视。一、对机器进行生产效率检查：在工厂里转一转，观察下面的情况：有没有某个机器造成生产的迟缓或阻塞有没有由于填料或退料装置造成机器运作缓慢有没有由于机器或生产过程的危险造成害怕或踌躇的现象有没有机器防护装置被改动、拆掉或毁坏的现象如果上述任何一个问题的回答是“是”，那么下面的指导则有助于处理机器生产率和安全的问题。消除危险，安装保护装置，使用个人保护用具，一定要按这个顺序做：发生事故不仅会带来经济损失，还会给人带来痛苦。然而，在存在许多危险因素的车间做到避免事故发生并非易事。这需要有一个良好的策略。其中最好的办法莫过于完全消除危险。这是否听起来不大可能？本文有许多彻底消除危险又能提高生产效率的办法。如果不能消除危险，那么就在周围安装防护装置。但防护装置必须要认真设计，否则将会起妨碍作用。把个人保护用具的使用作为最后办法，如果必须要使用个人保护用具，那么就要保证时常监督用具的正确使用：用危险小的机器或生产过程取代危险机器在危险区周围设立防护设施在危险未消除或未设防护设施之前，提供个人保护用具购买安全的机器：当订购新机器时，要特别注意强调机器结构的安全性。危险加工件应该处于不会伤害工人的位置，特别是操作点必须没有危险。自动化和机械化的填料和退料装置不仅能够消除危险，而且会大大提高生产效率。要考虑订购带有这样装置的机器，因为这些装置相对整个机器的成本来讲并不昂贵，但有与没有在提高生产效率方面大不一样。采用填料和退料装置提高生产率，减少机器的危险：许多机器（特别是冲压机）的生产效率受限于加工部件填入机器和从作业点退料的速度。在不采用填料和退料装置的情况下，生产速度将受工人操作的速度、工作的复杂性、防护装置的种类和必要的安全措施的影响。如果没有自动化和机械化的填料和退料装置，工人就会摆放冲压件，然后让手躲开机器的循环运转，之后把冲压件拿开，最后从机器上取出冲压成品。实际上，只有20%的时间用于机器生产。这还没有考虑在不是自动化的情况下作业点防护装置的拆卸和替换。有填料和退料装置，生产效率就会大大提高。然而，填料和退料方式也要考虑安全措施。在进行设备安装的过程中，在安装和调试填料和退料系统时，在清除填料故障和做润滑和保养的过程中，都要极度小心。另外对新的填料和退料装置进行评估时，也要非常仔细认真。有很多不同种类的填料和退料装置，如下：手工填料的基本办法就是不让操作工人的手进入危险区域这种填料装置最简单的形式就是推滑板（如图3-1）。推滑板上有一个模孔，工人在冲压点外部把冲压件放在模槽内。当推滑板推进冲压点时，机器就做冲压循环。转盘填料装置（3-2）根据同样原理，可以同样大大提高冲压操作工人的生产效率。如图，填料板有一个自动退料和收集冲压成品的装置。自由落体填料桶（3-3）可以结合到滑板装置上，这样省去了每一工作循环放入新的被冲压件。显而异见，采用滑板加装填料比手工填塞被冲压件要节省很多时间。还有许多从“自由”重力获益的其他办法。如图（3-4），简单的倾斜坡面填料装置可以用来把被冲压件送到默子里。注意：图中防护网把工人与冲压点隔开。图3-5，这个容易制作的重力填料装置可以使研磨作业的生产效率提高一倍。填料和退料系统的形式多种多样，从非常简单的到非常复杂的。有些容易设计和安装的是采用半自动化，以及采用斜面、转盘和滑板。退料作业可以采用机械装置、压缩气体或是填料系统整体的一部分装置。采用正确的防护种类：机器的运动形式不同，生产要求也不一样。同样，也存在不同种类的危险。转轴、轮子、滚轮、滑轮及齿轮会卷衣服和皮肤，也会毫不含糊的把工人拉进机器里。如果一部分与另一部分相对旋转（3-6），就会产生一个“夹点”，而这个夹点就可能抓住工人的手或衣服。机器来回移动，往复运动、移动的机器或毛坯与固定物体之间部分就可能抓住工人。直接进行切割、冲压、裁剪和板金操作也会导致严重事故发生。因此，根据机器要求或针对具体的危险特点配置相应的防护措施。1机器防护的种类：固定的防护装置比较简单，它直接固定在机器上或一稳固面上，比如墙壁或地面。制作它们的材料应该结实，可以起抵御飞削的作用。作业点固定防护装置应该与填料与退料系统相结合，从而不妨碍生产。固定防护装置应该只能用工具拆卸。互锁定防护装置要不时的与固定防护装置和防护盖结合使用。当防护装置或防护盖打开或移走时，这种装置就断开机器的供电或机器的工作（3-7）。它们还可以在机器工作周期之前挡住进入作业点（3-8），但必须特别当心的是：当作业过程有惯性时，要等机器停下以后再打开护盖。可调式防护装置可以调整同送入作业点的加工件尺寸相吻合，同时又有很好的保护作用（3-9、3-10）。双手控制。事故发生的原因往往是工人在用一只手填料加工件的同时，另一只手在操作开关，如果双手在时间上配合不当，那么当一只手还在机器里时，机器就开始工作了。一个可能的解决办法时，把控制机关设计成两个开关或搬动杆必须用两只手同时操作才有效。这样，工人的手必须移出机器后机器才能工作。但是，设计不用一只手操作的控制按钮或机关并不容易，当有粘连或卡住的情况时，要用漆盖顶住或采取防护机器继续运转的措施。因此，非常有必要监督这些装置使用是否使用恰当。2设计机器防护装置：制作防护装置的前提条件通常是低成本，有助于设计和制作防护装置的要点如下：完全封闭的防护要优于部分封闭的防护。应该避免部分防护防护装置尽可能的装在危险区域。但要当心：移动的皮带同固定的防护装置产生摩擦不仅会损坏皮带，而且还可能产生火花的潜在危险防护装置在可能的情况下要与填料和退料系统相结合为了不降低生产效率，防护装置在做简单维修时必须能快速打开和关上设计用来断开作业的合叶式防护如果打开进行维修的话，建议在重新作业之前更换新的防护装置应该保证作业区清晰可见任何妨碍生产效率的防护装置应该设计适当的维修机器：即便装了防护装置，机器维护不良同样会造成危险。他还会产生较多的停机和质量问题。适当的维护并非损失生产时间，而是一种提高生产率、降低修理费用的投入。维护应该包括防护装置。工人通常把护板拆下来清理机器，更换模具，或者进行维修加油。再有，防护板还常常由于降低了工人监督作业区产品质量的视线而被拆卸下来。这些防护装置应该检查、清洗、甚至更换，如果必要，应该重新设计。由合格人员对机器进行维护保养将会减少修理的频率，减少工人拆卸防护装置的次数。维护项目还包括必要监督质量的视线部分进行日常清洗。当对机器进行修理或保养时，机器的控制开关必须锁上，并贴上标签说明“危险勿动！”如果没有其他保护措施，采用个人保护用具：如果不能消除危险，也没有适当的防护装置保护工人，那么最后的办法就是采用个人防护用具。但有时危险在于个人防护用具会日益磨损老化。如果必须使用个人保护设备，那么只能根据具体危险情况采用国家认定的符合标准的用具。个人防护用具也必须根据根据的说明进行维护保养，当损坏或破旧时必须更换。使机器更安全，更有效的规则对机器进行效率检查消除危险或安装防护或在别无办法的情况下使用个人保护用具，一定要按这个顺序进行购买安全的机器采用填料或退料装置以提高生产率，减少机器的危险使用正确种类的防护装置适当的维修机器如果没有其他防护措施，使用个人保护用具第二章危险品的控制在中小企业里随处可见各种各样的危险品。对于大多数危险问题进行控制的方法既简单又廉价，改进后既可节约成本又可提高生产率。污染的环境往往对生产不利。大量的尘埃、油料、油漆和其他涂料都会影响生产效率、增加监测和清洁工作，还会污染原料和产品。因此，控制危险品对提高生产率和产品质量有很大潜力。大量化学物质的暴露会导致疲倦、头痛、头晕，并会对眼睛和呼吸道产生刺激，从而引起生产率和产品质量以及工人出勤率的下降。如果条件得到改善，劳动生产率就会大大提高。许多化学品危害问题不用花很多钱便可得到解决，同时还可以减少原材料和能源的消耗。用危险小的物品代替危险大的物品：很多企业用有机溶剂清晰金属和其他材料。有机溶剂既危险又昂贵。可以用危险小的物品取代有机溶剂，比如清洁油污一般用皂类清洁剂。含量为5-10%的肥皂水比有机溶剂既便宜又安全（4-1）。用碱性物质如氢氧化钠或氢氧化钙可以更彻底的去除油污。用这种办法清洗金属表面一般比较彻底，再上油漆不成问题。为防止生锈，洗净后的每一个加工件要干燥后分别保存。用类和碱类物质进行清洁和出垢比用有机溶剂更安全、更节省，而且不需要更多通风设备和氧化。但在碱性物质时要用防护眼睛或手套。用盖子、容器和隔离措施控制危害，减少损失：物质从开口的容器里蒸发或泄漏的容器或管道流失会导致直接损失、增加成本，且危害工人的身体。有机溶剂、油漆或胶的烟雾可用特殊盖子防止有害物质和废物的蒸发（4-2）。加盖隔离可减少机器如车床、刨床、磨床或搅拌器产生的污染。机器保养不善也会导致不必要的危险和损失。车床或其他机器的漏油滴在工人的皮肤和衣服上可导致湿疹和皮癌。这样的泄漏也浪费昂贵的机油。把某些工序限制在隔离的房间或空间可有效地减少危险。对少数接触污染物的工人应加以特别保护。节约用于加热化学物品的能源：高温液体蒸发可能会在空气中形成密集的危险物质。在装有有机溶剂、碱和酸的容器加热清洗或去除油垢时，胶质或带有涂层材料熔化时，在进行电镀时都有可能出现这样的问题。做这些工作往往不了解最高温度限制。温度过高增加了工人接触危险化学物品的机会，同时意味着在蒸发过程浪费更多的能源。使用恒温器可以获得最低的适当温度（4-3），如果没有恒温器，可用电流表加人工温度控制。正确地进行清洁，不要散落尘削：尘削来自磨、锯、钻、包装、搅拌和其他生产过程。尘削颗粒的大小及其危险程度各不相同。尘削增加对机器的磨损，从而增加对保养的要求，同时还可能影响原材料和成品的质量。尘削进入呼吸系统还会对工人的肺部造成危害。有些尘削还会通过皮肤吸收而导致疾病。尘削必须经常打扫。多数尘削须在产生时用排风装置彻底清除（如在磨床和圆形锯床接口处）。残留的尘削必须每天清除，除此之外还须定期对墙壁、储藏架和其他容易积尘的地方进行大扫除。玻璃窗、墙壁和台灯上的灰尘会很大程度的影响照明。注意：不要扫地扬尘。用压缩空气吹除地面或工作台尘削是常用的方法，但这种方法既危险又徒劳无益。灰尘里的细小尘埃在扫地后不会立即落下。0.001mm的尘埃在静态空气中下降1米需要3.5小时，也就是说尘埃在一天的工作时间里大多游离于空中，且容易被人吸入（4-4）。通常，灰尘在空中不易被人察觉，但第二天灰尘就会落得到处都是。有效控制尘削的办法包括使用真空吸尘器和喷雾器。尘削湿润后很容易用扫帚或水冲干净。有效利用现场通风设备：现场通风设备只能是在其他方法无效时用于降低化学危险的办法。很多企业花大量的钱安装通风设备来控制危险物质，但往往事倍功半。其实有很多方法既经济又高效，有些非常实用。1正确使用风扇：风扇用于清除工作场所的危险物质。有害空气从通风处吹走（4-5）。有些情况下灰尘可吹入一个积尘器（4-6），注意以下几点：风扇和通风口之间不能有障碍物，任何障碍物都会极大的影响效果风扇转速要低以减少震荡污染空气排除不能途径工人从工作场所排除的有害气体不能危害企业以外的人员2使用排风和进风通风设备：用排风机、风扇或通风设施排出有害气体功效有限（4-7）