重金属扩散模型

重金属扩散模型
重金属扩散模型

重金属系指密度4.0以上约60种元素或密度在5.0以上的45种元素。砷、硒是非金属,但是它的毒性及某些性质与重金属相似,所以将砷、硒列入重金属污染物范围内。环境污染方面所指的重金属主要是指生物毒性显著的汞、镉、铅、铬以及类金属砷,还包括具有毒性的重金属锌、铜、钴、镍、锡、钒等污染物。随着全球经济化的迅速发展,含重金属的污染物通过各种途径进入土壤,造成土壤严重污染。土壤重金属污染可影响农作物产量和质量的下降,并可通过食物链危害人类的健康,也可以导致大气和水环境质量的进一步恶化。因此引起世界各国的广泛重视。目前,世界各国土壤存在不同程度的重金属污染,全世界平均每年排放Hg约1.5万t、Cu为340万t、Pb为500万t、Mn为1500万t、Ni 为100万t[1]。中国北方大城市的蔬菜基地和部分商品粮基地也存在着不同程度的重金属污染,如北京、天津、西安、沈阳、济南、长春、郑州等地;。

南方相对较轻,如福州、宁波、上海、武汉、成都等地。土壤重金属污染将会造成生态系统的严重破坏。从中国土壤资源状况看,到2000年底中国人均耕地仅为0.1 hm2,而且随着今后中国经济社会的发展如生态退耕、农业结构调整及自然灾害损毁等,土壤资源将进一步减少。因而如何有效地控制及治理土壤重金属的污染,改良土壤质量,将成为生态环境保护工作中十分重要的一项内容。

本文主要从土壤中重金属污染物来源与分布、土壤中重金属污染物的现行治理方法入手,提出土壤中重金属污染物防治的环境矿物学新方法。旨在保护环境,提高土壤的环境质量。

1 土壤中重金属污染物来源与分布

土壤中重金属的来源是多途径的,首先是成土母质本身含有重金属,不同的母质、成土过程所形成的土壤含有重金属量差异很大。此外,人类工农业生产活动,也造成重金属对大气、水体和土壤的污染。

1.1 大气中重金属沉降

大气中的重金属主要来源于工业生产、汽车尾气排放及汽车轮胎磨损产生的大量含重金属的有害气体和粉尘等。它们主要分布在工矿的周围和公路、铁路的两侧。大气中的大多数重金属是经自然沉降[2]和雨淋沉降进入土壤的。如瑞典中部Falun市区的铅污染[3],它主要来自于市区铜矿工业厂、硫酸厂、油漆厂、采矿和化学工业产生大量废物,由于风的输送,这些细微颗粒的铅,从工业废物堆扩散至周围地区。南京某生产铬的重工业厂[4]铬污染叠加已超过当地背景值4.4倍,污染以车间烟囱为中心,范围达1.5 km2,污染范围最大延伸下限1.38 km。俄罗斯的一个硫酸生产厂[5]也是由工厂烟囱排放造成S、V、As的污染。

公路、铁路两侧土壤中的重金属污染,主要是Pb、Zn、Cd、Cr、Co、Cu的污染为主。它们来自于含铅汽油的燃烧,汽车轮胎磨损产生的含锌粉尘等。它们成条带状分布,以公路、铁路为轴向两侧重金属污染强度逐渐减弱;随着时间的推移,公路、铁路土壤重金属污染具有很强的叠加性。在宁—杭公路南京段[6]两侧的土壤形成Pb、Cr、Co污染晕带,且沿公路延长方向分布,自公路向两侧污染强度减弱。在宁—连一级公路淮阴段[7]两侧的土壤铅含量增高,向两侧含量逐渐降低,且在地表0~30 cm铅的含量较高。在法国索洛涅地区A71号高速公路[8]沿途严重污染重金属Pb、Zn、Cd,其沉降粒子浓度超过当地土壤背景值

2~8倍,而公路旁重金属浓度比沉降粒子中高7~26倍。在斯洛文尼亚[9]从居波加到扎各瑞波公路两侧,铅除了分布在公路两侧以外,还受阶地地貌和盛行风的影响,高铅出现在低地,公路顺风一侧铅含量较高。

经过自然沉降和雨淋沉降进入土壤的重金属污染,主要以工矿烟囱、废物堆和公路为中心,向四周及两侧扩散;由城市—郊区—农区,随距城市的距离加大而降低,特别是城市的郊区污染较为严重。此外,还与城市的人口密度、城市土地利用率、机动车密度成正相关;重工业越发达,污染相对就越严重。

此外,大气汞的干湿沉降[10~12]也可以引起土壤中汞的含量增高。大气汞通过干湿沉降进入土壤后,被土壤中的粘土矿物和有机物的吸附或固定,富集于土壤表层,或为植物吸收而转入土壤,造成土壤汞的浓度的升高。

1.2 农药、化肥和塑料薄膜使用

施用含有铅、汞、镉、砷等的农药和不合理地施用化肥,都可以导致土壤中重金属的污染。一般过磷酸盐中含有较多的重金属Hg、Cd、As、Zn、Pb,磷肥次之,,氮肥和钾肥含量较低,但氮肥中铅含量较高,其中As和Cd污染严重[13]。经过对上海地区菜园土地、粮棉地的研究[14],施肥后,Cd的含量从0.134 mg/kg 升到0.316 mg/kg,Hg的含量从0.22 mg/kg升到0.39 mg/kg,Cu、Zn 增长

2/3。通过新西兰[15]50 a前和现今同一地点58个土样分析,自施用磷肥后,镉从0.39 mg/kg升至0.85 mg/kg。在阿根廷[16]由于传统无机磷肥的施入,进而导致土壤重金属Cd、Cr、Cu、Zn、Ni、Pb的污染。

农用塑料薄膜生产应用的热稳定剂中含有Cd、Pb,在大量使用塑料大棚和地膜过程中都可以造成土壤重金属的污染。

1.3 污水灌溉

污水灌溉一般指使用经过一定处理的城市污水灌溉农田、森林和草地。城市污水包括生活污水、商业污水和工业废水。由于城市工业化的迅速发展,大量的工业废水涌入河道,使城市污水中含有的许多重金属离子,随着污水灌溉而进入土壤。在分布上,往往是靠近污染源头和城市工业区土壤污染严重,远离污染源头和城市工业区,土壤几乎不污染[17]。近年来污水灌溉已成为农业灌溉用水的重要组成部分,中国自60年代至今,污灌面积迅速扩大,以北方旱作地区污灌最为普遍,约占全国污灌面积的90%以上。南方地区的污灌面积仅占6%,其余在西北和青藏[18]。污灌导致土壤重金属Hg、Cd、Cr、As、Cu、Zn、Pb等含量的增加。淮阳污灌区自污灌以来,金属Hg、Cd、Cr、Pb、As等就逐渐增高,1995~1997年已超过警戒级[19]。太原污灌区的重金属Pb、Cd、Cr含量远远超过其当地背景值,且积累量逐年增高[20]。

1.4 污泥施肥

污泥中含有大量的有机质和氮、磷、钾等营养元素,但同时污泥中也含有大量的重金属,随着大量的市政污泥进入农田,使农田中的重金属的含量在不断增高。污泥施肥可导致土壤中Cd、Hg、Cr、Cu、Zn、Ni、Pb含量的增加,且污泥施用越多,污染就越严重,Cd、、Cu、Zn引起水稻、蔬菜的污染;Cd、Hg可引起小麦、玉米的污染;污泥增加,青菜中的Cd、Cu、Zn、Ni、Pb也增加[21]。Anthony[22]研究表明,用城市污水、污泥改良土壤,重金属Hg、Cd、Pb等的含量也明显增加。

1.5 含重金属废弃物堆积

含重金属废弃物种类繁多,不同种类其危害方式和污染程度都不一样。污染的范围一般以废弃堆为中心向四周扩散。通过对武汉市垃圾堆放场[23]、杭州某铬渣堆存区[24]、城市生活垃圾场[25]及车辆废弃场[26]附近土壤中的重金属污染的研究,这些区域的重金属Cd、Hg、Cr、Cu、Zn、Ni、Pb、As、Sb、V、Co、Mn的含量高于当地土壤背景值,重金属在土壤中的含量和形态分布特征受其垃

圾中释放率的影响,且随距离的加大重金属的含量而降低。由于废弃物种类不同,各重金属污染程度也不尽相同,如铬渣堆存区的Cd、Hg、Pb为重度污染,Zn 为中度污染,Cr、Cu为轻度污染。

1.6 金属矿山酸性废水污染

金属矿山的开采、冶炼、重金属尾矿、冶炼废渣和矿渣堆放等,可以被酸溶出含重金属离子的矿山酸性废水,随着矿山排水和降雨使之带入水环境(如河流等)或直接进入土壤,都可以间接或直接地造成土壤重金属污染。1989年我国有色冶金工业向环境中排放重金属Hg为56 t,Cd为88 t,As为173 t,Pb为226 t[27]。矿山酸性废水重金属污染的范围一般在矿山的周围或河流的下游,在河流中不同河段的重金属污染往往受污染源(矿山)控制,河流同一污染源的下段自上游到下游,由于金属元素迁移能力减弱和水体自净化能力的适度恢复,金属化学污染强度逐渐降低。江西乐安江沽口—中洲[28]由于遭受德兴铜矿的污染,水体及土壤中的重金属Cu、Pb、Zn、Cr含量增高,至鄱阳湖段重金属含量逐渐降低。美国科罗拉多州罗拉多流域[29]受采矿的影响,重金属元素Cd、Zn、Pb、As

的浓度,以污染源为最高,之后随着与污染源距离延长而逐渐降低。莱安河[30]重金属污染,来自一个大型铜矿,导致重金属浓度远远超过当地背景值。流域重金属污染随季节变化而异,枯水期重金属的含量明显高于丰水期[31]。河流流速减缓可以导致该流段重金属含量增加[32]。

同一区域土壤中重金属污染物的来源途径可以是单一的,也可以是多途径的。胡永定[33]通过研究徐州荆马河区域土壤重金属污染的成因中指出:Cr、Cu、Zn、Pb是由垃圾施用引起的,As是由农灌引起的,Cd是由农灌和垃圾施用引起的,Hg是各种途径都具备。王文祥[34]通过对山东省耕地重金属元素污染状况的研究说明,工业快速发展地区铅高于农业环境,铅与距公路远近有关。乡镇企业技术、设备落后,原材料利用率低,造成其周边土壤重金属污染相当严重。据贵州1986年的统计,全省乡镇排放汞14.7万kg,土壤中有的地方达56.64 mg/kg,超过未污染土壤的84.5倍。要引起高度重视。

总的来说:工业化程度越高的地区污染越严重,市区高于远郊和农村,地表高于地下,污染区污染时间越长重金属积累就越多,以大气传播媒介土壤重金属污染土壤的具有很强的叠加性,熟化程度越高重金属含量越高。

2 土壤中重金属污染物现行治理方法

关于土壤重金属污染物的研究,国外始于20世纪60~70年代,如澳大利亚、美国、德国等国家对土壤重金属较深入,尤其澳大利亚。我国在1983年对主要类型的土壤环境容量作过初步研究,如提出研究土壤重金属的生态效应、临界含量地带性分异规律和分区等。

当前,世界各国很重视对重金属污染治理方法研究,并开展广泛的研究工作[35~39]。总的来说,目前大致有以下四种治理措施:

2.1 工程治理方法

工程治理是指用物理或物理化学的原理来治理土壤重金属污染。主要有:客土是在污染的土壤上加入未污染的新土;换土是将以污染的土壤移去,换上未污染的新土;翻土是将污染的表土翻至下层;去表土是将污染的表土移去等。如日本富士县神通川流域的痛痛病发源地,就是由于长期食用含镉的稻米而引发的,他们通过研究,去表土15 cm,并压实心土,在连续淹水的条件下,稻米中镉的含量小于0.4 mg/kg;去表土后再客土20 cm,间歇灌溉稻米中镉的含量也不超标,

客土超过30 cm,其效果更佳。此外淋洗法是用淋洗液来淋洗污染的土壤;热处理法是将污染土壤加热,使土壤中的挥发性污染物(Hg)挥发并收集起来进行回收或处理;电解法是使土壤中重金属在电解、电迁移、电渗和电泳等的作用下在阳极或阴极被移走。

以上措施具有效果彻底、稳定等优点,但实施复杂、治理费用高和易引起土壤肥力降低等缺点。

2.2 生物治理方法

生物治理是指利用生物的某些习性来适应、抑制和改良重金属污染。主要有:动物治理是利用土壤中的某些低等动物蚯蚓、鼠类等吸收土壤中的重金属;微生物治理是利用土壤中的某些微生物等对重金属具有吸收、沉淀、氧化和还原等作用,降低土壤中重金属的毒性如Citrobacter sp产生的酶能使U、Pb、Cd形成难溶磷酸盐;原核生物(细菌、放线菌)比真核生物(真菌)对重金属更敏感,格兰氏阳性菌可吸收Cd、Cu、Ni、Pb等[44]。植物治理是利用某些植物能忍耐和超量积累某种重金属的特性来清除土壤中的重金属;重金属的植物吸收、淋溶和无效态数量将只依赖于它们的有效态的多少,重金属溶液浓度和它们的土壤的有效态之间关系遵循Freundlich吸附方程[41];超积累植物可吸收积累大量的重金属,目前已发现400多种,超积累植物积累Cr、Co、Ni、Cu、Pb的含量一般在0.1%以上,积累Mn、Zn含量一般在1%以上[40];印度芥菜(Brassica juncea)可吸收Zn、Cd、Cu、Pb等,在Cu为250 mg/kg,Pb为500 mg/kg、Zn为500 mg/kg条件下能生长,在Cd为200 mg/kg出现黄化现象[42];印度芥菜(Brassica juncea)可对Cr6+、Cd、Ni、Zn、Cu富集分别为58,52,31,17和7倍[45];高杆牧草(Agropyron elongatum)能吸收Cu等[43];英国的高山莹属类等,可吸收高浓度的Cu、Co、Mn、Pb、Se、Cd、Zn等。

生物治理措施的优点是实施较简便、投资较少和对环境破坏小,缺点是治理效果不显著。

2.3 化学治理方法

化学治理就是向污染土壤投入改良剂、抑制剂,增加土壤有机质、阳离子代换量和粘粒的含量,改变pH、Eh和电导等理化性质,使土壤重金属发生氧化、还原、沉淀、吸附、抑制和拮抗等作用,以降低重金属的生物有效性。其中沉淀法是指土壤溶液中金属阳离子在介质发生改变(pH值、OH-、SO42-等)时,形成金属沉淀物而降低土壤重金属的污染;如向土壤中投放钢渣,它在土壤中易被氧化成铁的氧化物,对Cd、Ni、Zn的离子有吸附和共沉淀作用,从而使金属固定。在沈阳张士污灌区进行的大面积石灰改良实验表明,每公顷施石灰1500~1875 kg籽实含镉量下降50%[18]。有机质法是指有机质中的腐殖酸能络合重金属离子生成难溶的络合物,而减轻土壤重金属的污染;吸附法是指重金属离子能被膨润土、沸石、粘土矿物等吸附固定,从而降低土壤重金属的污染。

化学治理措施优点是治理效果和费用都适中,缺点是容易再度活化。

2.4 农业治理方法

农业治理是因地制宜的改变一些耕作管理制度来减轻重金属的危害,在污染土壤上种植不进入食物链的植物。主要有:控制土壤水分是指通过控制土壤水分来调节其氧化还原电位(Eh),达到降低重金属污染的目的;选择化肥是指在不影响土壤供肥的情况下,选择最能降低土壤重金属污染的化肥;增施有机肥是指有机肥能够固定土壤中多种重金属以降低土壤重金属污染的措施;选择农作物品种是指选择抗污染的植物和不要在重金属污染的土壤上种植进入食物链的植物;如

在含镉100 mg/kg的土壤上改种苎麻,五年后,土壤镉含镉平均降低27.6%[46];因地制宜地种植玉米、水稻、大豆、小麦等,水稻根系吸收重金属的含量占整个作物吸收量的[35]58%~99%,玉米茎叶吸收重金属的含量占整个作物吸收量的20%~40%,玉米籽实吸收量最少,重金属在作物体内分配规律是根>茎叶>籽实[47]。土壤重金属污染也是导致生态系统破坏的重要因素。合理的利用农业生态系统工程措施,也可以保持土壤的肥力,改良和防治土壤重金属污染,提高土壤质量,并能与自然生态循环和系统协调运作。如可以在污染区公路两侧尽可能种树、种花、种草或经济作物(如蓖麻),种植草皮或观赏树木,移栽繁殖,不但可以美化环境,还可以净化土壤;蓖麻可用作肥皂的原料。也可以进行农业改良,即在污染区繁育种子(水稻、玉米),之后在非污染区种植;或种植非食用作物(高梁、玉米),收获后从秸秆提取酒精,残渣压制纤维板,并提取糠醛,或将残渣制作沼气作能源。

农业治理措施的优点是易操作、费用较低,缺点是周期长、效果不显著。

3 土壤中天然矿物治理重金属污染物新方法

土壤的主要矿物组成除粘土矿物外,还存在大量的天然铁锰铝氧化物及氢氧化物、硅氧化物、碳酸盐、有机质硫化物等天然矿物。在国内外关于土壤重金属污染物防治途径研究中,人们一直强调土壤自身的净化能力,但土壤自净化能力离不开土壤中矿物种对重金属的吸附与解吸作用、固定与释放作用,土壤中具体矿物的净化能力才真正体现土壤自身的净化能力和容纳能力。土壤中有毒有害元素含量的高低,并不是直接判定土壤环境质量优劣乃至土壤生态效应的唯一标志,关键问题是要揭示这些重金属在土壤中与各种无机物之间具有怎样的环境平衡

关系。在国内外为寻求地下水和土壤有机污染的修复方法而直接对土壤中多种粘土矿物进行改性研究,即利用有机表面活性剂去置换天然粘土矿物中存在着的大量可交换的无机阳离子,以形成有机粘土矿物,可有效截住或固定有机污染物,阻止地下水的进一步污染,限制有机污染物在土壤环境中迁移扩散。但特别需要指出的是,在粘土矿物改性过程中,其中的固定态重金属也一并被置换出来,导致土壤系统中业已建立环境平衡被打破,使得土壤环境中解吸释放态重金属污染物总量大大增加。至此,土壤中重金属污染物既来源于土壤中活动态的重金属,又来源于改性粘土矿物时被置换释放出来的重金属。以本实验室正在开展研究的环境矿物材料[48]—天然铁锰铝氧化物及氢氧化物为例[49, 50],其中磁铁矿、赤铁矿、针铁矿、软锰矿、硬锰矿与铝土矿等也正在成为国际上关于天然矿物净化污染方法研究方面的重点对象之一[51]。我们认为天然铁锰铝氧化物及氢氧化物的表面具有明显的化学吸附性特征,锰氧化物与氢氧化物还具有较完善的孔道特征,尤其是Fe、Mn为自然界中少数的但属于常见的变价元素,其氧化物和氢氧化物化合物往往可表现出一定的氧化还原作用。所以说天然铁锰铝氧化物及氢氧化物具有潜在的净化重金属污染物的功能,能成为土壤环境中吸附固定态重金属污染物的有效物质。

综上所述,国内外对土壤重金属污染现状与治理,取得了一定的成绩,也存在一些理论上和技术上的问题,如土壤中重金属与土壤中矿物之间的吸附与解吸、固定与释放的平衡关系的研究,土壤中重金属形态特征、转化与迁移规律的系统研究,土壤中二次污染物的及时处理等。

土壤重金属污染首先应从源头抓起,控制污染源,土壤重金属的污染已经达到相当严重的程度,要充分认识土壤重金属污染的长期性、隐匿性、不可逆性以及不

能完全被分解或消逝的特点。土壤质量问题是经济可持续发展和社会全面进步的战略问题,它直接影响土壤质别、水质状况、作物生长、农业产量、农产品品质等,并通过食物链对人体健康造成危害。对工业生产中排放的污染物尚未得到较彻底控制,尤其在农业生产中大量而盲目使用化学肥料和农药的今天,江河湖海、地下水及陆地中无机和有机污染物积累总量与日俱增,使土地环境质量变得极其脆弱。一旦土壤对这些污染物尤其是重金属的消纳容量达到饱和,这些污染物对耕地生产能力的潜在毁灭性破坏便有可能一触即发,有人已形象地称之为农业生产中的“定时炸弹”。从这个意义上来讲,土地管理与保护工作不仅是对耕地总量的监管,还应该加强对耕地质量的保护与改善。对土壤质量的保护便是对耕地生产能力的保护,更是提高土地利用效率的强有力措施之一。对于我国这样一个人口众多的农业大国,开展国土质量调查评价,对土壤重金属污染物进行试验研究,开发耕地污染的治理方法和技术,显得更为必要和迫切

大气污染物扩散模式

第四章 大气扩散浓度估算模式 第一节 湍流扩散的基本理论 一 湍流 1.定义:大气的无规则运动 风速的脉动 风向的摆动 2.类型: 按形成原因 热力湍流:温度垂直分布不均(不稳定)引起,取决于大气稳定度 机械湍流:垂直方向风速分布不均匀及地面粗糙度引起 3.扩散的要素 风:平流输送为主,风大则湍流大 湍流:扩散比分子扩散快105~106倍 二 湍流扩散理论(主要阐述湍流与烟流传播及湍流与物质浓度衰减的关系) 1.梯度输送理论 通过与菲克扩散理论类比建立起来的(菲克定律:单位时间内通过单位断面上的物质的数量与浓度梯 度呈正比) 类比于分子扩散,污染物的扩散速率与负浓度梯度成正比 x C k F ??-= 式中,F — 污染物的输送通量 k — 湍流扩散系数 C — 污染物的浓度 X — 与扩散截面垂直的空间坐标(扩散过程的长度) x C ??— 浓度梯度 要求得各种条件下某污染物的时、空分布,由于边界条件往往很复杂,不能求出严格的分析解,只能是在特定的条件下求出近似解,再根据实际情况进行修正。 2.湍流统计理论 泰勒首先将统计理论应用在湍流扩散上 图4-1显示:从原点O 放出的粒子,在风沿着x 方向吹的湍流大气中扩散。粒子的位置用y 表示,则结论为: ①y 随时间变化,但其变化的平均值为零 ②若从原点放出很多粒子,则在x 轴上粒子的浓度最高,浓席分布以x 轴为对称轴,并符合正态分布。 萨顿实用模式:解决污染物在大气中扩散的实用模式 高斯模式:应用湍流统计理论得出正态分布假设下的扩散模式 3.相似理论 第二节 高斯扩散模式 一 坐标系的建立—右手坐标系

1.原点O :无界点源或地面源,O 为污染物的排放点 高架源,O 为污染物的排放点在地面上的投影点 补充:点源 高架源 连续源 固定源 线源 地面源 间歇源 流动源 面源 2.x 轴:正向为平均风向,烟流中心线与x 轴重合 3.y 轴:垂直于x 轴 4.z 轴:垂直于xoy 平面 二 高斯模式的有关假定 1.污染物浓度在y 、z 轴上的分布为正态分布; )2exp(21 )(22 y y y y f σπ σ-= )2exp(21 )(22 z z z z f σπ σ-= y σ,z σ— 分别为污染物在y 和z 方向上分布的标准差,m 2.全部高度风速均匀稳定,即风速u 为常数; 3.源强是连续均匀稳定的,源强Q 为定值; 4.扩散中污染物是守恒的,不考虑转化,即烟云在扩散过程中没有沉降、化合、分解及地面吸收、吸附作用发生; 0=??t C 5.在x 方向上,输送作用远远大于扩散作用,即 )(x C k x x C u x ????>>??; 6.地面足够平坦。

气体泄漏及扩散计算

学号: 07412225 常州大学 毕业设计(论文) (2011届) 题目重气泄漏扩散模拟及应急救援 学生薛云龙 学院环境与安全工程学院专业班级安全072班 校内指导教师王新颖专业技术职务讲师 校外指导老师专业技术职务 二○一一年六月

重气泄漏扩散模拟及应急救援 摘要:重气泄漏扩散事故是经常发生且危害较大的一种事故形式,由于重气的密度大于空气,因此重气往往沿地面扩散,泄放物质进入人体将引起中毒事故,若泄放物质被点燃或引爆将引起大规模的燃烧爆炸事故。虽然人们对重气泄漏扩散所造成的危害十分重视,但由于缺乏足够有效的数据来提供人们作风险评估及预防改善措施,因此采用数学模型进行模拟是必要的。应在生产过程中,加强管理,强化生产者的安全生产教育。分析了泄漏扩散事故的七大影响因素,提取并建立了泄漏事故模式,并对各种事故模式的泄漏机理和发生条件进行了研究分析。通过试验研究得出在实际环境中大气主导风的风速,泄漏方向对气体扩散浓度分布有重大的影响,泄漏气体在下风向扩散的最快。静风时,随着时间的增加,空间各点的浓度有升高的趋势;在稳定风流中,空间各点的浓度随时间的变化不明显,可以认为是稳态的。泄漏的气体在下风向扩散的最快,在现场一旦发生天燃气泄漏,应综合考虑泄漏源的方向和该点当时的风向,风速等因素,及时准确预测泄漏气体可能扩散到危险区域,做好应对措施。 关键词:相似理论;泄漏模型;泄漏扩散模式;示踪法;重气;应急救援;

Heavy gas leak dispersion modeling and emergency rescue Abstract : As it is well-known, many industrial and domestic gases are toxic and flammable are stored in highly-pressurized vessels at liquefied state with ambient temperature. If there is by chance a sudden release, it often forms heavy-than-air vapour. The accident release and dispersion of toxic and flammable heavy gas can present a serious ris k to the public’s safety and to the environment. Disease may be caused when the flammable heavy gases are lit. Although great attention has been paid to the hazard of heavy gas dispersion, effective data of filed experiments are still insufficient to make risk assessment and precaution. Through the statistical analysis, draw a conclusion that chemical system in production, transportation and storage process, should first consideration and control of hazardous chemicals, and summarizes the characteristics of the leak diffusion process performance. Subjective factors, equipment inherent defect caused by leakage on China's chemical system is the main reason of the accident. In the process of production, should be strengthen management, strengthen the education of production safety producer. Analysis of the seven factors affecting diffusion of leakage accident, to extract and established the patterns of the leakage accident, and various and leakage accident modes mechanism and the conditions were studied and analyzed. Through the experimental study on practical environment atmosphere that dominated the wind, the wind of gas leakage direction spread concentration distribution, has enormous influence on the spread of gas leakage next wind fastest. Static, as time flies, the space increased concentration of the each point of the trend. In the stable romantic, space the concentration of each point does not change significantly over time, can be considered a steady. Leak gas diffusion next wind fastest, on the site once produce natural gas leak, should be taken into account in the direction and point source leaking the wind direction, wind speed at factors such as timely and accurate prediction leakage, gas may be spread to dangerous area, completes the countermeasures. Key words:Theory of similarity; Leakage model;Leakage diffusion mode;Trace method; heavy gas;Emergency rescue

2011年数学建模-重金属污染模型(特等奖作品)

兰州理工大学数学建模面朝大海 2011高教社杯全国大学生数学建模竞赛 承诺书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮 件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问 题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他 公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正 文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反 竞赛规则的行为,我们将受到严肃处理。 我们参赛选择的题号是(从A/B/C/D中选择一项填写): 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名):兰州理工大学 参赛队员 (打印并签名) :1. 2. 3. 指导教师或指导教师组负责人 (打印并签名): 日期: 2011 年 9 月 11 日赛区评阅编号(由赛区组委会评阅前进行编号):

2011高教社杯全国大学生数学建模竞赛 编号专用页 赛区评阅编号(由赛区组委会评阅前进行编号): 全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):

城市表层土壤重金属污染分析 摘要 随着城市经济的快速发展和城市人口的不断增加,人类活动对城市环境质量的影响日显突出。本文根据大气中重金属污染物的传播特征,将高斯扩散模型应用于土壤重金属污染问题,建立数学模型并进行求解。 问题一中,利用题目所提供的重金属元素样本值,借助matlab软件,形象地绘制出了城市地形图以及八种重金属在该城市五个功能区的二维等高线分布图。对所收集的土壤样品分别采用单因子污染指数法与内梅罗(N.L.Nemerow)综合污染指数法对该城市分区进行污染程度评价。 问题二中,通过对所给数据的分析,利用主成分分析法,推断出了重金属污染的主要原因。 问题三中,通过数据分析,发现了金属污染物的传播特征,由此建立重金属高斯扩散模型。利用数形结合的方法最终确定了污染源的位置。 问题四中,通过对建模过程中考虑因素的全面性以及所忽略因素的得当性分析,评估了所建模型的优缺点,确定了更好地研究城市地质环境的演变模式所需要搜集的信息,并给出了更合理建立数学模型的方法。 关键词:重金属污染主成分分析法单因子指数法内梅罗综合污染指数法高斯模型污染源

点污染源空气污染扩散模型

8 点、中午12 点、晚上9 点都没有排放气体,该怎么算,是不是需要找到一个关于时间t的函数,来计算多长时间之后污染还剩下多少 c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./((sigz+ eps).^2))); 这个函数对吗?该调用什么函数? 问题: 建立单污染源空气污染扩散模型,描述其对周围空气污染的动态影响规律。 现有河北境内某一工厂废气排放烟囱高50m,主要排放物为氮氧化物。早上9 点至下午 3 点期间的排放浓度为406.92mg/m3,排放速度为1200m3 /h;晚上10 点-凌晨4 点期间 的排放浓度为1160mg/m3,排放速度为5700m3 /h;通过你的扩散模型求解该工厂方圆51 公里分别在早上浓度8 点、中午12 点、晚上9 点空气污染分布和空气质量等级。 源代码 clear all clc [x,y]=meshgrid(0:20:5100,0:20:5100); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.^0.865014; sigz=0.0757182*x.^1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./((sigz+ eps).^2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel('C'), clear all clc [x,y]=meshgrid(-5100:20:5100,-5100:20:5100); Q=1836.7; z=1.5; H=50; u=1.7; sigy=0.3914238*x.^0.865014; sigz=0.0757182*x.^1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./((sigz+ eps).^2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel('C'), 分享到: 2015-05-29 16:32 提问者采纳 clear all [x,y]=meshgrid(-51000:100:51000,-51000:100:51000); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.^0.865014;

大气污染扩散模型

第一节大气污染物的扩散 一、湍流与湍流扩散理论 1. 湍流 低层大气中的风向是不断地变化,上下左右出现摆动;同时,风速也是时强时弱,形成迅速的阵风起伏。风的这种强度与方向随时间不规则的变化形成的空气运动称为大气湍流。湍流运动是由无数结构紧密的流体微团——湍涡组成,其特征量的时间与空间分布都具有随机性,但它们的统计平均值仍然遵循一定的规律。大气湍流的流动特征尺度一般取离地面的高度,比流体在管道内流动时要大得多,湍涡的大小及其发展基本不受空间的限制,因此在较小的平均风速下就能有很高的雷诺数,从而达到湍流状态。所以近地层的大气始终处于湍流状态,尤其在大气边界层内,气流受下垫面影响,湍流运动更为剧烈。大气湍流造成流场各部分强烈混合,能使局部的污染气体或微粒迅速扩散。烟团在大气的湍流混合作用下,由湍涡不断把烟气推向周围空气中,同时又将周围的空气卷入烟团,从而形成烟气的快速扩散稀释过程。 烟气在大气中的扩散特征取决于是否存在 湍流以及湍涡的尺度(直径),如图5-7所示。 图5-7(a)为无湍流时,烟团仅仅依靠分子 扩散使烟团长大,烟团的扩散速率非常缓慢, 其扩散速率比湍流扩散小5~6个数量级;图5 -7(b)为烟团在远小于其尺度的湍涡中扩散, 由于烟团边缘受到小湍涡的扰动,逐渐与周边 空气混合而缓慢膨胀,浓度逐渐降低,烟流几乎呈直线向下风运动;图5-7(c)为烟团在与其尺度接近的湍涡中扩散,在湍涡的切入卷出作用下烟团被迅速撕裂,大幅度变形,横截面快速膨胀,因而扩散较快,烟流呈小摆幅曲线向下风运动;图5-7(d)为烟团在远大于其尺度的湍涡中扩散,烟团受大湍涡的卷吸扰动影响较弱,其本身膨胀有限,烟团在大湍涡的夹带下作较大摆幅的蛇形曲线运动。实际上烟云的扩散过程通常不是仅由上述单一情况所完成,因为大气中同时并存的湍涡具有各种不同的尺度。 根据湍流的形成与发展趋势,大气湍流可分为机械湍流和热力湍流两种形式。机械湍流是因地面的摩擦力使风在垂直方向产生速度梯度,或者由于地面障碍物(如山丘、树木与建筑物等)导致风向与风速的突然改变而造成的。热力湍流主要是由于地表受热不均匀,或因大气温度层结不稳定,在垂直方向产生温度梯度而造成的。一般近地面的大气湍流总是机械湍流和热力湍流的共同作用,其发展、结构特征及强弱决定于风速的大小、地面障碍物形成的粗糙度和低层大气的温度层结状况。 2. 湍流扩散与正态分布的基本理论 气体污染物进入大气后,一面随大气整体飘移,同时由于湍流混合,使污染物从高浓度区向低浓度区扩散稀释,其扩散程度取决于大气湍流的强度。大气污染的形成及其危害程度在于有害物质的浓度及其持续时间,大气扩散理论就是用数理方法来模拟各种大气污染源在

扩散模型

2 扩散模型 2.1 高斯模型 燃气泄漏后会在泄漏源附近形成气团,气团在大气中的扩散计算通常采用高斯模型。高斯模型的基本形式是在如下的假设条件下推导出来的[1、9]:假定燃气在扩散的过程中没有沉降、化合、分解及地面吸收的发生;燃气连续均匀地排放;扩散空间的风速、大气稳定度都均匀、稳定;在水平和垂直方向上都服从正态分布。 泄漏燃气相对密度小于或接近1的连续泄漏采用高斯烟羽模型。以泄漏点为原点,风向方向为x轴的空间坐标系中的某一点(x,y,z)处的质量浓度计算公式如下[9]: 平均风速>1m/s时: 平均风速=0.5~1m/s时: 平均风速<0.5m/s时,假设气团围绕泄漏点浓度均匀分布,则距离泄漏点r 处的燃气质量浓度为: 式中ρ d (x,y,z)——扩散燃气在点(x,y,z)处的质量浓度,kg/m3 x、y、z——x、y、z方向上距泄漏点的距离,m u a ——平均风速,m/s δ x 、δ y 、δ z ——x、y、z方向的扩散系数,m

h——泄漏点高度,m ρ (r)——距离泄漏点r处的燃气质量浓度,kg/m3 d r——空间内任意一点到泄漏点的距离,m a、b——扩散系数,m t——静风持续时间,s,取3600的整数倍 扩散系数可查HJ/T 2.2—93《环境影响评价技术导则大气环境》得到。2.2 重气扩散模型 液化石油气密度比空气密度大,属于重气。该类气体泄漏时在重力的作用下会下沉,这时使用高斯模型计算的结果会使泄漏燃气扩散速度偏大,泄漏源附近的浓度偏小。为了解决这个问题,可以引入最早由Van Ulden提出,并由Manju Mohan等发展的箱式模型[1]。箱式模型分为两个阶段:泄漏后的重气扩散阶段和重气效应消失后的被动气体扩散阶段。 重气泄漏后首先是重气扩散阶段。在这个阶段,重气云团由于重力作用逐渐下沉并不断卷吸周围的空气,在卷吸空气的同时,气云受热,最终当重气云团与空气的密度差<0.001kg/m3时,可认为气云转变成中性状态。 随着重气的继续扩散,气云所受的重力不再是影响扩散的主要因素,而大气湍流扩散逐渐占主要地位,这时便是被动气体扩散阶段,可以应用高斯模型计算泄漏燃气的扩散。 3 结论 使用泄漏模型可以计算出燃气泄漏的理论量,此量为扩散计算提供基础数据,可以依据此量分析泄漏后的扩散范围以及预测评价事故后果。使用扩散模型可以对燃气泄漏后的危险区域进行预测。泄漏模型和扩散模型都有各自的适用条件和范围,应该根据泄漏扩散的具体情况分析选择相应模型。

重金属污染

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。 我们参赛选择的题号是(从A/B/C/D中选择一项填写): 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名): 参赛队员(打印并签名) :1. 2. 3. 指导教师或指导教师组负责人(打印并签名): 日期: 2012 年 8 月15 日赛区评阅编号(由赛区组委会评阅前进行编号):

2010高教社杯全国大学生数学建模竞赛 编号专用页 赛区评阅编号(由赛区组委会评阅前进行编号): 全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):

城市表层土壤重金属污染分析 摘要 本文运用Matlab、SPSS及Excel软件,综合分析该城区各区重金属污染程度及其重金属污染的主要原因,并根据重金属的传播特征建立数学模型,寻找出污染源,同时更进一步修改模型以更好研究城市地质环境的演变模式。 针对问题一,利用Matlab作出各重金属浓度空间分布图,可看出重金属浓度大多呈高浓度向外逐步扩散现象。同时,利用单因素污染指数公式及污染等级规则,通过模糊综合评价模型,得出各类区域重金属污染程度如下: 针对问题二,利用SPSS软件对各重金属变量进行R型聚类分析,将重金属污染分为四类:Cr、Ni、Cu为一类,Cd、Pb、Zn为一类,而As、Hg各为一类。并计算出各区域各单金属污染指数进行对照分析,说明该城区重金属污染主要来源于工业活动及其生产品和交通主干道污染。 针对问题三,我们根据重金属物质在土壤中扩散的特点,利用有衰减的扩散模型,构造出浓度与坐标的关系,并利用matlab编程,多次试验,合理筛选数据,通过回归分析的方法,求解参数,从而得到不同重金属的不同污染源。 针对问题四,由于问题三求解的局限性,我们在模型三的基础上,增加了时间、水文、人类活动等因素,得出一个更为具体并符合实际的模型。有利于计算土壤中重金属的浓度。 关键词单因素污染指数模糊综合评价 R型聚类分析有衰减的扩散模型扩散规律

重金属污染源的数学建模

承诺书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。 我们参赛选择的题号是(从A/B/C/D中选择一项填写): A 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名):中国人民解放军国防科学技术大学参赛队员(打印并签名) :1. 李腾骥 2. 袁雪强 3. 赵泉朴 指导教师或指导教师组负责人(打印并签名):杨文强 日期:2011 年9 月9 日赛区评阅编号(由赛区组委会评阅前进行编号):

编号专用页 赛区评阅编号(由赛区组委会评阅前进行编号): 全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):

基于MCMC的城市表层土壤重金属污染分析 摘要:本文针对某城市城区表层土壤重金属污染分析的问题,从各分区重金属元素污染的空间分布、污染程度、污染原因及传播途径等多方面进行了全面的统计与研究,同时采用贝叶斯推理和马尔科夫链蒙特卡洛法(MCMC),利用MATLAB进行编程,对所设计的模型进行求解,得到了各重金属污染源的地理位置坐标,最后对城市地质环境演变模式的问题进行了探究,得出了初步的结论。 首先,根据所给数据信息,利用MATLAB及EXCEL作图,得到了该地区的地势图及各重金属污染程度的分布图。 对于问题1,通过统计各重金属的重污染点在各区域的分布和计算不同区域各中金属污染程度之间的相关系数,得到了各重金属的大致空间分布,并采用改进的内梅罗指数法,得出不同区域重金属的污染程度。 对于问题2,在问题1所得数据及结论的基础上,结合实际情况,对各重金属污染的主要原因进行了较为全面的分析。 对于问题3,通过对污染的传播特征的研究,基于城市的地质特性,得到空气传播为污染主要传播方式的结论,并由此建立关于气体传播的概率模型,采用贝叶斯推理和马尔科夫链蒙特卡洛法,利用MATLAB进行50000次抽样,得到收敛的坐标值,从而计算出各重金属污染源的地理位置坐标分别为:As (2742,7294)、Cd(2382,3693)、Cr(3470,2309)、Cu(2707,2295)、Hg(2708,2294)、Ni(1647,2729)、Pb(2882,3617)、Zn(4152,2299)。 对于问题4,以题中城市土壤重金属污染的演变规律为例进行了一定深度的探究,并对所需信息及模型设计方案进行了初步的讨论。 关键词:内梅罗指数法相关系数贝叶斯推理MCMC法MH算法 一、问题的提出与分析 (一)背景介绍 随着城市经济的快速发展和城市人口的不断增加,人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证,以及如何应用查证获得的海量数据资料开展城市环境质量评价,研究人类活动影响下城市地质环境的演变模式,日益成为人们关注的焦点。 按照功能划分,城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等,分别记为1类区、2类区、……、5类区,不同的区域环境受人类活动影响的程度不同。 现对某城市城区土壤地质环境进行调查。为此,将所考察的城区划分为间距1公里左右的网格子区域,按照每平方公里1个采样点对表层土(0~10 厘米深度)进行取样、编号,并用GPS记录采样点的位置。应用专门仪器测试分析,获得了每个样本所含的多种化学元素的浓度数据。另一方面,按照2公里的间距在那些远离人群及工业活动的自然区取样,将其作为该城区表层土壤中元素的背景值。

大气污染物扩散高斯模型模拟

大气污染物扩散的高斯模型模拟:可视化模拟点源大气污染的扩散Gaussian Atmospheric Dispersion Model 突发性大气污染事故时有发生,对大气污染扩散进行模拟和分析,有利于减小事故的危害,减轻人员伤亡和财产损失。高斯扩散模型是国际原子能机构(IAEA)推荐使用于重气云扩散模拟的数学模型,该模型在非重气云扩散的应用日益广泛。高斯扩散模型是描述大气对有害气体的输移、扩散和稀释作用的物理或数学模型,是进行灾害预测和救援指挥的有力手段之一。 高斯扩散模型 高斯模型又分为高斯烟团模型和高斯烟羽模型。大气污染物泄漏分为瞬时泄漏和连续泄漏,瞬时泄漏是指污染物泄放的时间相对于污染物扩散的时间较短如突发泄漏等的情形,连续泄漏则是指污染物泄放的时间较长的情形。瞬时泄漏采用高斯烟团模型模拟,而连续泄漏采用高斯模型烟羽模型模拟。高斯模型适用于非重气云气体,包括轻气云和中性气云气体。要求气体在扩散过程中,风速均匀稳定。 在高斯烟团模型中,选择风向建立坐标系统,即取泄漏源为坐标原点,x轴指向风向,y轴表示在水平面内与风向垂直的方向,z轴则指向与水平面垂直的方向,具体公式见式: (mg/s); x、y、z轴上的扩散系数,需根据大气稳定度选择参数计算得到(m);x、y、z表示x、y、z上的坐标值(m);u 表示平均风速(m/s);t表示扩散时间(s);H 表示泄漏源的高度(m)。 同理,高斯烟羽模型的表达式如: 技术方法 若用高斯模型算出空间每一个点在一个时刻的污染浓度,这个计算量是很大的。因此所设计的系统一般都是采用先进行图层网格化,由高斯模型计算出有限个网格点的上的污染物浓度,在进行空间内插得到面上每一个点的污染物浓度,并由此得到污染物浓度的等值线。整个过程的示意图如图所示

污染空气的扩散模型

放射性气体扩散的预估模型 摘要:由于放射性气体泄漏造成惨重损失的报道在国际屡见不鲜,近日日本福岛核电站的放射性气体的泄漏事件更让我们关注放射性气体泄漏时在环境中的浓度问题,为了今后事故发生后提供积极的补救措施, 所以对放射性气体的扩散作深入的研究是很有必要的。本文结合高斯烟羽模型、线性拟合,以及微分方程模型,运用MA TLAB软件,分析了泄漏源强度、风速、大气稳定度参数、地面粗糙度参数和计算精确度等的因素对放射性气体扩散的影响,预测了放射性气体浓度在不同时间,不同地区的浓度变化,并且本文模型中的数据可以根据不同的实际情况而加以改变,因而使本文的应用范围大大增加,可以适用于具有较强的应用性。文章首先在第一问中利用MA TLAB软件对数据进行线性拟合,采用微分方程模型得到核电站周边放射性气体在不同地区,不同时间段的浓度变化,得出随着离泄漏源距离的延伸,最终放射性物质的浓度越来越小,趋近于零,即当L趋向无穷是,C(x,y,z,t)趋向于零;当时间趋于无穷时,C(x,y,z,t)也趋于无穷。问题二,问题三中,建立以核电站周边不同地区得距离以及风速为因变量,设置各个主要因素的参考数据,同时,利用高斯烟羽模型对核电站周边地区的浓度进行预测,然后,利用MATLAB软件,将相关数据代入程序,我们得到核电站周边地区的浓度分布的等高曲线。问题四中,通过实际收集数据,集合核电站周边地区的浓度等高曲线,可以直观的看出日本福岛核电站对我国东海岸以及美国西海岸的影响。 一.问题的提出 1.1背景的介绍 目前,核电的发展给国家带来了巨大的经济效益和社会效益,但核电正常运行以及发生泄露时不可避免的会有气载放射性核素排出,这样就给周围的环境产生了一定的影响,因此,正确的测出大气中放射性物质的浓度在环境检测以及安全评估中具有重要意义。 1.2需要解决的问题 的放射性气体以匀速排出,设有一座核电站遇自然灾害发生泄漏,浓度为p 速度为m kg/s,在无风的情况下,匀速在大气中向四周扩散, 速度为s m/s. (1)请你建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。 (2)当风速为k m/s时,给出核电站周边放射性物质浓度的变化情况。 (3)当风速为k m/s时,分别给出上风和下风L公里处,放射性物质浓度的预测模型。

最新城市表层土壤重金属污染分析

城市表层土壤重金属 污染分析

本科生毕业设计(论城市表层土壤重金属污染分析 二级学院:院 专业:数学 年级: 学号:2 作者姓名: 指导教师: 完成日期:2013年5月3日

目录 1 引言 (2) 1.1 问题由来 (2) 1.2 相关信息 (2) 1.3 问题的提出 (2) 2 基本假设 (2) 3 符号说明 (3) 4 问题分析 (4) 5 模型的准备 (4) 6 模型的建立与求解 (5) 6.1 问题1的分析与求解 (5) 6.1.1 模型的建立 (5) 6.1.2 模型的结果 (5) 6.1.3 污染程度的单项污染指数评价 (7) 6.1.4 综合污染程度的评价模型 (9) 6.2 问题2的分析与求解 (10) 6.2.1 模型的建立 (10) 6.2.2 模型的结果 (11) 6.3 问题3的分析与求解 (12) 6.3.1 土壤重金属污染物的来源 (12) 6.3.2 重金属污染物的传播特征 (13) 6.3.3基于matlab的数据处理与三维数据插值模型 (13) 6.4 问题4的分析与求解 (19) 6.4.1 模型的优点 (19) 6.4.2 数据的收集 (20) 6.4.3 模型的建立与求解 (20) 7 模型的推广 (21) 8 结语 (21) 参考文献 (22)

附录 (23)

城市表层土壤重金属污染分析 摘要:本文主要对某城市表层土壤地质环境重金属污染的情况进行分析.首先,利用MATLAB绘制出8种重金属元素在该城区的空间分布,通过内梅罗指数评价法来分析该城区内不同区域重金属的污染程度.接着,通过SPSS,用相关性分析和主成分分析,来分析重金属污染的主要原因.然后,建立三次多项式插值模型和高斯模型,确定污染源的位置.最后,运用灰色预测模型对土壤中的重金属的变化进行预测,较好地研究了城市地质环境的演变模式. 关键词:重金属污染;内梅罗指数评价法;相关性分析;主成分分析;高斯扩散模型;灰色预测模型 中图分类号:X53 Analysis of Heavy Metal Pollution of Urban Topsoil Abstract: This paper analyses the heavy metal pollution situation of urban topsoil geological environment. First,using MATLAB to draw the spatial distribution of eight heavy metals in the city, and using nemerow pollution index method for the analysis of the heavy metal pollution of the city zone in different areas. Then,by using SPSS, pollution causes of heavy metals are analyzed with correlation analysis and principal component analysis. After that, establishing cubic polynomial interpolation model and Gaussian diffusion model to confirm the pollution sources. At last, grey prediction model is set up to forecast the changes of heavy metals in the soil, preferably forecasting the evolution of the geological environment model. Key words: Pollution of heavy metals; Nemerow pollution index method ; Correlation analysis; Principal component analysis; Gaussian diffusion model;Grey prediction model.

SLAB用户手册模拟重气体泄漏的空气扩散模型中文简要

SLAB 用户手册:模拟重气体泄漏的空气扩散模型 中文简要用户使用手册 环境保护部环境工程评估中心 国家环境保护环境影响评价数值模拟重点实验室

手册说明 本用户手册基于《USER’S MANUAL FOR SLAB: AN ATMOSPHERIC DISPERSION MODEL FOR DENSER-THAN-AIR RELEASES》(1990.06)编写,仅对美国EPA网站所提供的模拟重气体泄漏的空气扩散模型SLAB的使用方法提供中文版简要说明,更详细的程序使用说明请查阅相关的软件手册及文档,或采用带图形界面版的商业软件。 本手册由环境保护部环境工程评估中心国家环境保护环境影响评价数值模拟重点实验室负责编写,参与人员包括:易爱华、陈陆霞、胡翠娟、梁昊、杨晔、丁峰等。 本手册版权所有,转载及印刷请与环境保护部环境工程评估中心联系。 本手册所涉及的模型系统及本手册电子版本下载地址:

一、SLAB简介 SLAB是用于模拟重气体泄漏的空气扩散模型。该模型最初基于Zeman于1982年提出的关于重气体云的空气卷吸和重力扩散的理念而开发。SLAB早期相关工作由美国能源署支持。SLAB的进一步开发由USAF工程和服务中心(1986年开始)和美国石油学会(1987开始)共同提供支持。 现行的SLAB版本可以模拟连续的、限时的和瞬时的物质泄漏,泄漏源包括以下4种:地面液池蒸发、高于地面的水平射流、烟囱或高于地面的垂直射流,以及瞬时释放。 SLAB除可以用于模拟重气体的扩散,还可以模拟中性浮力气体的烟云扩散,以及烟云轻于空气时的上升过程。 泄漏时的空气扩散过程可以通过求解质量、动量、能力和物质的守恒方程来计算,如图1所示。为了简化守恒方程的求解过程,方程可以通过将烟云作为稳态烟羽或瞬时烟团在空间上进行平均。连续排放(持续时间非常长的排放源)可以作为稳态烟羽。有限时间的排放采用稳态烟羽模式描述最初烟云的扩散,而且在该排放源持续泄漏的时间段内,可以一直使用稳态烟羽模式。释放一旦终止,烟团被视为瞬时烟团,之后的扩散采用瞬时烟团模式来计算。对于瞬时泄漏的排放源,整个过程都均使用瞬时烟团扩散模式。 二、理论介绍 2.1重气体扩散模型简介 重气体泄漏的空气扩散模型受到多种物理现象的影响,这些物理现象在中性或浮力气体泄漏中可能不会发生或者即便发生也不是很重要。这些物理现象包括:重气体烟云的稳定密度分层导致的湍流衰减;由于重力流和初始排放源动量导致的环境速率场的改变;由于液滴形成和挥发以及在过热或低温液体排放情况下的地面加热对烟云温度、浮力和湍流的热力学效应;此外,我们所关注的某种特定的重气体的浓度可能和典型大气污染物关注的累积浓度差别很大。例入,对于易燃气体,关注的是瞬间浓度;而对于有毒气体,关注的则是几分钟到几小时的浓度,以及累积浓度。因此,为了能够更好的预测出重气体泄漏时有毒浓度区的大小和持续时间,所有重要的物力现象都需要进行考虑,而且预测过程中要使用最合理的浓度平均时间。 为了满足重气体泄漏情形的要求,SLAB模型以质量、动量、能量和物质守恒方程的平均形式为起点,在该理论框架的基础上进行开发(如图1所示)。这些方程用于计算扩散烟团的空间平均性质,并且以两种方式来代表两种不同的扩散模式:稳态烟羽扩散模式和瞬时烟团扩散模式。

重金属污染物扩散FICK 模型

FICK 模型: <一>,Fick 第二定律简介: 22x C D T C ??=?? C T C =),0( 0)0,(C x C = 0),(C T C =∞ 符号说明: ),(T x C : T 时间后距离扩散源x 米处气体浓度。 0C : 初始浓度。 C : 气体浓度。 D : 扩散系数。 <二>,模型假设: (1)假设扩散系数为常数。 (2)扩散过程中没有特殊外界条件影响(刮风或下雨等)。 (3)扩散过程中外界温度恒定。 (4)扩散过程中气体没有与外界发生化学反应。 <三>,模型求解: 假设气体发生泄漏(T=0)时刻,所以气体泄漏中心周围处浓度C 与时间和距离(r )的关系可用Fick 第二定律求解,公式如下: 22r C D T C ??=?? (1) 令T r =λ代入(1)式可得: T C D r C D 1.222λ??=?? (2) 可得出(1)式为: λλλ .222??=??-C C D (3) 若n=2,D 41=α,代入化简,积分可得:

?+-=λλλ02)4exp(B d D A C (4) 令DT r D 22==λ β,由高斯误差积分公式可求得: ? ∞=-022)exp(πββd (5) 根据初始化条件可知,02C A π- =,0C B = 最终求得的模型公式为: )2()exp(2000200Dt r erf C C d C C C -=--=?β ββπ (6) 不妨取D=0.00001,发生泄漏时初始气体浓度1000=C 量纲,扩散时间T=1.0*10 8秒,利用Matlab 编程模拟,得到模拟图像如下: 模型评价: Fick 第二定律模拟出了理想状态下气体源扩散的过程,从仿真图(地平面点浓度)可以直观的看出,在气体源位置浓度最大,然后均匀向四周浓度逐步减少扩散,大致符合实际气体扩散过程,模型缺点为,条件过于苛刻,现实气体扩散中,往往收到外界条件影响,比如风向等等。

气体扩散模型

放射气体模型的预估模型 摘要 本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的。且结合了高斯烟羽模型、线性拟合,以及微分方程模型,运用MATLAB软件,分析泄漏源强度、风速、大气稳定度参数等因素对放射性气体扩散的影响,预测了放射性气体浓度在不同时间,不同地区的浓度变化,并且本文模型中数据可以根据不同的实际情况而加以改变,因而是本文的应用范围大大增加,可以适用于具有较强的应用型。 对于问题一,讨论在无风的情况下,放射性气体以s m/s的匀速在大气中向四周扩散。本问中由于不考虑风力的影响,且扩散出来的气体匀速向四周散开,这样经过任意时刻t,扩散的气体围成一个半径为st的球,且距球心位置不同的地方浓度值不同。采用列数列的表现方法,设定相同时间段t,把条件进行整理,并经过简单计算得出每段时间所预测得到的扩散距离r和浓度C。利用MATLAB 软件对数据进行线性拟合,采用微分方程模型得到核电站周边放射性气体在不同地区,不同时间段的浓度变化,得出随着离泄漏源距离的延伸,最后放射性物质的浓度越来越小,趋近于零,即当x趋向无穷时,C(x,y,z,t)趋向于零;当时间趋向于无穷时,C(x,y,z,t)也趋于无穷。 对于问题二,要探究风速对放射性物质浓度分布的影响。风速的处理是问题的中心,采用大气污染的经典高斯扩散模型,实现了高斯烟团气体扩散模型的动态预测,分析计算了气体扩散过程中的各关键参数。 对于问题三,本文在问题二的基础上,结合考虑风速和放射性物质扩散速度在空间中的矢量运算,将在上风和在下风不同情况下与传播速度s之间的比较的分析,利用高斯烟羽模型对核电站周边地区的浓度进行预测,然后,利用MATLAB软件,将相关数据代入程序,我们得到核电站周边地区的浓度分布的等高曲线。 对于问题四,本文参阅整理大量气象、地理、新闻资料,选择我国东海岸典型地域---山东半岛作为研究对象,综合考虑对应海域平均风速及风向、地理距离、海水对放射性物质扩散的部分反射系数等因素,集合核电站周边的浓度等高线,可。 关键词:放射性气体扩散浓度变化高斯修正模型预测

城市表层土壤重金属污染分析模型数学建模精品

城市表层土壤重金属污染分析模型数学建模精品

城市表层土壤重金属污染分析模型

高教社杯全国大学生数学建模竞赛 承诺书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮 件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的 问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他 公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正 文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反 竞赛规则的行为,我们将受到严肃处理。 我们参赛选择的题号是(从A/B/C/D中选择一项填写): A 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名):楚雄师范学院 参赛队员 (打印并签名) :1. 马银梅 2. 王华 3. 邓小娇 指导教师或指导教师组负责人 (打印并签名): 日期:年月 日

2011高教社杯全国大学生数学建模竞赛 编号专用页 赛区评阅编号(由赛区组委会评阅前进行编号): 全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):

城市表层土壤重金属污染分析模型 【摘要】: 随着城市经济的快速发展,城市人口也不断增加,工业化进程日趋加快,重金属污染日益严重,因此人们也越来越关注城市环境质量问题。 为剖析城市城区环境质量问题,对某市表层土壤重金属污染进行分析。根据城市土壤中重金属元素(As,Cd,Cr,Cu,Hg,Ni,Pb,Zn)的含量数据及背景值进行分段,并应用Spss及Matlab软件得到金属元素的空间分布。结合图形并用主成分分析法分析出该城区内不同区域重金属的污染程度,其传播特征是。并且明确了污染源的位置。通过模型分析明确了该地区的主要污染是:其主要来源于。研究结果充分表明:在各个不同的功能区内主要的重金属污染都有可能是不同的重金属,因此我们对不同区域不同的主要重金属污染进行分析。最后分析了模型的优缺点,指出为了更好地研究城市地质环境的演变模式,还应收集当地主要工业种类,城市交通图,土壤分布图,大气质量监测数据,气象数据(降水、风),水文观测数据等信息。 【关键词】 曲线拟合、空间分布、传播特征、主成分分析法、

相关文档
最新文档