版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、第四章第四章 地表水环境影响评价地表水环境影响评价第一节第一节 地表水体的污染和自净地表水体的污染和自净一、地表水资源一、地表水资源二、水体污染二、水体污染 点污染源点污染源 点污染源排放的废水量和污染物可以从管点污染源排放的废水量和污染物可以从管道或沟渠中直接量测流量和采样分析组分浓度道或沟渠中直接量测流量和采样分析组分浓度确定,在经费和其他条件有限制时,常采用排确定,在经费和其他条件有限制时,常采用排污指标污指标(例如排放系数例如排放系数)推算的方法推算的方法。居住区生活污水量计算式,式中:居住区生活污水量计算式,式中: QS居住区生活污水量,居住区生活污水量,L/s; q每人每日的排水定
2、额,每人每日的排水定额, L/(人人d); N设计人口数,人;设计人口数,人; Ks总变化系数总变化系数(1.51.7)。3600ism M KQt工业废水量计算式,式中:工业废水量计算式,式中: m单位产品废水量,单位产品废水量,L/t; M该产品的日产量,该产品的日产量,t; Ki总变化系数,根据工艺或总变化系数,根据工艺或 经验决定;经验决定; t 工厂每日工作时数,工厂每日工作时数,h。86400ssqNKQ 2. 2. 非点污染源非点污染源非点污染源:非点污染源:非点污染源又称面源,是指分散非点污染源又称面源,是指分散或均匀地通过岸线进入水体的废水和自然降水或均匀地通过岸线进入水体的
3、废水和自然降水通过沟渠进入水体的废水通过沟渠进入水体的废水。 主要包括城镇排水、农田排水和农村生活主要包括城镇排水、农田排水和农村生活废水、矿山废水、分散的小型禽畜饲养场废水,废水、矿山废水、分散的小型禽畜饲养场废水,以及大气污染物通过重力沉降和降水过程进入以及大气污染物通过重力沉降和降水过程进入水体等所造成的污染废水水体等所造成的污染废水。 非点源污染情况复杂,其污染影响较难定非点源污染情况复杂,其污染影响较难定量,但又不能忽视,特别是对点源已进行有效量,但又不能忽视,特别是对点源已进行有效控制后,非点源污染会日益突出。控制后,非点源污染会日益突出。(1)城市非点污染源负荷估计:城市非点污染
4、源负荷估计:城市非点污染源负荷来源:城市非点污染源负荷来源:城市雨水下水道及城市雨水下水道及合流制下水道的溢流。污染物自城市街道经排合流制下水道的溢流。污染物自城市街道经排水系统进入受纳水体。水系统进入受纳水体。城市非点源污染物被暴雨冲刷到接受水体的负城市非点源污染物被暴雨冲刷到接受水体的负荷的计算:荷的计算:基本程序:基本程序:首先估计暴雨事件中暴雨径流的大首先估计暴雨事件中暴雨径流的大小(径流深度和径流面积的乘积),从而确定小(径流深度和径流面积的乘积),从而确定暴雨的冲刷率,进而估计径流冲刷到受纳水体暴雨的冲刷率,进而估计径流冲刷到受纳水体的沉积物负荷,然后根据沉积物中污染物浓度的沉积物
5、负荷,然后根据沉积物中污染物浓度计算污染物负荷,或者根据固体废物与污染物计算污染物负荷,或者根据固体废物与污染物的统计相关关系计算污染物负荷。的统计相关关系计算污染物负荷。暴雨径流深度的估计:暴雨径流深度的估计: RCRPDs 式中:式中: R 总暴雨径流深度,总暴雨径流深度,cm; CR 总径流系数;总径流系数; P 降雨量,降雨量,cm; Ds 洼地存水,洼地存水,Cm。 总径流系数的估算方法:总径流系数的估算方法:粗略估算式:粗略估算式:式中:式中:I不透水区百分数;不透水区百分数; 按照不同坡度计算的不透水区按照不同坡度计算的不透水区(指指屋面、沥青和水泥路面或广场、庭院等屋面、沥青和
6、水泥路面或广场、庭院等)的径流的径流系数系数 。0.15 1100100RIIC准确计算式:准确计算式:式中:式中:Fi各种类型地区所占的面积;各种类型地区所占的面积; i对应的径流系数。对应的径流系数。iiRiFCF洼地存水洼地存水Ds的粗略估计:的粗略估计:0.630.48100sID 径流中冲刷到接受水体的颗粒物负荷:在总径流中冲刷到接受水体的颗粒物负荷:在总暴雨径流估算出来后,可估算暴雨冲刷率。一暴雨径流估算出来后,可估算暴雨冲刷率。一般认为般认为1 h内总径流为内总径流为1.27 cm时,可冲走时,可冲走90的的街道表面颗粒物(沉积物)。街道表面颗粒物(沉积物)。暴雨径流中冲刷的固体
7、负荷:暴雨径流中冲刷的固体负荷:式中:式中:Ysw暴雨冲刷到受纳水体的颗粒物负荷;暴雨冲刷到受纳水体的颗粒物负荷; te 等效的累积天数,等效的累积天数,d;Ysu街道表面颗粒物日负荷量,街道表面颗粒物日负荷量,kgd。swesuYt YPC 式中:式中: tr从最后一次暴雨事件算起的天数,从最后一次暴雨事件算起的天数,d; ts从最后一次清扫街道算起的天数,从最后一次清扫街道算起的天数,d;s街道清扫频率。街道清扫频率。 1erssstttt式中:式中:Lsu颗粒物日负荷率,颗粒物日负荷率,kg(kmd);Lst街道边沟长,约等于街道边沟长,约等于2倍的街道长,倍的街道长,km。susust
8、YLL 街道表面颗粒物日负荷取决于多种因素,街道表面颗粒物日负荷取决于多种因素,如交通强度、区域地表覆盖物的形式、径流量如交通强度、区域地表覆盖物的形式、径流量和降雨强度、灰尘沉降量、前期干旱时间、城和降雨强度、灰尘沉降量、前期干旱时间、城市街道清扫频率和清扫质量等。市街道清扫频率和清扫质量等。4-1径流中冲刷到受纳水体的有机污染负荷:径流中冲刷到受纳水体的有机污染负荷:用颗粒固体负荷乘上浓度因子计算有机物负荷:用颗粒固体负荷乘上浓度因子计算有机物负荷:式中:式中:You有机污染物的日负荷量,有机污染物的日负荷量,kgd; 单位转换因子,单位转换因子,10-6; Ysu总颗粒物固体日负荷量,总
9、颗粒物固体日负荷量,kgd; Cou有机污染物在颗粒物中的浓度,有机污染物在颗粒物中的浓度,gg。ousuouYYC 城市降雨径流问题是个十分复杂的问题,与城市降雨径流问题是个十分复杂的问题,与水分循环的每一个环节都有关系,并与多种因素水分循环的每一个环节都有关系,并与多种因素相关,如降水过程、大气污染、土地使用、人类相关,如降水过程、大气污染、土地使用、人类污染特征、自然特点等。由于变化性大、随机性污染特征、自然特点等。由于变化性大、随机性强、偶然因素多,尚未掌握其规律性。强、偶然因素多,尚未掌握其规律性。(2)农田径流污染负荷估算:农田径流污染负荷估算:第一种方法:避开污染物在农田表面实际
10、迁移第一种方法:避开污染物在农田表面实际迁移过程的变化,仅通过采集和分析各个集水区的过程的变化,仅通过采集和分析各个集水区的径流水样计算进入某一水环境中某种污染物总径流水样计算进入某一水环境中某种污染物总量,其公式如下:量,其公式如下:式中:式中: M某种污染物输出总量,某种污染物输出总量,kg;i第第i小时的该种污染物浓度,小时的该种污染物浓度,kgm3; Qi第第i小时的径流量,小时的径流量,m3; n观测的总时数,观测的总时数,h; j第第j个农田集水区;个农田集水区; m集水区总数。集水区总数。11mniijiMQ第二节污染物质在河流中的混合与扩散第二节污染物质在河流中的混合与扩散一、
11、污染物质在河流中的混合一、污染物质在河流中的混合 废水排入水体后,最先发生的过程是废水排入水体后,最先发生的过程是混合稀混合稀释释。对大多数保守污染物混合稀释是它们迁移的。对大多数保守污染物混合稀释是它们迁移的主要方式之一。对易降解的污染物混合稀释也是主要方式之一。对易降解的污染物混合稀释也是它们迁移的重要方式之一。水体的混合稀释、扩它们迁移的重要方式之一。水体的混合稀释、扩散能力,与其水体的水文特征密切相关。散能力,与其水体的水文特征密切相关。1河流的混合稀释模型河流的混合稀释模型 当废水进入河流后,便不断地与河水发生当废水进入河流后,便不断地与河水发生混合交换作用,使保守污染物浓度沿流程逐
12、渐混合交换作用,使保守污染物浓度沿流程逐渐降低,这一过程称为混合稀释过程。降低,这一过程称为混合稀释过程。 污水排入河流的入河口称为污水排入河流的入河口称为污水注入点污水注入点。污水。污水注入点以下的河段,污染物在断面上的浓度分注入点以下的河段,污染物在断面上的浓度分布是不均匀的,靠入河口一侧的岸边浓度高,布是不均匀的,靠入河口一侧的岸边浓度高,远离排放口对岸的浓度低。随着河水的流逝,远离排放口对岸的浓度低。随着河水的流逝,污染物在整个断面上的分布逐渐均匀。污染物在整个断面上的分布逐渐均匀。污染物浓度在整个断面上变为均匀一致的断面,污染物浓度在整个断面上变为均匀一致的断面,称为称为水质完全混合
13、断面水质完全混合断面。最早出现水质完全混合断面的位置称为最早出现水质完全混合断面的位置称为完全混完全混合点合点。污水注入点的上游称为污水注入点的上游称为初始段初始段,或,或背景河段背景河段;污水注入点到完全混合点之间的河段称为污水注入点到完全混合点之间的河段称为非均非均匀混合段匀混合段;完全混合点的下游河段称为完全混合点的下游河段称为均匀混合段均匀混合段。 设河水流量为设河水流量为Q,水质完全混合断面以前,水质完全混合断面以前,任一非均匀混合断面上参与和废水混合的河水任一非均匀混合断面上参与和废水混合的河水流量为流量为Qi。把参与和废水混合的河水流量。把参与和废水混合的河水流量Qi与与该断面河
14、水流量该断面河水流量Q的比值定义为的比值定义为混合系数混合系数,以,以表示:表示: 把参与和废水混合的河水流量把参与和废水混合的河水流量Qi与废水流与废水流量量q的比值定义为的比值定义为稀释比稀释比,以,以n表示:表示:iQQiQQnqq非均匀混合断面上的污染物平均浓度计算公式:非均匀混合断面上的污染物平均浓度计算公式:式中:式中: Q河流的流量,河流的流量,m3s;1排污口上游河流中污染物浓度,排污口上游河流中污染物浓度,mgL; q排人河流的废水流量,排人河流的废水流量,m3s;2废水中的污染物浓度,废水中的污染物浓度,mgL。在水质完全混合断面以下的任一断面的在水质完全混合断面以下的任一
15、断面的、n和和i均为常数。均为常数。1212iiiQqQqQqQq12iQqQq 当废水在岸边排入河流时,废水靠岸边向当废水在岸边排入河流时,废水靠岸边向下游流去,经过相当长的距离才能达到完全混下游流去,经过相当长的距离才能达到完全混合。在非均匀混合段的废水排入一侧的岸边形合。在非均匀混合段的废水排入一侧的岸边形成一个污染带。当完全混合距离成一个污染带。当完全混合距离Ln无实测数据无实测数据时,可参考下表确定。表中列举出了许多河流时,可参考下表确定。表中列举出了许多河流在岸边集中排入废水时,污水与河水达到完全在岸边集中排入废水时,污水与河水达到完全混合所需的时间。混合所需的时间。从下表中查取所
16、需时间与河从下表中查取所需时间与河水实际流速的乘积为完全混合距离水实际流速的乘积为完全混合距离。二、污染物质在河流中的扩散二、污染物质在河流中的扩散 污染物质在河流中的迁移总起来可分为两污染物质在河流中的迁移总起来可分为两类,即类,即推流推流和和扩散扩散。推流也称平流、随流输移。推流也称平流、随流输移。推流是指污染物质随水质点的流动一起移到新推流是指污染物质随水质点的流动一起移到新的位置。扩散可分为分子扩散、湍流扩散、剪的位置。扩散可分为分子扩散、湍流扩散、剪切流离散(弥散)和对流扩散。切流离散(弥散)和对流扩散。1分子扩散分子扩散 分子扩散是指物质分子的随机运动(即布分子扩散是指物质分子的随
17、机运动(即布朗运动)而引起的物质迁移或分散现象朗运动)而引起的物质迁移或分散现象。当水。当水体中污染物质浓度分布不均匀时,污染物质将体中污染物质浓度分布不均匀时,污染物质将会从浓度高的地方向浓度低的会从浓度高的地方向浓度低的地方移动。分子扩散过程服从地方移动。分子扩散过程服从费克第一定律。费克第一定律。xMJEx 即以扩散方式通过单位截面积的质量流量与扩即以扩散方式通过单位截面积的质量流量与扩散物质的浓度梯度成正比。散物质的浓度梯度成正比。 分子扩散系数一般很小。分子扩散引起的分子扩散系数一般很小。分子扩散引起的物质迁移与其它因素引起物质迁移相比,分子物质迁移与其它因素引起物质迁移相比,分子扩
18、散在水环境影响评价中往往被忽略。扩散在水环境影响评价中往往被忽略。2湍流扩散湍流扩散 当河流做湍流运动时,随机的湍流作用引当河流做湍流运动时,随机的湍流作用引起污染物的扩散,称为湍流扩散起污染物的扩散,称为湍流扩散。 湍流扩散所引起的污染物质量通量与浓度湍流扩散所引起的污染物质量通量与浓度梯度成正比。湍流扩散系数比梯度成正比。湍流扩散系数比分子扩散系数大分子扩散系数大78个数量级。个数量级。因此,在河流中污染物的迁移因此,在河流中污染物的迁移是以湍流为主的。是以湍流为主的。xxIEx 3剪切流离散剪切流离散 当垂直于流动方向的横断面上流速分布不当垂直于流动方向的横断面上流速分布不均匀或者说有流
19、速梯度存在的流动称为剪切流。均匀或者说有流速梯度存在的流动称为剪切流。剪切流离散又称弥散,它是由于横断面上各点剪切流离散又称弥散,它是由于横断面上各点的实际流速不等而引起的。的实际流速不等而引起的。 剪切流离散同样可以类比分子扩散,其引剪切流离散同样可以类比分子扩散,其引起的质量通量可用下式表示:起的质量通量可用下式表示:式中式中 :Dx剪切流离散系剪切流离散系数,或称弥漫系数,数,或称弥漫系数,m2s。xxJDx 4对流扩散对流扩散 对流扩散指由于温度差或密度分层不稳定对流扩散指由于温度差或密度分层不稳定性面引起的铅直方向对流运动所伴随的污染物性面引起的铅直方向对流运动所伴随的污染物迁移迁移
20、。 在自然界的水体中,各种形式扩散常常交在自然界的水体中,各种形式扩散常常交织在一起发生,除上述污染物几种主要迁移方织在一起发生,除上述污染物几种主要迁移方式以外,还存在着冲刷、淤积和悬浮等多种形式以外,还存在着冲刷、淤积和悬浮等多种形式式。除分子扩散外,所有各种迁移方式都和水。除分子扩散外,所有各种迁移方式都和水体流动特性有密切联系,因此,要研究物质的体流动特性有密切联系,因此,要研究物质的扩散输移规律应和研究水体的流动特性紧紧联扩散输移规律应和研究水体的流动特性紧紧联系在一起。系在一起。5移流扩散方程移流扩散方程 从流动的水体中,取一微分六面体。按照从流动的水体中,取一微分六面体。按照物质
21、守恒原理,从微分六面体流进与流出的污物质守恒原理,从微分六面体流进与流出的污染物质量之差应当等于同时段内微分六面体内染物质量之差应当等于同时段内微分六面体内质量的增量,从而导出三维的移流扩散方程为:质量的增量,从而导出三维的移流扩散方程为:对于二维问题,移流扩散方程为:对于二维问题,移流扩散方程为:xyuvwtxyzEEsxxyyzzxyuvtxyEEsxxyy对于一维问题,移流扩散方程为:对于一维问题,移流扩散方程为: 基本模型的求解因环境问题的复杂,往往求解基本模型的求解因环境问题的复杂,往往求解起来很困难,通常是利用有限差分法和有限元起来很困难,通常是利用有限差分法和有限元法求其数值解。
22、法求其数值解。xuEstxxx第三节第三节 河流和河口水质模型河流和河口水质模型 河流是沿地表的线形低凹部分集中的经常河流是沿地表的线形低凹部分集中的经常性或周期性水流。较大的叫河(或江),较小性或周期性水流。较大的叫河(或江),较小的叫溪。河口是河流注入海洋、湖泊或其他河的叫溪。河口是河流注入海洋、湖泊或其他河流的河段,可以分为入海河口、入湖河口及支流的河段,可以分为入海河口、入湖河口及支流河口。流河口。 应用水质模型预测河流水质时,应用水质模型预测河流水质时,常假设该常假设该河段内无支流,在预测时期内河段的水力条件河段内无支流,在预测时期内河段的水力条件是稳态的和只在河流的起点有恒定浓度和
23、流量是稳态的和只在河流的起点有恒定浓度和流量的废水(或污染物)排入的废水(或污染物)排入。 如果在如果在河段内有支流汇入,而且沿河有多河段内有支流汇入,而且沿河有多个污染源,这时应将河流划分为多个河段采用个污染源,这时应将河流划分为多个河段采用多河段模型多河段模型。2河流水质模型河流水质模型 河流水质模型是描述水体中污染物随时间河流水质模型是描述水体中污染物随时间和空间迁移转化规律的数学方程。和空间迁移转化规律的数学方程。1、水质模型的分类:、水质模型的分类:n按时间特性分:按时间特性分:分为动态模型和静态模型分为动态模型和静态模型。 描写水体中水质组分的浓度随时间变化的描写水体中水质组分的浓
24、度随时间变化的水质模型称为水质模型称为动态模型动态模型。 描述水体中水质组分的浓度不随时间变化描述水体中水质组分的浓度不随时间变化的水质模型称为的水质模型称为静态模型静态模型。 n按水质模型的空间维数分:按水质模型的空间维数分:分为零维、一维、二维、分为零维、一维、二维、三维水质模型三维水质模型。 当把所考察的水体看成是一个完全混合反应当把所考察的水体看成是一个完全混合反应器时,即水体中水质组分的浓度是均匀分布的,描述器时,即水体中水质组分的浓度是均匀分布的,描述这种情况的水质模型称为零维的水质模型。这种情况的水质模型称为零维的水质模型。 描述水质组分的迁移变化在一个方向上是重要描述水质组分的
25、迁移变化在一个方向上是重要的,另外两个方向上是均匀分布的,这种水质模型称的,另外两个方向上是均匀分布的,这种水质模型称为一维水质模型。为一维水质模型。 描述水质组分的迁移变化在两个方向上是重描述水质组分的迁移变化在两个方向上是重要的,在另外的一个方向上是均匀分布的,这种水质要的,在另外的一个方向上是均匀分布的,这种水质模型称为两维水质模型。模型称为两维水质模型。 描述水质组分迁移变化在三个方向进行的水描述水质组分迁移变化在三个方向进行的水质模型称为三维水质模型。质模型称为三维水质模型。n按描述水质组分的多少分:按描述水质组分的多少分:分为单一组分和分为单一组分和多组分的水质模型多组分的水质模型
26、。 水体中某一组分的迁移转化与其它组分水体中某一组分的迁移转化与其它组分没没有关系,描述这种组分迁移转化的水质模型称有关系,描述这种组分迁移转化的水质模型称为单一组分的水质模型。为单一组分的水质模型。 水体中一组分的迁移转化与另一组分(或水体中一组分的迁移转化与另一组分(或几个组分)的迁移转化是相互联系、相互影响几个组分)的迁移转化是相互联系、相互影响的,描述这种情况的水质模型称为多组分的水的,描述这种情况的水质模型称为多组分的水质模型。质模型。n按水体的类型可分为:按水体的类型可分为:河流水质模型、河口水河流水质模型、河口水质模型(受潮汐影响)、湖泊水质模型、水库质模型(受潮汐影响)、湖泊水
27、质模型、水库水质模型和海湾水质模型等水质模型和海湾水质模型等。河流、河口水质。河流、河口水质模型比较成熟,湖、海湾水质模型比较复杂,模型比较成熟,湖、海湾水质模型比较复杂,可靠性小。可靠性小。n按水质组分可分为:按水质组分可分为:耗氧有机物模型(耗氧有机物模型(BODBODDODO模型)模型) ,无机盐、悬浮物、放射性物质等,无机盐、悬浮物、放射性物质等单一组分的水质模型,难降解有机物水质模型,单一组分的水质模型,难降解有机物水质模型,重金属迁移转化水质模型重金属迁移转化水质模型。 水质模型的选择:水质模型的选择:选择水质模型必须对所研究的水选择水质模型必须对所研究的水质组分的迁移转化规律有清
28、楚地了解。因为水质组分质组分的迁移转化规律有清楚地了解。因为水质组分的迁移的迁移(扩散和平流扩散和平流)取决于水体的水文特性和水动力学取决于水体的水文特性和水动力学特性。在流动的河流中,平流迁移往往占主导地位,特性。在流动的河流中,平流迁移往往占主导地位,对某些组分可以忽略扩散项;在受潮汐影响的河口中,对某些组分可以忽略扩散项;在受潮汐影响的河口中,扩散是主导的迁移现象,扩散项必须考虑而不能忽略。扩散是主导的迁移现象,扩散项必须考虑而不能忽略。对这两者选择的模型就不应一样。对河床规整,断面对这两者选择的模型就不应一样。对河床规整,断面不变,污染物排人量不变的水体,可选用静态模型。不变,污染物排
29、人量不变的水体,可选用静态模型。为了减少模型的复杂性和减少所需的资料,对河流系为了减少模型的复杂性和减少所需的资料,对河流系统的水质模型往往选用静态的。但这种选择不能充分统的水质模型往往选用静态的。但这种选择不能充分评价时便输入对河流系统的影响。评价时便输入对河流系统的影响。 选择的水质模型必须反映所研究的水质组分,应选择的水质模型必须反映所研究的水质组分,应用条件和现实条件接近。用条件和现实条件接近。2、污染物在均匀流场中的扩散水质模型、污染物在均匀流场中的扩散水质模型 进入环境的污染物可以分为两大类:守恒进入环境的污染物可以分为两大类:守恒污染物(惰性污染物)和非守恒污染物。污染物(惰性污
30、染物)和非守恒污染物。守恒污染物守恒污染物:污染物进入环境以后,随着介质:污染物进入环境以后,随着介质的运动不断地变换所处的空间位置,还由于分的运动不断地变换所处的空间位置,还由于分散作用不断向周围扩散而降低其初始浓度,但散作用不断向周围扩散而降低其初始浓度,但它不会因此而改变总量,不发生衰减。这种污它不会因此而改变总量,不发生衰减。这种污染物称为守恒污染物。如重金属、很多高分子染物称为守恒污染物。如重金属、很多高分子有机化合物等。有机化合物等。非守恒污染物非守恒污染物:污染物进入环境以后,除了随:污染物进入环境以后,除了随着环境介质流动而改变位置,并不断扩散而降着环境介质流动而改变位置,并不
31、断扩散而降低浓度外,还因自身的衰减而加速浓度的下降。低浓度外,还因自身的衰减而加速浓度的下降。这种污染物称为非守恒污染物。这种污染物称为非守恒污染物。非守恒物质的衰减有两种方式:非守恒物质的衰减有两种方式:一种是由其自一种是由其自身的运动变化规律决定的身的运动变化规律决定的;如放射性物质的蜕;如放射性物质的蜕变;变;另一种是在环境因素的作用下,由于化学另一种是在环境因素的作用下,由于化学的或生物化学的反应而不断衰减的的或生物化学的反应而不断衰减的,如可生化,如可生化降解的有机物在水体中微生物作用下的氧化分降解的有机物在水体中微生物作用下的氧化分解过程。解过程。守恒污染物在均匀流场中的扩散方程守
32、恒污染物在均匀流场中的扩散方程 对于守恒污染物在运动过程中不发生衰减,对于守恒污染物在运动过程中不发生衰减,在上述移流扩散方程中应有在上述移流扩散方程中应有S0。在均匀流场。在均匀流场中,流速应为常数,扩散参数也应为常数。因中,流速应为常数,扩散参数也应为常数。因此,移流扩散方程式有下列形式:此,移流扩散方程式有下列形式:二维空间扩散方程式为:二维空间扩散方程式为:一维空间扩散方程式为:一维空间扩散方程式为:扩散方程的解扩散方程的解对于瞬时点源,守恒污染物在均匀无限大流场对于瞬时点源,守恒污染物在均匀无限大流场中,污染物浓度呈高斯分布。若设坐标原点在中,污染物浓度呈高斯分布。若设坐标原点在污染
33、物排放点,则有:污染物排放点,则有:二维扩散方程的解:二维扩散方程的解:2222dEdtxy22dEdtx22exp44xutyMEtEt一维扩散方程的解为:一维扩散方程的解为:对于守恒污染物,实际应用中,在不需要考虑对于守恒污染物,实际应用中,在不需要考虑其横向均匀混合时间的情况下,通常假设其可其横向均匀混合时间的情况下,通常假设其可以瞬间混合完毕,而采用完全混合公式来计算以瞬间混合完毕,而采用完全混合公式来计算河流断面的污染物浓度。河流断面的污染物浓度。2()exp44MxutEtEt 对非守恒污染物,在河流的流量和其他水对非守恒污染物,在河流的流量和其他水文条件不变的稳态条件下,可以采用
34、一维模型文条件不变的稳态条件下,可以采用一维模型进行污染物浓度预测。进行污染物浓度预测。220 xxEuKxx 对于非持久性或可降解污染物,若给定对于非持久性或可降解污染物,若给定x0,0,上式解为:,上式解为: 对于一般条件下的河流,推流形成的污染物迁移作对于一般条件下的河流,推流形成的污染物迁移作用要比弥散作用大得多,在稳态条件下,弥散作用用要比弥散作用大得多,在稳态条件下,弥散作用可以忽略,则有:可以忽略,则有:式中:式中: ux河流的平均流速,河流的平均流速,md或或ms;Ex废水与河水的纵向混合系数,废水与河水的纵向混合系数,m2d或或m2s; K污染物的衰减系数,污染物的衰减系数,
35、1d或或1s; x河水河水(从排放口从排放口)向下游流经的距离,向下游流经的距离,m。04exp112xxxxu xKEEu0expxxKu例例1 一个改扩建工程拟向河流排放废水,废水一个改扩建工程拟向河流排放废水,废水量量q0.15m3s,苯酚浓度为,苯酚浓度为30gL,河流流,河流流量量Q5.5m3s,流速,流速u0.3ms,苯酚背景浓,苯酚背景浓度为度为 0.5 g L,苯酚的降解(衰减)系数,苯酚的降解(衰减)系数K0.2d-1,纵向弥散系数,纵向弥散系数Ex10m2s。求排放点。求排放点下游下游10km处的苯酚浓度。处的苯酚浓度。解解 计算起始点处完全混合后的初始浓度:计算起始点处完
36、全混合后的初始浓度:(1)考虑纵向弥散条件下的下游)考虑纵向弥散条件下的下游10km处的浓处的浓度:度:00.15305.50.51.28/5.50.15g L24 0.2/86400 100.3 100001.28 exp111.19/2 100.3g L(2)忽略纵向弥散时的下游)忽略纵向弥散时的下游10km处的浓度:处的浓度: 0.2 100001.281.19/0.3 86400g L 由此看出,由此看出,在稳态条件下,忽略纵向弥散在稳态条件下,忽略纵向弥散系数与考虑纵向弥散系数的差异可以忽略。系数与考虑纵向弥散系数的差异可以忽略。 对水面宽阔的河流受纳污对水面宽阔的河流受纳污(废废)
37、水后的混合水后的混合过程和污染物的衰减可用二维模型预测;对于过程和污染物的衰减可用二维模型预测;对于水面又宽又深和流态复杂的河流水质预测宜采水面又宽又深和流态复杂的河流水质预测宜采用三维模型。用三维模型。 3污染物与河水完全混合所需距离污染物与河水完全混合所需距离 污染物从排污口排出后要与河水完全混合需污染物从排污口排出后要与河水完全混合需一定的纵向距离,这段距离称为混合过程段。一定的纵向距离,这段距离称为混合过程段。 当某一断面上任意点的浓度与断面平均浓当某一断面上任意点的浓度与断面平均浓度之比介于度之比介于0.95 至至1.05 之间时,称该断面已达之间时,称该断面已达到横向混合,由排放点
38、至完成横向断面混合的到横向混合,由排放点至完成横向断面混合的距离称为完成横向混合所需的距离。距离称为完成横向混合所需的距离。 当采用河中心排放时所需的完成横向混合当采用河中心排放时所需的完成横向混合的距离为:的距离为:在岸边上排时:在岸边上排时:20.1xyu BxE20.4xyu BxE二、二、BODDO耦合模型耦合模型 河水中溶解氧浓度河水中溶解氧浓度 (DO)是决定水质洁净程是决定水质洁净程度的重要参数之一,而排入河流的度的重要参数之一,而排入河流的 BOD在衰减在衰减过程中将不断消耗过程中将不断消耗DO,与此同时空气中的氧气,与此同时空气中的氧气又不断溶解到河水中。又不断溶解到河水中。
39、 描述一维河流中描述一维河流中BOD 和和DO消长变化规律消长变化规律的模型的模型(SP模型模型)。建立。建立SP模型的基本假设模型的基本假设如下:如下:河流中的河流中的BOD的衰减和溶解氧的复氧都是一的衰减和溶解氧的复氧都是一级反应;级反应;反应速度是定常的;反应速度是定常的;河流中的耗氧是由河流中的耗氧是由BOD衰减引起的,而河流衰减引起的,而河流中的溶解氧来源则是大气复氧。中的溶解氧来源则是大气复氧。S P方程:方程:临界氧亏发生的时间:临界氧亏发生的时间: 该方程是应用最广的河流水质中该方程是应用最广的河流水质中BODDO预测模型。预测模型。01220121ssDODODBODK tK
40、 tK tDODKeeeKK0021221111ln1DcBODKKKtKKKKSP模型的修正模型:模型的修正模型: SP模型的假设是不模型的假设是不完全符合实际的。为了计算河流水质的某些特完全符合实际的。为了计算河流水质的某些特殊问题,人们在殊问题,人们在 SP 模型的基础上附加了一些模型的基础上附加了一些新的假设,推导出了一些新的模型。新的假设,推导出了一些新的模型。(1)托马斯)托马斯(Thomas)模型模型 对一维静态河流,在对一维静态河流,在SP模型的基础上,模型的基础上,为了考虑沉淀、絮凝、冲刷和再悬浮过程对为了考虑沉淀、絮凝、冲刷和再悬浮过程对BOD去除的影响,引入了去除的影响,
41、引入了BOD沉浮系数沉浮系数k3,BOD变化速度为变化速度为k3L。由以下的基本方程组。由以下的基本方程组(忽忽略扩散项略扩散项):1312BODBODDBODDdkkdtdkkdt 解得:解得:(2)多宾斯)多宾斯坎普坎普(DobbinsCamp)模型模型 对一维静态河流,在托马斯模型的基础上,对一维静态河流,在托马斯模型的基础上,多宾斯多宾斯坎普提出了两条新的假设:坎普提出了两条新的假设: 考虑地面径流和底泥释放考虑地面径流和底泥释放BOD所引起的所引起的BOD变化速率,该速率以变化速率,该速率以R表示。表示。 考虑藻类光合作用和呼吸作用以及地面径考虑藻类光合作用和呼吸作用以及地面径流所引
42、起的溶解氧变化速率,该速率以流所引起的溶解氧变化速率,该速率以P表示。表示。多宾斯多宾斯坎普采用以下基本方程组:坎普采用以下基本方程组:1300132201132ekktBODBODBODkktk tk tDDkeeekkk三、污染物在河口中的混合和衰减模型三、污染物在河口中的混合和衰减模型 入海河口受海洋潮汐和上游河流来水双重入海河口受海洋潮汐和上游河流来水双重作用。海潮上溯与上游下泄的水流相汇形成强作用。海潮上溯与上游下泄的水流相汇形成强烈的混合作用。烈的混合作用。 一般污染比较严重的河口都是工业集中的一般污染比较严重的河口都是工业集中的城市或水陆交通枢纽。在无组织排放的条件下,城市或水陆
43、交通枢纽。在无组织排放的条件下,河口将受纳许多排放口废水。在通航的河口,河口将受纳许多排放口废水。在通航的河口,其宽度一般都较大,也比较深,污染物要完成其宽度一般都较大,也比较深,污染物要完成横向混合仍需要经过较长的距离。横向混合仍需要经过较长的距离。 当只需了解污染物在当只需了解污染物在个潮汐周期内的个潮汐周期内的平均浓度时,可以采用本节中介绍的河流相应平均浓度时,可以采用本节中介绍的河流相应情况的模型,其混合系数情况的模型,其混合系数Ey可以采用式可以采用式(467)的泰勒公式。的泰勒公式。 如果要求污染物与河口水混合过程中浓度如果要求污染物与河口水混合过程中浓度随时间变化情况,则应采用二
44、维动态混合数值随时间变化情况,则应采用二维动态混合数值模型预测:首先通过实测得到断面上各测点流模型预测:首先通过实测得到断面上各测点流速与断面平均流速的相关关系,同时用一维非速与断面平均流速的相关关系,同时用一维非恒定流方程数值模型计算出沿程各断面平均流恒定流方程数值模型计算出沿程各断面平均流速,这样就可得到河口的流场分布。速,这样就可得到河口的流场分布。二维动态混合物数值模型的微分方程见式:二维动态混合物数值模型的微分方程见式:2222xyxyuuEEKtxyxy四、河口和河网水质模型四、河口和河网水质模型 河口是入海河流受潮汐作用影响明显的河河口是入海河流受潮汐作用影响明显的河段。段。潮汐
45、对河口水质的双重影响:潮汐对河口水质的双重影响:上游下泄的水流相汇,形成强烈的混合作用,上游下泄的水流相汇,形成强烈的混合作用,使污染物的分布趋于均匀;使污染物的分布趋于均匀;由于潮流的顶托作用,延长了污染物在河口由于潮流的顶托作用,延长了污染物在河口的停留时间,有机物的降解会进一步消耗水中的停留时间,有机物的降解会进一步消耗水中的溶解氧,使水质下降。的溶解氧,使水质下降。此外,潮汐也使河口的含盐量增加。此外,潮汐也使河口的含盐量增加。 河口模型比河流模型复杂,求解也比较困河口模型比河流模型复杂,求解也比较困难。对河口水质有重大影响的评价项目,需要难。对河口水质有重大影响的评价项目,需要预测污
46、染物浓度随时间的变化。这时应采用水预测污染物浓度随时间的变化。这时应采用水力学中的非恒定流的数值模型,以差分法计算力学中的非恒定流的数值模型,以差分法计算流场,再采用动态水质模型,预测河口任意时流场,再采用动态水质模型,预测河口任意时刻的水质。当排放口的废水能在断面上与河水刻的水质。当排放口的废水能在断面上与河水迅速充分混合,则也可用一维非恒定流数值模迅速充分混合,则也可用一维非恒定流数值模型计算流场,再用一维动态水质模型预测任意型计算流场,再用一维动态水质模型预测任意时刻的水质。对河口水质有重大影响,但只需时刻的水质。对河口水质有重大影响,但只需预测污染物在一个潮汐周期内的平均浓度,这预测污
47、染物在一个潮汐周期内的平均浓度,这时可以用一维潮周平均模型预测。时可以用一维潮周平均模型预测。 一维(潮周平均)河口水质模型如下:一维(潮周平均)河口水质模型如下:式中:式中:r污染物的衰减速率,污染物的衰减速率,g/(m3.d); s系统外输入污染物的速率,系统外输入污染物的速率,g/(m3.d); ux不考虑潮汐作用,由上游来水不考虑潮汐作用,由上游来水(净泄净泄量量)产生的流速,产生的流速,m/s。假定假定s0和和rK1,对排放点上游(对排放点上游(x0)对排放点下游(对排放点下游(x 0)0 xxdddEursdxdxdx 1204exp112xxxxuK EEu1204exp112x
48、xxxuK EEu01241xxWK EQu 第三节第三节 湖泊湖泊(水库水库)水质数学模型水质数学模型 湖泊(水库)水流状态分为前进和振动两类。湖泊(水库)水流状态分为前进和振动两类。前者指湖流和混合作用,后者指波动和波漾。前者指湖流和混合作用,后者指波动和波漾。 (1)湖流:指湖水在水力坡度、密度梯度和风力等湖流:指湖水在水力坡度、密度梯度和风力等作用下产生沿一定方向的缓慢流动。湖流经常呈水作用下产生沿一定方向的缓慢流动。湖流经常呈水平环状运动(多出现在湖水较浅的场合)和垂直环平环状运动(多出现在湖水较浅的场合)和垂直环状运动状运动(湖水较深时湖水较深时)。 (2)混合:指在风力和水力坡度
49、作用下产生的湍流混合:指在风力和水力坡度作用下产生的湍流混合和由湖水密度差引混合和由湖水密度差引 起的对流混合作用。起的对流混合作用。 (3)波动:主要由风引起的,又称风浪。波动:主要由风引起的,又称风浪。 (4)波漾:是在复杂的外力作用下,湖中水位有节波漾:是在复杂的外力作用下,湖中水位有节奏的升降变化。奏的升降变化。湖泊湖泊(水库水库)的水质特征:的水质特征:水的停留时间较长(可达数月至数年),属于水的停留时间较长(可达数月至数年),属于缓流水域,其中的化学和生物学过程保持一个比缓流水域,其中的化学和生物学过程保持一个比较稳定的状态。较稳定的状态。进入湖泊和水库中的营养物质在其中容易不断进
50、入湖泊和水库中的营养物质在其中容易不断积累,致使水质发生富营养化。积累,致使水质发生富营养化。在水深较大的湖、库中,水温和水质是竖向分在水深较大的湖、库中,水温和水质是竖向分层的。层的。 湖泊水质模型分为描述湖、库营养状况的湖泊水质模型分为描述湖、库营养状况的箱箱式模型式模型、分层箱式模型分层箱式模型和描述温度与水质竖向分和描述温度与水质竖向分布的布的分层模型分层模型。一、完全混合模型一、完全混合模型 完全混合模型属箱式模型,也称沃兰伟德完全混合模型属箱式模型,也称沃兰伟德(Vollenwelder)模型。模型。 对于停留时间很长、水质基本处于稳定状态对于停留时间很长、水质基本处于稳定状态的中
51、小型湖泊和水库,可以简化为一个均匀混合的中小型湖泊和水库,可以简化为一个均匀混合的水体。沃兰伟德假定,湖泊中某种营养物的浓的水体。沃兰伟德假定,湖泊中某种营养物的浓度随时间的变化率,是输入、输出和在湖泊内沉度随时间的变化率,是输入、输出和在湖泊内沉积的该种营养物量的函数,可以用质量平衡方程积的该种营养物量的函数,可以用质量平衡方程表示:表示: 1污染物污染物(营养物营养物)混合和降解模型混合和降解模型1dVWQKVdt式中:式中: V湖、库的容积,湖、库的容积,m3; 污染物或水质参数的浓度,污染物或水质参数的浓度,mgL; 污染物或水质参数的平均排入量,污染物或水质参数的平均排入量,mgs;
52、 t时间,时间,s; Q出入湖、库流量,出入湖、库流量,m3s;K1-污染物或水质参数浓度衰减速率系数污染物或水质参数浓度衰减速率系数1s。积分上式得:积分上式得:式中:式中: W0现有污染物排入量,现有污染物排入量,mgs;W11exptWQKtQKVV0pWWq 拟建项目废水中污染物浓度,拟建项目废水中污染物浓度,mgL; q废水排放量,废水排放量,m3s。而而 式中:式中: 湖、库中污染物起始浓度,湖、库中污染物起始浓度,mgL。则:。则:对于持久性污染物对于持久性污染物K10,则:,则:当时间足够长,湖、库中污染物当时间足够长,湖、库中污染物(营养物营养物)浓度达浓度达到平衡时,到平衡
53、时, 。则平衡时浓度为:。则平衡时浓度为: p10WQK V0011 tttWQeeKVVQV0ddteWV2求湖、库中污染物达到一指定求湖、库中污染物达到一指定t所需时间所需时间t0。设设t/p=,则:,则:3无污染物输入无污染物输入(W0)时浓度随时间变化为时浓度随时间变化为 这时,可以求出污染物这时,可以求出污染物(营养物营养物)浓度达到初始浓浓度达到初始浓度之比为度之比为即即t 0 时,所需时间:时,所需时间: 1ln(1)VtQKV1(/)00Q VKttee11lnt4溶解氧模型溶解氧模型式中:式中:K2大气复氧系数,大气复氧系数,1/d或或1/s;DO0溶解氧起始浓度,溶解氧起始
54、浓度,mg/L;R湖库的生物和非生物因素耗氧总量,湖库的生物和非生物因素耗氧总量,mg/(m3.d)或或mg/(m3s); R=rABA养鱼密度,养鱼密度,kg/m3; r鱼类耗氧速率,鱼类耗氧速率,mg/(kg.d)或或mg/(kgs);DOs饱和溶解氧浓度,饱和溶解氧浓度,mg/L;B其他因素耗氧量,其他因素耗氧量,mg/(m3.d)或或mg/(m3s);02sDODODODODOdQKRdtV第四节第四节 水质模型的标定水质模型的标定 上述地表水体水质模型标定的目的是确定上述地表水体水质模型标定的目的是确定模型中各个系数,包括模型中各个系数,包括K1、K2、Ex、Ey等的值,等的值,这是
55、决定预测结果准确性的关键之一,其估值这是决定预测结果准确性的关键之一,其估值可以单个进行,也可用同时估值法;单个估值可以单个进行,也可用同时估值法;单个估值可以实测、应用经验式计算或借用类似水体的可以实测、应用经验式计算或借用类似水体的经验数据。经验数据。一、工业建设项目一、工业建设项目1建设期影响建设期影响工业建设项目在建设期工业建设项目在建设期(施工阶段施工阶段)的共同影响:的共同影响:(1)施工队伍大批进入现场,排放的生活污水和施工队伍大批进入现场,排放的生活污水和垃圾的污染。垃圾的污染。(2)施工机械运作、清洗、漏油等排放的含油和施工机械运作、清洗、漏油等排放的含油和悬浮物废水。悬浮物
56、废水。第五节第五节 开发行动对地表水影响的识别开发行动对地表水影响的识别 各种类型的人类开发行动如建设项目、区各种类型的人类开发行动如建设项目、区域和流域开发等都会对地表水环境的域和流域开发等都会对地表水环境的水量、水水量、水质、水生生物或底部沉积物质、水生生物或底部沉积物产生影响。建设项产生影响。建设项目和区域或流域开发行动在其目和区域或流域开发行动在其建设期、运行期建设期、运行期和服务期满和服务期满都会都会有不同性质和程度的影响有不同性质和程度的影响。(3)(3)基坑开挖和降低地下水位等操作排放含泥砂基坑开挖和降低地下水位等操作排放含泥砂废水。废水。(4)(4)施工场地清理和开辟施工机械通
57、行道路常大施工场地清理和开辟施工机械通行道路常大片破坏地面植被,造成裸土。在降雨片破坏地面植被,造成裸土。在降雨( (特别是暴特别是暴雨雨) )时,造成土壤侵蚀,使地表水中泥砂含量陡时,造成土壤侵蚀,使地表水中泥砂含量陡增,严重时造成河道阻塞。如果地表受过污染,增,严重时造成河道阻塞。如果地表受过污染,则污染物随雨水进入河道。则污染物随雨水进入河道。2 2运行期影响运行期影响 任何工业建设项目都有其特殊性,所以,任何工业建设项目都有其特殊性,所以,必须针对具体项目开展深入细微的工程分析才能必须针对具体项目开展深入细微的工程分析才能全面而有重点地识别出具体影响。全面而有重点地识别出具体影响。(1
58、)(1)石油炼制工业:一个炼油厂有四种主要操作:石油炼制工业:一个炼油厂有四种主要操作:分离、转化、精制和调和。分离、转化、精制和调和。废水主要来自:废水主要来自:含油废水主要来自油罐区和操作区的雨水、含油废水主要来自油罐区和操作区的雨水、油罐排水、冷却水排污、冲洗和清洗水及原油油罐排水、冷却水排污、冲洗和清洗水及原油脱盐等场所和工序;脱盐等场所和工序;苯酚、苯和有机酸等有机物以及硫化铵、金苯酚、苯和有机酸等有机物以及硫化铵、金属盐、无机盐等无机物来自汽提、原油裂解、属盐、无机盐等无机物来自汽提、原油裂解、洗涤、洗涤、油的化学处理、原油脱盐、催化裂解等工艺过油的化学处理、原油脱盐、催化裂解等工
59、艺过程;程;高温水高温水( (非污染水非污染水) )来自锅炉排污、冷却水排来自锅炉排污、冷却水排放等。放等。 炼油厂用水量和废水排放量都很大。炼油厂用水量和废水排放量都很大。(2)(2)钢铁工业:铁和钢制造钢铁工业:铁和钢制造般有六种操作:般有六种操作:焦炭制造和副产品回收;焦炭制造和副产品回收;铁矿石制备;铁矿石制备;高炉炼铁;高炉炼铁;转炉、电炉或平炉炼钢;转炉、电炉或平炉炼钢;铸、轧机操作:铸、轧机操作:精整操作。精整操作。废水主要来源如下:废水主要来源如下:a a焦炭生产和副产品回收过程的工艺用水和冷焦炭生产和副产品回收过程的工艺用水和冷却水如熄焦废水、酸洗废水和氨蒸馏废液中含高却水如
60、熄焦废水、酸洗废水和氨蒸馏废液中含高浓度酚、氰化物、硫氰酸盐和硫化物、氯化物;浓度酚、氰化物、硫氰酸盐和硫化物、氯化物;b b高炉炼铁的废水主要由排气洗涤和高炉炉渣高炉炼铁的废水主要由排气洗涤和高炉炉渣用水淬熄时排出的;用水淬熄时排出的;c c铸造和轧机操作主要排大量冷却水铸造和轧机操作主要排大量冷却水( (一般循环一般循环回用回用) ),轧机操作中产生的含铁碎屑和油滴的废,轧机操作中产生的含铁碎屑和油滴的废水;水; d d精整操作采用酸、碱浸渍去除锈和磷皮将排精整操作采用酸、碱浸渍去除锈和磷皮将排出含铁盐的酸性废水,含皂化油的碱性废水等。出含铁盐的酸性废水,含皂化油的碱性废水等。(3)(3)
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 五金件供货合同(2026版)
- 学校物业服务合同(2026版)
- 2023年6月卫生资格考试呼吸内科肺癌预防测试题及答案
- 庆阳环县高中学校选调专任教师笔试真题2025
- 赣州市面向大学生乡村医生专项计划招聘笔试真题2025
- 社群基础及运营 5
- 2026年秋季校园安全教育专题课 防诈骗教育
- 2026 年初中秋季开学第一课厄尔尼诺引发暴雨洪涝科普教育
- 2026年感染性疾病科细菌性感染护理培训
- 2026年秋道德与法治一年级上册教学工作计划
- 2022 低环境温度空气源多联式热泵(空调)机组
- 2023版MRI临床应用安全专家共识课件
- 2024年陕西工业职业技术学院教师招聘笔试真题
- 2011桂林奥林匹克花园项目发展战略及2012年推广报告2
- 钢结构厂房的施工方案
- 公司售电业务管理制度
- JBT 7784-2024 隐极同步发电机用交流励磁机 技术规范(正式版)
- 《风电场工程规划报告编制规程》(NB-T 31098-2016)
- 园林绿化修剪培训课件
- 人教版七年级数学下册尖子生培优必刷题专题5.8平行线的性质与判定大题专项提升训练(基础篇重难点培优30题)(原卷版+解析)
- 山东省广播电视有线网络安全播出规章制度样本
评论
0/150
提交评论