Для детальных расчетов D определяется по формуле
D |
qHср Vср |
, |
(4) |
|
|||
|
MC |
|
|
где Нср – средняя глубина на рассматриваемом участке, м; Vср – средняя скорость течения водотока на участке, м/с; q – ускорение свободного падения, м/с2;
С – коэффициент Шези, м1/2/c;
М – коэффициент, зависящий от С.
При условии: 10 < С < 60М = 0,7С + 6, при С ≥ 60M = 48 = const.
Произведение МС имеет размерность м/с2.
При переменных гидравлических условиях на отдельных участках распространения сточных вод до расчетного створа коэффициент турбулентной диффузии определяется для каждого участка по выражению (2), а затем для всей расчетной длины по соотношению
D |
D1L1 D2L2 ... DnLn |
, |
(5) |
|
|||
|
L1 L2 ... Ln |
|
|
где D1, D2, …, Dn – коэффициенты турбулентной диффузии отдельных участков,
L1, L2, …, Ln – протяженность отдельных участков. Распространение растворенных и взвешенных веществ в естест-
венных потоках вычисляются уравнениями (6) и (7)
c |
|
|
|
2c |
|
|
2c |
|
|
2c |
|
|
|
c |
|
|
|||||||||
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
U |
|
, |
(6) |
|||||
|
|
|
|
2 |
|
2 |
|
|
2 |
|
|||||||||||||||
t |
|
|
|
x |
|
|
y |
|
|
z |
|
|
|
|
y |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
c |
|
c |
V |
|
|
c |
V |
|
|
c |
V |
|
c |
, |
|
|
(7) |
||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||
t |
|
t |
|
x x |
|
|
y y |
|
z z |
|
|
|
|||||||||||||
где С – концентрация загрязняющего вещества в воде, мг/л; t – время, с;
U – гидравлическая крупность взвешенных веществ, м/с (для растворенных веществ U = 0 и уравнение записывается без последнего члена);
Vx, Vy, Vz — компоненты скорости течения (м/с) относительно коорди-
нат x, y, z (м).
Ось х направлена по течению потока, ось у – от поверхности ко дну, ось z – по ширине потока.
D– коэффициент турбулентной диффузии (м2/с).
Вслучае рассмотрения задачи распространения загрязнения в воде в одной какой-либо плоскости уравнение (6) может быть упрощено.
11
Например, в горизонтальной, тогда уравнение запишется в виде уравнения
V |
c |
D |
2c |
. |
(8) |
|
x |
z2 |
|||||
x |
|
|
|
Эта задача возникает, в частности, в том случае, когда по вертикальной оси перемешивание загрязняющих веществ с водой происходит очень быстро, поперечные течения отсутствуют и интересно проследить распространение загрязняющих веществ по ширине потока на различных расстояниях от места сброса сточных вод.
Используемые при расчете граничные условия основаны на принципе сохранения вещества, учитывающие, что перенос загрязняющего вещества через поверхность, ограничивающую поток жидкости, равен нулю, т. е. для береговой черты потока граничные условия запишутся в виде
D |
c |
0. |
(9) |
|
|||
|
x |
|
|
Поскольку коэффициент D не может равняться нулю (D ≠ 0), то выражение (9) приобретает вид
c |
0. |
(10) |
|
||
z |
|
|
Начальные условия могут задаваться либо в виде распределения концентрации загрязняющего вещества на начальном поперечнике, либо в виде расхода и концентрации поступающего в водный объект загрязняющего вещества с указанием места его поступления.
Уравнение (8) можно записать в форме конечных разностей. Содержащиеся в нем дифференциалы дс, дx, дz заменяются конечными приращениями ∆С, ∆X, ∆Z. Уравнение (8) приобретает вид:
XC |
|
Dср |
|
2Z |
. |
(11) |
|
X |
V |
Z2 |
|||||
|
|
|
|
||||
|
|
ср |
|
|
|
|
Для расчета турбулентной диффузии всю расчетную область потока делим плоскостями, параллельными координатным, на расчетные клетки – элементы. По оси X таких элементов k, по оси Z – m. Каждому элементу присваивается свой индекс по соответствующим осям координат.
Изменение индекса на единицу показывает переход от одного элемента к соседнему. Значениям концентрации в каждом элементе присваиваются те же индексы (рис. 2).
12
|
K |
K+1 |
K+2 |
|
|
|
|
|
|
Экстраполяция |
CK1экстр. |
|
|
|
|
|
|
|
|
Боковая стенка 1 |
CKp1 |
CK,m+1,1 |
|
X |
|
|
|
|
|
2 |
CK, 2 |
|
|
|
|
|
|
|
|
m-2 |
CK,m-2 |
|
|
|
|
|
|
|
|
m-1 |
CK,m-1 |
|
|
|
|
|
|
|
|
m |
CK,m |
CK,m+1,m |
|
|
|
|
|
|
|
m+1 |
CK,m+1 |
|
|
|
|
|
|
|
|
m+2 |
CK,m+2 |
∆X |
∆Z |
|
|
|
|
|
|
|
Z |
|
|
|
Рис 2. Сетка к расчету турбулентной диффузии для условий плоской задачи
Расчетное уравнение, позволяющее определять распространение концентрации загрязняющего вещества по длине и ширине потока, т. е. для условий плоской задачи, записывается в виде уравнения
Ck 1,m 0,5 Ck 1,m Ck,m 1 . |
(12) |
|||
Значения ∆Х и ∆Z связаны зависимостью (13) |
|
|||
V |
Z2 |
|
||
X |
ср |
|
. |
(13) |
|
|
|||
2 D
Коэффициент турбулентной диффузии определяется по формуле
(4).
Когда загрязняющее вещество достигнет граничных поверхностей потока, для расчета диффузии необходимо учитывать особые условия у стенок, уравнение (10), которое в конечных разностях запишется как
|
C |
0. |
(14) |
|
|
|
|
||
|
||||
|
Z |
граничная |
|
|
|
|
|
поверхность |
|
Поле концентраций и расчетную сетку можно условно распространить и за пределы потока (за стенку), т. е. проэкстраполировать концентрацию за ограничивающую водный поток поверхность.
При этом экстраполяционное значение концентрации Ck, экстр в элементе, примыкающем к внешней поверхности стенки, и значение концентрации Ck,1 в элементе, находящемся в потоке и примыкающем к внутренней поверхности стенки, должны удовлетворять условию (14), что возможно только когда соблюдается равенство
13
Ck,экстр Ck,1. |
(15) |
Соотношение (15) определяет правило экстраполяции концентрации загрязняющего вещества. При определении диффузии экстраполяционные значения концентрации используют как действительные. При выполнении расчета на плане водного объекта обозначают место поступления сточных вод в водный объект (начальный створ). Ниже по течению поток схематизируется и делится на расчетные элементы. Скорость поступления сточных вод в водный объект в месте сброса Vст принимается равной скорости течения водотока Vср.
Вычисляется условная площадь поперечного сечения притока δ в месте его впадения по формуле
|
qср |
. |
(16) |
|
|||
Vср |
|
||
Определение ширины загрязненной струи потока b в нулевом створе производят по формуле
b |
|
|
qcn |
. |
(17) |
Hcp |
|
||||
|
|
Hcp Vcp |
|
||
В соответствии с величиной b назначается ширина расчетного элемента ∆Z.
Наиболее допустимая величина ∆Z при береговом сбросе сточных вод находится из соотношения
b |
|
|
qcn |
. |
Hcp |
|
|||
|
|
Hcp Vcp |
||
При выпуске сточных вод на некотором расстоянии от берега или в фарватере величина ∆Z вычисляется из соотношения
Z |
b |
|
qcn |
. |
(18) |
|
|
2 2 Vcp Hcp
Необходимо соблюдать условие, при котором при назначении величины ∆Z выполнялось бы неравенство
Z |
1 |
B, |
(19) |
|
|||
10 |
|
|
|
где B – ширина водотока.
Таким образом, при расчете турбулентной диффузии весь участок потока (от нулевого створа до расчетного или створа, который нас интересует по условиям решаемой задачи) делят на клетки со сторонами ∆Х и ∆Z, получая расчетную сетку.
Клетки, попавшие в водоток со сточными водами в начальном
14
поперечнике (нулевой створ), заполняются числами, выражающими концентрацию загрязняющего вещества в сточной воде, остальные клетки заполняются числами, отражающими естественную концентрацию загрязняющего вещества в водотоке (в частном случае это может быть нулевая концентрация).
Если протяженность интересующего участка водотока велика, а размеры клеток малы, то расчет ведут до определенного створа, после чего клетки в сечении объединяют (укрупняют), получая новые средние значения концентрации загрязняющего вещества и новые линейные параметры клетки. Значения концентраций загрязняющего вещества получают как среднее арифметическое из суммы концентраций объединенных клеток.
Укрупнение клеток можно повторять несколько раз, начиная с определенного раствора.
Пример расчета
Определить максимальную концентрацию загрязняющего вещества в водотоке на расстоянии 700 м от места выпуска сточных вод по схеме плоской задачи. Выпуск сточных вод – береговой.
Расход сточных вод qст = 50,6 м3/с. Водоток характеризуется следующими показателями:
– |
средняя скорость течения |
Vср = 2,42 |
м/с; |
– |
средняя глубина водотока |
Нср = 2,37 |
м; |
– |
ширина водотока |
В = 26,5 м. |
|
Коэффициент турбулентной диффузии D = 0,073 м2/с. Для упрощения расчетов примем, что фоновое загрязнение водотока отсутству-
ет, т. е. Св = 0, а концентрация загрязняющего вещества в сточной воде
Сст = 100 г/м3.
Решение. Определим начальное сечение струи:
qср 50,6 20,9м2. Vср 2,42
Определим ширину загрязненной части водотока:
b |
|
|
20,9 |
8,8м. |
|
|
Hcp 2,37
Выбираем ширину расчетной клетки, соблюдая условие
Z 1,3м.
В этом случае число клеток по ширине потока, занятых загрязненной водой в результате выпуска сточных вод:
15