基于GIS的海上泄漏化学品行为数值模拟研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
随着海上化学品运输、转存量的不断增加,化学品泄漏风险也日益增加,对海上泄漏化学品行为进行合理、可信的数值模拟并将其结果直观地显示是目前应急辅助决策和风险评估的关键技术。因此利用地理信息系统(GIS)强大的空间数据分析显示功能,进行基于GIS的海上溢漏化学品行为数值模拟研究具有重要的理论意义和实用价值。
     在分析了国内外潮流模拟、近海水域环境污染扩散数值模拟以及GIS在海洋科学等领域的研究应用现状基础上,根据海上溢漏化学品的行为特点,采用Monte—Carlo法分别建立了不同类型的化学品在大气中挥发扩散、水面漂移扩展和水下溶解扩散的数值模拟模型,并通过对海气相互作用的研究,确立海洋—大气耦合模型,从而能够全方位地计算模拟溢漏化学品在海洋和大气中的行为动态。
     为了给数值模拟过程迅速提供必要的流场和风场动力条件,本文分别建立了近海多开边界三维潮流快速预报模块和风场插值预报模块。首先应用伴随算子法从岸边验潮站资料估计外海开边界条件,然后用ADI法进行三维潮流计算,并用干湿法处理动边界,最终将计算结果以调和常数形式存储成图层文件用以预报潮流场;风场预报是采用插值法将监测站位风场资料插值生成海面10米风场,外推到垂向各层的风场,并诊断形成三维风场。
     为提高数值模拟的应用性,本文采用GIS组件技术(MapX)和面向对象开发环境(Visual Basic),将上述数值计算模型与地理信息系统进行了嵌入式紧密集成,初步形成了应用软件系统。实现了GIS支持下的坐标和投影变换、计算网格划分、水深和栖息地数据库的调用、等值线绘制等功能,以及对模拟计算结果时空分布的可视化表达。最后通过实际案例应用证实了系统的准确性和实用性。
With the increasing transportation of liquid chemicals in bulk by ships, the risk of the chemicals spill is getting higher and higher. So the reasonable numerical modeling of the spilled chemicals at sea is a key technique for the risk assessment and contingency decision-making. By taking advantage of the GIS's capacities of analysis and-visual expression on spatial data, the paper develops the numerical simulation system to calculate and display the fates and behaviors of spilled chemicals at sea.
    Based on the analysis of the State-of-the-art numerical modeling on the tidal current, spreading of pollutant in the water column, and the application of the GIS in the marine science, the Monte-Carlo method is used to established the individual numerical model of emission in the air, advection on the water surface and dilution in the water column for different type of chemicals according to the chemical and physical character of each type of spilled chemicals. And a water-air interactive model is studied to consider the physical transfer between the water and air. So the fate and behaviors of the spilled chemicals at the sea and in the air can be simulated simultaneously.
    To provide the prevail conditions of the tidal current and the wind, the exoteric boundary condition using the data from exoteric station is established by conciliatory method, in which the movable boundary is dealt with by "dry-wet" method. Then the simulation models are computed by ADI method. The unison constants are saved in the form of layer databases for the aim to predict tidal current quickly. The wind field at the
    height of 10 meters can be obtained by inserting the observatory data, and the wind field at the other height can be extrapolated using the diagnostic method.
    To promote the application of the numerical models, the GIS component (MapX)
    
    
    
    and an Object-oriented Programming Environment (Visual Basic ) are used to integrate the numerical models and GIS as a software system. With the help of the GIS module, the conversions of projection and coordinate, the formation of calculation grid and isoclines, and other functions can be approved. And the simulation results of the fate and behaviors of the spilled chemicals in space and time can be visually expressed.
    In the end, the numerical modeling system is applied to simulate the "3.24" oil spill accident happened in the Pearl River Estuary. The correspondence between the simulation results and the fact indicates that the numerical models and the system are feasible to simulate the behaviors of the spilled chemicals.
引文
[1] 《1990国际油污防备、反应和合作公约》(OPRC公约)
    [2] 《1996年国际海上运输有害有毒物质损害责任和赔偿公约》(HNS公约)
    [3] 《国际散装运输危险化学品船舶构造和设备规则》(IBC规则)
    [4] 《经1978年议定书修订的1973年国际防止船舶造成污染公约》(MARPOL73/78公约)
    [5] 《中华人民共和国海洋环境保护法》
    [6] 杜海涛,大连湾溶解、轻质、保守液体化学品溢漏后的归宿,大连海事大学,硕士毕业论文,2000年
    [7] 朱国林,大连湾易溶保守液体化学品污染扩散的数值计算,大连海事大学,硕士毕业论文,1998年
    [8] IMO MEPC, MANUAL ON CHEMICAL POLLUTION
    [9] 刘斌、李凯,化学品船安全知识与操作,大连海事大学出版社,1999年
    [10] 交通部海事局,沿海船舶载运有毒液体物质污染事故应急对策研究,2002年
    [11] 国家环境保护总局监督管理司,环境影响评价教材,海洋环境影响评价
    [12] 江文胜,孙文心,渤海悬浮颗粒物的三维输运模式,《海洋与湖沼》,2001年第31卷第6期
    [13] 傅国伟编著.河流水质数学模型及其模拟计算.北京:中国环境出版社,1987
    [14] Lin P.N,Shen H W.2_D flow with sediment by characteristies method.Hyer Engrg ASCE. 1984,110(5)615-626
    [15] 赵隶华.戚晨,庚维德,徐稼华,裴中平.平面二维水流一水质有限体积法及黎曼近似解模型。水环境科学进展:Vol11(4)368—374
    [16] O.B.Anderson,P.L.Woodworth, R.A.Flather, Interecomprision of recent ocean tide models,J.Geophys,Res,1998.13(8):261-282
    [17] 何少苓,龚振瀛等,隐式破开算子法在二维潮流计算中的应用,海洋学报,1985,
    
    7(2):225—232
    [18] 李孟国.潮州港三维潮流模型研究.海洋通报:1995,4(6)1
    [19] 江毓武,陈宗团。厦门港三维潮流数值模拟。海洋工程:1998,16(4)66—72
    [20] 雷孝平,朱国林,陆凤梧。大连湾潮流和潮位的数值计算。大连海事大学学报:1998:24(3)40-44
    [21] 陈满春。厦门港潮流场的数值模拟。海洋通报,1997:16(3)13-20
    [22] 陈玉霞.厦门西港潮汐潮流的数值计算。台湾海峡,1985:4(2)111-118
    [23] 江毓武,张潞平,陈宗团。闽江口水动力和污染扩散数值模拟.台湾海峡:1997:16(2)169-173
    [24] 吴巍等。渤海局部海域风暴潮漫滩计算模式。海洋与湖沼,1995,27(2)149
    [25] 温生辉,陈季良。厦门附近海域潮流场的数值模拟.海洋学报,1996:18(2)15-25
    [26] 王泽良,陶建华。缓坡潮间带海湾的污染物输移扩散数值模拟方法研究。海洋学报,1998:20(3)120-127
    [27] 王泽良,陶建华。渤海湾流场及污染物分布的数值模拟研究。海洋与湖沼,1999:30(2)224-230
    [28] 孙琪,孙效功,李瑞杰。窄缝法在河口区悬沙输运数值计算中的应用.海洋与湖沼.2001:32(3)298-30]
    [29] Lynch and Gray.Finite element simulation of fluid in deforming region. J.Compute. Phys. 1980,39,135-153
    [30] Johns B.SK Dube,P C Sinha, U C Mohanty and A D Rao.The simulation of a continuously deforming lateral boundary to problems involving the shallow water equation.Comput,Fluids,1982,10,105-116
    [31] 史峰岩,林文心,巍更生。运动侧边界海洋问题的自适应网格模拟方法。海洋学报.1997,19(3)1-9
    [32] Takisawa A, S Koshizuka and S Kondo. Generalization of physical component boundary fitted coordinate(PCBFC) method for the analysis of free-surface flow. Int.J.Numer;method Fluids; 1992,15:1213-1237
    [33] 叶安乐,李凤歧编著。物理海洋学。青岛海洋大学出版社,1992
    
    
    [34] 魏皓、武建平、张平。海洋湍流模式应用研究。青岛海洋大学学报。2001:31(1)7—12.
    [35] 罗家海,潘南明,汪道明.汕头港湾附近水域潮流特征和污染物扩散的数值计算.海洋环境科学.1997,16(2):46-50.
    [36] 管卫兵,王丽娅,潘建明等。POM模式在河口湾污染物质输运过程模拟中的应用.海洋学报.2002,24(3):9-17.
    [37] 陈永平,刘家驹,喻国华.潮流数值模拟中紊动粘性系数的研究。河海大学学报.2002,30(1);39-43.
    [38] 黄思训,伍荣生编著.大气科学中的数学物理问题.气象出版社,北京.2001
    [39] 宇德明.易燃、易爆、有毒危险品储运过程定量风险评价.中国铁道出版社,北京.2000
    [40] 王宏模.大气扩散模式及其应用.包钢科技,2000,26(4):81-83
    [41] Anderson, E. L...The OILMAPWin/WOSM oil spill model: Application to hindcast a river spill. In Proceeding of the 18th Arctic and Marine Oil Spill Program, Technical Seminar, Edmonton, Alberta, Canada: 1995. 793—817
    [42] Leech, M., et al.., OSIS: A Windows Oil Spill Information System, In Proceeding of the 16th Arctic and Marine Oil Spill Program, Technical Seminar, Calgary, Alberta, Canada.,1983,5(1):27-30
    [43] Reed, M., N. Ekrol, H. Rye, and L. Turner. Oil Spill Contingency and Response(OSCAR) Analysis in Support of Environmental Impact Offshore Namibia. Spill Science and Technology Bulletin, 1999, 5(1): 29--38.
    [44] 娄安刚,黄祖珂.海面溢油轨迹的分析与预报.青岛海洋大学学报,1994,24(4):477—454.
    [45] 张波.中文Windows环境下的海上溢油预报系统。海洋环境科学,1997,16(1):
    
    37—41
    [46] 刘桦,何有声.河口三维流动数学模型研究进展.海洋工程,1999,18(2):87-93
    [47] 李孟国,曹祖德.海岸河口潮流数值模拟的研究与进展.海洋学报,1999,21(1):111-126
    [48] 吴冠,张波,张砚峰等.溢油应急系统中的潮流快速预报,海洋环境科学,1997,16(1):42-44
    [49] Backhaus, J. O..A semi-implicit scheme for the shallow water equations for application to shelf sea modeling. Continental Shelf Research, 1983, 4(2): 243-254
    [50] Backhaus, J. O..On the atmospherically induced variability of the circulation of the northwest European shelf sea and related phenomena—a model experiment. Modeling Marine System, 1990, 1: 94-134
    [51] 吕咸青,方国洪.渤海开边界潮汐的伴随法反演.海洋与湖沼,2002,33(2):113-120
    [52] 朱江,曾庆存,郭冬建等.利用伴随算子法从岸边潮位站资料估计近岸模式的开边界条件.中国科学(D辑),1997,27(5):462—468
    [53] 娄安刚,王学昌,于宜法等.蒙特卡罗方法在海洋溢油扩展预测中的应用研究.海洋科学,2000,24(5):7-10
    [54] 黄立文,王献孚,俞齐清等.多开边界群岛海域溢油模拟的水动力模式.武汉交通科技大学学报,2000,24(4):367-371
    [55] 严世强,熊德琪.潮流数值模拟研究进展.浙江海洋学院学报.2001,21(3):263-268.
    [56] 雷孝恩.大气污染扩散的数值模拟,1999.北京:气象出版社.105
    [57] 桑建国,温市耕.大气扩散的数值计算,1990.北京:气象出版社.5
    [58] Mcqueen J T, Draxler R R, Rolph G D,. Influence of grid size and Terrain Resolution on Wind Fields Prediction from an Operation Mesoscale Model. J .Appl.
    
    Met.1993.34:2166
    [59] 程麟生.中尺度大气数值模式和模拟,1994.北京:气象出版社.102~110
    [60] 桑建国主编.北京大学大气环境模式用户手册.北京:北京大学地球物理系,北京大学环境科学中心,992
    [61] 谷清,李云生编著.大气环境模式计算方法.2002:北京,气象出版社,5-14
    [62] 雷孝恩,张美根,韩志伟编著.大气污染数值预报基础和模式,1998:北京,气象出版社,59-66
    [63] 钟中,苏炳凯,赵鸣.大气数值模式中有效粗糙度计算的一种新方法.自然科学进展.2002,12(5):519-523
    [64] 杨洪斌,马雁军.辽东湾地区风温廓线及摩擦速度研究.辽宁气象,2002(4):17-18
    [65] 余兴,刘贵华,李兆元.三维风场调整和诊断模式的实用研究.陕西气象,1997(1):15-17
    [66] 黄鹏,施宗城,陈伟等.大气边界层风场模拟及测试技术的研究.同济大学学报,2001,29(1):40-44
    [67] 农尚尧,徐银薪.层结性对大气风场的影响.大气科学,1994,14(1):37-45
    [68] 康凌,陈家宜,蔡旭晖.核事故应急条件下风场的实时预报.辐射防护,2003,23(4):193-209
    [69] 廖木星.海面风场预报的技术研究.中国航海.2003,56(3):65-67
    [70] 冯凯,王嵘.边界层风场调整的分析.南京气象学院学报.1994,17(4):508-511
    [71] 余琦,刘原中.复杂地形上的风场内插方法.辐射防护.2001,21(4):213-218
    [72] 单铁良,董羽仑,黄金萍.卡尔曼滤波方法预报因子的选择.河南气象,2001,1:16-17
    [73] 刘敦训,金伟福.卡尔曼滤波方法在风的短期预报中的应用.山东气象,1998,18(1):29-32
    
    
    [74] 孙心文,江文胜,李磊编著.近海环境流体动力学数值模型.2003,北京,科学出版社,396-402
    [75] 王红梅,朱振海.海洋地理信息系统国内外研究进展.遥感技术与应用,1999,14(3):49-55.
    [76] 何明元,石汉青,孙立明.地图投影在气象卫星资料处理中的应用.气象仪器设备,2002,(1):18-20.
    [77] 钱辉,董火根,陈楚铭.极射赤平投影中的数学方法.高校地质学报,1997,3(3):328-338.
    [78] 苏志刚.雷达视频地图的坐标变换.中国民航学院学报,1999,17(1):16-18.
    [79] 林建国,殷佩海,徐洪磊.预测海上溢油动态的快速算法.航海技术,船舶机电专业委员会论文集专刊.2000:13—23
    [80] 严世强,熊德琪,滕广伟等.地图投影及其在基于GIS的海洋数值预报中的应用.浙江海洋学院学报(自然科学版),2003,22(4):321-326
    [81] 赵冬泉,贾海峰.基于Mapx的等值线绘制分析.清华大学学报(自然科学版),2002,42(8):1023-1025
    [82] Mackay D. et al. Oil spill process and model. Environmental Impact Control Directorate, Canada, Report, 1980: EE-8.
    [83] Mackay D. and Tebeau P. A.: Development and calibration of an oil sill behavior model, Report No' CG-D. 1982:27-83.
    [84] G. A. Robilliard. Evaluation of Compensation Formulate to Measure Natural Resource Damages, 1993 Oil Spill Conference Proceedings.
    [85] T. A. Grigalunas. Marine and Shelf Resources Damage Assessment Model, 1985 Oil Spill Conference Proceedings.
    
    
    [86] 苏奋振 周成虎等,海洋渔业地理信息系统的发展、应用与前景。水产学报,2002,26(2):169-174
    [87] Trimanes J A, Cotos J M, Tobar A, et al.A Geographic information System for Operational Use in Pelagic Fisheries-FIS.IEEE Oceans'94,1994,3:532-535
    [98] Haddad, Kenneth D, MacAulay, McGarry G, et al. GIS and Fisheries Management. Coastal Zone'93 Coastal Zone: Proceedings of the Symposium on Coastal and Ocean Management,1993,1:68-78
    [89] 邵全琴,周成虎,张明金等.海洋渔业电子地图系统软件设计与实现.水产学报 2001,25(4):367-372
    [90] Rongxing Li, Saxena N K. Maxine Geographic Information System: Approaches Applications and Trends.IEEE Oceans' 93,1993,3:49-54
    [91] Basu A. Development of Marine Information System for Various Data Analysis and Data Integration in the Hawrtisan Exclusive Economic Zone. IEEE Ocean'95,1995,1:146-153
    [92] 赵玲,赵冬至.基于GIS的海域环境质量评价模型研究.遥感技术与应用,1998,13(3):61-65
    [93] 罗云启,罗毅.数据化地理信息系统Mapinfo应用大全,1版,2001.北京希望电子出版社,2-20
    [94] 郑江玲.地理信息系统(GIS)软件技术,四川测绘,2001,25(1):18-20
    [95] 林孝松.基于GIS的重庆地质灾害信息管理系统.中国地质灾害与防治报.2001,12(3):74-80
    [96] 王风林,毕彤.沈阳突发性环境污染事故应急地理信息系统.环境保护科学,2000,26(2):27-32
    
    
    [97] 庄学强,王炜,熊德琪.船运液体化学品泄漏污染大气应急地理信息系统.中国航海.2003,1:23-26
    [98] 庄学强,王炜,熊德琪.基于GIS的船运液体化学品泄漏污染大气扩散模拟系统.交通环保.2003,1:14-18
    [99] 张存智.黄海北部海域三维潮流数值模型.海洋预报.2000,17(1):1-13
    [100] 徐钟济.蒙特卡罗方法.1985上海:上海科学技术出版社,.
    [101] 常清英.扩散方程的蒙特卡罗法误差估计.中南民族学院学报(自然科学版).1994,13(1):72-75
    [102] 卢贵武,李荣.蒙特卡罗计算机模拟方法及其在流体理论研究中的应用.石油大学学报(自然科学版).1999,23(3):112-116
    [103] 吉庆丰,刘超,郑邦民.计算复杂边界渗流的蒙特卡罗方法.计算力学学报.2000,17(4):390-393
    [104] 吉庆丰,郑邦民.不规则游动网格的蒙特卡罗方法.水动力学研究与进展.2000,15(2):177-181
    [105] 吉庆丰,郑邦民.抛物型方程的蒙特卡罗解法.江苏农学院学报,1992,(3).
    [106] 吉庆丰,郑邦民.蒙特卡罗解法解堰闸流动.水利学报,1998,(9).
    [107] 徐希孺,王子荣.用蒙特卡罗方法计算大气点扩散函数.遥感学报.1999,3(4):268-278
    [108] 田瑞明.Monte-Carlo法模拟复杂地形对扩散的影响.大气科学.1994,18(1):37-42
    [109] 姚仁太.排放高度对建筑物尾流扩散影响的分析与蒙特卡罗模拟.辐射防护 2001,21(4):209-214
    [110] 姚仁太,夏益华,赵宏圣.铀矿环境质量评价中进风井对大气扩散影响的随机行走模拟.辐射防护.1999,19(5):354-360
    
    
    [111] 李军.地球空间数据集成的基础理论与应用研究.中国科学院地理研究所博士论文,1998.
    [112] 张健挺.地理信息系统集成若干问题探讨.遥感信息,1998,(1):14-18.
    [113] 张犁.GIS系统集成的理论与实践.地理学报,1996,51(4):306-313.
    [114] 张明金.核心地理信息系统的设计研究.地理学报,1995,50(supplement):44-53.
    [115] Goodchild M F, Haining R P, Wise, S, ct al. Integrating GIS and spatial analysis: Problems and possibilities. International Journal of Geographical Information Systems,1992,6(5): 407-423.
    [116] Yoon, J. Watershed-Scale Nonpoint Source Pollution Management Based on Spatiotcmporal Pararnctcr Model and GIS Linkage. In: Pascolo P, Brcbbia C A. GIS Technologies and their Envionmcntal Applications. Southamton, UK: Computational Mechanics Publication, 1998: 3-12.
    [117] Wu F L. SimLand: A prototype to simulate land conversion through the integrated GIS and CA with AHP-derived transition rules. International Journal of Geographical Information Systems, 1998,12(1):363-82.
    [118] Todini E, Marsigl M, Zamboni L. The Role of GIS in the Decision Support System OADSSEI(Open Architecture Decision Support System for Environmental Impact Assessment). In: Pascolo P, Brebbia C A. GIS Technologies and their Envionmental Applications. Southamton, UK: Computational Mechanics Publication, 1998: 83-98.
    [119] 王晓蓉.环境化学,第一版.南京大学出版社,南京,1993:103
    [120] 庄学强.基于GIS的水运气体化学品泄漏大气扩散系统研究.大连海事大学硕士学位论文,2003

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700