Об утверждении Положения о повышении точности прогностических оценок радиационных характеристик радиоактивного загрязнения окружающей среды и дозовых нагрузок на персонал и население (утратил силу на основании приказа Ростехнадзора от 08.06.2010 N 465)
Приложение N 4
к Положению о повышении
точности прогностических оценок
радиационных характеристик
радиоактивного загрязнения
окружающей среды
и дозовых нагрузок на
персонал и население,
утвержденному приказом
Федеральной службы
по экологическому, технологическому
и атомному надзору
от 15 января 2010 года N 11
Методика обработки градиентных наблюдений
Существующие методики определения метеорологических параметров довольно сложны, поскольку требуют весьма точных измерений профиля скорости ветра, что возможно только на хорошо оборудованных метеостанциях. Рассматриваемая методика основана на нелинейной модели приземного слоя и является менее требовательной к точности градиентных наблюдений. Вычисления параметров проводят следующим образом [1]. Измеряют на двух уровнях скорость ветра и температуру, используя, как правило, значения z= 2 м и z
= 0,5 м. Находят разность Du = u(z
) - u(z
), D
=
(z
) - (z
).
Используя формулу ([2]) и выражения для скорости ветра и температуры через безразмерные величины (и = 
, = -
, 
), получаем 
/, Du =
.
| (1) |
= 
L
,где u,
- табулированные значения универсальных функций, вычисленных для различных z
(z
= z/L), a Du
, D
- их разность, k - постоянная Кармана. Поскольку Du, D
- измеряемые величины, a Du
, D
зависят от L, то выражение (1) есть неявная функция L. Для нахождения L задаются некоторым значением L
и варьируют его, например, L
= DLi, i = 1, 2, 3, ... N, DL = L
/N до тех пор, пока разность или относительная погрешность
| (2) |
, 
не будет минимальной (в пределе = 0). Найденное значение L*, при котором
минимальна, и определит искомое значение L : L = DL
*. Определив L и пересчитав z
при фиксированных z
и z
, т.е., таким образом, пересчитав D
, Du
найдем v
v | (3) |
или
v | (4) |
При стремлении
0, v
v
. Подобный метод расчета наиболее целесообразен при расчете метеопараметров на ЭВМ. Поскольку параметр L может быть как L > 0, так и L < 0 (при L = 0 режим движения теряет турбулентный характер ([1], стр.74)), то всевозможные вариации L
должны проводиться по формуле: L
= DL(N + L - i), i = 1, 2, 3, ..., N, N + 1, N + 2, ... ,2N + 1. Последнее позволит учесть различную стратификацию слоя атмосферы, задаваемую температурным режимом. Для расчета u
(z
), k
(z
) при найденном L целесообразно пользоваться не таблицами, а аналитическим значением у как функцией z
[2]. Выбор u
,
по заданному z
осуществляется следующим образом при известном z
находят y, по которому из таблицы находят значения соответствующие u
или
. Аналогично находят значения для другого значения z
(другого уровня), вычисляя затем разности
u
,
. После определения параметров L, v
значения u(z), k(z) находят по формулам (9), (10) приложения N 5 к настоящему Положению. Постоянную с
в [2] находят при z = z
и u(z)
= 0.
|
|
Рис.1. Зависимость скорости приземного ветра U(z) от высоты от подстилающей поверхности (модель приземного слоя атмосферы): | Рис.2. Зависимость коэффициента турбулентной диффузии K(z) от высоты от подстилающей поверхности при различных состояниях устойчивости атмосферы: |
Рассчитанные значения u(z), k(z) для двух случаев -L > 0, L < 0 приводятся в виде графиков на рис.1, 2.
Следует отметить, что для применения методики градиентных наблюдений над скоростью ветра и температурой можно использовать не только стандартные, но и любые другие уровни, на которых размещаются датчики для измерения метеопараметров [3, 4]. Более того, в условиях, когда уровень шероховатости подстилающей поверхности нельзя считать однородным на достаточно протяженной области X ~ 1,5 - 2,0 км, для наблюдения метеопараметров целесообразно выбирать уровни размещения датчиков не ниже 20 м [5].