集合变分法在大气边界层参数反演中应用研究.pdfVIP

集合变分法在大气边界层参数反演中的应用研究i 张晓辉，韩月琪l王云峰1高凡2 (1解放军旺大学气辅洋学院．甫囊211101；292853部队司令部气象台．麓芦岛125106) 摘要 希j用集合交分法(Ensemble．basedVar-哺fioua]Method，以下简称En．Vat)研究了大气边界层EK模式 (考虑了水平气压梯度力、科氏力和湍流粘性力三力平衡的大气边界层模式)中湍流粘性系数的反演问题。 首先利用集合变分法推导了目标涵数关于反演参数的梯度表达式然后绘出了反演计算方塞，最后在簸反 演湍流粘性系数是进行了一系列理想数值试验。数值试验结果表明：零j用集合变分法对EK模式中的湍流 牯性系数的及演比较成功。通过适当调整集合数，切始扰动均方差及初始值的大小，可以进一步提高集合 变分法对端流糨性系数k的反演精度。斑集合变分法应用于大气边界层参数的反演是一种值得研究的计算 方案。 关键词： 集合变分法反演EK模式湍流粘性系数 1引言 由于人类活动主要位子大气边界层中，而且边界层中湍流输运作用对地一气物质能量 循环起着至关熏要的作用，因此各种时空尺度的大气模式和污染物扩散模式都必须要对大 气边界层湍流问题进行足够精确的刻画。丈气湍流问题是一个远来解决的自然科学难题fI_ ：】，所以大气边界层数值模拟和业务数值预报模式以及浩如烟海的研究工作对于湍流摩擦这 一物理过程均采用了经验参数化的方法，即根据火量试验给出人气湍流系数的某种经验形 式。而事实上湍流系数是一个很难确定的物理参数，它与大气的动力热力结构密切相关 【31。因此给出火气湍流系数合理的计算方法无论对于发展理论研究还是改进业务数值天气 预报、污染预报都是非常迫切的。 20世纪80年代在大气科学界提出从反问题的角度，利用最优控制论思想研究数值预 报初值化14l。根据这一提示，可以把确定边界层湍流系数提成～个反问题，即在知道部分 风场观测资料的前提下，根据一定的边界层模式反演湍流系数15l。因此本文利用集合变分 法[61研究了大气边界层EK模式(考虑了水平压梯度力、科氏力和湍流粘性力三力平衡的大 气边界层模式)中湍流粘性系数的反演闯题。 象海洋学院基础理论研究基金． 作者简介：张晓辉，男，1989年生，在读研究生。主要从事大气动力学与数值模拟。 l：zxhuil9891003#163．com 电话F,．mai 论文已在C气象科技》发表 l 考虑到i觜时空密度大气资料获取的困难，本文研究主要是从理论方面构造大气湍流系 数的反演方案。因而也选取大气边界层理论模式作为研究对象。选取的模式是考虑水平气 压梯度力、科氏力和湍流粘性力三力平衡的EK模式，与经典EK模式不同的是湍流系数 随高度变化。 2集合变分梯度算法的推导 集合变分梯度算法的推导流程如图1所汞，_3l。 圈l集合变分梯度算法的流程图 Flowchartoftheensemble-ba．矧：lvariational Fig．1 gradientalgorithm 2 EK模式下湍流系数豇的反演计算方案 2．1EKman方程 定常状态下大气边界层Ekman方程为： 七 +f(v—Vg)=0， a一瑟a一如 ／，●，●●＼／，，，，●一， 七 鼬一赴加一瑟 一f(u—Hg)=0， 边界条件为： (”，V)L=(o，o)，(“，V)|蒯=(“多，V。H)， (2) 2 其中k为湍流系数，“、1，为水平风场纬向、经向分量·“g、％为地转风纬向、经向 分量，f为科氏参数，“多、《为边界层项地转风纬向、经向分量。 2．2离散格式 肚黔水TUi+I--Ui一砖刊。拇r”‰n脚Ⅲ]邝， I丢(七老)i=去(‰监产一t生≯)=吉[t％一(南+‰)K+岛‰]一