Методика автоматизированной комплексной интерпретации гравиметрических данных для сложных плотностных моделей: на примере Тимано-Печорской нефтегазоносной провинции тема диссертации и автореферата по ВАК РФ 04.00.12, кандидат геолого-минералогических наук Моисеенкова, Светлана Владиславовна

• Моисеенкова, Светлана Владиславовна
• кандидат геолого-минералогических наук
• 1998

Оглавление диссертации кандидат геолого-минералогических наук Моисеенкова, Светлана Владиславовна

ОГЛАВЛЕНИЕ

ВВЕДЕНИЕ

ГЛАВА 1. Интерпретационные возможности гравиметрии и их эффективность в ТПП

1.1 Классификация методов интерпретации данных гравиразведки

1.2 Анализ эффективности методов интерпретации гравиметрических данных в Тимано-Печорской провинции

ГЛАВА 2. Теоретические основы методов детального описания структурно-плотностного разреза с использованием данных гравиметрии

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

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

2.3 Постановка и решение обратной задачи гравиразведки для нерегулярных сетей наблюдения

ГЛАВА 3. Организация вычислительного процесса

3.1 Алгоритмическая вычислительная схема

3.2 Характеристика программного обеспечения

ГЛАВА 4. Формирование представлений о типичных геолого-геофизических моделях строения осадочного чехла ТПП

4.1 Грабенообразные прогибы. Рифогенные образования. Тектоническое строение

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

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

4.4 Зона клиновидных вдвигов. Тектоническое строение. Проявление в сейсмическом и гравитационном полях

ГЛАВА 5. Методика автоматизированной интерпретации

5.1 Проведение вычислительных экспериментов

5.2 Интерпретация результатов

5.3 Методические рекомендации к автоматизированной интерпретации сейсмогравиметрических данных для условий Тимано-

Печорской провинции

ЗАКЛЮЧНИЕ

ЛИТЕРАТУРА

Введение диссертации (часть автореферата) на тему «Методика автоматизированной комплексной интерпретации гравиметрических данных для сложных плотностных моделей: на примере Тимано-Печорской нефтегазоносной провинции»

ВВЕДЕНИЕ

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

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

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

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

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

При этом в работе решаются следующие задачи:

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

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

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

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

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

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

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

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

Автор выносит на защиту:

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

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

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

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

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

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

Апробация работы. Основные результаты диссертационной работы докладывались на Всесоюзных семинарах им. Д.Г.Успенского "Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей" (г. Воронеж, 1996г.; г.Москва, 1997г.; г.Ухта, 1998г.), на Тринадцатой Коми республиканской молодежной научной конференции (г.Сыктывкар, 1997г.), на научных семинарах кафедры геофизики Ухтинского индустриального института.

Публикации. По теме диссертации опубликовано 14 печатных работ.

Структура и объем работы. Диссертация состоит из введения, _5_глав, заключения, содержит _108_ страниц текста, _60_ рисунков и список литературы из 181_наименования.

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

Глава 1. Интерпретационные возможности гравиметрии их эффективность в Тимано-Печорской провинции

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

1.1. Классификация методов интерпретации данных гравиразведки

Известно, что основной задачей интерпретации геофизических данных принято считать их геологическое истолкование. Непосредственно геологическое истолкование гравиметрических данных предполагает решение двух вопросов: построение плотностной модели изучаемой геологической среды и последующая геологическая интерпретация полученной плотностной модели. Различные методы интерпретации геофизических (и, в частности, гравиметрических) данных, призванные решать такого рода задачи, в вычислительной основе своей содержат решение соответствующих обратных задач. Решение обратной задачи гравиразведки является наиболее важным этапом в процессе построения плотностной модели среды. Для постановки обратной задачи необходимо записать зависимость гравитационного эффекта от вызывающего его распределения плотностных неоднородностей в земной коре. Если выбрать прямоугольную систему координат, расположив плоскость ХОУ параллельно дневной поверхности и направить ось Ъ вниз, то распределение плотности, представленной в виде функции координат с(х,у,г) в некоторой области V нижнего полупространства, связано с вертикальной производной гравитационного потенциала Щхо, уо, г0) в точке (х0,уо, г0) соотношением:

(1.1)

где у - гравитационная постоянная.

Выражение (1.1) представляет собой оператор прямой задачи гравираз-ведки или нахождение элементов гравитационного поля по заданному распределению источников. Под решением обратной задачи понимают нахождение распределения источников при введенных о них модельных представлениях по заданным наблюдаемой и оператору прямой задачи. Известно, что такая задача является некорректно поставленной [164]. Всякую задачу называют корректной (по Адамару), если ее решение: 1) существует; 2) единственно; 3) устойчиво. Задача называется некорректной, если не выполняется одно из этих условий. Решение вопросов о корректности обратных задач занимает в теории интерпретации геофизических данных центральное место. Так, трудами А.Н.Тихонова, М.М.Лаврентьева, В.К.Иванова, В.Н.Страхова и их последователей развита теория некорректных задач [18, 81, 28, 117, 121, 123, 159]. В основе этой теории лежит идея привлечения в постановку обратной задачи дополнительной информации об искомом решении. Непосредственно поиск решения обратной задач осуществляется с помощью регуляризующих алгоритмов [135, 150, 164 и др.]. Рассмотрению вопросов о проблеме единственности обратной задачи гравиразведки посвящены работы П.С.Новикова, А.Х.Остромогильского, В.Н.Страхова, В.И.Старостенко, В.Г.Филатова, В.Г.Чередниченко и др. [104, 126, 139, 166, 169, 171 и др.]. Результаты исследований в этой области породили так называемые "теоремы единственности". Доказательство этих теорем возможно при введении некоторых ограничений или упрощений, позволяющих выделять такие типы наборов элементарных тел (классы единственности), параметры которых однозначно определяются по гравитационному полю. Проблема неединственности обратной задачи сопровождается ее неустойчивостью: малые изменения во входных данных порождают различные результаты решения [164]. А так как реальные гравиметрические данные всегда осложнены ошибками, то решение обратной задачи гравиразведки будет неустойчивым. Снижение этого эффекта позволили осуществить методы регуляризующих процедур.

На сегодняшний день множество утвердившихся и разрабатываемых методов решения обратных задач гравиметрии применяются в различных геологических ситуациях. Если рассматривать историю развития теории и методов интерпретации гравиметрических (и в целом геофизических) данных, то можно выделить четыре типа задач, каждый из которых последовательно становился приоритетным на соответствующем этапе интерпретационного обеспечения гравиметрии [138, 144, 146, 153, 155]:

1. Обнаружение геологических объектов.

2. Пространственная локализация геологических объектов.

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

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

Основоположником теоретических методов интерпретации гравитационных аномалий по праву можно считать известного русского ученого прошлого века Ф.А.Слудского. В своих работах еще в 1863 году он сформулировал идею моделирования, описал метод подбора, в 1893 - предложил приемы определения глубины центра тяжести по характерным точкам поля. Уровень работ Ф.А.Слудского был достигнут спустя несколько десятилетий - в 30-е годы нашего столетия - в период становления теории интерпретации потенциальных полей. На раннем этапе развития гравиметрии наиболее характерной задачей являлась задача обнаружения геологических объектов, состоящая в принятии решения о наличии либо отсутствии изучаемых объектов. Величина и форма гравитационной аномалии подвергались качественному описательному анализу, на основе чего делались соответствующие выводы. Так, впервые, в 20-х годах А.Д.Архангельский для изучения строения кристаллического фундамента Русской платформы применял гравиметрические карты в сопоставлении их с геологическими [4]. Подобный подход наблюдается и в современных методах интерпретации гравиметрических данных, например, в процедурах анализа трансформаций гравитационного поля. Анализу с целью обнаружения объектов подвергается не само гравитационное поле, а его преобразование - трансформация, полученная с помощью тех либо

и

иных вычислительных процедур. В настоящее время большое количество разработанных методов трансформаций гравитационных полей позволяет осуществлять различные цели: подавление (фильтрацию) случайных погрешностей [101]; усиление локализации аномальной компоненты поля; исключение регионального фона [2]; построение вертикального разреза аномального поля и пр. Это методы аналитического продолжения аномалий [137] и их пересчета в производные различных порядков [86-88], корреляционно-статистические [19, 20, 24, 173] и методы условных трансформаций [1, 11]. Результаты некоторых трансформаций поля могут рассматриваться как формальные решения обратной задачи гравиразведки (например, аналитическое продолжение поля в нижнее полупространство) или, во всяком случае, являются средством извлечения дополнительной информации, уточняющей в последствии решение обратной задачи.

Следующий этап развития теории интерпретации характеризуется преобладанием задач пространственной локализации геологического объекта, состоящих в первую очередь в оценке пространственного расположения, а также поиске его характерных параметров, таких как центра тяжести, суммарной массы, размеров объекта, элементов залегания. Первые работы по решению подобных задач появляются, приблизительно, начиная с 1927 г. Это, например, работы по изучению конфигурации одной контактной поверхности Б.В.Нумерова [109]. Исследования по определению характеристик возмущающих тел различными методами проводились так же Д.С.Миковым, Г.А.Гамбурцевым и др.[17, 97].

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

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

симости от результатов сопоставления. В рамках этого направления созданы эффективные приемы решения прямых задач на основе как палеток (Г.А.Гамбурцев, К.Ф.Тяпкин, В.Н.Страхов и др.) так и атласов рассчитанных теоретических кривых [96 и др.]. Получившие на этом этапе развитие так называемые «ручные» методы долгое время являлись наиболее распространенными и в настоящее время они так же применяются уже в более совершенном виде, с использованием современной вычислительной техники.

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

А(с*х)=и(х), (1.2)

где с1={с11, сЬ,...с1п} - параметры тела произвольной формы, определяющие пространственные и плотностные его характеристики и аналитически связаны с интерпретируемой компонентой и(х) гравитационного поля. С целью решения уравнения (1.2) созданы эффективные методы характерных точек гравитационного поля, графические методы [98]. Разработанные в этот период методы определения интегральных характеристик возмущающих тел так же позволяют определить основные элементы трехмерных и двухмерных тел, например, гармонические моменты масс [26], определение избыточной массы тела, вертикальной координаты центра тяжести и т.д. [134]. В силу упрощенных моделей изучамых сред вопросы единственности решения обратной задачи в то время не возникали, хотя уже были известны работы П.С.Новикова по этому вопросу [104]. Позже, рассмотрение и построение эффективных схем решения задачи об одной контактной поверхности привело к необходимости изучения вопросов единственности результатов решения обратной задачи. Так, в своих работах [141 и др.] В.Н.Страхов показал, что обратная задача для контактной поверхности обладает двухпараметрической неединственностью, в отличие от однопараметрической для задачи нахождения плотности изолированных объектов. Еще большая неединственность в решении

обратной задачи возникает при переходе к нескольким контактным поверхностям. В этом случае исследователи столкнулись с проблемой выделения из суммарного поля гравитационных эффектов от каждой поверхности. По трудности решения эта задача оказалась близкой с исходной обратной задачей. Исследования в этой области отражены в работах В.Н.Страхова, И.А.Непомнящих, В.Г.Филатова, Ю.В.Антонова и др. [2, 151, 167 и др.].

Результаты в области разработки методов решения прямых задач, теоретические исследования по эквивалентности в обратных задачах теории потенциала, формирование идей теории регуляризации некорректно поставленных задач предопределили появление формализованных методов решения обратных задач с целью расчленения геологического разреза. Начало этому направлению было положено работами Булаха Е.Г., Юнькова A.A., Шалаева C.B. [15, 172, 174, 175], а разработанные в этом направлении различные подходы к решению обраных задач получили название "методы подбора". Исследования, посвященные методам подбора, отражены так же в трудах В.Н.Страхова, В.И.Старостенко, С.М.Оганесяна, К.Ф.Тяпкина, Г.Я.Голиздры, А.И.Кобрунова и многих других [16, 40, 42, 135, 165 и др.]. В рамках этих методов рассматриваются три главных компонента: модельный класс, критерий подбора и алгоритм минимизации.

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

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

3. Алгоритм решения обратной задачи - представляет собой способ осуществления минимизации функционала J(uH, A(m)), в результате чего определяется искомая плотностная модель т\

В рамках методов подбора выделяется группа технологий, применяющих в качестве модельного класса К конструкции призматического типа [5, 8, 15, 19, 82, 84, 135, 136, 140, 161, 162]. В качестве минимизируемого функционала применяется квадратичная метрика. При необходимости вводятся веса - например, при неравноточности задания гравитационного поля, и высшие производные гравитационного потенциала - при подборе верхней части разреза. В зависимости от того, как задано гравитационное поле и с какой точностью, а так же от таких факторов, как наличие регионального фона, влияния пород, залегающих вне области поиска решения и т.п. принимается различный вид функционала J(uH, A(m)).

Следующая группа технологий в рамках методов подбора состоит в нахождении некоторых формальных источников наблюдаемого гравитационного поля с последующим их эквивалентным перераспределением. Такие источники находятся либо в результате осуществляемой аналитической аппроксимации поля с последующим пересчетом параметров подобранной аппроксимации в параметры источников, либо в результате непосредственно решения обратной задачи как задачи минимизации функционала J(uH, A(m)) с использованием алгоритмов регуляризации [27, 83, 99, 100, 102, 103, 110, 142, 147, 160, 168, 170 и др.].

В следующем направлении в рамках методов подбора использовалась полиномиальная аппроксимация [6,7] для нахождения плотности, а описание конфигурации плотностной границы производилось отрезком ряда Фурье [9, 10]. В качестве модельного класса к вводилась латерально изменчивая плотность [93-95, 105-108].

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

кают проблемы устойчивости и единственности. Как это показал в своих работах А.И.Кобрунов, при решении обратных задач такого типа возникает "эффект скрытой эквивалентности" [40, 42, 53, 57, 58, 67], который состоит в следующем. На исходном широком классе может найтись некоторый бессмысленный с геологической точки зрения элемент сГ(у) (например, гармоническая функция): А(о(у))=А(<тл(у)), который может быть аппроксимирован некоторым элементом пг из класса к лучше, чем те аппроксимирует а(у). Наиболее этот эффект проявляется в сеточных методах при решении линейных обратных задач гравиразведки, в которых искомыми являются значения плотности в узлах сетки, а сами сетки имеют внутренние ячейки [57, 62]. Актуальным становится вопрос о том, какое решение из области эквивалентности получается при решении обратной задачи методом подбора. Таким образом, из-за проявления эффектов эквивалентности следует внимательнее относиться к результатам методов подбора, особенно если таковой производился для сложнопостроенных сред.

Необходимо отметить развитие отдельного направления в теории интерпретации гравиметрических данных - статистических методов интерпретации, развитых в работах Г.И.Каратаева, В.И.Шрайбмана, М.С.Жданова, Ф.М.Гольцмана, Т.Б.Калининой и др. [20, 24, 29, 173 и др.]. Общая схема интерпретации геофизических аномалий корреляционно-статистическими методами следующая. Для расчетов выбираются эталонный и исследуемый участки, близкие по геологическому строению. По данным различных геолого-геофизических исследований на эталонном участке устанавливается корреляционная связь между глубиной плотностной границы и какой-либо компонентой гравитационного поля. Составленное уравнение регрессии применяется к исследуемому участку. То есть, результаты решений во многом зависят от степени соответствия эталонной и исследуемой площадей, а так же от полученной тесноты связи на эталонной площади. В статистических методах проблемным вопросом является выделение аномального гравитационного поля, наилучшим образом связанного с конкретной плотностной границей.

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

Развитие методов интерпретации гравитационных аномалий на рассматриваемом этапе с привлечением идей и методов вычислительной математики привело к возникновению и становлению такого направления как вычислительная геофизика [155]. А появление и внедрение вычислительной техники определило разработку и массовое применение в комплексе геологоразведочных работ различных систем автоматизированной интерпретации гравиметрических данных (В.В.Ломтадзе, В.И.Старостенко, Е.Г.Булах, [84, 136 и др.]).

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

В середине семидесятых - начале восьмидесятых годов появляются работы А.И.Кобрунова, развивающие методы решения обратных задач гравиметрии в рамках так называемого критериального подхода [47, 49, 50, 76]. Получение решения обратной задачи - единственного из множества эквива-

лентных - производится с помощью минимизации некоторого функционала (критерия оптимальности J), сконструированного на основе дополнительной априорной информации о параметрах среды. Основные положения данного подхода могут быть отнесены в операторной форме к решению обратных задач других геофизических методов, что позволяет строить удобные вычислительные схемы комплексной интерпретации для различных геофизических методов [23, 25, 78, 127].

Вопросы, связанные с такой постановкой обратной задачи, а так же их решение широко отражены в трудах В.Н.Страхова [143, 145, 148, 149, 152, 156, 157 и др.], А.И.Кобрунова [30-33, 35-38, 40, 42, 43, 46-56, 59-66, 69-78],

A.С.Маргулиса, М.В.Новоселицкого [90-92], С.М.Оганесяна,

B.И.Старостенко [111-124]. Подобная постановка известна в работах Backus G. и Gilbert F., разработавших линейную теорию решения обратных задач при использовании ограниченного набора данных [177-181], которая широко применяется зарубежом для решения различных обратных задач гравиметрии, сейсмологии и др. Подробная характеристика этих методов дана в [176].

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

1.2. Анализ эффективности методов интерпретации гравиметрических данных в Тимано-Печорской провинции

Первые геофизические исследования в Коми крае начались в 1931 г. -тогда под руководством Л.П.Смирнова и Л.В.Сорокина была проведена ва-риометрическая съемка по маршруту вдоль реки Ухты. В том же году под руководством В.М.Озерецкого опытные работы с вариометром проведены на реке Печоре. Позднее региональной маршрутной гравиметрической съемкой были охвачены территории вдоль тракта Усть-Вымь-Чибью, рек Ижма, Печора (В.И.Арест, В.В.Федынский, 1933 г.), Уса (Г.И.Разбегаев, 1937 г.). Первые материалы, полученные по результатам гравиметрических съемок, позволили выделить основные черты тектонического строения территории Тимано-Печорской провинции.

В течение 1931 - 1947 гг. гравиметрические исследования (вариомет-рическая, а затем и гравиметрическая съемка) в комплексе с магниторазведкой были проведены на территории Юго-Западного Притиманья, Тимана, Печорской гряды, Предуральского прогиба, а по их результатам удалось выявить характер сочленения надпорядковых тектонических элементов.

В 1950 - х годах точность проводимых гравиметрических наблюдений позволила осуществлять картирование масштабами 1:200000 и 1:100000 с целью геотектонического районирования и выявления перспективных структур. Съемкой масштаба 1:500000 была охвачена северная тундровая часть провинции. Благодаря этим исследованиям была построена гравиметрическая карта по всей территории ТПП.

Первое геотектоническое районирование Тимано-Печорской провинции было осуществлено благодаря гравиметрическим исследованиям, проведенным в комплексе с аэромагнитной съемкой - в 1955 г. Тиманской экспедицией ВНИГРИ (О.А.Калинина) составлена "Тектоническая схема Тимано-Печорской газонефтеносной области по данным геофизики", основой которой была гравиметрическая карта.

В шестидесятые - семидесятые годы гравиметрическая съемка масштабов 1:100000 и 1:50000 проводилась в районах перспективных на нефть и газ: Верхне-Печорской впадине, Печоро-Кожвинском мегавале, южной части Усино-Кожвинского мегавала и Хорейверской впадины, Косью-Роговской впадине, северо-восточном склоне Тимана и Средне-Печорском поперечном поднятии. Материалы гравиметрических исследований позволили уточнить границы тектонических элементов по фундаменту и осадочному чехлу, а так же выявлять локальные структуры в осадочном чехле. При интерпретации гравиметрических данных использовались в основном методы трансформаций поля силы тяжести (расчет производных гравитационного потенциала, аналитическое продолжение поля силы тяжести), а так же метод нормированного градиента В.М.Березкина. Так, в 1963 г. впервые была сделана попытка пересчета гравитационного поля в нижнее полупространство (К.А. Кривцов [80]). В 1964-65 гг. (М.А.Осада [125]) с использованием гравиметрических данных построена первая карта поверхности фундамента. В 1967 г. в рамках разведочных работ по Вуктылу (О.И.Плескачев, [129]) вручную (из-за отсутствия необходимой вычислительной техники) было рассчитано по одному профилю распределение нормированного градиента (Gn), так же был произведен пересчет аномального поля на глубины. Полученные результаты не были оценены авторами положительно. В 1970 г. в рамках работ по теме "Обобщение материалов детальных гравиметрических съемок Тимано-Печорской провинции. Составление единой карты гравитационного поля и его производных с применением ЭЦВМ. Составление схемы глубинного строения Тимано-Печорской провинции" (К.А. Кривцов, 1970г. [79]) начали применяться программы на ЭВМ для вычисления трансформированных гравитационного и магнитного полей. Как это отметили авторы, трансформированные аномалии были отягощены ошибками (особенно магнитные). Но несмотря на это, пересчитанные вверх остаточные аномалии позволили провести тектоническое районирование территории, выделить крупные блоки фун-

дамента, а пересчет аномального поля вниз позволил более точно проводить границы между сопряженными блоками фундамента, залегающего на различных гипсометрических уровнях или обладающих различной плотностью. В этой работе с целью решения поставленных задач авторы попытались установить корреляционную связь между глубинами до фундамента и значениями остаточной аномалии. С 1971 г. стали проводиться тематические работы с целью построения тектонической карты провинции (А.С.Бушуев, Ф.Н. Снисарь, Ю.Ф.Кузьмин, В.П.Сергеев, Е.С.Подловилин и др.) по данным комплекса геофизических методов разведки и бурения. В рамках этих работ велись исследования по изучению глубинного строения провинции и поискам локальных структур в осадочном чехле. Таким образом, в шестидесятых годах начинают применяться различные методы интерпретации гравиметрических данных (в том числе и с применением вычислительной техники) - метод полного нормированного градиента В.М.Березкина (О.И.Плескачев, 1967 г., С.И.Максимова, 1968 г.), методы пересчета гравитационных аномалий (К.А.Кривцов, Н.В.Вихорнов, 1963-64 гг., К.А.Кривцов, 1970 г. и др.).

В восьмидесятых годах гравиметрическая съемка уже проводилась в комплексе детальных геофизических исследований с целью: структурно-литологического картирования (масштабы 1:50000, 1:25000, 1:10000) территорий Среднего и Южного Тимана (Э.М.Репин, Э.М.Ходжаев, В.В.Мартынов, И.М.Плякина); изучения особенностей строения осадочного чехла нефтеперспективных районов (например, прослеживание верхнедевонских барьерных рифов на Западно-Тэбукской структуре и к северо-востоку от нее, С.И.Максимова, 1976-1982 гг.); детальное изучение строения фундамента и осадочного чехла, выявление барьерных рифов (масштабом 1:50000, 1:25000).

В период до 1994 г. детальная гравиметрическая съемка масштабом 1:25000 была проведена на территориях Ижма-Печорской впадины, Печоро-Кожвинского мегавала, гряды Чернышева, Хорейверской впадины и Вычегодского прогиба (см. рис. 1.1). Остальные территории остались покрытыми

Рис. 1.1 Схема изученности ТПП методом гравиразведки

гравиметрической съемкой меньшего масштаба - вплоть до 1:200000 (некоторые территории Ижма-Печорской впадины, северные тундровые части провинции). Интерпретация результатов детальных гравиметрических исследований показала высокую информативность последних в районах с неуверенными сейсмическими материалами: в прибортовых частях ВерхнеПечорской и Косью-Роговской впадин, гряде Чернышева. Интерпретация полученных гравиметрических данных, а так же переинтерпретация материалов прошлых лет в период восьмидесятых - девяностых годов представляла собой опробованные ранее методы интерпродолжения аномального гравитационного поля; расчеты вертикальных производных различного порядка и составление карт их распределения; применение метода нормированного градиента [11, 12]. Так, в рамках комплексных работ «Обобщение результатов грави- электроразведочных работ совместно с данными сейсморазведки и бурения на отдельных площадях с целью повышения эффективности геофизических исследований в ТИП» (А.Н.Мальцева, 1988 г.) был применен метод GRAD Е.А.Мудрецовой, основанный на регуляризованном аналитическом продолжении полей в сторону источников аномалий. В результате работ было проведено тектоническое районирование - поверхность фундамента картировалась по положению особых точек и зон максимальных градиентов продолженных аномалий силы тяжести [89]. Другой пример - применение метода нормированного градиента (с расчетом его фазовых характеристик) при интерпретации гравиметрических данных по Когельской площади (Т.И.Гончаренко, [22]) позволило определить положение тела рифа, которое было приурочено к максимумам Gn.

В начале девяностых годов в единичных случаях в практику начинают внедряться новые методы интерпретации гравиметрических данных, основанные на применении критериального подхода. В рамках этих работ (С.И.Максимова, Н.С.Борисов, 1990 г. [13, 14], 1994 г.) были получены качественно новые результаты - а именно, в трехмерной постановке решены обратные задачи гравиразведки в классе распределения плотностей и построения плотностных границ, в результате чего было уточнено строение Ярвож-

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

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

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

Заключение диссертации по теме «Геофизические методы поисков и разведки месторождений полезных ископаемых», Моисеенкова, Светлана Владиславовна

ЗАКЛЮЧЕНИЕ

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

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

В работе решены следующие задачи:

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

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

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

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

В результате проведенных исследований впервые были предложены:

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

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

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

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

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

Список литературы диссертационного исследования кандидат геолого-минералогических наук Моисеенкова, Светлана Владиславовна, 1998 год

Литература

1. Андреев Б.А., Клушин И.Г. Геологическое истолкование гравитационных аномалий,- Л.: Недра, 1965,- 495 с.

2. Антонов Ю.В. Разделение сложных аномальных полей силы тяжести,- Воронеж: ВГУ, 1985,- 214 с.

3. Аронов В.И. Обработка на ЭВМ значений силы тяжести при произвольном рельефе поверхности наблюдений. М.: Недра, 1976,- 131 с.

4. Архангельский А.Д. О соотношениях между геологической структурой и аномалиями силы тяжести Европейской части СССР// Докл. 7-й Балтийской геодез. комиссии. - Л., 1934. - Вып. 6. - С. 3-11.

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

6. Банк П.И. К аналитическому решению трехмерной обратной задачи теории потенциала// Докл. АН СССР, 1973,- т. 213, №4,- С. 830 -832.

7. Балк П.И., Балк Т.В. Трехмерная прямая и обратная задача гравиразведки при полиномиальной аппроксимации плотностных неоднородностей// Геология и геофизика, 1986,- № 5,- С. 120-124.

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

9. Бережная Л.Т., Телепин М.А. Определение плотности по гравиметрическим данным// Прикладная геофизика. - М.: Недра, 1966. -Вып.47- С.130-139.

10.Бережная Л.Т., Телепин М.А. Решение обобщенной обратной задачи гравиметрии для контактной поверхности// Прикладная геофизика.-М.: Недра, 1971. - Вып.64. - С.110-123.

П.Березкин В.М. Применение гравиразведки для поисков месторождений нефти и газа,- М.: Недра, 1973,- 164 с.

12.Березкин В.M. Применение метода полного нормированного градиента при геофизической разведке. М.: Недра, 1988,- 188 с.

13.Борисов Н.С. Отчет Средне-Кочмесской сейсмопартии 50489 о результатах региональных работ масштаба 1:200000 в 1988-1990 гг. ИГЭ, ОМЭ ПГО "Печорагеофизика",- Ухта, 1990,- 230 с.

14.Борисов Н.С. Комплексное обобщение результатов геолого-

геофизических работ с целью изучения геологического строения осадочного

чехла Косью-Роговской впадины// Отчет по теме 9043 за 1988-1990 гг. Книга \ ...

I. Текст отчета. ОМЭ ПГО "Печорагеофизика",- Ухта, 1990,- 290 с.

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

16.Булах Е.Г., Маркова М.Н. Решение обратных задач гравиметрии методом подбора//Геофизич. журнал,- 1992,- 14,- №4,- С.9-19.

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

18.Гласко В.Б. Обратные задачи математической физики,- М.: Изд-во МГУ, 1984,- 112 с.

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

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

21.Гончаренко Т.И. Отчет Шарью-Заостренской гравиразведочной партии 40694 о результатах поисковых работ масштаба 1:25000 в 1993-94 гг./ Книга I. Текст отчета. ОАО "Севергеофизика".- Ухта, 1997,- 60 с.

22.Гончаренко Т.Н. Отчет Когельской партии 41489 о результатах гравиразведочных работ масштаба 1:50000, проведенных в 1989 году на Когельской площади/ Книга I. Текст отчета. ПГО "Печорагеофизика", СГЭ №10. - Ухта, 1990. - 78 с.

23. Даниленко А.Н. Теория и методы построения термо-электроизображений пластов по данным геофизических исследований скважин. Автореф. дисс...канд. геол.- мин. наук,- Ухта, 1996,- 22 с.

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

25.Журавлева О.И. Методы решения обратных задач структурной гравиметрии на основе равномерных критериев оптимальности. Автореф. дисс...канд.физ,- мат. наук,- Киев, 1996,- 25 с.

26.3аморев A.A. Об определении производных гравитационного потенциала и соотношений между моментами возмущающих масс по производной, заданной на плоскости// Изв. АН СССР. Сер. геогр. и геофиз., 1939,- № 3 - С.275-286.

27.Иванов В.К. О разрешимости обратной задачи логарифмического потенциала в конечном виде// Докл. АН СССР, 1956,- Т.106, №4,- С.598-600.

28.Иванов В.К., Васин В.В., Танана В.П. Теория линейных некорректных задач и ее приложения,- М.: Наука, 1978,- 206 с.

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

30.Кобрунов А.И. The Theory Of Interpretation Of Gravity Data For Media Of Intricate Structure// Academic Express. Geophysical Express.- Lviv, 1994,-Vol. 1 (1).- 20 p.

31.Кобрунов А.И. Алгоритм решения систем линейных уравнений большой размерности// Геофизический журнал, 1984.-Т.6, №1,- С.81-85.

32.Кобрунов А.И. Варфоломеев В.А. К построению пространственных распределений масс с нулевым гравитационным эффектом// Геофизический журнал, 1983,- Т.5, №3. - С.53-56.

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

34.Кобрунов А.И. и др. Разработка методики выявления структурно-тектонических ловушек нефти и газа и ее применение в южной части Ижма-Печорской синеклизы. Книга I. Текст отчета. КРО РАЕН,- Ухта, 1997,- 231с.

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

36.Кобрунов А.И. К вопросу об интерпретации аномальных гравитационных полей методом оптимизации// Изв. АН СССР. Физика Земли, 1979,- №10,- С. 67-78.

37.Кобрунов А.И. К вопросу об эквивалентных перераспределениях в классе плотностных границ// Теория и методика интерпретации гравимагнитных полей,- Киев: Наукова Думка, 1981.- С.97-105.

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

39.Кобрунов А.И. К теории комплексной интерпретации// Геофизический журнал, 1980,- Т.2, №2,- С.31-38.

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

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

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

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

44.Кобрунов А.И. О детерминистическом подходе к теории комплексной интерпретации геофизических полей// Геофизические исследования глубинного строения земной коры,- Киев: Наукова Думка, 1979,- С.60-65.

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

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

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

48.Кобрунов А.И. О методе поиска оптимальных решений обратной задачи гравиразведки в классе распределений плотности// Геофизические исследования глубинного строения земной коры,- Киев: Наукова Думка, 1979,- С.65-70.

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

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

51.Кобрунов А.И. Об одной постановке задачи оптимизации, возникающей при интерпретации комплекса геофизических данных// Геофизический журнал.-1982,- Т.4,- №3.- С.50-56.

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

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

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

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

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

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

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

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

60.Кобрунов А.И. Теория интерпретации данных гравиметрии для сложнопостроенных сред: Учебное пособие.- Киев: УМК ВО, 1989.-100 с.

61.Кобрунов А.И. Теория интерпретации данных гравиметрии для сложнопостроенных сред// Геофизический журнал.-1995. - Т. 17. №1,- С.3-12.

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

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

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

65.Кобрунов А.И., Аникеев С.А., Денисюк Р.П., Степанюк В.П. Комплекс программ решения обратной задачи гравиразведки в классе распределения масс в площадном варианте "Масса-3"// Гос. ФАП СССР № 5087000081 от 22.02.87.

66.Кобрунов А.И., Аникеев С.Г., Благий И.И., Денисюк Р.П. Комплекс программ решения обратной задачи гравиразведки в классе плотностных границ в площадном варианте// Гос. ФАП СССР № 50880001415, 1988.

67.Кобрунов А.И., Варфоломеев В.А. О распределении масс с нулевым гравитационным эффектом// Разведка и разработка нефтяных и газовых месторождений.-1981,- Вып. 18,- С.27-33.

68.Кобрунов А.И, Варфоломеев В.А. Об одном методе г-эквивалентных перераспределений и его практическом использовании при интерпретации гравитационных полей// Изв. АН СССР. Физика Земли.-1981,- №10,- С. 25-44.

69.Кобрунов А.И., Варфоломеев В.А., Денисюк Р.П. Комплекс методов для интерпретации данных гравиразведки// Геофизический журнал.-1983,-Т.5.- №5,- С.3-13.

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

71.Кобрунов А.И., Денисюк Р.П. Об одном эвристическом алгоритме решения обратной задачи гравиметрии в классе плотностных границ// Разведка и разработка нефтяных и газовых месторождений,-1981,- Вып. 18,-С.33-39.

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

73.Кобрунов А.И., Журавлева О.И., Использование спектральных представлений для решения обратной задачи гравиразведки структурного типа (равномерная оптимизация)// Изв. АН СССР. Физика Земли.- 1991,-№5,- С. 47-58.

74.Кобрунов А.И., Журавлева О.И. Схема решения нелинейной структурной задачи гравиметрии на основе спектральных представлений (равномерная оптимизация)// Докл. АН УССР, Б.-1990,- №11,- С.18-21.

75.Кобрунов А.И., Ман О. ОЬгаэепа и1оЬа ^аутеШе-зеБ^еш пегоНка 1ш81о1шс11 гохЫаЫ// Межд. конференция "Современные проблемы гравиметрии". Прага (Либерцы), 1984,- С. 7.

Типовые программы обработки геофизических данных на ЦВМ Минск-2,-Л.,ВИРГ, 1966,- 123 с.

86.Маловичко А.К. Вычисление гравитационного аномального эффекта, создаваемого структурами третьего порядка// Вопросы обработки и интерпретации геофизических наблюдений. Вып.З,- Пермь, 1962. - С. 3-10.

87.Маловичко А.К., Костицын В.И. Гравиразведка,- М.: 1992,- 357 с.

88.Маловичко А.К., Костицын В.И., Тарунина О.Л. Детальная гравиразведка на нефть и газ,- М.: Недра, 1989,- 224 с.

89.Мальцева А.Н. Обобщение результатов грави-, электровазведочных работ совместно с данными сейсморазведки и бурения на отдельных площадях с целью повышения эффективности геофизических исследований в ТПП. - Ухта. ПГО "Печорагеофизика", СГЭ №10, 1988. - 193 с.

90.Маргулис A.C. Гармонические плотности и обратные задачи потенциала. - В кн.: Теория и методика интерпретации гравимагнитных полей. К.: Наукова думка, 1981,- С.130-136.

91.Маргулис A.C. К теории потенциала в классах Lp(W)// Изв. ВУЗов, Математика, 1982,- № 1,- С.33-41.

92.Маргулис A.C. О единственности решения обратной задачи гравиметрии для структурных моделей// Докл. АН СССР, 1984,- т.285, № 2,-С.242-246.

93.Маргулис A.C. Рудные и структурные обратные задачи гравиметрии. Нормальные решения и их приложения. Автореф. дисс... канд. физ,- мат. наук. - М.: ИФЗ АН СССР, 1984. - 18с.

94.Маргулис A.C., Новоселицкий В.М. Прогнозирование слоистого плотностного разреза по гравитационным аномалиям// Докл. АН УССР. Сер. Б,- 1982,-№9,- С.10-13.

95.Маргулис A.C., Новоселицкий В.М. Фурье-аналогия в обратных задачах гравиметрии для некоторых "нефтяных" и "планетарных" плотностных моделей // Изв. АН СССР, Физика Земли, 1982,- № 4,- С.242-246.

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

97.Миков Д.С. Практические приемы интерпретации гравитационных и магнитных аномалий прямыми методами// Разведка недр.-1954,- №5,- С.25-36.

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

99.Мудрецова Е.А., Филатов В.Г. Определение глубины залегания, формы, избыточной плотности и участка модуляции контактной поверхности// Прикладная геофиз. Вып. 78. - М.: Недра, 1975,- С.153-158.

100. Непомнящих A.A., Овчаренко A.B., Ли B.C., Соколов Л.В. Интерпретация гравитационных аномалий на основе пространственного изучения и разделения полей. - Алма-Ата: MB и ССО Каз. ССР, 1978. - 87 с.

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

102. Никонова Ф.И. Использование одного класса потенциалов для приближенного решения обратной задачи гравиразведки и магниторазведки/ Прикладная геофизика. - М.: Недра, 1978,- вып.93,- С.153-156.

103. Никонова Ф.И., Цирульский A.B. Интерпретация гравимагнитных аномалий на основе классов потенциалов, для которых обратная задача разрешима в конечном виде// Изв. АН СССР, Физика Земли, 1978,- № 2,- С.74-85.

104. Новиков П. С. О единственности решения обратной задачи потенциала// Докл. АН СССР, 1938,- т.18, № 3,- С.165-168.

105. Новоселицкий В.М. Интерпретация гравитационных аномалий в условиях латерального изменения плотности осадочных толщ. Автореф. дисс... докт. геол.- мин. наук. - Пермь, ПГУ, 1975. - 34 с.

106. Новоселицкий В.М. К теории определения изменения плотности в горизонтальном пласте по аномалиям силы тяжести// Изв. АН СССР, Физика Земли, 1965,- № 5,- С.25-32.

107. Новоселицкий В.М. Соотношения между переменными пластовыми, средними и эффективными плотностями при решении задач гравиразведки. - В кн.: Прикладная геофизика. - М.: Недра, 1971,- вып.62,-С.151-158.

108. Новоселицкий В.М., Губайдуллин М.Г. Некоторые эквивалентные представления градиентно-слоистых сред в задачах гравиразведки. - В кн.: Теория и методика интерпретации гравимагнитных полей. - К.: Наукова думка, 1981,- С.347-352.

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

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

111. Оганесян С.М. О связи одного итерационного процесса с вариационным методом А.Н. Тихонова - В кн.: Методы решения некорректных задач и их приложение/ Под ред. А.Н.Тихонова и М.М. Лаврентьева. - Новосибирск: 1982.- С.238-239.

112. Оганесян С.М. Об одном условии стабилизатора для решения многомерных обратных задач геофизики// Докл. АН Арм. ССР, 1987.-№3,-С.83

113. Оганесян С.М. Обратная задача гравиметрии в Lp(S) и вопросы приближения элементов из Wp(S). - Докл. АН УССР. Сер. Б,- 1987,- №2,-с.85

114. Оганесян С.М. Обратная задача гравиметрии при заданном носителе масс. - В кн.: Теория и практика геологической интерпретации гравитационных и магнитных аномалий. - ПГО "Казгеофизика", Алма-Ата, 1984,- С.47-49.

115. Оганесян С.М. Обратная задача теории потенциала в пространстве Lp(S). - В кн.: Теория и методы решения некорректно

поставленных задач и их приложения/ Под ред. акад. А.Н.Тихонова. - Изд-во Саратовского университета, 1985,- С.112-114.

116. Оганесян С.М. Решение линейных некорректных задач гравиметрии двойственным методом// Докл. АН УССР. Сер. Б,- 1982,- № 9,-С.13-18.

117. Оганесян С.М. Решение обратной задачи гравиметрии в классе Lp (S) распределения плотностей// Докл. АН УССР. Сер. Б,- 1981,- № 6,-С.39-43.

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

119. Оганесян С.М. Теория и численные методы решения трехмерных задач гравиметрии. Автореф. дисс... доктора физ.-мат. наук. Ленинакан, 1986,- 35с.

120. Оганесян С.М., Оганесян М.Г. Решение модельных обратных задач гравиразведки двойственным методом. - В кн.: Всесоюзное совещание по разработке и комплексированию геофизических методов при детальной и эксплутационной разведке рудных месторождений. - Ереван: Изд-во АН Арм.ССР, 1982,- С.104-105.

121. Оганесян С.М., Старостенко В.И. Двойственный метод решения линейных уравнений некорректных задач, использующих параметрический функционал Лагранжа и вариационный способ А.Н.Тихонова// Докл. АН СССР, 1982,- т.263, № 2,- С.297-301.

122. Оганесян С.М., Старостенко В.И. Некоторые свойства L -псевдорешения и элемента, минимизирующий обобщенный функционал А.Н.Тихонова// Докл. АН СССР, 1982,- т.263, № 3,- С.540-542.

123. Оганесян С.М., Старостенко В.И. О корректности постановки задач геофизики, представленных в виде систем уравнений и итерационных методов // Изв. АН СССР,- Сер. Физика Земли,- 1978,- №8.-С.54-64.

124. Оганесян С.М., Старостенко В.И., Оганесян М.Г. Двойственный метод решения линейных некорректных задач геофизики// Изв. АН СССР, Физика Земли, 1984,- № 6,- С.64-78.

125. Осада М.А. О работах тематической партии 203/64-65 "По обобщению результатов гравиразведочных работ, состоянию сводных гравиразведочных карт на определенные районы". УТГУ Ухтинской геофизической конторы,- Ухта, 1965,- 69 с.

126. Остромогильский А.Х. О единственности решения обратной задачи теории потенциала//ЖВМ и МФ.- 1969,- Т.9., №5.-С.1189-1191.

127. Петровский А.П. Детерминированные модели обратных задач комплексной интерпретации геофизических данных на примере задач сейсмогравиметрии. Автореф. дисс...канд. физ,- мат.наук,- Киев, 1989,- 20 с.

128. Пильник Л.Ф. Отчет Тальбейской с/п 20992 о результатах поисковых работ масштаба 1:100 000 в 1991-93 гг. ПГО "Печорагеофизика", Воркутинская геофизическая экспедиция.-Ухта, 1994,- 88 с.

129. Плескачев О.И. О работах Вуктыльской разведочной партии 298/67-68 в Печорском, Троицко-Печорском районах. - Ухта. УТГУ, 1968. -175 с.

130. Попов Ж.П. Обобщение материалов по физическим свойствам горных пород по Предуральскому прогибу и определение физических параметров по новым площадям. Отчет по теме 308/69-71. Кн.1. Текст отчета. СГЭ №10 УТГУ,- Ухта, 1971,- 281с.

131. Репин Э.М. Отчет Омра-Сойвинской партии 41487 о результатах гравиразведочных работ масштаба 1:25000, проведенных в 1987-88 гг. на Омра-Сойвинской площади. Книга I. Текст отчета. СГЭ №10 "Печорагеофизика",- Ухта, 1988,- 142 с.

132. Соборнов К.О. Зоны клиновидных вдвигов: строение и перспективы нефтегазоносности/ Геология нефти и газа.- 1990,- №8.-С.2-6.

133. Соборнов К.О., Пильник Л.Ф. Гряда Чернышева - фронт вдвиговой пластины?/ Докл. АН СССР,- 1991, т. 317,- №2,- С.430-433.

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

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

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

137. Страхов В.Н. Аналитические продолжение и решение обратной задачи гравиметрии// Изв. АН СССР. Физика Земли,- 1989,- № 3,- С.34-50.

138. Страхов В.Н. Гольдшмидт В.И., Калинина Г.Б., Старостенко В.И. Состояние и перспективы развития в СССР теории интерпретации гравитационных и магнитных полей/ Изв. АН СССР, Физика Земли, 1982.-№ 5,- С.11-30.

139. Страхов В.Н. К вопросу о единственности решения плоской обратной задачи теории потенциала // Изв.АН СССР.- Сер.Физика Земли,-1972.-№2,- С.38-49.

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

141. Страхов В.Н. К теории обратной задачи логарифмического потенциала для контактной поверхности// Изв. АН СССР. Сер. Физика Земли,- 1974,- №6,- С.39-60.

142. Страхов В.Н. К теории плоской задачи гравиметрии и магнитометрии - "аналитический мир", порождаемый выметанием Пуанкаре// Изв. АН СССР, сер. Физика Земли,- 1978,- №2,- С.47-73.

143. Страхов В.Н. Корреляционный метод решения линейной обратной задачи гравиметрии// Докл. АН СССР.- 1990,- 311, т.1.- С.63-66.

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

145. Страхов В.Н. О линейной обратной задаче гравиметрии// Межвузовский сборник научных трудов. - Пенз. политехнический институт, 1990,- № 9 - С.54-67.

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

147. Страхов В.Н. О подходе к решению обратных задач гравиразведки, основанном на теории эквивалентных перераспределений масс// Докл. АН СССР,- 1977,- т.236,- № 3.- С.571-574.

148. Страхов В.Н. О решении линейных обратных задач гравиметрии и магнитометрии// Докл. АН СССР,- 1990,- 310,- №6,- С.1348-1352.

149. Страхов В.Н. О решении линейных обратных задач гравиметрии и магнитометрии при наличии фоновых полей// Докл. АН СССР.- 1990,- 311, № 3,- С.586-589.

150. Страхов В.Н. О решении некорректных задач магнито- и гравиметрии, представленных интегральными уравнениями типа свертки// Изв. АН СССР. Сер. Физика Земли,- 1967,- №4,- С.36-54.

151. Страхов В.Н. Об обратной задаче логарифмического потенциала для контактной поверхности// Изв. АН СССР. Физика Земли. - 1974. - № 2. -С. 43-65.

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

153. Страхов В.Н. Основные идеи и методы извлечения информации из данных гравитационных и магнитных наблюдений. - В кн.: Теория и методы интерпретации гравитационных и магнитных аномалий. М.: Наука, 1979,- С. 146 -269.

154. Страхов В.Н. Основные положения теории конфигурационных и сеточных приближенных решений плоской обратной задачи гравиметрии. -Геология и геофизика,- 1980,- №2,- С. 88-93.

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

156. Страхов В.Н. Решение линейной обратной задачи гравиметрии с учетом структурной информации об искомом решении// Докл. АН СССР,-1990,-311,-№2,- С.331-334.

157. Страхов В.Н. Решение линейных задач геофизики при мультипликативно-аддитивных помехах// Докл. АН СССР.- 1991,- 319, № 3,-С. 559-603.

158. Страхов В.Н., Тетерин Д.Е. Линейные трансформации гравитационных и магнитных аномалий в случае многоэлементных съемок при произвольных сетях наблюдений// Докл. АН СССР,- 1991,- 318, № 3,- С. 572-576.

159. Страхов В.Н. Теория приближенного решения линейных некорректных задач в гильбертовом пространстве и ее использование в разведочной геофизике // Изв. АН СССР. Физика Земли. - 1969. - № 8. - С. 3053.

160. Страхов В.Н. Эквивалентность в обратной задаче гравиметрии и возможности ее практического использования при интерпретации гравитационных аномалий// Изв. АН СССР. Физика Земли,- 1980,- № 2,-С.44-64; 1980.- № 9,- С.38-69.

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

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

163. Тимонин И.И. История развития Печорской плиты в палеозое. Дисс...докт. геол.- мин. наук. Сыктывкар, 1998.- 470 с.

164. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач,- М.: Наука, 1979,- 288 с.

165. Тяпкин К.Ф. О методологических проблемах геологической интерпретации гравитационных и магнитных аномалий// Геоф.сб,- Киев: Наукова думка, 1971,- Вып.61,- С.29-38.

166. Филатов В.Г. О единственности решения некоторых обратных задач гравиразведки// Изв. АН СССР.- Сер.Физика Земли,- 1974,-№11. -С.97-101.

167. Филатов В.Г., Антонов Ю.В. Решение обратной задачи гравиметрии для двух контактных поверхностей// Прикладная геофизика,-М.: Недра, 1979,-Вып. 94.-С. 136-140.

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

169. Цирульский A.B. О решении прямой и обратной задачи гравиметрии// Изв. АН СССР. Сер. Физика Земли.-1974,- №7,- С.84-90.

170. Цирульский A.B., Никонова Ф.И. О разрешимости обратной задачи логарифмического потенциала в конечном виде/ Изв. АН СССР, Физика Земли, 1975,- № 5,- С.37-46.

171. Цирульский A.B., Никонова Ф.И., Федорова Н.В. Метод интерпретации гравитационных и магнитных аномалий с построением эквивалентных семейств решений. - Свердловск: УНЦ АН СССР, 1980. - 136 с.

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

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

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

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

176. Яновская Т.Б., Порохова JI.H. Обратные задачи геофизики. Л., ЛГУ, 1983.-209 с.

177. Backus G., Gilbert F. Numerical applications of a formalist for geophysical inverse problems. - Geophysical Journal of the Royal Astronomical Society, 1967. v.13, p.247-276.

178. Backus G., Gilbert F. The resolving power of gross earth data. -Geophysical Journal of the Royal Astronomical Society, 1968, v.16, p.169-205.

179. Backus G., Gilbert F. Uniqueness in the inversion of inaccurate gross Earth data/ - Philosophical Transactions of the Royal Society of London, 6470, A 266, p.123-142.

180. Jones P. Geologic studies of Cordilleran thrust belt. Denver: Rocky Mount. Ass. Of Geol- 1982,- vol.1.- P.61-75.

181. Week N. Inverse Probleme der Potential theorie. - Appl.analys, 1972, v/2, № 2, p. 105-204.