河流检测matlab,基于MATLAB的河流突发性污染事故应急监测布点研究

基于MATLAB的河流突发性污染事故应急监测布点研究

2 0 1 0年 7月 第 2 7卷第 7期 海洋开发与管理 J u l y 2 0 1 0 V0 1 ．2 7 NO． 7 基于 MA T L A B的河流突发性污染事故 应 急监测布点研 究 胡 越 ，张 戈 ( 辽宁师范大学城市与环境学 院大连 1 1 6 0 2 9 ) 摘 要 ：文章研究了MATL AB在突发性水污染事故应急监测布 点中的技术支持。 以松花 江污 染 事故为例 ，用 MATL A B语 言编 写程序 ，并在 MATL AB上实现地 表水 环境 系 统模 拟 ，根 据预 测结 果结合 布 点原 则布 设 合 理 的监 测 点 ，为 污 染 事 故 的处 理 工作 提 供 科 学 依 据 。 关 键 词 ：突 发性 环境 污染事 故 ；应急监 测 ；MATL AB语 言 近 年 来 ，我 国伴 有 严 重 安 全 与 环 境 后 果 的 重 大环境 污 染 事 故 和公 害 事 件 数 量 呈 持 续 上 升 趋势，直接经济损失 每年数亿元l 】 ] 。突发性环 境 污染事 故 是 指社 会 生 产 和 生 活 中使 用 的危 险 品在其产生 、运输 、使用 、储存 和处置的整个 生命周期中，由于 自然灾 害、人为 的疏忽或错 误 操作 ，在 瞬 间或短 时 间 内 大量 排 放 污染 物 质 ， 对环境造成严重污染和破坏 ，给人 民的生命和 国家 财 产 造 成 重 大 损 失 的恶 性 事 故[ 2 ] 。 因 此 ， 对 突发 性 环 境 污 染 事 故 的控 制 ， 已经 到 了刻 不 容缓 的地步 ，应 引 起 我 国 政 府 的 高度 重 视 ，同 时环境 污染 事故 应 急监 测 也 作 为一 个新 的课 题 ， 摆在 了环境 工作 者 的面前 。 1 基 于 MA T L A B地 表 水环 境 系统 模 拟 模型简介 MATL AB诞 生 于 2 0世 纪 7 0年 代 ，它 的 编 写者是 C l e v e Mo l e r 博士和他的同事。当时，他 们 利 用 F o r t r a n开 发 了 两 个 子 程 序 库— — E I s — P A C K 和 L I NP AC K。这 两个子 程序库 是 求解 先 行方 程 的程 序 库 。C l e v e Mo l e给 E I S P AC K 和 L I NP AC K的接 口程序取名为 MA TL AB ，意为 矩 阵( Ma t r i x ) 和 实 验 室( L a b o r a t o r y) 的 组合 。 对于一段长河流而言 ，在这里我们只研究其 纵 向扩 散忽略横 向扩 散 。因此 ，一维 模拟 模 型更 加适用 ，更能反 映出污染物在河 流中的变化情况 。 由于突发性污染 事故 中污染 源均 是 瞬间排 放 ，应 用 瞬时点源模 型 比较 合适 _ d ] 。下 面介绍 两种 适用 于河 流突发性污染事 故 的一维扩散 模型[ 5 ] 。 1 ． 1静 止 水体 瞬 时点 源模 型 该模型用于描述 瞬时释放的一个点污染源 到静态水体 中的污染 物扩散规律，此 时污染物 浓度随时间 t 和空间位置3 2发生变化 ] 。 一D ～ 丝 2 Dt Dx一 ( 一 o o 0 ) 初始条件： c： ( ， ￡ ) J 。一 C o ( ) ( 一 o o 0) ( 1 ) 式 中 ： ( ) 一 j z ∽ ； C o 为 瞬 时 释 放 式 中 ： ‘ 一 { o ( z≠ o) 为 瞬时释 放 的污染 源 在 单 位 面 积 上 形 成 的浓 度 ，可 表 示 为 C o =M／ A，M 和A 分 别为 瞬 时释放 质量为 M 的 污 染物 在单位 深度 的 面积 A 上 。 应 用 F o u r i e r 变 换 ，可 以得 到该 模 型 的解 析 解 为 C( x, t )= M e x pt ( 一 1I ( 2 ) 2 A D ＼ 4 z ￡ 式中：D 表示 方向的扩散系数，I n 。s 。 1 ． 2 有水 流瞬时点源 的模型 对于河流而言，水体经 常有一个主导流 向， 如果在 一0的地 方 瞬 时投 放 污染 源 ，在 投 放点 第 7期 胡越 ，等 ：基 于 MAT I AB的河流突发性污染事故应急监测布点研究4 5 形成 浓度 C 。 ，则 描 述 其 下 游 污 染 物 浓 度 时 空 变 化 的数学 模型 可 以表述 为 0 C — D OC — KC ( 1z 0， t 0 ) 初始条件： C ( x ， ) J 一 一C o ( z ) ( 3 ) ( z≥ 0 ) 边 界条 件 ： l i mC( x， ￡ )一 0 ( t O) 式 中 ：K 为污染 物 降解速 率 常数 ，S ～ 。 同样应 用 F o u r i e r 变 换 ，可 以得 到该 模 型 的 解析 解为 c c ， 一 e [ 一 ] e x p ( -- K t 一 exp[ 一 ] e x p ( -- K t2~A ／ 2 D - 2 t ) 一 ——。 l 一 — J 。 ( 4) 2河流应 急监 测布点原则 突发性环 境 污 染 事 故 应 急 监 测 是 指 在 发 生 环 境 污染 事 故 的紧 急情 况 下 ，为 发 现 和 查 明环 境 污染 情 况( 污 染 物 种 类 、污 染 范 围 和 污 染 程 度 等 )而 进行 的 环 境 监 测 ，包 括 现 场 定 点 监 测 和动态 监 测 。根 据 了解 的 事 故 信 息 及 事 故 发 生 地点 河 流周 围 的 地 形 地 貌 ，选 用 适 用 的 污染 物 泄漏模型，结合 当地 的气 象水文 等参数 ，根 据 MA TL AB预测 事故 范 围 ，并 在 电子 地 图上 显 示 事故 可 能影 响 的 范 围 ；在 影 响 范 围 内标 识 出环 境敏 感点 。 环境敏感 点是指 事故发生后 ，需要特别保 护的 区 域 ，如 人 群 稠 密 区 、学 校 、 医 院 、水 源 地等。河流监测一般应 布设对 照断 面、控制断 面 和削减 断 面 。设 置 监 测 断 面 后 ，应 根 据 水 面 的宽度 确定 断 面 上 的采 样 垂 线 ，再 根 据 采 样 垂 线处 水深确 定采 样点 的数 目和位 置 。 3 MA T L A B在松 花江 污染事故中的应用 2 0 0 5年 1 1月 1 3日 1 3时 4 0分 左 右 ，中石 油吉林石化公 司 1 0 1双苯厂发生爆 炸。事故 造 成 1 0 0 t 余笨及硝基苯泄露 。自发生爆炸事故 ， 到 2 3日中央有关部委各部 门开始组建专家组讨 论 并相 继赶 到哈 尔滨 处 理 水污 染 事 件 ，共 l O d ， 就是在 这 几 天 里 ，随 着 时 间 的 推 移 和 上 游 来 水 的稀释混合作 用，污染带 由原来 的数 千米延长 到 8 0 k m 余 。 该 问题 分 别利 用 式 ( 1 )和 式 ( 3 )求 解 ， 其中水体 的