WWW.KNIGA.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА - Книги, пособия, учебники, издания, публикации

 

Pages:   || 2 |

«Бисеркин Игорь Алексеевич Технология интерпретации данных потенциальных полей при изучении строения земной коры (на примере Байкитской антеклизы) Специальность 25.00.10 ...»

-- [ Страница 1 ] --

РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГЕОЛОГОРАЗВЕДОЧНЫЙ

УНИВЕРСИТЕТ им. Серго Орджоникидзе

Бисеркин Игорь Алексеевич

Технология интерпретации данных потенциальных

полей при изучении строения земной коры (на примере

Байкитской антеклизы)

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

ископаемых

ДИССЕРТАЦИЯ

на соискание ученой степени кандидата технических наук Научный руководитель кандидат технических наук, профессор А.М. Лобанов Москва, 2013г.

ОГЛАВЛЕНИЕ

ВВЕДЕНИЕ …………………………………………………………………………………… ГЛАВА 1. Современное состояние технологий построения 3D плотностных моделей……………………………………………………………………………………….. 1.1. Современное состояние методов решения обратных задач грави- и магниторазведки……………………………………………………………………………… 1.2. Общая схема решения обратной задачи. …………………………………..… 1.3. Компьютерные технологии по решению прямых и обратных задач……….. 1.4. Основные тенденции при моделировании потенциальных полей………….. ГЛАВА 2. Метод, методика и технология моделирования 3D плотностных моделей……………………………………………………………………………………… 2.1. Метод построения 3D плотностной модели………………………………….. 2.2. Процесс подбора 3D сеточной плотностной модели слоисто-блоковой среды…………………………………………………………................ 2.3. Технологические средства ГИС ИНТЕГРО………………………………….. ГЛАВА 3. Детализация плотностной модели земной коры на основе многопризнаковой фильтрации в среде ГИС-ИНТЕГРО…………………………… 3.1. Алгоритм многопризнаковой фильтрации………………………………….... 3.2. Методика применения алгоритма фильтрации для построения плотностной модели земной коры…………………………………………………………………………………………... ГЛАВА 4. Применение технологии построения плотностной слоисто-блоковой 3D модели на примере Байкитской антеклизы …………………………………………... 4.1. Построение плотностной модели с помощью алгоритма многопризнаковой фильтрации………………………………………………………………………………....... 4.2.Детализация плотностной модели с помощью алгоритма многопризнаковой фильтрации…………………………………………………………………………………... Заключение …………………………………………………………………………………... Список литературы………………………………………………………………………….. Список рисунков…………………………………………………………………………....... Введение Актуальность. Проблема изучения глубинного строения земной коры и верхней мантии является одной из стратегических направлений геофизических исследований, обеспечивающих развитие наук о Земле. При этом гравиразведка является одним из основных методов изучения строения земной коры.




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

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

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

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

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

Задачи исследований.

1. Анализ современного состояния методик и компьютерных технологий 3D моделирования среды.

2. Создание технологии интерактивного подбора плотностных 3D моделей для блоково-слоистой геосреды на основе использования ГИС ИНТЕГРО.

3. Разработка детализации плотностной модели земной коры на основе многопризнаковой фильтрации.

4. Апробация созданной технологии на территории Восточной Сибири.

Практическая значимость.

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

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





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

Научная новизна исследований.

1. Предложен метод формирования априорной 3D сеточной модели слоисто-блоковой среды.

2. Предложен способ изменения параметров сеточных моделей для реализации интерактивного 3D подбора.

3. Разработана технология, обеспечивающая подбор модели с использованием сеточного распределения плотности и векторно-заданных объектов.

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

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

1. Предложенный метод формирования априорной 3D сеточной модели слоистоблоковой среды на основе использования комплекта структурных карт контактных поверхностей и распределения плотностей в каждом слое, обеспечивает эффективное решение прямой задачи гравиразведки.

2. Технология построения (интерактивного подбора) плотностной 3D модели земной коры в среде ГИС ИНТЕГРО, включающая метод формирования 3D сеточных моделей слоисто-блоковой среды, а также путем сочетания сеточного и векторного представления геолого-геофизических данных обеспечивает оптимальное решение обратной задачи гравиразведки, адекватное реальным средам.

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

Апробация и публикации.

Основные положения и результаты, представленные в диссертационной работе, докладывались на международной научной конференции «Новые идеи в науках о Земле»(Москва 2011, 2012), 18-ой Международной конференции «Геологическая среда, минерагенические и сейсмотектонические процессы»(Воронеж 2012), 13-ой международной научно-практической конференции по проблемам комплексной интерпретации геологогеофизических данных при геологическом моделировании месторождений углеводородов «Геомодель». (Геленджик 2011).

Диссертация основана на теоретических, методических и экспериментальных исследованиях, выполненных автором в 2009-2013гг. Основные методические и технологические результаты получены непосредственно автором. По результатам выполненных исследований опубликовано 4 печатные работы, из них 2 статьи в журналах входящие в перечень ВАК и были представлены в материалах 4 международных конференциях.

Объем и структура работы.

Диссертация состоит из введения, четырех глав, заключения, списка литературы, содержит 86 страниц машинописного текста.

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

Глава 1. Современное состояние технологий построения 3D плотностных моделей 1.1.Современное состояние методов решения обратных задач грави- и магниторазведки Проблема обработки и интерпретации геофизических данных при изучении глубинного строения земной коры связана с необходимостью поиска новых подходов к геокартированию, выделению крупных нефтегазоносных провинций и рудных поясов.

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

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

Традиционно принято выделять четыре крупных класса задач, определивших эволюцию методов интерпретации гравиметрических и магниторазведочных данных:

1. Обнаружение объектов.

2. Пространственная локализация.

3. Расчленение геологического разреза.

4. Детальное описание (прогнозирование геологического разреза).

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

Учитывая это необходимо согласиться и принять к руководству вывод В.Н. Страхова о неадекватности современного состояния теории и методов гравиметрии реальной геофизической практике [34,44] и необходимости создания адекватных этой геофизической практике концепций и следующих из нее конструктивных решений. При этом основным резервом повышения эффективности геофизического интерпретационного обеспечения ГРР, является создание теории и методов, основанных на вовлечение в активный процесс извлечения информации геологических и тектонофизических концепций, моделей происхождения и эволюции изучаемых объектов, а в некоторых случаях и подчинение самих этих методов и порождаемых ими интерпретационных технологий этим концепциям. Именно в этом должен состоять диалект методов извлечения информации из геофизических данных [17,26,28]. Именно это направление объективно и неизбежно будет основным в настоящем и ближайшей перспективе, поскольку отражает реальные потребности практики ГРР и направлено на обеспечение ее потребностей.

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

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

Применяемые подходы, методы и средства для решения обратных задач отличаются большим разнообразием.

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

В рамках детерминированного подхода приемом приближенного решения обычно некорректно поставленных задач является метод регуляризации А. Н. Тихонова, для которого широкое распространение получили способы, основанные на вариационном принципе, т.е. на использовании сглаживающего параметрического функционала Тихонова [46,47]. Этот метод требует решения систем линейных или нелинейных уравнений достаточно большого размера.

Так же в рамках данного подхода, относятся методы аналитического продолжения, а также условно корректных трансформаций, направленные на определение положения особых точек поля. Методы аналитического продолжения получили практическое использование в виде модифицированного метода Б.А. Андреева в компьютерной системе КОСКАД-3Д и в виде спектрального разложения поля в алгоритмах И.И. Приезжев, данные алгоритмы так же были реализованы в ГИС-Интегро Геофизика. Среди трансформаций по определению положения особых точек наиболее известны методы А.В. Цирульского, В.М. Березкина[52].

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

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

В рамках вероятностно-статистического подхода решение обратных задач потенциальных полей базируется на теории статистических гипотез, разработанной Ф.М. Гольцманом и Т.Б. Калининой. Задача нахождения вектора параметров геологического разреза сводится к вычислению апостериорной вероятности ( / g) вектора при условии вектора наблюденных значений поля, например, g.

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

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

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

Смешанный вероятностно-детерминистский подход к решению линейных обратных задач грави- и магниторазведки предложен А.С. Долгалем и П.И. Балком. При этом подходе накладываются условия на адаптивную помеху, принимая ее медианные значениями нулю, а для параметров модели задается равномерное распределение в заданных пределах их изменения. Решения вычислительных экспериментов по методу минимизации эмпирического риска с решением по минимуму среднеквадратической невязки наблюденного и модельного полей (для детерминированного подхода) показали повышение точности в 2,5 раза оценки параметров аномалиеобразующих объектов.

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

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

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

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

1.2. Общая схема решения обратной задачи В общем случае схема решения обратной задачи гравиразведки включает:

1) выбор модели для решения прямой задачи, которая должна:

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

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

в) допускать свободные изменение параметров в процессе подбора;

г) позволять получение вертикальных или горизонтальных сечений, построение погоризонтных планов.

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

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

Модель первого приближения считается удачной, если гравитационный эффект от нее грубо отражает все основные особенности аномального поля;

3) уточнение параметров изучаемой геосреды. Изменения в значениях параметров осуществляется в тех частях модели, где данные отсутствуют, а физические и геометрические параметры изменяются в заданных, разумных пределах.

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

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

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

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

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

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

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

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

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

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

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

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

На данный момент выделяются следующие подходы к решению поставленных задач. В рамках вероятностно-статистического подхода выделяются технологии:

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

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

Так же в данном программном комплексе реализован алгоритм оценки параметров аномалиеобразующих объектов на основе модификаций метода Андреева.

Метод построен на использовании алгоритма фильтрации наблюдений в окне живой формы. В основе лежит предложенная Б.А. Андреевым[3] идея об оценке частотных составляющих поля с помощью фильтрации данных в окнах увеличивающегося размера.

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

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

-ПАНГЕЯ (ЗАО Пангея), имеющая много общего с технологией КОСКАД-3Д (ПАНГЕЯ – 10 лет. Специальный выпуск НТЖ ГЕОФИЗИКА, 2004, 106 с.) ;

Разработанная и опробованная в ряде регионов страны автоматизированная система комплексной обработки и интерпретации геолого-геофизической информации «ПАНГЕЯ», до настоящего времени была адаптирована на решение, в основном, задач нефтяной геологии.

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

Система обеспечивает на вероятностном уровне решение широкого круга геологических задач, в частности задач, направленных на:

1) структурно-тектоническое районирование площадей по комплексу геофизических признаков методами таксономической классификации;

2) классификацию ловушек углеводородов по степени нефтегазоперспективности с ранжированием их по очередности опоискования на основе кластерного анализа с обучением;

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

4) оконтуривание залежей нефти и газа с определением положения ВНК и ГНК на основе количественного прогнозирования эффективных мощностей;

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

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

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

В обобщенном виде схема обработки геолого-геофизических данных по системе «ПАНГЕЯ» при прогнозировании месторождений углеводородов включает в себя следующие основные этапы:

1) постановка задачи и создание физико-геологических моделей объектов поиска и разведки, выбор эталонов различных геофизических кластеров;

2) формирование многомерного признакового пространства, расчет дополнительных параметров физических полей, анализ законов распределений и изменчивости признаков с оценкой их статистических связей;

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

4) собственно прогнозирование на ЭВМ по различным программно-алгоритмическим комплексам на основе различных математических моделей с использованием кластерного, факторного и регрессионного анализов и теории информации методами классификации и распознавания образов;

5) интерпретация и геологическое истолкование результатов прогнозирования.

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

Комплексная интерпретация геолого-геофизических данных и прогнозирование на основе системы «ПАНГЕЯ» могут быть реализованы на всех стадиях ГРР.

Успешное применение системы «ПАНГЕЯ» предусматривает соблюдение определенной последовательности операций вычислительного процесса, что обеспечивает получение устойчивых решений. Подобный результат достигается с помощью рационального графа обработки многопараметровой информации как применительно к различным этапам исследований, так и при решении отдельных задач прогноза.

- Mult-Alt, созданная Д.Ф. Калининым в ОАО «Геотехника», и используемая для потенциальных полей при решении задач рудной геофизики на основе байесовского подхода Ф.М. Гольцмана, Т.Б. Калининой[21]; Технология «MultAlt» использует оптимальные статистические критерии и алгоритмы принятия решений по комплексу количественных и качественных признаков (геофизических, геохимических, геологических и др.). В процессе формализованных построений определяются численные оценки достоверности и качества решений об искомых альтернативных объектах на любых этапах комплексирования.

Компьютерная технология «MultAlt» имеет уникальную возможность контроля надежности и информативности решений с выявлением наиболее рационального комплекса признаков до проведения автоматизированной комплексной интерпретации. На различных стадиях работ технология предусматривает учет эвристических мнений экспертов-геологов - ГРАВИПАК и КОМПАК, два пакета, созданные В.В. Ломтадзе в ОАО «Иркутскгеология», в последнем реализуются методы факторного анализа и метод главных компонент по атрибутам потенциальных полей[32].

В рамках детерминированного подхода наиболее известными технологиями являются:

- СИГМА-3Д, созданная под руководством Ю.И. Блоха[4] и используемая в ОАО «Аэрогеофизика» при построении согласованных плотностных и магнитных моделей кристаллического фундамента по территории Восточно-Европейской платформы. Основными комплексами в данной технологии являются программа REIST и DVOP.

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

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

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

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

Число точек магнитной съемки, как правило, значительно больше числа элементарных призм.

Поэтому часть из них относится программным комплексом к активным, по которым ведется коррекция модели в процессе последовательных приближений. Другая часть точек наблюдений считается пассивной, для них вычисляется поле подбираемой модели, по ним же оцениваются погрешности подбора при каждой итерации. Применительно к данным гравиразведки описанный итерационный процесс реализует решение линейной обратной задачи, в то время как для магнитных аномалий осуществляется решение нелинейной обратной задачи. Это обстоятельство отличает технологию СИГМА-3D от других аналогичных компьютерных систем, для которых применяются линейное, чаще всего гармоническое приближение для аномалий.

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

Подобный процесс интерпретации полностью устраняет влияние на результаты постоянного регионального фона.

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

В технологию СИГМА-3D включен пакет DVOP, данная программа разработана для вычисления объемного распределения эффективной плотности и эффективной намагниченности. С помощью технологии DVOP можно выполнять изучение характера пространственной изменчивости магнитоактивных (гравитирующих) элементов строения нижнего полупространства, что позволяет судить об их морфологии, элементах залегания и взаимоотношениях. В программе реализуется одно из наиболее перспективных направлений в изучении потенциальных полей Земли - интерпретационная томография, то есть система исследования геологических структур по гравитационному либо магнитному полю, позволяющая получать их послойное отображение. Теоретические основы интерпретационной томографии к настоящему времени проработаны крайне слабо. Очевидно, в общем случае невозможно однозначно разделить вклады в аномальные гравитационные и магнитные поля от разных слоев разреза. Однако это ни в коей мере не является препятствием для разработки методов решения этой задачи в частных, но практически важных случаях.

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

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

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

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

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

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

Несомненным достоинством технологии СИГМА-3D является то, что исходные гравитационное и магнитное поля задаются в реальных точках наблюдений, т.е. при неравномерной сети с учетом высоты над поверхностью геоида. Кроме этого представляет интерес корректность использования в процессе моделирования сложных по своей природе магнитометрических наблюдений.

Компьютерная технологи комплексной интерпретации данных сейсморазведки и гравиразведки GCIS. Технология количественной комплексной интерпретации данных сейсморазведки и гравиразведки GCIS-авторы Кобрунов, Петровский[27]. Предназначена для решения прогнозно-поисковых задач на нефть и газ. Технология основана на совместном решении обратных задач сейсморазведки и гравиразведки. При этом изучаемые объект задается в виде модели слоистой среды, с априори заданными скоростными и плотностными (атрибутами). Интерпретация сейсмических данных осуществляется в рамках кинематической задачи. Несмотря на то, что скоростные и плотностные границы не могут быть тождественными, но достаточно близки между собой, поскольку скоростная слоистая и плотностная слоистая модели являются образами одного и того же геологического объекта.

Построение согласованной сейсмоплотной (согласование по скорости и плотности) модели среды осуществляется путем минимизации обобщенного функционала (,) на основе интерактивного подбора. При этом предполагается, что минимум такого функционала (,) реализуется при условии соответствия с одной стороны, скоростной модели волновому полю отражений ОГТ, а с другой стороны, плотностной модели по гравитационному полю наблюдений. При каждом фиксированном векторе параметров среды индекс означает, что «сейсмические границы зафиксированы и вошли в конструкцию функционала. Аналогично при каждом фиксированном функционал (,) представляет собой критерий оптимальности на множестве эквивалентности для плотностных границ. Итерационный процесс состоит в последовательном нахождении решений, при котором результат, полученный при интерпретации данных гравиразведки, передается для конструирования функционала по данным сейсморазведки и последующего нахождения решения, далее найденное решение для сейсморазведки передается для конструирования функционала по данным гравиразведки и так далее. Описанный итерационный процесс сходится. Его результатом является пара слоистых моделей- скоростная и плотностная, максимально близких между собой насколько это допустимо согласованностью волнового и гравитационного полей. Важным элементом приведенной технологии является параметризация модели среды, определяющая допустимые виды моделей в данной конкретной геологической ситуации, выражаемые через варьируемые в процессе решения обратной задачи параметры среды. Параметром служит одна единственная функция-функция Лагранжа, через которую и критерий оптимальности выражаются все границы среды. В рамках параметризации используется аппроксимация среды наборами призм, что обеспечивает расчет прямой задачи. Технология GCIS развивается путем математического моделирования динамики формирования изучаемого объекта, результаты которого является последовательность статистических моделей, а параметризация заключается в оценке параметров по управлению динамикой процесса формирования углеводородов.

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

Технология GMSYS на основе канадского программного комплекса Geosoft.

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

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

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

Доступна возможность построения трехмерной сети гравитационного или магнитного поля, рассчитанной для 2 3 4 D модели для изучения аномалий, находящихся вне линии профиля.

Вычисление эффекта может быть произведено от отдельных блоков модели.

GMSYS-3D представляет собой программный пакет трехмерного моделирования гравитационных и магнитных аномалий, предназначенный для модели слоистой среды.

Модель формируется из имеющегося набора сеточных моделей структурных горизонтов (грид-файлов, совместных со средой Geosoft), разделяющих слои с заданными физическими параметрами (плотности и магнитной восприимчивости). В процессе решения обратной структурной задачи(structural inversion) программа может деформировать и двигать вверх или вниз структурный горизонт(не изменяя при этом значения плотности вышележащего и нижележащего слоев) для наилучшего подбора формы поверхности для данной модели. Для получения более полного представления о геологическом строении производятся интерпретация аномалий полного тензора градиента гравитационного поля (Full Tesor Gravity Gradient Joint Inversion). Представлена возможность выполнять интерпретацию с использованием любой компоненты тензора градиента гравитационного поля совместно с магнитными данными и полным градиентом гравитационного поля, указав весовые коэффициенты для каждого из восьми возможных ограничивающих параметров. Программа отображает среднеквадратическое расхождение на каждом этапе подбора. Таким образом, пользователь имеет возможность наблюдать за результатами инверсии и в любой момент остановить процесс подбора. Поверхности слоев удобны для редактирования. Свойства слоев могут задаваться постоянными значениями или изменяться по латерали и с глубиной.

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

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

Корректная вставка двумерных разрезов в трехмерную модель.

Для сравнения подобранного поля с отфильтрованными значениями наблюденного поля в программе доступна возможность применить фильтры к кривым подобранного поля.

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

Доступ к Voxi Еarth Modelling осуществляется из Oasis montaj модели с использованием мощности облачных вычислений для быстрого создания трехмерных моделей крупных массивов разведочных данных.

Технология моделирования геологической среды система VECTOR.

Интерпретационная система «ВЕКТОР», разработанная в 80 – 90 годы XX века коллективом лаборатории геопотенциальных полей Горного института УрО РАН (научный руководитель – В.М.Новоселицкий).

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

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

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

Предложены пути комплексирования сейсмических исследований 3D и современной гравиметрии для повышения информативности геофизических работ при детальном изучении месторождений углеводородов. Результаты интерпретации гравиметрических данных позволяют выявлять плотностные (и скоростные) неоднородности среды как в верхней части разреза, так и в более глубоких отложениях, и, в итоге, построить детальную геологическую модель месторождения. Решение прямой задачи гравиразведки в рамках принципа контактных поверхностей, программный модуль позволяет вычислять эффект от набора поверхностей, образующих модель среды, и подбирать (корректировать) геометрию и плотностные аномалии границ. Так же в программе реализовано Решение трехмерной нелинейной обратной задачи гравиметрии монтажным методом. Монтажный подход к решению обратной задачи снимает проблему неустойчивости в ее классическом понимании, так как конечная размерность модели и природные ограничения на ее плотностные и геометрические параметры изначально приводят к компактному множеству возможных решений, а учет определенного объема априорной информации обеспечивает высокую геологическую содержательность конечного результата. Степень достоверности результатов количественной интерпретации гравитационных аномалий сводится к построению монтажным методом системы вложенных множеств допустимых решений обратной задачи (отвечающих различным уровням неопределенности в априорных данных). Сопоставление объемов информации, содержащейся в этих множествах, позволяет получить сравнительную оценку надежности обнаружения возмущающих масс в тех или иных фрагментах изучаемого объема геологической среды. Предложенный метод может существенно повысить вероятность вскрытия искомых аномалиеобразующих объектов в заданных интервалах глубин поисковыми и разведочными скважинами, рекомендованными по гравиметрическим данным. Вероятностно-детерминистский подход к решению линейных обратных задач в гравиметрии, магнитометрии и электрометрии. Предложен принципиально новый метод решения линейных обратных задач геофизики, базирующийся на синтезе функциональноаналитического и вероятностно-статистического подходов к интерпретации. Отличительными особенностями метода являются возможность подавления интенсивных знакопеременных помех, закон распределения которых может существенно отличаться от нормального, а также приближенная оценка точности определения физических параметров геологических объектов.

Технология моделирование геологической среды 3D GeoModeller 3D GeoModeller представляет собой программный пакет, обеспечивающий трехмерное моделирование геологической среды и изначально ключающий в себя логику геологического моделирования на основе первичных полевых наблюдений и общих представлений о геологическом строении района работ. Особенность 3D GeoModeller заключается в том, что он изначально разрабатывался в рамках проекта «Виртуальная Земля» («Terre Virtuelle») для решения геологических задач, поэтому в самом программном обеспечении заложены основы геологической логики(взаимоотношения разновозрастных объектов, разрывных нарушений, основные принципы интерполяции данных, оригинальный метод решения прямой задачи гравиразведки, и др.), а так же возможность построения объемных объектов. Кроме собственно геологического моделирования 3D GeoModeller обеспечивает решения в пространстве прямой и обратной задач гравии- и магниторазведки. Технология 3D GeoModeller отличает от рассмотренных выше изначальная ориентация на логику геологического моделирования, на основе первичных полевых наблюдений и общих представлений о геологическом строении района работ. Использование геофизических наблюдений в процессе моделирования подчинено построениям, базирующимся на геологической информации.

- АДМ-3Д, реализованная В.А. Кочневым, И.В. Гозом для трехуровневых магнитных съемок; технология решения прямых и обратных задач 3D гравиметрии и магнитометрии.[29,30] Основой представляемых технологий являются пакеты программ ADG-3D и ADM-3D и рекомендации по их применению.

В основу разработки положена слоисто-блочная модель среды, состоящая из набора прямоугольных параллелепипедов с заданными размерами по осям X и Y. Центры параллелепипедов образуют равномерную сетку, в узлах которой находятся аномальные значения гравитационного или магнитного поля. Высота параллелепипедов равна толщине слоя в узле сетки. Число таких параллелепипедов равно N=nx*ny*nz, где nx и ny – число узлов по осям X и Y, а nz – число слоев. Каждый блок имеет свою плотность или избыточную плотность.

При известной избыточной плотности в блоках прямая задача для любой точки пространства решается однозначно. Определенные трудности возникают при решении обратной задачи, в которой по значениям избыточного ускорения свободного падения g ij определяются ijk. Как правило, число неизвестных значительно больше числа уравнений, следовательно, задача имеет бесконечное множество решений. Для преодоления этой и других трудностей разработан адаптивный метод, позволяющий найти решение, ближайшее к начальному.

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

Для решения 3D задач магнитометрии использована та же слоисто-блочная модель, но вместо плотностей каждый блок имеет свою магнитную восприимчивость.

Внешним магнитом, намагничивающим блоки, является магнитное поле Земли с компонентами Hx, Hy, Hz, которые являются входными параметрами. Аналитические формулы для параллелепипедов более сложные, чем в гравиметрической задаче. Здесь приходится рассчитывать три компоненты аномального поля и получать из них вектор или скаляр аномалии полного магнитного поля Ta. При известных параметрах модели прямая задача расчета всех компонент Hx, Hy, Hz и Ta решается однозначно. Трудности здесь возникают при исключении аномалий на краях модели. Краевой эффект был преодолен путем продолжения в бесконечность блоков, окаймляющих модель. Основные трудности возникают при решении обратной задачи. Они преодолеваются тем же путем, что и в гравиметрии.

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

- Программный комплекс института геофизики УрО РАН, созданный П.С. Мартышко и И.Л. Пруткиным; В данном комплексе были реализованы новые подходы к 2D и 3D моделированию среды по данным магнитного поля, предложенные в последние годы.

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

- продолжение поля в верхнее полупространство на уровень Н для устранения поверхностных источников;

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

- пересчитанное поле вверх продолжается в нижнее полупространство на глубину Н.

При этом формула Пуассона (1):

используется в качестве интегрального уравнения в котором функция Ux, y, H задана, а в его ядре величина Н заменяется на значение 2Н. Искомым решением становится функция Ux, y,-H.

Решение уравнения (1) представляет собой некорректную задачу, что требует:

- регуляризации, которая осуществляется по схеме М.М. Лаврентьева путем решения уравнения K E W U.

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

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

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

- Программный комплекс С.А. Тихоцкого, реализованный в отделе гравиметрии ИФЗ РАН[48]; С.А. Тихоцким с соавторами предложен алгоритм, обеспечивающий восстановление гармонического аномального поля по данным об аномалиях T. Суть алгоритма реализуется в следующей схеме. За U обозначается скалярный потенциал, соответствующий аномальному полю T. За нулевое приближение гармонической компоненты T0 U t, где t, принимается измеренное поле T. Далее используется вспомогательное распределение магнитных масс, которое порождает скалярный магнитный потенциал V 0, соответствующей на исходной совокупности точек наблюдений гармонической компоненте и совпадающей с T:

Негармоническая часть поля T потенциала V 0 будет равна:

В качестве следующего (первого) приближения гармонической компоненты можно принять:

Это приближение T01 вновь аппроксимируется полем вспомогательного распределения источников, т.е. находится распределение магнитных масс, порождающее скалярный магнитный потенциал V 1 такой, что на исходной совокупности точек наблюдений V 1t T01. После повторения операций (2-3) получается следующее (второе) приближение гармонической компоненты.

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

где U – потенциал аномального поля T 0, а V i - потенциалы, порождаемые вспомогательным распределением масс и определяемые на каждом шаге итерационного процесса из условия:

Доказано С.А. Тихоцким утверждение: если итерационный процесс, определяемый формулами (4-5) сходится, то последовательность T0 k сходится к T0, т.е. оценка гармонической компоненты является состоятельной. Моделированием показано, что алгоритм сходится в диапазоне значений магнитного поля вплоть до 15 000 нТл. Независимо от соотношения направлений нормального поля и порождающей аномальное поле намагниченности.

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

Так, на основе алгоритма (1-5) решается задача выделения низкоамплитудного аномального поля, порождаемого намагниченными горизонтами в осадочном чехле платформ (например, траппами) на фоне высокоамплитудного поля пород кристаллического фундамента, определения положения особых точек магнитовозмущающего объекта по данным аномалий T. При этом требуется предварительно восстановить гармоническое аномального поля по рассмотренному алгоритму.

Технология решения прямых и обратных задач гравиразведки с применением аппарат спектральных преобразований ГИС-Интегро ГИС ИНТЕГРО имеет 3D редактор, позволяющий визуализировать и редактировать данные в пространстве.

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

В общем случае, выражение для спектра физического параметра на глубине z выглядит следующим образом:

В случае потенциала поля силы тяжести формула принимает вид:

Это означает, что, используя в качестве Z (1, 2, z ) произвольную функцию, зависящую от глубины, можно получить такое распределение физического параметра G(1, 2, z ), которое будет строго удовлетворять измеренному полю (1, 2,0). Существенным и единственным теоретическим ограничением является существование обратного преобразования Фурье для G (1, 2, z ). Разработанная в ГИС-Интегро методика подразумевает использование данной формулы для определения распределения физического параметра на основе априорно заданной плотностной модели. Соответственно функцию Z (1, 2, z ) можно задавать в виде куба или разреза априорно заданных атрибутов, или других геофизических атрибутов, имеющих связь с плотностью или намагниченностью.

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

Реализация алгоритма позволяет уверенно обнаруживать аномалии в поле силы тяжести, и обеспечивает переход к локализации тел, порождающих данные аномалии. При этом, уровень помех на порядок ниже, по сравнению с технологией решения обратной задачи, использующей аппарат Приезжева в программном комплексе «Коскад». Опробование на модельных данных показало, что амплитуда помех составляет 3-5% от амплитуды полезного сигнала.

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

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

Этот метод основан на том, что формула колебания для точечной массы I достаточно хорошо приближается суммой трх гауссианов:

s1=0. s2=1. s3=4. a1=0. a2=0. a3=0. проекции центра точечной массы на поверхность. Таким образом, можно получить приближнное аналитическое решение для оценки положения, массы и глубины точечной массы для вейвлет функции g ( x).Для g ( x) же было получено численное решение.

Для создания 3-х мерной эффективной плотности ( x0, h) измеренной гравиметрикой f используется пространство вейвлет разложения W (, s ) f, где x0 2 – смещение по горизонтали, а s Ch – смещение по глубине, C 1.38846642104819 :

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

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

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

Основные тенденции при моделировании потенциальных полей определяются:

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

- использованием генетических алгоритмов при решении обратных задач (И.И. Приезжев);

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

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

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

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

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

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

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

Построение согласованных по плотности и магнитной восприимчивости 3D моделей геосреды достаточно оправдано при изучении интрузивных массивов, магматических и метаморфических комплексов пород. Построение таких объектов эффективно при картировании пород кристаллического фундамента. С этой целью используется программный комплекс СИГМА-3D (Ю.И. Блох), с помощью которого реализовано построение обновленных карт фундамента на территории Русской (Восточно-Европейской) платформы.

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

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

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

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

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

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

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

4. Имеется хорошее аналитическое обеспечение решения прямой и обратной 3д задачи как для тел, так и сеточных моделей, но не реализован интерактивной 3D подбор.

5. В тоже время современные компьютерные технологии визуализации и редактирования пространственных объектов позволяют работать с довольно сложной средой.

В одних системах реализовано хорошее аналитическое обеспечение решения задач гравики, в других есть возможности визуализации, но основными нужными функциями и программами обладает только ГИС ИНТЕГРО, поэтому технология будет строиться на базе ГИС ИНТЕГРО.

Таким образом, существуют предпосылки для реализации решения обратной задачи гравиразведки методом интерактивного 3D подбора.

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

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

Технология реализована в среде ГИС ИНТЕГРО, в которой имеется аналитическое обеспечение для реализации технологии и широкие возможности визуализации и оперирования 3D данными.

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

2.1. Метод построения 3D плотностной модели Метод построения 3D плотностной модели реализуется для следующих моделей среды (рис.1):

-слоисто-блоковой среды, -слоисто-блоковой градиентной среды.

Исходными данными для построения модели в зависимости от типа среды являются (рис.2):

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

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

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

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

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

Слоистая среда, задаваемая для моделирования может иметь (рис.4):

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

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

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

Рассмотрим основные этапы проведения подбора.

2.2. Процесс подбора 3D сеточной плотностной модели слоисто-блоковой среды Процесс подбора 3D сеточной плотностной модели слоисто-блоковой среды реализуется по стандартной схеме и включает в себя 3 этапа (см. рис.6):

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

Рис.2. Исходные данные для построения 3D плотностной сеточной модели:

а) карты глубины залегания поверхностей в земной коре, б) карты скорости Vp на кровле глубинных поверхностей, в) пространственное векторное представление исходных данных (цветом на поверхностях показано распределение скорости) Рис.3. Технология построения 3D сеточных плотностных моделей разных типов:

А- для слоистой среды (плотность в слоях постоянная), Рис.3. Технология построения 3D сеточных плотностных моделей разных типов:

Б - для блоково-слоистой среды (плотность в слое изменяется по латерали), Рис.3. Технология построения 3D сеточных плотностных моделей разных типов:

В - для слоистой среды с градиентным изменением плотности Рис.4. Особенности исходных данных для технологии построения 3D плотностной сеточной модели:каждый нижний слой не пересекает верхний (А), слои могут быть субпараллельны (Б),слои могут выклиниваться или иметь прерывание мощности в небольшой области (В) Рис.5. Сеточная 3D плотностная блоково-слоистая модель Этап 1. Подготовка исходных данных и построение начальной модели. На этом этапе осуществляется подготовка исходных данных для построения 3D плотностной модели.

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

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

Источником данных о глубине и плотности слоев могут быть скважинные данные.

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

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

Исходные данные могут быть получены из результатов исследований площади другими геофизическими методами. Например, сейсморазведкой. При изучении глубинного строения земной коры используются в первую очередь данные ГСЗ о структуре коры и распределении скорости в ней, положении границы кора-мантия, а также корреляционная зависимость =f(Vp), позволяющая в первом приближении оценить плотностной разрез. В данном случае исходные данные о глубине и скорости слоев будут привязаны к профилям.

Для построения начальной модели исходными данными могут служить результаты решения 2D обратной задачи гравиразведки или магниторазведки. Например, полученные в результате определения параметров отдельных аномалообразующих тел. В этом случае данные о глубине залегания верхних кромок тел и плотности также будут распределены по неравномерной сети точек. Кроме того, исходными данными могут быть результаты подбора плотности по серии 2D разрезов.

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

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

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

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

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

Этап 2. Проведение подбора.

На этом этапе подбор осуществляется в следующей последовательности:

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

Наблюденное гравитационное поле Рис.6 Схема основных этапов проведения 3D моделирования Затем для построенной начальной 3D сеточной модели проводится решение прямой задачи гравиразведки методом Приезжева И И., в ходе которой рассчитывается двумерное гравитационное поле по регулярной 2D сети.

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

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

- редактирование исходных карт рельефа поверхностей, карт распределения плотности (рис.7), - проводя оценку параметров модели на 2D разрезах, - путем добавления в модель локальных неоднородностей.

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

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

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

3 Этап. Интерпретация полученной модели.

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

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

В ее аналитическое обеспечение в том числе включены: программы решения прямой задачи для 3D сеточных моделей (алгоритм Приезжева И.И.) и для отдельных тел (призм, многогранников, пластин, сфер, уступов и др.), а также для решения обратной задачи как в 2D варианте, так и с получением 3D распределения плотности (статистическое и корреляционное зондирование, вейвлет-анализ, метод Андреева, Приезжева И.И. и др.). Кроме того, в технологии используются программы интерполяции, разработанный специальный модуль для формирования 3D сети по набору двумерных поверхностей, программа пересчета объемных векторных объектов в сеточную форму и др.

В ГИС ИНТЕГРО имеются технологические средства обеспечивающие визуализацию и анализ 3D моделей, совместный анализ трехмерных данных с 2D моделями, а также их сопоставление с геологическими картами, скважинными данными, данными о глубинных поверхностях и т.д.

2.3.Технологические средства ГИС ИНТЕГРО Технологические средства ГИС ИНТЕГРО осуществляют:

- представление 3D сеточных моделей в виде кубов, горизонтальных и вертикальных сечений - визуализацию 3D моделей в виде отдельных тел, поверхностей - геометризацию моделей в вертикальных, горизонтальных плоскостях и в объеме - сопоставление 3D моделей полученных разными методами - сопоставление 3D моделей (представленных сеточными моделями, поверхностями, телами) и 2D моделей - отображение физических свойств на поверхности слоев. Данные, представляющие собой пространственную информацию по поверхностям и разрезам (сечениям) могут быть изображены в 3D визуализаторе как композиция физических свойств и геометрии. Например, подобные композиции физических свойств могут отображаться не плоско по сечениям, а укладываться непосредственно на 3D поверхность соответствующих слоев.

- сечение 3D модели заданной поверхностью Визуализатор 3D позволяет представление моделей залитыми поверхностями, точками, каркасом, объемом, сечениями, осуществлять управление шкалой раскраски, изменение масштабов по x, y, z, управление прозрачностью 3D объектов.

Вывод: Предложенная технология формирования 3D блоково-слоистых моделей осуществляемая с использованием ГИС ИНТЕГРО позволяет проводить моделирование плотностного разреза с использованием 3D сеточных моделей с привлечением векторных объектов, что позволяет с наименьшими затратами труда задавать и корректировать параметры модели.

Глава 3. Детализация плотностной модели земной коры на основе многопризнаковой фильтрации в среде ГИС-ИНТЕГРО.



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

«ВЕСТНИК НОВЫХ МЕДИЦИНСКИХ ТЕХНОЛОГИЙ – 2010 – N 1 Электронное издание УДК 612.017+615.399 АНТИОКСИДАТНОЕ И ИММУНОПРОТЕКТОРНОЕ ДЕЙСТВИЕ ЭКСТРАКТА ЛИЧИНОК ВОСКОВОЙ МОЛИ ПРИ ОКИСЛИТЕЛЬНОМ СТРЕССЕ У КРЫС, ВЫЗВАННОМ ПОТРЕБЛЕНИЕМ КОРМА, ОБОГАЩЕННОГО ЖЕЛЕЗОМ А.А. ОВСЕПЯН, Н.И. ВЕНЕДИКТОВА, М. В. ЗАХАРЧЕНКО, Р.Е. КАЗАКОВ, М.Н. КОНДРАШОВА, Е.Г. ЛИТВИНОВА, И.Р. СААКЯН, Т.В. СИРОТА, И.Г. СТАВРОВСКАЯ, П.М. ШВАРЦБУРД Институт теоретической и экспериментальной биофизики РАН E-mail: litvinova@rambler.ru,...»

«С.В. Шевчук РУДОЛЬФ ЭРНСТОВИЧ ТРАУТФЕТТЕР (1809–1889) — БОТАНИК, ЗНАТОК ФЛОРЫ РОССИИ Рудольф Эрнстович Траутфеттер (Ernst Rudolf von Trautvetter) родился в г. Митава Курляндской губернии 20 февраля (по новому стилю) 1809 г. в семье старшего преподавателя местной гимназии Эрнста Кристофа фон Траутфеттера — потомственного дворянина. Митава (или Митау) — это старое название Елгавы, второго по величине города Латвийской Республики. Этот красивый старинный и в то же время современный город находится...»

«П.Н. Николаев Михаил Васильевич Ломоносов и развитие физики в Московском университете Москва 2013 2 Михаил Васильевич Ломоносов и развитие физики в Московском университете Николаев П.Н. (Павел Николаевич) Михаил Васильевич Ломоносов и развитие физики в Московском университете. М., 2013. - 125 с. Изучается влияние идей и замыслов М.В. Ломоносова на развитие физики в Московском университете на протяжении всей его истории. Исследуются закономерности в осознании научным сообществом и обществом в...»

«Федеральное государственное автономное образовательное учреждение высшего профессионального образования Северный (Арктический) федеральный университет имени М.В. Ломоносова ПОПОВА Л. Ф. ХИМИЧЕСКОЕ ЗАГРЯЗНЕНИЕ УРБОЭКОСИСТЕМЫ АРХАНГЕЛЬСКА Монография Архангельск – 2014 Дается биогеохимическая оценка экологического состояния и степени загрязнения тяжлыми металлами компонентов природной среды г. Архангельска. Рассматриваются изменение физико-химических свойств почвы и особенности кумуляции,...»

«Библиография новых поступлений отдела абонемента за февраль 2012 г. Естественные науки в целом (20) 20.1 С 17 Самые красивые места мира / [отв. ред. А. Русакова]. — Москва : Мир энциклопедий Аванта+ : Астрель : Полиграфиздат, 2011. — 30, [2] с. : цв. ил. — (Знаменитые достопримечательности планеты). Величественные постройки древности, горные пики и глубочайшие пещеры, клокочущие вулканы и бурлящие водопады, гладь озер, морей и рек разных стран и континентов расположились на страницах этого...»

«Сергей Новиков: Мои Истории. История первая: Семья Новиковых-Келдышей и 20 век. Содержание. Введение: Мои Истории. Мой отец. Происхождение и юность. Ранний большевизм. Мать отца–старообрядка. Начало научной деятельности отца и матери. Стабилизация. Образование семьи Новиковых–Келдышей. Семья Келдышей. Берия и выборы в Академию. Эволюция МГУ в позднем СССР. Берия и Сталин. Перед войной. Троцкий. Война и последующий период. Сталин и создание Израиля. Гулаг. Конец Сталина и убийство Берии....»

«МИНИСТЕРСТВО ГЕОЛОГИИ СССР ПО Южморгеология Научно-исследовательский и проектный институт геофизических методов разведки океана (НИПИОКЕАНГЕОФИЗИКА) ТЕКТОНИКА МЕЗОКАИНОЗОИСКИХ ОТЛОЖЕНИЙ ЧЕРНОМОРСКОЙ ВПАДИНЫ -S8МОСКВА „НЕДРА 1385 УДК 551.24(262.5):551.77 Тектоника мезокайнозойских отложений Черноморской впадины. М., Нед­ ра, 1985, 215 с, ил. (Мингео СССР, ПО Южморгеология, НИПИокеангеофизика). Разрез Черноморской впадины в последние годы изучен сейсморазведкой М О П до глубин 14—15 км, и Черное...»

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования Тверской государственный университет УТВЕРЖДАЮ Декан физико-технического факультета Б.Б. Педько 2012 г. Учебно-методический комплекс по дисциплине ОБЩИЙ ФИЗИЧЕСКИЙ ПРАКТИКУМ. ОПТИКА для студентов 2 курса очной формы обучения направления 010700.62 Физика, специальностей 010801.65 Радиофизика и электроника, 010704.65 Физика конденсированного состояния вещества...»

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования Воронежский государственный педагогический университет УЧЕБНО-МЕТОДИЧЕСКИЙ КОМПЛЕКС ПО ДИСЦИПЛИНЕ Математическая физика Направление подготовки бакалавриат Специальность 050200.62 Физико-математическое образование профиль Информатика Форма обучения очная Срок освоения ООП 4 года Кафедра общей физики Разработчик: Доцент кафедры общей физики...»

«Федеральное агентство по образованию Государственное образовательное учреждение высшего профессионального образования Ухтинский государственный технический университет УТВЕРЖДАЮ Первый проректор _В.К. Хегай УЧЕБНО-МЕТОДИЧЕСКИЙ КОМПЛЕКС по дисциплине Концепции современного естествознания для студентов очной и заочной формы обучения, обучающихся по направлению 0306000 Связи с общественностью и специальности 030602 Связи с общественностью. Ухта Федеральное агентство по образованию...»

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ “ВОРОНЕЖСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ” (ФГБОУ ВПО ВГУ) УТВЕРЖДАЮ Ректор ВГУ _ Д.А. Ендовицкий.2013 ДОПОЛНИТЕЛЬНАЯ ПРОФЕССИОНАЛЬНАЯ ОБРАЗОВАТЕЛЬНАЯ ПРОГРАММА повышения квалификации научно-педагогических работников федеральных государственных учреждений высшего профессионального образования, подведомственных Министерству образования и...»

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

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ ГОСУДАРСТВЕННАЯ КОРПОРАЦИЯ ПО АТОМНОЙ ЭНЕРГИИ РОСАТОМ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ЯДЕРНЫЙ УНИВЕРСИТЕТ МИФИ СНЕЖИНСКИЙ ФИЗИКО-ТЕХНИЧЕСКИЙ ИНСТИТУТ НИЯУ МИФИ НАУЧНАЯ СЕССИЯ НИЯУ МИФИ – 2013 Сборник научных трудов ТРЕТЬЕ ЗАСЕДАНИЕ ТЕМАТИЧЕСКИХ СЕКЦИЙ ПО НАПРАВЛЕНИЮ ИННОВАЦИОННЫЕ ЯДЕРНЫЕ ТЕХНОЛОГИИ 4-6 февраля 2013 г., Снежинск Москва   УДК 001(06) ББК 621.039 Н34 НАУЧНАЯ СЕССИЯ НИЯУ МИФИ-2013. Сборник научных трудов. Третье заседание...»

«Формирование и характеристики пылеугольной летучей золы А.Н. Алехнович Уральская теплотехническая лаборатория, Челябинск Шлакование, загрязнение, абразивный износ поверхностей котлов с огневой стороны, а также надежность и эффективность работы систем газоочистки наряду с параметрами газового потока зависят от физикохимических свойств летучей золы. Для процессов, происходящих в котле в первую очередь важны свойства относительно крупных частиц (микроны). В последние годы вс большее внимание...»

«ДЕТЕКТОРЫ ИОНИЗИРУЮЩЕГО ИЗЛУЧЕНИЯ Гребенщиков В.В., СПбГТУ, кафедра Экспериментальной ядерной физики [Rev.PrB – 04/2002] Электронная версия статьи подготовлена фирмой АВТЭКС Санкт-Петербург, http://www.autex.spb.ru, E-mail: info@autex.spb.ru ДЕТЕКТОРЫ ИОНИЗИРУЮЩЕГО ИЗЛУЧЕНИЯ СОДЕРЖАНИЕ 1. Сцинтилляционные детекторы электромагнитного излучения 1.1. Сцинтилляционные блоки детектирования на основе фотоэлектронных умножителей. 1.2. Сцинтилляционные детекторы с полупроводниковыми фотоприемниками 2....»

«Федеральное агентство по образованию МОСКОВСКИЙ ИНЖЕНЕРНО-ФИЗИЧЕСКИЙ ИНСТИТУТ (ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ) ФИЗИЧЕСКОЕ МАТЕРИАЛОВЕДЕНИЕ В шести томах Под общей редакцией Б. А. Калина Том 2 Основы материаловедения Рекомендовано ИМЕТ РАН в качестве учебника для студентов высших учебных заведений, обучающихся по направлению Ядерные физика и технологии. Регистрационный номер рецензии 012 от 17 декабря 2007 года МГУП Москва 2007 УДК 620.22(075) ББК 30.3я7 К17 ФИЗИЧЕСКОЕ МАТЕРИАЛОВЕДЕНИЕ: Учебник для...»

«5 Биогенное почвообразование под минеральными кочками Вестник Томского государственного университета. Биология 2010 № 4 (12) АГРОХИМИЯ И ПОЧВОВЕДЕНИЕ УДК 574.2:631.484 А.В. Захарченко, Л.А. Изерская, Л.К. Цыцарева, В.Н. Тучак Обособленное структурное подразделение Научно-исследовательский институт биологии и биофизики Томского государственного университета (г. Томск) БИОГЕННОЕ ПОЧВООБРАЗОВАНИЕ ПОД МИНЕРАЛЬНЫМИ КОЧКАМИ ЕMPETRUM NIGRUM В ПСЕВДОФИБРОВЫХ ПОДЗОЛАХ В супесчаных почвах Кольского...»

«Министерство образования Республики Беларусь Учреждение образования Международный государственный экологический университет имени А. Д. Сахарова Факультет экологической медицины Кафедра биохимии и биофизики Е. И. Квасюк С. Б. Бокуть МЕТОДЫ СТРУКТУРНОЙ БИОЛОГИИ В ИЗУЧЕНИИ ФУНКЦИИ БЕЛКОВ ПРАКТИКУМ Минск 2009 УДК 577(075.8) ББК 28.072я73 К32 Рекомендовано к изданию НМС МГЭУ им. А. Д. Сахарова (протокол № 9 от 27 мая 2009 г.) А в то р ы : профессор кафедры биохимии и биофизики УО МГЭУ им. А. Д....»

«ПОДГОТОВКА НАУЧНЫХ КАДРОВ ЦФ РАН ведет активную деятельность по многоуровневой подготовке молодых специалистов - от бакалавриата до аспирантуры. В ЦФ РАН в 1998 году была создана базовая кафедра физики супрамолекулярных систем и нанофотоники Московского физико-технического института (Государственный университет) (МФТИ), созданная согласно приказу ректора по МФТИ от 09.02.1998 г за №26-1, на которой ведется подготовка кадров по современным направлениям исследований в области фотоники...»

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ ГОСУДАРСТВЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ ВОРОНЕЖСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ БИОЛОГО-ПОЧВЕННЫЙ ФАКУЛЬТЕТ ВСЕРОССИЙСКАЯ БИОЛОГИЧЕСКАЯ АССОЦИАЦИЯ СИМБИОЗ РОССИЯ СИМБИОЗ РОССИЯ 2011 Материалы IV Всероссийского с международным участием конгресса студентов и аспирантов-биологов Воронеж, 23–27 мая 2011 Том II Издательско-полиграфический центр Воронежского государственного университета 2011 УДК 57(063) ББК 28.0Я431 С37...»





Загрузка...



 
© 2014 www.kniga.seluk.ru - «Бесплатная электронная библиотека - Книги, пособия, учебники, издания, публикации»

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