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
- data
- pointwise
- mean filter
- radon
- 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 (4)
| 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 | 核工业北京地质研究院 | 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法 |
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 (7)
| 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 | 核工业北京地质研究院 | 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103576208B (zh) | 2016-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Breitenbach et al. | Constructing proxy records from age models (COPRA) | |
| US7266041B1 (en) | Multi-attribute background relative scanning of 3D geophysical datasets for locally anomaluous data points | |
| US10677046B2 (en) | Leakage detection using smart field technology | |
| CN103913781B (zh) | 一种预测铀矿体埋深的组合探测方法 | |
| WO2014138191A1 (en) | Geological model analysis incorporating cross-well electromagnetic measurements | |
| CN103576208A (zh) | 一种面向铀矿床定位的瞬时测氡数据异常提取方法 | |
| CN103678778B (zh) | 一种放射性物化探信息综合方法 | |
| CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
| CN105116452A (zh) | 一种确定地质异常体电阻率和极化率的方法和装置 | |
| CN105114067A (zh) | 一种岩性测井相的方法 | |
| CN105093297A (zh) | 一种微地震定位精度的质量控制方法 | |
| CN105301638B (zh) | 一种提取风化层底界面的方法和装置 | |
| CN103852789A (zh) | 用于地震数据的非线性层析方法及其装置 | |
| CN108008456A (zh) | 一种圈定热液型铀矿深部三维重点铀成矿有利靶区的方法 | |
| Selzer et al. | RaDeCC Reader: Fast, accurate and automated data processing for Radium Delayed Coincidence Counting systems | |
| Pang et al. | 3D multi-scale resistivity inversion method applied in the tunnel face to borehole observations for tunnel-ahead prospecting | |
| KR100978838B1 (ko) | 라돈 농도 측정을 통한 지하 단층 조사 방법 | |
| CN116045903A (zh) | 一种采煤区地面形变识别与评价方法 | |
| CN113031076B (zh) | 一种在电阻率法超前探测中考虑围岩各向异性影响的方法 | |
| CN104809745A (zh) | 一种改进低空间分辨率遥感影像提取的湖面面积的方法 | |
| Robinson et al. | Portable gamma ray spectrometry for archaeological prospection: A preliminary investigation at Silchester Roman Town | |
| Zielke et al. | Bedrock fault roughness resolves slip increments of large earthquakes: Case studies from Central Italy | |
| CN110927789B (zh) | 一种基于欠失数据预测页岩平面分布的方法及装置 | |
| Desherevskii et al. | Estimating the number of blasts in the Dushanbe-Vakhsh earthquake catalog | |
| CN110764163A (zh) | 基于γ总量及高密度电法测量下的伟晶岩锂矿的圈定方法 |
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 |