长江口粗粉砂和极细砂输移特性研究

1、长江口粗粉砂和极细砂输移特性研究摘要：在长江口区主航槽中利用悬浮物浓度测定系统(AP1)第一次获得了近底紊流“猝发(burst)和床面泥沙喷发(ejetin)的连续时间序列可视图像和数据，同时利用高分辨率流速仪以及旁侧声纳、热敏式双频测深仪等手段获得了大量底形沙波的实测图像和数据。这些图像和数据结果说明：长江口南支南港之间粗粉砂至极细砂在涨落潮流的作用下，呈现静止细颗粒泥沙起、悬扬粗颗粒喷发运移的同时进展大尺度底形运动细颗粒起扬、悬扬静止的输移循环形式，即当U50/s，泥沙在底床上根本不运动；50/sU65/s，小尺度紊流猝发发生，泥沙开场小尺度起扬，并有幕式猝发现象；65/sU100/s时，

2、小尺度紊流猝发强度和持续宽度急剧增大，床面细颗粒泥沙开场喷发导致床面开场冲淤，形成中小尺度沙波；100/sU110/s时，大尺度紊流形成，床面粗颗粒泥沙开场喷发；到了U110/s时，宏观紊流规模更大，床面粗颗粒泥沙的喷发强度和持续宽度急剧增大，导致床面形成大尺度沙波运动。由于U110/s的落潮流速远远大于涨潮流速，且落潮最大流速达137/s，落急(流速100/s)历时(3小时16分)长，因此发育了大量陡坡朝向下游的大尺度底形沙波，尤其是波长10的大尺度沙波都是朝向落潮流方向；而涨潮最大流速仅113/s，涨急(流速100/s)历时仅约40分钟，因此没有发育朝向上游方向的沙波。关键词：粗粉砂极细砂

3、幕式再悬浮猝发喷发底形运动输移形式1引言与波流共同作用下的砂质河口13相比，长江河口是一个悬沙浓度较高、底质分布较复杂的潮汐河口，有泥、粉砂、极细砂和浮泥、细砂等，且各河槽底床表层泥沙的平面分布存在较大差异(图1)4。中值粒径d50在0.1250.03的粗粉砂和极细砂在长江河口区分布广，如南支下游、南港、北港，北支及滩槽交界处，水深均在1015左右5。根据现代普遍采用的美国地球物理学会(AGU)泥沙分类标准6和国际沉积学会(SEP)的底形分类标准7，粒级为0.030.125的粗粉砂至极细砂由于其粒度太细，粘滞力较大，一般以冲泻质、平床方式运动，仅低流速时有波高0.03的沙纹运动，与悬浮泥沙根本

4、无交换8。甚至在KarandREiners(1978)的图解中被归到了粘性斯托克斯域(visusStkesregie)9内，亦无大尺度底形发育。图1长江口底质沉积物类型分布图DistributinfbedsedienttypesinthehangjiangEstuary周程熹等曾报道在长江口该粒径分布区洪季有大尺度底形(bedfr)沙波运动1012，因此多年来一直为人们所关注，但由于受观测方法和仪器条件的限制，迄今为止其存在的可能性仍受到疑心1315。而且粗粉砂至极细砂被认为因颗粒间的粘滞力明显增大，泥沙颗粒所受到的起动切力大，只有在足够大的涨急和落急水流条件下，细颗粒周围的粘土被悬浮后才能以

5、跳跃或半悬浮形式在近床面上运动12，16。作者于1997年12月3日至5日沿长江口南支下游至南港上游河段进展了走航式旁侧声纳、热敏式双频测深仪、Gehirp浅地层剖面仪测量，并在14小时的一个连续涨、落潮过程中进展了连续定点悬沙剖面(ASS)、END海流计、悬沙含量等新仪器的测量，获得了明晰、完好、系统的高分辨率三维底形沙波和实时同步高分辨率底流速实测资料、14小时连续时间序列的底沙变化运动全过程的数据和可视图象。这些研究成果证实了长江口粗粉砂至极细砂在洪、枯季都有底形运动的发育，对长江口泥沙的全沙输移模型、港口航道的疏浚等研究将具有重要意义，填补了国际沉积学界关于粗粉砂至极细砂粒级范围内大尺

6、度底形运动研究的空白。本文特向同行报道这一研究成果，并对长江口粗粉砂至极细砂的输移形式作一初步讨论。2研究区域长江是流域面积达1.8106k2，年平均流量293003/s，最大流量可达926003/s，最小流量46203/s，平均潮差3.12，每年向东海输送4.78108t的泥沙，其中有0.8108t的较粗颗粒的泥沙沉积在河口段7，与外海输入泥沙一起组成河口区底沙(bedaterial)。长江出徐六泾节点之后，首先被崇明岛分隔为南北两支，继而南支又被长兴、横沙岛分隔为南北两港，节点到口门的间隔 和口门的宽度相当，均为90k。长江口的平均粒径为0.160.002，中值粒径d50=0.170.00

7、2，即从细砂到粘土均有。平均而言，长江口水域沉积物较潮滩粗，水域中，口内河道较口外海滨粗；河道中，南支粗于北支，北港粗于南港。河道呈现从上游向下游变细的趋势：如北支上、下段z平均为0.06和0.045。由于水下地形滩槽交织，沉积物横向变化频繁剧烈，一般规律是自深槽向浅滩变细。考虑到南支对下游南、北港和南、北槽发育的控制性影响和前人的双频测深资料2，3，研究者选择了床面泥沙平均粒径在0.1770.031之间的南支下段至南港上段，展开了由吴淞口至下游(剖面AA、测线长18k)和由下游南港至上游南支(剖面BB、测线长21.07k)等两条剖面的实测工作(图2)，并在南支下游宝山河段选择一个实测点进展定

8、点定位测量(DZ1站)。在南支与南港两河段之间有黄浦江自吴淞入海口，吴淞年平均净径流量约3453/s，最大进潮流量可达121003/s，最大落潮流量为60003/s，最大涨潮平均流速1.8/s，最大落潮平均流速1.5/s，涨潮流时间约5小时，落潮流时间约7小时。图2长江口中的测线和测站DetetedlineandstatinfbednfiguratininthehanggjiangEstuary转贴于论文联盟.ll.论文联盟.LL.编辑。3研究方法调查时使用80t木质渔船，船速3节，所有观测仪器除了探测器外，其它设备均安装在船舱内。航线和测点定位采用美国DelNrteIn.公司的1008/58

9、6D.GPS(实时定位)。河槽床面形态的纵剖面特征由测深仪和浅地层剖面仪获得，前者为美国InnerspaeIn.公司的热敏式测深仪(449TheralDepthSunderRerder)沿两条测线(剖面AA、BB)展开，后者为英国GeAustisLiited公司的GehirpSubBttPrfilerSyste沿测线2展开，这两种仪器的探测器垂直悬挂在船右舷水面以下1处。河槽床面形态的平面分布那么利用英国UltraEletrnis公司的旁侧声纳(3050LsanSidesanSnar)进展测量，拖鱼状探测器放置在船尾水面漂福1997年12月4日至5日14小时连续定点定位测量主要有两项内容：一为

10、底沙输移随时间的变化过程；二为底沙上部水体流速流向随时间的变化过程。底沙随时间的输移变化过程采用中国科学院东海研究站设计出品的AP1声学浓度剖面仪测定，探测器(发射/接收共用换能器)与加重铅鱼垂直悬挂在船右舷水面以下6。底沙上部水体随时间的变化(流速、流向)过程那么利用美国Ende/YSIIn公司的流速仪(174SSurrernteter)进展测量，将其垂直悬挂在调查船左舷离水面以下1处；同时进展的传统流速、流向测量仪与加重铅鱼垂直悬挂在左舷床面以上1。水样的采集点为0.6H，0.8H，1.0H(H为水体深度)，采集时间间隔为30in,其将作为悬浮物的浓度和粒度分析之用，前者采用0.45定性滤

11、纸过滤法，后者采用TAulter计数器分析。底质粒度分析采用标准筛析法。4结果4.1底沙粒径组成根据粒度分析结果，本次调查区域的底沙在枯季为粉砂至细砂，以粗粉砂至极细砂为主(图3)，z=0.1250.031，部分采样点那么以细砂为主，z=0.083，且粗粉砂至极细砂中仍含有2030%的细砂。因此长江口南支南港沉积物粒度组成有较宽的粒度分布区间。表1长江口南支-南港底沙沉积物粒度参数GrainsizeparaetersfbedsedientinSuthBranhandSuthhannelinthehangjiangEstuary分级名称粒度区间()比例(%)比例(%)中砂0.50.250.050

12、.25细砂0.250.1254.0020.50极细砂0.1250.0632.6451.52粗粉砂0.060.0332.3216.17中粉砂0.030.0169.132.72细粉砂0.0160.0089.884.52极细粉砂0.0080.0045.072.11粘土0.0046.912.21图3长江口南支粗粉砂和极细砂质底沙的累计频率曲线(a)和频率曲线(b)Auulative(a)andfrequeny(b)urvesfbedsedientpsedfarsesiltandveryfinesandinSuthBranhinthehangjiangEstuary4.2底形特征本区测线水深在1017，

13、两条测线上绝大部分声纳图像均显现出河床外表呈波状起伏的韵律形态。根据双频测深图谱上的波长、波高大小和声纳图谱上沙波脊的宽度、形态、连续性等(图4、5，表2)可将本区沙波划分为以下几种类型：图4长江口南支南港粗粉砂和极细砂质分布区底形沙波的热敏式测深仪记录(测量时间：1997年12月3日至5日)其中A：小型对称和不对称沙波，波长510;B：对称沙波，波长1015；：不对称和叠置波波，前者波长1015，后者波长小于5；D：大型孤立沙波，波长大于30PTheralDepthRerderrerdsfbedfrspsedfaresesiltandveryfinesandinSuthhannelandSu

14、thBranh,hangjiangEstuaryduringDeeber35,1997图5长江口粗粉砂和极细砂质底沙大尺度底形动波的旁侧声纳记录(测量时间1997年12月3日至5日)A：直脊状小型沙波B：弯曲小型沙波：新月型沙波D：左弯曲沙波；右新月型沙波idesanSnarSidesanrerdsfbedfrspsedfarsesiltandveryfinesandinthehangjiangEstuarydurngDeeber35,1997表2长江口高分辨率底形沙波双频测深和旁侧声纳记录TheTheralDepthRerderandSnarSidesanrerdsfhighreslutin

15、bedfrsinthehangjiangEstuary名称波长()波高()波/高比陡坡(角度和方向)双频测深、声纳图像百分比(%)直脊状小型对称沙波5100.20.5253035下游直脊线状脊、谷密集均匀图4A、5A26.97直脊状小型不对称沙波5120.10.45080311下游为主直脊线状脊、谷不均匀图4A、5A34.25弯曲型对称沙波10150.30.7214057，下游为主脊弯曲或分叉，图4B、5B5.60弯曲型不对称沙波10150.20.750100712，下游为主脊线断续，图4、5B、5D左10.32新月型对称沙波15200.41.0204056，下游为主脊线断续、呈新月型图4、5

16、、5D右11.36新月型不对称沙波150.51.23050512，下游为主同上图4、5、5D右11.00孤立状沙丘301.030图4D0.54.3高分辨率流速、流向与悬浮物浓度变化过程定点定位流速、流向和ASS悬浮物浓度剖面仪测量图像和结果(图6、7和表3)显示：研究区内落潮历时大于涨潮历时约2小时，落潮流速大于涨潮流速，落潮最大流速达137/s，涨潮最大流速达113/s，落急(流速100/s)历时(3小时16分)亦长于涨急(流速100/s)历时(仅约40分钟，图6)；流速可分成以下几个区间，且长江南支南港由粗粉砂极细砂细砂组成的河槽床面形态和床沙的再悬浮等特征与流速的变化直接相关：U50/s

17、时，床面根本保持平整，无床沙起、悬扬；50U60/s时，床面稍有侵蚀和堆积地貌，床沙有掀动，起扬高程达640，仅有波长为数的沙纹；60/sU100/s时，床沙再悬浮(或喷发，ejetin)作用加强，形成大型沙纹和和小尺度沙丘，持续时间较短；100U110/s时，底沙再悬浮作用急剧加强，持续时间长且喷发高度大；U110/s时，底沙再悬浮作用很强，且能形成大尺度沙丘(图7和表3)。图6长江口南支1998年12月4日至5日连续14小时流速和流向变化过程Flutuatinfurrentvelityanddiretininntinuus14hursinSuthBranhinthehangjiangEst

18、uaryduringDeeber45,1997图7长江口粗粉砂和极细砂质底沙随涨落潮的流速变化周期发生底形运动和再悬浮输移变化过程的声学悬浮物浓度剖面仪(AP)记录(测量时间1997年12月4日至5日)A：伴随着沙纹的运动细颗粒的中等强度再悬浮(黄色)和粗颗粒的弱再悬浮(红紫色)B：伴随着小尺度沙波运动细颗粒的强再悬浮和粗颗粒的中等强度再悬浮：伴随着中等尺度沙波运动粗细颗粒均有强再悬浮D：伴随着大尺度沙波运动粗细颗粒均亦有很强的再悬浮AustiSuspendednentratinPrfilererdspsedfarsesiltandveryfinesandiththeurrentvelityf

19、lutuatinfebbandfldinthehangjiangEstuaryduringDeeber45，1997表31997年12月4至5日长江口南支南港枯季流速、流向和再悬浮变化特征(图6、7)urrentVelity,diretinandresuspensinhangeinSuthBranhandSuthhannelinthehangjiangEstuaryindryseasnduringDeeber45,1997时间流速(/s)流向床面形态特征再悬浮特征(起扬高程和持续宽度)1614163054.063.8303319少量掀沙和微量冲淤660130s1631171248.717.63

20、01321平整无171317453.319.6159311，方向经常变化大型沙波无1746180917.849.2183151平整无1810190050.9997.96137149少量掀沙、微量冲淤640115s19011906100.2107.5144148较大掀沙和冲淤2600160s1907202073.698.5139142同上同上20212337100137138151大强度掀沙和冲淤同上233800399.253.8141151中等强度掀沙和冲淤6600130s00402544.216.3106186平整无0261029.062.3203310，方向经常变化大型沙波无1032006

21、3.898.7311329少量掀沙和冲淤240115s201237100.2113.5314319较大掀沙和冲淤660130s23840452.096.2312317中等强度掀沙和冲淤2400160s40544413.649.2300332平整无4450.550168326，方向经常变化大型沙波无转贴于论文联盟.ll.5长江口粗粉砂至极细砂输移形式讨论通过上述对长江河口南支下游南港上游主槽及其附近的粗粉砂和极细砂质底床床面形态以及其上部水体运动进展三维立体观测时所获得大量可靠底沙再悬浮和大尺度底形沙波数据和图像资料研究结果，可以提出这样一个关于长江口粉砂和极细砂的输移形式(图8)，即在涨落潮流

22、的作用下，呈现静止细颗粒泥沙起、悬扬粗颗粒喷发运移的。同时进展大尺度底形运动细颗粒起扬、悬扬静止的输移循环形式，即当U50/s，泥沙在底床上根本不运动；50/sU65/s，微观紊流猝发(irturbulentburst)发生，泥沙开场小尺度起扬，并有幕式猝发现象；65/sU100s时，小尺度紊流猝发强度和持续宽度急剧增大，床面细颗粒泥沙开场喷发(ejetin)导致床面开场冲淤，形成中小尺度沙波；100/sU110/s时，大尺度紊流形成，床面粗颗粒泥沙开场喷发(ejetinrklks)；到了U110/s时，大尺度紊流规模更大，床面粗颗粒泥沙的喷发强度和持续宽度急剧增大，导致床面形成大尺度沙波运动

23、。图8长江口粗粉砂和极细砂输移概化形式neptualdelfarsesiltandveryfinesandtransprtinthehangjiangEstuary表4底流速(U)、再悬浮和底形特征(图6、7)urrentvelity(U),resuspensinandbedfrharateristis底流速(/s)再悬浮特征(浓度，持续时间和喷发高度)底形AP图像特征(图7)50无无5060少而低沙纹A60100较多而高大型沙纹和小型沙波B100110较强，持续时间长且高中等尺度沙波110140很强，长且高大尺度沙波D由于U110/s的落潮流速远远大于涨潮流速，且落潮最大流速达137/s，落

24、急(流速100/s)历时(3小时16分)长，因此发育了大量陡坡朝向下游的大尺度底形沙波，尤其是波长10的大尺度沙波都是朝向落潮流方向；而涨潮最大流速仅113/s，涨急(流速100/s)历时仅约40分钟，因此没有发育朝向上游方向的沙波。显然上述有关粗粉砂至极细砂质大尺度沙波形成和运动成因机制与目前地学界一般看法很不一致，后者认为大尺度底形沙波是由细砂级以上的泥沙颗粒作推移质运动而成，成因假说有四种：从部分扰动引起的近壁层流层的紊动特性出发1718；从床面外表作用力的平衡条件考虑床面的不稳定性和水面波的诱导19；前两者的结合20；从水流底部的紊动构造出发，结合运动底沙摆动特点的“准共振界面波21。

25、由于这四种假说针对细砂级以上粒径而言，均不能解释本文对长江口粗粉砂至极细砂分布区所开展的野外水体底部小尺度紊流和泥沙运动显影观测结果。大量实验室水槽中的水流显影研究(如Ra22，Suer等23)将紊动构造和猝发特征与床面附近的泥沙运动联络起来，获得了一些有意义的结论，即当紊动猝发自床面附近上升时，也夹带了那里的泥沙，假如泥沙的沉速较大，那么在抬升过程中会很快回到河床床面，这样的泥沙就属于跃移质或推移质。假设泥沙沉速较小，那么进入悬浮状态，在悬移质和床面泥沙或推移质之间存在着不断的交换(ffen,1975)24。HeathershaandThrne(1985)25、Thrne等(1989)26和

26、Lapinte(1992)27曾通过安置在河流和海洋环境中近底床含砾砂的流速仪观测到了“象猝发那样的信号。而作者在世界上第一次利用声学悬浮泥沙浓度测定系统和流速仪获得了河口区连续时间序列的可视“紊流猝发(图7)，说明长江口区水深较小，大尺度底形沙波与高流速和含沙量条件下，大、小尺度紊流的猝发特征(时间间隔和历时、频率和强度)及其对粗粉砂和极细砂的泥沙特性(起扬速度和临界运动)直接相关，可能与细砂级大尺度底形沙波与大尺度紊流的关系28相反。致谢：本文曾得到陈吉余院士、王宝灿教授的指导和帮助，还得到李九发教授、何青讲师的帮助，在此一并致谢。参考文献1Kstashuk,R.A.,etal.,1996

27、,Flandsedienttransprtverlargesubaqueusdunes:FraserRiver,anada,Sedientlgy,43,849-863.2VanDenBerg,J.H.,1987.Sedientlgy,34,681-698.3Berne,S.,etal.,1993.J.Sedi.Petrl.,63(5),780-793.4陈吉余主编。上海海岛资源调查。上海文献出版社，P60.1987.5杨世伦。长江口沉积物粒度参数的统计规律及其沉积动力学解释。泥沙研究，1994，(3).6SubitteentheSedientTerinlgy,Aer.Gephys.Unin,“

28、ReprtntheSubitteenSedientTerinlgy,AerianGephysialUnin，Trans.,Aeri.Gephys.Unin,Vl.28,N.6,1947,pp.936-938.7Ashley,G.,etal.,1990,lassifiatinflargesalesubaqueusbedfrs:anelkatanldprble.J.Sedi.Petrl.,Vl.60,N.1,160-172.8VandenBerg,J.H.,etal.,1993.Preditinfsuspendedbedaterialtransprtinflsversiltandveryfines

29、and.aterResuresRes.,Vl.29,N.5,1393-1404.9Kar,P.D.andREiers,.e.,1978.Grainshapeeffetsnsettlingrates.J.Gel.,86,193-209.10王谷谦，乐嘉钻。上海航道科技，1992，总第53期，1-9.11Zhuhengxietal.,1993.nSandaveinYashahannelftheYangtzeEstuary,PreedingsftheSendInternatinalSypsiufRiverSedientatin,aterResuresandEletriPerPress(BEIjing

30、),650-661.12李九发etal.长江口底沙运动规律。海洋与湖沼，Vl.7,N.4,441-450,1995.13hengHeqin,etal.,1998.BedfrsandepisdiresuspensinfarsersiltandveryfinesandinthehangjiangEstuary,hina.PublishedasaSuppleenttES,AerianGephysialUnin,H41E-7.14程和琴等。现代盐水楔河口湾底沙推移速度的估算方法。海洋科学，1998，(1)，27-29.15程和琴等。1996，波、流结合作用下的近岸沙波稳定性研究进展。地球科学进展，Vl.11,N.4.16Shi,Z.,etal.,1997,AustiiagingfhesivesedientresuspensinandreentrainentinthehangjiangEstuary,EasthinaSea,GearineLetters,17,162-168.17Liu,H.K.,ehanisfsedientripplefratin,J.Hyd.Div.,Pr.Aer.S.ivilEngrs.,Vl.83,N.HY2,1957,23.18Nelsn,J.,etal.,1989.ehanisfflverripplesanddunes,J