基于卷积神经网络的Bolund岛流场重构.pdf

文章编号：0258-1825(2023)09-0038-10基于卷积神经网络的 Bolund 岛流场重构李怡锋1，2，薛明2，吴定聪1，2，张伟1，倪章松2，*(1.西南科技大学土木工程与建筑学院，绵阳621010；2.成都流体动力创新中心，成都610072)摘要：基于卷积神经网络的流场重构方法是风资源评估领域的一种新兴方法，具有成本低、周期短、精度高等优点，在风电场建设评估和风资源利用等领域具有较强的学术价值和工程应用背景。以风资源评估中常用的 Bolund 岛为研究对象，采用 CFD 数值模拟的输入风廓线和输出流场为样本标签，构建了以多层融合卷积神经网络为基础的流场重构方法。以缩短计算时间为目标，研究了样本量对流场重构方法精度的影响，同时对重构方法的外插能力进行了分析。研究结果表明：建立的基于卷积神经网络的流场重构方法可以较好地还原流场信息；当重构目标处于训练样本范围内，全场最大重构误差不高于 2%，且在样本量从 100 减小到 25 时，全场最大重构误差仍在 5%以内；构建的重构方法具有一定的外插能力，但其能力有限，当重构目标超出训练样本范围，全场最大重构误差增加到 20%以上。关键词：风资源评估；复杂地形；风廓线；深度学习；卷积神经网络；流场重构中图分类号：TP183；P425.3文献标识码：Adoi：10.7638/kqdlxxb-2023.0012Flow field reconstruction for Bolund Island based on CNNLIYifeng1，2，XUEMing2，WUDingcong1，2，ZHANGWei1，NIZhangsong2，*(1.School of Civil Engineering and Architecture,Southwest University of Science and Technology,Mianyang621010,China；2.Chengdu Fluid Dynamics Innovation Center,Chengdu610031,China)Abstract:Flowfieldreconstructionbasedontheconvolutionalneuralnetwork(CNN)isanewmethodinthefieldofwindresourceassessment.Ithastheadvantagesoflowcost,shortevaluationandhighprecisionforthewindfarmconstructionandwindresourceutilization,whichprovidesvaluesinbothacademicresearchandengineeringpractice.ByselectingBolundIslandasacasestudy,theflowfieldreconstructionmethodbasedonamulti-layerfusionconvolutionneuralnetworkhasbeendeveloped,whichadoptssamplelabelsoftheinputwindprofileandtheoutputflowfieldasusedintheCFD.Aimingatshorteningthecomputationaltime,thisstudyexaminestheinfluenceofthesamplenumberontheaccuracyoftheflowfieldreconstructionmethod,andanalyzesitsextrapolationabilityinthemeantime.TheresultssuggestthattheflowfieldreconstructionmethodbasedonCNNcanrestoretheflowfieldinformationwell.Themaximumreconstructionerrorislessthan2%whenthereconstructiontargetiswithinthetrainingsamplerange,anditkeepswithin5%whenthenumberofsamplesisreducedfrom100to25.Thereconstructionmethodhassomeextrapolationability,butislimitedwithin some range.When the reconstruction target is beyond the training sample range,the maximumreconstructionerrorrisesfromlessthan5%tomorethan20%.Keywords:windresourceassessment；complexterrain；windprofile；deeplearning；convolutionalneuralnetwork；flowfieldreconstruction收稿日期：2023-02-06；修订日期：2023-03-14；录用日期：2023-03-16；网络出版时间：2023-04-07基金项目：国家自然科学基金（11902278）作者简介：李怡锋（1997），男，四川攀枝花人，硕士研究生，研究方向：风工程.E-mail：通信作者：倪章松*（1973），男，安徽太湖人，博士，正高级工程师，研究方向：低速空气动力学，飞机结冰与防除冰技术.E-mail：引用格式:李怡锋,薛明,吴定聪,等.基于卷积神经网络的 Bolund 岛流场重构J.空气动力学学报,2023,41(9):3847.LIYF,XUEM,WUDC,etal.FlowfieldreconstructionforBolundIslandbasedonCNNJ.ActaAerodynamicaSinica,2023,41(9):3847(inChinese).doi:10.7638/kqdlxxb-2023.0012第41卷第9期空气动力学学报Vol.41，No.92023年9月ACTA AERODYNAMICA SINICASept.，2023 0 引言风能作为可再生的清洁能源，受到世界各国的重视，并已在全球大规模开发利用。“大力提升风电规模与基础建设”作为一项重要内容被写入“十四五”规划和 2035 年远景目标纲要1。据中商产业研究院预测，2025 年我国风力发电量将接近十万亿千瓦时2。在风电开发中，风电场选址与风资源评估尤为重要，其决定着风电场建设的经济效益与风资源利用率3-5。传统的风资源评估主要有三种：基于气象历史观测资料的评估，基于测风塔观测数据的评估，基于数值计算的风资源评估6。受限于当时的技术手段，早期的研究只能采用基于气象历史观测资料的评估方法来进行风能资源评估。随着气象观测手段的丰富，基于测风塔观测数据的评估方法得到了发展，能够较为精确地计算当地的风资源。但这两种评估方法受限于测量器材的覆盖范围和风资源的季节波动，都会对风能资源评估结果的准确度产生一定的影响7-8。基于数值计算的风资源评估方法可以较准确地获得计算区域内的风能资源分布趋势，但模拟的风速在数值上会有系统性偏差，可以利用测风塔观测数据和气象站观测数据进行修正，能够获得较准确的区域风能资源分布9。但这三种方法都还有一定的缺陷，前两种评估方法需要进行长年累月的观测数据的收集，而后一种方法会耗费大量的计算资源，因此需要开展快速实现风资源评估的方法研究。随着计算机技术的发展，基于深度学习、神经网络等方法的流场重构技术在风资源评估中表现出了巨大潜力。利用测风塔观测、CFD 数值模拟等海量数据进行流场重构，在风电场选址和风资源评估过程中可以有效节约成本与周期，成为风资源评估的新兴研究方向之一。卷积神经网络（convolutionalneuralnetwork,CNN）作为一种深度学习模块，具有学习变量之间复杂非线性关系的能力10，研究人员以此开展了流场重构工作。针对风资源评估的研究领域，Li 等11提出了一种双边卷积神经网络，基于高可信度数值模拟得到的风电场尾迹数据集，建立了时间与尾迹流场的映射关系，研究表明，训练完成后的双边卷积神经网络不仅能重构风电场尾迹流场信息，还可以实现准确的尾迹预测，其利用卷积神经网络对风速的总体预测误差为 3.7%。Li 等12还利用卷积神经网络在较短的时间内实现了速度场的重构，能实时实现风电场尾迹的准确预测，所开发的模型对自由风速的总体预测误差为 3.7%，研究工作可为风电场布局优化提供指导。Zhang 等13提出了一种基于物理驱动的深度学习模型，针对风力涡轮机前三维时空风场，建立了激光雷达测得的风速与时间同步对应的数据集，研究结果表明，训练好的卷积神经网络可以实现三维时空风场的重建，流向、横向和垂直速度场的预测误差分别为0.263m/s、0.397m/s 和 0.361m/s，该项研究可应用在风力涡轮机控制和监测、风资源评估等领域。从以上学者的研究，可以发现基于 CNN 的流场重构方法在风资源评估的应用中极具潜力。本文基于卷积神经网络开展了 Bolund 岛的流场重构方法研究，技术路线如图 1 所示。以数值模拟的风廓线输入和输出流场作为样本标签，利用 CNN 建立的模型进行训练学习，建立了基于风廓线和流场的映射关系模型的流场重构方法，并研究了样本量对流场重构方法精度的影响，同时对重构方法的外插能力进行了研究分析。所构建的基于卷积神经网络的流场重构方法可有效缩短风资源评估周期，节约风资源开发利用成本，可为未来深度学习和神经网络方法应用于风资源评估提供一定参考。基于卷积神经网络的流场重构方法生成训练输出输入建立改变样本数量映射关系模型风廓线CFD 计算速度场重构精度外插能力边界条件速度卷积神经网络图 1 技术路线图Fig.1 Technique roadmap 1 风场的数值计算方法基于卷积神经网络构建的映射关系模型是流场重构方法的基础。训练模型所需数据通过数值模拟获得，因此数值计算方法的正确性是流场重构的重要组成部分。本小节介绍所采用的数值计算方法，并通过 Bolund 岛数值计算结果与测风塔实测数据对比，验证数值模拟的正确性。1.1 数值计算方法测风塔在实际测量时通常采用时均数据。为更好地模拟真实风场环境及减少额外干扰，选取在中性大气层结的条件下开展研究，即大气中的气块运动处于平衡状态。采用雷诺平均的 N-S 方程（RANS）14，第9期李怡锋等：基于卷积神经网络的 Bolund 岛流场重构39k-湍流模型为标准湍流模型和标准壁面函数。运输方程如下所示：(k)t+(ujk)t=xj(+tPrk)kxj+Gk(1)()t+(uj)t=xj(+tPr)xj+G1kGkG22k(2)kPrk=1.0Pr=1.32 G1=1.44G2=1.92 G22/k其中：k 为湍流动能；为耗散率；为由平均速度梯度引起的湍动能生成项；Gk为平均速度梯度引起的湍 动 能；经 验 常 数、；及为耗散项。空间半离散格式采用有限体积法，对流项离散格式为一阶迎风格式。1.2 边界条件在开展 Bolund 绕流数值模拟时，边界条件根据文献 15 设置，其中入口为速度入口边界条件，出口为完全发展自由出流边界条件，侧面与顶面采用自由滑移边界条件，底面为无滑移边界条件。uk按照公式（35）设置速度入口边界条件。入口速度、湍动能 和耗散率 分别为：u=ukvln(zzz0)(3)k=u2CC1ln(zzz0)+C2(4)=u3k(z+z0)C1ln(zzz0)+C2(5)uz0=0.000 3z=0.75 mC=0.09kv=0.4其中：摩擦速度（近地面层湍流运动的一个参考速度）依据相关情况进行设定；地表粗糙长度（定义为粗糙地表上方风速为零的点到地面的距离）；地面海拔高度；经验系数；冯卡门常数。C1C2和为湍流动能及其耗散率廓线的模型常数，采用以下公式确定：C1ln(zzz0)+C2=1,zz=0C1ln(zzz0)+C2=0.000 1,zz=H(6)zz=0u2/Czz=H其中：当时，湍流动能取值为；当时（H 为大气边界层厚度，这里为计算域高度），湍流动能趋近于 0，这里取值为 0.0001。2 基于 Bolund 的数值计算方法验证 2.1 Bolund 模型介绍=0=31Bolund 岛（图 2）地理位置处于哥本哈根厄勒海峡，是风资源评估研究中常用的标准模型。小岛整体为东西走向，东西长 130m，南北长 75m，海拔高度为 12m。小岛西侧为陡峭悬崖，除悬崖壁面之外的地表均被矮草覆盖16。Risoe 实验室的长期研究表明，Bolund 岛常年面临西南风和西风17，他们在岛上布置了十座测风塔（如图 3 所示）。为便于表述，定义风向角 为 x 轴正方向与风速之间的夹角，则西风和西南风可分别表示为和。本文主要针对两个风向进行数值模拟。模拟时，设置岛屿中心处测风塔 M3 为坐标原点，正东方向为 x 轴正方向，正北方向为 y 轴正方向。图 2 Bolund 岛实况图Fig.2 Overall view of Bolund Island100500501001502001000 x/mxy=31=0M0M9M8M4M3M5M6M7M1M2y/m100200300图 3 测风塔布设图17Fig.3 Location diagram of the wind towers 17 2.2 网格划分按照实际地形 11 建模，建立以测风塔 M3 为中心的计算域，计算域大小设置为 800m600m500m。计算域及边界条件设置如图 4 所示。z=0.75 m采用结构化网格，根据文献 18-19确定第一层网格高度为 1.07m，对 Bolund 岛附近局部区域进行加密处理。以处截面为例的网格示意图见图 5。2.3 网格无关性验证在开展 CFD 数值模拟之前，首先进行网格无关性验证。共设置了四套网格，所有网格中边界层第一层高度保持一致。四套网格具体尺寸如表 1 所示。40空气动力学学报第41卷u=0V=u2+v2uv选取摩擦速度=0.4m/s、风向角，分别采用以上四套网格进行数值模拟，计算结果如图 6 所示。图中横坐标为合速度，其中 和 分别为 x 和 y 方向的速度，纵坐标为竖直高度 z。由图可见，在风廓线的发展过程中，模拟结果与实测结果中风廓线的发展趋势一致。024681003691215Grid1,M3Grid2,M3Grid3,M3Grid4,M3实测数据,M3Grid1,M1Grid2,M1Grid3,M1Grid4,M1实测数据,M1V/(ms1)z/m图 6 测风塔 M1 和 M3 位置的风廓线Fig.6 Wind profiles of the wind towers M1 and M3z3 mz 3 m选择测风塔 M3 处的风廓线进行分析。当高度时，四套网格的模拟结果与实测数据的误差分别 为 0.027%、0.0132%、0.0409%、0.0473%；当 高 度时，风廓线数值与实测数据基本一致。在保证计算精度和考虑计算量的条件下，后续研究选取Grid2 的网格设置方案。2.4 数值计算结果与实验对比分析风向角=0时，计算得到图 7 所示的速度云图。可以看出，由于岛体的阻塞效应，岛前悬崖处加速效应明显，悬崖周围边界层变薄，悬崖位置加速特征明显，与文献 20 的研究结论相符。5V:678910 11 12 13 14 15(a)计算域内部速度云图5V:678910 11 12 13 14 15(b)Bolund 位置局部放大图图 7 风向角=0时的速度云图Fig.7 Velocity contours under the flow angle =0=0为了进一步验证数值模拟的可靠性，将实测数据与数值模拟结果进行对比。图 8 给出了风向角时不同高度处的速度分布对比，其中图 8（a）中标注了 Bolund 岛的轮廓线及测风塔位置。从图中可以看出，数值模拟得到的速度分布与测风塔实测数据在变化规律和数值上基本保持一致，仅在某些局部点有一定的波动，其原因可能是在数值计算时速度入口条件进行了简化后与实际速度入口条件有一定偏差。3 基于 CNN 的流场重构方法 3.1 映射关系模型流场重构的目的是建立风廓线到 Bolund 岛绕流OutflowWallBolund zxyVelocity-inletSymmetry图 4 计算域及边界条件设置Fig.4 Computational domain and boundary conditions0.75 m 高度处的水平面网格中心区域加密网格前端断崖处网格z=0.75 m图 5 截面网格示意图z=0.75 mFig.5 Grid diagram at 表 1 网格尺寸设置方案Table 1 Setting of the grid size网格方案名网格数量/万中心网格尺寸/mGrid1546.78250.58 0.29Grid2294.37200.74 0.37Grid3147.02840.90 0.45Grid474.32811.06 0.53第9期李怡锋等：基于卷积神经网络的 Bolund 岛流场重构41速度场的映射关系，其可以表示为一个非线性函数：V=f(U,)(7)VUf其中：为 Bolund 岛绕流速度场；为入口处的速度廓线；为可训练的网络参数；是以 CNN 为基础构建的映射关系模型。当 CNN 完成训练时，便会确定最优的网络参数，即建立了风廓线到 Bolund 岛绕流速度场的映射关系模型。3.2 卷积神经网络设置目前，卷积神经网络在流场重构研究中已得到广泛应用，有大量的网络结构出现。虽然卷积神经网络结构多变，但其始终由卷积、激活、全连接等基本元素组成。CNN 的训练过程就是一个参数优化问题，通过定义损失函数作为网络收敛的判据，网络不断地学习流场特征，最终建立映射关系模型。在 CNN 中设置了 3 层线性层和 9 层卷积层，利用线性层进行从稀疏数据到密集数据的上采样操作，在卷积层中设置了尺寸为 33 的卷积核提取特征。网络结构如图 9 所示。线性层和卷积层的具体设置见表 2。卷积层和线性层的输出都采用 ELU（exponentiallinearunits）激活函数，其表达式如下：f(x)=x,x 0ex1,x0(8)在卷积神经网络采用 ELU 激活函数不仅能够加速收敛，还能避免梯度饱和以及梯度消失的问题21。为了加速模型学习速度和缓解梯度消失22，在最后一层线性层中设置 Sigmoid 函数进行数据归一化处理。Smooth L1用损失函数衡量真实值和预测值之间的差距，函数的表达式如下：Smooth L1(Y|Y)=0.5(YY)2,?YY?1?YY?0.5,otherwise(9)YY式中：为真实值，为预测值。文中将数值模拟得到的流场数据作为真实值，CNN 输出的流场数据作为预测值。通过正向传播输出预测值，利用损失函数计算出预测值和真实值之间的损失值。以最小的损失值为优化目标，模型通过梯度下降法进行反向传播，采用自适应矩估计算法（Adam）23更新网络中的参数，来降低真实值与预测值之间的损失，使得模型生成的预测值向真实值靠拢，从而达到学习的目的。Yj=wijXi+bijXiYjwijbijij本文的卷积神经网络架构中，在路径融合前后都使用了线性层，在线性层中依据函数进行拟合。其中：为输入流场信息；为输出流场 信 息；为 权 重；为 偏 置；的 取 值 范 围 为1510；的取值范围为 11024。具有卷积操作的卷积层是基于 CNN 的流场重构141210V/(ms1)8642001000(a)z=5 m的速度分布z=5 mM7 M6M3M8速度实测数据100200141210V/(ms1)8642001000(b)z=2 m的速度分布z=2 m速度实测数据10020030z/m201002001000(c)Bolund岛在=0的直面截图100200 x/mx/mx/m图 8 不同高度处的速度分布与实测数据对比图Fig.8 Computed velocity distributions at different heightscompared with the measured data表 2 线性层和卷积层的设置Table 2 Settings for linear and convolutional layers类型卷积核尺寸/步长（或注释）输入尺寸输出尺寸线性层21.ELU2.ELU5101102411024116001卷积层933,14040158581线性层ELU33641122881201510z/m5005 10Input:入口处风速廓线Output:平面内的速度场线性层线性层卷积层15V/(ms1)20 25 30图 9 卷积神经网络结构图Fig.9 Structure of the convolutional neural network architecture42空气动力学学报第41卷方法的核心部分。在卷积层中，上一层输出的特征会在卷积层进行卷积运算，用滤波器将特征进行处理，输出新的特征。卷积层中的输出为：Z(l)u,v=li=0lj=0X(l1)i,jK(l)i,j+b(l)(10)x(i,j)=1,0 i,j n0,others(11)lX(l1)K(l)b(l)当输入层是第 层时，输入的流场特征是（mm），特征对应的卷积核是（nn）。为了增加 CNN 的数据拟合能力，输出的流场特征加入偏置单元进行修正。4 Bolund 岛的流场重构研究本节基于数值模拟得到的数据集，在流场重构方法中建立以风廓线为输入、流场为输出的映射关系模型，研究样本量对流场重构方法精度的影响，并探讨样本量对重构方法外插能力的影响。4.1 数据集样本u在样本工况中设置了不同的风向角和摩擦速度。其中摩擦速度可作为地形边界层的特征速度来表征地表摩擦作用，并且不同的摩擦速度对应不同的入口边界条件，具体可参考公式（3）。训练样本及测试样本的设置如表 3 所示。表 3 数据集样本设置Table 3 Setup of the dataset sample风向角/()摩擦速度u*/(ms1)训练样本0，310.01，0.02，0.99，1.00测试样本0，310.01，0.02，1.20，1.21uu训练样本设置为=0.011.0m/s，样本数量 N=100。为了验证 CNN 的外插能力，设置部分测试样本在训练样本范围之外，即=0.011.21m/s。4.2 重构结果分析采取以上样本数据进行映射关系模型的训练，训练过程中的损失值如图 10 所示。在整个训练过程中损失值表现出平稳下降的趋势。在 3000 个训练回合后，训练基本达到收敛状态，建立了流场参数之间的映射关系模型。u=0.75=0u=0.75通过样本训练学习，构建了从输入风廓线到输出流场的映射关系模型，选取m/s（样本范围0.011.21m/s）进行流场重构。图 11 给出了、m/s 时的速度云图及流线图。图中黑色线条为速度流线，白色线条为速度等值线。重构流场和数值模拟的速度流线及速度等值线基本一致，说明重构的速度云图保留了原始流场的信息特征。05001000150020002500300000.0050.0100.0150.0200.0250.030Epoch=0=31Smooth L1_loss图 10 CNN 的训练损失Fig.10 Training loss of CNN20015010050200150100500(a)=0,u*=0.75 m/s,CFD501001502002:4681012141618202224(b)=0,u*=0.75 m/s,CNN20015010050200150100500501001502002V/(ms1)V/(ms1):4681012141618202224z/mz/mx/mx/m图 11 流场重构结果与数值模拟结果对比图Fig.11 Comparison between the reconstructed resultand the CFD result在采用速度绝对值的误差评估流场重构精度的过程中，风廓线最大重构误差及全场最大重构误差为：第9期李怡锋等：基于卷积神经网络的 Bolund 岛流场重构43x=max(|V|V0|V0|)(x,z)(12)=max(x)(13)V=V(x,z)V0=V0(x,z)x其中：为重构流场的速度矢量；为数值模拟流场的速度矢量；为空间坐标 x 固定时流场相对误差的最大值，即风廓线最大重构误差；为全场最大重构误差。x针对不同的风向角选取不同的测风塔位置进行对比，图 12 给出了流场重构结果与数值模拟风廓线对比图，其中 u*=0.75m/s。由图可知，流场重构与数值模拟得到的风廓线十分接近，且风廓线最大重构误差在 1%以内，验证了本文流场重构方法的有效性。05101520250408012016004080120160CNN,M1CNN,M3CFD,M1CFD,M3CNN,M3CNN,M7CFD,M3CFD,M7(a)u*=0.75 m/s,=00510152025(b)u*=0.75 m/s,=31z/mz/mV/(ms1)V/(ms1)图 12 流场重构结果与数值模拟结果的风廓线对比Fig.12 Comparison of the wind profile betweenthe reconstructed result and the CFD result 4.3 外插能力分析uuu=1.11对构建的流场重构方法的外插能力进行检验。选取处于训练样本范围外，即=1.011.21m/s，检验外插能力。图 13 中展示了m/s 时的速度云图及流线图。从流场重构结果中可以发现，速度等值线（白色实线）出现一定的波动，如图 13（b）所示，表明在重构过程中有噪点的产生。但整体上来看，流场重构与数值模拟的速度云图及速度流线十分接近，表明本文的流场重构方法具备一定的外插能力，能够基本实现流场信息的还原，其重构精度还需进一步验证。20015010050200150100500(a)=0,u*=1.11 m/s,CFD501001502002V/(ms1)V/(ms1):4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 362:4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36(b)=0,u*=1.11 m/s,CNN20015010050200150100500 x/mx/m50100150200z/mz/m图 13 流场重构结果与数值模拟结果的速度场对比图Fig.13 Comparison of the velocity contour betweenthe reconstructed result and the CFD resultu=0.011.0u=0.75xu=1.11x为进一步检验外插能力，将流场重构和数值计算结果在不同测风塔处的风廓线进行对比分析，如图 14 所示。图中红色虚线与红色圆形标识分别为M3 位置处流场重构及数值计算得到的风廓线。可以发现，流场重构得到的风廓线与数值模拟的规律基本一致：在近地面，随着高度增加，速度迅速增大；当高度达到一定值后，速度变化趋于平缓。对比图 14 与图 12 可以发现，构建的流场重构方法可以在训练样本范围外进行流场重构，具有一定的外插能力，但外插流场的重构精度有所下降。在训练样本集为m/s 时，样本集范围内m/s 时获得的风廓线最大重构误差小于 1%，而m/s 时迅速上升到 11%。44空气动力学学报第41卷0102030400408012016004080120160V/(ms1)010203040V/(ms1)(a)u*=1.11 m/s,=0(b)u*=1.11 m/s,=31CNN,M1CNN,M3CFD,M1CFD,M3CNN,M7CNN,M3CFD,M7CFD,M3z/mz/m图 14 流场重构结果与数值模拟结果的风廓线对比图Fig.14 Comparison of the wind profile betweenthe reconstructed result and the CFD resultuuuu图 15 给出了全场最大重构误差随摩擦速度的变化趋势。当处于训练样本集范围内，最大重构误差 保持在 2%以内；当重构目标在训练样本集范围外，=1.011.10m/s 时，全场最大重构误差 约为10%，表明在小范围内，构建的流场重构方法在流场重构的过程中能够保持一定的精度，因此具备一定精度的外插能力。从 1.01m/s 增加至 1.21m/s 的过程u=0=31中，全场最大重构误差 快速增大。例如=1.21m/s时，和获得的 分别为 27.29%和 25.77%，表明训练得到的映射关系模型的外插能力在这个条件下已经有一定的减弱，重构精度也随之下降。4.4 样本数量对于重构精度的影响=0样本数量是影响流场重构的重要因素之一。选取风向角时的流场数据，研究样本数量对重构精度的影响。设置样本数量N=100、50、25，如表4 所示。表 4 样本参数Table 4 Sample parameters样本数量风向角/()摩擦速度u*/(ms-1)N=10000.01，0.02，0.99，1.00N=5000.02，0.04，0.98，1.00N=2500.04，0.08，0.96，1.00由表 5 中 CNN 训练得到的损失值，可以看出训练损失与样本数量有极大的关系。在本文研究的样本数量范围内，样本数量从 100 减少至 25 时，模型的训练损失明显变大。表 5 CNN 训练得到的损失值Table 5 Training loss of CNNN=100N=50N=25=02.51 1035.46 10343.3 103u=0.6为了探讨不同样本数量对于重构精度的影响，选取m/s，对比分析重构流场和数值计算得到的速度云图及流线图。由图 16 可知，重构得到的速度流线与数值结果基本一致。但是随着样本数量的减少，速度等值线的波动变得剧烈，表明重构图像的噪点逐渐增多，由此可见样本的数量对于速度云图的重构产生了巨大的影响。u将 CNN 对以上三种不同数量的训练样本学习后建立的映射关系模型进行分析，全场最大重构误差随增长而变化的趋势由图 17 所示。u从图中可知，在训练样本范围内（即=0.011.0m/s）进行流场重构时，即使样本数量从 N=100 减少到 N=25，全场最大重构误差 仍在 5%以内。在训练样本范围外进行流场重构，在同样的样本数量变化范围内，均迅速增大。这与图 15 得到的结论一致。uu在训练样本范围外（即=1.011.21m/s），在=1.21m/s 时最大重构误差达到最大值，样本数量 N=25、50、100 时的分别为 43.75%、32.34%和 27.29%。说明样本数量对于外插能力有着直接的影响。样本数量越少，映射关系模型的外插能力越弱，全场最大重构误差 增加越显著。因此在进行流场重构时，样00.20.40.60.81.01.21.45051015202530/%u*/(ms1)=0=31图 15 最大重构误差随摩擦速度变化的曲线图Fig.15 Graph of maximum reconstruction error varieswith friction velocity第9期李怡锋等：基于卷积神经网络的 Bolund 岛流场重构45本数量的选择要在优先满足最大重构误差的条件下，根据计算资源酌情考虑。5 结论本文基于数值模拟的数据，利用卷积神经网络开展了流场重构方法的研究。通过建立流场参数的映射关系模型，实现了从输入风廓线快速重构输出Bolund 岛表面绕流速度场的目标。根据流场重构结果，对重构精度和外插能力进行了分析，并且基于不同数量样本的学习情况，分析了样本数量对于外插能力的影响。研究得出以下主要结论：u1）构建的流场重构方法能够快速提取流场特征。重构目标处于训练样本范围内（即=0.011.0m/s）、样本数量充足时，全场最大重构误差 不超过 2%，在此范围内能较好地反应真实流场信息。uuu=312）重构目标处于训练样本范围外（即=1.011.21m/s），最大重构误差 快速增大，虽然=1.011.10m/s 的全场最大重构误差 仍保持在 10%左右，但在=1.21m/s 时，=0和的全场最大重构误差 分别为 27.29%和 25.77%。表明在训练样本范围外，通过外插能力能够实现对未知工况的流场重构，但外插能力向外延伸的范围有限。u3）样本数量对于重构精度有一定的影响。当重构目标处于训练样本范围内（即=0.011.0m/s）20015010050200150100500(a)CFDx/mx/mx/mx/m501001502001:2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 181:2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 181:2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 181:2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18(c)N=5020015010050200150100500z/mz/mz/mz/m5010015020020015010050200150100500(b)N=10050100150200(d)N=252001501005020015010050050100150200V/(ms1)V/(ms1)V/(ms1)V/(ms1)图 16 不同数量样本训练得到的流场重构结果与数值模拟结果的速度云图对比图Fig.16 Velocity contours reconstructed using different numbers of samples in comparison with the CFD result00.20.40.60.81.01.21.4010203040N=100N=50N=25/%u*/(ms1)u图 17 全场最大重构误差随增长的变化曲线uFig.17 Variation diagram of the maximum reconstruction errorin the full field with 46空气动力学学报第41卷u时，即使样本数量从 100 减少至 25，构建的流场重构方法仍保持着较好的精度，最大重构误差仍在 5%以内。当 重 构 目 标 处 于 训 练 样 本 范 围 外（即=1.011.21m/s）时，样本数量从 100 减少至 25，最大重构误差从 27.29%迅速增大至 43.75%，重构精度明显降低。需要指出的是，本文选择 Bolund 岛作为研究对象，其地形、地貌等条件已明确，入口风廓线与目标位置的流场具有特定的关联，研究结论仅适用于Bolund 岛的场景。当研究对象改变时，需采用对应的样本数据建立映射关系方可适用。目前本文所构建的流场重构方法与已有的流场重构方法一样，不具备普适性，在物理机制上存在有限度的适用性和合理性。参 考 文 献：中华人民共和国国民经济和社会发展第十四个五年规划和2035年远景目标纲要EB/OL.2021-03-13.http:/ evaluation for a highly complex terrain using ComputationalFluidDynamics(CFD)J.RenewableEnergy,2017,101:19.doi:10.1016/j.renene.2016.08.0324SHARMA P K,WARUDKAR V,AHMED S.Application of a newmethodtodevelopaCFDmodeltoanalyzewi