2025年南信工数值天气预报试题C及答案

2025年南信工数值天气预报试题C及答案一、名词解释（每题4分，共20分）1.正压模式：假设大气是正压的，即密度仅为气压的函数（ρ=ρ(p)），水平方向上气块的温度梯度为零，垂直速度由水平辐散辐合决定。该模式忽略斜压性，适用于大尺度天气系统的简化模拟，常用于研究水平风场与气压场的平衡关系。2.半隐式格式：时间积分方案中对慢变项（如平流项）采用显式计算，对快变项（如重力波、声波）采用隐式处理的混合格式。其通过隐式部分抑制快波的高频振荡，扩大稳定条件限制，同时保持显式部分的计算效率，是中尺度和全球模式中常用的时间离散方法。3.资料同化：将观测资料与模式预报结果结合，通过最优估计方法提供初始场的过程。核心是利用统计方法（如三维变分、四维变分、集合卡尔曼滤波）最小化观测误差与模式背景误差的加权距离，提高初始场的准确性，从而改善预报效果。4.B-grid：一种水平网格配置方案，其中水平速度（u、v）定义在网格元的边中点，位势高度、温度、气压等标量定义在网格元中心，垂直速度（ω或w）定义在网格元顶点。该网格通过交错布局减少计算噪声，平衡离散误差，广泛应用于WRF、ARPS等模式。5.地形追随坐标：垂直坐标设计中，坐标面随地形起伏变化（如σ坐标：σ=(p-p_top)/(p_s-p_top)，p为气压，p_s为地面气压，p_top为模式顶气压）。其优势是将复杂地形边界转化为坐标面的水平边界，简化数值计算；缺点是地形陡峭时可能引发坐标面倾斜，导致计算误差增大（如“地形波动”）。二、简答题（每题8分，共40分）1.数值模式中时间积分方案的选择需考虑哪些关键因素？答：时间积分方案的选择需综合以下因素：（1）稳定性：快波（如重力波、声波）的传播速度决定CFL条件，隐式或半隐式格式可放宽稳定限制；（2）精度：高阶格式（如Runge-Kutta）能减少时间截断误差，但计算成本高；（3）计算效率：显式格式计算量小但步长受限，隐式格式需解线性方程组（如三对角矩阵），需权衡时间步长与迭代次数；（4）守恒性：能量、质量等物理量的守恒对长期积分至关重要，如平流项的离散需保持正定性；（5）模式动力框架特性：全球模式因科里奥利参数随纬度变化（β平面近似），需考虑时间方案与空间离散的匹配性。2.为何数值模式中需要引入次网格尺度参数化？其核心假设是什么？答：次网格尺度参数化的必要性源于模式分辨率限制：实际大气中存在比网格尺度（Δx）更小的湍流、积云、边界层等过程，无法直接解析。核心假设是“次网格尺度过程对网格平均量的影响可通过网格平均量的函数表示”（如用网格平均温度、湿度计算积云对流的垂直输送）。参数化需满足：（1）物理合理性：反映次网格过程的主要机制（如积云的热量、水汽垂直输送）；（2）计算可行性：参数化方案需基于可解析的网格变量（如垂直速度、温度梯度）；（3）普适性：在不同天气背景（如对流性降水、层状云）下均能合理描述次网格效应。3.水平分辨率提高对数值模式模拟的影响主要体现在哪些方面？需注意哪些潜在问题？答：影响：（1）更好解析中小尺度系统（如雷暴、中尺度涡旋），改善降水、风场等要素的精细化预报；（2）减少次网格参数化的依赖（如积云参数化可在高分辨率下关闭，直接模拟深对流）；（3）更准确刻画地形、下垫面（如城市、海岸线）的动力与热力作用。潜在问题：（1）计算成本指数增长（如水平格点数约为(1/Δx)²，总计算量约为(1/Δx)^3）；（2）高分辨率下快波（如声波）传播速度加快，需减小时间步长（Δt∝Δx），进一步增加计算负担；（3）次网格参数化方案可能失效（如边界层参数化在Δx<1km时需调整，因湍流尺度与网格尺度接近）；（4）观测资料分辨率不足可能导致初始场误差增大（“过分辨率”问题）。4.初始场误差对数值预报的影响机制是什么？可预报性理论如何描述这一过程？答：初始场误差通过动力不稳定机制（如正压/斜压不稳定、非线性相互作用）随时间增长。具体机制：（1）线性阶段：误差增长率由大气中的不稳定模态（如斜压波）决定，增长率与扰动波长、垂直切变等有关；（2）非线性阶段：误差通过能量级串向不同尺度传播（如小尺度误差激发大尺度波动，或大尺度误差分裂为小尺度扰动），最终导致预报失效。可预报性理论（如Lorenz提出的“蝴蝶效应”）用误差增长的e折叠时间（τ）描述：若初始误差为ε₀，预报时间t后误差为ε(t)=ε₀·e^(t/τ)，当ε(t)接近气候噪声水平时，预报失去意义。全球模式中，中期预报（10天左右）的可预报性主要受斜压波不稳定控制，τ约为2-3天；而热带地区因斜压性弱，可预报性时间更长（约2周）。5.有限差分法与谱方法在数值模式动力框架中的优缺点比较。答：有限差分法（FDM）：优点：（1）计算局部性强（仅需相邻网格点），适合并行计算；（2）格式灵活，可构造高阶精度（如WENO格式）或自适应网格；（3）对复杂地形（如σ坐标）的离散更易实现。缺点：（1）截断误差随分辨率提高而减小，但高阶格式可能引发数值振荡（需加入人工耗散）；（2）全局守恒性（如角动量守恒）难以严格满足，需通过格式设计补偿。谱方法（SM）：优点：（1）全局高精度（谱精度，误差随模态数指数衰减）；（2）严格保持物理量的守恒性（如能量、角动量），适合长期气候模拟；（3）自然满足周期性边界条件（如全球模式的经向周期）。缺点：（1）计算全局性强（需计算所有模态的相互作用），并行效率低；（2）对复杂地形或非均匀网格（如区域模式）的适应性差；（3）非线性项（如平流项）的计算需通过快速傅里叶变换（FFT）转换到物理空间，计算成本高。三、计算题（每题10分，共40分）1.考虑一维平流方程∂u/∂t+c∂u/∂x=0（c>0为常数），采用向前时间、向后空间的差分格式（FTBS）：u_j^(n+1)=u_j^ncΔt/Δx(u_j^nu_{j-1}^n)（1）判断该格式是否守恒质量（即∑u_j^(n+1)Δx=∑u_j^nΔx）；（2）用VonNeumann稳定性分析推导其稳定条件；（3）计算截断误差的阶数。答：（1）质量守恒性：对空间求和∑u_j^(n+1)Δx=∑[u_j^ncΔt/Δx(u_j^nu_{j-1}^n)]Δx=∑u_j^nΔxcΔt∑(u_j^nu_{j-1}^n)。由于求和项为望远镜级数，∑(u_j^nu_{j-1}^n)=u_N^nu_0^n（假设网格点j=0到N），若边界条件为周期性（u_N=u_0），则第二项为0，质量守恒；否则不严格守恒。通常假设周期性边界，故该格式守恒质量。（2）稳定性分析：设扰动为u_j^n=A^ne^(ikjΔx)（k为波数），代入格式得：A=1cΔt/Δx(1e^(-ikΔx))=1α(1e^(-iθ))（α=cΔt/Δx，θ=kΔx）振幅因子|A|²=[1α(1cosθ)]²+[αsinθ]^2=12α(1cosθ)+α²(12cosθ+cos²θ+sin²θ)=12α(1cosθ)+2α²(1cosθ)稳定条件为|A|²≤1，即-2α(1cosθ)+2α²(1cosθ)≤0→2(1cosθ)(α²α)≤0。因1cosθ≥0，故需α²α≤0→α(α1)≤0。因α>0（c>0），故稳定条件为0<α≤1，即cΔt/Δx≤1（CFL条件）。（3）截断误差：将u(x,t)在(x_j,t_n)处泰勒展开：u_j^(n+1)=u+Δtu_t+(Δt²/2)u_tt+O(Δt³)u_j^ncΔt/Δx(u_j^nu_{j-1}^n)=ucΔt/Δx[u(uΔxu_x+(Δx²/2)u_xxO(Δx³))]=ucΔt/Δx(Δxu_x(Δx²/2)u_xx)=ucΔtu_x+(cΔtΔx/2)u_xx原方程u_t=-cu_x，故u_tt=-cu_xt=c²u_xx（对t求导）。代入左边得：左边=u+Δt(-cu_x)+(Δt²/2)(c²u_xx)+O(Δt³)右边=ucΔtu_x+(cΔtΔx/2)u_xx截断误差T=左边-右边=(c²Δt²/2cΔxΔt/2)u_xx+O(Δt³,ΔxΔt²)=(cΔt/2)(cΔtΔx)u_xx+O(Δt³,ΔxΔt²)。当Δt∝Δx（CFL条件下α=O(1)），则T=O(Δt,Δx)，即一阶截断误差。2.推导浅水方程在σ坐标（σ=z/H，H为静止流体深度）下的转换过程，写出转换后的水平动量方程（忽略科里奥利力）。答：浅水方程原始形式（z坐标）：∂u/∂t+u∂u/∂x+v∂u/∂y+w∂u/∂z=-g∂η/∂x∂v/∂t+u∂v/∂x+v∂v/∂y+w∂u/∂z=-g∂η/∂y∂η/∂t+∂(uH)/∂x+∂(vH)/∂y=0（η为自由面高度，H为静止深度，总深度h=H+η）σ坐标定义：σ=z/H，故z=Hσ，dz=Hdσ。垂直速度w=dz/dt=Hdσ/dt+σdH/dt（因H可能随时间变化？不，浅水模式中H为常数，故w=Hdσ/dt→dσ/dt=w/H。变量转换：∂/∂z=(1/H)∂/∂σ，∂/∂t=∂/∂t+w∂/∂z=∂/∂t+(w/H)∂/∂σ（随体导数）。水平动量方程中的垂直输送项w∂u/∂z=(w/H)∂u/∂σ=dσ/dt∂u/∂σ。因此，x方向动量方程转换为：∂u/∂t+u∂u/∂x+v∂u/∂y+dσ/dt∂u/∂σ=-g∂η/∂x由于浅水模式中u、v与z无关（正压假设），∂u/∂σ=0，垂直输送项消失，方程简化为：∂u/∂t+u∂u/∂x+v∂u/∂y=-g∂η/∂x（与原z坐标方程形式一致）。注：若考虑非正压（u随z变化），则需保留垂直输送项，但浅水模式通常假设水平速度不随高度变化，故σ坐标转换后方程形式不变。3.二维线性波动方程∂²φ/∂t²=c²(∂²φ/∂x²+∂²φ/∂y²)，采用中心差分格式：(φ_j,k^(n+1)2φ_j,k^n+φ_j,k^(n-1))/Δt²=c²[(φ_{j+1,k}^n2φ_j,k^n+φ_{j-1,k}^n)/Δx²+(φ_{j,k+1}^n2φ_j,k^n+φ_{j,k-1}^n)/Δy²]设Δx=Δy=h，求CFL数的稳定条件。答：VonNeumann分析：设φ_j,k^n=A^ne^(i(kxjh+kykh))（kx、ky为x、y方向波数），代入格式得：(A²2A+1)/Δt²=c²[(-4sin²(kxh/2)/h²)+(-4sin²(kyh/2)/h²)]整理得：A²2A+1=-4c²Δt²/h²[sin²(kxh/2)+sin²(kyh/2)]令α=cΔt/h（CFL数），θx=kxh，θy=kyh，则：A²2A+1=-4α²[sin²(θx/2)+sin²(θy/2)]方程化为A²2A+[1+4α²(sin²(θx/2)+sin²(θy/2))]=0特征根A=[2±√(44(1+4α²(sin²(θx/2)+sin²(θy/2))))]/2=1±i2α√(sin²(θx/2)+sin²(θy/2))振幅因子|A|=√[1²+(2α√(sin²(θx/2)+sin²(θy/2)))²]=√[1+4α²(sin²(θx/2)+sin²(θy/2))]≥1，这意味着格式无条件不稳定？但实际二维波动方程的中心差分格式稳定条件需满足CFL数α≤1/√2（当Δx=Δy时）。修正分析：正确的二维波动方程中心差分格式稳定条件推导应为：波动方程可分解为两个一维平流方程的组合（∂/∂t±c∇）φ=0，因此稳定条件需同时满足两个方向的CFL条件。对于二维情况，最大波数对应θx=θy=π（最短波长2h），此时sin²(π/2)=1，故|A|²=1+4α²(1+1)=1+8α²。为使|A|≤1，需8α²≤0，矛盾，说明中心差分格式对二维波动方程显式积分不稳定，需采用隐式格式或限制α≤1/√2（通过能量法推导）。正确稳定条件为cΔt/h≤1/√2。4.半隐式格式处理快波的核心思想是什么？以线性化的浅水方程为例，推导其时间离散方案（假设时间步长为Δt，显式处理平流项，隐式处理重力波项）。答：核心思想：将方程中的快变项（如重力波的气压梯度项）隐式离散，慢变项（如平流项）显式离散，通过隐式部分消除快波的高频振荡，放宽CFL条件限制。线性浅水方程（忽略平流项）：∂u/∂t=-g∂η/∂x∂η/∂t=-H∂u/∂x半隐式格式中，平流项（此处为0）显式，重力波项隐式。时间离散采用中心差分（n为当前时间层，n+1为未来层，n-1为过去层）：(u_j^(n+1)u_j^(n-1))/(2Δt)=-g(η_{j+1}^nη_{j-1}^n)/(2Δx)（显式处理？不，半隐式应隐式处理重力波项，即η用n+1层）正确离散：将重力波项（-g∂η/∂x）在n+1层计算（隐式），平流项（此处无）在n层计算（显式）。采用向前-向后格式：u_j^(n+1)=u_j^nΔtg(η_{j+1}^(n+1)η_{j-1}^(n+1))/(2Δx)（隐式，u的更新依赖η的n+1层）η_j^(n+1)=η_j^nΔtH(u_{j+1}^nu_{j-1}^n)/(2Δx)（显式，η的更新依赖u的n层）将η的显式表达式代入u的隐式方程，消去η^(n+1)，得到关于u^(n+1)的方程组，通过迭代或矩阵求解。该格式通过隐式处理重力波项，将稳定条件从显式的Δt≤Δx/√(gH)放宽为Δt≤O(Δx)（与水平平流项的CFL条件一致），从而允许更大的时间步长。四、综合分析题（每题10分，共30分）1.比较WRF模式与ECMWF全球模式在动力框架设计上的主要差异，并分析其对预报性能的影响。答：WRF（WeatherResearchandForecasting）为区域模式，ECMWF（欧洲中期天气预报中心）模式为全球模式，动力框架差异主要体现在：（1）坐标系统：WRF常用地形追随坐标（如η坐标）适应复杂地形；ECMWF采用混合坐标（Hybridσ-pressure坐标），高层为气压坐标（避免地形陡峭时的坐标倾斜），低层为σ坐标以解析边界层。（2）水平网格：WRF多采用ArakawaC-grid（u、v定义在网格边中点，标量在中心）减少计算噪声；ECMWF早期用谱方法（T模式），现逐步转向有限体积法（如IFS模式的FV3动力核），兼顾全局守恒与区域高分辨率嵌套。（3）时间积分：WRF对浅对流采用显式方案（Δt=10-30秒），对深对流关闭积云参数化（Δx<3km）；ECMWF因全球范围需大时间步长（Δt=15-30分钟），采用半隐式-半拉格朗日（SISL）格式，隐式处理重力波，拉格朗日方法减少平流误差。（4）守恒性：ECMWF谱模式严格守恒角动量、能量，适合长期预报；WRF有限差分格式通过通量修正（如MUSCL方案）近似守恒，更适应中小尺度强梯度（如锋面）。对预报性能的影响：WRF在区域高分辨率（Δx=1-10km）下能更好解析地形强迫（如地形降水）和中小尺度系统（如雷暴），但受限于侧边界条件误差；ECMWF全球模式通过谱方法或有限体积法的全局一致性，提供更准确的大尺度环流背景（如西风带、副热带高压），其半拉格朗日方法减少了平流过程的虚假扩散，有利于中期（10天）预报的稳定性。两者结合（如WRF嵌套于ECMWF分析场）可兼顾大尺度背景与区域细节。2.初始扰动的增长率与数值天气预报的可预报性有何联系？结合Lorenz63模型与实际大气，说明可预报性上限的主要影响因素。答：初始扰动增长率决定了误差从不可观测尺度（~10⁻⁶K）增长到可影响预报的尺度（~1K）的时间，即可预报性上限。Lorenz63模型（简化的对流系统）中，非线性项导致扰动指数增长（λ≈0.9天⁻¹，e折叠时间τ≈0.7天），这解释了短期预报（1-3天）的准确性和中期预报（10天）的不确定性。实际大气中，扰动增长率受