RU2501047C2 - Способ прогнозирования глубокозалегающих горизонтов на акваториях по результатам тренд-анализа магнитных и гравитационных аномалий - Google Patents
Способ прогнозирования глубокозалегающих горизонтов на акваториях по результатам тренд-анализа магнитных и гравитационных аномалий Download PDFInfo
- Publication number
- RU2501047C2 RU2501047C2 RU2011149904/28A RU2011149904A RU2501047C2 RU 2501047 C2 RU2501047 C2 RU 2501047C2 RU 2011149904/28 A RU2011149904/28 A RU 2011149904/28A RU 2011149904 A RU2011149904 A RU 2011149904A RU 2501047 C2 RU2501047 C2 RU 2501047C2
- Authority
- RU
- Russia
- Prior art keywords
- anomalies
- trend
- magnetic
- horizons
- gravitational
- Prior art date
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Изобретение относится к области геофизики и может быть использовано для построения структурных планов на акваториях: от фундамента до границы М. Для реализации способа используют магнитные, гравитационные поля и рельеф дна моря. Из магнитных и гравитационных аномалий поочередно исключаются локальные аномалии. По остаточным аномалиям полей и рельефу дна оценивается тренд. Находят линию максимальной изменчивости тренда, на которой проводят сейсмические исследования и выделяют глубинные горизонты. Оцениваются многомерные зависимости с глубинами до этих горизонтов. Полученные зависимости распространяются на ближайшее окружение эталонного профиля. Допускается интерполяция и экстраполяция зависимостей. Технический результат: повышение точности разведочных данных.
Description
Известен способ построения геофизических разрезов на акваториях, основанный на объединении результатов интерпретации данных магниторазведки в виде локальных объектов на нескольких соседних параллельных маршрутах [1]. Магнитоактивные горизонты в этом способе выделяются по преобладающему (наиболее вероятному) положению магнитоактивных объектов (верхних кромок, центров масс и т.п.).
Однако для набора достаточной статистики этот способ требует объединения большого числа маршрутов. При изменении глубины исследуемого горизонта это приводит к выбору границ не всегда отвечающих среднему маршруту. Если число суммируемых разрезов уменьшать, то число магнитоактивных объектов тоже уменьшиться до такого уровня, что нельзя будет выделить закономерность в их распределении как вдоль профиля, так и по глубине. Однако для акваторий, где обычно нет глубокого бурения и редкая сеть сейсмических профилей, а сейсмические границы к тому же не всегда отвечают магнитным, иногда и этот способ является основным [1].
Для выполнения надежной интерпретации требуется разделить наблюденные аномалии. Для разделения аномалий обычно используют различные трансформации.
Трансформации аномалий силы тяжести и магнитного поля Земли заключаются в осреднении, сглаживании или пересчете в некоторые функции и служат для разделения наблюденных аномалий на отдельные аномалии или группы аномалий, обусловленных конкретными объектами. Разделение аномалий основано на различии определенных характеристик этих аномалий и на свойстве аддитивности, которое присуще потенциальным полям. При интерпретации магнитных и гравитационных аномалий на акваториях можно предположить, что аномалии, созданные плавными глубинными границами раздела в земной коре, создают независимые от локальных аномалий эффекты, которые при измерении на поверхности Земли складываются (суммируются), а при интерпретации должны быть разделены на составляющие.
Очевидно, что локальные объекты могут присутствовать на больших глубинах, но находятся ли они на границе структурных горизонтов или между ними оценить трудно. Поэтому нельзя их использовать для построения таких структурных горизонтов, которые отличаются по магнитным и плотностным свойствам.
Наша задача заключается в выделении различных аномалий и применении к ним надежных способов интерпретации.
Если некоторые объекты возмущающие аномальное поле находятся на поверхности земли или на дне океана, или изучены другими способами, в том числе бурением и другими горными работами, то естественно при разделении аномалий их необходимо устранить (вычесть) и получить оставшиеся (остаточные) аномалии.
Задача в том и состоит, чтобы в остаточных аномалиях присутствовали лишь аномалии, отвечающие реальному распределению намагниченности или плотности, а аномалии, порожденные неустойчивостью вычислительной операции были бы пренебрежимо малы или случайны [2].
В настоящее время, когда аэромагнитная съемка стала высокоточной, а в комплексе аэро- и надводной магнитных съемок выполняется гравиметрическая съемка, необходимо пересмотреть способы интерпретации и, в первую очередь, - построение структурного плана - глубокозалегающих горизонтов. При высоких точностях съемок требуется по иному подойти к вопросу обнаружения, выделения и разделения с использованием для этой цели локальных аномалий. Для этого пригоден метод вычитания известных аномалий [3]. Тогда остаток (остаточные аномалии), например, в аномалиях магнитного поля Земли, будет являться эффектом от более глубоких слоев (границ), который при использовании стандартных способов интерпретации обычно относился к фону. Эти аномалии должны быть использованы для построения структурных планов глубокозалегающих горизонтов.
Метод вычитания теоретических и/или экспериментальных (к примеру, рассчитанных по результатам бурения или других геофизических методов) локальных аномалий может быть предложен в следующем виде.
Чтобы получить конкретные численные оценки локальных аномалий выполняется подбор магнитных масс, наилучшим образом описывающий магнитные и/или гравитационные аномалии. В качестве магнитоактивных и гравитирующих объектов рассматриваются шары, цилиндры и др. тела. Полученные теоретические локальные аномалии вычитаются из наблюденных.
В полученных (после вычитания локальных аномалий) остаточных аномалиях содержатся региональные аномалии, случайные ошибки и помехи за счет вычитания теоретических аномалий из экспериментальных. Зато в них не появляются составляющие иного (градиенты и др.) физического смысла [3]. Такие аномалии содержат информацию о плавных поверхностях, эффекты от которых могут быть оценены по систематической составляющей (тренд) остаточных аномалий. Систематическая составляющая оценивается путем сглаживания остаточных аномалий по скользящему окну или полиномом заданной степени. В обоих случаях сглаживание производится с учетом весов [4]. Если в исследуемом поле отсутствуют градиентные крупные локальные аномалии, то для разделения наблюденных аномалий на составляющие может быть использован метод оценки систематической составляющей - метод тренда.
Тренд-анализ наблюденных значений уже получил широкое распространение в геологии. Метод основан на аппроксимации исследуемых значений поверхностью двумерной регрессии на координаты точек наблюдений. При этом каждая наблюденная величина f(x, y) рассматривается в виде суммы систематической составляющей φ0(x, y) (тренд) и случайной составляющей (погрешности измерения, погрешности увязки и локальные аномалии) - ε(x, y). При этом наблюденные значения f(x, y) могут быть выражены простым уравнением:
Можно попытаться найти φ(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.
Claims (1)
- Способ прогнозирования глубокозалегающих горизонтов на акваториях по результатам тренд-анализа магнитных и гравитационных аномалий, включающий магниторазведку, гравиразведку и батиметрию, и интерпретацию локальных магнитных и гравитационных аномалий, отличающийся тем, что вычисляют остаточные аномалии по разности исходных и теоретических аномалий, созданных локальными магнитными объектами, производят трансформации остаточных магнитных, и гравитационных аномалий, и рельефа дна моря путем оценки поверхностей тренда на исследуемой площади, по картам тренда находят линию наибольшей изменчивости тренда исследуемых признаков и на этой линии-профиле проводят глубинные сейсмические исследования, выделяют сейсмические горизонты, выбирают коррелирующиеся с ними поверхности тренда различных порядков магнитного поля, гравитационного поля и рельефа дна, строят по ним многомерные зависимости между глубинами конкретных горизонтов и коррелирующихся с ними поверхностей тренда, и по найденным зависимостям прогнозируют глубокозалегающие структурные горизонты на площади исследований.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2011149904/28A RU2501047C2 (ru) | 2011-12-07 | 2011-12-07 | Способ прогнозирования глубокозалегающих горизонтов на акваториях по результатам тренд-анализа магнитных и гравитационных аномалий |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2011149904/28A RU2501047C2 (ru) | 2011-12-07 | 2011-12-07 | Способ прогнозирования глубокозалегающих горизонтов на акваториях по результатам тренд-анализа магнитных и гравитационных аномалий |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2011149904A RU2011149904A (ru) | 2013-06-20 |
RU2501047C2 true RU2501047C2 (ru) | 2013-12-10 |
Family
ID=48784960
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2011149904/28A RU2501047C2 (ru) | 2011-12-07 | 2011-12-07 | Способ прогнозирования глубокозалегающих горизонтов на акваториях по результатам тренд-анализа магнитных и гравитационных аномалий |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2501047C2 (ru) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107229594A (zh) * | 2017-05-05 | 2017-10-03 | 深圳市建设综合勘察设计院有限公司 | 一种多波束测深数据趋势面滤波方法及系统 |
WO2019104912A1 (zh) * | 2017-11-28 | 2019-06-06 | 中国科学院地质与地球物理研究所 | 一种获得剩余重力异常的方法及装置 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112485841A (zh) * | 2019-09-11 | 2021-03-12 | 中国石油天然气集团有限公司 | 深层地层岩性识别方法及装置 |
CN112666607B (zh) * | 2019-10-16 | 2024-06-25 | 中国石油天然气集团有限公司 | 重力反演黄土层厚度分布的方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU693319A1 (ru) * | 1977-08-10 | 1979-10-25 | Научно-Производственное Объединение "Геофизика" | Магнитометр с устройством компенсации магнитных помех от носител |
SU811179A1 (ru) * | 1979-03-28 | 1981-03-07 | Научно-Производственное Объедине-Ние "Геофизика" | Автоматический компенсатор магнит-НыХ пОМЕХ |
RU2008146577A (ru) * | 2008-11-25 | 2010-05-27 | Федеральное государственное унитарное предприятие "ВНИИОкеангеология" (RU) | Способ прогнозирования структурного плана глубокозалегающих горизонтов на акваториях по магнитным и гравитационным аномалиям |
US20100237870A1 (en) * | 2009-03-17 | 2010-09-23 | Geo Equipment Manufacturing Limited | Geophysical Prospecting Using Electric And Magnetic Components Of Natural Electromagnetic Fields |
RU2433429C2 (ru) * | 2010-01-29 | 2011-11-10 | Сергей Борисович Курсин | Способ определения стационарного геомагнитного поля при проведении морской магнитной съемки |
RU2433427C1 (ru) * | 2010-03-16 | 2011-11-10 | Сергей Петрович Алексеев | Способ определения стационарного геомагнитного поля при проведении морской магнитной съемки |
-
2011
- 2011-12-07 RU RU2011149904/28A patent/RU2501047C2/ru not_active IP Right Cessation
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU693319A1 (ru) * | 1977-08-10 | 1979-10-25 | Научно-Производственное Объединение "Геофизика" | Магнитометр с устройством компенсации магнитных помех от носител |
SU811179A1 (ru) * | 1979-03-28 | 1981-03-07 | Научно-Производственное Объедине-Ние "Геофизика" | Автоматический компенсатор магнит-НыХ пОМЕХ |
RU2008146577A (ru) * | 2008-11-25 | 2010-05-27 | Федеральное государственное унитарное предприятие "ВНИИОкеангеология" (RU) | Способ прогнозирования структурного плана глубокозалегающих горизонтов на акваториях по магнитным и гравитационным аномалиям |
US20100237870A1 (en) * | 2009-03-17 | 2010-09-23 | Geo Equipment Manufacturing Limited | Geophysical Prospecting Using Electric And Magnetic Components Of Natural Electromagnetic Fields |
RU2433429C2 (ru) * | 2010-01-29 | 2011-11-10 | Сергей Борисович Курсин | Способ определения стационарного геомагнитного поля при проведении морской магнитной съемки |
RU2433427C1 (ru) * | 2010-03-16 | 2011-11-10 | Сергей Петрович Алексеев | Способ определения стационарного геомагнитного поля при проведении морской магнитной съемки |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107229594A (zh) * | 2017-05-05 | 2017-10-03 | 深圳市建设综合勘察设计院有限公司 | 一种多波束测深数据趋势面滤波方法及系统 |
WO2019104912A1 (zh) * | 2017-11-28 | 2019-06-06 | 中国科学院地质与地球物理研究所 | 一种获得剩余重力异常的方法及装置 |
US11409020B2 (en) | 2017-11-28 | 2022-08-09 | Beijing Csgg Energy Technology Co. Ltd | Method and apparatus for obtaining residual gravity anomaly |
Also Published As
Publication number | Publication date |
---|---|
RU2011149904A (ru) | 2013-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wellmann et al. | Towards incorporating uncertainty of structural data in 3D geological inversion | |
US10329903B2 (en) | Methods of characterizing earth formations using physiochemical model | |
Gallardo et al. | Structure‐coupled multiphysics imaging in geophysical sciences | |
Linde et al. | Joint inversion of crosshole radar and seismic traveltimes acquired at the South Oyster Bacterial Transport Site | |
US8185313B2 (en) | Classifying potential hydrocarbon reservoirs using electromagnetic survey information | |
Pasquet et al. | SWIP: An integrated workflow for surface-wave dispersion inversion and profiling | |
Zhang et al. | Research on the architecture of submarine-fan lobes in the Niger Delta Basin, offshore West Africa | |
Backé et al. | Basin geometry and salt diapirs in the Flinders Ranges, South Australia: Insights gained from geologically-constrained modelling of potential field data | |
Al-Mudhafar | Bayesian kriging for reproducing reservoir heterogeneity in a tidal depositional environment of a sandstone formation | |
US10534100B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on time-lapse seismic data | |
RU2501047C2 (ru) | Способ прогнозирования глубокозалегающих горизонтов на акваториях по результатам тренд-анализа магнитных и гравитационных аномалий | |
Romain et al. | Reconstruction of three-dimensional eolian dune architecture from one-dimensional core data through adoption of analog data from outcrop | |
Ardid et al. | Bayesian magnetotelluric inversion using methylene blue structural priors for imaging shallow conductors in geothermal fields | |
Frahm et al. | Misinterpretation of velocity pull‐ups caused by high‐velocity infill of tunnel valleys in the southern Baltic Sea | |
US20160377752A1 (en) | Method of Digitally Identifying Structural Traps | |
León et al. | GIS-based mapping for marine geohazards in seabed fluid leakage areas (Gulf of Cadiz, Spain) | |
Trainor-Guitton et al. | A methodology for quantifying the value of spatial information for dynamic Earth problems | |
Orellana et al. | Influence of variograms in 3D reservoir-modeling outcomes: An example | |
Gyulai et al. | On the application of combined geoelectric weighted inversion in environmental exploration | |
Darijani et al. | Joint and constrained inversion of magnetic and gravity data: A case history from the McArthur River area, Canada | |
CN105308481B (zh) | 确定地表下温度的方法 | |
Mehdipour et al. | The Best Scenario for Geostatistical Modeling of Porosity in the Sarvak Reservoir in an Iranian Oil Field, Using Electrofacies, Seismic Facies, and Seismic Attributes | |
RU2711178C1 (ru) | Способ 3D прогнозирования свойств и строения геологических объектов на основе компьютерного анализа марковских свойств поверхностных геолого-геофизических полей | |
Zhan et al. | Full elastic finite-difference modeling and interpretation of karst system in a subsalt carbonate reservoir | |
US20190146108A1 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20141208 |