CN103576208A - 一种面向铀矿床定位的瞬时测氡数据异常提取方法 - Google Patents
一种面向铀矿床定位的瞬时测氡数据异常提取方法 Download PDFInfo
- Publication number
- CN103576208A CN103576208A CN201310551918.0A CN201310551918A CN103576208A CN 103576208 A CN103576208 A CN 103576208A CN 201310551918 A CN201310551918 A CN 201310551918A CN 103576208 A CN103576208 A CN 103576208A
- Authority
- CN
- China
- Prior art keywords
- point
- value
- data
- radon
- result
- 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.)
- Granted
Links
Landscapes
- Measurement Of Radiation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种面向铀矿床定位的瞬时测氡数据异常提取方法,依次包括:1、每个瞬时测氡布点采集一个测量值,所有测量值组成二位离散数据A;2、面向二位离散数据A建立一个均值滤波器B;3、对测量数据A采用均值滤波器B进行滤波处理,结果记录为C;4、对C每一点的数值逐点求平方,结果记录为D;5、对测量数据A的测量数值逐点求平方,并采均值滤波器B进行滤波处理,结果记录为E;6、对E和D逐点求差,并对差值逐点开平方,结果记录为F;7、用F数值逐点除以对应的C点数值,获得新的异常提取结果G。本发明利用这种不确定性来提取氡的异常分布,采取本发明具体的步骤及参数,可在砂岩型铀矿床定位中获得更高的定位精度。
Description
技术领域
本发明涉及一种面向铀矿床定位的瞬时测氡数据异常提取方法,特别是涉及一种通过分析瞬时测氡数据中由于覆盖层厚度对实测数据带来的不确定影响程度的面向铀矿床定位的瞬时测氡数据异常提取方法。
背景技术
氡是铀的直接衰变物,也是自然界唯一的放射性的气体,由镭衰变产生,镭又是由放射性元素铀衰变而来,铀起了一个氡的永久源作用,因此,在铀矿地质中,土壤和地表瞬时测氡是一种重要的铀矿床攻深找盲方法。大量矿床实例证明,深达几百米的铀矿石可以在地表形成明显的氡异常。
在数据处理和异常提取过程中,受到地质覆盖层岩性及其厚度和测量仪器的误差影响,一般对瞬时测氡数据只进行滤波或者小波变换等噪声滤除处理,数值较大的测量点即被视为异常点,毗邻的多个异常点形成异常区,通过绘制等值线进行分析应用。
不同地区的铀矿埋深是不同的,地质覆盖层的岩性、产状也不尽相同,地貌景观也有所差别,现有数据处理方法对此考虑不全面,只是简单的通过数值大小来判断异常区的位置。
因此亟需提供一种新型的面向铀矿床定位的瞬时测氡数据异常提取方法。
发明内容
本发明要解决的技术问题是提供一种在砂岩型铀矿床定位中获得更高的定位精度的面向铀矿床定位的瞬时测氡数据异常提取方法。
为解决上述技术问题,本发明一种面向铀矿床定位的瞬时测氡数据异常提取方法,依次包括以下步骤:
步骤1、确定待分析的瞬时测氡数据;
每个瞬时测氡布点采集一个测量值,所有测量值组成二位离散数据A;
步骤2、设定均值滤波器;
面向二位离散数据A建立一个均值滤波器B,所述均值滤波器的窗口尺寸为奇数;
步骤3、对测量数据A采用均值滤波器B进行滤波处理,结果记录为C;
步骤4、对C每一点的数值逐点求平方,结果记录为D;
步骤5、对测量数据A的测量数值逐点求平方,并采均值滤波器B进行滤波处理,结果记录为E;
步骤6、对E和D逐点求差,并对差值逐点开平方,结果记录为F;
步骤7、用F数值逐点除以对应的C点数值,获得新的异常提取结果G。
均值滤波器的窗口尺寸为5。
受到地表覆盖层的影响,氡和氡的母体和子体元素在由地下向地表扩散的过程中,其扩散路径受到周围地质体、岩性、产状、厚度、含水性和温度等因素的影响会不断发生变化。与当前仅直接通过测量数值大小确定氡异常的方法相比,本发明考虑到这一影响因素,并利用这种不确定性来提取氡的异常分布,采取本发明具体的步骤及参数,可在砂岩型铀矿床定位中获得更高的定位精度。
具体实施方式
步骤1、确定待分析的瞬时测氡数据;
通常应用于地质矿床定位的瞬时测氡布点构成一个二维矩阵,布点的位置为了便于分析通常沿行和列方向等距,非等距情况下可以通过数据重采样获得,每个布点记录一个测量值,在实际处理中也可以被认为是一幅离散图像并记录为A;
步骤2、设定均值滤波器;
面向二位离散数据A建立一个均值滤波器B,所述均值滤波器的窗口尺寸为奇数;窗口尺寸的大小可根据实测数据和区域地质情况确定,当地质覆盖层比较均匀时,滤波器的窗口尺寸可以适当取大一些;反之,则滤波器的窗口尺寸可以适当取小一些;优选取值5;
步骤3、对测量数据A采用均值滤波器B进行滤波处理,结果记录为C;
步骤4、对C每一点的数值逐点求平方,结果记录为D;
步骤5、对测量数据A的测量数值逐点求平方,并采均值滤波器B进行滤波处理,结果记录为E;
步骤6、对E和D逐点求差,并对差值逐点开平方,结果记录为F;
步骤7、用F数值逐点除以对应的C点数值,获得新的异常提取结果G。
Claims (2)
1.一种面向铀矿床定位的瞬时测氡数据异常提取方法,依次包括以下步骤:
步骤1、确定待分析的瞬时测氡数据;
每个瞬时测氡布点采集一个测量值,所有测量值组成二位离散数据A;
步骤2、设定均值滤波器;
面向二位离散数据A建立一个均值滤波器B,所述均值滤波器的窗口尺寸为奇数;
步骤3、对测量数据A采用均值滤波器B进行滤波处理,结果记录为C;
步骤4、对C每一点的数值逐点求平方,结果记录为D;
步骤5、对测量数据A的测量数值逐点求平方,并采均值滤波器B进行滤波处理,结果记录为E;
步骤6、对E和D逐点求差,并对差值逐点开平方,结果记录为F;
步骤7、用F数值逐点除以对应的C点数值,获得新的异常提取结果G。
2.根据权利要求1所述的一种面向铀矿床定位的瞬时测氡数据异常提取方法,其特征在于:所述均值滤波器的窗口尺寸为5。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310551918.0A CN103576208B (zh) | 2013-11-07 | 2013-11-07 | 一种面向铀矿床定位的瞬时测氡数据异常提取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310551918.0A CN103576208B (zh) | 2013-11-07 | 2013-11-07 | 一种面向铀矿床定位的瞬时测氡数据异常提取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103576208A true CN103576208A (zh) | 2014-02-12 |
| CN103576208B CN103576208B (zh) | 2016-05-18 |
Family
ID=50048370
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310551918.0A Active CN103576208B (zh) | 2013-11-07 | 2013-11-07 | 一种面向铀矿床定位的瞬时测氡数据异常提取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103576208B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103837908A (zh) * | 2014-03-05 | 2014-06-04 | 核工业北京地质研究院 | 一种适用于隐伏砂岩型铀矿快速找矿定位方法 |
| CN106324691A (zh) * | 2015-06-30 | 2017-01-11 | 核工业北京地质研究院 | 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 |
| CN108241176A (zh) * | 2016-12-26 | 2018-07-03 | 核工业北京地质研究院 | 一种放射性综合信息异常远景区圈定方法 |
| CN111045114A (zh) * | 2019-12-31 | 2020-04-21 | 核工业北京地质研究院 | 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法 |
| CN114609683A (zh) * | 2022-03-15 | 2022-06-10 | 核工业北京地质研究院 | 一种砂岩型铀矿砂体有利区确定方法及系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4468558A (en) * | 1980-10-20 | 1984-08-28 | Boliden Aktiebolag | Method and apparatus for measuring radon beneath and above the level of the ground |
| US4920263A (en) * | 1988-01-26 | 1990-04-24 | Gemini Research, Inc. | Radon detection system |
-
2013
- 2013-11-07 CN CN201310551918.0A patent/CN103576208B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4468558A (en) * | 1980-10-20 | 1984-08-28 | Boliden Aktiebolag | Method and apparatus for measuring radon beneath and above the level of the ground |
| US4920263A (en) * | 1988-01-26 | 1990-04-24 | Gemini Research, Inc. | Radon detection system |
Non-Patent Citations (2)
| Title |
|---|
| 刘敦旺: "基于小波分析的活性炭测氡数据处理与应用", 《中国优秀硕士学位论文全文数据库(工程科技I辑)》 * |
| 杨镜明 等: "小波变换在砂岩型铀矿地面氡气测量中的应用", 《铀矿地质》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103837908A (zh) * | 2014-03-05 | 2014-06-04 | 核工业北京地质研究院 | 一种适用于隐伏砂岩型铀矿快速找矿定位方法 |
| CN106324691A (zh) * | 2015-06-30 | 2017-01-11 | 核工业北京地质研究院 | 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 |
| CN106324691B (zh) * | 2015-06-30 | 2018-06-12 | 核工业北京地质研究院 | 一种用于铀矿勘探的土壤氡浓度数据调平校正方法 |
| CN108241176A (zh) * | 2016-12-26 | 2018-07-03 | 核工业北京地质研究院 | 一种放射性综合信息异常远景区圈定方法 |
| CN108241176B (zh) * | 2016-12-26 | 2019-08-13 | 核工业北京地质研究院 | 一种放射性综合信息异常远景区圈定方法 |
| CN111045114A (zh) * | 2019-12-31 | 2020-04-21 | 核工业北京地质研究院 | 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法 |
| CN111045114B (zh) * | 2019-12-31 | 2021-11-30 | 核工业北京地质研究院 | 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法 |
| CN114609683A (zh) * | 2022-03-15 | 2022-06-10 | 核工业北京地质研究院 | 一种砂岩型铀矿砂体有利区确定方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103576208B (zh) | 2016-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103576208B (zh) | 一种面向铀矿床定位的瞬时测氡数据异常提取方法 | |
| CN103809217B (zh) | 一种计算高、低伽马交互储层泥质含量的方法 | |
| CN110886604B (zh) | 一种基于计算机模拟技术的高效地热资源勘察方法 | |
| US20120223235A1 (en) | Systems and Methods for Determining Fluid Mobility in Rock Samples | |
| CN108019204B (zh) | 基于物性和地层水矿化度校正的中子伽马气层识别方法 | |
| CN106846475B (zh) | 一种河流历史形态数字复原方法 | |
| Olierook et al. | Paleodrainage and fault development in the southern Perth Basin, Western Australia during and after the breakup of Gondwana from 3D modelling of the Bunbury Basalt | |
| Kuo et al. | Anomalous decrease in groundwater radon before 2016 Mw 6.4 Meinong earthquake and its application in Taiwan | |
| Kulongoski et al. | Applications of groundwater helium | |
| CN106650192A (zh) | 一种火山岩型铀矿床磁性界面反演方法 | |
| CN103678778A (zh) | 一种放射性物化探信息综合方法 | |
| CN107272080B (zh) | 一种地下水渗漏通道勘察的场标准化分析方法 | |
| Schillig et al. | Upscaling point velocity measurements to characterize a glacial outwash aquifer | |
| CN108008456B (zh) | 一种圈定热液型铀矿深部三维重点铀成矿有利靶区的方法 | |
| CN114764150B (zh) | 一种断层活动强度和活动速率的计算方法及系统 | |
| CN107478392B (zh) | 基于氡同位素的沿海进排水闸设施渗漏检测方法 | |
| Nielsen et al. | New age for ferromanganese crust 109D‐C and implications for isotopic records of lead, neodymium, hafnium, and thallium in the Pliocene Indian Ocean | |
| CN117890997B (zh) | 一种基于土壤氢的快速识别硫化物多金属矿的方法 | |
| CN112666096A (zh) | 一种砂岩型铀矿土壤找矿信息高光谱提取方法 | |
| CN113836870B (zh) | 陆棚沉积亚相确定方法及装置 | |
| Zschornack et al. | Hydraulic profiling with the direct-push permeameter: assessment of probe configuration and analysis methodology | |
| CN108756869B (zh) | 碳酸盐岩低级序不整合面的识别方法 | |
| CN119308671B (zh) | 一种页岩气保存条件评价参数获取方法、装置和电子设备 | |
| Wagner | New developments in electrical resistivity imaging with applications to geological CO₂ storage | |
| Morales et al. | Electrical resistivity monitoring of CO2 injection at the Mont Terri underground laboratory, Switzerland |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |