近岸波生流运动三维数值模拟及验证 被引量：16

1

作者 解鸣晓 张玮 《水科学进展》 EI CAS CSCD 北大核心 2011年第3期391-399,共9页

开发建立了近岸波生流运动三维数值计算模式。模式中,引入了三维时均剩余动量、破波表面水滚、波浪水平与垂向紊动作为主要驱动力,同时考虑了波流共同作用的底部剪切力。推导了可综合反映底坡、能量传递率和密度影响的水滚能量传输方程;... 开发建立了近岸波生流运动三维数值计算模式。模式中,引入了三维时均剩余动量、破波表面水滚、波浪水平与垂向紊动作为主要驱动力,同时考虑了波流共同作用的底部剪切力。推导了可综合反映底坡、能量传递率和密度影响的水滚能量传输方程;将Larson-Kraus的二维波浪水平紊动系数表达式拓展至三维。采用大量实测数据和文献资料测试验证了所建模式,表明所建模式可有效模拟波浪增减水、底部离岸流、沿岸流、裂流、堤后环流等不同维度的波生流现象。此外,研究也表明破波水滚效应可解释波生流峰值向岸推移的物理现象,从而在模拟中不能忽略;破波带内沿岸流速垂向较为均匀的现象与波浪附加垂向紊动有关。 展开更多

关键词 三维波生流 时均剩余动量 破波水滚 波浪紊动 模式验证

在线阅读 下载PDF 职称材料

建闸河口潮波变形概化模型研究 被引量：5

2

作者 徐雪松 窦希萍 +2 位作者 赵晓冬 王向明 张新周 《水运工程》 北大核心 2011年第6期6-10,共5页

概述潮汐河口潮波变形的普遍性及建闸前后影响变形的不同因素。建立了概化河口闸下河道物理模型,说明了模型设计的理论依据及模型布置方案。通过概化模型对建闸河口不同引河长度的潮波变形特性进行模拟研究,得到不同闸址建闸后潮波变形... 概述潮汐河口潮波变形的普遍性及建闸前后影响变形的不同因素。建立了概化河口闸下河道物理模型,说明了模型设计的理论依据及模型布置方案。通过概化模型对建闸河口不同引河长度的潮波变形特性进行模拟研究,得到不同闸址建闸后潮波变形的共性和差异。闸下引河长度对潮波变形的影响明显,随着闸址距河口距离的增加,沿程潮差的增减趋势发生变化;建闸后涨潮平均流速与落潮平均流速的比值变化趋势与引河长度紧密相关。 展开更多

关键词 建闸河口 概化物理模型 引河长度 潮波变形

在线阅读 下载PDF 职称材料

近岸波生环流的三维数值模拟研究 被引量：3

3

作者 王平 张宁川 《哈尔滨工程大学学报》 EI CAS CSCD 北大核心 2015年第1期34-40,共7页

为研究近岸波浪破碎引起的复杂波生环流问题,本文基于三维水动力模型,引入随水深变化的辐射应力项,同时考虑波浪对水体紊动影响以及波流联合作用下的底部切应力,得到了波生流三维数值模型;模型基于非结构化网格拟合复杂的岸线变化,同时... 为研究近岸波浪破碎引起的复杂波生环流问题,本文基于三维水动力模型,引入随水深变化的辐射应力项,同时考虑波浪对水体紊动影响以及波流联合作用下的底部切应力,得到了波生流三维数值模型;模型基于非结构化网格拟合复杂的岸线变化,同时采用有限体积法进行数值离散。模拟了波浪浅化和破碎引起的斜坡地形上的垂向环流、人工岛后的平面环流以及波状地形下的裂流和沿岸流现象。结果表明:波致辐射应力的二维和深度分布不均会驱动平面和垂向的波生环流;同时与实验的比较,验证了本模式在计算近岸波生三维环流上的精确性和有效性。 展开更多

关键词 波生环流 三维数值模拟 非结构化网格 辐射应力 波浪紊动

在线阅读 下载PDF 职称材料

基于小波分析的珠江河口区盐度与流速、潮位的相关分析 被引量：2

4

作者 肖莞生 陈子燊 《热带地理》 北大核心 2013年第1期34-39,共6页

利用小波分析方法研究了珠江河口区盐度、流速与潮位变化的周期特征，进而利用交叉小波分析盐度与流速、潮位变化周期之间的相位关系。结果表明：珠江河口区盐度、流速及潮位变化存在着12．4h的半日周期和22．4h的日周期变化；其中盐度... 利用小波分析方法研究了珠江河口区盐度、流速与潮位变化的周期特征，进而利用交叉小波分析盐度与流速、潮位变化周期之间的相位关系。结果表明：珠江河口区盐度、流速及潮位变化存在着12．4h的半日周期和22．4h的日周期变化；其中盐度与流速在变化周期上存在显著相关关系，位相差约为85°，各站点超前和滞后特征各不相同；盐度与潮位在变化周期上存在显著相关关系，潮位变化超前于盐度变化约30°-60°。 展开更多

关键词 交叉谱 小波分析 盐度 流速 潮位 珠江河口

在线阅读 下载PDF 职称材料

乐清湾水域纳潮量演变分析 被引量：7

5

作者 王诚超 潘国富 +1 位作者 许雪峰 陈培雄 《海洋科学》 CAS CSCD 北大核心 2017年第8期76-85,共10页

基于1965年以来的不同时期水深地形数据和卫星遥感影像资料,对乐清湾岸线和不同特征值水深所围水域面积的历史变化进行研究,根据GIS技术计算乐清湾近50年的不同水域面积的演变特征。利用实测地形资料和水文数据建立乐清湾海域二维潮流... 基于1965年以来的不同时期水深地形数据和卫星遥感影像资料,对乐清湾岸线和不同特征值水深所围水域面积的历史变化进行研究,根据GIS技术计算乐清湾近50年的不同水域面积的演变特征。利用实测地形资料和水文数据建立乐清湾海域二维潮流数学模型分析乐清湾水域纳潮量的演变情况。结果表明:(1)乐清湾海域纳潮量近50年减少3.16×10~8 m^3,2013年较1965年减少17.69%,年均递减速率由0.06×10~8 m^3/a增加至近年的1.96×10~8 m^3/a;(2)海湾不同区域围填海造成相应海区水域面积的缩减,但纳潮量对水域面积改变的响应程度有显著差别。内湾滩涂围垦对乐清湾不同湾区纳潮量影响十分有限,外湾围垦对纳潮量的影响由外湾向内湾明显递减。漩门湾二期工程显著改变乐清湾的潮流形态,造成各个湾区纳潮量均出现大幅度的调整。本研究可以为海湾的生态环境保护和预测海湾的发展趋势提供量化的参考数据。 展开更多

关键词 乐清湾 水域面积 纳潮量 演变分析

在线阅读 下载PDF 职称材料

柳州市柳江河道景观工程对河道水流的影响研究 被引量：8

6

作者 方神光 黄胜伟 《中国农村水利水电》 北大核心 2010年第10期55-58,共4页

针对目前柳州市柳江河道上拟建景观对河道的影响,分别建立了柳州市柳江河道大范围河段和拟建景观工程局部范围曲线坐标系下的二维水流数学模型,并采用纯隐格式的混合有限分析法用来离散和求解该数学模型。建立的大范围的柳州市柳江河道... 针对目前柳州市柳江河道上拟建景观对河道的影响,分别建立了柳州市柳江河道大范围河段和拟建景观工程局部范围曲线坐标系下的二维水流数学模型,并采用纯隐格式的混合有限分析法用来离散和求解该数学模型。建立的大范围的柳州市柳江河道二维数学模型计算成果一方面用来进行验证,另一方面为景观工程局部水流数学模型提供初始和边界条件。结果分析表明,与工程前相比,景观工程局部河道地形有挖有填,对局部水位和流态存在一定的影响;由于对工程区堤防外的房屋进行了拆迁,一定程度改善了所处位置河道左岸的过流条件。 展开更多

关键词 混合有限分析法 二维水流数学模型 景观工程 壅水 柳江

在线阅读 下载PDF 职称材料

珠江蕉门分汊河口分流比水槽实验研究 被引量：4

7

作者 姚鹏 余志斌 +2 位作者 苏敏 安欣禧 周曾 《海洋学报》 CAS CSCD 北大核心 2022年第6期106-115,共10页

作为珠江主要的泄洪排沙口门,珠江蕉门主要承泄西北江的径流和泥沙,与虎门通过横向汊道凫洲水道连通,其分水分沙的比例对粤港澳大湾区防洪及通航安全至关重要。针对凫洲水道分流比问题,本文基于蕉门分汊河口的形态特点及地貌特征,设计... 作为珠江主要的泄洪排沙口门,珠江蕉门主要承泄西北江的径流和泥沙,与虎门通过横向汊道凫洲水道连通,其分水分沙的比例对粤港澳大湾区防洪及通航安全至关重要。针对凫洲水道分流比问题,本文基于蕉门分汊河口的形态特点及地貌特征,设计系列水槽实验,探究凫洲水道分流比、分流比影响因子及其变化趋势等。结果表明:在保证宽度比、深度比和流量比均相似的情况下,该水槽实验可较高精度地模拟蕉门分汊河口分流比。凫洲水道分流比受到蕉门与虎门共同影响:虎门与蕉门的相对流速比越大,凫洲水道分流比越大,其主汊地位越明显。在现有的地貌特征条件下,凫洲水道主支汊转变的阈值是虎门与蕉门相对流速比,约为0.35。研究成果可为蕉门整治、保障其行洪通航安全等提供参考依据。 展开更多

关键词 蕉门 凫洲水道 分汊河口 分流比 水槽实验

在线阅读 下载PDF 职称材料

伶仃洋河口水域纳潮特性分析 被引量：4

8

作者 方神光 陈文龙 崔丽琴 《海洋科学》 CAS CSCD 北大核心 2012年第6期105-112,共8页

为阐明伶仃洋河口水域纳潮现状和特性,分别采用经验公式和建立贴体曲线坐标系下的珠江河口二维数学模型的方法,对河口伶仃洋水域的纳潮特性进行了分析和探讨。数学模型的离散和求解采用纯隐格式的混合有限分析法,根据实测资料对其进行... 为阐明伶仃洋河口水域纳潮现状和特性,分别采用经验公式和建立贴体曲线坐标系下的珠江河口二维数学模型的方法,对河口伶仃洋水域的纳潮特性进行了分析和探讨。数学模型的离散和求解采用纯隐格式的混合有限分析法,根据实测资料对其进行了验证。比较显示,采用实测资料和经验公式计算到的伶仃洋水域纳潮量与数学模型的计算结果相差在5%以内,从而证实了数学模型的可靠性。结果表明,伶仃洋海域的主要纳潮口门为伶仃洋口门,约占总纳潮量的87.7%,通过香港暗士顿水道进出的潮量仅占12.3%左右,显示伶仃洋口门为该海域的主要径潮通道,对该河口的防洪和生态具有重要的作用。 展开更多

关键词 伶仃洋 纳潮量 经验公式 数学模型

在线阅读 下载PDF 职称材料

咸潮上溯期磨刀门河口垂向环流动力机制 被引量：3

9

作者 彭石 何用 +1 位作者 刘诚 邹华志 《水动力学研究与进展（A辑）》 CSCD 北大核心 2018年第5期650-657,共8页

为研究垂向环流在磨刀门咸潮上溯过程中的作用,该文引入淡水输运速度分析潮周期平均水流特征,获得了河口垂向环流结构在咸潮半月周期波动过程中的时空变化。通过改进环流强度与密度弗劳德数的计算,分析了两者的内在关联,且探讨了径流量... 为研究垂向环流在磨刀门咸潮上溯过程中的作用,该文引入淡水输运速度分析潮周期平均水流特征,获得了河口垂向环流结构在咸潮半月周期波动过程中的时空变化。通过改进环流强度与密度弗劳德数的计算,分析了两者的内在关联,且探讨了径流量、垂向盐度分层和地形等动力因素影响垂向环流的机理。结果表明,小潮期间垂向环流强于大潮期间,环流强度的增减与盐水头部(2‰等盐度线)的上移和回落同步;垂向环流形成的动力学条件为密度弗劳德数大于等于1;垂线平均淡水输运速度代表的惯性力及表底层盐度差和水深代表的密度梯度力,对磨刀门水道垂向环流结构起决定性作用;这两种作用力的相对优势转换,是引起咸潮上溯周期变化的主要原因。 展开更多

关键词 垂向环流 咸潮 径流 密度弗劳德数 磨刀门

原文传递

闸控潮汐河网区水环境调度模型研究 被引量：7

10

作者 贺新春 王翠婷 +1 位作者 汝向文 丁波 《华北水利水电大学学报（自然科学版）》 2018年第2期86-92,共7页

针对闸控潮汐河网区的水系条件和水动力特征,立足于水网区水环境改善的需求,构建了闸控潮汐河网区水环境调度模型。该模型由河网区一维水动力模型、河网一维水质模型和水闸调控模拟模型3个模型构成,采用实测资料对模型进行了率定和验证... 针对闸控潮汐河网区的水系条件和水动力特征,立足于水网区水环境改善的需求,构建了闸控潮汐河网区水环境调度模型。该模型由河网区一维水动力模型、河网一维水质模型和水闸调控模拟模型3个模型构成,采用实测资料对模型进行了率定和验证。选择珠江三角洲"中珠联围"作为研究区开展实证研究,提出了联围整体调度与局部重点改善相结合的调度方法。结果表明,统筹考虑中山和珠海的调度需求,实施"中珠联围"统一调度,形成"西北进、东南出"的水流格局,可有效地加快联围内河涌的水体更新速度,降低污染物平均浓度,对于"中珠联围"水环境的改善及供水水质安全的保障具有重要意义。 展开更多

关键词 潮汐河网 闸群联控 水环境 调度模型 “中珠联围”

在线阅读 下载PDF 职称材料

混合有限分析法在桥墩水流数值模拟中的应用 被引量：2

11

作者 方神光 黄胜伟 崔丽琴 《人民长江》 北大核心 2010年第7期75-78,共4页

在河道管理范围内布置桥墩时,较为关注桥墩对壅水和流态的影响,以柳州市局部柳江河道和拟建的维义大桥为例,分别建立大范围的河道和桥墩局部二维水动力数学模型,采用混合有限分析法进行了计算和求解。结果分析表明,维义桥的建设使局部... 在河道管理范围内布置桥墩时,较为关注桥墩对壅水和流态的影响,以柳州市局部柳江河道和拟建的维义大桥为例,分别建立大范围的河道和桥墩局部二维水动力数学模型,采用混合有限分析法进行了计算和求解。结果分析表明,维义桥的建设使局部河道产生壅水现象并对河势存在一定的影响。同时也证实了该方法在工程实际应用中的准确性和有效性。 展开更多

关键词 混合有限分析法 二维水流数学模型 桥墩 壅水

在线阅读 下载PDF 职称材料

澳门地区东侧水域夏季余流及通量的时空变化特征

12

作者 喻丰华 方神光 +1 位作者 何用 何青 《华东师范大学学报（自然科学版）》 CAS CSCD 北大核心 2023年第3期147-157,共11页

基于澳门地区东侧水域两座河口原型观测试验站夏季15 d的连续观测资料,分析了潮周期垂向平均欧拉余流、拉格朗日余流和斯托克斯余流时空变化特征的同时,将其概念推广到垂向潮汐动力影响下的分层余流分析中.研究显示,夏季澳门地区东侧南... 基于澳门地区东侧水域两座河口原型观测试验站夏季15 d的连续观测资料,分析了潮周期垂向平均欧拉余流、拉格朗日余流和斯托克斯余流时空变化特征的同时,将其概念推广到垂向潮汐动力影响下的分层余流分析中.研究显示,夏季澳门地区东侧南、北端水域整体涨落潮流态一致,南侧水域潮动力显著大于北侧;半月时段内,北侧水域潮流物质整体输移方向指向西北,南侧指向东南,拉格朗日余流流速分别为2.2 cm/s和5.1 cm/s,略小于欧拉余流.受海面西南季风影响,南、北侧水域表层欧拉余流、拉格朗日余流和斯托克斯余流流向均为东北向.潮周期垂向平均余流在洪水期指向外海侧,非洪水期北侧水域指向近岸,南侧水域指向伶仃洋河口东侧.海面风主导东南侧水域的斯托克斯漂移强度,对东北侧水域影响不明显,洪水期径流动力增强会削弱海面风引起的漂移作用;但南、北侧水域表层潮流物质输运强度和方向都与海面风密切相关.基于余流流态结合净潮通量分析显示,夏季澳门地区水道出口以下水域存在较为稳定的逆时针余环流,初步分析其原因主要是:珠江河口外海侧强劲的东北向沿岸流引起澳门地区东侧水域形成自东向西的补偿流所致.夏季,该独特的动力结构会截获上游东四口门下泄的部分水沙进入澳门地区水道,导致澳门地区水道滩槽淤积和水体交换不畅. 展开更多

关键词 余流 净潮通量 潮动力 海面风 澳门地区水域

在线阅读 下载PDF 职称材料

混合有限分析法在桥墩综合壅水叠加模拟中的应用 被引量：1

13

作者 方神光 余顺超 黄胜伟 《长沙理工大学学报（自然科学版）》 CAS 2010年第2期53-57,共5页

将一种纯隐格式-混合有限分析法推广到实际工程应用中.采用该方法对曲线坐标系下的二维水流数学模型进行了离散,将常用的干湿节点判别标准与冻结法结合起来,有效提高了模拟计算的效率和稳定性.采用建立的二维水流数学模型和该计算方法,... 将一种纯隐格式-混合有限分析法推广到实际工程应用中.采用该方法对曲线坐标系下的二维水流数学模型进行了离散,将常用的干湿节点判别标准与冻结法结合起来,有效提高了模拟计算的效率和稳定性.采用建立的二维水流数学模型和该计算方法,有效地模拟和重现了柳州市柳江河道上最近发生的洪水.计算和分析显示,柳州市柳江河道上当前已有桥梁由于间距太近,形成了显著的壅水叠加现象.采用经验公式计算了各单个桥梁的壅水高度,与数值模拟计算结果较为吻合. 展开更多

关键词 混合有限分析法 二维水流数学模型 桥墩 阻力概化

在线阅读 下载PDF 职称材料

低洼地区农田淹没成因分析方法及应用(Ⅱ)——洪水反演和分析 被引量：1

14

作者 姚志坚 方神光 陈文龙 《中国农村水利水电》 北大核心 2012年第3期154-156,共3页

为探求从化市水南村部分农田被淹没的主要原因,对水南村所在小流域的降雨资料和涌容关系进行了分析。在对流溪河洪水特性和牛心岭电站操作对淹没位置河段水位影响分析的基础上,真实反演出了农田被淹没的过程。结果显示,在不考虑淹没区... 为探求从化市水南村部分农田被淹没的主要原因,对水南村所在小流域的降雨资料和涌容关系进行了分析。在对流溪河洪水特性和牛心岭电站操作对淹没位置河段水位影响分析的基础上,真实反演出了农田被淹没的过程。结果显示,在不考虑淹没区涵洞排水的条件下,采用大坳站和太平场站的平均雨量资料,本次降雨造成淹没区的水位上涨幅度可达2.74m;淹没区水位经历了两次上涨,第1次主要由产汇流太急造成涵洞排水不及造成,第2次主要是流溪河洪水上涨造成涵洞无法排水所致。未封堵暗渠导致了流溪河洪水倒灌入淹没区,造成最大淹没深度为0.23m;牛心岭电站操作对淹没区影响较小。 展开更多

关键词 二维水流数学模型 混合有限分析法 洪水反演 淹没

在线阅读 下载PDF 职称材料

港珠澳大桥对伶仃洋河口水域纳潮影响分析 被引量：7

15

作者 方神光 陈文龙 《人民珠江》 2011年第1期18-22,72,共6页

为分析港珠澳大桥对珠江河口水域纳潮量的影响,分别采用经验公式和建立贴体曲线坐标系下的珠江河口二维数学模型的方法,对河口伶仃洋水域的纳潮特性进行了计算和分析。数学模型的离散和求解采用纯隐格式的混合有限分析法,根据实测资料... 为分析港珠澳大桥对珠江河口水域纳潮量的影响,分别采用经验公式和建立贴体曲线坐标系下的珠江河口二维数学模型的方法,对河口伶仃洋水域的纳潮特性进行了计算和分析。数学模型的离散和求解采用纯隐格式的混合有限分析法,根据实测资料对其进行了验证。比较显示,采用实测资料和经验公式计算到的伶仃洋水域纳潮量与数学模型的计算结果相差在5%以内,从而证实了数学模型的可靠性。 展开更多

关键词 伶仃洋 数学模型 纳潮量 经验公式 港珠澳大桥

在线阅读 下载PDF 职称材料

惠阳淡水河流域年底前关停11家污染企业

16

《给水排水》 CSCD 北大核心 2008年第11期55-55,共1页

关键词 重污染企业 淡水河 河流域 整治工作 流域污染 污染治理 惠州市

在线阅读 下载PDF 职称材料

ECOMSED模型动边界处理技术研究 被引量：2

17

作者 赵海亮 张艳飞 《水道港口》 2010年第5期395-398,共4页

在ECOMSED数值模式的水动力模块基础上,引入了动边界处理技术,建立了数值试验,并采用单一斜坡地形和圆形浅滩斜坡地形对该技术进行了验证。数值试验结果显示,改进后的水动力计算模式可成功再现由涨落潮引起的潮间带漫滩现象,而且模拟得... 在ECOMSED数值模式的水动力模块基础上,引入了动边界处理技术,建立了数值试验,并采用单一斜坡地形和圆形浅滩斜坡地形对该技术进行了验证。数值试验结果显示,改进后的水动力计算模式可成功再现由涨落潮引起的潮间带漫滩现象,而且模拟得到的水面线较为光滑,且未产生零水深和负水深。 展开更多

关键词 ECOMSED 动边界 数值模拟

在线阅读 下载PDF 职称材料

低洼地区农田淹没成因分析方法及应用——(I)数学模型验证及洪水特性分析

18

作者 方神光 姚志坚 陈文龙 《中国农村水利水电》 北大核心 2012年第2期153-155,共3页

为探求从化市水南村部分农田被淹没的原因,分别采用现场调查、水文地形资料分析、数学模型、调洪演算等方法,对本次流溪河洪水和下游牛心岭电站对淹没位置河段水位的影响进行了分析和探讨。结果表明,洪水从大拗坝到达水南村淹没位置河... 为探求从化市水南村部分农田被淹没的原因,分别采用现场调查、水文地形资料分析、数学模型、调洪演算等方法,对本次流溪河洪水和下游牛心岭电站对淹没位置河段水位的影响进行了分析和探讨。结果表明,洪水从大拗坝到达水南村淹没位置河道的时间约为2.6h;从5月7日凌晨4时至12时,流溪河水位明显高于淹没区地形,极容易导致外江洪水从未封堵暗渠进入淹没区;牛心岭电站只用4个闸孔泄洪,在洪峰时刻引起淹没区外侧河道水位抬高约0.05m。 展开更多

关键词 二维水流数学模型 混合有限分析法 洪水 农田淹没

在线阅读 下载PDF 职称材料

南粤绿讯

19

《环境》 2013年第3期60-61,共2页

省人大常委会召开关于深化淡水河、石马河流域整治工作座谈会 2月6日上午，关于深化淡水河、石马河流域整治工作座谈会在东莞召开，省人大常委会主任黄龙云、省人大常委会副主任肖志恒、雷于蓝、陈小川、秘书长陈逸葵，以及省环保厅厅... 省人大常委会召开关于深化淡水河、石马河流域整治工作座谈会 2月6日上午，关于深化淡水河、石马河流域整治工作座谈会在东莞召开，省人大常委会主任黄龙云、省人大常委会副主任肖志恒、雷于蓝、陈小川、秘书长陈逸葵，以及省环保厅厅长李清、副厅长陈光荣，省人大环资委、省国土厅、省水利厅和深圳、惠州、东莞市相关负责人出席会议。 展开更多

关键词 人大常委会 整治工作 淡水河 河流域 副主任 东莞市 秘书长 水利厅

在线阅读 下载PDF 职称材料

“地毯式”核查淡水河和观澜河污染

20

《环境》 2009年第7期43-44,共2页

6月10—15日，广东省环保局组织120余名工作人员，分成13个小组，对淡水河和观澜河流域重污染行业进行了“地毯式”核查。行动除按照污染源普查工作提供的1970家企业名单逐一开展核查外，还对近期新上建设项目以及违法建设项目进行突击... 6月10—15日，广东省环保局组织120余名工作人员，分成13个小组，对淡水河和观澜河流域重污染行业进行了“地毯式”核查。行动除按照污染源普查工作提供的1970家企业名单逐一开展核查外，还对近期新上建设项目以及违法建设项目进行突击检查，重点检查电镀、印染、化工、制革、畜禽养殖等重污染企业。 展开更多

关键词 污染行业 淡水河 地毯式 核查 观澜河 建设项目 重污染企业 工作人员

在线阅读 下载PDF 职称材料