悬浮物输运运动规律研究.docxVIP

悬浮物输运运动规律研究 悬浮物是水环境研究的重要参数，主要由沉积物、藻类、微生物和动植物碎屑组成。悬浮物的存在直接影响了水体透明度和溶解氧,影响大型沉水植物生长及鱼类等水生动物的生存,使水质变差。 镇江是长江三角洲典型的滨江城市,位于长江与京杭运河2条黄金水道交汇之处。镇江内湖位于长江镇扬河段中段的世业洲汊道下段、六圩弯道上段右岸征润洲边滩内,水面面积约为8.8km2,是镇江的主要水体。内湖上游通过引航道与长江相通,下游在焦山东侧与长江汇聚,在潮位变化下与长江进行水体交换。由于长江水体中大量泥沙的涌入,致使目前内湖水体悬浮物含量较高,透明度只有15～18cm,水体明显浑浊、泛黄、严重影响了内湖水生生态系统的平衡,使水质进一步恶化。为研究内湖悬浮物的输运扩散及运动规律,建立了二维悬浮物输移的数学模型。模型考虑了风场及潮流作用下悬浮物的输运,同时还考虑了悬浮物的沉降和再悬浮机制,采用切应力概念确定悬浮物的起悬量和沉降量。该模型基于无结构网格有限体积法,能较好地拟合计算区域的复杂地形和边界;采用高性能的Osher格式计算水量、动量及悬浮物输运等通量,提高了数值模拟的精度。同时,2次野外的水文和悬浮物含量的同步监测为模型的验证提供了详细准确的基础资料。 1 内部的同步监控 1.1 潮洪时期掳航道 根据长江镇江段近期潮汐变化情况,确定2003年4月21日(农历2月3日,大潮汛)和2003年8月18日(农历7月3日,大潮汛)的全日潮为枯季和洪季的2次监测时间。 在引航道上、下断面,运粮河,焦南闸布设了4个监测断面,每个断面均设3条垂线;在内湖主体水域布设了7个监测点,每点设一条垂线,同步监测水文、含沙量、悬浮物浓度等项目,具体点位见图1。 1.2 潮前后潮水文化观测 (1)水位观测。在引航道上、下断面、运粮河、焦南闸4处断面进行水位观测,其中焦南闸设置自记水位记实施观测,引航道上、下断面、运粮河采用人工现场观测即每1h观测1次,高(低)潮前后加密观测,时间为26h(同步观测)。(2)流速流向监测。在各垂线处进行定点流速流向测量。垂线定点流速、流向测验运用SL25-1流速流向仪,涨潮每逢整点小时前后施测1次,落潮每2h施测1次。各条垂线全潮定点流速、流向测验均同步进行。(3)悬浮物浓度监测。引航道上、下断面、运粮河、焦南闸4处断面,每个断面设3个采样点,内湖7个监测点各设一个采样点。各采样点均采集表层(水面下50cm)水样,样品采集后按照规范要求分批及时送回实验室,根据《水和废水监测分析方法》及《湖泊富营养化调查规范》(第二版)分析测定悬浮物浓度。 1.3 引航道断面和水域流速 2次监测的潮位特征值见表1,各垂线全潮流速特征值见表2,内湖水域垂线悬浮颗粒含量见表3。 由上表可知,内湖潮位日变化明显,每日出现2次高低潮。枯季(4月21日)潮差较大,实测最大潮差1.55m(运粮河断面),最小潮差1.21m(引航道上、下断面),最大涨潮历时3h 45min(引航道下断面),最大落潮历时9h20min(焦南闸断面)。洪季(8月18日)实测最大潮差1.16m(引航道下断面),最小潮差0.71m(运粮河断面),最大涨潮历时3h30min(引航道上、下断面),最大落潮历时9h40min(焦南闸断面)。洪季内湖水域流速普遍大于枯季流速,靠近深泓的14#、17#、18#、19#流速较大,枯季一般为0.1m/s～0.15m/s左右,洪季为0.3m/s～0.7m/s左右;浅滩处的16#流速较小,低于0.1m/s。由于枯季内湖有污染源排入,故枯季内湖水域悬浮物浓度大于洪季,且污染源附近的15#、17#、18#、19#浓度较大。洪季内湖水域悬浮物浓度为0.045kg/m3～0.09kg/m3之间。 2 悬浮物输运的数学模型 2.1 控制方程式 2.1.1 y-sfy-hfy1 {?h?t+?(hu)?x+?(hv)?y=0?(hu)?t+?(hu2+gh2/2)?x+?(huv)?y=gh(s0x-sfx)+hfv+hFx?(hv)?t+?(huv)?x+?(hv2+gh2/2)?y=gh(s0y-sfy)-hfu+hFy(1) 其中: {sfx=ρu√u2+v2hc2=ρn2u√u2+v2h4/3sfy=ρv√u2+v2hc2=ρn2v√u2+v2h4/3?Fx=1ρhρaCDua√u2a+v2a?Fy=1ρhρaCDva√u2a+v2a?s0y=-?Ζb?y?s0x=-?Ζb?x(2) 式中:sfx为x向的摩阻底坡;sfy为y向的摩阻底坡;s0x为x向的河底底坡;s0y为y向的河底底坡;Fx为摩擦力在x方向上的分量;Fy为摩擦力在y方向上的分量;风应力就是通过Fx、Fy而起作用;ρ、ρa分别为水和空气的密度;h为水深;CD为风拖曳系数;ua、va分别为风速在x