RU2664488C1 - Способ автоматического получения структурного строения из данных потенциального поля - Google Patents

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

Info

Publication number
RU2664488C1
RU2664488C1 RU2017133223A RU2017133223A RU2664488C1 RU 2664488 C1 RU2664488 C1 RU 2664488C1 RU 2017133223 A RU2017133223 A RU 2017133223A RU 2017133223 A RU2017133223 A RU 2017133223A RU 2664488 C1 RU2664488 C1 RU 2664488C1
Authority
RU
Russia
Prior art keywords
data
magnetic
structural
structural structure
order
Prior art date
Application number
RU2017133223A
Other languages
English (en)
Inventor
Дяньхуа ЦАО
Original Assignee
Инститьют Оф Минерал Рисорсиз, Чайниз Акедеми Оф Джиолоджикал Сайенсиз
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from CN201510096981.9A external-priority patent/CN104658037B/zh
Priority claimed from CN201510303316.2A external-priority patent/CN104965232B/zh
Application filed by Инститьют Оф Минерал Рисорсиз, Чайниз Акедеми Оф Джиолоджикал Сайенсиз filed Critical Инститьют Оф Минерал Рисорсиз, Чайниз Акедеми Оф Джиолоджикал Сайенсиз
Application granted granted Critical
Publication of RU2664488C1 publication Critical patent/RU2664488C1/ru

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38Processing data, e.g. for analysis, for interpretation, for correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
    • G01V11/002Details, e.g. power supply systems for logging instruments, transmitting or recording data, specially adapted for well logging, also if the prospecting method is irrelevant
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/081Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V7/00Measuring gravitational fields or waves; Gravimetric prospecting or detecting
    • G01V7/02Details
    • G01V7/06Analysis or interpretation of gravimetric records
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20016Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20036Morphological image processing
    • G06T2207/20044Skeletonization; Medial axis transform
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20064Wavelet transform [DWT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20172Image enhancement details
    • G06T2207/20182Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Electromagnetism (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • Computer Graphics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

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

Description

Область техники, к которой относится изобретение
Настоящее изобретение относится к способу автоматического получения структурного строения из данных гравитационного потенциального поля и данных магнитного потенциального поля, который представляет собой методику обнаружения геологической структуры на основе данных гравитационных аномалий и данных магнитных аномалий. В частности, настоящее изобретение относится к областям вейвлет-анализа, обработки изображений, геофизики, геологии и разведки полезных ископаемых. Способ настоящего изобретения может быть непосредственно применен в областях разведки полезных ископаемых и родственной геологической разведки.
Предпосылки создания изобретения
С развитием технологии воздушной геофизической разведки способы гравитационных и магнитных измерений (ниже называемые «гравитационно-магнитными») приобрели преимущества экономичности, оперативности и охвата труднодостижимых ландшафтных зон, и они играют все более и боле важную роль в разведке и оценке месторождений металлических руд. В особенности с развитием методики высокоточной аэромагнитной съемки, способы обнаружения геологических структур, обусловливающих образование рудных месторождений, на основании данных магнитных аномалий являются весьма важными на всех этапах разведки, от выбора района до выбора объектов. В настоящее время в области разведки полезных ископаемых гравитационно-магнитная разведка обычно используется для непосредственного обнаружения оруденелых объектов, обладающих сильными магнитными свойствами или высокой плотностью, и для интерпретации и реконструкции геологических структур с сильными аномалиями. Большинство месторождений руд цветных металлов и редких металлов не могут непосредственно образовывать заметные гравитационно-магнитные аномалии. Однако геологические структуры, обусловливающие образование рудных месторождений, могут порождать аномальные отклонения. Поэтому обнаружение геологических структур, обусловливающих образование рудных месторождений, на основе данных гравитационного потенциального поля и данных магнитного потенциального поля является важным для определения местоположения объектов разведки месторождений металлических руд.
В настоящее время способы автоматической идентификации и получения сведений о структуре из гравитационно-магнитных потенциальных полей, главным образом, включают способ аналитического сигнала, способ деконволюции Эйлера, способ многоуровневого обнаружения краев потенциального поля и т. д. Эти способы обладают недостатками, связанными с тем, что они являются нечувствительными к сведениям об ориентации и не могут получать полные и точные местоположения границ аномалий. В заявке на патент Китая №200810006676.6, озаглавленной «Multidirectional and Multiscale Edge Detection Method of Potential Field (Способ многонаправленного и многоуровневого обнаружения краев потенциального поля)» заявителя настоящего изобретения раскрыт способ многонаправленного и многоуровневого обнаружения краев потенциального поля, улучшающий сведения об ориентации при помощи направленного вейвлет-преобразования, получающий сведения о границах источников аномалий в разных направлениях, выполняющий автоматическое извлечение структурного строения, преодолевающий недостаток, связанный с нечувствительностью к сведениям об ориентации в способе аналитического сигнала, способе деконволюции Эйлера и способе многоуровневого обнаружения краев потенциального поля, и представляющий собой техническое решение для оперативной реконструкции трехмерных структур поверхностных корок на малых глубинах на основании данных потенциальных полей. Однако данный способ по-прежнему обладает следующими недостатками: (1) края, получаемые посредством вычисления, не обладают шириной точек в один пиксел, фактические ареалы, соответствующие краям и пересечениям краев в разных направлениях, являются большими, что приводит к невысокой точности результата анализа; (2) трехмерная структурная карта поверхностных корок на малых глубинах получается путем искусственной векторизации на разные уровни краев и принятия в качестве границ аномалий линий, соединяющих центры краев, что приводит к низкой производительности; (3) трехмерная структурная карта поверхностных корок на малых глубинах показывает лишь сведения о структурах на разных глубинах, не может отражать изменения литологического состава структурных поясов и с двух сторон, и не может указывать интенсивность деформации и активности структур; и (4) уровни определяются нечетко, т. е. признакам уровня не придается ни одно особое геофизическое свойство. Таким образом, существует потребность в способе, пригодном для точного получения структурного строения из потенциального поля на основе собранных данных потенциального поля.
Кроме того, по сравнению с гравитационными аномалиями, магнитные аномалии являются более сложными из-за влияния ориентации намагниченности. Обычно полная напряженность магнитного поля (Total magnetic intensity, TMI), регистрируемая современным магнитометром, представляет собой составляющую в направлении, параллельном направлению постоянного магнитного поля Земли. Аномалии TMI можно получить путем выполнения коррекции TMI по Международному эталону геомагнитного поля. Из-за влияния косой намагниченности аномалии TMI обладают недостатками, связанными, например, с боковым отклонением, деформацией картины, а также с чередованием положительных и отрицательных значений. Обычно влияние вышеуказанных факторов исключается при помощи обработки данных аномалий TMI методом приведения к полюсу. При приведении к полюсу наблюдаемые аномалии TMI преобразуются в вертикальные магнитные аномалии в состоянии вертикальной намагниченности, т. е. наблюдаемые аномалии TMI преобразуются в аномалии, которые могут быть измерены на северном магнитном полюсе, так что магнитные аномалии перемещаются в местоположения непосредственно над их источниками, что способствует геологической интерпретации магнитных аномалий.
Однако под влиянием малого магнитного наклонения и шума трудно получить надежные данные магнитных аномалий с вертикальной намагниченностью путем выполнения приведения к полюсу для данных аномалий TMI в низкоширотных зонах (которые обычно относятся к магнитному наклонению в диапазоне ±20°). Современные способы автоматической идентификации и получения сведений о магнитной структуре из данных аномалий TMI после приведения к полюсу, такие как способ аналитического сигнала, способ деконволюции Эйлера, симметрия фаз и способ многоуровневого обнаружения краев данных потенциального поля, не приспособлены для автоматической идентификации сведений о магнитной структуре в низкоширотных зонах.
Поэтому существует потребность в создании способа, способного точно получать геологическое структурное строение на основе данных магнитных аномалий в низкоширотных зонах.
Сущность изобретения
Целью настоящего изобретения является создание способа автоматического получения структурного строения из данных гравитационно-магнитных потенциальных полей с целью оперативного получения сведений о геологических структурах, обусловливающих образование рудных месторождений, с получением определения местоположения объектов разведки — месторождений металлических руд. Настоящее изобретение дополнительно предусматривает способ автоматического получения магнитного структурного строения в низкоширотных зонах, что позволяет получать не только линеаменты, но также и кольцевые структуры. В настоящем изобретении полученная таким образом структура определяется как структурное строение.
Для достижения данной цели настоящее изобретение предусматривает способ, включающий получение гравитационного структурного строения или магнитного структурного строения путем применения многонаправленного и многоуровневого обнаружения краев к измеренным данным гравитационного потенциального поля и магнитного потенциального поля, утончение полученного структурного строения на каждом уровне до однопиксельной ширины с использованием алгоритма определения морфологического скелета; выделение сведений о структурах на разных глубинах путем представления структурных строений на разных уровнях разными постепенно изменяющимися цветами и их наложения одно на другое с целью создания составной карты структурного строения; и выделение сведений об интенсивности изменения плотности и изменения напряженности магнитного поля путем представления значений структурной интенсивности, полученных из модуля градиента, соответствующего каждой краевой точке на структурных строениях на разных уровнях, разными постепенно изменяющимися цветами и их наложения одно на другое с целью создания составной карты структурного строения с отображением интенсивности. Таким образом, решение настоящего изобретения может получать сведения о гравитационно-магнитных аномалиях на разных глубинах в зоне исследований, сведения о распределении структурных строений, представляющих геологические структуры на разных глубинах, и сведения об интенсивности изменения плотности и изменения напряженности магнитного поля структурных строений на разных глубинах с целью идентификации и качественной интерпретации геологических структур, обусловливающих образование рудных месторождений, и может определять типы потенциальных рудных месторождений и свойства структур, обусловливающих образование рудных месторождений, в соответствии с априорными знаниями о зоне исследований, для отбора разных типов структурных строений с целью определения местоположения объектов разведки — месторождений металлических руд.
Настоящее изобретение предусматривает способ автоматического получения структурного строения потенциального поля, включающий следующие этапы:
1) Выполнение предварительной обработки гравитационных данных или магнитных данных, измеренных в зоне исследований.
Предварительная обработка включает приведение магнитных данных к полюсу с целью получения приведенных к полюсу магнитных аномалий или псевдогравитационное преобразование магнитных данных в псевдогравитационные аномалии; или предварительную обработку гравитационных данных с целью получения гравитационных аномалий Буге.
2) Применение многонаправленного обнаружения краев данных потенциального поля, соответственно, в отношении предварительно обработанных данных гравитационного потенциального поля или данных магнитного потенциального поля на множестве уровней, которое включает продолжение вверх данных потенциального поля на множество предварительно определенных высот с целью получения данных гравитационного потенциального поля или данных магнитного потенциального поля на множестве уровней, и применение многонаправленного обнаружения краев, соответственно, к полученным данным потенциального поля на каждом уровне с целью получения краев потенциального поля уровня.
За счет выбора разных направлений a для применения обнаружения краев для каждого уровня, продолженного вверх на множество предварительно определенных высот, можно выделить сведения о краях в разных направлениях. Для достижения полного покрытия значение a направления полагают равным kπ/(2n-1), где k=0, 1, 2…(2n-1), и n — целое число, большее или равное 2. В каждом направлении точки с локальным максимумом модуля двумерного направленного вейвлет-преобразования предварительно обработанных гравитационных данных или магнитных данных соединяют в направлении, перпендикулярном градиенту, с целью образования кривых и создания краев. Получают края, вычисленные на множестве разных направлений на одном уровне, в качестве краев на уровне принимают объединение всех полученных краев, и, таким образом, может быть получено многонаправленное обнаружение краев данных потенциального поля на всех уровнях.
3) Утончение вычисленных краев на уровнях до однопиксельной ширины, соответственно, с использованием алгоритма определения морфологического скелета с целью получения структурных строений на уровнях.
4) Наложение вычисленных структурных строений уровней одно на другое с целью создания составной карты структурного строения.
Полученные структурные строения на всех уровнях накладывают одно на другое с целью создания составной карты структурного строения, отражающей сведения на разных глубинах. Боковое отклонение краев на разных уровнях карты отражает сведения об залегании структурных строений.
Путем приведения краев, полученных из потенциальных полей, продолженных вверх на разные высоты, в соответствие структурам на разных глубинах, где глубина составляет половину высоты после продолжения вверх (см. случай А в публикации «Upward Continuation as a Standard Separation Filter for Potential-Field Maps (Продолжение вверх как стандартный фильтр разделения для карт потенциальных полей)», Jacobsen, B.H., Geophysics, т. 52, №. 8, 1987), можно получить карты структурных строений, представляющие разные глубины.
5) Принятие модуля градиента в каждой краевой точке на вычисленных краях в качестве значения интенсивности этой краевой точки на структурном строении. Величина значений интенсивности краевых точек на структурном строении, которым приписаны значения интенсивности, отражает величину изменений литологического состава структурного пояса и двух сторон, и отражает интенсивность деформации и активности структуры. Путем создания карт строения с отображением интенсивности разных уровней, отражающих сведения на разных уровнях, и путем наложения карт строения с отображением интенсивности разных уровней с целью создания составной карты структурного строения с отображением интенсивности, раскрывается региональное геологическое тектоническое строение.
Кроме того, на полные горизонтальные производные наклонных производных (tilt derivatives, TDR) и аналитические сигналы данных аномалий TMI не оказывает влияние магнитное наклонение, результаты вычислений не связаны с величиной магнитного наклонения, и вычисление полных производных TDR и аналитических сигналов может быть выполнено непосредственно в отношении данных аномалий TMI без обработки приведения к полюсу. Однако результаты вычисления по способу аналитического сигнала могут увеличивать диапазон магнитных аномалий, терять сведения о залегании геологической структуры и сведения о структурном зонировании, и являться нечувствительными в отношении идентификации геологических структур. Однако текущие результаты полных горизонтальных производных наклонных производных выражены в форме изображений с нанесением на координатную сетку или контурных карт, и полезные сведения для геологической интерпретации и разведки полезных ископаемых, такие как глубина структуры, отношение первичных и вторичных структур, отношение врезания и напряженность магнитного поля, представлены быть не могут. Традиционный способ обнаружения краев на основе горизонтального градиента без учета сведений о направленности из данных не позволяет получать полные и точные местоположения границ объектов-источников магнитных аномалий. Путем применения многонаправленного обнаружения краев на основе общих горизонтальных производных наклонных производных на множестве уровней к предварительно обработанным данным аномалий TMI настоящее изобретение может эффективно идентифицировать и устанавливать магнитные структурные строения в низкоширотных зонах. Следовательно, способ автоматического получения структурного строения в соответствии с настоящим изобретением в особенности является применимым для автоматического получения структурных строений из данных магнитных измерений в низкоширотной зоне.
В соответствии с одним аспектом настоящего изобретения настоящее изобретение предусматривает способ автоматического получения структурного строения из данных потенциального поля, включающий следующие этапы:
предварительную обработку данных потенциального поля из зоны исследований;
продолжение вверх предварительно обработанных данных потенциального поля на множество предварительно определенных высот с целью получения множества данных потенциального поля на соответствующих уровнях;
применение многонаправленного обнаружения краев, соответственно, к множеству данных потенциального поля на уровнях с целью получения множества краев потенциального поля на соответствующих уровнях; и
утончение полученных краев потенциального поля на уровнях, соответственно, до однопиксельной ширины с использованием алгоритма определения морфологического скелета с целью получения множества карт структурного строения на соответствующих уровнях.
Предпочтительно способ дополнительно включает наложение одна на другую полученного множества карт структурного строения на соответствующих уровнях для создания составной карты структурного строения.
Предпочтительно способ дополнительно включает использование модуля градиента каждой краевой точки на карте структурного строения на каждом уровне в качестве значения интенсивности краевой точки на карте структурного строения на уровне с целью получения множества карт структурного строения с отображением интенсивности на соответствующих уровнях.
Предпочтительно способ дополнительно включает наложение множества карт структурного строения с отображением интенсивности на соответствующих уровнях одна на другую с целью создания составной карты структурного строения с отображением интенсивности.
Предпочтительно данные потенциального поля представляют собой данные гравитационного потенциального поля или данные магнитного потенциального поля, и этап предварительной обработки дополнительно включает: предварительную обработку гравитационных данных с целью получения гравитационных аномалий Буге; или приведение магнитных данных к полюсу с целью получения приведенных к полюсу магнитных аномалий или применение к магнитным данным псевдогравитационного преобразования с целью получения псевдогравитационных аномалий.
Предпочтительно этап применения многонаправленного обнаружения краев к данным потенциального поля на каждом уровне включает следующие этапы:
в предположении, что уровень s=z/z0, z>z0, z0 — это высота измерения, и z — высота, увеличивающаяся вверх, — определение гравитационной аномалии или магнитной аномалии в местоположении (x, y) на нулевой высоте как f0(x, y),
определение сглаживающей функции на уровне s как:
Figure 00000001
,
где
Figure 00000002
, k(x, y, z) — функция Грина,
определение вейвлет-функции в направлении a как:
Figure 00000003
,
где D представляет производную первого порядка;
для уровня s и местоположения (x, y), — определение вейвлет-преобразования гравитационной аномалии или магнитной аномалии, f0(x, y), в направлении a как:
Figure 00000004
где * представляет операцию свертки,
при этом известно, что
Figure 00000005
в соответствии с формулой продолжения потенциального поля вверх, и
fz(x, y) — гравитационная аномалия или магнитная аномалия посредством продолжения вверх f0(x, y) от нулевой высоты до высоты z=sz0 и полученная путем продолжения вверх гравитационной аномалии или магнитной аномалии, fz0(x, y), измеренной на высоте z0 измерения, на высоту z-z0,
отсюда
Figure 00000006
Figure 00000007
;
кроме того, для уровня s и местоположения (x, y), — определение вейвлет-преобразования гравитационной аномалии или магнитной аномалии, f0(x, y), в направлении a+
Figure 00000008
как:
Figure 00000009
,
поэтому двумерное вейвлет-преобразование f0(x, y) записывается как в градиенте:
Figure 00000010
где
Figure 00000011
обозначает двумерный градиент, и
при этом для местоположения (x, y), уровня s и направления a двумерное направленное вейвлет-преобразование
Figure 00000012
f0(x, y) прямо пропорционально градиенту
Figure 00000013
fz(x, y), и двумерное направленное вейвлет-преобразование
Figure 00000012
f0(x, y) можно представить градиентом
Figure 00000013
fz(x, y),
для высоты z, — определение модуля градиента
Figure 00000014
как:
Figure 00000015
,
где соответствующий аргумент градиента вдоль горизонтальной оси имеет вид:
Figure 00000016
, и
определение точки, в которой модуль M
Figure 00000017
представляет собой локальный максимум вдоль направления градиента
Figure 00000018
, как краевой точки,
формирование кривой путем соединения точек с локальным максимумом модуля градиента в направлении, перпендикулярном градиенту, с целью создания края, и
для той же высоты вычисление краев на множестве разных направлений a и принятие объединения вычисленных краев как краев потенциального поля на соответствующих уровнях.
Предпочтительно, для того чтобы полностью покрыть двумерную плоскость, этап вычисления краев на множестве разных направлений a для одной и той же высоты дополнительно включает: принятие значения a направления равным kπ/(2n-1), где k=0, 1, 2…(2n-1), и n — целое число, большее или равное 2.
Предпочтительно края, полученные из потенциальных полей, продолженных вверх на предварительно определенные высоты, соответствуют структурам на разных глубинах, и составную карту структурного строения, отражающую сведения на разных глубинах, получают путем наложения карт структурного строения уровней одна на другую.
Предпочтительно для создания составной карты структурного строения высоты представляют с использованием постепенно изменяющиеся цвета.
Предпочтительно для создания составной карты структурного строения с отображением интенсивности величины значений интенсивности представляют с использованием постепенно изменяющихся цветов.
В соответствии с другим аспектом настоящего изобретения настоящее изобретение предусматривает способ автоматического получения структурного строения, включающий следующие этапы:
предварительную обработку магнитных данных или данных гравитационных измерений из зоны исследований с целью получения данных аномалий полной напряженности магнитного поля (TMI) или данных гравитационных аномалий Буге;
нанесение на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге и продолжение вверх нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге на множество предварительно определенных высот с целью получения множества нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге, Th, на соответствующих уровнях, где h представляет высоту после продолжения вверх;
соответствующее вычисление наклонной производной TDRh для данных аномалий TMI или данных гравитационных аномалий Буге на каждом уровне из нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге, Th, на уровне;
применение многонаправленного обнаружения краев на основе горизонтального градиента, соответственно, к наклонной производной нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге на каждом уровне с целью получения краев объектов-источников магнитных или гравитационных аномалий на уровнях; и
утончение полученных краев объектов-источников магнитных или гравитационных аномалий на уровнях, соответственно, до однопиксельной ширины с использованием алгоритма определения морфологического скелета с целью получения множества карт структурного строения на соответствующих уровнях.
Предпочтительно способ дополнительно включает наложение множества карт структурного строения на соответствующих уровнях с целью создания составной карты структурного строения.
Предпочтительно способ дополнительно включает приведение краев, полученных из нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге, продолженных вверх на множество определенных высот, в соответствие структурам на разных глубинах, и наложение полученных карт структурного строения на глубинах с целью получения составной карты структурного строения, отражающей сведения на разных глубинах врезания.
Предпочтительно этап применения многонаправленного обнаружения краев на основе горизонтального градиента, соответственно, к наклонной производной нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге на каждом уровне включает следующие этапы:
соответственно определение направленных производных наклонной производной TDRh в направлениях a и a+
Figure 00000008
как:
Figure 00000019
Figure 00000020
где D представляет производную первого порядка;
для высоты h и направления a, — представление горизонтального градиента наклонной производной TDRh как:
Figure 00000021
где
Figure 00000011
— горизонтальный градиент;
определение модуля горизонтального градиента
Figure 00000022
как:
Figure 00000023
где аргумент горизонтального градиента имеет вид:
Figure 00000024
, и
затем определение точек, в которых модуль
Figure 00000025
представляет собой локальный максимум вдоль направления
Figure 00000026
градиента, как краевых точек объектов-источников магнитных или гравитационных аномалий на высоте h в направлении a;
в каждом направлении a формирование кривой путем соединения точек с локальным максимумом модуля горизонтального градиента наклонной производной TDRh в направлении, перпендикулярном градиенту, с целью создания края и
для той же высоты вычисление краев на множестве разных направлений a и принятие объединения краев в качестве краев объектов-источников магнитных или гравитационных аномалий на соответствующих уровнях,
где для того, чтобы полностью покрыть двумерную плоскость, множество разных значений a направления соответственно полагают равными kπ/(2n-1), где k=0, 1, 2…(2n-1), и n — целое число, большее или равное 2.
Предпочтительно способ дополнительно включает, соответственно, использование модуля горизонтального градиента в каждой краевой точке на карте структурного строения на каждом уровне для представления структурной глубины залегания в краевой точке на карте структурного строения на уровне с целью получения множества карт структурного строения, представляющих структурные глубины залегания на соответствующих уровнях.
Предпочтительно способ дополнительно включает наложение множества карт структурного строения, представляющих структурные глубины залегания на соответствующих уровнях, одна на другую с целью создания составной карты структурного строения с отражением глубины залегания.
Предпочтительно способ дополнительно включает, соответственно, вычисление трехмерного аналитического сигнала ASh на основе нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге, Th, на каждом уровне для получения значения ASh каждой краевой точки с целью получения множества карт структурного строения, представляющих напряженность магнитного поля или интенсивность плотности в краевых точках на соответствующих уровнях.
Предпочтительно способ дополнительно включает наложение множества карт структурного строения, представляющих напряженность магнитного поля или интенсивность плотности в краевых точках на соответствующих уровнях, одна на другую с целью создания составной карты структурного строения с отображением напряженности магнитного поля или составной карты структурного строения с отображением интенсивности плотности.
Предпочтительно способ дополнительно включает выполнение обработки, связанной с подавлением шума, в отношении наклонных производных перед обнаружением краев.
Предпочтительно способ применим для автоматического получения магнитного структурного строения в низкоширотной зоне; предпочтительно способ применим к данным магнитных измерений в зоне с магнитным наклонением в диапазоне ±30°; и более предпочтительно способ применим к данным магнитных измерений в зоне с магнитным наклонением в диапазоне ±20°.
В соответствии со способами, предусмотренными настоящим изобретением, могут быть получены сведения о гравитационно-магнитных аномалиях на разных глубинах в зоне исследований, сведения о распределении структурных строений, представляющих геологические структуры на разных глубинах, и сведения об изменении напряженности магнитного поля и интенсивности плотности структурных строений на разных глубинах, достигается идентификация и качественная интерпретация геологических структур, обусловливающих образование рудных месторождений, в соответствии с априорными знаниями о зоне исследований определяются типы потенциальных рудных месторождений и свойства структур, обусловливающих образование рудных месторождений, выполняется отбор разных типов структурных строений, и, таким образом, определяется местоположение объектов разведки — месторождений металлических руд.
Способ автоматического получения структурных строений из данных магнитных измерений в низкоширотных зонах в настоящем изобретении решает существующую на предыдущем уровне техники проблему, связанную с трудностью точного получения структурных сведений путем использования данных магнитных измерений в низкоширотных зонах. По сравнению с существующими изображениями с координатной сеткой или контурными картами карты структурного строения, полученные в соответствии с настоящим изобретением, интуитивно представляют такие структурные сведения, как глубина врезания, глубина залегания, отношение первичных и вторичных структур, отношение врезания и напряженность магнитного поля, имеющие большое значение для геологической интерпретации и разведки полезных ископаемых.
Способ автоматического получения магнитного структурного строения в низкоширотной зоне, предусмотренный настоящим изобретением, также является применимым для автоматического получения магнитных структурных строений в высокоширотных зонах, а также является применимым для автоматического получения структурных строений гравитационных потенциальных полей.
Путем применения способов в соответствии с настоящим изобретением расширяется диапазон зон со структурными строениями, которые анализируют и получают при помощи магнитного способа, повышается точность автоматического получения структурных строений, может быть получена идентификация и качественная интерпретация геологических структур, обусловливающих образование рудных месторождений, в соответствии с априорными знаниями о зоне исследований определяют типы потенциальных рудных месторождений и свойства структур, обусловливающих образование рудных месторождений, выполняется отбор разных типов структурных строений, и, таким образом, достигается определение местоположений объектов разведки — месторождений металлических руд.
Описание графических материалов
Дополнительные признаки настоящего изобретения станут очевидными из нижеследующего описания примерных вариантов осуществления с отсылкой к графическим материалам. Описываемые в настоящем документе графические материалы используются для обеспечения более глубокого понимания изобретения и составляют часть настоящей заявки. Схематические варианты осуществления настоящего изобретения и их описание используются для разъяснения изобретения, но не для ненадлежащего ограничения изобретения. В графических материалах:
На фиг. 1 проиллюстрирована блок-схема способа автоматического получения структурного строения из данных потенциального поля в соответствии с Вариантом осуществления 1 настоящего изобретения.
На фиг. 2 проиллюстрирована карта структурного строения однопиксельной ширины в соответствии с Примером 1 настоящего изобретения.
На фиг. 3 проиллюстрирована составная карта структурного строения в соответствии с Примером 1 настоящего изобретения.
На фиг. 4 проиллюстрирована карта структурного строения с отображением интенсивности в соответствии с Примером 1 настоящего изобретения.
На фиг. 5 проиллюстрирована составная карта структурного строения с отображением интенсивности в соответствии с Примером 1 настоящего изобретения.
На фиг. 6 проиллюстрирована блок-схема способа автоматического выделения магнитного структурного строения в низкоширотной зоне в соответствии с Вариантом осуществления 2 настоящего изобретения.
На фиг. 7 проиллюстрирована карта структурного строения однопиксельной ширины в соответствии с Примером 2 настоящего изобретения.
На фиг. 8 проиллюстрирована составная карта структурного строения в соответствии с Примером 2 настоящего изобретения.
На фиг. 9 проиллюстрирована карта структурного строения, отражающая глубину залегания, в соответствии с Примером 2 настоящего изобретения.
На фиг. 10 проиллюстрирована составная карта структурного строения с отражением глубины залегания в соответствии с Примером 2 настоящего изобретения.
На фиг. 11 проиллюстрирована карта структурного строения, отражающая напряженность магнитного поля, в соответствии с Примером 2 настоящего изобретения.
На фиг. 12 проиллюстрирована составная карта структурного строения с отражением напряженности магнитного поля в соответствии с Примером 2 настоящего изобретения.
Описание вариантов осуществления
Для лучшего понимания технического решения настоящего изобретения способы реализации настоящего изобретения будут дополнительно описаны ниже с отсылкой к графическим материалам в сочетании с вариантами осуществления, которые, однако, не следует рассматривать как ограничения настоящего изобретения.
Вариант осуществления 1
На фиг. 1 проиллюстрирована блок-схема способа автоматического получения структурного строения из данных потенциального поля в соответствии с Вариантом осуществления 1 настоящего изобретения. Как проиллюстрировано на фиг. 1, способ включает следующие этапы:
На этапе 101 выполняется предварительная обработка измеренных данных гравитационного потенциального поля или данных магнитного потенциального поля.
В отношении магнитных данных, с целью получения приведенных к полюсу магнитных аномалий, выполняется вычисление приведения к полюсу, или, с целью получения псевдогравитационных аномалий, применяется псевдогравитационное преобразование.
Предварительная обработка выполняется в отношении гравитационных данных с целью получения гравитационных аномалий Буге.
На этапе 102 к предварительно обработанным данным гравитационного потенциального поля или данным магнитного потенциального поля с целью получения краев на множестве уровней применяется способ многонаправленного и многоуровневого обнаружения краев потенциального поля.
Способ многонаправленного и многоуровневого обнаружения краев потенциального поля включает этапы продолжения вверх предварительно обработанных данных гравитационного потенциального поля или данных магнитного потенциального поля на множество предварительно определенных высот с целью получения множества данных потенциального поля на соответствующих уровнях и применения многонаправленного обнаружения краев, соответственно, к данным потенциального поля на каждом уровне с целью получения края потенциального поля на указанном уровне.
Способ применения многонаправленного обнаружения краев соответственно к данным потенциального поля на каждом уровне включает следующие этапы:
Гравитационная аномалия или магнитная аномалия в точке (x, y) поверхности Земли на нулевой высоте определяется как f0(x, y).
Пусть уровень s=z/z0, и z>z0, где z0 представляет высоту измерения, а z представляет продолженную вверх высоту.
Сглаживающая функция для уровня s определяется как:
Figure 00000001
,
где
Figure 00000002
, k(x, y, z) — функция Грина.
Вейвлет-функция в направлении a определяется как:
Figure 00000003
,
где D представляет производную первого порядка.
Для уровня s и местоположения (x, y) вейвлет-преобразование гравитационной аномалии или магнитной аномалии f0(x, y) в направлении a определяется как:
Figure 00000004
где * представляет операцию свертки,
в соответствии с формулой продолжения потенциального поля вверх, известно, что:
Figure 00000005
, и
fz(x, y) — гравитационная аномалия или магнитная аномалия посредством продолжения вверх f0(x, y) от нулевой высоты до высоты z=sz0, и полученная путем продолжения вверх гравитационной аномалии или магнитной аномалии, fz0(x, y), измеренной на высоте измерения z0, для высоты z-z0,
отсюда
Figure 00000006
Figure 00000007
.
Аналогично, для уровня s и местоположения (x, y) вейвлет-преобразование гравитационной аномалии или магнитной аномалии f0(x, y) в направлении a+
Figure 00000008
определяется как:
Figure 00000009
.
Градиент можно использовать для представления двумерного вейвлет-преобразования f0(x, y) как:
Figure 00000010
где
Figure 00000011
— двумерный градиент.
Приведенное выше уравнение устанавливает отношение между горизонтальным градиентом гравитационной аномалии или магнитной аномалии fz(x, y) на любой высоте z (z>z0) и двумерным направленным вейвлет-преобразованием гравитационной аномалии или магнитной аномалии f0(x, y) на нулевой высоте.
Для местоположения (x, y), уровня s и направления a двумерное вейвлет-преобразование
Figure 00000012
f0(x, y) прямо пропорционально градиенту
Figure 00000013
fz(x, y), и двумерное вейвлет-преобразование
Figure 00000012
f0(x, y) можно представить градиентом
Figure 00000013
fz(x, y).
Для высоты z после продолжения вверх модуль градиента
Figure 00000014
определяется как:
Figure 00000015
,
где соответствующий аргумент градиента по горизонтальной оси имеет вид:
Figure 00000016
, и
точка, в которой модуль M
Figure 00000017
представляет собой локальный максимум вдоль направления
Figure 00000018
аргумента, определяется как краевая точка.
Для каждой высоты z после продолжения вверх для вычисления выбираются разные направления a, и могут быть выделены сведения о краях в разных направлениях. Для достижения полного покрытия двумерной плоскости, значение a направления может быть принято равным kπ/(2n-1), где k=0, 1, 2…(2n-1), и n — целое число, большее или равное 2. В каждом направлении, путем соединения точек с локальным максимумом модуля градиента в направлении, перпендикулярном градиенту, формируется кривая с целью создания края. Для той же высоты после продолжения вверх вычисляются края в направлениях, полностью охватывающих двумерную плоскость, и для получения края на уровне берется объединение всех вычисленных краев.
Путем выполнения вышеописанного вычисления в отношении данных потенциального поля на каждом уровне может быть достигнуто многонаправленное обнаружение краев данных потенциального поля на каждом уровне.
На этапе 103 полученные края на каждом уровне подвергают утончению до однопиксельной ширины с использованием алгоритма определения морфологического скелета с целью получения карты строения на каждом уровне.
Для полученных изображений краев края подвергают утончению до однопиксельной ширины с использованием алгоритма утончения строения, описанного в публикации «Thinning Methodologies-A Comprehensive Survey», Lam, L., Seong-Whan Lee and Ching Y. Suen, IEEE Transactions on Pattern Analysis and Machine Intelligence, т. 14, №. 9, сентябрь 1992.
За счет утончения полученных краев фактический ареал, соответствующий краям и пересечениям краев в разных направлениях, значительно сокращается по сравнению с предыдущим уровнем техники, так что полученное структурное строение потенциального поля является более близким к признакам фактической структуры, идентифицированной с использованием геологического картирования. С одной стороны, сведения карты являются четкими, и ее удобочитаемость улучшается. С другой стороны, также облегчается геологическая интерпретация.
На этапе 104 полученные структурные строения на всех уровнях накладываются одно на другое с целью создания составной карты структурного строения.
Края на всех уровнях, выделенные из данных потенциального поля после продолжения вверх на разные высоты, накладывают друг на друга с целью получения составной карты структурного строения потенциального поля зоны исследований.
Приведение разных высот после продолжения потенциального поля вверх в соответствие разным глубинам объектов-источников (см. случай А для продолжения вверх в качестве стандартного разделительного фильтра для карт потенциальных полей, Jacobsen, B.H., Geophysics, т. 52, №.8, 1987) позволяет получать структуры, соответствующие разным глубинам Земли. Строения на соответствующих разных уровнях, или на разных глубинах, представляют при помощи постепенного изменения цветов, и структурные строения на всех уровнях, представляемые постепенно изменяющимися цветами, накладывают одно на другое с целью выделения сведений о структурах на разных глубинах. Путем такого наложения получают составную карту структурного строения, отражающую сведения на разных глубинах.
На этапе 105 модуль градиента каждой краевой точки в строении на каждом уровне используют для отражения напряженности края.
Для разных направлений a модуль
Figure 00000027
градиента каждой краевой точки не изменяется и представляет собой постоянное значение. Модуль градиента краевой точки принимается в качестве значения структурной интенсивности в данной краевой точке с целью установления карт структурного строения с отображением интенсивности на разных уровнях, или разных глубинах.
На этапе 106 карты структурного строения с отображением интенсивности на множестве уровней накладывают одну на другую с целью создания составной карты структурного строения с отображением интенсивности.
Путем представления значений интенсивности краев на уровнях, или глубинах, разными постепенно изменяющимися цветами и наложения могут быть выделены сведения об изменении напряженности магнитного поля или плотности.
Пример 1
Техническое решение настоящего изобретения будет разъяснено на примере автоматического получения структурного строения потенциального поля из аэромагнитных данных для западного Юньнаня.
В первую очередь выполняют приведение к полюсу в отношении аэромагнитных данных в блоках, обработанные данные соединяют с формированием сеточного файла потенциального поля, и шаг сетки составляет 500 м. Аэромагнитные данные в данном примере представляют собой исторические данные, многократно измеренные и собранные в прошлом, высота измерений находится в диапазоне 800—1200 м.
Затем в отношении объединенного сеточного файла потенциального поля соответственно выполняют обработку продолжения вверх с получением данных потенциального поля после продолжения вверх, при этом высоты продолжения вверх составляют, соответственно, 1000 м, 1500 м, 2000 м, 2500 м, 3000 м, 4000 м, 5000 м, 10000 м, 15000 м, 20000 м, 25000 м и 30000 м, и к каждой продолженной вверх высоте применяют обнаружение краев в тридцати двух направлениях, при этом каждое направление a принимает значение kπ/(2n-1), где k=0, 1, 2,…, (2n-1), и n=5, , и, таким образом, получают края на каждом уровне.
Затем, с целью получения карты структурного строения на каждом уровне, в отношении полученных в результате краев на каждом уровне выполняют обработку утончения с использованием алгоритма определения скелета.
На фиг. 2 проиллюстрировано структурное строение однопиксельной ширины, полученное на уровне, соответствующем высоте, продолженной вверх на 5000 м, с краями, подвергнутыми утончению с использованием алгоритма определения скелета. Видно, что структурное строение потенциального поля однопиксельной ширины, автоматически полученное при помощи способа настоящего изобретения, является более близким к признакам фактической структуры, идентифицированной с использованием геологического картирования, и, таким образом, облегчается его геологическая интерпретация. В дополнение, сведения карты являются более четкими, и анализ наложения строений на разных уровнях облегчается.
На фиг. 3 проиллюстрирована составная карта структурного строения, сформированная путем представления при помощи постепенно изменяющихся цветов и наложения строений на всех уровнях, полученных путем обнаружения краев в тридцати двух направлениях в отношении данных потенциального поля посредством продолжения вверх для высот, соответственно, 1000 м, 1500 м, 2000 м, 2500 м, 3000 м, 4000 м, 5000 м, 10000 м, 15000 м, 20000 м, 25000 м и 30000 м. Составная карта структурного строения, соответствующая высотам, полученным посредством продолжения вверх до соответствующих глубин источников, может использоваться для представления сведений о структурных строениях на разных глубинах зоны исследований. Путем использования цветов, постепенно изменяющихся от серо-белого до черного, для представления продолжения вверх высот от низких до высоких, или глубин от малых до больших, составная карта структурного строения отражает сведения о структурах на разных глубинах.
На фиг. 4 проиллюстрирована структурная интенсивность, отражаемая модулем градиента в краевых точках строения, полученного после продолжения вверх для высоты 5000 м. Цвета, изменяющиеся от серо-белого до черного, отражают то, что значения структурной интенсивности, представленные модулем градиента, увеличиваются постепенно. Данная карта отражает изменение структурной интенсивности, соответствующее одному уровню.
На фиг. 5 проиллюстрирована составная карта структурного строения с отображением интенсивности, полученная путем наложения краевых интенсивностей, полученных путем продолжения данных потенциального поля вверх для высот 1000 м, 1500 м, 2000 м, 2500 м, 3000 м, 4000 м, 5000 м, 10000 м, 15000 м, 20000 м, 25000 м и 30000 м. Как видно на данной карте, на разных уровнях основные структурные пояса представляют высокую структурную интенсивность, и соответствующие магнитные аномалии изменяются явным образом, что указывает на то, что все основные структурные пояса представляют собой пояса резких изменений магнитных аномалий, при этом соответствующие глубины являются очень большими, и отражает то, что структурные пояса обусловливают процесс магматического образования руд на большой глубине.
Кроме того, на фиг. 3 и фиг. 5, карты структурного строения, содержащие сведения о глубине и интенсивности, могут в четко отражать первичные структуры на большой глубине и с большой протяженностью в зоне, вторичные структуры на малой глубине и с малой протяженностью, а также обоюдные отношения врезания. Следовательно, карта структурного строения, полученная способом в соответствии с настоящим изобретением, может помочь специалисту в данной области техники распознать структурное строение зоны изысканий.
При сравнении фиг. 3 и фиг. 5 с полученной наземным картированием геологической структурой данной зоны пространственное местоположение и диапазон очень хорошо согласуются, что указывает на точность и эффективность способа автоматического получения в соответствии с настоящим изобретением. В дополнение, в сравнении с полученной наземным картированием геологической структурой данной зоны, структурные сведения о трехмерной протяженности, интенсивности и структуре содержатся на фиг. 3 и фиг. 5 настоящего изобретения, что может способствовать идентификации скрытых структурных поясов.
Зона вблизи структурных поясов на больших глубинах и с большой структурной интенсивностью, пересекающиеся участки структур разных направлений, а также изогнутые и криволинейные участки структурных поясов являются важными местоположениями, где могут быть разведаны потенциальные месторождения металлических руд. С помощью способа в соответствии с настоящим изобретением карты структурного строения, содержащие сведения о глубине и интенсивности, могут быть оперативно и точно получены из аэромагнитных данных и гравитационных данных и помочь геологоразведчикам точно и оперативно исследовать потенциальные месторождения металлических руд.
Вариант осуществления 2
Ниже на примере данных магнитных измерений будет конкретно описана блок-схема способа автоматического получения магнитного структурного строения низкоширотной зоны в соответствии с Вариантом осуществления 2 настоящего изобретения. Специалисту в данной области техники будет понятно, что способ, предусматриваемый настоящим изобретением, не ограничен применением для автоматического получения структурных строений из данных магнитной разведки, и также является применимым для автоматического получения структурных строений из данных гравитационных измерений.
На фиг. 6 проиллюстрирована блок-схема способа автоматического получения магнитного структурного строения из низкоширотной зоны в соответствии с настоящим изобретением. Как проиллюстрировано на фиг. 6, способ включает следующие этапы:
На этапе 601 выполняют предварительную обработку в отношении измеренных значений, полученных посредством магнитной разведки, и для получения данных аномалий полной напряженности магнитного поля (TMI) выполняют коррекцию по нормальному геомагнитному полю (Международному эталону геомагнитного поля, IGRF).
В отношении данных аномалий TMI выполняют нанесение на координатную сетку, и шаг координатной сетки полагают равным от 1/8 до 1/4 расстояния между линиями съемки или минимальному расстоянию между точками.
На этапе 602 нанесенные на координатную сетку данные T аномалий TMI продолжают вверх на множество предварительно определенных высот с получением множества нанесенных на координатную сетку данных Th аномалий TMI на соответствующих уровнях, где h представляет высоту после продолжения вверх.
На этапе 603 наклонную производную TDRh данных аномалий TMI на каждом уровне вычисляют с использованием нанесенных на координатную сетку данных Th аномалий TMI на соответствующем уровне в соответствии со следующей формулой:
Figure 00000028
,
где VDRh и THDRh, соответственно, представляют вертикальную производную первого порядка и полную горизонтальную производную нанесенных на координатную сетку данных Th аномалий TMI:
Figure 00000029
,
Figure 00000030
.
На этапе 604 к наклонной производной TDRh на каждом уровне, соответственно, применяют многонаправленное обнаружение краев на основе горизонтального градиента.
Направленные производные наклонной производной TDRh в направлениях a и a+
Figure 00000008
, соответственно, определяют как:
Figure 00000019
Figure 00000020
где D представляет производную первого порядка.
Горизонтальный градиент наклонной производной TDRh для высоты h и направления a представляют как:
Figure 00000021
где
Figure 00000011
— горизонтальный градиент.
Модуль горизонтального градиента
Figure 00000022
определяют как:
Figure 00000023
где соответствующий аргумент градиента имеет вид:
Figure 00000024
,
полная горизонтальная производная наклонных производных (сокращенно: TDR_THDR) данных аномалий TMI имеет вид:
Figure 00000031
, и
амплитуда модуля
Figure 00000032
не связана с величиной магнитного наклонения.
Таким образом, точки, в которых модуль
Figure 00000025
представляет собой локальный максимум вдоль направления аргумента
Figure 00000026
, представляют собой краевые точки объектов-источников магнитных аномалий на высоте h в направлении a.
В каждом направлении a точки с локальным максимумом модуля градиента соединяют в направлении, перпендикулярном градиенту, для получения кривой с целью создания края.
Края на множестве разных направлений a вычисляют на одной и той же высоте h, и, для получения краев объектов-источников магнитных аномалий на соответствующих уровнях, берут объединение всех созданных краев.
Для достижения полного покрытия двумерной плоскости, разные значения a направления соответственно полагают равными kπ/(2n-1), где k=0, 1, 2…(2n-1), и n — целое число, большее или равное 2.
Выполнение вычисления TDR в отношении нанесенных на координатную сетку данных аномалий TMI может вносить шум в результаты вычисления, а последующее многонаправленное вычисление краев является очень чувствительным к шуму. Для предотвращения влияния шума на многонаправленное обнаружение краев, в отношении данных TDRh с высоким уровнем шума перед выполнением многонаправленного обнаружения краев предпочтительно выполняют обработку, связанную с подавлением шума, такую как гауссова фильтрация.
Поскольку полные горизонтальные производные наклонных производных данных аномалий TMI не зависят от магнитного наклонения, результаты вычисления не связаны с величиной магнитного наклонения, магнитное наклонение не оказывает влияния на края объектов-источников магнитных аномалий, полученные на вышеописанных этапах, т. е. не оказывают влияния широтные местоположения зоны, подлежащей анализу, и структурное строение низкоширотной зоны может быть представлено точно.
На этапе 605 вычисленные края на каждом уровне соответственно подвергают утончению до однопиксельной ширины с использованием алгоритма определения морфологического скелета с целью получения карт структурного строения на множестве уровней.
Для вычисленных изображений краев края подвергают уточнению до однопиксельной ширины с использованием алгоритма для утончения строения, описанного в публикации «Thinning Methodologies-A Comprehensive Survey (Методики утончения, Исчерпывающий обзор)», Lam, L., Seong-Whan Lee and Ching Y. Suen, IEEE Transactions on Pattern Analysis and Machine Intelligence, т. 14, №. 9, сентябрь 1992.
Посредством утончения полученных краев фактический ареал, соответствующий краям и пересечениям краев в разных направлениях, значительно сокращается по сравнению с предыдущим уровнем техники, так что полученное структурное строение является более близким к признакам фактической структуры, идентифицированной с использованием геологического картирования. С одной стороны, сведения карты являются четкими, и ее удобочитаемость улучшается. С другой стороны, также облегчается геологическая интерпретация.
На этапе 606 полученные карты структурного строения на всех уровнях накладывают одну на другую с целью создания составной карты структурного строения.
Полученные карты структурного строения на всех уровнях накладывают одну на другую с целью создания составной карты структурного строения, отражающей сведения на разных глубинах, и боковое отклонение краев на разных уровнях на карте отражает сведения об ориентации структурных строений.
Путем приведения краев, полученных из нанесенных на координатную сетку данных аномалий TMI после продолжения вверх для разных высот, в соответствие структурам на разных глубинах, где глубина составляет половину высоты после продолжения вверх (см. случай А в публикации «Upward Continuation as a Standard Separation Filter for Potential-Field Maps (Продолжение вверх как стандартный фильтр разделения для карт потенциальных полей)», Jacobsen, B.H., Geophysics, т.52, №.8, 1987), могут быть получены карты структурного строения, представляющие разные глубины врезания. Путем наложения одна на другую полученных карт структурного строения на всех глубинах можно получить составную карту структурного строения, отражающую сведения на разных глубинах врезания.
На этапе 607 модуль горизонтального градиента TDRh в каждой краевой точке на карте структурного строения на каждом уровне используют для представления структурной глубины залегания с целью получения карт структурного строения, представляющих структурные глубины залегания на множестве уровней.
Для разных направлений a модуль
Figure 00000033
горизонтального градиента в каждой краевой точке не изменяется и представляет собой постоянное значение.
Из-за свойства значения арктангенса тригонометрической функции, вне зависимости от того, насколько велики амплитуды VDRh и THDRh, амплитуда TDRh ограничена промежутком между -π/2 и +π/2. Поэтому модуль
Figure 00000034
горизонтального градиента в каждой краевой точке имеет малую корреляцию с амплитудой аномалий TMI, значение отражает глубины залегания объектов-источников, величина значения обратно пропорциональна глубине залегания, и чем больше значение, тем меньше глубина залегания (см. публикацию «New Insights into Magnetic Derivatives for Structural Mapping (Новые представления о магнитных производных для структурного картирования)», Bruno Verduzco и др., The Leading Edge, т. 23, №. 2, 2004).
Структурная глубина залегания, отражаемая модулем горизонтального градиента наклонной производной TDRh, представляет толщину покровного слоя над структурным поясом, и разные глубины, соответствующие продолжению вверх на множество предварительно определенных высот, отражают глубину врезания структурного пояса, которая представляет глубину протяженности структурного пояса вниз.
Модуль
Figure 00000035
градиента каждой краевой точки используют для представления связанной с краем глубины залегания краевой точки на уровне с целью установления карт структурного строения, отражающих связанную со структурой глубину залегания на разных уровнях, или глубинах. Изменения в значениях модуля градиента одного и того же структурного пояса отражают глубины залегания разных участков данного структурного пояса.
На этапе 608 карты структурного строения, представляющие глубины залегания на множестве уровней, накладывают одну на другую с целью получения составной карты структурного строения с отражением глубины залегания.
Значения модуля горизонтального градиента производной TDRh краевых точек на уровнях, или глубинах, накладывают одно на другое и представляют разными постепенно изменяющимися цветами с целью выделения сведений об изменениях глубины залегания структурных строений разных диапазонов глубин врезания.
На этапе 609 для нанесенных на координатную сетку данных Th аномалий TMI на каждой высоте вычисляют трехмерный аналитический сигнал ASh в соответствии с формулой:
Figure 00000036
Амплитуда получаемого в результате аналитического сигнала ASh обладает сильной корреляцией с амплитудой аномалий TMI, однако она не зависит от величины магнитного наклонения и может быть использована для указания местоположений и напряженности магнитного поля объектов-источников магнитных аномалий.
На этапе 610 значение ASh каждой краевой точки строения на каждом уровне используют для указания напряженности магнитного поля структуры.
Без учета глубин залегания объектов-источников магнитных аномалий амплитуда трехмерного аналитического сигнала ASh указывает напряженность магнитного поля источников магнитных аномалий, а структуры с сильными магнитными свойствами обычно имеют тесную связь с образованием руд.
На этапе 611 карты структурного строения с отображением напряженности магнитного поля на множестве уровней накладывают одну на другую с целью создания составной карты структурного строения с отображением напряженности магнитного поля.
Значения ASh краевых точек на уровнях, или глубинах, представляют постепенно изменяющимися цветами и накладывают одно на другое с целью выделения сведений об изменениях напряженности магнитного поля структурных строений в диапазоне разных глубин.
Пример 2
Техническое решение настоящего изобретения будет разъяснено ниже на примере автоматического получения магнитного структурного строения из высокоточных аэромагнитных данных Гвианского плоскогорья в низкоширотной зоне Южной Америки.
Масштаб съемок аэромагнитных данных в данном примере составляет 1:25000, а высота измерений находится в диапазоне 70—120 м.
В первую очередь выполняют предварительную обработку значений измерений, полученных посредством аэромагнитной съемки зоны исследований, выполняют коррекцию по нормальному геомагнитному полю (Международному эталону геомагнитного поля, IGRF) для получения данных аномалий полной напряженности магнитного поля (TMI) и выполняют нанесение на координатную сетку данных аномалий TMI с шагом координатной сетки приблизительно 10 м.
Затем соответственно выполняют обработку продолжения вверх в отношении нанесенных на координатную сетку данных аномалий TMI с целью получения нанесенных на координатную сетку данных Th аномалий TMI на множестве уровней, при этом высоты продолжения вверх составляют, соответственно, 100 м, 200 м, 300 м, 400 м и 500 м.
В отношении нанесенных на координатную сетку данных Th аномалий TMI на каждом уровне выполняют вычисление TDR.
Для каждой высоты продолжения вверх в отношении вычисленных TDR на каждой высоте после продолжения вверх выполняют многонаправленное обнаружение краев на основе горизонтального градиента. Для получения краев на каждом уровне обнаружение краев выполняют в шестидесяти четырех направлениях для каждого результата вычисления TDR, и разные значения a направления соответственно полагают равными kπ/(2n-1), где k=0, 1, 2…(2n-1), и n=6.
Затем в отношении получаемых в результате краев на каждом уровне выполняют обработку утончения с использованием алгоритма определения скелета с целью получения карты структурного строения на каждом уровне.
На фиг. 7 проиллюстрировано магнитное структурное строение однопиксельной ширины на соответствующем уровне, полученное путем извлечения из нанесенных на координатную сетку данных аномалий TMI после продолжения вверх на высоту 300 м и утончения краев с использованием алгоритма определения скелета. Видно, что магнитное структурное строение однопиксельной ширины, автоматически полученное при помощи способа настоящего изобретения, более близко к признакам фактической структуры, идентифицированной с использованием геологического картирования, и, таким образом, облегчается геологическая интерпретация. В дополнение, сведения карты являются более четкими, и анализ наложения структурных строений на разных уровнях облегчается.
На фиг. 8 проиллюстрирована составная карта структурного строения, сформированная путем представления постепенно изменяющимися цветами и наложения строений на всех уровнях, полученных при помощи обнаружения краев в шестидесяти четырех направлениях в отношении нанесенных на координатную сетку данных аномалий TMI после продолжения вверх для высот 100 м, 200 м, 300 м, 400 м и 500 м. За счет приведения высот, полученных посредством продолжения вверх, в соответствие глубинам объектов-источников, составная карта структурного строения может быть использована для представления сведений о структурных строениях на разных глубинах зоны исследований. За счет использовании цветов, постепенно изменяющихся от серо-белого до черного, для представления продолженных вверх высот от низких до высоких или глубин от мелких до глубоких, составная карта структурного строения отражает сведения о структурах на разных глубинах врезания.
Затем, для представления структурной глубины залегания, в каждой краевой точке строения на каждом уровне вычисляют модуль градиента TDRh. На фиг. 9 проиллюстрированы структурные глубины залегания, представленные модулем градиента в краевых точках строения, продолженного вверх на высоту 300 м, при этом цвета, постепенно изменяющиеся от серо-белого до черного, отражают то, что структурные глубины залегания, представленные модулем горизонтального градиента, постепенно уменьшаются. Данная карта отражает изменения структурной глубины залегания на соответствующем уровне.
Затем карты структурного строения на множестве уровней, которые представляют относительные глубины залегания, накладывают одну на другую с целью создания составной карты структурного строения с отражением глубины залегания. На фиг. 10 проиллюстрирована составная карта структурного строения с отражением глубин залегания, полученная путем наложения полученных карт структурного строения, представляющих относительные глубины залегания после продолжения вверх нанесенных на координатную сетку данных аномалий TMI для высот 100 м, 200 м, 300 м, 400 м и 500 м. Из карты видно, что на разных уровнях глубины залегания основных структурных поясов заметно изменяются, и обычно структуры с большими глубинами залегания, т. е. с более толстыми покровными слоями, имеют меньшие значения аномалий TMI, и структуры представляют собой скрытые структуры, и следовательно данный способ облегчает идентификацию скрытых структур.
Затем для нанесенных на координатную сетку данных аномалий TMI на каждом уровне вычисляют трехмерный аналитический сигнал ASh. Значение получаемого в результате аналитического сигнала ASh не зависит от величины магнитного наклонения, и его можно использовать для указания местоположений и напряженности магнитного поля объектов-источников магнитных аномалий.
В данном примере значение ASh каждой краевой точки строения на каждом уровне используется для указания напряженности магнитного поля структуры. На фиг. 11 проиллюстрирована карта структурного строения, на которой значения ASh указывают напряженность магнитного поля краевых точек структуры, продолженной вверх на высоту 300 м, при этом цвета, постепенно изменяющиеся от серо-белого до черного, представляют то, что магнитные свойства структурных поясов, указываемые значениями аналитического сигнала ASh, постепенно усиливаются.
Затем карты структурного строения с отображением напряженности магнитного поля на множестве уровней накладывают одну на другую с целью создания составной карты структурного строения с отображением напряженности магнитного поля. Как проиллюстрировано на фиг. 12, значения ASh краев на всех уровнях, или глубинах, накладывают одно на другое и представляют цветами, постепенно изменяющимися от серо-белого до черного, с целью выделения сведений об изменениях магнитных свойств структурных строений на разных глубинах.
Без учета глубин залегания объектов-источников амплитуда трехмерного аналитического сигнала ASh указывает напряженность магнитного поля источников магнитных аномалий, а структуры с сильными магнитными свойствами обычно имеют тесную связь с образованием руд.
В дополнение, при рассмотрении фиг. 8 и фиг. 10, карты структурного строения, содержащие сведения о глубине врезания и напряженности, могут четко отражать первичные структуры на больших глубинах и с длительной протяженностью в зоне, вторичные структуры малых глубин и краткой протяженности, глубины залегания структур, а также взаимные отношения врезания. Следовательно, карта структурного строения, полученная способом в соответствии с настоящим изобретением, может помочь специалисту в данной области техники распознать структурное строение зоны изысканий.
Важными местоположениями, в которых могут быть разведаны потенциальные месторождения металлических руд, являются: зона вблизи структурных поясов с сильными магнитными свойствами на больших глубинах, участки пересечения структур разных направлений, а также изогнутые и криволинейные участки структурных поясов. С использованием способа в соответствии с настоящим изобретением можно оперативно и точно получать из данных магнитных измерений карты структурного строения, содержащие сведения о глубине врезания, глубине залегания, напряженности магнитного поля, отношении первичных и вторичных структур, отношении врезания и т. п., что может содействовать геологоразведчикам в проведении точной и быстрой разведки потенциальных месторождений металлических руд.
Настоящее изобретение было подробно описано выше в комбинации с предпочтительными вариантами осуществления. Однако настоящее изобретение ими не ограничивается. Специалисты в данной области техники могут осуществлять разнообразные модификации в соответствии с принципами настоящего изобретения. Поэтому модификации, осуществленные в соответствии с принципами настоящего изобретения, следует понимать как включаемые в объем правовой защиты настоящего изобретения.

Claims (75)

1. Способ автоматического получения структурного строения из данных потенциального поля, включающий этапы:
предварительной обработки данных потенциального поля из зоны исследований;
продолжения вверх предварительно обработанных данных потенциального поля на множество предварительно определенных высот с целью получения множества данных потенциального поля на соответствующих уровнях;
применения многонаправленного обнаружения краев соответственно к данным потенциального поля на множестве уровней с целью получения множества краев потенциального поля на соответствующих уровнях; и
утончения полученных краев потенциального поля на уровнях соответственно до однопиксельной ширины с использованием алгоритма определения морфологического скелета с целью получения множества карт структурного строения на соответствующих уровнях.
2. Способ по п. 1, отличающийся тем, что способ дополнительно включает наложение множества карт структурного строения на соответствующих уровнях с целью создания составной карты структурного строения.
3. Способ по п. 1, отличающийся тем, что способ дополнительно включает применение модуля градиента каждой краевой точки на карте структурного строения на каждом уровне в качестве значения интенсивности краевой точки на карте структурного строения на уровне с целью получения множества карт структурного строения с отображением интенсивности на соответствующих уровнях.
4. Способ по п. 3, отличающийся тем, что способ дополнительно включает наложение множества карт структурного строения с отображением интенсивности на соответствующих уровнях с целью создания составной карты структурного строения с отображением интенсивности.
5. Способ по п. 1, отличающийся тем, что данные потенциального поля представляют собой данные гравитационного потенциального поля или данные магнитного потенциального поля, и этап предварительной обработки дополнительно включает:
предварительную обработку гравитационных данных с целью получения гравитационных аномалий Буге; или
приведение магнитных данных к полюсу с целью получения приведенных к полюсу магнитных аномалий или применение к магнитным данным псевдогравитационного преобразования с целью получения псевдогравитационных аномалий.
6. Способ по п. 1, отличающийся тем, что:
этап применения многонаправленного обнаружения краев к данным потенциального поля на каждом уровне включает следующие этапы:
в предположении, что уровень s=z/z0, z>z0, z0 — это высота измерения, и z — высота, увеличивающаяся вверх, — определение гравитационной аномалии или магнитной аномалии в местоположении (x, y) на нулевой высоте как f0(x, y),
определение сглаживающей функции на уровне s как:
Figure 00000037
,
где
Figure 00000038
, k(x, y, z) — функция Грина,
определение вейвлет-функции в направлении a как:
Figure 00000039
,
где D представляет производную первого порядка;
для уровня s и местоположения (x, y) — определение вейвлет-преобразования гравитационной аномалии или магнитной аномалии f0(x, y) в направлении a как:
Figure 00000040
где * представляет операцию свертки,
при этом известно, что
Figure 00000041
в соответствии с формулой продолжения вверх потенциального поля, и
fz(x, y) — гравитационная аномалия или магнитная аномалия посредством продолжения вверх f0(x, y) от нулевой высоты до высоты z=sz0, полученная путем продолжения вверх гравитационной аномалии или магнитной аномалии, fz0(x, y), измеренной на высоте z0 измерения, на высоту z-z0,
тогда
Figure 00000042
Figure 00000043
;
кроме того, для уровня s и местоположения (x, y), — определение вейвлет-преобразования гравитационной аномалии или магнитной аномалии f0(x, y) в направлении a+
Figure 00000044
как:
Figure 00000045
,
тогда двумерное направленное вейвлет-преобразование f0(x, y) записывают как градиент:
Figure 00000046
где
Figure 00000047
обозначает двумерный градиент,
для высоты z — определение модуля градиента
Figure 00000048
как:
Figure 00000049
,
где соответствующий аргумент градиента вдоль горизонтальной оси имеет вид:
Figure 00000050
, и
определение точки, в которой модуль M
Figure 00000051
представляет собой локальный максимум вдоль направления
Figure 00000052
аргумента, как краевой точки,
в каждом направлении a формирование кривой путем соединения точек с локальным максимумом модуля градиента в направлении, перпендикулярном градиенту, с целью создания края и
для той же высоты, — вычисление краев на множестве разных направлений a и принятие объединения краев в качестве краев потенциального поля на соответствующих уровнях.
7. Способ по п. 6, отличающийся тем, что этап вычисления краев на множестве разных направлений a для одной высоты дополнительно включает: принятие значения a направления равным kπ/(2n-1), где k=0, 1, 2…(2n-1), и n — целое число, большее или равное 2, для полного покрытия двумерной плоскости.
8. Способ по п. 6, отличающийся тем, что дополнительно включает приведение краев, полученных из потенциальных полей, продолженных вверх на предварительно определенные высоты, в соответствие структурам на разных глубинах и наложение карт структурного строения на уровнях одна на другую с целью получения составной карты структурного строения, отражающей сведения на разных глубинах.
9. Способ по п.2 или 8, отличающийся тем, что для представления высот с целью формирования составной карты структурного строения используют постепенно изменяющиеся цвета.
10. Способ по п.4, отличающийся тем, что величины значений интенсивности представляют с использованием постепенно изменяющихся цветов с целью формирования составной карты структурного строения с отображением интенсивности.
11. Способ автоматического получения структурного строения, включающий следующие этапы:
предварительную обработку магнитных данных или данных гравитационных изысканий из зоны исследований с целью получения данных аномалий полной напряженности магнитного поля (total magnetic intensity, TMI) или данных гравитационных аномалий Буге;
нанесение на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге и продолжение вверх нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге на множество предварительно определенных высот с целью получения множества нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге, Th на соответствующих уровнях, где h представляет высоту после продолжения вверх;
соответственно, вычисление наклонной производной TDRh данных аномалий TMI или данных гравитационных аномалий Буге на каждом уровне из нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге, Th, на каждом уровне;
применение многонаправленного обнаружения краев на основе горизонтального градиента соответственно к наклонной производной нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге на каждом уровне с целью получения краев объектов-источников магнитных или гравитационных аномалий на уровнях; и
утончение полученных краев объектов-источников магнитных или гравитационных аномалий на уровнях соответственно до однопиксельной ширины с использованием алгоритма определения морфологического скелета с целью получения множества карт структурного строения на соответствующих уровнях.
12. Способ по п. 11, отличающийся тем, что способ дополнительно включает наложение множества карт структурного строения на соответствующих уровнях с целью создания составной карты структурного строения.
13. Способ по п. 12, отличающийся тем, что способ дополнительно включает приведение краев, полученных из нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге, продолженных вверх на множество определенных высот, в соответствие структурам на разных глубинах, и наложение полученных карт структурного строения на глубинах с целью получения составной карты структурного строения, отражающей сведения о разных глубинах врезания.
14. Способ по п. 11, отличающийся тем, что этап применения многонаправленного обнаружения краев на основе горизонтального градиента соответственно к наклонной производной нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге на каждом уровне включает следующие этапы:
Соответственно, определение направленных производных наклонной производной TDRh в направлениях a и a+
Figure 00000044
как:
Figure 00000053
Figure 00000054
где D представляет производную первого порядка;
для высоты h и направления a — представление горизонтального градиента наклонной производной TDRh как:
Figure 00000055
где
Figure 00000047
— горизонтальный градиент;
определение модуля горизонтального градиента
Figure 00000056
как:
Figure 00000057
где аргумент горизонтального градиента имеет вид:
Figure 00000058
, и
затем определение точек, в которых модуль
Figure 00000059
представляет собой локальный максимум вдоль направления
Figure 00000060
градиента, как краевых точек объектов-источников магнитных или гравитационных аномалий на высоте h в направлении a;
в каждом направлении a формирование кривой путем соединения точек с локальным максимумом модуля горизонтального градиента наклонной производной TDRh в направлении, перпендикулярном градиенту, с целью создания края и
для той же высоты вычисление краев на множестве разных направлений a и принятие объединения краев в качестве краев объектов-источников магнитных или гравитационных аномалий на соответствующих уровнях,
где множество разных значений a направления соответственно полагают равными kπ/(2n-1), где k=0, 1, 2…(2n-1), и n — целое число, большее или равное 2.
15. Способ по п. 14, отличающийся тем, что способ дополнительно включает:
соответственно использование модуля горизонтального градиента в каждой краевой точке на карте структурного строения на каждом уровне для представления структурной глубины залегания в краевой точке на карте структурного строения на уровне с целью получения множества карт структурного строения, представляющих структурные глубины залегания на соответствующих уровнях.
16. Способ по п. 15, отличающийся тем, что способ дополнительно включает наложение множества карт структурного строения, представляющих структурные глубины залегания на соответствующих уровнях, с целью создания составной карты структурного строения с отражением глубин залегания.
17. Способ по п. 11, отличающийся тем, что способ дополнительно включает:
соответственно, вычисление трехмерного аналитического сигнала ASh на основе нанесенных на координатную сетку данных аномалий TMI или данных гравитационных аномалий Буге, Th, на каждом уровне для получения значения ASh каждой краевой точки с целью получения множества карт структурного строения, представляющих напряженность магнитного поля или интенсивность плотности краевых точек на соответствующих уровнях.
18. Способ по п. 17, отличающийся тем, что способ дополнительно включает наложение множества карт структурного строения, представляющих напряженность магнитного поля или интенсивность плотности в краевых точках на соответствующих уровнях, с целью создания составной карты структурного строения с отображением напряженности магнитного поля или составной карты структурного строения с отображением интенсивности плотности.
19. Способ по п. 11, отличающийся тем, что способ дополнительно включает выполнение обработки, связанной с подавлением шума, в отношении наклонных производных перед обнаружением краев.
20. Способ по п. 11, отличающийся тем, что способ применим для автоматического получения магнитного структурного строения в низкоширотной зоне; предпочтительно способ применим к данным магнитных измерений зоны с магнитным наклонением в диапазоне ±30°; и более предпочтительно способ применим к данным магнитных измерений зоны с магнитным наклонением в диапазоне ±20°.
RU2017133223A 2015-03-04 2016-03-04 Способ автоматического получения структурного строения из данных потенциального поля RU2664488C1 (ru)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
CN201510096981.9A CN104658037B (zh) 2015-03-04 2015-03-04 一种位场构造格架自动提取方法
CN201510096981.9 2015-03-04
CN201510303316.2A CN104965232B (zh) 2015-06-04 2015-06-04 低纬度地区磁构造格架自动提取方法
CN201510303316.2 2015-06-04
PCT/CN2016/075626 WO2016138874A1 (zh) 2015-03-04 2016-03-04 位场构造格架自动提取方法

Publications (1)

Publication Number Publication Date
RU2664488C1 true RU2664488C1 (ru) 2018-08-17

Family

ID=56848301

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2017133223A RU2664488C1 (ru) 2015-03-04 2016-03-04 Способ автоматического получения структурного строения из данных потенциального поля

Country Status (5)

Country Link
US (1) US10884161B2 (ru)
AU (1) AU2016228027B2 (ru)
CA (1) CA2978500C (ru)
RU (1) RU2664488C1 (ru)
WO (1) WO2016138874A1 (ru)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108008459B (zh) * 2017-11-28 2019-11-01 北京中科地物能源技术有限公司 一种获得剩余重力异常的方法及装置
CN110135506B (zh) * 2019-05-20 2021-02-09 杭州电子科技大学 一种应用于web的七类皮肤肿瘤检测方法
CN110414060B (zh) * 2019-06-28 2023-01-03 中国地质大学(武汉) 一种基于四阶谱矩的位场边界识别方法
CN110703347B (zh) * 2019-10-24 2021-04-16 中国石油化工股份有限公司 基于构造背景的重力断裂影像识别方法
CN111325766B (zh) * 2020-02-20 2023-08-25 腾讯科技(深圳)有限公司 三维边缘检测方法、装置、存储介质和计算机设备
CN112328955B (zh) * 2020-10-16 2024-02-27 中国地质调查局沈阳地质调查中心 重磁数据的处理方法、存储介质及装置
CN113111500A (zh) * 2021-03-31 2021-07-13 中国地质大学(北京) 一种基于重磁物理场使用二维经验模态分解异常分析方法
CN113421194B (zh) * 2021-06-04 2022-07-15 贵州省地质矿产勘查开发局 一种根据布格重力异常图像提取隐伏断层的方法
CN113779013B (zh) * 2021-09-18 2024-01-30 核工业航测遥感中心 基于地质约束的小范围磁场数据补缺处理方法
US11815647B1 (en) * 2022-04-20 2023-11-14 Chinese Academy Of Geological Sciences Gravity inversion method and system based on meshfree method
CN115236755B (zh) * 2022-07-25 2023-10-03 中国自然资源航空物探遥感中心 基于张量特征值的航磁异常边界检测方法、装置
CN117312898B (zh) * 2023-11-27 2024-03-15 山东省煤田地质规划勘察研究院 一种基于多重k均值聚类分析的找矿预测方法及系统
CN118051740B (zh) * 2024-02-27 2024-08-30 四川省地球物理调查研究所 融合多数据源的地质构造分析方法及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2169384C1 (ru) * 1999-12-17 2001-06-20 Закрытое акционерное общество "Петербургская геофизическая компания" Способ поиска нефтегазовых месторождений
US20030169253A1 (en) * 2002-03-07 2003-09-11 Samsung Electronic Co., Ltd. Method of wavelets-based multiresolution representation of three-dimensional image object
CN101256676A (zh) * 2008-01-31 2008-09-03 中国地质科学院矿产资源研究所 位场多方向多尺度边缘检测方法
RU2401443C2 (ru) * 2008-03-17 2010-10-10 Общество с ограниченной ответственностью "Научно-производственное предприятие "Нейво" Способ обнаружения и отображения фигуры газонефтяной лог-трубки
CN102236108A (zh) * 2010-05-06 2011-11-09 中国石油天然气集团公司 一种磁性地表三维地形改正方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7486811B2 (en) * 1996-09-16 2009-02-03 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual examination of objects, such as internal organs
US7127104B2 (en) * 2004-07-07 2006-10-24 The Regents Of The University Of California Vectorized image segmentation via trixel agglomeration
US8041141B2 (en) * 2006-06-30 2011-10-18 The University Of Louisville Research Foundation, Inc. Method and software for shape representation with curve skeletons
CN101520518B (zh) 2008-02-25 2011-01-12 中国石油集团东方地球物理勘探有限责任公司 一种利用重磁电异常的组合特征识别火成岩岩性的方法
CA2721008A1 (en) * 2008-04-11 2009-10-15 Terraspark Geosciences, Llc Visulation of geologic features using data representations thereof
MY164574A (en) * 2008-05-22 2018-01-15 Exxonmobil Upstream Res Co Seismic horizon skeletonization
US8126247B2 (en) * 2009-05-19 2012-02-28 National Tsing Hua University Image preprocessing system for 3D image database construction
US9558588B2 (en) * 2012-06-26 2017-01-31 Schlumberger Technology Corporation Method for building a 3D model of a rock sample
CN102937725B (zh) * 2012-11-12 2016-04-20 中国科学院地质与地球物理研究所 一种基于过渡区与相叠合的位场异常边缘增强方法
CN103955007A (zh) 2014-05-20 2014-07-30 中国石油化工股份有限公司胜利油田分公司西部新区研究院 复杂山前构造带的综合建模方法及建立的地质结构模型
CN104965232B (zh) * 2015-06-04 2017-11-10 中国地质科学院矿产资源研究所 低纬度地区磁构造格架自动提取方法
CN104658037B (zh) * 2015-03-04 2017-04-05 中国地质科学院矿产资源研究所 一种位场构造格架自动提取方法
US10136869B2 (en) * 2016-03-25 2018-11-27 Perkinelmer Health Sciences, Inc. Systems and methods for characterizing a central axis of a bone from a 3D anatomical image

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2169384C1 (ru) * 1999-12-17 2001-06-20 Закрытое акционерное общество "Петербургская геофизическая компания" Способ поиска нефтегазовых месторождений
US20030169253A1 (en) * 2002-03-07 2003-09-11 Samsung Electronic Co., Ltd. Method of wavelets-based multiresolution representation of three-dimensional image object
CN101256676A (zh) * 2008-01-31 2008-09-03 中国地质科学院矿产资源研究所 位场多方向多尺度边缘检测方法
RU2401443C2 (ru) * 2008-03-17 2010-10-10 Общество с ограниченной ответственностью "Научно-производственное предприятие "Нейво" Способ обнаружения и отображения фигуры газонефтяной лог-трубки
CN102236108A (zh) * 2010-05-06 2011-11-09 中国石油天然气集团公司 一种磁性地表三维地形改正方法

Also Published As

Publication number Publication date
CA2978500C (en) 2022-07-05
US10884161B2 (en) 2021-01-05
CA2978500A1 (en) 2016-09-09
US20180052251A1 (en) 2018-02-22
AU2016228027B2 (en) 2018-11-22
WO2016138874A1 (zh) 2016-09-09
AU2016228027A1 (en) 2017-09-21

Similar Documents

Publication Publication Date Title
RU2664488C1 (ru) Способ автоматического получения структурного строения из данных потенциального поля
Gunn et al. Interpreting aeromagnetic data in areas of limited outcrop
US5724309A (en) Method for geophysical processing and interpretation using instantaneous phase and its derivatives and their derivatives
CN104965232B (zh) 低纬度地区磁构造格架自动提取方法
Boniger et al. Subsurface utility extraction and characterization: Combining GPR symmetry and polarization attributes
Qi et al. A workflow to skeletonize faults and stratigraphic features
Zahra et al. Application of high-pass filtering techniques on gravity and magnetic data of the eastern Qattara Depression area, Western Desert, Egypt
Ansari et al. A new edge detection method based on the analytic signal of tilt angle (ASTA) for magnetic and gravity anomalies
Martinez et al. Lithologic characterization using airborne gravity gradient and aeromagnetic data for mineral exploration: A case study in the Quadrilátero Ferrífero, Brazil
Zhang et al. NAV-Edge: Edge detection of potential-field sources using normalized anisotropy variance
Sun et al. Edge enhancement of potential field data using spectral moments
Pham et al. Analysis of gravity data for extracting structural features of the northern region of the Central Indian Ridge
CN104658037B (zh) 一种位场构造格架自动提取方法
Koyan et al. 3D ground-penetrating radar attributes to generate classified facies models: A case study from a dune island
Pham A stable method for detecting the edges of potential field sources
Pei et al. Multi-scale edge detection method for potential field data based on two-dimensional variation mode decomposition and mathematical morphology
Karbalaali et al. Identification of shallow geohazard channels based on discontinuity seismic attributes in the South Caspian Sea
CN114609683B (zh) 一种砂岩型铀矿砂体有利区确定方法及系统
Amonpantang et al. Multiattribute analysis for channel element discrimination in the Taranaki Basin, offshore New Zealand
Zhang et al. Automatic lineament extraction from potential-field images using the Radon transform and gradient calculation
Eshaghzadeh et al. Magnetic field interpretation using singular value decomposition method based on correlation coefficient of eigenimages
Zahra et al. Tectonic and structural setting of the northeastern central Gulf of Suez area using aeromagnetic data
Shi et al. New enhancement filters for geological mapping
Koyan et al. 3D ground-penetrating radar data analysis and interpretation using attributes based on the gradient structure tensor
Guo et al. Reduction-to-the-pole method for aeromagnetic three-component data in different latitudes