CN110853114B - 一种利用重力异常线对海底热液硫化物矿点识别方法 - Google Patents
一种利用重力异常线对海底热液硫化物矿点识别方法 Download PDFInfo
- Publication number
- CN110853114B CN110853114B CN201910973195.0A CN201910973195A CN110853114B CN 110853114 B CN110853114 B CN 110853114B CN 201910973195 A CN201910973195 A CN 201910973195A CN 110853114 B CN110853114 B CN 110853114B
- Authority
- CN
- China
- Prior art keywords
- grid
- gravity anomaly
- layer
- anomaly
- gravity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/56—Information retrieval; Database structures therefor; File system structures therefor of still image data having vectorial format
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种利用重力异常线对海底热液硫化物矿点识别方法,导入重力异常原始网格数据,根据研究需要转换成相应需要精度的栅格图层;计算重力异常梯度方向,计算网格梯度方向累计的次数,计算域内的所有网格都将得到一个累计次数,并形成新的栅格图层;用窗口从第一个网格单元开始进行梯度计算,边界网格以窗口内现有网格数目计算结果为准,记录每个网格的最大梯度方向,并记录该网格最大梯度方向的网格编号;计算每个网格单元被其它网格单元梯度指向的次数,提取0值点,其连接成的线即为重力异常线。本发明的有益效果是能够识别洋中脊热液硫化物矿点空间分布。
Description
技术领域
本发明属于海底热液成矿地质技术领域,涉及一种利用自由空气重力异常线对海底热液硫化物矿点识别方法。
背景技术
在大西洋,前期的诸多研究(Rona et al.,1980,2001;Bostrom et al.,1972;Batuyev et al.,1994;German et al.,1998;Bates,2005等)对热液矿床的地质构造背景、矿体规模、矿化作用过程等方面,进行了研究和推理。目前热液调查中所用到的AUV、ROV、TVG、声学、光学、地球物理等热液信息捕捉和获取手段,只能通过现场原位测量,若仅采取这样的方式寻找热液点,就失去了经济探矿的意义。因此,从数学方法尤其是特征提取方法上改变思路和策略,充分挖掘隐含的微异常信息,使其能从不同角度反映成矿特征,从而开展硫化物资源空间分布预测和评价是硫化物资源调查和研究中的重要任务。
发明内容
本发明的目的在于提供一种利用重力异常线对海底热液硫化物矿点识别方法,本发明的有益效果是能够识别洋中脊热液硫化物矿点空间赋存位置区域,据此规划大洋中冀热液硫化物成矿远景区域。
本发明所采用的技术方案是按照以下步骤进行:
步骤1:利用ArcGIS软件,导入重力异常原始网格数据,根据研究需要转换成相应需要精度的栅格图层;
步骤2:将栅格图层,用ArcGIS的ArcToolbox工具箱的Flow Direction方法计算重力异常梯度方向,之后用Flow Accumulation计算网格梯度方向累计的次数,计算域内的所有网格都将得到一个累计次数,并形成新的栅格图层;
步骤3:用Raster Calculator栅格计算器求取上述梯度方向累计次数值为0的网格,形成新的栅格图层,即得到表征重力异常线性异常的栅格图层;
步骤4:利用ArcGIS软件的ArcToolbox的重分类工具对步骤3栅格图层进行数据筛选分类,调整结果;
步骤5:用ArcGIS软件建立线性图层,添加字段,并打开到编辑状态;用ArcGIS软件的ArcScan模块,在Vectorization中选择要写入的图层和线宽度值存放的字段等,将栅格数据转化为矢量图层;
步骤6:将矢量线性图层以1倍于重力异常数据精度的半径建立缓冲区;
步骤7:将重力异常低线性异常和重力异常高线性异常建立缓冲区之后的面状图层转化为栅格图层;
步骤8:将上述两个栅格图层进行图层计算,提取交叉区域和临近的线性异常高缓冲区域,此即重力异常线异常的最终提取结果。
步骤9:步骤8结果与已发现热液矿点之间的具有密切的空间关系,据此,可划分洋中脊区热液成矿远景区。
附图说明
图1是梯度方向;
图2是梯度方向累计数值;
图3是重力异常线性异常;
图4是重力异常线性异常矢量图;
图5是重力异常线处理结果。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
步骤1:利用ArcGIS软件,导入重力异常原始数据(一般为网格数据),根据研究需要转换成相应需要精度的栅格图层;
步骤2:将上述步骤1的栅格图层,用ArcGIS的ArcToolbox工具箱的FlowDirection方法计算重力异常梯度方向,之后用Flow Accumulation计算网格梯度方向累计的次数,计算域内的所有网格都将得到一个累计次数,并形成新的栅格图层。
Flow Direction方法计算中的重力异常梯度方向如图1,表示某重力网格单元(g5)及周围网格单元分布,g1、g2、g3、g4、g5、g6、g7、g8、g9为重力异常数据在网格单元中心点取值。
则网格单元g5的差分计算为:
Gg5=(g5-gi)/di i=1,2,3,4,6,……
其中,Gg5为示例窗口中心点网格梯度,d为网格中心点间距离。
Flow Accumulation计算中的梯度方向用图2所示方式,将上述g5网格单元的方向标记为相应的值。在网格单元经过梯度计算后有如左图所示的方向标记的情况下,统计每个单元被梯度指向的次数,每个单元的梯度指向累计次数值为右图所示,0值点构成的区域则为相对高异常区域。
步骤3:用Raster Calculator栅格计算器求取上述梯度方向累计次数值为0的网格,形成新的栅格图层,即得到表征重力异常线性异常的栅格图层。
(1)自由空气重力异常“脊”提取:
输入网格化处理好的重力异常图层,用Flow Direction计算重力异常梯度方向,之后用Flow Accumulation计算网格梯度方向累计的次数,计算域内的所有网格都将得到一个累计次数,并形成新的栅格图层,用Raster Calculator栅格计算器求取梯度方向累计次数为0的网格,形成新的栅格图层,即得到重力异常“脊”栅格图层。
(2)自由空气重力异常“谷”提取:
输入网格化处理好的重力异常图层,用Raster Calculator将上述输入图层每个栅格的重力异常值取反,之后重复上述(1)中的操作,即得到重力异常“谷”栅格图层。
步骤4:利用ArcGIS软件的ArcToolbox的重分类工具对步骤3栅格图层进行数据筛选分类,调整结果,得到重力异常线性异常栅格图层;
步骤5:用ArcGIS软件的ArcScan进行重力异常线性栅格图层矢量化,并将矢量化的重力异常线性高和线性低图层分别以1倍于重力异常数据精度的半径建立缓冲区,之后再分别转化为栅格图层,将两个栅格图层用Raster Calculator进行图层计算,提取交叉区域及临近的线性异常高缓冲区域,此即重力异常线的最终提取结果。
步骤6:前述提取结果中的重力异常高线性异常和重力异常低线性异常,这里形象地称为重力异常脊和谷,重力异常高的区域这里称为重力异常顶,提取结果与已发现硫化物矿点空间对比分析可以看到,热液点基本分布在重力异常脊与谷的交叉区域及其临近的重力异常脊附近。如图3所示为大西洋脊区重力异常线提取结果,其中的黑点为已发现热液点。
第一种情况,即热液点附近有不规则的重力异常顶特征,热液点位于重力异常脊与重力异常谷的交界区域。典型代表如TAG热液区。TAG热液区东侧有平行于脊轴的重力异常顶,西侧位于洋中脊断裂带上,且一条重力异常脊由南东方向向北西向穿过TAG热液区,表明此区有物质盈余异常,代表断裂构造的重力异常谷与代表物质盈余异常的重力异常脊交汇部位即是容易发现热液点的区域。大西洋脊80%的已发现热液点均明显分布在具有上述两种特征的位置。
第二种情况的典型代表是北大西洋的LostCity热液区,其发育区域重力异常脊在重力异常顶处的最高值部位,有表现为放射状脊线,或可说明地壳内岩浆物质的积聚和外散而导致的重力异常变化,这两个特征符合物质的积聚、盈余和散溢的特征。
以重力异常上述特征为推断依据,大西洋中脊尚有数个可能存在但尚未发现需要重点探测的热液区域。第一类热液区,其成矿体系比较封闭,面积和体积都较大,更加值得关注,类似于这样的目标区域,据统计南大西洋中脊附近约有12个,北大西洋中脊附近约30个。第二类脊、谷交叉特征的热液硫化物远景区在南大西洋至少有上百个。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (1)
1.一种利用重力异常线对海底热液硫化物矿点识别方法,其特征在于按照以下步骤进行:
步骤1:利用ArcGIS软件,导入重力异常原始网格数据,根据研究需要转换成相应需要精度的栅格图层;
步骤2:将栅格图层,用ArcGIS的ArcToolbox工具箱的Flow Direction方法计算重力异常梯度方向,之后用Flow Accumulation计算网格梯度方向累计的次数,计算域内的所有网格都将得到一个累计次数,并形成新的栅格图层;
步骤3:用Raster Calculator栅格计算器求取上述梯度方向累计次数值为0的网格,形成新的栅格图层,即得到表征重力异常线性异常的栅格图层;
步骤4:利用ArcGIS软件的ArcToolbox的重分类工具对步骤3栅格图层进行数据筛选分类,调整结果;
步骤5:用ArcGIS软件建立线性图层,添加字段,并打开到编辑状态;用ArcGIS软件的ArcScan模块,在Vectorization中选择要写入的图层和线宽度值存放的字段,将栅格数据转化为矢量图层;
步骤6:将矢量线性图层以1倍于重力异常数据精度的半径建立缓冲区;
步骤7:将重力异常低线性异常和重力异常高线性异常建立缓冲区之后的面状图层转化为栅格图层;
步骤8:将上述两个栅格图层进行图层计算,提取交叉区域和临近的线性异常高缓冲区域,此即重力异常线异常的最终提取结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910973195.0A CN110853114B (zh) | 2019-10-14 | 2019-10-14 | 一种利用重力异常线对海底热液硫化物矿点识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910973195.0A CN110853114B (zh) | 2019-10-14 | 2019-10-14 | 一种利用重力异常线对海底热液硫化物矿点识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110853114A CN110853114A (zh) | 2020-02-28 |
CN110853114B true CN110853114B (zh) | 2023-08-22 |
Family
ID=69596262
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910973195.0A Active CN110853114B (zh) | 2019-10-14 | 2019-10-14 | 一种利用重力异常线对海底热液硫化物矿点识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110853114B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011242231A (ja) * | 2010-05-18 | 2011-12-01 | Tokyo Univ Of Marine Science & Technology | 地下密度構造の推定方法 |
CN103605168A (zh) * | 2013-10-12 | 2014-02-26 | 国家海洋局第二海洋研究所 | 一种海底多金属硫化物综合信息快速找矿方法 |
CN108090929A (zh) * | 2017-12-04 | 2018-05-29 | 国家海洋局第海洋研究所 | 矿区线性异常分析提取新型方法 |
CN108133480A (zh) * | 2017-12-04 | 2018-06-08 | 国家海洋局第海洋研究所 | 一种热液喷口异常地形提取方法 |
-
2019
- 2019-10-14 CN CN201910973195.0A patent/CN110853114B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011242231A (ja) * | 2010-05-18 | 2011-12-01 | Tokyo Univ Of Marine Science & Technology | 地下密度構造の推定方法 |
CN103605168A (zh) * | 2013-10-12 | 2014-02-26 | 国家海洋局第二海洋研究所 | 一种海底多金属硫化物综合信息快速找矿方法 |
CN108090929A (zh) * | 2017-12-04 | 2018-05-29 | 国家海洋局第海洋研究所 | 矿区线性异常分析提取新型方法 |
CN108133480A (zh) * | 2017-12-04 | 2018-06-08 | 国家海洋局第海洋研究所 | 一种热液喷口异常地形提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110853114A (zh) | 2020-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kuhn et al. | Lithologic mapping using Random Forests applied to geophysical and remote-sensing data: A demonstration study from the Eastern Goldfields of Australia | |
Agterberg | Geomathematics: theoretical foundations, applications and future developments | |
Cheng | Spatial and scaling modelling for geochemical anomaly separation | |
Englund et al. | GEO-EAS (Geostatistical environmental assessment software) user's guide | |
Koike et al. | Neural network-based estimation of principal metal contents in the Hokuroku district, northern Japan, for exploring Kuroko-type deposits | |
CN101504659B (zh) | 一种基于广义自相似性原理提取地学空间信息的方法 | |
Liu et al. | Compositional balance analysis: an elegant method of geochemical pattern recognition and anomaly mapping for mineral exploration | |
Yang et al. | Quantification of uncertainty associated with evidence layers in mineral prospectivity mapping using direct sampling and convolutional neural network | |
CN116665067B (zh) | 基于图神经网络的找矿靶区优选系统及其方法 | |
Bruce et al. | Unconformity-type uranium systems: A comparative review and predictive modelling of critical genetic factors | |
Leung et al. | Sample truncation strategies for outlier removal in geochemical data: the MCD robust distance approach versus t-SNE ensemble clustering | |
Behera et al. | Gold prospectivity mapping in the Sonakhan Greenstone Belt, Central India: A knowledge-driven guide for target delineation in a region of low exploration maturity | |
Alqahtani et al. | Machine learning techniques in predicting bottom hole temperature and remote sensing for assessment of geothermal potential in the Kingdom of Saudi Arabia | |
Ferrer et al. | Modeling the uncertainty in the layout of geological units by implicit boundary simulation accounting for a preexisting interpretive geological model | |
CN110853114B (zh) | 一种利用重力异常线对海底热液硫化物矿点识别方法 | |
Feng et al. | Comparative Study on Three Autoencoder‐Based Deep Learning Algorithms for Geochemical Anomaly Identification | |
CN110717618A (zh) | 基于多指标综合要素的海底热液硫化物资源评价预测方法 | |
Williams et al. | A predictive geospatial exploration model for Mississippi valley type Pb–Zn mineralization in the southeast Missouri lead district | |
Ribeiro et al. | Geochemical multifractal modeling of soil and stream sediment data applied to gold prospectivity mapping of the Pitangui Greenstone Belt, northwest of Quadrilátero Ferrífero, Brazil | |
NL2028764B1 (en) | Method for identifying submarine hydrothermal sulphide ore spots by utilizing gravity anomaly lines | |
Louloudis et al. | Multi-seam coal deposit modeling via principal component analysis & GIS | |
Partington | Exploration targeting using GIS: More than a digital light table | |
CN113111928A (zh) | 一种基于地学数据库的半监督学习矿产资源定量预测方法 | |
Stensgaard et al. | Assessment of mineral potential using cross-validation techniques and statistical analysis: a case study from the Paleoproterozoic of West Greenland | |
NL2028908B1 (en) | Method for predicting submarine hydrothermal metallogenic prospective area based on the features of multi-index elements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |