1、环境质量评价与系统分析环境质量评价与系统分析安徽工业大学建工学院安徽工业大学建工学院二二四年六月八日四年六月八日8.环境系统最优化8.1 环境规划和系统最优化8.2 线性规划的概念8.3 图解法解二维线性规划问题8.4 单纯形法解LP问题8.5 对偶线性规划模型8.6 Excel的规划求解8.7 规划求解在大气污染控制中应用8.1 环境规划和系统最优化8.1.1 城市环境规划n城市开发规划概要城市开发规划概要n(1)工业规划(2)自然环境改变(3)人口变化n土地利用规划土地利用规划n(1)总体规划(2)工业区划(3)居住区和商业区划(4)农业、林业和畜牧业等区划(5)其他依据和标准n水资源管理
2、规划水资源管理规划n(1)用水规划，水资源保护规划(水质、水量)(3)水面利用规划n城市能源规划城市能源规划n(1)能源利用规划(2)能源环境影响预测(3)能源环境管理规划n工业污染源控制规划工业污染源控制规划n(1)工业污染源环境影响预测 (2)控制规划n大气污染综合防治规划及其他大气污染综合防治规划及其他n(1)大气环境质量预测(2)大气污染防治(3)固体废物，化学品、噪声污染预测及防治n城市交通规划城市交通规划n城市绿化和建立生态调节区特殊保护区城市绿化和建立生态调节区特殊保护区8.1.2 环境系统最优化nMin Z=f(,)nS.t.G(,)=0),.,()(21nXXXFZMinMa
3、xmnmnnbXXXgbXXXgbXXXg或或或,),.,(.,),.,(,),.,(2122121211最优化模型可以写成更易于理解的一般形式：S.t 城市污水排水系统优化例)()()()(ttpppppupuppACACACACMinZ)2()1(0iiiitiQQQSiiCC minmaxminVVVDDDhHHDppminmaxmin0min22maxmin11max00)(VVVDDDDHDHFFzHEzzHEzqQADiiuuus.t.污水处理厂约束条件（总流量为各分流量之和，处理效率在工艺相应限制下）；水质约束条件（控制水质符合环境标准）；压力输水管约束条件（输水总扬程=输水净扬
4、程+水头损失,设计管径属于标准管径系列，最小管径限制，管中水流流速在最小允许流速和最大允许流速之间）。重力流污水管约束条件:(1)水量连续方程，Q为各管的设计流量,g本段流量；(2)(3)管段的上下游地面标高与管顶标高的差，在允许最大管顶覆土和最小覆土厚度之间；(4)水流最大充盈度限制；(5)相邻的上游管段的管底高程高于下游。其余同压力输水管相应约束条件 例例8-1o某金属冶炼厂,每生产 1kg 金属产生 0.3 kg废物，这些废物随废水排放,浓度为 2 kg/m3，废水经部分处理，排入附近河流。政府对废物实行总量控制，为10 kg/d。工厂最大生产能力为 5500kg/d,售价为$13/kg
5、，生产成本为$9/kg,废水处理设施的废水处理能力为700m3/d，处理费用是$2/m3,废水处理效率与污染物的负荷有关,以 Q 表示废水处理量,单位为(100m3/d)，处理效率为1-0.06Q，试对该问题建立最优化模型，并求解。0.03Y2工 厂污水处理厂0.3X0.3X-YY河流图图 8-1 污染物的发生与产量、处理量的关系。oMax Z400X-100YoS.t.0.3X-Y+0.03Y 2 10；X55；Y14；0.3X-Y0,X0,Y0 024681012141618200204060XY=14X=550.3X-Y+0.03Y2=10可 行 域ABCY0.3X-Y=0可行域曲线上的
6、目标值XYZ（$）3301333340215640454176275061929355921193图图 8-2 废水管理问题的可行域8.2 线性规划的概念o8.2.1 线性规划问题例例8-2 在上节讨论优化问题时,以水处理方案为例建立了最优化模型。该例中,污水处理效率与负荷有关,所以可行域边界线有一段为曲线.将例8-1的问题稍作修改,如果污水处理厂的处理效率与废水处理量无关，始终为=0.85，其他条件仍相同，该如何进行选择。解：设X:工厂的金属产量(100 kg/d)；Y:送往废水处理设施处理的污染物量(100 kg/d)；建立的最优化模型成为：Max Z 400X-100YS.t.0.3XY
7、+(1-0.85）Y 10；X55；Y14；0.3X-Y0,X0,Y0；oo8.2.2 线性规划问题的标准形式例例8-3农药管理问题。o一个容积为 100000m3的湖泊,湖水的平均停留时间为6个月,周围有1000ha 农田,农作物上施加的一部分农药会流失到湖中,并危害到吃鱼的鹰。环保部门想知道如何管理农田才不致对鹰造成危害，生物学的研究证明湖水中的农药在食物链中被富集,并按几何级数增长。设湖水中的农药浓度为 C 1(ppm),湖水中的藻类中的农药浓度为C 2(ppm),食藻鱼体内浓度为C 3(ppm),食鱼的鹰体内浓度为C 4(ppm),鹰的最大耐药浓度为100ppm。在1000ha农田上种
8、植两种农作物,它们具有不同的收益和农药施加量具体数据如下：oMax Z140X1+100X2oS.t.0.9X1+0.5X2 632.5oX1+X2 1000oX1,X20 作物农药施加量(kg/ha)农药流失率%作物收入$/ha作物费用$/ha蔬菜615300160粮食2.520150508.2.2 线性规划问题的标准形式).)(2211nnXcXcXcZMinMaxmnmnmmnnnnbXaXaXabXaXaXabXaXaXaLP,.,.,.)(221122222121112121110,.,21nXXXS.t.如果将不等式约束条件，全部使用“”表示，称为线性规划问题的典则形典则形式式。我
9、们还可以将一般形式转化为线性规划问题的标准形式标准形式。线性规划问题的标准形式可采用如下的矩阵表达式：0.XBAXCXtsMaxZ),.,(),.,(),.,(212121TmTnnbbbXXXcccBXCmnmmnnaaaaaaaaa.212222111211A其中8.3 图解法解二维线性规划问题o在线性规划问题中，如果只含有两个变量时，称为二维线性规划问题，就可以用图解法求解。o8.3.1 可行域和目标线可行域和目标线o线性规划问题图解法过程：o根据线性规划问题的约束条件，画出约束条件函数线，围出满足全部约束条件的解的可行域；o根据线性规划问题的目标函数，对确定的 Z 值（目标值可任意给定
10、），画出目标函数的投影线。变动Z 值，确定目标函数增大或减小的方向；o根据线性规划问题目标函数极大化或极小化要求，在线性规划问题解的可行域上平行移动目标函数投影线，找到平行线与可行域相接的最终边际点，确定问题的最优解。02468101214161820010203040506070XY=14X=550.3X-0.85Y=10可行域ABCY0.3X-Y=0z=13000z=21235DMax Z 400X-100YS.t.0.3XY+(1-0.85）Y 10；X55；Y14；0.3X-Y0,X0,Y0；Max Z140X1+100X2S.t.0.9X1+0.5X2 632.5 X1+X2 100
11、0X1,X2 0 020040060080010001200020040060080010001200AB(331.25,668.75)CX1+X2=10000.9X1+0.5X2=632.5z=110500z=102500z=113250X1X28.3.2 灵敏度分析由于实际问题中模型参数的不确定性，系统分析关心模型中选择的参数将会对最优解产生影响的灵敏程度，即参数在何等范围内变化时，原解仍然是合理的。020040060080010001200020040060080010001200ABC(702.78,0)X1+X2=10000.9X1+0.5X2=632.5X1X2z=50X1+100
12、X2=100000z=220X1+100X2=1546108.4 单纯形法解LP问题Max Z=140X1+100X2+0S1+0S2 632,5=0.9X1+0.5X2+S1+0S2 1000=X1+X2+0S1+S2 1.化LP问题为线性规划问题的标准形式(LP)LP:MaxZ=140X1+100X2S.t.0.9X1+0.5X2 632.5 X1+X2 1000 X1,X20引入松驰变量 S1,S2;得LP如下:2.选择初始极点,写出未经迭代的单纯形表(零级单纯形表)Z 0140X1 100X2 S1632.50.9X1 0.5X2 S21000X1X2 3.迭代o 使零级单纯形表中的一
13、个零变量增值，以改进目标函数，变量的选择以高效为原则，即该变量的增值，会使目标函数得到最迅速的增加。o 4.重复迭代过程 Z98388.8889-155.555556S1+22.22222X2X1702.77778-1.111111S1-0.555556X2S2297.22222+1.111111S1-0.444444X2Z113250.147-99.99957S1-50.00044S2X1668.75668+2.50002S1-2.25002S2X2331.24332-2.50002S1+1.25002S25.判别新解是否为最优解o 由目标函数的系数进行判别:S1和S2的系数均小于零,说明这
14、些变量的增大不会改善Z值，目标函数已经是最优解。因此，对于单纯形表，如果目标函数关系式中的全部变量系数(检验数)均小于或等于零时，目标函数已经是最优解，整个迭代过程结束。o 从单纯形表可解读出 LP 的解,o S2=0,S1=0,X1=331.24,X2=668.76 和 Z=113250.15。8.5 对偶线性规划模型o 8.5.1 由算例认识对偶问题1.原模型原模型LP:Max Z=140X1+100X2S.t.0.9X1+0.5X2 632.5 X1+X2 1000X1,X202.用边际值定义新变量用边际值定义新变量设 Y1 为改变农药限制条件的边际值($/kg)，定义Y2为总种植面积的
15、边际值($/ha)，3.农药管理问题的对偶问题农药管理问题的对偶问题MinZ632.5Y1+1000Y2 0.9Y1+1.0Y2 140 0.5Y1+1.0Y2 100 Y1,Y20 农药管理对偶问题的目标线和可行域 020406080100120140160050100150200Y1Y20.9Y1+Y2=1400.5Y1+Y2=100Z=632.5Y1+1000Y2=113250(100,50)可行域MinZ632.5Y1+1000Y2 0.9Y1+1.0Y2 140 0.5Y1+1.0Y2 100 Y1,Y20 两种不同目标函数下的原始和对偶模型方案原 始 模 型对 偶 模 型1Max
16、Z=220X1+100X2S.t.0.9X1+0.5X2 632.5X1+X2 1000X1,X20Z154610X1702.78,X20MinZ632.5Y1+1000Y20.9Y1+1.0Y2 2200.5Y1+1.0Y2 100Y1,Y20Z154610Y1244.44,Y202Max Z=50X1+100X2S.t.0.9X1+0.5X2 632.5X1+X2 1000X1,X20Z100000X10,X21000MinZ632.5Y1+1000Y20.9Y1+1.0Y2 1400.5Y1+1.0Y2 100Y1,Y20Z100000 Y10 ,Y2100线性规划的对偶模型 3.线性规
17、划原问题和对偶问题的基本性质线性规划原问题和对偶问题的基本性质原型问题是一个极大化规划，则对偶问题是一个极小化规划；在典则形式中，若极大化的原型问题是“”约束条件，则在极小化对偶规划中是“”约束条件；原型问题约束方程的个数等于对偶变量的个数，反之亦然；原型问题约束方程右边的常数分别为对偶问题目标函数的系数，反之亦然；线性规划对偶问题的对偶是原问题；若是线性规划原问题的可行解,是对偶问题的可行解,则，等于的情况出现在原问题和对偶问题的最优解；对偶问题的最优解可由原设问题获得时，此最优解可在单纯形表的松弛变量中的检验数中得到，反之，原设问题的最优解可由对偶问题的剩余变量各列中的检验效中得到；对偶定
18、理：若原问题有最优解,那么对偶问题也有最优解,且其目标函数值相等。研究线性规划的对偶问题具有实用意义。通过线性规划对偶问题的研究，能够更加深入认识线性规划问题变量和约束的内容；对变量少但约束多的线性规划问题可通过对偶转换简化求解；应用对偶模型还可以讨论线性规划问题中参数的灵敏度。8.5.3 影子价格o 在典则形式的对偶问题中，对偶LP 模型的变量Y1,Y2,.Ym0,称为对偶变量。o 对偶变量有其经济解释，被称为影子价格。对偶变量Yi 是相应于第 i个原始约束条件的边际值，它是第 i 种资源每一个单位对目标值的贡献。o 影子价格并不是市场价格，它是在特定最优解条件下某一资源的潜在价格，反映这种
19、资源在实现最优解中的作用和紧迫程度。8.6 Excel的规划求解o例例8-6 用Excel的规划求解，解下列非线性问题。Max Z 8X1+12X2+4X3 S.t.X13+4X2+3X3327X1-X22+3X32 X1，X2，X3 0例例8-5 用Excel的规划求解，解农药管理问题。解：(1)由原模型LP:Max Z=140X1+100X2S.t.0.9X1+0.5X2 632.5 X1+X2 1000 X1,X20例例8-7 有九头鸟、鸡、兔同笼，上有152头，下有374脚，问笼中九头鸟、鸡、兔各多少？Excel的规划求解操作过程1在“工具”菜单中，单击“规划求解”命令。如果“规划求解
20、”命令没有出现在“工具”菜单中，则需要安装“规划求解”加载宏。2在“目标单元格”编辑框中，键入单元格引用或目标单元格的名称。目标单元格必须包含公式。3目标单元格中数值可选“最大值”、“最小值”或指定“目标值”。4在“可变单元格”编辑框中，键入每个可变单元格的名称或引用，用逗号分隔不相邻的引用。可变单元格必须直接或间接与目标单元格相联系。5在“约束”列表框中，输入相应的约束条件。约束条件是指“规划求解”问题中设置的限制条件。6单击“求解”按钮。7如果要在工作表中保存求解后的数值，请在“规划求解结果”对话框中，单击“保存规划求解结果”。8报告（1）运算结果报告（2）敏感性报告（3）极限值报告 8.
21、7 规划求解在大气污染控制中应用o 8.7.1 污染物排放控制n例例8-9 某区域内有三个排放总悬浮颗粒物(TSP)的点源，其中两个是燃煤发电厂，另一个是水泥厂的窑炉。发电厂每烧一吨煤排放 95kg的TSP，水泥厂每生产一吨水泥排放85kg的TSP。水泥厂的产量为250000 t/a水泥，两个燃煤发电厂的燃煤量分别是 400000 和 300000t/a。目前三个点源均无控制措施，环保机构希望将该地区的TSP削减 80%。点源去除TSP的可行方法去除效率和相应费用列于表，现需通过成本-效果分析，以最小费用达到环境目标。污染控制方法和费用控制方法去除效率电厂1($/t)电厂2($/t)水泥厂($
22、/t)隔板沉淀槽591.01.41.1多级除尘器74-1.2长锥除尘器84-1.5喷雾洗涤器942.02.23.0静电除尘器972.83.0-决策变量定义方法j 控制方法1电厂12电厂23水泥厂不控制隔板沉淀槽多级除尘器长锥除尘器喷雾洗涤器静电除尘器X10X11-X14X15X20X21-X24X25X30X31X32X33X34-大气污染控制的目的是使 TSP 总排放量削减 80%，用燃煤生产排放 TSP系数 95 kg/t 和水泥生产排放 TSP系数 85kg/t求出各污染源未加控制的TSP 的总排放量：污染源1:40000095=38000000 kg/a污染源2:30000095=28
23、500000 kg/a污染源3:25000085=21250000 kg/a总计 :87750000 kg/a最终限量:17550000 kg/a控制后污染源1：95X10+95(0.41)X11+95(0.06)X14+95(0.03)X15=95X10+39X11+5.7X14+2.9X15 不考虑生产水平的改变:X10+X11+X14+X15400000X20+X21+X24+X25300000X30+X31+X32+X33+X34250000 大气污染控制的规划模型Min Z1.0X11+2.0X14+2.8X15+1.4X21+2.2X24+3.0X25+1.1X31+1.2X32+
24、1.5X33+3.0X34S.t.X10+X11+X14+X15400000X20+X21+X24+X25300000X30+X31+X32+X33+X3425000095X10+39X11+5.7X14+2.9X15+95X20+39X21+5.7X24+2.9X25+85X30+34.9X31+22.1X32+13.6X33+5.1X3417600000Xij0,i=1,2,3；j=0,1,2,3,4,5Excel规划求解ABCDEFGHIJKLMNOP1X10X11X14X15X20X21X24X25X30X31X32X33X341042可变量02416001290002230Z3目标值
25、02432001.463.800026.44.501524源102416000000000040405源20000012900000030306源30000000000223025257TSP控093691003916500048640.80175917608目标值0122.801.42.2301.11.21.539源11111000000000系10源20000111100000数11源30000000011111矩12TSP控95395.72.995395.72.9853522.113.65.1阵138.7.2 大气质量管理的污染物迁移模型)22(exp)0,(2222zyzyHyuQyx
26、CC(x,y,0):点(x,y)的大气污染物地面浓度(g/m3)；y:水平方向上(Y)烟羽浓度的标准差（m）；z:垂直方向上(Z)烟羽浓度的标准差（m）；Q:污染源源强 g/s）；u:平均风速（m/s）；H:有效烟囱高度（m）对于单位源强，大气污染物的地面浓度仅仅是位置关系和气象条件的函数。污染源 i 在下风向任一点 k 处造成的大气污染物的地面浓度)22(exp12222ziyikzyikHyut当气象条件已知，污染源和采样考核点的位置相对固定时，tik 成为常数。多源在接受点k 处的最终污染物浓度，是每个污染源单独作用情况的叠加。由 m个污染源采用 n 种控制方法，在接受器 k 处形成的污
27、染物总浓度 为：minjijijpikkpXbtC11通用大气质量管理的污染迁移模型 Xij:污染源 i 采用控制方法 j 的产量或燃料量；Cij:污染源 i 采用控制方法 j 的年费用；aij:为一时表示控制方法 j 对污染源 i可行,否则为零；Si:污染源 i 应达到的总产量或燃料量；bijp:污染源 i,采用控制方法 j 时,污染物 p 的排放系数；tik:污染源 i 达到接受器 k 处的污染物迁移因子（单位源强浓度系数）；C0pk:接受器 k 处,污染物 p 的大气环境质量浓度标准。这是一个线性规划模型，共具有 m 个污染源，n 种控制方法,q 种污染物,且在r个接受器处，应达到大气环
28、境质量浓度标准。；miSXatsXCMinZinjijijminjijij,.,2,1.111 ,0r.,1,2,.k.,.,2,1011；jiXijqpCXbtpkminjijijpik第八章环境系统最优化学习要点o 本章主要讲述环境系统最优化方法、环境问题的线性规划求解算法等内容。1.应用系统分析方法解决环境问题的显著特点是通过模型化和最优化来协调环境系统中各要素之间的关系，实现经济效益、环境效益和社会效益的统一。常用的最优化方法有线性规划、动态规划与网络分析等。2.认识线性规划问题的一般形式、典则形式和标准形式。掌握线性规划问题的基本概念：目标函数、约束条件、可行域、目标线、灵敏度、搜索
29、策略、方案和目标值等。3.掌握线性规划问题的图解法、单纯形法等基本解法。4.掌握线性规划问题的对偶模型和灵敏度分析，并能利用影子价格认识其在原线性规划问题中所具有的意义。5.掌握运用 Excel求解规划问题的方法。根据规划问题的基本概念，正确解读出用 Excel生成的运算结果报告、敏感性报告和极限值报告等。熟练使用 Excel 规划求解解决大气污染控制的排污问题和污染物扩散问题。难点重点重点9.附录n 9.1 环境评价的法规与条例名录n 9.2 环境影响评价技术导则n 9.3 环境空气质量标准n 9.4 地表水环境质量标准n 9.5 大气污染物综合排放标准n 9.6 污水综合排放标准n 参考文献
