WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

О ВОЗМОЖНОСТИ ИЗУЧЕНИЯ ВНУТРЕННЕГО СТРОЕНИЯ

КРИСТАЛЛИЧЕСКОГО МАССИВА ГЕОЭЛЕКРИЧЕСКИМИ МЕТОДАМИ

С.П. Левашов1,2, Н.А. Якимчук1,2, И.Н. Корчагин3, Ю.М. Пищаный2

1

Институт прикладных проблем экологии, геофизики и геохимии, Киев, Украина;

2

Центр менеджмента и маркетинга в области наук о Земле при ИГН НАНУ, Киев;

3

Институт геофизики НАН Украины, Киев, Украина

Введение. Компьютеризированные экспресс-технологии геофизических исследований, разработанные на базе геоэлектрических, сейсмоакустических и георадарных методов, позволяют оперативно и эффективно решать широкий круг экологических, инженерно-геологических, геолого-геофизических и гидрогеологических задач. В частности, с помощью технологии [1-3] успешно решались задачи картирования зон загрязнения нефтепродуктами и повышенного увлажнения грунтов, изучения оползневых процессов и провалов, обнаружения и локализации техногенных вытоков из подземных водных коммуникаций, и т.д. Ниже представлены результаты применения экспресстехнологии геоэлектрических исследований, включающей методы становления короткоимпульсного поля (СКИП) и вертикального электрорезонансного зондирования (ВЭРЗ) [1-3], для изучения внутренней структуры кристаллического массива.

Общие сведения. В 2004-2005 г.г., в рамках проекта УНТЦ № «Обоснование возможности захоронения радиоактивных отходов в восточной части Коростенского кристаллического массива (полевые и модельные исследования)” геоэлектрические работы с использованием технологии СКИПВЭРЗ проведены на Вереснянской площади.

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

Экспресс-технология СКИП-ВЭРЗ [1-3] базируются на изучении геоэлектрических параметров среды в импульсных неустановившихся геоэлектрических полях, а также квазистационарного электрического поля Земли и его спектральных характеристик над исследуемыми объектами.

Методом СКИП регистрируется процесс становления поля короткого электрического импульса в малогабаритных дипольных ферритовых антеннах.

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

Зондирование ВЭРЗ основано на изучении процессов естественной поляризации среды и спектральных характеристик естественного электрического поля над объектом (залежью). Для горизонтально-слоистых разрезов эта составляющая технологии дает возможность эффективно разделять разрез в точке зондирования на отдельные стратиграфические подразделы и с высокой точностью определять глубину их залегания. Эффективность зондирования значительно повышается при наличии в районах исследований задокументированных параметрических буровых скважин, которые предоставляют возможность соответствующим образом “калибровать” диаграмму зондирования – связать отдельные интервалы значений на диаграмме с соответствующими стратиграфическими разновидностями разреза.

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

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

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

Методика полевых работ. При проведении съёмки методом СКИП в автомобильном варианте скорость движения автомобиля составляет 10- км/час. Профили съёмки прокладывались по дорогам в пределах Вереснянского участка. В процессе съемки, каждые 0.5 секунд автоматически производился опрос GPS-приёмника и при удалении от предыдущей точки съёмки более чем 50 м компьютер запускал генератор электромагнитных импульсов. С приемной антенны в цифровом виде осуществлялась регистрация процесса становления электромагнитного поля станцией сбора и накопления геофизической информации “GEMA”. При проведении съемки в пешеходном варианте момент генерации электромагнитных импульсов задавался вручную, через каждые 10м по профилю. Синхронно с регистрацией процесса становления выполнялась первичная обработка данных, определялись точки, принадлежащие аномальным зонам.



Для построения вертикальных геоэлектрических разрезов и карты кровли кристаллических пород выполнено зондирование ВЭРЗ по двум взаимно перпендикулярным профилям субширотного и субмеридионального направления. Пункты зондирования располагались через 500 м.

Результаты работ. Геоэлектрическая съёмка методом СКИП. По данным съемки построена карта аномальных значений СКИП для Вереснянского участка (рис. 1). В центральной части участка исследований выделяется зона относительно слабых аномальных значений. В зоне минимальных значений поля СКИП выделяется наиболее стабильный блок кристаллических пород.

Приблизительные размеры этого блока – 42 км.

Рис. 1. Карта геоэлектрических аномалий СКИП интенсивности аномальных значений; 2 – кристаллического фундамента граница участка; 3 – линии сейсмических (мощности осадочного чехла).

профилей; 4 – линии геоэлектрических разрезов Наименьшая глубина залегания по данным ВЭРЗ; 5 – зоны тектонических кровли кристаллических пород нарушений, выделенные по геоэлектрическим выявлена в центральной части данным; 6 – пункты съёмки СКИП; 7 – точки ВЭРЗ V0-V26; 8 – точки ВЭРЗ для определения участка – 165 м. В северной и глубины залегания фундамента.

величина мощности осадочного комплекса равна 180-185 м.

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

Результаты зондирования ВЭРЗ вдоль профилей № 1-1а, 2-2а. Вдоль линий профилей № 1 и 2 (рис. 1) с шагом 500 м проведено вертикальное зондирование в интервале глубин 0-5000 м в 27 пунктах. Шаг зондирования по глубине равнялся 2 м. По данным зондирования построены диаграммы ВЭРЗ, на которых выделены участки аномально поляризованных зон положительной и отрицательной поляризации АПЗ(+) и АПЗ(-). Зоны отрицательной поляризации связаны с участками относительного понижения геоэлектрического сопротивления, положительной поляризации – с участками повышенного сопротивления. Зоны нарушений пород фундамента выделялись по диаграммам зондирования на тех интервалах глубин, где амплитуда резонансной кривой превышает 40 mV.

Для пунктов зондирования № 4, 7, 8, 9, 22, 13, 17, 18 АПЗ практически не выявлены. Эти пункты расположены над ненарушенными блоками кристаллического массива. Для остальных точек зондирования на различных интервалах глубин выделены АПЗ, которые связаны с нарушенными зонами фундамента.

Построена также карта суммарных значений интенсивности абсолютных величин АПЗ. Минимальные значения суммарной интенсивности АПЗ фиксируются над ненарушенными блоками кристаллического фундамента.

Такая зона выделяется в центральной части Вереснянского участка. По предварительной оценке горизонтальные размеры блока – 42 км.

Вертикальные разрезы по данным ВЭРЗ. По данным зондирования построены вертикальные геоэлектрические разрезы № 1-1а и № 2-2а (рис. 1, 2).

На профиле №1 (рис. 2) массивный (ненарушенный) блок кристаллических пород выделен в интервале между пунктами зондирования № 23 и № 6, Рис. 2. Вертикальный разрез по профилю I-Ia через Вереснянскую площадь. 1 – осадочный чехол; 2 – ненарушенные зоны фундамента; 3 – зоны тектонических нарушений; 4 – точки зондирований ВЭРЗ № 0-9, 22-26; 5 – трещиноватые зоны фундамента.

определены на глубинах 2000м и 2500м. В северо-восточной части нарушенные зоны выделены на глубинах 1000, 2000, 2500 м и на интервале от 3000 до На профиле № 2 массивный блок выделен в интервале между пунктами зондирования № 11-№ 13 (протяженность блока 1 км, интервал глубин от до 2000 м), второй блок – между пунктами № 17 – № 18 (протяженность блока 0.5 км, интервал глубин 300-3000 км.) В юго-восточной части участка выделены зоны нарушений, расположенные в интервалах глубин 500-1000 м и 3000-4000 м.

Выводы. По результатам геоэлектрических исследований на Вереснянской площади Коростеньского плутона выделен и частично оконтурен блок ненарушенных пород кристаллического фундамента с горизонтальными размерами 42 км. Блок расположен в центральной части участка. Глубина кровли массивного блока варьирует от 300 до 500 м. На уровне поверхности кристаллического фундамента блок проявляется в виде небольшого поднятия с перепадом высот до 30 м. Подошва блока зафиксирована в интервале глубин от 2000 до 4000 м. С северной части блок ограничен субширотным тектоническим нарушением, с запада и востока – зонами дробления.

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

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

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

Комплексирование технологии СКИП-ВЭРЗ с гравиметрическим и магнитометрическим методами (и другими!) исследований будет способствовать существенному повышению эффективности геофизических методов при решении экологических, инженерно-геологических и геолого-геофизических задач в целом.

1. Левашов С.П., Якимчук Н.А., Корчагин И.Н. Электрорезонансное зондирование и его использование для решения задач экологии и инженерной геологии // Геологический журнал. – 2003. – № 4. – С. 24-28.

2. Левашов С.П., Самсонов А.И., Якимчук Н.А., Корчагин И.Н., Таскинбаев К.М.

Использование геоэлектрических методов при проведении рекогносцировочных исследований на нефть в Западном Казахстане // Геоинформатика. – 2004. – № 1. – С. 21-31.

3. Левашов С.П., Якимчук Н.А.., Пищаный Ю.М., Корчагин И.Н. Аппаратурный комплекс «GEMA» комплексных геоэлектрических исследований и его использование для поисков скоплений углеводородов. Сборник научных трудов НГА Украины № 13, Том 4. – Днепропетровск: РИК НГА. – 2002. – С. 78-83.

ИЗУЧЕНИЕ ПЛОЩАДЕЙ РАСПРОСТРАНЕНИЯ МИНЕРАЛЬНЫХ ВОД

ГЕОЭЛЕКТРИЧЕСКИМИ МЕТОДАМИ

С.П. Левашов1,2, Н.А. Якимчук1,2, И.Н. Корчагин3, Ю.М. Пищаный Институт прикладных проблем экологии, геофизики и геохимии, Киев, Украина;

Центр менеджмента и маркетинга в области наук о Земле при ИГН НАНУ, Киев;

Введение. На протяжении 2001-2005 г.г. экспресс-технология геоэлектрических исследований СКИП-ВЭРЗ широко применялась для решения различных экологических, инженерно-геологических и геологогеофизических задач. В частности, технология применялась для поисков и картирования водо-насыщенных коллекторов [1-3]. Ниже представлены результаты ее применения для картирования зон распространения минеральной воды в районе города-курорта Моршин в Западной Украине.

Общие сведения. В сентябре 2005 г. в пределах Моршинского месторождения минеральных вод выполнены геофизические исследования на трех участках. Участок № 1 (3.0 х 1.0 км) расположен в лесной зоне на правом берегу речки Бережница севернее г. Моршин. В пределах участка расположены источники минеральной воды № 1, № 2, № 3, № 4. Участок № 2 (1.5 х 0.8 км) расположен западнее г. Моршин, в его пределах имеется буровая скважина, в которой установлено наличие минеральной воды. Участок № 3 (1.2 х 0.5 км) расположен на южнее с. Лисовичи, на правом берегу речки Сукель.

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

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

Во время работ проводились геоэлектрические и сейсмоакустические исследования методами становления короткоимпульсного поля (СКИП), вертикального электрорезонансного зондирования (ВЭРЗ) и сейсмоакустического зондирования. Съемка СКИП использовалась для построения карт распространения минеральной воды. Методы зондирования применялись для определения глубин залегания водонасыщенных горизонтов и построения вертикальных колонок и разрезов.

Методика полевых работ. Полевая съемка методом СКИП проводились вдоль отдельных профилей с шагом 20-50 м. При выявлении аномальной зоны повышенной увлажненности грунтов по одному из профилей, шаг измерений уменьшался, а число профилей увеличивалось.

Через аномальные зоны, выявленные методом СКИП, проложены профили ВЭРЗ и сейсмоакустического зондирования. Сейсмоакустические исследования выполнялись с шагом 5 м. Глубина исследования ВЭРЗ и сейсмозондирования составляла 40 м. В отдельных точках ВЭРЗ выполнялось до 500 м. Координаты точек начала и конца профилей съемки и точек ВЭРЗ определялись при помощи GPS-приемника.

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

На участке № 1 аномальная зона повышенного содержания воды (рис. 1) залегает вдоль правого берега р. Бережница в виде полосы шириной от 100 до 250 м. Рабочие каптажи № 1, 2, 3 расположены в центральной части зоны, Рис. 1. Карта зоны распространения минеральной воды сильного увлажнения на участке № 1 месторождения. 1 – границы участка ; грунта. Для строительства 2 – тектонические нарушения; 3 – направление миграции каптажей необходимо подземных водных потоков; 4 – линии вертикальных выбирать места, которые разрезов; 5 – пункты съемки СКИП; 6 – станции ВЭРЗ; расположены на левых 7 – каптажи источников № 1-4; 8 – естественный источник. сторонах оврагов и балок, В северо-восточной части участка (точка К8) выявлен интенсивный подземный источник воды.

Поперек зоны повышенного обводнения пород проложены профили ВЭРЗ, вдоль которых построены схематические геологические разрезы (№ 1-8).

Профили проходят через зоны каптированных источников № 1, 2, 3, 4, а также через места, рекомендованные для сооружения новых каптажей. Во всех точках зондирования минеральная вода фиксируется в гравийно-галечных отложениях.

Профиль 1-1а, источник № 1 (К1). Каптаж расположен вблизи максимума аномальной зоны. Ширина зоны повышенного увлажнения – 150 м. Данные зондирования в районе каптажа К1: суглинки – 0-2 м; галечник – 2-9 м;

минеральная вода – 4-9 м. Зоны повышенной увлажненности в нарушенных песчано-глинистых отложениях установлены в интервалах 60-64 м и 449-463 м.

Профиль 2-2а, источник № 2 (К2). Ширина зоны обводнения – 250 м.

Каптаж расположен на левой стороне балки, в центре аномалии. Профиль зондирования проложен по дну балки. Данные зондирования в районе каптажа К2: суглинки – 0-2 м; галечник – 2-11 м; поверхностная вода – 2-3 м;

минеральная – 3-11 м (8 м). Зона повышенной увлажненности в нарушенных песчано-глинистых отложениях установлена в интервале 11-462 м.

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

Профиль 3-3а проложен по дну балки. Ширина зоны обводнения – 150 м.

Точки ВЄРЗ № 4 и № 5 размещены в центральной части зоны аномалии.

Данные ВЭРЗ для точки № 4: суглинки – 0-2.5 м; галечник – 2.5-10 м;

поверхностная вода 2.5-4.5 м; минеральная – 4.5-10 м. Зоны увлажнения в песчано-глинистых отложениях установлены в интервалах 62-64 м и 449-461 м.

Профиль 4-4а привязан к точке зондирования № 6. Ширина зоны обводнение – 180 м. Данные ВЭРЗ для точки № 6: суглинки – 0-3 м; галечник 3-12 м; поверхностная вода – 2-2.5 м; минеральная вода – 2.5-12 м.

Профиль 6-6а. Источники № 3, 4. Ширина аномальной зоны – 70 м.

Данные зондирования в районе каптажа К3: суглинки – 0-4 м; галечник – 4- м; поверхностная вода – 4-4.5 м; минеральная – 4.5-15 м. Данные ВЭРЗ в районе каптажа К4: суглинки – 0-4 г; галечник – 4-15 м; поверхностная вода – 4-7 м; минеральная вода – 8-15 м. Для нормального функционирования каптажа № 4 отбор воды должен быть оборудован на глубине 8 м и ниже.

Профиль 7-7а проложен по склону горы, вне овражной части. Ширина аномальной зоны – 80 м. Данные ВЭРЗ для точки № 9: суглинки – 0-10.5 м;

галечник – 10.5-16.5 м; поверхностная вода – 7.5-10.5 м; минеральная вода – 10.5-16.5 м. Зоны повышенной увлажненности в нарушенных песочноглинистых отложениях установлены в интервалах 60-66 м и 458-483 м.

Профиль 8-8а проложен вдоль ручья, который вытекает из интенсивного источника. Ширина аномалии – 130 м. Данные ВЭРЗ для точки № 9: суглинки – 0-1 м; галечник – 1-11 м; поверхностная вода – 0-1 м; минеральная вода 1-11 м.

Участок № 2 пересекают две балки. В борту одной из них пробурена скважина, которая вскрыла пласт с минеральной водой. Съемкой СКИП на площади выделена аномальная зона повышенного водосодержания (рис. 2), которая протягивается с юго-запада на северо-восток. Миграция подземных вод происходит в северо-восточном направлении вдоль тектонически ослабленной зоны. Эрозионные врезы балок пересекают зону, вскрывая гравийно-галечный горизонт, из которого часть воды поступает в ручьи. Через аномалию построены два вертикальных разреза по профилям № 1-1а и № 2-2а.

Результаты ВЭРЗ в районе скважины. Основные горизонты по данным ВЭРЗ: суглинки – 0-9 м; галечник – 9-21 м; поверхностная вода – 7-15 м;

Рис. 2. Карта зоны распространения минеральной воды вдоль профиля – 600 м.

на участке № 2 месторождения. 1 – границы участка Точки ВЭРЗ 11 и исследований; 2 – тектонические нарушения; 3 – размещены вблизи профили вертикальных разрезов 1-1а, 2-2а по данным максимума аномалии. В ВЭРЗ; 4 – направление миграции подземных водных этом месте наибольшая потоков; 5 – пункты съемки СКИП; 6 – станции ВЭРЗ; – контур зоны распространение минеральной воды; 8 – рекомендованные участки для заложения буровых горизонта: суглинки – 0- скважин на воду; 9 – направление миграции подземных м; галечник – 13- 25 м;

водных потоков.

13-25 м (12 м).

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

Участок № 3 расположен вдоль правого берега г. Сукель. Аномальная зона состоит здесь из двух расходящихся рукавов. Основной участок зоны простирается в северном направлении вдоль тектонически ослабленной зоны. В этом же направления происходит миграция воды подземных водных потоков.

Вдоль берега речки по профилю 1-1а проведено сейсмоакустическое и электрорезонансное зондирование. В точках ВЭРЗ № 114 и № зарегистрированы локальные экстремумы аномалии. Геологические границы для этих точек: ВЭРЗ № 114 – суглинки – 0-4.5 м; галечник – 4.5-11 м; поверхностная вода – 3-4.5 м; минеральная – 4.5-11 м; ВЭРЗ № 84 – суглинки – 0-4.5 м; галечник – 4.5-11 м; поверхностная вода – 4-7.5 м; минеральная – 7.5-11 м (3.5 м). Район точек ВЭРЗ № 114 и № 84 является оптимальным для заложения добывающих буровых скважин.

Общие результаты геофизических работ. В результате проведенных исследований на трех участках установлено, что Моршинськие минеральные воды залегают в основном вдоль тектонически нарушенных зон. Гравийногалечные отложения в верхней части разреза являются своеобразным накопительным резервуаром, в который минеральные воды поступают за счет вертикальной миграции из нижних горизонтов. Участки вертикальной миграции установлены в районе каптажа № 2 (К2, участок 1) и в центральной части участка № 2. Вертикальная миграция происходит в зонах максимального дробления пород и пересечения тектонических нарушений.

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

Выделены три основные нарушения, вдоль которых происходит миграция минеральной воды. Особый интерес представляет субмеридиональное нарушение, которое проходит восточнее от г. Моршин, через с. Лисовичи, лесную зону и через зону источника № 2 на участке № 1. Вдоль линии тектонического нарушения образовался целый ряд балок. Эрозионные врезы балок пересекают линию тектонически нарушенной зоны и могут вскрывать гравийно-галечный горизонт, по которому осуществляется миграция минеральных вод. Данная зона представляет интерес для проведения работ по поиску и разведке источников минеральной воды в районе г. Моршин.

Выводы. По данным геофизических исследований на участках № 1- выделены и закартированы зоны распространения минеральных вод.

Определены места максимального водонасыщения отложений гравийногалечного горизонта. Установлено, что основные залежи минеральной воды расположены вдоль зон тектонических нарушений. Гравийно-галечный пласт является накопительным резервуаром минеральных вод, которые поступают с глубин 450-460 м. В пределах участков № 1 и № 2 определено несколько мест вертикальной миграции воды из нижних горизонтов.

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

1. Левашов С.П., Якимчук Н.А., Корчагин И.Н. Использование геоэлектрических методов для поиска и картирования водоносных коллекторов / Проблемы водных ресурсов, геотермии и геоэкологии. Материалы Международной научной конференции, посвященной 100-летию со дня рождения академика Г.В.

Богомолова, Минск 1 – 3 июня 2005 г. В двух томах. Том 1. – С. 137-139.

2. Levashov S.P., Yakymchuk N.A., Korchagin I.N., Pyschaniy Ju.M., Yakymchuk Ju.N. The searching and mapping the water-saturated rock by geoelectric methods. 67nd EAGE Conference and Technical Exhibition. Madrid, Spain, 13 – 16 June 2005. CD-ROM Abstracts volume, 4 pages.

3. Levashov S.P., Yakymchuk N.A., Korchagin I.N., Pyschaniy Ju.M., Yakymchuk Ju.N. The water-saturated rock searching and mapping by geoelectric methods. The 4th Balkan Geophysical Congress – International Conference of Applied Geophysics and Earth Physics, Bucharest Romania, 9-12 Oct. 2005. Extended Abstracts compact disk. 4 pages.

О ВОЗМОЖНОСТИ РАЗДЕЛЕНИЯ ПО ГЛУБИНЕ ИСТОЧНИКОВ

МАГНИТНОГО ПОЛЯ

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

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

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

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

При пересчете вверх остаточного поля использовалась формула Пуассона, дающая решение задачи Дирихле для полупространства:

решения плоской задачи Дирихле на дневной поверхности, левая часть найденные значения поля на плоскости z = H в верхнем полупространстве.

Чтобы окончательно избавиться от источников в горизонтальном слое от дневной поверхности до глубины H, пересчитанное вверх поле было затем продолжено вниз на глубину H. При этом соотношение (1) используется в качестве интегрального уравнения, в котором теперь функция U (x, y, H ) задана, в ядре H заменяется величиной 2H, а искомым решением является функция U ( x, y,- H ). Решение уравнения 1-го рода (1) представляет собой некорректную задачу, что делает необходимым использование одного из методов регуляризации. Оператор уравнения (1) является положительно определенным и самосопряженным, поэтому можно использовать схему М. М. Лаврентьева [1] и вместо уравнения Kw = U, которое является операторной записью соотношения (1), решать регуляризованное уравнение (K + aE )w = U при подходящем выборе параметра регуляризации. Особенности найденной функции U ( x, y,- H ) лежат глубже H, выше этого уровня она является гармонической, в связи с чем можно трактовать ее как поле глубинных источников. Сумму этого поля, пересчитанного на дневную поверхность, и предварительно вычтенного решения задачи Дирихле мы использовали как поле боковых и глубинных источников. Разность измеренных значений гравитационного поля и полученного поля сторонних источников использовалась нами в качестве гравитационного эффекта от источников, лежащих в горизонтальном слое от дневной поверхности до глубины H.

Обсудим проблемы, связанные с применением методики для обработки магнитных данных. Как известно, поле точечного источника в магнитном случае имеет более сложный характер, чем в гравитационном: при удалении от максимума поле Z переходит через 0, а затем медленно возрастает от отрицательных значений к нулю. Возможно наложение на минимум максимума от соседнего источника. Обычно измеренное магнитное поле является гораздо более градиентным, чем гравитационное. Поле глубинных источников часто практически неразличимо из-за большого количества аномалий от приповерхностных локальных магнитных объектов. С другой стороны, магнитный потенциал и его производные являются гармоническими в верхнем полупространстве, поэтому при пересчете магнитного поля вверх также может использоваться формула Пуассона (1).

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

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

Применение методики привело к целому ряду любопытных результатов.

При малых значениях H магнитные и гравитационные поля от источников в полупространстве ниже глубины H значительно отличались. По-видимому, это позволяет сделать вывод о различной природе приповерхностных источников гравитационных и магнитных аномалий. При использовании значения H = 5 км обнаружилось, что поле от глубинных объектов, как в гравитационном, так и в магнитном случае имеет одинаковый характер: наблюдаются два максимума, положение которых качественно совпадает (см. рис. 1 и рис. 2 –значения в мгл и нТл, соответственно). Возможно, поля обусловлены глубинной границей раздела, на которой наблюдается скачок как плотности, так и магнитной восприимчивости.

Весьма интересным представляется последовательное устранение источников при увеличении значения H. Если источник присутствует в поле от полупространства ниже глубины H1 и практически исчезает при использовании значения H2, где H2 > H1, то можно предполагать, что он локализован в горизонтальном слое между глубинами H1 и H2. Это обстоятельство позволяет разделить по глубине источники поля.

1. Мартышко П.С., Пруткин И.Л. Технология разделения источников гравитационного поля по глубине // Геофизический журнал,т. 25, № 3.2003. С.159-168.

ГРАВИМЕТРИЧЕСКОЕ ИЗОБРАЖЕНИЕ СОЛЯНОКУПОЛЬНОГО

СТРУКТУРНОГО ЭТАЖА ПРИКАСПИЙСКОЙ ВПАДИНЫ

РГП «Специализированное гравиметрическое предприятие», Алматы, Казахстан Еще на заре применения трансформаций гравитационного поля для его геологической интерпретации действие площадного осреднения в пределах круга с радиусом R считалось равносильным «раскатыванию» аномальных масс в плоский диск, радиусом и высотой, равными R (А.Н. Тихонов, Ю.Д. Буланже, 1945 г.). Это на первый взгляд вполне логичное утверждение строго не доказано, более того нетрудно показать, что для многих моделей оно вообще ошибочно. Однако для некоторых геолого-геофизических условий, по крайней мере, в первом приближении оно может выполняться. Для этого класса моделей будем считать, что в результате осреднения аномалии, связанные с неоднородностями до глубины R в значительной степени «подавлены», то есть амплитуды их снижены. Было установлено, что эффект подавления аномалий («раскатывания» масс) увеличивается, при многократном осреднении поля в скользящем окне.

Выполним n последовательных осреднений гравитационного поля с увеличением размеров скользящего окна на постоянную величину 2·dH, начиная с 2·dH, каждый раз сглаживая поле Ui, полученное в предыдущем цикле осреднения. В результативной серии составляющих гравитационного поля будут подавлены аномалии, обусловленные неоднородностями верхней толщи разреза мощностью i·.dH, где i – порядковый номер осреднений. В разностных полях dUI=UI–UI+1 подавлены аномалии от неоднородностей выше i·dH и глубже (i+1)·dH и, как следствие, преобладают аномалии от (i+1)-ого слоя мощностью dH. Не трудно заметить, что в результате выполнена полосовая фильтрация исходного поля. Если каждую из вычисленных трансформант отнести к горизонтальной плоскости на уровне глубин i·dH, то полученный набор числовых матриц можно рассматривать, как объемную числовую модель гравитационного поля.

На рисунке 1а показан разрез объемного гравитационного поля по профилю, расположенному в центральной части Прикаспийской впадины (dH=0.5 км). На разрезе четко выделяется зона чередующихся интенсивных минимумов и максимумов, которые в нижней части разреза переходят в относительно спокойное поле. Это однозначно соответствует разделению геологического разреза на солянокупольный структурно-тектонический этаж и подсолевой комплекс отложений, поэтому рассматриваемые локальные аномалии назовем аномалиями солянокупольного типа. На разрезе вынесены данные сейсморазведки - подсолевой горизонт П1, и кровля соли в сводах куполов (на склонах и в межкупольных зонах отражения от кровли соли прослеживаются неуверенно). Минимумы хорошо согласуются с соляными куполами, однако глубина до кровли соли в сводах куполов не всегда удовлетворительно совпадает с верхней границей отрицательного поля.

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

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

Представление результатов фильтрации в виде увязанных по глубине вертикальных разрезов, облегчающее геологическую интерпретацию гравитационного поля, правомерно назвать гравиметрическим изображением геологической среды. Интересно, что такой подход вполне приемлем и для магнитного поля. В связи с этим рассмотрим один пример. Аномалии магнитного поля Прикаспийской впадины в основном связаны с неоднородностями кристаллического фундамента. По мнению многих исследователей, они обусловлены интрузивными телами, которые частично могли внедряться в нижние части разреза осадочного чехла. Наиболее глубокий горизонт, из прослеживаемых сейсморазведкой в осадочных отложениях Прикаспийской впадины, стратифицируется как девонский (П3). В некоторых частях региона горизонт П3 приближается и сливается с поверхностью фундамента. Учитывая это, результаты объемного преобразования магнитного поля в первом приближении можно увязывать с горизонтом П3. На рисунке показано грави-магнитное изображение геологической среды по одному из разрезов в южной части Прикаспийской впадины. Верхняя часть разреза, между подошвой юрских отложений (горизонт V) и подсолевым ложем (П1), совмещена с объемным гравитационным полем, которое отображает строение солянокупольного структурного этажа. В низах разреза (ниже горизонта П3, включая фундамент) показано объемное магнитное поле, полученное в результате преобразований данных аэромагнитной съемки по методике, разработанной нами для гравитационного поля. Объем статьи не позволяет углубляться в геологическую интерпретацию приведенного разреза. Поэтому здесь лишь показываем удобство интерпретации гравимагнитных материалов, представленных в виде разрезов полей, увязанных по глубине.

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

ОТОБРАЖЕНИЕ СОЛЯНЫХ КУПОЛОВ ПРИКАСПИЙСКОЙ ВПАДИНЫ

В МАГНИТНОМ ПОЛЕ И ПЕРСПЕКТИВЫ ЕГО ИСПОЛЬЗОВАНИЯ ДЛЯ

ИЗУЧЕНИЯ СОЛЯНОКУПОЛЬНОГО СТРУКТУРНОГО ЭТАЖА

РГП «Специализированное гравиметрическое предприятие», Алматы, Казахстан гравитационном поле интенсивными минимумами общеизвестно. Вот уже три четверти века гравиразведка используется для изучения строения солянокупольных структур. Аномалии магнитного поля связаны, в основном, с неоднородностями кристаллического фундамента, которые отображаются аномалиями интенсивностью десятки - сотни нТл. По единодушному мнению многих исследователей эти аномалии обусловлены интрузивными телами в фундаменте, которые могли частично внедряться в нижние части разреза осадочного чехла. В восточной и северо-восточной частях Прикаспийской впадины в результате детального анализа материалов аэромагнитных съемок были выделены аномалии в 2 - 10 нТл, контролирующие краевые части соляных куполов (В.А. Боронаев). Для этого выполнялась высокочастотная фильтрация поля в скользящем окне по маршрутам аэромагнитных съемок, что позволило избежать влияния погрешностей их взаимной увязки. Обязательным условием выделения таких аномалий было близкое к ортогональному пересечение границ куполов аэромаршрутами. При уменьшении углов пересечения до 50О и менее контуры куполов не выделялись.

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

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

Отображение соляных куполов в магнитном поле не столь очевидно, как проявление их в гравитационном поле. Сопоставление локальных минимумов магнитного поля и контуров солянокупольных структур в обобщенном виде показано на схеме (рис. 1). Из рисунка видно насколько неоднозначно соотношение между сводами куполов и аномалиями магнитного поля. Так на площади крупных соляных массивов Жекендысай, Сауырбай, Егизкара локальные аномалии не выделяются, с узким соляным телом купола Сауле достаточно точно совпадает контур минимума магнитного поля. На многих куполах наблюдается частичное перекрытие магнитных минимумов и сводов поднятий соли. Для объективной оценки с помощью средств системы СПОГМ выполнена статистическая оценка совпадений локальных минимумов со сводами куполов. Так как из-за неравномерной густоты сейсмических профилей своды куполов не всегда полностью оконтурены, то это было восполнено сопоставлением контуров минимумов магнитного и гравитационного полей. В результате установлено, что 68% площади сводов куполов совпадет с локальными минимумами магнитного поля. Примерно на столько же совпадает площадь гравитационных минимумов с магнитными (70%). Совпадение площади сводов куполов с гравитационными минимумами более значительное - 86%. В то же время только 44% площади магнитных минимумов совпадает со сводами куполов и 62% с гравитационными минимумами. Это объясняется тем, что общая площадь отрицательных аномалий магнитного поля превышает площадь сводов куполов на 50%, и больше площади гравитационных минимумов на 12%.

Из приведенных статистических данных ясно, что взаимосвязь между локальным магнитным полем и соляными куполами весьма сложная. Создается впечатление о невысокой перспективе использования данных магниторазведки для изучения солянокупольных структур. Более детальный анализ позволил нам сделать более оптимистические выводы. Мы уже отмечали, что крупные соляные массивы в магнитном поле не проявляются. Исключение их может существенно повысить статистические показатели отображения куполов в магнитном поле. В то же время узкие протяженные соляные гряды в магнитном поле выражаются довольно четко. Более того, в магнитном поле проявляются небольшие изометричные обособленные соляные штоки, расположенные во внутренних частях крупных межкупольных зон (Каратау, Бекет, Кзылжарский).

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

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

Отображение группы соляных куполов в магнитном и гравитационном полях а) магнитное поле; б) соляные купола (данные сейсморазведки);

в) локальные минимумы магнитного поля (площадная узкополосная фильтрация);

г)локальные минимумы гравитационного поля; (в и г на фоне контуров соляных куполов) Ярким примером отображения соляных куполов в магнитном поле может служить группа соляных куполов Тобекудук – Сауле (рис.2). Здесь локальные минимумы магнитного поля образуют две кольцевые зоны, с высокой точностью совпадающие с соляными грядами и локальными гравитационными минимумами. Даже небольшие своды соли отображаются иногда более контрастно, чем в гравитационном поле (например, свод на перешейке Тобекудук - Тобекудук Сев., купол Кособа Зап.).

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

1) Граница купола (сейсморазведка) совпадает с контуром минимума гравитационного и магнитного полей. 2) граница купола, совпадая с контуром магнитного минимума, выходит за пределы гравитационного минимума. 3) Небольшой обособленный купол, зафиксированный сейсморазведкой, четко отображен в магнитном поле, но слабо выражен в гравитационном. 4) На площади совпадающих, небольших по площади минимумов магнитного и гравитационного полей купол сейсморазведкой не оконтурен. 5) Купол в гравитационном поле отображен интенсивным гравитационным минимумом но локального минимума магнитного поля выделить не удается. Есть основания ожидать, что глубокий анализ совокупности этих данных может привести к важным геологическим выводам. Например, вполне вероятно, что ситуация связана с наличием соляного карниза, а в случае 4 купол сейсморазведкой пропущен (редкая сеть профилей или плохое качество материала).

УСТОЙЧИВОСТЬ МАК-ВОЛН В ЗЕМНОМ ЯДРЕ ЗЕМЛИ ДЛЯ

МОДЕЛИ СЖИМАЕМОЙ СРЕДЫ.

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

Характерные величины периодов колебаний близки к значениям: Т~360, 600, 900, 1200, 180, 2700, 3600, 5400 и 9000 лет. Согласно [Брагинский, 1970] средние частоты основного спектра, порядка тысячи лет, представляют собой распространяющиеся в азимутальном направлении МАК-волны, которые отражают принципиальную неустойчивость процесса генерации. Устойчивость МАК-волн обеспечивается балансом между магнитными, архимедовыми и кориолисовыми силами. Именно с наличием этих волн, которые возникают за счет действия архимедовых сил, связано появление асимметричных скоростей, необходимых для генерации полоидального поля в динамо.

В развитие результатов [Брагинский, 1967; Braginsky, 1980] рассмотрена модель МАК-волн с учетом сжимаемости вещества внешнего ядра Земли.

Система уравнений магнитной гидродинамики сжимаемой жидкости в отсутствии диссипации магнитного поля, в пренебрежении, меридиональными составляющими магнитной индукции и скорости для малых возмущений v ', B ', C ', p ' аналогично [Брагинский, 1967], может быть записано в следующем виде:

где P = ( p ' / r0 ) + Bj Bj' / m0 r0, r0 – средняя плотность во внешнем ядре Земли, r1 = r1 ( s, z ) – безразмерная функция, характеризующая равновесное меридиональное распределение относительное превышение локальной плотности r над гидростатической r0 r1.

В уравнении (1.1) отброшены инерциальные члены, а в (1.5) диссипативные члены.

где z – локальная угловая скорость. Далее для переменных v’ и В’ аналогично [Брагинский, 1967 ] вводятся вектора w и y таким образом, что справедливо divw = divy =0 и с использованием (1.2) и (1.3) имеем:

Решения для этих переменных ищутся в виде:

Соотношение для w и y выглядит следующим образом [Миндубаев, 2003]:

2. Для модели течения в торе с бесконечным радиусом аналитически получено условие устойчивости для распространения МАК-волн. В явном виде получена зависимость частот волн от основных физических параметров.

Показано, что частота МАК-волн определяется азимутальной скоростью в тех областях, где нарушается условие устойчивости. В данной модели это соответствует областям внешнего ядра Земли с пониженным значением азимутального магнитного поля. Качественно показано, что для вариаций, обладающих свойством кратности частот, характерна одна и та же скорость дрейфа, что соответствует наблюденным данным. Наличие по наблюденным данным в спектре вариаций наряду с характерным западным дрейфом восточного дрейфа, объясняется знакопеременностью геострофической скорости. Для западного направления дрейфа – отрицательное значение скорости, соответственно для восточного – положительное. Картина распространения МАК-волн в нашей модели с учетом сжимаемости вещества внешнего ядра Земли приводит локализации области возмущений у внешней границы.

Получены численные результаты распространения МАК-волн в модели течения в торе с конечным радиусом при различных распределениях магнитного поля. Рассматривается течение проводящей сжимаемой жидкости в торе прямоугольного сечения, вращающемся вокруг оси z. –z1 < z < z1, 0 < x < L, где x = (s - 0.35)/0.65, L=R1 – R2; R1 и R2 = 0.35R1 - соответственно размеры внешнего и внутреннего радиусов тора. Уравнения для малых возмущений магнитного поля, скорости, концентрации и давления описываются системой (1). Векторы g и Cp принимаются параллельными и направленными вдоль оси s.

Для распределения азимутального магнитного поля примем Bj ( x) = Bj 0 B1 ( x).

Решения уравнения ищутся в виде:

где k z = nzp / 2 z1, nz - целое. При принятых условиях задача сводится к одномерной и решается линейная система на собственные значения l и собственные функции yx(x) и yz(x):

Линейные вещественные операторы Lik имеют вид:

где за единицу длины принято R1, l = 2Ww R12 / ca 0 m, u ( x) = 2Wz 1R12 / ca 0, B = B12 ( x) / r1.

Среднее значение квадрата альфвеновской скорости в области и плотности принимаются соответственно gx(x) = gx0g1(x), Cx(x) = C0C1(x), x(x) = 01(x).

Граничные условия задаются в виде:

Для численного решения системы (2) в области 0 x 1 использовалась однородная пространственная сетка с постоянным шагом h = xi+1 – xi = 1/36, производные заменялись центральными разностями с квадратичной аппроксимацией. Полученная в результате система на собственные значения и собственные функции решалась с использованием математической библиотеки IMSL Fortran 90 MP.

Для безразмерных распределений ускорения силы тяжести, градиента концентрации, плотности и геострофической скорости принято:

1(x) = (1.1-0.2x), u(x) = u0(1 – x2)((0.51)1/2 – (1 – x2)1/2/(1 – (0.51)1/2), где D – коэффициент, u0 = 100.

Влияние архимедовой силы проявляется, если она имеет поперечную составляющую z. Для этого параллельные вектора g и Cp повернуты на малый положительный угол /36. Тогда с допустимой точностью в рассматриваемой задаче, можно принять xn x и z x (/36). Рассчитаны две модели с учетом z для распределений магнитного поля B1(x):

Полю В1 соответствует распределение магнитного поля с непроводящими границами на R1 и R2, а В2 с проводящей границей на R2. Значения коэффициентов выбраны из условия нормировки. Собственные значения i неустойчивых мод для распределения В1 и В2, полученные при различных значениях m, kz и D для несжимаемой и сжимаемой жидкости и помещены в таблицу.

При учете сжимаемости для собственных значений с положительной действительной частью, частоты волн уменьшаются, что соответствует эффекту локализации возмущений вблизи внешней границы области. Для частот с отрицательной действительной частью, учет сжимаемости не приводит к заметному изменению частот и росту возмущений. Заметное совпадение в собственных значениях проявляется для m=1 и kz=; m=2 и kz=2 и m=3 и kz=3.

kz= -7.6+1.2 i -7.4+1.1 i -12.6+4.6 i -12.6+4.7 i -11.6+2.3 i -11.6+2.4 i kz=2 -9.0+6.2 i -8.7+6.0 i -9.8+6.9 i -9.6+6.6 i -11.7+0.2 i -11.9+0.6 i kz=2 -7.7+1.2 i -6.9+1.3 i -8.2+1.5 i -8.1+1.3 i -11.7+0.5 i -11.9+0.6 i -10.5+8.1 i -10.3+7.6 i -11.0+8.5 i -10.8+8.0 i -12.2+0.6 i -12.8+0.3 i На рисунках 1 и 2 представлены распределения yz ( x) для азимутального волнового числа m=1 и k z = p распределений магнитного поля B1 и B соответственно.

Рис.1. Решения Re yz ( x) (сплошная) и Im y z ( x) (пунктир) для l1 = 7.7 + 2.8i (кривые 1, 1`), l2 = 6.2 + 3.0i (кривые 2, 2`), l3 = -12.6 + 4.6i (кривые 3, 3`) и l4 = -12.6 + 4.7i (кривые 4, 4`).

Собственные значения l1 и l3 для r1 = 1, l2 и l4 для r1 = (1.1 - 0.2 x). r1 - безразмерное распределение плотности по слою.

Рис.2. Решения Re yz ( x) (сплошная) и Im y z ( x) (пунктир) для l1 = 7.7 + 2.8i (кривые 1, 1`), l2 = 6.2 + 3.0i (кривые 2, 2`), l3 = -11.6 + 2.3i (кривые 3, 3`) и l4 = -11.6 + 2.4i (кривые 4, 4`).

Собственные значения l1 и l3 для r1 = 1, l2 и l4 для r1 = (1.1 - 0.2 x). r1 - безразмерное распределение плотности по слою.

Распределения yz ( x) для собственных значений с отрицательной



Похожие работы:

«Отчет Главы Любинского муниципального района о результатах деятельности за 2011 год Приоритетной задачей Администрации Любинского муниципального района на 2011 год в соответствии с Планом действий в реализации социально-экономической политики являлось развитие устойчиво функционирующей экономики, обеспечивающей благосостояние граждан и сохранение социальной стабильности в Любинском районе. Кроме того, План действий на 2011 год был направлен на реализацию на территории Любинского муниципального...»

«Практика получения лицензий ФСТЭК и ФСБ в Московском регионе www.bis-expert.ru МОСКВА 2013 Ассоциация по вопросам защиты информации www.bis-expert.ru Информация об авторе: Андрей Валерьевич Семёнов, эксперт по информационной безопасности ООО Эдвансед Трансформейшн Консалтинг (AT Consulting) Опыт работы в области ИБ (более 10 лет): 2011 г. – по н/в ООО Эдвансед Трансформейшн Консалтинг Эксперт по ИБ 2004 г. – 2011 г. Банк Возрождение (ОАО) Главный специалист Службы информационной безопасности...»

«1 от 14 мая 2014 года № 265-р г. Горно-Алтайск О проекте закона Республики Алтай Об исполнении республиканского бюджета Республики Алтай за 2013 год На основании пункта 3 статьи 121 Конституции Республики Алтай: Одобрить разработанный проект закона Республики Алтай Об исполнении республиканского бюджета Республики Алтай за 2013 год для внесения его в установленном порядке на рассмотрение Государственным Собранием – Эл Курултай Республики Алтай. Исполняющий обязанности Председателя Правительства...»

«Инвестиционное предложение по развитию комплекса Аэроград Ноябрь 2010 Содержание 2 Предложение 3 План комплекса 4 План расположения объектов 5 Финансовые параметры проекта 6 Сводная информация по проектам 7 Создание производственной базы 8 Продвижение бизнеса и культуры 10 Оптовая и розничная торговля 16 Инфраструктура 18 Описание земельного участка 24 Украина: краткий обзор 25 Украина: информация о стране Украина: государственный строй Макроэкономический обзор Торговое и экономическое...»

«Стратегия развития морской портовой инфраструктуры России до 2030 года Москва 2012   1 Паспорт Стратегии развития морской портовой инфраструктуры России до 2030 года Наименование Стратегия развития морской портовой инфраструктуры Стратегии России до 2030 года Приказ Министерства транспорта Российской Основание Федерации от 30 июля 2010 года № 167 О создании для разработки рабочей группы по разработке Стратегии развития морских портов Российской Федерации Головной разработчик -...»

«УДК 55: 502:64 Барабошкина Т. А. Природно-ресурсный потенциал Юго-Западного Крыма Московский государственный университет имени М. В. Ломоносова, м. Москва e-mail: [email protected] Аннотация. На рубеже тысячелетий в странах СНГ одним из следствий спада интенсивности производства явилось сокращение техногенного прессинга на компоненты экосистем. Данный социальноэкономический феномен позволил решить задачу идентифицикации ведущих геологических факторов риска, оказывающих значимое влияние на...»

«6/2008 Официальное издание Федеральной таможенной службы Таможенные ведомости бюллетень таможенной информации В НОМЕРЕ: Регламент организации Инструкция о совершении законопроектной работы отдельных таможенных операций Федеральной таможенной службы при использовании таможенного режима переработки вне таможенной территории Обзор практики рассмотрения жалоб на решения, действия или бездействие О местах доставки товаров, перемещаемых таможенных органов в сфере таможенного дела железнодорожным...»

«ОТЧЕТ ГОДОВОМУ СОБРАНИЮ АКЦИОНЕРОВ ОАО Аэропорт Мурманск за 2006 год Общие сведения об Открытом акционерном обществе Аэропорт Мурманск Открытое акционерное общество “Аэропорт Мурманск. Свидетельство о государственной регистрации № 000426 от 25.01.1996г. Юридический адрес: 184355, Мурманская область, п. Мурмаши, Аэропорт Почтовый адрес: 184355, Мурманская область, п. Мурмаши, Аэропорт Контактный телефон: (815-2) 28-12-54, 28-12-59 Факс: (815-2) 28-15-41 Адрес электронной почты: airport@onlinе.ru...»

«  В.В. Максимов Генеральный директор ООО ТРАНСПРОЕКТ www.transproekt.ru Государственно-частное партнерство в транспортном секторе. Зарубежный опыт. В последние годы вопросам модернизации транспортной инфраструктуры уделяется особое внимание. Современные транспортные узлы, артерии – это одно из определяющих условий роста экономики. Именно транспортная инфраструктура определяет интенсивность хозяйственных связей, мобильность рабочей силы, товаров и услуг. Негативный опыт многих стран показывает,...»

«ПРЕДВАРИТЕЛЬНО УТВЕРЖДЕН УТВЕРЖДЕН решением Единственного акционера решением АО МНК КазМунайТениз Совета директоров (протокол заседания АО МНК КазМунайТениз Правления (протокол АО НК КазМунайГаз от 07.04.2011г. №04/11-о) от 28.04.2011 года №68) Годовой отчет АО МНК КазМунайТениз за 2010 год г. Астана, 2011 г. СОДЕРЖАНИЕ ОБРАЩЕНИЕ ГЕНЕРАЛЬНОГО ДИРЕКТОРА 1. АО МНК КАЗМУНАЙТЕНИЗ МАРАБАЕВА Е.Н.3 НОВЫЕ АКТИВЫ..5 2. ОПЕРАЦИОННАЯ ДЕЯТЕЛЬНОСТЬ.. 3. СОЦИАЛЬНЫЕ ОБЯЗАТЕЛЬСТВА, ПОДГОТОВКА КАЗАХСТАНСКИХ 4....»

«ЧЕТЫРЕ КНИГИ ОБ А Р Х И Т Е К Т У Р Е АНДРЕА ПАЛЛАДИО В КОИХ,ПОСЛЕ КРАТКОГО ТРАКТАТА О П Я Т И ОРДЕРАХ И НАСТАВЛЕНИЙ НАИБОЛЕЕ НЕОБХОДИМЫХ ДЛЯ СТРОИТЕЛЬСТВА, ТРАКТУЕТСЯ О МАСТНЫХ ДОМАХ, Д О Р ОГАХ, МОСТАХ, ПЛОЩАДЯХ, КСИСТАХ И ХРАМАХ. В ПЕРЕВОДЕ АКАДЕМИКА АРХИТЕКТУРЫ И.В.ЖОАТОВСКОГО * * # ИЗДАТЕЛЬСТВО ВСЕСОЮЗНОЙ АКАДЕМИИ АРХИТЕКТУРЫ...»

«ЗАКОНОДАТЕЛЬНОЕ СОБРАНИЕ ПЕРМСКОГО КРАЯ ПРЕСС-РЕЛИЗЫ Проект закона Пермского края О порядке ведения органами местного самоуправления учета граждан в качестве нуждающихся в жилых помещениях, предоставляемых по договорам социального найма (первое чтение, инициатива Главы города Перми) Законопроектом предлагается внести изменения в Закон Пермского края от 01.12.2011 № 2694-601 О порядке ведения органами местного самоуправления учета граждан в качестве нуждающихся в жилых помещениях,...»

«Памятные места историко-культурного наследия Ивановской области Правительство Ивановской области Департамент образования Ивановской области Департамент культуры и культурного наследия Ивановской области Департамент спорта и туризма Ивановской области ГТРК Ивтелерадио Благотворительный фонд Гордость Отечества Памятные места историко-культурного наследия Ивановской области Дорогие друзья! Проект Путешествие по губернии посвящен 95-летию Иваново-Вознесенской губернии и призван стать одним из...»

«Приложение № 3 к Протоколу заседания Совета директоров ОАО Полюс Золото № 04-09/СД от 13 апреля 2009 года Предварительно утвержден Советом директоров ОАО Полюс Золото Протокол заседания Совета директоров № 04-09/СД от 13 апреля 2009 года Утвержден годовым Общим собранием акционеров ОАО Полюс Золото Протокол Общего собрания № 01-09/ОСА от 26 мая 2009 года ГОДОВОЙ ОТЧЕТ ОАО ПОЛЮС ЗОЛОТО ЗА 2008 год Генеральный директор ОАО Полюс Золото Е.И. Иванов 07 апреля 2009 года Главный бухгалтер ОАО Полюс...»

«6.1.Состояние материально-технической базы АНО ВПО КИУ.191 6.2.Состояние и развитие учебно-лабораторной базы, уровень ее оснащения..192 6.3. Социально-бытовые условия АНО ВПО КИУ.195 ЗАКЛЮЧЕНИЕ..194 Приложение 1. ПОКАЗАТЕЛИ ДЕЯТЕЛЬНОСТИ ОБРАЗОВАТЕЛЬНОЙ ОРГАНИЗАЦИИ ВЫСШНГО ОБРАЗОВАНИЯ, ПОДЛЕЖАЩЕЙ САМООБСЛЕДОВАНИЮ Автономная некоммерческая организация высшего профессионального образования Калининградский институт управления 1. Общие сведения об образовательной организации Институт – Автономная...»

«ВВЕДЕНИЕ Гидротехнические мелиорации представляют собой комплекс мероприятий, направленных на регулирование водного режима почв и улучшения режима питания путем осушения избыточно увлажненных земель и орошения земель с недостаточным увлажнением. Чаще всего конечными целями гидротехнических мелиорации в лесном хозяйстве являются увеличение продуктивности лесных земель, повышение санитарно-гигиенических и рекреационных функций лесов и устойчивости их к неблагоприятным воздействиям. В сочетании с...»

«МЕЖОТРАСЛЕВОЙ КАТАЛОГ ПРОМЫШЛЕННАЯ ЭКОЛОГИЯ-2012 Межотраслевой каталог для специалистов промышленных предприятий поможет найти актуальную информацию о ведущих производителях экологического оборудования, ООО ИНТЕХЭКО проектных и инжиниринговых компаниях, имеющих www.intecheco.ru реальный опыт решения различных экологических задач в промышленности, газоочистки, пылеулавливания, водоподготовки и водоочистки, переработки отходов, экологического мониторинга и охраны окружающей среды. МЕЖОТРАСЛЕВОЙ...»

«Структура Проекта перспективного развития МБОУ СОШ №7 городского округа г. Стерлитамак РБ на 2011-2015 гг. Паспорт Проекта перспективного развития... -3Введение.... -4Постановка проблемы 2.1.Анализ развития МБОУ СОШ №7 городского округа г. Стерлитамак РБ 2.1.1.Текущее ресурсное обеспечение... -5Режим работы... -7Внешние связи... -8Анализ состояния образовательного процесса... -9Кадровый состав... -9Контингент обучающихся... -16Воспитательная работа... -19Проблемное поле для...»

«ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ Байкальский государственный университет экономики и права Институт национального развития при администрации президента Монголии и МАН Иркутское отделение Вольного экономического общества Модернизация социально-экономического развития региона Сборник научных трудов Иркутск Издательство БГУЭП 2009 PDF created with pdfFactory Pro trial version www.pdffactory.com УДК 332.1(57) ББК 65.04 М Печатается по решению редакционно-издательского совета Байкальского...»

«АНАЛИЗ ПОДЗАКОННЫХ АКТОВ РЕСПУБЛИКИ ТАДЖИКИСТАН В ОБЛАСТИ ЛИЦЕНЗИРОВАНИЯ 20 июня 2008 г. Данный анализ опубликован благодаря помощи американского народа, предоставленной Агентством США по международному развитию (USAID). Анализ был подготовлен Нигиной Салибаевой, кандидатом юридических наук, доцентом кафедры международного права ТГНУ и Проектом USAID по улучшению бизнес среды. АНАЛИЗ ПОДЗАКОННЫХ АКТОВ РЕСПУБЛИКИ ТАДЖИКИСТАН В ОБЛАСТИ ЛИЦЕНЗИРОВАНИЯ ОГОВОРКА Мнение автора, высказанное в данной...»






 
2014 www.av.disus.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Программы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.