УДК 551.467; 551.321.62
DOI: 10.33764/2618-981Х-2019-2-2-187-194
ОПРЕДЕЛЕНИЕ ТОЛЩИНЫ ЛЕДЯНОГО ПОКРОВА ПО АКУСТИЧЕСКИМ ШУМАМ
Константин Владимирович Федин
Институт нефтегазовой геологии и геофизики им. А. А. Трофимука СО РАН, 630090, Россия, г. Новосибирск, пр. Академика Коптюга, 3, кандидат технических наук, научный сотрудник лаборатории динамических проблем сейсмики, тел. (383)333-34-19, e-mail: FedinKV@ipgg.sbras.ru; Новосибирский национальный исследовательский государственный университет, 630090, Россия, г. Новосибирск, ул. Пирогова, 2, старший преподаватель кафедры геофизики; Новосибирский государственный технический университет, 630073, Россия, г. Новосибирск, пр. К. Маркса, 20, доцент кафедры геофизических систем
Юрий Иванович Колесников
Институт нефтегазовой геологии и геофизики им. А. А. Трофимука СО РАН, 630090, Россия, г. Новосибирск, пр. Академика Коптюга, 3, доктор технических наук, доцент, главный научный сотрудник лаборатории динамических проблем сейсмики, тел. (383)333-31-38, e-mail: KolesnikovYI@ipgg.sbras.ru; Сейсмологический филиал Федерального исследовательского центра «Единая геофизическая служба РАН», 630090, Россия, г. Новосибирск, пр. Академика Коптюга, 3, ведущий геофизик
Приведены результаты натурных экспериментов, направленных на оценку возможностей определения толщины ледяного покрова водоемов по зарегистрированным на поверхности льда акустическим шумам. Показано, что по частотам генерируемых шумами в слое льда стоячих волн вертикального сжатия-расширения уверенно определяются как толщина ледяного покрова, так и тип подстилающей среды (вода или мерзлый грунт).
Ключевые слова: ледяной покров водоемов, определение толщины льда, акустические шумы, стоячие волны вертикального сжатия-растяжения.
DETERMINATION OF ICE COVER THICKNESS USING ACOUSTIC NOISE
Konstantin V. Fedin
Trofimuk Institute of Petroleum Geology and Geophysics SB RAS, 3, Prospect Аkademik Koptyug St., Novosibirsk, 630090, Russia, Ph. D., Researcher of Laboratory of Seismic Dynamic Analysis, phone: (383)333-34-19, e-mail: FedinKV@ipgg.sbras.ru; Novosibirsk National Research State University, 2, Pirogova St., Novosibirsk, 630073, Russia, Senior Lecturer of Geophysics Department; Novosibirsk State Technical University, 20, Prospect K. Marx St., Novosibirsk, 630073, Russia, Associate Professor of Department of Geophysical Systems
Yury I. Kolesnikov
Trofimuk Institute of Petroleum Geology and Geophysics SB RAS, 3, Prospect Аkademik Koptyug St., Novosibirsk, 630090, Russia, D. Sc., Associate Professor, Chief Researcher of Laboratory of Seismic Dynamic Analysis, phone: (383)333-31-38, e-mail: KolesnikovYI@ipgg.sbras.ru; Seis-mological Branch of Federal Research Center «United Geophysical Survey RAS», 3, Prospect Аkademik Koptyug, Novosibirsk, 630090, Russia, Lead Geophysicist
The results of field experiments aimed at assessing the possibilities of determining the thickness of the ice cover of water bodies using acoustic noise recorded on the ice surface are presented.
It is shown that the frequencies of vertical compressional standing waves generated by acoustic noise in the ice layer can be used for determining the thickness of the ice cover and the type of underlying medium (water or frozen soil).
Key words: ice cover of water bodies, determination of ice thickness, acoustic noise, vertical compressional standing waves.
Необходимость измерения толщины, образующегося над водной поверхностью при отрицательных температурах ледяного покрова может возникать при решении множества практических задач - обеспечения безопасного передвижения людей по замерзшим водоемам и транспортных средств по ледовым переправам и автозимникам, прогнозирования паводков, оценки качества ледовых карт, создаваемых на основе спутниковых методов дистанционного мониторинга ледовой обстановки и т.д. Кроме водоемов, иногда представляет интерес и оценка толщины льда, образовавшегося на суше, например, для оценки запасов ледниковой воды [11].
Среди методов определения толщины ледяного покрова самым точным, но и самым трудоемким, является прямой метод, предполагающий бурение лунок и измерение толщины льда ледомерной рейкой. При мониторинговых работах такие измерения проводят дистанционными методами, например, с использованием авиации [5] или космических аппаратов [3], но точность этих методов относительно невысока. При движении судоходного транспорта во льдах применяются визуальные или телевизионные наблюдения за обломками льдин, занимающих возле судна положение близкое к вертикальному, что дает возможность оценить толщину их бокового скола [4].
Наиболее распространенный метод определения толщины льда основан на электромагнитном зондировании ледяного покрова георадаром [11, 1, 7]. Этот метод характеризуется высокой производительностью и достаточно хорошей точностью, но при большой изменчивости электрических параметров льда (например, в соленой воде) он может требовать периодической калибровки аппаратуры.
В последнее время предлагаются и новые методы, основанные на анализе интерференционных диаграмм, создаваемых сигналами навигационных спутников ГЛОНАСС и GPS [2], на оптических эффектах (с использованием лида-ра) [10] или на анализе высоты звука от возбуждаемой в ледовой пластине волны Лэмба [12].
Акустические методы также применяются для определения толщины льда. Как правило, это активные методы, например, эхолокация с поверхности льда [8] или из-под воды с помощью сонара [13]. Ниже представлены результаты применения для определения толщины льда пассивного акустического метода, основанного на выделении стоячих волн из записанного на поверхности льда акустического шума.
Этот метод уже неоднократно применялся нами ранее, например, при проведении физического моделирования и натурных экспериментов по определению пустот в грунтовых отложениях [9]. Метод сводится к регистрации акусти-
ческих шумов на поверхности исследуемого ограниченного объекта и накоплению амплитудных спектров большого числа шумовых записей. Это позволяет выделять из шума, формирующиеся под его воздействием в объекте, стоячие волны.
В нашем случае таким ограниченным объектом является слой льда, лежащий на поверхности воды или мерзлого грунта. В таком слое при формировании стоячих волн между его нижней и верхней границами должно укладываться, в зависимости от условий отражения на них, либо целое число полудлин, либо нечетное число четвертей длин стоячих волн (подобно стоячим волнам в незакрепленном или закрепленном с одного конца стержнях [6]).
Частоты стоячих волн вертикального сжатия-растяжения в слое (собственные частоты слоя) в этих двух случаях определяются, соответственно, формулами
/. = £ (1)
для льда, под которым находится вода или воздух, или
к = (2)
для льда, лежащего поверх мерзлого грунта. Здесь п - номер моды стоячих волн, Ур - скорость продольных волн, к - расстояние между границами слоя.
В обоих случаях интервал между соседними собственными частотами равен ДF = (/п+1 —/п) = Ур/2И. На осредненных в результате накопления амплитудных спектрах это проявляется в появлении регулярной последовательности резонансных пиков на частотах стоячих волн. Таким образом, при известной скорости Ур и частот каких-либо двух соседних (по порядку моды) стоячих волн вертикального сжатия-растяжения толщина слоя льда может быть определена по формуле
К = ^ . (3)
2ДF 4 '
Заметим также, что для формулы (1) ДF = Д, а для формулы (2) ДF = 2Д. Поэтому вычислять толщину слоя можно и по частоте низшей моды, имея ввиду вид зависимости, которой соответствует распределение резонансных пиков на оси частот. То есть, над водой (или над воздухом в случае зависания льда) интервал ДF равен частоте низшей моды, а над грунтом - ее удвоенной величине.
Для оценки возможностей применения для определения толщины льда пассивного акустического метода, основанного на выделении стоячих волн из акустического шума, в конце января - начале февраля 2019 г. были проведены натурные эксперименты на двух водоемах в районе новосибирского Академгородка - на р. Зырянка и на пляже «Звезда» на Обском водохранилище (рис. 1).
Регистрация шумовых записей проводилась на линейных профилях с шагом 1 м. На реке профиль длиной 20 м был ориентирован вдоль берегов и находился примерно в ее средней части, причем заведомо над водой, так как из-подо льда отчетливо было слышно ее журчание. На пляже «Звезда» профиль наблюдений был примерно перпендикулярен к береговой линии и располагался частично над водой, частично над мерзлым песком.
Рис. 1. Положение участков наблюдений на р. Зырянка (1) и на пляже «Звезда» на Обском водохранилище (2)
Регистрация акустических шумов производилась с помощью цифрового осциллографа В-423 с частотой дискретизации 100 кГц. В качестве приемника использовался широкополосный пьезокерамический датчик поршневого типа с вертикально направленной осью максимальной чувствительности, который при измерениях устанавливался непосредственно на очищенную поверхность льда. Общая длительность шумовых записей в каждой точке составляла 30 секунд при наблюдениях на р. Зырянка и 2 минуты на пляже «Звезда», где течение практически отсутствует, и уровень шума, соответственно, был значительно ниже.
При обработке записи разбивались на фрагменты длительностью по 8192 отсчетов (примерно 82 мс), после чего проводилось накопление амплитудных спектров этих фрагментов. Частоты выделяемых на осредненных амплитудных спектрах резонансных пиков использовались для определения по формуле (3) толщины ледяного покрова в точках наблюдений. При этом предполагалось, что скорость продольных волн Ур для льда известна и в разных точках наблюдений меняется незначительно. Поэтому при расчетах использовалось значение Ур = 4090 м/с, определенное импульсным методом на керне льда, полученном при бурении контрольной лунки в районе пляжа «Звезда».
На рис. 2 приведены примеры осредненных спектров шумовых записей, зарегистрированных на двух участках наблюдений. Как можно видеть, даже
при относительно малой длительности регистрации шумов на их спектрах можно уверенно выделить несколько регулярных резонансных пиков. Регулярность пиков на амплитудных спектрах и согласие с формулами (1) и (2) позволяет идентифицировать эти пики как резонансы на частотах стоячих волн.
1 i
ч
и 0.8
° 0.6 га
j? 0 4 li.
I 0 2
IvMwi
10 20 30 40 50 60 70 Частота, кГц
а)
10 20 30 40 50 60 70 Частота, кГц
б)
0 10 20 30 40 50 60 70 Частота, кГц
Рис. 2. Примеры нормированных осредненных амплитудных спектров шумовых записей, зарегистрированных на ледяном покрове р. Зырянка (а) и возле береговой линии Обского водохранилища на льду над водой (б)
и над мерзлым песком (в)
Действительно, если лед покрывает более низкоскоростную относительно него воду (рис. 1, а-Ь), то эти пики расположены на оси частот с высокой точностью в соответствии с формулой (1). Например, на рис. 2, а частоты первых четырех мод стоячих волн равны 17.4 кГц, 34.8 кГц, 52.1 кГц и 69.5 кГц. В то же время над мерзлым песком (рис. 2, в) распределение пиков согласуется с формулой (2), что говорит о его большей акустической жесткости в сравнении со льдом (в основном, по-видимому, из-за более высокой плотности минеральных зерен).
Соответствие выделенных регулярных пиков стоячим волнам вертикального сжатия-растяжения ледового слоя, а не стоячим волнам других типов (например, сдвиговым), обусловлено применением при измерениях датчика, реги-
стрирующего преимущественно вертикальную компоненту акустического шума. Подтверждение такой природы стоячих волн дает сопоставление результатов прямого измерения толщины льда в контрольной лунке ледомерной рейкой и расчет толщины по формуле (3). Измеренная рейкой толщина льда над водой составила примерно 10.2 см. Определенная по записанному рядом с лункой шуму частота низшей моды равна ^ =19.92 кГц, что при измеренной скорости Ур = 4090 м/с в соответствии с формулой (3) дает практически такую же толщину 10.27 см. Этот пример, кроме того, подтверждает корректность применения рассматриваемого метода для определения толщины ледяного покрова водоемов.
Частоты трех низших мод для всех точек наблюдений двух экспериментов приведены на рис. 3. На этом рисунке хорошо видно, что для наблюдений на р. Зырянка частоты второй и третьей мод превышают частоту первой моды в 2 и 3 раза, соответственно, что согласуется с формулой (1). Такая же закономерность наблюдается для измерений, выполненных на пляже «Звезда», но только для участка профиля 9-15 м. На части профиля от 0 до 8 м частотный интервал между соседними модами равен удвоенной частоте низшей моды, что согласуется с формулой (2). Из этого следует, что от 0 до 8 м профиль проходит над мерзлым песком, от 9 до 15 м - над водой, а между отметками 8 и 9 м находится граница этих зон.
8 12
Расстояние, м
6 9 Расстояние, М
а)
б)
Рис. 3. Частоты первой (квадратные маркеры), второй (круглые маркеры) и третьей (треугольные маркеры) мод стоячих волн сжатия-растяжения для профилей наблюдений на р. Зырянка (а) и на пляже «Звезда» (б)
И, наконец, на рис. 4 показано распределение толщины льда вдоль профилей, на которых были выполнены шумовые измерения. Как можно видеть, на р. Зырянка толщина льда составила в среднем около 11 см, изменения на профиле в основном не превышают 3 см. На Обском водохранилище (возле пляжа «Звезда») толщина льда над водой примерно такая же - порядка 10-11 см. Но в зоне береговой линии наблюдается ее резкое снижение, связанное, по-видимому, с резким уменьшением наклона дна.
а)
б)
Рис. 4. Распределение толщины льда вдоль профилей наблюдений на р. Зырянка (а) и на пляже «Звезда» (б)
Резюмируя, можно отметить, что проведенные эксперименты показали эффективность использования акустических шумов для определения толщины ледяного покрова водоемов, а также для оценки типа подстилающей среды (вода или мерзлый грунт), на которой лежит слой льда.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Кулижников А.М. Георадиолокационные обследования ледовых переправ // Инженерные изыскания. - 2016. - № 8. - С. 22-29.
2. Измерение толщины льда на пресноводном пруде и реке с использованием сигналов ГЛОНАСС и GPS / М.И. Михайлов, К.В. Музалевский, В.Л. Миронов // Современные проблемы дистанционного зондирования Земли из космоса. - 2017. - Т. 14, № 2. - С. 167-174.
3. Спутниковый мониторинг морских льдов / В.Г. Смирнов, А.В. Бушуев, И.А. Бычкова, Н.Ю. Захваткина, В.С. Лощилов // Проблемы Арктики и Антарктики. - 2010. - № 2 (85). -С. 62-76.
4. Судовой телевизионный комплекс - реализация автоматизированной системы натурных измерений толщины морского льда / С.С. Сероветников, С.В. Фролов, А.Э. Клейн // Российская Арктика. - 2018. - №2. - С. 41-55.
5. Радиолокационные аэроледомерные съемки рек, озер, водохранилищ / М.И. Фин-кельштейн, Э.И. Лазарев, А.Н. Чижов - Л.: Гидрометеоиздат, 1984. - 112 с.
6. Хайкин С.Э. Физические основы механики. - М.: Наука, 1971. - 752 с.
7. Double-frequency ground penetrating radar for measurement of ice thickness and water depth in rivers and canals: Development, verification and application / H. Fu, Z. Liu, X. Guo, H. Cui // Cold Regions Science and Technology. - 2018. - Vol. 154. - P. 85-94.
8. Kirby M.S., Hansman R.J., Jr. Method and Apparatus for Measurement of Ice Thickness Employing Ultrasonic Pulse Echo Technique. - United States Patent, Number: 4,628,736, Date of Patent: Dec. 16, 1986.
9. Kolesnikov Yu.I., Fedin K.V. Detecting underground cavities using microtremor data: physical modelling and field experiment // Geophysical Prospecting. - 2018. - Vol. 66, No. 2. -P.342-353.
10. Ice thickness measurements by Raman scattering / S.M. Pershin, V.N. Lednev, V.K. Klinkov, R.N. Yulmetov, A.F. Bunkin // Optics Letters. - 2014. -Vol. 39, No. 9. -P. 2573-2575.
11. Estimation of glacier ice thickness using Ground Penetrating Radar in the Himalayan region / S.K. Singh, B.P. Rathore, I.M. Bahuguna, A.L. Ramnathan // Current Science. - 2012. -Vol. 103, No. 1. - P. 68-73.
12. Shaw A., Declercq N.F. Determination of Ice Thickness on Lakes and Canals Using Audio Sound // Acta Acustica United with Acustica. - 2015. - Vol. 101. - P. 484-493.
13. Arctic sea ice thickness characteristics in winter 2004 and 2007 from submarine sonar transects / P. Wadhams, N. Hughes, J. Rodrigues // Journal of Geophysical Research. - 2011. -Vol. 116. - C00E02.
© К. В. Федин, Ю. И. Колесников, 2019