CN111045110B - 圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法 - Google Patents
圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法 Download PDFInfo
- 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
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)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (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):
式中,
i表示测线的数目;
j表示每条测线上的测点数目;
m为整数,取值为1,2,3,……,i;
n为整数,取值为1,2,3,……,j;
步骤1.3.替换细粒级土壤化探铀含量数据中的异常高值和低值
步骤1.4.计算细粒级土壤化探铀含量异常下限值AT(UFS)
基于步骤1.3.后的全部数据,利用下式计算细粒级土壤化探铀含量异常下限值AT(UFS):
式中,
AT(UFS)表示细粒级土壤化探铀含量异常下限值,单位为×10-6g/g;
i表示测线的数目;
j表示每条测线上的测点数目;
m为整数,取值为1,2,3,……,i;
n为整数,取值为1,2,3,……,j;
步骤1.5.圈定细粒级土壤化探铀异常范围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):
式中,
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):
式中,
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):
式中,
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):
式中,
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):
式中,
i表示测线的数目;
j表示每条测线上的测点数目;
m为整数,取值为1,2,3,……,i;
n为整数,取值为1,2,3,……,j;
步骤1.3.替换细粒级土壤化探铀含量数据中的异常高值和低值
步骤1.4.计算细粒级土壤化探铀含量异常下限值AT(UFS)
基于步骤1.3.后的全部数据,利用下式计算细粒级土壤化探铀含量异常下限值AT(UFS):
式中,
AT(UFS)表示细粒级土壤化探铀含量异常下限值,单位为×10-6g/g。
i表示测线的数目;
j表示每条测线上的测点数目;
m为整数,取值为1,2,3,……,i;
n为整数,取值为1,2,3,……,j。
步骤1.5.圈定细粒级土壤化探铀异常范围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):
式中,
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):
式中,
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):
式中,
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):
式中,
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):
式中,
i表示测线的数目;
j表示每条测线上的测点数目;
m为整数,取值为1,2,3,……,i;
n为整数,取值为1,2,3,……,j;
步骤(1.3).替换细粒级土壤化探铀含量数据中的异常高值和低值
步骤(1.4).计算细粒级土壤化探铀含量异常下限值AT(UFS)
基于步骤(1.3).后的全部数据,利用下式计算细粒级土壤化探铀含量异常下限值AT(UFS):
式中,
AT(UFS)表示细粒级土壤化探铀含量异常下限值,单位为×10-6g/g;
i表示测线的数目;
j表示每条测线上的测点数目;
m为整数,取值为1,2,3,……,i;
n为整数,取值为1,2,3,……,j;
步骤(1.5).圈定细粒级土壤化探铀异常范围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):
式中,
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):
式中,
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):
式中,
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):
式中,
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)。
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 (9)
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 | 核工业北京地质研究院 | 一种三维铀成矿有利部位的圈定方法 |
CN114609683B (zh) * | 2022-03-15 | 2024-07-12 | 核工业北京地质研究院 | 一种砂岩型铀矿砂体有利区确定方法及系统 |
CN117275601B (zh) * | 2023-11-17 | 2024-02-20 | 核工业北京地质研究院 | 砂岩型铀矿异常信息的确定方法 |
CN117723579B (zh) * | 2024-02-07 | 2024-05-28 | 核工业北京地质研究院 | 通过矿物组合确定勘测区的砂岩型铀矿矿化部位的方法 |
Family Cites Families (7)
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 | 核工业北京地质研究院 | 一种花岗岩型铀矿深部矿化信息提取方法 |
-
2019
- 2019-12-17 CN CN201911299715.0A patent/CN111045110B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111045110A (zh) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111045110B (zh) | 圈定砂岩型铀矿深部三维铀成矿靶区的综合物化探方法 | |
CN103913776B (zh) | 盾构施工中球状孤石的探测方法 | |
MX2014010955A (es) | Tecnicas de correlacion para prospeccion electrosismica y sismoelectrica pasiva. | |
CN111679342B (zh) | 快速寻找盆地隐伏砂岩型铀矿方法 | |
CN109782360A (zh) | 一种低阻覆盖区深部砂体探测方法 | |
AU2013380971A1 (en) | Sensors for passive electroseismic and seismoelectric surveying | |
CN112782773B (zh) | 一种侵入岩型铀钍铌钽矿隐伏资源预测评价方法 | |
Bryman et al. | Muon geotomography—bringing new physics to orebody imaging | |
CN105093314B (zh) | 一种测定微地震震源的方法 | |
CN115327663A (zh) | 深部矿产资源勘查用空-地-井立体地球物理探测方法 | |
Wang et al. | Coal mine goaf interpretation: survey, passive electromagnetic methods and case study | |
CN102914797B (zh) | 一种获得地层各向异性系数的方法及装置 | |
CN108008456A (zh) | 一种圈定热液型铀矿深部三维重点铀成矿有利靶区的方法 | |
Essa et al. | Magnetic inversion approach for modeling data acquired across faults: various environmental cases studies | |
Zhang et al. | Thickness Identification of Tunnel Lining Structure by Time-Energy Density Analysis based on Wavelet Transform. | |
CN113850770B (zh) | 一种盾构隧道和隔震隧道的影像ai数字化检测方法及装置 | |
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) | 一种三维铀成矿有利部位的圈定方法 | |
CN107797160A (zh) | 弹性波和电磁波ct勘测数据联合分析系统及方法 | |
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) | 一种天地孔协同一体的瞬变电磁探测采空区的方法 | |
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 |