CN106324691B - 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 - Google Patents
一种用于铀矿勘探的土壤氡浓度数据调平校正方法 Download PDFInfo
- Publication number
- CN106324691B CN106324691B CN201510385512.9A CN201510385512A CN106324691B CN 106324691 B CN106324691 B CN 106324691B CN 201510385512 A CN201510385512 A CN 201510385512A CN 106324691 B CN106324691 B CN 106324691B
- Authority
- CN
- China
- Prior art keywords
- radon consistence
- radon
- fen
- areas
- original
- 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
- Measurement Of Radiation (AREA)
Abstract
本发明属于放射性物探数据校正技术领域,具体涉及一种用于铀矿勘探的土壤氡浓度数据调平校正方法。本发明包括以下步骤:记录原始测量时的测量数据和原始测点位置,在原始测点中选择部分测点作为重复测点开展重复测量;对每个分测区计算重复测量时氡浓度背景值和原始测量时氡浓度背景值;根据同一分测区重复测量时氡浓度背景值和原始测量时氡浓度背景值相等的关系、重复测点测得数据与同位置经调平校正后的原始测点测得数据之间的误差为最小的关系,计算调平校正系数;采用调平校正系数对分测区氡浓度数据进行调平校正。本发明能够对不同年度相邻测区土壤中氡浓度数据进行有效调平校正,获取整个勘查区的土壤中氡浓度分布特征,指导铀矿勘探工作。
Description
技术领域
本发明属于放射性物探数据校正技术领域,具体涉及一种用于铀矿勘探的土壤氡浓度数据调平校正方法。
背景技术
目前,在铀矿勘查领域中,普遍采用瞬时测氡或累积测氡来获取土壤中的氡浓度,进而通过圈定氡浓度异常区来指导划定铀成矿靶区。
当氡气勘查区面积较大时,需要将勘查区划分为几个分测区,并在不同年度来完成测氡任务。由于土壤中氡浓度受气候条件的影响较大,导致不同年度相邻分测区的氡浓度数据差别较大,因此需要对各分测区的氡浓度数据进行调平校正,从而有效地去除气候条件对氡浓度值的影响,获取整个勘查区的氡浓度分布特征。
发明内容
本发明需要解决的技术问题为:现有技术难以对不同年度相邻测区的土壤氡浓度数据进行有效调平校正。
本发明的技术方案如下所述:
一种用于铀矿勘探的土壤氡浓度数据调平校正方法,包括以下步骤:
步骤1.氡浓度重复测量:将不同年度在分测区进行氡浓度原始测量时设置的测点作为原始测点,之后在同一年度对所有分测区进行重复氡浓度测量,重复氡浓度测量过程中设置的测点作为重复测点,每个分测区中的重复测点均为该分测区原始测点中的一部分,即在该分测区部分原始测点同位置处设置重复测点,重复测点在该分测区均匀分布;
步骤2.计算氡浓度背景值:分别对每个分测区计算重复测量时氡浓度背景值和原始测量时氡浓度背景值;
步骤3.计算氡浓度调平校正系数:根据同一分测区重复测量时氡浓度背景值和原始测量时氡浓度背景值相等的关系、重复测点测得数据与同位置经调平校正后的原始测点测得数据之间的误差为最小的关系,计算调平校正系数;
步骤4.对分测区氡浓度数据进行调平校正:采用调平校正系数对分测区氡浓度数据进行调平校正。
作为优选方案:步骤2中,采用“逐步替代剔除法”分别对每个分测区计算重复测量时氡浓度背景值和原始测量时氡浓度背景值,具体包括以下步骤:
步骤2.1.计算重复测量时氡浓度背景值
采用下式计算重复测点氡浓度平均值、标准偏差和异常下限值:
式中,
MY(k)为第k次迭代得到的重复测点氡浓度平均值;
SY(k)为第k次迭代得到的重复测点氡浓度标准偏差;
AY(k)为第k次迭代得到的重复测点氡浓度异常下限值;
Yi(k)为第k次迭代采用的第i个重复测点的氡浓度值,其中,Yi(1)为第i个重复测点测得的氡浓度值;
n为该分测区重复测点个数;
式(1)计算中,迭代计算终止时的重复测点氡浓度平均值即为该分测区重复测量时氡浓度背景值BY;
步骤2.2.计算原始测量时氡浓度背景值
采用下式计算原始测点氡浓度平均值、标准偏差和异常下限值:
式中,
MX(k)为第k次迭代得到的原始测点氡浓度平均值;
SX(k)为第k次迭代得到的原始测点氡浓度标准偏差;
AX(k)为第k次迭代得到的原始测点氡浓度异常下限值;
Xi(k)为第k次迭代采用的与第i个重复测点同位置的原始测点的氡浓度值,其中,Xi(1)为第i个重复测点同位置的原始测点测得的氡浓度值;
式(2)计算中,迭代计算终止时的原始测点氡浓度平均值即为该分测区原始测量时氡浓度背景值BX。
作为优选方案:步骤3具体包括以下步骤:
根据原始测量时氡浓度背景值与重复测量时氡浓度背景值相等的关系建立函数关系式:
令BY=aBX+b (3)
式中,
a、b均为调平校正系数;
采用式(4)对与第i个重复测点同位置的原始测点测得的氡浓度值Xi(1)进行校正:
Xci=aXi(1)+b i=1,2,...,n (4)
式中,
Xci为校正后的与第i个重复测点同位置的原始测点测得的氡浓度值;
采用式(5)计算Xci与第i个重复测点测得的氡浓度值Yi(1)的误差E:
将式(3)、式(4)代入式(5)得到式(6):
根据重复测点测得数据与经调平校正后同位置的原始测点测得数据之间的误差为最小的关系,计算校正系数:
根据式(7)计算调平校正系数a,将a值代入式(3)计算调平校正系数b。
作为优选方案:步骤4具体包括以下步骤:根据式(4)计算出为校正后的与第i个重复测点同位置的原始测点测得的氡浓度值即完成了对该分测区氡浓度数据的调平校正。
作为优选方案:步骤1中,对于每个分测区,该分测区重复测点数量均大于等于该分测区原始测点数量的25%。
作为优选方案:步骤1中,对所有分测区的重复氡浓度测量均采用同一种测量方法,例如:均采用瞬时测氡法,或均采用累积测氡法。
本发明的有益效果为:
本发明的一种用于铀矿勘探的土壤氡浓度数据调平校正方法,能够对不同年度相邻测区土壤中氡浓度数据进行有效调平校正,获取整个勘查区的土壤中氡浓度分布特征,指导铀矿勘探工作。
具体实施方式
下面结合实施例对本发明的一种用于铀矿勘探的土壤氡浓度数据调平校正方法进行详细说明。
一种用于铀矿勘探的土壤氡浓度数据调平校正方法,包括以下步骤:
步骤1.氡浓度重复测量:将现有技术中不同年度在分测区进行氡浓度原始测量时设置的测点作为原始测点,之后在同一年度对所有分测区采用现有技术中的测量方法进行重复氡浓度测量,重复氡浓度测量过程中设置的测点作为重复测点,每个分测区中的重复测点均为该分测区原始测点中的一部分,即在该分测区部分原始测点同位置处设置重复测点,重复测点在该分测区均匀分布;
步骤2.计算氡浓度背景值:采用“逐步替代剔除法”分别对每个分测区计算重复测量时氡浓度背景值和原始测量时氡浓度背景值;
步骤3.计算氡浓度调平校正系数:根据相同分测区重复测量时氡浓度背景值和原始测量时氡浓度背景值相等的关系、重复测点测得数据与经调平校正后的原始测点测得数据之间的误差为最小的关系,计算调平校正系数;
步骤4.对分测区氡浓度数据进行调平校正:采用校正系数对分测区氡浓度数据进行调平校正。
实施例1
步骤1.氡浓度重复测量
在同一年度对所有分测区采用现有技术中的测量方法进行重复氡浓度测量:对所有分测区的重复氡浓度测量均采用同一种测量方法,如均采用瞬时测氡法或均采用累积测氡法;对于每个分测区,重复测点数量均需大于等于该分测区原始测点数量的25%。
步骤2.计算氡浓度背景值
针对每一个分测区,分别计算该分测区内所有重复测点和其同位置的原始测点的背景值。
步骤2.1.计算重复测量时氡浓度背景值
采用下式计算重复测点氡浓度平均值、标准偏差和异常下限值:
式中,
MY(k)为第k次迭代得到的重复测点氡浓度平均值;
SY(k)为第k次迭代得到的重复测点氡浓度标准偏差;
AY(k)为第k次迭代得到的重复测点氡浓度异常下限值;
Yi(k)为第k次迭代采用的第i个重复测点的氡浓度值,其中,Yi(1)为第i个重复测点测得的氡浓度值;
n为该分测区重复测点个数。
式(1)计算中,当k=p时,不再存在大于AY(p)的Yi(p)数据,迭代计算终止。迭代计算终止时的氡浓度平均值MY(p)即为该分测区重复测量时氡浓度背景值BY。
步骤2.2.计算原始测量时氡浓度背景值
采用下式计算原始测点氡浓度平均值、标准偏差和异常下限值:
式中,
MX(k)为第k次迭代得到的原始测点氡浓度平均值;
SX(k)为第k次迭代得到的原始测点氡浓度标准偏差;
AX(k)为第k次迭代得到的原始测点氡浓度异常下限值;
Xi(k)为第k次迭代采用的与第i个重复测点同位置的原始测点的氡浓度值,其中,Xi(1)为第i个重复测点同位置的原始测点测得的氡浓度值.
式(2)计算中,当k=q时,不再存在大于AX(q)的Xi(q)数据,迭代计算终止。迭代计算终止时的氡浓度平均值MX(q)即为该分测区原始测量时氡浓度背景值BX。
步骤3.计算氡浓度调平校正系数
根据原始测量时氡浓度背景值与重复测量时氡浓度背景值相等的关系建立函数关系式:
令BY=aBX+b (3)
式中,
a、b均为调平校正系数。
采用式(4)对与第i个重复测点同位置的原始测点测得的氡浓度值Xi(1)进行校正:
Xci=aXi(1)+b i=1,2,...,n (4)
式中,
Xci为校正后的与第i个重复测点同位置的原始测点测得的氡浓度值。
采用式(5)计算Xci与第i个重复测点测得的氡浓度值Yi(1)的误差E:
将式(3)、式(4)代入式(5)得到式(6):
根据重复测点测得数据与经调平校正后同位置的原始测点测得数据之间的误差为最小的关系,计算校正系数:
根据式(7)计算调平校正系数a,将a值代入式(3)计算调平校正系数b。
步骤4.对分测区氡浓度数据进行调平校正
根据式(4)计算出校正后的与第i个重复测点同位置的原始测点测得的氡浓度值即完成了对该分测区氡浓度数据的调平校正。
其他分测区采用与上述方法逐一进行调平校正。
Claims (7)
1.一种用于铀矿勘探的土壤氡浓度数据调平校正方法,其特征在于:包括以下步骤:
步骤1.氡浓度重复测量:将不同年度在分测区进行氡浓度原始测量时设置的测点作为原始测点,之后在同一年度对所有分测区进行重复氡浓度测量,重复氡浓度测量过程中设置的测点作为重复测点,每个分测区中的重复测点均为该分测区原始测点中的一部分,即在该分测区部分原始测点同位置处设置重复测点,重复测点在该分测区均匀分布;
步骤2.计算氡浓度背景值:分别对每个分测区计算重复测量时氡浓度背景值和原始测量时氡浓度背景值;
步骤3.计算氡浓度调平校正系数:根据同一分测区重复测量时氡浓度背景值和原始测量时氡浓度背景值相等的关系、重复测点测得数据与同位置经调平校正后的原始测点测得数据之间的误差为最小的关系,计算调平校正系数;
步骤4.对分测区氡浓度数据进行调平校正:采用调平校正系数对分测区氡浓度数据进行调平校正。
2.根据权利要求1所述的一种用于铀矿勘探的土壤氡浓度数据调平校正方法,其特征在于:步骤2中,采用“逐步替代剔除法”分别对每个分测区计算重复测量时氡浓度背景值和原始测量时氡浓度背景值,具体包括以下步骤:
步骤2.1.计算重复测量时氡浓度背景值
采用下式计算重复测点氡浓度平均值、氡浓度标准偏差和氡浓度异常下限值:
式中,
MY(k)为第k次迭代得到的重复测点氡浓度平均值;
SY(k)为第k次迭代得到的重复测点氡浓度标准偏差;
AY(k)为第k次迭代得到的重复测点氡浓度异常下限值;
Yi(k)为第k次迭代采用的第i个重复测点的氡浓度值,其中,Yi(1)为第i个重复测点测得的氡浓度值;
n为该分测区重复测点个数;
式(1)计算中,迭代计算终止时的重复测点氡浓度平均值即为该分测区重复测量时氡浓度背景值BY;
步骤2.2.计算原始测量时氡浓度背景值
采用下式计算原始测点氡浓度平均值、氡浓度标准偏差和氡浓度异常下限值:
式中,
MX(k)为第k次迭代得到的原始测点氡浓度平均值;
SX(k)为第k次迭代得到的原始测点氡浓度标准偏差;
AX(k)为第k次迭代得到的原始测点氡浓度异常下限值;
Xi(k)为第k次迭代采用的与第i个重复测点同位置的原始测点的氡浓度值,其中,Xi(1)为第i个重复测点同位置的原始测点测得的氡浓度值;
式(2)计算中,迭代计算终止时的原始测点氡浓度平均值即为该分测区原始测量时氡浓度背景值BX。
3.根据权利要求2所述的一种用于铀矿勘探的土壤氡浓度数据调平校正方法,其特征在于:步骤3具体包括以下步骤:
根据原始测量时氡浓度背景值与重复测量时氡浓度背景值相等的关系建立函数关系式:
令BY=aBX+b (3)
式中,
a、b均为调平校正系数;
采用式(4)对与第i个重复测点同位置的原始测点测得的氡浓度值Xi(1)进行校正:
式中,
Xci为校正后的与第i个重复测点同位置的原始测点测得的氡浓度值;
采用式(5)计算Xci与第i个重复测点测得的氡浓度值Yi(1)的误差E:
将式(3)、式(4)代入式(5)得到式(6):
根据重复测点测得数据与经调平校正后同位置的原始测点测得数据之间的误差为最小的关系,计算调平校正系数:
令
根据式(7)计算调平校正系数a,将a值代入式(3)计算调平校正系数b。
4.根据权利要求3所述的一种用于铀矿勘探的土壤氡浓度数据调平校正方法,其特征在于:步骤4具体包括以下步骤:根据式(4)计算出校正后的与第i个重复测点同位置的原始测点测得的氡浓度值即完成了对该分测区氡浓度数据的调平校正。
5.根据权利要求1或2或3或4所述的一种用于铀矿勘探的土壤氡浓度数据调平校正方法,其特征在于:步骤1中,对于每个分测区,该分测区重复测点数量均大于等于该分测区原始测点数量的25%。
6.根据权利要求1或2或3或4所述的一种用于铀矿勘探的土壤氡浓度数据调平校正方法,其特征在于:步骤1中,对所有分测区的重复氡浓度测量均采用同一种测量方法。
7.根据权利要求6所述的一种用于铀矿勘探的土壤氡浓度数据调平校正方法,其特征在于:步骤1中,对所有分测区的重复氡浓度测量均采用瞬时测氡法或均采用累积测氡法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510385512.9A CN106324691B (zh) | 2015-06-30 | 2015-06-30 | 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510385512.9A CN106324691B (zh) | 2015-06-30 | 2015-06-30 | 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106324691A CN106324691A (zh) | 2017-01-11 |
CN106324691B true CN106324691B (zh) | 2018-06-12 |
Family
ID=57728492
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510385512.9A Active CN106324691B (zh) | 2015-06-30 | 2015-06-30 | 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106324691B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107632137B (zh) * | 2017-08-25 | 2019-12-20 | 核工业北京地质研究院 | 一种铀矿勘查中土壤氡气瞬时测量数据校正方法 |
CN111815896A (zh) * | 2020-06-19 | 2020-10-23 | 淮阴师范学院 | 一种基于物联网的报警通信管理系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5833395A (en) * | 1997-04-24 | 1998-11-10 | Atlantic Richfield Company | Method for reducing the concentration of contaminating radioactive material in contaminated soil |
CN103091472A (zh) * | 2013-01-25 | 2013-05-08 | 南华大学 | 一种快速测量土壤潜势氡浓度的方法 |
CN103529481A (zh) * | 2013-08-13 | 2014-01-22 | 核工业北京地质研究院 | 一种深部铀矿化信息探测方法 |
CN103576208A (zh) * | 2013-11-07 | 2014-02-12 | 核工业北京地质研究院 | 一种面向铀矿床定位的瞬时测氡数据异常提取方法 |
CN103837908A (zh) * | 2014-03-05 | 2014-06-04 | 核工业北京地质研究院 | 一种适用于隐伏砂岩型铀矿快速找矿定位方法 |
CN103913781A (zh) * | 2014-03-18 | 2014-07-09 | 核工业北京地质研究院 | 一种预测铀矿体埋深的组合探测方法 |
-
2015
- 2015-06-30 CN CN201510385512.9A patent/CN106324691B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5833395A (en) * | 1997-04-24 | 1998-11-10 | Atlantic Richfield Company | Method for reducing the concentration of contaminating radioactive material in contaminated soil |
CN103091472A (zh) * | 2013-01-25 | 2013-05-08 | 南华大学 | 一种快速测量土壤潜势氡浓度的方法 |
CN103529481A (zh) * | 2013-08-13 | 2014-01-22 | 核工业北京地质研究院 | 一种深部铀矿化信息探测方法 |
CN103576208A (zh) * | 2013-11-07 | 2014-02-12 | 核工业北京地质研究院 | 一种面向铀矿床定位的瞬时测氡数据异常提取方法 |
CN103837908A (zh) * | 2014-03-05 | 2014-06-04 | 核工业北京地质研究院 | 一种适用于隐伏砂岩型铀矿快速找矿定位方法 |
CN103913781A (zh) * | 2014-03-18 | 2014-07-09 | 核工业北京地质研究院 | 一种预测铀矿体埋深的组合探测方法 |
Non-Patent Citations (2)
Title |
---|
《铀矿井下氡及子体浓度10年变化规律分析》;苏世标等;《第四届全国环境与职业医学研究生学术研讨会论文集》;20051231;264 * |
《铀矿采场氡浓度模拟初步探讨》;冀东等;《环境科学与技术》;20140630;第37卷(第6N期);304-307 * |
Also Published As
Publication number | Publication date |
---|---|
CN106324691A (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104133245B (zh) | 一种地震资料的静校正方法及系统 | |
CN105301639B (zh) | 基于vsp旅行时双加权层析反演速度场的方法及其装置 | |
WO2016141598A1 (zh) | 共姿态道集地表一致性振幅补偿方法及装置 | |
CN103529481B (zh) | 一种深部铀矿化信息探测方法 | |
CN107870001B (zh) | 一种基于椭球拟合的磁力计校正方法 | |
WO2009117784A1 (en) | System, apparatus and method for measuring soil moisture content | |
CN106324691B (zh) | 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 | |
CN110333543B (zh) | 基于反射系数分析的低阻体解释及成像方法与系统 | |
CN104776827B (zh) | Gps高程异常数据的粗差探测方法 | |
CN104360396B (zh) | 一种海上井间tti介质三种初至波走时层析成像方法 | |
CN104406715A (zh) | 一种遥感估算地表感热/潜热通量的精度评价方法及系统 | |
CN104363649A (zh) | 带有约束条件的ukf的wsn节点定位方法 | |
RU2009101300A (ru) | Коррекция влияния отклонения прибора при каротаже плотности в процессе бурения | |
CN106324692B (zh) | 一种地面伽玛能谱测量铀含量数据调平校正方法 | |
CN106802164A (zh) | 一种适用于粮情监测的绝对水势计算方法及云图生成方法 | |
CN108008456B (zh) | 一种圈定热液型铀矿深部三维重点铀成矿有利靶区的方法 | |
CN111721831A (zh) | 基于电刺激的三维层析成像堤坝隐伏渗漏通道扫描方法 | |
US10579756B2 (en) | Simulation method of surface water flow movement process in surface irrigation | |
CN106569259A (zh) | 基于构造倾角的正则化层析速度反演方法和装置 | |
CN105717538B (zh) | 起伏地表地震数据偏移基准面转换方法及装置 | |
CN113449254B (zh) | 任意网型变形监测稳定性分析方法及监控点位置确定方法 | |
CN105093318A (zh) | 一种自适应波动方程波场延拓静校正方法 | |
CN106526078A (zh) | 一种用于提取深部铀成矿信息的氡气‑地气联合测量方法 | |
CN117092702A (zh) | 孔-隧激发极化探水结构的施工方法及反演探水方法 | |
CN108387194B (zh) | 一种井下闭合导线坐标测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |