论文部分内容阅读
钱塘江水源地咸水入侵频繁,氯度超标严重,经常造成沿岸城市的水质性缺水。为确保杭州市用水安全,研究钱塘江咸水入侵规律,预报不同径流和潮汐条件下重要取水口水域的盐度变化,为科学调度上游水库下泄流量和水厂取水提供依据。根据研究需要设钱塘江咸水入侵实时数据监测系统,对闻家堰、闸口、七堡、仓前、盐官、澉浦、乍浦等地的氯度以及潮位进行全天候的实时监控,由此获得的实测数据给进一步的研究工作提供数据支撑。但由于野外环境复杂多变及现有仪器存在一定的缺陷,使得氯度实时监测数据经常出现无规律的缺失或者异常,对后期的研究工作造成很大影响,因此根据实时监测的需要开发相关程序。分析了氯度实时监测异常点及其原因,对不同特点的异常数据作了科学的分类处理。研究了氯度变化的时域特征,据此对数据进行高通滤波和限幅滤波等预处理,初步判断出有效数据和总体质量。通过短时离散傅里叶变换加汉宁窗的方法分析了氯度的频域特征,研究了观测站点氨度与天文潮的相关关系,据此设计滤波函数作频域滤波。然后用通过S-G方法判别观察点是否异常,并根据异常点分布特点诊断数据质量并报警,最后在保证合理的原始数据不失真的同时尝试通过最小二乘拟合作数据修复。研究氯度的频谱特征发现,氯度变化受天文潮影响,盐官以下受潮流影响大的地区特别明显。氯度的频域能量集中分布在全日、半日、四分之一日和六分之一日潮周期处,受多个分潮控制,其中K1、J1、S2、2SM2、S4和S6分潮的影响最大。在全日和半日潮周期的频段,氯度和天文潮的频域能量分布比较相似,在三分之一日处氯度频域能量最集中处的频率值要比分潮大1°/h左右,在四分之一日、五分之一日和六分之一日处氯度频域能量最集中处频率值要比分潮大20/h左右。本文使用C语言编写源代码,开发的程序异常点识别准确率高,运算速度快,对钱塘江仓前下游测站氯度异常数据识别率在99.5%以上,报警延迟在一个半小时以内。