湖泊科学   2022, Vol. 34 Issue (3): 994-1005.  DOI: 10.18307/2022.0323
0

研究论文——流域水文与水资源安全

引用本文 [复制中英文]

高昂, 吴修锋, 吴时强, 戴江玉, 王芳芳, 吴晨晖, 有限水域近水面气流演化特征试验研究. 湖泊科学, 2022, 34(3): 994-1005. DOI: 10.18307/2022.0323
[复制中文]
Gao Ang, Wu Xiufeng, Wu Shiqiang, Dai Jiangyu, Wang Fangfang, Wu Chenhui. Experimental study on the evolution characteristics of near-surface airflow in limited waters. Journal of Lake Sciences, 2022, 34(3): 994-1005. DOI: 10.18307/2022.0323
[复制英文]

基金项目

国家重点研发计划项目(2018YFC0407200)、国家自然科学基金项目(51909169)、江苏省基础研究计划(自然科学基金)项目(SBK2019042181)和中央级公益性科研院所基本科研业务费专项资金项目(Y121006, Y120010)联合资助

通信作者

高昂, E-mail: agao@nhri.cn

文章历史

2021-06-28 收稿
2021-08-20 收修改稿

码上扫一扫

有限水域近水面气流演化特征试验研究
高昂1,2 , 吴修锋1,2 , 吴时强1,2 , 戴江玉1,2 , 王芳芳1,2 , 吴晨晖1,2     
(1: 南京水利科学研究院, 水文水资源与水利工程科学国家重点实验室, 南京 210029)
(2: 水利部太湖流域水治理重点实验室, 南京 210029)
摘要:气流作为湖泊、湿地等有限水域的主要驱动力之一, 其演化特征是研究水-气间能量与物质传递的基础, 决定了水域水环境与水生态格局, 具有重要研究意义. 采用室内风洞水槽, 通过设置不同试验风速, 研究近水面气流与风速和吹程的响应关系与演化特征, 并将光滑壁面条件下的气流特征作为试验对照组, 分析时间域内风速垂向分布形式、摩阻风速、特征粗糙度等变化规律, 剖析频率域内风速的主频变化与波致气流边界层厚度等参数的分布特征. 结果表明: 近水面时均风速符合对数分布, 摩阻风速与特征粗糙度均与风速呈正相关, 且同等风速与吹程条件下两者均大于光滑壁面条件, 波致气流边界层厚度与风速和吹程均呈正相关. 分析了采用风速作为气流动力参数存在不妥的原因, 认为采用摩阻风速或切应力作为气流动力参数更具合理性, 基于前人研究数据的统计分析, 给出了特征粗糙度与风速、吹程和水深的回归关系式. 探讨了气流-风浪“同频率”的互馈模式, 发现波致气流边界层厚度不只与波长有关, 还受波高、波周期和水深等多因素影响, 建立了考虑风浪非线性特征的波致气流边界层厚度关系式.
关键词有限水域    演化特征    特征粗糙度    气流边界层    互馈模式    
Experimental study on the evolution characteristics of near-surface airflow in limited waters
Gao Ang1,2 , Wu Xiufeng1,2 , Wu Shiqiang1,2 , Dai Jiangyu1,2 , Wang Fangfang1,2 , Wu Chenhui1,2     
(1: State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Nanjing Hydraulic Research Institute, Nanjing 210029, P.R.China)
(2: Key Laboratory of Taihu Basin Water Management, Ministry of Water Resources, Nanjing 210029, P.R.China)
Abstract: As one of the main driving forces of limited water areas such as lakes and wetlands, the evolution characteristics of airflow are the basis for the study of water-air energy and material transfer, which determines the water environment as well as water ecological pattern of water area, and has important research significance. Using wind tunnel flume, by setting different wind speeds, the response relationship and evolution law of near-surface airflow with wind speed and blowing distance were studied, particularly, flow characteristics under the condition of smooth wall were taken as the test control group, vertical distributions of wind speed, friction wind speed, roughness and so on were analyzed in time domain, and domain frequency change of wind speed and boundary layer thickness of wave-induced airflow are analyzed in the frequency. Results show that the average wind speed near the water surface conforms to the logarithmic distribution, the friction wind speed and characteristic roughness are all positively correlated with the wind speed, and both of them are greater than the smooth wall condition under the same wind speed and blowing distance, and the wave-induced airflow boundary layer thickness is positively correlated with the wind speed and blowing distance. The reasons for the improper use of wind speed as the aerodynamic dynamic parameter are analyzed, and it is considered that it is more reasonable to use friction wind speed or shear stress as the aerodynamic dynamic parameter, and the regression relationship between characteristic roughness and wind speed, blowing distance and water depth is given. The mutual feed model of airflow and wind wave with "same frequency" is discussed, it is found that the thickness of the wave-induced airflow boundary layer is not only related to the wavelength, but also affected by wave height, wave period and water depth, and the relation formula of the wave-induced airflow boundary thickness considering the nonlinear characteristics of wind and wave is established.
Keywords: Limited waters    evolution characteristics    characteristic roughness    airflow boundary layer    mutual feed model    

全球71个典型湖泊研究结果表明,自1980s以来全球湖泊蓝藻“暴发”整体呈加剧趋势,水环境污染及水生态退化是湖泊面临的普遍问题[1]. 秦伯强[2]在回顾太湖治理历程的基础上,发现目前湖沼学研究与湖泊物理环境和过程的结合度不够,指出未来应深化空气动力学、水动力学与生物地球化学等多学科交叉研究. 对于大多数湖泊而言,气流是湖泊水体的主要驱动力之一,水气相互作用对大气与湖泊之间能量与物质传递至关重要[3]. 比如,当风速大于某一临界值时,湖泊水面的气体(如CO2、CH4等)传输速率与风速呈显著正相关[4-5],然而,不同研究得到的风速-气体传输速率之间的关系存在较大差异,造成这种差异的主要原因之一是不同研究使用的风速测点高度不同(同一测线不同高度处风速不同). 为了减小这种差异,部分研究采用风速廓线公式外延计算得到水面上方10 m高度处风速u10,并建立u10与气体传输速率的关系. 然而,对于面积有限的湖泊而言,摩阻风速u*和特征粗糙度z0与风速和吹程有关[3],导致使用风速廓线公式过程中选取u*z0时存在经验性,仍会造成风速-气体传输速率关系的不确定性. 可见目前对湖泊水域气流动力参数选取方面还未达成一致,需要通过近水面气流演化特征的研究进一步明确.

湖泊水域面积有限,属于有限水域,水体在气流驱动下产生风浪和风生流[6],两者随风速和吹程增加逐渐演化,如风浪尺度与动力强度的增加[7],风生流强度与风生环流模式的改变[8]等. 相应地,气流在与水流耦合作用过程中也逐渐演变,如风速随吹程的增加[9],气流边界层的逐渐发展[10-11]等. 近水面气流演化特征研究属于边界层流动问题,一直是海洋、湖泊领域关注的热点[12-13],但由于风浪的非恒定特征,近水面气流特征比平板气流更加复杂,其研究更具挑战性. 湖面风场的不均、不稳特点使得现场观测多具短历时、点测量特征,获取的数据系统性较差,制约了对湖面气流演化特征的深入研究. 风洞水槽因具有精细可控的优势,一直是水气作用规律研究的重要手段[14]. Liberzon和Shemer[9]、Longo[10]、Shaikh和Siddiqui[15]等通过试验研究了水气相互作用规律,阐述了水-气相互作用与气-平板相互作用之间的异同,发现风浪的存在增加了气流下垫面的粗糙度,剖析了不同风浪尺度对气流特征的差异影响,发现近水面一定高度内时均风速随着距离水面降低而快速减小[15],即近水面气流存在明显的风浪扰动区(波致气流边界层),该区域内水气作用强烈,是能量与物质交换最强烈的区域,也是水气作用研究中重点关注的区域,其中,波致气流边界层厚度δwind的确定是主要研究内容之一. 然而,但不同研究对δwind的判定存在差异,Sullivan等[16]试验研究发现小风速时δwind为主频波长Ld的0.16倍,而Liberzon和Shemer[9]的试验结果表明δwindLd的0.32倍,Zavadsky和Shemer[17]发现大风速时δwind可增大至Ld的0.48倍. 可见,目前多认为波致气流边界层厚度是主频波长的单一函数,但不同研究结果却未达成一致. 这就会给气象与环境监测带来难题,可能导致在湖面架设气象、环境等监测仪器时存在架设高度的不确定性,如测点架设在边界层以上,测量的数据可能不具代表性,且具有误导性. 因此,需理清水面波致气流边界层厚度与风浪尺度的关系,以指导湖泊、湿地等有限水域气象与环境监测.

有限水域气流动力参数选取不明与波致气流边界层厚度不清的现状彰显了目前湖泊水环境与湖泊动力的交叉研究还存在短板,其中,湖面近水面气流作为湖泊主要动力条件,研究近水面气流演化特征是补充上述短板的关键前提. 基于此,本研究以湖泊、湿地等有限水域为研究对象,在风洞水槽中模拟研究有限水域近水面气流演化特征. 首先,从时间域和频率域两个维度对风速垂向分布、摩阻风速、特征粗糙度、气流主频及边界层厚度等基本特征进行剖析. 在此基础上,量化分析气体传输率关系建立中存在不确定性的原因,提出更合理的水面气流动力参数,为水面气体传输速率参数选取提供理论依据; 探讨气流与风浪的互馈模式,明晰波致气流边界层厚度与风浪间的关系,为有限水域原位监测中确定仪器布设高度提供科学指导.

1 试验设计 1.1 风洞水槽介绍

试验在南京水利科学研究院风洞水槽中开展(图 1),该水槽中已开展了多项有限水域风、浪、流的理论与应用基础研究[18-20]. 水槽有效段长22.5 m,高1.2 m,宽1 m,顶部为弧形盖板(高0.4 m),进口及出口设置渐变段,边壁为钢化玻璃,底坡为平坡. 风洞进口布置了蜂窝管,以增加气流的紊动程度与均匀性,气流在经过蜂窝管整流后基本为均匀流. 考虑到天然湖泊近岸存在的人工消浪设施或水生植物的消浪功能,水槽尾部设置消浪板. 吸风机设置在风洞尾部,通过橡胶减震片与风洞相连,避免风机震动干扰水体运动. 风洞内可产生最大风速约20 m/s.

图 1 风洞水槽及仪器布设 (F4.5~F19.5表征吹程F为4.5~19.5 m) Fig.1 Wind tunnel flume and instrument layout
1.2 数据测量与数据处理

风洞内沿程布置6个测量断面,断面间距3 m. 在每个测量断面中线布设一条风速测线,每个测线沿垂向布置16支热线风速传感器,传感器间距3 cm,采集频率10 Hz,精度1 cm/s,可同步采集水面上方约50 cm高度的瞬时风速. 对应地,在每个测线下风向10 cm处沿风向平行布置2根(1对)电容式波高仪,共计12根,波高仪间距5 cm,采集频率75 Hz,精度1 mm.

风浪主频波长Ld基于1对波高仪同步测量的水面波动,采用互相关分析方法得到[20];气流与水面波动的频谱图通过快速傅里叶变换(FFT)方法得到,主频为频谱图中能量峰值对应的频率.

1.3 试验方案

本研究为针对有限水域气流演化特征的应用基础研究,试验方案设计在考虑天然水域实际条件的同时,结合研究内容综合考虑. 首先进行无水时(玻璃底壁)气流特征试验研究,而后开展水面气流特征研究. 考虑到湖泊的浅水特征,试验水深d为0.15 m,水温为(18±0.5)℃;共设置7级参考风速(为风洞过风断面的平均风速,由无级变频器精细控制、数字表盘显示),包含了典型湖泊盛行的风速区间,见表 1. 每级风速进行2组试验,试验结果的相对误差均在5 % 以内,采用2次结果的平均值为最终结果. 根据预试验,风、波、流、雍水等动力现象一般在15 min基本达到稳定,因此,每组试验时长为20 min,最后5 min内对风速与水面波动进行同步采集.

表 1 试验组次设计 Tab. 1 Test group design
2 结果分析 2.1 时均风速分布

从半对数坐标系中玻璃壁面条件时均风速分布(图 2a)可见,靠近壁面处气流受边壁影响形成气流边界层,同一测线不同高度处风速存在差异,风速随着测点高度的增加逐渐增加并存在风速极大值;图中还绘出了近壁面测点风速的拟合线,发现数据符合对数分布. 采用式(1)可回归得到摩阻风速u*与特征粗糙度z0

$ u_{z}=\frac{u_{*}}{\kappa} \ln \left(\frac{z}{z_{0}}\right)=\frac{u_{*}}{\kappa}\left(\ln z-\ln z_{0}\right) $ (1)
图 2 近玻璃壁面(a)与近水面(b)风速的垂向分布 (F=19.5 m) Fig.2 Vertical distribution of wind speed near glass wall (a) and near water surface (b) (F= 19.5 m)

式中,z为水面上方高度,uzz处平均风速. 由图 2a可见拟合线斜率(κ/u*κ为卡门常数,取0.4)随风速增加稍有减小,表明摩阻风速与切应力(τ=ρa ·u*2ρa为空气密度)随风速增加而增加. 从水面条件风速分布(图 2b)可见,时均风速分布与玻璃底壁条件相似,但相同吹程与风速条件下水面风速拟合线斜率大于玻璃壁面条件,表明水面条件下摩阻风速与切应力大于玻璃壁面条件.

2.2 摩阻风速与特征粗糙度

图 3a为两个界面条件下摩阻风速与风速的关系,图中也绘出了Reul等[21]、Liberzon和Shemer[9]与Buckley和Veron[11]水面条件试验结果,其中u10为采用式(1)计算得到的10 m高度处风速[22]. 可见水面条件下试验结果与前人结果较一致;无论是玻璃壁面还是水面条件,摩阻风速与风速均呈正相关,但水面摩阻风速整体大于玻璃壁面条件,小风速时两者差异较小,随着风速的增加,水面条件下摩阻风速随风速的增幅更快,表明水气相互作用后会增加摩阻风速与切应力,即同等风速条件下水气通量传输效率大于玻璃壁面条件下气流与玻璃壁面间的通量传输.

图 3 摩阻风速(a)和特征粗糙度(b)与风速的关系 Fig.3 Relationship of friction wind speed (a) and characteristic roughness (b) with wind speed

从风速与特征粗糙度的关系(图 3b)同样可见,水面特征粗糙度整体大于玻璃壁面条件;玻璃壁面条件下特征粗糙度随风速增加呈现了先增加后趋于平稳的趋势,而水面条件下两者呈正相关,这种差异体现了气-固、气-液的动力作用差异,反映了水气相互作用会增加界面的粗糙程度.

2.3 气流主频与波致气流边界层厚度

从7.7 m/s风速条件下风速能量(S)-频率(f)谱(简称能谱)分布与吹程的关系(图 4a)可见,不同吹程下能谱分布无显著差异,且平坦壁面下风速无主频特征. 从10.5 m吹程下风速能谱分布与风速的关系(图 4b)可见,能谱曲线谱形随着风速的增加整体上移,表明各频率风速分量均有所增加.

图 4 玻璃壁面条件下风速能谱分布与吹程(a)和风速(b)的关系 Fig.4 Relationship of wind speed spectrum distribution with blowing fetch (a) and wind speed (b) under glass wall condition

图 5为水面条件下风速能谱图. 图 5a为10.5 m吹程、12.3 m/s风速条件下,水面上方不同高度处能谱曲线,包括风速全频率能谱图与主频附近的局部能谱图(受风浪与风雍水影响,最低风速测点距时均水面高度为7.5 cm),由图可见气流受风浪的影响,靠近水面一定区域内存在1.6 Hz的主频,该测点处波致气流边界层厚度δwind介于16.5~19.5 cm之间,主频沿高度方向基本保持不变,越靠近水面主频能量越集中,气流受风浪的影响越显著. 需要说明的是δwind的确定是基于水面上方不同高度处风速能谱是否存在明显主频与能量集中现象,当某一高度的谱形没有主峰时,定义该点高度和相邻下方测点高度的平均值为波致气流边界层厚度,由于两风速测点间距3 cm,因此波致气流边界层厚度数据存在1.5 cm的误差.

图 5 水面条件下风速能谱 Fig.5 Wind speed energy spectrum under water surface condition

图 5b为19.5 m吹程、12.3 m/s风速条件下不同高度处风速能谱图,频谱随高度变化规律与图 5a相似,但该测点风速主频为1.2 Hz,小于F16.5测点主频(1.6 Hz),该测点下波致气流边界层厚度介于19.5~22.5 cm之间,大于F16.5测点边界层厚度(16.5~19.5 cm). 表明同等风速条件下随着吹程的增加,近水面气流主频降低而波致气流边界层厚度增加,是有限水域近水面气流沿程演化特征之一.

图 5c为19.5 m吹程、7.7 m/s风速条件下不同高度处风速能谱图. 对比图 5b可见,相同吹程时(F19.5),风速由12.3 m/s减低为7.7 m/s后,风速主频1.2 Hz增大至1.77 Hz,波致气流边界层厚度由19.5~22.5 cm降低至12.5~15.5 cm之间. 表明随着风速的减小,近水面气流主频增加而波致气流边界层厚度减小,体现了有限吹程条件下近水面风动力与风速的响应关系,是有限水域近水面气流又一演化特征.

图 6a为近水面风速主频fdwind与风速和吹程的响应关系,黑色数据点为实测值,彩色图为基于实测点的插值结果,由于仪器限制,部分数据缺失,但并不影响对整体趋势的判断,即近水面风速主频随吹程和风速的增加而减小. 从波致气流边界层厚度δwind与风速和吹程的关系(图 6b)可见,波致气流边界层厚度与风速和吹程均呈正相关,体现了水面气流边界层内风动力随风吹程及风速增加的演化特征.

图 6 气流主频(a)和边界层厚度(b)与风速和吹程的关系 Fig.6 Relationship of dominant frequency (a) and boundary layer thickness (b) of airflow with wind speed and blowing distance
3 讨论

湖泊系统涉及水-气和水-泥两个界面,气流是湖泊水体运动的主要驱动力之一,通过影响水体运动强度与模式来营造整个湖泊生态与环境系统. 对湖泊近水面气流特征与水气相互作用开展研究是认识和揭示水气界面物质、能量传递的基础,是物理湖泊学研究的主要内容之一. 然而,由于天然条件的复杂性,难以开展细致、可控的研究. 目前针对水气作用的研究多采用风洞水槽试验的方法,尽管风洞水槽与天然湖泊存在差异,但两个情境下水气相互作用的物理机制是一致的. 相对于野外条件的复杂性,风洞水槽具有的精细可控的优势为水气相互作用系统深入研究提供了有利条件,是天然条件下水气相互作用研究的重要且必要的补充. 下面针对前文水气基本演化特征,从水面气流动力参数的合理选取和气流-风浪互馈模式两方面展开讨论,分析了现有研究存在不足的原因,提出了更加合理的表征水面气流动力强度的参数和波致气流边界层厚度的确定方法.

3.1 水面气流动力参数的合理选取

气流动力条件对湖泊系统影响举足轻重,其中对气流动力参数选取不统一是水气传输速率不统一的主要原因之一. 根据前文试验结果(图 2),在一定的风速、吹程与水深条件下,同一测点处水面上方不同高度处风速存在差异,这就必然导致不同研究者选取风速测点高度时存在不确定性,即同一测点风速具有非唯一性,造成风速-气体传输速率间存在区别. 不少学者根据风速廓线式(1)推算标准的10 m高度处风速u10,并建立u10-气体传输速率间的关系,但这种看似标准化的做法依然存在不妥.

以实测水面上方3 m高度处风速u3为例,根据式(1)计算u10的关系式为(2),可见式(2)中存在未知量z0,该参数表征了水面粗糙程度,显然该值与气流和水面特征密切相关,也就表明z0在不同风速与测点条件下存在不同(参见图 4),即z0在湖面具有时空异质性. 这就解释了尽管许多学者建立了u10与气体传输速率的关系却依然与其他研究者存在差异的原因.

$ u_{10}=u_{3} \frac{\ln 10-\ln z_{0}}{\ln 3-\ln z_{0}} $ (2)

基于上述近水面气流特征与演化规律的分析,认为目前以风速作为气流动力参数构建与气体传输速率的关系存在不妥,以摩阻风速或切应力作为气流动力参数更具合理性. 因为它们表征了水气界面作用强度,更具力学意义. 且在一定的水深、吹程与风速条件下,摩阻风速和切应力均具有唯一性,因此,上述两个参数的使用会降低气体传输速率确定的不确定性. 摩阻风速与切应力的获取需要在近水面同时布设两个以上不同高度的风速测点,且两个测点均需在波致气流边界层内(具体波致气流边界层厚度参见3.2节),而后根据式(1)拟合得到摩阻风速,进一步计算得到切应力.

考虑到天然水域条件的复杂性与不确定性,确实存在不具备测量摩阻风速、切应力的情况. 因此,如能建立z0与风速和吹程的关系,进一步根据式(2)计算得到u10,可避免根据主观经验选取z0带来的不确定性. 在建立关系时为增加所选用数据的代表性,除选用本次试验数据外,数据集中还包括了室内风洞水槽试验,湖泊、海峡等有限水域以及外海3种尺度的数据,所有数据均为实际测量值(表 2),水深d介于0.15~59.00 m之间,风速u10介于5.0~24.4之间,吹程F介于4.5~6100000.0 m之间.

表 2 统计数据集 Tab. 2 Statistical data set

文献[22]指出水气相互作用还会还受到水深影响,因此在探寻z0的关系时考虑了风速、吹程和水深3个因素. 采用非线性回归分析方法,得到了ln z0u10Fd的回归关系,见式(3),相关系数为0.80. 图 7为ln z0实测值与根据式(3)得到的计算值的对比,可见不同研究者数据基本分布在45°线两侧,表明关系式(3)的合理性. 由式(3)可见ln z0与风速和吹程正相关,而与水深负相关,即风速和吹程越大,水面粗糙度越大;水深越小,水面粗糙度越大,这一定性结论与文献[22]一致,进一步佐证了式(3)的合理性.

图 7 计算ln z0与实测ln z0的对比 Fig.7 Comparison between the calculated ln z0 and the measured ln z0

在进行有限水域水面动力参数选取时,可首先确定研究点位的吹程和水深,将其代入式(3),而后将式(3)代入式(2),进一步将水面上方z高度处的实测风速uz代入式(2)(代替u3),最后采用试算法获得u10.

$ \ln z_{0}=2.303-\frac{12.649}{\left(0.143+0.184 u_{10}^{0.564} \cdot F^{0.064} \cdot d^{-0.037}\right)^{0.5}} $ (3)
3.2 气流-风浪互馈模式

有限吹程水域风浪随吹程逐渐演化,对应地,水面上方气流也进行着动力调整,两者在演进过程中相互作用并逐渐耦合. 图 8为吹程为19.5 m、风速为7.7 m/s条件下同步测量的风速与水面波动的变化过程(已减去时均值),可见两者存在一定的响应关系. 统计了相同风速及吹程条件下近水面气流主频与水面波动主频的比值fdwind/fdwave,建立与吹程的关系(图 9). 由图 9可见,fdwind/fdwave接近1,揭示了气流与风浪“同频率”的互馈模式,且这种模式与吹程和风速是无关的.

图 8 吹程19.5 m、风速7.7 m/s时风速与水面波动的变化过程 Fig.8 Change process of wind speed and water surface fluctuation when blowing fetch is 19.5 m and wind speed is 7.7 m/s
图 9 不同风速时fdwind/fdwave与吹程的关系 Fig.9 Relationship between fdwind/fdwave and blowing fetch at different wind speeds

气流与风浪相互作用除了体现在“同频率”外,还需考虑风浪直接影响的气流高度,并分析其与风浪尺度的关系,可为水面上方测控设备(如风速仪)布设高度的确定提供指导. 根据实测的波致气流边界层(图 6b)及对应主频波高Hd与主频波长Ld(HdLd的计算参考文献[20]),不同风速下δwind/Hdδwind/Ld随吹程变化见图 10,同时绘出各吹程平均值及误差线,结果表明δwind/Hdδwind/Ld与风速关系不显著,小吹程时δwind/Hdδwind/Ld均呈增加趋势,吹程大于10.5 m后基本趋于稳定,稳定后δwind约为Hd的3.3倍,约为Ld的0.28倍,与Liberzon和Shemer[9]的试验结果较接近(0.32倍Ld). 但本文在前人基础上,进一步揭示了气流边界层厚度与风浪尺度的关系并非一成不变,而是与吹程呈现一定关系,即在水气作用初期,波致气流边界层厚度增幅较风浪尺度快,当水气作用一段距离后,波致气流边界层厚度增幅与风浪尺度增幅相当,体现了气流与风浪作用的非线性特征.

图 10 不同风速下δwind/Hd (a) 和δwind/Ld (b) 与吹程的关系 Fig.10 Relationship of δwind/Hd (a) and δwind/Ld (b) with blowing distance at different wind speeds

尽管图 10在一定程度上体现了水气互馈过程中的气流边界层发展与风浪特征的关系,但整体而言数据离散程度较大,表明仅采用风浪的某一特征参数难以较全面地体现风浪特征. 考虑到非线性是风浪的重要特征之一,与风浪变形及破碎等现象密切相关,决定了风浪与气流及底壁的能量及物质交换强弱[30]. 因此,考虑采用风浪非线性参数表征风浪特征,并建立其与波致气流边界层厚度的关联. 无量纲数Fc(Fc=g1.25H0.5T2.5/d1.75HT分别表示平均波高和波周期)反映了波高、波周期和水深的综合影响,将其作为风浪的非线性特征参数[31]. 不同风速与吹程时风浪的Fc分布见图 11(Fc介于2~45之间),可见随着风速和吹程的增加,Fc值呈非线性增加趋势,吹程及风速较小时Fc增长速率较慢,吹程及风速较大时增幅较快,即风浪的非线性特征随着风速和吹程增加逐渐增强.

图 11 Fc与风速和吹程的关系 Fig.11 Relationship of Fc with wind speed and blowing distance

考虑水深对风浪尺度及气流边界层的潜在影响,为使试验结果不失一般性,采用水深dδwind进行无量纲化,建立无量纲的气流边界层厚度δwind/d,其与风浪非线性参数的关系见图 12,可见不同风速时数据点分布较紧凑(显著优于图 10),变化规律较一致,δwind/dFc增加呈现非线性增加趋势,在风浪非线性较弱时(Fc较小)增幅较快,随着风浪非线性增强,波致风浪边界层厚度的增幅逐渐放缓,有趋于饱和之势. 由于Fc表征风浪类型,当Fc小于10时风浪为深水波[31],此时底壁对风浪的约束作用可忽略,水气互馈作用不受底壁影响;而当Fc大于10并逐渐增加时,底壁对风浪的约束逐渐凸显,抑制风浪形态与动力的发展,进一步限制气流边界层发展,表明有限水域中水气互馈作用还会受到水深(底壁)影响. 此外,由图 12还可以看出,当Fc大于20时,δwind/d大于1,表明此时波致气流边界层厚度已大于水深.

图 12 无量纲波致气流边界层厚度与风浪非线性参数的关系 Fig.12 Relationship of dimensionless wave and induced airflow boundary layer thickness with wave nonlinear parameters

采用对数函数对数据点拟合得到的关系式(4),相关系数为0.95. 表明Liberzon和Shemer[9]、Sullivan等[16]、Zavadsky和Shemer[17]等仅采用Ld建立与δwind的关系具有片面性,本文提出采用非线性参数表征风浪特征,综合考虑了波高、波周期与水深等多个因素,在建立其与波致气流边界层厚度间的关系时更具合理性.

$ \delta_{\text {wind }} / d=0.368 \ln F_{\rm{c}}-0.069 $ (4)

基于上述波致气流边界层与风浪特征非线性关系的认识,在湖泊水面仪器布设时,应充分考虑当地风浪与水深条件,参考式(4)推算出波致气流边界层厚度,根据仪器功能确定其布设高度,如仪器需要测量水气作用较强烈区域的数据(如近水面风速),则需将其布设在波致气流边界层内.

4 结论

本文在风洞水槽中模拟研究了有限水域条件下近水面气流演化特征,得到主要结论如下:

1) 玻璃壁面与水面条件时近壁面气流均符合对数分布;水面条件下摩阻风速与特征粗糙度均与风速呈正相关,且同等风速条件下均大于光滑壁面时的值.

2) 风浪上方一定高度内气流具有主频特征,主频与风速和吹程呈负相关;波致气流边界层厚度与风速和吹程均呈正相关.

3) 分析指出了以风速作为气流动力参数构建其与气体传输速率关系存在不妥,阐述了以摩阻风速或切应力作为气流动力参数更具合理性的原因,建立了特征粗糙度与风速、吹程和水深的回归表达式,可为水气界面气体传输效率研究提供参考.

4) 近水面气流与风浪具有“同频率”的互馈模式;指出目前普遍采用波长(表征风浪尺度)建立与波致气流边界层厚度的关系具有差异性与片面性,认为采用综合考虑波高、波周期和水深影响的风浪非线性参数Fc建立与波致气流边界层厚度间的关系更具合理性,给出了δwind/dFc的经验关系式,可供应用参考.

5 参考文献

[1]
Ho JC, Michalak AM, Pahlevan N. Widespread global increase in intense lake phytoplankton blooms since the 1980s. Nature, 2019, 574(7780): 667-670. DOI:10.1038/s41586-019-1648-7
[2]
Qin BQ. Shallow lake limnology and control of eutrophication in LakeTaihu. J Lake Sci, 2020, 32(5): 1229-1243. [秦伯强. 浅水湖泊湖沼学与太湖富营养化控制研究. 湖泊科学, 2020, 32(5): 1229-1243. DOI:10.18307/2020.0501]
[3]
Jiménez PA, Dudhia J. On the need to modify the sea surface roughness formulation over shallow waters. Journal of Applied Meteorology and Climatology, 2018, 57(5): 1101-1110. DOI:10.1175/jamc-d-17-0137.1
[4]
Zhang YF, Yang P, Zhao GH et al. Diffusive CO2 flux across the water-air interface of reclaimed shrimp ponds in the Minjiang River estuary based on the TBL model. Environmental Science, 2019, 40(3): 1188-1197. [张逸飞, 杨平, 赵光辉等. 基于TBL模型的闽江口围垦养虾塘水-大气界面CO2扩散通量估算. 环境科学, 2019, 40(3): 1188-1197. DOI:10.13227/j.hjkx.201808025]
[5]
Xiao QT, Zhang M, Hu ZH et al. Estimate of gas transfer velocity between water-air interface in a large lake based on different models: A case study of Lake Taihu. J Lake Sci, 2018, 30(3): 790-801. [肖启涛, 张弥, 胡正华等. 基于不同模型的大型湖泊水-气界面气体传输速率估算. 湖泊科学, 2018, 30(3): 790-801. DOI:10.18307/2018.0321]
[6]
Karimpour A, Chen Q, Twilley RR. Wind wave behavior in fetch and depth limited estuaries. Scientific Reports, 2017, 7: 40654. DOI:10.1038/srep40654
[7]
Zhang YH, Hu WP. Wind wave characteristics of Lake Chaohu under the influence of typhoon. J Lake Sci, 2020, 32(1): 236-245. [张怡辉, 胡维平. 台风影响下巢湖风浪特征研究. 湖泊科学, 2020, 32(1): 236-245. DOI:10.18307/2020.0122]
[8]
Wu TF, Qin BQ, Ding WH et al. Field observation of different wind-induced basin-scale current field dynamics in a large, polymictic, eutrophic lake. Journal of Geophysical Research: Oceans, 2018, 123(9): 6945-6961. DOI:10.1029/2018jc014040
[9]
Liberzon D, Shemer L. Experimental study of the initial stages of wind waves' spatial evolution. Journal of Fluid Mechanics, 2011, 681: 462-498. DOI:10.1017/jfm.2011.208
[10]
Longo S. Wind-generated water waves in a wind tunnel: Free surface statistics, wind friction and mean air flow properties. Coastal Engineering, 2012, 61: 27-41. DOI:10.1016/j.coastaleng.2011.11.008
[11]
Buckley MP, Veron F. Structure of the airflow above surface waves. Journal of Physical Oceanography, 2016, 46(5): 1377-1397. DOI:10.1175/jpo-d-15-0135.1
[12]
Zhang T, Song JB. Effects of sea-surface waves and ocean spray on air-sea momentum fluxes. Advances in Atmospheric Sciences, 2018, 35(4): 469-478. DOI:10.1007/s00376-017-7101-7
[13]
Esters L, Breivik, Landwehr S et al. Turbulence scaling comparisons in the ocean surface boundary layer. Journal of Geophysical Research: Oceans, 2018, 123(3): 2172-2191. DOI:10.1002/2017jc013525
[14]
Barlow JB, Rae WH, Pope A eds. Low-speed wind tunnel testing. New Jersey: John Wiley & Sons Inc., 1999.
[15]
Shaikh N, Siddiqui K. Airside velocity measurements over the wind-sheared water surface using particle image velocimetry. Ocean Dynamics, 2008, 58(1): 65-79. DOI:10.1007/s10236-008-0132-y
[16]
Sullivan PP, McWilliams JC, Moeng CH. Simulation of turbulent flow over idealized water waves. Journal of Fluid Mechanics, 2000, 404: 47-85. DOI:10.1017/s0022112099006965
[17]
Zavadsky A, Shemer L. Characterization of turbulent airflow over evolving water-waves in a wind-wave tank. Journal of Geophysical Research: Oceans, 2012, 117(C11): 1-21. DOI:10.1029/2011jc007790
[18]
Wang Y, Wang H, Ji W et al. Experimental study on hydrodynamic characteristics of wind-driven flow. South-to-North Water Transfers and Water Science & Technology, 2018, 16(1): 122-126. [王勇, 王华, 纪伟等. 风生流水动力特性试验研究. 南水北调与水利科技, 2018, 16(1): 122-126.]
[19]
Gao A, Wu SQ, Chen L et al. Experimental study on the wave height distribution of wind-induced waves in the growth stage under finite water depth. Water Supply, 2019, 19(8): 2270-2279. DOI:10.2166/ws.2019.108
[20]
Gao A, Wu SQ, Wu XF et al. Experimental study on the influence of drift flow on the propagation characteristics of wind waves in shallow lakes. Advances in Water Science, 2020, 31(4): 592-600. [高昂, 吴时强, 吴修锋等. 滑移流对浅水湖泊风浪传播特性影响试验研究. 水科学进展, 2020, 31(4): 592-600.]
[21]
Reul N, Branger H, Giovanangeli JP. Air flow structure over short-gravity breaking water waves. Boundary-Layer Meteorology, 2008, 126(3): 477-505. DOI:10.1007/s10546-007-9240-3
[22]
Jones Ian SF, Toba Y eds. Wind stress over the ocean. Cambridge: Cambridge University Press, 2001.
[23]
Ding Y. The influence of wind sea state on wind stress[Dissertation]. Qingdao: Ocean University of China, 2005. [丁赟. 风浪状态对海面风应力的影响[学位论文]. 青岛: 中国海洋大学, 2005. ]
[24]
Caulliez G, Makin V, Kudryavtsev V. Drag of the water surface at very short fetches: Observations and modeling. Journal of Physical Oceanography, 2008, 38(9): 2038-2055. DOI:10.1175/2008jpo3893.1
[25]
Mahrt L, Vickers D, Howell J et al. Sea surface drag coefficients in the Risø Air Sea Experiment. Journal of Geophysical Research: Oceans, 1996, 101(C6): 14327-14335. DOI:10.1029/96jc00748
[26]
Johnson HK, Højstrup J, Vested HJ et al. On the dependence of sea surface roughness on wind waves. Journal of Physical Oceanography, 1998, 28(9): 1702-1716. DOI:10.1175/1520-0485(1998)0281702:otdoss>2.0.co;2
[27]
Donelan MA, Babanin AV, Young IR et al. Wave-follower field measurements of the wind-input spectral function. part Ⅱ: Parameterization of the wind input. Journal of Physical Oceanography, 2006, 36(8): 1672-1689. DOI:10.1175/jpo2933.1
[28]
Smith SD. Wind stress and heat flux over the ocean inGale force winds. Journal of Physical Oceanography, 1980, 10(5): 709-726. DOI:10.1175/1520-0485(1980)0100709:wsahfo>2.0.co;2
[29]
Geernaert GL, Larsen SE, Hansen F. Measurements of the wind stress, heat flux, and turbulence intensity during storm conditions over the North Sea. Journal of Geophysical Research Atmospheres, 1987, 92(C12): 13127. DOI:10.1029/jc092ic12p13127
[30]
Teixeira MAC, Argaín JL, Xu X. Reflection of nonlinear mountain waves by critical levels: Behaviour of the reflection coefficient. Quarterly Journal of the Royal Meteorological Society, 2020, 146(727): 1009-1025. DOI:10.1002/qj.3722
[31]
Nelson RC. Depth limited design wave heights in very flat regions. Coastal Engineering, 1994, 23(1/2): 43-59. DOI:10.1016/0378-3839(94)90014-0