СПОСОБ ПРОГНОЗИРОВАНИЯ ГЛУБОКОЗАЛЕГАЮЩИХ ГОРИЗОНТОВ НА АКВАТОРИЯХ ПО РЕЗУЛЬТАТАМ ТРЕНД-АНАЛИЗА МАГНИТНЫХ И ГРАВИТАЦИОННЫХ АНОМАЛИЙ Российский патент 2013 года по МПК G01V11/00 

Описание патента на изобретение RU2501047C2

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

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

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

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

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

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

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

Задача в том и состоит, чтобы в остаточных аномалиях присутствовали лишь аномалии, отвечающие реальному распределению намагниченности или плотности, а аномалии, порожденные неустойчивостью вычислительной операции были бы пренебрежимо малы или случайны [2].

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

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

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

В полученных (после вычитания локальных аномалий) остаточных аномалиях содержатся региональные аномалии, случайные ошибки и помехи за счет вычитания теоретических аномалий из экспериментальных. Зато в них не появляются составляющие иного (градиенты и др.) физического смысла [3]. Такие аномалии содержат информацию о плавных поверхностях, эффекты от которых могут быть оценены по систематической составляющей (тренд) остаточных аномалий. Систематическая составляющая оценивается путем сглаживания остаточных аномалий по скользящему окну или полиномом заданной степени. В обоих случаях сглаживание производится с учетом весов [4]. Если в исследуемом поле отсутствуют градиентные крупные локальные аномалии, то для разделения наблюденных аномалий на составляющие может быть использован метод оценки систематической составляющей - метод тренда.

Тренд-анализ наблюденных значений уже получил широкое распространение в геологии. Метод основан на аппроксимации исследуемых значений поверхностью двумерной регрессии на координаты точек наблюдений. При этом каждая наблюденная величина f(x, y) рассматривается в виде суммы систематической составляющей φ0(x, y) (тренд) и случайной составляющей (погрешности измерения, погрешности увязки и локальные аномалии) - ε(x, y). При этом наблюденные значения f(x, y) могут быть выражены простым уравнением:

f ( x , y ) = ϕ 0 ( x , y ) + ε ( x , y ) ( 1 )

Можно попытаться найти φ(x, y) близкую к φ0(x, y) и, тем самым, выделить закономерную составляющую - тренд.

Восстановление функции φ0(x, y) через φ(x, y) из f(x, y) может быть произведено путем подбора полинома степени р, который описывает поверхность тренда.

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

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

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

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

По картам тренда исследуемых полей по различным направлениям (например, через каждые 10°, от 0° до 170°) считывают значения по одному или двум основным (магнитное и гравитационное) полям и вычисляют их изменчивость: средние квадратические отклонения а тренда 2-го и более высоких порядков от тренда первого порядка.

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

Пример прогнозирования глубинных горизонтов в районе хребта Ломоносова

В качестве признаков были выбраны глубины моря, аномалии Фая, аномалии Буге с плотностью промежуточного слоя 2.3, 2.67 и переменной плотностью (р2), - оцененной методом Нетлтона, собственно сама плотность р2, и аномалии магнитного поля Земли (МПЗ).

Оказалось, что тренд аномалий Фая практически не информативен для оценки сейсмических границ, а тренд рельефа, аномалий Буге с переменной плотностью и тренд аномалии МПЗ - информативны. К примеру, аномалии Буге с р2 имеют корреляцию с сейсмическими границами (R=0,91-0,89). По информативным признакам строилось множественное уравнение регрессии для прогнозирования глубинных горизонтов по всей площади гравимагнитных съемок на площади примерно 100×700 км.

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

Коэффициенты корреляции сейсмических горизонтов с трендом МПЗ Т R1 R2 R3 R4 R5 R6 R7 R8 R9 Н1 0.48 0.49 0.64 0.63 0.63 0.62 0.54 0.50 0.46 0.45 Н2 0.39 0.84 0.72 0.93 0.82 0.85 0.77 0.72 0.64 0.59 Н3 0.33 0.82 0.87 0.94 0.95 0.91 0.86 0.78 0.71 0.65 Н4 -0.07 0.10 -0.48 -0.07 -0.47 -0.28 -0.46 -0.36 -0.42 -0.38 Н5 0.11 0.63 0.87 0.74 0.86 0.70 0.70 0.57 0.49 0.42 где Т - аномальное поле Земли; R1÷R9 - весовые тренды МПЗ с 1-го по 9-й порядок; H1-H5 - глубины сейсмических границ.

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

Таким образом, предложен способ прогнозирования глубокозалегающих структурных горизонтов на акваториях по результатам тренд-анализа магнитных и гравитационных аномалий, полученных по данным магниторазведки, гравиразведки и батиметрии, и результатам интерпретации локальных магнитных и гравитационных аномалий. По результатам интерпретации локальных аномалий вычисляются теоретические аномалии и строятся разности исходных и теоретических (вычисленных по результатам интерпретации) локальных аномалий - остаточные аномалии. Производят трансформации остаточных магнитных, гравитационных и рельефных аномалий путем оценки поверхностей тренда (уравнения двумерной регрессии на координаты точек наблюдения x и y) на исследуемой площади. По картам тренда в секторе от 0° до 180° находят линию (профиль) наибольшей изменчивости тренда исследуемых признаков и на этой линии-профиле проводят глубинные сейсмические исследования (зондирования) (ГСЗ) для оценки глубинных структурных горизонтов. Выделяют сейсмические горизонты, выбирают коррелирующиеся с ними поверхности тренда среди различных порядков и различных признаков (магнитное поле, гравитационное поле и рельеф дна). По выбранным для каждого структурного горизонта коррелирующимся признакам (трендам) строят многомерные зависимости между глубинами конкретных горизонтов и коррелирующихся с ними поверхностей тренда, и по найденным зависимостям прогнозируют глубоказалегающие структурные горизонты на сейсмическом (для оценки надежности) и пересекаемым им геофизическим (гравимагнитным) профилям, или на всей площади гравимагнитных исследований. Если район гравимагнитных и ГСЗ исследований имеет однородное геологическое строение, то найденные зависимости могут быть использованы для прогнозирования структурных горизонтов на больших территориях (акваториях). Во всяком случае, всегда имеется возможность провести в небольшом объеме дополнительные сейсмические исследования на границе исследуемого района для оценки надежности прогнозирования, или площадь прогнозирования «выйдет» на другой сейсмический профиль, где имеются данные о глубинном строении.

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

Список литературы

1. В.Э. Волк «Геофизические методы исследований Мирового океана», Ленинград, НИИГА, 1979, с.82.

2. Справочник геофизика, 5 т., Гравиразведка, Из-во: «Недра», 1968, с.512.

3. В.К. Паламарчук «Разделение магнитных аномалий путем исключения теоретических эффектов тел-помех», ж. «Геология и Геофизика» №2, Новосибирск, Наука, 1986.

4. В.К. Паламарчук «Опыт разделения аномалии методом тренда», ж. «Геология и Геофизика» №4, Новосибирск, Наука, 1972.

Похожие патенты RU2501047C2

название год авторы номер документа
СПОСОБ РАСШИРЕНИЯ РАЙОНОВ ПРИМЕНЕНИЯ НАВИГАЦИИ ПО МАГНИТНОМУ ПОЛЮ 2010
  • Паламарчук Василий Климентьевич
  • Бурдакова Елена Владиславовна
  • Глинская Надежда Викторовна
  • Мищенко Оксана Николаевна
  • Прялухина Любовь Александровна
RU2447405C2
СПОСОБ ПРОГНОЗИРОВАНИЯ СТРУКТУРНОГО ПЛАНА ГЛУБОКОЗАЛЕГАЮЩЕГО ОТРАЖАЮЩЕГО ГОРИЗОНТА ОСАДОЧНОГО ЧЕХЛА 1998
  • Щекин С.Н.
  • Нежданов А.А.
  • Туренков Н.А.
  • Миколаевский Э.Ю.
RU2145106C1
Способ картирования магнитных аномалий-индикаторов залежей углеводородов 2015
  • Паламарчук Василий Климентьевич
  • Глинская Надежда Викторовна
  • Мищенко Оксана Николаевна
  • Бурдакова Елена Владиславовна
  • Субботин Константин Петрович
RU2614351C2
Способ прогноза залежей углеводородов 2021
  • Нургалиев Данис Карлович
  • Хасанов Дамир Ирекович
  • Кузина Диляра Мтыгулловна
  • Зиганшин Эдуард Рашидович
RU2781752C1
СПОСОБ УЧЕТА ВАРИАЦИЙ ГЕОМАГНИТНОГО ПОЛЯ ПО ДОПОЛНИТЕЛЬНЫМ СЕКУЩИМ МАРШРУТАМ ПРИ ВЫПОЛНЕНИИ МАГНИТНЫХ СЪЕМОК НА АКВАТОРИЯХ 2011
  • Глинская Надежда Викторовна
  • Мищенко Оксана Николаевна
  • Паламарчук Василий Климентьевич
  • Бурдакова Елена Владиславовна
RU2539097C2
СПОСОБ ГЛУБИННОГО СТРУКТУРНОГО КАРТИРОВАНИЯ 1994
  • Керимов И.А.
  • Гайсумов М.Я.
  • Гайрабеков И.Г.
  • Моллаев З.Х.
RU2107314C1
СПОСОБ ПРОГНОЗА ЕМКОСТНЫХ ПАРАМЕТРОВ И ТИПА ФЛЮИДОНАСЫЩЕНИЯ КОЛЛЕКТОРОВ 2013
  • Тригубович Георгий Михайлович
  • Филатов Владимир Викторович
  • Багаева Татьяна Николаевна
  • Яковлев Андрей Георгиевич
  • Яковлев Денис Васильевич
  • Агафонов Юрий Александрович
  • Шарлов Максим Валерьевич
RU2540216C1
СПОСОБ ОБНАРУЖЕНИЯ И ОТОБРАЖЕНИЯ ФИГУРЫ ГАЗОНЕФТЯНОЙ ЛОГ-ТРУБКИ 2008
  • Пшеничный Геннадий Андреевич
  • Каминский Евгений Юрьевич
  • Матвеев Юрий Иванович
  • Свинин Михаил Васильевич
  • Бубнов Николай Геннадьевич
  • Крафт Якоб Соломонович
RU2401443C2
СПОСОБ МОРСКОЙ ВЫСОКОТОЧНОЙ МАГНИТНОЙ СЪЕМКИ 2015
  • Паламарчук Василий Климентьевич
  • Глинская Надежда Викторовна
  • Мищенко Оксана Николаевна
  • Бурдакова Елена Владиславовна
RU2665355C2
СПОСОБ ПРОГНОЗИРОВАНИЯ СТРУКТУРНОГО ПЛАНА ГЛУБОКОПОГРУЖЁННЫХ ГРАНИЦ РАЗДЕЛА ОСАДОЧНОГО ЧЕХЛА 2002
  • Малышев В.В.
RU2207602C1

Реферат патента 2013 года СПОСОБ ПРОГНОЗИРОВАНИЯ ГЛУБОКОЗАЛЕГАЮЩИХ ГОРИЗОНТОВ НА АКВАТОРИЯХ ПО РЕЗУЛЬТАТАМ ТРЕНД-АНАЛИЗА МАГНИТНЫХ И ГРАВИТАЦИОННЫХ АНОМАЛИЙ

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

Формула изобретения RU 2 501 047 C2

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

Документы, цитированные в отчете о поиске Патент 2013 года RU2501047C2

RU 2008146577 А, 27.05.2010
СПОСОБ ОПРЕДЕЛЕНИЯ СТАЦИОНАРНОГО ГЕОМАГНИТНОГО ПОЛЯ ПРИ ПРОВЕДЕНИИ МОРСКОЙ МАГНИТНОЙ СЪЕМКИ 2010
  • Курсин Сергей Борисович
  • Бродский Павел Григорьевич
  • Добротворский Александр Николаевич
  • Ставров Константин Георгиевич
  • Леньков Валерий Павлович
  • Жуков Юрий Николаевич
  • Ленькова Людмила Александровна
  • Чернявец Владимир Васильевич
  • Румянцев Юрий Владимирович
RU2433429C2
СПОСОБ ОПРЕДЕЛЕНИЯ СТАЦИОНАРНОГО ГЕОМАГНИТНОГО ПОЛЯ ПРИ ПРОВЕДЕНИИ МОРСКОЙ МАГНИТНОЙ СЪЕМКИ 2010
  • Алексеев Сергей Петрович
  • Курсин Сергей Борисович
  • Добротворский Александр Николаевич
  • Бродский Павел Григорьевич
  • Ставров Константин Георгиевич
  • Леньков Валерий Павлович
  • Чернявец Владимир Васильевич
  • Жуков Юрий Николаевич
  • Румянцев Юрий Владимирович
RU2433427C1
Магнитометр с устройством компенсации магнитных помех от носителя 1977
  • Хвостов Орион Павлович
  • Попов Михаил Сергеевич
  • Митяев Игорь Михайлович
  • Кондратьев Геннадий Яковлевич
SU693319A1
Автоматический компенсатор магнит-НыХ пОМЕХ 1979
  • Хвостов Орион Павлович
  • Аверкиев Владимир Витальевич
  • Бебекин Виктор Анатольевич
  • Горшков Владимир Викторович
  • Колосков Юрий Петрович
  • Попов Михаил Сергеевич
SU811179A1
Способ измерения магнитных помех носителя магнитометра 1991
  • Хвостов Орион Павлович
SU1803898A1
US 20100237870 A1, 23.09.2010.

RU 2 501 047 C2

Авторы

Глинская Надежда Викторовна

Мищенко Оксана Николаевна

Паламарчук Василий Климентьевич

Бурдакова Елена Владиславовна

Шарков Дмитрий Витальевич

Даты

2013-12-10Публикация

2011-12-07Подача