Построение моделей сложнопостроенных геологических сред по гравиметрическим данным

Тип работы:
Реферат
Предмет:
Общие и комплексные проблемы естественных и точных наук


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 550. 3: 004. 02
ПОСТРОЕНИЕ МОДЕЛЕЙ СЛОЖНОПОСТРОЕННЫХ ГЕОЛОГИЧЕСКИХ СРЕД ПО ГРАВИМЕТРИЧЕСКИМ ДАННЫМ
1Кобрунов А.И., 1Мотрюк Е.Н., 2Ломинский Д.О.
ФГБОУ ВПО «Ухтинский государственный технический университет»,
Ухта, e-mail: akobrunov@ugtu. net- 2ФГБОУ ВПО «Российский государственный университет нефти и газа имени И.М. Губкина», Москва
Для условий слабой изученности и сложного строения осадочных бассейнов выработаны модельные представления геологических сред по данным гравиразведки. При решении задач реконструкции геологической среды построенная модель распределения физического параметра должна соответствовать с заданной степенью приближения наблюдаемым физическим полям. При этом возникают три взаимоувязанных элемента, составляющих содержание понятия модели задачи: модель среды, модель поля и модель физического явления. Рассмотрены постановки прямой и обратной задач гравиразведки. Прямая задача состоит в нахождении модели поля по известным моделям среды и физического явления. Обратная задача заключается в реконструкции модели среды по модели поля и явления. Представлены различные способы представления информационной модели: одномерные объекты, двухмерные объекты и трёхмерные наблюдения, для которых применяются различные способы интерпретации данных.
Ключевые слова: гравиразведка, сложнопостроенные среды, модели среды, модели поля, прямые задачи гравиразведки, обратные задачи гравиразведки
CONSTRUCTION OF MODELS OF COMPLEX GEOLOGICAL ENVIRONMENTS OF GRAVIMETRIC DATA
1Kobrunov A.I., 1Motryuk E.N., 2Lominskiy D.O.
1FGBOU VPO «Ukhta State Technical University», Ukhta, e-mail: akobrunov@ugtu. net-
2FGBOU VPO «Russian State University of Oil and Gas named after Gubkin», Moscow
Model representations of geological environments according to the gravity prospecting for conditions and little knowledge of the complex structure of sedimentary basins are developed. In solving the problems of reconstruction of the geological environment constructed physical parameter distribution model must comply with the specified degree of approximation the observed physical fields. Thus there are three interrelated elements that make up the content of the concept model of the problem: a model of the environment, the model fields and model of a physical phenomenon. Consider the formulation of direct and inverse gravity. The direct problem consists in finding the model field from the known models of the environment and physical phenomena. The inverse problem is the reconstruction model on model fields and event. We present different ways of presenting the information model: one-dimensional objects, two-dimensional objects, and three-dimensional observation, for which there are different ways to interpret the data.
Keywords: gravity prospecting, complex structure environment, model of environment, model of the field, direct gravity problem, inverse problem of gravity
Для построения геолого-геофизических моделей осадочных бассейнов и их фрагментов эффективно применение ведущих современных технологий интерпретации гравиметрических данных. Для крупных плотностных структур типа осадочных бассейнов необходимо проводить постановку и решение интерпретационных задач с учётом условий слабой изученности и сложного строения среды. Это позволяет сделать технология, основанная на критериальном подходе к решению обратных задач гравиразведки [1−5, 7]. Она основана на комплексном анализе всей имеющейся геолого-геофизической информации. Возникает проблема выбора модели среды, которая должна быть положена в основу изучения имеющихся в осадочных бассейнах образований.
При решении задач реконструкции геологической среды построенная модель распределения физического параметра должна соответствовать с заданной степенью приближения наблюдаемым физическим полям. Три взаимоувязанных элемента, составляющих содержание понятия модели задачи: модель среды, модель поля и модель физического явления. Выделяются прямая и обратная задачи геофизики. Во введённых терминах прямая задача состоит в нахождении модели поля по известным моделям среды и физического явления. Обратная задача заключается в реконструкции модели среды по модели поля и явления. Исходную информацию, имеющуюся для изучения геологической среды, можно представить в виде совокупности трёх компонент [2, 5, 6, 8]:
одномерные объекты, например данные скважин- двухмерные объекты, т. е. некие профильные наблюдения, и трёхмерные наблюдения, в частности наблюдаемые геофизические поля.
Наиболее общей и полной в смысле количества априорной информации моделью среды является модель распределения плотности, или плотностная модель. Она может быть использована для решения задач локального прогнозирования, т. е. оценки локальных изменений физических свойств, в пределах отдельных блоков, пластов и разреза в целом. Эта модель является математической основой схем интерпретации на классе непрерывных функций или прогнозирования непрерывного геолого-геофизического разреза. Структурная геолого-геофизическая модель является упрощённым вариантом. Введение данной модели необходимо для решения задач, где априори вводится предположение о слоистом строении среды. В рамки таких предположений укладываются как классические структурные задачи — изучение поверхности кристаллического фундамента, структур синклинального и антисинклинального типов, надвиговых структур, так и задачи соляно-купольной тектоники.
Для описания моделей среды используются различные варианты задания данных и соответствующие им постановки прямых задач гравиразведки для двухмерного и трёхмерного случая. Суть двухмерной постановки состоит в следующем [2, 5, 6]. Часто при интерпретации гравиразведочных данных аномалии гравитационного поля носят ярко выраженный линейный характер. Тогда естественно ввести предположение о том, что источник поля имеет бесконечную протяжённость по какому-то направлению, а в каждой вертикальной плоскости, перпендикулярной к этому направлению, сечение носителя и распределение масс по этому сечению одни и те же. Такие тела принято называть двухмерными, а соответствующие им поля -двухмерными или плоскими.
Плотностная двухмерная модель (модель распределения плотности). Задание исходной информации производится по профилю. Используем прямоугольную систему координат с осью Ог, направленной вниз к массам, и осью Ох, совмещённой с дневной поверхностью и направленной вдоль линии задания данных. Модель среды конструируется в области S, лежащей в нижнем полупространстве. Для каждой точки (х, г) ставим в соответствие значения параметра плотности о = о (х, г).
Структурная двухмерная модель. Для представления структурной модели каж-
дой точке х сопоставляется глубина залегания г = ?(х) соответствующей границы к = 0, 1, …, N лежащей в области S, ограниченной горизонтальной полосой П. Каждый из пластов характеризуется своим неизменным по вертикали параметром плотности ок (х).
Аппроксимационные двухмерные модели, являясь тоже разновидностью плот-ностных, представляют собой разбиение области S сеткой на ячейки S. с плотностью о. так, что uSJ =?- пЯ- =0-. = 1, …, J конфигурация которых подбирается таким образом, чтобы изучаемый объект можно было ими хорошо приблизить.
Исходные данные в условиях недоопре-деленности задаются системой профилей Г, Г = {ГД, I = 1, …, Р, I — геолого-геофизические разрезы. Систему координат рассматриваем с осью Ог, направленной вниз и плоскостью хОу — совмещённой с дневной поверхностью. Модель среды конструируется в области V, лежащей в нижнем полупространстве. Для всех разрезов имеются координаты начала и конца в общей системе координат.
Плотностная объёмная модель. Для представления модели данного вида каждой точке V = (х, у, г) пространства V сопоставляется значение параметра плотности о = о^).
Структурная объёмная модель. Каждой точке 5 = {х, у} пространства V, ограниченного полосой П, сопоставляется глубина залегания г = ^(я) границы к = 0, 1, …, N. Каждый из пластов характеризуется своим параметром плотности ок (5), неизменным по вертикали в пределах каждого из пластов. Среди структурных моделей можно выделить:
1) модели, где плотность ок (5) = ок постоянна в пределах каждого из N + 1 пластов-
2) модели, где плотность ок (5) может меняться не только в пределах к-го пласта, но и в пределах q блоков, находящихся в составе данного пласта.
Модель первого вида является идеализированной, поэтому наибольший интерес представляет модель, где ок (5) Ф ок, т. е. плотность пласта не является постоянной. В этом случае модель представляет собой следующее:
• границы, ограничивающие пласты, представляют собой однозначные функции пространственных координат г = УХ5), к = 0, 1, …, N + 1-
• плотность пласта, заключённого между к-й и к + 1-й границами, есть функция горизонтальных координат ок+1 = ок+1(5), к = 0, 1, …, Л-
• нулевая граница есть горизонтальная пластина с глубиной/0 и плотностью пласта с0-
• плотность среды ниже границы с номером N есть сЛ+1-
• величины ?, с, г, с" считаются посто-
•& gt- 0'- N+1'- 0
янными- ^
• объёмная модель {/, Дб} имеет в качестве своего следа на поверхности {/} -проходящей через линию Г1 нормально к дневной поверхности двухмерную модель {МВД.
Для объёмной структурной модели введём краткое обозначение:
{Г, Ад}, / = {/0,/,… /"}, Д5= {Дсц, Дс^,… АстЛГ},
= ~ °к+" °0 = СТЛГ+1 = °& gt-
где Док — контрастность к-го пласта.
Аппроксимационные объемные модели представляют собой разбиение области V сеткой на подобласти V с плотностью о.
так, что иГ. = V- V. п V. = 0:. = 1, …, J.
Классические аппроксимационные модели характеризуются призматическим видом подобластей. Модели этого вида могут иметь в качестве аппроксимирующих тел также и другие тела: шар, пирамиду и др. Такие конструкции ориентированы на задачи типа рудных.
Иерархическая схема объёмных моделей среды, используемых гравиразведкой
и в дальнейшем обозначаемых x (v), приведена на рис. 1.
В зависимости от вида модели среды модель явления, которая суть оператор прямой задачи A, может описывать плотност-ные и структурные задачи гравиразведки.
Ax = U. (1)
Моделью поля в этих задачах можно считать значения вертикальной составляющей гравитационного потенциала в конечном числе точек. В случае двухмерных задач это будет двухмерный массив (x0, u), в случае трёхмерной — трёхмерный (x0, y, u). Под решением обратной задачи [1,2, 8] по -нимают нахождение распределения источников x (плотность o (v) либо конфигурация границ/s)) при введённых о них модельных представлениях по заданным наблюдаемому полю U и оператору прямой задачи A (x).
Оператор прямой двухмерной плот-ностной задачи гравиразведки. Вертикальная производная гравитационного потенциала u (x0, z0), где точка (x0) регистрируется на поверхности в E+(z & gt- 0) с уравнением z0 = y (x0), задана следующим соотношением между uz (x0, z0) и g = o (x, z):
s [(х-х0) + (z-|/(*0)) J
(2)
где у = 6,673−10 5 (поле в миллигалах, расстояния в метрах, плотность в г/см3) — гравитационная постоянная.
Рис. 1. Иерархическая схема объёмных моделей среды
Оператор прямой двухмерной струк- циала где v0 = (х0, у0, г0) и точка (х0,
турной задачи гравиразведки. Для струк- у0) = 50 е Ё регистрируется на поверхности
000
с уравнением г = у (50), задана следующим
турной модели среды, где система гра-
Л- / г /& quot-^л! соотношением между и (О и о = о^):
х)-11х)] имеет плотности ^ ^ 0 у '-
а (-с)={а1. (л-)}, связь с вертикальной производной гравитационного поля иг (х0, г0) определена соотношением
411,
Л (о (у))=и2 (я0,20=у (*0))= уа () (г —))(Ы, ус1г
пЗ/2
^(7(х), а (х))=мг (х0,г0 =|/(х0))= к -Хо)2 + -у0/ + (г-))2]
=11
• (4)
уАо,. (х) с1х
(3)
«[(*-*0)2 + (у-(х)-х)/(х0))2|2
Оператор прямой трёхмерной плот-ностной задачи гравиразведки. Вертикальная производная гравитационного потен- ется по формулег:
Оператор прямой трёхмерной структурной задачи гравиразведки. Для структурной модели среды для простоты плотность каждого пласта примем постоянной. Вертикальная производная гравитационного потенциала и0), в точке А0), вычисля-
•уа (/ (я) — г0) йхйф
Формула (5 а) для объёмной структурной модели среды в случае, когда г0 = у (50):
(5а)


[(х-Х0)2 +(у-у0)2 +(/(5)-|/(50))2]
(5Ь)
Вычислив в (5Ь) интеграл по переменной г и представив результат в виде суммы по N + 1 границе, получаем
С/г (50,70=У (50)) = ^({7,Ад})=ХЛ
уАа ?(я) с! х& lt-}у
/=0 5
[(х-х0)2 + Су-у0)2 +
,½ • (5с)
Оператор прямой задачи гравиразведки
Ах ~ и
Плотностная задача
Л (о (V)) = гф, = ИЛ))) =
Фк-ьг+ь-уУ+ь-уъо"2!& quot-
^ ч
/ Структурная задача
^угг_уйо^'-Шу
Прямая зад ача для Ч
аппроксныацнониой модели среды
ув/у& amp-г- ^?^УЬсОусЬ
-ГО
5 к*-^+0'--^+(*- к-иРГ
Рис. 2. Оператор объёмной прямой задачи гравиразведки для разных видов моделей среды
Оператор прямой трёхмерной задачи гравиразведки для аппроксимационной модели. Для каждого элемента V. задаётся значение плотности Oy., а гравитационное поле рассчитывается по формуле

j=1
где
ai- ttJi — fff yo (v)(z-y (s0))dxdydz
A^yut (vaJ=JJJr-----ту. (6)
Vj [(x-X0) +(y-y0) +(z-|/(50)) J
Виды модели явления для прямой задачи гравиразведки можно представить в виде схемы (рис. 2).
Наиболее полной и общей моделью среды, используемой гравиразведкой и позволяющей описывать любые неоднородности, является модель функции плотности как функции пространственных координат, однако она обладает серьёзными недостатками. Объясняется это тем, что по гравиметрическим данным не может быть восстановлено распределение плотности в нижнем полупространстве. Поэтому для изучения строения осадочных бассейнов следует использовать модели структурного типа.
Список литературы
1. Аминов Л. З., Кобрунов А. И., Моисеенкова С. В., Мо-трюк Е.Н., Шилова С. В., Мужикова А. В. Методика интегрированной интерпретации гравиметрических данных в условиях слабой изученности с целью построения объемных региональных плотностных моделей седиментационных бассейнов // Геология и минеральные ресурсы Европейского Северо-востока России: материалы XIV геологического съезда Республики Коми. — Том IV. — Сыктывкар: Геопринт, 2004. — С. 79−81/
2. Кобрунов А. И. Математические методы моделирования в прикладной геофизике (избранные главы) в 2-х частях. Ч.1 Функционально-аналитические основы // Международный журнал экспериментального образования. — 2014. -№ 11. — 78 с. Ч. 2 Системный анализ и моделирование в условиях неопределенности. — Ухта: УГТУ, 2014. — 154 с.
3. Кобрунов А. И. Скрытая эквивалентность и эффективность интерпретации гравиметрических данных // Изв. РАН сер Физика земли. — 2014. — № 2. — С. 53−62
4. Кобрунов А. И., Урбан А. В. О проблеме скрытой эквивалентности при реконструкции моделей геологических сред // Геофизика. — 2009. -№ 3. — С. 41−48.
5. Кобрунов А. И. Математические основы теории интерпретации геофизических данных: учеб. пособие. — Ухта: УГТУ, 2007. — 286 с.: ил.
6. Маловичко А. К., Костицын В. И. Гравиразведка. -М.: 1992. — 357 с.
7. Мотрюк Е. Н. Разработка технологии и методики моделирования глубинного строения земной коры на основе комплексной интерпретации геолого-геофизических данных // Вопросы теории и практики геологической интерпретации геофизических полей: материалы 38 сессии Международного семинара им. Д. Г. Успенского. — Пермь: ГИ УрО РАН, 2011. — С. 201−204.
8. Петровский А. П. Информационное обеспечение и модельные представления интегральной интерпретации геолого-геофизических данных при изучении нефтегазоносных структур // Геофизический журнал. — 2004. — Т.5. -№ 3. — С. 12−13.
References
1. Aminov L.Z., Kobrunov A.I., Moiseenkova S.V., Mo-trjuk E.N., Shilova S.V., Muzhikova A.V. Metodika integri-rovannoj interpretacii gravimetricheskih dannyh v uslovijah slaboj izuchennosti s celju postroenija ob#emnyh regionalnyh plotnostnyh modelej sedimentacionnyh bassejnov // Geologija i mineralnye resursy Evropejskogo Severo-vostoka Rossii: mate-rialy XIV geologicheskogo s#ezda Respubliki Komi. Tom IV. Syktyvkar: Geoprint, 2004. рр. 79−81.
2. Kobrunov A.I. Matematicheskie metody modelirovanija v prikladnoj geofizike (izbrannye glavy) v 2-h chastjah. Ch.1 Funkcionalno-analiticheskie osnovy // Mezhdunarodnyj zhur-nal jeksperimentalnogo obrazovanija. 2014. no. 11. 78 р. Ch. 2 Sistemnyj analiz i modelirovanie v uslovijah neopredelennosti. Uhta: UGTU, 2014. 154 р.
3. Kobrunov A.I. Skrytaja jekvivalentnost i jeffektivnost interpretacii gravimetricheskih dannyh. /Izv. RAN ser Fizika zemli. 2014. no. 2. рр. 53−62
4. Kobrunov A.I., Urban A.V. O probleme skrytoj jek-vivalentnosti pri rekonstrukcii modelej geologicheskih sred // Geofizika. 2009. -no. 3. рр. 41−48.
5. Kobrunov A.I. Matematicheskie osnovy teorii interpre-tacii geofizicheskih dannyh: ucheb. posobie. Uhta: UGTU, 2007. 286 р.: il.
6. Malovichko A.K., Kosticyn V.I. Gravirazvedka. M.: 1992. 357 р.
7. Motrjuk E.N. Razrabotka tehnologii i metodiki mod-elirovanija glubinnogo stroenija zemnoj kory na osnove kom-pleksnoj interpretacii geologo-geofizicheskih dannyh // Voprosy teorii i praktiki geologicheskoj interpretacii geofizicheskih polej: materialy 38 sessii Mezhdunarodnogo seminara im. D.G. Uspen-skogo. Perm: GI UrO RAN, 2011. рр. 201−204.
8. Petrovskij A.P. Informacionnoe obespechenie i model-nye predstavlenija integralnoj interpretacii geologo-geofizich-eskih dannyh pri izuchenii neftegazonosnyh struktur // Geofi-zicheskij zhurnal. 2004. T.5. no. 3. рр. 12−13.
Рецензенты:
Бурмистрова О. Н., д.т.н., заведующая кафедрой технологии и машин лесозаготовок, ФГБОУ ВПО «Ухтинский государственный технический университет», г. Ухта-
Андронов И. Н., д.т.н., заведующий кафедрой сопротивления материалов и деталей машин, ФГБОУ ВПО «Ухтинский государственный технический университет», г. Ухта.

ПоказатьСвернуть
Заполнить форму текущей работой