Транспорт наносов внутренними волнами

In the Boussinesq approximation, the mean flows induced by a wave through nonlinearities in the presence of eddy viscosity are determined for internal waves. The boundary layer solutions, decrement of fading of a wave, and tangential stress at the bottom are determined. The bottom concentration of w...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2007
Hauptverfasser: Слепышев, A.A., Носова, А.В.
Format: Artikel
Sprache:Russian
Veröffentlicht: Видавничий дім "Академперіодика" НАН України 2007
Schlagworte:
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/3140
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Транспорт наносов внутренними волнами / А.А. Слепышев, А.В. Носова // Доп. НАН України. — 2007. — № 9. — С. 95-101. — Бібліогр.: 9 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-3140
record_format dspace
spelling irk-123456789-31402009-07-02T12:00:22Z Транспорт наносов внутренними волнами Слепышев, A.A. Носова, А.В. Науки про Землю In the Boussinesq approximation, the mean flows induced by a wave through nonlinearities in the presence of eddy viscosity are determined for internal waves. The boundary layer solutions, decrement of fading of a wave, and tangential stress at the bottom are determined. The bottom concentration of wave-suspended sediments is obtained when the bottom stress exceeds a critical value. The vertical distribution of the concentration of sediments is determined in the diffusion approximation. 2007 Article Транспорт наносов внутренними волнами / А.А. Слепышев, А.В. Носова // Доп. НАН України. — 2007. — № 9. — С. 95-101. — Бібліогр.: 9 назв. — рос. 1025-6415 http://dspace.nbuv.gov.ua/handle/123456789/3140 551.466.8 ru Видавничий дім "Академперіодика" НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
topic Науки про Землю
Науки про Землю
spellingShingle Науки про Землю
Науки про Землю
Слепышев, A.A.
Носова, А.В.
Транспорт наносов внутренними волнами
description In the Boussinesq approximation, the mean flows induced by a wave through nonlinearities in the presence of eddy viscosity are determined for internal waves. The boundary layer solutions, decrement of fading of a wave, and tangential stress at the bottom are determined. The bottom concentration of wave-suspended sediments is obtained when the bottom stress exceeds a critical value. The vertical distribution of the concentration of sediments is determined in the diffusion approximation.
format Article
author Слепышев, A.A.
Носова, А.В.
author_facet Слепышев, A.A.
Носова, А.В.
author_sort Слепышев, A.A.
title Транспорт наносов внутренними волнами
title_short Транспорт наносов внутренними волнами
title_full Транспорт наносов внутренними волнами
title_fullStr Транспорт наносов внутренними волнами
title_full_unstemmed Транспорт наносов внутренними волнами
title_sort транспорт наносов внутренними волнами
publisher Видавничий дім "Академперіодика" НАН України
publishDate 2007
topic_facet Науки про Землю
url http://dspace.nbuv.gov.ua/handle/123456789/3140
citation_txt Транспорт наносов внутренними волнами / А.А. Слепышев, А.В. Носова // Доп. НАН України. — 2007. — № 9. — С. 95-101. — Бібліогр.: 9 назв. — рос.
work_keys_str_mv AT slepyševaa transportnanosovvnutrennimivolnami
AT nosovaav transportnanosovvnutrennimivolnami
first_indexed 2025-07-02T06:25:29Z
last_indexed 2025-07-02T06:25:29Z
_version_ 1836515358790385664
fulltext 11. Братусь М.Д., Даниш В.В., Сворень Й.М. Вуглеводневi сполуки гiдротермальних утворень Кар- пат // Там само. – 1981. – № 7. – С. 3–6. 12. Дудок И.В. Минералого-геохимические особенности жильных образований флишевых отложений Украинских Карпат (в связи с нефтегазоносностью) : Автореф. дис. . . . канд. геол.-мин. наук / АН УССР. Ин-т геологии и геохимии горюч. ископаемых – Львов, 1991. – 19 с. 13. Чекалюк Э.Б., Филяс Ю.И. Водо-нефтяные растворы. – Киев: Наук. думка, 1977. – 128 с. 14. Ладыженский Н. Р. К вопросу о времени формирования нефтяных месторождений Карпат // Геол. сб. Львов. геол. о-ва. – 1961. – 7./8. – С. 79–88. 15. Сворень Й.М., Наумко I.М. Нова теорiя синтезу i генезису природних вуглеводнiв: абiогенно-бiоген- ний дуалiзм // Доп. НАН України. – 2006. – № 2. – С. 111–116. Надiйшло до редакцiї 12.03.2007Iнститут геологiї i геохiмiї горючих копалин НАН України, Львiв УДК 551.466.8 © 2007 A.A. Слепышев, А.В. Носова Транспорт наносов внутренними волнами (Представлено членом-корреспондентом НАН Украины Л. В. Черкесовым) In the Boussinesq approximation, the mean flows induced by a wave through nonlinearities in the presence of eddy viscosity are determined for internal waves. The boundary layer solutions, decrement of fading of a wave, and tangential stress at the bottom are determined. The bottom concentration of wave-suspended sediments is obtained when the bottom stress exceeds a critical value. The vertical distribution of the concentration of sediments is determined in the diffusion approximation. Транспорт наносов обычно связывают с поверхностными волнами, однако их влияние рас- пространяется до глубин, которые составляют половину длины волны [1], т. е. 20–30 м в Чер- ном море. На больших глубинах отмечается сильное проявление внутренних волн. Опре- делению среднего течения, индуцируемого пакетом внутренних волн за счет нелинейности без учета турбулентной вязкости и диффузии посвящена статья [2]. Вертикальные потоки тепла, соли, импульса с учетом турбулентной вязкости находим в работе [3]. В настоящем сообщении описано определение горизонтальных средних течений, индуцированных волной за счет нелинейности с учетом турбулентной вязкости и диффузии; в диффузионном при- ближении — нахождение вертикального распределения концентрации наносов, взвешенных волной при превышении тангенциальными напряжениями критических значений, соответ- ствующих началу движения наносов. Постановка задачи. Рассматриваются свободные внутренние волны в приближении Буссинеска при реальной стратификации и постоянных коэффициентах турбулентного об- мена. Исходная нелинейная система уравнений гидродинамики для волновых возмущений с учетом турбулентной вязкости и диффузии решается асимптотическим методом много- масштабных разложений [1]. В первом порядке малости по крутизне волны находим реше- ние линейного приближения и дисперсионное соотношение. Во втором порядке малости по ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2007, №9 95 крутизне волны решаем краевую задачу по определению вертикальной структуры неосцил- лирующей поправки к функции тока, обуловленной нелинейностью, т. е. находим среднее течение, индуцированное волной. Принимая в качестве исходных уравнений для волновых возмущений уравнения Навье–Стокса для неоднородной жидкости, введем безразмерные переменные: t̃ = t ω∗ , k̃ = k H , ω̃ = ω∗ω, ũ1 = u1Hω∗, ũ3 = u3Hω∗, P̃ = ρ0H 2ω2 ∗ P, ρ̃ = ρ0ω 2 ∗ Hρ g , x̃i = Hxi, K̃i = Kiµ, M̃i = Miµ (i = 1, 3), где t̃, k̃, ω̃, ũ, P̃ , ρ̃, x̃, K̃, M̃ — размерные величины. Систему уравнений гидродинамики для волновых возмущений в безразмерных перемен- ных в приближении Буссинеска представим так: ∂u1 ∂t + ui ∂u1 ∂xi = − ∂P ∂x1 + ǫ22K1 ∂2u1 ∂x2 1 + ǫ22 ∂K1 ∂x3 ∂u3 ∂x1 + ǫ22K3 ∂2u1 ∂x2 3 + ǫ22 ∂K3 ∂x3 ∂u1 ∂x3 , ∂u3 ∂t + ui ∂u3 ∂xi = − ∂P ∂x3 + ǫ22K1 ∂2u3 ∂x2 1 + ǫ22 K3∂ 2u3 ∂x2 3 + 2ǫ22 ∂K3 ∂x3 ∂u3 ∂x3 − ρ, ∂ρ ∂t + ui ∂ρ ∂xi = ǫ22M1 ∂2ρ ∂x2 1 + ǫ22 ∂ ∂x3 ( M3 ∂ρ ∂x3 ) − u3 dρ0 dx3 , ∂ui ∂xi = 0. (1) Здесь (общие для приведенных формул обозначения) g — ускорение силы тяжести; x1, x3 — горизонтальная и вертикальная координаты соответственно, вертикальная ось направлена вверх; ρ и P — волновые возмущения плотности и давления; ρ0 — характерная средняя плот- ность воды; u1 и u3 — горизонтальная и вертикальная компоненты волновых возмущений скорости; K1, K3, M1, M3 — горизонтальные и вертикальные коэффициенты турбулентной вязкости и диффузии соответственно; H — глубина моря; ω∗ — характерная частота волны; ǫ22 = µ/(H2ω∗) — малый параметр, пропорциональный значению горизонтальной турбулент- ной вязкости. В качестве граничных условий на свободной поверхности используем кинематическое и динамическое условия: −P + ζ3g1 + 2ǫ22K3 ∂u3 ∂x3 = 0, (2) K3 ∂u1 ∂x3 +K1 ∂u3 ∂x1 = 0, dζ3 dt = u3, (3) где ζ3 — вертикальное смещение свободной поверхности; g1 = g/(ω2 ∗ H). Первые два условия определяют отсутствие нормальных и тангенциальных напряжений на свободной поверх- ности. На дне примем условия прилипания: u3 = u1 = 0 при x3 = −1. (4) Граничные условия по плотности опишем таким образом: ρ+ ζ3 ∂ρ0 ∂x3 + ζi ∂ρ ∂xi = 0 при x3 = 0, (5) ρ = 0, x3 = −1. (6) 96 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2007, №9 Пусть ψ(x1, x3, t) — функция тока, которая определяет поле волновых орбитальных ско- ростей. Волновые возмущения компонент скорости выразим через ψ: u1 = ∂ψ ∂x3 , u3 = − ∂ψ ∂x1 . (7) Систему уравнений (1) после данной подстановки преобразуем к виду ∂∆ψ ∂t + Jx1x3 (∆ψ,ψ) = ∂ρ ∂x1 + ǫ22 ∂ ∂x1 [ ∂ ∂x1 ( K1 ∂2ψ ∂x2 1 ) + ∂ ∂x3 ( K3 ∂2ψ ∂x1∂x3 )] + + ǫ22 ∂ ∂x3 [ ∂ ∂x1 ( K1 ∂2ψ ∂x1∂x3 ) + ∂ ∂x3 ( K3 ∂2ψ ∂x2 3 )] , (8) ∂ρ ∂t + Jx1x3 (ρ, ψ) − ǫ22 ∂ ∂x1 ( M1 ∂ρ ∂x1 ) − ǫ22 ∂ ∂x3 ( M3 ∂ρ ∂x3 ) − dρ0 dx3 ∂ψ ∂x1 = 0. (9) где Jx1x3 (a, b) = ∂a ∂x1 ∂b ∂x3 − ∂a ∂x3 ∂b ∂x1 — якобиан по переменным x1, x3. Данную систему будем решать асимптотическим методом многомасштабных разложе- ний [1] с введением “медленных” переменных ζ, τ (ζ = ǫ2x1, τ = ǫ2t) и “быстрой” переменной θ-фазы волны (k = ∂θ/∂x1, ω = −∂θ/∂t): ψ = ∑ n=1 ǫnψn(ζ, τ, x3, θ), ρ = ∑ n=1 ǫnρn(ζ, τ, x3, θ). (10) Решение уравнений линейного приближения будем искать в виде: ψ1 = Aϕ1(x3)e iθ + k.c., ρ1 = An1(x3)e iθ + k.c., (11) где k.c. — комплектно сопряженное слагаемое. Подставляя функции (11) в (10) и уравнения (8), (9), получим уравнения для ϕ1 и связь между n1(x3) и ϕ1(x3): k2 dρ0 dx3 ϕ1 = [ −ωi+ ǫ22M1k 2 − ǫ22 d dx3 ( M3d/dx3 )][ ωi ( d2ϕ1 dx2 3 − k2ϕ1 ) − − ǫ22 d dx3 [ K1k 2 dϕ1 dx3 − d dx 3 ( K3 d2ϕ1 dx2 3 )] + ǫ22k [ K1k 3ϕ1 −K3k d2ϕ1 dx2 3 ]] , (12) ( iω − k2M1ǫ 2 2 + ǫ22 d dx3 ( M3 d dx3 )) n1 = −ik dρ0 dx3 ϕ1. (13) Уравнения (12) и (13) следует дополнить граничными условиями, которые вытекают из условий (2)–(6): kg1 ω ϕ1 + ω k dϕ1 dx3 − ikK1ǫ 2 2 dϕ1 dx3 + i k ǫ22 d dx3 ( K3 d2ϕ1 dx2 3 ) − 2ikǫ22K3 dϕ1 dx3 = 0 при x3 = 0, (14) K3 d2ϕ1 dx2 3 +K1k 2ϕ1 = 0 при x3 = 0, (15) ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2007, №9 97 ϕ1 = dϕ1 dx3 = 0 при x3 = −1, (16) n1 + k ω dρ0 dx3 ϕ1 = 0 при x3 = 0, (17) n1 = 0 при x3 = −1. (18) Уравнение (12) при малом ǫ2 будем решать асимптотическим методом Люстерника– Вишика [4, 5], разлагая ϕ1, n1, ω в асимптотические ряды: ϕ1(x3) = ∑ i=0 ϕ1i(x3)ǫ i 2 + ǫ2 ∑ i=0 ǫi2ν 1 i ( x3 + 1 ǫ2 ) + ǫ22 ∑ i=0 ǫi2ν 0 i ( x3 ǫ2 ) , (19) n1(x3) = ∑ i=0 n1iǫ i 2 + ǫ2 ∑ i=0 ǫi2w 1 i ( x3 + 1 ǫ2 ) + ǫ22 ∑ i=0 ǫi2w 0 i ( x3 ǫ2 ) , (20) ω = ω01 + ∑ i=1 ǫi2ωi1. (21) Здесь ν1 i ((x3 +1)/ǫ2) — погранслойные решения в окрестности дна; ν0 i (x3/ǫ2) — погранслой- ные решения в окрестности свободной поверхности. Введение погранслойных решений не- обходимо для удовлетворения граничных условий (2)–(6). Подставляя разложения (19)–(21) в уравнение (12), получим краевую задачу для ϕ10, определяющую вертикальную структуру моды в линейном приближении, краевую задачу для ϕ12, из условия разрешимости которой находится декремент затухания волны на турбулентности δω = ǫ22ω21/i [3] и погранслойные решения в окрестности верхней и нижней границы [3]. Уравнение для неосциллирующей на временном масштабе волны поправки к функции тока находим из уравнения второго приближения по параметру ǫ, осредненного по периоду волны. Как и при отсутствии тур- булентности, неосциллирующую поправку к функции тока C(x3, τ, ǫ) следует искать в виде C = c(x3)A1A ∗ 1, где A1 = A exp(δω · t). Функция c(x3) удовлетворяет краевой задаче ǫ22 d2 dx2 3 ( K3 d2c dx2 3 ) = ki d dx3 (( −k2 + d2 dx2 3 ) ϕϕ∗ 1 ) + k.c., ikϕ∗ 1 d2ϕ1 dx2 3 + k. c. = ǫ22 d dx3 ( K3 d2c dx2 3 ) при x3 = 0, d2c dx2 3 = 0 при x3 = 0, dc dx3 = c = 0 при x3 = −1. (22) Горизонтальная компонента скорости индуцированного волной среднего течения опре- делялась через функцию c(x3) : uинд = |ǫA1| 2dc/dx3. Для оценки коэффициента верти- кального турбулентного обмена применим формулу Озмидова [6] для скорости диссипа- ции турбулентной энергии ǫt; K3 выражается через ǫt и критический масштаб плавучести L0 : ǫ 1/2 t = L0N 3/2, K3 = c1ǫ 1/3 t L 4/3 0 . Для оценки критического масштаба плавучести L0 применим формулу Бэлла [7]: L0 = N∗/(Nβ0), отсюда следует зависимость ǫt и K3 от N : ǫt = N2 ∗ N β0 , K3 = c1 N2 ∗ β2 0 N , (23) 98 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2007, №9 где N∗ = 3 цикл/ч; β0 = 1 м−1; c1 = 0,1. Будем полагать, что теряемая волной энергия цели- ком переходит в турбулентность и далее расходуется на работу турбулентности против сил плавучести и на диссипацию в тепло, т. е. скорость диссипации волновой энергии, проинтег- рированная по глубине, равна интегральной величине работы турбулентности против сил плавучести M3N 2 и скорости диссипации турбулентной энергии ǫt, т. е. 2|δω| H∫ 0 Edz = H∫ 0 (M3N 2 + ǫt) dz. (24) Здесь E = ǫ2A1A ∗ 1 [ k2 ( 1+2 N2 ω2 ) ϕ1ϕ ∗ 1+ dϕ1 dz dϕ∗ 1 dz ] — плотность энергии волны, z = (x3+1)H. Уравнение (24) позволяет найти коэффициент горизонтального турбулентного обмена K1. Определим осредненное за период волны тангенциальное напряжение у дна: τ = ρ0K3 √( ∂u1 ∂z )2 . Если тангенциальное напряжение у дна превышает критическое значение τ0, соответствую- щее началу движения наносов, то волна взмучивает наносы, осуществляя их горизонталь- ный перенос. В стационарном и горизонтально-однородном случае уравнение вертикальной диффузии для средней концентрации наносов n(z) имеет вид [8]: ∂ ∂z ((w − wg)n(z)) = ∂ ∂z ( M3 ∂n ∂z ) , (25) где wg — гидравлическая крупность наносов; w — вертикальная компонента скорости те- чения, индуцированного волной за счет нелинейности. Решение уравнения (25), затухающее при удалении от дна, приведем к виду n(z) = n0 exp ( z∫ 0 w − wg M3 dz ) . (26) Здесь n0 = χ(τ − τ0)/wg [9]. Расход наносов определяется по формуле G = H∫ 0 n(uинд(z) + u1s) dz, (27) где u1s — скорость стоксова дрейфа [3]. Результаты расчетов. Cредние течения, индуцированные волной за счет нелиней- ности рассчитывали на Северо-Западном шельфе Черного моря при стратификации, по- казанной на рис. 1. Волновое число у получасовых внутренних волн низшей моды k = = 1,44 ·10−2 м−1, декремент затухания волны δw = −1,34 ·10−4 рад/с. Коэффициент верти- кальной турбулентной вязкости определяли по формуле (23), коэффициент горизонтальной турбулентной вязкости — из уравнения (24) при условии, что M3 = 0,5K3, M1 = 0,5K1. Нор- мирующий множитель ǫA1 находился по известной максимальной амплитуде вертикальных ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2007, №9 99 Рис. 1. Средний профиль частоты Брента–Вяй- сяля Рис. 2. Вертикальный профиль горизонтальной эйлеровой скорости среднего течения, индуциро- ванного волной Рис. 3. Вертикальный профиль концентрации наносов, взвешенных волной Рис. 4. Зависимость критической амплитуды волны, соответствующей началу движения нано- сов от частоты волны смещений ζ0. Вертикальная структура индуцированного течения для волны с w = 2 цикл/ч, ζ0 = 5 м показана на рис. 2. Средневзвешенный диаметр донных осадков рассматриваемой акватории Черного моря составляет d = 0,03 мм [9]. Гидравлическая крупность наносов wg со средневзвешенным диаметром d = 0,03 мм равна 0,078 см/с, коэффициент χ в форму- ле (26) для донной концентрации наносов равен 5·10−6 с/см [8]. Критическое тангенциальное напряжение, соответствующее началу движения наносов, τ0 = 1 дин/см2 [8]. Максималь- ная амплитуда волны, при которой тангенциальное напряжение (23) равно критическому, ζ0кр = 4,5 м. Если максимальная амплитуда волны ζ0 = 5 м, то величина тангенциального напряжения у дна будет τ = 1,12 дин/см2. Донная концентрация взвешенных волной на- носов при ζ0 = 5 м составляет n0 = 7,4 мг/л. Коэффициент горизонтальной турбулентной вязкости K1 = 0,86 м2/c. Вертикальное распределение концентрации наносов, взвешенных волной, показано на рис. 3. Зависимость критической амплитуды волны низшей моды, при которой начинается движение наносов от частоты волны изображена на рис. 4. Величи- на расхода наносов (27) составляет G = 5,3 · 10−5 кг/(м · c). Интегральный поток наносов положительный, т. е. сонаправлен с горизонтальным волновым вектором и направлением распрастранения волны. 100 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2007, №9 1. Ястребов В. С., Парамонов А.Н. и др. Исследование природного слоя буксируемыми аппаратами. – Москва: Ин-т океанологии АН СССР, 1989. – 128 с. 2. Борисенко Ю.Д., Воронович А.Г., Леонов А.И., Миропольский Ю.З. К теории нестационарных сла- бонелинейных внутренних волн в стратифицированной жидкости // Изв. АН СССР. Физика атмо- сферы и океана. – 1976. – 12, № 3. – С. 293–301. 3. Пантелеев Н.А., Слепышев А.А. Тепломассоперенос слабонелинейными внутренними волнами при наличии турбулентности // Мор. гидрофиз. журн. – 1995. – № 4. – С. 3–23. 4. Задорожный А.И. Затухание длинных волн в экспоненциально стратифицированном море // Мор. гидрофиз. исследования. – 1975. – № 43. – С. 96–110. 5. Черкесов Л. В. Гидродинамика волн. – Киев: Наук. думка, 1980. – 259 с. 6. Озмидов Р.В. О турбулентном обмене в устойчиво стратифицированном море // Изв. АН СССР. Физика атмосферы и океана. – 1965. – 1, № 8. – С. 853–860. 7. Bell I. H. Internal wave-turbulence interpretation of ocean fine structure // Geophys. Res. Lett. – 1974. – No 6. – P. 253–255. 8. Шапиро Г.И., Аквис Т.М., Пыхов Н.В., Анциферов С.М. Перенос мелкодисперсного осадочного материала мезомасштабными течениями в шельфово-склоновой зоне моря // Океанология. – 2000. – 40, № 3. – С. 333–339. 9. Щербаков Ф.А., Куприн П.Н., Потапова Л.И., Поляков А.С., Забелина Э.К., Сорокин В.М. Осад- конакопление на континентальной окраине Черного моря. – Москва: Наука, 1978. – 210 с. Поступило в редакцию 23.03.2007Морской гидрофизический институт НАН Украины, Севастополь УДК 563.6:551.735(477.8) © 2007 В.Ф. Шульга, В.В. Огарь Первые находки коралловых построек в раннем карбоне Львовского палеозойского прогиба (Представлено членом-корреспондентом НАН Украины А.Ю. Митропольским) Organic buildups (biostromes and bioherms) in the coal formation of the Lvov Palaeozoic Trough are described for the first time. The primary framework builders were the colonial rugose corals. The coral buildups were formed in a shallow water environment and located wi- thin the tectonically active area. It is a narrow stripe of 10–12 km in width and more than 70 km in length. The late Viséan coral bioherms and biostromes under consideration resemble those from some regions of Western Europe. В раннекаменноугольную эпоху кораллы, наряду с водорослями, мшанками, криноидеями и др., являлись рифообразующими организмами [1, 2]. Однако, несмотря на значительное развитие в карбоне Львовского палеозойского прогиба (ЛПП) колониальных ругоз, состав- ляющих более половины коралловой фауны [3–5], а также широкое развитие двух мощных (до 140 м) известняковых толщ верхнего визе (олесковская и устилужская свиты), до са- мого последнего времени коралловые постройки на территории прогиба не были обнару- жены. По-нашему мнению, с одной стороны, это объясняется недостаточным вниманием к данному вопросу, а с другой — ограниченным количеством буровых скважин, вскрывших залегающие на больших глубинах вышеуказанные известняковые толщи. ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2007, №9 101