RU2664488C1 - Способ автоматического получения структурного строения из данных потенциального поля - Google Patents
Способ автоматического получения структурного строения из данных потенциального поля Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 104
- 238000003708 edge detection Methods 0.000 claims abstract description 23
- 238000013507 mapping Methods 0.000 claims abstract description 11
- 230000000877 morphologic effect Effects 0.000 claims abstract description 9
- 239000002131 composite material Substances 0.000 claims description 49
- 238000005259 measurement Methods 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 21
- 239000003086 colorant Substances 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000011160 research Methods 0.000 claims description 14
- 230000005484 gravity Effects 0.000 claims description 10
- 238000009499 grossing Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001629 suppression Effects 0.000 claims 1
- 239000002184 metal Substances 0.000 abstract description 11
- 229910052751 metal Inorganic materials 0.000 abstract description 11
- 230000000694 effects Effects 0.000 abstract description 4
- 230000003750 conditioning effect Effects 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 230000015572 biosynthetic process Effects 0.000 description 14
- 238000001514 detection method Methods 0.000 description 11
- 230000009467 reduction Effects 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 229910052500 inorganic mineral Inorganic materials 0.000 description 6
- 239000011707 mineral Substances 0.000 description 6
- 230000005358 geomagnetic field Effects 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 5
- 238000000926 separation method Methods 0.000 description 5
- 230000035515 penetration Effects 0.000 description 4
- 206010039509 Scab Diseases 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000005415 magnetization Effects 0.000 description 3
- 230000002547 anomalous effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- CWYNVVGOOAEACU-UHFFFAOYSA-N Fe2+ Chemical compound [Fe+2] CWYNVVGOOAEACU-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
- G01V11/002—Details, e.g. power supply systems for logging instruments, transmitting or recording data, specially adapted for well logging, also if the prospecting method is irrelevant
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric 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/081—Electric 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20036—Morphological image processing
- G06T2207/20044—Skeletonization; Medial axis transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth 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 как:
определение вейвлет-функции в направлении a как:
где D представляет производную первого порядка;
для уровня s и местоположения (x, y), — определение вейвлет-преобразования гравитационной аномалии или магнитной аномалии, f0(x, y), в направлении a как:
где * представляет операцию свертки,
fz(x, y) — гравитационная аномалия или магнитная аномалия посредством продолжения вверх f0(x, y) от нулевой высоты до высоты z=sz0 и полученная путем продолжения вверх гравитационной аномалии или магнитной аномалии, fz0(x, y), измеренной на высоте z0 измерения, на высоту z-z0,
кроме того, для уровня s и местоположения (x, y), — определение вейвлет-преобразования гравитационной аномалии или магнитной аномалии, f0(x, y), в направлении a+ как:
поэтому двумерное вейвлет-преобразование f0(x, y) записывается как в градиенте:
при этом для местоположения (x, y), уровня s и направления a двумерное направленное вейвлет-преобразование f0(x, y) прямо пропорционально градиенту fz(x, y), и двумерное направленное вейвлет-преобразование f0(x, y) можно представить градиентом fz(x, y),
где соответствующий аргумент градиента вдоль горизонтальной оси имеет вид:
определение точки, в которой модуль M представляет собой локальный максимум вдоль направления градиента , как краевой точки,
формирование кривой путем соединения точек с локальным максимумом модуля градиента в направлении, перпендикулярном градиенту, с целью создания края, и
для той же высоты вычисление краев на множестве разных направлений 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+ как:
где D представляет производную первого порядка;
для высоты h и направления a, — представление горизонтального градиента наклонной производной TDRh как:
где аргумент горизонтального градиента имеет вид:
затем определение точек, в которых модуль представляет собой локальный максимум вдоль направления градиента, как краевых точек объектов-источников магнитных или гравитационных аномалий на высоте 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 определяется как:
Вейвлет-функция в направлении a определяется как:
где D представляет производную первого порядка.
Для уровня s и местоположения (x, y) вейвлет-преобразование гравитационной аномалии или магнитной аномалии f0(x, y) в направлении a определяется как:
где * представляет операцию свертки,
в соответствии с формулой продолжения потенциального поля вверх, известно, что:
fz(x, y) — гравитационная аномалия или магнитная аномалия посредством продолжения вверх f0(x, y) от нулевой высоты до высоты z=sz0, и полученная путем продолжения вверх гравитационной аномалии или магнитной аномалии, fz0(x, y), измеренной на высоте измерения z0, для высоты z-z0,
Аналогично, для уровня s и местоположения (x, y) вейвлет-преобразование гравитационной аномалии или магнитной аномалии f0(x, y) в направлении a+ определяется как:
Градиент можно использовать для представления двумерного вейвлет-преобразования f0(x, y) как:
Приведенное выше уравнение устанавливает отношение между горизонтальным градиентом гравитационной аномалии или магнитной аномалии fz(x, y) на любой высоте z (z>z0) и двумерным направленным вейвлет-преобразованием гравитационной аномалии или магнитной аномалии f0(x, y) на нулевой высоте.
Для местоположения (x, y), уровня s и направления a двумерное вейвлет-преобразование f0(x, y) прямо пропорционально градиенту fz(x, y), и двумерное вейвлет-преобразование f0(x, y) можно представить градиентом fz(x, y).
где соответствующий аргумент градиента по горизонтальной оси имеет вид:
точка, в которой модуль M представляет собой локальный максимум вдоль направления аргумента, определяется как краевая точка.
Для каждой высоты 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 модуль градиента каждой краевой точки не изменяется и представляет собой постоянное значение. Модуль градиента краевой точки принимается в качестве значения структурной интенсивности в данной краевой точке с целью установления карт структурного строения с отображением интенсивности на разных уровнях, или разных глубинах.
На этапе 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 на соответствующем уровне в соответствии со следующей формулой:
где VDRh и THDRh, соответственно, представляют вертикальную производную первого порядка и полную горизонтальную производную нанесенных на координатную сетку данных Th аномалий TMI:
На этапе 604 к наклонной производной TDRh на каждом уровне, соответственно, применяют многонаправленное обнаружение краев на основе горизонтального градиента.
Направленные производные наклонной производной TDRh в направлениях a и a+, соответственно, определяют как:
где D представляет производную первого порядка.
Горизонтальный градиент наклонной производной TDRh для высоты h и направления a представляют как:
где соответствующий аргумент градиента имеет вид:
полная горизонтальная производная наклонных производных (сокращенно: TDR_THDR) данных аномалий TMI имеет вид:
Таким образом, точки, в которых модуль представляет собой локальный максимум вдоль направления аргумента , представляют собой краевые точки объектов-источников магнитных аномалий на высоте 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 модуль горизонтального градиента в каждой краевой точке не изменяется и представляет собой постоянное значение.
Из-за свойства значения арктангенса тригонометрической функции, вне зависимости от того, насколько велики амплитуды VDRh и THDRh, амплитуда TDRh ограничена промежутком между -π/2 и +π/2. Поэтому модуль горизонтального градиента в каждой краевой точке имеет малую корреляцию с амплитудой аномалий TMI, значение отражает глубины залегания объектов-источников, величина значения обратно пропорциональна глубине залегания, и чем больше значение, тем меньше глубина залегания (см. публикацию «New Insights into Magnetic Derivatives for Structural Mapping (Новые представления о магнитных производных для структурного картирования)», Bruno Verduzco и др., The Leading Edge, т. 23, №. 2, 2004).
Структурная глубина залегания, отражаемая модулем горизонтального градиента наклонной производной TDRh, представляет толщину покровного слоя над структурным поясом, и разные глубины, соответствующие продолжению вверх на множество предварительно определенных высот, отражают глубину врезания структурного пояса, которая представляет глубину протяженности структурного пояса вниз.
Модуль градиента каждой краевой точки используют для представления связанной с краем глубины залегания краевой точки на уровне с целью установления карт структурного строения, отражающих связанную со структурой глубину залегания на разных уровнях, или глубинах. Изменения в значениях модуля градиента одного и того же структурного пояса отражают глубины залегания разных участков данного структурного пояса.
На этапе 608 карты структурного строения, представляющие глубины залегания на множестве уровней, накладывают одну на другую с целью получения составной карты структурного строения с отражением глубины залегания.
Значения модуля горизонтального градиента производной TDRh краевых точек на уровнях, или глубинах, накладывают одно на другое и представляют разными постепенно изменяющимися цветами с целью выделения сведений об изменениях глубины залегания структурных строений разных диапазонов глубин врезания.
На этапе 609 для нанесенных на координатную сетку данных Th аномалий TMI на каждой высоте вычисляют трехмерный аналитический сигнал ASh в соответствии с формулой:
Амплитуда получаемого в результате аналитического сигнала 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 как:
определение вейвлет-функции в направлении a как:
где D представляет производную первого порядка;
для уровня s и местоположения (x, y) — определение вейвлет-преобразования гравитационной аномалии или магнитной аномалии f0(x, y) в направлении a как:
где * представляет операцию свертки,
fz(x, y) — гравитационная аномалия или магнитная аномалия посредством продолжения вверх f0(x, y) от нулевой высоты до высоты z=sz0, полученная путем продолжения вверх гравитационной аномалии или магнитной аномалии, fz0(x, y), измеренной на высоте z0 измерения, на высоту z-z0,
тогда двумерное направленное вейвлет-преобразование f0(x, y) записывают как градиент:
где соответствующий аргумент градиента вдоль горизонтальной оси имеет вид:
в каждом направлении 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 или данных гравитационных аномалий Буге на каждом уровне включает следующие этапы:
где D представляет производную первого порядка;
для высоты h и направления a — представление горизонтального градиента наклонной производной TDRh как:
где аргумент горизонтального градиента имеет вид:
в каждом направлении 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°.
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)
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)
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)
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 |
-
2016
- 2016-03-04 RU RU2017133223A patent/RU2664488C1/ru active
- 2016-03-04 WO PCT/CN2016/075626 patent/WO2016138874A1/zh active Application Filing
- 2016-03-04 CA CA2978500A patent/CA2978500C/en active Active
- 2016-03-04 AU AU2016228027A patent/AU2016228027B2/en active Active
- 2016-03-04 US US15/555,153 patent/US10884161B2/en active Active
Patent Citations (5)
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 |