CN106324692B - 一种地面伽玛能谱测量铀含量数据调平校正方法 - Google Patents
一种地面伽玛能谱测量铀含量数据调平校正方法 Download PDFInfo
- Publication number
- CN106324692B CN106324692B CN201510385514.8A CN201510385514A CN106324692B CN 106324692 B CN106324692 B CN 106324692B CN 201510385514 A CN201510385514 A CN 201510385514A CN 106324692 B CN106324692 B CN 106324692B
- Authority
- CN
- China
- Prior art keywords
- uranium content
- fen
- areas
- original
- uranium
- 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
- Other Investigation Or Analysis Of Materials By Electrical Means (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中，对于每个分测区，该分测区重复测点数量均大于等于该分测区原始测点数量的20％。
作为优选方案：步骤1中，对所有分测区的重复铀含量测量均采用地面伽玛能谱测量方法。
本发明的有益效果为：
本发明的一种地面伽玛能谱测量铀含量数据调平校正方法，能够对不同年度相邻测区地面伽玛能谱测量铀含量数据进行有效调平校正，获取整个勘查区的铀含量分布特征，指导铀矿勘探工作。
具体实施方式
下面结合实施例对本发明的一种地面伽玛能谱测量铀含量数据调平校正方法进行详细说明。
一种地面伽玛能谱测量铀含量数据调平校正方法，包括以下步骤：
步骤1.铀含量重复测量：将现有技术中不同年度在分测区进行铀含量原始测量时设置的测点作为原始测点，之后在同一年度对所有分测区采用现有技术中的测量方法进行重复铀含量测量，重复铀含量测量过程中设置的测点作为重复测点，每个分测区中的重复测点均为该分测区原始测点中的一部分，即在该分测区部分原始测点同位置处设置重复测点，重复测点在该分测区均匀分布；
步骤2.计算铀含量背景值：采用“逐步替代剔除法”分别对每个分测区计算重复测量时铀含量背景值和原始测量时铀含量背景值；
步骤3.计算铀含量调平校正系数：根据相同分测区重复测量时铀含量背景值和原始测量时铀含量背景值相等的关系、重复测点测得数据与经调平校正后的原始测点测得数据之间的误差为最小的关系，计算调平校正系数；
步骤4.对分测区铀含量数据进行调平校正：采用校正系数对分测区铀含量数据进行调平校正。
实施例1
步骤1.铀含量重复测量
在同一年度对所有分测区采用现有技术中的地面伽玛能谱测量方法进行重复铀含量测量，对于每个分测区，重复测点数量均需大于等于该分测区原始测点数量的20％。
步骤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 (6)
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＝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.根据权利要求3所述的一种地面伽玛能谱测量铀含量数据调平校正方法，其特征在于：步骤4具体包括以下步骤：根据式(4)计算出校正后的与第i个重复测点同位置的原始测点测得的铀含量值即完成了对该分测区铀含量数据的调平校正。
5.根据权利要求1或2或3或4所述的一种地面伽玛能谱测量铀含量数据调平校正方法，其特征在于：步骤1中，对于每个分测区，该分测区重复测点数量均大于等于该分测区原始测点数量的20％。
6.根据权利要求1或2或3或4所述的一种地面伽玛能谱测量铀含量数据调平校正方法，其特征在于：步骤1中，对所有分测区的重复铀含量测量均采用地面伽玛能谱测量方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510385514.8A CN106324692B (zh) | 2015-06-30 | 2015-06-30 | 一种地面伽玛能谱测量铀含量数据调平校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510385514.8A CN106324692B (zh) | 2015-06-30 | 2015-06-30 | 一种地面伽玛能谱测量铀含量数据调平校正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106324692A CN106324692A (zh) | 2017-01-11 |
| CN106324692B true CN106324692B (zh) | 2018-06-26 |
Family
ID=57728421
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510385514.8A Active CN106324692B (zh) | 2015-06-30 | 2015-06-30 | 一种地面伽玛能谱测量铀含量数据调平校正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106324692B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107728221A (zh) * | 2017-09-19 | 2018-02-23 | 中国地质调查局沈阳地质调查中心 | 一种放射性异常段的探测方法和系统 |
| CN111045097A (zh) * | 2019-12-30 | 2020-04-21 | 核工业北京地质研究院 | 一种钠交代型铀矿化的初步圈定方法 |
Citations (4)
| 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 |
| CN103913781A (zh) * | 2014-03-18 | 2014-07-09 | 核工业北京地质研究院 | 一种预测铀矿体埋深的组合探测方法 |
| CN103913780A (zh) * | 2014-03-18 | 2014-07-09 | 核工业北京地质研究院 | 一种用于勘查钙结岩型铀矿的方法 |
| CN104062307A (zh) * | 2014-07-01 | 2014-09-24 | 核工业北京地质研究院 | 一种车载伽玛能谱大气氡影响修正方法 |
-
2015
- 2015-06-30 CN CN201510385514.8A patent/CN106324692B/zh active Active
Patent Citations (4)
| 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 |
| CN103913781A (zh) * | 2014-03-18 | 2014-07-09 | 核工业北京地质研究院 | 一种预测铀矿体埋深的组合探测方法 |
| CN103913780A (zh) * | 2014-03-18 | 2014-07-09 | 核工业北京地质研究院 | 一种用于勘查钙结岩型铀矿的方法 |
| CN104062307A (zh) * | 2014-07-01 | 2014-09-24 | 核工业北京地质研究院 | 一种车载伽玛能谱大气氡影响修正方法 |
Non-Patent Citations (1)
| Title |
|---|
| 《铀矿采场氡浓度模拟初步探讨》;冀东等;《环境科学与技术》;20140630;第37卷(第6N期);304-307 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106324692A (zh) | 2017-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104133245B (zh) | 一种地震资料的静校正方法及系统 | |
| CN107239659B (zh) | 一种基于改进的k-rusle模型与成土速率计算土壤侵蚀的方法 | |
| CN111045072B (zh) | 一种适用于CeBr3探测器的γ能谱迭代解谱方法 | |
| CN106324692B (zh) | 一种地面伽玛能谱测量铀含量数据调平校正方法 | |
| CN102901985B (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
| CA2381107A1 (en) | Standoff compensation for nuclear measurements | |
| CN104090301B (zh) | 一种求取三维高频静校正量的方法 | |
| CN107132578A (zh) | 一种微地震地面监测速度模型校正算法 | |
| US20230384473A1 (en) | Method and system for advanced detection and optimization of tunnel resistivity based on depth resolution | |
| CN104360396B (zh) | 一种海上井间tti介质三种初至波走时层析成像方法 | |
| CN110333543A (zh) | 基于反射系数分析的低阻体解释及成像方法与系统 | |
| CN104406715A (zh) | 一种遥感估算地表感热/潜热通量的精度评价方法及系统 | |
| CN106324691B (zh) | 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 | |
| CN103678778A (zh) | 一种放射性物化探信息综合方法 | |
| CN108008456B (zh) | 一种圈定热液型铀矿深部三维重点铀成矿有利靶区的方法 | |
| CN107657119A (zh) | 提高伽马能谱数据品质的测量方法 | |
| CN116381772A (zh) | 一种用于剂量率仪的实时能量响应校正方法、系统及终端 | |
| CN105717538B (zh) | 起伏地表地震数据偏移基准面转换方法及装置 | |
| CN105093318A (zh) | 一种自适应波动方程波场延拓静校正方法 | |
| CN102819042A (zh) | 一种抽稀地震数据的近地表层析反演静校正方法 | |
| CN108387194B (zh) | 一种井下闭合导线坐标测量方法 | |
| CN106248740B (zh) | 一种土壤电阻率湿度校正方法 | |
| CN114236624B (zh) | 基于电磁法估算压裂改造空间体积的方法和系统 | |
| CN104504273A (zh) | 一种提高γ射线峰探测灵敏度的能谱处理方法 | |
| CN105445788A (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 |