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

Автореферат диссертации по теме "Развитие методов интегрированного анализа гравиметрических данных на основе эволюционно-динамических принципов"

На правах рукописи

КУЛЕШОВ ВЛАДИСЛАВ ЕВГЕНЬЕВИЧ

РАЗВИТИЕ МЕТОДОВ ИНТЕГРИРОВАННОГО АНАЛИЗА ГРАВИМЕТРИЧЕСКИХ ДАННЫХ НА ОСНОВЕ ЭВОЛЮЦИОННО-ДИНАМИЧЕСКИХ ПРИНЦИПОВ

Специальность 25 00 16 - Горнопромышленная и нефтегазопромысловая геология, геофизика, маркшейдерское дело и геометрия недр

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

02 онт гт

Ухта 2008

003448365

Работа выполнена в Ухтинском государственном техническом университете

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

Александр Иванович Кобрунов

Официальные оппоненты - доктор физико-математических наук, ведущий научный сотрудник Александр Сергеевич Долгаль

- кандидат физико-математических наук, доцент Смирнов Юрий Геннадьевич

Ведущее предприятие Филиал ООО «ВНИИГАЗ» - «Севернипигаз»,

г Ухта

Защита состоится «23» октября, в Ю00 на заседании диссертационного совета Д 212 29101 в Ухтинском государственном техническом университете по адресу 169300, Республика Коми, г Ухта, ул Первомайская, 13

С диссертацией можно ознакомиться в библиотеке УГТУ Автореферат разослан «19» сентября 2008 г

Ученый секретарь

диссертационного совета УляшеваН М

кандидат технических наук,

профессор

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

Актуальность работы

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

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

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

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

Их решение лежит в сфере создания и развития новых информационных технологий для методов разведки месторождений полезных ископаемых

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

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

Цель работы

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

Основные задачи исследований

1 анализ методов инверсии, направленных на реконструкцию сложно-построенных моделей сред,

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

3 анализ и классификация типовых тектоно-физических элементов Тимано-Печорской провинции (11111) в форме и содержании, обеспечивающих их включение в автоматизированную интерпретационную систему,

4 выполнение вычислительных экспериментов для типовых (эталонных) геолого-геофизических моделей,

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

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

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

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

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

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

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

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

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

Практическая ценность

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

Реализация результатов работы

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

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

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

Результаты исследований нашли свое отражение и представлены в научно-исследовательских отчётах по госбюджетной и хоздоговорной тематике УГТУ и КРО РАЕН, выполнены в т ч по заявкам ОАО «Газпром», ОАО «Газфлот» и других газо- и нефтедобывающих предприятий Результаты, полученные при выполнении работы,

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

Апробация работы

Основные положения диссертации докладывались на Международном семинаре им Д Г Успенского «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей» (Москва, 2007 г, Ухта 2008 г), Международной научно-практической конференции-выставке «Санкт-Петербург-2008» (Санкт-Петербург, 2008 г ), Межрегиональной научно-практической конференции «Освоение минеральных ресурсов Севера проблемы и решения» (Воркута, 2006 г), Международной научно-практической конференции «Геомодель-2006» (Геленджик, 2006 г ), Международной молодежной научной конференции «Севергеоэкотех» (Ухта, 2006 г), научно-технической конференции УГТУ (Ухта, 2006 г, 2007 г), научной конференции «Перспективы нефтегазоносности малоизученных территорий севера и северо-востока европейской части России» (Сыктывкар, 2007 г ), Региональной научно-практической конференции «Проблемы управления в XXI веке» (Ухта, 2007 г), научных семинарах Ухтинского государственного технического университета и КРО РАЕН

По теме диссертации опубликовано 15 работ

Структура и объем работы Диссертация состоит из введения, 5 глав, заключения, списка литературы из 190 наименований, содержит 115 страниц текста, включая 52 рисунка

Благодарности:

Автор выражает искреннюю благодарность своему научному руководителю, доктору физико-математических наук, профессору, академику Российской академии естественных наук А И Кобрунову за постановку задачи, оказание огромной помощи и постоянный контроль на всех этапах научно-исследовательской работы, а также кандидату технических наук, доценту С В Шиловой за внимание и помощь в практической реализации научных исследований

Автор признателен ректору УГТУ, профессору, доктору технических наук Н Д Цха-дая и зав кафедрой геологии нефти и газа, доценту А Н Смирнову за создание оптимальных условий, понимание и постоянное внимание к работе

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

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

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

В первой главе проведен комплексный анализ методов и технологий инверсии, прежде всего гравиметрических данных, направленных на реконструкцию сложно-построенных моделей геологических сред Основываясь на работах В Н Страхова (1970-2007) автор рассматривает этапы развития теории и практики интерпретации потенциальных полей Приведена их классификация, периодизация, подведены основные итоги Развитие информационного геофизического обеспечения и информационно-статистической теории методов инверсии прослежены по работам Ф М Гольцмана и Т Б Калининой (1997), А А Никитина (1999), Г И Каратаева (1996) Выполнен анализ методов инверсии, направленньсх на реконструкцию сложно-построенных моделей сред в условиях эквивалентности и скрытой эквивалентности (А И Кобрунов, 1999)

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

В соответствии с классификацией, предложенной, В Н Страховым, а так же изложенной в работах А. И Кобрунова (1978-2008), выделяется четыре типа задач, по мере развития ставившихся перед геофизическими методами и перед гравиразведкой в частности

1 обнаружение геологических объектов (новых объектов в предварительно сформированной геологической модели) - совсем простые модели да - нет,

2 пространственная локализация обнаруженных объектов - простые модели, эквивалентность отсутствует,

3 расчленение геологического разреза (или его фрагмента) - сложные модели, возникает явление скрытой эквивалентности по А И Кобрунову и в -эквивалентность,

4 детальное описание геологического разреза (или его фрагмента) - сложные модели с заведомо включенными в себя областями теоретической эквивалентности

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

В работе рассмотрены методы, которые на уровне постановок задач реализуют эти подходы

Прежде всего, следует отметить большое число методов инверсии, основанных на введении ограничений на искомые значения параметров, характеризующих модель среды Эти работы восходят к А А Юнькову (1958), Е Г Булаху (1973), С В Шалаеву (1972), В И Старостенко (1978), А А Непомнящих (1976) и продолжены большим числом авторов (см известную работу В Н Страхова «Общая схема и основные итоги развития теории и практики потенциальных полей в СССР и России в XX веке», опубликованную в трудах конференции «Развитие гравиметрии и магнитометрии в XX веке») Аргументом для ведения параметризаций и ограничений на значения параметров в этих методах, которые чаще всего называют методами призматической аппроксимации и аппроксимационным подходом, служит возможность достаточно хорошей аппроксимации моделью подразумеваемой среды Ограничение, которое тормозит развитие этих подходов, состоит в том, что в них проявляется

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

Альтернативный подход, который выше охарактеризован как определение допустимых комбинаций параметров модели, реализован в принципах использования критериев оптимальности для параметров среды при постановке задач инверсии Этот подход состоит в следующем Формируется функционал J(a), оценивающий качество того или иного из элементов выбранной многопараметрической модели среды, либо сам функционал J{cj) включает в себя конструкцию невязки полей Далее решается задача максимизации этого качества при условии, что поле от а соответствует заданному (с наперед заданной точностью) Эта задача на условный экстремум решается либо прямыми методами, что крайне неудобно, либо строится уравнение Эйлера (необходимые условия экстремума), которому решение должно удовлетворять Это уравнение и определяет те допустимые комбинации параметризаций, которые ограничивают область эквивалентности Далее необходимо из набора этих комбинаций всего лишь выделить тот элемент, который с требуемой точностью создает как раз необходимое поле Этот путь имеет, несмотря на кажущуюся простоту формулировки, очень много различных модификаций, связанных прежде всего с тем, как формируется функционал J(a), и какой смысл в него вкладывается Это также зависит от того, является ли исходный модельный класс классом единственности или эквивалентности К методам подобного типа относятся методы регуляризации при инверсии гравиметрических данных (В Н Страхов, В И Старостенко, В П Танана, В В Васин, А Н Тихонов, М М Лаврентьев и многие др (См статью В Н Страхова, о которой упоминается выше)) Центральной проблемой, которая возникает в рамках такого подхода, является проблема выбора критерия оптимальности J(a) Свойства решения предопределены этим выбором Оказывается, что для сложных и очень сложных моделей сред, казалось бы, разумные и хорошо аргументированные приемы выбора этого функционала (например, принцип наименьших квадратов или коэффицент корреляции) приводят к парадоксальным и бессмысленным в геологическом отношении результатам Свойства этих решений исследовались А И Кобруновым (1976-1990), В II Страховым (1977-1980), В И Старостенко (1976-1978), С М Оганесяном (1981-1987), А С Маргу-лисом (1984), В М Новосилицким (1981) Возникла необходимость создания методики выбора этих критериев как реальных носителей геолого-геофизической информации Здесь следует отметить работы А И Кобрунова (1979-2007), В А Варфоломеева (1981), Р П Де-нисюка (1981), А П Петровского (1997), С В Моисеенковой (1998), С В Шиловой (2000) и многих других Слабая управляемость геологической содержательностью результата, исклю-

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

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

лиг создать более полную структуру учета генезиса при анализе геофизических данных, состоящую из формализованных фрагментов, как элементов общей структуры Такой путь развивался, в частности, В П Мясниковым, В О Михайловым, В М Гординым, Е П Тимош-киной и др , (1979-2007) Моделирование динамических и эволюционно-динамических процессов опирается на известные параметры вязкости, которые входят, например, в уравнение Навье-Стокса Это известное уравнение вязких течений, получаемое из закона сохранения массы Решение этого уравнения неустойчиво относительно вариаций параметров вязкости В условиях их приближенного задания получаемые результаты чаще следует рассматривать как качественные, оценочные Однако, само эволюционное уравнение может служить ограничивающим фактором эквивалентности для задачи инверсии, обеспечивающим учет генезиса объекта на некотором эффективном «отрезке времени» В этом понимании введение эволюционно-динамических уравнений в процесс инверсии геофизических, в частности гравиметрических данных, аналогично принципу использования внешних критериев оптимальности относительно параметров среды Оно также дает характеристику «допустимого многообразия» решений, следующего на этот раз из эволюционных уравнений (а не уравнений Эйлера для вариационной задачи), на котором и следует выбрать элемент - предельный элемент «эволюции», характеризущийся тем, что поле от него есть требуемое Этот принцип также допускает многочисленные интерпретации Например, его можно интерпретировать, как введение для эволюционных уравнений краевых условий - соблюдение заданных геофизических полей (В И Старостенко, 2004) Все сказанное требует внимательного рассмотрения самих геодинамических принципов как признаков нефтегазоносности (гл 2), а так же процедур их включения в интерпретационный процесс (гл 4)

Во второй главе «Геодинамика и нефтегазоносность» рассмотрены принципы введения итогов геодинамического анализа при построении моделей нефтегазоносных бассейнов с использованием геофизических данных, в частности гравиметрических Ключевую роль в методах геодинамического анализа играет оценка сбалансированности по тем или иным принципам (физико-механическим, сейсмическим) физико-геологических моделей и внесение в нее минимальных балансировочных корректив (А И Кобрунов, С В Шилова, В Е Кулешов, 2005)

Обобщая проведенный в работе анализ методов решения обратных задач на основе эволюционно-динамических подходов, мы выделяем последовательность, составляющую их классификацию, которая служит конструктивной основой методов инверсии и выбора методики анализа данных (см. Рис 1 А И. Кобрунов, В Е. Кулешов, 2007)

Постановки задач инверсии, поясняющие эволюционно-динамические методы из приведенных на рис 1, даны в уравнениях (1-5) на рис 2 Пунктиром на рис 1 приведено

ясировакных

чЯШЖ

Рис. 1. Схема эволюционно-динамических подходов интегрированного анализа геофизических

данных.

Метод комплексных минимальных корректив

(1)

Метод ш иг спальных корректив

./(д(х)-77(х))-> тт.

./(с(х) - тш.

Метод дротпгческого продолжения

1ипА(,¿(/.у)) = 1Ця).

4'+1(У) = п' (V) + ос1а^срА'{,ч)1 77°(У) = С0(У) = ,7(У). <РВ{*) = Щ1(У))-У(.*)

(5)-

щ?

■ ¡-№1 гО

(V) = (у) + а, Су1 О)],

(V) = ??(У).

Рис. 2. Формулировки задач инверсии в эволюционно-динамических методов интегрированного анализа.

видимое дальнейшее развитие методов, которое не излагается, поскольку не включено в защищаемые положения Более подробное описание приведено в [1] Пояснения к обозначениям следующие г] - распределение физических параметров, для которых решается обратная задача, и, у - наблюдаемые геофизические поля, J - критерий оптимальности, Л(£(х)),В(;7(х)) - операторы решения прямой задачи (сейсморазведка, гравиразведка) для распределений параметров £(х)),(//(х)) ./( )- функционал, оценивающий меру удаления двух моделей Е[£(/,у)] = 0- эволюционное уравнение для параметра £(/,V)относительно эффективного «времени» I развития системы, которое начинается в момент /0 и продолжается до тех пор, пока не окажется соблюденным условие современного поля а - параметр релаксации, (? -сверочный оператор, связанный с производными от операторов А и В

В рамках идеи сбалансированных моделей мы выделяем три метода (см рис 1), третий из которых (МДП) включается также в идеологию моделирования динамики в соответствии со сформулированным утверждением о том, что принцип балансировки означает проверку согласованности с заданным физическим полем либо физическими законами Методы минимальных корректив состоят в нахождении минимально уклоняющегося заданного критерия оптимальности J от построенной модели >] распределения с(см рис 2) Этот принцип работает только для монометодов В том случае, когда работает система методов, мы вводим идею комплексных минимальных корректив, которая состоит во взаимной увязке моделей для разных методов Взаимная увязка имеет отличительную особенность, состоящую в том, что требуется нахождение таких распределений различных физических параметров, каждое из которых, соблюдая свое поле, менее всего уклоняется друг от друга Принцип наименьшего уклонения друг от друга позволяет отойти от поиска решения, одновременно удовлетворяющего и одному и другому полю Это положительное качество, в связи с тем, что такого решения просто не существует Метод динамического продолжения (комплексного динамического продолжения) по своей формульной записи использует эволюционные уравнения и реализует, как уже говорилось выше, построение такого распределения из семейства удовлетворяющих заданному уравнению динамики, которое в пределе обеспечивает сопротивление заданного поля Однако реализация метода сводится к итерационному процессу, аналогично тому, который использовался в методе минимальных корректив и в методе комплексных минимальных корректив Отличие состоит лишь в том, что оператор в, связанный с видом критерия оптимальности в последнем методе является динамически меняющимся от шага к шагу

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

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

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

В работах Хаина В Е, Клещева К А , Шеина В С, и др (1990) обоснован ряд новых теоретических положений, и объединены разработки в области применения теории тектоники плит в нефтегазовой геологии Это позволяет создать геодинамическую основу прогноза и поисков нефти и газа, разработать модели строения, эволюции НГБ, рассмотреть геологию и нефтегазоносность крупных регионов, например арктических акваторий (Гаврилов В П и др), юга бывшего СССР (Клещев К.А, Шеин ВС), территорий и акваторий Российской Федерации в целом Каждое из указанных направлений по-своему способствует повышению достоверности прогноза нефтегазоносности Практическое направление приближено к решению задач прогноза и поисков нефти и газа

Клещев К А и Шеин В С выделяют в своих работах составляющие современной геодинамической основы прогноза поисков и разведки нефти и газа

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

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

логическом районировании при оценке перспектив нефтегазоносности На основе такого районирования (Астафьев С В , 1997, Зоненшайн Л П, 1999, Игнатова В А, 2003, Пешкова И Н , Соколов Б. А, 1990 и др) выделеио пять крупных групп бассейнов, характеризующих их эволюцию I - континентальных и межконтинентальных палеорифтов и надрифтовых прогибов, депрессий, II - пассивных, трансформных континентальных палеоокраин и краевых прогибов, III - океанических рифтов, IV - орогенов столкновения плит, V - областей схождения плит (активных континентальных палеоокраин, островных дуг, окраинных морей)

Тимано-Псчорский НГБ относится ко П группе—пассивной окраины (Я-СО ВосточноЕвропейского папеоконтинента и краевого прогиба (Р-Т1), трансформированных столкновением с палеоконтинентом Баренция (V), с палеомикроконтинентами, островными дугами (С2> С,-Т) и изостазией (Т3-С>)

Анализ распределения УВ-ресурсов мира, образованных в разных геодинамических обстановках, осуществленный Геологической службой США (Ульмишек Г и др, 2004), показал, что наибольшими ресурсами обладают бассейны пассивных континентальных палеоокраин и предорогенных прогибов

Т о, в главе 2 дано обоснование привлечения информации геодинамического характера (эволюция развития НГБ) в оволюционно-динамическом подходе как альтернативы сложным математическим уравнениям, в частности уравнениям Навье-Стокса.

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

Геодинамическая эволюция истории ТПП рассматривает характер и последовательность развития во времени плитотектонических процессов и явлений в геодинамическом цикле и в стратиграфическом разрезе

Тимано-Печорская литосферная плита принимается в объеме от Западно-Тиманского разлома до Главного Уральского надвига на востоке (В В Юдин, В А Дедеев, 1987, В Н Пучков, 1997, Н И Тимонин, 1998, Н А. Малышев, 1998, А И Дьяконов, 1999 и др ) Прослежена эволюция Тимано-Печорской провинции и процессы, её сопровождающие Построены модели развития ТПП во времени. (Рис 3)

а) рифей - венд

б) кембрий

'■а.!?<4!х!ж кфушж»

Рис.3. Модели развития Тимано-Печорской провинции. (По материалам В. В. Юдина, В. А.

Дедеева, В. Н. Пучкова, Н. И. Тимонина, Н. А. Малышева, А. И. Дьяконова и др.)

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

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

В качестве примера локальной структуры из состава ТПП исследована схема формирования гряды Ченышева. Исходя из приведённых данных, сделаны выводы о структурной эволюции гряды по материалам А. Б. Кирмасова (2002).

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

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

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

Базовыми кинематическими типами глубинных региональных разломов ТПП являются сбросы и взбросы (Н А Малышев, А И Дьяконов) Например, на участке региональных геологических разрезов, проходящих через Припечорскую и Шапкинскую зоны разломов, наблюдаются сбросы и взбросы, форма и простирание которых повторяется и на локальном уровне в районе Шапкинского и Южно-Шапкинского месторождений

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

Другим примером из состава ТПП служат брахиантиклинальные структуры, которые выделяются на региональных профилях (22-РС) на участке Ижма-Печорской и ВерхнеПечорской впадины На локальном уровне их форма повторяется в Пашнинской, Мичаю-ской, Восточно- и Северо-Савиноборской структурах Аналогичная ситуация наблюдается для Западно-Тэбукского месторождения Здесь форма брахиантиклинальных структур также схожа с формой региональных структур Омра-Сойвинской и Тэбукской ступеней (Региональный профиль КМПВ 1-Ф) Однако в нижней части разреза Западно-Тэбукского месторождения форма антиклиналей изменяется и не совпадает с локальными структурами

Четласская горст-антиклиналь имеет северо-западное простирание Она образована двумя крупными разломами Блок выдвинут и смят в антиклинальную складку Эта антиклиналь отражается на крыльях аналогичными складками более высоких порядков северовосточного простирания до микроскладчатости (наблюдаются складки размерами до нескольких сантиметров) Сама складка является структурой первого порядка, и на обеих крыльях выделяется второй порядок и более Здесь мелкие складки подчеркивают крупную (А М Плякин, В М Пачуковский, 1968)

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

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

Эволюционно-динамический подход нашёл своё применение в программных продуктах GCIS, SIGGRA (включает в себя два независимых модуля SIGGRA(s) и S1GGRA (g)) и VIP. Совокупность и стыковка программ в единую технологию интегрированной интерпретации обеспечивается посредством форматов ввода-вывода. (Рис. 4). В главе сделан акцент на методический аспект получения, корректировки и балансировки моделей различного строения.

Автоматизированная система GCIS (Geophysical Complex Interpretation System - система комплексной интерпретации геофизических данных, авторы А. П. Петровский, А. И. Кобрунов и др.)

1) обеспечивает:

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

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

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

Рис. 4. Технология обработки ГГИ

• построение согласованных сейсмо-гравиметрических моделей 2) включает в себя процедуры

• решение эффективной кинематической задачи сейсморазведки на оптимизационных принципах,

• построения оптимальных структурно-плотностных моделей,

• построения локальных плотностных неоднородносгей по некомпенсированной структурной модели компоненты гравитационного поля,

• согласования скоростных и плотностных моделей

Более подробное изложение и состав системы приведены в работе [12] Система VIP (авторы Д Г Григорьев, А И Кобрунов) в рамках технологии является интегрирующим модулем и обеспечивает импорт данных и параметров из системы GCIS, визуализацию и редактирование данных

Программный комплекс SIGGRA (автор А И Кобрунов, реконструкция С С Антонов, В Е Кулешов) включает в себя два независимых модуля SIGGRA(s) и SIGGRA (g) которые ориентированы на реконструкцию среды в терминах модели распределений плотности (SIGGRA (s)) и модели среды в терминах системы плотностных границ (SIGGRA (g)), на основе принципов гравитационной балансировки Их комбинация - интерактивное исполнение, которое позволяет изучать модели, содержащие как структурные элементы, так и зоны локальных плотностных аномалий Некоторые функции инверсии, предусмотренные в комплексе GSIS, дублируются в комплексе SIGGRA Однако в отличие от GSIS, ориентированного на совместное решение обратных задач сейсмогравиметрии, комплекс SIGGRA ориентирован на методы гравитационной инверсии с включением эволюционно-динамических принципов Сращивание этих комплексов в единое целое при одновременном возможном включении внешних математических библиотек для выполнения процедур специальной обработки данных выполнено за счет программного комплекса VIP Программные комплексы SIGGRA обеспечивают

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

• учет влияния рельефа,

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

(gran).

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

• возможность создания "синтезированной" структурно-плотностной модели, полученной из результатов решения обратных задач в классе плотностных моделей (БККЖА^)) и структурных (ЭЮОНА

• возможность учёта при интерпретации эволюционно-динамических принципов.

После выделения элементов, характерных для строения и эволюции региона, эта информация используется в качестве критерия оптимальности в эволюционно-динамическом подходе. С помощью программных комплексов ЗКХЖА (в) и 5ЮСЖА (¿) в исходную модель вносятся коррективы. Границы пошагово изменяются согласно введённым данным. Пример внесения пошаговых изменений модели приведён на рис.5. Вначале задаётся модификатор (имитатор модели) в соответствии с типовым элементом строения региона, к которому относится изучаемая модель. Затем выделяется участок модели и изменяется согласно заданному модификатору.

Выбор необходимого модификатора

р| €

¡¡ш

1-' .........".................~~ 1———-—1 1

.....:> 1!

Выделение участка границы для внесения изменений

Изменение границы согласно выбранному модификатору

Варианты типовых элементов, характерных для строения региона (модификаторы)

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

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

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

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

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

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

Аналогичным образом с использованием метода МКМК проведено моделирование участка месторождения «Геофизическое» в районе Карского моря

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

Заключение

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

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

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

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

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

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

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

Список опубликованных работ:

1 Кулешов В Е, Кобрунов CAO методе балансировки при поисках мелких и слабоконтрастных месторождений нефти и газа // 4-я межрегиональная научно-практическая конференция «Освоение минеральных ресурсов Севера проблемы и решения» (12-14 марта 2006 г, Воркута) Труды конференции - Сыктывкар Изд-во СЛИ, 2006 г - С 608-613

2 Кулешов В Е Возможности применения технологий, основанных на комплексной интерпретации данных, для решения задач инженерной геофизики // Вторая международная научно-практическая конференция «Инженерная геофизика -2006» (17-22 апреля 2006 г , Геленджик) Тезисы докладов - Москва Изд-во МГУ, 2006 г - С 42-43

3 Кулешов В Е, Кобрунов А И, Петровский А П, Шилова С В, и др Технологии интегрированного эволюционно-динамического анализа гравиметрических данных И УИ-ая международная научно-практическая конференция «Геомодель-2006» (17-22 сентября 2006 г, Геленджик) Тезисы докладов -Москва Изд-во МГУ, 2006 г - С 239-240

4 Кулешов В Е, Кобрунов Л И, Петровский А П, Шилова С В, и др Теория и методы интегрированного эволюционно-динамического анализа гравиметрических данных // УН-ая международная научно-практическая конференция «Геомодель-2006» (17-22 сентября 2006 I , Геленджик) Тезисы докладов -Москва Изд-во МГУ, 2006 г -С 175-176

5 Кулешов В Е, Кобрунов А И, Петровский А П.Шилова С В.идр Принципы и подходы гравитационной балансировки при интегрированном эволюционно-динамическом анализе гравиметрических данных // УИ-ая международная научно-практическая конференция «Геомодель-2006» (1722 сентября 2006 г, Геленджик) Тезисы докладов -Москва Изд-во МГУ, 2006 г - С 176-178

6 Кучешов В Е Интегрированные технологии интерпретации и возможности их взаимодействия // Международная молодежная научная конференция «Севергеоэкотех - 2006» (17-19 марта 2006 г, Ухта) Материалы конференции - Ухта Изд-во УГТУ, 2006 г С 99-104

7 Кулешов В Е Кобрунов А И, Шилова С В, Григорьев Д Г Результаты предварительной интерпретации сейсмогравиметрических данных с применением интегрированного комплекса программ (на примере шельфа Карского моря) // VII научно-техническая конференция (18-26 апреля

2006 г, Ухта) Сборник научных трудов -Ухта Изд-во УГТУ, 2006 г С 45-49

8 ВелътистоваО М.Шилова С В, Кучешов В Е, Урбан А В Моделирование сейсмогеоло-гических разрезов Лемвинской екпадчато-надвиговой зоны на основе комплексной интерпретации геолого-геофизических данных // 34-я сессия Международного семинара им Д Г Успенского «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей» (29 февраля - 3 февраля 2007 г , Москва) Материалы семинара - Москва Изд-во ИФЗ РАН,

2007 г С 59-62

9 Кулешов В Е, Кобрунов А И, Антонов С С Интегрированные технологии, обеспечивающие эволюционно-динамический анализ геолого-геофизических данных // 34-я сессия Международного семинара им Д Г Успенского «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей» (29 февраля - 3 февраля 2007 г, Москва) Материалы семинара -Москва Изд-во ИФЗ РАН, 2007 г С 137-140

10 Кулешов В Е, Шилова С В, Коксов В В Результаты объемного сейсмогравитационного моделирования в пределах акватории Обской губы Карского моря // 34-я сессия Международного семинара им Д Г Успенского «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей» (29 февраля - 3 февраля 2007 г, Москва) Материалы семинара - Москва Изд-во ИФЗ РАН, 2007 г С 286-289

11 Кулешов В Е. Кобрунов А И, Шилова С В Технология и методика интегрированной интерпретации сейсмогравиметрических данных // Научная конференция «Перспективы нефтегазоност-ности малоизученных территорий севера и северо-востока европейской части России» (9-10 октября 2007 г, Сыктывкар) Тезисы докладов - Москва Изд-во ИГиРГИ, 2007 г - С 69-71

12 Кулешов В Е, Шилова С В, Урбан А В Функциональные особенности программного обеспечения автоматизированного комплекса ССК // Региональная научно-практическая конференция «Проблемы управления в XXI веке»(16-18 октября 2007 г, Ухта) Материалы конференции -Ухта Изд-во Института управления, информации и бизнеса, 2007 г - С 25-29

13 Кулешов В Е, Кобрунов А И, Шилова С В Повышение достоверности построения геолого-геофизических моделей на поисково-разведочном этапе работ на нефть и газ // Нефтяное хозяйство -2007 -№10 С 64-67

14 Кулешов В Е Развитие новых подходов в эволюциошю-динамических методах и технологиях интегрированного анализа ГГИ // 35-я сессия Международного семинара им Д Г Успенского «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей» (29 февраля - 3 февраля 2008 г, Ухта) Материалы семинара - Ухта 2008 г С 169177

Отпечатано в отделе оперативной полиграфии Ухтинского государственного технического университета Республика коми, г Ухтаул Октябрьская, 13 Уел печ л 1,4, уел изд л 1,1, тираж 150 экз Сдано в печать 15 09 2008г Заявка № 783

Содержание диссертации, кандидата технических наук, Кулешов, Владислав Евгеньевич

ВВЕДЕНИЕ

1 РАЗВИТИЕ ИНТЕРПРЕТАЦИОННОГО ГЕОФИЗИЧЕСКОГО ОБЕСПЕЧЕНИЯ

1.1 Развитие теории интерпретации геофизических данных

1.2 Развитие информационно-статистической теории

1.3 Эволюция задач, стоящих перед геофизическими методами

1.4 Комплексная(интегрированная) интерпретация геофизических данных '

1.5 Эволюционно-динамические подходы интегрированного анализа геофизических данных

2 ГЕОДИНАМИКА И НЕФТЕГАЗОНОСНОСТЬ

2.1 Задачи и методы геодипамического анализа

2.2 Применение геодинамического анализа па различных этапах ГРР

2.3 Нефтегазогеологическое районирование па основе геодинамики

3 ГЕОДИНАМИЧЕСКИЕ МОДЕЛИ РАЗВИТИЯ РЕГИОНОВ

3.1 Модели геодинамического развития ТПП

3.2 Геодинамическая история развития ТПП

3.3 Характеристика и эволюция разломных структур ТГ1П

3.4 Структура формирования гряды Чернышева

3.5 Геодинамические модели Балтийского щита

3.6 Использование идей фрактальности в геологии при формировании модельных представлений и систематизации материала

4 ПРОГРАММНО-МЕТОДИЧЕСКОЕ ОБЕСПЕЧЕНИЕ ИНВЕРСИИ С ИСПОЛЬЗОВАНИЕМ ЭВОЛЮЦИОННО-ДИНАМИЧЕСКИХ ПРИНЦИПОВ

4.1 Характеристика программного обеспечения эволюционно-динамического подхода и выбор параметров управляющих моделированием

4.2 Автоматизированная система GSIS

4.3 Система VIP '

4.4 Программный комплекс SIGGRA

4.5 Формирование типовых динамических моделей

5 ПРИМЕРЫ ИСПОЛЬЗОВАНИЯ ТЕХНОЛОГИЙ, ОСНОВАННЫХ НА ЭВОЛЮЦИОННО-ДИНАМИЧЕСКИХ ПРИНЦИПАХ ДЛЯ ИЗУЧЕНИЯ НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ

5.1 Геофизическое нефтегазовое месторождение

5.2 Хамботейско-Бухаринский участок

5.3 Юньяхинская площадь

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

Актуальность работы

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

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

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

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

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

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

Цель работы

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

Основные задачи исследований

1. анализ методов инверсии, направленных на реконструкцию сложно-построенных моделей сред;

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

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

4. выполнение вычислительных экспериментов для типовых (эталонных) геолого-геофизических моделей;

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

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

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

2. сформулированы и обоснованы принципы гравитационной балансировки сложно-построенных геологических сред;

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

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

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

3. включение принципов эволюционно-динамического моделирования в процесс инверсии повышает достоверность построенных моделей.

Практическая ценность

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

Реализация результатов работы

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

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

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

Результаты исследований нашли своё отражение и представлены в научно-исследовательских отчётах по госбюджетной и хоздоговорной тематике УГТУ и КРО РАЕН, выполнены в т. ч. по заявкам ОАО «Газпром», ОАО «Газфлот» и других газо- и нефтедобывающих предприятий. Результаты, полученные при выполнении работы, использованы в учебном процессе при подготовке специалистов и магистров по полевой нефтегазовой геофизике и геологии нефти и газа.

Апробация работы

Основные положения диссертации докладывались па Международном семинаре им. Д. Г. Успенского «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей» (Москва, 2007 г., Ухта 2008 г.), Международной научно-практической конференции-выставке «Сапкт-Петербург-2008» (Санкт-Петербург, 2008 г.), Межрегиональной научно-практической конференции «Освоение минеральных ресурсов Севера: проблемы и решения» (Воркута, 2006 г.), Международной научно-практической конференции «Геомодель-2006» (Геленджик, 2006 г.), Международной молодежной научной конференции «Севергеоэкотсх» (Ухта, 2006 г), научно-технической конференции УГТУ (Ухта, 2006 г., 2007 г.), научной конференции «Перспективы нефтегазоносности малоизученных территории севера и северо-востока европейской части России» (Сыктывкар, 2007 г.), Региональной научно-практической конференции «Проблемы управления в XXI веке» (Ухта, 2007 г), научных семинарах Ухтинского государственного технического университета и КРО РАЕН. По теме диссертации опубликовано 14 работ.

Структура и объем работы. Диссертация состоит из введения, 5 глав, заключения, списка литературы из 185 наименований, содержит 141 страницу текста, включая 57 рисунков.

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

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

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

Библиография Диссертация по наукам о земле, кандидата технических наук, Кулешов, Владислав Евгеньевич, Ухта

1. Актуальные научно-технические проблемы развития геолого-геофизических, промысловых и поисково-разведочных работ в Республике Коми: монография. - Ухта: КРО РАЕН, 2001. - 372 св.: ил.

2. Андреев Б. А., Кулишин И. Г. Геологическое истолкование гравитационных аномалий. М.: Гостоптехиздат, 1962.

3. Аронов В. И., Ланда Т. И. Система программ обработки на ЭВМ гравитационных и магнитных аномалий, заданных на плоскости или произвольной поверхности. М.: ВИЭМС, 1972. (Сер. регион., развед., промысл, геофизика).

4. Балагапский В. В., Глазнев В. II., Осипенко Л. Г. Раннепротерозойекая эволюция северо-востока Балтийского щита: анализ тектоно-стратиграфических террейнов //Геотектоника. 1998. 2. С. 16-28.

5. Балк П. И. Детерминистические модели интерпретации гравитационных полей. Диссертация на соискание уч. Ст. Доктора физ. мат. наук. Киев, 1989 г.

6. Балк П. И., Балк Т. В., Горчаков И. В. Об устойчивости решения обратной задачи гравиметрии для группы точечных источников // Геология и геофизика, 1982,- № 10.- С.118-126

7. Банковский В. И., Захаров В. С. Динамические процессы в геологии: первое знакомство с нелинейными системами: Электронная версия учебного пособия, (http://www.dinamo.geol.msu.ru/DPG.html).

8. Бонгард М. М. и др. Решение геологических задач с помощью программы распознавания // Сов. Геология. 1963. № 6.

9. Булах Е. Г. Автоматизированная система интерпретации гравитационных аномалий. К.: Наукова думка, 1973. - 202с.

10. Булах Е. Г., Ржаницын В. А., Маркова М. Н. Применение метода минимизации для решения задач структурной геологии по данным разведки. Киев: Наук, думка, 1976.

11. Воронин Ю. А. Применение в геологоразведке математики па базе ЭВМ: СССР и США // Препр. Новосибирск: ВЦ СО АН СССР, 1987. 26 с.

12. Гаврилов В. П. Геодинамика и нсфтегазоносность Арктики. М.: Недра, 1993 -323 с.

13. Гаврилов В. П. Геодинамическая модель нефтегазообразования в литосфере // Геология нефти и газа. 1998.- №6. - С. 2-12.

14. Гамбурцев Г. А. Об одном способе определения расположения подземных масс на основании магнитных и гравитационных наблюдений // Журн. прикл. физики. 1930. №7, вып.2.

15. Геодинамические основы прогноза и поисков нефти и газа и их внедрение в практику геологоразведочных работ / К. А. Клещев, В. С. Шеип, В. Е. Хаип и др. -М.: ВИЭМС, 1990.

16. Гзовский М. В. Основы текгонофизики. -М.: Наука, 1975. 536 с.

17. Гладкий К. В. Частотный анализ при обработке и интерпретации данных маг-нито- и гравиразведки и возможности автоматического решения задач обработки и интерпретации. Гостоптехиздат, 1960. (Тр. МИНХиГП; вып. 31).

18. Глазнев В. II. Комплексные геофизические модели литосферы Фепноскандии. Апатиты. ЗАО «КаэМ». 2003. - 252 с.

19. Гольдшмидт В. И. Оптимизация процесса количественной интерпретации данных гравиразведки. М.: Недра, 1984. - 185 с.

20. Гольцман Ф. М. Основы теории интерференционного приема регулярных волн. М.: Наука, 1964. 283 с.

21. Гольцман Ф. М. Проблемные вопросы информационно-статистической теории геофизических наблюдений // Изв. АН СССР. Физика Земли. 1977. № 12. С. 75-86.

22. Гольцман Ф. М. Развитие информационно-статистических методов решения обратных задач при интерпретации геоданпых в 20 веке // Развитие гравиметрии и магнитометрии в 20 веке М.: ОИФЗ РАН, 1997. - 234 с.

23. Гольцман Ф. М. Состояние и перспективы развития статистических методов интерпретации геофизических наблюдений // Разведочная геофизика на рубеже 70-х годов. М.: Недра, 1974.

24. Гольцман Ф. М. Статистические модели интерпретации. М.: Наука, 1971. 374 с.

25. Гольцман Ф. М., Калинина Т. Б. Априорная неопределенность выбора гсолого-геофизических моделей и эффективность интерпретации гравитационных и магнитных полей // Геология и геофизика. 1991. № 6. С. 127-134.

26. Гольцман Ф. М., Калинина Т. Б. Комплексировапис качественных и количественных признаков в задачах геолого-геофизической интерпретации //Изв. РАН. Физика Земли. 1993. №8. С. 88-93.

27. Гольцман Ф. М., Калинина Т. Б. Моделирование трехмерных геологических объектов по комплексу геоданпых в условиях априорной неопределенности // Изв. РАН. Физика Земли. 1996. №2, с. 82-89.

28. Гольцман Ф. М., Калинина Т. Б. Простые методы частотного анализа и синтеза сигналов и их применения к решению некоторых геофизических задач // Прикл. геофизика. 1958. Вып. 21.

29. Гольцман Ф. М., Калинина Т. Б. Статистическая интерпретация магнитных и гравитационных аномалий. Л.: Недра, 1983. 247 с.

30. Гольцман Ф. М., Калинина Т. Б. Статистический алгоритм интерпретации и его применение для решения обратной задачи магниторазведки // Изв. АН СССР. Физика Земли. 1967. №4.

31. Гольцман Ф. М., Калинина Т. Б., Калинин Д. Ф. Статистическая методология построения моделей геолого-геофизических объектов по комплексу геоданных // Рос. геофиз. жури. 1994. №3/4. С. 61-66.

32. Губермап Ш. А. Использование обучающих программ для решения геологических задач // Комплексная интерпретация геологических и геофизических данных па вычислительных машинах. М.: Гостехиздат, 1966. (Тр. МИНХиГП: вып. 62).

33. Дьяконов А. И. Геодинамика. Учеб. пособие. Ухта: УГТУ, 2000. - 67 с.

34. Дьяконов А. И., Белый Н. И. Тектонодинамичские критерии раздельного прогноза нефтсгазоносности. М.: ВЫИИОЭНГ, 1993 - 125 с.

35. Жданов М. С., Шрайбман В. И. Корреляционный метод разделения геофизических аномалий. М.: Недра, 1973.

36. Зоненштайн Jl. П., Кузмин М. И., Натагюв Л. М. Тектоника литосферных плит территории СССР. М.: Недра, 1990.

37. Зоненштайн Л. П., Савостин Л. А. Введение в геодинамику. — М.: Недра, 1979.-310с.

38. Иванов В. К. О линейных некорректных задачах // Докл. АН СССР. 1962. Т. 145, №2.

39. Казанский А. П. Определение элементов залегания полезного ископаемого по данным магниторазведки // Изв. АН СССР. Сер. геофиз. 1938. 138 с.

40. Калинина Т. Б., Гольцман Ф. М. Развитие информационно-статистических методов решения обратных задач в гравиразведке и магниторазведке в 20 веке // Развитие гравиметрии и магнитометрии в 20 веке — М.: ОИФЗ РАН, 1997. 234 с.

41. Каратаев Г. И. Корреляционная схема геологической интерпретации гравитационных и магнитных аномалий. Новосибирск: Паука, 1966. 135 е.

42. Карта пефтегазогеологического районирования СССР. Под ред. Г. А. Габри-элянца. 1990.

43. Кесарев В. В. Эволюция вещества вселенной. — М.: Атомиздат, 1976. 86 с.

44. Кирмасов А. Б. Структурная позиция мезозойских трубок щелочно-ультраосновных пород гряды Чернышева, Республика Коми, (http://geo.web.ru/geolab/publ/paper/rudyl/index.html).

45. Клещев К. А. Геодинамика и новые типы природных резервуаров для нефти и газа / К. А. Клещев, А. И. Петров, В. С. Шеин М.: Недра, 1995.

46. Клещев К. А. Плитотектонические модели нефтегазоносных бессейнов России /К. А. Клещев, В. С. Шеин // Геология нефти и газа. 2004. №1.

47. Клещев К. А. Современное состояние геодинамических основ прогноза, поисков и разведки нефти и газа. /К. А. Клещев, В. С. Шеин // Геология нефти и газа. — 2002. №4.

48. Клещев К. А., Петров А. И., Шеин В. С. Геодинамика и новые типы природных резервуаров нефти и газа. М.: Недра, 1995.

49. Кобруиов А. И. Заметки к истории развития методов решения обратной задачи гравиразведки в XX веке / А.И. Кобруиов // Развитие гравиметрии и магнитометрии в XX веке: Труды конференции. Москва, 23-25 сентября 1996 г. -М.: ОИФЗ РАН, 1997. -С. 188-201.

50. Кобрунов А. И. К анализу линейных приближений обратной задачи структурной гравиметрии//Докл. АН УССР,Б.-1982.- N9.-C.7-9.

51. Кобрунов А. И. К теории интерпретации данных гравиметрии для слоистых сред (равномерная оптимизация)//Изв. АН СССР. Физика Земли,- 1988,- N8.-C. 33-34.

52. Кобрунов А. И. К теории методов подбора//Геофизический журнал.- 1983.-Т.5.-N4.-С.34-43.

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

54. Кобрунов А. И. Метод расщепления решений//Докл. АН УССР, B.-1981.-N7.-С. 19-22.

55. Кобрунов А. И. О введении ограничений типа неравенств на значения плотности при интерпретации гравиметрических данных // Изв.вузов. Геология и разведка.-1981.-N 12.-С.75-81.

56. Кобрунов А. И. О методе оптимизации при решении обратной задачи гравиразведки// Изв. АН СССР. Физика Земли,- 1978,- N8.-C.73-78.

57. Кобрунов А. И. О методе поиска оптимальных решений обратной задачи гравиметрии. Дисс. На соискание ученой степени кандидата наук. Киев, 1978 г. 156 стр.

58. Кобрунов А. И. О построении решений обратной задачи гравиразведки в классе распределений плотности//Докл. АН УССР,B.-1977.-N12.-C. 1077-1080.

59. Кобрунов А. И. О построении решений обратной задачи гравиразведки в классе распределений плотности//Разведка и разработка нефтяных и газовых месторождений,-1978.-Вьтп.15.-С.48-50.

60. Кобрунов А. И. Об одном подходе к решению обратной задачи гравиметрии в плотпостных грапицах//Разведка и разработка нефтяных и газовых месторождений.-1979.-Вып.16.-С.29-33.

61. Кобрунов А. И. Равномерные критерии оптимальности в задачах гравимет-рии//Докл. АН УССР, B.-1986.-N11.-С. 11-14.

62. Кобрунов А. И. Разрешимость и эквивалентность в обратной задаче гравиразведки для нескольких плотностных границ//Изв. АН СССР. Физика Земли,- 1983.-N5.-C. 67-75.

63. Кобрунов А. И. Спектральные представления общего вида решения обратной задачи структурной гравиметрии//Докл. АН УССР, Б.-1988.- N9.-С. 18-21.

64. Кобрунов А. И. Теоретические основы критериального подхода к анализу геофизических данных (на примере задач гравиметрии). Ивано-Франковск, 1985,-229с. Деп в УкрНИИНТИ, 1280 Ук 86.

65. Кобрунов А. И. Теоретические основы решения обратных задач геофизики: Учебное пособие.- Ухта: УИИ, 1995.-226 с.

66. Кобрунов А. И. Теория и методы автоматизированной интерпретации гравиметрических данных для сложнопостроенпых сред / А.И. Кобрунов // Разведочная геофизика: Обзор/МГП 'Теоинформарк". -М., 1993. 51 с.

67. Кобрунов А. И. Экстремальные классы в задачах гравиметрии и их использование для построения плотпостных моделей геологических сред. Дисс. На соискание уч. Степени доктора физ. мат наук. Ивано-Франковск, 1983 г. С. 439.

68. Кобрунов А. И., Аникеев С. А., Варфоломеев В. А., Деиисюк Р. П., Степашок В.П. Комплекс программ решения обратной задачи гравиразведки в классе распределений масс в профильном варианте "Масса-2"// Гос.ФАП СССР N 5087000003 от 5.03.86.

69. Кобрунов А. И., Аникеев С. А., Дениеюк Р. П., Сгепанюк В. П. Комплекс программ решения обратной задачи гравиразведки в классе плотностных границ в профильном варианте "Граница-2" //Гос.ФАП СССР N 5087000035 от 5.03.86.

70. Кобрунов А. И., Войпова О. В. Результаты исследований оптимальных в равномерной метрике решений обратной задачи гравиметрии // Разведка и разработка нефтяных и газовых месторождений. -Львов, 1989.- Вып.26.-С.29-32.

71. Кобрунов А. И., Дениеюк Р. П. Решение обратной задачи гравиразведки в классе плотностных границ с переменной плотностью па контакте // Изв. вузов. Геология и разведка.-1982.-N9.-C. 108-117.

72. Кобруиов А. И., Панасенко В. Н. О вариационном методе решения обратной задачи гравиразведки// Разведка и разработка нефтяных и газовых месторождений.-1976.-Вып. 13.-С.47-51.

73. Кобрунов А.И. Дениеюк Р. П. О выборе параметра релаксации при решении обратной задачи гравиразведки// Геофизический журнал.-1983.-T.5.-N2.-C.63-68.

74. Кобрунов А.И. Некоторые особенности методов подбора в геофизических задачах// Докл. АН УССР,Б.-1984.- N4.-C. 10-13.

75. Кобрунов А. И. О классах оптимальности решения обратной задачи гравиразведки //Изв. АН СССР. Физика Земли.- 1982,- N2.-C. 100-107.

76. Кобрунов А. И. Оценки эквивалентности в методах подбора// Разведка и разработка нефтяных и газовых месторождений.-1983.- Вып. 17.-С.51-54.

77. Коваль Л. А., Гольдшмидт В. И. Автоматизированная система обработки данных магниторазведки с применением ЭВМ (АСОМАМ). Алма-Ата: Изд-во КазВИРГ, 1973.

78. Ковтун А. И., Порохова Л. Н., Гольцман Ф. М. Оптимальная машинная интерпретация кривых МТЗ // Изв. АН СССР. Физика Земли. 1969. №4.

79. Кольская сверхглубокая. Научные результаты и опыт исследований: Ред. В. П. орлов, Н. П. Лаверов. -М.: Технонефтегаз, 1998. 260 с.

80. Костюченко С. JI. Структура и тектоническая модель земной коры Тимапо-Печорского бассейна по результатам комплексного геолого-геофизического изучения // Тектоника и магматизм Восточно-Европейской платформы. М., КМК, 1994. С. 121 -133.

81. Кузнецов О. JL, Никитин А. А. Гсоинформатика. М.: Недра, 1992. 302 с.

82. Кулешов В. Е. Интегрированные технологии интерпретации и возможности их взаимодействия // Международная молодежная научная конференция «Севергеоэко-тех — 2006» (17-19 марта 2006 г.; Ухта): Материалы конференции. Ухта: Изд-во УГ-ТУ, 2006 г. С. 99-104.

83. Кулешов В. Е., Кобрунов А. И., Шилова С. В. Повышение достоверности построения геолого-геофизических моделей па поисково-разведочном этапе работ на нефть и газ // Нефтяное хозяйство. 2007. -№10. С. 64 -67.

84. Лаврентьев М. М. О задаче Коши для уравнения Лапласа // Докл. АН СССР. 1955. Т. 102, №2.

85. Леонов А. С. Об устойчивом решении обратной задачи гравиметрии на классе выпуклых тел. Изв. АН СССР, сер. "Физика Земли", 1976, № 7.- С.55-65.

86. Линник Ю. В. Способ наименьших квадратов и основы теории обработки наблюдений. М.: Физматгиз, 1958. 255 с.

87. Логачев А. А., Захаров В. П. Магниторазведка. М.: Недра, 1973. 347 с.

88. Ломтадзе В. В. Интерпретация гравитационных аномалий с помощью цифровых вычислительных машин. -/ Вопр. развед. геофизики, 1967, вып.6.- С.61-65.

89. Ломтадзе В. В. Принципы организации автоматизированной системы обработки геофизических данных // Прикл. геофизика. 1973. Вып. 72.

90. Ломтадзе В. В. Программное обеспечение обработки геофизических данных. Л.: Недра, 1982, 280 с.

91. Луговенко В. Н. Статистический анализ аномального магнитного поля. М.: Наука, 1974. 200 с.

92. Максимов С. П., Шсип В. С. Тектоника литосферных плит — теоретическая основа научного прогресса в геологии нефти и газа // Геология нефти и газа. 1986. -№9. - С. 2-14.

93. Малышев II. А. Разломы Европейского Северо-Востока СССР в связи с нефте-газоносностыо. Л.: Наука, 1986. - 110 с.

94. Марченко В. В. Человеко-машинные методы геологического прогнозирования. М.: Недра, 1988. 232 с.

95. Миков Д. С. Атлас теоретических кривых для интерпретации магнитных и гравитационных аномалий. Томск: Том. политехи, ин-т, 1955.

96. Миков Д. С. Атлас теоретических кривых для интерпретации магнитных и гравитационных аномалий. М., Госгеолтехиздат, 1956.

97. Минц М. В., Глазнев В. Н., Конилов и др. Ранний докембрий северо-востока балтийского щита: Палеогеодинамика, строение и эволюция континентальной коры. -М: Научный мир, 1996. 287 с.

98. Моделирование в рудной магниторазведке (сб. алгоритм и программ): Метод, рук-во. М.: СНИИГГИМС, 1973.

99. Мотрюк Е. II. Развитие теории и методов объёмной реконструкции плотностных моделей сложнопостроенных геологических сред. /Е. Н. Мотрюк// Дис.канд. техн. наук. Ухта, 2004 г. — 150 с.

100. Мудрецова Е. А. Определение элементов залегания крутопадающих пластов по кривой аномалии силы тяжести. Разведка и охрана недр, №5. М. Госгеолтехиздат, 1956.

101. Нефтегазоносность и геолого-геофизическая изученность ТПП: история, современность, перспективы: Моногр./А. И. Дьяконов, Г. В. Рассохин, Н. Д. Цхадая и др. Коми региональное отд-нис РАЕН, 1999.

102. Николаевский В. Н. Геомеханика и флюидодинамика. ~ М.: Недра, 1996. 447 с.

103. Нумеров Б. В. Интерпретация гравитационных наблюдений в случае одной контактной поверхности// Докл.АН СССР, 1930.- А, № 21.- С.566-574.

104. Овруцкий И. Г. Применение методов минимизации негладких функционалов для решения обратных задач геофизики/ Автореф.дисс. на соиск.уч.ст.канд.физ.-мат. наук. К.: ИГ АН УССР, 1983. - 20с.

105. Оганесян С. М. Решение обратных задач гравиметрии интеграционными методами// Изв.АН АрмССР, Науки о Земле, 1981.- № 5. С.68-74.

106. Отчёт о производстве опытной групповой геологической съёмки М 1 : 50 ООО. на среднем Тимане в 1972 1975 г. Текст. / ТГФ; исполн.: В. М. Пачуковский, В. И. Граф, А. М. Плякин. - Ухта, 1978.

107. Петров А. И., Шеип В. С. О необходимости учета современной геодинамики при оценке и пересчете промышленных запасов нефти и газа // Геология нефти и газа. -2001. -№3.

108. Пудовкин И. М. Применение проблемы Неймана к решению некоторых вопросов прикладной магнитометрии. Гидрометиздат, 1950. (Тр. НИИЗМ; Вып. 5 (15)).

109. Пучков В. Н. Образование Урало-Новоземельского складчатого пояса результат неравномерной косоориентированной коллизии континентов // Геотектоника. 1996. №5. с. 66 - 75.

110. Разработка способов сейсмической разведки и интерпретации геофизических полей в сложных сейсмологических условиях Крыма: Отчет о НИР / Фонды п/о «Крымморгеология», НИС ИФИНГ. Кобрунов А. И., Панасенко В. Н. и др. 1975. -105 с.

111. Родионов Д. А. Статистические методы разграничения геологических объектов по комплексу признаков. М.: Недра, 1968, 158 е.

112. Рудаков А. Г. Экспериментальные исследования групп источников на основе частотной теории фильтрации скоростей // Вопросы динамической теории распространения сейсмических волн. Л.: Изд-во ЛГУ, 1962. Вып. 4.

113. Седиментационные бассейны подвижных поясов /В. Н. Данилов, Ю. Б. Силантьев, Л. 3. Аминов и др. М.: Изд-во Академии горных наук, 2001. - 272 е.: ил.

114. Соколов Б. А. и др. Флюидодинамическая модель нефтегазообразования. -М.: изд-во МГУ, 1998.-115с.

115. Соколов Б. А. Эволюционно-динамические критерии оценки нефгегазоносно-сти недр. М.: Недра, 1984.

116. Сорокин JI. В. Гравиметрия и гравиметрическая разведка. М.: Госгоптехиздат, 1953,-484 с.

117. Сорохтин О. Г. Геодинамика, геофизика океана. — М.: Наука, 1979.

118. Сорохтин О. Г., Ушаков С. А. Глобальная эволюция Земли. М.: изд-во МГУ, 1991.-224 с.

119. Старостенко В. И., Заворотько А. Н. Методика и комплекс программ решения обратной линейной задачи гравиметрии на ЭВМ "Минск-22". К.: Наукова думка, 1976. - 64с.

120. Старостенко В. И. Устойчивые численные методы в задачах гравиметрии. К.: Наукова думка, 1978. - 228с.

121. Старостенко В. И., Дядгора В. А., Заворотько А. Н. Об интерпретации гравитационного поля методом подбора // Изв. АН СССР. Физика Земли. 1975. № 4. С.78-85.

122. Статистическая интерпретация геофизических данных. Л.: Изд-во ЛГУ, 1981, 253 с.

123. Статистические методы интерпретации геофизических наблюдений (коллективная монография). Л.: Изд-во Лепингр. ун-та, 1971. 75 с.

124. Страхов В. Н. Будущее теории интерпретации гравитационных и магнитных аномалий // Комплексные исследования по физике Земли. М.: Наука, 1989.

125. Страхов В. Н. Некоторое применение функционально-аналитических методов в математической теории интерпретации гравитационных и магнитных аномалий: Ав-тореф. дис. . д-ра физ.-мат. Наук. М.: ИФЗ АН СССР, 1972.

126. Страхов В. Н. О задачах, решаемых в рамках второй парадигмы в теории интерпретации гравитационных и магнитных полей // Изв. АН СССР. Физика Земли. 1986. №3.

127. Страхов В. Н. О новом этапе в развитии теории интерпретации гравитационных и магнитных аномалий // Изв. АН СССР. Физика Земли. 1977. №12.

128. Страхов В. Н. О проблеме параметризации в обратной задаче гравиметрии// Изв. АН СССР. Физика Земли. 1978. - №6. -С.39-50.

129. Страхов В. Н. О состоянии и задачах математической теории интерпретации магнитных и гравитационных аномалий // Изв. АН СССР. Физика Земли. 1970. №5.

130. Страхов В. Н. О состоянии и проблемах геологической интерпретации данных гравитационных и магнитных наблюдений // Разведочная геофизика СССР па рубеже 70-х годов. М.: Недра, 1974. С. 113-121.

131. Страхов В. Н. От вычислительной геофизики к геофизической кибернетике. //Изв. Вузов Геология и разведка, 1977,№5, С. 110 119.

132. Страхов В. Н., Лапина М. И. О монтажном принципе построения решения обратной задачи гравиметрии/ Геофиз. сб. АН УССР, вып.74, 1976.- С.3-19.

133. Страхов В. Н., Лапина М. И. Прямая и обратная задачи гравиметрии и магнитометрии для произвольных однородных многогранников. В кн.: Теория и практика интерпретации гравитационных и магнитных полей в СССР. - К.: Наукова думка, 1983.-C.3-87.

134. Страхов В. II. О задачах, решаемых в рамках второй парадигмы в теории интерпретации гравитационных и магнитных аномалий// Докл. АН ССР. Физика Земли. 1987,- № 3.- С.56-68.

135. Страхов В. Н. Об общих решениях обратной задачи гравиметрии и магнитометрии. // Изв. Вузов. Геология и разведка, 1978, №4, С. 104 117.

136. Судо М. М. Современная геология. М.: Знание, 1981.-62 с.

137. Тарбаев Б. И. Тектоника гряды Чернышева // Геология и нефтегазоноспость северо-востока Европейской части СССР. Вып. IV. Сыктывкар: Коми книж. из-во, 1977. С. 14-15.

138. Тафеев 10. П. Палетки для определения элементов залегания вертикального пласта по магнитной аномалии. Тр. ВИРГ, вып. 2. М., Госгеолиздат, 1950.

139. Тимонин Н. И. Тектоника гряды Чернышева. Л.: Наука, 1975.

140. Тихонов А. Н. О решении некорректно поставленных задач и методе регуляризации //Докл. АН СССР. 1963. Т. 151, № 3.

141. Тихонов А. II. Об устойчивости обратных задач // Докл. АН СССР. 1943. Т.39, №5.

142. Троян В. Н. Статистические методы обработки сейсмической информации при исследовании сложных сред. М.: Недра, 1982. 182 с.

143. Хаин В. Е. Глобальные закономерности пефтегазонакоплепия: современное состояние проблемы // Глобальные тектонические закономерности нефтегазонакопле-ния. М„ 1985.-С. 5-14.

144. Хаин В. Е. Общая геотектоника. М.: Недра, 1973. - 574 с.

145. Хаин В. Е. Тектоника континентов и океанов. М.: Научный мир, 2001.

146. Химмельблау Д. Прикладное нелинейное программирование. М., Мир, 1975.

147. Шалаев С. В. Геологическое истолкование геофизических аномалий. Л.: Недра, 1972.

148. Шалаев С.В. Геологическое истолкование геофизических аномалий с помощью линейного программирования. Л.: Недра, 1972. - 142с.

149. Шеин B.C. Геодинамический анализ нефтегазоносных территорий и акваторий в связи с поисками месторождений нефти и газа. // Геология нефти и газа. 2007. №2.

150. Шеин В. С: Геология и нефтегазоносность России. М.: Изд-во ВНИГНИ, 2006.

151. Шеин В. С. Условия нефтегазонакопления на платформах и в складчатых областях (с позиции теорий тектоники плит) / В. С. Шеин, К. А. Клещев // Геология нефти и газа. 1984. - №3.

152. Шеин В. С., Клещев К. А. Новые теоретические предпосылки оценки перспектив нефтегазоносностыо. М., 1984. 58 с.: (Обзор / ВНИИИОЭНГ. Сер. «Нефтегазовая геология и геофизика»).

153. Юдин В. В., Дедеев В. А. Геодинамическая модель Печорской плиты. Сыктывкар, 1987. 12 с.

154. Юдин В.В. Орогенез севера Урала и Пай-Хоя. Екатеринбург: УИФ «Наука», 1994.

155. Юньков А. А., Булах Е. Г. О точности определения аномальных масс методом сеток. ДАН УССР, № 11, 1958 б.

156. Юньков А. А., Булах Е. Г. Возможности использования метода сеток для интерпретации аномалий горизонтального градиента силы тяжести. Тр. Ип-та геологических наук АН УССР, сер. геофиз., вып. 2, 1958г.

157. АН Т. Ismail-Zadeh, Sergey L. Kostyuchcnko and Boris M. Naimark The Timan-Pechora Basin (northeastern European Russia): tectonic subsidence analysis and a model offormation mechanism / №283 1997), Tectonophysics, pg.205-218.

158. Bott M. H. P. Inverse methods in the interpretation of magnetic and gravity anomalies // Geophys. Prospect. 1973. Vol.21, N 3.

159. Grant F. S. Review of data processing interpretation methods in gravity and magnetics 1964-1971 //Geophysics, 1972. Vol. 37, N 4.

160. Jin D. V. Two dimentional mass distributions from gravity anomalies: a computer method // Geophys. Prosp. 1974. Vol. 22, N 4.