一种核设施废气择时排放的设计方法技术

本发明专利技术涉及一种核设施废气择时排放的设计方法，包括：(1)选取核设施厂址气象塔常规气象观测资料和部分时段的三维湍流资料作为数据基础；(2)建立核设施厂址气象塔常规气象观测资料和三维湍流资料的相似关系，并计算厂址粗糙度；(3)计算厂址逐时的湍流脉动速度，同时计算对应时段的理查森数以及对应的混合长，进而计算相应的拉格朗日时间尺度；(4)将步骤(3)中的逐时数据代入泰勒公式计算逐时扩散参数；(5)根据得到的扩散参数和当前时刻风速计算造成的最大浓度，最大浓度低于设定限值时进行废气排放。本发明专利技术将气象塔常规气象观测数据用于核设施排放时机决策，拓展了核设施气象塔数据的应用价值。的应用价值。

【技术实现步骤摘要】
一种核设施废气择时排放的设计方法


[0001]本专利技术涉及核环境大气扩散参数的估算技术，具体涉及一种核设施废气择时排放的设计方法。

技术介绍

[0002]在核设施的环境评价、设施设计以及日常生产的排放工作中，需要考虑厂址所在区域的大气稀释能力。扩散参数是评价厂址扩散稀释能力的综合指标，扩散参数对核设施排放时机具有重要的指导意义。
[0003]在实际工作中，扩散参数的估算一般采用以下方法：1)大气扩散示踪实验，在厂址释放大气示踪剂，依据采样结果估算主导风向下的大气扩散参数，该方法贴近实际，但实验周期长且费用较高，只能依据部分时段内的扩散情况对大气扩散参数进行估算，无法反映厂址全年不同气象条件下的扩散参数特征。2) 利用三维湍流观测设备，通过观测三维湍流脉动量估算扩散参数，该方法可以对一个完整的时间序列进行观测，因此可以得到完整的各类稳定度下的扩散参数。该方法相对大气扩散示踪实验更易实现，但三维湍流脉动观测设备价格高且易损坏，实际中也只用于数十天或几个月的观测，无法进行长时间的部署。大气边界层湍流脉动量还可采用常规铁塔气象观测资料通过湍流参数化方法计算。考虑到核设施一般配备常规气象观测塔，可以考虑采用气象塔观测资料估算实时扩散参数，进而指导核设施实时排放的系统。

技术实现思路

[0004]本专利技术的目的在于提供一种核设施废气择时排放的设计方法，采用核设施常规气象资料以及部分三维湍流观测资料，通过相似理论建立气象塔常规气象观测资料与三维湍流观测资料的相似性关系，对核设施所在区域的实时扩散参数进行估算，进而计算估算排放中实时的最大落地浓度，与预定限值比较后确定核设施是否进行排放。
[0005]本专利技术的技术方案如下：一种核设施废气择时排放的设计方法，包括如下步骤：
[0006](1)选取核设施厂址气象塔常规气象观测资料和部分时段的三维湍流资料作为数据基础；
[0007](2)建立核设施厂址气象塔常规气象观测资料和三维湍流资料的相似关系，并计算厂址粗糙度；
[0008](3)以建立的相似关系计算厂址逐时的湍流脉动速度，同时计算对应时段的理查森数以及对应的混合长，进而计算相应的拉格朗日时间尺度；
[0009](4)将步骤(3)中的逐时数据代入泰勒公式计算逐时扩散参数；
[0010](5)根据得到的扩散参数和当前时刻风速计算造成的最大浓度，最大浓度低于设定限值时进行废气排放。
[0011]进一步，如上所述的核设施废气择时排放的设计方法，步骤(1)中所述的气象塔常规气象观测资料包括气象塔各层的风速、风向、温度，三维湍流脉动资料包括三维风速、温
度脉动。
[0012]进一步，如上所述的核设施废气择时排放的设计方法，步骤(2)中核设施厂址气象塔常规气象观测资料和三维湍流资料的相似关系建立方法如下：
[0013]根据大气边界层湍流的相似性理论，湍流脉动速度满足如下的相似性关系：
[0014][0015]其中，σ
u,v,w
表示σ
u
或σ
v
或σ
w
，分别为纵向、侧向和垂向湍流脉动速度，U
*
表示特征速度，根据不同情况应取不同参数，例如地面摩擦速度，F
u,v,w
(z/L,h/L) 为纵向、侧向和垂向三个方向对应的无因次函数，z为高度，h为大气边界层的特征尺度，L为莫宁—奥布霍夫长度：
[0016][0017]其中，U
*
为特征速度，k为开尔曼常数，取0.38，θ为位温，w为垂向速度，θ
’
为温度脉动量，w
’
为垂向速度脉动量，
[0018]只考虑核设施近地面层相似性关系函数如下：
[0019][0020]这里U
*
取地面摩擦速度，对中性层结有：
[0021][0022]湍流扩散强弱可以用三个方向的湍流度I
u
、I
v
、I
w
来度量，它们分别为：
[0023][0024]其中，U是平均风速，如果风速廓线也遵守莫宁—奥布霍夫相似性，那么湍流度也遵守相似性，
[0025]通过部分湍流观测数据以及常规气象观测资料，确定相似性关系的函数形式以及中性条件下A、B、C的取值。
[0026]进一步，如上所述的核设施废气择时排放的设计方法，步骤(2)中，假定100米高度以下气层遵守近地面层规律：选取全年中性稳定度的风廓线资料，采用最小二乘法估算地面粗糙度。u2、u1分别为高度z2、 z1对应的风速，z0为地面粗糙度。
[0027]进一步，如上所述的核设施废气择时排放的设计方法，步骤(3)中计算理查森数的方法如下：
[0028]分别表示垂直方向上的风速梯度，表示位温梯度，可以由气象塔常规气象观测资料计算得到，稳定度函数S
h
、S
m
可表示为通量理查森数 R
f
的函数：
[0029][0030][0031]其中，通量查森数R
f
由梯度理查森数R
i
求得：
[0032][0033]其中，C
h
、C
m
、R
f1
、R
f2
为常数，R
fc
表示临界通量理查森数，R
i
≥0.191时， R
ifc
取常数0.2。
[0034]进一步，如上所述的核设施废气择时排放的设计方法，步骤(3)中计算拉格朗日时间尺度的方法如下：
[0035]l表征了影响污染物扩散湍涡尺度的大小，其计算结果为离地面高度z 与稳定度的变量：
[0036][0037]其中，k＝0.38，S
m
为稳定度函数，
[0038]格朗日时间尺度可以表示为，c
t
为常数：
[0039][0040]T
Lu
、T
Lv
、T
Lw
分别为水平u、水平v和w方向的拉格朗日时间尺度，u
’
、v
’
、 w
’
分别表示三个方向的湍流脉动速度。
[0041]进一步，如上所述的核设施废气择时排放的设计方法，步骤(4)中计算扩散参数的方法如下：
[0042][0043]式中，σ
i
为水平或垂向扩散参数，σ
i
为水平或垂向湍流脉动速度，t
L
为拉格朗日时间尺度，T为扩散时间。
[0044]进一步，如上所述的核设施废气择时排放的设计方法，步骤(5)中所述最大浓度的计算方法如下：
[0045][0046]其中，q
m
为最大浓度，Q为释放率，u为风速，H为释放高度，e为自然常数，σ
z
、σ
y
分别为垂向和水平扩散参数。
...

【技术保护点】

【技术特征摘要】
1.一种核设施废气择时排放的设计方法，其特征在于，包括如下步骤：(1)选取核设施厂址气象塔常规气象观测资料和部分时段的三维湍流资料作为数据基础；(2)建立核设施厂址气象塔常规气象观测资料和三维湍流资料的相似关系，并计算厂址粗糙度；(3)以建立的相似关系计算厂址逐时的湍流脉动速度，同时计算对应时段的理查森数以及对应的混合长，进而计算相应的拉格朗日时间尺度；(4)将步骤(3)中的逐时数据代入泰勒公式计算逐时扩散参数；(5)根据得到的扩散参数和当前时刻风速计算造成的最大浓度，最大浓度低于设定限值时进行废气排放。2.如权利要求1所述的核设施废气择时排放的设计方法，其特征在于，步骤(1)中所述的气象塔常规气象观测资料包括气象塔各层的风速、风向、温度，三维湍流脉动资料包括三维风速、温度脉动。3.如权利要求1所述的核设施废气择时排放的设计方法，其特征在于，步骤(2)中核设施厂址气象塔常规气象观测资料和三维湍流资料的相似关系建立方法如下：根据大气边界层湍流的相似性理论，湍流脉动速度满足如下的相似性关系：其中，σ
u,v,w
表示σ
u
或σ
v
或σ
w
，分别为纵向、侧向和垂向湍流脉动速度，U
*
表示特征速度，F
u,v,w
(z/L,h/L)为纵向、侧向和垂向三个方向对应的无因次函数，z为高度，h为大气边界层的特征尺度，L为莫宁—奥布霍夫长度：其中，U
*
为特征速度，k为开尔曼常数，θ为位温，θ
’
为温度脉动量，w
’
为垂向脉动速度，只考虑核设施近地面层相似性关系函数如下：这里U
*
取地面摩擦速度，对中性层结有：湍流扩散强弱可以用三个方向的湍流度I
u
、I
v
、I
w
来度量，它们分别为：其中，U是平均风速，如果风速廓线也遵守莫宁—奥布霍夫相似性，那么湍流度也遵守相似性，通过部分湍流观测数据以及常规气象观测资料，确定相似性关系的函数形式以及中性条件下A、B、C的取值。4.如权利要求1所述的核设施废气择时排放的设计方法，其特征在于，步骤(2)中，假定
100米高度以下气层遵守近地面层规律：其中，u2、u1分别...

【专利技术属性】
技术研发人员：赵多新，张俊芳，廉冰，陈龙泉，姚仁太，黄莎，吕明华，李云鹏，
申请(专利权)人：中国辐射防护研究院，
类型：发明
国别省市：