Компьютерная технология построения плотностных и магнитных моделей земной коры по данным грави- и магниторазведки

Зиновкин, Сергей Владимирович

Бесплатный автореферат и диссертация по наукам о земле на тему
Компьютерная технология построения плотностных и магнитных моделей земной коры по данным грави- и магниторазведки
ВАК РФ 25.00.10, Геофизика, геофизические методы поисков полезных ископаемых

Автореферат диссертации по теме "Компьютерная технология построения плотностных и магнитных моделей земной коры по данным грави- и магниторазведки"

На правах рукописи УДК 550.834

ЗИНОВКИН СЕРГЕЙ ВЛАДИМИРОВИЧ

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

Специальность 25.00.10 Геофизика, геофизические методы поисков полезных ископаемых

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

МОСКВА 2011

4850789

Работа выполнена на кафедре геофизики Российского Государственного Геологоразведочного Университета им. Серго Оржоникидзе (РГГРУ)

Научный руководитель: доктор физико-математических наук,

профессор Петров Алексей Владимирович

Официальные оппоненты:

доктор технических наук Галуев

Владимир Иванович

(ФГУП ГНЦ РФ ВНИИГеосистем)

кандидат технических наук Трусов Алексей Андреевич (ЗАО «ГНПП Аэрогеофизика»)

Ведущая организация:

ФГУНПП "Аэрогеология"

Защита состоится « 19 » мая в _1_5 час. 00 мин. на заседании диссертационного совета Д 212.121.07 при Российском Государственном Геологоразведочном Университете по адресу: 117997, Москва, ул. Миклухо-Маклая, 23, РГГРУ, ауд. 6-38.

С диссертацией можно ознакомиться в научной библиотеке Российского Государственного Геологоразведочного Университета.

Автореферат разослан « 15 » апреля 2011г.

Ученый секретарь диссертационного совета Д 212.121.07, доктор

физико-математических наук, профессор А.Д. Карински"

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

Построение плотностных и магнитных моделей земной коры представляет особую актуальность в связи с тем, что данные гравиразведки и магниторазведки в масштабе 1:200 ООО получены практически по всей территории Российской Федерации, а для ряда площадей получены в масштабе 1:50 000. Кроме этого, плотностные модели представляют основу прогноза поисковых исследований на углеводородное сырье и твердые полезные ископаемые при изучении осадочного чехла и кристаллического фундамента. Эффективное решение задачи расчленения пород кристаллического фундамента невозможно без оценки распределения магнитной восприимчивости с глубиной.

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

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

Цели и задачи исследования. Создание компьютерной технологии

построения плотностных и магнитных моделей земной коры в двухмерной (2Б) и трехмерной пространствах (21АТ) и ЗЭ) является основной целью исследований. Достижение цели базируется на решении следующих задач:

- построение оптимальной структуры базы данных для процесса моделирования и соответствующей программной реализации;

- разработка оригинальных алгоритмов редактирования геометрических и физических характеристик модели;

- разработка приемов решения прямых задач и методов оценки параметров анамалиеобразующих объектов на основе сверточной модели в пространстве и в спектральной области;

- использование методов кластерного анализа и распознавания образов в процессе моделирования.

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

Научная новизна исследований определяется:

- разработкой оригинальных алгоритмов решения прямой задачи грави- и магниторазведки на основе сверточной модели для точечных масс;

- включением в процесс построения плотностных и магнитных моделей методов автоматической и интерактивной классификации;

- программной реализацией процедуры построения плотностных и магнитных моделей;

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

Защищаемые положения:

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

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

создания согласованных физико-геологических моделей земной коры.

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

Личный вклад. Все положения, выносимые на защиту, выполнены автором или при его непосредственном участии. Автором проведены исследования по разработке оригинальных алгоритмов решения прямых задач гравиразведки и магниторазведки для сеточных моделей. Разработана методика построения плотностных и магнитных моделей земной коры. Создано программное обеспечение для редактирования модельных построений. Разработан оригинальный интерфейс и база данных для сопровождения модельных построений и визуализации результатов моделирования.

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

Практическая ценность работы состоит в создании компьютерной технологии по созданию плотностных и магнитных моделей, обеспечивающей обработку данных грави- и магниторазведки в режиме реального времени и внедрением программно-алгоритмического обеспечения во ВСЕГЕИ и РГГРУ.

Апробация работы. Основные результаты проводимых исследований, изложенных в работе, докладывались на 33-й, 34-й, 35-й и 37-й сессиях Международного семинара им. Д.Г. Успенского (2006-2008 и 2010гт.), X Международной конференции Новые идеи в науках о Земле (Москва, 2011г.). По теме диссертации опубликовано 16 печатных работ, в том числе 3 в реферируемых журналах.

Структура и объем диссертации. Диссертация состоит из введения, четырех глав и заключения. Общий объем составляет 122 страниц, в том числе 45 рисунков и 2 таблицы. Список литературы включает 78 наименований.

Благодарности. Автор выражает глубокую благодарность научному руководителю, доктору физико-математических наук, профессору Петрову Алексею Владимировичу и доктору физико-математических, заведующему кафедры геофизики, профессору Никитину Алексею Алексеевичу за внимание, помощь и поддержку, советы и ценные замечания, оказываемые за годы совместной работы; специалистам ВСЕГЕИ Мильштейн Е.Д. и Каличевой Т.И.

за консультации и конструктивную критику; сотрудникам геофизического факультета РГГРУ за моральную поддержку.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Методика обработки гравитационных и магнитных полей в настоящее время достаточно проработана и, по существу, сведена к стандартному графу [Никитин A.A., 2007]. При этом учитывается существенно нестационарный характер потенциальных полей, когда среднее значение поля, его дисперсия и корреляционные (спектральные) свойства претерпевают изменения, как вдоль одиночных профилей, так и по площади съемки.

Для оценки изменения статистических и градиентных атрибутов по глубине разреза земной коры A.B. Петровым предложена методика их статистического зондирования. Результатом такого зондирования является создание трехмерных сетей, состоящих из отдельных слоев, полученных путем расчета того или иного атрибута при конкретном размере скользящего окна.

Количественная интерпретация данных гравитационных и магнитных полей осуществляется с целью получения распределения в плоскости разреза земной коры эффективных значений плотности и намагниченности. Для этого также используется построение их объемных моделей. Трехмерные модели распределения эффективных значений плотности и намагниченности строятся путем вычисления трансформант гравитационного и магнитного полей, по алгоритмам И.И. Приезжаева (по пересчету полей в нижнее полупространство в спектральной области) [Приезжев И.И., 2005], Ю.И. Блоха (построение 3D-распределение эффективных источников магнитных масс) [Блох Ю.И., 2003], A.B. Петрова (модифицированный метод Б.А. Андреева) [Петров A.B., 2007].

Приведенными выше моделями описываются распределения избыточных физических свойств изучаемой геосреды в земной коре. Для получения непосредственно плотностной модели земной коры проводится нормирование эффективных значений плотности и построение плотностного разреза с учетом изменения плотности с глубиной. Обычно при этом задается линейный тренд плотностной характеристики среды при увеличении плотности

с глубиной [Галуев В.И. 2009].

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

В работе приведен обзор компьютерных систем построения согласованных моделей по геофизическим данным: ГИС ИНТЕГРО ГЕОФИЗИКА, вСК, СИГМА-ЗО, 0\МЗео; а также монометодных систем интерактивного построения гравимагнитных моделей ОеоБой вМ-БУБ, В1ШМ ЗБ СеоМос1е11ег и КОСКАД ЗБ.

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

Рассмотрим интегральное представление потенциала гравитационных масс, распределенных в области V:

Эта формула является основной при решении прямых задач для различных трехмерных геологических объектов. Дифференцируя ее можно получать интегральные соотношения для любых элементов гравитационного поля. Наибольшее практическое применение при интерпретации гравитационных аномалий имеет аномалия ускорения силы тяжести Лg, которая отождествляется с их\

Возможности современной вычислительной техники позволяют отказаться от представления модели в виде набора многоугольников или призм, позволяющих сократить время решения прямых задач, но усложняющих процедуру редактирования моделей, и перейти к более эффективным и удобным для редактирования - равномерным сеточным моделям. В этом случае исследуемое пространство разбивается на элементарные объемы, каждый из которых наделяется плотностными свойствами. В общем случае требуется решение прямой задачи от источников, расположенных в узлах трехмерной регулярной сети. В качестве источников можно использовать точечные источники, стержни или прямоугольные параллелепипеды. Для модели на площади 100x100км с шагом 0.1км и 10км в глубину требуется =1014 кратное

решение прямой задачи (в этом случае число пикетов и=1001, число профилей ш=1001, число слоев /г.г=101):

Д^ЁЕЁТ--j,i = l...nJ = L..m.

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

Значительный прирост скорости можно получить за счет послойного решения прямой задачи. Для источников, расположенных на одной глубине z0, форма аномалии одинаковая, и следовательно, гравитационный эффект слоя A?(x>jO|z=«rei можно выразить в виде свертки плотностей a(x,y,z0) этого слоя и поля точечного источника K(x,y,z-z0), расположенного на глубине слоя: &g(x, у) | z=consl = а(х, y,z0)* K(x,y,z - z0) = Ц cr(x,y,z0)K(x - t,y - q,z - z0 )dtdq.

Для регулярной сети, содержащей т профилей по п пикетов, получим:

т п

bgi,j\z=COns, = о).

р=\г=\

где i = \...m, j = \...п. При этом количество сложных вычислительных операций (расчет матриц К), требующих значительного машинного времени, сокращается до ~108 (для всего куба), а количество простых операций -сложение и умножение остается «1014. Схема вычислений проиллюстрирована на рис. 1.

Еще больший прирост скорости можно получить, если выполнять свертки в частотной области. Здесь потребуется «3*102 быстрых двумерных преобразований Фурье и »108 операций умножения.

Прямые задачи в спектральной области. Прямая задача в разведочной геофизике в общей постановке для статических полей на плоскости z0=const определяется следующим образом:

<р(хо ,JVzo = const) = Р JJJ а(х, у, z)k(x-x0,y-y0,z-z0 )dxdydz,

где <p(x0,y0,zQ = const) - наблюдаемая функция физического поля;

a(x,y,z) - функция физического параметра горных пород, определенная в

нижнем полупространстве;

K(x-x0,y-y0,z-z0) - ядро интеграла или функция, связывающая распределение физического параметра среды и наблюдаемое поле, или функция реакции среды на единичный импульс (поле от точечного источника);

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

Рис. 1. Схема вычислений послойного решения прямой задачи путем свертки плотности с полем точечного источника. Применим преобразование Фурье по осям X,Y и на основе теоремы о свертке получим:

®(w„w2,z0 = const) = JC(vv[, w2,z)K(w,,w2,z- z0)dz,

2=0

где W], w2 имеют смысл волновых чисел (пространственных или волночисловых частот) и именуемых ниже как частота.

Используя преобразование Ханкеля Н0{ }, можно убедиться, что двумерное преобразование ядра для гравитационного поля Ag имеет вид:

// J_Z~Z0_I = g-w(z-Zo),

_1 [(* -xof + (у ~Уо)2+iz~ zof 1')

где w = д/й),2 + со22 - радиальная частота.

Используя теоремы о производной и первообразной, формулы, приведенные Страховым В.Н. для вычисления спектра сопряженных гармонических функций, а также формулы, приведенные Серкеровым С.А. для вычисления спектра косо намагниченных объектов и AT, можно определить

спектр основных в геофизике ядер для статических полей в виде:

K(Wl ,w2,z - z0)= M(wi,w2)e~wi2'Zo), где M{ú\,co2) - в общем случае комплексная функция частоты, определяющая тип ядра. В работе приведена сводка видов функции для различных типов ядра [Приезжев И.И., 2005].

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

Программное обеспечение для формирования плотностных и магнитных моделей функционирует с индексно-последовательной, базой данных компьютерной технологии КОСКАД-ПРОФИЛЬ. Эта технология ориентирована на анализ и обработку профильных данных, с использованием как хорошо зарекомендовавших себя алгоритмов, составляющих функциональное наполнение технологии КОСКАД 3D, так и новых оригинальных алгоритмов.

Исходной информацией для технологии КОСКАД-ПРОФИЛЬ является сеть наблюдений, которая соответствует площади исследований, ограниченной прямоугольными координатами Хмин, Хмакс, YMUH, YMaKC, и содержит определенное число профилей. При этом каждый из профилей сети характеризуется количеством точек наблюдения, координатами X, Y, Z каждой точки профиля и значениями набора исследуемых признаков. Если по какому-либо признаку сети на отдельных профилях информация отсутствует, то соответствующие поля в БД заполняются кодом отсутствия измерения. Кроме этого, все профиля сети могут содержать несколько слоев, расположенных друг от друга на постоянном расстоянии dZ параллельно поверхности рельефа. Расстояние между соседними слоями dZ, как и количество слоев, является постоянным для всех профилей сети. Введение понятия слоя позволяет эффективно решать задачи трехмерного моделирования, при этом постоянство расстояния между отдельными слоями повышает эффективность решения большого числа задач, возникающих в процессе моделирования. В отдельно выбранной директории можно хранить 99 различных профильных сетей [Петров A.B., Юдин Д.Б, 2010].

В технологию гравитационного и магнитного моделирования входят:

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

- программы интерактивной и автоматической классификации;

- программы решения прямых 3D задач грави- и магниторазведки.

При моделировании плотностных и магнитных разрезов модель представляет собой совокупность точечных источников, организованную в двухмерную регулярную сеть. По горизонтальной оси откладываются координаты профиля, а по вертикальной (направленной вниз) глубины. Каждый узел сети наделяется физическими свойствами плотности или магнитной восприимчивости. Эти свойства далее будем назвать массами (гравитационные массы или магнитные массы).

Основными возможностями системы являются:

- моделирование слоистых сред с изменением физических свойств по латерали, в том числе и моделирование разломов;

- моделирование блоковых сред и геологических тел;

- расчет прямых гравитационных и магнитных эффектов модели и сопоставление наблюденного поля с результатом решения прямой задачи;

- интерактивный процесс построения модели с широкими возможностями графического интерфейса пользователя.

Для формирования плотностного или магнитного разреза слоистой среды используется набор полилиний. Полилиния - это заданная в плоскости разреза ломаная линия, в каждом узле которой определены "массы слева" и "массы справа". Задание различных масс слева и справа дает возможность моделировать разрывные нарушения слоистых сред (разломы). Моделирование осуществляется следующим образом:

1. Узлы сети модели очищаются (заполняются кодом "пропуска измерения").

2. В плоскости разреза задается одна или несколько ломаных линий - полилиний, в узлах которых задаются магнитные или гравитационные "массы слева" и "массы справа" (эти массы могут и совпадать). Вдоль ломаной эти массы интерполируются при помощи линейной интерполяции.

3. Массы, каждой полилинии, переносятся в сеть модели, причем "массы слева" переносятся в узлы, лежащие слева от ломаной линии, а "массы справа" переносятся в узлы, лежащие справа.

4. Для каждого слоя модели Zj=Zo+(j-l)dZ, интерполяции, определены.

Для формирования плотностного или магнитного разреза блоковой

261 2 60 259 2.58 2.57 2.55 255 2 54 2.53 2.52 2.51 2.50 2.49 2.48 2.47 2 46 2 45 2 44 2.43 2.42

2 37 2.36 2.35 2.34 2.33 2 32 2.31 2 30

при помощи линеинои производится восполнение узлов, массы в которых не

среды используется набор примыкающих друг к другу, горизонтально расположенных цилиндрических тел, сечениями которых являются многоугольники или полигоны. Простирание цилиндров перпендикулярно плоскости разреза. Для двухмерных моделей простирание цилиндров не ограничено, а в 2'АО случае простирание каждого из цилиндров ограничено плоскостями у=у° и у=у". В пределах каждого тела задаются значения масс, которые могут быть постоянными или изменятся по линейному закону.

Полигон - это заданный в плоскости разреза многоугольник, который является основанием (сечением) цилиндрического тела. Для каждого полигона задаются местоположение и размер по оси У - у" и у" и массы в двух точах. Масса каждого тела изменяется линейно относительно координатных осей плоскости разреза: М(х,г)=Ах+Вг+С. Константы А, В, С вычисляются автоматически в зависимости от расположения в плоскости разреза пары точек (на рис. 2 точки, соединенные отрезком) и соответствующих им масс, которые задаются в таблице свойств. Массы каждого полигона переносятся на сеть модели, заменяя массы слоистой модели, или суммируются с ними (если для полигона в контекстном меню установлено свойство "Относительная плотность"). Каждый полигон может быть включен в процесс моделирования

* | Обрабатывать

Относительная плотность

X 2 Э1д2 Х2. 22

19.7898312 29 2 92

41 60165 20 2452431 46 25572 69 3835754

89.66201 21.4883728 22 775898 22.775898

92 51144 25 2177582

67 2465057 25 839323

28 6842442 25 3065548

0 24 68499

► У0

Уп 20

Рис. 2. Моделирование геологических тел.

Учет свойств вмещающей среды дает возможность задавать объектам моделирования (полигонам и полилиниям) абсолютные (реальные) гравитационные и магнитные массы. Вместе с тем, если свойства вмещающей среды не заданы, то они считаются нулевыми. В этом случае объектам моделирования следует задавать избыточные (эффективные) массы. Свойства вмещающей среды определяются для каждого горизонтального слоя модели Zj=const при помощи таблицы, состоящей из двух колонок. В первой задаются глубины, а во второй соответствующие им массы.

Расчеты прямых задач производятся со следующими основными параметрами:

• тип поля - гравитационное или магнитное;

• постановка задачи - двухмерная или трехмерная;

• способ решения прямой задачи - в пространстве или в спектральной области;

• углы нормального магнитного склонения И и наклонения I;

• высота пересчета прямой задачи.

• режим восполнения узлов модели с неопределенными физическими свойствами. Восполнение производится методом ближайшего соседа;

• коэффициенты линейного преобразования рассчитанного поля для увязки с наблюденным полем.

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

К дополнительным возможностям системы относятся:

1. Построение графиков наблюденного и рассчитанного полей, и их разности.

2. Отображение в плоскости разреза подложки.

3. Настройка параметров визуализации сети и векторных объектов.

4. Масштабирование изображения, просмотр разреза с нарушением пропорций.

5. Отмена последних действий и повтор ранее отмененных.

6. Режим вывода эффективных плотностей модели.

7. Другие возможности (вывести сетку и узлы, выбрать масштаб графика разности и пр.).

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

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

- Запускается вычислительный процесс, состоящий в послойном приведении исходной сети к регулярной и вычислении прямого гравитационного или магнитного эффекта последней.

- В зависимости от выбранных параметров регулярная сеть может быть сохранена для дальнейшего использования в комплексе КОСКАД ЗЭ.

- Завершаются вычисления интерполяцией результатов решения прямой задачи в точки исходной профильной сети и сохранением результата в заданную сеть.

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

Программа интерактивной классификации многопризнаковых данных КЛаэБ входит в модуль "комплекс" компьютерной технологии КОСКАД ЗБ и предназначена для проведения классификаций в интерактивном режиме. Программа позволяет:

■ Организовать рабочую область (проект), содержащую исходные, промежуточные и выходные данные.

■ Проводить классификацию в интерактивном режиме.

■ Осуществлять совместный просмотр исходных данных и результата классификации в виде наложенных друг на друга карт изолиний с широкими сервисными возможностями.

При работе с программой предлагается следующая схема действий:

1. Определить исходные данные - признаковое пространство.

2. Задать начальное количество классов.

3. Разбить Л'-мерное признаковое пространство на непересекающиеся № мерные полуинтервалы.

4. Задать соответствие между каждым классом и элементами разбиения признакового пространства.

5. Построить карту классификации.

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

7. Для точек, не вошедших ни в один класс, определить классовую принадлежность. Построить карту классификации.

8. При необходимости вернуться к п. 6.

9. Сохранить результат классификации.

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

Методика построения плотностных разрезов при создании ЗБ моделей земной коры предполагает следующие этапы:

1. Этапы моделирования, реализуемые в программном комплексе спектрально-корреляционного анализа КОСКАД ЗЭ [Мильштейн Е.Д., 2007]:

■ Получение трехмерной модели относительного распределения плотностей методом Андреева.

■ Корректировка глубин источников, уточнение шкалы глубин.

■ Нормирование, переход к модели эффективной плотности.

■ Получение плотностной модели путем преобразования, обеспечивающего градиентно-слоистое нарастание плотности.

2. Подготовка данных для программы моделирования комплекса КОСКАД-ПРОФИЛЬ.

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

■ Векторизация плотностного разреза. Векторизация - это преобразование растровых изображений в элементы векторной графики (точки, линии, полигоны, надписи и др.). Плотностной разрез и другая априорная информация позволяет интерпретатору синтезировать геометрию модели, разбить разрез на слои и блоки. Эти действия могут быть выполнены в программе моделирования вручную. Далее рассмотрим интерактивный, полуавтоматический способ геометризации модели:

- В системе КОСКАД 30 из плотностного куба получаем разрез вдоль профиля.

- Выполняем кластеризацию полученного разреза. Для этого можно использовать программы классификации модуля "КОМПЛЕКС" системы КОСКАД 3D - метод динамических сгущений, классификация по Петрову A.B. и интерактивная классификация. Строим карту классификации и сохраняем изображение.

- Векторизуем полученное изображение. Для этого можно воспользоваться утилитой PowerTRACE графического приложения CorelDRAW. В результате получим набор полигонов.

- Выполняем проецирование полигонов в координаты разреза и сохранение векторной информации в формате, пригодном для системы моделирования. Для этого можно воспользоваться специальной макро-процедурой, разработанной автором для графического приложения CorelDRAW или приложением Golden Software Surfer.

3. Основные этапы, выполняемые в предложенной системе моделирования.

• Подготовка таблицы для учета свойств вмещающей среды.

■ Подбор геометрических и физических характеристик профильных моделей.

■ Построение трехмерной модели на основе послойной интерполяции с использованием методов геостатистики.

В четвертой главе представлены примеры практического использовании предлагаемой технологии моделирования. На рис. 3 представлена модель плотностного разреза по фрагменту опорного регионального профиля 1-ЕВ (ПК 850-1200).

На рис. 4 представлена модель плотностного разреза в осадочном чехле, построенная на основе геологического разреза (показанного в нижней части рисунка), поученного при проведении геологических изысканий на уран в пределах Урюмкано-Уровского района Забайкалья.

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

Рис. 3. Онежский участок регионального опорного профиля!ЕВ (ПК 850-1200).

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

ОСНОВНЫЕ ВЫВОДЫ И РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ

■ создание программно-алгоритмических решений и математического аппарата для оценки адекватности и увязки модельных построений по данным гравиразведки и магниторазведки;

■ дальнейшее совершенствование технологии создание магнитных моделей, учитывающей весь спектр свойств объектов, индуцирующих магнитное поле;

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

ОПУБЛИКОВАННЫЕ РАБОТЫ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Решение обратной линейной задачи гравиразведки на основе метода регуляризации для системы многоугольников. Соавторы Филатов В.Г., Гласко Ю.В., Овсепян М.Л., Сташевский В.Е. Геофизика №5, 2005, с. 58.

2. Адаптивные приемы выделения неоднородных геологических объектов в геофизических полях. Соавторы Никитин A.A., Петров A.B. Изв. ВУЗов. Геология и разведка №3,2006г. с. 50 - 56.

3. Развитие статистических приемов обработки и интерпретации геофизических полей в компьютерной технологии КОСКАД-ЗД. Соавторы Никитин A.A., Петров A.B. Изв. ВУЗов. Геология и разведка №6, 2007г.

4. Применение регуляризации и оптимальной фильтрации геофизических данных при поисках месторождений углеводородов (Глава 2). Соавторы Никитин A.A., Петров A.B., Мегеря В.М., Филатов В.Г., Лобанов A.M.. -Учебное пособие для ВУЗов, М.: Изд. РГГРУ, 2011, с. 29-78.

5. Оптимальная фильтрация и интропродолжение геополей с учетом вторичного магнитоминералообразования в нефтегазоразведке (Главы 3,4). Соавторы Никитин A.A., Петров A.B., Мегеря В.М., Старостенко В.И., Филатов В.Г., Лобанов A.M.. - Учебное пособие для ВУЗов, М.: Изд. РГГРУ, 2011, с. 48-98.

6. Интропродолжение и эпигенетическое магнитоминералообразование в нефтегазоразведке. Соавторы Лобанов A.M., Филатов В.Г., Петров A.B. — Учебное пособие для ВУЗов, М.: Изд. РГГРУ, 2009, 79с.

7. Компьютерная система геолого-математического моделирования рудных тел. Соавторы Сикорский В.А. - VI Международная конференция "Новые идеи в науках о Земле", Москва, РГГРУ, 2003, том 4, с. 229.

8. Возможности компьютерной технологии "Коскад 3Dt" при обработке и интерпретации данных глубинной сейсморазведки. Соавторы Петров A.B., Пискун П.В. - Тезисы докладов на V международной конференции "Крым 2003", Крым, Гурзуф, 8-13 сентября 2003г., с. 285.

9. Решение обратных задач гравимагниторазведки на основе методов регуляризации, оптимальной фильтрации и полного нормированного градиента при поисках месторождений нефти и газа. Соавторы Филатов

1. В.Г., Овсепян M.JI., Булычев Е.В., Александров Г.В., Пискун П.В. -Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей: Материалы 31-й сессии Международного семинара им. Д.Г. Успенского, Москва, 26-29 января 2004г. - М.: ОИФЗ РАН, 2004. - 86 с.

2. Компьютерная технология "КОСКАД СЕЙС" комплексного анализа сейсмической и скважинной информации. Соавторы Петров A.B., Пискун П.В. - Тезисы докладов VII-ой международной научно-практической конференции Геомодель-2005 г. Геленжик, 11-17 сентября 2005г. 144 с.

3. Возможности моделирования по гравимагнитным и сейсмическим данным в компьютерной технологии «КОСКАД 3D». - Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей: Материалы 33-й сессии Международного семинара им. Д.Г. Успенского, Екатеринбург, 30 января - 3 февраля 2006г. -Екатеринбург.: Институт геофизики УрО РАН, 2006. - 467 с.

4. Комплекс регуляризированных алгоритмов разделения потенциальных полей и решения обратных задач гравимагниторазведки при поисках нефтегазовых месторождений. Соавторы Филатов В.Г. и др. VII Международная конференция Новые идеи в науках о Земле, Москва, РГГРУ, 2005, с. 26.

5. Развитие компьютерной системы КОСКАД-ЗД. Соавтор Петров A.B. -Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей: Материалы 34-й сессии Международного семинара им. Д.Г. Успенского, Москва, 29 января - 3 февраля 2007г. - М.: ОИФЗ РАН, 2007. - 301 с.

6. Развитие статистических приемов обработки и интерпретации геофизических полей в компьютерной технологии КОСКАД-ЗД. Соавтор Петров A.B. Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей: материалы 35-й сессии Международного семинара им. Д.Г. Успенского, г. Ухта УГТУ, 2008г.-с. 241-242.

7. Компьютерная технология построения гравимагнитных моделей. Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей: Материалы 37-й сессии Международного семинара им. Д.Г. Успенского, Москва, 25 января — 29 января 2010 г., с. 149.

8. Возможности моделирования гравимагнитных полей в компьютерной системе КОСКАД-ПРОФИЛЬ. - X Международная конференция Новые идеи в науках о Земле, Москва, РГГРУ, 2011, Том 2, с. 25.

Подпись автора у!_/ C.B. Зиновкин /

Формат 60x90/16. Заказ 1405. Тираж 100 экз. Подписано в печать 14.04.2011 г.

Печать офсетная. Бумага для множительных аппаратов.

Отпечатано в ООО "ФЭД+", Москва, ул. Кедрова, д. 15, тел. 774-26-96

Содержание диссертации, кандидата технических наук, Зиновкин, Сергей Владимирович

ВВЕДЕНИЕ

ГЛАВА I. СОВРЕМЕННЫЕ КОМПЬЮТЕРНЫЕ ТЕХНОЛОГИИ СОЗДАНИЯ ФИЗИКО-ГЕОЛОГИЧЕСКИХ МОДЕЛЕЙ ЗЕМНОЙ КОРЫ

1.1. Компьютерные технологии построения согласованных физико-геологических моделей

1.2. Компьютерные системы построения несогласованных моделей

ГЛАВА II. АЛГОРИТМИЧЕСКОЕ ОБЕСПЕЧЕНИЕ КОМПЬЮТЕРНОЙ ТЕХНОЛОГИИ ПОСТРОЕНИЯ ПЛОТНОСТНЫХ И МАГНИТНЫХ МОДЕЛЕЙ

2.1. Способы решения прямых задач гравиразведки

2.2. Алгоритмы решения прямых и обратных задач в спектральной области

2.3. Модифицированный метод Б.А. Андреева

2.4. Классификационные алгоритмы

2.5. Интерполяционные алгоритмы

ГЛАВА III. ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ ТЕХНОЛОГИИ МОДЕЛИРОВАНИЯ ПЛОТНОСТНЫХ И МАГНИТНЫХ РАЗРЕЗОВ

3.1. Интерфейс компьютерной технологии гравитационного и магнитного моделирования

3.2. Создание проекта

3.3. Методика построения моделей плотностных и магнитных разрезов слоистой среды

3.4. Методика построения моделей плотностных и магнитных разрезов блоковой среды

3.5. Учет свойств вмещающей среды

3.6. Расчет прямой задачи

3.7. Сохранение результатов

3.8. Дополнительные возможности программного комплекса

3.9. Программы решения прямых задач

3.10. Программы интерактивной классификации

3.11. Методика построения моделей

ГЛАВА IV. ПРИМЕРЫ ПРАКТИЧЕСКОГО ИСПОЛЬЗОВАНИЯ ТЕХНОЛОГИИ ПОСТРОЕНИЯ ПЛОТНОСТНЫХ И МАГНИТНЫХ МОДЕЛЕЙ ЗЕМНОЙ КОРЫ

4.1. Плотностные модели по опорным региональным профилям

4.2. Пример построения плотностных моделей по геологическим разрезам

Введение Диссертация по наукам о земле, на тему "Компьютерная технология построения плотностных и магнитных моделей земной коры по данным грави- и магниторазведки"

Построение плотностных и магнитных моделей разрезов земной коры представляет особую актуальность в связи с тем, что данные гравиразведки и магниторазведки в масштабе 1:200 ООО получены практически по всей территории Российской Федерации, а для ряда площадей получены в масштабе 1:50 000. Кроме этого, плотностные модели представляют основу прогноза поисковых исследований на углеводородное сырье и твердые полезные ископаемые при изучении осадочного чехла и кристаллического фундамента. Эффективное решение задачи расчленения пород кристаллического фундамента невозможно без оценки распределения магнитной восприимчивости с глубиной.

Цели и задачи исследования. Создание компьютерной технологии' построения плотностных и магнитных моделей земной коры в двухмерной (2X5) и трехмерной пространствах (2УШ и ЗО), является основной целью исследований.

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

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

Научная новизна исследований определяется: разработкой оригинальных алгоритмов решения прямой задачи грави- и магниторазведки на основе сверточной модели для точечных масс; включением в процесс построения плотностных и магнитных моделей методов автоматической и интерактивной классификации; программной реализацией процедуры построения плотностных и магнитных моделей; созданием плотностных моделей глубинного строения земной коры по опорным региональным профилям, адекватных реальным средам.

Защищаемые положения:

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

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

Практическая г^енностъ работы состоит в создании компьютерной технологии по созданию плотностных и магнитных моделей, обеспечивающей обработку данных грави- и магниторазведки в режиме реального времени и внедрением программно-алгоритмического обеспечения во ВСЕГЕИ и РГГРУ.

Заключение Диссертация по теме "Геофизика, геофизические методы поисков полезных ископаемых", Зиновкин, Сергей Владимирович

ЗАКЛЮЧЕНИЕ

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

Библиография Диссертация по наукам о земле, кандидата технических наук, Зиновкин, Сергей Владимирович, Москва

1. Андреев Б А., Закашанский М.С. и др. Курс гравитационной разведки. М.:; Л.: Геолитиздат, 1941, 432с.

2. Андреев Б.А., Клушин И.Г. Геологическое истолкование гравитационных аномалий. Ленинград, ГОСТОПТЕХИЗДАТ, 1962, 496с.

3. Аронов В.И. Трехмерная аппроксимация как проблема обработки, моделирования и интерпретации геофизических и геологических данных. // Геофизика, 2000, №4. с.21-25.

4. Бабаянц П.С., Блох Ю.И., Трусов A.A. Изучение рельефа кристаллического основания платформенных областей по данным магниторазведки и гравиразведки. Геофизика, 2003, №6, с.55-58.

5. Белобородов М.А., Боков K.M. и др. ГИС ПАРК (версия 10.0): руководство пользователя. -М.: МФ ВСЕГЕИ, 2009, 333с.

6. Березкин В.М. Метод полного градиента при геофизической разведке. — М.: Недра, 1988, 188с.

7. Березкин В.М., Жбанков Ю.В.,Филатов В.Г., Трайнин П.Н., Булычев Е.В. Методические рекомендации по технологии обработки и интерпретации геофизических данных (ТПОИГД). -М.: изд. Нефтегеофизики, 1992, 80с.

8. Билибин С.И., Перепечкин М.В., Юканова Е.А. Технология построения геологических моделей залежей углеводородов в программном комплексе DV-Geo //Геофизика, 2007, №4, с. 191-194.

9. Блох Ю.И. Интерпретация гравитационных и магнитных аномалий. Учебное пособие. М., 2009. http://sigma3d.eom/content/view/24/2/.

10. Блох Ю.И. Решение прямых задач гравиразведки и магниторазведки. Учебное пособие. -М.: ММГА, 1993, 79с.

11. Булах Е.Г. Прямые и обратные задачи гравиметрии и магнитометрии. Математические методы геологической интерпретации гравиметрических и магнитометрических данных. — Изд. "Наукова думка", Киев, 2010, 464с.

12. Бычков С.Г. Методы обработки и интерпретации гравиметрических наблюдений при решении задач нефтегазовой геологии. Екатеринбург: УрО РАН, 2010, 187с.

13. Вахромеев Г.С., Давыденко А.Ю. Моделирование в разведочной геофизике. -М.: Недра, 1987, 192с.

14. Галуев В.И. ГИС ИНТЕГРО-ГЕОФИЗИКА геоинформационная система интегрированной интерпретации геофизических данных для изучения глубинного строения земной коры//Геоинформатика, 2006, №1, с. 3-9.

15. Галуев В.И., Каплан С.А., Никитин A.A. Технология создания физико-геологических моделей земной коры по опорным профилям на основе геоинформационных систем. — М.: ВНИИгеосистем, 2009, 236с.

16. Геоинформационная система ПАРК (версия 6.01). Введение в систему и технологию, ЛАНЭКО, 2000, 98с.

17. Геологическое строение и закономерности размещения полезных ископаемых. Сибирская платформа. Т. 4 // JL: Недра, 1987. 448 с.

18. Голиздра Г.Я. Комплексная интерпретация геофизических полей при изучении глубинного строения земной коры. — М.: Недра, 1988, 212с.

19. Гольцман Ф.М., Калинин Д.Ф., Калинина Т.Б. Компьютерная технология Mult Alt альтернативной классификации и прогноза по комплексу геоданных. // Российский геофизический журнал, 2000, № 17-18, с. 64-70.

20. Гравиразведка. Справочник геофизика. / Под ред. Е.А. Мудрецовой, К.Е. Веселова, 2-е изд. -М.: Недра, 1990, 607с.

21. Дуда Р., Харт П. Распознавание образов и анализ сцен: Пер. с англ. — М.: Мир, 1976,511с.

22. Зиновкин C.B., Сикорский В.А. Компьютерная система геолого-математического моделирования рудных тел. VI Международная конференция Новые идеи в науках о Земле 2003г., Москва РГГРУ, 2003, том 4, с. 229.

23. Зиновкин C.B., Филатов В.Г., Гласко Ю.В., Овсепян М.Л., Сташевский В.Е. Решение обратной линейной задачи гравиразведки на основе метода регуляризации для системы многоугольников. Геофизика №5, 2005, с. 58.

24. Зиновкин C.B. Интропродолжение и эпигенетическое магнито-минералообразование в нефтегазоразведке: Учеб. пособие для ВУЗов / Соавторы Лобанов A.M., Филатов В.Г., Петров A.B. и др., М.: РГГРУ, 2009, 79с.

25. Игнатьева Т.С. Петрофизическая характеристика разреза среднего протерозоя Онежского полуострова. / В кн.: Геофизические ипетрофизические исследования в Карелии, Петрозаводск, Карельский филиал АН СССР, Институт геологии, 1978.

26. Иксанов И .Я., Каплан С.А., Козлов A.C., Шпекторов A.JI., Хаустов М.Г. Сейсморазведка МОВ-ОСТ в интегрированном анализе и комплексной интерпретации геоинформации // Геоинформатика, 2006, № 3, с. 66-74.

27. Инструкция по гравиразведке. М.: Недра, 1980, 80с.

28. Картвелишвили K.M. Планетарная плотностная модель и нормальное гравитационное поле Земли. — М.: Наука, 1982.

29. Кашик A.C., Кириллов С.А., Ческис B.JI. Решение геологических задач в четырехмерном многопараметровом пространстве // Геофизика: Специальный выпуск «Технологии сейсморазведки-Ii», 2003, с. 3-16.

30. Кобрунов А.И. Математические основы теории интерпретации геофизических данных. М.: ЦентрЛитНефтеГаз, 2008, 288 стр.

31. Комплексирование методов разведочной геофизики: Справочник геофизика. / Под редакцией В.В. Бородового, A.A. Никитина. М.: Недра, 1984, 384с.

32. Кузнецов О.Л., Никитин A.A., Черемисина E.H. Геоинформационные системы. М.: ИЦ ВНИИгеосистем, 2005, 346с.

33. Кутрунов В.Н., Дмитриевский М.В. Аналог интерполяционного метода крайгинга без геостатистического обоснования. — Вестник ТюмГУ, №3, 2001, с. 208-216

34. Ларин Г.В. Эффективный компьютерный инструментарий геолога и геофизика при изучении нефтегазовых месторождений // Геофизика, 2010, № 3, с. 3-15.

35. Ломтадзе В.В. Программное и информационное обеспечение геофизических исследований. -М.: Недра, 1993, 268с.

36. Матусевич A.B. Объемное моделирование геологических объектов на ЭВМ. -М.: Недра, 1988, 184с.

37. Никитин A.A. Теоретические основы обработки геофизической информации. -М.: Недра, 1986. 342с.

38. Никитин A.A., Петров A.A. Основные процедуры обработки и интерпретации нестационарных геофизических полей. Геофизика №3, 2007, с. 63-70.

39. Никитин A.A., Петров A.B., Зиновкин C.B., Адаптивные процедуры выделения геологических объектов в геофизических полях. Из. Вузов. Геология и разведка, 2006, № 3, с. 50-56.

40. Никитин A.A., Хмелевский В.К. Комплексирование геофизических методов. Изд. ГЕРС, г. Тверь, 2004, 294с.

41. ПАНГЕЯ- 10 лет. Специальный выпуск. Геофизика, 2004, 106с.

42. Петров A.B., Никитин A.A., Зиновкин C.B. Развитие статистических приемов обработки и интерпретации геофизических полей в компьютерной технологии «КОСКАД 3D». Из. Вузов. Геология и разведка 2007, № 7 с.68-74.

43. Петров A.B., Трусов A.A. Компьютерная технология статистического и спектрально-корреляционного анализа трехмерной геоинформации КОСКАД 3D // Геофизика, 2000, № 4, с. 29-33.

44. Приезжев И.И. Информационные технологии интерпретации комплекса геолого-геофизических данных для геологического моделирования. Диссертация на соискание ученой степени доктора технических наук. Москва, РГГРУ, 2010, 212с.

45. Приезжев И.И. Построение распределений физических параметров среды по данным гравиразведки, магнитометрии и сейсморазведки. Геофизика №3, 2005, стр. 46.

46. Приезжев И.И. Уточнение геологической модели по данным гравитационного поля на основе критериальных методов решения обратных задач геофизики. Геофизика №1, 2010 стр. 65-68.

47. Романюк Т.В. Некоторые вопросы комплексной интерпретации данных ГСЗ и гравиметрии. Теория и практика интерпретации гравитационных аномалий. М. ИФЗ АН СССР, 1982, с. 26-59.

48. Савинский И.Д. О решении некорректной задачи при пересчете потенциального поля на нижние уровни. Изв. АН СССР, Физика Земли №6, 1967, с. 72-92.

49. Серкеров С.А. Гравиразведка и магниторазведка: Учеб. для вузов. — М.: ОАО Издательство "Недра", 1999, 437с.

50. Серкеров С.А. Корреляционные методы анализа в гравиразведке и магниторазведке. -М.: Недра, 1986, 247с.

51. Серкеров С.А. Спектральный анализ гравитационных и магнитных аномалий. -М.: ОАО Издательство "Недра", 2006.

52. Страхов В.Н. Использование методов теории функций комплексного переменного в плоских задачах гравиразведки и магниторазведки. Учебное пособие. М.: АН СССР, Институт Физики Земли им. О.Ю. Шмидта, 1982, 120с.

53. Страхов В.Н. Основные идеи и методы извлечения информации из данных гравитационных и магнитных наблюдений. — М.: ИФЗ АН СССР, 1979, с. 146-269.

54. Страхов В.Н. Основные направления теории и методологии интерпретации геофизических данных на рубеже XXI столетия. Ч. I // Геофизика №3, 1995, с. 9-18.

55. Хемминг Р.В. Численные методы. М.: Наука, 1968, 400с.

56. Черемисина Е.Н., Митракова О.В., Финкелылтейн M .Я. ГИС-ИНТЕГРО -инструмент постановки и решения природопользовательских задач. Информационный бюллетень № 3, 1988, ГИС-Ассоциация.

57. Черемисина Е.Н., Никитин А.А. Геоинформационные системы и технологии. Учебник для вузов. М.: ГНЦ РФ ВНИИгеосистем, 2011, 376с.

58. Baranov V. A new method for interpretation of aeromagnetic maps pseudo gravimetric anomalies. Geophysics, 1957, vol. 22 №2, p. 359-383.

59. Bott M.P., Smith R.A. The estimation of the limiting depth of graviting bodies. Geophysical processing, 1958, vol. 6, №1, p. 1-10.

60. Fairhead J.D., Odegard M.E. Advances in gravity survey resolution. // The Leading Edge, 2002, №1, p. 36-37.

61. Fedi M., Rapolla A. 3-D inversion of gravity and magnetic data with depth resolution. Geophysics, 1999, vol. 64, №2, p.452-460.

62. Hammer S. Deep gravity interpretation by stripping. Geophysics, 1963, vol. 28, №3, p. 369-378.

63. Hansen R.O. Gravity and magnetic methods at the turn of the millennium. Geophysics, 2001, vol. 66, №1, p. 36-37.

64. Moraes A.V.R., Hansen R.O. Constrained inversion of gravity fields for complex 3-D structures. Geophysics, 2001, vol. 66, №2, p. 501-510.

65. Silva J.B.C., Medeiros E.W., Barbosa V.C.F. Potential-field inversion: Choosing the appropriate technique to solve a geologic problem. Geophysics, 2001, vol. 66, №2, p. 511-520.

66. Yang Z., Wei Y. The gravity & seismic data joined formation separation technique for deep structure study. // SEG Technical Program Expanded Abstracts, 2005, p. 635-638.

Информация о работе

Зиновкин, Сергей Владимирович
кандидата технических наук
Москва, 2011
ВАК 25.00.10

Диссертация

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

Автореферат

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