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

Тип работы:
Реферат
Предмет:
Физико-математические науки


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

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

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

УДК 550. 3:004. 02
ОСОБЕННОСТИ ТЕХНОЛОГИЙ ПОСТРОЕНИЯ МОДЕЛЕЙ СЛОЖНОПОСТРОЕННЫХ ГЕОЛОГИЧЕСКИХ СРЕД ПО ГРАВИМЕТРИЧЕСКИМ ДАННЫМ
1Мотрюк Е.Н., 1Кобрунов А.И., 2Ломинский Д.О.
ФГБОУ ВПО «Ухтинский государственный технический университет», Ухта, e-mail: akobrunov@ugtu. net-
2ФГБОУ ВПО «Российский государственный университет нефти и газа имени И.М. Губкина»,
Москва
Моделирование геологических сред по данным гравиразведки в условиях слабой изученности и сложного строения осадочных бассейнов требует особого подхода. В условиях эквивалентности одно из важных направлений повышения информационной мощности интерпретационных процедур есть развитие технологий извлечения информации из геолого-геофизических данных. Исследования информационной характеристики геологоразведочных работ показали, что в условиях недостаточной изученности сложнопостроенных сред технологические решения должны быть сконцентрированы на активном управлении вводимой дополнительной информацией об исходных данных и информацией, управляющей формированием модельных представлений, а также динамикой их трансформации. Активное использование в технологиях решения обратных задач априорной информации означает вовлечение в процесс решения обратных задач комплекса методов. Комплексирование геофизических методов нашло отражение в интегрированной пассивной и активной интерпретации информации, в результате применения которой построенная модель среды описана согласованными данными различных методов.
Ключевые слова: гравиразведка, сложнопостроенные среды, пассивная и активная комплексная интерпретация, информационная характеристика геологоразведочных работ
FEATURES OF TECHNOLOGY CONSTRUCTION OF MODELS OF COMPLEX GEOLOGICAL ENVIRONMENTS FOR GRAVIMETRIC DATA
1Motryuk E.N., 1Kobrunov A.I., 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
Earth modeling of gravity data in a poor knowledge and the complex structure of sedimentary basins requires a special approach. Under the conditions of equivalence of input models one of the most important ways of increasing the power of information of interpretation procedures is the development of processes to extract information from geological and geophysical data. Research informational characteristics geological exploration works revealed that under the conditions of insufficient study of structurally complex environments technology solutions should be concentrated on the active management of the additional input information on the source data and the information on the formation of model representations as well as the dynamics of their transformation. Active use of technology in the solution of inverse problems of a priori information indicates involvement in the process of solution of inverse problems of the complex methods. Integration of geophysical methods is reflected in the integrated of passive and active integrated interpretation of information, as a result of which the constructed environment model describes the consistent data of various methods.
Keywords: gravity prospecting, complex structure environment, passive and active integrated interpretation, informational characteristic of geological exploration works
Для построения геолого-геофизических моделей осадочных бассейнов и их фрагментов эффективно применение ведущих современных технологий интерпретации гравиметрических данных. Для крупных плотностных структур типа осадочных бассейнов необходимо проводить постановку и решение интерпретационных задач с учётом условий слабой изученности и сложного строения среды. Это позволяет сделать технология, основанная на критериальном подходе к решению обратных задач гравиразведки [1−5, 7]. Она основана на комплексном анализе всей имеющейся гео-
лого-геофизической информации, где для построения отдельной модели среды используются данные других методов.
Рассмотрим подробно, что представляет собой слабая изученность и сложное строение геологической среды. Понятие слабая изученность территорий раскрывается как соотношение между объёмом требуемой информации для решения поставленной геологической задачи и реально имеющейся. Сложность объекта количественно характеризуется объёмом требуемой для его описания информации. Если соотношение требуемой и имеющейся информации
не в пользу первой, то говорят об условиях недостатка информации или, что-то же самое, о слабой изученности данного объекта. Реально имеющаяся информация определена исходными данными, полученными в результатах геологоразведочных работ (ГРР). Если этих данных недостаточно, то говорят об условиях сложного строения, подразумевая соотношение между информационным уровнем, имеющимся в распоряжении геолога, и требуемым. Таким образом, сложность строения и недостаток информации — взаимозависимые и взаимоопределяющие понятия.
Недостаток информации порожден двумя причинами: первая является объективной компонентой, состоящей в нехватке кондиционных данных, вторая характеризуется технико-экономической нецелесообразностью их наращивания за счет проведения дополнительных ГРР. Этот недостаток может быть устранен за счет усовершенствования используемых технологий извлечения информации из имеющегося объёма данных. В этой связи становится очевидной значимость исследований по наращиванию информационной мощности используемых технологий, которая, в свою очередь, имеет и объективное ограничение. В конечном итоге мы приходим к ситуации, когда между информацией, поданной на вход обрабатывающей интерпретационной технологии, и той, в терминах которой описывается требуемый результат, имеется существенный разрыв.
Информационной характеристикой комплекса ГРР [2, 5] служит соотношение объективной информации J заложенной в получаемых данных, объёма информации J извлекаемой из J0 с помощью используемой обрабатывающей интерпретационной технологии, и объёма информации J содержащейся в построенной модели объекта за счёт преобразования данных предшествующего этапа. Для обработки информации используются две технологии. Во-первых, это технология преобразования первичных данных J0 в некоторые подлежащие последующей геологической интерпретации геофизические параметры Jv Во-вторых, это технологии преобразования геофизических параметров, составляющих содержание информации J в некоторую геологическую модель, характеризующуюся объёмом информации J По сути J2 — это тот объем информации, который достигается выполненными с требуемой детальностью построениями.
Величина J1 является функцией J0 и многокомпонентного параметра характеризующего особенности всего цикла извлечения
информации J0. Очевидно, что J0 & gt- J1(J0, ?). Воспользовавшись терминологией, используемой при характеристике информационной модели геофизических исследований, можно утверждать, что переход от J0 к J1 реализуется за счёт процедур геофизической интерпретации — перехода от наблюдаемых к геофизическим параметрам. Переход от информационного уровня J1 к J2 обеспечивается, прежде всего, процедурами геологической интерпретации и состоит в переходе от геофизических к геологическим параметрам. Характеристику перехода технологии можно записать в виде многокомпонентного параметра п: J1(J0, 4) & gt- J2(J1, п). Переход от наблюдаемых к геофизическим параметрам включает в себя процедурную часть — фильтрацию и решение обратных задач и, кроме того, введение некоторой дополнительной к J1, как принято говорить, априорной информации Д/1. Эта часть информации, присутствующая в результатах геофизической интерпретации, может вводиться явно внешним образом либо присутствовать скрыто в параметрах используемых вычислительных процедур. Будем обозначать = А/г В этом случае У'- может превышать J0. Аналогичная ситуация с объёмом информации на втором шаге — построении геологической модели. Наряду с информацией J получаемой за счёт применения технологии п к информации У присутствует компонента Д/2, реализующая дополнение J до требуемого объёма. /'-. Тогда
5 =
dJ2{J,-п)
а/,
ал
¦Д/0+Д/,
+ А/,
Мультипликативная компонента
aj,
dJn
может быть встроена в технологию переработки информации. Аддитивная компонента ДJ2 стоит вне технологий извлечения информации, смысл которой состоит в явном привнесении опыта, интуиции интерпретатора, его неформализованных знаний. Величина Ы2 — это и фантазии интерпретатора, навеянные размышлениями над генезисом изучаемых объектов, соображениями подобия с чем-то ему хорошо знакомым и прочими. Как только они найдут своё обоснование, эта часть тут же войдёт в технологическую мультипликативную компоненту. Поэтому основным направлением развития технологий извлечения информации при изучении сложнопостроенных, слабоизученных сред является минимизация компоненты Д^, развитие технологий, обеспечивающих ак-
тивное управление той частью объективно необходимой дополнительной информации, которая входит в мультипликативную компоненту. В этой связи необходимо определить те методы и принципы, которые могут лежать в основе такой схемы.
Прежде всего, необходимо определить понятие интегрированной интерпретации, которая использует для построения модели среды в переинтерпретированной форме данные других методов. Синонимом термина интегрированная интерпретация, принятого в зарубежных публикациях, в отечественной литературе служит понятие комплексная интерпретация. Рассматривая модели интегрированной интерпретации [2, 5, 8], ограничимся лишь теми из них, которые основаны на принципах решения обратных задач в детерминированной постановке и направлены на реконструкцию физико-геологической модели изучаемой территории. Это соответствует введению понятия обратных задач интегрированной интерпретации как дальнейшему развитию понятия обратных задач геофизики.
В постановке обратных задач комплексной интерпретации выделяются её пассивная и активная формы. Пассивная форма состоит в том, что данные других методов вводятся в процедуры решения обратных задач гравиметрии — построения плотностной модели среды — в уже проинтерпретированной форме. В этом есть как положительные стороны: использование существующих наработок из инструментария анализа данных других методов — так и отрицательные: невозможность оперативного критического пересмотра введённых данных в процессе построения плотностной модели. Активная форма основана на создании методов и технологий совместного решения нескольких различных обратных задач, например совместного решения обратных задач сейсморазведки и гравиразведки. Здесь также есть свои положительные стороны (возможность динамической корректировки результатов нескольких методов) и отрицательные (сложности с введением самого понятия решения и согласованных модельных представлений).
Типичным примером модели пассивной комплексной интерпретации служит схема, состоящая в нахождении оптимальной относительно вводимого критерия плотност-ной модели среди множества других, удовлетворяющих тому же полю.
Рассмотрим особенности получения решения задачи сначала в классе распределения плотности. Вертикальная производная гравитационного потенциала и (у
где V = (V Уo, 2 о) и точка (хо& gt- У0) = Эр г& lt-Е Ео регистрируется на поверхности в Е (р & gt- 0) с уравнением z0) = для плотностной
модели, задана следующим соотношением между мг (у0) и о = о (у):

A (a (v))=u, (s0,z0=y (s0))= 7u (v) (z — \f (s0))dxdydz [(x-x0)2 + (y-y0)2 + (z- |f (s0))2 J
r. (1)
(2)
Весь комплекс априорных данных об искомом распределении плотности может быть выражен в виде принятия критерия оптимальности [1, 2] относительно o (v), имеющего вид требования минимума некоторого функционала J (o (v)):
y (a (v))=|F (o (v)-o°(v)| = = J|f (a (v) — a° (v))|26?v min.
V
Решение задачи (1), (2):
G (v) = a°(v) + F-1F*-1/CK) — (3)
Здесь Z (s0) — функция, параметризующая класс оптимальных решений уравнения (3) — F — линейный замкнутый оператор, имеющий обратный- ц (^0) — мера на поверхности задания поля, учитывающая способ его задания.
Подставив (3) в (1) и решив соответствующее уравнение относительно Z (s0), легко построить искомое оптимальное распределение. Однако есть более простой и вычислительно эффективный путь. Введём итерационный процесс:
G& quot-+1 (v)= G& quot- (v)+ а"^ф& quot-(50), n = 0, 1, 2, …, o°(v) = Oo (v) —

Г 2 ~|3/2 '-
[(x-*0)2 +(y-yor +(z-|/(50)) J
(4)
(5)
ФЛК)=ДО
(z-V (J0))O-(V)A-

[(x-x0)2+(y-yQf+(z- |/O0))2 J2 У
(6)
a" =--
'-0_
(7)
Критерий оптимальности, в определённом смысле эквивалентный (2):
supliT1 [ct (v)-g0 (v)]| min, veF. (8)
Формирование критерия основано на выборе его компонент: аддитивной — нулевого приближения, имеющего смысл модели среды, к которой ищется наилучшее приближение- мультипликативной — оператора К, определяющего смысл понятия минимального уклонения. Этот оператор многокомпонентен, имеет вид свёртки и включает в себя информацию о степени достоверности построения различных компонент нулевого приближения, данные о корреляционной связи между различными компонентами.
В операторной форме приведенная задача, как и множество подобных, может быть записана в виде унифицированной схемы: искомая модель среды хДу) должна соответствовать наблюдаемому полю и1(я), что выражается в операторном равенстве А (х^ = и1(я), и среди множества таких моделей должна обеспечивать оптимальность критерию, выражаемому в виде требования минимума уклонения от заданного фиксированного элемента
г. е. J (Xi — x0)
min.
A (xl) = u1(s), b{x2) = u2{s), — Xj) -& gt- min.
(9)
Решать эту задачу с учётом приближённости задания правых частей можно итерационным процессом. Для унифицикации записи введём обозначения: А = А1- В = А
Тогда процедура решения состоит в следующем:
фТ= А (& lt-)-"-i>-
х*=зГ1-
— + X) Ф2'-
ф 1=А2(хп2)-и2,
(10)
х2 — ,
а =
Ф1а1(х")К,(Х", Х№
Переход к активной форме задачи комплексной интерпретации обеспечивается введением дополнительного требования на модель х состоящую в её соответствии некоторому другому полю и (я): В (х2) = и (я). Нулевое приближение из фиксированного параметра превращается в подбираемый, и сама задача трансформируется в подбор пары распределений, каждое из которых удовлетворяет своему полю, а сами эти распределения — ближайшие друг к другу среди всех иных, допускаемых полями пар:
Здесь К. — должным образом подобранные операторы, обеспечивающие выполнение условия оптимальности в (9).
В приведённой модели активной ком -плексной интерпретации на этапе постановки задачи происходит отказ от попыток построения единой физической модели, соответствующей различным физическим полям. Вместо этого для каждого поля вводится своя модель, а интегрирующим фактором служит требование их наибольшей близости в заданном смысле. Модель активной комплексной интерпретации для двух методов разделяется на взаимоувязанные сопряжённые схемы. Это две схемы пассивной комплексной интерпретации с динамически меняющейся в процессе счета аддитивной и неизменной мультипликативной компонентами критерия оптимальности. В пассивных схемах неизменны и аддитивная, и мультипликативная компоненты. После завершения цикла вычислений и анализа полученного результата параметры модели и критерия оптимальности, включая его аддитивную и мультипликативную части, могут быть изменены и весь процесс повторен, но это уже другой вычислительный цикл. Для вышеизложенных основных принципов активной комплексной интерпретации геофизических данных разработаны вычислительные схемы. Решение обратных задач пассивной комплексной интерпретации реализовано в постоянно совершенствующемся программном ком -плексе ОС18. Применением такого подхода к анализу геофизических данных может служить работа А. П. Петровского и коллектива авторов [6] по Карпатскому региону Украины.
Список литературы
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 obemnyh regionalnyh plotnostnyh modelej sedimentacionnyh bassejnov // Geologija i mineralnye resursy Evropejskogo Severo-vostoka Rossii: ma-terialy XIV geologicheskogo s#ezda Respubliki Komi.- T. IV. Syktyvkar: Geoprint, 2004. pp. 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 p. Ch. 2 Sistemnyj analiz i modelirovanie v uslovijah neopredelennosti. Uhta: UGTU, 2014. 154 p.
3. Kobrunov A.I. Skrytaja jekvivalentnost i jeffektivnost interpretacii gravimetricheskih dannyh // Izv. RAN ser Fizika zemli. 2014. no. 2. pp. 53−62.
4. Kobrunov A.I., Urban A.V. O probleme skrytoj jekviva-lentnosti pri rekonstrukcii modelej geologicheskih sred // Zhur-nal Geofizika. 2009. no. 3. pp. 41−48.
5. Kobrunov A.I. Matematicheskie osnovy teorii interpretacii geofizicheskih dannyh: ucheb. posobie. Uhta: UGTU, 2007. 286 p.: il.
6. Malovichko A.K., Kosticyn V.I. Gravirazvedka. M., 1992. 357 p.
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. Us-penskogo. Perm, 2011. pp. 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. pp. 12−13.
Рецензенты:
Бурмистрова О. Н., д.т.н., заведующая кафедрой технологии и машин лесозаготовок, ФГБОУ ВПО «Ухтинский государственный технический университет», г. Ухта-
Андронов И. Н., д.т.н., заведующий кафедрой сопротивления материалов и деталей машин, ФГБОУ ВПО «Ухтинский государственный технический университет», г. Ухта.

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