Электронная библиотека » Вадим Романов » » онлайн чтение - страница 10


  • Текст добавлен: 6 июля 2014, 11:37


Автор книги: Вадим Романов


Жанр: Прочая образовательная литература, Наука и Образование


Возрастные ограничения: +12

сообщить о неприемлемом содержимом

Текущая страница: 10 (всего у книги 17 страниц)

Шрифт:
- 100% +

Оценки, проведенные в работе [92], показывают, что отклонение собственной скорость вихрей от скорости потока не превышает нескольких процентов, и ее влиянием на процессы рассеивания загрязняющих примесей можно пренебречь.

3.7. Измерения геометрических и динамических характеристик выбросов

Исследование физических процессов при авариях и их последствий в различных средах, как и во многих других областях науки и практики, приводит к необходимости построения моделей различных процессов и явлений, отражающих реальность. Эти задачи настолько сложны, что для успешного их решения надо умело сочетать теоретические представления в данной области знаний с использованием экспериментальных или статистических данных, относящихся к конкретному явлению.

В области протекания процессов горения детонации и взрыва в настоящее время глубоко разработаны теоретические основы этих явлений [60, 61, 77, 80, 82, 103], заранее известна структура модели и основные зависимости. Экспериментальные данные могут служить здесь лишь для определения отдельных уточняющих параметров модели.

Совершенно иная картина в области изучения формирования и развития кратковременных выбросов в реальной атмосфере. В этой области теоретические знания далеки от точных количественных представлений, а немногочисленные эксперименты [13, 33, 48-52, 133], проведенные в узких диапазонах изменений определяющих параметров, как правило, не дают общей картины явления. Следует отметить, что и в тех областях, где сравнительно хорошо разработана количественная теория (например, подъем термина в стратифицированной атмосфере), исследователь часто сталкивается с задачами, в которых по экспериментальным или статистическим данным требуется не просто определить отдельные параметры модели, но и в существенной мере восстановить общую картину явления, которая заранее может быть ясна лишь в очень грубом приближении или совсем неясна.

К таким задачам относятся, в нашем случае, задачи интерпретации геометрических измерений, когда по данным дальномерных съемок взрывных выбросов требуется определить картину формирования первичного газопылевого выброса, а также построить модель разлета твердой фазы взрыва из источника и подъема перегретого газовоздушного объема, насыщенного твердой фазой различного состава и дисперсности.

Сложная задача построения модели явления в целом, как правило, может быть разложена на этапы и фазы, на каждой из которых с успехом могут быть использованы формализованные методы обработки данных. Такими задачами являются задачи восстановления зависимостей. При решении задачи восстановления неизвестной зависимости по эмпирическим данным первым желанием бывает [84] искать зависимость как можно более общего вида, привлекая как можно большее число аргументов. Однако такой путь неизбежно сталкивается с ограниченностью экспериментального или статистического материала, которым располагает исследователь и ограниченностью вычислительных возможностей. Это противоречие разрешается при использовании структурной минимизации задачи. Она состоит в поиске вначале предельно грубой модели, а затем эта модель постепенно усложняется до достижения оптимального соотношения между точностью аппроксимации эмпирического материала и надежностью результата в условиях ограниченного объема данных.

Важным этапом построения математических моделей физического явления является его верификация, использующая данные экспериментов.

Измерения выбросов загрязняющих веществ в реальной атмосфере, кроме организационных и финансовых затрат, имеет ряд трудностей, связанных с методическими вопросами. Сравнение результатов измерений геометрических характеристик формирующихся в атмосфере объемов, полученных разными методами, является не всегда корректным и может привести к большим различиям.

В первую очередь это связано с тем, что при проведении натурных опытов с использованием измерительных устройств теодолитного типа определяются горизонтальные и вертикальные углы при наведении на граничные контуры выбросов. При этом могут возникнуть большие погрешности, связанные с частичным затенением одних частей объема другими.

Рисунок 3.11 иллюстрирует и частично объясняет противоречивость и нестыковку многих литературных данных о размерах формирующихся в атмосфере выбросов от эквивалентных источников. Кроме неизвестных различий в турбулентной активности воздуха при различных замерах геометрии выбросов существенное значение имеет и фактор места регистрации объекта. Он проявляется в форме зависимости полученных результатов от удаления измерительной аппаратуры от объекта и от его абсолютных размеров.

Из Рис. 3.11 видно, что, например, высотные размеры выбросов измеряемые в точках А и В, дают существенно разные результаты ZAА2 – НА1 и ZB– НВ2 – Нт , причем ZА ≠ Zв ≠ Z ,

где ZА и Zв – высоты выбросов при замерах в соответствующих точках; Z – истинный высотный его размер.

Такая же ситуация возникает при измерениях горизонтального размера выброса X.

Наряду с теодолитными замерами в некоторых натурных экспериментах проводятся самолетные и вертолетные зондировки, фотографирования дымовых выбросов и стереофотограмметрическая съемка [157].

Отметим, что измерения, производимые с летательных аппаратов, также не лишены недостатков. Самые существенные из них – влияние воздушных потоков на выброс и временной сдвиг в процессе измерения.



Рис. 3.11. Схема зависимостей визуальных размеров выброса от места измерений на поверхности земли: 1 – выброс; А, В – точки размещения измерительной аппаратуры; X, Z – истинные продольный и высотный размеры выброса.


Стереофотограмметрия весьма трудоемка и дорога, а ее результаты часто не оправдывают возложенных надежд.

Ошибки в определении характеристик выбросов

Основными натурными измерениями характеристик атмосферных выбросов, имеющими научный интерес в части изучения их динамики, являются фото– и киносъемки, а также регистрации геометрических размеров выбросов с помощью дальномера или высотомера. Такие измерения обычно проводятся с одновременным измерением метеорологических характеристик приземного слоя атмосферного воздуха, и после статистической обработки они могут иметь представительный характер.

Под измерением понимается [155] сравнение измеряемой величины с другой величиной, принятой за единицу измерения. При кинофоторегистрации выбросов и их измерении дальномером или высотомером определяемыми характеристиками выбросов являются угловые значения и удаления верхней и нижней, а также боковых границ наблюдаемого объема. Эти характеристики могут быть сравнены с единицей измерения длины и угла при помощи измерительных приборов, проградуированных в соответствующих единицах. Такие измерения изучаемых объектов могут отсчитываться по шкале прибора. Они называются прямыми.

При косвенных измерениях измеряемая величина вычисляется при использовании результатов прямых измерений других величин, которые связаны с измеряемой величиной определенной функциональной зависимостью. Примерами косвенных измерений при проведении работ по изучению атмосферных выбросов могут служить измерения высот верхней и нижней границ объема, его поперечного размера, координат центра масс в системе координат, связанной с точкой взрыва или пожара, используя линейные и угловые координаты соответствующих частей взрывного выброса. Динамические характеристики изучаемого объема определяются при использовании измерений длин пройденного пути и соответствующих промежутков времени.

Известно, что при измерении любой величины никогда не получают истинного значения этой величины – результат измерения дает лишь приближенное значение. Это объясняется как принципиально ограниченной возможностью точности измерения, так и природой самих измеряемых объектов.

Погрешности результата измерений определяются разностью измеряемой и истинной величин и зависят от многих причин. Погрешности (ошибки) разделяются на систематические и случайные. Систематические погрешности связаны с ограниченной точностью изготовления прибора (погрешностью прибора), неправильным выбором метода измерения, неправильной установкой прибора. Они также появляются, если пренебречь действием некоторых внешних факторов.

Таким фактором может быть неучет сноса ветром облака загрязняющих веществ при определении его координат в горизонтальной плоскости и пренебрежение высотным разворотом ветра при определении горизонтального размера облака. Неправильная установка начальных углов дальномера в горизонтальной и вертикальной плоскостях приводит к появлению систематических ошибок в измерениях линейных и угловых размеров выброса.

Таким образом, систематические погрешности вызываются вполне определенными причинами, их величина при всех повторных измерениях либо остается постоянной (как в случаях округления или смещения нуля шкалы прибора и т.п.), либо изменяется по определенному закону (в случае изменения ветра с высотой). Так как причины, вызывающие систематические погрешности, в большинстве случаев известны, то эти погрешности, в принципе, могут быть исключены изменением метода измерений, введением поправок к показаниям приборов, а также учетом систематического влияния внешних факторов.

Случайные погрешности вызываются большим числом случайных причин, действие которых на каждое измерение различно и не может быть заранее учтено. Сюда могут быть отнесено различие механических и термодинамических свойств взрываемых или сгорающих изделий, наличие флуктуации метеорологических параметров за время проведения серии экспериментов, различие состава и прочностных характеристик грунта в месте проведения работ.

Кроме того, ошибки случайного характера обуславливаются качеством и чувствительностью фотоматериалов и используемых кино– фотокамер. Появление случайных ошибок связано с техникой запечатления и увеличения изображения на фотопленке и фотобумаге. На точность определения пространственных координат точек облака загрязнений влияют ошибки измерения снимков и ошибки определения элементов внутреннего и внешнего ориентирования снимков [156]. В настоящее время ошибки в снятии плоских координат X, Z соизмеримы с разрешающей способностью фотоматериалов [157], поэтому дальнейшее повышение точности измерений снимков будет бесполезным, если не повысить разрешающую способность фотоматериалов и не снизить их деформацию.

Для того, чтобы правильно истолковать и использовать полученные в экспериментах результаты, необходимо знание достоверности проведенных расчетов, выявить ошибки, которые могут повлиять на точность вычислений, оценить погрешность полученного результата.

Источниками погрешностей и сшибок при вычислении могут быть недостаточно точное отображение реальных явлений их математической моделью (погрешности модели), приближенное знание величин, входящих в условие задачи (погрешности исходных данных), применение приближенных расчетных формул вместо точных (методические погрешности), особенности работы операционного устройства вычислительной машины (операционные погрешности). Дадим краткое описание этих погрешностей.

Погрешности модели.

Известно, что совокупность факторов реального физического явления невозможно в полной мере отразить в его математической модели. Чем больше факторов учитывается, тем сложнее анализ явления и расчет по его математической модели. Поэтому целесообразно принять различные условия и допущения, упрощающие решение задачи, но не искажающие ее физическое содержание. Необходимо помнить, что сколь бы ни было точным математическое решение задачи, оно не может быть точнее тех приближенных предпосылок, на которых основано.

Основные предпосылки при решении задачи формирования и движения антропогенных выбросов в атмосфере являются следующие:

– выброс представляется в виде правильного геометрического тела. В случае взрывного выброса – это полусфера. В случае газообразного облака – это либо сфера, либо эллипсоид;

– твердая фаза выброса представляется в виде частиц сферической формы;

– вещество выброса считается равномерно распределенным по его объему, все его газодинамические характеристики осреднены по объему, а центр приложения массовых сил совпадает с геометрическим центром выброса.

Погрешности исходных данных.

При проведении экспериментов, обработке статистических данных и округлении результатов предыдущих расчетов приходится работать с приближенными числами. Присущие им погрешности обуславливают погрешности результата, причем для каждой операций по-своему. Анализируя точность задания исходных данных, можно сделать заключение о возможней точности результата вычислений.

Основными исходными данными задачи изучения атмосферных выбросов являются характеристика аварийного объекта и метео данные. Массовые характеристики горящих объектов или взрываемых изделий могут быть определены с точностью до 4^-5%; термодинамические и теплофизические характеристики продуктов горения или топлив – с точностью до 10%. Точность линейных размеров + 0,1 м. Метеорологические параметры определяются со следующей погрешностью:

– плотность воздуха + 0,01 кг/мЗ;

– скорость ветра ± 1м/с;

– направление ветра I град;

– градиенты скорости ветра и температуры атмосферного воздуха ± 0,01.

Методические погрешности

Методические погрешности при решении задачи возникают, при использовании приближенных формул вместо точных. Этот прием используется при замене функций их разложением в ряды ( при интерполяции скорости ветра и температуры атмосферного воздуха в виде полиномов по Z ), замене производной – разностью (при замене дифференциальных уравнений – конечноразностными ), при заменах интегралов – суммами и т.д. Используя приближенные формулы, применяемые многократно вместо сложной точной формулы, удается в процессе последовательных приближений найти эффективное решение задачи. При этом всегда надо знать диапазон аргумента, для которого они применимы, и точность вычисления по ним. В общем случае при вычислениях на ЭВМ погрешность метода целесообразно выбирать такой, чтобы она была в несколько раз меньше погрешности исходных данных.

Погрешности измерений дальномером

При наблюдении с поверхности земли геометрических характеристик, газообразных аварийных выбросов, а также скоростей их перемещения в пространстве как целого из-за погрешностей визуального определения границ их объемов могут быть допущены существенные ошибки в определении указанных величин. Подобные ошибки неизбежно появляются из-за кучевой неоднородной структуры поверхности газообразных выбросов, когда их отдельные выступающие части заслоняются другими [133, 158, 159].

Оценим ошибки в определении характеристик некоторых модельных выбросов.

Первичный полусферический выброс.

Случай приземного полусферического выброса реализуется в первые мгновения после взрывной аварии. Исследователь, наблюдающий взрывной выброс радиуса R под углом а к горизонту (Рис. 3.12), фиксирует видимый его размер г = R · Cos α .

Из геометрических построений этого рисунка легко получить для дефекта радиуса следующее выражение


ΔR = R– г = R(1 – Cos α)


На графике этого рисунка приводится зависимость относительной ошибки в определении радиуса выброса  от угла α. Как следует из этого графика  резко увеличивается при увеличении угла наблюдения α ≈ 18 ÷ 20°. Приемлемые (менее 5%) ошибки будут при α < 20°.

Сферический выброс.

Кратковременный выброс после отрыва от земли хорошо моделируется сферой. Он характеризуется радиусом R, геометрическими координатами центра масс X и (Z+R) (см. Рис. 3.13), а также углами α и β . Эти углы определяются видимыми наблюдателю из точки М частями поверхности выброса – дугами К А и КД. Как следует из этого рисунка наблюдатель видит выброс уменьшенным в вертикальном направлении на величину (hB +hH), а в горизонтальном направлении на величину (Хв + Хн ). Найдем эти характеристики из геометрических построений рисунка. Из рассмотрения прямоугольников МАО, АВО, СДО и ДОМ находим центральные углы АО В и СОД:


∠AОВ = 2α + β;

∠СОД = β.

Из прямоугольных треугольников АОВ и СОД получаем связь геометрических параметров в виде:



Откуда для наблюдаемых дефектов вертикальных и горизонтальных размеров сферического выброса находим:



Таким образом, вертикальный и горизонтальный видимые размеры сферического приподнятого выброса описываются следующими соотношениями:



Найдем относительные ошибки в визуальном определении геометрических характеристик выброса [133].

Из соотношений (3.47) – (3.50) находим относительную ошибку в определении вертикального  и горизонтального  размеров сферического приподнятого выброса. Получаем:



Рис. 3.12. Зависимость относительной ошибки в определении радиуса первичного выброса от угла наблюдений.



Найдем теперь относительную (относительно геометрической высоты) ошибку в определении верхней границы облака. Из соотношения (3.47) получаем:



Из прямоугольных треугольников МОЕ находим



Подставив соотношение (3.58) в формулу (3.55), получаем окончательное выражение для относительной ошибки в определении верхней граница сферического приподнятого выброса:



Аналогично находится относительная ошибка в определении нижней границы выброса. Получаем:



На графике рис. 3.14 представлено изменение относительной ошибки в определении вертикального размера приподнятого сферического выброса в зависимости от угла β для различных полууглов наблюдения α .

Как следует из этих графиков угол возвышения выброса над местностью β существенно влияет на величину  . Еще большее влияние на ошибку в определении высоты выброса оказывает величина полуугла его наблюдения α . Для α > 200 h превышает 20% для любых значений угла β .



Рис. 3.13. Схема для расчета ненаблюдаемых с земли и реальных размеров сферического приподнятого выброса.


Оптимальное соотношение углов α и β , обеспечивающее возможность четкого наблюдения деталей выброса и не более чем десятипроцентную ошибку в определении его высотного размера является β ≈ 15° и α ≈10°. Для нахождения видимых геометрических параметров наземного сферического выброса следует в полученных нами выражениях перейти к пределу при β → 0 . Получаем для дефектов вертикальных видимых размеров следующие выражения, (т.к. tg α = R / X):



Как следует из формулы (3.62) для наземного выброса его нижняя граница определяется точно.

Видимый вертикальный размер облака Hназемнвид запишется так:




Рис. 3.14. Изменение относительной ошибки в определении вертикального размера приподнятого сферического выброса в зависимости от угла над местностью β для различных полууглов наблюдения выброса α.



Рис. 3.15. Изменение относительной ошибки в определении вертикального размера выброса в зависимости от полуугола его наблюдения.


Относительную ошибку в определении вертикального размера облака находим из соотношения



Как видно из анализа соотношения (3.55) при β → 0 относительная ошибка в определении верхней границы облака Hназемнв для наземного сферического выброса в точности равна относительной ошибке визуального определения вертикального размера облака, т.е.



Как следует из этого графика при α > 200 ошибка в определении высоты выброса превосходит 10% и резко возрастает с увеличением полуугла наблюдения α . На практике для уменьшения ошибки в определении вертикального размера облака следует его наблюдение вести на относительно больших удалениях.

3.8. Высоты подъемов выбросов в атмосфере

Как было показано в главе 1, одним из основных параметров в рамках любой математической моде-|ли расчета концентраций загрязняющей примеси является высота вторичного атмосферного источника – фактически высота выброса в месте потери им динамической индивидуальности.

В большинстве современных разработок авторы пытаются использовать аналитические выражения для этого параметра, однако практика применения подобных формул имеет слишком малую область корректного использования в отношении как к тепловой мощности источника, так и к метеопараметрам.

Кроме того, часто путают динамический подъем выбросов с тепловым всплытием их разрушившихся объемов. Ошибочно считают, что тепловой подъем дает искомый результат, после чего наступает фаза атмосферной диффузии.

За границу струи, например, предлагается [136] принять изолинию однопроцентной относительной избыточной температуры.

Не всегда имеются и достаточно точные определения самого понятия подъема выброса. Например, применительно к струям факельного типа за такую высоту принимается [137,138] высоту струи, когда угол касательной к траектории ее наветренной части в сносящем ветре равен 100, в других работах за такую высоту предлагается считать подъем выброса на фиксированном расстоянии от трубы или его подъем за фиксированное время и т.д.

Некоторые авторы считают, что «потолок» выброса достигается, где он еще различим с помощью измерительной или фотографической аппаратуры.

Считается, что в случае когда радиоактивные или химические опасные вещества поступают в атмосферу посредством взрыва, можно пользоваться результатами работы Бриггса [139]. Однако результаты вычислений по приведенным там формулам также имеют весьма ограниченный диапазон применения. Поэтому рекомендуется, если это возможно, эффективную высоту источника загрязнений определять натуральными измерениями или оценкой.

Бриггс в зависимости от метеорологических условий предлагает проводить расчет подъема струи Δh по одной из нескольких модельных формул. Приведем их. Для устойчивого равновесия атмосферы предлагается выражение:



Значения параметра р, входящего в эту формулу, в зависимости от класса устойчивости атмосферы представлены в таблице 3.6.


Таблица № 3.6.

Скоростной параметр р в зависимости от устойчивости атмосферы и типа местности (по данным [162])



Для высот Z > 200 м следует брать постоянную скорость ветра на высоте 200 м.

Для условий слабых ветров подъема факела на завершающем этапе подъема предлагается находить по формуле:

Δh = 5,3 • F1/4 -S3/8 -R0,

где R0 – радиус дымовой трубы.

Для условий, близких к нейтральным, при которых параметр S приближается к нулю, Бриггс предлагает следующую формулу расчета конечного подъема факела:

Δh = l,54(F/U SU2X)2/3-hs1/3

где Uх – скорость трения; hs – высота дымовой трубы.

В работе [22] предлагается формула, объединяющая начальный поток количества движения

М0 = W20R20 и плавучий поток F:


Δh = 3,75M01/2 /U(10м)+ 5F/U(10м)3


где U(10м) – скорости ветра, измеренная на высоте флюгера.

Одной из первых формул, при составлении которой сделана попытка учесть вклад динамического и теплового подъема выброса, является формула Холланда (Holland) [116]:



где скорость примеси в выходном сечении W0 в м/с; диаметр устройства выброса D0 в метрах.

Эта формула получена в результате обработки наблюдений за шлейфами, выброшенными из трубы на высотах, не превышающих 50 м.

В ряде работ [137], полный подъем выброса в атмосфере предлагается разделить на динамический (газодинамический) Δhr и тепловой Δhm – за счет перегрева вещества выброса. Полный подъем загрязняющей примеси Δh при этом определяется сложением динамической и тепловой составляющих:

Δh = Δhm + Δhr

В частности, для модельных экспериментов для условий тепловых электростанций получены значения Δhm и Δhr в следующем виде [157]:



где Тг – температура отходящих газов;


ΔТ = Тr—Те;


εy – интенсивность турбулентности в горизонтальной плоскости; S = 0,05 ÷ 0,17

β – угол касательной к траектории движения факела; остальные обозначения те же, что и ранее.

Отмечается [137], что если данные по исследованию динамической составляющей подъема факела у разных авторов по характеру влияния на подъем параметров сносящего потока и движущейся в ней струи совпадает, то при изучении тепловой составляющей такое единообразие отсутствует. Оно проявляется в различиях в значениях показателей в формуле для Δhr, которое приводит к большим отличиям в вычисленных значениях теплового всплытия.

Такое различие вычисляемых Δhr объясняется

различными теоретическими предпосылками при определении этой характеристики, трудностями исследования теплового подъема на моделях в лабораторных условиях, отсутствием опытных данных по влиянию различных факторов на подъем выброса.

Что касается ограничений для применения аналитических формул, то они не пригодны в случаях сильно стратифицированной атмосферы или при сильном сдвиге ветра.

Рассмотрим теперь некоторые литературные данные по высотам подъемов кратковременных выбросов. По зарубежным литературным источникам, обобщенным в работе [151] для наземных, приземных и воздушных ядерных взрывов высота центра облака после стабилизации может быть найдена по формуле:

hц=1070 ·q02

где q – мощность ядерного заряда;

[q]=T; [h]=M

Модификация этой формулы относительно верхней he и нижней hH границ взрывного облака дает следующие значения:



где a = (3 + 0,131gq)-1; e = (2,6 +0,4lgq)-1.

Независимая обработка данных по высотам 60 ядерных взрывов привела к появлению формулы, справедливой в диапазоне q от долей тонны до 105т (с надежной статистикой лишь в диапазоне 1-100кт), определяющую высоту подъема hц в виде [151]:


hц =1600·q0,21


Для наземных подрывов взрывчатки справедливо соотношение


hц = 284·q022 ·В-1/3 -0,36u · В-1, (3.68)

где u – средняя скорость ветра в слое 0·hц .

 – параметр, определяющий устойчивость атмосферы.

В случае изотермий, когда  и U = 5 м/с, а также устойчивых, но близких к нейтральным условиях формула (3.68) принимает вид [151]:

Δhц = 1400 • q0,22 – 52 (3.69)

Эта формула неприменима при , т.е. для неустойчивой стратификации атмосферы, при которой подъем выброса за счет сил плавучести теоретически ничем не ограничен.

При взрывах химических ВВ в серийных взрывах программы «Хардхет» в умеренноустойчивых условиях высота подъема облака оценивается в следующем виде [151]: hц ≈ 700 • qn, (3.70)

где n – принимает значения от 0,2 до 0,25.

Для более полного описания геометрии атмосферного источника при ядерных взрывах целесообразно привести формулы геометрических характеристик подобных источников. Для диаметра облака Добл, вертикального его размера ΔН и диаметра «ножки» облака Дн можно пользоваться следующими оценками:

Добл. = 1600 • q0117

ΔН = 1430 V0246

Дн= 1420 V0134

Эти же значения параметров, очевидно, могут быть применены для инженерных оценок выбросов при взрывах химических ВВ и авариях взрывного характера на АЭС и других энергоемких объектов. В любом случае после взрыва формируется универсальный по форме атмосферный источник, отличающийся лишь характером поступления примесей и их составом.

В разделе книги, посвященном рассеиванию примесей из вторичных источников, приводится пример, как используя стандартные модели рассеяния, можно получить суммарное поле приземных концентраций в виде суперпозиции концентраций двух источников: облака и «ножки».

Отметим, что приведенные в этом разделе формулы пригодны только для весьма грубых инженерных оценок в соответствующих диапазонах параметров атмосферы и источника загрязнений.

Наилучшие и наиболее корректные результаты в процедуре получения высот подъема выбросов в атмосфере дает, на наш взгляд, решение полной системы дифференциальных уравнений движения выброса при задании пульсационных характеристик атмосферы.

При этом для прерывания расчетного процесса необходимо использовать критерий потери выбросом динамической индивидуальности на фоне пульсационного движения окружающего воздуха.


Страницы книги >> Предыдущая | 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | Следующая
  • 0 Оценок: 0

Правообладателям!

Это произведение, предположительно, находится в статусе 'public domain'. Если это не так и размещение материала нарушает чьи-либо права, то сообщите нам об этом.


Популярные книги за неделю


Рекомендации