
книги из ГПНТБ / Бызова, Н. Л. Рассеяние примеси в пограничном слое атмосферы
.pdfи при малых временах диффузии отсутствует. Таким образом, сле дует ожидать его проявления на тех стадиях процесса, когда вре мя диффузии уже не мало по отношению к лагранжеву масштабу; следовательно, этот эффект должен сказываться на вертикальной диффузии заметнее, чем на горизонтальной.
Параметром, который определяет степень влияния гравитаци
онного оседания, является величина v', определяемая |
выражением |
|
(2.29). Связь параметра у' с параметром |
v, который |
используется |
в полуэмпирической теории, определяется |
выражением |
|
v' = 2p£v. |
|
(5.38) |
Поскольку В при уменьшении неустойчивости и переходе к устой чивым состояниям уменьшается, а $ растет обратно пропорцио
нально |
интенсивности той или иной пульсационной компоненты, |
||
можно ожидать, что величина ВВ при изменении |
характера устой |
||
чивости |
меняется мало. Для состояний, близких |
к безразличному, |
|
с отклонением в сторону слабой неустойчивости |
в приземном |
слое |
|
атмосферы В составляет около 5, при этом 28В приближенно |
мож |
||
но положить равным 0,2. |
|
|
Из-за большого разнообразия метеорологических условий экс периментальное выявление гравитационного оседания возможно только при синхронном проведении опытов, когда все условия (вы
сота источника, время его действия) |
совпадают |
и |
различаются |
только скорости гравитационного оседания. Таких |
парных опытов |
||
в рассматриваемом массиве имеется |
20; из них в |
12 |
случаях вы |
сота источника была более или равна 25 м, синхронность соблю далась неполностью (в части случаев выпуски частиц разной крупности проводились не одновременно, а последовательно), вре
мя действия |
источника |
немного |
различалось. |
Результаты этих |
||
12 пар не позволили |
выявить |
никакой тенденции |
зависимости |
|||
параметров диффузии Ьу |
и В от параметра v. |
более |
тщательно. |
|||
Остальные |
восемь |
опытов были поставлены |
Здесь высота источника составляла 2 м, что позволило использо
вать достаточно |
развитую пробоотборную сеть, |
опыты |
проводи |
|||
лись - синхронно. Метеорологические условия |
в |
четырех |
|
случаях |
||
соответствовали |
безразличной |
стратификации |
(класс 4), |
в |
осталь |
|
ных случаях — неустойчивой |
(классы 3 и 2). |
|
|
|
|
|
Результаты всех парных опытов показаны |
на |
рис. 5.12 |
и 5.13, |
где даны отношения характеристик диффузии В и Ъу к их соответ ствующим средним для данного класса устойчивости и высоты значениям в зависимости от параметра v. Из рис. 5.12 следует, что в пределах использованных времен диффузии и параметра v вли яние гравитационного оседания на горизонтальную диффузию практически отсутствует. Параметр же вертикальной диффузии В, который линейно связан с коэффициентом турбулентной диффу зии, обнаруживает явную зависимость от параметра v.
На рис. 5.13 приведена |
также |
шкала v' в предположении, что |
265 = 0,20 и экспериментальные |
точки, полученные Дунским и |
|
Оникулом [1968] по опытам |
Дунского и соавторов {il966]; нанесена |
также теоретическая зависимость, полученная в гл. 2. Эксперимен тальные данные позволяют уловить влияние гравитационного осе-
|
х |
• |
|
X |
|
|
|
|
|
to |
|
|
|
|
|
|
|
|
|
| х х х |
• |
• |
• |
|
|
|
|
|
|
|
£*х х х |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
* 2 |
|
Рис. |
5.12. |
Зависимость |
значений 6 у , |
|
|||||
нормированных |
на |
средние |
значения |
|
|||||
для |
данной |
устойчивости, |
от |
парамет |
|
||||
|
|
|
ра |
v—w/BU: |
|
|
|
|
|
/ — специальные парные |
опыты, |
Я= 2 м; 2 — |
|
||||||
остальные |
парные |
опыты, |
Н>25 м |
|
|||||
Б |
|
|
|
|
|
|
|
|
|
х |
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
• |
1 |
|
|
|
|
|
|
|
|
||
1.0 |
|
|
|
|
|
|
|
х |
г |
|
|
|
|
X |
|
|
|
о |
J |
Х |
X |
о . |
|
|
|
|
|
|
|
О • |
1 |
!>' |
РИС. 5.13. Зависимость значений В, нормирован
ных .на средние значения для данной устойчи- *
•вости, от параметра v=w BU и v'=*$By:
I — специальные |
парные |
опыты, И =2 м; 2 — остальные |
парные опыты, |
Я?-25 |
м; 3 — результаты Дунского |
(1958); |
кривая |
рассчитана по (2.26) |
дания на вертикальное рассеяние примеси, и оно оказывается близким к расчетному.
ЗАКЛЮЧЕНИЕ
Рассеяние примеси в пограничном слое атмосферы определяет ся свойствами самой примеси и источника, состоянием турбулент ности атмосферы и особенностями среднего потока, характером подстилающей поверхности и, наконец, взаимодействием приме си и среды.
Задачи расчетов рассеяния примеси от локальных источников можно разделить на три типа:
1. Расчеты концентрации или загрязнения при стационарных или квазистацпонарных условиях, когда предметом расчета яв ляются параметры загрязненной зоны, осредненные за время от нескольких минут до нескольких часов.
2. Те же параметры находятся при -осреднениях за большие промежутки времени — сутки, месяцы, годы; определяются экс тремальные условия, способствующие наиболее благоприятным или
неблагоприятным с точки зрения практики |
результатам. |
3. Требуется определять мгновенные |
значения концентраций |
или поведение клубов примеси от источников мгновенного дей ствия.
Не затрагивая свойства примеси и источников, которые в на стоящей работе практически не рассматривались (кроме гравита ционного оседания частиц), отметим те области, в которых наши
знания о свойствах атмосферы, определяющих |
рассеяние при |
меси, достаточны с точки зрения практических |
применений. |
Наилучшее состояние здесь соответствует |
задачам первого |
типа. Можно считать, что |
имеющиеся модели диффузии обеспечи |
||
вают удовлетворительную |
для |
практических |
применений точность, |
а установленные связи между |
различными |
метеорологическими |
параметрами и параметрами диффузии для нижней части погра ничного слоя атмосферы достаточны для выбора этих последних на основе измерения сравнительно несложных метеорологических характеристик. Влияние подстилающей поверхности, если только она не представляет собой сложный рельеф или город, также под дается учету.
Для верхней части пограничного слоя атмосферы полученные связи существенно менее надежны. Кроме того, на этих высотах может сказываться вертикальная составляющая атмосферных дви жений синоптического происхождения. Недостаточно также дан ных о взаимодействии очень мелких частиц с подстилающей по верхностью, а также сведений о поведении примеси в условиях
сложного рельефа и городской застройки. Мало данных о мезомасштабных закономерностях, определяющих рассеяние примеси на большие расстояния. Определенная осторожность требуется при расчетах диффузии во время нестационарных метеорологических условий.
Для задач третьего типа, т. е, расчета концентраций и рассея ния клубов примеси от действия мгновенного источника к недос таточно известным факторам надо добавить также влияние сдвига ветра и учет конечности скорости диффузии.
В зависимости от конкретных целей тех или иных оценок и расчетов рассеяния примеси в атмосфере необходимо определять требуемую степень точности, надежности и строгости подхода. В ряде случаев с помощью более тщательного подбора модели и обеспечения большей информации о состоянии атмосферы можно достигнуть значительного улучшения результатов, но в других слу чаях это принципиально невозможно из-за естественного разнооб разия и случайного характера, присущих факторам, определяю щим рассеяние примеси в атмосфере (пространственная и вре менная изменчивость метеорологических условий, неоднородность подстилающей поверхности).
Что касается расчетов второго типа то, как нам кажется, для их проведения, по-видимому, нет принципиальных затруднений, но требуется достаточная климатологическая и синоптическая инфор
мация |
и |
разработка |
методов |
ее |
применения. |
В |
этом |
направле |
|||||||||||||||||||
нии |
многое |
уже сделано |
как за рубежом, |
|
так |
и |
в |
Советском |
|||||||||||||||||||
Союзе. |
|
|
|
|
|
|
|
|
CONCLUSION |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
The |
diffusion |
in |
(he |
surface |
boundary |
layer |
depends |
on |
(he |
source and |
|||||||||||||||||
substance |
|
peculiarities, |
the |
turbulence |
state |
of |
atmosphere, |
|
the |
underlying |
|||||||||||||||||
surface and at last the substance |
interaction |
with |
the |
environment. |
|
|
|
||||||||||||||||||||
The |
problems |
of |
calculating |
the |
local |
|
source |
diffusion |
|
may |
be |
of |
three |
||||||||||||||
types: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1. Concentration or pollution calculations in stationary and quasi-stationary |
|||||||||||||||||||||||||||
conditions, when the contaminated area parameters |
averaged |
for |
over |
a |
period |
||||||||||||||||||||||
from several minutes to several hours, are needed. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
2. The |
same |
parameters |
are |
calculated |
for |
large time |
|
intervals: |
days, |
||||||||||||||||||
months, |
years; |
extremum |
conditions, |
most |
favourable |
and |
unfavourable |
From |
|||||||||||||||||||
the practical point of view, are found. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
3. |
The |
instantaneous |
puff |
concentration |
|
values |
or |
the |
behaviour |
estimates |
|||||||||||||||||
of puffs from the instataneous sources are required. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Leaving |
out |
the |
substance |
and |
source |
characteristics |
not |
|
considered |
in |
|||||||||||||||||
our book |
|
at |
all |
(besides |
the |
gravitational |
|
deposition |
of particles) |
we |
shall |
||||||||||||||||
mention |
those |
areas |
where |
our |
knowledge |
of the |
|
atmospheric |
diffusion |
is |
|||||||||||||||||
sufficient |
to be applied in practice. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
The best state in this case will correspond to the |
first |
|
type |
problems. One |
|||||||||||||||||||||||
may |
assume |
|
that the diffusion |
models |
available |
are |
accurate |
enough |
for |
||||||||||||||||||
practical |
applications, |
and the |
established |
correlations |
between |
meteorological |
|||||||||||||||||||||
and diffusion parameters Tor lower boundary |
layer |
are |
sufficient |
for |
estimating |
||||||||||||||||||||||
the latter on the basis of comparatively simple meteorological |
parameters. The |
||||||||||||||||||||||||||
underlying surface effect may be also taken into account, |
if |
the |
underlying |
sur |
|||||||||||||||||||||||
face is neither |
a complex relief nor a town |
building system. |
|
|
|
|
|
|
|
|
|
||||||||||||||||
For the upper boundary layer the correlations obtained |
are |
substantially |
less |
||||||||||||||||||||||||
reliable. |
Moreover, the vertical |
component |
of synoptical |
origin |
can |
count |
at |
||||||||||||||||||||
these |
altitudes. |
The |
scarecity |
|
of |
data |
on |
small |
particle |
interaction |
with |
the |
underlying |
surface |
and |
on the |
dispersion |
in |
complex |
relief |
and |
town |
|||||||
building |
conditions is |
also felt. There are very |
little data on rhesoscale cliffuson |
|||||||||||||
as well. |
One |
should |
be |
careful |
in calculating |
diffusion during |
nonstationary |
|||||||||
meteorological |
conditions. |
|
|
|
|
|
|
|
|
|
|
|
||||
In case |
|
of |
the |
third |
|
type problem, i. e. |
puff |
concentrations and |
their |
|||||||
dispersion, |
the |
wind |
shear |
effect |
and |
finite diffusion |
rate |
are |
waiting |
their |
||||||
further |
investigation. |
|
|
|
|
|
|
|
|
|
|
|
|
|||
One |
should |
determine |
the |
desired |
degree |
of |
accuracy, reliability and |
|||||||||
rigidness |
of |
approach as |
a |
function of particular uses |
of |
diffusion |
estimates. In |
a number ef cases it is possible to obtain better data with the help of a more thorough model choice and wider information on atmosphere; nonetheless this is impossible in some other cases due to the natural variety and randomness of
factors |
determining |
diffusion |
(space |
and |
time variations |
of |
meteorological |
|||
conditions, nonuniformity of |
the underlying surface). |
|
|
|
|
|||||
As |
for the second |
type |
calculations |
these |
might |
seem |
possible, |
no |
||
principal |
difficulties |
being encountered, |
though |
sufficient |
elimatological |
end |
||||
synoptical data and well-elaborated technique of their application |
are required. |
|||||||||
Much has been done in this field both abroad and in our country. |
|
|
ПРИЛОЖЕНИЕ
ПРИЛОЖЕНИЕ
ФОРМУЛЫ И ПАРАМЕТРЫ ДЛЯ РАСЧЕТА ХАРАКТЕРИСТИК ПРИЗЕМНОЙ КОНЦЕНТРАЦИИ И ВЫПАДЕНИИ
Формулы и параметры для расчета характеристик приземной концентрации и плотности выпадений составлены на основе изло женных в книге моделей и экспериментальных результатов. Они дополняют аналогичные материалы, основанные на применении других моделей и приведенные в [МАЭ, 1955; МАЭ, 1968].
Обозначения те же, что в главах 2 и 5 (см. список обозначе ний) .
А. Типизация метеорологических |
условий (см. |
гл. 3) |
|
|
• |
Состояние нижнего слоя атмосферы |
определяется |
классом ус |
тойчивости. При обычных условиях по классу устойчивости непо средственно определяются параметры диффузии. При особых усло виях точность определения падает, иногда требуются дополни
тельные |
поправки. |
|
|
|
|
|
||
К особым условиям относятся: |
|
|
|
|
||||
а) штили и слабые ветры; |
|
|
|
|
до Н (по |
|||
б) аномальные профили ветра в пределах от земли |
||||||||
вороты на 90° и более, резко |
выраженный |
максимум скорости); |
||||||
в) |
приподнятые инверсии; |
|
|
|
|
|||
г) |
переходное время суток |
(за |
1—2 |
ч до захода солнца, 1—3 ч |
||||
после восхода при ясной или почти ясной погоде); |
|
|||||||
д) |
осадки, |
туманы. |
|
|
|
|
|
|
Принято следующее соответствие между номером класса и ха |
||||||||
рактером |
устойчивости: |
|
|
|
|
|
||
сильная неустойчивость |
— |
класс 1; |
|
|
||||
умеренная |
» |
— |
» |
2; |
|
|
||
слабая |
» |
— |
» |
3; |
|
|
||
нейтральные условия |
— |
» |
4; |
|
|
|||
слабая устойчивость |
|
» |
5; |
|
|
|||
умеренная |
» |
|
» |
6; |
|
|
||
сильная |
» |
|
» |
7. |
|
|
||
Класс устойчивости можно определить различными способами. |
||||||||
а." По |
градиентным измерениям в |
приземном слое |
атмосферы. |
|||||
Соответствие |
с параметром |
устойчивости |
Б (формула |
3.9) — см. |
||||
табл. |
3.8 |
и 3.21. |
|
|
|
|
|
б. По градиентным измерениям в слое от 0 до Н. Класс |
устой |
||||||
чивости определяется по табл. 3.22 и 3.23. |
|
|
|
|
|||
в. По метеорологическим измерениям сетевого типа. Класс ус |
|||||||
тойчивости |
определяется |
по табл. 3.1—3.3 |
или 3.14—3.16 |
(более |
|||
точно). |
|
|
|
|
|
|
I |
Способ |
определения |
устойчивости выбирается |
в зависимости |
||||
от Н и требуемой |
точности. |
|
|
|
|
||
Б. Формулы |
для |
расчета характеристик |
приземной |
|
|||
|
|
|
концентрации |
|
|
|
|
Для постоянно |
действующего источника |
q — |
концентрация, |
||||
р — вертикальный |
поток |
на уровне земли, |
Q |
— |
производитель |
||
ность источника, Т — время осреднения при измерении. Для ис |
точника короткого времени действия q — интеграл от концентра |
||||||||
ции, р — плотность выпадений, |
Q |
— общее |
количество выпущен |
|||||
ного вещества, Т — время действия |
источника. |
|
||||||
|
Коэффициенты и параметры в формулах |
безразмерны, так что |
||||||
можно пользоваться любой системой единиц. |
|
|
||||||
|
Легкая (газообразная) |
примесь, не взаимодействующая с по |
||||||
верхностью земли (гл. 2, разд. 2.2.1). |
|
|
|
|||||
|
Наземная концентрация от точечного источника |
(формулы 2.1, |
||||||
2.37—2.40): |
|
|
|
|
|
|
|
|
|
|
Q ехр |
X |
|
ч |
|
|
|
|
q(x,y) |
|
|
|
(Б.1) |
|||
|
|
Г(1+р)1/2тсо-у £/ |
|
|
||||
где |
|
|
|
СН |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
В |
|
|
|
|
|
То же для линейного нормального к ветру |
источника |
||||||
|
|
Q ехр |
R |
|
(К |
|
|
|
|
|
|
|
|
|
|||
|
q{x) = |
|
|
|
|
|
(Б.2) |
|
|
Г(1 + р) UH |
\ х |
|
|||||
|
|
|
|
|||||
|
Расстояние до зоны максимальной наземной |
концентрации |
||||||
(формулы 2.49) |
|
|
|
|
|
|
|
|
|
|
хй |
— Лг |
Н_ |
|
|
(Б.З) |
|
|
|
|
|
В |
|
|
|
|
(г'=1 для линейного источника, 1 = 2 для точечного). |
|
|||||||
|
Максимальная наземная концентрация для линейного источни |
|||||||
ка |
(2.49 и 2.50) |
|
|
|
|
|
|
|
|
|
с7о = |
А т |
~ |
, |
|
(Б.4) |
|
|
|
|
UH |
|
|
|
|
|
для |
точечного |
|
QB |
|
|
|
|
|
|
|
|
|
|
|
(Б.5) |
||
|
|
|
|
|
|
|
|
urn*
Поперечное рассеяние (для легкой и оседающей примеси оди наково — гл. 5, разд. 5.3)
oy = by(H, Т)х, |
(Б.6) |
где |
|
/ Т |
\0'2 |
|
Ьу{Н, П = * у о ( 4 ) |
- |
(Б.7) |
|||
где Т — в минутах. |
|
|
|
|
|
|
Или |
оу=ауха. |
|
|
|
|
(Б.8) |
|
|
|
|
|
||
Оседающая примесь (гл. 2, разд. 2.2, гл. 5, разд. 5.1). |
||||||
Вспомогательные параметры |
(формулы 2.59, 2,63, 2.66) |
|||||
|
"> = -z—, |
xw |
= |
|
, |
(Б.9) |
|
BU |
|
|
w |
|
|
где |
|
|
|
|
|
|
|
5 = 2(1 + m)B** |
2В. |
|
(Б.10) |
||
Максимальная приземная концентрация, точечный источник |
||||||
Чо = 1 ^ - Ф ( У ) |
= ^ч<У), |
(Б. 11) |
||||
линейный источник • |
|
|
|
|
|
|
|
?о = - ^ - Ф ( > ) - |
|
(Б.12) |
|||
|
UN |
|
|
|
|
|
Максимальная |
плотность |
осадка |
(формулы |
2.61, 1.64—1.67, |
||
табл. 2.3) |
|
|
|
|
|
|
|
pQ = w q0 при |
v > 0 , l , |
(Б.13) |
|||
|
|
|||||
или (точечный источник, формула 5.36) |
|
|
||||
|
" ' - - Л - ^ Ш ' |
• |
< Б Л 4 : |
|||
|
при 0,005 < — < 0 , 2 5 . |
|
||||
|
|
|
U |
|
|
|
Расстояние до точки максимальной плотности осадка или кон |
||||||
центрации вдоль оси х |
|
|
|
|
|
|
Х О = |
J ^ L _ = XTfM |
|
= |
X E / 2 ( V ) , |
(Б-15 |
|
|
xw + xT |
|
|
|
|
,\ •.! |
где для точечного |
источника |
|
|
' |
|
••• |
|
j c r = - ^ - , |
|
|
(Б. 16 |
||
|
|
2В |
|
|
|
для линейного |
|
Н |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
(Б.17) |
|
|
|
|
хт |
———. |
|
|
|
||
|
|
|
|
|
|
|
|||
|
|
|
|
В |
|
|
|
|
|
Длина |
зоны |
загрязнения |
координата |
вдоль |
оси х, где |
||||
<7 —0,1 а0 |
(формула 5.9), |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
(Б.18) |
Значения А,г при б ф 0,1 определяются с помощью |
уравнения |
||||||||
(2.68) |
и рис. 5.1. |
|
|
|
|
|
|
|
|
Доля осевшей примеси до расстояния пх0 |
|
|
|
||||||
|
|
|
1 ГГр{х, |
у) |
dxdy |
|
|
(Б.19) |
|
|
|
|
Q |
|
|
|
|
|
|
определяется с помощью рис. 5.2. |
|
|
|
|
|
||||
Доля осевшей примеси внутри площади |
|
|
|
||||||
|
|
S = 4,28 оу (л'0 ) (л - А-0) = |
(v) |
|
|
(Б.20) |
|||
|
|
|
|
|
|
В |
|
|
|
составляет г{2 и определяется с помощью табл. В.З. |
|
|
|||||||
В. |
Определение |
безразмерных |
параметров |
для |
расчетных |
||||
|
|
|
|
формул |
|
|
|
|
|
Параметр вертикальной |
диффузии |
в зависимости |
от И и клас |
са устойчивости — табл. В.1, а также табл. 5.5.—5.7. Поправки на
слабый ветер — табл. 5.23 и 5.24; |
для приподнятых |
инверсий — |
||||||
табл. 5.8. |
|
|
|
|
|
|
Таблица В.1 |
|
|
|
|
|
|
|
|
||
Параметр |
В в зависимости от Я и класса |
устойчивости |
||||||
|
|
|
|
Класс |
устойчивости |
|
|
|
Н ы |
1 |
|
2 |
3 |
|
4 |
5 |
6 |
|
|
|
||||||
2 |
0,060 |
0,030 |
0,022 |
0,015 |
0,006 |
0,003 |
||
25—70 |
0,050 |
0,030 |
0,020 |
0,012 |
0,005 |
0,002 |
||
100-200 |
0,050 |
0.030 |
0,019 |
0,011 |
0,003 |
0,0015 |
||
200—300 |
0,050 |
0,028 |
0,015 |
0,010 |
0,003 |
0,0010 |
||
500 |
0,050 |
0,010 |
0,010 |
0,08* |
0,003* |
0,001* |
||
1000 |
0,005 |
0,005* |
0,005 |
0,004* |
|
|
||
* Менее надежные данные. |
|
|
|
|
||||
Параметры |
поперечной диффузии: |
Ьу0 — табл. |
5.17 и 5.18; |
|||||
'У и а — табл. 5.13—5.16. |
р, С, Dh |
D2, Ai, А2 |
для формул |
|||||
Безразмерные |
параметры |
|||||||
Б.1) — (Б.5) —табл. В.2. |
|
|
|
|
|
|||
Безразмерные функции v для формул |
(5.11) — (5.15) —рис. 5.3 |
|||||||
точечный источник) |
и табл. В.З, В.4. |
|
|
|