九龙江污染物通量估算.docxVIP

九龙江污染物通量估算 河口海岸的海水污染是陆地和海洋相互作用的研究内容之一，也是全球环境变化的研究重点。海洋中的水质和污染物的运输能力已成为研究的重点。我国是一个河流水系相对发达的国家,但是我国在海陆相互作用方而的研究还很少,仅在少数河流进行了相关的研究。 关于河流污染物通量的估算研究,尽管有许多种估算方法可用来进行估算,但由于不同估算方法计算出来的污染物年通量值相差悬殊,应该非常谨慎地选定估算方法。河流污染物入海通量通常以年为时段进行估算,与断面瞬间通量相比,估算时段通量的估算误差相对较为困难,时段的跨度越大,时段通量波动变化的方式也就越复杂,而用有限的瞬间实测数据进行估算误差也就越大。通量估算的实测误差来源包括:测流误差、水质采样误差、水质分析误差、断面离散采样的代表性不强、采样频率带来的误差等。估算流域河流污染物入海年时段通量,其误差主要由两点引起:点源与非点源差异引起的水质监测数据代表性不强;水文水质数据不同步引起的时均离散误差。针对这些存在的问题,Webb等人采用分时段通量和、时段平均浓度与时段水量之积两类估算方法构造了5种时段通量的计算公式,国内也有学者应用这些方法进行河流污染物通量估算。本文根据九龙江下游的水质水文监测数据,针对九龙江流域水污染的特点,采用上述合适的计算方法,进行九龙江污染物入海通量的初步估算。 1 九龙河大海流量分析 1.1 农业非点源污染 九龙江位于福建省南部,是福建省仅次于闽江的第二大河流,由北溪和西溪两大支流构成,以北溪为主流。其中北溪流域面积9 803 km2,干流河长285 km,西溪流域面积3 964 km2,干流河长172 km,两溪在龙海市境内长洲附近汇合流入厦门附近海域。九龙江流经福建省龙岩和漳州境内8个市县区,工业点源污染影响较大,同时农业非点源污染问题突出。九龙江河口和厦门附近海域的海水水质以及富营养化等问题受九龙江来水水质的影响较大。上世纪末,由于流域经济快速发展、人口急剧膨胀、城镇化进程加快等原因,流域水污染问题严重,从1999年开始,福建省实施了“九龙江流域水污染与生态破坏综合整治方案”,目前已经遏制了水污染问题的进一步恶化。 1.2 计算量数据,满足计算要求 九龙江北溪和西溪下游分别设有浦南和郑店水文监测站,其径流量可近似代表九龙江的入海水量。获取的可使用水文数据包括近10 a(1995~2004年)的逐日流量数据,可满足各种计算要求。九龙江北溪和西溪靠近下游入海口处分别设有江东桥和水头两个水质监测省控断面,水质监测数据始于1996年,代表性的水质监测数据有高锰酸钾指数(CODMn)、NH3-N等,但NH3-N仅从2000年开始监测。从1996~2001年,每年仅有代表丰、平、枯3个水期的水质监测数据,其后每逢单月进行水质监测,即每年有6个水质监测数据。由此可见,九龙江污染物入海通量估算的主要问题是水质监测数据的不完善,因此选择合适的通量计算方法显得尤为重要。 2 通量计算 根据九龙江的水文水质监测数据,分别采用目前流行的河流污染物通量5种估算方法计算其入海污染物的通量。河流入海污染物通量估算一般是指以年为时段的通量估算,河流污染物通量估算的方法并不确定,通常都会存在一定的误差。5种估算方法及其应用特点分析见公式(1)~(5)所示。各公式中符号含义:W为估算时间段的污染物通量;K为不同估算时间段的转换系数;n为估算时间段内的采样次数;ci为样品i的浓度值。 W=ΚnΣi=1cinnΣi=1Qin(1)W=KΣi=1ncinΣi=1nQin(1) 公式(1)中通量采用瞬时浓度ci的平均值与瞬时流量Qi的平均值的乘积进行通量计算,该方法忽略径流量时均变化的影响; W=Κ[nΣi=1cin]ˉQ(2)W=K[Σi=1ncin]Qˉˉˉ(2) 公式(2)中通量采用瞬时浓度ci的平均值与时段平均流量 的乘积进行通量计算,该方法突出径流量时均变化的影响; W=ΚnΣi=1cQin(3)W=KΣi=1ncQin(3) 公式(3)中通量采用求瞬时浓度ci与瞬时流量Qi乘积的平均值进行通量计算,该方法忽略径流量时均变化的影响,但突出了点源污染; W=ΚnΣi=1ciˉQ(4)W=KΣi=1nciQˉˉˉ(4) 公式(4)中通量采用瞬时浓度ci与代表时段平均流量ˉQQˉˉˉ乘积的平均值进行通量计算,该方法突出径流量时均变化的影响,同时也突出了非点源污染; W=ΚnΣi=1ciQinΣi=1QiˉQ(5)W=KΣi=1nciQiΣi=1nQiQˉˉˉ(5) 公式(5)中通量采用时段通量平均浓度与时段平均流量ˉQQˉˉˉ的乘积进行通量计算,该方法突出总径流量时均变化的影响,同时也突出了非点源污染。 3 龙江1年生时期污染物进出通量 选择水质监测数据较多的2003年对上述5种污染物通量的计