CN111045110B - 圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法 - Google Patents

圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法 Download PDF

Info

Publication number
CN111045110B
CN111045110B CN201911299715.0A CN201911299715A CN111045110B CN 111045110 B CN111045110 B CN 111045110B CN 201911299715 A CN201911299715 A CN 201911299715A CN 111045110 B CN111045110 B CN 111045110B
Authority
CN
China
Prior art keywords
uranium
deep
exploration
mineralization
dimensional
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
Application number
CN201911299715.0A
Other languages
English (en)
Other versions
CN111045110A (zh
Inventor
赵丹
蔡煜琦
吴国东
刘武生
易超
乔宝平
张字龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Research Institute of Uranium Geology
Original Assignee
Beijing Research Institute of Uranium Geology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Research Institute of Uranium Geology filed Critical Beijing Research Institute of Uranium Geology
Priority to CN201911299715.0A priority Critical patent/CN111045110B/zh
Publication of CN111045110A publication Critical patent/CN111045110A/zh
Application granted granted Critical
Publication of CN111045110B publication Critical patent/CN111045110B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00—Sampling; Preparing specimens for investigation
    • G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24—Earth materials

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明属于铀矿勘探技术领域，具体涉及一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法。本发明包括以下步骤：步骤1、圈定细粒级土壤铀含量异常范围；步骤2、开展航空瞬变电磁测量工作，圈定铀成矿有利砂体范围；步骤3、圈定铀成矿重点勘查区；步骤4、开展土壤氡气测量工作，圈定土壤氡气浓度异常范围；步骤5、开展可控源音频大地电磁测深工作，圈定深部铀成矿有利砂体范围；步骤6、圈定深部铀成矿有利地段；步骤7、圈定地电化学铀含量异常范围；步骤8、圈定深部三维铀成矿有利砂体范围；步骤9、圈定深部三维铀成矿靶区。本发明能够解决现有技术难以准确地圈定砂岩型铀矿深部三维铀成矿靶区的问题。

Description

圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法
技术领域
本发明属于铀矿勘探技术领域，具体涉及一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法。
背景技术
目前，在铀矿勘查领域中，通过放射性测量方法，埋藏较浅的砂岩型铀矿床已发现殆尽，砂岩型铀矿勘探的重点转入深部勘探，而深部铀矿体在地表的成矿信息比较微弱，需要通过适合的物化探方法和相应的数据处理技术来获取有效的深部铀成矿信息。
为圈定砂岩型铀矿深部铀成矿靶区，锁定深部铀矿体大体空间位置，采取的关键的物化探方法和相应的数据处理技术是本领域技术人员亟需解决的关键问题。
发明内容
本发明解决的技术问题：
本发明提供一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，能够解决现有技术难以准确地圈定砂岩型铀矿深部三维铀成矿靶区的问题。
本发明采用的技术方案：
一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，包括以下步骤：步骤1、在铀成矿勘查区S(U)内开展细粒度土壤化探工作，圈定细粒级土壤铀含量异常范围AS(UFS)；步骤2、在铀成矿勘查区S(U)内开展航空瞬变电磁测量工作，圈定铀成矿有利砂体范围SS(MAEM)；步骤3、圈定铀成矿重点勘查区IS(U)；步骤4、在铀成矿重点勘查区IS(U)内，开展土壤氡气测量工作，圈定土壤氡气浓度异常范围AS(Rn)；步骤5、在铀成矿重点勘查区IS(U)内开展可控源音频大地电磁测深(CSAMT)工作，圈定深部铀成矿有利砂体范围SS(MCSAMT)；步骤6、圈定深部铀成矿有利地段DS(U)；步骤7、在深部铀成矿有利地段DS(U)内开展地电化学测量工作，圈定地电化学铀含量异常范围AS(UEG)；步骤8、在深部铀成矿有利地段DS(U)内开展三维地震勘探工作，圈定深部三维铀成矿有利砂体范围SS(M3DSE)；步骤9、圈定深部三维铀成矿靶区TS(U)。
所述步骤1包括如下步骤，
步骤1.1.细粒级土壤化探数据采集
在铀成矿勘查S(U)区内开展细粒级土壤化探工作，采集浅地表土壤样品，获取各测点样品中的铀含量C(UFS)；
步骤1.2.计算细粒级土壤化探铀含量异常剔除上限值ADH(UFS)和下限值ADL(UFS)
利用下式计算出细粒级土壤化探铀含量异常剔除上限值ADH(UFS)和下限值ADL(UFS)：
Figure BDA0002321550450000031
式中，
Figure BDA0002321550450000032
表示(m，n)测点的原始细粒级土壤化探铀含量，单位为×10-6g/g；
Figure BDA0002321550450000033
表示所有测点的
Figure BDA0002321550450000034
的平均值；
Figure BDA0002321550450000035
表示所有测点的
Figure BDA0002321550450000036
的标准偏差；
i表示测线的数目；
j表示每条测线上的测点数目；
m为整数，取值为1，2，3，……，i；
n为整数，取值为1，2，3，……，j；
Figure BDA0002321550450000037
表示异常剔除上限值，单位为×10-6g/g；
Figure BDA0002321550450000038
表示异常剔除下限值，单位为×10-6g/g；
步骤1.3.替换细粒级土壤化探铀含量数据中的异常高值和低值
将原始细粒级土壤化探铀含量数据中大于
Figure BDA0002321550450000039
的值用ATH(UFS)代替；将小于
Figure BDA00023215504500000310
的数据值用ATL(UFS)代替。
步骤1.4.计算细粒级土壤化探铀含量异常下限值AT(UFS)
基于步骤1.3.后的全部数据，利用下式计算细粒级土壤化探铀含量异常下限值AT(UFS)：
Figure BDA0002321550450000041
式中，
Figure BDA0002321550450000042
表示(m，n)测点的经替换处理后的细粒级土壤化探铀含量，单位为×10- 6g/g；
Figure BDA0002321550450000043
表示所有测点的
Figure BDA0002321550450000044
的平均值；
Figure BDA0002321550450000045
表示所有测点的
Figure BDA0002321550450000046
的标准偏差；
AT(UFS)表示细粒级土壤化探铀含量异常下限值，单位为×10-6g/g；
i表示测线的数目；
j表示每条测线上的测点数目；
m为整数，取值为1，2，3，……，i；
n为整数，取值为1，2，3，……，j；
步骤1.5.圈定细粒级土壤化探铀异常范围AS(UFS)
采用“径向基函数插值方法”对铀成矿远景区内的
Figure BDA0002321550450000047
进行插值计算得到网格化数据，对网格化数据进行“低通滤波”处理，圈定大于等于AT(UFS)的数据范围，定义该范围为AS(UFS)。
所述步骤2包括如下步骤，
步骤2.1.航空瞬变电磁测量数据采集；
步骤2.2.获取航空瞬变电磁反演视电阻率断面图
采用奥康姆一维反演方法对ATEM数据进行反演计算，得到各条ATEM测线的电阻率断面图；
步骤2.3.圈定铀成矿有利砂体范围SS(MATEM)
根据ATEM测线的电阻率断面图，结合铀成矿勘查区S(U)内砂体电阻率特征和铀成矿地质特征，在平面上圈定成矿有利砂体范围SS(MATEM)。
所述步骤3中，将AS(UFS)与SS(MATEM)重叠的区域圈定为铀成矿重点勘查区IS(U)。
所述步骤4包括如下步骤，
步骤4.1.土壤氡气测量数据采集
在铀成矿重点勘查区IS(U)范围内获取测点处土壤中氡气浓度值C(Rn)；
步骤4.2.计算土壤氡气浓度剔除下限值T(Rn)
利用下式计算出氡浓度剔除下限值T(Rn)：
Figure BDA0002321550450000051
式中，
Cmn(Rn)表示测点(p,q)处的土壤氡气浓度，单位为Bq/m3；
M(Rn)表示所有测点的ln[Cpq(Rn)]的平均值；
S(Rn)表示所有测点的ln[Cpq(Rn)]的标准偏差；
T(Rn)表示土壤氡气浓度剔除下限值，单位为Bq/m3；
k表示测线的数目；
l表示每条测线上的测点数目；
p为整数，取值为1，2，3，……，k；
q为整数，取值为1，2，3，……，l；
步骤4.3.计算土壤氡气异常下限值A(Rn)
首先，采用“迭代剔除法”，逐次剔除大于等于T(Rn)的土壤氡气浓度数据，直到未有大于等于T(Rn)的氡浓度数据为止；然后，利用下式计算出剩余土壤氡气浓度数据的平均值Mr(Rn)、标准偏差Sr(Rn)和异常下限A(Rn)：
Figure BDA0002321550450000061
式中，
Cr(Rn)表示剔除后剩余测点的氡浓度，单位为Bq/m3；
Mr(Rn)表示剔除后剩余测点的ln[Cw(Rn)]的平均值；
Sr(Rn)表示剔除后剩余测点的ln[Cw(Rn)]的标准偏差；
t表示剔除后剩余测点的数目；
w为整数，取值为1，2，3，……，t；
A(Rn)表示氡浓度异常下限值，单位为Bq/m3；
步骤4.4.圈定土壤氡气浓度异常范围AS(Rn)
采用“反距离加权插值法”对测区的氡气浓度值C(Rn)进行插值计算得到网格化数据，然后圈定大于等于A(Rn)的数据范围，定义该范围为AS(Rn)。
所述步骤5包括如下步骤，
步骤5.1.可控源音频大地电磁测量数据采集
在铀成矿重点勘查区IS(U)范围内，开展可控源音频大地电磁测量工作；
步骤5.2.获取可控源音频大地电磁测深反演视电阻率断面图
采用“非线性共轭梯度”二维反演方法对可控源音频大地电磁数据进行反演计算，得到各条可控源音频大地电磁测线的电阻率断面图；
步骤5.3.圈定视电阻率断面图上深部铀成矿有利砂体范围SS(MCSAMT)
根据可控源音频大地电磁测线的电阻率断面图，结合IS(U)范围内铀成矿有利砂体电阻率特征和铀成矿地质特征，圈定深部铀成矿有利砂体范围SS(MCSAMT)。
所述步骤6中，将AS(Rn)与SS(MCSAMT)重叠的区域圈定为深部铀成矿有利地段DS(U)。
所述步骤7包括如下步骤，
步骤7.1.地电化学测量数据采集
在深部铀成矿有利地段DS(U)范围内，获取测点处土壤中地电化学铀含量值C(UEG)；
步骤7.2.计算地电化学铀含量异常剔除下限值T(UEG)
利用下式计算地电化学铀含量异常剔除下限值T(UEG)：
Figure BDA0002321550450000081
式中，
Cxy(UEG)表示测点(x，y)处的地电化学铀含量值，单位为×10-9g/g；
M(UEG)表示所有测点的ln[Cxy(UEG)]的平均值；
SD(UEG)表示所有测点的ln[Cxy(UEG)]的标准偏差；
T(UEG)表示地电化学铀含量剔除下限值，单位为×10-9g/g；
e表示测线的数目；
f表示每条测线上的测点数目；
x为整数，取值为1，2，3，……，e；
y为整数，取值为1，2，3，……，f；
步骤7.3.计算地电化学铀含量异常下限A(UEG)
采用“迭代剔除法”，逐次剔除大于等于T(UEG)的地电化学铀含量数据，直到未有大于等于T(UEG)的地电化学铀含量数据为止，然后，利用下式计算出剩余地电化学铀含量数据的平均值Mr(UEG)、标准偏差Sr(UEG)和异常下限A(UEG)：
Figure BDA0002321550450000091
式中，
Cu(UEG)表示剔除后剩余测点的地电化学铀含量值，单位为×10-9g/g；
Mr(UEG)表示剔除后剩余测点的ln[Cu(UEG)]的平均值；
Sr(UEG)表示剔除后剩余测点的ln[Cu(UEG)]的标准偏差；
w表示剔除后剩余测点的数目；
u为整数，取值为1，2，3，……，w；
A(UEG)表示地电化学铀异常下限值，单位为×10-9g/g；
步骤7.4.圈定地电化学铀异常范围AS(UEG)
采用“最小曲率插值法”对测区的地电化学铀含量值C(UEG)进行插值计算得到网格化数据，然后圈定大于等于A(UEG)的数据范围，定义该范围为AS(UEG)。
所述步骤8包括如下步骤，
步骤8.1.三维地震勘探数据采集
在深部铀成矿有利地段DS(U)范围内，开展三维地震勘探工作；
步骤8.2.获取三维地震勘探数据体
基于步骤8.1中采集的三维地震勘探数据，进行数据处理，获取三维地震勘探数据体；
步骤8.3.圈定深部三维铀成矿有利砂体范围SS(M3DSE)
基于步骤8.2中的三维地震勘探数据体，采用“拟声波地震反演”方法识别砂体，结合DS(U)范围内铀成矿砂岩层展布特征和铀成矿地质特征，圈定深部三维铀成矿有利砂体范围SS(M3DSE)。
所述步骤9中，将AS(UEG)与SS(M3DSE)重叠的区域圈定为深部三维铀成矿靶区TS(U)。
本发明的有益效果：
(1)本发明提供的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，能够快速、准确地圈定砂岩型铀矿深部三维铀成矿靶区，解决了砂岩型铀矿深部铀矿体赋存位置难以圈定的难题，切实保障了砂岩型铀矿深部勘查的技术可行性；
(2)本发明提供的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，为砂岩型铀矿深部勘查的钻孔布设提供了可靠依据，能够大量节约用于深部钻探前期勘查所需的人力、物力、财力；
(3)本发明提供的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，有效提升了深部砂岩型铀矿勘查的经济性，为我国深部砂岩型铀矿勘查提供了经济可行的技术手段，从而有效地助力于我国军用和民用铀资源保障。
附图说明
图1为本发明提供的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法流程图。
具体实施方式
下面结合附图和具体实施例对本发明的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法进行详细说明。
如图1所示，本发明的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，包括以下步骤：
步骤1.圈定细粒级土壤铀含量异常范围AS(UFS)
步骤1.1.细粒级土壤化探数据采集
在铀成矿勘查S(U)区内开展测网为1千米×1千米的细粒级土壤化探工作，采集浅地表土壤样品，晾干后筛至-400目，过筛后每个样品质量不小于10克，送至分析测试机构，获取各测点样品中的铀含量C(UFS)。
步骤1.2.计算细粒级土壤化探铀含量异常剔除上限值ADH(UFS)和下限值ADL(UFS)
利用下式计算出细粒级土壤化探铀含量异常剔除上限值ADH(UFS)和下限值ADL(UFS)：
Figure BDA0002321550450000111
式中，
Figure BDA0002321550450000121
表示(m，n)测点的原始细粒级土壤化探铀含量，单位为×10-6g/g；
Figure BDA0002321550450000122
表示所有测点的
Figure BDA0002321550450000123
的平均值；
Figure BDA0002321550450000124
表示所有测点的
Figure BDA0002321550450000125
的标准偏差；
i表示测线的数目；
j表示每条测线上的测点数目；
m为整数，取值为1，2，3，……，i；
n为整数，取值为1，2，3，……，j；
Figure BDA0002321550450000126
表示异常剔除上限值，单位为×10-6g/g；
Figure BDA0002321550450000127
表示异常剔除下限值，单位为×10-6g/g。
步骤1.3.替换细粒级土壤化探铀含量数据中的异常高值和低值
首先，将原始细粒级土壤化探铀含量数据中大于
Figure BDA0002321550450000128
的值用ATH(UFS)代替；然后，将小于
Figure BDA0002321550450000129
的数据值用ATL(UFS)代替。
步骤1.4.计算细粒级土壤化探铀含量异常下限值AT(UFS)
基于步骤1.3.后的全部数据，利用下式计算细粒级土壤化探铀含量异常下限值AT(UFS)：
Figure BDA00023215504500001210
式中，
Figure BDA0002321550450000131
表示(m，n)测点的经替换处理后的细粒级土壤化探铀含量，单位为×10- 6g/g；
Figure BDA0002321550450000132
表示所有测点的
Figure BDA0002321550450000133
的平均值；
Figure BDA0002321550450000134
表示所有测点的
Figure BDA0002321550450000135
的标准偏差；
AT(UFS)表示细粒级土壤化探铀含量异常下限值，单位为×10-6g/g。
i表示测线的数目；
j表示每条测线上的测点数目；
m为整数，取值为1，2，3，……，i；
n为整数，取值为1，2，3，……，j。
步骤1.5.圈定细粒级土壤化探铀异常范围AS(UFS)
采用“径向基函数插值方法”对铀成矿远景区内的
Figure BDA0002321550450000136
进行插值计算得到网格化数据，然后对网格化数据进行“低通滤波”处理，然后圈定大于等于AT(UFS)的数据范围，定义该范围为AS(UFS)。
步骤2.圈定铀成矿有利砂体范围SS(MATEM)
步骤2.1.航空瞬变电磁(ATEM)测量数据采集
在铀成矿勘查区S(U)内开展航空瞬变电磁(ATEM)测量工作，采用固定翼机载测量方式，测线布设方向垂直于S(U)区内主要构造方向，飞行线距500米，飞行高度≤100米，飞行速度在80千米/小时至100千米/小时之间。
步骤2.2.获取航空瞬变电磁(ATEM)反演视电阻率断面图
采用奥康姆(Occam)一维反演方法对ATEM数据进行反演计算，得到各条ATEM测线的电阻率断面图。
步骤2.3.圈定铀成矿有利砂体范围SS(MATEM)
根据ATEM测线的电阻率断面图，结合铀成矿勘查区S(U)内砂体电阻率特征和铀成矿地质特征，在平面上圈定成矿有利砂体范围SS(MATEM)。
步骤3.圈定铀成矿重点勘查区IS(U)
将AS(UFS)与SS(MATEM)重叠的区域圈定为铀成矿重点勘查区IS(U)。
步骤4.圈定土壤氡气浓度异常范围AS(Rn)
步骤4.1.土壤氡气测量数据采集
在铀成矿重点勘查区IS(U)范围内开展测网为100米×100米的土壤氡气测量工作，采用基于金硅面垒型探测器的α能谱仪测量土壤中的氡气浓度，取气深度为1米，获取测点处土壤中氡气浓度值C(Rn)。
步骤4.2.计算土壤氡气浓度剔除下限值T(Rn)
利用下式计算出氡浓度剔除下限值T(Rn)：
Figure BDA0002321550450000141
式中，
Cmn(Rn)表示测点(p,q)处的土壤氡气浓度，单位为Bq/m3；
M(Rn)表示所有测点的ln[Cpq(Rn)]的平均值；
S(Rn)表示所有测点的ln[Cpq(Rn)]的标准偏差；
T(Rn)表示土壤氡气浓度剔除下限值，单位为Bq/m3；
k表示测线的数目；
l表示每条测线上的测点数目；
p为整数，取值为1，2，3，……，k；
q为整数，取值为1，2，3，……，l。
步骤4.3.计算土壤氡气异常下限值A(Rn)
首先，采用“迭代剔除法”，逐次剔除大于等于T(Rn)的土壤氡气浓度数据，直到未有大于等于T(Rn)的氡浓度数据为止；然后，利用下式计算出剩余土壤氡气浓度数据的平均值Mr(Rn)、标准偏差Sr(Rn)和异常下限A(Rn)：
Figure BDA0002321550450000151
式中，
Cr(Rn)表示剔除后剩余测点的氡浓度，单位为Bq/m3；
Mr(Rn)表示剔除后剩余测点的ln[Cw(Rn)]的平均值；
Sr(Rn)表示剔除后剩余测点的ln[Cw(Rn)]的标准偏差；
t表示剔除后剩余测点的数目；
w为整数，取值为1，2，3，……，t；
A(Rn)表示氡浓度异常下限值，单位为Bq/m3。
步骤4.4.圈定土壤氡气浓度异常范围AS(Rn)
采用“反距离加权插值法”对测区的氡气浓度值C(Rn)进行插值计算得到网格化数据，然后圈定大于等于A(Rn)的数据范围，定义该范围为AS(Rn)。
步骤5.圈定深部铀成矿有利砂体范围SS(MCSAMT)
步骤5.1.可控源音频大地电磁(CSAMT)测量数据采集
在铀成矿重点勘查区IS(U)范围内，开展可控源音频大地电磁(CSAMT)测量工作，线距200米，点距50米，测线布设方向垂直于SI(U)范围内主要断裂构造方向。
步骤5.2.获取可控源音频大地电磁测深(CSAMT)反演视电阻率断面图
采用“非线性共轭梯度”二维反演方法对CSAMT数据进行反演计算，得到各条CSAMT测线的电阻率断面图。
步骤5.3.圈定视电阻率断面图上深部铀成矿有利砂体范围SS(MCSAMT)
根据CSAMT测线的电阻率断面图，结合IS(U)范围内铀成矿有利砂体电阻率特征和铀成矿地质特征，圈定深部铀成矿有利砂体范围SS(MCSAMT)。
步骤6.圈定深部铀成矿有利地段DS(U)
将AS(Rn)与SS(MCSAMT)重叠的区域圈定为深部铀成矿有利地段DS(U)。
步骤7.圈定地电化学铀含量异常范围AS(UEG)
步骤7.1.地电化学测量数据采集
在深部铀成矿有利地段DS(U)范围内，开展测网为20米×20米的地电化学测量工作，获取测点处土壤中地电化学铀含量值C(UEG)。
步骤7.2.计算地电化学铀含量异常剔除下限值T(UEG)
利用下式计算地电化学铀含量异常剔除下限值T(UEG)：
Figure BDA0002321550450000171
式中，
Cxy(UEG)表示测点(x，y)处的地电化学铀含量值，单位为×10-9g/g；
M(UEG)表示所有测点的ln[Cxy(UEG)]的平均值；
SD(UEG)表示所有测点的ln[Cxy(UEG)]的标准偏差；
T(UEG)表示地电化学铀含量剔除下限值，单位为×10-9g/g；
e表示测线的数目；
f表示每条测线上的测点数目；
x为整数，取值为1，2，3，……，e；
y为整数，取值为1，2，3，……，f。
步骤7.3.计算地电化学铀含量异常下限A(UEG)
采用“迭代剔除法”，逐次剔除大于等于T(UEG)的地电化学铀含量数据，直到未有大于等于T(UEG)的地电化学铀含量数据为止，然后，利用下式计算出剩余地电化学铀含量数据的平均值Mr(UEG)、标准偏差Sr(UEG)和异常下限A(UEG)：
Figure BDA0002321550450000181
式中，
Cu(UEG)表示剔除后剩余测点的地电化学铀含量值，单位为×10-9g/g；
Mr(UEG)表示剔除后剩余测点的ln[Cu(UEG)]的平均值；
Sr(UEG)表示剔除后剩余测点的ln[Cu(UEG)]的标准偏差；
w表示剔除后剩余测点的数目；
u为整数，取值为1，2，3，……，w；
A(UEG)表示地电化学铀异常下限值，单位为×10-9g/g。
步骤7.4.圈定地电化学铀异常范围AS(UEG)
采用“最小曲率插值法”对测区的地电化学铀含量值C(UEG)进行插值计算得到网格化数据，然后圈定大于等于A(UEG)的数据范围，定义该范围为AS(UEG)。
步骤8.圈定深部三维铀成矿有利砂体范围SS(M3DSE)
步骤8.1.三维地震勘探数据采集
在深部铀成矿有利地段DS(U)范围内，开展三维地震勘探工作，观测系统参数为16L×4S×128R，道间距为10米，炮点距为20米，接收线距为80米，炮线距为40米，纵向炮检距为635米-5-10-10-635米，震源采用可控震源车。
步骤8.2.获取三维地震勘探数据体
基于步骤8.1中采集的三维地震勘探数据，进行数据处理，获取三维地震勘探数据体。
步骤8.3.圈定深部三维铀成矿有利砂体范围SS(M3DSE)
基于步骤8.2中的三维地震勘探数据体，采用“拟声波地震反演”方法识别砂体，结合DS(U)范围内铀成矿砂岩层展布特征和铀成矿地质特征，圈定深部三维铀成矿有利砂体范围SS(M3DSE)。
步骤9.圈定深部三维铀成矿靶区TS(U)
将AS(UEG)与SS(M3DSE)重叠的区域圈定为深部三维铀成矿靶区TS(U)。

Claims (10)

1.一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：包括以下步骤：步骤(1)、在铀成矿勘查区S(U)内开展细粒度土壤化探工作，圈定细粒级土壤铀含量异常范围AS(UFS)；步骤(2)、在铀成矿勘查区S(U)内开展航空瞬变电磁测量工作，圈定铀成矿有利砂体范围SS(MAEM)；步骤(3)、圈定铀成矿重点勘查区IS(U)；步骤(4)、在铀成矿重点勘查区IS(U)内，开展土壤氡气测量工作，圈定土壤氡气浓度异常范围AS(Rn)；步骤(5)、在铀成矿重点勘查区IS(U)内开展可控源音频大地电磁测深工作，圈定深部铀成矿有利砂体范围SS(MCSAMT)；步骤(6)、圈定深部铀成矿有利地段DS(U)；步骤(7)、在深部铀成矿有利地段DS(U)内开展地电化学测量工作，圈定地电化学铀含量异常范围AS(UEG)；步骤(8)、在深部铀成矿有利地段DS(U)内开展三维地震勘探工作，圈定深部三维铀成矿有利砂体范围SS(M3DSE)；步骤(9)、圈定深部三维铀成矿靶区TS(U)。
2.根据权利要求1所述的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：所述步骤(1)包括如下步骤，
步骤(1.1).细粒级土壤化探数据采集
在铀成矿勘查S(U)区内开展细粒级土壤化探工作，采集浅地表土壤样品，获取各测点样品中的铀含量C(UFS)；
步骤(1.2).计算细粒级土壤化探铀含量异常剔除上限值ADH(UFS)和下限值ADL(UFS)
利用下式计算出细粒级土壤化探铀含量异常剔除上限值ADH(UFS)和下限值ADL(UFS)：
Figure FDA0003512818330000021
式中，
Figure FDA0003512818330000022
表示(m，n)测点的原始细粒级土壤化探铀含量，单位为×10-6g/g；
Figure FDA0003512818330000023
表示所有测点的
Figure FDA0003512818330000024
的平均值；
Figure FDA0003512818330000025
表示所有测点的
Figure FDA0003512818330000026
的标准偏差；
i表示测线的数目；
j表示每条测线上的测点数目；
m为整数，取值为1，2，3，……，i；
n为整数，取值为1，2，3，……，j；
Figure FDA0003512818330000027
表示异常剔除上限值，单位为×10-6g/g；
Figure FDA0003512818330000028
表示异常剔除下限值，单位为×10-6g/g；
步骤(1.3).替换细粒级土壤化探铀含量数据中的异常高值和低值
将原始细粒级土壤化探铀含量数据中大于
Figure FDA0003512818330000029
的值用ATH(UFS)代替；将小于
Figure FDA00035128183300000210
的数据值用ATL(UFS)代替；
步骤(1.4).计算细粒级土壤化探铀含量异常下限值AT(UFS)
基于步骤(1.3).后的全部数据，利用下式计算细粒级土壤化探铀含量异常下限值AT(UFS)：
Figure FDA0003512818330000031
式中，
Figure FDA0003512818330000032
表示(m，n)测点的经替换处理后的细粒级土壤化探铀含量，单位为×10-6g/g；
Figure FDA0003512818330000033
表示所有测点的
Figure FDA0003512818330000034
的平均值；
Figure FDA0003512818330000035
表示所有测点的
Figure FDA0003512818330000036
的标准偏差；
AT(UFS)表示细粒级土壤化探铀含量异常下限值，单位为×10-6g/g；
i表示测线的数目；
j表示每条测线上的测点数目；
m为整数，取值为1，2，3，……，i；
n为整数，取值为1，2，3，……，j；
步骤(1.5).圈定细粒级土壤化探铀异常范围AS(UFS)
采用“径向基函数插值方法”对铀成矿远景区内的
Figure FDA0003512818330000037
进行插值计算得到网格化数据，对网格化数据进行“低通滤波”处理，圈定大于等于AT(UFS)的数据范围，定义该范围为AS(UFS)。
3.根据权利要求2所述的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：所述步骤(2)包括如下步骤，
步骤(2.1).航空瞬变电磁测量数据采集；
步骤(2.2).获取航空瞬变电磁反演视电阻率断面图
采用奥康姆一维反演方法对ATEM数据进行反演计算，得到各条ATEM测线的电阻率断面图；
步骤(2.3).圈定铀成矿有利砂体范围SS(MATEM)
根据ATEM测线的电阻率断面图，结合铀成矿勘查区S(U)内砂体电阻率特征和铀成矿地质特征，在平面上圈定成矿有利砂体范围SS(MATEM)。
4.根据权利要求3所述的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：所述步骤(3)中，将AS(UFS)与SS(MATEM)重叠的区域圈定为铀成矿重点勘查区IS(U)。
5.根据权利要求4所述的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：所述步骤(4)包括如下步骤，
步骤(4.1).土壤氡气测量数据采集
在铀成矿重点勘查区IS(U)范围内获取测点处土壤中氡气浓度值C(Rn)；
步骤(4.2).计算土壤氡气浓度剔除下限值T(Rn)
利用下式计算出氡浓度剔除下限值T(Rn)：
Figure FDA0003512818330000051
式中，
Cpq(Rn)表示测点(p,q)处的土壤氡气浓度，单位为Bq/m3；
M(Rn)表示所有测点的ln[Cpq(Rn)]的平均值；
S(Rn)表示所有测点的ln[Cpq(Rn)]的标准偏差；
T(Rn)表示土壤氡气浓度剔除下限值，单位为Bq/m3；
k表示测线的数目；
l表示每条测线上的测点数目；
p为整数，取值为1，2，3，……，k；
q为整数，取值为1，2，3，……，l；
步骤(4.3).计算土壤氡气异常下限值A(Rn)
首先，采用“迭代剔除法”，逐次剔除大于等于T(Rn)的土壤氡气浓度数据，直到未有大于等于T(Rn)的氡浓度数据为止；然后，利用下式计算出剩余土壤氡气浓度数据的平均值Mr(Rn)、标准偏差Sr(Rn)和异常下限A(Rn)：
Figure FDA0003512818330000061
式中，
Cr(Rn)表示剔除后剩余测点的氡浓度，单位为Bq/m3；
Mr(Rn)表示剔除后剩余测点的ln[Cr(Rn)]的平均值；
Sr(Rn)表示剔除后剩余测点的ln[Cr(Rn)]的标准偏差；
t表示剔除后剩余测点的数目；
r为整数，取值为1，2，3，……，t；
A(Rn)表示氡浓度异常下限值，单位为Bq/m3；
步骤(4.4).圈定土壤氡气浓度异常范围AS(Rn)
采用“反距离加权插值法”对测区的氡气浓度值C(Rn)进行插值计算得到网格化数据，然后圈定大于等于A(Rn)的数据范围，定义该范围为AS(Rn)。
6.根据权利要求5所述的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：所述步骤(5)包括如下步骤，
步骤(5.1).可控源音频大地电磁测量数据采集
在铀成矿重点勘查区IS(U)范围内，开展可控源音频大地电磁测量工作；
步骤(5.2).获取可控源音频大地电磁测深反演视电阻率断面图
采用“非线性共轭梯度”二维反演方法对可控源音频大地电磁数据进行反演计算，得到各条可控源音频大地电磁测线的电阻率断面图；
步骤(5.3).圈定视电阻率断面图上深部铀成矿有利砂体范围SS(MCSAMT)
根据可控源音频大地电磁测线的电阻率断面图，结合IS(U)范围内铀成矿有利砂体电阻率特征和铀成矿地质特征，圈定深部铀成矿有利砂体范围SS(MCSAMT)。
7.根据权利要求6所述的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：所述步骤(6)中，将AS(Rn)与SS(MCSAMT)重叠的区域圈定为深部铀成矿有利地段DS(U)。
8.根据权利要求7所述的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：所述步骤(7)包括如下步骤，
步骤(7.1).地电化学测量数据采集
在深部铀成矿有利地段DS(U)范围内，获取测点处土壤中地电化学铀含量值C(UEG)；
步骤(7.2).计算地电化学铀含量异常剔除下限值T(UEG)
利用下式计算地电化学铀含量异常剔除下限值T(UEG)：
Figure FDA0003512818330000071
式中，
Cxy(UEG)表示测点(x，y)处的地电化学铀含量值，单位为×10-9g/g；
M(UEG)表示所有测点的ln[Cxy(UEG)]的平均值；
SD(UEG)表示所有测点的ln[Cxy(UEG)]的标准偏差；
T(UEG)表示地电化学铀含量剔除下限值，单位为×10-9g/g；
e表示测线的数目；
f表示每条测线上的测点数目；
x为整数，取值为1，2，3，……，e；
y为整数，取值为1，2，3，……，f；
步骤(7.3).计算地电化学铀含量异常下限A(UEG)
采用“迭代剔除法”，逐次剔除大于等于T(UEG)的地电化学铀含量数据，直到未有大于等于T(UEG)的地电化学铀含量数据为止，然后，利用下式计算出剩余地电化学铀含量数据的平均值Mr(UEG)、标准偏差Sr(UEG)和异常下限A(UEG)：
Figure FDA0003512818330000081
式中，
Cu(UEG)表示剔除后剩余测点的地电化学铀含量值，单位为×10-9g/g；
Mr(UEG)表示剔除后剩余测点的ln[Cu(UEG)]的平均值；
Sr(UEG)表示剔除后剩余测点的ln[Cu(UEG)]的标准偏差；
w表示剔除后剩余测点的数目；
u为整数，取值为1，2，3，……，w；
A(UEG)表示地电化学铀异常下限值，单位为×10-9g/g；
步骤(7.4).圈定地电化学铀异常范围AS(UEG)
采用“最小曲率插值法”对测区的地电化学铀含量值C(UEG)进行插值计算得到网格化数据，然后圈定大于等于A(UEG)的数据范围，定义该范围为AS(UEG)。
9.根据权利要求8所述的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：所述步骤(8)包括如下步骤，
步骤(8.1).三维地震勘探数据采集
在深部铀成矿有利地段DS(U)范围内，开展三维地震勘探工作；
步骤(8.2).获取三维地震勘探数据体
基于步骤(8.1)中采集的三维地震勘探数据，进行数据处理，获取三维地震勘探数据体；
步骤(8.3).圈定深部三维铀成矿有利砂体范围SS(M3DSE)
基于步骤(8.2)中的三维地震勘探数据体，采用“拟声波地震反演”方法识别砂体，结合DS(U)范围内铀成矿砂岩层展布特征和铀成矿地质特征，圈定深部三维铀成矿有利砂体范围SS(M3DSE)。
10.根据权利要求9所述的一种圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法，其特征在于：所述步骤(9)中，将AS(UEG)与SS(M3DSE)重叠的区域圈定为深部三维铀成矿靶区TS(U)。
CN201911299715.0A 2019-12-17 2019-12-17 圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法 Active CN111045110B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911299715.0A CN111045110B (zh) 2019-12-17 2019-12-17 圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911299715.0A CN111045110B (zh) 2019-12-17 2019-12-17 圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法

Publications (2)

Publication Number Publication Date
CN111045110A CN111045110A (zh) 2020-04-21
CN111045110B true CN111045110B (zh) 2022-04-19

Family

ID=70235128

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911299715.0A Active CN111045110B (zh) 2019-12-17 2019-12-17 圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法

Country Status (1)

Country Link
CN (1) CN111045110B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111505734B (zh) * 2020-05-19 2022-08-05 核工业二三0研究所 一种基于电磁声子理论的深部铀资源探测技术方法
CN111681124B (zh) * 2020-07-06 2022-06-03 核工业北京地质研究院 一种深部砂岩型铀矿化信息三维氡异常识别方法及系统
CN111679342B (zh) * 2020-07-21 2023-04-28 核工业航测遥感中心 快速寻找盆地隐伏砂岩型铀矿方法
CN112731557B (zh) * 2020-11-27 2023-04-14 核工业北京地质研究院 一种预测砂岩铀矿成矿有利区的综合地球物理方法
CN112764123B (zh) * 2020-12-30 2021-11-30 核工业北京地质研究院 一种火山岩型铀矿三维空间定位预测方法
CN112882124B (zh) * 2021-01-13 2024-02-20 核工业北京地质研究院 一种三维铀成矿有利部位的圈定方法
CN117275601B (zh) * 2023-11-17 2024-02-20 核工业北京地质研究院 砂岩型铀矿异常信息的确定方法
CN117723579B (zh) * 2024-02-07 2024-05-28 核工业北京地质研究院 通过矿物组合确定勘测区的砂岩型铀矿矿化部位的方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4081675A (en) * 1976-11-08 1978-03-28 Phillips Petroleum Company Geophysical and geochemical exploration
SE8105923L (sv) * 1981-10-07 1983-04-08 Boliden Ab Sett att indikera dolda fyndigheter
CN103678778B (zh) * 2013-11-22 2017-02-01 核工业北京地质研究院 一种放射性物化探信息综合方法
CN103852807A (zh) * 2014-03-20 2014-06-11 青海省地质矿产研究所 一种高海拔冻土区斑岩型铜多金属矿勘查技术组合方法
CN105549113B (zh) * 2014-10-31 2018-03-16 核工业北京地质研究院 火山岩型铀矿田科学深钻场址遴选方法
CN108008456B (zh) * 2016-10-27 2019-08-13 核工业北京地质研究院 一种圈定热液型铀矿深部三维重点铀成矿有利靶区的方法
CN109270588B (zh) * 2018-09-30 2020-02-21 核工业北京地质研究院 一种花岗岩型铀矿深部矿化信息提取方法

Also Published As

Publication number Publication date
CN111045110A (zh) 2020-04-21

Similar Documents

Publication Publication Date Title
CN111045110B (zh) 圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法
CN111679342B (zh) 快速寻找盆地隐伏砂岩型铀矿方法
CN103913776A (zh) 盾构施工中球状孤石的探测方法
CN109782360A (zh) 一种低阻覆盖区深部砂体探测方法
CN105093314B (zh) 一种测定微地震震源的方法
Bryman et al. Muon geotomography—bringing new physics to orebody imaging
CN107807413A (zh) 地质模型确定方法和装置
CN102914797B (zh) 一种获得地层各向异性系数的方法及装置
CN108008456A (zh) 一种圈定热液型铀矿深部三维重点铀成矿有利靶区的方法
Wang et al. Coal mine goaf interpretation: survey, passive electromagnetic methods and case study
CN112782773B (zh) 一种侵入岩型铀钍铌钽矿隐伏资源预测评价方法
Takahashi ISRM suggested methods for land geophysics in rock engineering
Essa et al. Magnetic inversion approach for modeling data acquired across faults: various environmental cases studies
Asfahani Multifractal approach for delineating uranium anomalies related to phosphatic deposits in Area-3, Northern Palmyrides, Syria
Alumbaugh et al. Monitoring infiltration within the vadose zone using cross borehole ground penetrating radar
CN112882124B (zh) 一种三维铀成矿有利部位的圈定方法
CN115327663A (zh) 深部矿产资源勘查用空-地-井立体地球物理探测方法
CN107797160A (zh) 弹性波和电磁波ct勘测数据联合分析系统及方法
Zhang et al. Thickness Identification of Tunnel Lining Structure by Time-Energy Density Analysis based on Wavelet Transform.
Barale et al. Estimation of natural asbestos content in rocks by fracture network modeling and petrographic characterization
Wang et al. The application of aerogeophysical prospecting technology in the investigation of ancient water conservancy projects—a case analysis of Li Bing's Gaojingguan, Shifang water conservancy project
Liang et al. A review on sustainable development of geological exploration technology and risk management
CN114019570B (zh) 一种天地孔协同一体的瞬变电磁探测采空区的方法
CN113850770B (zh) 一种盾构隧道和隔震隧道的影像ai数字化检测方法及装置
Airo et al. Aerogeophysical approach to ductile and brittle structures in the densely populated urban Helsinki area, southern Finland

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