CN104880726A - 测量铀矿回风井氡排放量的方法 - Google Patents
测量铀矿回风井氡排放量的方法 Download PDFInfo
- Publication number
- CN104880726A CN104880726A CN201510320774.7A CN201510320774A CN104880726A CN 104880726 A CN104880726 A CN 104880726A CN 201510320774 A CN201510320774 A CN 201510320774A CN 104880726 A CN104880726 A CN 104880726A
- Authority
- CN
- China
- Prior art keywords
- returnairshaft
- radon
- uranium ore
- formula
- unit
- 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)
Abstract
一种测量铀矿回风井氡排放量的方法,首先采用多点法对铀矿回风井井口特定样点进行氡浓度测量,然后利用插值积分法计算出回风井井口径向特定区间的平均出口风速和平均氡浓度,再将各区间平均风速、氡浓度平均值及特定区间面积三者的乘积进行加权计算出回风井的氡释放率。对于某特定回风井,其释放率与某截面井中心风速、井中心氡浓度及该截面面积之间存在恒定的数值关系,利用该规律和上述测量数据计算得到该回风井的转换系数。在铀矿回风井氡释放率的后续测量中,只需测得回风井某截面井中心处氡浓度及风速并乘以转换系数和回风井出口截面积即能够计算出回风井的氡释放率,最后再乘以回风井氡的排放时间即能够计算出铀矿回风井氡的排放量。
Description
技术领域
本发明涉及铀矿冶工业技术领域,特别是一种测量铀矿回风井氡排放量的方法。
背景技术
在铀矿开采、勘测及矿石加工的过程中,铀矿石中析出的氡及其子体对铀矿工作人员的健康造成危害且污染铀矿周边环境。目前,在铀矿开采过程中,均采用机械通风将矿井内析出的氡经由回风井排放到大气中,为保障铀矿工作人员的健康并确保对铀矿附近环境空气附加氡浓度的有效控制,回风井氡排放量的监测成为了各铀矿开采期间必要的工作。目前对于铀矿回风井氡排放量的测量,都是采用常规测量方法,常规测量方法是采用抓取样方法测量回风井出口的氡浓度,再将该氡浓度乘以回风井出口截面积、风速及排放时间计算出氡的排放量。由于铀矿回风井氡的排放速率随四季而变化,且还与铀矿井下工况相关,为及时掌控铀矿回风井氡的排放量,就需连续或频繁地对回风井出口进行氡浓度监测。然而目前普遍应用的常规测量方法未能同时解决可靠性和便捷性等方面的问题,还不能满足铀矿辐射防护与评价的要求。因此开发出一种新的铀矿回风井氡排放量的测量方法是十分必要的,对于准确评估回风井流出氡的环境影响、保障铀矿矿工和公众健康、保护生态环境具有重大意义。
发明内容
本发明的目的是克服现有技术的上述不足而提供一种测量铀矿回风井氡排放量的方法。
本发明的技术方案是:一种测量铀矿回风井氡排放量的方法,它是采用多点插值积分转换方法进行测量,首先采用多点法对铀矿回风井井口特定样点进行氡浓度测量,然后利用插值积分法计算出回风井井口径向特定区间的平均出口风速和平均氡浓度,再将各区间平均风速、氡浓度平均值及特定区间面积三者的乘积进行加权计算出回风井的氡释放率;
由于回风井口内流场和氡浓度分布极不均匀,然而对于某特定回风井,其排放率与某截面井中心风速、井中心氡浓度及该截面面积之间存在恒定的数值关系,利用该规律和上述测量数据能够计算得到该回风井的转换系数,因此在铀矿回风井氡释放率的后续测量中,只需测得回风井某截面井中心处氡浓度及风速并乘以转换系数和回风井出口截面积即能够计算出回风井的氡释放率,最后再乘以回风井氡的排放时间即能够计算出铀矿回风井氡的排放量。
其具体测量方法如下:
A、转换系数K的测量:首先采用多点法测得回风井井口内某截面i个取样点处的风速Vi及氡浓度Ci,其中i个取样点与井口截面中心点对称分布;
根据测得的i个取样点处风速值,利用插值法推导出风速V与回风井径向距离x的表达公式:
(1)
式中V是回风井径向距离为x处的风速,单位为m/s;x是取样点i到i+1之间任意一点到回风井前端井壁之间的径向距离,单位为m;xi是第i个取样点在回风井径向方向上与前端井壁的距离,单位为m;Vi是第i个取样点对应的风速值,单位为m/s;
对各点风速在回风井径向上不同区间进行积分求平均风速值Vi,积:
(2)
式中Vi,积是第i个取样点积分区间内的平均风速值,单位m/s:
同理,根据测得的i个样点处的氡浓度值Ci,利用插值法推导出氡浓度C与回风井径向距离x的表达公式:
(3)
式中C是回风井径向距离为x处的氡浓度,单位为Bq/m3;Ci是第i个样点所对应的氡浓度值,单位为Bq/m3;
对各点氡浓度在回风井径向上不同区间进行积分求平均氡浓度值Ci,积:
(4)
式中Ci,积是第i个取样点积分区间内的平均氡浓度,单位为Bq/m3;
根据风速与氡浓度积分区域面积公式计算出不同积分区域的面积Si:
(5)
式中Si是第i个取样点处积分区域的面积,单位为m2;ri、ri+1是同一截面上不同区间边界到回风井中心的距离,单位为m;
利用(2)、(4)式计算出的平均风速及平均氡浓度,根据公式:
(6)
计算出铀矿回风井氡释放率,式中P是铀矿回风井氡的释放率,单位为Bq/s;
研究表明,铀矿回风井圆心处氡浓度、风速及井口截面积三者的乘积与回风井氡释放率之间存在着恒定的数值关系:
(7)
式中K为转换系数;C0是铀矿回风井圆心处氡浓度,单位为Bq/m3;V0是铀矿回风井圆心处风速,单位为m/s;S为回风井出口截面积,单位为m2;
B、铀矿回风井氡的排放量的测量:对于铀矿回风井同一截面,其转换系数K为固定常数,因此在该铀矿后续测量过程中,只需测定回风井圆心一点处的风速、氡浓度及已知的转换系数K值即能够利用公式(7)计算出铀矿回风井在该时间段内的氡释放率,最后根据公式:
(8)
即能够计算出铀矿回风井氡的排放量,从而达到快速、方便计算出铀矿回风井氡排放量的目的,式中Q是铀矿回风井氡的排放量,单位为Bq;t是排放时间,单位为s。
本发明与现有技术相比具有如下特点:
(1)本发明通过首先测量计算转换系数K,再利用转换系数K实现对铀矿回风井氡排放量的可靠便捷测量,很好的实现了测量准确性、可靠性与便捷性的有机统一,为铀矿回风井流出物氡的长期可靠监测提供了技术手段。
(2)本发明提供的方法在对铀矿回风井流出物氡的长期监测过程中,只需测量回风井井中心的风速值与氡浓度值,与目前单点或多点测量计算方法相比,具有更小取样不确定性和统计涨落,因此该方法具有更高的准确度
以下结合附图和具体实施方式对本发明的详细结构作进一步描述。
附图说明
附图1为九点法取样点分布示意图,图中i1~i9为取样点,x1~x9为分别为各取样点到回风井井壁之间的径向距离。
具体实施方式
一种测量铀矿回风井氡排放量的方法,其具体测量方法如下:
本实施例采用九点法测量,首先采用九点法测得回风井1井口内某截面各点处的风速Vi及氡浓度Ci,其中九个取样点离回风井1井壁距离分别为:x1为3%d;x2为10%d;x3为20%d;x4为30%d;x5为50%d;x6为70%d;x7为80%d;x8为90%d;x9为97%d;其中d为回风井1出口直径,单位为m。
根据测得的九个取样点处风速值,利用插值法推导出的风速V与回风井1径向距离x的表达公式:
(1)
得到:
(1-1) (1-2)
(1-3)
(1-4)
(1-5)
(1-6)
(1-7)
(1-8)
式中V是回风井1径向距离为x处的风速,单位为m/s;x是待测点在回风井1径向方向上与前端井壁的距离,单位为m; V1、V2、V3、V4、V5、V6、V7、V8、V9分别是利用九点法测得九个样点对应的风速值,单位为m/s。
对各点风速在回风井1径向上不同区间进行积分求平均风速值Vi,积:
(2)
同理,根据测得的九个样点处的氡浓度值Ci,利用插值法推导出的氡浓度C与回风井1径向距离x的表达公式:
(3)
得到:
(3-1)
(3-2)
(3-3)
(3-4)
(3-5)
(3-6)
(3-7)
(3-8)
式中C是回风井1径向距离为x处的氡浓度,单位为Bq/m3;C1、C2、C3、C4、C5、C6、C7、C8、C9分别为九点法测得九个样点所对应的氡浓度值,单位为Bq/m3。
对各点氡浓度在回风井1径向上不同区间进行积分求平均氡浓度值Ci,积。
(4)
根据风速与氡浓度积分区域面积公式计算出不同积分区域的面积Si:
(5)
式中Si是i点处积分区域的面积,单位为m2;ri、ri+1是同一截面上不同区间边界到回风井1中心的距离,单位为m。
ri分别为:r1为0.5d;r2为0.4d;r3为0.3d;r4为0.2d;r5为0;r6为0.2d;r7为0.3d;r8为0.4d;r9为0.5d。
利用(2)、(4)式计算出的平均风速及平均氡浓度,根据公式:
(6)
计算出铀矿回风井1氡释放率,式中P是铀矿回风井1氡的释放率,单位为Bq/s。
研究表明,铀矿回风井1圆心处氡浓度、风速及井口截面积三者的乘积与回风井1氡释放率之间存在着恒定的数值关系:
(7)
式中K为稳定转换系数;C0是铀矿回风井1圆心处氡浓度,单位为Bq/m3;V0是铀矿回风井1圆心处风速,单位为m/s;S为回风井1出口截面积,单位为m2,对于铀矿回风井1同一截面,其转换系数K为固定常数,因此在该铀矿后续测量过程中,只需测定回风井1圆心一点处的风速、氡浓度及已知的转换系数K值即能够利用公式(7)计算出铀矿回风井1在该时间段内的氡释放率,最后根据公式:
(8)
即能够计算出铀矿回风井1氡的排放量,从而达到快速、方便计算出铀矿回风井氡排放量的目的。
常规方法与多点插值积分转换方法对比试验:
在某铀矿A回风井进行了多点插值积分转换方法现场应用,该回风井井径为3.5米,回风井的入口风速为8m/s,试验分别采用常规方法与多点插值积分转换方法对铀矿A回风井氡排放率进行多批次测量,并将两种方法的测量结果进行对比分析,验证多点插值积分转换方法的可靠性。
第一步骤:
利用九点监测法对某铀矿A回风井某一截面上各取样点进行风速及氡浓度的测量,风速测量结果如表一所示:
表一 铀矿A回风井某截面上各取样点风速测量结果
位置 | i1 | i2 | i3 | i4 | i5 | i6 | i7 | i8 | i9 |
风速(m/s) | 6.88 | 7.22 | 7.70 | 8.83 | 9.66 | 8.59 | 7.23 | 6.65 | 6.39 |
回风井各取样点氡浓度测量结果如表二所示:
表二 铀矿A回风井某截面各取样点氡浓度测量结果
位置 | i1 | i2 | i3 | i4 | i5 | i6 | i7 | i8 | i9 |
氡浓度(Bq/m3) | 50450 | 23500 | 23800 | 42100 | 52015 | 46100 | 25200 | 26500 | 46300 |
首先将表一中各取样点风速值代入公式:
(1)
计算出V的表达式,并将该表达式按照不同积分区间分别代入积分公式:
(2)
即可计算出不同积分区间的平均风速Vi,积;
同理将表二中各取样点的氡浓度代入公式:
(3)
计算出C的表达式,并将该表达式按照不同区间分别代入积分公式:
(4)
即可计算出不同积分区间的平均氡浓度Ci,积;
计算出不同积分区间边界到同一截面中心点的距离,ri取值分别为:r1为1.75;r2为1.4;r3为1.05;r4为0.7;r5为0;r6为0.7;r7为1.05;r8为1.4;r9为1.75,将不同的ri分别代入公式:
(5)
计算出该截面上不同积分区间的面积Si;
将计算出的Vi,积和Ci,积及Si代入公式:
(6)
计算出氡的释放率为P=3.07×106Bq/s,在九个取样点中i5为截面中心点,分别由表一、表二可知该截面中心点处的风速值和氡浓度值,并且由回风井直径可计算出该截面的面积;
将中心点风速值V0=9.66m/s、中心点氡浓度值C0=52015Bq/m3、截面面积S=9.6m2及氡释放率P=3.07×106Bq/s代入公式:
(7)
即计算出转换系数K=0.637 ;
第二步骤:
对该铀矿A回风井某截面风速及氡浓度进行第二批次、第三批次及第四批次的测量,获得中心点及其他八个取样点处的风速和氡浓度,其中该截面各取样点风速的测量结果如表三所示:
表三 同一截面不同批次各取样点的风速值
该截面不同批次各取样点氡浓度测量结果如下表四所示:
表四 同一截面不同批次各取样点氡浓度值
采用本发明提供的方法和A回风井排放量的转换系数K,将表三、表四中回风井中心点风速、氡浓度代入公式(7),得到:
P=K×C0×V0×S
分别计算出不同批次铀矿A回风井氡的释放率,同时采用常规测量方法,根据表三、表四中各取样点的风速及氡浓度值进行铀矿回风井氡释放率的计算,两种方法的测量结果和相对标准偏差如表五所示:
表五 不同方法下氡释放率的测量结果和相对标准偏差
由表五可知采用多点插值积分转换法所测量的结果与常规方法测量的结果之间具有较小的相对标准偏差,说明多点插值积分转换方法能实现对铀矿回风井氡排放量的快捷、准确估算。
Claims (2)
1.一种测量铀矿回风井氡排放量的方法,其特征是:它是采用多点插值积分转换方法进行测量,首先采用多点法对铀矿回风井井口特定样点进行氡浓度测量,然后利用插值积分法计算出回风井井口径向特定区间的平均出口风速和平均氡浓度,再将各区间平均风速、氡浓度平均值及特定区间面积三者的乘积进行加权计算出回风井的氡释放率;
研究表明:回风井口内流场和氡浓度分布极不均匀,然而对于某特定回风井,其释放率与某截面井中心风速、井中心氡浓度及该截面面积之间存在恒定的数值关系,利用该规律和上述测量数据能够计算得到该回风井的转换系数,因此在铀矿回风井氡释放率的后续测量中,只需测得回风井某截面井中心处氡浓度及风速并乘以转换系数和回风井出口截面积即能够计算出回风井的氡释放率,最后再乘以回风井氡的排放时间即能够计算出铀矿回风井氡的排放量。
2.如权利要求1所述的一种测量铀矿回风井氡排放量的方法,其特征是:其具体测量方法如下:
A、转换系数K的测量:首先采用多点法测得回风井井口内某截面i个取样点处的风速Vi及氡浓度Ci,其中i个取样点与井口截面中心点对称分布;
根据测得的i个取样点处风速值,利用插值法推导出风速V与回风井径向距离x的表达公式:
(1)
式中V是回风井径向距离为x处的风速,单位为m/s;x是取样点i到i+1之间任意一点到回风井前端井壁之间的径向距离,单位为m;xi是第i个取样点在回风井径向方向上与前端井壁的距离,单位为m;Vi是第i个取样点对应的风速值,单位为m/s;
对各点风速在回风井径向上不同区间进行积分求平均风速值Vi,积:
(2)
式中Vi,积是第i个取样点积分区间内的平均风速值,单位m/s:
同理,根据测得的i个样点处的氡浓度值Ci,利用插值法推导出氡浓度C与回风井径向距离x的表达公式:
(3)
式中C是回风井径向距离为x处的氡浓度,单位为Bq/m3;Ci是第i个样点所对应的氡浓度值,单位为Bq/m3;
对各点氡浓度在回风井径向上不同区间进行积分求平均氡浓度值Ci,积:
(4)
式中Ci,积是第i个取样点积分区间内的平均氡浓度,单位为Bq/m3;
根据风速与氡浓度积分区域面积公式计算出不同积分区域的面积Si:
(5)
式中Si是第i个取样点处积分区域的面积,单位为m2;ri、ri+1是同一截面上不同区间边界到回风井中心的距离,单位为m;
利用(2)、(4)式计算出的平均风速及平均氡浓度,根据公式:
(6)
计算出铀矿回风井氡释放率,式中P是铀矿回风井氡的释放率,单位为Bq/s;
研究表明,铀矿回风井圆心处氡浓度、风速及井口截面积三者的乘积与回风井氡释放率之间存在着恒定的数值关系:
(7)
式中K为转换系数;C0是铀矿回风井圆心处氡浓度,单位为Bq/m3;V0是铀矿回风井圆心处风速,单位为m/s;S为回风井出口截面积,单位为m2;
B、铀矿回风井氡的排放量的测量:对于铀矿回风井同一截面,其转换系数K为固定常数,因此在该铀矿后续测量过程中,只需测定回风井圆心一点处的风速、氡浓度及已知的转换系数K值即能够利用公式(7)计算出铀矿回风井在该时间段内的氡释放率,最后根据公式:
(8)
即能够计算出铀矿回风井氡的排放量,从而达到快速、方便计算出铀矿回风井氡排放量的目的,式中Q是铀矿回风井氡的排放量,单位为Bq;t是排放时间,单位为s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510320774.7A CN104880726B (zh) | 2015-06-12 | 2015-06-12 | 测量铀矿回风井氡排放量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510320774.7A CN104880726B (zh) | 2015-06-12 | 2015-06-12 | 测量铀矿回风井氡排放量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104880726A true CN104880726A (zh) | 2015-09-02 |
CN104880726B CN104880726B (zh) | 2017-07-11 |
Family
ID=53948292
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510320774.7A Active CN104880726B (zh) | 2015-06-12 | 2015-06-12 | 测量铀矿回风井氡排放量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104880726B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114526118A (zh) * | 2022-02-18 | 2022-05-24 | 中国矿业大学(北京) | 巷道断面瓦斯流量测定方法和装置、存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201876550U (zh) * | 2010-12-03 | 2011-06-22 | 深圳市清华环科检测技术有限公司 | 采氡仪 |
JP2011209080A (ja) * | 2010-03-29 | 2011-10-20 | Central Res Inst Of Electric Power Ind | ラドンフラックスの測定方法 |
JP2011209079A (ja) * | 2010-03-29 | 2011-10-20 | Central Res Inst Of Electric Power Ind | ラドンによる岩石内部の間隙表面積測定方法及びラドンによる岩石内部の間隙表面積測定装置 |
JP2013047683A (ja) * | 2012-10-01 | 2013-03-07 | Mitsubishi Electric Corp | 放射性ヨウ素モニタ装置 |
RU2519079C1 (ru) * | 2013-01-09 | 2014-06-10 | Айрат Муратович Сафаров | Способ экспресс-определения загрязнения участков почв и подземных вод нефтью и нефтепродуктами |
-
2015
- 2015-06-12 CN CN201510320774.7A patent/CN104880726B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011209080A (ja) * | 2010-03-29 | 2011-10-20 | Central Res Inst Of Electric Power Ind | ラドンフラックスの測定方法 |
JP2011209079A (ja) * | 2010-03-29 | 2011-10-20 | Central Res Inst Of Electric Power Ind | ラドンによる岩石内部の間隙表面積測定方法及びラドンによる岩石内部の間隙表面積測定装置 |
CN201876550U (zh) * | 2010-12-03 | 2011-06-22 | 深圳市清华环科检测技术有限公司 | 采氡仪 |
JP2013047683A (ja) * | 2012-10-01 | 2013-03-07 | Mitsubishi Electric Corp | 放射性ヨウ素モニタ装置 |
RU2519079C1 (ru) * | 2013-01-09 | 2014-06-10 | Айрат Муратович Сафаров | Способ экспресс-определения загрязнения участков почв и подземных вод нефтью и нефтепродуктами |
Non-Patent Citations (3)
Title |
---|
DONG XIE ET AL: "Radon dispersion modeling and dose assessment for uranium mine", 《JOURNAL OF ENVIRONMENTAL RADIOACTIVITY》 * |
冀东 等: "铀矿采场氡浓度模拟初步探讨", 《环境科学与技术》 * |
王志勇 等: "铀矿区放射性核素氡在大气中的迁移研究", 《安全与环境学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114526118A (zh) * | 2022-02-18 | 2022-05-24 | 中国矿业大学(北京) | 巷道断面瓦斯流量测定方法和装置、存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN104880726B (zh) | 2017-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103643997B (zh) | 一种小风速掘进工作面的瓦斯涌出量的计算方法 | |
CN105545284B (zh) | 一种随钻伽马成像数据处理方法 | |
CN109680669B (zh) | 一种基于地下水动力加载的基坑稳定性监测预警方法 | |
CN102305612A (zh) | 一种位移/挠度测量系统与方法 | |
GB2575940A (en) | Apply multi-physics principle for well integrity evaluation in a multi-string configuration | |
CN102645385A (zh) | 一种不同参数对材料疲劳裂纹扩展性能影响的检测方法 | |
CN109507743B (zh) | 一种高精度扫描探雷装置及扫描检测方法 | |
WO2013155075A1 (en) | Heap monitoring | |
CN103778480A (zh) | 一种基于敏感度分析的裂隙带高度预测方法 | |
CN106446445A (zh) | 一种巷道平均风速的单点统计测量方法 | |
CN108151766A (zh) | 磁钉的定位方法、磁钉定位导航误差修正方法及定位装置 | |
CN106991509A (zh) | 基于径向基函数神经网络模型的测井曲线预测方法 | |
CN108470095A (zh) | 基于数据驱动径向基函数模型的tbm推进力预测方法 | |
CN107016620A (zh) | 一种基于层次分析的断层含水性评价方法 | |
CN110145301A (zh) | 一种随钻密度成像测井仪环境校正方法 | |
CN105651311A (zh) | 农机作业卫星导航自动驾驶精度的测试方法 | |
CN103984017B (zh) | 一种微震震源自动定位方法 | |
CN103954987A (zh) | 一种利用γ相机与γ谱仪联用测量放射性物质质量的方法 | |
CN104880726A (zh) | 测量铀矿回风井氡排放量的方法 | |
CN106777556B (zh) | 一种评估边坡开挖期稳定状态的空间分析方法 | |
CN206311626U (zh) | 一种地下水流速和流向探测装置 | |
CN103323571B (zh) | 粮仓粮食储藏的粮堆中虫霉活动定位的方法 | |
CN103033274B (zh) | 一种日照温度概率密度的测定方法 | |
CN103954197B (zh) | 一种巷道表面位移、深部位移同位测量方法 | |
CN103325065A (zh) | 农产品质量安全检测抽样决策方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |