Электромагнитные зондирования плато Гиза (Египет)

Построение 3D модели электрических сопротивлений горных пород до глубины 150 метров на основе результатов импульсных электромагнитных зондирований территории Сфинкс. Описание обработки данных и дизайна модели среды как результата 3D блоковой инверсии.

Рубрика Геология, гидрология и геодезия
Вид статья
Язык русский
Дата добавления 02.02.2019
Размер файла 366,5 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

УДК 550.837

ЭЛЕКТРОМАГНИТНЫЕ ЗОНДИРОВАНИЯ ПЛАТО ГИЗА (ЕГИПЕТ)

П.О. Барсуков

Э.Б. Файнберг

Об авторах:

1 - Центр Геоэлектромагнитных Исследований Института Физики Земли РАН, Москва

В 1994 г. в рамках совместного проекта НРИАГ (Египет) и ИФЗ РАН (Россия) [1] были проведены исследования территории Сфинкс - п. Хефрена импульсным электромагнитным методом. Расположение пирамид, Сфинкса и рельеф местности показаны на рис. 1. Целью исследований была оценка конфигурации водоносных горизонтов и резервуаров грунтовых вод и комбинирование полученных данных с результатами ранее проведённых НРИАГ вертикальных электрических зондирований (ВЭЗ) в окрестностях Сфинкса. Для возбуждения и измерения полей использовались квадратные совмещённые антенны размером 50Ч50 м и технология TEM-FAST [2] зондирований. Отклики среды на импульсное возбуждение поля измерялись в диапазоне 4-1000 µs, что обеспечивало исследование электрического сопротивления пород в интервале глубин 5-150 м.

Рис. 1. Территория исследований.

А - изолинии - рельеф земной поверхности (в метрах), точечные линии - подземные входы в пирамиды. Б - пункты зондирований показанные на рис. 1В и профиль зондирования (ломаная линия) вдоль которого строился разрез. B - переходные отклики в виде кажущегося сопротивления

1. Технология TEM-FAST зондирований

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

Отметим особенности технологии TEM-FAST зондирований. Обычно, после того, как импульсный отклик измерен, его трансформируют в зависимость кажущегося сопротивления от времени a(t), используя известную асимптотическую формулу становления поля для поздних времён. Примеры таких откликов приведены на рис. 1. Однако для построения 3D моделей среды такой нормировки недостаточно. В технологии TEM-FAST применяется трансформация откликов в зависимости кажущегося сопротивления от глубины исследования (h). Алгоритм трансформации содержит три последовательных этапа: 1) аппроксимация измеренных откликов функцией E(t), представляющей собой сумму затухающих во времени экспоненциальных членов [8]: 2) вычисление сопротивления однородного полупространства (t), отклик от которого совпадает с E(t) на всех временах t; 3) вычисление параметрически заданной функции (h) для каждого параметра t.

Пример трансформаций (t) и (h) иллюстрируют рис. 2А и Б, соответственно. горный порода электромагнитный зондирование

Трансформации сглаживают кусочно-однородное распределение сопротивлений среды, несколько искажая значения и глубины h. Тем не менее, основные особенности модели сохраняются. Кроме того, важно, что при таком подходе для дизайна устойчивой трансформации не требуется никакой априорной информации об изучаемой среде; локальные 1D трансформации (h), «сшитые» в 2D или 3D структуры, создают электромагнитный образ среды, который в 2D варианте называется псевдо-разрезом сопротивлений.

Рис. 2. А - импульсные отклики в виде кажущегося сопротивления; Б - ?а(t) и ?(t) - зависимость кажущегося сопротивления от глубины исследования. Кусочно-однородная кривая - результат 1D инверсии отклика, полученного при зондировании площадки под Сфинксом

Следующим шагом построения 3D модели среды является проецирование трансформаций на вертикальную сетку модели, т.е. замена непрерывной функции (h) (рис. 3а) её ступенчатой, кусочно-однородной аппроксимацией. Количество слоёв с разными сопротивлениями может достигать нескольких десятков в зависимости от необходимой детальности исследований. Сетка анализируемой здесь 3D модели состояла из одинаковых ячеек размером x =y=25 м, z= 2 м. На рис. 3 А показано сечение 3D модели, «сшитой» из локальных 1D моделей. Положение профиля показано на рис. 1Б, а экспериментальные данные и рассчитанные вдоль профиля 3D отклики показаны на рис. 3Б,В. Видно, что невязка (rms) между наблюдёнными и модельными данными недопустима высока.

На следующем шаге производится 1D инверсия откликов в рамках многослойных локальных моделей в каждой точке измерения. В качестве стартовых моделей используются ранее построенные на основе трансформаций (h) ступенчатые структуры. В процессе 1D инверсии подбираются сопротивления слоистых структур, при которых минимизируется невязка между наблюдёнными и расчётными откликами в каждой точке измерения. Дополнительным ограничением на значения сопротивлений слоёв служит дискретная шкала сопротивлений, содержащая ограниченное количество градаций (в данном случае 10). Результат 1D инверсии вдоль профиля показан на рис. 3Б. На верхней панели показаны экспериментальные данные и рассчитанные 3D отклики. Видно, что невязка уменьшилась, однако по-прежнему оставляет желать лучшего.

На заключительном этапе производится 3D инверсия данных в рамках построенной блоковой структуры, т.е. минимизируется невязка между наблюдёнными и модельными откликами, при этом подбираются сопротивлениями блоков без изменения их границ. В данном случае инверсия осуществлялась относительно 10 сопротивлений в рамках шкалы 1-1000 Ом-м. Результат показан на рис. 3В. Невязка между наблюдёнными и модельными данными оказалась в пределах 12% для t=250 мкс. Общая невязка по всему временному диапазону rms ~8%., что в целом соответствует погрешности данной съёмки.

Рис. 3. Разрезы сопротивлений 3D модели по профилю, показанному на рис. 1Б.

A - структура, созданная на основе трансформаций (h);

Б - структура, созданная на основе локальных 1D инверсий;

В - результат 3D инверсии

Анализируя результаты инверсии на разных этапах можно заметить, что главные особенности модели проявляются уже на начальном этапе трансформаций (рис. 3A) - это три блока с высокими сопротивлениями, разбитые субвертикальными низкоомными швами. Разрезы, представленные на панелях Б и В, принципиально не отличаются, однако их 3D отклики более контрастны и лучше согласуются с наблюдёнными данными.

2. Анализ геоэлектрической модели

Созданная геоэлектрическая 3D модель показывает следующее. Сопротивление поверхностного слоя (первые 5 м) меняется в пределах от 10 до 500 Ом-м, достигая минимума восточнее Сфинкса. Этот результат полностью согласуется с данными ВЭЗ, полученными в окрестности (100 м) монумента. Здесь также при разносе AB/2 < 5 м сопротивление не превышает 20 Ом-м. Глубже идёт возрастание сопротивления пород до глубины ~10 м от поверхности, что также подтверждается ВЭЗ [5] . На горизонте 0 м (отсчитывается от уровня моря) центр и южное крыло территории имеет сопротивление более 200 Ом-м, в то время как северное крыло имеет сопротивление ~10-20 Ом-м. В центральной части прослеживается низкоомный шов мощностью 30-50 м SW-NE направления. На горизонте -50 м общее сопротивление пород падает до 5-10 Ом-м, за исключением двух областей. На глубоком горизонте -100 м виден единственный блок у южной грани п. Хефрен сопротивлением более 200 Ом-м, на фоне практически однородной среды с сопротивлением менее 10 Ом-м. На вертикальном разрезе модели (рис. 3В) видны три блока пород с сопротивлением более 200 Ом-м, разделённые субвертикальными швами сопротивлением менее 10 Ом-м. Все блоки подстилаются толщей пород сопротивлением менее 10 Ом-м. Восточный, самый мощный блок, на котором расположен Сфинкс, прослеживается до горизонта - 50 м.

3. Геологическая интерпретация.

Как известно, плато пирамид Гиза сложено Средне-Эоценовыми известняками формации Моккатам; породы этого горизонта в регионе нигде, кроме плато, не встречаются [7]. Известняки моноклиинально падают в SE направлении под углом 5-7o. Толща Моккатам несогласно залегает на позднемеловых карбонатах, которые интенсивно смяты в результате позднемеловой тектонической активности. В работе [6] приведён фрагмент геологической карты, на которой видна сложная складчатая структура Абу-Раваш (9 км NW от плато Пирамид), сложенная Меловыми породами, в пределах которой формация Моккатам отсутствует. Меловые отложения перекрыты более молодыми породами Верхне-Эоценовой Маади формацией. Мощность формации Моккатам не определена. В NW направлении подошва Моккатам обнажается на горизонте +100 м в 2.5 км от плато Пирамид. Если предположить, что подошва толщи плоская, можно оценить максимальную мощность Моккатам в основании Сфинкса (+25 м) - Hmax~190 - 280 м при углах падения 5-70 SE, соответственно. Сопротивление пород в созданной 3D модели изменяется почти на 3 порядка, что явно противоречит предположению о том, что мы имеем дело только с латерально однородной толщей известняков. Столь резкая вертикальная неоднородность сопротивлений в модели объясняется сложным рельефом кровли подстилающей толщи Меловых отложений. Так, например, в структуре Абу-Раваш пачки пород изменяют угол и вектор падения с 150 SE на 150 NW на расстоянии ~2 км, образуя “мини депрессию” глубиной ~250 м. Эта же складчатая структура подстилает и плато Пирамид [3,4].

Таким образом, высокоомные блоки модели сложены плотными известняками формации Моккатам и разбиты обводнённой разломной зоной SW-NE направления. Максимальная мощность плотных известняков под южной гранью п. Хефрен (+65 м) - не менее 160 м, а вблизи Сфинкса (+25 м) ~80 м. Низкоомные подстилающие горизонты представляют собой смятые в складки трещиноватые меловые карбонаты, насыщенные минерализованными водами. В Западной части плато (0.5-1 км западнее п. Хефрен), где были сделаны «калибровочные» зондирования, импульсные отклики аналогичны sp24 (рис. 2), т.е. сопротивление поверхностного слоя (25 м) > 500 Ом-м. Далее наблюдается 100-кратное падение сопротивления до 5 Ом-м на глубине 45 м и до 1 Ом-м на глубине 70 м. Верхний слой соответствует плотным, слабопроницаемым известнякам, а подстилающая толща представлена насыщенными минерализованной водой трещиноватыми породами. Таким образом, граница высокоомных и низкоомных образований в созданной модели соответствует границе меловых и палеогеновых отложений. Схематический геологический разрез плато Пирамид по линии Сфинкс-Хефрен показан на рис. 5.

Рис. 4. Схематический геологический разрез плато Пирамид. Пунктир - разломная зона

Отметим, что соответствие данных «калибровочных» измерений и откликов, полученных на основной территории исследований, свидетельствует об отсутствии существенных погрешностей ЭМ съёмки, обусловленных обилием рукотворных объектов на плато пирамид. Поверхностный (~5м) слой Моккатам известняков характеризуется весьма неоднородным электрическим сопротивлением., что подтверждается ВЭЗами [6] В работах [3,4] отмечается, что нижняя пачка формации Моккатам (основание Сфинкса) представлена глинистыми известняками с вкраплением галита, а верхняя пачка (п. Хефрен) - массивными рифовыми известняки (биогермы). По-видимому, этим можно объяснить вариации сопротивления в приповерхностном слое плато пирамид Гиза.

Библиографический список

1. Barsukov P.O, Hassaneen A.Gh., Svetov B.S., Osman S.Sh. Geoelectrical Study in the Vicinity of Giza Pyramids, Egypt //13-th Workshop on Electromagnetic Induction in the Earth, 1996. Japan

2. Barsukov P.O., Fainberg E.B., Khabensky E.O. Shallow investigation by TEM-FAST technique: methodology and case histories // Methods of geochemistry and geophysics, 40, Elsevier (Ed. Spichak V.V.), 2007. C. 55-77

3. Gauri K.L. Geologic Study of the Sphinx. ARCE Newsletter, No. 127, 1984. C. 24-43

4. Guirand R., Bosworth W. Phanerozoic geodynamic evolution of northeastern Africa and the northwest Arabian Platform //Tectonophysics 315 (1-4),1999. C. 73-108

5. NRIAG (National Research Institute of Astronomy and Geophysics). Report of the first six months of underground water monitoring at Sphinx area to Egyptian Authority of Monuments, 1993, Egypt

6. Raynaud S., Boisse H., Makroum F.M., Bertho J. Geological and Geomorphological study of the original hill at the base of Fourth Dynasty Egyptian monuments // Bulletin de la Societe Geologique de France, V. 181, no. 3, 2010. C. 279-290

7. Said R. The Geology of Egypt. Elsevier Publishing Co., Amsterdam-New York, 1962. 377 C.

8. Svetov B. S., Barsukov P. O. Transformation of the quasi-stationary transient processes in the geoelectrically equivalent waves // Physics of the Earth, V.8, 1984. C. 29-37

Аннотация

На основе результатов импульсных электромагнитных зондирований территории Сфинкс - Пирамида Хефрена (Гиза, Египет) построена 3D модель электрических сопротивлений горных пород до глубины 150 м. Описана технология зондирований, обработки данных и дизайна модели среды как результата 3D блоковой инверсии. На основе геологической интерпретации модели впервые определена граница палеогеновых известняков, слагающих плато пирамид, и подстилающей складчатой меловой структуры.

Ключевые слова: электромагнитные зондирования, TEM-FAST технология, 3D инверсия, геоэлектрический разрез, плато Гиза, Пирамиды, Сфинкс.

Based on the results of pulsed electromagnetic sounding of the territory of the Sphinx - the Pyramid (Giza, Egypt) built a 3D model of electrical resistance of rocks to a depth of 150 m. the technique of intubation, data and design models the environment as a result of 3D block inversion. Based on the geological interpretation of the model is identification of the boundary of the Paleogene limestones composing the plateau of the pyramids, and underlain by folded Cretaceous structure.

Keywords: electromagnetic sounding, TEM-FAST technology, 3D inversion, the geoelectric, Giza plateau, great Pyramids, Sphinx.

Размещено на Allbest.ru

...

Подобные документы

  • Электромагнитные свойства минералов и горных пород. Электрохимическая активность и поляризуемость. Пьезоэлектрические модули у турмалина. Особенности применения электроразведки. Майко-Лебедское рудное поле. Методы полевых электромагнитных зондирований.

    презентация [1,6 M], добавлен 30.10.2013

  • Исторический образ, обзор первобытной обработки камня. Залегания горных пород и их внешний вид. Структура, текстура горных пород Южного Урала. Способы и оборудование для механической обработки природного камня. Физико-механические свойства горных пород.

    курсовая работа [66,9 K], добавлен 26.03.2011

  • Подготовка горных пород к выемке. Параметры взрывных работ. Определение парка буровых станков карьера. Выбор модели экскаватора-мехлопаты (для экскавации полезного ископаемого). Транспортировка горной массы. Выбор модели бульдозера, фронта разгрузки.

    курсовая работа [486,7 K], добавлен 21.12.2011

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

    курсовая работа [1,4 M], добавлен 24.06.2011

  • Проверка крайних значений вариационных рядов по проходке интервала от 400 до 2100 метров. Проверка однородности пачки одинаковой буримости. Выбор типа буровых долот по механическим свойствам горных пород и порядок определения осевой нагрузки на долото.

    курсовая работа [61,2 K], добавлен 24.10.2012

  • Общее описание и характерные черты осадочных горных пород, их основные свойства и разновидности. Типы слоистости осадочных горных пород и структура. Содержание и элементы обломочных пород. Характеристика и пути образования химических, органогенных пород.

    реферат [267,1 K], добавлен 21.10.2009

  • Классификация горных пород по происхождению. Особенности строения и образования магматических, метаморфических и осадочных горных пород. Процесс диагенеза. Осадочная оболочка Земли. Известняки, доломиты и мергели. Текстура обломочных пород. Глины-пелиты.

    презентация [949,2 K], добавлен 13.11.2011

  • Сущность интрузивного магматизма. Формы залегания магматических и близких к ним метасоматических пород. Классификация хемогенных осадочных пород. Понятие о текстуре горных пород, примеры текстур метаморфических пород. Геологическая деятельность рек.

    реферат [210,6 K], добавлен 09.04.2012

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

    курсовая работа [417,1 K], добавлен 11.12.2012

  • Методы определения возраста горных пород, слагающих Землю. Возраст пород слоя Базальт Карденас в восточной части Большого Каньона. Геологическая “блоковая" схема расположения пластов горных пород Большого Каньона. Ошибки радиологического датирования.

    реферат [1,4 M], добавлен 03.06.2010

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

    курсовая работа [878,7 K], добавлен 18.04.2015

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

    контрольная работа [555,6 K], добавлен 29.11.2012

  • Исследование особенностей осадочных и метафорических горных пород. Характеристика роли газов в образовании магмы. Изучение химического и минералогического состава магматических горных пород. Описания основных видов и текстур магматических горных пород.

    лекция [15,3 K], добавлен 13.10.2013

  • Декриптометрические методы исследования минералов, пород и руд, их распространение. Типизация вакуумных декриптограмм пород гранитоидного ряда. Обработка и интерпретация результатов вакуумно-декриптометрических анализов метасоматически измененных пород.

    контрольная работа [702,3 K], добавлен 21.06.2016

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

    реферат [25,0 K], добавлен 19.03.2009

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

    курсовая работа [2,6 M], добавлен 08.02.2016

  • Характеристика структуры, изучение строения и определение размеров пор горных пород. Исследование зависимости проницаемости и пористости горных пород. Расчет факторов проницаемости и методов определения содержания в пористой среде пор различного размера.

    курсовая работа [730,4 K], добавлен 11.08.2012

  • Проведение на электронных вычислительных машинах имитационных лабораторных испытаний горных пород и определение их механических свойств (пределов прочности, модуля упругости и коэффициента Пуассона). Теории определения прочности горных пород Кулона-Мора.

    курсовая работа [3,8 M], добавлен 27.06.2014

  • Основные стадии процесса добычи полезного ископаемого. Предел прочности горных пород при растяжении, методы и схемы определения, количественная оценка. Деформация твердого тела. Методы определения хрупкости горных пород. Хрупкое разрушение материала.

    реферат [303,3 K], добавлен 14.02.2014

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

    контрольная работа [956,1 K], добавлен 10.01.2014

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.