涌潮河口底摩擦系数伴随资料同化研究 Study of bottom friction coefficient in tidal bore affected estuaries using the adjoin method of data assimilation.pdfVIP

第22卷第3期 水科学进展 V01．22，No．3 2011年5月 ADVANCESINWATERSCIENCE May，2011 DOI：CNKI：32—1309／P1449．015 涌潮河口底摩擦系数伴随资料同化研究 陈耀登1’2，闵锦忠1’2，高玉芳1’3 (1．南京信息工程大学气象灾害省部共建教育部重点实验室，江苏南京210044； 2．南京信息工程大学大气科学学院，江苏南京210044；3．淮河流域气象中心，安徽蚌埠233000) 摘要：为更合理选取涌潮河口的底摩擦系数，基于无结构三角网格和有限体积法，建立了适用于反演潮汐涌潮河口 底摩擦系数随空间和涨、落潮变化的伴随资料同化模型。一系列“孪生”实验均较精确地反演出了给定底摩擦系 数。孪生实验表明底摩擦系数的空间分布会影响反演的精度，也表明充足的观测数据资料能够提高反演的准确性。 149和0．001520(相当于曼 利用实测数据进行了钱塘江河口伴随同化实验，得到涨、落潮底摩擦系数分别为0．000 宁系数为0．012206和0．038987)，利用该值模拟出了潮波变形和涌潮形成等现象，验证了前人在钱塘江河口数值 模拟时关于底摩擦系数选取经验的合理性，反映了钱塘江河口底摩擦“涨小落大”的实质。 关键词：伴随方法；资料同化；底摩擦系数；涌潮；河口 中图分类号：TVl31．2文献标志码：A 文章编号：1001-6791(2011)03-0400-07 在水动力数值模拟计算中，底部切应力项中的不确定性参数——底摩擦系数作为一个敏感参数，其数值 的合理计算和选取对数值模拟的结果有着直接影响，如何方便合理地确定这个参数，一直是数值计算研究人 员重视的问题…。传统上采用人工试算的方法来率定此参数，工作量往往十分巨大。基于最优控制和变分 理论的伴随数据同化方法，能够利用伴随方程方便地计算出控制参数下降的梯度，并将控制参数的确定变成 依赖于区域内观测值的数值迭代过程，成为参数估计的有力工具被应用于参数优化选取的研究中嵋…。前人 的研究多基于传统有限差分法，研究的区域为大范围海域或是比较顺直的河道。事实上在对河口等区域进行 数值模拟时，物理变量可能在较小区域内变化较大，甚至出现涌潮等间断现象，且此类岸线曲折复杂，在此 类区域数值模拟时对所用的数值模型的计算方法和网格分辨率要求较高¨…。而且研究还表明¨。：当潮波自 深海向近岸传播时，由于地形变浅，潮流的加速作用使运动的阻力增大，而退潮时则明显减小。在潮汐河口 的数值模拟研究，常通过判别潮流的运动状态分别给定不同的底摩擦系数呻…。因此，在伴随同化过程中需 要判别潮流的运动状态，以更合理反演该系数。 基于此，本文通过建立适用于潮汐涌潮河口底摩擦系数随空间和涨落潮变化的底摩擦校正关系式，利用 具有较高空间分辨率且适应间断解问题的无结构三角网格有限体积法的正模型和伴随数据同化模型，进行了 概化河口和钱塘江河口的底摩擦系数反演的伴随同化实验研究。 1 基本模型 1．1 正模型 深度平均的非线性守恒型浅水方程为 收稿Et期：2010-06-10：网络出版时间：2011-05-14 网络出版地址：http：／／www．cnki．net／kcms／detail／32．1309．P1449．015．html 基金项目：国家自然科学基金资助项目；淮河流域气象中心开放基金资助项目(HRM200903) 作者简介：陈耀登(1980一)，男，福建沙县人，博士，副教授，主要从事数值模拟与资料同化研究。 E—mail：keyu@nuist．edu．cn 万方数据 第3期 陈耀登，等：涌潮河口底摩擦系数伴随资料同化研究 40l —OH+型+型：o Ot Ox Oy