澳大利亚奥伯龙湖风致去层理化的数值评估及其对水质管理的启示

《Journal of Water Process Engineering》：Numerical assessment of wind-induced destratification and its implications for water-quality management in Lake Oberon, Australia

　　

强风事件使奥伯龙湖夏季分层减弱最高达1.28 °C。风持续时间解释了分层响应变化的约87%。纵向风通过更大的有效吹程（effective fetch）增强了温跃层位移。Wedderburn数（W）关联了风强迫、分层强度与温跃层响应。基于情景推导的阈值支持风后监测与适应性水库运行。热分层是控制湖泊与水库水动力循环、垂直输运、生物地球化学转化及水质动态的基本过程。暖季太阳辐射建立暖上层（epilimnion）覆盖冷下层（hypolimnion），由变温层（metalimnion）密度梯度抑制垂直交换。下层（hypolimnion）长期隔离限制氧补给，促进底—水界面还原条件，可能增强铁（Fe）、锰（Mn）、无机磷（P）及铵态氮（NH4+–N）释放。在饮用水水库中，这些过程增加氧化、过滤、曝气与化学加药需求，并可能提升蓝藻增殖风险。风作为分层水体主要能量来源，通过界面动量传递产生表面输运、剪切、湍流与流域尺度环流，持续定向强迫引起风涌（wind setup）与温跃层倾斜，风松弛后以内波（internal seiche）振荡重分配动量与势能。风响应取决于风速、持续时间、初始分层强度、有效吹程、盆地形态及风事件时长与内波时间尺度关系。已有观测与数值研究多针对单一事件，定量风—响应关系在许多运行水库仍受限，澳大利亚饮用水水库尤甚。研究人员以新南威尔士州温带饮用水水库奥伯龙湖为对象，其细长且形态异质盆地导致有效吹程方向差异显著。研究采用校准并独立验证的三维水生态模型（AEM3D）系统量化风速、持续时间与方向对热响应的调节，整合常规热指标、Schmidt稳定性（St）、Wedderburn数（W）与基于情景的敏感性分析，将水动力机制转化为工程与运行决策指标。目标包括：（i）校准验证AEM3D；（ii）量化风速与时长对温跃层位移、表底温差（ΔT）与水柱稳定性的影响；（iii）确定风向、有效吹程与盆地形态对空间异质混合的调控；（iv）用W与内波时间尺度解释响应；（v）将响应范围转为风后监测与适应性管理的情景阈值。结果表明强风使夏季分层减弱达1.28 °C；风持续时间解释约87%分层响应变异；纵向风借有效吹程增强温跃层位移；W联结风强迫、分层与温跃层响应；情景阈值支持风后监测与适应性运行。

研究背景方面，热分层是湖泊与水库中控制水动力循环、垂直输运、生物地球化学转化及水质动态的基本过程。暖季太阳辐射形成暖上层覆盖冷下层，变温层密度梯度抑制垂直交换，下层长期隔离会限制氧补给并促进铁、锰、无机磷与铵态氮等还原释放，在饮用水水库中增加处理难度与蓝藻风险。风是分层水体主要能量来源，通过界面动量传递产生表面输运、剪切与流域尺度环流，持续定向强迫引起风涌与温跃层倾斜，风松弛后以内波振荡重分配能量。然而风响应受风速、持续时间、初始分层、有效吹程、盆地形态共同控制，既有研究多针对单一事件，许多运行水库缺乏定量化风—响应关系，澳大利亚饮用水水库尤为不足。奥伯龙湖位于新南威尔士州中央高地，1946年筑坝形成，库容约45000 ML，面积410 ha，最大深度近30 m，Fish River贡献约87%年入流，细长且形态异质导致有效吹程方向差异显著，是检验风强迫—温跃层位移—内波响应—分层减弱相互作用的合适系统。研究人员开展本研究旨在用校准并独立验证的三维水生态模型系统量化风速、持续时间与方向对热响应的调节，并将水动力机制转为运行管理指标。

关键技术方法方面，研究人员以2023年8月至2025年7月为模拟期，前接30天自旋_up，初始温度场来自观测剖面，流速场设静水。水平网格50 m覆盖约4.1×4.1 km流域，垂向在表层与变温层加密。气象强迫含气温、风速风向、湿度、气压、辐射与降雨，边界含河流入流（流量与温度）、出流与无法向流侧边界，底床摩擦与绝热边界。模型为AEM3D，采用Boussinesq与静力近似，垂速由连续方程诊断。校准期2023年8月至2024年2月调垂直混合与表面热交换参数，验证期2024年3月至2025年7月用未参与率定观测，以RMSE、Nash–Sutcliffe效率（NSE）、偏差与Pearson相关系数评价。风数值实验在夏季峰值分层后施加恒定风：18 km h^?1下分1、3、5天与南、北、东、西风，36 km h^?1仅南、北风1天匹配对照。指标用表底温差（ΔT）、混合效率（ME）、Schmidt稳定性（S_t）、去层化效率（DE）与Wedderburn数（W= g′h_e²/(u^*2L)，g′为折减重力，h_e为上层深度，u^*为水侧摩擦速度，L为有效吹程）。

研究结果方面，3.1节基线分层显示暖季表温18–20 °C、底温低约5 °C即ΔT≈5 °C，温跃层由12月初约5 m加深至2月末10 m，下游深盆热稳定性最高、下层隔离最强，冬季混合由冷却削弱浮力梯度实现，模型对季节演变与垂向梯度复现良好但内波与垂速仅属模型派生机制。3.2节风速响应中，18 km h^?1单日南、北风引起可测温跃层位移与部分向下传热，东、西风较弱；36 km h^?1匹配单日南、北风相对18 km h^?1使ΔT降低与ME增32.2%，S_t降低与DE增30.4%，但短时高强度仍不足以全混，说明风速须与时长方向联评。3.3节风持续时间表明其解释热响应指标变异约86.6–86.7%，远高于风向10.6–10.8%，1天强迫仅扰上层，3–5天持续使下层增温、ΔT与S_t进一步降，内部seiche与斜坡耗散促局部混合；南、北向（纵向）因有效吹程大位移最强，东、西向弱。W值低表示风强迫相对浮力强、温跃层倾斜潜力大（W<1非全混），与内波时间尺度匹配时响应放大。研究人员由此给出情景阈值：中等偏强18 km h^?1需≥3天纵向风方达中等混合效率（10–20%），5天可入强（20–30%）；36 km h^?1单日纵向风接近强混合但未均质；混合效率分级为弱<10%、中等10–20%、强20–30%、很强>30%。

讨论与结论部分，研究人员指出定量风—响应关系具强场地特异性，奥伯龙湖形态使纵向风最有效，风持续时间而非瞬时风速主导去层化变异，W提供无量纲桥接风强迫、上层深度与有效吹程的解释框架，内波周期匹配增强响应。因缺溶解氧连续观测，水质含义经分层—下层缺氧—氧化还原—内源释养链条间接推断，不作直接DO验证。研究意义在于把已识机制转为水库专属运行指标：风后监测表底温差、S_t与W，可预警下层缺氧与营养盐释放，指导饮用水水库在强风事件后调整取水层、加药与曝气，支持适应性运行。论文发表于《Journal of Water Process Engineering》。结论可译为：强风事件通过风持续时间主导的机制使奥伯龙湖夏季热分层减弱最高1.28 °C；风持续时间解释约87%分层响应变异，风速由18增至36 km h^?1使匹配单日响应增约30–32%；纵向风借更大有效吹程增强温跃层位移；Wedderburn数联结风强迫、分层强度与温跃层响应；基于情景阈值支持风后监测与适应性饮用水水库运行。