CN103529486B - 一种地球化学异常圈定方法 - Google Patents
一种地球化学异常圈定方法 Download PDFInfo
- Publication number
- CN103529486B CN103529486B CN201310274760.7A CN201310274760A CN103529486B CN 103529486 B CN103529486 B CN 103529486B CN 201310274760 A CN201310274760 A CN 201310274760A CN 103529486 B CN103529486 B CN 103529486B
- Authority
- CN
- China
- Prior art keywords
- subarea
- value
- area
- district
- survey
- 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
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于地球化学异常圈定技术领域,具体涉及一种能够快速、有效圈定地球化学异常的方法。包括以下步骤:根据测区实际地质情况,确定子区的形状;针对子区形状,确定子区的大小;以子区在横向和纵向上轴线长度的1/3至1/2为横向间距和纵向间距进行滑动,划分测区;针对每个子区,采用含量与面积分形模型来确定该子区的异常下限值;将异常下限值作为子区中心位置处的背景值,采用克里格方法对整个测区范围内的子区背景值进行空间插值,得到整个测区的地球化学背景值;将整个测区的地球化学数据减去地球化学背景值,得到整个测区的地球化学异常值。本发明技术方案能够反映实际地质情况下地球化学异常的分布特征,对地质找矿工作具有重要的实际意义。
Description
技术领域
本发明属于地球化学异常圈定技术领域,具体涉及一种能够快速、有效圈定地球化学异常的方法。
背景技术
地球化学异常是矿产资源勘查中的重要决策依据之一,对其进行准确地圈定对地质找矿工作具有重要的实际意义。传统的地球化学异常圈定的方法,以整个测区地球化学数据的平均值加1至3倍标准差作为异常下限,圈定地球化学异常范围。在实际工作中,测区的地质情况往往比较复杂,且岩性的分布也不可能十分均一,不同的岩性和地层具有不同的地球化学异常下限值。如果整个测区采用同一数值作为下限,不但难以反映局部地段之间的地球化学场差异,而且容易遗失对地质找矿工作具有重要意义的异常信息。
在长期的地质过程中,多种地质作用相互叠加,使地球化学场表现为复杂的叠加场。为了有效地圈定地球化学异常,需要将叠加的地球化学场分解为正常场和异常场。分形理论指出,如果某一随机过程可以用各种等级的空间尺度等几率去描述,那么由该过程产生的现象往往具有分形特征。地球化学场便是由这种随机过程形成的复杂景观,因此用分形几何方法去刻划地球化学场的特征是可行的。
因此,亟需开发一种符合测区实际地质情况且能够有效圈定地球化学异常范围的方法,从而对地质找矿工作的进一步发展提供重要支撑。
发明内容
本发明的目的在于提供一种能够快速、有效圈定地球化学异常的方法。
为解决上述技术问题,本发明采取的技术方案是:
一种地球化学异常圈定方法,包括以下步骤:
(1)根据测区实际地质情况,确定子区的形状:
在沉积盆地地区,子区的形状为矩形;
在南方山区地区,无侵入岩体存在时,子区的形状为圆形;
在南方山区地区,有侵入岩体存在时,子区的形状为椭圆形;
(2)针对步骤(1)中选择的子区形状,确定子区的大小:
在沉积盆地地区,矩形子区的长度≥10km,宽度≥10km;
在南方山区地区,圆的直径长度、椭圆的长轴和短轴长度都∈1~10km;
(3)根据步骤(1)和(2)中确定的子区形状和大小,以子区在横向和纵向上轴线长度的1/3至1/2为横向间距和纵向间距进行滑动,划分测区;
(4)针对每个子区,采用含量与面积分形模型来确定该子区的异常下限值;
含量与面积分形模型由式(1)确定:
A(ρ>υ)∝ρ-α(1)
式中,A(ρ>υ)表示元素含量大于某一值υ的等值线所围成的面积,在元素含量的最大值与最小值之间均匀设定n个不同的阈值υn,统计含量值大于υn的所有元素所围成的面积An;
在logA和logυ的双对数坐标图上,利用最小二乘法拟合出两段直线,将两段直线的交叉点所对应的元素含量值作为异常下限值;
(5)将步骤(4)中得到的异常下限值作为子区中心位置处的背景值,采用克里格方法对整个测区范围内的子区背景值进行空间插值,得到整个测区的地球化学背景值;
(6)将整个测区的地球化学数据减去步骤(5)中得到的地球化学背景值,得到整个测区的地球化学异常值。
进一步的,如上所述的一种地球化学异常圈定方法,步骤(3)中,划分测区时的具体原则为:
子区形状为矩形时,分别以矩形的长度和宽度的1/3至1/2为横向间距和纵向间距进行滑动,将整个测区划分为若干个子区;
子区的形状为圆形时,分别以圆形直径长度的1/3至1/2为横向间距和纵向间距进行滑动,将整个测区划分为若干个子区;
子区的形状为椭圆形时,分别以椭圆长轴和短轴长度的1/3至1/2为横向间距和纵向间距进行滑动,将整个测区划分为若干个子区。
进一步的,如上所述的一种地球化学异常圈定方法,步骤(4)中,在元素含量的最大值与最小值之间均匀设定15~20个不同的阈值υn。
利用本发明技术方案后,对地球化学异常进行圈定的异常能够反映实际地质情况下地球化学异常的分布特征,对地质找矿工作具有重要的实际意义。
附图说明
图1是本发明所提供的一种地球化学异常圈定方法流程图。
具体实施方式
下面结合附图对本发明技术方案进行详细说明。
如图1所示,本发明一种地球化学异常圈定方法,包括以下步骤:
(1)根据测区实际地质情况,确定子区的形状:
在沉积盆地等地质条件相对均一的地区,子区的形状为矩形;
在南方山区等地形和地质条件较为复杂的地区,无侵入岩体存在时,子区的形状为圆形;有侵入岩体存在时,子区的形状为椭圆形;
(2)针对步骤(1)中选择的子区形状,确定子区的大小:
在沉积盆地地区等地质条件相对均一的地区,矩形子区的长度≥10km,宽度≥10km;
在南方山区地区等地质条件相对复杂的地区,圆的直径长度、椭圆的长轴和短轴长度都∈1~10km;
(3)根据步骤(1)和(2)中确定的子区形状和大小,以子区在横向和纵向上轴线长度的1/3至1/2为横向间距和纵向间距进行滑动,划分测区;划分测区时的具体原则为:
子区形状为矩形时,分别以矩形的长度和宽度的1/3至1/2为横向间距和纵向间距进行滑动,将整个测区划分为若干个子区;
子区的形状为圆形时,分别以圆形直径长度的1/3至1/2为横向间距和纵向间距进行滑动,将整个测区划分为若干个子区;
子区的形状为椭圆形时,分别以椭圆长轴和短轴长度的1/3至1/2为横向间距和纵向间距进行滑动,将整个测区划分为若干个子区。
(4)针对每个子区,采用含量与面积分形模型来确定该子区的异常下限值;
含量与面积分形模型由式(1)确定:
A(ρ>υ)∝ρ-α(1)
式中,A(ρ>υ)表示元素含量大于某一值υ的等值线所围成的面积,随着υ值的提高,A(ρ>υ)总是相应减小。A随υ变化的规律取决于指数α的大小,在背景值和异常值范围,这种变化对应不同的α值。在元素含量的最大值与最小值之间均匀设定15~20个不同的阈值υn,统计含量值大于υn的所有元素所围成的面积An;
在logA和logυ的双对数坐标图上,利用最小二乘法拟合出两段直线,将两段直线的交叉点所对应的元素含量值作为异常下限值;
(5)将步骤(4)中得到的异常下限值作为子区中心位置处的背景值,采用克里格方法对整个测区范围内的子区背景值进行空间插值,得到整个测区的地球化学背景值;
(6)将整个测区的地球化学数据减去步骤(5)中得到的地球化学背景值,得到整个测区的地球化学异常值。
Claims (2)
1.一种地球化学异常圈定方法,其特征在于:包括以下步骤:
(1)根据测区实际地质情况,确定子区的形状:
在沉积盆地地区,子区的形状为矩形;
在南方山区地区,无侵入岩体存在时,子区的形状为圆形;
在南方山区地区,有侵入岩体存在时,子区的形状为椭圆形;
(2)针对步骤(1)中选择的子区形状,确定子区的大小:
在沉积盆地地区,矩形子区的长度≥10km,宽度≥10km;
在南方山区地区,圆的直径长度、椭圆的长轴和短轴长度都∈1~10km;
(3)根据步骤(1)和(2)中确定的子区形状和大小,以子区在横向和纵向上轴线长度的1/3至1/2为横向间距和纵向间距进行滑动,划分测区;
(4)针对每个子区,采用含量与面积分形模型来确定该子区的异常下限值;
含量与面积分形模型由式(1)确定:
A(ρ>υ)∝ρ-α(1)
式中,A(ρ>υ)表示元素含量大于某一值υ的等值线所围成的面积,在元素含量的最大值与最小值之间均匀设定n个不同的阈值υn,统计含量值大于υn的所有元素所围成的面积An;
在logA和logυ的双对数坐标图上,利用最小二乘法拟合出两段直线,将两段直线的交叉点所对应的元素含量值作为异常下限值;
(5)将步骤(4)中得到的异常下限值作为子区中心位置处的背景值,采用克里格方法对整个测区范围内的子区背景值进行空间插值,得到整个测区的地球化学背景值;
(6)将整个测区的地球化学数据减去步骤(5)中得到的地球化学背景值,得到整个测区的地球化学异常值;
划分测区时的具体原则为:
子区形状为矩形时,分别以矩形的长度和宽度的1/3至1/2为横向间距和纵向间距进行滑动,将整个测区划分为若干个子区;
子区的形状为圆形时,分别以圆形直径长度的1/3至1/2为横向间距和纵向间距进行滑动,将整个测区划分为若干个子区;
子区的形状为椭圆形时,分别以椭圆长轴和短轴长度的1/3至1/2为横向间距和纵向间距进行滑动,将整个测区划分为若干个子区。
2.如权利要求1所述的一种地球化学异常圈定方法,其特征在于:步骤(4)中,在元素含量的最大值与最小值之间均匀设定15~20个不同的阈值υn。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310274760.7A CN103529486B (zh) | 2013-07-02 | 2013-07-02 | 一种地球化学异常圈定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310274760.7A CN103529486B (zh) | 2013-07-02 | 2013-07-02 | 一种地球化学异常圈定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103529486A CN103529486A (zh) | 2014-01-22 |
CN103529486B true CN103529486B (zh) | 2016-05-18 |
Family
ID=49931639
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310274760.7A Active CN103529486B (zh) | 2013-07-02 | 2013-07-02 | 一种地球化学异常圈定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103529486B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103869053B (zh) * | 2014-03-24 | 2015-07-15 | 焦振志 | 区域地球化学调查样品分析异常点抽检方法 |
CN106126882B (zh) * | 2016-06-15 | 2019-01-29 | 中国地质大学(北京) | 基于肯德尔和谐系数排序的元素地球化学数据优选方法 |
CN107807221B (zh) * | 2017-09-22 | 2019-12-10 | 中国石油天然气集团公司 | 一种地球化学普查实验室样品分析的异常点抽查检查方法 |
CN108133480A (zh) * | 2017-12-04 | 2018-06-08 | 国家海洋局第海洋研究所 | 一种热液喷口异常地形提取方法 |
CN110633310B (zh) * | 2019-09-12 | 2023-03-21 | 甘肃省地质调查院 | “直线逼近”异常下限确定方法 |
CN112380308B (zh) * | 2020-11-17 | 2021-06-29 | 中国地质科学院矿产资源研究所 | 一种基于数据正则化的地球化学异常圈定方法及系统 |
CN112990322A (zh) * | 2021-03-23 | 2021-06-18 | 中国地质调查局发展研究中心 | 一种确定地球化学异常识别临界值的方法、装置及设备 |
CN113408019B (zh) * | 2021-06-23 | 2022-04-22 | 河北地质大学 | 一种基于bme-gwr的水系沉积物地球化学异常制图方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1716469A1 (ru) * | 1989-07-25 | 1992-02-28 | Производственное Геологическое Объединение "Аэрогеология" | Способ дистанционного зондировани земной поверхности |
CN102254095A (zh) * | 2011-07-01 | 2011-11-23 | 合肥工业大学 | 基于多维分形克里格方法的成矿异常提取方法 |
-
2013
- 2013-07-02 CN CN201310274760.7A patent/CN103529486B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1716469A1 (ru) * | 1989-07-25 | 1992-02-28 | Производственное Геологическое Объединение "Аэрогеология" | Способ дистанционного зондировани земной поверхности |
CN102254095A (zh) * | 2011-07-01 | 2011-11-23 | 合肥工业大学 | 基于多维分形克里格方法的成矿异常提取方法 |
Non-Patent Citations (2)
Title |
---|
基于MAPGIS平台下分形理论在地球化学异常圈定中的应用;郭科等人;《成都理工大学学报(自然科学版)》;20060801;第33卷(第4期);全文 * |
庞西垌地区地球化学异常圈定方法讨论;王琨等人;《地质学刊》;20120201;第36卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103529486A (zh) | 2014-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103529486B (zh) | 一种地球化学异常圈定方法 | |
CN104330828B (zh) | 甜点储层的预测方法及预测装置 | |
Witter et al. | 3D geophysical inversion modeling of gravity data to test the 3D geologic model of the Bradys geothermal area, Nevada, USA | |
CN104695862B (zh) | 一种致密砂岩岩性气藏水平井钻井地质设计方法 | |
CN103824133A (zh) | 一种花岗岩型铀矿田远景区综合预测方法 | |
CN103064124B (zh) | 一种校正电磁勘探中地形影响的比值方法 | |
CN103868503B (zh) | 一种基于insar制作dem的水域高程处理方法 | |
CN109407143A (zh) | 基于压缩感知的地震勘探非规则观测系统设计方法 | |
CN103995281A (zh) | 一种三维观测系统快速滚动布设的方法 | |
CN110286416B (zh) | 一种基于物性函数的快速二维密度反演方法 | |
CN108121009B (zh) | 开发后期复杂断块变速构造成图方法 | |
CN103824510A (zh) | 一种基于Voronoi图的筛选电子地图点状要素的方法 | |
Ni et al. | Lineament length and density analyses based on the segment tracing algorithm: a case study of the gaosong field in gejiu tin mine, China | |
Jonah et al. | A survey for groundwater at a lot at the Dan Zaria Academic Estate, Federal University of Technology, Minna, Central Nigeria | |
CN107092995A (zh) | 地热资源评价方法和装置 | |
CN105116447A (zh) | 一种基于曲率异常条带的地质河道方向判别方法 | |
CN104459800A (zh) | 一种预测砂体尖灭的方法和装置 | |
CN104020508A (zh) | 一种用于地质雷达层析探测的直射线追踪算法 | |
CN104570070A (zh) | 一种建立二维近地表地质模型的方法及设备 | |
CN102736104A (zh) | 近地表表层模型构建方法 | |
CN102914790B (zh) | 二维观测系统和三维观测系统一次采集的观测系统方法 | |
CN106842334A (zh) | 一种电磁勘探方法和装置 | |
Salahuddin et al. | Innovative integration of subsurface data and history matching validation to characterize and model complex carbonate reservoir with high permeability streaks and low resistivity pay issues, onshore Abu Dhabi | |
CN103105629B (zh) | 一种等视电阻率曲面图电磁勘探资料解释方法 | |
Guo et al. | A new method for identifying fractures in tight sandstone of a gentle structural area using well logs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |