利用测风塔实测风速数据确定大气稳定度的方法和装置与流程

1.本发明涉及风力发电领域技术领域，具体涉及利用测风塔实测风速数据确定大气稳定度的方法和装置。


背景技术：

2.在大气边界层中，大气稳定度是衡量大气湍流运动强弱的一个重要指标。不同大气稳定度条件下，大气的湍流强度及风切变指数有显著变化，因此准确计算当地环境条件下的大气稳定度对于风能资源评估的准确性以及风力机运行特性的研究具有重要意义。
3.当前，大气稳定度的确定方法主要有莫宁-奥布霍夫长度法、梯度理查森数法、总体理查森数法等，每种方法均需要利用不同高度层的温度和速度等数据来计算大气稳定度。风电场建设过程中一般会设立测风塔，通过在不同高度层安装风速测量装置和温度测量装置来得到风速和温度数据，进而计算大气稳定度。然而，以上几种方法所需计算参数较多，需要同时处理风速和温度数据，计算也较为复杂。


技术实现要素：

4.因此，本发明提供了利用测风塔实测风速数据确定大气稳定度的方法和装置，在确定大气稳定度时仅使用了测风塔的不同高度层风速数据，所需计算参数少，计算方法简单快捷，以解决上述背景技术中提出的问题。
5.为达到上述目的，本发明提供如下技术方案：
6.第一方面，本发明实施例提供利用测风塔实测风速数据确定大气稳定度的方法，包括：
7.提取测风塔预设时间范围内预设高度层的实测风速数据，记为第一时序集合，并基于预设选取方法对其进行选取得到大气环境中性状态的风速数据，记为第二时序集合；
8.对第一时序集合和第二时序集合分别进行预设数据处理，对应得到第一时均对数集合和第二时均对数集合；
9.根据大气边界层不同位置的风速与高度对应的关系式，对第二时均对数集合进行拟合，确定拟合参数；根据拟合参数和关系式，对第一时均对数集合进行拟合，确定风廓线稳定度修正参数；
10.根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式计算莫宁-奥布洛夫长度，并根据莫宁-奥布洛夫长度以及稳定度判别标准确定大气稳定度的级别。
11.优选地，基于预设选取方法对其进行选取得到大气环境中性状态的风速数据的过程，包括：
12.在第一时序集合中，将每个高度层每个时刻的风速数据按照风速大小进行排序，从中选取预设比例的高风速样本数据，得到大气环境中性状态的风速数据。
13.优选地，预设数据处理包括：求时均值处理和取对数处理。
14.优选地，大气边界层不同位置的风速与高度对应的关系式为：
[0015][0016]
其中，u为大气边界层不同位置的风速，z为每个高度层的高度，u
*
为摩擦速度，k为冯卡门常数，z0为地表粗糙度，ψm为风廓线稳定度修正参数。
[0017]
优选地，确定拟合参数的过程，包括：
[0018]
当大气环境处于中性状态时，风廓线稳定度修正参数ψm＝0，则大气边界层不同位置的风速与高度对应的关系式为：
[0019]
将第二时均对数集合带入上述关系式得到拟合参数，拟合参数为和ln z0。
[0020]
优选地，根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式计算莫宁-奥布洛夫长度的过程，包括：
[0021]
根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式，得到莫宁-奥布洛夫长度l的公式为：
[0022][0023]
其中，z为每个高度层的高度。
[0024]
优选地，稳定度判别标准为：
[0025]
当-100≤l《0时，大气环境处于强不稳定状态；
[0026]
当-500≤l《-100时，大气环境处于不稳定状态；
[0027]
当|l|》500时，大气环境处于中性状态；
[0028]
当100≤l《500时，大气环境处于稳定状态；
[0029]
当0≤l《100时，大气环境处于强稳定状态。
[0030]
第二方面，本发明实施例提供利用测风塔实测风速数据确定大气稳定度的装置，包括：
[0031]
数据获取单元，用于提取测风塔预设时间范围内预设高度层的实测风速数据，记为第一时序集合，并基于预设选取方法对其进行选取得到大气环境中性状态的风速数据，记为第二时序集合；
[0032]
数据处理单元，用于对第一时序集合和第二时序集合分别进行预设数据处理，对应得到第一时均对数集合和第二时均对数集合；
[0033]
参数确定单元，用于根据大气边界层不同位置的风速与高度对应的关系式，对第二时均对数集合进行拟合，确定拟合参数；根据拟合参数和关系式，对第一时均对数集合进行拟合，确定风廓线稳定度修正参数；
[0034]
大气稳定度确定单元，用于根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式计算莫宁-奥布洛夫长度，并根据莫宁-奥布洛夫长度以及稳定度判别标准确定大气稳定度的级别。
[0035]
第三方面，本发明实施例提供一种计算机设备，包括：至少一个处理器，以及与至少一个处理器通信连接的存储器，其中，存储器存储有可被至少一个处理器执行的指令，指令被至少一个处理器执行，以使至少一个处理器执行本发明实施例第一方面的利用测风塔实测风速数据确定大气稳定度的方法。
[0036]
第四方面，本发明实施例提供一种计算机可读存储介质，计算机可读存储介质存储有计算机指令，计算机指令用于使计算机执行本发明实施例第一方面的利用测风塔实测风速数据确定大气稳定度的方法。
[0037]
本发明技术方案，具有如下优点：
[0038]
本发明提供的利用测风塔实测风速数据确定大气稳定度的方法和装置，方法包括：提取测风塔预设时间范围内预设高度层的实测风速数据，记为第一时序集合，并基于预设选取方法对其进行选取得到大气环境中性状态的风速数据，记为第二时序集合；对第一时序集合和第二时序集合分别进行预设数据处理对应得到第一时均对数集合和第二时均对数集合；根据大气边界层不同位置的风速与高度对应的关系式，对第二时均对数集合进行拟合，确定拟合参数；根据拟合参数和关系式，对第一时均对数集合进行拟合，确定风廓线稳定度修正参数；根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式计算莫宁-奥布洛夫长度，并根据莫宁-奥布洛夫长度以及稳定度判别标准确定大气稳定度的级别。通过本发明提供的确定大气稳定度的方法和装置，仅使用了测风塔的不同高度层的风速数据，所需计算参数少，计算方法简单快捷。
附图说明
[0039]
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0040]
图1本发明实施例中提供的利用测风塔实测风速数据确定大气稳定度的方法的流程示意图；
[0041]
图2本发明实施例中提供的第一时均集合和第二时均集合的分布图；
[0042]
图3本发明实施例中提供的第二时均对数集合的线性拟合图；
[0043]
图4本发明实施例中提供的第一时均对数集合的线性拟合图；
[0044]
图5本发明实施例中提供的利用测风塔实测风速数据确定大气稳定度的装置的组成图；
[0045]
图6本发明实施例中提供的计算机设备一个具体示例的组成图。
具体实施方式
[0046]
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分的实施例，不是全部的实施例，而并非要限制本发明公开的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要的混淆本发明公开的概念。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实
施例，都应当属于本发明保护的范围。
[0047]
此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
[0048]
实施例1
[0049]
本发明实施例提供利用测风塔实测风速数据确定大气稳定度的方法，如图1所示，该方法包括：
[0050]
步骤s1：提取测风塔预设时间范围内预设高度层的实测风速数据，记为第一时序集合，并基于预设选取方法对其进行选取得到大气环境中性状态的风速数据，记为第二时序集合。
[0051]
在本实施例中，预设时间范围内为一年时间范围内，仅作为举例说明，不以此为限制。
[0052]
实际应用中，测风塔一般会在n个高度层设置风速测量装置，用于实时测量风速，通常n≥5，故在本实施例中，预设高度层的取值范围为大于等于5，具体数据不做限制，依据实际应用做适应性修改。
[0053]
在本实施例中，基于预设选取方法对其进行选取得到大气环境中性状态的风速数据的过程，包括：
[0054]
在第一时序集合中，将每个高度层每个时刻的风速数据按照风速大小进行排序，从中选取预设比例的高风速样本数据，得到大气环境中性状态的风速数据。具体地，预设比例为经验参数，根据环境条件选取，其取值范围为小于30％。需要说明的是，风速越大对应的大气稳定度越趋近于中性，因此认为第二时序集合所对应的大气环境为中性状态。
[0055]
步骤s2：对第一时序集合和第二时序集合分别进行预设数据处理，对应得到第一时均对数集合和第二时均对数集合。
[0056]
在本实施例中，预设数据处理包括：求时均值处理和取对数处理。具体地，对第一时序集合和第二时序集合的风速数据进行求时均值处理，对处理后的第一时序集合和第二时序集合的高度数据进行取对数处理。
[0057]
步骤s3：根据大气边界层不同位置的风速与高度对应的关系式，对第二时均对数集合进行拟合，确定拟合参数；根据拟合参数和关系式，对第一时均对数集合进行拟合，确定风廓线稳定度修正参数。
[0058]
在本实施例中，大气边界层不同位置的风速与高度对应的关系式为：
[0059][0060]
其中，u为大气边界层不同位置的风速，z为每个高度层的高度，u
*
为摩擦速度，k为冯卡门常数，z0为地表粗糙度，ψm为风廓线稳定度修正参数。
[0061]
当大气环境处于中性状态时，风廓线稳定度修正参数ψm＝0，则大气边界层不同位置的风速与高度对应的关系式为：
[0062]
将第二时均对数集合带入上述关系式得到拟合参数，拟合参数为和ln z0。
[0063]
由于第一时均对数集合所对应的大气环境状态未知，即风廓线稳定度修正参数ψm为未知的待求解参数，故将得到的拟合参数为和ln z0以及第一时均对数集合带入大气边界层不同位置的风速与高度对应的关系式，求得风廓线稳定度修正参数ψm。
[0064]
步骤s4：根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式计算莫宁-奥布洛夫长度，并根据莫宁-奥布洛夫长度以及稳定度判别标准确定大气稳定度的级别。
[0065]
在本实施例中，风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式为：
[0066][0067]
根据上述关系式，得到莫宁-奥布洛夫长度l的公式为：
[0068][0069]
其中，z为每个高度层的高度。
[0070]
在本实施例中，将步骤s3中得到的风廓线稳定度修正参数ψm带入上述关系式，求得对应的莫宁-奥布洛夫长度l，根据稳定度判别标准对莫宁-奥布洛夫长度l进行大气稳定度的判别。具体地，稳定度判别标准如下：
[0071]
当-100≤l《0时，大气环境处于强不稳定状态；
[0072]
当-500≤l《-100时，大气环境处于不稳定状态；
[0073]
当|l|》500时，大气环境处于中性状态；
[0074]
当100≤l《500时，大气环境处于稳定状态；
[0075]
当0≤l《100时，大气环境处于强稳定状态。
[0076]
一具体实施例中，对测风塔的不同高度层的实测风速数据进行大气稳定度的确定，具体过程包括：
[0077]
1、提取某一测风塔一年时间范围内5个高度层的风速数据，作为第一时序集合。
[0078]
2、将第一时序集合中的风速数据按照风速大小进行排序，从中选取10％的高风速样本数据，作为第二时序集合。
[0079]
3、对第一时序集合和第二时序集合中的风速数据分别进行求时均值的处理，对应得到第一时均集合和第二时均集合，如图2所示。
[0080]
具体地，第一时均集合记为{(z1,u1),(z2,u2),...,(zn,un)}，第二时均集合记为{(z
h1
,u
h1
),(z
h2
,u
h2
),...,(z
hn
,u
hn
)}，其中，z为每个高度层的高度，u为每个高度层在一年时间范围内风速的时均值，n为高度层。在本实施例中，n＝5，可知：
[0081]
第一时均集合为：
[0082]
{(z1,u1),(z2,u2),(z3,u3),(z4,u4)(z5,u5)}＝{(20,5.38),(50,6.25),(70,6.64),(90,6.88),(100,6.94)}；
[0083]
第二时均集合为：
[0084]
{(z
h1
,u
h1
),(z
h2
,u
h2
),(z
h3
,u
h3
),(z
h4
,u
h4
)(z
h5
,u
h5
)}＝{(20,5.72),(50,6.54),(70,6.90),(90,7.15),(100,7.21)}。
[0085]
4、对第一时均集合和第二时均集合中的高度数据分别进行取对数处理，对应得到第一时均对数集合和第二时均对数集合。
[0086]
具体地，第一时均对数集合记为{(p1,u1),(pi,ui),...,(p n
,un)}，第二时均对数集合记为{(p
h1
,u
h1
),(p
hi
,u
hi
),...,(p
hn
,u
hn
)}，其中，pi＝ln zi,i＝1,2,
…
,n，p
hi
＝ln z
hi
,i＝1,2,
…
,n。在本实施例中，n＝5，可知：
[0087]
第一时均对数集合为：
[0088]
{(p1,u1),(p2,u2),(p3,u3),(p4,u4)(p5,u5)}＝{(3.00,5.38),(3.91,6.25),(4.25,6.64),(4.50,6.88),(4.61,6.94)}；
[0089]
第二时均对数集合为：
[0090]
{(p
h1
,u
h1
),(p
h2
,u
h2
),(p
h3
,u
h3
),(p
h4
,u
h4
)(p
h5
,u
h5
)}＝{(3.00,5.72),(3.91,6.54),(4.25,6.90),(4.50,7.15),(4.61,7.21)}。
[0091]
5、对第二时均对数集合进行线性拟合，确定拟合参数和ln z0。
[0092]
具体地，由于第二时均对数集合对应的大气环境为中性状态，故风廓线稳定度修正参数ψm＝0，对应的大气边界层不同位置的风速与高度对应的关系式为：
[0093]
根据p
hi
＝ln z
hi
,i＝1,2,
…
,n，上述关系式变形为：
[0094]uhi
＝ap
hi-a
×b[0095]
其中，拟合参数为a和b，b＝ln z0。
[0096]
将第二时均对数集合带入上述变形后的关系式进行拟合，确定拟合参数a和b。拟合曲线如图3所示，由图可知：a＝0.941，b＝-3.069。
[0097]
6、根据拟合参数a和b，对第一时均对数集合进行线性拟合，确定风廓线稳定度修正参数ψm。
[0098]
本实施例中，由于第一时均对数集合所对应的大气环境状态未知，根据拟合参数a和b以及pi＝ln zi,i＝1,2,
…
,n，得到对应的大气边界层不同位置的风速与高度对应的关系式为：ui＝ap
i-a(b+ψm)。
[0099]
基于a＝0.941，b＝-3.069和上述关系式，对第一时均对数集合进行线性拟合，确定风廓线稳定度修正参数ψm。拟合曲线如图4所示，由图可知：ψm＝0.303。
[0100]
7、根据莫宁-奥布洛夫长度l的公式：
[0101][0102]
当每个高度层的高度z＝90m时，将风廓线稳定度修正参数ψm＝0.303代入求解，得到莫宁-奥布洛夫长度l＝-785。
[0103]
8、根据莫宁-奥布洛夫长度l＝-785，参照稳定度判别标准可知，当前的大气环境
处于中性状态。
[0104]
综上，本发明提供的利用测风塔实测风速数据确定大气稳定度的方法，仅需用到测风塔的不同高度层的风速数据，所需计算参数少，计算方法简单快捷。
[0105]
实施例2
[0106]
本发明实施例提供利用测风塔实测风速数据确定大气稳定度的装置，如图5所示，包括：
[0107]
数据获取单元，用于提取测风塔预设时间范围内预设高度层的实测风速数据，记为第一时序集合，并基于预设选取方法对其进行选取得到大气环境中性状态的风速数据，记为第二时序集合；此模块执行实施例1中的步骤s1所描述的方法，在此不再赘述。
[0108]
数据处理单元，用于对第一时序集合和第二时序集合分别进行预设数据处理，对应得到第一时均对数集合和第二时均对数集合；此模块执行实施例1中的步骤s2所描述的方法，在此不再赘述。
[0109]
参数确定单元，用于根据大气边界层不同位置的风速与高度对应的关系式，对第二时均对数集合进行拟合，确定拟合参数；根据拟合参数和关系式，对第一时均对数集合进行拟合，确定风廓线稳定度修正参数；此模块执行实施例1中的步骤s3所描述的方法，在此不再赘述。
[0110]
大气稳定度确定单元，用于根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式计算莫宁-奥布洛夫长度，并根据莫宁-奥布洛夫长度以及稳定度判别标准确定大气稳定度的级别；此模块执行实施例1中的步骤s4所描述的方法，在此不再赘述。
[0111]
通过本发明提供的利用测风塔实测风速数据确定大气稳定度的装置，在确定大气稳定度时仅需用到测风塔的不同高度层的风速数据，所需计算参数少，计算方法简单快捷。
[0112]
实施例3
[0113]
本发明实施例提供一种计算机设备，如图6所示，包括：至少一个处理器601，至少一个通信接口603，存储器604和至少一个通信总线602。其中，通信总线602用于实现这些组件之间的连接通信，通信接口603可以包括显示屏和键盘，可选通信接口603还可以包括标准的有线接口、无线接口。存储器604可以是高速易挥发性随机存取存储器，也可以是非不稳定的存储器，还可以是至少一个位于远离前述处理器601的存储装置。其中处理器601可以执行实施例1的利用测风塔实测风速数据确定大气稳定度的方法。存储器604中存储一组程序代码，且处理器601调用存储器604中存储的程序代码，以用于执行实施例1的利用测风塔实测风速数据确定大气稳定度的方法。
[0114]
其中，通信总线602可以是外设部件互连标准(peripheral component interconnect，简称pci)总线或扩展工业标准结构(extended industry standard architecture，简称eisa)总线等。通信总线602可以分为地址总线、数据总线、控制总线等。为便于表示，图6中仅用一条线表示，但并不表示仅有一根总线或一种类型的总线。
[0115]
其中，存储器604可以包括易失性存储器(volatile memory)，例如随机存取存储器(random access memory，简称ram)；存储器也可以包括非易失性存储器(non-volatile memory)，例如快闪存储器(flash memory)，硬盘(hard disk drive，简称hdd)或固态硬盘(solid-state drive，简称ssd)；存储器604还可以包括上述种类的存储器的组合。
[0116]
其中，处理器601可以是中央处理器(central processing unit，简称cpu)，网络
处理器(network processor，简称np)或者cpu和np的组合。
[0117]
其中，处理器601还可以进一步包括硬件芯片。上述硬件芯片可以是专用集成电路(application-specific integrated circuit，简称asic)，可编程逻辑器件(programmable logic device，简称pld)或其组合。上述pld可以是复杂可编程逻辑器件(complex programmable logic device，简称cpld)、现场可编程逻辑门阵列(field programmable gate array，简称fpga)、通用阵列逻辑(generic array logic，简称gal)或其任意组合。
[0118]
可选地，存储器604还用于存储程序指令。处理器601可以调用程序指令，实现如本发明执行实施例1中的利用测风塔实测风速数据确定大气稳定度的方法。
[0119]
本发明实施例还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机可执行指令，该计算机可执行指令可执行实施例1的利用测风塔实测风速数据确定大气稳定度的方法。其中，所述存储介质可为磁碟、光盘、只读存储记忆体(read only memory，简称rom)、随机存储记忆体(random access memory，简称ram)、快闪存储器(flash memory)、硬盘(hard disk drive，简称hdd)或固态硬盘(solid state drive，简称ssd)等；所述存储介质还可以包括上述种类的存储器的组合。
[0120]
显然，上述实施例仅是为清楚地说明所作的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

技术特征：
1.一种利用测风塔实测风速数据确定大气稳定度的方法，其特征在于，包括：提取测风塔预设时间范围内预设高度层的实测风速数据，记为第一时序集合，并基于预设选取方法对其进行选取得到大气环境中性状态的风速数据，记为第二时序集合；对所述第一时序集合和第二时序集合分别进行预设数据处理，对应得到第一时均对数集合和第二时均对数集合；根据大气边界层不同位置的风速与高度对应的关系式，对所述第二时均对数集合进行拟合，确定拟合参数；根据所述拟合参数和所述关系式，对所述第一时均对数集合进行拟合，确定风廓线稳定度修正参数；根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式计算莫宁-奥布洛夫长度，并根据莫宁-奥布洛夫长度以及稳定度判别标准确定大气稳定度的级别。2.根据权利要求1所述的利用测风塔实测风速数据确定大气稳定度的方法，其特征在于，所述基于预设选取方法对其进行选取得到大气环境中性状态的风速数据的过程，包括：在第一时序集合中，将每个高度层每个时刻的风速数据按照风速大小进行排序，从中选取预设比例的高风速样本数据，得到大气环境中性状态的风速数据。3.根据权利要求1所述的利用测风塔实测风速数据确定大气稳定度的方法，其特征在于，所述预设数据处理包括：求时均值处理和取对数处理。4.根据权利要求1所述的利用测风塔实测风速数据确定大气稳定度的方法，其特征在于，所述大气边界层不同位置的风速与高度对应的关系式为：其中，u为大气边界层不同位置的风速，z为每个高度层的高度，u
*
为摩擦速度，k为冯卡门常数，z0为地表粗糙度，ψ
m
为风廓线稳定度修正参数。5.根据权利要求4所述的利用测风塔实测风速数据确定大气稳定度的方法，其特征在于，所述确定拟合参数的过程，包括：当大气环境处于中性状态时，风廓线稳定度修正参数ψ
m
＝0，则大气边界层不同位置的风速与高度对应的关系式为：将第二时均对数集合带入上述关系式得到拟合参数，所述拟合参数为和lnz0。6.根据权利要求1所述的利用测风塔实测风速数据确定大气稳定度的方法，其特征在于，所述根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式计算莫宁-奥布洛夫长度的过程，包括：根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式，得到莫宁-奥布洛夫长度l的公式为：
其中，z为每个高度层的高度。7.根据权利要求6所述的利用测风塔实测风速数据确定大气稳定度的方法，其特征在于，所述稳定度判别标准为：当-100≤l<0时，大气环境处于强不稳定状态；当-500≤l<-100时，大气环境处于不稳定状态；当|l|>500时，大气环境处于中性状态；当100≤l<500时，大气环境处于稳定状态；当0≤l<100时，大气环境处于强稳定状态。8.一种利用测风塔实测风速数据确定大气稳定度的装置，其特征在于，包括：数据获取单元，用于提取测风塔预设时间范围内预设高度层的实测风速数据，记为第一时序集合，并基于预设选取方法对其进行选取得到大气环境中性状态的风速数据，记为第二时序集合；数据处理单元，用于对所述第一时序集合和第二时序集合分别进行预设数据处理，对应得到第一时均对数集合和第二时均对数集合；参数确定单元，用于根据大气边界层不同位置的风速与高度对应的关系式，对所述第二时均对数集合进行拟合，确定拟合参数；根据所述拟合参数和所述关系式，对所述第一时均对数集合进行拟合，确定风廓线稳定度修正参数；大气稳定度确定单元，用于根据风廓线稳定度修正参数与莫宁-奥布洛夫长度的关系式计算莫宁-奥布洛夫长度，并根据莫宁-奥布洛夫长度以及稳定度判别标准确定大气稳定度的级别。9.一种计算机设备，其特征在于，包括：至少一个处理器，以及与所述至少一个处理器通信连接的存储器，其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器执行权利要求1-7中任一所述的利用测风塔实测风速数据确定大气稳定度的方法。10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使所述计算机执行权利要求1-7中任一所述的利用测风塔实测风速数据确定大气稳定度的方法。

技术总结
本发明公开了利用测风塔实测风速数据确定大气稳定度的方法和装置，方法包括：提取测风塔预设时间范围内预设高度层的实测风速数据，记为第一时序集合，对其基于预设选取方法得到大气环境中性状态的风速数据，记为第二时序集合；对第一和第二时序集合分别进行预设数据处理得到第一和第二时均对数集合；根据大气边界层不同位置的风速与高度对应的关系式对第二时均对数集合进行拟合确定拟合参数，再对第一时均对数集合进行拟合确定风廓线稳定度修正参数；根据莫宁-奥布洛夫长度关系式计算莫宁-奥布洛夫长度，并结合稳定度判别标准用于确定大气稳定度的级别。本发明提供的确定大气稳定度的方法和装置具有计算参数少以及计算简单快捷的优势。算简单快捷的优势。算简单快捷的优势。


技术研发人员：张子良 文仁强 杜梦蛟 易侃 张皓 王浩
受保护的技术使用者：中国长江三峡集团有限公司
技术研发日：2023.03.08
技术公布日：2023/7/19